Description: We present here the first results, for the preindustrial and mid-Holocene climatological periods, of the newly developed isotope-enhanced version of the fully coupled Earth system model MPI-ESM, called hereafter MPI-ESM-wiso. The water stable isotopes H216O, H218O and HDO have been implemented into all components of the coupled model setup. The mid-Holocene provides the opportunity to evaluate the model response to changes in the seasonal and latitudinal distribution of insolation induced by different orbital forcing conditions. The results of our equilibrium simulations allow us to evaluate the performance of the isotopic model in simulating the spatial and temporal variations of water isotopes in the different compartments of the hydrological system for warm climates. For the preindustrial climate, MPI-ESM-wiso reproduces very well the observed spatial distribution of the isotopic content in precipitation linked to the spatial variations in temperature and precipitation rate. We also find a good model-data agreement with the observed distribution of isotopic composition in surface seawater but a bias with the presence of surface seawater that is too 18O-depleted in the Arctic Ocean. All these results are improved compared to the previous model version ECHAM5/MPIOM. The spatial relationships of water isotopic composition with temperature, precipitation rate and salinity are consistent with observational data. For the preindustrial climate, the interannual relationships of water isotopes with temperature and salinity are globally lower than the spatial ones, consistent with previous studies. Simulated results under mid-Holocene conditions are in fair agreement with the isotopic measurements from ice cores and continental speleothems. MPI-ESM-wiso simulates a decrease in the isotopic composition of precipitation from North Africa to the Tibetan Plateau via India due to the enhanced monsoons during the mid-Holocene. Over Greenland, our simulation indicates a higher isotopic composition of precipitation linked to higher summer temperature and a reduction in sea ice, shown by positive isotope-temperature gradient. For the Antarctic continent, the model simulates lower isotopic values over the East Antarctic plateau, linked to the lower temperatures during the mid-Holocene period, while similar or higher isotopic values are modeled over the rest of the continent. While variations of isotopic contents in precipitation over West Antarctica between mid-Holocene and preindustrial periods are partly controlled by changes in temperature, the transport of relatively 18O-rich water vapor near the coast to the western ice core sites could play a role in the final isotopic composition. So, more caution has to be taken about the reconstruction of past temperature variations during warm periods over this area. The coupling of such a model with an ice sheet model or the use of a zoomed grid centered on this region could help to better describe the role of the water vapor transport and sea ice around West Antarctica. The reconstruction of past salinity through isotopic content in sea surface waters can be complicated for regions with strong ocean dynamics, variations in sea ice regimes or significant changes in freshwater budget, giving an extremely variable relationship between the isotopic content and salinity of ocean surface waters over small spatial scales. These complicating factors demonstrate the complexity of interpreting water isotopes as past climate signals of warm periods like the mid-Holocene. A systematic isotope model intercomparison study for further insights on the model dependency of these results would be beneficial.

Description: This data set contains model values from 4 simulations produced with the isotope-enabled atmosphere GCMs ECHAM5-wiso and ECHAM6-wiso for the period 1979-2018. These simulations been performed at different spatial resolutions (T63 and T127) or with different reanalyses for the nudging (ERA5 and ERA-Interim). A complete description can be found in Cauquoin and Werner (J. Adv. Model. Earth Syst., 2021).

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 investigate the impact of spatial resolution on the discrepancy between simulated δ18O and observed δ18O in Greenland ice cores, regional climate simulations are performed with the isotope-enabled regional climate model (RCM) COSMO_iso. For this purpose, isotope-enabled general circulation model (GCM) simulations with the ECHAM5-wiso general circulation model (GCM) under present-day conditions and the MPI-ESM-wiso GCM under mid-Holocene conditions are dynamically downscaled with COSMO_iso for the Arctic region. The capability of COSMO_iso to reproduce observed isotopic ratios in Greenland ice cores for these two periods is investigated by comparing the simulation results to measured δ18O ratios from snow pit samples, Global Network of Isotopes in Precipitation (GNIP) stations and ice cores. To our knowledge, this is the first time that a mid-Holocene isotope-enabled RCM simulation is performed for the Arctic region. Under present-day conditions, a dynamical downscaling of ECHAM5-wiso (1.1∘×1.1∘ ) with COSMO_iso to a spatial resolution of 50 km improves the agreement with the measured δ18O ratios for 14 of 19 observational data sets. A further increase in the spatial resolution to 7 km does not yield substantial improvements except for the coastal areas with its complex terrain. For the mid-Holocene, a fully coupled MPI-ESM-wiso time slice simulation is downscaled with COSMO_iso to a spatial resolution of 50 km. In the mid-Holocene, MPI-ESM-wiso already agrees well with observations in Greenland and a downscaling with COSMO_iso does not further improve the model–data agreement. Despite this lack of improvement in model biases, the study shows that in both periods, observed δ18O values at measurement sites constitute isotope ratios which are mainly within the subgrid-scale variability of the global ECHAM5-wiso and MPI-ESM-wiso simulation results. The correct δ18O ratios are consequently not resolved in the GCM simulation results and need to be extracted by a refinement with an RCM. In this context, the RCM simulations provide a spatial δ18O distribution by which the effects of local uncertainties can be taken into account in the comparison between point measurements and model outputs. Thus, an isotope-enabled GCM–RCM model chain with realistically implemented fractionating processes constitutes a useful supplement to reconstruct regional paleo-climate conditions during the mid-Holocene in Greenland. Such model chains might also be applied to reveal the full potential of GCMs in other regions and climate periods, in which large deviations relative to observed isotope ratios are simulated.

Description: Stable water isotopes in polar ice cores are widely used to reconstruct past temperature variations over several orbital climatic cycles. One way to calibrate the isotope–temperature relationship is to apply the present-day spatial relationship as a surrogate for the temporal one. However, this method leads to large uncertainties because several factors like the sea surface conditions or the origin and transport of water vapor influence the isotope–temperature temporal slope. In this study, we investigate how the sea surface temperature (SST), the sea ice extent, and the strength of the Atlantic Meridional Overturning Circulation (AMOC) affect these temporal slopes in Greenland and Antarctica for Last Glacial Maximum (LGM, ∼ 21 000 years ago) to preindustrial climate change. For that, we use the isotope-enabled atmosphere climate model ECHAM6-wiso, forced with a set of sea surface boundary condition datasets based on reconstructions (e.g., GLOMAP) or MIROC 4m simulation outputs. We found that the isotope–temperature temporal slopes in East Antarctic coastal areas are mainly controlled by the sea ice extent, while the sea surface temperature cooling affects the temporal slope values inland more. On the other hand, ECHAM6-wiso simulates the impact of sea ice extent on the EPICA Dome C (EDC) and Vostok sites through the contribution of water vapor from lower latitudes. Effects of sea surface boundary condition changes on modeled isotope–temperature temporal slopes are variable in West Antarctica. This is partly due to the transport of water vapor from the Southern Ocean to this area that can dampen the influence of local temperature on the changes in the isotopic composition of precipitation and snow. In the Greenland area, the isotope–temperature temporal slopes are influenced by the sea surface temperatures near the coasts of the continent. The greater the LGM cooling off the coast of southeastern Greenland, the greater the transport of water vapor from the North Atlantic, and the larger the temporal slopes. The presence or absence of sea ice very near the coast has a large influence in Baffin Bay and the Greenland Sea and influences the slopes at some inland ice core stations. The extent of the sea ice far south slightly influences the temporal slopes in Greenland through the transport of more depleted water vapor from lower latitudes to this area. The seasonal variations of sea ice distribution, especially its retreat in summer, influence the isotopic composition of the water vapor in this region and the modeled isotope–temperature temporal slopes in the eastern part of Greenland. A stronger LGM AMOC decreases LGM-to-preindustrial isotopic anomalies in precipitation in Greenland, degrading the isotopic model–data agreement. The AMOC strength modifies the temporal slopes over inner Greenland slightly and by a little on the coasts along the Greenland Sea where the changes in surface temperature and sea ice distribution due to the AMOC strength mainly occur.

Description: Stable hydrogen isotopes in monsoonal precipitation (δDp) at three sites (Port Blair, Barisal and Darjeeling) reveal the factors governing δDp variations over a south-north gradient across the Bay of Bengal. We found that the δDp at each site continuously decreases from May to September and these trends become more pronounced from south to north. The decreasing trends of downstream δDp closely follow the decreasing trends of upstream stable hydrogen isotopes in water vapor (δDv), which indicates that upstream δDv properties shape initial spatiotemporal patterns of the downstream δDp (“shaping effect”). Additionally, our results demonstrate that, during moisture transport, upstream vertical air motions (convection and downward motion) and topographic relief magnify the amplitude of the decreasing trends of downstream δD (“magnifying effect”). Our findings imply that upstream δD properties and relevant atmospheric and pv topographical conditions along the moisture transport pathway need to be considered collectively to better interpret paleoclimate records.

Description: The understanding of the relationship between the variation of precipitation stable oxygen isotope ratio (δ18Op) and monsoon activity in the Asian monsoon region is crucial for an in-depth comprehension of the regional hydrological cycle processes and for reconstructing the history of Asian paleomonsoon changes. Based on the 1979–2017 summer δ18Op output by two isotope-enabled atmospheric general circulation models nudged to climate reanalysis data, this study explores the associations of the Indian summer monsoon (IM) and western North Pacific summer monsoon (WNPM) intensities with the interannual variations of the regional δ18Op and their possible physical mechanisms. Statistical analyses demonstrate that the East Asian δ18Op is negatively correlated with the IM intensity while the Indian δ18Op is positively correlated with the WNPM intensity. Moreover, the underlying mechanisms linking the monsoon and δ18Op vary in different regions. In strong IM years, with the intensified convection and increased precipitation near the Indian peninsula, the water vapor isotope ratio (δ18Ov) transported to East Asia has lower values, resulting in the depletion of δ18Op there. The opposite is true for weak IM years. In years of strong WNPM, the intensified convection over the tropical western Pacific and the suppressed convection over the western Indian Ocean may be linked to a Walker-type circulation anomaly, accompanied by the enlarging of the vertical wind shear between the western Pacific and the western Indian Ocean. Accordingly, the decreasing of convection and precipitation over the Arabian Sea results in higher δ18Ov values in the upstream area of India, which ultimately increases δ18Op values in the Indian peninsula through the monsoonal moisture transport; and vice versa.

Description: Numerical simulations employing prognostic sta- ble water isotopes can not only facilitate our understanding of hydrological processes and climate change but also al- low 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 We- gener 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) simula- tions. The model reproduces the observed PI isotope compo- sitions in both precipitation and seawater well and captures their major differences from the MH conditions. The sim- ulated relationship between the isotope composition in precipitation (d18Op) and surface air temperature is very similar between the PI and MH conditions, and it is largely consis- tent 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 pre- cipitation 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 dis- tinction between pre-monsoon and monsoon seasons. Com- bining 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 with- drawal 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: table water isotopes in polar ice cores are widely used to reconstruct past temperature variations over several orbital climatic cycles. One way to calibrate the isotope– temperature relationship is to apply the present-day spatial relationship as a surrogate for the temporal one. However, this method leads to large uncertainties because several fac- tors like the sea surface conditions or the origin and transport of water vapor influence the isotope–temperature temporal slope. In this study, we investigate how the sea surface tem- perature (SST), the sea ice extent, and the strength of the Atlantic Meridional Overturning Circulation (AMOC) affect these temporal slopes in Greenland and Antarctica for Last Glacial Maximum (LGM, ∼ 21 000 years ago) to preindus- trial climate change. For that, we use the isotope-enabled at- mosphere climate model ECHAM6-wiso, forced with a set of sea surface boundary condition datasets based on recon- structions (e.g., GLOMAP) or MIROC 4m simulation out- puts. We found that the isotope–temperature temporal slopes in East Antarctic coastal areas are mainly controlled by the sea ice extent, while the sea surface temperature cooling af- fects the temporal slope values inland more. On the other hand, ECHAM6-wiso simulates the impact of sea ice extent on the EPICA Dome C (EDC) and Vostok sites through the contribution of water vapor from lower latitudes. Effects of sea surface boundary condition changes on modeled isotope-emperature temporal slopes are variable in West Antarctica. This is partly due to the transport of water vapor from the Southern Ocean to this area that can dampen the influence of local temperature on the changes in the isotopic compo- sition of precipitation and snow. In the Greenland area, the isotope–temperature temporal slopes are influenced by the sea surface temperatures near the coasts of the continent. The greater the LGM cooling off the coast of southeastern Green- land, the greater the transport of water vapor from the North Atlantic, and the larger the temporal slopes. The presence or absence of sea ice very near the coast has a large influ- ence in Baffin Bay and the Greenland Sea and influences the slopes at some inland ice core stations. The extent of the sea ice far south slightly influences the temporal slopes in Greenland through the transport of more depleted water vapor from lower latitudes to this area. The seasonal vari- ations of sea ice distribution, especially its retreat in sum- mer, influence the isotopic composition of the water vapor in this region and the modeled isotope–temperature tempo- ral slopes in the eastern part of Greenland. A stronger LGM AMOC decreases LGM-to-preindustrial isotopic anomalies in precipitation in Greenland, degrading the isotopic model– data agreement. The AMOC strength modifies the temporal slopes over inner Greenland slightly and by a little on the coasts along the Greenland Sea where the changes in surface temperature and sea ice distribution due to the AMOC strength mainly occur.

Description: In order to complement the picture of the atmospheric water cycle in the Southern Ocean, we have continuously monitored water vapor isotopes since January 2020 on Amsterdam Island in the Indian Ocean. We present here the first 2-year long water vapor isotopic record at this site. We show that the water vapor isotopic composition largely follows the water vapor mixing ratio, as expected in marine boundary layers. However, we detect 11 periods of a few days where there is a strong loss of correlation between water vapor δ18O and water vapor mixing ratio as well as abrupt negative excursions of water vapor δ18O. These excursions often occur toward the end of precipitation events. Six of these events show a decrease in gaseous elemental mercury, suggesting subsidence of air from a higher altitude. Our study aims to further explore the mechanism driving these negative excursions in water vapor δ18O. We used two different models to provide a data–model comparison over this 2-year period. While the European Centre Hamburg model (ECHAM6-wiso) at 0.9° was able to reproduce most of the sharp negative water vapor δ18O excursions, hence validating the physics process and isotopic implementation in this model, the Laboratoire de Météorologie Dynamique Zoom model (LMDZ-iso) at 2° (3°) resolution was only able to reproduce seven (one) of the negative excursions, highlighting the possible influence of the model resolution for the study of such abrupt isotopic events. Based on our detailed model–data comparison, we conclude that the most plausible explanations for such isotopic excursions are rain–vapor interactions associated with subsidence at the rear of a precipitation event.