maintainer: Anita C. Risch

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

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

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

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

format: XLS

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

Paper Citation:

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

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

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

ammonification

display_name: ammonification

climate

display_name: climate

fertilization

display_name: fertilization

global-change

display_name: global-change

grassland

display_name: grassland

herbivore

display_name: herbivore

mineralization

display_name: mineralization

nitrification

display_name: nitrification

nitrogen

display_name: nitrogen

nutrient-network

display_name: nutrient-network

soil

display_name: soil

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

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

en: Anthropogenic change and soil net N mineralization


maintainer: Amt für Abfall, Wasser, Energie und Luft des Kantons Zürich, Abteilung Luft, Klima und Strahlung - Sektion Monitoring

de: lufttemperaturen luft

label: LoRa-Datenbeschrieb

state: active

type: dataset

de: LoRa-Sensor-Messwerte 2021/02

media_type: text/csv

format: CSV

mimetype: text/csv

de: LoRa-Sensor-Messwerte 2020/04

de: LoRa-Sensor-Messwerte 2022/08

de: LoRa-Sensor-Messwerte 2022/11

de: LoRa-Sensor-Messwerte 2021/04

de: LoRa-Sensor-Messwerte 2019/10

de: LoRa-Sensor-Messwerte 2020/05

de: LoRa-Sensor-Messwerte 2021/11

de: LoRa-Sensor-Messwerte 2020/10

de: LoRa-Sensor-Messwerte 2022/09

de: LoRa-Sensor-Messwerte 2022/10

de: LoRa-Sensor-Messwerte 2022/04

de: LoRa-Sensor-Messwerte 2020/07

de: LoRa-Sensor-Messwerte 2022/12

de:
Seit Frühjahr 2019 betreibt das AWEL des Kantons Zürich ein Messnetz aus rund 40 Lufttemperatur und -feuchte Sensoren (LoRa). Die Sensor-Standorte befinden sich hauptsächlich im städtischen Raum, verdichtet in der Stadt Zürich, aber auch in weiteren Gemeinden im Kanton Zürich. Ziel des Messnetzes ist es, Messdaten für Stadtklima-Analysen bereitzustellen, damit Aspekte des Lokalklimas und Klimawandels stärker sichtbar werden und in der Raumplanung und Stadtentwicklung berücksichtigt werden können. Mit Hilfe dieser Daten lassen sich u.a. die Auswirkungen von Hitzewellen im städtischen Raum beschreiben. Das Messnetz wird durch das AWEL betrieben, unter Beteiligung der Stadt Zürich; die weiteren Standort-Gemeinden sind über die Messaktivität informiert.
Der hier verfügbare Datensatz wird täglich aktualisiert und umfasst alle, für die Standorte gültigen Messwerte im 10 Minuten Mittel, inklusive entsprechender Standortinformationen. Weitere Angaben zu Datensatz-Struktur und verwendeten Sensoren sind im Datensatz-Beschrieb zu finden.

awel-kanton-zuerich

de: Amt für Abfall, Wasser, Energie und Luft Kanton Zürich

de: Das Amt für Abfall, Wasser, Energie und Luft (AWEL) des Kantons Zürich hat den Auftrag, der Bevölkerung im Kanton Zürich eine gesunde, sichere und attraktive Umwelt zu bieten. Zu den Themen Abfall, Wasser, Energie und Luft stellt das AWEL Fachleuten und Interessierten Informationen zur Verfügung.

de: Lufttemperatur und Luftfeuchte LoRa-Sensor-Messwerte


maintainer: EnviDat Support

en: picea-abies forest nitrogen co2 face soil

en: Gross N mineralisation and nitrification

en: Gross rates of N mineralisation and nitrification (15N dilution)

en: Collection of soil solution

en: Soil solution collected by suction cups

en: Water analyses

en: Analyses of water, especially soil solution

en: Soil C and N

en: Analyses of C, N and their stable isotopes in the soil, and soil pH

en: Collection points

en: Collection points under the trees

en: Resin bags

en: Ammonium, nitrate and their 15N in resin bags

en: Data obtained in the free-air CO2 enrichment (FACE) experiment at Hofstetten, NW Switzerland, between 2009 and 2016. This dataset contains analyses of the soil solution throughout the experiment, especially for nitrate, as well as different analyses done at the end of the experiment: ammonium and nitrate captured by ion-exchange resin bags and extracted from soil cores, gross N mineralisation and nitrification measured by isotope dilution.

face

display_name: face

forest

display_name: forest

picea-abies

display_name: picea-abies

education

fr: Education, science

de: Bildung, Wissenschaft

it: Formazione e scienza

envidat

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

it: EnviDat: il portale dei dati ambientali

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

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

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

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester

en: Nitrogen availability under trees exposed to CO2 enrichment (FACE)


en: ice-layers snow-stratigraphy snow-temperature water-front weissfluhjoch snowmicropen

en: README

media_type: text/plain

format: TXT

mimetype: text/plain

en: WFJ_ICE_LAYERS dataset

media_type: application/zip

format: ZIP

mimetype: application/zip

en: The WFJ_ICE_LAYERS dataset contains multi-instrument snowpack measurements at high temporal resolution, which enable to monitor the formation of deep ice layers due to preferential water flow, at the Weissfluhjoch research site, Davos, Switzerland. It covers the winter 2016/2017, with a focus on the early melting season. This dataset includes traditional snowpack profiles (weekly resolution, 15/11/2016-29/05/2017), SnowMicroPen penetration resistance profiles (daily resolution, 01/02/2017-19/04/2017), snow temperatures measured at different heights in the snowpack (half-hourly resolution, 01/03/2017-15/04/2017) and the water front height derived from an upward-looking ground penetrating radar (3-hour resolution, 04/03/2017-08/04/2017). The measurements are complemented by initialization files for SNOWPACK model simulations with the ice reservoir parameterization at Weissfluhjoch for the winter 2016/2017.

ice-layers

display_name: ice-layers

snow-stratigraphy

display_name: snow-stratigraphy

snow-temperature

display_name: snow-temperature

snowmicropen

display_name: snowmicropen

water-front

display_name: water-front

weissfluhjoch

display_name: weissfluhjoch

education

en: Education and science

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: WFJ_ICE_LAYERS: Multi-instrument data for monitoring deep ice layer formation in an alpine snowpack


en: qgis-plugin swiss-geodata wms-data digital-planning digital-elevation-model geoinformation-system seilaplan open-source cable-yarding

en: View video on youtube!

format: SERVICE

en: WMS-Layer als Hintergrundkarten laden

en:
In order to digitally plan a cable line using the QGIS plugin ‘Seilaplan’, maps with various background information are helpful. In this tutorial we show you how to obtain maps that are helpful for cable line planning, for example a national map of Switzerland at different scales, the NFI vegetation height model or the NFI forest mix rate. For this we explain what WMS datasets are and how to integrate them into QGIS. No download of large data is needed for this, only a good internet connection. Please note that the tutorial language is German!

Link for the integration of WMS data: https://wms.geo.admin.ch/

Link to the description on the Swisstopo website:
https://www.geo.admin.ch/en/geo-services/geo-services/portrayal-services-web-mapping/web-map-services-wms.html

Link to the Seilaplan website: https://seilaplan.wsl.ch

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

Für die Verwendung des QGIS Plugins Seilaplan zur digitalen Seillinienplanung sind verschiedene Hintergrundkarten hilfreich. In diesem Tutorialvideo zeigen wir, was WMS Daten sind und wie man diese in QGIS einbinden kann. Dafür müssen die Daten nicht heruntergeladen werden. Es braucht lediglich eine gute Internetverbindung. Für die Seillinienplanung hilfreiche Karten sind bspw. die Landeskarte der Schweiz in verschiedenen Massstäben, das Vegetationshöhenmodell LFI oder der Waldmischungsgrad LFI.

Link zur Einbindung der WMS Daten: https://wms.geo.admin.ch/

Link zur Beschreibung auf der Swisstopo Webseite: https://www.geo.admin.ch/de/geo-dienstleistungen/geodienste/darstellungsdienste-webmapping-webgis-anwendungen/web-map-services-wms.html

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

cable-yarding

display_name: cable-yarding

digital-elevation-model

display_name: digital-elevation-model

digital-planning

display_name: digital-planning

geoinformation-system

display_name: geoinformation-system

open-source

display_name: open-source

qgis-plugin

display_name: qgis-plugin

seilaplan

display_name: seilaplan

swiss-geodata

display_name: swiss-geodata

wms-data

display_name: wms-data

en: Seilaplan Tutorial: Load WMS layers as background maps


maintainer: Alexandre Buttler

en: silvopastoral-systems biogeochemisty soil climate-change transplantation-experiment

en: Dataset

en: Zip archive of all datasets detailed in study.

en:
Silvopastoral systems are highly productive and combine long-term wood production with annual provision of forage for livestock. In the Swiss Jura Mountains these systems are a key component of the landscape. As in other cold biomes, climate change can potentially accelerate landscape change within these historically sustainable systems.

In order to anticipate the evolution of subalpine wooded pasture ecosystems under future climate and land-use changes, this project focused on the interplay between soil, vegetation and climate. It was aimed at providing experimental evidence for chief ecosystem processes, with emphasis on the quality of the ecosystem services provided. The main interest was placed on vegetation turf resistance to climate change along an unwooded – sparsely wooded - densely wooded pasture gradient (land-use intensity), where plant productivity, diversity and succession along with rates of carbon cycling and microbial activity provided measures of ecosystem functioning at both plot and landscape level.

Experimental transplantation of monolith soil turfs to lower altitudes allowed to simulate soil warming and reduced annual precipitation. In order to simulate a year-round warmer and drier climate the natural climate variation along an altitudinal gradient was used as a proxy. The aim was to simulate realistic climate change scenarios for the second half of the 21st century predicted by the IPCC report and downscaled for Switzerland providing regionalized interpolated projections integrating therein trends for temperature increase and precipitation decrease. By using permanent meteorological stations within the network of the Federal Office of Meteorology and Climatology (MeteoSwiss), we obtained high resolution regional data on the variation of mean annual temperature (MAT) and mean annual precipitation (MAP) in relation to altitude in the Swiss Jura Mountains. We observed a general increase of +0.5 K in MAT and a decrease of -20 % MAP for each 100 m decrease in altitude along the SE slope of the Swiss Jura Mountains. These relationships served for the selection of the transplantation sites such that in comparison to a control site at 1350 m a.s.l. (Combe des Amburnex, N 46°54’, E 6°23’) a +2 K MAT and -20 % MAP was achieved at 1010 m a.s.l. (Saint-George, N 46°52’, E 6°26’), a +4 K MAT and -40 % MAP at 570 m a.s.l., (Arboretum d’Aubonne, N 46°51’, E 6°37’), and a +5 K MAT and -50 % MAP at 395 m a.s.l. (Les Bois Chamblard, N 46°47’, E 6°41’). The two stations at 1010 m a.s.l. and 570 m a.s.l. corresponded to the IPCC scenario A1B for a moderate increase in greenhouse gas emissions and to scenario A2 for a high increase in greenhouse gas emissions, respectively. The station at 395 m a.s.l. was chosen to represent an extreme scenario with climate variables lying at the positive tail distribution of model predictions under the A2 scenario.

Soil cores were assembled into rectangular PVC boxes of 60  80 cm2 size and of 35 cm height. All mesocosms were dug down to surface level into previously prepared trenches in the ground thus preventing lateral heat exchange with the atmosphere. Since at each site the mesocosms were placed in a common garden with no light interception, mesocosms with turfs from the two wooded pastures were shaded from direct sun light to simulate the natural light conditions in the corresponding habitats. Each mesocosm was equipped with a drainage system and was connected to a water tank thus representing a zero potential lysimeter collecting soil solution and precipitation/snowmelt runoff. ECH2O EC-TM sensor probes coupled to Em50 data-loggers (Decagon Devices, Inc., USA) recorded soil temperature and volumetric water content in each mesocosm at the top-soil (0 to -3 cm) every minute and data were averaged over one hour intervals. Climate parameters at each transplantation site were monitored continuously throughout the experiment by means of automated weather stations (Sensor Scope Sàrl, Switzerland), measuring rain precipitation (non-heated tipping bucket gauges) and air temperature and humidity 2 m above the ground surface at one minute intervals.

A list of above- and belowground variables were measured to assess the resilience of biogeochemical processes, plant productivity, tree regeneration, and carbon sequestration for each respective land-use practice. Furthermore, the experimental data were used to improve on (parameterization) the existing spatially explicit, dynamic model WoodPaM and refine the modelʼs climatic and land-use variables so that different scenarios of climate change and land use change could be simulated. Natural and management induced disturbance patterns were incorporated into the model.
The data have been made available within the project CCES Mounted. The climate stations Sensorscope are still in use within the project CLIMARBRE (Wald und Klimawandel, WSL/BAFU).

#References

1. Puissant, J., Cécillon, L., Mills, R.T.E., Robroek, B.J.M. Gavazov, K., De Danieli, S., Spiegelberger, T., Buttler, A., Brun, J.J. 2015. Seasonal influence of climate manipulation on microbial community structure and function in mountain soils. Soil Biology and Biochemistry 80: 296–305.
2. Mills, R., K. Gavazov, T. Spiegelberger, D. Johnson and A. Buttler 2014. Diminished soil functions occur under simulated climate change in a sup-alpine pasture, but heterotrophic temperature sensitivity indicates microbial resilience. Science of the Total Environment, vol. 473–474(0): 465-472.
3. Gavazov, K., Spiegelberger, T. and Buttler, A. 2014. Transplantation of subalpine wood-pasture turfs along a natural climatic gradient reveals lower resistance of unwooded pastures to climate change compared to wooded ones. Oecologia (174) : 1425-1435.
4. Peringer A., Siehoff S., Chételat J., Spiegelberger T., Buttler A. & Gillet F. 2013. Past and future landscape dynamics in pasture-woodlands of the Swiss Jura Mountains under climate change. Ecology and Society, 18, 3: 11. DOI: 10.5751/ES-05600-180311. [online] URL: http://www.ecologyandsociety.org/vol18/iss3/art11/
5. Gavazov, K. S., A. Peringer, A. Buttler, F. Gillet and T. Spiegelberger. 2013. Dynamics of Forage Production in Pasture-woodlands of the Swiss Jura Mountains under Projected Climate Change Scenarios. Ecology and Society 18 (1): 38. [online] URL: http://www.ecologyandsociety.org/vol18/iss1/art38/

biogeochemisty

display_name: biogeochemisty

climate-change

display_name: climate-change

silvopastoral-systems

display_name: silvopastoral-systems

soil

display_name: soil

transplantation-experiment

display_name: transplantation-experiment

en: Biogeochemical data from a transplantation experiment of monolith soil turfs along an altitudinal gradient to simulate climate change scenarios


en: reflection snow photovoltaic irradiance mountain energy albedo

en: PV Data

en: Windows users use the 7-zip app to extract the zip-file.

en:
 Overview
The SUNWELL Modelling Environment is a combination of data and code that models electricity production from satellite-derived irradiance data and other spatial data sets for all of Switzerland. This ensemble accompanies the publication "The bright side of PV production in snow-covered mountains", published in the Proceedings of the National Academy of Science and reproduces all results and figures of. Code and resources are in their original form (with documentation). A new version with a more generalized application to PV modelling and with more flexibility in terms of input and output formats will be released in the coming months.

 Format
All code is written and has to be executed in Matlab. The input and output data sets are also in the Matlab-specific .mat format. Whenever publicly available, the original data is provided as geotif, .xlsx or other common format.
This is the case for:

- Digital Elevation Model (InputsFromMatlab/MSG/OriginalData/ASTERDEM),
- Landsurface cover type (InputsFromMatlab/MSG/OriginalData/CORINE),
- Population Density (InputsFromMatlab/MSG/OriginalData/popdensRaster,
- Electricity production from three of our validation sites (/Validation/WSL),
- Measured irradiance for two validation sites (/Validation/ASRB)
The ‘Metadata’ documents in the respective folders provide further information about the data sources and processing. Figures are produced either in .pdf or .png format.

 Structure
The central level of the SUNWELL environment holds the 5 Mains, which run the different modelling aspects of the paper; each code is documented separately. Additional code is located in the ‘DataProcessing’ and the ‘functions’ folder. Functions are called in the different Mains.

‘InputsFromMatlab’ contains the radiation and albedo input data sets in separate subfolders (SIS/SISDIR/ALB). The original data is not publicly available, but can be requested for research purposes free of charge. We provide a processed subset of the data set that was used to run the SUNWELL simulations. The MSG subfolder contains additional spatial input data sets.

‘Outputs’ contains the output files from the different mains (matching names, Main_CHallpixels.m  Prod_CHallpixels)

‘Publication_figures’ contains all individual figures from the PNAS publication, as well as the generating code (/code_plot) and the power point figures (/ppts) that provide the combined final figures.

‘Validation’ contains the data sets used in the model validation:

- Electricity production from three of our validation sites (/WSL),
- Measured irradiance for two validation sites (/ASRB)

Electricity production from a validation site at Lac des Toules in Wallis (/LDT), this data set was provided under an NDA and cannot be made publicly available.

Paper Citation:
 > _Annelen Kahl; Jérôme Dujardin; Michael Lehning (2018). Dataset on PV Production in Snow Covered Mountains. PNAS - Proceedings of the National Academy of Sciences. (in press)_

albedo

display_name: albedo

energy

display_name: energy

irradiance

display_name: irradiance

mountain

display_name: mountain

photovoltaic

display_name: photovoltaic

reflection

display_name: reflection

snow

display_name: snow


en: stillberg afforestation davos mountain-forest treeline-ecosystems foresty

en: Site photo

media_type: JPEG

en: Dataset V1

media_type: GDB

en: This dataset has been depricated in favour of an updated version.

format: GDB

en: Documentation

media_type: application/msword

mimetype: application/msword

en:
#Background Information

A history of many centuries of seasonal livestock grazing of alpine pastures has significantly depressed the treeline in many European mountain ranges. Therefore, observations of treeline patterns are more likely to be obscured by historical land use than in regions with less human influence and might not provide relevant information about environmental factors that influence mortality and growth patterns at the natural high-elevation limit of tree existence. To overcome this, we used a large afforestation experiment in the Swiss Alps, "Stillberg", to address the following questions:

Can trees survive above the current treeline in the European Alps? Are environmental factors that influence spatial patterns of mortality the same as those that influence height growth? Does the relative importance of these environmental variables change over the first 40 years after planting as seedlings?

#Dataset
The 5 ha Stillberg long-term afforestation research area is located near Davos (Central Alps, Switzerland). The bottom of the afforestation is approximately at the current treeline and it covers an altitudinal gradient of more than 150 m (2075 to 2230 m a.s.l.; see picture above). Approximately 92'000 seedlings of three treeline species (Larix decidua, Pinus mugo ssp. uncinata and Pinus cembra) were systematically planted in 1975, and mortality and height growth were closely monitored during the following 40 years.

Publications
Barbeito, I., Dawes, M. A., Rixen, C., Senn, J. and Bebi, P. (2012), Factors driving mortality and growth at treeline: a 30-year experiment of 92 000 conifers. Ecology, 93: 389–401. http://dx.doi.org/10.1890/11-0384.1

afforestation

display_name: afforestation

davos

display_name: davos

foresty

display_name: foresty

mountain-forest

display_name: mountain-forest

stillberg

display_name: stillberg

treeline-ecosystems

display_name: treeline-ecosystems

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

en: Long-term treeline research dataset at Stillberg, Davos


maintainer: Lisa Hülsmann

en: tree-growth forest-reserves inventory-data tree-mortality empirical-mortality-models

en: Calibration data

en:
A set of three consecutive inventories was used to generate records for the calibration of mortality models based on trees that were alive in the first and second inventory and either dead or alive in the third inventory.

 Meta data
Reserve: Name of the forest reserve

Plot: Permanent plot number

 Tree data
Species: Latin species name (or species group, cf. _Quercus_ spp.)

y: Status of the tree with 0 = alive and 1 = dead

dt: Length of the mortality period (years)

DBH: Diameter at breast height (mm) of the second inventory

relBAI: Annual relative basal area increment (-) calculated based on the first and the second DBH measurement as the compound annual growth rate of the trees basal area

 Stand data
_All stand variables were calculated for the second inventory._

mDBH: Mean diameter at breast height (mm)

iqrDBH: Interquartile range of diameter at breast height (mm)

BA: Basal area per hectare (m2 ha-1)

N: Number of trees per hectare (ha-1)

 Climate data
_Both climate variables were calculated for the period between the second and the third inventory._

P: Mean annual precipitation sum (mm)

mT: Mean annual air temperature (°C)

en:
The dataset comprises > 90 000 records from inventories in 54 strict forest reserves in [Switzerland](https://www.wsl.ch/de/wald/biodiversitaet-naturschutz-urwald/naturwaldreservate.html) and [Lower Saxony / Germany](http://naturwaelder.de/) along a considerable environmental gradient. It was used to develop parsimonious, species-specific mortality models for 18 European tree species based on tree size and growth as well as additional covariates on stand structure and climate.

 Inventory data
Measurements had been conducted repeatedly on up to 14 permanent plots per reserve for up to 60 years with re-measurement intervals of 4 - 27 years. The permanent plots vary in size between 0.03 and 3.47 ha. The inventories provide diameter measurements at breast height (DBH) and information on the species and status (alive or dead) of trees with DBH ≥ 4 cm for Switzerland and ≥ 7 cm for Germany.

 Data selection
We excluded three permanent plots where at least 80 % of the trees died during an interval of 10 years, and mortality could be clearly assigned to a disturbance agent. Mortality in the remaining stands was rather low, with a mean annual mortality rate of 1.5 % and strong variation between plots from 0 to 6.5 % (assessed for trees of all species with DBH ≥ 7 cm).
We only used data from permanent plots with at least 20 trees per species to obtain reliable plot-level mortality rates even for species with low mortality rates (about 5 % during 10 years), and selected tree species occurring on at least 10 plots to cover sufficient ecological gradients. This led to a dataset of 197 permanent plots and 18 tree or shrub species: _Abies alba_ Mill., _Acer campestre_ L., _Acer pseudoplatanus_ L., _Alnus incana_ Moench., _Betula pendula_ Roth, _Carpinus betulus_ L., _Cornus mas_ L., _Corylus avellana_ L., _Fagus sylvatica_ L., _Fraxinus excelsior_ L., _Picea abies_ (L.) Karst, _Pinus mugo_ Turra, _Pinus sylvestris_ L., _Quercus pubescens_ Willd., _Quercus_ spp. (_Q. petraea_ Liebl. and _Q. robur_ L.; not properly differentiated in the Swiss inventories), _Sorbus aria_ Crantz, _Tilia cordata_ Mill. and _Ulmus glabra_ Huds..

 Predictors of tree mortality
We considered tree size and growth as key indicators for mortality risk. Radial stem growth between the first and second inventory and DBH at the second inventory were used to predict tree status (alive or dead) at the third inventory. To this end, the annual relative basal area increment (relBAI) was calculated as the compound annual growth rate of the trees basal area. Additional covariates on stand structure and climate comprise mean annual precipitation sum (P), mean annual air temperature (mT), the mean and the interquartile range of DBH (mDBH, iqrDBH), basal area (BA) and the number of trees (N) per hectare.

 Further information
For further information, refer to Hülsmann _et al_. (in press) How to kill a tree – Empirical mortality models for eighteen species and their performance in a dynamic forest model. _Ecological Applications_.

empirical-mortality-models

display_name: empirical-mortality-models

forest-reserves

display_name: forest-reserves

inventory-data

display_name: inventory-data

tree-growth

display_name: tree-growth

tree-mortality

display_name: tree-mortality

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: Calibration data for empirical mortality models of 18 European tree species


maintainer: James W. Kirchner

en: groundwater sap-flow snowmelt diel-cycle evapotranspiration diurnal-cycle ecohydrology streamflow

en: Subcatchment polygons

en: Subcatchment polygons for use in the Google Earth Engine scripts.

en: Metadata

media_type: application/pdf

en: Column descriptions for both data files, and map coordinates for all measurement locations.

format: PDF

mimetype: application/pdf

en: MODIS NDSI script

en: Google Earth Engine script for computing MODIS normalized difference snow index (NDSI) values and averaging them over subcatchment areas.

en: 30-minute time series

en: 30-minute hydrometeorological and ecohydrological time series.

en: Daily time series

en: Daily hydrometeorological and ecohydrological time series

en: MODIS EVI2 script

en: Google Earth Engine script for computing MODIS two-band enhanced vegetation index (EVI2) values and averaging them over subcatchments.

en:
Hydrometerological and ecohydrological time series from Sagehen Creek and Independence Creek, Sierra Nevada, USA, illustrating hydrological responses to daily cycles in snowmelt and evapotranspiration forcing.

Data include 30-minute time series of

- weather variables,

- sap flow fluxes,

- groundwater levels (in two riparian transects of shallow groundwater wells),

- and stream stages (at 12 sites spanning a 500-meter elevation gradient),

and daily time series of

- temperature, precipitation, and snow water equivalent at three nearby snow telemetry stations

- diel cycle index values for groundwater levels and stream stages,

- and MODIS normalized difference snow index (NDSI) and enhanced vegetation index (EVI2) values averaged over selected subcatchments.

Google Earth Engine scripts for extracting the MODIS data are also provided.

diel-cycle

display_name: diel-cycle

diurnal-cycle

display_name: diurnal-cycle

ecohydrology

display_name: ecohydrology

evapotranspiration

display_name: evapotranspiration

groundwater

display_name: groundwater

sap-flow

display_name: sap-flow

snowmelt

display_name: snowmelt

streamflow

display_name: streamflow

en: Daily cycles in solar flux, snowmelt, transpiration, groundwater, and streamflow at Sagehen and Independence Creeks, Sierra Nevada, USA


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

en: CORE Dataset (80 GB)

media_type: CORE

media_type: .sh

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

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

media_type: HTML

format: HTML

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

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

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

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

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

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

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

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

CORE Specifications:

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

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

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

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

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

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

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

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

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

Major issues to be solved:

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

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

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

Example

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

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

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

core

display_name: core

geodata

display_name: geodata

gridded-data

display_name: gridded-data

raster

display_name: raster

research-data-management

display_name: research-data-management

time-series

display_name: time-series

video-encoding

display_name: video-encoding

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

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


maintainer: Marcia Phillips

en: temperature gemsstock permos permafrost borehole

en: DATASET & PLOTS

media_type: URL

en: Data visualiser and download service provided by PERMOS.

en: PERMOS DATA POLICY

media_type: PDF

en: Map of PERMOS borehole network

en:
Processed ground temperature measurements at the Gemsstock permafrost borehole in canton Uri, Switzerland. The borehole is located at 2940 m asl on a steep (50°) North-West slope (315°). The surface material is bedrock and borehole depth is 40 m. Thermistors used YSI 44008. Year of drilling 2006. This borehole is part of the Swiss Permafrost network, PERMOS (www.permos.ch). Contact phillips@slf.ch for details of processing applied.

Publications

1. A. Haberkorn, M. Phillips, R. Kenner, H. Rhyner, M. Bavay, S.P. Galos, M. Hoelzle. Thermal regime of rock and its relation to snow cover in steep Alpine rock walls: Gemsstock, central Swiss Alps. 2015. Geografiska Annaler: Series A, Physical Geography. Volume 97. Issue 3. 579–597. http://dx.doi.org/10.1111/geoa.12101. 10.1111/geoa.12101.

2. R. Kenner, M. Phillips, C. Danioth, C. Denier, P. Thee, A. Zgraggen. Investigation of rock and ice loss in a recently deglaciated mountain rock wall using terrestrial laser scanning: Gemsstock, Swiss Alps. 2011. Cold Regions Science and Technology. Volume 67. Issue 3. 157–164. http://dx.doi.org/10.1016/j.coldregions.2011.04.006. 10.1016/j.coldregions.2011.04.006.

borehole

display_name: borehole

gemsstock

display_name: gemsstock

permafrost

display_name: permafrost

permos

display_name: permos

temperature

display_name: temperature

en: Processed permafrost borehole data (2940 m asl), Gemsstock, Switzerland


maintainer: Dirk Nikolaus Karger

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

en: CHELSA V1.2 (deprecated)

media_type: geotiff

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

format: GeoTIFF

en: CHELSA V2.1 (current)

en: CHELSA Version 2.1 based on ERA5 forcings and CMIP6

en: CHELSA Technical Documentation

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

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

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

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

bioclim

display_name: bioclim

climate

display_name: climate

climate-change

display_name: climate-change

climate-layers

display_name: climate-layers

climatology

display_name: climatology

free-climate-data

display_name: free-climate-data

precipitation

display_name: precipitation


envidat-supports-open-science

maintainer: Ionut Iosifescu Enescu

title_for_slug: envidat-supports-open-science

en: juypter-notebooks envidat open-science research-data-management wsl

en: Article as PDF

en:
The article "EnviDat Supports Open Science" originally appeared in WSLintern No. 3 (2020), page 14-15 and it is republished here with permission from the WSLintern editorial team. It contains guidelines for WSL scientists about the main issues behind Open Science and how to pragmatically approach the complexities of doing Open Science with EnviDat’s support.

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

envidat

display_name: envidat

juypter-notebooks

display_name: juypter-notebooks

open-science

display_name: open-science

wsl

display_name: wsl

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

en: EnviDat Supports Open Science


maintainer: Christoph Marty

en: Dataset

en: Datasets are grouped under experiment "Aare" and "Graubunden".Data is also made available in netCDF format. A range of tools are available to [read netCDF](http://www.unidata.ucar.edu/software/netcdf/software.html).

en: Graubunden input

media_type: application/x-tar

en: Grid input data for the 'Graubunden' CH2014 simulation.

mimetype: application/x-tar

en: Aare input

en: Grid input data for the 'Aare' CH2014 simulation.

en: CH2014 REPORT

en:
CH2014-Impacts: Toward Quantitative Scenarios of Climate Change Impacts in Switzerland

The CH2014-Impacts report presents the results of a concerted effort of 15 Swiss research institutions to describe impacts of climate change in Switzerland quantitatively. This dataset contributes to the Crosphere/Snowcover section.

en:
 Overview

The CH2014-Impacts initiative is a concerted national effort to describe impacts of climate change in Switzerland quantitatively, drawing on the scientific resources available in Switzerland today. The initiative links the recently developed Swiss Climate Change Scenarios CH2011 with an evolving base of quantitative impact models. The use of a common climate data set across disciplines and research groups sets a high standard of consistency and comparability of results. Impact studies explore the wide range of climatic changes in temperature and precipitation projected in CH2011 for the 21st century, which vary with the assumed global level of greenhouse gases, the time horizon, the underlying climate model, and the geographical region within Switzerland. The differences among climate projections are considered using three greenhouse gas scenarios, three future time periods in the 21st century, and three climate uncertainty levels (Figure 1). Impacts are shown with respect to the reference period 1980-2009 of CH2011, and add to any impacts that have already emerged as a result of earlier climate change.

 Experimental Setup

Future snow cover changes are simulated with the physics-based model Alpine3D (Lehning et al., 2006). It is applied to two regions: The canton of Graubünden and the Aare catchment. These domains are modeled with a Digital Elevation Model (DEM) with a resolution of 200 m × 200 m. This defines the simulation grid that has to be filled with land cover data and downscaled meteorological input data for each cell for the time period of interest at hourly resolution. The reference data set consists of automatic weather station data. All meteorological input parameters are spatially interpolated to the simulation grid. The reference period comprises only thirteen years (1999–2012), because the number of available high elevation weather stations for earlier times is not sufficient to achieve unbiased distribution of the observations with elevation. The model uses projected temperature and precipitation changes for all greenhouse gas scenarios (A1B, A2, and RCP3PD) and CH2011 time periods (2035, 2060, and 2085).

 Data

Snow cover changes are projected to be relatively small in the near term (2035) (Figure 5.1 top), in particular at higher elevations above 2000 m asl. As shown by Bavay et al. (2013) the spread in projected snow cover for this period is greater between different climate model chains (Chapter 3) than between the reference period and the model chain exhibiting the most moderate change. In the 2085 period much larger changes with the potential to fundamentally transform the snow dominated alpine area become apparent (Figure 5.1 bottom). These changes include a shortening of the snow season by 5–9 weeks for the A1B scenario. This is roughly equivalent to an elevation shift of 400–800 m. The slight increase of winter precipitation and therefore snow fall projected in the CH2011 scenarios (with high associated uncertainty) can no longer compensate for the effect of increasing winter temperatures even at high elevations. In terms of Snow Water Equivalents (SWE), the projected reduction is up to two thirds toward the end of the century (2085). A continuous snow cover will be restricted to a shorter time period and/or to regions at increasingly high elevation. In Bern, for example, the number of days per year with at least 5 cm snow depth will decrease by 90% from now 20 days to only 2 days on average.

snow-depth

display_name: snow-depth

snow-water-equivalent

display_name: snow-water-equivalent

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: Alpine3D simulations of future climate scenarios CH2014


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

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

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

Paper Citation:

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

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

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

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

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

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

bacteria

display_name: bacteria

biogeochemical-cycling

display_name: biogeochemical-cycling

carcass

display_name: carcass

carrion

display_name: carrion

decomposition

display_name: decomposition

detritus

display_name: detritus

fungi

display_name: fungi

soil-biodiversity

display_name: soil-biodiversity

wolf

display_name: wolf

yellowstone

display_name: yellowstone

en: Elk and bison carcasses in Yellowstone, USA