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: Meridional atmospheric transport is an important process in the climate system and has implications for the availability of heat and moisture at high latitudes. Near-surface cold and warm temperature advection over the ocean in the context of extratropical cyclones additionally leads to important air–sea exchange. In this paper, we investigate the impact of these air–sea fluxes on the stable water isotope (SWI) composition of water vapour in the Southern Ocean’s atmospheric boundary layer. SWIs serve as a tool to trace phase change processes involved in the atmospheric water cycle and, thus, provide important insight into moist atmospheric processes associated with extratropical cyclones. Here we combine a 3-month ship-based SWI measurement data set around Antarctica with a series of regional high-resolution numerical model simulations from the isotope-enabled numerical weather prediction model COSMOiso. We objectively identify atmospheric cold and warm temperature advection associated with the cold and warm sector of extratropical cyclones, respectively, based on the air–sea temperature difference applied to the measurement and the simulation data sets. A Lagrangian composite analysis of temperature advection based on the COSMOiso simulation data is compiled to identify the main processes affecting the observed variability of the isotopic signal in marine boundary layer water vapour in the region from 35 to 70◦ S. This analysis shows that the cold and warm sectors of extratropical cyclones are associated with contrasting SWI signals. Specifically, the measurements show that the median values of δ18O and δ2H in the atmospheric water vapour are 3.8 ‰ and 27.9 ‰ higher during warm than during cold advection. The median value of the second-order isotope variable deuterium excess d, which can be used as a measure of non-equilibrium processes during phase changes, is 6.4 ‰ lower during warm than during cold advection. These characteristic isotope signals during cold and warm advection reflect the opposite air–sea fluxes associated with these large-scale transport events. The trajectory-based analysis reveals that the SWI signals in the cold sector are mainly shaped by ocean evaporation. In the warm sector, the air masses experience a net loss of moisture due to dew deposition as they are advected over the relatively colder ocean, which leads to the observed low d. We show that additionally the formation of clouds and precipitation in moist adiabatically ascending warm air parcels can decrease d in boundary layer water vapour. These findings illustrate the highly variable isotopic composition in water vapour due to contrasting air–sea interactions during cold and warm advection, respectively, induced by the circulation associated with extratropical cyclones. SWIs can thus potentially be useful as tracers for meridional air advection and other characteristics associated with the dynamics of the storm tracks over interannual timescales.

Description: Winter (October to March) precipitation δ18OP and δDP values in central Europe correlate with the winter North Atlantic Oscillation index (wNAOi), but the causal mechanisms remain poorly understood. Here we analyse the relationships between precipitation-weighted δ18OP and δDP datasets (δ18Opw and δDpw) from European GNIP and ANIP stations and the wNAOi, with a focus on isotope gradients. We demonstrate that longitudinal δ18Opw and δDpw gradients across Europe (“continental effect”) depend on the wNAOi state, with steeper gradients associated with more negative wNAOi states. Changing gradients reflect a combination of air temperature and variable amounts of precipitable water as a function of the wNAOi. The relationships between the wNAOi, δ18Opw and δDpw can provide additional information from palaeoclimate archives such as European speleothems that primarily record winter δ18Opw. Comparisons between present-day and past European longitudinal δ18O gradients inferred from Holocene speleothems suggest that atmospheric pressure configurations akin to negative wNAO modes dominated the early Holocene, whereas patterns resembling positive wNAO modes were more common in the late Holocene, possibly caused by persistent shifts in the relative locations of the Azores High and the Icelandic Low.

Description: The project MUSICA (MUlti-platform remote Sensing of Isotopologues for investigating the Cycle of Atmospheric water) has shown that the sensor IASI aboard the satellite MetOp can measure the free tropospheric {H2O,δD} pair distribution twice per day on a quasi-global scale. Such data are very promising for investigating tropospheric moisture pathways, however, the complex data characteristics compromise their usage in the context of model evaluation studies. Here we present a tool that allows for simulating MUSICA MetOp/IASI {H2O,δD} pair remote sensing data for a given model atmosphere, thereby creating model data that have the remote sensing data characteristics assimilated. This model data can then be compared to the MUSICA data. The retrieval simulation method is based on the physical principles of radiative transfer and we show that the uncertainty of the simulations is within the uncertainty of the MUSICA MetOp/IASI products, i.e. the retrieval simulations are reliable enough. We demonstrate the working principle of the simulator by applying it to ECHAM5-wiso model data. The few case studies clearly reveal the large potential of the MUSICA MetOp/IASI {H2O,δD} data pairs for evaluating modelled moisture pathways. The tool is made freely available in form of MATLAB and Python routines and can be easily connected to any atmospheric water vapour isotopologue model.

Description: Variations in oxygen isotope compositions (δ18O) provide insight into modern climate and past changes in climate and topography. In addition, in regions such as Tibet, geologic archives of isotope ratios record climate change driven by plateau uplift and therefore also provide information about the surface uplift history. A good understanding of modern-day controls on δ18O is crucial for interpreting geologic δ18O in this context. We use the ECHAM5-wiso global atmospheric general circulation model to calculate δ18O in precipitation (δ18Op) for the present-day climate. In the region of the Tibetan Plateau, spatial variations of monthly means of δ18Op are statistically related to spatial variations of 2 m air temperature and precipitation rate, as well as to topography. The size and location of investigated regions are varied in our study to capture regional differences in these relationships and the processes governing the modern δ18Op. In addition to correlation analyses, a cross-validated stepwise multiple regression is carried out using δ18Op as the predictand, and topography and atmospheric variables (temperature and precipitation amount) as predictors. The 2 m air temperature and topography yield the highest spatial correlation coefficients of 〉0.9 and 〈 0.9, respectively, throughout most of the year. Particularly high correlation coefficients are calculated for the region along the Himalayan orogeny and parts of western China. The predictors explain 〉90% of the δ18Op spatial variance in the same regions. The 2 m air temperature is the dominant predictor and contributes 93.6% to the total explained spatial variance on average. The results demonstrate that most of the δ18Op pattern on and around the Tibetan Plateau can be explained by variation in 2 m air temperature and altitude. Correlation of the dependent predictors indicate that in low-altitude regions where topography does not determine temperature variability, the high correlation of temperature and δ18Op may partially be explained by variations in precipitation rates.

Description: The skills of isotope-enabled general circulation models are evaluated against atmospheric water vapor isotopes. We have combined in situ observations of surface water vapor isotopes spanning multiple field seasons (2010, 2011, and 2012) from the top of the Greenland Ice Sheet (NEEM site: 77.45°N, 51.05°W, 2484 m above sea level) with observations from the marine boundary layer of the North Atlantic and Arctic Ocean (Bermuda Islands 32.26°N, 64.88°W, year: 2012; south coast of Iceland 63.83°N, 21.47°W, year: 2012; South Greenland 61.21°N, 47.17°W, year: 2012; Svalbard 78.92°N, 11.92°E, year: 2014). This allows us to benchmark the ability to simulate the daily water vapor isotope variations from five different simulations using isotope-enabled general circulation models. Our model-data comparison documents clear isotope biases both on top of the Greenland Ice Sheet (1–11‰ for δ18O and 4–19‰ for d-excess depending on model and season) and in the marine boundary layer (maximum differences for the following: Bermuda δ18O = ~1‰, d-excess = ~3‰; South coast of Iceland δ18O = ~2‰, d-excess = ~ 5‰; South Greenland δ18O = ~4‰, d-excess = ~7‰; Svalbard δ18O = ~2‰, d-excess = ~7‰). We find that the simulated isotope biases are not just explained by simulated biases in temperature and humidity. Instead, we argue that these isotope biases are related to a poor simulation of the spatial structure of the marine boundary layer water vapor isotopic composition. Furthermore, we specifically show that the marine boundary layer water vapor isotopes of the Baffin Bay region show strong influence on the water vapor isotopes at the NEEM deep ice core-drilling site in northwest Greenland. Our evaluation of the simulations using isotope-enabled general circulation models also documents wide intermodel spatial variability in the Arctic. This stresses the importance of a coordinated water vapor isotope-monitoring network in order to discriminate amongst these model behaviors

Description: This modelling study aims at an improved understanding of the processes that determine the water vapour budget in the stratosphere by means of the investigation of water isotope ratios. An additional (and separate from the actual) hydrological cycle has been introduced into the chemistry–climate model EMAC, including the water isotopologues HDO and H2 18 O and their physical fractionation processes. Additionally an explicit computation of the contribution of methane oxidation to H2O and HDO has been incorporated. The model expansions allow detailed analyses of water vapour and its isotope ratio with respect to deuterium throughout the stratosphere and in the transition region to the troposphere. In order to assure the correct representation of the water isotopologues in the model’s hydrological cycle, the expanded system has been evaluated in several steps. The physical fractionation effects have been evaluated by comparison of the simulated isotopic composition of precipitation with measurements from a ground-based network (GNIP) and with the results from the isotopologue-enabled general circulation model ECHAM5-wiso. The model’s representation of the chemical HDO precursor CH3D in the stratosphere has been confirmed by a comparison with chemical transport models (1-D, CHEM2D) and measurements from radiosonde flights. Finally, the simulated stratospheric HDO and the isotopic composition of water vapour have been evaluated, with respect to retrievals from three different satellite instruments (MIPAS, ACE-FTS, SMR). Discrepancies in stratospheric water vapour isotope ratios between two of the three satellite retrievals can now partly be explained.

Description: Reconstructions of Quaternary climate are often based on the isotopic content of paleo-precipitation preserved in proxy records. While many paleo-precipitation isotope records are available, few studies have synthesized these dispersed records to explore spatial patterns of late-glacial precipitation δ18O. Here we present a synthesis of 86 globally distributed groundwater (n = 59), 
cave calcite (n = 15) and ice core (n = 12)
 isotope records spanning the late-glacial (defined as
 ~ 50000 to ∼ 20000 years ago) to the late-Holocene (within the past ∼5000 years). We show that precipitation δ18O changes from the late-glacial to the late-Holocene range from −7.1 ‰ (δ18Olate-Holocene 〉 δ18Olate-glacial) to +1.7 ‰ (δ18Olate-glacial 〉 δ18Olate-Holocene), with the majority (77 %) of records having lower late-glacial δ18O than late-Holocene δ18O values. High-magnitude, negative precipitation δ18O shifts are common at high latitudes, high altitudes and continental interiors (δ18Olate-Holocene 〉 δ18Olate-glacial by more than 3‰). Conversely, low-magnitude, positive precipitation δ18O shifts are concentrated along tropical and subtropical coasts (δ18Olate-glacial 〉 δ18Olate-Holocene by less than 2 ‰). Broad, global patterns of late-glacial to late-Holocene precipitation δ18O shifts suggest that stronger-than-modern isotopic distillation of air masses prevailed during the late-glacial, likely impacted by larger global temperature differences between the tropics and the poles. Further, to test how well general circulation models reproduce global precipitation δ18O shifts, we compiled simulated precipitation δ18O shifts from five isotope-enabled general circulation models simulated under recent and last glacial maximum climate states. Climate simulations generally show better inter-model and model-measurement agreement in temperate regions than in the tropics, highlighting a need for further research to better understand how inter-model spread in convective rainout, seawater δ18O and glacial topography parameterizations impact simulated precipitation δ18O. Future research on paleo-precipitation δ18O records can use the global maps of measured and simulated late-glacial precipitation isotope compositions to target and prioritize field sites.

Description: The oxygen isotopes of water (H218O and H216O) are tracers widely used for the investigation of Earth science problems. The tracer applications are based on the premise that the 18O/16O ratio of open-water evaporation (δ18O ) can be calculated from environmental conditions. A long-standing issue concerns the role of kinetic fractionation, or diffusion transport, in the evaporation process. Here we deployed an optical instrument at a large lake (area 2,400 km2) to make in situ measurement of δ18O and δD of atmospheric vapor, then determined δ18O and δD of open-water evaporation using the gradient-diffusion method. Our results show a much weaker kinetic effect than suggested by the kinetic factor εk adopted in some previous studies of lake hydrology (14.2‰). By incorporating into the H218O isotopic mass balance of the lake a lower εk value (about 6.2‰) used for ocean evaporation in global climate models, we obtain an annual lake evaporation rate that agrees with an independent eddy-covariance observation, but the rate is 72% higher than if the commonly used lake εk value of 14.2‰ is applied. The applicability of this results to small lakes is uncertain and in need of field-based assessment.