High-resolution record of Northern Hemisphere climate extending into the last interglacial period

Earthquake supercycles and fault interaction over the past 32 ka in the Lake Ebinur area, Xinjiang, China

Earthquake prediction and disaster assessment in tectonically active regions require a continuous and complete regional seismic archive, which is commonly difficult to obtain, especially for prehistoric records. Here, high-resolution analysis of the sedimentary sequence from Lake Ebinur in Xinjiang revealed a detailed history of environment evolution since 32 ka ago. Both the Cl content and ultrafine proportion revealed the changing climate: the climate was relatively dry with low lake-water volumes from 32 to 12 ka, while the climate became warmer and wetter since 12 ka. In addition, eight earthquakes were identified by comprehensive analysis of grain size and geochemical element proxies, showing more than two seismic supercycles, with gaps of ∼10.4 ka; these gaps are much larger than those inferred previously (∼4-7 ka). Notably, these seismic events exhibited a pattern of mutual transmittance between the BoA and Jinghenan faults. Such fault interaction can occur in the Lake Ebinur area because it is dominated by weak lithosphere in which strain is easily accumulated and released; the interaction can also be attributed to the unique spatial distribution and immature nature of both faults. Combined with trenching investigations, our high-resolution analysis of lacustrine sediments can reveal a complete history of tectonic activity, which can efficiently serve regional earthquake prediction and disaster assessment.

Legacies of millennial-scale climate oscillations in contemporary biodiversity in eastern North America

The Atlantic meridional overturning circulation (AMOC) has caused significant climate changes over the past 90 000 years. Prior work has hypothesized that these millennial-scale climate variations effected past and contemporary biodiversity, but the effects are understudied. Moreover, few biogeographic models have accounted for uncertainties in palaeoclimatic simulations of millennial-scale variability. We examine whether refuges from millennial-scale climate oscillations have left detectable legacies in the patterns of contemporary species richness in eastern North America. We analyse 13 palaeoclimate estimates from climate simulations and proxy-based reconstructions as predictors for the contemporary richness of amphibians, passerine birds, mammals, reptiles and trees. Results suggest that past climate changes owing to AMOC variations have left weak but detectable imprints on the contemporary richness of mammals and trees. High temperature stability, precipitation increase, and an apparent climate fulcrum in the southeastern United States across millennial-scale climate oscillations aligns with high biodiversity in the region. These findings support the hypothesis that the southeastern United States may have acted as a biodiversity refuge. However, for some taxa, the strength and direction of palaeoclimate-richness relationships varies among different palaeoclimate estimates, pointing to the importance of palaeoclimatic ensembles and the need for caution when basing biogeographic interpretations on individual palaeoclimate simulations. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

Using ancient sedimentary DNA to forecast ecosystem trajectories under climate change

Ecosystem response to climate change is complex. In order to forecast ecosystem dynamics, we need high-quality data on changes in past species abundance that can inform process-based models. Sedimentary ancient DNA (sedaDNA) has revolutionised our ability to document past ecosystems' dynamics. It provides time series of increased taxonomic resolution compared to microfossils (pollen, spores), and can often give species-level information, especially for past vascular plant and mammal abundances. Time series are much richer in information than contemporary spatial distribution information, which have been traditionally used to train models for predicting biodiversity and ecosystem responses to climate change. Here, we outline the potential contribution of sedaDNA to forecast ecosystem changes. We showcase how species-level time series may allow quantification of the effect of biotic interactions in ecosystem dynamics, and be used to estimate dispersal rates when a dense network of sites is available. By combining palaeo-time series, process-based models, and inverse modelling, we can recover the biotic and abiotic processes underlying ecosystem dynamics, which are traditionally very challenging to characterise. Dynamic models informed by sedaDNA can further be used to extrapolate beyond current dynamics and provide robust forecasts of ecosystem responses to future climate change. This article is part of the theme issue 'Ecological novelty and planetary stewardship: biodiversity dynamics in a transforming biosphere'.

A marine record of Patagonian ice sheet changes over the past 140,000 years

Terrestrial glacial records from the Patagonian Andes and New Zealand Alps document quasi-synchronous Southern Hemisphere-wide glacier advances during the late Quaternary. However, these records are inherently incomplete. Here, we provide a continuous marine record of western-central Patagonian ice sheet (PIS) extent over a complete glacial-interglacial cycle back into the penultimate glacial (~140 ka). Sediment core MR16-09 PC03, located at 46°S and ~150 km offshore Chile, received high terrestrial sediment and meltwater input when the central PIS extended westward. We use biomarkers, foraminiferal oxygen isotopes, and major elemental data to reconstruct terrestrial sediment and freshwater input related to PIS variations. Our sediment record documents three intervals of general PIS marginal fluctuations, during Marine Isotope Stage (MIS) 6 (140 to 135 ka), MIS 4 (~70 to 60 ka), and late MIS 3 to MIS 2 (~40 to 18 ka). These higher terrigenous input intervals occurred during sea-level low stands, when the western PIS covered most of the Chilean fjords, which today retain glaciofluvial sediments. During these intervals, high-amplitude phases of enhanced sediment supply occur at millennial timescales, reflecting increased ice discharge most likely due to a growing PIS. We assign the late MIS 3 to MIS 2 phases and, by inference, older advances to Antarctic cold stages. We conclude that the increased sediment/meltwater release during Southern Hemisphere millennial-scale cold phases was likely related to higher precipitation caused by enhanced westerly winds at the northwestern margin of the PIS. Our records complement terrestrial archives and provide evidence for PIS climate sensitivity.

Dansgaard-Oeschger cycles of the penultimate and last glacial period recorded in stalagmites from Türkiye

The last glacial period is characterized by abrupt climate oscillations, also known as Dansgaard-Oeschger (D-O) cycles. However, D-O cycles remain poorly documented in climate proxy records covering the penultimate glacial period. Here we present highly resolved and precisely dated speleothem time series from Sofular Cave in northern Türkiye to provide clear evidence for D-O cycles during Marine Isotope Stage (MIS) 6 as well as MIS 2-4. D-O cycles are most clearly expressed in the Sofular carbon isotope time series, which correlate inversely with regional sea surface temperature (SST) records from the Black Sea. The pacing of D-O cycles is almost twice as long during MIS 6 compared to MIS 2-4, and could be related to a weaker Atlantic Meridional Overturning Circulation (AMOC) and a different mean climate during MIS 6 compared to MIS 2-4, leading most likely to a higher threshold for the occurrence of D-O cycles.

Exploring different methodological approaches to unlock paleobiodiversity in peat profiles using ancient DNA

Natural and human-induced environmental changes deeply affected terrestrial ecosystems throughout the Holocene. Paleoenvironmental reconstructions provide information about the past and allow us to predict/model future scenarios. Among potential records, peat bogs are widely used because they present a precise stratigraphy and act as natural archives of highly diverse organic remains. Over the decades, several techniques have been developed to identify debris occurring in peat, including their morphological description. However, this is strongly constrained by the researcher's ability to distinguish residues at the species level, which typically requires many years of experience. In addition, potential contamination hampers using these techniques to obtain information from organisms such as fungi or bacteria. Environmental DNA metabarcoding and shotgun metagenome sequencing could represent a solution to detect specific groups of organisms without any a priori knowledge of their characteristics and/or to identify organisms that have rarely been considered in previous investigations. Moreover, shotgun metagenomics may allow the identification of bacteria and fungi (including both yeast and filamentous life forms), ensuring discrimination between ancient and modern organisms through the study of deamination/damage patterns. In the present review, we aim to i) present the state-of-the-art methodologies in paleoecological and paleoclimatic studies focusing on peat core analyses, proposing alternative approaches to the classical morphological identification of plant residues, and ii) suggest biomolecular approaches that will allow the use of proxies such as invertebrates, fungi, and bacteria, which are rarely employed in paleoenvironmental reconstructions.

Stable-isotope ratios (δ18O and δD) in a firn core from West Antarctica

A 22.48 m long ice core (BR-IC-4) was collected in the West Antarctic Ice Sheet (at 83°58'59.4" S, 80°07'01.4" W, 1,295 m above the sea level) during the Austral summer of 2004-2005, as a contribution to the International Trans-Antarctic Expedition program. The isotopic composition (δD and δ18O) of 599 samples, corresponding to the upper 12.98 m of the ice core, was determined by gas source mass spectrometry and cavity ring-down spectroscopy. Relative dating was based on the isotopic ratios and major ions (MS-, Na+, nssSO4 2-) and trace elements (Na, S, Sr) concentrations. The record covers approximately 13 years - from 1990 to 2003. The mean accumulation rate of 0.48 ± 0.09 m water equivalent per year (m eq H2O a1) is relatively high for the geographical area and possibly results from snowdrifting from near areas, as attested by ice glaze surfaces in other sites in the region. The stable isotope δD content varies between -367.90‰ and 256.30‰ (mean -314.42 ± 19.01‰); and δ18O ranges from -44.96‰ to 35.08‰ (mean -39.95 ± 2.05‰). Deuterium excess values (mean 3.70 ± 1.54‰) indicate episodic intense oceanic evaporation and high relative humidity in the moisture sources.

The firing temperatures of burnt clay from the Chinese neolithic cultural relics and its paleoenvironmental imprints

The impact of human-induced fires on the surrounding environment has been particularly significant during the Anthropocene epoch. Neolithic burnt clay, which is widely distributed in the archaeological relics of ancient civilizations across the Eurasian continent, provides pivotal information about the ancient firing stories. However, understanding of the paleoenvironmental imprints of fire in burnt clay has been largely limited by the lack of sufficient analytical data on comprehensive knowledge of ancient firing conditions. In this study, a detailed magnetic analysis was conducted on burnt clay materials from a Neolithic site in Fujian, southeastern China, which presented a burnt clay-based record of the Neolithic firing temperature in relation to paleoenvironmental conditions. Based on magnetic analysis, the ancient firing temperature was determined to be approximately 620 °C, which is comparable with other records from Eurasian Neolithic sites. Frequency-dependent magnetic susceptibility provides alternative indicative information for ancient firing conditions in addition to conventional magnetic susceptibility. Furthermore, magnetic properties of burnt clay may decode the in-situ source characteristics with respect to geological background. In addition, a potential link between temporal variations in ancient firing temperatures in burnt clay and surrounding paleoenvironmental changes is tentatively interpreted by local environmental feedback of temperature-moisture conditions and anthropogenic activity. This study further confirms the archaeological potential of thermomagnetic properties as useful indicators in paleoenvironmental studies. More work combining paleoenvironmental and archaeological archives is critically essential to understand ancient firing history in the context of environmental conditions.

Climate amelioration, abrupt vegetation recovery, and the dispersal of Homo sapiens in Baikal Siberia

The dispersal of Homo sapiens in Siberia and Mongolia occurred by 45 to 40 thousand years (ka) ago; however, the climatic and environmental context of this event remains poorly understood. We reconstruct a detailed vegetation history for the Last Glacial period based on pollen spectra from Lake Baikal. While herb and shrub taxa including Artemisia and Alnus dominated throughout most of this period, coniferous forests rapidly expanded during Dansgaard-Oeschger (D-O) events 14 (55 ka ago) and 12 to 10 (48 to 41 ka ago), with the latter presenting the strongest signal for coniferous forest expansion and Picea trees, indicating remarkably humid conditions. These abrupt forestation events are consistent with obliquity maxima, so that we interpret last glacial vegetation changes in southern Siberia as being driven by obliquity change. Likewise, we posit that major climate amelioration and pronounced forestation precipitated H. sapiens dispersal into Baikal Siberia 45 ka ago, as chronicled by the appearance of the Initial Upper Paleolithic.

Antarctic evidence for an abrupt northward shift of the Southern Hemisphere westerlies at 32 ka BP

High-resolution ice core records from coastal Antarctica are particularly useful to inform our understanding of environmental changes and their drivers. Here, we present a decadally resolved record of sea-salt sodium (a proxy for open-ocean area) and non-sea salt calcium (a proxy for continental dust) from the well-dated Roosevelt Island Climate Evolution (RICE) core, focusing on the time period between 40-26 ka BP. The RICE dust record exhibits an abrupt shift towards a higher mean dust concentration at 32 ka BP. Investigating existing ice-core records, we find this shift is a prominent feature across Antarctica. We propose that this shift is linked to an equatorward displacement of Southern Hemisphere westerly winds. Subsequent to the wind shift, data suggest a weakening of Southern Ocean upwelling and a decline of atmospheric CO2 to lower glacial values, hence making this shift an important glacial climate event with potentially important insights for future projections.

Last glacial hydroclimate variability in the Yucatán Peninsula not just driven by ITCZ shifts

We reconstructed hydroclimate variability in the Yucatán Peninsula (YP) based on stalagmite oxygen and carbon isotope records from a well-studied cave system located in the northeastern YP, a region strongly influenced by Caribbean climate dynamics. The new stalagmite isotopic records span the time interval between 43 and 26.6 ka BP, extending a previously published record from the same cave system covering the interval between 26.5 and 23.2 ka BP. Stalagmite stable isotope records show dominant decadal and multidecadal variability, and weaker variability on millennial timescales. These records suggest significant precipitation declines in the broader Caribbean region during Heinrich events 4 and 3 of ice-rafted discharge into the North Atlantic, in agreement with the antiphase pattern of precipitation variability across the equator suggested by previous studies. On millennial timescales, the stalagmite isotope records do not show the distinctive saw-tooth pattern of climate variability observed in Greenland during Dansgaard-Oeschger (DO) events, but a pattern similar to North Atlantic sea surface temperature (SST) variability. We propose that shifts in the mean position of the Intertropical Convergence Zone (ITCZ), per se, are not the dominant driver of last glacial hydroclimate variability in the YP on millennial timescales but instead that North Atlantic SSTs played a dominant role. Our results support a negative climate feedback mechanism whereby large low latitude precipitation deficits resulting from AMOC slowdown would lead to elevated salinity in the Caribbean and ultimately help reactivate AMOC and Caribbean precipitation. However, because of the unique drivers of future climate in the region, predicted twenty-first century YP precipitation reductions are unlikely to be modulated by this negative feedback mechanism.

Pre-Younger Dryas megafaunal extirpation at Rancho La Brea linked to fire-driven state shift

The cause, or causes, of the Pleistocene megafaunal extinctions have been difficult to establish, in part because poor spatiotemporal resolution in the fossil record hinders alignment of species disappearances with archeological and environmental data. We obtained 172 new radiocarbon dates on megafauna from Rancho La Brea in California spanning 15.6 to 10.0 thousand calendar years before present (ka). Seven species of extinct megafauna disappeared by 12.9 ka, before the onset of the Younger Dryas. Comparison with high-resolution regional datasets revealed that these disappearances coincided with an ecological state shift that followed aridification and vegetation changes during the Bølling-Allerød (14.69 to 12.89 ka). Time-series modeling implicates large-scale fires as the primary cause of the extirpations, and the catalyst of this state shift may have been mounting human impacts in a drying, warming, and increasingly fire-prone ecosystem.

The spatiotemporal extent of the Green Sahara during the last glacial period

The Sahara Desert, one of today's most inhospitable environments, has known periods of enhanced precipitation that supported pre-historic humans. However, the Green Sahara timing and moisture sources are not well known due to limited paleoclimate information. Here, we present a multi-proxy (δ¹⁸O, δ¹³C, Δ¹⁷O, and trace elements) speleothem-based climate record from Northwest (NW) Africa. Our data document two Green Sahara periods during Marine Isotope Stage (MIS) 5a and the Early to Mid-Holocene. Consistency with paleoclimate records across North Africa highlights the east-west geographical extent of the Green Sahara, whereas millennial-scale North Atlantic cooling (Heinrich) events consistently resulted in drier conditions. We demonstrate that an increase in westerly-originating winter precipitation during MIS5a resulted in favorable environmental conditions. The comparison of paleoclimate data with local archaeological sequences highlights the abrupt climate deterioration and the decline in human density in NW Africa during the MIS5-4 transition, which suggests climate-forced dispersals of populations, with possible implications for pathways into Eurasia.

Dynamic and thermodynamic influences on precipitation in Northeast Mexico on orbital to millennial timescales

The timing and mechanisms of past hydroclimate change in northeast Mexico are poorly constrained, limiting our ability to evaluate climate model performance. To address this, we present a multiproxy speleothem record of past hydroclimate variability spanning 62.5 to 5.1 ka from Tamaulipas, Mexico. Here we show a strong influence of Atlantic and Pacific sea surface temperatures on orbital and millennial scale precipitation changes in the region. Multiple proxies show no clear response to insolation forcing, but strong evidence for dry conditions during Heinrich Stadials. While these trends are consistent with other records from across Mesoamerica and the Caribbean, the relative importance of thermodynamic and dynamic controls in driving this response is debated. An isotope-enabled climate model shows that cool Atlantic SSTs and stronger easterlies drive a strong inter-basin sea surface temperature gradient and a southward shift in moisture convergence, causing drying in this region.

Multi-proxy constraints on Atlantic circulation dynamics since the last ice age

Uncertainties persist in the understanding of the Atlantic meridional overturning circulation and its response to external perturbations such as freshwater or radiative forcing. Abrupt reduction of the Atlantic circulation is considered a climate tipping point that may have been crossed when Earth's climate was propelled out of the last ice age. However, the evolution of the circulation since the Last Glacial Maximum (22-18 thousand years ago) remains insufficiently constrained due to model and proxy limitations. Here we leverage information from both a compilation of proxy records that track various aspects of the circulation and climate model simulations to constrain the Atlantic circulation over the past 20,000 years. We find a coherent picture of a shallow and weak Atlantic overturning circulation during the Last Glacial Maximum that reconciles apparently conflicting proxy evidence. Model-data comparison of the last deglaciation-starting from this new, multiple constrained glacial state-indicates a muted response during Heinrich Stadial 1 and that water mass geometry did not fully adjust to the strong reduction in overturning circulation during the comparably short Younger Dryas period. This demonstrates that the relationship between freshwater forcing and Atlantic overturning strength is strongly dependent on the climatic and oceanic background state.

Widespread glacier advances across the Tian Shan during Marine Isotope Stage 3 not supported by climate-glaciation simulations

Whether there were more extensive glaciations during the Marine Isotope Stage (MIS) 3 relative to MIS 2 across the Tian Shan in Central Asia is intensely debated because of the uncertainty in chronological data and fully understanding the driving mechanisms. To help resolve the ongoing debate, we assess the climate sensitivity of the glaciers and reconstruct the extent of glaciation during MIS 2 and 3 across the Tian Shan, using a glacier-resolving (250 × 250 m) ice sheet model asynchronously coupled with a global climate model. Our results demonstrate that the equilibrium-line altitude (ELA) over the Tian Shan decreases by ∼180 m for every 1 °C cooling under a modern precipitation regime, but precipitation reduction greatly lowers the sensitivity of the glaciers to temperature decrease (e.g., the effect of 2 °C cooling is broadly offset by a 50% decrease in precipitation). Under the modeled colder/drier-than-present climate, the model predicts an ELA depression (∆ELA) of ∼75 m (162 m) over the Tian Shan during MIS 3 interstadials (stadials). The extent of MIS 3 glaciation is much smaller than that during MIS 2 (i.e., ∆ELA = ∼726 m). The more extensive glaciation during MIS 2 is largely attributed to the enhanced summer cooling. Furthermore, through a site-to-site model-data comparison, we find that the closest match between the modeled glacier margin and the locations of the glacial deposits previously argued to be MIS 3 is generally achieved under MIS 2 climatic conditions. These results suggest more extensive glacier advances over the Tian Shan during MIS 2 than MIS 3 on a regional scale, although MIS 3 glaciation may still occur in individual glacier catchments. This pattern suggests general synchronicity with the timing of maximum Northern Hemisphere ice sheets during the last glacial, which should be further tested in a multimodel framework in the future.

Hydroclimate variability in the Caribbean during North Atlantic Heinrich cooling events (H8 and H9)

We present a speleothem record from western Cuba, spanning the period 98.7-84.9 ka BP. Our record shows two distinctive periods of high δ¹⁸O corresponding to dry and/or cold periods during 85-87.6 and 90.2-93.1 ka BP, synchronous with Heinrich events 8 and 9 (H8 and H9). Hence, we provide the first proxy evidence of the local Caribbean climate response to H8 and H9. Interestingly, H8 is more pronounced compared to H9, which may be a local response to lower temperatures in the North Atlantic resulting in a weak AMOC and reduced deep water formation, therefore a stronger south shift of the ITCZ. Our data complement existing speleothem records from western Cuba which, collectively, provide a nearly continuous paleoclimate time-series spanning the last 100 ka BP, indicating a consistent response to millennial-scale events as dry and/or cooler conditions. The comparison with regional paleoclimate records reveals an anti-phased relationship with South America, caused by the southern movements of the ITCZ during millennial-scale events which lead to dry conditions in the Caribbean and a stronger South American Monsoon System.

The extreme yet transient nature of glacial erosion

Ice can sculpt extraordinary landscapes, yet the efficacy of, and controls governing, glacial erosion on geological timescales remain poorly understood and contended, particularly across Polar continental shields. Here, we assimilate geophysical data with modelling of the Eurasian Ice Sheet - the third largest Quaternary ice mass that spanned 49°N to 82°N - to decipher its erosional footprint during the entire last ~100 ka glacial cycle. Our results demonstrate extreme spatial and temporal heterogeneity in subglacial erosion, with rates ranging from 0 to 5 mm a^-1 and a net volume equating to ~130,000 km³ of bedrock excavated to depths of ~190 m. A hierarchy of environmental controls ostensibly underpins this complex signature: lithology, topography and climate, though it is basal thermodynamics that ultimately regulates erosion, which can be variously protective, pervasive, or, highly selective. Our analysis highlights the remarkable yet fickle nature of glacial erosion - critically modulated by transient ice-sheet dynamics - with its capacity to impart a profound but piecemeal geological legacy across mid- and high latitudes.

Sea ice fluctuations in the Baffin Bay and the Labrador Sea during glacial abrupt climate changes

Sea ice decline in the North Atlantic and Nordic Seas has been proposed to contribute to the repeated abrupt atmospheric warmings recorded in Greenland ice cores during the last glacial period, known as Dansgaard-Oeschger (D-O) events. However, the understanding of how sea ice changes were coupled with abrupt climate changes during D-O events has remained incomplete due to a lack of suitable high-resolution sea ice proxy records from northwestern North Atlantic regions. Here, we present a subdecadal-scale bromine enrichment (Br enr ) record from the NEEM ice core (Northwest Greenland) and sediment core biomarker records to reconstruct the variability of seasonal sea ice in the Baffin Bay and Labrador Sea over a suite of D-O events between 34 and 42 ka. Our results reveal repeated shifts between stable, multiyear sea ice (MYSI) conditions during cold stadials and unstable, seasonal sea ice conditions during warmer interstadials. The shift from stadial to interstadial sea ice conditions occurred rapidly and synchronously with the atmospheric warming over Greenland, while the amplitude of high-frequency sea ice fluctuations increased through interstadials. Our findings suggest that the rapid replacement of widespread MYSI with seasonal sea ice amplified the abrupt climate warming over the course of D-O events and highlight the role of feedbacks associated with late-interstadial seasonal sea ice expansion in driving the North Atlantic ocean–climate system back to stadial conditions.

Giant clam (Tridacna) distribution in the Gulf of Oman in relation to past and future climate

The Oman upwelling zone (OUZ) creates an unfavorable environment and a major biogeographic barrier for many coral reef species, such as giant clams, thus promoting and maintaining faunal differences among reefs on the east and west side of the Arabian Peninsula. We record the former existence of Tridacna in the Gulf of Oman and review its stratigraphic distribution in the Persian Gulf to provide new insights on the connectivity of coral reef habitats around southern Arabia under changing climate and ocean conditions. Fossil shells were carbon-14 dated and employed as sclerochronological proxy archives. This reveals that the Omani population represents a last glacial colonization event during the Marine Isotope Stage 3 interstadial under colder-than-present temperatures and variable upwelling intensity linked to Dansgaard-Oeschger climate oscillations. It was favored by temperatures just above the lower threshold for the habitat-forming reef coral communities and instability of the upwelling barrier. We conclude that the distribution of Tridacna in the northern Arabian Sea is generally limited by either strong upwelling or cool sea surface temperature under gradually changing climate conditions at the interglacial-glacial scale. Opportunities for dispersal and temporary colonization existed only when there was a simultaneous attenuation of both limiting factors due to high-frequency climate variability. The OUZ will unlikely become a future climate change refuge for giant clams because they will be exposed either to thermal stress by rapid anthropogenic Indian Ocean warming or to unfavorable upwelling conditions.

Multiple carbon cycle mechanisms associated with the glaciation of Marine Isotope Stage 4

Here we use high-precision carbon isotope data (δ¹³C-CO₂) to show atmospheric CO₂ during Marine Isotope Stage 4 (MIS 4, ~70.5-59 ka) was controlled by a succession of millennial-scale processes. Enriched δ¹³C-CO₂ during peak glaciation suggests increased ocean carbon storage. Variations in δ¹³C-CO₂ in early MIS 4 suggest multiple processes were active during CO₂ drawdown, potentially including decreased land carbon and decreased Southern Ocean air-sea gas exchange superposed on increased ocean carbon storage. CO₂ remained low during MIS 4 while δ¹³C-CO₂ fluctuations suggest changes in Southern Ocean and North Atlantic air-sea gas exchange. A 7 ppm increase in CO₂ at the onset of Dansgaard-Oeschger event 19 (72.1 ka) and 27 ppm increase in CO₂ during late MIS 4 (Heinrich Stadial 6, ~63.5-60 ka) involved additions of isotopically light carbon to the atmosphere. The terrestrial biosphere and Southern Ocean air-sea gas exchange are possible sources, with the latter event also involving decreased ocean carbon storage.

Summary for general audience: We used carbon stable isotope data from an Antarctic ice core to evaluate which mechanisms caused changes in atmospheric carbon dioxide 74-59 thousand years ago, including a ~40 ppm decrease at the beginning of the last ice age.

Subglacial precipitates record Antarctic ice sheet response to late Pleistocene millennial climate cycles

Ice cores and offshore sedimentary records demonstrate enhanced ice loss along Antarctic coastal margins during millennial-scale warm intervals within the last glacial termination. However, the distal location and short temporal coverage of these records leads to uncertainty in both the spatial footprint of ice loss, and whether millennial-scale ice response occurs outside of glacial terminations. Here we present a >100kyr archive of periodic transitions in subglacial precipitate mineralogy that are synchronous with Late Pleistocene millennial-scale climate cycles. Geochemical and geochronologic data provide evidence for opal formation during cold periods via cryoconcentration of subglacial brine, and calcite formation during warm periods through the addition of subglacial meltwater originating from the ice sheet interior. These freeze-flush cycles represent cyclic changes in subglacial hydrologic-connectivity driven by ice sheet velocity fluctuations. Our findings imply that oscillating Southern Ocean temperatures drive a dynamic response in the Antarctic ice sheet on millennial timescales, regardless of the background climate state.

Piccione et al find evidence for Antarctic ice sheet instability driven by millennial cycles in Southern Ocean temperature, providing clues for the mechanisms that link climate change and rapid Antarctic ice loss events.

Evolution of tropical land temperature across the last glacial termination

The tropical West Pacific hosts the warmest part of the surface ocean and has a considerable impact on the global climate system. Reconstructions of past temperature in this region can elucidate climate connections between the tropics and poles and the sensitivity of tropical temperature to greenhouse forcing. However, existing data are equivocal and reliable information from terrestrial archives is particularly sparse. Here we constrain the magnitude and timing of land temperature change in the tropical West Pacific across the last deglaciation using an exceptionally precise paleothermometer applied to a well-dated stalagmite from Northern Borneo. We show that the cave temperature increased by 4.4 ± 0.3 °C (2 SEM) from the Last Glacial Maximum to the Holocene, amounting to 3.6 ± 0.3 °C (2 SEM) when correcting for sea-level induced cave altitude change. The warming closely follows atmospheric CO₂ and Southern Hemisphere warming. This contrasts with hydroclimate, as reflected by drip water δ¹⁸O, which responds to Northern Hemisphere cooling events in the form of prominent drying, while temperature was rising. Our results thus show a close response of tropical temperature to greenhouse forcing, independent of shifts in the tropical circulation patterns.

Ancient drip water in a Borneo stalagmite reveals a strong land temperature rise across the last glacial termination in close correspondence with atmospheric CO2, and an intriguing decoupling between tropical temperature and hydroclimate.

Past, present, and future predictions on the suitable habitat of the Slender racer (Orientocoluber spinalis) using species distribution models

Species distribution models (SDMs) across past, present, and future timelines provide insights into the current distribution of these species and their reaction to climate change. Specifically, if a species is threatened or not well-known, the information may be critical to understand that species. In this study, we computed SDMs for Orientocoluber spinalis, a monotypic snake genus found in central and northeast Asia, across the past (last interglacial, last glacial maximum, and mid-Holocene), present, and future (2070s). The goal of the study was to understand the shifts in distribution across time, and the climatic factors primarily affecting the distribution of the species. We found the suitable habitat of O. spinalis to be persistently located in cold-dry winter and hot summer climatic areas where annual mean temperature, isothermality, and annual mean precipitation were important for suitable habitat conditions. Since the last glacial maximum, the suitable habitat of the species has consistently shifted northward. Despite the increase in suitable habitat, the rapid alterations in weather regimes because of climate change in the near future are likely to greatly threaten the southern populations of O. spinalis, especially in South Korea and China. To cope with such potential future threats, understanding the ecological requirements of the species and developing conservation plans are urgently needed.

Human activities have reduced plant diversity in eastern China over the last two millennia

Understanding the history and regional singularities of human impact on vegetation is key to developing strategies for sustainable ecosystem management. In this study, fossil and modern pollen datasets from China are employed to investigate temporal changes in pollen composition, analogue quality, and pollen diversity during the Holocene. Anthropogenic disturbance and vegetation's responses are also assessed. Results reveal that pollen assemblages from non-forest communities fail to provide evidence of human impact for the western part of China (annual precipitation less than 400 mm and/or elevation more than 3000 m.a.s.l.), as inferred from the stable quality of modern analogues, principal components, and diversity of species and communities throughout the Holocene. For the eastern part of China, the proportion of fossil pollen spectra with good modern analogues increases from ca. 50% to ca. 80% during the last 2 millennia, indicating an enhanced intensity of anthropogenic disturbance on vegetation. This disturbance has caused the pollen spectra to become taxonomically less diverse over space (reduced abundances of arboreal taxa and increased abundances of herbaceous taxa), highlighting a reduced south-north differentiation and divergence from past vegetation between regions in the eastern part of China. We recommend that care is taken in eastern China when basing the development of ecosystem management strategies on vegetation changes in the region during the last 2000 years, since humans have significantly disturbed the vegetation during this period.

Subsurface ocean warming preceded Heinrich Events

Although the global environmental impact of Laurentide Ice-Sheet destabilizations on glacial climate during Heinrich Events is well-documented, the mechanism driving these ice-sheet instabilities remains elusive. Here we report foraminifera-based subsurface (~150 m water depth) ocean temperature and salinity reconstructions from a sediment core collected in the western subpolar North Atlantic, showing a consistent pattern of rapid subsurface ocean warming preceding the transition into each Heinrich Event identified in the same core of the last 27,000 years. These results provide the first solid evidence for the massive accumulation of ocean heat near the critical depth to trigger melting of marine-terminating portions of the Laurentide Ice Sheet around Labrador Sea followed by Heinrich Events. The repeated build-up of a subsurface heat reservoir in the subpolar Atlantic closely corresponds to times of weakened Atlantic Meridional Overturning Circulation, indicating a precursor role of ocean circulation changes for initiating abrupt ice-sheet instabilities during Heinrich Events. We infer that a weaker ocean circulation in future may result in accelerated interior-ocean warming of the subpolar Atlantic, which could be critical for the stability of modern, marine-terminating Arctic glaciers and the freshwater budget of the North Atlantic.

Late quaternary biotic homogenization of North American mammalian faunas

Biotic homogenization—increasing similarity of species composition among ecological communities—has been linked to anthropogenic processes operating over the last century. Fossil evidence, however, suggests that humans have had impacts on ecosystems for millennia. We quantify biotic homogenization of North American mammalian assemblages during the late Pleistocene through Holocene (~30,000 ybp to recent), a timespan encompassing increased evidence of humans on the landscape (~20,000–14,000 ybp). From ~10,000 ybp to recent, assemblages became significantly more homogenous (>100% increase in Jaccard similarity), a pattern that cannot be explained by changes in fossil record sampling. Homogenization was most pronounced among mammals larger than 1 kg and occurred in two phases. The first followed the megafaunal extinction at ~10,000 ybp. The second, more rapid phase began during human population growth and early agricultural intensification (~2,000–1,000 ybp). We show that North American ecosystems were homogenizing for millennia, extending human impacts back ~10,000 years.

Biotic homogenization, which is increased similarity in the composition of species among communities, is rising due to human activities. Using North American mammal fossil records from the past 30,000 years, this study shows that this phenomenon is ancient, beginning between 12,000 and 10,000 years ago with the extinction of the mammal megafauna.

Southern hemisphere forced millennial scale Indian summer monsoon variability during the late Pleistocene

Peninsular India hosts the initial rain-down of the Indian Summer Monsoon (ISM) after which winds travel further east inwards into Asia. Stalagmite oxygen isotope composition from this region, such as those from Belum Cave, preserve the vital signals of the past ISM variability. These archives experience a single wet season with a single dominant moisture source annually. Here we present high-resolution δ¹⁸O, δ¹³C and trace element (Mg/Ca, Sr/Ca, Ba/Ca, Mn/Ca) time series from a Belum Cave stalagmite spanning glacial MIS-6 (from ~ 183 to ~ 175 kyr) and interglacial substages MIS-5c-5a (~ 104 kyr to ~ 82 kyr). With most paleomonsoon reconstructions reporting coherent evolution of northern hemisphere summer insolation and ISM variability on orbital timescale, we focus on understanding the mechanisms behind millennial scale variability. Finding that the two are decoupled over millennial timescales, we address the role of the Southern Hemisphere processes in modulating monsoon strength as a part of the Hadley circulation. We identify several strong and weak episodes of ISM intensity during 104–82 kyr. Some of the weak episodes correspond to warming in the southern hemisphere associated with weak cross-equatorial winds. We show that during the MIS-5 substages, ISM strength gradually declined with millennial scale variability linked to Southern Hemisphere temperature changes which in turn modulate the strength of the Mascarene High.

Environmental and Hydrological Changes of Lake Coatetelco in Central Mesoamerica (Southwest Mexico) Over the Holocene and Comparison With Climatic Forcing

Elemental composition of the inorganic fraction, carbon isotopes, and C/N of organic matter from a new radiocarbon-dated sedimentary sequence collected from Lake Coatetelco (960 m asl) extend information about the environmental and hydrological conditions of central Mesoamerica from the earliest Holocene to the interval of first human settlements in the lake vicinity and nearby streams. Proxy-based reconstructions of erosion/runoff (Ti), water column salinity (CaCO3), sediment–water interactions (PIA), and provenance of organics (δ13Corg and C/N) showed similarities with the summer insolation modulated ITCZ position between ∼11.5 and 4.2 cal ka BP, and more frequent ENSO between ∼4.2 and 2.1 cal ka BP. After a possible depositional hiatus between ∼11.2 and 10.2 cal ka BP, the moderate-to-extremely altered sediments were deposited with enhanced erosion/runoff during the wetter ∼10.2–6 cal ka BP and the organic matter was dominantly autochthonous. Comparison of δ13Corg and C/N suggested that the contribution of C4 plants to organic matter increased over the drier ∼6–4.2 cal ka BP. Sediments representing this middle- Holocene drought-like condition showed geochemical similarity with sediments of the Post-Classic drought (∼1–0.4 cal ka BP), coeval with abandonment of the Xochicalco culture. Variation in seasonal insolation at orbital scales might have forced frequent droughts between ∼6 and 4.2 cal ka BP and the ENSO related short-lived storms possibly led to an unstable hydroclimate after ∼4.2 cal ka BP when the first Olmec settlements commenced in the region. Dissimilarity between this lacustrine archive and speleothems from southwest Mexico for the later part of the Holocene reflected different sensitivities of the geological records to seasonal and-annual precipitation.

Pollen-Based Holocene Thawing-History of Permafrost in Northern Asia and Its Potential Impacts on Climate Change

As the recent permafrost thawing of northern Asia proceeds due to anthropogenic climate change, precise and detailed palaeoecological records from past warm periods are essential to anticipate the extent of future permafrost variations. Here, based on the modern relationship between permafrost and vegetation (represented by pollen assemblages), we trained a Random Forest model using pollen and permafrost data and verified its reliability to reconstruct the history of permafrost in northern Asia during the Holocene. An early Holocene (12–8 cal ka BP) strong thawing trend, a middle-to-late Holocene (8–2 cal ka BP) relatively slow thawing trend, and a late Holocene freezing trend of permafrost in northern Asia are consistent with climatic proxies such as summer solar radiation and Northern Hemisphere temperature. The extensive distribution of permafrost in northern Asia inhibited the spread of evergreen coniferous trees during the early Holocene warming and might have decelerated the enhancement of the East Asian summer monsoon (EASM) by altering hydrological processes and albedo. Based on these findings, we suggest that studies of the EASM should consider more the state of permafrost and vegetation in northern Asia, which are often overlooked and may have a profound impact on climate change in this region.

Maritime-oriented foragers during the Late Pleistocene on the eastern costa del sol (Southeast Iberia): Cueva Victoria (Málaga, Spain)

The Mediterranean coast of Spain is marked by several clusters of Palaeolithic sites: to the south of the Pyrenees, in the area around the Ebro River, in the central part, and on the south coast, one of the southernmost regions in Europe. The number of sites is small compared with northern Iberia, but like that region, the Palaeolithic occupations are accompanied by several rock art ensembles. The archaeological material (both biotic and abiotic resources) and radiocarbon dates presented here were obtained during archaeological fieldwork of professor J. Fortea in the Late Pleistocene deposits in Cueva Victoria, located near the modern coastline and about 150 km north of the Strait of Gibraltar. In the three occupation phases, marine resources were acquired by shell-fishing (focusing almost exclusively on the clam Ruditapes decussatus), fishing, and the use of beached marine mammals. This contrasts with the limited data about the exploitation of terrestrial resources by hunting and gathering animals and plants. The study is completed by the study of artefacts (lithic and bone industry and objects of adornment) that help to understand the subsistence strategies of the cave occupants and enable a comparison with other groups inhabiting the Mediterranean coasts of the Iberian Peninsula during Greenland Interstadial 1, between ca. 15.1 and 13.6 cal BP.

Higher sea surface temperature in the Indian Ocean during the Last Interglacial weakened the South Asian monsoon

Understanding the drivers of South Asian monsoon intensity is pivotal for improving climate forecasting under global warming scenarios. Solar insolation is assumed to be the dominant driver of monsoon variability in warm climate regimes, but this has not been verified by proxy data. We report a South Asian monsoon rainfall record spanning the last ∼130 kyr in the Ganges–Brahmaputra–Meghna river catchment. Our multiproxy data reveal that the South Asian monsoon was weaker during the Last Interglacial (130 to 115 ka)—despite higher insolation—than during the Holocene (11.6 ka to present), thus questioning the widely accepted model assumption. Our work implies that Indian Ocean warming may increase the occurrence of severe monsoon failures in South Asia.

Addressing and anticipating future South Asian monsoon changes under continuing global warming is of critical importance for the food security and socioeconomic well-being of one-quarter of the world’s population. However, climate model projections show discrepancies in future monsoon variability in South Asian monsoon domains, largely due to our still limited understanding of the monsoon response to warm climate change scenarios. Particularly, climate models are largely based on the assumption that higher solar insolation causes higher rainfall during similar warm climatic regimes, but this has not been verified by proxy data for different interglacial periods. Here, we compare Indian summer monsoon (ISM) variability during the Last Interglacial and Holocene using a sedimentary leaf wax δD and δ¹³C record from the northern Bay of Bengal, representing the Ganges–Brahmaputra–Meghna (G-B-M) river catchment. In combination with a seawater salinity record, our results show that ISM intensity broadly follows summer insolation on orbital scales, but ISM intensity during the Last Interglacial was lower than during the Holocene despite higher summer insolation and greenhouse gas concentrations. We argue that sustained warmer sea surface temperature in the equatorial and tropical Indian Ocean during the Last Interglacial increased convective rainfall above the ocean but dampened ISM intensity on land. Our study demonstrates that besides solar insolation, internal climatic feedbacks also play an important role for South Asian monsoon variability during warm climate states. This work can help to improve future climate model projections and highlights the importance of understanding controls of monsoonal rainfall under interglacial boundary conditions.

A Jurassic record encodes an analogous Dansgaard-Oeschger climate periodicity

Earth’s past climate exhibits short-term (1500-year) pronounced fluctuations during the last glacial period, called Dansgaard–Oeschger (DO) glacial events, which have never been detected in pre-Quaternary times. The record of DO equivalent climate variability in Mesozoic strata can provide constraints on understanding these events. Here we highlight a prominent 1500-year cyclicity in a Jurassic (~ 155 Ma) ice-free sedimentary record from the Tethyan Basin. This Jurassic 1500-year cyclicity is encoded in high-resolution magnetic susceptibility (MS) proxy data reflecting detrital variations, and expressed as marl-limestone couplets. Additionally, MS data detect the modulation of these DO-scale couplets by supercouplet sets, reflecting the precession and its harmonics. We suggest that this Jurassic DO-like cyclicity may originate from paleo-monsoon-like system, analogous to the record of DO events in the Pleistocene East Asian monsoon archives. Paleogeographic reconstructions and atmosphere–ocean simulations further support the potential existence of strong, ancient monsoon circulations in the Tethyan Basin during the Jurassic.

Impact of paleoclimate on present and future evolution of the Greenland Ice Sheet

Using transient climate forcing based on simulations from the Alfred Wegener Institute Earth System Model (AWI-ESM), we simulate the evolution of the Greenland Ice Sheet (GrIS) from the last interglacial (125 ka, kiloyear before present) to 2100 AD with the Parallel Ice Sheet Model (PISM). The impact of paleoclimate, especially Holocene climate, on the present and future evolution of the GrIS is explored. Our simulations of the past show close agreement with reconstructions with respect to the recent timing of the peaks in ice volume and the climate of Greenland. The maximum and minimum ice volume at around 18-17 ka and 6-5 ka lag the respective extremes in climate by several thousand years, implying that the ice volume response of the GrIS strongly lags climatic changes. Given that Greenland's climate was getting colder from the Holocene Thermal Maximum (i.e., 8 ka) to the Pre-Industrial era, our simulation implies that the GrIS experienced growth from the mid-Holocene to the industrial era. Due to this background trend, the GrIS still gains mass until the second half of the 20th century, even though anthropogenic warming begins around 1850 AD. This is also in agreement with observational evidence showing mass loss of the GrIS does not begin earlier than the late 20th century. Our results highlight that the present evolution of the GrIS is not only controlled by the recent climate changes, but is also affected by paleoclimate, especially the relatively warm Holocene climate. We propose that the GrIS was not in equilibrium throughout the entire Holocene and that the slow response to Holocene climate needs to be represented in ice sheet simulations in order to predict ice mass loss, and therefore sea level rise, accurately.

Climate changes modulated the history of Arctic iodine during the Last Glacial Cycle

Iodine has a significant impact on promoting the formation of new ultrafine aerosol particles and accelerating tropospheric ozone loss, thereby affecting radiative forcing and climate. Therefore, understanding the long-term natural evolution of iodine, and its coupling with climate variability, is key to adequately assess its effect on climate on centennial to millennial timescales. Here, using two Greenland ice cores (NEEM and RECAP), we report the Arctic iodine variability during the last 127,000 years. We find the highest and lowest iodine levels recorded during interglacial and glacial periods, respectively, modulated by ocean bioproductivity and sea ice dynamics. Our sub-decadal resolution measurements reveal that high frequency iodine emission variability occurred in pace with Dansgaard/Oeschger events, highlighting the rapid Arctic ocean-ice-atmosphere iodine exchange response to abrupt climate changes. Finally, we discuss if iodine levels during past warmer-than-present climate phases can serve as analogues of future scenarios under an expected ice-free Arctic Ocean. We argue that the combination of natural biogenic ocean iodine release (boosted by ongoing Arctic warming and sea ice retreat) and anthropogenic ozone-induced iodine emissions may lead to a near future scenario with the highest iodine levels of the last 127,000 years.

Late Quaternary dynamics of Arctic biota from ancient environmental genomics

During the last glacial-interglacial cycle, Arctic biotas experienced substantial climatic changes, yet the nature, extent and rate of their responses are not fully understood1-8. Here we report a large-scale environmental DNA metagenomic study of ancient plant and mammal communities, analysing 535 permafrost and lake sediment samples from across the Arctic spanning the past 50,000 years. Furthermore, we present 1,541 contemporary plant genome assemblies that were generated as reference sequences. Our study provides several insights into the long-term dynamics of the Arctic biota at the circumpolar and regional scales. Our key findings include: (1) a relatively homogeneous steppe-tundra flora dominated the Arctic during the Last Glacial Maximum, followed by regional divergence of vegetation during the Holocene epoch; (2) certain grazing animals consistently co-occurred in space and time; (3) humans appear to have been a minor factor in driving animal distributions; (4) higher effective precipitation, as well as an increase in the proportion of wetland plants, show negative effects on animal diversity; (5) the persistence of the steppe-tundra vegetation in northern Siberia enabled the late survival of several now-extinct megafauna species, including the woolly mammoth until 3.9 ± 0.2 thousand years ago (ka) and the woolly rhinoceros until 9.8 ± 0.2 ka; and (6) phylogenetic analysis of mammoth environmental DNA reveals a previously unsampled mitochondrial lineage. Our findings highlight the power of ancient environmental metagenomics analyses to advance understanding of population histories and long-term ecological dynamics.

Globally resolved surface temperatures since the Last Glacial Maximum

Climate changes across the past 24,000 years provide key insights into Earth system responses to external forcing. Climate model simulations^1,2 and proxy data^3-8 have independently allowed for study of this crucial interval; however, they have at times yielded disparate conclusions. Here, we leverage both types of information using paleoclimate data assimilation^9,10 to produce the first proxy-constrained, full-field reanalysis of surface temperature change spanning the Last Glacial Maximum to present at 200-year resolution. We demonstrate that temperature variability across the past 24 thousand years was linked to two primary climatic mechanisms: radiative forcing from ice sheets and greenhouse gases; and a superposition of changes in the ocean overturning circulation and seasonal insolation. In contrast with previous proxy-based reconstructions^6,7 our results show that global mean temperature has slightly but steadily warmed, by ~0.5 °C, since the early Holocene (around 9 thousand years ago). When compared with recent temperature changes¹¹, our reanalysis indicates that both the rate and magnitude of modern warming are unusual relative to the changes of the past 24 thousand years.

Roles of insolation forcing and CO2 forcing on Late Pleistocene seasonal sea surface temperatures

Late Pleistocene changes in insolation, greenhouse gas concentrations, and ice sheets have different spatially and seasonally modulated climatic fingerprints. By exploring the seasonality of paleoclimate proxy data, we gain deeper insight into the drivers of climate changes. Here, we investigate changes in alkenone-based annual mean and Globigerinoides ruber Mg/Ca-based summer sea surface temperatures in the East China Sea and their linkages to climate forcing over the past 400,000 years. During interglacial-glacial cycles, there are phase differences between annual mean and seasonal (summer and winter) temperatures, which relate to seasonal insolation changes. These phase differences are most evident during interglacials. During glacial terminations, temperature changes were strongly affected by CO₂. Early temperature minima, ~20,000 years before glacial terminations, except the last glacial period, coincide with the largest temperature differences between summer and winter, and with the timing of the lowest atmospheric CO₂ concentration. These findings imply the need to consider proxy seasonality and seasonal climate variability to estimate climate sensitivity.

How temperatures at different seasons differ in response to different forcings is not well known. Here, the authors reconstruct annual and seasonal sea surface temperatures in the East China Sea and show that they react differently to CO2 and insolation forcing on glacial-interglacial timescales.

Adaptive capacity and flexibility of the Neanderthals at Heidenschmiede (Swabian Jura) with regard to core reduction strategies

The branched reduction system at the Heidenschmiede described here is hitherto exceptional for the Middle Paleolithic of the Swabian Jura. By means of refits and supporting objects, we are able to describe a superordinate reduction system that combines several individual reduction concepts, such as Levallois and blade production, within one volume. In the Middle Paleolithic of the Swabian Jura, blade technology has thus far played a rather minor role. On the one hand, it is possible to split a selected volume (nodule) into three parts, which are reduced separately according to individual concepts. On the other hand, it is also possible to reduce parts of a volume with one concept first and then with another. The hypothetical reduction system can be branched or linear, thus emphasizing the technological flexibility in core reduction, which requires a high degree of cognitive skills of three-dimensional imagination.

Insolation triggered abrupt weakening of Atlantic circulation at the end of interglacials

Abrupt cooling is observed at the end of interglacials in many paleoclimate records, but the mechanism responsible remains unclear. Using model simulations, we demonstrate that there exists a threshold in the level of astronomically induced insolation below which abrupt changes at the end of interglacials of the past 800,000 years occur. When decreasing insolation reaches the critical value, it triggers a strong, abrupt weakening of the Atlantic meridional overturning circulation and a cooler mean climate state accompanied by high-amplitude variations lasting for several thousand years. The mechanism involves sea ice feedbacks in the Nordic and Labrador Seas. The ubiquity of this threshold suggests its fundamental role in terminating the warm climate conditions at the end of interglacials.

Speleothem Records of the Hydroclimate Variability throughout the Last Glacial Cycle from Manita peć Cave (Velebit Mountain, Croatia)

We present stable carbon (δ13C) and oxygen (δ18O) isotope records from two partially coeval speleothems from Manita peć Cave, Croatia. The cave is located close to the Adriatic coast (3.7 km) at an elevation of 570 m a.s.l. The site experienced competing Mediterranean and continental climate influences throughout the last glacial cycle and was situated close to the ice limit during the glacial phases. U-Th dating constrains the growth history from Marine Isotope Stage (MIS) 5 to MIS 3 and the transition from MIS 2 to MIS 1. 14C dating was used to estimate the age of the youngest part of one stalagmite found to be rich in detrital thorium and thus undatable by U-Th. On a millennial scale, δ18O variations partly mimic the Dansgaard–Oeschger interstadials recorded in Greenland ice cores (Greenland Interstadials, GI) from GI 22 to GI 13. We interpret our δ18O record as a proxy for variations in precipitation amount and/or moisture sources, and the δ13C record is interpreted as a proxy for changes in soil bioproductivity. The latter indicates a generally reduced vegetation cover towards MIS 3–MIS 4, with shifts of ~8‰ and approaching values close to those of the host rock. However, even during the coldest phases, when a periglacial setting and enhanced aridity sustained long-residence-time groundwater, carbonic-acid dissolution remains the driving force of the karstification processes. Speleothem morphology follows changes in environmental conditions and complements regional results of submerged speleothems findings. Specifically, narrow sections of light porous spelaean calcite precipitated during the glacial/stadial sea-level lowstands, while the warmer and wetter conditions were marked with compact calcite and hiatuses in submerged speleothems due to sea-level highstands. Presumably, the transformation of this littoral site to a continental one with somewhat higher amounts of orographic precipitation was a site-specific effect that masked regional environmental changes.

Ancient Faunal History Revealed by Interdisciplinary Biomolecular Approaches

Starting four decades ago, studies have examined the ecology and evolutionary dynamics of populations and species using short mitochondrial DNA fragments and stable isotopes. Through technological and analytical advances, the methods and biomolecules at our disposal have increased significantly to now include lipids, whole genomes, proteomes, and even epigenomes. At an unprecedented resolution, the study of ancient biomolecules has made it possible for us to disentangle the complex processes that shaped the ancient faunal diversity across millennia, with the potential to aid in implicating probable causes of species extinction and how humans impacted the genetics and ecology of wild and domestic species. However, even now, few studies explore interdisciplinary biomolecular approaches to reveal ancient faunal diversity dynamics in relation to environmental and anthropogenic impact. This review will approach how biomolecules have been implemented in a broad variety of topics and species, from the extinct Pleistocene megafauna to ancient wild and domestic stocks, as well as how their future use has the potential to offer an enhanced understanding of drivers of past faunal diversity on Earth.

Time reversal symmetry and the difference between relaxations and building-up periods

Autocorrelations in stationary systems are time-symmetric, irrespective of the signal's properties. Linear dynamics is usually associated with signals which are statistically time inversion invariant. This is known to be broken for non-Gaussian models. In this paper, we develop a theoretical framework of time reversibility of linear models based on noncontinuous driving. We identify the inverse decay exponent of the autocorrelation function as either a characteristic time for the building-up of extreme events in the time series or of relaxations after these extreme events. If the characteristic time is known, the dynamics can be inverted in both directions in time and the residuals can be compared, which gives a criterion for the type of time inversion asymmetry. The method is applied to two time series from atmospheric science with different behavior.

A review of current and emerging approaches for Quaternary marine sediment dating

Dating methodologies for Quaternary marine sediments play increasingly important roles in the reconstruction of paleoenvironments and paleoclimate in (paleo)oceanography. Previous reviews or studies have focused mainly on one or two methodologies, and their applications in one specific environment. With the continuing technological and methodological advances in different methods over the past few decades, an up-to-date comparison of the pros and cons of each dating methodology is needed to clearly understand their applications in marine geoscience research. In this review, we first briefly summarized the common methods of absolute dating and relative dating. These are (1) radioisotope dating with different half-lives using natural nuclides of ²³⁴Th, ²¹⁰Pb, ²³⁰Th, and ²²⁶Ra, cosmogenic nuclides of ⁷Be, ¹⁴C, ¹⁰Be, ³²Si, ²⁶Al, ³⁶Cl and ²¹Ne, and the artificial radionuclides of ¹³⁷Cs, ^{239, 240}Pu, ²⁴¹Am and ¹²⁹I that have been induced by atmospheric nuclear tests, accidents in nuclear plants, and discharges of radioactive wastes; (2) radiation exposure dating of luminescence and electron paramagnetic resonance (ESR) dating; and (3) stratigraphic dating of δ¹⁸O and paleomagnetic sequence. Applications and limitations from the marine terraces, estuaries, to hadal trenches have been summarized to each technique in the study of Quaternary marine geoscience extending from the Anthropocene through the Pleistocene. Finally, we introduced some emerging event dating methods, namely the arrivals of microplastics, mercury isotopes, and organic pollutant deposition that all appeared after the industrial resolution in our now changing ocean influenced by acidification, global warming, and anthropogenic activities. We ended by discussing future perspectives for reliable and high-resolution chronology by interdisciplinary methods including computer programming to better understand the natural geological evolution and predict the future changes in earth science.

Addressing alpine plant phylogeography using integrative distributional, demographic and coalescent modeling

Phylogeographic studies of alpine plants have evolved considerably in the last two decades from ad hoc interpretations of genetic data to statistical model-based approaches. In this review we outline the developments in alpine plant phylogeography focusing on the recent approach of integrative distributional, demographic and coalescent (iDDC) modeling. By integrating distributional data with spatially explicit demographic modeling and subsequent coalescent simulations, the history of alpine species can be inferred and long-standing hypotheses, such as species-specific responses to climate change or survival on nunataks during the last glacial maximum, can be efficiently tested as exemplified by available case studies. We also discuss future prospects and improvements of iDDC.

Sedimentary ancient DNA shows terrestrial plant richness continuously increased over the Holocene in northern Fennoscandia

The effects of climate change on species richness are debated but can be informed by the past. Here, we generated a sedimentary ancient DNA dataset covering 10 lakes and applied novel methods for data harmonization. We assessed the impact of Holocene climate changes and nutrients on terrestrial plant richness in northern Fennoscandia. We find that richness increased steeply during the rapidly warming Early Holocene. In contrast to findings from most pollen studies, we show that richness continued to increase thereafter, although the climate was stable, with richness and the regional species pool only stabilizing during the past three millennia. Furthermore, overall increases in richness were greater in catchments with higher soil nutrient availability. We suggest that richness will increase with ongoing warming, especially at localities with high nutrient availability and assuming that human activity remains low in the region, although lags of millennia may be expected.

Temporal dynamics of teleost populations during the Pleistocene: a report from publicly available genome data

Background

Global climate oscillation, as a selection dynamic, is an ecologically important element resulting in global biodiversity. During the glacial geological periods, most organisms suffered detrimental selection pressures (such as food shortage and habitat loss) and went through population declines. However, during the mild interglacial periods, many species re-flourished. These temporal dynamics of effective population sizes (N_e) provide essential information for understanding and predicting evolutionary outcomes during historical and ongoing global climate changes.

Results

Using high-quality genome assemblies and corresponding sequencing data, we applied the Pairwise Sequentially Markovian Coalescent (PSMC) method to quantify N_e changes of twelve representative teleost species from approximately 10 million years ago (mya) to 10 thousand years ago (kya). These results revealed multiple rounds of population contraction and expansion in most of the examined teleost species during the Neogene and the Quaternary periods. We observed that 83% (10/12) of the examined teleosts had experienced a drastic decline in N_e before the last glacial period (LGP, 110–12 kya), slightly earlier than the reported pattern of N_e changes in 38 avian species. In comparison with the peaks, almost all of the examined teleosts maintained long-term lower N_e values during the last few million years. This is consistent with increasingly dramatic glaciation during this period.

Conclusion

In summary, these findings provide a more comprehensive understanding of the historical N_e changes in teleosts. Results presented here could lead to the development of appropriate strategies to protect species in light of ongoing global climate changes.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07816-7.