R code and datasets to support "Stability and changes of global climate diversity by 2100"

Abstract

Climate Diversity Calculation — Code and Data This repository accompanies the manuscript "Stability and changes of global climate diversity by 2100". It provides a reproducible example of the climate diversity (CD) calculation on a focal study region, using the actual 15 standardised bioclimatic variables and seven moving-window sizes. --- 1. System requirements Software dependencies - R (version 4.4.2)- RStudio (2026.08.0, recommended)- R packages (install once): ```rinstall.packages(c("terra"))``` Operating system - Windows 10 (tested). Tested versions - R 4.4.2 on Windows 10 with RStudio 2026.08.0. Non-standard hardware None. A standard desktop or laptop computer is sufficient. 2. Installation guide No installation is required beyond R and the packages above. 1. Download and unzip the archive.2. Open the project in RStudio.3. Set the working directory to the unzipped folder.4. Install missing R packages if prompted: ```rinstall.packages(c("terra"))``` 3. Repository structure upload-zenodo/├── cd-calculate_CUR_multi_mw.R # Main CD calculation script├── cur-std.R # mask input variables│ ├── cur_bio_1.tif ... cur_bio_17.tif├── cur-result/ # CD results for seven window sizes (GeoTIFF)│ ├── disc_CUR_mw3.tif│ ├── disc_CUR_mw5.tif│ ├── disc_CUR_mw7.tif│ ├── disc_CUR_mw9.tif│ ├── disc_CUR_mw11.tif│ └── disc_CUR_mw13.tif└── study_area/ └── study_area.tif # Study area Folder descriptions - `cur-std/` — The 15 standardised bioclimatic variables (bio1–bio7, bio10–bio17; bio8/9/18/19 are excluded), already clipped to the study area and aligned to a common 5 km grid in the EPSG:6933 equal-area projection.- `cur-result/` — Pre-computed climate diversity GeoTIFFs for the current period under seven moving-window sizes (mw = 3, 5, 7, 9, 11, 13 pixels).- `study_area/` — The study region (`study_area.tif`), aligned pixel-by-pixel with the variables in `cur-std/`.- `cd-calculate_CUR_multi_mw.R` — The main script that reads the 15 standardised variables, computes CD for each window size, and writes outputs to `cur-result/`. 4. Demo We provide a ready-to-run example on a focal study region across seven moving-window sizes . Full recomputation 1. Open `cd-calculate_CUR_multi_mw.R`.2. Update the file paths and `window_sizes` vector at the top of the script.3. Run the script. It reads the 15 variables in `cur-std/`, computes CD for each window size, and writes `disc_CUR_mw*.tif` to `cur-result/`. Expected output - For each window size, a GeoTIFF `disc_CUR_mw*.tif` in `cur-result/`.- Values are non-negative; cells outside the study area are `NA`. Expected run time - Full recomputation for the focal region (seven window sizes): approximately 5–10 minutes on a standard Windows desktop (Intel i7, 32 GB RAM). 5. Applying the code to your own data Prepare input layers. Standardise each variable (z-score) before computing CD — the equal-weight Euclidean distance in the formula requires comparable units across variables. Update paths and window sizes. In `cd-calculate_CUR_multi_mw.R`, point the input directory to your standardised variables and set `window_sizes` as needed (e.g. `c(3, 5, 7, 9, 11, 13)`). Run the script. Source the script.