Tree-ring width and glacier mass balance are two highly sensitive climatic proxies which are often used as indicators of biological and geophysical changes in high-altitude ecosystems. Tree-ring data have been widely used to reconstruct past temperatures and also to reconstruct past glacier mass balance. Here we show that tree-ring chronologies from a high-altitude Pinus cembra L. dendroclimatic network and glaciers from the same region in the European Alps have non-stationary responses to air temperature, and have also been responding non-proportionally to temperature extremes in recent decades. Both ring-width chronologies and the mass-balance series of some glaciers from the same region have shown an increasing sensitivity to summer (JJA) temperatures. Our results demonstrate that the sensitivity to climate of tree-ring chronologies and glacier mass balance may change over time and has been increasing in recent decades, posing some limitations to tree-ring-based glacier mass-balance reconstruction. Moreover, we found these reconstructions in the European Alps are more reliable for large rather than for small glaciers, and may not be able to reveal years of extreme ablation that could have occurred in the past.
The impacts of climate change on physical and biological components of ecosystems are not homogeneous worldwide, and some regions may be more disproportionately impacted than others. Changes in temperature and precipitation regimes (e.g., Jones et al., 1999; Mokhov et al., 2006) and a higher frequency of extreme meteorological events (e.g., Beniston, 2007) have been recorded in many regions worldwide (Intergovernmental Panel on Climate Change, 2007). Changes in the physical and biological components of ecosystems due to climate change are especially evident in those regions where climate factors greatly limit the productivity of local ecosystems (Hughes, 2000; Huber et al., 2005; Barnett et al., 2005). Among the indicators used to understand the ongoing climatic changes, glacier mass balance and tree rings are two of the more sensitive climatic proxies (e.g., Oerlemans, 2001; Hughes, 2002). Glacier net mass balance is the variation in ice mass, expressed in m w.eq. (water equivalent), and is usually calculated for mid- and high-latitude glaciers as the sum of the winter (accumulation) and summer (ablation) balance at the end of the hydrological year. Net mass balance is influenced by local weather conditions, mostly (summer) temperature and (snow) precipitation. As such it is considered an indicator of climate that reflects large-scale atmospheric circulation patterns (Bitz and Battisti, 1999; McCabe et al., 2000; Diolaiuti et al., 2002). In the interior portion of the European Alps, where climatic conditions are moderately continental, mass balance is mostly influenced by the summer temperature regime (Dyurgerov, 2006; Haeberli et al., 2007). How glaciers respond to climatic changes depends, however, on many factors, which are mostly related to the glaciers' size and geometry, but also to their surface characteristics influencing their albedo. For example, debris-covered glaciers are less sensitive to climatic inputs because the rocky layer influences the surface energy exchanges (Nakawo and Young, 1982; Mihalcea et al., 2006). Changes in mass, and consequently in equilibrium line altitude (ELA) of debris-free glaciers, on the other hand, are directly driven by climate (Haeberli, 1990; Pelfini and Smiraglia, 1997; Dyurgerov, 2002). As a result, the fluctuation of the glacier terminus, induced by changes in ice mass, have different response times (Oerlemans, 2001), varying from a few years for small, steep mountain glaciers up to several decades or even more than a century for large, flat valley glaciers (Haeberli and Hoelzle, 1995; Hoelzle et al., 2003).
Climate changes can directly impact animals and plants. Their adaptive capacities are known to be species-specific and depend on many factors often directly connected to their physiology (Parmesan, 1996; Hughes, 2000). Trees living in climatically extreme environments record information on the climatic conditions of a region in their annual radial growth (Fritts, 1976; Hughes, 2002). At high altitudes (and high latitudes), trees growing close to the treeline are usually temperature-limited and the derived tree-ring chronologies may be used for proxy-based reconstructions of past temperatures (e.g., Briffa et al., 1990; Esper et al., 2002). Long tree-ring chronologies may allow for millennial-long reconstructions of climatic conditions during the growth period of vegetation, providing reliable information on interannual variability and long-term trends.
Both glacier mass balance and tree rings are very sensitive to climate (e.g., Theurillat and Guisan, 2001; Hughes, 2002; Beniston, 2003; Huber et al., 2005; Oerlemans, 2005; Haeberli et al., 2007), and some efforts have been made to reconstruct past glacier mass balance on the basis of dendrochronological data. The main task in such proxy-based reconstructions is the assessment of the past natural variability of a climatic parameter for time periods prior to instrumental measurements. Successful tree-ring-based reconstructions of past glacier mass balance have been produced for the European Alps and the Western Cordillera of North America. In the European Alps, the best results have been obtained using tree-ring data of Pinus cembra L. (Nicolussi, 1994; Leonelli et al, 2008). Based on a regression model, past mass-balance estimates for the Hintereis Glacier cover about five centuries (Nicolussi and Patzelt, 1996). In the Western Cordillera and the Rocky Mountains, tree-ring based reconstructions have been performed to estimate the past mass balance of four glaciers in a region extending over ∼300 km (Lewis and Smith, 2004), as well as the summer and winter balances of several glaciers (Larocque and Smith, 2005) and of the Peyto Glacier (Watson and Luckman, 2004).
All proxy-based reconstructions are based on the uniformitarian principle, which relies on the stability of relationships between climatic records and proxies. The stability of relationships is established over a common period of record overlap before full reconstructions of past climate are generated. This implies that the sensitivity to climate of both glacier mass balance and tree rings is assumed to be stable over time. However, recent studies, have shown that tree-ring sensitivity to climate may vary over time (Büntgen et al., 2006a; Carrer and Urbinati, 2006; Carrer et al., 2007; Leonelli et al., 2009) and that some divergence between tree-ring chronologies and some climatic records is a global issue (D'Arrigo et al., 2008). This has implications for all tree-ring-based reconstructions. Another problem with long-term reconstructions of a glacier mass balance is that statistical relationships found over the calibration period may be biased since glaciers of the late-20th century are considerably smaller than they were few hundred years earlier, under different climate conditions (e.g., during the Little Ice Age) (Meier at al., 2003). The corresponding effect could indeed amount to a factor of two for the time period since the maximum glacier extent during the Little Ice Age (around 1850 in the European Alps; Nemec et al., 2009).
In this paper, on a yearly resolved base, we test the temporal stability of glacier and tree-ring sensitivity to climate over a common period, with a focus on the potential for tree-ring widths to reconstruct glacier mass balance in the European Alps. These two climatic proxies are compared by analyzing net mass balances from debris-free glaciers of different size, and tree-ring chronologies from a temperature-sensitive high-altitude P. cembra network already used in a previous dendroclimatic study (Leonelli et al., 2009). Choosing the same region allowed us to investigate how the same climatic conditions affected these two climatic proxies in an Alpine environment. Particular attention was given to the characterization of the climatic signal recorded in the tree-ring chronologies that show the strongest and weakest correlations with the mass-balance series. This was done to define an approach to select tree-ring chronologies a priori for a mass-balance reconstruction in the region.
Data and Methods
The aim of our analysis is to test the temporal stability of tree-ring and glacier mass-balance responses to the same climatic input. The tree-ring data came from a high-altitude P. cembra tree-ring network composed of 13 chronologies from six valleys of the Ortles-Cevedale Group (46°28′N, 10°34′E) in the central European Alps (Fig. 1). These were compared with three glacier mass-balance series from the same region. The central European Alps are characterized by several mountain groups with many peaks over 3000 m a.s.l. The Ortles-Cevedale Group, on the Italian side of the Alpine chain, covers an area extending over 30 km for which a considerable wealth of glaciological, dendrochronological, and meteorological data are available. The climate in the region is moderately continental, with the highest temperatures and precipitation during summer.
We considered the net mass-balance series from a valley glacier in the Oetztal Alps, on the Austrian side of the Alpine chain (the Hintereis Glacier, 7.7 km long, with a surface of about 9.5 km2, mean elevation 3020 m a.s.l., aspect N/NE; National Snow and Ice Data Center, 2007), and the series of a small mountain glacier in the Ortles-Cevedale Group, on the Italian side (the Careser Glacier, with a surface of 4.7 km2, mean elevation 3090 m a.s.l., aspect S). Additionally, we calculated a regional net mass-balance series for three mountain glaciers from the same mountain group (Fontana Bianca: 0.7 km2 at 3158 m a.s.l., E; Sforzellina: 0.39 km2 at 2925 m a.s.l., NW; and Careser: see above for details) after standardizing each series (z-scores) and calculating the mean departures. The three mass-balance series of the Hintereis (HEF; 1953–2003), of the Careser (CAR; 1967–2003), and of the regional group (COMP; 1967–2003) were then standardized, resulting in three mass-balance index (MBI) series (Fig. 2). These were then used in further analyses. The selected series are from debris-free glaciers and are from glaciers at a mean elevation of 3060 m a.s.l. During the common period 1967–2003, they show a mean mass-balance value of −0.68 m w.eq. and a mean autocorrelation of 0.38 (Table 1). All net mass-balance series were obtained from the World Glacier Monitoring Service (W.G.M.S.) (1959–2000 and online values: http://www.geo.unizh.ch/wgms/).
The main characteristics of the site and the dendrochronological/glaciological series. AR1 = first order autocorrelation referred to the regional curve (RC) and mass-balance index (MBI) series. The first letter of the site chronology codes corresponds to the first letter of the valley names (F = Forni; M = Martello; P = Peio; T = Trafoi; U = Ultimo; Z = Zebrù) and the middle letter identifies the site in the valley. The last letter corresponds to the site slope aspect. HEF = Hintereis MBI; CAR = Careser MBI; COMP = Composed MBI.
The P. cembra ring-width chronologies are from 13 high-altitude sites in 6 valleys of the Ortles-Cevedale Group. We prepared the dendrochronological data by removing the non-climatic growth-trend due to tree age or stand disturbances in the tree-ring chronologies (Fritts, 1976), preserving low-frequency variability in the resulting chronologies. For each site, raw measurements were transformed into ring-width index series by dividing each value by the one predicted by the regional curve (RC; Esper et al., 2003). The RC was calculated by aligning all the measurement series according to their cambial age and smoothing the mean values with a cubic spline 10% of the series length (Büntgen et al., 2006b). A biweight robust mean was applied to all the series, finally resulting in the site chronologies. A mean chronology was also prepared by averaging all 13 site chronologies. Moreover, six subsets of the tree-ring chronologies were created by averaging the four chronologies with the strongest (Schr_HEF, Schr_CAR, and Schr_COMP) or the weakest (Wchr_HEF, Wchr_CAR and Wchr_COMP) absolute correlation values with the three different MBI series (Fig. 2 and Table 2). The tree-ring chronologies are from sites with mean elevations ranging from 2100 to 2250 m (mean 2190), mainly from N-facing slopes. The mean ring widths range from 0.7 to 1.6 mm (mean 1.2) and the series have a first-order autocorrelation ranging from 0.7 to 0.9 (mean 0.8) (Table 1).
Constituents of the Schr and Wchr chronologies. The two groups of chronologies, Schr and Wchr, are defined as subsets of the site chronologies presenting the four strongest (Schr) and the four weakest (Wchr) correlation values versus the three MBI considered. Only the chronologies of the Schr group were used to compare the dendrochronological series with the MBI series. The groupings for the CAR and the COMP MBI series are identical.
Meteorological gridded data (monthly and seasonal anomalies) from 1865 to 2003 for the region comprising the study area were taken from the HISTALP database (Auer et al., 2007), grid point 46°N, 10°E. Anomalies refer to the 20th century means (1901–2000) and were derived only from homogenized data, both for the seasonal and monthly temperature series (high-altitude record, >1400 m a.s.l.). July and summer temperatures over the region show an increasing trend over the last century and recent decades (Fig. 2), whereas precipitation does not show any particular trend for either monthly and seasonal totals (not shown).
The selected glacier and tree-ring series were compared with the monthly and seasonal climatic variables during the same periods using classical and moving correlation analyses. Monthly and seasonal climatic analyses were performed on the Schr and Wchr series to check the sensitivity of the tree-ring growth to climate in the two groups of chronologies. Tree-ring chronologies were analyzed for the period 1968–2003, selecting the temperature and precipitation variables from October of the year prior to growth to September of the year of growth, as these generally represent the hydrological year (Leonelli et al., 2008) and from spring (MAM-1) of the year prior to growth up to autumn (SON) of the year of growth. The months prior to growth are commonly selected in dendroclimatic analysis because annual growth is influenced by important biological lags (e.g., needle retention, bud formation, carbohydrate storage) which are the reason for autocorrelation in the tree-ring series (Fritts, 1976). The same method was used to evaluate the mass-balance series relationships to climate.
For the moving correlation analysis, we analyzed 30-year periods progressively shifted by one year from 1953 to 1982 (1967 to 1993 for CAR and COMP series) up to 1974–2003. We calculated the bootstrapped Pearson's correlation coefficients (1000 iterations) between both the dendrochronological and the glaciological series, and summer (JJA) temperatures, spring (MAM) and summer (JJA) precipitation. Seasonal precipitation variables were chosen according to the statistically significant results obtained from the monthly analysis of mass balance and of the chronologies, respectively. Correlation coefficients calculated over 30-year periods may be highly variable. However, a period of 30 years is a standard commonly used in climatology to assess climate normals and it also adapts well to the selected mass-balance series lengths. The significance of the correlation coefficients for the overlapping intervals was assessed using the Bonferroni correction method, which adjusts confidence levels when multiple statistical tests are evaluated simultaneously (Snedecor and Cochran, 1989). Statistical comparisons between tree rings and glacier mass balance were performed using only the three Schr chronologies (Schr_HEF, _CAR, and _COMP).
A linear regression of the MBI series on the Schr series was also performed for the longest time period available and on the common period 1967–2003. To investigate the influence of non-proportional impacts arising from extreme temperature, a second linear regression was calculated after excluding the year 2003, which is an outlier in the temperature record due to the an extremely hot summer in Europe.
The correlations between all the dendrochronological and glaciological series, calculated for the common period 1967–2003 showed that the glacier series were most strongly correlated to one another (r ≥ 0.90, p < 0.001; Table 3). The mean value of the correlation coefficients calculated in all the comparisons between the pairs of tree-ring chronologies is 0.56. Many of the coefficients are significant at the p < 0.001 level (Table 3). The FCN chronology correlated least with the other chronologies, with up to three non-significant correlations.
Correlation coefficients calculated over the common period 1967–2003 between the 13 RC site chronologies and between the three MBI series selected for this study. Significant values at the p < 0.001 level are written in bold; non-significant values (p > 0.05) are written in italics. Abbreviations as in Table 1.
The climatic analysis on the Schr and Wchr series, performed for the period 1968–2003, revealed that the chronologies correlating better with the MBI series (the Schr) are those containing the strongest J and JJA temperatures signals (Figs. 3a and 3b, respectively). These chronologies also showed that JJA-1 temperature had a positive influence on tree-growth and that March precipitation (springtime snow) had a negative influence on tree-ring width. The Wchr chronologies had, in contrast, a negative correlation with the previous winter (D-1JF) and February temperatures. In all the chronologies the previous autumn (SON-1) temperature had a positive influence on tree-ring width and summer precipitation (JJA P) a negative influence.
The three glaciological series correlated with the climatic inputs in an almost opposite way to the tree-ring chronologies (especially with respect to the seasonal variables). The summer (JJA) temperature had the strongest influence on net mass balance and all summer months were significantly and negatively correlated (Figs. 3c and 3d). Positive correlations were found with previous Winter (D-1JF), Spring (MAM), and summer (JJA) precipitation.
The sensitivity of the Pinus cembra mean chronology to JJA temperatures varied over time, and the correlation coefficient values calculated for 30-year periods tend to increase (Fig. 4a). The values up to r = 0.64 for the 1973–2002 period are evidence of a recent trend towards a greater sensitivity. A trend towards a higher sensitivity to JJA precipitation was also found for the mean chronology in the recent period (Fig. 4a).
In general, the Hintereis series (HEF) showed stronger correlation values with summer temperatures than the mass-balance series of the smaller glaciers. However, the correlations were always significant (p < 0.001) for all the series and in all the intervals considered (Fig. 4b). Correlation coefficients reveal a generally increasing sensitivity of the glaciers to summer temperatures (this is especially evident for the Hintereis). A marked trend of a higher influence of MAM precipitation on glacier mass balance is especially evident for the Hintereis Glacier since the period 1967–1996 up to the latest 30-year periods analyzed (Fig. 4b).
Correlations between the three MBI series and the three Schr tree-ring chronologies (Schr-HEF, -CAR, and -COMP) in general reveal a convergence between these two climatic proxies in recent decades, especially since the period 1965–1994 (Fig. 5). The correlations were generally stronger for the large glacier and always significant, except in the period 1968–1997 for the COMP series of mass balance. Since the period 1970–1999, the Careser series has also shown stronger correlations with tree-ring data, nearly reaching the Hintereis values.
A linear regression of the MBI series HEF and CAR on the respective Schr series was performed for the period 1953–2003 and for the common period 1967–2003 (Fig. 6). Regressions were calculated using complete series and series without the year 2003, which is an outlier in the mass-balance series. The regression of MBI series on the tree-ring chronologies consistently showed negative slopes. Tree-ring chronologies could explain higher percentages of variations in the mass balance of the Hintereis Glacier (21%, and 26% in the 1967–2003 period) than of those of the Careser Glacier (16%). After removing year 2003 from the analysis, higher proportions of the mass-balance variations were explained by the tree-ring chronologies. Tree-ring chronologies explained up to 25% (33% in the 1967–2003 period) of the variations in the Hintereis series and up to 20% in the Careser series. The COMP series always showed the lowest values for the coefficient of determination and non-significant regressions (data not shown).
Recent research in the European Alps indicates that tree-ring sensitivity to climate may change over time. Büntgen et al. (2006a) found a variation in sensitivity to climate in a high-altitude Picea abies tree-ring network, Carrer and Urbinati (2006) in a Larix decidua network, and Carrer et al. (2007) in a P. cembra network. The climate sensitivity of the P. cembra dendroclimatic network used in this study also appeared to change over time, depending partly on site slope aspect (Leonelli et al., 2009). Chronologies from N-facing sites showed almost stable relationships with JJA temperature over time, whereas chronologies from S- and W-facing sites showed a dramatic decrease in sensitivity. In the present study the moving correlation analysis revealed an increasing sensitivity of the mean chronology to JJA temperatures (Fig. 4a). However, the analysis was performed for 30-year periods, which is half the length of the time period used in Leonelli et al. (2009). Using this time scale may induce higher variability in the correlation coefficients than using longer time periods.
The synchronicity and the similar patterns of the mass-balance series reflect a regional climatic signal. In general, glaciers within the same region typically present similar patterns in the ice-mass variation, and the similarity of glaciological series tends to decrease with an increasing distance between glaciers (Cogley and Adams, 1998; Diolaiuti et al., 2002; Letréguilly and Reynaud, 1990). However, the Hintereis mass-balance series generally showed a stronger correlation with the JJA temperature, revealing a different response to climate from those of the other two series from smaller glaciers (Fig. 3d).
Our findings also indicate that the same climatic conditions over a region may result in non-proportional variations in a glacier mass balance and in tree-ring growth, with the years presenting extreme variations in the mass balance being the least likely to be correctly estimated by these relationships. A clear example of this is summer 2003, when a heat wave occurred over Europe and the Alps. Glaciers then underwent extreme melting, whereas the tree-ring width increment at high altitudes was not particularly enhanced (Fig. 2). Glaciers in the Alps seem to be increasingly showing signs of downwasting and collapse due to rapid warming rather than “retreat” (Haeberli et al., 2007). The influence of local factors like topographic features and the increasing areas covered by rock outcrops in the glacier bodies largely depends on glacier size: under similar climatic conditions, glacier fragmentation and melting rates are especially high for small and thin glaciers (Paul et al., 2004). This could also explain the weaker correlations of the mass balance of small glaciers with summer temperature (Figs. 3c and 3d). Smaller glaciers, like the ones selected for this study, are the most likely to disappear first from the alpine landscape. In fact, some of them, including the Careser Glacier, are already close to full disintegration (Paul et al., 2004; Carturan and Seppi, 2007). Moreover, positive feedbacks accelerate glacier melt, such as a decrease in albedo due to firn area loss and dust/soot deposition on ice/firn surfaces (Oerlemans et al., 2009; Paul et al., 2005).
A convergence between glacier mass change and temperature records has also been reported on the global scale (Dyurgerov, 2006). The mass balance of glaciers may be considered an integrator of large-scale atmospheric circulation and surface climatic factors, and may thus be interpreted as the first indicator of an undelayed climate signal (Solomina et al., 2008). The sensitivity of a glacier's mass balance has a regional/climatic component (continentality, mass-balance gradient) as well as an individual/topographic component (hypsography: distribution of glacier area with altitude). The tongues of large valley glaciers may be especially sensitive to temperature change since a large portion of their bodies typically extend to low altitudes, where summer temperatures generally result in less snowfall and more melting. This could explain the higher sensitivity of glaciers, especially for the Hintereis Glacier, to spring (MAM) snowfall in recent decades (Fig. 4b).
The general convergence noticed between the three MBI series and the three tree-ring chronologies (Schr-HEF, -CAR, and -COMP; Fig. 5) is probably due to the increasing sensitivity of tree rings and glaciers to summer temperatures (Figs. 4a and 4b). The variation in the strength of the linear relationships over time between the MBI series of both large and small glaciers, and the tree-ring chronologies poses some limitations on the glacier mass-balance reconstructions. The reconstruction of past mass balances of the small glaciers selected for this study appears to be especially problematic since the series are short and correlations are generally weaker and also more variable over time. Non-stationary relationships with climatic elements and between variations in climatic proxies (both physical and biological ones) have significant implications for all proxy-based reconstructions.
Comparisons between the two climatic proxies, mass-balance series and tree-ring chronologies, have led to successful mass-balance reconstructions in Europe and North America (Nicolussi and Patzelt, 1996; Lewis and Smith, 2004; Larocque and Smith, 2005; Watson and Luckman, 2004). Upon close examination of our dendrochronological data for the central European Alps, however, we found stronger correlations for the larger Hintereis Glacier series, whereas relationships to smaller glaciers closer to the dendroclimatic network were generally weaker. The tree-ring chronologies that correlated better with the MBI series (Schr) are those containing a stronger July and JJA temperature signal and a negative influence of March snow (Figs. 3a and 3b). Those chronologies presenting weaker correlations with the MBI (Wchr) also had weaker July and JJA temperature signals.
The regression of the MBI series on tree-ring chronologies always showed negative slopes, but the tree-ring chronologies always presented higher percentage of mass-balance variations of the Hintereis glacier than of the small glaciers. After removing the year 2003, considerably higher proportions of mass-balance variations were explained by tree-ring chronologies. These findings emphasize how the same climatic conditions drive non-proportional variations in glaciers' mass variations and tree-ring growth, especially in the years when extreme meteorological conditions induce strong variations. The different responses to climate could be related to the fact that these two climatic proxies respond to common but also different climate variables and to the fact that high-altitude trees have threshold temperatures above which they may not respond linearly (D'Arrigo et al., 2004). In particular, we found that tree rings mainly responded to July temperature and glacier mass balance to JJA temperature (Figs. 3a and 3c). Also the influence of precipitation variables slightly differs, with tree rings being mainly influenced by JJA precipitation and glaciers by January to May precipitation (Figs. 3b and 3d). Since precipitation is a controlling factor modulating both tree-ring growth and glacier mass balance, another possible improvement for mass-balance reconstruction in this area therefore could be the inclusion of precipitation-sensitive chronologies and the consideration of winter and summer mass balances (Watson and Luckman, 2004).
Both glacier mass balance and tree rings from high-altitude sites show a non-stationary sensitivity to climate over time. The recent convergence between the glacier MBI series and ring-width chronologies was probably driven by the generally increasing sensitivity of the glacier mass balance to JJA temperatures and of the tree rings especially to J temperatures for the selected chronologies. This convergence between the two proxies may be mostly ascribed to the modern temperature regime over the region that shows a marked increasing trend both for JJA and J temperature. In the first portion of the considered period of analysis the two proxies were, instead, less correlated: changes in the climatic inputs occurred over different months and the generally cooler climatic regime may have caused different responses in the two proxies over time. This fact may pose some limitations on tree-ring-based glacier mass-balance reconstructions in the region studied, which may also apply in other mountain regions. The mass balance of the large Hintereis valley glacier generally correlated better with summer (JJA) temperatures than did those of small glaciers. These results indicate that tree-ring-based glacier mass-balance reconstructions in the European Alps are probably more reliable for large than for small glaciers.
We found that summers with extremely warm conditions, which are likely to become more frequent in the European Alps (Beniston, 2007), are inducing strong mass variation in glaciers but they are not correspondingly enhancing ring widths. While tree rings may provide a fairly reliable mass-balance reconstruction of large glaciers, such a reconstruction cannot reveal years of strong ablation that could have occurred in the past. For the Central Alps we suggest that first tree-ring data should be used to reconstruct the mass balance of large valley glaciers. Non-stationary relationships with climatic elements over time also pose some problems for tree-ring based reconstructions, because typically only short calibration periods are available. Given our results, before applying the uniformitarian principle to reconstruct past glacier mass-balance changes, we recommend rigorous testing of both the physical and the biological climatic proxies over common and maximum periods to determine what climatic variable they mostly respond to and to find out how stationary the influence of climate is over time.
This research was funded by the MIUR-COFIN 2005 and PRIN 2009 projects. The authors thank the reviewers for their useful comments and help, Dr. Reinhard Böhm (ZAMG, Vienna, Austria) for the HISTALP database, Dr. Silvia Dingwall for revising the English, Josée Frenette (Université du Québec AT) for fieldwork assistance, and the Stelvio National Park for tree-sampling permits.