The Gravity Recovery and Climate Experiment (GRACE) has been successfully used to monitor variations in terrestrial water storage (GRACE_TWS) and groundwater storage (GRACE_GWS) 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 (GRACE_TWS time series); three clusters were identified (cluster 1: 13,700 km²; cluster 2: 59,200 km²; cluster 3: 33,100 km²; 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, GRACE_TWS 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 GRACE_TWS solutions on a coarser scale (13,700–33,100 km²), and were used to estimate GRACE_TWS at a spatial resolution matching that of the fine-scale (0.125° × 0.125° or 120 km²) inputs. The statistical models were evaluated by comparing the observed and modeled GRACE_TWS 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 GRACE_GWS were extracted using outputs of land surface models and those of the optimum downscaling methodology (downscaled GRACE_TWS) (Step IV). Findings demonstrate that (1) consideration should be given to the cluster-based extreme gradient boosting technique in downscaling GRACE_TWS 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.

Persistent scatterer interferometric analyses were conducted on a stack of 84 Environmental Satellite's Advanced Synthetic Aperture Radar scenes spanning 7 years (2004 to 2010) over the entire Nile Delta of Egypt and surroundings (area: 40,416 km²) to monitor the ongoing spatial and temporal land deformation, identify the factors controlling the deformation, and model the interplay between sea level rise and land subsidence to identify areas and populations threatened by sea encroachment by the end of the 21st century. Findings include the following: (1) general patterns of subsidence (average rate: −2.4 mm/year) in the northern delta, near‐steady to slight subsidence in the southern delta (average rate: 0.4 mm/year), separated by a previously mapped flexure zone (minimum width: 20–40 km) undergoing uplift (average rate: 2.5 mm/year); (2) high subsidence rates (up to −8.9 mm/year) over the north central and northeastern delta (area: ~4,815 km²), possibly due to compaction of recent (5 m) silt and clay‐rich Holocene sediments; (3) high subsidence rates (up to −9.7 mm/year) in areas where the highest groundwater extraction rates were reported in southern delta (Menoufia governorate) and in reclaimed desert land in the western delta (Beheira governorate); (4) high subsidence rates (up to −9.7 mm/year) over onshore gas fields, notably the Abu Madi gas field, where high gas extraction rates have been recorded; and (5) using extracted deformation rates, high‐resolution TanDEM‐X digital elevation model, a eustatic sea level rise of 0.44 m, and applying a bathtub inundation model, an estimated 2,660 km² in northern delta will be inundated by year 2100.

An integrated approach [field, Interferometric Synthetic Aperture Radar (InSAR), hydrogeology, geodesy, and spatial analysis] was adopted to identify the nature, intensity, and spatial distribution of deformational features (sinkholes, fissures, differential settling) reported over fossil aquifers in arid lands, their controlling factors, and possible remedies. The Lower Mega Aquifer System (area 2 × 10⁶ km²) in central and northern Arabia was used as a test site. Findings suggest that excessive groundwater extraction from the fossil aquifer is the main cause of deformation: (1) deformational features correlated spatially and/or temporally with increased agricultural development and groundwater extraction, and with a decline in water levels and groundwater storage (− 3.7 ± 0.6 km³/year); (2) earthquake events (years 1985–2016; magnitude 1–5) are largely (65% of reported earthquakes) shallow (1–5 km) and increased from 1 event/year in the early 1980s (extraction 1 km³/year), up to 13 events/year in the 1990s (average annual extraction > 6.4 km³). Results indicate that faults played a role in localizing deformation given that deformational sites and InSAR-based high subsidence rates (− 4 to − 15 mm/year) were largely found within, but not outside of, NW–SE-trending grabens bound by the Kahf fault system. Findings from the analysis of Gravity Recovery and Climate Experiment solutions indicate that sustainable extraction could be attained if groundwater extraction was reduced by 3.5–4 km³/year. This study provides replicable and cost-effective methodologies for optimum utilization of fossil aquifers and for minimizing deformation associated with their use.

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.

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.

The normalized difference vegetation index (NDVI) has been frequently used to map hail damage to vegetation, especially in agricultural areas, but observations can be blocked by cloud cover during the growing season. Here, the European Space Agency’s Sentinel-1A/1B C-band synthetic aperture radar (SAR) imagery in co- and cross polarization is used to identify changes in backscatter of corn and soybeans damaged by hail during intense thunderstorm events in the early and late growing season. Following a June event, hail-damaged areas produced a lower mean backscatter when compared with surrounding, unaffected pixels [vertical–vertical (VV): −1.1 dB; vertical–horizontal (VH): −1.5 dB]. Later, another event in August produced an increase in co- and cross-polarized backscatter (VV: 0.7 dB; VH: 1.7 dB) that is hypothesized to result from the combined effects of crop growth, change in structure of damaged crops, and soil moisture conditions. Hail damage regions inferred from changes in backscatter were further assessed through coherence change detections to support changes in the structure of crops damaged within the hail swath. While studies using NDVI have routinely concluded a decrease in NDVI is associated with damage, the cause of change with respect to the damaged areas in SAR backscatter values is more complex. Influences of environmental variables, such as vegetation structure, vegetation maturity, and soil moisture conditions, need to be considered when interpreting SAR backscatter and will vary throughout the growing season.