# The role of plant functional groups mediating climate impacts on carbon and biodiversity of alpine grasslands

Vigdis Vandvik, Inge H J Althuizen, Francesca Jaroszynska, Linn C Krüger, Hanna Lee, Deborah E Goldberg, Kari Klanderud, Siri L Olsen, Richard J Telford, Silje a H Östman, et al.

► To cite this version:

Vigdis Vandvik, Inge H J Althuizen, Francesca Jaroszynska, Linn C Krüger, Hanna Lee, et al. The role of plant functional groups mediating climate impacts on carbon and biodiversity of alpine grasslands. Scientific Data, 2022, 9 (1), pp.451. 10.1038/s41597-022-01559-0. hal-04037535

HAL Id: hal-04037535

https://hal.science/hal-04037535

Submitted on 20 Mar 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
scientific **data**

<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

DATA DESCRIPTOR

# The role of plant functional groups mediating climate impacts on carbon and biodiversity of alpine grasslands

Vigdis Vandvik ID 1,2✉, Inge H. J. Althuizen ID 1,2,3, Francesca Jaroszynska 1,4, Linn C. Krüger 1, Hanna Lee 2,3,5, Deborah E. Goldberg 6, Kari Klanderud ID 7, Siri L. Olsen 7,8, Richard J. Telford ID 1,2, Silje A. H. Östman 1, Sara Busca 1, Ingrid J. Dahle 1, Dagmar D. Egelkraut 1,2, Sonya R. Geange ID 1,2, Ragnhild Gya ID 1,2, Josh S. Lynn 1,2, Eric Meineri 1,9, Sherry Young 10 & Aud H. Halbritter ID 1,2

Plant removal experiments allow assessment of the role of biotic interactions among species or functional groups in community assembly and ecosystem functioning. When replicated along climate gradients, they can assess changes in interactions among species or functional groups with climate. Across twelve sites in the Vestland Climate Grid (VCG) spanning 4 °C in growing season temperature and 2000 mm in mean annual precipitation across boreal and alpine regions of Western Norway, we conducted a fully factorial plant functional group removal experiment (graminoids, forbs, bryophytes). Over six years, we recorded biomass removed, soil microclimate, plant community composition and structure, seedling recruitment, ecosystem carbon fluxes, and reflectance in 384 experimental and control plots. The dataset consists of 5,412 biomass records, 360 species-level biomass records, 1,084,970 soil temperature records, 4,771 soil moisture records, 17,181 plant records covering 206 taxa, 16,656 seedling records, 3,696 ecosystem carbon flux measurements, and 1,244 reflectance measurements. The data can be combined with longer-term climate data and plant population, community, ecosystem, and functional trait data collected within the VCG.

## Background & Summary

Climate change poses a threat to biodiversity and ecosystem functioning of alpine ecosystems by altering plant, animal and microbial distributions<fnref n="1" />,², community composition³, and food webs⁴; as well as affecting biotic interactions between organisms and functional groups⁵⁻⁸. These impacts operate through both direct and indirect pathways. Individual organism's physiological rates such as photosynthesis and respiration are directly regulated by climatic factors such as temperature and precipitation⁹,¹⁰. In turn, these impacts will aggregate to affect ecosystem-level processes such as productivity, decomposition, and carbon and nutrient fluxes in response to climate change¹¹,¹². Because biodiversity and ecosystem functioning are linked through biotic interactions within and among trophic levels and functional groups⁵,¹³, indirect impacts of climate, operating on and via these interactions, are also of critical importance. Disentangling the direct and indirect consequences of climate on ecosystems is key to understanding and predicting the full spectrum of biodiversity and functional responses to a changing climate.

Plants are key players in ecosystems: they make up the majority of terrestrial biomass globally, are the main source of photosynthetically fixed carbon, and serve as habitat, food sources, and hosts for other organismal

<footnote n="1">¹Department of Biological Sciences, University of Bergen, Bergen, Norway. ²Bjerknes Centre for Climate Research, University of Bergen, Bergen, Norway. ³Norwegian Research Centre, Bjerknes Centre for Climate Research, Bergen, Norway. ⁴Office Français de la Biodiversité, Pérols, France. ⁵Department of Biology, Norwegian University of Science and Technology, Trondheim, Norway. ⁶Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, USA. ⁷Faculty of Environmental Sciences and Natural Resource Management, Norwegian University of Life Sciences, Ås, Norway. ⁸Norwegian Institute for Nature Research, Oslo, Norway. ⁹Aix Marseille Univ, Avignon Univ, CNRS, IRD, IMBE, Marseille, France. ¹⁰Simon Fraser University, Vancouver, Canada. ✉e-mail: vigdis.vandvik@uib.no</footnote>
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Study area, site selection, experimental design, and field sampling overview for the FunCaB plant functional group experiment. (a) Location of the study area and the 12 study sites in Vestland county, Western Norway. (b) The Vestland Climate Grid is established across independent broad-scale biogeographic gradients in summer temperature (mean of four warmest months, °C), and annual precipitation (mean annual precipitation, mm). Temperature levels are given as boreal, sub-alpine, and alpine; precipitation levels range from 1 (dry) to 4 (wet). (c) Experimental design, with the timeline (species composition recordings [dataset v], seedling recruitment monitoring [vi], ecosystem carbon flux monitoring [vii], and functional group removals [i] indicated, see legend for dataset symbols, for other datasets see text), an overview of the eight factorial removal treatments, and plot layout with subplots used for the community and seedling recording. The removal treatments are described by the functional groups removed from the respective plots, note that in the figure the 'Intact' plot refers to the no removal controls whereas the 'Gap' are plots with all functional groups removed. For detailed description of treatments and their abbreviations as used in the datasets, see text.</figcaption></figure>

groups, while also competing with them for resources¹³⁻¹⁵. A changing climate will drive predictable shifts in plant species composition; warmer and wetter climates generally favour taller plants, more resource-acquisitive traits, and species that are better competitors for light and nutrients¹⁶. Previous experiments and observational studies from arctic and alpine systems show that plant functional groups encapsulate much of this variation, and also respond differently to climate and climate change, with graminoids and shrubs typically increasing¹⁷⁻¹⁹, and bryophytes and lichens decreasing in response to warming¹⁷,¹⁸,²⁰,²¹. These shifts have knock-on effects on ecosystem productivity, moisture availability, mineralization, and decomposition processes²²,²³.

Removal experiments are a 'gold standard' approach for disentangling the roles, interactions, and effects of co-occurring (groups of) plants (see¹³,²⁴). Here we report on a fully factorial removal experiment to assess the roles of and interactions between the major different plant functional groups in alpine grasslands - graminoids, forbs, and bryophytes - in controlling alpine biodiversity, microclimate, and ecosystem functioning. Our semi-natural alpine grassland study system is of high biodiversity and conservation value²⁵, yet structurally simple and fine-scaled enough to allow whole-ecosystem (plants and soil community) experimental manipulation. Our experiment was replicated across broad-scale temperature and precipitation gradients in Western Norway in a 'climate grid' design that allows disentangling effects of temperature and precipitation change. This is useful, as temperature and precipitation are often confounded, both along elevational gradients and in warming experiments (changes in moisture availability are a side-effect of standard experimental warming approaches, see²⁶). We use this study design to disentangle the direct and indirect effects of climate on plant taxonomic and functional diversity, microclimate, and carbon dynamics of alpine grasslands by comparing treatments with plant functional groups removed (direct effects of climate) with treatments with one or more functional groups present (net of direct and indirect climate effects) across sites.

Our experiment is replicated across a macroecological context of twelve calcareous semi-natural grassland sites in south-western Norway (Fig. 1a). The fjords of western Norway offer large variation in temperature and precipitation, and we used this variability to select sites in a 'climate grid' where three temperature levels are crossed with four precipitation levels (Fig. 1b). The Vestland Climate Grid (VCG), was established in 2008, and a host of biotic and abiotic data are available from previous research, as detailed below. Within each VCG site, we established the FunCaB plant functional group removal experiment from which we measured carbon and biodiversity responses over 7 years (Fig. 1c). Below, we first describe the site selection and provide basic site information, then outline the experimental design and setup, before we detail variables related to (i) functional group biomass and the removal treatments, (ii) species-level biomass, (iii-iv) soil microclimate, (v) plant community composition, (vi) seedling recruitment, (vii) ecosystem carbon fluxes, and (viii) ecosystem reflectance (Table 1). Note that these data vary in spatial and temporal coverage, as described below. We provide raw and clean datasets with consistent structure and variable naming²⁷, and associated code for cleaning and combining data²⁸. By documenting and communicating data structures and qualities we aim to facilitate data reuse and combination for new applications in the future.
<table><caption>**Table 1.** Description and location of the datasets in the FunCaB plant functional group experiment. This table summarises information on dataset number, response variable(s), number of observations, temporal range of the data, and location of the primary data, the final published data, and the code for extracting and cleaning data from the primary data.</caption><thead><tr><th>Dataset</th><th>Response variable</th><th>Number of observations</th><th>Temporal range</th><th>Citation information for raw data, clean data, and code</th></tr></thead><tbody><tr><td>i</td><td>Functional group removal</td><td>5,412</td><td>2015 – 2021</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>ii</td><td>Species-level biomass</td><td>360</td><td>2016</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>iii</td><td>Soil temperature</td><td>1,084,970</td><td>2015 – 2016</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>iv</td><td>Soil moisture</td><td>4,771</td><td>2015 – 2019</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>v</td><td>Plant community composition</td><td>plant records: 17,181; plant taxa: 206</td><td>2015 – 2019</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>vi</td><td>Seedling recruitment</td><td>16,656</td><td>2018 – 2019</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>vii</td><td>Ecosystem carbon fluxes</td><td>3,696</td><td>2015 – 2017</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr><tr><td>viii</td><td>Reflectance</td><td>1,244</td><td>2019, 2021</td><td>Raw data<sup>27</sup>, clean data<sup>27</sup>, code<sup>28</sup></td></tr></tbody></table>

<figure><img src="image_4.png" /><figcaption>**Fig. 2** The data collection and management workflow of the FunCaB project. Reproducibility throughout the research process is assured as follows: Experimental design and data collection was based on best-practice community methods and protocols, adapted for the projects' needs. Measurements were digitalized and the raw data stored in the project Open Science Foundation (OSF) repository before the raw data were cleaned and managed through code-based data curation, with version control secured via GitHub. The clean data are stored at the OSF repository, and a time-stamped version of the code to retrieve and clean data is provided through Zenodo. This data paper describes and documents the data collection and workflow, and describes how to access and use clean data, raw data, and code.</figcaption></figure>

## Methods

**Data management and workflows.** We adopt best-practice approaches for open and reproducible research planning, execution, reporting, and management throughout the project (see e.g.<sup>29–32</sup>). Specifically, we use community-approved standards for experimental design and data collection<sup>29</sup>, and clean and manage the data using a fully scripted and reproducible data workflow, with data and code deposited at open repositories (Fig. 2).

**Research site selection and basic information, and general study setup.** Site selection. Our study is conducted across the twelve calcareous grassland experimental sites in the Vestland Climate Grid (VCG), in south-western Norway (Fig. 1a). The VCG sites were chosen to fit within a climate grid reflecting a fully factorial design encompassing the major bioclimatic variation in Norway. Potential sites were identified using a combination of topographic maps, geological maps (NGU) and interpolated maps of summer temperature and annual precipitation using the 1960–1990 climate normal (100 m resolution gridded data, met.no; see<sup>33</sup> and references therein). The three temperature levels (alpine, sub-alpine, boreal) and four levels of precipitation in the climate grid (Fig. 1b) were selected to reflect a difference in mean growing season temperature of ca. 2 °C between three temperature levels (alpine = 6.5 °C, sub-alpine = 8.5 °C, boreal = 10.5 °C mean temperature of the four warmest months of the year) and a difference in mean annual precipitation of 700 mm between four precipitation levels (precipitation levels 1 – 4 representing 700 mm, 1400 mm, 2100 mm, and 2800 mm, respectively). Climate data for the site selection was based on 100-m resolution downscaled data using the 1960–1990 climate normal from met.no. The final sites were selected from approximately 200 potential sites visited and surveyed in the summer of 2008, with selection criteria set to ensure that other factors such as grazing regime and history, bedrock, vegetation type and structure, slope and exposure were kept as constant as possible among the selected sites<sup>34</sup>. Geographical distance between sites is on average 15 km and ranges from 175 km to 650 m.
<table><caption>**Table 2.** Data dictionary for the FunCaB functional group-level biomass (dataset i). Data dictionary with column descriptions for dataset i – the biomass removed from 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2021, and total biomass in 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway.</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2015 – 2021</td><td>yyyy</td><td>defined</td></tr><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2015-07-22 – 2021-07-30</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>round</td><td>Round of sampling</td><td>numeric</td><td>1 – 2</td><td>count</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust-Vikesland</td><td></td><td>defined</td></tr><tr><td>temperature_level</td><td>Temperature level with the climate grid</td><td>categorical</td><td>boreal, sub-alpine, alpine</td><td></td><td>defined</td></tr><tr><td>precipitation_level</td><td>Precipitation level within the climate grid (low to high)</td><td>categorical</td><td>1, 2, 3, 4</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 – Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B – Vik5GF</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes, C = control and XC = extra control</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C, XC</td><td></td><td>defined</td></tr><tr><td>removedFG</td><td>Removed functional group, where F = forbs, B = bryophytes, G = graminoids. For extra controls also L = litter, P = pteridophytes, LI = lichens, and C = cryptograms</td><td>categorical</td><td>For experimental plots: F, G, B, For extra controls: F, G, B, L, P, LI, C</td><td></td><td>defined</td></tr><tr><td>functional_group</td><td>Removed functional group, including forbs, bryophytes, graminoids. For extra controls also litter, pteridophytes, lichens, and cryptograms</td><td>categorical</td><td>bryophytes-pteridophytes</td><td></td><td>defined</td></tr><tr><td>biomass</td><td>Dry weight of removed functional_group</td><td>numeric</td><td>0 – 41.31</td><td>g</td><td>measured</td></tr><tr><td>name</td><td>Name of data collector</td><td>categorical</td><td>AB – William</td><td></td><td>recorded</td></tr><tr><td>remark</td><td>Remarks</td><td>categorical</td><td></td><td></td><td></td></tr></tbody></table>

*Study system and experimental area selection within sites.* At each site, we selected an experimental area of ca. 75 – 200 m², targeting a homogeneous and representative part of the target grassland vegetation at large at that site. The experimental areas were placed on southerly-facing slopes, avoiding depressions and concave areas in the landscape and other features such as big rocks or formations that may affect light conditions, hydrology and/or snowdrift. The target vegetation type was forb-rich semi-natural upland grassland vegetation³⁴,³⁵, within the plant sociological association Potentillo-Festucetum ovinae tending towards Potentillo-Poligonium vivipara in the alpine sites and Nardo-Agrostion tenuis in some lowland sites³⁶. The most common vascular plants across sites, based on sum of covers, are the graminoids *Agrostis capillaris*, *Festuca rubra*, *Avenella flexuosa*, *Anthoxanthum odoratum*, and *Nardus stricta* and the forbs *Leucazentem vulgare*, *Hypericum maculatum*, *Silene acaulis*, *Alchemilla alpina*, and *Lotus corniculatus*. Common bryophytes are *Pleurotium schreberi*, *Hylocomium splendens*, *Polythricum spp*, *Racomitrium lanuginosum*, *R. fasciculare*, and *Dicranum* spp. All sites were moderately grazed prior to the study by sheep, cattle, goats, reindeer, deer, moose, and/or horses; and the experimental areas were fenced for the duration of the study to prevent animal and human disturbance of the experimental infrastructure. The fenced area was lightly mowed at the end of each growing season to mimic past grazing pressure and minimize fence effects. For further description of the sites, see³⁴ and for access to and further description of site-level data, see³⁵.

*Block and experimental plot setup.* Within these study areas we established four blocks, with a distance between the blocks ranging from one up to (in rare cases) 50 meters. Blocks were selectively placed in homogenous grassland vegetation, avoiding rocks, depressions, and other features as described above. Each block accommodates eight 25 × 25 cm plots, with at least 25 cm between adjacent plots. If a plot contained more than 10% bare rock, shrubs, or other non-grassland features, they were rejected or moved slightly to avoid these features. The plots were permanently marked with four aluminium 10-cm long pipes in the soil in the outer corners of all the 25 × 25 cm treatment plots, ensuring the pipes to fit the corners of a standardized vegetation analysis frame (aluminium frame demarking a 25 × 25 cm inner area, with poles fixed in the corners that fit into the aluminium tubes used for plot demarcation in the field). The upslope left corner tube was marked with a colour-coded waterproof tape. Note that in 31 out of 48 cases (12 sites × 4 blocks), the blocks were located within larger experimental blocks in the VCG sites, and control plots and various block-level data are then shared with other experiments in these larger blocks. Linking keys are described in the FunCaB data dictionaries below (see Fig. 3 and data records iii-vii below). For some datasets, additional plots within blocks were needed. These are described as needed below.
<table><caption>**Table 3.** Data dictionary for the FunCaB species-level biomass (dataset ii). Data dictionary with column descriptions for dataset ii – the total forb biomass per species in the 48 extra controls in 2016 at twelve sites in the Vestland Climate Grid, Vestland county, Norway.</caption><thead><tr><td>Variable name</td><td>Description</td><td>Variable type</td><td>Variable range or levels</td><td>Units</td><td>How measured</td></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2016 - 2016</td><td>yyyy</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust - Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 - Vik4</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1XC - Vik4XC</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes and C = control, and XC = extra control</td><td>categorical</td><td>XC - XC</td><td></td><td>defined</td></tr><tr><td>functional_group</td><td>Removed functional group, including forbs, bryophytes, graminoids, litter, pteridophytes, lichens, and cryptograms</td><td>categorical</td><td>forb - forb</td><td></td><td>defined</td></tr><tr><td>species</td><td>Species abbreviation</td><td>categorical</td><td>Ach.mil - Vio.sp</td><td></td><td>defined</td></tr><tr><td>biomass</td><td>Dry weight of removed functional group</td><td>numeric</td><td>0.01 - 12.95</td><td>g</td><td>measured</td></tr><tr><td>sorted_by</td><td>Person that sorted the biomass</td><td>categorical</td><td>AP - PP</td><td></td><td>recorded</td></tr></tbody></table>

**Background abiotic and biotic data from the Vestland Climate Grid.** The Vestland Climate Grid field sites were established in 2008, and from a series of research projects within the grid over the years we have collected a broad range of datasets on the climate and environment, soils, land-use and environment, vegetation, and ecosystems, along with basic descriptive data of the 12 sites, as described in³⁴. All these datasets are available from the previous projects through the VCG OSF (Open Science Framework) repository³⁵, and the results are presented in associated papers, see for example³⁴,³⁷⁻⁴⁵. The overall data structure, and the most relevant datasets from the VCG for the FunCaB project is laid out in Fig. 3, and briefly described below. Code to download and link these data to the FunCaB experimental data and sites are provided in the FunCaB github repository²⁸ (see R/download_VCG_data).

A new research project, 'FUNDER - Direct and indirect climate impacts on the biodiversity and FUNCTIONing of the UNDERground ecosystem' funded by the Norwegian Research Council KLIMAFORSK programme (project number 315249, 2021 – 2025) will augment the FunCaB experiment with data on the belowground components of the plant-soil ecosystem, including roots, mesofauna, fungi and microbes. These upcoming data will all link with the FunCaB and VCG project based on the given experimental, site and organismal keys, as indicated in Fig. 3.

**VCG Basic site-level attributes.** Basic descriptive data on the 12 sites include latitude, longitude, elevation, geology, land-use, soils, and their position within the climate grid (precipitation and temperature levels). These data are described in³⁴,⁴⁰, provided in³⁵, and can be downloaded using²⁸ (see R/download_VCG_data). For convenience, the climate grid information is also provided in the biomass dataset (see below).

**VCG Site-level climate data.** Temperature was measured continuously at each of the 12 VCG sites at four heights (2 m and 30 cm above ground, at ground level, and 5 cm below ground), soil moisture was measured continuously with two replicate loggers ca. 5 cm below ground, and precipitation was measured at each site during the snow-free season. For these measurements, we used Delta T GP1 loggers (Delta T devices, Cambridge, UK) equipped with two temperature probes, two SM200 moisture sensors which were later replaced as necessary with SM300 and SM150T loggers, and an ARG 100 tipping bucket (EML LTD, North Shields, UK) from 2009 onwards. UTL-3 version 3.0 temperature loggers (GEOTEST AG, Zollikofen, Switzerland) were used for measuring the 2 m and 30 cm temperatures. Soil moisture was measured as the mean of four measurements taken along each side of the turf, several times during the growing season using a Delta T HH2 version 2.3 Moisture Meter with the same probes as for the GP1 logger (SM200, SM150T). These data are described in³⁴,⁴⁰, provided in³⁵, and can be downloaded using²⁸ (see folder R/download_VCG_data).

**VCG Soil chemical and structural data.** Over the years, various soil chemical variables have been measured at the block level within each of the 12 VCG sites, including soil pH (2009) and % Loss-On-Ignition (2009, 2013), and available N, as sum of N available as NH4-N and NO3-N (available N per deployment period, 2010 & 2012). Soil pH was measured after adding 50 ml distilled water to 25 g soil and mixing for two hours. Loss-on-ignition (LOI), was measured by weighing dry soil (105 °C for 24, one hour in a desiccator), and burnt soil (six hours at 550 °C, one hour in the desiccator) and calculating LOI as the (burnt soil mass/dry soil mass) × 100. NH4-N and NO3-N were measured using in-situ ion exchange resin bags (IERBs) were used to measure the amount of plant-available nutrients in the soil. These data are partially described in³⁴,⁴⁰, and the full data are provided in³⁵.

**VCG Litter decomposition data.** Decomposition has been assessed at each of the 12 VCG sites using local plant litter and the Tea Bag Index method (Keuskamp et al., 2013). Local litter (dead leaves detached from live plants) was collected at each site in 2013 or 2014, with the specific timing of the collections at each site tuned to ensure that litter was present, not covered by snow, and not decomposed. In practice, this necessitated litter collection after snowmelt in spring in many sites. The litter was washed, dried, and stored in dark, dry, cool conditions. In 2016,
<table><caption>**Table 4.** Data dictionary for the FunCaB soil temperatures (dataset ii). Data dictionary with column descriptions for dataset iii – soil temperature from 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2016 at twelve sites in the Vestland Climate Grid, Vestland county, Norway.</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2015 – 2016</td><td>yyyy</td><td>defined</td></tr><tr><td>date_time</td><td>Date and time of sampling (UTC + 2)</td><td>date_time</td><td>2015-07-12 17:58:01 – 2016-06-27 23:05:01</td><td>yyyy-mm-dd hh:mm:ss</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust – Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 – Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B – Vik5GF</td><td></td><td>defined</td></tr><tr><td>iButtonID</td><td>Unique iButton ID</td><td>categorical</td><td>003E3B5C41_2016 – F83E3E2A41_2016</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes, and C = control</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C</td><td></td><td>defined</td></tr><tr><td>soiltemperature</td><td>Soil temperature measurement per plot</td><td>numeric</td><td>-8.75 – 32.8</td><td>°C</td><td>recorded</td></tr><tr><td>comments</td><td>Additional comments</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>turfID</td><td>Unique turf ID from SeedClim transplant experiment</td><td>categorical</td><td>101 TTC – 73 TTC</td><td></td><td>defined</td></tr></tbody></table>

five replicate litter bags containing 1 g of graminoid litter were buried at each site, and collected at four points in time after burial (1, 2, 3 and 12 months). Harvested litter bags were cleaned (soil and roots removed), dried for 48 h at 60 °C and weighed. The Tea Bag Index method⁴⁶ was used in 2014, 2015 and 2016 to measure decomposition at all sites of the climate grid. At each site, 10 replicates of each tea type were buried pair-wise, 8 cm below ground and with at least 10 cm between the tea bags. For a couple of sites, the number of replicate tea bag pairs was higher in 2015 (12 replicates at the site Gudmedalen and 16 replicates at Låvisdalen). After collection, adhering soil particles and roots were removed and the tea bags were dried (48 h at 60 °C) and weighed. These data are partially described in⁴⁷, and the full data are provided in³⁵ and can be downloaded using²⁸ (see folder R/download_VCG_data).

**VCG Species-level cover, biomass, and performance data.** A variety of plant species and community composition, cover, biomass, fitness, and reproductive data exists for the sites and blocks in the VCG from 2008 to 2021. These data are described in e.g³⁴,³⁷,³⁸,⁴¹,⁴³–⁴⁵,⁴⁸–⁵⁰, and provided in³⁵.

**VCG Site-level plant functional traits.** In 2016 and 2017, we measured 11 leaf functional traits that are related to potential physiological growth rates and environmental tolerance of plants, following the standardized protocols in Pérez-Harguindeguy et al.⁵¹: leaf area (LA, cm²), leaf thickness (LT, mm), leaf dry matter content (LDMC, g/g), specific leaf area (SLA, cm²/g), carbon (C, %), nitrogen (N, %), phosphorus (P, %), carbon nitrogen ratio (C:N), nitrogen phosphorus ratio (N:P), carbon13 isotope ratio (δ13C, ‰), and nitrogen15 isotope ratio (δ15N, ‰). Trait data are available at the site level for species making up at least 80% of the vegetation cover in the control plots at each of the 12 VCG sites. The plants were collected outside of the experimental plots and within a 50 m perimeter from the blocks, and we aimed to collect up to five individuals from each species in each site. To avoid repeated sampling from a single clone, we selected individuals that were visibly separated from other ramets of that species. The sampled plant individuals were labelled, put in plastic bags with moist paper towels, and stored in darkness at 4 °C until processing, which was done as soon as possible and always within 4 days. These data are described in⁵², provided in³⁵, and can be downloaded using²⁸ (see folder R/download_VCG_data).

**Experimental design.** The functional group removal experiment was designed to examine the impact of aboveground interactions among the major plant functional groups - graminoids, forbs and bryophytes - on the performance and functioning of other components of the vegetation and ecosystem. The experiment consists of eight 25 × 25 cm plots per site and block, with a fully factorial combination of removals of three plant functional groups, with treatments randomized within each block. The general experimental design, with the different removal treatments detailed, are provided as an insert to the timeline in Fig. 1c. The functional groups are delineated and abbreviated in the various datasets as follows: G = graminoids (including grasses, sedges and rushes), F = forbs (including herbaceous forbs, pteridophytes, dwarf-shrubs, and small individuals of trees and shrubs), B = Bryophytes (including mosses, liverworts, and hornworts). Note that all species are also coded by their respective functional group into which they were classified in the FunCaB taxon table. The experimental treatments are coded by functional group removed so that FGB = bare-ground gaps with all plants removed, FB = only graminoids remaining, GB = only forbs remaining, GF = only bryophytes remaining, B = graminoids and forbs remaining, F = bryophytes and graminoids remaining, G = bryophytes and forbs remaining, and C = intact vegetation controls with no vegetation removed. In 2016, four extra control (XC) plots were marked per site for aboveground biomass harvest and ecosystem carbon flux measurements. This sampling regime gave a total of 384 plots in the core FunCaB experiment, plus the additional 48 controls in 2016.

Functional group removals were done once in 2015 (at peak growing season due to late snowmelt), twice per year in 2016 and 2017 (after the spring growth and at peak growing season) and annually from 2018 to 2021 (at
<table><caption>**Table 5.** Data dictionary for the FunCaB soil moisture (dataset iv). Data dictionary with column descriptions for dataset iv – soil moisture from 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2016 at twelve sites in the Vestland Climate Grid, Vestland county, Norway.</caption><thead><tr><td>Variable name</td><td>Description</td><td>Variable type</td><td>Variable range or levels</td><td>Units</td><td>How measured</td></tr></thead><tbody><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2015-06-02 – 2019-08-09</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust – Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr – Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B – Vik5GF</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes, and C = control</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C</td><td></td><td>defined</td></tr><tr><td>soilmoisture</td><td>Soil moisture measurement per plot</td><td>numeric</td><td>0 – 100</td><td>(m³ water × m⁻³ soil) × 100</td><td>recorded</td></tr><tr><td>weather</td><td>Weather conditions during data collection</td><td>categorical</td><td>Blue sky - windy, cloudy</td><td></td><td>recorded</td></tr><tr><td>recorder</td><td>Data collector</td><td>categorical</td><td>? – Vojta</td><td></td><td>recorded</td></tr><tr><td>turfID</td><td>Unique turf ID from SeedClim transplant experiment</td><td>categorical</td><td>101 TTC – TTC 281</td><td></td><td>defined</td></tr></tbody></table>

peak growing season) as regrowth had declined (see below) and biannual removals were no longer necessary. At each sampling, all above-ground biomass of the relevant plant functional group was removed from each plot as follows: for each plot, all the above-ground parts of the relevant functional group(s) were removed using scissors and tweezers to cut the plants at the ground layer (i.e., the soil-vegetation interface). Roots and other below-ground parts were not removed, and non-target plant functional groups and litter were left intact.

**Species identification, taxonomy, and flora.** All vascular plant species were identified to the species level in the field, with nomenclature following Lid and Lid⁵³. Exceptions are sterile specimens of species that are not possible to identify without reproductive parts, and where flowers are either too rare or individuals too short-lived for comparisons of the position of individuals within the plots over years to be used to ascertain identifications (For example, *Alchemilla* spp. excluding *A. alpina*, and the annual *Euphrasia* spp.). Species identifications were confirmed by comparing records over time as described below. All unidentified specimens are included and flagged in the dataset, as described in Data Records below. The full taxon names are provided in the taxon table on OSF (Fig. 3).

**Dataset collection methods.** *Datasets (i–ii): Biomass and functional group removal.* As described above, functional group removals were done once in 2015 at peak growing season, and twice per year in 2016 and 2017 (after the spring growth; at peak growing season) and annually at peak growing season from 2018 to 2021. For each removal plot and occasion, a picture was taken of the plot pre-removal, the biomass to be removed was collected in separate pre-marked paper bags for each functional group (graminoids, forbs and bryophytes), and a picture was taken post-removal. The collected biomass was then dried at 60°C for 48 hours and weighed to the nearest 0.01 g (Model LPG-1002, VWR). From the four extra control (XC) plots in 2016, total above-ground biomass as well as litter (defined as dead biomass detached from live plants, see²⁸) was collected at peak growing season. From these plots, biomass was sorted into functional groups as described above, except the forb functional group, which was sorted into species. The graminoid and bryophyte functional groups, each forb species, and litter were individually dried and weighed as described above. The data is available as (i) a biomass dataset, consisting of the removed biomass per plot, date, removal treatment, and functional group for all treatment plots, and the total biomass per functional group plus litter for the extra control plots in 2016, and (ii) a species-level forb biomass dataset from the extra control plots in 2016 (Fig. 3, Table 1).

**Datasets (iii–iv) – Soil microclimate.** We measured soil temperature 3–5 cm below the soil surface for each plot using iButton temperature sensors (DS1922L, Manufacturer reports temperature accuracy of ±0.5°C, Maxim Integrated INC., San Jose, CA, USA). The data are reported with a resolution of 0.0625 at 140 min intervals from June 2015 to July 2016. We measured soil moisture as volumetric soil moisture; expressed as % water volume per soil volume ((m³ water /m³ soil) × 100). These measurements were done c. 3–5 times during the growing seasons from 2015–2019, usually in connection with the flux and vegetation measurements, by taking the average of four measurements, one at each side of each plot (SM300, Manufacturer reports accuracy ±2.5% vol over 0 to 50% vol and 0–60°C, Delta-T Devices, Cambridge, UK). The data is available as (iii) temperature and (iv) volumetric soil moisture % per plot and time point (temperature) or date (moisture) (Fig. 3, Table 1).

**Dataset (v): Vascular plant community composition and vegetation structure.** We recorded the full vascular plant species composition of all experimental plots in 2015 (pre treatment), and the control plots plus the extra control plots in 2016. In 2017, 2018, and 2019, we recorded the community composition in controls and in the functional groups that remained in the experimental plots according to the plot's treatment. At each analysis, each plot was sub-divided into 25 5 × 5 cm subplots, using a subplot overlay. We first recorded all species of vascular plants in the central five subplots, (i.e., the central + shaped area of each plot, Fig. 1c) noting the subplot cover of each species present in each of the five subplots (1 – 25% = 1, 26 – 50% = 2, 51 – 75% = 3, >76% = 4). Additionally, we noted if the individual was fertile (records circled if buds, flowers, or fruits were present). The five subplots were recorded and numbered (1-5) by row, and from left to right, starting from the
<figure><img src="image_5.png" /><figcaption>**Fig. 3** Data structure for the FunCaB functional group removal experiment and associated Vestland Climate Grid (VCG) and FUNDER project data. Within each of the three projects, boxes represent data tables. The FunCaB project data tables include biomass of functional groups removed and forbs species-level biomass (datasets i, ii), soil temperature and moisture (datasets iii, iv) plant community composition and the associated taxon table (dataset v), seedling recruitment (dataset vi), ecosystem carbon fluxes (dataset vii) and reflectance (dataset viii). Names of individual data tables are given in the coloured title area, and a selection of the main variables available within tables in the internal lists. For full sets of variables for each FunCaB dataset, see Tables 3–9. The lines linking three of the boxes exemplify links using **species** as keys across tables, note that all bold variables are shared between several tables and can be used as keys to join them. Keys can also be used to link to/from data from other projects in the VCG (for general VCG project keys, see top right hatched outline box, for keys between the FunCaB and FUNDER projects see the bottom right hatched outline box (both including an example value for each variable on the right). The (other) datasets* boxes refer to extensive datasets on plant community composition, cover, biomass, fitness, and reproduction available from previous projects in the VCG²⁷ and upcoming datasets in the FUNDER project.</figcaption></figure>

top up-slope subplot. For the entire 25 × 25 cm plot, any additional species not present in one of the central sub-plots were recorded and their fertility noted. We then visually estimated the percentage cover of each vascular plant species in the whole plot to the nearest 1% and measured vegetation height in mm at four points within the plot. Note that the total coverage in each plot can exceed 100% due to layering of the vegetation. The vascular plant vegetation data is available as percentage cover and fertility status (sterile or fertile) per species per subplot and plot per sampling date, and vegetation height in mm per plot per sampling date (Fig. 3, Table 1).
<table><caption>**Table 6.** Data dictionary for the FunCaB plant community composition (dataset v). Data dictionary with column descriptions for dataset v – the plant community composition of 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2019, and in 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway.</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2015 - 2019</td><td>yyyy</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust - Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 - Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B - Vik5GF</td><td></td><td>defined</td></tr><tr><td>removal</td><td>Pre (2015) or post (&gt;2015) removal treatment. Pre removal is related to the start of the experiment.</td><td>categorical</td><td>post - pre</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = Graminoids, and B = Bryophytes. C = Controls and XC = extra controls</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C, XC</td><td></td><td>defined</td></tr><tr><td>total_graminoids</td><td>Cover of graminoids</td><td>numeric</td><td>0 - 100</td><td>percentage</td><td>recorded</td></tr><tr><td>total_forbs</td><td>Cover of forbs</td><td>numeric</td><td>0 - 100</td><td>percentage</td><td>recorded</td></tr><tr><td>total_bryophytes</td><td>Cover of bryophytes</td><td>numeric</td><td>0 - 120</td><td>percentage</td><td>recorded</td></tr><tr><td>vegetation_height</td><td>Height of vegetation</td><td>numeric</td><td>0 - 350</td><td>mm</td><td>measured</td></tr><tr><td>moss_height</td><td>Height of bryophytes</td><td>numeric</td><td>0 - 120</td><td>mm</td><td>measured</td></tr><tr><td>litter</td><td>Cover of litter</td><td>numeric</td><td>1 - 100</td><td></td><td>recorded</td></tr><tr><td>species</td><td>Species abbreviation</td><td>categorical</td><td>Ach.mil - Vio.tri or NA</td><td></td><td>defined</td></tr><tr><td>cover</td><td>Individual species cover</td><td>numeric</td><td>1 - 98</td><td>percentage</td><td>recorded</td></tr><tr><td>functional_group</td><td>Removed functional group, including forbs, bryophytes, graminoids. For extra controls also litter, pteridophytes, lichens, and cryptograms</td><td>categorical</td><td>forb - graminoid</td><td></td><td>defined</td></tr><tr><td>sumcover</td><td>Total sum of species cover</td><td>numeric</td><td>1 - 176</td><td>percentage</td><td>recorded</td></tr><tr><td>recorder</td><td>Data collector</td><td>categorical</td><td>Aud - W</td><td></td><td>recorded</td></tr><tr><td>turfID</td><td>Unique turf ID from SeedClim transplant experiment</td><td>categorical</td><td>101 TTC - 73 TTC</td><td></td><td>defined</td></tr></tbody></table>

Other variables that were measured were percentage cover of bryophytes, litter, bare ground, and rock (measured per plot and per subplot) and moss layer depth in mm (mean of 4 measurements/plot), date of analysis, recorder/scribe (if any), and free-text comments. These data are available as % cover, depth in mm, date (year, month.day) and text strings per subplot and /or plot per sampling date (Fig. 3, Table 1).

**Dataset (vi): Seedling recruitment.** The total number of forb seedlings that emerged in the plots was recorded in 2018 and 2019. At peak growing season in 2018 (round 1, July-August, depending on site), all dicotyledonous seedlings were marked with wooden toothpicks and their x and y coordinates in the plot (mm, recorded from the bottom left hand-corner of the plot, Fig. 1c) and tentative species identity noted. Toward the end of the growing season (round 2, August-September, depending on site), each plot was revisited, seedling survival established, and any further seedlings marked. Survival (recorded when a seedling was present in subsequent surveys; recorded as mortality if absent) and new seedling emergence were followed up in the same manner in 2019 (rounds 3 and 4, respectively). Species identification was (re)assessed at all censuses and corrected if needed as the seedlings grew and identification uncertainty decreased. New seedlings were differentiated from emergent clonal ramets by looking for cotyledons or signs of above- or below-ground ramet connections. These data are available as talleys of seedlings, each with a status (dead or alive) and species identity (or NA when not identifiable), per subplot and /or plot per sampling round (Fig. 3, Table 1).

**Dataset (vii): Ecosystem carbon flux data and flux calculations.** Carbon flux measurements. Ecosystem CO₂ fluxes were measured to estimate net ecosystem exchange (NEE), ecosystem respiration (Reco) and gross primary production (GPP). The dataset covers the years 2015, 2016 and 2017, and individual plots have multiple measurements for ecosystem carbon flux per year as detailed below. At peak growing season in 2015, a median of 2 sets of paired carbon flux measurements were measured pre-removal for all plots, where a paired set consists of a light and a dark flux measurement of an individual plot. In 2016, a median of 8 sets of paired measurements were made for all control plots, and a median of 7 for the 4 extra controls (see experimental design above). In the data files, some additional measurements exist for other experiments in the VCG sites (a median of 7 paired sets of measurements for controls (TTC) and graminoid removal plots (RTC), see⁴² for a presentation of this experiment and³⁵ for technical details). In 2017, a median of 5 paired sets of measurements were made for all treated plots in nine of the sites, excluding the second wettest precipitation level (sites Gudmedalen, Rambera,
<table><caption>**Table 7.** Data dictionary for the FunCaB seedling recruitment (dataset v). Data dictionary with column descriptions for dataset v – the plant community composition of 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2019, and in 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway.</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2018 - 2019</td><td>yyyy</td><td>defined</td></tr><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2018-07-01 - 2019-08-15</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust - Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 - Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B - Vik5GF</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes and C = control, and XC = extra control</td><td>categorical</td><td>B - GF</td><td></td><td>defined</td></tr><tr><td>seedID</td><td>Unique seedling ID</td><td>categorical</td><td>1_Alr3B - z_Ves3GB</td><td></td><td>defined</td></tr><tr><td>round</td><td>Round of sampling</td><td>numeric</td><td>1 - 4</td><td>count</td><td>defined</td></tr><tr><td>species</td><td>Species abbreviation</td><td>categorical</td><td>Ach.mil - Vio.sp</td><td></td><td>defined</td></tr><tr><td>presence</td><td>Presence (1) or absence (0) of a seedling per census</td><td>numeric</td><td>0 - 1</td><td></td><td>recorded</td></tr><tr><td>x</td><td>x coordinate in the plot</td><td>numeric</td><td>0 - 250</td><td>mm</td><td>recorded</td></tr><tr><td>y</td><td>y coordinate in the plot</td><td>numeric</td><td>0 - 250</td><td>mm</td><td>recorded</td></tr><tr><td>comment</td><td>Comment on measurement</td><td>categorical</td><td>? - x220, y113</td><td></td><td>defined</td></tr><tr><td>turfID</td><td>Unique turf ID from SeedClim transplant experiment</td><td>categorical</td><td>11 TTC - 73 TTC</td><td></td><td>defined</td></tr><tr><td>functional_group</td><td>Removed functional group, including forbs, bryophytes, graminoids, litter pteridophytes, lichens, and cryptograms</td><td>categorical</td><td>forb - forb</td><td></td><td>defined</td></tr><tr><td>round</td><td>Round of sampling; round 1-2 correspond to 2018, and round 3-4 to 2019</td><td>numeric</td><td>1 - 4</td><td></td><td>defined</td></tr></tbody></table>

and Arhelleren). These measurements were made ca. 1 week after the first round of plant functional group removals in that season.

At each sampling occasion, a clear chamber (25 × 25 × 40 cm) equipped with two fans for air circulation and connected to an infrared gas analyzer (Li-840; Manufacturer reports accuracy within 1.5% of the reading value; LI-COR Biosciences, Lincoln, NE, USA) was used to measure CO₂ fluxes at all plots. To prevent cutting of roots and disruption of water flow in the plots by installing collars, we instead attached a windshield to the bottom of the chamber and weighed it down on the ground by a heavy chain to prevent wind-air mixing. At each sampling occasion we made paired measurements of fluxes under light and dark conditions, covering the chamber with a fitted light-excluding cover for the dark measurements.

NEE was estimated from measurements of CO₂ flux under ambient light and dark conditions: NEElight = GPP - Reco, NEEdark = (-) Reco. We define NEE such that negative values reflect CO₂ uptake in the ecosystem, and positive values reflect CO₂ release from the ecosystem to the atmosphere. For each measurement, CO₂ concentration was recorded at 5 s intervals over a period of 90–120 s. NEE was calculated from the temporal change of CO₂ concentration within the closed chamber according to the following formula:

$$NEE = \frac{\delta CO_2}{\delta t} \times \frac{PV}{R \times A \times (T + 273.15)}$$

where $\delta CO_2$ is the slope of the CO₂ concentration against time (µmol mol⁻¹ s⁻¹), P is the atmospheric pressure (kPa), R is the gasconstante (8.314 kPa m³ K⁻¹ mol⁻¹), T is the air temperature inside the chamber (°C), V is the chamber volume (m³) and A is the surface area (m²).

Light intensity was measured as photosynthetically active radiation (PAR, µmol m⁻² s⁻¹) using a quantum sensor (Li-190; Manufacturer reports absolute calibration accuracy of ±5%; LI-COR Biosciences, Lincoln, NE, USA) placed inside the chamber. Temperature inside the chamber was measured using an iButton temperature logger (DS1922L, Manufacturer reports temperature accuracy of ±0.5 °C, Maxim Integrated, San Jose, CA, USA). Volumetric soil moisture content (m³ water/m³ soil) × 100 was measured by calculating the average of four measurements with a soil moisture sensor (SM300, Manufacturer reports moisture accuracy of ±2.5%, Delta-T Devices, Cambridge, UK), taken at each side of a plot.

Data management and calculations. Data from the LiCOR data logger and iButton was downloaded in the field and stored. The information from the field data sheets (metadata of CO₂ measurements and plot soil moisture) was manually entered into digital worksheets, manually proof-read and stored. Data from the data logger (PAR
<table><caption>**Table 8.** Data dictionary for the FunCaB ecosystem light measurement carbon fluxes (dataset vii). Data dictionary with column descriptions for dataset vii – the ecosystem light measurement carbon fluxes from 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2017, and in 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway. Note that the dataset also contains additional measurements from another graminoid removal experiment in the VCD and the corresponding dark measurements for each plot and measurement time (Table 9).</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2015 – 2017</td><td>yyyy</td><td>defined</td></tr><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2015-06-30 – 2017-08-01</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust – Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 – Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B – Vik5RTC</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes, C = control, XC = extra control, RTC = VCG graminoid removal</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C, XC, RTC</td><td></td><td>defined</td></tr><tr><td>starttime</td><td>Start time of light measurement (UTC + 2)</td><td>date_time</td><td>2015-06-30 08:10:30 – 2017-08-01 14:17:10</td><td>yyyy-mm-dd hh:mm:ss</td><td>recorded</td></tr><tr><td>stoptime</td><td>End time of light measurement (UTC + 2)</td><td>date_time</td><td>2015-06-30 08:12:30 – 2017-08-01 14:19:10</td><td>yyyy-mm-dd hh:mm:ss</td><td>recorded</td></tr><tr><td>PAR</td><td>PAR value of light measurement</td><td>numeric</td><td>4.659 – 2555</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>recorded</td></tr><tr><td>soiltemp</td><td>Soil temperature value of light measurement</td><td>numeric</td><td>3 – 23.2</td><td>°C</td><td>recorded</td></tr><tr><td>soilmoisture</td><td>Soil moisture measurement per plot</td><td>numeric</td><td>0 – 76.775</td><td>m<sup>3</sup> water × m<sup>-3</sup> soil × 100</td><td>recorded</td></tr><tr><td>tempK</td><td>Air temperature in the chamber during light measurement</td><td>numeric</td><td>279.225 – 312.784</td><td>Kelvin</td><td>recorded</td></tr><tr><td>vegHeight</td><td>Vegetation height</td><td>numeric</td><td>0.01 – 27</td><td>mm</td><td>measured</td></tr><tr><td>nee</td><td>Net ecosystem exchange</td><td>numeric</td><td>-15.392 – 17.885</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>calculated</td></tr><tr><td>rsqd</td><td>R<sup>2</sup> squared of slope of linear regression fitting the CO<sub>2</sub> concentration versus time</td><td>numeric</td><td>0 – 1</td><td></td><td>calculated</td></tr><tr><td>chamber</td><td>Carbon flux chamber</td><td>numeric</td><td>1 – 2</td><td></td><td>defined</td></tr><tr><td>removal</td><td>Pre (2015) or post (>2015) removal treatment. Pre removal is related to the start of the experiment.</td><td>categorical</td><td>post – pre</td><td></td><td>defined</td></tr><tr><td>weather</td><td>Weather conditions during data collection</td><td>categorical</td><td>cloud – Windy</td><td></td><td>recorded</td></tr><tr><td>flag</td><td>Flag for data quality (x or DROP)</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>comment</td><td>Comment on measurement</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>time</td><td>Time in seconds for light measurement</td><td>numeric</td><td>0 – 195</td><td>second</td><td>recorded</td></tr></tbody></table>

and CO₂) and the iButton temperature logger were linked based on information from the field data sheets. All measurements were first visually evaluated for quality and only measurements that showed a consistent linear relationship between CO₂ over a time for a period of at least 60 s were used for NEE calculations. A second inclusion criterion was that this relationship had R² ≤ 0.2 or R² ≥ 0.8 for NEE measurement in light conditions and R² ≥ 0.8 for NEE dark measurements (R<sub>eco</sub>). Measurements of NEE in light conditions with R² ≤ 0.2 ensures representation of measurements with equal rates for R<sub>eco</sub> and GPP. Third, paired measurements that were more than 2 h apart were excluded. These data are available as raw fluxes and as GPP and R<sub>eco</sub> per plot per measurement (Fig. 3, Table 1).

**Dataset (vi): Reflectance.** Reflectance measures of Normalized Difference Vegetation Index (NDVI) were taken for each plot during the 2019 (post functional group removal) and 2021 (pre and post removal) field seasons (July-August), using a Trimble Greenseeker Handheld Crop Sensor (Trimble Inc., Sunnydale, CA, USA). As the sensor measures an elliptical plane, two measures perpendicular to each other were taken for each sub-plot (25 × 25 cm plot), with the centre of each ellipse being the centre of the subplot. Care was taken to ensure that sampling quadrat frames were not within the sensor range when conducting measurements (see methods Dataset ii). Measures of NDVI were taken at 60 cm above the surface where possible. Height was measured perpendicular to the sampled ground surface. These data are available as reflectance per plot per sampling date (Fig. 3, Table 1).

## Data Records

This paper reports on data from a plant functional group removal experiment replicated across the twelve sites in the Vestland Climate Grid along broad-scale bioclimatic gradients in boreal and alpine grasslands in western Norway, conducted from the 2015 growing season onwards. The datasets include the treatment variable, biomass removed, and a number of response variables reflecting microclimate, plant community composition, seedling recruitment, ecosystem carbon fluxes, and vegetation reflectance collected from 2015 (the pre-treatment year) through 2021, with variable number of years of data and temporal resolution between datasets (Fig. 1c, Table 1). Additional information and covariables reflecting site climate and environment, soils, litter decomposition, and plant functional traits in the Vestland Climate Grid can be obtained from the VCG data descriptor³⁵ and linked
<table><caption>**Table 9.** Data dictionary for the FunCaB ecosystem dark measurement carbon fluxes (dataset vii). Data dictionary with column descriptions for dataset vii – the ecosystem dark measurement carbon fluxes from 384 25 × 25 cm experimental plant functional group removal plots from 2015 – 2017, and in 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway. Note that the dataset also contains additional measurements from another graminoid removal experiment in the VCD and the corresponding light measurements for each plot and measurement time (Table 8). The light measurement data are used to calculate gpp (gross primary production) and delta (time difference between light and dark measurements).</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>year</td><td>Year of sampling</td><td>numeric</td><td>2015 - 2017</td><td>yyyy</td><td>defined</td></tr><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2015-06-30 - 2017 - 08-01</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust - Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 - Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B - Vik5RTC</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes, C = control, XC = extra control, RTC = VCG graminoid removal</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C, XC, RTC</td><td></td><td>defined</td></tr><tr><td>vegHeight</td><td>Vegetation height</td><td>numeric</td><td>0.01 - 27</td><td>mm</td><td>measured</td></tr><tr><td>gpp</td><td>Gross primary production calculated from nee - Reco</td><td>numeric</td><td>-24.342 - 0.004</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>calculated</td></tr><tr><td>starttime_Reco</td><td>Start time of dark measurement (UTC + 2)</td><td>date_time</td><td>2015-06-30 08:14:30 - 2017-08-01 14:19:55</td><td>yyyy-mm-dd hh:mm:ss</td><td>recorded</td></tr><tr><td>stoptime_Reco</td><td>End time of dark measurement (UTC + 2)</td><td>date_time</td><td>2015-06-30 08:16:00 - 2017-08-01 14:21:55</td><td>yyyy-mm-dd hh:mm:ss</td><td>recorded</td></tr><tr><td>time_Reco</td><td>Time in seconds for dark measurement</td><td>numeric</td><td>5 - 180</td><td>second</td><td>recorded</td></tr><tr><td>PAR_Reco</td><td>PAR value of dark measurement</td><td>numeric</td><td>0 - 140.131</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>recorded</td></tr><tr><td>soiltemp_Reco</td><td>Soil temperature value of dark measurement</td><td>numeric</td><td>3 - 23.2</td><td>°C</td><td>recorded</td></tr><tr><td>tempK_Reco</td><td>Air temperature in the chamber during dark measurement</td><td>numeric</td><td>279.225 - 311.977</td><td>Kelvin</td><td>recorded</td></tr><tr><td>Reco</td><td>Ecosystem respiration</td><td>numeric</td><td>0.549 - 30.354</td><td>µmol m<sup>-2</sup> s<sup>-1</sup></td><td>recorded</td></tr><tr><td>rsqd_Reco</td><td>R squared of slope of linear regression fitting the CO2 concentration versus time</td><td>numeric</td><td>0.8 - 1</td><td></td><td>calculated</td></tr><tr><td>chamber_Reco</td><td>Carbon flux chamber</td><td>numeric</td><td>1 - 3</td><td></td><td>defined</td></tr><tr><td>delta</td><td>Time difference between dark and light measurement</td><td>numeric</td><td>65 - 7050</td><td>second</td><td>calculated</td></tr><tr><td>removal</td><td>Pre (2015) or post (&gt;2015) removal treatment. Pre removal is related to the start of the experiment.</td><td>categorical</td><td>post - pre</td><td></td><td>defined</td></tr><tr><td>weather</td><td>Weather conditions during data collection</td><td>categorical</td><td>cloud - Windy</td><td></td><td>recorded</td></tr><tr><td>flag</td><td>Flag for data quality (x or DROP)</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>comment</td><td>Comment on measurement</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>time</td><td>Time in seconds for light measurement</td><td>numeric</td><td>0 - 195</td><td>second</td><td>recorded</td></tr></tbody></table>

to the FunCaB data using shared site, block, plot, treatment and taxonomic variables as keys (Fig. 3). Code to download and link selected VCG data to the FunCaB experiment is provided in the open FunCaB GitHub repository, and in a time-stamped version on Zenodo²⁸.

Data outputs consist of 6 datasets, (i) the functional group biomass removed from all treatment plots from 2015 through 2021, (ii) species-level forb biomass from the extra control plots in 2016, (iii) soil temperature in 2015 and 2016, (iv) soil moisture from 2015 through 2019, (v) species composition and structure of the experimental plots from 2015 through 2019, (vi) seedling recruitment in 2018 and 2019, (vii) ecosystem fluxes from 2015 through 2017, and (viii) reflectance in 2019 and 2021 (Table 1). The structure and relationship between datasets are described in Fig. 3. These data were checked and cleaned according to the procedures described in Fig. 2, see Technical Validation for details, before final data files and associated metadata were produced.

The final data files (see Table 1 for an overview) and all raw data, including raw flux measurements and pictures of the plots pre-and post-removals, are available at the FunCaB Open Science Framework (OSF) data repository²⁷. To ensure reproducibility and open workflows, the code necessary to access the raw data and produce cleaned datasets, along with explanations of the various data cleaning steps, issues, and outcomes, are available in the open FunCaB GitHub code repository, in the file R/data/download_clean_data.R, with a
<table><caption>**Table 10.** Data dictionary for the FunCaB vegetation reflectance (dataset viii). Data dictionary with column descriptions for dataset vii – the vegetation reflectance from 384 25 × 25 cm experimental plant functional group removal plots from 2019 and 2021, at twelve sites in the Vestland Climate Grid, Vestland County, Norway.</caption><thead><tr><th>Variable name</th><th>Description</th><th>Variable type</th><th>Variable range or levels</th><th>Units</th><th>How measured</th></tr></thead><tbody><tr><td>date</td><td>Date of sampling</td><td>date</td><td>2019-07-19 – 2021-07-30</td><td>yyyy-mm-dd</td><td>defined</td></tr><tr><td>siteID</td><td>Unique site ID</td><td>categorical</td><td>Alrust – Vikesland</td><td></td><td>defined</td></tr><tr><td>blockID</td><td>Unique block ID a combination of site and block</td><td>categorical</td><td>Alr1 – Vik5</td><td></td><td>defined</td></tr><tr><td>plotID</td><td>Unique plot ID is a combination of site, block and treatment</td><td>categorical</td><td>Alr1B – Vik5GF</td><td></td><td>defined</td></tr><tr><td>treatment</td><td>Plant functional groups removed, where F = forbs, G = graminoids, B = bryophytes and C = control, extra controls</td><td>categorical</td><td>FGB, FG, FB, GB, G, F, B, C, XC</td><td></td><td>defined</td></tr><tr><td>pre_post cut</td><td>Measurement was taken before or after the cut</td><td>categorical</td><td>post – pre</td><td></td><td>recorded</td></tr><tr><td>ndvi</td><td>NDVI value</td><td>numeric</td><td>0.315 – 4.265</td><td>percentage</td><td>measured</td></tr><tr><td>notes</td><td>Notes</td><td>categorical</td><td></td><td></td><td></td></tr><tr><td>turfID</td><td>Unique turf ID from SeedClim transplant experiment</td><td>categorical</td><td>TTC 101 – TTC 69</td><td></td><td>defined</td></tr><tr><td>weather</td><td>Weather conditions during data collection</td><td>categorical</td><td>Cloudy – Sunny</td><td></td><td>recorded</td></tr><tr><td>time</td><td>Time of sampling (UTC +2)</td><td>time</td><td>08:00 – 19:55</td><td>hh:mm</td><td>defined</td></tr></tbody></table>

time-stamped versioned copy archived in Zenodo²⁸. The reader is referred to the code and the detailed coding, data cleaning, and data accuracy comments and the associated raw and cleaned data and metadata tables available in these repositories for detailed information about the data cleaning process. Figure 2 and the Usage Notes section in this paper summarises the data accuracy and data cleaning procedures, see the latter for caveats regarding data quality and our advice on ‘best practice’ data usage.

**Dataset (i): Plant functional group removed and total biomass.** This dataset contains the biomass removed per experimental plot and functional group per year, and the total biomass in the extra controls in 2016 (Tables 1, 2, Fig. 4). For convenience, this dataset also contains basic information on the temperature and precipitation levels in the VCG climate grid (more details available in³⁵ and²⁸). From the experimental plots we removed a total of 17,839 g biomass over the project period; 5,881 g graminoids, 6,039 g forbs, and 5,918 g bryophytes. Biomass removed decreased over time, 3,742 g, 6,396 g, 1,956 g, 1,095 g, 1,701 g, 1,502 g, and 1,446 g per year from 2015 through 2021, suggesting that the removed functional groups were not able to recover from treatments and gradually were excluded from the treated plots as the experiment progressed. The lowest biomass was removed in 2018, a drought year.

The extra control plots in 2016 contained a total of 1,495 g total above-ground biomass; 357.0 g grasses, 300.0 g forbs, 322.0 g bryophytes, 7.9 g lichens, 0.9 g ferns, 1.9 g cryptogams, and 506.0 g litter. This corresponds to overall average biomass across the twelve sites of 119.0, 100.0, 107.3 and 168.7 g/m² for the three functional groups and litter, respectively. For an overview over the clean dataset see Table 2. The raw data is provided on the FunCaB OSF²⁷, the code to download and clean the data can be found in the FunCaB GitHub repository²⁸ in the file R/biomass/read_biomass.R.

**Dataset (ii): Species-level biomass in extra controls in 2016.** This dataset contains the biomass removed per plot and species for forbs in the 48 extra control plots in 2016 (Tables 1, 3). Species-level biomass was recorded for 79 forbs species, of which *Potentilla erecta* (24.3 g), *Achillea millefolium* (16.6 g), and *Alchemilla alpina* (16.3 g) were the most abundant overall. Abundances varied across sites, in the alpine *Alchemilla alpina* (12.9 g), *Alchemilla sp.* (10.0 g), and *Salix herbacea* (5.97 g) were the most abundant species, whereas in boreal sites *Achillea millefolium* (10.2 g), *Potentilla erecta* (6.88 g), and *Hieracium pilosella* (5.51 g), dominated. For an overview over the clean dataset see Table 3. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data dic/download_clean_data.R.

**Dataset (iii): Soil temperature.** The iButton dataset contains soil temperature data per plot from July 2015 to July 2016 (Tables 1, 4). This dataset has data from 322 loggers with a total of 1,084,970 observations. Soil temperatures ranged from -8.25 to 32.8 °C. During the growing season months, May through August, the average daily mean temperature was lowest in plots where only forbs were present, 9.47 ± 0.16 °C (GB treatment) and highest in the bare ground plots 10.40 ± 0.14 °C (FGB treatment). During the same period, the daily max temperatures were also lowest where only forbs were present, 17.80 ± 0.16 °C (GB treatment) and highest in the bare ground plots 20.40 ± 0.14 °C (FGB treatment). For an overview over the clean dataset see Table 4. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data dic/download_clean_data.R.

**Dataset (iv): Soil moisture.** The soil moisture dataset contains volumetric soil moisture data per plot from June 2015 to August 2018 (Tables 1, 5). The database contains 4,771 soil moisture measurements from
<figure><img src="image_6.png" /><figcaption>**Fig. 4** Biomass over time in the FunCaB plant functional group removal experiment. Biomass removed per site, treatment, and plant functional group over time from 384 25 × 25 cm experimental plant functional group removal plots from 2015–2021, and total biomass from 48 extra controls in 2016, at twelve sites in the Vestland Climate Grid, Vestland County, Norway. Sites are coded by vegetation zone (boreal, sub-alpine, alpine) and precipitation level (low = 1 to high = 4). B = bryophytes, F = forbs, G = graminoids. Removal treatments are coded by the plant functional group removed (i.e., in the B treatments, F and G are not removed).</figcaption></figure>

2015 to 2019, with the median number of measurements per plot per year ranging from 3 – 5, except in 2018 where the median was one measurement per plot. Average soil moisture was 20.2 ± 0.4% in 2015, 26.1 ± 0.5% in 2016, 29.6 ± 0.6% in 2017, 19.0 ± 0.9% in 2018, and 13.7 ± 0.5% in 2019, which was a drought year. Plots with forbs only and all functional groups removed were the driest, with an average soil moisture of 21.8 ± 0.7% (F) and 22.1 ± 0.7 (FGB) whereas the intact control plots were the wettest with 24.1 ± 0.6% (C). For an overview over the clean dataset see Table 5. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data//download_clean_data.R.

**Dataset (v): Vascular plant community composition and structure.** The plot-level plant community dataset contains a total of 206 taxa and 17,181 observations (taxa × plots × years) (Tables 2, 6). Overall, 299 observations (1.74%) were not identified to taxon. In 2015, the pre-treatment vegetation had an average species richness of 5.9 graminoids (ranging from 3.2 to 9.3 across sites) and 10.8 forbs (ranging from 3.8 to 14.9 across sites) per plot. The highest species richness was found in alpine sites. Sum of covers of the non-removed functional groups varied over sites and years (Fig. 5). In many cases, biomass of the remaining functional group increases over time in plots where other functional groups were removed, relative to controls, suggesting release from competitive interactions may be occurring. However, in some cases, notably forbs in colder and wetter sites, responses are neutral to negative, which is consistent with earlier work in the system suggesting that competitive plant-plant interactions are less important and facilitative interactions more important under these conditions¹². Similarly, there were generally lower abundances and also lower relative biomass gains in removals in 2018, which was a drought year, suggesting competitive release was less important under this climate extreme event (Fig. 5). Note that some of the recorded data (e.g., sub-plot level presence data, life-history stage/fertility data, acrocarp/pleurocarp bryophyte cover data) are available in the raw data tables, but these are not processed or presented in the final data table. For an overview over the clean dataset see Table 6. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data and the taxon table can be found in the FunCaB GitHub repository²⁸ in the file R/data//download_clean_data.R.

**Dataset (vi): Seedling recruitment.** The seedling recruitment dataset has data on a total of 16,656 seedlings of 44 species (Tables 2, 7). Overall, 61.2% of the seedlings were identified to species, the remaining 38.8% were either unidentifiable for the duration of the study or they died before they could be identified. Seedling densities were highest in the bare soil and forb only plots (FGB and GB treatments) with 1,297 and 1,451 seedlings m⁻², respectively, and lowest in the graminoid and bryophyte and graminoid only plots (F and FB treatments) with 228 and 224 seedlings m⁻², respectively. The intact vegetation controls (C) had intermediate seedling densities, with 764 seedlings m⁻². Seedling emergence differed between the two study years, with 5,888 seedlings recorded in the drought year 2018 compared to 10,768 in 2019. Overall, 5,870 seedlings survived till the last census. For an overview over the clean dataset see Table 7. The final cleaned data together with a data dictionary
is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data下载_clean_data.R.

**Dataset (vii)- Ecosystem carbon fluxes.** The plot-level ecosystem carbon flux dataset has a total of 3,696 paired light (GPP) and dark (Reco) measurements between 2015 and 2017 (Tables 1, 8, 9). A total of 926 paired fluxes were recorded in 2015, 1,286 in 2016 and 1,484 in 2017. The median number of measurements per plot and year in 2015, 8 in 2016 (when we measured control plots only) and 5 in 2017. The GPP values ranged from -24.30 ± 0.06 to -0.01 ± 0.06 µmol mol⁻¹ s⁻¹ CO₂, the PAR between 4.66 ± 8.46 and 2555.00 ± 8.46 (µmol m⁻² s⁻¹), the soil moisture between 0 ± 0.28 and 76.78 ± 0.28 and the air temperature within the chamber varied between 6.08 ± 0.10 and 39.63 ± 0.10 °C. The respiration measurements (Rₑco) ranged from 0.55 ± 0.05 to 30.35 ± 0.05 µmol mol⁻¹ s⁻¹ CO₂, the soil moisture between 0 ± 0.28 and 76.78 ± 0.28 and the associated air temperature varied within the chamber measurements between 6.08 ± 0.10 and 38.83 ± 0.10 °C. Note that light and dark measurements have very similar soil moisture and air temperature, as would be expected given that they are both measured on the same plots within a short period. Average GPP varied among treatments, and on average, the highest CO₂ uptake was found in intact vegetation (C; -8.28 ± 0.19 µmol mol⁻¹ s⁻¹ CO₂) while the lowest uptake was found in bare ground plots (FGB treatment; -4.95 ± 0.27 µmol mol⁻¹ s⁻¹ CO₂). Respiration did not vary much across treatments, but for GPP, the highest value 6.63 ± 0.15 µmol mol⁻¹ s⁻¹ CO₂ was again found in the intact vegetation (C) and the lowest value 5.04 ± 0.17 µmol mol⁻¹ s⁻¹ CO₂ in the bare ground plots (FGB treatment). For an overview over the clean dataset see Tables 8, 9. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data dic/download_clean_data.R.

**Dataset (xi) - Reflectance.** The plot-level community reflectance dataset has a total of 1,244 observations. In 2019 we made 480 post-treatment measurements, and in 2021 we made 382 pre-treatment measurements and 382 post-treatment measurements (Tables 2, 10). Average reflectance post-treatment was 0.634 ± 0.004 (mean ± SE) in 2019 and 0.679 ± 0.006 in 2021, whereas the reflectance pre-treatment in 2021 was 0.760 ± 0.010. In both years, the highest reflectance post cut was found in the intact vegetation (C), 0.695 ± 0.007 in 2019 and 0.774 ± 0.009 in 2021, whereas the lowest reflectance was found in the bare ground treatment (FGB) 0.524 ± 0.008 in 2019 and 0.507 ± 0.015 in 2021. For an overview over the clean dataset see Table 10. The final cleaned data together with a data dictionary is provided on the FunCaB OSF²⁷, and the code to download the final cleaned data can be found in the FunCaB GitHub repository²⁸ in the file R/data dic/download_clean_data.R.

**Technical Validation**

**Experimental validation.** The plant functional group experiments and field measurements (before-after treatment photographs, soil moisture measurements, site-level VCD climate data) were maintained over seven years under the oversight of LCK, supported by IHJA and SAHÖ in some years and field campaigns. Consistent leadership of the campaigns ensured consistency in experimental treatments and data collection, which was critical as a large number of students, technicians and interns were involved in the plant functional group removal work, typically 4 - 6 people each year (see Acknowledgements). Errors made during field work are noted in the data sheets and entered as comments in the plant functional group removal biomass data tables. The photographs of each plot pre- and post-removal for each field campaign (with some gaps, and photos are lacking from 2020) can be used for further documentation and validation²⁷. Biomass removed decreased over time, especially after the second year (Fig. 5, suggesting the communities stabilised in the experimentally defined functional group composition, with relatively little regrowth of the removed functional group after this point.

**Taxonomic validation.** During the 7-year data collection period a number of people were involved in community and seedling data collection (see notes in dataset for recorder information), which introduces a risk of observer errors. In particular, species can be misidentified (e.g., sterile graminoids) or might be overlooked in one of the community censuses. Such errors, if uncorrected, generate pseudo-turnover in the plant community and seedling data. To detect and correct errors in the community data, we compared each recorded species in each plot or subplot over the time-series. We used these comparisons to assign unidentified or missing species to an existing identified species in that plot if it was likely that this could be the correct ID, and if there was a record of that species in the plot in the previous and following year. Further, we imputed species covers in cases where cover was either missing or clearly too low or high to be realistic when comparing with the total sum of covers and covers of the same species in adjacent years in the time-series, replacing such erroneous values with the mean cover from the previous and following year. We further checked and corrected 'botanist effects' in cover estimation between observers, in height estimations (using cm vs mm), and ensured taxonomic consistency. The data-checking code and outcomes for these various procedures is documented in²⁸. In addition, 23 forbs and six graminoid taxa are identified to the genus level in the community dataset and one graminoid, two forbs and seven seedling observations are unidentified (these are named NID.herb, NID.gram, and NID.seedling in the datasets) (datasets v, vi). A full species list of all identified taxa across datasets is also available in the taxon table in the OSF repository (see Table 2). This table also includes the functional group (F, G, or B) to which each taxon was assigned in the experimental treatments. Note that the sub-plot-level presences and life-history stage/fertility data are available in the OSF, but these data have not been processed or cleaned.

**Ecosystem C flux validation.** All measurements were first visually evaluated for data quality. If there were clear outliers (values far outside the normal range of CO₂ concentrations within a measurement), these values were flagged and removed before flux calculation. The timeframes of measurements were adjusted if the linear
<figure><img src="image_7.png" /></figure>

**Fig. 5** Functional group cover in response to removal treatments in the FunCaB plant functional group experiment. Difference from unmanipulated control plots (C treatment) over time in the sum of covers of graminoids, forbs, and bryophytes in plots where this functional group grows alone (i.e., the FB, GB, and GF removal treatments, respectively, see Table 6 for abbreviations). Data from 192 25 × 25 cm experimental plant functional group removal plots from 2015 – 2019 at twelve sites in the Vestland Climate Grid, Vestland County, Norway.

relationship over time was affected by either change in cloud cover or there was indication of chamber leakage (CO₂ increase/decrease levelling off after steady linear relationship). Resulting fluxes were only included if they showed a consistent linear relationship between CO₂ and time for a period of at least 60 seconds. Measurements of poor quality due to erratic CO₂ concentrations over time were flagged and removed. An additional criterion was that this relationship had R² ≤ 0.2 or R² ≥ 0.8 for NEE measurement in light conditions and R² ≥ 0.8 for NEE dark measurements (R eco; 52, 150 and 466 measurements removed in 2015, 2016 and 2017, respectively). Measurements of NEE in light conditions with R² ≤ 0.2 ensure inclusion of measurements with equal rates for R eco and GPP. Dark measurements with resulting NEE < 0 were removed, as they would indicate uptake CO₂ uptake (49 measurements removed in 2017). Sets of L and D measurement per plot were then matched together based on date and time of day for calculation of GPP but were not paired if they were more than 2 hours apart from each other. Finally, GPP fluxes with values > 0 were removed, as they would represent CO₂ release, not uptake/primary productivity (0 measurements were thus removed in 2015, 18 in 2016, and 191 in 2017). The PAR sensor for one of the chambers was malfunctioning between 03 July 2017 and 6 July 2017. PAR values for these measurements were corrected by using an estimated value based on recorded PAR values in the metadata from the other par sensor.

**Usage Notes**

**Data use best practice.** We suggest that all FunCaB project data presented here and available in the FunCaB OSF²⁷ be cited to this data paper. Other VCG data available through the VCG OSF³⁵ should be cited to the relevant publications, specifically, the general description of the climate and environment, soils, land-use and environment, vegetation, and ecosystems along with basic descriptive data of the 12 sites are presented in³⁴, whereas specific biotic and abiotic response data are presented in e.g.³⁷⁻⁴⁸. Note that major VCG datasets will be made
available in a forthcoming data paper, and they can henceforth be cited to that paper. See Fig. 3 for a description of the interrelations and linkage keys among datasets.

The data are available for use under a CC-BY licence. We appreciate being contacted for advice or collaboration, if relevant, by users of our data. In cases where our data make up >10% of the data used in a downstream publication, we expect that appropriately acknowledging our contributions would require an invitation for collaboration.

**Combination and augmentation with existing and new data from the VCG.** We note that for plant functional group removal biomass (dataset i), soil moisture (dataset iv), plant community (dataset v), carbon flux (dataset vii) and reflectance (dataset viii), a final census planned for 2022. These data will be added to²⁷ and²⁸.

The data presented here relates to a large amount of site and block level data from the same sites³⁵, the most relevant of which are briefly explained under 'Background abiotic and biotic data from the Vestland Climate Grid' above. Code is provided for linking and downloading the most relevant datasets in the FunCaB github repository, with a time-stamper version on Zenodo²⁸. See Fig. 2 for a conceptual representation of how these datasets are linked via shared variables/keys.

In connection with a new project, FUNDER (Norwegian Research Council project number 315249, 2021 – 2025), the FunCaB experiment will be augmented with data on the belowground components of the plant-soil ecosystem, including roots, mesofauna, fungi and microbes. Plot-level slope, aspect, and soil depth will also be measured. These data will link with the FunCaB and VCG project based on the given experimental, site and organismal keys, as indicated in Fig. 3.

The VCG and FUNDER data will be made available in forthcoming data papers.

**Data quality comments and options.** The procedures and consequences for data collection, management, and cleaning are detailed in this paper, and the associated code²⁸. The code²⁸ describes and implements our suggested data cleaning and checking procedures that result in producing what we consider the clean and 'best practice' final datasets²⁷ and the various 'flag', 'comment' or 'notes' columns in the different dataset tables (Tables 2–9) that indicate additional specific data points that could be removed to create even more robust datasets. Users who might prefer stricter or more inclusive data handling strategies should check the flags in the raw data sets and adjust the data cleaning accordingly.

Note that unidentified taxa are not harmonized across datasets (relevant for datasets ii, v, vi; Table 1).

Also note that dataset ii species-level forb biomass has some blocks where data at the species level is lacking (11 cases, out of 48). For these blocks, the forb biomass is provided at the functional group level.

Dataset vi, seedling recruitment, does not contain any information on the plots that did not have any seedlings. If the user wants a complete dataset with all plots per site, for example for calculating overall seedling densities per treatment or site, these missing plots need to be extracted from another dataset, e.g. dataset v, plant community composition.

For dataset vii, ecosystem carbon fluxes, users might want to set more or less restrictive data exclusion thresholds regarding linear fits of fluxes to include in analysis, and can do so by altering the cleaning code. For example, many paired light-dark measurements were deleted due to the quality control of the light measurements. Users who are interested in dark measurements only might be interested in accessing the full dark measurements dataset available from the raw data available in²⁷. If users want to visually inspect individual measurements and set their own timeframes the code for running through individual combinations of metadata, logger data and temperature data are provided in the OSF. Note that in 2016, some control plots shared with another experiment (treatment = TTC) and other graminoid removal plots (treatment = RTC) from another project³⁵, see also⁴², were also measured. These measurements are found in the FunCaB flux dataset.

**Code availability**

The code used for checking, cleaning and analysing the data is available in the open GitHub repository https://github.com/Between-the-Fjords/funcab_data, of which a versioned copy is available at Zenodo⁵⁴.

Received: 21 December 2021; Accepted: 12 July 2022;

Published online: 28 July 2022

**References**

1. Steinbauer, M. J. et al. Accelerated increase in plant species richness on mountain summits is linked to warming. *Nature* **556**, 231–234 (2018).

2. Walther, G.-R. et al. Ecological responses to recent climate change. *Nature* **416**, 389–395 (2002).

3. Elmendorf, S. C. et al. Experiment, monitoring, and gradient methods used to infer climate change effects on plant communities yield consistent patterns. *Proc. Natl. Acad. Sci. USA* **112**, 448–452 (2015).

4. Bardgett, R. D. & Wardle, D. A. *Aboveground-belowground linkages: biotic interactions, ecosystem processes, and global change*. (Oxford University Press, 2010).

5. Tylianakis, J. M., Didham, R. K., Bascompte, J. & Wardle, D. A. Global change and species interactions in terrestrial ecosystems. *Ecol. Lett.* **11**, 1351–1363 (2008).

6. Berg, M. P. et al. Adapt or disperse: understanding species persistence in a changing world. *Global Change Biology* **16**, 587–598 (2010).

7. Roslin, T. et al. Higher predation risk for insect prey at low latitudes and elevations. *Science* **356**, 742–744 (2017).

8. Alexander, J. M., Diez, J. M. & Levine, J. M. Novel competitors shape species' responses to climate change. *Nature* **525**, 515–518 (2015).

9. Scafaro, A. P. et al. Responses of leaf respiration to heatwaves. *Plant Cell Environ.* **44**, 2090–2101 (2021).

10. Lin, Y.-S., Medlyn, B. E. & Ellsworth, D. S. Temperature responses of leaf net photosynthesis: the role of component processes. *Tree Physiol.* **32**, 219–231 (2012).
11. Wu, Z., Dijkstra, P., Koch, G. W., Peñuelas, J. & Hungate, B. A. Responses of terrestrial ecosystems to temperature and precipitation change: a meta-analysis of experimental manipulation. *Glob. Chang. Biol.* **17**, 927–942 (2011).

12. Sundqvist, M. K. et al. Responses of tundra plant community carbon flux to experimental warming, dominant species removal and elevation. *Funct. Ecol.* **34**, 1497–1506 (2020).

13. Fanin, N. et al. Effects of plant functional group removal on structure and function of soil communities across contrasting ecosystems. *Ecology Letters Volume 22*, 1095–1103 (2019).

14. Bar-On, Y. M., Phillips, R. & Milo, R. The biomass distribution on Earth. *Proc. Natl. Acad. Sci. USA* **115**, 6506–6511 (2018).

15. Beerling, D. *The Emerald Planet: How Plants Changed Earth's History*. (Oxford University Press, 2017).

16. Bjorkman, A. D. et al. Plant functional trait change across a warming tundra biome. *Nature* **562**, 57–62 (2018).

17. Elmendorf, S. C. et al. Global assessment of experimental climate warming on tundra vegetation: heterogeneity over space and time. *Ecol. Lett.* **15**, 164–175 (2012).

18. Elmendorf, S. C. et al. Plot-scale evidence of tundra vegetation change and links to recent summer warming. *Nat. Clim. Chang.* **2**, 453–457 (2012).

19. Ross, L. C., Woodin, S. J., Hester, A. J., Thompson, D. B. A. & Birks, H. J. B. Biotic homogenization of upland vegetation: patterns and drivers at multiple spatial scales over five decades. *Journal of Vegetation Science* **23**, 755–770 (2010).

20. Klanderud, K. & Totland, Ø. Simulated climate change altered dominance hierarchies and diversity of an alpine biodiversity hotspot. *Ecology* **86**, 2047–2054 (2005).

21. Michel, P., Burritt, D. J. & Lee, W. G. Bryophytes display allelopathic interactions with tree species in native forest ecosystems. *Oikos* **120**, 1272–1280 (2011).

22. Soudzilovskaja, N. A. et al. How do bryophytes govern generative recruitment of vascular plants? *New Phytol.* **190**, 1019–1031 (2011).

23. Wardle, D. A. et al. Linking vegetation change, carbon sequestration and biodiversity: insights from island ecosystems in a long-term natural experiment. *J. Ecol.* **100**, 16–30 (2012).

24. Guido, A. & Pillar, V. D. Are removal experiments effective tools for assessing plant community resistance and recovery from invasion? *J. Veg. Sci.* **26**, 608–613 (2015).

25. Artsdatabanken. Norsk rodliste for naturtyper. https://www.artsdatabanken.no/rodlistefornaturtyper (2018).

26. Bokhorst, S. et al. Variable temperature effects of Open Top Chambers at polar and alpine sites explained by irradiance and snow depth. *Global Change Biology* **19**, 64–74 (2013).

27. Vandvik, V. & the FunCaB Project Team. FunCaB - The role of functional group interactions in mediating climate change impacts on the Carbon dynamics and biodiversity of alpine ecosystems, *Open Science Framework*, https://doi.org/10.17605/OSF.IO/4C5V2 (2022).

28. Halbritter A. H. and the FunCaB Data Team. GitHub repository for FunCaB - The role of functional group interactions in mediating climate change impacts on the Carbon dynamics and biodiversity of alpine ecosystems. Version 2.1. Available at: github.com.

29. Halbritter, A. H. et al. The handbook for standardized field and laboratory measurements in terrestrial climate change experiments and observational studies (ClimEx). *Methods Ecol. Evol.* **11**, 22–37 (2019).

30. Wilkinson, M. D. et al. The FAIR Guiding Principles for scientific data management and stewardship. *Sci Data* **3**, 160018 (2016).

31. Alston, J. M. & Rick, J. A. A beginner's guide to conducting reproducible research. *Bull. Ecol. Soc. Am.* **102**(2), e01801 (2021).

32. Hampton, S. E. et al. The Tao of open science for ecology. *Ecosphere* **6**, art120 (2015).

33. Tveito, O. E., Bjordal, I., Skjelvåg, A. O. & Aune, B. A GIS-based agro-ecological decision system based on gridded climatology. *Meteorol. Appl.* **12**, 57–68 (2005).

34. Klanderud, K., Vandvik, V. & Goldberg, D. The Importance of Biotic vs. Abiotic Drivers of Local Plant Community Composition Along Regional Bioclimatic Gradients. *PLoS One* **10**, e0130205 (2015).

35. Vandvik, V. & The Vestland Climate Grid and SeedClim Project Teams. Vestland Climate Grid (VCG). *Open Science Framework* https://doi.org/10.17605/OSF.IO/NPFA9 (2022).

36. Fremstad, E. *Vegetasjonstyper i Norge*. (NINA, Norsk Institutt for Naturforskning, 1998).

37. Klanderud, K. et al. Vital rates in early life history underlie shifts in biotic interactions along bioclimatic gradients: An experimental test of the Stress Gradient Hypothesis. *J. Veg. Sci.* **32**, e13006 (2021).

38. Vandvik, V., Klanderud, K., Meineri, E., Mären, I. E. & Töpper, J. Seed banks are biodiversity reservoirs: species-area relationships above versus below ground. *Oikos* **125**, 218–228 (2016).

39. Lynn, J. S., Klanderud, K., Telford, R. J., Goldberg, D. & Vandvik, V. Macroecological context predicts species' responses to climate warming. *Glob. Chang. Biol.* **27**, 2088–2101 (2021).

40. Vandvik, V. et al. Biotic rescaling reveals importance of species interactions for variation in biodiversity responses to climate change. *Proc. Natl. Acad. Sci. USA* (2020).

41. Guittar, J. et al. Quantifying the roles of seed dispersal, filtering, and climate on regional patterns of grassland biodiversity. *Ecology* e03061 (2020).

42. Olsen, S. L., Töpper, J. P. & Skarpaas, O. From facilitation to competition: temperature-driven shift in dominant plant interactions affects population dynamics in seminatural grasslands. *Glob. Chang. Biol.* **22**, 1915–1926 (2016).

43. Klanderud, K., Meineri, E., Töpper, J., Michel, P. & Vandvik, V. Biotic interaction effects on seedling recruitment along bioclimatic gradients: testing the stress-gradient hypothesis. *J. Veg. Sci.* **28**, 347–356 (2017).

44. Töpper, J. P. et al. The devil is in the detail: Nonadditive and context-dependent plant population responses to increasing temperature and precipitation. *Glob. Chang. Biol.* **24**, 4657–4666 (2018).

45. Meineri, E., Klanderud, K., Guittar, J., Goldberg, D. E. & Vandvik, V. Functional traits, not productivity, predict alpine plant community openness to seedling recruitment under climatic warming. *Oikos* **129**, 13–23 (2019).

46. Keuskamp, J. A., Dingemans, B. J. J., Lehtinen, T., Sarneel, J. M. & Hefting, M. M. Tea Bag Index: a novel approach to collect uniform decomposition data across ecosystems. *Methods Ecol. Evol.* **4**, 1070–1075 (2013).

47. Althuizen, I. H. J., Lee, H., Sarneel, J. M. & Vandvik, V. Long-Term Climate Regime Modulates the Impact of Short-Term Climate Variability on Decomposition in Alpine Grassland Soils. *Ecosystems* **21**, 1580–1592 (2018).

48. Skarpaas, O., Meineri, E., Bargmann, T. & Pötsch, C. Biomass partitioning in grassland plants along independent gradients in temperature and precipitation. *Perspect. Plant Ecol. Evol. Syst.* **19**, 1–11 (2016).

49. Meineri, E., Skarpaas, O. & Spindelböck, J. Direct and size-dependent effects of climate on flowering performance in alpine and lowland herbaceous species. *J. Veg. Sci.* **25**, 275–286 (2014).

50. Meineri, E., Spindelböck, J. & Vandvik, V. Seedling emergence responds to both seed source and recruitment site climates: a climate change experiment combining transplant and gradient approaches. *Plant Ecol.* **214**, 607–619 (2013).

51. Pérez-Harguindeguy, N. et al. New handbook for standardised measurement of plant functional traits worldwide. *Aust. J. Bot.* **61**, 167–234 (2013).

52. Gya, R. The role of intraspecific variability in driving community trait shifts along temperature and precipitation gradients in alpine and boreal semi-natural grasslands. (The University of Bergen, 2017).

53. Lid, J. & Lid, D. T. Norsk flora (Norwegian flora). Oslo: Samlaget (2005).

54. Halbritter, A. H. & The FunCaB Project Team. Between-the-Fjords/funcab_data: Publish FunCaB data paper, Zenodo https://doi.org/10.5281/zenodo.6752847 (2022).

55. CRediT - Contributor Roles Taxonomy. https://casrai.org/credit/ (2019).
## Acknowledgements

This research was conducted in the 12 sites of the Vestland Climate Grid, and we would like to thank the landowners for permission to established and maintain the experiments on their lands. The **FunCaB** project 'The role of Functional group interactions in mediating climate change impacts on the Carbon dynamics and Biodiversity of alpine ecosystems' was funded by the Research Council of Norway KLIMAFORSK program, project number 244525 (2015 - 2019). The research sites and VCG data collection has been funded by the Research Council of Norway under the projects **SEEDCLIM**: NORKLIMA project number 184912 (2008 - 2015), **INCLINE** FRIMEDBIO project number 274712 (2018 - 2022), and **EMERALD**: KLIMAFORSK project number 294948 (2019 - 2022). Student internships were funded under the Research Council of Norway projects **RECITE INTPART** project number 274831 (2018 - 2021), **ExperTS**: INTPART project number 287784 (2019 - 2022), and the Centre for Internationalisation in Education (SIU) project **TraitTrain**: SIU project HNP-2015/10037 (2016 - 2018). We acknowledge a covid-19 infection (v. Delta) affecting VV and thereby providing the isolated time necessary to write this manuscript, in lieu of the otherwise required long-distance flight. The following people assisted with field and lab work: Terezie Novakova, Vojtěch Škara, Tabea Galusser, Alba Tousfandos, Linn Vassvik, Kari Elva Aarvik, Dan from America, Amund Skogrand, Maja Ugrin, Emilia Meneses-Reyes, Aynhoa Gomez, Arild Breistøl, Olav V. Skogrand, Ingrid Toresen, Ardian Høgøy Abaz, Nataniel Nesse, Frida Klubb, Kristina Svare Toby Santamaria, Peter Petrik, David Escobar Cruz, Anja Petek, Loushu He, Elisabeth Bergman, William Erazo Garcia, Jesslyn Tjendra, Lucely Vilca Bustamante, Sára Nagy, Kengo Mieda, Betsabe Kuan, Daiki Koga, Vincent Zimmermann, Philippine Surer, Eira Rønneberg, Alexander Vågenes, Leire Campos, Emma Little, Maryam Kazemi, Morgane Demeaux, Regine Gulbrandsen, Eliska Suchoparova, Julie Eriksen, and Jana Rüthers.

## Author contributions

We follow the CreDiT taxonomy⁵⁵, and recognize the following author contributions, Conceptualization (co), Data curation (da), Formal analysis (fo), Funding acquisition (fu), Investigation (in), Methodology (me), Project administration (pr), Resources (re), Software (so), Supervision (su), Validation (va), Visualization (vi), Writing - original draft (wo), and Writing - review & editing (wr), as follows: V.V. (co, fu, in, me, pr, re, su, wo), F.J. and I.H.J.A. (da, fo, in, me, va, vi, wr), D.E.G., H.L., K.K. and S.L.O. (co, fu, me, re, su, wr), R.J.T. (da, fo, me, re, su, va, wr), L.C.K. (in, pr, va, wr), S.R.G. (da, in, me, wr), D.D.E. and S.A.H. Ö. (in, pr, wr), E.M., I.D., J.S.L., R.G., S.B. and S.Y. (in, wr), and A.H.H. (da, fo, me, pr, re, su, va, vi, wo).

## Funding

Open access funding provided by University of Bergen.

## Competing interests

The authors declare no competing interests.

## Additional information

**Correspondence** and requests for materials should be addressed to V.V.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_8.png" /></figure>

**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 license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license 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 license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022