Impact model: BioScen1.5-SDM-GBM

Sector
Terrestrial biodiversity
Region
global

The BioScen1.5-SDM-GBM models are species distribution models of 15496 terrestrial vertebrate species (2964 amphibians, 8493 birds and 4039 mammals ) using Generalized Boosted Models as model algorithm. BioScen1.5-SDM-GBM are statistical models, which were calibrated using a 30-year average (1980-2009) of the EWEMBI observed climate input data set and IUCN and BirdLife range maps of each vertebrate species. The models are then projected to simulated historical, current and future 30-year average climate data from ISIMIP2b, resulting in probabilities of occurrence and summed probabilities of occurrence for the three different vertebrate taxa.

Information for the model BioScen1.5-SDM-GBM is provided for the simulation rounds shown in the tabs below. Click on the appropriate tab to get the information for the simulation round you are interested in.

Person responsible for model simulations in this simulation round
Matthias Biber: matthias.biber@tum.de, 0000-0002-7726-988X, Technical University of Munich (Germany)
Christian Hof: christian.hof@tum.de, 0000-0002-7763-1885, Technical University of Munich (Germany)
Output Data
Experiments: II, III
Climate Drivers: None
Date: 2018-11-23
Basic information
Model Version: BioScen15-SDM-GBM-v1
Model Output License: CC BY 4.0
Reference Paper: Main Reference: Hof C, Voskamp A, Biber M, Böhning-Gaese K, Engelhardt E, Niamir A, Willis S, Hickler T et al. Bioenergy cropland expansion may offset positive effects of climate change mitigation for global vertebrate diversity. Proceedings of the National Academy of Sciences,115,13294-13299,2018
Reference Paper: Other References:
Resolution
Spatial aggregation: regular grid
Horizontal resolution: 0.5°x0.5°
Input data
Simulated atmospheric climate data sets used: IPSL-CM5A-LR, HadGEM2-ES, GFDL-ESM2M, MIROC5
Observed atmospheric climate data sets used: EWEMBI
Other data sets used: Land-sea mask
Additional input data sets: Present biodiversity data for amphibians, terrestrial mammals and terrestrial birds were obtained from expert range maps provided by the International Union for Conservation of Nature (2011) and BirdLife International (2014) and gridded to a resolution of 0.5°x0.5°.
Climate variables: tasmax, tasmin, pr
Additional information about input variables: 30-year monthly means of minimum temperature (tasmin), maximum temperature (tasmax) and total precipitation (pr) were calculated and used to derive 19 bioclimatic variables (see Hijmans et al. 2005) . The best combination of 4 variables was selected for each taxon (amphibians, birds, mammals) by comparing the performance of all uncorrelated combinations of 4 variables. As a result, Bio4 (temperature seasonality), Bio5 (max temperature of warmest month), Bio12 (annual precipitation) and Bio15 (precipitation seasonality), were used for modelling birds and mammals and Bio4, Bio5, Bio18 (precipitation of warmest quarter) and Bio19 (precipitation of coldest quarter) were used for modelling amphibians.
Model specifications
Model algorithm: GBM
Explanatory variables: Bio4 (temperature seasonality), Bio5 (max temperature of warmest month), Bio12 (annual precipitation) and Bio15 (precipitation seasonality) were used for modelling birds and mammals. Bio4, Bio5, Bio18 (precipitation of warmest quarter) and Bio19 (precipitation of coldest quarter) were used for modelling amphibians.
Response variable: absence/presence of species
Additional information about response variable: Bio4 (temperature seasonality), Bio5 (max temperature of warmest month), Bio12 (annual precipitation) and Bio15 (precipitation seasonality) were used for modelling birds and mammals. Bio4, Bio5, Bio18 (precipitation of warmest quarter) and Bio19 (precipitation of coldest quarter) were used for modelling amphibians.
Parameters: n.trees=10500, n.minobsinnode = 10, interaction.depth=c(1,2,3), shrinkage=c(0.01, 0.001), bag.fraction=0.5
Additional information about parameters: Model performance was evaluated using cross-validation and the model parameters (combination of interaction.depth and shrinkage) that resulted in the lowest cross-validation error were used for the final model. Only models with an AUC of > 0.7 were kept as results.
Software function: gbm()
Software package: gbm
Software program: R
Additional information about model output: The model output is the probability of occurrence of a species, which can vary between 0 and 1 (see for example amphibian-prob output files). Note: Probability of occurrence is projected to the present and all neighbouring realms of a species and so sort of represents the unlimited dispersal of a species in the future. The output (probability of occurrence) of all species is then stacked (added up) to a summed probability of occurrence (a proxy of species richness) for (i) all species, (ii) all endemic species and (iii) all threatened species of the three modelled vertebrate taxa (see for example amphibian-sum-prob output files). Summed probability of occurrence is split into different dispersal scenarios (no dispersal, d/2, d, 2*d, full dispersal). Full dispersal represents the sum of the probability of occurrence output files. No dispersal assumes that species can only be present where they are actually present according to the IUCN and BirdLife range maps. The other three dispersal scenarios consider species-specific dispersal buffers added to the present range, where d is the largest diameter of the original range of the species.
Note of caution
Note of caution about applicability of model use due to modeling concept: Projecting models to pre-industrial and extended future time periods was only performed to allow for comparability with the other modeling sectors and should be used with caution. Historical projections are limited to the species that were modeled and do not include species which went extinct over the last decades or centuries (which we did not model). Historical richness is thus very unlikely to be a realistic estimate of past richness as it is probably highly underestimated. Future and especially extended future projections should also be used with caution, as species interactions and other processes are not considered when running individual models even though they do influence species distributions and may change in the future, altering the future distribution of a species. Furthermore, the distribution data used to calibrate the models are range maps drawn by experts based on current and historical distribution records as well as additional information and expert knowledge. Thus, they do not represent the distribution at a specific point in time, but a rough outline of a species current distribution. Therefore, precise projections of individual years are not sensible for this modelling approach.
Person responsible for model simulations in this simulation round
Matthias Biber: matthias.biber@tum.de, 0000-0002-7726-988X, Technical University of Munich (Germany)
Christian Hof: christian.hof@tum.de, 0000-0002-7763-1885, Technical University of Munich (Germany)
Output Data
Experiments: historical
Climate Drivers: None
Date: 2018-11-23
Basic information
Model Version: BioScen15-SDM-GBM-v1
Model Output License: CC BY 4.0
Reference Paper: Main Reference: Hof C, Voskamp A, Biber M, Böhning-Gaese K, Engelhardt E, Niamir A, Willis S, Hickler T et al. Bioenergy cropland expansion may offset positive effects of climate change mitigation for global vertebrate diversity. Proceedings of the National Academy of Sciences,115,13294-13299,2018
Reference Paper: Other References:
Resolution
Spatial aggregation: regular grid
Horizontal resolution: 0.5°x0.5°
Input data
Observed atmospheric climate data sets used: EWEMBI
Other data sets used: Land-sea mask
Climate variables: tasmax, tasmin, pr
Model specifications
Model algorithm: GBM
Explanatory variables: Bio4 (temperature seasonality), Bio5 (max temperature of warmest month), Bio12 (annual precipitation) and Bio15 (precipitation seasonality) were used for modelling birds and mammals. Bio4, Bio5, Bio18 (precipitation of warmest quarter) and Bio19 (precipitation of coldest quarter) were used for modelling amphibians.
Response variable: absence/presence of species
Additional information about response variable: Range maps only give us information about the presence of a species, absence data was created by using a distance kernel and random selection of at least 1000 points.
Distribution of response variable: Binomial
Parameters: n.trees=10500, n.minobsinnode = 10, interaction.depth=c(1,2,3), shrinkage=c(0.01, 0.001), bag.fraction=0.5
Additional information about parameters: Model performance was evaluated using cross-validation and the model parameters (combination of interaction.depth and shrinkage) that resulted in the lowest cross-validation error were used for the final model. Only models with an AUC of > 0.7 were kept as results.
Software function: gbm()
Software package: gbm
Software program: R
Additional information about model output: The model output is the probability of occurrence of a species, which can vary between 0 and 1 (see for example amphibian-prob output files). Note: Probability of occurrence is projected to the present and all neighbouring realms of a species and so sort of represents the unlimited dispersal of a species in the future. The output (probability of occurrence) of all species is then stacked (added up) to a summed probability of occurrence (a proxy of species richness) for (i) all species, (ii) all endemic species and (iii) all threatened species of the three modelled vertebrate taxa (see for example amphibian-sum-prob output files). Summed probability of occurrence is split into different dispersal scenarios (no dispersal, d/2, d, 2*d, full dispersal). Full dispersal represents the sum of the probability of occurrence output files. No dispersal assumes that species can only be present where they are actually present according to the IUCN and BirdLife range maps. The other three dispersal scenarios consider species-specific dispersal buffers added to the present range, where d is the largest diameter of the original range of the species.