Dimensionality reduction for groundwater forecasting under drought and intensive irrigation with neural networks

Abstract

Study Region: This study focuses on the Berrechid aquifer system in northern Morocco. Study Focus: The research explores Principal Component Analysis (PCA) for optimizing input selection in groundwater level forecasting using neural networks. PCA efficiently reduces input dimensionality while preserving critical information, making it beneficial for neural network modelling of natural systems with extensive input variables in a low-resource scenarios requiring feature engineering. A Long Short-Term Memory (LSTM) model predicted groundwater levels in six monitoring bores using four hydro-climatic variables, precipitation, land surface temperature (LST), actual evapotranspiration (AET), and the normalized difference vegetation index (NDVI). Model performance was compared using two approaches: the LSTM-XGB model with the best-selected input features and the LSTM-PC1 model based on the first principal component (PC1). New Hydrological Insights for the Region: Results showed that NDVI, AET, and LST were the dominant inputs across different monitoring bores. On average, PC1 accounted for 68.3 % of the variance in hydro-climatic variables, with an eigenvalue of 2.75, surpassing the combined variance of two individual hydro-climatic variables. Both models performed effectively, achieving R² values of 0.982–0.999 during training and 0.885–0.999 during validation. The models successfully captured groundwater fluctuations and the declining trend during drought. LSTM-XGB slightly outperformed LSTM-PC1 in certain cases, but the differences were minimal. The use of PC1 not only mitigates overfitting risks but also allows for generalized predictions across multiple monitoring sites, making it a practical choice for large datasets. • PC1, with an average eigenvalue of 2.75, explains more variance than two individual hydro-climatic variables. • PC1 captures on average 68.3 % of the total variance in hydro-climatic variables. • The LSTM-PC1 model exhibited exceptional performance, achieving R² values as high as 0.999. • PCA reduces overfitting risks and enhances model generalization across different monitoring wells. • PCA-driven modelling provides an efficient and scalable approach for groundwater forecasting in data-limited environments.