maintainer: Anita C. Risch

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

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

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: 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: Anthropogenic change and soil net N mineralization


maintainer: Matthias Haeni

en: meteorology relative-humidity meteo par precipitation climate-change wind-speed swiss-forest-lab temperature long-term

state: active

type: dataset

en: Site aerial photo (Nationalpark)

en: Site aerial photo from both Nationalpark stations. Source: Google. 2016. Site aerial photo from both Nationalpark stations. Retrieved November 11, 2016, from https://www.google.ch/maps

format: SERVICE

mimetype: image/jpeg

en: Forest station data (NAB)

media_type: text/csv

en: Data file.

format: CSV

mimetype: text/csv

en: Data plot of the forest station (NAB)

media_type: application/pdf

en: Data overview for air temperature, relative humidity, precipitation, PAR, wind speed.

format: PDF

mimetype: application/pdf

en: Site photo of the field station (NAF)

en: Site photo of the field station form the side.

en: Field station data (NAF)

en: Site photo of the forest station (NAB)

en: Site photo of the forest station from the side.

en: Data plot of the field station (NAF)

en: High quality meteorological data are needed for long-term forest ecosystem research, particularly in the light of global change. The long-term data series published here comprises almost 20 years of measurements for two meteorological stations in Nationalpark in Switzerland where one station is located within a natural coniferous forest stand (NAB) with mountain pine (_Pinus mugo_; 210 yrs) as dominant tree species. A second station is situated in the very vicinity outside of the forest (field station, NAF). The meteorological time series are presented in hourly time resolution of air temperature, relative humidity, precipitation, photosynthetically active radiation (PAR) and wind speed. Nationalpark is part of the Long-term Forest Ecosystem Research Programme (LWF) established and maintained by the Swiss Federal Research Institute WSL.

climate-change

display_name: climate-change

long-term

display_name: long-term

meteo

display_name: meteo

meteorology

display_name: meteorology

par

display_name: par

precipitation

display_name: precipitation

relative-humidity

display_name: relative-humidity

swiss-forest-lab

display_name: swiss-forest-lab

temperature

display_name: temperature

wind-speed

display_name: wind-speed

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

it: Istruzione, cultura e sport

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

it: EnviDat: il portale dei dati ambientali

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

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

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

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester

en: Nationalpark, Switzerland: Long-term forest meteorological data from the Long-term Forest Ecosystem Research Programme (LWF), from 1997 onwards


maintainer: Ionuț Iosifescu Enescu

en: map 0.2GB

media_type: image/tiff

format: TIFF

mimetype: image/tiff

en: map preview 0.5MB

en: large map 10GB

en: Bulk Data Access

en: Link to the S3 bucket where all the data resources can be downloaded in bulk

en: map 1GB

en: huge map 100GB

en: map 2GB

en: map 0.5GB

en:
The attached data are some large GIS raster files (GeoTIFFs) made with Natural Earth data. Natural Earth is a free vector and raster map data @ naturalearthdata.com. The data used for creating these large files was the "Cross Blended Hypso with Shaded Relief and Water". Data was concatenated to achieve larger and larger files. Internal pyramids were created, in order that the files can be opened easily in a GIS software such as QGIS or by a (future) GIS data visualisation module integrated in EnviDat.

Made with Natural Earth. Free vector and raster map data @ naturalearthdata.com

gis

display_name: gis

large-files

display_name: large-files

natural-earth

display_name: natural-earth

qgis

display_name: qgis

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: Large GIS raster data derived from Natural Earth Data (Cross Blended Hypso with Shaded Relief and Water)


maintainer: Ionut Iosifescu Enescu

en: hpc open-data research-data-management open-science open-source jupyter-notebooks reproducibility

en: Open Science Support @ WSL (Poster)

en:
This poster was originally created for the swissuniversities Open Science Action Plan: Kick-Off Forum, and showed to the audience on 17.10.2019. It illustrates how the environmental data portal EnviDat provides the tools for fostering Open Science and Reproducibility of scientific research at WSL.

Supporting open science is a highly relevant user requirement for EnviDat and for implementing FAIR (Findability, Accessibility, Interoperability and Reusability) principles at dataset level.
EnviDat encourages WSL scientists to complement data publication with a complete description of research methods and the inclusion of the open source software, code or scripts used for processing the dataset or for obtaining the published results. By openly publishing open software (e.g. as Jupyter notebooks) alongside research data sets, researchers can contribute to mitigate reproducibility issues.

EnviDat also promotes and supports, where possible and practical, the publication of software as Jupyter notebooks. Jupyter notebooks provide a solution for improved documentation and interactive execution of open code in a wide range of programming languages (Python, R, Octave/Matlab, Java or Scala). These programming languages are widely used in environmental research at WSL and well supported by the Jupyter-compatible kernels. We have sucessfully interfaced EnviDat-hosted notebooks with the WSL High-Performance Computing (HPC) Linux Cluster through a JupyterHub/JuypterLab beta installation on the HPC cluster implemented in close collaboration with the WSL IT-Services.
For existing software that cannot be easily migrated to Jupyter Notebooks, the Open Science and Reproducibility is assisted by containerisation. We have proven that several Singularity containers can successfully run on WSL's HPC cluster.

Finally, the researchers can upload the data/results complemented by code (e.g. as Jupyter Notebooks, or Singularity containers) and any additional documentation in EnviDat. Consequently, they will receive a DOI for the entire dataset, which they can reference in their science paper in order to publish a more reproducible research.

_License_: This poster is released by WSL and the EnviDat team to the public domain under a Creative Commons 4.0 CC0 "No Rights Reserved" international license. You can reuse this poster in any way you want, for any purposes and without restrictions.

hpc

display_name: hpc

jupyter-notebooks

display_name: jupyter-notebooks

open-data

display_name: open-data

open-science

display_name: open-science

open-source

display_name: open-source

reproducibility

display_name: reproducibility

research-data-management

display_name: research-data-management

en: Open Science Support at the Swiss Federal Research Institute WSL. The EnviDat Concept


maintainer: EnviDat Support

en: emma stable-isotopes tracers end-member-mixing-analysis eemma mixing-models

en: eemma vignette script

en: An R script for the vignettes (simple data analysis demonstrations) shown in the user's guide

format: OCTET

mimetype: application/octet-stream

en: An R script that performs ensemble end-member mixing analysis

en: eemma user's guide

en: A brief user's guide to EEMMA, with illustrations using a simple data set

en: vignette data

media_type: text/plain

en: A simple data set for the vignette script

format: TXT

mimetype: text/plain

en:
The R script eemma.R, which implements Ensemble End-Member Mixing Analysis (EEMMA) to estimate source fractions in mixtures, exploiting information contained in time-series correlations among tracer time series.

A brief user's guide, a demonstration script, and a demonstration data set are also provided, to accompany Kirchner, J.W., Mixing models with multiple, overlapping, or incomplete end-members, quantified using time series of a single tracer, Geophysical Research Letters, 2023.

The user's guide is available for public use under Creative Commons CC-BY-SA. Public use of the scripts is permitted under GNU General Public License 3 (GPL3); for details see https://www.gnu.org/licenses/

eemma

display_name: eemma

emma

display_name: emma

end-member-mixing-analysis

display_name: end-member-mixing-analysis

mixing-models

display_name: mixing-models

stable-isotopes

display_name: stable-isotopes

tracers

display_name: tracers

en: eemma.R, an R script for Ensemble End-Member Mixing Analysis


en: environmental-dna edna fish-composition teleo coral-reef marine indian-ocean laperouse

en: Data access

en: Please use cyberduck to download for the moment

en:
Fish environmental DNA data set collected in the shallow seamount of La Pérouse located between Madagascar and Réunion Island

In 2019, we sampled the shallow seamount of La Pérouse, located between Madagascar and Réunion Island in the WIO. La Pérouse is an extinct volcano with a summit depth at 60 m below the sea surface. The plateau of the summit has a maximum length of 12 km and a maximum width of 4 km. One side of the seamount might have collapsed in the past, leading to a less common crescent-shaped summit. The abyssal plains that surround this isolated pinnacle are at a depth of 5000 m. Despite its proximity to Réunion Island (160 km northwest), La Pérouse and its surroundings are still poorly studied.
To sample eDNA in this location we used a submersible pump (Subspace, Geneva, Switzerland; nominal flow of ~ 1.0 L/min) to collect samples at various depths during close circuit rebreather dives. The pump has an internal battery and can be activated underwater allowing in situ filtration directly at the targeted habitat and depth. The filtration lasted for 30 min to collect a total water volume of ~ 30 L. During close circuit rebreather dives, water samples from various depths, ranging from 60 to 140 m below the sea surface, were collected using only the Submersible method (eight stations: no replicates). Out of the eight water samples, six were collected during the day, while two were collected at night (IS2_1 at 62 m and IS7_1 at 100–125 m). We applied this filtration protocols to collect and filter water throughout a VigiDNA® 0.2 μM cross-flow filtration capsule (SPYGEN, le Bourget du Lac, France) using disposable sterile tubing for each filtration capsule. After the filtration process, the remaining water in the capsules was emptied. The capsules were then filled with 80 ml of lysis conservation buffer (CL1 buffer SPYGEN, le Bourget du Lac, France) and stored at room temperature in the dark. At La Pérouse, a total of eight water samples were collected during the MONT LA PÉROUSE expedition (https://gombessa-expeditions.com/mont-la-perouse/) on board the ship La Curieuse to the seamount La Pérouse from October 27th to November 5th 2019. We followed a strict contamination control protocol in both field and laboratory stages. Each water sample processing included the use of disposable gloves and single-use filtration equipment to avoid any risk of contamination. Libraries were prepared with ligation using the MetaFast protocol (Fasteris).
Data content:
	* rawdata/: contains the raw reads for each individual sample. One archive contains the paired-end reads specified by the _R1 or _R2 suffix as well as individually tagged PCR replicates (if available) together with an archive containing all extraction and PCR blank samples of the library. Reads have been demultiplexed using cutadapt but not trimmed, individual demultiplexing tags and primers remain present in the sequences.
	* taxadata/: contains the table with all detected taxonomy for each sample after bioinformatic processing (see Polanco et al. 2020 for details; https://doi.org/10.1002/edn3.140) and associated field metadata.
	* metadata/: contains two metadata files, one related to the data collected in the field for each filter, and the second related to the sequencing process in the lab (including the tag sequence, library name, and marker information for each sample)

coral-reef

display_name: coral-reef

edna

display_name: edna

environmental-dna

display_name: environmental-dna

fish-composition

display_name: fish-composition

indian-ocean

display_name: indian-ocean

laperouse

display_name: laperouse

marine

display_name: marine

teleo

display_name: teleo

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Environmental DNA Marine France LaPerouse 2019


en: biogeochemical-cycling yellowstone wolf soil-biodiversity detritus bacteria carrion carcass fungi decomposition

en:
Site	ID for different carcasses across the Yellowstone Landscape
Location	Sampling location, either under carcass or control plot
Species	Indicating whether carcass was from bison or elk
Gender	Gener of carcass; F = female, M = male, ? = unknown
Lat	Latitude
Long	Longitude
DSK	Days since kill; in days
Soil.N	Soil N concentration (%)
Soil.C	Soil C concentration (%)
Veg.C	Vegetation C concentration (%); dominant grasses sampled
Veg.N	Vegetation N concentration (%); dominant grasses sampled
Veg.CtoN	Vegetation C to N ratio; dominant grasses sampled
Soil.resp	Soil respiration (g m-2 h-1)
soil.temp	Soil temperature (°C)
Soil.moist	Soil moisture (volumentric; %)
CRR	Cotton rotting rate
bacteria.richness	Bacteria richness
bacteria.diversity	Bacteria diverstiy
fungi.richness	Fungi richness
fungi.diversity	Fungi diversity
bacteria.abundance	Bacteria abundance
fungi.abundance	Fungi abundance

format: XLS

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

Paper Citation:

Risch, A.C., Frossard, A., Schütz, M., Frey, B., Morris, A.W., Bump, J.K. (accepted) Effects of elk and bison carcasses on soil microbial communities and ecosystem functions in Yellowstone, USA. (accepted). Functional Ecology
doi: ...

Methods
Study area and study sites
This study was conducted in YNP’s Northern Range (NR), located in north-western Wyoming and south-western Montana, USA (~44.9163° N, 110.4169° W). The NR expands over ~1000 km2 and features long cold winters and short dry summers. Grasslands and shrublands dominate the NR that is the home of large migratory herds of bison (winter counts 2017: ~3919 individuals; Geremia, Wallen, & White, 2017) and elk (~5349 individuals) as well as their main predators, approximately five packs of wolves with a total of 33 individuals (Smith et al., 2017). As part of a long-term research program within YNP, wolf predation has been studied since their reintroduction in 1995.
For our study, we received ground-truthed coordinates of bison and elk carcasses from winter 2016/17 (November 2016 through April 2017) from the YNP Wolf Project. Between June 20 and July 1, 2017, we visited 24 carcasses in total. At five sites, we could not sample as the carcasses were no longer found. In total we located remains (hairmats, rumen content, bones, teeth) of 19 adult male and female carcasses (7 bison, 12 elk; Supplementary Table 1). Live body weights of adult bison and elk are approximately 730 kg (male bison), 450 kg (female bison), 330 kg (male elk), and 235 kg (female elk, Meagher, 1973; Quimby & Johnson, 1951).
The kills and subsequent consumption happened between 34 and 173 days prior to our sampling (hereafter “days since kill”, DSK), for which we accounted in our statistics. Note that wolves and other scavengers consumed the soft tissue of the carcasses quickly, hence, there is close to no soft tissue left for decomposition as compared to an intact body left on the soil surface. The 19 carcass sites covered the extent of YNP’s NR, with both bison and elk carcasses showing similar distributions; elevation ranged from 1703 to 2884 m a.s.l. (Supplementary Fig 1 & Supplementary Table 1). The carcasses were all located in grassland or sage-brush shrubland, with or without sparsely scattered trees, and both bison and elk carcasses showed the same distribution of DSK. At each study site, we selected a reference plot (hereafter “control”) that was of comparable size, slope aspect and vegetation to the carcass location (hereafter “carcass”). The control was at least 10 m away (Danell, Berteaux, & Brathen, 2002; Melis et al., 2007) from the carcass itself to ensure the absence of potential direct and indirect carcass effects (paired design; (Bump, Webster, et al., 2009; Bump, Peterson, et al., 2009).

Ecosystem functions and soil properties
We randomly collected 50 g of mineral soil from three locations on both control and carcass plots to a depth of 5 cm with sterile techniques and gently mixed the material to obtain a composite sample. Half the soil sample was immediately bagged in plastic bags (whirl packs), stored in a cooler with ice packs (~5 ºC), sieved (2-mm) and frozen within 4-6 hours of collection to assess soil microbial communities. For this purpose, we extracted total genomic DNA from 0.5 g soil using the PowerSoil DNA Isolation Kit (Qiagen, Hilden, Germany). DNA concentrations were measured using PicoGreen (Molecular Probes, Eugene, OR, USA). PCR amplifications of partial bacterial small-subunit ribosomal RNA genes (region V3–V4 of 16S rRNA) and fungal ribosomal internal transcribed spacers (region ITS2) were performed as described previously (Frey et al., 2016). Each sample consisting of 40 ng DNA was amplified in triplicate and pooled before purification with Agencourt AMPure XP beads (Beckman Colter, Berea, CA, USA) and quantified with the Qubit 2.0 fluorometric system (Life Technologies, Paisley, UK). Amplicons were sent to the Genome Quebec Innovation Center (Montreal, Canada) for barcoding using the Fluidigm Access Array technology and paired-end sequencing on the Illumina MiSeq v3 platform (Illumina Inc., San Diego, CA, USA).
Quality control of bacterial and fungal reads was performed using a customized pipeline (Supplementary Table 2; Frey et al., 2016). Paired-ends reads were matched with USEARCH (Edgar & Flyvbjerg, 2015), substitution errors were corrected using Bayeshammer (Nikolenko, Korobeynikov, & Alekseyev, 2013) and PCR primers were trimmed (allowing for 1 mismatch, read length >300 bp for 16S and >200 bp for ITS primers) using Cutadapt (M. Martin, 2011). Sequences were dereplicated and singleton reads removed prior to clustering into operational taxonomic units (OTUs) at 97% identity using USEARCH (Edgar, 2013). The remaining centroid sequences were tested for the presence of ribosomal signatures using Metaxa2 (Bengtsson-Palme et al., 2015) or ITSx (Bengtsson-Palme et al., 2013). Taxonomic assignments of the OTUs were obtained using Bayesian classifier (Wang, Garrity, Tiedje, & Cole, 2007) with a minimum bootstrap support of 60% implemented in mothur (Schloss et al., 2009) by querying the bacterial and fungal reads against the SILVA Release 128 (Quast et al., 2013) and UNITE 8.0 (Abarenkov et al., 2010) reference databases for 16S and ITS OTUs, respectively.
Abundances of the bacterial 16S rRNA gene and fungal ITS amplicon were determined by quantitative real-time PCR (qPCR) on an ABI7500 Fast Real-Time PCR system (Applied Biosystems, Foster City, CA, USA) as described previously (Frossard et al., 2018). The same primers (without barcodes) and cycling conditions as for the sequencing approach were used for the 16S and ITS qPCR. Three standard curves per target region were obtained using tenfold serial dilutions of plasmids generated from cloned targets (Frey, Niklaus, Kremer, Lüscher, & Zimmermann, 2011). Data were converted to represent mean copy number of targets per gram of soil (dry weight).
The other half of the soil sample was bagged in paper, dried to constant weight at 60°C, passed through a 2 mm sieve and analyzed for total C and N concentration with a CE Instruments NC 2100 soil analyzer (CE Elantech Inc., Lakewood NJ, USA). We also collected 20 mature and undamaged leaves of the dominant grass species growing on control and carcass sites, but taxa were not recorded. The plant material was dried at 60°C, finely ground till homogenized and also analyzed to obtain total C and N concentrations. Soil temperature (10 cm depth) was measured with a waterproof digital thermometer (Barnstead International, Dubuque IA, USA) at three locations each at the control and carcass site. Soil moisture (0 – 10 cm depth) was measured with time domain reflectometry (Field-Scout TDR-100; Spectrum Technologies, Plainfield IL, USA) at five randomly chosen points on control and carcass sites. We measured soil respiration at five randomly chosen points at both control and carcass sites with a PP-Systems SRC-1 soil respiration chamber (closed circuit) attached to a PP-Systems EGM-4 infrared gas analyzer (PP-Systems, Amesbury, MA, USA). For each measurement the soil chamber (15 cm high; 10 cm diameter) was tightly placed on the soil surface, after clipping plants to avoid measuring plant respiration or photosynthesis. Measurements were conducted over 120 s.
In addition, we assessed the decomposition rates of standardized OM using the cotton strip assay (Latter & Howson, 1977; Latter & Walton, 1988). Cotton cloth tensile strength loss (CTSL) is a measure of decomposition, and an index to express the combined effect of soil microclimatic, physical, chemical and biological properties on decomposition while accounting for OM quality (Latter & Walton, 1988; Risch, Jurgensen, & Frank, 2007; Withington & Sanford Jr., 2007). We placed five 20 cm wide x 13 cm long sheets of 100% unbleached cotton cloth (American Type SM 1/18’’, Warp: 34/1, Weft: 20/1, Weave plain, 29.5 picks/cm warp, 22 picks/cm weft, 237 g/m2; Daniel Jenny & Co., Switzerland;) at each carcass and control site vertically into the soil by making slits with a flat spade to a depth of 12 cm. We inserted each cloth with the spade, and then pushed the slit closed to assure tight contact with the soil. The cloths were retrieved after 18 to 27 days. After retrieval, the cloths were air-dried, remaining soil gently removed by hand, and 1.5 cm wide strips were cut at the 3.5-5.0 cm (top) and the 9-10.5 cm (bottom) soil depth. The strips were equilibrated at 50 % relative humidity and 20°C for 48 hours (climate chamber) prior to strength testing (Scanpro Awetron TH-1 tensile strength tester; AB Lorentzen and Wettre, Kista, Sweden). Cotton rotting rate (CRR) = (CTScontrol - CTSfinal/CTSfinal)1/3 * (365/t), where CTScontrol is the cotton tensile strength of a control cloth and CTSfinal the cotton tensile strength of the incubated sample, t is the incubation period in days. Control cloths were inserted into the ground and immediately retrieved to account for tensile strength loss associated with cloth insertion. We averaged the CRR of top and bottom strips for further analyses as no difference was found between the two. All sampling and cloth insertion took place between June 20 and July 1, 2017, cloths were retrieved between July 17 and 20, 2017. Soil respiration, average CRR, vegetation N concentration and vegetation C:N ratio are defined as ecosystem functions, soil C and N concentration, soil temperature and moisture as soil abiotic properties, and bacterial and fungal richness (number of taxa), diversity (Shannon) and abundance as soil biotic properties.

Statistical analyses
Univariate analyses for ecosystem functions, soil biotic and abiotic properties
We tested whether individual ecosystem functions, soil biotic and abiotic properties differed between carcass and control (“Location”), bison and elk (“Species”) and days since kill (“DSK”). For this purpose, we used linear mixed effect models (LMM, “nlme” package v 3.1 – 131.1 in R v 3.4.4; Pinheiro, Bates, DebRoy, & Sarkar, 2018; R Core Team, 2019) with Location, Species, Location x Species and DSK as fixed effects. Site was included as random effect to account for the paired design. We developed a separate model for all dependent variables. All but bacterial richness, fungal richness, fungal diversity and vegetation N concentration were natural-log transformed to meet model assumptions. For each LMM, we calculated contrasts to assess the specific comparisons we were interested in with the “lsmeans” package v 2.27-62 (Lenth & Love, 2018): 1) carcass vs control, 2) carcass bison vs control bison, and 3) carcass elk vs control elk. We also tested whether we had differences between bison and elk carcasses or the sites where bison and elk were killed and included contrasts 4) carcass bison vs carcass elk and 5) control bison vs control elk.
We calculated the log response ratio (LRR = ln[carcass/control]) to obtain carcass effects for all variables for both species separately. LRR < 0 indicates higher value at control compared to carcass, LRR > 0 indicates higher values at carcass compared control. We used LRRs for visualization and to assess spatial patterns in carcass effects across YNP. For this purpose we calculated the Moran’s I statistic for each ecosystem function, soil biotic and abiotic property based on a latitude-longitude matrix with the “moran.test” function in the “spdep” package version 1.1-3 (Bivand et al., 2019).
Multivariate analyses
Rare OTUs, defined as OTUs with a low abundance of reads, were retained in multivariate methods because they only marginally influence these analyses (Gobet, Quince, & Ramette, 2010). Bray–Curtis dissimilarity matrices were generated based on square-root-transformed matrices. We used Principal Coordinate Analyses (PCoA) to assess how soil bacterial and fungal communities differed between control and carcass of bison and elk (“vegan” package v 2.5-4, Oksanen et al., 2019). We then extracted PCoA axes scores 1 and 2 and used LMM (“nlme” package) with Location, Species, Location x Species and DSK as fixed effects. Site was, again, included as random effect. We again calculated the contrasts as described above using the “lsmeans” package. We also assessed how ecosystem functions, and soil abiotic and biotic properties were related to the soil bacteria and fungi community structure associated with bison and elk control and carcasses using the “envfit” function in the “vegan” package (Oksanen et al., 2019).
Indicator species analyses were performed using the multipatt function implemented in the “indicspecies” package version 1.7.6 with 100000 permutations (De Caceres & Jansen, 2016). This step allowed to identify OTUs that led to changes in multivariate patterns between control and carcass of both bison and elk separately (De Cáceres, Legendre, & Moretti, 2010). The multipatt function uses a point biserial correlation coefficient statistical test. Indicator OTUs were defined as bacterial and fungal OTUs with more than 50 sequences, i.e., removing rare taxa and taxa with low abundances containing little indicator information (Rime et al., 2015) and that were significantly correlated with Location (p < 0.05, correlation coefficient > 0.3). A heatmap of these OTUs were generated with the vegan and ggplot2 packages. The indicator analyses were performed in R version 3.3.3 (R Core Team, 2017).

References
Abarenkov, K., Henrik Nilsson, R., Larsson, K.-H., Alexander, I. J., Eberhardt, U., Erland, S., … Kõljalg, U. (2010). The UNITE database for molecular identification of fungi – recent updates and future perspectives. New Phytologist, 186(2), 281–285. doi:10.1111/j.1469-8137.2009.03160.x
Bengtsson-Palme, J., Hartmann, M., Eriksson, K. M., Pal, C., Thorell, K., Larsson, D. G. J., & Nilsson, R. H. (2015). metaxa2: improved identification and taxonomic classification of small and large subunit rRNA in metagenomic data. Molecular Ecology Resources, 15(6), 1403–1414. doi:10.1111/1755-0998.12399
Bengtsson-Palme, J., Ryberg, M., Hartmann, M., Branco, S., Wang, Z., Godhe, A., … Nilsson, R. H. (2013). Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods in Ecology and Evolution, 4(10), 914–919. doi:10.1111/2041-210X.12073
Bivand, R., Altman, M., Anselin, L., Assuncao, R., Berke, O., Blanchet, G., … Yu, D. (2019). spdep: Spatial dependence, weighthing schemes, statistics. R package version 1.1-3.
Bump, J. K., Peterson, R. O., & Vucetich, J. A. (2009). Wolves modulate soil nutrient heterogeneity and foliar nitrogen by configuring the distribution of ungulate carcasses. Ecology, 90(11), 3159–3167.
Bump, J. K., Webster, C. R., Vucetich, J. A., Peterson, R. O., Shields, J. M., & Powers, M. D. (2009). Ungulate carcasses perforate ecological filters and create biogeochemical hotspots in forest herbaceous layers allowing trees a competitive advantage. Ecosystems, 12(6), 996–1007. doi:10.1007/s10021-009-9274-0
Danell, K., Berteaux, D., & Brathen, K. A. (2002). Effect of muskox carcasses on nitrogen concentration in tundra vegetation. Arctic, 55(4), 389392.
De Caceres, M., & Jansen, F. (2016). indicspecies: relationship between species and groups of species. R package version 1.7.6.
De Cáceres, M., Legendre, P., & Moretti, M. (2010). Improving indicator species analysis by combining groups of sites. Oikos, 119(10), 1674–1684. doi:10.1111/j.1600-0706.2010.18334.x
Edgar, R. C. (2013). UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature Methods, 10, 996.
Edgar, R. C., & Flyvbjerg, H. (2015). Error filtering, pair assembly and error correction for next-generation sequencing reads. Bioinformatics, 31(21), 3476–3482. doi:10.1093/bioinformatics/btv401
Frey, B., Niklaus, P. A., Kremer, J., Lüscher, P., & Zimmermann, S. (2011). Heavy-machinery traffic impacts methane emissions as well as methanogen abundance and community structure in oxic forest soils. Applied and Environmental Microbiology, 77(17), 6060–6068. doi:10.1128/AEM.05206-11
Frey, B., Rime, T., Phillips, M., Stierli, B., Hajdas, I., Widmer, F., & Hartmann, M. (2016). Microbial diversity in European alpine permafrost and active layers. FEMS Microbial Ecology, 92(3), fiw018.
Frossard, A., Donhauser, J., Mestrot, A., Gygax, S., Bååth, E., & Frey, B. (2018). Long- and short-term effects of mercury pollution on the soil microbiome. Soil Biology and Biochemistry, 120, 191–199. doi:https://doi.org/10.1016/j.soilbio.2018.01.028
Geremia, C., Wallen, R., & White, P. J. (2017). Status report of the Yellowstone bison population, September 2017. Yellowstone National Park, Mammoth, WY, USA: National Park Service, Yellowstone Center for Resources.
Gobet, A., Quince, C., & Ramette, A. (2010). Multivariate cutoff level analysis (MultiCoLA) of large community data sets. Nucleic Acids Research, 38(15), e155–e155. doi:10.1093/nar/gkq545
Latter, P., & Howson, G. (1977). The use of cotton strips to indicate cellulose decomposition in the field. Pedobiologia, (17), 145–155.
Latter, P., & Walton, D. (1988). The cotton strip assay for cellulose decomposition studies in soil: history of the assay and development. In Cotton strip assay: an index for decomposition in soils (pp. 7–9). ITE Symposium, Institute of Terrestrial Ecology, Natural Environment Research Council, UK.
Lenth, R., & Love, J. (2018). lsmeans: least-squares means. R package version 2.27-62.
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet.Journal, 17(1), 10–12.
Meagher, M. M. (1973). The bison of Yellowstone National Park. NPS Scientific Monograph (Vol. 1). National Park Service, Yellowstone Center for Resources.
Melis, C., Selva, N., Teurlings, I., Skarpe, C., Linnell, J. D. C., & Andersen, R. (2007). Soil and vegetation nutrient response to bison carcasses in Białowieża Primeval Forest, Poland. Ecological Research, 22(5), 807–813. doi:10.1007/s11284-006-0321-4
Nikolenko, S. I., Korobeynikov, A. I., & Alekseyev, M. A. (2013). BayesHammer: Bayesian clustering for error correction in single-cell sequencing. BMC Genomics, 14(1), S7. doi:10.1186/1471-2164-14-S1-S7
Oksanen, J., Blanchet, F. G., Friendly, M., Kindt, R., Legendre, P., McGlinn, D., … Wagner, H. H. (2019). vegan: community ecology package. R package version 2.5-4.
Pinheiro, J., Bates, D., DebRoy, S., & Sarkar, D. (2018). nlme: Linear and nonlinear mixed effect models. R package version 3.1-131.1.
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., … Glöckner, F. O. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Research, 41(Database issue), D590–D596. doi:10.1093/nar/gks1219
Quimby, D. C., & Johnson, D. E. (1951). Weights and measurements of Rocky Mountain elk. Journal of Wildlife Management, 15, 57–62.
R Core Team. (2017). R: A language and environment for statistical computing. R Foundation for Statistical Computing, Zurich, Switzerland.
R Core Team. (2019). R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing.
Rime, T., Hartmann, M., Brunner, I., Widmer, F., Zeyer, J., & Frey, B. (2015). Vertical distribution of the soil microbiota along a successional gradient in a glacier forefield. Molecular Ecology, 24(5), 1091–1108. doi:10.1111/mec.13051
Risch, A. C., Jurgensen, M. F., & Frank, D. A. (2007). Effects of grazing and soil micro-climate on decomposition rates in a spatio-temporally heterogeneous grassland. Plant and Soil, 298(1–2), 191–201. doi:10.1007/s11104-007-9354-x
Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., … Weber, C. F. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Applied and Environmental Microbiology, 75(23), 7537–7541. doi:10.1128/AEM.01541-09
Smith, D., Stahler, D., Cassidy, K., Stahler, E., Metz, M., Cassidy, B., … Cato, E. (2018). Yellowstone National Park wolf project annual report 2017. Yellowstone National Park, Mammoth, WY, USA: National Park Service, Yellowstone Center of Resources.
Wang, Q., Garrity, G. M., Tiedje, J. M., & Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Applied and Environmental Microbiology, 73(16), 5261–5267. doi:10.1128/AEM.00062-07
Withington, C., & Sanford Jr., R. (2007). Decomposition rates of buried substances increase with altitude in a forest-alpine tundra ecotone. Soil Biology and Biochemistry, (39), 68–75.

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

bacteria

display_name: bacteria

biogeochemical-cycling

display_name: biogeochemical-cycling

carcass

display_name: carcass

carrion

display_name: carrion

decomposition

display_name: decomposition

detritus

display_name: detritus

fungi

display_name: fungi

soil-biodiversity

display_name: soil-biodiversity

wolf

display_name: wolf

yellowstone

display_name: yellowstone

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: Elk and bison carcasses in Yellowstone, USA


maintainer: Peter Bebi

en: soil-temperature stillberg experiment climate-change vegetation carbon-dioxide

en: Project Description FACE Stillberg

media_type: text/html

en: Link to FACE project description on the WSL website.

format: HTML

en: Metadata FACE Stillberg Data File Description

en: List of descriptions for each of the .csv FACE data files.

en: Data FACE Stillberg

media_type: application/zip

en: Folder in .zip format containing the FACE datasets.

format: ZIP

mimetype: application/zip

en: Metadata FACE Stillberg Experimental Design

en: Schematic illustration and tabular listing of the FACE experimental design.

en:
 Background information

High elevation ecosystems are important in research about environmental change because shifts in climate associated with anthropogenic greenhouse gas emissions are predicted to be more pronounced in these areas compared to most other regions of the world. This project involved a Free Air CO2 Enrichment (FACE) and soil warming experiment located in a natural treeline environment near Davos, Switzerland (Stillberg, 2200 m a.s.l.). Elevated atmospheric CO2 concentrations (+200 ppm) were applied from 2001 until 2009, and a soil warming treatment (+4 °C) was applied from 2007 until 2012. The combined CO2 enrichment and warming treatment reflects conditions expected to occur in this region in approximately 2050. A broad range of ecological and biogeochemical research was carried out as part of this environmental change project.

 Experimental design

The experiment consisted of 40 hexagonal 1.1 m² plots, 20 with a *Pinus mugo* ssp. *uncinata* (mountain pine, evergreen) individual in the centre and 20 with a *Larix decidua* (European larch, deciduous) individual in the centre. A dense cover of understorey vegetation surrounded the tree in each plot, including the dominant dwarf shrub species *Vaccinium myrtillus* (bilberry), *Vaccinium gaultherioides* (group *V. uliginosum agg.*, northern bilberry) and *Empetrum nigrum* ssp. *hermaphroditum* (crowberry) plus several herbaceous and non-vascular species.
At the beginning of the experimental period, the 40 plots were assigned to ten groups of four neighbouring plots (two larch and two pine trees per group) in order to facilitate the logistics of CO2 distribution and regulation. Half of these groups were randomly assigned to an elevated CO2 treatment, while the remaining groups served as controls and received no additional CO2. In spring 2007, one plot of each tree species identity was randomly selected from each of the 10 CO2 treatment groups and assigned a soil warming treatment, yielding a balanced design with a replication of five individual plots for each combination of CO2 level, warming treatment and tree species.

 Data description

Soil and air conditions have been monitored closely throughout the study period, with most measurements made during the combined CO2 x warming experiment (2007-2009). The data comprise of air temperature, soil temperature, soil moisture, sapflow, tree diameter and CO2 measurements.

carbon-dioxide

display_name: carbon-dioxide

climate-change

display_name: climate-change

experiment

display_name: experiment

soil-temperature

display_name: soil-temperature

stillberg

display_name: stillberg

vegetation

display_name: vegetation

en: FACE: Stillberg CO2 enrichment and soil warming study


maintainer: Giulia Francesca Azzurra Donati

en: circuit-theory integrated-modelling urban-spatial-planning biodiversity ecological-corridors human-dominated-landscapes

en: Amphibian breeding sites of national importance - cantons of Aargau and Zürich

en: Amphibian breeding sites of national importance (focal nodes) in the cantons of Aargau and Zürich, 30x30 m resolution, binarised ascii format.

en: Amphibian whole life cycle predictors - 10m resolution

en:
Amphibian whole life cycle predictors include 6 variable types:
(1)	Topographic – slope [°]: derived from the Swiss Topographical Landscape model (swissTLM 3D 2015) at 2 m resolution (elevation uncertainty ± 0.5m) and, subsequently, aggregated to a median value;
(2)	Hydrologic: two predictors including nearest distance to water [m] (all water-related objects and rasterized lotic, lentic and reedbeds from the regional cadastral maps – AV: Daten der Amtliche Vermessung) and runoff coefficient [C] (between 0 and 1 relating amount of runoff to amount of precipitation, mapped AV: Daten der Amtliche Vermessung data against literature values for different land covers);
(3)	Edaphic – soil moisture variability (not available): can gathered from the spatial modelling of Ecological Indicator Values in Switzerland (Descombes et al., 2020) and resampled from 100m to 10m;
(4)	Vegetation: we calculated the median and standard deviation of the Normalized Difference Vegetation Index (NDVI) from all available cloud free Sentinel-2 pixels available between April and October for the years 2016 to 2019 (Drusch et al., 2012) and the median vegetation height [m] (not available) can gathered from Ginzler and Hobi (2015);
(5)	Land-use derived: includes three predictors, grassland density [grassland/cell] (not available) - can be gathered from Pazur et al. (2021), total grasslands area evaluated using a spatial weight matrix and a 250m neighborhood radius around each cell; urbanization [buildings/cell] – average building density for a 250m neighborhood around each cell based on the AV: Daten der Amtliche Vermessung ; and traffic intensity [vehicles/day] – average daily vehicle count (Nationales Personenverkehrsmodells);
(6)	Movement-ecology related: includes three predictors of proximity – nearest forest distance [m]; nearest rock-gravel and sandy area (RGS) distance [m] and nearest road distance [m] (considering any road type) – all calculated based on the AV: Daten der Amtliche Vermessung maps.

en: Amphibian species specific landscape resistance maps.

en: For each amphibian species we used the weighted average Ensemble model suitability map and converted it to a resistance map using an exponential decay function following (Duflot et al., 2018) using the species-specific mpa threshold previously computed (Preau et al., 2020). The function assigns a resistance value of 1,000 when habitat suitability equals zero, and a value of one when habitat suitability is greater than or equal to the binarization threshold. The exponential shape strengthens the barrier effect of less favorable areas.

en: Binarized "green" and "grey" landcover data (Cantons of Aargau and Zürich)

en: "Green" landcover refers to lawns and vegetated urban surfaces such as, for example private gardens, public parks, cemeteries, green stripes along roadsides, etc. While "grey" landcover refers to impervious grounds such as roads, troittoirs, other impervious. Binarized maps, 30x30m resolution.

en: Amphibian monitoring data

en: Processed amphibian occurrence data (2017-2019) in the cantons of Aargau and Zürich. The dataset includes information for the following species: _Alytes obstetricans_,_Bombina variegata_, _Bufo bufo_, _Hyla arborea_, _Ichtyosaura alpestris_, _Lissotriton helveticus_, _Lissotriton vulgaris_, _Rana temporaria_, _Salamandra salamandra_, and _Triturus cristatus_.

en: The data includes (1) amphibian occurrence data (2017-2019) for ten species across the cantons of Aargau and Zürich gathered from the Coordination Center for the Protection of Amphibians and Reptiles of Switzerland (http://www.karch.ch), (2) amphibian whole-life cycle environmental predictors (i.e. topographic, hydrologic, edaphic, vegetation, land-use derived, movement-ecology related), and (3) local urban "green" and "grey" landcover data which can be used to identify opportunities for Blue-Green Infrastructure (through green or grey transitions) in support of regional landscape connectivity.

biodiversity

display_name: biodiversity

circuit-theory

display_name: circuit-theory

ecological-corridors

display_name: ecological-corridors

human-dominated-landscapes

display_name: human-dominated-landscapes

integrated-modelling

display_name: integrated-modelling

urban-spatial-planning

display_name: urban-spatial-planning

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Amphibian and urban-rural landscape data Swiss Lowlands


en: drone uav glacier alps dem dsm orthoimage gletscher griesgletscher orthophoto findelengletscher sankt-annafirn

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Griesgletscher on the 26th of September 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 19th of April 2016.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 21th of July 2015.

en: This .pdf contains an orthoimage for each dataset as well as a Table which displays the number of flights made, image acquired, number of GCP set and the total area covered.

en: Bulk Data Access

en: Link to the S3 bucket where all the data resources can be downloaded in bulk

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 15th of September 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 22th of September 2016.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Griesgletscher on the 3rd of August 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Griesgletscher on the 15th of August 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Sankt Annafirn on the 15th of July 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Sankt Annafirn on the 15th of April 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 2nd of July 2016.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Sankt Annafirn on the 28th of September 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Griesgletscher on the 13th of September 2016.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 09th of April 2015.

en:
This .zip folder contains:

 - An orthoimage (oi) of 0.1m resolution

- A Digital Surface Model (dsm) of 0.1m resolution

- The Agisoft Photoscan automatic processing report

for the UAV field campaign conducted on the Findelengletscher on the 30th of August 2016.

en: This folder contains all UAV flights paths that we created for the Sankt Annafirn, the Findelengletscher and Griesgletscher. They have been created with the software eMotion2, that was developped for the eBee platform from Sensefly.

en:
 UAV-derived DSMs and orthoimages

Unmanned Aerial Vehicle (UAV) surveys were conducted between 2015 and 2016 on the Sankt Annafirn, Findelen- and Griesgletscher, situated in the Swiss Alps. Three surveys at the Sankt Annafirn allowed for a full glacier coverage, four surveys at Griesgletscher allowed an almost full glacier coverage and seven surveys at Findelengletscher allowed for a partial coverage of the glacier tongue (see individual datasets for exact extent). For each survey, a high resolution orthoimage and a Digital Surface Model (DSM) was created.

 UAV surveys:

Prior flight, Ground Control Points (GCPs) were deployed on the glacier surface and measured with a differential GPS (Trimble R7 or Leica GPS 1200). They allowed precise georeferencing of the UAV-derived datasets. UAV flight plans were planned with the software *eMotion 2* and a SenseFly eBee was used as surveying platform. The images were then processed with the software Agisoft Photoscan Pro 1.1.6 .

The location and dates of each survey can be found in the table together with the number of flights performed (Nflights), the number of acquired images (Nimages), the number of GCPs set (NGCPs) and the surveyed area.

A folder for each dataset is available (see folder name in table), which contains:

- An orthoimage *glacier_date_photoscan_oi_CH1903+_LV95_0.1m.tif*

- A Digital Surface Model *glacier_date_photoscan_dsm_CH1903+_LV95_0.1m.tif*

- The Agisoft Photoscan automatic processing report  *glacier_date_photoscan_report.pdf*

where:

- *glacier* is the name of the surveyed glacier

- *date* is the date of the UAV image acquisition

- *photoscan* is the name of the photogrammetric software

- *oi* or *dsm* the type of dataset

- *CH1903+_LV95* is the coordinate system and datum of the dataset

- *0.1m* is the resolution of the dataset in meter

- *.tif* is the extention of the dataset

Details about the UAV surveys, the image processing and the accuracy of the UAV-derived products can be found in this publication below.

Paper Citation:
 > _Gindraux et al. 2017. Accuracy Assessment of Digital Surface Models from Unmanned Aerial Vehicles’Imagery on Glaciers, Remote Sensing, 9, 186, 1-15, [doi: 10.3390/rs9020186](https://doi.org/10.3390/rs9020186)._

The folder UAV_flight_paths.zip contains all UAV flights performed on the Sankt Annafirn, Findelengletscher and Griesgletscher. The flights were planned with the software eMotion2 and have the .afp extention.

alps

display_name: alps

dem

display_name: dem

drone

display_name: drone

dsm

display_name: dsm

findelengletscher

display_name: findelengletscher

glacier

display_name: glacier

gletscher

display_name: gletscher

griesgletscher

display_name: griesgletscher

orthoimage

display_name: orthoimage

orthophoto

display_name: orthophoto

sankt-annafirn

display_name: sankt-annafirn

uav

display_name: uav

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: UAV-derived Digital Surface Models and orthoimages for three alpine glaciers


disdrometer-data-davos-wolfgang

title_for_slug: disdrometer-data-davos-wolfgang

en: precipitation fall-velocity disdrometer davos raclets

en: Disdrometer Data Davos Wolfgang

en:
The zip folder contains one netCDF file per day for Davos Wolfgang from February 8 to March 19 2019 with a temporal resolution of 30 seconds (time in UTC) and includes:

- characteristics of hydrometeors (e.g. diameter, fall velocity, amount per diameter class,...)
- precipitation rate
- radar reflectivity
- visibility range
- weather codes
- instrument performance

en: The dataset contains information on precipitation amount and type for Davos Wolfgang (LON: 9.853594, LAT: 46.835577) from February 8 to March 19 2019. It includes: characteristics of hydrometeors (e.g. diameter, fall velocity, amount per diameter class,...), precipitation rate, radar reflectivity, visibility range, weather codes and instrument performance.

davos

display_name: davos

disdrometer

display_name: disdrometer

fall-velocity

display_name: fall-velocity

precipitation

display_name: precipitation

raclets

display_name: raclets

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal


maintainer: Markus Schlegel

en: metabarcoding spore-dispersal environmental-dna monitoring air-sampling

en: Data files and analysis scripts (Version 2)

en: This ZIP archive contains the data and source code necessary to reproduce the analyses for the manuscript "A case study on the application of spore sampling for the monitoring of macrofungi". Please refer to README.pdf inside the archive for more information.

en: Data files and analysis scripts (Version 1)

en: This ZIP archive contains the data and source code necessary to reproduce the analyses for the manuscript "A case study on the application of spore sampling for the monitoring of fungi". Please refer to README.pdf inside the archive for more information.

en:
Resources associated with the manuscript *A case study on the application of spore sampling for the monitoring of macrofungi* ([Schlegel et al. 2024](https://doi.org/10.1111/1755-0998.13941)).
Airborne environmental DNA (eDNA) monitoring of fungi was evaluated on a species-rich grassland site. Extensive fruiting body surveys yielded 29 waxcap and clavarioid species, 19 of which were also detectable in air eDNA. The analysis indicates that spores of locally occurring fruiting bodies were detected, in addition to a large diversity of common fungal species, while rare and threatened species were under-represented in the air.
*Version 1* of the dataset was submitted during the review phase. *Version 2* contains minor changes and clarifications, and is up to date with the [published supplement](https://onlinelibrary.wiley.com/action/downloadSupplement?doi=10.1111%2F1755-0998.13941&file=men13941-sup-0001-AppendixS1.pdf) of the study.
The contents of the ZIP archive:
* analysis scripts in [R-Markdown](https://rmarkdown.rstudio.com) format, including the rendered PDF output with figures and tables (*analysis* directory)
* eDNA sample metadata (*amplicon_meta*) and results of the amplicon clustering (*amplicon_data*) and the source code of the pipeline (*amplicon_pipeline*)
* positions of all fruiting body groups found in the surveys (3838 GPS points) (*surveys*)
* positions of the spore traps, taxonomy of studied grassland fungi and red-list species, and many more resources
Please refer to `README.pdf` inside the archive for a complete description of the file structure and information on how to reproduce the analyses.

air-sampling

display_name: air-sampling

environmental-dna

display_name: environmental-dna

metabarcoding

display_name: metabarcoding

monitoring

display_name: monitoring

spore-dispersal

display_name: spore-dispersal

en: Supporting information for case study on the application of spore sampling for the monitoring of macrofungi


distributed-acoustic-sensing-brienz

title_for_slug: distributed-acoustic-sensing-brienz

en: natural-hazard rockfall landslide seismic early-warning machine-learning

en: Data access

en:
This dataset contains the Distributed Acoustic Sensing (DAS), radar detection data used for training and result analysis in the GRL paper titled `Automatic Monitoring of Rock-Slope Failures Using Distributed Acoustic Sensing and Semi-Supervised Learning`.
The DAS dataset (both waveform and cross-spectral density matrices), extracted features, labeled dataset, two trained models (feature extraction model and xgboost classification model), scripts to reproduce the whole training and classification processes, and a notebook to replicate the result analysis part are provided under the MIT license. To provide a reasonable data size, we chunked the raw data to a few hundred channels which we used in our project.
Abstract:
Effective use of the wealth of information provided by Distributed Acoustic Sensing (DAS) for mass movement monitoring remains a challenge. We propose a semi-supervised neural network tailored to screen DAS data related to a series of rock collapses leading to a major failure of approximately 1.2 million cubic meters on 15 June 2023 in Brienz, Eastern Switzerland. Besides DAS, the dataset from 16 May to 30 June 2023 includes Doppler radar data for partially ground-truth labeling. The proposed algorithm is capable of distinguishing between rock-slope failures and background noise, including road and train traffic, with a detection precision of over 95%. It identifies hundreds of precursory failures and shows sustained detection hours before and during the major collapse. Event size and signal-to-noise ratio (SNR) are the key performance dependencies. As a critical part of our algorithm operates unsupervised, we suggest that it is suitable for general monitoring of natural hazards.

early-warning

display_name: early-warning

landslide

display_name: landslide

machine-learning

display_name: machine-learning

natural-hazard

display_name: natural-hazard

rockfall

display_name: rockfall

seismic

display_name: seismic

en: Distributed Acoustic Sensing Brienz


maintainer: Joan Casanelles Abella

en: tartu floral-resources floral-inventories urban antwerp plants angiospermae paris zurich poznan

en: List of plant species sampled in the five cities (Antwerp, Paris, Poznan, Tartu and Zurich) during the different sampling perionds (four for Zurich, three for the other cities) with their frequency (oc1-oc16) within a 100 m radius

en: Floral units, mean number of flowers per floral unit and floral abundance of the different plants counted in quadrats in the study sites in Zurich during the four sampling periods

en: Trait values for the sampled plants

en: Data of a survey of flowering plants in 80 sites in five European cities and urban agglomerations (Antwerp, Belgium; greater Paris, France; Poznan, Poland; Tartu, Estonia; and Zurich, Switzerland) sampled between April and July 2018.

angiospermae

display_name: angiospermae

antwerp

display_name: antwerp

floral-inventories

display_name: floral-inventories

floral-resources

display_name: floral-resources

paris

display_name: paris

plants

display_name: plants

poznan

display_name: poznan

tartu

display_name: tartu

urban

display_name: urban

zurich

display_name: zurich

envidat

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

en: Flowering Plants (Angiospermae) in Urban Green Areas in five European Cities


en: ecology switzerland environmental--policy environmental---history historical-carthography wetlands

en: Wetlands 1850 - 2010 (shapefiles)

en:
 History of wetlands - GIS-Layers
 Location of wetlands based on signatures in topographic maps
* Basis 1850: original measuring table sheets Dufour created between the years 1834 and 1863.
* Basis 1900/ Basis 1950: The Siegfried maps were updated and published at regular intervals between 1870 and 1949. This long period of time made it possible to reconstruct wetlands for two periods based on the Siegfried Maps: the first edition (hereinafter referred to as "1900") and the last edition ("1950").
* Basis 2010: Map sheets from the National Map from 2006 to 2011 were used for the most recent time period of the reconstruction.
 Location of wetlands due to reconstructions
The reconstruction **must not** be interpreted in the sense of a surface-accurate separation, but shows without claim to location accuracy the approximate extent of the areas that would have corresponded to the mapping criteria for wetlands used in today's topographic maps around 1850, 1900 and 1950. The basics can be found in Stuber & Bürgi 2019, chapter 5.4.

en:
Naturally, large parts of the Swiss Plateau are characterised by wetlands
and meandering rivers. That this is no longer the case today is the result
of centuries of efforts to obtain dry land. But how did this process take
place? What were the relevant actors and what were their motivations? And
what can be said about the ecological consequences of this development?

In a research project on the history of wetlands in Switzerland since 1700,
we conducted (a) a historical analysis of the development of land use in
wetlands and the actors involved, (b) a historical-cartographic
reconstruction of wetland extent since 1850 and (c) an evaluation of
ecological effects of changes in wetlands on various organisms groups.

The series of GIS layers on wetland history stem from the second part of
the project. The area reconstruction is based on digitized and homogenized
signatures from national map series, as they have been available since
about 1850. Details about the digitalization process and the homogenization
procedures applied ("Rekonstruktionen") are included in Stuber & Bürgi
2019.

Book Citation:
 > Stuber M, Bürgi M (2019) Vom «eroberten Land» zum Renaturierungsprojekt. Geschichte der Feuchtgebiete in der Schweiz seit 1700. "Bristol Schriftenreihe", Band 59. Haupt Verlag, Bern, Stuttgart, Wien. 262 Seiten.

ecology

display_name: ecology

environmental---history

display_name: environmental---history

environmental--policy

display_name: environmental--policy

historical-carthography

display_name: historical-carthography

switzerland

display_name: switzerland

wetlands

display_name: wetlands

educ

en: Education, culture and sport

en: EnviDat: the Environmental Data Portal

en: History of wetlands in Switzerland since 1850


maintainer: Virginie Marques

en: fish marine dna-metabarcoding environmental-dna plankton biodiversity fjord

en: Metabarcoding data

en: Raw metabarcoding data for both primers (teleo and 18S) and corresponding file to perform the demultiplexing and primer trimming

en: Code, data and figures

en: Here there are all the data and code to reproduce the figures and analysis presented in the paper, with a README.md file

en:
This dataset contains the raw environmental DNA data associated with the publication *Environmental drivers of eukaryotic plankton and fish biodiversity in an Arctic fjord* in the journal Polar Biology (2023).

 Methods

**Sampling**
We sampled the Lilliehöök fjord on the west coast of Spitsbergen (Svalbard, Norway) over 3 days from 3 to 5 of August 2021. Samples were taken from the glacier front up to the fjord mouth of the Krossfjorden system, around 30 km long, after the Lilliehöök fjord merged with the mouth of Möller fjord. The fjord’s maximum depth has been recorded at 373 m (Svendsen et al. 2002) and has no sill at its entrance, thereby facilitating water exchange with the open ocean of the West Spitsbergen Current. We used a research vessel to sample 5 sites for a total of 15 samples, sampling 3 depths per site (3-m, chlorophyll a maximum and 85-m, unless sea floor was shallower). Shallow and intermediate samples between 3-m and 12-m represent ~35-L of water filtered in-situ using long tubing and a peristaltic pump, and all other deeper samples were taken from a total of 3 Niskin bottles (General Oceanics), representing 22-L of water sampled per sample. Water was filtered through a VigiDNA filtration capsule (SPYGEN) with a 0.20-µm pore size using an Athena peristaltic pump (Proactive Environmental Products, Bradenton, Florida) with a flow rate of ~1-L/min. Each sample was handled with single use tubing and gloves.

**Molecular**
To perform the amplification, we used two sets of primers: teleo (forward: ACACCGCCCGTCACTCT, reverse: CTTCCGGTACACTTACCATG; Valentini et al. 2016) and the universal eukaryotic 1389F/1510R primer pair, amplifying the V9-18S rDNA gene (Amaral-Zettler et al. 2009) (forward: TTGTACACACCGCCC, reverse: CCTTCYGCAGGTTCACCTAC).

 Data content:

+ Metabarcoding data: This zip file contains the 2 sequencing libraries filtered to only retain the samples used in the present study.
+ Code, data and figure: This zip file contains all data and code to reproduce the figures and the analysis in the study, with an associated README explaining the content of each folder.

 Additional informations
For more details, please see the Methods in the associated publication: DOI: 10.1007/s00300-023-03187-9.

biodiversity

display_name: biodiversity

dna-metabarcoding

display_name: dna-metabarcoding

environmental-dna

display_name: environmental-dna

fish

display_name: fish

fjord

display_name: fjord

marine

display_name: marine

plankton

display_name: plankton

en: Data: Environmental drivers of eukaryotic plankton and fish biodiversity in an Arctic fjord


en: modeling urbanization statistical-model land-use landscape-services land-use-change switzerland

en: Scripts for Assessing potential landscape service trade-offs driven by urbanization in Switzerland

en:
Scripts ordered by their main purpose within sequential folders in this ReadMe marked with

Within the folders the capital letters refer to the parallel or sequential use of the scripts; scripts here in ReadMe are marked with

 00_Function Scripts
 Functions.r
Script containing several functions for analyzing statistical models: cohen.kappa, meva.table, D2, max-Tss, plot-TSS

 Service_functions.R
Script containing functions for calculating all services and returning sums (and median for bd and rec) as well as rasters

 01_Data Preparation
 A_AS_raster_preparation.R
Preparation of the land use rasters for the years 1985, 1997, and 2009 based on the original Arealstatistik (AS/ASCH); Assigning our categories and assining new values to all "construction sites" for calibrating the models
For model projection starting point (2009) and for service calculations (all ASCH years) the "other" categories were assigned to the land use categories based on dominant neighbor, so that there were no empty cells
Creation of water & mountainous area raster that stayed stable throughout the projections
Creation of MFH raster (2009), to have as base for all new projection years, as this category never becames any other

 B_Factors_preparation.R
Preparation of the factor rasters that were used to calculate service values per cell (see Table 2): LWfactor, Transportation, Vol, Rivers, Lakes, slope, TWW, Hazard_risk

 B_Predictors_preparation.R
Preparation of the raster of the predictors used for calibrating (ASCH years) and projecting (stable) the transition models that do not depend on the land use categories: Temperature (ASCH years), solar radiation (stable), Public transportation access (stable), Elevation (stable), Workers density (ASCH years), Population change (ASCH transitions)

 B_Regions_preparation.R
Preparation of bioregions as raster to be used to mask asxxxx later to calculated services and land use area just for one region

 C_Pred_proj.R
Projecting the dynamic predictors into the future for projecting the land use transitions with the model: Linear extrapolations for workers densities and population; more detailed extrapolation of temperature (based on data for all stations that matched between the temp table we used for the calibration period and the projections from CS2M, and then making a spatial interpolation for all years based on the interpolation of temperature for all years per station)

 02_Model calibration
 A_Bioreg_pred_corr.R
Calculating correlation between the predictors for each bioregion and land use category, Output are those with too high correlation so that one predictor can be removed from transitions were both predictors occur

 A_trans_numbers_bioreg.R
Calculate transition numbers for each bioregion to decide which transitions to primarily consider

 B_Bioreg_ext_eval.R
Calculating TSS based on the external evaluation for the bioregion transitions that happened often enough to try to model them (calibration data: 85-97, evaluation data: 97-09)

 B_Bioreg_int_eval.R
Calculating TSS based on the internal evaluation for the bioregion transitions that happened often enough to try to model them (calibration data: 70%, evaluation data: 30%)

 C_Bioreg_TSS_analysis.R
Extracting TSS numbers from internal and external evaluation, calculating mean and median for each transition for each region for each evaluation method and saving as csv per region to be able to select only the transitions with high enough model performance

 D_Bioreg_Models_for_prediction.R
Calibration of the models that were then used to project land use transitions into the future; calibrated separately for each bioregion

 E_Residuals_trans_models.R
Calculation of Moran's I statistics to analyze the residuals of the land use transition models

 03_Projections
 A_Bioreg_Projections_run.R
Master script for: Extrapolaitng into the future with the regional models and then assembling for each time step to calculate new AS predictors; Several Outputs: txt file to test for NAs in Predictors, txt file with random seed number for sampling, txt file with times for calculation steps (check progress of script while running), csv sheet with projected transition frequencies, tif-rasters for every region and time step, pdf with raster plot for every region and time step

 04_Service calculation
 A_Bioreg_as_services.R
Calculate services per bioregion for years 1985, 1997, 2009

 A_Service_calculation_bioregions.R
Extracting service numbers(sums) per bioregions for all projection runs and storing them in a list of matrixes to be available for further plotting; can be run for swiss wide projections or regional ones

 05_Analysis
 A_Bioreg_ASCH_development.R
Extracting the occurrence of each category per region for all runs (first part), then for ASCH years, and then plotting them nicely; only works for regional projections

 A_Bioreg_Predictors_importanceD2.R
Calculating D2 for each transition and each considered predictor from bivariate models to get value of importance

 B_Bioreg_Pred_imp_analysis.R
Using calculated D2 to sort Predictors according to their importance for each transition

 A_Bioreg_Spatial_results_analysis.R
Plotting randomly sampled runs for ASCH and services to analyses their spatial development

 A_Bioreg_Transition_numbers_check_v1.R
Analyzing the frequency for all transitions for all regions and producing the point plots for each transition separately

 A_Bioreg_Transition_numbers_check_v2.R
Analyzing percentage and numbers of all transitions but plotting them per region and originating ASCH category

 B_Transition_Analysis.R
Further Analysis of transition numbers

 A_Service_analysis.R
Plotting of the service development for bioregions in several ways for analyzing

en: Data for Assessing potential landscape service trade-offs driven by urbanization in Switzerland

en:
The provided folder structure follows the set up of the scripts.
Data is only provided when it has been fully produced by us, otherwise dummies are given.
For data that will be produced by one of the scripts, no dummies are given.
The folders may contain sub-folders, which would be filled if the scripts were all run

 00_Analysis_Data
Data produced throughout the process of the calibration and projection that is later needed for analysis

 01_Raw_Data
Primary data from external sources that were used to compile the rasters and numbers used during the calibration and projection of landscape transitions as well as service calculation

 02_Prepared_Data
Data that will be prepared for further use by scripts are put in here
Additionally it contains the parameter sheets for the transition predictors that we compiled

 03_Model_Data
 Model outputs
Folder were the projection models statistics will be saved
 Models
Folder into which the actual transition models will be saved
 Residuals
Folder into which data on models residuals will be saved for future use

 04_Projected_Data
The output of the projections will be saved into this folder

 05_Service_Data
Any service calculations will be saved in here

 06_Results_Analysis
Any Analysis plots and tables as well as intermediate data will be saved here in the respective sub-folder

en:
Data and scripts of publication:

Madleina Gerecke, Oskar Hagen, Janine Bolliger, Anna M. Hersperger, Felix Kienast, Bronwyn Price, Loïc Pellissier (2019) Assessing potential landscape service trade-offs driven by urbanization in Switzerland. Palgrave communications.

Contains land use projections for Switzerland and scripts and data for these projections as well as the calculation of landscape services.

Data Folder:
Contains sub-folder with the data necessary for this study (provided were no copyright issues, otherwise placeholders with descriptions), and folders where produced data may be stored

Scripts Folder:
Contains scripts organized into subfolders depending on their purpose

Note: Some abbreviations within the scripts and data are derived from German words and not English.

land-use

display_name: land-use

land-use-change

display_name: land-use-change

landscape-services

display_name: landscape-services

modeling

display_name: modeling

statistical-model

display_name: statistical-model

urbanization

display_name: urbanization

en: Land use projections and services for Switzerland


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

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

complex-terrain

display_name: complex-terrain

cosmo-wrf

display_name: cosmo-wrf

grayzone

display_name: grayzone

numerical-simulation

display_name: numerical-simulation

wrf

display_name: wrf

educ

en: Education, culture and sport

envidat

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

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


en: wagital snow-water-equivalent meteodat-gmbh gcos snow

en:
Total water reserves of the snow cover [mio m3] for Wägital catchment, Switzerland, for reference date April 1. Data is separated in 2 elevation zones 900m-1500m asl and 1500m-2300m asl. Time period 1943-2024, status 2024-04-22.

Funded currently or in the past by:

- Federal Office of Meteorology and Climatology MeteoSwiss in the context of GCOS Switzerland
- Meteodat GmbH
- Institute of Geography, University of Zurich
- WSL Institute for Snow and Avalanche Research SLF
- Institute of Geography, ETH Zurich (IAC ETH Zurich)
- AG Kraftwerk Wägital (AXPO and EWZ)

See also https://www.meteodat.ch/waegital.html

gcos

display_name: gcos

meteodat-gmbh

display_name: meteodat-gmbh

snow

display_name: snow

snow-water-equivalent

display_name: snow-water-equivalent

wagital

display_name: wagital

en: EnviDat: the Environmental Data Portal

en: Snow water equivalent for reference date April 1 for Wägital catchment, starting 1943


en: x-ray-tomography microstructure fem snow metamorphism time-series

en: Simulation parameters and outputs for "A rigorous approach to the specific surface area evolution in snow during temperature gradient metamorphism"

format: ODS

en:
In the associated study [1], two time-lapse temperature gradient metamorphism series of three-dimensional micro-computed tomography images of snow (obtained by [2]) have been used to model the decrease of specific surface area (SSA) over time based on the pore-scale physics. We conducted finite element simulations of one-way coupled heat and mass diffusion in order to estimation the spatial pattern of water vapor deposition and sublimation, which controls the evolution of the SSA over time. We notably studied the influence of the condensation coefficient, a key but poorly constrained physical parameter.
This dataset provides the parameters used for the mesh generation and the finite element simulations. It also includes the ice fraction, specific surface area per unit volume and surface average of mean curvature and vapor field obtained as outputs from the mesh process and the simulations.
[1]: Braun, A., Fourteau, K., and Löwe, H.: A rigorous approach to the specific surface area evolution in snow during temperature gradient metamorphism, EGUsphere [preprint], https://doi.org/10.5194/egusphere-2023-1947, 2023.
[2]: Pinzer, B. R., Schneebeli, M., and Kaempfer, T. U.: Vapor flux and recrystallization during dry snow metamorphism under a steady temperature gradient as observed by time-lapse micro-tomography, The Cryosphere, 6, 1141–1155, https://doi.org/10.5194/tc-6-1141-2012, 2012.

fem

display_name: fem

metamorphism

display_name: metamorphism

microstructure

display_name: microstructure

time-series

display_name: time-series

x-ray-tomography

display_name: x-ray-tomography

en: Simulation parameters and outputs for a rigorous approach to the specific surface area evolution in snow during temperature gradient metamorphism


en: fish-composition fish teleo marine curacao caribbean-region edna coral-reefs environmental-dna

en: Data access

en:
Fish taxonomic, functional and phylogenetic composition variations recovered from eDNA metabarcoding along the coast of Curacao

In February 2020, we collected a total of 20 water samples, from 10 stations, with two filtration replicates per station, in the two investigated coastal areas. Each station consisted of a transect of 2 km at an overall constant distance from the coast. We recorded the GPS coordinates at the start and end of the transect, which we used to map the transect positions. We conducted eDNA sampling by using a filtration device composed of an Athena® peristaltic pump (Proactive Environmental Products LLC, Bradenton, Florida, USA; nominal flow of 1.0 L/min), a VigiDNA® 0.20 μM cross-flow filtration capsule (SPYGEN, le Bourget du Lac, France) and disposable sterile tubing for each filtration capsule. We performed two filtration replicates in parallel on each side of a boat, at each station, for 30 min corresponding to a volume of ~30 L of water filtered by each capsule. At the end of each filtration, the water inside the capsules was emptied, and we filled the capsules with 80 ml of CL1 Conservation buffer (SPYGEN, le Bourget du Lac, France) and stored at room temperature. We followed a strict contamination control protocol in both field and laboratory stages (Valentini et al., 2016). Each water sample processing included the use of disposable gloves and single-use filtration equipment to avoid any risk of contamination. Libraries were prepared with ligation using the MetaFast protocol (Fasteris).

caribbean-region

display_name: caribbean-region

coral-reefs

display_name: coral-reefs

curacao

display_name: curacao

edna

display_name: edna

environmental-dna

display_name: environmental-dna

fish

display_name: fish

fish-composition

display_name: fish-composition

marine

display_name: marine

teleo

display_name: teleo

en: Environmental DNA Marine Curacao 2020