restoring-grassland-multifunctionality

title_for_slug: restoring-grassland-multifunctionality

en: restoration plants soil-biodiversity bacteria nematodes biodiversity ecosystem-multifunctionality mineralization grassland fungi

en: Functional diversity above_belowground fauna

en: Please note that C.sto and Corg have been corrected.

en: Please note that EMF, MF.below, MF.abio and Soil.hetero have been corrected.

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

Resch, M. C., Schütz, M., Buchmann, N., Frey, B., Graf, U., van der Putten, W. H., Zimmermann, S., Risch, A. C. 2021. Evaluating long-term success in grassland restoration – an ecosystem multifunctionality approach. Ecological Applications 31, e02271.

 Study area

The study was conducted in the Canton of Zurich, Switzerland, in and around two nature reserves Eigental and Altläufe der Glatt (47°27’ to 47°29’ N, 8°37’ to 8°32’ E, 417 to 572 m a.s.l.). All studied grasslands were located with a radius of approximately 4 km. Average monthly temperatures range from 0.7 ± 2.0 °C (January) to 19.0 ± 1.5 °C (July), and monthly precipitation range from 60 ± 42 mm (January) to 118 ± 46 mm (July [maxima]; 1989-2017; MeteoSchweiz 2018). In our study, we focused on semi-dry and semi-wet oligo- to mesotrophic grasslands characterized by high plant species richness and groundwater fluctuations throughout the year (Delarze et al. 2015, see also Resch et al. 2019).

 Experimental design and sampling

A large-scale restoration experiment to expand and reconnect isolated remnants of species-rich grasslands was initiated in the nature reserve Eigental in 1990. Twenty hectares of adjacent intensive grasslands were chosen for restoration. In 1995, three restoration methods of increasing intervention intensities were implemented. The goal of all three methods was to lower the availability of soil nutrients and hence, facilitate ecosystem development towards the targeted nutrient-poor grasslands. These methods were: Harvest only (hay harvest twice a year), Topsoil (removal of the nutrient-rich topsoil), and Topsoil+Propagules (topsoil removal combined with the application of hay from target vegetation). Plant biomass harvest (once a year in late summer/early autumn) commenced in Topsoil and Topsoil+Propagules five years after the soils were removed and is still ongoing today. We measured restoration success by comparing the three restoration methods with intensively managed (Initial) and semi-natural grasslands (Target) 22 years after restoration. Initial grassland sites share the same agricultural history as the restored sites: mowing and subsequent fertilizing (manure) up to five times a year, as well as different tillage regimes (Resch et al. 2019). Target sites were the sites from which hay for seeding the Topsoil+Propagules sites was collected. Soil conditions (i.e., soil types, soil texture) were comparable to those found in the restored grasslands (Resch et al. 2019). Additionally, Target sites were selected to represent a variety of semi-natural grasslands, including semi-dry to semi-wet conditions. In Target grasslands, biomass is harvested once a year in late summer or early autumn. Eleven 5 m x 5 m (25 m2) plots were randomly established in each of the five treatments (in total 55 plots; for a detailed map see Neff et al. 2020). An additional 2 m x 2 m (4 m2) subplot was randomly established at least 2 m away from each 25 m2 plot for destructive sampling. Data sampling took place between June and September 2017.

Vegetation properties

All plant species were identified within the 25 m2 plots (nomenclature: Lauber and Wagner 1996) in mid-June 2017 (in total 250 species). Vegetation structure and plant biomass were assessed diagonally on a transect of 2 m x 10 cm within the 25 m2 plot in early July 2017. We measured the maximum and mean height of the vegetation at the start, middle and end of the transect and calculated the standard deviation of these measures to describe vegetation structural heterogeneity (Schuldt et al. 2019). Thereafter, biomass was clipped on the entire transect to 1 cm height, sorted into five functional groups (graminoids, forbs, legumes, litter, and woody species), dried at 60 °C for 48 h, and weighed (Meyer et al. 2015).

 Aboveground arthropods

Aboveground arthropods were sampled at two locations in each 25 m2 plot in early July 2017 (see also Neff et al. 2020). Briefly, two cylindrical baskets (50 cm diameter, 67 cm height; woven fabric) were thrown simultaneously from outside the plot into two opposite corners. A closable mosquito mesh sleeve was mounted to the top of the baskets and an integrated metal ring at the bottom was fixed to the ground with metal stakes to assure that insects could not escape. A suction sampler (Vortis, Burkhard Manufacturing Co. Ltd., Hertfordshire, England) was then inserted into one of the baskets through the opening of the sleeve and the plot was “vacuumed" twice for 105 seconds with a 30 seconds break. The collected animals were immediately transferred into 70% ethanol. Arthropods were sorted and assigned to 23 taxonomic groups. Holometabolic larvae were lumped into one category while hemimetabolic larvae were grouped separately from adults in the respective taxonomic rank. We used mean values of individuals per plot for total abundance. Aboveground arthropod richness was defined by the number of different taxa to lowest taxonomic level (in total 23 taxa). All taxa were assigned to one of five trophic levels: 1) primary producers, 2) primary consumers, 3) secondary consumers, 4) tertiary consumers, and 5) quaternary consumers.

 Belowground fauna

Sampling of all belowground fauna took place in mid-July 2017. Earthworms were sampled in two 30 cm x 30 cm x 20 cm soil monoliths at two opposite corners of the 25 m2 plot (opposite to aboveground arthropod sampling). The excavated soil monolith was broken by hand, all earthworms collected and immediately transferred in a 4% formaldehyde solution. Thereafter, earthworm individuals were identified to species level (in total 10 taxa; Christian and Zicsi 1999) and species assigned to three functional groups (Bouché 1977).
To assess soil arthropod communities, we randomly collected one undisturbed soil core (5 cm diameter, 12 cm depth) in each 4 m2 subplot with a slide hammer corer lined with a plastic sleeve (AMS Samplers, American Falls, Idaho, USA). Soil arthropods were extracted using Berlese-Tullgren funnels (3 mm mesh), starting the day of sampling and lasting 14 days. Individuals were stored in 70% ethanol. Soil arthropods were assigned to 41 taxonomic groups and 4 feeding types. Holometabolic and hemimetabolic larvae were treated as previously described for aboveground arthropods. Belowground arthropod richness refers to the 41 taxonomic groups.

For soil nematode sampling, we randomly collected eight soil cores of 2.2 cm diameter (Giddings Machine Company, Windsor, CO, USA) within each 4 m2 subplot to a depth of 12 cm. The eight cores were combined, gently homogenized, placed in coolers, kept at 4 °C and transported to the laboratory at NIOO in Wageningen (NL) within one week after collection. Free-living nematodes were extracted from 200 g of fresh soil using Oostenbrink elutriator (Oostenbrink 1960) and prepared for morphological identification and quantification as described by Resch et al. (2019). Nematodes were identified to family level (39 taxa) according to Bongers (1988), assigned to 17 functional groups, 5 feeding types and 5 colonizer-persister (C-P) classes (Yeates et al. 1993, Bongers 1990, Resch et al. 2019).

We randomly collected two more soil cores (2.2 cm diameter x 12 cm depth) within each 4 m2 subplot to determine soil microbial communities. Again, the soil cores were combined, homogenized, placed in coolers and transported to the laboratory at WSL in Birmensdorf (Switzerland) where the metagenomic DNA was extracted from 8 g sieved soil (2 mm) using the DNeasy PowerMax Soil Kit (Quiagen, Hilden, NRW, GER) according to the manufacturer`s instructions. PCR amplification of the V3-V4 region of the prokaryotic small-subunit (16S) and the ribosomal internal transcribed spacer region (ITS2) of eukaryotes was performed with 1 ng of template DNA utilizing PCR primers and conditions as previously described (Frey et al., 2016). PCRs were run in triplicates and pooled. The pooled amplicons were sent to the Genome Quebec Innovation Centre (Montreal, QC, Canada) for barcoding using the Fluidigm Access Array technology (Fluidigm) and paired-end sequencing on the Illumina MiSeq v3 platform (Illumina Inc., San Diego, CA, USA). Quality filtering, clustering into operational taxonomic units (OTUs) and taxonomic assignment were performed as described by Frey et al. (2016) and Adamczyk et al. (2019). We used a customised pipeline largely based on UPARSE (Edgar 2013) implemented in USEARCH v. 9.2 (Edgar 2010). After discarding singletons of dereplicated sequences, clustering into OTUs with 97% sequence similarity was performed (Edgar 2013). Quality-filtered reads were mapped on the filtered set of centroid sequences. Taxonomic classification of prokaryotic and fungal sequences was conducted querying against most recent versions of SILVA (v.132, Quast et al. 2013) and UNITE (v.8, Nilsson et al. 2018). Only taxonomic assignments with confidence rankings equal or higher than 0.8 were accepted (assignments below 0.8 set to unclassified). Prokaryotic OTUs assigned to mitochondria or chloroplasts as well as eukaryotic OTUs assigned other than fungi were removed prior to data analysis. In addition, prokaryotic and fungal datasets were filtered to discard singletons and doubletons. Thereafter, OTU abundance matrices were rarefied to the lowest number of sequences per community, to normalize the total number of reads and achieve parity between samples (Prokaryota: 29,843 reads; Fungi: 26,690 reads). Finally, prokaryotic and fungal observed richness (number of OTUs) were estimated (Prokaryota: 14,010 OTUs; Fungi: 5,813 OTUs). For prokaryotes, we distinguished five and for fungi six functional types based on lowest taxonomic resolution (Nguyen et al. 2016, Tedersoo et al. 2014). Belowground taxon richness was defined by the total number of earthworm, arthropod, nematode, fungi, and prokaryote taxa assigned to lowest taxonomic level. Finally, all belowground taxa were assigned to the same five trophic levels as the aboveground arthropods.

 Soil chemical and physical properties, soil nitrogen mineralization

We randomly collected three 5 cm diameter x 12 cm depth soil samples in each 4 m2 subplot with a slide hammer corer (AMS Samplers, American Falls, Idaho, USA), pooled them and then made two subsamples. One was field-fresh and stored at 3 °C until analysis, the other was dried for 48 h at 60 °C and passed through a 4 mm mesh. From the dried sample, we measured soil pH potentiometrically in 0.01 M CaCl2 (soil:solution ratio=1:2; 30 minutes equilibration time). Total and organic carbon content were measured on fine-ground samples (≤ 0.5 mm) by dry combustion using a CN analyzer NC 2500 (CE Instruments, Wigan, United Kingdom). Inorganic carbon of samples with a pH > 6.5 was removed with acid vapor prior to analysis of organic carbon (Walthert et al. 2010). We calculated total soil carbon (C) storage after correcting its content for soil depth, stone content and density of fine earth (see below). Exchangeable cations were determined on another 5 g dry soil sample with 50 mL unbuffered 1 M NH4Cl solution (soil:solution ratio=1:10, end-over-end shaker for 1.5 hours) and measured by an ICP-OES (Optima 7300 DV, Perkin-Elmer, Waltham, Massachusetts, USA). Thereafter, cation exchange capacity (CEC) was calculated as the sum of exchangeable cations and protons (and expressed as mmolc per 1 kg soil) and used to describe nutrient retention capacity in our plots. Concentrations of exchangeable protons were calculated as the difference between total and Al-induced exchangeable acidity as determined by the KCl-method (Thomas 1982).

Ammonium (NH4+) and nitrate (NO3−) were extracted from a 20 g fresh subsample with 80 mL 1M KCl for 1.5 hours on an end-over-end shaker and filtered through ashless folded filter paper (DF 5895 150, ALBET LabScience, Hahnemühle FineArt GmbH, Dassel, Germany). NH4+ concentrations were determined colorimetrically by automated flow injection analysis (FIAS 300, Perkin-Elmer, Waltham, Massachusetts, USA). NO3− concentrations were measured colorimetrically according to Norman and Stucki (1981). Potential soil net nitrogen (N) mineralization was assessed during an 8-week incubation period under controlled moisture (60% of field capacity), temperature (20 °C) and light conditions (dark) in the laboratory. We weighed duplicate samples of fresh soil equivalent to 8 g dry soil (24 h at 104 °C) into 50 mL Falcon tubes. Soil samples were extracted for NH4+ and NO3− at the beginning and after eight weeks as described above. Soil net N mineralization was calculated as the difference between the inorganic nitrogen (NH4+ and NO3−) before and after the incubation (Hart et al. 1994), corrected for the total incubation time and represented per day values expressed as mg N kg-1 soil d-1.

To assess soil physical properties, we randomly collected one undisturbed soil core per 4 m2 subplot (5 cm diameter, 12 cm depth) in a steel cylinder that fit into the slide hammer (AMS Samplers, American Falls, Idaho, USA). The cylinder was capped in the field to avoid disturbance. We then measured field capacity in the laboratory. For this purpose, the cylinder and soil therein were saturated in a water bath and drained on a sand/silt-bed with a suction corresponding to 60 cm hydrostatic head. The moist soil was dried at 105 °C to constant weight. Field capacity was calculated by dividing the mass of water by the total mass of wet soil contained at 60 cm hydrostatic head and used to describe water holding capacity. Thereafter, samples were passed through a 4 mm mesh. Fine-earth and skeleton fractions were weighed separately to assess bulk soil density (fine-earth plus skeleton), density of fine earth, and proportion of skeleton. Particle density was determined with the pycnometer method (Blake and Hartge 1986), and total porosity and proportion of fine pores were calculated (Danielson and Sutherland 1986). Clay, silt, and sand contents were quantified with the sediment method (Gee and Bauder 1986).

Surface and soil temperature (12 cm depth, water-resistant digital pocket thermometer; IP65, H-B Instrument, Trappe, Pennsylvania, USA) as well as volumetric soil moisture content (12 cm depth, time domain reflectometry; Field-Scout TDR 300, Spectrum Technologies, Aurora, Illinois, USA) were measured at five random locations within the 4 m2 subplots every month from June to September. We calculated the standard deviation of each temperature and moisture measure over four months to describe seasonal variations. Slope inclination was determined at plot-level via GPS measurements (GPS 1200, Leica Geosystem, Heerbrugg, Switzerland) and categorized into slope gradient classes according to FAO standards (1990). Thickness of the topsoil horizon (equivalent to Ah or Aa horizon) was determined at one soil monolith (30 x 30 x 30 cm3) per 4 m2 subplot in cm and rounded to next integer.

 References

Adamczyk, M., F. Hagedorn, S. Wipf, J. Donhauser, P. Vittoz, C. Rixen, A. Frossard, J. Theurillat, and B. Frey. 2019. The soil microbiome of GLORIA mountain summits in the Swiss Alps. Frontiers in Microbiology 10:1080-1101.

Blake, G.R., and K. H. Hartge. 1986. Particle Density. Pages 377-382 in A. Klute, editor. Methods of soil analysis: Part 1—Physical and mineralogical methods. Soil Science Society of America (SSSA) Inc., Madison.

Bongers, T. 1988. De nematoden van Nederland. Stichting Uitgeverij van de Koniklijke Nederlandse Natuurhistorische Verenigung (KNNV), Utrecht.

Bongers, T. 1990. The maturity index: an ecological measure of environmental disturbance based on nematode species composition. Oecologia 83:14-19.

Bouché, M. B. 1977. Strategies lombriciennes. Ecological Bulletins 25:122-132.

Christian, E., and A. Zicsi. 1999. Ein synoptischer Bestimmungsschlüssel der Regenwürmer Österreichs (Oligochaeta: Lumbricidae). Die Bodenkultur 50:121-131.

Danielson, R. E., and P. L. Sutherland. 1986. Porosity. Pages 443-461 in A. Klute, editor. Methods of soil analysis: Part 1—Physical and mineralogical methods. Soil Science Society of America (SSSA) Inc., Madison.

Delarze, R., Y. Gonseth, S. Eggenberg, and M. Vust. 2015. Lebensräume der Schweiz: Ökologie ‐ Gefährdung ‐ Kennarten, Ott Verlag, Bern.

Edgar, R. C. 2010. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26:2460–2461.

Edgar, R. C. 2013. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nature Methods 10:996–998.

FAO. 1990. Guidelines for soil description, third ed. Land and Water Development Division at the Food and Agriculture Organization of the United Nations (FAO), Rome.

Frey, B., T. Rime, M. Phillips, B. Stierli, I. Hajdas, F. Widmer, and M. Hartmann. 2016. Microbial diversity in European alpine permafrost and active layers. FEMS Microbiology Ecology 92:fiw018.

Gee, G.W., and J. W. Bauder. 1986. Particle-size analysis. Pages 383-411 in A. Klute, editor. Methods of soil analysis: Part 1—Physical and mineralogical methods. Soil Science Society of America (SSSA) Inc., Madison.

Hart, S. C, J. M. Stark, E. A. Davidson, and M. K. Firestone. 1994. Nitrogen mineralization, immobilization, and nitrification. Pages 985-1016 in R. W. Weaver, S. Angle, P. Bottomley, D. Bezdicek, S. Smith, A. Tabatabai, and A. Wollum, editors. Methods of soil analysis: Part 2—Microbiological and biochemical properties. Soil Science Society of America (SSSA) Inc., Madison.

Lauber, K., Wagner, G., 1996. Flora Helvetica. Flora der Schweiz. Haupt Verlag, Bern.

MeteoSchweiz, 2018. Klimabulletin Jahr 2017. MeteoSchweiz, Zürich.

Meyer, S. T., C. Koch, and W. W. Weisser. 2015. Towards a standardized rapid ecosystem function assessment (REFA). Trends in Ecology and Evolution 30:390-397.

Neff, F., M. C. Resch, A. Marty, J. Rolley, M. Schütz, A. C. Risch, and M. M. Gossner. 2020. Long-term restoration success of insect herbivore communities in semi-natural grasslands – a functional approach. Ecological Applications 0:e02133

Nguyen, N.H., Z. W. Song, S. T. Bates, S. Branco, L. Tedersoo, J. Menke, J. S. Schilling, and P. G. Kennedy. 2016. FUNGuild: an open annotation tool for parsing fungal community datasets by ecological guild. Fungal Ecology 20:241–248.

Nilsson, R. H., K.-H. Larsson, A. F. S. Taylor, J. Bengtsson-Palme, T. S. Jeppesen, D. Schigel, P. G. Kennedy, K. Picard, F. O. Glöckner, L. Tedersoo, I. Saar, U. Kõljalg, and K. Abarenkov. 2018. The UNITE database for molecular identification of fungi: handling dark taxa and parallel taxonomic classifications. Nucleic Acids Research 47:259–264.

Norman, R. J., and J. W. Stucki. 1981. The Determination of Nitrate and Nitrite in Soil Extracts by Ultraviolet Spectrophotometry. Soil Science Society of America Journal 45:347-353.

Oostenbrink, M. 1960. Estimating nematode populations by some selected methods. Pages 81-101 in J. J. Sasser, and W. R. Jenkins, editors. Nematology. Univ. of North Carolina Press, Chapel Hill.

Quast, C., E. Pruesse, P. Yilmaz, J. Gerken, T. Schweer, P. Yarza, J. Peplies, and F. O. Glöckner. 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Research 41:590-596.

Resch, M. C., M. Schütz, U. Graf, R. Wagenaar, W.H. van der Putten, and A. C. Risch. 2019. Does topsoil removal in grassland restoration benefit both soil nematode and plant communities? Journal of Applied Ecology 56:1782-1793.

Schuldt, A., A. Ebeling, M. Kunz, M. Staab, C. Guimarães-Steinicke, D. Bachmann, N. Buchmann, W. Durka, A. Fichtner, F. Fornoff, W. Härdtle, L. R. Hertzog, A-N. Klein, C. Roscher, J. Schaller, von G. Oheimb, A. Weigelt, W. Weisser, C.Wirth, J. Zhang, H. Bruelheide, and N. Eisenhauer. 2019. Multiple plant diversity components drive consumer communities across ecosystems. Nature Communications 10:1460.

Tedersoo, L., M. Bahram, S. Põlme, U. Kõljalg, N. S. Yorou, R. Wijesundera, L. Villarreal Ruiz, A. M. Vasco-Palacios, P. Q. Thu, A. Suija, M. E. Smith, C. Sharp, E. Saluveer, A. Saitta, M. Rosas, T. Riit, D. Ratkowsky, K. Pritsch, K. Põldmaa, M. Piepenbring, C. Phosri, M. Peterson, K. Parts, K. Pärtel, E. Otsing, E. Nouhra, A. L. Njouonkou, R. H. Nilsson, L. N. Morgado, J. Mayor, T. W. May, L. Majuakim, D. J. Lodge, S. See Lee, K.-H. Larsson, P. Kohout, K. Hosaka, I. Hiiesalu, T. W. Henkel, H. Harend, L.-D. Guo, A. Greslebin, G. Grelet, J. Geml, G. Gates, W. Dunstan, C. Dunk, R. Frenkhan, L. Dearnaley, A. De Kesel, T. Dang, X. Chen, F. Buegger, F. Q. Brearley, G. Bonito, S. Anslan, S. Abell, and K. Abarenkov. 2014. Global diversity and geography of soil fungi. Science 346:1256688.

Thomas, G.W. 1982. Exchangeable cations. Pages 159-165 in A. L. Page, R. H. Miller, and D. R. Keeney, editors. Methods of Soil Analysis: Part 2—Chemical and microbiological properties. Soil Science Society of America (SSSA) Inc., Madison.

Walthert, L., U. Graf, A. Kammer, J. Luster, D. Pezzotta, S. Zimmermann, and F. Hagedorn. 2010. Determination of organic and inorganic carbon, δ13C, and nitrogen in soils containing carbonates after acid fumigation with HCl. Journal of Plant Nutrition and Soil Sciences 173:207-216.

Yeates, G. W., T. Bongers, R. G. M. de Goede, D. W. Freckman, and S. S. Georgieva. 1993. Feeding habits in soil nematode families and genera – an outline for soil ecologists. Journal of Nematology 25:315-331.

bacteria

display_name: bacteria

ecosystem-multifunctionality

display_name: ecosystem-multifunctionality

fungi

display_name: fungi

grassland

display_name: grassland

mineralization

display_name: mineralization

nematodes

display_name: nematodes

plants

display_name: plants

restoration

display_name: restoration

soil-biodiversity

display_name: soil-biodiversity

en: Restoring grassland multifunctionality


maintainer: EnviDat Support

en: density buffer jupyter gis roads geospatial

state: active

type: dataset

en: Road Density Calculation (Python Notebook) online

format: SERVICE

media_type: application/zip

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

format: ZIP

mimetype: application/zip

en: Plots

en: Random plot locations

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

buffer

display_name: buffer

density

display_name: density

geospatial

display_name: geospatial

gis

display_name: gis

jupyter

display_name: jupyter

roads

display_name: roads

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

en: Education, culture and 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: Dataset for OGRS 2018 publication


en: backward-trajectories clouds raclets ice-nucleation aerosol

media_type: text/csv

format: CSV

mimetype: text/csv

en: This dataset contains number concentrations of ice-nucleating particles active at -15 °C observed at Weissfluhjoch during February and March 2019, as well as complementary data (measured aerosol number concentrations and modelled total precipitation along air mass trajectories). This data formed the basis of our paper with the title “Towards parameterising atmospheric concentrations of ice-nucleating particles active at moderate supercooling”.

aerosol

display_name: aerosol

backward-trajectories

display_name: backward-trajectories

clouds

display_name: clouds

ice-nucleation

display_name: ice-nucleation

raclets

display_name: raclets

en: Ice nucleating particle concentrations active at -15 °C at Weissfluhjoch


en: spatial-modelling avalanche-forecast avalanche-release-probability machine-learning forecast-models

en: Model predictions (interpolated)

en: The zip file contains csv-files the interpolated model predictions for three models, in forecast and nowcast-mode, for the three different data sets, as described in the publication.

en: GPS data and human-triggered avalanches

en: Contains (1) the subset of points extracted from GPS tracks recorded during backcountry touring and uploaded to the website skitourenguru.ch and (2) the corresponding subset of human-triggered avalanches. Data (1) was provided by Skitourenguru GmbH and data (2) was collected by SLF. - Variables contained are date, Swiss coordinates (x, y), aspect (N = 0, NE = 45, ...), elevation (ele), and event (event = 1, non-event = 0).

format: XLS

en: Human-triggered avalanches 2022/2023

en: The file contains the reported human-triggered avalanches. The variables are as follows: id, release date (date), name, Swiss coordinates (x,y), elevation (ele), remark, reported avalanche size expressed using avalanche activity index (sizeAAI), aspect, event type (type), event (yes = 1, no = 0)

en: Natural avalanches 2023/2024

en: Natural avalanche observation (dry avalanches). 2023/2024 forecasting season.

en: Randomly sampled subset of grid points from 1 km DEM

en: Contains the randomly sampled subset of grid points from 1 km DEM. Three variables: coordinates (x, y) and elevation (elev).

en: Snowline estimates 2023/2024

en: Daily snowline estimates provided by observers for North (aspect = 0) and South (aspect = 180) slopes for the 2023/2024 season. For aspects in between North and South, values are interpolated. Contains observation date (date), coordinates (Swiss coordinates x, y), snowline elevation (snowline), and aspect.

en: Model predictions (raw data)

en: Model predictions (raw data) for three models, in nowcast and forecast mode. Not all models are available for both seasons. - These data are the basis for the interpolated predictions, which can be found in the other zip-folder

en: Avalanche forecast

en: Data contains forecasts relating to dry-snow avalanche conditions only. Forecast issued at 17.00 local time. Contains variables for the warning region (sector_id), publication date (date), danger level - sub-level combination expressed using numeric values (integer value of danger level -0.33 for sub-level minus, -0 for sub-level neutral, +0.33 for sub-level plus, altitude as indicated in the bulletin, and whether aspect (North = 0, Northeast = 45, ...) was indicated in the bulletin (no: 0, yes: 1).

en: Snowline estimates 2022/2023

en: Daily snowline estimates provided by observers for North (aspect = 0) and South (aspect = 180) slopes for the 2022/2023 season. For aspects in between North and South, values are interpolated. Contains observation date (date), coordinates (Swiss coordinates x, y), snowline elevation (snowline), and aspect.

en: Natural avalanches 2022/2023

en: Natural avalanche observation (dry avalanches). 2022-2023 forecasting season.

en: Human-triggered avalanches 2023/2024

en:
This data set was used in the analysis by Techel et al. *Forecasting avalanche danger: human-made forecasts vs. fully automated model-driven predictions*, submitted to *Natural Hazards Earth System Sciences* on 20 Aug 2024.
The repository contains data from two avalanche forecasting seasons (2022/2023, 2023/2024) in Switzerland.
**Interpolated predictions**
- The .zip file contains the interpolated predictions for the three models in nowcast- and forecast- mode. This data is needed to reproduce the figures and tables in the submitted preprint.
The other data are the **raw data** underlying the interpolations:
- Avalanche forecast by WSL Institute for Snow and Avalanche Research SLF, published at 17.00 local time, valid for the following 24 hours and relating to dry snow avalanche conditions.
- Model predictions in *nowcast*- and *forecast*-mode for three models (*danger level*, *instability*, *natural avalanche*), valid for 12.00 local time
- Subset of points extracted from GPS tracks (courtesy of Skitourenguru GmbH)
- Avalanche observations - natural avalanches and human-triggered avalanches
- Estimates of the snowline
- Randomly chosen subset of grid points used for generating reference distributions
For details regarding the data sets refer to the publication.

avalanche-forecast

display_name: avalanche-forecast

avalanche-release-probability

display_name: avalanche-release-probability

forecast-models

display_name: forecast-models

machine-learning

display_name: machine-learning

spatial-modelling

display_name: spatial-modelling

en: Evaluating the predictive performance of human avalanche forecasts and model predictions in Switzerland


maintainer: Anita C. Risch

en: non-native-plant herbivory weed-management grassland native-plant

en: Correction: Data file native and non-native plant interactions in Australian grasslands

en:
taxa tab:	list and overview of all the recorded species

* family:                 family of the plant

* genus:                genus of the plant

* species:                plant species name

* origin:                indicates whether the plant is native to Australia or not (non-native)

* common name:	        common name of the plant species

* functional group:	        functional group of the plant species: forb, graminoid, fern, tree

* functional group detailed:	gives more detail about the type of graminoid and forb: legume, forb (non-leguminous),
                        grass, sedge, rush, fern, tree

* lifeform:                annual or perennial lifeform (species which can have both of these lifeforms where
                        considered as annuals)

plot tab:	list of species found in each plot and their cover (%)

* sampling.session:        provides information about the sampling time:
                        May2020 = data that was collected in May 2020
                        Nov2020 = data that was collected in November 2020

* farm.name:	            name of the farms where the sampling took place:
                        SF = farm of the Southern group
                        NF = farm of the Northern group
                        numbers 1-3 = randomly allocated to the Southern group of farms
                        numbers 4-6 = randomly allocated to the Northern group of farms

* site.type:	                KG: Kangaroo grass dominated site
                        ALG: ALG+KG site (Kangarooo grass and African lovegrass co-dominated site)

* plot.no:                plot number: 9-12

* scientific.name:	        scientific name of the plant species found in a certain plot

* functional group:	        indicates whether the species found is a tree, forb, graminoid or fern

* functional group detailed:	indicates whether the species found is a tree, legume, forb (non-leguminous),
                        grass, sedge, rush, fern

* lifeform:                specifies the lifeform of the species: annual or perennial lifeform
                        (species which can have both of these lifeforms where considered as annuals)

corrections tab:	Detailed description of corrections made in plot.data tab

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

Paper Citation:

_Schlierenzauer, C., Risch, A.C., Schütz, M., Firn, J. 2021. Non-native Eragrostis curvula reduces plant species diversity in pastures of South-eastern Australia even when native Themeda triandra remains co-dominant. Plants 10, 596._

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

Study area
The study was conducted in the lowland grassy woodlands of the Bega Valley Region, which is located in the south-east corner of New South Wales, Australia. Embedded between the Pacific Ocean and the Australian Alps, the lowland grassy woodlands are mostly located on granitic substrates and reach elevations of roughly 500 m above sea level. Typically, these grassy woodlands receive less precipitation (mean annual precipitation between 700-1100 mm) compared to the more elevated areas that surround them (NSW Government - Office of Environment and Heritage 2017).
The vegetation is dominated by an open tree canopy layer consisting of Eucalyptus tereticornis Sm, Angophora floribunda Sm. (Sweet) and a range of other eucalypt species. Sometimes shrub or small trees are also present, whereas grasses and forbs form the ground-cover. In areas without intensive agricultural history, this layer is dominated by perennial, tussock grasses such as Themeda triandra Forssk, Microlaena stipoides R.Br (Weeping Grass), Eragrostis leptostachya Steud. (Paddock Lovegrass) and Echinopogon ovatus P.Beauv (Forest Hedgehog Grass). The remaining inter-tussock spaces are occupied by a diversity of growth-restricted grasses and herbaceous forbs (NSW - Department of Planing, Industries and Environment 2019; NSW Government - Office of Environment and Heritage 2017).
Clearing, pasture sowing, fertilizer application and livestock grazing resulted in a dramatic decrease in the extent of these natural woodlands, with less than five percent within conservation reserves and overall, with only about 20% of their original extent in New South Wales still existing (Tozer et al. 2010). The remaining areas outside of reserves are threatened by altered fire frequencies, habitat clearing, livestock grazing and especially by non-native plant invasion, particularly Eragrostis curvula (Schrad.) Nees. For this reason, the grassy woodlands are listed as an endangered ecological community in the NSW state legislation. Additionally, they are considered as critically endangered by the Commonwealth of Australia (Threatened Species Scientific Committee (TSSC) 2013).

Experimental design and sampling
The study was conducted on six farms and in each of them two sites were chosen, representing a paired design. One of the sites at each farm is dominated by native Themeda triandra, the other one co-dominated by non-native Eragrostis curvula and Themeda triandra. All farms are within a radius of approximately 10 km from the town Candelo. Three of the farms are located North (36°40’ to 36°42’ S and 149°38’ to 149°42’ E) and three of them are located South (36°51’ to 36°49’ S and 149°38’ to 149°42’ E) of Candelo.
Non-native herbivores (mainly cattle, sheep and rabbits) and native herbivorous marsupials (mainly kangaroos, wallabies and wombats) are present in the area of these sites. On each site, data was collected within four plots (each 1 x 1 m) in May and November 2020. All plant species found within a plot were recorded and their relative abundance was estimated.

References
NSW - Department of Planing, Industries and Environment. 2019. “Lowland Grassy Woodland in the South East Corner Bioregion - Endangered Ecological Community Listing.” https://www.environment.nsw.gov.au/topics/animals-and-plants/threatened-species/nsw-threatened-species-scientific-committee/determinations/final-determinations/2004-2007/lowland-grassy-woodland-south-east-corner-bioregion-endangered-ecological-community-l (February 18, 2021).
NSW Government - Office of Environment and Heritage. 2017. “Lowland Grassy Woodland in the South East Corner Bioregion - Profile.” https://www.environment.nsw.gov.au/threatenedSpeciesApp/profile.aspx?id=20070 (January 31, 2021).
Threatened Species Scientific Committee (TSSC). 2013. Environment Protection and Biodiversity Conservation Act 1999 (EPBC Act) Conservation Advice for Lowland Grassy Woodland in the South East Corner Bioregion. http://www.environment.gov.au/biodiversity/threatened/communities/pubs/82-conservation-advice.pdf.
Tozer, Mark et al. 2010. “Native Vegetation of Southeast NSW: A Revised Classification and Map for the Coast and Eastern Tablelands.” Cunninghamia : a journal of plant ecology for eastern Australia 11(3): 359–406.

grassland

display_name: grassland

herbivory

display_name: herbivory

native-plant

display_name: native-plant

non-native-plant

display_name: non-native-plant

weed-management

display_name: weed-management

educ

en: Education, culture and sport

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: Native and no-native plant interactions in Australian grasslands


maintainer: Giulia Francesca Azzurra Donati

en: ecosystem-services urban-planning network-analysis ecological-corridors biodiversity circuit-theory

en: Amphibian movement resistance layers and scripts for connectivity modeling - part 1

en: The zipped folder includes the script to aggregate the resistance layers into squared and hexagonal grids. The subfolders “Boundaries”, “Layers”, “Reference_grids”, and “Resistance_from_SDMs” contain files either called by the code or created in the process.The subfolder “Inputs_for_Circuitscape” contains the inputs to run Circuitscape, after doing some subsequent processing steps in QGIS (listed in the R script) with the files created with the script.

en: Amphibian movement resistance layers and scripts for connectivity modelling - part 2

en: The zipped folder includes the script to create and analyse the networks from Circuitscape outputs. The subfolders “Boundaries”, “Layers”, and “Reference_net” contain files either called by the code or created in the process. The subfolder “Outputs_from_Circuitscape” contains the outputs from Circuitscape including node and branch cumulative currents in a network format.

en: This data includes (1) Scripts to aggregate landscape resistance layers into squared and hexagonal grids (i.e., different representations and resolutions), (2) Input resistance layers and focal nodes in .txt format to run in Circuitscape (Python implementation v4.0.5). Circuitscape is a software tool for modeling and analyzing landscape connectivity, which simulates movement of organisms across landscapes by estimating resistance to movement across each point of the landscape. (3) Scripts for the ecological network analysis, and (4) environmental predictors for amphibian whole-life cycle habitats used to describe the local environment for BGI design (i.e. topographic, hydrologic, edaphic, vegetation, land-use derived, movement-ecology related).

biodiversity

display_name: biodiversity

circuit-theory

display_name: circuit-theory

ecological-corridors

display_name: ecological-corridors

ecosystem-services

display_name: ecosystem-services

network-analysis

display_name: network-analysis

urban-planning

display_name: urban-planning

en: Multi-Scale Prioritization framework for Urban Blue-Green Infrastructure Planning to Support Biodiversity: Data & Codes


maintainer: Thorsten Bartels-Rausch

en: alpine metamorphism brine ozone chemistry interface laboratory-study reaction arctic snow solution antarctica heterogeneous

en: Matlab Script to view and plot data

en: Data Figure 1

en: Data of published Figure

en: Readme file

media_type: application/pdf

en: This file shows data, code, and plots.

format: PDF

mimetype: application/pdf

en: Matlab Data

en: All data as Matlab mat file

en: data list

en: Lists file names, snow sample properties, and flow settings for each data-set.

en: Ozone data

en: Trend of ozone mixing ratio ( _O3mixingRatio_ppb_ with time ( _time_s_ ) as measured downstream of the snow sample during the experiments. Time "0" denotes the time when data acquisition with ozone analyser was started. And, ozone (_nOzone_molec_per_s_ ) loss for each experiment with duration of ozone exposure ( _duration_s_ ) of ozone exposure, as derived from the ozone trends shown above. Time "0" denotes the time when the carrier gas with ozone was passed via the snow sample.

en: Acetone data

en:
intensity of the acetone-water-cluster and water-cluster (as derived by a chemical ionisation mass spectrometer) with time in the carrier gas downstream of packed bed flow tube holding the snow sample or the by-pass, respectively. Time "0" denotes the time when the mass spectrometer data acquisition was started. When the carrier
gas with acetone was passed via the packed bed flow tube, that is switched from by-pass to pass the snow, acetone concentration in the carrier gas rapidly dropped due to surface adsorption and slowly recovered as the adsorption sites saturate. From this, the surface area can be obtained.

en: Earth’s snow cover is very dynamic on diurnal time scales. The changes to the snow structure during this metamorphism have wide ranging impacts such as on avalanche formation and on the capacity of surface snow to exchange trace gases with the atmosphere. Here, we investigate the influence of dry metamorphism, which involves fluxes of water vapor, on the chemical reactivity of bromide in the snow. For this, the heterogeneous reactive loss of ozone at a concentration of 5-6E12 molecules cm-3 is investigated in artificial, shock-frozen snow samples doped with 6.2 uM sodium bromide and with varying metamorphism history. The oxidation of bromide in snow is one reaction initiating polar bromine releases and ozone depletions.

alpine

display_name: alpine

antarctica

display_name: antarctica

arctic

display_name: arctic

brine

display_name: brine

chemistry

display_name: chemistry

heterogeneous

display_name: heterogeneous

interface

display_name: interface

laboratory-study

display_name: laboratory-study

metamorphism

display_name: metamorphism

ozone

display_name: ozone

reaction

display_name: reaction

snow

display_name: snow

solution

display_name: solution

en: Data set on bromide oxidation by ozone in snow during metamorphism from laboratory study


maintainer: Marco Moretti

en: urbanization wild-bees honeybees zurich bettergardens beekeeping pollinator-experiment resource-availability urban-biodiversity gardens city4bees urban-ecosystems

en: Metadata for the predictors and response dataset

en: Metadata for the predictor and response dataset

en: Response and predictors used for the SEMs

en:
Mean, standard deviation , minimum and maximum values of six functional traits. The data has the following structure:
- Bee_species: bee scientific name
- mean_trait: mean value
- sd_trait: standard deviation
_min_trait: minimum value
- max_trait: maximum value
- trait: the six traits used. ITD = Intertegular distance (mm); Pheno.duration = duration of the phenology; Pheno.start = start of the phenology; Pheno.end = end of the phenology; Rel.tongue.length= relative tongue length (mm); Daytime.activity = daily activity (h)
DistNichePart: functional dissimilarity in a multidimentional trait space defined by five traits between wild bee individuals of a wild bee species and the honeybee individuals. Values are at the species level.

en:
The data contains the following columns:

  SiteID: numeric, the identity of the garden where the bee was sampled.
  Taxonname: character, species name
  Lecty: numeric discreet, a mesure of specialisation degree. Levels: 1 = monolectic; 2=oliolectic; 3=polilectic
  Pheno.start= numeric discreet, the starting week of the flying activity.
  Pheno.end= numeric discreet, the ending week of the flying activity.
  Pheno.duration=numeric discreet, Pheno.end - Pheno.start
  Daytime.activity = numeric continuous, individual-level measure of the day hour when the bee was sampled
  ITD = numeric continuous, individual-level measure of the InterTegular Distance, a metric of mobility in mm.
  Rel.tongue.length =numeric continuous, individual-level measure of the relative tongue length (Tongue length / ITD), a metric of specilaisation degree, in mm.

en:
Cities are socio-ecological systems that filter and select species, thus establishing unique species assemblages and biotic interactions. Urban ecosystems can host richer wild bee communities than highly intensified agricultural areas, specifically in resource-rich urban green spaces such as allotment and family gardens. At the same time, urban beekeeping has boomed in many European cities, raising concerns that the fast addition of a large number of managed bees could deplete the existing floral resources, triggering competition between wild bees and honeybees. The data has been used to investigated the interplay between resource availability and the number of honeybees at local and landscape scales and how this relationship influences wild bee diversity.

This dataset contains the raw and processed data supporting the findings from the paper: "Low resource availability drives feeding niche partitioning between wild bees and honeybees in a European city".

The data contains:
1. Bee trait measurements at the species and individual-level of five functional traits.
2. The values of the feeding niche partitioning (functional dissimilarity to honeybees)
3. The predictors of resource availability and beekeeping intensity at local and landscape scales used in the modelling of the paper for the 23 experimental sites.

beekeeping

display_name: beekeeping

bettergardens

display_name: bettergardens

gardens

display_name: gardens

honeybees

display_name: honeybees

pollinator-experiment

display_name: pollinator-experiment

resource-availability

display_name: resource-availability

urban-biodiversity

display_name: urban-biodiversity

urban-ecosystems

display_name: urban-ecosystems

urbanization

display_name: urbanization

wild-bees

display_name: wild-bees

zurich

display_name: zurich

educ

en: Education, culture and sport

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: Niche partitioning between wild bees and honeybees


maintainer: Frank Techel

en: avalanche-statistics avalanche-fatalities long-term-trends avalanche-accident-statistics avalanche-accident

en: Data avalanche fatalities 1969/1970 - 2014/2015

media_type: text/plain

en: The data comprises avalanche fatality data from the Alpine regions in the seven countries Austria, France, Germany, Italy, Liechtenstein, Slovenia, and Switzerland.

format: TXT

mimetype: text/plain

en:
During the last 45 years, about 100 people lost their lives in avalanches in the European Alps each year. Avalanche fatalities in settlements and on transportation corridors have considerably decreased since the 1970s. In contrast, the number of avalanche fatalities during recreational activities away from avalanche-secured terrain doubled between the 1960s and 1980s and has remained relatively stable since, despite a continuing strong increase in winter backcountry recreational activities.

Data complementing Figure 2 in:
_"Avalanche fatalities in the European Alps: long-term trends and statistics"_, by Techel, F., Jarry, F., Kronthaler, G., Mitterer, S., Nairz, P., Pavšek, M., Valt, M., and Darms, G.

Data description:
please refer to section 2 (Data and Methods) in the mentioned publication

avalanche-accident

display_name: avalanche-accident

avalanche-accident-statistics

display_name: avalanche-accident-statistics

avalanche-fatalities

display_name: avalanche-fatalities

avalanche-statistics

display_name: avalanche-statistics

long-term-trends

display_name: long-term-trends

en: Avalanche fatalities in the European Alps (1969/1970 - 2014/2015)


maintainer: Mathias Bavay

en: air-pressure wind-speed treeline stillberg air-temperature wind-direction dew-point-temperature shortwave-radiation precipitation-sum snow-depth soil-temperature relative-humidity

en: The Davos Stillberg meteorological data (Stillberg_meteo_v1.1.zip) consists of timeseries from five meteorological stations operated at the Stillberg treeline ecological research site. Metadata are provided in the folder 'Stillberg_meteo_metadata' and the data in the folder 'Stillberg_meteo_data' with subfolders for the three quality levels 'raw_standardized', ‘raw_edited’, and ‘filtered’.

en: The Davos Stillberg meteorological data (Stillberg_meteo_v2.0.zip) consists of timeseries from five meteorological stations operated at the Stillberg treeline ecological resea...

en: The Davos Stillberg meteorological data (Stillberg_meteo_v1.zip) consists of timeseries from five meteorological stations operated at the Stillberg treeline ecological research site. Metadata are provided in the folder 'Stillberg_meteo_metadata' and the data in the folder 'Stillberg_meteo_data' with subfolders for the three quality levels 'raw_standardized', ‘raw_edited’, and ‘filtered’.

en:
 Background information
The Stillberg ecological treeline research site is located in the transition zone between the relatively humid climate of the Northern Alps and the continental climate of the Central Alps. In 1975, 92,000 seedlings of the high-elevation conifer species *Larix decidua* Mill. (European larch), *Pinus cembra* L. (Cembran pine), and *Pinus mugo* ssp. *uncinata* (DC.) Domin (mountain pine) were systematically planted across an area of 5 hectares along an elevation gradient of about 150 metres, with the aim to develop ecologically, technically, and economically sustainable afforestation techniques at the treeline to reduce the risk of snow avalanches. In the course of time, additional research aspects gained importance, such as the ecology of the treeline ecotone under global change. Alongside the ecological long-term monitoring of the afforestation, several meteorological stations have recorded local meteorological conditions at the Stillberg research site. Here, we provide the Davos Stillberg meteorological timeseries of five stations from 1975 (01-01-1975), the year of the afforestation establishment, until the end of the year 2022 (31-12-2022).

 Station description
The five meteorological stations were all installed at the same location (46°46′25.015″N 9°52′01.792″E) at 2090 m a.s.l., in the lower part of the afforestation area. In general, the five stations were operated sequentially (Stillberg_meteo_metadata_stations_v1.csv). However, there are some overlapping time periods when more than one station was operated in parallel. The stations have recorded environmental parameters, such as air and soil temperature, dew point temperature, air pressure, relative humidity, wind direction and velocity, radiation, precipitation, and snow depth (Stillberg_meteo_metadata_parameters_v1.csv). The meteorological measurements were recorded hourly from 1975 until 1996 and have been recorded in 10-minute intervals since 1997.

 Data description
We processed the Davos Stillberg meteorological timeseries with the MeteoIO meteorological data pre-processing library (Bavay & Egger, 2014). Data files are provided for each station and quality level separately and named according to the station (see ‘Stillberg_meteo_metadata_stations_v1.csv’). From the raw data in their original formats, we generated three data quality levels: raw standardized (folder ‘raw_standardized’), edited (folder ‘raw_edited’) and filtered (folder ‘filtered’). The processing level is indicated in the headers of the data files. The whole processing protocol is described in a set of human-readable configuration files that are used by MeteoIO to generate the required data quality levels. This improves long-term reproducibility (Bavay et al., 2022), as the data could be regenerated in the future, even using a completely different software, to account for additional data points or to introduce new data corrections. The first quality level (raw standardized) is generated by parsing the original data files and interpreting them in order to convert all data points to a common format and meteorological parameter naming scheme, while excluding unreadable or duplicated data lines. The generated data files are derivatives of CSV files, with a standardised header that contains the metadata that are necessary to interpret and use the data (use metadata) and to populate a data index (search metadata). The latter is a textual implementation of the Attribute Convention for Data Discovery (ACDD) metadata standard (Attribute Convention for Data Discovery 1-3, 2022). The second quality level (edited) builds on the raw data by performing low-level data editing, such as removing some data periods that are known to be unusable (often based on maintenance records or anecdotal evidence) or applying undocumented calibration factors (for example, when there seems to be an obvious offset on a measured parameter for a period between two documented maintenance operations). The third quality level is generated by applying statistical filters on the data (per station and per meteorological parameter) to exclude presumably wrong values. We did not perform gap filling, as no single strategy could be relied upon that would work best for all possible data usage scenarios.

air-pressure

display_name: air-pressure

air-temperature

display_name: air-temperature

dew-point-temperature

display_name: dew-point-temperature

precipitation-sum

display_name: precipitation-sum

relative-humidity

display_name: relative-humidity

shortwave-radiation

display_name: shortwave-radiation

snow-depth

display_name: snow-depth

soil-temperature

display_name: soil-temperature

stillberg

display_name: stillberg

treeline

display_name: treeline

wind-direction

display_name: wind-direction

wind-speed

display_name: wind-speed

en: Long-term meteorological station Stillberg, Davos, Switzerland at 2090 m a.s.l.


en: natural-climate-variability snow-water-equivalent rain-on-snow cryosphere uncertainty climate-change-impact

en: climate realizations

en: The dataset contains for three variables (snow water equivalent, surface water input and liquid precipitation) 50 realizations of current and future climate periods for two time horizons (mid end end of century), two emission senarions (RCP 4.5 and 8.5) and 10 climate model chains (all EUR11 chains within CH2018). To quantify natural climate variability for projections of snow conditions and resulting rain-on-snow (ROS) flood events, a weather generator was applied to simulate inherently consistent climate variables for multiple realizations of current and future climates at 100 m spatial and hourly temporal resolution over a 12 x 12 km high-altitude study area in the Swiss Alps. The output of the weather generator was used as input for subsequent simulations with an energy balance snow model. The data was extracted in 2021 from original model output.

climate-change-impact

display_name: climate-change-impact

cryosphere

display_name: cryosphere

natural-climate-variability

display_name: natural-climate-variability

rain-on-snow

display_name: rain-on-snow

snow-water-equivalent

display_name: snow-water-equivalent

uncertainty

display_name: uncertainty

en: Multiple realizations of daily snow water equivalent, surface water input and liquid precipitation projections for mid- and late-century


en: specific-surface-area microct intercomparison icecube snow

en: The file contains all data for the plot in figure 1 of the co-located density and SSA measurement at Summit, Greenland in May 2015, conducted by Schneebeli, Proksch, Schmid.

en: Influence of black carbon on the measurement of the specific surface of snow with the near infrared method - Philipp Baumann

en: These experiments investigate different possible errors by measuring the specific surface area (SSA) of snow with the near infrared method. The influence of black carbon on this measurement was disproved, the theory was therefore validated. It shows that sunlight has no influence on the measurement with the IceCube. Wind drift snow as thin layer on the sample can increase the measured value of the specific surface by a few percent. This could be part of the error, observed in former measurements at SLF. As reference for the measured SSA values with the near infrared method, computer tomography was used to verify the SSA values.

en: Grain survey of snow particles

en: The file contains all measurements of the artificial particles for snow types A, B, C, D and E, including the artificial particle mass in the Petri dish.

en: SMP measurement files for snow types A, B, C and D.

en: artificial_snow_data set

en: The file contains all micro-CT and IC measurements of the artificial snow data set. The data set was originally produced to investigate the influence of black carbon on SSA measurements.

en: TARTES simulations

en: Jupyter notebook with TARTES simulations.

en: IceCube and MicroCT data

en: The file contains all micro-CT and IC measurements for snow types A, B, C, D and E used in our study. All samples were processed (and stored) in a cold laboratory at -15°C.

en: Statistics of samples

en: The specific surface area (SSA) of different snow types were measured with the IceCube instrument and the Scanco Medical microCT 40. In addition, the snow particles created during the preparation of IceCube samples were counted. The difference in SSA between these instruments is explained by the formation of the surface particles. A numerical simulation using TARTES simulation support the observations.

icecube

display_name: icecube

intercomparison

display_name: intercomparison

microct

display_name: microct

snow

display_name: snow

specific-surface-area

display_name: specific-surface-area

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

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


maintainer: Matthias Haeni

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

en: Site aerial photo (Beatenberg)

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

mimetype: image/jpeg

en: Site photo of the field station (BAF)

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

en: Data plot of the forest station (BAB)

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

en: Forest station data (BAB)

en: Data file.

en: Data plot of the field station (BAF)

en: Field station data (BAF)

en: Site photo of the forest station (BAB)

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

en:
High quality meteorological data are needed for long-term forest ecosystem research, particularly in the light of global change. The long-term data series published here comprises almost 20 years of measurements for two meteorological stations in Beatenberg in Switzerland where one station is located within a natural coniferous forest (BAB) with Norway spruce (_Picea abies_; 190-210 yrs) as dominant tree species. A second station is situated in the very vicinity outside of the forest (field station, BAF). The meteorological time series are presented in hourly time resolution of air temperature, relative humidity, precipitation, photosynthetically active radiation (PAR) and wind speed. Beatenberg is part of the Long-term Forest Ecosystem Research Programme (LWF) established and maintained by the
Swiss Federal Research Institute WSL.

climate-change

display_name: climate-change

long-term

display_name: long-term

meteo

display_name: meteo

meteorology

display_name: meteorology

par

display_name: par

precipitation

display_name: precipitation

relative-humidity

display_name: relative-humidity

swiss-forest-lab

display_name: swiss-forest-lab

temperature

display_name: temperature

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


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

en: Details about the different folders

en: Stable water isotopes in vapor

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

en: Weather

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

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

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

en: Stable water isotopes in snow

en: Stable water isotopes in snow measurements

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

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

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

latent-heat

display_name: latent-heat

natural-tracer

display_name: natural-tracer

paleoclimate

display_name: paleoclimate

post-depositional-processes

display_name: post-depositional-processes

stable-water-isotopes

display_name: stable-water-isotopes

sublimation

display_name: sublimation

weissfluhjoch

display_name: weissfluhjoch

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Stable Water Isotopes in snow and vapor on the Weissfluhjoch


causal-effect-of-lup

title_for_slug: causal-effect-of-lup

en: annual-effect plan-evaluation land-use-planning land-use-change psm-did

en: causal-effect-of-lup

en:
Data.dbf contains the varibles of 1662 villages in Zhangzhou Cities
PSM_DID.R can re-produce the results

en:
Title: Does zoning contain built-up land expansion? Causal evidence from Zhangzhou City, China.

Research objective: Built-up land zoning is an imporatant policy measure of land use planning (LUP) to contain built-up land expansion in China. We used a difference-indifference model with propensity score matching to estimate the average and annual effect of built-up land zoning on built-up land expansion in Zhangzhou City, China between 2010 and 2020.

Data: Data.dbf contains the varibles of 1662 villages in Zhangzhou Cities in 1995, 2000, 2005, 2010, 2013, 2015, 2018, and 2020.
XZQDM2 is villages' unique administrative ID;
Area is the land area of village i;
Dis2water is the Euclidean distance from village i to the nearest waterbody;
Dis2coastl is the Euclidean distance from village i to the nearest coastline;
Dis2city is the the Euclidean distance from village i to the city center;
Dis2county is the the Euclidean distance from village i to the nearest county center;
Elevation is the the average elevation within village i;
Dis2road is the the Euclidean distance from village i to the nearest road;
Nei_Built_ is the the area of built-up land (Nei Built.upit) in the neighboring villages of village i in year t;
Treated is a binary variable, Treated = 1 to the villages that were partially or entirely located inside the development-permitted zones, and Treated = 0 to the villages that were entirely located outside the development-permitted zones;
Intensity is the percentage of land that was assigned to the development-permitted zones in village i;
Year represent the year in 1995, 2000, 2005, 2010, 2013, 2015, 2018, and 2020;
BuLE is the dependent variable, representing built-up land expansion in village i in year t;
Town is town' unique administrative ID.

Method: First, we employed propensity score matching to overcome the selection bias and satisfy the parallel trend
assumption. Second, we built four Difference-in-Difference models to estimate the average and annual effect.

annual-effect

display_name: annual-effect

land-use-change

display_name: land-use-change

land-use-planning

display_name: land-use-planning

plan-evaluation

display_name: plan-evaluation

psm-did

display_name: psm-did

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: Causal effect of LUP


en: snow-depth snow-water-equivalent model-validation-data model-calibration-data snow-density

en: swe2hs calibration and validation data

en:
 Contents of this archive

 README.md

 A readme file containing a description of the other files.

 automatic_stations_validation_data.csv

 Daily data from the automatic weather stations used for validation.

 Data sources and providers are listed in the SWE2HS paper, Table 1.

 Columns:

  * date: timestamp of the format YYYY-MM-DD
  * HS_[m]: snow depth in [m]
  * SWE_[m]: snow water equivalent in [m]
  * site_id: code for the measurement site
  - HS_interpolated: boolean flag for interpolated values
  - SWE_interpolated: : boolean flag for interpolated values

 coordinates_swe2hs_paper_stations.csv

 Coordinates of all measurement sites used in the SWE2HS paper.

 Columns:

  * site_id: code for the measurement site
  * location_name: location name of the measurement site
  * lon_[wgs84]: longitude in WGS84
  * lat_[wgs84]: latitude in WGS84
  * elevation_[m]: atitude of the measurement site in [m]
  * dataset: dataset in which the station is included

 manual_station_calibration_data.csv

 Calibration data from the Swiss manual observer stations. Long station records
 trmmed to 15 random water-years.

 Columns:

  * date: timestamp of the format YYYY-MM-DD
  * HS_[m]: snow depth in [m]
  * SWE_[m]: snow water equivalent in [m]
  * SWE_from_profile_[m]: biweekly measurement from the SWE in the measurement
   pit. Corrected with the bulk density in order to match the snow depth from
   the measurement stake. Please refer to the paper for further details.
  * site_id: code for the measurement site

 manual_stations_validation_data.csv

 Validation data from the Swiss manual observer stations. Remaining water-years
 from the long sttaion record not used for the calibration data.

 Columns:

  * date: timestamp of the format YYYY-MM-DD
  * HS_[m]: snow depth in [m]
  * SWE_[m]: snow water equivalent in [m]
  * SWE_from_profile_[m]: biweekly measurement from the SWE in the measurement
   pit. Corrected with the bulk density in order to match the snow depth from
   the measurement stake. Please refer to the paper for further details.
  * site_id: code for the measurement site

 manual_stations_used_deltasnow_parameters.csv

 Site specific DeltaSNOW model parameters used for the calculation of the daily
 SWE series in 'manual_stations_calibration_data.csv' and
 'manual_stations_validation_data.csv'.

 Parameters were optimized by reducing RMSE to the biweekly snow pit SWE measurements.
 Please refer to the paper for further details.

 A python implementation of the DeltaSNOW model was used for the calculation:
 pydeltasnow (https://pypi.org/project/pydeltasnow/)

 For a detailed description of the DeltaSNOW model parameters please see the
 respective publication:

 Winkler, M., Schellander, H., & Gruber, S. (2021). Snow water equivalents
 exclusively from snow depths and their temporal changes: the Δ snow model.
 Hydrology and Earth System Sciences, 25(3), 1165-1187.

 Columns:

  * site_id: code for the measurement site
  * rho_max: Maximum density of an individual snow layer produced by the DeltaSNOW model in [kg/m3]
  * rho_null: Fresh snow density for a newly created layer [kg/m3
  * c_ov: Overburden factor due to fresh snow [-]
  * k_ov: Defines the impact of the individual layer density on the compaction due to overburden [-]
  * k: Exponent of the exponential-law compaction [m3/kg]
  * tau: Uncertainty bound [m]
  * eta_null: Effective compactive viscosity of snow for “zero-density” [Pa s]

en:
The data in this repository was used for the calibration and validation of the
SWE2HS model in the following publication:
Aschauer, J., Michel, A., Jonas, T., & Marty, C. (2023). An empirical model to calculate snow depth from daily snow water equivalent: SWE2HS 1.0. Geoscientific Model Development Discussions, 1-19. https://doi.org/10.5194/gmd-2022-258
Contains daily snow water equivalent and snow depth timeseries from stations in the European Alps.

model-calibration-data

display_name: model-calibration-data

model-validation-data

display_name: model-validation-data

snow-density

display_name: snow-density

snow-depth

display_name: snow-depth

snow-water-equivalent

display_name: snow-water-equivalent

en: SWE2HS model calibration and validation data


swiss-fluxnet-site-lageren

title_for_slug: swiss-fluxnet-site-lageren

en: forest flux swiss-fluxnet eddy-covariance meteorological-measurements greenhouse-gas swiss-forest-lab

en: Fluxnet2015 dataset

en: The flux and meteo data of the Swiss FluxNet site Lägeren are also part of the Fluxnet2015 dataset.

en: Site description

en: Data access

en: The flux and meteo data of the Swiss FluxNet site Lägeren are openly available from the European Fluxes Database Cluster.

en:
The Swiss FluxNet Site Lägeren is a managed mixed deciduous mountain forest located on the steep Lägeren mountain (NW of Zurich, Swiss Plateau). The forest is highly diverse, dominated by beech, but also including ash, maple, spruce and fir trees. Eddy covariance flux measurements were started in April 2004. The site was part of the international CarboEurope IP network and the National Air Pollution Monitoring Network (NABEL). In addition to Swiss FluxNet, the site is part of the Long-term Forest Ecosystem Research (LWF) of WSL and the biological drought and growth indicator network (TreeNet) of WSL.

Measurements

- Ecosystem flux measurements of CO2, H2O vapour are performed with the eddy-covariance method. This method is based on measurements of trace gas mixing ratios, using infrared gas analyzers (for CO2, H2O vapor), combined with wind speed and wind direction measurements, using 3D sonic anemometers. To resolve the short-term turbulent fluctuations in the atmosphere, very fast measurements are needed: we measure at 10-20 Hz, i.e., 10-20 times per second. To assess the energy budget of each ecosystem, also radiation sensors and soil climate profiles are installed at the site.
- Sub-canopy eddy fluxes (CO2, H2O), soil respiration campaigns
- Continuous CO2 profile measurements.
- Auxiliary micrometeorology and soil climate measurements.

Data availability

All data are available from the European Fluxes Database Cluster, but are also part of Fluxnet2015 dataset.

Data policy

ICOS data license:
[https://www.icos-cp.eu/data-services/about-data-portal/data-license](https://www.icos-cp.eu/data-services/about-data-portal/data-license)

Detailed site info:
[https://www.swissfluxnet.ethz.ch/index.php/sites/ch-lae-laegeren/site-info-ch-lae//](https://www.swissfluxnet.ethz.ch/index.php/sites/ch-lae-laegeren/site-info-ch-lae/)

eddy-covariance

display_name: eddy-covariance

flux

display_name: flux

forest

display_name: forest

greenhouse-gas

display_name: greenhouse-gas

meteorological-measurements

display_name: meteorological-measurements

swiss-fluxnet

display_name: swiss-fluxnet

swiss-forest-lab

display_name: swiss-forest-lab

en: Swiss FluxNet Site Lägeren


maintainer: Giulia Mazzotti

en: snow-energy-balance forest-snow snow-modelling snow-distribution forest-structure

en: Canopy data

en: Includes all canopy structure metrics computed at all surveyed locations

en: Snow depth data forest sites

en: Contains all snow depth data acquired at all forest sites during winters 2012/13, 2013714 and 2014/15.

en: Coordinates data

en: Includes coordinate information of each surveyed point (7 x 276 = 1932 forest points, 2 open areas)

en: Grid point nomenclature

en: The grid structure of each forest site and the respective nomenclature applied to label individual points within each grid are outlined here.  Each point marks the location of snow depth measurements and co-located canopy structure parameters.

en: Snow depth data open sites

en: Contains all snow detph data accquired at open sites nearby the forest stands during the same period.

en: Meteo data to drive FSM2 model

en: Includes meteorological driving data used to run the forest snow model FSM2, processed and parsed in text files that match the format required by FSM2. Separate files are available for each field area (Laret, Drusatscha and Ischlag) and season (2012/13, 2013/14, 2014/15). Each file includes data for the period from October 1st to May 31st at hourly resolution.

en:
This dataset contains all snow, canopy and meteorological data presented and used in the publication: Mazzotti, G., Essery, R., Moeser, D. & Jonas T. (2020) 'Resolving spatial variability of forest snow using an energy-balance model with a 1-layer canopy'. Water Resources Research, https://doi.org/10.1029/2019WR026129.
This publication must be cited when using this dataset.

forest-snow

display_name: forest-snow

forest-structure

display_name: forest-structure

snow-distribution

display_name: snow-distribution

snow-energy-balance

display_name: snow-energy-balance

snow-modelling

display_name: snow-modelling

educ

en: Education, culture and sport

envidat

en: EnviDat: the Environmental Data Portal

en: Snow depth, canopy structure and meterorological datasets from the Davos area, Switzerland, Winters 2012/13-2014/15, used for high-resolution forest snow modelling


maintainer: Jana von Freyberg

en: catchment hydrology precipitation stable-water-isotope stream-water alptal

en: Single comma-separated txt-file with daily data of stable water isotopes measured in streamwater and precipitation in the Swiss Alptal catchment and two of its tributaries (Erlenbach and Vogelbach), as well as daily data of key hydrometeorological variables from the Alptal valley

en: Data Set Documentation: alptal_Isotopes_daily_2015-2019_v2

en: Documentations contains meta data of the data columns and information about the study catchments.

en:
Version 2 of the dataset contains additional information on:
- line-conditioned (lc) excess
- precipitation sample volumes
- isotope data quality

en: This dataset contains four years of daily measurements of the natural isotopic composition (2H, 18O) of precipitation and stream water at the Alp catchment (area 47 km2) in Central Switzerland and two of its tributaries (0.73 km2 and 1.55 km2). In addition, the dataset contains daily measurements of key hydrometeorological variables.

alptal

display_name: alptal

catchment

display_name: catchment

hydrology

display_name: hydrology

precipitation

display_name: precipitation

stable-water-isotope

display_name: stable-water-isotope

stream-water

display_name: stream-water

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: Four years of daily stable water isotope data in stream water and precipitation from three Swiss catchments


maintainer: Anita C Risch

en: soil-properties microbial-respiration microbial-biomass-carbon herbivore-exclusion nutrient-addition anthropogenic-management 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

climate

display_name: climate

herbivore-exclusion

display_name: herbivore-exclusion

microbial-biomass-carbon

display_name: microbial-biomass-carbon

microbial-respiration

display_name: microbial-respiration

nutrient-addition

display_name: nutrient-addition

soil-properties

display_name: soil-properties