Sea Surface Temperature Variability Research Articles

Interhemispheric synchrony of Mid-Late Holocene SST in the Indo-Pacific warm pool linked to upper water column dynamics

The sea surface temperature (SST) of the Indo-Pacific warm pool (IPWP) plays a crucial role in global climate system. Despite its importance, knowledge of its SST changes over century to millennial timescales remains limited and controversial, hindering our comprehensive understanding of climate change mechanisms. Our study bridges this gap with a novel SST reconstruction spanning approximately 7000 years, derived by aggregating 75 coral Sr/Ca-SST records from the South China Sea. These records reveal a notable synchronization in SST variations with foraminiferal reconstructions across the IPWP, reflecting a prevalent pattern of interhemispheric cooling and multi-century variability in the IPWP during the Mid-Late Holocene. Further analysis reveals that these (sub)millennial SST variations are closely linked to subsurface water temperature and ocean heat content (OHC) changes, supporting the theory that the release of deeper OHC acted as a heat source for low-latitude processes during the Holocene.

Seasonal evolution modes of tropical sea surface temperature anomalies and their links to antarctic sea ice anomalies

Abstract Previous studies have explored the teleconnections between variability of Antarctic sea ice cover and tropical sea surface temperature (SST) across the Pacific, Atlantic, and Indian Ocean basins, typically focusing on each basin individually. However, there has been limited investigation into the impact of tropical SSTs—particularly from a seasonal evolution perspective—on Antarctic sea ice cover. In this study, we employ the self-organizing map method to identify and analyze the primary modes of seasonal SST evolution in the tropical oceans from 1854 to 2022. We also project changes in the frequency of these modes through the 21st century. Moreover, we examine the seasonal variability of Antarctic sea ice concentration in relation to these tropical SST modes over the past four decades. Our results reveal that tropical SST anomalies display both uniform and shifted seasonal evolution patterns. Notably, the frequency of switched modes—namely, transitions from La Niña to El Niño (node 8) and from El Niño to La Niña (node 3)—is expected to increase in future climate. Interestingly, nearly mirrored SST seasonal evolution patterns do not lead to entirely opposite atmospheric circulation anomalies in the southern mid-high latitudes, nor do they result in completely inverse Antarctic sea ice cover anomalies.

Projected increase in the frequency of extremely active Atlantic hurricane seasons.

Future changes to the year-to-year swings between active and inactive North Atlantic tropical cyclone (TC) seasons have received little attention, yet may have great societal implications in areas prone to hurricane landfalls. This work investigates past and future changes in North Atlantic TC activity, focusing on interannual variability and evaluating the contributions from anthropogenic forcing. We show that interannual variability of Atlantic TC activity has already increased, evidenced by an increase in the occurrence of both extremely active and inactive TC seasons. TC-resolving general circulation models project a 36% increase in the variance of North Atlantic TC activity, measured by accumulated cyclone energy, by the middle of the 21st century. These changes are the result of increased variability in vertical wind shear and atmospheric stability, in response to enhanced Pacific-to-Atlantic interbasin sea surface temperature variations. Robust anthropogenic-forced intensification in the variability of Atlantic TC activity will continue in the future, with important implications for emergency planning and societal preparedness.

Role of the Europe–China Pattern Teleconnection in the Interdecadal Autumn Dry–Wet Fluctuations in Central China

Based on statistical analyses of long-term reanalysis data, we have investigated the interdecadal variations of autumn precipitation in central China (APC-d) and the associated atmospheric teleconnection. It reveals that the increased autumn rainfall in central China during the last decade is a portion of the APC-d, which exhibits a high correlation coefficient of 0.7 with the interdecadal variations of the Europe–China pattern (EC-d pattern) teleconnection. The EC-d pattern teleconnection presents in a “+-+” structure over Eurasia, putting central China into the periphery of a quasi-barotropic anticyclonic high-pressure anomaly. Driven by positive vorticity advection and the inflow of warmer and moist air from the south, central China experiences enhanced ascending motion and abundant water vapor supply, resulting in increased rainfall. Further analysis suggests that the EC-d pattern originates from the exit of the North Atlantic jet and propagates eastward. It is captured by the Asian westerly jet stream and proceeds towards East Asia through the wave–mean flow interaction. The wave train acquires effective potential energy from the mean flow by the baroclinic energy conversion and simultaneously obtains kinetic energy from the basic westerly jet zones across the North Atlantic and the East Asian coasts. The interdecadal variation of the mid-latitude North Atlantic sea surface temperature (MAT-d) exhibits a significant negative relationship with EC-d, serving as a modulating factor for the EC-d pattern teleconnection. Experiments with CMIP6 models predict that the interdecadal variations in APC-d, EC-d, and MAT-d will maintain stable high correlations for the rest of the 21st century. These findings may contribute to forecasting the interdecadal autumn dry–wet conditions in central China.

Dakar Niño under global warming investigated by a high-resolution regionally coupled model

Abstract. In this study, we investigated interannual variability in sea surface temperature (SST) along the northwestern African coast, focusing on strong Dakar Niño and Niña events and their potential alterations under the RCP8.5 emission scenario for global warming, using a high-resolution regional coupled model. Our model accurately reproduces the SST seasonal cycle along the northwestern African coast, including its interannual variability in terms of amplitude, timing, and the position of maximum variability. Comparing Dakar Niño variability between the 1980–2010 and 2069–2099 periods, we found that it intensifies under a warmer climate without changing its location and timing. The intensification is more pronounced during Dakar Niñas (cold SST events) than during Dakar Niños (warm SST events). In the future, SST variability will be correlated with ocean temperature and vertical motion at deeper layers. The increase in Dakar Niño variability can be explained by the larger variability in meridional wind stresses, which is likely to be amplified in the future by enhanced land–sea thermal contrast and associated sea-level-pressure anomalies extending from the Iberian Mediterranean area. A heat budget analysis of the mixed layer suggests that surface heat flux and horizontal-advection anomalies are comparably important for Dakar Niño and Niña events in the present climate. However, the future intensification of Dakar Niños and Niñas is likely to be driven by surface heat flux (latent heat flux and shortwave radiation). While horizontal- and vertical-advection anomalies also contribute to the intensification, their roles are secondary.

How do model biases affect large-scale teleconnections that control Southwest US precipitation? Part II: seasonal models

Abstract We explore the skill in predicting Southwest United States (SWUS) October to March precipitation and associated large-scale teleconnections in an ensemble of hindcasts from seasonal prediction systems. We identify key model biases that degrade the models’ capability to predict SWUS precipitation. The subtropical jet in the Pacific sector is generally too zonal and elongated. This is reflected in the models’ North Pacific ENSO teleconnections that are generally too weak with exaggerated northwest-southeast tilt, compared to observations. Also, the models are too dependent on tropical, El Niño-like, wave train anomalies for producing high seasonal SWUS precipitation, when in observations there is a larger influence of zonal Rossby wave trains such as the one observed in 2016/17. Overall, this is consistent with biases in the basic flow inducing errors in the propagation of zonal wave trains in the North Pacific, which affects SWUS precipitation downstream. Although higher skill may be gained from reducing mean flow biases in the models, a case study of the 2016/17 winter illustrates the great challenge behind skillful seasonal prediction of SWUS precipitation. Unsurprisingly, the almost record-breaking precipitation observed that year in the absence of ENSO is not predicted in the hindcasts, and model perturbation experiments suggest that even a perfect prediction of tropical sea surface temperature and tropical atmospheric variability would not have sufficed to produce a reasonable seasonal precipitation prediction. On a more positive note, our perturbation experiments suggest a potential role for Arctic variability that supports findings from prior studies and suggests re-examining high-latitude drivers of SWUS precipitation.

Widening of wind-stress anomalies amplifies ENSO in a warming climate

Abstract Climate change simulations generally indicate the strengthening of El Niño Southern Oscillation (ENSO) sea surface temperature (SST) variability through the 21st century, yet a robust physical mechanism explaining this change across different models is still lacking, and the projections for ENSO amplitude exhibit a large spread. Most commonly, changes in the background state of the tropical Pacific are invoked to explain these changes of ENSO. Here we show that changes in the structure of wind-stress anomalies associated with ENSO are potentially as important as these background state changes. Specifically, changes in the magnitude, meridional width, and zonal structure of wind-stress anomalies can explain approximately 53% of the inter-model variance in the projected change of ENSO magnitude through the 21st century as well as 43% in ENSO periodicity changes. Among these changes in the wind structure, the meridional widening of wind anomalies plays the most important role. To demonstrate that these changes are indeed critical, we develop a hybrid model of ENSO based on the Community Earth System Model version 2, which incorporates a dynamical ocean coupled to a simplified statistical atmosphere within the tropical Pacific. In the absence of external forcing and corresponding mean-state changes, the imposed changes in wind stress anomalies in this hybrid model result in an increase of ENSO amplitudes of nearly 10% along the equator. Our results are also theoretically supported by a recharge-oscillator model that incorporates the meridional wind structure. Thus, changes in the structure of wind-stress anomalies, together with changes in the mean state, likely play a critical role in the projected strengthening of ENSO.

Comparative Evaluation of Niño1+2 and Niño3.4 Indices in Terms of ENSO Effects Over the Euro‐Mediterranean Region

ABSTRACTGlobal or regional impacts of El Niño Southern Oscillation (ENSO) have predominantly been investigated through the Niño3.4 index, representing the sea surface temperature (SST) variations in the central Tropical Pacific. In this study, we comparatively evaluated the usefulness of Niño1+2, a relatively less utilised index that represents SST variability in the Eastern Tropical Pacific. In our analyses, we focused on ENSO impacts on Euro‐Mediterranean (EM) climate variability during boreal winter, using data from the NCEP/NCAR Reanalysis. The correlation analysis involving Niño1+2 depicts more distinct temperature and precipitation patterns over the EM region. Amongst the SST‐based Niño indices, it has the highest correlation with the East Atlantic index (0.47, statistically significant at > 99% confidence level), a prominent regional teleconnection associated primarily with the strength of the East Atlantic ridge, which produces dipole‐type climate patterns between East Atlantic/Western Europe and Central/Eastern Mediterranean. Moreover, its lagged correlations with the following spring (0.39), summer (0.31), and autumn (0.36) are all statistically significant at ≥ 99% confidence levels. The composite analysis shows that different Niño regions have distinct effects on atmospheric circulation and climate in the EM region. The Niño1+2 index is particularly helpful in identifying the years when warm SST anomalies of El Niño extend to the Eastern Equatorial Pacific, which results in a reversal of temperatures across the EM region. Thus, this study suggests that Niño1+2 is a useful index for studying climate variability and predictability in the EM region, especially when used in conjunction with other Niño indices, as it captures some ENSO features that they may not encompass.

Reconstructing tropical monthly sea surface temperature variability by assimilating coral proxy datasets

Coral reconstruction often serves as a major proxy of high-resolution sea surface temperature (SST) variability beyond the instrumental era. However, coral reconstructions are sparse and are usually studied for interannual variability, with few studies on the monthly features. In this study, we reconstruct the monthly SST spatial field by applying the paleoclimate data assimilation method to the coral records of the latest CoralHydro2k data set for the instrument period of 1880–2000. A comparison with observed SST variability shows that our assimilated tropical SST variability performs reasonably well for the seasonal cycle and monthly ENSO characteristics, notably the phase-locking and onset timing, and more realistic spatial fields relative to the model simulations. This study suggests the feasibility of applying paleoclimate data assimilation to reconstruct the monthly SST in the historical period.

A coordination attention residual U-Net model for enhanced short and mid-term sea surface temperature prediction

Sea surface temperature (SST) is crucial for studying global oceans and evaluating ecosystems. Accurately predicting short and mid-term daily SST has been a significant challenge in oceanography. Traditional deep learning methods can handle temporal data and spatial features but often struggle with long-range spatiotemporal dependencies. To address this, we propose a coordination attention residual U-Net(CResU-Net) model designed to better capture the dynamic spatiotemporal correlations of high-resolution SST. The model integrates coordinate attention mechanisms, multiple residual modules, and depthwise separable convolutions to enhance prediction capabilities. The spatiotemporal variations of SST across different areas of the South China Sea are complex, making accurate predictions challenging. Experiments across various regions of the South China Sea show the model’s effectiveness and robust generalization in predicting high-resolution daily SST. For a 10-day forecast period, the model achieves approximately 0.3 °C in RMSE, outperforming several advanced models.

Persistently active El Niño–Southern Oscillation since the Mesozoic

The El Niño-Southern Oscillation (ENSO), originating in the central and eastern equatorial Pacific, is a defining mode of interannual climate variability with profound impact on global climate and ecosystems. However, an understanding of how the ENSO might have evolved over geological timescales is still lacking, despite a well-accepted recognition that such an understanding has direct implications for constraining human-induced future ENSO changes. Here, using climate simulations, we show that ENSO has been a leading mode of tropical sea surface temperature (SST) variability in the past 250 My but with substantial variations in amplitude across geological periods. We show this result by performing and analyzing a series of coupled time-slice climate simulations forced by paleogeography, atmospheric CO2 concentrations, and solar radiation for the past 250 My, in 10-My intervals. The variations in ENSO amplitude across geological periods are little related to mean equatorial zonal SST gradient or global mean surface temperature of the respective periods but are primarily determined by interperiod difference in the background thermocline depth, according to a linear stability analysis. In addition, variations in atmospheric noise serve as an independent contributing factor to ENSO variations across intergeological periods. The two factors together explain about 76% of the interperiod variations in ENSO amplitude over the past 250 My. Our findings support the importance of changing ocean vertical thermal structure and atmospheric noise in influencing projected future ENSO change and its uncertainty.

Hydroclimate modulation of central-eastern Mexico by the North Atlantic subtropical high since the little ice age

We reconstructed the hydroclimate of central-eastern Mexico over the last 700 common era (CE) based on inferences from multi-proxies from a stalagmite (K-Inc) collected at Karmidas cave, eastern México. Projections on hydroclimate variability in Mexico raise concerns about possible future occurrences of severe droughts and seasonal water balance fluctuations related to increased global temperatures caused by anthropogenic climate change (Murray-Tortarolo, 2021). The eastern region influences the production and supply of food to Mexico. Simulations of past climates, validated by paleoclimate records, yield valuable perspectives on climate change and enhance our understanding of future projections. However, the paucity of paleoclimatic records hinders understanding past hydroclimatic variations and their climatic mechanisms in eastern Mexico. Our record covers the Little Ice Age (LIA) through the Historical Interval (HI), a crucial period for understanding the climate repercussions spanning the transition from Earth's climatic history to the post-industrial era. The reduced intensity of the North Atlantic Subtropical High (NASH) during the LIA enabled a prominent negative phase of NAO-like variability from 1600 CE until the end of LIA. Consequently, preferent meridional airflow within the continent fosters the encounter of moisture-laden intrusions with the increased frequency of cold surges as the occurrence of frontal rain in eastern Mexico, impairing the amount effect on the K-Inc δ18O record. However, after the artificial opening of the cave in 1910 CE, the δ18O records of K-Inc began to exhibit a ∼20-year oscillatory periodicity. In this context, the trace elements of K-Inc help elucidate the climatic conditions that governed the precipitation regime during the investigated period. The visual alignment between the zonal sea surface temperature (SST) variability in the eastern equatorial Pacific (EEP) and the trace elements (Sr/Ca and Ba/Ca) of K-Inc reveals their relationship. Warm zonal SST in the EEP appears to be associated with changes in the length of the winter and summer seasons in eastern Mexico during the LIA. In contrast, over the HI, the trace elements of K-Inc show an anti-phase response to Warm zonal SST in the EEP, denoting wetter climate conditions at the vicinities of Karmidas Cave. This configuration led to questioning the influence of SST zonal variability in the EEP during the HI, which was probably masked by more relevant climate forcing. Our findings enabled us to draw climate scenarios by addressing the main climate drivers in our records.

A global poleward shift of atmospheric rivers.

Atmospheric rivers (ARs) are key agents in distributing extratropical precipitation and transporting moisture poleward. Climate models forced by historical anthropogenic forcing suggest an increase in AR activity in the extratropics over the past four decades. However, reanalyses indicate a ~6° to 10° poleward shift of ARs during boreal winter in both hemispheres, featuring a rise along 50°N and 50°S and a decrease along 30°N and 30°S. Our analysis demonstrates that low-frequency sea surface temperature variability in the tropical eastern Pacific exhibits a cooling tendency since 2000 that plays a key role in driving this global AR shift, mostly over extratropical oceans, through a tropical-driven eddy-mean flow feedback. This mechanism also operates on interannual timescales, controlled by the El Niño-Southern Oscillation, and is less pronounced over the Southern Ocean due to weaker eddy activity during austral summer. These highlight the sensitivity of ARs to large-scale circulation changes driven by both internal variability and external forcing in current and upcoming decades.

Planktonic foraminiferal assemblages as tracers of paleoceanographic changes within the northern Benguela current system since the Early Pleistocene

Abstract. The Benguela Upwelling System (BUS), located in the southeastern Atlantic Ocean, represents one of the world's most productive regions. This system is delimited to the south by the Agulhas retroflection region. The northern boundary of the BUS is, instead, represented by the Angola–Benguela Front (ABF), which is a thermal feature separating warm waters of the Angola Basin (including the South Atlantic Central Water; SACW) from the cooler Benguela Oceanic Current (BOC). We performed statistical analyses on planktonic foraminiferal assemblages in 94 samples from Holes U1575A and U1576A, cored during International Ocean Discovery Program (IODP) Expedition 391. Drilled sites are located along the Tristan–Gough–Walvis Ridge (TGW) seamount track in the northern sector of the BUS (offshore the Namibian continental margin). The analyzed stratigraphic intervals span the Early–Late Pleistocene, marked by the Early–Middle Pleistocene transition (EMPT; 1.40–0.40 Myr), during which important glacial–interglacial sea surface temperature (SST) variabilities occurred. This work provides novel insights on the local paleoceanographic evolution of the northern BUS and associated thermocline variability based on the ecological significance of the foraminiferal assemblages. Specifically, variations in the assemblage content allowed us to characterize the different water masses (BOC, SACW, and Agulhas waters) and reconstruct their interactions during the Quaternary. The interplay of the previously mentioned water masses induced perturbations in the BUS (ABF latitudinal shifts and input of tropical waters from the Agulhas retroflection region). Furthermore, we investigated the possible link between changes in the paleoceanographic conditions and climatic events (e.g., Benguela Niño-/Niña-like phases and deglaciation stages) recorded since the EMPT.

Two-decadal variation of satellite-derived chlorophyll-a concentration and sea surface temperature in the Arabian Sea and Bay of Bengal

Two-decadal variation of satellite-derived chlorophyll-a concentration and sea surface temperature in the Arabian Sea and Bay of Bengal

Influence of mesoscale sea‐surface temperature structures on the Mediterranean cyclone Ianos in convection‐permitting simulations: Contributions of surface warming and cold wakes

AbstractThis study examined the influence of mesoscale sea‐surface temperature (SST) structures and warming on the tropical‐like cyclone Ianos. The controlling mechanisms were investigated using a set of sensitivity experiments with the non‐hydrostatic weather research and forecasting model. A simulation forced by a high‐resolution SST field from the mesoscale eddy‐resolving HYCOM reanalysis (CTRLH) performed better than a simulation forced by a coarse‐resolution SST field from Optimum Interpolation Sea Surface Temperature (OISST) in terms of precipitation, wind intensity, and landfall location. The suppression of mesoscale SST variability affected the cyclone's intensity, its path, and associated processes. This effect depends on the nature of the anomalies (cold/warm). The suppression of mesoscale SST forcing intensified wind and precipitation in the presence of cold anomalies, leading to a northwestward shift of the storm's path during its mature phase. The magnitude of these changes varied with the scale of smoothing (mesoscale SST). Compared to the CTRLH simulation, a simulation in which a 60‐km smoothing filter was applied to the SST field resulted in a cyclone centered slightly to the west (or delayed by 3–6 hr), which was further shifted to the west (or delayed) when the smoothing filter was increased to 150 km. The intensification of wind and precipitation may be associated with the reduction in the storm's cold wake, which would lead to enhanced turbulent fluxes and increased atmospheric water vapor. These results highlight the role of small‐scale air–sea interaction processes that could have a strong impact on Mediterranean cyclones' path and intensity.

Signature of the western boundary currents in local climate variability.

The western boundary currents are characterized by narrow, intense ocean jets and are among the most energetic phenomena in the world ocean. The importance of the western boundary currents to the mean climate is well established: they transport vast quantities of heat from the subtropics to the midlatitudes1, and they govern the structure of the climatological mean surface winds2-6, precipitation4-6 and extratropical storm tracks7-13. Their importance to climate variability is much less clear, as the tropospheric response to extratropical sea surface temperature (SST) variability is generally modest relative to the internal variability in the midlatitude atmosphere12-14. Here we exploit novel local analyses based on high-spatial-resolution data to demonstrate that SST variability in the western boundary currents has a more robust signature in climate variability than has been indicated in previous work. Our results indicate that warm SST anomalies in the major boundary currents of both hemispheres are associated with a distinct signature of locally enhanced precipitation and rising motion anomalies that extend throughout the depth of the troposphere. The tropospheric signature closely mirrors that of ocean dynamical processes in the boundary currents. Thus, the findings indicate a distinct and robust pathway through which extratropical ocean dynamical processes influence local climate variability. The observational relationships are also reproducible in Earth system model simulations but only when the simulations are run at high spatial resolution.

SST trends in certain areas of the Barents Sea in the winter season and mechanisms of their formation

The climate changes observed over the past few decades are most clearly manifested in the Arctic Ocean. Sea surface temperature (SST) is one of the most reliable indicators of climate change. In this paper we analyze the changes of winter SST for the western, northeastern and southeastern regions of the Barents Sea and examine the relationship of the emerging STS trends with the influence of various external factors. The working data set is represented by average monthly SST values taken from the ERA-5 reanalysis for the period 1949–2023 with a spatial resolution of 0.25×0.25° and average water temperature values on the Kola Meridian section in the 0–50 m layer. Additionally, the Arctic Oscillation (AO), Arctic Dipole (AD) and Atlantic Multidecadal Oscillation (AMO) indices were used as external factors that may affect SST variability. The time series analyzed was divided into three periods: 1949–1969, 1970–1990, 1991–2023, where the variability of the analyzed parameters was different. Thus, in the first period the trend in SST changes was negative, for the second period it was slightly negative or neutral, and for the third period it was positive. It is shown that SST in all the regions of the Barents Sea has undergone significant changes, which were most noticeable in the “warm” period of 1991–2023, when the rate of SST increasing was up to 10·10-2 °C/year in areas under the warm Atlantic water influence. The analysis of SST variability in the Barents Sea shows that the positive anomalies observed in the recent years are most likely associated with the changes in the atmospheric circulation. The Wavelet coherence analysis showed the closest agreement between the changes in the sea surface temperature and the AD index in the winter season, and with the AMO index.

Distinct seasonal changes and precession forcing of surface and subsurface temperatures in the mid-latitudinal North Atlantic during the onset of the Late Pliocene

Abstract. The Late Pliocene marks the intensification of Northern Hemisphere glaciation (iNHG), offering a unique opportunity to study climate evolution and ice-sheet-related feedback mechanisms. In this study, we present high-resolution Mg / Ca-based sea surface temperatures (SSTs) and subsurface temperatures (SubTs) derived from the foraminiferal species Globigerinoides ruber and Globorotalia hirsuta, respectively, at the Integrated Ocean Drilling Program (IODP) Expedition 306 Site U1313 in the mid-latitudinal North Atlantic during the early Late Pliocene, 3.65–3.37 million years ago (Ma). We find distinct differences between our new G. ruber Mg / Ca-based SST record and previously published alkenone-based SST records from the same location. These discrepancies in both absolute values and variations highlight distinctly different seasonal influences on the proxies. The G. ruber Mg / Ca-based SST data were primarily influenced by local summer insolation, showing a dominant precession cycle. Conversely, the variations in alkenone-based SST, dominated by the obliquity and lacking the precession cycle, are found to be more indicative of cold-season changes, despite previous interpretations of these records as reflecting annual mean temperatures. A simultaneous decline in Mg / Ca-based SST and SubT records from 3.65 to 3.5 Ma suggests a diminished poleward oceanic heat transport, implying a weakening of the North Atlantic Current (NAC). A comparison with Early Pleistocene G. ruber Mg / Ca-based SST records shows a shift in the dominant climatic cycle from precession to obliquity, alongside a marked increase in amplitude, indicating an enhanced influence of obliquity cycles correlated with the expansion of Northern Hemisphere ice sheets.

Sensitive response of erosion and weathering to the Indian Summer Monsoon changes in South Asia during Dansgaard-Oeschger oscillations

Continental weathering plays a key role in regulating the long-term climate stability via negative feedback. However, whether and how erosion and weathering respond to rapid climate change (e.g millennial scale) in large river basins remains unclear, partly due to the lack of archives with robust age models and high sampling resolution. As one of the largest sediment source-to-sink systems, the Himalaya-Ganga/Brahmaputra River-Bay of Bengal system is considered as an ideal laboratory for examining the weathering-climate relationships over the past. Here we present the elemental and mineralogical data from well age-constrained core sediments in the Bay of Bengal, which document the erosion and weathering conditions in South Asia during the last glacial period. All proxies generally show larger values during glacial interstadials than the stadials, consistent with the millennial variabilities of regional sea surface temperature (SST) and the Indian Summer Monsoon (ISM) strength. We confirm that the erosion and chemical weathering signals in Ganga/Brahmaputra River basins, which sensitively respond to the ISM strength and/or temperature change, could be propagated through the large river systems and faithfully preserved in the South Asian marginal seas, without noticeable time lag on a millennial scale. Our findings provide valuable data for response of erosion and weathering dynamics in a large river basin to rapid climate changes, and highlight the immense potential of monsoon-driven continental silicate weathering in floodplains as a substantial short-term carbon sink, thus underscoring a pivotal natural carbon sequestration mechanism amidst the escalating threat of global warming.