Future abundance and distribution of key bird species for pathogen transmission in the Netherlands
Data files
Oct 25, 2024 version files 57.32 MB
-
climate_ranges.csv
437 B
-
modelInputs.zip
52.27 MB
-
outputs.zip
5.04 MB
-
README.md
5.52 KB
Abstract
We have modelled the current and future abundance and distribution of three different bird species in the Netherlands: blackbirds, mallards, and house sparrows. These are all species which are widespread in north-western Europe and which are likely to play a role in the transmission of several pathogens. We used random forest models created using the SDMaps R package. Climate, land use, and vegetative cover were all found to be important predictors. We found that the abundance of mallards and house sparrows is expected to increase in the future, while the abundance of blackbirds shows little change. These changes will likely have consequences for disease risk. By providing these future abundance maps, we are enabling detailed modelling of the future risk of disease outbreaks.
https://doi.org/10.5061/dryad.r2280gbmc
Description of the data and file structure
The following data are included. All raster files are on a 1km grid with CRS EPSG: 28992.
- modelInputs.zip
- predictors
- climate: All climatic predictor rasters used for the modelling, divided by scenario:
- H: high emissions, used for SSPs 3 and 5
- M: medium emissions, used for SSP4
- L: low emissions, used for SSP1
- r: reference scenario
- Predictors have ‘breeding’ and ‘winter’ variants. Breeding means the breeding season and is different for each species (see the accompanying paper), winter is December 1st to February 28th (or 29th in leap years) and is the same for all species. The predictors are:
- cold_days: number of days with minimum temperature below 0C
- prec_covar: precipitation coefficient of variation
- prec_dry: precipitation (mm) in driest 4-week period
- prec_tot: total precipitation (mm)
- temp_max: average maximum daily temperature (C)
- temp_mean: mean daily temperature (C)
- temp_min: average minimum daily temperature (C)
- temp_range: average temperature range over the time period (C)
- ref: All non-climate predictors rasters used for the reference scenario. The numbers in the file names indicate the buffer size
- The predictors are:
- altitude (m)
- arable: proportion of arable farmland
- CBC: Common Bird Census dataset (kept as zero)
- clay: clay percentage in soil
- forest: proportion of forested area
- France: French dataset (kept as zero)
- grass: proportion of grass area
- industry_commerce_transport: proportion of industry, commerce or transport area
- MAS: Meetnet Agrarische Soorten dataset (kept as zero)
- mine_dump_construction: proportion of mine, dump or construction area
- MUS: Meetnet Urbane Soorten (kept as zero)
- open_unvegetated: proportion of open invegetated area
- other_rural: proportion of rural area which does not fall under other rural types
- pasture: proportion of pasture area
- permanent crops: proportion of permanent cropland
- sand: sand percentage in soil
- tree: tree density
- urban: proportion of urban area
- urban_green: proportion of urban green area
- vegetation: proportion of vegetated area
- water_inland: proportion of inland water area
- water_wetland_marine: proportion of marine water or wetland area
- wetland_inland: proportion of inland wetland area
- The predictors are:
- SSP1/3/4/5: All non-climate predictors rasters used for each SSP scenario. The numbers in the file names indicate the buffer size
- The predictors are the same as for the reference scenario
- climate: All climatic predictor rasters used for the modelling, divided by scenario:
- predictors
- Outputs
- predictions: raster files with predicted relative abundance for each species and scenario
- uncertainty: raster files with the confidence interval for the predictions for each species for the reference period
- variableImportance: csv files with the variable importance of each predictor for each species
- anovaResults.csv: results of pairwise comparisons between predictions to test for significant differences. Results shown are afer the application of the Bonferroni multiple comparison test. The numbers after the species indicate the SSP scenario and ‘R’ indicates the reference prediction
- climate_ranges.csv: Climate ranges used in determining areas of France where the current climate is similar to the future Dutch climate. These values were taken from summarised statistics from future climate scenarios produced by the Dutch Meteorological Association: https://cdn.knmi.nl/system/data_center_publications/files/000/071/901/original/KNMI23_klimaatscenarios_gebruikersrapport_23-03.pdf?1696579253
Code
All code was written in R, using R version 4.0.4 (the SDMaps package does not run in newer versions). The model generation was performed using the SDMaps package v0.15-6.
- makePredictors.R: Generates the current and future predictors (as raster files). This includes the derivation of future land use, as described in the document Deriving future and use.pdf
- runScenarios.R: Model generation and making reference and future predictions
- uncertainty.R: Runs bootstrapped model and calculates confidence interval of predictions
- anova.R: Performs anova on model predictions to check for significant differences
Supplementary material
- Data summary.pdf: A summary of the bird abundance data used for the modelling. The full dataset is available on request from the Dutch Centre for Field Ornithology (https://www.sovon.nl/)) and the Common Bird Monitoring Scheme database for France (https://www.vigienature.fr/fr/suivi-temporel-des-oiseaux-communs-stoc).
- Deriving future land use.pdf: Descibes the land use classifications used and the procedure for deriving future land use, grass and tree cover maps.
- Partial dependence plots.pdf: Shows the predictor maps and partial dependence plots for each species for all model covariates.
We used bird abundance data from the Netherlands and France to train random forest models, using the SDMaps package in R (Sierdsema et al., n.d.). Future scenario data was taken from future climate scenarios produced by the KNMI (Dutch Meteorological Organisation) ((KNMI, 2023) and the Dutch One Health SSPs (Dellar et al., n.d., 2024).
The reference or 'current' time period used was 1991-2020. The future period was 2036-2065.
Full details can be found in the accompanying article.
Dellar, M., Geerling, G., Kok, K., van Bodegom, P., Schrama, S., & Boelee, E. (2024). Creating the Dutch One Health Shared Socio-economic Pathways (SSPs). Regional Environmental Change, 24, 16. https://doi.org/https://doi.org/10.1007/s10113-023-02169-1
Dellar, M., Geerling, G., Kok, K., van Bodegom, P., van der Schrier, G., Schrama, M., & Boelee, E. (n.d.). Future land use maps for the Netherlands based on the Dutch One Health Shared Socio-economic Pathways. Submitted.
KNMI. (2023). KNMI National Climate Scenarios 2023 for the Netherlands. knmi.nl/klimaatscenarios
Sierdsema, H., Kampichler, C., & Hallmann, C. (n.d.). Creating distribution maps from monitoring data and casual observations. Retrieved March 21, 2022, from https://www.sovon.nl/nl/content/trimmaps