Description: Numerical simulations employing prognostic stable water isotopes can not only facilitate our understanding of hydrological processes and climate change but also allow for a direct comparison between isotope signals obtained from models and various archives. In the current work, we describe the performance and explore the potential of a new version of the Earth system model AWI-ESM (Alfred Wegener Institute Earth System Model), labeled AWI-ESM-2.1-wiso, in which we incorporated three isotope tracers into all relevant components of the water cycle. We present here the results of pre-industrial (PI) and mid-Holocene (MH) simulations. The model reproduces the observed PI isotope compositions in both precipitation and seawater well and captures their major differences from the MH conditions. The simulated relationship between the isotope composition in precipitation (δ18Op) and surface air temperature is very similar between the PI and MH conditions, and it is largely consistent with modern observations despite some regional model biases. The ratio of the MH–PI difference in δ18Op to the MH–PI difference in surface air temperature is comparable to proxy records over Greenland and Antarctica only when summertime air temperature is considered. An amount effect is evident over the North African monsoon domain, where a negative correlation between δ18Op and the amount of precipitation is simulated. As an example of model applications, we studied the onset and withdrawal date of the MH West African summer monsoon (WASM) using daily variables. We find that defining the WASM onset based on precipitation alone may yield erroneous results due to the substantial daily variations in precipitation, which may obscure the distinction between pre-monsoon and monsoon seasons. Combining precipitation and isotope indicators, we suggest in this work a novel method for identifying the commencement of the WASM. Moreover, we do not find an obvious difference between the MH and PI periods in terms of the mean onset of the WASM. However, an advancement in the WASM withdrawal is found in the MH compared to the PI period due to an earlier decline in insolation over the northern location of Intertropical Convergence Zone (ITCZ).

Description: In order to quantify the relative importance of individual boundary conditions and forcings, including greenhouse gases, ice sheets, and Earth's orbital parameters, on determining Last Glacial Maximum (LGM) climate, we have performed a series of LGM experiments using a state-of-the-art climate model AWI-ESM, in which different combinations of boundary conditions and forcings have been applied following the protocol of Paleoclimate Modelling Intercomparison Project phase 4 (PMIP4). In good agreement with observational proxy records, a general colder and drier climate is simulated in our full-forced LGM experiment as compared to the present-day simulation. Our simulated results from non-full-forced sensitivity simulations reveal that both the greenhouse gases and ice sheets play a major role in defining the anomalous LGM surface temperature compared to today. Decreased greenhouse gases in LGM as compared to present day leads to a non-uniform global cooling with polar amplification effect. The presence of LGM ice sheets favors a warming over the Arctic and northern Atlantic oceans in boreal winter, as well as a cooling over regions with the presence of ice sheets. The former is induced by a strengthening in the Atlantic meridional overturning circulation (AMOC) transporting more heat to high latitudes, whilst the latter is due to the increased surface albedo and elevation of ice sheets. We find that the Northern Hemisphere monsoon precipitation is influenced by the opposing effects of LGM greenhouse gases and ice sheets. Specifically, the presence of ice sheets leads to significant drying in the Northern Hemisphere monsoon regions, while a reduction in greenhouse gases results in increased monsoon rainfall. Based on our model results, continental ice sheets exert a major control on atmospheric dynamics and the variability of El Niño–Southern Oscillation (ENSO). Moreover, our analysis also implies a nonlinearity in climate response to LGM boundary conditions and forcings.

Description: During the Middle Miocene climate transition about 14 million years ago, the Antarctic ice sheet expanded to near-modern volume. Surprisingly, this ice sheet growth was accompanied by a warming in the surface waters of the Southern Ocean, whereas a slight deep-water temperature increase was delayed by more than 200 thousand years. Here we use a coupled atmosphere–ocean model to assess the relative effects of changes in atmospheric CO2 concentration and ice sheet growth on regional and global temperatures. In the simulations, changes in the wind field associated with the growth of the ice sheet induce changes in ocean circulation, deep-water formation and sea-ice cover that result in sea surface warming and deep-water cooling in large swaths of the Atlantic and Indian ocean sectors of the Southern Ocean. We interpret these changes as the dominant ocean surface response to a 100-thousand-year phase of massive ice growth in Antarctica. A rise in global annual mean temperatures is also seen in response to increased Antarctic ice surface elevation. In contrast, the longer-term surface and deep-water temperature trends are dominated by changes in atmospheric CO2 concentration. We therefore conclude that the climatic and oceanographic impacts of the Miocene expansion of the Antarctic ice sheet are governed by a complex interplay between wind field, ocean circulation and the sea-ice system.

Description: During glacial periods of the Late Pleistocene, an abundance of proxy data demonstrates the existence of large and repeated millennial-scale warming episodes, known as Dansgaard–Oeschger (DO) events1. This ubiquitous feature of rapid glacial climate change can be extended back as far as 800,000 years before present (BP) in the ice core record2, and has drawn broad attention within the science and policy-making communities alike3. Many studies have been dedicated to investigating the underlying causes of these changes, but no coherent mechanism has yet been identified3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15. Here we show, by using a comprehensive fully coupled model16, that gradual changes in the height of the Northern Hemisphere ice sheets (NHISs) can alter the coupled atmosphere–ocean system and cause rapid glacial climate shifts closely resembling DO events. The simulated global climate responses—including abrupt warming in the North Atlantic, a northward shift of the tropical rainbelts, and Southern Hemisphere cooling related to the bipolar seesaw—are generally consistent with empirical evidence1, 3, 17. As a result of the coexistence of two glacial ocean circulation states at intermediate heights of the ice sheets, minor changes in the height of the NHISs and the amount of atmospheric CO2 can trigger the rapid climate transitions via a local positive atmosphere–ocean–sea-ice feedback in the North Atlantic. Our results, although based on a single model, thus provide a coherent concept for understanding the recorded millennial-scale variability and abrupt climate changes in the coupled atmosphere–ocean system, as well as their linkages to the volume of the intermediate ice sheets during glacials.

Description: Understanding how the Antarctic ice sheet will respond to global warming relies on knowledge of how it has behaved in the past. The use of numerical models, the only means to quantitatively predict the future, is hindered by limitations to topographic data both now and in the past, and in knowledge of how subsurface oceanic, glaciological and hydrological processes interact. Incorporating the variety and interplay of such processes, operating at multiple spatio-temporal scales, is critical to modeling the Antarctic’s system evolution and requires direct observations in challenging locations. As these processes do not observe disciplinary boundaries neither should our future research.

Description: The equilibrium climate sensitivity (ECS) of climate models is calculated as the equilibrium global mean surface air warming resulting from a simulated doubling of the atmospheric CO2 concentration. In these simulations, long-term processes in the climate system, such as land ice changes, are not incorporated. Hence, climate sensitivity derived from paleodata has to be compensated for these processes, when comparing it to the ECS of climate models. Several recent studies found that the impact these long-term processes have on global temperature cannot be quantified directly through the global radiative forcing they induce. This renders the prevailing approach of deconvoluting paleotemperatures through a partitioning based on radiative forcings inaccurate. Here, we therefore implement an efficacy factor ε[LI] that relates the impact of land ice changes on global temperature to that of CO2 changes in our calculation of climate sensitivity from paleodata. We apply our refined approach to a proxy-inferred paleoclimate dataset, using ε[LI]=0.45+0.34−0.20 based on a multi-model assemblage of simulated relative influences of land ice changes on the Last Glacial Maximum temperature anomaly. The implemented ε[LI] is smaller than unity, meaning that per unit of radiative, forcing the impact on global temperature is less strong for land ice changes than for CO2 changes. Consequently, our obtained ECS estimate of 5.8±1.3 K, where the uncertainty reflects the implemented range in ε[LI], is ∼50 % higher than when differences in efficacy are not considered.

Description: The evolution of past global ice sheets is highly uncertain. One example is the missing ice problem during the Last Glacial Maximum (LGM, 26 000-19 000 years before present) – an apparent 8-28 m discrepancy between far-field sea level indicators and modelled sea level from ice sheet reconstructions. In the absence of ice sheet reconstructions, researchers often use marine δ 18 O proxy records to infer ice volume prior to the LGM. We present a global ice sheet reconstruction for the past 80 000 years, called PaleoMIST 1.0, constructed inde- pendently of far-field sea level and δ 18 O proxy records. Our reconstruction is compatible with LGM far-field sea-level records without requiring extra ice volume, thus solving the missing ice problem. However, for Marine Isotope Stage 3 (57 000-29 000 years before present) - a pre-LGM period - our reconstruction does not match proxy-based sea level reconstructions, indicating the relationship between marine δ 18 O and sea level may be more complex than assumed.

Description: Understanding the dynamics of warm climate states has gained increasing importance in the face of anthropogenic climate change. During the Last Interglacial (LIG, ∼128 to 116 ka), greenhouse gas concentrations and high latitude insolation were higher than pre-industrial levels, causing a high-latitude warming (Turney and Jones, 2010; Pfeiffer and Lohmann, 2016). We present a suite of climate model results (COSMOS, MPI-ESM, AWI-CM, EC-Earth) to evaluate the patterns and compare the simulations with the above-mentioned surface temperature reconstructions, seasonal archives (Felis et al., 2015; Brocas et al., 2017), and sea ice reconstructions (Stein et al., 2017). As a result of this modestly warmer climate, polar ice sheets were smaller and estimates report that the global mean sea level was 6-9 meters higher than today (Dutton et al., 2015). The sensitivity of the Antarctic Ice sheet is related to the local temperature around the West Antarctic Ice Sheet (WAIS) (Sutter et al., 2016). Our ice sheet model experiments indicate that a 2-3°C local warming causes already a partially collapsed, irreversible WAIS. A pronounced subsurface oceanic warming can destabilize the WAIS, resulting in an oceanic gateway between the Ross and Weddell Seas. A sensitivity study using the new oceanic gateway between the Atlantic and Pacific Oceans as a bathymetrical boundary condition indicates that this region would be covered by sea ice. Mixing due to sea-ice formation prevents a pronounced warming around the WAIS and would stabilize the WAIS. Thus, the disintegration of the WAIS is probably related to non-local influences like in Hellmer et al. (2017) where the shelves of West Antarctica are warmed from below by Circumpolar Deep Water.

Description: Stable water isotope records from Antarctica are key for our understanding of Quaternary climate variations. However, the exact quantitative interpretation of these important climate proxy records in terms of surface temperature, ice sheet height and other climatic changes is still a matter of debate. Here we report results obtained with an atmospheric general circulation model equipped with water isotopes, run at a high-spatial horizontal resolution of one-by-one degree. Comparing different glacial maximum ice sheet reconstructions, a best model data match is achieved for the PMIP3 reconstruction. Reduced West Antarctic elevation changes between 400 and 800 m lead to further improved agreement with ice core data. Our modern and glacial climate simulations support the validity of the isotopic paleothermometer approach based on the use of present-day observations and reveal that a glacial ocean state as displayed in the GLAMAP reconstruction is suitable for capturing the observed glacial isotope changes in Antarctic ice cores.

Description: In this manuscript we describe the experimental procedure employed at the Alfred Wegener Institute in Germany in the preparation of the simulations for the Pliocene Model Intercomparison Project (PlioMIP). We present a description of the utilized Community Earth System Models (COSMOS, version: COSMOS-landveg r2413, 2009) and document the procedures that we applied to transfer the Pliocene Research, Interpretation and Synoptic Mapping (PRISM) Project mid-Pliocene reconstruction into model forcing fields. The model setup and spin-up procedure are described for both the paleo- and preindustrial (PI) time slices of PlioMIP experiments 1 and 2, and general results that depict the performance of our model setup for mid-Pliocene conditions are presented. The mid-Pliocene, as simulated with our COSMOS setup and PRISM boundary conditions, is both warmer and wetter in the global mean than the PI. The globally averaged annual mean surface air temperature in the mid-Pliocene standalone atmosphere (fully coupled atmosphere-ocean) simulation is 17.35 °C (17.82 °C), which implies a warming of 2.23 °C (3.40 °C) relative to the respective PI control simulation.