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

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

Evaluation of coupled ocean-atmosphere simulations of northern hemisphere extratropical climates in the mid-Holocene

We have used the BIOME4 biogeography–biochemistry model and comparison with palaeovegetation data to evaluate the response of six ocean–atmosphere general circulation models to mid-Holocene changes in orbital forcing in the mid- to high-latitudes of the northern hemisphere. All the models produce: (a) a northward shift of the northern limit of boreal forest, in response to simulated summer warming in high-latitudes. The northward shift is markedly asymmetric, with larger shifts in Eurasia than in North America; (b) an expansion of xerophytic vegetation in mid-continental North America and Eurasia, in response to increased temperatures during the growing season; (c) a northward expansion of temperate forests in eastern North America, in response to simulated winter warming. The northward shift of the northern limit of boreal forest and the northward expansion of temperate forests in North America are supported by palaeovegetation data. The expansion of xerophytic vegetation in mid-continental North America is consistent with palaeodata, although the extent may be over-estimated. The simulated expansion of xerophytic vegetation in Eurasia is not supported by the data. Analysis of an asynchronous coupling of one model to an equilibrium-vegetation model suggests vegetation feedback exacerbates this mid-continental drying and produces conditions more unlike the observations. Not all features of the simulations are robust: some models produce winter warming over Europe while others produce winter cooling. As a result, some models show a northward shift of temperate forests (consistent with, though less marked than, the expansion shown by data) and others produce a reduction in temperate forests. Elucidation of the cause of such differences is a focus of the current phase of the Palaeoclimate Modelling Intercomparison Project

Simulating the mid-Pliocene Warm Period with the CCSM4 model

This paper describes the experimental design and model results from a 500 yr fully coupled Community Climate System, version 4, simulation of the mid-Pliocene Warm Period (mPWP) (ca. 3.3–3.0 Ma). We simulate the mPWP using the "alternate" protocol prescribed by the Pliocene Model Intercomparison Project (PlioMIP) for the AOGCM simulation (Experiment 2). Results from the CCSM4 mPWP simulation show a 1.9 °C increase in global mean annual temperature compared to the 1850 preindustrial control, with a polar amplification of ~3 times the global warming. Global precipitation increases slightly by 0.09 mm day−1 and the monsoon rainfall is enhanced, particularly in the Northern Hemisphere (NH). Areal sea ice extent decreases in both hemispheres but persists through the summers. The model simulates a relaxation of the zonal sea surface temperature (SST) gradient in the tropical Pacific, with the El Niño–Southern Oscillation (Niño3.4) ~20% weaker than the preindustrial and exhibiting extended periods of quiescence of up to 150 yr. The maximum Atlantic meridional overturning circulation and northward Atlantic oceanic heat transport are indistinguishable from the control. As compared to PRISM3, CCSM4 overestimates Southern Hemisphere (SH) sea surface temperatures, but underestimates NH warming, particularly in the North Atlantic, suggesting that an increase in northward ocean heat transport would bring CCSM4 SSTs into better alignment with proxy data

The Connected Isotopic Water Cycle in the Community Earth System Model Version 1

Because of the pervasive role of water in the Earth system, the relative abundances of stable isotopologues of water are valuable for understanding atmospheric, oceanic, and biospheric processes, and for interpreting paleoclimate proxy reconstructions. Isotopologues are transported by both largeâ scale and turbulent flows, and the ratio of heavy to light isotopologues changes due to fractionation that can accompany condensation and evaporation processes. Correctly predicting the isotopic distributions requires resolving the relationships between largeâ scale ocean and atmospheric circulation and smallerâ scale hydrological processes, which can be accomplished within a coupled climate modeling framework. Here we present the water isotopeâ enabled version of the Community Earth System Model version 1 (iCESM1), which simulates global variations in water isotopic ratios in the atmosphere, land, ocean, and sea ice. In a transient Last Millennium simulation covering the 850â 2005 period, iCESM1 correctly captures the lateâ twentiethâ century structure of δ18O and δD over the global oceans, with more limited accuracy over land. The relationship between salinity and seawater δ18O is also well represented over the observational period, including interbasin variations. We illustrate the utility of coupled, isotopeâ enabled simulations using both Last Millennium simulations and freshwater hosing experiments with iCESM1. Closing the isotopic mass balance between all components of the coupled model provides new confidence in the underlying depiction of the water cycle in CESM, while also highlighting areas where the underlying hydrologic balance can be improved. The iCESM1 is poised to be a vital community resource for ongoing model development with both modern and paleoclimate applications.Key PointsAn isotopeâ enabled version of the Community Earth System Model (iCESM1) is now publicly availableiCESM1 simulates the major observed features of δ18O and δD over the late twentieth centuryiCESM1 is useful for both modern climate and paleoclimate applicationsPeer Reviewedhttps://deepblue.lib.umich.edu/bitstream/2027.42/151857/1/jame20931.pdfhttps://deepblue.lib.umich.edu/bitstream/2027.42/151857/2/jame20931_am.pd

Pliocene and Eocene provide best analogs for near-future climates

As the world warms due to rising greenhouse gas concentrations, the Earth system moves toward climate states without societal precedent, challenging adaptation. Past Earth system states offer possible model systems for the warming world of the coming decades. These include the climate states of the Early Eocene (ca. 50 Ma), the Mid-Pliocene (3.3–3.0 Ma), the Last Interglacial (129–116 ka), the Mid-Holocene (6 ka), preindustrial (ca. 1850 CE), and the 20th century. Here, we quantitatively assess the similarity of future projected climate states to these six geohistorical benchmarks using simulations from the Hadley Centre Coupled Model Version 3 (HadCM3), the Goddard Institute for Space Studies Model E2-R (GISS), and the Community Climate System Model, Versions 3 and 4 (CCSM) Earth system models. Under the Representative Concentration Pathway 8.5 (RCP8.5) emission scenario, by 2030 CE, future climates most closely resemble Mid-Pliocene climates, and by 2150 CE, they most closely resemble Eocene climates. Under RCP4.5, climate stabilizes at Pliocene-like conditions by 2040 CE. Pliocene-like and Eocene-like climates emerge first in continental interiors and then expand outward. Geologically novel climates are uncommon in RCP4.5 (<1%) but reach 8.7% of the globe under RCP8.5, characterized by high temperatures and precipitation. Hence, RCP4.5 is roughly equivalent to stabilizing at Pliocene-like climates, while unmitigated emission trajectories, such as RCP8.5, are similar to reversing millions of years of long-term cooling on the scale of a few human generations. Both the emergence of geologically novel climates and the rapid reversion to Eocene-like climates may be outside the range of evolutionary adaptive capacity

Hydroclimate footprint of pan-Asian monsoon water isotope during the last deglaciation

Oxygen isotope speleothem records exhibit coherent variability over the pan-Asian summer monsoon (AM) region. The hydroclimatic representation of these oxygen isotope records for the AM, however, has remained poorly understood. Here, combining an isotope-enabled Earth system model in transient experiments with proxy records, we show that the widespread AM delta O-18(c) signal during the last deglaciation (20 to 11 thousand years ago) is accompanied by a continental-scale, coherent hydroclimate footprint, with spatially opposite signs in rainfall. This footprint is generated as a dynamically coherent response of the AM system primarily to meltwater forcing and secondarily to insolation forcing and is further reinforced by atmospheric teleconnection. Hence, widespread delta O-18(p) depletion in the AM region is accompanied by a northward migration of the westerly jet and enhanced southwesterly monsoon wind, as well as increased rainfall from South Asia (India) to northern China but decreased rainfall in southeast China

The Pliocene Model Intercomparison Project (PlioMIP) Phase 2: Scientific Objectives and Experimental Design

The Pliocene Model Intercomparison Project (PlioMIP) is a co-ordinated international climate modelling initiative to study and understand climate and environments of the Late Pliocene, and their potential relevance in the context of future climate change. PlioMIP examines the consistency of model predictions in simulating Pliocene climate, and their ability to reproduce climate signals preserved by geological climate archives. Here we provide a description of the aim and objectives of the next phase of the model intercomparison project (PlioMIP Phase 2), and we present the experimental design and boundary conditions that will be utilised for climate model experiments in Phase 2. Following on from PlioMIP Phase 1, Phase 2 will continue to be a mechanism for sampling structural uncertainty within climate models. However, Phase 1 demonstrated the requirement to better understand boundary condition uncertainties as well as uncertainty in the methodologies used for data-model comparison. Therefore, our strategy for Phase 2 is to utilise state-of-the-art boundary conditions that have emerged over the last 5 years. These include a new palaeogeographic reconstruction, detailing ocean bathymetry and land/ice surface topography. The ice surface topography is built upon the lessons learned from offline ice sheet modelling studies. Land surface cover has been enhanced by recent additions of Pliocene soils and lakes. Atmospheric reconstructions of palaeo-CO2 are emerging on orbital timescales and these are also incorporated into PlioMIP Phase 2. New records of surface and sea surface temperature change are being produced that will be more temporally consistent with the boundary conditions and forcings used within models. Finally we have designed a suite of prioritized experiments that tackle issues surrounding the basic understanding of the Pliocene and its relevance in the context of future climate change in a discrete way

No Consistent Simulated Trends in the Atlantic Meridional Overturning Circulation for the Past 6,000 Years

The Atlantic Meridional Overturning Circulation (AMOC) is a key feature of the North Atlantic with global ocean impacts. The AMOC's response to past changes in forcings during the Holocene provides important context for the coming centuries. Here, we investigate AMOC trends using an emerging set of transient simulations using multiple global climate models for the past 6,000 years. Although some models show changes, no consistent trend in overall AMOC strength during the mid-to-late Holocene emerges from the ensemble. We interpret this result to suggest no overall change in AMOC, which fits with our assessment of available proxy reconstructions. The decadal variability of the AMOC does not change in ensemble during the mid- and late-Holocene. There are interesting AMOC changes seen in the early Holocene, but their nature depends a lot on which inputs are used to drive the experiment

Arctic sea ice simulation in the PlioMIP ensemble

Eight general circulation models have simulated the mid-Pliocene warm period (mid-Pliocene, 3.264 to 3.025 Ma) as part of the Pliocene Modelling Intercomparison Project (PlioMIP). Here, we analyse and compare their simulation of Arctic sea ice for both the pre-industrial period and the mid-Pliocene. Mid-Pliocene sea ice thickness and extent is reduced, and the model spread of extent is more than twice the pre-industrial spread in some summer months. Half of the PlioMIP models simulate ice-free conditions in the mid-Pliocene. This spread amongst the ensemble is in line with the uncertainties amongst proxy reconstructions for mid-Pliocene sea ice extent. Correlations between mid-Pliocene Arctic temperatures and sea ice extents are almost twice as strong as the equivalent correlations for the pre-industrial simulations. The need for more comprehensive sea ice proxy data is highlighted, in order to better compare model performances