Use of Geophysical and Remote Sensing Data for Assessment of Aquifer Depletion and Related Land Deformation

Intellligent sustainable agricultural water practice using multi sensor spatiotemporal evolution

The amount of water taken from non-renewable resources such as aquifers to fulfill irrigation requirements is rarely monitored, putting sustainable agriculture under threat in the face of changing climate. In the present research, an attempt was made to apply multi-sensor (Landsat suite, GRACE, GRACE-FO) satellite data to monitor spatiotemporal evolution of agriculture for the Al-Qassim region, Kingdom of Saudi Arabia (KSA). For this purpose, time series of NDVI (Normalized Difference Vegetation Index), SAVI (Soil-Adjusted Vegetation Index), and MSAVI2 (Modified Soil-Adjusted Vegetation 2) was utilized to assess vegetation pattern change in the study area. The present investigation used High-resolution Planetscope (PS) nanosatellite data to validate the vegetation results. Mann Kendall trend analysis and linear regression were performed to study the temporal pattern, and the relationship between vegetation, GRACE, and climate variables was performed from 1984 to 2020. Water extraction based on the averaged value of JPL GWS and CSR GWS showed a decreasing trend of -10.24 ± 1.4 mm/year from 2003-2020. The annual rainfall showed a decreasing trend, while the annual temperature showed an increasing trend from 1982-2020. The correlation of vegetation indices with rainfall of one-month lag showed a significantly better relationship of 0.74, 0.74, and 0.75, respectively, for NDVI, SAVI, and MSAVI2. The correlation between temperature and all three vegetation indices is a strong negative correlation: -0.85 for NDVI and -0.9 for SAVI and MSAVI.

Rapid groundwater decline and some cases of recovery in aquifers globally

Groundwater resources are vital to ecosystems and livelihoods. Excessive groundwater withdrawals can cause groundwater levels to decline1-10, resulting in seawater intrusion11, land subsidence12,13, streamflow depletion14-16 and wells running dry17. However, the global pace and prevalence of local groundwater declines are poorly constrained, because in situ groundwater levels have not been synthesized at the global scale. Here we analyse in situ groundwater-level trends for 170,000 monitoring wells and 1,693 aquifer systems in countries that encompass approximately 75% of global groundwater withdrawals18. We show that rapid groundwater-level declines (>0.5 m year-1) are widespread in the twenty-first century, especially in dry regions with extensive croplands. Critically, we also show that groundwater-level declines have accelerated over the past four decades in 30% of the world's regional aquifers. This widespread acceleration in groundwater-level deepening highlights an urgent need for more effective measures to address groundwater depletion. Our analysis also reveals specific cases in which depletion trends have reversed following policy changes, managed aquifer recharge and surface-water diversions, demonstrating the potential for depleted aquifer systems to recover.

Remote Sensing Application for Landslide Detection, Monitoring along Eastern Lake Michigan (Miami Park, MI)

We assessed the nature and spatial and temporal patterns of deformation over the Miami Park bluffs on the eastern margin of Lake Michigan and investigated the factors controlling its observed deformation. Our approach involved the following steps: (1) extracting bluff deformation rates (velocities along the line of sight of the satellite) using a stack of Sentinel-1A radar imagery in ascending acquisition geometry acquired between 2017 and 2021 and applying the Intermittent Small Baseline Subset (ISBAS) InSAR time series analysis method; (2) generating high-resolution (5 cm) elevation models and orthophotos from temporal unmanned aerial vehicle (UAV) surveys acquired in 2017, 2019, and 2021; and (3) comparing the temporal variations in mass wasting events to other relevant datasets including the ISBAS-based bluff deformation time series, lake level (LL) variations, and local glacial stratigraphy. We identified areas witnessing high line-of-sight (LOS) deformation rates (up to −21 mm/year) along the bluff from the ISBAS analysis and seasonal deformation patterns associated with freeze-thaw cycles, suggesting a causal effect. The acceleration of slope failures detected from field and UAV acquisitions correlated with high LLs and intensified onshore wave energy in 2020. The adopted methodology successfully predicts landslides caused by freezes and thaws of the slope face by identifying prolonged slow deformation preceding slope failures, but it does not predict the catastrophic landslides preceded by short-lived LOS deformation related to LL rise.

Using InSAR Time Series to Monitor Surface Fractures and Fissures in the Al-Yutamah Valley, Western Arabia

Western Arabia routinely experiences geophysical phenomena that deform the surface of the earth in a variety of ways. These phenomena include earthquakes, volcanic eruptions, sinkholes, and earth fissuring and fracturing. We perform a time-series analysis of interferometric synthetic aperture radar (InSAR) observations derived from the ESA Sentinel-1 radar satellite constellation to map regional surface displacements in western Arabia as a function of time. We rely on InSAR products generated by the JPL-Caltech ARIA project to detect regions with short wavelength anomalies, and then manually reprocess InSAR products at a higher resolution for these regions to maximize spatial and temporal coverage. We post-process InSAR products using MintPy workflows to develop the InSAR time series. We report short wavelength anomalies localized within alluvial valleys across western Arabia and find a 5 cm/year line-of-sight surface displacement within the Al-Yutamah Valley. Part of the observed subsidence is correlated with surface fractures that developed in conjunction with severe rainfall events in regions characterized mainly by alluvial sediments at the surface. Regions of observed subsidence that are not associated with any surface fractures or fissures are correlated with the presence of basalt layers at the surface. Both regions are subject to groundwater exploitation. The observed subsidence is inferred to be driven by groundwater withdrawal perhaps modulated by the presence of a preexisting depositional environment (e.g., paleo-lake deposits) that promotes unconsolidated soil compaction.

Land Subsidence in the Texas Coastal Bend: Locations, Rates, Triggers, and Consequences

Land subsidence and sea level rise are well-known, ongoing problems that are negatively impacting the entire Texas coast. Although ground-based monitoring techniques using long-term global navigation satellite systems (GNSS) records provide accurate subsidence rates, they are labor intensive, expensive, time-consuming, and spatially limited. In this study, interferometric synthetic aperture radar (InSAR) data and techniques were used to map the locations and quantify rates of land subsidence in the Texas Coastal Bend region during the period from October 2016 to July 2019. InSAR-derived land subsidence rates were then validated and calibrated against GNSS-derived rates. The factors controlling the observed land subsidence rates and locations were investigated. The consequences of spatial variability in land subsidence rates in Coastal Bend were also examined. The results indicated that: (1) land subsidence rates in the Texas Coastal Bend exhibited spatial variability, (2) InSAR-derived land subsidence rates were consistent with GNSS-derived deformation rates, (3) land subsidence in the Texas Coastal Bend could be attributed mainly to hydrocarbon and groundwater extraction as well as vertical movements along growth faults, and (4) land subsidence increased both flood frequency and severity in the Texas Coastal Bend. Our results provide valuable information regarding not only land deformation rates in the Texas Coastal Bend region, but also the effectiveness of interferometric techniques for other coastal rural areas around the globe.

Land Subsidence Induced by Rapid Urbanization in Arid Environments: A Remote Sensing-Based Investigation

The rapid increase in the population of many of the older major cities within the countries of the Saharan-Arabian Desert is steering vast and disorganized urban expansion and in many cases introducing adverse environmental impacts such as soil erosion, rise in groundwater levels, and contamination of shallow aquifers, as well as development of deformational features including land subsidence. Using the rapidly growing city of Riyadh (1992: 467 km2; 2018: 980 km2), the capital of the Kingdom of Saudi Arabia as a test site, we utilized Small Baseline Subset (SBAS) interferometric analyses of 2016 to 2018 Sentinel-1 images together with multi-temporal high-resolution images viewable on Google Earth, GPS, field, land use land cover (LULC), and geological data to assess the distribution and rates of land subsidence and their causal effects. Three main causes of subsidence were identified and assessed: (1) discharge of wastewater effluents from septic systems in newly urbanized areas that lead to an increase in soil moisture, rise in groundwater levels, waterlogging, and wetting and hydrocompaction of dry alluvium loose sediments causing land subsidence (up to −20 mm/y) in wadis and lowlands; (2) the subsurface dissolution of karst formation by wastewater effluents and the collapse of voids and cavities at depth under stresses introduced by heavy construction machinery, causing sagging and land subsidence (up to −5 mm/y); and (3) leveling, compaction, and degradation of municipal and building waste materials in organized landfills and disorganized dump sites that resulted in significant land subsidence (up to −21 mm/y) and differential settling that could jeopardize the stability of structures erected over these sites. Our findings highlight the potential use of the advocated integrated approach to assess the nature and extent of land deformation associated with rapid urban growth in arid lands, and to identify areas most impacted for the purpose of directing and prioritizing remediation efforts.

Countrywide Monitoring of Ground Deformation Using InSAR Time Series: A Case Study from Qatar

Over the past few decades the country of Qatar has been one of the fastest growing economies in the Middle East; it has witnessed a rapid increase in its population, growth of its urban centers, and development of its natural resources. These anthropogenic activities compounded with natural forcings (e.g., climate change) will most likely introduce environmental effects that should be assessed. In this manuscript, we identify and assess one of these effects, namely, ground deformation over the entire country of Qatar. We use the Small Baseline Subset (SBAS) InSAR time series approach in conjunction with ALOS Palsar-1 (January 2007 to March 2011) and Sentinel-1 (March 2017 to December 2019) synthetic aperture radar (SAR) datasets to assess ground deformation and conduct spatial and temporal correlations between the observed deformation with relevant datasets to identify the controlling factors. The findings indicate: (1) the deformation products revealed areas of subsidence and uplift with high vertical velocities of up to 35 mm/yr; (2) the deformation rates were consistent with those extracted from the continuously operating reference GPS stations of Qatar; (3) many inland and coastal sabkhas (salt flats) showed evidence for uplift (up to 35 mm/yr) due to the continuous evaporation of the saline waters within the sabkhas and the deposition of the evaporites in the surficial and near-surficial sabkha sediments; (4) the increased precipitation during Sentinel-1 period compared to the ALOS Palsar-1 period led to a rise in groundwater levels and an increase in the areas occupied by surface water within the sabkhas, which in turn increased the rate of deposition of the evaporitic sediments; (5) high subsidence rates (up to 14 mm/yr) were detected over landfills and dumpsites, caused by mechanical compaction and biochemical processes; and (6) the deformation rates over areas surrounding known sinkhole locations were low (+/−2 mm/yr). We suggest that this study can pave the way to similar countrywide studies over the remaining Arabian Peninsula countries and to the development of a ground motion monitoring system for the entire Arabian Peninsula.

Statistical Applications to Downscale GRACE-Derived Terrestrial Water Storage Data and to Fill Temporal Gaps

The Gravity Recovery and Climate Experiment (GRACE) has been successfully used to monitor variations in terrestrial water storage (GRACETWS) and groundwater storage (GRACEGWS) across the globe, yet such applications are hindered on local scales by the limited spatial resolution of GRACE data. Using the Lower Peninsula of Michigan as a test site, we developed optimum procedures to downscale GRACE Release-06 monthly mascon solutions. A four-fold exercise was conducted. Cluster analysis was performed to identify the optimum number and distribution of clusters (areas) of contiguous pixels of similar geophysical signals (GRACETWS time series); three clusters were identified (cluster 1: 13,700 km2; cluster 2: 59,200 km2; cluster 3: 33,100 km2; Step I). Variables (total precipitation, normalized difference vegetation index (NDVI), snow cover, streamflow, Lake Michigan level, Lake Huron level, land surface temperature, soil moisture, air temperature, and evapotranspiration (ET)), which could potentially contribute to, or correlate with, GRACETWS over the test site were identified, and the dataset was randomly partitioned into training (80%) and testing (20%) datasets (Step II). Multivariate regression, artificial neural network, and extreme gradient boosting techniques were applied on the training dataset for each of the identified clusters to extract relationships between the identified hydro-climatic variables and GRACETWS solutions on a coarser scale (13,700–33,100 km2), and were used to estimate GRACETWS at a spatial resolution matching that of the fine-scale (0.125° × 0.125° or 120 km2) inputs. The statistical models were evaluated by comparing the observed and modeled GRACETWS values using the R-squared, the Nash–Sutcliffe model efficiency coefficient (NSE), and the normalized root-mean-square error (NRMSE; Step III). Lastly, temporal variations in GRACEGWS were extracted using outputs of land surface models and those of the optimum downscaling methodology (downscaled GRACETWS) (Step IV). Findings demonstrate that (1) consideration should be given to the cluster-based extreme gradient boosting technique in downscaling GRACETWS for local applications given their apparent enhanced performance (average value: R-squared: 0.86; NRMSE 0.37; NSE 0.86) over the multivariate regression (R-squared: 0.74; NRMSE 0.56; NSE 0.64) and artificial neural network (R-squared: 0.76; NRMSE 0.5; NSE 0.37) methods; and (2) identifying local hydrologic variables and the optimum downscaling approach for individual clusters is critical to implementing this method. The adopted method could potentially be used for groundwater management purposes on local scales in the study area and in similar settings elsewhere.

Response of deep aquifers to climate variability

There is a general agreement that deep aquifers experience significant lag time in their response to climatic variations. Analysis of Temporal Gravity Recovery and Climate Experiment (GRACE), Soil Moisture and Ocean Salinity mission (SMOS), satellite altimetry, stable isotopic composition of groundwater, and precipitation and static global geopotential models over the Nubian Sandstone Aquifer System (NSAS) revealed rapid aquifer response to climate variability. Findings include: (1) The recharge areas of the NSAS (Northern Sudan Platform subbasin) witnessed a dry period (2002-2012), where average annual precipitation (AAP) was modest (85 mm) followed by a wet period (2013-2016; AAP: 107 mm), and during both periods the AAP remained negligible (<10 mm) over the northern parts of the NSAS (Dakhla subbasin); (2) the secular trends in terrestrial water storage (TWS) over the Dakhla subbasin were estimated at -3.8 ± 1.3 mm/yr and + 7.8 ± 1 mm/yr for the dry and wet periods, respectively; (3) spatial variations in TWS values and phase are consistent with rapid groundwater flow from the Northern Sudan Platform subbasin and Lake Nasser towards the Dakhla subbasin during the wet period and from the lake during the dry period; and (4) networks of densely fractured and karstified bedrocks provide preferential pathways for groundwater flow. The proposed model is supported by (1) rapid response in groundwater levels in distant wells (>280 km from source areas) and in soil moisture content in areas with shallow (<2 m) groundwater levels to fluctuations in Lake Nasser surface water, and (2) the isotopic composition (O, H) of groundwater along the preferred pathways, consistent with mixing of enriched (Lake Nasser water or precipitation over Sudan) and depleted (NSAS fossil water) endmembers. Findings provide new insights into the response of large, deep aquifers to climate variability and address the sustainability of the NSAS and similar fossil aquifers worldwide.

Agriculture Sprawl Assessment Using Multi-Temporal Remote Sensing Images and Its Environmental Impact; Al-Jouf, KSA

In this paper, multispectral and multi-temporal satellite data were used to assess the spatial and temporal evolution of the agriculture activities in the Al-Jouf region, Kingdom of Saudi Arabia (KSA). In the current study, an attempt was made to map the agriculture sprawl from 1987 to 2017 using temporal Landsat images in a geographic information system (GIS) environment for better decision-making and sustainable agriculture expansion. Our findings indicated that the agriculture activities developed through two crucial stages: high and low rise stages. Low rise stages occurred during three sub-stages from April 1987 to April 1988, from September 1993 to August 1998, and from April 2008 to May 2015, with overall change rates of 37.9, 44.4, and 30.5 km2/year, respectively. High rise stages occurred during three sub-stages from April 1988 to February 1993, from September 2000 to March 2006, and from April 2016 to August 2017, with overall change rates of 132.4, 159.1, and 119.5 km2/year, respectively. Different environmental problems due to uncontrolled agriculture activities were observed in the area, including substantial depletion of the groundwater table. Another environmental impact observed was the appearance of sinkholes that occurred suddenly with no warning signs. These environmental impacts will increase in the future if no regulated restrictions are implemented by decision-makers.