
==== Front
Clim Dyn
Clim Dyn
Climate Dynamics
0930-7575
1432-0894
Springer Berlin Heidelberg Berlin/Heidelberg

7195
10.1007/s00382-024-07195-1
Article
Revising Alpine summer temperatures since 881 CE
http://orcid.org/0000-0002-1246-6030
Kuhl Eileen eikuhl@uni-mainz.de

1
http://orcid.org/0000-0003-3919-014X
Esper Jan 17
Schneider Lea 2
http://orcid.org/0000-0002-2683-8704
Trouet Valerie 3
http://orcid.org/0009-0001-7116-1114
Kunz Marcel 1
Klippel Lara 4
http://orcid.org/0000-0002-3821-0818
Büntgen Ulf 5789
http://orcid.org/0000-0001-9492-4674
Hartl Claudia 6
1 https://ror.org/023b0x485 grid.5802.f 0000 0001 1941 7111 Department of Geography, Johannes Gutenberg University, Johann-Joachim-Becher Weg 32, 55128 Mainz, Germany
2 https://ror.org/033eqas34 grid.8664.c 0000 0001 2165 8627 Department of Geography, Justus-Liebig-University, Gießen, Germany
3 https://ror.org/03m2x1q45 grid.134563.6 0000 0001 2168 186X Laboratory of Tree-Ring Research, University of Arizona, Tucson, USA
4 https://ror.org/02nrqs528 grid.38275.3b 0000 0001 2321 7956 Deutscher Wetterdienst, Offenbach, Germany
5 https://ror.org/013meh722 grid.5335.0 0000 0001 2188 5934 Department of Geography, University of Cambridge, Cambridge, UK
6 Nature Rings - Environmental Research and Education, Mainz, Germany
7 Global Change Research Centre (CzechGlobe), Brno, Czech Republic
8 grid.419754.a 0000 0001 2259 5533 Swiss Federal Research Institute (WSL), Birmensdorf, Switzerland
9 https://ror.org/02j46qs45 grid.10267.32 0000 0001 2194 0956 Department of Geography, Masaryk University, Brno, Czech Republic
22 3 2024
22 3 2024
2024
62 7 61416157
30 11 2023
10 3 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by/4.0/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Europe experienced severe heat waves during the last decade, which impacted ecological and societal systems and are likely to increase under projected global warming. A better understanding of pre-industrial warm-season changes is needed to contextualize these recent trends and extremes. Here, we introduce a network of 352 living and relict larch trees (Larix decidua Mill.) from the Matter and Simplon valleys in the Swiss Alps to develop a maximum latewood density (MXD) chronology calibrating at r = 0.8 (p > 0.05, 1901–2017 CE) against May–August temperatures over Western Europe. Machine learning is applied to identify historical wood samples aligning with growth characteristics of sites from elevations above 1900 m asl to extend the modern part of the chronology back to 881 CE. The new Alpine record reveals warmer conditions in the tenth century, followed by an extended cold period during the late Medieval times, a less-pronounced Little Ice Age culminating in the 1810s, and prolonged anthropogenic warming until present. The Samalas eruption likely triggered the coldest reconstructed summer in Western Europe in 1258 CE (-2.32 °C), which is in line with a recently published MXD-based reconstruction from the Spanish Pyrenees. Whereas the new Alpine reconstruction is potentially constrained in the lowest frequency, centennial timescale domain, it overcomes variance biases in existing state-of-the-art reconstructions and sets a new standard in site-control of historical samples and calibration/ verification statistics.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00382-024-07195-1.

Keywords

Climate reconstruction
Climate change
Dendrochronology
Europe
European Alps
Tree rings
http://dx.doi.org/10.13039/501100001659 Deutsche Forschungsgemeinschaft HA 8048/1-1 Hartl Claudia http://dx.doi.org/10.13039/501100000781 European Research Council AdG 882727 AdG 882727 AdG 882727 Esper Jan Kunz Marcel Büntgen Ulf http://dx.doi.org/10.13039/501100001823 Ministerstvo Školství, Mládeže a Tělovýchovy CZ.02.01.01/00/22_008/0004635 CZ.02.01.01/00/22_008/0004635 Esper Jan Büntgen Ulf Johannes Gutenberg-Universität Mainz (1030)Open Access funding enabled and organized by Projekt DEAL.

issue-copyright-statement© Springer-Verlag GmbH Germany, part of Springer Nature 2024
==== Body
pmcIntroduction

There is a rising urgency to improve our understanding of natural and anthropogenic drivers of temperature variations and reduce uncertainties in models predictions (Eyring et al. 2019). Analysing past climate variability, identifying spatial patterns and updating paleoclimate models is indispensable to establish the basic conditions for these models (Esper and Büntgen 2021). Common paleoclimate archives are tree-rings from high-elevation and high-latitude sites, which are used to reconstruct past temperature variability and place anthropogenic warming into historical context (Schweingruber 1988). This is not only true for long-term changes in climate history but for annual fluctuations, for example, exceptionally warm/cold and wet/dry years. With information on extreme events, we can set recent extreme events (and the likelihood of future events) in a historical context regarding their timing, occurrence and strength (e.g., Qin et al. 2011; Lyu et al. 2016; Borkotoky et al. 2021).

The European Alps have been of interest for such studies for decades as the regional warming already exceeded 2 °C from 1864 to 2017 CE (Allgaier Leuch et al. 2017). Tree-ring width (TRW) and maximum latewood density (MXD) have been used as proxies to develop eleven Alpine temperature reconstructions (Schweingruber et al. 1987, 1988; Büntgen et al. 2005, 2006; Frank et al. 2005; Frank and Esper 2005a; Esper et al. 2007b; Corona et al. 2010, 2011; Trachsel et al. 2012; Coppola et al. 2013; Leonelli et al. 2016; Table 1). Mainly, these publications used larch (Larix decidua Mill.) tree-ring series with one exception, which is a spruce (Picea abies L.) reconstruction (Esper et al. 2007b). Table 1 Existing TRW- or MXD-based temperature reconstructions based on from the European Alps in order of publication year. The table includes the number of included sites, the elevation, the used proxies, the range of years covered, as well as the species (Larix decidua Mill. = LADE, Abies alba Mill. = ABAL, Pinus cembra L. = PICE, Picea abies L. = PIAB) and the temperature signal used for reconstructing (letters denote to the first letter of the corresponding months, e.g. JJA June, July, August)

No	Publication	# of Sites	Elevation [m asl]	Proxy	Period [CE]	Species	Temperature Signal	
1	Schweingruber et al. 1987, 1988h	16	100–1700, 1800	MXD	982–1976	PIAB, ABAL	JJAS	
2	Büntgen et al. 2005h	5	 > 1500	TRW	951–2002	LADE, PICE	JJA	
3	Frank and Esper 2005a, b	53	 > 1500	TRW / MXD	1600–1988/ 1650–1987	PIAB, ABAL, LADE, PICE	JJA / AMJJAS	
4	Frank et al. 2005	53	 > 1500	TRW	1760–1990	PIAB, ABAL, LADE, PICE	JJA	
5	Büntgen et al. 2006h	4	 > 1900	MXD	755–2004	LADE	JJAS	
6	Esper et al. 2007bh	4	-	TRW /MXD	1028–2003	PIAB	JJ / AS	
7	Corona et al. 2010h	38	 > 1725	TRW & MXD	1000–2000	LADE, PICE	JJA	
8	Corona et al. 2011h	34	 > 1600	TRW	751–2003	LADE	JJA	
9	Coppola et al. 2013	4	 > 1910	TRW	1610–2008	LADE	JJA	
10	Trachsel et al. 2012h	4	-	TRW, MXD, BI, Chrionomid/ biogenic Silica	755–2004	LADE, PIAB, PICE	JJA	
11	Leonelli et al. 2016	42	 > 1800	TRW	1470–2010	LADE, PIAB, PICE	JJA	
h = used historical series from buildings

Most analyses targeted summer temperatures from June to August (JJA), whereas a few MXD chronologies were used for longer seasons: Schweingruber et al. (1987, 1988) and Büntgen et al. (2006) calibrated their MXD chronologies to June to September (JJAS) temperatures, while Frank and Esper (2005a, b) reconstructed April to September temperatures. One exception is Esper et al. (2007b), which found best correlations with temperatures from June-July for TRW and August–September for MXD. Seven of the eleven reconstructions extend back to the first millennium CE (Schweingruber et al. 1987, 1988; Büntgen et al. 2005, 2006; Esper et al. 2007b; Corona et al. 2010, 2011; Trachsel et al. 2012). MXD-based chronologies reached correlations up to 0.69 with instrumental JJAS temperatures (Büntgen et al. 2006). The most recent reconstruction reported correlations up to 0.78 with JJA from 1763 to 2014 CE (Leonelli et al. 2016). Out of all mentioned publications, Büntgen et al. (2006) received by far the most citations for their reconstruction from the Swiss Lötschental (hereafter referred to as Lötschental) implying the data often being used in other studies and the findings having a generally high impact on the scientific community. Above all, it is the longest existing MXD-based reconstruction from the Alps based on a single species and serves as comparison for this study.

The mentioned millennium-long chronologies all include wood samples from historical buildings (e.g. Corona et al. 2010; Trachsel et al. 2012) as high-elevation forests in the Swiss Alps have been heavily influenced by anthropogenic forest use since the Middle Ages (Conedera et al. 2017). Therefore, the likelihood of living trees spanning earlier ages is low and in situ dead wood in the Alps is rare. The long history of human activity, however, enables the use of construction timber of old, high-elevation settlements like Zermatt in Switzerland (Coolidge 1912). The recurring issue with historical tree-ring samples in all of the existing reconstructions spanning a millennium is the limited site control and validation that these woods originate from elevations near treeline and contain strong temperature signals (Riechelmann et al. 2020). Multiple studies showed that the temperature signal in tree-ring parameters weakens with decreasing elevation (Neuwirth et al. 2004; Affolter et al. 2010; King et al. 2013; Hartl-Meier et al. 2014a, b, 2017b; Salzer et al. 2014; Zhang et al. 2015; Hartl et al. 2022), which can consequently influence the temporal robustness of a reconstruction. Kuhl et al. (2023) recently introduced a method to mitigate these biasing effects by fitting tree-ring parameters and growth characteristics to a classification model to sort historical series to living stands with different distances to treeline. This provenance method was tested on data from the Swiss Simplon Valley in Kuhl et al. (2023) and is now applied to the Swiss Matter Valley data as well to exclude historical series with growth behaviours similar to sites further away from the current treeline. With this approach, we aim to build an improved millennium-length temperature reconstruction by the selection of historical material based on the information we gain from the provenance models.

We present a new MXD-based temperature reconstruction from 352 living and historic larch trees (Larix decidua Mill.) sampled in the southwestern Swiss Alps. We analyse the influences of detrending and variance stabilisation on reconstruction variability and trends, assess and mitigate the effects of larch budmoth infestation events, and revise Alpine temperature history since 881 CE including severe changes compared to existing reconstructions from the European Alps. The record is compared with other high-resolution warm season temperature reconstructions, and common trends and extremes are discussed to evaluate large-scale climate forcings.

Methods

Study sites and data distribution

The study sites in the canton Valais in the South-Western Swiss Alps include a total number of 352 series from living trees and historical buildings from valleys south of the Rhone valley (Fig. 1a-c). In the Simplon (SV) and Matter (MV) valleys (Hartl et al. 2022; Kunz et al. 2023), 146 living series of larch (Larix decidua Mill.) from high elevations (hereafter SV1-3 and MV1-3) were sampled. Likewise, 99 historical series from buildings in Simplon village and 206 historical series from buildings in Zermatt and Zmutt villages (Schmidhalter, M., Riechelmann et al. 2013, 2020) were sampled (Fig. 1b-d). From these samples, X-ray densitometric measurements were conducted using a Walesch2003 (WALESCH, Electronic GmbH, Switzerland) as described in Björklund et al. (2019). The distribution of historical series in time was evaluated to prevent biases arising from an unbalanced dataset (Fig. 2, Esper et al. 2016).Fig. 1 a Study site in the Swiss Alps close to the Italian border, b and locations of the sampling sites for living trees (turquoise triangles) and historical buildings (dark blue quarters) in the Matter Valley (MV) and the Simplon Valley (SV) c The historical village of Zmutt in the Matter Valley (credits: C. Hartl) d Information about the elevation [m asl] and the number of samples (increment cores and discs) per site

Fig. 2 Temporally balanced distribution of the 146 living (grey) and 206 relict (coloured) samples, sorted by their innermost rings in a bar plot. Pie charts show the even distribution of samples per house (top) and the total number of yearly measurement points per building (bottom). Both, pie charts and bar plot support a balanced ratio between the number of historical and living measurements

Provenancing of the historical material

To determine the provenance of the historical series, two Machine Learning (ML) models were trained based on densitometric measurements and statistical data from transect sites in both valleys. For SV, the procedure is described in Kuhl et al. (2023) and was likewise applied to MV. Compared to SV, MV had only a limited number of available sites from different elevations (MV1 and MV2/MV3). Hence, samples showing a prediction probability lower than 0.8 in the ML model output were excluded from the analysis. In rare cases, when A and B samples were predicted to origin from different elevation classes, these were manually assigned to the most likely elevation. From all measured historical samples, 36 (out of 99) from SV and likewise 170 (out of 206) from MV were used for this study. The final datasets were produced by including the historical series to the corresponding living series collectives (Fig. 1c, Table S1).

Detrending

After provenancing, the sites were mean-adjusted valley by valley to the MXD level of the highest sites (Fig. S1). Differences in mean levels were estimated by calculating the offset between regional curves (RCs) over a period of replication ≥ 15 series (100 cambial years SV, 250 cambial years MV) for each site. The data from the two valleys were then combined and likewise adjusted in their mean levels over the first 300 years of cambial age meeting > 15 series. Missing rings were in-filled using Arstan (Version 41d, Cook 1985), series were power transformed, pruned (Briffa et al. 2001) beyond 300 years of cambial age (considering pith offsets) and standardized in R 4.2 (R Core Team 2021). To remove age-related density changes (Bräker 1981) the series were detrended as residuals (Cook and Peters 1997) using regional curve standardization (RCS, Esper et al. 2003) with a 67% mean cambial age smoothing filter (frequency response = 0.5). In addition, the signal free (SF) approach (Melvin and Briffa 2008) on RCS and age-dependent spline detrending (Melvin et al. 2007),as well as a 300-year smoothing spline (Cook and Peters 1981) and Hugershoff detrending (Cook et al. 1990) were applied for comparison (see Fig. S2). SF detrending was calculated using the program Signal Free (Version 45_v2b, Cook et al. 2017).

Larch budmoth treatment

Larch trees in these elevations and regions of the European Alps are affected by recurring larch budmoth (Zeiraphera griseana Hb., LBM) mass outbreaks, leading to defoliation, and consequently to reduced MXD values and radial growth in the event year (Esper et al. 2007a; Baltensweiler et al. 2008; Hartl-Meier et al. 2016, 2017b; Kunz et al. 2023). These distinct declines in MXD are independent of weather conditions and need to be removed when reconstructing temperature. Several methods were introduced by Büntgen et al. (2009) or Kunz et al. (2023) for LBM outbreak detection. From these methods, impulse indicator saturation (Pretis et al. 2018, IIS) is the most effective one for our purpose as it not only detects LBM events but also gives a correction factor for these years. The raw MXD series were detrended prior to IIS using a 30-year smoothing spline. This algorithm identifies rapid breaks (Fig. S3a) as negative outliers and uses these as correction coefficients to be subtracted from the RCS detrended series (Fig. S3b). The method was initially introduced to accurately detect volcanic events in simulated temperature time series and has proven its potential for LBM detection (Pretis et al. 2017; Schneider et al. 2017; Kunz et al. 2023). When applying IIS, a non-host chronology can be included as a regressor to ensure that the algorithm does not detect outliers that show climate induced declines (e.g. volcanic events) (Kunz et al. 2023). Between 1616 and 2017 CE a non-host Swiss stone pine (Pinus cembra L.) site from MV (Kunz et al. 2023) and a non-host larch (Larix decidua Mill.) site from the Northern Alps (Hartl-Meier et al. 2014a) were included in the algorithm (Fig. S3c). As the non-host chronologies only extend back to 1616 CE, prior comparison to the host chronology reaching back to 881 CE was prevented. Consequently, the algorithm detected known volcanic events as LBM events in this period. To avoid the correction of climatic induced growth declines, stratospheric aerosol optical depth (SAOD) years from Sigl et al. (2021) were averaged over the Westerlies zone (~ 30–60°N). In years, where SAOD values exceeded 0.03 and matched an IIS detected year, the LBM correction value was not subtracted from the host chronology. The valleys were treated individually to address the lag of outbreak years between the valleys (Johnson et al. 2004; Saulnier et al. 2017; Kunz et al. 2023). Comparison of LBM detections of the IIS algorithm corrected chronology and the non-host chronology (Fig. S3c) confirms that the algorithm can detect non-climate-related growth declines. The prevention of false detection of volcanic cooling events as LBM outliers by using the SAOD values is necessary to preserve climatic signals in the chronology, but it can lead to overestimated cooling when LBM events and volcanic induced cooling years align.

Variance stabilization

The final chronology was constructed using Tukey’s bi-weight robust mean. To account for changes in variance, chronology variance was stabilized using the method detailed in Osborn et al. (1997). On trial, different window lengths were tested during stabilization. Results showed shorter windows (e.g. 31 years) stabilize variance better than larger ones. However, applying a narrow window for variance stabilization increases the risk to eradicate low-frequency trends (Frank et al. 2007). We approached this risk by calculating residuals between the chronology and a 100-year low-pass Butterworth filter. The variance of the high-frequency residuals was then stabilized in a 31-year window, whereas the low-pass chronology was stabilized in a 201-year window (Osborn et al. 1997)). Afterwards, the stabilized sub-chronologies were added back together to the final chronology from 881 to 2017 CE. Calculations were performed in R 4.2 (R Core Team 2021) using the package “dplR” (Bunn 2010).

Climate signals

Pearson’s correlations between instrumental temperatures (TS 4.06; Harris et al. 2020) and the MXD chronology were calculated from 1901–2017 CE in a classical bootstrap approach (bootn = 1000, α = 0.95) using the R package “treeclim” (Zang and Biondi 2015). Correlations were computed for monthly January–September (J-S) and JJA and May–September (MJJAS) seasonal means. For every grid from -10°E to 20°W and 35–65°N, we likewise calculated correlations between MJJAS mean temperature anomalies and the tree-ring chronology. We averaged the temperature data of all grids exceeding r = 0.7. Static and running correlations (31-year window) were performed using this averaged temperature record.

Calibration, verification, and reconstruction

Three different reconstruction approaches (linear regression, polynomial regression, and scaling) were tested using a two-fold cross-validation from 1901–1959 and 1960–2017 CE against MJJAS grid-averaged temperature anomalies. The correlation between predicted and observed data (r), explained variance (R2) and root mean squared error (RMSE) were used to assess model accuracies. For the final reconstruction, the best performing model was chosen and then calibrated over the full period from 1901–2017 CE.

Uncertainty bands for the reconstruction were estimated by averaging RMSEs of sub-chronology reconstructions for different, stepwise increasing sample depths from n = 5–50 series in five steps. For each of these replication steps, bootstrapped sub-chronologies of n randomly sampled series were built (bootn = 1000). These individual sub-chronologies were then used to reconstruct MJJAS temperature anomalies on the calibration period 1914–2002 CE. For each reconstruction r, R2 and RMSE were calculated. The mean RMSE for each replication step (RMSEn) was then used to estimate the error of the reconstruction depending on its sample depth n (error = 2*RMSEn).

Standard reconstruction approach

A reconstruction including all existing historical series was additionally modelled neglecting provenance information. The data were gap-filled, power transformed and pruned (300-years). In contrast to the new approach, no mean adjustment was performed. All series were detrended, LBM corrected, and variance stabilized, just as the main chronology of this paper. A linear model with the same temperature data was used to reconstruct the MJJAS mean temperature variability (averaged grids), and the results compared with the novel provenance considering reconstruction.

Superposed epoch analysis and extreme year analysis

Superposed epoch analysis (Lough and Fritts 1987, SEA) was applied to evaluate the significance of major volcanic eruptions (Table S2) considering a period of five years prior (reference period) and five years after the event using the package “dplR” (Bunn 2010). For this, a 50-year high-pass filter was run on the new reconstruction and the Lötschental record (Büntgen et al. 2006). From bootstrap resampling (bootn = 10,000), 99% confidence intervals were calculated to estimate significant deviations after volcanic event year period. In addition, the coldest 20 years of the 50-year high-pass filtered records were compared against volcanic event years.

Results

Chronology development

The developed chronology covers the period from 881–2017 CE including 1137 years and thus has the potential to revise climate variability back to the early medieval period. The minimum sample depth n = 5 spans the period 881 to 901 CE and the maximum is reached in 1963 CE with n = 93 series. The average correlation between the series is 0.54, with a mean segment length of 146 years due to pruning of the data and a lag-1 autocorrelation of 0.3.

The removal of non-climate related biases in tree-growth can help to increase the robustness of a chronology. This improvement can be accomplished by selecting historical series, which have been sorted to high-elevation microclimate growth characteristics by the provenance model and are therefore more likely to have a stronger temperature signal. Considering this information of the series allows for exploration of differences in MXD mean level offsets between the sites (Fig. S1). If not addressed, these differences can lead to deflated or inflated values when RCS detrending is applied (Zhang et al. 2015; Römer et al. 2021; Hartl et al. 2022), especially if a period is mainly covered with series from one elevational growth characteristics. RCS mainly requires a great sample depth (here: 352 series), a homogeneous distribution of tree-ring series in time (Fig. 2), heterogenous tree age and homogeneity in the provenance of historical and living series (Briffa et al. 1992; Düthorn et al. 2015; Esper et al. 2003; Melvin et al. 2013), which is improved by the provenance approach of Kuhl et al. (2023). The RCS detrended chronology shows more variability and a stronger negative trend towards the early nineteenth century compared to 300-year spline and Hugershoff detrending (Fig. S2a), indicating a stronger preservation of the low-frequency (Briffa et al. 1992). The standard deviation of the 100-year smoothed chronology is almost twice the size compared to other detrending methods when we apply SF age-dependent spline detrending (Fig. S2a). The SF RCS chronology, however, does not show significant differences to the classical RCS approach (Fig. S2b).

Spectrum analysis emphasizes that a smaller window of 31 years can stabilize the variance better than a wider window of 201 years. However, stabilization on narrow window can remove low-frequency signals which are related to temperature (Fig. S4a). While the 31-year window stabilization has a reduced power spectrum moving towards lower frequencies, greater windows (e.g. 201 years) preserve more low-frequency information but in turn do not stabilize the variance as well (Fig. S4b). Our new approach detangles low- from high-frequency bands by stabilizing the variance individually on both, the high- and the low-frequency sub-chronologies. Treating low- and high-frequency on different window lengths reduces the loss of low-frequency signals when applying a short window but at the same time improves stabilization compared to using a wider window for the full chronology.

Climate-growth relationship

Correlations between temperature anomalies and the MXD chronology reveal a positive relationship between tree-growth and temperature. Spatial correlations for MJJAS and the chronology between 1901–2017 CE show the significantly high (p < 0.05) response between Western European temperature variability and MXD indices (Fig. 3a). Grid points over France, Switzerland, Southern Germany, Northern Italy, and Western Austria reached highest correlations of r ≥ 0.7. Altogether, correlations do not drop below 0.5 in Western Europe. The ability of the Alpine chronology to represent past temperature variability in Western Europe is strengthened by calculating static correlations from 1901–2017 CE between temperatures of the grids exceeding 0.7 for MJJAS in Fig. 3a and the MXD indices (Fig. 3b). From March onwards, significant correlations (p < 0.05) with temperature are observed with highest monthly correlations for August (r = 0.71). The coherence between the proxy and instrumental temperatures increases when considering seasonal means up to r = 0.8 (MJJAS). Moving window correlations indicate a temporally stable relationship over the full instrumental period from 1901–2017 CE (Fig. 3c). Although correlations slightly decline in recent years, r values do not fall below 0.6.Fig. 3 a Correlations between the final (combined and variance stabilized) MXD chronology and gridded mean May–September (MJJAS) temperature data (CRU TS 4.06 0.5°) from 1901–2017 CE b Temperature data of all grids with spatial MJJAS correlations ≥ 0.7 (a) are averaged to calculate monthly correlations between January and September, as well as using two different summer means (June–August JJA and May–September MJJAS). Grey colouring denotes non-significant correlations (p < 0.05), orange shows significant correlations and red illustrates the MJJAS window that is chosen for further analysis c Moving 31-year correlations between MJJAS mean temperature and the MXD chronology

Reconstruction of Alpine summer temperature anomalies

For reconstruction, a linear regression model was chosen as these models perform best on both cross-validation periods (1901–1959 CE, 1960–2017 CE) with correlations of 0.8 on the entire period 1901–2017 CE between observed and predicted time series (Fig. 4a). These cross-validation models explain up to 70% of the variance and are stable in time with the validation periods reaching r > 0.82 (Table S2). Although both models exhibit a high level of performance, they predict underexaggerated values for the latest decades. While the model trained on the early period fails to capture the recent warming from 1995 CE onwards, the model for the later period cannot accurately predict the increase in temperatures since 2010 CE.Fig. 4 a Linear Regression Models between the MXD chronology and May–September mean temperature anomalies [w.r.t. 1961–1990 CE] (MJJAS Temp) of calibration periods (cp.) 1901–1959 CE (blue) and 1960–2017 CE (red). Detailed performance measures of the cross-validation are listed in Table S2 b The effect of replication on regression model accuracy: Correlations r between MJJAS Temp and instrumental data between 1914–2002 CE depending on different sample replications (5–50). Black dots represent correlation coefficient between the individual model runs (1000 per replication step) and instrumental data. Red dots represent the medians of the 1000 runs (rmedian). The black line shows r for the regression performed on the full period 1901–2017 (min. replication 71 samples) c Model fit of the linear regression between the MXD indexes and instrumental data on the full period 1901–2017 CE, which is used for the reconstruction d Final MJJAS Temp reconstruction: the blue line shows the reconstruction with a 100-yr smoothing filter (dark blue). Grey shading presents the error range of the reconstruction and orange the instrumental data. Light blue denotes to the sample depth n

Uncertainties of the reconstruction are highly related to changes in replication. With decreased replication further back in time, model uncertainties are likely increasing. To test the reliability of the linear regression model depending on replication, randomly built sub-chronologies of different replications were fitted in linear regression models over a period from 1914–2002 CE, where the replication is constantly over 50 series. This smaller window was chosen to test the replication effect from 5 up to 50 series per sub-chronology. Results show a high correlation of the sub-chronologies to the instrumental data until the replication drops below 15 series (rmedian ≥ 0.8) (Fig. 4b). When the sample replication equals ten series or lower, the linear models lose accuracy as the confidence intervals increase towards lower correlation values. Still, rmedian does not fall below 0.7.

The robust temperature signals and the stability of the linear model with respect to replication fluctuations enable the construction of a millennium-length MJJAS temperature reconstruction. The reconstruction model fitted to 1901–2017 CE explains 64% of the variance with r = 0.8 and is used to predict and revise Alpine temperature anomalies until 881 CE (Fig. 4c, d).

The temperature reconstruction (Fig. 4d) from the European Alps displays that long-term temperature trends decrease during the early and high medieval times from 881–1200 CE. Temperatures rise again until late medieval ages in the co-called medieval climate anomaly (MCA) and peak in the 1360 s, which are on average 0.76 °C (uncertainty range: -0.08 to 1.6 °C) warmer than the reference period 1961–1990 CE. This maximum is followed by a slow long-term decrease in temperatures during the Little Ice Age (LIA) until the early nineteenth century, when the Tambora eruption caused the “year without summer” of 1816 CE (Stothers 1984; Oppenheimer 2003). This decade is the coldest decade of the record and marks the peak of the LIA. Here, temperature anomalies drop to -1.28 °C (uncertainty range -2.12 to -0.44 °C). Between the end of the MCA and the peak of the LIA, temperatures decrease in total by 2.04 °C. Afterwards temperatures rise again from 1900 until 2017 CE. This increase is shortly disrupted between 1950 and 1980 CE.

The most prominent long-term periods of cooling are between 1000–1300 and 1750–1870 CE, while warmer phases are found between 910–1000, 1300–1400, 1475–1575 and 1900–2017 CE. The 971–980 CE period represents the warmest decade of the record with a decadal average of + 1.17 °C (uncertainty range: 0.25 to 2.09 °C), followed by the 960 s, 1500 s, 980 s and 2000s. Besides the 1810s, coldest decades of the record are the 1180 s, 1040 s, 1030 s and 1170 s. The absolute warmest and coldest years of the record are not necessarily located within these decades. Warmest years are recorded in 970 (2.72 °C), 1483 (2.07 °C), 2003 (2.06 °C), 1100 (1.82 °C) and 980 CE (1.80 °C), while the coldest years are 1180 CE (-2.9 °C), 1258 (-2.86 °C), 1181 (-2.8 °C), 1816 (-2.43 °C) and 1821 CE (-2.29 °C).

Discussion

The new approach: eliminating interference signals from elevation and insects

The here presented new approaches aim to make reconstructions, which depend on historical series, more robust. Elevation and the potential weak temperature signals in historical series have been neglected in the past when building chronologies. Using a ML based provenance model helps to overcome this problem (Kuhl et al. 2023). The provenance of the historical series allows the selection of series, that are likely to be most sensitive to temperature, and therefore improves the temporal stability of the climate sensitivity in periods beyond the lifespans of the living trees.

As MXD increases with decreasing elevations (Zhang et al. 2015; Hartl et al. 2022), mean adjusting the series can reduce elevational biases while detrending with RCS. Lower elevation series have higher MXD levels on cambial age scale compared to higher elevations and can therefore potentially alter the regional curve towards higher mean values. This can be true for elevational differences of 100 m (Fig. S1). Microsite differences like this have already been observed and critically viewed in multiple studies in terms of regional curves biases in RCS detrending (Gunnarson et al. 2011; Düthorn et al. 2013; Hartl et al. 2021, 2022). In previous reconstructions, these potential offsets between historical series have been neglected, which is mainly due to the lack of information on the elevational origin of historical tree-ring series. Simultaneously, differing MXD values of the series can also result from temperature variability and changes in treeline (Büntgen et al. 2022a) and there is a likelihood that mean adjustment causes a loss of temperature information. However, site control of historical series offers the possibility to improve RCS detrending of living and historical series to reduce visible biases in the data due to elevational non-temperature related MXD offsets as presented in Fig. S1.

The correction of negative outliers is indispensable when chronologies are based on tree species with recurring insect outbreaks like larch (Büntgen et al. 2006; Aryal et al. 2023). As the outliers in a chronology can potentially result from volcanic cooling, the inclusion of a non-host chronology in combination with SAOD-analysis present necessary pre-steps before applying correction coefficients. This can ensure that high-frequency temperature signals are not removed from the data during this process.

The chronology depicts a different course of the temperature history when provenance information of historical series is neglected in the reconstruction process (Fig. S5). This reconstruction profits from more available samples to extend the record further until 729 (n = 5). Compared to the main reconstruction, this record shows higher temperature estimates in the medieval period until 1000 CE and between the twelfth-fifteenth century. A visible decrease in temperatures during the LIA is almost 150 years later than in the main reconstruction but shows a stronger decline due to the increased values from the twelfth-fifteenth century. Higher temperature estimates in this period most likely result from the inclusion of historical series from lower elevation growth characteristics with higher MXD levels and from the neglection of potential level offsets (e.g., see Hartl et al. 2022 for SV).

Mitigated low frequency signals

For almost two decades, the Lötschental reconstruction has been the state of the art record for European summer temperatures and was used in several larger network studies to present Western European temperature variability (Esper et al. 2007c; Brázdil et al. 2010; Luterbacher et al. 2016; Wilson et al. 2016; Anchukaitis et al. 2017; Torbenson et al. 2023; Wang et al. 2023). The comparison between the new reconstruction and the Lötschental record highlights some major differences between 1000–1900 CE (Fig. 5a). In contrast to the Lötschental record, the new reconstruction shows less long-term variability. Although a long-term temperature decrease during the LIA is observed, it is not as pronounced as in the Lötschental. There, the onset of the cold period begins after a peak at 1200 CE whereas the cooling in the new reconstruction starts approximately 160 years later in the mid-fourteenth century. In the Lötschental record, the drop in temperatures (w.r.t. 1901–2000 CE) between the warmest pre-LIA decade in the 1200 s (0.66 °C) and the 1810s (-2.31 °C) is 2.97 °C and the increase in temperatures between the coldest decade (1810s) and the warmest following decade of the 1940s is 3.08 °C. The reconstructed temperature decrease since the end of the MCA is almost a degree less in the new record and could hint that the European Little Ice Age temperature decrease was not as strong as the Lötschental record estimated and that the reconstruction implies that summer temperatures are not as cool as in the period between 1000 and 1250 CE. However, there is evidence that the LIA must have existed in Europe, for example by historical records or glacier advances in the Swiss Alps (Grove 2001), which challenges the reconstruction’s skill in low-frequency variation.Fig. 5 a Comparison between the Lötschental (green) and our record (blue) and b comparison between the Pyrenees (coral) and our record (blue). Dark lines in a and b present the 30-year smoothed data. r equals the correlation between both records in their common period, respectively c Differences in standard deviation of the Pyrenees (coral), the Lötschental (green) and the new Alpine (blue) records (100-year running window)

We tested and eliminated various potential sources of error during the development of the reconstruction (see Table S3 for more detailed information). Sources, which can alter the low-frequency signals, include the distribution of tree age in time, the detrending method, the selection of series or an inadequate classification by the ML models (Kuhl et al. 2023). For example, the model for MV uses data of only two different elevations as reference and despite the model performance score (f1-score) of 0.99, series might align more with patterns from lower elevations but will be sorted to either of the two elevations. We tackle this problem by selecting only those series that are sorted to an elevational class with a prediction probability greater than 0.8. As mentioned before, the mean adjustment might alter climate signals, when tree-ring series are adjusted to higher elevation MXD levels. However, excluding these steps (provenancing and mean adjustment) does not increase the low-frequency signal, especially not during the LIA (Fig. S5).

Besides the difference in the low-frequency signal in both records, the variance of the new record is more stable in time (Fig. 5c). In the Lötschental record, however, the variance is not as stable and shows stronger temporal fluctuations, especially before the thirteenth century CE. As variance changes are likely artificial and an artefact of alternating replication or varying ecological habitats, the increase in variance pre 1200 CE might result from a change in replication and the different origin of this historical tail being exclusively composed of the Simplon valley (Büntgen et al. 2006; Esper et al. 2016). Despite these differences, the correlation in the common period (881–2004 CE) between both records is 0.54 and 0.62, when a 500-year high-pass filter is applied to the Lötschental record. It indicates that albeit low-frequency signals are mitigated in the new reconstruction higher frequency signals are more coherent.

A recently developed reconstruction from the Pyrenees (Büntgen et al. 2024), however, coincides with the new Alpine record (Fig. 5b). Both reconstructions show a similar low-frequency variability. They correlate with 0.4 but correlations reach 0.74 when both chronologies are 300-year low-pass filtered. Although the LIA period is weakly pronounced in both reconstructions, the Pyrenees record does not record a strong cooling in the 1810s, which is prominent in most Alpine records (Table 1) (Schweingruber et al. 1987, 1988; Büntgen et al. 2005, 2006; Frank and Esper 2005b; Corona et al. 2010, 2011; Trachsel et al. 2012; Coppola et al. 2013; Leonelli et al. 2016). The standard deviation of the Pyrenees record in Fig. 5c is stable with lower variation than in the Lötschental record or the record of this study.

A focus on the warm and cold decades of the Lötschental and the new record reveals that while both reconstructions agree on the 1810s as the coldest decade, the warmest decade in the new record is not in the last decades but earlier in the 970 s. The most recent decades, which are the absolute warmest decades according to the Lötschental record, are not the warmest in the new record anymore. Here, the warmest decades are mainly between 960 – 990 CE. Nonetheless, both records show the 2000s as one of the five warmest decades. Surprisingly, the new record does not estimate the 2010s to be among the warmest decades, which is not in line with the instrumental data.

Revising cold and warm extremes in Alpine temperature history

A stable and unbiased variance is crucial for robust extreme year analysis (Frank et al. 2007). Artificial fluctuations of variance can lead to falsely detected extreme events and consequently alter interpretations of extreme year distributions. This influence becomes visible when 30-year high-pass filter reconstructions and their top 20 extreme years are compared (Fig. 6a and c, Table 2). The asterisks in the plots denote the 20 warmest and coldest years of the records, respectively. If variance is not stable in time (Fig. 5b), extreme years like volcanic cooling events are accumulated in periods with increased variance, for example the 11-twelfth century in the Lötschental record. As these volcanic induced cooling years are recorded by tree-ring proxies like MXD, they are frequently used to detect volcanic events (Esper et al. 2017) and reveal information on their impact on temperatures and on their link to societal developments like mass migrations, famines or cultural heydays (Büntgen et al. 2020). A bias in variance can thus alter volcanic event detection and interpretation. Although the Lötschental record detects some of the strongest 20 volcanic SAOD events during the last millennium, a SEA of these events reveals that a highly significant (p < 0.01) cooling in the year of the event and year after the event is not detectable in this record (Fig. 6d). In comparison, the new reconstruction has a temporally more stable variance (Fig. 5b), and the event year as well as following years are found to be significantly cooler than the five years prior to the event (Fig. 6b). As most volcanic events of the SEA fall within periods of lower variance in the Lötschental record, the averaged cooling of these events does not exceed the non-volcanic induced declines in periods when the variance is higher. Further, the growth decline relative to the years prior to the event is less pronounced than in the new record. With fluctuating variance, the volcanic induced cooling years in the Lötschental record are either no outliers in periods of low variance or surrounded by outliers in periods of higher variance. The residuals between the event year temperatures and the reference period are consequently lower.Fig. 6 Extreme Year Analysis: 30-year high-pass filtered reconstructions from this study in a and Lötschental in c. Asterisks denote to the 20 warmest and the 20 coldest years, respectively. Triangles show 20 strongest volcanic events between 880–2000 AD (see Table S4). For these 20 events, panels b and d show the Superposed Epoch Analysis (lag = 5, residuals from the 5 years prior to event) with mean (black line) and 99% confidence intervals (grey) after bootstrap resampling (n = 10,000)

Table 2 Top 20 minimum and maximum extreme years (descending order) of the new Alpine record between 881–2004 CE. Highlighted entries align with an eruption year or the year after using the dates of Table S4

An example of the influence of variance on the interpretation of extreme events is given by focusing on the Samalas eruption in 1257 CE, which is the largest volcanic eruption of the last millennium (Lavigne et al. 2013). The following volcanic cooling is recorded in both reconstructions with -2.32 °C in the new record and -1.18 °C in the Lötschental record (w.r.t. 1252–1256 CE), respectively. While the Samalas induced cooling in 1258 CE presents the strongest negative extreme event in the new record, it does not fall under the extreme events detected in the Lötschental record. This might be explained by variance fluctuations in the high-frequency domain. The resulting loss of climate information might explain why this eruption is not imprinted in the top 20 cold events of the record. Guillet et al. (2017) present historical evidence for bad (grape) harvest, increased prices, heavy precipitation and floodings in the regions of Western European region, where temperatures correlate high with the new record (Fig. 3a). A strong cooling in 1258 CE can be observed in the Pyrenees record being -3.26 °C colder than the average of the previous five years (see Fig. 5b, Büntgen et al. 2024). Other NH and Western Europe reconstructions record decreases between -0.55 to -1.8 °C compared to the ten year mean prior the events (here: -2.52 °C) (Esper et al. 2013; Schneider et al. 2015; Stoffel et al. 2015; Wilson et al. 2016; Hartl-Meier et al. 2017a). Less pronounced cooling in these records compared to the Pyrenees or the Alps is likely due to the observed spatially diverse intensity of the post volcanic cooling (Guillet et al. 2017; Büntgen et al. 2022b). For the new Pyrenees and Alpine reconstructions, the 1258 CE cooling is unrivalled in the last millennium (Figs. 5b, S7). In both, the cooling following the Samalas eruption is stronger than the cooling observed for 1816 CE after Tambora, which is less prominent in the Pyrenees record compared to the Alpine region. The new record agrees with previous studies that showed that summer temperatures in the 1800/1810s have been relatively cold already before the event compared to the situation in the 1250 s (e.g., Frank and Esper 2005b; Corona et al. 2010; Trachsel et al. 2012; Schneider et al. 2015). Although the volcanic induced decrease of 1.01 °C compared to the five years prior to the event appears less severe than the Samalas induced temperature decline, the impact of the 1810s CE cooling on humanity was strong which is proven by historical evidence (e.g., Brázdil et al. 2017; Büntgen et al. 2020). We are aware that the LBM correction can bias the strength of volcanic cooling which is mediated by the reconstruction. In years where LBM events and volcanic events align, LBM induced density decline can amplify the reconstructed compared to the true volcanic cooling. However, as the extreme year 1258 CE is also detected in the new Pyrenees record by Büntgen et al. (2024), which is based on LBM-unaffected pine trees (Pinus uncinata Ramond s. str.), it supports our findings that people in Western Europe experienced a year of unusual cold summer temperatures in 1258 CE.

The top 20 summer heat extremes in the new record are relatively balanced throughout the covered period with clusters in the late medieval times and during LIA. Although long-term temperatures during the LIA are decreasing, exceptionally warm summers relative to surrounding decades are observed in a recent study (Wanner et al. 2022). For example, the summer in 1473 CE was extremely hot and dry in Europe causing early (grape) harvests or harvest failures, loss of transportation routes (rivers), water shortages and wildfires (Camenisch et al. 2020). Likely due to a decrease in variance in this period, this year is not observed as an extreme year in the Lötschental record. Similarly, the 2003 CE heatwave appears as an extreme year when long-term trends are included (Büntgen et al. 2006), but not when the record is high-pass filtered. From 2003–2017 CE no cool extremes are observed in the new Alpine record but no heat extremes either which does not agree with instrumental data. Resuming the fit of the reconstruction model (Fig. 4a) this might result from the models’ limitation to fully capture the instrumental warming after 2010 CE.

Is there a divergence problem in the Alps?

Although the fit between instrumental and proxy data (r = 0.8) is very high, the dataset shows signs of a beginning divergence from 2010 CE onwards (Fig. S7). The divergence phenomenon describes the weakening of the linear relationship between temperature records and tree-ring proxies (TRW and MXD) during the twentieth century (D’Arrigo et al. 2008). This phenomenon has been observed in multiple sites all over the Norther Hemisphere (Wilmking et al. 2005; Wilson et al. 2007; D’Arrigo et al. 2008; Andreu-Hayles et al. 2011; Shah and Shah 2015; St. George and Esper 2019) and in low elevation sites in the European Alps (Corona et al. 2010). However, divergence has not been detected in high elevation sites in this region during the twentieth century for European larch (Larix decidua Mill.) (Büntgen et al. 2008). Here, we find first indications of a developing divergence phenomenon in the Swiss Alps in larch. From 2010 CE onwards MXD does not follow the increased warming observed in instrumental data. This is underlined by a decline in the 30-year running correlations between temperature and MXD indices from 1990 CE until present in Fig. 3c. The linear relationship is weakened and consequently, the prediction performance of linear regression likely decreases with increasing divergence. Correlations with other climate parameters, for example drought or precipitation, do not show a shift in signal in the Alps yet, but might in the future as it has been observed in other regions, for example Corsica (Römer et al. 2021). There, drought, and consequently summer precipitation become more important for limiting tree growth. Precipitation signals are, however, difficult to assess in network datasets like the new Alpine chronology since precipitation patterns in mountainous regions can vary strongly between valleys and elevations (Sevruk 1997; Auer et al. 2001; Schmidli et al. 2002).

Conclusion

In this study, we present a new temperature reconstruction based on tree-ring data from the Swiss Alps. We introduce a new approach to improve temperature reconstructions and their temporal reliability by considering elevational differences among data from historical buildings, removing LBM signals with improved methods, and applying novel approaches to stabilize variance. Albeit low-frequency trends are muted in our new reconstruction compared to previous records, high-frequency temperature variability shows strong agreements with large volcanic eruptions which is in line with historical evidence. We show how artificial variance increases may alter extreme year analysis including the 1257 CE Samalas eruption, which led to stronger regional cooling than the 1815 CE Tambora eruption. Our study thereby emphasizes how decision making, data selection and treatment in dendroclimatology influences the outcome of a reconstruction approach (Büntgen et al. 2021). The different applied methods, in particular the (1) data selection (of proxy and target), (2) mean adjustment, (3) LBM correction, (4) variance stabilization, and (5) the chosen method of reconstruction, alter the resulting record and interpretation of temperature variability towards its extreme events and long-term temperature signals. This study highlights the urgency to continue the development and testing of new methods to improve the temporal stability and reliability of climate reconstructions. It stresses the importance that with new methods revising existing records to improve our knowledge about past climate variability is indispensable.

Supplementary Information

Below is the link to the electronic supplementary material.Supplementary file1 (PDF 1278 KB)

Acknowledgements

We thank Markus Kochbeck, Frederick Reinig, Philipp Römer, Christian Gnaneswaran, Benedikt Lang, Philipp Schulz, Lara Meurer, Sophie Spelsberg and Jannes Fischer for help in the laboratory and field. We thank Martin Schmidhalter from Dendrolabor Valais for the provided data.

Author contributions

All authors contributed significantly to the study conception and design. Material preparation, data collection and analysis were performed by Eileen Kuhl, Claudia Hartl, Jan Esper, Ulf Büntgen, Lea Schneider, Marcel Kunz, Lara Klippel and Valerie Trouet. The first draft of the manuscript was written by Eileen Kuhl and all authors commented on previous versions of the manuscript. All authors read and approved the final manuscript.

Funding

Open Access funding enabled and organized by Projekt DEAL. This research was funded by the German Research Foundation [grant number HA 8048/1–1]. J.E. and U.B. received funding from the ERC Advanced project Monostar (AdG 882727) and the project AdAgrif (CZ.02.01.01/00/22_008/0004635). M.K. received funding from the ERC Advanced project Monostar (AdG 882727).

Data availability

Raw tree-ring measurements will be submitted to the International Tree Ring Database (ITRDB).

Declarations

Competing interests

The authors declare that they have no relevant financial or non-financial conflict of interest.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
==== Refs
References

Affolter P, Büntgen U, Esper J, Rigling A, Weber P, Luterbacher J, Frank D (2010) Inner Alpine conifer response to 20th century drought swings. Eur J For Res 129:2889–298. 10.1007/s10342-009-0327-x
Allgaier Leuch B Streit K Brang P Der Schweizer Wald im Klimawandel: welche Entwicklungen kommen auf uns zu? Merkblatt Für Die Praxis 2017 59 1 12
Allgaier Leuch B, Streit K, Brang P (2017) Der Schweizer Wald im Klimawandel: welche Entwicklungen kommen auf uns zu? Merkblatt Für Die Praxis 59:1–12
Anchukaitis KJ Wilson R Briffa KR Büntgen U Cook ER D’Arrigo R Davi N Esper J Frank D Gunnarson BE Hegerl G Helama S Klesse S Krusic PJ Linderholm HW Myglan V Osborn TJ Zhang P Rydval M Schneider L Schurer A Wiles G Zorita E Last millennium Northern Hemisphere summer temperatures from tree rings: part II, spatially resolved reconstructions Quat Sci Rev 2017 163 1 22 10.1016/j.quascirev.2017.02.020
Anchukaitis KJ, Wilson R, Briffa KR, Büntgen U, Cook ER, D’Arrigo R, Davi N, Esper J, Frank D, Gunnarson BE, Hegerl G, Helama S, Klesse S, Krusic PJ, Linderholm HW, Myglan V, Osborn TJ, Zhang P, Rydval M, Schneider L, Schurer A, Wiles G, Zorita E (2017) Last millennium Northern Hemisphere summer temperatures from tree rings: part II, spatially resolved reconstructions. Quat Sci Rev 163:1–22. 10.1016/j.quascirev.2017.02.020
Andreu-Hayles L, D’Arrigo R, Anchukaitis K, Beck P, Frank D, Goetz S (2011) Varying boreal forest response to Arctic environmental change at the Firth River, Alaska. Environ Res Lett. 10.1088/1748-9326/6/4/045503
Aryal S Grießinger J Arsalani M Meier WJ-H Fu P-L Fan Z-X Bräuning A Insect infestations have an impact on the quality of climate reconstructions using Larix ring-width chronologies from the Tibetan plateau Ecol Indic 2023 148 110124 10.1016/j.ecolind.2023.110124
Aryal S, Grießinger J, Arsalani M, Meier WJ-H, Fu P-L, Fan Z-X, Bräuning A (2023) Insect infestations have an impact on the quality of climate reconstructions using Larix ring-width chronologies from the Tibetan plateau. Ecol Indic 148:110124. 10.1016/j.ecolind.2023.110124
Auer I Böhm R Maugeri M A new long-term gridded precipitation data-set for the alps and its application for Map and Alpclim Phys Chem Earth B 2001 26 5–6 421 424 10.1016/S1464-1909(01)00029-6
Auer I, Böhm R, Maugeri M (2001) A new long-term gridded precipitation data-set for the alps and its application for Map and Alpclim. Phys Chem Earth B 26(5–6):421–424. 10.1016/S1464-1909(01)00029-6
Baltensweiler W Weber UM Cherubini P Tracing the influence of larch-bud-moth insect outbreaks and weather conditions on larch tree-ring growth in Engadine (Switzerland) Oikos 2008 117 2 161 172 10.1111/j.2007.0030-1299.16117.x
Baltensweiler W, Weber UM, Cherubini P (2008) Tracing the influence of larch-bud-moth insect outbreaks and weather conditions on larch tree-ring growth in Engadine (Switzerland). Oikos 117(2):161–172. 10.1111/j.2007.0030-1299.16117.x
Björklund J, von Arx G, Nievergelt D, Wilson RJS, Van den Bulcke J, Günther B, Loader NJ, Rydval M, Fonti P, Scharnweber T, Andreu‐Hayles L, Büntgen U, D’Arrigo R, Davi N, De Mil T, Esper J, Gärtner H, Geary J, Gunnarson BE, Hartl C, Hevia A, Song H, Janecka K, Kaczka RJ, Kirdyanov AV, Kochbeck M, Liu Y, Meko M, Mundo I, Nicolussi K, Oelkers R, Pichler T, Sánchez‐Salguero R, Schneider L, Schweingruber F, Timonen M, Trouet V, Van Acker J, Verstege A, Villalba R, Wilmking M, Frank D (2019) Scientific merits and analytical challenges of tree‐ring densitometry. Rev Geophys 57(4). 10.1029/2019RG000642
Borkotoky SS, Williams AP, Cook ER, Steinschneider S (2021) Reconstructing extreme precipitation in the Sacramento River watershed using tree‐ring based proxies of cold‐season precipitation. Water Resour Res 57:e2020WR028824. 10.1029/2020WR028824
Bräker OU (1981) Der Alterstrend bei Jahrringdichten und Jahrringbreiten von Nadelhölzern und sein Ausgleich. Mitt Forstl Bundes-Vers.anst Wien 142:75–102
Brázdil R Dobrovolný P Luterbacher J Moberg A Pfister C Wheeler D Zorita E European climate of the past 500 years: new challenges for historical climatology Clim Change 2010 101 1–2 7 40 10.1007/s10584-009-9783-z
Brázdil R, Dobrovolný P, Luterbacher J, Moberg A, Pfister C, Wheeler D, Zorita E (2010) European climate of the past 500 years: new challenges for historical climatology. Clim Change 101(1–2):7–40. 10.1007/s10584-009-9783-z
Brázdil R, Reznícková L, Valásek H, Dolák L, Oldrich (2017) Climatic and other responses to the Lakagígar 1783 and Tambora 1815 volcanic eruptions in the Czech Lands. Geografie 122:147–168. 10.37040/geografie2017122020147
Briffa KR Jones PD Bartholin TS Eckstein D Schweingruber F Karlén W Zetterberg P Eronen M Fennoscandian summers from ad 500: temperature changes on short and long timescales Clim Dyn 1992 7 3 111 119 10.1007/BF00211153
Briffa KR, Jones PD, Bartholin TS, Eckstein D, Schweingruber F, Karlén W, Zetterberg P, Eronen M (1992) Fennoscandian summers from ad 500: temperature changes on short and long timescales. Clim Dyn 7(3):111–119. 10.1007/BF00211153
Briffa KR Osborn TJ Schweingruber F Harris IC Jones PD Shiyatov SG Vaganov EA Low-frequency temperature variations from a northern tree ring density network J Geophys Res 2001 106 D3 2929 2941 10.1029/2000JD900617
Briffa KR, Osborn TJ, Schweingruber F, Harris IC, Jones PD, Shiyatov SG, Vaganov EA (2001) Low-frequency temperature variations from a northern tree ring density network. J Geophys Res 106(D3):2929–2941. 10.1029/2000JD900617
Bunn AG Statistical and visual crossdating in R using the dplR library Dendrochronologia 2010 28 4 251 258 10.1016/j.dendro.2009.12.001
Bunn AG (2010) Statistical and visual crossdating in R using the dplR library. Dendrochronologia 28(4):251–258. 10.1016/j.dendro.2009.12.001
Büntgen U Esper J Frank D Nicolussi K Schmidhalter M A 1052-year tree-ring proxy for Alpine summer temperatures Clim Dyn 2005 25 2–3 141 153 10.1007/s00382-005-0028-1
Büntgen U, Esper J, Frank D, Nicolussi K, Schmidhalter M (2005) A 1052-year tree-ring proxy for Alpine summer temperatures. Clim Dyn 25(2–3):141–153. 10.1007/s00382-005-0028-1
Büntgen U Frank D Wilson RJS Carrer M Urbinati C Esper J Testing for tree-ring divergence in the European Alps: tree-ring divergence in the European Alps Glob Change Biol 2008 14 10 2443 2453 10.1111/j.1365-2486.2008.01640.x
Büntgen U, Frank D, Wilson RJS, Carrer M, Urbinati C, Esper J (2008) Testing for tree-ring divergence in the European Alps: tree-ring divergence in the European Alps. Glob Change Biol 14(10):2443–2453. 10.1111/j.1365-2486.2008.01640.x
Büntgen U Frank D Liebhold A Johnson D Carrer M Urbinati C Grabner M Nicolussi K Levanic T Esper J Three centuries of insect outbreaks across the European Alps New Phytol 2009 182 4 929 941 10.1111/j.1469-8137.2009.02825.x 19383093
Büntgen U, Frank D, Liebhold A, Johnson D, Carrer M, Urbinati C, Grabner M, Nicolussi K, Levanic T, Esper J (2009) Three centuries of insect outbreaks across the European Alps. New Phytol 182(4):929–941. 10.1111/j.1469-8137.2009.02825.x19383093
Büntgen U Piermattei A Crivellaro A Reinig F Krusic PJ Trnka M Torbenson M Esper J Common Era treeline fluctuations and their implications for climate reconstructions Glob Planet Change 2022 219 103979 10.1016/j.gloplacha.2022.103979
Büntgen U, Piermattei A, Crivellaro A, Reinig F, Krusic PJ, Trnka M, Torbenson M, Esper J (2022a) Common Era treeline fluctuations and their implications for climate reconstructions. Glob Planet Change 219:103979. 10.1016/j.gloplacha.2022.103979
Büntgen U Smith SH Wagner S Krusic P Esper J Piermattei A Crivellaro A Reinig F Tegel W Kirdyanov A Trnka M Oppenheimer C Global tree-ring response and inferred climate variation following the mid-thirteenth century Samalas eruption Clim Dyn 2022 59 1–2 531 546 10.1007/s00382-022-06141-3
Büntgen U, Smith SH, Wagner S, Krusic P, Esper J, Piermattei A, Crivellaro A, Reinig F, Tegel W, Kirdyanov A, Trnka M, Oppenheimer C (2022b) Global tree-ring response and inferred climate variation following the mid-thirteenth century Samalas eruption. Clim Dyn 59(1–2):531–546. 10.1007/s00382-022-06141-3
Büntgen U Reinig F Verstege A Piermattei A Kunz M Krusic P Slavin P Stěpanék P Torbenson M Martinez del Castillo E Arosio T Kirdyanov A Oppenheimer C Trnka M Palosse A Bebchuk T Camarero JJ Esper J Recent summer warming over the western Mediterranean region is unprecedented since medieval times Glob Planet Change 2024 232 104336 10.1016/j.gloplacha.2023.104336
Büntgen U, Reinig F, Verstege A, Piermattei A, Kunz M, Krusic P, Slavin P, Stěpanék P, Torbenson M, Martinez del Castillo E, Arosio T, Kirdyanov A, Oppenheimer C, Trnka M, Palosse A, Bebchuk T, Camarero JJ, Esper J (2024) Recent summer warming over the western Mediterranean region is unprecedented since medieval times. Glob Planet Change 232:104336. 10.1016/j.gloplacha.2023.104336
Büntgen U, Frank D, Nievergelt D, Esper J (2006) Summer Temperature Variations in the European Alps, a.d. 755–2004. J Clim 19(21):5606–5623. 10.1175/JCLI3917.1
Büntgen U, Arseneault D, Boucher E, Churakova (Sidorova) OV, Gennaretti F, Crivellaro A, Hughes MK, Kirdyanov AV, Klippel L, Krusic PJ, Linderholm HW, Ljungqvist FC, Ludescher J, McCormick M, Myglan VS, Nicolussi K, Piermattei A, Oppenheimer C, Reinig F, Sigl M, Vaganov EA, Esper J (2020) Prominent role of volcanism in Common Era climate variability and human history. Dendrochronologia 64:125757. 10.1016/j.dendro.2020.125757
Büntgen U, Allen K, Anchukaitis KJ, Arseneault D, Boucher E, Bräuning A, Chatterjee S, Cherubini P, Churakova OV, Corona C, Gennaretti F, Grießinger J, Guillet S, Guiot J, Gunnarson B, Helama S, Hochreuther P, Hughes MK, Huybers P, Kirdyanov AV, Krusic PJ, Ludescher J, Meier WJ-H, Myglan VS, Nicolussi K, Oppenheimer C, Reinig F, Salzer MW, Seftigen K, Stine AR, Stoffel M, St. George S, Tejedor E, Trevino A, Trouet V, Wang J, Wilson R, Yang B, Xu G, Esper J (2021) The influence of decision-making in tree ring-based climate reconstructions. Nat Commun 12(1):3411. 10.1038/s41467-021-23627-6
Camenisch C Brázdil R Kiss A Pfister C Wetter O Rohr C Contino A Retsö D Extreme heat and drought in 1473 and their impacts in Europe in the context of the early 1470s Reg Environ Change 2020 20 1 19 10.1007/s10113-020-01601-0
Camenisch C, Brázdil R, Kiss A, Pfister C, Wetter O, Rohr C, Contino A, Retsö D (2020) Extreme heat and drought in 1473 and their impacts in Europe in the context of the early 1470s. Reg Environ Change 20(1):19. 10.1007/s10113-020-01601-0
Conedera M Colombaroli D Tinner W Krebs P Whitlock C Insights about past forest dynamics as a tool for present and future forest management in Switzerland For Ecol Manag 2017 388 100 112 10.1016/j.foreco.2016.10.027
Conedera M, Colombaroli D, Tinner W, Krebs P, Whitlock C (2017) Insights about past forest dynamics as a tool for present and future forest management in Switzerland. For Ecol Manag 388:100–112. 10.1016/j.foreco.2016.10.027
Cook ER Peters K The smoothing spline: a new approach to standardizing forest interior tree-ring width series for dendroclimatic studies Tree-Ring Bull 1981 41 45 53
Cook ER, Peters K (1981) The smoothing spline: a new approach to standardizing forest interior tree-ring width series for dendroclimatic studies. Tree-Ring Bull 41:45–53
Cook ER Peters K Calculating unbiased tree-ring indices for the study of climatic and environmental change Holocene 1997 7 3 361 370 10.1177/095968369700700314
Cook ER, Peters K (1997) Calculating unbiased tree-ring indices for the study of climatic and environmental change. Holocene 7(3):361–370. 10.1177/095968369700700314
Cook ER Briffa KR Shiyatov SG Mazepa V Cook ER Kairiukstis LA Tree-ring standardization and growth-trend estimation Methods of dendrochronology: applications in the environmental sciences 1990 1 Dordrecht Kluwer Academic Publishers 104 123
Cook ER, Briffa KR, Shiyatov SG, Mazepa V (1990) Tree-ring standardization and growth-trend estimation. In: Cook ER, Kairiukstis LA (eds) Methods of dendrochronology: applications in the environmental sciences, 1st edn. Kluwer Academic Publishers, Dordrecht, pp 104–123
Cook ER, Krusic PJ, Peters K, Melvin TM (2017) Program signal free (Version 45_v2b). RCS signal free tree-ring standardization program. Tree-ring laboratory of lamont-doherty earth observatory. https://www.geog.cam.ac.uk/research/projects/dendrosoftware/
Cook ER (1985) A time series analysis approach to tree ring standardization. Dissertation, University of Arizona
Coolidge WAB (1912) The names of Zermatt. Engl Hist Rev XXVII(CVII):522–530. 10.1093/ehr/XXVII.CVII.522
Coppola A Leonelli G Salvatore MC Pelfini M Baroni C Tree-ring–based summer mean temperature variations in the Adamello-Presanella Group (Italian Central Alps), 1610–2008 AD Clim Past 2013 9 1 211 221 10.5194/cp-9-211-2013
Coppola A, Leonelli G, Salvatore MC, Pelfini M, Baroni C (2013) Tree-ring–based summer mean temperature variations in the Adamello-Presanella Group (Italian Central Alps), 1610–2008 AD. Clim Past 9(1):211–221. 10.5194/cp-9-211-2013
Corona C Guiot J Edouard JL Chalié F Büntgen U Nola P Urbinati C Millennium-long summer temperature variations in the European Alps as reconstructed from tree rings Clim Past 2010 6 3 379 400 10.5194/cp-6-379-2010
Corona C, Guiot J, Edouard JL, Chalié F, Büntgen U, Nola P, Urbinati C (2010) Millennium-long summer temperature variations in the European Alps as reconstructed from tree rings. Clim Past 6(3):379–400. 10.5194/cp-6-379-2010
Corona C Edouard J-L Guibal F Guiot J Bernard S Thomas A Denelle N Long-term summer (AD751–2008) temperature fluctuation in the French Alps based on tree-ring data Boreas 2011 40 351 366 10.1111/j.1502-3885.2010.00185.x
Corona C, Edouard J-L, Guibal F, Guiot J, Bernard S, Thomas A, Denelle N (2011) Long-term summer (AD751–2008) temperature fluctuation in the French Alps based on tree-ring data. Boreas 40:351–366. 10.1111/j.1502-3885.2010.00185.x
D’Arrigo R Wilson RJS Liepert B Cherubini P On the ‘Divergence Problem’ in Northern Forests: a review of the tree-ring evidence and possible causes Glob Planet Change 2008 60 3–4 289 305 10.1016/j.gloplacha.2007.03.004
D’Arrigo R, Wilson RJS, Liepert B, Cherubini P (2008) On the ‘Divergence Problem’ in Northern Forests: a review of the tree-ring evidence and possible causes. Glob Planet Change 60(3–4):289–305. 10.1016/j.gloplacha.2007.03.004
Düthorn E Holzkämper S Timonen M Esper J Influence of micro-site conditions on tree-ring climate signals and trends in central and northern Sweden Trees 2013 27 5 1395 1404 10.1007/s00468-013-0887-8
Düthorn E, Holzkämper S, Timonen M, Esper J (2013) Influence of micro-site conditions on tree-ring climate signals and trends in central and northern Sweden. Trees 27(5):1395–1404. 10.1007/s00468-013-0887-8
Düthorn E, Schneider L, Konter O, Schön P, Timonen M, Esper J (2015) On the hidden significance of differing micro-sites on tree-ring based climate reconstructions. Silva Fenn 49(1). 10.14214/sf.1220
Esper J Cook ER Krusic PJ Peters K Schweingruber F Tests of the RCS method for preserving low-frequency variability in long tree-ring chronologies Tree-Ring Res 2003 59 2 81 89
Esper J, Cook ER, Krusic PJ, Peters K, Schweingruber F (2003) Tests of the RCS method for preserving low-frequency variability in long tree-ring chronologies. Tree-Ring Res 59(2):81–89
Esper J Büntgen U Frank D Nievergelt D Liebhold A 1200 years of regular outbreaks in alpine insects Proc R Soc B 2007 274 1610 671 679 10.1098/rspb.2006.0191
Esper J, Büntgen U, Frank D, Nievergelt D, Liebhold A (2007a) 1200 years of regular outbreaks in alpine insects. Proc R Soc B 274(1610):671–679. 10.1098/rspb.2006.0191
Esper J Frank D Büntgen U Verstege A Luterbacher J Xoplaki E Long-term drought severity variations in Morocco Geophys Res Lett 2007 34 17 L17702 10.1029/2007GL030844
Esper J, Frank D, Büntgen U, Verstege A, Luterbacher J, Xoplaki E (2007c) Long-term drought severity variations in Morocco. Geophys Res Lett 34(17):L17702. 10.1029/2007GL030844
Esper J Schneider L Krusic PJ Luterbacher J Büntgen U Timonen M Sirocko F Zorita E European summer temperature response to annually dated volcanic eruptions over the past nine centuries Bull Volcanol 2013 75 7 736 10.1007/s00445-013-0736-z
Esper J, Schneider L, Krusic PJ, Luterbacher J, Büntgen U, Timonen M, Sirocko F, Zorita E (2013) European summer temperature response to annually dated volcanic eruptions over the past nine centuries. Bull Volcanol 75(7):736. 10.1007/s00445-013-0736-z
Esper J Krusic PJ Ljungqvist FC Luterbacher J Carrer M Cook ER Davi N Hartl-Meier C Kirdyanov AV Konter O Myglan VS Timonen M Treydte K Trouet V Villalba R Yang B Büntgen U Ranking of tree-ring based temperature reconstructions of the past millennium Quat Sci Rev 2016 145 134 151 10.1016/j.quascirev.2016.05.009
Esper J, Krusic PJ, Ljungqvist FC, Luterbacher J, Carrer M, Cook ER, Davi N, Hartl-Meier C, Kirdyanov AV, Konter O, Myglan VS, Timonen M, Treydte K, Trouet V, Villalba R, Yang B, Büntgen U (2016) Ranking of tree-ring based temperature reconstructions of the past millennium. Quat Sci Rev 145:134–151. 10.1016/j.quascirev.2016.05.009
Esper J Büntgen U Hartl-Meier C Oppenheimer C Schneider L Northern Hemisphere temperature anomalies during the 1450s period of ambiguous volcanic forcing Bull Volcanol 2017 79 6 41 10.1007/s00445-017-1125-9
Esper J, Büntgen U, Hartl-Meier C, Oppenheimer C, Schneider L (2017) Northern Hemisphere temperature anomalies during the 1450s period of ambiguous volcanic forcing. Bull Volcanol 79(6):41. 10.1007/s00445-017-1125-9
Esper J, Büntgen U (2021) The future of paleoclimate. Clim Res 83:57–59.10.3354/cr01636
Esper J, Büntgen U, Frank D, Pichler T, Nicolussi K (2007b) Updating the Tyrol tree-ring dataset. In: Haneca K, Verheyden A, Beeckman H, Gärtner H, Helle G, Schleser G (eds) Schirften des Forschungszentrums Jülich. Forschungszentrum Jülich, Tervuren, Belgium, pp 80–84
Eyring V Cox PM Flato GM Gleckler PJ Abramowitz G Caldwell P Collins WD Gier BK Hall AD Hoffman FM Hurtt GC Jahn A Jones CD Klein SA Krasting JP Kwiatkowski L Lorenz R Maloney E Meehl GA Pendergrass AG Pincus R Ruane AC Russell JL Sanderson BM Santer BD Sherwood SC Simpson IR Stouffer RJ Williamson MS Taking climate model evaluation to the next level Nat Clim Change 2019 9 2 102 110 10.1038/s41558-018-0355-y
Eyring V, Cox PM, Flato GM, Gleckler PJ, Abramowitz G, Caldwell P, Collins WD, Gier BK, Hall AD, Hoffman FM, Hurtt GC, Jahn A, Jones CD, Klein SA, Krasting JP, Kwiatkowski L, Lorenz R, Maloney E, Meehl GA, Pendergrass AG, Pincus R, Ruane AC, Russell JL, Sanderson BM, Santer BD, Sherwood SC, Simpson IR, Stouffer RJ, Williamson MS (2019) Taking climate model evaluation to the next level. Nat Clim Change 9(2):102–110. 10.1038/s41558-018-0355-y
Frank D Esper J Characterization and climate response patterns of a high-elevation, multi-species tree-ring network in the European Alps Dendrochronologia 2005 22 2 107 121 10.1016/j.dendro.2005.02.004
Frank D, Esper J (2005a) Characterization and climate response patterns of a high-elevation, multi-species tree-ring network in the European Alps. Dendrochronologia 22(2):107–121. 10.1016/j.dendro.2005.02.004
Frank D Esper J Temperature reconstructions and comparisons with instrumental data from a tree-ring network for the European Alps Int J Climatol 2005 25 11 1437 1454 10.1002/joc.1210
Frank D, Esper J (2005b) Temperature reconstructions and comparisons with instrumental data from a tree-ring network for the European Alps. Int J Climatol 25(11):1437–1454. 10.1002/joc.1210
Frank D Wilson R Esper J Synchronous variability changes in Alpine temperature and tree-ring data over the past two centuries Boreas 2005 34 498 505 10.1080/03009480500231443
Frank D, Wilson R, Esper J (2005) Synchronous variability changes in Alpine temperature and tree-ring data over the past two centuries. Boreas 34:498–505. 10.1080/03009480500231443
Frank D, Esper J, Cook ER (2007) Adjustment for proxy number and coherence in a large-scale temperature reconstruction. Geophys Res Lett 34(16). 10.1029/2007GL030571
Grove JM The initiation of the" Little Ice Age" in regions round the North Atlantic Clim Change 2001 48 53 82 10.1023/A:1005662822136
Grove JM (2001) The initiation of the" Little Ice Age" in regions round the North Atlantic. Clim Change 48:53–82. 10.1023/A:1005662822136
Guillet S, Corona C, Stoffel M, Khodri M, Lavigne F, Ortega P, Eckert N, Sielenou PD, Daux V, Churakova (Sidorova) OV, Davi N, Edouard J-L, Zhang Y, Luckman BH, Myglan VS, Guiot J, Beniston M, Masson-Delmotte V, Oppenheimer C (2017) Climate response to the Samalas volcanic eruption in 1257 revealed by proxy records. Nat Geosci 10(2):123–128. 10.1038/ngeo2875
Gunnarson BE Linderholm HW Moberg A Improving a tree-ring reconstruction from west-central Scandinavia: 900 years of warm-season temperatures Clim Dyn 2011 36 1–2 97 108 10.1007/s00382-010-0783-5
Gunnarson BE, Linderholm HW, Moberg A (2011) Improving a tree-ring reconstruction from west-central Scandinavia: 900 years of warm-season temperatures. Clim Dyn 36(1–2):97–108. 10.1007/s00382-010-0783-5
Harris I Osborn TJ Jones P Lister D Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset Sci Data 2020 7 1 109 10.1038/s41597-020-0453-3 32246091
Harris I, Osborn TJ, Jones P, Lister D (2020) Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset. Sci Data 7(1):109. 10.1038/s41597-020-0453-332246091
Hartl C Düthorn E Tejedor E Kirchhefer AJ Timonen M Holzkämper S Büntgen U Esper J Micro-site conditions affect Fennoscandian forest growth Dendrochronologia 2021 65 125787 10.1016/j.dendro.2020.125787
Hartl C, Düthorn E, Tejedor E, Kirchhefer AJ, Timonen M, Holzkämper S, Büntgen U, Esper J (2021) Micro-site conditions affect Fennoscandian forest growth. Dendrochronologia 65:125787. 10.1016/j.dendro.2020.125787
Hartl C Schneider L Riechelmann DFC Kuhl E Kochbeck M Klippel L Büntgen U Esper J The temperature sensitivity along elevational gradients is more stable in maximum latewood density than tree-ring width Dendrochronologia 2022 73 125958 10.1016/j.dendro.2022.125958
Hartl C, Schneider L, Riechelmann DFC, Kuhl E, Kochbeck M, Klippel L, Büntgen U, Esper J (2022) The temperature sensitivity along elevational gradients is more stable in maximum latewood density than tree-ring width. Dendrochronologia 73:125958. 10.1016/j.dendro.2022.125958
Hartl-Meier C Dittmar C Zang C Rothe A Mountain forest growth response to climate change in the Northern Limestone Alps Trees 2014 28 3 819 829 10.1007/s00468-014-0994-1
Hartl-Meier C, Dittmar C, Zang C, Rothe A (2014a) Mountain forest growth response to climate change in the Northern Limestone Alps. Trees 28(3):819–829. 10.1007/s00468-014-0994-1
Hartl-Meier C Zang C Dittmar C Esper J Göttlein A Rothe A Vulnerability of Norway spruce to climate change in mountain forests of the European Alps Clim Res 2014 60 2 119 132 10.3354/cr01226
Hartl-Meier C, Zang C, Dittmar C, Esper J, Göttlein A, Rothe A (2014b) Vulnerability of Norway spruce to climate change in mountain forests of the European Alps. Clim Res 60(2):119–132. 10.3354/cr01226
Hartl-Meier C Büntgen U Smerdon JE Zorita E Krusic PJ Ljungqvist FC Schneider L Esper J Temperature covariance in tree ring reconstructions and model simulations over the past millennium: last-millennium temperature covariance Geophys Res Lett 2017 44 18 9458 9469 10.1002/2017GL073239
Hartl-Meier C, Büntgen U, Smerdon JE, Zorita E, Krusic PJ, Ljungqvist FC, Schneider L, Esper J (2017a) Temperature covariance in tree ring reconstructions and model simulations over the past millennium: last-millennium temperature covariance. Geophys Res Lett 44(18):9458–9469. 10.1002/2017GL073239
Hartl-Meier C Esper J Liebhold A Konter O Rothe A Büntgen U Effects of host abundance on larch budmoth outbreaks in the European Alps Agr for Entomol 2017 19 4 376 387 10.1111/afe.12216
Hartl-Meier C, Esper J, Liebhold A, Konter O, Rothe A, Büntgen U (2017b) Effects of host abundance on larch budmoth outbreaks in the European Alps. Agr for Entomol 19(4):376–387. 10.1111/afe.12216
Hartl-Meier C, Büntgen U, Esper J (2016) On the occurrence of cyclic larch budmoth outbreaks beyond its geographical hotspots. In: Hevia A, Sánchez-Salguero R, Linares JC, Olano JM, Camarero JJ, Gutiérrez E, Helle G, Gärtner H (eds) Proceedings of the DENDROSYMPOSIUM 2015: May 20th - 23rd, 2015 in Sevilla, Spain, (Scientific Technical Report ; 16/04), 14th TRACE conference (Tree rings in archaeology, climatology and ecology). Deutsches GeoForschungsZentrum GFZ, (Sevilla, Spain), Potsdam, pp 86–92
Johnson DM Bjørnstad ON Liebhold AM Landscape geometry and travelling waves in the larch budmoth: landscape geometry and travelling waves in the LBM Ecol Lett 2004 7 10 967 974 10.1111/j.1461-0248.2004.00659.x
Johnson DM, Bjørnstad ON, Liebhold AM (2004) Landscape geometry and travelling waves in the larch budmoth: landscape geometry and travelling waves in the LBM. Ecol Lett 7(10):967–974. 10.1111/j.1461-0248.2004.00659.x
King GM Gugerli F Fonti P Frank D Tree growth response along an elevational gradient: climate or genetics? Oecologia 2013 173 4 1587 1600 10.1007/s00442-013-2696-6 23771802
King GM, Gugerli F, Fonti P, Frank D (2013) Tree growth response along an elevational gradient: climate or genetics? Oecologia 173(4):1587–1600. 10.1007/s00442-013-2696-623771802
Kuhl E, Zang C, Esper J, Riechelmann DFC, Büntgen U, Briesch M, Reinig F, Römer P, Konter O, Schmidhalter M, Hartl C (2023) Using machine learning on tree-ring data to determine the geographical provenance of historical construction timbers. Ecosphere 1–14. 10.1002/ecs2.4453
Kunz M Esper J Kuhl E Schneider L Büntgen U Hartl C Combining tree-ring width and density to separate the effects of climate 2 variation and insect defoliation Forests 2023 14 1478 10.3390/f14071478
Kunz M, Esper J, Kuhl E, Schneider L, Büntgen U, Hartl C (2023) Combining tree-ring width and density to separate the effects of climate 2 variation and insect defoliation. Forests 14:1478. 10.3390/f14071478
Lavigne F, Degeai J-P, Komorowski J-C, Guillet S, Robert V, Lahitte P, Oppenheimer C, Stoffel M, Vidal CM, Surono, Pratomo I, Wassmer P, Hajdas I, Hadmoko DS, de Belizal E (2013) Source of the great A.D. 1257 mystery eruption unveiled, Samalas volcano, Rinjani Volcanic Complex, Indonesia. Proc Natl Acad Sci USA 110(42):16742–16747. 10.1073/pnas.1307520110
Leonelli G Coppola A Baroni C Salvatore MC Maugeri M Brunetti M Pelfini M Multispecies dendroclimatic reconstructions of summer temperature in the European Alps enhanced by trees highly sensitive to temperature Clim Change 2016 137 1–2 275 291 10.1007/s10584-016-1658-5
Leonelli G, Coppola A, Baroni C, Salvatore MC, Maugeri M, Brunetti M, Pelfini M (2016) Multispecies dendroclimatic reconstructions of summer temperature in the European Alps enhanced by trees highly sensitive to temperature. Clim Change 137(1–2):275–291. 10.1007/s10584-016-1658-5
Lough JM, Fritts HC (1987) An assessment of the possible effects of volcanic eruptions on North American climate using tree-ring data, 1602 to 1900 A.D. Clim Change 10(3):219–239. 10.1007/BF00143903
Luterbacher J Werner JP Smerdon JE Fernández-Donado L González-Rouco FJ Barriopedro D Ljungqvist FC Büntgen U Zorita E Wagner S Esper J McCarroll D Toreti A Frank D Jungclaus JH Barriendos M Bertolin C Bothe O Brázdil R Camuffo D Dobrovolný P Gagen M García-Bustamante E Ge Q Gómez-Navarro JJ Guiot J Hao Z Hegerl GC Holmgren K Klimenko VV Martín-Chivelet J Pfister C Roberts N Schindler A Schurer A Solomina O von Gunten L Wahl E Wanner H Wetter O Xoplaki E Yuan N Zanchettin D Zhang H Zerefos C European summer temperatures since Roman times Environ Res Lett 2016 11 2 024001 10.1088/1748-9326/11/2/024001
Luterbacher J, Werner JP, Smerdon JE, Fernández-Donado L, González-Rouco FJ, Barriopedro D, Ljungqvist FC, Büntgen U, Zorita E, Wagner S, Esper J, McCarroll D, Toreti A, Frank D, Jungclaus JH, Barriendos M, Bertolin C, Bothe O, Brázdil R, Camuffo D, Dobrovolný P, Gagen M, García-Bustamante E, Ge Q, Gómez-Navarro JJ, Guiot J, Hao Z, Hegerl GC, Holmgren K, Klimenko VV, Martín-Chivelet J, Pfister C, Roberts N, Schindler A, Schurer A, Solomina O, von Gunten L, Wahl E, Wanner H, Wetter O, Xoplaki E, Yuan N, Zanchettin D, Zhang H, Zerefos C (2016) European summer temperatures since Roman times. Environ Res Lett 11(2):024001. 10.1088/1748-9326/11/2/024001
Lyu S Li Z Zhang Y Wang X A 414-year tree-ring-based April–July minimum temperature reconstruction and its implications for the extreme climate events, northeast China Clim past 2016 12 1879 1888 10.5194/cp-12-1879-2016
Lyu S, Li Z, Zhang Y, Wang X (2016) A 414-year tree-ring-based April–July minimum temperature reconstruction and its implications for the extreme climate events, northeast China. Clim past 12:1879–1888. 10.5194/cp-12-1879-2016
Melvin TM Briffa KR A “signal-free” approach to dendroclimatic standardisation Dendrochronologia 2008 26 2 71 86 10.1016/j.dendro.2007.12.001
Melvin TM, Briffa KR (2008) A “signal-free” approach to dendroclimatic standardisation. Dendrochronologia 26(2):71–86. 10.1016/j.dendro.2007.12.001
Melvin TM Briffa KR Nicolussi K Grabner M Time-varying-response smoothing Dendrochronologia 2007 25 1 65 69 10.1016/j.dendro.2007.01.004
Melvin TM, Briffa KR, Nicolussi K, Grabner M (2007) Time-varying-response smoothing. Dendrochronologia 25(1):65–69. 10.1016/j.dendro.2007.01.004
Melvin TM Grudd H Briffa KR Potential bias in ‘updating’ tree-ring chronologies using regional curve standardisation: re-processing 1500 years of Torneträsk density and ring-width data Holocene 2013 23 3 364 373 10.1177/0959683612460791
Melvin TM, Grudd H, Briffa KR (2013) Potential bias in ‘updating’ tree-ring chronologies using regional curve standardisation: re-processing 1500 years of Torneträsk density and ring-width data. Holocene 23(3):364–373. 10.1177/0959683612460791
Neuwirth B Esper J Schweingruber FH Winiger M Site ecological differences to the climatic forcing of spruce pointer years from the LoÈ tschental, Switzerland Dendrochronologia 2004 21 2 69 78 10.1078/1125-7865-00040
Neuwirth B, Esper J, Schweingruber FH, Winiger M (2004) Site ecological differences to the climatic forcing of spruce pointer years from the LoÈ tschental, Switzerland. Dendrochronologia 21(2):69–78. 10.1078/1125-7865-00040
Oppenheimer C Climatic, environmental and human consequences of the largest known historic eruption: Tambora volcano (Indonesia) 1815 Prog Phys Geogr: Earth Environ 2003 27 2 230 259 10.1191/0309133303pp379ra
Oppenheimer C (2003) Climatic, environmental and human consequences of the largest known historic eruption: Tambora volcano (Indonesia) 1815. Prog Phys Geogr: Earth Environ 27(2):230–259. 10.1191/0309133303pp379ra
Osborn TJ Briffa KR Jones PD Adjusting variance for sample-size in tree-ring chronologies and other regional mean timeseries Dendrochronologia 1997 15 89 99
Osborn TJ, Briffa KR, Jones PD (1997) Adjusting variance for sample-size in tree-ring chronologies and other regional mean timeseries. Dendrochronologia 15:89–99
Pretis F, Reade J, Sucarrat G (2018) Automated general-to-specific (GETS) regression modeling and indicator saturation for outliers and structural breaks. J Stat Softw 86(3):1–44. 10.18637/jss.v086.i03
Pretis F, Schneider L, Smerdon JE, Hendry DF (2017) Detecting volcanic eruptions in temperature reconstructions by designed break-indicator saturation. J Econ Surv 30:403–429. 10.1111/joes.12148
Qin C Yang B Bräuning A Sonechkin DM Huang K Regional extreme climate events on the northeastern Tibetan Plateau since AD 1450 inferred from tree rings Glob Planet Change 2011 75 143 154 10.1016/j.gloplacha.2010.10.013
Qin C, Yang B, Bräuning A, Sonechkin DM, Huang K (2011) Regional extreme climate events on the northeastern Tibetan Plateau since AD 1450 inferred from tree rings. Glob Planet Change 75:143–154. 10.1016/j.gloplacha.2010.10.013
R Core Team (2021) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/
Riechelmann DFC, Schmidhalter M, Büntgen U, Esper J (2013) Extending the high elevation larch ring width chronology from the Simplon region in the Swiss Alps over the past millennium. In: Helle G, Gärtner H, Beck W, Heinrich I, Heussner K-U, Müller A, Sanders T (eds) Proceedings of the DENDROSYMPOSIUM 2012: May 8th - 12th, 2012 in Potsdam and Eberswalde, Germany, (Scientific Technical Report ; 13/05), 11th TRACE conference (Tree Rings in Archaeology, Climatology and Ecology). Deutsches GeoForschungsZentrum GFZ, Potsdam, Eberswalde, pp 103–107
Riechelmann DFC, Hartl C, Esper J (2020) The effect of provenance of historical timber on tree-ring based temperature reconstructions in the Western Central Alps. iForest 13(1):351–359. 10.3832/ifor3412-013
Römer P Hartl C Schneider L Bräuning A Szymczak S Huneau F Lebre S Reinig F Büntgen U Esper J Reduced temperature sensitivity of maximum latewood density formation in high-elevation corsican pines under recent warming Atmosphere 2021 12 7 804 10.3390/atmos12070804
Römer P, Hartl C, Schneider L, Bräuning A, Szymczak S, Huneau F, Lebre S, Reinig F, Büntgen U, Esper J (2021) Reduced temperature sensitivity of maximum latewood density formation in high-elevation corsican pines under recent warming. Atmosphere 12(7):804. 10.3390/atmos12070804
Salzer MW Larson ER Bunn AG Hughes MK Changing climate response in near-treeline bristlecone pine with elevation and aspect Environ Res Lett 2014 9 11 114007 10.1088/1748-9326/9/11/114007
Salzer MW, Larson ER, Bunn AG, Hughes MK (2014) Changing climate response in near-treeline bristlecone pine with elevation and aspect. Environ Res Lett 9(11):114007. 10.1088/1748-9326/9/11/114007
Saulnier M Roques A Guibal F Rozenberg P Saracco G Corona C Edouard J-L Spatiotemporal heterogeneity of larch budmoth outbreaks in the French Alps over the last 500 years Can J for Res 2017 47 5 667 680 10.1139/cjfr-2016-0211
Saulnier M, Roques A, Guibal F, Rozenberg P, Saracco G, Corona C, Edouard J-L (2017) Spatiotemporal heterogeneity of larch budmoth outbreaks in the French Alps over the last 500 years. Can J for Res 47(5):667–680. 10.1139/cjfr-2016-0211
Schmidli J Schmutz C Frei C Wanner H Schär C Mesoscale precipitation variability in the region of the European Alps during the 20th century Int J Climatol 2002 22 9 1049 1074 10.1002/joc.769
Schmidli J, Schmutz C, Frei C, Wanner H, Schär C (2002) Mesoscale precipitation variability in the region of the European Alps during the 20th century. Int J Climatol 22(9):1049–1074. 10.1002/joc.769
Schneider L Smerdon JE Pretis F Hartl-Meier C Esper J A new archive of large volcanic events over the past millennium derived from reconstructed summer temperatures Environ Res Lett 2017 12 9 094005 10.1088/1748-9326/aa7a1b
Schneider L, Smerdon JE, Pretis F, Hartl-Meier C, Esper J (2017) A new archive of large volcanic events over the past millennium derived from reconstructed summer temperatures. Environ Res Lett 12(9):094005. 10.1088/1748-9326/aa7a1b
Schneider L, Smerdon JE, Büntgen U, Wilson RJS, Myglan VS, Kirdyanov AV, Esper J (2015) Revising midlatitude summer temperatures back to A.D. 600 based on a wood density network. Geophys Res Lett 42(11):4556–4562. 10.1002/2015GL063956
Schweingruber F Tree rings 1988 Bern, Stuttgart Paul Haupt
Schweingruber F (1988) Tree rings. Paul Haupt, Bern, Stuttgart
Schweingruber FH Bräker OU Schar E Temperature information from a European dendroclimatological sampling network Dendrochronologia 1987 5 9 33
Schweingruber FH, Bräker OU, Schar E (1987) Temperature information from a European dendroclimatological sampling network. Dendrochronologia 5:9–33
Schweingruber FH Bartholin T Schar E Briffa KR Radiodensitometric-dendroclimatological conifer chronologies from Lapland (Scandinavia) and the Alps (Switzerland) Boreas 1988 17 559 566 10.1111/j.1502-3885.1988.tb00569.x
Schweingruber FH, Bartholin T, Schar E, Briffa KR (1988) Radiodensitometric-dendroclimatological conifer chronologies from Lapland (Scandinavia) and the Alps (Switzerland). Boreas 17:559–566. 10.1111/j.1502-3885.1988.tb00569.x
Sevruk B Regional dependency of precipitation-altitude relationship in the Swiss Alps Clim Change 1997 36 355 369 10.1023/A:1005302626066
Sevruk B (1997) Regional dependency of precipitation-altitude relationship in the Swiss Alps. Clim Change 36:355–369. 10.1023/A:1005302626066
Shah S, Shah C (2015) Tree rings for the assessment of the potential impact of climate change on forest growth. Appl Ecol Environ Res 13(1):277–288. 10.15666/aeer/1301_277288
Sigl M, Toohey M, McConnell JR, Cole-Dai J, Severi M (2021) HolVol: reconstructed volcanic stratospheric sulfur injections and aerosol optical depth for the Holocene (9500 BCE to 1900 CE). PANGAEA. 10.1594/PANGAEA.928646
St. George S, Esper J (2019) Concord and discord among Northern Hemisphere paleotemperature reconstructions from tree rings. Quat Sci Rev 203:278–281. 10.1016/j.quascirev.2018.11.013
Stoffel M Khodri M Corona C Guillet S Poulain V Bekki S Guiot J Luckman BH Oppenheimer C Lebas N Beniston M Masson-Delmotte V Estimates of volcanic-induced cooling in the Northern Hemisphere over the past 1,500 years Nat Geosci 2015 8 10 784 788 10.1038/ngeo2526
Stoffel M, Khodri M, Corona C, Guillet S, Poulain V, Bekki S, Guiot J, Luckman BH, Oppenheimer C, Lebas N, Beniston M, Masson-Delmotte V (2015) Estimates of volcanic-induced cooling in the Northern Hemisphere over the past 1,500 years. Nat Geosci 8(10):784–788. 10.1038/ngeo2526
Stothers RB The Great Tambora Eruption in 1815 and its aftermath Science 1984 224 1191 1198 10.1126/science.224.4654.1191 17819476
Stothers RB (1984) The Great Tambora Eruption in 1815 and its aftermath. Science 224:1191–119817819476
Torbenson MCA Büntgen U Esper J Urban O Balek J Reinig F Krusic PJ Martinez del Castillo E Brázdil R Semerádová D Štěpánek P Pernicová N Kolář T Rybníček M Koňasová E Arbelaez J Trnka M Central European Agroclimate over the past 2000 years J Clim 2023 36 13 4429 4441 10.1175/JCLI-D-22-0831.1
Torbenson MCA, Büntgen U, Esper J, Urban O, Balek J, Reinig F, Krusic PJ, Martinez del Castillo E, Brázdil R, Semerádová D, Štěpánek P, Pernicová N, Kolář T, Rybníček M, Koňasová E, Arbelaez J, Trnka M (2023) Central European Agroclimate over the past 2000 years. J Clim 36(13):4429–4441. 10.1175/JCLI-D-22-0831.1
Trachsel M Kamenik C Grosjean M McCarroll D Moberg A Brázdil R Büntgen U Dobrovolný P Esper J Frank DC Friedrich M Glaser R Larocque-Tobler I Nicolussi K Riemann D Multi-archive summer temperature reconstruction for the European Alps, AD 1053–1996 Quat Sci Rev 2012 46 66 79 10.1016/j.quascirev.2012.04.021
Trachsel M, Kamenik C, Grosjean M, McCarroll D, Moberg A, Brázdil R, Büntgen U, Dobrovolný P, Esper J, Frank DC, Friedrich M, Glaser R, Larocque-Tobler I, Nicolussi K, Riemann D (2012) Multi-archive summer temperature reconstruction for the European Alps, AD 1053–1996. Quat Sci Rev 46:66–79. 10.1016/j.quascirev.2012.04.021
Wang J Yang B Fang M Wang Z Liu J Kang S Synchronization of summer peak temperatures in the medieval climate anomaly and little ice age across the Northern Hemisphere varies with space and time scales Clim Dyn 2023 60 11–12 3455 3470 10.1007/s00382-022-06524-6
Wang J, Yang B, Fang M, Wang Z, Liu J, Kang S (2023) Synchronization of summer peak temperatures in the medieval climate anomaly and little ice age across the Northern Hemisphere varies with space and time scales. Clim Dyn 60(11–12):3455–3470. 10.1007/s00382-022-06524-6
Wanner H Pfister C Neukom R The variable European Little Ice Age Quat Sci Rev 2022 287 107531 10.1016/j.quascirev.2022.107531
Wanner H, Pfister C, Neukom R (2022) The variable European Little Ice Age. Quat Sci Rev 287:107531. 10.1016/j.quascirev.2022.107531
Wilmking M D’Arrigo RD Jacoby GC Juday GP Increased temperature sensitivity and divergent growth trends in circumpolar boreal forests Geophys Res Lett 2005 32 15 L15715 10.1029/2005GL023331
Wilmking M, D’Arrigo RD, Jacoby GC, Juday GP (2005) Increased temperature sensitivity and divergent growth trends in circumpolar boreal forests. Geophys Res Lett 32(15):L15715. 10.1029/2005GL023331
Wilson RJS D’Arrigo R Buckley B Büntgen U Esper J Frank D Luckman B Payette S Vose R Youngblut D A matter of divergence: tracking recent warming at hemispheric scales using tree ring data J Geophys Res 2007 112 D17 D17103 10.1029/2006JD008318
Wilson RJS, D’Arrigo R, Buckley B, Büntgen U, Esper J, Frank D, Luckman B, Payette S, Vose R, Youngblut D (2007) A matter of divergence: tracking recent warming at hemispheric scales using tree ring data. J Geophys Res 112(D17):D17103. 10.1029/2006JD008318
Wilson RJS Anchukaitis KJ Briffa KR Büntgen U Cook ER D’Arrigo R Davi N Esper J Frank D Gunnarson BE Hegerl GC Helama S Klesse S Krusic PJ Linderholm HW Myglan V Osborn TJ Rydval M Schneider L Schurer A Wiles G Zhang P Zorita E Last millennium northern hemisphere summer temperatures from tree rings: part I: the long term context Quat Sci Rev 2016 134 1 18 10.1016/j.quascirev.2015.12.005
Wilson RJS, Anchukaitis KJ, Briffa KR, Büntgen U, Cook ER, D’Arrigo R, Davi N, Esper J, Frank D, Gunnarson BE, Hegerl GC, Helama S, Klesse S, Krusic PJ, Linderholm HW, Myglan V, Osborn TJ, Rydval M, Schneider L, Schurer A, Wiles G, Zhang P, Zorita E (2016) Last millennium northern hemisphere summer temperatures from tree rings: part I: the long term context. Quat Sci Rev 134:1–18. 10.1016/j.quascirev.2015.12.005
Zang C Biondi F treeclim: an R package for the numerical calibration of proxy-climate relationships Ecography 2015 38 4 431 436 10.1111/ecog.01335
Zang C, Biondi F (2015) treeclim: an R package for the numerical calibration of proxy-climate relationships. Ecography 38(4):431–436. 10.1111/ecog.01335
Zhang P Björklund J Linderholm HW The influence of elevational differences in absolute maximum density values on regional climate reconstructions Trees 2015 29 4 1259 1271 10.1007/s00468-015-1205-4
Zhang P, Björklund J, Linderholm HW (2015) The influence of elevational differences in absolute maximum density values on regional climate reconstructions. Trees 29(4):1259–1271. 10.1007/s00468-015-1205-4
