Articles | Volume 14, issue 4
Research article
08 Apr 2020
Research article |  | 08 Apr 2020

How useful is snow accumulation in reconstructing surface air temperature in Antarctica? A study combining ice core records and climate models

Quentin Dalaiden, Hugues Goosse, François Klein, Jan T. M. Lenaerts, Max Holloway, Louise Sime, and Elizabeth R. Thomas

Improving our knowledge of the temporal and spatial variability of the Antarctic Ice Sheet (AIS) surface mass balance (SMB) is crucial to reduce the uncertainties of past, present, and future Antarctic contributions to sea level rise. An examination of the surface air temperature–SMB relationship in model simulations demonstrates a strong link between the two. Reconstructions based on ice cores display a weaker relationship, indicating a model–data discrepancy that may be due to model biases or to the non-climatic noise present in the records. We find that, on the regional scale, the modeled relationship between surface air temperature and SMB is often stronger than between temperature and δ18O. This suggests that SMB data can be used to reconstruct past surface air temperature. Using this finding, we assimilate isotope-enabled SMB and δ18O model output with ice core observations to generate a new surface air temperature reconstruction. Although an independent evaluation of the skill is difficult because of the short observational time series, this new reconstruction outperforms the previous reconstructions for the continental-mean temperature that were based on δ18O alone. The improvement is most significant for the East Antarctic region, where the uncertainties are particularly large. Finally, using the same data assimilation method as for the surface air temperature reconstruction, we provide a spatial SMB reconstruction for the AIS over the last 2 centuries, showing large variability in SMB trends at a regional scale, with an increase (0.82 Gt yr−2) in West Antarctica over 1957–2000 and a decrease in East Antarctica during the same period (−0.13 Gt yr−2). As expected, this is consistent with the recent reconstruction used as a constraint in the data assimilation.

1 Introduction

The spatial coverage of climate observations in Antarctica and the Southern Ocean is sparse (e.g., Jones et al.2016; Neukom et al.2018). Consequently, the climate dynamics of the high southern latitudes are still poorly understood, leading to large uncertainties in the processes governing climate variability (Church et al.2013). Since around 1995, the contribution to the global sea level rise from the ice sheets – Greenland Ice Sheet (GrIS) and the Antarctic Ice Sheet (AIS) – has strongly increased and is slowly outpacing the contributions from mountain glaciers and ocean thermal expansion (Shepherd et al.2018). The GrIS has been dominating the ice sheet contribution so far (Rignot et al.2019), but AIS mass loss increased 5-fold in 2012–2017 relative to 1992–1997, with current AIS mass loss values that approach those of the GrIS.

The (grounded) AIS mass balance (MB) is the difference between the surface mass balance (SMB) and the solid ice discharge (Lenaerts et al.2019; Fyke et al.2018). Reliable estimates of AIS MB and its relationship with internal climate variability and transient climate forcing are needed to constrain future climate and sea level projections (Bamber et al.2018). The current AIS MB is negative (Rignot et al.2019) because of large values of ice discharge (Shepherd et al.2018).

The SMB is defined as the difference between the incoming and outgoing mass at the surface of the ice sheet. In Antarctica, the main source term of the SMB, and its interannual variations, is precipitation in the form of snow (e.g., Lenaerts et al.2012; Agosta et al.2018). Unlike Greenland, AIS surface melt is small, and most surface melt water refreezes in place, not contributing to SMB (Trusel et al.2015; Kuipers Munneke et al.2012). As a result, the surface sublimation and sublimation of blowing snow are the main sink terms of the AIS SMB (e.g., Frezzotti et al.2013; van Wessem et al.2018).

Ice cores provide information on past changes in surface air temperature (SAT) and SMB across Antarctica on timescales of centuries to millennia (e.g., Stenni et al.2017a; Thomas et al.2017). In particular, it has become standard to reconstruct past temperature changes from water stable isotopes, and in particular δ18O (e.g., Jouzel2003; Masson-Delmotte et al.2006). However, ice core studies suffer from several limitations: (1) the ice core network is still relatively sparse, despite recent coordinated international drilling efforts (Thomas et al.2017; Stenni et al.2017a); (2) annually resolved surface air temperature and SMB records are not available from extremely dry areas, such as the East Antarctic Plateau; (3) changes in precipitation seasonality (e.g., Sime et al.2008), moisture origin (e.g., Holloway et al.2016a), and other processes can modify the expected relationship between δ18O and surface air temperature (e.g., Jouzel et al.1997; Klein et al.2019). Combined, these factors lead to large uncertainties in the reconstruction of surface air temperatures.

Until recently, AIS SMB had been considered to display no significant trends since the mid-twentieth century (Monaghan et al.2006; Frezzotti et al.2013). Based on recent work, this hypothesis has been revised: using a larger ice core network (PAGES2k database), Thomas et al. (2017) and Medley and Thomas (2019) have shown that AIS SMB has increased significantly since 1900, albeit with important regional differences. The Antarctic Peninsula has witnessed a considerable SMB increase during the twentieth century (e.g., Thomas et al.2015; Goodwin et al.2016), as well as some regions of Dronning Maud Land (e.g., Philippe et al.2016; Lenaerts et al.2013; Medley et al.2018; Shepherd et al.2012). In contrast, other regions of Dronning Maud Land have been subject to a SMB decrease over the recent past (Schlosser et al.2014; Altnau et al.2015). All these studies point out the need to densify the ice core network over Antarctica, but also to gain more insight into what is driving the trends in AIS SMB and its spatial signatures. For the latter, output of climate model simulations can be very useful (e.g., Lenaerts et al.2018).

In the last decade, output of several climate model simulations that cover the last millennium has become available (Schmidt et al.2011). Thus far, model evaluation has been mainly focused on surface air temperature (PAGES 2k-PMIP3 group2015). These results have shown discrepancies in AIS surface air temperature between climate model simulations and reconstructions during the last millennium. In contrast to climate model results, surface air temperature reconstructions show no clear warming over the 20th century at the continental scale (Goosse et al.2012; Stenni et al.2017a; PAGES 2k-PMIP3 group2015; Neukom et al.2018). This mismatch can be explained by an overestimation of the response of climate models to external forcing, the strong natural variability occurring in Antarctica, underestimation of the signal in proxy-based reconstructions, or a combination of all of these (Jones et al.2016; Neukom et al.2018). Unlike temperature changes, modeled AIS SMB variations over the past millennium are poorly documented.

In a warmer climate, AIS SMB is expected to increase due to higher snowfall associated with the greater moisture holding capacity at higher air temperature (e.g., Lenaerts et al.2016). Taken alone, this straightforward thermodynamical effect would mitigate the sea level rise (Huybrechts et al.2004; Krinner et al.2007; Frieler et al.2015). According to Frieler et al. (2015), the observed sensitivity of Antarctic snowfall accumulation to surface air temperature was about 5 % K−1 during the 1960–1999 period. Based on climate model simulations, this sensitivity is expected to increase in the future with an estimated conversion value of 7.4 % K−1 for the end of the 21st century (2080–2099; Palerme et al.2017). The link between surface air temperature and SMB has been confirmed for small regions at the centennial timescale (200 years; e.g., Oerter et al.2000; Medley et al.2018) and on longer timescales (glacial–interglacial; Frieler et al.2015) for the full AIS using climate models and ice cores. However, some studies using surface air temperature reconstructions based on δ18O data (Fudge et al.2016; Altnau et al.2015; Philippe et al.2016; Goursaud et al.2019) suggest that this SMB–surface air temperature relationship is not always positive and varies spatially and temporally. These results suggest that in some regions, especially along the AIS coasts, the SMB variability is dominated by large-scale atmospheric circulation rather than by thermodynamic processes (such as the Clausius–Clapeyron relation), limiting the correlation with δ18O.

The first goal of this study is to document the relationship between surface air temperature and SMB in Antarctica on a regional scale using climate models and ice core records over the past 2 centuries and over the last millennium. The final goal is to use the covariance between both variables to reconstruct past changes over the last 2 centuries by using a data assimilation procedure. All reconstruction methods depend on the number and quality of the input data. However, while the statistical methods classically used to infer past surface air temperature (see for instance Stenni et al.2017a) rely on the length of the calibration period, on the quality of the record during this period, and on the stationarity of the link between the proxy and the variable of interest, which can be strong assumptions in the case of the δ18O–temperature relationship (Klein et al.2019), data assimilation does not. In recent years, data assimilation has become a standard procedure in paleoclimatology to optimally combine the information from model results and proxies and to provide estimates of past climate states (e.g., Hakim et al.2016; Widmann et al.2010; Goosse et al.2010; Matsikaris et al.2015; Steiger et al.2014). Nevertheless, Antarctic SMB, to the best of our knowledge, has never been assimilated in a climate model. The biggest advantage of using data assimilation is that it takes into account information brought by both SMB and δ18O without making the strong assumptions that the statistical methods do. Additionally, using the covariance between them might lead to better estimates of past changes in the two variables, particularly when proxy records are scarce and few instrumental data are available, which is the case for Antarctica. The resulting reconstructions will have the benefit of being compatible with the physics of the climate system as represented by the models.

2 Data: model simulations and observations

2.1 Global climate model simulations

The climate model simulations selected for this study are those for which the required variables (i.e., precipitation and sublimation–evaporation) are available for the last millennium from the PMIP3–CMIP5 database (Otto-Bliesner et al.2009; Taylor et al.2012). In addition to these simulations, the CESM1–CAM5 model simulations covering the last millennium (Otto-Bliesner et al.2015) are also used. The characteristics and references of each model are described in Table A1. All these general circulation models (GCMs) use the GMTED2010 elevation dataset (Danielson and Gesch2011) as topography, adapted to their spatial horizontal resolution. The simulations are driven by both natural (orbital, solar and volcanic) and anthropogenic (greenhouse gases, land use, aerosol, and ozone) forcings through the last millennium (Schmidt et al.2011, 2012). Except for CESM1–CAM5, CSIRO Mk3L-1-2, and MPI-ESM-P, the simulations do not cover the entire millennium. Historical simulations covering 1851–2005 CE were launched independently of simulations covering 850–1850 CE (referred to as the past1000 experiment). In order to obtain results over the full millennium, we adopt the approach from Klein and Goosse (2018) and merge the first ensemble members (r1i1p1) of the past1000 experiment with the corresponding ensemble members of the historical experiment. Although not continuous, there is no large discrepancy between the two merged simulations (e.g., Klein and Goosse2018).

Simulations performed with the isotope-enabled climate models, ECHAM5–MPI-OM (Sjolte et al.2018), ECHAM5-wiso (Steiger et al.2017), and iHadCM3 (Tindall et al.2009; Holloway et al.2016b) are also analyzed. These simulations allow for a direct comparison with observed water isotope content. ECHAM5–MPI-OM is a fully coupled general circulation model (GCM). The simulation used here covers the period 800–2000 CE forced by natural and anthropogenic forcing (Sjolte et al.2018). The horizontal resolution of the atmospheric model is 3.75×3.75. The simulation of ECHAM5-wiso, which only includes an atmospheric component, was performed by Steiger et al. (2017) and covers the period 1871–2011 CE at  1 resolution. The model is driven by the sea surface temperature and sea ice from the Rayner et al. (2003) dataset. Due to a lack of Antarctic sea ice data before 1973, this dataset is based on historical climatologies of sea ice concentration for the period 1871–1973 CE, with no interannual variability. Finally, iHadCM3 is the version of HadCM3 (fully coupled climate model; Turner et al.2016) which has an explicit representation of the water isotopes. The resolution of the atmospheric model is 3.75×2.5. While only one simulation is available for ECHAM5–MPI-OM and ECHAM5-wiso, we have an ensemble of seven iHadCM3 simulations spanning the industrial period from 1851 to 2003 CE. The initial conditions for each of these simulations correspond to different years in the pre-industrial control simulation of the iHadCM3 model. Comparisons of the results of these three isotope-enabled models with modern δ18O observations indicate that they all reproduce the main characteristics of the spatial distribution of the isotopic composition of precipitation over Antarctica including the latitudinal distribution (negative δ18O gradient from the coast to the plateau). According to Tindall et al. (2009) and Sime et al. (2008), the small biases in δ18O (for example, an underestimation of the spatial δ18O variability in rugged areas) in the iHadCM3 simulation mainly come from the coarse horizontal resolution of the model and not from the isotopic model itself. ECHAM5-wiso and ECHAM5–MPI-OM display an overall underestimation of δ18O in Antarctica but reproduce the general Antarctic δ18O pattern well (Goursaud et al.2018; Klein et al.2019, see reference of each model for more details).

Klein et al. (2019) recently described an evaluation of Antarctic surface air temperature in reconstructions and model simulations over the last millennium. In accordance with Abram et al. (2016), they highlighted the early onset of industrial warming simulated by the PMIP and CMIP models, which is not observed in the δ18O-based temperature reconstructions of Stenni et al. (2017a). This suggests that the Antarctic surface air temperatures simulated by the models are too sensitive to the anthropogenic forcing.

2.2 The regional climate model RACMO2 simulation

The evaluation of AIS SMB simulated by GCMs for the present period (1979–2005) is mainly based on the results of the regional atmospheric climate model RACMO2.3p_2 (RACMO2 hereafter) covering the entire AIS over 1979–2016 (van Wessem et al.2018). This is because (1) the SMB observations are very sparse on the AIS (Favier et al.2013); (2) the interannual (year-to-year) variability is different between observations and GCMs given that the latter are freely evolving coupled models. Consequently, the comparison can be only made on multi-decadal timescales (> 20 years), which drastically reduces the availability of observations; (3) unlike observations, RACMO2 provides a complete SMB field over the entire AIS; and, finally, (4) RACMO2 has been extensively evaluated against available measurements and displays a very good agreement (e.g., van Wessem et al.2018; Lenaerts et al.2012). In an intercomparison of AIS SMB from reanalysis, atmospheric models, and observations, Wang et al. (2016) showed that the RACMO2 model best fits the recent AIS SMB observations compared to all other available datasets.

RACMO2 combines the physics package of the European Centre for Medium-Range Weather Forecasts (ECMWF2008) Integrated Forecast System and the hydrostatic dynamics of the High Resolution Limited Area Model (HIRLAM, Unden et al.2002). RACMO2 is specially adapted to polar regions since it includes the interactions between the atmosphere and the multilayered snow model that calculates physical processes occurring in the firn: meltwater production, percolation, runoff, refreezing, and snow grain size and resulting snow albedo (Greuell and Thomas1994; Ettema et al.2010). RACMO2 also includes a drifting snow scheme simulating the interactions between the near-surface air and snow (Lenaerts et al.2010). All the SMB components are explicitly calculated by this regional model on a 27 km resolution grid. The digital elevation model of Bamber et al. (2009) is taken as the reference of the Antarctic topography. ERA-Interim reanalysis data (Dee et al.2011) are used to force the regional model at its lateral boundaries. For more details on RACMO2, see van Wessem et al. (2018).

2.3 Snow accumulation database from Antarctica2k

The annually resolved Antarctica2k (Ant2k) snow accumulation database (Thomas et al.2017) is used for the evaluation of AIS SMB simulated by GCMs before 1979. The estimate of the SMB from ice cores is based on the physical distance between suitable age markers within the ice core. The age markers used depend on the timescale of interest, ranging from glacial cycles (e.g., bulk changes in isotopic compositions) to seasonal variations reflected by changes in stable water isotopes, while volcanic eruptions can provide information on decadal to millennial timescales (Dansgaard and Johnsen1969). Once the age markers are identified, since the firn density generally increases with depth in the ice core, it is necessary to consider those variations to convert the age and depth to mass (Van Den Broeke et al.2008). Doing so, SMB is converted to meters of water equivalent based on measured density and corrected for the vertical strain rate effect – the differential vertical velocity with depth leading to layer thinning with depth (Thomas et al.2017).

Figure 1Antarctic regions used in this study. The definitions of the regions are those of Thomas et al. (2017).

This database is composed of 79 records that are assigned to seven geographical regions (Fig. 1) with distinctly different climates. East Antarctica above 2000 m elevation constitutes the East Antarctica Plateau (EAP). West Antarctica is separated into two parts: the Antarctic Peninsula (AP) and the West Antarctic Ice Sheet (WAIS), with a division at 88 W. The coastal region of East Antarctica is divided into four regions: Victoria Land (VL; 150–170 E), the Wilkes Land coast (WL; 70–150 E), Dronning Maud Land (DML; 15 W–150 E), and the Weddell Sea coast (WS; 15–60 E). For each region, this database covers the past 1000 years except for EAP, AP, and DML, for which the period covered is 1240–2005 CE, 1703–2010 CE, and 1737–2010 CE, respectively. Hereafter, West Antarctica is composed of WAIS and AP, while East Antarctica comprises all of the other regions. Since some Antarctic regions lack long-term data, the SMB reconstruction for the whole Antarctic ice sheet is only available from 1737 CE onward. This regional SMB reconstruction has been compared to RACMO2, concluding that the reconstruction captures a large proportion of the regional spatial SMB variability as defined by RACMO2 for the 1979–2010 period (Thomas et al.2017).

2.4 Water stable isotope records and surface air temperature reconstructions from Antarctica2k

Stenni et al. (2017a) built δ18O regional composites from 112 individual ice cores compiled in the framework of the PAGES Antarctica2k working group for seven similar Antarctic regions as in Thomas et al. (2017, see Sect. 2.3) over the last 2 millennia. Based on those δ18O composites, they reconstructed regional surface air temperatures over the last 2 millennia based on the statistical relationship between δ18O and surface air temperature. Three methods have been used to scale the δ18O composites. The second reconstruction (borehole reconstruction) is used throughout this study for two reasons: (1) this is not based on surface air temperature observations, which are used here to estimate the skill of the reconstructions which would have led to a bias; (2) because it is based on more information, the borehole reconstruction is expected to be more accurate (see Sect. S1 of the Supplement for details). The temporal resolution is the same as for the δ18O composites: 10 years over 0–1800 and 5 years over 1800–2010.

3 Methods: reconstructing SMB and surface air temperatures using data assimilation

3.1 Data assimilation method: a particle filter using fixed ensembles

Data assimilation optimally combines observations (proxy data in our case) and climate model states. Two types of data assimilation methods are usually applied in paleoclimatology. First, online methods follow standard sequential data assimilation approaches, in which the analysis at a single time step depends on the state at the previous step. Information is thus propagated forward in time. However, because data assimilation requires a large ensemble of model simulations (tens to hundreds), for paleoclimate reconstructions, performing online data assimilation at high spatial climate model resolution (e.g., CMIP5 class as used here) becomes impractical. Second, when working with so-called offline methods, ensemble members are constructed from existing model simulations, which is of great interest in terms of computation time compared to online methods. Here, ensemble members are constructed by individual years and not by independent model simulations. Therefore, in contrast to online methods, offline methods do not maintain temporal consistency. However, when the predictability on interannual timescales is limited, for surface air temperature or precipitation because of the dominant role of their chaotic nature for example, online methods do not outperform offline ones (Matsikaris et al.2015). Indeed, offline methods have provided skillful data-assimilation-based reconstructions for various types of data (e.g., Steiger et al.2017; Klein and Goosse2018; Hakim et al.2016). Nevertheless, the online approach is preferred when focusing on ocean dynamics because of the long ocean memory (e.g., Goosse2017; Pendergrass et al.2012).

The offline data assimilation method applied in this study is based on a particle filter (e.g., van Leeuwen2009; Dubinkina et al.2011) using fixed ensembles from climate model outputs. The implementation described in Dubinkina et al. (2011) is identical to that of previous studies (e.g., Klein and Goosse2018). Hence, only a brief description of the methodology will be given here. At each time step of the data assimilation procedure (yearly; see Sect. 3.2), each ensemble member, called particle, is compared to the proxy-based reconstruction by computing its likelihood, assumed here to be Gaussian, taking into account data uncertainties (see Dubinkina et al.2011, for details). Depending on its likelihood, each particle receives a weight. Then, all the weights are multiplied by the number of particles and rounded to the nearest integer toward negative infinity by ensuring that the sum of the weights equals the number of particles throughout the data assimilation process (see Dubinkina et al.2011, for details). Considering all particle weights, we can compute a weighted average, providing a reconstruction for this time step. In this study, the ensemble members are derived from the three isotope-enabled climate model outputs ECHAM5–MPI-OM, ECHAM5-wiso, and iHadCM3, presented in Sect. 2.1. These models have been chosen because they explicitly simulate δ18O. Because iHadCM3 offers an ensemble of seven simulations, while the other isotope-enabled models have only a single realization, we mainly focus on the iHadCM3 outputs in the main text. The results from the other models are shown in the Supplement and provide similar conclusions to the ones obtained with iHadCM3. Using an ensemble allows us to estimate the contribution of internal variability over the last century, and the range provided by this ensemble is larger than the one given by the other two models, increasing the probability of finding model results close to the assimilated records during the data assimilation process. Additionally, the iHadCM3 ensemble provides a more realistic Antarctic surface air temperature and snow accumulation over the recent past than the ECHAM5–MPI-OM simulation (Fig. S8; Klein et al.2019).

3.2 Experiment design

Data assimilation is used in this study to reconstruct surface air temperature and SMB by taking advantage of the covariance between these variables. They are assimilated together as well as separately in three different experiments. In the first experiment, the seven region composites of δ18O data (Stenni et al.2017a) are used to constrain model results. Assimilating δ18O instead of surface air temperature potentially accounts for the nonstationarity and the nonlinearity of the stable oxygen ratios–surface air temperature link (Masson-Delmotte et al.2008; Klein et al.2019). For the second experiment, the SMB reconstruction for the seven regions (Thomas et al.2017) is used in the data assimilation process. Finally, both δ18O and SMB are taken into account together in the last experiment. This allows us to estimate independently the consistency of the reconstructed SMB and surface air temperature between the various records and model results. In addition, our experiments allow us to assess the information acquired on surface air temperature by assimilating SMB and the information on SMB by assimilating δ18O. In all the experiments, we assimilate annual-mean proxies. All modeled δ18O values are precipitation-weighted as this quantity is most similar to the one measured in ice cores.

Since the number of ice cores is limited before 1800 CE (both for δ18O and for SMB), which drastically decreases the quality of the regional composites (Thomas et al.2017), the experiments are performed on the 1800–2010 period. Contrary to the SMB composites, which have an annual resolution, the composites of δ18O are 5-year averages. Consequently, the δ18O data have been interpolated linearly over the studied period to match the temporal resolution of the SMB reconstruction. However, as recommended by Stenni et al. (2017a), the results are analyzed only for the 5-year averages. This temporal averaging reduces uncertainties in dating and the noise induced by non-climatic processes (e.g., Laepple et al.2018; Fan et al.2014).

In order to assess the skill of our data-assimilation-based surface air temperature reconstructions, we evaluate them at first with the reconstructions of Stenni et al. (2017a). But this is biased since they are only based on δ18O and we cannot thus evaluate the added value acquired from SMB data and model physics in the data assimilation experiments. Therefore, independent data are needed to properly assess the potential of SMB and δ18O in reconstructing surface air temperature. This is done here using the surface air temperature reconstruction from Nicolas and Bromwich (2014a), which is based on surface air temperature records and not on δ18O data, over the 1958–2010 period. SMB estimates are also available for the last decades (e.g., Medley et al.2014), but they cover a too short period or have a too small spatial coverage to provide an independent validation of our reconstruction. It is thus not possible to estimate if the assimilation of SMB and δ18O measurements provides an improvement on SMB estimates.

Figure 2Surface mass balance anomalies (Gt yr−1) simulated by the GCMs (Table A1) and snow accumulation reconstructions (Thomas et al.2017) from 1000 to 2005 and from 1800 to 2005 for West Antarctica, East Antarctica, and Antarctica as a whole. Anomalies are relative to the 1871–2000 period. The shaded area corresponds to the range of the CESM1–CAM5 simulations. For visibility, data have been smoothed with a 100-year moving average for the last millennium and a 30-year moving average for the last 200 years. The equivalent for the seven regions is given in Fig. S1 in the Supplement.


4 Results

4.1 Reconstructed and simulated SMB changes over the last centuries

The AIS SMB over the last millennium has been estimated for each GCM by computing the difference between precipitation and sublimation–evaporation. Runoff is assumed to be negligible as surface meltwater generally refreezes in the cold firn (Magand et al.2008; Kuipers Munneke et al.2012). Our short evaluation of SMB simulated by GCMs over the present day (see Sect. S3) suggests that the selected GCMs (including the isotope-enabled models) display reasonable SMB climatology when compared to RACMO outputs.

Before the 19th century, all GCM simulations are characterized by large decadal variability, but no long-term trend (Fig. 2). A positive trend, albeit initiated at different times, is shown at the end of the simulation (around 1950 CE). All models agree on an AIS SMB increase from  1975 onwards, which is consistent with the SMB reconstruction of Thomas et al. (2017). However, the contrast in the SMB trends between East Antarctica and West Antarctica is clearly stronger in the reconstruction based on ice cores than in GCMs on average. Indeed, over the last decades (1950–2000), the ice core SMB reconstruction shows a large increase for West Antarctica (25.6 Gt yr−1 per decade) and a small decrease (−3.6 Gt yr−1 per decade) for East Antarctica, while, on average, the GCMs simulate a strong SMB increase in both regions (8.9±9.2 Gt yr−1 per decade and 14.2±13.5 Gt yr−1 per decade, respectively; Figs. 2 and 3 and Table S1 in the Supplement). When analyzing the ensemble of simulations performed with CESM1–CAM5, the ensemble mean also shows a relatively homogeneous increase, but some simulations display a contrast between East Antarctica and West Antarctica close to the one observed in the reconstruction (Fig. 3). This suggests that internal variability has a dominant contribution in the current Antarctic SMB changes and might explain why the observed contrast between East and West Antarctica is only present in a few simulations.

Figure 3(a) Comparison between the reconstructed and the simulated SMB trends (mm w.e. 100 yr−2) over the period 1950–2000 CE in West Antarctica (y axis) and East Antarctica (x axis). (b) As on the left but for SMB–SAT sensitivity factors (% K−1). For the reconstruction, data from Thomas et al. (2017) and Nicolas and Bromwich (2014a) are used.


4.2 Relationship between SMB and surface air temperatures in Antarctica

Averaged across all GCMs, the relationship between SMB and surface air temperature is positive for each Antarctic region (Fig. 4). A very similar result is obtained when the annual-mean surface air temperature and SMB derived from the RACMO2 simulation over the recent period (1979–2016) are used. The regional correlations are much weaker for the reconstructions based on ice cores than those obtained from model outputs (Fig. 4). These results are also true for detrended times series, indicating that this modeled link is valid at the interannual timescale (not shown).

Figure 4The 5-yearly correlations (r) between SMB and surface air temperature for the seven Antarctic regions (see Fig. 1 for geographical definitions) for GCMs over 1850–2000 CE (a), for RACMO2 over 1979–2016 CE (b), and for ice core reconstructions (Thomas et al.2017; Stenni et al.2017a) for 1850–2000 CE (c). For the CMIP5 models and RACMO2, their correlations are all statistically significant (p value < 0.05). For the reconstructions, the statistically significant (p value < 0.05) correlations are obtained for the Antarctic Peninsula and Dronning Maud Land coast. See Fig. S2 for the correlations for each CMIP5 model.

Figure 5SMB sensitivity to surface air temperature (SAT) over the 850–1850 and the 1850–2005 periods for each Antarctic region (see Fig. 1 for geographical definitions) for GCM outputs. Additionally, the SMB–δ18O sensitivity for ice-core-based reconstructions (i.e., Thomas et al.2017; Stenni et al.2017a) over 1850–2005 is represented by a solid black vertical line while a solid red vertical line represents the SMB-observed surface air temperature sensitivity (i.e., Thomas et al.2017; Nicolas and Bromwich2014a) over 1960–2010. For the CESM1–CAM5 model, the 12 simulations are plotted as grey points.


To quantify more precisely the link between surface air temperature and SMB in model outputs and reconstructions, the SMB sensitivity to temperature – defined as the slope of the linear fit between surface air temperature and SMB – has been calculated. Firstly, the GCMs and reconstructions (i.e., Thomas et al.2017; Nicolas and Bromwich2014a) suggest that this sensitivity is similar for both West Antarctica and East Antarctica over the 1950–2000 period (Fig. 3). Secondly, on average over the entire continent, this sensitivity reaches 3.6 % K−1 in ice-core-based reconstructions for the 1850–1949 period. According to these reconstructions, the SMB sensitivity to temperature has increased a lot for the recent period (1950–2005; 15.52 % K−1), confirming the findings of Frieler et al. (2015). However, Frieler et al. (2015) do not obtain such an increase in SMB sensitivity (only +40 %). Additionally, this recent increase found here in the reconstructions is not represented by the GCMs: on average, the simulated sensitivity of SMB to surface air temperature is 5.0±1.1 % K−1 over 1850–1949 and 5.4±2.0 % K−1 over 1950–2005. When looking at the regional scale over 1850–2005, the average SMB sensitivity over all models for West Antarctica (6.8 % K−1) is in good agreement with the one deduced from the reconstructions (8.0 % K−1; Fig. 5), while for East Antarctica, the sensitivity of the model mean is higher than the one obtained from the reconstructions (6.2 % K−1 and 2.1 % K−1, respectively). The very low SMB sensitivity in the reconstructions for East Antarctica, especially on the Antarctic Plateau (0.5 % K−1), is somewhat unexpected, given that this region is continental and thus less affected by synoptic activities than coastal areas (Monaghan and Bromwich2008). Actually, when using the observed surface air temperatures (e.g., Nicolas and Bromwich2014a) instead of the reconstructed ones of Stenni et al. (2017a), the Antarctic SMB sensitivity to temperature is strongly reduced (4.02 % K−1 for the 1958–2010 period) and thus closer to the resulting sensitivity found in the GCMs (5.4±2.0 % K−1 for the 1950–2005 period).

Figure 6The 5-year correlations between SMB and δ18O, surface air temperature and δ18O, and SMB and surface air temperature for the seven Antarctic regions over the 1850–1995 period from the iHadCM3 outputs. The error bars correspond to the range (maximum and minimum) of the iHadCM3 simulations while the dot is the mean of the simulation ensemble is shown. In black circles, the correlation between the SMB ice core reconstructions from Thomas et al. (2017) and the δ18O of Antarctic ice cores aggregated for the seven Antarctic regions (Stenni et al.2017a). In black squares, the correlation between the SMB reconstructions from Thomas et al. (2017) and the observed surface air temperatures aggregated for the seven Antarctic regions (Nicolas and Bromwich2014a) is shown. This latter dataset only covers the 1960–2010 period (5-year averages). Non-significant correlations (p value >= 0.05) are shown in pale.


In the study of Neukom et al. (2018), the authors argue that the data sampling, the noise in proxy data, and the deficiencies in the reconstruction methods can partly explain the discrepancy between models and reconstructions for the surface air temperature during the last millennium, especially for the Southern Hemisphere. The spatial coverage of the surface air temperature and SMB reconstructions based on ice cores is poor, in particular for East Antarctica (Stenni et al.2017a; Thomas et al.2017). Moreover, due to the low snow accumulation in some regions, the uncertainties of the reconstructions are large for both surface air temperatures and SMB, leading to noise in the time series (Stenni et al.2017a; Thomas et al.2017; Frezzotti et al.2007). Since the SMB reconstruction is only based on direct snow accumulation measurements, this is expected to be more accurate than the δ18O-based temperature reconstruction, which is built by assuming a stationary link between δ18O and surface air temperature. Because a lot of processes (such as precipitation seasonality or moisture origin) can significantly modify this relationship over time (e.g., Jouzel et al.1997; Sime et al.2008), this is computed over a short calibration period, but this might be too short to be representative (Klein et al.2019). Consequently, all these processes could explain the large differences between models and proxy-based reconstructions in the estimation of SMB sensitivity to surface air temperature.

4.3 SMB and surface air temperature reconstructions from data assimilation

The high correlation values obtained between SMB and surface air temperatures in GCMs suggest that we can potentially use SMB to reconstruct Antarctic surface air temperature. The analysis of isotope-enabled model results reinforces this hypothesis (Fig. 6): the iHadCM3 outputs show high correlations between these two variables. For most regions, the link between surface air temperature and SMB (r=0.70 on average over the seven regions for the 1850–2000 period) is higher than that between surface air temperatures and δ18O (r=0.55 on average over the seven regions for the 1850–2000 period). This is consistent with the observations: the regional correlations between SMB from ice cores (e.g., Thomas et al.2017) and the observed surface air temperatures (i.e., Nicolas and Bromwich2014a) are high for several regions over the 1960–2010 period (using 5-year averages as for Stenni et al.2017a). In particular, this correlation for East Antarctica is 0.82 (statistically significant). The results with the outputs of ECHAM5-wiso and ECHAM5–MPI-OM are a bit more nuanced than those from iHadCM3 (Fig. S3). The results of ECHAM5-wiso and ECHAM5–MPI-OM confirm this strong link between SMB and temperature, but, in contrast to iHadCM3, the correlations are not systematically higher than between δ18O and temperature. When analyzing the long ECHAM5–MPI-OM simulation (800–2000 CE), the relationship between SMB and surface air temperature is generally higher than between δ18O and surface air temperature but the difference is small. For some regions, the SMB–surface air temperature link is much higher than the δ18O–surface air temperature link but it is weaker for other regions. Compared to the δ18O–surface air temperature link, the SMB–surface air temperature is also less spatially variable (minimum regional correlation is 0.54 compared with 0.07 for the δ18O–surface air temperature link).

Figure 7Reconstructed temperatures (5-year mean) for West Antarctica, East Antarctica, and Antarctica as a whole from the data assimilation experiment (red) using the iHadCM3 outputs and the δ18O (Stenni et al.2017a) and SMB reconstructions (Thomas et al.2017) as constraints in the data assimilation process. The period covered is 1800–2010. The surface air temperature reconstructions of Stenni et al. (2017a) are represented in black and those from Nicolas and Bromwich (2014a) are in blue. DA δ18O (a–c) is the data assimilation experiment using only the δ18O data to constrain the model while DA SMB (d–f) uses the SMB reconstruction and DA δ18O and SMB (g–i) use both. For each experiment and each region, the correlation (r) between the reconstruction based on ice cores (in black) and that based on data assimilation is computed (in red). The shaded areas represent ±1 standard deviation of the model particles. Asterisks represent the statistically significant correlations (p value < 0.05).


Stenni et al. (2017a)

Table 1The 5-year mean correlations between the three surface air temperature reconstructions from data assimilation experiments using the ECHAM5–MPI-OM, ECHAM5-wiso, and iHadCM3 outputs and the statistical reconstruction of Stenni et al. (2017a), with the surface air temperature reconstructions from Nicolas and Bromwich (2014a) over the 1958–2010 period for West Antarctica, East Antarctica, and Antarctica as a whole. Asterisks represent statistically significant correlations (p value < 0.10). DA stands for data assimilation.

Download Print Version | Download XLSX

4.3.1 Surface air temperature reconstruction

When constraining the model with the SMB reconstruction of Thomas et al. (2017), the obtained surface air temperature reconstruction is less well correlated with the reconstruction of Stenni et al. (2017a) than for the data assimilation reconstruction constrained by only the δ18O (Fig. 7). However, the difference is relatively small, despite the fact that SMB and surface air temperatures are more strongly correlated in models than in the ice core reconstruction (0.86 for iHadCM3 compared with 0.16 for ice cores; Fig. 6). When compared to observed surface air temperature over the 1958–2010 period (i.e., Nicolas and Bromwich2014a), the surface air temperature reconstruction of Stenni et al. (2017a) as well as the reconstruction when only δ18O is assimilated are in good agreement with the observed surface air temperatures for West Antarctica (Table 1, correlation coefficients are 0.79 and 0.69, respectively, both statistically significant) but not for East Antarctica (coefficient correlations are 0.10 and 0.13, respectively, both not statistically significant).

In contrast to the data assimilation experiment in which only δ18O is assimilated, the skill of the surface air temperature reconstruction is almost identical for both regions in the data assimilation experiment where only SMB is assimilated: r=0.55 (p value < 0.1) for West Antarctica and r=0.60 (p value < 0.1) for East Antarctica. Assimilating SMB thus provides a more spatially robust temperature reconstruction than when assimilating δ18O. When both δ18O and SMB are taken into account in the data assimilation process, the skill of the surface air temperature reconstructions for the two sub-Antarctic continents is higher (r=0.72 and 0.61 for West Antarctica and for East Antarctica, respectively, both significant) than when separately assimilating the δ18O or the SMB. Moreover, the only reconstruction that provides statistically significant results for all the sub-Antarctic continents (West, East, and Antarctica as a whole; p value < 0.1) is when both δ18O and SMB are assimilated, implying that assimilating both proxies offers more robust results than only assimilating one of them.

When looking at the linearly detrended time series, our final reconstruction (i.e., when δ18O and SMB are assimilated) displays a null correlation with observed surface air temperature (p value = 0.99) for West Antarctica, but the correlation remains high for East Antarctica (r=0.60; p value = 0.07). During the 1958–2012 period, a significant warming is observed in West Antarctica while no significant change is noticed for East Antarctica (Nicolas and Bromwich2014a). Consequently, data assimilation tends to reproduce the warming for West Antarctica and the interannual variability for East Antarctica, explaining our different results between the original and detrended time series. Additionally, as well as our reconstruction based on only δ18O, the correlation of the detrended δ18O-based temperature reconstruction of Stenni et al. (2017a) with the observed one for East Antarctica is non-significant and negative, suggesting that SMB constitutes a better proxy than δ18O for surface air temperatures, at least at the interannual timescale (see Table S2).

Table 2SMB trends over grounded (i.e., excluding ice shelves) West Antarctica, East Antarctica, and Antarctica as a whole from (1) our reconstruction based on data assimilation using iHadCM3 outputs and SMB and δ18O data in the data assimilation procedure; (2) Medley and Thomas (2019); (3) RACMO2 outputs for various time intervals (Gt yr−2). Asterisks stand for statistically significant trends at the 5 % level.

Download Print Version | Download XLSX

Regarding surface air temperature trends over the last 2 centuries, our reconstruction displays an increase of 0.02 C per decade for West Antarctica and 0.023 C per decade for East Antarctica, which finally leads to an increase of 0.022 C per decade for Antarctica as a whole (all statistically significant). For the 1961–2010 period, our reconstruction is able to simulate the observed contrast between West and East Antarctica (0.22 C per decade (significant) and 0.053 C per decade (not significant), respectively, for Nicolas and Bromwich (2014a) compared with 0.1 C per decade and 0.06 C per decade, respectively, for our reconstruction, both significant). The resulting contrast in our reconstruction is thus not as large as observed (see Table S3 for details). However, because of the short time period considered, these values can be highly variable depending on the time interval chosen (not shown).

Figure 8Reconstructed SMB (5-year mean) for West Antarctica, East Antarctica, and Antarctica as a whole from the data assimilation experiment using the iHadCM3 outputs and δ18O (Stenni et al.2017a) and SMB reconstruction (in black; Thomas et al.2017) as constraints in the data assimilation process. The period is 1800–2010. DA δ18O (a–c) is the data assimilation experiment using only the δ18O data to constrain the model while DA SMB (d–f) uses the SMB reconstruction and DA δ18O and SMB (g–i) uses both. For each experiment and each region, the correlation (r) between the reconstruction based on ice cores (in black) and that based on data assimilation is computed (in red). The shaded areas represent ±1 standard deviation of the model particles. Asterisks represent the statistically significant correlations (p value < 0.05).


4.3.2 SMB reconstruction

Constraining the model with the δ18O data leads to a poor SMB reconstruction, especially for West Antarctica (correlation coefficient of 0.29; Fig. 8). Moreover, the constraint derived from observed δ18O on SMB is weak as illustrated by the large error band of the reconstruction (estimated by the weighted variance of the particles with nonzero weight). When assimilating both δ18O and SMB, the SMB reconstruction is in good agreement with the reconstruction of Thomas et al. (2017) as expected.

According to this data-assimilation-based SMB reconstruction, the grounded AIS SMB (i.e., SMB over the AIS excluding ice shelves) increased at a pace of 0.26 Gt yr−2 (p value < 0.05) during the 1801–2000 period (Table 2) and 0.70 Gt yr−2 (p value > 0.05) for the 1957–2000 period. Over this latter period, West Antarctica witnessed a significant increase of 0.82 Gt yr−2 while East Antarctica displayed a decrease of 0.13 Gt yr−2 (p value > 0.05). Unlike West Antarctica, the non-statistical significance of the SMB trend for East Antarctica might imply that internal variability currently plays a large role in the SMB variability there (e.g., Jones et al.2016). However, if we focus on the shorter 1979–2000 period, a significant decrease is obtained for East Antarctica (−3.3 Gt yr−2; p value < 0.05) while it is still positive for West Antarctica (1.6 Gt yr−2; p value > 0.05), which is consistent with RACMO2 outputs (−3.7 Gt yr−2 for East Antarctica and 2.0 Gt yr−2 for West Antarctica, both not significant).

5 Discussion and conclusions

This paper discusses the AIS SMB over the last 2 centuries and its links with surface air temperature in reconstructions and model simulations. The analysis of the relationship between SMB and surface air temperature in models and in ice core reconstructions highlighted the covariance between both variables that can potentially be used to reconstruct past changes. The relevance of SMB in the reconstruction of surface air temperature in Antarctica is based on a relatively simple concept: Antarctic precipitation originates mainly from lower latitudes, in the form of warm and wet air masses (Goodwin et al.2016; Turner et al.2016; Clem et al.2018). Nevertheless, δ18O also provides useful temperature-related information that can be used to complement the information provided by SMB, such as changes in moisture origin (e.g., Holloway et al.2016a). Our analyses pointed out significant model–data discrepancies in the SMB–surface air temperature relationship. On the one hand, models show a strong correlation between δ18O and SMB for many Antarctic regions, and, on the other hand, the reconstructions based on ice cores display a weak relationship. Furthermore, unlike previous studies (e.g., Frieler et al.2015) that suggest an increase in the SMB sensitivity to surface air temperature for the future in Antarctica ( 40 %), we show that the current sensitivity is not exceptionally high compared to the last 200 years, according to CMIP5 models.

These large discrepancies between model results and reconstructions can be explained by different factors. The GCMs may have biases in the simulated temperature changes. For example, as shown by Klein et al. (2019), GCMs display on average a homogeneous warming over Antarctica during the last decades while observations mainly show warming for West Antarctica with no significant change for East Antarctica. Additionally, climate model simulations generally display a warming starting in the 19th century in Antarctica while it begins much later in proxy-based reconstructions (Abram et al.2016). This suggests that reconstructions underestimate the response to anthropogenic forcing or that climate models overestimate it. In this latter case, this may contribute to an overestimation of the contribution of the simple thermodynamic link between temperature and precipitation and thus snow accumulation while it underestimates the role of changes in atmospheric circulation variability (Abram et al.2016; Klein et al.2019; PAGES 2k-PMIP3 group2015). Nevertheless, by removing the linear trend of time series, we obtained similar results. Models may also neglect processes such as blowing snow that can reduce the correlation between temperature and SMB. On the other hand, RACMO2, which includes a simple representation of blowing snow and is nudged by observed temperature and large-scale circulation changes, displays similar correlations to those of the GCMs. Another hypothesis is that the differences could rather arise from uncertainties in the reconstructions. To understand the potential origin of the disagreements between model results and reconstructions over the last millennium, Neukom et al. (2018) used pseudoproxy experiments. They found that uncertainties in the reconstructions and the data sampling could be an explanation for many observed discrepancies between models and reconstructions.

By analyzing isotope-enabled climate models, we showed that on average over the models the relationship between SMB and surface air temperature is often higher (or at least equivalent) and more stable than the one between surface air temperature and δ18O. This is true on both the continental and regional scales. Unlike SMB, δ18O can be subject to large uncertainties linked to precipitation seasonality (Sime et al.2008) or changes in moisture origins (Holloway et al.2016a), which can explain the weaker correlations.

Our data assimilation experiments confirm the benefits of using both proxies – SMB and δ18O – to reconstruct surface air temperature. When assimilating both δ18O and SMB data, the resulting reconstruction shows a higher correlation with observed surface air temperature over the period 1958–2010 (i.e., Nicolas and Bromwich2014a) for the entire Antarctic continent (r=0.73) than the one obtained with the reconstruction based on the statistical method of Stenni et al. (2017ar=0.57). The difference is larger for East Antarctica, where the reconstruction skill is enhanced by incorporating SMB data (r=0.61 for our reconstruction compared with 0.10 for the reconstruction of Stenni et al.2017a). For West Antarctica, our reconstruction is very similar to the statistical method of Stenni et al. (2017a). This improvement can be explained by the large uncertainties in δ18O data for East Antarctica, probably because of the low number of ice cores and low snow accumulation in those areas. In comparison to Stenni et al. (2017a) and Klein et al. (2019), who obtain a higher surface air temperature trend over the last 2 centuries for East Antarctica (0.03 and 0.018 C per decade, respectively, both significant) than for West Antarctica (0.011 and 0.01 C per decade, respectively, both not significant), our data-assimilation-based reconstruction reveals similar surface air temperature trends for both regions (0.02 and 0.023 C per decade, respectively, both significant). However, over the entire continent, the trend is almost the same between the different datasets (0.022 C per decade in this study (significant), 0.019 C per decade for Stenni et al.2017a, significant, and 0.016 C per decade for Klein et al.2019, not significant). Over the last decades (1961–2010), all the reconstructions are able to reproduce the observed contrast between West Antarctica (large warming) and East Antarctica (weak warming), but overall they underestimate this regional contrast (see Table S3 for details). Our reconstruction displays smaller variance in time than the reconstruction from Stenni et al. (2017a), which is a standard characteristic of estimates based on data assimilation using only a few uncertain data (e.g., Goosse et al.2010).

Regarding changes in SMB over the last 2 centuries, our reconstruction shows large regional differences in SMB trends, both in magnitude and in signs, in accordance with Medley and Thomas (2019, Fig. S4), who used the same ice core dataset but a different method. While they obtain a statistically significant SMB increase of 0.4 Gt yr−2 over the grounded AIS for 1801–2000, our result suggests a weaker increase (0.26 Gt yr−2; p value < 0.001; see Table 2 for details). A similar underestimation is noticed for the 1957–2000 period, (1.0 Gt yr−2 for Medley and Thomas2019, not significant, and 0.70 Gt yr−2 for our reconstruction, p value = 0.130). Over the last decades (1979–2000), both Medley and Thomas (2019) and our results reveal that grounded West Antarctica gains mass at its surface (1.6 Gt yr−2 in this study and 1.7 Gt yr−2 for Medley and Thomas2019, both not significant) while grounded East Antarctica has experienced a very large SMB decrease (−3.3 and −4.5 Gt yr−2, respectively, both significant), which is consistent with the value from the RACMO2 outputs (2.0 Gt yr−2 for West Antarctica and −3.7 Gt yr−2 for East Antarctica, both not significant).

More generally, in contrast to statistical methods, data assimilation ensures that reconstructions are compatible with the physics of the system as represented in the models chosen. Although it is not possible to independently evaluate our SMB reconstruction, our good results regarding surface air temperatures and SMB reconstructions suggest that the strong simulated correlation between surface air temperatures and SMB in GCMs is not a model artifact. This is supported by a strong link between these two variables in observations when using snow accumulation data from Thomas et al. (2017) and surface air temperatures from Nicolas and Bromwich (2014a), in particular for East Antarctica (r=0.82, statistically significant). Therefore, our study shows that SMB records seem to be a relevant proxy in reconstructing surface air temperature as a complement to δ18O records. Since only a few records are available before the instrumental period over Antarctica, any relevant record to reconstruct the Antarctic climate and more specifically surface air temperature is welcome. Additionally, data assimilation appears particularly well adapted for reconstructing surface air temperatures as the covariance between variables is obtained directly from climate models that explicitly include physical processes while statistical approaches restrict the problem to empirical linear relationships. By using data assimilation, no assumption such as stationarity or long calibration periods is required to estimate the link between variables, assumptions whose validity can strongly vary in time and space (Klein et al.2019). However, to get a skillful data assimilation-based reconstruction, it is essential that the selected climate models have an adequate representation of climate variability and that good uncertainty estimates are available for the chosen datasets.

Code and data availability

The resulting Antarctic SMB and surface air temperature reconstructions are available upon request to Quentin Dalaiden ( All CMIP5–PMIP3 model simulations can be directly downloaded on (last access: 10 February 2018). iHadCM3 data are available by request to Max Holloway ( ECHAM5-wiso data covering the 1871–2011 period can be downloaded from (Steiger2018). Products from the ECHAM5–MPI-OM model simulation are available by request to Jesper Sjolte ( RACMO2 data are available by request to Jan Lenaerts (jan.lenaerts@Colorado.EDU). δ18O, surface air temperature, and SMB reconstructions are stored at the UK Polar Data Centre and at the NOAA World Data Center for Paleoclimatology (, last access: 18 March 2018, Stenni et al.2017b), or they can be requested from Elizabeth R. Thomas ( Antarctic observed surface air temperatures are available at (last access: 4 July 2018, Nicolas and Bromwich2014b).


The supplement related to this article is available online at:

Author contributions

This study is part of QD's thesis under the supervision of HG. QD and HG designed the study. QD performed the analysis and made the figures. FK contributed to the preparation, launch, and interpretation of the data assimilation experiments. JTML provided RACMO data and advice on their use. MH and LS provided iHadCM3 outputs and contributed to their interpretation in the framework of this study. ERT provided the PAGES2k databases and advice on the climate proxy interpretation. QD led the writing of the manuscript with contributions from all authors.

Competing interests

The authors declare no competing interests.


We would like to thank Nathan Steiger and Jesper Sjolte for the isotopic model outputs. We would also like to thank Marie Cavitte for her feedback on our study. We acknowledge the World Climate Research Programme Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modeling groups for producing and making available their model outputs. We sincerely thank the two reviewers and the editor for their constructive comments, which have improved our manuscript. This work was supported by the Belgian Research Action through Interdisciplinary Networks (BRAIN-be) from Belgian Science Policy Office in the framework of the project “East Antarctic surface mass balance in the Anthropocene: observations and multiscale modelling (Mass2Ant)” (contract no. BR/165/A2/Mass2Ant). Quentin Dalaiden is a research fellow with the Fonds pour la formation á la Recherche dans l'Industrie et dans l'Agronomie (FRIA-Belgium) and Hugues Goosse is the research director within the F.R.S.-FNRS.

Financial support

This research has been supported by the BRAIN-be (grant no. BR/165/A2/Mass2Ant).

Review statement

This paper was edited by Kenichi Matsuoka and reviewed by two anonymous referees.


Abram, N. J., McGregor, H. V., Tierney, J. E., Evans, M. N., McKay, N. P., Kaufman, D. S., Thirumalai, K., Martrat, B., Goosse, H., Phipps, S. J., Steig, E. J., Kilbourne, K. H., Saenger, C. P., Zinke, J., Leduc, G., Addison, J. A., Mortyn, P. G., Seidenkrantz, M. S., Sicre, M. A., Selvaraj, K., Filipsson, H. L., Neukom, R., Gergis, J., Curran, M. A., and Von Gunten, L.: Early onset of industrial-era warming across the oceans and continents, Nature, 536, 411–418,, 2016. a, b, c

Agosta, C., Amory, C., Kittel, C., Orsi, A., Favier, V., Gallée, H., van den Broeke, M. R., Lenaerts, J. T. M., van Wessem, J. M., van de Berg, W. J., and Fettweis, X.: Estimation of the Antarctic surface mass balance using the regional climate model MAR (1979–2015) and identification of dominant processes, The Cryosphere, 13, 281–296,, 2019. a

Altnau, S., Schlosser, E., Isaksson, E., and Divine, D.: Climatic signals from 76 shallow firn cores in Dronning Maud Land, East Antarctica, The Cryosphere, 9, 925–944,, 2015. a, b

Bamber, J. L., Gomez-Dans, J. L., and Griggs, J. A.: A new 1 km digital elevation model of the Antarctic derived from combined satellite radar and laser data – Part 1: Data and methods, The Cryosphere, 3, 101–111,, 2009. a

Bamber, J. L., Westaway, R. M., Marzeion, B., and Wouters, B.: The land ice contribution to sea level during the satellite era, Environ. Res. Lett., 13, 063008,, 2018. a

Church, J., Clark, P., Cazenave, A., Gregory, J., Jevrejeva, S., Levermann, A., Merrifield, M., Milne, G., Nerem, R., Nunn, P., Payne, A., Pfeffer, W., Stammer, D., and Unnikrishnan, A. S.: The Physical Science Basism Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Stocker, T. F., Qin, D., Plattner, G.-K., Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex V., and Midg, P. M., Tech. rep., Cambridge University Press, Cambridge, UK, New York, NY, USA, 2013. a

Clem, K. R., Renwick, J. A., and McGregor, J.: Autumn Cooling of Western East Antarctica Linked to the Tropical Pacific, J. Geophys. Res.-Atmos., 123, 89–107,, 2018. a

Danielson, J. and Gesch, D.: Global Multi-resolution Terrain Elevation Data 2010(GMTED2010), U.S. Geological Survey, 2010, 26 pp., 2011. a

Dansgaard, W. and Johnsen, S. J.: A Flow Model and a Time Scale for the Ice Core from Camp Century, Greenland, J. Glaciol., 8, 215–223,, 1969. a

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, I., Biblot, J., Bormann, N., Delsol, C., Dragani, R., Fuentes, M., Greer, A. J., Haimberger, L., Healy, S. B., Hersbach, H., Holm, E. V., Isaksen, L., Kallberg, P., Kohler, M., Matricardi, M., McNally, A. P., Mong-Sanz, B. M., Morcette, J.-J., Park, B.-K., Peubey, C., de Rosnay, P., Tavolato, C., Thepaut, J. N., and Vitart, F.: The ERA-Interim reanalysis: Configuration and performance of the data assimilation system, Q. J. Roy. Meteorol. Soc., 137, 553–597,, 2011. a

Dubinkina, S., Goosse, H., Sallaz-Damaz, Y., Crespin, E., and Crucifix, M.: Testing a Particle Filter To Reconstruct Climate Changes Over the Past Centuries, Int. J. Bifurcat. Chaos, 21, 3611–3618,, 2011. a, b, c, d

Dufresne, J.-L., Foujols, M.-A., Denvil, S., Caubel, A., Marti, O., Aumont, O., Balkanski, Y., Bekki, S., Bellenger, H., Benshila, R., Bony, S., Bopp, L., Braconnot, P., Brockmann, P., Cadule, P., Cheruy, F., Fairhead, L., Fichefet, T., Codron, F., Cozic, A., Cugnet, D., Noblet, N. D., Ethe, C., Flavoni, S., Friedlingstein, P., Guez, L., Guilyardi, E., Hauglustaine, D., Hourdin, F., Idelkadi, A., Ghattas, J., Joussaume, S., Kageyama, M., Krinner, G., Labetoulle, S., Lahellec, A., Lefevre, F., Levy, C., Li, Z. X., Lloyd, J., Lott, F., Madec, G., Mancip, M., Marchand, M., Masson, S., Meurdesoif, Y., Mignot, J., Musat, I., Parouty, S., Polcher, J., Rio, C., Schulz, M., Swingedouw, D., Szopa, S., Talandier, C., Terray, P., Viovy, N., and Vuichard, N.: Climate change projections using the IPSL-CM5 Earth System Model: from CMIP3 to CMIP5, Clim. Dynam., 40, 2123–2165,, 2013. a

ECMWF: Part IV : Physical Processes (CY33R1), Tech. rep., available at: (last access: 3 April 2019), 2008. a

Ettema, J., van den Broeke, M. R., van Meijgaard, E., van de Berg, W. J., Box, J. E., and Steffen, K.: Climate of the Greenland ice sheet using a high-resolution climate model – Part 1: Evaluation, The Cryosphere, 4, 511–527,, 2010. a

Fan, T., Deser, C., and Schneider, D.: Recent Antarctic sea ice trends in thecontext of Southern Ocean surface climate variations since 1950, Geophys. Res. Lett., 41, 2419–2426,, 2014. a

Favier, V., Agosta, C., Parouty, S., Durand, G., Delaygue, G., Gallée, H., Drouet, A.-S., Trouvilliez, A., and Krinner, G.: An updated and quality controlled surface mass balance dataset for Antarctica, The Cryosphere, 7, 583–597,, 2013. a

Frezzotti, M., Urbini, S., Proposito, M., Scarchilli, C., and Gandolfi, S.: Spatial and temporal variability of surface mass balance near Talos Dome, East Antarctica, J. Geophys. Res.-Earth, 112, F02032,, 2007. a

Frezzotti, M., Scarchilli, C., Becagli, S., Proposito, M., and Urbini, S.: A synthesis of the Antarctic surface mass balance during the last 800 yr, The Cryosphere, 7, 303–319,, 2013. a, b

Frieler, K., Clark, P. U., He, F., Buizert, C., Reese, R., Ligtenberg, S. R., Van Den Broeke, M. R., Winkelmann, R., and Levermann, A.: Consistent evidence of increasing Antarctic accumulation with warming, Nat. Clim. Change, 5, 348–352,, 2015. a, b, c, d, e, f

Fudge, T. J., Markle1, B. R., Cuffey, K. M., Buizert, C., Taylor4, K. C., Steig, E. J., Waddington, E. D., Conway, H., and Koutnik, M.: Variable relationship between accumulation and temperature in West Antarctica for the past 31,000 years, Geophys. Res. Lett., 43, 3795–3803,, 2016. a

Fyke, J., Sergienko, O., Löfverström, M., Price, S., and Lenaerts, J. T.: An Overview of Interactions and Feedbacks Between Ice Sheets and the Earth System, Rev. Geophys., 56, 361–408,, 2018. a

Gent, P. R., Danabasoglu, G., Donner, L. J., Holland, M. M., Hunke, E. C., Jayne, S. R., Lawrence, D. M., Neale, R. B., Rasch, P. J., Vertenstein, M., Worley, P. H., Yang, Z. L., and Zhang, M.: The community climate system model version 4, J. Climate, 24, 4973–4991,, 2011. a

Goodwin, B. P., Mosley-Thompson, E., Wilson, A. B., Porter, S. E., and Roxana Sierra-Hernandez, M.: Accumulation variability in the Antarctic Peninsula: The role of large-scale atmospheric oscillations and their interactions, J. Climate, 29, 2579–2596,, 2016. a, b

Goosse, H.: Reconstructed and simulated temperature asymmetry between continents in both hemispheres over the last centuries, Clim. Dynam., 48, 1483–1501,, 2017. a

Goosse, H., Crespin, E., de Montety, A., Mann, M. E., Renssen, H., and Timmermann, A.: Reconstructing surface temperature changes over the past 600 years using climate model simulations with data assimilation, J. Geophys. Res., 115, D09108,, 2010. a, b

Goosse, H., Braida, M., Crosta, X., Mairesse, A., Masson-Delmotte, V., Mathiot, P., Neukom, R., Oerter, H., Philippon, G., Renssen, H., Stenni, B., van Ommen, T., and Verleyen, E.: Antarctic temperature changes during the last millennium: evaluation of simulations and reconstructions, Quatern. Sci. Rev., 55, 75–90,, 2012. a

Goursaud, S., Masson-Delmotte, V., Favier, V., Orsi, A., and Werner, M.: Water stable isotope spatio-temporal variability in Antarctica in 1960–2013: observations and simulations from the ECHAM5-wiso atmospheric general circulation model, Clim. Past, 14, 923–946,, 2018. a

Goursaud, S., Masson-Delmotte, V., Favier, V., Preunkert, S., Legrand, M., Minster, B., and Werner, M.: Challenges associated with the climatic interpretation of water stable isotope records from a highly resolved firn core from Adélie Land, coastal Antarctica, The Cryosphere, 13, 1297–1324,, 2019. a

Greuell, W. and Thomas, K.: Numerical modelling of the energy balance and the englacial temperature of the Greenland Ice Sheet. Calculations for the ETH-Camp location (West Greenland, 1155 m a.s.l.), Global Planet. Change, 9, 91–114, 1994. a

Hakim, G. J., Emile-Geay, J., Steig, E. J., Noone, D., Anderson, D. M., Tardif, R., Steiger, N., and Perkins, W. A.: The last millennium climate reanalysis project: Framework and first results, J. Geophys. Res., 121, 6745–6764,, 2016. a, b

Holloway, M. D., Sime, L. C., Singarayer, J. S., Tindall, J. C., Bunch, P., and Valdes, P. J.: Antarctic last interglacial isotope peak in response to sea ice retreat not ice-sheet collapse, Nat. Commun., 7, 12293,, 2016a. a, b, c

Holloway, M. D., Sime, L. C., Singarayer, J. S., Tindall, J. C., and Valdes, P. J.: Reconstructing paleosalinity from δ18O: Coupled model simulations of the Last Glacial Maximum, Last Interglacial and Late Holocene, Quatern. Sci. Rev., 131, 350–364,, 2016b. a

Huybrechts, P., Gregory, J., Janssens, I., and Wild, M.: Modelling Antarctic and Greenland volume changes during the 20th and 21st centuries forced by GCM time slice integrations, Global Planet. Change, 42, 83–105,, 2004. a

Jones, J. M., Gille, S. T., Goosse, H., Abram, N. J., Canziani, P. O., Charman, D. J., Clem, K. R., Crosta, X., De Lavergne, C., Eisenman, I., England, M. H., Fogt, R. L., Frankcombe, L. M., Marshall, G. J., Masson-Delmotte, V., Morrison, A. K., Orsi, A. J., Raphael, M. N., Renwick, J. A., Schneider, D. P., Simpkins, G. R., Steig, E. J., Stenni, B., Swingedouw, D., and Vance, T. R.: Assessing recent trends in high-latitude Southern Hemisphere surface climate, Nat. Clim. Change, 6, 917–926,, 2016. a, b, c

Jouzel, J.: Magnitude of isotope/temperature scaling for interpretation of central Antarctic ice cores, J. Geophys. Res., 108, 4361,, 2003. a

Jouzel, J., Alley, R. B., Cuffey, K. M., Dansgaard, W., Grootes, P., Hoffmann, G., Johnsen, S. J., Koster, R. D., Peel, D., Shuman, C. A., Stievenard, M., Stuiver, M., and White, J.: Validity of the temperature reconstruction from ice cores, J. Geophys. Res., 102, 26471–26487, 1997. a, b

Klein, F. and Goosse, H.: Reconstructing East African rainfall and Indian Ocean sea surface temperatures over the last centuries using data assimilation, Clim. Dynam., 50, 3909–3929,, 2018. a, b, c, d

Klein, F., Abram, N. J., Curran, M. A. J., Goosse, H., Goursaud, S., Masson-Delmotte, V., Moy, A., Neukom, R., Orsi, A., Sjolte, J., Steiger, N., Stenni, B., and Werner, M.: Assessing the robustness of Antarctic temperature reconstructions over the past 2 millennia using pseudoproxy and data assimilation experiments, Clim. Past, 15, 661–684,, 2019. a, b, c, d, e, f, g, h, i, j, k, l

Krinner, G., Magand, O., Simmonds, I., Genthon, C., Dufresne, J.-l., Krinner, G., Magand, O., Simmonds, I., Genthon, C., and Simulated, J.-l. D.: Simulated antarctic precipitation and surface mass balance of the end of the 20th and 21st centuries To cite this version : HAL Id : hal-00184741 Simulated Antarctic precipitation and surface mass balance at, Clim. Dynam., 28, 215–230,, 2007. a

Kuipers Munneke, P., Picard, G., Van Den Broeke, M. R., Lenaerts, J. T., and Van Meijgaard, E.: Insignificant change in Antarctic snowmelt volume since 1979, Geophys. Res. Lett., 39, 6–10,, 2012. a, b

Laepple, T., Münch, T., Casado, M., Hoerhold, M., Landais, A., and Kipfstuhl, S.: On the similarity and apparent cycles of isotopic variations in East Antarctic snow pits, The Cryosphere, 12, 169–187,, 2018. a

Lenaerts, J. T., Van Meijgaard, E., Van Den Broeke, M. R., Ligtenberg, S. R., Horwath, M., and Isaksson, E.: Recent snowfall anomalies in Dronning Maud Land, East Antarctica, in a historical and future climate perspective, Geophys. Res. Lett., 40, 2684–2688,, 2013. a

Lenaerts, J. T., Vizcaino, M., Fyke, J., van Kampenhout, L., and van den Broeke, M. R.: Present-day and future Antarctic ice sheet climate and surface mass balance in the Community Earth System Model, Clim. Dynam., 47, 1367–1381,, 2016. a

Lenaerts, J. T., Fyke, J. G., and Medley, B.: The signature of ozone depletion in recent Antarctic precipitation change: a study with the Community Earth System Model, Geophys. Res. Lett., 45, 12931–12939,, 2018. a

Lenaerts, J. T., Medley, B., Broeke, M. R., and Wouters, B.: Observing and Modeling Ice‐Sheet Surface Mass Balance, Rev. Geophys., 57, 376–420,, 2019. a

Lenaerts, J. T. M., van den Broeke, M. R., Déry, S. J., König-Langlo, G., Ettema, J., and Munneke, P. K.: Modelling snowdrift sublimation on an Antarctic ice shelf, The Cryosphere, 4, 179–190,, 2010. a

Lenaerts, J. T. M., Van Den Broeke, M. R., Van De Berg, W. J., Van Meijgaard, E., and Kuipers Munneke, P.: A new, high-resolution surface mass balance map of Antarctica (1979–2010) based on regional atmospheric climate modeling, Geophys. Res. Lett., 39, 1–5,, 2012. a, b

Magand, O., Picard, G., Brucker, L., Fily, M., and Genthon, C.: Snow melting bias in microwave mapping of Antarctic snow accumulation, The Cryosphere, 2, 109–115,, 2008. a

Masson-Delmotte, V., Kageyama, M., Braconnot, P., Charbit, S., Krinner, G., Ritz, C., Guilyardi, E., Jouzel, J., Abe-Ouchi, A., Crucifix, M., Gladstone, R. M., Hewitt, C. D., Kitoh, A., LeGrande, A. N., Marti, O., Merkel, U., Motoi, T., Ohgaito, R., Otto-Bliesner, B., Peltier, W. R., Ross, I., Valdes, P. J., Vettoretti, G., Weber, S. L., Wolk, F., and Yu, Y.: Past and future polar amplification of climate change: Climate model intercomparisons and ice-core constraints, Clim. Dynam., 26, 513–529,, 2006. a

Masson-Delmotte, V., Hou, S., Ekaykin, A., Jouzel, J., Aristarain, A., Bernardo, R. T., Bromwich, D., Cattani, O., Delmotte, M. M., Falourd, S., Frezzotti, M., Gallée, H., Genoni, L., Isaksson, E., Landais, A., Helsen, M. M., Hoffmann, G., Lopez, J., Morgan, V., Motoyama, H., Noone, D., Oerter, H., Petit, J. R., Royer, A., Uemura, R., Schmidt, G. A., Schlosser, E., Simões, J. C., Steig, E. J., Stenni, B., Stievenard, M., Van Den Broeke, M. R., Van De Wal, R. S., Van De Berg, W. J., Vimeux, F., and White, J. W.: A review of antarctic surface snow isotopic composition: Observations, atmospheric circulation, and isotopic modeling, J. Climate, 21, 3359–3387,, 2008. a

Matsikaris, A., Widmann, M., and Jungclaus, J.: On-line and off-line data assimilation in palaeoclimatology: a case study, Clim. Past, 11, 81–93,, 2015. a, b

Medley, B. and Thomas, E. R.: Increased snowfall over the Antarctic Ice Sheet mitigated twentieth-century sea-level rise, Nat. Clim. Change, 9, 34–39,, 2019. a, b, c, d, e, f

Medley, B., Joughin, I., Smith, B. E., Das, S. B., Steig, E. J., Conway, H., Gogineni, S., Lewis, C., Criscitiello, A. S., McConnell, J. R., van den Broeke, M. R., Lenaerts, J. T. M., Bromwich, D. H., Nicolas, J. P., and Leuschen, C.: Constraining the recent mass balance of Pine Island and Thwaites glaciers, West Antarctica, with airborne observations of snow accumulation, The Cryosphere, 8, 1375–1392,, 2014. a

Medley, B., McConnell, J. R., Neumann, T. A., Reijmer, C. H., Chellman, N., Sigl, M., and Kipfstuhl, S.: Temperature and Snowfall in Western Queen Maud Land Increasing Faster Than Climate Model Projections, Geophys. Res. Lett., 45, 1472–1480,, 2018. a, b

Monaghan, A. J. and Bromwich, D. H.: Advances in describing recent Antarctic climate variability, B. Am. Meteorol. Soc., 89, 1295–1306,, 2008. a

Monaghan, A. J., Bromwich, D. H., Fogt, R. L., Wang, S.-H., Mayewski, P. A., Dixon, D. A., Ekaykin, A., Frezzotti, M., Goodwin, I., Isaksson, E., Kaspari, S. D., Morgan, V. I., Oerter, H., Van Ommen, T. D., Van der Veen, C. J., and Wen, J.: Insignificant Change in Antarctic Snowfall Since the International Geophysical Year, Science, 313, 827–831,, 2006. a

Neukom, R., Schurer, A. P., Steiger, N. J., and Hegerl, G. C.: Possible causes of data model discrepancy in the temperature history of the last Millennium, Sci. Rep., 8, 1–15,, 2018. a, b, c, d, e

Nicolas, J. P. and Bromwich, D. H.: New reconstruction of antarctic near-surface temperatures: Multidecadal trends and reliability of global reanalyses, J. Climate, 27, 8070–8093,, 2014a. a, b, c, d, e, f, g, h, i, j, k, l, m, n

Nicolas, J. P. and Bromwich, D. H.: Reconstruction of Antarctic near-surface temperatures (1958–2012), Data set, Polar Meterology Group, available at: (last access: 4 July 2018), 2014b. a

Oerter, H., Wilhelms, F., Jung-Rothenhausler, F., Goktas, F., Miller, H., Graf, W., and Sommer, S.: Accumulation rates in Dronning Maud Land, Antarctica, as revealed by dielectric-profiling measurements of shallow firn cores, Ann. Glaciol., 30, 27–34,, 2000. a

Otto-Bliesner, B. L., Joussaume, S., Braconnot, P., Harrison, S. P., and Abe-Ouchi, A.: Modeling and Data Syntheses of Past Climates: Paleoclimate Modelling Intercomparison Project Phase II Workshop, Eos T. Am. Geophys. Un., 90, p. 93, 2009. a

Otto-Bliesner, B. L., Brady, E. C., Fasullo, J., Jahn, A., Landrum, L., Stevenson, S., Rosenbloom, N., Mai, A., and Strand, G.: Climate Variability and Change since 850 CE: An Ensemble Approach with the Community Earth System Model, B. Am. Meteorol. Soc., 97, 735–754,, 2015. a, b

PAGES 2k-PMIP3 group: Continental-scale temperature variability in PMIP3 simulations and PAGES 2k regional temperature reconstructions over the past millennium, Clim. Past, 11, 1673–1699,, 2015. a, b, c

Palerme, C., Genthon, C., Claud, C., Kay, J. E., Wood, N. B., and L'Ecuyer, T.: Evaluation of current and projected Antarctic precipitation in CMIP5 models, Clim. Dynam., 48, 225–239,, 2017. a

Pendergrass, A. G., Hakim, G. J., Battisti, D. S., and Roe, G.: Coupled air-mixed layer temperature predictability for climate reconstruction, J. Climate, 25, 459–472,, 2012. a

Philippe, M., Tison, J.-L., Fjøsne, K., Hubbard, B., Kjær, H. A., Lenaerts, J. T. M., Drews, R., Sheldon, S. G., De Bondt, K., Claeys, P., and Pattyn, F.: Ice core evidence for a 20th century increase in surface mass balance in coastal Dronning Maud Land, East Antarctica, The Cryosphere, 10, 2501–2516,, 2016. a, b

Rayner, N. A., Parker, E. B., Horton, C. K., Holland, C. K., Alexander, L. V., Rowell, D. P., Kent, E. C., and Kaplan, A.: Global analyses of sea surface temperature, sea ice, and night marine air temperature since the late nineteenth century, J. Geophys. Res., 108, 4407,, 2003. a

Rignot, E., Mouginot, J., Scheuchl, B., van den Broeke, M. R., van Wessem, M. J., and Morlighem, M.: Four decades of Antarctic Ice Sheet mass balance from 1979–2017, P. Natl. Acad. Sci. USA, 116, 1095–1103,, 2019. a, b

Rotstayn, L. D., Collier, M. A., Dix, M. R., Feng, Y., Gordon, H. B., O'Farrell, S. P., Smith, I. N., and Syktus, J.: Improved simulation of Australian climate and ENSO-related rainfall variability in a global climate model with an interactive aerosol treatment, Int. J. Climatol., 30, 1067–1088,, 2010. a

Schlosser, E., Anschütz, H., Divine, D., Martma, T., Sinisalo, A., Altnau, S., and Isaksson, E.: Recent climate tendencies on an East Antarctic ice shelf inferred from a shallow firn core network, J. Geophys. Res.-Atmos., 119, 6578–6595,, 2014. a

Schmidt, G. A., Jungclaus, J. H., Ammann, C. M., Bard, E., Braconnot, P., Crowley, T. J., Delaygue, G., Joos, F., Krivova, N. A., Muscheler, R., Otto-Bliesner, B. L., Pongratz, J., Shindell, D. T., Solanki, S. K., Steinhilber, F., and Vieira, L. E. A.: Climate forcing reconstructions for use in PMIP simulations of the Last Millennium (v1.1), Geosci. Model Dev., 5, 185–191,, 2012. a, b

Schmidt, G. A., Jungclaus, J. H., Ammann, C. M., Bard, E., Braconnot, P., Crowley, T. J., Delaygue, G., Joos, F., Krivova, N. A., Muscheler, R., Otto-Bliesner, B. L., Pongratz, J., Shindell, D. T., Solanki, S. K., Steinhilber, F., and Vieira, L. E. A.: Climate forcing reconstructions for use in PMIP simulations of the Last Millennium (v1.1), Geosci. Model Dev., 5, 185–191,, 2012. a

Schmidt, G. A., Kelley, M., Nazarenko, L., Ruedy, R., Russell, G. L., Aleinov, I., Bauer, M., Bauer, S. E., Bhat, M. K., Bleck, R., Canuto, V., Chen, Y.-h., Cheng, Y., Clune, T. L., Genio, A. D., Fainchtein, R. D., Faluvegi, G., Hansen, J. E., Healy, R. J., Kiang, N. Y., Koch, D., Lacis, A. A., Legrande, A. N., Lerner, J., Lo, K. K., Matthews, E. E., Menon, S., Miller, R. L., Oinas, V., and Oloso, A. O.: Journal of Advances in Modeling Earth Systems contributions to the CMIP5 archive, J. Adv. Model. Earth Syst., 6, 141–184,, 2014. a

Shepherd, A., Ivins, E. R., A, G., Barletta, V. R., Bentley, M. J., Bettadpur, S., Briggs, K. H., Bromwich, D. H., Forsberg, R., Galin, N., Horwath, M., Jacobs, S., Joughin, I., King, M. A., Lenaerts, J. T. M., Li, J., Ligtenberg, S. R. M., Luckman, A., Luthcke, S. B., McMillan, M., Meister, R., Milne, G., Mouginot, J., Muir, A., Nicolas, J. P., Paden, J., Payne, A. J., Pritchard, H., Rignot, E., Rott, H., Sørensen, L. S., Scambos, T. A., Scheuchl, B., Schrama, E. J. O., Smith, B., Sundal, A. V., van Angelen, J. H., van de Berg, W. J., van den Broeke, M. R., Vaughan, D. G., Velicogna, I., Wahr, J., Whitehouse, P. L., Wingham, D. J., Yi, D., Young, D., and Zwally, H. J.: A Reconciled Estimate of Ice-Sheet Mass Balance, Science, 338, 1183–1189,, 2012. a

Shepherd, A., Ivins, E., Rignot, E., Smith, B., van den Broeke, M., Velicogna, I., Whitehouse, P., Briggs, K., Joughin, I., Krinner, G., Nowicki, S., Payne, T., Scambos, T., Schlegel, N., A, G., Agosta, C., Ahlstrøm, A., Babonis, G., Barletta, V., Blazquez, A., Bonin, J., Csatho, B., Cullather, R., Felikson, D., Fettweis, X., Forsberg, R., Gallee, H., Gardner, A., Gilbert, L., Groh, A., Gunter, B., Hanna, E., Harig, C., Helm, V., Horvath, A., Horwath, M., Khan, S., Kjeldsen, K. K., Konrad, H., Langen, P., Lecavalier, B., Loomis, B., Luthcke, S., McMillan, M., Melini, D., Mernild, S., Mohajerani, Y., Moore, P., Mouginot, J., Moyano, G., Muir, A., Nagler, T., Nield, G., Nilsson, J., Noel, B., Otosaka, I., Pattle, M. E., Peltier, W. R., Pie, N., Rietbroek, R., Rott, H., Sandberg-Sørensen, L., Sasgen, I., Save, H., Scheuchl, B., Schrama, E., Schröder, L., Seo, K.-W., Simonsen, S., Slater, T., Spada, G., Sutterley, T., Talpe, M., Tarasov, L., van de Berg, W. J., van der Wal, W., van Wessem, M., Vishwakarma, B. D., Wiese, D., Wouters, B., and the IMBIE team: Mass balance of the Antarctic Ice Sheet from 1992 to 2017, Nature, 558, 219–222,, 2018. a, b

Sime, L. C., Tindall, J. C., Wolff, E. W., Connolley, W. M., and Valdes, P. J.: Antarctic isotopic thermometer during a CO2 forced warming event, J. Geophys. Res.-Atmos., 113, D24119,, 2008. a, b, c, d

Sjolte, J., Sturm, C., Adolphi, F., Vinther, B. M., Werner, M., Lohmann, G., and Muscheler, R.: Solar and volcanic forcing of North Atlantic climate inferred from a process-based reconstruction, Clim. Past, 14, 1179–1194,, 2018. a, b

Steiger, N. J.: Historical climate model output of ECHAM5-wiso from 1871–2011 at T106 resolution, Data set, Zenodo,, 2018. a

Steiger, N. J., Hakim, G. J., Steig, E. J., Battisti, D. S., and Roe, G. H.: Assimilation of time-averaged pseudoproxies for climate reconstruction, J. Climate, 27, 426–441,, 2014. a

Steiger, N. J., Steig, E. J., Dee, S. G., Roe, G. H., and Hakim, G. J.: Climate reconstruction using data assimilation of water isotope ratios from ice cores, J. Geophys. Res.-Atmos., 122, 1545–1568,, 2017. a, b, c

Stenni, B., Curran, M. A. J., Abram, N. J., Orsi, A., Goursaud, S., Masson-Delmotte, V., Neukom, R., Goosse, H., Divine, D., van Ommen, T., Steig, E. J., Dixon, D. A., Thomas, E. R., Bertler, N. A. N., Isaksson, E., Ekaykin, A., Werner, M., and Frezzotti, M.: Antarctic climate variability on regional and continental scales over the last 2000 years, Clim. Past, 13, 1609–1634,, 2017a. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa, ab, ac, ad

Stenni, B., Curran, M. A. J., Abram, N. J., Orsi, A. J., Goursaud, S., Masson-Delmotte, V., Neukom, R., Goosse, H., Divine, D. V., van Ommen, T. D., Steig, E. J., Dixon, D. A., Thomas, E. R., Bertler, N. A. N., Isaksson, E., Ekaykin, A. A., Werner, M., and Frezzotti, M.: PAGES Antarctica2k Temperature Reconstructions, Data set, NOAA, available at: (last access: 18 March 2018), 2017b. a

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 Sys., 5, 146–172,, 2013. a

Taylor, K. E., Stouffer, R. J., and Meehl, G. A.: An overview of CMIP5 and the experiment design, B. Am. Meteorol. Soc., 93, 485–498,, 2012. a

Thomas, E. R., Hosking, J. S., Tuckwell, R. R., Warren, R. A., and Ludlow, E. C.: Twentieth century increase in snowfall in coastal West Antarctica, Geophys. Res. Lett., 42, 9387–9393,, 2015. a

Thomas, E. R., van Wessem, J. M., Roberts, J., Isaksson, E., Schlosser, E., Fudge, T. J., Vallelonga, P., Medley, B., Lenaerts, J., Bertler, N., van den Broeke, M. R., Dixon, D. A., Frezzotti, M., Stenni, B., Curran, M., and Ekaykin, A. A.: Regional Antarctic snow accumulation over the past 1000 years, Clim. Past, 13, 1491–1513,, 2017. a, b, c, d, e, f, g, h, i, j, k, l, m, n, o, p, q, r, s, t, u, v, w, x, y, z, aa

Tindall, J. C., Valdes, P. J., and Sime, L. C.: Stable water isotopes in HadCM3: Isotopic signature of El Niño-Southern Oscillation and the tropical amount effect, J. Geophys. Res.-Atmospheres, 114, D04111,, 2009. a, b

Trusel, L. D., Frey, K. E., Das, S. B., Karnauskas, K. B., Kuipers Munneke, P., Van Meijgaard, E., and Van Den Broeke, M. R.: Divergent trajectories of Antarctic surface melt under two twenty-first-century climate scenarios, Nat. Geosci., 8, 927–932,, 2015. a

Turner, J., Connolley, W. M., and Marshall, G. J.: The performance of the hadley centre climate model (HadCM3) in high southern latitudes, Int. J. Climatol., 112, 91–112,, 2006. a

Turner, J., Lu, H., White, I., King, J. C., Phillips, T., Hosking, J. S., Bracegirdle, T. J., Marshall, G. J., Mulvaney, R., and Deb, P.: Absence of 21st century warming on Antarctic Peninsula consistent with natural variability, Nature, 535, 411–415,, 2016. a, b

Unden, P., Rontu, L., Järvinen, H., Lynch, P., Calvo, J., Cats, G., Cuxart, J., Eerola, K., Fortelius, C., Garcia-Moya, J. A., Jones, C., Lenderlink, G., McDonald, A., McGrath, R., Navascues, B., Woetman Nielsen, N., Odegaard, V., Rodriguez, E., Rummukainen, M., Room, R., Sattler, K., Sass, B. H., Savijarvi, O., Wichers Schreur, B., Sigg, R., The, H., and Tijm, A.: HIRLAM-5 Scientific Documentation, Environ. Geol., 43, page238,, 2002. a

Van Den Broeke, M., Van De Berg, W. J., and Van Meijgaard, E.: Firn depth correction along the Antarctic grounding line, Antarct. Sci., 20, 513–517,, 2008. a

van Leeuwen, P. J.: Particle Filtering in Geophysical Systems, Mon. Weather Rev., 137, 4089–4114,, 2009. a

van Wessem, J. M., van de Berg, W. J., Noël, B. P. Y., van Meijgaard, E., Amory, C., Birnbaum, G., Jakobs, C. L., Krüger, K., Lenaerts, J. T. M., Lhermitte, S., Ligtenberg, S. R. M., Medley, B., Reijmer, C. H., van Tricht, K., Trusel, L. D., van Ulft, L. H., Wouters, B., Wuite, J., and van den Broeke, M. R.: Modelling the climate and surface mass balance of polar ice sheets using RACMO2 – Part 2: Antarctica (1979–2016), The Cryosphere, 12, 1479–1498,, 2018. a, b, c, d

Wang, Y., Ding, M., van Wessem, J. M., Schlosser, E., Altnau, S., van den Broeke, M. R., Lenaerts, J. T., Thomas, E. R., Isaksson, E., Wang, J., and Sun, W.: A comparison of Antarctic ice sheet surface mass balance from atmospheric climate models and in situ observations, J. Climate, 29, 5317–5337,, 2016. a

Widmann, M., Goosse, H., van der Schrier, G., Schnur, R., and Barkmeijer, J.: Using data assimilation to study extratropical Northern Hemisphere climate over the last millennium, Clim. Past, 6, 627–644,, 2010. a

Wu, T., Song, L., Li, W., Wang, Z., Zhang, H., Xin, X., Zhang, Y., Zhang, L., Wu, F., Y, L., Zhang, F., Shi, X., Chu, M., Zhang, J., Fang, Y., Wang, F., Lu, Y., Liu, X., Wei, M., Liu, Q., Zhou, W., Dong, M., Zhao, Q., Ji, J., Laurent, L., and Zhou, M.: An overview of BCC climate system model development and application for climate change studies, J. Meteorol. Res., 28, 34–56,, 2014.  a

Yukimoto, S., Adachi, Y., Hosaka, M., Sakami, T., Yoshimura, H., Hirabara, M., Tanaka, T. Y., Shindo, E., Tsujino, H., Deushi, M., Mizuta, R., Yabu, S., Obata, A., Nakano, H., Tsuyoshi, K., Ose, T., and Kitoh, A.: A New Global Climate Model of the Meteorological Research Institute: MRI-CGCM3 – Model Description and Basic Performance, J. Meteorol. Soc. Jpn., 90, 23–64,, 2012. a

Short summary
Large uncertainties remain in Antarctic surface temperature reconstructions over the last millennium. Here, the analysis of climate model outputs reveals that snow accumulation is a more relevant proxy for surface temperature reconstructions than δ18O. We use this finding in data assimilation experiments to compare to observed surface temperatures. We show that our continental temperature reconstruction outperforms reconstructions based on δ18O, especially for East Antarctica.