# Chinese stalagmite δ<SUP>18</SUP>/<SUP>18</SUP>O records reveal the diverse moisture trajectories during the middle to late last glacial period

Huihui Yang, Yu-Min Chou, Xiuyang Jiang, Wei Zheng, Yaoqi He, Yogaraj Banerjee, Chuan-Chou Shen, Tsai-Luen Yu, Yi Zhong, F Humbert, et al.

► To cite this version:

Huihui Yang, Yu-Min Chou, Xiuyang Jiang, Wei Zheng, Yaoqi He, et al. Chinese stalagmite δ<SUP>18</SUP>/<SUP>18</SUP>O records reveal the diverse moisture trajectories during the middle to late last glacial period. Geological Magazine, 2024, 10.1017/S0016756824000013. hal-04577754

HAL Id: hal-04577754

https://hal.science/hal-04577754

Submitted on 16 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
## Original Article

**Cite this article:** Yang H, Chou Y-M, Jiang X, Zheng W, He Y, Banerjee Y, Shen C-C, Yu T-L, Zhong Y, Humbert F, and Liu Q. Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period. *Geological Magazine* https://doi.org/10.1017/S0016756824000013

Received: 30 June 2023
Revised: 4 December 2023
Accepted: 3 January 2024

**Keywords:**

stalagmite δ¹⁸O; Heinrich stadials; moisture trajectory; East Asian Summer Monsoon; Indian Summer Monsoon

**Corresponding authors:**

Yu-Min Chou; Email: chouym@sustech.edu.cn; Xiuyang Jiang; Email: xyjiang@fjnu.edu.cn

# Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period

Huihui Yang¹,² ID, Yu-Min Chou², Xiuyang Jiang³, Wei Zheng², Yaoqi He⁴, Yogaraj Banerjee⁵,⁶, Chuan-Chou Shen⁵,⁶, Tsai-Luen Yu⁷, Yi Zhong², Fabien Humbert²,⁸ ID and Qingsong Liu²

¹School of Environment, Harbin Institute of Technology, Harbin, 150006, China; ²Center for Marine Magnetism (CM²), Department of Ocean Science and Engineering, Southern University of Science and Technology, Shenzhen, 518055, China; ³Key Laboratory of Humid Subtropical Eco-Geographical Processes, College of Geography Science, Fujian Normal University, Fuzhou, 350117, China; ⁴College of Tourism and Air Service, Guizhou Minzu University, Guiyang, 550025, China; ⁵High-Precision Mass Spectrometry and Environment Change Laboratory (HISPEC), Department of Geosciences, National Taiwan University, Taipei, 10617, Taiwan ROC; ⁶Marine Industry and Engineering Research Center, National Academy of Marine Research, Taipei, 10617, Taiwan ROC; ⁷Research Center for Future Earth, National Taiwan University, Taipei 10617, Taiwan ROC and ⁸University of Rennes, CNRS, Géosciences Rennes, UMR 6118, Rennes, France

## Abstract

Based on 30 high-resolution U-Th dating controls, we reconstruct stalagmite δ¹⁸O records from 45 to 15 thousand years ago (ka B.P., before AD 1950) from the Shizhu Cave, which is located in southwestern China under the influence of both the Indian Summer Monsoon (ISM) and the East Asian Summer Monsoon (EASM). By integrating with the other stalagmite δ¹⁸O records in Asia during the middle to late last glacial, our results reveal two main moisture trajectories: one from the Indian Ocean, through the Shizhu Cave towards central China, and the other from the Pacific Ocean to central and northern China. The systematic decrease of the average values of stalagmite δ¹⁸O records from oceans to inland China reveals a spatial pattern of water vapour fractionation and moisture trajectory during the middle to late last glacial. In contrast, the variation amplitude, which is defined as the departures apart from the background δ¹⁸O records during Heinrich stadials 1 to 4 (HS1-HS4), show an increasing trend from the coastal oceans to mid-latitude inland China, presenting a 'coastal-inland' pattern, which can be interpreted by the enhanced East Asian Winter Monsoon (EAWM) and the weakened EASM. More specifically, the enriched stalagmite δ¹⁸O records in the EASM region during HS1 to HS4 are caused by the decreased summer rainfall amount or/and the increased proportion of summer moisture resources from the Pacific Ocean. These new observations deepen our understanding of the complicated stalagmite δ¹⁸O records in the EASM region.

## 1. Introduction

Stalagmite δ¹⁸O has been one of the most important proxies to track climate evolution during the late Quaternary in recent decades (e.g., Cheng et al. 2016, 2019). However, the interpretation of stalagmite δ¹⁸O records from regions impacted by the Asian summer monsoon remains controversial. Wang et al. (2001) suggested that the evolution of δ¹⁸O records in the east China stalagmite is driven predominantly by ratios of the summer and winter precipitation. Alternatively, Yuan et al. (2004) argued that variations in stalagmite δ¹⁸O in the Dongge Cave were controlled mostly by the fraction of water vapour between the moisture sources and cave sites. Other views also suggested the contribution of various degrees of different moisture sources from the Indian and Pacific Oceans on the δ¹⁸O signature (Maher & Thompson, 2012; Baker et al. 2015). Some studies noted that reconciling those factors may better explain the δ¹⁸O of precipitation and stalagmites (Liu et al. 2014).

Pausata et al. (2011) suggested that Chinese stalagmite δ¹⁸O was controlled by fluctuations in the Indian Summer Monsoon (ISM) intensity using climate model simulations for Heinrich stadial 1. They indicated that decreased intensity of the ISM results in enriched stalagmite δ¹⁸O in eastern China during Heinrich stadial 1. Several geological archives have also recorded the contribution of the Indian Ocean to the δ¹⁸O~P~ (δ¹⁸O of precipitation) in eastern China (Yang et al. 2014; Chen et al. 2015; Liu et al. 2015). However, this viewpoint was challenged by other geological records (Tan, 2009; Tan, 2013; Li et al. 2019; Liang et al. 2020). For example, Tan et al. argued that moisture sources of the monsoon regions of China including both the Indian Ocean and the Pacific Ocean, with depleted and enriched δ¹⁸O signal in precipitation (Tan, 2009, 2013). They found that the δ¹⁸O signal in stalagmite reflects the intensity contrast between two

© The Author(s), 2024. Published by Cambridge University Press. This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.

<figure><img src="image_1.png" /></figure>

CAMBRIDGE
UNIVERSITY PRESS
## Original Article

**Cite this article:** Yang H, Chou Y-M, Jiang X, Zheng W, He Y, Banerjee Y, Shen C-C, Yu T-L, Zhong Y, Humbert F, and Liu Q. Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period. *Geological Magazine* https://doi.org/10.1017/S0016756824000013

Received: 30 June 2023
Revised: 4 December 2023
Accepted: 3 January 2024

**Keywords:**
stalagmite δ¹⁸O; Heinrich stadials; moisture trajectory; East Asian Summer Monsoon; Indian Summer Monsoon

**Corresponding authors:**
Yu-Min Chou; Email: chouym@sustech.edu.cn; Xiuyang Jiang; Email: xyjiang@fjnu.edu.cn

# Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period

Huihui Yang<fnref n="1" />,² ID, Yu-Min Chou², Xiuyang Jiang³, Wei Zheng², Yaoqi He⁴, Yogaraj Banerjee⁵,⁶, Chuan-Chou Shen⁵,⁶, Tsai-Luen Yu⁷, Yi Zhong², Fabien Humbert²,⁸ ID and Qingsong Liu²

## Abstract

Based on 30 high-resolution U-Th dating controls, we reconstruct stalagmite δ¹⁸O records from 45 to 15 thousand years ago (ka B.P., before AD 1950) from the Shizhu Cave, which is located in southwestern China under the influence of both the Indian Summer Monsoon (ISM) and the East Asian Summer Monsoon (EASM). By integrating with the other stalagmite δ¹⁸O records in Asia during the middle to late last glacial, our results reveal two main moisture trajectories: one from the Indian Ocean, through the Shizhu Cave towards central China, and the other from the Pacific Ocean to central and northern China. The systematic decrease of the average values of stalagmite δ¹⁸O records from oceans to inland China reveals a spatial pattern of water vapour fractionation and moisture trajectory during the middle to late last glacial. In contrast, the variation amplitude, which is defined as the departures apart from the background δ¹⁸O records during Heinrich stadials 1 to 4 (HS1-HS4), show an increasing trend from the coastal oceans to mid-latitude inland China, presenting a ‘coastal-inland’ pattern, which can be interpreted by the enhanced East Asian Winter Monsoon (EAWM) and the weakened EASM. More specifically, the enriched stalagmite δ¹⁸O records in the EASM region during HS1 to HS4 are caused by the decreased summer rainfall amount or/and the increased proportion of summer moisture resources from the Pacific Ocean. These new observations deepen our understanding of the complicated stalagmite δ¹⁸O records in the EASM region.

## 1. Introduction

Stalagmite δ¹⁸O has been one of the most important proxies to track climate evolution during the late Quaternary in recent decades (e.g., Cheng et al. 2016, 2019). However, the interpretation of stalagmite δ¹⁸O records from regions impacted by the Asian summer monsoon remains controversial. Wang et al. (2001) suggested that the evolution of δ¹⁸O records in the east China stalagmite is driven predominantly by ratios of the summer and winter precipitation. Alternatively, Yuan et al. (2004) argued that variations in stalagmite δ¹⁸O in the Dongge Cave were controlled mostly by the fraction of water vapour between the moisture sources and cave sites. Other views also suggested the contribution of various degrees of different moisture sources from the Indian and Pacific Oceans on the δ¹⁸O signature (Maher & Thompson, 2012; Baker et al. 2015). Some studies noted that reconciling those factors may better explain the δ¹⁸O of precipitation and stalagmites (Liu et al. 2014).

Pausata et al. (2011) suggested that Chinese stalagmite δ¹⁸O was controlled by fluctuations in the Indian Summer Monsoon (ISM) intensity using climate model simulations for Heinrich stadial 1. They indicated that decreased intensity of the ISM results in enriched stalagmite δ¹⁸O in eastern China during Heinrich stadial 1. Several geological archives have also recorded the contribution of the Indian Ocean to the δ¹⁸O~P~ (δ¹⁸O of precipitation) in eastern China (Yang et al. 2014; Chen et al. 2015; Liu et al. 2015). However, this viewpoint was challenged by other geological records (Tan, 2009; Tan, 2013; Li et al. 2019; Liang et al. 2020). For example, Tan et al. argued that moisture sources of the monsoon regions of China including both the Indian Ocean and the Pacific Ocean, with depleted and enriched δ¹⁸O signal in precipitation (Tan, 2009, 2013). They found that the δ¹⁸O signal in stalagmite reflects the intensity contrast between two

<footnote n="1">¹School of Environment, Harbin Institute of Technology, Harbin, 150006, China; ²Center for Marine Magnetism (CM²), Department of Ocean Science and Engineering, Southern University of Science and Technology, Shenzhen, 518055, China; ³Key Laboratory of Humid Subtropical Eco-Geographical Processes, College of Geography Science, Fujian Normal University, Fuzhou, 350117, China; ⁴College of Tourism and Air Service, Guizhou Minzu University, Guiyang, 550025, China; ⁵High-Precision Mass Spectrometry and Environment Change Laboratory (HISPEC), Department of Geosciences, National Taiwan University, Taipei, 10617, Taiwan ROC; ⁶Marine Industry and Engineering Research Center, National Academy of Marine Research, Taipei, 10617, Taiwan ROC; ⁷Research Center for Future Earth, National Taiwan University, Taipei 10617, Taiwan ROC and ⁸University of Rennes, CNRS, Géosciences Rennes, UMR 6118, Rennes, France</footnote>
<figure><img src="image_2.png" /><figcaption>**Figure 1.** (Colour online) Cave locations in Asia. The detailed information for every cave is as follows: Mawmluh Cave: 25°15′N, 91°52′E, 1290 m a.s.l. (Dutt *et al.* 2015); Shizhu Cave: 26°02′N, 107°16′E, 950 m a.s.l.; Xianyun Cave: 116°59′E, 25°33′N; 970 m a.s.l. (Zhang *et al.* 2021; Qiu *et al.* 2022); Hulu Cave: 32°32′N, 119°10′E, 90 m a.s.l. (Wang *et al.* 2001); Yongxing Cave: 31°35′N, 111°14′E, 750 m a.s.l. (Chen *et al.* 2016); Songjia Cave: 32°24′N, 107°10′E; 680 m a.s.l. (Zhou *et al.* 2008). The little black arrows represent the present-day averaged June–July–August (JJA) wind at 850 hPa (NCEP/NCAR, 1994). Four arrows represent the ISM, EASM, Westerly and Winter monsoon, respectively.</figcaption></figure>

summer monsoons, the ISM and the East Asian Summer Monsoon (EASM), which is further dependent on the variation of the Western Pacific Subtropical High (WPSH). Li *et al.* (2019) subtracted the typical $\delta^{18}$O signal of the ISM from a stalagmite record of eastern China during the Heinrich stadial 1 to isolate the Pacific’s $\delta^{18}$O signal. The latter signal shows a clear negative correlation with the Sea Surface Temperature (SST) of the Western Pacific Warm Pool.

Moreover, the position of the westerly jet in summer is also important in affecting stalagmite $\delta^{18}$O records in China. The shifts of the westerly jet change the geographic position and $\delta^{18}$O$_P$ of the monsoon rainband (Sampe & Xie, 2010; Baker *et al.* 2015; Chiang *et al.* 2015, 2017, 2020). The rapid millennial-scale climatic changes during the last glacial occurred in the high-latitude Northern Hemisphere and transduced to East Asia, which is realized through the northern westerlies (e.g., Porter & An, 1995; Guo *et al.* 1996; Nagashima *et al.* 2011; Sun *et al.* 2012).

We suggest that systematic comparisons of the stalagmite $\delta^{18}$O records at different spatial-temporal scales can help to resolve these disputes. The last glacial is the most recent period when global climate changed dramatically, with expanding ice sheet, lowering sea levels, exposing shelves (Siddall *et al.* 2003, 2008), dropping in SST (Kucera *et al.* 2005) and reorganizing in atmospheric-oceanic circulation (Bohm *et al.* 2015; Markle *et al.* 2016; Lynch-Stieglitz, 2017; Dong *et al.* 2022). In the Asian region, atmospheric-oceanic circulation changed, with a southward shift of the Intertropical Convergence Zone (ITCZ) in the South China Sea (Zhao *et al.* 2018; Kaboth-Bahr *et al.* 2021), weakened Asian summer monsoon (Cheng *et al.* 2012, 2016) and intensified winter monsoon (Sun *et al.* 2006, 2012).

Here, we established a high-resolution $\delta^{18}$O record from two stalagmites of the Shizhu Cave from 45 to 15 ka B.P. in

southwestern China (Fig. 1). The hydroclimate there is dominated by the ISM but also partly influenced by the EASM and the westerlies (Wang & Chen, 2012; Cai *et al.* 2015). Thus, stalagmites of the Shizhu Cave are suitable materials to better constrain the moisture trajectory of ISM in the past. Further, we compared $\delta^{18}$O records of the Shizhu Cave with other available records that were widely deposited in the ISM and EASM regions during the middle to late last glacial to investigate the possible moisture trajectories of the Asian summer monsoon during this period. In this study, first, we explore the possible moisture trajectories over the ISM and EASM regions during the middle to late last glacial. Second, we obtain amplitudes of stalagmite $\delta^{18}$O records from caves in the Asian region during Heinrich stadial 1 to 4 (HS1 to HS4) and investigate their distribution pattern to further discuss the controlling factors behind the pattern.

## 2. Materials and methods

### 2.a. Setting and sampling

The Shizhu Cave is located in the southern Guizhou Province (26°02′N, 107°28′E; 1089 m above sea level), southwestern China (Fig. 1). The mean annual temperature and precipitation are 16.1 °C and 1431.1 mm, respectively (http://data.cma.cn/dataService). Warm and humid air from the Indian Ocean and the equatorial Pacific Ocean generates heavy and continuous rain in this region from May to October, constituting about 80% of the yearly precipitation (Dykoski *et al.* 2005). Two stalagmites, SZ-1 and SZ-3 (Fig. 2), were sampled from the Shizhu Cave at a chamber approximately 50 m below the surface. SZ-1 and SZ-3 are 430 mm and 542.5 mm in length, respectively (Fig. 3).
<figure><img src="image_2.png" /><figcaption>**Figure 1.** (Colour online) Cave locations in Asia. The detailed information for every cave is as follows: Mawmluh Cave: 25°15′N, 91°52′E, 1290 m a.s.l. (Dutt *et al.* 2015); Shizhu Cave: 26°02′N, 107°16′E, 950 m a.s.l.; Xianyun Cave: 116°59′E, 25°33′N; 970 m a.s.l. (Zhang *et al.* 2021; Qiu *et al.* 2022); Hulu Cave: 32°32′N, 119°10′E, 90 m a.s.l. (Wang *et al.* 2001); Yongxing Cave: 31°35′N, 111°14′E, 750 m a.s.l. (Chen *et al.* 2016); Songjia Cave: 32°24′N, 107°10′E; 680 m a.s.l. (Zhou *et al.* 2008). The little black arrows represent the present-day averaged June–July–August (JJA) wind at 850 hPa (NCEP/NCAR, 1994). Four arrows represent the ISM, EASM, Westerly and Winter monsoon, respectively.</figcaption></figure>

summer monsoons, the ISM and the East Asian Summer Monsoon (EASM), which is further dependent on the variation of the Western Pacific Subtropical High (WPSH). Li *et al.* (2019) subtracted the typical $\delta^{18}$O signal of the ISM from a stalagmite record of eastern China during the Heinrich stadial 1 to isolate the Pacific's $\delta^{18}$O signal. The latter signal shows a clear negative correlation with the Sea Surface Temperature (SST) of the Western Pacific Warm Pool.

Moreover, the position of the westerly jet in summer is also important in affecting stalagmite $\delta^{18}$O records in China. The shifts of the westerly jet change the geographic position and $\delta^{18}$O$_P$ of the monsoon rainband (Sampe & Xie, 2010; Baker *et al.* 2015; Chiang *et al.* 2015, 2017, 2020). The rapid millennial-scale climatic changes during the last glacial occurred in the high-latitude Northern Hemisphere and transduced to East Asia, which is realized through the northern westerlies (e.g., Porter & An, 1995; Guo *et al.* 1996; Nagashima *et al.* 2011; Sun *et al.* 2012).

We suggest that systematic comparisons of the stalagmite $\delta^{18}$O records at different spatial-temporal scales can help to resolve these disputes. The last glacial is the most recent period when global climate changed dramatically, with expanding ice sheet, lowering sea levels, exposing shelves (Siddall *et al.* 2003, 2008), dropping in SST (Kucera *et al.* 2005) and reorganizing in atmospheric-oceanic circulation (Bohm *et al.* 2015; Markle *et al.* 2016; Lynch-Stieglitz, 2017; Dong *et al.* 2022). In the Asian region, atmospheric-oceanic circulation changed, with a southward shift of the Intertropical Convergence Zone (ITCZ) in the South China Sea (Zhao *et al.* 2018; Kaboth-Bahr *et al.* 2021), weakened Asian summer monsoon (Cheng *et al.* 2012, 2016) and intensified winter monsoon (Sun *et al.* 2006, 2012).

Here, we established a high-resolution $\delta^{18}$O record from two stalagmites of the Shizhu Cave from 45 to 15 ka B.P. in

southwestern China (Fig. 1). The hydroclimate there is dominated by the ISM but also partly influenced by the EASM and the westerlies (Wang & Chen, 2012; Cai *et al.* 2015). Thus, stalagmites of the Shizhu Cave are suitable materials to better constrain the moisture trajectory of ISM in the past. Further, we compared $\delta^{18}$O records of the Shizhu Cave with other available records that were widely deposited in the ISM and EASM regions during the middle to late last glacial to investigate the possible moisture trajectories of the Asian summer monsoon during this period. In this study, first, we explore the possible moisture trajectories over the ISM and EASM regions during the middle to late last glacial. Second, we obtain amplitudes of stalagmite $\delta^{18}$O records from caves in the Asian region during Heinrich stadial 1 to 4 (HS1 to HS4) and investigate their distribution pattern to further discuss the controlling factors behind the pattern.

## 2. Materials and methods

### 2.a. Setting and sampling

The Shizhu Cave is located in the southern Guizhou Province (26°02′N, 107°28′E; 1089 m above sea level), southwestern China (Fig. 1). The mean annual temperature and precipitation are 16.1 °C and 1431.1 mm, respectively (http://data.cma.cn/dataService). Warm and humid air from the Indian Ocean and the equatorial Pacific Ocean generates heavy and continuous rain in this region from May to October, constituting about 80% of the yearly precipitation (Dykoski *et al.* 2005). Two stalagmites, SZ-1 and SZ-3 (Fig. 2), were sampled from the Shizhu Cave at a chamber approximately 50 m below the surface. SZ-1 and SZ-3 are 430 mm and 542.5 mm in length, respectively (Fig. 3).
<figure><img src="image_3.png" /><figcaption>**Figure 2.** (Colour online) Age models and growth rates of the stalagmites SZ-1 and SZ-3 from Shizhu Cave. Age–depth relationship for SZ-1 (green line and axis) and SZ-3 (blue line and axis). Chronologies are established by U-Th dating and the StalAge model; error bars are indicated for every dating in red.</figcaption></figure>

## 2.b. U-Th dating

The U-Th dating method is one of the most important advantages of stalagmite paleoclimate research because of its high degree of accuracy in absolute timing (Cheng et al. 2019). Thirty subsamples with 20 mm in length, 10 mm in width and 2.5 mm in depth were collected discontinuously from SZ-1 and SZ-3 along their growth axis for U-Th dating (Table 1). Uranium and thorium isotopic composition and content were measured using a multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) in the High-Precision Mass Spectrometry and Environment Change Laboratory (HISPEC), Department of Geoscience, National Taiwan University (Shen et al. 2012). Uncertainties of the isotope data and U-Th dates, relative to 1950 AD, are given at the two-sigma (2s) level or two standard deviations of the mean (2s_m) unless otherwise noted. Decay constants used for age calculation are available in Cheng et al. (2013).

## 2.c. Oxygen isotopic measurements

The measurement technology of oxygen isotopes has been already quite mature after years of development. Here, 406 and 276 subsamples (cylinder of 0.3 mm diameter wide) were drilled,

respectively, from the SZ-1 and SZ-3 with carbide dental burrs at a spacing of 2 mm along their growth axes. Acidic hydrolysis of CaCO₃ was then performed from the extracted material, to which 0.2 mg of calcite was added to each subsample, using H₃PO₄ at a concentration close to 100 % at a constant controlled temperature of ~72 °C to extract the CO₂. Released CO₂ was carried with helium (He) flow to an accessory Finnigan Gas Bench, where water vapour was separated through gas chromatography. The resultant CO₂ was ionized and accelerated to mass separation in an automated system. Measurements of oxygen and carbon stable isotopes were performed on a Gasbench attached to a Thermo MAT 253 mass spectrometer at Fujian Normal University with an external 1σ precision of ± 0.06‰. Oxygen isotope values were presented as δ¹⁸O (δ¹⁸O = [(¹⁸O/¹⁶O)sample/($^{18}$O/$^{16}$O)standard - 1] × 1000‰) with the Vienna Pee Dee Belemnite (VPDB) standard.

## 2.d. X-ray diffraction (XRD) analysis

Calcite is more stable than other carbonate materials like aragonite in recording stable isotope information (Zhang et al. 2014). Thus, it is necessary to detect the properties of the research materials. Eight subsamples from the sectioned surfaces of the stalagmite SZ-1 and SZ-3 were collected for XRD analysis. For each subsample, 2 mg powder was loaded on stages with a silicon tablet background and analysed by a Rigaku MiniFlex benchtop X-ray diffraction (XRD) instrument (Cu-Kα radiation) in the Center for Social Science, Southern University of Science and Technology. The XRD spectra were obtained from 2° to 100° 2θ. The major phase identification and semi-quantification were conducted using the PDXL2 software.

## 2.e. Data collecting and statistical analysis

The stack curve between stalagmite δ¹⁸O records of SZ-1 and SZ-3 in Figure 4 and the stack curve among five caves in Figure 7b were both obtained using the bootstrap resampling method, which can reserve the general trend of different curves as much as possible. Before stacking SZ-1 and SZ-3, we calculated the Pearson correlation coefficient between them in Fig. S1. The stack curve of SZ-1 and SZ-3 and the error bars are shown in Figure 4.

# 3. Results

## 3.a. Geochronology

The measured isotope ratios of U and Th, decay constants and calculated ages are listed in Table 1. Most of the 2σ uncertainties are better than 200 years for SZ-1 and 100 years for SZ-3, respectively. The StalAge chronology was used to establish the age models. Ages versus depths show that stalagmites SZ-1 and SZ-3 continuously grew from 45.0 to 22.0 ka B.P. and 32.8 to 14.4 ka B.P., respectively (Figs. 2 and 3). The growth rates vary from 1 to 7 cm/ka B.P. for SZ-1 and from 0.4 to 2.4 cm/ka B.P. for SZ-3. It is noted that samples before 22.28 ka B.P. of SZ-1 were lost before obtaining U-Th results, corresponding to a depth of about 0–9.75 cm. Thus, the δ¹⁸O records before 22.28 ka are unusable.

## 3.b. XRD results

The XRD results (Fig. S2) indicate that samples from SZ-1 and SZ-3 are predominantly composed of calcite, which makes them suitable for studying oxygen isotope variations. It is noted that ‘unclear laminations’ occurring during the period 22.28 to 24.64 ka B.P. were observed in stalagmite SZ-1 (Fig. 3); they may be formed by the recrystallization of the primary calcite crystals, alter the
<figure><img src="image_3.png" /><figcaption>**Figure 2.** (Colour online) Age models and growth rates of the stalagmites SZ-1 and SZ-3 from Shizhu Cave. Age–depth relationship for SZ-1 (green line and axis) and SZ-3 (blue line and axis). Chronologies are established by U-Th dating and the StalAge model; error bars are indicated for every dating in red.</figcaption></figure>

## 2.b. U-Th dating

The U-Th dating method is one of the most important advantages of stalagmite paleoclimate research because of its high degree of accuracy in absolute timing (Cheng *et al.* 2019). Thirty subsamples with 20 mm in length, 10 mm in width and 2.5 mm in depth were collected discontinuously from SZ-1 and SZ-3 along their growth axis for U-Th dating (Table 1). Uranium and thorium isotopic composition and content were measured using a multi-collector inductively coupled plasma mass spectrometer (MC-ICP-MS) in the High-Precision Mass Spectrometry and Environment Change Laboratory (HISPEC), Department of Geoscience, National Taiwan University (Shen *et al.* 2012). Uncertainties of the isotope data and U-Th dates, relative to 1950 AD, are given at the two-sigma (2s) level or two standard deviations of the mean (2s_m) unless otherwise noted. Decay constants used for age calculation are available in Cheng *et al.* (2013).

## 2.c. Oxygen isotopic measurements

The measurement technology of oxygen isotopes has been already quite mature after years of development. Here, 406 and 276 subsamples (cylinder of 0.3 mm diameter wide) were drilled,

respectively, from the SZ-1 and SZ-3 with carbide dental burrs at a spacing of 2 mm along their growth axes. Acidic hydrolysis of CaCO₃ was then performed from the extracted material, to which 0.2 mg of calcite was added to each subsample, using H₃PO₄ at a concentration close to 100 % at a constant controlled temperature of ~72 °C to extract the CO₂. Released CO₂ was carried with helium (He) flow to an accessory Finnigan Gas Bench, where water vapour was separated through gas chromatography. The resultant CO₂ was ionized and accelerated to mass separation in an automated system. Measurements of oxygen and carbon stable isotopes were performed on a Gasbench attached to a Thermo MAT 253 mass spectrometer at Fujian Normal University with an external 1σ precision of ± 0.06‰. Oxygen isotope values were presented as δ¹⁸O (δ¹⁸O = [(¹⁸O/¹⁶O)sample/($^{18}$O/$^{16}$O)standard - 1] × 1000‰) with the Vienna Pee Dee Belemnite (VPDB) standard.

## 2.d. X-ray diffraction (XRD) analysis

Calcite is more stable than other carbonate materials like aragonite in recording stable isotope information (Zhang *et al.* 2014). Thus, it is necessary to detect the properties of the research materials. Eight subsamples from the sectioned surfaces of the stalagmite SZ-1 and SZ-3 were collected for XRD analysis. For each subsample, 2 mg powder was loaded on stages with a silicon tablet background and analysed by a Rigaku MiniFlex benchtop X-ray diffraction (XRD) instrument (Cu-Kα radiation) in the Center for Social Science, Southern University of Science and Technology. The XRD spectra were obtained from 2° to 100° 2θ. The major phase identification and semi-quantification were conducted using the PDXL2 software.

## 2.e. Data collecting and statistical analysis

The stack curve between stalagmite δ¹⁸O records of SZ-1 and SZ-3 in Figure 4 and the stack curve among five caves in Figure 7b were both obtained using the bootstrap resampling method, which can reserve the general trend of different curves as much as possible. Before stacking SZ-1 and SZ-3, we calculated the Pearson correlation coefficient between them in Fig. S1. The stack curve of SZ-1 and SZ-3 and the error bars are shown in Figure 4.

# 3. Results

## 3.a. Geochronology

The measured isotope ratios of U and Th, decay constants and calculated ages are listed in Table 1. Most of the 2σ uncertainties are better than 200 years for SZ-1 and 100 years for SZ-3, respectively. The StalAge chronology was used to establish the age models. Ages versus depths show that stalagmites SZ-1 and SZ-3 continuously grew from 45.0 to 22.0 ka B.P. and 32.8 to 14.4 ka B.P., respectively (Figs. 2 and 3). The growth rates vary from 1 to 7 cm/ka B.P. for SZ-1 and from 0.4 to 2.4 cm/ka B.P. for SZ-3. It is noted that samples before 22.28 ka B.P. of SZ-1 were lost before obtaining U-Th results, corresponding to a depth of about 0–9.75 cm. Thus, the δ¹⁸O records before 22.28 ka are unusable.

## 3.b. XRD results

The XRD results (Fig. S2) indicate that samples from SZ-1 and SZ-3 are predominantly composed of calcite, which makes them suitable for studying oxygen isotope variations. It is noted that ‘unclear laminations’ occurring during the period 22.28 to 24.64 ka B.P. were observed in stalagmite SZ-1 (Fig. 3); they may be formed by the recrystallization of the primary calcite crystals, alter the
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (Colour online) Cross sections of the stalagmites SZ-1 (a) and SZ-3 (b); the scale is the same for both stalagmites. The yellow rectangle in SZ-1 highlights the area where recrystallization might have occurred. Dating results and related dated areas of each stalagmite are indicated by the red rectangles. The white bars are remnants of the original marks and can be ignored. The oxygen isotope sampling track was performed directly on the left of the central axis of each stalagmite, which is marked by the orange lines.</figcaption></figure>

primary stable isotope composition (Zhang *et al.* 2014; Lachniet, 2015) and thus have been excluded for further paleoclimate interpretations (Fig. 4a).

### 3.c. Combination of curve SZ-1-3

The δ¹⁸O values of SZ-1 range between -2.31‰ and -7.87‰ with an average value of -5.50‰, while the δ¹⁸O values of SZ-3 varied from -2.38‰ to -7.47‰ with an average value of -5.60‰. Centennial- to decadal-scale variability of SZ-1 and SZ-3 can be identified (Fig. 4), with averages of temporal resolution of 128 years (SZ-1) and 66 years (SZ-3). Although the evolution of δ¹⁸O between SZ-1 and SZ-3 (Fig. 4c) could not overlap perfectly due to different growth rates and sampling resolutions during 32.65 to 24.64 ka B.P., their overall trend, amplitudes and the average value of variations are comparable. The result of the Pearson correlation coefficient in Fig. S1 is 0.56, indicating a medium correlation and implying that SZ-1 and SZ-3 are available for further stack analysis. We stacked these two curves deposited between 24.64 ka and 32.65 ka to composite SZ-1-3 with the error curves in Figure 4d. The centennial-scale amplitudes of SZ-1-3 are 2.53 ‰, 1.05 ‰, 0.96 ‰ and 2.86 ‰ at nearly 16 ka B.P., 24 ka B.P., 31 ka B.P. and 39 ka B.P., respectively, indicating the existence of HS1 to HS4 (Wang *et al.* 2001; Andersen, 2004). Additionally, the average values of SZ-1-3 during Marine Isotope Stages 2 and 3 (MIS 2 and MIS 3) are -5.16‰ and -5.56‰, respectively.

## 4. Discussion

### 4.a. The Asian stalagmite δ¹⁸O records reveal moisture trajectories during the middle to late glacial

Previous studies have shown that variations in the stalagmite δ¹⁸O indicate the mean state of summer monsoon intensity and/or the

integrated moisture transports (Cheng *et al.* 2012, 2016). However, the exact mechanisms controlling the stalagmite δ¹⁸O values are still widely debated, including the 'upstream depletion mechanism' (Yuan *et al.* 2004), the ratio of summer to winter rainfall (Wang *et al.* 2001) and the ratio of water vapour between the Pacific Ocean and the Indian Ocean (Tan, 2013). We suggest that systematic comparison at different spatial-temporal scales may partially settle those disputes. Then, we selected relatively continuous stalagmites δ¹⁸O records from caves that are far enough to capture spatial variation in the Asian region during the middle to late last glacial for comparison.

The stalagmite δ¹⁸O records in the monsoon region of China are controlled by δ¹⁸O of yearly precipitation (Qiu *et al.* 2021), of which almost 80% is contributed by the EASM and the last 20% is dominated by the East Asian Winter Monsoon (EAWM) in the current climate setting. While, under the background of weakened Asian summer monsoon and intensified winter monsoon during the last glacial (Sun *et al.* 2006), this ratio may change partially. However, in this discussion, we do not explore the changed proportion emphasized by Wang *et al.* (2001) but focus on the averages of these stalagmite δ¹⁸O records over the last glacial period.

Wang *et al.* (2001) explained the changed stalagmite δ¹⁸O records of the Hulu Cave as the ratio between the East Asia summer and winter monsoon, which is very reasonable. Whereas, when comparing averages of stalagmite δ¹⁸O records from different regions over the last glacial period, the enhanced winter winds had a relatively common effect on these caves, although the influences on these caves may differ. It is conceivable that summer precipitation still constitutes the majority of yearly precipitation (Jiang *et al.* 2011; Cai *et al.* 2015) and dominates the composition of stalagmite δ¹⁸O records. Thus, the comparison among averages of stalagmite δ¹⁸O records reflects summer moisture sources and
<figure><img src="image_4.png" /><figcaption>**Figure 3.** (Colour online) Cross sections of the stalagmites SZ-1 (a) and SZ-3 (b); the scale is the same for both stalagmites. The yellow rectangle in SZ-1 highlights the area where recrystallization might have occurred. Dating results and related dated areas of each stalagmite are indicated by the red rectangles. The white bars are remnants of the original marks and can be ignored. The oxygen isotope sampling track was performed directly on the left of the central axis of each stalagmite, which is marked by the orange lines.</figcaption></figure>

primary stable isotope composition (Zhang *et al.* 2014; Lachniet, 2015) and thus have been excluded for further paleoclimate interpretations (Fig. 4a).

### 3.c. Combination of curve SZ-1-3

The δ¹⁸O values of SZ-1 range between -2.31‰ and -7.87‰ with an average value of -5.50‰, while the δ¹⁸O values of SZ-3 varied from -2.38‰ to -7.47‰ with an average value of -5.60‰. Centennial- to decadal-scale variability of SZ-1 and SZ-3 can be identified (Fig. 4), with averages of temporal resolution of 128 years (SZ-1) and 66 years (SZ-3). Although the evolution of δ¹⁸O between SZ-1 and SZ-3 (Fig. 4c) could not overlap perfectly due to different growth rates and sampling resolutions during 32.65 to 24.64 ka B.P., their overall trend, amplitudes and the average value of variations are comparable. The result of the Pearson correlation coefficient in Fig. S1 is 0.56, indicating a medium correlation and implying that SZ-1 and SZ-3 are available for further stack analysis. We stacked these two curves deposited between 24.64 ka and 32.65 ka to composite SZ-1-3 with the error curves in Figure 4d. The centennial-scale amplitudes of SZ-1-3 are 2.53 ‰, 1.05 ‰, 0.96 ‰ and 2.86 ‰ at nearly 16 ka B.P., 24 ka B.P., 31 ka B.P. and 39 ka B.P., respectively, indicating the existence of HS1 to HS4 (Wang *et al.* 2001; Andersen, 2004). Additionally, the average values of SZ-1-3 during Marine Isotope Stages 2 and 3 (MIS 2 and MIS 3) are -5.16‰ and -5.56‰, respectively.

## 4. Discussion

### 4.a. The Asian stalagmite δ¹⁸O records reveal moisture trajectories during the middle to late glacial

Previous studies have shown that variations in the stalagmite δ¹⁸O indicate the mean state of summer monsoon intensity and/or the

integrated moisture transports (Cheng *et al.* 2012, 2016). However, the exact mechanisms controlling the stalagmite δ¹⁸O values are still widely debated, including the 'upstream depletion mechanism' (Yuan *et al.* 2004), the ratio of summer to winter rainfall (Wang *et al.* 2001) and the ratio of water vapour between the Pacific Ocean and the Indian Ocean (Tan, 2013). We suggest that systematic comparison at different spatial-temporal scales may partially settle those disputes. Then, we selected relatively continuous stalagmites δ¹⁸O records from caves that are far enough to capture spatial variation in the Asian region during the middle to late last glacial for comparison.

The stalagmite δ¹⁸O records in the monsoon region of China are controlled by δ¹⁸O of yearly precipitation (Qiu *et al.* 2021), of which almost 80% is contributed by the EASM and the last 20% is dominated by the East Asian Winter Monsoon (EAWM) in the current climate setting. While, under the background of weakened Asian summer monsoon and intensified winter monsoon during the last glacial (Sun *et al.* 2006), this ratio may change partially. However, in this discussion, we do not explore the changed proportion emphasized by Wang *et al.* (2001) but focus on the averages of these stalagmite δ¹⁸O records over the last glacial period.

Wang *et al.* (2001) explained the changed stalagmite δ¹⁸O records of the Hulu Cave as the ratio between the East Asia summer and winter monsoon, which is very reasonable. Whereas, when comparing averages of stalagmite δ¹⁸O records from different regions over the last glacial period, the enhanced winter winds had a relatively common effect on these caves, although the influences on these caves may differ. It is conceivable that summer precipitation still constitutes the majority of yearly precipitation (Jiang *et al.* 2011; Cai *et al.* 2015) and dominates the composition of stalagmite δ¹⁸O records. Thus, the comparison among averages of stalagmite δ¹⁸O records reflects summer moisture sources and
<table><caption>**Table 1.** U-Th isotopic measurements for Shizhu Cave stalagmite SZ-1 and SZ-3, on MC-ICP-MS</caption><thead><tr><th>Sample</th><th colspan="2"><sup>238</sup>U</th><th colspan="2"><sup>232</sup>Th</th><th>d<sup>234</sup>U</th><th colspan="2"><sup>230</sup>Th/<sup>238</sup>U</th><th colspan="2"><sup>230</sup>Th/<sup>232</sup>Th</th><th>Age (ka ago)</th><th colspan="2">Age (ka ago)</th><th>d<sup>234</sup>U<sub>initial</sub></th></tr><tr><th>Length(cm)</th><th>10<sup>-9</sup>g/g<sup>a</sup></th><th>±10</th><th>10<sup>-9</sup>g/g</th><th>±10</th><th>Measured<sup>a</sup></th><th>Activity<sup>c</sup></th><th>Atomic (x 10<sup>-6</sup>)</th><th>Uncorrected</th><th>Corrected<sup>c,d</sup></th><th>Uncorrected</th><th>Corrected</th><th>Corrected<sup>b</sup></th></tr></thead><tbody><tr><td>SZ1-9.75</td><td>2265.5</td><td>±2.7</td><td>7891.</td><td>±21</td><td>420.3</td><td>±1.6</td><td>0.2656</td><td>±0.0010</td><td>1257.2</td><td>±5.8</td><td>22.34</td><td>±0.10</td><td>22.28</td><td>±0.11</td><td>447.6</td><td>±1.7</td></tr><tr><td>SZ1-13.25</td><td>222.575</td><td>±0.042</td><td>6132.4</td><td>±4.1</td><td>412.00</td><td>±0.34</td><td>0.2943</td><td>±0.0014</td><td>176.11</td><td>±0.82</td><td>25.20</td><td>±0.13</td><td>24.64</td><td>±0.31</td><td>441.67</td><td>±0.53</td></tr><tr><td>SZ1-16.75</td><td>419.104</td><td>±0.049</td><td>6835.8</td><td>±5.2</td><td>401.65</td><td>±0.27</td><td>0.31982</td><td>±0.00051</td><td>323.30</td><td>±0.57</td><td>27.890</td><td>±0.050</td><td>27.56</td><td>±0.17</td><td>434.13</td><td>±0.36</td></tr><tr><td>SZ1-19.25</td><td>222.95</td><td>±0.26</td><td>5719</td><td>±16</td><td>413.1</td><td>±1.7</td><td>0.3398</td><td>±0.0022</td><td>218.4</td><td>±1.6</td><td>29.59</td><td>±0.23</td><td>29.12</td><td>±0.33</td><td>448.5</td><td>±1.9</td></tr><tr><td>SZ1-21.25</td><td>453.404</td><td>±0.060</td><td>3350.5</td><td>±2.5</td><td>419.95</td><td>±0.37</td><td>0.36618</td><td>±0.00041</td><td>817.0</td><td>±1.1</td><td>32.034</td><td>±0.043</td><td>31.887</td><td>±0.085</td><td>459.49</td><td>±0.42</td></tr><tr><td>SZ1-24.75</td><td>319.51</td><td>±0.41</td><td>148.1</td><td>±4.2</td><td>459.7</td><td>±2.0</td><td>0.39521</td><td>±0.00083</td><td>14062</td><td>±403</td><td>33.844</td><td>±0.099</td><td>33.835</td><td>±0.099</td><td>505.8</td><td>±2.2</td></tr><tr><td>SZ1-27.25</td><td>268.801</td><td>±0.062</td><td>1460.6</td><td>±2.3</td><td>425.61</td><td>±0.51</td><td>0.39759</td><td>±0.00062</td><td>1206.4</td><td>±2.7</td><td>35.056</td><td>±0.065</td><td>34.948</td><td>±0.084</td><td>469.73</td><td>±0.58</td></tr><tr><td>SZ1-28.5</td><td>247.687</td><td>±0.022</td><td>1345.4</td><td>±2</td><td>420.88</td><td>±0.33</td><td>0.4024</td><td>±0.0011</td><td>1221.5</td><td>±3.7</td><td>35.69</td><td>±0.11</td><td>35.60</td><td>±0.12</td><td>465.36</td><td>±0.40</td></tr><tr><td>SZ1-30.0</td><td>236.640</td><td>±0.045</td><td>6872.5</td><td>±7.3</td><td>422.37</td><td>±0.83</td><td>0.4166</td><td>±0.0021</td><td>236.5</td><td>±1.2</td><td>37.12</td><td>±0.22</td><td>36.60</td><td>±0.34</td><td>468.3</td><td>±1.0</td></tr><tr><td>SZ1-33.25</td><td>295.01</td><td>±0.37</td><td>285.9</td><td>±2.1</td><td>431.27</td><td>±1.95</td><td>0.42960</td><td>±0.00106</td><td>7309</td><td>±55</td><td>38,191</td><td>±128.04</td><td>38,174</td><td>±128</td><td>480.32</td><td>±2.18</td></tr><tr><td>SZ1-36.75</td><td>253.70</td><td>±0.32</td><td>571.5</td><td>±2.2</td><td>439.72</td><td>±2.02</td><td>0.45151</td><td>±0.00125</td><td>3305</td><td>±15</td><td>40,210</td><td>±149.38</td><td>40,170</td><td>±151</td><td>492.51</td><td>±2.27</td></tr><tr><td>SZ1-39</td><td>3459.60</td><td>±0.93</td><td>27919.0</td><td>±33.8</td><td>435.56</td><td>±0.51</td><td>0.46926</td><td>±0.00047</td><td>959</td><td>±1</td><td>42,253</td><td>±53.81</td><td>42,110</td><td>±90</td><td>490.52</td><td>±0.59</td></tr><tr><td>SZ1-40.5</td><td>234.75</td><td>±0.05</td><td>221.6</td><td>±2.2</td><td>443.43</td><td>±0.51</td><td>0.48017</td><td>±0.00166</td><td>8386</td><td>±87</td><td>43,138</td><td>±180.13</td><td>43,121</td><td>±180</td><td>500.81</td><td>±0.63</td></tr><tr><td>SZ1-43.25</td><td>280.79</td><td>±0.40</td><td>1907.6</td><td>±4.0</td><td>451.75</td><td>±2.36</td><td>0.49576</td><td>±0.00187</td><td>1203</td><td>±5</td><td>44,508</td><td>±221.92</td><td>44,389</td><td>±229</td><td>512.04</td><td>±2.70</td></tr><tr><td>SZ3-1.75</td><td>672.65</td><td>±0.97</td><td>4027.0</td><td>±8.3</td><td>795.7</td><td>±3.1</td><td>0.2386</td><td>±0.0012</td><td>657.2</td><td>±3.4</td><td>15.391</td><td>±0.085</td><td>15.235</td><td>±0.095</td><td>830.9</td><td>±3.2</td></tr><tr><td>SZ3-5.25</td><td>676.3</td><td>±1.0</td><td>5737.0</td><td>±9.9</td><td>778.5</td><td>±3.0</td><td>0.2639</td><td>±0.0013</td><td>513.0</td><td>±2.6</td><td>17.316</td><td>±0.099</td><td>17.12</td><td>±0.12</td><td>817.2</td><td>±3.1</td></tr><tr><td>SZ3-8.75</td><td>844.1</td><td>±1.1</td><td>2227.7</td><td>±4.2</td><td>798.8</td><td>±2.7</td><td>0.28348</td><td>±0.00087</td><td>1771.1</td><td>±6.0</td><td>18.470</td><td>±0.068</td><td>18.363</td><td>±0.071</td><td>841.5</td><td>±2.8</td></tr><tr><td>SZ3-14.0</td><td>869.8</td><td>±1.2</td><td>1129.1</td><td>±3.0</td><td>742.2</td><td>±3.2</td><td>0.30798</td><td>±0.00078</td><td>3912</td><td>±13</td><td>20.925</td><td>±0.072</td><td>20.837</td><td>±0.072</td><td>787.3</td><td>±3.4</td></tr><tr><td>SZ3-21.0</td><td>911.7</td><td>±1.4</td><td>220.6</td><td>±2.0</td><td>722.7</td><td>±3.1</td><td>0.31831</td><td>±0.00063</td><td>21696</td><td>±195</td><td>21.964</td><td>±0.065</td><td>21.891</td><td>±0.065</td><td>768.9</td><td>±3.3</td></tr><tr><td>SZ3-24.5</td><td>910.6</td><td>±1.0</td><td>2003.7</td><td>±4.0</td><td>714.1</td><td>±2.1</td><td>0.32590</td><td>±0.00083</td><td>2442.0</td><td>±7.5</td><td>22.664</td><td>±0.070</td><td>22.562</td><td>±0.072</td><td>761.2</td><td>±2.2</td></tr><tr><td>SZ3-28.0</td><td>933.1</td><td>±1.3</td><td>169.6</td><td>±2.2</td><td>705.6</td><td>±2.5</td><td>0.33723</td><td>±0.00058</td><td>30588</td><td>±406</td><td>23.663</td><td>±0.059</td><td>23.591</td><td>±0.059</td><td>754.3</td><td>±2.7</td></tr><tr><td>SZ3-35.0</td><td>1104.2</td><td>±1.1</td><td>281.6</td><td>±2.1</td><td>687.8</td><td>±1.7</td><td>0.35266</td><td>±0.00049</td><td>22802</td><td>±169</td><td>25.155</td><td>±0.048</td><td>25.082</td><td>±0.048</td><td>738.4</td><td>±1.9</td></tr><tr><td>SZ3-38.5</td><td>891.05</td><td>±0.96</td><td>5684</td><td>±10</td><td>690.5</td><td>±2.0</td><td>0.3656</td><td>±0.0015</td><td>945.0</td><td>±4.1</td><td>26.13</td><td>±0.12</td><td>25.97</td><td>±0.13</td><td>743.2</td><td>±2.1</td></tr><tr><td>SZ3-39.25</td><td>828.71</td><td>±0.13</td><td>9661</td><td>±10</td><td>686.81</td><td>±0.41</td><td>0.36758</td><td>±0.00094</td><td>519.9</td><td>±1.4</td><td>26.358</td><td>±0.075</td><td>26.25</td><td>±0.12</td><td>739.48</td><td>±0.50</td></tr><tr><td>SZ3-41.0</td><td>704.98</td><td>±0.10</td><td>7170.3</td><td>±9.4</td><td>627.78</td><td>±0.33</td><td>0.37541</td><td>±0.00018</td><td>608.57</td><td>±0.84</td><td>28.107</td><td>±0.016</td><td>28.016</td><td>±0.082</td><td>679.30</td><td>±0.39</td></tr><tr><td>SZ3-42.0</td><td>788.20</td><td>±0.84</td><td>3592.0</td><td>±5.7</td><td>617.1</td><td>±2.0</td><td>0.3834</td><td>±0.0012</td><td>1387.0</td><td>±4.8</td><td>29.00</td><td>±0.11</td><td>28.85</td><td>±0.12</td><td>669.6</td><td>±2.1</td></tr><tr><td>SZ3-44.5</td><td>778.91</td><td>±0.81</td><td>3513.7</td><td>±6.4</td><td>641.8</td><td>±1.8</td><td>0.3980</td><td>±0.0013</td><td>1454.8</td><td>±5.4</td><td>29.73</td><td>±0.12</td><td>29.59</td><td>±0.12</td><td>697.8</td><td>±2.0</td></tr><tr><td>SZ3-49.0</td><td>889.45</td><td>±0.98</td><td>116.8</td><td>±2.5</td><td>653.4</td><td>±1.9</td><td>0.41952</td><td>±0.00070</td><td>52681</td><td>±1117</td><td>31.292</td><td>±0.073</td><td>31.221</td><td>±0.073</td><td>713.7</td><td>±2.1</td></tr><tr><td>SZ3-54.0</td><td>950.4</td><td>±1.7</td><td>558.1</td><td>±2.2</td><td>651.1</td><td>±3.3</td><td>0.4362</td><td>±0.0011</td><td>12248</td><td>±54</td><td>32.77</td><td>±0.12</td><td>32.69</td><td>±0.12</td><td>714.1</td><td>±3.6</td></tr></tbody></table>

Analytical errors are 2s of the mean.

a[238U] = [235U] × 137.818 (±0.65‰) (Hiess et al. 2012); d[234U] = [(234U/238U)activity - 1] × 1000.

b d[234U]initial corrected was calculated based on $^{230}$Th age (T), i.e., d[234U]initial = d[234U]measured × e$^{234}$T, and T is corrected age.

c[230Th/238U]activity = 1 - e$^{-12307}$ + (d[234U]measured/1000)[(230/(230 - 234))][1 - e$^{-230 - 234}$]T, where T is the age.

Decay constants are 9.1705 × 10⁻⁶ yr⁻¹ for $^{230}$Th, 2.8221 × 10⁻⁶ yr⁻¹ for $^{234}$U (Cheng et al. 2013) and 1.55125 × 10⁻¹⁰ yr⁻¹ for $^{238}$U (Jaffey et al. 1971).

dAge corrections, relative to chemistry date, were calculated using an estimated atomic $^{230}$Th/$^{232}$Th ratio of 4 (± 2) × 10⁻⁶.

Those are the values for material at secular equilibrium, with the crustal $^{232}$Th/$^{238}$U value of 3.8. The errors are arbitrarily assumed to be 50%.
<table><caption>**Table 1.** U-Th isotopic measurements for Shizhu Cave stalagmite SZ-1 and SZ-3, on MC-ICP-MS</caption><thead><tr><th>Sample</th><th colspan="2"><sup>238</sup>U</th><th colspan="2"><sup>232</sup>Th</th><th>d<sup>234</sup>U</th><th colspan="2"><sup>230</sup>Th/<sup>238</sup>U</th><th colspan="2"><sup>230</sup>Th/<sup>232</sup>Th</th><th>Age (ka ago)</th><th colspan="2">Age (ka ago)</th><th>d<sup>234</sup>U<sub>initial</sub></th></tr><tr><th>Length(cm)</th><th>10<sup>-9</sup>g/g<sup>a</sup></th><th></th><th>10<sup>-9</sup>g/g</th><th></th><th>Measured<sup>a</sup></th><th>Activity<sup>c</sup></th><th></th><th>Atomic (x 10<sup>-6</sup>)</th><th>Uncorrected</th><th></th><th>Corrected<sup>c,d</sup></th><th></th><th>Corrected<sup>b</sup></th></tr></thead><tbody><tr><td>SZ1-9.75</td><td>2265.5</td><td>±2.7</td><td>7891.</td><td>±21</td><td>420.3</td><td>±1.6</td><td>0.2656</td><td>±0.0010</td><td>1257.2</td><td>±5.8</td><td>22.34</td><td>±0.10</td><td>22.28</td><td>±0.11</td><td>447.6</td><td>±1.7</td></tr><tr><td>SZ1-13.25</td><td>222.575</td><td>±0.042</td><td>6132.4</td><td>±4.1</td><td>412.00</td><td>±0.34</td><td>0.2943</td><td>±0.0014</td><td>176.11</td><td>±0.82</td><td>25.20</td><td>±0.13</td><td>24.64</td><td>±0.31</td><td>441.67</td><td>±0.53</td></tr><tr><td>SZ1-16.75</td><td>419.104</td><td>±0.049</td><td>6835.8</td><td>±5.2</td><td>401.65</td><td>±0.27</td><td>0.31982</td><td>±0.00051</td><td>323.30</td><td>±0.57</td><td>27.890</td><td>±0.050</td><td>27.56</td><td>±0.17</td><td>434.13</td><td>±0.36</td></tr><tr><td>SZ1-19.25</td><td>222.95</td><td>±0.26</td><td>5719</td><td>±16</td><td>413.1</td><td>±1.7</td><td>0.3398</td><td>±0.0022</td><td>218.4</td><td>±1.6</td><td>29.59</td><td>±0.23</td><td>29.12</td><td>±0.33</td><td>448.5</td><td>±1.9</td></tr><tr><td>SZ1-21.25</td><td>453.404</td><td>±0.060</td><td>3350.5</td><td>±2.5</td><td>419.95</td><td>±0.37</td><td>0.36618</td><td>±0.00041</td><td>817.0</td><td>±1.1</td><td>32.034</td><td>±0.043</td><td>31.887</td><td>±0.085</td><td>459.49</td><td>±0.42</td></tr><tr><td>SZ1-24.75</td><td>319.51</td><td>±0.41</td><td>148.1</td><td>±4.2</td><td>459.7</td><td>±2.0</td><td>0.39521</td><td>±0.00083</td><td>14062</td><td>±403</td><td>33.844</td><td>±0.099</td><td>33.835</td><td>±0.099</td><td>505.8</td><td>±2.2</td></tr><tr><td>SZ1-27.25</td><td>268.801</td><td>±0.062</td><td>1460.6</td><td>±2.3</td><td>425.61</td><td>±0.51</td><td>0.39759</td><td>±0.00062</td><td>1206.4</td><td>±2.7</td><td>35.056</td><td>±0.065</td><td>34.948</td><td>±0.084</td><td>469.73</td><td>±0.58</td></tr><tr><td>SZ1-28.5</td><td>247.687</td><td>±0.022</td><td>1345.4</td><td>±2</td><td>420.88</td><td>±0.33</td><td>0.4024</td><td>±0.0011</td><td>1221.5</td><td>±3.7</td><td>35.69</td><td>±0.11</td><td>35.60</td><td>±0.12</td><td>465.36</td><td>±0.40</td></tr><tr><td>SZ1-30.0</td><td>236.640</td><td>±0.045</td><td>6872.5</td><td>±7.3</td><td>422.37</td><td>±0.83</td><td>0.4166</td><td>±0.0021</td><td>236.5</td><td>±1.2</td><td>37.12</td><td>±0.22</td><td>36.60</td><td>±0.34</td><td>468.3</td><td>±1.0</td></tr><tr><td>SZ1-33.25</td><td>295.01</td><td>±0.37</td><td>285.9</td><td>±2.1</td><td>431.27</td><td>±1.95</td><td>0.42960</td><td>±0.00106</td><td>7309</td><td>±55</td><td>38,191</td><td>±128.04</td><td>38,174</td><td>±128</td><td>480.32</td><td>±2.18</td></tr><tr><td>SZ1-36.75</td><td>253.70</td><td>±0.32</td><td>571.5</td><td>±2.2</td><td>439.72</td><td>±2.02</td><td>0.45151</td><td>±0.00125</td><td>3305</td><td>±15</td><td>40,210</td><td>±149.38</td><td>40,170</td><td>±151</td><td>492.51</td><td>±2.27</td></tr><tr><td>SZ1-39</td><td>3459.60</td><td>±0.93</td><td>27919.0</td><td>±33.8</td><td>435.56</td><td>±0.51</td><td>0.46926</td><td>±0.00047</td><td>959</td><td>±1</td><td>42,253</td><td>±53.81</td><td>42,110</td><td>±90</td><td>490.52</td><td>±0.59</td></tr><tr><td>SZ1-40.5</td><td>234.75</td><td>±0.05</td><td>221.6</td><td>±2.2</td><td>443.43</td><td>±0.51</td><td>0.48017</td><td>±0.00166</td><td>8386</td><td>±87</td><td>43,138</td><td>±180.13</td><td>43,121</td><td>±180</td><td>500.81</td><td>±0.63</td></tr><tr><td>SZ1-43.25</td><td>280.79</td><td>±0.40</td><td>1907.6</td><td>±4.0</td><td>451.75</td><td>±2.36</td><td>0.49576</td><td>±0.00187</td><td>1203</td><td>±5</td><td>44,508</td><td>±221.92</td><td>44,389</td><td>±229</td><td>512.04</td><td>±2.70</td></tr><tr><td>SZ3-1.75</td><td>672.65</td><td>±0.97</td><td>4027.0</td><td>±8.3</td><td>795.7</td><td>±3.1</td><td>0.2386</td><td>±0.0012</td><td>657.2</td><td>±3.4</td><td>15.391</td><td>±0.085</td><td>15.235</td><td>±0.095</td><td>830.9</td><td>±3.2</td></tr><tr><td>SZ3-5.25</td><td>676.3</td><td>±1.0</td><td>5737.0</td><td>±9.9</td><td>778.5</td><td>±3.0</td><td>0.2639</td><td>±0.0013</td><td>513.0</td><td>±2.6</td><td>17.316</td><td>±0.099</td><td>17.12</td><td>±0.12</td><td>817.2</td><td>±3.1</td></tr><tr><td>SZ3-8.75</td><td>844.1</td><td>±1.1</td><td>2227.7</td><td>±4.2</td><td>798.8</td><td>±2.7</td><td>0.28348</td><td>±0.00087</td><td>1771.1</td><td>±6.0</td><td>18.470</td><td>±0.068</td><td>18.363</td><td>±0.071</td><td>841.5</td><td>±2.8</td></tr><tr><td>SZ3-14.0</td><td>869.8</td><td>±1.2</td><td>1129.1</td><td>±3.0</td><td>742.2</td><td>±3.2</td><td>0.30798</td><td>±0.00078</td><td>3912</td><td>±13</td><td>20.925</td><td>±0.072</td><td>20.837</td><td>±0.072</td><td>787.3</td><td>±3.4</td></tr><tr><td>SZ3-21.0</td><td>911.7</td><td>±1.4</td><td>220.6</td><td>±2.0</td><td>722.7</td><td>±3.1</td><td>0.31831</td><td>±0.00063</td><td>21696</td><td>±195</td><td>21.964</td><td>±0.065</td><td>21.891</td><td>±0.065</td><td>768.9</td><td>±3.3</td></tr><tr><td>SZ3-24.5</td><td>910.6</td><td>±1.0</td><td>2003.7</td><td>±4.0</td><td>714.1</td><td>±2.1</td><td>0.32590</td><td>±0.00083</td><td>2442.0</td><td>±7.5</td><td>22.664</td><td>±0.070</td><td>22.562</td><td>±0.072</td><td>761.2</td><td>±2.2</td></tr><tr><td>SZ3-28.0</td><td>933.1</td><td>±1.3</td><td>169.6</td><td>±2.2</td><td>705.6</td><td>±2.5</td><td>0.33723</td><td>±0.00058</td><td>30588</td><td>±406</td><td>23.663</td><td>±0.059</td><td>23.591</td><td>±0.059</td><td>754.3</td><td>±2.7</td></tr><tr><td>SZ3-35.0</td><td>1104.2</td><td>±1.1</td><td>281.6</td><td>±2.1</td><td>687.8</td><td>±1.7</td><td>0.35266</td><td>±0.00049</td><td>22802</td><td>±169</td><td>25.155</td><td>±0.048</td><td>25.082</td><td>±0.048</td><td>738.4</td><td>±1.9</td></tr><tr><td>SZ3-38.5</td><td>891.05</td><td>±0.96</td><td>5684</td><td>±10</td><td>690.5</td><td>±2.0</td><td>0.3656</td><td>±0.0015</td><td>945.0</td><td>±4.1</td><td>26.13</td><td>±0.12</td><td>25.97</td><td>±0.13</td><td>743.2</td><td>±2.1</td></tr><tr><td>SZ3-39.25</td><td>828.71</td><td>±0.13</td><td>9661</td><td>±10</td><td>686.81</td><td>±0.41</td><td>0.36758</td><td>±0.00094</td><td>519.9</td><td>±1.4</td><td>26.358</td><td>±0.075</td><td>26.25</td><td>±0.12</td><td>739.48</td><td>±0.50</td></tr><tr><td>SZ3-41.0</td><td>704.98</td><td>±0.10</td><td>7170.3</td><td>±9.4</td><td>627.78</td><td>±0.33</td><td>0.37541</td><td>±0.00018</td><td>608.57</td><td>±0.84</td><td>28.107</td><td>±0.016</td><td>28.016</td><td>±0.082</td><td>679.30</td><td>±0.39</td></tr><tr><td>SZ3-42.0</td><td>788.20</td><td>±0.84</td><td>3592.0</td><td>±5.7</td><td>617.1</td><td>±2.0</td><td>0.3834</td><td>±0.0012</td><td>1387.0</td><td>±4.8</td><td>29.00</td><td>±0.11</td><td>28.85</td><td>±0.12</td><td>669.6</td><td>±2.1</td></tr><tr><td>SZ3-44.5</td><td>778.91</td><td>±0.81</td><td>3513.7</td><td>±6.4</td><td>641.8</td><td>±1.8</td><td>0.3980</td><td>±0.0013</td><td>1454.8</td><td>±5.4</td><td>29.73</td><td>±0.12</td><td>29.59</td><td>±0.12</td><td>697.8</td><td>±2.0</td></tr><tr><td>SZ3-49.0</td><td>889.45</td><td>±0.98</td><td>116.8</td><td>±2.5</td><td>653.4</td><td>±1.9</td><td>0.41952</td><td>±0.00070</td><td>52681</td><td>±1117</td><td>31.292</td><td>±0.073</td><td>31.221</td><td>±0.073</td><td>713.7</td><td>±2.1</td></tr><tr><td>SZ3-54.0</td><td>950.4</td><td>±1.7</td><td>558.1</td><td>±2.2</td><td>651.1</td><td>±3.3</td><td>0.4362</td><td>±0.0011</td><td>12248</td><td>±54</td><td>32.77</td><td>±0.12</td><td>32.69</td><td>±0.12</td><td>714.1</td><td>±3.6</td></tr></tbody></table>

Analytical errors are 2s of the mean.

a[238U] = [235U] × 137.818 (±0.65‰) (Hiess et al. 2012); d[234U] = ([234U]/[238U]activity - 1) × 1000.

b[234U]$_{initial}$ corrected was calculated based on $^{230}$Th age (T), i.e., d[234U]$_{initial}$ = d[234U]$_{measured}$ × e[234T], and T is corrected age.

c[230Th/238U]$_{activity}$ = 1 - e[-12307 + (d[234U]$_{measured}$/1000)][[230]/([230 - [234]])(1 - e[-(230 - [234])T]), where T is the age.

Decay constants are 9.1705 × 10⁻⁶ yr⁻¹ for $^{230}$Th, 2.8221 × 10⁻⁶ yr⁻¹ for $^{234}$U (Cheng et al. 2013) and 1.55125 × 10⁻¹⁰ yr⁻¹ for $^{238}$U (Jaffey et al. 1971).

dAge corrections, relative to chemistry date, were calculated using an estimated atomic $^{230}$Th/$^{232}$Th ratio of 4 (± 2) × 10⁻⁶.

Those are the values for material at secular equilibrium, with the crustal $^{232}$Th/$^{238}$U value of 3.8. The errors are arbitrarily assumed to be 50%.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** (Colour online) δ¹⁸O records from the stalagmites SZ-1 and SZ-3 investigated in this study. The original δ¹⁸O records from (a) SZ-1 (forest green) and (b) SZ-3 (purple). (d) Composite Shizhu δ¹⁸O record SZ-1-3 (navy blue). Error bars of the ²³⁰Th dating points are shown at the top for both stalagmites. The light yellow bar in (a) highlights the stage lacking dating results before 22.28 ka and the possible recrystallization during 22.28 ka to 24.64 ka B.P. of SZ-1. The black lines in (d) are error curves of stack analysis. The light pink and green rectangles present the periods of MIS 2 and MIS 3.</figcaption></figure>

trajectories during the middle to late last glacial period based on Rayleigh fractionation (Dansgaard, 1964).

Hu *et al.* (2008) compared stalagmite δ¹⁸O records from two caves along the same moisture transport path to quantify the history of rainfall during the Holocene, verifying the effectiveness of comparing records from different regions based on spatially coherent variability. Further, Li *et al.* (2017) demonstrated the different moisture sources and paths of southwest and east China by comparing stalagmite δ¹⁸O records from different regions. In the ISM region, the spatial variation of stalagmite δ¹⁸O records from coast to inland has been demonstrated on a subregional scale in south and Southeast Asia by Liu *et al.* (2020) over the past 40 ka.

In the EASM region, Zhang *et al.* (2021) demonstrated that the Yongxing Cave was located downstream of the EASM moisture path during the last glacial. These previous studies inspire us but lack systematic comparison on larger spatial scales.

The graph in Figure 5 illustrates how the average isotope values change over time. The Shizhu Cave is located at the main path of the ISM moisture trajectory whether in the present day or during the last glacial (Fig. 1; Fig. 6). Therefore, SZ-1-3 records are used to compare with those from sites along the ISM moisture trajectory (i.e., Mawmluh and Yongxing Caves) between 45 and 15 ka B.P. (Fig. 5). Their average values in Figure 5a show systematically decrease from the Indian Ocean to inland China. Specifically, the
<figure><img src="image_5.png" /><figcaption>**Figure 4.** (Colour online) δ¹⁸O records from the stalagmites SZ-1 and SZ-3 investigated in this study. The original δ¹⁸O records from (a) SZ-1 (forest green) and (b) SZ-3 (purple). (d) Composite Shizhu δ¹⁸O record SZ-1-3 (navy blue). Error bars of the ²³⁰Th dating points are shown at the top for both stalagmites. The light yellow bar in (a) highlights the stage lacking dating results before 22.28 ka and the possible recrystallization during 22.28 ka to 24.64 ka B.P. of SZ-1. The black lines in (d) are error curves of stack analysis. The light pink and green rectangles present the periods of MIS 2 and MIS 3.</figcaption></figure>

trajectories during the middle to late last glacial period based on Rayleigh fractionation (Dansgaard, 1964).

Hu *et al.* (2008) compared stalagmite δ¹⁸O records from two caves along the same moisture transport path to quantify the history of rainfall during the Holocene, verifying the effectiveness of comparing records from different regions based on spatially coherent variability. Further, Li *et al.* (2017) demonstrated the different moisture sources and paths of southwest and east China by comparing stalagmite δ¹⁸O records from different regions. In the ISM region, the spatial variation of stalagmite δ¹⁸O records from coast to inland has been demonstrated on a subregional scale in south and Southeast Asia by Liu *et al.* (2020) over the past 40 ka.

In the EASM region, Zhang *et al.* (2021) demonstrated that the Yongxing Cave was located downstream of the EASM moisture path during the last glacial. These previous studies inspire us but lack systematic comparison on larger spatial scales.

The graph in Figure 5 illustrates how the average isotope values change over time. The Shizhu Cave is located at the main path of the ISM moisture trajectory whether in the present day or during the last glacial (Fig. 1; Fig. 6). Therefore, SZ-1-3 records are used to compare with those from sites along the ISM moisture trajectory (i.e., Mawmluh and Yongxing Caves) between 45 and 15 ka B.P. (Fig. 5). Their average values in Figure 5a show systematically decrease from the Indian Ocean to inland China. Specifically, the
<figure><img src="image_6.png" /><figcaption>**Figure 5.** (Colour online) Comparison of δ¹⁸O records on ISM and EASM moisture trajectories. (a) Mawmluh (pink), Shizhu (navy blue; this study) and Yongxing (purple) Cave records. (b) Hulu (green), Yongxing (purple) and Songjia (black) Cave records. The light yellow bars highlight the four weak monsoon periods in the Asian monsoon regions, coeval with the Heinrich stadials 1 to 4 (Wang *et al.* 2001). The light pink and green rectangles are the same as in Figure 4.</figcaption></figure>

stalagmite δ¹⁸O records of Mawmluh Cave are significantly enriched because it is close to the moisture source, the Indian Ocean (Dutt *et al.* 2015). In contrast, the δ¹⁸O records in the Yongxing Cave are considerably lower than those of the Shizhu and Mawmluh Caves because of its midstream or downstream location along the moisture trajectory of the ISM (yellow arrows in Figure 6; Hu *et al.* 2008; Cai *et al.* 2015; Chen *et al.* 2016). We demonstrate here that water vapour fractionation along this coastal-inland moisture transport route from the Indian Ocean to the coastal Mawmluh Cave, then to the Shizhu Cave, and probably extended to the Yongxing Cave during 45 and 15 ka (Fig. 6).

In the eastern monsoonal region of China, the coastal Xianyun and Hulu Caves received more enriched δ¹⁸O of precipitation from the adjacent tropical and western Pacific Oceans (Wang *et al.* 2001;

Zhang *et al.* 2021; Qiu *et al.* 2022). From the coastal Hulu/Xianyun Caves to the western inland Yongxing and then Songjia Caves, the average stalagmite δ¹⁸O values gradually decrease (Fig. 5b). We suggest that water vapour fractionation along this coastal-inland moisture transport route from the Pacific Ocean to the Yongxing Cave can be extended to the Songjia Cave during 45 and 15 ka (orange arrows in Fig. 6). Additionally, from coastal Hulu/Xianyun Caves to Dragon records in northern China (Dong *et al.* 2018) probably imply a northern moisture trajectory. Overall, we propose that the two moisture trajectories from the Indian Ocean and the Pacific Ocean may exist during 45–15 ka B.P., just as similar to the climatology of modelled summer water vapour transport during the last glacial period (Cai *et al.* 2015; Fig. 6).
<figure><img src="image_6.png" /><figcaption>**Figure 5.** (Colour online) Comparison of δ¹⁸O records on ISM and EASM moisture trajectories. (a) Mawmluh (pink), Shizhu (navy blue; this study) and Yongxing (purple) Cave records. (b) Hulu (green), Yongxing (purple) and Songjia (black) Cave records. The light yellow bars highlight the four weak monsoon periods in the Asian monsoon regions, coinciding with the Heinrich stadials 1 to 4 (Wang *et al.* 2001). The light pink and green rectangles are the same as in Figure 4.</figcaption></figure>

stalagmite δ¹⁸O records of Mawmluh Cave are significantly enriched because it is close to the moisture source, the Indian Ocean (Dutt *et al.* 2015). In contrast, the δ¹⁸O records in the Yongxing Cave are considerably lower than those of the Shizhu and Mawmluh Caves because of its midstream or downstream location along the moisture trajectory of the ISM (yellow arrows in Figure 6; Hu *et al.* 2008; Cai *et al.* 2015; Chen *et al.* 2016). We demonstrate here that water vapour fractionation along this coastal-inland moisture transport route from the Indian Ocean to the coastal Mawmluh Cave, then to the Shizhu Cave, and probably extended to the Yongxing Cave during 45 and 15 ka (Fig. 6).

In the eastern monsoonal region of China, the coastal Xianyun and Hulu Caves received more enriched δ¹⁸O of precipitation from the adjacent tropical and western Pacific Oceans (Wang *et al.* 2001;

Zhang *et al.* 2021; Qiu *et al.* 2022). From the coastal Hulu/Xianyun Caves to the western inland Yongxing and then Songjia Caves, the average stalagmite δ¹⁸O values gradually decrease (Fig. 5b). We suggest that water vapour fractionation along this coastal-inland moisture transport route from the Pacific Ocean to the Yongxing Cave can be extended to the Songjia Cave during 45 and 15 ka (orange arrows in Fig. 6). Additionally, from coastal Hulu/Xianyun Caves to Dragon records in northern China (Dong *et al.* 2018) probably imply a northern moisture trajectory. Overall, we propose that the two moisture trajectories from the Indian Ocean and the Pacific Ocean may exist during 45–15 ka B.P., just as similar to the climatology of modelled summer water vapour transport during the last glacial period (Cai *et al.* 2015; Fig. 6).
<figure><img src="image_7.png" /><figcaption>**Figure 6.** (Colour online) Possible ISM (yellow arrows) and EASM (orange arrows) moisture trajectories during the middle to late last glacial. The base map is modified from Cai *et al.* (2015), showing modelled June-August vapour transport (arrows, $m^{-1}skgm^{-2}$) and isotopic composition of column integrated vapour (colour-shading, ‰) during the last glacial maximum. Dark grey lines indicate the present coastline; the blue lines and hatched area indicate the coastline and the exposed continental shelf during the LGM when the sea level was ~120 m lower than the present day. Legends in the left bottom indicate distribution intervals of the averages of stalagmite $\delta^{18}$O records during the middle to late last glacial.</figcaption></figure>

### 4.b. Spatial pattern of millennial-scale variations of Asian stalagmite $\delta^{18}$O records

Millennial-scale abrupt climate changes also occur during 45–15 ka B.P., that is, Heinrich stadials. To highlight these Heinrich stadials, the Hulu, Yongxing, Dragon, Mawmluh and Shizhu $\delta^{18}$O records in the Asian region in Figure 7a were selected according to their data availability during the middle to late last glacial and far distances among them. We highlighted their changes during HS1 to HS4 after minus their respective averages of $\delta^{18}$O records in Figure 7a. We further stacked them using the bootstrap resampling method to avoid visual clutter in Figure 7b. The stacked curve shows rather consistent changes during HS1 to HS4. This consistent strongly suggests that they are dominated simultaneously by the same mechanism – integral changes in atmospheric-oceanic circulation patterns of the East Asian Monsoon (Wang *et al.* 2001), which is further linked to the paleoclimatic changes in the North Atlantic region (Fig. 7c, d).

Abrupt cooling at the North Atlantic region during the Heinrich stadials presents a north-south propagation of climatic signal initiated by the intrusion of vast amounts of freshwater into the Atlantic from melting continental ice sheets, which have been proved by Ice Raft Debris (IRD in Fig. 7c). This causes a disruption of the meridional overturning circulation (Bohm *et al.* 2015; Markle *et al.* 2016; Lynch-Stieglitz, 2017; Dong *et al.* 2022) and further a reduction of the northward heat transport, leading to a cooling in the middle and high-latitude North Atlantic (Fig. 7d) and warming of the Southern Hemisphere (Zhang & Delworth, 2005). This distinct change in heat balance between both hemispheres causes anomalously southward migration of the ITCZ (Chiang & Bitz, 2005; Broccoli *et al.* 2006), severely altering the spatial precipitation pattern in the tropics and subtropics (Jacobel *et al.* 2016), including the Asian summer monsoon region.

Specifically, the Asian summer monsoon consists of tropical ISM and hybrid EASM, which in response to Heinrich stadials different in detail. From the perspective of ISM, the ITCZ shifts southward, leading to a relatively dry and cool climate over the northern Indian Ocean (Dahl *et al.* 2005; Clement & Peterson,

2008; Mohtadi *et al.* 2014; Zorzi *et al.* 2015; Tierney *et al.* 2016; Bradley *et al.* 2017) and decreasing rainfall and river runoff in the south Asian subcontinent further (Mohtadi *et al.* 2014; Deplazes *et al.* 2014). Simulations also show that the sudden increase in North Atlantic sea-ice extent in the HS1 impacts the Indian Ocean climate and causes a sudden increase in $\delta^{18}$O$_P$ over the Indian subcontinent, as well as isotopically heavier water vapour exported eastward into Southeast Asia (Pausata *et al.* 2011).

Overall, it is commonly agreed that climate changes in the Indian Ocean are the primary cause of enriched stalagmite $\delta^{18}$O records in the ISM regions from HS1 to HS4. The change in water vapour in the Indian Ocean and the Indian subcontinent to what extent affects the EASM region is still controversial. As described above, Pausata *et al.* (2011) proposed that the changed water vapour in the Indian Ocean leads to enriched $\delta^{18}$O$_P$ in the ISM region. These two mutually account for $\delta^{18}$O$_P$ of the downstream EASM region. That is to say, variations of stalagmite $\delta^{18}$O records in Chinese caves during Heinrich stadials are almost entirely affected by the ISM rather than the EASM.

However, recent studies suggest that the meridional position of the westerlies also affects the precipitation seasonality, and consequently, the $\delta^{18}$O of precipitation and stalagmite records in the EASM region (Chiang *et al.* 2015, 2017, 2020; Liang *et al.* 2020). For example, Liang *et al.* (2020) suggested that 28°N is a boundary line for stalagmite $\delta^{18}$O records, with higher amplitude in the north (>1.5‰) and lower amplitude in the south (~1‰) during HS1. The impediment of the westerly jet stream caused a small fraction of maritime moisture transported into inland China (less $\delta^{18}$O-depleted precipitation), but a large fraction (more $\delta^{18}$O-depleted precipitation) can still reach the lower latitude regions under modern settings and the Heinrich stadials. Thus, a mixture of yearly rainfall could cause small amplitude changes in calcite $\delta^{18}$O and regional differences in the north and south of the boundary line.

Furtherly, to detect factors affecting stalagmite $\delta^{18}$O records in the EASM region, we took amplitudes, which are defined as the departures apart from the background, of existing stalagmite $\delta^{18}$O
<figure><img src="image_7.png" /><figcaption>**Figure 6.** (Colour online) Possible ISM (yellow arrows) and EASM (orange arrows) moisture trajectories during the middle to late last glacial. The base map is modified from Cai *et al.* (2015), showing modelled June-August vapour transport (arrows, $m^{-1}skgm^{-2}$) and isotopic composition of column integrated vapour (colour-shading, ‰) during the last glacial maximum. Dark grey lines indicate the present coastline; the blue lines and hatched area indicate the coastline and the exposed continental shelf during the LGM when the sea level was ~120 m lower than the present day. Legends in the left bottom indicate distribution intervals of the averages of stalagmite $\delta^{18}$O records during the middle to late last glacial.</figcaption></figure>

### 4.b. Spatial pattern of millennial-scale variations of Asian stalagmite $\delta^{18}$O records

Millennial-scale abrupt climate changes also occur during 45–15 ka B.P., that is, Heinrich stadials. To highlight these Heinrich stadials, the Hulu, Yongxing, Dragon, Mawmluh and Shizhu $\delta^{18}$O records in the Asian region in Figure 7a were selected according to their data availability during the middle to late last glacial and far distances among them. We highlighted their changes during HS1 to HS4 after minus their respective averages of $\delta^{18}$O records in Figure 7a. We further stacked them using the bootstrap resampling method to avoid visual clutter in Figure 7b. The stacked curve shows rather consistent changes during HS1 to HS4. This consistent strongly suggests that they are dominated simultaneously by the same mechanism – integral changes in atmospheric-oceanic circulation patterns of the East Asian Monsoon (Wang *et al.* 2001), which is further linked to the paleoclimatic changes in the North Atlantic region (Fig. 7c, d).

Abrupt cooling at the North Atlantic region during the Heinrich stadials presents a north-south propagation of climatic signal initiated by the intrusion of vast amounts of freshwater into the Atlantic from melting continental ice sheets, which have been proved by Ice Raft Debris (IRD in Fig. 7c). This causes a disruption of the meridional overturning circulation (Bohm *et al.* 2015; Markle *et al.* 2016; Lynch-Stieglitz, 2017; Dong *et al.* 2022) and further a reduction of the northward heat transport, leading to a cooling in the middle and high-latitude North Atlantic (Fig. 7d) and warming of the Southern Hemisphere (Zhang & Delworth, 2005). This distinct change in heat balance between both hemispheres causes anomalously southward migration of the ITCZ (Chiang & Bitz, 2005; Broccoli *et al.* 2006), severely altering the spatial precipitation pattern in the tropics and subtropics (Jacobel *et al.* 2016), including the Asian summer monsoon region.

Specifically, the Asian summer monsoon consists of tropical ISM and hybrid EASM, which in response to Heinrich stadials different in detail. From the perspective of ISM, the ITCZ shifts southward, leading to a relatively dry and cool climate over the northern Indian Ocean (Dahl *et al.* 2005; Clement & Peterson,

2008; Mohtadi *et al.* 2014; Zorzi *et al.* 2015; Tierney *et al.* 2016; Bradley *et al.* 2017) and decreasing rainfall and river runoff in the south Asian subcontinent further (Mohtadi *et al.* 2014; Deplazes *et al.* 2014). Simulations also show that the sudden increase in North Atlantic sea-ice extent in the HS1 impacts the Indian Ocean climate and causes a sudden increase in $\delta^{18}$O$_P$ over the Indian subcontinent, as well as isotopically heavier water vapour exported eastward into Southeast Asia (Pausata *et al.* 2011).

Overall, it is commonly agreed that climate changes in the Indian Ocean are the primary cause of enriched stalagmite $\delta^{18}$O records in the ISM regions from HS1 to HS4. The change in water vapour in the Indian Ocean and the Indian subcontinent to what extent affects the EASM region is still controversial. As described above, Pausata *et al.* (2011) proposed that the changed water vapour in the Indian Ocean leads to enriched $\delta^{18}$O$_P$ in the ISM region. These two mutually account for $\delta^{18}$O$_P$ of the downstream EASM region. That is to say, variations of stalagmite $\delta^{18}$O records in Chinese caves during Heinrich stadials are almost entirely affected by the ISM rather than the EASM.

However, recent studies suggest that the meridional position of the westerlies also affects the precipitation seasonality, and consequently, the $\delta^{18}$O of precipitation and stalagmite records in the EASM region (Chiang *et al.* 2015, 2017, 2020; Liang *et al.* 2020). For example, Liang *et al.* (2020) suggested that 28°N is a boundary line for stalagmite $\delta^{18}$O records, with higher amplitude in the north (>1.5‰) and lower amplitude in the south (~1‰) during HS1. The impediment of the westerly jet stream caused a small fraction of maritime moisture transported into inland China (less $\delta^{18}$O-depleted precipitation), but a large fraction (more $\delta^{18}$O-depleted precipitation) can still reach the lower latitude regions under modern settings and the Heinrich stadials. Thus, a mixture of yearly rainfall could cause small amplitude changes in calcite $\delta^{18}$O and regional differences in the north and south of the boundary line.

Furtherly, to detect factors affecting stalagmite $\delta^{18}$O records in the EASM region, we took amplitudes, which are defined as the departures apart from the background, of existing stalagmite $\delta^{18}$O
<figure><img src="image_8.png" /><figcaption>**Figure 7.** (Colour online) Response of stalagmite δ¹⁸O records to IRD and SST in North Atlantic. (a) Records from the Mawmluh, Shizhu, Yongxing, Hulu and Dragon Caves after minus their respective average values; (b) the stacked curve of these five caves which is calculated by bootstrap resampling method; (c–d) IRD and SST records are referred to Lisiecki & Stern (2016).</figcaption></figure>

records (shown in Fig. S3) widely in Asian summer monsoon during HS1 to HS4 for comparison by the method in Fig. S4. As shown in Tab. S1 and Figure 8, amplitudes of stalagmite δ¹⁸O records from different caves increase from coastal oceans to mid-latitude inland China, presenting a ‘coastal-inland’ pattern, which is probably associated with two types of explanations.

First, this is possibly interpreted by the increased continental climate in central China, which concerns decreased summer rainfall amounts. Changing North Atlantic surface water strengthens the EAWM and weakens the summer monsoon probably through the northern westerlies (Porter & An, 1995; Sun et al. 2012). The southern shifts of the ITCZ in the tropical ocean (Mohtadi et al. 2014) and the South China Sea (Kaboth-Bahr et al. 2021) together with exposed shelves in the Asian region caused less maritime moisture transported into inland China, which means less δ¹⁸O-depleted precipitation in a year as Wang et al. (2001)

proposed. Chiang *et al.* (2015) also suggested that a southward movement of the westerly jet could prevent the inland intrusion of the low-level monsoonal flow and Asian summer monsoon moisture. Eventually, the summer seasonality is gradually reducing in central China with increasing distance to oceans, resulting in increasing amplitudes of stalagmite δ¹⁸O records during Heinrich stadials.

Another possibility is the increased proportion of summer moisture resources from the Pacific Ocean. Tan (2009, 2011, 2013, 2016) proposed that the position of the WPSH affects this proportion in summer monsoonal China. The WPSH is further associated with the global atmospheric and oceanic circulation. Specifically, a strong WPSH at a more southwestern location may weaken the ISM over the Bay of Bengal and strengthen the EASM over the South China Sea and the Western North Pacific, and vice versa (e.g., Gong & Ho, 2002; Zhou et al. 2009; Cao et al. 2012,
<figure><img src="image_8.png" /><figcaption>**Figure 7.** (Colour online) Response of stalagmite δ¹⁸O records to IRD and SST in North Atlantic. (a) Records from the Mawmluh, Shizhu, Yongxing, Hulu and Dragon Caves after minus their respective average values; (b) the stacked curve of these five caves which is calculated by bootstrap resampling method; (c–d) IRD and SST records are referred to Lisiecki & Stern (2016).</figcaption></figure>

records (shown in Fig. S3) widely in Asian summer monsoon during HS1 to HS4 for comparison by the method in Fig. S4. As shown in Tab. S1 and Figure 8, amplitudes of stalagmite δ¹⁸O records from different caves increase from coastal oceans to mid-latitude inland China, presenting a ‘coastal-inland’ pattern, which is probably associated with two types of explanations.

First, this is possibly interpreted by the increased continental climate in central China, which concerns decreased summer rainfall amounts. Changing North Atlantic surface water strengthens the EAWM and weakens the summer monsoon probably through the northern westerlies (Porter & An, 1995; Sun et al. 2012). The southern shifts of the ITCZ in the tropical ocean (Mohtadi et al. 2014) and the South China Sea (Kaboth-Bahr et al. 2021) together with exposed shelves in the Asian region caused less maritime moisture transported into inland China, which means less δ¹⁸O-depleted precipitation in a year as Wang et al. (2001)

proposed. Chiang *et al.* (2015) also suggested that a southward movement of the westerly jet could prevent the inland intrusion of the low-level monsoonal flow and Asian summer monsoon moisture. Eventually, the summer seasonality is gradually reducing in central China with increasing distance to oceans, resulting in increasing amplitudes of stalagmite δ¹⁸O records during Heinrich stadials.

Another possibility is the increased proportion of summer moisture resources from the Pacific Ocean. Tan (2009, 2011, 2013, 2016) proposed that the position of the WPSH affects this proportion in summer monsoonal China. The WPSH is further associated with the global atmospheric and oceanic circulation. Specifically, a strong WPSH at a more southwestern location may weaken the ISM over the Bay of Bengal and strengthen the EASM over the South China Sea and the Western North Pacific, and vice versa (e.g., Gong & Ho, 2002; Zhou et al. 2009; Cao et al. 2012,
<figure><img src="image_9.png" /><figcaption>**Figure 8.** (Colour online) Changing amplitudes of caves during Heinrich stadials 1 to 4. The pink and orange arrows indicate the positions of ISM and EASM. The information on the base map is the same as in Figure 6.</figcaption></figure>

2016). The enriched stalagmite δ¹⁸O records of caves in southeast China during HS1 to HS4 were possibly caused by an increased proportion of moisture from the Pacific Ocean. As to the ‘coastal-inland’ pattern, Zhao *et al.* (2019) proposed that caves in other monsoon regions of China with dominant moisture sources from either the Pacific or the Indian Ocean are not significantly correlated with WPSH variations. Considering this, we suggest that the caves located in central China where influenced by both ISM and EASM are most affected by the changed proportion, showing the biggest amplitudes in the southeast China region during Heinrich stadials.

Additionally, the distribution pattern of HS2 (Fig. 8b) seems to differ from the last three, without an obvious ‘coastal-inland’ pattern. While records in the EASM region in Figure 8b show a less obvious ‘south-north’ pattern in mid-low latitudes. Some researchers suggest that this pattern is possibly caused by the migration of the westerlies during Heinrich stadials directly (Chiang *et al.* 2015, 2017, 2020; Nagashima *et al.* 2011). The specific mechanisms have been described above by the Liang *et al.* (2020). Besides, it is noted that the Dragon Cave in northern China in Figure 8d presents smaller amplitudes, compared with these caves in central China. A previous study of backward trajectory analysis in northern China reveals the importance of westerly

borne and near-sourced moisture, dominating 77% of the summertime rainfall (Draxler & Hess, 1998). This probably leads to a weaker source effect and smaller amplitudes of the Dragon Cave during Heinrich stadials compared with caves in central China.

## Conclusion

To further explore the moisture transportation process from oceans to the continent over China in the past, we present stalagmite δ¹⁸O records between 45 and 15 ka B.P. from the Shizhu Cave, southwestern China. After comparing it with other Asian δ¹⁸O records, we propose two main moisture trajectories: one from the Indian Ocean, through the Shizhu Cave to central China, and one from the Pacific Ocean to central China. This gradually systematic decreasing average values of stalagmite δ¹⁸O from oceans to inland China present a spatial pattern of water vapour fractionation during the middle to late last glacial. Further, to explore factors affecting stalagmite δ¹⁸O records in the EASM region during Heinrich stadials, we took amplitudes of stalagmite δ¹⁸O records widely in Asian summer monsoon during HS1 to HS4 for comparison. The amplitude results show an increasing trend from coastal oceans to mid-latitude inland China, presenting a
<figure><img src="image_9.png" /><figcaption>**Figure 8.** (Colour online) Changing amplitudes of caves during Heinrich stadials 1 to 4. The pink and orange arrows indicate the positions of ISM and EASM. The information on the base map is the same as in Figure 6.</figcaption></figure>

2016). The enriched stalagmite δ¹⁸O records of caves in southeast China during HS1 to HS4 were possibly caused by an increased proportion of moisture from the Pacific Ocean. As to the ‘coastal-inland’ pattern, Zhao *et al.* (2019) proposed that caves in other monsoon regions of China with dominant moisture sources from either the Pacific or the Indian Ocean are not significantly correlated with WPSH variations. Considering this, we suggest that the caves located in central China where influenced by both ISM and EASM are most affected by the changed proportion, showing the biggest amplitudes in the southeast China region during Heinrich stadials.

Additionally, the distribution pattern of HS2 (Fig. 8b) seems to differ from the last three, without an obvious ‘coastal-inland’ pattern. While records in the EASM region in Figure 8b show a less obvious ‘south-north’ pattern in mid-low latitudes. Some researchers suggest that this pattern is possibly caused by the migration of the westerlies during Heinrich stadials directly (Chiang *et al.* 2015, 2017, 2020; Nagashima *et al.* 2011). The specific mechanisms have been described above by the Liang *et al.* (2020). Besides, it is noted that the Dragon Cave in northern China in Figure 8d presents smaller amplitudes, compared with these caves in central China. A previous study of backward trajectory analysis in northern China reveals the importance of westerly

borne and near-sourced moisture, dominating 77% of the summertime rainfall (Draxler & Hess, 1998). This probably leads to a weaker source effect and smaller amplitudes of the Dragon Cave during Heinrich stadials compared with caves in central China.

## Conclusion

To further explore the moisture transportation process from oceans to the continent over China in the past, we present stalagmite δ¹⁸O records between 45 and 15 ka B.P. from the Shizhu Cave, southwestern China. After comparing it with other Asian δ¹⁸O records, we propose two main moisture trajectories: one from the Indian Ocean, through the Shizhu Cave to central China, and one from the Pacific Ocean to central China. This gradually systematic decreasing average values of stalagmite δ¹⁸O from oceans to inland China present a spatial pattern of water vapour fractionation during the middle to late last glacial. Further, to explore factors affecting stalagmite δ¹⁸O records in the EASM region during Heinrich stadials, we took amplitudes of stalagmite δ¹⁸O records widely in Asian summer monsoon during HS1 to HS4 for comparison. The amplitude results show an increasing trend from coastal oceans to mid-latitude inland China, presenting a
‘coastal-inland’ pattern, which is probably associated with the decreased summer rainfall amount in inland areas or/and the increased proportion of summer moisture resources from the Pacific Ocean.

**Supplementary material.** To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756824000013

**Acknowledgements.** This research was supported by grants from the National Natural Science Foundation of China (42074071, 42071106, 41874078, 92158208, 42274094), Shenzhen Science and Technology Program (grant no. KQTD20170810111725321), Science and Technology Innovation Committee of Shenzhen Municipality (ZDSYS201802081843490), Technology and Innovation Commission of Shenzhen Municipality (20200925154739001), Southern University of Science and Technology (Y01316111), and Key Special Project for Introduced Talents Team of Southern Marine Science, the International Cooperation Project of Fujian (202210013 to X. J) and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0210). U-Th dating was partially supported by grants from the National Science and Technology Council (111-2116-M-002-022-MY3), the National Taiwan University (112L894202 to C.-C.S.), and the Higher Education Sprout Project of the Ministry of Education (110L901001 and 110L8907). Shanghai Sheshan National Geophysical Observatory (SSKP202102).

**Competing interests.** We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work, and there is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, the manuscript entitled, ‘Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period’.

## References

* Andersen KK (2004) High-resolution record of northern hemisphere climate extending into the last interglacial period. *Nature* **431**, 147–51.

* Baker AJ, Sodemann H, Baldini JUL, Breitenbach SFM, Johnson KR, Hunen JV and Zhang PZ (2015) Seasonality of westerly moisture transport in the East Asian summer monsoon and its implications for interpreting precipitation δ¹⁸O. *Journal of Geophysical Research: Atmospheres* **120**, 5850–62.

* Bohm E, Lippold J, Gutjahr M, Frank M, Blaser P, Antz B, Fohlmeister J, Frank N, Andersen MB and Deininger M (2015) Strong and deep Atlantic meridional overturning circulation during the last glacial cycle. *Nature* **517**, 73–6.

* Bradley RM, Eric JS, Christo B, Spruce WS, Cecilia MB, Fudge TJ, Joel BP, Ding QH, Tyler RJ, Jones WCW and Todd S (2017) Global atmospheric teleconnections during Dansgaard-Oeschger events. *Nature Geoscience* **10**, 36–40.

* Broccoli AJ, Dahl KA and Stouffer RJ (2006) Response of the ITCZ to Northern Hemisphere cooling. *Geophysical Research Letters* **33**, L01702.

* Cai YJ, Fung IY, Edwards RL, An ZS, Cheng H, Lee JE, Tan LC, Shen C-C, Wang XF, Day JA, Zhou WJ, Kelly MJ and Chiang JC (2015) Variability of stalagmite-inferred Indian monsoon precipitation over the past 252,000 y. *Proceedings of the National Academy of Sciences of the United States of America* **112**, 2954–9.

* Cao J, Gui S, Su Q and Yang Y (2016) The variability of the Indian-East Asian summer monsoon interface in relation to the spring seesaw mode between the Indian Ocean and the central-western Pacific. *Journal of Climate* **29**, 5027–40.

* Cao J, Hu J and Tao Y (2012) An index for the interface between the Indian summer monsoon and the East Asian summer monsoon. *Journal of Geophysical Research: Atmospheres* **117**, 108.

* Chen FH, Chen XM, Chen JH, Zhou AF, Wu D, Tang LY, Zhang XJ, Huang XZ and Yu JQ (2015) Holocene vegetation history, precipitation changes, and Indian summer monsoon evolution documented from sediments of Xingyun Lake, southwest China. *Journal of Quaternary Science* **29**, 661–74.

* Chen ST, Wang YJ, Cheng H, Edwards RL, Wang XF, Kong XG and Liu DB (2016) Strong coupling of Asian Monsoon and Antarctic climates on suborbital timescales. *Scientific Report* **8**, 32995.

* Cheng H, Edwards RL, Shen C-C, Polyak VJ, Asmerom Y, Woodhead J, Hellstrom J, Wang YJ, Kong XG, Spötl C, Wang XF and Alexander JEC (2013) Improvements in $^{230}$Th dating, $^{230}$Th and $^{234}$U half-life values, and U-Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry. *Earth and Planetary Science Letters* **371–372**, 82–91.

* Cheng H, Edwards RL, Sinha A, Spötl C, Yi L, Chen ST, Kelly M, Kathayat G, Wang XF, Li XL, Kong XG, Wang YJ, Ning YF and Zhang HW (2016) The Asian monsoon over the past 640,000 years and ice age terminations. *Nature* **534**, 640–6.

* Cheng H, Sinha A, Wang XF, Cruz FW and Edwards RL (2012) The Global Paleomonsoon is seen through speleothem records from Asia and the Americas. *Climate Dynamics* **39**, 1045–62.

* Cheng H, Zhang HW, Zhao JY, Li HY, Ning YF and Kathayat G (2019) Chinese stalagmite paleoclimate researches: a review and perspective. *Science China Earth Sciences* **62**, 1–25.

* Chiang JCH and Bitz CM (2005) Influence of high latitude ice cover on the marine Intertropical Convergence Zone. *Climate Dynamics* **25**, 477–96.

* Chiang JCH, Fung IY, Wu CH, Cai YJ, Edman JP, Liu YW, Day JA, Bhattacharya T, Mondal Y and Labrousse CA (2015) Role of seasonal transitions and westerly jets in East Asian paleoclimate. *Quaternary Science Reviews* **108**, 111–29.

* Chiang JCH, Herman MJ, Yoshimura K and Fung IY (2020) Enriched East Asian oxygen isotope of precipitation indicates reduced summer seasonality in regional climate and westerlies. *Proceedings of the National Academy of Sciences of the United States of America* **117**, 14745–50.

* Chiang JCH, Swenson LM and Kong W (2017) Role of seasonal transitions and the westerlies in the interannual variability of the East Asian summer monsoon precipitation. *Geophysical Research Letters* **44**, 3788–95.

* Clement AC and Peterson LC (2008) Mechanisms of abrupt climate change of the last glacial period. *Reviews of Geophysics* **46**, RG4002, 1–39.

* Dahl KA, Broccoli A and Stouffer RJ (2005) Assessing the role of North Atlantic freshwater forcing in millennial scale climate variability: a tropical Atlantic perspective. *Climate Dynamics* **24**, 325–46.

* Dansgaard W (1964) Stable isotopes in precipitation. *Tellus* **16**, 436–68.

* Deplazes G, Luckge A, Stuut J-BW, Patzold J, Kuhlmann H, Husson D, Fant M and Haug GH (2014) Weakening and strengthening of the Indian monsoon during Heinrich events and Dansgaard-Oeschger oscillations. *Paleoceanography* **29**, 99–144.

* Dong JG, Shen C-C, Kong XG, Wang YJ and Duan FC (2018) Asian monsoon dynamics at Dansgaard/Oeschger events 14-8 and Heinrich events 5-4 in northern China. *Quaternary Geochronology* **47**, 72–80.

* Dong XY, Kathayat G, Rasmussen SO, Svensson A, Severinghaus JP, Li HY, Sinha A, Xu Y, Zhang HW, Shi ZG, Cai YJ, Pérez-Mejías C, Baker J, Zhao JY, Spötl C, Columbu A, Ning YF, Strikis NM, Chen ST, Wang XF, Gupta AK, Dutt S, Zhang F, Cruz FW, An ZS, Edwards RL and Cheng H (2022) Coupled atmosphere-ice-ocean dynamics during Heinrich Stadial 2. *Nature Communications* **13**, 5867.

* Draxler RR and Hess GD (1998) An overview of the HYSPLIT_4 modelling system for trajectories, dispersion, and deposition. *Australian Meteorology magazine* **47**, 295–308.

* Dutt S, Gupta AK, Clemens SC, Cheng H, Singh RK, Kathayat G and Edwards RL (2015) Abrupt changes in Indian summer monsoon strength during 33,800 to 5500 years B.P. *Geophysical Research Letters* **42**, 5526–32.

* Dykoski C, Edwards RL, Cheng H, Yuan DX, Cai YJ, Zhang ML, Lin YS, Qing JM, An ZS and Revenaugh J (2005) A high-resolution, absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave, China. *Earth And Planetary Science Letters* **233**, 71–86.

* Gong DY and Ho CH (2002) Shift in the summer rainfall over the Yangtze River valley in the late 1970s. *Geophysical Research Letters* **29**, 71–8.

* Guo ZT, Liu T-S, Guiot J, Wu N, Lü HY and Han JM (1996) High frequency pulses of East Asia monsoon climate in the last two glaciations: link with the North Atlantic. *Climate Dynamics* **12**, 701–9.

* Hiess J, Condon DJ, McLean N, Noble SR (2012) 238U/235U Systematics in terrestrial uranium- bearing minerals. *Science* **335**, 1610–1614.
‘coastal-inland’ pattern, which is probably associated with the decreased summer rainfall amount in inland areas or/and the increased proportion of summer moisture resources from the Pacific Ocean.

**Supplementary material.** To view supplementary material for this article, please visit https://doi.org/10.1017/S0016756824000013

**Acknowledgements.** This research was supported by grants from the National Natural Science Foundation of China (42074071, 42071106, 41874078, 92158208, 42274094), Shenzhen Science and Technology Program (grant no. KQTD20170810111725321), Science and Technology Innovation Committee of Shenzhen Municipality (ZDSYS201802081843490), Technology and Innovation Commission of Shenzhen Municipality (20200925154739001), Southern University of Science and Technology (Y01316111), and Key Special Project for Introduced Talents Team of Southern Marine Science, the International Cooperation Project of Fujian (202210013 to X. J) and Engineering Guangdong Laboratory (Guangzhou) (GML2019ZD0210). U-Th dating was partially supported by grants from the National Science and Technology Council (111-2116-M-002-022-MY3), the National Taiwan University (112L894202 to C.-C.S.), and the Higher Education Sprout Project of the Ministry of Education (110L901001 and 110L8907). Shanghai Sheshan National Geophysical Observatory (SSKP202102).

**Competing interests.** We declare that we have no financial and personal relationships with other people or organizations that can inappropriately influence our work, and there is no professional or other personal interest of any nature or kind in any product, service, and/or company that could be construed as influencing the position presented in, or the review of, the manuscript entitled, ‘Chinese stalagmite δ¹⁸O records reveal the diverse moisture trajectories during the middle to late last glacial period’.

## References

* Andersen KK (2004) High-resolution record of northern hemisphere climate extending into the last interglacial period. *Nature* **431**, 147–51.

* Baker AJ, Sodemann H, Baldini JUL, Breitenbach SFM, Johnson KR, Hunen JV and Zhang PZ (2015) Seasonality of westerly moisture transport in the East Asian summer monsoon and its implications for interpreting precipitation δ¹⁸O. *Journal of Geophysical Research: Atmospheres* **120**, 5850–62.

* Bohm E, Lippold J, Gutjahr M, Frank M, Blaser P, Antz B, Fohlmeister J, Frank N, Andersen MB and Deininger M (2015) Strong and deep Atlantic meridional overturning circulation during the last glacial cycle. *Nature* **517**, 73–6.

* Bradley RM, Eric JS, Christo B, Spruce WS, Cecilia MB, Fudge TJ, Joel BP, Ding QH, Tyler RJ, Jones WCW and Todd S (2017) Global atmospheric teleconnections during Dansgaard-Oeschger events. *Nature Geoscience* **10**, 36–40.

* Broccoli AJ, Dahl KA and Stouffer RJ (2006) Response of the ITCZ to Northern Hemisphere cooling. *Geophysical Research Letters* **33**, L01702.

* Cai YJ, Fung IY, Edwards RL, An ZS, Cheng H, Lee JE, Tan LC, Shen C-C, Wang XF, Day JA, Zhou WJ, Kelly MJ and Chiang JC (2015) Variability of stalagmite-inferred Indian monsoon precipitation over the past 252,000 y. *Proceedings of the National Academy of Sciences of the United States of America* **112**, 2954–9.

* Cao J, Gui S, Su Q and Yang Y (2016) The variability of the Indian-East Asian summer monsoon interface in relation to the spring seesaw mode between the Indian Ocean and the central-western Pacific. *Journal of Climate* **29**, 5027–40.

* Cao J, Hu J and Tao Y (2012) An index for the interface between the Indian summer monsoon and the East Asian summer monsoon. *Journal of Geophysical Research: Atmospheres* **117**, 108.

* Chen FH, Chen XM, Chen JH, Zhou AF, Wu D, Tang LY, Zhang XJ, Huang XZ and Yu JQ (2015) Holocene vegetation history, precipitation changes, and Indian summer monsoon evolution documented from sediments of Xingyun Lake, southwest China. *Journal of Quaternary Science* **29**, 661–74.

* Chen ST, Wang YJ, Cheng H, Edwards RL, Wang XF, Kong XG and Liu DB (2016) Strong coupling of Asian Monsoon and Antarctic climates on suborbital timescales. *Scientific Report* **8**, 32995.

* Cheng H, Edwards RL, Shen C-C, Polyak VJ, Asmerom Y, Woodhead J, Hellstrom J, Wang YJ, Kong XG, Spötl C, Wang XF and Alexander JEC (2013) Improvements in $^{230}$Th dating, $^{230}$Th and $^{234}$U half-life values, and U-Th isotopic measurements by multi-collector inductively coupled plasma mass spectrometry. *Earth and Planetary Science Letters* **371–372**, 82–91.

* Cheng H, Edwards RL, Sinha A, Spötl C, Yi L, Chen ST, Kelly M, Kathayat G, Wang XF, Li XL, Kong XG, Wang YJ, Ning YF and Zhang HW (2016) The Asian monsoon over the past 640,000 years and ice age terminations. *Nature* **534**, 640–6.

* Cheng H, Sinha A, Wang XF, Cruz FW and Edwards RL (2012) The Global Paleomonsoon is seen through speleothem records from Asia and the Americas. *Climate Dynamics* **39**, 1045–62.

* Cheng H, Zhang HW, Zhao JY, Li HY, Ning YF and Kathayat G (2019) Chinese stalagmite paleoclimate researches: a review and perspective. *Science China Earth Sciences* **62**, 1–25.

* Chiang JCH and Bitz CM (2005) Influence of high latitude ice cover on the marine Intertropical Convergence Zone. *Climate Dynamics* **25**, 477–96.

* Chiang JCH, Fung IY, Wu CH, Cai YJ, Edman JP, Liu YW, Day JA, Bhattacharya T, Mondal Y and Labrousse CA (2015) Role of seasonal transitions and westerly jets in East Asian paleoclimate. *Quaternary Science Reviews* **108**, 111–29.

* Chiang JCH, Herman MJ, Yoshimura K and Fung IY (2020) Enriched East Asian oxygen isotope of precipitation indicates reduced summer seasonality in regional climate and westerlies. *Proceedings of the National Academy of Sciences of the United States of America* **117**, 14745–50.

* Chiang JCH, Swenson LM and Kong W (2017) Role of seasonal transitions and the westerlies in the interannual variability of the East Asian summer monsoon precipitation. *Geophysical Research Letters* **44**, 3788–95.

* Clement AC and Peterson LC (2008) Mechanisms of abrupt climate change of the last glacial period. *Reviews of Geophysics* **46**, RG4002, 1–39.

* Dahl KA, Broccoli A and Stouffer RJ (2005) Assessing the role of North Atlantic freshwater forcing in millennial scale climate variability: a tropical Atlantic perspective. *Climate Dynamics* **24**, 325–46.

* Dansgaard W (1964) Stable isotopes in precipitation. *Tellus* **16**, 436–68.

* Deplazes G, Luckge A, Stuut J-BW, Patzold J, Kuhlmann H, Husson D, Fant M and Haug GH (2014) Weakening and strengthening of the Indian monsoon during Heinrich events and Dansgaard-Oeschger oscillations. *Paleoceanography* **29**, 99–144.

* Dong JG, Shen C-C, Kong XG, Wang YJ and Duan FC (2018) Asian monsoon dynamics at Dansgaard/Oeschger events 14-8 and Heinrich events 5-4 in northern China. *Quaternary Geochronology* **47**, 72–80.

* Dong XY, Kathayat G, Rasmussen SO, Svensson A, Severinghaus JP, Li HY, Sinha A, Xu Y, Zhang HW, Shi ZG, Cai YJ, Pérez-Mejías C, Baker J, Zhao JY, Spötl C, Columbu A, Ning YF, Strikis NM, Chen ST, Wang XF, Gupta AK, Dutt S, Zhang F, Cruz FW, An ZS, Edwards RL and Cheng H (2022) Coupled atmosphere-ice-ocean dynamics during Heinrich Stadial 2. *Nature Communications* **13**, 5867.

* Draxler RR and Hess GD (1998) An overview of the HYSPLIT_4 modelling system for trajectories, dispersion, and deposition. *Australian Meteorology magazine* **47**, 295–308.

* Dutt S, Gupta AK, Clemens SC, Cheng H, Singh RK, Kathayat G and Edwards RL (2015) Abrupt changes in Indian summer monsoon strength during 33,800 to 5500 years B.P. *Geophysical Research Letters* **42**, 5526–32.

* Dykoski C, Edwards RL, Cheng H, Yuan DX, Cai YJ, Zhang ML, Lin YS, Qing JM, An ZS and Revenaugh J (2005) A high-resolution, absolute-dated Holocene and deglacial Asian monsoon record from Dongge Cave, China. *Earth And Planetary Science Letters* **233**, 71–86.

* Gong DY and Ho CH (2002) Shift in the summer rainfall over the Yangtze River valley in the late 1970s. *Geophysical Research Letters* **29**, 71–8.

* Guo ZT, Liu T-S, Guiot J, Wu N, Lü HY and Han JM (1996) High frequency pulses of East Asia monsoon climate in the last two glaciations: link with the North Atlantic. *Climate Dynamics* **12**, 701–9.

* Hiess J, Condon DJ, McLean N, Noble SR (2012) 238U/235U Systematics in terrestrial uranium- bearing minerals. *Science* **335**, 1610–1614.
Hu CY, Henderson GM, Huang JH, Xie SC, Sun Y and Johnson KR (2008) Quantification of Holocene Asian monsoon rainfall from spatially separated cave records. *Earth and Planetary Science Letters* **266**, 221–32.

Jacobel AW, McManus JF, Anderson RF and Winckler G (2016) Large deglacial shifts of the Pacific intertropical convergence zone. *Nature Communications* **7**, 10449.

Jaffey AH, Flynn KF, Glendenin LE, Bentley WC, Essling AM (1971) Precision measurement of half-lives and specific activities of 235U and 238U. *Physical Review C* **4**, 1889–1906.

Jiang DB, Lang XM, Tian ZP and Guo DL (2011) Last glacial maximum climate over China from PMIP simulations. *Palaeogeography Palaeoclimatology Palaeoecology* **309**, 347–57.

Kaboth-Bahr S, Bahr A, Yamoah KA, Chuang K-C, Li H-C, Su C-C and Wei K-Y (2021) Rapid humidity changes across the Northern South China Sea during the last ~40 kyrs. *Marine Geology* **440**, 106579.

Kucera M, Weinelt M, Kiefer T, Pflaumann U, Hayes A, Weinelt M, Chen M-T, Mix AC, Borrows TT, Cortijo E, Duprat J, Juggins S and Waelbroeck C (2005) Reconstruction of sea-surface temperatures from assemblages of planktonic foraminifera: multi-technique approach based on geographically constrained calibration data sets and its application to glacial Atlantic and Pacific Oceans. *Quaternary Science Reviews* **24**, 951–98.

Lachniet MS (2015) Are aragonite stalagmites reliable paleoclimate proxies? Tests for oxygen isotope time-series replication and equilibrium. *Geological Society of America Bulletin* **50**, 489–510.

Li D, Tan LC, Cai YJ, Jiang XY, Ma L, Cheng H, Edwards RL, Zhang HW, Gao YL and An ZS (2019) Is Chinese stalagmite δ¹⁸O solely controlled by the Indian summer monsoon? *Climate Dynamics* **53**, 2969–83.

Li YX, Rao ZG, Cao JT, Jiang H and Gao YL (2017) Highly negative oxygen isotopes in precipitation in southwest China and their significance in paleoclimatic studies. *Quaternary International* **440**, 64–71.

Liang YJ, Zhao K, Edwards RL, Wang YJ, Shao QF, Zhang ZQ, Zhao B, Wang Q, Cheng H and Kong XG (2020) East Asian monsoon changes early in the last deglaciation and insights into the interpretation of oxygen isotope changes in the Chinese stalagmite record. *Quaternary Science Reviews* **250**, 106699.

Lisiecki LE and Stern JV (2016) Regional and global benthic δ¹⁸O stacks for the last glacial cycle. *Paleoceanography* **31**, 1368–94.

Liu G, Li X, Chiang HW, Yuan S, Chawchai S, He S, Lu Y, Aung LT, Maung PM, Tun WN, Oo KM and Wang XF (2020) On the glacial-interglacial variability of the Asian monsoon in speleothem δ¹⁸O records. *Science Advances* **6**, eaay8189, 1–10.

Liu JB, Chen JH, Zhang XJ, Li Y, Rao ZG and Chen FH (2015) Holocene East Asian summer monsoon records in northern China and their inconsistency with Chinese stalagmite δ¹⁸O records. *Earth Science Reviews* **148**, 194–208.

Liu ZY, Wen XY, Brady EC, Otto-Bliesner B, Yu G, Lu HY, Cheng H, Wang YJ, Zheng WP, Ding YH, Edwards RL, Cheng J, Liu W and Yang H (2014) Chinese cave records and the East Asia summer monsoon. *Quaternary Science Reviews* **83**, 115–28.

Lynch-Stieglitz J, 2017. The Atlantic meridional overturning circulation and abrupt climate change. *Annual Review of Marine Science* **9**, 83–104.

Maher BA and Thompson R (2012) Oxygen isotopes from Chinese caves: records not of monsoon rainfall but of circulation regime. *Journal of Quaternary Science* **27**, 615–24.

Markle BR, Steig EJ, Buizert C, Schoenemann SW, Bitz CM, Fudge TJ, Pedro JB, Ding QH, Jones TR, White JWC and Sowers T (2016) Global atmospheric teleconnections during Dansgaard-Oeschger events. *Nature Geoscience* **10**, 36–40.

Mohtadi M, Prange M, Oppo DW, Pol-Holz RD, Merkel U, Zhang X, Steinke S and Lucke A (2014) North Atlantic forcing of tropical Indian Ocean climate. *Nature* **509**, 76–80.

Nagashima K, Tada R, Tani A, Sun Y, Isozaki Y, Toyoda S and Hasegawa H (2011) Millennial-scale oscillations of the westerly jet path during the last glacial period. *Journal of Asian Earth Sciences* **40**, 1214–20.

Pausata FSR, Battisti DS, Nisancioglu KH and Bitz CM (2011) Chinese stalagmite δ¹⁸O controlled by changes in the Indian monsoon during a simulated Heinrich event. *Nature Geoscience* **4**, 474–80.

Porter S and An ZS (1995) Correlation between climate events in the North-Atlantic and China during the last glaciation. *Nature* **375**, 305–8.

Qiu H-Y, Li TY, Chen C-J, Huang R, Wang T, Wu Y, Xiao S-Y, Xu Y-Z, Huang Y-Y, Zhang J, Yang Y and Li JY (2021) Significance of active speleothem δ¹⁸O at annual-decadal timescale: a case study from monitoring in Furong Cave. *Applied Geochemistry* **126**, 104873.

Qiu WY, Zhang X, Jiang XY, Hu HM, Ma L, Xiao HY, Cai BG and Shen C-C (2022) Double-plunge structure of the East Asian summer monsoon during Heinrich Stadial 1 recorded in Xianyun Cave, southeastern China. *Quaternary Science Reviews* **282**, 107442.

Sampe T and Xie SP (2010) Large-scale dynamics of the Meiyu-Baiu rainband: environmental forcing by the westerly jet. *Journal of Climate* **23**, 113–34.

Shen C-C, Wu C-C, Cheng H, Edwards RL, Hsieh Y-T, Gallet S, Chang C-C, Li T-Y, Lam DD, Kano A, Hori M and Spötl C (2012) High-precision and high-resolution carbonate ²³⁰Th dating by MC-ICP-MS with SEM protocols. *Geochimica et Cosmochimica Acta* **99**, 71–86.

Siddall M, Rohling EJ, Almogl-Labin A, Hemleben Ch, Meischner D, Schmelzer I and Smeed DA (2003) Sea-level fluctuations during the last glacial cycle. *Nature* **432**, 853–8.

Siddall M, Rohling EJ, Thompson WG and Waelbroeck C (2008) Maine isotope stage 3 sea level fluctuations: data synthesis and new outlook. *Reviews of Geophysics* **46**, RG4003.

Sun YB, Chen J, Clemens SC, Liu QS, Ji JF and Tada R (2006) East Asian monsoon variability over the last seven glacial cycles recorded by a loess sequence from the northwestern Chinese Loess Plateau. *Geochemistry Geophysics and Geosystems* **7**, 1–16.

Sun YB, Clemens SC, Morrill C, Lin XP, Wang XL and An ZS (2012). Influence of Atlantic meridional overturning circulation on the East Asian winter monsoon. *Nature Geoscience* **5**, 46–9.

Tan M (2009) Circulation effect: climatic significance of the short term variability of the oxygen isotopes stalagmites form monsoonal China—dialogue between paleoclimate records and modern climate research. *Quaternary Science* **29**, 851–62 (in Chinese with English abstract).

Tan M (2011) Trade-wind driven inverse coupling between stalagmite δ¹⁸O from monsoon region of China and large scale temperature. *Quaternary Science* **31**, 1086–97 (in Chinese with English abstract).

Tan M (2013) Circulation effect: response of precipitation δ¹⁸O to the ENSO cycle in monsoon regions of China. *Climate Dynamics* **42**, 1067–77.

Tan M (2016) Circulation background of climate patterns in the past millennium: uncertainty analysis and re-reconstruction of ENSO-like state. *Science China Earth Science* **59**, 1225–41.

Tierney JE, Pausata FSR and DeMenocal P (2016) Deglacial Indian monsoon failure and North Atlantic Stadials linked by Indian Ocean surface cooling. *Nature Geoscience* **9**, 46–50.

Wang HJ and Chen HP (2012) Climate control for southeastern China moisture and precipitation: Indian or East Asian monsoon? *Journal of Geophysical Research: Atmospheres* **117**, 1–9.

Wang YJ, Cheng H, Edwards RL, An ZS, Wu JY, Shen C-C and Dorale JA (2001) A high-resolution absolute-dated Late Pleistocene monsoon record from Hulu Cave, China. *Science* **294**, 2345–8.

Yang XL, Liu JB, Liang FY, Yuan DX, Yang Y, Lu YB and Chen FH (2014) Holocene stalagmite δ¹⁸O records in the East Asian monsoon region and their correlation with those in the Indian monsoon region. *Holocene* **24**, 1657–64.

Yuan DX, Cheng H, Edwards RL, Dykoski CA, Kelly MJ, Zhang ML, Qing JM, Lin YS, Wang YJ, Wu JY, Dorale JA, An ZS and Cai YJ (2004) Timing, duration, and transitions of the last interglacial Asian monsoon. *Science* **304**, 575–8.

Zhang HW, Cai YJ, Tan LC, Qin SJ and An ZS (2014) Stable isotope composition alteration produced by the aragonite-to-calcite transformation in speleothems and implications for paleoclimate reconstructions. *Sedimentary Geology* **309**, 1–14.

Zhang R and Delworth TL (2005) Simulated tropical response to a substantial weakening of the Atlantic thermohaline circulation. *Journal of Climate* **18**, 1853–60.

Zhang X, Qiu WY, Jiang XY, Hu HM, Xiao HY, Cai BG and Shen C-C (2021) Three-phase structure of the East Asia summer monsoon during Heinrich
Hu CY, Henderson GM, Huang JH, Xie SC, Sun Y and Johnson KR (2008) Quantification of Holocene Asian monsoon rainfall from spatially separated cave records. *Earth and Planetary Science Letters* **266**, 221–32.

Jacobel AW, McManus JF, Anderson RF and Winckler G (2016) Large deglacial shifts of the Pacific intertropical convergence zone. *Nature Communications* **7**, 10449.

Jaffey AH, Flynn KF, Glendenin LE, Bentley WC, Essling AM (1971) Precision measurement of half-lives and specific activities of 235U and 238U. *Physical Review C* **4**, 1889-1906.

Jiang DB, Lang XM, Tian ZP and Guo DL (2011) Last glacial maximum climate over China from PMIP simulations. *Palaeogeography Palaeoclimatology Palaeocology* **309**, 347–57.

Kaboth-Bahr S, Bahr A, Yamoah KA, Chuang K-C, Li H-C, Su C-C and Wei K-Y (2021) Rapid humidity changes across the Northern South China Sea during the last ~40 kyrs. *Marine Geology* **440**, 106579.

Kucera M, Weinelt M, Kiefer T, Pflaumann U, Hayes A, Weinelt M, Chen M-T, Mix AC, Borrows TT, Cortijo E, Duprat J, Juggins S and Waelbroeck C (2005) Reconstruction of sea-surface temperatures from assemblages of planktonic foraminifera: multi-technique approach based on geographically constrained calibration data sets and its application to glacial Atlantic and Pacific Oceans. *Quaternary Science Reviews* **24**, 951–98.

Lachniet MS (2015) Are aragonite stalagmites reliable paleoclimate proxies? Tests for oxygen isotope time-series replication and equilibrium. *Geological Society of America Bulletin* **50**, 489–510.

Li D, Tan LC, Cai YJ, Jiang XY, Ma L, Cheng H, Edwards RL, Zhang HW, Gao YL and An ZS (2019) Is Chinese stalagmite δ¹⁸O solely controlled by the Indian summer monsoon? *Climate Dynamics* **53**, 2969–83.

Li YX, Rao ZG, Cao JT, Jiang H and Gao YL (2017) Highly negative oxygen isotopes in precipitation in southwest China and their significance in paleoclimatic studies. *Quaternary International* **440**, 64–71.

Liang YJ, Zhao K, Edwards RL, Wang YJ, Shao QF, Zhang ZQ, Zhao B, Wang Q, Cheng H and Kong XG (2020) East Asian monsoon changes early in the last deglaciation and insights into the interpretation of oxygen isotope changes in the Chinese stalagmite record. *Quaternary Science Reviews* **250**, 106699.

Lisiecki LE and Stern JV (2016) Regional and global benthic δ¹⁸O stacks for the last glacial cycle. *Paleoceanography* **31**, 1368–94.

Liu G, Li X, Chiang HW, Yuan S, Chawchai S, He S, Lu Y, Aung LT, Maung PM, Tun WN, Oo KM and Wang XF (2020) On the glacial-interglacial variability of the Asian monsoon in speleothem δ¹⁸O records. *Science Advances* **6**, eaay8189, 1–10.

Liu JB, Chen JH, Zhang XJ, Li Y, Rao ZG and Chen FH (2015) Holocene East Asian summer monsoon records in northern China and their inconsistency with Chinese stalagmite δ¹⁸O records. *Earth Science Reviews* **148**, 194–208.

Liu ZY, Wen XY, Brady EC, Otto-Bliesner B, Yu G, Lu HY, Cheng H, Wang YJ, Zheng WP, Ding YH, Edwards RL, Cheng J, Liu W and Yang H (2014) Chinese cave records and the East Asia summer monsoon. *Quaternary Science Reviews* **83**, 115–28.

Lynch-Stieglitz J, 2017. The Atlantic meridional overturning circulation and abrupt climate change. *Annual Review of Marine Science* **9**, 83–104.

Maher BA and Thompson R (2012) Oxygen isotopes from Chinese caves: records not of monsoon rainfall but of circulation regime. *Journal of Quaternary Science* **27**, 615–24.

Markle BR, Steig EJ, Buizert C, Schoenemann SW, Bitz CM, Fudge TJ, Pedro JB, Ding QH, Jones TR, White JWC and Sowers T (2016) Global atmospheric teleconnections during Dansgaard-Oeschger events. *Nature Geoscience* **10**, 36–40.

Mohtadi M, Prange M, Oppo DW, Pol-Holz RD, Merkel U, Zhang X, Steinke S and Lucke A (2014) North Atlantic forcing of tropical Indian Ocean climate. *Nature* **509**, 76–80.

Nagashima K, Tada R, Tani A, Sun Y, Isozaki Y, Toyoda S and Hasegawa H (2011) Millennial-scale oscillations of the westerly jet path during the last glacial period. *Journal of Asian Earth Sciences* **40**, 1214–20.

Pausata FSR, Battisti DS, Nisancioglu KH and Bitz CM (2011) Chinese stalagmite δ¹⁸O controlled by changes in the Indian monsoon during a simulated Heinrich event. *Nature Geoscience* **4**, 474–80.

Porter S and An ZS (1995) Correlation between climate events in the North-Atlantic and China during the last glaciation. *Nature* **375**, 305–8.

Qiu H-Y, Li TY, Chen C-J, Huang R, Wang T, Wu Y, Xiao S-Y, Xu Y-Z, Huang Y-Y, Zhang J, Yang Y and Li JY (2021) Significance of active speleothem δ¹⁸O at annual-decadal timescale: a case study from monitoring in Furong Cave. *Applied Geochemistry* **126**, 104873.

Qiu WY, Zhang X, Jiang XY, Hu HM, Ma L, Xiao HY, Cai BG and Shen C-C (2022) Double-plunge structure of the East Asian summer monsoon during Heinrich Stadial 1 recorded in Xianyun Cave, southeastern China. *Quaternary Science Reviews* **282**, 107442.

Sampe T and Xie SP (2010) Large-scale dynamics of the Meiyu-Baiu rainband: environmental forcing by the westerly jet. *Journal of Climate* **23**, 113–34.

Shen C-C, Wu C-C, Cheng H, Edwards RL, Hsieh Y-T, Gallet S, Chang C-C, Li T-Y, Lam DD, Kano A, Hori M and Spötl C (2012) High-precision and high-resolution carbonate ²³⁰Th dating by MC-ICP-MS with SEM protocols. *Geochimica et Cosmochimica Acta* **99**, 71–86.

Siddall M, Rohling EJ, Almogl-Labin A, Hemleben Ch, Meischner D, Schmelzer I and Smeed DA (2003) Sea-level fluctuations during the last glacial cycle. *Nature* **432**, 853–8.

Siddall M, Rohling EJ, Thompson WG and Waelbroeck C (2008) Maine isotope stage 3 sea level fluctuations: data synthesis and new outlook. *Reviews of Geophysics* **46**, RG4003.

Sun YB, Chen J, Clemens SC, Liu QS, Ji JF and Tada R (2006) East Asian monsoon variability over the last seven glacial cycles recorded by a loess sequence from the northwestern Chinese Loess Plateau. *Geochemistry Geophysics and Geosystems* **7**, 1–16.

Sun YB, Clemens SC, Morrill C, Lin XP, Wang XL and An ZS (2012). Influence of Atlantic meridional overturning circulation on the East Asian winter monsoon. *Nature Geoscience* **5**, 46–9.

Tan M (2009) Circulation effect: climatic significance of the short term variability of the oxygen isotopes stalagmites form monsoonal China—dialogue between paleoclimate records and modern climate research. *Quaternary Science* **29**, 851–62 (in Chinese with English abstract).

Tan M (2011) Trade-wind driven inverse coupling between stalagmite δ¹⁸O from monsoon region of China and large scale temperature. *Quaternary Science* **31**, 1086–97 (in Chinese with English abstract).

Tan M (2013) Circulation effect: response of precipitation δ¹⁸O to the ENSO cycle in monsoon regions of China. *Climate Dynamics* **42**, 1067–77.

Tan M (2016) Circulation background of climate patterns in the past millennium: uncertainty analysis and re-reconstruction of ENSO-like state. *Science China Earth Science* **59**, 1225–41.

Tierney JE, Pausata FSR and DeMenocal P (2016) Deglacial Indian monsoon failure and North Atlantic Stadials linked by Indian Ocean surface cooling. *Nature Geoscience* **9**, 46–50.

Wang HJ and Chen HP (2012) Climate control for southeastern China moisture and precipitation: Indian or East Asian monsoon? *Journal of Geophysical Research: Atmospheres* **117**, 1–9.

Wang YJ, Cheng H, Edwards RL, An ZS, Wu JY, Shen C-C and Dorale JA (2001) A high-resolution absolute-dated Late Pleistocene monsoon record from Hulu Cave, China. *Science* **294**, 2345–8.

Yang XL, Liu JB, Liang FY, Yuan DX, Yang Y, Lu YB and Chen FH (2014) Holocene stalagmite δ¹⁸O records in the East Asian monsoon region and their correlation with those in the Indian monsoon region. *Holocene* **24**, 1657–64.

Yuan DX, Cheng H, Edwards RL, Dykoski CA, Kelly MJ, Zhang ML, Qing JM, Lin YS, Wang YJ, Wu JY, Dorale JA, An ZS and Cai YJ (2004) Timing, duration, and transitions of the last interglacial Asian monsoon. *Science* **304**, 575–8.

Zhang HW, Cai YJ, Tan LC, Qin SJ and An ZS (2014) Stable isotope composition alteration produced by the aragonite-to-calcite transformation in speleothems and implications for paleoclimate reconstructions. *Sedimentary Geology* **309**, 1–14.

Zhang R and Delworth TL (2005) Simulated tropical response to a substantial weakening of the Atlantic thermohaline circulation. *Journal of Climate* **18**, 1853–60.

Zhang X, Qiu WY, Jiang XY, Hu HM, Xiao HY, Cai BG and Shen C-C (2021) Three-phase structure of the East Asia summer monsoon during Heinrich
Stadial 4 recorded in Xianyun Cave, southeastern China. *Quaternary Science Reviews* **274**, 107267.

Zhao JY, Cheng H, Yang Y, Tan LC, Spötl C, Ning YF, Zhang HW, Cheng X, Sun Z, Li XL, Li HY, Liu W and Edwards RL (2019) Reconstructing the Western boundary variability of the Western Pacific subtropical high over the past 200 years via Chinese cave oxygen isotope records. *Climate Dynamics* **52**, 3741–57.

Zhao SH, Liu ZF, Colin C, Zhao YL, Wang XX and Jian ZM (2018) Responses of the East Asian summer monsoon in the low-latitude South China Sea to high-latitude millennial-scale climatic changes during the last glaciation: evidence from a high-resolution clay mineralogical record. *Paleoceanography and Paleoclimatology* **33**, 2017PA003235.

Zhou HY, Zhao JX, Feng YX, Gagan MK, Zhou GQ and Yan J (2008) Distinct climate change synchronous with Heinrich event one, recorded by stable oxygen and carbon isotopic compositions in stalagmites from China. *Quaternary Research* **69**, 306–15.

Zhou TJ, Yu R, Zhang J, Drange H, Cassou C, Deser C, Hodson DLR, Sanchez-Gomez E, Li J, Keenlyside N, Xin XG and Okumura Y (2009) Why the Western Pacific subtropical high has extended westward since the late 1970s. *Journal of Climate* **22**, 2199–215.

Zorzi C, Goñi MFS, Anupama K, Prasad S, Hanquiez V, Johnson J and Giosan L (2015) Indian monsoon variations during three contrasting climatic periods: the Holocene, Heinrich Stadial 2, and the last interglacial-glacial transition. *Quaternary Science Reviews* **125**, 50–60.
Stadial 4 recorded in Xianyun Cave, southeastern China. *Quaternary Science Reviews* **274**, 107267.

Zhao JY, Cheng H, Yang Y, Tan LC, Spötl C, Ning YF, Zhang HW, Cheng X, Sun Z, Li XL, Li HY, Liu W and Edwards RL (2019) Reconstructing the Western boundary variability of the Western Pacific subtropical high over the past 200 years via Chinese cave oxygen isotope records. *Climate Dynamics* **52**, 3741–57.

Zhao SH, Liu ZF, Colin C, Zhao YL, Wang XX and Jian ZM (2018) Responses of the East Asian summer monsoon in the low-latitude South China Sea to high-latitude millennial-scale climatic changes during the last glaciation: evidence from a high-resolution clay mineralogical record. *Paleoceanography and Paleoclimatology* **33**, 2017PA003235.

Zhou HY, Zhao JX, Feng YX, Gagan MK, Zhou GQ and Yan J (2008) Distinct climate change synchronous with Heinrich event one, recorded by stable oxygen and carbon isotopic compositions in stalagmites from China. *Quaternary Research* **69**, 306–15.

Zhou TJ, Yu R, Zhang J, Drange H, Cassou C, Deser C, Hodson DLR, Sanchez-Gomez E, Li J, Keenlyside N, Xin XG and Okumura Y (2009) Why the Western Pacific subtropical high has extended westward since the late 1970s. *Journal of Climate* **22**, 2199–215.

Zorzi C, Goñi MFS, Anupama K, Prasad S, Hanquiez V, Johnson J and Giosan L (2015) Indian monsoon variations during three contrasting climatic periods: the Holocene, Heinrich Stadial 2, and the last interglacial-glacial transition. *Quaternary Science Reviews* **125**, 50–60.