East Asian summer monsoon precipitation variability since the last deglaciation.

The lack of a precisely-dated, unequivocal climate proxy from northern China, where precipitation variability is traditionally considered as an East Asian summer monsoon (EASM) indicator, impedes our understanding of the behaviour and dynamics of the EASM. Here we present a well-dated, pollen-based, ~20-yr-resolution quantitative precipitation reconstruction (derived using a transfer function) from an alpine lake in North China, which provides for the first time a direct record of EASM evolution since 14.7 ka (ka = thousands of years before present, where the "present" is defined as the year AD 1950). Our record reveals a gradually intensifying monsoon from 14.7-7.0 ka, a maximum monsoon (30% higher precipitation than present) from ~7.8-5.3 ka, and a rapid decline since ~3.3 ka. These insolation-driven EASM trends were punctuated by two millennial-scale weakening events which occurred synchronously to the cold Younger Dryas and at ~9.5-8.5 ka, and by two centennial-scale intervals of enhanced (weakened) monsoon during the Medieval Warm Period (Little Ice Age). Our precipitation reconstruction, consistent with temperature changes but quite different from the prevailing view of EASM evolution, points to strong internal feedback processes driving the EASM, and may aid our understanding of future monsoon behaviour under ongoing anthropogenic climate change.

Whether the Asian summer monsoon responds directly on an orbital time scale to Northern Hemisphere summer insolation forcing without phase lag12, or whether the response is significantly delayed by internal feedback forcings345, has been intensively debated. This is mainly due to the lack of reliable palaeoclimatic records with a combination of a robust chronology and an unequivocal monsoon proxy678. As a major component of the Asian summer monsoon system, EASM variability controls almost all aspects of hydrology and ecology of East Asia910. Through its anomalous precipitation behaviour causing severe floods or droughts, the EASM considerably influences the economic and societal activities of one-third of the world’s population1112. In the last two decades, numerous studies have focused on the evolution of the EASM during the late Quaternary by investigating diverse natural archives such as aeolian sequences1314, lake sediments15, peat deposits16, and cave deposits21718. Among these studies, oxygen-isotope (δ18O) records from speleothems have been viewed as one of the most robust EASM records mainly because of their very precise age controls1718. However, the interpretation of speleothem δ18O in China remains highly controversial. Recent observational analysis67 and isotope modelling studies81920 suggested that speleothem δ18O in southern China does not reflect the local rainfall variability in the cave region of southern China. Rather, it is essentially a signal of precipitation isotopic composition, which is determined mainly by rainfall variability in the upstream region over the Indian Ocean and Indian Monsoon region21. Therefore, high-resolution reconstructions combining both a robust chronology and an unambiguous proxy which can help resolve the discrepancies between the different views of EASM evolution during the late Quaternary222 are urgently needed. In particular, the timing of the Holocene monsoon maximum, fundamental to revealing the phase relationship between the summer monsoon and insolation, is very contentious within the various reconstructions15172324. This impedes our full understanding of EASM behaviour and dynamics and limits our ability to predict EASM variability in the future, especially under a global warming scenario.

The EASM intensity can be directly represented by precipitation in North China: stronger (weaker) EASM circulation carries more (less) water vapour from the tropical Pacific and Indian Oceans, resulting in higher (lower) precipitation over North China2526. Here we use a high-resolution, pollen-based quantitative precipitation reconstruction from a well-dated sediment core retrieved from an alpine lake in North China to assess EASM variability and to characterize the underlying dynamical mechanisms since the last deglaciation, during which major drivers of Earth’s climate system (e.g. insolation27, Atlantic Meridional Overturning Circulation (AMOC)2829, El Niño/Southern Oscillation (ENSO)30) experienced dramatic changes from glacial to inter-glacial conditions.

Gonghai Lake (38°54' N, 112°14' E; 1,860 m above mean sea level), a hydrologically-closed alpine lake, is located on the NE margin of the Chinese Loess Plateau (Fig. 1), a typical EASM dominated region13142324. This montane area is well-suited for documenting natural climatic variability given the relatively limited human impact31, in contrast to the lowland plains and river valleys which have a long history of human activity. The current climate is dominated by the EASM with on average 77% of the annual 445 mm of precipitation occurring from June to September (Supplementary Fig. S1). Instrumental data indicate that precipitation at the site is representative of much of North China (Supplementary Fig. S2), and thus the site is ideal for assessing EASM variability. The regional vegetation is a transitional forest-steppe whose ecological variability is highly sensitive to precipitation changes, and past vegetation changes can reliably be reconstructed using pollen analysis. Our investigation of pollen assemblages and modern vegetation distribution suggests that the fossil pollen in the Gonghai Lake sediments mainly represents the vegetation composition of the area surrounding the lake, with only a minor contribution (<10%) from vegetation farther than 5 km from the lake (Supplementary Text S1). Given both the dominance of regional zonal vegetation (very few intra-zonal taxa) and the significant altitudinal zonation around the lake, the stratigraphic pollen data can be used to reflect regional vegetation succession and climate changes. Furthermore, results from a statistical approach for evaluating biotic-assemblage-based quantitative reconstructions32 demonstrate that the pollen record of Gonghai Lake is highly suitable for precipitation reconstruction (Supplementary Fig. S8).

Figure 1

Summer (June-July-August, JJA) mean 700 hPa streamline based on NCEP/NCAR Reanalysis during 1971–2000.

Blue dot indicates the location of Gonghai Lake. ‘EASM’, ‘ISM’, and ‘Westerlies’ denote the regions mainly influenced by the East Asian Summer Monsoon, the Indian Summer Monsoon, and the Westerlies, respectively. The modern Asian summer monsoon limit is shown by a green dashed line. Areas above 3000 m (above sea level) are shaded in grey. The Chinese Loess Plateau is shaded in light yellow. Daihai, Qinghai and Tianchi Lakes (red dots) and Dongge Cave (red cross) are additional key sites mentioned in the text. The figure was generated using GrADS v1.5.1.1269.

We obtained a sediment core (GH09B) from the central part of Gonghai Lake and focus here on the upper 9.42-m which consist of fine-grained lacustrine sediments (see Methods and Supplementary Fig. S3). The age model is based on a combination of 210Pb/137Cs dating of the uppermost 0.35 m (Supplementary Fig. S4) and twenty-five accelerator mass spectrometry (AMS) 14C dates of terrestrial plant macrofossils (Supplementary Table S1 and Supplementary Fig. S5), all converted to calendar years using the calibration of Reimer et al.33 (see Methods and Supplementary Table S1). The Bayesian age-depth model indicates that sedimentation was continuous, with a mean sediment accumulation rate of 64.3 cm/ka spanning the last ca. 14.7 ka. We also validated the stratigraphy and chronology of core GH09B by comparing the lithology and dating results of parallel core GH09C from the same location (Supplementary Fig. S3; also see Methods).

Core GH09B was sub-sampled at a 1-cm interval (20-yr resolution on average) and analysed for pollen, magnetic susceptibility (MS) and calcium carbonate (CaCO3) content using standard procedures (see Methods). Pollen analysis was employed to document vegetation changes and quantitatively to reconstruct precipitation variability. The MS values of the Gonghai Lake sediments mainly reflect catchment surface stability generally controlled by regional environmental conditions34. CaCO3 content is considered as an indicator of changes in water column chemistry with respect to lake status, given the dearth of carbonate minerals in the Gonghai Lake catchment.

Pollen assemblages in core GH09B are dominated by temperate-forest (Pinus, Picea, Betula, Quercus and Ulmus) and temperate-steppe (Artemisia, Chenopodiaceae and Poaceae) pollen types (Supplementary Fig. S6 and Supplementary Text S2). The abundance of arboreal and herbaceous taxa exhibits large (~60–80%) and abrupt changes over the last 14.7 ka (Fig. 2a,b). The percentages of tree pollen were low and herb pollen percentages high during the last deglacial period, and vice versa during the early-mid Holocene. Broad-leaved tree pollen reached maximum abundance from ~7.8–5.3 ka (Fig. 2c). After 3.3 ka, tree pollen decreased continuously, whereas herb pollen increased. Using a calibration-function method35, the pollen data were transferred quantitatively into annual precipitation (PANN). For this purpose, we used a modern surface pollen dataset consisting of 509 samples distributed along a large precipitation gradient (PANN: 0–1000 mm; PANN at nearby Ningwu Station is 445 mm) and a small mean annual temperature gradient (TANN: 3.6–8.6°C; TANN at Ningwu station is 6.1°C) to ensure reliable estimates of the precipitation optima and tolerances of the major pollen taxa, and hence a robust PANN reconstruction. In addition to statistical validation of the calibration function (the two component weighted-averaging partial least square (WAPLS) model featuring a coefficient of determination of 0.84 and a root mean square error of prediction (RMSEP) of 86 mm, see Methods and Supplementary Table S2), the close similarity of the reconstructed PANN to the instrumental PANN record for the most recent past (AD 1962–2008) (Supplementary Fig. S9) further demonstrates the robustness of our PANN reconstruction.

Figure 2

Variation in the reconstructed annual precipitation and other environmental proxies since 14.7 ka from Gonghai Lake, North China.

(a) Percentage of tree pollen (%). (b) Percentage of herb pollen (%). (c) Percentage of broad leaf tree pollen (%). (d) Annual precipitation reconstructed by a calibration function applied to the fossil pollen assemblages. (e) Magnetic susceptibility. (f) Carbonate content. Shaded area from ~11.5–3.3 ka denotes an interval of stable lake status and maximum vegetable cover in the Gonghai Lake catchment.

The time series of the pollen-based PANN portrays the evolution of the EASM since the last deglaciation (Fig. 2d). Compared to modern instrumental PANN (445 mm), a gradual 50% increase in precipitation from 380 mm at 14.7 ka to 570 mm at 7.0 ka indicates a continuously intensifying EASM across the late deglacial-Holocene transition. The maximum monsoon intensity, with an average PANN of 574 mm, ~30% higher than the modern value, occurred from 7.8 to 5.3 ka. A two-step gradual decrease in monsoon intensity commenced at 5.3 ka and 3.3 ka, with a more prominent decline at 3.3 ka (Supplementary Fig. S12) which continued until the end of the Little Ice Age (LIA). These long-term monsoon trends inferred from our PANN model are generally consistent with the MS and CaCO3 records (Fig. 2e,f), which suggest a stable lake and surrounding catchment with a high lake level during the early and mid-Holocene (~11.5–3.3 ka).

Our PANN reconstruction documents the major sub-orbital scale global climatic events since the last deglaciation [Fig. 3g (blue curve)]. EASM precipitation generally decreased during cold climatic episodes, and increased during warm events. For example, following the warm Bølling/Allerød interstadial, a ~100 mm decrease in precipitation from ca. 12.8–11.7 ka indicates a weakened EASM [Fig. 3g (blue curve)] during the Younger Dryas (YD) cold event within the limits of the available chronology, when the Northern Hemisphere temperature was significantly reduced36 [Fig. 3e (magenta curve)]. A further precipitation decrease from 9.5 to 8.5 ka in the early Holocene corresponds to the widespread cold episode inferred from marine proxies associated with weakened thermohaline circulation37. The weakest Holocene EASM occurred during the LIA, when the monsoon precipitation was about 100 mm less (~20%) than during the Medieval Warm Period (MWP) [Fig. 3g (blue curve)]. The increased precipitation during the last ~150 years [Fig. 3g (blue curve)], since the end of the LIA, which is evident in our reconstructed precipitation record and is supported by a recent tree-ring precipitation reconstruction38, is consistent with the Northern Hemisphere warming trends inferred from proxy-based surface temperatures39. The synchroneity of these EASM responses and global climatic events further confirms the sensitivity of Gonghai Lake records to EASM variability in North China.

Figure 3

Comparison of Gonghai Lake EASM record with various other regional and global environmental signals.

(a) Dongge cave speleothem δ18O records1754. (b) Tree pollen percentages from Qinghai Lake41. (c) Tree pollen percentages from Daihai Lake15. (d) Frequency distribution of Chinese Loess Plateau palaeosol dates40. (e) Synthesized Northern Hemisphere (30°–90°N) temperature record during the last deglaciation36 (magenta line) and Holocene55 (red line). (f) Western subtropical Atlantic 231Pa/230Th record28 (blue circles) and northeast Atlantic sortable silt record29 (cyan triangles), both of which may indicate AMOC strength, and synthesized meltwater flux in the Northern Hemisphere26 (blue line) and the rate of sea-level rise from the Laurentide Ice Sheet44 (purple line), demonstrating continuous freshwater input during the last deglaciation and early Holocene, with the intervals of rapid melting during the YD and from 9.5 to 8.5 ka. (g) Pollen-based annual precipitation (PANN) reconstructed from Gonghai Lake (blue line, this study) and similar reconstruction for the past 6 ka from the nearby Tianchi Lake42 (green curve) together with 65°N summer insolation27 (orange line). Black dots at the bottom are twenty-five AMS 14C dates from terrestrial-plant macrofossils with an uncertainty interval of 1σ. The shaded blue bars indicate periods of significantly decreased precipitation. The dashed lines indicate changes in the factors forcing EASM variability.

Previous proxy records from northern China may be unable fully to constrain the timing and amplitude of EASM variability in a detailed and precise way, given that many of them suffer from problems of limited temporal resolution, possible dating errors or ambiguous proxies47824. However, those records can be used to carry out general comparisons with our PANN reconstruction on a multi-millennial timescale. A stronger EASM from ~8 to 3 ka, as indicated by the highest and the lowest contents of tree and herb pollen, respectively (Fig. 2a,b), with around 550 mm of reconstructed PANN [Fig. 3g (blue curve)], is supported by the increased frequency of palaeosol development from ~8.6–3.2 ka in the Chinese Loess Plateau (Fig. 3d), and decreased aeolian-sand activity from ~8.6–3.2 ka in the four main sandlands located to the west and NE of Gonghai Lake40. In addition, a vegetation record from Daihai Lake (~200 km to the north of Gonghai Lake; Fig. 1) clearly demonstrates that the wettest climate and strongest summer monsoon occurred during the mid-Holocene15 (Fig. 3c). Further evidence for a mid-Holocene EASM maximum comes from Qinghai Lake, ~1000 km to the west of Gonghai Lake (Fig. 1), where pollen data suggest a mid-Holocene climatic optimum41 (Fig. 3b). A generally decreasing trend of precipitation since ~6 ka is also documented in Tianchi Lake on the Chinese Loess Plateau42 [Fig. 3g (green curve)]. Furthermore, the mid-Holocene maximum EASM inferred from our reconstruction from Gonghai Lake corresponds to significant Neolithic cultural developments in northern China from ~7–5 ka, the period when the Yangshao Culture attained maximum prosperity compared to other Neolithic cultures before and after (Supplementary Fig. S10 and Supplementary Text S3). This suggests that the high precipitation [~160 mm higher vs. modern value; Fig. 3g (blue curve)] associated with an intense EASM during the mid-Holocene may have provided a favourable environment for cultural development in semi-arid North China. These various lines of evidence support the contention that our quantitative precipitation reconstruction indeed represents precipitation change over a large area of northern China, and that it reliably documents the Holocene EASM evolution with maximum monsoon precipitation occurring in the mid-Holocene.

The foregoing highlights the complexity of monsoon dynamics since the last deglaciation, because the EASM was not driven simply by changing insolation, as inferred from speleothem isotope records2. Although the EASM intensity tracks the insolation curve, the deviation between them during the early-Holocene suggests a ~4 ka delayed response to the summer insolation peak at ~10 ka, which we attribute to the remnant melting Laurentide ice sheet which delivered continuous freshwater input to the North Atlantic until ~7 ka26434445. The linkage between North Atlantic cold events and a weak Asian summer monsoon on centennial-millennial time scales since the last glacial has been suggested by various proxy records174647. As our simulation results demonstrate, an anomalous freshening of the North Atlantic and the resultant weakened AMOC could have strongly depressed the insolation-driven EASM (Fig. 4a) through the strengthening of westerly and northerly winds (Supplementary Fig. S11). This would have caused the EASM to be weaker during the early Holocene than during the mid-Holocene (Fig. 4b), but still to be stronger than at present (Fig. 4c). Such a mechanical connection is also supported by present-day climate simulations showing that a reduction in the AMOC causes northerly surface wind anomalies over East Asia48. It is noteworthy that two intervals of significantly weakened millennial-scale EASM precipitation [Fig. 3g (blue curve)], during the YD and from 9.5–8.5 ka, correspond well to a weakened AMOC2829.

Figure 4

Summer precipitation changes (in mm/day) in sensitivity experiments (refer to Table S4 for detailed information).

The results show the weakening of the EASM at 8.5 ka in northern China when considering glacial boundary conditions (a). The weakened EASM at 8.5 ka was weaker than at 6 ka (b) but stronger than the present (c). Blue dot indicates the location of Gonghai Lake. Dotted areas indicate changes which are significant at the 95% confidence level.

After around 7 ka, the EASM strength changed broadly in response to orbitally-induced insolation forcing27 [Fig. 3g (orange curve)] because the glacial boundary conditions can be regarded as similar to the modern state, given both the dramatically decelerated rates of eustatic sea-level rise and the generally stable global ice volume after that time4950. However, different from the gradually decreasing insolation during the late Holocene, the EASM experienced an abrupt decline at 3.3 ka [Fig. 3g (blue curve) and Supplementary Fig. S12], which may be partly attributed to the influence of tropical ocean conditions. During El Niño events, the tropical eastern Pacific Ocean warming increases tropospheric temperature and induces an eastward propagating Kelvin wave, leading to an increase in the tropospheric temperature in the western Pacific Ocean5152, which weakens the land–sea thermal contrast and thus the EASM53. Both a paleo-ENSO reconstruction30 and our simulation indicate that the late Holocene exhibited an increased frequency of El Niño events (Supplementary Fig. S13), favouring a weak EASM (r = −0.50).

Our demonstration of a ~4 ka delay in the response of the maximum monsoon intensity to the Northern Hemisphere summer insolation maximum indicates that the prevailing view of an early Holocene EASM maximum inferred from speleothem oxygen isotope records in southern China1754 (Fig. 3a) should be significantly revised. Our comparative study of the EASM behaviour in the context of coeval North Atlantic and tropical ocean climate records indicates a complex relationship between the EASM, the external forcing factor (Northern Hemisphere summer insolation), and internal feedback mechanisms which involve both high-latitude (AMOC) and low-latitude (ENSO) processes. We anticipate that our findings may improve the predictive ability of future model-based precipitation simulations. In addition, the overall positive relationship demonstrated here between temperature3655 (Fig. 3e) and EASM precipitation [Fig. 3g (blue curve)], both on orbital and suborbital timescales, has significant implications for EASM behaviour under ongoing global warming, although the different responses of EASM to natural and anthropogenic forcing56 warrant further study.

Methods

Study area

Gonghai Lake (38°54′N, 112°14′E; 1860 m above mean sea level) is located in Ningwu County, Shanxi Province, on the northern margin of the Chinese Loess Plateau (Fig. 1). It is a freshwater alpine lake formed on a plateau of the watershed between the Sanggan and Fenhe rivers. The basin is hydrologically-closed and the water source is mainly from precipitation. The surface area is ~0.36 km2 and the maximum water depth is ~10 m; the lake floor has a flat topography. The area is located in a transitional zone between semi-arid and semi-humid conditions, on the fringe of the modern Asian summer monsoon. The mean annual precipitation in Ningwu is 445 mm, of which, 77% occurs from June to September. The regional vegetation in the surrounding Lvliang Mountains is dominated by Larix principis-rupprechtii, Pinus tabulaeformis and Populus davidiana forest, while on the plateau, Hippophae rhamnoides scrub, Bothriochloa ischaemum grassland and Carex spp. are widely distributed. Documentary evidence suggests that intensive deforestation in this region did not occur until ~600 cal yr BP (Ming Dynasty)31. The exposed bedrock of the zonal ground surfaces was mainly formed during the Archaean to Cenozoic, and is comprised mainly of sandstone.

Sediment lithology

Core GH09B extended to impenetrable substrate at a depth of 9.80 m. The interval 9.80–9.42 m consists of sand and gravel. The upper 9.42 m consists of lake sediment which can be divided into three lithological units: Unit 1 (9.42–7.70 m) consists of alternating layers of dark, silty clay and light–dark silt with the occasional occurrence of yellow sand. Unit 2 (7.70–2.60 m) consists of dark organic mud with light grey mud between 4.70 m and 2.60 m. Unit 3 (2.60–0.00 m) is dominated by silty clay with occasional plant macrofossils.

Analytical methods

The bulk magnetic susceptibility (MS) of cores GH09B and GH09C was measured at a 1-cm interval using a Bartington MS2 magnetic susceptibility meter and MS2B sensor (Supplementary Figs. S3D and S3E). 771 samples, at ~1-cm interval, were analysed for pollen. ~1 g of sediment was treated using a modified HCl-NaOH-HF procedure57. One tablet of Lycopodium spores was added to each sample in order to estimate the pollen concentration. More than 500 terrestrial pollen grains were counted for each sample, and 81 pollen taxa were identified, including aquatic pollen grains, and spores. The sum of total terrestrial pollen was used as the denominator when calculating pollen percentages. Identifications are based on the Pollen Flora of China58 assisted by modern reference collections from Hebei Normal University. Pollen diagrams were plotted using Tilia v2.0.b.459 and pollen-assemblage zones were constructed using stratigraphically-constrained cluster analysis (CONISS)60.

Chronology

The chronology of the top 35 cm of core GH09B was obtained using a combination of 210Pb and 137Cs dating (Supplementary Fig. S4). Eighteen terrestrial plant macrofossil samples from core GH09B and seven terrestrial plant-macrofossil samples from core GH09C were measured for radiocarbon (14C) at the Xi’an AMS Laboratory, China, and at Beta Analytic Inc., Florida, USA. The MS records of these two cores correlate well and therefore the 14C sample depths of core GH09C were transferred to the corresponding depths of core GH09B, using magnetic susceptibility correlations. Thus we developed for the first time an exceptionally high-resolution age model for northern Chinese lakes. All conventional 14C ages were calibrated to calendar years using the IntCal09 calibration curve33. Supplementary Table S1 lists all measured and calibrated ages. The ages are expressed in years before present (BP) where “present” is defined as AD 1950. Bayesian age-depth modelling for Core GH09B (Supplementary Fig. S5) was performed using OxCal v4.2.2 and a Poisson-process (P-sequence) single depositional model at 1-cm increment with a K value of 10061.

Quantitative reconstruction of precipitation

A modern calibration set of 1,860 surface pollen samples was used after excluding samples from tropical and arid regions, and those from sites potentially influenced by human activity, from the total of 2,689 samples62636465. The reconstruction was refined by limiting the mean annual temperature (TANN) range to obtain a more robust pollen-based PANN calibration function, given that temperature and precipitation are the most significant variables controlling regional vegetation and that the target variable of this reconstruction is precipitation. A subset of 509 samples (Supplementary Fig. S7) within the TANN range 3.6–8.6 °C, which were ±2.5 °C from TANN of the nearest Ningwu Station, was used. The range of PANN is relatively large, from 0 to 1,000 mm. Long compositional gradient lengths (axis 1 = 3.4 standard deviation (SD) units of compositional turnover, axis 2 = 3.6 SD units) in the pollen data set suggested a unimodal-based method for regression and calibration66. A two component WAPLS model exhibiting the strongest statistical performance (r2 = 0.84, RMSEP = 86 mm) (Supplementary Table S2) was chosen for quantitative precipitation reconstruction66. In addition, a novel method was employed to assess the statistical significance of the reconstruction using both the modern and fossil data32. This showed that the PANN reconstruction based on the WAPLS-2 model and stratigraphical pollen data from core GH09B explains more of the variation in the fossil data than 99% of reconstructions (p = 0.001) derived from calibration-functions trained on random environmental data, whereas the TANN reconstruction is not statistically significant (p = 0.434) (Supplementary Fig. S8). In addition, this method was employed to reconstruct PANN changes using pollen percentage data from Tianchi Lake over the past 6 ka42, plotted in Fig. 3g (green curve). The best-performing two-component WAPLS model based on 525 modern samples has an r2 of 0.87 and a RMSEP of 74.7 mm [99% confidence level (p = 0.01)].

EASM simulation experiments

The Kiel Climate Model (KCM), a coupled atmosphere–ocean–sea-ice general circulation model (AOGCM), was used to model the EASM from 9.5 ka to 0 ka with 10-times acceleration, as detailed in Lorenz and Lohmann67. Prior to the transient simulation (HT), two 1,000-yr control simulations for the pre-Industrial period (AD 1800, H0K) and the early Holocene (9.5 ka, H9K) were performed. The HT simulation was then started at the end of the H9K equilibrium run with transient orbital parameters. Supplementary Table S3 summarizes the boundary conditions adopted for these three simulations. A series of sensitivity simulations (hereafter PI, Exp_6 ka, Exp_8.5 ka, Exp_8.5 kaICE, and Exp_8.5 kaMELTICE) were performed following the description in Jin et al.68 with the Community Climate System Model 3 (CCSM3), which is also an AOGCM. PI, Exp_6 ka and Exp_8.5 ka are control runs for AD 1800, 6 ka and 8.5 ka, with corresponding orbital parameters and greenhouse gas concentrations. Additional forces were implanted in Exp_8.5 kaICE, and Exp_8.5 kaMELTICE. Supplementary Table S4 summarizes the boundary conditions adopted for these five simulations.

Additional Information

How to cite this article: Chen, F. et al. East Asian summer monsoon precipitation variability since the last deglaciation. Sci. Rep. 5, 11186; doi: 10.1038/srep11186 (2015).