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

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

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

Paper Citation:

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

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

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

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

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

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

bacteria

display_name: bacteria

biogeochemical-cycling

display_name: biogeochemical-cycling

carcass

display_name: carcass

carrion

display_name: carrion

decomposition

display_name: decomposition

detritus

display_name: detritus

fungi

display_name: fungi

soil-biodiversity

display_name: soil-biodiversity

wolf

display_name: wolf

yellowstone

display_name: yellowstone

educ

en: Education, culture and sport

en: Elk and bison carcasses in Yellowstone, USA


maintainer: Henning Löwe

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

state: active

type: dataset

en: Data Structure (Python Notebook - view online)

format: SERVICE

en: Data Structure (Python Notebook)

en: Data Access

media_type: text/plain

format: TXT

mimetype: text/plain

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

discrete-element-modelling

display_name: discrete-element-modelling

microstructure

display_name: microstructure

snow-mechanics

display_name: snow-mechanics

spherical-model-snow

display_name: spherical-model-snow

unconfined-compression-experiments

display_name: unconfined-compression-experiments

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

it: Istruzione, cultura e sport

envidat

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

en: EnviDat: the Environmental Data Portal

it: EnviDat: il portale dei dati ambientali

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

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

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

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

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester

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


photogrammetric-drone-data-dorfberg

maintainer: Yves Bühler

title_for_slug: photogrammetric-drone-data-dorfberg

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

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

All products are also available as a cloud optimized geotiff.

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

drone-data

display_name: drone-data

dsm

display_name: dsm

orthophoto

display_name: orthophoto

photogrammetry

display_name: photogrammetry

snow

display_name: snow

snow-depth-mapping

display_name: snow-depth-mapping

snow-distribution

display_name: snow-distribution

uav

display_name: uav

en: Photogrammetric Drone Data Dorfberg


maintainer: EnviDat Support

en: remote-sensing swe snow-depth snow uas

en: .tfw file

en: coordinates for the .tif raster

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

media_type: image/tiff

en:
10 cm resolution
15 cm accuracy

format: TIFF

mimetype: image/tiff

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

A detailed description is given here:

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

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

remote-sensing

display_name: remote-sensing

snow-depth

display_name: snow-depth

swe

display_name: swe

uas

display_name: uas

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


maintainer: Sofia Pagliarin

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

media_type: text/csv

en: Information on data sources

format: CSV

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

en:
Dataset description

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

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

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

 Dataset composition

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

_IMD15_BCN_MI_Sources.csv_: Information on data sources

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

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

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

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

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

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

 Further information on the Dataset

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

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

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

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

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

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

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

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

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

barcelona

display_name: barcelona

copernicus

display_name: copernicus

imperviousness

display_name: imperviousness

milan

display_name: milan

soil-sealing

display_name: soil-sealing

urban-atlas

display_name: urban-atlas

educ

en: Education, culture and sport

en: Soil sealing Barcelona and Milan different territorial levels


maintainer: Chunhong Zhao

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

en: Land use dynamics in Austin since 1992

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

austin

display_name: austin

land-dynamics

display_name: land-dynamics

land-use

display_name: land-use

landscape

display_name: landscape

urban-morphology

display_name: urban-morphology

envidat

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


maintainer: Anita C. Risch

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

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


maintainer: Anita C Risch

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

en: Risch et al. GEB_data.xlsx

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

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

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

Paper Citation:

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

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

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

anthropogenic-management

display_name: anthropogenic-management

herbivore-exclusion

display_name: herbivore-exclusion

microbial-biomass-carbon

display_name: microbial-biomass-carbon

microbial-respiration

display_name: microbial-respiration

nutrient-addition

display_name: nutrient-addition

soil-properties

display_name: soil-properties


maintainer: Marco Moretti

en: remote-sensing urbanization honeybees sustainability pollination urban-biodiversity

en: Dataset 3

en: Aggregated data on number hives, beekeeping locations and ratios

en: Dataset 1

en: Dataset 1: Processed data for urban beekeeping 2012-2018

en: Dataset 2

en: Dataset 2: Processed data for modelling balance urban greenspace

en: Data on: (1) spatial distribution of urban beekeeping (number of hives and number of beekeeping locations) in 14 Swiss cities (Geneva, Lausanne, Biel, Neuchatel, Basel, Zurich, Chur, Luzern, St. Gallen, Winterthur, Bern, Lugano, Bellinzona, Thun) for the period 2012-2018; (2) aggregated data to model the sustainability of urban beekeeping.

honeybees

display_name: honeybees

pollination

display_name: pollination

remote-sensing

display_name: remote-sensing

sustainability

display_name: sustainability

urban-biodiversity

display_name: urban-biodiversity

urbanization

display_name: urbanization

educ

en: Education, culture and sport

en: Challenging the sustainability of urban beekeeping: evidence from Swiss cities


en: bioenergy review switzerland biogas manure

en: Production de biogaz à partir d’engrais de ferme en Suisse

media_type: application/pdf

en: L'objectif de ce livre blanc est de fournir aux décideurs, aux administrations et aux parties prenantes les résultats de recherche les plus récents afin de promouvoir l'utilisation optimale de la bioénergie issue des engrais de ferme dans la transition énergétique suisse. A cette fin, les résultats du centre de compétence suisse pour la recherche en bioénergie - SCCER BIOSWEET - sont résumés et présentés dans un contexte plus large. Si rien d'autre n'est mentionné, les résultats se réfèrent à la Suisse et, dans le cas de la matière première, aux potentiels nationaux de biomasse.

format: PDF

bioenergy

display_name: bioenergy

biogas

display_name: biogas

manure

display_name: manure

review

display_name: review

switzerland

display_name: switzerland

envidat

en: EnviDat: the Environmental Data Portal

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


en: groundwater switzerland alptal studibach time-series

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

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

alptal

display_name: alptal

groundwater

display_name: groundwater

studibach

display_name: studibach

time-series

display_name: time-series

en: Groundwater time series Studibach (Rinderer et al., 2019, WRR)


en: shortwave-radiation longwave-radiation shading forest fractional-snow

en: Sky-view fraction along the cablecar transect

en: Data is used and described in the following paper: Malle et al. (2019). Shading by Trees and Fractional Snow Cover control the Subcanopy Radiation Budget. Journal of Geophysical Research. DOI:10.1029/2018JD029908.

en: Radiation Data 25.4.2018

en:
Spatiotemporal incoming/outgoing shortwave and longwave radiation as well as sunlit and shaded-snow view fractions below an evergreen forest stand in the Swiss alps. Data was collected on 25th of April 2018 along a 48m long, heterogeneous forest stand under partial snow cover conditions.

The presented data is used and further described in the following paper: Malle et al. (2019). Shading by Trees and Fractional Snow Cover control the Subcanopy Radiation Budget. Journal of Geophysical Research Atmospheres. DOI:10.1029/2018JD029908.

en: Radiation Data 18.04.2019

en:
Spatiotemporal incoming/outgoing shortwave and longwave radiation as well as sunlit and shaded-snow view fractions below an evergreen forest stand in the Swiss alps. Data was collected on 18th of April 2018 along a 48m long, heterogeneous forest stand under full snow cover conditions.

The presented data is used and further described in the following paper: Malle et al. (2019). Shading by Trees and Fractional Snow Cover control the Subcanopy Radiation Budget. Journal of Geophysical Research Atmospheres. DOI:10.1029/2018JD029908.

en: Seasonal sub-canopy radiation data

en:
Seasonal, spatially aggregated incoming and outgoing shortwave radiation below an evergreen forest stand collected with a cable car system.

The presented data is used and further described in the following paper: Malle et al. (2019). Shading by Trees and Fractional Snow Cover Control the Subcanopy Radiation Budget. Journal of Geophysical Research Atmospheres. DOI:10.1029/2018JD029908.

en: This data set consists of incoming and outgoing short- and longwave radiation as well as sunlit-snow-view-fraction as described in the JGR-Atmospheres paper "Shading by trees and fractional snow cover control the sub-canopy radiation budget", by Malle et al. (2019). Data was collected along a 48m long, heterogeneous forest transect between January and June 2018 close to Davos, Switzerland.

forest

display_name: forest

fractional-snow

display_name: fractional-snow

longwave-radiation

display_name: longwave-radiation

shading

display_name: shading

shortwave-radiation

display_name: shortwave-radiation


en: strategic-projects europe land-use-transformations urban-transformations spatial-planning land-use-change strategic-spatial-planning large-scale-urban-development-projects

en: QCA R script large-scale urban development projects Europe

en:
For the set up of this R script, we used the “QCA” (Duşa, 2019) and “SetMethods” (Oana and Schneider, 2018) R packages.

The following further resources were also used:
Dusa, A. (2019). QCA with R: A comprehensive resource. Springer: Berlin.
Schneider, C. Q., & Wagemann, C. (2012). Set-theoretic methods for the social sciences. A guide to qualitative comparative analysis. Cambridge: Cambridge University Press.
Thomann, E., Oana, I.-E., & Wittwer, S. (2018). Performing fuzzy- and crisp set QCA with R. A user-oriented beginner's guide. Retrieved from http://www.evathomann.com/links/qca-r-manual

en: QCA data set large-scale urban development projects

en: Please refer to the sections: 3. The data set lsUDPs: 3.1 Selection of cases and data collection; 3.2 Data management and operationalization in the Metadata for further details on this resource.

en:
Table of Content: 1. General context of the data set "lsUDPs" ; 2. Background and aims of the study using the data set lsUDPs; 3. The data set lsUDPs: 3.1 Selection of cases and data collection; 3.2 Data management and operationalisation

1. General context of the data set "lsUDPs"
The data set "lsUDPs" has been generated as part of the CONCUR research project (https://www.wsl.ch/en/projects/concur.html) led by Dr. Anna M. Hersperger and funded by the Swiss National Science Foundation (ERC TBS Consolidator Grant (ID: BSCGIO 157789) for the period 2016-2020. The CONCUR research project is interdisciplinary and aims to develop a scientific basis for adequately integrating spatial policies (in this case, strategic spatial plans) into quantitative land-change modelling approaches at the urban regional level. The first stage (2016-2017) of the CONCUR project focussed on 21 urban regions in Western Europe. The urban regions were selected through a multi-stage strategy for empirical research (see Hersperger, A. M., Grădinaru, S., Oliveira, E., Pagliarin, S., & Palka, G. (2019). Understanding strategic spatial planning to effectively guide development of urban regions. Cities, 94, 96–105. https://doi.org/10.1016/j.cities.2019.05.032 ).

2. Background and aims of the study using the data set lsUDPs
As part of the CONCUR project, a specific task was to examine the relationship between strategic spatial plans and the formulation and implementation (i.e. urban land change) of large-scale urban development projects in Western Europe.
Strategic urban projects are typically large-scale, prominent urban transformations implemented locally with the aim to stimulate urban growth, for instance in the form of urban renewals of deprived neighborhoods, waterfront renewals and transport infrastructures. While strategic urban projects are referred to in the literature with multiple terms, in the CONCOR project we call them large-scale urban development projects (lsUDPs).
Previous studies acknowledged both local and supra-local (or structural) factors impacting the context-specific implementation of lsUDPs. Local governance factors, such as institutional capacity, coordination among public and private actors and political leadership, intertwine with supra-local conditions, such as state re-scaling processes and devolution of state competencies in spatial planning, de-industrialisation and increasing social inequality. Hence, in implementing lsUDPs, multi-scalar factors act in combination.
Because the formulation and implementation of lsUDPs require multi-scalar coordination among coalitions of public and private actors over an extended period of time, they are generally linked to strategic spatial plans (SSPs). Strategic spatial plans convey collective visions and horizons of action negotiated among public and private actors at the local and/or regional level to steer future urban development, and can contain legally binding dispositions, but also indicative guidelines.
The key question remains as to what extent large-scale urban development projects and strategic spatial plans can be regarded as aligned. By alignment, or “concordance”, we mean that strategic projects are formulated and implemented as part of the strategic planning process (“high concordance”), or that the strategic role of projects is reconfirmed in (subsequent) strategic plans (“moderate concordance”). Lack of concordance is found when lsUDPs have been limitedly (or not at all) acknowledged in strategic spatial plans.
We assume that certain local and supra-local factors, characterising the development of the projects, foster (but not strictly “cause”) the degree of alignment between lsUPDs and SSPs. In this study, we empirically examine how, and to what extent, the concordance between 38 European large-scale urban development projects and strategic plans (outcome: CONCOR) has been enabled by five multi-scalar factors (or conditions): (i) the role of the national state (STATE), (ii) the role of (inter)national private actors (PRIVATE), (iii) the occurrence of supra-regional external events (EVENTS), (iv) the degree of transport connectivity (TRANSP), and (v) local resistance from civil society (RESIST). We adopted a (multi-data) case-based qualitative strategy for empirical research and applied the formalised procedure of within- and cross-case comparison offered by fuzzy-set Qualitative Comparative Analysis appropriate for the goal of this study. Based on set theory, QCA formally integrates contextual sensitivity to case specificities (within-case knowledge) with systematic comparative analysis (across-case knowledge).
The research question the data set has been created to reply to is the following: which conditions, and combinations of conditions, enable the concordance between large-scale urban development projects and strategic spatial plans?
The conditions (“independent variables”) considered are. STATE: the set of large-scale urban projects characterized by a high degree of state intervention and support in their formulation and implementation, PRIVATE: the set of large-scale urban projects characterized by a high degree of involvement of (inter)national private actors in their formulation and implementation, EVENTS: the set of large-scale strategic projects whose formulation and implementation have been strongly affected by unforeseen international events and/or global trends, TRANSP: the set of large-scale strategic projects with a high degree of road and/or transit connectivity, and RESIST: set of large-scale strategic projects whose realization has been characterized by resistances that have substantially delayed or modified the project implementation.
The outcome (“dependent variable”) under analysis is CONCOR: the set of large-scale urban projects having a high degree of concordance/alignment/integration with strategic spatial plans

3. The data set lsUDPs

   3.1 Selection of cases and data collection
To generate the current data set on large-scale urban development projects in European urban regions (data set "lsUDPs"), we identified 35 large-scale urban development projects in a sample of the 21 Western urban regions considered in the CONCUR project (see supra, Hersperger et al. 2019): Amsterdam, Barcelona, Copenhagen, Hamburg, Lyon, Manchester, Milan, Stockholm, Stuttgart. The criteria we followed to identify the 35 large-scale urban development projects are: geographical location, size (large-scale), site (located either in the city core or in the larger urban region) and urban function (e.g. housing, transportation infrastructures, service and knowledge economic functions). Employing these criteria facilitated the selection of diverse large-scale urban development projects while still ensuring sufficient comparability.
In 2016, we performed 47 in-depth interviews with experts in urban and regional planning and large-scale strategic projects and infrastructure (i.e. academics and practitioners) about the formulation, implementation and development (1990s–2010s) of each project in each of the 9 selected urban regions. On average, each interviewee answered questions on 3.1 large-scale urban development projects. Three cases were subdivided into two cases because a clear differentiation between specific implementation stages was identified by the interviewees (expansion of the Barcelona airport, cases “bcn_airport80-90” and “bcn_airport00-16”; realisation of Lyon Part-Dieu, cases “lyo_partdieu70-90” and “lyo_partdieu00-16”; MediaCityUK, cases “man_salfordquays80-00” and “man_mediacityuk00-16”). Therefore, from the initial 35 cases, the final number of analysed cases in the lsUDPs dataset is 38.

    3.2 The data set lsUDPs: Data management and operationalisation
Interviews were fully transcribed and analysed through MAXQDA (version 12.3, VERBI GmbH, Berlin, Germany), and intercoder agreement was evaluated on a sample of nine interviews. We also compiled “synthetic case descriptions” (SCD) for each case (totalling more than 160 SCDs) to spot potential inconsistencies among interviewees’ accounts and to facilitate completion of the “calibration table” for each case (see below). An online expert survey distributed to the interviewees (response rate 78%) helped systematise the information collected during the interviews. We also consulted both academic and gray literature on the case studies to check for possible ambiguity and inconsistencies in the interview data, and to solve discrepancies between our assigned set membership scores and questionnaire values. Site visits were also carried out to retrieve additional information on the selected cases.
For each case (i.e. each of the 38 selected large-scale urban development projects), we operationalised each condition (i.e. STATE, PRIVATE, EVENTS, TRANSP, RESIST) and the outcome (CONCOR) in terms of sets, for subsequent application of Qualitative Comparative Analysis. This process is called “calibration”; we used a number of indicators for each condition to qualitatively assess each large-scale project across the conditions. The case-based qualitative assessment was then transformed into fuzzy-set membership values. Fuzzy-set membership values range from 0 to 1, and should be conceived as “fundamentally interpretative tools” that “operationalize theoretical concepts in a way that enhances the dialogue between ideas and evidence” (Ragin 2000:162, in “Fuzzy-set Social Science”. Chicago: University Press). We employed a four-value fuzzy-set scale (0, 0.33, 0.67, 1) to “quantify” into set membership scores the individual histories of cases retrieved from interview data. Only the condition TRANSP was calibrated as a crisp-set (0, 1).
The translation of qualitative case-based information into numerical fuzzy-set membership values was iteratively performed by populating a calibration table following standard practices recently emerged in QCA when dealing with qualitative (interview) data.

europe

display_name: europe

land-use-change

display_name: land-use-change

land-use-transformations

display_name: land-use-transformations

large-scale-urban-development-projects

display_name: large-scale-urban-development-projects

spatial-planning

display_name: spatial-planning

strategic-projects

display_name: strategic-projects

strategic-spatial-planning

display_name: strategic-spatial-planning

urban-transformations

display_name: urban-transformations

educ

en: Education, culture and sport

en: lsUDPS Large-scale urban development projects in European urban regions


maintainer: Jürg Schweizer

en: avalanche-size avalanche-activity snow-avalanche danger-level avalanche-forcasting

en: Avalanche activity data

en: This data file includes data on the avalanche activity per day aggregated from the avalanche data.

mimetype: text/csv

en: Snow avalanche data

en: This data file includes data on the avalanches mapped in the region of Davos.

en: Readme

en: The Readme file contains the description of the two data files.

mimetype: application/pdf

en:
These data include all avalanches that were mapped in the region of Davos, Switzerland during the winters 1998-1999 to 2018-2019 (21 years), in total 13,918 avalanches, and the corresponding forecast danger level valid on the day of avalanche occurrence, 3533 days and danger ratings in total.

This avalanche activity data set was analysed and results published by Schweizer et al. (2020). They found that the
number of avalanches per day strongly increased with increasing danger level, but avalanche size was poorly related to avalanche danger level.

The data are provided in two files: the first includes the avalanche data (13,918 records); the second includes the avalanche activity per day (3533 records).
Please refer to the Read-me file for further details on the data.

These data are the basis of the following publication:
Schweizer, J., Mitterer, C., Techel, F., Stoffel, A. and Reuter, B., 2020. On the relation between avalanche occurrence and avalanche danger level. The Cryosphere, 14, 737-750, https://doi.org/10.5194/tc-14-737-2020.

avalanche-activity

display_name: avalanche-activity

avalanche-forcasting

display_name: avalanche-forcasting

avalanche-size

display_name: avalanche-size

danger-level

display_name: danger-level

snow-avalanche

display_name: snow-avalanche

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: Snow avalanche data Davos, Switzerland, 1999-2019


en: terrestrial-laser-scans sea-ice-elevation tls polar sea-ice

en: Link to post-processed TLS’s of Antarctic sea ice

media_type: application/zip

format: ZIP

mimetype: application/zip

en: Data processing documentation

en: Post Processed TLS's

en: Compressed dat files

en:
This dataset comprises of a post-processed set of terrestrial laser scans (TLS’s) of Antarctic sea ice obtained during the Sea Ice Physics and Ecosystem Experiment-2 (SIPEX-2, http://seaice.acecrc.org.au/sipex2012/) in September-November 2012. The post-processing steps include the registration of the individual scans into a single 3-dimensional point cloud, the removal of unwanted noise caused by particles in the air (i.e., snow crystals), and the final generation of surface grids based on the cleaned individual point returns. The final product includes the ‘xyz’ coordinates of the individual point measurements, and gridded surfaces covering study areas of 100m x 100 m, and at resolutions of 0.01 m, 0.1 m, 0.25 m, 0.5 m and 1 m for each of the survey dates. Additionally, subgrid statistics that include the mean elevation, standard deviation, minimum and maximum elevations, range, and number of point returns in each gridcell are generated. The final product is provided in space-delimited text files, with the surface grids provided in Digital Terrain Model (DTM) format ready for visualization in any GIS software.

How to cite:
Please also cite the original publication when using this data set.: Trujillo, E., K. Leonard, T. Maksym, and M. Lehning (2016), Changes in snow distribution and surface topography following a snowstorm on Antarctic sea ice, J. Geophys. Res. Earth Surf., 121, doi:[10.1002/2016JF003893](https://dx.doi.org/10.1002/2016JF003893).

polar

display_name: polar

sea-ice

display_name: sea-ice

sea-ice-elevation

display_name: sea-ice-elevation

terrestrial-laser-scans

display_name: terrestrial-laser-scans

tls

display_name: tls

en: High resolution sea ice surface topography from the SIPEX-2 expedition, East Antarctica, 2012


maintainer: Martin Schneebeli

en: rounded-snow snow 3d-model snow-cast weak-layer new-snow depth-hoar 3d-print

en: Description of snowtypes, filenames and other metadata

en: This file contains information about the stl-files of the different snow models. Snowtype (Schneeart), International Snow Classification, preview picture, notes about the origin of the snow, filename of the stl-data, and physical size of the model are given.

en: Snow models in STL-format

en: The zipped file contains all data described in the file 3d_snow_models_metadata. The stl-files should be ready for additive manufacturing.

en: The dataset contains several snow models in the Standard Tesselated Geometry File Format (stl) for 3D visualization, printing and additive manufacturing. Different snow types are available (new snow, rounded snow, depth hoar, buried surface hoar, graupel).

depth-hoar

display_name: depth-hoar

new-snow

display_name: new-snow

rounded-snow

display_name: rounded-snow

snow

display_name: snow

snow-cast

display_name: snow-cast

weak-layer

display_name: weak-layer

educ

en: Education, culture and sport