Abstract
The reconstruction of past climates and vegetation can provide valuable information for understanding the mechanisms of the variability of the East Asian summer monsoon in eastern China. In this study, organic geochemical evidence from the compositions of sedimentary leaf-wax stable isotopes and n-alkane parameters investigates the changes in vegetation and climate over the last 1200 years in the Xiyaohu peatland, Jiangxi Province, southeast China. Combined with temperature records, three climatic periods are presented: (a) a warm and humid period with an increase in C4 plants from 900 to 1450 AD, which coincides with the Medieval Warm Period (MWP); (b) a cool and dry period with the expansion of C3 plants from 1450 to 1800 AD, coinciding with the Little Ice Age (LIA); and (c) the Present Warm Period (PWP) from 1800 AD until the present, with warm and wet conditions. The sub-stages within the MWP and LIA intervals are also presented. The earlier MWP stage (900–1125 AD) was drier than the latter one (1125–1450 AD), and the earlier LIA stage (1450–1650 AD) was drier than the late LIA (1650–1800 AD). Increased solar irradiance and enhanced El Niño activities are related to the warm and humid climate during the MWP and PWP, whereas reduced solar irradiance and La Niña activities correspond to the cool and dry climate during the LIA. The present results provide insights into paleoclimatic changes in eastern monsoonal China and provide an understanding of centennial-scale climatic fluctuations and their driving factors.
Introduction
The Asian monsoon is a key component of the global climate system. Summer monsoon-derived rainfall is crucial to terrestrial ecosystems and human societies in East Asia (An, 2000). The last 1200 years include the Medieval Warm Period (MWP) and Little Ice Age (LIA), which are multi-centennial-scale global climatic events (Mann et al., 1999). The MWP prevailed in various parts of the world, and historical documents and paleoclimate reconstructions show that the temperature increased during the MWP. Continental and oceanic temperature reconstruction records demonstrate that the LIA is the nearest mountain-glacier expansion episode on a centennial time scale and featured colder conditions than those of medieval times (D’Arrigo et al., 2006; Holzhauser et al., 2005; Kayastha and Harrison, 2008; Lorrey et al., 2014; Matthews and Briffa, 2016). Paleoclimate reconstructions and simulations show that precipitation during the MWP and LIA presents different patterns in China. The wet–dry pattern of the eastern monsoon region of China is unstable and controversial (Chen et al., 2015). The middle Yangtze River in the East Asia monsoon area is typically recognized as one of the origins of Chinese civilization and supports prosperous rice cultivation, which depends primarily on water availability and mild temperatures (Normile, 1997). Climatic variability during the MWP and LIA profoundly affects southeast China (Zhao and Piperno, 2000), therefore, a clear understanding of past monsoonal climate changes and their drivers, particularly during key periods such as the MWP and LIA, is essential to study the response of terrestrial ecosystems to climate change and its influence on human development.
n-alkanes are biomarkers that are composed of a chain of methylene groups and epidermal cells by the elongation of C16 fatty acids. In the absence of functional groups, they are thermally stable and relatively resistant to degradation (Naafs et al., 2019). n-alkane distribution and stable isotopic composition (δD and δ13C) have been widely used to reconstruct various paleoenvironmental changes, such as temperature variability, paleo-vegetation evolution, and hydrological changes (Feakins et al., 2019; Ficken et al., 2000; Wang et al., 2017).
Ombrotrophic peatlands ideal archives for high-resolution paleoclimate studies owing to their high accumulation rates and sensitivity to climate change (Gałka et al., 2017; Schittek et al., 2016). Studies have been conducted in South China using lipid biomarker approaches (Huang et al., 2018; Zheng et al., 2009; Zhou et al., 2005). The hydrogen isotope ratios of n-alkanes from the Poyang Lake sequence were used to reconstruct precipitation changes (Yao et al., 2015). n-alkane distributions and other proxies from a peat deposit from southern China were applied to investigate paleoclimate changes from 18 ka (Zhou et al., 2005). However, high-resolution studies of n-alkane δD and δ13C values in peat from South China are relatively few.
Xiyaohu (XYH) is a relatively elevated peatland in Jiangxi Province in the south of the middle Yangtze River catchment, and multiple proxies from the peatland, including pollen, phytolith, and biomarkers have been used to reconstruct vegetation changes (Cui et al., 2018; Zhang et al., 2019; Pang et al., under review). A previous study conducted a grain-size analysis to study climate change in the past 2000 years at Xiyaohu (Li et al., 2017). In this study, sediment leaf-wax n-alkane distributions, δD, and δ13C, from the XYH peat sequence for the last 1200 years were analyzed. Leaf-wax δD records δD precipitation in the middle Yangtze River, and δ13C typically indicates photosynthetic pathways for C3 and C4 plants (Collister et al., 1994). This study aimed to detect the changes in the botanical sources of the peatland and to reconstruct paleoclimate changes during the past 1200 years. Considering the integration of human activities and climatic changes, the results were compared with previously published proxies and historical documents to analyze different anthropogenic and natural influencing factors.
Study area
The XYH peatland (28.44°N, 115.40°E; Figure 1b) is located on the northwest main peak of Mount Xishan in Jiangxi province, in the middle reaches of the Yangtze River. The bottom of the peatland is flat and covered by a silty clay layer and heavily weathered exposures of granite bedrock (Figure 1c). Rainfall is the primary water source, and hygrophilous vegetation is abundant because of recharge from slope runoff. This provides an ideal situation for the formation of ombrotrophic peatland. The XYH peatland is the largest and thickest at Mount Xishan, with a surface area of 12,000 m2 and a maximum depth of 3.8 m.

(a) Monthly temperature and precipitation from the nearby Nanchang meteorological station. (b) The Location of XYH (red star) and another site discussed in the study (yellow star). The dashed line represents the boundary of the modern summer monsoon, and the arrows represent two storm trajectories (Chen et al., 2010). The pink shadow represents the Jiang-Nan area (Zheng et al., 2006). The base map is provided from the ArcGIS Hub. (c) The General cross-sectional view of the studied peatland in Xishan Mountain (modified from Zhang et al., 2019).
The study area is influenced by the humid subtropical East Asian summer monsoon and experiences summer-dominated precipitation, with maximum precipitation from April to June (Figure 1a). The total annual precipitation is 1400 mm, and the annual mean temperature is 17°C–18°C with a short spring and autumn but a long summer and winter (Bureau, 2015).
The vegetation on Mount Xishan changes with elevation. The foothills contain economic crops, namely two-season croplands, rapeseed, and tea-oil tree forests. The hillsides are dominated by widespread scrubs (i.e. Loropetalum chinense, Caccinium bracteatum, and Rhododendron simsii scrub), grass-forb communities (Arundinella setosa and Miscanthus sinensis), and sporadic forests (Phyllostachys pubescens, Pinus massoniana, Loropetalum chinense, Rhododendron simsii, and Cunninghamia lanceolate). Areas with elevation above 700 m feature the grass group, comprising swamps (Miscanthus sacchariflorus, Phragmites communis, Carex spp., and Juncus effuses), meadows (Carex cinerascens), and scrub (Salix spp. and Lespedeza bicolor). Vegetation in the peatland area includes widespread grasses, sporadic scrubs, meadows, and swamps) (Vegetation Atlas of China, 2001). The dominant vegetation in the valley includes Scirpus cyperius, Salix brachypoda, and Sphagnum assemblages.
Materials and methods
Core sampling and stratigraphy
An XYH sediment core was drilled from the center of the XYH peatland using a manual sampling peat drill with a PVC tube with a length of 300 cm and an inner diameter of 60 mm. The total depth of the sample (300 cm) was divided into three units based on lithology. The top (0–120 cm) of the core consisted gray-black clay with abundant plant residue. From 120 to 215 cm, gray–yellow and dark gray muddy sand is the primary component. The bottom (215–300 cm) changed from taupe to gray-black peat soil (Figure 2).

Lithology of the XYH peatland and age-depth model. Gray dots indicate the 95% probability intervals. Average linear sedimentation rates (sr) are also shown.
Core chronology
All procedures were conducted in the Pollen and Paleoecology Lab, Department of Geography and Ocean Science, Nanjing University. Radiocarbon dating of sporopollen concentrates, plant residues, and charcoal samples (total of 12) was conducted by acceleration mass spectrometry (AMS) at the Beta Analytic Laboratory (Florida USA; Table 1). The Bayesian software Bacon “Winbacon2.2” (Blaauw and Christeny, 2011) and R statistical software were used to develop an age model using 12 AMS radiocarbon dates (Blaauw and Christeny, 2011). CALIB Rev. 5.0.2 was used to calibrate the radiocarbon age to calendar years. Generally, the section from 86 to 120 cm displayed notably lower sedimentation rates than the one from 153 to 300 cm. The top section from 0 to 86 cm has a medium sedimentation rate.
AMS 14C dating and correct results of XYH core.
Lipid extraction
The core was selected at 5- and 4-cm intervals through the top and bottom 150 cm, respectively; 68 samples (5 g/sample) were freeze-dried and homogenized. Sediments were solvent-extracted with 9:1 v/v dichloromethane (DCM) to methanol (MeOH) at 30°C for 30 min and 1500 psi for 5 min using an accelerated solvent extraction system. Total lipid extracts were separated using column chromatography through silica gels into the hexane and neutral fractions, using hexane and 9:1 v/v DCM: MeOH, respectively. The hexane and neutral fractions contain n-alkanes and the remaining parts of the lipid, such as long-chain ketone, respectively.
n-alkanes analysis
Compound identification and the molecular abundance distribution of long-chain n-alkanes were detected by gas chromatography (GC, Agilent 7890B) equipped with a flame ionization detector (FID) for separation and quantification. The FID temperature was 310°C, the H2 speed was 30 mL/min, the airflow speed was 400 mL/min, and each sample tested was 1 µL. Individual n-alkane identification was performed by comparing the retention times of n-alkanes in samples with the aforementioned standards. We report the average chain length (ACL) (Poynter et al., 1989), carbon preference index (CPI) (Bray and Evans, 1961; Ortiz et al., 2010), and Paq (Ficken et al., 2000):
Hydrogen and carbon isotope analysis
δD and δ13C values of long-chain n-alkanes were measured with a gas chromatographer (Trace GC Ultra) connected to a Thermo Fisher Delta V Advantage isotope ratio mass spectrometer. δD values were normalized to the Vienna Standard Mean Ocean Water (VSMOW) hydrogen isotopic scale. We measured the δ13C values of leaf-wax n-alkanes. Analytical conditions were the same as the abovementioned, and the δ13C values were normalized to the Vienna Pee Dee Belemnite (VPDB) carbon isotopic scale.
Results
n-alkane distributions
The carbon numbers in the peat samples ranged from C17 to C35. The relative abundances of n-alkanes in the peat samples show the dominance of mid- to long-chain n-alkanes with variable major peaks. Figure 3b and c show two distinct n-alkane distribution patterns. Pattern A (a single-peak type) features the highest peak at C29 or C31, and pattern B (a double-peak type) features the highest peak at C25, C27, C29, or C31. Pattern B has considerably more mid-chain n-alkanes than pattern A.

Results of n-alkanes and its parameters values. (a) n-alkane distributions and Paq, ACL, and CPI values of the XYH core. (b) and (c) n-alkane distribution patterns. P-A and P-B in (a) are the acronyms of pattern A and pattern B, respectively.
Figure 3a shows the variations in the n-alkane distributions and parameters; ACL, Paq, and CPI values. The ACL values generally fluctuated from 28.0 to 30.7 with a mean of 29.0, from 813 to 2010 AD. The values remained low from 813 to 1025 AD then rapidly increased to a mean of 29.5 at the end of 1125 AD. The mean ACL value fell to 28.9 from 1125 to 1450 AD before growing significantly and remaining stable for more than 200 years. From 1650 to 1925 AD, ACL values reached a minimum of 28.4 and then recovered slightly to 28.9. Paq values ranged from 0.11 to 0.40 with a mean of 0.24. The average Paq values had a relatively low mean of 0.19 from 813 to 1025 AD, increased to 0.26, from 1125 to 1440 AD, varied for a decade, and increased to a peak of approximately 0.32 between 1650 and 1925 AD. CPI values ranged from 1.24 to 6.01 with a mean of 2.99. The values remained high from 813 to 900 AD with a mean of 3.48, and then sharply decreased to 1.36 at approximately 1000 AD. After a rapid increase from 1025 to 1125 AD, the values decreased, with a mean of 2.01 from 1125 to 1250 AD and a mean of 2.26 from 1250 to 1450 AD. The CPI values then increased considerably from 1450 to 1800 AD with a mean of 3.68, followed by a decrease in CPI from 1800 to 1925 AD.
Hydrogen and carbon isotope compositions of n-alkanes
The δD values of C27, C29, and C31 n-alkanes (Figure 4a) ranged from −178.00‰ to −199.00‰, −179.00‰ to −200.95‰, and −182.00‰ to −233.00‰, respectively. The C27 n-alkanes exhibited less covariance with C29 (r2 = 0.2809, p = 0.011) and C31 n-alkanes (r2 = 0.175, p = 0.0124), whereas C29 and C31 n-alkanes exhibited similar covariances (r2 = 0.5225, p < 0.0001). The δDC29 values decreased from 813 to 1125 AD, increased by 9‰ from 1125 to 1300 AD and 1300 to 1650 AD, and thereafter decreased significantly from −189‰ to −199‰ and fluctuated at −199‰ for approximately one decade. The δDC29 values increased again and fluctuated at −195‰ after 1650 AD. The weighted mean values of δDC29 and δDC31 were calculated to minimize the influence from vegetation and ecosystem type variation. The δ13C values of C27, C29, and C31 n-alkanes (Figure 4b) ranged from −24.9‰ to −33.2‰, −25.9‰ to −33.3‰, and −26.9‰ to −33.9‰, respectively, and exhibited significant covariation over time. The weighted average δ13Cmean values of C27, C29, and C31 n-alkane δ13C values, which ranged from −26.0‰ to −33.5‰ with a mean value of −30.3‰, were presented. The lowest δ13Cmean values (−33.47‰) were reported at approximately 820 and 1700–1800 AD, and the largest δ13Cmean values (−26.01‰) occurred at approximately 1400 AD. A growth of 4.4‰ occurred from 813 to 1025 AD, followed by a rapid decrease. The δ13Cmean values gradually increased by 2.2‰ at 1440 AD and generally decreased by 4.28‰ from 1440 to 1650 AD, with several short intervals of increase, before slowly growing and stabilizing at approximately −32.9‰.

Results of hydrogen and carbon isotope composition of n-alkanes. (a) δD values of C27, C29, and C31 n-alkanes and δDmean values. (b) δ13C values of C27, C29, and C31 n-alkanes, and δ13Cmean values.
Discussion
Origin of sedimentary n-alkanes and proxies interpretation
Modern studies have illustrated that short-chain n-alkanes with 15–19 carbon atoms are primarily produced by algae and photosynthetic bacteria. Moss and aquatic macrophytes, including submerged, floating, and emerged plants, biosynthesize multiple mid-chain n-alkanes with 21–27 carbon atoms. Terrestrial plants predominantly produce odd-numbered long-chain n-alkanes with 25–37 carbon atoms (Barnes and Barnes, 1978; Cranwell, 1984; Cranwell et al., 1987; Eglinton and Hamilton, 1967). In this study, the proportions of mid- to- long-chain n-alkanes (C21–C31) accounted for approximately 95% of the total. Although the proportional n-alkanes input from diverse plants does not accurately indicate these amounts, terrestrial plants, moss, and aquatic plants with little input of microorganisms or diatoms are considered the dominant contributors to sedimentary n-alkanes. The present results indicate a prominent change in the proportion of aquatic and terrestrial plant inputs. Paq was defined to reconstruct the input of aquatic macrophytes in peat (Ficken et al., 2000). ACL23–33 values are the indicators of temperature, precipitation, and effective moisture depending on the dominant regional climatic control factors (Naafs et al., 2019). In this study, ACL23–33 values exhibited an opposite trend to that of Paq change, indicating that hydrological changes caused a decrease in aquatic macrophytes and an increase in terrestrial plants. The two n-alkane distribution patterns alternated through the sedimentary core (Figure 3a). During the time intervals of 900–1125 AD, 1450–1650 AD, and 1800–2000 AD, C29 and C31 n-alkanes were the most abundant homologues and exhibited pattern A one-peak mode. From 813 to 900 AD, 1125 to 1450 AD, and 1650 to 1800 AD, pattern B that represented two-peak variants had higher mid-chain n-alkanes C25 than pattern A. Pattern A represents the period in which terrestrial plants thrived and provided the principal sediment contribution. However, pattern B indicates the combined input of aquatic macrophytes and terrestrial plants. This inference is consistent with the variation in Paq. When n-alkanes were distributed as pattern A, the Paq values and C25 n-alkane concentrations were remarkably low, that is, from 1450 to 1650 AD. Pollen records in the same study core exhibited a conspicuous decline in wetland herbs and increase in Poaceae. When n-alkanes were distributed as pattern B, the Paq values were higher as the C29 and C31 n-alkane concentrations increased. During these periods, the abundance of wetland-herb pollen increased significantly, whereas that of Poaceae, birch, and alder pollen decreased. Therefore, the primary n-alkane sources were the peatland and surrounding hills, as recorded by regional paleoclimate changes (Cui et al., 2018).
The δ13C values of n-alkanes from sedimentary leaf wax lipids are typically used as the indicators of the relative contribution of the different photosynthetic pathways, C3 and C4. Almost all trees, cool-season grasses, and cool-season sedges use the C3 photosynthetic pathway, whereas the C4 photosynthetic pathway is used by warm-season grasses and sedges (Castañeda et al., 2009). Submerged and surface aquatic plants typically use the carbon in CO2 to maintain regular photosynthesis. C3 and C4 plants have different CO2 fixation modes, whereby, C3 plants incorporate CO2 with the Calvin cycle and show enriched 12C (depleted δ13C values range from −34‰ to −24‰). C4 plants fix CO2 with the Hatch–Slack cycle and exhibit enriched 13C (enriched δ13C values ranging from −19‰ to −6‰) (Lichtfouse et al., 1994). The δ13C values from the Beta lab showed a significant difference from the one from our n-alkanes. We noticed that most of the AMS δ13C values ranged from −18‰ to −22‰, but some of them have an obvious deviation (samples at the depth of 51, 86, 245, and 300 cm). According to Beta Analytic White Paper (www.radiocarbon.com/accelerator-mass-spectrometry.htm), the AMS δ13C values were measured on accelerator mass spectrometry with a series of pretreatments such as the combustion of the samples, the CO2 graphitization reaction, and the AMS deflection during the testing. It is widely accepted that isotopic fractionation during the processes occurs, which could lead to 0‰–20‰ anomaly from the real δ13C from the samples tested on gas chromatographer in our lab. In addition, the AMS δ13C values could be used to 14C isotopic correction and calculation, and cannot represent the condition of the original samples. The differences that occurred from the table might come from the different materials’ pretreatment and the AMS testing processes.
Relative humidity and vegetation type influence leaf-wax δD. δD values in vascular plant leaf waxes, such as C29 n-alkanes, have been used to identify the δDprecipitation source water (in the present study area, the source water is primarily precipitation) values during growth periods (Sachse et al., 2006). Nevertheless, other factors, such as leaf-water evaporation and species distribution, can affect leaf-wax δD. Previous studies have shown that εwax/p (the apparent fractionation between the leaf-wax n-alkane δD and precipitation δDprecipitation) is affected by multiple factors, such as the growing season length, plant life form, precipitation amount, and evapotranspiration (Sachse et al., 2012). However, ombrotrophic peatlands are singular in vegetation and environment in some ways, minimizing the influence of leaf-water evaporation and species (Nichols et al., 2009). Leaf water evaporation is primarily influenced by the leaf-to-air vapor difference, which depends on relative humidity and stomatal conductance (Hou et al., 2008). The relative humidity in ombrotrophic peatlands fluctuates between 75% and 85%, and therefore, the effect of variations in relative humidity on the δD of leaf wax in the peatland could be minimized. A previous study showed that the impact of relative humidity on plant leaf-wax δD was rather small in growth experiments, with only a 5.7‰ deviation for a relative humidity difference of 50% (Hou et al., 2008). The stratigraphic pollen record shows that the primary assemblages in the area are Pinus, Castanopsis, Poaceae, Artemisia, Cyperaceae, and Apiaceae (Cui et al., 2018); 40% of the pollen assemblage was dominated by upland herbs. Although regional vegetation changes substantially throughout the period, the site remains a peatland with relatively low diversity. Most of the vegetation changes in peatland are between vascular plants and the Sphagnum, and the sediment leaf wax originates from peatland plants in situ. The differences between species can be minimized. Herein, the δD of vascular plants was used to interpret the δDprecipitation. The δD of stored peatland water represents the mean annual precipitation and is more negative with abundant precipitation.
Aquatic plants are dominated by mid-chain n-alkanes, whereas terrestrial plants are dominated by long-chain n-alkanes. The C23/C31 n-alkane ratio was used to reconstruct the input of Sphagnum in modern peat (Nott et al., 2000). However, the C23/C31 n-alkane ratio can also be affected by environmental and biosynthetic processes. The Paq index indicates the relative contributions of aquatic macrophytes to peatlands and the changes in peat hydrology (Nichols et al., 2006; Routh et al., 2014; Zhou et al., 2005). High Paq values (>0.7) indicate a significant contribution of aquatic macrophytes, whereas low Paq values (<0.1) indicate a contribution from terrestrial plants. However, Paq values cannot be applied to distinguish the similar lipid distributions of Sphagnum and aquatic macrophytes (Baas et al., 2000; Ficken et al., 2000), particularly when plant macrofossils are absent or poorly preserved (Ortiz et al., 2010). In this study, the Paq values complement the C23/C31 n-alkane ratio to reflect the Sphagnum (aquatic macrophyte) input relative to the terrestrial plants in peats. The increase in C23/C31 n-alkane ratio and Paq indicates a large amount of aquatic macrophytes, which typically flourish under wet and humid conditions, and vice versa.
The CPI values of n-alkanes indicate the predominance of odd-numbered over even-numbered n-alkanes. Although CPI from fresh plant is high (e.g. 4–40), it decreases with diagenesis and microbial input, reaching 1. The mean CPI value in this study was 2.99, which is most likely caused by the bacterial diagenesis of n-alkanes and potential production from bacterial long-chain n-alkanes. This biosynthesis process often occurs at elevated temperatures. In this circumstance, CPI values cannot be considered as a direct proxy to infer the changes in paleoclimate but can support other temperature records in the nearby area.
Reconstruction of the paleoclimate and vegetation
Based on the n-alkanes and stable isotope records from the XYH peat core, we divided the climate evolution into five stages (Figure 5).

Comparison of different climatic records. (a) Paq record (this study), (b) C23/C31 n-alkane ratio record (this study). (c) δDmean records (this study). (d) Wetland-herb pollen concentration of core XYH (Cui et al., 2018). (e) Arid index Iph phytolith records from XYH (Zhang et al., 2019). (f) Speleothem δ18O record from Cave Heshang (Hu et al., 2008). (g) Dry-Wet index in Jiang-Nan area, approximately 25°N–31°N (Jingyun et al., 2006). (h) δDC31 record from Poyang Lake (Yao et al., 2015). (i) Warm index Iw phytolith records from XYH (Zhang et al., 2019). (j) δ13Cmean record (this study). (k) Mean air temperature anomaly in East China over the past 3000 years (Liu, 2004 after Zhang et al., 2000).
Stage I (900–1125 AD) (red bar in Figures 5 and 6): The sediment was primarily dark peat with approximately 10 cm of silty mud. The sedimentation rates increased from 0.27 to 0.37 cm/year. The low Paq and C23/C31 ratios, compared to the later period, revealed a smaller amount of aquatic plant input, which was correlated with relatively less negative δDmean values. The overall decrease in δDmean during this period implied a wetting trend. The overall dry conditions weakened the growth of aquatic plants, which is supported by the relatively low pollen concentration of wetland herbs (Cui et al., 2018) and phytolith index (Iph) values (Zhang et al., 2019). The δ13Cmean values were relatively negative but increased, as C3 biomass was the dominant and C4 plants increased. Studies have shown significantly negative correlations between δ13C values and mean annual precipitation for C3 plants and a significantly positive correlation between δ13C values and mean annual temperature for C4 plants (Liu and An, 2020). From the compiled sedimentary and historical records of temperature in east China (Liu, 2004; Zhang et al., 2000), the temperature increased significantly compared to that from 813 to 900 AD. This increase in temperature corresponded with the high phytolith warm index Iw values and likely implied warming conditions. The δ13Cmean values increased as the number of C4 plants increased, and therefore, it is suggested that higher temperatures benefit C4 plant growth.

(a) Paq record (this study). (b) C23/C31 ratio record (this study). (c) δDmean record (this study). (d) δ13Cmean record (this study). (e) El Niño 3.4 index (Emile-Geay et al., 2013). (f) Solar irradiance record (Bard et al., 2000). (g) Sea surface temperature (SST) of the eastern equatorial Pacific (Conroy et al., 2010). (h) Historic population change records from the historical documents in the Library of Nanjing University.
Stage II (1125–1450 AD) (yellow bar in Figures 5 and 6): The sediment was primarily silty mud with lower sedimentation rates. The Paq and C23/C31 ratios increased markedly, which indicates an increased amount of aquatic plant input compared to the last stage. The n-alkane distribution exhibited a two-peak pattern, which supports the Paq and C23/C31 changes. The more negative δDmean values implied wetter conditions, which corresponded with the sharp increase in the aquatic type of phytolith (Zhang et al., 2019). The wet climate during this period may have favored the growth of aquatic plants. The finer sedimentary conditions of silty mud also indicated the wet climate (Li et al., 2017). The more negative δ13Cmean values from the first ten decades implied that the C3 plants initially grew rapidly, and continuously less negative δ13Cmean values demonstrated increasing C4 plants. The concentration of Cyperaceae pollen peaked during this humid period (Cui et al., 2018), indicating the growth of C4 plants. The temperature remained high and the warm period continued, where Iw values exhibited a similar trend (Zhang et al., 2019). Herein, it is assumed that this period featured warm and wet conditions. The wet condition did not restrict the growth of C3 and C4 plants; however, the long-lasting high temperature decreased C3 plant growth and contributed to C4 plant recovery. C4 plants gradually surpassed C3 plants before peaking.
Stage III (1440–1650 AD) (dark-blue bar in Figures 5 and 6): The sediment was dark peat with considerably lower sedimentation rates of approximately 0.17 cm/yr. The remarkably increasing δDmean values indicate a drier condition, and the low Paq and C23/C31 ratio suggest a decrease in aquatic plant input, which could also be reflected by the low pollen concentration of the wetland herbs (Cui et al., 2018). The one-peak mode of n-alkane distribution reveals that terrestrial higher plants were dominant. The δ13Cmean values peaked at approximately 1450 AD and then gradually decreased during this stage. C3 plants thrived and gradually dominated, whereas C4 plants suffered. The temperature records showed a significant cooling trend (Liu, 2004; Zhang et al., 2000), which was supported by the decreasing Iw values (Zhang et al., 2019). C4 plants are significantly influenced by temperature than by precipitation, and therefore, the cooling condition depressed the growth of C4 plants, although dry condition was suitable for C4 plants.
Stage IV (1650–1800 AD) (light-blue bar in Figures 5 and 6): The dark peat during this period was deposited at a higher rate. The markedly decreasing δDmean values reveal the significant transition to a wet condition. The increasing Paq and C23/C31 ratios indicate the growing amount of aquatic plant input, which was supported by the larger pollen percentage of wetland herbs (Cui et al., 2018). The significantly decreasing δ13Cmean values demonstrate a similar trend to that of Stage III, such as the continuous growth of C3 plants. According to the temperature records, after reaching the lowest anomaly point, the condition remained at a cool level, which was also indicated by the Iw values (Zhang et al., 2019). It is assumed that the low temperature weakened the C4 plants, whereas the wet conditions favored the thriving of C3 plants. C3 plants were dominant again after this period. A wet and cool environment results in coarse particles.
Stage V (1800–2010 AD) (green bar in Figures 5 and 6): The dark peat during this period accumulated at a slightly higher stable rate than Stage IV. The δDmean values remained negative at a level similar to that of Stage IV. Increasing Paq and C23/C31 ratios suggested a slight shift to wet conditions, with more aquatic plants in the peat. The lower concentrations of wetland herbs and Iph values exhibit a reverse trend (Cui et al., 2018; Zhang et al., 2019). This contradiction is likely caused by human activities because of the increasing population from 1800 AD. The wetland herbs on the hillsides might have been cleared for rice or tea cultivation, which would have probably decreased original wetland herbs. The δ13Cmean values increased slightly when the temperature records implied a shift to warm conditions. The influencing factors of C3 and C4 plants cannot be explained because human activities may have begun to influence the evolution of vegetation since approximately 1800 AD.
Combined with temperature records, the curves described herein represent five primary climatic states. These are the warming climate and overall dry condition from 900 to 1125 AD, the warm and wet climate from 1125 to 1450 AD, the cool and dry condition from 1450 to 1650 AD, the cool and wet climate from 1650 to 1800 AD, and the warm and wet climate from 1800 to 2010 AD. Previous studies have shown several climate intervals in the Northern Hemisphere, namely MWP (900–1300 AD), LIA (1350–1850 AD), and PWP (1900 AD to now). In the present study, n-alkane records showed that Stage I and Stage II from 900 to 1450 AD could be a response to MWP. The MWP could be divided into two periods, and the distributions of n-alkanes and δDmean values revealed a dry-to-wet transformation. Overall, the records exhibited a wetter MWP. Climate changes in Stages III and IV from 1450 to 1800 AD respond to the LIA period, which could also be divided into two sub-climatic conditions, with the former characterized by a drier and cooling condition, and the latter being a considerably wetter and cooling condition. After the LIA, the later part of Stage V beginning in 1800 AD corresponds to the PWP, with an overall warm and wet climate.
Compared with the LIA, the core sections of the MWP and the PWP contain more fine particles (Li et al., 2017), which also supports the warmer and wetter conditions than those of the LIA. The interval between 1450 and 1650 AD had a significantly lower sedimentation rate (Figure 2). Although the n-alkane proxies, indicating a cooler and drier climate during this period would constrain microbial degradation and facilitate the accumulation of organic matter, vegetation types could significantly affect the degradation. The decay rate of vascular peat-forming plants is generally higher than that of Sphagnum species (Schellekens et al., 2015). The C23/C31 ratio and Paq were low during 1450–1650 AD, indicating the low input of aquatic plants, and the less negative δ13Cmean values indicated the increased growth of C4 sedges. It is assumed that the vegetation type change is the primary cause of the remarkably low sedimentation rate of 1450–1650 AD.
Climatic intercomparison and possible controlling factors
Comparison of paleoclimate records
The changes in climate and vegetation from the current findings, particularly during the universal MWP, LIA, and PWP periods (Figure 5a–c and j), were recorded at other sites or historical records in the East Asia monsoon area. The XYH peatland δDmean record and C23/C31 ratio indicate centennial-scale fluctuations in the hydroclimate, with an overall wetter MWP, drier LIA, and wetter PWP (Figure 5c). There exist sub-stages in both the MWP and LIA, such as the drier-to-wetter transfer during the MWP and the drier-to-wetter transfer during the LIA. The D–W index (Figure 5g) (Jingyun et al., 2006) calculated from official historical records in the Jian-Nan area (approximately 25°N–31°N) suggested the same trend. However, there are differences, such as the degree of wetter conditions during the early LIA from 1450 to 1650 AD and the level of aridity during the later LIA from 1650 to 1800 AD. This discrepancy is acceptable for the deviation between historical and natural records. This is not the only hydrology reconstruction in the middle–lower Yangtze River area, and although it is a rare high-resolution δD reconstruction from ombrotrophic peatlands covering the MWP and LIA period, we compared the values to proxy archives in the region. The late MWP and late LIA trended toward more depleted speleothem δ18O from Heshang Cave, and therefore, the early MWP and LIA trended toward more enriched values of δ18O (Figure 5f) (Hu et al., 2008). Although there are some anomalies because of dating error and resolution effects, the C23/C31 ratio, δDmean, and δ18O generally exhibited a similar hydroclimate history. The Poyang Lake δDC31 of plant wax from lake sediments exhibited a discrepancy, where the lake δDC31 values (Figure 5h) indicated a wetter condition in the early LIA and a drier condition in the late LIA (Yao et al., 2015). As discussed earlier, δD from plant wax could be affected by the “amount effect,” which would be recorded in the sedimentary n-alkane δD values (Sachse et al., 2012). The study site is nearby Poyang Lake. However, peatland characteristics, such as high relative humidity and relatively singular plant group, deplete the “amount effect” and the effect of vegetation lifeforms, which may lead to the short time difference. In general, plant lipids from sedimentary accumulations tend to average and alleviate these effects (Sachse et al., 2012), and therefore, the robust δD records of peatland and lake cores reveal a wetter MWP and a drier LIA.
Possible controlling factors
Previous studies show that the δ18O values of stalagmites are synchronized with the intensity of the East Asian Summer Monsoon (EASM), which controls the amount of precipitation (An, 2000; Ding et al., 2008; Zhang et al., 2011; Zheng et al., 2017). The humidity variations illustrated by n-alkane proxies and δD correspond to δ18O, and therefore, δD could be an indicator of EASM intensity, which is also significantly influenced by the El Niño Southern Oscillation (ENSO) (An, 2000; Chen et al., 2010, 2015; Oppo et al., 2009;). The West Pacific Subtropical High moves toward the southwest causing abundant precipitation during El Niño years and moves toward the northeast during La Niña years. Cold sea surface Temperatures (SSTs) in the eastern equatorial Pacific (La Niña) are correlated with excessive precipitation in Southeast Asia, whereas a warm eastern tropical Pacific (El Niño) causes drought conditions (Chen et al., 2015; Xu et al., 2016). In this study, the East Equatorial Pacific (EEP) SST (Figure 6g) was low during the early MWP, then swiftly turned high, and later returned to low values. The n-alkane data identified progressively wetter El Niño conditions from the MWP and drier La Niña conditions at the end of the MWP. The cool and dry EEP conditions associated with La Niña since 1450 AD are shown by n-alkane data, and El Niño conditions since 1800 AD later caused a slightly wetter climate. The SST and n-alkanes data showed that El Niño gradually predominated since 1800 AD. The generally warm-humid MWP period corresponds to El Niño-like conditions (relatively higher yet cooling SST of the EEP), whereas the cool-drier LIA period corresponds to La Niña-like conditions (relatively colder SST of the EEP). The Niño 3.4 index emphasizes the El Niño-like conditions in the MWP and the PWP and conspicuous La Niña-like conditions in the LIA.
During the warm period, the C3 plant-dominated scenario had fewer graminoids and woody angiosperms, which typically have lower CPIs. Furthermore, the more productive bacterial diagenesis under warmer conditions could also lead to the lower CPIs (Naafs et al., 2019). These two factors lead to low CPIs. In contrast, during the cool period, the abundance of C4 plants brings more graminoids and woody angiosperm plants, and the less productive diagenesis leads to higher CPIs. The CPI curve indicates that the centennial warm MWP and PWP experienced partially stronger solar irradiance, whereas the cooler LIA experienced weaker irradiance.
The Mount Xishan study area was not affected by human activities until the Song Dynasty (960–1279 AD), following which time war caused increased migration (Xu, 2000). North–south differences in the environment (northern drought and humid south) encouraged invaders from Mongolia to move toward southern China and compete for resources. Migration contributed to the population of Jiang Xi Province (Nan Chang County Annals) and lasted until the Yuan Dynasty (1271–1368 AD) when the southern China experienced a warm period (Yan et al., 2015). The warm and humid climate in the middle Yangtze River catchment encouraged population development. Although the stable society and agriculture-oriented policies supported by the Yuan Dynasty were suitable for development, the population decreased substantially, likely because of the cool and dry climate during the LIA. This trend may have led to the fall of the Ming Dynasty (1368–1644 AD). After the LIA, the warmer and more humid climate enabled a significant increase in population after 1750 AD.
Conclusions
In this study, a biomarker approach was applied to investigate the past 1200 years of climate variability from a peat core in the XYH peatland in southeast China. The input changes of the aquatic plants from n-alkane distribution patterns and the parameters (Paq and C23/C31 ratios) were reconstructed. Furthermore, the paleohydrological changes based on the δD from plant wax with the n-alkane changes and inferred vegetation changes from δ13C were qualitatively reconstructed. During the MWP (900–1450 AD), Paq, C23/C31 ratios, and δD generally inferred that the hydroclimate in the middle Yangtze River is wet and C4 vegetation increased under the dominance of C3 vegetation. In contrast, during the LIA (1450–1800 AD), a shift to relatively drier conditions occurred, and C3 vegetation generally increased. In the PWP (1800–2010 AD), the condition again became wetter compared to the LIA. This study also revealed that there are two intervals between the MWP and the LIA. The results are generally consistent with previous paleoclimatic reconstructions nearby and in situ but exhibit different hydroclimatic patterns to those from the southwest and northeast China in the monsoon areas (Chen et al., 2015).
In this study, it is proposed that the changes in ENSO and the EASM may have played key roles in producing hydroclimatic contrasts during the MWP and the LIA. The drier LIA was likely caused by the weaker EASM and La Niña conditions, whereas the wetter MWP and PWP resulted from the opposite conditions. Climatic conditions play an important role in social development. It is suggested that the up-and-down trend of the population during the MWP and the LIA was affected by the different water-heat conditions. With the development of instruments and fertilizers, the impact of human activities might have affected the local vegetation and microclimate from 1800 AD. However, further research is required to confirm this hypothesis
Research Data
sj-xls-1-hol-10.1177_09596836211025966 – Research Data for Biogeochemical evidence for environmental and vegetation changes in peatlands from the middle Yangtze river catchment during the medieval warm period and little ice Age
Research Data, sj-xls-1-hol-10.1177_09596836211025966 for Biogeochemical evidence for environmental and vegetation changes in peatlands from the middle Yangtze river catchment during the medieval warm period and little ice Age by Jia Sun, Chunmei Ma, Bin Zhou, Jiawei Jiang and Cheng Zhao in The Holocene
Footnotes
Acknowledgements
We thank Chaohao Ling and Guangjiu Ling for their assistance with field sampling work and the Nanjing Institute of Geography & Limnology, Chinese Academy of Sciences for the isotope test.
Availability of data and material
Datasets archiving for this research is underway. We have uploaded our data on Pangaea.de, yet they are being checked and processed. A copy of this data could be found in the “Supporting Information” file.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was supported by the National Natural Science Foundation of China (grant numbers 41371202, 41671196, and 41977378), Major Social Science Project (20&ZD247), and the Global Change Program of the National Key Research and Development Program of China (grant number 2016 YFA0600501).
References
Supplementary Material
Please find the following supplemental material available below.
For Open Access articles published under a Creative Commons License, all supplemental material carries the same license as the article it is associated with.
For non-Open Access articles published, all supplemental material carries a non-exclusive license, and permission requests for re-use of supplemental material or any part of supplemental material shall be sent directly to the copyright owner as specified in the copyright notice associated with the article.
