Functional Ecology. 2024;38:219–232. �  | 219wileyonlinelibrary.com/journal/fec Received: 31 March 2023  | Accepted: 9 October 2023 DOI: 10.1111/1365-2435.14466 R E S E A R C H A R T I C L E Higher vascular plant abundance associated with decreased ecosystem respiration after 20 years of warming in the forest–tundra ecotone Eero Myrsky1,2,3  | Juha Mikola4  | Elina Kaarlejärvi3  | Johan Olofsson5  | Sofie Sjögersten6  | Boris Tupek4  | Minna K. Männistö2  | Sari Stark1 1Arctic Centre, University of Lapland, Rovaniemi, Finland; 2Natural Resources Institute Finland (Luke), Natural Resources Unit, Rovaniemi, Finland; 3Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland; 4Natural Resources Institute Finland (Luke), Helsinki, Finland; 5Department of Ecology and Environmental Sciences, Umeå University, Umeå, Sweden and 6School of Biosciences, University of Nottingham, Loughborough, UK Correspondence Eero Myrsky Email: eero.myrsky@helsinki.fi Funding information Academy of Finland, Grant/Award Number: 310776; Lapin Rahasto, Grant/ Award Number: 40221879 Handling Editor: Ji Chen Abstract 1. The on-going climate warming is promoting shrub abundance in high latitudes, but the effect of this phenomenon on ecosystem functioning is expected to depend on whether deciduous or evergreen species increase in response to warming. 2. To explore effects of long-term warming on shrubs and further on ecosystem functioning, we analysed vegetation and ecosystem CO2 exchange after 20 years of warming in the forest–tundra ecotone in subarctic Sweden. A previous study conducted 9 years earlier had found increased evergreen Empetrum nigrum ssp. hermaphroditum in the forest and increased deciduous Betula nana in the tundra. 3. Following current understanding, we expected continued increase in shrub abun- dance that would be stronger in tundra than in forest. We expected warming to in- crease ecosystem respiration (Re) and gross primary productivity (GPP), with a greater increase in Re in tundra due to increased deciduous shrub abundance, leading to a less negative net ecosystem exchange and reduced ecosystem C sink strength. 4. As predicted, vascular plant abundances were higher in the warmed plots with a stronger response in tundra than in forest. However, whereas B. nana had in- creased in abundance since the last survey, E. hermaphroditum abundance had declined due to several moth and rodent outbreaks during the past decade. In contrast to predictions, Re was significantly lower in the warmed plots irrespec- tive of habitat, and GPP increased marginally only in the forest. The lower Re and a higher GPP under warming in the forest together led to increased net C sink. Re was negatively associated with the total vascular plant abundance. 5. Our results highlight the importance of disturbance regimes for vegetation re- sponses to warming. Climate warming may promote species with both a high capacity to grow under warmer conditions and a resilience towards herbivore outbreaks. Negative correlation between Re and total vascular plant abundance This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited. © 2023 The Authors. Functional Ecology published by John Wiley & Sons Ltd on behalf of British Ecological Society. www.wileyonlinelibrary.com/journal/fec mailto: https://orcid.org/0000-0001-6128-2580 https://orcid.org/0000-0002-4336-2648 https://orcid.org/0000-0003-0014-0073 https://orcid.org/0000-0002-6943-1218 https://orcid.org/0000-0003-4493-1790 https://orcid.org/0000-0003-3466-0237 https://orcid.org/0000-0001-9390-1104 https://orcid.org/0000-0003-4845-6536 mailto:eero.myrsky@helsinki.fi http://creativecommons.org/licenses/by/4.0/ http://crossmark.crossref.org/dialog/?doi=10.1111%2F1365-2435.14466&domain=pdf&date_stamp=2023-11-22 220  |    MYRSKY et al. 1  |  INTRODUC TION The high latitudes of Earth are warming as much as four times faster than the global average (Rantanen et al., 2022), which is drastically changing the plant community structure and functioning of ecosys- tems in these areas (Wookey et al., 2009). The most important plant group responding to warming are the shrubs, which have increased in height as well as expanded their distribution along both latitudes and altitudes—the process generally known as Arctic shrubification or greening (Elmendorf, Henry, Hollister, Björk, Bjorkman, et al., 2012; Elmendorf, Henry, Hollister, Björk, Boulanger-Lapointe, et al., 2012; Myers-Smith et al., 2011, 2020). Vegetation composition plays a major role in ecosystem pro- cesses (Bråthen et al., 2017; Wookey et al., 2009). Increased growing season temperatures together with faster nutrient cycles facilitate shrub growth (Bjorkman et al., 2020; Zamin & Grogan, 2012), leading into greater gross primary productivity (GPP) and carbon uptake from the atmosphere (Cahoon et  al.,  2012; Hobbie & Chapin III,  1998; Shaver et  al., 2007). In the long term, however, warming-induced vegetation change can modify ecosystem CO2 fluxes through several indirect mechanisms over and above the direct effects (Weintraub & Schimel, 2005; Wookey et al., 2009). Taller and darker shrub vege- tation reduces albedo (Blok et al., 2011) and increases accumulation of snowpacks (Way & Lapalme, 2021), which both can lead to higher wintertime soil temperatures (Leffler & Welker, 2013), N mineraliza- tion (DeMarco et al., 2011) and ecosystem respiration (Re) (Nobrega & Grogan, 2007). These effects may accelerate decomposition of soil organic matter (SOM) and lead to increased rates of Re (Nobrega & Grogan, 2007), thus contributing to a greater CO2 release from tun- dra (Biasi et al., 2008; Cahoon et al., 2012; Väisänen et al., 2014). If climate warming accelerates Re more than GPP, net ecosystem ex- change (NEE) is reduced, which induces more carbon losses from the ecosystem (Bradford et al., 2016). However, the amplitude of shru- bification and the dynamics between vegetation and carbon sink are not straightforward but rather a sum of complex site- and time-de- pendent interactions (Mekonnen et al., 2021). For example, although shrubification may increase soil temperatures, the taller and denser vegetation cover may also buffer the soil against warming during peak summer season through increased shading (Blok et al., 2010), and it is not yet fully understood whether shrubification eventu- ally results in positive or negative temperature feedbacks (Way & Lapalme, 2021). It has been hypothesized that one key predictor in the effect of Arctic shrubification on carbon cycling is whether warming increases deciduous or evergreen dwarf shrubs (Vowles & Björk, 2018). The expansion of shrubs has often been attributed to an increase in de- ciduous species such as willows (Salix sp.) and dwarf birch, Betula nana (Elmendorf, Henry, Hollister, Björk, Bjorkman, et  al.,  2012; Elmendorf, Henry, Hollister, Björk, Boulanger-Lapointe, et al., 2012; Myers-Smith et al., 2011). These species can enhance soil respiration through producing easily decomposing litter and having ectomycor- rhizal fungal symbionts that efficiently decompose SOM (Parker et  al., 2021). In contrast, evergreen dwarf shrubs have lower GPP and respiration rates, store greater amount of nutrients in long-lived shoots, leaves and roots and thus decelerate carbon cycles (Vowles & Björk, 2018 and references therein). In particular, the evergreen dwarf shrub Empetrum nigrum ssp. hermaphroditum (hereafter E. hermaphroditum), one of the most common plant species found in Fennoscandian subarctic environments, produces slowly decompos- ing litter with high concentrations of allelopathic compounds that can decelerate soil nutrient and carbon cycles (Bråthen et al., 2017). Increased abundance of E. hermaphroditum in response to climate warming (Stark et al., 2021; Vuorinen et al., 2017) could thus lead to slower rather than faster carbon cycles (Vowles & Björk, 2018). Here, we examine whether longer-term changes in vegetation differ from shorter-term changes, and whether the effects of warm- ing on CO2 exchange differ between tundra and forest due to the different response of deciduous and evergreen dwarf shrubs to warming. For this, we report the responses of plant species com- position and ecosystem CO2 fluxes to 20 years of simulated climate warming, carried out using open top chambers (hereafter OTCs) at tundra and mountain birch forest sites located in the forest–tundra ecotone in northernmost Sweden. Results from the same study sites 9 years earlier showed that although E. hermaphroditum is the domi- nant species in both forest and tundra, its abundance had increased in response to warming in forest only, while B. nana had increased in tundra (Kaarlejärvi et al., 2012), making this an optimal setting to test the influence of these species on the response of CO2 fluxes to warming (Figure 1). We predict that (H1) the effect of warming on dwarf shrub growth will be generally stronger in tundra than forest because the deciduous B. nana is able to effectively respond to improved growing conditions (Bret-Harte et al., 2001) and because shrub growth is po- tentially more temperature limited in tundra than forest (Sundqvist further indicate that the indirect impacts of increased plants on soil microclimate may become increasingly important for ecosystem CO2 exchange in the long run, which adds to the different mechanisms that link warming and CO2 fluxes in northern ecosystems. K E Y W O R D S arctic greening, climate change, CO2 exchange, deciduous dwarf shrubs, evergreen dwarf shrubs, moth outbreaks 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense     |  221MYRSKY et al. et  al.,  2020). As found in earlier warming experiments (Biasi et al., 2008, Väisänen et al., 2014,), we predict that (H2) warming will increase GPP and Re both in the forest and the tundra, but owing to the positive effects of increased B. nana abundance on Re (Vowles & Björk, 2018), the increase in Re will be stronger in tundra than forest. Finally, we predict that (H3) in tundra the effect of warming on Re will be stronger than on GPP, thus leading to a reduced ecosystem C sink strength (manifested as less negative NEE), whereas in forest, no difference, or even increased sink strength will be found. 2  |  MATERIAL S AND METHODS 2.1  |  Study area and experimental design The study area is in the forest–tundra ecotone—a 3- to 4-km wide mosaic of mountain birch forest patches extending to tundra. It is located in the oroarctic region (Virtanen et al., 2016), a few kilome- tres south of the Abisko Scientific Research Station in northernmost Sweden (68°21′ N, 18°49′ E). In the area, the forest and tundra sites are located at altitudes of 523 and 616 m a.s.l. respectively. The mean annual temperature in the period 1960–1990 was −0.2°C (Callaghan et al., 2013) and in the period 1991–2018 was 0.7°C (Abisko Scientific Research Station, meteorological data from the Abisko Observatory; monthly sums for 1991–2018). Mean annual (covering years 1914– 2013) precipitation in the area is 307 mm (Abisko Scientific Research Station). In 2018, the year of the latest vegetation survey and CO2 measurements, the annual mean temperature was −0.2°C and the precipitation 287 mm. The mean temperature of July in 2018 was 15°C, the second warmest in the recorded history and around 1.5°C above the long-term average. In the area, mountain birch (Betula pubescens ssp. czerepanovii, hereafter B. czerepenovii) forms the treeline. The forest understo- rey is dominated by the dwarf shrubs E. hermaphroditum, Vaccinium myrtillus and Vaccinium uliginosum and the grass Deschampsia flexu- osa. Bryophytes are common, but only a few species of lichens are found (Kaarlejärvi et al., 2012). In tundra, the field layer is dominated by the dwarf shrubs E. hermaphroditum, V. vitis-idaea, B. nana, V. uliginosum and the bottom layer by lichens such as Peltigera aph- tosa and Nephroma arcticum and bryophytes such as Ptilidium cili- are (Kaarlejärvi et al., 2012). In the area, the vegetation is affected by background herbivory as well as cyclic population outbreaks of insects and rodents (Callaghan et al., 2013; Kristensen et al., 2020; Olofsson et al., 2013). Several rodent and geometric moth outbreaks have occurred while the current experiment has been ongoing and some of them have caused severe damage to the vegetation. Most notable is the severe geometric moth outbreak by Epirrita autumnata and Operophthera brumara in 2012 (Olofsson et al., 2013). The forest and tundra sites used in this study belong to a net- work of experiments established in mountain birch forests and tun- dra heath patches in 1998 (Sjögersten et al., 2003). By the time of the study site establishment, no permission for the fieldwork was required and the nature of the experiment does not contravene pro- tections currently in place in the area. For this study, we included seven control plots (1 m2) and seven experimentally warmed plots, randomly set out at the mountain birch forest site, and eight control plots and eight warmed plots at the tundra site (Table 1). Warming was induced using International Tundra Experiment (ITEX) hexagonal F I G U R E 1  Location of the study area on the map and photos of two study sites: Tundra and Forest. Photo credit: Minna Männistö. Scale of inference Scale at which the factor of interest is applied Number of replicates at the appropriate scale Plots Square or hexagonal field plot 7 or 8 of each combination, except 4 or 5 for comparison between years TA B L E 1  Replication statement. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense 222  |    MYRSKY et al. OTCs with a maximum basal diameter of 146 cm. Vegetation in OTCs as well as in the control plots only includes species from the field and ground layers and no B. czerepenovii trees grow in the plots. At the early phases of the experiment, OTCs increased July air tem- peratures by 0.8 and 2.5°C in the forest and the tundra sites respec- tively (Sjögersten et al., 2003). No air temperature is available for the time of the present investigation, but the average soil temperature, measured every 2 h at 3 cm depth for 8 June–20 August 2018 using EasyLog EL-USB-1 data loggers (Lascar Electronics, Whiteparish) was 9.1 ± 0.4°C and 9.4 ± 0.2°C in the control and OTC plots in the forest site, and 8.7 ± 0.3 and 8.1 ± 0.5°C in the control and OTC plots in the tundra site respectively. Increased plant abundances insulated the ground to an extent that the OTCs experienced soil tempera- tures close to those in the control plots in the forest, and even lower temperatures than control plots in the tundra site (Stark et al., 2023). 2.2  |  Vegetation analyses The plant community composition was earlier recorded in 1999 and 2009 in five control plots and five OTCs in both habitats (Kaarlejärvi et al., 2012). We used the same plots during the present investiga- tion and analysed the composition of vegetation in July 2018 with the point intercept method: in OTCs, a total of 87 pins was systemat- ically distributed among three diagonals of the hexagons, 29 pins per diagonal (as in Kaarlejärvi in 2012). For each pin, the total number of hits as well as the height of the highest hit were recorded for each plant. Only one hit for each species was counted at the ground layer for each pin. The same method was applied to control plots. Later the total number of hits was normalized to hits per 100 pins. Data from one forest plot were discarded because of poor plot condition. 2.3  |  Ecosystem carbon flux analyses For the ecosystem carbon flux analyses, we included a few additional plots to have seven plots per treatment in the forest and eight plots per treatment in the tundra. The fluxes were analysed at 2-week intervals throughout the growing season 2018 (from 5 June to 19 August) using a closed system composed of a custom-built acrylic chamber (diameter 146 cm, height 60 cm) coupled to a Vaisala Carbon Dioxide Probe GMP343, Vaisala Humidity and Temperature Probe HMP75 and Vaisala Measurement Indicator MI70. Measurements included four consecutive measures of gradually changing light in- tensity: ambient light, 35% and 60% shading, and darkness to reveal ecosystem respiration, Re. Shading was implemented using hoods made of single- and double-layer white mosquito nets while darkness was obtained by covering the chamber with an opaque white hood. The chamber was vented before each measurement and placed carefully on top of the study plot so that the leakage of air from beneath the chamber was minimized (as in Väisänen et  al.,  2014). Photosynthetically active radiation (PAR) within the chamber was recorded using an HD 9021 Quantum-Photo-Radiometer. The CO2 concentration, temperature and humidity within chambers were logged at 5-s intervals for 90 s. The CO2 flux was calculated using CO2 and the chamber microclimate data and corrected for changes in temperature and water vapour pressure (Hooper et al., 2002). The net CO2 flux with light intensity above zero was regarded as NEE. For NEE, negative fluxes indicate a net uptake of CO2 from the at- mosphere, whereas positive fluxes indicate a net release of CO2 into the atmosphere. For the comparison of daily CO2 flux measurements between the treatments and control plots we normalized GPP to the PAR level of 600 μmol m−2 s−1. The GPP was calculated from the NEE and Re as: Daily plot-specific GPP values were fitted to their corresponding PAR levels using the nonlinear least squares (nls) function from stats pack- age in R software environment as (R Core Team, 2021): where i stands for ith plot and j for jth date, Amax is the maximum GPP rate when saturated to light (mg CO2 m−2 h−1) and k is the half-satura- tion light constant (μmol m−2 s−1). Subsequently, the GPP600 was calcu- lated for each plot and day at the light level of 600 μmol m−2 s−1 using Equation (2) with PAR set to 600 (Figure S1). 2.4  |  Statistical analyses To test the effects of OTC, habitat, time (year or date) and their inter- actions on the abundance of total vegetation, plant species and plant functional types, and on NEE600, GPP600 and Re, we used a repeated measures ANOVA. When ANOVA revealed a significant interaction, pairwise t-test was used to test the effect of year, habitat or OTC in the split data. The normality of data distribution was evaluated using the Shapiro–Wilk test and histograms, and the homogene- ity of variances using the Levene's test. To meet the assumptions of ANOVA, B. nana, V. uliginosum and graminoid abundances were square root transformed. The abundance of B. nana and V. myrtil- lus were tested only for warming and year effects, since they grew almost exclusively in one of the habitats only. The abundances of forbs and pteridophytes were not tested individually as they grew only on few plots. To investigate the links between CO2 fluxes and total vascular plant, deciduous and evergreen shrub abundances, residual flux data devoid of habitat effect were first produced using ANOVA with habitat as a predictor. Pearson correlation was then used to test the significance of associations between plant abundances and residual CO2 fluxes. If correlations were calculated with origi- nal data they would be dominated by mean habitat differences and would not reflect the within-habitat association of variables across study plots. Finally, the vascular vegetation from all plots, where the CO2 exchange was measured, was fitted to an ordination diagram with non-metric multidimensional scaling using the vegan package (1)GPP = NEE − Re . (2)GPPij = AmaxPAR k + PAR , 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense     |  223MYRSKY et al. (Oksanen et al., 2022). All statistical analyses were performed and figures drawn with R software version 4.1.0 for statistical computing (R Core Team, 2021). 3  |  RESULTS 3.1  |  Vegetation development We observed a significant treatment × habitat × year interaction ef- fect on total vascular plant abundance (Table 2; Figure 2). A signifi- cant positive effect of OTCs was found in both forest (p = 0.006) and tundra (p = 0.026) in 2009 (with an especially strong effect in forest where vascular plants were 126% more abundant in OTCs than in controls), and a positive OTC effect was found in tundra (p = 0.003) but not in forest (p = 0.101) in 2018 (Figure 2). In control plots, the total abundance of vascular plants did not change significantly from 1999 to 2009 (p = 0.674 in forest and p = 0.123 in tundra) or from 2009 to 2018 (p = 0.704 in forest and p = 0.164 in tundra) (Figure 2). In contrast, the abundance increased by 118% in forest (p = 0.04) and by 103% in tundra (p < 0.001) (Figure 2) in OTCs from 1999 to 2009, with no further statistically significant changes in the abundance from 2009 to 2018 in either forest (p = 0.211) or tundra (p = 0.724). Overall, the vascular plants were 71% and 59% more abundant in OTCs than control plots in years 2009 and 2018 respectively. Throughout the years, a major part of the vegetation consisted of dwarf shrubs, both deciduous and evergreen. The overall per- centage of dwarf shrubs of the total vascular plant abundance was 79%, 95% and 90% in the years 1999, 2009 and 2018 respectively. In both habitats, evergreen dwarf shrubs dominated over decid- uous ones. Evergreen dwarf shrubs were on average 51% more abundant in OTCs than controls in both 2009 (p = 0.006) and 2018 (p = 0.002), but not in 1999 (p = 0.436) (Figure  3). The abundance increased on average by 133% from 1999 to 2009 (p < 0.001), but then decreased by 30% from 2009 to 2018 (p = 0.003). Empetrum hermaphroditum dominated the evergreen dwarf shrubs: on average, it formed 83% of the functional group abundance and 56% of the total vascular vegetation, but the abundance differed between the habitats (Table 2), being 143% higher in tundra than forest (Figure 3). Response of evergreen shrubs to warming was largely driven by E. hermaphroditum (Figure  3). Another evergreen dwarf shrub spe- cies, V. vitis-idaea, formed on average 15% of total abundance (see Figure S2 and Table S1 for V. vitis-idaea results). The deciduous dwarf shrubs made up on average 22% of the total vascular vegetation. Their abundance was on average 72% higher in OTCs than controls (Table  2; Figure  3). There was a significant year × habitat effect (Table 2) because deciduous shrubs were 557% more abundant in forest than in tundra in 1999 (p < 0.001), while no significant habitat effect was observed later in 2009 (p = 0.639) and 2018 (p = 0.584). Accordingly, the abundance of deciduous dwarf shrubs increased by 757% from 1999 to 2018 in tundra (p = 0.005), but not in forest (p = 0.883) (Figure 3). Betula nana grew almost ex- clusively in tundra, where it was on average 131% more abundant TA B LE 2   F- a nd p -s ta tis tic s of A N O VA o f t he e ff ec ts o f t re at m en t ( co nt ro l a nd O TC p lo ts ), ye ar (1 99 9, 2 00 9 an d 20 18 ) a nd h ab ita t ( fo re st a nd tu nd ra ) o n to ta l v as cu la r p la nt , f un ct io na l g ro up an d sp ec ie s ab un da nc es . Tr ea tm en t Ye ar H ab ita t Tr ea tm en t ×  ye ar Tr ea tm en t ×  h ab ita t Ye ar  ×  h ab ita t Tr ea tm en t ×  h ab ita t ×  ye ar F p F p F p F p F p F p F p Va sc ul ar p la nt s 25 .0 27 < 0. 00 1 14 .7 00 < 0. 00 1 21 .1 58 0. 00 1 5. 79 1 0. 01 0 0. 07 4 0. 78 8 2. 22 8 0. 13 3 6. 30 1 0. 00 7 Ev er gr ee ns 11 .0 48 0. 00 3 17 .8 76 < 0. 00 1 54 .1 57 < 0. 00 1 4. 80 9 0. 01 9 0. 91 5 0. 35 0 3. 09 5 0. 06 6 2. 61 7 0. 09 7 E. h er m ap hr od itu m 8. 27 0 0. 00 9 13 .0 68 < 0. 00 1 42 .9 08 < 0. 00 1 3. 40 0 0. 05 3 0. 47 9 0. 49 6 1. 54 6 0. 23 6 2. 22 6 0. 13 3 D ec id uo us 19 .1 56 < 0. 00 1 5. 72 3 0. 01 0 2. 48 9 0. 13 0 3. 36 6 0. 05 4 2. 01 8 0. 17 0 6. 45 5 0. 00 7 2. 57 5 0. 10 0 B. n an aa 6. 26 7 0. 02 1 9. 37 3 0. 00 1 — — 0. 96 8 0. 39 0 — — — — — — G ra m in oi ds a 1. 20 2 2. 28 5 6. 92 1 0. 00 5 25 .7 67 < 0. 00 1 4. 88 9 0. 01 8 0. 42 8 0. 52 0 17 .6 37 < 0. 00 1 0. 24 0 0. 78 9 Br yo ph yt es 4. 23 6 0. 05 2 3. 20 5 0. 06 1 11 .4 82 0. 00 3 2. 49 7 0. 10 6 0. 24 2 0. 62 8 8. 86 6 0. 00 2 0. 36 5 0. 69 9 N ot e: S ig ni fic an ce fo r p ≤ 0 .0 5 is in di ca te d by b ol d. a Sq ua re ro ot tr an sf or m ed . 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense 224  |    MYRSKY et al. in OTCs than controls (Table 2; Figure 3) and where its abundance increased through the years in both treatments, first by 258% from 1999 to 2009 (p = 0.036) and then further by 157% from 2009 to 2018 (p = 0.002) (Figure 3). The two other common deciduous dwarf shrub species were V. myrtillus and V. uliginosum. Of these, V. uligi- nosum was common in both habitats, while V. myrtillus grew almost exclusively in forest (see Figure S2 and Table S1 for V. myrtillus and V. uliginosum results). For graminoids the year × habitat interaction effect (Table  2) was due to graminoids being more abundant in forest than tundra in 1999 (p < 0.001), but not in 2009 (p = 0.397) or 2018 (p = 0.335) (Figure  4). Accordingly, their abundance decreased by 70% from F I G U R E 2  Vascular plant abundance (mean ± SD, n = 4–5) in control and OTC plots from 1999 to 2018 in forest and tundra habitats. F I G U R E 3  Abundance (mean ± SD, n = 4–5) of all evergreen and deciduous dwarf shrubs, Empetrum hermaphroditum and Betula nana in control and OTC plots from 1999 to 2018 in forest and tundra habitats. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense     |  225MYRSKY et al. 1999 to 2009 in forest (p = < 0.001), but not in tundra (p = 0.352) (Figure 4). The year × treatment effect (Table 2) appeared because graminoids were significantly more abundant in controls than in OTCs in 2018 (p = 0.049), but not earlier. For bryophytes, the hab- itat × year interaction effect (Table 2) was because bryophytes were more abundant in forest than tundra in 1999 (p = 0.005) and 2018 (p = 0.018), but not in 2009 (p = 0.2). OTCs decreased bryophyte abundance on average by 24% (Figure 4), but this effect remained marginally significant (Table 2). 3.2  |  CO₂ exchange NEE600, GPP600 and Re varied across the growing season, and for all three variables the habitat and treatment effects depended on the measurement date (Table  3). The three-way treatment × habi- tat × date interaction was statistically significant for NEE600 (Table 3): no significant OTC effect was observed in tundra (p = 0.122), whereas in forest, OTCs increased NEE600 (p < 0.001) from mid-June to the end of July, but not in the first and last measurement (Figure 5a). Tundra acted as a C sink throughout the growing season and NEE600 remained relatively stable from July to August in both treatments (Figure 5a). In forest, NEE600 varied more in time and control plots shifted temporarily from sink to source in late July (Figure 5a). The mean growing season NEE600 was on average 64% more negative in OTCs than controls (Figure 5b; Table 4). Treatment × habitat interaction effect on mean growing season GPP600 (Figure 5d; Table 4) was because OTCs increased GPP600 by 25% in forest (p = 0.070), while no statistically significant effect ap- peared in tundra (p = 0.118). A treatment × habitat × date interaction effect was also observed on GPP600 although the effect was only marginally significant (Table 3; Figure 5c). In tundra, GPP600 was on average 91% higher in control than OTCs in the fourth and sixth measurements (p = 0.036 and p = 0.013 respectively), while in forest, GPP600 was 161% higher in OTCs than controls in the second mea- surement (p = 0.029) (Figure 5c). The mean growing season Re differed between habitats and treatments (Table 4): Re was on average 28% greater in forest than tundra and on average 74% greater in controls than OTCs (Figure 5f). For Re, the treatment × date interaction effect (Table 3) was because Re was on average 48% lower in OTCs than controls in the last three measurements (p = 0.001, p = 0.009, p < 0.001 respectively), but not yet in June (Figure 5e). The habitat × date interaction effect appeared because Re was similar or higher in tundra than forest until mid-July, but lower in tundra than forest in later measurements (Figure 5e). 3.3  |  Correlations between plant group abundances and CO₂ fluxes When the abundance of all vascular plants and evergreen and de- ciduous dwarf shrubs were contrasted with the residuals of CO2 fluxes devoid of the habitat effect, the only statistically significant correlation that appeared was between residual Re and total vascular plant abundance (Figure 6; Table 5). This correlation was negative, as was the marginally significant correlation between evergreen dwarf shrub abundance and residual NEE600 (Figure 6; Table 5). 3.4  |  Multivariate analysis of treatment and habitat effects on vegetation The non-metric multidimensional scale analysis (NMDS) of vascular plant abundances from 2018 clearly separated forest and tundra habitats (Figure 7). Within the habitats, the treatments overlapped but the treatment effect was still clear and OTCs from the two habi- tats resembled each other more than the controls did. Fitted vectors with a significant p-value show that both deciduous and evergreen F I G U R E 4  Abundance (mean ± SD, n = 4–5) of graminoids and bryophytes in control and OTC plots from 1999 to 2018 in forest and tundra habitats. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense 226  |    MYRSKY et al. dwarf shrubs increased towards tundra (r2 = 0.26, p = 0.009 and r2 = 0.18, p = 0.049 respectively), with the deciduous shrubs increas- ing especially towards OTCs. None of the CO2 fluxes gained a statis- tically significant value and appeared in the NMDS diagram. 4  |  DISCUSSION The effects of warming on shrub abundance followed in many parts our predictions and earlier observations, whereas the effects on ecosystem CO2 exchange contrasted our predictions. In particular, warming effects on Re were negative, not positive, which calls for a reconsideration of mechanisms that link warming and CO2 fluxes in northern ecosystems. We predicted that after 20 years, (H1) warm- ing should have a stronger positive effect on dwarf shrub growth in tundra than forest because plant growth may be more sensitive to changes in temperature in tundra (Sundqvist et al., 2020), and an earlier survey had found increased abundance of B. nana (Kaarlejärvi et al., 2012). In line with our first hypothesis, OTCs had significantly higher total abundance of vascular plants than control plots in tun- dra, but not in forest in 2018 (Figure 2), and as expected, B. nana was driving this difference as other shrubs responded to warming similarly in both habitats (Figure  3). In particular, E. hermaphrodi- tum abundance was higher in OTCs than controls in both habitats (Figure 3). Our prediction that (H2) OTCs will increase GPP and Re in both habitats was clearly refuted: warming decreased Re in both habitats and increased GPP only in forest (Figure 5f). These results also refuted our prediction that warming should increase Re more in tundra than forest due to the positive effect of B. nana on carbon cycling, and that (H3) the greater increase in Re under warming in tundra should outweigh the increase in GPP and result in reduced ecosystem C sink strength (i.e. a less negative NEE). In contrast, ecosystem C sink strength was significantly stronger in OTCs than controls in both habitats (Figure 5b). For forest, we predicted (H3) a neutral or slightly positive warming effect on C sink strength as warming was earlier shown to increase the abundances of evergreen dwarf shrubs in forest (Kaarlejärvi et al., 2012). This prediction was supported in terms that the positive effect of warming on C sink strength was stronger in forest than tundra (Figure 5a,b). 4.1  |  Vegetation development Our results add to the cumulating evidence that woody shrubs are driving the Arctic greening under climate warming through both in- creasing production and expanding distribution (Elmendorf, Henry, Hollister, Björk, Bjorkman, et al., 2012; Elmendorf, Henry, Hollister, Björk, Boulanger-Lapointe, et al., 2012): evergreen shrubs were on average 51% more abundant in OTCs than controls in the 2009 and 2018 surveys (Figure 3) and deciduous shrubs on average 72% more abundant in OTCs than controls (Figure 3). Shrub expansion has been associated with increasing abundance of deciduous spe- cies (e.g. Elmendorf, Henry, Hollister, Björk, Bjorkman, et al., 2012; TA B LE 3   F- a nd p -s ta tis tic s of A N O VA o f t he e ff ec ts o f t re at m en t ( co nt ro l a nd O TC p lo ts ), ha bi ta t ( fo re st a nd tu nd ra ) a nd d at e (ti m e of th e gr ow in g se as on ) o n ec os ys te m N EE 60 0, G PP 60 0 an d R e. Tr ea tm en t H ab ita t D at e Tr ea tm en t ×  h ab ita t Tr ea tm en t ×  d at e H ab ita t ×  d at e Tr ea tm en t ×  h ab ita t ×  d at e F p F p F p F p F p F p F p N EE ₆₀ ₀ 13 .5 09 < 0. 00 1 0. 36 2 0. 54 9 5. 28 4 < 0. 00 1 5. 09 9 0. 02 7 2. 20 4 0. 06 30 7. 75 4 < 0. 00 1 2. 37 6 0. 04 7 G PP ₆₀ ₀ 0. 00 6 0. 93 9 1. 05 9 0. 30 7 4. 93 6 < 0. 00 1 6. 40 8 0. 01 4 3. 43 7 0. 00 8 2. 75 9 0. 02 5 2. 09 9 0. 07 5 Rₑ 33 .2 81 < 0. 00 1 9. 13 5 0. 00 3 35 .0 65 < 0. 00 1 1. 90 4 0. 17 2 5. 41 1 < 0. 00 1 4. 31 2 0. 00 2 0. 64 7 0. 66 5 N ot e: S ig ni fic an ce fo r p ≤ 0 .0 5 is in di ca te d by b ol d. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense     |  227MYRSKY et al. Elmendorf, Henry, Hollister, Björk, Boulanger-Lapointe, et al., 2012; Myers-Smith et al., 2011), but recent studies have also found strong increases in the abundance of evergreen species (Stark et al., 2021; Vuorinen et al., 2017). In our study, both evergreen and deciduous shrubs increased under warming, but the greater percentual in- crease in deciduous than evergreen shrubs suggest that, also in our study area, the deciduous shrubs benefit of warming relatively more than the evergreen shrubs. Supporting the idea that shrub expansion often occurs at the expense of graminoids (e.g. Myers-Smith et  al.,  2011) and bryo- phytes (Alatalo et al., 2020), we found that the abundances of gram- inoids and bryophytes were lower in OTCs than controls (Figure 4). Considering that graminoids are the plant group in Arctic environ- ments whose growth most closely follows soil nutrient availability (Croll et  al.,  2005), the negative warming effect may be mediated by the decelerated soil nutrient mineralization, supposedly following the accumulation of E. hermaphroditum litter (Bråthen et al., 2017; Vowles & Björk, 2018). Bryophytes in turn likely decrease when the dominant shrubs become so tall that they limit light availability at the ground layer (Alatalo et al., 2020). Although the OTC effect on plant groups in our study largely followed our predictions, changes along the years show that the OTC effect was not cumulative. Most notably, only B. nana of all shrub species had higher abundance in OTCs in 2018 than in 2009. The likely reason for this finding is a series of herbivore outbreaks that took place in the study area between 2009 and 2018. The geometric moth outbreaks by E. autumnata and O. bru- mata in 2012 and 2013, a lemming outbreak in 2011 (Callaghan et  al.,  2013; Olofsson et  al.,  2013) and another moth outbreak in 2017 (S. Sjögersten, personal observation) had drastic effects F I G U R E 5  Measurement (a, c, e) and growing season (b, d, f) means (± SD) of NEE600, GPP600 and Re in forest and tundra habitats in 2018. All values are expressed as mg CO2 m−2 h−1. Measurement 1 = 6–7 June, 2 = 17–19 June, 3 = 30 June–1 July, 4 = 13–16 July, 5 = 27–28 July and 6 = 18–19 August. Forest Tundra 2 4 6 2 4 6 −900 −600 −300 0 300 N EE 60 0 m g C O 2 m −2 h −1 Control OTC(a) Forest Tundra −400 −200 0 Control OTC(b) Forest Tundra 2 4 6 2 4 6 −1000 −750 −500 −250 0 G PP 60 0 m g C O 2 m −2 h −1 (c) Forest Tundra −600 −400 −200 0 (d) Forest Tundra 2 4 6 2 4 6 0 250 500 MeasurementR e m g C O 2 m −2 h −1 (e) Forest Tundra 0 100 200 300 400 (f) Treatment Habitat Treatment × habitat F p F p F p NEE₆₀₀ 10.323 0.003 0.626 0.436 3.454 0.074 GPP₆₀₀ 0.005 0.946 0.375 0.546 5.417 0.028 Rₑ 26.702 <0.001 5.823 0.023 1.341 0.257 Note: Significance for p ≤ 0.05 is indicated by bold. TA B L E 4  F- and p-statistics of ANOVA of the effects of treatment (control and OTC plots) and habitat (forest and tundra) on growing season averages of ecosystem NEE600, GPP600 and Re. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense 228  |    MYRSKY et al. on the mountain birches and the ground vegetation, and a study at nearby plots reported 34%–76% reductions of biomass for dwarf shrub species after the outbreaks (Olofsson et  al.,  2013). Microtine rodents can significantly reduce the abundance dwarf shrubs, including the lower palatability species like E. hermaphro- ditum (Olofsson et al., 2014; Tuomi et al., 2019), as well as bryo- phytes (Olofsson et  al.,  2014). The rodent outbreak thus likely explains the decline in bryophytes in tundra between 2009 and 2018. Both moth and rodent outbreaks have regular cycles in sub- arctic Fennoscandia and can be seen as natural drivers of vege- tation (Callaghan et  al.,  2013; Olofsson et  al.,  2013, 2014). Our results from a long-term experiment provide a novel perspective into how strongly these drivers can control vegetation trajectories under warming Arctic. Assuming that herbivore outbreaks caused comparable damage to all shrub species in our study plots, B. nana appears as an excep- tion in two ways. First, the positive OTC effect seemed to accumu- late in B. nana abundance over years despite herbivory. Second, B. nana abundance also increased in tundra controls over years, thus suggesting that, despite herbivory outbreaks, B. nana could also ben- efit of the naturally increasing temperatures over the study period (Rantanen et al., 2022). These results are in line with the high plas- ticity of B. nana growth strategy (Bret-Harte et al., 2001). Altogether our results support the framework that herbivory outbreaks may nullify the cumulation of positive warming effects on plant growth (Kaarlejäarvi et al., 2015), and further highlight that herbivory may promote species such as B. nana that both significantly benefit of warming and can quickly recover from damage. 4.2  |  Growing season CO2 exchange While the vegetation mostly responded to warming as predicted, this was not the case for CO2 exchange: instead, in contrast to several earlier warming experiments, which report increasing Re under warming (Biasi et al., 2008; Oberbauer et al., 2007; Väisänen et al., 2014; Ylänne et al., 2015), we observed a clear decrease in Re in both forest and tundra. This finding not only opposed our predic- tion, but also showed that the consequence of warming on Re did not depend on whether deciduous or evergreen dwarf shrubs in- crease in abundance under warming, as suggested by Vowles and Björk  (2018). We further predicted that warming will increase the F I G U R E 6  Association of the abundance of all vascular plants and evergreen and deciduous dwarf shrubs with CO2 flux residuals devoid of habitat effect. Residual NEE₆₀₀ Residual GPP₆₀₀ Residual Re r p r p r p Vascular plants −0.25 0.180 0.03 0.890 −0.39 0.032 Evergreens −0.34 0.068 −0.12 0.520 −0.31 0.099 Deciduous −0.01 0.950 0.14 0.450 −0.22 0.250 Note: Significance for p ≤ 0.05 is indicated by bold. TA B L E 5  Correlation coefficients and p-values of Person correlation between plant groups and CO2 flux residuals devoid of habitat effect. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense     |  229MYRSKY et al. GPP in both habitats, but found a marginal positive effect only in forest. The effect of warming on the GPP was strongest during early growing season, a time period where light availability for the understorey vegetation is not yet limited by the tree leaves (Ťupek et  al.,  2008). Studies in the Arctic have indicated that the effect of warming on Re rather than on GPP drives the pattern in NEE (Cahoon et al., 2012; Väisänen et al., 2014). In our study, both lower Re and higher GPP contributed to increased C sink strength in forest, whereas in tundra, where no warming effect on GPP was detected, the lower Re was the sole driver. Our findings that GPP was not affected and Re decreased under warming in tundra contrast earlier studies, which have suggested that increase in deciduous dwarf shrubs accelerates carbon cycling (Parker et al., 2021; Vowles & Björk, 2018) and that B. nana has a key role in this process even in sites where it accounts for a minor pro- portion of the total above-ground vascular plant biomass (Cahoon et al., 2016; Metcalfe & Olofsson, 2015). Although our result was un- expected, a number of studies have also found CO2 fluxes in tundra to be relatively unresponsive to changes in plant community compo- sition (Sundqvist et al., 2020; Ylänne et al., 2015). It therefore seems likely that under some conditions, other mechanisms outweigh the importance of plant species on ecosystem CO2 exchange. Noteworthy, we found a significant negative correlation be- tween Re and the total vascular plant abundance, indicating that the increase in total amount of vascular plants rather than a change in a single plant functional group was responsible for the lower Re in OTCs. Previous studies have demonstrated that taller and denser vegetation cover buffers the soil against warmer air temperatures (e.g. Way & Lapalme, 2021; Weintraub & Schimel, 2005), and that the effect of shrubs on CO2 exchange in tundra ecosystems depends on soil temperatures (Cahoon et al., 2012). Increased insulation of soil through increased vascular plant cover could thus significantly contribute to the lower Re under warming. Why this outcome has not been observed in other studies may derive from the longevity of our experiment: it may take decades until shrubs increase to the extent that their cooling impact on the growing season soil microcli- mate outweighs the direct effects of plant groups on CO2 exchange. However, not even in our case can the indirect effects through soil microclimate solely explain decreasing Re as we detected lower mean soil temperature and daily temperature maximums inside the OTCs in tundra, but not in the forest (Stark et al., 2023). In the forest, shading by B. czerepanovii trees likely diminishes the role of ground vegetation in soil microclimate. Results from the same experiment also suggest that nutrient competition between plants and soil mi- croorganisms increase under warming, thus negating the effects of warming on microbial activity in C decomposition (Stark et al., 2023). It appears that several different mechanisms could contribute to de- creasing Re. The growing season 2018 was especially hot and dry during the mid and late growing season. This led to a decrease in GPP and Re in other parts of Fennoscandia (Silfver et al., 2020), but we found an increasing trend of GPP during the same phase of the growing season. We also found the highest Re rates in mid and late July, sug- gesting that Re followed the same pattern as temperature and was not negatively affected by high temperatures. These contrasts likely arise from differences in climate and other environmental factors between the research sites. The Silfver et  al.  (2020) study site in northern Finland has subcontinental climate, acidic bedrock and lower vegetation cover while our site has a more maritime climate with higher precipitation, more nutrient-rich bedrock and a thicker moss and dwarf shrub cover, which all likely buffer against severe moisture stress. F I G U R E 7  Two-dimensional non- metric multidimensional scale (NMDS) ordination of plant species composition within different treatments in 2018 (stress = 0.17). The ordination shows the vectors of species and functional groups with a p-value less than 0.05. Ellipses represent 95% confidence intervals for treatment groups. 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense 230  |    MYRSKY et al. 4.3  |  Implications Our findings of higher shrub abundances after 20 years of warming agree with the notion of Arctic shrubification (Elmendorf, Henry, Hollister, Björk, Bjorkman, et al., 2012; Elmendorf, Henry, Hollister, Björk, Boulanger-Lapointe, et  al.,  2012; Myers-Smith et  al.,  2020). However, owing to several population outbreaks of geometric moths and rodents during the experiment (Callaghan et al., 2013; Olofsson et al., 2013), the effect of warming on shrubs was not cumulative over time. The only exception was B. nana, which continued to in- crease in abundance over time in both controls and OTCs. Given that moth outbreaks will likely increase in frequency under warming climate (Jepsen et al., 2013), we suggest that climate warming may promote species that can combine efficient compensatory growth with a capacity of benefiting from improved growing conditions. Furthermore, in the long term, the indirect impact of increased total vascular plant abundance on soil microclimate may become increas- ingly important for ecosystem CO2 exchange. Given the slow re- sponses of dwarf shrub vegetation to changing temperatures and the cyclicity of herbivore outbreaks in the Arctic, only long-term field studies may be able to reveal the impact of climate warming on the dynamics of vegetation and the ecosystem C sink strength. AUTHOR CONTRIBUTIONS Sofie Sjögersten and Johan Olofsson established and maintained the field experiment. Sari Stark and Minna K. Männistö planned the pre- sent study with input from Johan Olofsson. Eero Myrsky performed new data collection. Eero Myrsky conducted statistical analyses to- gether with Juha Mikola. Elina Kaarlejärvi, Johan Olofsson and Sofie Sjögersten contributed with earlier data. Boris Tupek processed the light standardization of the CO2 exchange data. Eero Myrsky led the writing of the paper to which all co-authors contributed with discus- sion and text and gave the final approval of the paper. ACKNO​WLE​DG E​MENTS We thank Jonas Gustafsson for helping with vegetation analyses, and Manoj Kumar and Juho Haveri-Heikkilä for helping with CO2 ex- change measurements and data processing. This project was funded by the Academy of Finland (decision number 310776 to Minna K. Männistö). Eero Myrsky also worked with a personal grant from Finnish Cultural Foundations Lapland regional fund (decision num- ber 40221879). CONFLIC T OF INTERE S T S TATEMENT The authors state that they have no conflict of interest. DATA AVAIL ABILIT Y S TATEMENT Data available from the Dryad Digital Repository (Myrsky et  al., 2023): https://​doi.​org/​10.​5061/​dryad.​612jm​649d. ORCID Eero Myrsky  https://orcid.org/0000-0001-6128-2580 Juha Mikola  https://orcid.org/0000-0002-4336-2648 Elina Kaarlejärvi  https://orcid.org/0000-0003-0014-0073 Johan Olofsson  https://orcid.org/0000-0002-6943-1218 Sofie Sjögersten  https://orcid.org/0000-0003-4493-1790 Boris Tupek  https://orcid.org/0000-0003-3466-0237 Minna K. Männistö  https://orcid.org/0000-0001-9390-1104 Sari Stark  https://orcid.org/0000-0003-4845-6536 R E FE R E N C E S Alatalo, J. M., Jägerbrand, A. K., Erfanian, M. B., Chen, S., Sun, S.-Q., & Molau, U. (2020). Bryophyte cover and richness decline after 18 years of experimental warming in alpine Sweden. AoB Plants, 12, plaa061. https://​doi.​org/​10.​1093/​aobpla/​plaa061 Biasi, C., Meyer, H., Rusalimova, O., Hämmerle, R., Kaiser, C., Baranyi, C., Daims, H., Lashchinsky, N., Barsukov, P., & Richter, A. (2008). Initial effects of experimental warming on carbon exchange rates, plant growth and microbial dynamics of a lichen-rich dwarf shrub tundra in Siberia. Plant and Soil, 307, 191–205. Bjorkman, A. D., García Criado, M., Myers-Smith, I. H., Ravolainen, V., Jónsdóttir, I. S., Westergaard, K. B., Lawler, J. P., Aronsson, M., Bennett, B., Gardfjell, H., Heiðmarsson, S., Stewart, L., & Normand, S. (2020). Status and trends in Arctic vegetation: Evidence from ex- perimental warming and long-term monitoring. Ambio, 49(3), 678– 692. https://​doi.​org/​10.​1007/​s13280-​019-​01161-​6 Blok, D., Heijmans, M. P. D., Schaepman-Strub, G., Kononov, A. V., Maximov, T. C., & Berendse, F. (2010). Shrub expansion may reduce summer permafrost thaw in Siberian tundra. Global Change Biology, 16, 1296–1305. Blok, D., Sass-Klaassen, U., Schaepman-Strub, G., Heijmans, M. M. P. D., Sauren, P., & Berendse, F. (2011). What are the main climate drivers for shrub growth in northeastern Siberian tundra? Biogeosciences, 8, 1169–1179. Bradford, M. A., Wieder, W. R., Bonan, G. B., Fierer, N., Raymond, P. A., & Crowther, T. W. (2016). Managing uncertainty in soil carbon feed- backs to climate change. Nature Climate Change, 6, 751–757. Bråthen, K. A., Gonzales, V. T., & Yoccoz, N. G. (2017). Gatekeepers to the effects of climate warming? Niche construction restricts plant community changes along a temperature gradient. Perspectives in Plant Ecology, Evolution and Systematics, 30, 71–81. Bret-Harte, M. S., Shaver, G. R., Zoerner, J. P., Johnstone, J. F., Wagner, J. L., Chavez, A. S., Gunkelman, R. F. I. V., Lippert, S. C., & Laundre, J. A. (2001). Developmental plasticity allows Betula nana to dominate tundra subjected to an altered environment. Ecology, 82, 18–32. Cahoon, S. M. P., Sullivan, P. F., & Post, E. (2016). Greater abundance of Betula nana and early onset of the growing season increase eco- system CO2 uptake in West Greenland. Ecosystems, 19, 1149–1163. Cahoon, S. M. P., Sullivan, P. F., Shaver, G. R., Welker, J. M., & Post, E. (2012). Interactions among shrub cover and the soil microclimate may deter- mine future Arctic carbon budgets. Ecology Letters, 15, 1415–1422. Callaghan, T. V., Jonasson, C., Thierfelder, T., Yang, Z., Hedenås, H., Johansson, M., Molau, U., Van Bogaert, R., Michelsen, A., Olofsson, J., Gwynn-Jones, D., Bokhorst, S., Phoenix, G., Bjerke, J. W., Tømmervik, H., Christensen, T. R., Hanna, E., Koller, E. K., & Sloan, V. L. (2013). Ecosystem change and stability over multiple decades in the Swedish subarctic: Complex processes and multiple drivers. Philosophical Transactions of the Royal Society B: Biological Sciences, 368, 20120488. Croll, D. A., Maron, J. L., Estes, J. A., Danner, E. M., & Byrd, G. V. (2005). Introduced predators transform subarctic islands from grassland to tundra. Science, 307, 1959–1961. DeMarco, J., Mack, M. C., & Bret-Harte, M. S. (2011). The effects of snow, soil microenvironment, and soil organic matter quality on N availability in three Alaskan arctic plant communities. Ecosystems, 14, 804–817. https://​doi.​org/​10.​1007/​s10021-​011-​9447-​5 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense https://doi.org/10.5061/dryad.612jm649d https://orcid.org/0000-0001-6128-2580 https://orcid.org/0000-0001-6128-2580 https://orcid.org/0000-0002-4336-2648 https://orcid.org/0000-0002-4336-2648 https://orcid.org/0000-0003-0014-0073 https://orcid.org/0000-0003-0014-0073 https://orcid.org/0000-0002-6943-1218 https://orcid.org/0000-0002-6943-1218 https://orcid.org/0000-0003-4493-1790 https://orcid.org/0000-0003-4493-1790 https://orcid.org/0000-0003-3466-0237 https://orcid.org/0000-0003-3466-0237 https://orcid.org/0000-0001-9390-1104 https://orcid.org/0000-0001-9390-1104 https://orcid.org/0000-0003-4845-6536 https://orcid.org/0000-0003-4845-6536 https://doi.org/10.1093/aobpla/plaa061 https://doi.org/10.1007/s13280-019-01161-6 https://doi.org/10.1007/s10021-011-9447-5     |  231MYRSKY et al. Elmendorf, S., Henry, G. H. R., Hollister, R. D., Björk, R. G., Boulanger- Lapointe, N., Cooper, E. J., Cornelissen, J. H. C., Day, T. A., Dorrepaal, E., Elumeeva, T. G., Gill, M., Gould, W. A., Harte, J., Hik, D. S., Hofgaard, A., Johnson, D. R., Johnstone, J. F., Jónsdóttir, I. S., Jorgenson, J. C., … Wipf, S. (2012). Plot-scale evidence of tun- dra vegetation change and links to recent summer warming. Nature Climate Change, 2, 453–457. Elmendorf, S. C., Henry, G. H. R., Hollister, R. D., Björk, R. G., Bjorkman, A. D., Callaghan, T. V., Collier, L. S., Cooper, E. J., Cornelissen, J. H. C., Day, T. A., Fosaa, A. M., Gould, W. A., Grétarsdóttir, J., Harte, J., Hermanutz, L., Hik, D. S., Hofgaard, A., Jarrad, F., Jónsdóttir, I. S., … Wookey, P. A. (2012). Global assessment of experimental climate warming on tundra vegetation: Heterogeneity over space and time. Ecology Letters, 15, 164–175. Hobbie, S. E., & Chapin, F. S., III. (1998). The response of tundra plant biomass, aboveground production, nitrogen, and CO2 flux to exper- imental warming. Ecology, 79, 1526–1544. Hooper, D. U., Cardon, Z. G., Chapin, F., & Durant, M. (2002). Corrected calculations for soil and ecosystem measurements of CO2 flux using the LI-COR 6200 portable photosynthesis system. Oecologia, 132, 1–11. Jepsen, J. U., Biuw, M., Ims, R. A., Kapari, L., Schott, T., Vindstad, O. P. L., & Hagen, S. B. (2013). Ecosystem impacts of a range expand- ing forest defoliator at the forest-tundra ecotone. Ecosystems, 16, 561–575. Kaarlejärvi, E., Baxter, R., Hofgaard, A., Hytteborn, H., Khitun, O., Molau, U., Sjögersten, S., Wookey, P. A., & Olofsson, J. (2012). Effects of warming on shrub abundance and chemistry drive ecosystem-level changes in a forest-tundra ecotone. Ecosystems, 15, 1219–1233. Kaarlejäarvi, E., Hoset, K. S., & Olofsson, J. (2015). Mammalian herbi- vores confer resilience of Arctic shrub-dominated ecosystems to changing climate. Global Change Biology, 21, 3379–3388. Kristensen, J. A., Michelsen, A., & Metcalfe, D. B. (2020). Background insect herbivory increases with local elevation but makes minor contribution to element cycling along natural gradients in the sub- arctic. Ecology and Evolution, 10, 11684–11698. Leffler, A. J., & Welker, J. M. (2013). Long-term increases in snow pack elevate leaf N and photosynthesis in Salix arctica: Responses to a snow fence experiment in the High Arctic of NW Greenland. Environmental Research Letters, 8(2), 025023. Mekonnen, Z. A., Riley, W. J., Berner, L. T., Bouskill, N. J., Torn, M. S., Iwahana, G., Breen, A. L., Myers-Smith, I. H., Criado, M. G., Liu, Y., Euskirchen, E. S., Goetz, S. J., Mack, M. C., & Grant, R. F. (2021). Arctic tundra shrubification: A review of mechanisms and impacts on eco- system carbon balance. Environmental Research Letters, 16, 053001. Metcalfe, D. B., & Olofsson, J. (2015). Distinct impacts of different mammalian herbivore assemblages on arctic tundra CO2 exchange during the peak of the growing season. Oikos, 124, 1632–1638. Myers-Smith, I. H., Forbes, B. C., Wilmking, M., Hallinger, M., Lantz, T., Blok, D., Tape, K. D., Macias Fauria, M., Sass-Klaassen, U., Lévesque, E., Boudreau, S., Ropars, P., Hermanutz, L., Trant, A., Collier, L. S., Weijers, S., Rozema, J., Rayback, S. A., Schmidt, N. M., … Hik, D. S. (2011). Shrub expansion in tundra ecosystems: Dynamics, impacts and research priorities. Environmental Research Letters, 6, 045509. Myers-Smith, I. H., Kerby, J. T., Phoenix, G. K., Bjerke, J. W., Epstein, H. E., Assmann, J. J., John, C., Andreu-Hayles, L., Angers-Blondin, S., Beck, P. S. A., Berner, L. T., Bhatt, U. S., Bjorkman, A. D., Blok, D., Bryn, A., Christiansen, C. T., Cornelissen, J. H. C., Cunliffe, A. M., Elmendorf, S. C., … Wipf, S. (2020). Complexity revealed in the greening of the Arctic. Nature Climate Change, 10, 106–117. Myrsky, E., Mikola, J., Kaarlejärvi, E., Olofsson, J., Sjögersten, S., Tupek, B., Männistö, M., & Stark, S. (2023). Higher vascular plant abun- dance associated with decreased ecosystem respiration after 20 years of warming in the forest-tundra -ecotone [dataset]. Dryad Digital Repository, https://​doi.​org/​10.​5061/​dryad.​612jm​649d Nobrega, S., & Grogan, P. (2007). Deeper snow enhances winter respira- tion from both plant-associated and bulk soil carbon pools in birch hummock tundra. Ecosystems, 10, 419–431. Oberbauer, S. F., Tweedie, C. E., Welker, J. M., Fahnestock, J. T., Henry, G. H. R., Webber, P. J., Hollister, R. D., Walker, M. D., Kuchy, A., Elmore, E., & Starr, G. (2007). Tundra CO2 fluxes in response to experimen- tal warming across latitudinal and moisture gradients. Ecological Monographs, 77, 221–238. Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., Simpson, G. L., Solymos, P., Henry, M., Stevens, H., & Wagner, H. (2022). Vegan: Community ecology package. http://​cran.​r-​proje​ct.​ org/​packa​ge=​vegan​ Olofsson, J., Oksanen, L., Oksanen, T., Tuomi, M., Hoset, K. S., Virtanen, R., & Kyrö, K. (2014). Long-term experiments reveal strong interac- tions between lemmings and plants in the Fennoscandian highland tundra. Ecosystems, 17, 606–615. https://​doi.​org/​10.​1007/​s1002​ 1-​013-​9740-​6 Olofsson, J., te Beest, M., & Ericson, L. (2013). Complex biotic interac- tions drive long-term vegetation dynamics in a subarctic ecosys- tem. Philosophical Transactions of the Royal Society, B: Biological Sciences, 368, 20120486. Parker, T. C., Thurston, A. M., Raundrup, K., Subke, J.-A., Wookey, P. A., & Hartley, I. P. (2021). Shrub expansion in the Arctic may induce large-scale carbon losses due to changes in plant–soil interactions. Plant and Soil, 463, 643–665. R Core Team. (2021). R: A language and environment for statistical comput- ing. R Foundation for Statistical Computing. Rantanen, M., Karpechko, A. Y., Lipponen, A., Nordling, K., Hyvärinen, O., Ruosteenoja, K., Vihma, T., & Laaksonen, A. (2022). The Arctic has warmed nearly four times faster than the globe since 1979. Communications Earth & Environment, 3, 168. https://​doi.​org/​10.​ 1038/​s43247-​022-​00498-​3 Shaver, G. R., Street, L. E., Rastetter, E. B., Van Wijk, M. T., & Williams, M. (2007). Functional convergence in regulation of net CO2 flux in heterogeneous tundra landscapes in Alaska and Sweden. Journal of Ecology, 95(4), 802–817. http://www.jstor.org/stable/4496035 Silfver, T., Heiskanen, L., Aurela, M., Myller, K., Karhu, K., Meyer, N., Tuovinen, J.-P., Oksanen, E., Rousi, M., & Mikola, J. (2020). Insect herbivory dampens subarctic birch forest C sink response to warm- ing. Nature Communications, 11, 2529. Sjögersten, S., Turner, B. L., Mahieu, N., Condron, L. M., & Wookey, P. A. (2003). Soil organic matter biochemistry and potential suscep- tibility to climatic change across the forest-tundra ecotone in the Fennoscandian mountains. Global Change Biology, 9, 759–772. Stark, S., Kumar, M., Myrsky, E., Vuorinen, J., Kantola, A., Telkki, V.-V., Sjögersten, S., Olofsson, J., & Männistö, M. (2023). Decreased soil microbial nitrogen under vegetation ‘Shrubification’ in the subarc- tic forest–tundra ecotone: The potential role of increasing nutrient competition between plants and soil microorganisms. Ecosystems, 1–20. https://​doi.​org/​10.​1007/​s1002​1-​023-​00847​-​z Stark, S., Ylänne, H., & Kumpula, J. (2021). Recent changes in mountain birch forest structure and understory vegetation depend on the seasonal timing of reindeer grazing. Journal of Applied Ecology, 58, 941–952. Sundqvist, M. K., Sanders, N. J., Dorrepaal, E., Lindén, E., Metcalfe, D. B., Newman, G. S., Olofsson, J., Wardle, D. A., & Classen, A. T. (2020). Responses of tundra plant community carbon flux to experimen- tal warming, dominant species removal and elevation. Functional Ecology, 34, 1497–1506. Tuomi, M., Stark, S., Hoset, K. S., Väisänen, M., Oksanen, L., Murguzur, F. J. A., Tuomisto, H., Dahlgren, J., & Bråthen, K. A. (2019). Herbivore effects on ecosystem process rates in a low-productive system. Ecosystems, 22, 827–843. Ťupek, B., Minkkinen, K., Kolari, P., Starr, M., Chan, T., Alm, J., Vesala, T., Laine, J., & Nikinmaa, E. (2008). Forest floor versus ecosystem 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense https://doi.org/10.5061/dryad.612jm649d http://cran.r-project.org/package=vegan http://cran.r-project.org/package=vegan https://doi.org/10.1007/s10021-013-9740-6 https://doi.org/10.1007/s10021-013-9740-6 https://doi.org/10.1038/s43247-022-00498-3 https://doi.org/10.1038/s43247-022-00498-3 http://www.jstor.org/stable/4496035 https://doi.org/10.1007/s10021-023-00847-z 232  |    MYRSKY et al. CO2 exchange along boreal ecotone between upland forest and lowland mire. Tellus B: Chemical and Physical Meteorology, 60(2), 153–166. Väisänen, M., Ylänne, H., Kaarlejärvi, E., Sjögersten, S., Olofsson, J., Crout, N., & Stark, S. (2014). Consequences of warming on tundra carbon balance determined by reindeer grazing history. Nature Climate Change, 4, 384–388. Virtanen, R., Oksanen, L., Oksanen, T., Cohen, J., Forbes, B. C., Johansen, B., Käyhkö, J., Olofsson, J., Pulliainen, J., & Tommervik, H. (2016). Where do the treeless tundra areas of northern highlands fit in the global biome system: Toward an ecologically natural subvision of the tundra biome. Ecology and Evolution, 6, 143–158. Vowles, T., & Björk, R. G. (2018). Implications of evergreen shrub expan- sion in the Arctic. Journal of Ecology, 107, 650–655. Vuorinen, K. E. M., Oksanen, L., Oksanen, T., Pyykönen, A., Olofsson, J., & Virtanen, R. (2017). Open tundra persist, but arctic features de- cline—Vegetation changes in the warming Fennoscandian tundra. Global Change Biology, 23, 3794–3807. Way, R. G., & Lapalme, C. M. (2021). Does tall vegetation warm or cool the ground surface? Constraining the ground thermal impacts of upright vegetation in northern environments. Environmental Research Letters, 16, 054077. Weintraub, M. N., & Schimel, J. P. (2005). Nitrogen cycling and the spread of shrubs control changes in the carbon balance of arctic tundra ecosystems. BioScience, 55, 408–415. Wookey, P. A., Aerts, R., Bardgett, R. D., Baptist, F., Bråthen, K. A., Cornelissen, J. H. C., Gough, L., Hartley, I. P., Hopkins, D. W., Lavorel, S., & Shaver, G. R. (2009). Ecosystem feedbacks and cas- cade processes: Understanding their role in the responses of Arctic and alpine ecosystems to environmental change. Global Change Biology, 15, 1153–1172. Ylänne, H., Stark, S., & Tolvanen, A. (2015). Vegetation shift from decid- uous to evergreen dwarf shrubs in response to selective herbivory offsets carbon losses: Evidence from 19 years of warming and sim- ulated herbivory in the subarctic tundra. Global Change Biology, 21, 3696–3711. Zamin, T. J., & Grogan, P. (2012). Birch shrub growth in the low Arctic: The relative importance of experimental warming, enhanced nu- trient availability, snow depth and caribou exclusion. Environmental Research Letters, 7, 034027. SUPPORTING INFORMATION Additional supporting information can be found online in the Supporting Information section at the end of this article. Figure S1: Comparison between the observed and modelled GPP (mg CO2 m−2 h−1) and NEE (mg CO2 m−2 h−1) for control (a, d), tundra (b, e), and forest ground vegetation (c, f) study plots. Figure S2: The significant year × habitat interaction effect on V. vitis-idaea abundance (Table S1) was because V. vitis-idaea was 231% more abundant in tundra than forest in 1999 (p = 0.021) and 59 % more abundant in forest than tundra in 2018 (p = 0.021), while no significant difference was found in 2009 (p = 0.516). Table S1: F-and p-statistics of ANOVA of the effects of treatment (control and OTC plots), year (1999, 2009 and 2018) and habitat (forest and tundra) on total vascular plant, functional group and species abundances. How to cite this article: Myrsky, E., Mikola, J., Kaarlejärvi, E., Olofsson, J., Sjögersten, S., Tupek, B., Männistö, M. K., & Stark, S. (2024). Higher vascular plant abundance associated with decreased ecosystem respiration after 20 years of warming in the forest–tundra ecotone. Functional Ecology, 38, 219–232. https://doi.org/10.1111/1365-2435.14466 13652435, 2024, 1, D ow nloaded from https://besjournals.onlinelibrary.w iley.com /doi/10.1111/1365-2435.14466 by L uonnonvarakeskus, W iley O nline L ibrary on [08/07/2024]. See the T erm s and C onditions (https://onlinelibrary.w iley.com /term s-and-conditions) on W iley O nline L ibrary for rules of use; O A articles are governed by the applicable C reative C om m ons L icense https://doi.org/10.1111/1365-2435.14466 Higher vascular plant abundance associated with decreased ecosystem respiration after 20 years of warming in the forest–tundra ecotone Abstract 1|INTRODUCTION 2|MATERIALS AND METHODS 2.1|Study area and experimental design 2.2|Vegetation analyses 2.3|Ecosystem carbon flux analyses 2.4|Statistical analyses 3|RESULTS 3.1|Vegetation development 3.2|CO₂ exchange 3.3|Correlations between plant group abundances and CO₂ fluxes 3.4|Multivariate analysis of treatment and habitat effects on vegetation 4|DISCUSSION 4.1|Vegetation development 4.2|Growing season CO2 exchange 4.3|Implications AUTHOR CONTRIBUTIONS ACKNO​WLE​DGE​MENTS CONFLICT OF INTEREST STATEMENT DATA AVAILABILITY STATEMENT REFERENCES