A quasi-elastic aquifer deformational behavior: Madrid aquifer case study

The purpose of this paper is to analyze the quasi-elastic deformational behavior that has been induced by groundwater withdrawal of the Tertiary detrital aquifer of Madrid (Spain). The spatial and temporal evolution of ground surface displacement was estimated by processing two datasets of radar satellite images (SAR) using Persistent Scatterer Interferometry (PSI). The first SAR dataset was acquired between April 1992 and November 2000 by ERS-1 and ERS-2 satellites, and the second one by the ENVISAT satellite between August 2002 and September 2010. The spatial distribution of PSI measurements reveals that the magnitude of the displacement increases gradually towards the center of the well field area, where approximately 80 mm of maximum cumulated displacement is registered. The correlation analysis made between displacement and piezometric time series provides a correlation coefficient greater than 85% for all the wells. The elastic and inelastic components of measured displacements were separated, observing that the elastic component is, on average, more than 4 times the inelastic component for the studied period. Moreover, the hysteresis loops on the stress–strain plots indicate that the response is in the elastic range. These results demonstrate the quasi-elastic behavior of the aquifer. During the aquifer recovery phase ground surface uplift almost recovers from the subsidence experienced during the preceding extraction phase. Taking into account this unique aquifer system, a one dimensional elastic model was calibrated in the period 1997–2000. Subsequently, the model was used to predict the ground surface movements during the period 1992–2010. Modeled displacements were validated with PSI displacement measurements, exhibiting an error of 13% on average, related with the inelastic component of deformation occurring as a long-term trend in low permeability fine-grained units. This result further demonstrates the quasi-elastic deformational behavior of this unique aquifer system.This work was developed during Pablo Ezquerro research stay within the Geohazards InSAR laboratory and Modeling group of the Instituto Geológico y Minero de España in the framework of DORIS project (Ground Deformation Risk Scenarios: an Advanced Assessment Service) funded by the EC-GMES-FP7 initiative (Grant Agreement nº 242212). This work has been also supported by the Spanish Ministry of Science and Research (MICINN) under project TEC2011-28201-C02-02 and EU FEDER. Additional funding was obtained from Spanish Research Program through the project ESP2013-47780-C2-2-R

Analysis of regional large-gradient land subsidence in the Alto Guadalentín Basin (Spain) using open-access aerial LiDAR datasets

Land subsidence associated with groundwater overexploitation in the Alto Guadalentín Basin (Spain) aquifer system has been detected during the last decades. In this work, for the first time, we propose a new point cloud differencing methodology to detect land subsidence at basin scale, based on the multiscale model-to-model cloud comparison (M3C2) algorithm. This method is applied to two open-access airborne LiDAR datasets acquired in 2009 and 2016, respectively. First the internal edge connection errors in the different flight lines were addressed by means of a smoothing point cloud method. LiDAR datasets capture information from ground and non-ground points. Therefore, a method combining gradient filtering and cloth simulation filtering (CSF) algorithms was applied to remove non-ground points. The iterative closest point (ICP) algorithm was used for point cloud registration of both point clouds exhibiting a very stable and robust performance. The results show that vertical deformation rates are up to −14 cm/year in the basin from 2009 to 2016, in agreement with the displacement reported by previous studies. LiDAR results have been compared to the velocity measured by continuous GNSS stations and an InSAR dataset. For the GNSS-LiDAR and InSAR-LiDAR comparison, we computed a common 100 × 100 m grid in order to assess any similarities and discrepancies. The results show a good agreement between the vertical displacements obtained from the three different surveying techniques. Furthermore, LiDAR results were compared with the distribution of compressible soil thickness showing a clear relationship. The study underlines the potential of open-access and non-customized LiDAR to monitor the distribution and magnitude of vertical deformations in areas prone to be affected by groundwater-withdrawal-induced land subsidence.This research was funded by the ESA-MOST China DRAGON-5 project (ref. 59339) and by a Chinese Scholarship Council studentship awarded to Liuru Hu (Ref. 202004180062). María I. Navarro-Hernández and Guadalupe Bru are funded by the PRIMA programme supported by the European Union under grant agreement No 1924, project RESERVOIR

We lose ground: Global assessment of land subsidence impact extent

Depletion of groundwater aquifers along with all of the associated quality and quantity problems which affect profitability of direct agricultural and urban users and linked groundwater-ecosystems have been recognized globally. During recent years, attention has been devoted to land subsidence—the loss of land elevation that occurs in areas with certain geological characteristics associated with aquifer exploitation. Despite the large socioeconomic impacts of land subsidence most of these effects are still not well analyzed and not properly recognized and quantified globally. In this paper we developed a land subsidence impact extent (LSIE) index that is based on 10 land subsidence attributes, and applied it to 113 sites located around the world with reported land subsidence effects. We used statistical means to map physical, human, and policy variables to the regions affected by land subsidence and quantified their impact on the index. Our main findings suggest that LSIE increases between 0.1 and 6.5% by changes in natural processes, regulatory policy interventions, and groundwater usage, while holding all other variables unchanged. Effectiveness of regulatory policy interventions vary depending on the lithology of the aquifer system, in particular its stiffness. Our findings suggest also that developing countries are more prone to land subsidence due to lower performance of their existing water governance and institutions.Partial funding was provided by the Giannini Foundation of Agricultural Economics Minigrant Program. Dinar would like to acknowledge support from the W4190 Multistate NIFA-USDA-funded Project, “Management and Policy Challenges in a Water-Scarce World.” Esteban, Calvo, and Albiac would like to acknowledge support from the project INIA RTA2017-00082-00-00 by the Spanish Ministry of Science and Innovation, and support by funding to the research group ECONATURA from the Government of Aragon. Tomás would like to acknowledge support from the Spanish Ministry of Economy and Competitiveness, the State Agency of Research and the European Funds for Regional Development under project TEC2017-85244-C2-1-P. Tomás, Herrera, Ezquerro, and Teatini acknowledge the European Union support from the RESERVOIR project (GA nº 1924) developed in the framework of the PRIMA program

Evaluation of the potential of InSAR time series to study the spatio-temporal evolution of piezometric levels in the Madrid aquifer

The Tertiary detritic aquifer of Madrid (TDAM), with an average thickness of 1500 m and a heterogeneous, anisotropic structure, supplies water to Madrid, the most populated city of Spain (3.2 million inhabitants in the metropolitan area). Besides its complex structure, a previous work focused in the north-northwest of Madrid city showed that the aquifer behaves quasi elastically trough extraction/recovery cycles and ground uplifting during recovery periods compensates most of the ground subsidence measured during previous extraction periods (Ezquerro et al., 2014). Therefore, the relationship between ground deformation and groundwater level through time can be simulated using simple elastic models. In this work, we model the temporal evolution of the piezometric level in 19 wells of the TDAM in the period 1997–2010. Using InSAR and piezometric time series spanning the studied period, we first estimate the elastic storage coefficient (Ske) for every well. Both, the Ske of each well and the average Ske of all wells, are used to predict hydraulic heads at the different well locations during the study period and compared against the measured hydraulic heads, leading to very similar errors when using the Ske of each well and the average Ske of all wells: 14 and 16 % on average respectively. This result suggests that an average Ske can be used to estimate piezometric level variations in all the points where ground deformation has been measured by InSAR, thus allowing production of piezometric level maps for the different extraction/recovery cycles in the TDAM.This work is supported by the Spanish Ministry of Economy and Competitiveness and EU FEDER funds under projects TEC2011-28201-C02-02, TIN2014-55413-C2-2-P and ESP2013-47780-C2-2-R, by the Ministry of Education, Culture and Sport trough the project PRX14/00100 and by the project 15224/PI/10 from the Regional Agency of Science and Technology in Murcia

Vulnerability Assessment of Buildings due to Land Subsidence using InSAR Data in the Ancient Historical City of Pistoia (Italy)

The launch of the medium resolution Synthetic Aperture Radar (SAR) Sentinel-1 constellation in 2014 has allowed public and private organizations to introduce SAR interferometry (InSAR) products as a valuable option in their monitoring systems. The massive stacks of displacement data resulting from the processing of large C-B and radar images can be used to highlight temporal and spatial deformation anomalies, and their detailed analysis and postprocessing to generate operative products for final users. In this work, the wide-area mapping capability of Sentinel-1 was used in synergy with the COSMO-SkyMed high resolution SAR data to characterize ground subsidence affecting the urban fabric of the city of Pistoia (Tuscany Region, central Italy). Line of sight velocities were decomposed on vertical and E–W components, observing slight horizontal movements towards the center of the subsidence area. Vertical displacements and damage field surveys allowed for the calculation of the probability of damage depending on the displacement velocity by means of fragility curves. Finally, these data were translated to damage probability and potential loss maps. These products are useful for urban planning and geohazard management, focusing on the identification of the most hazardous areas on which to concentrate efforts and resources.This work was supported by the Spanish Ministry of Economy, Industry and Competitiveness (MINECO), the State Agency of Research (AEI) and European Funds for Regional Development (FEDER) under projects AQUARISK (ESP2013-47780-C2-2-R) and TEMUSA (TEC2017-85244-C2-1-P) and STAR-EO (TIN2014-55413-C2-2-P). The first author shows gratitude for the PhD student contract BES-2014-069076. The work was conceived during the research stay of P. Ezquerro and R. Tomás in the Università degli Studi di Firenze and the research stay of G. Herrera in the IGOT Lisbon University, supported by the Spanish Ministry of Education, Culture and Sport under fellowships EEBB-I-18-13014, PRX17/00439 and PRX19/00065, respectively. The S-1 monitoring activity is funded and supported by the Tuscany Region under the agreement “Monitoring ground deformation in the Tuscany Region with satellite radar data.” The authors also gratefully acknowledge TRE ALTAMIRA for having processed the S-1 data. The project was carried out using CSK® Products, © ASI (Italian Space Agency), delivered under the ASI Project Id Science 678 - “High resolution Subsidence investigation in the urban area of Pistoia (Tuscany Region, central Italy). The work is under the framework of the e-shape project, which has received funding from the European Union’s Horizon 2020 research and innovation program under grant agreement 820852. This paper is also supported by the PRIMA programme under grant agreement No 1924, project RESERVOIR. The PRIMA programme is supported by the European Union

Improving multi-technique monitoring using Sentinel-1 and Cosmo-SkyMed data and upgrading groundwater model capabilities

Aquifer-systems have become a strategic source of fresh water in the present climatic conditions, especially under stress in arid regions like the Iberian Mediterranean Arc. Understanding the behavior of groundwater reservoirs is crucial to their well-management and mitigation of adverse consequences of overexploitation. In this work, we use space geodetic measurements from satellite interferometric synthetic aperture radar (InSAR) and Global Positioning System (GPS) data, covering the period 2011–2017, to predict and validate the ground surface displacement over the fastest subsiding basin due to groundwater withdrawal in Europe (>10 cm/year). The 2D decomposition of InSAR displacements from Cosmo-SkyMed and Sentinel-1 satellites allows us to detect horizontal deformation towards the basin center, with a maximum displacement of 1.5 cm/year. InSAR results were introduced in a newly developed methodology for aquifer system management to estimate unknown pumping rates for the 2012–2017 period. This study illustrates how the combination of InSAR data, groundwater flow and deformation models can be used to improve the aquifer-systems sustainable management.This work was supported by the Spanish Ministry of Science, Innovation and Universities (MICINN), Spain; the State Agency of Research (AEI), Spain; and European Funds for Regional Development (FEDER), under projects AQUARISK (ESP2013-47780-C2-2-R), TEMUSA (TEC2017-85244-C2-1-P) and STAR-EO (TIN2014-55413-C2-2-P). First author shows gratitude for PhD student contract BES-2014-069076. A first version of this work was written during the research stay of first and second authors in the Università degli Studi di Firenze supported by the Spanish Ministry of Education, Culture and Sport, under fellowships EEBB-I-18-13014 and PRX17/00439, respectively

Imaging land subsidence in the Guadalentín River Basin (SE Spain) using Advanced Differential SAR Interferometry

Aquifer overexploitation can lead to the irreversible loss of groundwater storage caused by the compaction or consolidation of unconsolidated fine-grained sediments resulting in land subsidence. Advanced Differential SAR Interferometry (A-DINSAR) is particularly efficient to monitor progressive ground movements, making it an appropriate method to study depleting aquifers undergoing overexploitation and land subsidence. The Guadalentín River Basin (Murcia, Spain) is a widely recognized subsiding area that exhibits the highest rates of groundwater-related land subsidence recorded in Europe (>10 cm/yr). The basin covers an extension of more than 500 km2 and is underlain by an overexploited aquifer-system formed by two contiguous hydraulically connected units (Alto Guadalentín and Bajo Guadalentín). Although during the last years the piezometric levels have partially stabilized, the ongoing aquifer-system deformation is evident and significant, as revealed by the A-DInSAR analysis presented. In this work, we submit the first vertical and horizontal (E-W) decomposition results of the LOS velocity and displacement time series of the whole Guadalentín Basin obtained from two datasets of Sentinel-1 SAR acquisitions in ascending and descending modes. The images cover the period from 2015 to 2021 and they were processed using the Parallel Small BAseline Subset (P-SBAS) implemented by CNRIREA in the Geohazards Exploitation Platform (GEP) on-demand web tool, which is funded by the European Space Agency. The output ascending and descending measurement points of P-SBAS lie on the same regular grid, which is particularly suited for the geometrical decomposition. Time series displacements are compared to a permanent GNSS station located in the Bajo Guadalentín basin.This study has received funding in framework of the RESERVOIR project (Sustainable groundwater RESources managEment by integrating eaRth observation deriVed monitoring and flOw modelIng Results), funded by the Partnership for Research and Innovation in the Mediterranean Area (PRIMA) programme supported by the European Union (Grant Agreement 1924; https://reservoir-prima.org/). The study has also been supported by the Grant FPU19/03929 (funded by MCIN/AEI/10.13039/501100011033 and by “FSE invests in your future”); the Project CGL2017-83931-C3-3-P (funded by MCIN/ AEI/10.13039/501100011033 and by “ERDF A way of making Europe”); the ESA-MOST China DRAGON-5 Project (ref. 59339) and the SARAI Project PID2020-116540RB-C22 (funded by MCIN/AEI/10.13039/501100011033). Copernicus Sentinel-1 IW SAR data were provided and processed in ESA’s Geohazards Exploitation Platform (GEP), in the framework of the GEP Early Adopters Programme

3D groundwater flow and deformation modelling of Madrid aquifer

A novel methodological approach to calibrate and validate three-dimensional (3D) finite element (FE) groundwater flow and geomechanical models has been implemented using Advanced Differential Interferometric SAR (A-DInSAR) data. In particular, we show how A-DInSAR data can be effectively used to (1) constrain the model set-up in evaluating the areal influence of the wellfield and (2) characterise the aquifer system, specifically the storage coefficient values, which represents a fundamental step in managing groundwater resources. The procedure has been tested to reconstruct the surface vertical and horizontal movements caused by the Manzanares-Jarama wellfield located northwest of Madrid (Spain). The wellfield was used to supply freshwater during major droughts over the period between 1994 and 2010. Previous A-DInSAR outcomes obtained by ERS-1/2 and ENVISAT acquisitions clearly revealed the seasonality of the land displacements associated to the withdrawal and recovery cycles that characterized the wellfield development. A time-lag of about one month, which is in the order of the time span between two SAR acquisitions, between the hydraulic head changes and the displacements has been detected in this site by a wavelet analysis of A-DInSAR and piezometer time series. The negligible delay between the forcing factor and the system response and the complete subsidence recovery when piezometric head recovers supported the understanding of a minor role played by the pore pressure propagation within clay layers and the almost perfectly elastic behavior of the system (viscosity is negligible), respectively. The developed geomechanical model satisfactorily reproduces the pumping-induced deformations with a Root Mean Square Error (RMSE) between observed and simulated land displacements in the order of 0.1-0.3 mm. The results give insights about the approach benefits in deeply understanding the spatio-temporal aquifer-system response to the management of this strategic water resource for Madrid.The activity has been developed within the scientific collaboration established in the framework of the UNESCO Land Subsidence International Initiative (LaSII - https://www.landsubsidence-unesco.org/). The research was funded by University of Pavia in the framework of a research grant award “assegno di tipo A premiale” for research activities at the Dept. of Earth and Environmental Sciences, within the research project entitled “Sustainable groundwater resources management by integrating A-DInSAR derived monitoring and flow modeling results” assigned to Roberta Bonì in March 2019. This research was partially funded by the Spanish Ministry of Economy, Industry and Competitiveness (MINECO), the State Agency of Research (AEI), and the European Funds for Regional Development (FEDER) under project TEC2017-85244-C2-1-P

Understanding the dynamic behaviour for the Madrid aquifer (Spain): insights from the integration of A-DInSAR and 3-D groundwater flow and geomechanical models

Advanced Differential Interferometric Synthetic Aperture Radar (A-DInSAR) techniques and 3-D groundwater flow and geomechanical models are integrated to improve our knowledge about the Tertiary detritic aquifer of Madrid (TDAM). In particular, the attention is focused on the Manzanares-Jarama well field, located to the northwest of Madrid, which experienced five cycles of extensive groundwater withdrawal followed by natural recovery, to cope with the droughts occurred in summer 1995, 1999, 2002, 2006, and 2009. Piezometric records and A-DInSAR data acquired by ERS-1/2 and ENVISAT satellites during the periods 1992–2000 and 2002–2010, respectively, have been used to calibrate the groundwater flow and the geomechanical models. A time-lag of about one month between the hydraulic head changes and the displacements of the land surface has been detected by a joint wavelet analysis of A-DInSAR and piezometer head time series. Overall, the results show the effectiveness of the proposed integrated approach composed of A-DInSAR and 3-D geomechanical model to characterize the aquifer-system response during and after the groundwater withdrawal.This research has been supported by the Università degli Studi di Pavia (Assegno di tipo A premiale: “Sustainable groundwater resources management by integrating A-DInSAR derived monitoring and flow modeling results” assigned to Roberta Bonì in March 2019) and the Spanish Ministry of Economy, Industry and Competitiveness (MINECO), the State Agency of Research (AEI), and the European Funds for Regional Development (FEDER) (grant no. TEC2017-85244-C2-1-P)

Wavelet analysis of land subsidence time-series: Madrid Tertiary aquifer case study

Interpretation of land subsidence time-series to understand the evolution of the phenomenon and the existing relationships between triggers and measured displacements is a great challenge. Continuous wavelet transform (CWT) is a powerful signal processing method mainly suitable for the analysis of individual nonstationary time-series. CWT expands time-series into the time-frequency space allowing identification of localized nonstationary periodicities. Complementarily, Cross Wavelet Transform (XWT) and Wavelet Coherence (WTC) methods allow the comparison of two time-series that may be expected to be related in order to identify regions in the time-frequency domain that exhibit large common cross-power and wavelet coherence, respectively, and therefore are evocative of causality. In this work we use CWT, XWT and WTC to analyze piezometric and InSAR (interferometric synthetic aperture radar) time-series from the Tertiary aquifer of Madrid (Spain) to illustrate their capabilities for interpreting land subsidence and piezometric time-series information.This research has been supported by the Spanish Ministry of Economy and Competitiveness, the State Agency of Research and the European Funds for Regional Development (grant no. TEC2017-85244-C2-1-P)