Tree-ring δ18O from an Alpine catchment reveals changes in glacier stream water inputs between 1980 and 2010

ABSTRACT The tree-ring stable carbon and oxygen isotope chronologies from two forest sites located in the Forni Glacier forefield (Italy)—one along the glacier stream (GL) and the other toward the valley slope (SL)—were analyzed with the aim of disentangling the precipitation and glacier meltwater inputs in source water δ18O, as reflected by the tree-ring cellulose δ18O. The cellulose δ18O from the GL trees has a negative correlation with winter and summer temperatures, whereas the cellulose δ18O from the SL trees has a positive correlation with precipitation δ18O. The isotopic signature of the source water at the GL site is also influenced by waters of glacial origin, as confirmed by the 18O-depleted glacier meltwater inputs (GMWI_δ18O) estimated by means of an isotope model. The GMWI_δ18O values are consistent with the mean difference measured between the δ18O in the glacier stream and in the precipitation and the winter and summer temperature explains up to 37 percent of the GMWI_δ18O variance. Our results show an increasing influence of glacier meltwater throughout the past decade for the GL site. Our analysis opens new opportunities to reconstruct changes in water regimes of the glacier streams by means of the tree-ring cellulose δ18O.


Introduction
Stable carbon and stable oxygen isotopes in relation to climate and climate-related dynamics have often been applied both in forest ecology, for assessing climatechange impacts (e.g., Battipaglia et al. 2013;Wieser et al. 2016), and in paleoclimate studies, for reconstructing past climatic conditions on a yearly basis (e.g., Kress et al. 2010;Sidorova et al. 2013;Hartl-Meier et al. 2015). The impacts of a globally relatively fast-changing climate involve several forest environments, and in their upper altitudinal limits toward the treeline, forests are probably facing the most exceptional changes. Among these changes, longer and warmer growing seasons, higher concentrations of carbon dioxide in the atmosphere, higher N-deposition, and the transformation of forest use (in the European Alps mainly the result of land abandonment) represent the most important factors acting at both the regional and global scales.
One of the most widespread impacts of the increasing mean air temperature at the global scale is the cryosphere shrinkage, which is especially evidenced by the ongoing prevalent trends of increasing ice-melting rates and negative glacier mass balances in high mountain environments (Fountain et al. 2012;Huss et al. 2017), the extensive glacier front retreats, and the increasing frequency of glacier extinctions. Additionally, the annual number of days with snow cover is generally decreasing at high latitudes (e.g., Dye and Tucker 2003) and at high altitudes in the Alps (e.g., Hantel and Lucia-Maria 2007;Marty et al. 2017). Thus, this decrease in snow cover also contributes to changes in the hydrological cycle in mountain catchments in terms of thaw anticipation, high glacier ablation rates, and higher water discharges in the glacier streams. Changes in the water discharge of Alpine streams and rivers because of changes in glacier meltwaters may indeed drive important changes in mountain ecosystems (Slemmons, Saros, and Simon 2013) and potentially impact tree growth in mountain forests. Because of glacier and permafrost melt, a higher discharge in glacier streams is already ongoing; however, an extensive glacier runoff reduction is expected in the mid-to long term because the glacier mass is dramatically reducing and because other processes, such as sublimation and evaporation, will probably have a higher impact (Kaser 2001). Likely, small glaciers will disappear first from the Alpine landscape. Larger glaciers, instead, will still likely show increasing water discharge amounts during the summer ablation period because of the presence of large amounts of ice. Medium-sized glaciers will decrease their water discharge in the near future; for example, the Corbassière Glacier in Valais will increase its water discharge only until 2020, and a subsequent constant declining trend in glacier stocks is foreseen (Terrier et al. 2011). In addition to the climaticecological dynamics involving the upper portions of the forests and the treelines, the forests that colonize and grow in glacier forefields will face exceptional changes after the retreat of the glacier tongues. Trees in glacier forefields typically grow on initial and very young soils formed on glacial deposits and till. These recently deglaciated areas at elevations below the plant species' limit and below the treeline-which is the case for most of the large valley glaciers in the Alps, including our study site -typically undergo primary successions and forest invasion (Huggett 1998;Walker et al. 2010;Garbarino et al. 2010;D'Amico et al. 2015;D'Agata et al. 2019).
The understanding of the mid-to long-term spatiotemporal hydrological changes that occur in glacier streams that are still fed by glaciers is therefore a critical issue for assessing both changes in water resources in the Alpine environment and the impacts of glacier meltwater changes on forest ecosystems. Treering stable isotopes are highly sensitive proxies for reconstructing and dating past and ongoing climatic or environmental changes (e.g., Battipaglia et al. 2008;Treydte et al. 2009;Daux et al. 2011). Tree photosynthetic capacity and stomatal conductance mainly drive the cellulose δ 13 C variability (Stuvier and Braziunas 1987;Farquhar, Ehleringer, and Hubick 1989), and even at moist sites, significant relationships are found with environmental parameters, such as summer irradiance, temperature, and drought (e.g., McCarroll and Loader 2004;Seibt et al. 2008;Kress et al. 2014). The δ 13 C is further directly impacted by the fossil-fuel emission rates in the atmosphere at least since the nineteenth century, which can be well corrected for (Treydte et al. 2009). Several other factors, such as the internal mixture and storage of carbon with different turnover at the leaf or the stem levels or the external combined effects of environmental factors, may alter the climatic signals (Seibt et al. 2008;Leonelli et al. 2012).
The analysis of the tree-ring stable isotopes from glacial environments, in the presence of young and shallow soils, has the potential to make mid-to longterm spatiotemporal hydrological reconstructions in the vicinity of glacier streams, especially focusing on cellulose δ 18 O. In fact, trees growing in the glacial environments of the Alps are highly sensitive to the δ 18 O changes in the soil waters induced by glacier meltwater or by snowmelt and may be used to date and reconstruct past glacier runoff events  or to reconstruct past winter and summer precipitation δ 18 O, respectively (Leonelli et al. 2017). An important assumption for the present study is that in the study area, the cellulose δ 18 O may be influenced by the δ 18 O of winter snow and summer precipitation (Leonelli et al. 2017). In the intermountain west of the United States the use of hydrogen isotope ratios to track the source waters of trees has underlined differences between young and mature trees growing along or far from a stream (Dawson and Ehleringer 1991): young trees were fed by precipitation (far from the stream) or by stream waters (along the stream), whereas mature trees along the stream used little or none of the stream water, being mainly fed by waters from deeper soil strata, thus underlining a possible change in the source water used by trees while growing and developing a more mature root system. Combined analyses of the stable carbon and oxygen isotope ratios in the annual tree rings provide useful information about the water-use efficiency and physiological processes of the trees and may provide highresolution information on past and ongoing trends in hydrological changes because the variations in source water δ 18 O are often more important for explaining cellulose-δ 18 O compared to the needle-water enrichments (Treydte et al. 2014).
With the aim of characterizing the impacts of both the climate and the ongoing glacio-hydrological changes on trees colonizing and growing in the neodeglaciated forefield of the Forni Glacier, one of the largest Italian valley glaciers, we (1) analyze the climate impacts on tree-ring growth at two sites where we assume the different importance of glacier meltwater: one site, the GL site, was chosen along the glacier stream, where trees grow on the stream banks very close to the running waters that come directly from the Forni Glacier; the other site, the SL site, was chosen as a control site toward the left valley slope, where the trees are mainly fed by snow meltwaters and precipitation and by secondary, small, slope tributaries. (2) We assess the climate influences on tree-ring stable carbon and stable oxygen isotopes. (3) We estimate the soilwater δ 18 O signature at both the GL and SL sites by means of an isotope model. (4) We evaluate the input of the glacier meltwater δ 18 O signatures on the treering cellulose δ 18 O over time. Our hypotheses are that the tree physiology and growth at the GL site may be influenced by the vicinity of the dynamic environment of the glacier stream, with its daily and seasonal changes in water discharge, but that the strongest signals are found in the δ 18 O signature, with trees at the GL site being more depleted in δ 18 O than those at the SL site.

Study sites
Two forest sites of European larch (Larix decidua Mill.) were selected in the summer of 2011 in the Forni Glacier forefield (central Italian Alps, lat. 46°25ʹ6ʺN; long. 10°34ʹ00ʺE) at an altitude of approximately 2,195 m a.s.l.: the SL site, toward the valley slope, and the GL site, along the glacier stream left bank with trees growing very close to the stream waters ( Figure 1). At each site, five larch trees of the same age and size classes were selected, and two cores per tree were extracted from the stems at a 1.3 m height, which totalled twenty samples.
The Forni Valley is located in the northwestern end of the Valtellina, and its water catchment area is highly glacierized and surrounded by several high-altitude mountains with peaks taller than 3,000 m. The valley at its end hosts the tongue of the the Forni Glacier (lat. 46°24ʹ1ʺN; long. 10°35ʹ10ʺE) that, in 2010, extended approximately 20 km 2 from an altitude of 3,678 m a.s.l. to 2,510 m a.s.l. (WGMS 2012) and is north-northwest oriented. South of the study site on the northeastoriented valley flanks there are another two small glaciers, the San Giacomo and the northwest-San Giacomo (lat. 46°24ʹ8ʺN; long 10°33ʹ20ʺE), which together total approximately 0.16 km 2 and reach, in their lower portions, 3,020 m a.s.l. and 2,920 m a.s.l., respectively. The Forni Glacier tongue retreated from the study area in the valley bottom in 1914-1926 and left several morainic deposits that formed arches in the glacier forefield and on the valley slopes . The soils in the valley bottom are poorly developed; they present only A and C horizons and are 15-30 cm deep. Gravel abundance is high (40-50%) and the pH is acidic (5-5.9). The lithology is characterized by mica-schist and paragneiss of the Pejo unit (Argenton et al. 1980). Currently, the glacier forefield, where the GL and SL sites are located, is mainly colonized by Norway spruce (Picea abies Karst.) and European larch (Larix decidua Mill.). Swiss stone pine (Pinus cembra L.) is nearly absent from the glacier forefield, whereas outside the frontal Holocene moraines and along the lower Forni Valley it forms a pure forest up to the treeline, reaching approximately 2,300 m a.s.l.

Construction of ring-width and stable isotope chronologies
All wood cores were prepared cutting transversal surfaces with a blade and then the ring widths were counted and measured with a 1/100 mm accuracy using a measuring table and the TSAP software (both RINNTECH, Heidelberg, Germany). The twenty individual growth series were cross-dated visually against each other and statistically against the respective mean chronology (using the software COFECHA; Grissino-Mayer 2001), and the age-dependent growth trends were removed from the raw series by fitting a cubic spline with cutoff at two-thirds the series length. Growth indices were obtained as ratios of measured versus modeled values, and a robust bi-weight mean was then calculated (Cook, Shiyatov, and Mazepa 1990), finally obtaining the GL and SL standard chronologies.
To construct the chronologies of stable carbon and oxygen isotope ratios, dated tree rings for the period 1980-2010 were split from the cores using a razor blade. For each site and year, the tree rings were pooled and milled together to a powder less than 0.05 mm, from which the α-cellulose was extracted (method as in Loader et al. 1997). A mass of 0.6 ± 0.1 mg of cellulose per year was weighed in silver capsules to determine δ 18 O by using an isotope ratio mass spectrometer, treating the samples with high-temperature pyrolysis. For δ 13 C, a mass of 0.7 ± 0.1 mg per year was weighed in tin capsules, and the samples were combusted under an excess of oxygen. Stable isotope measurements were done at the IRMS Lab of the Swiss Federal Institute for Forest, Snow and Landscape Research WSL, and crosstests on some cellulose samples of the SL site were performed at the IRMS Lab of the University of Campania. At the WSL, the samples were weighed in silver capsules (aliquots of 0.8-1.0 mg) and pyrolyzed at 1,450°C (PYRO-cube, Elementar, Hanau, Germany) to determine the δ 18 O values of the obtained CO by a Delta Plus XP isotope ratio mass spectrometer (ThermoFinnigan MAT, Bremen, Germany) via a ConFlo III interface (ThermoFinnigan MAT). Isotope ratios were referenced to the international standard VSMOW with a reproducibility of 0.2 per mill (Weigt et al. 2015). At the University of Campania, the carbon and oxygen stable isotope compositions were measured by continuous-flow isotope ratio mass spectrometry (Delta V Advantage, Thermo Scientific, Bremen, Germany) using 0.06 mg and 0.3 mg for 13 C and 18 O determinations, respectively. Isotopic compositions are expressed in delta notation (‰) relative to accepted reference standards: Vienna Pee Dee Belemnite for carbon and Vienna Standard Mean Ocean Water for oxygen. The standard deviation for the repeated analysis of an internal standard (commercial cellulose) was lower than 0.1 per mill for C and lower than 0.5 per mill for O. Both the δ 13 C and δ 18 O chronologies at the SL site have already been presented by Leonelli et al. (2017), as were the related climatic analyses. A correction to the δ 13 C raw series at both sites was applied by means of linear interpolation to compensate for the decreasing δ 13 C in the organic matter related to the fossil-fuel combustion and increasing atmospheric concentration (Francey et al. 1999;McCarroll and Loader 2004).
δ 18 O signature of precipitation, glacier meltwater, and soil Water samples from the glacier stream directly fed by the Forni Glacier site were collected monthly from July to September in 2013, 2014, and again in 2016 (missing 2015 for operational reasons), and the δ 18 O of these water samples was measured. Concurrent with the water sampling, soil-water samples were collected at the rooting depth of approximately 20 cm (soil gravelly C horizon), from July to September 2016, by means of two suction lysimeters equipped with ceramic cups: one of the lysimeters was located at the SL site, and the other was at the GL site ( Figure 1). The suction cup lysimeters (5 cm diameter, 30 cm length) consist of a porous ceramic cup (two-bar air entry pressure) epoxy resin-bonded to a PVC body. Nylon compression fittings threaded into the top cap were used to attach lengths of polyethylene tubing to the sample. A manual pressure-vacuum pump was used to create a depression inside the cup and to recover the sample, once an equilibrium between cup and soil pressures was reached. The two-bar air entry pressure is the pressure required to force air through a saturated ceramic pore, this means that it is suitable to sample in the range of 0-2 bar of soil suction, which is from saturation much beyond the soil field capacity (0.33 bar).
From the Global Network of Isotopes in Precipitation (GNIP, IAEA/WMO 2018) we obtained monthly data of precipitation δ 18 O at Pontresina (50 km west from the study sites) that was completed back in time with the series from the Guttannen, Grimsel (Switzerland), and Villacher Alpe (Austria) stations. All sites in the database were within a range of 200 km from our study site and were at an altitude of more than 1,000 m a.s.l. (Leonelli et al. 2017 (Bowen 2018). We calculated the average differences between the predicted δ 18 O monthly values of the Forni site and the predicted δ 18 O monthly values of the Pontresina station: the mean difference (i.e., the ALG) is −0.68 ± 0.10 per mill; the bar above the term ALG, and above some of the terms for the following equations, means that it is an average value.
Because winter and summer precipitation δ 18 O (W&S_δ 18 O; i.e., the mean δ 18 O of the period from October to March and from May to September) is the main driver of δ 18 O variability at the SL site (Leonelli et al. 2017 The soil water δ 18 O at the SL site for each year i was then linked to W&S_δ 18 O Forni as follows: where f is an unknown fractionation factor. As previously described, soils in the area are very thin and the precipitation waters are expected to enter in the profile and replenish the soil waters almost without fractionations, with the exception of soil-surface evaporative processes: we here assume that the same evaporative processes and the following soil fractionation of minor entities occur at two sites in the same way, changing the δ 18 O of precipitation of the same f factor. The glacier meltwater input δ 18 O (GMWI_δ 18 O) to the cellulose δ 18 O at the GL site was finally estimated for each year i throughout 1980-2010 by a procedure based on the following steps: (1) We first calculated the overall mean tree-ring cellulose δ 18 O enrichment at the SL site () during 1980-2010, considering only the input (precipitation) and output (cellulose) δ 18 O values for each year i and disregarding all the fractionation passages from the soil surface to deeper soil horizons, to the leaves and stem: (2) Because both SL and GL sites are under the same climatic and pedologic conditions but not under the same geomorphologic settings (given the vicinity of the GL site to the glacier stream), we expect that the same fractionation mechanisms in the trees, because of climate and related physiological responses occurring at the SL site, also apply to the GL site. At this site, however, the soil-water δ 18 O may also be influenced by the abundance of glacier meltwaters running very close to the trees. The expected soil water δ 18 O signature at the GL site (equation 4) was therefore estimated by subtracting the mean cellulose δ 18 O enrichment at the SL site (calculated from winter and summer precipitation δ 18 O, before surface fractionation) from the GL cellulose δ 18 O and by adding the same unknown f factor of equation 2 (thus assuming that the same physiological conditions occur at both the SL and GL trees, as well as the same surface fractionations):

Climate data
Meteorological data were collected at the Forni station, located at the Small Dam (in Figure 1), and at the Bormio and Santa Caterina Valfurva stations, within a distance of 16.5 km from the study sites. These stations are maintained by the Regional Environmental Protection Agency (ARPA) of Lombardy, which manages the Lombardy stations of the network of the former Italian Hydrological Service. Their monthly temperature and precipitation records were then subjected to the homogenization procedure described in Brunetti et al. (2006), using also other reference series from the database of Italian temperature and precipitation records set up by the Institute of Atmospheric Sciences and Climate of the National Research Council (ISAC/CNR) and by the University of Milan (Brunetti et al. 2006;Simolo et al. 2010;Crespi et al. 2018). Finally, to set up complete temperature and precipitation records for the study site, the Bormio and Santa Caterina data were used to fill the gaps in the Forni observations. Gap filling was performed by means of the anomaly method, as described in Crespi et al. (2018).

Climate sensitivity
The climate signals in the ring-width and stable-isotope chronologies of the GL and SL sites were detected by correlation analyses using the monthly and seasonal data of the temperature, precipitation, and precipitation δ 18 O from June of the year prior to growth to the current September. The climatic analysis on the SL tree-ring records is the only part of the data in the current research that has already been presented in Leonelli et al. (2017), although with different, dendroclimatic purposes. The analysis was performed during the period 1980-2010, and simple linear regressions between the tree-ring δ 18 O and seasonal temperature of winter (ONDJ), summer (AMJJ), and their aggregations were performed. Finally, the same climatic analyses were applied to the GMWI_δ 18 O series.

Results
The two standard chronologies of European larch constructed at the SL and GL sites show a good agreement in growth patterns during the analyzed period 1980-2010 ( Figure 2). In particular, the similarity index Glk (Gleichläufigkeit; Schweingruber 1988) indicates a statistically significant agreement of the interval trends between the two standard chronologies (with 77 percent of upward or downward concordant trends) and a correlation coefficient greater than 0.5 (Table 1). Throughout the considered period, both series show a maximum growth peak in 1994 and a very pronounced decrease in tree-ring growth variability since 1996. In 1986 and 1987, tree growth at the GL site shows a marked decrease in contrast to the control site SL.
The chronologies of stable carbon isotopes corrected for the decline in atmospheric δ 13 C still show a slightly negative trend during 1980-2010 (Figure 3a). The two series show the highest coherence in terms of Glk and correlation coefficients (Table 1), both of which show highly significant values (p < .001). Both series show a maximum enrichment in δ 13 C in 1992, and the difference between the GL and the SL chronologies shows no trend throughout the considered period ( Figure 3c). The two chronologies of the stable oxygen isotopes show low similarities in interval trends (Figure 3b), but the considered indices still show statistically significant values ( Table 1). The GL site chronology shows a marked negative trend of cellulose δ 18 O depletion since 2001, whereas in the previous period it shows values that are rather similar to those of the SL site and higher values in both 1980 and 1981. The difference between the GL and SL sites shows a marked and a statistically highly significant negative trend over time, with the lowest δ 18 O concentrations measured in 2001 and 2009.
The climatic analysis performed on the standardized ring-width chronologies and the δ 13 C chronologies revealed common responses to the selected climatic variables (Figure 4a, b). Tree-ring growth at both sites is enhanced by the summer months' temperatures, particularly by the July temperature, and by the late summer precipitation (August temperature has an influence only at the SL site). Early summer precipitation up to July, instead, exerts a negative influence on tree-ring growth. June temperature is negatively correlated with δ 13 C at both sites. For δ 18 O, the SL and GL chronologies show almost opposite responses both to temperature and precipitation (Figure 4c, d). The δ 18 O at the GL site is negatively correlated with the monthly and seasonal temperature variables from October of the year prior to growth to January and from April to July, whereas the δ 18 O at the SL site is positively correlated with the August temperature and negatively correlated with March precipitation. Additionally, the response to the precipitation δ 18 O is different at the two sites, with only the SL chronology being positively correlated with previous winter and summer precipitation δ 18 O ( Figure 4e); indeed, the treering cellulose δ 18 O at this site is strongly driven by the δ 18 O from November of the previous year to March (i.e., of winter snowfall), and by the August precipitation δ 18 O (for further details on the signals recorded in the SL chronology, refer to Leonelli et al. [2017]).
The regression of the tree-ring δ 18 O as a function of the winter (ONDJ) and summer (AMJJ) temperatures and their aggregation (Figure 5a-c, respectively) shows the significant distribution of GL points along a line of negative slope, whereas for SL, no significant relationships with these variables are found.
The fractionation processes involved in the treering stable isotope signature can be detected by the observation of the stable isotope distribution of water Table 1. Statistics calculated between standard chronologies (RWI = ring-width index), δ 13 C, and δ 18 O chronologies at the SL and SL sites, during the period 1980-2010 (n = 31): Glk = Gleichläufigkeit (Glk), correlation coefficient (CC), and t value. Significance level: * = p < .05; ** = p < .01; *** = p < .001. GL (Table 2), the δ 18 O general mean from July to August is −9.26 ± 1.06 per mill in precipitation (estimated values from the Pontresina series and equation 1; arithmetic unweighted means of isotope ratios) and −13.75 ± 0.69 per mill in the glacier stream (mean of measured values), with a mean difference of −4.49 ± 1.10 per mill ( Figure 6 and Table 2). The two lysimeters show almost the same values of δ 18 O at both sites. These values are slightly higher than the precipitation signature; however, they show the same trend of the precipitation δ 18 O. shows several values negatively exceeding the soil-water δ 18 O at the SL site and its variability given by ± 1 standard deviation (Figure 7a). The glacier meltwater input (GMWI_δ 18 O), given by the differences in soil-water δ 18 O at the two sites (equation 5) during the considered period is especially evident since 2001 (Figure 7a), when the strongly negative input to the δ 18 O in the treering cellulose at the GL site reached values exceeding −4 per mill. This input is consistent with the glacier meltwater signature, which has more depleted values of approximately −4.49 per mill with respect to the precipitation δ 18 O of summer months ( Table 2).
The regression of GMWI_δ 18 O as a function of winter and summer temperature (W&S temperature; i.e., the mean temperature of the period from October to March and from May to September; Figure 7b), shows that the relation with this temperature variable is stronger than the relation found for the GL cellulose δ 18 O. The W&S temperature explains up to 37 percent of GMWI_δ 18 O variance and only up to 33 percent of GL_δ 18 O cellulose variance (Figure 5c).

Discussion
Our results underline the impacts of water regime changes on the tree-ring cellulose δ 18 O signature of a high-altitude forest colonizing the Forni Glacier forefield. This glacier is still retreating after the short-term re-advance that occurred in the mid-1980s. Its total retreat amounts to more than 700 m during the period 1980-2010 (Figure 8), witnessing the dramatic transformation occurring in these environments that is primarily driven by climate change and warmer temperatures.
The tree-ring growth at the SL and GL sites present similar patterns as expected, because they have trees of similar age and size and the two sites are rather close to each other. Trees at the two sites grow under the same climatic, pedologic, and physiological conditions, with the only exception being the glacier stream influence at the GL site. The two δ 13 C chronologies show a common negative trend in δ 13 C, probably because of the increase of atmospheric CO 2 affecting tree Figure 5. Linear regression between the cellulose δ 18 O at the GL and SL sites (black dots and gray points, respectively) and (a) the mean temperature of the months from October prior to growth (−1) to the current January, (b) the months from April to July, and (c) the aggregation of these two variables. The linear regression equation is reported for all regressions, as well as the corresponding equations, standard errors of the regression parameters, and determination coefficients. Statistically significant regressions: *p < .01; ***p < .001. physiology and stomatal conductance (McCarroll et al. 2009). These chronologies are highly correlated at the two sites, thus highlighting also the similar ecophysiological conditions of the trees growing in the glacier forefield. These similar patterns in tree-ring growth and in δ 13 C are also reflected in the similar responses to climate. Summer temperature, particularly in July, enhances tree-ring growth, whereas precipitation in this month negatively influences tree-ring growth. These opposite response patterns to temperature and precipitation are typical, especially when these climatic factors are negatively correlated with each other, even if larch at high altitudes usually mainly responds to June temperatures (e.g., Coppola et al. 2012). A positive influence of precipitation on treering growth is instead recorded in August (SL) and September (GL). The uncommon negative correlation between δ 13 C and June temperature could be because of the presence of snow on the ground in this month (as also suggested by the winter-precipitation signal of δ 18 O recorded at the SL site); therefore, snowmelt induced by higher temperatures could be a water source for the trees. In the context of higher water supply, a decrease in δ 13 C is then expected. Both the SL and GL sites present this atypical response to temperature, underlining similar environmental conditions. The tree-ring growth of larch in the glacier forefield at both sites seems to be delayed by one month with respect to what is usually found in the literature, thus underlining a possible negative effect exerted by snow meltwaters replenishing the glacier forefield water table, or by cold air drainage off the glacier, pooling in the valley.
In our model for estimating the δ 18 O of glacier meltwater inputs (equation 5), precipitation, snowmelt, and the surface fractionation factor f (equations 2 and 4) are not measured and are assumed as factors of the same magnitude at both sites, because trees grow under the same climatic, altitudinal, and eco-physiological conditions, with the only exception being the glacier influence at the GL site. About the climatic influence on δ 18 O, the nearly opposite responses to temperature and precipitation at the SL and GL sites, respectively, clearly show the influence of different water sources at the two sites. In fact, only at the GL site do we find a negative correlation between previous winter and summer temperatures on the cellulose δ 18 O. This relationship is explainable by the double effect that higher (or lower) temperatures may exert on glacier meltwaters: higher temperatures induce higher ablation rates in the glacier, which indicates, first, higher discharges of the glacier stream and, second, higher ratios of δ 18 O-depleted waters in the glacier stream. Considering the strength of the interannual climatic variability on the cellulose δ 18 O at the GL site (likely linked to glacier melt processes, differently from the SL site), a possible influence of riverbank erosion that could occur along the stream seems a less important factor in determining the recent changes occurring in the soil-water signature at the GL site.
The climatic signals of previous winter and summer precipitation δ 18 O were found only at the SL site, thus underlining the different only meteoric origins of the soil water at this site. Indeed, at the SL site the climatic signal of precipitation δ 18 O is well recorded in the cellulose δ 18 O (Leonelli et al. 2017), while at the GL site, the cellulose δ 18 O is also driven by the glacier inputs. Considering the glacier inputs over time, it is especially since 1999 that major changes were recorded at the GL site, whereas the two sites previously showed almost no differences in soil waters (with the exception Figure 7. (a) Soil-water δ 18 O at the SL site is depicted with a grey line together with ± 1 standard deviation lines (SDEV), and the expected soil water δ 18 O at the GL site during the 1980-2010 period is depicted with a dashed line: in both cases (equations 2 and 4, respectively) the soil surface fractionation factor was set to zero, f = 0. In the bottom, the blue histogram bars depict the annual GMWI δ 18 O calculated with equation 5 and the green line depicts the values of ONDJ-AMJJ mean temperatures expressed as z-scores (−1 refers to months prior to growth). (b) Regression of the GMWI δ 18 O as function of ONDJ-AMJJ mean temperatures is expressed in°C; the regression line also reports the equations as well as the standard errors of the regression parameters and the determination coefficient.
of 1980 and 1981). Interestingly, W&S temperature correlates better with the calculated glacier meltwater inputs (GMWI_δ 18 O) than with the cellulose δ 18 O of the GL site. Linear regressions underline that this climatic variable explains as much as 37 percent of GMWI_δ 18 O variance (Figure 7b) and as much as 33 percent of cellulose δ 18 O variance (Figure 5c). At the SL site, we find a purer climatic signal of precipitation δ 18 O than at the GL site, whose cellulose δ 18 O is also influenced by the glacier meltwater inputs.
As found in east England with oaks (Robertson et al. 2001) or in the Swiss Alps (Treydte et al. 2009), treering cellulose δ 18 O may reflect the precipitation δ 18 O. Additionally, Danis et al. (2006) found that European larch cellulose reflects the signature of previous winter precipitation recharging soil and ground waters in the French Alps. Indeed, in high mountains the isotopic signature of snow may be predominant, as found in an undisturbed subalpine forest of the Rocky Mountain, where Hu et al. (2010) underlined that the xylem water by mid-August may reflect the isotopic signature of snowmelt from 57 percent to 68 percent. Our explanation is that during summer at the Forni Valley study site, the slope water tables are replenished by snowmelt and precipitation, whereas along the glacier stream the isotopic signature of the glacier meltwaters may be predominant on the snowmelt and precipitation signals. Our finding indicates that trees along the glacier stream are fed also by the glacier stream waters and may record mid-to long-term changes in the water regimes of the glacier stream. Tracking source water δ 18 O changes by analyzing the tree-ring cellulose δ 18 O and by applying the proposed isotope model may help in reconstructing past and ongoing changes in the environment (e.g., Ellsworth et al. 2013;Bose et al. 2016) especially in remote sites where few on-site stable isotope data are available and long series are lacking.

Conclusions
We have demonstrated that in the past decade of analysis trees in the glacier forefield of a highly glacierized Alpine catchment growing along the glacier stream are mainly fed by meltwaters of glacial origins, whereas trees growing in the glacier forefield toward the valley slope are mainly fed by waters of meteoric origins, as reflected in their respective cellulose δ 18 O. The negative correlation between (previous) winter and summer temperatures with the δ 18 O of glacier meltwater inputs in soil waters at the GL site clearly shows the glacier influence on treering isotopic signature. The tree-ring cellulose δ 18 O isotopic signature of the trees growing along the glacier stream underlines a recent period of stronger δ 18 O depletion that may be linked to either a higher glacier discharge or a higher concentration of glacier meltwaters in the stream. Both scenarios are compatible with the ongoing trends of higher atmospheric temperatures inducing higher ablation rates of the Forni Glacier. Further field measurements could investigate the inner mechanisms explaining the more δ 18 O-depleted soil waters along the stream banks; however, at this stage we can state that these changes were recorded since 1999.
Because of the ongoing temperature warming, glaciers are retreating almost everywhere on the Alps (WGMS 2017); consequently, glacier forefields are rapidly expanding and are a place of forest succession if they are at an altitude below the treeline. These forests, even if they are marginal in terms of land cover and productivity, can be seen as sentinels of the ongoing climatic and hydrologic changes that involve the European Alps and highmountain areas in general. Monitoring the hydrological changes in these sensitive forests in the mid-to long term will allow the reconstruction of the ongoing glaciohydrologic dynamics and the assessment of climatechange impacts on mountain forests for what concerns changes in glacier meltwater regimes. Our tree-ring sampling approach, together with the proposed isotope model for estimating δ 18 O meltwater inputs, opens new opportunities to reconstruct mid-to long-term changes in glacier meltwater streams in remote sites.
Roberto Comolli for suggestions on soil-water sampling; Dr. Giuliano Dreossi for measuring stable isotopes in water samples from 2016 at the Università Ca' Foscari, Venezia; Loic Schneider and Stella Fusco for the lab analyses at WSL; and Aurora Merlini for her help in the water-sample collection during the summer of 2014.