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

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: 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: 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: 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: 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: 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: 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:
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

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


wind-lidar-davos-wolfgang

maintainer: EnviDat Support

title_for_slug: wind-lidar-davos-wolfgang

en: profiles wolfgangpass wind raclets lidar

state: active

type: dataset

en: Wind LIDAR data Wolfgang

media_type: application/zip

en: A Leosphere WidCube 100S was deployed at Davos Wolfang, measuring wind profiles as well as performing PPI and RHI scans. Time format is in UTC with highest time resolution of 5 seconds.

format: ZIP

mimetype: application/zip

en: Scanning wind Lidar from Meteoswiss was installed at Davos Wolfgang (LON: 9.853594, LAT: 46.835577) and measured from 200 m above ground to 8100 m. The time resolution is up to 5 seconds. The Lidar was measuring wind profiles but also performed plan position indicator (PPI) and range height indicator (RHI) scans.

lidar

display_name: lidar

profiles

display_name: profiles

raclets

display_name: raclets

wind

display_name: wind

wolfgangpass

display_name: wolfgangpass

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

en: Education, culture and sport

it: Istruzione, cultura e sport

envidat

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.

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.

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: Wind LIDAR Davos Wolfgang


en: icos dbh atmospheric-conditions dendrometer seehornwald davos carbon-storage open-data subalpine carbon-cycle carbon-dioxide eddy-covariance-measurements climate-change climate-growth sap-flow icp-forests tree-stem treenet flux eddy-covariance forest-characteristics automated-monitoring co2 ecosystem-functioning forest forest-function forest-biomass forest-monitoring picea-abies stem-growth forest-ecosystems metadata carbon-sequestration lwf ecophysiology forest-canopy forest-dynamics ilter biomass water-relations twd forest-area

en: Collection of files with information about the research site Seehornwald Davos CH-DAV

en: Readme

en: Readme file contains the list of material that comes with the download

format: SERVICE

en: Site description of the Seehornwald Davos (CH-DAV) research site. The material provided offers information on the biomass and growth of trees in CH-DAV, including supporting materials such as maps, geolocations and analyses. CH-DAV is part of networks such as ICOS, LWF, TreeNet, ICPForests and is dedicated to forest monitoring and research.

atmospheric-conditions

display_name: atmospheric-conditions

automated-monitoring

display_name: automated-monitoring

biomass

display_name: biomass

carbon-cycle

display_name: carbon-cycle

carbon-dioxide

display_name: carbon-dioxide

carbon-sequestration

display_name: carbon-sequestration

carbon-storage

display_name: carbon-storage

climate-change

display_name: climate-change

climate-growth

display_name: climate-growth

davos

display_name: davos

dbh

display_name: dbh

dendrometer

display_name: dendrometer

ecophysiology

display_name: ecophysiology

ecosystem-functioning

display_name: ecosystem-functioning

eddy-covariance

display_name: eddy-covariance

eddy-covariance-measurements

display_name: eddy-covariance-measurements

flux

display_name: flux

forest

display_name: forest

forest-area

display_name: forest-area

forest-biomass

display_name: forest-biomass

forest-canopy

display_name: forest-canopy

forest-characteristics

display_name: forest-characteristics

forest-dynamics

display_name: forest-dynamics

forest-ecosystems

display_name: forest-ecosystems

forest-function

display_name: forest-function

forest-monitoring

display_name: forest-monitoring

icos

display_name: icos

icp-forests

display_name: icp-forests

ilter

display_name: ilter

lwf

display_name: lwf

metadata

display_name: metadata

open-data

display_name: open-data

picea-abies

display_name: picea-abies

sap-flow

display_name: sap-flow

seehornwald

display_name: seehornwald

stem-growth

display_name: stem-growth

subalpine

display_name: subalpine

tree-stem

display_name: tree-stem

treenet

display_name: treenet

twd

display_name: twd

water-relations

display_name: water-relations

en: Site description Davos Seehornwald CH-DAV


maintainer: Elisabeth Hafner

en: natural-hazard avalanche satellite 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

snow

display_name: snow

en: SPOT6 Avalanche outlines 16 January 2019


en: tree-species switzerland time-series sentinel-2 remote-sensing

media_type: text/csv

en: Tree species profiles for 2019, SWIR-2 band

format: CSV

en: Tree species profiles for 2019, SWIR-1 band

en: Tree species profiles for 2019, Red-Edge-1 band

en: Tree species profiles for 2020, CCI index

en: Tree species profiles for 2019, NDMI/NDWI index

en: Tree species profiles for 2019, NDVI index

en: Tree species profiles for 2020, Green band

en: Tree species profiles for 2019, Red-Edge-2 band

en: Tree species profiles for 2020, Red-Edge-1 band

en: Tree species profiles for 2019, Blue band

en: Tree species profiles for 2019, CCI index

en: Tree species profiles for 2020, SWIR-1 band

en: Tree species profiles for 2019, NIR band

en: Tree species profiles for 2020, NDMI index

en: Tree species profiles for 2019, Green band

en: Tree species profiles for 2019, EVI index

en: Tree species profiles for 2020, SWIR-2 band

en: Tree species profiles for 2019, Red band

en: Tree species profiles for 2020, CIre index

en: Tree species profiles for 2019, CIre Index

en: Tree species profiles for 2019, Red-Edge-3 band

en: Tree species profiles for 2020, NIR band

en: Tree species profiles for 2020, Red-Edge-3 band

en: Tree species profiles for 2020, Red-Edge-2 band

en: Tree species profiles for 2020, EVI index

en: Tree species profiles for 2020, NDVI index

en: Tree species profiles for 2020, Blue band

en: Tree species profiles for 2020, Red band

en:
 Tree species profiles from Sentinel-2 time series
We obtained tree species-specific Sentinel-2 time series for the years 2019 and 2020. The Sentinel-2 time series originated from [here](https://www.envidat.ch/#/metadata/sentinel-2-time-series-of-switzerland), and the tree species information was derived from pure plots of the Swiss National Forest Inventory, spread all over Switzerland.
In addition to the 10 and 20 m bands from Sentinel-2 (Blue, Green, Red, Red-Edge-1, Red-Edge-2, Red-Edge-3, NIR, SWIR-1, SWIR-2), we calculated indices (CCI, CIre, EVI, NDMI/NDWI, NDVI).
In total, tree species profiles of 763 pixels for 2020 and 562 pixels for 2019 are available for seven tree species (*Abies alba*, *Castanea sativa*, *Fagus sylvatica*, *Fraxinus excelsior*, *Larix* spp., *Picea abies*, *Pinus sylvestris*).

remote-sensing

display_name: remote-sensing

switzerland

display_name: switzerland

time-series

display_name: time-series

tree-species

display_name: tree-species


en: raster-processing sentinel-2 time-series switzerland remote-sensing force

media_type: image/tiff

format: TIFF

en:
We processed Sentinel-2 image time series from 2017 to 2023 for Switzerland with the Software FORCE (Frantz 2019) on the basis of [Sentinel-2 images](https://envidat.ch/#/metadata/sentinel-2-imagery-of-switzerland). The respective parameter files can be found here: [Github](https://github.com/TLKoch/Sentinel-2_CH). All the image time series consist of several TB and therefore access will be granted upon request.
The available bands (in spatial reference system EPSG 3035) are the following:
Red, Green, Blue, NIR, Red-Edge-1, Red-Edge-2, Red-Edge-3, SWIR-1, SWIR-2
The available indices (in spatial reference system EPSG 3035) are the following:
CCI, CIRE, NDWI/NDMI, NDVI, EVI
Further indices might be provided upon request.

 Processing
On the basis of processed Sentinel-2 images for the 14 Sentinel-2 tiles covering Switzerland (T31TGN, T32TLT, T32TMT, T32UMU, T32TNT, T32TPT, T31TGM, T32TLS, T32TMS, T32TNS, T32TPS, T32TLR, T32TMR, T32TNR), we processed the image time series further with FORCE v. 3.7.8-12.
We generated interpolated Sentinel-2 time series with a 5-day interval, corresponding to the theoretical revisit time of the Sentinel-2 satellites. It's important to note that the 5-day time series consist of interpolated and smoothed composites, not the original images. We used the radial basis convolutional filtering (RBF) available in the FORCE time series analysis (TSA) submodule (Schwieder et al. 2016). The RBF is similar to a spatial moving window average approach over time (Schwieder et al. 2016). We applied kernel width values of 10, 20, 30, and 50 days. We spectrally adjusted all the images to match Sentinel-2A, and we removed curve outliers and pixels that failed the quality checks for clouds and their shadows, snow, saturation, and limited illumination. The processed image time series are available in tiles of 30 by 30 km.

 Example images
Uploaded is an example of the index EVI for one of the generated 30 by 30 km tiles located around the city of Zürich. The values are multiplied by 10.000. The time series spans the month of July from 2018.

force

display_name: force

raster-processing

display_name: raster-processing

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: Sentinel-2 time series of Switzerland


greenland-shrubs-and-microclimate

maintainer: Nathalie Chardon

title_for_slug: greenland-shrubs-and-microclimate

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

en: GitLab Repository

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

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: 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


en: growing-degree-days growing-season-length precipitation relative-humidity temperature soil-water-balance solar-radiation

en: CHELSA-BIOCLIM+

en: Global extented set of bioclimatic variables at different temporal (e.g. monthly, annual, long term normals) resolution in Geotiff format.

format: GeoTIFF

en: CHELSA-BIOCLIM+ Technical Documentation

media_type: application/pdf

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

format: PDF

en: A multitude of physical and biological processes on which ecosystems and human societies depend are governed by climatic conditions. Understanding how these processes are altered by climate change is central to mitigation efforts. Based on mechanistically downscaled climate data, we developed a set of climate-related variables at yet unprecedented spatiotemporal detail as a basis for environmental and ecological analyses. We created gridded data for near-surface relative humidity (hurs), cloud area fraction (clt), near-surface wind speed (sfcWind), vapour pressure deficit (vpd), surface downwelling shortwave radiation (rsds), potential evapotranspiration (pet), climate moisture index (cmi), and site water balance (swb), at a monthly temporal and 30 arcsec spatial resolution globally starting 1980 until 2018. At the same spatial resolution, we further estimated climatological normals of frost change frequency (fcf), snow cover days (scd), potential net primary productivity (npp), growing degree days (gdd), and growing season characteristics for the periods 1981-2010, 2011-2040, 2041-2070, and 2071-2100, considering three shared socioeconomic pathways (SSP126, SSP370, SSP585) and five Earth system models. Time-series variables showed high accuracy when validated against observations from meteorological stations. Climatological normals were also highly correlated to observations although some variables showed notable biases, e.g., snow cover days (scd). Together, the data sets presented here allow improving our understanding of patterns and processes that are governed by climate, including the impact of recent and future climate changes on the world’s ecosystems and associated services to societies.

growing-degree-days

display_name: growing-degree-days

growing-season-length

display_name: growing-season-length

precipitation

display_name: precipitation

relative-humidity

display_name: relative-humidity

soil-water-balance

display_name: soil-water-balance

solar-radiation

display_name: solar-radiation

temperature

display_name: temperature

en: CHELSA-BIOCLIM+ A novel set of global climate-related predictors at kilometre-resolution


weather-station-davos-wolfgang

maintainer: Jörg Wieder

title_for_slug: weather-station-davos-wolfgang

en: raclets temperature pressure weather relative-humidity

en: Weather Station Davos Wolfgang

en: The data file contains date and time (YYYY-MM-DD hh:mm:ss) in UTC, temperature (°C), relative humidity (%), pressure (hPa) at 10 m above ground and above sea level. The data was averaged over 10 minutes.

mimetype: text/csv

en: The dataset contains weather parameters measured at Davos Wolfgang (LON: 9.853594, LAT: 46.835577).

pressure

display_name: pressure

raclets

display_name: raclets

weather

display_name: weather

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.


maintainer: Dirk Nikolaus Karger

en: temperature climatology climate-change precipitation bioclim climate free-climate-data climate-layers

en: CHELSA V1.2 (deprecated)

en: CHELSA Version 1.2 based on ERA-Interim forcings and CMIP5

en: CHELSA V2.1 (current)

en: CHELSA Version 2.1 based on ERA5 forcings and CMIP6

en: CHELSA Technical Documentation

en:
High-resolution information on climatic conditions is essential to many applications in environmental and ecological sciences. Here we present the CHELSA (Climatologies at high resolution for the earth’s land surface areas) data of downscaled temperature and precipitation to a high resolution of 30 arc sec. The temperature algorithm is based on statistical downscaling of atmospheric temperatures. The precipitation algorithm incorporates orographic predictors including wind fields, valley exposition, and boundary layer height, with a subsequent bias correction.

CHELSA data published in EnviDat includes the deprecated version 1.2 (originally published under 10.5061/dryad.kd1d4). Please use the current 2.1 version.

Paper Citation:
 > _Karger DN. et al. Climatologies at high resolution for the earth’s land surface areas, Scientific Data, 4, 170122 (2017) [doi: 10.1038/sdata.2017.122](https://doi.org/10.1038/sdata.2017.122)._

bioclim

display_name: bioclim

climate-change

display_name: climate-change

climate-layers

display_name: climate-layers

climatology

display_name: climatology

free-climate-data

display_name: free-climate-data


maintainer: Frank Techel

en: avalanche-accident-statistics avalanche-accident avalanche-fatalities long-term-trends avalanche-statistics

en: Data avalanche fatalities 1969/1970 - 2014/2015

media_type: text/plain

en: The data comprises avalanche fatality data from the Alpine regions in the seven countries Austria, France, Germany, Italy, Liechtenstein, Slovenia, and Switzerland.

format: TXT

mimetype: text/plain

en:
During the last 45 years, about 100 people lost their lives in avalanches in the European Alps each year. Avalanche fatalities in settlements and on transportation corridors have considerably decreased since the 1970s. In contrast, the number of avalanche fatalities during recreational activities away from avalanche-secured terrain doubled between the 1960s and 1980s and has remained relatively stable since, despite a continuing strong increase in winter backcountry recreational activities.

Data complementing Figure 2 in:
_"Avalanche fatalities in the European Alps: long-term trends and statistics"_, by Techel, F., Jarry, F., Kronthaler, G., Mitterer, S., Nairz, P., Pavšek, M., Valt, M., and Darms, G.

Data description:
please refer to section 2 (Data and Methods) in the mentioned publication

avalanche-accident

display_name: avalanche-accident

avalanche-accident-statistics

display_name: avalanche-accident-statistics

avalanche-fatalities

display_name: avalanche-fatalities

avalanche-statistics

display_name: avalanche-statistics

long-term-trends

display_name: long-term-trends

en: Avalanche fatalities in the European Alps (1969/1970 - 2014/2015)


maintainer: Michael Lehning

en: snow-cover snow-distribution snow-hydrology snow alpine-precipitation alpine3d-simulations alpine3d

en: ADS datasets and simulation inputs

en: The file contains the ADS datasets of snow depth and input data used for simulations.

en:
Dataset (Model input, snow distribution and validation) for the precipitation scaling paper, which should be cited along with the data set citation.

This data is useful for distributed hydrological modelling or other tasks that involve the study of snow distribution and precipitation in the high Alpine.

The format of the data is for Alpine3D (models.slf.ch) model runs but other models could be used, too.

Please cite:

_Vögeli, C., Lehning, M., Wever, N., Bavay M., 2016: Scaling Precipitation Input to Spatially Distributed Hydrological Models by Measured Snow Distribution., Front. Earth Sci. 4: 108. doi: 10.3389/feart.2016.00108._

Dataset is provided as a single zip file.
The archive contains two directories, the valuable distributed snow depth maps for the landscape Davos and the simulation input.
The archive also contains the file: "ReadMeMetadataDataSetPrecipitationScaling" which explains the data structure.

alpine-precipitation

display_name: alpine-precipitation

snow

display_name: snow

snow-cover

display_name: snow-cover

snow-distribution

display_name: snow-distribution

snow-hydrology

display_name: snow-hydrology

en: Precipitation Scaling Data Set (Vögeli et al., Frontiers)


en: aridity precipitation frost-change-velocity ground-temperature sdms soil

en: Soil organic carbon content modified from SoilGrids

en: Soil pH modified from SoilGrids

en: Frost change frequency modified from CHELSA

en: Mean annual temperature modified from CHELSA

en: Precipitation seasonality from CHELSA

en: Precipitation of driest quarter modified from CHELSA

en: Soil clay content modified from SoilGrids

en: Mean diurnal range modified from CHELSA

en: The dataset contains seven environmental layers (average annual temperature, aridity [annual precipitation divided by annual potential evapotranspiration], frost change frequency, precipitation in the driest quarter, mean diurnal temperature range, and precipitation seasonality) modified from CHELSA (https://chelsa-climate.org/) and three soil layers (soil organic matter content, pH water, and clay content) modified from SoilGrids (https://soilgrids.org/).

aridity

display_name: aridity

frost-change-velocity

display_name: frost-change-velocity

ground-temperature

display_name: ground-temperature

precipitation

display_name: precipitation

sdms

display_name: sdms

soil

display_name: soil

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: Environmental layers for SDM simulations (GDPlants)


maintainer: M. Carol Resch

en: topsoil-removal vegetation long-term-experiment grassland nematodes restoration diversity

en: Nematode feeding types

en:
Plot: Plot ID combining treatment code and plot number (in total 55 plots)

Treatment: Three restoration treatments and two reference treatments. Initial=intensively managed grasslands; Harvest only= multiple hay harvest (2-3 times) per years; Topsoil= Topsoil removal (upper 10-30 cm); Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Target= species-rich ancient grasslands

He.ecto.a: Herbivorous ectoparasite nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

He.endo.a: Herbivorous endoparasite nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Ba.a: Bacterivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Fu.a: Fungivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Om.a: Omnivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Ca.a: Carnivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

format: XLS

en: Coordinates of one research plot were wrong. This was corrected in the current version of the file.

en: Nematode colonizer-persister classes

en:
Plot: Plot ID combining treatment code and plot number (in total 55 plots)

Treatment: Three restoration treatments and two reference treatments. Initial=intensively managed grasslands; Harvest only= multiple hay harvest (2-3 times) per years; Topsoil= Topsoil removal (upper 10-30 cm); Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Target= species-rich ancient grasslands

CP1.a: C-P (colonizer-persister) class 1 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990; equivalent to enrichment indicators according to Ferris, Bongers, & de Goede 2001
CP2.a: C-P (colonizer-persister) class 2 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990; equivalent to stress indicators according to Ferris, Bongers, & de Goede 2001
CP3.a: C-P (colonizer-persister) class 3 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990
CP4.a: C-P (colonizer-persister) class 4 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990
CP5.a: C-P (colonizer-persister) class 5 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990

en: Plant community composition (plant species)

en:
Plot: Plot ID combining treatment code and plot number (in total 55 plots)

Treatment: Three restoration treatments and two reference treatments. Initial=intensively managed grasslands; Harvest only= multiple hay harvest (2-3 times) per years; Topsoil= Topsoil removal (upper 10-30 cm); Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Target= species-rich ancient grasslands

column 3-252: Plant abundance per species (in total 252 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: Nematode families

en:
Plot: Plot ID combining treatment code and plot number (in total 55 plots)

Treatment: Three restoration treatments and two reference treatments. Initial=intensively managed grasslands; Harvest only= multiple hay harvest (2-3 times) per years; Topsoil= Topsoil removal (upper 10-30 cm); Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Target= species-rich ancient grasslands

column 3-41: Nematode abundance (individuals per 100g dry soil) per family; nomenclature according to Bongers 1988

en: Nemadotes and vegetation

en:
Plot: Plot ID combining treatment code and plot number (in total 55 plots)

Treatment: Three restoration treatments and two reference treatments. Initial=intensively managed grasslands; Harvest only= multiple hay harvest (2-3 times) per years; Topsoil= Topsoil removal (upper 10-30 cm); Topsoil+Propagules= Topsoil removal (upper 10-30 cm) combined with introduction of target plant species; Target= species-rich ancient grasslands

X: X-Coordinate

Y: Y-Coordinate

Ri.veg: Total number of plant species

Div.veg: Shannon diversity of plant species

Ri.nem: Total number of nematode families

Div.nem: Shannon diversity of nematode families

abund: Total nematode abundance (individuals per 100g dry soil)

MI: Maturity index according to Bongers 1990; calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

PPI: Plant parasite index according to Bongers 1990; calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

PPI.MI: Plant parasite to Maturity index according to Bongers, van der Meulen, & Korthals 1997; calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

CI: Channel index according to Ferris, Bongers, & de Goede 2001 (%); calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

EI: Enrichment index according to Ferris, Bongers, & de Goede 2001 (%); calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

SI: Structure index according to Ferris, Bongers, & de Goede 2001 (%); calculations done in online tool NINJA (Sieriebriennikov, Ferris, & de Goede, 2014)

He.a: Herbivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

He.ecto.a: Herbivorous ectoparasite nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

He.endo.a: Herbivorous endoparasite nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Ba.a: Bacterivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Fu.a: Fungivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Om.a: Omnivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Ca.a: Carnivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014

Om.Ca.a	Sum of omnivorous and carnivorous nematode abundance (individuals per 100g dry soil); classification according to Sieriebriennikov, Ferris, & de Goede, 2014
CP1.a: C-P (colonizer-persister) class 1 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990; equivalent to enrichment indicators according to Ferris, Bongers, & de Goede 2001

CP2.a: C-P (colonizer-persister) class 2 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990; equivalent to stress indicators according to Ferris, Bongers, & de Goede 2001

CP3.a: C-P (colonizer-persister) class 3 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990

CP4.a: C-P (colonizer-persister) class 4 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990

CP5.a: C-P (colonizer-persister) class 5 nematode abundance (individuals per 100g dry soil); C-P classification according to Bongers 1990

CP3-5.a: Sum of C-P3, C-P4 and C-P5 nematode abundance (individuals per 100g dry soil) equivalent to structure indicators according to Ferris, Bongers, & de Goede 2001

protect.sp.pro: plant species categorized into species of concern based on their red list status in Switzerland as well as their protection status in Switzerland and the Canton Zurich (Moser, Gygax, Bäumler, Wyler, & Palese, 2002)

tar.sp: Number of target species (all species recorded in the eleven target plots plus propagules of species applied by hand, resulting in a total of 143 species)

Indi.moist: Plant indicator values for soil moisture according to Landolt et al. 2010

Indi.nutri: Plant indicator values for soil nutrients according to Landolt et al. 2010

Management: level of mowing and manuring frequency (35-40 kg N ha-1 application-1); intensity1 = mowing once and no manure; intensity2 = 2-3 times mowing and no manure; intensity3 = 2-3 times mowing and manuring; intensity4 = 4-5 times mowing and manuring

Tillage: last time of tillage

Land.use.history: for Initial plots = recent agricultural use; for all restored sites = agricultural use before restoration
Mowing	earliest date of biomass harvesting (excluding “Initial”) based on contracts between farmers and the Nature Protection Agency of Zurich

Soil.type: based on soil mapping survey 1991 (before restoration)

Soil.texture: determined by the sediment method (Gee & Bauder 1986) and classified into soil textural classes according to FAO (1990)

pH: soil pH (CaCl2)

Slope: slope gradient classes according to FAO (1990; %)

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

Paper Citation:
 > Resch, M.C., Schütz, M., Graf, U., Wagenaar, R., van der Putten, W.H., Risch, A.C. 2019. Does topsoil removal in grassland restoration benefit both soil nematode and plant communities? Journal of Applied Ecology 56: 1782-1793.

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

 Methods
 Study area and experimental settings
The study was conducted in a nature reserve (Eigental: 47° 27’ to 47° 29’ N, 8° 37’ E, 461 to 507 m a.s.l.) that is located on the Swiss Central plateau close to Zurich airport (Canton Zurich, Switzerland). The mean annual temperature in this area ranges from 8.9 to 10.6 °C, mean annual precipitation from 910 to 1260 mm [10-year average (2007-2017); MeteoSchweiz, 2018]. The main soil types are calcaric to gleyic Cambisol and Gleysols. The reserve was established in 1967 to protect small remnants of oligotrophic semi-natural grasslands (roughly 12 ha). The plant community can be characterized as Molinion and Mesobromion (semi-wet to semi-dry), depending on the site-specific groundwater level and slope inclination (Delarze, Gonseth, Eggenberg, & Vust, 2015). These remnants represent species-rich islands in an otherwise intensively managed agricultural landscape. Semi-natural grasslands covered an area of 60,000 ha in the Canton Zurich in 1939, however, by 2005 only roughly 600 ha remained (Baudirektion Kanton Zürich, 2007). In 1990, the government of Canton Zurich decided to enlarge the nature reserve Eigental. The goal was to incorporate eleven patches of 20 ha adjacent intensively farmed land and transform these patches into semi-natural grasslands. The patches had a different agricultural history, ranging from permanent (no tillage for >50 years) to temporary grassland (as part of crop rotation; last tillage <5 years). On all freshly integrated patches fertilization was stopped in 1992 and from then on biomass was harvested three times a year and removed. After 5 years without noticeable effects on vegetation composition, the Nature Conservation Agency of Canton Zurich decided to increase the restoration efforts. In 1995, a large-scale experiment was initialized to evaluate if certain treatments can facilitate restoration within a reasonable timeframe of 5 to 10 years after treatment implementation.
The three restoration treatments used were:
i. “Harvest only”: Plots are being mowed two to three times a year and the biomass is removed.
ii. “Topsoil”: Topsoil was removed to a depth of 10 to 20 cm, depending on the depth of the O and A horizon, in four randomly selected areas within each of the eleven patches in late autumn 1995. The size of each topsoil removal area depended on individual patch size and was between 2700 and 7000 m2.
iii. “Topsoil+Propagules”: Propagules from target vegetation were added on half of the area where topsoil was removed, using fresh, seed-containing hay originating from a mixture of semi-dry to semi-wet species-rich grasslands of local provenance (within a radius of 7 km). Hay applications were conducted twice in 1995 and 1996. Repeated applications were chosen to account for the low quantity of available plant material per transfer, since area ratio between receptor and donor sites was roughly 1:1. In addition, hand-collected propagules from 15 selected target species of regional provenance (within a radius of 30 km) were equally applied in 1996 and 1997.
“Topsoil” and “Topsoil+Propagules” plots are mowed once a year, and the biomass is removed. Mowing on these plots started five years after the treatment was implemented.
Eleven permanent plots of 5 m x 5 m were randomly established in each treatment to monitor the vegetation development. The experiment was complemented with 11 control plots that represent the initial state of intensively managed grasslands, further referred to as “Initial”, and 11 control plots that represent the targeted state of donor sites for “Topsoil+Propagules”, further referred to as “Target”. Consequently, the experiment consists of 55 plots (5 treatments x 11 replicates). Management of intensively used grasslands includes mowing and fertilizing (manure) between two to five times a year, as well as different tillage regimes (no tillage for >50 years; last time of tillage <5 years).

 Nematode and plant sampling
Soil nematodes were sampled in 2 m x 2 m plots, randomly established at least 2 m away from the vegetation plots. We collected eight soil cores with a 2.2 cm diameter soil core sampler (Giddings Machine Company, Windsor, CO, USA) to a depth of 12 cm (representing the majority of the plant rooting system) in each plot at the beginning of July 2017. The eight cores within each replicate plot were combined, gently homogenized, placed in coolers and transported to the laboratory of NIOO in Wageningen, the Netherlands, within one week. Free-living nematodes were extracted from 200 g of fresh soil using Oostenbrink elutriator (Oostenbrink, 1960) and concentrated, resulting in 6 mL nematode solution. The nematode solution was subdivided into three subsamples, two for morphological identification and quantification, and one for molecular work (not used in this study). For morphological identification and quantification, nematodes were heat-killed at 90 °C and fixed in 4 % formaldehyde solution (final volume 10 mL per subsample). All nematodes in 1 mL of formaldehyde solution were counted, and a minimum of 150 individuals per 1 mL sample (or all if less nematodes were present) were identified to family level using Bongers (1988). We then extrapolated the numbers of each nematode taxa identified to the entire sample and expressed them per 100 g dry soil for further analyses.
We calculated number of nematode taxa and Shannon diversity and assessed nematode community composition. In addition, we classified the nematode taxa into feeding types (herbivores, bacterivores, fungivores, omni-carnivores), structural and functional guilds (Table S4). Structural guilds assign nematode taxa according to life-history traits into five colonizer-persister (C-P) classes, ranging from one (early colonizers of new resources) to five (persisters in undisturbed habitats; Bongers 1990). C-P classes can be categorized as indicators for nutrient-enriched (C-P1), stressed (C-P2) and structured (C-P3 + C-P4 + C-P5) soil conditions (Ferris, Bongers, & de Goede, 2001). Functional guilds assign nematode taxa according to their C-P classification combined with their feeding habits (Ferris, Bongers, & de Goede, 2001). Based on the structural and functional guild classification we calculated five additional indices to assess soil nutrient status, disturbance and food web characteristics using NINJA (Sieriebriennikov, Ferris, & de Goede, 2014). 1) The Maturity index indicates the degree of different environmental perturbations (e.g., tillage, nutrient enrichment, pollution) and is used to monitor colonization and subsequent succession after disturbances (Bongers, 1990). 2) The ratio between the Plant Parasite (C-P of herbivorous nematodes only) to Maturity index is used to monitor the recovery of disturbed habitats incorporating information of life-history traits for all feeding types (Bongers, van der Meulen, & Korthals, 1997). 3) The Enrichment index indicates nutrient-enriched soils and agricultural management practices (Ferris, Bongers, & de Goede, 2001). 4) The Structure index provides information about the succession stage of the soil food web and therefore correlates with the degree of maturity of an ecosystem (Ferris, Bongers, & de Goede, 2001). 5) The Channel index provides information about the predominant decomposition pathways, where higher values stand for a higher proportion of energy transformed through the slow fungal decomposition channel (Ferris, Bongers, & de Goede, 2001). In addition, the Structure and Enrichment indices can be displayed in a biplot where nematode assemblages are plotted along a structure (x-axis) and enrichment (y-axis) trajectory (increasing index values). Each biplot quadrat reflects different levels of disturbance, soil nutrient pools and decomposition pathways (Ferris, Bongers, & de Goede, 2001).
The plant surveys were conducted on the 25 m2 permanent plots in June 2017. Plant species cover was visually assessed according to the semi-quantitative cover-abundance scale of Braun-Blanquet (1964; nomenclature: Lauber & Wagner, 1996). We calculated number of species and Shannon diversity, and assessed plant community composition. We also counted the number of target species (all species recorded in the eleven target plots plus propagules of species applied by hand, resulting in a total of 143 species) and categorized plant species into species of concern based on their red list status in Switzerland as well as their protection status in Switzerland and the Canton Zurich (Moser, Gygax, Bäumler, Wyler, & Palese, 2002). Furthermore, we calculated indicator values for soil moisture and soil nutrients for each species according to Landolt et al. (2010).

 References
Baudirektion Kanton Zürich (2007). 10 Jahre Naturschutz-Gesamtkonzept für den Kanton Zürich 1995-2005 – Stand der Umsetzung. Zürich: Baudirektion Kanton Zürich.
Bongers, T. (1988). De nematoden van Nederland. Utrecht: Stichting Uitgeverij Koninklijke Nederlandse Natuurhistorische Vereniging.
Bongers, T. (1990). The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia, 83, 14-19. doi:10.1007/BF00324627
Bongers, T., van der Meulen, H., & Korthals, G. (1997). Inverse relationship between the nematode maturity index and plant parasite index under enriched nutrient conditions. Applied Soil Ecology, 6, 195-199. doi:10.1016/S0929-1393(96)00136-9
Braun-Blanquet, J. (1964). Pflanzensoziologie, Grundzüge der Vegetationskunde (3rd ed.). Wien: Springer.
Delarze, R., Gonseth, Y., Eggenberg, S., & Vust, M. (2015). Lebensräume der Schweiz: Ökologie - Gefährdung - Kennarten (3rd ed.). Bern: Ott.
Ferris, H., Bongers, T., & de Goede, R.G.M. (2001). A framework for soil food web diagnostics: extension of the nematode faunal analysis concept. Applied Soil Ecology, 18, 13-29. doi:10.1016/S0929-1393(01)00152-4
Landolt, E., Bäumler, B., Erhardt, A., Hegg, O., Klötzli, F., Lämmler, W., … Wohlgemuth, T. (2010). Flora indicativa. Ecological indicator values and biological attributes of the Flora of Switzerland and the Alps (2nd ed.). Bern: Haupt.
Lauber, K., & Wagner, G. (1996). Flora Helvetica. Flora der Schweiz. Bern: Haupt.
MeteoSchweiz (2018). Klimabulletin Jahr 2017, Zürich: MeteoSchweiz.
Moser, D., Gygax, A., Bäumler, B., Wyler, N., & Palese, R. (2002) Rote Liste der gefährteten Farn- und Blütenpflanzen der Schweiz. Bern: BUWAL.
Oostenbrink, M. (1960). Estimating nematode populations by some selected methods. In N.J. Sasser & W.R. Jenkins (Eds.), Nematology (pp. 85-101). Chapel Hill: University of North Carolina Press.
Sieriebriennikov, B., Ferris, H., & de Goede, R.G.M (2014). NINJA: An automated calculation system for nematode-based biological monitoring. European Journal of Soil Biology, 61, 90-93. doi:10.1016/j.ejsobi.2014.02.004

diversity

display_name: diversity

grassland

display_name: grassland

long-term-experiment

display_name: long-term-experiment

nematodes

display_name: nematodes

restoration

display_name: restoration

topsoil-removal

display_name: topsoil-removal

vegetation

display_name: vegetation

en: Grassland restoration: nematodes and plant communities


en: shortwave-radiation davos outgoing longwave-radiation incoming

en: Parameter file

en: DATA ACCESS

en: Processing doc

mimetype: application/pdf

en: Metadata

en:
Incoming and outgoing shortwave and longwave 2 min radiation measurements in Davos Dorf, CH.

 References
1. Marty, C., Philipona, R., Frohlich, C., Ohmura, A.. Altitude dependence of surface radiation fluxes and cloud forcing in the alps: results from the alpine surface radiation budget network. 2002. Theoretical and Applied Climatology. Volume 72. Issue 3-4. 137-155. http://dx.doi.org/10.1007/s007040200019. 10.1007/s007040200019.
2. Christoph Marty. Surface Radiation, Cloud Forcing and Greenhouse Effect in the Alps. 2000. Institute fuer Klimaforschung ETH. Zuercher Klima-Schriften. Volume 79. http://e-collection.library.ethz.ch/eserv/eth:23491/eth-23491-01.pdf.

davos

display_name: davos

incoming

display_name: incoming

longwave-radiation

display_name: longwave-radiation

outgoing

display_name: outgoing

shortwave-radiation

display_name: shortwave-radiation

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: ASRB_DAV: Shortwave and longwave radiation measurements (2 min) in Davos Dorf


automatic-detection-of-avalanches

maintainer: Matthias Heck

title_for_slug: automatic-detection-of-avalanches

en: monitoring avalanche-activity hidden-markov-model seismic snow

en: Results and Scripts

en: This zip package contains the classification results, coherence results, reference data set, python scripts used for the classification and the start and end times for the pre-processed continuous seismic data. A further readme file is also include for more detailed information.

en: This dataset contains the results obtained by an automatic classification using hidden Markov models of a continuous seismic dataset. To avoid long computational times, we reduced the seismic data using pre-processing step. The start and end times of the windows used for the classification are also included in this dataset. Furthermore, an avalanche reference data set is included and the python scripts used to perform the processing steps and the classification.

avalanche-activity

display_name: avalanche-activity

hidden-markov-model

display_name: hidden-markov-model

monitoring

display_name: monitoring

seismic

display_name: seismic

snow

display_name: snow

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: Automatic detection of avalanches


en: lake-level-regulation switzerland perialpine-lakes climate-change hydrologic-and-hydrodynamic-modelling

en: Lake outflow scenarios CH2018

en: Lake outflow scenarios for Lake Walen, Lake Zurich, Lake Brienz and Lake Thun.

en: Lake water level scenarios CH2018

en: Lake water level scenarios for Lake Walen, Lake Zurich, Lake Brienz and Lake Thun.

en: Lake characteristics

en: Lake characteristics: stage-area relation for Lake Walen, Lake Zurich, Lake Brienz and Lake Thun, stage-discharge relations for the unregulated Lake Walen and semi-regulated Lake Brienz and, lake level management rules for the two regulated lakes Zurich and Thun.

en:
The dataset "Lake_climate_change_scenarios_CH2018" provides simulation-based climate change impact scenarios for perialpine lakes in Switzerland. These transient future scenarios were produced by combining the hydrologic model PREVAH with the hydrodynamic model MIKE11 to simulate daily lake water level (Lake_water_level_scenarios_CH2018.xls) and outflow scenarios (Lake_outflow_scenarios_CH2018.xls) from 1981 to 2099, using the Swiss Climate Change Scenarios CH2018. The future scenarios contain a total of 39 model members for three Representative Concentration Pathways, RCP2.6 (concerted mitigation efforts), RCP4.5 (limited climate mitigation) and RCP8.5 (no climate mitigation measures). These scenarios result from the study titled "Lower summer lake levels in regulated perialpine lakes, caused by climate change," authored by Wechsler et al. in 2023.

The dataset emphasizes four specific Swiss lakes, each subject to different degrees of lake level management: an unregulated lake (Lake Walen), a semi-regulated lake (Lake Brienz), and two regulated lakes (Lake Zurich and Lake Thun).

In addition, the file (Lake_characteristics.xlsx) includes data used in the modeling process, encompassing the stage-area relation for the four lakes, stage-discharge relations for the unregulated and semi-regulated lakes, and lake level management rules for the two regulated lakes.

climate-change

display_name: climate-change

hydrologic-and-hydrodynamic-modelling

display_name: hydrologic-and-hydrodynamic-modelling

lake-level-regulation

display_name: lake-level-regulation

perialpine-lakes

display_name: perialpine-lakes

switzerland

display_name: switzerland

en: Lake climate change scenarios CH2018


maintainer: Anita C. Risch

en: mineralization grassland global-change ammonification climate nutrient-network soil nitrogen herbivore fertilization nitrification

en: Risch et al. 2020 - Raw data for publication "Global impacts of fertilization and herbivore removal on soil net nitrogen mineralization are modulated by local climate and soil properties"

en:
we assessed how the removal of mammalian herbivores (Fence) and fertilization with growth-limiting nutrients (N, P, K, plus nine essential macro- and micronutrients; NPK) individually, and in combination (NPK+Fence), affected potential and realized soil net Nmin across 22 natural and semi-natural grasslands on five continents. Our sites spanned a comprehensive range of climatic and edaphic conditions found across the grassland biome. We focused on grasslands, because they cover 40-50% of the ice-free land surface and provide vital ecosystem functions and services. They are particularly important for forage production and C sequestration. Worldwide, grasslands store approximately 20-30% of the Earth’s terrestrial C, most of it in the soil (Schimel, 1995; White et al., 2000).

Study sites and experimental design
The 22 sites contributing to this project are part of the Nutrient Network Global Research Cooperative (NutNet, https://nutnet.umn.edu/). Mean annual temperature across our 22 sites ranged from -4 to 22°C, mean annual precipitation from 252 to 1,592 mm, and elevations from 6 to 4,261 m above sea level (Fig. 1, Supplementary Table S1). Soil organic C varied from 0.8 to 7.8%, soil total N from 0.1 to 0.6%, and the soil C:N ratio from 9.1 to 21.5. Soil clay content spanned from 3.0 to 35%, and soil pH from 3.4 to 7.6 (Supplementary Table S2). Thus, the sites covered a wide range of environments in which grasslands occur (Fig. 1, Supplementary Table S1 & S2).
At each site, the effects of nutrient addition and herbivore removal were tested via a randomized-block design (Borer et al., 2014; Supplementary Fig. S1a). Three replicate blocks with 10 treatment plots each were established at each site, with the exception of the site at bldr.us, where only two blocks were established (Supplementary Fig. S1a). The 10 plots were randomly assigned to a nutrient or fencing treatment, but only a subset of four plots was used in the current study, each with a different treatment (see below; Supplementary Fig. S1a). All plots were 5 x 5 m and divided into four 2.5 x 2.5 m subplots (Supplementary Fig. S1b). 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; Supplementary Fig. S1b). Plots were separated by at least 1 m wide walkways.
In this study, we collected data from the following four treatments: (i) untreated control plots (Control), (ii) herbivore removal 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 removal (NPK+Fence; Supplementary Fig. S1a). The number of years of treatment differed among sites (2 – 9 years since start of treatment; Supplementary Table S1). For the nutrient additions, all sites applied 10 g N m-2 yr-1 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.
The vertebrate herbivore removal treatment (Fence) was established by fencing two plots, one control and one NPK plot, within each of the blocks (Supplementary Fig. S1a). We designed the fences so that they would effectively exclude aboveground mammalian herbivores with a body mass of over 50 g (Borer et al., 2014). At the majority of sites, the height of the fences was 180 cm, and the fence design included wire mesh (1 cm holes) on the first 90 cm along 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. While most sites only had native herbivores, a few sites (4) were also grazed by domestic animals (Supplementary Table S1).
Potential and realized soil net N mineralization, ammonification, nitrification and other soil properties
Each site participating in the study received a package containing identical material from the Swiss Federal Institute for Forest, Snow and Landscape Research (WSL) to be used for sampling and on-site N incubations. For the field incubation, we followed the protocol by Risch et al. (2015, 2019). Briefly, we drove a 5 x 15 cm (diameter x depth) steel cylinder 13.5 cm deep into the soil after clipping the vegetation at randomized locations in each plot. The top 1.5 cm of the cylinder remained empty to capture incoming N from run-off or deposition with a polyester mesh bag (mesh-size 250 µm) filled with 13.2 ± 0.9 g of acidic and alkaline exchanger resin (1:1 mixture; ion-exchanger I KA/ion-exchanger III AA, Merck AG, Darmstadt). The bag was fixed in place with a metal Seeger ring (Bruetsch-Rüegger Holding, Urdorf, Switzerland). We then removed 1.5 cm soil at the bottom of the cylinder and placed another resin bag to capture N leached from the soil column. The exchange resin was saturated with H+ and Cl- prior to filling the bags by stirring the mixture in 1.2M HCl for 1 h and then rinsing it with demineralized water until the electrical conductivity of the water reached 5 μS/cm. The cylinders were then re-inserted into the cored hole, level with the soil surface, and incubated for an average of 42 days (range 40 to 57days). The individual site coordinators chose the timing of incubation to start approximately six weeks prior to peak plant biomass production. All incubations were completed between February 2015 and January 2016 accounting for differences in growing season between northern and southern hemispheres. At the end of the incubation, the cylinders were re-collected and immediately shipped to WSL in an insulated box together with cold packs to halt further mineralization. Gloves were worn at all times to avoid contamination of the samples. Upon arrival at WSL, we extracted the resin bags and a 20 g subsample of sieved soil (4 mm) separately in a 100 ml PE-bottle with 80 ml 1 M KCl for 1.5 h on an end-over-end shaker and filtered through ashless folded filter paper (DF 5895 150, ALBET LabScience). We measured NO3- (colorimetrically; Norman & Stucki, 1981) and NH4+concentrations (flow injection analysis; FIAS 300, Perkin Elmer) on these filtrates.
At the start of the field incubation, we additionally collected two soil cores of 5 x 12 cm (diameter x depth) in each sampling plot and composited them to measure potential soil net Nmin, soil chemical and biological properties (see below). We also collected an additional sample (5 x 12 cm) to assess soil physical properties, which remained within the steel cylinder. Both ends were tightly closed with plastic caps. Cores were carefully packed to avoid further disturbance, and together with the composited soil samples, were shipped to the laboratory at WSL within a few days after collection.
From the composited samples, we extracted an equivalent of 20 g dry soil with KCl, as described above, and measured NO3- and NH4+concentrations. Realized soil net Nmin was then calculated as the difference between the inorganic N content of samples collected at the end of the incubation (plus N extracted from the bottom resin bag) and the N content at the beginning of the incubation; values were scaled to represent daily mineralization rates (mg N kg-1soil d-1; Risch et al., 2015). Realized soil net Nmin values represent an average period of 42 days prior to peak biomass, typically the highest period of biological activity, and not the entire year (Risch et al., 2019).
A second subsample of the composited sample was used to determine potential soil net Nmin in the laboratory (Risch et al., 2019). Briefly, we weighed duplicate samples (8 g dry soil) into 50-ml Falcon tubes. Soil moisture was brought to 60% of the field capacity of each plot, the Falcon tubes tightly closed and then incubated at 20°C for 42 days in a dark room. Every week the Falcon tubes were opened and ventilated. At the end of the incubation, the soil samples were extracted the same way as described above and NO3- and NH4+ was determined. Potential soil net Nmin was calculated as the difference between the N content before and after the incubation and scaled to represent daily values (mg N kg-1soil d-1). Using our NO3- and NH4+ measures we also calculated potential and realized soil net nitrification and soil net ammonification to be able to better understand the drivers of fertilization and herbivore removal effects on potential and realized soil net Nmin.
A third subsample of the composite soil sample was sieved (2 mm mesh) and microbial biomass (μg Cmic g-1 soil dry weight) was estimated 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 (Anderson & Domsch, 1978). The rest of the composited sample was dried at 65°C for 48 h, ground and sieved (2 mm mesh) to assess a series of soil chemical properties (Risch et al., 2019). We measured the percentage of clay as an indicator of soil texture (Gee & Bauder, 1986; Risch et al., 2019).
Statistical analyses
Potential and realized soil net Nmin were square root transformed to account for a highly skewed data distribution (yt = sign(y)*sqrt|y|; negative values in the data set impeded log transformation). To assess treatment effects on potential and realized soil net Nmin, we used linear mixed effects models (LMMs) fitted by maximum likelihood using the lme function from the nlme package (version 3.131.1; Pinheiro et al., 2016), R version 3.6.1; R Foundation for Statistical Computing. Treatment (Control, Fence, NPK, NPK+Fence) was a fixed factor, with site and block as random factors, where block was nested within site. We also tested for effects of time since start of treatments in preliminary analyses by adding total treatment years as an additional fixed factor. We did not find a significant effect of years of treatment, and thus dropped this variable from the models. The LMMs were corrected using varIdent if the homogeneity of variance criterion was not met (Pinheiro et al., 2016). To visualize our results, we calculated treatment effects using Cohens’ d statistic (Cohen, 1977; Koricheva et al., 2013). Note that calculating response ratios (or log response ratios) was not possible with our data, as we have both negative and positive values. We also fitted LMMs for potential and realized soil net ammonification and nitrification to gain more insight into how global change affects the processes underpinning potential and realized soil net Nmin. We also sqrt-transformed (see above) these dependent variables. Treatment was included as a fixed factor with random factors as described above. In addition, we assessed how potential and realized soil net Nmin were related to potential and realized soil net ammonification and nitrification, respectively. For this, we calculated site averages for each treatment separately. We then ran LMMs, with potential and realized soil net Nmin as the dependent variable, potential and realized soil net ammonification/nitrification as the independent ones. Site was included as a random factor.
Based on our previous work (Risch et al., 2019) and the existing literature (Schimel & Bennett, 2004; Liu et al., 2017), we developed a priori causal conceptual models of relationships among treatments, environmental drivers, and potential and realized soil net Nmin (Supplementary Fig. S2) to test with structural equation modelling (SEM) using a d-sep approach (Shipley, 2009; Lefcheck, 2016). The variables included in the model were long-term climatic conditions, specifically, site-level mean annual precipitation (MAP) and temperature of the wettest quarter (T.q.wet), plot-level soil texture (clay content) and soil microbial biomass. Mean annual precipitation and T.q.wet were obtained from WorldClim (Hijmans et al., 2005) (http://www.worldclim.org/) and together with the experimental treatments were predicted to directly affect soil properties and soil net Nmin (Supplementary Fig. S2). Soil clay content was, in turn, predicted to affect microbial biomass and soil net Nmin. Because we determined microbial biomass prior to incubating the samples in the laboratory or field, we assumed that the abundance of these microbes would be responsible for N process rates and not vice versa (Supplementary Fig. S2). We tested our conceptual model (Supplementary Fig. S2) 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 (Shipley, 2009; Lefcheck, 2016). 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 covariates and experimental treatments detected in the SEMs, we calculated treatment effect sizes, i.e. the differences in potential or realized soil net Nmin between Control and treatments (Fence, NPK, NPK+Fence) and plotted these values against the climate or soil covariates. Finally, we fitted LMMs for the soil variables included in our SEMs, with treatment as the fixed factor, and with site and block as random factors, where block was nested within site.

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

Paper Citation:

Risch Anita C., Zimmermann, Stefan, Moser, Barbara, Schütz, Martin, Hagedorn, Frank, Firn, Jennifer, Fay, Philip A., Adler, Peter B., Biederman, Lori A., Blair, John M., Borer, Elizabeth T., Broadbent, Arthur A.D., Brown, Cynthia S., Cadotte, Marc W., Caldeira, Maria C., Davies, Kendi F., di Virgilio, Augustina, Eisenhauer, Nico, Eskelinen, Anu, Knops, Johannes M.H., MacDougall, Andrew S., McCulley, Rebecca L., Melbourne, Brett A., Moore, Joslin L., Power, Sally A., Prober, Suzanne M., Seabloom, Eric W., Siebert, Julia, Silveira, Maria L. , Speziale, Karina L., Stevens, Carly J., Tognetti, Pedro M., Virtanen, Risto, Yahdjian, Laura, Ochoa-Hueso, Raul (accepted). Global impacts of fertilization and herbivore removal on soil net nitrogen mineralization are modulated by local climate and soil properties. Global Change Biology

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

We assessed how the removal of mammalian herbivores (Fence) and fertilization with growth-limiting nutrients (N, P, K, plus nine essential macro- and micronutrients; NPK) individually, and in combination (NPK+Fence), affected potential and realized soil net Nmin across 22 natural and semi-natural grasslands on five continents. Our sites spanned a comprehensive range of climatic and edaphic conditions found across the grassland biome. We focused on grasslands, because they cover 40-50% of the ice-free land surface and provide vital ecosystem functions and services. They are particularly important for forage production and C sequestration. Worldwide, grasslands store approximately 20-30% of the Earth’s terrestrial C, most of it in the soil (Schimel, 1995; White et al., 2000).

ammonification

display_name: ammonification

climate

display_name: climate

fertilization

display_name: fertilization

global-change

display_name: global-change

grassland

display_name: grassland

herbivore

display_name: herbivore

mineralization

display_name: mineralization

nitrification

display_name: nitrification

nitrogen

display_name: nitrogen

nutrient-network

display_name: nutrient-network

soil

display_name: soil

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Anthropogenic change and soil net N mineralization