ALBERT

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: 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: 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: 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: Central Asia is one of the largest arid regions in the world, however, multiple lakes have existed here since the Neogene. These lakes were able to sustain themselves despite the aridification trend in Asia through the PlioPleistocene. For example, long-term geological multiproxy records, including carbonate δ18O, from lake sediments of the Qaidam, Gaxun Nur, and Orog Nuur Basins indicate multiple changes in the hydrological cycle of the region with alternate phases of prevailing evaporation and precipitation. These changes are attributed either to Neogene global climate change or regional tectonic events. In this study, we use the isotope-equipped atmospheric general circulation model ECHAM5-wiso for modeling of Asia climate evolution and associated changes in precipitation δ18O during key periods of the Neogene. High-resolution simulations (T159L31, ca. 0.8°x0.8° and 31 vertical levels, 6 hour output frequency) with Mid-Holocene, Pleistocene, Pliocene and Miocene boundary conditions allow us to estimate the contributions of global climate change into the hydrological budget over the Central Asia. We complement this work with a Lagrangian Trajectory analysis (wind back-trajectories) applied to the ECHAM5-wiso outputs to trace changes in the origin of precipitation-producing air masses. We show that in addition to precipitation amount variations associated with changes in large-scale atmosphere dynamics, considerable changes in moisture sources between the time slices considered contribute to the isotopic signature of precipitation within the Qaidam, Gaxun Nur, and Orog Nuur Basins. Finally, comparison of simulated δ18O results to wind backtrajectory analysis suggests that local process, such as moisture recycling, exert an increasing control for more recent time periods.

Description: A 22.4 m-long shallow firn core was extracted during the 2006/2007 field season from coastal Adélie Land. Annual layer counting based on subannual analyses of δ18O and major chemical components was combined with 5 reference years associated with nuclear tests and non-retreat of summer sea ice to build the initial ice-core chronology (1946– 2006), stressing uncertain counting for 8 years. We focus here on the resulting δ18O and accumulation records. With an average value of 21.8 ± 6.9 cm w.e. yr−1 , local accumulation shows multi-decadal variations peaking in the 1980s, but no long-term trend. Similar results are obtained for δ18O, also characterised by a remarkably low and variable amplitude of the seasonal cycle. The ice-core records are compared with regional records of temperature, stake area accumulation measurements and variations in sea-ice extent, and outputs from two models nudged to ERA (European Reanalysis) atmospheric reanalyses: the high-resolution atmospheric general circulation model (AGCM), including stable water isotopes ECHAM5-wiso (European Centre Hamburg model), and the regional atmospheric model Modèle Atmosphérique Régional ( AR). A significant linear correlation is identified between decadal variations in δ18O and regional temperature. No significant relationship appears with regional sea-ice extent. A weak and significant correlation appears with Dumont d’Urville wind speed, increasing after 1979. The model-data comparison highlights the inadequacy of ECHAM5-wiso simulations prior to 1979, possibly due to the lack of data assimilation to constrain atmospheric reanalyses. Systematic biases are identified in the ECHAM5-wiso simulation, such as an overestimation of the mean accumulation rate and its interannual variability, a strong cold bias and an underestimation of the mean δ18O value and its interannual variability. As a result, relationships between simulated δ18O and temperature are weaker than observed. Such systematic precipitation and temperature biases are not displayed by MAR, suggesting that the model resolution plays a key role along the Antarctic ice sheet coastal topography. Interannual variations in ECHAM5-wiso temperature and precipitation accurately capture signals from meteorological data and stake observations and are used to refine the initial ice-core chronology within 2 years. After this adjustment, remarkable positive (negative) δ18O anomalies are identified in the ice-core record and the ECHAM5-wiso simulation in 1986 and 2002 (1998–1999), respectively. Despite uncertainties associated with post-deposition processes and signal-to-noise issues, in one single coastal ice-core record, we conclude that the S1C1 core can correctly capture major annual anomalies in δ18O as well as multi-decadal variations. These findings highlight the importance of improving the network of coastal high-resolution ice-core records, and stress the skills and limitations of atmospheric models for accumulation and δ18O in coastal Antarctic areas. This is particularly important for the overall East Antarctic ice sheet mass balance.

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.