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

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

educ

en: Education, culture and sport

en: Elk and bison carcasses in Yellowstone, USA


maintainer: EnviDat Support

en: avalanche-accidents number-of-caught-persons number-of-dead-persons avalanche-accident-statistics avalanche-fatalities

state: active

type: dataset

en: Fatal avalanche accidents in Switzerland since 1995-1996

media_type: text/csv

format: CSV

mimetype: text/csv

en:
Attention: this data is not updated after 2022 anymore.
This data collection contains information concerning all accidents by snow avalanches causing at least one fatality in Switzerland. The data set commences on 01/10/1995. After the completion of a hydrological year, the new data is added.
The following information is provided:
 * avalanche identifier
 * date of the accident
 * accuracy of the date in range of days before and after
 * canton
 * name of the locality
 * start zone of the avalanche
  * coordinates (Swiss coordinate system, approximately in middle of start zone)
  * accuracy of the coordinates in meters
 * elevation (in meteres above sea level, app. in middle of start zone)
 * slope aspect (main orientation of start zone)
 * slope inclination (in degree, steepest point within start zone)
 * number of dead persons
 * number of caught persons
 * number of fully buried persons
 * forecasted avalanche danger level
 * activity/location of the accident party at the time of the incident

avalanche-accident-statistics

display_name: avalanche-accident-statistics

avalanche-accidents

display_name: avalanche-accidents

avalanche-fatalities

display_name: avalanche-fatalities

number-of-caught-persons

display_name: number-of-caught-persons

number-of-dead-persons

display_name: number-of-dead-persons

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

it: Istruzione, cultura e sport

envidat

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

en: EnviDat: the Environmental Data Portal

it: EnviDat: il portale dei dati ambientali

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

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

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

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

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester


en: snow-stability-evaluation avalanche snowpack avalanche-forecast snow

en: README

media_type: text/plain

en: README file explaining data structure and content

format: TXT

mimetype: text/plain

en: SNOWPACK simulations

media_type: application/zip

en: This dataset contains all SNOWPACK simulations, including input files, setup and output files.

format: ZIP

mimetype: application/zip

en: Stability tests

en: Stability tests including csv-files with CT, ECT and PST experiments for each year and each field site.

en: Extracted data

en: Extracted data contain validation and results published in Richter et al. (2019).

en: Manually observed snow profiles

en: Manually observed snow profiles as csv-files containing depth [cm], thickness [cm], density [kg m-3], hand hardness, first grain type, second grain type, first grain size [mm], and second grain size [mm] of each layer.

en:
To validate the critical crack length as implemented in the snow cover model SNOWPACK, PST experiments were conducted for three winter seasons (2015-2017) at two field site above Davos, Switzerland. This dataset contains manually observed snow profiles and stability tests. Furthermore, corresponding SNOWPACK simulations are included.
These data were analyzed and results were published in Richter et al. (2019).

Please refer to the Readme file for further details on the data.

These data are the basis of the following publication:
Richter, B., Schweizer, J., Rotach, M. W., and van Herwijnen, A.: Validating modeled critical crack length for crack propagation in the snow cover model SNOWPACK, The Cryosphere, 13, 3353–3366, https://doi.org/10.5194/tc-13-3353-2019, 2019.

avalanche

display_name: avalanche

avalanche-forecast

display_name: avalanche-forecast

snow

display_name: snow

snow-stability-evaluation

display_name: snow-stability-evaluation

snowpack

display_name: snowpack

en: Validating and improving the critical crack length in SNOWPACK


en: land-use-change spatial-planning digital-zoning-plans land-use-simulations land-use-transformations

format: SERVICE

en: Legend simulated maps

en:
Table of Content: 1. General context of the data set "LUSzoning”; 2. Background and aims of the study using the data set LUSzoning; 3. The data set LUSzoning.
1. General context of the data set "LUSzoning".
The data set "LUSzoning" stands for Land-use simulations integrating zoning regulations in Spanish functional urban areas. The data set has been generated as part of the CONCUR research project (https://www.wsl.ch/en/projects/concur.html) led by Dr. Anna M. Hersperger and funded by the Swiss National Science Foundation (ERC TBS Consolidator Grant (ID: BSCGIO 157789) for the period 2016-2021. The CONCUR research project is interdisciplinary and aims to develop a scientific basis for adequately integrating spatial policies (in this case, digital zoning plans) into quantitative land-change modelling approaches at the urban regional level.
2.	Background and aims of the study using the data set “LUSzoning”.
As part of the CONCUR project, a specific task was to integrate planning spatial policies in land-change modelling. Planning can be implemented in modelling using either hard or gradual restrictions. Different studies have addressed the inclusion of spatial planning policies in land-use change modelling. However, the integration of zoning constraints is generally established as hard or Boolean-based restrictions (e.g., whether urban development is allowed or not), while not accounting for the spatial heterogeneity or gradual characteristics within planning zones (e.g., whether planning regulations allow low, medium or high urban density), though these could improve real patterns simulations in urban areas. We assume Spanish General Zoning plans were suitable to explore the integration of planning into land-change modelling as soft constrains because they define land-use intensities in the buildable zoning areas. In light of the above considerations, the overall aim of the study was to model urban land-use changes using a multi-scenario approach that integrates digitized zoning plans for the Functional Urban Areas (FUAs) of Madrid, Barcelona, Valencia, and Zaragoza. The following specific objectives were addressed: i) to analyse the role of planning by defining three future scenarios that integrate digitized zoning plans and one scenario that assumes almost no planning intervention; ii) to introduce zoning constraints that reflect different degrees of urban densities; iii) to generate a transferable spatially-explicit modelling framework to integrate planning into land-use change simulations.
Four future land-use demands scenarios were defined for the FUAs. Storylines were created considering probable development scenarios related to zoning plans, current Spanish legislation and sustainability goals defined along two axes: a high market-oriented vs. high planning-intervention axis, and an axis of short-term economic growth vs. long-term sustainable growth. The sustainable development scenario (S1) is characterized by low gross floor area (GFA) growth that is limited to areas that are currently under development according to zoning plans. The business-as-usual scenario (S2) is characterized by medium GFA growth in the range of on-going trends. The strong development scenario (S3) is characterized by high GFA growth rates. Growth is restricted to buildable areas without urbanization project designated in zoning plans. The unrestricted development scenario (S4) prioritizes a high degree of market liberalization characterized by high GFA growth that surpasses population demands. S4 follows a rapid economic growth pattern with almost no planning intervention.
3.	 The data set “LUSzoning”.
The dataset includes 16 .asc raster layers providing the simulated land-uses under four defined scenarios for Barcelona, Madrid, Valencia and Zaragoza Functional Urban Areas (FUAs) for 2030. The simulated raster layers were created using CLUMondo simulation framework and have a spatial resolution of 30m. The .asc layers name include the name of the FUA and scenario number. For example, the output from simulating the urban growth for the city of Zaragoza under Scenario 2 is named “Zaragoza_S2.tif”. Furthermore, a .txt file named “Legend.txt” includes the numeric value of the land-use and the category of land-use that represents to interpret the .asc raster layers. The name of the land-use classes is a reclassification of the Urban Atlas 2012 land-use classes within the four Spanish FUAs analyzed.

digital-zoning-plans

display_name: digital-zoning-plans

land-use-change

display_name: land-use-change

land-use-simulations

display_name: land-use-simulations

land-use-transformations

display_name: land-use-transformations

spatial-planning

display_name: spatial-planning

educ

en: Education, culture and sport

en: LUSzoning: Land-use simulations integrating zoning regulations in Spanish functional urban areas


en: switzerland alpine3d-simulations snow climate-change graubunden

en: Configuration file

en: Simulation data

en: Reference data as well as the 2050 and 2099 scenarios for ARPEGE and ETHZ other RCMs available on demand)

en:
This is the simulation dataset from _"Response of snow cover and runoff to climate change in high Alpine catchments of Eastern Switzerland"_, M. Bavay, T. Grünewald, M. Lehning, Advances in Water Resources 55, 4-16, 2013

A model study on the impact of climate change on snow cover and runoff has been conducted for the Swiss Canton of Graubünden. The model Alpine3D has been forced with the data from 35 Automatic Weather Stations in order to investigate snow and runoff dynamics for the current climate. The data set has then been modified to reflect climate change as predicted for the 2021-2050 and 2070-2095 periods from an ensemble of regional climate models.

The predicted changes in snow cover will be moderate for 2021-2050 and become drastic in the second half of the century. Towards the end of the century the snow cover changes will roughly be equivalent to an elevation shift of 800 m. Seasonal snow water equivalents will decrease by one to two thirds and snow seasons will be shortened by five to nine weeks in 2095.

Small, higher elevation catchments will show more winter runoff, earlier spring melt peaks and reduced summer runoff. Where glacierized areas exist, the transitional increase in glacier melt will initially offset losses from snow melt. Larger catchments, which reach lower elevations will show much smaller changes since they are already dominated by summer precipitation.

climate-change

display_name: climate-change

graubunden

display_name: graubunden

switzerland

display_name: switzerland

envidat

en: EnviDat: the Environmental Data Portal

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

en: Alpine3D simulations of future climate scenarios for Graubunden


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

media_type: image/tiff

en: Frost change frequency modified from CHELSA

format: TIFF

en: Mean annual temperature modified from CHELSA

en: Soil pH modified from SoilGrids

en: Mean diurnal range modified from CHELSA

en: Precipitation seasonality from CHELSA

en: Precipitation of driest quarter modified from CHELSA

en: Soil clay content modified from SoilGrids

en: Soil organic carbon content modified from SoilGrids

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

aridity

display_name: aridity

frost-change-velocity

display_name: frost-change-velocity

ground-temperature

display_name: ground-temperature

precipitation

display_name: precipitation

sdms

display_name: sdms

soil

display_name: soil

en: Environmental layers for SDM simulations (GDPlants)


maintainer: Athanasios Nenes

en: dmt-ccn-chamber weissfluhjoch sigma-w hygroscopicity raclets cloud-modeling ccn supersaturation cloud-droplets

en: Predicted cloud droplet number concentrations Weissfluhjoch

en:
In-situ measurements of aerosol size distributions along with updraft velocities collected from a wind Doppler Lidar are coupled with a state-of-the art droplet activation parameterization (Morales and Nenes, 2014) to determine the cloud droplet number and maximum in-cloud supersaturation that can potentially form over the Weissfluhjoch station between February 24 and March 8 2019. For this, a hygroscopicity parameter of 0.25 is assumed, according to the CCN-derived values. The data file contains:

* Date and time (DD/MM/YYYY hh:mm:ss) in UTC
* The wind-lidar derived standard deviation of updraft velocities (sigma-w) in m/s, used to apply the parameterization employing a single characteristic velocity (Morales and Nenes, 2010)
* The maximum predicted cloud parcel supersaturation (Smax%)
* The total in-situ aerosol number concentration (Naer) in cmˉ³, measured by the SMPS
* The predicted cloud droplet numbers (Nd) in cmˉ³
* The activation fraction (%)

en: Hygroscopicity parameter Weissfluhjoch

en:
The aerosol hygroscopicity parameter is derived combining the CCN and SMPS size distribution data collected between February 24 and March 8 2019 at Weissfluhjoch. The data file contains:

* Date and time (DD/MM/YYYY hh:mm:ss) in UTC
* The instrument supersaturation (SS%)
* The calculated hygroscopicity parameter (kappa)
* The calculated critical dry diameter of the particles (Dcr) in nm.

Missing values are indicated by "-999".

en: CCN data Weissfluhjoch

en:
In-situ CCN measurements were conducted between February 24 and March 8 2019 at Weissfluhjoch. The data file contains:

* Date and time (DD/MM/YYYY hh:mm:ss) in UTC
* Filtered CCN number concentrations (in cmˉ³) measured at 6 discrete supersaturations (SS) ranging from 0.09% to 0.74%. The name of each CCN column is indicative of the SS developed within the instrument.
* Instrument characteristics such as sample temperature (in °C) and pressure (in mbar)
* A code that refers to the reason of filtering.

Missing or filtered CCN values are indicated by "-999", while a declaration of each "Error_code" can be found in the following list:

_Error_code_:

0000: Measurements without error

1000: Period of transient operation after a supersaturation change

2000: Sharp change (decrease) in the room temperature housing the instrument, which induces a reset in column temperature

3000: Irregular values in the room temperature housing the instrument

9999: Missing data

en:
Cloud Condensation Nuclei (CCN) data: A Droplet Measurement Technologies (DMT) single-column continuous-flow streamwise thermal gradient chamber (CFSTGC; Roberts and Nenes, 2005) was deployed at the measurement site Weissfluhjoch (2700 m a.s.l., LON: 9.806475, LAT: 46.832964) to record the in-situ CCN number concentrations between February 24 and March 8 2019 for different supersaturations (SS). To account for the difference between the ambient (~735 mbar) and the calibration pressure (~800 mbar), the SS reported by the instrument is adjusted by a factor of 0.92. The CFSTGC was cycled between 6 discrete SS values ranging from 0.09% to 0.74%, producing a full CCN spectrum every hour. The raw CCN measurements are filtered to discount periods of transient operation and whenever the room temperature housing the instrument changed sufficiently to induce a reset in column temperature. Additional information can be found in Section 2.1.2 [here](https://acp.copernicus.org/preprints/acp-2020-1036/).

Hygroscopicity data: The CCN number concentration measurements were directly related to the size distribution and total aerosol concentration data measured by the Scanning Mobility Particle Size Spectrometer (SMPS) instrument at the same station (https://www.envidat.ch/dataset/aerosol-data-weissfluhjoch) to infer the particles hygroscopicity parameter (kappa). For each SMPS scan, the particles critical dry diameter (Dcr) is estimated by integrating backward the SMPS size distribution, until the aerosol number matches the CCN concentration observed for the same time period as the SMPS scan. Assuming the particle chemical composition is internally mixed, the kappa is determined from Dcr and SS, applying Köhler theory. Additional information can be found in Section 2.2 [here](https://acp.copernicus.org/preprints/acp-2020-1036/).

Predicted cloud droplet numbers: Droplet calculations are carried out with the physically based aerosol activation parameterization of Morales and Nenes (2014), employing the “characteristic velocity” approach of Morales and Nenes (2010). Aerosol size distribution observations required to predict the cloud droplet numbers and maximum in-cloud supersaturation are obtained from the SMPS instrument deployed at Weissfluhjoch. The required vertical velocity measurements are derived from the wind Doppler Lidar (https://www.envidat.ch/dataset/lidar-wind-profiler-data) deployed at Davos Wolfgang and are extracted for the altitude of interest, being 1100 m above ground level for Weissfluhjoch. Additional information can be found in Section 2.3 [here](https://acp.copernicus.org/preprints/acp-2020-1036/).

ccn

display_name: ccn

cloud-droplets

display_name: cloud-droplets

cloud-modeling

display_name: cloud-modeling

dmt-ccn-chamber

display_name: dmt-ccn-chamber

hygroscopicity

display_name: hygroscopicity

raclets

display_name: raclets

sigma-w

display_name: sigma-w

supersaturation

display_name: supersaturation

weissfluhjoch

display_name: weissfluhjoch

educ

en: Education, culture and sport

en: CCN, hygroscopicity, predicted cloud droplet numbers Weissfluhjoch


maintainer: Axel Volkwein

en: trajectory field-data field-tests experiment experiments rockfall mass-movements field-test natural-hazards

en: Overview of all tests

en:
This table gives an overview on all tests performed including
- block used
- start & end position
- number of jumps
- etc.
-

en: These data contain the point clouds obtained from laser scanning the test slope and the six test blocks used for the experiments. The coordinates of the slope are given in CH-1903 format. The coordinates of the block point clouds are kept in a lokal x/y/z coordinate system.

en: Block trajectories recorded from local positioning system

en: This file contains more than 80 rockfall trajectories projected onto the terrain, i.e. the positions of the single blocks over time during the single tests that were recorded by the local positioning system with about 10Hz. It contains the detected positions as well as a smoothened path (see Volkwein et al., Geosciences, 2018). The table also contains the derived velocities and accelerations along the slope (attention: these accelerations are the once that occur during an impact).

en: Jumps and impacts

en: This table contains data of more than 2200 of each jumps and impacts of the test blocks on a grassy slope. It gives there durations in combination with acting accerlartions and the rotational speed of the blocks. The data were retrieved using the block internal accelerometer and rotation sensors. The combination of these data with the measurement from the local positioning system allows the deduction of jump lengths, slope gradient etc.

en:
In summer 2014, 6 rock blocks between 20 and 80kg have been thrown in total 111 times down a slope at the Swiss Oberalppass close to the village Tschamut. The slope was mainly covered by grass and its lower part was flat and large enough to provide natural runouts of the single trajectories.

An extensive measurement program has been set up to measure the block trajectories: With surveyor's instruments the slope and the six used rock blocks were scanned and the start and end positions of each test were recorded. During the single events two cameras filmed the trajectories. A special sensor device located within the blocks recorded the acting accelerations and rotational speeds over time. Further, the device emitted a Wifi signal that got detected from eight receivers around the slope. Based on this signal the block position has been recorded over time.

The dataset contains all data that were gathered through above field campaign.

experiment

display_name: experiment

experiments

display_name: experiments

field-data

display_name: field-data

field-test

display_name: field-test

field-tests

display_name: field-tests

mass-movements

display_name: mass-movements

natural-hazards

display_name: natural-hazards

rockfall

display_name: rockfall

trajectory

display_name: trajectory

envidat

en: EnviDat: the Environmental Data Portal

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

en: Repetitive trajectory testing in Tschamut 2014


en: drifting-snow snow-saltation experiments water-isotopes metamorphism wind-tunnel cold-lab blowing-snow

en: Wind tunnel water vapour isotopic composition

en: The water vapour and its calibrated isotopic composition sampled inside the wind tunnel at 46cm height during the experiments. The vapour was pumped through a heated copper tube at a rate of 0.2 l/min to the CRDS instrument located in a temperate room outside the cold laboratory. The data has 1 second resolution with the timestamp in UTC (CRDS reference time). The measured variables are the calibrated d18O (‰), dD (‰) and d-excess (‰) signals and the uncalibrated ("H2Oraw", ppmv) and calibrated ("H2Ocali",ppmv) humidity readings of the CRDS and converted to specific humidity ("spechum", g/kg). The isotope data was corrected for its humidity-isotope dependency and referenced against the VSMOW-SLAP scale. Calibration runs were performed daily and details can be found in Wahl et al. (2024) including the Supplementary Materials.

en: Snow samples - microCT measurements

en: The data of the snow samples that were analysed in a Scanco® CT-40 scanner using a sample holder of 15 mm diameter which allowed for a high-resolution analysis with voxel size of 8 μm. The timestamp (local time, Davos, Switzerland) of snow sampling is given as well as the type of snow sample (see Snow samples - stable water isotopic composition dataset), the corresponding wind tunnel experiment number ("ExpNo", Wahl et al. (2024), Table 1) and the time since the snow introduction to the wind tunnel. The measured SSA is given in units 1/mm. The column "hist_filename" lists the name of the corresponding .txt file which contains the sphere size statistics of the snow samples. The sphere size statistics files are compressed in a .zip file.

en:
The dataset contains the stable water isotopic composition (d18O, dD. d-excess) of the snow samples taken during the ring wind tunnel experiments in January and May 2023 at the SLF cold laboratory facilities in Davos. The snow samples were collected in air-tight WhirlPaks and kept frozen at -20°C until they were melted in the closed sampling bags, filtered and 1.8ml were transferred into a 2ml glass vial and shipped to the Paul Scherrer Institut in Villigen, Switzerland where they were analysed following established laboratory procedures. The combined uncertainty on the individual isotopic composition is <0.5‰ for dD and <0.1‰ for d18O which yields an uncertainty of 0.9‰ on d-excess.
Timestamps are given for local Davos time (LT) and synchronised to the CRDS instrument which is in UTC. The "ExpNo" column refers to the experiment numbers given in Wahl et al. (2024), Table 1. The "type" of the snow samples is: "ini" = initial snow before introduction to the wind tunnel, "airborne" = snow sampled out of the air, "acc_in" = inner part of accumulated snow, "acc_out" = outer surface layer of accumulated snow.

en: Meteorological variables measured inside the wind tunnel during the January 2023 experiments

en: The dataset contains the 1-minute averages of the calibrated meteorological variables: temperature (T, °C), relative humidity (RH, %), relative humidity referenced against saturation over ice (RHice, %), specific humidity (q, g/kg) and wind speed (WS, m/s) as measured in the observation section of the wind tunnel at different heights. The measurement heights correspond to: "Top"=42.5cm, "Mid"=30cm, "Bot"=15cm, "PT100_Ta"=3cm. "WS" was measured at 36.5cm height. The timestamp is in local time (UTC+1) and indicates the start of the averaging period.

en: Meteorological variables measured inside the wind tunnel during the May 2023 experiments

en: The dataset contains the 1-minute averages of the calibrated meteorological variables: temperature (T, °C), relative humidity (RH, %), relative humidity referenced against saturation over ice (RHice, %), specific humidity (q, g/kg) and wind speed (WS, m/s) as measured in the observation section of the wind tunnel at different heights and the room temperature of the laboratory measured outside the wind tunnel (Lab_Ta, °C). The measurement heights correspond to: "Top"=42.5cm, "Mid"=30cm, "Bot"=15cm, "PT100_Ta"=3cm. "WS" was measured at 36.5cm height. The time was synchronised to the concurrent CRDS measurements (which are UTC) and adjusted for the time lag resulting from pumping vapour from inside the wind tunnel to the CRDS. The timestamps indicate the beginning of the averaging period.

en: Snow samples - sphere size statistics files

en: The individual measurement reports of the microCT scan for each snow sample as combined zip file. The file names and the sample number indicate which snow sample it is, see "Snow samples - microCT measurements". Each file contains the statistics and the count of fitted sphere size diameters and the respective share (in-%) in the frequency distribution. Immediately after snow sampling, a small portion of snow was transferred to the cylindrical μCT sample holder with a diameter of 15 mm and analysed at -20°C using a Scanco® CT-40 scanner. The analysis followed Hildebrand and Rüegsegger (1997).

en: This dataset collection contains all datasets collected during the ring wind tunnel experiments in January and May 2023 at the SLF cold laboratory facilities in Davos. A full description of the experiment set-up can be found in Wahl et al. (2024). The collection contains data from 19 experiments. The dataset collection entails measurements of the stable water isotopic composition of snow samples and the water vapour inside the wind tunnel, measurements of the meteorological variables inside the wind tunnel and snow sample properties as measured with microCT measurements (sphere size distribution and specific surface area (SSA)).

blowing-snow

display_name: blowing-snow

cold-lab

display_name: cold-lab

drifting-snow

display_name: drifting-snow

metamorphism

display_name: metamorphism

snow-saltation

display_name: snow-saltation

water-isotopes

display_name: water-isotopes

wind-tunnel

display_name: wind-tunnel

en: Ring wind tunnel experiments - airborne snow metamorphism and stable water isotopes


maintainer: Matthias Haeni

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

en: Forest station data (OTB)

en: Data file.

en: Site photo of the field station (OTF)

en: Site photo of the field station from below (viewing upwards the slope).

mimetype: image/jpeg

en: Site aerial photo (Othmarsingen)

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

en: Site photo of the forest station (OTB)

en: Site photo of the forest station from above (viewing downwards the slope).

en: Data plot of the field station (OTF)

media_type: application/pdf

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

format: PDF

mimetype: application/pdf

en: Field station data (OTF)

en: Data plot of the forest station (OTB)

en: High quality meteorological data are needed for long-term forest ecosystem research, particularly in the light of global change. The long-term data series published here comprises almost 20 years of measurements for two meteorological stations in Othmarsingen in Switzerland where one station is located within a natural deciduous forest stand (OTB) with European beech (_Fagus sylvatica_; 120-140 yrs) and lime trees (_Tilia sp._; 120-140 yrs) as dominant tree species. A second station is situated in the very vicinity outside of the forest (field station, OTF). The meteorological time series are presented in hourly time resolution of air temperature, relative humidity, precipitation, photosynthetically active radiation (PAR) and wind speed. Othmarsingen is part of the Long-term Forest Ecosystem Research Programme (LWF) established and maintained by the Swiss Federal Research Institute WSL.

climate-change

display_name: climate-change

long-term

display_name: long-term

meteo

display_name: meteo

meteorology

display_name: meteorology

par

display_name: par

precipitation

display_name: precipitation

relative-humidity

display_name: relative-humidity

swiss-forest-lab

display_name: swiss-forest-lab

temperature

display_name: temperature

wind-speed

display_name: wind-speed

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

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


en: cloud-droplets pmcamx-uf total-number modeling nucleation

en: data_JRL_New particle formation events can reduce cloud droplets in boundary layer clouds at the continental scale.Paper_2023GL106182.zip

en: This repository contains PMCMx-UF model outputs used for the paper: D. Patoulias K. Florou S. N. Pandis and A. Nenes: New particle formation events can reduce cloud droplets in boundary layer clouds at the continental scale, Geophysical Research Letters, in review, 2023.

cloud-droplets

display_name: cloud-droplets

modeling

display_name: modeling

nucleation

display_name: nucleation

pmcamx-uf

display_name: pmcamx-uf

total-number

display_name: total-number

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: Dataset for the publication New particle formation events can reduce cloud droplets in boundary layer clouds at the continental scale


maintainer: Zhao

en: land-dynamics driving-factors urban-land-use planning austin

en: Potential driving factors of urban transformations of Austin over 25 years

en: In this study, the Austin metropolitan area, Texas, U.S., one of the fastest urban transformations and transformations regions, is selected to test the hypothesis that spatial planning and policies are important factors of urban transformations. Despite ample previous work in understanding the interactions between human and urban form transformation at specific areas, the actual interventions and outcomes of planning and policies (e.g., ‘smart growth’) on urban forms have been poorly measured. In this study, the potential influencing factors of urban transformations of Austin over 25 years were selected and collected.

austin

display_name: austin

driving-factors

display_name: driving-factors

land-dynamics

display_name: land-dynamics

planning

display_name: planning

urban-land-use

display_name: urban-land-use


maintainer: Paraskevi Georgakaki

en: wrf liquid-water-content secondary-ice-production ice-crystal-number-concentration swiss-alps blowing-snow jungfraujoch ice-water-content

en: Observational data

en:
The .zip file contains 6 .csv files including all microphysical measurements taken at the Jungfraujoch station during the two periods of interest:

* 25-28 January 2014 : northwesterly-winds (NW)
* 29 January - 1 February 2014 : southeasterly-winds (SE)

The total Ice Crystal Number Concentration (ICNC; reported in L⁻¹) is measured by the two-dimensional stereo hydrometeor spectrometer (2D-S; Lawson et al., 2006) over the size range 10-1280 μm. The total ice water content (IWC; reported in gm⁻³) is also derived by the 2D-S data using the Brown and Francis (1995) mass-diameter relationship. Finally, the quantification of the liquid water content (LWC; reported in gm⁻³) is based on the liquid droplet size distribution data derived from a DMT cloud droplet probe (CDP; Lance et al., 2010) over the size range 2-50 μm. The frequency of the measurements in the .csv files is 30 min, to match the output frequency of the model. Datetime is also reported in the files (DD/MM/YYYY h:mm) in UTC.

Acknowledgment: Collection and pre-processing of the microphysical measurements has been performed by Gary Lloyd (gary.lloyd@manchester.ac.uk). More details about the instrumentation during the CLACE2014 campaign can be found in Lloyd et al. (2015).

en: WRF outputs

en:
The .zip file contains all the NetCDF outputs produced by the mesoscale Weather Research and Forecasting (WRF) model (version: 4.0.1). Two case studies were simulated in this study, between 25-28 January 2014 (00:00 UTC) and 29 January - 1 February 2014 (00:00 UTC), respectively. The former is denoted as "NW" and the latter as "SE", owing to the prevailing wind conditions at the Jungfraujoch (JFJ) station.

The outputs come from the innermost domain of the simulations with a spatial resolution of 1 km and are extracted from the nearest to JFJ model grid point (46.55°N 7.98°E). The dimensions of each model variable is [145,100], which corresponds to time and height, respectively. The frequency of the outputs is 30 min. Note that the first 24 hours of each simulation (i.e., the first 48 outputs) were considered spin up time.

The following sets of sensitivity simulations were presented in this study:

* CNTRL_NW.nc, CNTRL_SE.nc : Sensitivity simulations performed with the default version of WRF
* DS_PHILL_NW.nc, DS_PHILL_SE.nc : Sensitivity simulations including a description of the droplet shattering process following Phillips et al. (2018)
* DS_SULL_NW.nc, DS_SULL_SE.nc : Sensitivity simulations including a description of the droplet shattering process following Sullivan et al. (2018)
* PHILL_NW.nc, PHILL_SE.nc : Sensitivity simulations including a description of the ice-ice collisional break-up process following Phillips et al. (2017)
* TAKAH_NW.nc, TAKAH_SE.nc : Sensitivity simulations including a description of the ice-ice collisional break-up process following Sullivan et al. (2018) based on Takahashi et al. (1995)
* TAKAHrim_NW.nc, TAKAHrim_SE.nc : Same set-up as before but applied only to collisions between rimed ice particles
* BIPS_Xk0_DY_NW.nc, BIPS_Xk0_DY_SE.nc : Sensitivity simulations including a source of Xk0 (either 10.0 or 100.0 L⁻¹) ice crystals with sizes of DY (either 10 or 100 μm) in the first atmospheric level of the model
* BIPS100_PHILL_NW.nc, BIPS100_PHILL_SE.nc : Sensitivity simulations where the combined effect of blowing snow and ice-ice collisional break-up is examined

en:
This repository contains all WRF model outputs and observational data sets used for the paper:
Georgakaki, P., Sotiropoulou, G., Vignon, É., Billault-Roux, A.-C., Berne, A., and Nenes, A.: Secondary ice production processes in wintertime alpine mixed-phase clouds, Atmos. Chem. Phys. Discuss., https://doi.org/10.5194/acp-2021-760, in review, 2021.

blowing-snow

display_name: blowing-snow

ice-crystal-number-concentration

display_name: ice-crystal-number-concentration

ice-water-content

display_name: ice-water-content

jungfraujoch

display_name: jungfraujoch

liquid-water-content

display_name: liquid-water-content

secondary-ice-production

display_name: secondary-ice-production

swiss-alps

display_name: swiss-alps

wrf

display_name: wrf

en: Data for the publication "Secondary ice production processes in wintertime alpine mixed-phase clouds"


en: snow-density weisfluhjoch snow-stratigraphy snow-ice-temperature snow-water-equivalent

en: WFJ Snow Profiles 1999-2014

en:
Biweekly manual snow profiles from Weissfluhjoch, Davos, Switzerland.
Version: 2015-09-14T16:25

en: README

en: README file giving use metadata.

en: Dataset of manual bi-weekly snow profiles from Weissfluhjoch, Davos, Switzerland. Typical snow profile measurements and observations are included (temperature, density, grain size, grain type, hardness, wetness), following the guidelines of the The International Classification for Seasonal Snow on the Ground (ICSSG) [Fierz, C., Armstrong, R.L., Durand, Y., Etchevers, P., Greene, E., McClung, D.M., Nishimura, K., Satyawali, P.K. and Sokratov, S.A. 2009. The International Classification for Seasonal Snow on the Ground. IHP-VII Technical Documents in Hydrology N°83, IACS Contribution N°1, UNESCO-IHP, Paris].

snow-density

display_name: snow-density

snow-ice-temperature

display_name: snow-ice-temperature

snow-stratigraphy

display_name: snow-stratigraphy

snow-water-equivalent

display_name: snow-water-equivalent

weisfluhjoch

display_name: weisfluhjoch

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Manual bi-weekly snow profiles from Weissfluhjoch, Davos, Switzerland


en: weissfluhjoch latent-heat stable-water-isotopes natural-tracer snow sublimation post-depositional-processes paleoclimate

en: Folder containing the scripts and data used in "Surface snow isotopic composition: post-depositional processes in an alpine snowpack" (Sadowski et al., 2022)

en: Stable water isotopes in snow

en: Stable water isotopes in snow measurements

en: Stable water isotopes in vapor

en:
Stable water isotopes in vapor measurements. Processed (calibrated) and raw data for February - June 2017. Raw data for January and February 2020.
For the time in between July 2017 and April 2020, raw vapor isotopic measurements are available, but no calibrations were performed during that time. For data delivery, please contact one of the authors.
This folder also contains the post-processing scripts for stable water isotopes in vapor measurements.

en: Weather

en: Meteorological data, snow height measurement, and SNOWPACK simulations covering the different field campaigns.

en:
Details on the data collected during the different field campaigns and available in this dataset.
Information on the sampling procedure of the snow isotopes samples.

en: Details about the different folders

en:
Notice: Changes to the dataset are still possible. Please do not use this dataset until the final publication with a DOI. Contact the authors if you have questions about this.

This dataset contains measurements of stable water isotopes in snow and vapor on the Weissfluhjoch from different field campaigns (Winter 2017 (Trachsel, 2019), January 2020, December 2020, and March 2021 (Sadowski et al., 2022).
Snow profiles and surface samples are available at different frequencies for each campaign. Please see "Data_description.pdf" for details.

Scripts and SNOWPACK simulations used in (Trachsel, 2019) and (Sadowski et al., 2022) are also provided.

latent-heat

display_name: latent-heat

natural-tracer

display_name: natural-tracer

paleoclimate

display_name: paleoclimate

post-depositional-processes

display_name: post-depositional-processes

snow

display_name: snow

stable-water-isotopes

display_name: stable-water-isotopes

sublimation

display_name: sublimation

weissfluhjoch

display_name: weissfluhjoch

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: Stable Water Isotopes in snow and vapor on the Weissfluhjoch


en: occupancy-model monitoring sessile-species swisslichens detectability inter-observer-difference bayesian-inference lichens detection-probability biodiversity observer-error cryptogams

en: Repository containing the full data, code, README and workflow files.

en:
The available lichen data consists of detection/nondetection data (1/0) of 373 tree-inhabiting (*corticolous*) lichen species from 416 plots surveyed 1-2 times. The lichen data were originally collected for the purpose of the Red List of epiphytic lichen species in Switzerland ([Scheidegger *et al.* 2002](https://www.bafu.admin.ch/bafu/de/home/themen/biodiversitaet/publikationen-studien/publikationen/rote-liste-gefaehrdete-arten-baum-erdbewohnende-flechten.html)), but updated to recent nomenclature for the purpose of this study.
This repository contains all the supporting data and R code for the paper:
von Hirschheydt, G., Kéry, M., Ekman, S., Stofer, S., Dietrich, M., Keller, C., Scheidegger, C. (2024) **Occupancy model reveals limited detectability of lichens in a standardised large-scale monitoring**. *Journal of Vegetation Science*.
Results and figures presented in the manuscript should be reproducible (with small differences in the latter digits due to stochasticity of the MCMC sampler) with the provided data and code.
The downloadable `.zip` folder has the following structure:
* `0_data/`
* `1_code/`
* `2_output/`
* `lichen_detectability.Rproj`
* `README.txt`
* `workflow.html`
- `workflow.Rmd`

The main folder and the three subordinate folders each have their own `README*.txt` file. These describe each available file in detail and should be consulted prior to using the data or running any code. The file `lichen_detectability.Rproj` stores the information about the R project. The user can open the project by clicking/double-clicking on this file which will automatically define the repository as working directory for the R session. If the user does not use RStudio/Posit, they may have to set the working directory manually to the stored location in the R files. The files `workflow.*` guide the user through the analysis (`1_code/*.R`) in the correct order so that they can:
- bundle the cleaned data into a data list readable for JAGS
- fit the multi-species occupancy model to the data and store the output
- assess the goodness-of-fit of the model to the data
- conduct a prior sensitivity analysis with 2 additional sets of priors
- extract the summary statistics reported in the manuscript and supplementary materials
- generate the figures shown in the manuscript and supplementary materials

bayesian-inference

display_name: bayesian-inference

biodiversity

display_name: biodiversity

cryptogams

display_name: cryptogams

detectability

display_name: detectability

detection-probability

display_name: detection-probability

inter-observer-difference

display_name: inter-observer-difference

lichens

display_name: lichens

monitoring

display_name: monitoring

observer-error

display_name: observer-error

occupancy-model

display_name: occupancy-model

sessile-species

display_name: sessile-species

swisslichens

display_name: swisslichens

en: Repeated detection-nondetection data of corticolous lichens from a standardised monitoring across Switzerland


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

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

format: ODS

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

fem

display_name: fem

metamorphism

display_name: metamorphism

microstructure

display_name: microstructure

time-series

display_name: time-series

x-ray-tomography

display_name: x-ray-tomography

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

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


en: invertebarates trophic-networks vascular-plants switzerland vertebrates food-webs metaweb

en: trophiCH

en: We provide seven datasets: 1) the metaweb, 2) the taxa checklist, 3) the data source meta-dataset, 4) the list of generalist basal and predator families and polylectic species with citations and 5) a dataset with citations for the inferences of missing predators, 6) a dataset with citations for the parallel inference of diets from similar species and 7) a list of other existing metawebs. The empirical dataset is 00_empirical_metaweb.csv. To reproduce the final metaweb (01_final_metaweb.csv), start the R Project file: trophich.Rproj. We provide five scripts, accompanying functions, and the raw data required to run these scripts to reproduce the taxonomic expansion and validation of the datasets. In the first script (01_inferring_interactions.R), we infer interactions using genus and family level interactions and for basal feeding groups (see Methods: Taxonomic expansion). In the second script (02_inferring_interactions_special_cases.R), we infer further interactions for a few special cases with detailed explanations. In the third script (03_metaweb_comparisons.Rmd), we provide the statistical comparisons between our metaweb and other empirical metawebs as an R Markdown document. We additionally provide a a Python Jupyter Notebook document, outlining the error validation of the data extraction process (04_error_validation.ipynb and an accompanying .html file). Finally, we provide a script to reproduce figure 1 from the data paper (05_metaweb_summary_figure_1.R).

en:
Environmental pressures on species can cascade within food webs and even extend beyond individual ecosystems to interconnected systems at large spatial scales. To facilitate the exploration of these dynamics, we construct a data-based national trophic meta-food web (henceforth metaweb), that includes well-documented vertebrates, invertebrates, and vascular plants within Switzerland's national boundaries, and compiles 160 years of ecological knowledge. We additionally use a combination of taxonomic and geographic information to infer further species-level interactions. Our comprehensive dataset catalogues 1,112,073 trophic interactions involving 23,151 species and 125 feeding guilds (e.g. detritivores, fungivores, etc). While explorations of large-scale food webs in space and time have often relied on modelling approaches due to limited data availability, this empirically based metaweb paves the way for data-driven large-scale studies of real-world food webs. By integrating the metaweb with knowledge of local species assemblages, future studies can gain insights into the impact of global change drivers, including climate change and land-use intensification, on food web structures across spatial scales.
**Please note that this dataset is under embargo until publication of the paper.**

food-webs

display_name: food-webs

invertebarates

display_name: invertebarates

metaweb

display_name: metaweb

switzerland

display_name: switzerland

trophic-networks

display_name: trophic-networks

vascular-plants

display_name: vascular-plants

vertebrates

display_name: vertebrates

en: trophiCH v1 - a food web for Switzerland


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

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

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: An ice oxygen K-edge NEXAFS spectroscopy data set on gas-phase processing


en: incoming weissfluhjoch shortwave-radiation outgoing longwave-radiation

en: Parameter file

en: DATA ACCESS

en: Processing doc

en: Metadata

en:
Corrected incoming and outgoing shortwave and longwave 2 min radiation measurements at the Weissfluhjoch summit, Davos, CH.

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

incoming

display_name: incoming

longwave-radiation

display_name: longwave-radiation

outgoing

display_name: outgoing

shortwave-radiation

display_name: shortwave-radiation

weissfluhjoch

display_name: weissfluhjoch

en: ASRB_WFJ: Shortwave and longwave radiation measurements (2 min) at the Weissfluhjoch research site, Davos


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

en: Site photo of the field station (SCF)

en: Site photo of the field station from above (viewing downwards the slope).

en: Forest station data (SCB)

en: Data file.

en: Data plot of the forest station (SCB)

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

en: Site photo of the forest station (SCB)

en: Site photo of the forest station from below (viewing upwards the slope).

en: Further details on the data repository

en: Meta data documentation.

en: Field station data (SCF)

en: Site aerial photo (Schänis)

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

en: Data plot of the field station (SCF)

en: High quality meteorological data are needed for long-term forest ecosystem research, particularly in the light of global change. The long-term data series published here comprises almost 20 years of measurements for two meteorological stations in Schänis in Switzerland where one station is located within a natural mixed forest stand (SCB) with European beech (_Fagus sylvatica_; 130-150 yrs), European silver fir (_Abies alba_; 130-150 yrs) and European ash (_Fraxinus excelsior_; 130-150 yrs) as dominant tree species. A second station is situated in the very vicinity outside of the forest (field station, SCF). The meteorological time series are presented in hourly time resolution of air temperature, relative humidity, precipitation, photosynthetically active radiation (PAR) and wind speed. Schänis is part of the Long-term Forest Ecosystem Research Programme (LWF) established and maintained by the Swiss Federal Research Institute WSL.

climate-change

display_name: climate-change

long-term

display_name: long-term

meteo

display_name: meteo

meteorology

display_name: meteorology

par

display_name: par

precipitation

display_name: precipitation

relative-humidity

display_name: relative-humidity

swiss-forest-lab

display_name: swiss-forest-lab

temperature

display_name: temperature

wind-speed

display_name: wind-speed

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Schänis, Switzerland: Long-term forest meteorological data from the Long-term Forest Ecosystem Research Programme (LWF), from 1998 onwards


maintainer: Yohann Chauvier

en: species-distribution synonyms filtering r-package occurrence-records landsca large-scale range-maps accepted-name ecoregions hard-limit flags conservation species-observations taxonomy gbif

en: GitHub URL

en: gbif.range R package version 0.2 (deprecated)

en: gbif.range R package version 1.1 (current)

en: wsl.gbif R package version 0.1 (deprecated)

en: Test version, not to be used.

en: Although species range may be obtained using expert maps or modeling methods, expert data is often species-limited and statistical models need more technical expertise as well as many species observations. When unavailable, such information may be extracted from the Global Biodiversity Information facility (GBIF), the largest public data repository inventorying georeferenced species observations worldwide. However, retrieving GBIF records at large scale may be tedious if users are unaware of specific tools and functions that need to be employed. Here we present *gbif.range*, an R library that contains automated methods to generate species range maps from scratch using in-house ecoregions shapefiles and an easy-to-use GBIF download wrapper. Finally, this library also offers a set of additional very useful parameters and functions for large GBIF datasets (generate doi, extract GBIF taxonomy, records filtering...). [gbif.range R project](https://github.com/8Ginette8/gbif.range)

accepted-name

display_name: accepted-name

conservation

display_name: conservation

ecoregions

display_name: ecoregions

filtering

display_name: filtering

flags

display_name: flags

gbif

display_name: gbif

hard-limit

display_name: hard-limit

landsca

display_name: landsca

large-scale

display_name: large-scale

occurrence-records

display_name: occurrence-records

r-package

display_name: r-package

range-maps

display_name: range-maps

species-distribution

display_name: species-distribution

species-observations

display_name: species-observations

synonyms

display_name: synonyms

taxonomy

display_name: taxonomy

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: gbif.range - An R package to generate species range maps based on ecoregions and a user-friendly GBIF wrapper


maintainer: Anita C Risch

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

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