Abstract
In tundra ecosystems, snow cover protects plants from low temperatures in winter and buffers temperature fluctuations in spring. Climate change may lead to reduced snowfall and earlier snowmelt, potentially exposing plants to more frequent and more severe frosts in the future. Frost can cause cell damage and, in combination with high solar irradiance, reduce the photochemical yield of photosystem II (ΦPSII). Little is known about the natural variation in frost exposure within individual habitats of tundra plant populations and the populations’ resilience to this climatic variation. Here, we assessed how natural differences in snowmelt timing affect microclimatic variability of frost exposure in habitats of the evergreen Vaccinium vitis-idaea in sub-Arctic alpine Finland and whether this variability affects the extent of cell damage and reduction in ΦPSII. Plants in early melting plots were exposed to more frequent and more severe frost events, and exhibited a more pronounced decrease in ΦPSII, during winter and spring compared to plants in late-melting plots. Snowmelt timing did not have a clear effect on the degree of cell damage as assessed by relative electrolyte leakage. Our results show that sub-Arctic alpine V. vitis-idaea is currently exposed to strong climatic variation on a small spatial scale, similar to that projected to be caused by climate change, without significant resultant damage. We conclude that V. vitis-idaea is effective in mitigating the effects of large variations in frost exposure caused by differences in snowmelt timing. This suggests that V. vitis-idaea will be resilient to the ongoing climate change.
Similar content being viewed by others
Introduction
Climate change in northern regions will lead to rising temperatures in winter and changes in the amount of snowfall and the timing of snowmelt (Callaghan et al. 2011; AMAP 2017). In tundra ecosystems, snow cover protects plants from low and fluctuating temperatures, so an advance in snowmelt timing caused by climatic warming potentially exposes plants to more frequent early season frost events. Frost (sub-freezing temperatures), especially during the period of dehardening and growth initiation in spring, can cause cell damage to vegetative and flower buds and stress to the photosystem with consequences at the ecosystem level for pollinators, seed recruitment, growth and carbon cycling (Ensminger et al. 2004; Taulavuori et al. 2004; Inouye 2008; Wipf et al. 2009).
Snow manipulation experiments have been used to study the effect of simulated changes in snowmelt timing on plant physiology (Starr 2001; Starr et al. 2008; Lundell et al. 2010; Taulavuori et al. 2011). These effects can be studied alternatively along natural gradients in snowmelt timing that occur on the scale of a few metres in sub-Arctic alpine ecosystems, due to the effect of wind blow and the uneven topography on the distribution of snow. (Billings and Bliss 1959; McKay and Gray 1981; Ögren 1996; Bokhorst et al. 2010, 2017). Spatial heterogeneity in snow cover on a small local scale implies that plants growing in early melting vegetation patches may be exposed to early season frost events, whilst thicker snow cover on late-melting patches insulates plants from low air temperatures until later in the season when temperatures are milder (Kimball et al. 1973). This variation in plant exposure is not reflected in those air temperature measurements used as standard in species distribution models (SDMs). This is because (1) SDMs usually use global climate data with a coarse resolution which does not capture small-scale spatial heterogeneity and (2) such climate data are usually recorded 2 m above the ground, which does not account for the temperature in the soil or near the ground, especially under the snow (Scherrer and Körner 2010, 2011; Scherrer et al. 2011; Graae et al. 2012; Aalto et al. 2017). Heterogeneous landscapes can harbour future refugia or stepping stones which allow plants to persist in a changing climate by providing suitable microhabitats for survival and migration as the macroclimatic conditions change (Niskanen et al. 2017a, b; Graae et al. 2018). Furthermore, it has been suggested that plant populations exposed to highly variable environments are likely to hold large adaptive potential, which may help them to adapt to new environments in the future (Crawford 2008). Hence, by overlooking microclimatic variability, SDMs risk underestimating the role of phenotypic plasticity for the persistence of a species in a changing climate (Lenoir et al. 2017; Lembrechts et al. 2018).
The risk of frost exposure and the persistence of snow are important determinants of local plant performance in cold ecosystems, as well as influencing regional plant distribution patterns (Inouye 2000; Körner 2003; Niittynen and Luoto 2017). Despite this, small-scale spatial variation in frost exposure at the plant level is not commonly reported in snow gradient studies, even though it can vary significantly within microhabitats. In order to assess the likelihood that certain plant species will be able to persist in their current habitats under future climates, or persist locally in nearby niches, we need information on how much environmental variation is present in their current habitats and how plant populations cope with this variation.
Some species, like the evergreen dwarf shrub Vaccinium vitis-idaea L., exist in a wide range of microclimates and hence a single population may experience a multitude of climatic conditions, which is reflected in its high phenotypic plasticity (Ritchie 1955). With its distribution across polar, boreal and alpine ecosystems, V. vitis-idaea is an important species of ecosystems that are largely snow-covered in winter. Wintertime photosynthesis beneath the snow is well-documented for evergreen dwarf shrubs such as V. vitis-idaea, Ledum palustre and Cassiope tetragona (Starr and Oberbauer 2003; Lundell et al. 2008, 2010; Saarinen and Lundell 2010). Underneath the snow, the conditions can be suitable for photosynthesis when sufficient light penetrates to the leaves and the temperature and moisture conditions are favourable (Neuner et al. 1999; Tahkokorpi et al. 2007; Lundell et al. 2008). Wintertime photosynthesis can even compensate for respiratory losses, as the light compensation point in winter is low (Starr and Oberbauer 2003; Lundell et al. 2008; Saarinen and Lundell 2010; Saarinen et al. 2015). However, compared to the subnivean environment, upon snowmelt in spring, plants are suddenly exposed to high solar radiation and low temperature (Gerland et al. 2000). Due to the low temperature inhibition of photosynthetic electron transport and carbon dioxide uptake, this can create an imbalance between the temperature-insensitive absorption of light energy and the downstream photosynthetic reactions (Öquist and Huner 2003). When more light energy is absorbed by chlorophyll than can be passed down the electron transport chain, the excess energy can cause photooxidative damage (Adams III et al. 2002; Öquist and Huner 2003). Plants possess several mechanisms and structures to avert this risk during unfavourable winter conditions, such as the thermal dissipation (non-photochemical quenching, NPQ) of the photon energy by the photosystem II protein complex and light-harvesting antennae aided by zeaxanthin (Li et al. 2000). This mechanism can be activated transiently upon exposure to light during unfavourable conditions (e.g. low temperatures) or it develops during winter hardening (sustained NPQ) and is released in spring upon return to favourable weather conditions (Adams III et al. 2004). Damage to the photosystem II and NPQ both result in a reduced photochemical yield of photosystem II (ΦPSII) (Adams III et al. 2002, 2004). Lack of snow cover in winter combined with both low temperatures and high irradiance causes a decrease in ΦPSII in sub-Arctic V. vitis-idaea, compared to when it is growing under snow cover (Saarinen and Lundell 2010; Taulavuori et al. 2011). However, we do not know whether the degree of inhibited ΦPSII varies amongst individuals of V. vitis-idaea growing along a natural snowmelt gradient where the plants are exposed to differences in low temperature and high irradiance.
Temperatures below the threshold for freezing resistance can directly damage plant cell membranes. Frost hardiness, or frost survival, in plants changes seasonally and depends on the preceding temperature and moisture conditions, but winter snow cover reduces hardiness (Palacio et al. 2015; Sklenář 2017). The degree of sub-lethal cell damage can be measured as the extent of electrolytes leaking from damaged cells (Bajji et al. 2001). It is not known whether cell damage differs amongst individuals of V. vitis-idaea from adjacent microclimate regimes.
In this study, we assessed the extent of spatial and temporal microclimate variability of snowmelt timing and frost exposure and its effect on the leaf photosynthetic function and cell damage in V. vitis-idaea in a sub-Arctic environment above the treeline in northern Finland. We hypothesised that
-
(1)
Frost intensity, in terms of the duration and minimum temperature of frost events, increases with earlier day of year (DOY) of snowmelt,
-
(2)
V. vitis-idaea experiences more cell damage and reduced ΦPSII in spring shortly after snowmelt compare to summer and
-
(3)
Cell damage and reduction in ΦPSII in spring increase with earlier DOY of snowmelt.
We also evaluate whether frost intensity can explain these ecophysiological responses and addressed these hypotheses under consideration of the potentially confounding effects of mountain slope aspect and year of the measurements.
Materials and methods
Study site
The study was carried out in sub-Arctic alpine tundra on Saana mountain in Kilpisjärvi, Finland (69° 02′ 37″ N, 20° 51′ 22″ E) in 2015 and 2016 (Online Resource 1a). For the period 1981–2010, mean monthly air temperatures range from − 12.9 °C (January) to 11.2 °C (June) and the mean annual precipitation is 487 mm (Pirinen et al. 2012). Snowmelt begins in April–May but can last until July in patches with thick snow cover and snow depth can vary by up to 256 cm within a few metres (Aalto et al. 2018). The vegetation in the study area consists primarily of woody dwarf shrubs such as Empetrum nigrum, Betula nana, Salix spp. and Vaccinium spp. On two mountain slope aspects (‘Sites’: north-west facing, termed ‘North’ and south-east facing, termed ‘South’), we identified areas with different timing of snowmelt in spring and classified them into three groups based on their sequential melt-out time with respect to each other to cover a wide range of snowmelt timing (Online Resource 1b). These snowmelt categories were termed ‘Early’, ‘Mid’ and ‘Late’ and located at approximately 700–770 m a.s.l. Within each site, four plots of 2 m × 2 m were marked per snowmelt group in 2014 for an earlier study (Gehrmann et al. 2018). In total, 21 sample plots were used in the present study with 12 plots on the North site (four in Early, four in Mid and four in Late) and 9 plots on the South site (four in Early, four in Mid and one in Late). We could only select one plot in South-Late as V. vitis-idaea did not occur in the other plots. Plots within the same Site were between 3 and 200 m apart. A detailed description of the site and study design can be found in Gehrmann et al. (2018).
Temperature measurements and frost severity
Air temperature at approximately 15 cm above the ground, soil surface temperature and temperature 10 cm below the ground (rooting zone) were recorded at hourly intervals with one or two iButtons® (DS1922L-F5 thermochrons; Homechip Ltd., UK) in each plot, so that each plot had between three and six temperature loggers. The temperatures were recorded continuously from spring 2014 until the end of summer in 2016. The air temperature loggers were protected from direct sunlight with cylindrical radiation shields (height = 15 cm, diameter = 6 cm) made of white PVC plastic.
The timing of snowmelt was estimated as the day on which the range of daily fluctuations in soil surface temperatures reached at least 5 °C. Precise measurements of snow depth were not possible due to the risk of avalanches in the area, but early melting plots had the least amount of snow and late-melting plots the most.
In order to assess frost intensity, for each plot, each continuous sequence of hourly temperatures below 0 °C was extracted and constituted one freezing event in the analyses. Multiple events per day were possible. For each freezing event, the duration (number of hours) of the event and its minimum temperature were determined.
Photosynthetic photo flux density
Photosynthetic photon flux density (PPFD, quantified as the amount of photons in µmol in the range of 400–700 nm reaching a surface of 1 m2 per second) for Kilpisjärvi was obtained from the STRÅNG model of the Swedish Meteorological and Hydrological Institute (https://strang.smhi.se/), as locally measured data were not available. The model produces hourly PPFD data at a resolution of 2.5 km × 2.5 km based on calculations of clear sky spectral irradiance coupled with a cloud effect function. The input data for the model are modelled meteorological parameters and satellite ozone observations. For each freezing event, we extracted the highest PPFD (PPFDmax) modelled during that event. The PPFDmax during each freezing event was determined to find out whether the number of high irradiance and low temperature events varied along the snowmelt gradient. Based on prior experiments (Lundell et al. 2008), we considered PPFDmax > 500 µmol m−2 s−1 during spring freezing events to be a potentially stressful high-light exposure.
ΦPSII measured using chlorophyll fluorescence
The photochemical quantum yield of photosystem II (dark-adapted ΦPSII = Fv/Fm) in V. vitis-idaea was assessed by measuring the ratio of variable (Fv) to maximum fluorescence (Fm) in dark-adapted leaves. This ratio is a measure of photosystem II stress. The measurements were made in 2015 and 2016 once shortly after snowmelt (‘Spring’) and once in late July (‘Summer’, Table 1). In each plot, 24 healthy looking leaves of different shoots from the previous growing season (the second or third leaf from the tip of the shoot) were dark-adapted in leaf clips for 20 min before measuring Fv/Fm with a portable chlorophyll fluorometer (PAM-2100, Walz GmbH). A few measurements where the chlorophyll fluorescence signal was too noisy to confidently determine the Fv/Fm were excluded.
Freezing damage estimated by relative electrolyte leakage
Relative electrolyte leakage is a measure of cell damage (Sutinen et al. 1992; Luoranen et al. 2004). The relative electrolyte leakage of V. vitis-idaea leaves was measured twice both in 2015 and 2016 (Table 1). The first set of measurements was carried out within two weeks after snowmelt (‘Spring’), except in early melting plots in 2015 and late-melting plots in 2016 in the South, where the measurements were made three to four weeks after snowmelt. In each year, the second set of measurements was carried out in July (‘Summer’) on the same or consecutive days within the same site. In late-melting plots, no second set of measurements was made in 2015.
From each plot, we collected 72 healthy looking leaves from the previous growing season (the second or third leaf from the tip of the shoot; different leaves to those used for chlorophyll fluorescence) evenly across the plot in sample boxes containing moist tissue paper and stored them at 4 °C until the measurements were carried out on the same or the following day. We pooled the leaves into samples of three leaves (i.e. 24 samples per plot) and cut each leaf into six evenly sized pieces. Each sample was placed in a separate 50 mL Eppendorf tube and rinsed with approximately 15 mL distilled water. The leaves were then incubated in 20 mL of fresh distilled water and shaken at 100 rpm for 14 h. Subsequently, the tubes were boiled for 60 min to kill the plant cells. The conductivity of each sample was measured with a conductivity probe (Hanna Instruments, Inc.) before and after boiling, the latter after the samples had reached room temperature. To account for spontaneous ionisation of pure water, the conductivity of control samples containing only 20 mL of deionised water was subtracted from each sample measurement. The relative electrolyte leakage was calculated as the ratio of conductivity before boiling to that after boiling and converted to percentage. Further details of the method can be found in Lundell et al. (2008).
Statistical analyses
For the following statistical analyses, we used the actual date of snowmelt timing (SnowmeltDOY) of each plot instead of the category of snowmelt (Early, Mid, Late) to achieve a better representation of the continuous gradient in snowmelt timing at the study site.
We tested whether the duration and minimum temperature of each freezing event varied significantly between years (2015, 2016), sites (North, South), seasons (spring, summer) and the DOY of snowmelt timing (SnowmeltDOY) using a linear mixed model with Plot as a random effect and a weighted variance for SnowmeltDOY, the latter to improve homogeneity of residuals. We used the log-likelihoods to compare nested models and only retained variables in the final model which significantly (p < 0.05) improved the likelihood ratio. The significance of individual covariates was assessed with an F-test applied to the restricted maximum likelihood estimates of the model.
Akaike’s information criterion (AIC) was used to evaluate whether Season (Spring and Summer), Site (North and South) or SnowmeltDOY or a combination of those, were good predictors of relative electrolyte leakage and of Fv/Fm. For this, we fitted mixed effect maximum likelihood models with either of the three factors, or interactions amongst them, as explanatory variables. Relative electrolyte leakage (REL [%]) or Fv/Fm, was the response variable and Plot was a random effect in the model. The model with the lowest AIC value was considered to have the best fit. Models were fitted to each year separately due to different timing of the measurements during the two years. An F-test was applied to the best-fit model (with restricted maximum likelihood estimates) to test the significance (p < 0.05) of the effects of the individual variables and their interactions. Pairwise comparisons were carried out with the emmeans package (Lenth 2020, version 1.4.6). The assumptions of normality, independence and homogeneity of variance of the residuals was confirmed with visual assessment of standard residual plots. A weighted variance function was applied to one or several factors if these assumptions were not met. We used the lme function in the nlme package version 3.1–147 (Pinheiro et al. 2020) for fitting the models. All analyses and figures were made in R version 4.0.0 (Pinheiro et al. 2020; R Core Team 2020).
Results
Snowmelt timing and temperature conditions at the plant level
In 2015, the early melting plots became snow-free between 22 April and 26 May, the mid-melting plots between 26 May and 06 June and the late-melting plots between 08 and 29 June (Table 1). In 2016, an unusually warm period at the end of April caused plots classified as early or mid melting to melt out at about the same time, between 28 April and 06 May (Table 1). Late-melting plots became snow-free between 15 May and 11 Jun in 2016. South plots melted out before North plots within the same snowmelt timing.
In early and mid-melting plots, air temperatures at 15 cm above the ground ranged between − 5 and − 20 °C and fluctuated by more than ten degrees within a single day during the winter months (January–May; Fig. 1, only data from North site shown). Late-melting plots retained air temperatures at 15 cm above the ground between 0 and − 5 °C from January until the snow melted, with the exception of a few colder days (Fig. 1). In comparison, air temperature measured at the weather station in Kilpisjärvi (ca. 1.5 km from the study site and 500 m a.s.l.) fluctuated more than those at our plots in the winter, reaching a minimum of − 37 °C, and weather station temperature was generally lower in the summer than the canopy-level air temperature measurements (Fig. 1).
In early and mid-melting plots, the soil surface temperature followed the fluctuations of the air temperature during January–May, with the coldest temperature being − 13 °C (Fig. 1). In the late-melting plots, the soil surface temperature was always above − 2 °C, apart from the cold days in February where it was − 3.7 °C at the coldest. Soil surface temperature in late-melting plots was decoupled from the air temperature fluctuations until snowmelt.
In late-melting plots, especially in the South, the rooting zone temperature was often above freezing already before the snow melted and it never dropped below − 5 °C after January (data from South site not shown). In all of the other plots, the rooting zone temperature fluctuated considerably during the winter and reached down to − 14 °C, with lower temperatures occurring usually in plots on the North site. In all plots, rooting zone temperatures remained around 0 °C during melt-out due to the constant heat transfer between meltwater and the frozen soil water, termed the zero curtain effect (Outcalt et al. 1990). A fast rise in the rooting zone temperatures above 0 °C, alongside its increased fluctuation above freezing, took place simultaneously with the estimated dates of snowmelt timing (Fig. 1 and Table 1).
Frost intensity
In line with our first hypothesis, frost intensity decreased significantly towards later snowmelt dates. From earliest to latest snowmelt DOY, the mean number of frost events decreased from 43 to 0 (data not shown), the frost event duration from 17 to 0 h (SnowmeltDOY F1,710 = 133.096, p < 0.0001, Fig. 2a) and the mean minimum temperature of frost events from − 11.3 to − 0.2 °C (SnowmeltDOY F1,711 = 239.0232, p < 0.0001, Fig. 2b). Freezing duration was significantly higher in the South site compared to the North site, although the difference was small and only present in spring (effect size 1.3 h, p = 0.0178). In the South site, freezing duration in summer was 1.5 h less than in spring (p < 0.0001). Apart from SnowmeltDOY, minimum freezing temperature was only significantly affected by season where it was marginally higher in the summer (effect size 0.3 °C, p = 0.0001).
Exceptionally long and cold freezing events occurred on the south-facing slope in early melting plots in 2015 shortly after snowmelt (Fig. 3). Freezing events there lasted up to 17 h with a minimum temperature of − 11.3 °C. In 2016, early and mid-melting plots were exposed to much more similar freezing events than in 2015, as both melted out on similar dates (Fig. 3).
Concomitant high irradiance and low temperature events
During the freezing events occurring after snowmelt, the highest received irradiance (PPFDmax) was 685 µmol m−2 s−1 (Online Resource 2). Some freezing temperatures occurred during the night early in spring, so the lowest PPFDmax was 0 µmol m−2 s−1. Low temperature events with concomitant high irradiance were only present in early melting plots in the South site in 2015 (marked with a red asterisk in Fig. 3).
Chlorophyll fluorescence
Chlorophyll fluorescence was best explained by the interaction of Season with Site and Season with SnowmeltDOY in 2015 and by the interaction of Season, Site and SnowmeltDOY in 2016 (Online Resource 3). In accordance with our hypothesis 2, the combined mean Fv/Fm (chlorophyll fluorescence of photosystem II) in V. vitis-idaea was higher in the summer (0.75 ± 0.01, n = 42) than in spring (0.40 ± 0.01, n = 42) (Fig. 4). This difference was statistically significant in both 2015 and 2016 as a main effect (Table 2) and within each Site (Online Resource 3). In 2015, Fv/Fm increased by 0.004 per day of delay in snowmelt in spring and increased correspondingly by 0.0002 per day of delay in snowmelt in summer, which was according to our hypothesis 3. In 2016, Fv/Fm also increased with later snowmelt in spring, by 0.002 and 0.012 per day of later snowmelt in the North and South site, respectively, but in the South site in summer it decreased by 0.002 per day of later snowmelt delay (Table 2, Fig. 4).
Relative electrolyte leakage
Relative electrolyte leakage (REL) was best explained by the two-way interactions of Season, Site and SnowmeltDOY in 2015 and by Season in 2016 (Online Resource 3). Electrolyte leakage was always significantly higher in summer compared to spring, but only by 2.4% (F1,875 = 4.4255, p = 0.0357) and 0.7% (F1,736 = 309.2752, p = < 0.0001) in 2015 and 2016, respectively (Fig. 5, Table 2, Online Resource 3). The effect of snowmelt timing (SnowmeltDOY) found in 2015 was opposite in the two sites. Relative electrolyte leakage decreased with later snowmelt timing in the North site and increased with later snowmelt timing in the South site. These results only partially agree with our hypotheses 2 and 3, that REL would be higher in spring than in summer (the measured difference was very small) and higher in spring after earlier compared to later snowmelt (this was only observed in 2015 in the North site).
Discussion
Microclimate determines frost exposure on a small spatial scale
Our small-scale temperature measurements along a snowmelt gradient show that spatial variability in the microclimate is high in this tundra ecosystem and is decoupled from the regional conditions. This leads to a large discrepancy between the air temperatures in the winter near the plant canopy under snow and air temperatures measured at the weather station. The microtopography of the landscape causes large gradients in snowmelt timing on the scale of a few metres, which exposes plants to a spatially highly variable environment (Scherrer and Körner 2010; Aalto et al. 2017; Graae et al. 2018). The earliest melting plots melted out up to seven weeks before the latest melting plots, which is similar to the 20% decrease in the length of snow cover duration predicted for this region in 2071–2100 (Jylhä et al. 2009). Along the snowmelt gradient, the number of freezing events after snowmelt decreased by 95% from the earliest to the latest melting plots, with some plots being exposed to no freezing events at all after snowmelt (Fig. 3). This difference in freezing events is considerably higher than the 30% reduction in the number of freezing days predicted for the period 2071 to 2100 in northern Finland (Jylhä et al. 2009). In a snow manipulation experiment in sub-Arctic Alaska, the number of frost days decreased by only 45% in late-melting plots compared to control plots (Wipf et al. 2006), which is not similar to the natural variation we found here. Our study also shows that the exposure to frost largely depends on the timing of snowmelt, which is highly dependent on microtopography, snowfall and spring temperature. Changes in the timing of snowmelt are difficult to predict, especially due to the uncertainty in expected cloud cover (Rosenfeld et al. 2014), and because microtopography is highly heterogeneous in alpine areas, the actual future frost exposure of sub-Arctic alpine plants may not correspond to the general regional predictions. This highlights the need for detailed quantitative data of microclimatic variation along local environmental gradients (Graae et al. 2012), such as we have reported here.
The risk of frost exposure is an important determinant of plant distribution in cold ecosystems. Plants occurring in these ecosystems have acquired traits which either confer avoidance or tolerance of frost damage (Inouye 2000; Körner 2003). Vaccinium vitis-idaea has successfully established itself along the whole snowmelt gradient in our study site, so it must be well adapted to cope with the present variation in frost risk. In this species, the phenological development starts late compared to co-occurring species, and this may represent a strategy allowing the evergreen leaves to cope with early season frost events and its new leaves and flowers to emerge only once frost events are rarer and less severe (Gehrmann et al. 2018). Actively differentiating and young tissues are usually more sensitive to frost than older tissues, so delaying growth until conditions are benign can be an advantage (Inouye 2000).
Apart from a stable, topographically controlled microclimate, sudden weather changes can also temporarily affect the microclimate at a particular site. Warming events during winter are common across Arctic and sub-Arctic areas and are predicted to become more frequent and more severe in the future (see Vikhamar-Schuler et al. 2016 and references therein). During our study period in 2016, one of these unusual winter warming events occurred. The air temperature rose above 5 °C between the end of April and the beginning of May and reached 8.9 °C during this period (Fig. 1). In comparison, the average maximum April temperature during the normal period (1981–2010) is 0 to 2 °C for the region (Pirinen et al. 2012) and during the same period in April 2015 the temperature remained below 4 °C. The warming event caused rapid snowmelt, which led our mid-melting plots to melt out unusually early in 2016, creating a similar microclimate as in the early melting plots. Warming events like this can lead to the sudden exposure of plants to freezing temperatures when the usual seasonal cold returns, as was the case in the early and mid-melting plots in our study. In vascular plants, the effect of such events varies with species, but in general, plants show reduced flower abundance, reduced reproductive success, reduced shoot and root growth and delayed bud burst (Bokhorst et al. 2009, 2011; Semenchuk et al. 2013). In one study, V. vitis-idaea suffered 51% shoot mortality after simulated winter warming events, but its growth, reproduction and phenology were not affected (Bokhorst et al. 2011). In our study, V. vitis-idaea usually had delayed vegetative and reproductive bud burst in mid- compared to early-melting plots, but the extreme warming event in 2016 advanced the development in mid-melting plots compared to previous years, leading to synchronised bud burst across these plots in 2016 (Gehrmann et al. 2018). Despite the advanced phenological development in the mid-melting plots in 2016, we did not see any visible damage to V. vitis-idaea in these plots.
Dynamics of ΦPSII after snowmelt
Springtime decrease in ΦPSII in V. vitis-idaea was evident along the entire snowmelt gradient as shown by the low Fv/Fm ratio shortly after snowmelt in spring compared to the Fv/Fm ratio in July (Fig. 4). In V. vitis-idaea, the photochemical yield of PSII is typically high in leaves under snow cover in winter and starts declining when the snow depth has decreased to less than 20 cm, indicating that there is no sustained wintertime NPQ in response to frost hardening (Lundell et al. 2008). Upon snowmelt, the photochemical yield of PSII drops substantially and exhibits a slow recovery over several weeks, which is thought to be associated with increased thermal dissipation in response to cold temperatures and high-light exposure (Lundell et al. 2008; Solanki et al. 2019). Therefore, the springtime drop in ΦPSII, which we observed a few days up to a couple of weeks after snowmelt along the whole gradient, likely reflects a strategy of V. vitis-idaea to adjust the photosynthetic apparatus to the early spring conditions (Huner et al. 1993). This indicates that V. vitis-idaea successfully mitigates the exposure to unfavourable spring conditions regardless of snowmelt timing, allowing full recovery of the potential yield of PSII to c. 0.80 by the summer (Maxwell and Johnson 2000; Starr 2001; Lundell et al. 2008).
In a review of photoinhibition (one of the components of NPQ) in plants from different environments, Míguez et al. (2015) found a temperature dependency of cold-induced photoinhibition during winter, suggesting that it is determined by minimum temperatures specifically rather than freezing days in general. In shrubs in particular, ΦPSII was reduced by 75% when the mean winter temperature was below 0 °C. This is comparable to the approximately 50% reduction in ΦPSII which we found in the early melting plots of the snowmelt gradient in the spring in connection with the occurrence of freezing events, as compared to the ΦPSII measured in summer, when temperatures were already above 0 °C. In the later melting plots, the high Fv/Fm in spring was related to few or no freezing events so that only a 40% reduction compared to summer values was observed. In the South site in 2015, late-melting plots had the lowest Fv/Fm in spring. Lundell et al. (2010) have shown that photoinhibition in V. vitis-idaea in the same region already begins when the snow cover becomes thinner due to melting and this decrease in Fv/Fm continues during the first four days after snowmelt. This was the likely reason also for our observations in the late-melting plots in the South in 2015. The higher Fv/Fm measured in mid-melting plots in the South site in Spring 2015 may be due to the low number of freezing events prior to the measurements, which have enabled faster recovery. This contrasts with occasions where additional freezing events have slowed down the recovery of ΦPSII during spring in the boreal Scots pine and in V. vitis-idaea in southern Finland (Ensminger et al. 2004, 2008; Lundell et al. 2010).
Freezing events during which PPFD was higher than 500 µmol m−2 s−1 did not appear to determine the observed differences in ΦPSII as these events only occurred in early melting plots in the South site in 2015. However, based on Lundell et al. (2008, 2010), the rate of photosynthesis in V. vitis-idaea can become saturated at light intensities much lower than 500 µmol m−2 s−1, so this threshold can be considered a conservative estimate of light intensities which can decrease ΦPSII during freezing temperatures.
Leaf damage
Electrolyte leakage in the North was approximately 10% higher in summer compared to spring of 2015, but otherwise it was similar amongst sites (c. 30%) to that found in other such studies throughout the spring (Fig. 5) (Saarinen and Lundell 2010; Preece and Phoenix 2013). This was contrary to our expectation that electrolyte leakage would be higher in spring than in summer due to the occurrence of spring freezing events, which could cause cell damage. This suggests that approximately 30% is a natural amount of electrolyte leakage in V. vitis-idaea in the sub-Arctic tundra in Fenno-Scandinavia during spring. (Saarinen and Lundell 2010; Preece and Phoenix 2013).
Vaccinium vitis-idaea retains frost resistance to approximately -6 °C even in a dehardened state at the end of August (Taschler and Neuner 2004; Sklenář 2017). However, the temperature of frost resistance is often defined as the temperature at which 50% mortality occurs, so sub-lethal damage is not accounted for in measures of frost resistance. This may be the reason why the electrolyte leakage values we found in the summer were higher than we expected based on the frost resistance estimated from mortality. Furthermore, V. vitis-idaea has been suggested to be frost resistant through the tolerance of extracellular ice and mechanisms involving the C-repeat binding factor (CBF), rather than by modifying the fatty acid composition of cell membranes to increase their rigidity which is a common frost resistance mechanism in other plants, including the evergreen Empetrum nigrum (Bokhorst et al. 2017; Sklenář 2017). Therefore, the high electrolyte leakage values we found in V. vitis-idaea in both spring and summer may be due to the lack of a temperature response of membrane fatty acid composition. Other mechanisms likely act to keep frost damage at a similar level in spring, especially after early snowmelt, compared to summer, in order to avert the negative effects of increased frost events early in the spring.
The potential of Vaccinium vitis-idaea to resist the effects of changes in snowmelt timing
Winter conditions are predicted to become harsher for plants in northern Finland, with less snow, shorter duration of snow cover and more frost exposure by 2099 compared to the normal period (1981–2010, Jylhä et al. 2009; Räisänen and Eklund 2012). The speed of these changes requires that local plant populations either migrate to nearby suitable climate niches or are acclimatised to strongly varying climate conditions. In this study, we addressed these mechanisms simultaneously and show that V. vitis-idaea benefits from a combination of the large microclimate heterogeneity in its habitat and its local acclimation to microclimatic variation.
Our study and other published results reveal that sub-Arctic and alpine habitats comprise large microclimate variability which create a buffer against regional climate change (Scherrer and Körner 2010, 2011; Ford et al. 2013). In a study of grassland species in the UK, such heterogeneity in the microclimate allowed the preservation of sufficiently cold habitats during a rise in temperature over 30 years which allowed plants to persist in the colder microclimates rather than go extinct (Maclean et al. 2015). We found that the differences in snow conditions and frost exposure, which exist at a scale of a few metres in this heterogeneous sub-Arctic alpine landscape, correspond in magnitude to projected climate changes. This suggests that even if snow conditions change dramatically in the future, patches with current day snow conditions might persist, allowing V. vitis-idaea to migrate to a suitable microclimate over a short distance.
Many cold-climate plant species have been subjected to repeated periods of glaciation and warming over their evolutionary history, which led them to retain high intra-specific genetic diversity and phenotypic plasticity. This allows populations and species as a whole to be pre-adapted to different climates (Crawford and Abbott 1994). Our results suggest that V. vitis-idaea does not suffer detrimental effects from the exposure to large gradients in snowmelt timing and frost, as shown by its persistence along the microclimatic gradient and the low variation in cell damage. We also showed that the amount of reduction in potential yield of PSII in the springtime is related to the timing of snowmelt and fully reversible, indicating a large degree of phenotypic plasticity in this population and a high capacity to cope with the prevailing climate conditions.
However, other sub-Arctic plant species with a higher affinity for snow may be more sensitive to early snowmelt than V. vitis-idaea, which survives across a wide range of snow conditions (Wipf et al. 2009; Palacio et al. 2015). Extreme weather events, such as sudden winter warming, pose a greater risk to plants than gradual climate changes, because the fast changes only allow for phenotypic plastic responses instead of evolutionary genetic adaptation. The effect of extreme weather events on plants can range from simple stress responses, such as increased light stress on photosystem II, which we observed in mid-melting plots in the warm spring in 2016, to increased mortality. These responses can lead to decreased biomass and productivity, collectively termed “Arctic browning”, which can eventually result in a severe decrease in net ecosystem exchange (Treharne et al. 2019). Despite the high microclimatic heterogeneity in our study site and the persistence of V. vitis-idaea in this environment, extreme weather events could still have implications for CO2 exchange if they increase in severity and frequency. At the same time, too high plasticity can become maladaptive if it prevents genetic adaptation to a new climate (Benito Garzón et al. 2019). More insight is therefore needed into the plasticity of stress responses, especially along heterogeneous environmental gradients.
References
Aalto J, Riihimäki H, Meineri E, Hylander K, Luoto M (2017) Revealing topoclimatic heterogeneity using meteorological station data. Int J Climatol 37:544–556. https://doi.org/10.1002/joc.5020
Aalto J, Scherrer D, Lenoir J, Guisan A, Luoto M (2018) Biogeophysical controls on soil-atmosphere thermal differences: implications on warming Arctic ecosystems. Environ Res Lett 13:074003. https://doi.org/10.1088/1748-9326/aac83e
Adams WW III, Demmig-Adams B, Rosenstiel TN, Brightwell AK, Ebbert V (2002) Photosynthesis and photoprotection in overwintering plants. Plant Biol 4:545–557. https://doi.org/10.1055/s-2002-35434
Adams WW III, Zarter CR, Ebbert V, Demmig-Adams B (2004) Photoprotective strategies of overwintering evergreens. Bioscience 54:41–49. https://doi.org/10.1641/0006-3568(2004)054[0041:PSOOE]2.0.CO;2
AMAP (2017) Snow, Water, Ice and Permafrost in the Arctic (SWIPA) (2017). Arctic Monitoring and Assessment Programme (AMAP), Oslo
Bajji M, Kinet J-M, Lutts S (2001) The use of the electrolyte leakage method for assessing cell membrane stability as a water stress tolerance test in durum wheat. Plant Growth Regul 00:1–10. https://doi.org/10.1023/A:1014732714549
Benito Garzón M, Robson TM, Hampe A (2019) ΔTraitSDM: species distribution models that account for local adaptation and phenotypic plasticity. New Phytol. https://doi.org/10.1111/nph.15716
Billings WD, Bliss LC (1959) An alpine snowbank environment and its effects on vegetation, plant development and productivity. Ecology 40:388–397
Bokhorst SF, Bjerke JW, Tømmervik H, Callaghan TV, Phoenix GK (2009) Winter warming events damage sub-Arctic vegetation: consistent evidence from an experimental manipulation and a natural event. J Ecol 97:1408–1415. https://doi.org/10.1111/j.1365-2745.2009.01554.x
Bokhorst S, Bjerke JW, Davey MP, Taulavuori K, Taulavuori E, Laine K, Callaghan TV, Phoenix GK (2010) Impacts of extreme winter warming events on plant physiology in a sub-Arctic heath community. Physiol Plant 140:128–140. https://doi.org/10.1111/j.1399-3054.2010.01386.x
Bokhorst S, Bjerke JW, Street LE, Callaghan TV, Phoenix GK (2011) Impacts of multiple extreme winter warming events on sub-Arctic heathland: phenology, reproduction, growth, and CO2 flux responses. Glob Chang Biol 17:2817–2830. https://doi.org/10.1111/j.1365-2486.2011.02424.x
Bokhorst S, Jaakola L, Karppinen K, Edvinsen GK, Mæhre HK, Bjerke JW (2017) Contrasting survival and physiological responses of sub-Arctic plant types to extreme winter warming and nitrogen. Planta 247:635–648. https://doi.org/10.1007/s00425-017-2813-6
Callaghan TV, Johansson M, Brown RD, Groisman PY, Labba N, Radionov V, Barry RG, Bulygina ON, Essery RLH, Frolov DM, Golubev VN, Grenfell TC, Petrushina MN, Razuvaev VN, Robinson DA, Romanov P, Shindell D, Shmakin AB, Sokratov SA, Warren S, Yang D (2011) The changing face of Arctic snow cover: a synthesis of observed and projected changes. Ambio 40:17–31. https://doi.org/10.1007/s13280-011-0212-y
Crawford RMM (2008) Cold climate plants in a warmer world. Plant Ecol Divers 1:285–297. https://doi.org/10.1080/17550870802407332
Crawford RMM, Abbott RJ (1994) Pre-adaptation of Arctic plants to climate change. Bot Acta 107:271–278. https://doi.org/10.1111/j.1438-8677.1994.tb00796.x
Ensminger I, Sveshnikov D, Campbell DA, Funk C, Jansson S, Lloyd J, Shibistova O, Öquist G (2004) Intermittent low temperatures constrain spring recovery of photosynthesis in boreal Scots pine forests. Glob Chang Biol 10:995–1008. https://doi.org/10.1111/j.1365-2486.2004.00781.x
Ensminger I, Schmidt L, Lloyd J (2008) Soil temperature and intermittent frost modulate the rate of recovery of photosynthesis in Scots pine under simulated spring conditions. New Phytol 177:428–442. https://doi.org/10.1111/j.1469-8137.2007.02273.x
Ford KR, Ettinger AK, Lundquist JD, Raleigh MS, Hille Ris Lambers J (2013) Spatial heterogeneity in ecologically important climate variables at coarse and fine scales in a high-snow mountain landscape. PLoS ONE 8:e65008. https://doi.org/10.1371/journal.pone.0065008
Gehrmann F, Hänninen H, Liu C, Saarinen T (2018) Phenological responses to small-scale spatial variation in snowmelt timing reveal compensatory and conservative strategies in subarctic-alpine plants. Plant Ecol Divers 10:453–468. https://doi.org/10.1080/17550874.2018.1428693
Gerland S, Liston GE, Winther J-G, Ørbaek JB, Ivanov BV (2000) Attenuation of solar radiation in Arctic snow: field observations and modelling. Ann Glaciol 31:364–368. https://doi.org/10.3189/172756400781820444
Graae BJ, De Frenne P, Kolb A, Brunet J, Chabrerie O, Verheyen K, Pepin N, Heinken T, Zobel M, Shevtsova A, Nijs I, Milbau A (2012) On the use of weather data in ecological studies along altitudinal and latitudinal gradients. Oikos 121:3–19. https://doi.org/10.1111/j.1600-0706.2011.19694.x
Graae BJ, Vandvik V, Armbruster WS, Eiserhardt WL, Svenning J-C, Hylander K, Ehrlén J, Speed JDM, Klanderud K, Bråthen KA, Milbau A, Opedal ØH, Alsos IG, Ejrnæs R, Bruun HH, Birks HJB, Westergaard KB, Birks HH, Lenoir J (2018) Stay or go – how topographic complexity influences alpine plant population and community responses to climate change. Perspect Plant Ecol Evol Syst 30:41–50. https://doi.org/10.1016/j.ppees.2017.09.008
Huner NPA, Öquist G, Hurry VM, Krol M, Falk S, Griffith M (1993) Photosynthesis, photoinhibition and low temperature acclimation in cold tolerant plants. Photosynth Res 37:19–39. https://doi.org/10.1007/BF02185436
Inouye DW (2000) The ecological and evolutionary significance of frost in the context of climate change. Ecol Lett 3:457–463. https://doi.org/10.1046/j.1461-0248.2000.00165.x
Inouye DW (2008) Effects of climate change on phenology, frost damage, and floral abundance of montane wildflowers. Ecology 89:353–362. https://doi.org/10.1890/06-2128.1
Jylhä K, Ruosteenoja K, Räisänen J, Venäläinen A, Tuomenvirta H, Ruokolainen L, Saku S, Seitola T (2009) Arvioita Suomen muuttuvasta ilmastosta sopeutumistutkimuksia varten - ACCLIM-hankkeen raportti 2009. (The changing climate in Finland: estimates for adaptation studies, ACCLIM project report 2009.). Finnish Meterological Institute, Helsinki
Kimball SL, Bennett BD, Salisbury FB (1973) The growth and development of montane species at near-freezing temperatures. Ecology 54:168–173
Körner C (2003) Alpine plant life: functional plant ecology of high mountain ecosystems. Springer, Berlin
Lembrechts JJ, Nijs I, Lenoir J (2018) Incorporating microclimate into species distribution models. Ecography 42:1–13. https://doi.org/10.1111/ecog.03947
Lenoir J, Hattab T, Pierre G (2017) Climatic microrefugia under anthropogenic climate change: implications for species redistribution. Ecography 40:253–266. https://doi.org/10.1111/ecog.02788
Lenth R (2020) Emmeans: estimated marginal means, aka least-squares means. R package version 1(4):6
Li X-P, Björkman O, Shih C, Grossman AR, Rosenquist M, Jansson S, Niyogi KK (2000) A pigment-binding protein essential for regulation of photosynthetic light harvesting. Nature 403:391–395. https://doi.org/10.1038/35000131
Lundell R, Saarinen T, Åström H, Hänninen H (2008) The boreal dwarf shrub Vaccinium vitis-idaea retains its capacity for photosynthesis through the winter. Botany 86:491–500. https://doi.org/10.1139/B08-022
Lundell R, Saarinen T, Hänninen H (2010) Effects of snowmelt on the springtime photosynthesis of the evergreen dwarf shrub Vaccinium vitis-idaea. Plant Ecol Divers 3:121–130. https://doi.org/10.1080/17550874.2010.497195
Luoranen J, Repo T, Lappi J (2004) Assessment of the frost hardiness of shoots of silver birch (Betula pendula) seedlings with and without controlled exposure to freezing. Can J For Res 34:1108–1118. https://doi.org/10.1139/x03-285
Maclean IMD, Hopkins JJ, Bennie J, Lawson CR, Wilson RJ (2015) Microclimates buffer the responses of plant communities to climate change. Glob Ecol Biogeogr 24:1340–1350. https://doi.org/10.1111/geb.12359
Maxwell K, Johnson GN (2000) Chlorophyll fluorescence: a practical guide. J Exp Bot 51:659–668. https://doi.org/10.1093/jxb/51.345.659
McKay GA, Gray DM (1981) The distribution of snowcover. Principles, processses, management and use. In: Gray DM, Male DH (eds) Handbook of snow. Pergamon Press Toronto, Elmsford, pp 153–190
Míguez F, Fernández-Marín B, Becerril JM, García-Plazaola JI (2015) Activation of photoprotective winter photoinhibition in plants from different environments: A literature compilation and meta-analysis. Physiol Plant 155:414–423. https://doi.org/10.1111/ppl.12329
Neuner G, Ambach D, Aichner K (1999) Impact of snow cover on photoinhibition and winter desiccation in evergreen Rhododendron ferrugineum leaves during subalpine winter. Tree Physiol 19:725–732. https://doi.org/10.1093/treephys/19.11.725
Niittynen P, Luoto M (2017) The importance of snow in species distribution models of arctic vegetation. Ecography 40:1–13. https://doi.org/10.1111/ecog.03348
Niskanen AKJ, Heikkinen RK, Mod HK, Väre H, Luoto M (2017a) Improving forecasts of arctic-alpine refugia persistence with landscape-scale variables. Geogr Ann Ser A 99:2–14. https://doi.org/10.1080/04353676.2016.1256746
Niskanen AKJ, Luoto M, Väre H, Heikkinen RK (2017b) Models of Arctic-alpine refugia highlight importance of climate and local topography. Polar Biol 40:489–502. https://doi.org/10.1007/s00300-016-1973-3
Ögren E (1996) Premature dehardening in Vaccinium myrtillus during a mild winter: a cause for winter dieback? Funct Ecol 10:724–732
Öquist G, Huner NPA (2003) Photosynthesis of overwintering evergreen plants. Annu Rev Plant Biol 54:329–355. https://doi.org/10.1146/annurev.arplant.54.072402.115741
Outcalt SI, Nelson FE, Hinkel KM (1990) The zero-curtain effect: Heat and mass transfer across an isothermal region in freezing soil. Water Resour Res 26:1509–1516. https://doi.org/10.1029/WR026i007p01509
Palacio S, Lenz A, Wipf S, Hoch G, Rixen C (2015) Bud freezing resistance in alpine shrubs across snow depth gradients. Environ Exp Bot 118:95–101. https://doi.org/10.1016/j.envexpbot.2015.06.007
Pinheiro J, Bates D, DebRoy S, Sarkar D, R Core Team (2020) nlme: Linear and Nonlinear Mixed Effects Models. R package version 3.1–147
Pirinen P, Simola H, Aalto J, Kaukoranta J-P, Karlsson P, Ruuhela R (2012) Climatological statistics of Finland 1981–2010. Finnish Meterological Institute, Helsinki
Preece C, Phoenix GK (2013) Responses of sub-arctic dwarf shrubs to low oxygen and high carbon dioxide conditions. Environ Exp Bot 85:7–15. https://doi.org/10.1016/j.envexpbot.2012.06.003
R Core Team (2020) R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna
Räisänen J, Eklund J (2012) 21st Century changes in snow climate in Northern Europe: a high-resolution view from ENSEMBLES regional climate models. Clim Dyn 38:2575–2591. https://doi.org/10.1007/s00382-011-1076-3
Ritchie AJC (1955) Vaccinium vitis-idaea L. J Ecol 43:701–708
Rosenfeld D, Sherwood S, Wood R, Donner L (2014) Climate effects of aerosol-cloud interactions. Science 343:379–381. https://doi.org/10.1126/science.1247490
Saarinen T, Lundell R (2010) Overwintering of Vaccinium vitis-idaea in two sub-Arctic microhabitats: a reciprocal transplantation experiment. Polar Res 29:38–45. https://doi.org/10.1111/j.1751-8369.2010.00152.x
Saarinen T, Rasmus S, Lundell R, Kauppinen O-K, Hänninen H (2015) Photosynthetic and phenological responses of dwarf shrubs to the depth and properties of snow. Oikos 125:001–010. https://doi.org/10.1111/oik.02233
Scherrer D, Körner C (2010) Infra-red thermometry of alpine landscapes challenges climatic warming projections. Glob Chang Biol 16:2602–2613. https://doi.org/10.1111/j.1365-2486.2009.02122.x
Scherrer D, Körner C (2011) Topographically controlled thermal-habitat differentiation buffers alpine plant diversity against climate warming. J Biogeogr 38:406–416. https://doi.org/10.1111/j.1365-2699.2010.02407.x
Scherrer D, Schmid S, Körner C (2011) Elevational species shifts in a warmer climate are overestimated when based on weather station data. Int J Biometeorol 55:645–654. https://doi.org/10.1007/s00484-010-0364-7
Semenchuk PR, Elberling B, Cooper EJ (2013) Snow cover and extreme winter warming events control flower abundance of some, but not all species in high arctic Svalbard. Ecol Evol 3:2586–2599. https://doi.org/10.1002/ece3.648
Sklenář P (2017) Seasonal variation of freezing resistance mechanisms in north-temperate alpine plants. Alp Bot 127:31–39. https://doi.org/10.1007/s00035-016-0174-6
Solanki T, Aphalo PJ, Neimane S, Hartikainen SM, Pieristè M, Shapiguzov A, Porcar-Castell A, Atherton J, Heikkilä A, Robson TM (2019) UV-screening and springtime recovery of photosynthetic capacity in leaves of Vaccinium vitis-idaea above and below the snow pack. Plant Physiol Biochem 134:40–52. https://doi.org/10.1016/j.plaphy.2018.09.003
Starr G (2001) A multiscale assessment of physiological processes in arctic tundra plants under natural and simulated climate change scenarios. Dissertation, Florida International University
Starr G, Oberbauer SF (2003) Photosynthesis of arctic evergreens under snow: implications for tundra ecosystem carbon balance. Ecology 84:1415–1420. https://doi.org/10.1890/02-3154
Starr G, Oberbauer SF, Ahlquist LE (2008) The photosynthetic response of Alaskan tundra plants to increased season length and soil warming. Arctic Antarct Alp Res 40:181–191. https://doi.org/10.1657/1523-0430(06-015)
Sutinen M-L, Palta JP, Reich PB (1992) Seasonal differences in freezing stress resistance of needles of Pinus nigra and Pinus resinosa: evaluation of the electrolyte leakage method. Tree Physiol 11:241–254. https://doi.org/10.1093/TREEPHYS/11.3.241
Tahkokorpi M, Taulavuori K, Laine K, Taulavuori E (2007) After-effects of drought-related winter stress in previous and current year stems of Vaccinium myrtillus L. Environ Exp Bot 61:85–93. https://doi.org/10.1016/j.envexpbot.2007.03.003
Taschler D, 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. https://doi.org/10.1111/j.1365-3040.2004.01176.x
Taulavuori KMJ, Taulavuori EB, Skre O, Nilsen J, Igeland B, Laine KM (2004) Dehardening of mountain birch (Betula pubescens ssp. czerepanovii) ecotypes at elevated winter temperatures. New Phytol 162:427–436. https://doi.org/10.1111/j.1469-8137.2004.01042.x
Taulavuori KMJ, Bauer E, Taulavuori E (2011) Overwintering stress of Vaccinium vitis-idaea in the absence of snow cover. Environ Exp Bot 72:397–403. https://doi.org/10.1016/j.envexpbot.2011.02.009
Treharne R, Bjerke JW, Tømmervik H, Stendardi L, Phoenix GK (2019) Arctic browning: Impacts of extreme climatic events on heathland ecosystem CO2 fluxes. Glob Chang Biol 25:489–503. https://doi.org/10.1111/gcb.14500
Vikhamar-Schuler D, Isaksen K, Haugen JE, Tømmervik H, Luks B, Schuler TV, Bjerke JW (2016) Changes in winter warming events in the nordic arctic region. J Clim 29:6223–6244. https://doi.org/10.1175/JCLI-D-15-0763.1
Wipf S, Rixen C, Mulder CPH (2006) Advanced snowmelt causes shift towards positive neighbour interactions in a subarctic tundra community. Glob Chang Biol 12:1496–1506. https://doi.org/10.1111/j.1365-2486.2006.01185.x
Wipf S, Stoeckli V, Bebi P (2009) Winter climate change in alpine tundra: Plant responses to changes in snow depth and snowmelt timing. Clim Change 94:105–121. https://doi.org/10.1007/s10584-009-9546-x
Acknowledgements
Open access funding provided by Umeå University. The staff at Kilpisjärvi Biological Station provided practical support for the fieldwork and knowledge about the local flora. Che Liu, Sannimari Käärmelahti and Katja Vrabec assisted with the field observations. We thank T. Matthew Robson and Alexey Shapiguzov (Organismal and Evolutionary Biology Research Programme, University of Helsinki) for helpful comments on the manuscript. Miska Luoto (Department of Geosciences and Geography) facilitated the selection of the study sites and provided insightful advice on the project. Pedro J. Aphalo gave advice on statistical testing and R programming. STRÅNG data used here are from the Swedish Meteorological and Hydrological Institute (SMHI), and were produced with support from the Swedish Radiation Protection Authority and the Swedish Environmental Agency. This study was funded by grants from the Doctoral Programme in Plant Science (University of Helsinki), the Finnish Society of Forest Sciences, Nordenskiöld-samfundet, Societas pro Fauna et Flora Fennica and Oskar Öflunds Stiftelse to FG.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Conflict of interest
The authors declare that they have no conflict of interest.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Electronic supplementary material
Below is the link to the electronic supplementary material.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Gehrmann, F., Lehtimäki, IM., Hänninen, H. et al. Sub-Arctic alpine Vaccinium vitis-idaea exhibits resistance to strong variation in snowmelt timing and frost exposure, suggesting high resilience under climatic change. Polar Biol 43, 1453–1467 (2020). https://doi.org/10.1007/s00300-020-02721-3
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s00300-020-02721-3