en: snow-hydrology avalanches snow-depth-distribution davos dischma photogrammetry

en: Snowdepth maps

format: GeoTIFF

en:
The available datasets are snow depth maps with a spatial resolution of 0.5 m derived from images of the survey camera Vexcel Ultracam mounted on a piloted airplane. Image acquisition was carried out during the approximately peak of winter (time when the thickest snowpack is expected) in spring. The snow depth maps are calculated by the subtraction of a summer-DTM from the processed winter- DSM of the corresponding date. The summer-DTM used was derived from a point cloud of an airborne laser scanner from 2020.

Due to the occurrence of inaccuracies of the calculated snow depth values caused by the photogrammetric method, we applied different masks to significantly increase the reliability of the snow depth maps. We masked out settled areas, high-frequented streets and technical constructions, pixels with high vegetation (height > 0.5 m) , outliers and unrealistic snow depth values. In addition, we modified the snow depth values of snow-free pixels to 0.
The information on buildings and infrastructure comes from the exactly classified ALS point cloud and the TLM dataset from Swisstopo (https://www.swisstopo.admin.ch/de/geodata/landscape/tlm3d.html#links). High vegetation is also derived from the classification and the calculated object height from the point cloud.
Outliers and unrealistic snow depth values are defined as negative snow depth values and snow depths exceeding 10 m. The classification of each pixel of the corresponding orthophoto into snow-covered or snow-free is based on the application of a threshold of the NDSI or manually determined ratios of the RGB values.

An extensive accuracy assessment proves the high accuracy of the snow depth maps with a root mean square error of 0.25 m for the year 2017 and 0.15 m for the following snow depth maps.

The work is published in:

avalanches

display_name: avalanches

davos

display_name: davos

dischma

display_name: dischma

photogrammetry

display_name: photogrammetry

snow-depth-distribution

display_name: snow-depth-distribution

snow-hydrology

display_name: snow-hydrology

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: Snow depth mapping by airplane photogrammetry (2017 -  ongoing)


maintainer: James W. Kirchner

en: catchment chloride stable-water-isotopes tracer streamflow deuterium precipitation plynlimon

state: active

type: dataset

en: 7-hourly stable water istopes

media_type: text/plain

en: Data set of 7-hourly stable water isotope data measured at Carreg Wen (precipitation) and Upper Hafren (streamflow).

format: TXT

mimetype: text/plain

en: Data Set Documentation

media_type: application/pdf

en: This file contains a documentation of the data sets, including references to other helpful documentations.

format: PDF

mimetype: application/pdf

en: Weekly stable water isotopes

en: Data set of weekly stable water isotope data collected at Carreg Wen (precipitation) and Lower Hafren and Tanllwyth (streamflow).

en: The data base contains timeseries of stable water isotopes in precipitation and streamflow at Plynlimon, Wales, UK. One data set contains weekly stable water isotope data from the Lower Hafren and Tanllwyth catchments, and the other data set contains 7-hourly stable water isotope data from Upper Hafren. Both data sets also include chloride concentrations in precipitation and streamflow.

catchment

display_name: catchment

chloride

display_name: chloride

deuterium

display_name: deuterium

plynlimon

display_name: plynlimon

precipitation

display_name: precipitation

stable-water-isotopes

display_name: stable-water-isotopes

streamflow

display_name: streamflow

tracer

display_name: tracer

educ

fr: Éducation, culture et sport

de: Bildung, Kultur und Sport

it: Istruzione, cultura e sport

fr: EnviDat: le portail de données environnementales

de: EnviDat: das Umwelt-Datenportal

it: EnviDat: il portale dei dati ambientali

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

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

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

political_level: other

wsl

capacity: public

type: organization

approval_status: approved

key: guid

key: uri

value: Envidat Harvester

en: Stable water isotopes in precipitation and streamflow at Plynlimon, Wales, UK


maintainer: EnviDat Support

en: temperature growing-degree-days precipitation soil-water-balance growing-season-length relative-humidity solar-radiation

en: CHELSA-BIOCLIM+

en: Global extented set of bioclimatic variables at different temporal (e.g. monthly, annual, long term normals) resolution in Geotiff format.

en: CHELSA-BIOCLIM+ Technical Documentation

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

en: A multitude of physical and biological processes on which ecosystems and human societies depend are governed by climatic conditions. Understanding how these processes are altered by climate change is central to mitigation efforts. Based on mechanistically downscaled climate data, we developed a set of climate-related variables at yet unprecedented spatiotemporal detail as a basis for environmental and ecological analyses. We created gridded data for near-surface relative humidity (hurs), cloud area fraction (clt), near-surface wind speed (sfcWind), vapour pressure deficit (vpd), surface downwelling shortwave radiation (rsds), potential evapotranspiration (pet), climate moisture index (cmi), and site water balance (swb), at a monthly temporal and 30 arcsec spatial resolution globally starting 1980 until 2018. At the same spatial resolution, we further estimated climatological normals of frost change frequency (fcf), snow cover days (scd), potential net primary productivity (npp), growing degree days (gdd), and growing season characteristics for the periods 1981-2010, 2011-2040, 2041-2070, and 2071-2100, considering three shared socioeconomic pathways (SSP126, SSP370, SSP585) and five Earth system models. Time-series variables showed high accuracy when validated against observations from meteorological stations. Climatological normals were also highly correlated to observations although some variables showed notable biases, e.g., snow cover days (scd). Together, the data sets presented here allow improving our understanding of patterns and processes that are governed by climate, including the impact of recent and future climate changes on the world’s ecosystems and associated services to societies.

growing-degree-days

display_name: growing-degree-days

growing-season-length

display_name: growing-season-length

relative-humidity

display_name: relative-humidity

soil-water-balance

display_name: soil-water-balance

solar-radiation

display_name: solar-radiation

temperature

display_name: temperature

en: CHELSA-BIOCLIM+ A novel set of global climate-related predictors at kilometre-resolution


maintainer: Felix Kienast

en: photovoltaic energy landscape marxan optimization social-acceptance wind-power ecosystem-services spatial-planning

media_type: application/zip

en:
This is the compressed Shapefile (ESRI) from the MARXAN optimization with 2216 rows (objects) and 18 columns (variables) with sf (simple feature) and data frame classes.
Column description:
X01 = Planning_units (PUs)
X02 = incl.OS. number of times pu was selected in MARXAN when optimized for ecological costs
X03 = incl.OS. number of times pu was selected in MARXAN when optimized for social costs
X04 = incl.OS. ecological costs (ECUess) of pu when realizing the total energy (normalized)
X05 = incl.OS. social costs (ECUsoc) of pu when realizing the total energy (normalized)
X06 = incl.OS. total energy (MWh/a)
X07 = incl.OS. energy from roof-mounted PV infrastructure (MWh/a)
X08 = incl.OS. energy from wind turbines (MWh/a)
X09 = incl.OS. energy from ground-mounted PV infrastructure (MWh/a)

X10 = excl.OS. number of times pu was selected in MARXAN when optimized for ecological costs
X11 = excl.OS. number of times pu was selected in MARXAN when optimized for social costs
X12 = excl.OS. ecological costs (ECUess) of pu when realizing the total energy (normalized)
X13 = excl.OS. social costs (ECUsoc) of pu when realizing the total energy (normalized)
X14 = excl.OS. total energy (MWh/a)
X15 = excl.OS. energy from roof-mounted PV infrastructure (MWh/a)
X16 = excl.OS. energy from wind turbines (MWh/a)
X17 = excl.OS. energy from ground-mounted PV infrastructure (MWh/a)
geometry = Simple feature XY geometry (SFC_POINT) each representing the center of a 4x4km planning unit (PU) in EPSG 21781 (CRS CH1903 / LV03).

format: ZIP

en:
This study examines the siting scenarios for renewable energy infrastructure (REI) in Switzerland, incorporating the external costs of ecosystem services and, innovatively, social preferences. This approach challenges the prevalent techno-economic siting paradigm, which often overlooks these externalities. To minimize the external costs of the scenarios while maximizing energy yield, Marxan, an optimization software, was employed.
MARXAN was run for 2 versions: a) without ground-mounted open space PV infrastructure (excl. OS) and b) with ground-mounted open space PV infrastructure (incl. OS). In each version optimization was done using ecological costs (ECUess) or social costs (ECUsoc) in a regular grid of 4x4km (planning unit).

File:
PU_data.shp: compressed Shapefile (ESRI) from the MARXAN optimization with 2216 rows (objects) and 18 columns (variables) with sf (simple feature) and data frame classes.

Headers are described below:
X01 = Planning_units (PUs)

X02 = incl.OS. number of times pu was selected in MARXAN when optimized for ecological costs
X03 = incl.OS. number of times pu was selected in MARXAN when optimized for social costs
X04 = incl.OS. ecological costs (ECUess) of pu when realizing the total energy (normalized)
X05 = incl.OS. social costs (ECUsoc) of pu when realizing the total energy (normalized)
X06 = incl.OS. total energy (MWh/a)
X07 = incl.OS. energy from roof-mounted PV infrastructure (MWh/a)
X08 = incl.OS. energy from wind turbines (MWh/a)
X09 = incl.OS. energy from ground-mounted PV infrastructure (MWh/a)

X10 = excl.OS. number of times pu was selected in MARXAN when optimized for ecological costs
X11 = excl.OS. number of times pu was selected in MARXAN when optimized for social costs
X12 = excl.OS. ecological costs (ECUess) of pu when realizing the total energy (normalized)
X13 = excl.OS. social costs (ECUsoc) of pu when realizing the total energy (normalized)
X14 = excl.OS. total energy (MWh/a)
X15 = excl.OS. energy from roof-mounted PV infrastructure (MWh/a)
X16 = excl.OS. energy from wind turbines (MWh/a)
X17 = excl.OS. energy from ground-mounted PV infrastructure (MWh/a)

geometry = Simple feature XY geometry (SFC_POINT) each representing the center of a 4x4km planning unit (PU) in EPSG 21781 (CRS CH1903 / LV03).

For MARXAN optimizations (for meaning of standard files in MARXAN see https://scholar.google.ch/scholar_url?url=https://courses.washington.edu/cfr590/software/Marxan1810/marxan_manual_1_8_2.pdf&hl=de&sa=X&ei=NAWHZI6JJ8PFmAGMo47oDQ&scisig=AGlGAw-DBrsV4kUzkR6GkN1jtG66&oi=scholarr), pu_data.shp was used to generate MARXAN files input.dat, pu.dat, puvfeat.dat, spec.dat.
Coordinate Reference System: CH1903 / LV03. EPSG: 21781.

ecosystem-services

display_name: ecosystem-services

energy

display_name: energy

landscape

display_name: landscape

marxan

display_name: marxan

optimization

display_name: optimization

photovoltaic

display_name: photovoltaic

social-acceptance

display_name: social-acceptance

spatial-planning

display_name: spatial-planning

wind-power

display_name: wind-power

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: Optimizing renewable energy siting in the Swiss landscape


maintainer: Elisabeth Hafner

en: avalanche validation snow avalanche-hazard-mapping satellite

en: Satellite Avalanche Mapping Validation data and data description

en: 14 shapefiles to be used only together with the data description also provided as a pdf.

en:
Validation points, validation area, ground truth coverage, SPOT 6 avalanche outlines, Sentinel-1 avalanche outlines, Sentinel-2 avalanche outlines, Davos avalanche mapping (DAvalMap) avalanche outlines as shapefiles and a detailed attribute description (DataDescription_EvalSatMappingMethods.pdf).

Coordinate system: CH1903+_LV95

The generation of this dataset is described in detail in:

Hafner, E. D., Techel, F., Leinss, S., and Bühler, Y.: Mapping avalanches with satellites – evaluation of performance and completeness, The Cryosphere, https://doi.org/10.5194/tc-2020-272, 2021.

avalanche

display_name: avalanche

avalanche-hazard-mapping

display_name: avalanche-hazard-mapping

satellite

display_name: satellite

snow

display_name: snow

validation

display_name: validation

en: Satellite avalanche mapping validation data


en: field-data restoration traffic-noise urban-landscape forest

en: 07_Studyprotocol_Test site pictures.zip

en: This folder includes two pictures of each test site. The locations for the two pictures on each route were chosen to be as characteristic as possible in terms of the visual and acoustical impressions of the walk.

en: This file contains the informed consent document that participants signed to participate in the study.

en: This document includes the questionnaire for Papers 2 and 3, administered at the first measurement time point (before the walk began).

en: 06_Studyprotocol_GPS tracks of the walking routes.zip

en: This folder contains the GPS tracks of the walking routes in the test sites.

en: 03a_Studyprotocol_Data on evaluation of acoustic impression of the test settings_Studyprotocol.sav

en: This data file includes information for an additional post hoc classification of the sound environment of the test sites. This evaluation was conducted in the afternoons of the same week in 2023 at two locations within each test site by a civil servant who was not familiar with the study's research questions.

format: OCTET

en: 09_P2_Fotos Noise ROS paper.zip

en: This folder included the example pictures of the conditions included in paper 2 ("Psychophysiological effects of walking in forests and urban built environments with disparate road traffic noise exposure: A randomized controlled trial“).

en: 03b_Studyprotocol_Script_Analysis stationary noise evaluation.sps.sps

en: This SPSS syntax file contains a script for analyzing data related to the stationary evaluation of the acoustic impressions at the test sites, which is numbered 03a in this repository.

en: This document includes the questionnaire for Papers 2 and 3, administered at the second measurement time point (after the walk).

en: This SPSS data file includes the data used for analysis in paper 2 (Psychophysiological effects of walking in forests and urban built environments with disparate road traffic noise exposure: A randomized controlled trial) and paper 3 ("The effects of walking in urban forests and urban built environments with disparate road traffic noise exposure on repetitive negative thinking and connectedness with the non-human world: A randomized controlled trial").

en: 01_Studyprotocol_Locations of the starting points of the walking routes in the test settings.pdf

en: This document includes GPS coordinates of the starting points of the walking routes in the test sites.

en: This document includes the questionnaire for Papers 2 and 3, administered at the third measurement time point (10 days after participation on site).

en: 11_P4_S1_Instructions for individual mindful walking.pdf

en: This document contains the booklet with instructions for mindful walking, which participants in the mindfulness condition were encouraged to take with them and use during their individual walks over the 10 days following on-site participation.

en: 12_P4_S2_Description of mindful walking intervention on site.pdf

en: This file includes the guided mindful walking instructions provided to participants in the mindfulness condition during the on-site walks.

en: This file includes the information sheet provided to participants before they consented to take part in the study.

en: 02_Studyprotocol_Data on Laeq from noise recordings of walking routes_Studyprotocol.sav

en: This data file contains information on the LAeq levels recorded along the final walking routes between September 2021 and March 2022, prior to the commencement of data collection.

en: 13_P4_Data mindful walking.sav

en: This file contains data for the manuscript "Effects of a mindful walking intervention in urban forests for healthy adults".

en:
This repository contains data related to the field experiments of the RESTORE project. This project is a collaboration between the Swiss Federal Institute for Forest, Snow and Landscape Research (WSL) and the Swiss Federal Laboratories for Materials Science & Technology (EMPA). It has received funding from the Swiss National Science Foundation. The overall objective of the RESTORE project was to assess the effects of green spaces as facilitators and noise as impediment to recover from stress in people's daily environments across Switzerland.
We conducted a randomized, controlled field study to compare the psychophysiological benefits of exposure to forests and urban built environments with different levels of road traffic noise in healthy adults during 30-minute walks in Zürich, Switzerland to explore:
1) The psychophysiological effects of walking in forests and urban built environments with disparate road traffic noise exposure
2) The effects of walking in urban forests and urban built environments with disparate road traffic noise exposure on repetitive negative thinking and connectedness with the non-human world
3) The effects of a mindful walking intervention in urban forests for healthy adults
Accordingly, this repository includes data and documentation for the 3 scientific articles:
Scientific article 1. “Psychophysiological effects of walking in forests and urban built environments with disparate road traffic noise exposure: study protocol of a randomized controlled trial", BMC Psychology, 2024. https://doi.org/10.1186/s40359-024-01720-x
Scientific article 2. "Psychophysiological effects of walking in forests and urban built environments with disparate road traffic noise exposure: A randomized controlled trial“. Journal of Environmental Psychology, 2025.
Scientific article 3. “The effects of walking in urban forests and urban built environments with disparate road traffic noise exposure on repetitive negative thinking and connectedness with the non-human world: A randomized controlled trial". To be submitted to PLOS One, 2025.
Scientific article 4: "Effects of a mindful walking intervention in urban forests for healthy adults". To be submitted to Current Psychology, 2025.

field-data

display_name: field-data

forest

display_name: forest

restoration

display_name: restoration

traffic-noise

display_name: traffic-noise

urban-landscape

display_name: urban-landscape


en: drone sea-ice antarctica uav snow irra snow-temperature snow-roughness

en:
This dataset accompanies the publication "How Flat is Flat? Investigating the spatial variability of snow surface temperature and roughness on landfast sea ice using UAVs in McMurdo Sound, Antarctica".
The dataset consists of 1) processed orthomosaic of RGB values, snow topography (digital elevation model), and surface temperature, 2) topography-dependent irradiance map, and 3) point data of measured radiation and atmospheric parameters in the McMurdo Sound, October-December 2022. The data was collected as part of the New Zealand Marsden Fund Research Grant 21-VUW-103 "Can Snow Change the Fate of Antarctic Sea Ice?". Unprocessed data associated with this field campaign are published in 10.16904/envidat.633.
For the airborne measurements, we used a DJI Matrice 30T, a multi-rotor (quadcopter) UAV with self-heating batteries, suitable for operations in polar conditions (temperature range from -20 to +50 °C). The UAV has a wide-angle camera (12 megapixels) for RGB images and a thermal camera (Uncooled Vox Microbolometer, long-wave infrared spectrum 8-14 µm, 1.31 megapixels) for TIR images with an accuracy of ±2 °C or ± 2 %.
GNSS data includes high-precision location information (latitude, longitude, ellipsoidal height, receiver clock offset, and timestamps formatted in decimal hours and day-of-year notation).
WE calculated the distribution of received solar irradiance over the sea-ice surface using a terrain-corrected radiative transfer model and measured shortwave radiation components, accounting for surface roughness, aspect, slope, and solar position during the flight, while applying terrain shading and measured diffuse radiation.
The broadband radiation station data consists of shortwave (Sw) and longwave (Lw) radiation measurements (Kipp & Zonen CMP22 pyranometers and CGR4 pyrgeometers).
Surface temperature readings from cold, sediment, and hot targets were recorded using Apogee SI-121-SS sensors.
This dataset also contains automated weather station (AWS) data (air and snow surface temperatures, relative humidity, wind direction, and wind speed).
For detailed information, please use the README file.

antarctica

display_name: antarctica

drone

display_name: drone

irra

display_name: irra

sea-ice

display_name: sea-ice

snow-roughness

display_name: snow-roughness

snow-temperature

display_name: snow-temperature

uav

display_name: uav

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: Snow On Antarctic Sea Ice - McMurdoSound 2022 - UAV Retrievals of Snow Topography and Snow Surface Temperature


en: polar sea-ice drone-data spatial-variability snow antarctica

en:
Data of snow and sea ice in the McMurdo Sound, October-December 2022. The data was collected as part of the New Zealand Marsden Fund Research Grant 21-VUW-103 "Can Snow Change the Fate of Antarctic Sea Ice?"
The dataset includes raw data of the manual snow and sea ice measurements from snow pits and ice cores (temperature, density, salinity, dO18), measurements of snow water equivalent (SWE), spatial information of snow height (MagnaProbe) and sea ice thickness (EM-31), AWS (air temperature, wind speed, wind direction, relative humidity, pressure), radiations stations (shortwave, longwave, thermal IR, spectral shortwave), differential GPS data (3 fixed stations on different sea ice thicknesses, + 1 rover station for georeferencing UAV measurements), SIMBA buoy temperature (+heated temperature) data (3 buoys during November, 1 buoy for 15 months), UAV data: RGB, thermal IR, broadband albedo, spectral albedo, Chlorophyll-a from ice cores (bottom 10 cm), NIR reflectivity data of snow at 850 nm, and 940 nm (snow surface, profile, ice surface), photographs (1. overview of field sites, 2. for Structure from Motion for surface roughness, 3. macrophotos of snow) surface impurity concentrations, microCT data of snow microstructure, Denoth probe (density) and InfraSnow (specific surface area - SSA).
See the README file in each dataset for detailed information.

drone-data

display_name: drone-data

polar

display_name: polar

spatial-variability

display_name: spatial-variability

en: Snow on Antarctic Sea Ice - McMurdo Sound 2022


maintainer: Christoph Marty

en: snow climate-change-impacts snowfall hydrology water-cycle

en: Synthesis report chapter

en: DORA Persistent URL

format: SERVICE

en:
This report was prepared as one of the synthesis report chapters of the Hydro-CH2018 project of the Federal Office for the Environment (FOEN). An important feature of snow cover is the fact that its volume and duration is subject to large year-to-year fluctuations. As frozen precipitation, snow cover is nothing other than a natural water reservoir that delays precipitation to runoff and is thus of outstanding importance for the seasonal water balance in Switzerland. Over a whole year, approximately 40% (22 km3) of the annual runoff currently comes from snow melting and only 1% from glacier melting. Typically, the snow cover in the Alpine region builds up over the autumn and winter months, reaches its maximum between February and May, depending on the altitude, and dominates the runoff processes during melting in the following spring and summer months.
Due to the great dependence on minus temperatures and precipitation, the snow cover reacts sensitively to temperatures above 0° Celsius and more or less precipitation. Due to climate change and the associated warming, the proportion of precipitation that falls as snow decreases measurably. In addition to this reduction in snowfall, the warmer temperatures also cause the snow cover to melt more quickly. The decline in snowfall has so far mainly affected lower altitudes, where winter temperatures often reach positive levels.
As climate change progresses, this trend is likely to continue and above all affect higher zones. Even at higher altitudes, the snow cover will then start later, melt away earlier and is increasingly no longer permanently present. This development will also have an effect on the water bodies. Today nival regimes, i.e. regimes shaped by snow, are shifting towards pluvial regimes, i.e. regimes dominated by rain. Overall, winter runoff increases, summer runoff decreases. By the end of the century, the proportion of runoff from snowmelt will decrease throughout Switzerland, albeit to a lesser extent than the proportion from glacier melt.

climate-change-impacts

display_name: climate-change-impacts

hydrology

display_name: hydrology

snow

display_name: snow

snowfall

display_name: snowfall

water-cycle

display_name: water-cycle

en: Snow and the water cycle in a changing climate


maintainer: Henning Löwe

en: discrete-element-modelling angle-of-repose fast-sintering granular-snow grain-shape

en: Angle of repose and 3D CT images

en:
5 Snow types were examined and glass beads as a cohesion-less reference, named as follows:

RG: rounded grains (sieved with 0.7 - 1.4 mm)

FC: faceted grains (1.4 - 2 mm)

FCS: small faceted grains (0.7 - 1.4 mm)

IB: ice beads (4 mm)

IBS: small ice beads (2 - 2.2 mm)

GB: glass beads(2 mm)

en:
Angle of repose experiments were performed with different snow types at temperatures between -2 and -40°C. They were used to examine granular snow dynamics on the grain-scale with focus on the role of grain shape and cohesion. The angle of repose was observed by sieving snow onto a round, freestanding base until a stationary heap was formed.

This dataset consists of 1) the images of the experimental heaps that were taken to determine the angle of repose, 2) one binary 3D micro computed tomography image of each snow type. The CT images were taken with the SLF micro-CT40 to characterize snow properties and grain shape.

The experiments with natural snow types (rounded and faceted grains) and spherical model snow allowed for an examination of the differences in granular properties between natural grain shapes and spherical particles in view of Discrete Element Modelling. With the chosen temperatures, the effect of sintering could be observed that increases the angle of repose with increasing temperature.

angle-of-repose

display_name: angle-of-repose

discrete-element-modelling

display_name: discrete-element-modelling

fast-sintering

display_name: fast-sintering

grain-shape

display_name: grain-shape

granular-snow

display_name: granular-snow

en: Angle of repose experiments with natural and spherical snow


maintainer: Anita Risch

en: exclusion invertebrates grassland vertebrates subalpine grazing ecosystem-multifunctionality ecosystem-coupling swiss-national-park ecosystem-functions

en: Datafile

en:
ID: Plot ID composited of Fence & Treatment number (90 plots total)

Grassland: Name of the six subalpine grasslands where the exclosure setups were erected. In large grasslands (GRIM, MING, STAB), we established four hierarchical exclosure setups (two in short- and two in tall-grass vegetation), whereas in the small grasslands (BOTS, DADA, MARG) only one setup was erected per vegetation type. BOTS = Val dal Botsch, GRIM = Grimmels, MING = Val Mingér, STAB = Alp Stabelchod, DADA = Stabelchod da daint, MARG = Margunet

Fence: Number of exclosure setup (1 - 18). Note that uneven numbers (1,3,5,7,9,11,13,15,17) represent setups established on short-grass vegetation, even numbers (2,4,6,8,10,12,14,16,18) setups established in tall-grass vegetation
Vegetation	Vegetation type, 1 = short-grass vegetation, 2 = tall-grass vegetation; (9 replicate fences on each vegetation type)
Treatment	Progressive animal exclusion, numbers indicate which animals had access to the plot: 1 = Large/medium/small mammals, and invertebrates have access ("L/M/S/I"), 2 = Medium/small mammals, and invertebrates have access (“M/S/I”), 3 = Small mammals and invertebrates have access (“S/I”), 4 = Invertebrates have access (“I”), 5 = No animals have access (“None”)

Bulkden:	Fine fraction soil bulk density (g cm-3)

Soil.moist	: Soil moisture (%)

Soil.org.C: Rhizosphere soil C content (%)

pH: Soil pH

Nema.PC1: First principal component (PC1) of principal component analysis (PCA) of the soil nematode community composition

Nema.PC2: Second principal component (PC2) of principal component analysis (PCA) of the soil nematode community composition

Micro.PC1: 	First principal component (PC1) of principal component analysis (PCA) of the soil microbial community composition

Micro.PC2: Second principal component (PC2) of principal component analysis (PCA) of the soil microbial community composition

Abv.Inv.PC1: First principal component (PC1) of principal component analysis (PCA) of the aboveground-dwelling invertebrate community composition

Abv.Inv.PC2: Second principal component (PC2) of principal component analysis (PCA) of the aboveground-dwelling invertebrate community composition

Arthr.PC1: First principal component (PC1) of principal component analysis (PCA) of the soil arthropod community composition

Arthr.PC2: Second principal component (PC2) of principal component analysis (PCA) of the soil arthropod community composition

Plant.PC1: First principal component (PC1) of principal component analysis (PCA) of the plant community composition

Plant.PC2: Second principal component (PC2) of principal component analysis (PCA) of the plant community composition

Root.biomass: Root biomass (g m-2)

Micro.biomass	: Microbial biomass carbon (mg kg-1)

Plant.richness: Plant species richness

Plant.N: Plant tissue nitrogen content (%)

Soil.N.min: Soil net nitrogen mineralisation (kg ha-1 season)

Soil.resp: Soil carbon dioxide emissions (soil respiration; g m-2 h-1)

format: XLS

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

Paper Citation:
 > Risch AC, Ochoa-Hueso R, van der Putten WH, Bump JK, Busse MD, Frey B, Gwiazdowicz DJ, Page-Dumroese DS, Vandegehuchte ML, Zimmermann S, Schütz M. Size-dependent loss of aboveground animals differentially affects grassland ecosystem coupling and functions. 2018. Nature Communications 9: 3684. [doi: 10.1038/s41467-018-06105-4](https://doi.org/10.1038/s41467-018-06105-4).

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

#Methods
Study sites
The experimental exclosure setups were installed within the SNP (IUCN category Ia preserve; Dudley 2008), in south-eastern Switzerland. The park covers 172 km2 of forests and subalpine and alpine grasslands along with scattered rock outcrops and scree slopes. The entire area has been protected from human impact (no hunting, fishing, camping or off-trail hiking) since 1914. Large, fairly homogenous patches of short- and tall-grass vegetation, which originate from different historical management and grazing regimes, cover the park’s subalpine grasslands entirely. Short-grass vegetation developed in areas where cattle used to rest (nutrient input) prior to the park’s foundation (14th century to 1914) (Schütz and others 2003, 2006) and is dominated by lawn grass species such as Festuca rubra L., Briza media L. and Agrostis capillaris L. (Schütz and others 2003, 2006). Today, this vegetation type is intensively grazed by diverse vertebrate and invertebrate communities that inhabit the park and consume up to 60% of the available biomass (Risch and others 2013). Tall-grass vegetation developed where cattle formerly grazed, but did not rest, and is dominated by rather nutrient-poor tussocks of Carex sempervirens Vill. and Nardus stricta L. (Schütz and others 2003, 2006). This vegetation type receives considerably less grazing, with only roughly 20% of the biomass consumed (Risch and others 2013). Consequently, the two vegetation types together represent a long-term trajectory of changes in grazing regimes. Underlying bedrock of all grasslands is dolomite, which renders these grasslands rather poor in nutrients regardless of former and current land-use regimes.
Experimental design
To progressively exclude aboveground vertebrate and invertebrate animals, we established 18 size-selective exclosure setups (nine in short-grass, nine in tall-grass vegetation) distributed over six subalpine grasslands across the SNP (Risch and others 2013, 2015). Elevation differences of exclosure locations did not exceed 350 m (between 1975 and 2300 m a.s.l.). The exclosures were established immediately after snowmelt in spring 2009 and were left in place for five consecutive growing seasons (until end of 2013). They were, however, temporarily dismantled every fall (late October after first snowfall) to protect them from avalanches. They were re-established in the same location every spring immediately after snowmelt. Each size-selective exclosure setup consisted of five plots (2 x 3 m) that progressively excluded aboveground vertebrates and invertebrates from large to small. The plots are labelled according to the guilds that had access to them “L/M/S/I”, “M/S/I”, “S/I”, “I”, “None”; L = large mammals, M = medium mammals, S = small mammals, I = invertebrates, None = no animals had access. As we only had permission to have the experimental setup in place for five consecutive growing seasons, the experiment had to be completely dismantled in the late fall of 2013 and all material removed from the SNP.
Our exclosure design was aimed at excluding mammalian herbivores, but naturally also excluded the few medium and small mammalian predators, as well as the entire aboveground invertebrate food web. A total of 26 large to small mammal species can be found in the SNP, but large apex predators are missing (wolf, bear, lynx) . Reptiles, amphibians and birds are scarce to absent in the subalpine grasslands under study. Only two reptile species occur in the park and they are confined to rocky areas that warm up enough for them to survive. One frog species spawns in an isolated pond far from our grasslands. Only three bird species occasionally feed on the subalpine grasslands. Using game cameras (Moultrie 6MP Game Spy I-60 Infrared Digital Game Camera, Moultrie Feeders, Alabaster, AL, USA), we did observe that the medium- and small-sized mammals (marmot/hares and mice) were not afraid to enter the fences and feed on their designated plots. We never spotted reptiles, amphibians or birds on camera. We distinguished between 59 higher aboveground-dwelling invertebrate taxa that our size-selective exclosures excluded (see also methods for aboveground-dwelling invertebrates below).
The “L/M/S/I” plot (not fenced) was located at least 5 m from the 2.1 m tall and 7 x 9 m large main electrical fence that enclosed the other four plots. The bottom wire of this fence was mounted at 0.5 m height and was not electrified to enable safe access for medium and small mammals, while fencing out the large ones. Within each main fence, we randomly established four 2 x 3 m plots separated by 1-m wide walkways from one another and from the main fence line: 1) the “M/S/I” plots were unfenced, allowing access to all but the large mammals; 2) the “S/I” plots (10 x 10 cm electrical mesh fence) excluded all medium-sized mammals. Note that the bottom 10 cm of this fence remained non-electrified to enable safe access for small mammals; 3) the “I” plots (2 x 2 cm metal mesh fence) excluded all mammals. We double-folded the mesh at the bottom 50 cm to reduce the mesh size to smaller than 1 x 1 cm openings; and 4) the “None” plots were surrounded by a 1 m tall mosquito net (1.5 x 2 mm) to exclude all animals. The top of the plot was covered with a mosquito-meshed wooden frame mounted to the corner posts (roof). We treated these plots a few times with biocompatible insecticide (Clean kill original, Eco Belle GmbH, Waldshut-Tiengen, Germany) to remove insects that might have entered during data collection or that hatched from the soil, but amounts were negligible and did not impact soil moisture conditions within these plots.
To assess whether the design of the “None” exclosure (mesh and roof) affected the response variables within the plots and, therefore, influenced the results, we established an additional six “micro-climate control” exclosures (one in each of the six grasslands) (Risch and others 2013, 2015). These exclosures were built as the “None” exclosures but were open at the bottom (20 cm) of the 3-m side of the fence facing away from the prevailing wind direction to allow invertebrates to enter. A 20-cm high and 3-m long strip of metal mesh was used to block access to small mammals. Thus, this construction allowed a comparable micro-climate to the “None” plots, but also a comparable feeding pressure by invertebrates to the “I” plots. We compared various properties within these exclosures against one another to assess if our construction altered the conditions in the “None” plots. We showed that differences in plant (e.g., vegetation height, aboveground biomass) and soil properties (e.g., soil temperature, moisture) found between the “I” and the “None” treatments were not due to the construction of the “None” exclosure, but a function of animal exclusions, although the amount of UV light reaching the plant canopy was significantly reduced (Risch and others 2013).
Aboveground invertebrate sampling
Aboveground invertebrates were sampled with two different methods to capture both ground- and plant-dwelling organisms: 1) we randomly placed two pitfall traps (67 mm in diameter, covered with a roof) filled with 20% propylene glycol in one 1 x 1 m subplot of the 2 x 3 m treatment plots in spring 2013 (May) and emptied them every two weeks until late September 2013 (Vandegehuchte and others 2017b, 2018). A pitfall trap consisted of a plastic cylinder (13 cm depth, 6.75 cm diameter). Within each cylinder we placed a 100 ml plastic vial with outer diameter 6.70 cm and on top of the cylinder we placed a plastic funnel to guide the invertebrates into the vials. Each trap was cover with a cone-shaped and transparent plastic roof to protect the trap from rain (Vandegehuchte and others 2017b, 2018). Note that in the “None” plots only one trap was placed as control to check for effectiveness of the exclosure. 2) We vacuumed all invertebrates from a 60 x 60 cm area on another 1 x 1 m subplot with a suction sampler (Vortis, Burkhard manufacturing CO, Ltd., Rickmansworth, Hertfordshire, UK) every month from June to September 2013 (Vandegehuchte and others 2017b, 2018). For this purpose, we quickly placed a square plastic frame (60 x 60 x 40 cm) with a closable mosquito mesh sleeve attached to the top edge into the plot from the outside. The suction sample was then inserted into through the sleeve and operated for 45 s to collect the invertebrates (Vandegehuchte and others 2017b, 2018).
We sorted the ≈100 000 individuals collected with both methods by hand and identified each individual morphologically to the lowest taxonomic level feasible (59 taxa, including orders, suborders, subfamilies, families; phylum for Mollusca). These taxa belonged to the following feeding types: 19 herbivores, 16 detritivores, 9 predators, 8 mixed feeders, 5 omnivores and 2 non-classified feeders (or not feeding as adults) (Vandegehuchte and others 2017b). We summed the numbers from the two pitfall traps and the suction sampling over the course of the 2013 season to represent the aboveground invertebrate abundance and community composition of a plot. Note: we did not specifically attempt to catch flying invertebrates with e.g., sticky traps, thus a few flying insects may have been missed with our vacuum sampling approach.
Sampling of plant properties
The vascular plant species composition was assessed at peak biomass every summer (July) by estimating the frequency of occurrence of each species with the pin count method in each plot (Frank and McNaughton 1990). A total of 172 taxa occurred within our 90 plots and we calculated plant species richness for each plot separately. We used the 2013 data in this study. Plant quality was assessed every year in July and September; here we use plant quality at the end of the experiment (September 2013). Two 10 x 100 cm wide strips of vegetation per plot were clipped, combined, dried at 65°C, and ground (Pulverisette 16, Fritsch, Idar-Oberstein, Germany) to pass through a 0.5 mm sieve. Twenty randomly selected samples across all treatments were analysed for N (Leco TruSpec Analyser, Leco, St. Joseph, Michigan, USA) (Vandegehuchte and others 2015). Nitrogen concentrations of the other samples were then estimated from models established for the experiment and the entire SNP relating Fourier transform-near infrared reflectance (FT-NIR) spectra to the measured values of N using a multi-purpose FT-NIR spectrometer (Bruker Optics, Fällanden, Switzerland) (Vandegehuchte and others 2015). Root biomass was sampled every fall by collecting five 2.2 cm diameter x 10 cm deep soil samples (Giddings Machine Company, Windsor, CO, USA) per plot (450 samples year-1). The samples were dried at 30 °C and roots were sorted from the sample by hand. We sorted each sample for 1 h which allowed to retrieve over 90% of all roots present in the samples (Risch and others 2013). The roots were then dried at 65 °C for 48 and weighed to the nearest mg. We averaged the values per plot and used the 2013 data only in this study.
Sampling of edaphic communities
In 2009, 2010, and 2011 we collected three composited soil samples (5 cm diameter x 10 cm depth; AMS Samplers, American Falls, ID, USA) and assessed bacterial community structure using T-RFLP profiling (Liu and others 1997; Blackwood and others 2003; Hodel and others 2014). We detected a total of 89 operational taxonomic units (OTUs). These values are in accordance with other studies reporting OTU richness (Wirthner and others 2011; Zumsteg and others 2012; Meola and others 2014) using T-RFLP profiling, a method that detects the most abundant, and thus likely, the most relevant, taxa. We averaged the data over the three years of collections for our calculations. Microbial biomass carbon (MBC) was determined with the substrate-induced method (Anderson and Domsch 1978) every fall (September) between 2009 and 2013 by collecting three mineral soil samples (5 cm diameter × 10 cm mineral soil core, AMS Samplers, American Falls, ID, USA). The three samples were combined (90 samples for each sampling year), immediately put on ice, taken to the laboratory, passed through a 2-mm sieve and stored at 4°C. Again, we only used the 2013 data in this study.
Soil samples (5 cm diameter x 10 cm depth) to extract soil arthropods were collected in June, July, and August 2011 with a soil corer lined with a plastic sleeve to ensure an undisturbed sample (total of 270 samples). The plastic line core was immediately sealed on both ends using cling film and put into a cooler. All plots were sampled within three days and the extraction of arthropods started the evening of the sampling day using a high-gradient Tullgren funnel apparatus (Crossley and Blair 1991; Vandegehuchte and others 2015). Samples were kept in the extractor for four days and the soil arthropods were collected in 95% ethanol. All individuals were counted and each individual was identified morphologically to the lowest level feasible [76 taxa, including orders, suborders, subfamilies, families (Protura, Thysanoptera, Aphidina, Psylina, Coleoptera, Brachycera, Nematocera, Auchenorryncha, Heteroptera, Formicidae); sub-phylum for Myriapoda, for Acari and Collembola also including morpho-species). Note that we also included larval stages (nine of the 76 taxa) (Vandegehuchte and others 2015). All data were summed over the season. A detailed species list for mites and collembolans is published (Vandegehuchte and others 2017a) [https://doi.org/10.1371/journal.pone.0118679.s001]. Earthworms are rare in the SNP and therefore were not included. We collected eight random 2.2 cm diameter x 10 cm deep soil cores from each plot in September 2013 to determine the soil nematode community composition. The samples were mixed and the nematodes were extracted from 100 ml of fresh soil using Oostenbrink elutriators (Oostenbrink 1960). All nematodes in a 1 ml of the 10 ml extract were counted, a minimum of 150 individuals sample-1 were identified to genus or family level using (Bongers 1988), the numbers of all nematodes were extrapolated to the entire sample and expressed for a 100 g dry sample. In total we identified 63 genus or family levels (Vandegehuchte and others 2015). The list of all the nematodes found is published (Vandegehuchte and others 2015) [http://www.oikosjournal.org/appendix/oik-03341] or DOI: [doi: 10.1111/oik.03341].
We are aware that sampling soil microbes from 2009 to 2011 and soil arthropods in 2011 was not ideal, but we are positive that this does not bias the results. Most of the parameters measured in our experiment either already showed a treatment response after the first growing season (e.g., plant biomass) or did not respond over the entire time experiment (e.g., microbial biomass C). The microbial community composition (2009 – 2011) was highly influenced by inter-annual differences in temperature and precipitation, but did not differ between treatments or vegetation types (Hodel and others 2014). We therefore felt comfortable using the 2009 through 2011 data for describing the soil microbial community in our experimental treatments. Similarly, we are positive that our soil arthropod data are representative. We did assess soil arthropods in August 2012 and found no differences to the August 2011 data. However, we did not feel comfortable combining the 2011 June, July, August data with only August data for 2012 for our analyses.
Sampling of soil properties
We collected three soil samples (5 cm diameter x 10 cm depth) in each plot in September 2013 after removing the vegetation. First, we collected the top layer of mineral soil rich in organic matter, the surface organic layer or rhizosphere, typically 1 to 3 cm in depth with a soil corer (AMS Samples, American Falls, Idaho, USA). Second, we collected a 10 cm mineral soil core beneath this surface layer. The cores for each layer were composited, dried at 65 °C for 48 h and fine-ground to pass a 0.5 mm screen. We then analysed all samples for total C using a Leco TruSpec Analyser (Leco, St. Joseph, Michigan, USA). Mineral soil pH was measured potentiometrically in 1:2 soil:CaCl2 solution with an equilibration time of 30 min.
Soil net N mineralisation was assessed during the 2013 growing season (Risch and others 2015). For this purpose, we randomly collected a 5 cm diameter x10 cm deep soil sample with a soil corer (AMS Samples, American Falls, Idaho, USA) after clipping the vegetation in June 2013. After weighing and sieving (4 mm mesh) the soil, we extracted a 20 g subsample in 1 mol l-1 KCl for 1.5 h on an end-over-end shaker and thereafter filtered it through ashless folded filter paper (DF 5895 150, ALBET LabScience, Hahnenmühle FineArt GmbH, Dassel, Germany). From these filtrates NO3- concentrations were measured colorimetrically (Norman and Stucki 1981) and NH4+with flow injection analysis (FIAS 300, Perkin Elmer, Waltham Massachusetts, USA) (Risch and others 2015). We dried the rest of the sample 105 °C to constant mass to determine fine,fraction bulk density. A second soil sample was collected within each plot in June 2013 with a corer lined with a 5 x 13 cm aluminium cylinder. The corer was driven 11.5 cm deep into the soil so that the top 1.5 cm of the cylinder remained empty. Into this space we placed a polyester bag (250 µm) filled an ion-exchanger resin to capture the incoming N. The bag was filled with a 1:1 mixture of acidic and alkaline exchanger resin (ion-exchanger I KA/ion exchanger IIIAA, Merck AG, Darmstadt, Germany). We then removed 1.5 cm soil at the bottom of the cylinder and placed a second resin exchanger bag into this space to capture the N leached from the soil column. To assure that the exchange resin was saturated with H+ and Cl- prior to filling the bags, the mixture was stirred with 1.2 ml l-1 HCl for 1 h and then rinsed with demineralized water until the electrical conductivity of the water reached 5 µm cm-1. The cylinder with the resin bags in place was reinserted into the soil with the top flush to the soil surface and incubated for three months. We recollected the cylinders in September 2013. Each resin bag and 20 g of sieved soil (4 mm mesh) from each cylinder were then separately extracted with KCl and NO3- and NH4+ concentrations were measured. Nitrate and NH4+ concentrations of all samples were then converted to a content basis by multiplying their values with fine fraction bulk density. Net N mineralisation was thereafter calculated as the difference between the N content of the samples collected at the end of the three-month incubation (including the N extracted from the bottom resin bag) and the N content at the beginning of the incubation (Risch and others 2015).
Soil CO2 emissions were measured every two weeks between 0900 and 1700 hrs from early May through late September 2013 with a PP-Systems SRC-1 soil respiration chamber (15 cm high, 10 cm diameter; closed circuit) attached to a PP-Systems EGM-4 infrared gas analyser (PP-Systems, Amesbury, MA, USA) on two locations per plot (Risch and others 2013). The chamber was placed on randomly placed, permanently installed PVC collars (10 cm diameter) driven 5 cm into the soil at the beginning of the study (Risch and others 2013). Freshly germinated plants growing within the collars were removed prior to each measurement to avoid measuring plant respiration or photosynthesis. The two measurements collected per plot and sampling date were averaged.
Soil moisture (with time domain reflectometry; Field-Scout TDR-100, Spectrum Technologies, Plainfield, Illionois, USA) and temperature (with a waterproof digital pocket thermometer; Barnstead International, Dubuque, Iowa, USA) were measured at five random locations per plot every two weeks during the growing seasons during the experiment for the 0 to 10 cm depth (Risch and others 2013, 2015). As soil moisture and soil temperature were highly negatively correlated (Risch and others 2013), we only used soil moisture for this study. We used plot-level averages of all values available to capture soil moisture variability during the five years of the experiment. The results remained unchanged when we only used soil moisture from the 2013 growing season.
Numeral calculations and statistical analyses
Ecosystem coupling. We conducted principal component analyses (PCAs; unscaled) at the complete dataset level using the abundances of each taxonomical entity to describe each of the five different communities used in this study: aboveground-dwelling invertebrates, vascular plants, soil microorganisms, soil arthropods and soil nematodes. We retained the first two components (PCA axis 1 and PCA axis 2) of each analysis as we found them to adequately represent the temporal and spatial variability of our 90 treatment plots in previous studies55,67. Together they explained a total of 71.70% of the variation for aboveground invertebrates, 44.36% for plants, 44.85% for soil microorganisms, 61.85% for soil arthropods and 77.19% for soil nematodes. In addition, we used soil pH and soil organic C content as a proxy for soil chemical properties, soil bulk density as a proxy for soil physical properties and soil moisture (negatively correlated with soil temperature) as a proxy for soil micro-climatic conditions for an overall total of fourteen constituents.
We calculated ecosystem coupling9 for each exclosure treatment within each vegetation type (i.e., 2  5 treatment combinations in total) as an integrated measure of pairwise ecological interactions between ecosystem constituents representing ecological communities and the soil abiotic environment. These ecological interactions are defined by non-parametric Spearman rank correlation analyses between two constituents, excluding interactions involving two abiotic constituents (e.g., soil pH vs. soil moisture) and interactions between the first (PC1) and second (PC2) component of each community type, as these are orthogonal by definition. Interactions between abiotic constituents were excluded from the analyses because the focus of our study was on communities and how they interact with one another and their surrounding environment; therefore, including abiotic-abiotic interactions was not of interest here. Given that the effectiveness of our experimental design resulted in that no community composition data of aboveground-dwelling invertebrates was available for the “None” plots (all animals excluded), only thirteen instead of fourteen constituents were included in the ecosystem coupling calculations for this treatment. The complete absence of aboveground invertebrates represents the most extreme case of disturbance between aboveground animal communities and the rest of the ecosystem constituents. This may have resulted in a slight overestimation of ecosystem coupling for these plots.
	Average ecosystem coupling was calculated as follows:
Ecosystem coupling=
where Xi is the absolute Coupling was calculated value of the Spearman’s rho coefficient of the ith correlation for each treatment within each vegetation type (i.e., based on nine replicates each), considering and n is the number of pairwise comparisons considered (n = a total of 80; interactions (56 in the case of the “None” treatment). We considered a total of 40 biotic-biotic interactions (i.e., concerning two community-level principal components such as plants and microbes; 24 in the case of the “None” treatment) and 40 abiotic-biotic (i.e., concerning one community-level principal component and one abiotic factor, e.g., plant community and soil properties; 32 in the case of the “None” treatment).	Coupling was calculated for each treatment within each vegetation type (i.e., based on nine replicates each), considering a total of 80 interactions (56 in the case of the “None” treatment). We considered a total of 40 biotic-biotic interactions (i.e., concerning two community-level principal components such as plants and microbes; 24 in the case of the “None” treatment) and 40 abiotic-biotic (i.e., concerning one community-level principal component and one abiotic factor, e.g., plant community and soil properties; 32 in the case of the “None” treatment).
To establish whether constituents were significantly and positively coupled within treatments (i.e., the average of their correlation coefficients were greater than in a null model where correlation only happens by chance), we calculated one-tailed p-values based on permutation tests with 999 permutations.

We considered six ecosystem functions and process rates commonly used to assess ecosystem functioning (Meyer and others 2015; Manning and others 2018). Plant N content represents a measure of forage quality, while plant richness has been shown to stabilise biomass production, thus allowing the system to respond to changes in herbivory. Soil net N mineralisation, soil respiration, root biomass, and microbial biomass represent fluxes or stocks of energy. For all functions and processes higher values represent higher functioning (Manning and others 2018). All these variables were measured in the last year of the experiment (2013). We then quantified ecosystem multifunctionality using the multiple threshold approach (Byrnes and others 2014; Manning and others 2018), which considers the number of functions that are above a certain threshold, over a series of threshold values (typically 10-99%) that are defined based on the maximum value of each function. We weighted all our functions equally for these calculations (Manning and others 2018). The number of functions in a plot with values higher than a given threshold value for the respective function is summed up. The sum represents ecosystem multifunctionality for that plot. Given that choosing any particular threshold as a measure of ecosystem multifunctionality is arbitrary, we calculated the average of thresholds from 10-90% (in 10% intervals) as a more integrated representation of ecosystem multifunctionality.
We used Pearson correlations to explore the relationships between ecosystem coupling (all interactions, biotic-biotic interactions, abiotic-biotic interactions involving above- and belowground constituents, and all interactions, biotic-biotic interactions, abiotic-biotic interactions involving belowground constituents only) and ecosystem multifunctionality by calculating the slopes of all relationships between ecosystem coupling and multifunctionality for all thresholds between 10 and 99%. We also related ecosystem coupling with the average of multifunctionality at thresholds between 30-80% as explained before and considered this correlation as a robust indication of the type of association between these two variables. In addition, we explored the relationships between ecosystem coupling (all interactions, biotic-biotic interactions, abiotic-biotic interactions involving above- and belowground constituents, and all interactions, biotic-biotic interactions, abiotic-biotic interactions involving belowground constituents only) and individual ecosystem functions. The effects of exclosures and vegetation type on individual functions and multifunctionality were evaluated using linear mixed effects models ('lme' function of the nlme package), with exclosure and vegetation type as fixed effects and fence as a random factor. All statistical analyses and numerical calculations were done in R version 3.4.0 (R Core Team 2016).

#References

- Anderson J, Domsch K. 1978. A physiological method for the quantitative measurement of microbial biomass in soil. Soil Biol Biochem 10:215–21.
- Blackwood CB, Marsh T, Kim S-H, Paul EA. 2003. Terminal Restriction Fragment Length Polymorphism Data Analysis for Quantitative Comparison of Microbial Communities. Appl Environ Microbiol 69:926–32. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC143601/
- Bongers T. 1988. De nematoden von Nederland. Schoorl, The Netherlands: Pirola
- Byrnes JEK, Gamfeldt L, Isbell F, Lefcheck JS, Griffin JN, Hector A, Cardinale BJ, Hooper DU, Dee LE, Duffy JE. 2014. Investigating the relationship between biodiversity and ecosystem multifunctionality: Challenges and solutions. Methods Ecol Evol 5:111–24.
- Crossley DAJ, Blair JM. 1991. A high-efficiency low-technology Tulgren-type extractor for soil microarthopods. Agric Ecosyst Environ 34:187–92.
- Dudley N. 2008. Guidelines for applying protected area managment categories. Gland: IUCN
- Frank DA, McNaughton SJ. 1990. Aboveground biomass estimation with the canopy intercept method: A plant growth form caveat. Oikos 57:57–60.
- Haynes AG, Schütz M, Buchmann N, Page-Dumroese DS, Busse MD, Risch AC. 2014. Linkages between grazing history and herbivore exclusion on decomposition rates in mineral soils of subalpine grasslands. Plant Soil 374.
- Hodel M, Schütz M, Vandegehuchte ML, Frey B, Albrecht M, Busse MD, Risch AC. 2014. Does the aboveground herbivore assemblage influence soil bacterial community composition and richness in subalpine grasslands? Microb Ecol 68:584–95.
- Liu WT, Marsh TL, Cheng H, Forney LJ. 1997. Characterization of microbial diversity by determining terminal restriction fragment length polymorphisms of genes encoding 16S rRNA. Appl Environ Microbiol 63:4516–22. http://www.ncbi.nlm.nih.gov/pmc/articles/PMC168770/
- Manning P, van der Plas F, Soliveres S, Allan E, Maestre FT, Mace G, Whittingham MJ, Fischer M. 2018. Redefining ecosystem multifunctionality. Nat Ecol Evol 2:427–36. https://doi.org/10.1038/s41559-017-0461-7
- Meola M, Lazzaro A, Zeyer J. 2014. Diversity, resistance and resilience of the bacterial communities at two alpine glacier forefields after a reciprocal soil transplantation. Environ Microbiol 16:1918–34. https://onlinelibrary.wiley.com/doi/abs/10.1111/1462-2920.12435
- Meyer ST, Koch C, Weisser WW. 2015. Towards a standardized Rapid Ecosystem Function Assessment (REFA). Trends Ecol Evol 30:390–7. http://www.sciencedirect.com/science/article/pii/S0169534715000968
- Norman R., Stucki JW. 1981. The determination of nitrate and nitrite in soil extracts by ultraviolet spectrophotometry. Soil Sci Soc Am J 45:347–53.
- Ochoa-Hueso R. 2016. Non-linear disruption of ecological interactions in response to nitrogen deposition. Ecology 87:2802–2814.
- Oostenbrink M. 1960. Estimating nematode populations by some selected methods. In: Sasser NJ, Jenkins WR, editors. Nematology. Chapel Hill, NC, USA: University of North Carolina Press. pp 85–101.
- R Core Team. 2016. R: A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing
- Risch AC, Haynes AG, Busse MD, Filli F, Schütz M. 2013. The response of soil CO2 fluxes to progressively excluding vertebrate and invertebrate herbivores depends on ecosystem type. Ecosystems 16:1192–202.
- Risch AC, Schütz M, Vandegehuchte ML, Van Der Putten WH, Duyts H, Raschein U, Gwiazdowicz DJ, Busse MD, Page-Dumroese DS, Zimmermann S. 2015. Aboveground vertebrate and invertebrate herbivore impact on net N mineralization in subalpine grasslands. Ecology 96:3312–22.
- Schütz M, Risch AC, Achermann G, Thiel-Egenter C, Page-Dumroese DS, Jurgensen MF, Edwards PJ. 2006. Phosphorus translocation by red deer on a subalpine grassland in the Central European Alps. Ecosystems 9:624–633.
- Schütz M, Risch AC, Leuzinger E, Krüsi BO, Achermann G. 2003. Impact of herbivory by red deer (Cervus elaphus L.) on patterns and processes in subalpine grasslands in the Swiss National Park. For Ecol Manage 181:177–88.
- Vandegehuchte ML, van der Putten WH, Duyts H, Schütz M, Risch AC. 2017a. Aboveground mammal and invertebrate exclusions cause consistent changes in soil food webs of two subalpine grassland types, but mechanisms are system-speciﬁc. Oikos 126:212–23.
- Vandegehuchte ML, Raschein U, Schütz M, Gwiazdowicz DJ, Risch AC. 2015. Indirect short- and long-term effects of aboveground invertebrate and vertebrate herbivores on soil microarthropod communities. PLoS One 10:e0118679.
- Vandegehuchte ML, Schütz M, de Schaetzen F, Risch AC. 2017b. Mammal-induced trophic cascades in invertebrate food webs are modulated by grazing intensity in subalpine grassland. J Anim Ecol 86:1434–46.
- Vandegehuchte ML, Trivellone V, Schütz M, Firn J, de Schaetzen F, Risch AC. 2018. Mammalian herbivores affect leafhoppers associated with specific plant functional types at different timescales. Funct Ecol 32:545–55.
- Wirthner S, Frey B, Busse MD, Schütz M, Risch AC. 2011. Effects of wild boar (Sus scrofa L.) rooting on the bacterial community structure in mixed-hardwood forest soils in Switzerland. Eur J Soil Biol 47:296–302. http://dx.doi.org/10.1016/j.ejsobi.2011.07.003
- Zumsteg A, Luster J, Göransson H, Smittenberg RH, Brunner I, Bernasconi SM, Zeyer J, Frey B. 2012. Bacterial, Archaeal and Fungal Succession in the Forefield of a Receding Glacier. Microb Ecol 63:552–64. https://doi.org/10.1007/s00248-011-9991-8

ecosystem-coupling

display_name: ecosystem-coupling

ecosystem-functions

display_name: ecosystem-functions

ecosystem-multifunctionality

display_name: ecosystem-multifunctionality

exclusion

display_name: exclusion

grassland

display_name: grassland

grazing

display_name: grazing

invertebrates

display_name: invertebrates

subalpine

display_name: subalpine

swiss-national-park

display_name: swiss-national-park

vertebrates

display_name: vertebrates

educ

en: Education, culture and sport

en: Ecosystem coupling and multifunctionality - exclosure experiment


maintainer: Esther Frei

en: treeline mountain-forest stillberg afforestation long-term-experiment

mimetype: application/x-zip-compressed

en:
 Background information

The Stillberg ecological treeline research site in the Swiss Alps was established in 1975, with the aim to develop ecologically, technically, and economically sustainable reforestation 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. Long-term monitoring of the large-scale high-elevation afforestation has generated data about tree growth, survival, and vitality. In addition, detailed characteristics of the microsite conditions of the research were conducted. Besides providing a scientific basis and practical guidelines for high-elevation afforestation, this research has contributed to a comprehensive understanding of ecological processes in the treeline ecotone.

 Experiment description

The Stillberg afforestation experiment was established in 1975 by planting 92,000 seedlings of *Larix decidua*, *Pinus cembra* and *Pinus mugo* ssp. *uncinata* in the alpine treeline ecotone. The afforestation site is located on a northeast-facing slope with steep, topographically highly structured terrain and covers elevations from 2075 to 2230 m a.s.l. The afforestation site was divided into 4052 square plots of 3.5 × 3.5 m, arranged in a regular species-alternating pattern over the whole area. Each plot contained 25 trees of one species (1350 plots per species), and the seedlings were systematically planted 70 cm apart. The trees have been monitored since 1975. Specifically, tree mortality was assessed annually from 1975 until 1995 and has been documented every ten years since then, with surveys in 2005 and 2015 (the next survey is due in 2025). Height of the surviving trees was measured in 1975, 1979, 1982, 1985, 1990, 1995, 2005, and 2015. In 1995, 2005, and 2015, drivers of tree vitality were assessed for a subset of trees per plot. Additionally, an extensive set of environmental parameters characterizing microsite conditions of the afforestation area were recorded before and after the planting of the trees.

 Data description

The five datasets from the afforestation experiment comprise ecological and environmental data from the main afforestation experiment in five datasets with accompanying metadata (Stillberg_afforestation_all_metadata.xlsx). All data and metadata files are bundled in a ZIP-file (Stillberg_afforestation_v1.zip). In particular, a first dataset contains environmental data characterising microsite conditions of the 4000 plots with regard to soil, topography, vegetation and microclimatic conditions (Stillberg_afforestation_plot_data_v1.csv; Stillberg_afforestation_plot_metadata_v1.csv. In each plot, the natural tree regeneration was assessed by counting seedings of several tree species in 2005 and 2015 (Stillberg_afforestation_regeneration_data_v1.csv; Stillberg_afforestation_regeneration_metadata_v1.csv). Furthermore, specific information about each of the 92’000 planted trees of the tree species is available (Stillberg_afforestation_tree_parameter_data_v1.csv; Stillberg_afforestation_tree_parameter_metadata_v1.csv). Survival data for each of the 92’000 individual trees can be found in a separate dataset (Stillberg_afforestation_tree_survival_data_v1.csv; Stillberg_afforestation_tree_survival_metadata_v1.csv). Tree growth and vitality parameters are available for all trees from 1995, and for subsets of trees for 2005 and 2015 (Stillberg_afforestation_tree_measurements_data_v1.csv; Stillberg_afforestation_tree_measurements_metadata_v1.csv).

afforestation

display_name: afforestation

long-term-experiment

display_name: long-term-experiment

mountain-forest

display_name: mountain-forest

stillberg

display_name: stillberg

treeline

display_name: treeline

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: Long-term afforestation experiment at the Alpine treeline site Stillberg, Switzerland


tree-species-map-of-switzerland

maintainer: Tiziana Li Koch

title_for_slug: tree-species-map-of-switzerland

en: tree-species map remote-sensing forest switzerland

media_type: image/png

en: Image of the tree species map

format: PNG

en:
Dominant tree species map of Switzerland
We created a tree species map of Switzerland for the dominant tree species in the forested areas. The spatial resolution of the map is 10 m and the coordinate system is ETRS89-extended / LAEA Europe (EPSG 3035). The map comprises Sentinel-2 index time series from the year 2020, a digital elevation model and species reference data from the Swiss National Forest Inventory. The map is available as raster (.tif) or vector dataset (.gpkg).
**Access will be granted upon request.**
In total, the following 15 species were mapped:
*Abies alba*, *Acer pseudoplatanus*, *Alnus glutinosa*, *Alnus incana*, *Betula pendula*, *Castanea sativa*, *Fagus sylvatica*, *Fraxinus excelsior*, *Picea abies*, *Pinus cembra*, *Pinus mugo arborea*, *Pinus sylvestris*, *Quercus petraea*, *Quercus robur*, *Sorbus aucuparia*.

 Approach

 Data
- Swiss National Forest Inventory Data (stand species with > 60 % dominance in upper canopy; on at least more than 9 plots dominant)
- Sentinel-2 time series (2020, Indices: CCI, CIRE, NDMI, EVI, NDVI)
- Digital elevation model (DEM) (swissalti3d, 5 m)
- Biogeographical regions (Federal Office for the Environment FOEN)
- Forest mask 2017 (Approach: Waser et al., 2015)

 Modeling approach
We identified the most meaningful variables that led to separation of the respective groups by using random forest models with a forward feature selection (Meyer et al., 2018; Ververidis & Kotropoulos, 2005). In this approach, the final random forest model is solely built from the selected meaningful variables. By identifying meaningful variables, we can determine which variables might influence the grouping. Further, to avoid overfitting and overly optimistic results, we applied 10-fold spatial cross-validation and put all pixels from a plot in the same spatial fold.
The modeling was realized using the CAST package in R (Meyer et al., 2022), based on the well-known caret package (Kuhn, 2022). We used the ranger package in R (Wright & Ziegler, 2017) to implement the random forest models, due to its short computation time.

 Training data for modeling
- 295 Sentinel-2, DEM & Biogeographical variables
- 10525 tree species pixels

 Selected variables for final model
1. EVI of 2020.05.16
2. NDMI of 2020.03.12
3. CIRE of 2020.04.16
4. NDMI of 2020.07.05
5. CCI of 2020.05.11
6. dem
7. CCI of 2020.08.14
8. NDMI of 2020.08.24
9. CCI of 2020.12.22
10. NDMI of 2020.04.21
11. NDMI of 2020.11.17
12. NDMI of 2020.08.09
13. CIRE of 2020.03.22
14. CIRE of 2020.08.09
14. CCI of 2020.11.02
15. CIRE of 2020.06.10

 Overall Accuracy of final model
- 0.759

 Nationwide prediction
- Predicted throughout forest mask 2017 (Approach: Waser et al., 2015)
- Not applied on incomplete Sentinel-2 time series (own category in final map: incomplete_ts)
- Applied the Area of Applicability (Meyer 2022) to sort out pixels outside of the feature space; basically where the model had not the same values for pixels as in the available training data

 *Be aware that the map is only validated with the training data itself, an independent validation with other data sources remains missing*

 References
- Kuhn, M. (2022). Classification and Regression Training. 6.0-93.
- Meyer, H., Reudenbach, C., Hengl, T., Katurji, M., & Nauss, T. (2018). *Improving performance of spatio-temporal machine learning models using forward feature selection and target-oriented validation*. Environmental Modelling and Software, 101, 1-9. https://doi.org/10.1016/j.envsoft.2017.12.001
- Meyer, H., Milà, C., & Ludwig, M. (2022). *CAST: 'caret' Applications for Spatial-Temporal Models*. 0.7.0.
- Ververidis, D., & Kotropoulos, C. (2005). *Sequential forward feature selection with low computational cost*. 2005 13th European Signal Processing Conference.
- Waser, L., Fischer, C.,Wang, Z., & Ginzler, C. (2015). *Wall-to-Wall Forest Mapping Based on Digital Surface Models from Image-Based Point Clouds and a NFI Forest Definition*. Forests, 6, 12, 4510–4528.
- Wright, M. N., & Ziegler, A. (2017). *ranger: A Fast Implementation of Random Forests for High Dimensional Data in C++ and R*. Journal of Statistical Software, 77(1), 1-17. https://doi.org/doi:10.18637/jss.v077.i01

forest

display_name: forest

map

display_name: map

remote-sensing

display_name: remote-sensing

switzerland

display_name: switzerland

tree-species

display_name: tree-species

en: Tree species map of Switzerland


maintainer: Virginie Marques

en: fish environmental-dna plankton fjord marine biodiversity dna-metabarcoding

en: Metabarcoding data

en: Raw metabarcoding data for both primers (teleo and 18S) and corresponding file to perform the demultiplexing and primer trimming

en: Code, data and figures

en: Here there are all the data and code to reproduce the figures and analysis presented in the paper, with a README.md file

en:
This dataset contains the raw environmental DNA data associated with the publication *Environmental drivers of eukaryotic plankton and fish biodiversity in an Arctic fjord* in the journal Polar Biology (2023).

 Methods

**Sampling**
We sampled the Lilliehöök fjord on the west coast of Spitsbergen (Svalbard, Norway) over 3 days from 3 to 5 of August 2021. Samples were taken from the glacier front up to the fjord mouth of the Krossfjorden system, around 30 km long, after the Lilliehöök fjord merged with the mouth of Möller fjord. The fjord’s maximum depth has been recorded at 373 m (Svendsen et al. 2002) and has no sill at its entrance, thereby facilitating water exchange with the open ocean of the West Spitsbergen Current. We used a research vessel to sample 5 sites for a total of 15 samples, sampling 3 depths per site (3-m, chlorophyll a maximum and 85-m, unless sea floor was shallower). Shallow and intermediate samples between 3-m and 12-m represent ~35-L of water filtered in-situ using long tubing and a peristaltic pump, and all other deeper samples were taken from a total of 3 Niskin bottles (General Oceanics), representing 22-L of water sampled per sample. Water was filtered through a VigiDNA filtration capsule (SPYGEN) with a 0.20-µm pore size using an Athena peristaltic pump (Proactive Environmental Products, Bradenton, Florida) with a flow rate of ~1-L/min. Each sample was handled with single use tubing and gloves.

**Molecular**
To perform the amplification, we used two sets of primers: teleo (forward: ACACCGCCCGTCACTCT, reverse: CTTCCGGTACACTTACCATG; Valentini et al. 2016) and the universal eukaryotic 1389F/1510R primer pair, amplifying the V9-18S rDNA gene (Amaral-Zettler et al. 2009) (forward: TTGTACACACCGCCC, reverse: CCTTCYGCAGGTTCACCTAC).

 Data content:

+ Metabarcoding data: This zip file contains the 2 sequencing libraries filtered to only retain the samples used in the present study.
+ Code, data and figure: This zip file contains all data and code to reproduce the figures and the analysis in the study, with an associated README explaining the content of each folder.

 Additional informations
For more details, please see the Methods in the associated publication: DOI: 10.1007/s00300-023-03187-9.

biodiversity

display_name: biodiversity

dna-metabarcoding

display_name: dna-metabarcoding

environmental-dna

display_name: environmental-dna

fish

display_name: fish

fjord

display_name: fjord

marine

display_name: marine

plankton

display_name: plankton

en: Data: Environmental drivers of eukaryotic plankton and fish biodiversity in an Arctic fjord


en: temperate-forests forest-dynamics ecological-niche forest-inventory forest tree-mortality mortality-rate

en:
These datasets and script document the analyses and the plot generation for the article **“Trends in background mortality in unmanaged forests across Europe over the last century”**.
We used a comprehensive network of permanent plots with repeated censuses spanning from 1936 to 2020 in 299 unmanaged forest reserves across Europe. The database includes 1.5 million tree records over successive censuses covering a total sampled area of 853.7 ha.
It consists of the following files:
- **Mortality_script.Rmd**: R markdown document to reproduce the analysis and the plot generation.
- **plot_level.RDS**: This dataset contains annual mortality rate values for each permanent plot at each census, of plots with at least three censuses. This dataset is used to calculate the weighted mean annual background mortality rates by forest types, and the series of models to study the trends of annual mortality rates, and then to create Figure 1 and Figure 2.
- **tree_level.RDS**: This dataset contains tree census data for each permanent plot at each census, of plots with at least two censuses. This dataset is used to model species-specific size-dependent mortality probabilities of the studied tree species, and then to create Figure 3.

ecological-niche

display_name: ecological-niche

forest-dynamics

display_name: forest-dynamics

forest-inventory

display_name: forest-inventory

mortality-rate

display_name: mortality-rate

temperate-forests

display_name: temperate-forests

tree-mortality

display_name: tree-mortality

educ

en: Education, culture and sport


c-n-p-release-treeline

maintainer: Frank Hagedorn

title_for_slug: c-n-p-release-treeline

en: enzymatic-activities nutrient-fluxes biogeochemical-cycling treeline-ecosystems litter-decomposition climate-change carbon-cycle microbial-biomass

en: Litter layer characteristics, C-N-P-release during decomposition, microbial C-N-P, enzyme activity, mineralization of 13C-labeled glucose-6-phosphate from samples collected across two treelines

en:
1. Chemical characteristics of litter layer sampled along elevation gradient across treeline ecotones in the Khibiny mountains (birch) and South Urals (spruce)
2. Release of C-N-P during decomposition at 5°C and 15° studied in an incubation experiment with microlysimeters
3. Data on microbial biomass C-N-P, potential activity of C-N-P requiring enzymes in the litter layer
4. Mineralization of 13C labeled glucose-6 phosphate, microbial 13C uptake, and net P immobilization of released P to identify microbial processing of labile organic P

biogeochemical-cycling

display_name: biogeochemical-cycling

carbon-cycle

display_name: carbon-cycle

climate-change

display_name: climate-change

enzymatic-activities

display_name: enzymatic-activities

litter-decomposition

display_name: litter-decomposition

microbial-biomass

display_name: microbial-biomass

nutrient-fluxes

display_name: nutrient-fluxes

treeline-ecosystems

display_name: treeline-ecosystems

envidat

en: EnviDat: the Environmental Data Portal

en: C-N-P-Release Treeline


maintainer: Armin Sigmund

en: snow-sublimation backward-trajectories stable-water-isotopes lagrangian-analysis antarctica isotopic-distillation evaporation

en: Model code

en: Model Sublimation and Model Air Parcel (python scripts)

en: Python scripts for data visualization

en: Scripts for generating the figures of the related research article and generating the data files (csv) from the model output

en: Readme

en: Meta data

en: Data

en: Modeled data including the data shown in the figures of the related research article

en:
 Summary

This data set contains Python programming code and modeled data discussed in a related research article. We developed a simple isotope model to study the drivers of the particularly depleted vapour isotopic composition measured on the ship of the Antarctic Circumnavigation Expedition close to the outlet of the Mertz glacier, East Antarctica, in the 6-day period from 27 January 2017 to 1 February 2017. The model considers the stable water isotopologues H2(16O), H2(18O), and HD(16O). It uses data from the ERA5 reanalysis product with a spatial resolution of 0.25° x 0.25° (Hersbach et al., 2018) and 10-day backward trajectories for the location of the ship, published by Thurnherr et al. (2020a). Our data set includes the model code, Python scripts for visualizing the results, and data produced by the model including the results shown in the figures of the related research article. Here, we summarize the most important model characteristics while further details can be found in the readme.txt file and the related research article including its supporting information.

 Main model characteristics

The modeling approach consists of two steps called *Model Sublimation* and *Model Air Parcel*. The former estimates the isotopic compositions of the snow and sublimation flux across the Antarctic Ice Sheet using an Eulerian frame of reference while the latter models the vapour isotopic composition and specific humidity along air parcel trajectories using a Lagrangian frame of reference. The isotope effects of most phase changes are represented by equilibrium fractionation. Only for ocean evaporation, kinetic fractionation is additionally taken into account (original Craig-Gordon formula). For snow sublimation, two assumptions are tested: *Run E* assumes that sublimation is associated with equilibrium fractionation while *Run N* assumes that sublimation occurs without isotopic fractionation.

 Model Sublimation
Model Sublimation uses a simple one-dimensional mass-balance approach in each grid cell, considering snow accumulation due to snowfall and vapour deposition and snow ablation due to sublimation. The snowpack is represented by 100 layers of equal thickness (e.g., 1 cm) and density (350 kg m-3). The isotopic composition of snowfall is parameterized by generalizing a site-specific, empirical relationship between the daily mean air temperature and snowfall isotopic composition. In the case of vapour deposition, Model Sublimation assumes equilibrium fractionation and estimates the isotopic composition of the atmospheric vapour as the average value for two idealized situations: (i) locally sourced vapour which has the same isotopic composition as the sublimation flux; (ii) non-locally sourced vapour in isotopic equilibrium with snowfall. Model Sublimation is run with a time step of 1 h, independently of Model Air Parcel.

 Model Air Parcel
Every hour, an ensemble of trajectories arrives at different heights in the ABL above the ship. For each of these trajectories, we consider an air parcel with a constant volume of 1 x 1 x 1 m3. The air parcels are initialized at the first suitable time when the trajectories are located in the ABL, either over the ice-free ocean in conditions of evaporation or over snow (Antarctic Ice Sheet or sea ice). Subsequently, the masses of the water isotopologues in the air parcels are simulated with a time step of 3 h, considering vapour uptake or removal due to the moisture flux at the snow or liquid ocean surface (only if the parcel is in the ABL) and cloud/precipitation formation (if the saturation specific humidity is reached). Sea ice is taken into account in a very simplified way. We represent the sea ice by grid cells with a sea-ice cover of more than 90% and assume the isotopic composition of the sublimation flux to be identical to that in the nearest grid cell of the Antarctic Ice Sheet. The isotopic composition of the sublimation flux is taken from Model Sublimation whereas the isotopic composition of the vapour deposition flux (over snow) and condensation flux (over ice-free ocean) is simulated assuming an isotopic equilibrium with the air parcel. Isotope effects of cloud/precipitation formation are represented using the classic Rayleigh distillation model with equilibrium fractionation, where the cloud water is assumed to precipitate immediately after formation.

 References

Hersbach, H., Bell, B., Berrisford, P., Biavati, G., Horanyi, A., Munoz Sabater, J.,... others (2018). *ERA5 hourly data on single levels from 1979 to present. Copernicus Climate Change Service (C3S) Climate Data Store (CDS)*. doi: 10.24381/cds.bd0915c6

Thurnherr, I., Wernli, H., & Aemisegger, F. (2020a). *10-day backward trajectories from ECMWF analysis data along the ship track of the Antarctic Circumnavigation Expedition in austral summer 2016/2017*. Zenodo. doi: 10.5281/zenodo.4031705

antarctica

display_name: antarctica

backward-trajectories

display_name: backward-trajectories

evaporation

display_name: evaporation

isotopic-distillation

display_name: isotopic-distillation

lagrangian-analysis

display_name: lagrangian-analysis

snow-sublimation

display_name: snow-sublimation

stable-water-isotopes

display_name: stable-water-isotopes

envidat

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

en: Modeled Isotopic Composition of Water Vapour Along Air Parcel Trajectories in the Antarctic


sensitivity-of-modeled-snow-instability

maintainer: Bettina Richter

title_for_slug: sensitivity-of-modeled-snow-instability

en: global-sensitivity-analysis snowpack snow-instability numerical-avalanche-forecasting snow-cover-modeling

en: Reference run

en: Reference run for the sensitivity analysis at Weissfluhjoch field site for winter season 2016-2017. Meteorological data from the automatic weather stations are used. No biases were introduced to the input data.

mimetype: application/zip

en: Manually observed snow profiles

en: Manually observed snow profiles at Weissfluhjoch field site for winter season 2016-2017, including layer depth, hand hardness, grain size and grain shape for each layer.

en: Case SL

en: Extracted output data for case SL. Here biases were introduced to meteorological input data during the period of slab formation. Data include quasi random biases, which were introduced, extracted output of simulated snow profiles and sobol indices for each output variable.

en: Case WL

en: Extracted output data for case ALL. Here biases were introduced to meteorological input data during the period of weak layer formation. Data include quasi random biases, which were introduced, extracted output of simulated snow profiles and sobol indices for each output variable.

en: Case ALL

en: Extracted output data for case ALL. Here biases were introduced to meteorological input data during the whole simulation period. Data include quasi random biases, which were introduced, extracted output of simulated snow profiles and sobol indices for each output variable.

en: We investigated the sensitivity of modeled snow instability to meteorological input data using SNOWPACK. We therefore used input data from the automatic weather station at the Weissfluhjoch field site for the year 2016-2017. We investigated three scenarios and performed 14'000 simulations for each scenario. The dataset contains extracted output data from modeled SNOWPACK simulations, including setup files to reproduce the simulations. For further information read the README file.

global-sensitivity-analysis

display_name: global-sensitivity-analysis

numerical-avalanche-forecasting

display_name: numerical-avalanche-forecasting

snow-cover-modeling

display_name: snow-cover-modeling

snow-instability

display_name: snow-instability

snowpack

display_name: snowpack

educ

en: Education, culture and sport

en: Sensitivity of modeled snow instability


en: danger-level snow avalanche-danger snow-instability avalanche avalanche-forecasting

en: Read-me file

en: This read-me file contains the descriptions of the variables.

en: Field observations of snow instabilities

media_type: text/csv

format: CSV

en:
This data set includes 589 snow profile observations including a rutschblock test, observations of signs of instability and an assessment of the local avalanche danger level, mainly recorded in the region of Davos (eastern Swiss Alps) during the winter seasons 2001-2002 to 2018-2019.

These data were analyzed and results published by Schweizer et al. (2021). They characterized the avalanche danger levels based on signs of instability (whumpfs, shooting cracks, recent avalanches), snow stability class and new snow height.

The data are provided in a csv file (589 records); the variables are described in the corresponding read-me file.

These data are the basis of the following publication: Schweizer, J., Mitterer, C., Reuter, B., and Techel, F.: Avalanche danger level characteristics from field observations of snow instability, Cryosphere, 15, 3293-3315, https://doi.org/10.5194/tc-15-3293-2021, 2021.

 Acknowlegements
Many of the data were recorded by SLF observers and staff members, among those Roland Meister, Stephan Harvey, Lukas Dürr, Marcia Phillips, Christine Pielmeier and Thomas Stucki. Their contribution is gratefully acknowledged.

avalanche

display_name: avalanche

avalanche-danger

display_name: avalanche-danger

avalanche-forecasting

display_name: avalanche-forecasting

danger-level

display_name: danger-level

snow

display_name: snow

envidat

en: EnviDat: the Environmental Data Portal


en: solar-irradiance agri-pv plane-of-array speed2zero rooftop-pv alpine-pv sweet-edge

en: This document provides an in-depth description of the CH-POA300 dataset. It includes methodology and validation details.

en:
Contains seasonal and annual statistics (mean, sum) of the CH-POA300 dataset.
- "yyyy_agg_period_hourly_stat.nc": statistic <stat> for the period <agg_period> of the year <yyyy>
- "all_agg_period_stat_mean.nc": statistic <stat> for the period <agg_period> reduced on the entire dataset (all years).
\
yyyy: [2009, 2019]
agg_period: [winter, summer, annual]
stat: [mean, sum]

en:
CH-POA300 is a set of maps of hourly solar irradiance projected onto oriented tilted surfaces, similar to monofacial solar panels, at a spatial resolution of 300m over Switzerland. Two temporal extents are available depending on the number of computed projections, i.e. azimuths and tilts combinations:
 1. Long term
 - Temporal extent: 2009-2019
 - Azimuths x Tilts: [180°] x [0°,35°,60°,90°]
 2 .Typical year
 - Temporal extent: 2016
 - Azimuths x Tilts: [0°,90°,135°,180°,225°,270°] x [0°,3°,8°, 15°, 25°, 35°, 45°, 60°, 80°, 90]
NB: Azimuth: 0°: N, 90°: E, 180°: S, 270°: W. Tilt: 0°: vertical, 90°: horizontal.
Pre-computed seasonal and annual statistics of these maps are also provided in the dataset.
Details on the method and the quality assessment can be found in the specifications (CH-POA300-v#-specs.pdf).

agri-pv

display_name: agri-pv

alpine-pv

display_name: alpine-pv

plane-of-array

display_name: plane-of-array

rooftop-pv

display_name: rooftop-pv

solar-irradiance

display_name: solar-irradiance

sweet-edge

display_name: sweet-edge

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: CH-POA300 -- Hourly Plane-of-Array Solar Irradiance at High Resolution for Switzerland


maintainer: Yves Bühler

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

en: Data Access

en:
- DSM (10 cm resolution)
- Ortho (3 cm resolution)
- HS (10 cm resolution)
- original RGB drone images
- Agisoft Metashape Processing Report

All products are also available as a cloud optimized geotiff.

en: To map and assess snow depth on different dates, 9 flights were conducted in the winter season of 2020/21 at the Latschüelfurgga in Davos. The data was captured with a Sony RX1R II mounted on a Wingtra drone and was processed with the Agisoft Metashape software. High-resolution DSMs, orthomosaics and snow height rasters, as well as the original RGB images from each flight are available.

drone-data

display_name: drone-data

dsm

display_name: dsm

orthophoto

display_name: orthophoto

photogrammetry

display_name: photogrammetry

snow-depth-mapping

display_name: snow-depth-mapping

uav

display_name: uav

wingtra

display_name: wingtra

en: Photogrammetric Drone Data Latschuelfurgga


maintainer: Mahdi Jafari

en: numerical-simulation convection-in-snow snow-density-profile openfoam vapor-transport

en: snowpackBuoyantPimpleFoam

en:
This zipped folder includes (1) the source code files, (2) an example case for simulation, and (3) an instruction how to compile and use the new solver.

Note that the most updated version of snowpackBuoyantPimpleFoam can be found in the WSL/SLF GitLab repository: https://gitlabext.wsl.ch/snow-models/snowpackbuoyantpimplefoam

en: snowpackBuoyantPimpleFoam is a two-phase solver implemented to model convection of water vapor with phase change in snowpacks. This new solver is based on the standard solver of buoyantPimpleFoam in the open-source fluid dynamics software, OpenFOAM 5.0 (www.openfoam.org).

convection-in-snow

display_name: convection-in-snow

numerical-simulation

display_name: numerical-simulation

openfoam

display_name: openfoam

snow-density-profile

display_name: snow-density-profile

vapor-transport

display_name: vapor-transport

educ

en: Education, culture and sport

en: snowpackBuoyantPimpleFoam: an OpenFOAM Eulerian–Eulerian two-phase solver for modelling convection of water vapor in snowpacks


lwf-tea-bag-sites

maintainer: Marcus Schaub

title_for_slug: lwf-tea-bag-sites

en: litter green-tea roiboos decomposition modeling verification tea-bag forest carbon

en: Site Information

en: Descriptive information on sites

mimetype: text/csv

en:
Decomposition of plant litter is a key process for the transfer of carbon and nutrients in ecosystems. Carbon contained in the decaying biomass is released to the atmosphere as respired CO2, and may contribute to global warming. Litterbag studies have been used to improve our knowledge of the drivers of litter decomposition, but they lack comparability because litter quality is plant species-specific. The use of commercial tea bags as a standard substrate was suggested in order to harmonize studies, where green tea and rooibos represent more labile and more recalcitrant C compounds as surrogates of local litter.
The tea bag approach was implemented on eight sites of the Swiss long-term Forest Ecosystem Research (LWF) network (https://www.wsl.ch/LWF). This allowed us to take advantage from the existing infrastructure and data from a previous litterbag study with local litter. In Beatenberg and Schaenis, additional elevation transects were established (1200-1800 m and 540-1150 m, respectively) to examine particularly the effect of temperature on decomposition. In Pfynwald (https://www.wsl.ch/de/ueber-die-wsl/versuchsanlagen-und-labors/flaechen-im-wald/pfynwald.html) and Salgesch, infrastructure of running projects was used to examine the effect of drought and understory removal, respectively. In Novaggio, tea bags were incubated in summer and winter to study the effect of seasonality particularly precipitation. Tea bags are collected after 3, 12, 24, and 36 months; for the two time-shifted experiments additionally after 6 and 9 months.
The study has two primary objectives. Firstly, it contributes to TeaComposition initiative (http://teacomposition.org/) which aims at investigating long-term litter decomposition and its key drivers at present as well as under different future climate scenarios using a common protocol and standard litter (tea) across nine terrestrial biomes. Secondly, the data are used to further develop decomposition models such as Yasso (http://en.ilmatieteenlaitos.fi/yasso) which is used by several countries, including Switzerland to estimate the annual carbon fluxes in dead wood, litter, and soil for reporting in National Greenhouse Gas Inventories under the United Nations Framework Convention on Climate Change and the Kyoto Protocol.

carbon

display_name: carbon

decomposition

display_name: decomposition

forest

display_name: forest

green-tea

display_name: green-tea

litter

display_name: litter

modeling

display_name: modeling

roiboos

display_name: roiboos

tea-bag

display_name: tea-bag

verification

display_name: verification

envidat

en: EnviDat: the Environmental Data Portal

en: LWF-Tea bag sites