Hostname: page-component-7b9c58cd5d-9klzr Total loading time: 0 Render date: 2025-03-15T19:46:12.958Z Has data issue: false hasContentIssue false

East Asian summer monsoon climates and cave hydrological cycles over Dansgaard-Oeschger events 14 to 11 revealed by a new stalagmite record from Hulu Cave

Published online by Cambridge University Press:  05 August 2019

Yijia Liang
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
Yongjin Wang*
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
Quan Wang
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China School of Geography and Tourism, Chongqing Normal University, Chongqing 401331, China
Jiangying Wu
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
Qingfeng Shao
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
Zhenqiu Zhang
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China School of Life Sciences, Nanjing Normal University, Nanjing 210023, China
Shaohua Yang
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
Xinggong Kong
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Key Laboratory of Virtual Geographic Environment, Nanjing Normal University, Ministry of Education, Nanjing 210023, China
R. Lawrence Edwards
Affiliation:
School of Geography, Nanjing Normal University, Nanjing 210023, China Department of Earth Sciences, University of Minnesota, Minnesota 55455, USA
*
*Corresponding author. e-mail address: yjwang@njnu.edu.cn (Y. Wang).
Rights & Permissions [Opens in a new window]

Abstract

A 230Th/U-dated stalagmite from Hulu Cave was analyzed for δ18O, δ13C, and trace elements. A ~10-yr-resolution δ18O record, spanning 51.7–42.6 ka, revealed Dansgaard-Oeschger (DO) events 14 to 11. A similar rapid transition and synchronous timing of the onset of DO 12 is evident between the Greenland and Hulu Cave records, which suggests a common forcing mechanism of DO cycles in the North Atlantic and monsoonal region of Asia. Centennial-scale monsoonal oscillations in the cave δ18O record are indicative of hydroclimatic instability during interstadials. After removing the signals of remote moisture sources, the proportion of moisture from nearby sources is found to be higher during stadials than during interstadials. To explain this, we propose that the movement of the westerly jet is an important control on the balance of nearby and distant moisture sources in East Asia. In addition, the records of δ13C and trace element ratios, which are proxies of local environmental changes, resemble the δ18O record on the scale of DO cycles, as well as on even shorter timescales. This suggests that hydrological processes and biological activity at the cave site respond sensitively to the monsoonal changes.

Type
Research Article
Copyright
Copyright © University of Washington. Published by Cambridge University Press, 2019 

INTRODUCTION

During the last glacial period, roughly 120–10 ka, high-amplitude, relatively short-duration climatic oscillations are evident in Greenland ice-core records (Johnsen et al., Reference Johnsen, Clausen, Dansgaard, Fuhrer, Gundestrup, Hammer and Iversen1992; Dansgaard et al., Reference Dansgaard, Johnsen, Clausen, Dahl-Jensen, Gundestrup, Hammer and Hvidberg1993; Andersen et al., Reference Andersen, Azuma, Barnola, Bigler, Biscaye, Caillon and Chappellaz2004). These so-called Dansgaard-Oeschger (DO) events were rapid, decadal-scale transitions from cold conditions (Greenland stadials, GS) to warm conditions (Greenland interstadials, GI), followed by a slow return to cold stadial conditions, within an interval of centuries to millennia (Wolff et al., Reference Wolff, Chappellaz, Blunier, Rasmussen and Svensson2010). In addition to Greenland ice cores, DO cycles have been observed in numerous geological archives in the Northern Hemisphere, including terrestrial and marine sediments (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001; Voelker, Reference Voelker2002; Deplazes et al., Reference Deplazes, Lückge, Peterson, Timmermann, Hamann, Hughen and Röhl2013; Rousseau et al., Reference Rousseau, Svensson, Bigler, Sima, Steffensen and Boers2017).

The character of DO events in Chinese cave deposits differs from, and is more complex than, paleoclimatic records from northern, high latitudes. An annually laminated stalagmite record from northern China resembles the Greenland NGRIP (the North Greenland Ice Core Project) ice-core record in terms of the timing and duration of the abrupt transitions into DO 15.2 and 14 (Duan et al., Reference Duan, Cheng, Tan and Edwards2016). The onset of DO 12, however, is extended in stalagmite records from two caves in southwestern China, lasting from one to several millennia (Cai et al., Reference Cai, An, Cheng, Edwards, Kelly, Liu and Wang2006; Han et al., Reference Han, Li, Cheng, Edwards, Shen, Li and Huang2016). These regional differences may result from somewhat different forcing mechanisms. In northern China, monsoon climates are dominated by temperature changes in high, northern latitudes mediated via the westerlies (Duan et al., Reference Duan, Cheng, Tan and Edwards2016). Climates in southwestern China, however, are controlled by cross-equatorial airflows and temperature changes in the Southern Hemisphere (Cai et al., Reference Cai, An, Cheng, Edwards, Kelly, Liu and Wang2006; An et al., Reference An, Wu, Li, Sun, Liu, Zhou and Cai2015). In addition, statistical analysis of a single paleoclimatic record from Hulu Cave in eastern China showed that both Northern and Southern Hemisphere climates affected the pattern of monsoonal DO events, with different ratios of southern and northern climate signals evident (Rohling et al., Reference Rohling, Liu, Roberts, Stanford, Rasmussen, Langen and Siddall2009). Indeed, paleoclimatic records from caves have been regarded as correlating either with temperature changes in Greenland (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001, Reference Wang, Cheng, Edwards, Kong, Shao, Chen and Wu2008; Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003; Zhao et al., Reference Zhao, Wang, Edwards, Cheng and Liu2010; Deplazes et al., Reference Deplazes, Lückge, Peterson, Timmermann, Hamann, Hughen and Röhl2013) or in Antarctica (Cai et al., Reference Cai, An, Cheng, Edwards, Kelly, Liu and Wang2006; Chen et al., Reference Chen, Wang, Hai, Edwards, Wang, Kong and Liu2016; Han et al., Reference Han, Li, Cheng, Edwards, Shen, Li and Huang2016) on the millennial scale. Thus, further analysis of speleothem records from Hulu Cave can potentially provide insights into the regional nature of DO signals, since the site was likely influenced in part by moisture from the Indian monsoon subsystem during the last ice age (Pausata et al., Reference Pausata, Battisti, Nisancioglu and Bitz2011). In addition, due to the sparsity of high-resolution speleothem records, little evidence is available to conduct a detailed analysis of monsoon climates on the centennial scale during Marine Oxygen Isotope Stage (MIS) 3.

Whereas the broad anti-correlation between rainfall somewhere in the monsoon system and speleothem δ18O from Chinese caves is generally accepted, the specifics are controversial. Several studies have interpreted such records as resulting from differential Rayleigh fractionation of water vapor between tropical moisture sources and cave sites and/or changes in the seasonal duration or intensity of monsoon rainfall (perhaps related to the seasonal position of the subtropical jet), a set of processes often referred to as “monsoon intensity” (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001; Cai et al., Reference Cai, An, Cheng, Edwards, Kelly, Liu and Wang2006, Reference Cai, Tan, Cheng, An, Edwards, Kelly and Kong2010; Kelly et al., Reference Kelly, Edwards, Cheng, Yuan, Cai, Zhang and Lin2006; Duan et al., Reference Duan, Wu, Wang, Edwards, Cheng, Kong and Zhang2015; Orland et al., Reference Orland, Edwards, Cheng, Kozdon, Cross and Valley2015; Tan et al., Reference Tan, Cai, An, Cheng, Shen, Gao and Edwards2017, Reference Tan, Cai, Cheng, Edwards, Gao, Xu and Zhang2018). The extent to which these processes correlate with rainfall amount has been confirmed in some localities through independent quantitative lake-level records of past rainfall variability (Zhang et al., Reference Zhang, Jia, Lai, Long and Yang2011; Goldsmith et al., Reference Goldsmith, Broecker, Xu, Polissar, Demenocal, Porat and Lan2017). A study of the 10Be record of Chinese loess deposits, a proxy of summer monsoon rainfall, suggested that orbital-scale monsoonal changes were also forced by differences in the proportions of moisture from Indian and East Asian monsoon sources (Beck et al., Reference Beck, Zhou, Li, Wu, White, Xian and Kong2018). This argument is consistent with some climate model simulations and interpretations based upon atmospheric reanalysis data (Maher, Reference Maher2008; Dayem et al., Reference Dayem, Molnar, Battisti and Roe2010; Wu et al., Reference Wu, Zhang, Li, Li and Huang2015). In the present study, we developed a new strategy for determining how the proportions of water vapor derived from remote and nearby sources affect speleothem δ18O records from China.

Proxies such as δ13C and trace element content of speleothems are effective indicators of palaeohydrological conditions under appropriate circumstances (Fairchild and Treble, Reference Fairchild and Treble2009 and references therein). However, compared to δ18O records, they have been underutilized in paleoenvironmental studies of speleothems in China. Nevertheless, in a specific cave environment, δ13C and trace elements can be applied as hydrological tracers at a local scale (Genty et al., Reference Genty, Blamart, Ouahdi, Gilmour, Baker, Jouzel and Van-Exter2003; Sinclair et al., Reference Sinclair, Banner, Taylor, Partin, Jenson, Mylroie and Goddard2012; Huang et al., Reference Huang, Wang, Cheng, Edwards, Shen, Liu and Shao2016; Zhao et al., Reference Zhao, Wang, Edwards, Cheng, Liu, Kong and Ning2016; Wang et al., Reference Wang, Wang, Shao, Liang, Zhang and Kong2018). Moreover, a relationship between δ13C and trace element ratios in speleothems has been observed; for example, owing to the impact of prior calcite precipitation (PCP) and CO2 degassing under low-flow conditions, metal/Ca ratios and δ13C are elevated in dripwater and hence in stalagmites (Treble et al., Reference Treble, Shelley and Chappell2003; Johnson et al., Reference Johnson, Hu, Belshaw and Henderson2006; Fairchild and Treble, Reference Fairchild and Treble2009; Stoll et al., Reference Stoll, Müller and Prieto2012; Chen and Li, Reference Chen and Li2018). Trace element enrichment can also be associated with warm and humid climatic conditions due to enhanced chemical weathering, and thus it reflects local hydrological dynamics (Borsato et al., Reference Borsato, Frisia, Fairchild, Somogyi and Susini2007; Zhou et al., Reference Zhou, Chi, Michael, Zhao, Yan, Alan and Feng2008). Hellstrom and McCulloch (Reference Hellstrom and McCulloch2000) showed that concentrations of Sr and Ba in speleothems were positively correlated with changes in the overlying vegetation cover, which confirms the effectiveness of carbon isotopes and trace elements for reconstructing local environmental changes. A recent study of a 3000-yr annually laminated stalagmite from Hulu Cave demonstrated the covariation of records of δ18O and Sr/Ca ratio on centennial to multi-decadal scales, which enabled the authors to interpret the δ18O record in terms of the monsoonal rainfall (Duan et al., Reference Duan, Wu, Wang, Edwards, Cheng, Kong and Zhang2015).

In this context, we conducted a multi-proxy study of a new 230Th/U-dated stalagmite from Hulu Cave (HL161), spanning the interval of 51.7–42.6 ka (relative to the present, defined as AD 1950). The δ18O record has an average resolution of ~10 yr, which enables us to identify fine-scale monsoonal variations during DO cycles 14–11. Our study has three aims: (1) to determine the timing and structure of the monsoonal response to DO events; (2) to reconstruct local hydrological processes; and (3) to provide a better understanding of the origin of speleothem δ18O records from caves in China.

SITE, MATERIALS AND METHODS

Stalagmite sample HL161 was found in a naturally detached state on the northern slope in Hulu Cave, Nanjing, in eastern China (32°30′N, 119°10′E, 86 m above sea level). A detailed site description is given in Wang et al. (Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001) and Duan et al. (Reference Duan, Wu, Wang, Edwards, Cheng, Kong and Zhang2015). The stalagmite, which has a length of 243 mm and an asymmetric shape, was sliced along its growth axis and polished; no depositional hiatuses are evident (Fig. 1).

Figure 1. (Color online) Chronology of stalagmite HL161. (a) Scanned image and location of fifteen 230Th samples. (b) The linearly-interpolated age model, including error bars for 230Th ages.

For 230Th/U dating, 15 subsamples of up to 100 mg were extracted using a 0.9-mm-diameter carbide dental drill. The procedures used for U/Th chemical separation and isotopic measurements are detailed in Shao et al. (Reference Shao, Pons-Branchu, Zhu, Wang, Valladas and Fontugne2017). The carbonate samples were weighed and dissolved in 7N HNO3 in Teflon beakers containing a known quantity of a 229Th-233U-236U triple spike. U and Th were preconcentrated by coprecipitation with iron hydroxide and then separated from each other and from other cations by passing the sample solution through a U-TEVA resin column. The U/Th fractions were then dried and diluted in a mixture of 0.1N HNO3 and 0.01N HF for isotopic analysis by a multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS, Neptune) at the School of Geography, Nanjing Normal University, using the methods described in Shao et al. (Reference Shao, Li, Huang, Liao, Arps, Huang and Chou2019). The U isotopic data were acquired by two static sequences. The first sequence measured 233U, 235U, 236U, and 238U in Faraday cups and simultaneously 234U on the secondary electron multiplier (SEM); and the second sequence measured 236U on the SEM and the other isotopes are in Faraday cups. Thorium measurements were carried out immediately after U measurements for the same sample. The 229Th and 230Th were alternately measured on the SEM and 232Th in a Faraday cup. Mass fractionation was corrected by comparing the measured 238U/235U to the natural value of 137.760 for HU-1 and 137.818 for unknown samples (Hiess et al., Reference Hiess, Condon, McLean and Noble2012). Hydride interferences, machine abundance sensitivity, and amplifier gains were evaluated every day prior to measurements. 230Th/U ages were calculated using half-lives of 75,584 yr and 245,620 yr for 230Th and 234U, respectively (Cheng et al., Reference Cheng, Edwards, Shen, Polyak, Asmerom, Woodhead and Hellstrom2013). The 230Th/U age uncertainty was estimated by Monte-Carlo simulations (n = 106). All of the speleothem ages have 2-sigma analytical errors of roughly 0.2–0.4% (Table 1).

Table 1. 230Th dating results for sample HL161 stalagmite from Hulu Cave, China.

*Errors are 2-sigma analytical errors. Decay constant values are λ230 = 9.1705 × 10−6 yr−1, λ234 = 2.82206 × 10−6 yr−1, λ238 = 1.55125 × 10−10 yr−1. Corrected 230Th ages (ka before AD 1950) assume an initial 230Th/232Th atomic ratio of (4.4 ± 2.2) × 10−6.

Only the section from the top to 187-mm depth was used for stable isotope analyses, because the lower part was less densely crystallized. Powdered subsamples (each ~50 µg) were shaved at a resolution of 0.1 mm from along the growth axis with a knife. Every second sample (n = 969) was measured using a Finnigan-MAT 253 mass spectrometer coupled with a Kiel Carbonate Device at the School of Geography, Nanjing Normal University, China. All results are reported in parts per mil (‰) relative to the Vienna Pee Dee Belemnite. Repeated analyses of an international standard (NBS19) indicated long-term reproducibility, with precisions better than 0.06‰ for δ18O and 0.05‰ for δ13C at the 1-sigma level.

A total of 170 powdered samples were drilled with carbide dental burrs along the growth axis and used for trace element analyses at Chongqing Key Laboratory of Karst Environment, Southwest University, China. Each sample weighed 300 ± 50 µg and was dissolved in a solution of 3% HNO3 and 1% HF. Mg and Ca were analyzed using an inductively coupled plasma optical emission spectrometer (ICP-OES, Perkin-Elmer), and Sr and Ba were measured using a single-collector inductively coupled plasma mass spectrometer (SC-ICP-MS, Element XR). International standard SLRS-5 was used to determine the accuracy and precision of the analyses. The precisions are better than 3% for Ca, 3% for Mg, 5% for Ba, and 10% for Sr.

Trace element analyses of bedrock and dripwater samples were made using an inductively coupled plasma atomic emission spectrometry (ICP-AES, Perkin-Elmer Optima) in the Department of Civil and Structural Engineering, Hongkong Polytechnic University, China. Nine bedrock samples and three dripwater samples were collected from inside the cave. Powdered subsamples (~10 mg) were drilled in the unweathered part of each bedrock sample using carbide dental burrs. The subsamples were then digested with concentrated HNO3 and HClO4 (in a ratio of 4:1) and taken to complete dryness on a hot plate. The water samples were filtered through a Millipore membrane (0.45 µm) to remove fine particles, transferred to a polyethylene bottle, and then acidified with 0.2% HNO3. Details of the chemical extraction procedures and the ICP-AES analysis are given in Li et al. (Reference Li, Coles, Ramsey and Thornton1995). Reagent blanks, standard reference materials (NIST 1646a), and sample replicates were used to assess the accuracy and precision of the analyses. The precisions for the measured elements (Mg, Sr, and Ca) are generally <10%.

RESULTS

Chronology

U-Th isotopic compositions and 230Th ages for stalagmite HL161 are presented in Table 1. Measured 238U concentrations range from 155–465 ppb and 232Th concentrations range from 29–781 ppt. The corrections for initial detrital 230Th are negligible, as indicated by the high values of the 230Th/232Th activity ratios (>1000). The age model for stalagmite HL161 was constructed by linear interpolation of the 230Th/U dates (Fig. 1) and the resulting chronology spans the interval of 51.7–42.6 ka, which includes DO 14–11. The average temporal resolution is ~10 yr for stable isotopes (δ18O and δ13C) and ~50 yr for trace elements (Mg/Ca, Sr/Ca, and Ba/Ca; Fig. 2). Thus, the records are well-suited for determining the paleoenvironmental evolution of the monsoonal region of China on millennial to centennial scales.

Figure 2. Multi-proxy records for stalagmite HL161. (A) δ18O records for stalagmite HL161 (green curve; this study), MSL (blue curve; Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001), MSD (grey curve; Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001) and a new composite record (yellow curve; Cheng et al., Reference Cheng, Edwards, Sinha, Spötl, Yi, Chen and Kelly2016). 2-sigma error bars for the records are also shown with corresponding colors. The previously published Hulu records are plotted for comparison with the record of stalagmite HL161. (B–E) are other multi-proxy records from stalagmite HL161. (B) δ13C. (C) Mg/Ca × 10−3. (D) Sr/Ca × 10−3. (E) Ba/Ca × 10−3. The bold lines in (C–E) are 13-point running averages of the raw data and show the long-term trends. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Time series of the proxy records

The stable isotope records vary substantially during the studied period. δ18O varies from -8.3 to -5.9‰ (Fig. 2A) and δ13C from -5.8 to -3.2‰ (Fig. 2B). Three positive δ18O excursions, centered at 49.7, 47.5, and 43.4 ka, enable the entire profile to be divided into four negative phases. These four negative δ18O intervals correspond to GI 14-11 (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001). Rapid δ18O shifts (within <100 yr), with amplitudes >1‰, mark the onset of GI 13, 12, and 11. The changes in δ13C exhibit an antiphased relationship with δ18O on the millennial scale, with 13C enrichment corresponding to 18O depletion. However, this relationship is muted during GI 14, during which a negative period is evident in the δ13C record, with no corresponding shift in the δ18O record.

The “Hendy test” and “Replication test” are used to assess whether calcite precipitation is in equilibrium (Hendy, Reference Hendy1971; Dorale and Liu, Reference Dorale and Liu2009). The negative correlation between δ18O and δ13C in stalagmite HL161 (R = −0.37, P < 0.01) indicates insignificant kinetic isotopic effects (Hendy, Reference Hendy1971). Moreover, the record strongly resembles the pattern of millennial-scale fluctuations previously recorded from Hulu Cave, within the dating uncertainties (Fig. 2A; Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001; Cheng et al., Reference Cheng, Edwards, Sinha, Spötl, Yi, Chen and Kelly2016). Thus, the record passes the “Replication test” (Dorale and Liu, Reference Dorale and Liu2009). Under isotopic equilibrium conditions, the δ18O records of the stalagmites of Hulu Cave reliably reflect the δ18O signal of meteoric precipitation. Modern observations of dripwater and rainwater at the cave site support this interpretation (Wang et al., Reference Wang, Wang, Shao, Liang, Zhang and Kong2018). Following the reasoning of previous studies (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001, Reference Wang, Cheng, Edwards, Kong, Shao, Chen and Wu2008; Liu et al., Reference Liu, Wen, Brady, Otto-Bliesner, Yu, Lu and Cheng2014; Duan et al., Reference Duan, Wu, Wang, Edwards, Cheng, Kong and Zhang2015; Cheng et al., Reference Cheng, Edwards, Sinha, Spötl, Yi, Chen and Kelly2016; Wang et al., Reference Wang, Wang, Shao, Liang, Zhang and Kong2018), we conclude that the stalagmite δ18O record is a proxy of the East Asian summer monsoon (EASM), with intervals of δ18O depletion representing stronger monsoon. However, stalagmite δ18O records are likely to integrate the effects of complex hydroclimatic processes, including moisture sources, degree of Rayleigh fractionation, and seasonal duration and intensity of monsoon rainfall.

The concentrations of Mg, Sr, and Ba of stalagmite HL161 are expressed as ratios to Ca. The Mg/Ca, Sr/Ca and Ba/Ca ratios have the following ranges: 3.4 × 10−3 to 8 × 10−3, 0.1 × 10−3 to 0.4 × 10−3, and 0.2 × 10−3 to 0.3 × 10−3, respectively (Fig. 2C–E). The Sr/Ca and Ba/Ca ratios, in particular, are strongly positively correlated (R = 0.65, n = 170, P < 0.01). The trace element ratios exhibit a coherent pattern of long-term fluctuations (bold lines in Fig. 2C–E): relative stability from the beginning of the record until ~48 ka and then a slight increase followed by an overall decreasing trend until ~42.6 ka. Notably, the trace element ratios exhibit a long-term pattern of variation similar to that of δ13C, with higher values during GI periods, and vice versa.

DISCUSSION

Climatic teleconnection between the East Asian monsoon and North Atlantic regions

There are striking similarities between the δ18Ocalcite record from stalagmite HL161 and the Greenland NGRIP δ18Oice record (Fig. 3A and B; Andersen et al., Reference Andersen, Azuma, Barnola, Bigler, Biscaye, Caillon and Chappellaz2004). First, a stronger EASM indicated by depleted δ18Ocalcite values corresponds closely to Greenland interstadials and vice versa. This millennial-scale coupling of northern, high-latitude climates and low-latitude monsoons has been suggested by several other studies (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001; Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003; Carolin et al., Reference Carolin, Cobb, Adkins, Clark, Conroy, Lejau and Malang2013; Deplazes et al., Reference Deplazes, Lückge, Peterson, Timmermann, Hamann, Hughen and Röhl2013; Chen et al., Reference Chen, Wang, Hai, Edwards, Wang, Kong and Liu2016; Kathayat et al., Reference Kathayat, Cheng, Sinha, Spötl, Edwards, Zhang and Li2016). Second, abrupt monsoonal intensification closely tracks Greenland warming at the start of GI 13, 12, and 11. For example, the onset of GI 12 in the δ18Ocalcite record occurs over ~80 yr, which is close to the 60-yr transition of GI 12 in the δ18Oice record, according to Rampfit analysis of both records (Mudelsee, Reference Mudelsee2000). Third, the duration of GI 12 in the δ18Ocalcite record, which lasts ~2600 yr, is approximately equivalent to the duration of GI 12 estimated by Landais et al. (Reference Landais, Caillon, Goujon, Grachev, Barnola, Chappellaz and Jouzel2004) and Wolff et al. (Reference Wolff, Chappellaz, Blunier, Rasmussen and Svensson2010). In addition, one centennial-scale warming event which stands out from decadal oscillations in the Greenland δ18Oice record (referred to as a “rebound event” by Capron et al. [2010]) is also evident in low-latitude hydroclimatic records (dashed line in Fig. 3; see Deplazes et al. [2013]). These comparisons confirm that low-latitude monsoonal strengthening (weakening) corresponds to increased (decreased) Greenland temperature. Potentially, therefore, a common forcing mechanism exists for the DO cycles evident in both the North Atlantic and the Asian monsoonal regions. Potential mechanisms have been suggested to be associated with the Atlantic Meridional Ocean Circulation (Alley et al., Reference Alley2007) and/or astronomical forcing (Braun et al., Reference Braun, Christl, Rahmstorf, Ganopolski, Mangini, Kubatzki and Roth2005).

Figure 3. Comparison of δ18O records in the monsoonal regions with Greenland paleotemperature. (A) Greenland NGRIP δ18Oice record (Andersen et al., Reference Andersen, Azuma, Barnola, Bigler, Biscaye, Caillon and Chappellaz2004) plotted on the GICC05 timescale (Svensson et al., Reference Svensson, Andersen, Bigler, Clausen, Dahl-Jensen, Davies and Johnsen2008). (B) δ18Ocalcite record from stalagmite HL161 (this study). (C) δ18Ocalcite record from stalagmite M1-2 (Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003). (D) Reconstructed water vapor source. △18O is calculated from the HL161 and M1-2 records and expressed as z-scores. Inferred periods of nearby water vapor sources (z-score △18O > 0.53) are shaded with standard deviations denoted by red dashed lines. The grey dashed line indicates a “rebound event” detected in the NGRIP δ18O ice record (Capron et al., Reference Capron, Landais, Chappellaz and Schilt2010). The blue dashed line denotes an overall cooling trend during GI 12, and the green and orange dashed lines indicate stable monsoon climates. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Previous studies have shown that the GICC05 time scale for Greenland ice cores is systematically too young and that the age offsets increase further back in time (Wang et al., Reference Wang, Cheng, Edwards, An, Wu, Shen and Dorale2001; Fleitmann et al., Reference Fleitmann, Cheng, Badertscher, Edwards, Mudelsee, Göktürk and Fankhauser2009; Buizert et al., Reference Buizert, Cuffey, Severinghaus, Baggenstos, Fudge, Steig and Markel2014). The age offsets can be explained by the presence of uncertain annual layers (Andersen et al., Reference Andersen, Svensson, Johnsen, Rasmussen, Bigler, Rothlisberger and Ruth2005; Svensson et al., Reference Svensson, Andersen, Bigler, Clausen, Dahl-Jensen, Davies and Johnsen2008). Comparison with independent, well-dated time makers is an alternative means of evaluating the controversial GICC05 time scale. The onsets of GI 12 and 11 in the record from Hulu Cave are both in excellent agreement with equivalent GI events in the NGRIP record with nominal differences well within dating uncertainties. Our results are in accord with seven previous studies within dating errors (Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003; Lachniet et al., Reference Lachniet, Johnson, Asmerom, Burns, Polyak, Patterson and Burt2009; Wagner et al., Reference Wagner, Cole, Beck, Patchett, Henderson and Barnett2010; Carolin et al., Reference Carolin, Cobb, Adkins, Clark, Conroy, Lejau and Malang2013; Moseley et al., Reference Moseley, Spötl, Svensson, Cheng, Brandstatter and Edwards2014; Chen et al., Reference Chen, Wang, Hai, Edwards, Wang, Kong and Liu2016; Dong et al., Reference Dong, Shen, Kong, Wang and Duan2018) and thus support the accuracy of the GICC05 time scale within the range of DO 12 and 11.

However, there are noticeable differences in the character of these high- and low-latitude climatic records. For example, during GI 12, Greenland temperature exhibits a decreasing trend (Fig. 3A), while the monsoon remains strong (Fig. 3B and C). In addition, the NGRIP δ18Oice record is characterized by decadal oscillations (10–50 yr; Boers et al., Reference Boers2018), which are different from the pattern of centennial monsoonal variability (~250 yr) superimposed on GI events (Fig. 4). This periodicity of monsoon instability approximates that of the de Vries solar cycle (Wagner et al., Reference Wagner, Beer, Masarik, Muscheler, Kubik, Mende and Laj2001), which may suggest a sensitive monsoonal response to solar activity under warm climatic conditions (Ji et al., Reference Ji, Shen, Balsam, Chen, Liu and Liu2005; Wang et al., Reference Wang, Cheng, Edwards, He, Kong, An and Wu2005).

Figure 4. Results of wavelet analysis of the δ18O record from stalagmite HL161. The horizontal dashed line indicates the ~250-yr band prominent during interstadials. Spectral power is shown by colors ranging from deep blue (weak) to deep red (strong). The 95% significance level against red noise is shown by the thick solid line. The spectra were estimated using the method of Grinsted et al. (Reference Grinsted, Moore and Jevrejeva2004; http://noc.ac.uk/using-science/crosswavelet-wavelet-coherence). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

The monsoon system is driven by two primary mechanisms: differential sensible heating between land and sea, and variations in latent heat exported from the southern subtropical Indian Ocean via cross-equatorial airflows (An et al., Reference An, Wu, Li, Sun, Liu, Zhou and Cai2015). A temporally variable relationship under different climatic conditions between the EASM and the Northern Hemisphere “pull” factors and Southern Hemisphere “push” factors was determined by Rohling et al. (Reference Rohling, Liu, Roberts, Stanford, Rasmussen, Langen and Siddall2009). According to their study, a southern “push” effect has a stronger control on millennial-scale monsoon variability during glacial states, and a northern “pull” effect has a stronger control during deglacial and interglacial states (Rohling et al., Reference Rohling, Liu, Roberts, Stanford, Rasmussen, Langen and Siddall2009). The striking similarity between the NGRIP and HL161 records in terms of their general trends and the abruptness of DO onsets indicates a rapid transmission via an atmospheric “pull” effect. Thus, we suggest that the “push” and “pull” processes affecting the variability of the cave δ18O signals are different between southwestern and eastern/northern China, leading to different DO patterns (Cai et al., Reference Cai, An, Cheng, Edwards, Kelly, Liu and Wang2006; Han et al., Reference Han, Li, Cheng, Edwards, Shen, Li and Huang2016; Duan et al., Reference Duan, Cheng, Tan and Edwards2016).

Dynamic changes of the westerly jet over China may be an important component of the climatic teleconnection between northern high latitudes and the Asian monsoonal region (Nagashima et al., Reference Nagashima, Tada, Matsui, Irino, Tani and Toyoda2007, Reference Nagashima, Tada, Tani, Sun, Isozaki, Toyoda and Hasegawa2011; Molnar et al., Reference Molnar, Boos and Battisti2010; Orland et al., Reference Orland, Edwards, Cheng, Kozdon, Cross and Valley2015). Notably, changes in the timing of the EASM are consistent with the shift of the westerly jet from a location to the south of the Tibetan Plateau to a location to the north (Molnar et al., Reference Molnar, Boos and Battisti2010; Orland et al., Reference Orland, Edwards, Cheng, Kozdon, Cross and Valley2015). In addition, the jet responds sensitively to northern high-latitude climate, including storm-track movements, sea-ice extent, and sea surface temperature (Laîné et al., Reference Laîné, Kageyama, Salas-Mélia, Voldoire, Rivière, Ramstein and Planton2009). During cold GS phases, the westerly jet may be weaker and its northward movement is delayed, essentially maintaining East Asia in a prolonged dry pre-monsoonal state (Orland et al., Reference Orland, Edwards, Cheng, Kozdon, Cross and Valley2015). The southern location of the westerly jet prevents isotopically light water vapor from penetrating into East Asia, resulting in heavier mean δ18O of rainfall. During the transition from GS to GI climate states, following the abrupt warming in northern high latitudes, a rapid northward movement of the jet may allow rapid and large-scale, low-level monsoonal flow into the interior of East Asia, with a corresponding rapid decrease in δ18Ocalcite values. This mechanism is supported by changes in the water vapor supply of EASM rainfall, as discussed below.

The high degree of similarity of the HL161 δ18O record to other records from Chinese and Indian caves confirms that the δ18O changes are of regional extent (Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003; Chen et al., Reference Chen, Wang, Hai, Edwards, Wang, Kong and Liu2016; Kathayat et al., Reference Kathayat, Cheng, Sinha, Spötl, Edwards, Zhang and Li2016; Dong et al., Reference Dong, Shen, Kong, Wang and Duan2018). One possible contributor to δ18O variability is the water vapor supply from source regions (Dayem et al., Reference Dayem, Molnar, Battisti and Roe2010; Baker et al., Reference Baker, Sodemann, Baldini, Breitenbach, Johnson, van Hunen and Zhang2015; Wu et al., Reference Wu, Zhang, Li, Li and Huang2015). In this study, we removed the δ18O signals of remote water vapor sources from the HL161 record by comparing it with a high-resolution calcite δ18O record (M1-2) from Moomi Cave in Socotra Island. Socotra Island is located upstream of the EASM - the Indian Ocean (Baker et al., Reference Baker, Sodemann, Baldini, Breitenbach, Johnson, van Hunen and Zhang2015) and is solely influenced by the Indian summer monsoon (Fig. 3C; Burns et al., Reference Burns, Fleitmann, Matter, Kramers and Al-Subbary2003). The two caves are linked in that they are both under the control of large-scale monsoonal circulation (Wang et al., Reference Wang, Wang, Cheng, Fasullo, Guo, Kiefer and Liu2014). According to the reconstructed global gridded precipitation prediction map of δ18O (Terzer et al., Reference Terzer, Wassenaar, Araguásaraguás and Aggarwal2013), the δ18O values decrease progressively along the pathway of water vapor from Socotra Island (-2.9 to 0.0‰) to South China (-8.9 to -6.0‰). The difference between the average δ18O values at the two cave sites (6.8%) is well within the range of the δ18O difference (8.9 to 3.1‰) between the water vapor source (the Indian Ocean) and the site of rainout (Nanjing), which could have been caused by Rayleigh fractionation during transport (Dansgaard, Reference Dansgaard1964).

Accordingly, we subtracted the M1-2 δ18O record from the HL161 δ18O record, after tuning M1-2 to the HL161 chronology with the aid of the very distinctive GI 12. We then performed a z-score transform (standard deviation = 0.53) and the results are expressed as z-scored △δ18O (Fig. 3D). The △δ18O values represent changes in water vapor supply from relatively close sources, including the West Pacific Ocean and/or inland China. The record enables us to elucidate the control of the westerly jet over the water vapor source of the EASM rainfall (Fig. 3D). Three periods with positive values, centered at 49.7, 47.6, and 43.7 ka and corresponding respectively to GS 14, 13, and 12, are prominent within the record. The positive △δ18O values represent increased water vapor supply from nearby sources due to the delayed shift in the position of the westerly jet. We then applied Empirical Mode Decomposition (Huang et al., Reference Huang, Shen and Long1998; Huang and Wu, Reference Huang and Wu2008) to determine the water vapor contribution from the Indian Ocean and the West Pacific Ocean/inland China during GS 14 to 12 (excluding the incomplete events GI 14 and 11). The estimated average water vapor proportions derived from distant and sources near Hulu Cave are 56.7 and 43.3%, respectively. However, during GS states, the nearby water vapor source comprises a higher proportion, with values of 47.6, 52.9 and 69.0% for GS 14, 13, and 12, respectively.

Response of cave hydrological cycles to monsoonal changes

The high degree of consistency of the δ13C records from Hulu Cave (Fig. 5A and B) suggests that they reflect a common environmental process (Dorale and Liu, Reference Dorale and Liu2009). The antiphased relationship between δ13C and δ18O records of stalagmite HL161 on the millennial scale suggests that δ13C record is sensitive to monsoonal changes (Fig. 2A and B). A possible explanation for this correlation is a “damping model” which reflects local hydrology, with changes linked to rainfall (Kong et al., Reference Kong, Wang, Wu, Cheng, Edwards and Wang2005). A “damping model” requires a specific combination of the conditions of speleothem formation and climate factors: a thin soil cover and elevated precipitation (Baker et al., Reference Baker, Ito, Smart and Mcewan1997). Modern field observations show that the site of Hulu Cave has a thin soil cover, <30 cm deep, while rainfall is the main contributor to the recharge of the aquifer above the cave. The aquifer is also highly permeable, which allows rainwater to infiltrate rapidly within a few days after heavy rainfalls. This combination provides the necessary conditions for the “damping model.” Under heavy rainfall, equilibrium between soil water and soil CO2 cannot be attained due to the short residence time of the water, which results in a low degree of depletion of 13C in the speleothem (Baker et al., Reference Baker, Ito, Smart and Mcewan1997; Kong et al., Reference Kong, Wang, Wu, Cheng, Edwards and Wang2005). The inverse process operates under drier and colder climatic conditions. Notably, the similarity of speleothem δ13C and δ18O records is also demonstrated during the penultimate glacial period (Wang et al., Reference Wang, Wang, Shao, Liang, Zhang and Kong2018).

Figure 5. (A and B) Results of a replication test of δ13C records from Hulu Cave. (A) δ13C record for stalagmite HL161. (B) δ13C records from stalagmite MSL (brown) and stalagmite MSD (gold; Kong et al., Reference Kong, Wang, Wu, Cheng, Edwards and Wang2005). (C and D) Results of an analysis of trace element from Hulu Cave. (C) Relationship between ln(Mg/Ca) versus ln(Sr/Ca) for stalagmite HL161 with the slope of the best-fit regression line indicated. (D) Comparison of values of Sr/Ca and Mg/Ca for bedrock (black dots), cave water (red dots), and stalagmite HL161 (purple crosses). The cave water and stalagmite samples are depleted in Mg and Sr relative to Ca in comparison to bedrock. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

The similarity between the trace element ratios (Fig. 2C– E) on the millennial scale also suggests that they reflect a common set of local processes. Trends in stalagmite composition are ultimately derived from variations in the seepage water which are primarily controlled by local soil, regolith, and bedrock (Fairchild and Treble, Reference Fairchild and Treble2009). Processes such as PCP, differential dissolution of calcite and dolomite, selective leaching, and acids in water are potential controls of the trace element composition of speleothems (Fairchild et al., Reference Fairchild, Borsato, Tooth, Frisia, Hawkesworth, Huang and McDermott2000; Hellstrom and McCulloch, Reference Hellstrom and McCulloch2000; Fairchild and Baker, Reference Fairchild and Baker2012). Among them, PCP is one of the most-often cited explanations for the covariation of trace elements (Johnson et al., Reference Johnson, Hu, Belshaw and Henderson2006; Fairchild and Treble, Reference Fairchild and Treble2009; Sinclair et al., Reference Sinclair, Banner, Taylor, Partin, Jenson, Mylroie and Goddard2012; Stoll et al., Reference Stoll, Müller and Prieto2012). Sinclair et al. (Reference Sinclair, Banner, Taylor, Partin, Jenson, Mylroie and Goddard2012) presented a mathematical model of PCP and showed that the positive linear correlation between ln(Mg/Ca) versus ln(Sr/Ca) was given by $\displaystyle{{{\rm K}{\rm d}_{{\rm Sr}}-1} \over {{\rm K}{\rm d}_{{\rm Mg}}-1}}$. This is a universal property of PCP and can be applied even when there is no prior characterization of the dripwater or host rock (Sinclair et al., Reference Sinclair, Banner, Taylor, Partin, Jenson, Mylroie and Goddard2012). Here, we used this model to calculate the slope of the theoretical PCP gradient in Hulu Cave. At a mean annual cave temperature of ~15°C, the partition coefficients adopted here are KdMg = 0.019 (Huang and Fairchild, Reference Huang and Fairchild2001) and KdSr = 0.1 (Treble et al., Reference Treble, Fairchild, Griffiths, Baker, Meredith, Wood and McGuire2015). Assuming that the partition coefficients are constant and that PCP is the only process occurring in a closed water-rock system, the theoretical PCP gradient in Hulu Cave will be 0.92. However, the calculated slope for stalagmite HL161 is much smaller (0.21, Fig. 5C), thus indicating minor PCP effects during deposition.

At Hulu Cave, the congruent Ordovician limestone (Wang et al., Reference Wang, Wang, Shao, Liang, Zhang and Kong2018) indicates the absence of differential dissolution (Fairchild et al., Reference Fairchild, Borsato, Tooth, Frisia, Hawkesworth, Huang and McDermott2000). Selective leaching is also unlikely to occur. Compared with the bedrock mixing line, the concentrations of Sr and Mg relative to Ca in the dripwater and stalagmite are depleted (Fig. 5D; Fairchild and Treble, Reference Fairchild and Treble2009). Therefore, we exclude these factors as possible mechanisms for the covariation.

Our favored explanation relates to the acid concentrations in the soil water. The evidence for this is the case of Ba, which is a relatively immobile element due to the high cation exchange selectivity of Ba2+ (Mcbride, Reference McBride1994). Biological activity, however, can increase the weathering and leaching of Ba2+ by producing more concentrated carbonic acid, hydrogen ions, and organic exudates (Taylor et al., Reference Taylor, Leake, Quirk, Hardy, Banwart and Beerling2009). Under warm and wet conditions, these reactive species are supplied to the infiltrating water by elevated root respiration, litter decomposition, and other biological activities (Brook et al., Reference Brook, Folkoff and Box1983; Taylor et al., Reference Taylor, Leake, Quirk, Hardy, Banwart and Beerling2009; Plestenjak et al., Reference Plestenjak, Eler, Vodnik, Ferlan, Čater, Kanduč and Simončič2012), thus leading to increased concentrations of Ba2+ relative to Ca2+ (Hellstrom and McCulloch, Reference Hellstrom and McCulloch2000). The positive correlation between Sr/Ca and Ba/Ca ratios (R = 0.65, n = 170, P < 0.01) and the covariation of their long-term trends suggest that similar processes apply to Sr and Mg. The mobilized cations (Ba, Sr, and Mg) are thereby introduced into the cave environment from surficial soils and host rock by infiltrating water.

The variations in proxy δ13C and metal/Ca ratios strongly resemble the HL161 δ18O record. On the millennial scale, more negative δ18O values are consistent with the elevated δ13C values and Mg/Ca, Sr/Ca, and Ba/Ca ratios. This correlation even extends to the centennial scale, such as the episodes of monsoon strengthening at 51.2, 46.5, 45.8, and 44.3 ka. The general agreement between trace element ratios and stable isotope values indicates that both hydrological and biological activity at the local level responded sensitively to changes in monsoon.

CONCLUSIONS

We have produced a high-resolution and 230Th/U-dated multi-proxy record from a new stalagmite (HL161) from Hulu Cave. The records span the interval of 51.7–42.6 ka that encompasses DO 14–11. The records of the geochemical indicators of δ18O, δ13C, Mg/Ca, Sr/Ca, and Ba/Ca provide new perspectives on the EASM and cave hydrological cycles on millennial to centennial timescales.

The ~10-yr-resolution δ18O record is strikingly similar to that of the Greenland ice-core record in terms of the general pattern, suggesting a common forcing factor of DO cycles in both high and low latitudes of the Northern Hemisphere. The rapid transitions at the onset of GI events in our stalagmite record resemble the abrupt DO warming events in the NGRIP record and they are synchronous within dating errors. However, centennial-scale variations (~250 yr) during Chinese interstadials differ from the predominantly decadal-scale (10–50 yr) oscillations in the climate over Greenland. The surprisingly rapid response of the subtropical monsoon to Greenland temperature variations indicates a rapid atmospheric transmission mechanism. After removing the variance associated water vapor from the Indian Ocean from our δ18O record, the proportion of the nearby water vapor sources is higher during GS phases, implying that the westerly jet plays a key role in controlling the water vapor transport to eastern China. Therefore, we conclude that the calcite δ18O record of Hulu Cave can largely be interpreted as reflecting changes in monsoon intensity, which incorporates changes in water vapor sources.

The δ13C record of stalagmites from Hulu Cave reflects the effects of a “damping model,” which are linked with local rainfall. The Mg/Ca, Sr/Ca, and Ba/Ca ratios of stalagmite HL161 reflect the effects of biological activity within the soil above the cave during water infiltration. On the millennial and even the centennial scale, changes in δ13C and metal/Ca ratios resemble the δ18O record, with higher δ13C and metal/Ca ratios corresponding to lower δ18O values. This observation indicates that changes in local hydrology and ecosystem respond sensitively to variations in monsoon.

ACKNOWLEDGEMENTS

The authors are grateful to Jeffrey Dorale and two anonymous reviewers for their valuable comments, which substantially improved the manuscript. This research was supported by the National Nature Science Fund of China (awards 41571102, 41672164, 41372174, and 41130210), the Priority Academic Program Development of Jiangsu Higher Education Institutions (award 164320H116), the Jiangsu Center for Collaborative Innovation in Geographical Information Resource development and Application, the 111 program of China (approved Number: D19002), and United States National Science Foundation Grant 1702816.

References

REFERENCES

Alley, R.B., 2007. Wally was right: predictive ability of the North Atlantic “Conveyor Belt” Hypothesis for abrupt climate change. Annual Review of Earth and Planetary Sciences 35, 241272.Google Scholar
An, Z., Wu, G., Li, J., Sun, Y., Liu, Y., Zhou, W., Cai, Y., et al. , 2015. Global monsoon dynamics and climate change. Journal of Earth Environment 43, 2977.Google Scholar
Andersen, K.K., Azuma, N., Barnola, J.M., Bigler, M., Biscaye, P., Caillon, N., Chappellaz, E., et al. , 2004. High-resolution record of Northern Hemisphere climate extending into the last interglacial period. Nature 431, 147151.Google Scholar
Andersen, K.K., Svensson, A., Johnsen, S., Rasmussen, S.O., Bigler, M., Rothlisberger, R., Ruth, U., et al. , 2005. The Greenland ice core chronology 2005, 15–42 ka. Part 1: constructing the time scale. Quaternary Science Reviews 25, 32463257.Google Scholar
Baker, A., Ito, E., Smart, P.L., Mcewan, R.F., 1997. Elevated and variable values of δ13C in speleothems in a British cave system. Chemical Geology 136, 263270.Google Scholar
Baker, A.J., Sodemann, H., Baldini, J.U.L., Breitenbach, S.F.M., Johnson, K.R., van Hunen, J., Zhang, P., 2015. Seasonality of westerly moisture transport in the East Asian summer monsoon and its implications for interpreting precipitation δ18O. Journal of Geophysical Research Atmospheres 120, 58505862.Google Scholar
Beck, J.W., Zhou, W., Li, C., Wu, Z., White, L., Xian, F., Kong, X., et al. , 2018. A 550,000-year record of East Asian monsoon rainfall from 10Be in loess. Science 360, 877881.Google Scholar
Boers, N., 2018. Early-warning signals for Dansgaard-Oeschger events in a high-resolution ice core record. Nature Communications 9, 2556.Google Scholar
Borsato, A., Frisia, S., Fairchild, I.J., Somogyi, A., Susini, J., 2007. Trace element distribution in annual stalagmite laminae mapped by micrometer-resolution x-ray fluorescence: implications for incorporation of environmentally significant species. Geochimica et Cosmochimica Acta 71, 14941512.Google Scholar
Braun, H., Christl, M., Rahmstorf, S., Ganopolski, A., Mangini, A., Kubatzki, C., Roth, K., et al. , 2005. Possible solar origin of the 1,470-year glacial climate cycle demonstrated in a coupled model. Nature 438, 208211.Google Scholar
Brook, G.A., Folkoff, M.E., Box, E.O., 1983. A world model of soil carbon dioxide. Earth Surface Processes and Landforms 8, 7988.Google Scholar
Buizert, C., Cuffey, K.M., Severinghaus, J.P., Baggenstos, D., Fudge, T.J., Steig, E.J., Markel, B.R., et al. , 2014. The WAIS-Divide deep ice core WD2014 chronology-part 2: methane synchronization (68–31 ka BP) and the gas age-ice age difference. Climate of the Past 10, 153173.Google Scholar
Burns, S.J., Fleitmann, D., Matter, A., Kramers, J., Al-Subbary, A.A., 2003. Indian Ocean climate and an absolute chronology over Dansgaard/Oeschger events 9 to 13. Science 301, 13651367.Google Scholar
Cai, Y., An, Z., Cheng, H., Edwards, R.L., Kelly, M.J., Liu, W., Wang, X., et al. , 2006. High-resolution absolute-dated Indian monsoon record between 53 and 36 ka from Xiaobailong Cave, southwestern China. Geology 34, 621624.Google Scholar
Cai, Y., Tan, L., Cheng, H., An, Z., Edwards, R.L., Kelly, M.J., Kong, X., et al. , 2010. The variation of summer monsoon precipitation in central China since the last deglaciation. Earth and Planetary Science Letters 291, 2131.Google Scholar
Capron, E., Landais, A., Chappellaz, J., Schilt, A., 2010. Millennial and sub-millennial scale climatic variations recorded in polar ice cores over the last glacial period. Climate of the Past 6, 345365.Google Scholar
Carolin, S.A., Cobb, K.M., Adkins, J.F., Clark, B., Conroy, J.L., Lejau, S., Malang, J., et al. , 2013. Varied response of western pacific hydrology to climate forcings over the last glacial period. Science 340, 15641566.Google Scholar
Chen, C., Li, T., 2018. Geochemical characteristics of cave drip water respond to ENSO based on a 6-year monitoring work in Yangkou cave, southwest China. Journal of Hydrology 561, 896907.Google Scholar
Chen, S., Wang, Y., Hai, C., Edwards, R.L., Wang, X., Kong, X., Liu, D., 2016. Strong coupling of Asian monsoon and Antarctic climates on sub-orbital timescales. Scientific Reports 6, 32995.Google Scholar
Cheng, H., Edwards, R.L., Shen, C., Polyak, V.J., Asmerom, Y., Woodhead, J., Hellstrom, J., et al. , 2013. Improvements in 230Th dating, 230Th and 234U half-life values, and U-Th isotopic measurements by multi-collector inductively coupled plasma mass spectroscopy. Earth and Planetary Science Letters 371–372, 8291.Google Scholar
Cheng, H., Edwards, R.L., Sinha, A., Spötl, C., Yi, L., Chen, S., Kelly, M., et al. , 2016. The Asian monsoon over the past 640,000 years and ice age terminations. Nature 534, 640646.Google Scholar
Dansgaard, W., 1964. Stable isotopes in precipitation. Tellus 16, 436468.Google Scholar
Dansgaard, W., Johnsen, S.J., Clausen, H.B., Dahl-Jensen, D., Gundestrup, N.S., Hammer, C.U., Hvidberg, C.S., et al. , 1993. Evidence for general instability of past climate from a 250-kyr ice-core record. Nature 364, 218220.Google Scholar
Dayem, K.E., Molnar, P., Battisti, D.S., Roe, G. H., 2010. Lessons learned from oxygen isotopes in modern precipitation applied to interpretation of speleothem records of paleoclimate from eastern Asia. Earth and Planetary Science Letters 295, 219230.Google Scholar
Deplazes, G., Lückge, A., Peterson, L.C., Timmermann, A., Hamann, Y., Hughen, K.A., Röhl, U., et al. , 2013. Links between tropical rainfall and North Atlantic climate during the last glacial period. Nature Geoscience 6, 213217.Google Scholar
Dong, J., Shen, C., Kong, X., Wang, Y., Duan, F., 2018. Asian monsoon dynamics at Dansgaard/Oeschger events 14–8 and Heinrich events 5–4 in northern China. Quaternary Geochronology 47, 7280.Google Scholar
Dorale, J.A., Liu, Z., 2009. Limitations of Hendy test criteria in judging the paleoclimatic suitability of speleothems and the need for replication. Journal of Cave and Karst Studies 71, 7380.Google Scholar
Duan, F., Wu, J., Wang, Y., Edwards, R.L., Cheng, H., Kong, X., Zhang, W., 2015. A 3000-yr annually laminated stalagmite record of the last glacial maximum from Hulu Cave, China. Quaternary Research 83, 360369.Google Scholar
Duan, W., Cheng, H., Tan, M., Edwards, R.L., 2016. Onset and duration of transitions into Greenland interstadials 15.2 and 14 in northern China constrained by an annually laminated stalagmite. Scientific Reports 6, 20844Google Scholar
Fairchild, I.J., Baker, A., 2012. Speleothem science: from process to past environments. John Wiley & Sons, Chichester.Google Scholar
Fairchild, I.J., Borsato, A., Tooth, A.F., Frisia, S., Hawkesworth, C.J., Huang, Y., McDermott, F., et al. , 2000. Controls on trace element (Sr–Mg) compositions of carbonate cave waters: implications for speleothem climatic records. Chemical Geology 166, 255269.Google Scholar
Fairchild, I.J., Treble, P.C., 2009. Trace elements in speleothems as recorders of environmental change. Quaternary Science Reviews 28, 449468.Google Scholar
Fleitmann, D., Cheng, H., Badertscher, S., Edwards, R.L., Mudelsee, M., Göktürk, O.M., Fankhauser, A., et al. , 2009. Timing and climatic impact of Greenland interstadials recorded in stalagmites from northern Turkey. Geophysical Research Letters 36, L19707.Google Scholar
Genty, D., Blamart, D., Ouahdi, R., Gilmour, M., Baker, A., Jouzel, J., Van-Exter, S., 2003. Precise dating of Dansgaard-Oeschger climate oscillations in Western Europe from stalagmite data. Nature 421, 833837.Google Scholar
Goldsmith, Y., Broecker, W.S., Xu, H., Polissar, P.J., Demenocal, P.B., Porat, N., Lan, J., et al. , 2017. Northward extent of East Asian monsoon covaries with intensity on orbital and millennial timescales. Proceedings of the National Academy of Sciences of the United States of America 114, 18171821.Google Scholar
Grinsted, A., Moore, J.C., Jevrejeva, S., 2004. Application of the cross wavelet transform and wavelet coherence to geophysical time series. Nonlinear Processes in Geophysics 11, 561566.Google Scholar
Han, L., Li, T., Cheng, H., Edwards, R.L., Shen, C., Li, H., Huang, C., et al. , 2016. Potential influence of temperature changes in the Southern Hemisphere on the evolution of the Asian summer monsoon during the last glacial period. Quaternary International 392, 239250.Google Scholar
Hellstrom, J.C., McCulloch, M.T., 2000. Multi-proxy constraints on the climatic significance of trace element records from a New Zealand speleothem. Earth and Planetary Science Letters 179, 287297.Google Scholar
Hendy, C.H., 1971. The isotopic geochemistry of speleothems - I. The calculation of the effects of different modes of formation on the isotopic composition of speleothems and their applicability as palaeoclimatic indicators. Geochimica et Cosmochimica Acta 3, 801824.Google Scholar
Hiess, J., Condon, D.J., McLean, N., Noble, S.R., 2012. 238U/235U Systematics in terrestrial uranium-bearing minerals. Science 335, 16101614.Google Scholar
Huang, N., Shen, Z., Long, S., 1998. The empirical mode de-composition and the Hilbert spectrum for nonlinear and non-stationary time series analysis. Proceedings of the Royal Society A Mathematica 454, 903995.Google Scholar
Huang, N., Wu, Z., 2008. A review on Hilbert-Huang transform: method and its applications to geophysical studies. Reviews of Geophysics 46, RG2006. http://dx.doi.org/10.1029/2007rg000228Google Scholar
Huang, W., Wang, Y., Cheng, H., Edwards, R.L., Shen, C., Liu, D., Shao, Q., et al. , 2016. Multi-scale Holocene Asian monsoon variability deduced from a twin-stalagmite record in southwestern China. Quaternary Research 86, 3444.Google Scholar
Huang, Y., Fairchild, I.J., 2001. Partitioning of Sr2+ and Mg2+ into calcite under karst-analogue experimental conditions. Geochimica et Cosmochimica Acta 65, 4762.Google Scholar
Ji, J., Shen, J., Balsam, W., Chen, J., Liu, L., Liu, X., 2005. Asian monsoon oscillations in the northeastern Qinghai-Tibet Plateau since the late glacial as interpreted from visible reflectance of Qinghai Lake sediments. Earth and Planetary Science Letters 233, 6170.Google Scholar
Johnsen, S.J., Clausen, H.B., Dansgaard, W., Fuhrer, K., Gundestrup, N., Hammer, C.U., Iversen, P., et al. , 1992. Irregular glacial interstadials recorded in a new Greenland ice core. Nature 359, 311313.Google Scholar
Johnson, K.R., Hu, C., Belshaw, N.S., Henderson, G.M., 2006. Seasonal trace-element and stable-isotope variations in a Chinese speleothem: the potential for high-resolution paleomonsoon reconstruction. Earth and Planetary Science Letters 244, 394407.Google Scholar
Kathayat, G., Cheng, H., Sinha, A., Spötl, C., Edwards, R.L., Zhang, H., Li, X., et al. , 2016. Indian monsoon variability on millennial-orbital timescales. Scientific Reports 6, 24374.Google Scholar
Kelly, M., Edwards, R.L., Cheng, H., Yuan, D., Cai, Y., Zhang, M., Lin, Y., et al. , 2006. High resolution characterization of the Asian monsoon between 146,000 and 99,000 years BP from Dongge Cave, China and global correlation of events surrounding termination II. Palaeogeography Palaeoclimatology Palaeoecology 236, 2038.Google Scholar
Kong, X., Wang, Y., Wu, J., Cheng, H., Edwards, R.L., Wang, X., 2005. Complicated responses of stalagmite δ13C to climate change during the last glaciation from Hulu Cave, Nanjing, China. Science in China Series D-Earth Sciences 48, 21742181.Google Scholar
Lachniet, M.S., Johnson, L., Asmerom, Y., Burns, S.J., Polyak, V., Patterson, W. P., Burt, L., et al. , 2009. Late quaternary moisture export across Central America and to Greenland: evidence for tropical rainfall variability from Costa Rican stalagmites. Quaternary Science Reviews 28, 33483360.Google Scholar
Laîné, A., Kageyama, M., Salas-Mélia, D., Voldoire, A., Rivière, G., Ramstein, G., Planton, S., et al. , 2009. Northern Hemisphere storm tracks during the last glacial maximum in the PMIP2 ocean-atmosphere coupled models: energetic study, seasonal cycle, precipitation. Climate Dynamics 32, 593614.Google Scholar
Landais, A., Caillon, N., Goujon, C., Grachev, A.M., Barnola, J.M., Chappellaz, J., Jouzel, J., et al. , 2004. Quantification of rapid temperature change during DO event 12 and phasing with methane inferred from air isotopic measurements. Earth and Planetary Science Letters 225, 221232.Google Scholar
Li, X., Coles, B.J., Ramsey, M.H., Thornton, I., 1995. Sequential extraction of soils for multielement analysis by ICP-AES. Chemical Geology 124, 109123.Google Scholar
Liu, Z., Wen, X., Brady, E.C., Otto-Bliesner, B., Yu, G., Lu, H., Cheng, H., et al. , 2014. Chinese cave records and the East Asia Summer Monsoon. Quaternary Science Reviews 83, 115128.Google Scholar
Maher, B.A., 2008. Holocene variability of the East Asian summer monsoon from Chinese cave records: a re-assessment. Holocene 18, 861866.Google Scholar
McBride, M.B., 1994. Environmental Chemistry of Soils. Oxford University Press, New York.Google Scholar
Molnar, P., Boos, W.R., Battisti, D.S., 2010. Orographic controls on climate and paleoclimate of Asia: Thermal and mechanical roles for the Tibetan Plateau. Annual Review of Earth and Planetary Sciences 38, 77102.Google Scholar
Moseley, G.E., Spötl, C., Svensson, A., Cheng, H., Brandstatter, S., Edwards, R.L., 2014. Multi-speleothem record reveals tightly coupled climate between central Europe and Greenland during Marine Isotope Stage 3. Geology 42, 10431046.Google Scholar
Mudelsee, M., 2000. Ramp function regression: a tool for quantifying climate transitions. Computer and Geosciences 26, 293307.Google Scholar
Nagashima, K., Tada, R., Matsui, H., Irino, T., Tani, A., Toyoda, S., 2007. Orbital- and millenial-scale variations in Asian dust transport path to the Japan Sea. Palaeogeography, Palaeoclimatology, Palaeoecology 247, 144161.Google Scholar
Nagashima, K., Tada, R., Tani, A., Sun, Y., Isozaki, Y., Toyoda, S., Hasegawa, H., 2011. Millennial-scale oscillations of the westerly jet path during the last glacial period. Journal of Asian Earth Sciences 40, 12141220.Google Scholar
Orland, I.J., Edwards, R.L., Cheng, H., Kozdon, R., Cross, M., Valley, J.W., 2015. Direct measurements of deglacial monsoon strength in a Chinese stalagmite. Geology 43, 555558.Google Scholar
Pausata, F.S.R., Battisti, D.S., Nisancioglu, K.H., Bitz, C.M., 2011. Chinese stalagmite δ18O controlled by changes in the Indian monsoon during a simulated Heinrich event. Nature Geoscience 4, 474480.Google Scholar
Plestenjak, G., Eler, K., Vodnik, D., Ferlan, M., Čater, M., Kanduč, T., Simončič, P., et al. , 2012. Sources of soil CO2 in calcareous grassland with woody plant encroachment. Journal of Soils and Sediments 12, 13271338.Google Scholar
Rohling, E.J., Liu, Q., Roberts, A.P., Stanford, J.D., Rasmussen, S.O., Langen, P. L., Siddall, M., 2009. Controls on the East Asian monsoon during the last glacial cycle, based on comparison between Hulu Cave and polar ice-core records. Quaternary Science Reviews 28, 32913302.Google Scholar
Rousseau, D.D., Svensson, A., Bigler, M., Sima, A., Steffensen, J.P., Boers, N., 2017. Eurasian contribution to the last glacial dust cycles: how are loess sequences built? Climate of the Past 13, 11811197.Google Scholar
Shao, Q., Li, C., Huang, M., Liao, Z., Arps, J., Huang, C., Chou, Y., et al. , 2019. Interactive programs of MC-ICPMS data processing for 230Th/U geochronology. Quaternary Geochronology 51, 4352.Google Scholar
Shao, Q., Pons-Branchu, E., Zhu, Q., Wang, W., Valladas, H., Fontugne, M., 2017. High precision U/Th dating of the rock paintings at Mt. Huashan, Guangxi, southern China. Quaternary Research 88, 113.Google Scholar
Sinclair, D.J., Banner, J. L., Taylor, F.W., Partin, J., Jenson, J., Mylroie, J., Goddard, E., et al. , 2012. Magnesium and strontium systematics in tropical speleothems from the western pacific. Chemical Geology 294–295, 117.Google Scholar
Stoll, H.M., Müller, W., Prieto, M., 2012. I-STAL, a model for interpretation of Mg/ca, Sr/ca and Ba/ca variations in speleothems and its forward and inverse application on seasonal to millennial scales. Geochemistry, Geophysics, Geosystems 13, Q09004, doi:10.1029/2012GC004183.Google Scholar
Svensson, A., Andersen, K.K., Bigler, M., Clausen, H.B., Dahl-Jensen, D., Davies, S.M., Johnsen, S.J., et al. , 2008. A 60,000 year Greenland stratigraphic ice core chronology. Climate of the Past Discussions 3, 4757.Google Scholar
Tan, L., Cai, Y., An, Z., Cheng, H., Shen, C., Gao, Y., Edwards, R.L., 2017. Decreasing monsoon precipitation in southwest China during the last 240 years associated with the warming of tropical ocean. Climate Dynamics 48, 17691778.Google Scholar
Tan, L., Cai, Y., Cheng, H., Edwards, R.L., Gao, Y., Xu, H., Zhang, H., et al. , 2018. Centennial-to decadal-scale monsoon precipitation variations in the upper Hanjiang River region, China over the past 6650 years. Earth and Planetary Science Letters 482, 580590.Google Scholar
Taylor, L.L., Leake, J.R., Quirk, J., Hardy, K., Banwart, S.A., Beerling, D.J., 2009. Biological weathering and the long-term carbon cycle: integrating mycorrhizal evolution and function into the current paradigm. Geobiology 7, 171191.Google Scholar
Terzer, S., Wassenaar, L.L., Araguásaraguás, L.J., Aggarwal, P.K., 2013. Global isoscapes for δ18O and δ2H in precipitation: improved prediction using regionalized climatic regression models. Hydrology and Earth System Sciences 17, 47134728.Google Scholar
Treble, P., Fairchild, I.J., Griffiths, A., Baker, A., Meredith, K.T., Wood, A., McGuire, E., 2015. Impacts of cave air ventilation and in-cave prior calcite precipitation on Golgotha Cave dripwater chemistry, southwest Australia. Quaternary Science Reviews 127, 6172.Google Scholar
Treble, P., Shelley, J.M.G., Chappell, J., 2003. Comparison of high resolution sub-annual records of trace elements in a modern (1911–1992) speleothem with instrumental climate data from southwest Australia. Earth and Planetary Science Letters 216, 141153.Google Scholar
Voelker, A.H.L., 2002. Global distribution of centennial-scale records for Marine Isotope Stage (MIS) 3: a database. Quaternary Science Reviews 21, 11851212.Google Scholar
Wagner, G., Beer, J., Masarik, J., Muscheler, R., Kubik, P.W., Mende, W., Laj, C., et al. , 2001. Presence of the solar de Vries cycle (~205 years) during the last ice age. Geophysical Research Letters 28, 303306.Google Scholar
Wagner, J.D.M., Cole, J.E., Beck, J.W., Patchett, P.J., Henderson, G.M., Barnett, H.R., 2010. Moisture variability in the southwestern United States linked to abrupt glacial climate change. Nature Geoscience 3, 110113.Google Scholar
Wang, P., Wang, B., Cheng, H., Fasullo, J., Guo, Z., Kiefer, T., Liu, Z., 2014. The global monsoon across timescales: coherent variability of regional monsoons. Climate of the Past 10, 21632291.Google Scholar
Wang, Q., Wang, Y., Shao, Q., Liang, Y., Zhang, Z., Kong, X., 2018. Millennial-scale Asian monsoon variability during the late Marine Isotope Stage 6 from Hulu Cave, China. Quaternary Research 75, 112.Google Scholar
Wang, Y., Cheng, H., Edwards, R.L., An, Z., Wu, J., Shen, C., Dorale, J.A., 2001. A high-resolution absolute-dated late Pleistocene Monsoon record from Hulu Cave, China. Science 294, 23452348.Google Scholar
Wang, Y., Cheng, H., Edwards, R.L., He, Y., Kong, X., An, Z., Wu, J., et al. , 2005. The Holocene Asian monsoon: links to solar changes and North Atlantic climate. Science 308, 854857.Google Scholar
Wang, Y., Cheng, H., Edwards, R.L., Kong, X., Shao, X., Chen, S., Wu, J., et al. , 2008. Millennial- and orbital-scale changes in the East Asian monsoon over the past 224,000 years. Nature 451, 10901093.Google Scholar
Wolff, E.W., Chappellaz, J., Blunier, T., Rasmussen, S.O., Svensson, A., 2010. Millennial-scale variability during the last glacial: the ice core record. Quaternary Science Reviews 29, 28282838.Google Scholar
Wu, H., Zhang, X., Li, X., Li, G., Huang, Y., 2015. Seasonal variations of deuterium and oxygen-18 isotopes and their response to moisture source for precipitation events in the subtropical monsoon region. Hydrological Processes 29, 90102.Google Scholar
Zhang, J., Jia, Y., Lai, Z., Long, H., Yang, L.H., 2011. Holocene evolution of Huangqihai Lake in semi-arid northern China based on sedimentology and luminescence dating. Holocene 21, 12611268.Google Scholar
Zhao, K., Wang, Y., Edwards, R.L., Cheng, H., Liu, D., 2010. High-resolution stalagmite δ18O records of Asian monsoon changes in central and southern China spanning the MIS 3/2 transition. Earth and Planetary Science Letters 298, 191198.Google Scholar
Zhao, K., Wang, Y., Edwards, R.L., Cheng, H., Liu, D., Kong, X., Ning, Y., 2016. Contribution of ENSO variability to the East Asian summer monsoon in the late Holocene. Palaeogeography, Palaeoclimatology, Palaeoecology 449, 510519.Google Scholar
Zhou, H., Chi, B., Michael, L., Zhao, J., Yan, J., Alan, G., Feng, Y., 2008. High-resolution and precisely dated record of weathering and hydrological dynamics recorded by manganese and rare-earth elements in a stalagmite from central China. Quaternary Research 69, 438446.Google Scholar
Figure 0

Figure 1. (Color online) Chronology of stalagmite HL161. (a) Scanned image and location of fifteen 230Th samples. (b) The linearly-interpolated age model, including error bars for 230Th ages.

Figure 1

Table 1. 230Th dating results for sample HL161 stalagmite from Hulu Cave, China.

Figure 2

Figure 2. Multi-proxy records for stalagmite HL161. (A) δ18O records for stalagmite HL161 (green curve; this study), MSL (blue curve; Wang et al., 2001), MSD (grey curve; Wang et al., 2001) and a new composite record (yellow curve; Cheng et al., 2016). 2-sigma error bars for the records are also shown with corresponding colors. The previously published Hulu records are plotted for comparison with the record of stalagmite HL161. (B–E) are other multi-proxy records from stalagmite HL161. (B) δ13C. (C) Mg/Ca × 10−3. (D) Sr/Ca × 10−3. (E) Ba/Ca × 10−3. The bold lines in (C–E) are 13-point running averages of the raw data and show the long-term trends. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 3

Figure 3. Comparison of δ18O records in the monsoonal regions with Greenland paleotemperature. (A) Greenland NGRIP δ18Oice record (Andersen et al., 2004) plotted on the GICC05 timescale (Svensson et al., 2008). (B) δ18Ocalcite record from stalagmite HL161 (this study). (C) δ18Ocalcite record from stalagmite M1-2 (Burns et al., 2003). (D) Reconstructed water vapor source. △18O is calculated from the HL161 and M1-2 records and expressed as z-scores. Inferred periods of nearby water vapor sources (z-score △18O > 0.53) are shaded with standard deviations denoted by red dashed lines. The grey dashed line indicates a “rebound event” detected in the NGRIP δ18O ice record (Capron et al., 2010). The blue dashed line denotes an overall cooling trend during GI 12, and the green and orange dashed lines indicate stable monsoon climates. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 4

Figure 4. Results of wavelet analysis of the δ18O record from stalagmite HL161. The horizontal dashed line indicates the ~250-yr band prominent during interstadials. Spectral power is shown by colors ranging from deep blue (weak) to deep red (strong). The 95% significance level against red noise is shown by the thick solid line. The spectra were estimated using the method of Grinsted et al. (2004; http://noc.ac.uk/using-science/crosswavelet-wavelet-coherence). (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)

Figure 5

Figure 5. (A and B) Results of a replication test of δ13C records from Hulu Cave. (A) δ13C record for stalagmite HL161. (B) δ13C records from stalagmite MSL (brown) and stalagmite MSD (gold; Kong et al., 2005). (C and D) Results of an analysis of trace element from Hulu Cave. (C) Relationship between ln(Mg/Ca) versus ln(Sr/Ca) for stalagmite HL161 with the slope of the best-fit regression line indicated. (D) Comparison of values of Sr/Ca and Mg/Ca for bedrock (black dots), cave water (red dots), and stalagmite HL161 (purple crosses). The cave water and stalagmite samples are depleted in Mg and Sr relative to Ca in comparison to bedrock. (For interpretation of the references to color in this figure legend, the reader is referred to the web version of this article.)