Data required to run the analyses in 'Bumble bee probability of occurrence responds to interactions between local and landscape land use, climatic niche properties and climate change'

Contained here are all of the data needed to run the analyses in 'Bumblebee occupancy responds to complex interactions between local and landscape land use, climatic niche properties and climate change', as detailed in the analysis code on GitHub (https://github.com/timnewbold/BumblebeesLUClimStudy).
The datasets contained in this repository are:
BUMBLEBEE NICHE PROPERTIES
Estimates of the realized climatic niches of bumblebees, presented as the temperature (BaselineTEI_Spp) or precipitation niche position (BaselinePEI_Spp) across each species distribution, where a value of 0 represents the coldest or driest conditions and 1 the warmest or wettest conditions within the species realized niche. Also estimates of how these temperature and precipitation position values have changed under recent climate change (DeltaTEI_Spp_Period3 and DeltaPEI_Spp_Period3).
The creation of these estimates is described in Soroye et al. (2020), DOI: 10.1126/science.aax8591

LAND-USE HISTORY
Data on the year in which landscapes were first substantially converted to human land uses, where substantial use is defined either as 10% (yocPS15002005_1_0.1), 30% (yocPS15002005_1_0.3) or 50% (yocPS15002005_1_0.5) conversion to human land-use types (cropland, pasture or urban areas).
For more information on these estimates, see Newbold et al. (2015), DOI: 10.1038/nature14324

UN SUB-REGIONS MAP
A map of UN sub-regions (UN_subregion), for plotting purposes

Datasets that are required to run the analysis but omitted here are those for which we do not have permission to republish, but which are nevertheless available in other public repositories. Instructions for obtaining these additional datasets are as follows:
COVERAGE OF NATURAL HABITAT
Data on the coverage of natural (primary and secondary) habitats can be downloaded here:
https://data.csiro.au/collection/csiro:15276v3
These data then need to be pre-processed by multiplying by 1000 and converting to an integer value, with the resulting raster files saved as:
'pri_1km_int' and 'sec_1km_int'

FERTILIZERS
Global data on the density of application (kg/ha) of nitrogen, phosphorous and potassium fertilizer on 17 major crops can be downloaded here:
http://www.earthstat.org/nutrient-application-major-crops/
The downloaded data need to be in a folder 'FertilizerCropSpecific_Geotiff'
PESTICIDES
Global data on the density of application (kg/ha) of different types of pesticides on different agricultural crops can be downloaded (in GeoTIFF format) here:
https://sedac.ciesin.columbia.edu/data/set/ferman-v1-pest-chemgrids-v1-01/data-download
The downloaded data need to be in the following directory structure:
'PEST-CHEMGRIDS_v1_01_APR/GEOTIFF'

PREDICTS DATABASE
The PREDICTS database, documenting the structure of ecological assemblages across different land-use types and intensities, can be downloaded from the data portal of the Natural History Museum:
http://dx.doi.org/10.5519/0066354
You should download the RDS-format file, and save as database.rds

ELEVATION
A Digital Elevation Model at 30 arc-second spatial resolution can be downloaded from WorldClim:
https://www.worldclim.org/data/worldclim21.html
Filename should be wc2.1_30s_elev.tif