Abstract: The Piacenzian stage of the Pliocene (2.6 to 3.6 Ma) is the most recent past interval of sustained global warmth with mean global temperatures markedly higher (by ~2–3 °C) than today. Quantifying CO2 levels during the mid-Piacenzian Warm Period (mPWP) provides a means, therefore, to deepen our understanding of Earth System behaviour in a warm climate state. Here we present a new high-resolution record of atmospheric CO2 using the δ11B-pH proxy from 3.35 to 3.15 million years ago (Ma) at a temporal resolution of 1 sample per 3–6 thousand years (kyrs). Our study interval covers both the coolest marine isotope stage of the mPWP, M2 (~3.3 Ma) and the transition into its warmest phase including interglacial KM5c (centered on ~3.205 Ma) which has a similar orbital configuration to present. We find that CO2 ranged from $${389}_{-8}^{+38}$$ppm to $${331}_{-11}^{+13},$$ppm, with CO2 during the KM5c interglacial being $${371}_{-29}^{+32}\,$$ppm (at 95% confidence). Our findings corroborate the idea that changes in atmospheric CO2 levels played a distinct role in climate variability during the mPWP. They also facilitate ongoing data-model comparisons and suggest that, at present rates of human emissions, there will be more CO2 in Earth’s atmosphere by 2025 than at any time in at least the last 3.3 million years.

Abstract: Marine food-reliant subsistence systems such as those in the African Middle Stone Age (MSA) were not thought to exist in Europe until the much later Mesolithic. Whether this apparent lag reflects taphonomic biases or behavioral distinctions between archaic and modern humans remains much debated. Figueira Brava cave, in the Arrabida range (Portugal), provides an exceptionally well preserved record of Neandertal coastal resource exploitation on a comparable scale to the MSA and dated to ~86 to 106 thousand years ago. The breadth of the subsistence base-pine nuts, marine invertebrates, fish, marine birds and mammals, tortoises, waterfowl, and hoofed game-exceeds that of regional early Holocene sites. Fisher-hunter-gatherer economies are not the preserve of anatomically modern people; by the Last Interglacial, they were in place across the Old World in the appropriate settings.

Abstract: The Last Interglacial (LIG), a warmer period 130,000–116,000 years before present, is a potential analogue for future climate change. Stronger LIG summertime insolation at high northern latitudes drove Arctic land summer temperatures 4–5 °C higher than in the pre-industrial era. Climate model simulations have previously failed to capture these elevated temperatures, possibly because they were unable to correctly capture LIG sea-ice changes. Here, we show that the latest version of the fully coupled UK Hadley Center climate model (HadGEM3) simulates a more accurate Arctic LIG climate, including elevated temperatures. Improved model physics, including a sophisticated sea-ice melt-pond scheme, result in a complete simulated loss of Arctic sea ice in summer during the LIG, which has yet to be simulated in past generations of models. This ice-free Arctic yields a compelling solution to the long-standing puzzle of what drove LIG Arctic warmth and supports a fast retreat of future Arctic summer sea ice. Arctic climate in the Last Interglacial (LIG)—a warm period 130,000–116,000 years ago—is poorly simulated by modern climate models. A model with improved sea-ice melt-pond physics reproduces LIG Arctic temperatures, suggests an ice-free Arctic during this period and predicts the same by 2035.

Abstract: Previous studies have suggested that during the late Pleistocene ice ages, surface-deep exchange was somehow weakened in the Southern Ocean's Antarctic Zone, which reduced the leakage of deeply sequestered carbon dioxide and thus contributed to the lower atmospheric carbon dioxide levels of the ice ages. Here, high-resolution diatom-bound nitrogen isotope measurements from the Indian sector of the Antarctic Zone reveal three modes of change in Southern Westerly Wind-driven upwelling, each affecting atmospheric carbon dioxide. Two modes, related to global climate and the bipolar seesaw, have been proposed previously. The third mode-which arises from the meridional temperature gradient as affected by Earth's obliquity (axial tilt)-can explain the lag of atmospheric carbon dioxide behind climate during glacial inception and deglaciation. This obliquity-induced lag, in turn, makes carbon dioxide a delayed climate amplifier in the late Pleistocene glacial cycles.

Abstract: A climate/vegetation model simulates episodic wetter and drier periods at the 21,000-y precession period in eastern North Africa, the Arabian Peninsula, and the Levant over the past 140,000 y. Large orbitally forced wet/dry extremes occur during interglacial time, ∼130 to 80 ka, and conditions between these two extremes prevail during glacial time, ∼70 to 15 ka. Orbital precession causes high seasonality in Northern Hemisphere (NH) insolation at ∼125, 105, and 83 ka, with stronger and northward extended summer monsoon rains in North Africa and the Arabian Peninsula and increased winter rains in the Mediterranean Basin. The combined effects of these two seasonally distinct rainfall regimes increase vegetation and narrow the width of the Saharan-Arabian desert and semidesert zones. During the opposite phase of the precession cycle (∼115, 95, and 73 ka), NH seasonality is low, and decreased summer insolation and increased winter insolation cause monsoon and storm track rains to decrease and the width of the desert zone to increase. During glacial time (∼70 to 15 ka), forcing from large ice sheets and lowered greenhouse gas concentrations combine to increase winter Mediterranean storm track precipitation; the southward retreat of the northern limit of summer monsoon rains is relatively small, thereby limiting the expansion of deserts. The lowered greenhouse gas concentrations cause the near-equatorial zone to cool and reduce convection, causing drier climate with reduced forest cover. At most locations and times, the simulations agree with environmental observations. These changing regional patterns of climate/vegetation could have influenced the dispersal of early humans through expansions and contractions of well-watered corridors.

Abstract: Southeast Asia has emerged as an important region for understanding hominin and mammalian migrations and extinctions. High-profile discoveries have shown that Southeast Asia has been home to at least five members of the genus Homo1–3. Considerable turnover in Pleistocene megafauna has previously been linked with these hominins or with climate change4, although the region is often left out of discussions of megafauna extinctions. In the traditional hominin evolutionary core of Africa, attempts to establish the environmental context of hominin evolution and its association with faunal changes have long been informed by stable isotope methodologies5,6. However, such studies have largely been neglected in Southeast Asia. Here we present a large-scale dataset of stable isotope data for Southeast Asian mammals that spans the Quaternary period. Our results demonstrate that the forests of the Early Pleistocene had given way to savannahs by the Middle Pleistocene, which led to the spread of grazers and extinction of browsers—although geochronological limitations mean that not all samples can be resolved to glacial or interglacial periods. Savannahs retreated by the Late Pleistocene and had completely disappeared by the Holocene epoch, when they were replaced by highly stratified closed-canopy rainforest. This resulted in the ascendency of rainforest-adapted species as well as Homo sapiens—which has a unique adaptive plasticity among hominins—at the expense of savannah and woodland specialists, including Homo erectus. At present, megafauna are restricted to rainforests and are severely threatened by anthropogenic deforestation. Stable isotope data for Southeast Asian mammals across the Quaternary period shed light on environmental change from the Early Pleistocene to the Holocene epoch, contextualizing hominin evolution and megafauna extinction in the region.

Abstract: A range of future climate scenarios are projected for high atmospheric CO2 concentrations, given uncertainties over future human actions as well as potential environmental and climatic feedbacks. The geological record offers an opportunity to understand climate system response to a range of forcings and feedbacks which operate over multiple temporal and spatial scales. Here, we examine a single interglacial during the late Pliocene (KM5c, ca. 3.205±0.01 Ma) when atmospheric CO2 exceeded pre-industrial concentrations, but were similar to today and to the lowest emission scenarios for this century. As orbital forcing and continental configurations were almost identical to today, we are able to focus on equilibrium climate system response to modern and near-future CO2. Using proxy data from 32 sites, we demonstrate that global mean sea-surface temperatures were warmer than pre-industrial values, by ∼2.3 ∘C for the combined proxy data (foraminifera Mg∕Ca and alkenones), or by ∼3.2–3.4 ∘C (alkenones only). Compared to the pre-industrial period, reduced meridional gradients and enhanced warming in the North Atlantic are consistently reconstructed. There is broad agreement between data and models at the global scale, with regional differences reflecting ocean circulation and/or proxy signals. An uneven distribution of proxy data in time and space does, however, add uncertainty to our anomaly calculations. The reconstructed global mean sea-surface temperature anomaly for KM5c is warmer than all but three of the PlioMIP2 model outputs, and the reconstructed North Atlantic data tend to align with the warmest KM5c model values. Our results demonstrate that even under low-CO2 emission scenarios, surface ocean warming may be expected to exceed model projections and will be accentuated in the higher latitudes.

Abstract: We investigate hydrology during a past climate slightly warmer than the present: the last interglacial (LIG). With daily output of preindustrial and LIG simulations from eight new climate models we force hydrological model PCR-GLOBWB and in turn hydrodynamic model CaMa-Flood. Compared to preindustrial, annual mean LIG runoff, discharge, and 100-yr flood volume are considerably larger in the Northern Hemisphere, by 14%, 25%, and 82%, respectively. Anomalies are negative in the Southern Hemisphere. In some boreal regions, LIG runoff and discharge are lower despite higher precipitation, due to the higher temperatures and evaporation. LIG discharge is much higher for the Niger, Congo, Nile, Ganges, Irrawaddy, and Pearl and lower for the Mississippi, Saint Lawrence, Amazon, Parana, Orange, Zambesi, Danube, and Ob. Discharge is seasonally postponed in tropical rivers affected by monsoon changes. Results agree with published proxies on the sign of discharge anomaly in 15 of 23 sites where comparison is possible.

Abstract: While Asian monsoon (AM) changes have been clearly captured in Chinese speleothem oxygen isotope (δ18O) records, the lack of glacial-interglacial variability in the records remains puzzling. Here, we report speleothem δ18O records from three locations along the trajectory of the Indian summer monsoon (ISM), a major branch of the AM, and characterize AM rainfall over the past 180,000 years. We have found that the records close to the monsoon moisture source show large glacial-interglacial variability, which then decreases landward. These changes likely reflect a stronger oxygen isotope fractionation associated with progressive rainout of AM moisture during glacial periods, possibly due to a larger temperature gradient and suppressed plant transpiration. We term this effect, which counteracts the forcing of glacial boundary conditions, the moisture transport pathway effect.

Abstract: Pulse-like carbon dioxide release to the atmosphere on centennial time scales has only been identified for the most recent glacial and deglacial periods and is thought to be absent during warmer climate conditions. Here, we present a high-resolution carbon dioxide record from 330,000 to 450,000 years before present, revealing pronounced carbon dioxide jumps (CDJ) under cold and warm climate conditions. CDJ come in two varieties that we attribute to invigoration or weakening of the Atlantic meridional overturning circulation (AMOC) and associated northward and southward shifts of the intertropical convergence zone, respectively. We find that CDJ are pervasive features of the carbon cycle that can occur during interglacial climate conditions if land ice masses are sufficiently extended to be able to disturb the AMOC by freshwater input.

Abstract: The end of the Pleistocene in North America saw the extinction of 38 genera of mostly large mammals. As their disappearance seemingly coincided with the arrival of people in the Americas, their extinction is often attributed to human overkill, notwithstanding a dearth of archaeological evidence of human predation. Moreover, this period saw the extinction of other species, along with significant changes in many surviving taxa, suggesting a broader cause, notably, the ecological upheaval that occurred as Earth shifted from a glacial to an interglacial climate. But, overkill advocates ask, if extinctions were due to climate changes, why did these large mammals survive previous glacial-interglacial transitions, only to vanish at the one when human hunters were present? This question rests on two assumptions: that previous glacial-interglacial transitions were similar to the end of the Pleistocene, and that the large mammal genera survived unchanged over multiple such cycles. Neither is demonstrably correct. Resolving the cause of large mammal extinctions requires greater knowledge of individual species' histories and their adaptive tolerances, a fuller understanding of how past climatic and ecological changes impacted those animals and their biotic communities, and what changes occurred at the Pleistocene-Holocene boundary that might have led to those genera going extinct at that time. Then we will be able to ascertain whether the sole ecologically significant difference between previous glacial-interglacial transitions and the very last one was a human presence.

Abstract: Disrupting North Atlantic Deep Water (NADW) ventilation is a key concern in climate projections. We use (sub)centennially resolved bottom water δ13C records that span the interglacials of the last 0.5 million years to assess the frequency of and the climatic backgrounds capable of triggering large NADW reductions. Episodes of reduced NADW in the deep Atlantic, similar in magnitude to glacial events, have been relatively common and occasionally long-lasting features of interglacials. NADW reductions were triggered across the range of recent interglacial climate backgrounds, which demonstrates that catastrophic freshwater outburst floods were not a prerequisite for large perturbations. Our results argue that large NADW disruptions are more easily achieved than previously appreciated and that they occurred in past climate conditions similar to those we may soon face.

Abstract: Efforts to improve sea level forecasting on a warming planet have focused on determining the temperature, sea level and extent of polar ice sheets during Earth’s past interglacial warm periods1–3. About 400,000 years ago, during the interglacial period known as Marine Isotopic Stage 11 (MIS11), the global temperature was 1 to 2 degrees Celsius greater2 and sea level was 6 to 13 metres higher1,3. Sea level estimates in excess of about 10 metres, however, have been discounted because these require a contribution from the East Antarctic Ice Sheet3, which has been argued to have remained stable for millions of years before and includes MIS114,5. Here we show how the evolution of 234U enrichment within the subglacial waters of East Antarctica recorded the ice sheet’s response to MIS11 warming. Within the Wilkes Basin, subglacial chemical precipitates of opal and calcite record accumulation of 234U (the product of rock–water contact within an isolated subglacial reservoir) up to 20 times higher than that found in marine waters. The timescales of 234U enrichment place the inception of this reservoir at MIS11. Informed by the 234U cycling observed in the Laurentide Ice Sheet, where 234U accumulated during periods of ice stability6 and was flushed to global oceans in response to deglaciation7, we interpret our East Antarctic dataset to represent ice loss within the Wilkes Basin at MIS11. The 234U accumulation within the Wilkes Basin is also observed in the McMurdo Dry Valleys brines8–10, indicating11 that the brine originated beneath the adjacent East Antarctic Ice Sheet. The marine origin of brine salts10 and bacteria12 implies that MIS11 ice loss was coupled with marine flooding. Collectively, these data indicate that during one of the warmest Pleistocene interglacials, the ice sheet margin at the Wilkes Basin retreated to near the precipitate location, about 700 kilometres inland from the current position of the ice margin, which—assuming current ice volumes—would have contributed about 3 to 4 metres13 to global sea levels. Uranium isotopes in subglacial precipitates from the Wilkes Basin of the East Antarctic Ice Sheet reveal ice retreat during a warm Pleistocene interglacial period about 400,000 years ago.

Abstract: The nature of human dispersals out of Africa has remained elusive because of the poor resolution of paleoecological data in direct association with remains of the earliest non-African people. Here, we report hominin and non-hominin mammalian tracks from an ancient lake deposit in the Arabian Peninsula, dated within the last interglacial. The findings, it is argued, likely represent the oldest securely dated evidence for Homo sapiens in Arabia. The paleoecological evidence indicates a well-watered semi-arid grassland setting during human movements into the Nefud Desert of Saudi Arabia. We conclude that visitation to the lake was transient, likely serving as a place to drink and to forage, and that late Pleistocene human and mammalian migrations and landscape use patterns in Arabia were inexorably linked.

Abstract: Southerly wind in the lower troposphere is an essential feature of East Asian summer monsoon (EASM) circulation, which is reported to be enhanced under global warming scenarios and interglacial epochs. Based on an analysis of an ensemble of CMIP6 models, this study shows that the magnitude of intensification of the EASM circulation is much smaller under global warming scenarios than during interglacial epochs. Distinct changes in the western North Pacific subtropical high (WNPSH) are responsible for the different responses of the EASM circulation. The WNPSH is substantially enhanced during interglacial epochs, which acts to strengthen the southerly wind associated with the EASM on the western flank of the WNPSH. However, the change in the WNPSH is insignificant and cannot strengthen the EASM under global warming scenarios, and the weakly enhanced EASM circulation may be a direct response to intensified heating over the Tibetan Plateau. The land–ocean thermal contrast explains the different responses of the WNPSH. During interglacial epochs, the summertime surface warming over the subtropical North Pacific is much weaker than over Eurasia due to the large thermal inertia of the ocean to increased insolation, and the WNPSH is intensified as a response to the suppressed latent heating over the subtropical North Pacific. The fast response of the WNPSH to abrupt quadrupling of CO2 without sufficient ocean warming is an analog to the interglacial epochs, but it is offset by the effect of slow oceanic warming, resulting in an insignificant change of the WNPSH under global warming scenarios.

Abstract: The impact of the ongoing anthropogenic warming on the Arctic Ocean sea ice is ascertained and closely monitored. However, its long-term fate remains an open question as its natural variability on centennial to millennial timescales is not well documented. Here, we use marine sedimentary records to reconstruct Arctic sea-ice fluctuations. Cores collected along the Lomonosov Ridge that extends across the Arctic Ocean from northern Greenland to the Laptev Sea were radiocarbon dated and analyzed for their micropaleontological and palynological contents, both bearing information on the past sea-ice cover. Results demonstrate that multiyear pack ice remained a robust feature of the western and central Lomonosov Ridge and that perennial sea ice remained present throughout the present interglacial, even during the climate optimum of the middle Holocene that globally peaked ∼6,500 y ago. In contradistinction, the southeastern Lomonosov Ridge area experienced seasonally sea-ice-free conditions, at least, sporadically, until about 4,000 y ago. They were marked by relatively high phytoplanktonic productivity and organic carbon fluxes at the seafloor resulting in low biogenic carbonate preservation. These results point to contrasted west–east surface ocean conditions in the Arctic Ocean, not unlike those of the Arctic dipole linked to the recent loss of Arctic sea ice. Hence, our data suggest that seasonally ice-free conditions in the southeastern Arctic Ocean with a dominant Arctic dipolar pattern, may be a recurrent feature under “warm world” climate.

Abstract: The Antarctic Circumpolar Current (ACC) is the world's largest current system connecting all three major basins of the global ocean. Our knowledge of glacial‐interglacial changes in ACC dynamics in the southeast Pacific is not well constrained and presently only based on reconstructions covering the last glacial cycle. Here we use a combination of mean sortable silt grain size of the terrigenous sediment fraction (10–63 μm, "Sortable Silt") and X‐ray fluorescence scanner‐derived Zr/Rb ratios as flow strength proxies to examine ACC variations at the Pacific entrance to the Drake Passage (DP) in the vicinity of the Subantarctic Front. Our results indicate that at the DP entrance, ACC strength varied by ~6–16% on glacial‐interglacial time scales, yielding higher current speeds during interglacial times and reduced current speeds during glacials. We provide evidence that previous observations of a reduction in DP throughflow during the last glacial period are part of a consistent pattern extending for at least the last 1.3 Ma. The orbital‐scale cyclicity follows well‐known global climate changes from prevailing ca. 41‐kyr cycles in the early part of the record (1.3 Ma to 850 ka; marine isotope stage 21) across the mid‐Pleistocene transition into the middle and late Pleistocene 100‐kyr world. A comparison to a bottom water flow record from the deep western boundary current off New Zealand (Ocean Drilling Program Site 1123) reveals anti‐phased changes between the two sites. The enhanced supply of deep water along the DP and into the Atlantic Ocean during interglacials corresponds to a weakened flow of the SW Pacific deep western boundary current.

Abstract: . A valuable analogue for assessing Earth’s sensitivity to warming is the Last Interglacial (LIG; 129–116 kyr), when global temperatures (0−+2 °C) and mean sea level (+6–11 m) were higher than today. The direct contribution of warmer conditions to global sea level (thermosteric) are uncertain. We report here a global network of LIG sea surface temperatures (SST) obtained from various published temperature proxies (e.g. faunal/floral assemblages, Mg/Ca ratios of calcareous plankton, alkenone UK’37). Each reconstruction is averaged across the LIG (anomalies relative to 1981–2010), corrected for ocean drift and with varying seasonality (189 annual, 99 December-February, and 92 June–August records). We summarise the current limitations of SST reconstructions for the LIG and the spatial temperature features of a naturally warmer world. Because of local δ18O seawater changes, uncertainty in the age models of marine cores, and differences in sampling resolution and/or sedimentation rates, the reconstructions are restricted to mean conditions. To avoid bias towards individual LIG SSTs based on only a single (and potentially erroneous) measurement or a single interpolated data point, here we average across the entire LIG. To investigate the sensitivity of the reconstruction to high temperatures, we also report maximum values during the first 5 ka of the LIG (129–124 kyr). The global dataset provides a remarkably coherent pattern of higher SST increases at polar latitudes than in the tropics, with comparable estimates between different SST proxies. We report mean global annual SST anomalies of 0.2 ± 0.1 °C and a maximum of 0.9 ± 0.2 °C respectively. Using the reconstructed SSTs suggests a mean thermosteric sea level rise of 0.01 ± 0.1 m and a maximum of 0.13 ± 0.1 m respectively. The data provide an important natural baseline for a warmer world, constraining the contributions of Greenland and Antarctic ice sheets to global sea level during a geographically widespread expression of high sea level, and can be used to test the next inter-comparison of models for projecting future climate change. The dataset described in this paper, including summary temperature and thermosteric sea-level reconstructions, are available at https://doi.org/10.1594/PANGAEA.904381 (Turney et al., 2019).