Abstract
Understanding the features of the paleoclimate and paleoenvrionment from key areas is crucial for predicting the climate variation of the future. In this study, fossil charcoal coupled with high-accuracy accelerator mass spectrometry (AMS) 14C dating from the Guanzhong Basin is analyzed to reconstruct the paleoclimate. Here, the coexistence approach (CA) was applied, and the result showed that the mean annual temperature (MAT) was about 14.8°C, and the annual precipitation (AP) was about 831.1 mm in the Guanzhong Basin during 6200–5600 cal. a BP. Comparing the climate between the mid-Holocene and present in the Xi’an area, the MAT was about 1.1°C higher than today and the AP was about 278 mm higher than today, similar to the modern climate of the Hanzhong area in the southern Qinling Mountains.
Introduction
‘Global Warming’ has a profound impact on governments and people worldwide; however, the warming extent and its climatic–ecological effects are still uncertain. In order to understand and predict climate changes in the future, one approach is to study the features of the paleoclimate and paleoenvironment variation from key areas, identify historical analogues of predicted warming, and reconstruct the paleotemperature and paleoprecipitation from robust climatic proxies.
The early-Holocene experienced increased temperatures, and then, a warm and humid mid-Holocene was followed by a cooling trend of the late-Holocene (Marcott et al., 2013; Wang et al., 2005). The Holocene Megathermal occurred in the mid-Holocene (An et al., 2000; Renssen et al., 2009; Salonen et al., 2011; Wanner et al., 2008), when the mean global temperature was about 1–2°C higher than today and is similar to scenarios of future global warming (IPCC, 2013). The Megathermal Maximum in China occurred during 8–6.4 ka BP when the temperature was about 1.5°C higher (Fang and Hou, 2011) and the precipitation was about 50–300 mm higher than today (Fang et al., 2011), but these values are very generalized. The warming extent varies from eastern to western China (Ge et al., 2006), indicating regional differences (Steig, 1999). Therefore, it is important to obtain high-resolution climatic records from different areas, which can help us better understanding the processes and patterns of climate change (Broecker, 1997; Steig, 1999).
The Chinese Loess Plateau (CLP) is in a semi-arid and sub-humid area of middle and western China and the climate is strongly controlled by the East Asian monsoon. The ecosystems are fragile and sensitive to climate change (Gong and Wang, 2005), and thus, this area is significant for carrying out research on response to climate changes. To date, some achievements have been made. A pollen and lake record from the western CLP indicates that the wettest period in the Holocene appears during 7.5–5.9 cal. ka BP (An et al., 2003). Fossil charcoal from the Tianshui Basin shows that the mean annual temperature (MAT) was approximately 13.2°C, and the annual precipitation (AP) was approximately 778 mm between 5.2 and 4.3 cal. ka BP (Sun and Li, 2012). Magnetic susceptibility analysis from Luochuan, which is located in the middle of the CLP, suggests that the MAT reached a maximum of 12.2 ± 1.4°C and the AP peaked at 701 ± 74 mm during the Holocene (Lu et al., 1994). Organic stable carbon isotopic composition data from a loess-sand sequence which is located to the north of the CLP indicate that the AP was 435 mm after 8 ka (Chen et al., 2015). Nevertheless, the Loess Plateau covers a large area, the climate character varies from west to east and from north to south, and paleoclimate data from each part of the CLP will make a contribution toward a more thorough understanding of the climate change mechanisms in this region. The Guanzhong Basin is located on the southern CLP, and it covers about 300 km from west to east. In the west Guanzhong Basin, plant phytolith assemblies from Baoji indicate that the MAT was 14–16°C and the AP was 700–800 mm during the Holocene (Lu et al., 1996). In the east Guanzhong Basin, the geochemical transfer functions from Weinan show that the MAT was about 12–14°C and the AP was above 700 mm during the Holocene optimum (Sun et al., 1999a). However, data on paleoclimates from the hinterland were insufficient; this study will broaden the paleoclimate database of the CLP.
Botanical records are important proxies for paleoclimate studies; for example, pollen (Huntley, 1993; Seppä and Birks, 2001; Seppä et al., 2004), plant macrofossils (Yao et al., 2011), fruits (Grein et al., 2011), and fossil charcoal (Sun and Li, 2012) have been used. Pollen sites for analysis are rare in the CLP, and it has proven difficult to quantitatively reconstruct the paleoclimate of the region based on pollen. Fortunately, Neolithic culture and early agriculture flourished in the CLP, and many relics such as fossil charcoal are well preserved in Neolithic sites. Fossil charcoal comes from the incomplete burning of wood, and the anatomical characters of the original wood are retained (McGinnes et al., 1974); this raises the possibility of much greater precision in the level of taxonomic identification, overcoming the limitation of some plant microfossils, thus improves the reliability of paleoclimate reconstruction. Therefore, fossil charcoals are ideal material for this study.
The Holocene Megathermal refers to the warmest interglacial period after the last glacial, and in this period, agriculture-based Neolithic culture developed rapidly and the Yangshao culture reached a peak in Neolithic prosperity. In order to find out the possible connections between environment and human cultures, two Yangshao culture (mid-Holocene) sites from the Guanzhong Basin around Xi’an area were selected; the fossil charcoal, with high-precision accelerator mass spectrometry (AMS) 14C dating, was analyzed to quantify the paleoclimate for the period.
Study area
The Guanzhong Basin is located in the middle region of the Yellow River valley; it is bounded on the west by Baoji, the east by Tongguan, and the south by the Qinling Mountains. Because of its special geographic setting, a warm–humid climate developed in this area made it an ideal place for settlements. The Neolithic cultures include the Laoguantai, Yangshao, and the Longshan cultures, and they were well developed in the Guanzhong Basin, making the Guanzhong Basin recognition as the cradle of Chinese civilization. Xi’an is located in the hinterland of the Guanzhong Basin, with annual temperature of 13.7°C and AP of 553.3 mm, and two Yangshao culture sites including Huiduipo and Beiniu sites near Xi’an were selected for this study.
The Huiduipo site (34°34′4.1″N, 109°01′41.8″E) located in the Huiduipo village of Gaoling county, in Xi’an city (Figure 1), is a site that covers an area of 60,000 m2, and the culture layer is about 5 m thick. Relics such as red and black clay pottery are common in the site. A 300-cm-deep sediment profile with abundant fossil charcoals from the southern site was selected for this study, and three layers were divided according to the color and lithologic features of the sediment (Figure 2): above 20 cm is the modern cultivation layer, containing grass roots; between 20 and 240 cm, there is a brownish culture layer, containing abundant charcoal fragments and pottery pieces; below 240 cm is yellowish unconsolidated loess.

Study area.

Stratigraphic sections from the Huiduipo and the Beiniu sites.
The Beiniu site (34°28′14.5″N, 109°19′2.6″E) is located in Beiniu village of Lintong county, in Xi’an city, and is close to the southern bank of the Weihe River (Figure 1). Extensive excavation projects have been conducted at the Beiniu site by the Shaanxi Provincial Institute of Archaeology, and its area of 200,000 m2 has relics such as gray and red pottery (Sun et al., 2006). A 300-cm-deep sediment profile with abundant fossil charcoals from the eastern site was selected for this study, and eight layers were divided according to the color and lithologic features of the sediment (Figure 2): (1) 0–10 cm: modern cultivation layer, containing grass roots; (2) 10–50 cm: yellowish unconsolidated loess; (3) 50–70 cm: blackish cultural layer, containing many charcoal fragments and pottery pieces; (4) 70–100 cm: brownish cultural layer, containing a few charcoal fragments and pottery pieces; (5) 100–200 cm: blackish cultural layer; (6) 200–210 cm: brownish cultural layer; (7) 210–280 cm: blackish cultural layer; (8) 280–300 cm: yellowish unconsolidated loess.
Methods
Two sections with continuous and undisturbed cultural sediment in the Huiduipo and Beiniu sites were selected; the fossil charcoal and seeds within the sediment were recovered using the floatation method (Tsuyuzaki, 1994). In this study, nine AMS radiocarbon dates were measured at the Australian Nuclear Science and Technology Organization (ANSTO) based on recovered seeds, and the conventional 14C ages were then converted to calibrated ages using the IntCal09 dataset (Reimer et al., 2009) (Table 1). Other ages were inferred from interpolation between dated samples.
Accelerator mass spectrometry (AMS) radiocarbon dates from the Huiduipo (HDP) and the Beiniu (BN) sites.
Based on a statistical analysis of fossil charcoal from temperate regions (Keepax, 1988; Smart and Hoffman, 1988; Sun and Li, 2012), a minimum of about 100 pieces of fossil charcoal from each sample is considered suitable for reliable analysis. Fossil charcoal samples were selected at random, following standard procedures: first, pressure fractured charcoal fragments were prepared with a razor blade, in order to produce fresh clean surfaces to show transverse, radial, and tangential sections (Leney and Casteel, 1975). These were examined under a stereomicroscope, categorized, and one or two samples from each type were photographed under a scanning electron microscope (SEM). Identification of the taxa was carried out using wood anatomy atlases (Cheng et al., 1985; Schweingruber, 1990; Yao, 1988; Yao et al., 2002).
The coexistence approach (CA) is a method for quantitative reconstruction of paleoclimate (Mosbrugger and Utescher, 1997), and it performs well for lowland and mid-altitude areas (Grimm and Denk, 2012), the reliability of which has been validated by previous studies (Sun and Li, 2012; Yang et al., 2007). It is based on the assumption that fossil plants have similar climatic requirements as their nearest living relative species (NLRS), and while the climatic requirements for each NLRS may be different, the overlap interval of the climatic data from different fossil plants is supposed to be the paleoclimate when the fossil plants co-existed (Mosbrugger and Utescher, 1997). Since plants in the Holocene are the result of long-term natural selection, their ecological amplitude and climatic tolerance are assumed to be the same as their living plant counterparts; thus, the CA is well suited for obtaining quantitative information about past climate.
Most published CA analyses rely on climate data for NLRS stored in the paleoflora database (PFDB). However, this database is mainly from North American and European species and is incomplete for China which is a region of megadiversity. The distribution characteristic of modern plants in China is well known, and distribution maps for more than 11,000 species at the level of counties are documented in the ‘Atlas of woody plants in China’ (Fang et al., 2009). This provides information on both species and environmental variables, including MAT (°C), annual biotemperature (ABT, °C), warmth index (WI, °C), coldness index (CI, °C), AP (mm), and seven other climatic variables (Fang et al., 2009; Kou, 2005). Here, two main climatic factors are calculated: MAT and AP based on the database from the atlas.
Results
Huiduipo site
The AMS radiocarbon dates from the Huiduipo section show that the sediment was deposited between 6500 and 5600 cal. a BP (Table 1). In order to obtain sufficient charcoal fragments for quantitative analysis, six samples were selected from the sediment between 20 and 200 cm depth that was deposited between 6300 and 5600 cal. a BP.
A total of 516 charcoal fragments from six samples were examined, and 17 different taxa were identified with one unknown species (Table 2). The abundance of Quercus sp. (14.92%), Fargesia sp. (13.37%), Ulmus sp. (12.02%), Sassafras sp. (9.88%), and Salix sp. (9.11%) make up about 60% of the fossil charcoal assemblages, while Abies sp. (7.56%), Vitis sp. (6.59%), and Toxicodendron sp. (4.84%) have lower abundance. It is worth noting that five subtropical taxa Castanopsis sp., Fargesia sp., Sassafras sp., Stewartia sp., and Toxicodendron sp. appear in the charcoal assemblages.
The fossil charcoal assemblages from the Huiduipo site during 6300–5600 cal. a BP.
Beiniu site
The AMS radiocarbon dates and the interpolation dates from the Beiniu section show that the sediment was deposited between 6200 and 4700 cal. a BP (Table 1). A total of 704 charcoal fragments from five samples were examined, and 15 different taxa were identified (Table 3). Almost all the taxa appeared in the earlier three samples (493 charcoal fragments) dated between 6200 and 5600 cal. a BP, while many taxa disappeared in the later two samples (211 charcoal fragments) dated between 5600 and 4700 cal. a BP; here, our study will focus on the samples between 6200 and 5600 cal. a BP.
The fossil charcoal assemblages from the Beiniu site during 6200–5600 cal. a BP.
Table 3 shows that the abundance of Paliurus sp. (28.2%), Quercus sp. (27.18%), and Pinus sp. (11.56%) make up about 67% of the fossil charcoal assemblages, and Cyclobalanopsis sp. (5.68%), Platycarya sp. (5.68%), and Zelkova sp. (5.07%) are lower in abundance. Also, five subtropical taxa including Castanopsis sp., Cyclobalanopsis sp., Fargesia sp., Melia sp., and Platycarya sp. appear in the charcoal assemblages.
The Huiduipo and Beiniu sites are located in Xi’an – a hinterland region of the Guanzhong Basin, and they are in close proximity (Figure 1); as mentioned above, the AMS 14C dating of the samples from the two sites used here is 6300–5600 and 6200–5600 cal. a BP, respectively. For this reason, the climate and environment of the two sites are assumed to be similar during 6200–5600 cal. a BP, and the plant taxa found in the two sites is due to environmental preferences related to climate; therefore, the plant taxa from both the Huiduipo and Beiniu site will reflect the climate in the Guanzhong Basin between 6200 and 5600 cal. a BP.
After combining the fossil charcoal taxa from the two sites, we observe a total of 27 taxa (duplicate records were deleted) appeared in the study area during the mid-Holocene, and the climatic factors were obtained from the CA analysis based on these fossil charcoal taxa. The CA results show that the MAT was 9–15.8°C, and the AP was 760–1478 mm in the Guanzhong Basin during 6200–5600 cal. a BP (Figure 3).

Coexistence intervals for Xi’an during 6200–5600 cal. a BP.
Discussion
Evaluation and validation on the quantitative reconstructions
The Guanzhong Basin is located in the transition area between the Qinling Mountains and the CLP (Figure 1). The Qinling Mountains lie across central China and is an important dividing line of vegetation and climate, with subtropical climate in the south and warm temperate climate in the north; accordingly, subtropical evergreen broadleaf forest is developed on the south and temperate deciduous forest on the north slopes. The modern climate of the CLP is governed by the East Asian monsoon system, and it shows a regular pattern from north to south; the climate changes gradually from cold and dry to warm and wet, and this brings a corresponding change in the vegetation type from north to south. In order to find the relationship between vegetation and climate in the study area, and verify the validity of the CA in quantitative reconstruction of the paleoclimate, CA is applied to the modern vegetation types of different climatic zones to obtain the estimated value of climatic factors, which will be effectively compared with the actual climatic factors.
According to the climate character of the CLP, a belt around the Guanzhong Basin from the desert/loess transition zone to the south slopes of the Qinling Mountains was selected; five places including Shenmu (which is located in Yulin), Yan’an, Xi’an, Hanzhong, and Ankang (Figure 1) were chosen to carry out a vegetation–climate relationship study. In this paper, most of the charcoal fragments were identified to genus and a few to species level; hence, different genera of the local modern woody plants from the five places were selected according to local chronicles and the Flora of China (Ankang Local Chronicles Compilation Committee, 1989; Du, 2009; Flora of China Editorial Committee of Chinese Academy of Sciences, 2004; Hanzhong Local Chronicles Compilation Committee, 1994; Shenmu County Chronicles Compilation Committee, 1990; Xi’an Local Chronicles Compilation Committee, 1996; Yan’an Local Chronicles Compilation Committee, 1994), related modern climatic factors come from the Surface Meteorological Data of China (1971–2000, Table 4).
The modern climate and estimated climate from CA.
MAT: mean annual temperature; CA: coexistence approach; AP: annual precipitation.
Based on the climate information of each plant from the ‘Atlas of Woody Plants in China’ (Fang et al., 2009), CA analysis was applied to the modern woody plants in the five places to obtain the estimated climatic factors, including the lowest and the highest estimated value of MAT (MATmin, MATmax) and AP (APmin, APmax), and then the correlation between estimates and actual value (MAT, AP) was analyzed (using Statistical Package for the Social Science (SPSS)). Table 5 shows that the Pearson correlation coefficient is 0.986, which indicates that the correlation between the independent variables (MATmax, MATmin) and the dependent variable (MAT) is high; R2 is 0.973, thus about 97% of the variation in the dependent variable can be explained by the variation in the independent variable and thus estimated temperature is close to the actual temperature; the significance level of the F test is 0.027 (Table 6); therefore, the regression equation established here is strong. Based on the regression coefficients (Table 7), the equation showing the relationship is as follows:
Model summary a (for MAT).
MAT: mean annual temperature.
Dependent variable: MAT.
Predictors: constant, MATmax, MATmin.
ANOVA a (for MAT).
ANOVA: analysis of variance; MAT: mean annual temperature.
Dependent variable: MAT.
Predictors: constant, MATmax, and MATmin.
Coefficients a (for MAT).
MAT: mean annual temperature.
Dependent variable: MAT.
Table 8 shows that the Pearson correlation coefficient is 0.992, which indicates that the correlation between the independent variables (APmax, APmin) and the dependent variable (AP) is high; R2 is 0.984, indicating that about 98% of the variation in the dependent variable can be explained by the variation in the independent variable, thus estimated precipitation is close to the actual precipitation; the significance level of the F test is 0.016 (Table 9), hence the regression equation established here is also strong. Based on the regression coefficients (Table 10), the equation is established as follows:
Model summary a (for AP).
AP: annual precipitation.
Dependent variable: AP.
Predictors: constant, APmax, and APmin.
ANOVA a (for AP).
ANOVA: analysis of variance; AP: annual precipitation.
Dependent variable: AP.
Predictors: constant, APmax, and APmin.
Coefficients a (for AP).
AP: annual precipitation.
Dependent variable: AP.
Climate plays an important role in controlling vegetation distribution, and climate change will therefore cause vegetation shifts southwards or northwards (Ni et al., 2010; Prentice et al., 2000; Sun et al., 1999b; Yu et al., 1998, 2000); also, the shift can occur in the belt from the desert/loess transition zone in the north to the south slope of the Qinling Mountains. Thus, natural vegetation types around Xi’an can shift in sympathy with northern (e.g. Yan’an, Shenmu) or southern neighbors (e.g. Hanzhong, Ankang) depending on the climate changes. Since the regression equations applicable to each place in this belt today, we assume that it is applicable for Xi’an during the mid-Holocene.
The estimated climatic factors of the Xi’an area obtained from the CA analysis are as follows: MATmin = 9°C, MATmax = 15.8°C; APmin = 760 mm, APmax = 1478 mm; in order to find the most probable value of the MAT and AP in the study area, the estimated values are substituted into Eqs (1) and (2), respectively, and the result shows that the MAT was about 14.8°C and the AP was about 831.1 mm in Xi’an area between 6200 and 5600 cal. a BP. Thus, in the mid-Holocene, the MAT was about 1.1°C higher than today and the AP was about 278 mm higher than today for Xi’an, these are similar to the modern climate of the Hanzhong area in southern Shaanxi in the southern Qinling Mountains (Figure 1; Table 11).
The comparison of the climatic factors between Xi’an and Hanzhong.
MAT: mean annual temperature; MAP: mean annual precipitation.
Paleoclimate reconstructions from west (Baoji), hinterland (Xi’an), and east (Weinan) of the Guanzhong Basin during the Holocene now exist. After comparing the paleo and modern climate data (Table 12), we see that the temperature was about 1–2°C higher, and precipitation was about 200–300 mm higher than today across the Guanzhong Basin.
The comparison of the paleoclimate during the Holocene and modern climate in the Guanzhong Basin.
MAT: mean annual temperature; MAP: mean annual precipitation.
Implications of paleoclimatic reconstructions
Many studies show that water vapor in the atmosphere closely follows temperature in agreement with the Clausius–Clapeyron relationship. An increase in temperature by 1°C in the lower troposphere implies an increase in the vertical profile of water vapor by 6–7% (Bengtsson, 2010), this will cause an increase in precipitation in the East China Sea and South China Sea (Deng et al., 2014). And it is predicted that the precipitation will also increase where the rise in sea surface temperature exceeds the mean surface warming in the tropics (Huang et al., 2013; Vecchi and Soden, 2007; Xie et al., 2010).
The results of Climate simulations demonstrate that global warming would enhance the thermal depression over the Asian continent in summer, which will further enhance the East Asian Summer Monsoon (EASM). Simultaneously, the increase in water vapor in the tropics area will provide sufficient water vapor resource for the East Asia; thus, future warming will likely cause an increase in precipitation in the EASM region (Sun and Ding, 2010). Precipitation in the CLP is mainly controlled by the water vapor transported by the EASM (An et al., 2000; Chen et al., 1991; Tan et al., 2014); for this reason, increasing temperature may have a positive impact on the precipitation in the CLP.
Human activity and its response to climate change in mid-Holocene
In the long history, food is essential for human survival and development, and the origin of agriculture guarantees that people could have a reliable supply of it. China is the origin and development center of rain-fed millet (Lu et al., 2009; Zhao, 1998) and cultivated rice (Crawford, 2006; Yan, 1997). Rain-fed agriculture was distributed mainly in the Yellow River valley in North China (Barton et al., 2009), and the common millet first appeared around 10 ka BP and the foxtail millet came after that (Lu et al., 2009). The origin and domestication of the millet are considered to be closely related to the typical climate (rain hot during the same period) of the East Asian monsoon area. Because of the sufficient precipitation, the millet cultivation spread rapidly to all over the CLP by 7 ka BP, and subsequent development of millet agriculture resulted in the flourishing of the Yangshao Culture (Huang et al., 2004). The higher temperature and precipitation during the Yangshao period (this paper) created a warm and wet environment which was able to meet the needs of rice cultivation, and this is supposed to be one of the reasons that rice spread into the Guanzhong Basin during this period (Zhang et al., 2010). The emergence of both millet and rice in the Guanzhong Basin indicates that the human farming activity changed from millet cultivation to mixed millet–rice cultivation (dominantly cultivated millets and locally cultivated rice), and this would significantly change the proportion of food production in the Neolithic subsistence economy, thus establish a very strong basis for the formation of the social civilization.
Conclusion
The quantitative reconstruction of the paleoclimate in the hinterland of the Guanzhong Basin suggests that the MAT was 14.8°C, and the AP was about 831.1 mm. Higher temperature corresponds with the higher precipitation, which would favor the growth of subtropical plants taxa and then may cause vegetation shift northward. In this case, a warmer–wetter environment will develop in the south of the CLP if the temperature rises about 1–2°C in the future.
Footnotes
Funding
This work was supported by the Chinese National Natural Science Foundation (grant nos 41202131 and 41372175); Natural Science Basic Research Plan in Shaanxi Province of China (grant no. 2016JQ4026) the Special Fund for Basic Scientific Research of Central Colleges, Chang’an University, China (grant nos 2014G1271059 and 2014G3272013); and Strategic Priority Research Program of the Chinese Academy of Sciences (grant no. XDB03020303).
