skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Detecting Paleoclimate Transitions With Laplacian Eigenmaps of Recurrence Matrices (LERM)
Abstract Paleoclimate records can be considered low‐dimensional projections of the climate system that generated them. Understanding what these projections tell us about past climates, and changes in their dynamics, is a main goal of time series analysis on such records. Laplacian eigenmaps of recurrence matrices (LERM) is a novel technique using univariate paleoclimate time series data to indicate when notable shifts in dynamics have occurred. LERM leverages time delay embedding to construct a manifold that is mappable to the attractor of the climate system; this manifold can then be analyzed for significant dynamical transitions. Through numerical experiments with observed and synthetic data, LERM is applied to detect both gradual and abrupt regime transitions. Our paragon for gradual transitions is the Mid‐Pleistocene Transition (MPT). We show that LERM can robustly detect gradual MPT‐like transitions for sufficiently high signal‐to‐noise (S/N) ratios, though with a time lag related to the embedding process. Our paragon of abrupt transitions is the “8.2 ka” event; we find that LERM is generally robust at detecting 8.2 ka‐like transitions for sufficiently high S/N ratios, though edge effects become more influential. We conclude that LERM can usefully detect dynamical transitions in paleogeoscientific time series, with the caveat that false positive rates are high when dynamical transitions are not present, suggesting the importance of using multiple records to confirm the robustness of transitions. We share an open‐source Python package to facilitate the use of LERM in paleoclimatology and paleoceanography.  more » « less
Award ID(s):
2002556
PAR ID:
10537510
Author(s) / Creator(s):
; ; ;
Publisher / Repository:
American Geophysical Union
Date Published:
Journal Name:
Paleoceanography and Paleoclimatology
Volume:
39
Issue:
1
ISSN:
2572-4517
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The 8.2 ka event is the most significant global climate anomaly of the Holocene epoch, but a lack of records from Mainland Southeast Asia (MSEA) currently limits our understanding of the spatial and temporal extent of the climate response. A newly developed speleothem record from Tham Doun Mai Cave, Northern Laos provides the first high resolution record of this event in MSEA. Our multiproxy record (d18O, d13C, Mg/Ca, Sr/Ca, and petrographic data), anchored in time by 9 U-Th ages, reveals a significant reduction in local rainfall amount and weakening of the monsoon at the event onset at ~8.29 +/- 0.03 ka BP. This response lasts for a minimum of ~170 years, similar to event length estimates from other speleothem d18O monsoon records. Interestingly, however, our d13C and Mg/Ca data, proxies for local hydrology, show that abrupt changes to local rainfall amounts began decades earlier (~70 years) than registered in the d18O. Moreover, the d13C and Mg/Ca also show that reductions in rainfall continued for at least ~200 years longer than the weakening of the monsoon inferred from the d18O. Our interpretations suggest that drier conditions brought on by the 8.2 ka event in MSEA were felt beyond the temporal boundaries defined by d18O-inferred monsoon intensity, and an initial wet period (or precursor event) may have preceded the local drying. Most existing Asian Monsoon proxy records of the 8.2 ka event may lack the resolution and/or multiproxy information necessary to establish local and regional hydrological sensitivity to abrupt climate change. 
    more » « less
  2. null (Ed.)
    In this paper, a new decadal resolution stalagmite δ18O record covering 10.4–6.5 ka BP from Kulishu cave in Beijing, north China is presented in combination with the published stalagmite δ18 O record covering 10.4–14.0 ka BP in the same cave. Five significant monsoon collapses were identified around 11.5, 11.0, 10.0, 9.4, and 8.2 ka BP as well as three smaller ones around 10.3, 9.0, and 8.6 ka BP. The weak monsoon episodes around 8.6 and 8.2 ka BP form the two-step structure of the 8.2 ka event. All monsoon collapses, coeval with the cooling in northern high-latitude records, are correlated with Lakes Agassiz-Ojibway outbursts. Thus, our data support the idea of freshwater forcing of abrupt climate anomalies during the early Holocene. Nevertheless, the decreased irradiance together with freshwater outburst may account for the 9.2/9.3 ka event, which is expressed more significantly in low-latitude records. 
    more » « less
  3. Abstract The 8.2 ka event is the most significant global climate anomaly of the Holocene epoch, but a lack of records from Mainland Southeast Asia (MSEA) currently limits our understanding of the spatial and temporal extent of the climate response. A newly developed speleothem record from Tham Doun Mai Cave, Northern Laos provides the first high‐resolution record of this event in MSEA. Our multiproxy record (δ18O, δ13C, Mg/Ca, Sr/Ca, and petrographic data), anchored in time by 9 U‐Th ages, reveals a significant reduction in local rainfall amount and weakening of the monsoon at the event onset at ∼8.29 ± 0.03 ka BP. This response lasts for a minimum of ∼170 years, similar to event length estimates from other speleothem δ18O monsoon records. Interestingly, however, our δ13C and Mg/Ca data, proxies for local hydrology, show that abrupt changes to local rainfall amounts began decades earlier (∼70 years) than registered in the δ18O. Moreover, the δ13C and Mg/Ca also show that reductions in rainfall continued for at least ∼200 years longer than the weakening of the monsoon inferred from the δ18O. Our interpretations suggest that drier conditions brought on by the 8.2 ka event in MSEA were felt beyond the temporal boundaries defined by δ18O‐inferred monsoon intensity, and an initial wet period (or precursor event) may have preceded the local drying. Most existing Asian Monsoon proxy records of the 8.2 ka event may lack the resolution and/or multiproxy information necessary to establish local and regional hydrological sensitivity to abrupt climate change. 
    more » « less
  4. Glacial-interglacial transitions and abrupt millennial-scale events are the most prominent features in many paleoclimate records. Understanding these oscillations requires high-resolution time series from multiple locations to constrain the latitudinal response to forcings. Few high-resolution records exist from the Southern Hemisphere tropics that predate the last two glaciations. We present a high-resolution speleothem oxygen and carbon isotope record from Huagapo Cave in the Central Peruvian Andes covering Marine Isotope Stage (MIS) 8 glacial and MIS 9 interglacial (339 to 249 ka). Uranium-series dates on three stalagmites (n=18) with small age uncertainty ±1% allows us to resolve abrupt climate events similar in structure and duration to Dansgaard-Oescchger and Heinrich events. The South American Summer Monsoon (SASM) controls modern hydroclimate variability in the Andes, and previous records from Huagapo Cave have provided records of past SASM variability. Termination three (T-III) in our record has a steep increase in δ18O values of 5‰, punctuated by two stadial event decreases of ~3‰ (S8.1 and S8.2). This pattern is mirrored in the δ13C record, indicating that these millennial-scale events record hydroclimate and vegetation productivity changes. The same structure as our T-III record is found in other records globally, where they are noted to be Heinrich-like events. Frequency analysis indicates that the occurrence of these abrupt events changes between glacial cycles. Precession is weakly expressed in the δ18O record during MIS 8; similar to speleothem records from the region dating to the Last Glacial Maximum (LGM). Global ice cover and sea levels were similar in the LGM and MIS 8, but the Milankovitch insolation forcing differed. This change in SASM behavior is not observed in the East Asian monsoon, where the precession signal is dominant throughout. Interglacial precessional control is apparent during the latter half of MIS 9 and during Huagapo Cave intervals dating to MIS 6 and 7. These data indicate that the response to high-latitude forcing in the Southern Hemisphere tropics fluctuates through time, and potential explanations for low-latitude sensitivity to forcing factors are further explored. 
    more » « less
  5. Abstract. Reconstructions of past temperature and precipitation are fundamental to modeling the Greenland Ice Sheet and assessing its sensitivity to climate. Paleoclimate information is sourced from proxy records and climate-model simulations; however, the former are spatially incomplete while the latter are sensitive to model dynamics and boundary conditions. Efforts to combine these sources of information to reconstruct spatial patterns of Greenland climate over glacial–interglacial cycles have been limited by assumptions of fixed spatial patterns and a restricted use of proxy data. We avoid these limitations by using paleoclimate data assimilation to create independent reconstructions of mean-annual temperature and precipitation for the last 20 000 years. Our method uses oxygen isotope ratios of ice and accumulation rates from long ice-core records and extends this information to all locations across Greenland using spatial relationships derived from a transient climate-model simulation. Standard evaluation metrics for this method show that our results capture climate at locations without ice-core records. Our results differ from previous work in the reconstructed spatial pattern of temperature change during abrupt climate transitions; this indicates a need for additional proxy data and additional transient climate-model simulations. We investigate the relationship between precipitation and temperature, finding that it is frequency dependent and spatially variable, suggesting that thermodynamic scaling methods commonly used in ice-sheet modeling are overly simplistic. Our results demonstrate that paleoclimate data assimilation is a useful tool for reconstructing the spatial and temporal patterns of past climate on timescales relevant to ice sheets. 
    more » « less