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

en: Soil properties

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

en: Ecosystem coupling_interactions

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

en: Nematode data

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

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

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

en: Distance to target results

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

en: Vegetation data

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

en: Networks

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

en: Fungi data

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

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

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

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

Paper Citation:

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

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

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

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

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

bacteria

display_name: bacteria

biodiversity

display_name: biodiversity

fungi

display_name: fungi

grassland

display_name: grassland

nematodes

display_name: nematodes

plants

display_name: plants

restoration

display_name: restoration

soil

display_name: soil

soil-biodiversity

display_name: soil-biodiversity

soil-biology

display_name: soil-biology

soil-chemistry

display_name: soil-chemistry

topsoil-removal

display_name: topsoil-removal

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


maintainer: Sebastian Schlögl

en: snow atmosphere stable-boundary-layer atmospheric-turbulence turbulent-sensible-heat-fluxes

state: active

type: dataset

en: Turbulence data WFJ07

media_type: application/zip

en: lower sensor

format: ZIP

mimetype: application/zip

en: Turbulence data AA00

en: sensor in 25 m height

en: Turbulence data PM07

en: NWW

en: sensor in 1 m height

en: Turbulence data GR00

en: upper sensor

en: SEE

en: sensor in 0.2 m height

en: Turbulence data WFJ11

en: sensor in 2 m height

en:
We used five different atmospheric turbulence datasets from four test sites, with these sites showing differences in their topographical characteristics. We chose one typical alpine test site with high topographical complexity (Weissfluhjoch, Davos, Switzerland) and three test sites consisting of one glacier site (Plaine Morte, Crans-Montana, Switzerland) and two polar sites (Greenland and Antarctica) representing a quasi-ideal site with homogeneous surface and quasi infinite fetch in all directions.
The turbulent sensible heat flux was calculated using the eddy-covariance method. Note that the sonic temperature fluctuations have been converted into virtual temperature fluctuations.
Three-dimensional wind velocity and air temperature were processed using a linear detrending (Rannik and Vesala, 1999) and a planar fit approach (Massmann and Lee, 2002) to rotate the coordinate system. Air temperature, relative humidity and air pressure from weather stations were used to calculate air properties, which are required for the data processing. The weather stations are located in the immediate vicinity of the turbulence tower and are affected by the same air masses. Turbulence data were averaged to 30-min intervals, whilst changing to a 15-min time interval marginally affects the heat fluxes at the Weissfluhjoch test site (Mott et al., 2011). Note that we define a negative sensible heat flux as being directed towards the snow surface and a positive sensible heat flux as being directed upwards.
The selected datasets and corresponding test sites are briefly introduced below:

Weissfluhjoch 2007 (WFJ07): A vertical set-up of two three-dimensional ultrasonic anemometers (CSAT3, Campbell Scientific, Inc.) was used at the traditional field site Weissfluhjoch (2540 m asl.) to measure three-dimensional wind velocity and air temperature at a frequency of 20 Hz. The sensors were mounted 3 m and 5 m above the ground and provided reliable data for 50 days between 11 February 2007 and 24 April 2007. Further information on the field campaign can be found in Stössel et al. (2010) and Mott et al. (2011).

Weissfluhjoch 2011-13 (WFJ11): Three-dimensional wind velocity and air temperature were recorded at 5 m above the ground at a frequency of 10 Hz with a three-dimensional ultrasonic anemometer (CSAT3). The analysis was conducted for data obtained between February and March in the years 2011-13.

Plaine Morte 2007 (PM07): Two three-dimensional ultrasonic anemometers (CSAT3) were installed on a horizontal boom facing opposite directions (west-north-west vs. east-south-east) at 3.75 m above the ground to measure air temperature and three-dimensional wind velocity at 20 Hz. The data were collected at the almost flat field site on the Plaine Morte glacier (2750 m asl.) near Crans-Montana, Switzerland from February to April 2007. High quality meteorological data were additionally recorded and used to force the model. A detailed description about the set-up at the Plaine Morte glacier can be found in Huwald et al. (2009) and Bou-Zeid et al. (2010).

Greenland 2000 (GR00): High-frequency three-dimensional ultrasonic anemometer measurements (CSAT3) were recorded at 50 Hz at the Summit Camp (72.3 °N, 38.8 °W, 3208 m asl.) located on the northern dome of the Greenland ice sheet. Data were collected at 1 m and 2 m above the snow surface during summer in 2000 and 2001. Additionally, meteorological measurements were obtained for the post processing and used to force the model. More information about the field campaign can be found in Cullen et al. (2007, 2014).

Antarctica 2000 (AA00): A set-up of three vertical three-dimensional ultrasonic anemometers (DA-600, Kaijo Denki) were installed at Mizuho Station (70°42' S, 44°20' E, 2230 m asl.) in Eastern Antarctica at 0.2, 1 and 25 m and recorded turbulence data at a frequency of 100 Hz from October to November 2000. Longwave and shortwave radiation, relative humidity, air and snow surface temperature were additionally measured and used to force the model. More information about the field campaign can be found in Nishimura and Nemoto (2005).

atmosphere

display_name: atmosphere

atmospheric-turbulence

display_name: atmospheric-turbulence

snow

display_name: snow

stable-boundary-layer

display_name: stable-boundary-layer

turbulent-sensible-heat-fluxes

display_name: turbulent-sensible-heat-fluxes

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

en: Education, culture and sport

it: Istruzione, cultura e sport

envidat

fr: EnviDat: le portail de données environnementales

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

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

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

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

en: How do stability corrections perform in the stable boundary layer over snow?


maintainer: Adam Winstral

en: snow-depth snow jim snow-simulation jules-investigation-model

en: HS and SWE observations

en: Automated (HS) and manual (SWE) data used to evaluate model simulations

en: DEM

media_type: text/plain

en: 1-km resolution DEM used for downscaling gridded NWP data to stations and for distributing BDE evaluations at HS-stations to SWE-stations

format: TXT

mimetype: text/plain

en: Station forcing data

en: All forcing data need to run JIM at HS and SWE stations. These are data downscaled directly from COSMO NWP data (WY14, WY15) and further enhanced by additional station observations of snow depth (HS) and air temperature (Ta) in WY17. Further information can be found in README file. Note that in WY14 and WY15, shortwave radiation was not split into direct and diffuse (not necessary for JIM) and that the default method in JIM for splitting precipitation phase was used. In WY 17, direct and diffuse were combined as input to JIM and JIM was forced with pre-processed liquid and solid precipitation phases.

en: File/data descriptor

en: These data were used to drive and evaluate Jules Investigation Model (JIM) snow simulations. The data provided are the forcing data used for the "deterministic" runs as described in Winstral et al., 2019. The bias-detecting ensemble (Winstral et al., 2019) used observed snow depths (HS) to detect biases in these deterministic simulations related to precipitation and energy inputs to JIM. Simulations that included the BDE evaluations substantially improved JIM simulations.

jim

display_name: jim

jules-investigation-model

display_name: jules-investigation-model

snow-depth

display_name: snow-depth

snow-simulation

display_name: snow-simulation

de: EnviDat: das Umwelt-Datenportal

en: EnviDat: the Environmental Data Portal

it: EnviDat: il portale dei dati ambientali

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester

en: Meteorological data used to develop and validate the bias-detecting ensemble (BDE)


maintainer: EnviDat Support

en: air-temperature precipitation-amount orographic-precipitation climate-change cmip5

en: Download Server

format: SERVICE

en:
Predicting future climatic conditions at high spatial resolution is essential for many applications in science. Here we present data for monthly time series of precipitation and minimum and maximum temperature for four downscaled global circulation models. We used model output statistics in combination with mechanistic downscaling (the CHELSA algorithm) to calculate mean monthly maximum and minimum temperatures, as well as monthly precipitation sums at ~5km spatial resolution globally for the years 1850-2100. We validated the performance of the downscaling algorithm by comparing model output with observed climates for the years 1950-2069.

CHELSA_cmip5_ts is licensed under a Creative Commons Attribution 2.0 Generic (CC BY 2.0) license.

air-temperature

display_name: air-temperature

climate-change

display_name: climate-change

orographic-precipitation

display_name: orographic-precipitation

precipitation-amount

display_name: precipitation-amount

en: High resolution monthly precipitation and temperature timeseries for the period 2006-2100


en: urbanization snails insects species-richness-models arthropods urban-ecology diversity sdms species-abundance-distributions spiders birds

en: abundance and occurrence data of the 1446 species studied, used to investigate species-abundance distributions and species-occurrence distributions

en: contains the raw data to perform the species richness models of the rare species

en: contains the raw data to perform species distribution models, only for common species

en: contains the raw data to perform the species richness models of the hypercommon species

en: contains the raw data to perform the species richness models of all species

en:
Richness, site occurrence and abundance data of bees, beetles, birds, hoverflies, net-wingeds, true bugs, snails, spiders, milipides, wasps collected in the city of Zurich using different sampling techniques, and the environmental variables for each sampling site.
Data are provided on request to contact person against bilateral agreement.

arthropods

display_name: arthropods

birds

display_name: birds

diversity

display_name: diversity

insects

display_name: insects

sdms

display_name: sdms

snails

display_name: snails

species-abundance-distributions

display_name: species-abundance-distributions

species-richness-models

display_name: species-richness-models

spiders

display_name: spiders

urban-ecology

display_name: urban-ecology

urbanization

display_name: urbanization

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: Ecological properties of urban ecosystems. Biodiversity dataset of Zurich


causal-effect-of-mfoz

title_for_slug: causal-effect-of-mfoz

en: selection-bias built-up-land-change the-major-function-oriented-zone fujian-province causality propensity-score-matching

en: Data and Stata code for repeating the results

en:
Data.xlsx is the variables of 952 townships
Table 1-2 NN(1-1)_noreplace_caliper(0.01).do is STATA code to produce the results in Table 1 and 2
Figure 5 Distribution of the propensity scores.R is R code rto produce Figure 5
Figure 5 Standard mean difference.R is R code rto produce Figure 5

en:
Title: Closer to causality: How effective is spatial planning in governing built-up land expansion in Fujian Province, China?

Research objective: The Major Function Oriented Zone (MFOZ), the first strategic spatial plan in China, is developed to achieve a coordinated regional development, through spatial regulation and zoning of development. The MFOZ he MFOZ divided land into four major function-oriented zones: The development-optimized zone, the development-prioritised zone, the development-restricted zone, and the development-prohibited zone. We used propensity score marching to evaluate the effect of the MFOZ on built-up land expansion in Fujian Province over three time intervals (2013–2015, 2013–2018 and 2013–2020).

Data: Data.xlsx contains the variables of 954 towns in Fujian Province.
Town_ID is the town unique ID;
County_ID is the county unique ID;
City_ID is the city unique ID;
MFOZ is the the development-prioritised zone and the development-restricted zone (The development-optimized zone and the development-prohibited zone are excluded);
Builtup_13_15 is the built-up land expansion from 2013 to 2015;
Builtup_13_18 is the built-up land expansion from 2013 to 2018;
Builtup_13_20 is the built-up land expansion from 2013 to 2020;
Dis2water is the Euclidean distance from the town to the nearest waterbody;
Slope is the the average slope within the town;
GDP is the average GDP in 2010 within the town;
Pop is the average population in 2010 within the town;
Road is the average population in 2010 within the town;
Dis2city is the Euclidean distance from the town to the nearest prefectural city centre;
Nei_Arable, Nei_Forest, and Nei_Built.up are the area of arable land, forest land, and built-up land neighbouring town i in 2010.

Method: we used the propensity score matching to compare the changes in the amount of built-up land in the towns of the development-prioritised zone with the matched towns of the development-restricted zone. Additionally, we used three evaluation intervals (2013–2015, 2013–2018 and 2013–2020) to evaluate temporal variation in the causal effect of the MFOZ on built-up land expansion.

built-up-land-change

display_name: built-up-land-change

causality

display_name: causality

fujian-province

display_name: fujian-province

propensity-score-matching

display_name: propensity-score-matching

selection-bias

display_name: selection-bias

the-major-function-oriented-zone

display_name: the-major-function-oriented-zone

en: Causal effect of MFOZ


radiosondes

maintainer: Jörg Wieder

title_for_slug: radiosondes

en: temperature wind radiosonde wind-direction pressure holoballoon raclets relative-humidity

en: Radiosond Template Files

en:
The "template" files contain standard parameters, i.e. ambient pressure (mbar and Pa), temperature (°C), relative humidity (%), wind speed (m s‾ ¹) and wind direction (true deg), as well as additional information, e.g. time (UTC and local time), dew point (°C), mixing ration (g/kg), density (kg m³), altitude (m MSL), height above ground (ft and m AGL), and geographic coordinates (longitude, latitude). If the wind speed was zero, no wind direction is given.
The zip files contain ascents on the following days:
February 10 (morning and afternoon),
February 11 (morning),
February 16 (afternoon),
February 22 (morning),
March 1 (morning),
March 5 (morning),
March 11 (morning),
March 12 (morning),
March 13 (morning),
March 14 (morning),
March 17 (afternoon), and
March 18.
The data was managed with Windsond Receiver Application Demo Version 2.92, (C) Sparv Embedded AB 2016.

en: Radiosond Sounding Files

en:
The "soundings" radiosond data contains information on ambient pressure (mbar), temperature (C), relative humidity (%), wind speed (m s‾ ¹) and wind direction (true deg), and are ordered with respect to the height above ground (m AGL) with a step size of 20 m. The first line of each file idicates the geographic location and the start time (in UTC) of the ascent. If the wind speed was zero, no wind direction is given.
The zip files contain ascents on the following days:
February 10 (morning and afternoon),
February 11 (morning),
February 16 (afternoon),
February 22 (morning),
March 1 (morning),
March 5 (morning),
March 11 (morning),
March 12 (morning),
March 13 (morning),
March 14 (morning),
March 17 (afternoon), and
March 18.
The data was managed with Windsond Receiver Application Demo Version 2.92, (C) Sparv Embedded AB 2016.

en: Radiosondes (Windsond, Sparv Embedded AB) were started in Davos Wolfgang (LON: 9.853594, LAT: 46.835577) to report height profiles of pressure, relative humidity and temperature at specific days. In addition to regular launches of radiosondes, sondes were attached to [HoloBalloon](https://www.envidat.ch/group/clouds-in-situ-raclets) to report the ambient conditions of the in-situ measurements. Further profiles of meteorological measures were recorded at [HoloGondel](https://www.envidat.ch/group/clouds-in-situ-raclets) which was installed at the gondola moving between Gotschnaboden and Gotschnagrat at 2285 m a.s.l.

holoballoon

display_name: holoballoon

pressure

display_name: pressure

raclets

display_name: raclets

radiosonde

display_name: radiosonde

relative-humidity

display_name: relative-humidity

temperature

display_name: temperature

wind

display_name: wind

wind-direction

display_name: wind-direction

envidat

en: EnviDat: the Environmental Data Portal

en: Radiosondes


en: ice chemistry snow-atmosphere-interactions laboratory-study atmosphere

en:
The table lists all NEXAFS raw data with the following columns:

* name:						raw file name
* xRayEkinRange_eV: 		Smallest and largest X-ray photon energy in eV
* scientaPassEnergy_eV: 	Pass energy of electron analyzer in eV
* scientaDwellTime_ms: 		Dwell time of electron analyzer in ms
* regionName:				O K-edge for all data
* scientaAcquisitionMode: 	Fixed energy (snapshot) mode
* scientaEkinRange_eV: 		Smallest and largest energy of the kinetic energy window used for detection in this work
* scientaEkinStep_eV		Not defined in fixed energy mode
* scientaLensMode:			Tune table
* scientaRegionIterations: 	Iterations (sweeps)
* scientaSequenceIterations: Iterations (sweeps)
* sampleTemp_dC: 			Temperature as monitored at the sample holder in °C
* cellPressure_mbar: 		Total pressure in measurement cell in mbar
* iceTemp_dC: 				Temperature of ice sample calculated based on vapour pressure (assuming total pressure equals partial pressure of water) in °C
* smplX_mm: 				positon of measurement spot, x dimension in mm
* smplY_mm: 				positon of measurement spot, y dimension in mm
* smplZ_mm: 				positon of measurement spot, z dimension in mm
* sealingTemp:				temperature at base of sample manipulator
* lastModifiedDatestr: 		time and date when data aquisition was finished
* sample: 					sample description

en: Near Ambient Pressure partial Auger-Meitner electron yield X-ray absorption spectra and photon flux data. X-ray energy in eV in column 1, analyzer electron count intensity in arb. units in column 2.

en: Processed data as shown in Figure 4 of publication. X-ray energy in eV in column one, electron count intensity in arb. units in column 2.

en: Data are compiled that have been used to demonstrate the impact of high water partial pressure on X-ray absorption spectra of ice.

atmosphere

display_name: atmosphere

chemistry

display_name: chemistry

ice

display_name: ice

laboratory-study

display_name: laboratory-study

snow-atmosphere-interactions

display_name: snow-atmosphere-interactions

en: An ice oxygen K-edge NEXAFS spectroscopy data set on gas-phase processing


maintainer: Joan Casanelles Abella

en: floral-inventories plants tartu angiospermae urban paris antwerp zurich floral-resources 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

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: Flowering Plants (Angiospermae) in Urban Green Areas in five European Cities


bryophytes-of-europe-traits-bet-dataset

title_for_slug: bryophytes-of-europe-traits-bet-dataset

en: red-list iucn liverworts indicator-values mosses life-history-traits hornworts bioclimatic-variables

media_type: application/pdf

en: PDF including the (trait) explanations.

format: PDF

en: Text file including the BET dataset.

en: Excel file including the BET dataset and (trait) explanations.

format: XLS

en: The Bryophytes of Europe Traits (BET) dataset includes values for 65 biological and ecological traits and 25 bioclimatic variables for all 1816 bryophytes included in the European Red List (Hodgetts et al. 2019). The traits are compiled from several regional trait datasets and manually complemented using Floras, species-specific literature and expert knowledge. The bioclimatic variables are calculated using the European range of each species. Details regarding the trait compilation and extraction of bioclimatic variables can be found in the corresponding data paper (Van Zuijlen et al. 2023).

bioclimatic-variables

display_name: bioclimatic-variables

hornworts

display_name: hornworts

indicator-values

display_name: indicator-values

iucn

display_name: iucn

life-history-traits

display_name: life-history-traits

liverworts

display_name: liverworts

mosses

display_name: mosses

red-list

display_name: red-list


intratrait

maintainer: Christian Rixen

title_for_slug: intratrait

en: climate-change elevation-range neighbour-interactions cold-adapted-plants elevation-gradient plant-traits

media_type: text/csv

en: Plant trait data along elevation gradients on mountains

format: CSV

en:
This data set was used to test whether species specialized to high elevations or with narrow elevational ranges show more conservative (i.e. less variable) trait responses across their elevational distribution, or in response to neighbours, than species from lower elevations or with wider elevational ranges. We did so by studying intraspecific trait variation of 66 species along 40 elevational gradients in four countries (Switzerland, Australia, New Zealand, China) in both hemispheres. As an indication of potential neighbour interactions that could drive trait variation, we also analysed plant species’ height ratio, its height relative to its nearest neighbour.
The following traits and parameters were measured and are available in this data set:
As an indication of plant stature, we measured vegetative and generative height, where vegetative height was distance from soil to highest vegetative leaf and generative height was distance to the highest point on the reproductive shoot. As a measure of reproductive investment, we noted the presence of flowers on the randomly chosen individuals (see below). As a measure of individual and genet basal area, we measured individual plant and patch diameters, in two dimensions (along the largest diameter and perpendicular to it). In clonal plant species, plant diameter was equivalent to an individual rosette, whereas patch diameter referred to the whole genet and could represent the size of a tuft, tussock or cushion. For genera with more singular growth forms (e.g., some Gentiana species) plant and patch diameter were the same. The two diameter measurements were made at right angles, allowing estimates of patch and plant areas to be calculated as an ellipse (i.e., area = 0.5 a 0.5 b Π). All traits were measured on ten randomly selected individuals per site. Flower count data was considered in a binary fashion on a per individual basis (because for some species individuals only produce one flower when flowering) so that the presence or absence of flower(s) was a nominal value between 0 and 10 for each species at each site. We then collected at least three leaves (up to 30 for small and light leaves) from each of the first three individuals selected from each species for determination of leaf dry matter content (LDMC) and specific leaf area (SLA). For calculations of LDMC and SLA, fresh leaves were scanned on a flatbed scanner to determine leaf area. Leaves were then weighed on a balance to a precision of +/- 0.001g, prior to being air dried and reweighed with a balance to a precision of +/- 0.0001g. LDMC was calculated by dividing dry leaf mass by fresh leaf mass. SLA was calculated by dividing leaf area by dry leaf mass. Additionally, within an area of 10 cm diameter around the target individual, we determined the tallest neighbouring species and measured its vegetative and generative height, and estimated the percent cover of the target species, other vegetation, rock, and bare soil.
For more details see Rixen et al. 2022, Journal of Ecology.

climate-change

display_name: climate-change

cold-adapted-plants

display_name: cold-adapted-plants

elevation-gradient

display_name: elevation-gradient

elevation-range

display_name: elevation-range

neighbour-interactions

display_name: neighbour-interactions

plant-traits

display_name: plant-traits

envidat

en: EnviDat: the Environmental Data Portal

en: intratrait


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

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

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


maintainer: Sofia Pagliarin

en: barcelona imperviousness soil-sealing milan urban-atlas copernicus

en: Information on data sources

en: Shareable project in Esri ArcGIS format including the HRL IMD data in raster format for each of the territorial boundaries as specified in letter a)-h), see Dataset description

en:
Dataset description

This dataset is a recalculation of the Copernicus 2015 high resolution layer (HRL) of imperviousness density data (IMD) at different spatial/territorial scales for the case studies of Barcelona and Milan. The selected spatial/territorial scales are the following:

 * a)	Barcelona city boundaries
 * b)	Barcelona metropolitan area, Àrea Metropolitana de Barcelona (AMB)
 * c)	Barcelona greater city (Urban Atlas)
 * d)	Barcelona functional urban area (Urban Atlas)
 * e)	Milan city boundaries
 * f)	Milan metropolitan area, Piano Intercomunale Milanese (PIM)
 * g)	Milan greater city (Urban Atlas)
 * h)	Milan functional urban area (Urban Atlas)

In each of the spatial/territorial scales listed above, the number of 20x20mt cells corresponding to each of the 101 values of imperviousness (0-100% soil sealing: 0% means fully non-sealed area; 100% means fully sealed area) is provided, as well as the converted measure into squared kilometres (km2).

 Dataset composition

The dataset is provided in .csv format and is composed of:

_IMD15_BCN_MI_Sources.csv_: Information on data sources

 _IMD15_BCN.csv_: This file refers to the 2015 high resolution layer of imperviousness density (IMD) for the selected territorial/spatial scales in Barcelona:
 * a)	Barcelona city boundaries (label: bcn_city)
 * b)	Barcelona metropolitan area, Àrea metropolitana de Barcelona (AMB) (label: bcn_amb)
 * c)	Barcelona greater city (Urban Atlas) (label: bcn_grc)
 * d)	Barcelona functional urban area (Urban Atlas) (label: bcn_fua)

 _IMD15_MI.csv_: This file refers to the 2015 high resolution layer of imperviousness density (IMD) for the selected territorial/spatial scales in Milan:
 * e)	Milan city boundaries (label: mi_city)
 * f)	Milan metropolitan area, Piano intercomunale milanese (PIM) (label: mi_pim)
 * g)	Milan greater city (Urban Atlas) (label: mi_grc)
 * h)	Milan functional urban area (Urban Atlas) (label: mi_fua)

_IMD15_BCN_MI.mpk_: the shareable project in Esri ArcGIS format including the HRL IMD data in raster format for each of the territorial boundaries as specified in letter a)-h).

Regarding the territorial scale as per letter f), the list of municipalities included in the Milan metropolitan area in 2016 was provided to me in 2016 from a person working at the PIM.

In the IMD15_BCN.csv and IMD15_MI.csv, the following columns are included:

 * Level: the territorial level as defined above (a)-d) for Barcelona and e)-h) for Milan);
 * Value: the 101 values of imperviousness density expressed as a percentage of soil sealing (0-100%: 0% means fully non-sealed area; 100% means fully sealed area);
 * Count: the number of 20x20mt cells corresponding to a certain percentage of soil sealing or imperviousness;
 * Km2: the conversion of the 20x20mt cells into squared kilometres (km2) to facilitate the use of the dataset.

 Further information on the Dataset

This dataset is the result of a combination between different databases of different types and that have been downloaded from different sources. Below, I describe the main steps in data management that resulted in the production of the dataset in an Esri ArcGIS (ArcMap, Version 10.7) project.

 1. The high resolution layer (HRL) of the imperviousness density data (IMD) for 2015 has been downloaded from the official website of Copernicus. At the time of producing the dataset (April/May 2021), the 2018 version of the IMD HRL database was not yet validated, so the 2015 version was chosen instead. The type of this dataset is raster.

 2. For both Barcelona and Milan, shapefiles of their administrative boundaries have been downloaded from official sources, i.e. the ISTAT (Italian National Statistical Institute) and the ICGC (Catalan Institute for Cartography and Geology). These files have been reprojected to match the IMD HRL projection, i.e. ETRS 1989 LAEA.

 3. Urban Atlas (UA) boundaries for the Greater Cities (GRC) and Functional Urban Areas (FUA) of Barcelona and Milan have been checked and reconstructed in Esri ArcGIS from the administrative boundaries files by using a Eurostat correspondence table. This is because at the time of the dataset creation (April/May 2021), the 2018 Urban Atlas shapefiles for these two cities were not fully updated or validated on the Copernicus Urban Atlas website. Therefore, I had to re-create the GRC and FUA boundaries by using the Eurostat correspondence table as an alternative (but still official) data source. The use of the Eurostat correspondence table with the codes and names of municipalities was also useful to detect discrepancies, basically stemming from changes in municipality names and codes and that created inconsistent spatial features. When detected, these discrepancies have been checked with the ISTAT and ICGC offices in charge of producing Urban Atlas data before the final GRC and FUA boundaries were defined.

Steps 2) and 3) were the most time consuming, because they required other tools to be used in Esri ArcGIS, like spatial joins and geoprocessing tools for shapefiles (in particular dissolve and area re-calculator in editing sessions) for each of the spatial/territorial scales as indicated in letters a)-h).

Once the databases for both Barcelona and Milan as described in points 2) and 3) were ready (uploaded in Esri ArcGIS, reprojected and their correctness checked), they have been ‘crossed’ (i.e. clipped) with the IMD HRL as described in point 1) and a specific raster for each territorial level has been calculated. The procedure in Esri ArcGIS was the following:

 * Clipping: Arctoolbox > Data management tools > Raster > Raster Processing > Clip. The ‘input’ file is the HRL IMD raster file as described in point 1) and the ‘output’ file is each of the spatial/territorial files. The option "Use Input Features for Clipping Geometry (optional)” was selected for each of the clipping.
 * Delete and create raster attribute table: Once the clipping has been done, the raster has to be recalculated first through Arctoolbox > Data management tools > Raster > Raster properties > Delete Raster Attribute Table and then through Arctoolbox > Data management tools > Raster > Raster properties > Build Raster Attribute Table; the "overwrite" option has been selected.

 Other tools used for the raster files in Esri ArcGIS have been the spatial analyst tools (in particular, Zonal > Zonal Statistics). As an additional check, the colour scheme of each of the newly created raster for each of the spatial/territorial attributes as per letters a)-h) above has been changed to check the consistency of its overlay with the original HRL IMD file. However, a perfect match between the shapefiles as per letters a)-h) and the raster files could not be achieved since the raster files are composed of 20x20mt cells.

The newly created attribute tables of each of the raster files have been exported and saved as .txt files. These .txt files have then been copied in the excel corresponding to the final published dataset.

barcelona

display_name: barcelona

copernicus

display_name: copernicus

imperviousness

display_name: imperviousness

milan

display_name: milan

soil-sealing

display_name: soil-sealing

urban-atlas

display_name: urban-atlas

educ

en: Education, culture and sport

en: Soil sealing Barcelona and Milan different territorial levels


maintainer: Chunhong Zhao

en: land-dynamics land-use landscape urban-morphology austin

en: Land use dynamics in Austin since 1992

en: The present dataset is part of the published scientific paper Zhao C, Weng Q, Hersperger A M. Characterizing the 3-D urban morphology transformation to understand urban-form dynamics: a case study of Austin, Texas, USA. Landscape and urban planning, 2020, 203:103881. The overall objective of this paper is to understand urban form dynamics in the Austin metropolitan area for the periods 2006–2011 and 2011–2016. The study also aims to understand to what extent the changes in the built environment (in terms of ‘efficient growth’ versus ‘inefficient growth’) from the 1990s to 2016 in the Austin metropolitan area corresponded with ‘compact and efficient growth’ planning policy documents. The UMT distribution can be found in the paper. The area of transitioning UMT was provided in Table 2 and Table 3 can be found in the Appendix of the paper. A protocol was developed to perform the content analysis of the strategic plans and gather the data. The detailed list of protocol items can be found in Appendix B of the paper. This study demonstrates the advantage of applying Lidar data to characterize 3-D urban morphology type (UMT) transition and understand its dynamics, which helps develop a comprehensive understanding of the urbanization process and provides a tool for planning intentions and policies evaluation on urban development over time. The UMT maps can be found in Appendix A of the paper. The Lidar point datasets and the 30 × 30 m National Land Cover Database (NLCD) are the two main data sources of UMT mapping. Lidar datasets were gathered from different projects that had been conducted and collected by state agencies and other organizations between 2007 and 2017. Table A1 in the appendix in the paper shows the accuracies and acquisition parameters of the Lidar projects from 2007 to 2017. Land use/cover dynamics in Austin metropolitan area dataset provides Land use/cover patterns in the years 1992, 2001, 2004, 2006, 2008, 2011, 2013, 2016 with a spatial resolution of 30 meters. Since NLCD 1992 used a different classification system for the urban land classes, we first reclassified the NLCD 1992 using a customized Arcpy package.

austin

display_name: austin

land-dynamics

display_name: land-dynamics

land-use

display_name: land-use

landscape

display_name: landscape

urban-morphology

display_name: urban-morphology

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: Land use/cover dynamics in Austin metropolitan area since 1992


photogrammetric-drone-data-dorfberg

maintainer: Yves Bühler

title_for_slug: photogrammetric-drone-data-dorfberg

en: drone-data snow-depth-mapping snow-distribution photogrammetry orthophoto uav dsm snow

en: Data Access

en:
- 3cm-resolution Orthophoto
- 10cm-resolution DSM
- 10cm-resolution snow depth raster
- original RGB images
- Agisoft Metashape Report

All products are also available as a cloud optimized geotiff.

en: The data was collected with a Wingtra Gen II drone and a Sony RX1R II sensor. In total, 10 flights were conducted at different dates, both in summer and winter. A DSM, an orthophoto, a snow depth raster and the original drone images from every flight are available at a high resolution (10cm and 3cm, respectively).

drone-data

display_name: drone-data

dsm

display_name: dsm

orthophoto

display_name: orthophoto

photogrammetry

display_name: photogrammetry

snow

display_name: snow

snow-depth-mapping

display_name: snow-depth-mapping

snow-distribution

display_name: snow-distribution

uav

display_name: uav

en: EnviDat: the Environmental Data Portal

en: Photogrammetric Drone Data Dorfberg


en: remote-sensing swe snow-depth snow uas

en: .tfw file

en: coordinates for the .tif raster

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

media_type: image/tiff

en:
10 cm resolution
15 cm accuracy

format: TIFF

mimetype: image/tiff

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

A detailed description is given here:

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

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

remote-sensing

display_name: remote-sensing

snow-depth

display_name: snow-depth

swe

display_name: swe

uas

display_name: uas

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


maintainer: Henning Löwe

en: discrete-element-modelling spherical-model-snow unconfined-compression-experiments microstructure 3d-ct-images snow-mechanics

en: Data Structure (Python Notebook - view online)

en: Data Structure (Python Notebook)

en:
For the investigation of microstructural and mechanical properties of snow unconfined compression experiments and 3D computed tomography (CT) imaging were performed on sintered rounded grain snow and spherical model snow. The spherical model snow was generated to create geometrically simplified, well-defined microstructures for calibration of numerical models, such as discrete element models (DEM) in which the microstructure is represented by spherical particles.
In the experiments, microstructural variation was created by varying the sintering time (contact size) and the density of the ice sphere samples (number of contacts). The 3D CT images allow for a complete reconstruction of the entire experimental sample (cylindrical sample dimension: diameter = 33.6 mm; height = 14 mm).

discrete-element-modelling

display_name: discrete-element-modelling

microstructure

display_name: microstructure

snow-mechanics

display_name: snow-mechanics

spherical-model-snow

display_name: spherical-model-snow

unconfined-compression-experiments

display_name: unconfined-compression-experiments

educ

en: Education, culture and sport

en: Unconfined compression experiments and 3D CT images of spherical model snow and RG snow samples


maintainer: Anita C Risch

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

en: Risch et al. GEB_data.xlsx

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

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

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

Paper Citation:

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

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

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

anthropogenic-management

display_name: anthropogenic-management

climate

display_name: climate

herbivore-exclusion

display_name: herbivore-exclusion

microbial-biomass-carbon

display_name: microbial-biomass-carbon

microbial-respiration

display_name: microbial-respiration

nutrient-addition

display_name: nutrient-addition

soil-properties

display_name: soil-properties

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.