A twenty-year period of severe land subsidence evolution in the Alto Guadalentín Basin (southeast Spain) is monitored using multi-sensor SAR images, processed by advanced differential interferometric synthetic aperture radar (DInSAR) techniques. The SAR images used in this study consist of four datasets acquired by ERS-1/2, ENVISAT, ALOS and COSMO-SkyMed satellites between 1992 and 2012. The integration of ground surface displacement maps retrieved for different time periods allows us to quantify up to 2.50 m of cumulated displacements that occurred between 1992 and 2012 in the Alto Guadalentín Basin. DInSAR results were locally compared with global positioning system (GPS) data available for two continuous stations located in the study area, demonstrating the high consistency of local vertical motion measurements between the two different surveying techniques. An average absolute error of 4.6 ± 4 mm for the ALOS data and of 4.8 ± 3.5 mm for the COSMO-SkyMed data confirmed the reliability of the analysis. The spatial analysis of DInSAR ground surface displacement reveals a direct correlation with the thickness of the compressible alluvial deposits. Detected ground subsidence in the past 20 years is most likely a consequence of a 100–200 m groundwater level drop that has persisted since the 1970s due to the overexploitation of the Alto Guadalentín aquifer system. The negative gradient of the pore pressure is responsible for the extremely slow consolidation of a very thick (> 100 m) layer of fine-grained silt and clay layers with low vertical hydraulic permeability (approximately 50 mm/h) wherein the maximum settlement has still not been reached.This work was developed during Roberta Bonì's stay within the Geohazards InSAR Laboratory and Modeling Group of the Instituto Geológico y Minero de España in the framework of the DORIS project (Ground Deformation Risk Scenarios: an Advanced Assessment Service) funded by the EC-GMES-FP7 initiative (Grant Agreement 423 no. 242212). ALOS PALSAR images were provided by the project JAXA-1209. Part of this work is supported by the Spanish Government under project TEC2011-28201-C02 and by the project 15224/PI/10 from the Regional Agency of Science and Technology in Murcia. Additional funding was obtained from the Spanish Research Program through the projects AYA2010-17448, ESP2013-47780-C2-1-R and ESP2013-47780-C2-2-R and by the Ministry of Education, Culture and Sport through the project PRX14/00100

There is an increasing need for accurate groundwater-level (GWL) prediction to support effective seasonal water management. It is desirable for forecasting tools to be not only accurate but also accessible for decision-makers. We test the Prophet forecasting procedure, an open-source code released by Facebook, to address these challenges. It is based on an additive model considering non-periodic changes and periodic components in a Bayesian framework with easily-interpretable parameters. Predictions of daily GWL data in an area affected by pumping near a tourist complex in the Ramsar wetland area of Doñana (Spain) are compared to other forecasting methods. Prophet outperforms most methods in predicting GWL making it a fast and flexible forecasting tool for hydrologists and water managers. Furthermore, it allows gaining insight into the influence of each component of the forecast separately, helping to assess the hydrodynamic response to external drivers such as groundwater pumping.

Land subsidence resulting from groundwater extractions is a global phenomenon adversely affecting many regions worldwide. Understanding the governing processes and mitigating associated hazards require knowing the spatial distribution of the implicated factors (piezometric levels, lithology, ground deformation), usually only known at discrete locations. Here, we propose a methodology based on the Kriging with External Drift (KED) approach to interpolate sparse point measurements of variables influencing land subsidence using high density InSAR measurements. In our study, located in the Alto Guadalentín basin, SE Spain, these variables are GPS vertical velocities and the thickness of compressible soils. First, we estimate InSAR and GPS rates of subsidence covering the periods 2003–2010 and 2004–2013, respectively. Then, we apply the KED method to the discrete variables. The resulting continuous GPS velocity map shows maximum subsidence rates of 13 cm/year in the center of the basin, in agreement with previous studies. The compressible deposits thickness map is significantly improved. We also test the coherence of Sentinel-1 data in the study region and evaluate the applicability of this methodology with the new satellite, which will improve the monitoring of aquifer-related subsidence and the mapping of variables governing this phenomenon.

Groundwater resources are regularly the principal water supply in semiarid and arid climate areas. However, groundwater levels (GWL) in semiarid aquifers are suffering a general decrease because of anthropic exploitation of aquifers and the repercussions of climate change. Effective groundwater management strategies require a deep characterization of GWL fluctuations, in order to identify individual behaviors and triggering factors. In September 2019, the Guadalquivir River Basin Authority (CHG) declared that there was over-exploitation in three of the five groundwater bodies of the Almonte-Marismas aquifer, Southwest Spain. For that reason, it is critical to understand GWL dynamics in this aquifer before the new Spanish Water Resources Management Plans (2021–2027) are developed. The application of GWL series clustering in hydrogeology has grown over the past few years, as it is an extraordinary tool that promptly provides a GWL classification; each group can be related to different responses of a complex aquifer under any external change. In this work, GWL time series from 160 piezometers were analyzed for the period 1975 to 2016 and, after data pre-processing, 24 piezometers were selected for clustering with k-means (static) and time series (dynamic) clustering techniques. Six and seven groups (k) were chosen to apply k-means. Six characterized types of hydrodynamic behaviors were obtained with time series clustering (TSC). Number of clusters were related to diverse affections of water exploitation depending on soil uses and hydrogeological spatial distribution parameters. TSC enabled us to distinguish local areas with high hydrodynamic disturbance and to highlight a quantitative drop of GWL during the studied period.