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

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

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

Paper Citation:

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

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

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

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

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

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

bacteria

display_name: bacteria

biogeochemical-cycling

display_name: biogeochemical-cycling

carcass

display_name: carcass

carrion

display_name: carrion

decomposition

display_name: decomposition

detritus

display_name: detritus

fungi

display_name: fungi

soil-biodiversity

display_name: soil-biodiversity

wolf

display_name: wolf

yellowstone

display_name: yellowstone

en: Elk and bison carcasses in Yellowstone, USA


maintainer: Lisa Hülsmann

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

state: active

type: dataset

en: Calibration data

media_type: text/csv

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

 Meta data
Reserve: Name of the forest reserve

Plot: Permanent plot number

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

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

dt: Length of the mortality period (years)

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

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

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

mDBH: Mean diameter at breast height (mm)

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

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

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

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

P: Mean annual precipitation sum (mm)

mT: Mean annual air temperature (°C)

format: CSV

mimetype: text/csv

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

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

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

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

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

empirical-mortality-models

display_name: empirical-mortality-models

forest-reserves

display_name: forest-reserves

inventory-data

display_name: inventory-data

tree-growth

display_name: tree-growth

tree-mortality

display_name: tree-mortality

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


maintainer: Athanasios Nenes

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

en: CCN data Weissfluhjoch

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

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

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

_Error_code_:

0000: Measurements without error

1000: Period of transient operation after a supersaturation change

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

3000: Irregular values in the room temperature housing the instrument

9999: Missing data

en: Predicted cloud droplet number concentrations Weissfluhjoch

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

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

en: Hygroscopicity parameter Weissfluhjoch

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

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

Missing values are indicated by "-999".

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

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

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

ccn

display_name: ccn

cloud-droplets

display_name: cloud-droplets

cloud-modeling

display_name: cloud-modeling

dmt-ccn-chamber

display_name: dmt-ccn-chamber

hygroscopicity

display_name: hygroscopicity

raclets

display_name: raclets

sigma-w

display_name: sigma-w

supersaturation

display_name: supersaturation

weissfluhjoch

display_name: weissfluhjoch

en: CCN, hygroscopicity, predicted cloud droplet numbers Weissfluhjoch


en: primeval-forests forest-reserves tree-mortality fagus-sylvatica mortality-models forest-inventory

en: Swiss data full

en:
READ ME for

Hülsmann L, Bugmann H, Commarmot B, Meyer P, Zimmermann S, Brang P
Data from: Does one model fit all? Patterns of beech mortality in natural forests of three European regions.
Ecological Applications http://dx.doi.org/10.5061/dryad.h4s6t

Forest inventory data of beech for the development of individual-tree mortality models from strict forests reserves in Switzerland, Germany and Ukraine.

The data package provides the files
CH_reduced: inventory data of the reduced dataset from Switzerland that was used to assess the influence of additional environmental and stand characteristics on beech mortality
CH: inventory data from strict forest reserves in Switzerland
DE: inventory data from strict forest reserves in Germany
UA: inventory data from a primeval beech forest in Ukraine

A set of three consecutive inventories was used to generate records for the calibration of mortality models based on trees that were alive in the first and second inventory and either dead or alive in the third inventory.
Diameter at breast height: dbh (mm)
Annual relative basal area increment: relBAI (), calculated based on the first and the second dbh measurement as the compound annual growth rate of the trees basal area (Bigler and Bugmann 2004)
Status of the tree: y, 0 = alive, 1 = dead
Length of the mortality period: dt (years)

Additional environmental and stand characteristics were derived for the reduced Swiss dataset.
Note the varying seasonal level of the covariates considering the entire year, the growing season from April to September (G) and the winter month January to March (W).
All climatic covariates were calculated for both the ‘growth’ (2) and the ‘mortality’ period (3). The stand covariates were derived for the second inventory of each record (2).

Mean diameter at breast height: dbh.mean.2 (mm)
Interquartile range of diameter at breast height: dbh.IQR.2 (mm)
Median diameter at breast height: dbh.median.2 (mm)
Basal area per hectare: BA.2 (m2 ha-1)
Number of trees per hectare: N.2 (ha-1)
pH-value: pH ()
Available water capacity in the soil: AWC (mm)
Mean annual precipitation sum: P.2, P.3 (mm)
Mean precipitation sum during growing season: PG.2, PG.3 (mm)
Mean annual air temperature: mT.2, mT.3 (°C)
Mean air temperature during growing season: mTG.2, mTG.3 (°C)
Mean air temperature during winter months: mTW.2, mTW.3 (°C)
Mean growing degree-days: GDD.2, GDD.3 ()
Mean number of months during growing season with water deficit: mDEFdurG.2, mDEFdurG.3 ()
Mean annual maximum water deficit: mDEFmax.2, mDEFmax.3 (mm)
Overall maximum of the annual maximum water deficit: maxDEFmax.2, maxDEFmax.3 (mm)
Mean of cumulative water deficit during growing season: mDEFcumG.2, mDEFcumG.3 (mm)
Overall maximum of cumulative water deficitduring growing season: maxDEFcumG.2, maxDEFcumG.3 (mm)
Percentage of years with at least one month water deficit: percDEF.2, percDEF.3 (%)

en: German data

en: Swiss data reduced

en: Ukrainian data

en:
The datasets comprise nearly 19’000 trees of European beech (_Fagus sylvatica_ L.) from unmanaged forests in Switzerland, Germany / Lower Saxony and Ukraine. Tree death was modelled as a function of size and growth, i.e., stem diameter (DBH) and relative basal area increment (relBAI). To explain the spatial and temporal variability in mortality patterns, we considered a large set of environmental and stand characteristics.

 Inventory data
The strict forest reserves in Switzerland and Germany had been established in the period of 1961-1975 and 1971-1974, respectively. Every reserve included up to 10 permanent plots ranging from 0.09 to 1.8 ha in size, with slightly irregular re-measurement intervals. Permanent plots with pure or mixed beech stands were selected from the reserves of both networks. Reserves with considerable wind disturbance during the monitored intervals were excluded from the analysis. In addition to data from the Swiss and German reserves, data from a 10 ha plot in the primeval beech forest Uholka in Western Ukraine including three remeasurements were used. The inventory data provide diameter measurements at breast height (dbh) for revisited trees with a diameter of more than 4, 7 and 6 cm for Switzerland, Germany and Ukraine, respectively.

 Mortality predictors
A set of three consecutive inventories was used to generate records for the calibration of mortality models based on trees that were alive in the first and second inventory and either dead or alive in the third inventory. As an explanatory variable, the annual relative basal area increment (relBAI) was calculated based on the first and the second dbh measurement as the compound annual growth rate of the trees basal area. Tree dbh in the second inventory was used in addition to relBAI to model tree status (alive or dead) of the third inventory.

To increase the generality of the mortality models, we selected environmental variables that are known to have a considerable influence on growth and mortality of beech. We emphasized the effects of water availability using a large set of drought characteristics that were calculated based on the local site water balance. We also related beech mortality to soil pH, temperature, precipitation and growing degree-days. Additionally, we considered stand characteristics that reflect the development stage, competition and structure of the forests.

 Further information
For further information, refer to Hülsmann _et al_. (2016) Does one model fit all? patterns of beech mortality in natural forests of three European regions. _Ecological Applications_.

fagus-sylvatica

display_name: fagus-sylvatica

forest-inventory

display_name: forest-inventory

mortality-models

display_name: mortality-models

primeval-forests

display_name: primeval-forests

en: Data from: Does one model fit all? patterns of beech mortality in natural forests of three European regions


maintainer: Roksolana Petrovska

en: fagus lar regeneration dominance shade-tolerance acer mortality-rate

media_type: r script (flat-ascii text)

en: The R-script of survival analysis for accelerated failure time model

format: SERVICE

en: Mortality of regeneration

media_type: CSV

en: Acer platanoides, Acer pseudoplatanus and Fagus sylvatica living and dead juveniles (0-500 cm in height) grown in the deep shade of a prime beech-dominated forest.

en:
 One individual per species, vitality class (low and high) and height class (eight classes: 0–10, 11–20, 21–35, 36–60, 61–90, 91–130, 131–200 and 201–500 cm) was randomly selected and harvested in each of the six plots. This resulted in a sample of 82, 80 and 89 living individuals of A. platanoides, A. pseudoplatanus and F. sylvatica, respectively.

 Additionally stems of dead Acer spp. and F. sylvatica trees that had died within the last three years (2015–2018) were randomly harvested, matching the height classes of the harvested living trees wherever possible. In total, 179 dead young trees (60 A. platanoides, 72 A. pseudoplatanus and 47 F. sylvatica) were collected.

 Variables:
* species_code: a_pla - Acer platanoides, a_pse - Acer pseudoplatanus, f_syl - Fagus sylvatica
* species: as above
* dummy: 0 - living individual, 1 - dead individual
* LAR_cm2_g: leaf area ratio or ratio of leaf area to total plant biomass, [cm2/g]
* tree_age: in years
* avg_ring_micron: average width of the last 5 rings in tree life excluding the last ring
* dry_mass_g: aboveground and belowground biomass
* DLI: direct light index (measured only under living individuals)
* BLI: diffuse light index (measured only under living individuals)
* GLI: global light index

acer

display_name: acer

dominance

display_name: dominance

fagus

display_name: fagus

lar

display_name: lar

mortality-rate

display_name: mortality-rate

regeneration

display_name: regeneration

shade-tolerance

display_name: shade-tolerance

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: Mortality of regeneration: Acer spp. and Fagus sylvatica


greenland-shrubs-and-microclimate

maintainer: Nathalie Chardon

title_for_slug: greenland-shrubs-and-microclimate

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

en: Kangerluarsunnguaq

media_type: .csv

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

en: shortwave-radiation incoming outgoing weissfluhjoch longwave-radiation

en: Parameter file

media_type: text/plain

format: TXT

mimetype: text/plain

en: Metadata

media_type: application/pdf

format: PDF

mimetype: application/pdf

en: DATA ACCESS

en: Processing doc

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

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

incoming

display_name: incoming

longwave-radiation

display_name: longwave-radiation

outgoing

display_name: outgoing

shortwave-radiation

display_name: shortwave-radiation

weissfluhjoch

display_name: weissfluhjoch

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


maintainer: Frank Techel

en: snow avalanches alps observations shapefile

en: Selected wet snow avalanche activity data Davos, Switzerland (2011-2014)

media_type: application/zip

en: Zip archive with shapefiles for the selected wet snow avalanche polygons, including metadata.

format: ZIP

mimetype: application/zip

en:
Polygons of wet snow avalanches in the Davos area, as documented by the Swiss avalanche warning service. The georeferenced outlines of the avalanches contain both the release as well as the deposit area, but without separating between both.

The dataset is a subset of the total record of 1615 avalanches classified as wet snow avalanches from October 2011 - September 2014, containing those 255 avalanches exceeding 0.0125 km^2. Every polygon comes with meta data, including the date of occurrence.

This dataset is the underlying dataset to:
Wever, N., Vera Valero, C. and Techel, F. (2018) _Coupled snow cover and avalanche dynamics simulations to evaluate wet snow avalanche activity_. Submitted to J. Geophys. Res., in review.

alps

display_name: alps

avalanches

display_name: avalanches

observations

display_name: observations

shapefile

display_name: shapefile

snow

display_name: snow

education

en: Education and science

envidat

en: EnviDat: the Environmental Data Portal

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


maintainer: Axel Volkwein

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

en: Jumps and impacts

media_type: ascii, tab separated, also as excel table

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

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

en: Overview of all tests

media_type: ZIP

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

en: Block trajectories recorded from local positioning system

media_type: ascii, tab separated, also as Excel table

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

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

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

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

experiment

display_name: experiment

experiments

display_name: experiments

field-test

display_name: field-test

field-tests

display_name: field-tests

mass-movements

display_name: mass-movements

natural-hazards

display_name: natural-hazards

rockfall

display_name: rockfall

trajectory

display_name: trajectory

en: Repetitive trajectory testing in Tschamut 2014


en: weisfluhjoch longwave-radiation shortwave-radiation radiation davos

en: Parameter file

en: Description of ASRB parameters.

en: Metadata

en: DATA ACCESS

media_type: SMET

en: ASRB radiation station at the Weissfluhjoch versuchsfeld.

en: Processing doc

en: Description of processing code that was ported to R by P. Frei (2016).

en:
Incoming and outgoing shortwave and longwave 2 min radiation measurements at the Weissfluhjoch research site, Davos, CH. The experimental site at the Weissfluhjoch (WFJ, 46.83 N, 9.81 E) is located at an altitude of 2540 m in the Swiss Alps near Davos. During the winter months, almost all precipitation falls as snow at this altitude. As a consequence, a continuous seasonal snow cover builds up every winter, with a maximum snow height ranging from 153–366 cm over the period 1934–2012. The measurement site is located in an almost flat part of a southeast oriented slope.

 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

longwave-radiation

display_name: longwave-radiation

radiation

display_name: radiation

weisfluhjoch

display_name: weisfluhjoch

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


maintainer: Benjamin Reuter

en: avalanche-problem-type avalanche avalanche-activity crocus snowpack snow-cover-modeling snow-instability avalanche-danger weissfluhjoch snow

en: 2015-2016 Avalanche problem types from SNOWPACK

en:
Avalanche problem types were derived from a snow cover simulation with the model SNOWAPCK for the winter season 2015-2016 at the Weissfluhjoch study plot, Davos, CH.
The data cover the period 01.01.2016-30.04.2016. The binary indicator (0,1) tells whether an avalanche problem type was identified in the snow cover simulation (1) or not (0).
The data are presented in Figure 7 in Reuter et al. (2021).

en: 2015-2016 Avalanche problem types from Crocus

en:
Avalanche problem types were derived from a snow cover simulation with the model Crocus for the winter season 2015-2016 at the Weissfluhjoch study plot, Davos, CH.
The data cover the period 01.01.2016-30.04.2016. The binary indicator (0,1) tells whether an avalanche problem type was identified in the snow cover simulation (1) or not (0).
The data are presented in Figure 7 in Reuter et al. (2021).

en: 2015-2016 Seismic avalanche acitivity catalogue

en:
Avalanches were indentified in seismic data recorded by two seismic arrays, one deployed at Wannengrat (2,500 m) and another one deployed in the Dischma valley (2,000 m).
The data cover the period from 01.01.2016 - 30.04.2016 and are presented in Figure 7 in Reuter et al. (2021).

en: 1999-2017 Avalanche problem types from Crocus

en:
Avalanche problem types were derived from snow cover simulations with the model Crocus at the Weissfluhjoch study plot, Davos, CH.
The data cover the winter seasons between 1999-2017 and include the simulated annual frequencies of new snow, wind slab, persistent weak layer, aging-persistent weak layer and wet snow problem types.
The data are presented in Figure 8 in Reuter et al. (2021).

en: 1999-2017 Avalanche problem types from SNOWPACK

en:
Avalanche problem types were derived from snow cover simulations with the model SNOWPACK at the Weissfluhjoch study plot, Davos, CH.
The data cover the winter seasons between 1999-2017 and include the simulated annual frequencies of new snow, wind slab, persistent weak layer, aging-persistent weak layer and wet snow problem types.
The data are shown in Figure 8 in Reuter et al. (2021).

en:
Avalanche problem types were derived from snow cover simulations with the models Crocus and SNOWPACK at the Weissfluhjoch study plot, Davos, CH. The data include annual frequencies of avalanche problem types for the seasons 1999-2017 and daily presence of avalanche problem types for the period 01.01.2016 - 30.04.2016.
Avalanche activity was derived from two seismic sensor arrays deployed no further than 15 km from Weissfluhjoch, Davos, CH. The data cover the period 01.01.2016 - 30.04.2016.

avalanche

display_name: avalanche

avalanche-activity

display_name: avalanche-activity

avalanche-danger

display_name: avalanche-danger

avalanche-problem-type

display_name: avalanche-problem-type

crocus

display_name: crocus

snow-cover-modeling

display_name: snow-cover-modeling

snow-instability

display_name: snow-instability

snowpack

display_name: snowpack

weissfluhjoch

display_name: weissfluhjoch

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: Simulated avalanche problem types and seismic avalanche activity around Weissfluhjoch


en: climate maximum-temperature minimum-temperature precipitation century timeseries free-climate-data climate-layers air-temperature

en: CHELSA Technical Documentation

en: As the CHELSA project is continuously expanded, certain new parameters, time scales etc. are not included anymore in the original publication (Karger et al. 2017). We are therefore provide this document as a guideline of the current climatic parameters available and will expand it as new parameters or time periods will become available. Please refer to the version history for changes made in the document.

en: Data Access

media_type: geotiff

format: GeoTIFF

en:
CHELSAcruts is a delta change monthly climate dataset for the years 1901-2016 for mean monthly maximum temperatures, mean monthly minimum temperatures, and monthly precipitation sum. Here we use the delta change method by B-spline interpolation of anomalies (deltas) of the CRU TS 4.01 dataset. Anomalies were interpolated between all CRU TS grid cells and are then added (for temperature variables) or multiplied (in case of precipitation) to high resolution climate data from CHELSA V1.2 (Karger et al. 2017, Scientific Data). This method has the assumption that climate only varies on the scale of the coarser (CRU TS) dataset, and the spatial pattern (from CHELSA) is consistent over time. This is certainly a rather crude assumption, and for time periods for which more accurate data is available CHELSAcruts should be avoided if possible (e.g. use CHELSA V1.2 for 1979-2015). Different to CHELSA V1.2, CHELSAcruts does not take changing wind patterns, or temperature lapse rates into account, but rather expects them to be constant over time, and similar to the long term averages.

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

air-temperature

display_name: air-temperature

century

display_name: century

climate

display_name: climate

climate-layers

display_name: climate-layers

free-climate-data

display_name: free-climate-data

maximum-temperature

display_name: maximum-temperature

minimum-temperature

display_name: minimum-temperature

precipitation

display_name: precipitation

timeseries

display_name: timeseries

en: CHELSAcruts -  High resolution temperature and precipitation timeseries for the 20th century and beyond


backward-trajectories

maintainer: Michael Sprenger

title_for_slug: backward-trajectories

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

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


en: permafrost ground-temperature ritigraben borehole permos

en: DATA ACCESS

media_type: URL

en: Map of PERMOS borehole network

en: PERMOS DATA POLICY

media_type: PDF

en:
_ENVIDAT NOTE: Data currently unvailable and measures are being taken to recover and restore the data files._

Processed ground temperature measurements at the Ritigraben permafrost borehole (RIT_0102) in canton Valais, Switzerland. The borehole is located at 2690 m asl on a flat site. The surface material is coarse blocks and borehole depth is 30 m. Thermistors used YSI 44006. Year of drilling 2002. This borehole is part of the Swiss Permafrost network, PERMOS (www.permos.ch). Contact phillips@slf.ch for details of processing applied.

borehole

display_name: borehole

ground-temperature

display_name: ground-temperature

permafrost

display_name: permafrost

permos

display_name: permos

ritigraben

display_name: ritigraben

en: RIT1: Processed permafrost borehole data (2690  m asl), Ritigraben, Switzerland


en: particles turbulence hill dust atmosphere deposition snow

en: preferential deposition patterns of all simulations performed

en: The data consist of matlab files. Each files contains a 2D array with the values of deposition at the surface nodes. Each array has Nx x Ny values. Please refer to the original paper for additional information on the spatial resolution and grid size.

en:
Preferential deposition of snow and dust over complex terrain is responsible for a wide range of environmental processes, and accounts for a significant source of uncertainty in surface mass balances of cold and arid regions. Despite the growing body of literature on the subject, previous studies reported contradictory results on the location and magnitude of deposition maxima and minima. This study aims at unraveling the governing processes of preferential deposition in neutrally stable atmosphere and to reconcile seemingly inconsistent results of previous works. For this purpose, a comprehensive modeling approach is developed, based on large eddy simulations of the turbulent airflow, Lagrangian stochastic model of particle trajectories, and immersed-boundary method to represent the underlying topography. The model performance is tested against wind tunnel measurements of dust deposition around isolated and sequential hills. A scale analysis is then performed to investigate the dependence of snowfall deposition on the particle Froude and Stokes numbers, which fully account for the governing processes of inertia, flow advection, and gravity. Additional simulations are performed, to test whether the often used assumption of inertialess particles yields accurate deposition patterns. We finally show that our scale analysis provides qualitatively similar results for hills with different aspect ratios.

This dataset contains the results of the LES-LSM model. Each Matlab file contains a 2D array of deposition values (in kg/m2) in each surface node (ix, iy) of the Cartesian grid. The file names are consistent with the simulation numbers listed in the original paper. For additional information, please refer to "Preferential deposition of snow and dust over hills: governing processes and relevant scales" by F. Comola, M. G. Giometto, S. T. Salesky, M. B. Parlange, and M. Lehning, Journal of Geophysical Research: Atmospheres, 2019.

atmosphere

display_name: atmosphere

deposition

display_name: deposition

dust

display_name: dust

hill

display_name: hill

particles

display_name: particles

snow

display_name: snow

turbulence

display_name: turbulence

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: forest-persistence urban-morphology built-up-classification planning-outcome urban-compactness forest-gain nature-conservation urban-sprawl land-use-and-cover-change spatial-modelling basic-services topography

en: Increase in clean water provision service

en: Population density per pixel (change per year)

en: Urban persistence: Built-up 2005 = Built-up 2015

en: Human Development Index (HDI) per municipality (change per year)

en: slope

en: Slope °

media_type: .shp

en: Master plan Ubatuba (zones)

media_type: SHP

en: Master plan Ihabela (zones)

format: SHP

en: isolated

en: New built-up pixels classified as isolated

en: Permanent housing density increase

en: Increase in waste collection service

en: Variables description

en: Table of content. Variables description, source, and other important information

en: Ecological-Economic Zoning (EEZ)

en: Descricao dos dados em portiuguês

en: Portugues description

en: Mean income in Reais (R$) per census sector (change per year)

en: Axial

en: New built-up pixels classified as axial

en: Increase in Normalized Difference Vegetation Index

en: Native forest cover gains from 2005 to 2015

en: Native forest cover persistence from 2005 to 2015

en: infill

en: New built-up pixels classified as infill

en: Increase in sanitation service

en: Topographic Index Position (TPI)

en:
Basic education in % per census sector (change per
year)

en:
The present dataset is part of the published scientific paper entitled “The role of spatial planning in land change: An assessment of urban planning and nature conservation efficiency at the southeastern coast of Brazil” (Pierri Daunt, Inostroza and Hersperger, 2021). In this work, we evaluated the conformance of stated spatial planning goals and the outcomes in terms of urban compactness, basic services and housing provision, and nature conservation for different land-use strategies. We evaluate the 2005 Ecological-Economic Zoning (EEZ) and two municipal master plans from 2006 in a coastal region in São Paulo State, Brazil. We used Partial Least Squares Path Modelling (PLS-PM) to explain the relationship between the plan strategies and land-use change ten years after implementation in terms of urban compactness, basic services and housing increase, and nature conservation.

We acquired the data for the explanatory variables from different sources listed on Table 1. Since the model is spatially explicit, all input data were transformed to a 30 m resolution raster. Regarding the evaluated spatial plans, we acquired the zones limits from the São Paulo State Environmental Planning Division (CPLA-SP), Ilhabela and Ubatuba municipality.

1)	Land use and cover data: Urban persistence, Urban axial, Urban infill, Urban Isolates, Forest cover persistence, Forest cover gain, NDVI increase

We acquired two Landsat Collection 1 Higher-Level Surface Reflectance images distributed by the U.S. Geological Survey (USGS), covering the entire study area (paths 76 and 77, row 220, WRS-2 reference system, https://earthexplorer.usgs.gov/). We classified one image acquired by the Landsat 5 Thematic Mapper (TM) sensor on 2005-05-150, and one image from the Landsat 8 Operational Land Imager (OLI) sensor from 2015-08-15. We collected 100 samples for forest cover, 100 samples for built-up cover and 100 samples for other classes. We then classified these three classes of land cover at each image date using the Support Vector Machine (SVM) supervised algorithm (Hsu et al., 2003), using ENVI 5.0 software.

Land-use and land-cover changes from 2005 to 2015 were quantified using map algebra, by mathematically adding them together in pairs (10*LULC2015 + LULC2005). We reclassified the LULC data into forest gain (conversion of any 2005 LULC to forest cover in 2015); forest persistence (2005 forested pixels that remained forested in 2015); new built-up area (conversion of any 2005 LULC to built-up in 2015); and urban maintenance (2005 built-up pixels that remained built-up in 2015).

To describe the spatial configuration of the urban expansion, we classified the new built-up areas into axial, infill and isolated, following Inostroza et al. (2013) (For details, please refer to Supplementary Material I at the original publication).

The NDVI was obtained from the same source used for the LULC data. With the Google Engine platform, we used an annual average for the best pixels (without clouds) for 2005 and 2015, and we calculated the changes between dates. We used increases of > 0.2 NDVI to represent an improvement in forest quality.

2)	Federal Census data organization: Urban Basic Services and Housing indicator, socioeconomic and population:

The data used to infer the values of basic services provision, socioeconomic and population drivers was derived from the Brazilian National Census data (IBGE, 2000 and 2010). Population density, permanent housing unit density, mean income, basic education, and the percentage of houses receiving waste collection, sanitation and water provision services, called basic services in the context of this study, were calculated per 30 m pixel. The Human Development Index is only available at the municipality level. We attributed the HDI for the vector file with the municipality border, and we rasterized (30 m resolution) this file in QGIS. Annual rates of change were then calculated to allow comparability between LULC periods. To infer the BSH, we used only areas with an increase in permanent housing density and basic services provision (See Supplementary Material I at the original publication).

3)	Topographic drivers

To infer the values of the topographic driver, we used the slope data and the Topographic Index Position (TPI) based on the digital elevation model from SRTM (30 m resolution) produced by ALOS (freely available at eorc.jaxa.jp/ALOS/en/about/about_index.htm), and both variables were considered constant from 2005 to 2015 (See Supplementary Material I at the original publication).

basic-services

display_name: basic-services

en: Spatially explicit data to evaluate spatial planning outcomes in a coastal region in São Paulo State, Brazil


en: compression-experiments microstructure snow-mechanics snow x-ray-tomography

en: Snowmicroquakes Data

en: Compression experiments, X-ray images and simulation data with supplemented README file

en: README

media_type: TXT

en: README file explaining the organization of the data

en:
When snow is compressed with a certain speed, micro-snowquakes are triggered in the porous structure of bonded crystals. The present dataset covers uniaxial compression experiments of snow at different strain rates and concurrent X-ray tomography imaging documenting this feature.

The experiments were conducted in a micro-compression stage operated in the X-ray tomography scanner in the SLF cold laboratory. The dataset comprises the compression force data of 17 compression experiments, the 3D image data from 4 X-ray tomography scans and the results of numerical simulations.

compression-experiments

display_name: compression-experiments

microstructure

display_name: microstructure

snow-mechanics

display_name: snow-mechanics

x-ray-tomography

display_name: x-ray-tomography

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: Compressive stick slip and snow-micro-quakes


maintainer: Anita C. Risch

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

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

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

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

format: XLS

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

Paper Citation:

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

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

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

ammonification

display_name: ammonification

climate

display_name: climate

fertilization

display_name: fertilization

global-change

display_name: global-change

grassland

display_name: grassland

herbivore

display_name: herbivore

mineralization

display_name: mineralization

nitrification

display_name: nitrification

nitrogen

display_name: nitrogen

nutrient-network

display_name: nutrient-network

soil

display_name: soil

en: Anthropogenic change and soil net N mineralization