Nitrogen Addition and Understory Removal but Not Soil Warming Increased Radial Growth of Pinus cembra at Treeline in the Central Austrian Alps.

Beside low temperatures, limited tree growth at the alpine treeline may also be attributed to a lack of available soil nutrients and competition with understory vegetation. Although intra-annual stem growth of Pinus cembra has been studied intensively at the alpine treeline, the responses of radial growth to soil warming, soil fertilization, and below ground competition awaits clarification. In this study we quantified the effects of nitrogen (N) fertilization, soil warming, and understory removal on stem radial growth of P. cembra at treeline. Soil warming was achieved by roofing the forest floor with a transparent polyvinyl skin, while understory competition was prevented by shading the forest floor with a non-transparent foil around six trees each. Six trees received N- fertilization and six other trees served as controls. Stem growth was monitored with band dendrometers during the growing seasons 2012-2014. Our 3 years experiment showed that soil warming had no considerable effect on radial growth. Though understory removal through shading was accompanied by root-zone cooling, understory removal as well as N fertilization led to a significant increase in radial growth. Hardly affected was tree root biomass, while N-fertilization and understory removal significantly increased in 100-needle surface area and 100-needle dry mass, implying a higher amount of N stored in needles. Overall, our results demonstrate that beside low temperatures, tree growth at cold-climate boundaries may also be limited by root competition for nutrients between trees and understory vegetation. We conclude that tree understory interactions may also control treeline dynamics in a future changing environment.

Introduction

Although the alpine treeline has attracted the interest of researchers for more than one century (Brockmann-Jerosch, 1919; Däniker, 1923; Wardle, 1974; Tranquillini, 1979; Körner, 1998, 2012; Holtmeier, 2003; Wieser and Tausz, 2007; Piper et al., 2016), the causes to treeline formation are still under debate. Nevertheless, it has been hypothesized that low temperatures limit growth processes in meristematic tissues (growth limitation hypothesis; Körner, 1998, 2012; Hoch and Körner, 2003). A world-wide survey across natural high elevation treelines indicates that growing season mean root-zone temperature constrains tree growth in temperature limited ecosystems (Sveinbjörnsson, 2000; Körner and Hoch, 2006). In a global survey, Körner and Paulsen (2004) found that a growing season mean soil temperature of 6.4 ± 0.7°C in 10 cm soil depths matches the upper elevational limit of tree growth. Moreover, a growth decline at the alpine treeline may also be attributed to a lack of available soil nutrients (Tranquillini, 1979; Sveinbjörnsson et al., 1992), especially nitrogen (Hoch, 2013) and competition with understory vegetation (Elliott, 2011; Grau et al., 2012; Liang et al., 2016).

Given that alpine treelines are strongly temperature limited, they are considered to be sensitive to climate warming (Walther et al., 2005; Holtmeier and Broll, 2007; Wieser et al., 2009; Liang et al., 2016; Camarero et al., 2017). So far, only a few studies evaluated the effect of soil warming on tree growth in boreal forests (Strömgren and Linder, 2002; Dao et al., 2015) and at treeline in the Swiss (Dawes et al., 2015) and French Alps (Loranger et al., 2016) with results ranging from a strong growth stimulation to no stimulation in growth. In a pilot study at treeline in the Austrian Alps, Gruber et al. (2010) investigated the effect of root-zone warming and cooling on radial growth of Pinus cembra. Although not statistically significant, results of this study indicated that P. cembra responded to soil warming with an increase and to soil cooling with a decline in stem radial growth, when compared to control trees with soil temperature left unmanipulated. Moreover, differences in stem growth with respect to soil temperature may also be attributed to varying soil nutrient contents with respect to micro topography (Anschlag et al., 2008).

Beside low temperatures, the termination of tree growth their upper distribution limit may also be attributed to a low availability of soil nutrients (Sveinbjörnsson et al., 1992) and hence also to an insufficient nutrient uptake (Susiluoto et al., 2010). Particularly total and available nitrogen (N), which is a key nutrient, is known to limit plant productivity in numerous terrestrial ecosystems (LeBauer and Treseder, 2008). This also holds for the treeline ecotone (Tranquillini, 1979) where low soil temperatures limit N mineralization and decomposition (Nadelhoffer et al., 1991; Rustad et al., 2001; Melillo et al., 2002). Yet, soil fertilization studies in boreal treelines indicated enhanced tree growth after soil fertilization (Weih, 2000; Weih and Karlsson, 2001; Susiluoto et al., 2010).

Beside the two abiotic factors, low temperature and low nutrient availability, root competition for nutrients may also have a noticeable influence on tree growth (Nilsson and Wardle, 2005; Matsushima and Chang, 2006; Elliott et al., 2015; Du et al., 2016). There are some data indicating that the absence of understory (i.e., the removal of below-ground competition) may enhance tree (Platt et al., 2004; Song et al., 2016) and seedling growth (Okano and Bert-Hane, 2015) in subalpine forests and the treeline ecotone, respectively. To our knowledge, however, the role of root competition between adult trees and understory vegetation in regulating tree growth at treeline in the Central European Alps received no attention.

Although the intra-annual stem radial increment of P. cembra has been studied intensively at the alpine treeline (Loris, 1981; Rossi et al., 2006, 2007; Gruber et al., 2009a,b), the responses of radial growth to soil warming, soil fertilization, and below ground competition, respectively, still awaits clarification for conifers at the alpine treeline. In this unique study, we investigated how soil temperature, N fertilization, and understory removal influences stem radial growth of P. cembra at treeline. Additionally, we also estimated root biomass, needle and foliar nutrient concentrations as well as specific leaf area. We hypothesized that: (1) soil warming, (2) understory removal, and (3) N fertilization will enhance radial growth of P. cembra at treeline in the Central Tyrolean Alps. Soil warming and understory removal was incited by root-zone roofing throughout three consecutive growing seasons, while continuously monitoring intra-annual changes in stem radius with electronic band dendrometers (Gruber et al., 2010; Oberhuber, 2017). In addition, some trees received a soil nitrogen fertilization. Findings are expected to contribute to an increased understanding of the importance of root competition and soil N fertilization on tree growth at treeline in the Central Austrian Alps, where mean annual air temperature increased by 0.50°C per decade, during the past 35 years (Wieser et al., 2016).

Materials and Methods

Study Site and Experimental Design

The study was carried out in a south exposed P. cembra afforestation at treeline above Haggen near St. Sigmund in the Sellrain Valley, Tyrol, Austria (“Haggener Sonnberg,” 47°12′42″N, 11°05′04″E, 2150 m a.s.l.). Slope angel and aspect were 25° and SSW, respectively (Kronfuss, 1997). The long term mean annual temperature (1975–1994) at a nearby weather station at 1800 m a.s.l. (Kronfuss, 1997) was 3.2°C and the mean annual precipitation was 909 mm. For the same period the growing season (May through September) had a mean air temperature of 8.5°C and a mean precipitation of 537 mm. The soil at our study site is a podzolic cambisol (Neuwinger, 1972; World Base for soil Resources classification, FAO, 2008) which derived from gneisses and mica schist bedrock (Kronfuss and Havranek, 1999). The sandy loam subsoil was covered by a ca. 5 cm thick humus layer (Wieser et al., 2015). Hydraulic field capacity at -0.033 MPa (sensu Blume et al., 2010) of the top sub soil (0–25 cm) is 25% volume and the top soil is enriched by 8% of organic matter (Wieser et al., 2015).

The stand formed a sparsely open canopy permitting a dense understory of herbaceous species together with some dwarf shrubs. During the study period (2011–2014) the trees were ≈25 years old, with a stem diameter at breast height of 7.1 ± 1.2 cm, and an average height of 3.4 ± 0.4 m. Trees selected for the experiment were separated at least by a distance of 3–5 m. In early summer 2011 we established 24 quadratic 16 m2 plots with one P. cembra tree in the center, resulting in six plots each of the four experimental groups: (1) nitrogen (N) fertilization (N treatment), (2) an increase in soil temperature (warming treatment), (3) elimination of understory vegetation (understory removal treatment), and (4) controls, for with soil temperature left un-manipulated, the understory vegetation was present (competition), and trees were not fertilized (control treatment). In order to facilitate the logistics of soil temperature and soil moisture monitoring (see below), three adjacent plots per treatment were assembled into a block.

The nitrogen plots were fertilized following the recommendations of Kilian and Müller (1994). The plots were fertilized twice, in spring 2012 and 2013, with 10 g m-2 calcium ammonium nitrate (NAC 27 N; Borealis L.A.T., AT) containing 27% nitrogen with equal parts of NO3-N and NH4 and 12% CaO; the latter counteracting potential effects of N acidification. Soil warming was achieved by roofing the forest floor with a heat trapping 0.5-mm thick transparent polyvinyl skin (Wieser et al., 2015). In the no understory plots, herbaceous species and dwarf shrubs were completely removed in July 2011. This condition was maintained by shading the forest floor with a non-transparent foil. The respective foils were fixed on frames 15 cm above ground. The sides were open for allowing air circulation. Surface and slope run-off of water was allowed to penetrate the soil during and after rainfall (in total 67% of precipitation; Neuwinger et al., 1988; Wieser et al., 2015). Soil warming and understory removal through shading operated from end July throughout October 2011, and continued during the snow-free periods (May–October) of 2012, 2013, and 2014.

As expected, understory removal through shading was accompanied by root-zone cooling. Additionally, there is also evidence that artificial soil cooling by roofing the forest floor may have negative effects on stem radial growth of P. cembra growing in the krummholz (crippled trees) limit at 2180 m a.s.l. on Mt. Patscherkofel south of Innsbruck, Austria (Gruber et al., 2010). Therefore, in order to discriminate between potential antagonistic effects of understory removal and soil cooling on radial growth, we also removed the understory in all the control and warmed plots in early April 2014 and covered the ground with a permeable weed fleece to avoid any re-growth of the understory vegetation throughout the entire growing season of 2014. Detailed information on potential understory competition with respect to treatment is given in Table .

Environmental Measurements and Dendrometer Records

Air temperature (Tair) and relative humidity (CS215 Temperature and Relative Humidity Sensor), solar radiation (SP1110 Pyranometer Sensor), wind velocity (A100R Anemometer) and precipitation (ARG100 Rain Gauge, all sensors Campbell Scientific, Shepshed, United Kingdom) were measured continuously at 2 m height at the study site. In order to examine differences in the seasonal course of soil temperature (Tsoil) and volumetric soil water content (𝜃) between control, warming, and understory removal blocks, six soil temperature probes (T 107 Temperature Probe, Campbell Scientific, Shepshed, United Kingdom) and two soil moisture sensors (EC 5 Soil Moisture Sensor, Decagon Devices Inc. Pullman, WA, United States) were installed in each block close to the trees used for dendrometer records. While 𝜃 was measured at 10 cm soil depth solely, Tsoil was measured at 5, 10, and 20 cm soil depth (two probes per depth and block). To evaluate a potential influence of forest-floor roofing on stem temperature (Tstem), a T 107 Temperature Probe (Campbell Scientific, Shepshed, United Kingdom) was mounted at 50 cm stem height on the north facing side of all the trees in control, warmed and understory removal treatment, respectively. All the environmental data were recorded with two CR1000 data loggers (Campbell Scientific, Shepshed, United Kingdom) programmed to record 30-min averages of measurements taken every minute.

Intra-annual changes in stem radius of all the selected study trees were continually monitored using electronic band dendrometers with automatic temperature compensation (DC2, Ecomatik, Dachau, Germany) installed 0.5 m aboveground in early May 2011 to record intra-annual radial growth throughout the snow free periods (April–October) of 2012, 2013, and 2014. The measuring cable consisted of Invar-steel with a thermal expansion coefficient <1.4 × 10-6 K. Dead outermost layers (periderm) of the bark were slightly removed to reduce the influence of hygroscopic swelling and shrinkage of the bark on dendrometer traces and to ensure close contact with the stem (cf. Gruber et al., 2009a,b; Oberhuber, 2017). Data were recorded with a CR1000 data logger (Campbell Scientific, Shepshed, United Kingdom) programmed to record 30-min averages of measurements taken every minute. Circumference variations derived from band dendrometers were transformed to radial variations, and daily stem radius variations were determined by calculating the difference between mean values of two consecutive days (“daily mean approach,” Deslauriers et al., 2007). As for dendrometers no unambiguous date of growth onset can be determined (Downes et al., 1999; Deslauriers et al., 2003; Gruber et al., 2009b), we defined the onset of radial growth as the date when radial increment permanently exceeded the initial value after snow melt in May. To separate daily patterns of water movement from irreversible expansion growth, we applied a Gompertz function for describing the long-term development of radial growth over an entire growing season according to Rossi et al. (2006):

where y is the cumulative sum of growth, A is the upper asymptote of the total annual radial growth, β is the x-axis placement parameter, κ is the rate of change parameter, and t is the time computed in Julian days. The inflection point (Ip) corresponding to the maximum value of the radial growth rate was calculated as Ip = β/κ (Rossi et al., 2006). Finally, as dendrometer records are affected by stem water status, the end of stem radial growth was considered when 90% of A given by the Gompertz function was reached (Gruber et al., 2009b).

Additional Measurements

Tree characteristics including specific leaf area (SLA), 100 needle dry weight, and foliar nutrient concentrations of all the study trees, as well as root biomass at the plot level were determined in late fall 2013 (November 18, 2013). SLA (cm2 g-1) was calculated from measured needle dry weight and measured projected needle surface area. For nutrient analyses the needles were dried to constant weight at 60°C, ground and stored dry before analysis. The concentrations of nitrogen, phosphor and potassium were estimated according to EPA 3052[1]. For assessing root biomass mass (including tree and understory vegetation roots) at the plot level we sampled six soil cores (volume 125 cm3) per treatment in 0–10 and in 10–20 cm soil depth, respectively. Roots were rinsed with tap water in a sieving cascade to remove all soil particles while minimizing the loss of fine roots. Subsequently, the roots were divided into fine (diameter ≤ 1 mm) and coarse roots (diameter > 1 mm) and dried to constant weight at 75°C before determining dry weight. Finally, increment cores (5 mm in diameter) were taken at sensor height in fall 2014 for estimating radial increment in 2010, the year preceding the experiment.

Data Analysis

As daily mean Tsoil and 𝜃 values in the two control, the two warmed, and the two understory removal blocks, as well as Tstem of all the trees under study did not differ significantly [all P-values > 0.1; one-way analysis of variance (ANOVA)], we used repeated-measures ANOVA to test for differences in Tsoil, 𝜃, and Tstem between control, warmed, and understory removal blocks. Differences in the overall mean parameters of the Gompertz functions for modeling intra-annual radial growth (upper asymptote, inflection point, rate of change parameter, time when 90% of increment were produced) between the control, the warming, the understory removal, and the fertilization treatment during the growing seasons 2012, 2013, and 2014 were tested for significance using analysis of variance (ANOVA) and multiple comparison Tukey’s HSD tests.

We used binary logistic regressions (logit models; SPSS for Windows) to determine the probability of radial growth being active at a given air and soil temperature. Values of daily radial increment where binary coded as: no increment (value 0) or increment (value 1) during the period May 1–September 30. The model was fitted for each tree, treatment and year with the respective temperature series; i.e., daily mean air and soil temperature in 10 cm soil depth. Temperature thresholds were calculated when the probability of radial increment was 0.5 (cf. Rossi et al., 2007). Fitting verification included χ2 of the likelihood ratio, Wald’s χ2 for regression parameter and goodness of fit, and Hosmer–Lemeshow Ĉ for eventual lack of fit. Finally, thermal thresholds were averaged over the three snow free periods 2012, 2013, and 2014 and compared between treatments (control, warmed, understory removal ≈ cooling) using ANOVA.

ANOVA was also used for testing differences in SLA, needle area, needle mass, needle nutrient concentrations, and total fine root mass between treatments. A probability level of P < 0.05 was considered as statistically significant and P-values ≥ 0.05 and < 0.1 as marginally significant. Statistical analyses were made with the IBM SPSS Statistics 21 (IBM, NY, United States).

Results

Environmental Conditions

Growing season (May–September) mean Tair (Figure ) was 9.0°C in 2012, 8.1°C in 2013, and 7.3°C in 2014. Daily mean relative humidity varied between 50% during cloudless warm days and close to 100% on overcast rainy days (Figure ). Stem temperature at 50 cm stem height was not considerably altered by the experimental setup (data not shown; P-values from repeated measure analysis: P = 0.215, P = 0.418, and P = 0.390 for 2012, 2013, and 2014, respectively).

(A) Seasonal course of daily mean air temperature and daily mean soil temperatures in 10 cm soil depth in control, warming and understory removal blocks. (B) Relative air humidity. (C) Daily sum of precipitation and daily mean soil water content in 10 cm soil depth in control, warming and understory removal blocks from April to October 2012, 2013, and 2014. Arrows indicate roof closure.

Independent of treatment, Tsoil generally followed seasonal trends in Tair (Figure ). Warming caused Tsoil to average +1.5 ± 0.5, +1.3 ± 0.2, and +1.0 ± 0.2°C above the corresponding control levels at 5, 10, and 20 cm soil depth, respectively (all P < 0.05; Table ). In the understory removal blocks by contrast, shading caused Tsoil to be at average -3.0 ± 0.8, -2.6 ± 0.5, and -2.4 ± 0.4°C below the corresponding control levels at 5, 10, and 20 cm soil depth, respectively (all P < 0.05; Table ).

Precipitation during the growing seasons 2012, 2013, and 2014 amounted 940, 806, and 785 mm, respectively. Due to frequent precipitation over the entire growing seasons (Figure ) soil water content (𝜃) in 10 cm soil varied between 44 and 19% vol. (Figure ). Average 𝜃 over the growing seasons 2012, 2013, and 2014 did not differ significantly (all P > 0.64) between control (30.8 ± 0.3% vol.), warmed (31.1 ± 1.6% vol.) and understory removal (shaded) blocks (30.1 ± 1.7% vol.), respectively (Table ).

Treatment Effects on Growth

In 2010, the year preceding the experiment, annual radial stem increment did not differ significantly (all P > 0.45) between the treatments and averaged 3.5 ± 0.7, 3.4 ± 0.5, 3.9 ± 0.5, and 3.8 ± 0.7 mm in the control, the warming, the understory removal, and the fertilization treatment, respectively. The corresponding values for 2011 were 3.2 ± 0.7, 3.1 ± 1.0, 3.5 ± 0.4, and 3.4 ± 0.6 mm, respectively. Independent of treatment and year, radial growth started in the 2nd week of May (doy 128–133; Figure ). Based on the calculated Gompertz functions, maximum daily radial growth peaked in June around summer solstice (doy 167–173) and radial growth ended in early August (doy 218–227), with no treatment effects (Table ). Treatment, however, affected total annual stem radial increment (Table ), which corresponds to the upper asymptote of the Gompertz function (Figure ). Fertilization caused annual radial growth to increase significantly by 28 and 48% in 2012 and 2013, respectively, (both P < 0.004) above the level of the control treatment (Figure and Table ). Understory removal caused radial growth to increase marginally significant on average by 11% (P = 0.08) in 2012, and significantly on average by 22% (P = 0.04) in 2013 above the level of control trees (Figure and Table ). However, it has to be noted that understory removal through shading was accompanied by significant cooling of the rooting-zone (cf. Figure and Table ). Warmed and control trees by contrast, resembled each other with respect to radial growth (Figure ), in the absence of a warming effect on annual stem radial increment in 2012 and 2013 (Table ). Removing the understory in 2014 in all the control and warmed plots also failed to find significant temperature effects (both P > 0.35) on total annual stem radial increment between the control, the warming and the understory removal treatment (Table ).

Time series of mean daily Pinus cembra dendrometer records and standard deviations of the control, warming, understory removal and the N during growing seasons 2012, 2013, and 2014. Circles denote mean upper asymptotes ± SD of the Gompertz functions of each treatment.

Threshold Temperatures

The calculated threshold temperatures at which radial growth had a 0.5 probability of being active were averaged over the three growing seasons 2012, 2013, and 2014. Irrespective of treatment, the mean air temperature at which radial growth of P. cembra was active in control, warmed and understory removal plots varied between 5.1 and 5.9°C (Table ; all P > 0.05). The mean soil temperatures at which there was a 0.5 probability of radial growth were in general higher than the air temperature thresholds, being 7.8 ± 0.4, 9.1 ± 0.6 and 6.0 ± 0.6°C in the control, the warmed, and the understory removal treatment, respectively (Table ; all P < 0.05).

Specific Leaf Area, Foliar Nutrient Concentrations of the Study Trees and Plot-Level Fine Root Biomass

Treatment did not influence SLA (P = 0.98). At the end of the growing season 2013 current-year needle SLA averaged 43.5 ± 1.4 cm2 g-1 across all treatments. Current year 100-needle surface area and the corresponding 100-needle dry weight, however, were significantly higher in the understory removal and fertilization treatment than in the warming and the control treatment (Figure ). Treatment had no effect in foliar nutrient concentrations. Averaged across all the four treatments nitrogen, phosphor and potassium concentrations of current-year needles were 18.3 ± 1.3, 1.9 ± 01, and 5.7 ± 0.4 mg g-1, respectively, at the end of the growing season 2013.

Box plots showing the median, lower (25%) and upper (75%) quartile and minimum and maximum of 100 needle area and 100 needle mass of Pinus cembra trees in control (c), warmed (w), understory removal (ur) and N fertilized (N) plots at the end of the growing season 2013. Different letters indicate significant differences at P < 0.05.

The impact of understory removal is also reflected at the plot level. After three growing seasons total fine-root biomass mass (including tree and understory vegetation roots) in 0–10 cm soil depth was significantly lower in understory removal plots as compared to control and warmed (both P < 0.05) plots (Table ). Although statistically not significant, total fine-root biomass mass in fertilized plots was considerably reduced when compared to warmed, and control plots, but considerably higher than in understory removal plots (Table ). Treatment, however, had no effect on tree root biomass in 0–10 cm soil depth as well as on total fine-root biomass in the 10–20 cm soil layer (Table ).

Discussion

In this study, we tested the hypothesis that stem radial growth of P. cembra at treeline in the Central Tyrolean Alps is enhanced by soil warming, understory removal, as well as N-fertilization. Our experimental approach was appropriate to manipulate Tsoil, competition and N availability, enabling the clarification of radial growth in situ under a wide range of conditions. Although the roofing prevented 33% of growing season precipitation to reach the soil, treatment differences in daily mean 𝜃 between control, warmed and understory removal blocks stayed within the typical variation of at the study site (Neuwinger, 1972) which confirmed that the employed roofing system did not prevent any shortage in soil water availability (Figure and Table ). This is further corroborated by sap flow measurements, where roofing caused sap flow density of P. cembra to increase above the levels of control trees, whereas leaf level net CO2 uptake, conductance for water vapor, and water-use efficiency stayed unchanged (Wieser et al., 2015). Therefore, we assumed that our roofing induced 33% reduction of growing season precipitation reaching the soil did not considerably influence radial growth as a potential “side-effect” due to soil roofing.

After 3 years of treatment, our results indicate that soil warming had no considerable effect on radial growth, which allows hypothesis (1) to be rejected. Although understory removal through shading was accompanied by root-zone cooling (Figure and Table ), understory removal as well as soil N-fertilization significantly increased radial growth, confirming hypothesis (2) and (3). Despite the fact that growing season stem radial growth of trees in the understory removal and the N addition treatment was significantly larger than in the control and warming treatment (Figure and Table ), we observed no considerable differences in the dynamics of tree growth (Figure ).

Similar to our findings, Gruber et al. (2010) also failed to find a significant increase in radial growth of P. cembra in the kampfzone, i.e., the upper border line of the treeline ecotone (Körner, 2012), on Mt. Patscherkofel south of Innsbruck (2.180 m a.s.l.) after 1 year of 0.5°C soil warming. An increase in soil temperature of 4°C and as a consequence a 1-week earlier snow melt also did not affect growth and cambial phenology of Picea mariana in a boreal forest in Canada (Lupi et al., 2012). Dao et al. (2015) also failed to find substantial effects on xylem phenology and cell production of mature black spruce after 6 years of soil warming by 4°C in a boreal forest. As in our study, 6 years of soil warming by 3.2°C in 10 cm soil depth at treeline in the Swiss Alps (2.180 m a.s.l.) also resulted in no considerable increase in the growth of P. uncinata, and Larix decidua during the first 2 years of the experiment (Hagedorn et al., 2010), while in the following 4 years slightly improved shoot growth persisted in P. uncinata, but diminished in L. decidua (Dawes et al., 2017). These observed differences in the response to soil warming may be attributed to differences in the rooting depth between pine and larch, because coarse roots of larch can penetrate into deeper soil layers when compared to pine (Kutschera and Lichtenegger, 2002), and thus may probably have not experienced the warming treatment (Dawes et al., 2015). Strömgren and Linder (2002) reported an increase in stem production of boreal Picea abies trees in a 7-year study combining soil warming and nutrition.

The absence of significant growth stimulation due to soil warming in our study suggests that growth stimulation due to increasing soil temperatures may not be expected for P. cembra trees at treeline in the Central European Alps in the future. However, a potential effect of soil warming may occur under a longer period of experimentation (Dao et al., 2015), as there is evidence that P. cembra exhibits determinate shoot growth, in any year growth is primarily impelled by the bud formation occurring in the previous year (Tranquillini, 1979). Additionally, in alpine ecosystems, growth response to warming is mediated through a time lag by one or more years after initiation of warming (Danby and Hik, 2007; Hagedorn et al., 2010). Responses during the extraordinary warm year of 2003 support our result, as the growth of P. cembra at the treeline in the Austrian Alps was hardly affected by a 4°C warmer summer (Oberhuber et al., 2008).

As understory removal by shading was accompanied by a significant cooling of the rooting-zone (Figure and Table ), we were able to calculate threshold temperatures above which significant tree growth occurs (for a review see Körner, 2006). Körner and Paulsen (2004) found that a mean annual air temperature and 10 cm soil temperature threshold of 5.6–8.5°C and 6.4°C, respectively, define the upper elevational limit of tree growth worldwide. For P. cembra in the Italian Alps, Rossi et al. (2007) estimated a daily threshold mean air temperature when radial growth was active of 5.8–8.6°C. In our 3-year study period radial growth was active when mean air temperature was above 5.1–5.9°C (Table ). Hence, our results support the existence of a comparable air temperature threshold for radial growth in P. cembra at the alpine treeline. Soil temperature thresholds by contrast differed significantly between treatments and ranged from 6.0°C in trees experiencing understory removal (≈cooling) to 9.1°C in trees experiencing soil warming (Table ), indicating that soil temperature is not the main limiting factor for radial growth for P. cembra at treeline (Rossi et al., 2007; Lupi et al., 2012).

Beside low temperatures, low soil nutrient availability may also affect tree growth at high altitude and latitude (Tranquillini, 1979; Jarvis and Linder, 2000; Rossi et al., 2011) because low soil temperatures limit mineralization, decomposition, and N-cycling (Nadelhoffer et al., 1991; Rustad et al., 2001; Melillo et al., 2002). Thus, on the long term future climate warming may affect wood production, although no short term effect could be observed (Vaganov et al., 1999; Jarvis and Linder, 2000; Rustad et al., 2001; Rossi et al., 2007). Additionally, there is also evidence that fertilization effects may increase under conditions of experimental soil warming as shown for seedlings (Weih, 2000; Hoch, 2013) and mature forest trees (Rustad et al., 2001; Melillo et al., 2011; Lupi et al., 2012; but see Dao et al., 2015). In their soil heating study, Strömgren and Linder (2002) attributed the observed growth stimulation to an increase in N mineralization. The enhanced growth of P. uncinata at treeline in the Swiss Alps during the first 3 years of soil warming has also been attributed to a warming induced increase in soil N mineralization (Dawes et al., 2017). However, this was not investigated in our study.

As shown previously for boreal treeline sites (Sveinbjörnsson et al., 1992; Weih, 2000; Weih and Karlsson, 2001; Susiluoto et al., 2010), N-fertilization significantly increased radial growth of P. cembra at our study site (Figure and Table ). Similar to N-fertilization, understory removal, and thus removal of root competition, also significantly enhanced stem growth of P. cembra (Figure and Table ). Removal of below-ground competition also enhanced tree (Platt et al., 2004; Song et al., 2016) and seedling growth (Okano and Bert-Hane, 2015) in subalpine forests and in the treeline ecotone, respectively. In addition, there is also evidence that the existence of understory can limit seedling and tree establishment above the current treeline (Smith et al., 2009; Batllori et al., 2010; Elliott, 2011; Grau et al., 2012). Assessment of long term changes in species interaction, however, are still a matter of debate (Matsushima and Chang, 2006; Liang et al., 2016; Camarero et al., 2017), and await clarification for the Central European Alps.

Nevertheless, beside enhanced stem growth, N-fertilization, as well as understory removal led to a significant increase in 100-needle surface area and 100-needle dry mass (Figure ) which may be attributed to a considerable increase in N uptake, although needle N concentration did not vary between treatments. Because a greater needle size with equal N concentrations implies a higher amount of N stored in needles, N uptake seemed indeed to be a limiting factor at our study site as trees did profit from N addition and understory removal both, at the stem- and the canopy level. Moreover, at the plot level, understory removal significantly reduced total (=tree and understory vegetation) fine-root biomass in 0–10 cm soil depth, while tree root biomass was hardly affected when compared to control plots (Table ). Thus our findings suggest that at treeline in the Central Austrian Alps root competition for nutrients as well as low soil nutrient availability, rather than soil warming limits tree growth of P. cembra.

Conclusion

Results of this study on the effect of soil warming, N fertilization, and understory removal on stem radial growth of P. cembra at treeline clearly showed that warming failed to induce enhanced stem growth, which allowed hypothesis (1) to be rejected. Understory removal as well as soil N-fertilization by contrast, significantly increased radial growth, confirming hypothesis (2) and (3). As hypothesis (1) was rejected, results of this study strongly suggest that beside low air temperatures, the termination of tree growth at cold-climate tree life boundaries may also be attributed to root competition for nutrients between trees and understory vegetation (Nilsson and Wardle, 2005; Matsushima and Chang, 2006; Elliott et al., 2015; Du et al., 2016). Finally, the presented results underscore the importance of including both temperature manipulation and species interactions in future studies on tree growth within the treeline ecotone.

Author Contributions

AG, WO, and GW conceived and designed the experiments and analyzed the data. AG performed the experiments. GW, wrote the manuscript. 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.

Table 1

Understory characteristics in the control, the warming, the understory removal, and the N fertilization treatment during the growing seasons 2012, 2013, and 2014.

Treatment	2012	2013	2014
Control	Intact understory	Intact understory	No understory
Warming	Intact understory	Intact understory	No understory
Understory removal	No understory	No understory	No understory
N treatment	Intact understory	Intact understory

Table 2

Seasonal average soil temperature (Tsoil; °C) at 5, 10, and 20 cm soil depth and volumetric soil water content (𝜃) at 10 cm soil depth, in control, warmed and understory removal blocks for the periods May 1– September 30, 2012, 2013, and 2014.

Year	Treatment	Tsoil 5 cm	Tsoil 10 cm	Tsoil 20 cm	𝜃 10 cm
2012	Control	11.2 ± 0.2a	10.7 ± 0.2a	10.0 ± 0.2a	30.5 ± 0.2a
	Warming	12.4 ± 0.2b	11.9 ± 0.2b	11.1 ± 0.2b	30.0 ± 0.1a
	Understory removal	9.1 ± 0.2c	9.1 ± 0.2c	8.4 ± 0.2c	30.5 ± 0.2a
2013	Control	10.1 ± 0.3a	9.8 ± 0.3a	9.4 ± 0.2a	31.0 ± 0.2a
	Warming	12.1 ± 0.3b	11.3 ± 0.3b	10.5 ± 0.3b	30.5 ± 0.3a
	Understory removal	8.6 ± 0.3c	8.3 ± 0.2c	7.9 ± 0.2c	28.3 ± 0.4a
2014	Control	9.3 ± 0.2a	9.1 ± 0.2a	8.9 ± 0.2a	31.0 ± 0.2a
	Warming	10.5 ± 0.2b	10.2 ± 0.2b	9.7 ± 0.2b	32.9 ± 0.1a
	Understory removal	8.4 ± 0.2c	8.1 ± 0.2c	7.8 ± 0.2c	31.6 ± 0.1a

Table 3

Parameters of the Gompertz function (Mean ± SD) for intra-annual radial growth of Pinus cembra trees selected for dendrometer records in 2012, 2013, and 2014.

Year	Treatment	Ip (DOY)	κ	90% DOY	A (mm)
2012	Control	167 ± 0.3a	0.044 ± 0.003a	218 ± 3.5a	3.71 ± 0.41a
	Warming	168 ± 1.9a	0.043 ± 0.003a	221 ± 5.0a	3.80 ± 0.54a
	Understory removal	166 ± 2.4a	0.041 ± 0.005a	222 ± 8.0a	4.11 ± 0.37a∗
	N treatment	168 ± 2.2a	0.043 ± 0.003a	220 ± 4.8a	4.74 ± 0.31b
2013	Control	168 ± 3.0a	0.043 ± 0.003a	221 ± 5.6a	3.09 ± 0.20a
	Warming	173 ± 3.6a	0.041 ± 0.002a	227 ± 3.9a	3.49 ± 0.60ab
	Understory removal	171 ± 2.6a	0.041 ± 0.003a	223 ± 8.4a	3.83 ± 0.34b
	N treatment	171 ± 2.6a	0.044 ± 0.002a	222 ± 3.8a	4.59 ± 0.43c
2014	Control	167 ± 1.3a	0.042 ± 0.001a	220 ± 1.7a	3.08 ± 0.53a
	Warming	167 ± 2.6a	0.044 ± 0.003a	219 ± 6.0a	3.48 ± 0.46a
	Understory removal	169 ± 2.5a	0.040 ± 0.003a	227 ± 7.3a	3.54 ± 0.45a

Table 4

Mean threshold temperatures corresponding with the 0.5-probability of active radial growth of Pinus cembra in the control, the warmed and the understory removal (≈cooling) treatment estimated during the growing seasons 2012, 2013, and 2014.

Treatment	Tair (°C)	Tsoil (°C)
Control	5.9 ± 0.6a	7.8 ± 0.4a
Warmed	5.9 ± 0.8a	9.1 ± 0.6b
Understory removal	5.1 ± 0.9a	6.0 ± 0.6c

Table 5

Total (tree and understory) and tree fine-root biomass (g dw l-1 soil); (Mean ± SD) in 0–10 cm and 10–20 cm soil depth in control (c), warmed (w), understory removal (ur) and N fertilized (N) blocks at the end of the growing season 2013.

Soil depth	Root biomass	c	w	ur	N t
0–10 cm	Tree and understory1	9.4 ± 2.9a	9.8 ± 2.0a	4.5 ± 1.6b	7.6 ± 1.4ab
0–10 cm	Tree2	1.0 ± 0.9a	1.2 ± 1.6a	0.8 ± 0.6a	0.7 ± 0.7a
10–20 cm	Tree and understory1	3.1 ± 1.4a	3.5 ± 0.6a	3.0 ± 1.1a	3.1 ± 0.7a