High quality InSAR data linked to seasonal change in hydraulic head for an agricultural area in the San Luis Valley, Colorado

In the San Luis Valley (SLV), Colorado legislation passed in 2004 requires that hydraulic head levels in the confined aquifer system stay within the range experienced in the years 1978–2000. While some measurements of hydraulic head exist, greater spatial and temporal sampling would be very valuable in understanding the behavior of the system. Interferometric synthetic aperture radar (InSAR) data provide fine spatial resolution measurements of Earth surface deformation, which can be related to hydraulic head change in the confined aquifer system. However, change in cm-scale crop structure with time leads to signal decorrelation, resulting in low quality data. Here we apply small baseline subset (SBAS) analysis to InSAR data collected from 1992 to 2001. We are able to show high levels of correlation, denoting high quality data, in areas between the center pivot irrigation circles, where the lack of water results in little surface vegetation. At three well locations we see a seasonal variation in the InSAR data that mimics the hydraulic head data. We use measured values of the elastic skeletal storage coefficient to estimate hydraulic head from the InSAR data. In general the magnitude of estimated and measured head agree to within the calculated error. However, the errors are unacceptably large due to both errors in the InSAR data and uncertainty in the measured value of the elastic skeletal storage coefficient. We conclude that InSAR is capturing the seasonal head variation, but that further research is required to obtain accurate hydraulic head estimates from the InSAR deformation measurements

Modeling the two- and three-dimensional displacement field in Lorca, Spain, subsidence and the global implications

Land subsidence associated with overexploitation of aquifers is a hazard that commonly affects large areas worldwide. The Lorca area, located in southeast Spain, has undergone one of the highest subsidence rates in Europe as a direct consequence of long-term aquifer exploitation. Previous studies carried out on the region assumed that the ground deformation retrieved from satellite radar interferometry corresponds only to vertical displacement. Here we report, for the first time, the two- and three-dimensional displacement field over the study area using synthetic aperture radar (SAR) data from Sentinel-1A images and Global Navigation Satellite System (GNSS) observations. By modeling this displacement, we provide new insights on the spatial and temporal evolution of the subsidence processes and on the main governing mechanisms. Additionally, we also demonstrate the importance of knowing both the vertical and horizontal components of the displacement to properly characterize similar hazards. Based on these results, we propose some general guidelines for the sustainable management and monitoring of land subsidence related to anthropogenic activities.Peer ReviewedPostprint (published version

Remote Sensing and Modeling of Stressed Aquifer Systems and the Associated Hazards

abstract: Aquifers host the largest accessible freshwater resource in the world. However, groundwater reserves are declining in many places. Often coincident with drought, high extraction rates and inadequate replenishment result in groundwater overdraft and permanent land subsidence. Land subsidence is the cause of aquifer storage capacity reduction, altered topographic gradients which can exacerbate floods, and differential displacement that can lead to earth fissures and infrastructure damage. Improving understanding of the sources and mechanisms driving aquifer deformation is important for resource management planning and hazard mitigation. Poroelastic theory describes the coupling of differential stress, strain, and pore pressure, which are modulated by material properties. To model these relationships, displacement time series are estimated via satellite interferometry and hydraulic head levels from observation wells provide an in-situ dataset. In combination, the deconstruction and isolation of selected time-frequency components allow for estimating aquifer parameters, including the elastic and inelastic storage coefficients, compaction time constants, and vertical hydraulic conductivity. Together these parameters describe the storage response of an aquifer system to changes in hydraulic head and surface elevation. Understanding aquifer parameters is useful for the ongoing management of groundwater resources. Case studies in Phoenix and Tucson, Arizona, focus on land subsidence from groundwater withdrawal as well as distinct responses to artificial recharge efforts. In Christchurch, New Zealand, possible changes to aquifer properties due to earthquakes are investigated. In Houston, Texas, flood severity during Hurricane Harvey is linked to subsidence, which modifies base flood elevations and topographic gradients.Dissertation/ThesisDoctoral Dissertation Geological Sciences 201

Prediction of sustainable management and associated land subsidence features in the Toluca aquifer system

EbookRegional land subsidence accompanying groundwater abstraction in the Toluca aquifer-system is a challenge for managing groundwater resources and mitigating associated hazards. In order to improve this situation, groundwater management scenarios for the Toluca Valley are examined with a three dimensional groundwater flow model coupled to a one dimensional compaction module. Subsequently, the land subsidence evolution was investigated by integrating SAR interferometry and geological and hydrogeological data to shed insight on the underlying processes governing subsidence. The results indicate that continuing at current rates of water consumption will lead to subsidence of more than 1.6 m over a 40 year period (2010–2050). Completely stopping exports to Mexico City is not the most important factor in controlling subsidence because the pumping system is mostly located in regions with low clay content, where subsidence is lower. However, decreasing exports by half and relocating the pumping centres to low-clay-content areas does have a positive effect on the overall water budget and subsidence. From 2003 to 2016, groundwater level declines of up to 1.6 m/yr, land subsidence up to 77 mm/yr, and major infrastructure damages are observed. Groundwater level data show highly variable seasonal responses according to their connectivity to recharge areas. However, the trend of groundwater levels consistently range from −0.5 to −1.5 m/yr regardless of the well location and depth. By analysing the horizontal gradients of vertical land subsidence, we provide a potential ground fracture map to assist in future urban development planning in the Toluca Valley. The approach taken in this study could be applied to their locations with similar problems in order to determine the most viable option for water supply.CONACyT UAEM Institute for Groundwater Research, University of Guelp

Characterization of Ground Deformation Associated with Shallow Groundwater Processes Using Satellite Radar Interferometry

Shallow groundwater processes maylead to ground deformation and even geohazards. With the features of day-and-night accessibility and large-scale coverage, time-series interferometric synthetic aperture radar (InSAR) has proven a useful tool for mapping the deformation over various landscapes at cm to mm level with weekly to monthly updates. However, it has limitations such as, decorrelation,atmospheric artifacts, topographic errors, andunwrapping errors, in particular for the hilly, vegetated, and complicated deformation patterns. In this dissertation, I focus on characterizing the ground deformation over landslides, aquifer systems, and mine tailings impoundment, using the designed advanced time-series InSAR strategy, as well as theinterdisciplinary knowledge of geodesy, hydrology, geophysics, and geology. Northwestern USA has been exposed to extreme landslide hazards due to steep terrain, high precipitation, and loose root support after wildfire. I characterize the rainfall-triggered movements of Crescent Lake landslide, Washington State. The seasonal deformation at the lobe, with larger magnitudes than the downslope riverbank, suggests an amplified hydrological loading effect due to a thicker unconsolidated zone. High-temporal-resolution InSAR and GPS data reveal dynamic landslide motions. Threshold rainfall intensities and durations wet seasons have been associated with observed movement upon shearing: antecedent rainfall triggered precursory slope-normal subsidence, and the consequent increase in pore pressure at the basal surface reduces friction and instigates downslope slip over the course of less than one month. In addition, a quasi-three-dimensional deformation field is created using multiple spaceborne InSAR observations constrained by the topographical slope, and is further used to invert for the complex geometry of landslide basal surface based on mass conservation. Aquifer skeletons deform in response to hydraulic head changes with various time scales of delay and sensitivity. I investigate the spatio-temporal correlation among deformation, hydrological records and earthquake records over Salt Lake Valley, Utah State. A clear long-term and seasonal correlation exists between surface uplift/subsidence and groundwater recharge/discharge, allowing me to quantify hydrogeological properties. Long-term uplift reflects the net pore pressure increase associated with prolonged water recharge, probably decades ago. The distributions of previously and newly mapped faults suggest that the faultsdisrupt the groundwater flow andpartition hydrological units. Mine tailings gradual settle as the pore pressure dissipates and the terrain subsides, andtailings embankment failures can be extremely hazardous. I investigate the dynamics of consolidation settlement over the tailings impoundment in the vicinity of Great Salt Lake, Utah State, as well as its associated impacts to the surrounding infrastructures. Largest subsidence has been observed around the low-permeable decant pond clay at the northeast corner.The geotechnical consolidation model reveals and predicts the long-term exponentially decaying settlement process. My studies have demonstrated that InSAR methods can advance our understanding about the potential anthropogenic impacts and natural hydrological modulations on various geodynamic settings in geodetic time scale

Reservoir characterization in an underground gas storage field using joint inversion of flow and geodetic data

Characterization of reservoir properties like porosity and permeability in reservoir models typically relies on history matching of production data, well pressure data, and possibly other fluid-dynamical data. Calibrated (history-matched) reservoir models are then used for forecasting production and designing effective strategies for improved oil and gas recovery. Here, we perform assimilation of both flow and deformation data for joint inversion of reservoir properties. Given the coupled nature of subsurface flow and deformation processes, joint inversion requires efficient simulation tools of coupled reservoir flow and mechanical deformation. We apply our coupled simulation tool to a real underground gas storage field in Italy. We simulate the initial gas production period and several decades of seasonal natural gas storage and production. We perform a probabilistic estimation of rock properties by joint inversion of ground deformation data from geodetic measurements and fluid flow data from wells. Using an efficient implementation of the ensemble smoother as the estimator and our coupled multiphase flow and geomechanics simulator as the forward model, we show that incorporating deformation data leads to a significant reduction of uncertainty in the prior distributions of rock properties such as porosity, permeability, and pore compressibility.Eni S.p.A. (Firm

Urban Deformation Monitoring using Persistent Scatterer Interferometry and SAR tomography

This book focuses on remote sensing for urban deformation monitoring. In particular, it highlights how deformation monitoring in urban areas can be carried out using Persistent Scatterer Interferometry (PSI) and Synthetic Aperture Radar (SAR) Tomography (TomoSAR). Several contributions show the capabilities of Interferometric SAR (InSAR) and PSI techniques for urban deformation monitoring. Some of them show the advantages of TomoSAR in un-mixing multiple scatterers for urban mapping and monitoring. This book is dedicated to the technical and scientific community interested in urban applications. It is useful for choosing the appropriate technique and gaining an assessment of the expected performance. The book will also be useful to researchers, as it provides information on the state-of-the-art and new trends in this fiel

Development and Application of a 1D Compaction Model to Understand 65 Years of Subsidence in the San Joaquin Valley

High rates of land subsidence, caused by groundwater overdraft, are resulting in millions of dollars of infrastructure damage in California\u27s San Joaquin Valley (SJV). In recent years, the use of interferometric synthetic aperture radar (InSAR) has enabled us to substantially improve our understanding of this subsidence. However, only very occasionally have the InSAR data been integrated with a physical model of subsurface compaction. Here, we have used InSAR and other data to parameterize and calibrate a 1D compaction model. We applied our model to a study area in the SJV where we had access to the necessary information on hydraulic head to develop model inputs. Our model simulated subsidence in the three aquifer system layers over the period 1952–2017, and is the first 1D compaction model in the SJV to simulate multiple aquifer system layers from the 1950s to 2017. The results from our model suggest that previous studies have significantly underestimated the time constants governing the slow, residual compaction of subsurface clays. We suggest that residual compaction of clays is a process that continues for decades-to-centuries, indicating that to significantly reduce subsidence requires some recovery of head, not just a stabilization. We also show how compaction in the lower, confined aquifer has accounted for over 90% of subsidence in the past 20 years. Although our study area is small, our findings are likely representative of the subsiding regions of the SJV, and our methodology can be applied to unconsolidated aquifer systems exhibiting subsidence worldwide

Assessment of groundwater exploitation and land subsidence development in the Toluca aquifer system, Mexico

EbookRegional land subsidence accompanying groundwater abstraction in the Toluca aquifer-system is a challenge for managing groundwater resources and mitigating associated hazards. In order to improve this situation, groundwater management scenarios for the Toluca Valley are examined with a three dimensional groundwater flow model coupled to a one dimensional compaction module. Subsequently, the land subsidence evolution was investigated by integrating SAR interferometry and geological and hydrogeological data to shed insight on the underlying processes governing subsidence. The results indicate that continuing at current rates of water consumption will lead to subsidence of more than 1.6 m over a 40 year period (2010–2050). Completely stopping exports to Mexico City is not the most important factor in controlling subsidence because the pumping system is mostly located in regions with low clay content, where subsidence is lower. However, decreasing exports by half and relocating the pumping centres to low-clay-content areas does have a positive effect on the overall water budget and subsidence. From 2003 to 2016, groundwater level declines of up to 1.6 m/yr, land subsidence up to 77 mm/yr, and major infrastructure damages are observed. Groundwater level data show highly variable seasonal responses according to their connectivity to recharge areas. However, the trend of groundwater levels consistently range from −0.5 to −1.5 m/yr regardless of the well location and depth. By analysing the horizontal gradients of vertical land subsidence, we provide a potential ground fracture map to assist in future urban development planning in the Toluca Valley. The approach taken in this study could be applied to their locations with similar problems in order to determine the most viable option for water supply.CONACyT, Institute for Groundwater Research, University of Guelph UAE

Remote estimation of geologic composition using interferometric synthetic-aperture radar in California's Central Valley

California's Central Valley is the national agricultural center, producing 1/4 of the nation's food. However, land in the Central Valley is sinking at a rapid rate (as much as 20 cm per year) due to continued groundwater pumping. Land subsidence has a significant impact on infrastructure resilience and groundwater sustainability. In this study, we aim to identify specific regions with different temporal dynamics of land displacement and find relationships with underlying geological composition. Then, we aim to remotely estimate geologic composition using interferometric synthetic aperture radar (InSAR)-based land deformation temporal changes using machine learning techniques. We identified regions with different temporal characteristics of land displacement in that some areas (e.g., Helm) with coarser grain geologic compositions exhibited potentially reversible land deformation (elastic land compaction). We found a significant correlation between InSAR-based land deformation and geologic composition using random forest and deep neural network regression models. We also achieved significant accuracy with 1/4 sparse sampling to reduce any spatial correlations among data, suggesting that the model has the potential to be generalized to other regions for indirect estimation of geologic composition. Our results indicate that geologic composition can be estimated using InSAR-based land deformation data. In-situ measurements of geologic composition can be expensive and time consuming and may be impractical in some areas. The generalizability of the model sheds light on high spatial resolution geologic composition estimation utilizing existing measurements.Comment: 10 pages, 7 figures, NeurIPS 202