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: 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: 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: 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.

Description: The stable water isotopes H218O and HDO are incorporated as passive tracers into the oceanic general circulation model MPI-OM, and a control simulation under present-day climate conditions is analyzed in detail. Both δ18O and δD distributions at the ocean surface and deep ocean are generally consistent with available observations on the large scale. The modelled δD-δ 18O relations in surface waters slightly deviates from the slope of the global meteoric water line in most basins, and a much steeper slope is detected in Arctic Oceans. The simulated deuterium excess of ocean surface waters shows small variations between 80° S and 55° N, and a strong decrease north of 55° N. The model is also able to capture the quasi-linear relationship between δ18O and salinity S, as well as δD and S, as seen in observational data. Both in the model results and observations, the surface δ−S relations show a steeper slope in extra-tropical regions than in tropical regions, which indicates relatively more addition of isotopically depleted water at high latitudes.

Description: Interpreting stable oxygen isotope (δ18O) records from stalagmites is still one of the complex tasks in speleothem research. Here, we present a novel model-based approach, where we force a model describing the processes and modifications of δ18O from rain water to speleothem calcite (Oxygen isotope Drip water and Stalagmite Model – ODSM) with the results of a state-of-the-art atmospheric general circulation model enhanced by explicit isotope diagnostics (ECHAM5-wiso). The approach is neither climate nor cave-specific and allows an integrated assessment of the influence of different varying climate variables, e.g. temperature and precipitation amount, on the isotopic composition of drip water and speleothem calcite. First, we apply and evaluate this new approach under present-day climate conditions using observational data from seven caves from different geographical regions in Europe. Each of these caves provides measured δ18O values of drip water and speleothem calcite to which we compare our simulated isotope values. For six of the seven caves modeled δ18O values of drip water and speleothem calcite are in good agreement with observed values. The mismatch of the remaining caves might be caused by the complexity of the cave system, beyond the parameterizations included in our cave model. We then examine the response of the cave system to mid-Holocene (6000yr before present, 6ka) climate conditions by forcing the ODSM with ECHAM5-wiso results from 6ka simulations. For a set of twelve European caves, we compare the modeled mid-Holocene-to-modern difference in speleothem calcite δ18O to available measurements. We show that the general European changes are simulated well. However, local discrepancies are found, and might be explained either by a too low model resolution, complex local soil-atmosphere interactions affecting evapotranspiration or by cave specific factors such as non-equilibrium fractionation processes. The mid-Holocene experiment pronounces the potential of the presented approach to analyse δ18O variations on a spatially large (regional to global) scale. Modelled as well as measured European δ18O values of stalagmite samples suggest the presence of a strong, positive mode of the North Atlantic Oscillation at 6 ka before present, which is supported by the respective modelled climate parameters.

Description: Modelling the evolution of the Earth system on long timescales requires the knowledge and understanding of driving mechanisms, such as the hydrological cycle. This is dominant in all components of the Earth’s system, such as atmosphere, ocean, land surfaces/vegetation and the cryosphere. Observations and measurements of stable water isotopes in climate archives can help to decipher and reconstruct climate change and its regional variations. For the cryosphere, the δ18O cycle in the current generation of Earth-System-Models is missing and an efficient and accurate tracer transport scheme is required. We describe ISOPOLIS 1.0 a modular semi-Lagrangian transport scheme of second-order accuracy, which is coupled to the polythermal and thermomechanical ice sheet model SICOPOLIS (version 2.9). Model skill is demonstrated by experiments with a simplified ice sheet geometry and by comparisons of simulated ice cores with data from Greenland (GRIP) and Antarctica (Vostok). The presented method is a valuable tool to investigate the transport of any kind of passive tracer inside a polythermal ice sheet and is an important step to model the whole δ18O cycle.

Description: Sea ice in both polar regions is an important indicator of the expression of global climate change and its polar amplification. Consequently, broad interest exists on sea ice coverage, variability and long-term change. However, its predictability is complex and it depends strongly on different atmospheric and oceanic parameters. In order to provide insights into the potential development of a monthly/seasonal signal of sea ice evolution, we applied a robust statistical model based on different oceanic and atmospheric parameters to calculate an estimate of the September sea ice extent (SSIE) on a monthly timescale. Although previous statistical attempts of monthly/seasonal SSIE forecasts show a relatively reduced skill, when the trend is removed, we show here that the September sea ice extent has a high predictive skill, up to 4 months ahead, based on previous months’ oceanic and atmospheric conditions. Our statistical model skillfully captures the interannual variability of the SSIE and could provide a valuable tool for identifying relevant regions and oceanic and atmospheric parameters that are important for the sea ice development in the Arctic and for detecting sensitive/critical regions in global coupled climate models with a focus on sea ice formation.

Description: Holocene sea surface temperature trends and variability are underestimated in models compared to paleoclimate data. The idea is presented that the local trends and variability are related, which is elaborated in a conceptual framework of the stochastic climate model. The relation is a consequence of the fluctuation–dissipation theorem, connecting the linear response of a system to its statistical fluctuations. Consequently, the spectrum can be used to estimate the timescale-dependent climate response. The non-normality in the propagation operator introduces enhanced long-term variability related to nonequilibrium and/or Earth system sensitivity. The simple model can guide us to analyze comprehensive models' behavior.