Description: This dataset contains the time series (1241-1970 CE) of the first two principal components (PCs) of reconstructed winter (DJF) Sea Level Pressure (SLP) in the North Atlantic region (20-70N, 90W-40E). The SLP is reconstructed by matching modeled d18O patterns from an isotope enabled climate model simulation to the spatio-temporal variability of the winter d18O from Greenland ice cores, and extracting the model SLP from an ensemble of the best matching years. PC1 of reconstructed SLP corresponds to the North Atlantic Oscillation, and PC2 of reconstructed SLP corresponds to the Eastern Atlantic pattern. Both time series are normalized and centred. The sign of the data series is as shown in Figure 2 of Sjolte et al. (2018) with the corresponding spatial patterns shown in Supplementary Figure S3. For example, positive values of PC1 means positive NAO, while positive values of PC2 means negative pressure anomalies south of Greenland.

Description: Seasonal climate reconstructions of sea level pressure and temperature for the North Atlantic region 1241-1970. Details on the reconstruction method and validation of the data can be found in Sjolte et al. (2018, 2020). The data include reconstructed 2m temperature (t2m) (data 01) and sea level pressure (slp) (data 02). The summer reconstruction (May-Oct) is done in two versions, one with, and one without tree-ring data, while there is only one version for winter (DJF). The variables ending on ''tr'' include tree-ring data.

Description: Changes in solar activity have previously been proposed to cause decadal- to millennial-scale fluctuations in both the modern and Holocene climates1. Direct observational records of solar activity, such as sunspot numbers, exist for only the past few hundred years, so solar variability for earlier periods is typically reconstructed from measurements of cosmogenic radionuclides such as 10Be and 14C from ice cores and tree rings2, 3. Here we present a high-resolution 10Be record from the ice core collected from central Greenland by the Greenland Ice Core Project (GRIP). The record spans from 22,500 to 10,000 years ago, and is based on new and compiled data4, 5, 6. Using 14C records7, 8 to control for climate-related influences on 10Be deposition, we reconstruct centennial changes in solar activity. We find that during the Last Glacial Maximum, solar minima correlate with more negative δ18O values of ice and are accompanied by increased snow accumulation and sea-salt input over central Greenland. We suggest that solar minima could have induced changes in the stratosphere that favour the development of high-pressure blocking systems located to the south of Greenland, as has been found in observations and model simulations for recent climate9, 10. We conclude that the mechanism behind solar forcing of regional climate change may have been similar under both modern and Last Glacial Maximum climate conditions.

Description: The incorporation of water isotopologues into the hydrology of general circulation models (GCMs) facilitates the comparison between modeled and measured proxy data in paleoclimate archives. However, the variability and drivers of measured and modeled water isotopologues, as well as the diversity of their representation in different models, are not well constrained. Improving our understanding of this variability in past and present climates will help to better constrain future climate change projections and decrease their range of uncertainty. Speleothems are a precisely datable terrestrial paleoclimate archives and provide well-preserved (semi-)continuous multivariate isotope time series in the lower latitudes and mid-latitudes and are therefore well suited to assess climate and isotope variability on decadal and longer timescales. However, the relationships of speleothem oxygen and carbon isotopes to climate variables are influenced by site-specific parameters, and their comparison to GCMs is not always straightforward. Here we compare speleothem oxygen and carbon isotopic signatures from the Speleothem Isotopes Synthesis and Analysis database version 2 (SISALv2) to the output of five different water-isotope-enabled GCMs (ECHAM5-wiso, GISS-E2-R, iCESM, iHadCM3, and isoGSM) over the last millennium (850–1850 CE). We systematically evaluate differences and commonalities between the standardized model simulation outputs. The goal is to distinguish climatic drivers of variability for modeled isotopes and compare them to those of measured isotopes. We find strong regional differences in the oxygen isotope signatures between models that can partly be attributed to differences in modeled surface temperature. At low latitudes, precipitation amount is the dominant driver for stable water isotope variability; however, at cave locations the agreement between modeled temperature variability is higher than for precipitation variability. While modeled isotopic signatures at cave locations exhibited extreme events coinciding with changes in volcanic and solar forcing, such fingerprints are not apparent in the speleothem isotopes. This may be attributed to the lower temporal resolution of speleothem records compared to the events that are to be detected. Using spectral analysis, we can show that all models underestimate decadal and longer variability compared to speleothems (albeit to varying extents). We found that no model excels in all analyzed comparisons, although some perform better than the others in either mean or variability. Therefore, we advise a multi-model approach whenever comparing proxy data to modeled data. Considering karst and cave internal processes, e.g., through isotope-enabled karst models, may alter the variability in speleothem isotopes and play an important role in determining the most appropriate model. By exploring new ways of analyzing the relationship between the oxygen and carbon isotopes, their variability, and co-variability across timescales, we provide methods that may serve as a baseline for future studies with different models using, e.g., different isotopes, different climate archives, or different time periods.

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

Description: Equatorial volcanic eruptions are known to impact the atmospheric circulation on seasonal time scales through a strengthening of the stratospheric zonal winds followed by dynamic ocean-atmosphere coupling. This emerges as the positive phase of the North Atlantic Oscillation in the first 5years after an eruption. In the North Atlantic, other modes of atmospheric circulation contribute to the climate variability but their response to volcanic eruptions has been less studied. We address this by retrieving the stable water isotopic fingerprint of the four major atmospheric circulation modes over the North Atlantic (Atlantic Ridge, Scandinavian Blocking and the negative and positive phases of the North Atlantic Oscillation (NAO�and NAOþ)) by using monthly precipitation data from Global Network of Isotopes in Precipitation (GNIP) and 500 mb geo-potential height from the 20th Century Reanalysis. The simulated stable isotopic pattern of each atmospheric circulation mode is further used to assess the retrieved pattern. We test if changes in the atmospheric circulation as well as moisture source conditions as a result of volcanic eruptions can be identified by analyzing the winter climate response after both equatorial and high-latitude North Hemispheric volcanic eruptions in data, reanalysis and simulations. We report of an NAOþmode in the first two years after equatorial eruptions followed by NAO � in year 3 due to a decrease in the meridional temperature gradient as a result of volcanic surface cooling. This emerges in both GNIP data as well as reanalysis. Although the detected response is stronger after equatorial eruptions compared to high latitude eruptions, our results show that the response after high latitude eruptions tend to emerge as NAO-� in year 2 followed by NAO+ in year 3–4.

Description: The volcanic fingerprint on the winter North Atlantic atmospheric circulation and climate is analyzed in six ensemble runs of ECHAM5/MPI-OM covering 800–2000 CE, both for equatorial and Northern Hemisphere (NH) eruptions. Large volcanic eruptions influence climate on both annual and decadal time scales due to dynamic interactions of different climate components in the Earth's system. It is well known that the North Atlantic Oscillation (NAO) tends to shift towards its positive phase during winter in the first 1–2 years after large tropical volcanic eruptions, causing warming over Europe, but other North Atlantic weather regimes have received less attention. Here we investigate the four dominant weather regimes in the North Atlantic: The negative and positive phase of NAO as well as the Atlantic Ridge, Scandinavian blocking. The volcanic fingerprint is detected as a change in the frequency of occurrence and anomalies in the wind and temperature fields as well as in the sea ice cover. We observe a strong significant increase in the frequency of Atlantic Ridge in the second year after equatorial eruptions that precede the NAO+ detected in year 3–5 as a result of a strong zonal wind anomalies in year 1–2. Evidence for a stronger polar vortex is detected in years 12–14 where NAO+ is detected both as a frequency increase and in the wind and temperature fields. A short-term response is also detected 2–4 years after NH eruptions. The longterm signal after NH eruptions indicate a weak polar vortex around a decade after an eruption. Although the signal after NH eruptions is weaker our results stress the need for further studies. The simulated atmospheric response recorded in ECHAM5 after volcanic eruptions suggest a more dynamic response than previously thought. The methodology used can also be applied to other forcing scenario, for example for future climate projections where the aim is to search for a long-term climate signal.

Description: The effect of external forcings on atmospheric circulation is debated. Due to the short observational period, the analysis of the role of external forcings is hampered, making it difficult to assess the sensitivity of atmospheric circulation to external forcings, as well as persistence of the effects. In observations, the average response to tropical volcanic eruptions is a positive North Atlantic Oscillation (NAO) during the following winter. However, past major tropical eruptions exceeding the magnitude of eruptions during the instrumental era could have had more lasting effects. Decadal NAO variability has been suggested to follow the 11-year solar cycle, and linkages have been made between grand solar minima and negative NAO. However, the solar link to NAO found by modeling studies is not unequivocally supported by reconstructions, and is not consistently present in observations for the 20th century. Here we present a reconstruction of atmospheric winter circulation for the North Atlantic region covering the period 1241–1970 CE. Based on seasonally resolved Greenland ice core records and a 1200-year-long simulation with an isotope-enabled climate model, we reconstruct sea level pressure and temperature by matching the spatiotemporal variability in the modeled isotopic composition to that of the ice cores. This method allows us to capture the primary (NAO) and secondary mode (Eastern Atlantic Pattern) of atmospheric circulation in the North Atlantic region, while, contrary to previous reconstructions, preserving the amplitude of observed year-to-year atmospheric variability. Our results show five winters of positive NAO on average following major tropical volcanic eruptions, which is more persistent than previously suggested. In response to decadal minima of solar activity we find a high-pressure anomaly over northern Europe, while a reinforced opposite response in pressure emerges with a 5-year time lag. On centennial timescales we observe a similar response of circulation as for the 5-year time-lagged response, with a high-pressure anomaly across North America and south of Greenland. This response to solar forcing is correlated to the second mode of atmospheric circulation, the Eastern Atlantic Pattern. The response could be due to an increase in blocking frequency, possibly linked to a weakening of the subpolar gyre. The long-term anomalies of temperature during solar minima shows cooling across Greenland, Iceland and western Europe, resembling the cooling pattern during the Little Ice Age (1450–1850 CE). While our results show significant correlation between solar forcing and the secondary circulation pattern on decadal (r = 0.29, p 〈 0.01) and centennial timescales (r = 0.6, p 〈 0.01), we find no consistent relationship between solar forcing and NAO. We conclude that solar and volcanic forcing impacts different modes of our reconstructed atmospheric circulation, which can aid in separating the regional effects of forcings and understanding the underlying mechanisms.

Description: The Antarctic temperature changes over the past millennia remain more uncertain than in many other continental regions. This has several origins: (1) the number of high-resolution ice cores is small, in particular on the East Antarctic plateau and in some coastal areas in East Antarctica; (2) the short and spatially sparse instrumental records limit the calibration period for reconstructions and the assessment of the methodologies; (3) the link between isotope records from ice cores and local climate is usually complex and dependent on the spatial scales and timescales investigated. Here, we use climate model results, pseudo-proxy experiments and data assimilation experiments to assess the potential for reconstructing the Antarctic temperature over the last 2 millennia based on a new database of stable oxygen isotopes in ice cores compiled in the frame- work of Antarctica2k (Stenni et al., 2017). The well-known covariance between δ18O and temperature is reproduced in the two isotope-enabled models used (ECHAM5/MPI-OM and ECHAM5-wiso), but is generally weak over the different Antarctic regions, limiting the skill of the reconstructions. Furthermore, the strength of the link displays large variations over the past millennium, further affecting the potential skill of temperature reconstructions based on statistical methods which rely on the assumption that the last decades are a good estimate for longer temperature reconstructions. Using a data assimilation technique allows, in theory, for changes in the δ18O–temperature link through time and space to be taken into account. Pseudoproxy experiments confirm the benefits of using data assimilation methods instead of statistical methods that provide reconstructions with unrealistic variances in some Antarctic subregions. They also confirm that the relatively weak link between both variables leads to a limited potential for reconstructing temperature based on δ18O. However, the reconstruction skill is higher and more uniform among reconstruction methods when the reconstruction target is the Antarctic as a whole rather than smaller Antarctic subregions. This consistency between the methods at the large scale is also observed when reconstructing temperature based on the real δ18O regional composites of Stenni et al. (2017). In this case, temperature reconstructions based on data assimilation confirm the long-term cooling over Antarctica during the last millennium, and the later onset of anthropogenic warming compared with the simulations without data assimilation, which is especially visible in West Antarctica. Data assimilation also allows for models and direct observations to be reconciled by reproducing the east–west contrast in the recent temperature trends. This recent warming pattern is likely mostly driven by internal variability given the large spread of individual Paleoclimate Modelling Intercomparison Project (PMIP)/Coupled Model Intercomparison Project (CMIP) model realizations in simulating it. As in the pseudoproxy framework, the reconstruction methods perform differently at the subregional scale, especially in terms of the variance of the time series produced. While the potential benefits of using a data assimilation method instead of a statistical method have been highlighted in a pseudoproxy framework, the instrumental series are too short to confirm this in a realistic setup.