Understanding Long-Term Groundwater Storage Variability Using GRACE Data and Explainable Machine Learning

Abstract

The decline in groundwater storage (GWS) poses a critical threat to water security in semi-arid regions where increasing agricultural water demand and climate variability are increasing pressure on aquifers. This study presents a novel hybrid modeling framework integrating multi-source satellite and climate data (GRACE, GLDAS, TerraClimate, and MODIS) with machine learning and explanatory artificial intelligence techniques for the long-term assessment and interpretation of GWS anomalies in the data-poor Iğdır Basin. Three different modeling approaches were developed: XGBoost, Long Short-Term Memory (LSTM) networks, and their combined model, and interpreted using the Shapley Additive Explanations (SHAP) method. The results showed a significant long-term decreasing trend in groundwater storage anomalies at a rate of −0.87 mm per month during the 2002–2016 period, indicating continuous depletion. The LSTM model demonstrated the best performance with R2 of 0.59, RMSE of 19.5 mm, and MAE of 15.1 mm, revealing the dominant role of temporal dependencies in groundwater systems. SHAP analysis identified lagged groundwater anomalies (especially GWS_lag3) as the most effective predictors; this may reflect the memory effect and lagged response specific to semi-arid aquifer systems, but this interpretation needs to be validated in different study areas. Snow water equivalent and total water storage anomalies also emerged as significant determinants, while the direct effect of instantaneous precipitation was found to be limited. This study addresses significant gaps in the literature by combining sequence-based modeling with model interpretability in a semi-arid closed basin. The findings highlight the necessity of using system memory and explainable artificial intelligence together for reliable groundwater prediction. While the proposed hybrid approach has the potential for application in other semi-arid regions, its broader usability needs to be supported by independent validation studies under different hydrogeological and climatic conditions.