ORIGINAL RESEARCH article
Sec. Functional Plant Ecology
Stable Water Use Efficiency under Climate Change of Three Sympatric Conifer Species at the Alpine Treeline
- 1Department of Alpine Timberline Ecophysiology, Federal Research and Training Centre for Forests, Natural Hazards and Landscape, Innsbruck, Austria
- 2Institute of Botany, Leopold-Franzens-Universität Innsbruck, Innsbruck, Austria
- 3Ecophysiology of Plants, Department of Ecology and Ecosystem Management, Technische Universität München, Freising, Germany
The ability of treeline associated conifers in the Central Alps to cope with recent climate warming and increasing CO2 concentration is still poorly understood. We determined tree ring stable carbon and oxygen isotope ratios of Pinus cembra, Picea abies, and Larix decidua trees from 1975 to 2010. Stable isotope ratios were compared with leaf level gas exchange measurements carried out in situ between 1979 and 2007. Results indicate that tree ring derived intrinsic water-use efficiency (iWUE) of P. cembra, P. abies and L. decidua remained constant during the last 36 years despite climate warming and rising atmospheric CO2. Temporal patterns in Δ13C and Δ18O mirrored leaf level gas exchange assessments, suggesting parallel increases of CO2-fixation and stomatal conductance of treeline conifer species. As at the study site soil water availability was not a limiting factor iWUE remained largely stable throughout the study period. The stability in iWUE was accompanied by an increase in basal area increment (BAI) suggesting that treeline trees benefit from both recent climate warming and CO2 fertilization. Finally, our results suggest that iWUE may not change species composition at treeline in the Austrian Alps due to similar ecophysiological responses to climatic changes of the three sympatric study species.
High-altitude forest ecosystems at the timberline-treeline transition have raised concern as they may undergo significant alterations due to climate warming and changes in ground-level air chemistry (Holtmeier and Broll, 2007; Wieser et al., 2009). Dendroclimatological studies conducted within the treeline ecotone of the Central European Alps have shown radial stem growth to be limited by low summer temperature (Carrer and Urbinati, 2004; Oberhuber, 2004; Büntgen et al., 2005; Frank and Esper, 2005; Oberhuber et al., 2008). During recent decades, several authors report treeline-associated conifers to reflect increased radial growth, putatively related to climate warming (Graumlich et al., 1989; Peterson et al., 1990; Jacoby and D’Arrigo, 1997; Rolland and Florence-Schueller, 1998; Bunn et al., 2005). Moreover, increasing atmospheric CO2 concentration may act in concert with climate warming to increase carbon accumulation within the treeline ecotone (cf. Graumlich, 1991; Saurer et al., 1997; Duquesnay et al., 1998; Sidorova et al., 2009). Notwithstanding, reduced radial growth has been attributed to late-summer drought under increasing treeline temperature in the European Alps (Büntgen et al., 2006; Carrer and Urbinati, 2006; Oberhuber et al., 2008; Wieser et al., 2009).
Stable isotope ratios of carbon and oxygen, i.e., 13C/12C and 18O/16O, respectively, may serve as dendrochronological proxies that facilitate mechanistic understanding of climate-related influences on physiological processes such as leaf gas exchange and stem wood formation (Loader et al., 2007; Weigt et al., 2015; and further references therein). In plant organic matter, δ13C expressing the 13C/12C ratio in relation to an international standard (Pee Dee Belimnite) depends on variables such as leaf conductance for water vapor (gw) that modify the net CO2 uptake rate (A; Farquhar et al., 1989). In addition, δ13C of plant organic matter (δ13Cp) is a function of atmospheric δ13C, which is accounted for by calculating discrimination of photosynthesis against 13C (Δ13C; Farquhar and Richards, 1984) in relation to intrinsic water-use efficiency (iWUE), i.e., the ratio of A (rate of net CO2 fixation) versus gw (leaf conductance for water vapor). For a review see Brugnoli and Farquhar (2000).
When analyzing δ13Cp alone, the impacts of A (demand of CO2) and gw (supply of CO2) on iWUE are difficult to separate (Saurer et al., 2008). The oxygen isotope ratio (δ18O), however, may allow a distinction between biochemical and stomatal limitations of photosynthesis as it is not affected by the photosynthetic CO2 carboxylation but linked to gw (Barbour et al., 2000; Grams et al., 2007). It is, therefore, an ideal covariable to estimate to what degree photosynthesis and stomatal conductance modify δ13Cp (Scheidegger et al., 2000; Werner et al., 2012). Originally, this dual isotope approach was introduced for photosynthetic tissue and only recently tested conceptually for the interpretation of tree-ring data (Roden and Farquhar, 2012), although several critical points should be taken into account during interpretion. Among the most important issues that need to be considered are the facts that the δ18O of source and atmospheric water can vary spatially and temporarily and that post-photosynthetic and post-evaporative oxygen atom exchange processes could affect the initial leaf-level isotope signal (see below).
At the leaf level, δ18O of photoassimilates derive primarily from leaf water, typically being enriched in 18O compared to the source water (i.e., xylem water) through evaporative enrichment at the site of transpiration. This enrichment is counteracted by the so-called Péclet effect and transpiratory leaf cooling (for a review see Barbour, 2007), which may result in the above mentioned negative correlation between δ18Op and gw (e.g., Barbour et al., 2000; Grams et al., 2007). However, to an extent that may depend on species and site conditions, the signal is dampend by oxygen exchange with source water during biomass formation at the stem level (Gessler et al., 2013). This causes, at least partially, a decoupling between oxygen isotopic signatures of photoassimilates and the tree ring organic matter. However, in a recent report Weigt et al. (2015) confirmed that information may be exploited to relate ecophysiological responses of trees to environmental changes for both total sapwood organic matter and extracted cellulose. In any case, it appears advisable to confirm interpretation from the dual isotope approach by gas exchange assessments whenever possible.
Previous tree-ring carbon isotope studies carried out in tropical, arid, Mediterranean, temperate and boreal forest ecosystems have identified an increase in iWUE over the past 40 years in response to increasing atmospheric CO2 concentration (Penuelas et al., 2011; Saurer et al., 2014) and just recently supported by tree-ring δ13C dynamic global vegetation models comparisons (DGVMs; Frank et al., 2015). Tree growth on the contrary, remained stable or even declined, suggesting local site conditions to override a potential CO2-induced increase in growth (Penuelas et al., 2011; Silva and Anand, 2013; Levesque et al., 2014). Increasing iWUE accompanied by a reduced productivity has been attributed to the combined effect of elevated CO2 and climate change-induced soil drying (Penuelas et al., 2011; Saurer et al., 2014).
The impact of the steadily increasing CO2 level and concurrent climate change, however, still awaits clarification for treeline-associated conifers in the Central Austrian Alps where low temperatures limit tree growth (Oberhuber, 2007; Wieser et al., 2009). At treeline in the Central Austrian Alps ample precipitation during the growing season prevails every third to fourth day on average (Wieser, 2012), so that soil water limitation stays absent, allowing trees to meet their water demand (Tranquillini, 1979; Mayr, 2007; Matyssek et al., 2009). Hence, whole-tree conductance stays high and mainly depends on the evaporative demand in terms of irradiance and vapor pressure deficit (Wieser, 2012). Therefore, we hypothesize that treeline trees passively respond to the increasing atmospheric CO2 level (Ca), so that their leaf-intercellular CO2 concentration (Ci) rises in parallel, while iWUE remains unchanged. The hypothesis was evaluated by stable carbon and oxygen isotope sampling and radial growth analysis over the past 36 years (1975–2010) in stems of mature Pinus cembra, Picea abies and Larix decidua trees growing at the treeline of Mt. Patscherkofel in the Central Tyrolean Alps. Observed long-term trends in δ13Cp, δ18Op in tree rings and iWUE were at least for some years compared with in situ leaf-level gas exchange data, assessed at the same study site between 1979 and 2007 in adult P. cembra and L. decidua trees. Results are discussed in view of tree response to climate warming at the treeline ecotone.
Materials and Methods
Study Site, Climatic Data, and Tree Species
The study was conducted in a scattered stand at the lower edge of the treeline ecotone at 1950 m a.s.l. on Mt. Patscherkofel (47°12′37″ N, 11°27′07″ E), south of Innsbruck, Austria. The site is characterized by a cool subalpine climate, the possibility of frost during the entire year and a continuous snow cover from October through April. We used monthly mean temperatures and monthly total precipitation from 1975 to 2010 from a weather station nearby (Klimahaus Research Station and Alpengarten; 1950 m a.s.l.) for our analysis. Mean annual precipitation averaged 878 mm, with 58% falling during the growing season (May through September). Mean annual air temperature averaged 2.4°C, with summer maxima of up to 27°C and winter minima of -28°C.
The geology of Mt. Patscherkofel is dominated by gneisses and schist. The soil at the study site is a haplic podzol, being a typical soil type of the treeline ecotone in the Central Tyrolean Alps (Neuwinger, 1972). The water holding capacity of the soil (at 5–65 cm depth) at saturation (-0.001 MPa) averages 0.60 m3 m-3. Due to frequent precipitation during the growing season soil water potential rarely drops below -0.01 MPa, approximating soil water contents above 0.35 m3 m-3 (Wieser, 2012; including the dry summer of 2003, unpublished data).
The stand is composed of the dominant tree species Pinus cembra, accounting for 84% of the tree population, and accompanied by Larix decidua (9%) and Picea abies L. Karst (7%) at some locations. Trees grew either as isolated trees or in groups of four to five. The distance between single trees or tree groups was 20–30 m. From 20 P. cembra trees cored at the lower edge of the treeline ecotone Oberhuber et al. (2008) derived an expressed signal population (EPS) value of 0.94, reflecting a strong climate signal in the site chronology. From these trees we selected five trees which had the strongest correlation to the site specific mean tree-ring chronology, no missing rings, and regular ring boundaries. In addition we cored five dominant P. abies and L. decidua trees each to account for potential inter-specific differences of the three associated treeline species. In 2010 the trees were 69 ± 9 years old, with stem heights averaging 12 ± 1.3 m. The stem diameter at breast height (DBH) averaged 22 ± 3.2 cm.
Tree Ring and Basal Stem Area Increment
In fall 2010 we obtained two increment cores per trees at DBH using a 5-mm-diameter increment bore. For contrast enhancement of tree ring boundaries the cores were dried in the laboratory, non-permanently mounted on a holder, and the surface was prepared with a razor blade (Pilcher, 1990). Ring widths were measured to the nearest 1 μm using a reflecting microscope (Olympus SZ61) and the software package TSAP WIN Scientific. Ring widths of both cores from each sample tree were averaged and individual tree ring chronologies were then checked for dating accuracy using the COFECHA software (Holmes, 1994; Grissino-Mayer, 2001). As ring width may be biased by a negative correlation with the time course during tree maturation, ring width was converted to basal stem area increment (BAI) according to:
where R is stem radius inside tree bark and n is the year of tree ring formation (Fritts, 1976). Bark thickness was subtracted from stem radius. Finally BAI of each year were averaged over the five sample trees of each species.
Stable Isotope Analysis
δ13C and δ18O analyses for the years 1975–2010 were performed on the same cores as used for BAI assessment. Annual rings (early wood plus late wood) were cut exactly at ring boundaries by use of a scalpel and a reflecting microscope (Wild 308700). For each of the five study trees per species the two samples per tree ring were pooled and homogenized with a swing mill (Retsch MM301, Retsch Haan, Germany). In a subsample, we compared isotope signatures in bulk wood with those in cellulose for determining the necessity of cellulose extraction in our study trees. Cellulose extraction was performed using a modified version of the method of Brendel et al. (2000). The methodological comparison corroborated significant correlations in the cases of δ13C and δ18O (Figure 1) as reported earlier from coniferous and other tree species (Jaggi et al., 2002; Sohn et al., 2013) and is in accordance with a recent report (Weigt et al., 2015). On average, δ13C in cellulose was 1.0–1.1 ‰ higher than in bulk wood (Table 1), being smaller than 1.3–1.4 ‰ found in Picea abies by Borella et al. (1998) and Sohn et al. (2013). Mean δ18O in cellulose was 4.3–4.9 ‰ higher than in bulk wood (Table 1) being somewhat lower than 5.9 ‰ in bulk wood of Picea abies (Sohn et al., 2013). Based on these findings and in accordance with a recent methodological study (Weigt et al., 2015), we used bulk wood samples rather than extracted cellulose for isotope analysis.
FIGURE 1. Comparison of (A) δ13C and (B) δ18O in bulk wood and cellulose of annual growth rings of P. cembra (solid circle), P. abies (solid square) and L. decidua (open circle). Dashed lines reflect the one-to-one lines for comparison and the solid lines correspond to the linear regression analyses of the data points: δ13C: y = 1.03x + 1.86, r2= 0.92, P < 0.001, n = 36; δ18O; y = 0.92x + 6.65, r2= 0.73, P < 0.001, n = 18.
TABLE 1. δ13C and δ18O difference between cellulose and bulk wood in annual growth rings of Pinus cembra, Picea abies and Larix decidua.
Regarding δ13C, 2.0 ± 0.02 mg of homogenized samples were weighed into tin capsules each (3.5 × 5 mm, IVA Analysentechnik e.K., Meerbush, Germany) and combusted to CO2 in an elemental analyzer (Eurovector EA3000) connected to an isotope ratio mass spectrometer (Isoprime, Elementar, Hanau, Germany). For δ18O analysis 0.7 ± 0.05 mg were weighed into silver capsules each (3.5 × 5 mm, IVA Analysentechnik e.K., Meerbush, Germany) to obtain CO at 1,430°C in a high-temperature pyrolysis system (HTO, HekaTech, Wegberg, Germany) which was connected via an open-split interface (Conflo III; Finnigan MAT, Bremen, Germany) to an isotope ratio mass spectrometer (Delta Plus; Finnigan MAT, Bremen, Germany). Isotope abundances were expressed using the δ-notation in ‰ relative to the international standards:
where Rsample is the molar fraction of the 13C/12C or 18O/16O ratio of the sample, and Rstandard that of the international IAEA standards V-PDB for carbon and V-SMOW for O. The analytical precision was <0.12 ‰ and <0.28 ‰ regarding δ13C and δ18O, respectively (expressed as standard deviation of the internal laboratory standard at the same sample mass).
Isotope Discrimination and iWUE
Tree ring specific δ13Ctring were corrected for the progressive decline in atmospheric δ13Catm through calculating 13C discrimination (Δ13C):
To this end, δ13Catm with its nearly linear time course during 1980 through 20101 was extrapolated for the years 1975 through 1979. In a simplified model, Farquhar et al. (1982) related Δ13C through plant physiological processes during CO2 fixation in C3 plants with the ratio of intercellular to ambient CO2 partial pressure (Ci/Ca):
where a (=4.4 ‰) refers to the slower diffusivity of 13CO2 relative to 12CO2 in air and b (=27 ‰) is the isotopic fractionation caused by enzymatic C fixation. Ca was obtained from published data2. It should be noted that Δ13C is determined by the ratio of chloroplast to the ambient CO2 mole fraction (Cc/Ca) rather than Ci/Ca, as used in equation 3, making the here calculated value sensitive to mesophyll conductance (gm; Seibt et al., 2008). The latter varies in accordance to changes in environmental conditions such as temperature, irradiance, water and CO2 availability (Flexas et al., 2008). Consequently using Ca may be problematic if gm to CO2 is not constant (Seibt et al., 2008). However, as information on mesophyll conductance of the three conifers under study is not available and published means of gm would not improve results (Cernusak et al., 2013), we chose using the simplified linear model of Farquhar et al. (1982). Hence, iWUE, i.e., the ratio of the net carbon gain (A) versus leaf conductance for water vapor (gw), was calculated as follows:
where 1.6 is the ratio between the diffusivities of water vapor and CO2 in air.
Enrichment in 18O in tree rings over source water (Δ18O), resulting from incorporation of 18O-enriched photoassimilates into stem biomass, was calculated from d18O of tree ring organic matter (δ18Otring) and precipitation (δ18Oprep) according to:
Precipitation is the only source of water on Mt. Patscherkofel and thus was assumed to reflect source water of trees. Furthermore, as there is evidence from a treeline in the Central Swiss Alps that mean values of δ18Oprep and δ18O of soil water do not differ significantly, although soil water δ18O carries a distinct signal from snow melt water far into the growing season (Treydte et al., 2014), we used annual means of δ18Oprep sampled on top of Mt. Patscherkofel (2246 m a.s.l.) (Umweltbundesamt Austria; personal communication) approximately 300 m south of the selected study trees for our Δ18O calculation. We elevationally corrected δ18Oprepby a factor of 0.17 ‰ per 100 m of elevation (Umweltbundesamt Austria; personal communication) which is considerably lower than the global mean of 0.28 ‰ per 100 m of elevation (Poage and Chamberlain, 2001).
Leaf Level Gas Exchange Data
To illustrate long term trends in foliar CO2 and H2O gas exchange we compiled published gas exchange data of mature, field grown P. cembra and L. decidua trees carried out at our study site between 1979 and 2007 (Table 2). Maximum net photosynthetic capacity at ambient CO2 (Amax; sensu Larcher, 2001) was assessed of sun exposed twigs from the upper canopy. Employed in situ were thermoelectrically climate-controlled cuvettes (Walz, Effeltrich, Germany) or a portable exchange system (CIRAS 1, PP Systems, Hitchin, Hertfordshire, UK). For methodological details see publications given in Table 2.
TABLE 2. A comparison of published maximum net photosynthetic capacity (Amax) and leaf conductance for water vapor of sun exposed shoots of mature Pinus cembra and Larix decidua trees at the lower end of the treeline ecotone on Mt. Patscherkofel.
Temperature, precipitation, vapor pressure deficit, Δ13C, Δ18O, BAI, Ci, Ci/Ca and iWUE trends were calculated for the time period 1975–2010 by least-squares linear regression analysis. For a given variable, differences among trends (slopes) between P. cembra, P. abies and L. decidua were assessed by the two-slope comparison test (Zar, 1999). We used repeated measures ANOVA to detect significant differences in the mean values (1975–2010) of Δ13C, Δ18O, BAI, Ci, Ci/Ca and iWUE of P. cembra, P. abies and L. decidua. Following Kunter et al. (2004) we used multiple least-squares linear regression models to assess the influence of atmospheric CO2 concentration (Ca) and mean growing season (May-Sep) air temperature (Tveg) and their interactions (explanatory variables) on tree-ring variables. For assessing the climatic impact on tree ring variables (BAI, Δ13C, and Δ18O) statistical analyses were based on mean monthly air temperature (°C) and total monthly precipitation (mm) throughout the study period (1975–2010). For each species Pearson’s correlation coefficients between BAI, isotope chronologies and both climate variables were calculated from August of the year prior to growth to September of the growth year. All the statistical analysis were conducted by use of the SPSS 16 software package (SPSS. Inc. Chicago, IL, USA), and a probability level of P < 0.05 was considered as statistically significant. As suggested by Sarris et al. (2013) we did not remove any age related trend from our tree-ring chronologies by conventional detrending procedures, thus avoiding the risk of removing any environmental signal or trend captured by our tree-ring series.
Inter-annual Trends in Climate and Tree-Ring Indices
A warming trend is reflected at our treeline site during the growing seasons (0.50°C per decade P < 0.001) of 1975–2010, without concurrent trends in precipitation and vapor pressure deficit (Figure 2). During the whole study period Δ13C and Δ18O chronologies were synchronized between the three studied species. In each species Δ13C increased over time (Figure 3A; Table 3) whereas Δ18O decline (Figure 3B; Table 3). The increase in Δ13C was accompanied by rising of Amax for both P. cembra (1979–2007) and L. decidua (1980–1993) by about 50%. Likewise, gw increase by about 75% in L. decidua (Table 3).
FIGURE 2. Temporal variation in (A) growing season mean air temperature (Tair), (B) total growing season precipitation (P) and (C) growing season mean vapor pressure deficit (VPD) during the period 1975 throughout 2010. Data were fit by linear regression analysis: Tair: y = 0.053x-97.0, r2= 0.30, P < 0.001; P: y = 0.075x-365.3, r2= 0.00, P = 0.95; VPD: y = -0.055x + 13.3, r2= 0.01, P = 0.49.
FIGURE 3. Δ13C (A), Δ18O (B), and BAI (C) chronologies of P. cembra (solid circle, solid line), P. abies (solid square, dashed line) and L. decidua (open circle, dotted line) between 1975 and 2010. See Table 2 for regression information.
The mean tree-ring Δ13C was highest in L. decicua, although the increase was significantly higher in P. cembra and P. abies (Table 4). Temporal changes in Δ18O by contrast, did not differ significantly between the tree species (Table 4). P. cembra showed the highest Δ18O while P. abies presented the lowest Δ18O and L. decidua displayed an intermediate mean (Table 4). On average, growth of P. cembra was significantly higher than growth of P. abies and L. decidua (Figure 3C; Table 4). During the whole study period all three species showed an increase in growth expressed as BAI, being significantly lower in P. cembra than P. abies, and L. decidua (Figure 3C; Table 4).
TABLE 4. Tree-ring carbon isotope characteristics (Δ13C, Δ18O, BAI, Ci, Ci/Ca, and iWUE) in P. cembra, P. abies, and L. decidua during the period 1975–2010.
Paralleling atmospheric CO2 enhancement (Figure 4A), tree ring derived Ci increased from 1975 through 2010 from 180 to 234 μmol mol-1 in P. cembra, from 185 to 227 μmol mol-1 in P. abies and from 213 to 256 μmol mol-1 in L. decidua (Figure 4A; Table 2). Although species specific differences in the temporal change of Ci were not statistically significant different from each other, mean Ci was significantly lower in P. cembra and P. abies as compared to L. decidua (Table 4). Averaged over the study period P. cembra showed the lowest and L. decidua, the highest Ci/Ca, while the Ci/Ca of P. abies was intermediate (Table 4). The increase in Ci/Ca over time (Figure 4B; Table 3) was significantly higher in P. cembra and P. abies than in L. decidua (Table 4). In all the three species iWUE had remained stable during the study period (Figure 4C; Table 3). However, we observed statistically significant between species, with P. cembra showing the highest and Larix decidua showing the lowest iWUE averaged over the study period (Table 4).
FIGURE 4. (A) Ambient armospheric CO2 concentration (Ca, thin solid line) Ci, (B) Ci/Ca, and (C) iWUE chronologies of P. cembra (solid circle, solid line), P. abies (solid square, dashed line) and L. decidua (open circle, dotted line) between 1975 and 2010. See Table 2 for regression information.
Effects Ca and Tveg on Δ13C, Δ18O, BAI, Ci, Ci/Ca and iWUE
Multiple linear regression analysis show that Δ13C and Δ18O of all species significantly increased with increasing ambient CO2 concentration (Ca) while growing season mean air temperature (Tveg) had no effect on Δ13C and on Δ18O (Table 5). Growth (BAI) of P. cembra, and P. abies significantly increased with increasing Ca and Tveg. L. decidua presented a significant increase in BAI at increasing Ca without any response to Tveg (Table 5). For all species we found a significant increase in Ci and Ci/Ca at higher Ca but not at higher Tveg (Table 5). iWUE of P. cembra, P. abies, and L. decidua, however, did not significantly respond to increasing Ca and Tveg (Table 5).
TABLE 5. Summary of multiple linear regression models fitted to explain inter-annual changes (1975–2010) in Δ13C, Δ18O, BAI, Ci, Ci/Ca, and iWUE of P. cembra, P. abies, and L. decidua in response to atmospheric CO2 concentration (Ca) and mean growing season (May-Sep) air temperature (Tveg).
Δ13C, Δ18O, and BAI Response to Climate (Climate-Growth Relationships)
Climate-response relationships of Δ13C, Δ18O, and BAI differed both in time and in signal strength (Figures 5 and 6). In all three species Δ13C was significantly positive correlated with April throughout June temperatures (Figure 5A) and significantly negative correlated with January precipitation (Figure 6A). Previous-year August and October temperature also favored Δ13C in P. cembra and in P. abies, respectively (Figure 5A), whereas previous- and current-year August precipitation did so in P. abies (Figure 6A).
FIGURE 5. Pearson correlation coefficients between (A) Δ13C (B) Δ18O, and (C) BAI, and monthly mean temperature for P. cembra (solid circle, solid line), P. abies (solid square, dashed line) and L. decidua (open circle, dotted line) between 1975 and 2010. Horizontal gray lines indicates P < 0.05.
FIGURE 6. Pearson correlation coefficients between (A) Δ13C (B) Δ18O, and (C) BAI, and monthly total precipitation for P. cembra (solid circle, solid line), P. abies (solid square, dashed line) and L. decidua (open circle, dotted line) between 1975 and 2010. Horizontal gray lines indicates P < 0.05.
The effects of temperature and precipitation on Δ18O were clearly in opposite directions (Figures 5B and 6B). In all the three species tree-ring Δ18O was negative correlated to air temperature from April to June of the current year and significantly positive correlated to previous-year December temperature (Figure 5B) as well as to January precipitation (Figure 6B). From Figure 5B, previous- year August temperature had a negative correlation with L. decidua and current-year January temperature showed a negative correlation with L. decidua. Previous-year November and current-year March precipitation showed a negative correlation with Δ18O in L. decidua, as did June precipitation in P. cembra (Figure 6B).
We also found significant positive correlations between BAI and temperature during April, May, and June in P. cembra, P. abies, and L. decidua (Figure 5C). Current-year August temperature also favored BAI in L. decidua and P. abies, and also previous-year October temperature in P. abies (Figure 5C). The correlations between BAI and precipitation were weak, except for significant positive correlations in February and current-year August in P. cembra and a significant negative correlation in January in P. abies (Figure 6C).
Similar growth and Δ13C, and Δδ18O responses were found over time in P. cembra, P. abies and L. decidua at the lower edge of the treeline ecotone in the central Austrian Alps. From 1975 throughout 2010 the three species increased Δ13C and BAI, while Δ18O showed a declining trend. Apparently, underlying response mechanisms were similar across the three studied species.
Our observed correlations for temperature and precipitation for Δ13C and Δ18O (Figures 5 and 6) are consistent with results reported for oak and pine trees at temperate sites in Switzerland (Saurer et al., 2008). Weather conditions prevailing during April through June predominantly were responsible for variations in tree-ring Δ13C, Δ18O, and BAI of P. cembra, P. abies, and L. decidua. We found positive correlations between April to June temperatures and Δ13C. Δ13C is strongly affected by net CO2 uptake rates, which at treeline are governed by both photon flux density and temperature (Treydte et al., 2001; McCarrol and Pawellek, 2004; Kress et al., 2011) as well as enhanced plant transpiration (Liu et al., 2015). Moreover, at our study site in situ net photosynthetic capacity of sun exposed twigs from the upper canopy of mature P. cembra and L. decidua trees measured under clear summer days also tended to increase between 1979 and 2007 (Table 2), which might be attributed to both the observed increase in atmospheric Ca and Tair. Additionally, the temperature optimum of Amax for P. cembra increased from 12.5°C in 1956 (Pisek et al., 1969, 1973) to 15.0°C in 2002 (Wieser, 2004) and to 17.1°C in 2007 (Wieser et al., 2010), matching the observed increase in mean growing season air temperature of 0.9°C per decade (Figure 2). An increase in net photosynthetic rates under elevated CO2 was also observed in P. mugo and L. decidua after nine years of free-air CO2 enrichment at the Swiss treeline (Dawes et al., 2013; Streit et al., 2014). Three years of ecosystem warming also increased carbon uptake of Pinus cembra at treeline in the Austrian Alps (Wieser et al., 2015).
Our precipitation signals suggest that trees do not suffer from moisture stress. Indeed the observed declining trend Δ18O and the strong negative correlations between Δ18O and growing season precipitation is consistent with the physiological isotopic responses (Barbour and Farquhar, 2000), suggesting that stomatal conductance is increased during the study period. Although a leaf physiological signal in δ18O will be dampened at the level of tree rings due to oxygen exchange with source water during cellulose biosynthesis, impact of gw on 18O in tree rings may be still detectable, even in whole wood analyses (Weigt et al., 2015). Ecosystem warming accompanied by unchanged VPD also increased gw and hence also transpiration in boreal Picea abies (Bergh and Linder, 1999), Pinus sylvestris (Kellomäki and Wang, 1998), Picea mariana (Van Herk et al., 2011), Pinus cembra at treeline (Wieser et al., 2015), and Populus deltoides (Barron-Gafford et al., 2007). Thus, it seems that in cold environments under non-limiting water availability increasing temperatures counteract the diminishing effect of rising CO2 on leaf conductance (Saurer et al., 2014).
The observed positive correlations between BAI and April–June temperatures are also reflected in wood formation. At the study site wood formation of larch, spruce and pine generally starts in May, reaches its maximum in June, and terminates in August (Havranek, 1981; Loris, 1981; Gruber et al., 2009). Beside summer temperatures (Figure 5C) other climatic variable like winter and August precipitation (Figure 6C) are also known to influence radial growth at treeline as shown for P. cembra by Oberhuber (2004), reflecting minor soil water effects on tree growth at treeline (Tranquillini, 1979; Wieser, 2004). Although treeline trees are saturated with carbohydrates (Gruber et al., 2011), growth of trees at treeline is primarily affected by temperature dependent carbon sink activity during tissue formation (Hoch and Körner, 2003, 2012). The observed increase in BAI (Figure 3C) suggests that our treeline trees benefit from climate warming, although effects of CO2 fertilization on growth may not completely ruled out (Table 5). Four years of experimental air warming with open-top-chambers also stimulated radial growth of Picea glauca seedlings at the subarctic treeline in southwest Yukon, Canada (Danby and Hik, 2007). Thus, when growth is stimulated and there is plenty of water gw can increase as indicated by a decline in Δ18O (Figure 3B) along with increasing A, resulting in a constant iWUE (Figure 4C).
Elevated atmospheric CO2 is expected to affect plant carbon-water relationships, as a decline in stomatal conductance is often observed when plants are exposed to elevated CO2 (Battipaglia et al., 2013). If stomatal conductance declines under increasing CO2 in combination with an increased or unchanged carbon assimilation, this will decrease the Ci to Ca ratio and thus decrease Δ13C. Conversely, in all three study species Δ13C increased from 1975 throughout 2010, while tree-ring derived iWUE remained stable (Figure 4C) although ambient CO2 concentration increased by 60 μmol mol-1 (Figure 2B). The stability of iWUE resulted as Ci drifted upward paralleling the rise in Ca (Figure 2B). Likewise, in Picea schrenkiana at treeline in the western Tianshan Mountains in China iWUE remained also unchanged from 1985 to 2010 (Wu et al., 2015). No change in iWUE (i.e., homeostasis) over the last 100 years was also reported for three conifer species in the Selkirk Range (Rocky Mountains, Idaho, ID, USA) by Marshall and Monserud (1996). Other studies by contrast observed a 20% increase in iWUE from the 1960 throughout 2000 in mature trees in tropical, arid, Mediterranean, wet temperate and boreal forests distributed through Europe, Asia, Africa, America, and Oceania (Penuelas et al., 2011; Saurer et al., 2014; Frank et al., 2015). In these latter studies, increasing iWUE was attributed to the combined effect of increasing CO2 and climate change-induced soil drying that reduced stomatal aperture. Soil drought can be ruled out along the treeline ecotone of the Central Alps (Mayr, 2007; Wieser, 2012). Occurrence of soil drought strongly depends on site conditions such as precipitation patterns, water holding capacity of the soil, and evaporative demand. Ample precipitation and moderate evaporative demand in general cause soil water availability to be sufficiently high to meet the trees’ water demand at treeline in the Central Tyrolean Alps (Mayr, 2007; Wieser et al., 2009). As a consequence, treeline trees are rarely forced to restrict transpiration (Tranquillini, 1979; Benecke et al., 1981; Matyssek et al., 2009; Wieser and Leo, 2012; Wieser et al., 2014, 2015). Given the ample soil water availability whole-tree conductance of P. cembra, P. abies, and L. decidua remains high for CO2 uptake because leaf conductance for water vapor depends only on the evaporative demand driven by irradiance and vapor pressure deficit (Wieser, 2012).
Beside climate warming and increasing Ca, nitrogen deposition could also be important for explaining the observed increase in tree growth as increasing nitrogen deposition during the 1980ties (Smidt and Mutsch, 1993) has been suggested as a possible growth stimulator. However, there is evidence that nitrogen contents per needle dry mass are higher in trees at treeline as compared to trees growing at lower elevation sites (Körner, 1989; Birmann and Körner, 2009). Furthermore, since 1988, nitrogen deposition at treeline in the Tyrolean Alps is steadily declining (Amt der Tiroler Landesregierung, 2015), and a nitrogen fertilizer experiment at the alpine treeline in the Swiss Alps showed little or no growth stimulation (Keller, 1970). Thus, it seems that presently nitrogen deposition is insufficient to explain observed growth trends at treeline as reported previously by Tranquillini (1979) and Nicolussi et al. (1995).
Treeline trees respond to the increasing atmospheric CO2 level (Ca) in a way that their leaf-intercellular CO2 concentration (Ci) drifted upward paralleling the rise in atmospheric CO2 while iWUE remained stable over the last 36 years. The stability in iWUE was accompanied by an increase in BAI suggesting that treeline trees benefit from both recent climate warming) and CO2 fertilization. In addition, treeline trees are rarely forced to restrict transpiration due to ample soil water availability (Tranquillini, 1979; Matyssek et al., 2009; Wieser et al., 2015). A stable iWUE suggests an increase of both carbon gain and leaf conductance for water vapor as also indicated by stable C and O isotope analysis and direct gas exchange assessments. Furthermore, iWUE may not change species composition at treeline in the Austrian Alps due to similar ecophysiological responses to climatic changes of the three sympatric study species. Our finding that growth of treeline associated conifers increases with slowly rising ambient CO2 concentration and warming may be relevant for assessing complex growth models with empirical data, finally leading to model improvements and better estimations of forest-climate feedbacks.
GW, RM, and TG conceived and designed the experiment. WO, AG, and ML performed the experiment. GW, WO, AG, and ML analyzed the data. GW, RM, and TG wrote the manuscript and WO and AG provided editorial advice.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
This work was supported by the Austrian Science Fund Project No. FWF P22206-B16 “Transpiration of conifers in contrasting environments.” Thanks to the Umweltbundesamt Austria for providing us with annual mean δ18O data for Mt. Patscherkofel.
Barbour, M. M., and Farquhar, G. D. (2000). relative humidity and ABA-induced variation in carbon and oxygen isotope ratios of cotton leaves. Plant Cell Environ. 23, 473–485. doi: 10.1046/j.1365-3040.2000.00575.x
Barbour, M. M., Fischer, R. A., Sayre, K. D., and Farquhar, G. D. (2000). Oxygen isotope ratio of leaf and grain material correlates with stomatal conductance and grain yield in irrigated wheat. Austr. J. Plant Physiol. 27, 625–637.
Barron-Gafford, G. A., Grieve, K. A., and Murthy, R. (2007). Leal-and stand-level responses of a forested mesocosm to independent manipulations of temperature and vapour pressure deficit. New Phytol. 174, 614–625. doi: 10.1111/j.1469-8137.2007.02035.x
Battipaglia, G., Saurer, M., Cherubini, P., Calfapietra, C., McCarthy, H. R., Norby, R. J., et al. (2013). Elevated CO2 increases tree-level intrinsic water use efficiency: insights from carbon and oxygen isotope analyses in tree rings across three forest FACE sites. New Phytol. 197, 544–554. doi: 10.1111/nph.12044
Benecke, U., Schulze, E.-D., Matyssek, R., and Havranek, W. M. (1981). Environmental control of CO2-assimilation and leaf conductance in Larix decidua Mill. I. a comparison of contrasting natural environmets. Oecologia 50, 54–61. doi: 10.1007/BF00378793
Borella, S., Leuenberger, M., Saurer, M., and Siegwolf, R. (1998). Reducing uncertainties in δ13C analysis of tree rings: pooling, milling and cellulose extraction. J. Geophys. Res. 103, 519–526. doi: 10.1029/98JD01169
Brendel, O., Iannetta, P. P. M., and Stewart, D. (2000). Arapid and simple method to isolate pure alpha-cellulose. Phytochem. Anal. 11, 7–10. doi: 10.1002/(SICI)1099-1565(200001/02)11:1<7::AID-PCA488>3.0.CO;2-U
Brugnoli, E., and Farquhar, G. D. (2000). “Photosynthetic fractionation of carbon isotopes,” in Photosynthesis: Physiology and Metabolism, Vol. 9, eds R. C. Leegood, T. D. Sharkey, and S. von Caemmerer (Dordrecht: Kluwer Academic Publishers), 399–434.
Bunn, A. G., Graumlich, L. J., and Urban, D. L. (2005). Trends in twentieth century tree growth at high elevations in the Sierra Nevada and White Mountains, USA. Holocene 15, 481–488. doi: 10.1191/0959683605hl827rp
Cernusak, L. A., Ubierna, N., Winter, K., Holtum, J. A. M., Marshall, J. D., and Farquhar, G. D. (2013). Environmental and physiological determinants of carbon isotope discrimination in terrestrial plants. New Phytol. 200, 950–965. doi: 10.1111/nph.12423
Danby, R. K., and Hik, D. S. (2007). Response of white spruce (Pinus glauca) to experimental warming at a subarctic alpine treeline. Glob. Change Biol. 13, 437–451. doi: 10.1111/j.1365-2486.2006.01302.x
Dawes, M. A., Hagedorn, F., Handa, I. T., Streit, K., Ekbald, Q. A., Rixen, C., et al. (2013). An alpine treeline in a carbon dioxide-rich world: synthesis of a nine-year free-air carbon dioxide enrichment study. Oecologia 171, 623–637. doi: 10.1007/s00442-012-2576-5
Duquesnay, A., Breda, N., Stievenard, M., and Dupouey, J. L. (1998). Changes of tree-ring δ13C and water-use efficiency of beech (Fagus sylvatica L.) in north-eastern France during the past century. Plant Cell Environ. 21, 565–572. doi: 10.1046/j.1365-3040.1998.00304.x
Farquhar, G. D., O’Leary, M. H., and Berry, J. A. (1982). On the relationship between carbon isotope discrimination and the intercellular carbon dioxide concentration in leaves. Aust. J. Plant. Physiol. 9, 121–137. doi: 10.1071/pp9820121
Flexas, J., Ribas-Carbo, M., Diaz-Espejo, A., Glamés, J., and Medrano, H. (2008). Mesophyll conductance to CO2: current knowledge and future prospects. Plant Cell Environ. 31, 602–621. doi: 10.1111/j.1365-3040.2007.01757.x
Frank, D., and Esper, J. (2005). Temperature reconstructions and comparisons with instrumental data from a tree-ring network for the European Alps. Int. J. Climatol. 25, 1437–1454. doi: 10.1002/joc.1210
Frank, D. C., Poulter, B., Saurer, M., Esper, J., Huntingford, C., Helle, G., et al. (2015). Water-use efficiency and transpiration across European forests during the Anthropocene. Nat. Clim. Change 5, 579–583. doi: 10.1038/NCLIMATE2614
Gessler, A., Brandes, E., Keitel, C., Boda, S., Kayler, Z., Granier, A., et al. (2013). The oxygen isotope enrichment of leaf-exported assimilates—does it always reflect lamina leaf water enrichment? New Phytol. 200, 144–157. doi: 10.1111/nph.12359
Grams, T. E. E., Kozovits, A. R., Häberle, K.-H., Matyssek, R., and Dawson, T. E. (2007). Combining δ13C and δ18O analyses to unravel competition, CO2 and O3 effects on the physiological performance of different-aged trees. Plant Cell Environ. 30, 1023–1034. doi: 10.1111/j.1365-3040.2007.01696.x
Gruber, A., Pirkebner, D., Oberhuber, W., and Wieser, G. (2011). Spatial and seasonal variations in mobile carbohydrates in Pinus cembra in the timberline ecotone of the Central Austrian Alps. Eur. J. Forest Res. 130, 173–179. doi: 10.1007/s10342-010-0419-7
Gruber, A., Wieser, G., and Oberhuber, A. (2009). Intra-annual dynamics in stem CO2 eﬄux in relation to cambial activity and xylem development in Pinus cembra. Tree Physiol. 29, 641–649. doi: 10.1093/treephys/tpp001
Jaggi, M., Saurer, M., Fuhrer, J., and Siegwolf, R. (2002). The relationship between the stable carbon isotope composition of needle bulk material, starch, and tree rings in Picea abies. Oecologia 131, 325–332. doi: 10.1007/s00442-002-0881-0
Kellomäki, S., and Wang, K.-Y. (1998). Sap flow in Scots pines growing under conditions of year-round carbon dioxide enrichment and temperature elevation. Plant Cell Environ. 21, 969–981. doi: 10.1046/j.1365-3040.1998.00352.x
Kress, A., Saurer, M., Siegwolf, R. T. W., Frank, D., Esper, J., and Bugmann, H. (2011). A 350 year drought reconstruction from Alpine tree ring stable isotopes. Global Biogeochem. Cycles 24. doi: 10.1029/2009GB003613
Levesque, M., Siegwolf, R., Saurer, M., Eilmann, B., Brang, P., Bugmann, H., et al. (2014). Drought response of five conifer species under contrasting water availability suggests high vulnerability of Norway spruce and European larch. Glob. Change Biol. 19, 3184–3199. doi: 10.1111/gcb.12268
Liu, X., Wang, W., Xu, G., Zeng, X., Wu, G., Zhang, X., et al. (2015). Tree growth and intrinsic water-use efficiency of inland riparuian forests in norhwestern China: evaluation via δ13C and δ18O analysis of tree rings. Tree Physiol. 34, 966–980. doi: 10.1093/treephys/tpu067
Loader, N. J., McCarrol, D., Gagen, M., Robertson, I., and Jalkanen, R. (2007). “Extracting climatic information from stable isotopes in tree rings,” in Stable Isotopes as Indicators of Ecological Change, eds T. E. Dawson and R. T. W. Siegwolf (London: Academic Press), 27–48.
Matyssek, R., Wieser, G., Patzner, K., Blaschke, H., and Häberle, K.-H. (2009). Transpiration of forest trees and stands at different altitude: consistencies rather than contrasts? Eur. J. Forest Res. 128, 579–596. doi: 10.1007/s10342-008-0243-5
McCarrol, D., and Pawellek, E. (2004). Stable carbon isotope ratios of Pinus sylvestris from northern Finland and the potential for extracting a climate signal from long term Fennoscandian chronologies. Holocene 11, 517–526. doi: 10.1191/095968301680223477
Oberhuber, W., Kofler, W., Pfeifer, K., Seeber, A., Gruber, A., and Wieser, G. (2008). Long-term changes in tree-ring-climate relationships at Mt. Patscherkofel (Tyrol, Austria) since the mid 1980s. Trees 22, 31–40. doi: 10.1007/s00468-007-0166-7
Penuelas, J., Canadell, J. G., and Ogaya, R. (2011). Increased water-use efficiency during the 20th century did not translate into enhanced tree growth. Global Ecol. Biogeogr. 20, 597–608. doi: 10.1111/j.1466-8238.2010.00608.x
Peterson, D. L., Arbaugh, M. J., Robinson, L. J., and Dederian, B. I. (1990). Growth trends of whitebark pine and lodgepole pine in a subalpine Sierra Nevada forest, California, USA. Arct. Alp. Res. 22, 233–243. doi: 10.2307/1551586
Pilcher, J. R. (1990). “Sample preparation, cross-dating and measurement,” in Methods of Dendrochronology: Applications in the Environmental Sciences, eds E. R. Cook and L. A. Kairiukstis (Dordrecht: Kluwer), 40–51.
Pisek, A., Larcher, W., Moser, W., and Pack, I. (1969). Kardinale Temperaturbereiche der Photosynthese und Grenztemperaturen des Lebens der Blätter verschiedener Spermatophyten. III. Temperaturabhängigkeit und optimaler Temperaturbereich der Netto-Photosynthese. Fora Ab. B 158, 608–630.
Pisek, A., Larcher, W., Vegis, A., and Napp-Zinn, K. (1973). “The normal temperature range,” in Temperature and Life, eds H. Precht, J. Christophersen, H. Hensel, and W. Larcher (Berlin: Springer), 102–194.
Poage, M. A., and Chamberlain, C. P. (2001). Empirical relationships between elevation and stable isotope composition of precipitation and surface waters: considerations for studies of paleoelevation change. Am. J. Sci. 301, 1–15. doi: 10.2475/ajs.301.1.1
Roden, J. S., and Farquhar, G. D. (2012). A controlled test of the dual-isotope approach for the interpretation of stable carbon and oxygen isotope ratio variation in tree rings. Tree Physiol. 32, 490–503. doi: 10.1093/treephys/tps019
Sarris, D., Siegwolf, R., and Körner, C. (2013). Inter and intra-annual stable carbon and oxygen isotope signals in response to drought in Mediterranean pines. Agric. Forest Meteorol. 168, 59–68. doi: 10.1016/j.agrformet.2012.08.007
Saurer, M., Cherubini, P., Reynolds-Henne, C. E., Treydte, K. S., Anderson, W. T., and Siegwolf, R. T. W. (2008). An investigation of the common signal in tree ring stable isotope chronologies at temperate sites. J. Geophys. Res. 113:G04035, doi: 10.1029/2008JG000689
Saurer, M., Spahni, R., Frank, D. C., Joos, F., Leuenberger, M., Loader, N. J., et al. (2014). Spatial variability and temporal trends in water-se efficiency of European forest. Glob. Change Biol. 20, 332–336. doi: 10.1111/gcb.12717
Scheidegger, Y., Saurer, M., Bahn, M., and Siegwolf, R. (2000). Linking stable isotopes with stomatal conductance and photosynthetic capacity: a conceptual model. Oecologia 125, 350–357. doi: 10.1007/s004420000466
Sidorova, O. V., Siegwolf, R. T. W., Saurer, M., Shashkin, A. V., Knorre, A. A., Prokushkin, A. S., et al. (2009). Do centennial tree-ring and stable isotope trends of Larix gmelinii (Rupr.) indicate increasing water shortage in the Siberian north? Oecologia 161, 825–835. doi: 10.1007/s00442-009-1411-0
Silva, L. C. R., and Anand, M. (2013). Probing for the influence of atmospheric CO2 and climate change on forest ecosystems across biomes. Glob. Ecol. Biogeogr. 22, 83–92. doi: 10.1111/j.1466-8238.2012.00783.x
Smidt, S., and Mutsch, F. (1993). “Messungen der nassen Freilanddeposition an alpinen Höhenprofilen,” in Proceedings of the International Symposium ‘Stoffeinträge aus der Atmosphäre und Waldbodenbe-Lastung in der Ländern der ARGE ALP und ALPEN ADRIA’. GSF-Report, Neuherberg, 21–29.
Sohn, J. A., Gebhardt, T., Ammer, C., Bauhus, J., Häberle, K.-H., Matyssek, R., et al. (2013). Mitigation of drought by thinning: short-term and long-term effects on growth and physiological performance of Norway spruce (Picea abies). Forest Ecol. Manag. 308, 188–197. doi: 10.1016/j.foreco.2013.07.048
Streit, K., Siegwolf, R. T. W., Hagedorn, F., Schaub, M., and Buchmann, N. (2014). Lack of photosynthetic or stomatal regulation after 9 years of elevated [CO2] and 4 years of soil warming in two conifer species at the alpine treeline. Plant Cell Environ. 37, 315–326. doi: 10.1111/pce.12197
Treydte, K., Boda, S., Graf Pannatier, E., Fonti, P., Frank, D., Ullrich, B., et al. (2014). Seasonal transfer of oxygen isotopes from precipitation and soil to the tree ring: source water versus needle water enrichment. New Phytol. 202, 772–783. doi: 10.111/nph.12741
Treydte, K. S., Schleser, G. H., Schweingruber, F. H., and Winiger, M. (2001). The climatic significance of δ13C in subalpine spruces (Lötschental. Swiss Alps). Tellus Ser. B 55, 593–611. doi: 10.1034/j.1600-0889.2001.530505.x
Van Herk, I. G., Gower, S. T., Bronson, D. R., and Tanner, M. S. (2011). Effect of climate warming on canopy water dynamics of a boreal black spruce plantation. Can. J. For. Res. 41, 217–227. doi: 10.1139/X10-196
Weigt, R. B., Bräunlich, S., Zimmermann, L., Saurer, M., Grams, T. E. E., Dietrich, H.-P., et al. (2015). Comparison of δ13C and δ18O values between tree-ring whole wood and cellulose in five species growing under two different site conditions. Rapid Commun. Mass Spectr. 29, 2233–2244. doi: 10.1002/rcm.7388
Werner, C., Schnyder, H., Cuntz, M., Keitel, C., Zeeman, M. J., Dawson, T. E., et al. (2012). Progress and challenges in using stable isotopes to trace plant carbon and water relations across scales. Biogeosciences 9, 3083–3111. doi: 10.5194/bg-9-3083-2012
Wieser, G., Gigele, T., and Pausch, H. (2005). Seasonal and spatial variation of woody tissue respiration in a Pinus cembra tree at the alpine timberline in the Central European Alps. Eur. J. For. Res. 124, 1–8. doi: 10.1007/s10342-004-0050-6
Wieser, G., Grams, T. E. E., Matyssek, R., Oberhuber, W., and Gruber, A. (2015). Soil warming increased whole-tree water use of Pinus cembra at the treeline in the Central Tyrolean Alps. Tree Physiol. 35, 279–288. doi: 10.1093/treephys/tvp009
Wieser, G., Gruber, A., and Oberhuber, W. (2014). Sap flow characteristics and whole-tree water use of Pinus cembra across the treeline ecotone in the central Tyrolean Alps. Eur. J. For. Res. 133, 287–295. doi: 10.1007/s10342-013-0760-8
Wieser, G., Matyssek, R., Luzian, R., Zwerger, P., Pindur, P., Oberhuber, W., et al. (2009). Effects of atmospheric and climate change at the timberline of the Central European Alps. Ann. For. Sci. 66:402. doi: 10.1051/forest/2009023
Wieser, G., Oberhuber, W., Walder, L., Spieler, D., and Gruber, A. (2010). Photosynthetic temperature adaptation of Pinus cembra within the timberline ecotone of the Central Austrian Alps. Ann. For. Sci. 67:201. doi: 10.1051/forest/2009094
Wu, G., Liu, X., Chen, T., Xu, G., Wang, W., Zeng, X., et al. (2015). Elevation-dependent variations of tree growth and intrinsic water-use efficiency in Schrenk spruce (Picea schrenkiana) in the westen Tianshan Mountain, China. Front. Plant Sci. 6:309. doi: 10.3389/fpls.2015.00309
Keywords: stable isotopes, intrinsic water use efficiency, tree growth, climate change, treeline, Central Alps
Citation: Wieser G, Oberhuber W, Gruber A, Leo M, Matyssek R and Grams TEE (2016) Stable Water Use Efficiency under Climate Change of Three Sympatric Conifer Species at the Alpine Treeline. Front. Plant Sci. 7:799. doi: 10.3389/fpls.2016.00799
Received: 22 April 2016; Accepted: 22 May 2016;
Published: 07 June 2016.
Edited by:Boris Rewald, University of Natural Resources and Life Sciences, Vienna, Austria
Reviewed by:Eryuan Liang, Chinese Academy of Sciences, China
Cristina Nabais, University of Coimbra, Portugal
Copyright © 2016 Wieser, Oberhuber, Gruber, Leo, Matyssek and Grams. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Gerhard Wieser, email@example.com