Description: Understanding climate proxy records that preserve physical characteristics of past climate is a prerequisite to reconstruct long‐term climatic conditions. Water stable isotope ratios (δ18O) constitute a widely used proxy in ice cores to reconstruct temperature and climate. However, the original climate signal is altered between the formation of precipitation and the ice, especially in low‐accumulation areas such as the East Antarctic Plateau. Atmospheric conditions under which the isotopic signal is acquired at Aurora Basin North (ABN), East Antarctica, are characterized with the regional atmospheric model Modèle Atmosphérique Régional (MAR). The model shows that 50% of the snow is accumulated in less than 24 days/year. Snowfall occurs throughout the year and intensifies during winter, with 64% of total accumulation between April and September, leading to a cold bias of −0.86°C in temperatures above inversion compared to the annual mean of −29.7°C. Large snowfall events are associated with high‐pressure systems forcing warm oceanic air masses toward the Antarctic interior, which causes a warm bias of +2.83°C. The temperature‐δ18O relationship, assessed with the global atmospheric model ECHAM5‐wiso, is primarily constrained by the winter variability, but the observed slope is valid year‐round. Three snow δ18O records covering 2004–2014 indicate that the anomalies recorded in the ice core are attributable to the occurrence of warm winter storms bringing precipitation to ABN and support the interpretation of δ18O in this region as a marker of temperature changes related to large‐scale atmospheric conditions, particularly blocking events and variations in the Southern Annular Mode.

Description: The research area of climate field reconstructions has developed strongly during the past 20 years, motivated by the need to understand the complex dynamics of the earth system in a changing climate. Climate field reconstructions aim to build a consistent gridded climate reconstruction of different variables, often from a range of climate proxies, using either statistical tools or a climate model to fill the gaps between the locations of the proxy data. Commonly, large-scale climate field reconstructions covering more than 500 years are of annual resolution. In this method study, we investigate the potential of seasonally resolved climate field reconstructions based on oxygen isotope records from Greenland ice cores and an isotope-enabled climate model. Our analogue-type method matches modeled isotope patterns in Greenland precipitation to the patterns of ice core data from up to 14 ice core sites. In a second step, the climate variables of the best-matching model years are extracted, with the mean of the best-matching years comprising the reconstruction. We test a range of climate reconstructions, varying the definition of the seasons and the number of ice cores used. Our findings show that the optimal definition of the seasons depends on the variability in the target season. For winter, the vigorous variability is best captured when defining the season as December–February due to the dominance of large-scale patterns. For summer, which has weaker variability, albeit more persistent in time, the variability is better captured using a longer season of May–October. Motivated by the scarcity of seasonal data, we also test the use of annual data where the year is divided during summer, that is, not following the calendar year. This means that the winter variability is not split and that the annual data then can be used to reconstruct the winter variability. In particularly when reconstructing the sea level pressure and the corresponding main modes of variability, it is important to take seasonality into account, because of changes in the spatial patterns of the modes throughout the year. Targeting the annual mean sea level pressure for the reconstruction lowers the skill simply due to the seasonal geographical shift of the circulation modes. Our reconstructions based on ice core data also show skill for the North Atlantic sea surface temperatures, in particularly during winter for latitudes higher than 50◦ N. In addition, the main modes of the sea surface temperature variability are qualitatively captured by the reconstructions. When testing the skill of the reconstructions using 19 ice cores compared to the ones using eight ice cores, we do not find a clear advantage of using a larger data set. This could be due to a more even spatial distribution of the eight ice cores. However, including European tree-ring data to further constrain the summer temperature reconstruction clearly improves the skill for this season, which otherwise is more difficult to capture than the winter season.

Description: Russian-German project PLOT (Paleolimnological Transect) aims at investigating the regional responses of the quaternary climate and environment on external forcing and feedback mechanisms along a more than 6000 km long longitudinal transect crossing Northern Eurasia. The well-dated record from Lake El´gygytgyn used as reference site for comparison the local climatic and environmental histories. Seismic surveys and sediment coring up to 54 m below lake floor performed in the frame of the project on Ladoga Lake (North-West of Russia; 2013), Lake Bolshoye Shchuchye (Polar Ural; 2016), Lake Levinson-Lessing and Lake Taymyr (Taymyr Peninsula; 2016-2017), Lake Emanda (Verkhoyansk Range; 2017). Fieldwork at Polar Ural and Taymyr Peninsula was conducted in collaboration with the Russian-Norwegian CHASE (Climate History along the Arctic Seaboard of Eurasia) project. Here, we present the major results of the project obtained so far.

Description: We present a climatology of the near-sea-surface temperature (NSST) anomaly and the sea-ice extent during the Last Glacial Maximum (LGM, 23000–19000 years before present) mapped on a global regular 1◦ × 1◦ grid. It is an extension of the Glacial Atlantic Ocean Mapping (GLAMAP) reconstruction of the Atlantic NSST based on the faunal and floral assemblage data of the Multiproxy Approach for the Reconstruction of the Glacial Ocean Surface (MARGO) project and several recent estimates of the LGM sea-ice extent. Such a gridded climatology is highly useful for the visualization of the LGM climate, calculation of global and regional NSST averages, and estimation of the equilibrium climate sensitivity, as well as a boundary condition for atmospheric general circulation models. The gridding of the sparse NSST reconstruction was done in an optimal way using the Data-Interpolating Variational Analysis (DIVA) software, which takes into account the uncertainty in the reconstruction and includes the calculation of an error field. The resulting Glacial Ocean Map (GLOMAP) confirms the previous findings by the MARGO project regarding longitudinal and meridional NSST differences that were greater than today in all oceans. Taken at face value, the estimated global and tropical cooling would imply an equilibrium climate sensitivity at the lower end of the currently accepted range. However, because of anticipated changes in the seasonality and thermal structure of the upper ocean during the LGM as well as uneven spatial sampling, the estimated cooling and implied climate sensitivity are likely to be biased towards lower values.

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: Reconstruction of past changes in monsoon climate from speleothem oxygen isotope (δ18O) records is complex because δ18O signals can be influenced by multiple factors including changes in precipitation, precipitation recycling over land, temperature at the moisture source, and changes in the moisture source region and transport pathway. Here, we analyse 〉 150 speleothem records of the Speleothem Isotopes Synthesis and AnaLysis (SISAL) database to produce composite regional trends in δ18O in monsoon regions; compositing minimises the influence of site-specific karst and cave processes that can influence individual site records. We compare speleothem δ18O observations with isotope-enabled climate model simulations to investigate the specific climatic factors causing these regional trends. We focus on differences in δ18O signals between the mid-Holocene, the peak of the Last Interglacial (Marine Isotope Stage 5e) and the Last Glacial Maximum as well as on δ18O evolution through the Holocene. Differences in speleothem δ18O between the mid-Holocene and the Last Interglacial in the East Asian and Indian monsoons are small, despite the larger summer insolation values during the Last Interglacial. Last Glacial Maximum δ18O values are significantly less negative than interglacial values. Comparison with simulated glacial–interglacial δ18O shows that changes are principally driven by global shifts in temperature and regional precipitation. Holocene speleothem δ18O records show distinct and coherent regional trends. Trends are similar to summer insolation in India, China and southwestern South America, but they are different in the Indonesian–Australian region. Redundancy analysis shows that 37 % of Holocene variability can be accounted for by latitude and longitude, supporting the differentiation of records into individual monsoon regions. Regression analysis of simulated precipitation δ18O and climate variables show significant relationships between global Holocene monsoon δ18O trends and changes in precipitation, atmospheric circulation and (to a lesser extent) source area temperature, whereas precipitation recycling is non-significant. However, there are differences in regional-scale mechanisms: there are clear relationships between changes in precipitation and δ18O for India, southwestern South America and the Indonesian–Australian regions but not for the East Asian monsoon. Changes in atmospheric circulation contribute to δ18O trends in the East Asian, Indian and Indonesian–Australian monsoons, and a weak source area temperature effect is observed over southern and central America and Asia. Precipitation recycling is influential in southwestern South America and southern Africa. Overall, our analyses show that it is possible to differentiate the impacts of specific climatic mechanisms influencing precipitation δ18O and use this analysis to interpret changes in speleothem δ18O.

Description: The isotopic compositions of oxygen (δ18O) and hydrogen (δ2H) are widely used to locate the geo- graphical origin of biological remains or manufactured products. In this paper, we analyse the distributions of δ18O and δ2H in tap waters sampled across France, and in precipitation interpolated with the Online Isotopes in Precipitation Calculator and modelled with the isotope-enabled ECHAM6-wiso model. Our aim is to provide isoscapes usable in archaeology and forensics and evaluate whether the modelled data could be surrogates for measured ones. The δ18O and δ2H in the 396 tap waters sampled varied spatially within a range of 10‰ and 77‰, respectively. Their consistent distributions followed rules summarized by the effects of altitude and dis- tance from the coast. Their variations along the year were small. Therefore, the database provides a solid ref- erence for δ18O and δ2H of the water supply system at the regional scale. The areas with the most uncommon oxygen and hydrogen isotopic compositions (Atlantic coast south of Brittany and the highest elevations in the Alps) are the most accurately traceable areas in provenancing studies. The isotopic compositions of modelled precipitation have the same spatial distributions but different absolute values from those of tap waters. There- fore, our results favour the use of statistical isoscapes rather than general circulation model-based isoscapes in provenancing studies.

Description: In the context of the Arctic amplification of climate change affecting the regional atmospheric hydrological cycle, it is crucial to characterize the present-day moisture sources of the Arctic. The isotopic composition is an important tool to enhance our understanding of the drivers of the hydrological cycle due to the different molecular characteristics of water stable isotopes during phase change. This study introduces 2 years of continuous in situ water vapour and precipitation isotopic observations conducted since July 2015 in the eastern Siberian Lena delta at the research station on Samoylov Island. The vapour isotopic signals are dominated by variations at seasonal and synoptic timescales. Diurnal variations of the vapour isotopic signals are masked by synoptic variations, indicating low variations of the amplitude of local sources at the diurnal scale in winter, summer and autumn. Low-amplitude diurnal variations in spring may indicate exchange of moisture between the atmosphere and the snow-covered surface. Moisture source diagnostics based on semi-Lagrangian backward trajectories reveal that different air mass origins have contrasting contributions to the moisture budget of the Lena delta region. At the seasonal scale, the distance from the net moisture sources to the arrival site strongly varies. During the coldest months, no contribution from local secondary evaporation is observed. Variations of the vapour isotopic composition during the cold season on the synoptic timescale are strongly related to moisture source regions and variations in atmospheric transport: warm and isotopically enriched moist air is linked to fast transport from the Atlantic sector, while dry and cold air with isotopically depleted moisture is generally associated with air masses moving slowly over northern Eurasia.

Description: The paper presents oxygen and hydrogen isotopes of 284 precipitation event samples systematically collected in Irkutsk, in the Baikal region (southeast Siberia), between June 2011 and April 2017. This is the first high‐resolution dataset of stable isotopes of precipitation from this poorly studied region of continental Asia, which has a high potential for isotope‐based palaeoclimate research. The dataset revealed distinct seasonal variations: relatively high δ18O (up to −4‰) and δD (up to −40‰) values characterize summer air masses, and lighter isotope composition (−41‰ for δ18O and −322‰ for δD) is characteristic of winter precipitation. Our results show that air temperature mainly affects the isotope composition of precipitation, and no significant correlations were obtained for precipitation amount and relative humidity. A new temperature dependence was established for weighted mean monthly precipitation: +0.50‰/°C (r2 = 0.83; p 〈.01; n = 55) for δ18O and +3.8‰/°C (r2 = 0.83, p 〈 0.01; n = 55) for δD. Secondary fractionation processes (e.g., contribution of recycled moisture) were identified mainly in summer from low d excess. Backward trajectories assessed with the Hybrid Single‐Particle Lagrangian Integrated Trajectory (HYSPLIT) model indicate that precipitation with the lowest mean δ18O and δD values reaches Irkutsk in winter related to moisture transport from the Arctic. Precipitation originating from the west/southwest with the heaviest mean isotope composition reaches Irkutsk in summer, thus representing moisture transport across Eurasia. Generally, moisture transport from the west, that is, the Atlantic Ocean predominates throughout the year. A comparison of our new isotope dataset with simulation results using the European Centre/Hamburg version 5 (ECHAM5)‐wiso climate model reveals a good agreement of variations in δ18O (r2 = 0.87; p 〈.01; n = 55) and air temperature (r2 = 0.99; p 〈.01; n = 71). However, the ECHAM5‐wiso model fails to capture observed variations in d excess (r2 = 0.14; p 〈 0.01; n = 55). This disagreement can be partly explained by a model deficit of capturing regional hydrological processes associated with secondary moisture supply in summer.

Description: A long‐term precipitation stable isotope dataset covering every month of the year at least 3 times has been obtained for Irkutsk in southeast Siberia. Distinct seasonal variations are mainly due to changes in air temperature. A new reliable and representative temperature dependency was established for weighted mean monthly precipitation: +0.50‰/°C for δ18O and +3.8‰/°C for δD.

Description: BMBF http://dx.doi.org/10.13039/501100002347

Description: SSHRC http://dx.doi.org/10.13039/501100000155

Description: Ministry of Education and Science of the Russian Federation http://dx.doi.org/10.13039/501100003443

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 50km improves the agreement with the measured δ18O ratios for 14 of 19 observational data sets. A further increase in the spatial resolution to 7km 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 50km. 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.