dataset for "Climate-driven warming-herbicide synergies drive cryptic functional degradation of coastal oceans"

This repository contains reproducible R code and processed datasets accompanying manuscripts "Climate-driven warming-herbicide synergies drive cryptic functional degradation of coastal oceans" investigating climate-driven changes in coastal herbicide discharge and their cascading impacts on marine phytoplankton communities through warming-herbicide interactions. Repository contentsAnalysis Scripts1. herbicide_erosivity_model.R: XGBoost machine learning model projecting coastal herbicide discharge under RCP4.5 climate change scenarios2. herbicide_sst_phytoplankton_gam.R: Generalized Additive Models quantifying herbicide-SST interaction effects on nine phytoplankton functional groups3. phytoplankton_herbicide_meta_analysis.R: Meta-analysis of experimental phytoplankton herbicide sensitivity at environmentally relevant concentrations Processed datasets1. 50 km buffer spatially aggregated dataset for XGBoost analysis combining herbicide discharge, rainfall erosivity (ERA5-Land), application rates, and runoff coefficients herbicide_erosivity_50km_for_xgboost.csvherbicide_erosivity_50km_for_xgboost.rds2. 50 km buffer spatially aggregated dataset combining satellite-derived chlorophyll-a, chlorophyll-a related to phytoplankton functional types, herbicide discharge, sea surface temperature and nutrients for GAM analysis. herbicide_phyto_SST_integrated_50km_for_gam.csv herbicide_phyto_SST_integrated_50km_for_gam.rds Usage instructions# Load processed dataxgb_data <- readRDS("herbicide_erosivity_50km_for_xgboost.rds")# ORxgb_data <- read.csv("herbicide_erosivity_50km_for_xgboost.csv") # Run XGBoost analysissource("herbicide_erosivity_model.R") # Load GAM datagam_data <- readRDS("herbicide_phyto_SST_integrated_50km_for_gam.rds") # Run GAM analysissource("herbicide_sst_phytoplankton_gam.R") # Run meta-analysis (requires experimental data - included in supplementary files of the manuscript)source("phytoplankton_herbicide_meta_analysis.R")