Model sensitivity to North Atlantic freshwater forcing at 8.2 ka

We compared four simulations of the 8.2 ka event to assess climate model sensitivity and skill in responding to North Atlantic freshwater perturbations. All of the simulations used the same freshwater forcing, 2.5 Sv for one year, applied to either the Hudson Bay (northeastern Canada) or Labrador Sea (between Canada's Labrador coast and Greenland). This freshwater pulse induced a decadal-mean slowdown of 10–25% in the Atlantic Meridional Overturning Circulation (AMOC) of the models and caused a large-scale pattern of climate anomalies that matched proxy evidence for cooling in the Northern Hemisphere and a southward shift of the Intertropical Convergence Zone. The multi-model ensemble generated temperature anomalies that were just half as large as those from quantitative proxy reconstructions, however. Also, the duration of AMOC and climate anomalies in three of the simulations was only several decades, significantly shorter than the duration of ~150 yr in the paleoclimate record. Possible reasons for these discrepancies include incorrect representation of the early Holocene climate and ocean state in the North Atlantic and uncertainties in the freshwater forcing estimates

Climate response to large, high-latitude and low-latitude volcanic eruptions in the Community Climate System Model

Explosive volcanism is known to be a leading natural cause of climate change. The second half of the 13th century was likely the most volcanically perturbed half-century of the last 2000 years, although none of the major 13th century eruptions have been clearly attributed to specific volcanoes. This period was in general a time of transition from the relatively warm Medieval period to the colder Little Ice Age, but available proxy records are insufficient on their own to clearly assess whether this transition is associated with volcanism. This context motivates our investigation of the climate system sensitivity to high- and low-latitude volcanism using the fully coupled NCAR Community Climate System Model (CCSM3). We evaluate two sets of ensemble simulations, each containing four volcanic pulses, with the first set representing them as a sequence of tropical eruptions and the second representing eruptions occurring in the mid-high latitudes of both the Northern and Southern hemispheres. The short-term, direct radiative impacts of tropical and high- latitude eruptions include significant cooling over the continents in summer and cooling over regions of increased sea-ice concentration in Northern Hemisphere (NH) winter. A main dynamical impact of moderate tropical eruptions is a winter warming pattern across northern Eurasia. Furthermore, both ensembles show significant reductions in global precipitation, especially in the summer monsoon regions. The most important long-term impact is the cooling of the high-latitude NH produced by multiple tropical eruptions, suggesting that positive feedbacks associated with ice and snow cover could lead to long-term climate cooling in the Arctic

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 sub-thermocline-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

Paleoclimate constraints on the spatiotemporal character of past and future droughts

Author Posting. © American Meteorological Society, 2020. This article is posted here by permission of American Meteorological Society for personal use, not for redistribution. The definitive version was published in Journal of Climate 33(22), (2020): 9883-9903, https://doi.org/10.1175/JCLI-D-20-0004.1.Machine-learning-based methods that identify drought in three-dimensional space–time are applied to climate model simulations and tree-ring-based reconstructions of hydroclimate over the Northern Hemisphere extratropics for the past 1000 years, as well as twenty-first-century projections. Analyzing reconstructed and simulated drought in this context provides a paleoclimate constraint on the spatiotemporal characteristics of simulated droughts. Climate models project that there will be large increases in the persistence and severity of droughts over the coming century, but with little change in their spatial extent. Nevertheless, climate models exhibit biases in the spatiotemporal characteristics of persistent and severe droughts over parts of the Northern Hemisphere. We use the paleoclimate record and results from a linear inverse modeling-based framework to conclude that climate models underestimate the range of potential future hydroclimate states. Complicating this picture, however, are divergent changes in the characteristics of persistent and severe droughts when quantified using different hydroclimate metrics. Collectively our results imply that these divergent responses and the aforementioned biases must be better understood if we are to increase confidence in future hydroclimate projections. Importantly, the novel framework presented herein can be applied to other climate features to robustly describe their spatiotemporal characteristics and provide constraints on future changes to those characteristics.This material is based upon work supported by the National Center for Atmospheric Research, which is a major facility sponsored by the National Science Foundation under Cooperative Agreement 1852977. JAF was also supported by the Regional and Global Model Analysis (RGMA) component of the Earth and Environmental System Modeling Program of the U.S. Department of Energy's Office of Biological & Environmental Research (BER) via National Science Foundation IA 1844590. JS was supported in part by the U.S. National Science Foundation through Grants AGS-1602920 and AGS-1805490, and by the National Oceanic and Atmospheric Administration by Grant NA20OAR4310425. We acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups (listed in Table 1) for producing and making available their model output. For CMIP, the U.S. Department of Energy’s Program for Climate Model Diagnosis and Intercomparison provides coordinating support and led development of software infrastructure in partnership with the Global Organization for Earth System Science Portal. We thank the editor and two reviewers for comments that greatly improved the quality of this manuscript. This is SOEST Publication No. 11116 and LDEO Publication No. 8450.2021-04-1

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.Stuart C. Brown, Tom M. L. Wigley, Bette L. Otto-Bliesner and Damien A. Fordha

Modeling Northern Hemisphere ice-sheet distribution during MIS 5 and MIS 7 glacial inceptions

The present manuscript compares Marine Iso- tope Stage 5 (MIS 5, 125–115 kyr BP) and MIS 7 (236– 229 kyr BP) with the aim to investigate the origin of the difference in ice-sheet growth over the Northern Hemi- sphere high latitudes between these last two inceptions. Our approach combines a low resolution coupled atmosphere– ocean–sea-ice general circulation model and a 3-D thermo- mechanical ice-sheet model to simulate the state of the ice sheets associated with the inception climate states of MIS 5 and MIS 7. Our results show that external forcing (orbitals and GHG) and sea-ice albedo feedbacks are the main fac- tors responsible for the difference in the land-ice initial state between MIS 5 and MIS 7 and that our cold climate model bias impacts more during a cold inception, such as MIS 7, than during a warm inception, such as MIS 5. In addition, if proper ice-elevation and albedo feedbacks are not taken into consideration, the evolution towards glacial inception is hardly simulated, especially for MIS 7. Finally, results high- light that while simulated ice volumes for MIS 5 glacial in- ception almost fit with paleo-reconstructions, the lack of pre- cipitation over high latitudes, identified as a bias of our cli- mate model, does not allow for a proper simulation of MIS 7 glacial inception

Pliocene Warmth Consistent With Greenhouse Gas Forcing

With CO2 concentrations similar to today (410 ppm), the Pliocene Epoch offers insights into climate changes under a moderately warmer world. Previous work suggested a low zonal sea surface temperature (SST) gradient in the tropical Pacific during the Pliocene, the so‐called “permanent El Niño.” Here, we recalculate SSTs using the alkenone proxy and find moderate reductions in both the zonal and meridional SST gradients during the mid‐Piacenzian warm period. These reductions are captured by coupled climate model simulations of the Pliocene, especially those that simulate weaker Walker circulation. We also produce a spatial reconstruction of mid‐Piacenzian warm period Pacific SSTs that closely resembles both Pliocene and future, low‐emissions simulations, a pattern that is, to a first order, diagnostic of weaker Walker circulation. Therefore, Pliocene warmth does not require drastic changes in the climate system—rather, it supports the expectation that the Walker circulation will weaken in the future under higher CO2

LGM permafrost distribution: how well can the latest PMIP multi-model ensembles perform reconstruction?

Here, global-scale frozen ground distribution from the Last Glacial Maximum (LGM) has been reconstructed using multi-model ensembles of global climate models, and then compared with evidence-based knowledge and earlier numerical results. Modeled soil temperatures, taken from Paleoclimate Modelling Intercomparison Project phase III (PMIP3) simulations, were used to diagnose the subsurface thermal regime and determine underlying frozen ground types for the present day (pre-industrial; 0 kya) and the LGM (21 kya). This direct method was then compared to an earlier indirect method, which categorizes underlying frozen ground type from surface air temperature, applying to both the PMIP2 (phase II) and PMIP3 products. Both direct and indirect diagnoses for 0 kya showed strong agreement with the present-day observation-based map. The soil temperature ensemble showed a higher diversity around the border between permafrost and seasonally frozen ground among the models, partly due to varying subsurface processes, implementation, and settings. The area of continuous permafrost estimated by the PMIP3 multi-model analysis through the direct (indirect) method was 26.0 (17.7) million km2 for LGM, in contrast to 15.1 (11.2) million km2 for the pre-industrial control, whereas seasonally frozen ground decreased from 34.5 (26.6) million km2 to 18.1 (16.0) million km2. These changes in area resulted mainly from a cooler climate at LGM, but from other factors as well, such as the presence of huge land ice sheets and the consequent expansion of total land area due to sea-level change. LGM permafrost boundaries modeled by the PMIP3 ensemble-improved over those of the PMIP2 due to higher spatial resolutions and improved climatology-also compared better to previous knowledge derived from geomorphological and geocryological evidence. Combinatorial applications of coupled climate models and detailed stand-alone physical-ecological models for the cold-region terrestrial, paleo-, and modern climates will advance our understanding of the functionality and variability of the frozen ground subsystem in the global eco-climate system

Critical evaluation of climate syntheses to benchmark CMIP6/PMIP4 127 ka Last Interglacial simulations in the high-latitude regions

The Last Interglacial (LIG, ∼129-116 thousand years ago, ka) represents an excellent case study to investigate the response of sensitive components of the Earth System and mechanisms of high-lati tude amplification to a climate warmer than present-day. The Paleoclimate Model Intercomparison Project (Phase 4, hereafter referred as PMIP4) and the Coupled Model Intercomparison Project (Phase 6, hereafter referred as CMIP6) are coordinating the design of (1) a LIG Tier 1 equilibrium simulation to simulate the climate response at 127 ka, a time interval associated with a strong orbital forcing and greenhouse gas concentrations close to preindustrial levels and (2) associated Tier 2 sensitivity experiments to examine the role of the ocean, vegetation and dust feedbacks in modulating the response to this orbital forcing. Evaluating the capability of the CMIP6/PMIP4 models to reproduce the 127 ka polar and sub-polar climate will require appropriate data-based benchmarks which are currently missing. Based on a recent data synthesis that offers the first spatio-temporal representation of high-latitude (i.e. poleward of 40°N and 40°S) surface temperature evolution during the LIG, we produce a new 126–128 ka time slab, hereafter named 127 ka time slice. This 127 ka time slice represents surface temperature anomalies relative to preindustrial and is associated with quantitative estimates of the uncertainties related to relative dating and surface temperature reconstruction methods. It illustrates warmer-than-preindustrial conditions in the high-latitude regions of both hemispheres. In particular, summer sea surface temperatures (SST) in the North Atlantic region were on average 1.1 °C (with a standard error of the mean of 0.7 °C) warmer relative to preindustrial and 1.8 °C (with a standard error of the mean of 0.8 °C) in the Southern Ocean. In Antarctica, average 127 ka annual surface air temperature was 2.2 °C (with a standard error of the mean of 1.4 °C) warmer compared to preindustrial. We provide a critical evaluation of the latest LIG surface climate compilations that are available for evaluating LIG climate model experiments. We discuss in particular our new 127 ka time-slice in the context of existing LIG surface temperature time-slices. We also compare the 127 ka time slice with the ones published for the 125 and 130 ka time intervals and we discuss the potential and limits of a data-based time slice at 127 ka in the context of the upcoming coordinated modeling exercise. Finally we provide guidance on the use of the available LIG climate compilations for future model-data comparison exercises in the framework of the upcoming CMIP6/PMIP4 127 ka experiments. We do not recommend the use of LIG peak warmth-centered syntheses. Instead we promote the use of the most recent syntheses that are based on coherent chronologies between paleoclimatic records and provide spatio-temporal reconstruction of the LIG climate. In particular, we recommend using our new 127 ka data-based time slice in model-data comparison studies with a focus on the high-latitude climate.E. C. is funded by the European Union's Seventh Framework Programme for research and innovation under the Marie Skłodowska-Curie grant agreement no 600207. B. L. O-B is supported by the U.S. National Science Foundation (NSF) sponsorship of NCAR. R. F. acknowledges the funding of the NSF Arctic System Science. E.W.W. is supported by the Royal Society. This is LSCE contribution no 6117