Highly restricted near-surface permafrost extent during the mid-Pliocene warm period.

Accurate understanding of permafrost dynamics is critical for evaluating and mitigating impacts that may arise as permafrost degrades in the future; however, existing projections have large uncertainties. Studies of how permafrost responded historically during Earth's past warm periods are helpful in exploring potential future permafrost behavior and to evaluate the uncertainty of future permafrost change projections. Here, we combine a surface frost index model with outputs from the second phase of the Pliocene Model Intercomparison Project to simulate the near-surface (~3 to 4 m depth) permafrost state in the Northern Hemisphere during the mid-Pliocene warm period (mPWP, ~3.264 to 3.025 Ma). This period shares similarities with the projected future climate. Constrained by proxy-based surface air temperature records, our simulations demonstrate that near-surface permafrost was highly spatially restricted during the mPWP and was 93 ± 3% smaller than the preindustrial extent. Near-surface permafrost was present only in the eastern Siberian uplands, Canadian high Arctic Archipelago, and northernmost Greenland. The simulations are similar to near-surface permafrost changes projected for the end of this century under the SSP5-8.5 scenario and provide a perspective on the potential permafrost behavior that may be expected in a warmer world.

Role of eruption season in reconciling model and proxy responses to tropical volcanism.

The response of the El Niño/Southern Oscillation (ENSO) to tropical volcanic eruptions has important worldwide implications, but remains poorly constrained. Paleoclimate records suggest an "El Niño-like" warming 1 year following major eruptions [Adams JB, Mann ME, Ammann CM (2003) Nature 426:274-278] and "La Niña-like" cooling within the eruption year [Li J, et al. (2013) Nat Clim Chang 3:822-826]. However, climate models currently cannot capture all these responses. Many eruption characteristics are poorly constrained, which may contribute to uncertainties in model solutions-for example, the season of eruption occurrence is often unknown and assigned arbitrarily. Here we isolate the effect of eruption season using experiments with the Community Earth System Model (CESM), varying the starting month of two large tropical eruptions. The eruption-year atmospheric circulation response is strongly seasonally dependent, with effects on European winter warming, the Intertropical Convergence Zone, and the southeast Asian monsoon. This creates substantial variations in eruption-year hydroclimate patterns, which do sometimes exhibit La Niña-like features as in the proxy record. However, eruption-year equatorial Pacific cooling is not driven by La Niña dynamics, but strictly by transient radiative cooling. In contrast, equatorial warming the following year occurs for all starting months and operates dynamically like El Niño. Proxy reconstructions confirm these results: eruption-year cooling is insignificant, whereas warming in the following year is more robust. This implies that accounting for the event season may be necessary to describe the initial response to volcanic eruptions and that climate models may be more accurately simulating volcanic influences than previously thought.

Northern Hemisphere forcing of Southern Hemisphere climate during the last deglaciation.

According to the Milankovitch theory, changes in summer insolation in the high-latitude Northern Hemisphere caused glacial cycles through their impact on ice-sheet mass balance. Statistical analyses of long climate records supported this theory, but they also posed a substantial challenge by showing that changes in Southern Hemisphere climate were in phase with or led those in the north. Although an orbitally forced Northern Hemisphere signal may have been transmitted to the Southern Hemisphere, insolation forcing can also directly influence local Southern Hemisphere climate, potentially intensified by sea-ice feedback, suggesting that the hemispheres may have responded independently to different aspects of orbital forcing. Signal processing of climate records cannot distinguish between these conditions, however, because the proposed insolation forcings share essentially identical variability. Here we use transient simulations with a coupled atmosphere-ocean general circulation model to identify the impacts of forcing from changes in orbits, atmospheric CO(2) concentration, ice sheets and the Atlantic meridional overturning circulation (AMOC) on hemispheric temperatures during the first half of the last deglaciation (22-14.3 kyr BP). Although based on a single model, our transient simulation with only orbital changes supports the Milankovitch theory in showing that the last deglaciation was initiated by rising insolation during spring and summer in the mid-latitude to high-latitude Northern Hemisphere and by terrestrial snow-albedo feedback. The simulation with all forcings best reproduces the timing and magnitude of surface temperature evolution in the Southern Hemisphere in deglacial proxy records. AMOC changes associated with an orbitally induced retreat of Northern Hemisphere ice sheets is the most plausible explanation for the early Southern Hemisphere deglacial warming and its lead over Northern Hemisphere temperature; the ensuing rise in atmospheric CO(2) concentration provided the critical feedback on global deglaciation.

The Holocene temperature conundrum.

A recent temperature reconstruction of global annual temperature shows Early Holocene warmth followed by a cooling trend through the Middle to Late Holocene [Marcott SA, et al., 2013, Science 339(6124):1198-1201]. This global cooling is puzzling because it is opposite from the expected and simulated global warming trend due to the retreating ice sheets and rising atmospheric greenhouse gases. Our critical reexamination of this contradiction between the reconstructed cooling and the simulated warming points to potentially significant biases in both the seasonality of the proxy reconstruction and the climate sensitivity of current climate models.

Impact of abrupt deglacial climate change on tropical Atlantic subsurface temperatures.

Both instrumental data analyses and coupled ocean-atmosphere models indicate that Atlantic meridional overturning circulation (AMOC) variability is tightly linked to abrupt tropical North Atlantic (TNA) climate change through both atmospheric and oceanic processes. Although a slowdown of AMOC results in an atmospheric-induced surface cooling in the entire TNA, the subsurface experiences an even larger warming because of rapid reorganizations of ocean circulation patterns at intermediate water depths. Here, we reconstruct high-resolution temperature records using oxygen isotope values and Mg/Ca ratios in both surface- and subthermocline-dwelling planktonic foraminifera from a sediment core located in the TNA over the last 22 ky. Our results show significant changes in the vertical thermal gradient of the upper water column, with the warmest subsurface temperatures of the last deglacial transition corresponding to the onset of the Younger Dryas. Furthermore, we present new analyses of a climate model simulation forced with freshwater discharge into the North Atlantic under Last Glacial Maximum forcings and boundary conditions that reveal a maximum subsurface warming in the vicinity of the core site and a vertical thermal gradient change at the onset of AMOC weakening, consistent with the reconstructed record. Together, our proxy reconstructions and modeling results provide convincing evidence for a subsurface oceanic teleconnection linking high-latitude North Atlantic climate to the tropical Atlantic during periods of reduced AMOC across the last deglacial transition.

Younger Dryas cooling and the Greenland climate response to CO2.

Greenland ice-core δ(18)O-temperature reconstructions suggest a dramatic cooling during the Younger Dryas (YD; 12.9-11.7 ka), with temperatures being as cold as the earlier Oldest Dryas (OD; 18.0-14.6 ka) despite an approximately 50 ppm rise in atmospheric CO(2). Such YD cooling implies a muted Greenland climate response to atmospheric CO(2), contrary to physical predictions of an enhanced high-latitude response to future increases in CO(2). Here we show that North Atlantic sea surface temperature reconstructions as well as transient climate model simulations suggest that the YD over Greenland should be substantially warmer than the OD by approximately 5 °C in response to increased atmospheric CO(2). Additional experiments with an isotope-enabled model suggest that the apparent YD temperature reconstruction derived from the ice-core δ(18)O record is likely an artifact of an altered temperature-δ(18)O relationship due to changing deglacial atmospheric circulation. Our results thus suggest that Greenland climate was warmer during the YD relative to the OD in response to rising atmospheric CO(2), consistent with sea surface temperature reconstructions and physical predictions, and has a sensitivity approximately twice that found in climate models for current climate due to an enhanced albedo feedback during the last deglaciation.

Global climate evolution during the last deglaciation.

Deciphering the evolution of global climate from the end of the Last Glacial Maximum approximately 19 ka to the early Holocene 11 ka presents an outstanding opportunity for understanding the transient response of Earth's climate system to external and internal forcings. During this interval of global warming, the decay of ice sheets caused global mean sea level to rise by approximately 80 m; terrestrial and marine ecosystems experienced large disturbances and range shifts; perturbations to the carbon cycle resulted in a net release of the greenhouse gases CO(2) and CH(4) to the atmosphere; and changes in atmosphere and ocean circulation affected the global distribution and fluxes of water and heat. Here we summarize a major effort by the paleoclimate research community to characterize these changes through the development of well-dated, high-resolution records of the deep and intermediate ocean as well as surface climate. Our synthesis indicates that the superposition of two modes explains much of the variability in regional and global climate during the last deglaciation, with a strong association between the first mode and variations in greenhouse gases, and between the second mode and variations in the Atlantic meridional overturning circulation.

Ice-shelf collapse from subsurface warming as a trigger for Heinrich events.

Episodic iceberg-discharge events from the Hudson Strait Ice Stream (HSIS) of the Laurentide Ice Sheet, referred to as Heinrich events, are commonly attributed to internal ice-sheet instabilities, but their systematic occurrence at the culmination of a large reduction in the Atlantic meridional overturning circulation (AMOC) indicates a climate control. We report Mg/Ca data on benthic foraminifera from an intermediate-depth site in the northwest Atlantic and results from a climate-model simulation that reveal basin-wide subsurface warming at the same time as large reductions in the AMOC, with temperature increasing by approximately 2 °C over a 1-2 kyr interval prior to a Heinrich event. In simulations with an ocean model coupled to a thermodynamically active ice shelf, the increase in subsurface temperature increases basal melt rate under an ice shelf fronting the HSIS by a factor of approximately 6. By analogy with recent observations in Antarctica, the resulting ice-shelf loss and attendant HSIS acceleration would produce a Heinrich event.

North Atlantic cooling triggered a zonal mode over the Indian Ocean during Heinrich Stadial 1.

Abrupt changes in the Atlantic meridional overturning circulation (AMOC) are thought to affect tropical hydroclimate through adjustment of the latitudinal position of the intertropical convergence zone (ITCZ). Heinrich Stadial 1 (HS1) involves the largest AMOC reduction in recent geological time; however, over the tropical Indian Ocean (IO), proxy records suggest zonal anomalies featuring intense, widespread drought in tropical East Africa versus generally wet but heterogeneous conditions in the Maritime Continent. Here, we synthesize proxy data and an isotope-enabled transient deglacial simulation and show that the southward ITCZ shift over the eastern IO during HS1 strengthens IO Walker circulation, triggering an east-west precipitation dipole across the basin. This dipole reverses the zonal precipitation anomalies caused by the exposed Sunda and Sahul shelves due to glacial lower sea level. Our study illustrates how zonal modes of atmosphere-ocean circulation can amplify or reverse global climate anomalies, highlighting their importance for future climate change.

Past terrestrial hydroclimate sensitivity controlled by Earth system feedbacks.

Despite tectonic conditions and atmospheric CO2 levels (pCO2) similar to those of present-day, geological reconstructions from the mid-Pliocene (3.3-3.0 Ma) document high lake levels in the Sahel and mesic conditions in subtropical Eurasia, suggesting drastic reorganizations of subtropical terrestrial hydroclimate during this interval. Here, using a compilation of proxy data and multi-model paleoclimate simulations, we show that the mid-Pliocene hydroclimate state is not driven by direct CO2 radiative forcing but by a loss of northern high-latitude ice sheets and continental greening. These ice sheet and vegetation changes are long-term Earth system feedbacks to elevated pCO2. Further, the moist conditions in the Sahel and subtropical Eurasia during the mid-Pliocene are a product of enhanced tropospheric humidity and a stationary wave response to the surface warming pattern, which varies strongly with land cover changes. These findings highlight the potential for amplified terrestrial hydroclimate responses over long timescales to a sustained CO2 forcing.

African Hydroclimate During the Early Eocene From the DeepMIP Simulations.

The early Eocene (∼56-48 Myr ago) is characterized by high CO2 estimates (1,200-2,500 ppmv) and elevated global temperatures (∼10°C-16°C higher than modern). However, the response of the hydrological cycle during the early Eocene is poorly constrained, especially in regions with sparse data coverage (e.g., Africa). Here, we present a study of African hydroclimate during the early Eocene, as simulated by an ensemble of state-of-the-art climate models in the Deep-time Model Intercomparison Project (DeepMIP). A comparison between the DeepMIP pre-industrial simulations and modern observations suggests that model biases are model- and geographically dependent, however, these biases are reduced in the model ensemble mean. A comparison between the Eocene simulations and the pre-industrial suggests that there is no obvious wetting or drying trend as the CO2 increases. The results suggest that changes to the land sea mask (relative to modern) in the models may be responsible for the simulated increases in precipitation to the north of Eocene Africa. There is an increase in precipitation over equatorial and West Africa and associated drying over northern Africa as CO2 rises. There are also important dynamical changes, with evidence that anticyclonic low-level circulation is replaced by increased south-westerly flow at high CO2 levels. Lastly, a model-data comparison using newly compiled quantitative climate estimates from paleobotanical proxy data suggests a marginally better fit with the reconstructions at lower levels of CO2.

Abrupt Heinrich Stadial 1 cooling missing in Greenland oxygen isotopes.

Abrupt climate changes during the last deglaciation have been well preserved in proxy records across the globe. However, one long-standing puzzle is the apparent absence of the onset of the Heinrich Stadial 1 (HS1) cold event around 18 ka in Greenland ice core oxygen isotope δ18 O records, inconsistent with other proxies. Here, combining proxy records with an isotope-enabled transient deglacial simulation, we propose that a substantial HS1 cooling onset did indeed occur over the Arctic in winter. However, this cooling signal in the depleted oxygen isotopic composition is completely compensated by the enrichment because of the loss of winter precipitation in response to sea ice expansion associated with AMOC slowdown during extreme glacial climate. In contrast, the Arctic summer warmed during HS1 and YD because of increased insolation and greenhouse gases, consistent with snowline reconstructions. Our work suggests that Greenland δ18 O may substantially underestimate temperature variability during cold glacial conditions.

Deglacial variability of South China hydroclimate heavily contributed by autumn rainfall.

The deglacial hydroclimate in South China remains a long-standing topic of debate due to the lack of reliable moisture proxies and inconsistent model simulations. A recent hydroclimate proxy suggests that South China became wet in cold stadials during the last deglaciation, with the intensification proposed to be contributed mostly by the East Asian summer monsoon (EASM). Here, based on a deglacial simulation in a state-of-the-art climate model that well reproduces the evolution of EASM, winter monsoon (EAWM) and the associated water isotopes in East Asia, we propose that the intensified hydroclimate in South China is also contributed heavily by the rainfall in autumn, during the transition between EASM and EAWM. The excessive rainfall in autumn results from the convergence between anomalous northerly wind due to amplified land-sea thermal contrast and anomalous southerly wind associated with the anticyclone over Western North Pacific, both of which are, in turn, forced by the slowdown of the Atlantic thermohaline circulation. Regardless the rainfall change, however, the modeled δ18Op remains largely unchanged in autumn. Our results provide new insights to East Asia monsoon associated with climate change in the North Atlantic.

Antarctic surface temperature and elevation during the Last Glacial Maximum.

Water-stable isotopes in polar ice cores are a widely used temperature proxy in paleoclimate reconstruction, yet calibration remains challenging in East Antarctica. Here, we reconstruct the magnitude and spatial pattern of Last Glacial Maximum surface cooling in Antarctica using borehole thermometry and firn properties in seven ice cores. West Antarctic sites cooled ~10°C relative to the preindustrial period. East Antarctic sites show a range from ~4° to ~7°C cooling, which is consistent with the results of global climate models when the effects of topographic changes indicated with ice core air-content data are included, but less than those indicated with the use of water-stable isotopes calibrated against modern spatial gradients. An altered Antarctic temperature inversion during the glacial reconciles our estimates with water-isotope observations.

StableClim, continuous projections of climate stability from 21000 BP to 2100 CE at multiple spatial scales.

Paleoclimatic data are used in eco-evolutionary models to improve knowledge of biogeographical processes that drive patterns of biodiversity through time, opening windows into past climate-biodiversity dynamics. Applying these models to harmonised simulations of past and future climatic change can strengthen forecasts of biodiversity change. StableClim provides continuous estimates of climate stability from 21,000 years ago to 2100 C.E. for ocean and terrestrial realms at spatial scales that include biogeographic regions and climate zones. Climate stability is quantified using annual trends and variabilities in air temperature and precipitation, and associated signal-to-noise ratios. Thresholds of natural variability in trends in regional- and global-mean temperature allow periods in Earth's history when climatic conditions were warming and cooling rapidly (or slowly) to be identified and climate stability to be estimated locally (grid-cell) during these periods of accelerated change. Model simulations are validated against independent paleoclimate and observational data. Projections of climatic stability, accessed through StableClim, will improve understanding of the roles of climate in shaping past, present-day and future patterns of biodiversity.

An Efficient Ice Sheet/Earth System Model Spin-up Procedure for CESM2-CISM2: Description, Evaluation, and Broader Applicability.

Spinning up a highly complex, coupled Earth system model (ESM) is a time consuming and computationally demanding exercise. For models with interactive ice sheet components, this becomes a major challenge, as ice sheets are sensitive to bidirectional feedback processes and equilibrate over glacial timescales of up to many millennia. This work describes and demonstrates a computationally tractable, iterative procedure for spinning up a contemporary, highly complex ESM that includes an interactive ice sheet component. The procedure alternates between a computationally expensive coupled configuration and a computationally cheaper configuration where the atmospheric component is replaced by a data model. By periodically regenerating atmospheric forcing consistent with the coupled system, the data atmosphere remains adequately constrained to ensure that the broader model state evolves realistically. The applicability of the method is demonstrated by spinning up the preindustrial climate in the Community Earth System Model Version 2 (CESM2), coupled to the Community Ice Sheet Model Version 2 (CISM2) over Greenland. The equilibrium climate state is similar to the control climate from a coupled simulation with a prescribed Greenland ice sheet, indicating that the iterative procedure is consistent with a traditional spin-up approach without interactive ice sheets. These results suggest that the iterative method presented here provides a faster and computationally cheaper method for spinning up a highly complex ESM, with or without interactive ice sheet components. The method described here has been used to develop the climate/ice sheet initial conditions for transient, ice sheet-enabled simulations with CESM2-CISM2 in the Coupled Model Intercomparison Project Phase 6 (CMIP6).

Using paleo-archives to safeguard biodiversity under climate change.

Strategies for 21st-century environmental management and conservation under global change require a strong understanding of the biological mechanisms that mediate responses to climate- and human-driven change to successfully mitigate range contractions, extinctions, and the degradation of ecosystem services. Biodiversity responses to past rapid warming events can be followed in situ and over extended periods, using cross-disciplinary approaches that provide cost-effective and scalable information for species' conservation and the maintenance of resilient ecosystems in many bioregions. Beyond the intrinsic knowledge gain such integrative research will increasingly provide the context, tools, and relevant case studies to assist in mitigating climate-driven biodiversity losses in the 21st century and beyond.

Drier tropical and subtropical Southern Hemisphere in the mid-Pliocene Warm Period.

Thermodynamic arguments imply that global mean rainfall increases in a warmer atmosphere; however, dynamical effects may result in more significant diversity of regional precipitation change. Here we investigate rainfall changes in the mid-Pliocene Warm Period (~ 3 Ma), a time when temperatures were 2-3ºC warmer than the pre-industrial era, using output from the Pliocene Model Intercomparison Projects phases 1 and 2 and sensitivity climate model experiments. In the Mid-Pliocene simulations, the higher rates of warming in the northern hemisphere create an interhemispheric temperature gradient that enhances the southward cross-equatorial energy flux by up to 48%. This intensified energy flux reorganizes the atmospheric circulation leading to a northward shift of the Inter-Tropical Convergence Zone and a weakened and poleward displaced Southern Hemisphere Subtropical Convergences Zones. These changes result in drier-than-normal Southern Hemisphere tropics and subtropics. The evaluation of the mid-Pliocene adds a constraint to possible future warmer scenarios associated with differing rates of warming between hemispheres.

Past climates inform our future.

As the world warms, there is a profound need to improve projections of climate change. Although the latest Earth system models offer an unprecedented number of features, fundamental uncertainties continue to cloud our view of the future. Past climates provide the only opportunity to observe how the Earth system responds to high carbon dioxide, underlining a fundamental role for paleoclimatology in constraining future climate change. Here, we review the relevancy of paleoclimate information for climate prediction and discuss the prospects for emerging methodologies to further insights gained from past climates. Advances in proxy methods and interpretations pave the way for the use of past climates for model evaluation-a practice that we argue should be widely adopted.

Glacial changes in tropical climate amplified by the Indian Ocean.

The mechanisms driving glacial-interglacial changes in the climate of the Indo-Pacific warm pool are poorly understood. Here, we address this question by combining paleoclimate proxies with model simulations of the Last Glacial Maximum climate. We find evidence of two mechanisms explaining key patterns of ocean cooling and rainfall change interpreted from proxy data. Exposure of the Sahul shelf excites a positive ocean-atmosphere feedback involving a stronger surface temperature gradient along the equatorial Indian Ocean and a weaker Walker circulation-a response explaining the drier/wetter dipole across the basin. Northern Hemisphere cooling by ice sheet albedo drives a monsoonal retreat across Africa and the Arabian Peninsula-a response that triggers a weakening of the Indian monsoon via cooling of the Arabian Sea and associated reductions in moisture supply. These results demonstrate the importance of air-sea interactions in the Indian Ocean, amplifying externally forced climate changes over a large part of the tropics.