en: grassland soil topsoil-removal plants bacteria nematodes biodiversity soil-biology soil-chemistry soil-biodiversity restoration fungi

en: Networks

en:
Treatment	Treatment groups; Y.4-18 = restored grasslands; Initial = intensively managed grasslands; Target = species-rich semi-wet to semi-dry grasslands
Correlation	non-parametric Spearman rank correlations
Vertex.1+2	Comparison between network components; N.pool = total soil nitrogen pool [kg ha-1 12 cm-1]; Corg.pool = soil organic carbon pool [kg ha-1 12 cm-1], pH = soil pH [CaCl2]; BD = bulk density [g cm-3]; FC = field capacity [%]
	PCoA 1 + 2 biotic communities [Veg = vegetation; Nem = nematoda; Fun = fungi; Prok = prokaryota]
Type	Interactions type; A.A = abiotic.abiotic; A.B = abiotic.biotic; B.B = biotic.biotic

en: Fungi data

en:
Sample.data
Plot.ID	consecutive numbers for plots (in total 36 plots)
Treatment	One restoration treatment and two reference treatments. Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Initial=intensively managed grasslands; Target= species-rich semi-wet to semi-dry grasslands
treat.yr	Treatment abbreviation. For restored grasslands: Treatment abbreviation + time since restoration; Initial + Target = NA
yr.trt	year when topsoil was removed; Initial + Target = NA
Age	time since restoration (sampling year 2018 - yr.trt); Initial + Target = NA
Age.class	restored sites of similar age grouped into four age classes. Y.4 = Year 3+4 (recently restored); Y.18 = Year 17+19 (mid-term); Y.24 = Year 23+25 (long-term); Y.30 = Year 27+32 (long-term); Initial + Target = NA
Cluster	clustered spatial arrangment of sites in the field (see Supporting Information Fig. S1)
Site	clustered spatial arrangment (replicates) of plots per site (see Supporting Information Fig. S1)

Taxonomy
OTU.ID	consecutive numbers for operational taxonomic units (OTUe 1-7794)
Taxonomic rank	Domain, Phylum, Class, Order, Family, Genus, Species

Abundance
columns	Plot.ID (see Sample.data)
rows	OTU.ID (see Taxonomy)

en: Ecosystem coupling_interactions

en:
Treatment	Treatment groups; Y.4-18 = restored grasslands; Initial = intensively managed grasslands; Target = species-rich semi-wet to semi-dry grasslands
Type	Interactions type; A.A = abiotic.abiotic; A.B = abiotic.biotic; B.B = biotic.biotic
Correlation	non-parametric Spearman rank correlations between biotic and abiotic parameters (ecological interactions)

en: Distance to target results

en:
Community	biotic community (vegetation, nematoda, fungi, prokaryota)
Treatment	Treatment groups; Y.4-30= restored grasslands (time since restoration); Initial=intensively managed grasslands (I)
Comparison	within- or between-community Bray-Curtis dissimilarity distance comparison
Label	within-community distances of initial/restored grasslands: treatment abbreviation; between-communities distances of initial/restored-to-target grasslands: treatment abbreviation plus Target (T; species-rich grasslands)
mean	Mean distance to Target (6 plots per treatment group)
sd	Standard deviation of distances to Target (6 plots per treatment group)

en: Vegetation data

en:
Plot.ID	consecutive numbers for plots (in total 36 plots)
Treatment	One restoration treatment and two reference treatments. Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Initial=intensively managed grasslands; Target= species-rich semi-wet to semi-dry grasslands
treat.yr	Treatment abbreviation. For restored grasslands: Treatment abbreviation + time since restoration; Initial + Target = NA
yr.trt	year when topsoil was removed; Initial + Target = NA
Age	time since restoration (sampling year 2018 - yr.trt); Initial + Target = NA
Age.class	restored sites of similar age grouped into four age classes. Y.4 = Year 3+4 (recently restored); Y.18 = Year 17+19 (mid-term); Y.24 = Year 23+25 (long-term); Y.30 = Year 27+32 (long-term); Initial + Target = NA
Cluster	clustered spatial arrangment of sites in the field (see Supporting Information Fig. S1)
Site	clustered spatial arrangment (replicates) of plots per site (see Supporting Information Fig. S1)
column I-LW	Plant abundance per species (in total 327 plant species) expressed as cover percentage estimation according to van der Maarel (vdM 1979; 7 classes: 1-7); nomenclature according to Lauber & Wagner (1996)

en: Soil properties

en:
Plot.ID	consecutive numbers for plots (in total 36 plots)
Treatment	One restoration treatment and two reference treatments. Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Initial=intensively managed grasslands; Target= species-rich semi-wet to semi-dry grasslands
treat.yr	Treatment abbreviation. For restored grasslands: Treatment abbreviation + time since restoration; Initial + Target = NA
yr.trt	year when topsoil was removed; Initial + Target = NA
Age	time since restoration (sampling year 2018 - yr.trt); Initial + Target = NA
Age.class	restored sites of similar age grouped into four age classes. Y.4 = Year 3+4 (recently restored); Y.18 = Year 17+19 (mid-term); Y.24 = Year 23+25 (long-term); Y.30 = Year 27+32 (long-term); Initial + Target = NA
Cluster	clustered spatial arrangment of sites in the field (see Supporting Information Fig. S1)
Site	clustered spatial arrangment (replicates) of plots per site (see Supporting Information Fig. S1)
N.pool	Soil nitrogen pool (June 2017); corrected for bulk density and 12 cm sampling depth [kg ha-1 12cm-1]
Corg.pool	Soil organic carbon pool (June 2017); corrected for bulk density and 12 cm sampling depth [kg ha-1 12cm-1]
pH	soil pH [CaCl2]
BD	Bulk density [g cm-3]
FC	Field capacity [%]
Veg.1	principle coordinate axis 1 (vegetation)
Veg.1	principle coordinate axis 2 (vegetation)
Nem.1	principle coordinate axis 1 (nematoda)
Nem.2	principle coordinate axis 2 (nematoda)
Fun.1	principle coordinate axis 1 (fungi)
Fun.2	principle coordinate axis 2 (fungi)
Prok.1	principle coordinate axis 1 (prokaryota)
Prok.2	principle coordinate axis 2 (prokaryota)
Veg.H	Shannon diversity (vegetation)
Nem.H	Shannon diversity (nematoda)
Fun.H	Shannon diversity (fungi)
Prok.H	Shannon diversity (prokaryota)

en: Nematode data

en:
Sample.data
Plot.ID	consecutive numbers for plots (in total 36 plots)
Treatment	One restoration treatment and two reference treatments. Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Initial=intensively managed grasslands; Target= species-rich semi-wet to semi-dry grasslands
treat.yr	Treatment abbreviation. For restored grasslands: Treatment abbreviation + time since restoration; Initial + Target = NA
yr.trt	year when topsoil was removed; Initial + Target = NA
Age	time since restoration (sampling year 2018 - yr.trt); Initial + Target = NA
Age.class	restored sites of similar age grouped into four age classes. Y.4 = Year 3+4 (recently restored); Y.18 = Year 17+19 (mid-term); Y.24 = Year 23+25 (long-term); Y.30 = Year 27+32 (long-term); Initial + Target = NA
Cluster	clustered spatial arrangment of sites in the field (see Supporting Information Fig. S1)
Site	clustered spatial arrangment (replicates) of plots per site (see Supporting Information Fig. S1)

Taxonomy
ASV.ID	unit number.letter code per amplicon sequence variant (in total 904)
Taxonomic rank	Phylum, Class, Order, Family, Genus, Species

Abundance
columns	Plot.ID (see Sample.data)
rows	ASV.ID (see Taxonomy)

en:
This dataset contains all data, on which the following publication below is based.

Paper Citation:

_Resch, M.C., Schütz, M., Ochoa-Hueso, R., Buchmann, N., Frey, B., Graf, U., van der Putten, W.H., Zimmermann, S., Risch, A.C. (in review). Long-term recovery of above- and belowground interactions in restored grassland after topsoil removal and seed addition. Journal of Applied Ecology_

Please cite this paper together with the citation for the datafile.

Study area and experimental design
The study was conducted in and around two nature reserves, Eigental and Altläufe der Glatt, which were located approximately 5 km apart (47°27´ to 47°29´ N, 8°37´ to 8°32´ E, 417 to 572 m a.s.l., Canton of Zurich, Switzerland; Figure S1 and S2, Table S1). Mean annual temperature and precipitation are 9.8 ± 0.6 °C and 990 ± 168 mm (Kloten climate station 1988-2018; MeteoSchweiz, 2019). TFor this study, we used a space-for-time approach based on eight restoration sites that were between 3 and 32 years old. We measured recovery and restoration success by comparing the restored grasslands with intensively managed and semi-natural grasslands. Using a space-for-time approach requires high similarities in historical properties of the site, such as soil conditions and management regimes, to assure that temporal processes are appropriately represented by spatial patterns (Walker et al., 2010). This was the case in our study. The restored sites had similar soil conditions (i.e., soil type, structure, water availability) as the targeted semi-natural grasslands, while they shared the same agricultural legacy with intensively managed grasslands, i.e., biomass harvest and fertilization (manure and/or slurry) three to five times a year as well as tillage. We randomly established three 5 m x 5 m (25-m2) plots for plant identification and three 2 m x 2 m (4-m2) subplots for soil biotic and abiotic data collection at least 2 m away from the 25-m2 plots in each restoration site. Sites of similar age were grouped into four age classes: Y.4 (3 & 4 years after restoration), Y.18 (17 & 19 years), Y.24 (23 & 25 years), and Y.30 (27 & 32 years). Six intensively managed (Initial) and six semi-natural grassland (Target) sites complemented the experimental set-up, for a total of 36 plots. All plots were sampled under similar conditions, i.e., day of the year, air temperature, soil moisture, and time since last rain event, in June/July 2017 (intensively managed and semi-natural plots) and 2018 (restored plots).

Collection of plants and selected soil biota data
Plant species cover (in %) was visually estimated in each 25-m2 plot in mid-June (Braun-Blanquet, 1964; nomenclature: Lauber & Wagner, 1996). We calculated Shannon diversity and assessed plant community structure. We included soil microbial (fungi, procaryotes) and nematodes in our study as they represent the majority of soil biotic diversity and abundance (Bardgett & van der Putten, 2014), cover various trophic levels of the soil food web (Bongers & Ferris, 1999), and play key roles in soil functioning and ecosystem processes (Bardgett & van der Putten, 2014). In particular, soil nematodes were found to be well suited belowground indicators to evaluate recovery/development after restoration (e.g. Frouz, et al. 2008; Kardol et al., 2009; Resch et al., 2019). We randomly collected ten soil cores (2.2 cm diameter x 12 cm depths; sampler from Giddings Machine Company, Windsor, USA) in the 4-m2 subplots to assess soil nematode and microbial (fungal, prokaryotic) diversities and community structures. For soil nematodes, eight of the soil cores were combined and gently homogenized, placed in coolers and stored at 4 °C and transported to the laboratory (Netherlands Institute of Ecology, NIOO, Wageningen, Netherlands) within three days after collection. Free-living nematodes were extracted from 200 g of fresh soil using Oostenbrink elutriators (Oostenbrink, 1960). After extraction, each sample was divided into three subsamples, two for molecular identification and one to determine nematode abundance (see Resch et al., 2019). For the molecular work, two subsamples were stored in 70% ethanol (final volume 10 mL each) and transported to the laboratory at the Swiss Federal Research Institute WSL (Birmensdorf, Switzerland). Each subsample was reduced to roughly 200 μL by centrifugation and removal of the supernatant. The remaining ethanol was vaporized (65 °C for 3 h). Thereafter, 180 μL ATL buffer solution (Qiagen, Hilden, Germany) was immediately added and samples were stored at 4 °C until further processing. From these samples, nematode metagenomic DNA was extracted using the DNeasy Blood & Tissue Kit (Qiagen, Hilden, Germany) according to the manufacturer`s protocol, except for the incubation step which was run at 56 °C for 4 h. PCR amplification of the V6-V8 region of the eukaryotic small-subunit (18S) was performed with 7.5 μL of genomic DNA template (ca. 1 ng/μL) in 25 μL reactions containing 5 μL PCR reaction buffer, 2.5 mM MgCL2, 0.2 mM dNTPs, 0.8 μM of each primer (NemF: Sapkota & Nicolaisen, 2015; 18Sr2b: Porazinska et al., 2009), 0.5 μL BSA, and 0.25 μL GoTaq G2 Hot Start Polymerase (Promega Corporation, Madison, USA). Amplification was using an initial DNA denaturation step of 95 °C for 2 min, followed by 35 cycles at 94 °C for 40 sec, 58 °C for 40 sec, 72 °C for 1 min, and a final elongation step at 72 °C for 10 min. Filtering, dereplication, sample inference, chimera identification, and merging of paired-end reads was implemented using the DADA2 pipeline (v.1.12; Callahan et al., 2016) to finally assign amplicon sequence variants (ASVs) as taxonomic units.
We combined and homogenized the remaining two soil cores to assess soil microbes, placed them in coolers (4 °C) and transported them to the laboratory at WSL. Metagenomic DNA was extracted from 8 g sieved soil (2 mm) using the DNAeasy PowerMax Soil Kit (Qiagen, Hilden, Germany) according to the manufacturer´s protocol. PCR amplification of the V3-V4 region of the small-subunit (16S) of prokaryotes (i.e., bacteria and archaea) and the ribosomal internal transcribed spacer region (ITS2) of fungi was performed with 1 ng of template DNA using PCR primers and conditions as previously described (Frey et al., 2016). PCRs were run in triplicates, pooled and sent to the Genome Quebec Innovation Centre (Montreal, QC, Canada) for barcoding using the Fluidigm Access Array technology (Fluidigm) and paired-end sequencing on the Illumina MiSeq v3 platform (Illumina Inc., San Diego, USA). Quality filtering, clustering into operational taxonomic units (OTUs, 97% similarity cutoffs) and taxonomic assignment were performed as previously described (Resch et al., 2021).Taxonomic classification of nematode, prokaryotic and fungal sequences was conducted querying against the most recent versions of PR2 (v.4.11.1; Guillou et al., 2013), SILVA (v.132; Quast et al., 2013), and UNITE (v.8; Nilsson et al., 2019) reference sequence databases. Taxonomic assignment cutoffs were set to confidence rankings ≥ 0.8 (below ranked as unclassified). Prokaryotic OTUs assigned to mitochondria or chloroplasts as well as OTUs or ASVs assigned to other than Fungi or Nematoda were manually removed prior to data analysis. The three datasets were filtered to discard singletons and doubletons. Taxonomic abundance matrices were rarefied to the lowest number of sequences per community to achieve parity of the total number of reads between samples (Prokaryotes: 10,929 reads; Fungi: 18,337 reads; Nematodes: 6,662 reads). We calculated Shannon diversity and assessed community structures for soil nematodes, prokaryotes and fungi based on their relative abundances of ASV or OTU at the taxon level.

Collection of soil physical and chemical properties
We randomly collected one undisturbed soil core (5 cm diameter, 12 cm depth) per 4-m2 subplot using a steel cylinder that fit into the soil corer. The cylinders were capped to avoid disturbance during transport and used to measure field capacity, rock content and fine earth density as previously described (Resch et al., 2021). We randomly collected another three soil cores (5 cm diameter, 12 cm depths) in each 4-m2 subplot to determine soil chemical properties. The cores were pooled, dried at 60 °C for 48 h and passed through a 2 mm sieve. We measured soil pH (CaCl2) on dried samples, total nitrogen (N) and organic carbon (C) concentration on dried and fine-ground samples (≤ 0.5 mm; for details see Resch et al., 2021). We calculated total N and organic C pools after correcting its concentration for soil depth, rock content and fine earth density.

bacteria

display_name: bacteria

biodiversity

display_name: biodiversity

fungi

display_name: fungi

grassland

display_name: grassland

nematodes

display_name: nematodes

plants

display_name: plants

restoration

display_name: restoration

soil

display_name: soil

soil-biodiversity

display_name: soil-biodiversity

soil-biology

display_name: soil-biology

soil-chemistry

display_name: soil-chemistry

topsoil-removal

display_name: topsoil-removal

envidat

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: Long-term recovery of above-and belowground interactions in restored grasslands


maintainer: Anita C Risch

en: herbivore-exclusion climate microbial-respiration anthropogenic-management soil-properties nutrient-addition microbial-biomass-carbon

state: active

type: dataset

en: Risch et al. GEB_data.xlsx

en:
Study sites and experimental design
We collected data from 23 sites that are part of the Nutrient Network Global Research Cooperative (NutNet, https://nutnet.umn.edu/). The mean annual air temperature (MAT) across these sites ranged from -4 to 22°C, mean annual precipitation (MAP) from 252 to 1592 mm, and elevations from 6 to 4261 m above sea level (Fig 1a, Supplementary Table S1), hence cover a wide range of climatic conditions under which grasslands occur (Fig 1b). Soil organic C concentrations ranged between 0.8 to 7.8%, soil total N concentrations between 0.1 and 0.6%, and the soil C:N ratio between 9.1 and 21.5. Soil clay content spanned from 3.0 to 35%, and soil pH from 3.4 to 7.6 (Supplementary Table S2).
At each site, the effects of nutrient addition and herbivore exclusion were tested via a randomized-block design (Borer et al., 2014). Three blocks with 10 treatment plots each were established at each site, except for the site at bldr.us (only two blocks). Each of these 10 plots was randomly assigned to a nutrient or fencing treatment. An individual plot was 5 x 5 m, divided into four 2.5 x 2.5 m subplots. Each subplot was further divided into four 1 x 1 m square sampling plots, one of which was set aside for soil sampling (Borer et al., 2014). Plots were separated by at least 1 m wide walkways. We collected soil samples from four different treatments for this study: (i) untreated control plots (Control), (ii) herbivore exclusion plots (Fence), (iii) plots fertilized with N, P, K, plus nine essential macro and micronutrients (NPK), and (iv) plots with simultaneous fertilizer addition and herbivore exclusion (NPK+Fence). The experiments were established at different times in the past, with years of treatment different among sites (2 – 9 years since start of treatment; Supplementary Table S1). For the nutrient additions, all sites applied 10 g N m-2 each year as time-release urea; 10 g P m-2 yr-1 as triple-super phosphate; 10 g K m-2 yr-1 as potassium sulfate. A micro-nutrient mix (Fe, S, Mg, Mn, Cu, Zn, B, Mo, Ca) was applied at 100 g m-2 together with K in the first year of treatments but not thereafter.
We excluded large vertebrate herbivores (Fence) by fencing two plots, one with and one without NPK additions, within each block. The fences excluded all aboveground mammalian herbivores with a body mass of over 50 g (Borer et al., 2014). At most sites, the fences were 180 cm high, and the fence contained a wire mesh (1 cm holes) for the bottom 90 cm with a 30 cm outward-facing flange stapled to the ground to exclude burrowing animals. Climbing and subterranean animals may potentially still access these plots (Borer et al., 2014). For slight modifications in fence design at a few sites see Supplementary Table S3. Most sites only had wild herbivores, although four sites were also grazed by domestic animals (Supplementary Table S1).

Collection of soil samples, soil microbial respiration, microbial biomass, and other soil properties
Each of the 23 sites received a package containing identical material from the Swiss Federal Institute for Forest, Snow and Landscape Research WSL, Switzerland to be used for sampling (Risch et al., 2015; Risch et al., 2019). We collected two soil cores of 5 cm diameter and 12 cm depth in each sampling plot and composited them to measure MR, MBC, and soil chemical properties (see below). An additional sample (5 x 12 cm) was collected to assess soil physical properties. This sample remained within a steel sampling core after collection and both ends were tightly closed with plastic caps to avoid disturbance. All soils were shipped cooled to the laboratory at (Location will be disclosed after manuscript acceptance) within a few days after collection. Soils were sampled roughly 6 weeks prior to peak biomass at each site during 2015 and 2016.
To assess MR (CO2 production) in a laboratory incubation experiment we weighed duplicate soil samples (8 g dry soil equivalent) into 50-ml Falcon tubes. No additional substrate (glucose, sugar) was added to these samples. We adjusted the soil moisture of each sample to 60% field capacity. We then placed a 15 ml plastic test tube (Semadeni 1701A) containing 7.25 ml 0.05 M NaOH over each soil sample. The test tube was fixed with a plastic rod so that it was not in contact with the soil sample. The Falcon tubes were then sealed with a screw cap and placed in an incubator under completely dark conditions at 20°C. The CO2 produced by microbial respiration was absorbed by the 0.05 M NaOH. For five weeks we measured the decrease in conductivity within the 0.05 M NaOH solution on a weekly basis with a Multimeter WTW Multi 3410 (WTW GmbH, Germany) and replaced the 0.05 M NaOH with fresh solution. We included Falcon tubes without soil samples in each incubation run as blanks to test if tubes were tight and no CO2 could enter or escape. We calibrated the relationship between conductivity reduction and NaOH absorbed as follows: 400 ml 0.05 M NaOH was placed in a beaker and its conductivity was measured with the multimeter. While stirring, air containing CO2 was blown into the solution for approximately one minute, which reacts with NaOH to form Na2CO3. After this process, conductivity was measured again. We then transferred 7.25 ml of the solution into a smaller beaker and added 1 ml of 0.1 M BaCl2 to precipitate Na2CO3 and then titrated the solution with 0.05 M HCl to determine the remaining NaOH. We then repeated these steps with the remaining solution a total of nine times and plotted the conductivities (y-axis) against the NaOH consumed (x-axis, Supplementary Fig S1). This regression line was used to infer the consumption of NaOH from the conductivity reduction in the incubation experiments and to calculate CO2 evolution during incubation. In addition, we determined the optimum concentration for the NaOH solution in series of preliminary experiments, so that the concentration was not too high to become insensitive, but also not too low so that not all NaOH reacts during incubation. We then calculated MR (mg CO2-C kg dry soil-1 h-1) as total amount of CO2 released over the 5 weeks divided by the duration of the entire incubation in hrs.
Soil microbial biomass carbon (MBC; mg C kg soil-1 ) was measured at the beginning of the experiment by measuring the maximal respiratory response to the addition of glucose solution (4 mg glucose per g soil dry weight dissolved in distilled water; substrate-induced respiration method) on approximately 5.5 g of soil (J. P. E. Anderson & Domsch, 1978; Nico Eisenhauer et al., 2018; Scheu, 1992). For this purpose we used an O2-micro-compensation apparatus (Scheu, 1992). More specifically, substrate-induced respiration was calculated from the respiratory response to D-glucose for 10 hr at 20°C. Glucose was added according to preliminary studies to saturate the catabolic enzymes of microorganisms (4 mg g soil-1 dissolved in 400 ml deionized water). The mean of the lowest three readings within the first 10 hrs (between the initial peak caused by disturbing the soil and the peak caused by microbial growth) was taken as maximum initial respiratory response (MIRR; ml O2 kg soil-1 h-1) and microbial biomass (mg C kg soil-1) was calculated as 38 x MIRR (Beck et al., 1997; Cesarz et al., 2022; Thakur et al., 2015).
The rest of the composited sample was dried at 65°C for 48 h, ground and sieved (2 mm mesh) to assess the soil pH, mineral soil total C and N and C:N ratio, and mineral soil organic C (Risch et al., 2019). The undisturbed sample was used to assess water holding capacity (WHC), bulk density (BD), and soil texture [sand, silt, clay; methods in (Risch et al., 2019)]. We used the percentage of sand and clay as an indicator of soil texture in this study. MAT (°C), MAP (mm) and temperature of the wettest quarter (°C) were obtained from www.worldclim.com (Fick & Hijmans, 2017; Hijmans, Cameron, Parra, Jones, & Jarvis, 2005). These variables were selected as they were found to be drivers of soil nutrient processes across these sites in earlier studies (Risch et al., 2020; Risch et al., 2019). Mean annual soil temperatures (MAST; °C) for the 0 to 5 cm soil layer were obtained for each site from the SoilTemp maps (J. Lembrechts et al., 2021; J. J. Lembrechts et al., 2022), global gridded modelled products of soil bioclimatic variables for the 1979-2013 period at a 1-km² resolution, based on CHELSA, ERA5 and in-situ soil temperature measurements.
Numerical calculations and statistical analyses
We calculated MMQ as MR/MBC. We corrected this measure using the average soil temperature of each site (MMQsoil). This temperature correction is necessary as incubation temperatures are usually much higher than site mean annual soil temperatures (see Xu et al. 2017). MMQsoil = MMQ x Q10(MAST – 20)/10, where Q10 was assumed to be 2 (Xu et al. 2017). See Supplementary Fig S2 for comparison of air and soil temperatures across the 23 sites as well as the incubation temperature.
Some of the explanatory variables (clay, soil organic C, C:N ratio) were skewed and were thus log-transformed prior to analyses. All continuous explanatory variables were centred and scaled to have a mean of zero and variance of one. To avoid collinearity between them we filtered them using correlation analysis (Supplementary Fig S3). From the variables that were strongly correlated (Pearson’s |r| > 0.70) (Dormann et al., 2013), we selected the ones that allowed us to minimize the number of variables (Supplementary Fig S3). Specifically, soil total N concentration, soil total C concentration, soil sand content and soil bulk density were dropped from the dataset. We then assessed how these edaphoclimatic variables are related to MMQ across our global grasslands.
For this, we used linear mixed effects models (LMMs) fitted by maximum likelihood with the lme function in the nlme package (version 3.1-153) (Pinheiro, Bates, DebRoy, & Sarkar, 2021) in R version 3.6.3. (R Core Team, 2019). We used treatment as a fixed effect and plot nested in site as random effects to assess treatment differences in MMQsoil, as well as MR, and MBC. The number of years since the treatment started was included as a fixed effect in all the initial models but was not significant and therefore not retained in the models. To assess how differences in MMQsoil were affected by environmental factors (soil, climatic properties) we again used LMMs. Soil and climatic properties were included as fixed effects and plot nested in site as random effects. We did not include interactions between environmental variables. We then used the MuMin package (Barton, 2018) (version 1.42.1) to select the best models that explained the most variation based on Akaike’s information criterion (AIC; model.avg function). We used the corrected AIC (AICc) to account for our small sample size and selected the top models that fell within 4 AICc units (delta AICc < 4) (Burnham & Anderson, 2002; Johnson & Omland, 2004). We present all our top models rather than model averages. Conditional averages are provided in the Supplementary material.
Based on findings from analyses described above and the literature, we developed a conceptual model of direct and indirect relationships between both edaphoclimatic variables and experimental treatments (Supplement Figure S4) to obtain a more holistic approach in understanding how these properties affect MMQsoil. We had data from 23 sites with 272 observations. We tested this model using structural equation modelling based on a d-sep approach (Lefcheck, 2016; Shipley, 2009). We considered those environmental drivers that were included in our top LMMs, namely temperature of the wettest quarter (T.q.wet), soil pH, water holding capacity (WHC) and soil organic C (organic C; Supplementary Figure S4). These factors were allowed to directly affect MMQsoil, and via their interactions with treatments. In addition, treatments were allowed to directly affect MMQsoil. Treatments were included as dummy variables in the model. We tested our conceptual model (Supplementary Fig S4) using the piecewiseSEM package (version 2.0.2; Lefcheck, 2016) in R 3.4.0, in which a structured set of linear models are fitted individually. This approach allowed us to account for the nested experimental design, and overcome some of the limitations of standard structural equation models, such as small sample sizes (Lefcheck, 2016; Shipley, 2009). We used the lme function of the nlme package to model response variables, including site as a random factor. Good fit of the SEM was assumed when Fisher’s C values were non-significant (p > 0.05). For all significant interactions between soil or climate variables and treatments detected in the SEMs, we calculated treatment effect sizes, i.e., the differences in MMQsoil between Control and treatments as log response ratios (LRR) and plotted these values against the climate or soil factors. The LRR were defined as log(Control/Treatment), where treatment was either Fence, NPK or NPK+Fence. To assess which of the LRR-climate or soil property relationships were significant we again used LMMs, in which soil and climatic properties were included as fixed effects and plot nested in site as random effects.

format: SERVICE

en:
This dataset contains all data on which the following publication below is based.

Paper Citation:

Risch Anita C., Zimmermann, Stefan, Schütz, Martin, Borer, Elizabeth T., Broadbent, Arthur A.D., Caldeira, Maria C., Davies, Kendi F., Eisenhauer, Nico, Eskelinen, Anu, Fay, Philip A., Hagedorn, Frank, Knops, Johannes M.H., Lembrechts, Jonas, J., MacDougall, Andrew S., McCulley, Rebecca L., Melbourne, Brett A., Moore, Joslin L., Power, Sally A., Seabloom, Eric W., Silveira, Maria L., Virtanen, Risto, Yahdjian, Laura, Ochoa-Hueso, Raul (accepted). Drivers of the microbial metabolic quotient across global grasslands. Global Ecology and Biogeography

Please cite this paper together with the citation for the datafile.

The microbial metabolic quotient (MMQ; mg CO2-C mg MBC-1 h-1), defined as the amount of microbial CO2 respired (MR; mg CO2-C kg soil-1 h-1) per unit of microbial biomass C (MBC; mg C kg soil-1), is a key parameter for understanding the microbial regulation of the carbon (C) cycle, including soil C sequestration. Here, we experimentally tested hypotheses about the individual and interactive effects of multiple nutrient addition (NPK+micronutrients) and herbivore exclusion on MR, MBC, and MMQ across 23 sites (5 continents). Our sites encompassed a wide range of edaphoclimatic conditions, thus we assessed which edaphoclimatic variables affected MMQ the most and how they interacted with our treatments. Soils were collected in plots with established experimental treatments. MR was assessed in a five-week laboratory incubation without glucose addition, MBC via substrate-induced respiration. MMQ was calculated as MR/MBC and corrected for soil temperatures (MMQsoil). Using LMMs and SEMs, we analysed how edaphoclimatic characteristics and treatments interactively affected MMQsoil. MMQsoil was higher in locations with higher mean annual temperature, lower water holding capacity, and soil organic C concentration, but did not respond to our treatments across sites as neither MR nor MBC changed. We attributed this relative homeostasis to our treatments to the modulating influence of edaphoclimatic variables. For example, herbivore exclusion, regardless of fertilization, led to greater MMQsoil only at sites with lower soil organic C (<1.7%). Our results pinpoint the main variables related to MMQsoil across grasslands and emphasize the importance of the local edaphoclimatic conditions in controlling the response of the C cycle to anthropogenic stressors. By testing hypotheses about MMQsoil across global edaphoclimatic gradients, this work also helps to align the conflicting results of prior studies.

anthropogenic-management

display_name: anthropogenic-management

climate

display_name: climate

herbivore-exclusion

display_name: herbivore-exclusion

microbial-biomass-carbon

display_name: microbial-biomass-carbon

microbial-respiration

display_name: microbial-respiration

nutrient-addition

display_name: nutrient-addition

soil-properties

display_name: soil-properties

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

en: Education, culture and sport

it: Istruzione, cultura e sport

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

en: EnviDat: the Environmental Data Portal

it: EnviDat: il portale dei dati ambientali

fr: EnviDat est le portail de données du WSL; il offre un accès uniformisé et réglementé aux données de surveillance et de recherche en matière d’environnement. Ce portail est capable d’héberger et de publier des ensembles de données précieuses et uniques, qui sont partagées aux quatre coins de la planète. Ce faisant, EnviDat favorise la recherche et contribue à l’avancement des sciences environnementales.

de: EnviDat ist das Datenportal der WSL für einen einheitlichen und verwalteten Zugang zu Umweltdaten aus Monitoring- und Forschungsprojekten. Auf dem Portal können Datensätze gehostet und veröffentlicht werden. Durch die weltweite Bereitstellung einzigartiger und wertvoller Datensätze leistet EnviDat einen Beitrag zur Forschung im Bereich der Umweltwissenschaften.

it: EnviDat è il portale dati del WSL che fornisce un accesso unificato e gestito ai dati di monitoraggio ambientale e di ricerca. Il portale ha la capacità di ospitare e pubblicare serie di dati preziosi e unici che vengono condivisi in tutto il mondo. Così facendo, EnviDat promuove la ricerca e contribuisce al progresso della scienza ambientale.

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester


en: biogeochemical-cycling microbial-biomass plfa forest drought c13-labeling

en: Gao_Frank_C13 pulse labeling.xlsx

en:
Data of manuscript by Gao t al. : Drought resistance and resilience of rhizosphere communities in forest soils from the cellular to ecosystem scale – insights from 13C pulse labeling
Inludes data on 13C pulse-labeling experiment in young beech model ecosystems and a mature pine forest. Data of soil climate, total PLFA, 13C label strength in PLFA of pulse-labeling carried out under drought and after rewetting

format: XLS

en: Data of 13C-pulselabeling experiments in a young beech model ecosystem and a mature pine forest exposed to drought and rewetting. Data include microclimate during the experiment, total PLFA contents, d13C signatures in PLFA , microbial biomass, soil-respired CO2 measure during 30 days following 3-4 hour long pulse labeling with highly 13C-enriched CO2.

biogeochemical-cycling

display_name: biogeochemical-cycling

drought

display_name: drought

forest

display_name: forest

microbial-biomass

display_name: microbial-biomass

plfa

display_name: plfa

en: Gao_Drought resistance and resilience of rhizosphere communities_13C-PLFA


swiss-fluxnet-site-lageren

maintainer: EnviDat Support

title_for_slug: swiss-fluxnet-site-lageren

en: greenhouse-gas flux swiss-fluxnet swiss-forest-lab eddy-covariance meteorological-measurements forest

en: Data access

en: The flux and meteo data of the Swiss FluxNet site Lägeren are openly available from the European Fluxes Database Cluster.

en: Site description

en:
The Swiss FluxNet Site Lägeren is a managed mixed deciduous mountain forest located on the steep Lägeren mountain (NW of Zurich, Swiss Plateau). The forest is highly diverse, dominated by beech, but also including ash, maple, spruce and fir trees. Eddy covariance flux measurements were started in April 2004. The site was part of the international CarboEurope IP network and the National Air Pollution Monitoring Network (NABEL). In addition to Swiss FluxNet, the site is part of the Long-term Forest Ecosystem Research (LWF) of WSL and the biological drought and growth indicator network (TreeNet) of WSL.

Measurements

- Ecosystem flux measurements of CO2, H2O vapour are performed with the eddy-covariance method. This method is based on measurements of trace gas mixing ratios, using infrared gas analyzers (for CO2, H2O vapor), combined with wind speed and wind direction measurements, using 3D sonic anemometers. To resolve the short-term turbulent fluctuations in the atmosphere, very fast measurements are needed: we measure at 10-20 Hz, i.e., 10-20 times per second. To assess the energy budget of each ecosystem, also radiation sensors and soil climate profiles are installed at the site.
- Sub-canopy eddy fluxes (CO2, H2O), soil respiration campaigns
- Continuous CO2 profile measurements.
- Auxiliary micrometeorology and soil climate measurements.

Data availability

All data are available from the European Fluxes Database Cluster, but are also part of Fluxnet2015 dataset.

Data policy

ICOS data license:
[https://www.icos-cp.eu/data-services/about-data-portal/data-license](https://www.icos-cp.eu/data-services/about-data-portal/data-license)

Detailed site info:
[https://www.swissfluxnet.ethz.ch/index.php/sites/ch-lae-laegeren/site-info-ch-lae//](https://www.swissfluxnet.ethz.ch/index.php/sites/ch-lae-laegeren/site-info-ch-lae/)

eddy-covariance

display_name: eddy-covariance

flux

display_name: flux

greenhouse-gas

display_name: greenhouse-gas

meteorological-measurements

display_name: meteorological-measurements

swiss-fluxnet

display_name: swiss-fluxnet

swiss-forest-lab

display_name: swiss-forest-lab

en: Swiss FluxNet Site Lägeren


maintainer: Thorsten Bartels-Rausch

en: sea-salt-aerosol snow cryosphere antartic arctic polar laboratory-study

en: Experimental results of uptake coefficient vs. temperature and ozone concentration

media_type: text/csv

en:
Data to Figure 3: Measured uptake coefficients as a function of temperature at 5E12 molecules per cm3 ozone concentration.

Data to Figure 6: Measured uptake coefficients at different temperatures as a function of gas- phase ozone concentration.

format: CSV

mimetype: text/csv

en: Experimental Results: Uptake coefficients vs. temperature and ozone concentration

en:
* Final Data to Figure 3: Uptake coefficients as a function of temperature at 5E12 molecules per cm3 ozone concentration.
* Final Data to Figure 6: Uptake coefficients at different temperatures as a function of gas- phase ozone concentration.

The average ozone uptake coefficients of the three cycles between the bypass and the flow tube per repeat is taken as the mean uptake coefficient for each experiment. Gamma was corrected for gas phase diffusion limitation following the Cooney–Kim– Davis (CKD) method52 (see supporting material of publication for details).

en: The reaction of ozone with bromide in polar regions results in the formation of reactive bromide species with impacts on ozone budget and the oxidative capacity of the lower atmosphere. Here, we present a data investigating the temperature dependence of bromide oxidation by ozone using a coated wall flow tube reactor coated with an aqueous mixture of citric acid and sodium bromide, a proxy for sea salt aerosol in snow or the free troposphere. Thus study shows the effect of of organic species at relatively mild temperatures between the freezing point and eutectic temperature as typical for Earth's cryosphere.

antartic

display_name: antartic

arctic

display_name: arctic

cryosphere

display_name: cryosphere

laboratory-study

display_name: laboratory-study

polar

display_name: polar

sea-salt-aerosol

display_name: sea-salt-aerosol

snow

display_name: snow

envidat

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: Kinetic experiments on the oxidation of bromide by ozone from 289-245 K.


en: traffic-noise outdoor-restoration ppgis survey landscape-perception

en: Variable list for statistical analyses for article 1

media_type: application/pdf

en: Description of the sample stratification process for the invitation to participate in the survey

format: PDF

en: Processed survey data for window view for article 1

en: RESTORE WP4_Questionnaire_French.pdf

en: Questionnaire in French

en: Coding analysis first round for article 2

en: Coding analysis second round for article 2

en: Processed data for article 2

en: RESTORE WP4_Questionnaire_German.pdf

en: Questionnaire in German

en: Cleaned raw data from survey for article 1

en: Coding analysis third round for article 2

en: RESTORE WP4_Questionnaire_Italian.pdf

en: Questionnaire in Italian

en: Study protocol of the survey

en: Raw data for article 2

en: Variable list for article 2

en: Processed survey data of restorative locations for article 1

en: Processed data for coding analysis for article 2

en:
This repository contains data related to the Swiss-wide survey of the RESTORE project. This project is a collaboration between the Swiss Federal Institute for Forest, Snow and Landscape Research (WSL) and the Swiss Federal Laboratories for Materials Science & Technology (EMPA). It has received funding from the Swiss National Science Foundation. The objective of the RESTORE project was to assess the effects of green spaces as facilitators and noise as impediment to recover from stress in people's daily environments across Switzerland.
We conducted a nation-wide participatory mapping survey among Swiss residents to explore:
1.	The effects of greenness and noise exposure on psychological restoration when looking out of the window at home or when visiting nearby green spaces. It adopted a landscape approach to examine the interplay between environmental factors, individual perceptions, and personal traits – Scientific article 1
2.	The multiple aspects that constraint restoration when visiting nearby green spaces – Scientific article 2
Accordingly, this repository includes data and documentation for the following 2 scientific articles:
1.	Scientific article 1. “The role of greenness and road traffic noise for psychological restoration in everyday environments. A participatory mapping approach”. Landscape and Urban Planning, 2025
2.	Scientific article 2. “Perceived constraints for psychological restoration in nearby greenspaces. An exploratory and multi-dimensional approach”. Urban Forestry and Urban Greening, 2025

landscape-perception

display_name: landscape-perception

outdoor-restoration

display_name: outdoor-restoration

ppgis

display_name: ppgis

survey

display_name: survey

traffic-noise

display_name: traffic-noise

educ

en: Education, culture and sport

en: EnviDat: the Environmental Data Portal

en: Swiss-wide survey on restorative places and road traffic noise


en: herbivory climate-change drought-legacy phytochemistry fagus-sylvatica

en:
What are the research data files about:
Raw data on various beech (Fagus sylvatica) leaf traits. Beech leaf chemistry (primary and specialized metabolites), leaf damage measurements from various herbivore feeding guilds and crown health assessment

Which methods were used:
Field surveys on trees that showed either prematured defoliation due to drought-stress or showed not pre-mature defoliation. .

When and where was the data collected/analyzed:
Field data was collected in 2019, 2020
Data was analyzed in 2021-2023

climate-change

display_name: climate-change

drought-legacy

display_name: drought-legacy

fagus-sylvatica

display_name: fagus-sylvatica

herbivory

display_name: herbivory

phytochemistry

display_name: phytochemistry

envidat

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: Legacy effects of premature defoliation in response to an extreme drought event modulate phytochemical profiles with subtle consequences for leaf herbivory in European beech


en: radiation-transfer forest statistical-interpolation biodiversity ecology climate-change microclimate

media_type: application/zip

en: Predicted microclimate air temperatures at 100 cm above ground. 10 m spatial resolution raster files are provided for mean daily minimum (tmin), mean (tmean) and maximum (tmax) temperatures during each month from April (04) to October (10) over the period between 2012 to 2021.

format: ZIP

en: Predicted microclimate topsoil temperatures at 5 cm below ground. 10 m spatial resolution raster files are provided for mean daily minimum (tmin), mean (tmean) and maximum (tmax) temperatures during each month from April (04) to October (10) over the period between 2012 to 2021.

en: Predicted microclimate air temperatures at 5 cm above ground. 10 m spatial resolution raster files are provided for mean daily minimum (tmin), mean (tmean) and maximum (tmax) temperatures during each month from April (04) to October (10) over the period between 2012 to 2021.

en: Climate data matching the scales at which organisms experience climatic conditions are often missing. Yet, such data on microclimatic conditions are required to better understand climate-change impacts on biodiversity and ecosystem functioning. Here we combine a national network of microclimate temperature measurements with a novel radiative transfer model to map monthly minimum, mean and maximum temperatures during the vegetation period at a 10 meter spatial resolution across Switzerland. The temperature measurements took place in 107 sampling plots distributed across different habitat types, with 62 plots in forests, 22 below trees outside forests, and 23 in open grasslands. In each plot we measured temperature in the topsoil (-5cm), as well as in the air at 5cm and 100cm height above ground. Spatial interpolation was achieved by using a hybrid approach based on linear mixed effects models with input from detailed radiation estimates that account for topographic and vegetation shading, as well as other predictor variables related to the macroclimate, topography and vegetation height. Our data reveals strong horizontal and vertical variability in microclimate temperature, particularly for maximum temperatures at 5 cm above the ground and within the topsoil. Compared to macroclimate conditions as measured by weather stations outside forests, diurnal air and topsoil temperature ranges inside forests were reduced by up to 3.0 and 7.8 °C, respectively, while below trees outside forests, e.g. in hedges and below solitary trees, this buffering effect was 1.8 and 7.2 °C. We also found that in open grasslands, maximum temperatures at 5 cm above ground are on average 3.4 °C warmer than that of macroclimate, suggesting that in such habitats heat exposure close to the ground is often underestimated when using macroclimatic data. After accounting for macroclimate effects, microclimate patterns were primarily driven by radiation, with particularly strong effects on maximum temperatures. Results from spatial block cross-validation revealed predictive accuracies as measured by RSME’s ranging from 1.18 to 3.43 °C, with minimum temperatures generally being predicted more accurately than maximum temperatures. The microclimate maps presented here enable a more biologically relevant perspective when analysing climate-species interactions, which is expected to lead to a better understanding of biotic and ecosystem responses to climate and land use change.

biodiversity

display_name: biodiversity

ecology

display_name: ecology

forest

display_name: forest

microclimate

display_name: microclimate

radiation-transfer

display_name: radiation-transfer

statistical-interpolation

display_name: statistical-interpolation

en: Monthly topsoil and near surface microclimate temperature data for Switzerland


greenland-shrubs-and-microclimate

maintainer: Nathalie Chardon

title_for_slug: greenland-shrubs-and-microclimate

en: climate arctic validation accessible-shrub-forest field-data sdm

en: Kangerluarsunnguaq

en:
 Collected August 2020 (PI: Nathalie Isabelle Chardon)

These data collected August 2020 (PI: Nathalie Isabelle Chardon). Please contact N. I. Chardon (nathalie.chardon_at_gmail.com) for individual quadrat (i.e. 1 m x 1 m) data, which are averaged at the 90 m spatial grid here.

 COLUMN ABBREVIATIONS

 Metadata for 90 m spatial grid

id: sampling zone and unique 90 m grid identifier

long: WGS84 longitude

lat: WGS84 latitude

x_stere: longitude in Polar Stereographic (EPSG 4326)

y_stere: latitude in Polar Stereographic (EPSG 4326)

date: sampling day.month.year

 Species presence [1] or absence [0] at grid level

alcr: _Alnus crispa_

bena: _Betula nana_

dila: _Diapensia lapponica_

emni: _Empetrum nigrum_

hahy: _Harimanella hypoides_

lopr: _Loiseleuria procumbens_

oxpa: _Oxycoccus palustre_

phca: _Phyllodoce caerulea_

rhod: _Rhododendron_ sp.

saar: _Salix arctica_

sagl: _Salix glauca_

sahe: _Salix herbacea_

vaul: _Vaccinium uliginosum_

drin: _Dryas integrifolia_

thpr: _Thymus praecox_

 Habitat characteristics [% cover]

shrub: all shrub species

gram: graminoids

forb: forbs

bryo: bryophytes

lich: lichens

litt: litter

bare: bare ground

stone: stones (> 2 cm)

org: organic crust

 Shrub species height (maximum vegetative height) [mm] and cover (estimated in 1x1 m quadrat) [%]

alcr_c: _Alnus crispa_ cover

alcr_h: _Alnus crispa_ height

bena_c: _Betula nana_ cover

bena_h: _Betula nana_ height

dila_c: _Diapensia lapponica_ cover

dila_h: _Diapensia lapponica_ height

emni_c: _Empetrum nigrum_ cover

emni_h: _Empetrum nigrum_ height

hahy_c: _Harimanella hypoides_ cover

hahy_h: _Harimanella hypoides_ height

lopr_c: _Loiseleuria procumbens_ cover

lopr_h: _Loiseleuria procumbens_ height

oxpa_c: _Oxycoccus palustre_ cover

oxpa_h: _Oxycoccus palustre_ height

phca_c: _Phyllodoce caerulea_ cover

phca_height: _Phyllodoce caerulea_ height

rhod_c: _Rhododendron_ sp. cover

rhod_h: _Rhododendron_ sp. height

saar_c: _Salix arctica_ cover

saar_h: _Salix arctica_ height

sagl_c: _Salix glauca_ cover

sagl_h: _Salix glauca_ height

sahe_c: _Salix herbacea_ cover

sahe_h: _Salix herbacea_ height

vaul_c: _Vaccinium uliginosum_ cover

vaul_h: _Vaccinium uliginosum_ height

drin_c: _Dryas integrifolia_ cover

drin_h: _Dryas integrifolia_ height

 Community characteristics

rich_pa: shrub species richness at grid level

rich: shrub species richness averaged across 3 1x1 m sampling quadrats per grid

shan: Shannon diversity (‘vegan’ package; Oksanen et al. 2019)

comp_bena: shrub competition on _Betula nana_ (i.e. summed cover of non-_B. nana_ species) [% cover]

comp_sagl: shrub competition on _Salix glauca_ (i.e. summed cover of non-_S. glauca_ species) [% cover]

ht_cwm: maximum shrub canopy height (i.e. community weighted mean by abundance of maximum shrub height) [mm]

 Topography and terrain wetness (calculated from)

tcws: Tasseled Cap Wetness component (Sentinal images)

slope: slope angle (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

twi: SAGA wetness index (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015)

asp: aspect (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

 Downscaled climate (time series 1984 – 2013)

insoljja: average June-August insolation [MJ cm-2 yr-1]

precipmam: cumulative March-May precipitation [mm]

tempcont: yearly maximum – minimum temperature [ºC]

tempmin: yearly minimum temperature [ºC]

precipjja: cumulative June-August precipitation [mm]

tempjja: average maximum June-August temperature [ºC]

tempmax: yearly maximum temperature [ºC]

en: GitLab Repository

en: Zip drive of GitLab repository contents (https://gitlab.com/nathaliechardon/gl_microclim).

en: Nuup Kangerlua

en:
 Collected July-August 2011-2013 (PI: Jacob Nabe-Nielsen)

These data collected July-August 2011-2013 (PI: Jacob Nabe-Nielsen). Please contact J. Nabe-Nielsen (jnn_at_ecos.au.dk) to request access to data on all surveyed shrub species. Available here are only data used in Chardon et al. 2022.

 COLUMN ABBREVIATIONS

 Metadata for 90 m spatial grid

year: sampling year

alt: average altitude of grid as recorded during field surveys [m]

long: WGS84 longitude

lat: WGS84 latitude

x_stere: longitude in Polar Stereographic (EPSG 4326)

y_stere: latitude in Polar Stereographic (EPSG 4326)

 Species presence [1] or absence [0] at grid level

pa_bena: _Betula nana_

pa_sagl: _Salix glauca_

 Shrub species cover (calculated from pin hits in 0.7 x 0.7 m frame, based on total of 25 hits) [%]

cov_bena: _Betula nana_ cover

cov_sagl: _Salix glauca_ cover

 Downscaled climate (time series 1984 – 2013)

tempjja: average maximum June-August temperature [ºC]

tempmax: yearly maximum temperature [ºC]

tempmin: yearly minimum temperature [ºC]

tempcont: yearly maximum – minimum temperature [ºC]

precipmam: cumulative March-May precipitation [mm]

precipjja: cumulative June-August precipitation [mm]

insoljja: average June-August insolation [MJ cm-2 yr-1]

 Topography and terrain wetness (calculated from)

twi: SAGA wetness index (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015)

slope: slope angle (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

asp: aspect (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

tcws: Tasseled Cap Wetness component (Sentinal images)

 Community characteristics

rich: shrub species richness averaged across pin-point frames per grid

shan: Shannon diversity (‘vegan’ package; Oksanen et al. 2019)

comp_bena: shrub competition on Betula nana (i.e. summed cover of non-_B. nana_ species) [% cover]

comp_sagl: shrub competition on Salix glauca (i.e. summed cover of non-_S. glauca_ species) [% cover]

en:
 Study Aim

We collected these data to alternatively train and validate high resolution (~ 90 m) Species Distribution Models (SDMs) and Species Abundance Models (SAMs) for _Betula nana_ L. (dwarf birch, Betulaceae) and _Salix glauca_ L. (grey willow, Salicaceae) in Southwest Greenland to assess how well such models can predict local-scale patterns.

 Data Description

Individual (presence-absence, abundance, maximum vegetative height) and community (species composition, maximum canopy height) shrub data for two fjords near Nuuk, Southwest Greenland. Also provided are corresponding downscaled climate data as well as calculated topographic and terrain wetness indicator variables.

 Nuup Kangerlua (Godthåbsfjord)

_Betula nana_ and _Salix glauca_ presence-absence, abundance, community species richness

 Kangerluarsunnguaq (Kobbefjord)

Shrub presence-absence, abundance, maximum vegetative height, community composition, maximum shrub canopy height

 Methods

 Field survey in Nuup Kangerlua

We conducted a stratified systematic plant survey along the length of Nuup Kangerlua (NK) fjord in Soutwesth Greenland (Fig. 1 in Chardon et al. 2022; following Nabe-Nielsen et al., 2017). At five distinct sites, we sampled along elevational gradients to collect data on presences, absences, abundance, and species composition of all woody species using a 0.7 x 0.7 m pin-point frame (Fig. 1e in Chardon et al. 2022). For model training, we converted these pin-point data to percent cover estimates based on the number of pins dropped (n = 25 per plot) and averaged them across the 119 spatio-climatic grids (see next section) corresponding to the plot locations (for details see Appendix S2 in Chardon et al. 2022).

 Field survey in Kangerluarsunnguaq

We conducted a random stratified plant survey in Kangerluarsunnguaq (K) fjord in Southwest Greenland. We used a preliminary Species Abundance Model trained with summed pin counts of _Betula nana_ in NK fjord (see Fig. S1.3 in Chardon et al. 2022) to stratify the ~ 27 x 17 km fjord landscape into low, medium, and high abundances classes. We randomly selected 90 x 90 m spatio-climatic grids to survey in each class for a total of 200 grids, ensuring that they were accessible by foot or boat (for details see Appendix S2 in Chardon et al. 2022). Within each grid, we sampled within three 1 m2 quadrats arranged in a randomly rotated equilateral triangle centered on the mid-point of the cell. We used a gridded sampling quadrat with 1% delineations (Fig. 1h in Chardon et al. 2022) to record woody species presences, absences, and composition, estimated percent cover, and measured maximum shrub species vegetatitve height. At every plot, we also visually scanned the area in a 20 m radius from the plot and recorded the presence of any additional shrub species to estimate grid-level species richness. As in NK fjord, we averaged these data at the grid level (for details see Appendix S2 in Chardon et al. 2022).

 Biotic variables

We calculated biotic microscale variables from the plant survey data collected in NK and K fjords. We calculated shrub species richness, diversity, and competition (i.e. sum of non-B. nana or non-S. glauca pin hits or percent cover). In K fjord, we also calculated canopy height as the community weighted mean (by abundance) of maximum vegetative shrub height.

 Climate variables

We computed high resolution temperature, precipitation, and insolation for local scale data for the study area by statistically downscaling climate time series (1982 - 2013) from the monthly CHELSA data (Karger et al. 2017). We downscaled these data from 30 arc sec (~ 400 m at the latitude of our study) to our target grid size of ~ 90 m with geographic weighted regression and using the MEaSUREs Greenland Ice Mapping Project (GIMP) Digital Elevation Model (DEM) v. 1 (Howat et al., 2014, 2015). We then calculated 30-year averages of the climate parameters: average summer (June – August) maximum temperature, yearly maximum temperature, yearly minimum temperature, temperature continentality (yearly max. - min. temperatures), cumulative Spring (March – May) precipitation, cumulative summer precipitation, and average summer incident solar radiation (henceforth, insolation) (for calculation details see Appendices S2, S3 in Chardon et al. 2022 and Appendix S2 in von Oppen et al. 2021).

 Topography and terrain wetness indicator variables

We calculated several topographic and terrain wetness indices at a local scale. We derived slope, aspect, and the SAGA wetness index (hereafter TWI; Boehner et al., 2002; Boehner and Selige, 2006) from the GIMP DEM. TWI is a measure of how ‘wet’ an area is, based on water drainage from the surrounding landscape. We also calculated the tasseled cap wetness component (hereafter TCW, Crist and Cicone 1984) from satellite images (for details see Appendices S2, S3 in Chardon et al. 2022) as an alternative measure of wetness.

 Computer code

Attached as zip file and available on GitLab (https://gitlab.com/nathaliechardon/gl_microclim)

 Third-party data

Data used to calculate climate, topography, and terrain wetness indicator variables are publicly available (see Appendix S2 in Chardon et al. 2022 for all data references).

accessible-shrub-forest

display_name: accessible-shrub-forest

arctic

display_name: arctic

climate

display_name: climate

field-data

display_name: field-data

sdm

display_name: sdm

validation

display_name: validation

en: Greenland shrubs and microclimate


maintainer: Mahdi Jafari

en: numerical-simulation convection-in-snow openfoam vapor-transport snow-density-profile

en: snowpackBuoyantPimpleFoam

en:
This zipped folder includes (1) the source code files, (2) an example case for simulation, and (3) an instruction how to compile and use the new solver.

Note that the most updated version of snowpackBuoyantPimpleFoam can be found in the WSL/SLF GitLab repository: https://gitlabext.wsl.ch/snow-models/snowpackbuoyantpimplefoam

en: snowpackBuoyantPimpleFoam is a two-phase solver implemented to model convection of water vapor with phase change in snowpacks. This new solver is based on the standard solver of buoyantPimpleFoam in the open-source fluid dynamics software, OpenFOAM 5.0 (www.openfoam.org).

convection-in-snow

display_name: convection-in-snow

numerical-simulation

display_name: numerical-simulation

openfoam

display_name: openfoam

snow-density-profile

display_name: snow-density-profile

vapor-transport

display_name: vapor-transport

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: snowpackBuoyantPimpleFoam: an OpenFOAM Eulerian–Eulerian two-phase solver for modelling convection of water vapor in snowpacks


maintainer: Leo Gallus Bont

en: swiss-geodata digital-elevation-model qgis-plugin geoinformation-software cable-yarding seilaplan digital-planning swisstopo open-source

en: View video on youtube!

en: Herunterladen des DHM von swisstopo Website

en:
In order to use the QGIS plugin ‘Seilaplan’ for digital cable line planning, a digital terrain model (DTM) is required. As an alternative to using the ‘Swiss Geo Downloader’ plugin, the DTM can be obtained directly from Swisstopo. In this tutorial we explain step by step how to download the necessary DTM from the Swisstopo Website, and how to use it in QGIS for the digital planning of a cable line using the plugin ‘Seilaplan’. Please note that the tutorial language is German!

Link to the elevation model on the swisstopo website: https://www.swisstopo.admin.ch/de/geodata/height/alti3d.html#technische_details

Link to the rope map website: https://seilaplan.wsl.ch

********************

Für die Verwendung des QGIS Plugins Seilaplan zur digitalen Seillinienplanung ist ein digitales Höhenmodell (DHM) nötig. Als Alternative zum Swiss Geo Downloader erklären wir in diesem Tutorial Schritt für Schritt, wie man das nötige Höhenmodell von der Swisstopo Webseite herunterladen und in QGIS zur Seillinienplanung verwenden kann.

Link zum Höhenmodell auf der swisstopo Webseite: https://www.swisstopo.admin.ch/de/geodata/height/alti3d.html#technische_details

Link zur Seilaplan-Website: https://seilaplan.wsl.ch

cable-yarding

display_name: cable-yarding

digital-elevation-model

display_name: digital-elevation-model

digital-planning

display_name: digital-planning

geoinformation-software

display_name: geoinformation-software

open-source

display_name: open-source

qgis-plugin

display_name: qgis-plugin

seilaplan

display_name: seilaplan

swiss-geodata

display_name: swiss-geodata

swisstopo

display_name: swisstopo

en: Seilaplan Tutorial: DTM download from swisstopo website


carabid-beetles-in-forests

title_for_slug: carabid-beetles-in-forests

en: carabidae forest arthropoda running-beetles trapping coleoptera insects

media_type: text/plain

en: Data extract in TAB-delimted format. First row is header, first column is counter.

format: TXT

en: Tab-delimited data of the according sql query.

en: SQL query used to extract the data from fauna database.

en: Data query used on Fauna Database to retrieve the submitted data.

en:
Carabidae data from all historic up to the recent projects (21.10.2019) of WSL, collected with various methods in forests of different types.

Version 2 ('FIDO_global_extract 2019-11-22_18-11-24 WSL-Forest-Carabidae') contains additional data field PROJ_FALLENBEZEICHNUNG.

Data are provided on request to contact person against bilateral agreement.

arthropoda

display_name: arthropoda

carabidae

display_name: carabidae

coleoptera

display_name: coleoptera

forest

display_name: forest

insects

display_name: insects

running-beetles

display_name: running-beetles

trapping

display_name: trapping

en: Carabid beetles in forests


en: precipitation multi-year-drought chelsa extreme-events megadrought drought potential-evapotranspiration

en: Standardized Precipitation Index (SPI)

en: Global monthly layers of the Standardized Precipitation Index (SPI) a 30 arcsec (~1km) spatial resolution

en: Standardised Precipitation Evapotranspiration Index (SPEI)

en: Global monthly layers of the Standardized Precipitation Evapotranspiration Index (SPEI) a 30 arcsec (~1km) spatial resolution

en: Statistics of the multiyear meteorological droughts, containing the ID, duration, area, and empirical ranking of the droughts.

en: Multiyear meteorological droughts (MYMD) statistics - Description

en: Description of column names, and units of the fields contained in the dataset "Multiyear meteorological droughts (MYMD) statistics".

en: Technical documentation document including the variable names, units, scaling, offset, and background.

en: Global annual layers of the multiyear meteorological droughts (MYMD) at a 1/24° (~5km) spatial resolution

en:
Multiyear meteorological droughts (MYMD10) statistics
Statistics of the multiyear meteorological droughts, containing the ID, duration, area, and empirical ranking of the droughts, and the anomalies in kNDVI

en: Global annual layers of the multiyear meteorological droughts (MYMD) at a 1/12° (~10km) spatial resolution

en: Multiyear meteorological droughts (MYMD10) statistics - Description

en: Description of column names, and units of the fields contained in the dataset "Multiyear meteorological droughts (MYMD10) statistics".

en: Kernel Normalized Difference Vegetation Index Anomaly (QKNDVI)

en: Global annual layers of the the kernel normalized difference vegetation index anomaly (qkndvi) at a 1/12° (~10km) spatial resolution.

en:
The dataset consists of the standardized precipitation (SPI) and the standardized precipitation evapotranspiration index (SPEI) index with a 30 arcsec (~1km) horizontal resolution.
Data for precipitation rates (pr) and potential evapotranspiration (pet) are taken from CHELSA V2.1 (https://doi.org/10.16904/envidat.228.v2.1). Both data sets are monthly time series from 1980 to 2018 at 30 arcsec. and have the unit of kg m-2 month-1. SPEI is a more reliable measure of drought than SPI as it additionally considers the effect of potential evapotranspiration. Potential evapotranspiration (pet) is the amount of water per area per time that can evaporate at the soil surface or be transpired by plants in absence of water supply limitations. SPEI effectively estimates for a given location how the climatic water balance, i.e., the difference between precipitation rate (pr) and potential evapotranspiration (pet), relates to the long-term mean. Here, pet was calculated with the Penman-Montheith approach, assuming surface conductance of a reference crop of 12 cm height following the definition of the Food and Agriculture Organization of the United Nations (FAO), originating from CHELSA-BIOCLIM+ database (https://www.doi.org/10.16904/envidat.332). To obtain global time series of SPEI, in each pixel the frequency distribution of water balance estimates (pr-pet) is approximated with a log-logistic probability distribution function. We calculate monthly SPEI using a 12-month memory and approximate the parameters of the log-logistic probability distribution based on the time series from 1980-2018. Such a 12-month time window is suitable to capture long-term drought events that are long enough to impose substantial impacts on agriculture, hydrology, and ecosystems while ones at shorter time scales would be more appropriate to detect short-term drought. For each month, SPEI is calculated considering the respective water-balance value and the values of the eleven preceding months. Here we use SPEI-12 in December to represent the overall water-balance condition of the respective year. The resulting SPEI data has a horizontal resolution of 130 arcsec (~1km). The calculation of SPI and SPEI is done in R with the package SPEI (https://www.doi.org/10.32614/CRAN.package.SPEI).

chelsa

display_name: chelsa

drought

display_name: drought

extreme-events

display_name: extreme-events

megadrought

display_name: megadrought

multi-year-drought

display_name: multi-year-drought

potential-evapotranspiration

display_name: potential-evapotranspiration

precipitation

display_name: precipitation

en: High resolution global standardized drought indices


en: soil-moisture plants soil-mapping ecological-maps landolt ecological-indicator-values soil-ph switzerland

en: Description

en: Description of the predictors provided in this dataset

en: Predictors

en: Predictors generated and used in Descombes et al. (2020)

mimetype: application/zip

en:
Ecologically meaningful predictors are often neglected in plant distribution studies, resulting in incomplete niche quantification and low predictive power of species distribution models (SDMs). Because environmental data are rare and expensive to collect, and because their relationship with local climatic and topographic conditions are complex, mapping them over large geographic extents and at high spatial resolution remains a major challenge. Here, we derived environmental data layers by mapping ecological indicator values (EIVs) in space by using a large set of environmental predictors in Switzerland.

This dataset contains the predictors (raster layers) generated and used in the following publication (Descombes et al. 2020). Only predictors for which we have the rights to share them are provided. Other datasets and predictors can be accessed via the original data provider. Details on the predictors and sources are fully described in the publication. The predictors are provided as GeoTIFF files, at 93 m spatial resolution and Mercator projection ("+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +ellps=WGS84 +towgs84=0,0,0,0,0,0,0 +units=m +no_defs"). The excel file (xlsx) provides a short description of the raster layers.

Paper Citation:

Descombes, P. et al. (2020). Spatial modelling of ecological indicator values improves predictions of plant distributions in complex landscapes. Ecography. (accepted)

ecological-indicator-values

display_name: ecological-indicator-values

ecological-maps

display_name: ecological-maps

landolt

display_name: landolt

plants

display_name: plants

soil-mapping

display_name: soil-mapping

soil-moisture

display_name: soil-moisture

soil-ph

display_name: soil-ph

switzerland

display_name: switzerland

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.


en: arctic chemistry snow alpine solution brine antarctica interface heterogeneous ozone laboratory-study reaction metamorphism

en: Readme file

en: This file shows data, code, and plots.

mimetype: application/pdf

en: Matlab Data

en: All data as Matlab mat file

en: Ozone data

en: Trend of ozone mixing ratio ( _O3mixingRatio_ppb_ with time ( _time_s_ ) as measured downstream of the snow sample during the experiments. Time "0" denotes the time when data acquisition with ozone analyser was started. And, ozone (_nOzone_molec_per_s_ ) loss for each experiment with duration of ozone exposure ( _duration_s_ ) of ozone exposure, as derived from the ozone trends shown above. Time "0" denotes the time when the carrier gas with ozone was passed via the snow sample.

en: Matlab Script to view and plot data

en: Data Figure 1

en: Data of published Figure

en: data list

en: Lists file names, snow sample properties, and flow settings for each data-set.

en: Acetone data

en:
intensity of the acetone-water-cluster and water-cluster (as derived by a chemical ionisation mass spectrometer) with time in the carrier gas downstream of packed bed flow tube holding the snow sample or the by-pass, respectively. Time "0" denotes the time when the mass spectrometer data acquisition was started. When the carrier
gas with acetone was passed via the packed bed flow tube, that is switched from by-pass to pass the snow, acetone concentration in the carrier gas rapidly dropped due to surface adsorption and slowly recovered as the adsorption sites saturate. From this, the surface area can be obtained.

en: Earth’s snow cover is very dynamic on diurnal time scales. The changes to the snow structure during this metamorphism have wide ranging impacts such as on avalanche formation and on the capacity of surface snow to exchange trace gases with the atmosphere. Here, we investigate the influence of dry metamorphism, which involves fluxes of water vapor, on the chemical reactivity of bromide in the snow. For this, the heterogeneous reactive loss of ozone at a concentration of 5-6E12 molecules cm-3 is investigated in artificial, shock-frozen snow samples doped with 6.2 uM sodium bromide and with varying metamorphism history. The oxidation of bromide in snow is one reaction initiating polar bromine releases and ozone depletions.

alpine

display_name: alpine

antarctica

display_name: antarctica

arctic

display_name: arctic

brine

display_name: brine

chemistry

display_name: chemistry

heterogeneous

display_name: heterogeneous

interface

display_name: interface

laboratory-study

display_name: laboratory-study

metamorphism

display_name: metamorphism

ozone

display_name: ozone

reaction

display_name: reaction

snow

display_name: snow

solution

display_name: solution

en: Data set on bromide oxidation by ozone in snow during metamorphism from laboratory study


en: numerical-simulation wrf complex-terrain grayzone cosmo-wrf

en:
This dataset contains simulation data for *"Influence of Air Flow Features on Alpine Wind Energy Potential" (Kristianti et al., 2024).*
•	The WRF_dia2702 and WRF_dia1103 folders contain the namelist.wps, namelist.input, and output files from the simulation in the Diablerets region on 27/02/2021 and 11/03/2021, respectively.
•	The WRF_luk2310 and WRF_luk0912 folders contain the namelist.wps, namelist.input, and output files from the simulation in the Lukmanier region on 23/10/2020 and 09/12/2020, respectively.
•	The postprocessing file contains Python codes used to create figures in Kristianti et al. (2024).
**Notes:**
•	Topography input is provided by Gerber et al. (2021).
•	Simulations are run using the COSMO-WRF model developed by Gerber et al. (2018).
•	Visualization by Python code uses wrf-python provided by Ladwig et al. (2017).
*This dataset was produced by the EDGE consortium sponsored by the Swiss Federal Office of Energy's SWEET programme, Swiss National Science Foundation (SNSF): Grant 179130, and the Swiss National Supercomputing Centre (CSCS) projects s938, s1115 and s1242.*

complex-terrain

display_name: complex-terrain

cosmo-wrf

display_name: cosmo-wrf

grayzone

display_name: grayzone

numerical-simulation

display_name: numerical-simulation

wrf

display_name: wrf

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: COSMO-WRF Dataset for Swiss Alps Simulations in Gray-Zone Resolution


maintainer: Robert Pazúr

en: switzerland cropland habitat phenology grassland agriculture vegetation-indices remote-sensing sentinel-2 machine-learning

en: Cropland and grassland map of Switzerland  (v.1.0)

media_type: image/tiff

en: !!! Please check out the newer version in this repository !!! (1) cropland; (2) grassland; other LC classes are masked out from the map

format: TIFF

en: (current version) Cropland, grassland, and shrubland map of Switzerland (v.1.5)

format: GeoTIFF

en: (current version) Cropland and grassland map of Switzerland (v.1.5)

en: Cropland, grassland, and shrubland map of Switzerland (v1.0)

en: !!! Please check out the newer version in this repository !!! (1) cropland; (2) grassland; (3) shrubland; other LC classes are masked out from the map

en: Mask of unsealed surfaces of Switzerland based on 3 years Sentinel-2 NDVI 95th percentile threshold value: 0.6  (v1.0)

en:
!!! Please check out the newer version in the repository !!!
Mask of unsealed surfaces of Switzerland based on 3 years (2017-2019) Sentinel-2 NDVI 95th percentile threshold value: 0.6

en: (current version) Mask of unsealed surfaces of Switzerland based on 3 years Sentinel-2 NDVI 95th percentile threshold value: 0.7 (v1.5)

en: GEE (Google Earth Engine) script (JavaScript) to derived the output map  (v1.5)

en: The code used in Google Earth Engine to derive the cropland, grassland and shrubland map of Switzerland. Compare to the previous version the newer version includes the shadow masking procedure which might be important in the mountainous areas.

en: (current version) Mask of unsealed surfaces of Switzerland based on 3 years Sentinel-2 NDVI 95th percentile threshold value: 0.6 (v1.5)

en: Assuming that sealed surfaces have little or no phenological response we calculated a mask of unsealed surfaces based on the 3 years (2017-2019) Sentinel-2 NDVI 95th percentile threshold value: 0.6.

en: Mask of unsealed surfaces of Switzerland based on 3 years Sentinel-2 NDVI 95th percentile threshold value: 0.7 (v1.0)

en: !!! Please check out the newer version in this repository !!! Assuming that sealed surfaces have little or no phenological response we calculated a mask of unsealed surfaces based on the 3 years (2017-2019) Sentinel-2 NDVI 95th percentile threshold value: 0.7.

en: GEE (Google Earth Engine) script (JavaScript) to derived the output map  (v1.0)

en: !!! Please check out the newer version in the repository !!! The code used in Google Earth Engine to derive the cropland, grassland and shrubland map of Switzerland

en: We developed a map of cropland and grassland allocation for Switzerland based on several indices dominantly derived from Sentinel-2 satellite imagery captured over multiple growing seasons. The classification model was trained based on parcel-based data derived from landholder reporting. The mapping was conducted on Google Earth Engine platform using random forest classifier. Areas of high vegetation, shrubland, sealed surface and non-vegetated areas were masked out from the country-wide map. The resulting map has high accuracy in lowlands as well as mountainous areas.

agriculture

display_name: agriculture

cropland

display_name: cropland

grassland

display_name: grassland

habitat

display_name: habitat

machine-learning

display_name: machine-learning

phenology

display_name: phenology

remote-sensing

display_name: remote-sensing

switzerland

display_name: switzerland

vegetation-indices

display_name: vegetation-indices

en: Cropland and grassland map of Switzerland based on Sentinel-2 data


maintainer: Elisabeth Hafner

en: natural-hazard satellite avalanche hazard-mapping snow

en: SPOT6 Avalanche Outlines 16 January 2019

en: Shapefile of the manually mapped avalanche outlines. Please only use this data together with the example key document provided with this dataset and within the limitations of the dataset described in the publications.

format: SHP

en:
Outlines of 6'041 avalanches mapped from SPOT6 satellite data over the Swiss Alps on 16 January 2019. The dataset was acquired following a period with very high avalanche danger. The outlines have different attributes described in the data example key (ExampleKey_AvalMapping16012019.pdf)

The generation of the data is described in:

Bühler, Y., E. D. Hafner, B. Zweifel, M. Zesiger, and H. Heisig (2019), Where are the avalanches? Rapid SPOT6 satellite data acquisition to map an extreme avalanche period over the Swiss Alps, The Cryosphere, 13(12), 3225-3238, doi:10.5194/tc-13-3225-2019.

The data was comprehensivly validated in a subset area in

Hafner, E.D.; Techel, F.; Leinss, S.; Bühler, Y., 2021: Mapping avalanches with satellites - evaluation of performance and completeness. Cryosphere, 15, 2: 983-1004. doi: 10.5194/tc-15-983-2021

avalanche

display_name: avalanche

hazard-mapping

display_name: hazard-mapping

natural-hazard

display_name: natural-hazard

satellite

display_name: satellite

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: SPOT6 Avalanche outlines 16 January 2019


en: wind-direction anemometer wind-velocity wind-speed wind raclets

en: Sonic temperature

en: The data set contains values of temperature (°C) measured with the anemometer at a sampling frequency of 20 Hz.

en: Sonic time

en: The data set contains the time in YYYY-mm-dd HH:MM:SS.sss for the anemometer data. The time is given in UTC+1.

en: Sonic wind component u

en:
The data contains values for u direction of the wind in a sampling frequency of 20 Hz.
Erroneour spikes were filterd in two steps.

1. Define moving window uₓ = u(i-3:i+3)
2. If u(i)-mean(uₓ) > 2*std(uₓ):u(i) = NaN

Frozen detectors resulted in continuous zero values (not filtered).
The wind speed and direction were calculated as follows:

wind speed = sqrt(u² + v²)
wind direction = 180° + (180° / PI)* arctan2 (-u, -v)

en: Sonic wind component w

en:
The data contains values for w direction of the wind in a sampling frequency of 20 Hz.
Erroneour spikes were filterd in two steps.

1. Define moving window uₓ = u(i-3:i+3)
2. If u(i)-mean(uₓ) > 2*std(uₓ):u(i) = NaN

Frozen detectors resulted in continuous zero values (not filtered).
The wind speed and direction were calculated as follows:

wind speed = sqrt(u² + v²)
wind direction = 180° + (180° / PI)* arctan2 (-u, -v)

en: Sonic wind component v

en:
The data contains values for v direction of the wind in a sampling frequency of 20 Hz.
Erroneour spikes were filterd in two steps.

1. Define moving window uₓ = u(i-3:i+3)
2. If u(i)-mean(uₓ) > 2*std(uₓ):u(i) = NaN

Frozen detectors resulted in continuous zero values (not filtered).
The wind speed and direction were calculated as follows:

wind speed = sqrt(u² + v²)
wind direction = 180° + (180° / PI)* arctan2 (-u, -v)

en: A Young 81000 sonic anemomenter was deployed at Gotschnagrat (LON: 46.859 LAT: 9.849) to record three components of the wind velocity (u, v, w in [m s‾ ¹]) and air temperature (Ts in [°C]). The anemomenter was mounted in direction North at a height of 1.5 m above snow surface at the beginning. The time within each data set is given in UTC+1. Instrument specifications can be found [here](http://www.youngusa.com/Manuals/81000-90(I).pdf) .

anemometer

display_name: anemometer

raclets

display_name: raclets

wind

display_name: wind

wind-direction

display_name: wind-direction

wind-speed

display_name: wind-speed

wind-velocity

display_name: wind-velocity

en: Snow Drift Station - 3D Ultrasonic


en: snow-water-equivalent natural-climate-variability uncertainty climate-change-impact cryosphere rain-on-snow

en: climate realizations

en: The dataset contains for three variables (snow water equivalent, surface water input and liquid precipitation) 50 realizations of current and future climate periods for two time horizons (mid end end of century), two emission senarions (RCP 4.5 and 8.5) and 10 climate model chains (all EUR11 chains within CH2018). To quantify natural climate variability for projections of snow conditions and resulting rain-on-snow (ROS) flood events, a weather generator was applied to simulate inherently consistent climate variables for multiple realizations of current and future climates at 100 m spatial and hourly temporal resolution over a 12 x 12 km high-altitude study area in the Swiss Alps. The output of the weather generator was used as input for subsequent simulations with an energy balance snow model. The data was extracted in 2021 from original model output.

climate-change-impact

display_name: climate-change-impact

cryosphere

display_name: cryosphere

natural-climate-variability

display_name: natural-climate-variability

rain-on-snow

display_name: rain-on-snow

snow-water-equivalent

display_name: snow-water-equivalent

uncertainty

display_name: uncertainty

educ

en: Education, culture and sport

en: Multiple realizations of daily snow water equivalent, surface water input and liquid precipitation projections for mid- and late-century


en: snow uas snow-depth swe remote-sensing

en: .tfw file

en: coordinates for the .tif raster

en: HS map Brämabül, 14. April 2015

en:
10 cm resolution
15 cm accuracy

mimetype: image/tiff

en:
This snow depth map was generated 14 January 2015, close to peak of winter accumulation, applying Unmanned Aerial System digital surface models with a spatial resolution of 10 cm. The covered area is 285'000 m2 at the top of Brämabüel, 2490 m a.s.l. covering all expositions. Coordinate system: CH1903LV03.

A detailed description is given here:

Bühler, Y., Adams, M. S., Bösch, R., and Stoffel, A.: Mapping snow depth in alpine terrain with unmanned aerial systems (UASs): potential and limitations, The Cryosphere, 10, 1075-1088, 10.5194/tc-10-1075-2016, 2016.

Abstract:
Detailed information on the spatial and temporal distribution, and variability of snow depth (HS) is a crucial input for numerous applications in hydrology, climatology, ecology and avalanche research. Nowadays, snow depth distribution is usually estimated by combining point measurements from weather stations or observers in the field with spatial interpolation algorithms. However, even a dense measurement network is not able to capture the large spatial variability of snow depth in alpine terrain.
Remote sensing methods, such as laser scanning or digital photogrammetry, have recently been successfully applied to map snow depth variability at local and regional scales. However, such data acquisition is costly, if manned airplanes are involved. The effectiveness of ground-based measurements on the other hand, is often hindered by occlusions, due to the complex terrain or acute viewing angles. In this paper, we investigate the application of unmanned aerial systems (UAS), in combination with structure-from-motion photogrammetry, to map snow depth distribution. Such systems have the advantage that they are comparatively cost-effective and can be applied very flexibly to cover also otherwise inaccessible terrain. In this study we map snow depth at two different locations: a) a sheltered location at the bottom of the Flüela valley (1900 m a.s.l.) and b) an exposed location (2500 m a.s.l.) on a peak in the ski resort Jakobshorn, both in the vicinity of Davos, Switzerland. At the first test site, we monitor the ablation on three different dates. We validate the photogrammetric snow depth maps using simultaneously acquired manual snow depth measurements. The resulting snow depth values have a root mean square error (RMSE) better than 0.07 to 0.15 m on meadows and rocks and a RMSE better than 0.30 m on sections covered by bushes or tall grass. This new measurement technology opens the door for efficient, flexible, repeatable and cost effective snow depth monitoring for various applications, investigating the worlds cryosphere.

remote-sensing

display_name: remote-sensing

snow

display_name: snow

snow-depth

display_name: snow-depth

swe

display_name: swe

uas

display_name: uas

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat is the WSL data portal providing unified and managed access to environmental monitoring and research data. The portal has the capability to host and publish data sets. Through the sharing of unique and high-value data sets worldwide, EnviDat fosters research and contributes to advancing environmental science.

en: UAS based snow depth maps Brämabüel, Davos, CH


en: time-series sentinel-2 video-encoding geodata research-data-management core raster gridded-data

en: View publication in the Digital Object Repository at WSL (DORA)

media_type: text/html

format: HTML

en:
Shell script to extract data from the CORE format. You can 1) download and unzip the file, 2) make the script executable and 3) run it.
Disclaimer: only provided for peer review. Proper software will become available in time on https://github.com/EnviDat/?q=core (both Web/JS for Frontend and Python for Backend).

en:
Shell script to reproduce the CORE format. You can 1) download and unzip the file, 2) make the script executable and 3) run it.
Disclaimer: only provided for peer review. Proper software for working with CORE will become available in time on https://github.com/EnviDat/?q=core (both Web/JS for Frontend and Python for Backend).

en: CORE Dataset (80 GB)

en:
DISCLAIMER: CORE is still in development. Interested parties are warmly invited to join common development, to comment, discuss, find bugs, etc.

Acknowledgement: The CORE format was proudly inspired by the Cloud Optimized GeoTIFF ([COG](https://www.cogeo.org/)) format, by considering how to leverage the ability of clients issuing HTTP GET range requests for a time-series of remote sensing and aerial imagery (instead of just one image).

License: The Cloud Optimized Raster Encoding (CORE) specifications are released to the public domain under a Creative Commons 1.0 CC0 "No Rights Reserved" international license. You can reuse the information contained herein in any way you want, for any purposes and without restrictions.

-----------------------

Summary:
The Cloud Optimized Raster Encoding (CORE) format is being developed for the efficient storage and management of gridded data by applying video encoding algorithms. It is mainly designed for the exchange and preservation of large time series data in environmental data repositories, while in the same time enabling more efficient workflows on the cloud. It can be applied to any large number of similar (in pixel size and image dimensions) raster data layers. CORE is not designed to replace COG but to work together with COG for a collection of many layers (e.g. by offering a fast preview of layers when switching between layers of a time series).

WARNING: Currently only applicable to RGB/Byte imagery. The final CORE specifications may probably be very different from what is written herein or CORE may not ever become productive due to a myriad of reasons (see also 'Major issues to be solved').

With this early public sharing of the format we explicitly support the Open Science agenda, which implies "shifting from the standard practices of publishing research results in scientific publications towards sharing and using all available knowledge at an earlier stage in the research process" (quote from: European Commission, Directorate General for Research and Innovation, 2016. Open innovation, open science, open to the world).

CORE Specifications:

1) a MP4 or WebM video digital multimedia container format (or any future video container playable as HTML video in major browsers)

2) a free to use or open video compression codec such as H.264, VP9, or AV1 (or any future video codec that is open sourced or free to use for end users)

Note: H.264 is currently recommended because of the wide usage with support in all major browsers, fast encoding due to acceleration in hardware (which is currently not the case for AV1 or VP9) and the fact that MPEG LA has allowed the free use for streaming video that is free to the end users. However, please note that H.264 is restricted by patents and its use in proprietary or commercial software requires the payment of royalties to [MPEG LA](https://www.mpegla.com/programs/avc-h-264/). However, when AV1 matures and accelerated hardware encoding becomes available, AV1 is expected to offer 30% to 50% smaller file size in comparison with H.264, while retaining the [same quality](https://trac.ffmpeg.org/wiki/Encode/AV1).

3) the encoding frame rate should be of one frame per second (fps) with each layer segmented in internal tiles, similar to COG, ordered by the main use case when accessing the data: either layer contiguous or tile contiguous;

Note: The internal tile arrangement should support an easy navigation inside the CORE video format, depending on the use case.

4) a CORE file is optimised for streaming with the moov atom at the beginning of the file (e.g. with -movflags faststart) and optional additional optimisations depending on the codec used (e.g. -tune fastdecode -tune zerolatency for H.264)

5) metadata tags inside the moov atom for describing and using geographic image data (that are preferably compatible with the [OGC GeoTIFF standard](https://www.ogc.org/standards/geotiff) or any future standard accepted by the geospatial community) as well as list of original file names corresponding to each CORE layer

6) it needs to encode similar source rasters (such as time series of rasters with the same extent and resolution, or different tiles of the same product; each input raster should be having the same image and pixel size)

7) it provides a mechanism for addressing and requesting overviews (lower resolution data) for a fast display in web browser depending on the map scale (currently external overviews)

Major issues to be solved:

- Internal overviews (similar to COG), by chaining lower resolution videos in the same MP4 container for fast access to overviews first); Currently, overviews are kept as separate files, as external overviews.

- Metadata encoding (how to best encode spatial extent, layer names, and so on, for each of the layer inside the series, which may have a different geographical extent, etc...; Known issues: adding too many tags with FFmpeg which are not part of the standard MP4 moov atom; metadata tags have a limited string length.

- Applicability beyond RGB/Byte datasets; defining a standard way of converting cell values from Int16/UInt16/UInt32/Int32/Float32/Float64/ data types into multi-band Byte values (and reconstructing them back to the original data type within acceptable thresholds)

Example

Notice: The provided CORE (.mp4) examples contain modified Copernicus Sentinel data [2018-2021].

For generating the CORE examples provided, 50 original Sentinel 2 (S-2) TCI data images from an area located inside Switzerland were downloaded from www.copernicus.eu, and then transformed into CORE format using ffmpeg with H.264 encoding using the [x264 library](https://www.videolan.org/developers/x264.html).

For full reproducibility, we provide the original data set and results, as well scripts for data encoding and extraction (see resources).

core

display_name: core

geodata

display_name: geodata

gridded-data

display_name: gridded-data

raster

display_name: raster

research-data-management

display_name: research-data-management

time-series

display_name: time-series

video-encoding

display_name: video-encoding


en: plant-traits climate-change soil-nutrients plant-community treeline-ecosystems biomass-dynamics

en: Data on soils, tree biomass, plant nutrients, rock mineralogy, climate of treeline ecosystems in Khibiny mountains, Kola Penninsula, Russia

en: Data on soil nutrients and charcteristics, tree and understory biomass, foliage nutrients, mineralogy, soil temperature and moisture along an elevation gradient from boreal forest to mountain tundra in teh Khibiny mountains (Kola Penninsula, Russia). Samples were taken along two gradients at 7 elevation level either under tree canopy (tree) and in open areas (open).

biomass-dynamics

display_name: biomass-dynamics

climate-change

display_name: climate-change

plant-community

display_name: plant-community

plant-traits

display_name: plant-traits

soil-nutrients

display_name: soil-nutrients

treeline-ecosystems

display_name: treeline-ecosystems

educ

en: Education, culture and sport