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.
Introduction
Stable isotopes of water (HD16O and H218Oδ18O) are fractionated during any phase transition. This fractionating process depends on temperature (Dansgaard, 1953; Craig and Gordon, 1965; Jouzel and Merlivat,
1984), and thus water isotopic ratios (expressed here in the usual δ
notation, δD and δ18O with respect to the Vienna
Standard Mean Ocean Water, V-SMOW) reflect the atmospheric conditions under
which the fractionating process took place (Dansgaard, 1964; Merlivat and
Jouzel, 1979; Gat, 1996). This process is generally utilized to reconstruct
paleo-climate conditions such as past temperature changes, using isotopic
ratios stored in climate archives (Dansgaard et al., 1969; Masson-Delmotte
et al., 2005; Jouzel, 2013).
In Arctic regions like Greenland, ice cores constitute an exceptional
climate archive. Over thousands of years, accumulated snow was solidified to
ice, preserving the water isotopic ratios since the last
interglacial period at some locations. Climate reconstructions based on these ice cores show
that the climate conditions changed considerably in Greenland during the
Holocene (here defined as the period between present day and 12 ka; Marcott
et al., 2013). Between the early Holocene and the Holocene Thermal Maximum
in the mid-Holocene (6 ka), a pronounced warm phase took place. Since then,
temperatures steadily decreased until the late Holocene (Marcott et al., 2013; Moossen et al., 2015). In this context, the mid-Holocene is a period
of particular interest, as by that time an Arctic warming had taken place
due to orbital forcing variations and their related feedbacks on large-scale
climate variations, which exhibits similarities to the strong recent Arctic
warming. For Greenland, the mid-Holocene provides the opportunity to
investigate the processes leading to this warming in more detail and to
potentially obtain new insights about the future development of the Arctic
region (Yoshimori and Suzuki, 2019).
While general circulation models (GCMs) are generally able to reproduce the
direction and large-scale patterns of past climate changes (e.g., Timm and
Timmermann, 2007; Smith and Gregory, 2012), they often fail to reproduce the
magnitude of regional changes (Braconnot et al., 2012; Harrison et al., 2014) documented in various local climate archives. Thus, a scale gap might
exist between the measured point information and the large-scale climate
information generated by GCMs. The comparison of observational and GCM data
can therefore be subject to considerable uncertainties (Felzer and Thompson,
2001).
Especially for structured landscapes, the spatial resolution in GCMs is
often too coarse to resolve relevant local factors (Jost et al., 2005;
Fischer and Jungclaus, 2011). Important properties like topography and
surface conditions are consequently only represented in a generalized and
imprecise form in climate simulations. In most cases, this does not
adequately represent the complex characteristics of the land surface and its
associated interactions with the atmosphere. For stable water isotopes, key
physical processes of isotope fractionation are therefore not well resolved
in coarse resolution GCMs, leading to differences between simulated and
observational isotope data, especially in complex terrains (Sturm et al., 2005; Werner et al., 2011). Isotope-enabled GCMs are consequently not able
to quantitatively reproduce regional changes in isotope ratios (e.g., Risi et
al., 2010), and the simulated isotope ratios with GCMs exhibit in many cases
larger deviations relative to observed ratios than the results of
corresponding regional climate model (RCM) simulations. For instance, Sturm
et al. (2007) were able to reduce the bias of simulated isotope ratios in
precipitation through a regional downscaling of an isotope-enabled GCM run
in South America. Comparable results were achieved by Sjolte et al. (2011)
for isotope-enabled RCM simulations in Greenland.
Therefore, in the present study, isotope-enabled GCM simulation results for
the Arctic region are dynamically downscaled with an isotope-enabled RCM to
a higher temporal and spatial resolution. By means of such regional
simulations, the spatial and temporal variability of the isotopic ratios in
the Arctic is potentially increased, accounting for the heterogeneity of
local conditions at the different ice core locations and the associated
uncertainties. In this way, the impact of highly resolved local conditions
on the spatial and temporal variability of isotopic ratios is investigated,
and the impact of such small-scale variability on the discrepancy between
simulated and observed paleo-climate conditions in the Arctic region is
examined.
To explore this, the isotope-enabled version of the RCM COSMO-CLM (Rockel et
al., 2008), COSMO_iso (Pfahl et al., 2012; Christner et al., 2018), is used. In a first step, the general suitability of
COSMO_iso to be used for isotope applications in Greenland is
assessed. For this purpose, near-surface temperatures and precipitation
amounts simulated with the standard COSMO version are compared with
observations in the Arctic region. Subsequently, the capability of
COSMO_iso to simulate realistic water isotopic ratios for
Greenland is tested by downscaling a global present-day simulation with an
isotope-enabled GCM for the Arctic region. The GCM and RCM results are then
compared to measured water isotope ratios in precipitation and snow pit
samples. Afterwards, the tested isotope-enabled COSMO_iso
model system is used to downscale an isotope-enabled GCM simulation for a
mid-Holocene time slice. The simulated isotopic ratios are evaluated against
Greenland ice core data. Such a dynamical downscaling of global isotope
simulations for Greenland under mid-Holocene conditions is performed for the
first time in the framework of this study.
MethodsCOSMO_isoModel description
In this study, simulated stable water isotope concentrations of HD16O
and H218Oδ18O with isotope-enabled GCMs (Sect. 2.1.2), are regionally downscaled with COSMO_iso (Pfahl et al., 2012), an
isotope-enabled version of the numerical weather prediction model COSMO
(Consortium for Small-scale Modeling; Baldauf et al., 2011) (version 4.18).
For the purpose of long-term climate simulations, isotope routines of
COSMO_iso were implemented in COSMO-CLM (Rockel et al., 2008), the climate version of COSMO. In this context, the δD and
δ18O ratios in the soil water and the surface-layer snow are
simulated with TERRA_iso version 1 (Dütsch, 2017; Christner et al., 2018), the isotope-enabled version of the multi-layer Land Surface Model TERRA-ML (Schrodin and Heise, 2001) in COSMO. In several studies, COSMO_iso and TERRA_iso were successfully employed for the simulation of isotopic ratios in the midlatitudes (Pfahl et al., 2012; Aemisegger et al., 2015; Christner et al., 2018). In the present study, the model system will be applied to the
Arctic region. For this, some additional modifications regarding the
treatment of snow and ice had to be implemented in the model.
Snow albedo
The surface albedo of fresh snow is increased from 0.7 to 0.8 to improve the
model agreement with measured values of short-wave reflectance and 2 m
temperature at stations from the Cooperative Institute for Research in
Environmental Sciences at the University of Colorado Boulder (CIRES) in
central Greenland (Karremann and Schädler, 2021).
Snow layer thickness
In the standard configuration of COSMO, the Greenland ice sheet is treated
as a constant mass of ice that is covered by a single snow layer. But in
this model structure, dynamical processes within the ice sheet (flow, basal
melt) are not included. As a result, the depth of the snow layer is constantly increasing and thus so is its heat capacity. To avoid this
spurious model behavior, the snow layer depth is limited to 5 cm in this
study. Using this value, realistic diurnal cycles of the 2 m air temperature
could be simulated.
Marine regions with sea ice cover
To be able to simulate reasonable fractionation processes for marine regions
with sea ice cover, a snow layer is also implemented on top of the sea ice
(e.g., as suggested in Bonne et al., 2019). The isotopic composition of this
surface snow layer is in this case set to the isotopic composition of the
most recent precipitation.
Fractionation at snow covered surfaces
Isotope fractionation during sublimation from a surface snow layer is poorly
understood. Several different processes are suggested to be involved, which
are not yet taken into account in state-of-the-art isotope-enabled models
(see, e.g., discussion in Christner et al., 2017), such as non-fractionating
layer-by-layer sublimation (e.g., Ambach et al., 1968), kinetic fractionation
during sublimation into sub-saturated air, a diurnal cycle of sublimation
combined with fractionating vapor deposition on the snow (e.g., Steen-Larsen
et al., 2014) and fractionating melt water evaporation combined with
recrystallization of residual melt water (Gurney and Lawrence, 2004). To
approximate this complex interplay of different influencing factors, in this
study an equilibrium fractionation during sublimation from surface layer
snow and sea ice is assumed. However, the authors are aware that this is a
simplified description of isotope fractionation during sublimation.
Model simulation setup
The capability of COSMO_iso to realistically reproduce the
fractionating processes of stable water isotopes in Greenland is evaluated.
For this, the nudged simulation outputs (standard and isotopic) from an
isotope-enabled atmospheric model ECHAM5-wiso (Werner et al., 2011)
simulation are dynamically downscaled with COSMO_iso for the
whole Arctic region. The data from the same ECHAM5-wiso simulation have been
already used as boundary conditions for COSMO_iso simulations
over Europe by Christner et al. (2018). The simulation outputs from
ECHAM5-wiso are at a T106 horizontal spatial resolution (1.1∘×1.1∘) and on 31 atmospheric vertical levels. The dynamical fields
were nudged every 6 h towards ERA-Interim reanalysis data (Dee et al., 2011). Monthly varying sea surface temperatures and sea ice cover were
prescribed as lower boundaries over sea, also based on the ERA-Interim data.
The simulation period is 2008–2014. In order to guarantee a physically consistent transition between the
coarse model resolution of the GCM and the fine model resolution of the RCM at the boundaries
of the regional model domain,
the model resolution is increased stepwise. This procedure is called
nesting. In a first nesting step, the spatial resolution of
COSMO_iso is set to 0.44∘×0.44∘, corresponding to 50km×50km in rotated coordinates (COSMO_iso_50km). In a second nesting step, an additional
COSMO_iso simulation with a spatial resolution of
0.0625∘×0.0625∘ (corresponding to about 7km×7km) for Greenland (COSMO_iso_7km)
is nested in the COSMO_iso_50km simulation. This high-resolution simulation covers the year 2011. In the COSMO_iso runs, the horizontal wind fields above the 850 hPa level are spectrally nudged (von Storch et al., 2000) towards the
reanalysis-based dynamical fields of ECHAM5-wiso. This method ensures that
consistent atmospheric boundary conditions build the framework for the
fractionating processes simulated in COSMO_iso. The model
domains of both simulations is shown in Fig. 1.
COSMO_iso model domain of the 50 km (blue) and the
7 km (black) simulation.
The same model chain is applied to the mid-Holocene period. Atmospheric
fields have been retrieved from a mid-Holocene simulation of the fully coupled isotope-enabled Max Planck Institute Earth System Model
(MPI-ESM-wiso, Cauquoin et al., 2019), whose atmospheric component is
ECHAM6-wiso. The major ECHAM6 model changes compared to ECHAM5 include an
improved representation of radiative transfer in the solar part of the
spectrum, an improved representation of surface albedo, a new aerosol
climatology and an improved representation of the middle atmosphere (Stevens
et al., 2013). The ocean component is the Max Planck Institute Ocean Model
(MPIOM, Jungclaus et al., 2013). With COSMO_iso, only a representative time slice of 30 years is simulated for this climate period, as the regional COSMO_iso simulations are computationally very expensive. The greenhouse gas concentrations and the orbital parameters are adapted according to the Paleoclimate Modelling Intercomparison Project 4 experiment design (PMIP4, Kageyama et al., 2018). The model domain of the COSMO_iso simulations is identical to the present-day simulations.
Observations
The capability of the isotope-enabled regional climate model COSMO_iso to reproduce measured isotopic ratios in Greenland is evaluated by comparing the simulation results to observational data. The simulated isotopic composition in precipitation is assessed by comparing the model results in the Arctic region with observed monthly data from the Global Network of Isotopes in Precipitation (GNIP) of the International Atomic Energy Agency and the World Meteorology Organization (IAEA, 2018) over the period 2008–2014 (Table 1). Furthermore, simulated δ18O ratios are compared to snow pit samples collected during the North Greenland Traverse (Fischer et al., 1998; Weißbach et al., 2016a) and top core samples from five ice core locations (Renland, Vinther et al., 2008; Neem, Masson-Delmotte et al., 2015; GISP2, Grootes and Stuiver, 1997; Summit, Fischer, 2003; SE-Dome, Furukawa et al., 2017). The station numbers assigned to the respective samples within this study, as well as their locations and δ18O values are summarized in Table 2. Since all snow pit samples cover different time periods, the present-day δ18O values (black numbers in Table 2) are calculated as an average of all available δ18O values measured between 1940 and 2014. With this procedure, uncertainties in snow pit samples and top ice core samples, associated with post-depositional diffusion and wind erosion and the resulting constraints in analyzing annual and interannual top ice core data (e.g., Johnsen et al., 2000), can be neglected. However, further uncertainties in snow pit samples and ice core data remain regarding the timescale assignment (Steig et al., 2005) and the spatial variability (Weißbach et al., 2016b).
Description of the snow pit and ice core samples used in this study. The present-day δ18O values are calculated as an average of all available δ18O values measured in snow pit samples between 1940–2014. The mid-Holocene δ18O values are calculated as an average of the measured δ18O values in ice cores over the period 5.5–6.5 ka. Black numbers indicate present-day δ18O values and the respective standard deviations, and bolded numbers indicate mid-Holocene values and the corresponding standard deviations.
No.NameSampleLongitudeLatitudeδ18OσReference1Renlandtop core-26.7371.27-27.38 (-26.44)0.42 (0.31)Vinther et al. (2008)2NEEMtop core-51.0677.45-33.241.15Masson-Delmotte et al. (2015)3GISP2top core-38.4872.58-34.95 (-34.83)0.69 (0.33)Grootes and Stuiver (1997)4Summittop core-37.6473.03-36.461.11Fischer (2003)5B27_B28snow pit-46.4876.65-34.051.52Weißbach et al. (2016a)6NGT03C93snow pit-37.6273.94-37.021.06Weißbach et al. (2016a)7NGT06C93snow pit-37.6275.25-36.891.38Weißbach et al. (2016a)8NGT14C93snow pit-36.476.61-36.181.64Weißbach et al. (2016a)9NGT23C94snow pit-36.578.83-35.181.78Weißbach et al. (2016a)10NGT27C94snow pit-41.1380-34.011.55Weißbach et al. (2016a)11NGT30C94snow pit-45.9179.34-34.191.69Weißbach et al. (2016a)12NGT33C94snow pit-4478-36.131.62Weißbach et al. (2016a)13NGT37C95snow pit-49.2177.25-33.811.38Weißbach et al. (2016a)14NGT39C95snow pit-46.4876.65-34.951.4Weißbach et al. (2016a)15NGT42C95snow pit-43.4976-35.531.24Weißbach et al. (2016a)16NGT45C95snow pit-4275-35.331.32Weißbach et al. (2016a)17GRIPtop core-37.6472.58-34.730.32Vinther et al. (2006)18NGRIPtop core-42.3275.1-34.690.36Vinther et al. (2006)
Since both snow pit samples and top core samples from ice cores represent an
integrated signal of the isotopic composition in precipitation, we computed
modeled annual mean δ18O values in precipitation and compared
the multi-year 2008–2014 model mean to the observed values. For the
calculation of the yearly modeled mean values, the modeled δ18O in precipitation is weighted with accumulation rate, i.e., months with high precipitation amounts get a higher weight compared to months with small precipitation amounts.
Ice core samples are also used to evaluate the simulated isotopic ratios for
the mid-Holocene. Besides the already-mentioned Renland and GISP2 samples,
two more ice core samples, namely GRIP and NGRIP (Vinther et al., 2006), are
used for the model evaluation. The mid-Holocene δ18O values
(blue numbers in Table 2) are calculated as an average of the measured
δ18O values in ice cores over the period 5.5–6.5 ka.
In a first step, the general capability of the COSMO model to reproduce
observed standard climatological parameters in present-day simulations for
Greenland is assessed. For this purpose, the results of an ERA-Interim
reanalysis-driven (Dee et al., 2011) simulation with the standard COSMO
model (without isotope application) are compared with observations
collected by the Danish Meteorological Institute (DMI). Evaluated are the
yearly mean 2 m temperatures (Fig. 2a) and the yearly mean precipitation
sums (Fig. 2b) over the period 1995–2015. Both simulated 2 m temperatures and simulated precipitation amounts are in good agreement with the DMI observations. The simulated 2 m temperatures coincide especially well with the observed values. For the precipitation sums, the spread of simulated and observed values is higher than for the 2 m temperatures, a feature generally occurring in weather and climate simulations. Thus, the COSMO model is generally able to simulate reasonable climate conditions in Greenland and can therefore be used for isotope applications in this region. A detailed analysis of the standard COSMO model performance in the Arctic region is presented in Karremann and Schädler (2021).
Simulated yearly mean (a) 2 m temperatures and (b) precipitation sums of a standard COSMO simulation, driven with ERA-Interim, for Greenland over the period 1995–2015 compared to DMI observations. The solid black line is the 1:1 line in (a, b).
Comparison of simulated δ18O data to station data
Figure 3 shows the yearly mean δ18O values for the period
2008–2014 for Greenland, simulated with COSMO_iso_50km (a) and ECHAM5-wiso (b). Additionally, the locations and the observed δ18O values of the 19 snow pit samples used to assess the models' capability to reproduce observed δ18O ratios in Greenland are illustrated. In general, COSMO_iso in a 50km×50km spatial resolution is able to reproduce the observed isotopic ratios at the snow pit samples and improves the simulation results of ECHAM5-wiso. In both simulations, the δ18O ratios are high near the coastline and low in central Greenland. But in COSMO_iso_50km, the δ18O ratios decline more rapidly from the coastline to the inland plateau than in ECHAM5-wiso. The spatial δ18O differences are consequently more pronounced, and the general overestimation of δ18O ratios, which occurs in ECHAM5-wiso, is reduced in COSMO_iso_50km. As
a consequence, the regional simulation reaches a better agreement with the
observations. In COSMO_iso_50km the root-mean-squared error (RMSE) is reduced by 0.98 ‰ over all snow pit samples (the RMSE of ECHAM5-wiso is 2.42 ‰, and the RMSE of COSMO_iso_50km is 1.44 ‰). This RMSE reduction is significant at the 95 % level, assessed by performing a t test. Especially for the snow pit samples for which ECHAM5-wiso exhibits strong deviations from the observed δ18O values (1,3,4,6,7,8,9,10,16,17,19; see Table 2), a regional downscaling with COSMO_iso_50km reduces the bias considerably (Fig. 4). For these stations, the RMSE of ECHAM5-wiso of 3.09 ‰ is reduced by 1.65 ‰ to 1.44 ‰. But for snow pit samples at which ECHAM5-wiso already has a high agreement with the observations (2,5,11,13,14), COSMO_iso_50km increases the RMSE from 0.34 ‰ to 1.51 ‰.
Yearly mean δ18O values of COSMO_iso_50km (a) and ECHAM5-wiso (b, interpolated to the COSMO_iso_50km grid) for the period 2008–2014 and the corresponding observations for the 19 snow
pit samples (Table 2).
Differences (Δ) between simulated and observed δ18O values (model minus observation) for the model results of ECHAM5-wiso, COSMO_iso_50km, and COSMO_iso_7km and snow pit samples and top core
samples from ice cores from Greenland (averaging periods: COSMO_iso_50km: 2008–2014; COSMO_iso_7km: 2011; observation: 1940–2014). Numbers refer to the different snow pit locations shown in Fig. 3.
Figure 5 shows that these annual biases of the COSMO_iso_50km simulation are not caused by systematic seasonal biases, as for example reported by Sjolte et al. (2011) for RCM simulations in Greenland. Shown are the simulated monthly δ18O values with COSMO_iso_50km compared to observed monthly δ18O values in precipitation for the period 2008–2014, collected at Arctic stations of the GNIP data set (Table 1). In general, the modeled δ18O values for precipitation are in good agreement with the monthly GNIP data. But in the COSMO_iso_50km simulation no systematic overestimation or underestimation of
observed isotope ratios is simulated with the RCM. This is true for each
season. Neither in winter (low δ18O values) nor in summer
(high δ18O values) are systematic deviations to the
observations simulated. Thus, the seasonal variability in the
COSMO_iso_50km results has no systematic impact on the yearly mean δ18O values and is therefore not the reason for systematic differences between the coarse model results and observations.
Monthly δ18O simulated with
COSMO_iso_50km for the period 2008–2014
versus the corresponding observations for nine GNIP stations (Table 1).
As is visible in Fig. 3, these systematic differences are instead caused by a
southward shift of the area of low yearly mean δ18O values in
central northern Greenland in COSMO_iso_50km
relative to ECHAM5-wiso. As a result, the simulated δ18O values in central northern Greenland in COSMO_iso_50km are higher than in ECHAM5-wiso. Since ECHAM5-wiso already has a high agreement with the observed δ18O values, a model bias is introduced in COSMO_iso_50km, causing the deviations relative to the observations in northern Greenland. But COSMO_iso_50km in general yields an overall
improvement in simulating the yearly mean δ18O values compared
to ECHAM5-wiso.
A further downscaling with COSMO_iso to a spatial resolution of 7km×7km does not improve the simulation results further. The RMSE reduction with respect to ECHAM5-wiso is 0.83 ‰
(compared to 0.98 ‰ for COSMO_iso_50km). The only exception constitutes the snow pit sample from Renland (1). Here, a considerable model bias in ECHAM5-wiso and COSMO_iso_50km is strongly reduced in
COSMO_iso_7km. The coastal area of Renland is
characterized by complex terrain and constitutes a special case for
isotope-enabled modeling in Greenland. The snow pit sample is located in a
transition zone from the homogeneous inland glaciation to the rugged
coastline, where the glaciers calve into the sea. Thus, within short
distances large differences in altitude and land surface characteristics
occur in this region. The isotopic ratios in the snow pit sample are
therefore strongly affected by these heterogeneous local conditions, which
are insufficiently represented in the coarse model resolution of
ECHAM5-wiso. By increasing the spatial resolution with regional climate
modeling, the representation of the associated small-scale processes is
improved. This leads generally to an improved agreement of the simulation
results with observations, as seen for the COSMO_iso_7km run for Renland (Fig. 4).
However, an increase in spatial resolution is also associated with an
increased heterogeneity of the surface characteristics and the related
small-scale processes, especially in complex terrains. This is because the
GCM grid boxes are further divided in smaller RCM grid boxes and
consequently both higher and lower values (e.g., altitude) are now
included in the respective GCM grid boxes. As a consequence, an additional
spatial variability is introduced in the RCM simulations in comparison to
the GCM results. Due to uncertainties accompanied by model simulations, this
can potentially increase the RCM bias with respect to in situ point
measurements, which may actually be closer to the spatially averaged values
simulated by the coarse GCM model. This effect can be observed for the
SE-Dome ice core (19) in southeastern Greenland. Comparable to the Renland
ice core, SE-Dome is located near the coastline. But in contrast to the
Renland ice core, an increase in the spatial resolution to 7 km does not
further improve the RCM results for SE-Dome. On the contrary, the δ18O bias is even higher than in the ECHAM5-wiso simulation.
Despite the lack of improvement in the point to grid cell comparison, higher
resolved RCM simulations allow the subgrid-scale variability of δ18O within GCM grid boxes to be simulated and compared to observed δ18O values. In this way, the inherent uncertainty of in situ measurements, associated with a local micrometeorological variability, can be considered. Thus, in the following sections, snow pit samples are no longer solely compared to the model grid boxes covering the sample locations. Instead, it is investigated whether the δ18O range of all adjacent RCM grid boxes to a snow pit location is consistent with the observed δ18O value of the same site. For this, all RCM grid boxes located within the corresponding GCM grid box are included in the comparison with the observations.
δ18O variability
The spatial isotopic ratio variability of the COSMO_iso_50km grid boxes surrounding the 19 snow pit samples is shown as a box and whisker plot in Fig. 6a. The spatial isotopic ratio variability of the COSMO_iso_7km is shown in Fig. 6b. This spatial isotopic ratio variability includes the δ18O values of all COSMO_iso (50 km) and COSMO_iso (7 km) grid boxes within the ECHAM5_wiso grid box closest to the snow pit sample. For 14 of the 19 snow pit samples (1,2,3,4,5,6,7,8,14,15,16,17,18,19) the observed δ18O
values are within the range of the spatial COSMO_iso_50km grid box variability. But for 5 of the 19 snow pit samples (9–13) the spatial isotope range of the COSMO_iso_50km simulation does not fit with the observations. Since these stations are all located in the north of Greenland (Fig. 3), this is likely partially associated with the southward shift of the area of low yearly mean δ18O values in central northern Greenland in COSMO_iso_50km in comparison to ECHAM5-wiso, as already described in Sect. 3.1.2.
Present-day isotopic ratio variability of the COSMO_iso grid boxes surrounding the 16 snow pit samples for the (a) 50 km and (b) 7 km simulation. The black bar in the box and whisker plot represents the median of the isotope ratio distribution. The box comprises the upper and lower quartile, and the whiskers show the whole distribution. The MPI-ESM-wiso results are shown by the blue dots, and the observed δ18O values are shown by the red dots.
A downscaling to 7 km does slightly increase the spread of the
COSMO_iso results. However, the observed δ18O values from 5 of 19 snow pit samples still fall outside the range of the modeled COSMO_iso_7km grid box variability (Fig. 6b). Thus, a further downscaling to a spatial resolution of 7 km does not substantially change the simulated isotopic ratio spread within an ECHAM5-wiso grid box. In accordance with the missing benefits of the COSMO_iso_7km simulation and its increased computing time costs, only a COSMO_iso_50km simulation is performed for the mid-Holocene (Sect. 3.2).
The high spatial δ18O variability in COSMO_iso
simulations is also reflected in the spatial δ18O–temperature
slope of the COSMO_iso_50km run (Fig. 7c), a measure that is frequently used to analyze how strongly isotope ratios and surface temperatures are interrelated. The spatial isotope–temperature slope constitutes a linear fit between the simulated δ18O ratios and the surface temperatures at all COSMO_iso_50km grid boxes within the respective ECHAM5-wiso grid box. The spatial δ18O–temperature slope is high in central Greenland (1.5–2.0) and moderate at the coastline (∼1.0). In order to better understand these spatial δ18O–temperature patterns, both quantities affecting the spatial δ18O–temperature slope, i.e., the spatial δ18O variability and the spatial temperature
variability, are explicitly shown in Fig. 7a and b. These spatial
variabilities are calculated as the standard deviation of all COSMO_iso_50km grid boxes within the respective ECHAM5-wiso grid boxes.
Present-day spatial subgrid-scale variability (calculated as standard deviation) of (a)δ18O and (b) surface temperature, derived from the COSMO_iso_50km grid boxes within the respective ECHAM5-wiso grid boxes for all of Greenland. Present-day (c) spatial and (d) interannual temporal δ18O–temperature slope for Greenland, based on yearly mean values.
Figure 7a shows that at the coastline, the spatial δ18O
variability of COSMO_iso_50km is considerably increased within the ECHAM5-wiso grid boxes. In central Greenland the spatial isotopic ratio variability is lower. Thus, the simulated spatial δ18O variability is high in regions where large orographic differences occur within short distances, like the coastal areas of
Greenland, and lower for homogeneous terrain like the inland plateau.
Nevertheless, widespread areas with higher spatial isotopic variability also occur in the inland plateau of Greenland. This is not the case for the
spatial surface temperature variability. In central Greenland, surface
temperature variability is very low (Fig. 7b), but near the coastline the
spatial surface temperature variability is also high, highlighting how
important the land surface characteristics are for the regional temperature
variability. The moderate spatial δ18O–temperature slope at the coastline (Fig. 7c) is therefore a result of a high surface temperature
variability in this region counteracting the high δ18O
variability in the slope calculation. There, the correlation between both
quantities is consequently high (Fig. S1a in the Supplement, mainly between 0.7 and 0.99). In central Greenland, the spatial δ18O–temperature slope is further increased due to the relatively high δ18O variability compared with the surface temperature variability. Therefore, the spatial distribution of δ18O cannot be solely explained by land surface processes and the associated spatial temperature variability. The additional spatial δ18O variability consequently must be caused by dynamic atmospheric processes. In this way, isotopic ratios based on atmospheric fractionation processes along the trajectory of an air mass are transported to central Greenland and increase the isotopic variability there.
In order to investigate the temporal interrelations between the isotope
ratios and the surface temperature, the interannual temporal δ18O–temperature slope is calculated for the COSMO_iso_50km simulation, based on the yearly mean δ18O and surface temperature values (Fig. 7d). The interannual temporal δ18O–temperature slope is, in contrast to the spatial
δ18O–temperature slope, small all over Greenland, which is in
agreement with the results of Sjolte et al. (2011). That means that the
interannual δ18O variability is less pronounced than the
interannual surface temperature variability, and thus the correlation
between both quantities is also not as strong as for the spatial δ18O variability and the spatial temperature variability (Fig. S1b). The mean spatial correlation over Greenland is 0.61, compared to 0.25 for the mean temporal correlation. The impact of interannual surface temperature variations on the temporal δ18O variability in Greenland is therefore not as dominant as for the spatial δ18O variability.
Mid-HoloceneComparison of simulated δ18O data to ice core data
In contrast to the present-day simulations, for the mid-Holocene
COSMO_iso_50km is driven by MPI-ESM-wiso rather than ECHAM5-wiso. While in ECHAM5-wiso oceanic boundary conditions are prescribed by monthly varying sea surface temperatures and sea ice cover, ocean states are calculated internally in the fully coupled atmosphere–ocean Earth system model MPI-ESM-wiso. Systematic deviations between the COSMO_iso_50km simulations for the mid-Holocene and present day, caused by these different forcing approaches, therefore cannot be excluded. For this reason, a comparison of the mid-Holocene δ18O anomalies to the present-day conditions is omitted, and an analysis is performed for simulated mid-Holocene δ18O ratios with comparison to observed mid-Holocene δ18O values.
In Fig. 8a the absolute differences between the simulated MPI-ESM-wiso
(blue) and COSMO_iso_50km (green) grid box results and the observed δ18O ratios at the corresponding ice cores are presented for the mid-Holocene. As in Fig. 6, the spatial isotopic ratio variability of the COSMO_iso_50km grid boxes surrounding the four Greenland ice core samples is shown as a box and whisker plot. MPI-ESM-wiso properly reflects the isotopic ratios of the mid-Holocene from ice core data. For the inland ice cores (GRIP, GISP2, NGRIP), the simulated δ18O deviates only by about 1 ‰ relative to the observation, and at Renland the
deviation is about 3 ‰. For GRIP and GISP2 the MPI-ESM-wiso simulations slightly underestimate the δ18O ratios, and for NGRIP and Renland the δ18O values are slightly overestimated.
(a) Mid-Holocene isotopic ratio variability of the COSMO_iso_50km grid boxes surrounding four Greenland ice core samples. In each grid box, the observed ratios derived from the ice cores are subtracted from the simulated δ18O ratios. The black bar in the box and whisker plot represents the median of the isotope ratio difference distribution. The box comprises the upper and lower quartile, and the whiskers show the whole distribution. The MPI-ESM-wiso (blue dots) and COSMO_iso_50km (green dots) results for the grid points closest to the ice cores are also shown as differences with respect to the observed δ18O ratios. (b) Anomalies of the MPI-ESM-wiso simulation relative to pre-industrial (PI) conditions, based on an MPI-ESM-wiso PI reference simulation (Cauquoin et al., 2019) are shown as blue points, and the observed anomalies for the mid-Holocene relative to present-day are shown as red points.
For COSMO_iso_50km, the deviation of δ18O values relative to the observations are opposite in sign compared with MPI-ESM-wiso at all locations except Renland. That means that in GRIP and GISP2, the underestimated δ18O values in MPI-ESM-wiso are turned into overestimated δ18O values in COSMO_iso_50km; at NGRIP the overestimation is turned into an underestimation, but the net bias is not reduced. At Renland, the bias is even increased. Thus, by just looking at the absolute biases, the downscaling does not seem to bring an added value to the MPI-ESM-wiso results for mid-Holocene conditions. However, when the spatial isotopic ratio variability within MPI-ESM-wiso grid cells simulated by
COSMO_iso_50km isotopic ratios is taken into account, the model results are in agreement with the isotopic ratios of the ice core samples.
In Fig. 8b, the MPI-ESM-wiso model anomalies with reference to the pre-industrial period (PI) conditions, which are based on an MPI-ESM-wiso
PI reference simulation performed by Cauquoin et al. (2019), are compared to
the measured mid-Holocene-PI δ18O anomalies of the ice cores.
The positive δ18O anomalies between the mid-Holocene and PI values for
both ice core data and MPI-ESM-wiso model results are associated with higher
temperatures, especially during the summer, and a reduction in Arctic
sea ice during the mid-Holocene (Cauquoin et al., 2019). In Renland and NGRIP
simulated anomalies are slightly underestimated, and in GRIP and GISP2
anomalies are slightly overestimated, but overall the biases of the
MPI-ESM-wiso mid-Holocene–PI period model anomalies to the observed mid-Holocene–PI period
anomalies are small and statistically not significant at
the 95 % level for all ice cores (assessed by performing a t test).
δ18O variability
The fact that, in contrast to the present-day simulations, only four
observational data sets are available for the mid-Holocene, makes the
assessment of the simulation results difficult. Moreover, with the GRIP and GISP2
ice cores being located very close to each other (Fig. 9), only three
local isotope distributions that are clearly different from each other are available.
Therefore, in Fig. 9a, the spatial δ18O variability of the
COSMO_iso_50km simulation for the mid-Holocene is illustrated for all of Greenland, which is, in accordance to the analysis of the present-day simulation, again calculated as the standard deviation of all COSMO_iso_50km grid boxes within the respective GCM grid boxes. In general, the δ18O variability of COSMO_iso_50km in the mid-Holocene is high at the coastline, while it is lower in central Greenland. The Renland ice core is consequently located in an area of a high isotopic variability, and the GRIP and GISP2 ice cores are located in an area of low isotopic variability. But regions with increased isotopic variability also occur in the inland plateau of Greenland. The NGRIP ice core, for instance, is located in such an area of a moderate isotopic ratio variability. The four ice core drill sites are therefore located in three regions of Greenland with substantially different sub-grid isotopic ratio variabilities.
The same as Fig. 7 but for the mid-Holocene. The locations of
the ice core samples are shown in green.
The spatial surface temperature variability in the COSMO_iso_50km mid-Holocene simulation is shown in Fig. 9b. The mid-Holocene simulation shows a high spatial surface temperature variability near the coastline and almost no variability in central Greenland. As a consequence, the spatial δ18O–temperature slope is moderate at the coastline (∼1.0) and high in central Greenland (1.5–2.0; Fig. 9c). Moreover, the interannual δ18O–temperature slope is very small over Greenland in the mid-Holocene (the mean interannual slope over Greenland is 0.47), although in some regions high temporal slopes are simulated (Fig. 9d). However, the influence of interannual surface temperature variations on the temporal δ18O variability in the mid-Holocene is, in general, not as dominant as for the spatial variability.
In general, the results of the COSMO_iso_50km mid-Holocene simulation exhibit the same spatial characteristics as for the present-day simulation (Figs. 7 and 9). Comparable spatial patterns are simulated for the surface temperature variability (Figs. 7b and 9b) and the δ18O variability (Figs. 7a and 9a) within a GCM grid box, although regions of increased δ18O variability in central Greenland are more widely present in the mid-Holocene run than in the present-day one. In the inland plateau, regions of a high spatial δ18O–temperature slope are therefore more extensive in the mid-Holocene than under present-day conditions (Figs. 7c and 9c). Nevertheless, the spatial δ18O–temperature interrelations are comparable in both periods. This is also the case for the temporal variabilities of δ18O and the surface temperature in the COSMO_iso_50km simulation results for the mid-Holocene and the present day (Figs. 7d and 9d), although both simulations are driven by two different forcing approaches (ECHAM5-wiso nudged to ERA-Interim reanalysis with prescribed monthly varying oceanic boundary conditions vs. the fully-coupled atmosphere–ocean Earth system model MPI-ESM-wiso). This finding indicates that the spatial and interannual δ18O variability of COSMO_iso_50km within a GCM grid box over Greenland is not strongly
dependent on the oceanic boundary conditions.
Discussion and conclusions
The results of several global paleo-climate simulations exhibit considerable deviations from the observed regional climate patterns during the Holocene (Braconnot et al., 2012). In the present study, regional climate simulations with an isotope-enabled RCM are, for the first time, performed for Greenland to potentially improve the agreement with climate observations in this region for the mid-Holocene. In a first step, the capability of the isotope-enabled RCM COSMO_iso to reproduce observed isotopic ratios for Greenland is demonstrated.
The COSMO_iso simulation results show that a spatial
resolution of 50 km produces reasonable δ18O values. Especially
in regions where the global ECHAM5-wiso model, which has been used to derive
necessary forcing fields for the COSMO_iso simulations,
deviates strongly from the observed δ18O values, the RMSE is
significantly reduced by 1.65 ‰ for regional climate simulations with COSMO_iso. In complex terrain like the coastal areas of Greenland, the results can be further improved with an additional downscaling to a spatial resolution of 7 km. In such simulations with high spatial resolution, small-scale processes are described in more detail (e.g., Torma et al., 2015; Coppola et al., 2020), and thus the local characteristics at ice core sites are better taken into account (Sturm et al., 2005; Werner et al., 2011). But for the rest of Greenland, highly resolved regional climate simulations do not yield further improvements. For northern Greenland, regional climate simulations with COSMO_iso increase the bias with respect to observations. A comparison of simulated isotope ratios in precipitation with measured values at GNIP stations shows that such deviations between model results and observations are not caused by systematic seasonal biases in the RCM, as it was the case in a similar study by Sjolte et al. (2011) for Greenland. In central northern Greenland, a model bias is instead introduced, due to a southward shift of the area of low yearly mean δ18O values. All in all, the results of this study show that COSMO_iso is generally
able to provide reasonable isotopic ratios for Greenland and that the model can be applied for paleo-climate simulations.
For the mid-Holocene, MPI-ESM-wiso is in good agreement with observed ice
core data in Greenland, as already described by Cauquoin et al. (2019). The
model bias is, in this context, not further reduced by a downscaling with
COSMO_iso. But an increase in the spatial model resolution also leads to an increase in the models' degrees of freedom. This in turn can lead to additional noise and thus a deviating RCM behavior and even an increase in the absolute model bias, as seen for the Renland station.
Another consequence of these increased degrees of freedom in the
COSMO_iso simulation is that the spatial variability of the
simulated δ18O ratios is enhanced. This enhanced spatial
variability represents the subgrid-scale uncertainty of the driving GCM,
which can be derived in a physically consistent way by regional
downscaling. By analyzing this subgrid-scale variability, the spatial
uncertainties in the comparison between GCM data and point measurements can
be considered. In this way, it can be demonstrated that most of the observed
δ18O values lie within the local δ18O
uncertainties of the coarse GCM results. This applies for both the
present-day runs and the regional paleo-climate simulations for the
mid-Holocene in Greenland. The deviation between the more coarsely resolved GCM
results and the more finely resolved observations is therefore potentially caused
by the missing representation of important small-scale processes, which are
induced by, e.g., the surface conditions or orographic effects over Greenland.
Shi et al. (2020), for instance, were able to demonstrate that GCM
deficiencies in reproducing the observed water isotope variability in the
southeastern Tibetan Plateau are associated with the missing representation
of such small-scale processes in coarse GCM simulations.
As δ18O ratios are used as an indicator for temperatures in
past climates (Dansgaard et al., 1969; Masson-Delmotte et al., 2005; Jouzel,
2013), it is important to understand how the presented COSMO_iso simulations might be able to improve these isotope-based temperature
reconstructions. In general, the regional surface temperature variability
and the regional δ18O variability show similar patterns for
Greenland. In both cases the variability is high at the coast and low on the
inland plateau. Similar patterns to those in the mid-Holocene can also be seen for
the present-day simulations. These patterns of spatial variability in
δ18O and the surface temperature are in line with the results
of Sjolte et al. (2011) for RCM simulations under present-day conditions for
Greenland. Based on these patterns of variability, it can be derived that
the regional surface temperature variability highly depends on the surface
characteristics in Greenland. However, for the regional isotopic ratio
variability, this dependence appears to be less pronounced. At the coastline, a clear relationship between surface temperatures and measured δ18O ratios in ice cores can be deduced, while in central
Greenland this relation is weaker. These spatial differences might be
explained by the fact that isotope changes are an integrated signal of the
mesoscale variability of atmospheric processes (Dansgaard, 1964; Merlivat
and Jouzel, 1979; Gat, 1996), which might be partially decoupled from
surface temperature changes in homogeneous terrain.
Consistent patterns over Greenland are also modeled for the interannual
temporal δ18O–temperature slope in the mid-Holocene and the
present-day simulation. But in comparison to the spatial δ18O–temperature slope, the interannual temporal δ18O–temperature slope is rather small. This weaker interannual δ18O–temperature slope is again in line with the results of Sjolte et al. (2011).
The presented study demonstrates that the isotope-enabled MPI-ESM-wiso–COSMO_iso model chain with realistically implemented stable
water isotope fractionation processes constitutes a useful supplement to
reconstruct regional paleo-climate conditions during the mid-Holocene in
Greenland. By means of such an isotope-enabled GCM-RCM model chain, locally
measured isotope ratios in an ice core can be adequately linked to spatially
coarse climate model results, and conclusions on the underlying climatic
processes leading to these ratios can be drawn in a physically consistent
way. This approach might also be very helpful for other isotope-enabled GCMs
and understanding their deviations from observed isotope ratios in different
paleo-time periods and regions. Particularly in regions in which large
differences occur between simulated and observed δ18O ratios,
due to small-scale orographic variations, like parts of Europe and North
America (Cauquoin et al., 2019; Comas-Bru et al., 2019), an improved
representation of small-scale processes can potentially reduce these biases,
and consequently the reconstruction of regional paleo-climate patterns can
become more reliable. To test this hypothesis, in follow-up studies more
time slices will be simulated with the presented MPI-ESM-wiso–COSMO_iso model chain for different periods and different
regions.
Code availability
The isotope-enabled version COSMO_iso is available upon
request from Marcus Breil. The code of the isotopic version MPI-ESM-wiso is
available upon request on the AWI's GitLab repository (https://gitlab.awi.de/mwerner/mpi-esm-wiso, last access: 13 October 2019, Cauquoin et al., 2019, 10.5194/cp-15-1913-2019).
Data availability
The model data presented in the paper can
be accessed by contacting Marcus Breil. The data will be archived
at the DKRZ (German Climate Computing Centre).
The supplement related to this article is available online at: https://doi.org/10.5194/cp-17-1685-2021-supplement.
Author contributions
MB performed the COSMO_iso mid-Holocene simulations. EC performed the COSMO_iso present-day simulation. The ECHAM5-wiso simulations were performed by MW. The MPI-ESM-wiso simulations were performed by AC. MK performed the
comparison of the model results with present-day observations for the
standard climatological parameters. MB analyzed the presented δ18O model results and wrote the manuscript with contributions
from all co-authors.
Competing interests
The authors declare that they have no conflict of interest.
Disclaimer
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Special issue statement
This article is part of the special issue “Paleoclimate data synthesis and analysis of associated uncertainty (BG/CP/ESSD inter-journal SI)”. It is not associated with a conference.
Acknowledgements
All simulations were performed at the
German Climate Computing Center (DKRZ). Marcus Breil acknowledges support
from the KIT Publication Fund of the Karlsruhe Institute of Technology.
Financial support
This research has been supported by the German Federal Ministry of Education and Research (BMBF) (grant no. FKZ:01LP1511B).The article processing charges for this open-access publication were covered by the Karlsruhe Institute of Technology (KIT).
Review statement
This paper was edited by André Paul and reviewed by two anonymous referees.
ReferencesAemisegger, F., Spiegel, J. K., Pfahl, S., Sodemann, H., Eugster, W., and Wernli, H.: Isotope meteorology of cold front passages: A case study combining observations and modeling, Geophys. Res. Lett., 42, 5652–5660, 10.1002/2015GL063988, 2015.Ambach, W., Dansgaard, W., Eisner, H., and Moller, J.: The altitude effect
on the isotopic composition of precipitation and glacier ice in the Alps, Tellus, 20, 595–600, 10.1111/j.2153-3490.1968.tb00402.x, 1968.
Baldauf, M., Seifert, A., Förstner, J., Majewski, D., and Raschendorfer, M.: Operational convective-scale numerical weather prediction with the COSMO
model: description and sensitivities, Mon. Weather Rev., 139, 3887–3905, 2011.Bonne, J. L., Behrens, M., Meyer, H., Kipfstuhl, S., Rabe, B., Schönicke, L., Steen-Larsen, H. C., and Werner, M.: Resolving the controls of water vapour isotopes in the Atlantic sector, Nat. Commun., 10, 1632, 10.1038/s41467-019-09242-6, 2019.
Braconnot, P., Harrison, S. P., Kageyama, M., Bartlein, P. J., Masson-Delmotte, V., Abe-Ouchi, A., Otto-Bliesner, B., and Zhao, Y.: Evaluation of climate models using palaeoclimatic data, Nat. Clim. Change, 2, 417–424, 2012.Cauquoin, A., Werner, M., and Lohmann, G.: Water isotopes – climate relationships for the mid-Holocene and preindustrial period simulated with an isotope-enabled version of MPI-ESM, Clim. Past, 15, 1913–1937, 10.5194/cp-15-1913-2019, 2019 (data available at: https://gitlab.awi.de/mwerner/mpi-esm-wiso, last access: 19 October 2019).Christner, E., Kohler, M., and Schneider, M.: The influence of snow sublimation and meltwater evaporation on δD of water vapor in the atmospheric boundary layer of central Europe, Atmos. Chem. Phys., 17, 1207–1225, 10.5194/acp-17-1207-2017, 2017.Christner, E., Aemisegger, F., Pfahl, S., Werner, M., Cauquoin, A., Schneider, M., Hase, F., Barthlott, S., and Schädler, G.: The climatological impacts of continental surface evaporation, rainout, and subcloud processes on δD of water vapor and precipitation in Europe, J. Geophys. Res.-Atmos., 123, 4390–4409, 2018.Comas-Bru, L., Harrison, S. P., Werner, M., Rehfeld, K., Scroxton, N., Veiga-Pires, C., and SISAL working group members: Evaluating model outputs using integrated global speleothem records of climate change since the last glacial, Clim. Past, 15, 1557–1579, 10.5194/cp-15-1557-2019, 2019.Coppola, E., Sobolowski, S., Pichelli, E., et al.: A first-of-its-kind multi-model convection permitting ensemble for investigating convective phenomena over Europe and the Mediterranean, Clim. Dynam., 55, 3–34, 10.1007/s00382-018-4521-8, 2020.
Craig, H. and Gordon, L. I.: Deuterium and oxygen 18 variations
in the ocean and the marine atmosphere, in: Stable Isotopes in
Oceanographic Studies and Paleotemperatures, edited by: Ton-
giorgi, E., Spoleto, Italy, 9–130, 1965.
Dansgaard, W.: The abundance of O18 in atmospheric water and water vapour,
Tellus, 5, 461–469, 1953.
Dansgaard, W.: Stable isotopes in precipitation, Tellus, 16, 436–468, 1964.
Dansgaard, W., Johnsen, S. J., Møller, J., and Langway, C. C.: One
thousand centuries of climatic record from Camp Century on the Greenland ice
sheet, Science, 166, 377–380, 1969.
Dee, D. P., Uppala, S.M., Simmons, A.J., Berrisford, P., Poli, P., Kobayashi, S., Andrae, U., Balmaseda, M.A., Balsamo, G., Bauer, P., Bechtold, P., Beljaars, A.C.M., van de Berg, L., Bidlot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Geer, A.J., Haimberger, L., Healy, S.B., Hersbach, H., Hólm, E.V., Isaksen, L., Kållberg, P., Köhler, M., Matricardi, M., McNally, A.P., Monge-Sanz, B.M., Morcrette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thépaut, J.-N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteor. Soc., 137, 553–597, 2011.
Dütsch, M.: Stable water isotope fractionation processes
in weather systems and their influence on isotopic variability on different
time scales, Diss no. 23939, PhD thesis, ETH Zurich, 2017.
Felzer, B. and Thompson, S. L.: Evaluation of a regional climate model for
paleoclimate applications in the Arctic, J. Geophys. Res.-Atmos., 106, 27407–27424, 2001.Fischer, H.: Stable oxygen isotopes on snow pit ngt01C93 from the North
Greenland Traverse, PANGAEA [data set], 10.1594/PANGAEA.133399, 2003.
Fischer, H., Werner, M., Wagenbach, D., Schwager, M., Thorsteinnson, T., Wilhelms, F., Kipfstuhl, J., and Sommer, S.: Little ice age clearly recorded in northern Greenland ice cores, Geophys. Res. Lett., 25, 1749–1752, 1998.Fischer, N. and Jungclaus, J. H.: Evolution of the seasonal temperature cycle in a transient Holocene simulation: orbital forcing and sea-ice, Clim. Past, 7, 1139–1148, 10.5194/cp-7-1139-2011, 2011.
Furukawa, R., Uemura, R., Fujita, K., Sjolte, J., Yoshimura, K., Matoba, S., and Iizuka, Y.: Seasonal-Scale Dating of a Shallow Ice Core From Greenland
Using Oxygen Isotope Matching Between Data and Simulation, J. Geophys. Res.-Atmos., 122, 10873–10887, 2017.Gat, J. R.: Oxygen and hydrogen isotopes in the hydrological cycle, Annu. Rev. Earth Pl. Sc., 24, 225–262, 10.1146/annurev.earth.24.1.225, 1996.Grootes, P. M. and Stuiver, M.: Oxygen 18/16 variability in Greenland snow
and ice with 10-3- to 105-year time resolution, J. Geophys. Res., 102, 26455–26470, 1997.
Gurney, S. D. and Lawrence, D. S. L.: Seasonal trends in the stable
isotopic composition of snow and meltwater runoff in a subarctic catchment
at Okstindan, Norway, Hydrol. Res., 35, 119–137, 2004.
Harrison, S. P., Bartlein, P. J., Brewer, S., Prentice, I. C., Boyd, M., Hessler, I., Holmgren, K., Izumi, K., and Willis, K.: Climate model benchmarking with glacial and mid-Holocene climates, Clim. Dynam., 43, 671–688, 2014.IAEA: International Atomic Energy Agency/World Meteorological Organization, Global Network of Isotopes in Precipitation, The GNIP Database, available at: https://nucleus.iaea.org/wiser (last access: 18 December 2019), 2018.
Johnsen, S. J., Clausen, H. B., Cuffey, K. M., Hoffmann, G., Schwander, J., and Creyts, T.: Diffusion of stable isotopes in polar firn and ice: the
isotope effect in firn diffusion, in: Physics of ice core records, Hokkaido University Press, Sapporo, Japan, 121–140, 2000.
Jost, A., Lunt, D., Kageyama, M., Abe-Ouchi, A., Peyron, O., Valdes, P. J.,
and Ramstein, G.: High-resolution simulations of the last glacial maximum
climate over Europe: a solution to discrepancies with continental
palaeoclimatic reconstructions?, Clim. Dynam., 4, 577–590, 2005.Jouzel, J.: A brief history of ice core science over the last 50 yr, Clim. Past, 9, 2525–2547, 10.5194/cp-9-2525-2013, 2013.
Jouzel, J. and Merlivat, L.: Deuterium and oxygen 18 in precipitation:
Modeling of the isotopic effects during snow formation, J. Geophys. Res.-Atmos., 89, 11749–11757, 1984.Jungclaus, J. H., Fischer, N., Haak, H., Lohmann, K., Marotzke,J., Matei,
D., Mikolajewicz, U., Notz, D., and von Storch, J. S.: Characteristics of
the ocean simulations in the Max Planck Institute Ocean Model (MPIOM) the
ocean component of the MPI-Earth system model, J. Adv. Model. Earth Sy., 5,
422–446, 10.1002/jame.20023, 2013.Kageyama, M., Braconnot, P., Harrison, S. P., Haywood, A. M., Jungclaus, J. H., Otto-Bliesner, B. L., Peterschmitt, J.-Y., Abe-Ouchi, A., Albani, S., Bartlein, P. J., Brierley, C., Crucifix, M., Dolan, A., Fernandez-Donado, L., Fischer, H., Hopcroft, P. O., Ivanovic, R. F., Lambert, F., Lunt, D. J., Mahowald, N. M., Peltier, W. R., Phipps, S. J., Roche, D. M., Schmidt, G. A., Tarasov, L., Valdes, P. J., Zhang, Q., and Zhou, T.: The PMIP4 contribution to CMIP6 – Part 1: Overview and over-arching analysis plan, Geosci. Model Dev., 11, 1033–1057, 10.5194/gmd-11-1033-2018, 2018.
Karremann, M. K. and Schädler, G.: Parametrisation of variables having
an impact on the Greenland surface mass balance using regional climate model
simulations, Atmosphere, in preparation, 2021.
Marcott, S. A., Shakun, J. D., Clark, P. U., and Mix, A. C.: A
reconstruction of regional and global temperature for the past 11,300 years, Science, 339, 1198–1201, 2013.
Masson-Delmotte, V., Jouzel, J., Landais, A., Stievenard M., Johnsen, S. J., White, J. W. C., Werner, M., Sveinbjornsdottir A., and Fuhrer, K.: GRIP deuterium excess reveals rapid and orbital-scale changes in Greenland moisture origin, Science, 309, 118–121, 2005.Masson-Delmotte, V., Steen-Larsen, H. C., Ortega, P., Swingedouw, D., Popp, T., Vinther, B. M., Oerter, H., Sveinbjornsdottir, A. E., Gudlaugsdottir, H., Box, J. E., Falourd, S., Fettweis, X., Gallée, H., Garnier, E., Gkinis, V., Jouzel, J., Landais, A., Minster, B., Paradis, N., Orsi, A., Risi, C., Werner, M., and White, J. W. C.: Recent changes in north-west Greenland climate documented by NEEM shallow ice core data and simulations, and implications for past-temperature reconstructions, The Cryosphere, 9, 1481–1504, 10.5194/tc-9-1481-2015, 2015.
Merlivat, L. and Jouzel, J.: Global climatic interpretation of the
deuterium-oxygen 18 relationship for precipitation, J. Geophys. Res., 84, 5029–5033, 1979.
Moossen, H., Bendle, J., Seki, O., Quillmann, U., and Kawamura, K.: North
Atlantic Holocene climate evolution recorded by high-resolution terrestrial
and marine biomarker records, Quaternary Sci. Rev., 129, 111–127, 2015.Pfahl, S., Wernli, H., and Yoshimura, K.: The isotopic composition of precipitation from a winter storm – a case study with the limited-area model COSMOiso, Atmos. Chem. Phys., 12, 1629–1648, 10.5194/acp-12-1629-2012, 2012.Risi, C., Bony, S., Vimeux, F., and Jouzel, J.: Water-Stable Isotopes in the LMDZ4 General Circulation Model: Model Evaluation for Present-Day and Past Climates and Applications to Climatic Interpretations of Tropical Isotopic Records, J. Geophys. Res.-Atmos., 115, D12118, 10.1029/2009JD013255, 2010.
Rockel, B., Will, A., and Hense, A.: The regional climate model COSMO-CLM
(CCLM), Meteorol. Z., 17, 347–348, 2008.
Schrodin, R. and Heise, E.: The multi-layer version of the DWD soil model
TERRA-LM, Consortium for Small-Scale Modelling (COSMO), Tech. Rep., 2, 16 pp., 2001.Shi, X., Risi, C., Pu, T., Lacour, J.-L., Kong, Y., Wang, K., He, Y., and Xia, D.: Variability of isotope composition of precipitation in the Southeastern Tibetan Plateau from the synoptic to seasonal time scale, J. Geophys. Res.-Atmos., 125, e2019JD031751, 10.1029/2019JD031751, 2020.Sjolte, J., Hoffmann, G., Johnsen, S. J., Vinther, B. M., Masson-Delmotte,
V., and Sturm, C.: Modeling the water isotopes in Greenland precipitation
1959–2001 with the meso-scale model REMO-iso, J. Geophys. Res., 116, D18105, 10.1029/2010JD015287, 2011.
Smith, R. S. and Gregory, J.: The last glacial cycle: transient simulations
with an AOGCM, Clim. Dynam., 38, 1545–1559, 2012.Steen-Larsen, H. C., Masson-Delmotte, V., Hirabayashi, M., Winkler, R., Satow, K., Prié, F., Bayou, N., Brun, E., Cuffey, K. M., Dahl-Jensen, D., Dumont, M., Guillevic, M., Kipfstuhl, S., Landais, A., Popp, T., Risi, C., Steffen, K., Stenni, B., and Sveinbjörnsdottír, A. E.: What controls the isotopic composition of Greenland surface snow?, Clim. Past, 10, 377–392, 10.5194/cp-10-377-2014, 2014.
Steig, E. J., Mayewski, P. A., Dixon, D. A., Kaspari, S. D., Frey, M. M., Schneider, D. P., Arcone, S. A., Hamilton, G. S., Blue Spikes, V., Albert, M., Meese, D., Gow, A. J., Shuman, C. A., White, J. W. C., Sneed, S., Flaherty, J., and Wumkes, M.: High-resolution ice cores from US ITASE (West
Antarctica): development and validation of chronologies and determination of
precision and accuracy, Ann. Glaciol., 41, 77–84, 2005.
Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric component of the MPI-M earth system
model: ECHAM6, J. Adv. Model. Earth Sy., 5, 146–172, 2013.Sturm, K., Hoffmann, G., Langmann, B., and Stichler, W.: Simulation of
δ18O in precipitation by the regional circulation model REMOiso, Hydrol. Process., 19, 3425–3444, 2005.
Sturm, C., Hoffmann, G., and Langmann, B.: Simulation of the stable water
isotopes in precipitation over South America: Comparing regional to global
circulation models, J. Climate, 20, 3730–3750, 2007.
Timm, O. and Timmermann, A.: Simulation of the last 21 000 years using
accelerated transient boundary conditions, J. Climate, 20, 4377–4401, 2007.
Torma, C., Giorgi, F., and Coppola, E.: Added value of regional climate
modeling over areas characterized by complex terrain – Precipitation over
the Alps, J. Geophys. Res.-Atmos., 120, 3957–3972, 2015.Vinther, B. M., Clausen, H. B., Johnsen, S. J., Rasmussen, S. O., Andersen, K. K., Buchardt, S. L., Dahl-Jensen, D., Seierstad, I. K., Siggaard-Andersen, M.-L., Steffensen, J. P., Svensson, A., Olsen, J., and Heinemeier, J.: A synchronized dating of three Greenland ice cores throughout the Holocene, J. Geophys. Res., 111, D13102, 10.1029/2005JD006921, 2006.Vinther, B. M., Clausen, H. B., Fisher, D. A., Koerner, R. M., Johnsen, S. J., Andersen, K. K., Dahl-Jensen, D., Rasmussen, S. O., Steffensen, J. P., and Svensson, A. M.: Synchronizing ice cores from the Renland and
Agassiz ice caps to the Greenland Ice Core Chronology, J. Geophys. Res., 113, D08115, 10.1029/2007JD009143, 2008.von Storch, H., Langenberg, H., and Feser, F.: A Spectral Nudging Technique
for Dynamical Downscaling Purposes, Mon. Weather Rev., 128, 3664–3673,
10.1175/1520-0493(2000)128<3664:ASNTFD>2.0.CO;2, 2000.Weißbach, S., Wegner, A., Opel, T., Oerter, H., Vinther, B. M., and
Kipfstuhl, S.: Accumulation rate and stable oxygen isotope ratios of the ice
cores from the North Greenland Traverse, PANGAEA [data set], 10.1594/PANGAEA.849161, 2016a.Weißbach, S., Wegner, A., Opel, T., Oerter, H., Vinther, B. M., and Kipfstuhl, S.: Spatial and temporal oxygen isotope variability in northern Greenland – implications for a new climate record over the past millennium, Clim. Past, 12, 171–188, 10.5194/cp-12-171-2016, 2016b.Werner, M., Langebroek, P. M., Carlsen, T., Herold, M., and Lohmann, G.:
Stable water isotopes in the ECHAM5 general circulation model: Toward
high-resolution isotope modeling on a global scale, J. Geophys. Res., 116, D15109, 10.1029/2011JD015681, 2011.Yoshimori, M. and Suzuki, M.: The relevance of mid-Holocene Arctic warming to the future, Clim. Past, 15, 1375–1394, 10.5194/cp-15-1375-2019, 2019.