Lucrezia Unterholzner1*
Angela Luisa Prendin1,2
Raffaella Dibona1Roberto Menardi1
Valentino Casolo3Sara Gargiulo3,4
Francesco Boscutti3
Marco Carrer1
- 1Department of Land Environment Agriculture and Forestry, University of Padova, Legnaro, Italy
- 2Department of Biology, Ecoinformatics and Biodiversity, Aarhus University, Aarhus, Denmark
- 3Department of Agricultural Food Environmental Animal Sciences, University of Udine, Udine, Italy
- 4Department of Life Sciences, University of Trieste, Trieste, Italy
With the recent climate warming, tundra ecotones are facing a progressive acceleration of spring snowpack melting and extension of the growing season, with evident consequences to vegetation. Along with summer temperature, winter precipitation has been recently recognised as a crucial factor for tundra shrub growth and physiology. However, gaps of knowledge still exist on long-living plant responses to different snowpack duration, especially on how intra-specific and year-to-year variability together with multiple functional trait adjustments could influence the long-term responses. To fill this gap, we conducted a 3 years snow manipulation experiment above the Alpine treeline on the typical tundra species Juniperus communis, the conifer with the widest distributional range in the north emisphere. We tested shoot elongation, leaf area, stomatal density, leaf dry weight and leaf non-structural carbohydrate content of plants subjected to anticipated, natural and postponed snowpack duration. Anticipated snowpack melting enhanced new shoot elongation and increased stomatal density. However, plants under prolonged snow cover seemed to compensate for the shorter growing period, likely increasing carbon allocation to growth. In fact, these latter showed larger needles and low starch content at the beginning of the growing season. Variability between treatments slightly decreased over time, suggesting a progressive acclimation of juniper to new conditions. In the context of future warming scenarios, our results support the hypothesis of shrub biomass increase within the tundra biome. Yet, the picture is still far from being complete and further research should focus on transient and fading effects of changing conditions in the long term.
Introduction
Arctic and Alpine regions are currently experiencing higher warming rates than the global average and this trend is likely to increase in the near future, especially during winter periods (Hoegh-Guldberg et al., 2018). Temperature variations, together with decreasing frequencies of cold extremes and local alteration of the precipitation pattern (Allan et al., 2021), deeply affect the cryosphere. For example, with each decade from the 1980s, the Arctic has shown a reduction in spring (May) snow cover extent of ca. 10% and a reduction in April snow mass of 5% (Richter-Menge et al., 2016). This also induced an advance of the snow-free period by ca. 5 days per decade in the Northern Hemisphere (Dye, 2002; Bormann et al., 2018). As a result, there are already significant, cascading effects on plant communities. Species composition, distribution and relative abundance, and species traits, such as plant height, biomass allocation and canopy cover, have been altered (Myers-Smith et al., 2011), with a significant effect on plant diversity (Boscutti et al., 2018). These changes result in heterogeneous rates of Alpine and Arctic tundra greenness (Berner et al., 2020; Myers-Smith et al., 2020) and phenological responses, and in an overall shrubline advance around the globe observed in the last decades (Dial et al., 2015; Carlson et al., 2017; Myers-Smith and Hik, 2018; Steinbauer et al., 2018; Frost et al., 2020). Among the tundra systems, the European Alps at high elevation may be even more vulnerable than the Arctic, due to differences in the photoperiod and resource availability (Ernakovich et al., 2014). In parallel with northern hemisphere tendency, by the end of the twenty-first century, above the alpine treeline is predicted to see a reduction of 30–80% of snow water equivalent and an acceleration of spring snowpack melting by up to 50–60 days (Beniston et al., 2003; Steger et al., 2013). Nevertheless, the key mechanisms driving the response of plant communities to these shifting scenarios are still poorly understood, leaving uncertainties on the possible future tundra-vegetation feedback. This calls for the importance of deeply understanding the sensitivity of morphological and physiological responses of tundra species under varying climate conditions. This is especially true for shrubs species, the dominant long-lived organisms above the treeline, of which survival and fitness depend more on their capacity to acclimate/adapt to changing conditions, rather than to migrate to more fitting sites. Plants with slow regeneration times, e.g., Juniperus communis (J. communis) (García et al., 1999; Gruwez et al., 2017; Jacquemart et al., 2021) or scarce dispersal capacity might be forced to face the environmental changes, being hindered in migration (Midgley et al., 2006; Aitken et al., 2008; Sedlacek et al., 2015).
Tundra represents a typical heat-limited environment, where most of the investigations on climate change effects on vegetation focus on rising temperature, overlooking the role of changing snow cover dynamics (Dormann and Woodin, 2002; Walker et al., 2006; Elmendorf et al., 2012; Bjorkman et al., 2018). However, the date of snow melting is deemed a key factor driving the seasonal radial growth of trees (Kirdyanov et al., 2003). Recent studies demonstrated that, besides summer-temperature, winter precipitation can play a crucial role in shrub growth, opening new questions on the association between snow and shrub dynamics (Hallinger et al., 2010; Pellizzari et al., 2014; Martin et al., 2017; Carrer et al., 2019; Pandey et al., 2020). In fact, in heat limited environments, close-to-soil shrubs, with aboveground meristems below the boundary layer, may benefit from warmer microclimate conditions in respect to trees that have air-exposed meristems (Körner, 2012). Still, the influence of winter-precipitation on shrub performance is controversial due to the interaction of several effects: snow cover can provide protection against frost and soil nutrient losses, and ensure water availability but, at the same time, when abundant, can shorten the growing season, hindering growth or altering phenological dynamics (Hallinger et al., 2010; Rixen et al., 2010; Francon et al., 2017; Carrer et al., 2019; Gao et al., 2021).
So far, investigations targeting snow cover effects mainly explore single plant’s parameters as radial annual growth or shoot elongation (Rixen et al., 2010; Schmidt et al., 2010; Pellizzari et al., 2014; Francon et al., 2017; Carrer et al., 2019), while other key functional traits or physiological parameters crucial for plant fitness, productivity, and survival, as leaf traits or non-structural carbohydrate content, are mostly neglected, especially in shrub-related studies (Wipf et al., 2009; Sedlacek et al., 2015; Filippi et al., 2021). In fact, leaves and stomata features, driving photosynthesis and transpiration processes, reflect plant functioning in specific growing conditions (Chapin et al., 2002). Similarly, non-structural carbohydrates (NSCs) content, that is the pool of photosynthetic products suitable for growth, metabolism maintenance and defence, are pivotal to understanding plant response to environmental drivers (Kozlowski, 1992; Martínez-Vilalta et al., 2016). In particular, NSCs have been intensively studied at the treeline with regard to the growth limitation hypothesis (Hoch et al., 2002; Fajardo et al., 2012; Hoch and Körner, 2012; Fajardo and Piper, 2014). Recently, these investigations, on the role of NSC in response to elevational gradient, have been extended to shrubs (Dolezal et al., 2019) and dwarf shrubs (Casolo et al., 2020). However, only a few studies, with contrasting results, compared NSC content with snow depth and cover length (Palacio et al., 2015; Wheeler et al., 2016; Domisch et al., 2017).
Among the long-living shrubs, the common juniper (J. communis L.) is a cold and drought tolerant pioneer species with the widest distributional range across the northern hemisphere, representing one of the most typical tundra woody species at high latitude and elevation (Adams, 2014; Enescu et al., 2016). Hence, it represents a model species to investigate physiological responses of woody species to environmental variability above the treeline ecotone. To assess the responses of J. communis to snow cover duration, we established a snow manipulation experiment above the treeline in the Italian Alps. Primary growth, leaf traits and leaf NSC content were measured for three years to test the effects of different snowmelt timings on common juniper axial growth and carbon gain. Specifically, we expected: (i) a short-term and cumulative (across years) increasing effect of snow cover duration on shrub primary growth with the positive effect of early snowmelt and negative for late snowmelt as a response to longer/shorter growing season; (ii) a higher risk of early frost events in photosynthetic organs with earlier snowmelt, reflected by a higher concentration of soluble NSC in leaves and (iii) a significant influence on leaf traits (e.g., leaf size, dry weight and stomatal density) as the new needles tend to rapidly adjust to maximise light exposure and gas exchange for the different vegetation period.
Materials and Methods
Study Site and Experimental Design
The area is located above the treeline nearby Giau pass in the Dolomites, Italy (46°29′9′′N; 12° 3′21′′E, 2155 m a.s.l.; Figure 1A). The substrate is siliceous and soils are deep leptosol (ARPAV-Servizio Osservatorio Suolo e Bonifiche, 2011). We selected an east-facing gentle slope (ca. 20°) with morphological homogeneity and abundance of J. communis individuals growing on early snow-free mounds (Figure 1B). Within the area, we defined a plot of 30 m × 4 m, almost fully covered by the target species, and identified 15 shrubs clearly distinguished from each other. We delimited the area with a fence to avoid livestock grazing and trampling damage. We applied three different treatments, each on five randomly distributed individuals (experimental blocks) within the plot: the early snow melting, (E), from which we removed snow before the natural melting period (Figure 1C); the control (C), where snowpack remained untouched and snow melted naturally according to yearly precipitation and temperature regimes; and last, the late snow melting (L), on which we prolonged the snow cover as much as possible (Figure 1D).
Figure 1. (A) Location (46°29′9′′N; 12° 3′21′′E, 2,155 m a.s.l.) and (B) Overview of the study site. (C) Early snow melting treatment where snowpack had been removed and juniper was early air-exposed. (D) Late snow melting treatment with snowpack protected by thermal sheets. Sources: (A) Esri World Topo Map.
The day of snowpack removal on the E-plants was modulated year by year according to the meteorological conditions and snow abundance that was monitored through a webcam located within the study site and regular visits to the site. In contrast, to prolong the snow cover on the L-plants, we manually accumulated snow and applied reflective insulating thermal sheets (padding in aluminium and polyester) ca. 1.5 mm thick, to postpone the melting (Figure 1C).
The experiment was setup in early autumn 2017 and conducted continuously during the three following years 2018, 2019, and 2020.
Time Series of Climate and Microsite Condition
A meteorological station was set up within the study site to record air temperature at 2 m above the ground, through Type T Thermocouple (error ± 0.5°C). The volumetric soil water content was also measured every hour through Time Domain Reflectometry (TDR) method, with probes that were previously calibrated for the soil type (precision ± 2%; Campbell Scientific co. Mo. SC615). Snow depth records were automatically collected every 30 min from 1988 at the close Falzarego Pass (46°31′10′′N; 12° 0′33′′E, 2090 m a.s.l.) and manually every 1 or 2 weeks at Giau Pass station (46°29′9′′N; 12° 3′21′′E, 2155 m a.s.l.; ARPA Veneto—Centro Valanghe di Arabba). Since data from the two stations were significantly correlated (r = 0.99), and Falzarego Pass snow-depth data have higher resolution and longer historical record, we used the Falzarego Pass data.
Primary Growth Measurement
New shoot elongation of common juniper individuals was measured with a calliper every 7–10 days from the moment they sprouted from the buds (between May and June according to snow cover duration) until they reached the maximum length (September). Branches and shoots were randomly selected and labelled with coloured tags to be identified from one season to the next. Sample number increased over years from 1 branch × 1 shoot × individual (n = 15 shoots) in 2018 to 3 branch × 3 shoots × individual (n = 135) and finally to 6 branches × 3 shoots (n = 270) in 2020.
Leaf Traits
At the end of September 2020, we collected from each individual ca. 30 to 50 mature and healthy needles from the last 3 yearly internodes. We used 225 needles for the stomatal density analysis (5 needles × 3 years × 15 shrubs): we imprinted the lower side of the needles with nail varnish, then we peeled off the rubber with an adhesive tape, placed it on a glass slide and collected digital images of the stomatal imprints at 10X with a NIKON Eclipse 80i microscope. Then with an image analysis software (GIMP 2.10), we counted the number of stomata for a fixed 1 mm2 area.
The remaining needles were then used to define leaf area, using WinSEEDLE (V 5.0) image analysis system, and the dry weight. Dry weight was assessed on needles left 24 h in the oven at 50°C.
Finally, specific leaf area (SLA) was computed as the ratio between needle area and dry weight.
Non-structural Carbohydrates
To explore the variation of NSC content, we collected samples of leaves from 2- to 3-year-old branches of 15 plants at the beginning (June) and end (September) of each growing season. Sampling occurred between 11.00 and 12.00 a.m. on sunny days to avoid diurnal variability and plant material was immediately stored in a plastic bag inside a cool portable fridge, brought to the lab and microwaved for 3 min at 700 W to stop carbohydrate consumption. Then, we ground leaves (particles size less to 0.15 mm) with liquid nitrogen and dried them in an oven at 55°C overnight. Soluble NSC and starch content were extracted following the method proposed by Quentin et al. (2015) and Landhäusser et al. (2018). About 15 mg of ground samples were weighed into Eppendorf tubes (2 ml volume) and incubated in 0.5 ml of 80% (v/v) ethanol for 30 min at 80°C. Samples were centrifuged at 14,000 rpm for 3 min (with Mikro 120, Hettuch zentrifugen) and supernatants were transferred into new Eppendorf tubes (1.5 ml volume). A second incubation of pellets was performed with 0.5 ml ethanol. Then we added the resulting supernatant to the first and dried it in the oven until ethanol evaporation was complete. Crystallised soluble NSCs were suspended in 0.5 ml of 50 mM of Tris-HCl pH 7.5 and stored at −20°C until the application of anthrone assay (Yemm and Willis, 1954). This method allows for quantifying NSC content by means of spectrophotomertic measurements by converting absorbance at 620 nm, in the amount of glucose (g/g DW) through a titration curve obtained with known amounts of glucose.
For starch quantification, the pellet coming from sugar extraction was suspended in 1 ml of 0.2 M sodium acetate trihydrate (pH 4.6) and subjected to gelatinisation by boiling it for 1 h. Starch enzymatic hydrolysis to glucose was induced by means of 100 U of α-amylase and 25 U of amyloglucosidase per sample. The reaction was performed at 55°C in an oven overnight, then samples were boiled for five minutes to interrupt the enzymatic activity. Glucose derived from starch hydrolysis was evaluated as NADPH formation following the method proposed by Bergmeyer and Bernt (1974) by means of 0.3 U of Hexokinase and 0.5 U of Glucose-6-phosphate dehydrogenase, in buffer solution (50 mM Tris-HCl, 2 M MgCl2, 50 mM NADP+ and 0.4 M NaATP) at 32°C. NADPH produced was read at spectrophotometer at 340 nm and converted into glucose through a titration curve.
All the spectrophotometric measurements were performed with VICTOR3 Multilabel Counter Plate Reader (Perkin Elmer, Boston, MA, United States).
Statistical Analysis
Data of all variables were tested for normality and homoscedasticity, through the Shapiro–Wilk and Levene tests. To detect air temperature, soil water content and snow depth trends, we averaged data by the day of the year (doy) for 2018, 2019, and 2020 and interpolated them through a spline function.
In the case of non-normal distribution, variables were square-root transformed, except for shoot elongation for which requirements were not satisfied even after transformation. For all variables, except for new shoot elongation, differences between treatments for single years and for the average of the whole period were tested through linear mixed-effects models (nested ANOVA), accounting for the treatment as a fixed factor and experimental block as a random factor. We used the nlme R package (Pinheiro et al., 2021). Then, we run a pairwise comparison on the same models via the emmeans package in R (Lenth, 2022). Finally, to analyse shoot elongation, a generalised additive model (GAM) was used, through the mgcv R package (Wood, 2017), according to the following formula:
Where shoot elongation is modelled according to the interaction between the day of the year (doy) and treatment with year and specimen as covariates. Specimens (individual, stem and branch) are treated as nested parameters being one subgroup of the other; s represents the smoothing function that can be linear or non-linear (Wood, 2017). Also, for shoot elongation, both single years and the average for the period 2018–2020 were considered. The maximum elongation rate and the beginning time of elongation were derived from the GAM model, considering all single measurements together. All statistical analyses were performed using R software (R Core Team, 2021).
Results
Climate Data and Growing Season Length
Daily mean air temperature recorded at the experimental site followed the typical seasonal trend, peaking between mid-July and August with a value of 12.2°C and a minimum of −6.7°C in January. Daily soil water content was rather constant until mid-March and later showed two peaks, in mid-May (33.05%) and at the beginning of November (30.22%; Supplementary Figure 1A).
In the three years of the experiment, the snow melted 10–35 days earlier than usual respective to the average (1988–2021) melting date occurring at the beginning of June. Yet, snowpack depth was thicker than usual in 2018 and 2020, with 122 and 112 cm, respectively, while in 2019 it was 70 cm, equal to the average value of the whole monitoring 1988–2021 period (Supplementary Figure 1B).
Considering the average value of the three years, E-plants were air-exposed 62 days before the controls, while L plants were air-exposed 33 days after the controls (Figure 2 and Supplementary Table 1).
Figure 2. Day of the year (doy) with indication of months when common junipers for the three treatments (late-melt L, control C, and early-melt E) emerged from the snow cover in 2018, 2019, and 2020 together with the mean values for the 3 years.
Primary Growth Trend During the Growing Season
Shoot primary growth followed a sigmoidal trend for all treatments and all years (Figure 3 and Supplementary Figure 2). At the end of the growing season, the final length was 22.80 ± 3.23 mm (mean ± SE) for E-shrubs, 21.02 ± 2.33 mm for C-shrubs and 20.55 ± 2.07 for L-shrubs (Figure 3 and Supplementary Table 2).
Figure 3. GAM modelled elongation trend (mm) with 95% confidence intervals of Juniperus communis new shoots over the growing season for the control (green), early-melt (orange) and late-melt (blue) treatment.
All the predictive variables considered in the model had a significant effect on shoot elongation (Table 1).
Table 1. Outcome of GAM model applied to elongation as an independent variable in relation to year, shrub, stem, and snow treatment.
Maximum growth rate, computed as the 1st derivate of GAM-modelled elongation curves, was 0.55, 0.54, and 0.50, and it occurred on day 183, 182, and 196, for E-, C- and L-plants, respectively (Supplementary Table 3).
The beginning of shoot elongation did not reflect the number of days of snow removal between treatments. In fact, considering the three years together, despite the marked difference in the dates of snow disappearing between treatments (Figure 2 and Supplementary Table 1), E- and C-shoots started to grow in June with just a mean difference across the 3 years of 1 day (doy 170 for the E and 171 for the C) and C and L with a difference of 5 days (doy 176 for the L).
Comparing the elongation dynamic between the three years, while the control remained stable both in terms of maximum elongation rate and final length, the E- and L-plants did not feature any clear trend. However, for the manipulated groups, we observed a progressive reduction of the inter-individual and inter-group variability from the 1st to the 3rd year (Supplementary Figure 2 and Table 2).
Table 2. F-value, p-values of the linear mixed-effects models (nested ANOVA), and degrees of freedom (DF) of stomatal density, leaf area, leaf dry weight, leaf soluble sugars and leaf starch content in 2018, 2019, 2020, and averaging these 3 years (mean).
Treatment Effects on Leaf Traits
Stomatal density was 266 ± 4, 260 ± 4, and 254 ± 4 n°/mm2 for the E-, C- and L-groups, respectively. These values were not significantly different considering the single years and remained almost constant across the years for all the treatments. However, considering the whole period (2018–2020) stomatal density was significantly higher in E-plants compared to L-plants (Figure 4A, Table 2, and Supplementary Table 4).
Figure 4. Outcome of linear mixed-effects models (nested ANOVA) of leaf traits. Mean ± SE stomatal density (A; n°/mm2), needle area (B; mm2) and dry weight per leaf (C; mg) of early-melt (orange), control (green), and late-melt (blue) shrubs in average from 2018 to 2020. Different letters indicate significant differences between groups (p < 0.05).
Leaf area was on average 9.19 ± 0.09, 8.80 ± 0.09, and 9.78 ± 0.09 mm2 for the E-, C-, and L groups (Figure 4B and Supplementary Table 4) respectively, with significant difference between treatments considering both single years and all the three years together (Table 2). L-group had a larger leaf area than C for all time periods and a larger leaf area than E-plants in 2019 and the period 2018–2020. In addition, the leaf area did not present any trend over time but had a tendency to reduce the differences between treatments across the years and showed a step-like increase for all treatments from 2019.
Needle dry weight resulted 2.15 ± 0.11 mg for E-groups, 1.97 ± 0.08 mg for C and 2.25 ± 0.11 mg for the L (Figure 4C and Supplementary Table 4). This trait highlights no significant differences between treatments and no trend over time (Table 2).
Regarding SLA, no significant differences emerged between treatments (p-value = 0.82; Supplementary Figure 4).
Treatment Effects on Non-structural Carbohydrates
At the beginning of the growing season (June), soluble NSC content was on average 0.24 ± 0.02 g/g for E, 0.27 ± 0.01 g/g for C and 0.28 ± 0.01 g/g for the L-plants, while at the end of the vegetative period (September), it was 0.26 ± 0.01 g/g, 0.25 ± 0.03 g/g and 0.30 ± 0.02 g/g for the E, C and L, respectively (Figure 5B and Supplementary Table 6). The models revealed no significant differences between treatments neither at the beginning nor at the end of the growing season when considering the period 2018–2020 (Table 2). However, in June sugar content was slightly lower in E-shrubs in contrast to other treatments, even if not significantly (Figure 5B). In 2020, E-shrubs showed significantly lower soluble NSC content than the other two groups at the beginning of the season. Moreover, at the end of the growing season, C-plants and E-plants showed lower sugar content in 2019 and 2020, respectively (Table 2 and Supplementary Figure 5A).
Figure 5. Outcome of linear mixed-effects models (nested ANOVA) of NSC leaf content. Mean ± SE starch (g/g; A) and NSC content (g/g; B) of early-melt (E; orange), control (C; green), and late-melt (L; blue) shrubs, considering the mean from 2018 to 2020, measured in June and September. Different letters indicate significant differences between groups (p < 0.05).
On the contrary, starch content showed significant differences between treatments at the beginning of the season, but not at the end, for both single years and the whole period. The only exception is at the end of the 2020 season when L-plants have higher starch content than the other two groups. At the beginning of the season, L-plants had significantly less starch (Table 2 and Supplementary Figure 5A). Starch content at the beginning of the vegetative season was 0.08 ± 0.005 g/g for E, 0.08 ± 0.003 g/g for C and 0.02 ± 0.03 g/g for L, while at the end it was 0.01 ± 0.002 g/g for E, 0.01 ± 0.01 g/g for C and 0.02 ± 0.004 g/g for L (Figure 5A and Supplementary Table 5).
Non-structural carbohydrates between treatments did not show a clear trend over time or a specific response to treatments over years. However, variability of starch content increased within groups and decreased between treatments (Supplementary Figure 5 and Supplementary Tables 5, 6).
Discussion
With a manipulative 3-year experiment, we monitored the immediate and carryover effects of different snowpack duration on primary growth, leaf traits and NSC leaf content on J. communis growing above the treeline.
According to our initial hypothesis, a few traits were affected by snow cover duration already from the first year, namely primary growth, leaf area and leaf starch content. However, in contrast to our expectations, the differences between treatments for these traits seemed to decrease progressively, except for soluble NSC content. In fact, soluble NSC, together with stomatal density showed an overall response to different snow cover treatments only at the end of the experiment or looking at the whole period. Instead, other traits (growth rate, leaf dry weight and SLA) seemed rather unaffected by the varying persistence of the snow cover.
Above the alpine treeline, early frost is generally considered a crucial stressor, especially on early-exposure sites (Wheeler et al., 2014). This suggests that earlier exposure of plants to air compared to natural snowpack melting might increase the risk of damages due to freezing events, yet with the concurrent effect of accelerating the onset of the growing season (Björk and Molau, 2007). At this elevation, freezing temperatures can occur during the vegetative season, highlighting the lack of a frost-free period (Larcher and Wagner, 2009; Neuner, 2014). Common juniper living in harsh conditions is well-adapted to frost (Thomas et al., 2007), and freezing risk might not be considered as a crucial limiting factor for E-plants in this experiment. Additionally, juniper needles withstood temperatures of up to −10°C (Taschler and Neuner, 2004) also during the summer, when frost resistance is usually lower compared to winter or early spring (Neuner, 2014). Accordingly, we did not observe any significant sign of frost damage, like dead branches or needles, during the three years. The absence of frost damages on the morphology of the needles, besides the shrub resistance to frost, could also be related to the relatively low number of random events experienced by plants during the experiment, such as intense frost (Supplementary Figure 1A).
The low frost impact on E-plants was also confirmed by the fact that these started to elongate earlier, at a higher rate throughout the summer to finally reach longer extension, 8.5% more than control shoots, at the end of the season (Figure 3 and Supplementary Table 2). This may suggest that early snow-cover removal has a potential direct effect on the length of the vegetation period with E-plants benefiting from earlier exposure to daylight. Since stomata are the regulator of photosynthetic performance (Hetherington and Woodward, 2003) and higher stomatal density might improve plant productivity (Schlueter et al., 2003), we may also speculate that the observed enhanced elongation could result from the higher stomatal density detected on E-individuals. In fact, the stomatal density of developing leaves, besides CO2 concentration and plant hormones, is positively related to temperature (Schlueter et al., 2003) and irradiance received by mature leaves (P. W. Hetherington and Woodward, 2003; Thomas et al., 2004; Casson and Gray, 2008). E- and L-plants showed a slight but divergent response in terms of stomatal density when considering the whole period (Figure 4A). Accordingly, considering that light is one of the key varying factors between treatments, we may hypothesise that E-plants received overall more light than L-plants, letting more epidermal cells develop into stomata (Schoch et al., 1980), and vice versa.
Late melting-shoots, despite the delayed and lower peak of elongation rate, reached maximal length just 2.2% shorter than controls at the end of the vegetative period. The reduced growth is in line with the detrimental effect of snow cover persistence on long-term radial growth, that similarly to primary growth is linked to carbon gain, of common juniper across the Alps (Carrer et al., 2019). However, J. communis showed high plasticity in growth phenology that could potentially allow it to tune its growth rates during the growing seasons to better cope with the peculiar pace of climate variability, and modulate carbon investment accordingly (Tumajer et al., 2021). Plasticity in the growth rate could be the key attribute permitting L-plants to recover the growth lag in respect to control and E-plants. This growth acclimation is shared with other shrub species able to colonize a wide range of environments, such as Salix artica and Rhododendron ferrugineum growing in heat-limited environments, of which radial growth is negatively affected by snow persistence (Schmidt et al., 2010; Francon et al., 2020). On the contrary, early snowmelt did not boost the primary growth of other alpine species, such as Salix herbacea, Vaccinium myrtillus, Vaccinium uliginosum, and Loiseleuria procumbens (Wheeler et al., 2014). These differences could likely derive from different species and growing conditions, besides different experimental methodologies that can highlight or attenuate such a result. In general, the responses in terms of primary growth of shrubs to changing snow persistence seem rather heterogeneous (Wipf and Rixen, 2010), highlighting the need for species-specific studies and stressing the role of microclimatic conditions in promoting growth (Myers-Smith et al., 2015; Tumajer et al., 2021).
Nevertheless, the upswing of L-shoot growth may also be supported by larger leaves, that is larger photosynthetic area found in L-plants, as they increase the radiation absorption and therefore carbon input (Chapin et al., 2002). In fact, L-shrubs produced significantly larger leaves than both C- and E-shrubs (Figure 4B and Table 2), featuring a plastic response in this leaf trait in relation to altered snowmelt timings, in accordance with our hypothesis. The reason behind larger leaves could be connected to the mean warmer conditions experienced by L-shrubs that should favour larger leaf area. Despite the shorter growing season, L-shrubs emerged from the snowpack in correspondence with the peak mean air temperature (Figure 2 and Supplementary Figure 1A). Similar investigations on S. herbacea in the Alps (Sedlacek et al., 2015) and, more in general, across the tundra biome found a similar outcome, with a positive association between temperature and leaf area (Bjorkman et al., 2018). Likely, together with the larger photosynthetic area, also higher efficiency of light may have enhanced L-plants growth, in line with the response of Arctic tundra vegetation to shorter growing season caused by prolonged snow cover (Bosiö et al., 2014). In parallel, the leaf dry weight of L-plants was just slightly higher but not significantly affected by snow removal (Figure 4C and Table 2). This result, together with the larger L-needles, seems in accordance with the tendency of leaves to reduce their construction costs (represented by dry weight) and increase the photosynthetic rate with increasing latitude to maximally exploit the short growing season (Tian et al., 2016).
In general, NSC content in C-plants is consistent with seasonal data proposed by Rabska et al. (2020), even though, surprisingly, another study on the same species proposed a soluble NSC content 100-fold lower (Magaña Ugarte et al., 2021). The high level of soluble NSC we measured can be attributed to the time of collection that corresponds to the maximal diurnal photosynthetic rate (Bickford et al., 2009). Nonetheless, we consider that the observed differences in NSC should be attributed to the plant response to the distinct snowpack duration. In fact, the evergreen common juniper does not rely on new needles to perform photosynthesis when it is out from snow cover and its photosynthetic capacity does not dramatically vary seasonally (Rabska et al., 2021).
Plants exposed to abiotic stresses often accumulate starch that can be remobilized for energy and carbon supply, under conditions of photosynthesis limitation (Thalmann and Santelia, 2017). However, starch in leaves is a transitory form accumulated in the chloroplast during active photosynthesis. Focusing on the lower level of starch, observed in leaves of L-plants in late spring (Figure 5), we may notice that these plants likely developed a good level of photosynthetic rate; in fact, no significant differences were observed in the sugar pool. This suggests that transitory starch was poorly accumulated in the chloroplast because sucrose is promptly exported to the sink organ actively growing, e.g., buds (Ray and Savage, 2021). Consistently, L-plants showed a high growth rate in late spring, i.e., when NSCs were measured (see Figure 3 doy 200–210).
Moreover, in our case, the primary growth dynamic seemed significantly influenced not just by the treatment but also by other factors related to time (year) and the individual (Table 1). For example, selected stem and branch, besides shrub, highlighted a significant effect on primary growth (Table 1), suggesting that the intra-individual variability could be potentially as high, or even higher, than the intraspecific variability. This may reflect the intraspecific homogeneity, at least for some physiological traits, such as functional traits related to hydraulic efficiency, of this taxon (Unterholzner et al., 2020).
We found a remarkable converging tendency between treatments over years for primary growth dynamics, leaf area, leaf dry weight and, partially, also for starch leaf content. This may suggest a plant’s progressive acclimation to the new conditions across time. On one hand, L-plants seem to enhance their efficiency in exploiting the shorter growing seasons, on the other hand, E-plants seem to lose their advantage of experiencing longer growing seasons. In the long term, E-plants may be affected by the harsher soil conditions, in terms of microbial communities, when exposed to early snowmelt (Broadbent et al., 2021), and in parallel soil under prolonged snow cover may increase N soil availability (Semenchuk et al., 2015).
These results are even more remarkable when considering that, despite the difference of snow-removal days between treatments (Figure 3), shrubs started to grow within just a few days of each other (E/L shoots 1/5 days before/after the controls), when daily mean air temperature was around 10°C (Figure 2A). Processes needed for alpine plants to grow to require a certain temperature threshold, below which they do not activate (Körner, 2003). Since snowpack duration, in parallel with temperature increase, is predicted to drastically shrink in the next decades in both Alpine (Steger et al., 2013) and Arctic regions (Meredith et al., 2019), a substantial phytomass increase may be expected. This is consistent with the overall, even if heterogeneous, current tundra shrub encroachment and with the ca. 20% biomass increase detected in the Northern hemisphere in the recent decades (Epstein et al., 2012; Frost et al., 2020). Nonetheless, future predictions are far from being straightforward considering the high plasticity of this species and the potential parallel change of other external factors as a consequence of snow regime change. Finally, since changes in growing season length may also influence the interplay between plants and pathogens, pollinators or herbivores, we could expect potential reshaping in trophic interactions (Roy et al., 2004; Høye et al., 2013; Sedlacek et al., 2015). Also, even if individual performance may benefit from the future climate, predicting potential vegetation shifts would benefit from also considering climatic responses of shrub community. For instance, the benefit of warmer temperature for recruitment of Juniper species at high elevation may be hindered by water stress, highlighting the role of abiotic interaction on vegetation dynamics (Sigdel et al., 2021). Indeed, common juniper (but also other plants) growing on snow beds could benefit from the additional soil water available due to, e.g., late snowmelt. Still, additional studies are needed to disentangle the role of the different abiotic (environmental and climatic) and biotic (e.g., competition) components to get insight into shrubs’ responses. Therefore, we encourage future studies to focus not only on the individual but also on the community and ecosystem levels.
Data Availability Statement
The raw data supporting the conclusions of this article will be made available by the authors on request.
Author Contributions
MC conceived the experiment and designed the study with input from the co-authors. LU, ALP, RD, RM, and MC contributed to the experiment maintenance, collected samples, and performed fieldwork activities. LU measured stomatal density and leaf area, analysed the data with input from ALP and FB, and designed tables and figures. VC, SG, and FB contributed to run non-structuralcarbohydrates measurements. LU interpreted the results with inputs from MC, ALP, VC, SG, and FB and drafted the 1st version of the manuscript with inputs from MC and the contribution of all other co-authors. All authors contributed to the article and read and approved the submitted version.
Funding
ALP was supported by the 2017 BIRD Project of TeSAF Department University of Padova and partially by a Marie Sklodowska-Curie Individual Fellowship (IF) under contract number 895233.
Conflict of Interest
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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Vinicio Carraro, Tommaso Baggio, Filippo Biscontin, Elisabetta Candeago, Nicola Cappellari, Mattia Gioele Prendin, and Matteo Zava for their support in the experiment maintenance and data collection over the 3 years. We also thank Ashley Pearcy Buitenwerf for proofreading the manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.822901/full#supplementary-material
References
Adams, R. P. (2014). Junipers of the world: the genus Juniperus, 4th Edn. Bloomington, IN: Trafford Publishing Co.
Aitken, S. N., Yeaman, S., Holliday, J. A., Wang, T., and Curtis-McLane, S. (2008). Adaptation, migration or extirpation: climate change outcomes for tree populations. Evol. Appl. 1, 95–111. doi: 10.1111/j.1752-4571.2007.00013.x
Allan, R. P., Hawkins, E., Bellouin, N., and Collins, B. (2021). IPCC, 2021: summary for Policymakers. Cambridge: Cambridge University Press.
ARPAV-Servizio Osservatorio Suolo e Bonifiche (2011). Regione del Veneto – L.R. n. 28/76 – Formazione della Carta Tecnica Regionale. Available online at: https://idt2.regione.veneto.it/* (accessed March 8, 2022).
Beniston, M., Keller, F., Koffi, B., and Goyette, S. (2003). Estimates of snow accumulation and volume in the Swiss Alps under changing climatic conditions. Theor. Appl. Climatol. 76, 125–140. doi: 10.1007/s00704-003-0016-5
Bergmeyer, H. U., and Bernt, E. (1974). Methods of Enzymatic Analysis, 2nd Edn. Cambridge, Massachusetts: Academic Press, Inc. doi: 10.1016/B978-0-12-091302-2.X5001-4
Berner, L. T., Massey, R., Jantz, P., Forbes, B. C., Macias-Fauria, M., Myers-Smith, I., et al. (2020). Summer warming explains widespread but not uniform greening in the Arctic tundra biome. Nat. Commun. 11, 1–12. doi: 10.1038/s41467-020-18479-5
Bickford, C. P., McDowell, N. G., Erhardt, E. B., and Hanson, D. T. (2009). High-frequency field measurements of diurnal carbon isotope discrimination and internal conductance in a semi-arid species, Juniperus monosperma. Plant Cell Environ. 32, 796–810. doi: 10.1111/j.1365-3040.2009.01959.x
Björk, R. G., and Molau, U. (2007). Ecology of alpine snowbeds and the impact of global change. Arct. Antarct. Alp. Res. 39, 34–43. doi: 10.1657/1523-0430(2007)39[34:eoasat]2.0.co;2
Bjorkman, A. D., Myers-Smith, I. H., Elmendorf, S. C., Normand, S., Rüger, N., Beck, P. S. A. A., et al. (2018). Plant functional trait change across a warming tundra biome. Nature 562, 57–62. doi: 10.1038/s41586-018-0563-7
Bormann, K. J., Brown, R. D., Derksen, C., and Painter, T. H. (2018). Estimating snow-cover trends from space. Nat. Clim. Chang. 8, 924–928. doi: 10.1038/s41558-018-0318-3
Boscutti, F., Casolo, V., Beraldo, P., Braidot, E., Zancani, M., and Rixen, C. (2018). Shrub growth and plant diversity along an elevation gradient : evidence of indirect effects of climate on alpine ecosystems. PLoS One 13:e0196653. doi: 10.1371/journal.pone.0196653
Bosiö, J., Stiegler, C., Johansson, M., Mbufong, H. N., and Christensen, T. R. (2014). Increased photosynthesis compensates for shorter growing season in subarctic tundra—8 years of snow accumulation manipulations. Clim. Change 127, 321–334. doi: 10.1007/s10584-014-1247-4
Broadbent, A. A. D., Bahn, M., Pritchard, W. J., Newbold, L. K., Goodall, T., Guinta, A., et al. (2021). Shrub expansion modulates belowground impacts of changing snow conditions in alpine grasslands. Ecol. Lett. 25, 52–64. doi: 10.1111/ele.13903
Carlson, B. Z., Corona, M. C., Dentant, C., Bonet, R., Thuiller, W., and Choler, P. (2017). Observed long-term greening of alpine vegetation - A case study in the French Alps. Environ. Res. Lett. 12:114006. doi: 10.1088/1748-9326/aa84bd
Carrer, M., Pellizzari, E., Prendin, A. L., Pividori, M., and Brunetti, M. (2019). Winter precipitation - not summer temperature - is still the main driver for Alpine shrub growth. Sci. Total Environ. 682, 171–179. doi: 10.1016/j.scitotenv.2019.05.152
Casolo, V., Braidot, E., Petrussa, E., Zancani, M., Vianello, A., and Boscutti, F. (2020). Relationships between population traits, nonstructural carbohydrates, and elevation in alpine stands of Vaccinium myrtillus. Am. J. Bot. 107, 639–649. doi: 10.1002/ajb2.1458
Casson, S., and Gray, J. E. (2008). Influence of environmental factors on stomatal development. New Phytol. 178, 9–23. doi: 10.1111/J.1469-8137.2007.02351.X
Chapin, S. F. III, Maton, P. A., and Vitousek, P. M. (2002). Principles of Terrestrial Ecosystem Ecology, 2nd Edn. Berlin: Springer.
Dial, R. J., Scott Smeltz, T., Sullivan, P. F., Rinas, C. L., Timm, K., Geck, J. E., et al. (2015). Shrubline but not treeline advance matches climate velocity in montane ecosystems of south-central Alaska. Glob. Chang. Biol. 22, 1841–1856. doi: 10.1111/GCB.13207
Dolezal, J., Kopecky, M., Dvorsky, M., Macek, M., Rehakova, K., Capkova, K., et al. (2019). Sink limitation of plant growth determines tree line in the arid Himalayas. Funct. Ecol. 33, 553–565. doi: 10.1111/1365-2435.13284
Domisch, T., Martz, F., Repo, T., and Rautio, P. (2017). Winter survival of Scots pine seedlings under di fferent snow conditions. Tree Physiol. 38, 602–616. doi: 10.1093/treephys/tpx111
Dormann, C. F., and Woodin, S. J. (2002). Climate change in the Arctic: using plant functional types in a meta-analysis of field experiments. Funct. Ecol. 16, 4–17. doi: 10.1046/j.0269-8463.2001.00596.x
Dye, D. G. (2002). Variability and trends in the annual snow-cover cycle in Northern Hemisphere land areas, 1972-2000. Hydrol. Process. 16, 3065–3077. doi: 10.1002/hyp.1089
Elmendorf, S. C., Henry, G. H. R., Hollister, R. D., Björk, R. G., Boulanger-Lapointe, N., Cooper, E. J., et al. (2012). Plot-scale evidence of tundra vegetation change and links to recent summer warming. Nat. Clim. Chang. 2, 453–457. doi: 10.1038/nclimate1465
Enescu, C. M., Houston Durrant, T., Cadullo, G., and de Rigo, D. (2016). “Juniperus communis in Europe: distribution, habitat, usage and threats,” in European Atlas of Forest Tree Species, eds J. San-Miguel-Ayanz, D. de Rigo, G. Caudullo, T. Houston Durrant, and A. Mauri (Luxembourg: Publ. Off. EU).
Epstein, H. E., Raynolds, M. K., Walker, D. A., Bhatt, U. S., Tucker, C. J., and Pinzon, J. E. (2012). Dynamics of aboveground phytomass of the circumpolar Arctic tundra during the past three decades. Environ. Res. Lett. 7:015506. doi: 10.1088/1748-9326/7/1/015506
Ernakovich, J. G., Hopping, K. A., Berdanier, A. B., Simpson, R. T., Kachergis, E. J., Steltzer, H., et al. (2014). Predicted responses of arctic and alpine ecosystems to altered seasonality under climate change. Glob. Chang. Biol. 20, 3256–3269. doi: 10.1111/gcb.12568
Fajardo, A., and Piper, F. I. (2014). An experimental approach to explain the southern Andes elevational treeline. Am. J. Bot. 101, 788–795. doi: 10.3732/ajb.1400166
Fajardo, A., Piper, F. I., Pfund, L., Körner, C., and Hoch, G. (2012). Variation of mobile carbon reserves in trees at the alpine treeline ecotone is under environmental control. New Phytol. 195, 794–802. doi: 10.1111/j.1469-8137.2012.04214.x
Filippi, A., Braidot, E., Petrussa, E., Fabro, M., Vuerich, M., and Boscutti, F. (2021). Plant growth shapes the effects of elevation on the content and variability of flavonoids in subalpine bilberry stands. Plant Cell Environ. 23, 241–249. doi: 10.1111/plb.13194
Francon, L., Corona, C., Roussel, E., Lopez-Saez, J., Stoffel, M., Lopez Saez, J., et al. (2017). Warm summers and moderate winter precipitation boost Rhododendron ferrugineum L. growth in the Taillefer massif (French Alps). Sci. Total Environ. 586, 1020–1031. doi: 10.1016/j.scitotenv.2017.02.083
Francon, L., Corona, C., Till-Bottraud, I., Choler, P., Carlson, B. Z., Charrier, G., et al. (2020). Assessing the effects of earlier snow melt-out on alpine shrub growth: the sooner the better? Ecol. Indic. 115:106455. doi: 10.1016/j.ecolind.2020.106455
Frost, G. V., Bhatt, U. S., Epstein, H. E., Myers-Smith, I., Phoenix, G. K., Berner, L. T., et al. (2020). “Tundra Greenness,” in Arctic Report Card 2020, eds J. Richter-Menge, M. L. Druckenmiller, and M. Jeffries (Washington, DC: NOAA), doi: 10.25923/46rm-0w23
Gao, D., Bai, E., Yang, Y., Zong, S., and Hagedorn, F. (2021). A global meta-analysis on freeze-thaw effects on soil carbon and phosphorus cycling. Soil Biol. Biochem. 159:108283. doi: 10.1016/j.soilbio.2021.108283
García, D., Zamora, R., Hódar, J. A., and Gómez, J. M. (1999). Age structure of Juniperus communis L. in the Iberian peninsula: conservation of remnant populations in Mediterranean mountains. Biol. Conserv. 87, 215–220. doi: 10.1016/S0006-3207(98)00059-7
Gruwez, R., De Frenne, P., De Schrijver, A., Vangansbeke, P., and Verheyen, K. (2017). Climate warming and atmospheric deposition affect seed viability of common juniper (Juniperus communis) via their impact on the nutrient status of the plant. Ecol. Res. 32, 135–144. doi: 10.1007/s11284-016-1422-3
Hallinger, M., Manthey, M., and Wilmking, M. (2010). Establishing a missing link: warm summers and winter snow cover promote shrub expansion into alpine tundra in Scandinavia. New Phytol. 186, 890–899. doi: 10.1111/j.1469-8137.2010.03223.x
Hetherington, A. M., and Woodward, F. I. (2003). The role of stomata in sensing and driving environmental change. Nature 424, 901–908. doi: 10.1038/nature01843
Hoch, G., and Körner, C. (2012). Global patterns of mobile carbon stores in trees at the high-elevation tree line. Glob. Ecol. Biogeogr. 21, 861–871. doi: 10.1111/j.1466-8238.2011.00731.x
Hoch, G., Popp, M., and Körner, C. (2002). Altitudinal increase of mobile carbon pools in Pinus cembra suggests sink limitation of growth at the Swiss treeline. Oikos 98, 361–374. doi: 10.1034/j.1600-0706.2002.980301.x
Hoegh-Guldberg, O., Jacob, D., Taylor, M., Bindi, M., Brown, S., Camilloni, I., et al. (2018). “Impacts of 1.5°C global warming on natural and human systems,” in Global Warming of 1.5°C. An IPCC Special Report on the impacts of global warming of 1.5°C above pre-industrial levels and related global greenhouse gas emission pathways, in the context of strengthening the global response to the threat of climate change, sustainable development, and efforts to eradicate poverty, eds V. P. Masson-Delmotte, H.-O. Zhai, D. Pörtner, J. Roberts, P. R. Skea, A. Shukla, et al. (Rome, Italy: IPPC).
Høye, T. T., Post, E., Schmidt, N. M., Trøjelsgaard, K., and Forchhammer, M. C. (2013). Shorter flowering seasons and declining abundance of flower visitors in a warmer Arctic. Nat. Clim. Chang. 3, 759–763. doi: 10.1038/nclimate1909
Jacquemart, A. L., Buyens, C., Delescaille, L. M., and Van Rossum, F. (2021). Using genetic evaluation to guide conservation of remnant Juniperus communis (Cupressaceae) populations. Plant Biol. 23, 193–204. doi: 10.1111/plb.13188
Kirdyanov, A., Hughes, M., Vaganov, E., Schweingruber, F., and Silkin, P. (2003). The importance of early summer temperature and date of snow melt for tree growth in the Siberian Subarctic. Trees Struct. Funct. 17, 61–69. doi: 10.1007/s00468-002-0209-z
Körner, C. (2003). Alpine Plant Life: functional Plant Ecology of High Mountain Ecosystems, 2nd Edn. Berlin: Springer.
Körner, C. (2012). Treelines Will be Understood Once the Functional Difference Between a Tree and a Shrub Is. Ambio 41, 197–206. doi: 10.1007/s13280-012-0313-2
Kozlowski, T. T. (1992). Carbohydrate sources and sinks in woody plants. Bot. Rev. 58, 107–222. doi: 10.1007/BF02858600
Landhäusser, S. M., Chow, P. S., Turin Dickman, L., Furze, M. E., Kuhlman, I., Schmid, S., et al. (2018). Standardized protocols and procedures can precisely and accurately quantify non-structural carbohydrates. Tree Physiol. 38, 1764–1778. doi: 10.1093/treephys/tpy118
Larcher, W., and Wagner, J. (2009). High mountain bioclimate: temperatures near the ground recorded from the timberline to the nival zone in the Central Alps. Contrib. Nat. Hist. 12, 857–874.
Lenth, R. V. (2022). emmeans: estimated Marginal Means, aka Least-Squares Means. Available online at: https://CRAN.R-project.org/package=emmeans
Magaña Ugarte, R., Escudero, A., and Gavilán, R. G. (2021). Assessing the Role of Selected Osmolytes in Mediterranean High-Mountain Specialists. Front. Ecol. Evol. 9:576122. doi: 10.3389/fevo.2021.576122
Martin, A. C., Jeffers, E. S., Petrokofsky, G., Myers-Smith, I., and MacIas-Fauria, M. (2017). Shrub growth and expansion in the Arctic tundra: an assessment of controlling factors using an evidence-based approach. Environ. Res. Lett. 12:085007. doi: 10.1088/1748-9326/aa7989
Martínez-Vilalta, J., Sala, A., Asensio, D., Galiano, L., Hoch, G., Palacio, S., et al. (2016). Dynamics of non-structural carbohydrates in terrestrial plants: a global synthesis. Ecol. Monogr. 86, 495–516. doi: 10.1002/ecm.1231
Meredith, M., Sommerkorn, M., Cassotta, S., Derksen, C., Ekaykin, A., Hollowed, A., et al. (2019). ““Polar Regions,” in IPCC Special Report on the Ocean and Cryosphere in a Changing Climate, eds H.-O. Pörtner, D. C. Roberts, V. Masson-Delmotte, P. Zhai, M. Tignor, E. Poloczanska, et al. (Italy: IPPC), 203–320.
Midgley, G. F., Hughes, G. O., Thuiller, W., and Rebelo, A. G. (2006). Migration rate limitations on climate change-induced range shifts in Cape Proteaceae. Divers. Distrib. 12, 555–562. doi: 10.1111/j.1366-9516.2006.00273.x
Myers-Smith, I. H., Elmendorf, S. C., Beck, P. S. A., Wilmking, M., Hallinger, M., Blok, D., et al. (2015). Climate sensitivity of shrub growth across the tundra biome. Nat. Clim. Chang. 5, 887–891. doi: 10.1038/nclimate2697
Myers-Smith, I. H., Forbes, B. C., Wilmking, M., Hallinger, M., Lantz, T., Blok, D., et al. (2011). Shrub expansion in tundra ecosystems: dynamics, impacts and research priorities. Environ. Res. Lett. 6, 657–660. doi: 10.1088/1748-9326/6/4/045509
Myers-Smith, I. H., and Hik, D. S. (2018). Climate warming as a driver of tundra shrubline advance. J. Ecol. 106, 547–560. doi: 10.1111/1365-2745.12817
Myers-Smith, I. H., Kerby, J. T., Phoenix, G. K., Bjerke, J. W., Epstein, H. E., Assmann, J. J., et al. (2020). Complexity revealed in the greening of the Arctic. Nat. Clim. Chang. 10, 106–117. doi: 10.1038/s41558-019-0688-1
Neuner, G. (2014). Frost resistance in alpine woody plants. Front. Plant Sci. 5:654. doi: 10.3389/fpls.2014.00654
Palacio, S., Lenz, A., Wipf, S., Hoch, G., and Rixen, C. (2015). Bud freezing resistance in alpine shrubs across snow depth gradients. Environ. Exp. Bot. 118, 95–101. doi: 10.1016/j.envexpbot.2015.06.007
Pandey, J., Sigdel, S. R., Lu, X., Salerno, F., Dawadi, B., Liang, E., et al. (2020). Early growing-season precipitation drives radial growth of alpine juniper shrubs in the central Himalayas. Geogr. Ann. Ser. A Phys. Geogr. 102, 317–330. doi: 10.1080/04353676.2020.1761097/FORMAT/EPUB
Pellizzari, E., Pividori, M., and Carrer, M. (2014). Winter precipitation effect in a mid-latitude temperature-limited environment: the case of common juniper at high elevation in the Alps. Environ. Res. Lett. 9:104021. doi: 10.1088/1748-9326/9/10/104021
Pinheiro, J., Bates, D., DebRoy, S., Sarkar, D., and Team, R. C. (2021). nlme: linear and Nonlinear Mixed Effects Models. Available Online at: https://cran.r-project.org/package=nlme* (accessed March 8, 2022).
Quentin, A. G., Pinkard, E. A., Ryan, M. G., Tissue, D. T., Baggett, L. S., Adams, H. D., et al. (2015). Non-structural carbohydrates in woody plants compared among laboratories. Tree Physiol. 35, 1146–1165. doi: 10.1093/treephys/tpv073
R Core Team (2021). R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
Rabska, M., Pers-Kamczyc, E., Żytkowiak, R., Adamczyk, D., and Iszkuło, G. (2020). Sexual dimorphism in the chemical composition of male and female in the dioecious tree, Juniperus communis L., growing under different nutritional conditions. Int. J. Mol. Sci. 21, 1–16. doi: 10.3390/ijms21218094
Rabska, M., Robakowski, P., Ratajczak, E., Żytkowiak, R., Iszkuło, G., and Pers-Kamczyc, E. (2021). Photochemistry differs between male and female Juniperus communis L. independently of nutritional availability. Trees Struct. Funct. 35, 27–42. doi: 10.1007/s00468-020-02004-8
Ray, D. M., and Savage, J. A. (2021). Seasonal changes in temperate woody plant phloem anatomy and physiology: implications for long-distance transport. AoB Plants 13, 1–10. doi: 10.1093/aobpla/plab028
Richter-Menge, J., Overland, J. E., and Mathis, J. T. (2016). Arctic Report Card 2016: persistent warming trend and loss of sea ice are triggering extensive Arctic changes. Available Online at: http://www.arctic.noaa.gov/Report-Card* (accessed March 8, 2022).
Rixen, C., Schwoerer, C., and Wipf, S. (2010). Winter climate change at different temporal scales in Vaccinium myrtillus, an Arctic and alpine dwarf shrub. Polar Res. 29, 85–94. doi: 10.1111/j.1751-8369.2010.00155.x
Roy, B. A., Güsewell, S., and Harte, J. (2004). Response Of Plant Pathogens And Herbivores To A Warming Experiment. Ecology 85, 2570–2581. doi: 10.1890/03-0182
Schlueter, U., Muschak, M., Berger, D., and Altmann, T. (2003). Photosynthetic performance of an Arabidopsis mutant with elevated stomatal density (sdd1-1) under different light regimes. J. Exp. Bot. 54, 867–874. doi: 10.1093/jxb/erg087
Schmidt, N. M., Baittinger, C., Kollmann, J., and Forchhammer, M. C. (2010). Consistent Dendrochronological Response of the Dioecious Salix arctica to Variation in Local Snow Precipitation across Gender and Vegetation Types. Arct. Antarct. Alp. Res. 42, 471–475. doi: 10.1657/1938-4246-42.4.471
Schoch, P. G., Zinsou, C., and Sibi, M. (1980). Dependence of the Stomatal Index on Environmental Factors during Stomatal Differentiation in Leaves of Vigna sinensis. J. Exp. Bot. 31, 1211–1216. doi: 10.1093/JXB/31.5.1211
Sedlacek, J., Wheeler, J. A., Cortés, A. J., Bossdorf, O., Hoch, G. G., Lexer, C., et al. (2015). The response of the alpine dwarf Shrub Salix herbacea to altered snowmelt timing: lessons from a multi-site transplant experiment. PLoS One 10:e0122395. doi: 10.1371/journal.pone.0122395
Semenchuk, P. R., Elberling, B., Amtorp, C., Winkler, J., Rumpf, S., Michelsen, A., et al. (2015). Deeper snow alters soil nutrient availability and leaf nutrient status in high Arctic tundra. Biogeochemistry 124, 81–94. doi: 10.1007/s10533-015-0082-7
Sigdel, S. R., Pandey, J., Liang, E., Muhammad, S., Babst, F., Leavitt, S. W., et al. (2021). No benefits from warming even for subnival vegetation in the central Himalayas. Sci. Bull. 66, 1825–1829. doi: 10.1016/j.scib.2021.06.005
Steger, C., Kotlarski, S., Jonas, T., and Schär, C. (2013). Alpine snow cover in a changing climate: a regional climate model perspective. Clim. Dyn. 41, 735–754. doi: 10.1007/s00382-012-1545-3
Steinbauer, M. J., Grytnes, J.-A., Jurasinski, G., Kulonen, A., Lenoir, J., Pauli, H., et al. (2018). Accelerated increase in plant species richness on mountain summits is linked to warming. Nature 556, 231–234. doi: 10.1038/s41586-018-0005-6
Taschler, D., and Neuner, G. (2004). Summer frost resistance and freezing patterns measured in situ in leaves of major alpine plant growth forms in relation to their upper distribution boundary. Plant Cell Environ. 27, 737–746. doi: 10.1111/j.1365-3040.2004.01176.x
Thalmann, M., and Santelia, D. (2017). Starch as a determinant of plant fitness under abiotic stress. New Phytol. 214, 943–951. doi: 10.1111/nph.14491
Thomas, P. A., El-Barghathi, M., and Polwart, A. (2007). Biological Flora of the British Isles: juniperus communis L. J. Ecol. 95, 1404–1440. doi: 10.1111/j.1365-2745.2007.01308.x
Thomas, P. W., Woodward, F. I., and Quick, W. P. (2004). Systemic Irradiance Signalling in Tobacco. New Phytol. 161, 193–198. doi: 10.1046/j.1469-8137.2003.00954.x
Tian, M., Yu, G., He, N., and Hou, J. (2016). Leaf morphological and anatomical traits from tropical to temperate coniferous forests : mechanisms and influencing factors. Sci. Rep. 6:19703. doi: 10.1038/srep19703
Tumajer, J., Buras, A., Camarero, J. J., Carrer, M., Shetti, R., Wilmking, M., et al. (2021). Growing faster, longer or both? Modelling plastic response of Juniperus communis growth phenology to climate change. Glob. Ecol. Biogeogr. 30, 2229–2244. doi: 10.1111/geb.13377
Unterholzner, L., Carrer, M., Bär, A., Beikircher, B., Dämon, B., Losso, A., et al. (2020). Juniperus communis populations exhibit low variability in hydraulic safety and efficiency. Tree Physiol. 40, 1668–1679. doi: 10.1093/treephys/tpaa103
Walker, M. D., Wahren, H. C., Hollister, R. D., Henryd, G. H. R., Ahlquist, L. E., Alatalo, J. M., et al. (2006). Plant community responses to experimental warmingacross the tundra biome. Proc. Natl. Acad. Sci. U. S. A. 103, 1342–1346. doi: 10.1073/pnas.0503198103
Wheeler, J. A., Cortés, A. J., Sedlacek, J., Karrenberg, S., van Kleunen, M., Wipf, S., et al. (2016). The snow and the willows: earlier spring snowmelt reduces performance in the low-lying alpine shrub Salix herbacea. J. Ecol. 104, 1041–1050. doi: 10.1111/1365-2745.12579
Wheeler, J. A., Hoch, G., Cortés, A. J., Sedlacek, J., Wipf, S., and Rixen, C. (2014). Increased spring freezing vulnerability for alpine shrubs under early snowmelt. Oecologia 175, 219–229. doi: 10.1007/s00442-013-2872-8
Wipf, S., and Rixen, C. (2010). A review of snow manipulation experiments in Arctic and alpine tundra ecosystems. Polar Res. 29, 95–109. doi: 10.1111/j.1751-8369.2010.00153.x
Wipf, S., Stoeckli, V., and Bebi, P. (2009). Winter climate change in alpine tundra: plant responses to changes in snow depth and snowmelt timing. Clim. Change 94, 105–121. doi: 10.1007/s10584-009-9546-x
Wood, S. N. (2017). Generalized Additive Models: an introduction with R, 2nd Edn. London: Chapman and Hall/CRC.
Keywords: snowmelt, alpine tundra, shrub phenology, Juniperus communis, leaf traits, non-structural carbohydrates, primary growth, climate change
Citation: Unterholzner L, Prendin AL, Dibona R, Menardi R, Casolo V, Gargiulo S, Boscutti F and Carrer M (2022) Transient Effects of Snow Cover Duration on Primary Growth and Leaf Traits in a Tundra Shrub. Front. Plant Sci. 13:822901. doi: 10.3389/fpls.2022.822901
Received: 26 November 2021; Accepted: 14 February 2022;
Published: 05 April 2022.
Edited by:
Martina Pollastrini, University of Florence, ItalyReviewed by:
Cristina Nabais, University of Coimbra, PortugalShalik Ram Sigdel, Institute of Tibetan Plateau Research (CAS), China
Julia Wheeler, Agriculture and Agri-Food Canada (AAFC), Canada
Copyright © 2022 Unterholzner, Prendin, Dibona, Menardi, Casolo, Gargiulo, Boscutti and Carrer. 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) and the copyright owner(s) 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: Lucrezia Unterholzner, lucrezia.unterholzner@phd.unipd.it