Satellite-Based Groundwater Monitoring: History
Please note this is an old version of this entry, which may differ significantly from the current revision.
Subjects: Water Resources
Contributor:

Satellite is an effective technology to assess water storages remotely through remote sensing and geographic information system (RS and GIS), which are very useful tools in groundwater exploration mapping. They provide inputs to assess the total groundwater resources in a region and for the selection of suitable sites for drilling or artificial recharge.

  • Satellites
  • groundwater
  • GRACE

1. Most Common Satellites

The first application of satellite data for earth monitoring, including hydrological studies, started back in 1957 [61,62]. Satellite data of LANDSAT, MODIS, and other similar satellites (Sentinel-2, SPOT etc.) have been used for the assessment of cryosphere (snow and glaciers), soil water, surface water, groundwater, evaporation, and other constituents of the global hydrologic system [63,64,65,66].
The Moderate Resolution Imaging Spectrometer (MODIS) data are used to help improve the understanding of global dynamics and processes that occur on the land, in the oceans, and in the lower atmosphere [67]. These specifically include total precipitable water, cloud, atmospheric profiles, land cover, evapotranspiration, water mask, ocean products, and snow, glaciers, and sea ice cover [16,67,68,69]. MODIS is a sensor that operates on the two satellites, i.e., Terra and Aqua, which were launched by NASA in 1999 and 2002, respectively [70]. Terra passes around the earth from north to south across the equator in the morning, whereas the Aqua satellite passes from south to north in the afternoon. In this way, both satellites view the entire Earth’s surface every one to two days, obtaining data in 36 spectral bands. The spatial resolution of MODIS data is 250 m, 500 m, and 1 km [71]. Relatively coarse resolution is its limitation; however, its benefit is its open-access availability.
MODIS satellite does not assess ground water resources directly; however, it helps to improve the different model inputs for groundwater modeling. An example consists of using the Nation-wide Groundwater Recharge Model (NGRM) for New Zealand. It uses MODIS-derived evapotranspiration and vegetation satellite data, available rainfall dataset, elevation, soil and geology, etc., and this research shows that with minor model adjustments and use of improved input data, large-scale models and satellite data can be used to derive rainfall recharge estimates for local scales. The estimated recharge of the NGRM model compares well to most local and regional data and recharge models [72].
Landsat satellite is also a helpful tool for the groundwater assessment; however, it has high resolution (i.e., 30 m) [73] relative to MODIS. Its images are also freely available for research purposes, which is a great benefit. The Landsat Program is a continuous Earth-observing satellite mission, which is jointly managed by NASA and the U.S. Geological Survey. It started in 1972 with Landsat 1 and has been continued throughout the following four decades with Landsat 2, 3, 4, 5, 6, 7, and 8, and the last one (Landsat 9) launched in mid-2021 [74].
Landsat data collect spectral information from earth’s surface. Landsat data have evolved as a valuable resource of agriculture, land use, water resources, forestry, and natural resource investigation along with the understanding of existing circumstances of land changes in fresh water supplies [73,75,76,77].
The European Space Agency (ESA) under the Copernicus Program provides Sentinel-2 with temporal and spatial (10–20 m) resolution of 5 days and 10–20 m, respectively, which has opened new vistas for many applications for having higher resolution than both MODIS and Landsat [78,79,80,81,82,83,84,85]. Sentinel-2 is a multispectral operational imaging mission for worldwide land observation. It provides the data on soil, water, and vegetation cover for land, inland waterways, and coastal areas [86]. Dubois et al. [87] used Sentinel-2 data to estimate the groundwater withdrawals in the Maghreb area. Some other studies also estimated groundwater needs indirectly via determining irrigation requirements for agricultural production [81,88].
Interferometric Synthetic Aperture Radar (INSAR) makes high-density measurements over large areas by using radar signals through Earth-orbiting satellites for measuring the variations of land surface elevation at higher resolution and spatial detail [89]. With INSAR, the variations in land surface elevation, the spatial extent, and the magnitude of deformation that is linked with fluid abstraction and natural hazards (landslides, volcanoes, earthquakes) are measured effectively. By utilizing synthetic aperture radar satellite data through InSAR techniques, ground deformation is assessed [90,91]. For local or regional scales, InSAR maps millimeter-scale displacements of the Earth’s surface with radar satellite measurements. It helps to survey remote sites quickly by measuring areas (thousands of kilometers) and attaining a spatial resolution of around 20 m.
There can be a combination of multiples data sources and different techniques to find a particular parameter of groundwater. Srivastava and Bhattacharya [92] utilized IRS-LISS III (Indian Remote Sensing 1D Linear Imaging Self-Scanning Sensor 3) with TM/ETM + (Landsat Thematic Mapper/Enhanced TM) digital data and SRTM DEM (Shuttle Radar Topographic Mission Digital Elevation Models) to create multiple models. These include geomorphology, lithology, land uses, soil, lineament, river gradient, drainage density, and slope maps for detailed groundwater modeling.

2. GRACE Satellite Data

The direct measurement of water column remained unsolved till 2002, when Gravity Recovery and Climate Experiment (GRACE), consisting of two satellites, was launched. These two satellites were upgraded to GRACE-FO in 2018. The combined task of the U.S. and German space agencies aims for detailed measurements of Earth’s gravity field changes. These have revolutionized the explorations about Earth’s water reservoirs over land, ice, and oceans along with earthquakes and crustal deformations. The data of this project are free for research and management purposes. GRACE offers monthly, vertically combined estimates of Terrestrial Water Storage (TWS) anomalies (departures from the long-term mean) at a coarse spatial resolution of about 300 km. TWS is comprised of groundwater, soil water, surface water, snow, and ice. GRACE estimations are used to assess the groundwater-depletion rates throughout the world [93]. Although the resolution of this satellite is coarse, it is very useful for the large basins on a regional or global scale.

2.1. GRACE Products

The data obtained from the GRACE provide the information about the gravity anomalies. It needs extensive processing to extract groundwater-related direct information. The GRACE data can be classified into three stages/levels. The raw data, collected from satellites (Level 0), is calibrated and time tagged in a non-destructive (or reversible) sense and then labeled as Level 1A. Level 1A data products are not disseminated directly to the public. The data then undergo extensive and irreversible processing and are transformed to clean data products, labeled as Level 1B. These products include the inter-satellite range, range-rate, range-acceleration, the non-gravitational accelerations from each satellite, the pointing estimates, the orbits, etc. The Level 1B products are processed to produce the monthly gravity field estimates in form of spherical harmonic coefficients. These estimates are labeled Level 2. Occasionally, several months of data are combined to produce an estimate of the mean or static gravity field. After validation, all Level 2 and associated Level 1B products are released to the public.
GRACE Level 3 data represent monthly land-mass grids that possess terrestrial anomalies because of aquifers, snow covers, and river basins, etc. [94]. From the GRACE project data products, many GRACE users have put together the resources to generate and distribute value-added (or Level 3) products [95]. For the ease of users who would prefer to access GRACE data products as mass anomalies (for example water layer), some standardized products and tools are available from other sources, which are called Level 3 Data Products [96]. Figure 1 describes the GRACE data levels, processing, and the characteristics of these products.
Figure 1. GRACE Satellite Data Levels, characteristics, and processing.
Spherical harmonic (SH) (Level 2) and mascon (Level 3) solutions are the two main GRACE products used for the estimation of TWS. Both techniques have successfully been used in the retrieval of time-varying gravity fields from GRACE. For SH products, filtering and leakage corrections are required to reduce the noise, while no such correction is required for the mascon solutions. The capability of SH and mascon solutions to determine basin level mass variation is computed with an assessment of how the noise and errors are handled, which are inherent in GRACE solution. Although mascon solutions have an edge over the SH by avoiding leakage correction, the performance of SH for local scales is better than the mascon if corrections are applied correctly.
Leakage error is one of the many GRACE error constituents (for example, measurement and post-processing error). Generally, a leakage adjustment is needed to recover biases. The most extensively used approach is forward modeling, which uses a priori information (generally from global model simulations or in situ observations) to mimic the GRACE data processing techniques. Uncertainties of GRACE data is another main issue. Detailed analysis of uncertainties is important to achieve reliable and optimized products for obtaining TWS changes globally. For this purpose, different statistical techniques/approaches are used. For example, to enumerate uncertainties in the changes of TWS in GRACE observations, the three-cornered hat technique is in practice. Amid all other TWS change products, the BMA-based changes in TWS demonstrate the maximum consistency using with the WGHM output [97]. Empirical Orthogonal Function decomposition technique can be used to show detailed patterns of TWS variations annually and seasonally to see the long-term trends [97]. Additionally, Monte Carlo simulation framework is effective for reduction of uncertainty and improved parameter estimation while using GRACE in the calibration of Land Surface Models (LSM) parameters [22,98,99].

2.2. Estimation of Groundwater Storage (GWS)

Due to absence of long observational data and the sparse spread of monitoring wells, satellite-based observations of groundwater-storage variations through the GRACE satellite are often used to estimate the groundwater-storage changes [17,100,101]. Terrestrial water storage includes all landforms of water available to an area. It includes groundwater, surface water and canopy water, soil moisture, and snow and is a main component of the global and continental hydrologic cycle. Groundwater-storage (GWS) changes are acquired from TWS through deducting soil moisture, surface water, canopy water, and snow water equivalent data gained from GLDAS or with any other dataset [26,102,103,104]. Figure 2 explains the components of freshwater storage and processes. The same is explained in the following Equations:
Figure 2. GRACE based total water storage variations are a composite of various continental water storage compartments.
Since 2002, the GRACE has supplied a highly valued dataset, which allows the study of TWS over large river basins globally. Shamsudduha et al. [105] used GRACE datasets for a long time period and found that for the recent time period (2003 to 2007), groundwater-storage changes correlate well with in situ borehole records in the Bengal Basin in Bangladesh. For the Indus Basin, the variation in groundwater storage at numerous spatial scales has been studied by Iqbal et al. [106] with GRACE data. They found that the groundwater storage changes through GRACE-based estimation is reasonably sufficient for the provision of monthly groundwater-storage changes.
Variations in TWS can be acquired through the GRACE data. After the changes in groundwater, the variation of soil moisture is normally the highest factor in TWS changes and intermediate zone storage [107]. Strassberg et al. [108] found a robust link between the sum of soil moisture and groundwater storage with TWS by means of GRACE data. Before GRACE satellite, estimates of TWS changes were used to derive from the in situ data of soil moisture and snow. These approaches cover analyzing the variations in groundwater, soil water, surface water, and snow through utilizing a long time series [107].
Since the GRACE satellite was launched in 2002, a great deal of researchers has conducted the studies of groundwater using GRACE satellite data. According to the research articles published in various journals (source: https://scholar.google.com/, dated 18 January 2022), an annual frequency graph of articles has been developed to envisage the work carried out using GRACE data (refer to Figure 3).
Figure 3. Annual Article Frequency on GRACE.

2.3. GRACE Applications

GRACE provides useful information about global groundwater depletion. The shallow GWS declines faster than the deep GWS [109]. Globally, GRACE could be the only hope for the assessment of groundwater reduction [104]. GRACE data are useful for the estimation of regional groundwater monitoring on monthly to seasonal scale. Yeh et al. [110] confirmed the ability of GRACE to simulate regional groundwater-storage changes from monthly to seasonal scale for an area of about 200,000 km2 area in Illinois. It also provides reasonable results for the heterogeneous GWS variations as observed by Huang et al. [109]. Döll et al. [111] conducted groundwater storage monitoring for the Mississippi Basin. However, the study failed to obtain reliable results to monitor the total water storage due to water withdrawals at the whole scale of the basin. GRACE data were applied in global hydrological investigations, including water storage change evaluation, groundwater and evapotranspiration retrieving, drought analysis, and glacier response to global change [112]. In another study, complementary relationship (CR) approach was used to assess the evapotranspiration and a water balance method to estimate the runoff and base flow in the Conterminous United States (CONUS), which demonstrates how to use it with limited parameterization and datasets [113].
There are many GRACE-based studies of groundwater that showed substantial aquifer depletion for the large regions, for example, the North China, the northwest India, the Middle East, the Murray Darling Basin of Australia, and California’s Central Valley aquifers of the USA [17,101,114,115,116,117]. The shrinking phenomena of groundwater-fed lakes have also been monitored through GRACE and GLDAS [16]. Sun et al. [118] developed a method for estimation of aquifer storage parameters (specific yield and storability) utilizing the GRACE data for the connected aquifers of Edwards-Trinity Plateau and Pecos Valley in central Texas. Frappart and Ramillien [60] also found that these data were used to evaluate the parameters; the specific yield, for example, is used to relate groundwater level to storage or to describe the indices of groundwater depletion and stress. In another GRACE-GLDAS study, it was found that high consistency exists between GRACE observations and changes in water wells [119].

2.4. Additional Computational Tools

For groundwater management, the usage of GRACE data at a local scale has been inadequate due to in-built uncertainties in GRACE data and problems in disintegrating several TWS components.
Artificial Neural Network (ANN) models can be developed to foresee the variations of groundwater levels directly by means of a gridded GRACE product and other openly existing hydro-meteorological datasets. ANN is used as a statistical downscaling method, which is extensively used in streamflow projecting and management of water resources [120,121,122,123,124,125,126]. A trained ANN model is comparatively strong to data noise and is quite suitable to support real-time decision making [127,128]. Ensemble ANN models can be utilized to forecast monthly and seasonal groundwater-level changes for various wells situated in certain areas.
Because of the GRACE data, the ANN models forecast water-level variations instead of absolute water levels [128]. By combining GRACE, TRMM, and GLDAS, ANN predicts and evaluates groundwater drought [129]. Gemitzi et al. [130] used ANN to perform statistical downscaling of GRACE data using local meteorological data. A typical application process of ANN is shown in Figure 4.
Figure 4. Predicted downscaled GRACE-GLDAS GWS changes with ANN.
Srivastava and Bhattacharya integrated numerous thematic plots on the GIS platform based on Saaty’s Analytical Hierarchy Process (AHP) with the Groundwater Potential Index (GWPI) values (0.175 to 0.940) to produce a good match with the resistivity and pumping test data [92]. Similarly, AHP was used by Mallick et al. [41] to find the most important contributing features to delineate the groundwater potential zones [131].

3. Supporting Terrestrial Modeling Systems

To monitor the other constituents of terrestrial water storage, for example, soil water, surface water, snow and ice water, canopy water, and groundwater, separately, some system is required to address this aspect, which was fulfilled by upgrading the GRACE system. This can be done with Global Land Data Assimilation System (GLDAS) and WaterGAP Global Hydrological Model (WGHM), which help GRACE to disintegrate the TWS into its required water component. For example, TWS disintegrates with help of these models to give groundwater anomalies.

3.1. Global Land Data Assimilation System (GLDAS)

The Global Land Data Assimilation System (GLDAS) is a worldwide, high-resolution, offline terrestrial modeling system that integrates satellite and ground-based observations to produce optimum fields of land surface conditions and fluxes in near-real time.
It produces a sequences of worldwide land surface conditions (for example, soil moisture and surface temperature) and fluxes (for example, evaporation and sensible heat flux) products that were simulated by four Land Surface Models (LSM), that is, VIC, Mosaic, CLAM, and Noah.
Several current and proposed weather and climate prediction, water resources applications, and water-cycle investigations rely on GLDAS output. GLDAS runs globally and on user defined domains at resolutions ranging from 2.5° to 1 km. The GRACE and GLDAS simulations provide precise and reliable datasets that are used to describe the variations of groundwater storage regionally [132]. The assessed groundwater variations by GRACE-GLDAS are usually in line with in situ observations [22,97,98,109,110,132,133].

3.2. WaterGAP Global Hydrological Model (WGHM)

WGHM calculates timeseries of surface and subsurface runoff, groundwater recharge, and river flows along with water storage variations in canopy, snow, soil, groundwater, rivers, lakes, and wetlands. Therefore, it can quantify the whole water resources along with the groundwater resources. Institute of Physical Geography (IPG) developed WGHM to provide freshwater fluxes with storage data on the land surface globally [134,135]. The spatial resolution of the model is 0.5° × 0.5°, and computations are carried out on daily basis, whereas the output is analyzed on a monthly time step. The WGHM considers the human-induced water consumption along with the daily groundwater simulations [134]. Figure 5 explains the possible approaches to incorporated GLDAS and WGHM data to obtain GWS from the GRACE data. In this figure, it can be seen that TWS is obtained through GRACE satellite, whereas both GLDAS and WGHM give surface water, soil moisture, snow water equivalent, and groundwater storage. Final groundwater storage profiles are obtained through the interaction of GRACE with GLDAS and GRACE with WGHM, which can be compared with in situ groundwater measurements based upon the statistical selection criteria as given in the figure.
Figure 5. GRACE–GLDAS and GRACE–WGHM processing to get GWS from TWS.
This entry is offline, you can click here to edit this entry!
Video Production Service