As expected, our warming treatment increased average surface temperature by ca. 1°C, a warming similar to the 1.1°C increase of global surface temperature since 1850 (IPCC, 2023) and to the 0.8°C increase in air temperature in the last decade near our study site (Fig. S10). Additionally, it created summer, heatwave-like peaks in which surface temperatures inside OTCs were > 10°C higher than in controls, a common observation in ITEX sites (Bokhorst et al., 2013; Hollister et al. 2023). Our initial hypothesis of a warming-induced acceleration of ecological succession was based on the expected increase in average surface temperatures. However, our findings suggest that ecosystem responses at the CRUST experiment were driven to a large extent by the simulated summer heat pulses.
Heatwaves, which are becoming more frequent in the high north (Dobricic et al., 2020), can be a major source of stress for tundra plants like Salix polaris (Jónsdóttir et al., 2023) and close relatives (Marchand et al., 2005), possibly including S. herbacea, the dominant vascular plant in our study site. Contrary to our hypothesis of vascular plants growing on biocrust, we observed net decreases of vascular plant cover and increases of biocrust cover. It is important to notice that for the period 2018–2021, when we estimated cover using the quadrat method (Fig. S2), an increase of biocrust cover means an increase of the visible biocrust. It is possible that during this period the observed increase of biocrust cover was because of reduction in mosses and vascular plants, thus exposing more of the existing biocrust. In a transplant experiment in the Colorado Plateau (Antoninka et al., 2022), biocrusts from different sources had higher cover and community stability at the highest and coldest elevation. In this sense, the observed increase of biocrust cover in our experiment between 2018 and 2021 could be due to the resistance and/or resilience of Anthelia jurkatzkana biocrust to environmental change. On the other hand, the point intercept that we applied in 2022 accounts for biocrust covered by leaves. Therefore, our observation in 2022 of higher biocrust cover in OTCs (60.4 ± 2.7%) than in control plots (52.6 ± 3.6%) indicates a net increase in biocrust cover. It is not rare for biocrust be more resistant to environmental stress than vascular plants (e.g., Condon et al., 2024). In our case, this observation could be at least partially explained by the capacity of biocrust organisms to rapidly adjust their metabolic rates to sudden changes in the environment (Chongfeng et al., 2017).
Surface cover responses at the CRUST experiment may still be transient. In a similar experiment in northern Sweden, Alatalo et al., (2018) found that the effects of simulated warming on the relative abundance of bryophytes and lichens (two common components of biocrust) varied between years five and seven. In their experiment, warming reduced the relative abundance of bryophytes, with most of the change in the first five years. Warming also decreased the relative abundance of lichens, but in that case most of the reduction happened in the last two years of the experiment. In our case, too, the effects of OTC on surface cover differed between years and the most conspicuous differences between control and manipulated plots were at the end of the experiment, particularly for the cover of vascular plants. Future analyses at the CRUST experiment will show whether the effects of OTCs on surface cover will stabilize (as in the case of the bryophytes in Alatalo et al., 2018) or will intensify (as in the case of lichen in Alatalo et al., 2018, and as with the cover of vascular plants in our study). Also, in a laboratory experiment we found that as long as water is not limiting warming (between 0 and 25°C) can increase N fixation rates in Anthelia biocrust (Salazar et al., 2022). This suggests that as global temperatures continue rising, biocrusts from the highlands of Iceland and similar cold environments may fix more N2 into terrestrial ecosystems, which could, in the long-term, create favorable nutritional conditions for the establishment of vascular plants.
In addition to simulating warming, OTCs exclude grazers. Sheep grazing has transformed the landscape of North-Atlantic ecosystems for centuries (Ross et al., 2016). In the communal highland grazing areas of Iceland, where our experiment is located, sheep grazing is one of the main drivers of ecosystem degradation (Arnalds and Barkarson, 2003; Barrio et al., 2018; Marteinsdóttir et al., 2017). At the CRUST site, we have seen higher abundances of palatable species such as Bistorta vivipara in OTCs than in control plots (Gunnlaugsdóttir, 2022), which might be an effect of excluding selective grazing rather than of warming. However, instead of an increase in total vascular plant cover within OTCs that could have been explained in part by herbivory exclusion, in the last two years of measurements we observed substantial decreases of vascular plant cover. Because of this, we believe that the differences between OTC and control plots in our experiment were caused, to a larger extent, by the environmental manipulation (both warming and drying) rather than by the unintentional grazing exclusion.
Warming may also have changed the relative abundances of the organisms forming the biocrust. A 13 yrs manipulation experiment at the Colorado Plateau showed that an in situ warming of + 2–4°C with infrared lamps decreased the cover of mosses by 15% (which is in line with our findings) and differentially affected dark (i.e. late successional) and light (i.e. early successional) biocrusts (Phillips et al., 2022). Overall, in that experiment warming resulted in a stagnation of the system in an early successional state dominated by light-pigmented cyanobacteria such as Microcoleus vaginatus. We know that on a short-term time scale (hours to days), warming can alter the microbial community and function of the Anthelia jurkatzkana biocrust of our study site (Salazar et al., 2022). Similar analyses but at the multiyear scale remain to be done.
Regarding our analysis of biocrust function, we first need to acknowledge one caveat in our k measurements. For logistic reasons, we incubated the tea bags in the field for ca. 1 year, instead of 90 days as recommended by Keuskamp et al. (2013). This means that when we collected the tea bags the hydrolysable or labile fraction (i.e. the sum of non-polar extraction, water soluble and acid soluble fractions) could have been exhausted or almost exhausted. Since k is calculated based on a direct relationship with the hydrolysable fraction of the rooibos tea, it is possible that the observed lack of effect of the OTCs on k in our experiment was an artifact caused by an underestimation of k values. Also, k decreases with incubation time. Therefore, the differences in k values between years could reflect differences in incubation periods. In this sense, k values from 2018–2019, calculated based on a shorter incubation period than in 2021–2022, and which suggest an expected negative relationship between depth and k (similar to Miralles et al., 2020), are more reliable. In warm drylands, experimental warming can increase decomposition rates by degrading biocrust and therefore increasing soil-litter mixing (Chuckran et al., 2020). We did not see any evidence of biocrust degradation, which could at least partially explain why warming did not have any effect on k rates in our experiment.
In the case of S rates, which are calculated based on a constant hydrolysable fraction of green tea, we have high confidence in our estimates for both measuring periods. Moreover, all tea bags, in control and OTC plots, were under the same conditions for the same period of time, and therefore the observed differences in S between treatments are comparable. Although biocrust effects are often restricted to the top layer of the soil and decrease with depth (Miralles et al., 2020), in our study warming decreased litter mass loss and S both at 2 and 8 cm depth, suggesting a possible “cascading effect” of the biological activity in the biocrust down to deeper layers of the soil. Overall, our study suggests that warming within limits projected for the high north in the coming decades could increase mass loss in and under the biocrust from our site and similar ecosystems, but reduce the capacity of such biocrust to stabilize labile carbon into recalcitrant forms.
Our findings of higher respiration rates, as well as higher litter cover and mass losses, within the OTCs than in control plots, may indicate a warming-induced increase in turnover rates of carbon. Warming generally decreases biocrust respiration in warm drylands (García-Palacios et al., 2018). This is because of higher mortality of biocrust-forming organisms (e.g., lichens) and/or because of microbial thermal acclimation, as indicated by reductions in the activity of extracellular enzymes such as β-glucosidase (García-Palacios et al., 2018). However, in wetter conditions, such as wet years in semi-arid grasslands (Escolar et al., 2015) or in biocrust-dominated ecosystems regularly covered with snow (Darrouzet-Nardi et al., 2015; this study), warming can increase biocrust respiration rates and ultimately increase carbon losses. A previous study in our site suggests that surface temperature and respiration are unaffected by OTCs for most part of the snow-free season (Gunnlaugsdóttir, 2022). However, our interannual comparison of daytime (ca. 11:00–13:00 UTC) biocrust respiration in the peak of the summer (June-August) suggests a marginal, positive effect of the OTC-warming. Overall, this is consistent with a recent synthesis of 56 OTC in-situ experiments located at Arctic and alpine tundra sites that reports an average 30% increase in ecosystem respiration with warming (Maes et al., 2024).
At least part of our observations, such as the negative ΔRB recorded in the dry June of 2019, could have been caused by drying rather than by warming. Mean annual precipitation is projected to increase in Arctic and sub-Arctic regions in the coming decades (IPCC, 2023). However, biocrust-dominated ecosystems such as the one from our experiment will continue experiencing periods of dryness, especially at sites similar to ours where soils have limited water holding capacity and capillary water transfer from sub-surface horizons during droughts is hindered by near-surface coarse tephra layers (i.e. ash; Arnalds, 2015). More data is needed to better understand the ways in which changes in moisture and temperature could affect carbon dynamics in the vast areas of the high north covered with biocrust. Our findings point towards warming increasing carbon emissions from biocrust in the highlands of Iceland, which combined with our observations of higher litter mass loss and lower S values in OTCs than controls, could mean future depletion of soil carbon and possibly other soil nutrients (Ferrenberg et al., 2022).
Warming-induced biocrust development has been reported in many locations of the Arctic and sub-Arctic, often on bare ground after a glacier retreat (e.g., Breen and Levesque, 2006; Yoshitake et al., 2018). In our study, however, the bare ground cover was negligible. What we observed, was an orchestrated response of the different plant-soil components of an understudied type of biocrust habitat. An observation that fills an empty space in the puzzle that is understanding how tundra ecosystems could be further affected by climate change.