maintainer: Anita C. Risch

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

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

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

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

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

Paper Citation:

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

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

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

ammonification

display_name: ammonification

climate

display_name: climate

fertilization

display_name: fertilization

global-change

display_name: global-change

grassland

display_name: grassland

herbivore

display_name: herbivore

mineralization

display_name: mineralization

nitrification

display_name: nitrification

nitrogen

display_name: nitrogen

nutrient-network

display_name: nutrient-network

soil

display_name: soil

en: Anthropogenic change and soil net N mineralization


ceilometer-klosters

maintainer: EnviDat Support

title_for_slug: ceilometer-klosters

en: cloud-base-height klosters ceilometer raclets vertical-visibility

state: active

type: dataset

en: Ceilometer data Klosters

media_type: mat

en: Ceilometer installed at Klosters measuring cloud base height, cloud cover and backscatter profile. The time resolution is 30 seconds, with time given in number of days from January 0, 0000, UTC. More information on the variables in the manufacturer manual.

format: SERVICE

en: Cloud base height (m) and vertical visibility (m) were measured with the VAISALA Ceilometer CL31 in Klosters (LON: 9.880413, LAT: 46.869019). The CL31 is an instrument with constant reliability for all weather conditions and simultaneous detection of three cloud layers in heights up to 7.6 km.

ceilometer

display_name: ceilometer

cloud-base-height

display_name: cloud-base-height

klosters

display_name: klosters

raclets

display_name: raclets

vertical-visibility

display_name: vertical-visibility

education

fr: Education, science

de: Bildung, Wissenschaft

en: Education and science

it: Formazione e scienza

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


radiosondes

maintainer: Jörg Wieder

title_for_slug: radiosondes

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

en: Radiosond Template Files

media_type: application/zip

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

format: ZIP

mimetype: application/zip

en: Radiosond Sounding Files

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

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

holoballoon

display_name: holoballoon

pressure

display_name: pressure

radiosonde

display_name: radiosonde

relative-humidity

display_name: relative-humidity

temperature

display_name: temperature

wind

display_name: wind

wind-direction

display_name: wind-direction

en: Radiosondes


maintainer: Albrecht von Boetticher

en: interfoam debris-flow dberisintermixing murgang finite-volumes natural-hazard shallow-landslide hangmuren barrier-impact-simulation

en: debrisInterMixing for OF-Version 6

en:
This version of debrisInterMixing is based on OpenFOAM Version 4.x and was tested on Ubuntu 18.04.1 LTS.
Please follow the instructions given in ReadMe_debrisInterMixing4.pdf for installation.

en: Installation Instructions for debrisInterMixing-6

media_type: application/pdf

en: These are the installation instructions for debrisInterMixing4

format: PDF

mimetype: application/pdf

en:
Here the updated versions of debrisInterMixing are provided for download. The first OpenFoam-compatible Version 2.3.x are available as supplement to v. Boetticher, A., Turowski, J. M., McArdell,W. B., Rickenmann, D., Hürlimann, M., Scheidl, C., and Kirchner, J. W.: DebrisInterMixing-2.3: A Finite Volume solver for three dimensional debris flow simulations based on two calibration parameters. Part two: model validation with experiments. Geoscientific Model Development, 10, 11: 3963-3978. doi: 10.5194/gmd-10-3963-2017.
DebrisInterMixing is a Volume-of-Fluid based Finite Volume code that accounts for shear-thinning sensitive shares of fine sediment suspension together with pressure-sensitive components of the gravel grains within debris flow mixtures. All model properties can be derived from a material sample except for a grid-sensitive calibration parameter.
For more information, please contact albrecht.vonboetticher@wasserbau.ch.

barrier-impact-simulation

display_name: barrier-impact-simulation

dberisintermixing

display_name: dberisintermixing

debris-flow

display_name: debris-flow

finite-volumes

display_name: finite-volumes

hangmuren

display_name: hangmuren

interfoam

display_name: interfoam

murgang

display_name: murgang

natural-hazard

display_name: natural-hazard

shallow-landslide

display_name: shallow-landslide

en: Three-dimensional debris flow simulation tool debrisInterMixing


biogas-aus-hofdunger-in-der-schweiz

title_for_slug: biogas-aus-hofdunger-in-der-schweiz

en: bioenergy manure switzerland review biogas

en: Biogas aus Hofdünger in der Schweiz

media_type: PDF

en: Ziel dieses Whitepapers ist es, Entscheidungsträgern, Verwaltungen und Stakeholdern die aktuellsten Forschungsergebnisse zur Verfügung zu stellen, um die optimale Nutzung von Bioenergie aus Hofdünger in der Schweizer Energiewende zu fördern. Zu diesem Zweck werden die Ergebnisse des Schweizer Kompetenzzentrums für Bioenergieforschung - SCCER BIOSWEET - zusammengefasst und in einem breiteren Kontext dargestellt. Wenn nichts anderes erwähnt wird, beziehen sich die Ergebnisse auf die Schweiz und im Falle der Ressourcen auf die heimischen Biomassepotenziale.

bioenergy

display_name: bioenergy

biogas

display_name: biogas

manure

display_name: manure

review

display_name: review

switzerland

display_name: switzerland

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal


en: buffer gis geospatial jupyter density roads

en: Road Density Calculation (Python Notebook) online

media_type: url

en: Road data from OSM for the surroundings of Menzberg, Switzerland

en: Plots

en: Random plot locations

media_type: ipynb

en:
This dataset contains the road and plot data used for the geospatial analysis example showcased in "Fostering Open Science at WSL with the EnviDat Environmental Data Portal", a contribution to the 5th Open Source Geospatial Research and Education Symposium (OGRS), 2018.
The example uses Jupyter Notebook to calculate road densities in the neighbourhood of sample plot locations with Python. Road data were extracted from OpenStreetMap, while the point data (sample plots) were generated manually.

buffer

display_name: buffer

density

display_name: density

geospatial

display_name: geospatial

gis

display_name: gis

jupyter

display_name: jupyter

roads

display_name: roads

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: Dataset for OGRS 2018 publication


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

media_type: TIFF

en: Soil organic carbon content modified from SoilGrids

format: TIFF

en: Frost change frequency modified from CHELSA

en: Soil pH modified from SoilGrids

media_type: tif

en: Precipitation seasonality from CHELSA

en: Mean diurnal range modified from CHELSA

en: Precipitation of driest quarter modified from CHELSA

en: Soil clay content modified from SoilGrids

en: Mean annual temperature modified from CHELSA

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

aridity

display_name: aridity

frost-change-velocity

display_name: frost-change-velocity

ground-temperature

display_name: ground-temperature

precipitation

display_name: precipitation

sdms

display_name: sdms

soil

display_name: soil

en: Environmental layers for SDM simulations (GDPlants)


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

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

format: XLS

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

Paper Citation:

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

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

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

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

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

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

bacteria

display_name: bacteria

biogeochemical-cycling

display_name: biogeochemical-cycling

carcass

display_name: carcass

carrion

display_name: carrion

decomposition

display_name: decomposition

detritus

display_name: detritus

fungi

display_name: fungi

soil-biodiversity

display_name: soil-biodiversity

wolf

display_name: wolf

yellowstone

display_name: yellowstone

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

en: Elk and bison carcasses in Yellowstone, USA


en: picea-abies co2 nitrogen face soil forest

en: Resin bags

media_type: text/csv

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

format: CSV

mimetype: text/csv

en: Water analyses

en: Analyses of water, especially soil solution

en: Gross N mineralisation and nitrification

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

en: Collection points

en: Collection points under the trees

en: Soil C and N

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

en: Collection of soil solution

en: Soil solution collected by suction cups

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

nitrogen

display_name: nitrogen

picea-abies

display_name: picea-abies

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: Nitrogen availability under trees exposed to CO2 enrichment (FACE)


en: sense-of-place landscape place-attachment urban-ecology electrodermal-activity

en: 02. Questionnaire Data

media_type: ZIP

en: 01. EDA Processed Data

en: 04. Scripts

en: 00. EDA Raw Data

en: 05. Visuals

en: 03. Publication Ready Data

en: This repository contains data on EDA measurements of visitors with different cultural backgrounds in virtual urban park settings. The parks are a Persian garden (Shiraz, Iran) and a historical park in Zurich, Switzerland. The cultural background of the visitors is Persian and Central European. The repository contains raw data from EDA, processed time series and statistical procedures.

electrodermal-activity

display_name: electrodermal-activity

landscape

display_name: landscape

place-attachment

display_name: place-attachment

sense-of-place

display_name: sense-of-place

urban-ecology

display_name: urban-ecology

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

en: Electrodermal Activity (EDA) of Bi-cultural Visitors In Virtual Park Settings


maintainer: Leonie Schönbeck

en: vpd pinus-sylvestris photosynthesis gas-exchange acclimation fluorescence drought

en: Methods ACi curves

media_type: CSV

en:
Parameters and methods for CO2 response curves of Photosynthesis.

Metatdata format by
Ely et al (2021). A reporting format for leaf-level gas exchange data and metadata. Ecological Informatics. https://doi.org/10.1016/j.ecoinf.2021.101232

en: Methods fluorescence

en:
Parameters and methodology for Fluorescence measurements

Metadata format from
Ely et al (2021). A reporting format for leaf-level gas exchange data and metadata. Ecological Informatics. https://doi.org/10.1016/j.ecoinf.2021.101232

en: Methods light response

en: Methods and parameters for light response curves

en: Methods survey

media_type: .csv

en: Parameters and methods for point measurements / survey measurements

en: Methods diurnal measurements

en:
Parameters and methodology for diurnal measurements

Metadata format from
Ely et al (2021). A reporting format for leaf-level gas exchange data and metadata. Ecological Informatics. https://doi.org/10.1016/j.ecoinf.2021.101232

en: Variable names and units

media_type: XLSX

en: Variable names, description and units for the Pfynwald Gas exchange data

en: Instrument details

en:
Instrument details for all LiCor instruments used during collection of dataset.

Metatdata created following guidelines from
Ely et al (2021). A reporting format for leaf-level gas exchange data and metadata. Ecological Informatics. https://doi.org/10.1016/j.ecoinf.2021.101232

en: Pfynwald gas exchange

en:
Full dataset of Pfynwald gas exchange data taken in the period from June 2013 to October 2020.
All metdata can be found in attached documents.
Data uploaded after quality control using guidelines by Ely et al (2021). A reporting format for leaf-level gas exchange data and metadata. Ecological Informatics. https://doi.org/10.1016/j.ecoinf.2021.101232

en:
Gas exchange was measured on control, irrigated and irrigation-stop trees at the irrigation experiment Pfynwald, during the years 2013, 2014, 2016-2020. The measurement campaigns served different purposes, resulting in a large dataset containing survey data, CO2 response curves of photosynthesis, light response curves of photosynthesis, and fluorescence measurements. Measurements were done with LiCor 6400 and LiCor 6800 instruments. Until 2016, measurements were done on excised branches or branches lower in the canopy. From 2016 onwards, measurements were done in the top of the canopy using fixed installed scaffolds.
All metadata can be found in the attached documents.

acclimation

display_name: acclimation

drought

display_name: drought

fluorescence

display_name: fluorescence

gas-exchange

display_name: gas-exchange

photosynthesis

display_name: photosynthesis

pinus-sylvestris

display_name: pinus-sylvestris

vpd

display_name: vpd

en: 2013-2020 gas exchange at Pfynwald


maintainer: Rebecca Mott

en: patchy-snow-covers tls riegl-vz6000 snow-ablation-rates terrestrial-laser-scanning

en: Digital Elevation Modell 2015

en: Digital Elevation Modell 2014

en: Snowmask 2015

en: Snowmask 2016

en: Snowmask 2014

en: Gletschboden test site

en: Digital Elevation Modell 2016

en: We recorded snow ablation maps with a terrestrial laser scanner (TLS, Riegl-VZ6000) at the Gletschboden area. The TLS position is located approximately 30 vertical meters above the Gletschboden area at a northerly exposed slope. In total 44 TLS measurement sets have been conducted in three consecutive years 2014-2016 (2014: 13 measurements; 2015: 17 measurements; 2016: 14 measurements). The TLS system has a single-point measurement frequency of 300 kHz and a beam divergence of 0.007°. This set-up allows a horizontal resolution of approximately 0.01 m in 100 m distance to the TLS position. One scan of the Gletschboden area lasts approximately 15 minutes. The travel time from the laser scanner towards the surface is recorded and afterwards converted into a point cloud of distances. 5 reflectors located at the Gletschboden area and in the closer surroundings were additionally scanned during each measurement to transform the point cloud from the scanner own coordinate system into Swiss coordinates. Additionally, orthophotos have been created by using pictures recorded from the TLS in order to provide snow mask maps. Snow and bare ground can be distinguished by the RGB color information of the orthophoto. Cells with blue band information greater than 175 were categorized as snow and all cells with values smaller or equal 175 were categorized as bare ground.

patchy-snow-covers

display_name: patchy-snow-covers

snow-ablation-rates

display_name: snow-ablation-rates

terrestrial-laser-scanning

display_name: terrestrial-laser-scanning

tls

display_name: tls

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: DISCHMEX - High-resolution daily snow ablation rates in an Alpine environment


lidar-davos-wolfgang

maintainer: Patric Seifert

title_for_slug: lidar-davos-wolfgang

en: tropos polly lidar raclets clouds

en: Raman Lidar Davos

en: A portable Raman lidar system (Polly) from Leibnitz Institute for Tropospheric Research (Tropos) was deployed at Davos Wolfgang (LON: 9.853594, LAT: 46.835577). Please use this [link](http://polly.tropos.de/?p=lidarzeit&Ort=39), to be directly forwarded to the Davos location and select the date of interest from the calendar (bold numbers). The data can be requested directly at the Polly team.

clouds

display_name: clouds

lidar

display_name: lidar

polly

display_name: polly

raclets

display_name: raclets

tropos

display_name: tropos

en: LIDAR Davos Wolfgang


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

de: luft lufttemperaturen

label: LoRa-Datenbeschrieb

de: LoRa-Sensor-Messwerte 2023/02

de: LoRa-Sensor-Messwerte 2020/03

de: LoRa-Sensor-Messwerte 2023/03

de: LoRa-Sensor-Messwerte 2020/07

de: LoRa-Sensor-Messwerte 2021/10

de: LoRa-Sensor-Messwerte 2019/04

de: LoRa-Sensor-Messwerte 2023/06

de: LoRa-Sensor-Messwerte 2021/09

de: LoRa-Sensor-Messwerte 2020/05

de: LoRa-Sensor-Messwerte 2019/09

de: LoRa-Sensor-Messwerte 2019/05

de: LoRa-Sensor-Messwerte 2020/06

de: LoRa-Sensor-Messwerte 2022/09

de: LoRa-Sensor-Messwerte 2020/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: Matthias Haeni

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

en: Data plot of the field station (JUF)

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

en: Site aerial photo (Jussy)

media_type: image/jpeg

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

mimetype: image/jpeg

en: Data file.

en: Data plot of the forest station (JUB)

en: Forest station data (JUB)

en: Site photo of the field station (JUF)

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

en: Site photo of the forest station (JUB)

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

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 Jussy in Switzerland where one station is located within a natural broad-leaved forest stand (JUB) with sessile oak (_Quercus petrea_; 90 yrs), aspen (_Populus tremula_; 60 yrs) and European hornbeam (_Carpinus betulus_; 60 yrs) as dominant tree species. A second station is situated in the very vicinity outside of the forest (field station, JUF). The meteorological time series are presented in hourly time resolution of air temperature, relative humidity, precipitation, photosynthetically active radiation (PAR) and wind speed. Jussy 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

education

de: Bildung, Wissenschaft

envidat

de: EnviDat: das Umwelt-Datenportal

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.

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


en: groundwater switzerland alptal studibach time-series

media_type: text/plain

en:
Groundwater time series of the monitoring network in the Studibach (C7) catchment, Alptal, Switzerland published in
Rinderer M., van Meerveld I, McGlynn B. (2019): From points to patterns – Assessing runoff source area dynamics and hydrological connectivity using time series clustering. Water Resources Research, doi: 2018WR023886R

format: TXT

mimetype: text/plain

en: Groundwater time series between 2010 and 2014 of the distributed monitoring system in the Studibach (C7), Alptal, Switzerland. Data published in Rinderer M., van Meerveld I, McGlynn B. (2019): From points to patterns – Assessing runoff source area dynamics and hydrological connectivity using time series clustering. Water Resources Research, doi: 2018WR023886R

alptal

display_name: alptal

groundwater

display_name: groundwater

studibach

display_name: studibach

switzerland

display_name: switzerland

time-series

display_name: time-series

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: Groundwater time series Studibach (Rinderer et al., 2019, WRR)


disdrometer-data-gotschnagrat

title_for_slug: disdrometer-data-gotschnagrat

en: raclets visibility precipitation disdrometer hydrometeors

en: Disdrometer Data Gotschnagrat

en: The zip folder contains one netCDF file per day for Gotschnagrat from February 11 to March 27 2019 with a temporal resolution of 30 seconds (time in UTC). The naming convention for the files is Gotschna_YYYYMMDD.nc and the files include the following: characteristics of hydrometeors (e.g. diameter, fall velocity, amount per diameter class,...), precipitation rate, radar reflectivity, visibility range, weather codes and instrument performance.

en: A laser optical disdrometer (Parsivel² , OTT Hydromet) was deployed at Gotschnagrat (LON: 9.849, LAT: 46.859) to measure hydrometeors by extinction when passing a laser beam. The instrument can classify eight different kinds of precipitation, including rain, hail, snow, drizzle, and hybrid forms. The dataset contains information on precipitation amount and type for the period of February 11 to March 27 2019 at Gotschnagrat.

disdrometer

display_name: disdrometer

hydrometeors

display_name: hydrometeors

raclets

display_name: raclets

visibility

display_name: visibility


en: davos incoming outgoing longwave-radiation shortwave-radiation

en: Processing doc

en: Metadata

en: DATA ACCESS

media_type: SMET

en: Parameter file

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

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

davos

display_name: davos

incoming

display_name: incoming

longwave-radiation

display_name: longwave-radiation

outgoing

display_name: outgoing

shortwave-radiation

display_name: shortwave-radiation

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


greenland-shrubs-and-microclimate

maintainer: Nathalie Chardon

title_for_slug: greenland-shrubs-and-microclimate

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

en: Nuup Kangerlua

en:
 Collected July-August 2011-2013 (PI: Jacob Nabe-Nielsen)

These data collected July-August 2011-2013 (PI: Jacob Nabe-Nielsen). Please contact J. Nabe-Nielsen (jnn_at_ecos.au.dk) to request access to data on all surveyed shrub species. Available here are only data used in Chardon et al. 2022.

 COLUMN ABBREVIATIONS

 Metadata for 90 m spatial grid

year: sampling year

alt: average altitude of grid as recorded during field surveys [m]

long: WGS84 longitude

lat: WGS84 latitude

x_stere: longitude in Polar Stereographic (EPSG 4326)

y_stere: latitude in Polar Stereographic (EPSG 4326)

 Species presence [1] or absence [0] at grid level

pa_bena: _Betula nana_

pa_sagl: _Salix glauca_

 Shrub species cover (calculated from pin hits in 0.7 x 0.7 m frame, based on total of 25 hits) [%]

cov_bena: _Betula nana_ cover

cov_sagl: _Salix glauca_ cover

 Downscaled climate (time series 1984 – 2013)

tempjja: average maximum June-August temperature [ºC]

tempmax: yearly maximum temperature [ºC]

tempmin: yearly minimum temperature [ºC]

tempcont: yearly maximum – minimum temperature [ºC]

precipmam: cumulative March-May precipitation [mm]

precipjja: cumulative June-August precipitation [mm]

insoljja: average June-August insolation [MJ cm-2 yr-1]

 Topography and terrain wetness (calculated from)

twi: SAGA wetness index (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015)

slope: slope angle (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

asp: aspect (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

tcws: Tasseled Cap Wetness component (Sentinal images)

 Community characteristics

rich: shrub species richness averaged across pin-point frames per grid

shan: Shannon diversity (‘vegan’ package; Oksanen et al. 2019)

comp_bena: shrub competition on Betula nana (i.e. summed cover of non-_B. nana_ species) [% cover]

comp_sagl: shrub competition on Salix glauca (i.e. summed cover of non-_S. glauca_ species) [% cover]

en: GitLab Repository

media_type: .zip

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

en: Kangerluarsunnguaq

en:
 Collected August 2020 (PI: Nathalie Isabelle Chardon)

These data collected August 2020 (PI: Nathalie Isabelle Chardon). Please contact N. I. Chardon (nathalie.chardon_at_gmail.com) for individual quadrat (i.e. 1 m x 1 m) data, which are averaged at the 90 m spatial grid here.

 COLUMN ABBREVIATIONS

 Metadata for 90 m spatial grid

id: sampling zone and unique 90 m grid identifier

long: WGS84 longitude

lat: WGS84 latitude

x_stere: longitude in Polar Stereographic (EPSG 4326)

y_stere: latitude in Polar Stereographic (EPSG 4326)

date: sampling day.month.year

 Species presence [1] or absence [0] at grid level

alcr: _Alnus crispa_

bena: _Betula nana_

dila: _Diapensia lapponica_

emni: _Empetrum nigrum_

hahy: _Harimanella hypoides_

lopr: _Loiseleuria procumbens_

oxpa: _Oxycoccus palustre_

phca: _Phyllodoce caerulea_

rhod: _Rhododendron_ sp.

saar: _Salix arctica_

sagl: _Salix glauca_

sahe: _Salix herbacea_

vaul: _Vaccinium uliginosum_

drin: _Dryas integrifolia_

thpr: _Thymus praecox_

 Habitat characteristics [% cover]

shrub: all shrub species

gram: graminoids

forb: forbs

bryo: bryophytes

lich: lichens

litt: litter

bare: bare ground

stone: stones (> 2 cm)

org: organic crust

 Shrub species height (maximum vegetative height) [mm] and cover (estimated in 1x1 m quadrat) [%]

alcr_c: _Alnus crispa_ cover

alcr_h: _Alnus crispa_ height

bena_c: _Betula nana_ cover

bena_h: _Betula nana_ height

dila_c: _Diapensia lapponica_ cover

dila_h: _Diapensia lapponica_ height

emni_c: _Empetrum nigrum_ cover

emni_h: _Empetrum nigrum_ height

hahy_c: _Harimanella hypoides_ cover

hahy_h: _Harimanella hypoides_ height

lopr_c: _Loiseleuria procumbens_ cover

lopr_h: _Loiseleuria procumbens_ height

oxpa_c: _Oxycoccus palustre_ cover

oxpa_h: _Oxycoccus palustre_ height

phca_c: _Phyllodoce caerulea_ cover

phca_height: _Phyllodoce caerulea_ height

rhod_c: _Rhododendron_ sp. cover

rhod_h: _Rhododendron_ sp. height

saar_c: _Salix arctica_ cover

saar_h: _Salix arctica_ height

sagl_c: _Salix glauca_ cover

sagl_h: _Salix glauca_ height

sahe_c: _Salix herbacea_ cover

sahe_h: _Salix herbacea_ height

vaul_c: _Vaccinium uliginosum_ cover

vaul_h: _Vaccinium uliginosum_ height

drin_c: _Dryas integrifolia_ cover

drin_h: _Dryas integrifolia_ height

 Community characteristics

rich_pa: shrub species richness at grid level

rich: shrub species richness averaged across 3 1x1 m sampling quadrats per grid

shan: Shannon diversity (‘vegan’ package; Oksanen et al. 2019)

comp_bena: shrub competition on _Betula nana_ (i.e. summed cover of non-_B. nana_ species) [% cover]

comp_sagl: shrub competition on _Salix glauca_ (i.e. summed cover of non-_S. glauca_ species) [% cover]

ht_cwm: maximum shrub canopy height (i.e. community weighted mean by abundance of maximum shrub height) [mm]

 Topography and terrain wetness (calculated from)

tcws: Tasseled Cap Wetness component (Sentinal images)

slope: slope angle (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

twi: SAGA wetness index (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015)

asp: aspect (MEaSUREs Greenland Ice Mapping Project Digital Elevation Model v. 1; Howat et al., 2014, 2015) [º]

 Downscaled climate (time series 1984 – 2013)

insoljja: average June-August insolation [MJ cm-2 yr-1]

precipmam: cumulative March-May precipitation [mm]

tempcont: yearly maximum – minimum temperature [ºC]

tempmin: yearly minimum temperature [ºC]

precipjja: cumulative June-August precipitation [mm]

tempjja: average maximum June-August temperature [ºC]

tempmax: yearly maximum temperature [ºC]

en:
 Study Aim

We collected these data to alternatively train and validate high resolution (~ 90 m) Species Distribution Models (SDMs) and Species Abundance Models (SAMs) for _Betula nana_ L. (dwarf birch, Betulaceae) and _Salix glauca_ L. (grey willow, Salicaceae) in Southwest Greenland to assess how well such models can predict local-scale patterns.

 Data Description

Individual (presence-absence, abundance, maximum vegetative height) and community (species composition, maximum canopy height) shrub data for two fjords near Nuuk, Southwest Greenland. Also provided are corresponding downscaled climate data as well as calculated topographic and terrain wetness indicator variables.

 Nuup Kangerlua (Godthåbsfjord)

_Betula nana_ and _Salix glauca_ presence-absence, abundance, community species richness

 Kangerluarsunnguaq (Kobbefjord)

Shrub presence-absence, abundance, maximum vegetative height, community composition, maximum shrub canopy height

 Methods

 Field survey in Nuup Kangerlua

We conducted a stratified systematic plant survey along the length of Nuup Kangerlua (NK) fjord in Soutwesth Greenland (Fig. 1 in Chardon et al. 2022; following Nabe-Nielsen et al., 2017). At five distinct sites, we sampled along elevational gradients to collect data on presences, absences, abundance, and species composition of all woody species using a 0.7 x 0.7 m pin-point frame (Fig. 1e in Chardon et al. 2022). For model training, we converted these pin-point data to percent cover estimates based on the number of pins dropped (n = 25 per plot) and averaged them across the 119 spatio-climatic grids (see next section) corresponding to the plot locations (for details see Appendix S2 in Chardon et al. 2022).

 Field survey in Kangerluarsunnguaq

We conducted a random stratified plant survey in Kangerluarsunnguaq (K) fjord in Southwest Greenland. We used a preliminary Species Abundance Model trained with summed pin counts of _Betula nana_ in NK fjord (see Fig. S1.3 in Chardon et al. 2022) to stratify the ~ 27 x 17 km fjord landscape into low, medium, and high abundances classes. We randomly selected 90 x 90 m spatio-climatic grids to survey in each class for a total of 200 grids, ensuring that they were accessible by foot or boat (for details see Appendix S2 in Chardon et al. 2022). Within each grid, we sampled within three 1 m2 quadrats arranged in a randomly rotated equilateral triangle centered on the mid-point of the cell. We used a gridded sampling quadrat with 1% delineations (Fig. 1h in Chardon et al. 2022) to record woody species presences, absences, and composition, estimated percent cover, and measured maximum shrub species vegetatitve height. At every plot, we also visually scanned the area in a 20 m radius from the plot and recorded the presence of any additional shrub species to estimate grid-level species richness. As in NK fjord, we averaged these data at the grid level (for details see Appendix S2 in Chardon et al. 2022).

 Biotic variables

We calculated biotic microscale variables from the plant survey data collected in NK and K fjords. We calculated shrub species richness, diversity, and competition (i.e. sum of non-B. nana or non-S. glauca pin hits or percent cover). In K fjord, we also calculated canopy height as the community weighted mean (by abundance) of maximum vegetative shrub height.

 Climate variables

We computed high resolution temperature, precipitation, and insolation for local scale data for the study area by statistically downscaling climate time series (1982 - 2013) from the monthly CHELSA data (Karger et al. 2017). We downscaled these data from 30 arc sec (~ 400 m at the latitude of our study) to our target grid size of ~ 90 m with geographic weighted regression and using the MEaSUREs Greenland Ice Mapping Project (GIMP) Digital Elevation Model (DEM) v. 1 (Howat et al., 2014, 2015). We then calculated 30-year averages of the climate parameters: average summer (June – August) maximum temperature, yearly maximum temperature, yearly minimum temperature, temperature continentality (yearly max. - min. temperatures), cumulative Spring (March – May) precipitation, cumulative summer precipitation, and average summer incident solar radiation (henceforth, insolation) (for calculation details see Appendices S2, S3 in Chardon et al. 2022 and Appendix S2 in von Oppen et al. 2021).

 Topography and terrain wetness indicator variables

We calculated several topographic and terrain wetness indices at a local scale. We derived slope, aspect, and the SAGA wetness index (hereafter TWI; Boehner et al., 2002; Boehner and Selige, 2006) from the GIMP DEM. TWI is a measure of how ‘wet’ an area is, based on water drainage from the surrounding landscape. We also calculated the tasseled cap wetness component (hereafter TCW, Crist and Cicone 1984) from satellite images (for details see Appendices S2, S3 in Chardon et al. 2022) as an alternative measure of wetness.

 Computer code

Attached as zip file and available on GitLab (https://gitlab.com/nathaliechardon/gl_microclim)

 Third-party data

Data used to calculate climate, topography, and terrain wetness indicator variables are publicly available (see Appendix S2 in Chardon et al. 2022 for all data references).

accessible-shrub-forest

display_name: accessible-shrub-forest

arctic

display_name: arctic

climate

display_name: climate

field-data

display_name: field-data

sdm

display_name: sdm

validation

display_name: validation

en: Greenland shrubs and microclimate


maintainer: Varun Sharma

en: snow-atmosphere-interaction antarctica mesoscale-modelling blowing-snow alps

en: REPRODUCIBILITY_SIMULATION.zip contains all the extraneous files required to perform the simulations in associated manuscript (see Citations 3 in README.txt). The purpose of sharing these files apart from purposes of reproducibility is to also serve as templates for users to learn the workflow of CRYOWRF and perform simulations for their own purposes.

en:
REPRODUCIBILITY_POSTPROCESSING.zip contains data and scripts used to produce all figures from for the associated manuscript submitted for peer review ( see section on Citations below )

The folders stand for the three case studies presented in the manuscript. Two case studies (case Ia and case Ib) are on simulations in Antarctica and one case study
case II is for simulations in the Alps.

The associated data in the form of netcdf files or text files were produced using CRYOWRF v1.0

Apart from reproducibility, the intention of sharing these scripts is also to make it easy for new users to adopt CRYOWRF and be efficient, in particular in processing snowpack outputs.

The scripts were tested on multiple machines (different OSes) with the attached YAML file (wrf_analysis.yml) to setup the conda environment. All scripts were run using spyder as the IDE.

en:
This repository contains data required for reproducibility of the results to be published in the associated manuscript.

Apart from reproducibility, the attached datasets also serve as templates for new users to adopt CRYOWRF in their research.

The datasets consist of two folders organized in zip format:

1. REPRODUCIBILITY_SIMULATION:
  Consists of namelists for WPS, WRF and SNOWPACK to reproduce simulations published in the manuscript
  Additional files include datasets (from IMAU-FDM / RACMO, see "credits" below ) as well as helper python scripts to
  produce *.sno files
  which are used as initial conditions for SNOWPACK in CRYOWRF.

2. REPRODUCIBILITY_POSTPROCESSING:
   Includes outputs of CRYOWRF and python scripts used to prepare figures in the manuscript.

Each of the folders have their own readme files for more details.

 Code
citation: Varun Sharma. (2021, July 2). vsharma-next/CRYOWRF: CRYOWRF v1.0 (Version v1.0). Zenodo. http://doi.org/10.5281/zenodo.5060165

location: https://gitlabext.wsl.ch/atmospheric-models/CRYOWRF   (stable releases / institutional repo)
        https://github.com/vsharma-next/CRYOWRF            (dev branches / developer repo)

 Publication
**Introducing CRYOWRF v1.0: Multiscale atmospheric flow simulations with advanced snow cover modelling.**
Varun Sharma, Fraziska Gerber and Michael Lehning, Submitted to Geoscientific Model Development

 Acknowledgements
We thank Peter Kuipers Munneke (P.KuipersMunneke@uu.nl) for preparing and sharing outputs of IMAU-FDM and RACMO used for initial conditions for case Ia. The relevant citations for the methods through which these datasets were generated are:

* Kuipers Munneke, P., S. R. M. Ligtenberg, B. P. Y. Noël, I. M. Howat, J. E. Box, E. Mosley-Thompson, J. R. McConnell, K. Steffen, J. T. Harper, S. B. Das and M. R. van den Broeke. 2015. Elevation change of the Greenland ice sheet due to surface mass balance and firn processes, 1960-2014. The Cryosphere, 9, 2009-2025. doi:10.5194/tc-9-2009-2015

 * Ligtenberg, S. R. M., P. Kuipers Munneke, B. P. Y. Noël, and M. R. van den Broeke. 2018. Brief communication: Improved simulation of the present-day Greenland firn layer (1960-2016). The Cryosphere, 12, doi:10.5194/tc-12-1643-2018

alps

display_name: alps

antarctica

display_name: antarctica

blowing-snow

display_name: blowing-snow

mesoscale-modelling

display_name: mesoscale-modelling

snow-atmosphere-interaction

display_name: snow-atmosphere-interaction

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: Reproducibility Dataset for CRYOWRF v1.0


en: pressure-perturbations landslides tramm cces pore-water landslide-testsite

en: TRAMM Ruedlingen data

en: TRAMM project data from the Ruedlingen site (exp 1 + 2). Measured parameters include soil suction, water level, soil volumetric water content, rain intensity, and soil temperature.

en:
A landslide testsite dataset related to pore water pressure perturbations on the stability of unsaturated silty sand slopes leading to the initiation and propagation of the shear deformations and eventual rapid mass movements. This project was initiated and led by the Institute of Geotechnical Engineering (IGT) of the Swiss Federal Institute of Technology (ETH Zurich) and was incorporated in a Swiss national (TRAMM) and a European Union (SafeLand) multidisciplinary research project. Field site: The experimental slope is 7.5 m wide by 35 m long, located in the Swiss lowlands on an east facing slope over-looking the river Rhine, at an altitude of ~ 350 masl. Originally there were forestry covertures of circa 80%, heights of 5-20 m. Shrubs up to 1-5 m high and a free herb layer covered ~ 50% of the surface. The average gradient was determined to be from 38° to 43° with a slightly concave surface. The underlying rock consists mainly of Molasse, which is formed by alternate layers of sea deposits under the Tethys Sea (Seawater Molasse) and land deposits (Freshwater Molasse). Several augured samples, as well as an outcrop of the bedrock about 20 m above the selected field, revealed horizontal layering of fine grained sand- and marlstone at the test site. The sandstone was later proven to be highly permeable and fissured. Grain-size distributions were determined and the soil was classified as medium-low plasticity silty sand. Site instrumentation:Measurements of soil suction, groundwater level, soil volumetric water content, rain intensity and soil temperature were taken and combined with geophysical monitoring using Electrical Resistance Tomography (ERT) and investigations into subsurface flow by means of tracer experiments. Deformations were monitored during the experiment, both on the surface via photogrammetrical methods and within the soil mass, using a flexible probe equipped with strain gauges at different points and two axis inclinometers on the top and acoustic sensors. Instruments were installed mainly in three clusters at depths of 15, 30, 60, 90, 120, and 150 cm below the ground surface over the slope, including jet-fill tensiometers, TDRs, Decagon TDRs, piezometers, soil temperature sensors, deformation probes, earth pressure cells, acoustic sensors and rain gauges. A ring-net barrier (provided by Geobrugg AG) was set up at the foot of the slope to protect the road. Experiments: A sprinkling experiment was carried out in September 2008 to investigate the hydrological and mechanical response of the slope (Experiment 1), followed by a second one to trigger a landslide in March 2009 (Experiment 2).

Publications

1. Lehmann, P., F. Gambazzi, B. Suski, L. Baron, A. Askarinejad, S. M. Springman, K. Holliger, and D. Or (2013), Evolution of soil wetting patterns preceding a hydrologically induced landslide inferred from electrical resistivity survey and point measurements of volumetric water content and pore water pressure, Water Resour. Res., 49, 7992–8004, doi:[10.1002/2013WR014560](http://dx.doi.org/10.1002/2013WR014560).

2. Springman, S. M., Kienzler, P., Casini, F., & Askarinejad, A. (2009). Landslide triggering experiment in a steep forested slope in Switzerland. In 17th International Conference of Soil Mechanics and Geotechnical Engineering, Alexandria, Egypt (pp. 1698-1701). doi: [10.3233/978-1-60750-031-5-1698](http://dx.doi.org/10.3233/978-1-60750-031-5-1698)

cces

display_name: cces

landslide-testsite

display_name: landslide-testsite

landslides

display_name: landslides

pore-water

display_name: pore-water

pressure-perturbations

display_name: pressure-perturbations

tramm

display_name: tramm

en: TRAMM project Ruedlingen experimental landslide dataset, Switzerland


causal-effect-of-mfoz

title_for_slug: causal-effect-of-mfoz

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

en: Data and Stata code for repeating the results

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

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

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

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

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

built-up-land-change

display_name: built-up-land-change

causality

display_name: causality

fujian-province

display_name: fujian-province

propensity-score-matching

display_name: propensity-score-matching

selection-bias

display_name: selection-bias

the-major-function-oriented-zone

display_name: the-major-function-oriented-zone

en: Causal effect of MFOZ


backward-trajectories

maintainer: Michael Sprenger

title_for_slug: backward-trajectories

en: alps raclets cosmo atmospheric-dynamics lagrangian-analysis backward-trajectories complex-topopgraphy

en: Backward Trajectories

en: The trajectories are started from surface up to 5000 m altitude, with an equidistant vertical spacing of 200 m, and calculated 120 h backward in time. As long as the trajectories stay in the COSMO1 domain (essentially the Alpine region), as run operationally at MeteoSwiss, the COSMO analysis wind fields are used. When the trajectories leave this domain, they are extended further backward in time based on ECMWF analysis data. The following fields are traced along the trajectories: longitude, latitude, geopotential height, pressure, specific humidity, temperature and surface topography. Furthermore, it is indicated for each trajectory time instance on which model data, COSMO or ECMWF, the trajectory time is based. Finally, the trajectories are started every hour in the specified time period, and positions and fields are available every 10 min. Each day is saved as a separate file. The last 8 digits of the file names indicate the starting date of the trajectories.

en: Backward trajectories were calculated from two positions: Davos Wolfgang (LON: 9.85361, LAT: 46.83551) and Weissfluhjoch (LON: 9.80646 LAT: 46.83304) for the time period February 2 until March 27 2019 using COSMO or ECMWF, respectively.

alps

display_name: alps

atmospheric-dynamics

display_name: atmospheric-dynamics

backward-trajectories

display_name: backward-trajectories

complex-topopgraphy

display_name: complex-topopgraphy

cosmo

display_name: cosmo

lagrangian-analysis

display_name: lagrangian-analysis

raclets

display_name: raclets

en: Backward  Trajectories


vegetation-height-model-nfi

maintainer: Christian Ginzler

title_for_slug: vegetation-height-model-nfi

en: vegetation-height-model photogrammety point-cloud 3d forest

en: 2012 Vegetation Height Model NFI

media_type: geotiff

en:
Data 'Vegetation Height Model NFI 2012' available on request (christian.ginzler@wsl.ch).

A national vegetation height model was calculated for Switzerland using digital aerial images. We used the stereo aerial images acquired by the Federal Office of Topography swisstopo using the ADS80 sensor to first calculate a digital surface model (DSM) with a very high spatial resolution (1 × 1 m). The DSM was then normalized to obtain the actual vegetation heights using the digital terrain model (DTM) SwissAlti3D based on airborne laser data. Buildings are masked out using the NDVI information of the aerial images to produce a vegetation height model, and to produce a vegetation height model (VHM). Such a model will be calculated in the framework of the Swiss National Forest Inventory (NFI) with consistent methods and a very high level of detail. For covering the whole of Switzerland, we used summer aerial images from 2007 to 2012.

format: GeoTIFF

en: 2019 Vegetation Height Model NFI (current)

en:
Data 'Vegetation Height Model NFI 2019' available on request (christian.ginzler@wsl.ch).

A national vegetation height model was calculated for Switzerland using digital aerial images. We used the stereo aerial images acquired by the Federal Office of Topography swisstopo using the ADS80 sensor to first calculate a digital surface model (DSM) with a very high spatial resolution (1 × 1 m). The DSM was then normalized to obtain the actual vegetation heights using the digital terrain model (DTM) SwissAlti3D based on airborne laser data. Buildings are masked out using the NDVI information of the aerial images to produce a vegetation height model, and to produce a vegetation height model (VHM). Such a model will be calculated in the framework of the Swiss National Forest Inventory (NFI) with consistent methods and a very high level of detail. For covering the whole of Switzerland, we used summer aerial images from 2013 to 2019.

en: 2016 Vegetation Height Model NFI

en:
Data 'Vegetation Height Model NFI 2016' available on request (christian.ginzler@wsl.ch).

A national vegetation height model was calculated for Switzerland using digital aerial images. We used the stereo aerial images acquired by the Federal Office of Topography swisstopo using the ADS80 sensor to first calculate a digital surface model (DSM) with a very high spatial resolution (1 × 1 m). The DSM was then normalized to obtain the actual vegetation heights using the digital terrain model (DTM) SwissAlti3D based on airborne laser data. Buildings are masked out using the NDVI information of the aerial images to produce a vegetation height model, and to produce a vegetation height model (VHM). Such a model will be calculated in the framework of the Swiss National Forest Inventory (NFI) with consistent methods and a very high level of detail. For covering the whole of Switzerland, we used summer aerial images from 2011 to 2016.

en: 2015 Vegetation Height Model NFI

en:
Data 'Vegetation Height Model NFI 2015' available on request (christian.ginzler@wsl.ch).

A national vegetation height model was calculated for Switzerland using digital aerial images. We used the stereo aerial images acquired by the Federal Office of Topography swisstopo using the ADS80 sensor to first calculate a digital surface model (DSM) with a very high spatial resolution (1 × 1 m). The DSM was then normalized to obtain the actual vegetation heights using the digital terrain model (DTM) SwissAlti3D based on airborne laser data. Buildings are masked out using the NDVI information of the aerial images to produce a vegetation height model, and to produce a vegetation height model (VHM). Such a model will be calculated in the framework of the Swiss National Forest Inventory (NFI) with consistent methods and a very high level of detail. For covering the whole of Switzerland, we used summer aerial images from 2010 to 2015.

en: Countrywide Stereo-Image Matching for Updating Digital Surface Models in the Framework of the Swiss National Forest Inventory

en: Thumbnail

en:
A national vegetation height model was calculated for Switzerland using digital aerial images.
We used the stereo aerial images acquired by the Federal Office of Topography swisstopo using the ADS80 sensor to first calculate a digital surface model (DSM) with a very high spatial resolution
(1 × 1 m and 0.5 x 0.5 m). The DSM was then normalized to obtain the actual
vegetation heights using a digital terrain model (DTM)
based on laser data with the buildings masked out, and to
produce a vegetation height model (VHM). Such a model
will be calculated in the framework of the Swiss National Forest Inventory (NFI) with consistent methods and a very
high level of detail. For covering the whole of Switzerland,
we use summer aerial images from six years.

Latest version is from 2019.

forest

display_name: forest

photogrammety

display_name: photogrammety

point-cloud

display_name: point-cloud

vegetation-height-model

display_name: vegetation-height-model

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: Vegetation Height Model NFI