Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

An optimized chronology for a stalagmite using seasonal trace element cycles from Shihua Cave, Beijing, North China

An optimized chronology for a stalagmite using seasonal trace element cycles from Shihua Cave,... www.nature.com/scientificreports OPEN An optimized chronology for a stalagmite using seasonal trace element cycles from Shihua Cave, Received: 12 December 2017 Beijing, North China Accepted: 28 February 2018 Published: xx xx xxxx 1,2 2 3 4 5 Fengmei Ban , Andy Baker , Christopher E. Marjo , Wuhui Duan , Xianglei Li , 1 2 3 4 2 Jinxian Han , Katie Coleborn , Rabeya Akter , Ming Tan & Gurinder Nagra Stalagmites play an important role in paleoclimatic reconstructions from seasonal to orbital time scales as Th-dating can provide an accurate absolute age. Additionally, seasonal trace element and optical layers can provide a precise age. We analyzed the seasonal variability of multiple trace elements on a stalagmite (XMG) in Shihua Cave, Beijing and compared them with results from laminae counting. The results show that (1) the polished section of the topmost part of XMG has obvious bi-optical layers under a conventional transmission microscope, however, laminae are not observed using this method in the rest of the sample, and (2) The variations of P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca show seasonal cycles throughout the sample. The PC1 in the Principal Component Analysis (PCA) of five trace elements represents the annual cycle. This stalagmite was deposited over 150 ± 1 years through PC1 peak counting. This result corresponds well with the annual layers and U-Th dating. Trace element cyclicity of PC1 can increase the accuracy of stalagmite dating, especially in the absence of obvious laminae and are a powerful method to identify seasonal changes in a strongly contrasting wet-dry monsoon climate region. Stalagmites are increasingly important for paleoclimatic reconstruction because of their wide distribution, mul- tiple proxies, high resolution and precise dating. Especially, the timescale over which it is possible to build paleo- 1,2 climate reconstructions can be from seasonal to orbital timescales . Trace element variations of stalagmites have 1,3–11 been used to explore sub-annual climatic and cave environmental signals and for chronology building . With the development of in-situ measurements, high resolution sampling methods such as secondary ion mass spec- 1,3,8–10,12 4–6,11,13,14 trometry (SIMS) , laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) and 15,16 micro-X-ray uo fl rescence imaging have been used to measure trace elements to get annual and even seasonal signals, which can be used to explore sub-annual climatic, cave environmental signals and chronology building. Because of a seasonally changing climate, such as wet and dry seasons, and associated seasonal changes in 17–19 hydrological or ventilation conditions within the cave , some stalagmites show obvious annual laminae under visible-light or UV-u fl orescence, which can be used to obtain an accurate age by counting laminae and compar - 20–25 ison to radiometric methods . Some stalagmites fail to exhibit visible or UV-fluorescence laminae, or they lack laminae in some sections, but the same stalagmites may show cycles in their trace element composition . The impact of seasonality on speleothems would be expected to be strong enough to generate annual chemical 26,27 10 variations, whether or not conventional growth laminae are visible . In 2009, Smith et al. built an objective and automated technique for establishing trace element chronologies analysing individual elements in stalagmites from two Alpine caves. en Th Nagra et al. exploited seasonal variations in trace elements to construct chronolo- gies in a Mediterranean climate employing principal components analysis (PCA). These studies imply that cyclic- ity of trace elements can be used to obtain an accurate age for non-laminated speleothem samples. Faculty of Environmental Economics, Shanxi University of Finance & Economics, Taiyuan, 030006, China. 2 3 Connected Waters Initiative Research Centre, UNSW Sydney, Sydney, NSW, 2052, Australia. Mark Wainwright Analytical Centre, UNSW Sydney, Sydney, NSW, 2052, Australia. Key Laboratory of Cenozoic Geology and Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, 100029, Beijing, China. Institute of Global Environmental Change, Xi’an Jiaotong University, Xi’an, 710049, China. Correspondence and requests for materials should be addressed to F.B. (email: banfm@163.com) Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 1 www.nature.com/scientificreports/ Figure 1. Profile map of Shihua Cave (aer L ft iu Hong et al ., 2013) and the location of stalagmite XMG (red star). Figure used by permission from Cai Binggui. In Shihua Cave, stalagmites typically have obvious laminae observed under visible-light and UV-fluorescence , and the annual laminae thickness, the isotope of oxygen and carbon and trace elements have 25,29 been used to reconstruct past climate and paleoenvironmental signals , including a 2650-year annual resolu- tion record of warm season temperature . Observations of drip water and calcite deposition in this cave shows that opaque layers (the fluorescent layer) was formed by the pulses of soil organic matter in the rainy season, while 19,30 the thick calcite layer was formed in the dry season . Stalagmite-XMG was growing when it was collected in 2015. U-Th dating has demonstrated that it grew for around the last 145 years . Under the confocal fluorescent microscope, it shows clear fluorescent cycles throughout the stalagmite transect, and 143 ± 7 cycles were counted . Its oxygen isotopes have been interpreted as recording East Asian Summer Monsoon (EASM) variability over the last 145 years . In this study, we aim to establish the optimal chronology for stalagmite XMG by comparison between the PCA of multiple trace ele- ment ratios (P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca) and the stalagmite fluorescent laminae. These comparisons allow us to evaluate the accuracy of different dating techniques and test whether the trace elements can supply the similarly accurate chronology. It is important to build stalagmite chronology by using annual trace element variability of multiple trace elements in a monsoonal climate. Additionally, we investigate whether individual climatically-sensitive elements in the stalagmite, such as P and Sr, can relate to dry and wet signals and have the potential for paleoclimate reconstruction. Site Description Shihua Cave (115°56′E, 39°47′N, 251 m above sea level at the entrance) is located in Baihua Hill, Fangshan county, about 50 km away from southwest Beijing. Shihua Cave is situated in the Middle Ordovician carbonate rocks of the Maojiagou Formation, which mainly consists of limestone, with some interlayered dolomite. This limestone strikes in an east-west direction and is inclined towards the south at about 30°. The thickness of the bedrock over - lying the cave varies between 30 and 130 m. Coal layers are unevenly distributed in the overlying bedrock. Today, the natural vegetation above the cave is dominated by shrubs, while persimmon and walnut trees are grown in the valley. The overlying soil is brown soil and its thickness varies from 0.6 to 1.0 m in the valley and 0 to 0.5 m on the hill slopes. The soil above the cave consists of a 5–15 cm humic layer that has a high content of total organic −1 carbon (40–70 g·kg ). Situated within the East Asian monsoon zone, the Shihua Cave region typically has cold/dry winters and warm/wet summers. The East Asian monsoon is the major source of precipitation. In Beijing more than 70% of the total annual precipitation falls during the summer monsoon months (June to August), which are also the warmest months. The dry period lasts for at least 8 months from October to May of the next year. From 1971– 2000 CE, the mean annual temperature was 12.3 °C while mean annual precipitation was 572 mm. Shihua Cave has a multi-level and multi-branched structure. Based on the newest survey, the cave has five levels (Fig. 1). The passages open to visitors (first to third levels) are about 2500 m long. The XMG stalagmite, with a diameter of 40 mm and height of 36 mm, was located in the transition level between the level 2 and 3 (See Fig. 1). It was confirmed that it was actively growing when we collected it in October 2015 through the monitoring of calcite deposition on a glass plate for half a year. A total of 143 annual uo fl rescent layers have been previously counted with a seven year error at the bottom, consistent within the mar - 230 25 gin of error of the Th analyses . Results XMG stalagmite laminae. The XMG stalagmite comprises calcite that has obvious layers under conven- tional visible-light and mercury-source ultra-violet transmission microscopy, but only in the topmost part. From the top to 3.8–4 mm, the variability in depth is due to uneven growth layer thickness(Fig. 2).These layers are 21,28,29 similar to other stalagmites from the cave , where the laminae have been described as bi-optical (the inter- 21,28 faces seem dark under transmitted light and luminescent under reflected and ultra-violet light) . There are no layers apparent using visible-light or mercury-source ultra-violet transmission microscopy in the rest of XMG using the same method. There is also a significant difference in the character of the laminae in the topmost part of XMG and previously published stalagmites in Shihua Cave. Stalagmite XMG has a thick opaque layer and a thin transparent layer under conventional transmission microscopy, while previously published stalagmites records Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 2 www.nature.com/scientificreports/ Figure 2. Stalagmite XMG images collected during this study showing the sampling transects used for LA- ICP-MS analysis and (inset expanded on the right) micrograph of stalagmite XMG showing the contrast in topmost layers revealed by white-light illumination. 19,21,29 comprised thin dark layers and thick transparent layers . In most stalagmites from the cave, the thin inter- faces of the laminae consist of opaque calcite with the presence of organic matter. In stalagmite XMG, using mercury-source ultra-violet transmission microscopy, a total of 23 visible and flu- orescent layers were counted in the top 3.8–4.0 mm, with an accumulated error of 1 year. Assuming that these are annual, this indicates that the deposited calcite fabric on XMG has changed since the early 1990s. This could be related to the cave environment, as an artificial exit was created in 1987 when the cave was opened to tourists (see Fig. 1). Although no data exists to verify this hypothesis, we suggest this exit has likely changed ventilation and cave air CO concentrations. The XMG stalagmite was located on the transition between levels 2 and 3, where the passage is narrow. Ventilation patterns are likely to have changed here due to seasonal changes in ventilation between the two the entrances at different elevations induced by the chimney effect . Under the confocal fluorescence microscope, stimulated with a 40 mW laser (488 nm wavelength), and using a filter which isolated the emitted wavelengths at 515/30 nm—the green portion of the visible spectrum, it has been previously shown that clear fluorescent laminae are present throughout the whole stalagmite . This instrument has better resolution, image sectioning and signal-to-noise ratio than the traditional fluorescent microscope . It is this confocal fluorescence microscope data that is used in comparison with trace element data. Trace element chronology. Figure 3 shows the trace element time-series for P/Ca, Sr/Ca, Ba/Ca U/Ca, and Mg/Ca, using a 17-point Savitzky-Golay smoothing (Supplementary Fig. S1 includes raw data). We observe similar annual cycles in all the elements (Fig. 3). Table 1 shows that Sr/Ca is inversely correlated to P/Ca (Table 1). When Sr/Ca and P/Ca are compared with the visible-light lamina image, we can see that the peaks of Sr/Ca corre- spond with the white layers and the peaks of P/Ca match closely with the opaque layers (Fig. 4). In Shihua Cave, the white layer (low-organic calcite) deposits on the stalagmites in late winter and spring when the climate is dry, 19,30 then the opaque layer (organic-rich calcite) is formed in the high rainfall monsoon period . This can confirm P is a soil-derived element, which increases in the high rainfall monsoon period. From Table 2, we can see there is close agreement in the number of peaks between the transect 2 and 3 (155 and 153 peaks respectively). In these two transects, PC1 has main contributions from Sr/Ca and P/Ca (Table 3). In contrast, transect 1 has 169 peaks for PC1 representing the changes in U/Ca and P/Ca. Transect 1 also has a higher number of peaks of Sr/Ca, and appears to include more sub-annual peaks compared to the other transects. By ana- lyzing three tracks, we can judge that the PC1 in the transect 2 and 3 are comparable and can represent the annual cycle and be used to build the chronology. This result demonstrates the advantage of measuring at least three tran- sects to enable rejection of transects that are less representative of annual trace element periodicity. Ae ft r using the auto-counting method, further screening of the peaks was done by comparing the transects 2 and 3. Eight PC1 peaks in transect 2 and seven in transect 3 were removed manually (see Methods, Supplementary Fig. S2). Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 3 www.nature.com/scientificreports/ Figure 3. e va Th riation of P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca in the three transects. The similar annual cycles of smoothing lines are showed (17 point window Savitzky-Golay method). P/Ca Sr/Ca Ba/Ca Mg/Ca U/Ca Transect-1 P/Ca 1 Sr/Ca −0.36* 1 Ba/Ca −0.01 0.70* 1 Mg/Ca 0.45* −0.07* 0.21* 1 U/Ca 0.55* 0.02 0.35* 0.43* 1 Transect-2 P/Ca 1 Sr/Ca −0.55* 1 Ba/Ca −0.32* 0.64* 1 Mg/Ca 0.34* −0.33* 0.11* 1 U/Ca 0.34* −0.24* −0.07* 0.23* 1 Transect-3 P/Ca 1 Sr/Ca −0.59* 1 Ba/Ca −0.28* 0.58* 1 Mg/Ca 0.42* −0.39* 0.17* 1 U/Ca 0.51* −0.48* −0.13* 0.45* 1 Table 1. e P Th earson correlation coefficient of five trace elements in three transects. *Correlation is significant using 2-tailed test at the 0.05 level. Aer t ft hese peaks were removed, 147 peaks and 146 peaks are preserved respectively in transect 2 and 3. Because the first peak of trace elements measurement is on the fourth laminae from top (Fig.  5b), the XMG stalagmite grew for 150 years with a 1 year error based in its trace elements chronology To verify the chronology of the trace elements in the stalagmite, the peaks of PC1 in the transect 2 and 3 are compared with the fluorescent microscope image of the laminae from Li et al . (Fig. 5a). In the topmost part of the figure, we compare the transparent image with PC1 (Fig.  5b). Using microscopy, a total of 23 annual layers were counted with an error of 1 year. Due to the start position of the laser ablation transects, the first trace element peaks occur on the fourth lamina. The elemental chronology shows 19 and 18 peaks respectively in transects 2 and 3, which matches the microscopy layer counting very well. The absolute position of peaks in the two transects has a small dier ff ence because of uneven deposition rates across the surface (Fig.  5 purple ellipse). Through this comparison, the annual nature of the trace elements can be verified from their good agreement. Importantly, Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 4 www.nature.com/scientificreports/ Figure 4. Comparison between the peaks of Sr/Ca and P/Ca and layers in transect 3. The vertical dashed lines represent the transparent laminae formed in dry seasons (the blue dashed line indicates uncertain lamina). Note the y-axis of Sr/Ca is inverted, so that the peaks of Sr/Ca in the graph corresponds to the lower value of Sr/ Ca. Stalagmite XMG images were collected during this study. PC1 PC2 P/Ca Sr/Ca Transect-1 169 165 166 166 Transect-2 155 194 165 153 Transect-3 153 186 161 153 Table 2. e r Th esults of auto-counting peak based on three transects. P/Ca Sr/Ca Ba/Ca U/Ca Mg/Ca variance% PC1 0.772 −0.030 0.400 0.847 0.759 41.021 transect1 PC2 −0.422 0.940 0.846 0.070 −0.058 35.718 PC1 0.798 −0.875 −0.631 0.475 0.499 45.513 transect2 PC2 0.133 0.238 0.699 0.675 0.444 24.302 PC1 0.821 −0.865 −0.468 0.601 0.744 51.126 transect3 PC2 0.045 0.275 0.833 0.652 0.267 25.346 Table 3. e lo Th adings of each element onto the first two PCs from the PCA. PCs are listed which explain around 70% of the variance in each transect. the measurement of trace elements is consistent even in places without obvious laminae (Fig. 5 white square). er Th efore, the trace elements not only confirm the annual nature of the laminae but also provide a good alterna- tive for stalagmites without obvious laminae. Despite small variations in the absolute peak location determined by trace element peak counting and micros- copy, the distance between peaks can be used to determine an annual growth rate. In stalagmite XMG, the growth rate is faster below 25 mm from the top (0.26 mm/year) than above it (0.22 mm/year). The cumulative annual growth rate data can be used to derive age-depth relationships. Figure 6 shows that the age models of PC1 for the two transects are within the error of dating by U-Th and visible lamina counting. The PC1 variations in two transects are not synchronous, a result of uneven deposition perpendicular to the growth axis. Discussion The regular alternation of dry and wet season climate, and the seasonal flush of rain can control the hydro- logical and hydrochemical properties of the drip water that feeds stalagmites, and can cause the annual cycles of trace elements. Seasonal variations have previously been observed in the trace elements in stalagmites from 4,11,12 Mediterranean and Alpine climate areas, which also have a strong seasonality . In monsoonal climates, trace elements show seasonal cycles because of the dominant seasonal control of hydrology, which ae ff cts infiltration and prior calcite precipitation (PCP). During dry periods, the residence of recharge water in a karst aquifer can be long enough to result in high trace element-to-calcium ratios such as Sr and Ba because of PCP . In Shihua Cave, less than 30% annual rainfall is in the dry period. During this period, most drip water comes from water which was recharged to the unsaturated zone in the preceding rain season, and due to PCP has higher Sr [Ba]/Ca. This contrasts with tropical climates, such as at Gunung Mulu and Gunung Buda National Parks in northern Borneo, where stalagmite Mg [Sr]/Ca ratios are poor indicators of hydroclimate conditions, due to a wetter climate and less strong seasonality . Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 5 www.nature.com/scientificreports/ Figure 5. e co Th mparison between the laminae images (revised from Li et al . 2017 by recoloring and adjusting the brightness and contrast of image using the Microsoft Powerpoint 2010) and the two PC1. (Purple ellipses means the location of peaks are different between the two transects; white square area represent the locations without clear fluorescent layers). The numbers of b) are calendar age in which the layers were formed. Figure 6. Age model of stalagmite XMG. Phosphorous, as a nutrient element, is associated with soil and vegetation decay on the surface. It has been 12,26,36 shown to have an association with colloid-associated metals e.g. Zn, Pb in cave drip waters . At Shihua Cave, in the wet monsoon period, an increase in P/Ca is observed. The increase in P/Ca observed in stalagmite XMG, is indicative of a flush of soil-derived water in the summer monsoon. In Table  2 the P/Ca peaks are slightly greater than the number of PC1 peaks, suggesting that in some instances a sub-annual flush of soil derived material is possible. In Shihua Cave, field monitoring showed that intra-annual variability in drip discharge was significantly higher than inter-annual variability. The flushing peak of soil organic matter occurs above a threshold rainfall intensity, and can have several annual pulses corresponding with intense rain events . The fact that the amount of P/Ca peaks is more than the annual cycles can be explained by multiple flushes in a year. Therefore, P/Ca can be as the signal of the long-duration and/or high-magnitude precipitation events. Sr/Ca is strongly negatively correlated with the P/Ca in the XMG stalagmite and positively with Ba/Ca in 10,12 Shihua Cave (Table 1) and this has been observed in previous studies . e Th strong similarity between the var - iation of Sr/Ca and Ba/Ca, as also observed the Heshang Cave , suggests a common control mechanism, namely PCP. Sr and Ba come from the overlying limestone, and the unsaturated zone water residence time may also be important for these elements on inter-annual scales. Mg may have multiple sources in Shihua Cave. At the site the soil is derived both from bedrock weathering and wind-blown material. The latter includes clay which can be a significant source of Mg, as has been observed Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 6 www.nature.com/scientificreports/ in regions of wind-blown dust deposition . Soil mineral derived elements therefore have a complex origin that can change over time, as it can be influenced by the degree of clay weathering (temperature and humidity), its dissolution rate, as well as unsaturated zone water residence time (PCP). Mg/Ca, with a positive relation with U/ Ca, P/Ca, and an anti-correlation with Sr/Ca, shows it has a closer relation with soil than bedrock. In stalagmite XMG, U is complicated to interpret as it is organic-transported, but has to compete with other organic and colloidal bound elements. It may also come from both a soil and bedrock source. It is a complicated element for paleoclimate interpretation at our site to explain. The further monitoring field work should be carried to confirm its transferring mechanism. The strong annual variation of key trace elements, can provide chronological markers for high-resolution studies of the climate proxies. Regular and continuous annual variations in trace element composition for PC1 and Sr/Ca in stalagmite XMG are apparent. Sufficient mixing of infiltration waters in the karst before reaching the cave must occur, which smooth much of sub-annual hydrological variability but retains sufficient seasonal varia- bility. Stalagmite XMG, situated at a depth of ~60 m below the surface, has a hydrological connectivity sufficient to resolve this seasonality in trace elements. Seasonality in trace elements in drip waters and speleothems have 11 12 4 been observed at shallower sites e.g. ~4 m below surface , ~20 m below surface and ~30 m below surface . This contrasts to other deeper sites, where no annual variation of trace elements has been observed e.g. over 1000 m depth at Corchia Cave of central Italy . Overall, the variations of five trace element-to-calcium ratios (P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca) show the seasonal cycle from the XMG stalagmite in Shihua Cave, located in East Asian monsoon zone. The PCA of the v fi e trace elements showed that the cycles of PC1 represents the annual cycle. A PCA approach is beneficial for identifying the best elements sensitive to PCP and the annual hydrological cycle, and these can then be used as a record of annual laminae in the China monsoon area. We recommend P and Sr as the least ambiguous for use in paleoclimate in the Chinese monsoon regions. Sr cycles are closest to the annual cycles, and P provides a record of soil-connectivity, which may also include sub-annual climate information. Importantly the annual laminae of trace elements cover the shortage of unclear layer-counting. Therefore, the trace elements not only confirm the annual nature of the laminae but also provide a good alternative for stalagmites without obvious laminae. Methods Trace element measurements. e XM Th G stalagmite thick-section for trace elements analysis was cut with a wire line saw along its growth axis. The LA-ICP-MS analysis was completed at the Mark Wainwright Analytical Centre of the University of New South Wales, using NWR213 Laser Ablation sample accessory (ESI, USA) and a NexION 300D ICP-MS (PerkinElmer, USA). e Th laser wavelength was 213 nm operating at a frequency of 10 Hz at 70% power, equivalent to a fluence of 11.83 J/cm , focused down to 80 μm on the sample, with a scan speed of 10 μm/sec along the transect. Two standards, NIST612 (Sr: 78.4, Ba: 41.0, P: 39.1, Mg: 85.09, and U: 37.38ppm) and NIST610 (Sr: 515.5, Ba: 453, Mg: 457 and U: 461.5 ppm) were used to calibrate the data before and aer e ft ach scan transect. To account for uneven growth, three parallel transects separated by 1 mm (Fig. 2) along the central growth axis of the sample were chosen to reduce within stalagmite variability in our measurements. Trace elements cycles counting. All trace data were ratioed to the calcium signal and standardised. The Pauta criterion was used to remove the outliers that coincided with small holes on the sample surface. The missing data was then interpolated using the average of neighbouring data. Principal component analysis (PCA) can be used to identify components that explain the variability of all trace elements along the transect. Variation of the first component score along the transect should be related to the annual periodicity. The PCA/peak counting method for building chronologies will be most effective if only those elements that have the clearest cycles are included . We chose the five trace elements (Sr/Ca, P/Ca, Mg/ Ca, U/Ca and Ba/Ca) among the stalagmites which showed the annual variability then employed PCA using free software named Past3.x (download from https://folk.uio.no/ohammer/past/ ) in each transect to calculate the first principal component (PC1) that was associated with annual changes in trace element concentration. e s Th cores plots for PC1 along the transect was smoothed using the Savitzky-Golay method using a 17-point window then an automated peak counting program, ‘Peak Analyser’ in Originlab12 was used to locate peaks. A 25 point filter window was selected based the half of growth speed (~0.23 mm/year) to avoid missing annual peaks. To the automated peak counting process was refined manually by identifying cross-calibration between two transects by applying the criteria, peaks that only emerged in one transect with a height less 30% of the neigh- bouring relative peaks were removed. References 1. Orland, I. J. et al. Seasonal climate signals (1990–2008) in a modern Soreq Cave stalagmite as revealed by high-resolution geochemical analysis. Chem. Geol. 363, 322–333 (2014). 2. Cheng, H. et al. The Asian monsoon over the past 640,000 years and ice age terminations. Nature. 534, 640–646, https://doi. org/10.1038/nature18591 (2016). 3. Roberts, M. S., Smart, P. L. & Baker, A. Annual trace element variations in a Holocene speleothem. Earth Planet. Sci. Lett. 154, 237–246 (1998). 4. Treble, P., Shelley, J. M. G. & Chappell, J. Comparison of high resolution sub-annual records of trace elements in a modern (1911–1992) speleothem with instrumental climate data from southwest Australia. Earth Planet. Sci. Lett. 216, 141–153 (2003). 5. Treble, P., Chappell, J., Gagan, M., Mckeegan, K. & Harrison, T. In situ measurement of seasonal δ18O variations and analysis of isotopic trends in a modern speleothem from southwest Australia. Earth Planet. Sci. Lett. 233, 17–32 (2005). 6. Johnson, K. R., Hu, C., Belshaw, N. S. & Henderson, G. M. Seasonal trace-element and stable-isotope variations in a Chinese speleothem: The potential for high-resolution paleomonsoon reconstruction. Earth Planet. Sci. Lett. 244, 394–407 (2006). 7. Liu, D. et al. A detailed comparison of Asian Monsoon intensity and Greenland temperature during the Allerød and Younger Dryas events. Earth Planet. Sci. Lett. 272, 691–697 (2008). Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 7 www.nature.com/scientificreports/ 8. Orland, I. J. et al. Climate deterioration in the eastern Mediterranean as revealed by ion microprobe analysis of a speleothem that grew from 2.2 to 0.9 ka in Soreq Cave, Israel. Quat. Res. 71, 27–35 (2009). 9. Orland, I. J. et al. Seasonal resolution of Eastern Mediterranean climate change since 34 ka from a Soreq Cave speleothem. Geochim. Cosmochim. Acta. 89, 240–255 (2012). 10. Smith, C. L. et al. Chronology building using objective identification of annual signals in trace element profiles of stalagmites. Quatern. Geochronol. 4, 11–21 (2009). 11. Nagra, G. et al. Dating stalagmites in mediterranean climates using annual trace element cycles. Sci. Rep. 7, 1–12, https://doi. org/10.1038/s41598-017-00474-4 (2017). 12. Huang, Y. et al. Seasonal variations in Sr, Mg and P in modern speleothems. Chem. Geol. 175, 429–448 (2001). 13. Jochum, K. P. et al. Accurate trace element analysis of speleothems and biogenic calcium carbonates by LA-ICP-MS. Chem. Geol. 318–319, 31–44 (2012). 14. Finch, A. A., Shaw, P. A., Weedon, G. P. & Holmgren K. Trace element variation in speleothem aragonite: potential for palaeoenvironmental reconstruction. Earth Planet. Sci. Lett. 186, 255–267 (2001). 15. Frisia, S. et al. A re-evaluation of the palaeoclimatic significance of phosphorus variability in speleothems revealed by high- resolution synchrotron micro XRF mapping. Clim. Past 8, 2039–2051, https://doi.org/10.5194/cp-8-2039-2012 (2012). 16. Wu, J. Y. et al. Stable isotope and trace element investigation of two contemporaneous annually-laminated stalagmites from northeastern China surrounding the “8.2 ka event”. Clim. Past 8, 1497–1507, https://doi.org/10.5194/cp-8-1497-2012 (2012). 17. Frisia, S., Borsato, A., Fairchild, I. J. & McDermott, F. Calcite fabrics, growth mechanisms, and environments of formation in speleothems from the Italian Alps and South-western Ireland. J. Sediment. Res. 70, 1183–1196 (2000). 18. Spӧtl, C., Fairchild, I. J. & Tooth, A. F. Cave air control on dripwater geochemistry, Obir caves (Austria): implications for speleothem deposition in dynamically ventilated caves. Geochim. Cosmochim. Acta. 69, 2451–2468 (2005). 19. Cai, B., Zhu, J., Ban, F. & Tan, M. Intra-annual variation of the calcite deposition rate of drip water in Shihua Cave, Beijing, China and its implications for palaeoclimatic reconstructions. Boreas. 40(3), 525–535, https://doi.org/10.1111/j.1502-3885.2010.00201.x (2011). 20. Baker, A., Smart, P. L., Edwards, R. L. & Richards, D. Annual growth banding in a cave stalagmite. Nature. 364, 518–520 (1993). 21. Tan, M. et al. Cyclic rapid warming on centennial-scale revealed by a 2650-year stalagmite record of warm season temperature. Geophysical Research Letters. 30(12), 1617–1620, doi:10D1029/2003GL017352D2003 (2003). 22. Ma, Z. et al. Timing and structure of the Younger Dryas event in northern China. Quat. Sci. Rev. 41, 83–93 (2012). 23. Liu, Y. et al. Links between the East Asian monsoon and North Atlantic climate during the 8,200 year event. Nature Geosci. 6, 117–120 (2013). 24. Duan, W. et al. Onset and duration of transitions into Greenland Interstadials 15.2 and 14 in northern China constrained by an annually laminated stalagmite. Sci. Rep. 6, 20844, https://doi.org/10.1038/srep20844 (2016). 25. Li, X. et al. The East Asian summer monsoon variability over the last 145 years inferred from the Shihua Cave record, North China. Sci. Rep. 7, 7078, https://doi.org/10.1038/s41598-017-07251-3 (2017). 26. Fairchild, I. J. et al. Annual to sub-annual resolution of multiple trace-element trends in speleothems. J. Geol. Soc. 158, 831–841 (2001). 27. Fairchild, I. J., Frisia, S., Borsato, A. & Tooth, A. F. Speleothems. In Geochemical Sediments and Landscapes.(eds Nash, D. J., McLaren, S. J.) 200–245 (Blackwells, Oxford, 2007). 28. Tan, M. et al. Bioptical microcycles of laminated speleothems from China and their chronological significance. Chinese Science Bulletin. 44(17), 1604–1607 (1999). 29. Hou, J. Z., Tan, M., Cheng, H. & Liu, T. S. Stable isotope records of plant cover change and monsoon variation in the past 2200 years: evidence from laminated stalagmites in Beijing, China. Boreas. 32, 304–313 (2003). 30. Ban, F., Pan, G., Zhu, J., Cai, B. & Tan, M. Temporal and spatial variations in the discharge and dissolved organic carbon of drip waters in Beijing Shihua Cave, China. Hydrol. Process. 22, 3749–3758 (2008). 31. Liu, H. et al. Structure and evolution of the Shihua Cave in Beijing. Carsologica Sinica (in Chinese). 34(1), 27–34 (2015). 32. Covington, M. D. & Perne., M. Consider a cylindri-cal cave: A physicist’s view of cave and karst science. Acta Carsologic. 44/3, 363–380 (2015). 33. Sanderson, M. J., Smith, I., Parker, I. & Bootman, M. D. FluorescenceMicroscopy. Cold Spring Harb Protoc. 10, 1–36, https://doi. org/10.1101/pdb.top071795 (2014). 34. Fairchild, I. J. & Treble, P. C. Trace elements in speleothems as recorders of environmental change. Quat. Sci. Rev. 28, 449–468 (2009). 35. Partin, J. W., Cobb, K. M., Adkins, J. F., Tueu, A. A. & Clark, B. Trace metal and carbon isotopic variations in cave dripwater and stalagmite geochemistry from northern Borneo. Geochem. Geophy. Geosy. 14(9), 3567–3585, https://doi.org/10.1002/ggge.20215 (2013). 36. Borsato, A., Frisia, S., Fairchild, I. J., Somogyi, A. & Susini, J. Trace element distribution in annual stalagmite laminae mapped by micrometer-resolution X-ray fluorescence: implications for incorporation of environmentally significant species. Geochim. Cosmochim. Acta. 71, 1494–1512 (2007). 37. Rutlidge, H. et al. Dripwater organic matter and trace element geochemistry in a semi-arid karst environment: Implications for speleothem paleoclimatology. Geochim. Cosmochim. Acta. 135, 217–230 (2014). 38. Regattieri, E. et al. Lateglacial to Holocene trace element record (Ba, Mg, Sr.) from Corchia Cave (Apuan Alps, central Italy): paleoenvironmental implications. J. Quat. Sci. 29(4), 381–392, https://doi.org/10.1002/jqs.2712 (2014). Acknowledgements Many thanks to Binggui Cai for his permission for using map of Shihua Cave, to Haiying Sun, Haijian Chen and Linhua Li at the administrative committee of Shihua Cave for field assistance, and to our collaboration, Jie Guo for her helpful discussions. This work was supported by National Natural Science Foundation of China (NSFC) (Grant No. 41402161) Natural Science Foundation of Shanxi Province, China (Grant No. 2015011085) and the Mark Wainwright Analytical Centre at UNSW, Australia. Author Contributions F.M.B., A.B. and C.E.M. conceived the project and wrote the first draft of the paper. W.H.D. contributed to the manuscript revision at different stages. F.M.B., K.C. and R.A. contribute to trace elements measurements. F.M.B., X.L.L. and J.X.H. did the field and monitoring work. F.M.B. and X.L.L. collected the XMG stalagmite sample. G.N. contributes to guide suggestion of Peak Analyzer software. M.T. provided many suggestions of counting layers. Data collected in this study has been archived at the Mark Wainwright Analytical Centre at UNSW, Australia and is available on request from the authors. Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 8 www.nature.com/scientificreports/ Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-22839-z. Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 9 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

An optimized chronology for a stalagmite using seasonal trace element cycles from Shihua Cave, Beijing, North China

Loading next page...
 
/lp/springer_journal/an-optimized-chronology-for-a-stalagmite-using-seasonal-trace-element-citizZDzE6

References (72)

Publisher
Springer Journals
Copyright
Copyright © 2018 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
DOI
10.1038/s41598-018-22839-z
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN An optimized chronology for a stalagmite using seasonal trace element cycles from Shihua Cave, Received: 12 December 2017 Beijing, North China Accepted: 28 February 2018 Published: xx xx xxxx 1,2 2 3 4 5 Fengmei Ban , Andy Baker , Christopher E. Marjo , Wuhui Duan , Xianglei Li , 1 2 3 4 2 Jinxian Han , Katie Coleborn , Rabeya Akter , Ming Tan & Gurinder Nagra Stalagmites play an important role in paleoclimatic reconstructions from seasonal to orbital time scales as Th-dating can provide an accurate absolute age. Additionally, seasonal trace element and optical layers can provide a precise age. We analyzed the seasonal variability of multiple trace elements on a stalagmite (XMG) in Shihua Cave, Beijing and compared them with results from laminae counting. The results show that (1) the polished section of the topmost part of XMG has obvious bi-optical layers under a conventional transmission microscope, however, laminae are not observed using this method in the rest of the sample, and (2) The variations of P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca show seasonal cycles throughout the sample. The PC1 in the Principal Component Analysis (PCA) of five trace elements represents the annual cycle. This stalagmite was deposited over 150 ± 1 years through PC1 peak counting. This result corresponds well with the annual layers and U-Th dating. Trace element cyclicity of PC1 can increase the accuracy of stalagmite dating, especially in the absence of obvious laminae and are a powerful method to identify seasonal changes in a strongly contrasting wet-dry monsoon climate region. Stalagmites are increasingly important for paleoclimatic reconstruction because of their wide distribution, mul- tiple proxies, high resolution and precise dating. Especially, the timescale over which it is possible to build paleo- 1,2 climate reconstructions can be from seasonal to orbital timescales . Trace element variations of stalagmites have 1,3–11 been used to explore sub-annual climatic and cave environmental signals and for chronology building . With the development of in-situ measurements, high resolution sampling methods such as secondary ion mass spec- 1,3,8–10,12 4–6,11,13,14 trometry (SIMS) , laser-ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) and 15,16 micro-X-ray uo fl rescence imaging have been used to measure trace elements to get annual and even seasonal signals, which can be used to explore sub-annual climatic, cave environmental signals and chronology building. Because of a seasonally changing climate, such as wet and dry seasons, and associated seasonal changes in 17–19 hydrological or ventilation conditions within the cave , some stalagmites show obvious annual laminae under visible-light or UV-u fl orescence, which can be used to obtain an accurate age by counting laminae and compar - 20–25 ison to radiometric methods . Some stalagmites fail to exhibit visible or UV-fluorescence laminae, or they lack laminae in some sections, but the same stalagmites may show cycles in their trace element composition . The impact of seasonality on speleothems would be expected to be strong enough to generate annual chemical 26,27 10 variations, whether or not conventional growth laminae are visible . In 2009, Smith et al. built an objective and automated technique for establishing trace element chronologies analysing individual elements in stalagmites from two Alpine caves. en Th Nagra et al. exploited seasonal variations in trace elements to construct chronolo- gies in a Mediterranean climate employing principal components analysis (PCA). These studies imply that cyclic- ity of trace elements can be used to obtain an accurate age for non-laminated speleothem samples. Faculty of Environmental Economics, Shanxi University of Finance & Economics, Taiyuan, 030006, China. 2 3 Connected Waters Initiative Research Centre, UNSW Sydney, Sydney, NSW, 2052, Australia. Mark Wainwright Analytical Centre, UNSW Sydney, Sydney, NSW, 2052, Australia. Key Laboratory of Cenozoic Geology and Environment, Institute of Geology and Geophysics, Chinese Academy of Sciences, 100029, Beijing, China. Institute of Global Environmental Change, Xi’an Jiaotong University, Xi’an, 710049, China. Correspondence and requests for materials should be addressed to F.B. (email: banfm@163.com) Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 1 www.nature.com/scientificreports/ Figure 1. Profile map of Shihua Cave (aer L ft iu Hong et al ., 2013) and the location of stalagmite XMG (red star). Figure used by permission from Cai Binggui. In Shihua Cave, stalagmites typically have obvious laminae observed under visible-light and UV-fluorescence , and the annual laminae thickness, the isotope of oxygen and carbon and trace elements have 25,29 been used to reconstruct past climate and paleoenvironmental signals , including a 2650-year annual resolu- tion record of warm season temperature . Observations of drip water and calcite deposition in this cave shows that opaque layers (the fluorescent layer) was formed by the pulses of soil organic matter in the rainy season, while 19,30 the thick calcite layer was formed in the dry season . Stalagmite-XMG was growing when it was collected in 2015. U-Th dating has demonstrated that it grew for around the last 145 years . Under the confocal fluorescent microscope, it shows clear fluorescent cycles throughout the stalagmite transect, and 143 ± 7 cycles were counted . Its oxygen isotopes have been interpreted as recording East Asian Summer Monsoon (EASM) variability over the last 145 years . In this study, we aim to establish the optimal chronology for stalagmite XMG by comparison between the PCA of multiple trace ele- ment ratios (P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca) and the stalagmite fluorescent laminae. These comparisons allow us to evaluate the accuracy of different dating techniques and test whether the trace elements can supply the similarly accurate chronology. It is important to build stalagmite chronology by using annual trace element variability of multiple trace elements in a monsoonal climate. Additionally, we investigate whether individual climatically-sensitive elements in the stalagmite, such as P and Sr, can relate to dry and wet signals and have the potential for paleoclimate reconstruction. Site Description Shihua Cave (115°56′E, 39°47′N, 251 m above sea level at the entrance) is located in Baihua Hill, Fangshan county, about 50 km away from southwest Beijing. Shihua Cave is situated in the Middle Ordovician carbonate rocks of the Maojiagou Formation, which mainly consists of limestone, with some interlayered dolomite. This limestone strikes in an east-west direction and is inclined towards the south at about 30°. The thickness of the bedrock over - lying the cave varies between 30 and 130 m. Coal layers are unevenly distributed in the overlying bedrock. Today, the natural vegetation above the cave is dominated by shrubs, while persimmon and walnut trees are grown in the valley. The overlying soil is brown soil and its thickness varies from 0.6 to 1.0 m in the valley and 0 to 0.5 m on the hill slopes. The soil above the cave consists of a 5–15 cm humic layer that has a high content of total organic −1 carbon (40–70 g·kg ). Situated within the East Asian monsoon zone, the Shihua Cave region typically has cold/dry winters and warm/wet summers. The East Asian monsoon is the major source of precipitation. In Beijing more than 70% of the total annual precipitation falls during the summer monsoon months (June to August), which are also the warmest months. The dry period lasts for at least 8 months from October to May of the next year. From 1971– 2000 CE, the mean annual temperature was 12.3 °C while mean annual precipitation was 572 mm. Shihua Cave has a multi-level and multi-branched structure. Based on the newest survey, the cave has five levels (Fig. 1). The passages open to visitors (first to third levels) are about 2500 m long. The XMG stalagmite, with a diameter of 40 mm and height of 36 mm, was located in the transition level between the level 2 and 3 (See Fig. 1). It was confirmed that it was actively growing when we collected it in October 2015 through the monitoring of calcite deposition on a glass plate for half a year. A total of 143 annual uo fl rescent layers have been previously counted with a seven year error at the bottom, consistent within the mar - 230 25 gin of error of the Th analyses . Results XMG stalagmite laminae. The XMG stalagmite comprises calcite that has obvious layers under conven- tional visible-light and mercury-source ultra-violet transmission microscopy, but only in the topmost part. From the top to 3.8–4 mm, the variability in depth is due to uneven growth layer thickness(Fig. 2).These layers are 21,28,29 similar to other stalagmites from the cave , where the laminae have been described as bi-optical (the inter- 21,28 faces seem dark under transmitted light and luminescent under reflected and ultra-violet light) . There are no layers apparent using visible-light or mercury-source ultra-violet transmission microscopy in the rest of XMG using the same method. There is also a significant difference in the character of the laminae in the topmost part of XMG and previously published stalagmites in Shihua Cave. Stalagmite XMG has a thick opaque layer and a thin transparent layer under conventional transmission microscopy, while previously published stalagmites records Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 2 www.nature.com/scientificreports/ Figure 2. Stalagmite XMG images collected during this study showing the sampling transects used for LA- ICP-MS analysis and (inset expanded on the right) micrograph of stalagmite XMG showing the contrast in topmost layers revealed by white-light illumination. 19,21,29 comprised thin dark layers and thick transparent layers . In most stalagmites from the cave, the thin inter- faces of the laminae consist of opaque calcite with the presence of organic matter. In stalagmite XMG, using mercury-source ultra-violet transmission microscopy, a total of 23 visible and flu- orescent layers were counted in the top 3.8–4.0 mm, with an accumulated error of 1 year. Assuming that these are annual, this indicates that the deposited calcite fabric on XMG has changed since the early 1990s. This could be related to the cave environment, as an artificial exit was created in 1987 when the cave was opened to tourists (see Fig. 1). Although no data exists to verify this hypothesis, we suggest this exit has likely changed ventilation and cave air CO concentrations. The XMG stalagmite was located on the transition between levels 2 and 3, where the passage is narrow. Ventilation patterns are likely to have changed here due to seasonal changes in ventilation between the two the entrances at different elevations induced by the chimney effect . Under the confocal fluorescence microscope, stimulated with a 40 mW laser (488 nm wavelength), and using a filter which isolated the emitted wavelengths at 515/30 nm—the green portion of the visible spectrum, it has been previously shown that clear fluorescent laminae are present throughout the whole stalagmite . This instrument has better resolution, image sectioning and signal-to-noise ratio than the traditional fluorescent microscope . It is this confocal fluorescence microscope data that is used in comparison with trace element data. Trace element chronology. Figure 3 shows the trace element time-series for P/Ca, Sr/Ca, Ba/Ca U/Ca, and Mg/Ca, using a 17-point Savitzky-Golay smoothing (Supplementary Fig. S1 includes raw data). We observe similar annual cycles in all the elements (Fig. 3). Table 1 shows that Sr/Ca is inversely correlated to P/Ca (Table 1). When Sr/Ca and P/Ca are compared with the visible-light lamina image, we can see that the peaks of Sr/Ca corre- spond with the white layers and the peaks of P/Ca match closely with the opaque layers (Fig. 4). In Shihua Cave, the white layer (low-organic calcite) deposits on the stalagmites in late winter and spring when the climate is dry, 19,30 then the opaque layer (organic-rich calcite) is formed in the high rainfall monsoon period . This can confirm P is a soil-derived element, which increases in the high rainfall monsoon period. From Table 2, we can see there is close agreement in the number of peaks between the transect 2 and 3 (155 and 153 peaks respectively). In these two transects, PC1 has main contributions from Sr/Ca and P/Ca (Table 3). In contrast, transect 1 has 169 peaks for PC1 representing the changes in U/Ca and P/Ca. Transect 1 also has a higher number of peaks of Sr/Ca, and appears to include more sub-annual peaks compared to the other transects. By ana- lyzing three tracks, we can judge that the PC1 in the transect 2 and 3 are comparable and can represent the annual cycle and be used to build the chronology. This result demonstrates the advantage of measuring at least three tran- sects to enable rejection of transects that are less representative of annual trace element periodicity. Ae ft r using the auto-counting method, further screening of the peaks was done by comparing the transects 2 and 3. Eight PC1 peaks in transect 2 and seven in transect 3 were removed manually (see Methods, Supplementary Fig. S2). Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 3 www.nature.com/scientificreports/ Figure 3. e va Th riation of P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca in the three transects. The similar annual cycles of smoothing lines are showed (17 point window Savitzky-Golay method). P/Ca Sr/Ca Ba/Ca Mg/Ca U/Ca Transect-1 P/Ca 1 Sr/Ca −0.36* 1 Ba/Ca −0.01 0.70* 1 Mg/Ca 0.45* −0.07* 0.21* 1 U/Ca 0.55* 0.02 0.35* 0.43* 1 Transect-2 P/Ca 1 Sr/Ca −0.55* 1 Ba/Ca −0.32* 0.64* 1 Mg/Ca 0.34* −0.33* 0.11* 1 U/Ca 0.34* −0.24* −0.07* 0.23* 1 Transect-3 P/Ca 1 Sr/Ca −0.59* 1 Ba/Ca −0.28* 0.58* 1 Mg/Ca 0.42* −0.39* 0.17* 1 U/Ca 0.51* −0.48* −0.13* 0.45* 1 Table 1. e P Th earson correlation coefficient of five trace elements in three transects. *Correlation is significant using 2-tailed test at the 0.05 level. Aer t ft hese peaks were removed, 147 peaks and 146 peaks are preserved respectively in transect 2 and 3. Because the first peak of trace elements measurement is on the fourth laminae from top (Fig.  5b), the XMG stalagmite grew for 150 years with a 1 year error based in its trace elements chronology To verify the chronology of the trace elements in the stalagmite, the peaks of PC1 in the transect 2 and 3 are compared with the fluorescent microscope image of the laminae from Li et al . (Fig. 5a). In the topmost part of the figure, we compare the transparent image with PC1 (Fig.  5b). Using microscopy, a total of 23 annual layers were counted with an error of 1 year. Due to the start position of the laser ablation transects, the first trace element peaks occur on the fourth lamina. The elemental chronology shows 19 and 18 peaks respectively in transects 2 and 3, which matches the microscopy layer counting very well. The absolute position of peaks in the two transects has a small dier ff ence because of uneven deposition rates across the surface (Fig.  5 purple ellipse). Through this comparison, the annual nature of the trace elements can be verified from their good agreement. Importantly, Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 4 www.nature.com/scientificreports/ Figure 4. Comparison between the peaks of Sr/Ca and P/Ca and layers in transect 3. The vertical dashed lines represent the transparent laminae formed in dry seasons (the blue dashed line indicates uncertain lamina). Note the y-axis of Sr/Ca is inverted, so that the peaks of Sr/Ca in the graph corresponds to the lower value of Sr/ Ca. Stalagmite XMG images were collected during this study. PC1 PC2 P/Ca Sr/Ca Transect-1 169 165 166 166 Transect-2 155 194 165 153 Transect-3 153 186 161 153 Table 2. e r Th esults of auto-counting peak based on three transects. P/Ca Sr/Ca Ba/Ca U/Ca Mg/Ca variance% PC1 0.772 −0.030 0.400 0.847 0.759 41.021 transect1 PC2 −0.422 0.940 0.846 0.070 −0.058 35.718 PC1 0.798 −0.875 −0.631 0.475 0.499 45.513 transect2 PC2 0.133 0.238 0.699 0.675 0.444 24.302 PC1 0.821 −0.865 −0.468 0.601 0.744 51.126 transect3 PC2 0.045 0.275 0.833 0.652 0.267 25.346 Table 3. e lo Th adings of each element onto the first two PCs from the PCA. PCs are listed which explain around 70% of the variance in each transect. the measurement of trace elements is consistent even in places without obvious laminae (Fig. 5 white square). er Th efore, the trace elements not only confirm the annual nature of the laminae but also provide a good alterna- tive for stalagmites without obvious laminae. Despite small variations in the absolute peak location determined by trace element peak counting and micros- copy, the distance between peaks can be used to determine an annual growth rate. In stalagmite XMG, the growth rate is faster below 25 mm from the top (0.26 mm/year) than above it (0.22 mm/year). The cumulative annual growth rate data can be used to derive age-depth relationships. Figure 6 shows that the age models of PC1 for the two transects are within the error of dating by U-Th and visible lamina counting. The PC1 variations in two transects are not synchronous, a result of uneven deposition perpendicular to the growth axis. Discussion The regular alternation of dry and wet season climate, and the seasonal flush of rain can control the hydro- logical and hydrochemical properties of the drip water that feeds stalagmites, and can cause the annual cycles of trace elements. Seasonal variations have previously been observed in the trace elements in stalagmites from 4,11,12 Mediterranean and Alpine climate areas, which also have a strong seasonality . In monsoonal climates, trace elements show seasonal cycles because of the dominant seasonal control of hydrology, which ae ff cts infiltration and prior calcite precipitation (PCP). During dry periods, the residence of recharge water in a karst aquifer can be long enough to result in high trace element-to-calcium ratios such as Sr and Ba because of PCP . In Shihua Cave, less than 30% annual rainfall is in the dry period. During this period, most drip water comes from water which was recharged to the unsaturated zone in the preceding rain season, and due to PCP has higher Sr [Ba]/Ca. This contrasts with tropical climates, such as at Gunung Mulu and Gunung Buda National Parks in northern Borneo, where stalagmite Mg [Sr]/Ca ratios are poor indicators of hydroclimate conditions, due to a wetter climate and less strong seasonality . Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 5 www.nature.com/scientificreports/ Figure 5. e co Th mparison between the laminae images (revised from Li et al . 2017 by recoloring and adjusting the brightness and contrast of image using the Microsoft Powerpoint 2010) and the two PC1. (Purple ellipses means the location of peaks are different between the two transects; white square area represent the locations without clear fluorescent layers). The numbers of b) are calendar age in which the layers were formed. Figure 6. Age model of stalagmite XMG. Phosphorous, as a nutrient element, is associated with soil and vegetation decay on the surface. It has been 12,26,36 shown to have an association with colloid-associated metals e.g. Zn, Pb in cave drip waters . At Shihua Cave, in the wet monsoon period, an increase in P/Ca is observed. The increase in P/Ca observed in stalagmite XMG, is indicative of a flush of soil-derived water in the summer monsoon. In Table  2 the P/Ca peaks are slightly greater than the number of PC1 peaks, suggesting that in some instances a sub-annual flush of soil derived material is possible. In Shihua Cave, field monitoring showed that intra-annual variability in drip discharge was significantly higher than inter-annual variability. The flushing peak of soil organic matter occurs above a threshold rainfall intensity, and can have several annual pulses corresponding with intense rain events . The fact that the amount of P/Ca peaks is more than the annual cycles can be explained by multiple flushes in a year. Therefore, P/Ca can be as the signal of the long-duration and/or high-magnitude precipitation events. Sr/Ca is strongly negatively correlated with the P/Ca in the XMG stalagmite and positively with Ba/Ca in 10,12 Shihua Cave (Table 1) and this has been observed in previous studies . e Th strong similarity between the var - iation of Sr/Ca and Ba/Ca, as also observed the Heshang Cave , suggests a common control mechanism, namely PCP. Sr and Ba come from the overlying limestone, and the unsaturated zone water residence time may also be important for these elements on inter-annual scales. Mg may have multiple sources in Shihua Cave. At the site the soil is derived both from bedrock weathering and wind-blown material. The latter includes clay which can be a significant source of Mg, as has been observed Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 6 www.nature.com/scientificreports/ in regions of wind-blown dust deposition . Soil mineral derived elements therefore have a complex origin that can change over time, as it can be influenced by the degree of clay weathering (temperature and humidity), its dissolution rate, as well as unsaturated zone water residence time (PCP). Mg/Ca, with a positive relation with U/ Ca, P/Ca, and an anti-correlation with Sr/Ca, shows it has a closer relation with soil than bedrock. In stalagmite XMG, U is complicated to interpret as it is organic-transported, but has to compete with other organic and colloidal bound elements. It may also come from both a soil and bedrock source. It is a complicated element for paleoclimate interpretation at our site to explain. The further monitoring field work should be carried to confirm its transferring mechanism. The strong annual variation of key trace elements, can provide chronological markers for high-resolution studies of the climate proxies. Regular and continuous annual variations in trace element composition for PC1 and Sr/Ca in stalagmite XMG are apparent. Sufficient mixing of infiltration waters in the karst before reaching the cave must occur, which smooth much of sub-annual hydrological variability but retains sufficient seasonal varia- bility. Stalagmite XMG, situated at a depth of ~60 m below the surface, has a hydrological connectivity sufficient to resolve this seasonality in trace elements. Seasonality in trace elements in drip waters and speleothems have 11 12 4 been observed at shallower sites e.g. ~4 m below surface , ~20 m below surface and ~30 m below surface . This contrasts to other deeper sites, where no annual variation of trace elements has been observed e.g. over 1000 m depth at Corchia Cave of central Italy . Overall, the variations of five trace element-to-calcium ratios (P/Ca, Sr/Ca, Ba/Ca, U/Ca and Mg/Ca) show the seasonal cycle from the XMG stalagmite in Shihua Cave, located in East Asian monsoon zone. The PCA of the v fi e trace elements showed that the cycles of PC1 represents the annual cycle. A PCA approach is beneficial for identifying the best elements sensitive to PCP and the annual hydrological cycle, and these can then be used as a record of annual laminae in the China monsoon area. We recommend P and Sr as the least ambiguous for use in paleoclimate in the Chinese monsoon regions. Sr cycles are closest to the annual cycles, and P provides a record of soil-connectivity, which may also include sub-annual climate information. Importantly the annual laminae of trace elements cover the shortage of unclear layer-counting. Therefore, the trace elements not only confirm the annual nature of the laminae but also provide a good alternative for stalagmites without obvious laminae. Methods Trace element measurements. e XM Th G stalagmite thick-section for trace elements analysis was cut with a wire line saw along its growth axis. The LA-ICP-MS analysis was completed at the Mark Wainwright Analytical Centre of the University of New South Wales, using NWR213 Laser Ablation sample accessory (ESI, USA) and a NexION 300D ICP-MS (PerkinElmer, USA). e Th laser wavelength was 213 nm operating at a frequency of 10 Hz at 70% power, equivalent to a fluence of 11.83 J/cm , focused down to 80 μm on the sample, with a scan speed of 10 μm/sec along the transect. Two standards, NIST612 (Sr: 78.4, Ba: 41.0, P: 39.1, Mg: 85.09, and U: 37.38ppm) and NIST610 (Sr: 515.5, Ba: 453, Mg: 457 and U: 461.5 ppm) were used to calibrate the data before and aer e ft ach scan transect. To account for uneven growth, three parallel transects separated by 1 mm (Fig. 2) along the central growth axis of the sample were chosen to reduce within stalagmite variability in our measurements. Trace elements cycles counting. All trace data were ratioed to the calcium signal and standardised. The Pauta criterion was used to remove the outliers that coincided with small holes on the sample surface. The missing data was then interpolated using the average of neighbouring data. Principal component analysis (PCA) can be used to identify components that explain the variability of all trace elements along the transect. Variation of the first component score along the transect should be related to the annual periodicity. The PCA/peak counting method for building chronologies will be most effective if only those elements that have the clearest cycles are included . We chose the five trace elements (Sr/Ca, P/Ca, Mg/ Ca, U/Ca and Ba/Ca) among the stalagmites which showed the annual variability then employed PCA using free software named Past3.x (download from https://folk.uio.no/ohammer/past/ ) in each transect to calculate the first principal component (PC1) that was associated with annual changes in trace element concentration. e s Th cores plots for PC1 along the transect was smoothed using the Savitzky-Golay method using a 17-point window then an automated peak counting program, ‘Peak Analyser’ in Originlab12 was used to locate peaks. A 25 point filter window was selected based the half of growth speed (~0.23 mm/year) to avoid missing annual peaks. To the automated peak counting process was refined manually by identifying cross-calibration between two transects by applying the criteria, peaks that only emerged in one transect with a height less 30% of the neigh- bouring relative peaks were removed. References 1. Orland, I. J. et al. Seasonal climate signals (1990–2008) in a modern Soreq Cave stalagmite as revealed by high-resolution geochemical analysis. Chem. Geol. 363, 322–333 (2014). 2. Cheng, H. et al. The Asian monsoon over the past 640,000 years and ice age terminations. Nature. 534, 640–646, https://doi. org/10.1038/nature18591 (2016). 3. Roberts, M. S., Smart, P. L. & Baker, A. Annual trace element variations in a Holocene speleothem. Earth Planet. Sci. Lett. 154, 237–246 (1998). 4. Treble, P., Shelley, J. M. G. & Chappell, J. Comparison of high resolution sub-annual records of trace elements in a modern (1911–1992) speleothem with instrumental climate data from southwest Australia. Earth Planet. Sci. Lett. 216, 141–153 (2003). 5. Treble, P., Chappell, J., Gagan, M., Mckeegan, K. & Harrison, T. In situ measurement of seasonal δ18O variations and analysis of isotopic trends in a modern speleothem from southwest Australia. Earth Planet. Sci. Lett. 233, 17–32 (2005). 6. Johnson, K. R., Hu, C., Belshaw, N. S. & Henderson, G. M. Seasonal trace-element and stable-isotope variations in a Chinese speleothem: The potential for high-resolution paleomonsoon reconstruction. Earth Planet. Sci. Lett. 244, 394–407 (2006). 7. Liu, D. et al. A detailed comparison of Asian Monsoon intensity and Greenland temperature during the Allerød and Younger Dryas events. Earth Planet. Sci. Lett. 272, 691–697 (2008). Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 7 www.nature.com/scientificreports/ 8. Orland, I. J. et al. Climate deterioration in the eastern Mediterranean as revealed by ion microprobe analysis of a speleothem that grew from 2.2 to 0.9 ka in Soreq Cave, Israel. Quat. Res. 71, 27–35 (2009). 9. Orland, I. J. et al. Seasonal resolution of Eastern Mediterranean climate change since 34 ka from a Soreq Cave speleothem. Geochim. Cosmochim. Acta. 89, 240–255 (2012). 10. Smith, C. L. et al. Chronology building using objective identification of annual signals in trace element profiles of stalagmites. Quatern. Geochronol. 4, 11–21 (2009). 11. Nagra, G. et al. Dating stalagmites in mediterranean climates using annual trace element cycles. Sci. Rep. 7, 1–12, https://doi. org/10.1038/s41598-017-00474-4 (2017). 12. Huang, Y. et al. Seasonal variations in Sr, Mg and P in modern speleothems. Chem. Geol. 175, 429–448 (2001). 13. Jochum, K. P. et al. Accurate trace element analysis of speleothems and biogenic calcium carbonates by LA-ICP-MS. Chem. Geol. 318–319, 31–44 (2012). 14. Finch, A. A., Shaw, P. A., Weedon, G. P. & Holmgren K. Trace element variation in speleothem aragonite: potential for palaeoenvironmental reconstruction. Earth Planet. Sci. Lett. 186, 255–267 (2001). 15. Frisia, S. et al. A re-evaluation of the palaeoclimatic significance of phosphorus variability in speleothems revealed by high- resolution synchrotron micro XRF mapping. Clim. Past 8, 2039–2051, https://doi.org/10.5194/cp-8-2039-2012 (2012). 16. Wu, J. Y. et al. Stable isotope and trace element investigation of two contemporaneous annually-laminated stalagmites from northeastern China surrounding the “8.2 ka event”. Clim. Past 8, 1497–1507, https://doi.org/10.5194/cp-8-1497-2012 (2012). 17. Frisia, S., Borsato, A., Fairchild, I. J. & McDermott, F. Calcite fabrics, growth mechanisms, and environments of formation in speleothems from the Italian Alps and South-western Ireland. J. Sediment. Res. 70, 1183–1196 (2000). 18. Spӧtl, C., Fairchild, I. J. & Tooth, A. F. Cave air control on dripwater geochemistry, Obir caves (Austria): implications for speleothem deposition in dynamically ventilated caves. Geochim. Cosmochim. Acta. 69, 2451–2468 (2005). 19. Cai, B., Zhu, J., Ban, F. & Tan, M. Intra-annual variation of the calcite deposition rate of drip water in Shihua Cave, Beijing, China and its implications for palaeoclimatic reconstructions. Boreas. 40(3), 525–535, https://doi.org/10.1111/j.1502-3885.2010.00201.x (2011). 20. Baker, A., Smart, P. L., Edwards, R. L. & Richards, D. Annual growth banding in a cave stalagmite. Nature. 364, 518–520 (1993). 21. Tan, M. et al. Cyclic rapid warming on centennial-scale revealed by a 2650-year stalagmite record of warm season temperature. Geophysical Research Letters. 30(12), 1617–1620, doi:10D1029/2003GL017352D2003 (2003). 22. Ma, Z. et al. Timing and structure of the Younger Dryas event in northern China. Quat. Sci. Rev. 41, 83–93 (2012). 23. Liu, Y. et al. Links between the East Asian monsoon and North Atlantic climate during the 8,200 year event. Nature Geosci. 6, 117–120 (2013). 24. Duan, W. et al. Onset and duration of transitions into Greenland Interstadials 15.2 and 14 in northern China constrained by an annually laminated stalagmite. Sci. Rep. 6, 20844, https://doi.org/10.1038/srep20844 (2016). 25. Li, X. et al. The East Asian summer monsoon variability over the last 145 years inferred from the Shihua Cave record, North China. Sci. Rep. 7, 7078, https://doi.org/10.1038/s41598-017-07251-3 (2017). 26. Fairchild, I. J. et al. Annual to sub-annual resolution of multiple trace-element trends in speleothems. J. Geol. Soc. 158, 831–841 (2001). 27. Fairchild, I. J., Frisia, S., Borsato, A. & Tooth, A. F. Speleothems. In Geochemical Sediments and Landscapes.(eds Nash, D. J., McLaren, S. J.) 200–245 (Blackwells, Oxford, 2007). 28. Tan, M. et al. Bioptical microcycles of laminated speleothems from China and their chronological significance. Chinese Science Bulletin. 44(17), 1604–1607 (1999). 29. Hou, J. Z., Tan, M., Cheng, H. & Liu, T. S. Stable isotope records of plant cover change and monsoon variation in the past 2200 years: evidence from laminated stalagmites in Beijing, China. Boreas. 32, 304–313 (2003). 30. Ban, F., Pan, G., Zhu, J., Cai, B. & Tan, M. Temporal and spatial variations in the discharge and dissolved organic carbon of drip waters in Beijing Shihua Cave, China. Hydrol. Process. 22, 3749–3758 (2008). 31. Liu, H. et al. Structure and evolution of the Shihua Cave in Beijing. Carsologica Sinica (in Chinese). 34(1), 27–34 (2015). 32. Covington, M. D. & Perne., M. Consider a cylindri-cal cave: A physicist’s view of cave and karst science. Acta Carsologic. 44/3, 363–380 (2015). 33. Sanderson, M. J., Smith, I., Parker, I. & Bootman, M. D. FluorescenceMicroscopy. Cold Spring Harb Protoc. 10, 1–36, https://doi. org/10.1101/pdb.top071795 (2014). 34. Fairchild, I. J. & Treble, P. C. Trace elements in speleothems as recorders of environmental change. Quat. Sci. Rev. 28, 449–468 (2009). 35. Partin, J. W., Cobb, K. M., Adkins, J. F., Tueu, A. A. & Clark, B. Trace metal and carbon isotopic variations in cave dripwater and stalagmite geochemistry from northern Borneo. Geochem. Geophy. Geosy. 14(9), 3567–3585, https://doi.org/10.1002/ggge.20215 (2013). 36. Borsato, A., Frisia, S., Fairchild, I. J., Somogyi, A. & Susini, J. Trace element distribution in annual stalagmite laminae mapped by micrometer-resolution X-ray fluorescence: implications for incorporation of environmentally significant species. Geochim. Cosmochim. Acta. 71, 1494–1512 (2007). 37. Rutlidge, H. et al. Dripwater organic matter and trace element geochemistry in a semi-arid karst environment: Implications for speleothem paleoclimatology. Geochim. Cosmochim. Acta. 135, 217–230 (2014). 38. Regattieri, E. et al. Lateglacial to Holocene trace element record (Ba, Mg, Sr.) from Corchia Cave (Apuan Alps, central Italy): paleoenvironmental implications. J. Quat. Sci. 29(4), 381–392, https://doi.org/10.1002/jqs.2712 (2014). Acknowledgements Many thanks to Binggui Cai for his permission for using map of Shihua Cave, to Haiying Sun, Haijian Chen and Linhua Li at the administrative committee of Shihua Cave for field assistance, and to our collaboration, Jie Guo for her helpful discussions. This work was supported by National Natural Science Foundation of China (NSFC) (Grant No. 41402161) Natural Science Foundation of Shanxi Province, China (Grant No. 2015011085) and the Mark Wainwright Analytical Centre at UNSW, Australia. Author Contributions F.M.B., A.B. and C.E.M. conceived the project and wrote the first draft of the paper. W.H.D. contributed to the manuscript revision at different stages. F.M.B., K.C. and R.A. contribute to trace elements measurements. F.M.B., X.L.L. and J.X.H. did the field and monitoring work. F.M.B. and X.L.L. collected the XMG stalagmite sample. G.N. contributes to guide suggestion of Peak Analyzer software. M.T. provided many suggestions of counting layers. Data collected in this study has been archived at the Mark Wainwright Analytical Centre at UNSW, Australia and is available on request from the authors. Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 8 www.nature.com/scientificreports/ Additional Information Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-018-22839-z. Competing Interests: The authors declare no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2018 Scientific REpo R TS | (2018) 8:4551 | DOI:10.1038/s41598-018-22839-z 9

Journal

Scientific ReportsSpringer Journals

Published: Mar 14, 2018

There are no references for this article.