Endurance-Recovery-Resilience (ERR) Framework: Code and Data for Dynamic Groundwater Resilience Assessment

ERR Framework: Endurance-Recovery-Resilience for Groundwater Systems ===================================================================== Authors : J. Akhil, I. Bhatnagar, C. T. Dhanya Affiliation: Department of Civil Engineering, IIT Delhi WHAT THIS DOES -------------- Takes seasonal rainfall and groundwater level (GWL) time series as input and computes, per the manuscript's Methodology (Eq. 1-11): - Rainfall-induced groundwater fluctuations, h_rf(t) (Eq. 1-3) - Groundwater Response to Rainfall, GRR(t) (Eq. 4-5) - Response regimes and statistically validated transitions (Sec 4.3.1) - Mean Response Level (MRL) per response regime (Eq. 6) - Endurance (E) and Recovery (R), raw and normalized (Eq. 7-11) - Resilience class assignment (Critical to Dynamic, 3x3 ERR grid) FILES ----- err_framework.py Main script (run this) README.txt This file requirements.txt Python dependencies Data/gwl_ts.nc Seasonal GWL input (place here before running) Data/rain_ts.nc Seasonal rainfall input (place here before running) HOW TO RUN ---------- Step 1 - Install dependencies pip install -r requirements.txt Step 2 - Place input files in the Data/ folder Data/gwl_ts.nc NetCDF file Dimensions : basin_code x time Variable : Level (groundwater level, m bgl) Time : seasonal (4 seasons/year: Dec-Feb, Mar-May, Jun-Aug, Sep-Nov) Data/rain_ts.nc NetCDF file Dimensions : basin_code x time Variable : Rainfall (seasonal rainfall, mm) Step 3 - Run python err_framework.py Step 4 - Check outputs in Data/ ERR_Table_S3_Full_Diagnostics.csv MRL, CV_GRR, CV_Rain, Damping, E_raw/E, N_fail/N_trans/P, R_raw/R, Resilience_Class per basin and response regime (reproduces the manuscript's Table S3) ERR_Summary_Wide.csv One row per basin (both response regimes side by side) PARAMETERS (adjustable at the top of err_framework.py) ------------------------------------------------------- WAVELET Wavelet family default: db4 DECOMP_LEVEL Decomposition level default: 6 MAX_LAG Max lag (seasons) searched in the component-rainfall weighting (Eq. 2) default: 3 MIN_RR_STEPS Min response-regime duration default: 40 (= 10 years x 4 seasons) EXCLUDED Basin codes to skip default: BAR, BRM, IND DATA SOURCES ------------ Groundwater levels : Central Ground Water Board (CGWB) https://nwic.in/wris Rainfall : India Meteorological Department (IMD) Pai et al. (2014), Journal of Earth System Science https://imdpune.gov.in Note: Raw data are not included. Basin-aggregated NetCDF files used in this study are provided in the Zenodo archive.