Environmental Data: Amphibian and urban-rural landscape data Swiss Low...

Amphibian monitoring data

Processed amphibian occurrence data (2017-2019) in the cantons of Aargau and Zürich. The dataset includes information for the following species: _Alytes obstetricans_,_Bombina variegata_, _Bufo bufo_, _Hyla arborea_, _Ichtyosaura alpestris_, _Lissotriton helveticus_, _Lissotriton vulgaris_, _Rana temporaria_, _Salamandra salamandra_, and _Triturus cristatus_.

Amphibian breeding sites of national importance - cantons of Aargau and Zürich

Amphibian breeding sites of national importance (focal nodes) in the cantons of Aargau and Zürich, 30x30 m resolution, binarised ascii format.

Amphibian whole life cycle predictors - 10m resolution

Amphibian whole life cycle predictors include 6 variable types: (1) Topographic – slope [°]: derived from the Swiss Topographical Landscape model (swissTLM 3D 2015) at 2 m resolution (elevation uncertainty ± 0.5m) and, subsequently, aggregated to a median value; (2) Hydrologic: two predictors including nearest distance to water [m] (all water-related objects and rasterized lotic, lentic and reedbeds from the regional cadastral maps – AV: Daten der Amtliche Vermessung) and runoff coefficient [C] (between 0 and 1 relating amount of runoff to amount of precipitation, mapped AV: Daten der Amtliche Vermessung data against literature values for different land covers); (3) Edaphic – soil moisture variability (not available): can gathered from the spatial modelling of Ecological Indicator Values in Switzerland (Descombes et al., 2020) and resampled from 100m to 10m; (4) Vegetation: we calculated the median and standard deviation of the Normalized Difference Vegetation Index (NDVI) from all available cloud free Sentinel-2 pixels available between April and October for the years 2016 to 2019 (Drusch et al., 2012) and the median vegetation height [m] (not available) can gathered from Ginzler and Hobi (2015); (5) Land-use derived: includes three predictors, grassland density [grassland/cell] (not available) - can be gathered from Pazur et al. (2021), total grasslands area evaluated using a spatial weight matrix and a 250m neighborhood radius around each cell; urbanization [buildings/cell] – average building density for a 250m neighborhood around each cell based on the AV: Daten der Amtliche Vermessung ; and traffic intensity [vehicles/day] – average daily vehicle count (Nationales Personenverkehrsmodells); (6) Movement-ecology related: includes three predictors of proximity – nearest forest distance [m]; nearest rock-gravel and sandy area (RGS) distance [m] and nearest road distance [m] (considering any road type) – all calculated based on the AV: Daten der Amtliche Vermessung maps.

Binarized "green" and "grey" landcover data (Cantons of Aargau and Zürich)

"Green" landcover refers to lawns and vegetated urban surfaces such as, for example private gardens, public parks, cemeteries, green stripes along roadsides, etc. While "grey" landcover refers to impervious grounds such as roads, troittoirs, other impervious. Binarized maps, 30x30m resolution.

Amphibian species specific landscape resistance maps.

For each amphibian species we used the weighted average Ensemble model suitability map and converted it to a resistance map using an exponential decay function following (Duflot et al., 2018) using the species-specific mpa threshold previously computed (Preau et al., 2020). The function assigns a resistance value of 1,000 when habitat suitability equals zero, and a value of one when habitat suitability is greater than or equal to the binarization threshold. The exponential shape strengthens the barrier effect of less favorable areas.