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

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.