Habitat contractions to drive future planktonic foraminifera distribution and diversity change

Description This repository contains MATLAB code and data supporting the manuscript currently in review titled "Habitat contractions to drive future planktonic foraminifera population and diversity loss". The scripts and datasets included allow reproduction of abundance and diversity maps under various climate change scenarios (SSPs). This work includes the usage of two ecological niche models models (ESM and GRNN) to study the abundance distribution of 6 planktonic foraminifera species and biodiversity projections across past and future scenarios. Files and Descriptions Scripts esm_diversity_allspecies_script.m and grnn_diversity_allspecies_script.m Scripts to create abundance maps for all 41 species under all 3 SSP scenarios using a common set of model parameters derived from the 6 studied species. diversity_stackspeciesmaps_script.mCode to convert abundance maps into presence/absence maps with variable thresholds and generate species richness maps for each of 5 ecological groups. Requires prior execution of the above abundance map scripts. esm_6species_script.m and grnn_6species_script.mScripts to run model simulations individually for each of the 6 species and each SSP scenario for the ESM and GRNN models, respectively. Contains toggleable sections for running specific species-scenario combinations. Datasets esm_6species_allmaps.mat and grnn_6species_allmaps.matPast, future, and anomaly abundance maps for the 6 studied species for all SSP scenarios. esm_5ecogroups_allmaps.mat and grnn_5ecogroups_allmaps.matDiversity maps (past, future, anomaly) for each SSP scenario, categorized into 5 ecological groups: all species, common species, symbiont-bearing species, symbiont-barren species, and facultatively symbiont species. env_inputs_gridded.matEnvironmental variable datasets (chlorophyll, shortwave downward radiation, thermocline depth, temperature) gridded at 5°x5° spatial resolution and binned temporally by 5-year intervals, for past and future SSP scenarios. species_inputs.mat Ungridded species abundance data for the 6 studied species, split into pseudo-random training and validation datasets with columns: longitude, latitude, mean sample depth, year (binned), and abundance. diversity_inputs.matUngridded original abundance data divided by ecological groups, with columns: longitude, latitude, mean sample depth, year, and abundance. esm_41species_allmaps.mat and grnn_41species_allmaps.mat Abundance maps for all 41 species used for diversity mapping. Generated from esm_diversity_allspecies_script and grnn_diversity_allspecies_script using a common model parameter set. Keywords planktonic foraminifera, species distribution, climate change, SSP scenarios, ecological modeling, ENM, GRNN, ESM, biodiversity