Next Article in Journal
Assessing Leaf Biomass of Agave sisalana Using Sentinel-2 Vegetation Indices
Next Article in Special Issue
Groundwater Depletion Signals in the Beqaa Plain, Lebanon: Evidence from GRACE and Sentinel-1 Data
Previous Article in Journal
Sugarcane Yield Mapping Using High-Resolution Imagery Data and Machine Learning Technique
Previous Article in Special Issue
GRACE—Gravity Data for Understanding the Deep Earth’s Interior
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Assessment of Hydrologic- and Flood-Induced Land Deformation in Data-Sparse Regions Using GRACE/GRACE-FO Data Assimilation

by
Natthachet Tangdamrongsub
1,2,* and
Michal Šprlák
3
1
Earth System Science Interdisciplinary Center, University of Maryland, College Park, MD 20740, USA
2
Hydrological Sciences Laboratory, NASA Goddard Space Flight Center, Greenbelt, MD 20771, USA
3
NTIS—New Technologies for the Information Society, Faculty of Applied Sciences, University of West Bohemia, 306 14 Plzeň, Czech Republic
*
Author to whom correspondence should be addressed.
Submission received: 17 December 2020 / Revised: 3 January 2021 / Accepted: 8 January 2021 / Published: 12 January 2021
(This article belongs to the Special Issue GRACE Satellite Gravimetry for Geosciences)

Abstract

:
The vertical motion of the Earth’s surface is dominated by the hydrologic cycle on a seasonal scale. Accurate land deformation measurements can provide constructive insight into the regional geophysical process. Although the Global Positioning System (GPS) delivers relatively accurate measurements, GPS networks are not uniformly distributed across the globe, posing a challenge to obtaining accurate deformation information in data-sparse regions, e.g., Central South-East Asia (CSEA). Model simulations and gravity data (from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GRACE-FO)) have been successfully used to improve the spatial coverage. While combining model estimates and GRACE/GRACE-FO data via the GRACE/GRACE-FO data assimilation (DA) framework can potentially improve the accuracy and resolution of deformation estimates, the approach has rarely been considered or investigated thus far. This study assesses the performance of vertical displacement estimates from GRACE/GRACE-FO, the PCRaster Global Water Balance (PCR-GLOBWB) hydrology model, and the GRACE/GRACE-FO DA approach (assimilating GRACE/GRACE-FO into PCR-GLOBWB) in CSEA, where measurements from six GPS sites are available for validation. The results show that GRACE/GRACE-FO, PCR-GLOBWB, and GRACE/GRACE-FO DA accurately capture regional-scale hydrologic- and flood-induced vertical displacements, with the correlation value and RMS reduction relative to GPS measurements up to 0.89 and 53%, respectively. The analyses also confirm the GRACE/GRACE-FO DA’s effectiveness in providing vertical displacement estimates consistent with GRACE/GRACE-FO data while maintaining high-spatial details of the PCR-GLOBWB model, highlighting the benefits of GRACE/GRACE-FO DA in data-sparse regions.

Graphical Abstract

1. Introduction

Land deformation is a response of (1) geophysical processes governed by, e.g., hydrology, atmosphere, ocean mass transports associated with the earth system, and (2) anthropogenic impact associated with the human intervention [1,2]. Such processes cause the Earth’s surface to subside or uplift, severely affecting the scientific and industrial sectors [3,4]. Despite being affected by different Earth system drivers, land deformation is dominated by the hydrology cycle at annual and interannual scales [1,5,6,7]. For instance, vertical motions’ annual variations in Japan [8] and the Northwestern United States [9] are induced by snow loadings. Similarly, seasonal variations in California are explained by groundwater recharge and pumping activities [10]. In monsoon regions, significant terrestrial water storage variation (ΔTWS), influenced by exceptional annual precipitation (e.g., >1000 mm), is a primary driver of the land deformation’s annual variations (e.g., [11,12,13]). Monsoon flooding can also pose significant loads over inundated areas, causing rapid vertical displacement changes (or land subsidence) of a few mm [5,14]. Despite considering hydrologic loads as a signal or interference [1], accurate deformation measurements are needed to gain constructive insight into regional geophysical processes.
Measuring land deformation using Global Positioning System (GPS) networks has been used to achieve <1-cm measurement’s accuracy [6,15,16,17]. However, this technique requires moderately dense GPS sites to operate continuously. This is impractical in many parts of the world concerning, e.g., high operational/maintenance costs and limited accessibility (in remote areas), ending up with sparse GPS networks that can hardly be used to monitor regional (large-scale) land deformation rigorously. In data-sparse regions, satellite gravity data have been considered a primary measurement regarding multiple advantages such as continuous time series, extended ground coverage, and no operational cost required from users. The Gravity Recovery and Climate Experiment (GRACE [18]) and its successor, GRACE-Follow On (GRACE-FO [19]), satellite missions have provided accurate monthly earth gravity field measurements in various applications, including water resource, land deformation, and climate [17,20,21,22,23]. For convenience, the term GRACE is used for both GRACE and GRACE-FO throughout this paper. The successful measuring of large-scale land deformation by GRACE has been reported in various geographic and climate conditions, e.g., the Himalayas and India [17], Africa [11], and Australia [24]. The downside of GRACE data is the relatively coarse spatial resolution (i.e., ~300 × 300 km2) that may degrade the derived local deformation’s accuracy [6].
Estimating land deformation from land surface or hydrology models can be alternatively considered [5,25,26], particularly to achieve higher spatial resolution than GRACE or improved data availability under an absent gravity data’s scenario. Models implement sophisticated physics to represent the water cycle, which can provide an accurate TWS estimate comparable with GRACE results [21]. However, despite the overall positive impacts, model estimates commonly suffer from severe uncertainties associated with, e.g., inadequate physics representation, ineffective model configuration, or inaccurate metrological forcing [27], leading likely to inaccurate or underestimated TWS [25,28]. GRACE data assimilation (DA) can be used to improve model TWS estimates [27]. GRACE DA adopts a sequential DA approach to assimilate GRACE observations into a model simulation according to their relative errors [29], resulting in (updated) model TWS estimates consistent with the assimilated observations [30,31]. This approach has been successfully used in several hydrologic applications, e.g., drought assessment [31] and groundwater monitoring [32], while its benefit to the land deformation estimate has only been considered recently [26].
This study aims to assess the performance of GRACE DA-derived vertical displacement in data-sparse regions between April 2002 and December 2018. Central South-East Asia (CSEA) is used as a tested region, considering its monsoon climate and frequently experienced severe flooding that can be used to investigate the hydrologic- and flood-induced land deformation behavior [12,33]. Importantly, ground measurements from a few GPS stations in CSEA are available to validate the results. Data assimilation is performed by assimilating GRACE observations into the PCRaster Global Water Balance (PCR-GLOBWB [34]) model using the Ensemble Kalman Smoother (EnKS [30]) to improve the TWS and vertical displacement estimates. PCR-GLOBWB simulates the TWS at five-arcmins (~10 km) resolution, allowing much greater spatial details to be recovered. The evaluation is performed by comparing vertical displacement estimates (from GRACE, PCR-GLOBWB, and GRACE DA) with GPS measurements in terms of the correlation coefficient (ρ), root mean square difference (RMSD), and RMS reduction [35] to assess the benefits/limitations of different approaches.
The objectives of this study are threefold: (1) to estimate land deformation in data-sparse regions (here is CSEA) using GRACE, PCR-GLOBWB, GRACE DA, and GPS data; (2) to evaluate the GRACE DA’s performance (as well as GRACE and PCR-GLOBWB); and (3) to verify the primary driver of land deformation in CSEA and quantify the vertical displacement changes during flooding events. The presentation of this paper follows these objectives. Section 2 provides the details of the study area and data processing approaches. Section 3 describes the methods for vertical displacement computation and GRACE DA. The results are presented in Section 4, where the benefits of GRACE and GRACE DA in data-sparse regions are elaborated. Section 5 benchmarks the GRACE DA’s performance, and Section 6 concludes the findings of this paper.

2. Study Area and Data Processing

2.1. Study Area

This study focuses on Central South-East Asia (CSEA; 11° N–19° N, 97° E–106° E; Figure 1). The study area stretches over 800,000 km2, encompassing five different countries (i.e., Myanmar, Thailand, Laos, Cambodia, and Vietnam) and two major river basins (i.e., the Chao Phraya River and Tonlé Sap Basins). CSEA has a tropical wet climate with an average temperature and annual rainfall of approximately 25°C and 1800 mm, respectively. Essential sources of water supply in most countries are lakes and reservoirs, in which the surface water contributes significantly to regional water storage dynamics, e.g., >40% over the Tonlé Sap Basin [36]. The region has experienced several episodes of extreme monsoon rainfalls and flooding [37]. A recent severe flooding was seen in Thailand in 2011, where the inundated areas extended over 18,000 km2 [38]. Substantial hydrologic mass variations in CSEA induce a seasonal vertical motion that can likely be measured by a GPS network [33]. However, GPS measurements are very sparse in CSEA. Only six GPS sites’ measurements have collected data for longer than three years, and only two of them are still operational (see Table 1). This data-sparse problem limits the comprehensive spatiotemporal deformation details to be observed. The use of satellite observations or model estimates might help improve the deformation monitoring ability in CSEA.

2.2. GRACE Data

The GRACE and GRACE-FO Level 2 release 06 (RL06) products provided by the Center for Space Research (CSR) at the University of Texas at Austin (http://www2.csr.utexas.edu/grace/RL06.html and http://www.csr.utexas.edu/grace-fo) are used to derive the ΔTWSs and of vertical displacement variations (Δr) between April 2002 and December 2018. The products contain monthly (geopotential) spherical harmonic coefficients (SHC) up to a maximum degree and order of 96. The GRACE data are processed as follows. First, degree 1 coefficients are restored using values provided by Swenson et al. [39]. C20 and C30 coefficients are replaced by the value estimated from the satellite laser ranging [40]. Second, the long-term SHC mean (between April 2002 and December 2018) is computed and removed from each monthly product to obtain monthly SHC variations. Third, destriping [41] and 300-km radius Gaussian smoothing filters [42] are applied to the SHC variations to alleviate high-frequency errors. Fourth, the ΔTWS or Δr is computed from the filtered SHC variations (see also Section 3.1), and the Glacial Isostatic Adjustment correction provided by Caron et al. [43] is applied. Finally, a signal restoration approach [36,44] is performed to restore the damped signal caused by the applied filters. It is noteworthy that the signal restoration used in this study is applied to the TWS (not to displacements directly). The restoration approach searches for genuine GRACE signals through an iteration process, and the determined SHCs are used to compute ΔTWS and Δr (see Section 3 for more detail).
In contrast to the GPS measurements, the GRACE Level 2 data do not include short-term atmospheric and nontidal oceanic variations. Therefore, the de-aliasing atmospheric and oceanic signals provided by the GRACE Atmosphere and Ocean De-aliasing RL06 product [45] are added to the GRACE-derived vertical displacement when comparing Δr estimates with GPS measurements.

2.3. Hydrology Model

The global distributed hydrological model, PCR-GLOBWB Version 2 [34,46], is a grid-based hydrology model that can simulate various water storage components, e.g., soil moisture, groundwater, surface water, and snow. The model has spatial and temporal resolutions of five arcmins and one day, respectively. PCR-GLOBWB-simulated TWS comprises 27 different water storage components, i.e., eight soil moisture, four interceptions, eight snow, four inundated topwater, one surface water, and two groundwater. The complete details of PCR-GLOBWB can be found in, e.g., [34,46,47,48].
To perform a PCR-GLOBWB simulation, we use daily precipitation and 2-m air temperature from the Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG) [49] and the hourly land product of the European Centre for Medium-range Weather Forecasts Reanalysis Version 5 (ERA5; https://confluence.ecmwf.int/display/CKB/ERA5-Land), respectively. Both products have a spatial resolution of 0.1°, and they are resampled to five arcmins to overlay with the model grid space. The ERA5′s hourly data is resampled to a daily interval, consistent with the model temporal resolution. The daily reference potential evapotranspiration is calculated using the Hamon method [50].

2.4. GPS Measurements

GPS data are obtained from the Nevada Geodetic Laboratory at the University of Nevada, Reno (UNR) [51]. A final solution provides an hourly vertical motion associated with the IGS14 reference frame. Multiple GPS sites are located in CSEA (see http://geodesy.unr.edu/NGLStationPages/gpsnetmap/GPSNetMap.html) but most contain significant data gaps, i.e., only a few records available in the time series. In total, only six GPS sites have sufficiently long records (i.e., >3 years) for our analysis (Table 1). At each GPS site, the hourly records of the up (vertical) component are averaged into a monthly averaged time series, consistent with the GRACE temporal resolution. The long-term trend is removed from the time series to reduce the plate tectonic effects [7]. The trend is estimated by fitting the GPS time series using the least-squares adjustment associated with the offset, trend, annual, and semi-annual parameters. Then, the long-term mean is subtracted from the time series to obtain GPS position variations. Finally, the measurements of GPS sites inside the same model grid cell (i.e., five-arcmins grid box) are averaged before performing the validation. This averaging process results in three GPS time series associated with three independent model grid cells (S1–S3; see Table 1).

2.5. Flood Extent

The flood extent is derived from surface reflectance measurements from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensors onboard NASA’s Terra satellites. The MOD09A1 collection 6 product [52] provides surface reflectance in 7 different frequency bands at 500-m spatial resolution within an 8-day window. The combinations of specific frequency bands can be used to identify open water bodies larger than 0.25 km2. In this study, the multifrequency surface reflectance is extracted from two MODIS tiles, h27v07 and h28v07 (covering the entire study area), between April 2002 and December 2018. The pixels with filled values or severe cloud cover are masked, using data quality control flags provided with the MODIS product. The pixel associated with surface water can be classified using the normalized difference water index (NDWI) computed using measurements from green and near-infrared (NIR) channels as follows [53]:
NDWI   = green     NIR green   +   NIR
The NDWI values are between −1 and 1, where positive values represent open water, and zero or negative values represent clouds, vegetation canopy, or soil. The flooded extent is computed by multiplying the total number of open water pixels (NDWI > 0) by 0.25 km2 (the spatial resolution of the surface reflectance dataset).

3. Methods

The computation and analysis follow three steps: (1) data preparation (see Section 2), (2) vertical displacement computation, and (3) validation and accuracy assessment. Figure 2 illustrates the data processing diagram.

3.1. Computation of Terrestrial Water Storage and Vertical Displacement Variations

The GRACE monthly geopotential spherical harmonic coefficients (SHC), together with the load Love numbers [7], are used to compute the ΔTWS and Δr as follows [54]:
Δ TWS ( θ , λ ) = a ρ ave 3 ρ w n = 1 N max m = 0 n P ¯ nm ( cos θ ) ( Δ C ¯ nm cos ( m λ ) + Δ S ¯ nm sin ( m λ ) ) 2 n + 1 1 + k n
Δ r ( θ , λ ) = a n = 1 N max m = 0 n P ¯ nm ( cos θ ) ( Δ C ¯ nm cos ( m λ ) + Δ S ¯ nm sin ( m λ ) ) h n 1 + k n
where ( θ , λ ) are the colatitude and east longitude, ( k n , h n ) are the load Love numbers of degree n , P ¯ nm is the normalized associated Legendre function of the first kind of degree n and order m , ( ρ ave , ρ w ) are the average density of the Earth and freshwater, a is the radius of the Earth, and ( Δ C ¯ nm , Δ S ¯ nm ) are the SHC variation, computed by removing a long-term mean SHC (e.g., calculated between April 2002 and December 2018) from the monthly ( C ¯ nm , S ¯ nm ) . The load Love numbers computed from the Preliminary Reference Earth Model (PREM) [55] are used in this study.
Vertical displacement associated with model-simulated ΔTWS can also be computed using Equation (3). However, the model provides ΔTWS in a geographical gridded format (instead of SHC), and the conversion from the gridded data to the geopotential SHC is firstly performed as follows [56]:
{ Δ C ¯ nm Δ S ¯ nm } = 3 ρ w 4 π a ρ ave 1 + k n 2 n + 1 θ = 0 π λ = 0 2 π Δ TWS ( θ , λ ) P ¯ nm ( cos θ ) { cos m λ sin m λ } sin θ   d θ d λ
In this study, the maximum harmonic degree (Nmax) is approximated by N max = 180 ° / ψ 1 / 2 [57], where ψ 1 / 2 is a half-wavelength of a field (in degrees), i.e., a sampling distance of model-simulated ΔTWS. The atmospheric and oceanic signals from the atmospheric and oceanic de-aliasing (AOD) product (see Section 2.2) are then added to the vertical displacement estimate to obtain signals consistent with the GPS measurements.

3.2. GRACE Data Assimilation

GRACE-derived ΔTWS is assimilated into PCR-GLOBWB using the ensemble Kalman smoother (EnKS) [30]. Comprehensive details on the DA configurations are given in Appendix A, while this section describes the core processing concepts. The EnKS comprises two-step runs, a forward run, and a repeat run. In GRACE DA, a forward run is performed to obtain the monthly average TWS states (consistent with a GRACE or GRACE-FO time window), and a repeat run is used to distribute the update evenly throughout the month. A repeat run acts like a smoother to smooth the simulation’s trajectory and reduce the time series’ spurious discontinuity. The DA process begins with a model propagation () of each ensemble i (at time step t) as follows:
x t i = ( x t 1 i , f t i , u i )
where x is a vector containing the state variables, f is a vector of meteorological forcing data, and u represents the model parameters. The ensemble generation process is described in Appendix A. The number of ensembles used in this study is 100. A forward run is performed for approximately one month (T) to obtain the monthly TWS state variables. When GRACE observation is available, the model states are updated as follows:
X T + i = X T i + K T [ y T i H T X T i ]
where X contains 27 monthly average TWS variables (e.g., soil moisture, groundwater, and snow; see [30] for more details); y is an observation vector; H is a measurement operator relating the state to the observation vector, and K is a Kalman gain representing the weighting factor that determines the contribution of observation in the state variable’s adjustment. The subscript signs (−,+) indicate the update status: before (−) and after (+) the update, respectively. The Kalman gain is expressed as
K = C X H ( C H + C v ) 1
where C H is the covariance matrix of model TWS predictions ( H T X T i ), C v is the error covariance matrix of the GRACE observation, and C X H is the cross-covariance matrix computed between model TWS states and model TWS predictions. The covariance matrices can be computed empirically from the ensemble members [29]. The daily TWS update is then computed by dividing the monthly update values (from Equation (6)) by the number of days in that month. Next, a repeat run is performed by repeating a forward run (Equation (5)) but adding the daily TWS update to the initial states (associated with TWS variables) at the beginning of each day. The repeated run is skipped when GRACE or GRACE-FO observation is absent.

4. Results

4.1. The Estimations of TWS and Vertical Displacement Variations

The methods described in Section 3.1 are used to compute the monthly ΔTWS and Δr associated with the GRACE data, PCR-GLOBWB model, and GRACE DA estimate between April 2002 and December 2018. The estimated time series are least-squares-fitted and the ΔTWS annual amplitude estimates are shown in Figure 3. Note that only available GRACE/GRACE-FO data are used in the amplitude computation. The PCR-GLOBWB ΔTWS estimate (Figure 3b) exhibits much greater spatial details than the GRACE estimate (Figure 3a). Strong TWS variations are seen in river and lake areas, implying a significant contribution of surface water to the hydrological dynamics in CSEA. Although a direct comparison between GRACE and PCR-GLOBWB estimates is difficult due to the inconsistency in their spatial resolution, some distinct differences between them can still be seen, particularly in Central and Eastern Thailand, where PCR-GLOBWB has smaller annual amplitudes by approximately a factor of two (see Figure 3a vs. Figure 3b). The model errors associated with (but not limited to), e.g., the inadequate representation of physical processes, uncalibrated model parameters, or inaccurate forcing data, likely lead to underestimated ΔTWS [27]. After assimilating GRACE into PCR-GLOBWB, the ΔTWS annual amplitude estimate is closer to the GRACE result (Figure 3c). The GRACE DA’s application also preserves the high-spatial details of the ΔTWS (from PCR-GLOBWB), and the ΔTWS features over lakes/reservoirs are still clearly seen.
The GRACE Δr estimate’s annual amplitude exhibits a similar spatial pattern to the GRACE ΔTWS result, showing more significant variations in the northwestern and southeastern parts of CSEA (Figure 3d). Relative to the ΔTWS result (Figure 3a), the Δr spatial pattern is slightly shifted toward Cambodia’s coastal areas (see, e.g., the red oval in Figure 3d), attributed to significant atmospheric and oceanic loading effects in these semi-enclosed bay regions (see later discussions). The PCR-GLOBWB Δr estimate agrees with the GRACE result in, e.g., Northern Thailand and Laos (Figure 3e). However, the model underestimates Δr in Central Thailand and Cambodia, consistent with underestimated ΔTWS seen in Figure 3b. Significant annual Δr are also seen in the reservoir and lake areas (Figure 3e), confirming the large surface water seasonality in CSEA. In line with the ΔTWS estimate, applying GRACE DA leads to an improved agreement between the Δr and GRACE estimates and maintains the relatively high spatial resolution of the PCR-GLOBWB model (Figure 3f).
Land deformation is governed by multiple mass loading sources, i.e., land, atmosphere, and ocean. While the land component dominates the Δr in CSEA, atmospheric and oceanic loadings are found to play a crucial role in certain parts of the region, seen from their notable contributions to total vertical displacement (Figure 4). The signal contribution associated with different components is computed as the Δ r C / Δ r , where Δ r C is the annual amplitude of vertical displacement induced by an individual mass loading component (i.e., land and ocean/atmosphere), and Δ r is the amplitude of the total vertical displacement estimate (i.e., the integration of all components). Figure 4 suggests that the land water storage variation is the primary driver of Δr (see Figure 4a vs. Figure 4b), while the contribution of atmospheric and oceanic loadings is relatively small (~20%) in the inland areas (Figure 4c). However, large atmospheric and oceanic variations are found near Thailand and Cambodia coasts (see, e.g., a red oval in Figure 4c), influenced by a strong seasonality of ocean bottom pressure variations in the Gulf of Thailand [58], which contributes up to 60% to the Δr annual amplitude in the coastal regions.

4.2. The Impact of Flooding on Vertical Displacement Estimates

CSEA has experienced several flood events induced by extreme monsoon rainfalls. Figure 5 illustrates the MODIS-derived NDWI during the flood event in 2011 (Figure 5a; the most severe flood since 2002) compared with the average NDWI value computed between 2002 and 2018 (Figure 5b). The positive NDWI values reveal approximate surface water extents in the region, in which the lake, reservoir, and flooding features are clearly seen. During the flood event in 2011 (July–December 2011; Figure 5a), the surface water areas of CSEA were increased by a factor of 2.5 relative to the regionally average value (Figure 5b). The flooding severities were apparent in the Chao Phraya River Basin, Tonlé Sap Basin, and most of Southeastern Cambodia, where several major cities were inundated for more than one month [38].
The basin-averaged ΔTWS and Δr time series over the Chao Phraya River and Tonlé Sap Basins are analyzed together with the precipitation and flood extent (Figure 6). The PCR-GLOBWB ΔTWS and Δr estimates depict a significant agreement with the GRACE estimates (Figure 6a–d) in terms of seasonal variation (with ρ > 0.85; see Table 2). The discrepancy between PCR-GLOBWB and GRACE is seen in the annual amplitude estimates, particularly in the Tonlé Sap Basin, where the model underestimates the ΔTWS and Δr by >10% (Figure 6b; see also Table 2). After applying the GRACE DA, both the ΔTWS and Δr estimates move toward the GRACE observations, and the ρ and RMSD estimates (with respect to GRACE) are improved by up to 9% and 44%, respectively (see Table 2). Figure 6 also illustrates that the ΔTWS and Δr seasonal variations are governed mainly by the precipitation seasonality (Figure 6e,f). Exceptionally low or high ΔTWSs correspond to extreme precipitation and flooding (see, e.g., ΔTWS in 2011; Figure 6a,e,g), which results in noticeable land subsidence (negative Δr) as a response to the increased water mass.
The correlation analysis shown in Figure 7 reveals that the flood extent has a stronger connection with the ΔTWS than precipitation. The precipitation–flooding relationship is affected by several factors, e.g., climate conditions, land cover (urbanization), and catchment characteristics [59], which likely lessen the correlation between the two. On the other hand, the ΔTWS reflects a saturated soil and groundwater storage condition, which is an essential factor controlling flood occurrence. When soil or groundwater storage is saturated (e.g., a high ΔTWS), the throughfall cannot infiltrate through the soil but exits the system as surface runoff or a flood. This land process mechanism also leads to a phase delay of flood occurrence (relative to, e.g., the ΔTWS), where the maximum flood extent is observed after the ΔTWS peak by about one month (see, e.g., the Chao Phraya River Basin in 2011; Figure 6g). It is noteworthy that the Δr estimates’ anticorrelations reflect the nature of vertical displacement responses to water mass loads, which results in about one cycle (π) phase difference between Δr and ΔTWS (and, also, the flood extent) (see, e.g., Figure 6a vs. Figure 6c).

4.3. Validation with GPS Measurements

The Δr estimates are evaluated against the GPS measurements at the S1–S3 sites (Figure 8). The long-term mean of vertical displacement estimates (GRACE, PCR-GLOBWB, and GRACE DA) with respect to the GPS data period is removed from the time series prior to the comparison. The GRACE, PCR-GLOBWB, and GRACE DA Δr estimates depict consistent seasonal variability with the GPS time series in all three evaluated locations. Discrepancies are seen during extreme events (e.g., the flood event in 2011), where GPS measurements exhibit higher variability than that observed from the GRACE and model estimates (see, e.g., the blue ovals in Figure 8b). The higher Δr variations are likely associated with local signals that coarse spatial resolution satellite observations (~300 km) or model estimates (~10 km) cannot effectively capture [6,60]. At S2 and S3 (see Figure 8c,e), the excessive land subsidence (negative Δr) during the flood event in 2011 is clearly seen. The subsidence magnitude in July–December 2011 is increased by up to 90% compared with the average subsidence computed from all July–December Δr estimates (see, e.g., Figure 8d,f). The S1 location is located beyond the flood zone (see the S1 location in, e.g., Figure 5b), and no noticeable land subsidence feature is detected in 2011 (Figure 8a).
The notable agreement between Δr estimates and GPS measurements is also seen in the annual amplitude and phase estimates (Table 3). The annual amplitude estimates from GRACE, PCR-GLOBWB, and GRACE-DA agree with the GPS measurements within ~1 mm, which is reasonable, considering the 1-cm accuracy of the GPS data in the up component [16]. Similarly, the phase estimates agree within ~0.3 months, consistent with the evaluated time interval (one month). The GRACE results show the highest agreement with the GPS estimates, providing a ~17% closer annual amplitude than the PCR-GLOBWB results (e.g., 9.11 mm vs. 7.64 mm with respect to 8.67 mm; see Table 3). As seen in Section 4.2, the GRACE DA updates the ΔTWS estimates toward the GRACE observations, resulting in an improved agreement with the GPS data compared with the model estimates alone (e.g., the GRACE DA increases the agreement in the annual amplitude by 8%).
Validating against the GPS measurements, the GRACE DA Δr estimates provide ρ and RMSD values of approximately 0.7 and 5 mm, respectively (Table 4). The RMSD is smaller than the GRACE errors by a factor of two to three [61], denoting the effectiveness of GRACE data usage in data-sparse regions. Like previous analyses, the Δr estimates associated with GRACE observations (GRACE and GRACE DA) outperform the model estimate, i.e., by providing higher ρ and RMSD values by ~3%. The model, however, still presents one advantage in providing a continuous time series, in which Δr can still be estimated despite the absence of GRACE, GRACE-FO, or GPS data (see, e.g., Figure 8a).
Our final analysis investigates the impact of atmospheric and oceanic loadings on the vertical displacement accuracy. The assessment is performed by comparing two Δr estimate scenarios, including and excluding atmospheric and oceanic loadings in the Δr computation. Table 4 illustrates that excluding atmospheric and oceanic signals from Δr estimates worsens the ρ, RMSD, and RMS reduction estimates by ~10%. This analysis reveals that, although land hydrology explains about 90% of the Δr accuracy, the result is likely suboptimal when considering the land component alone. The inclusion of atmospheric and oceanic signals should be considered to achieve an accurate Δr estimate.

5. Discussion

The performances of GRACE, PCR-GLOBWB, and GRACE DA-derived Δr in data-sparse regions (here is CSEA) are assessed in this study. GRACE provides the most accurate Δr estimates, closest to the GPS measurements, with ρ and RMSD values of > 0.8 and ~5 mm, respectively. A similar agreement is seen in most monsoon regions [7,11,13,33]. GRACE fails to capture local signals, due likely to the gravity data’s limited spatial resolution (e.g., ~300 km), causing a noticeable mismatch between the GRACE and GPS time series [60]. However, the mismatch is negligible compared with GPS uncertainties [16] and does not significantly affect the GRACE’s overall performance in CSEA. A notable agreement in annual amplitude and phase estimates is also found. Previous studies (e.g., [11,33,62]) discussed the difficulty of aligning the amplitude or phase between GRACE and GPS data, which was caused by ineffective implementations of earth models. Such an issue is not observed here, likely due to the use of the PREM model, which was found to be the optimal earth model in the Southeast Asia region [33].
The PCR-GLOBWB model can also deliver an accurate Δr estimate, despite its inferior performance to GRACE. The model underestimates the Δr annual amplitudes but provides greater spatial details than the GRACE estimates. Underperformed Δr is attributed in part to incorrect model ΔTWS estimates associated with, e.g., imperfect model physics, ineffective parameter calibration, or inaccurate meteorological forcing input [27]. Assimilating GRACE into PCR-GLOBWB apparently improves Δr estimates (compared to PCR-GLOBWB results). It is noticed that GRACE DA does not outperform GRACE observations in terms of achieving a high Δr accuracy. This behavior can be explained by the GRACE DA’s nature that optimizes the contributions from the model and observation [29], resulting in the Δr accuracy falling between the model and GRACE results. However, GRACE DA estimates present one clear advantage (over that of GRACE observations) in providing high spatiotemporal resolution Δr estimates (e.g., ~10 km daily). In particular, accurate vertical deformation can always be obtained from GRACE DA, despite the GRACE observation absence. The GRACE DA’s effectiveness found from this study also supports the concept of Springer et al. [26] in utilizing DA estimates to isolate high-frequency hydrologic signals from GPS measurements.
The GPS, GRACE, PCR-GLOBWB, and GRACE DA results depicted elevated land subsidence during flood events, particularly in 2011, where the subsidence intensity was increased by up to 90%. During flooding, saturated soil and groundwater conditions appear to cause large-scale land subsidence that can be detected from ground, model, and satellite data. Additionally, we found that the CSEA flood extent had a much stronger relationship to the ΔTWS than precipitation. This finding is in line with Tangdamrongsub et al. [36], who utilized such a robust connection for predicting the ΔTWS (from the MODIS-derived surface water extent) in the Tonlé Sap Basin. This study’s analysis may underline the potential usage of ΔTWS estimates in large-scale flood monitoring applications.
It is important to include atmospheric and oceanic loadings in the Δr computation. Omitting such components leads to an accuracy reduction of about 10%. The impact of neglecting the ocean loading component may be intensified in the semi-enclosed bay regions, e.g., the Gulf of Thailand, where the ocean bottom pressure effect is significant [58,63]. Note that this study obtains the atmosphere and ocean signals from the AOD product consistent with GRACE Level-2 processing while using different de-aliasing products or geophysical models, e.g., [64], which may have a different impact on the Δr estimate. The performance comparison is beyond the scope of this study but can be considered in future analyses.

6. Conclusions

This study assessed large-scale land deformation in CSEA using GRACE/GRACE-FO, PCR-GLOBWB, GRACE/GRACE FO DA, and GPS data. It was found that the Δr estimate in CSEA is primarily governed by hydrologic variations and extreme events (i.e., monsoon flooding), which can be detected by the GRACE/GRAE-FO satellites, hydrology model, and ground measurements. The flooding analysis associated with MODIS data illustrated a significant impact of extreme weather conditions on the Δr estimates, where a degree of land subsidence can be intensified by almost 100% during flooding. Evaluating against the GPS data, the GRACE/GRACE-FO, PCR-GLOBWB, and GRACE/GRACE-FO DA effectively observed regional Δr annual variations (e.g., spatial resolution > 10 km), while only the GPS measurements could capture local Δr signals. The evaluation also exhibited the advantage of GRACE/GRACE-FO DA in improving Δr estimates and maintaining great spatiotemporal details. A similar GRACE DA performance can likely be expected in other data-sparse regions, particularly ones with similar monsoon climates like CSEA. Future experiments in different regions can be performed to confirm the GRACE DA’s effectiveness in Δr computations.

Author Contributions

N.T. designed and performed the experiments, analyzed the data, and wrote the paper. M.Š. edited the paper. All authors have read and agreed to the published version of the manuscript.

Funding

This paper received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

This study was performed based on public-access data. No new data were created in this study. Data sharing is not applicable to this article.

Acknowledgments

NT was supported by the NASA Earth Science Division in support of the National Climate Assessment. MŠ was supported by the project LO1506 of the Czech Ministry of Education, Youth and Sports.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. The Configuration of GRACE Data Assimilation

In the GRACE DA process, the ensemble sets of forcing data, model parameters, and observations are first generated by adding random noise to the nominal value. For forcing data, the precipitation is perturbed with lognormal noise with 30% of the nominal values range, while the temperature is perturbed with 2 °C Gaussian noise. For the model parameters, 14 TWS-related parameters are perturbed using Gaussian noise with 20% of the nominal values range. The parameters include the saturated conductivity of the upper and lower soil layers, the storage capacity of the upper and lower soil layers, the minimum and maximum soil depth fractions, the groundwater recession coefficient, and the minimum crop coefficient (see [30] for more details). The error magnitudes are selected through trial and error, mainly to allow ensemble growth between two consecutive updates in the DA process. For the GRACE observations, the temporal mean value of the TWS (April 2002–December 2018) from the PCR-GLOBWB estimate is first added to the GRACE-derived ΔTWS to obtain the absolute TWS prior to the assimilation process. This process is to alleviate the systematic mean difference between the model TWS estimate and the GRACE data. The GRACE data is perturbed with 2-cm (and a 3° spatial correlation error) Gaussian noise [61].

References

  1. Tregoning, P.; Watson, C.; Ramillien, G.; McQueen, H.; Zhang, J. Detecting hydrologic deformation using GRACE and GPS. Geophys. Res. Lett. 2009, 36. [Google Scholar] [CrossRef] [Green Version]
  2. Chaussard, E.; Amelung, F.; Abidin, H.; Hong, S. Sinking cities in Indonesia: ALOS PALSAR detects rapid subsidence due to groundwater and gas extraction. Remote Sens. Environ. 2013, 128, 150–161. [Google Scholar] [CrossRef]
  3. Phien-Wej, N.; Giao, P.; Nutalaya, P. Land subsidence in Bangkok, Thailand. Eng. Geol. 2006, 82, 187–201. [Google Scholar] [CrossRef]
  4. Schmidt, D.A.; Burgmann, R. Time-dependent land uplift and subsidence in the Santa Clara valley, California, from a large interferometric synthetic aperture radar data set. J. Geophys. Res. Space Phys. 2003, 108. [Google Scholar] [CrossRef] [Green Version]
  5. Dill, R.; Dobslaw, H. Numerical simulations of global-scale high-resolution hydrological crustal deformations. J. Geophys. Res. Solid Earth 2013, 118, 5008–5017. [Google Scholar] [CrossRef]
  6. Fu, Y.; Freymueller, J.T.; Jensen, T. Seasonal hydrological loading in southern Alaska observed by GPS and GRACE. Geophys. Res. Lett. 2012, 39. [Google Scholar] [CrossRef] [Green Version]
  7. Wang, L.; Chen, C.; Du, J.; Wang, T. Detecting seasonal and long-term vertical displacement in the North China Plain using GRACE and GPS. Hydrol. Earth Syst. Sci. 2017, 21, 2905–2922. [Google Scholar] [CrossRef] [Green Version]
  8. Heki, K. Seasonal Modulation of Interseismic Strain Buildup in Northeastern Japan Driven by Snow Loads. Science 2001, 293, 89–92. [Google Scholar] [CrossRef] [PubMed]
  9. Fu, Y.; Argus, D.F.; Landerer, F.W. GPS as an independent measurement to estimate terrestrial water storage variations in Washington and Oregon. J. Geophys. Res. Solid Earth 2015, 120, 552–566. [Google Scholar] [CrossRef]
  10. Argus, D.F.; Fu, Y.; Landerer, F.W. Seasonal variation in total water storage in California inferred from GPS observations of vertical land motion. Geophys. Res. Lett. 2014, 41, 1971–1980. [Google Scholar] [CrossRef]
  11. Nahmani, S.; Bock, O.; Bouin, M.-N.; Santamaría-Gómez, A.; Boy, J.-P.; Collilieux, X.; Métivier, L.; Panet, I.; Genthon, P.; De Linage, C.; et al. Hydrological deformation induced by the West African Monsoon: Comparison of GPS, GRACE and loading models. J. Geophys. Res. Space Phys. 2012, 117. [Google Scholar] [CrossRef] [Green Version]
  12. Fu, Y.; Argus, D.F.; Freymueller, J.T.; Heflin, M.B. Horizontal motion in elastic response to seasonal loading of rain water in the Amazon Basin and monsoon water in Southeast Asia observed by GPS and inferred from GRACE. Geophys. Res. Lett. 2013, 40, 6048–6053. [Google Scholar] [CrossRef]
  13. Steckler, M.; Nooner, S.L.; Akhter, S.H.; Chowdhury, S.K.; Bettadpur, S.; Seeber, L.; Kogan, M.G. Modeling Earth deformation from monsoonal flooding in Bangladesh using hydrographic, GPS, and Gravity Recovery and Climate Experiment (GRACE) data. J. Geophys. Res. Space Phys. 2010, 115. [Google Scholar] [CrossRef]
  14. Klein, E.; Duputel, Z.; Zigone, D.; Vigny, C.; Boy, J.-P.; Doubre, C.; Meneses, G. Deep Transient Slow Slip Detected by Survey GPS in the Region of Atacama, Chile. Geophys. Res. Lett. 2018, 45. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Tregoning, P.; Ramillien, G.; McQueen, H.; Zwartz, D. Glacial isostatic adjustment and nonstationary signals observed by GRACE. J. Geophys. Res. Space Phys. 2009, 114, 114. [Google Scholar] [CrossRef]
  16. Kouba, J. A Guide to Using International GPS Service (IGS) Products; IGS Central Bureau: Pasadena, CA, USA, 2003. [Google Scholar]
  17. Saji, A.P.; Sunil, P.; Sreejith, K.M.; Gautam, P.K.; Kumar, K.V.; Ponraj, M.; Amirtharaj, S.; Shaju, R.M.; Begum, S.K.; Reddy, C.D.; et al. Surface Deformation and Influence of Hydrological Mass Over Himalaya and North India Revealed from a Decade of Continuous GPS and GRACE Observations. J. Geophys. Res. Earth Surf. 2020, 125. [Google Scholar] [CrossRef]
  18. Tapley, B.D.; Bettadpur, S.V.; Ries, J.; Thompson, P.F.; Watkins, M.M. GRACE Measurements of Mass Variability in the Earth System. Science 2004, 305, 503–505. [Google Scholar] [CrossRef] [Green Version]
  19. Flechtner, F.; Neumayer, K.-H.; Dahle, C.; Dobslaw, H.; Fagiolini, E.; Raimondo, J.-C.; Güntner, A. What Can be Expected from the GRACE-FO Laser Ranging Interferometer for Earth Science Applications? In Remote Sensing and Water Resources; Cazenave, A., Champollion, N., Benveniste, J., Chen, J., Eds.; Springer: Cham, Switzerland, 2016; pp. 263–280. ISBN 978-3-319-32449-4. [Google Scholar]
  20. Tapley, B.D.; Watkins, M.M.; Flechtner, F.M.; Reigber, C.; Bettadpur, S.; Rodell, M.; Sasgen, I.; Famiglietti, J.S.; Landerer, F.W.; Chambers, D.P.; et al. Contributions of GRACE to understanding climate change. Nat. Clim. Chang. 2019, 9, 358–369. [Google Scholar] [CrossRef]
  21. Scanlon, B.R.; Zhang, Z.; Save, H.; Wiese, D.N.; Landerer, F.W.; Long, D.; Longuevergne, L.; Chen, J. Global evaluation of new GRACE mascon products for hydrologic applications. Water Resour. Res. 2016, 52, 9412–9429. [Google Scholar] [CrossRef]
  22. Ahmed, M. Sustainable management scenarios for northern Africa’s fossil aquifer systems. J. Hydrol. 2020, 589, 125196. [Google Scholar] [CrossRef]
  23. Swenson, S.C.; Lawrence, D.M. Assessing a dry surface layer-based soil resistance parameterization for the Community Land Model using GRACE and FLUXNET-MTE data. J. Geophys. Res. Atmos. 2014, 119. [Google Scholar] [CrossRef]
  24. King, M.; Moore, P.; Clarke, P.; Lavallée, D. Choice of optimal averaging radii for temporal GRACE gravity solutions, a comparison with GPS and satellite altimetry. Geophys. J. Int. 2006, 166, 1–11. [Google Scholar] [CrossRef] [Green Version]
  25. Ferreira, V.G.; Montecino, H.D.; Ndehedehe, C.E.; Del Rio, R.A.; Cuevas, A.; De Freitas, S.R.C. Determining seasonal displacements of Earth’s crust in South America using observations from space-borne geodetic sensors and surface-loading models. Earth Planets Space 2019, 71, 84. [Google Scholar] [CrossRef] [Green Version]
  26. Springer, A.; Karegar, M.A.; Kusche, J.; Keune, J.; Kurtz, W.; Kollet, S. Evidence of daily hydrological loading in GPS time series over Europe. J. Geod. 2019, 93, 2145–2153. [Google Scholar] [CrossRef]
  27. Tangdamrongsub, N.; Steele-Dunne, S.C.; Gunter, B.C.; Ditmar, P.G.; Weerts, A. Data assimilation of GRACE terrestrial water storage estimates into a regional hydrological model of the Rhine River basin. Hydrol. Earth Syst. Sci. 2015, 19, 2079–2100. [Google Scholar] [CrossRef] [Green Version]
  28. Scanlon, B.R.; Zhang, Z.; Save, H.; Sun, A.Y.; Schmied, H.M.; Van Beek, L.P.H.; Wiese, D.N.; Wada, Y.; Long, D.; Reedy, R.C.; et al. Global models underestimate large decadal declining and rising water storage trends relative to GRACE satellite data. Proc. Natl. Acad. Sci. USA 2018, 115, E1080–E1089. [Google Scholar] [CrossRef] [Green Version]
  29. Evensen, G. The Ensemble Kalman Filter: Theoretical formulation and practical implementation. Ocean Dyn. 2003, 53, 343–367. [Google Scholar] [CrossRef]
  30. Tangdamrongsub, N.; Steele-Dunne, S.C.; Gunter, B.C.; Ditmar, P.G.; Sutanudjaja, E.H.; Sun, Y.; Xia, T.; Wang, Z. Improving estimates of water resources in a semi-arid region by assimilating GRACE data into the PCR-GLOBWB hydrological model. Hydrol. Earth Syst. Sci. 2017, 21, 2053–2074. [Google Scholar] [CrossRef] [Green Version]
  31. Li, B.; Rodell, M.; Kumar, S.V.; Beaudoing, H.; Getirana, A.; Zaitchik, B.F.; De Goncalves, L.G.; Cossetin, C.; Bhanja, S.N.; Mukherjee, A.; et al. Global GRACE Data Assimilation for Groundwater and Drought Monitoring: Advances and Challenges. Water Resour. Res. 2019, 55, 7564–7586. [Google Scholar] [CrossRef] [Green Version]
  32. Yin, W.; Han, S.; Zheng, W.; Yeo, I.-Y.; Hu, L.; Tangdamrongsub, N.; Ghobadi-Far, K. Improved water storage estimates within the North China Plain by assimilating GRACE data into the CABLE model. J. Hydrol. 2020, 590, 125348. [Google Scholar] [CrossRef]
  33. Chanard, K.; Avouac, J.P.; Ramillien, G.; Genrich, J. Modeling deformation induced by seasonal variations of continental water in the Himalaya region: Sensitivity to Earth elastic structure. J. Geophys. Res. Solid Earth 2014, 119, 5097–5113. [Google Scholar] [CrossRef] [Green Version]
  34. Sutanudjaja, E.H.; Van Beek, L.P.H.; Wanders, N.; Wada, Y.; Bosmans, J.H.C.; Drost, N.; Van Der Ent, R.J.; De Graaf, I.E.M.; Hoch, J.M.; De Jong, K.; et al. PCR-GLOBWB 2: A 5 arcmin global hydrological and water resources model. Geosci. Model Dev. 2018, 11, 2429–2453. [Google Scholar] [CrossRef] [Green Version]
  35. Van Dam, T.; Wahr, J.M.; Lavallée, D. A comparison of annual vertical crustal displacements from GPS and Gravity Recovery and Climate Experiment (GRACE) over Europe. J. Geophys. Res. Space Phys. 2007, 112. [Google Scholar] [CrossRef] [Green Version]
  36. Tangdamrongsub, N.; Ditmar, P.; Steele-Dunne, S.; Gunter, B.; Sutanudjaja, E. Assessing total water storage and identifying flood events over Tonlé Sap basin in Cambodia using GRACE and MODIS satellite observations combined with hydrological models. Remote Sens. Environ. 2016, 181, 162–173. [Google Scholar] [CrossRef] [Green Version]
  37. Jamrussri, S.; Toda, Y. Simulating past severe flood events to evaluate the effectiveness of nonstructural flood countermeasures in the upper Chao Phraya River Basin, Thailand. J. Hydrol. Reg. Stud. 2017, 10, 82–94. [Google Scholar] [CrossRef]
  38. World Bank. Thai Flood 2011: Rapid Assessment for Resilient Recovery and Reconstruction Planning; World Bank: Washington, DC, USA, 2012. [Google Scholar]
  39. Swenson, S.; Chambers, D.P.; Wahr, J. Estimating geocenter variations from a combination of GRACE and ocean model output. J. Geophys. Res. Sol. Earth 2008, 113, 08410. [Google Scholar] [CrossRef] [Green Version]
  40. Loomis, B.D.; Rachlin, K.E.; Wiese, D.N.; Landerer, F.W.; Luthcke, S.B. Replacing GRACE/GRACE-FO With Satellite Laser Ranging: Impacts on Antarctic Ice Sheet Mass Change. Geophys. Res. Lett. 2020, 47. [Google Scholar] [CrossRef]
  41. Swenson, S.; Wahr, J. Post-processing removal of correlated errors in GRACE data. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef]
  42. Jekeli, C. Alternative Methods to Smooth the Earth’s Gravity Field; Scientific Report 327; The Ohio State University: Columbus, OH, USA, 1981. [Google Scholar]
  43. Caron, L.; Ivins, E.R.; Larour, E.; Adhikari, S.; Nilsson, J.; Blewitt, G. GIA Model Statistics for GRACE Hydrology, Cryosphere, and Ocean Science. Geophys. Res. Lett. 2018, 45, 2203–2212. [Google Scholar] [CrossRef]
  44. Zhou, H.; Luo, Z.; Tangdamrongsub, N.; Zhou, H.; He, L.; Xu, C.; Li, Q.; Wu, Y. Identifying Flood Events over the Poyang Lake Basin Using Multiple Satellite Remote Sensing Observations, Hydrological Models and In Situ Data. Remote Sens. 2018, 10, 713. [Google Scholar] [CrossRef] [Green Version]
  45. Dobslaw, H.; Bergmann-Wolf, I.; Dill, R.; Poropat, L.; Thomas, M.; Dahle, C.; Esselborn, S.; König, R.; Flechtner, F. A new high-resolution model of non-tidal atmosphere and ocean mass variability for de-aliasing of satellite gravity observations: AOD1B RL06. Geophys. J. Int. 2017, 211, 263–269. [Google Scholar] [CrossRef] [Green Version]
  46. Van Beek, L.P.H.; Wada, Y.; Bierkens, M.F.P. Global monthly water stress: 1. Water balance and water availability. Water Resour. Res. 2011, 47. [Google Scholar] [CrossRef]
  47. Sutanudjaja, E.H.; Van Beek, L.P.H.; De Jong, S.M.; Van Geer, F.C.; Bierkens, M.F.P. Large-scale groundwater modeling using global datasets: A test case for the Rhine-Meuse basin. Hydrol. Earth Syst. Sci. 2011, 15, 2913–2935. [Google Scholar] [CrossRef] [Green Version]
  48. Sutanudjaja, E.H.; Van Beek, L.P.H.; De Jong, S.M.; Van Geer, F.C.; Bierkens, M.F.P. Calibrating a large-extent high-resolution coupled groundwater-land surface model using soil moisture and discharge data. Water Resour. Res. 2014, 50, 687–705. [Google Scholar] [CrossRef]
  49. Huffman, G.J.; Bolvin, D.T.; Braithwaite, D.; Hsu, K.; Joyce, R.; Kidd, C.; Nelkin, E.J.; Sorooshian, S.; Tan, J.; Xie, P. NASA Global Precipitation Measurement (GPM) Integrated Multi-SatellitE Retrievals for GPM (IMERG); Algorithm Theoretical Basis Document (ATBD) Version 06; NASAGSFC: Greenbelt, MD, USA, 2019.
  50. Lu, J.; Sun, G.; McNulty, S.G.; Amatya, D.M. A Comparison of Six Potential Evapotranspiration Methods for Regional Use in the Southeastern United States. JAWRA J. Am. Water Resour. Assoc. 2005, 41, 621–633. [Google Scholar] [CrossRef]
  51. Blewitt, G.; Hammond, W.C.; Kreemer, C. Harnessing the GPS Data Explosion for Interdisciplinary Science. Available online: https://eos.org/science-updates/harnessing-the-gps-data-explosion-for-interdisciplinary-science (accessed on 30 September 2020).
  52. Roger, P.J.C.; Vermote, E.F.; Ray, J.P. MODIS Surface Reflectance User’s Guide; Collection 6; NASAGSFC: Greenbelt, MD, USA, 2015.
  53. McFeeters, S.K. The use of the Normalized Difference Water Index (NDWI) in the delineation of open water features. Int. J. Remote Sens. 1996, 17, 1425–1432. [Google Scholar] [CrossRef]
  54. Tangdamrongsub, N.; Han, S.-C.; Jasinski, M.F.; Šprlák, M. Quantifying water storage change and land subsidence induced by reservoir impoundment using GRACE, Landsat, and GPS data. Remote Sens. Environ. 2019, 233, 111385. [Google Scholar] [CrossRef]
  55. Wang, H.; Xiang, L.; Jia, L.; Jiang, L.; Wang, Z.; Hu, B.; Gao, P. Load Love numbers and Green’s functions for elastic Earth models PREM, iasp91, ak135, and modified models with refined crustal structure from Crust 2.0. Comput. Geosci. 2012, 49, 190–199. [Google Scholar] [CrossRef]
  56. Wahr, J.; Molenaar, M.; Bryan, F. Time variability of the Earth’s gravity field: Hydrological and oceanic effects and their possible detection using GRACE. J. Geophys. Res. Space Phys. 1998, 103, 30205–30229. [Google Scholar] [CrossRef]
  57. Vishwakarma, B.D.; Devaraju, B.; Sneeuw, N. What Is the Spatial Resolution of grace Satellite Products for Hydrology? Remote Sens. 2018, 10, 852. [Google Scholar] [CrossRef] [Green Version]
  58. Wouters, B.; Chambers, D.P. Analysis of seasonal ocean bottom pressure variability in the Gulf of Thailand from GRACE. Glob. Planet. Chang. 2010, 74, 76–81. [Google Scholar] [CrossRef]
  59. Sharma, A.; Wasko, C.; Lettenmaier, D.P. If Precipitation Extremes Are Increasing, Why Aren’t Floods? Water Resour. Res. 2018, 54, 8545–8551. [Google Scholar] [CrossRef]
  60. Ferreira, V.G.; Liu, Z.; Montecino, H.C.; Yuan, P.; Kelly, C.I.; Mohammed, A.S.; Han, L.Y. Reciprocal comparison of geodetically sensed and modeled vertical hydrological loading products. Acta Geod. Geophys. 2019, 55, 23–49. [Google Scholar] [CrossRef]
  61. Wahr, J.M.; Swenson, S.; Velicogna, I. Accuracy of GRACE mass estimates. Geophys. Res. Lett. 2006, 33. [Google Scholar] [CrossRef] [Green Version]
  62. Bettinelli, P.; Avouac, J.; Flouzat, M.; Bollinger, L.; Ramillien, G.; Rajaure, S.; Sapkota, S. Seasonal variations of seismicity and geodetic strain in the Himalaya induced by surface hydrology. Earth Planet. Sci. Lett. 2008, 266, 332–344. [Google Scholar] [CrossRef]
  63. Van Dam, T.; Collilieux, X.; Wuite, J.; Altamimi, Z.; Ray, J.R. Nontidal ocean loading: Amplitudes and potential effects in GPS height time series. J. Geod. 2012, 86, 1043–1057. [Google Scholar] [CrossRef] [Green Version]
  64. Mémin, A.; Boy, J.-P.; Santamaría-Gómez, A. Correcting GPS measurements for non-tidal loading. GPS Solut. 2020, 24, 45. [Google Scholar] [CrossRef]
Figure 1. The Central South-East Asia (CSEA) region encompassing two major river basins, i.e., the Chao Phraya River and Tonlé Sap Basins. The locations of six available Global Positioning System (GPS) sites are shown as green triangles. Characteristics of GPS sites can be found in Table 1.
Figure 1. The Central South-East Asia (CSEA) region encompassing two major river basins, i.e., the Chao Phraya River and Tonlé Sap Basins. The locations of six available Global Positioning System (GPS) sites are shown as green triangles. Characteristics of GPS sites can be found in Table 1.
Remotesensing 13 00235 g001
Figure 2. Data processing diagram used in this study. The given number in each processing step indicates the paper’s section where comprehensive details can be found. GRACE: Gravity Recovery and Climate Experiment and TWS: terrestrial water storage.
Figure 2. Data processing diagram used in this study. The given number in each processing step indicates the paper’s section where comprehensive details can be found. GRACE: Gravity Recovery and Climate Experiment and TWS: terrestrial water storage.
Remotesensing 13 00235 g002
Figure 3. Annual amplitude estimates of the terrestrial water storage variations (ΔTWS) (first row; ac) and vertical displacement variations (Δr) (second row; de) derived from the GRACE (first column; a,d), PCRaster Global Water Balance (PCR-GLOBWB) model (second column; b,e), and GRACE DA (data assimilation) estimates (third column; c,f). The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)). The red oval in (d) illustrates the location where the atmospheric and oceanic loadings are prominent.
Figure 3. Annual amplitude estimates of the terrestrial water storage variations (ΔTWS) (first row; ac) and vertical displacement variations (Δr) (second row; de) derived from the GRACE (first column; a,d), PCRaster Global Water Balance (PCR-GLOBWB) model (second column; b,e), and GRACE DA (data assimilation) estimates (third column; c,f). The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)). The red oval in (d) illustrates the location where the atmospheric and oceanic loadings are prominent.
Remotesensing 13 00235 g003
Figure 4. (a,b) The Δr estimate associated with land hydrology (a) and atmospheric and oceanic loadings (b). (c) The atmospheric and oceanic signals’ contribution (%) to the total Δr estimate. The contribution is computed as a ratio between the annual amplitude of the atmospheric and oceanic loadings and the annual amplitude of the total loadings. The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)). The red oval in (c) illustrates the location where the atmospheric and oceanic signal contributions are significant.
Figure 4. (a,b) The Δr estimate associated with land hydrology (a) and atmospheric and oceanic loadings (b). (c) The atmospheric and oceanic signals’ contribution (%) to the total Δr estimate. The contribution is computed as a ratio between the annual amplitude of the atmospheric and oceanic loadings and the annual amplitude of the total loadings. The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)). The red oval in (c) illustrates the location where the atmospheric and oceanic signal contributions are significant.
Remotesensing 13 00235 g004
Figure 5. The averaged Moderate Resolution Imaging Spectroradiometer (MODIS)-derived normalized difference water index (NDWI) during the monsoon season, July–December. (Left; a) The averaged NDWI during the flooding in 2011, July–December 2011. (Right; b) The averaged NDWI of all July–December data between 2002 and 2018. For visualization, the NDWI < 0.05 are shown in black. The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)).
Figure 5. The averaged Moderate Resolution Imaging Spectroradiometer (MODIS)-derived normalized difference water index (NDWI) during the monsoon season, July–December. (Left; a) The averaged NDWI during the flooding in 2011, July–December 2011. (Right; b) The averaged NDWI of all July–December data between 2002 and 2018. For visualization, the NDWI < 0.05 are shown in black. The evaluation sites (S1–S3) are shown as green triangles (see the site’s name in (a)).
Remotesensing 13 00235 g005
Figure 6. The basin-averaged ΔTWS (first row; a,b), Δr (second row; c,d), Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG)-derived precipitation (third row; e,f), and MODIS-derived surface water extent (fourth row; g,h) in the Chao Phraya River (first column) and Tonlé Sap (second column) Basins between April 2002 and December 2018. The gray background indicates the flood period of 2011.
Figure 6. The basin-averaged ΔTWS (first row; a,b), Δr (second row; c,d), Integrated Multi-satellite Retrievals for Global Precipitation Measurement Mission (IMERG)-derived precipitation (third row; e,f), and MODIS-derived surface water extent (fourth row; g,h) in the Chao Phraya River (first column) and Tonlé Sap (second column) Basins between April 2002 and December 2018. The gray background indicates the flood period of 2011.
Remotesensing 13 00235 g006
Figure 7. Correlation coefficients between the MODIS-derived surface water extent and different datasets; the ΔTWS (GRACE, PCR-GLOBWB, and GRACE DA); the Δr (GRACE, PCR-GLOBWB, and GRACE DA); and precipitation (IMERG) in the Chao Phraya River (left; a) and Tonlé Sap Basins (right; b).
Figure 7. Correlation coefficients between the MODIS-derived surface water extent and different datasets; the ΔTWS (GRACE, PCR-GLOBWB, and GRACE DA); the Δr (GRACE, PCR-GLOBWB, and GRACE DA); and precipitation (IMERG) in the Chao Phraya River (left; a) and Tonlé Sap Basins (right; b).
Remotesensing 13 00235 g007
Figure 8. (Left column) The Δr estimates from the GPS, GRACE, PCR-GLOBWB, and GRACE DA at the three evaluation sites: S1 (a), S2 (c), and S3 (e). Note that the evaluation period (horizontal axis) is different for different evaluated sites regarding GPS measurement availability. In (c), the blue ovals demonstrate the periods where the GPS Δr is more significant than the GRACE or model estimates. (Right column) The average land subsidence (negative Δr) computed during the great flood of 2011, July–December 2011 (orange) estimated from the GPS, GRACE, PCR-GLOBWB (PCR-G), and GRACE DA (DA) at three evaluation sites: S1 (b), S2 (d), and S3 (f). For comparison, the average land subsidence computed from all July–December (2002–2018) values is also shown (blue).
Figure 8. (Left column) The Δr estimates from the GPS, GRACE, PCR-GLOBWB, and GRACE DA at the three evaluation sites: S1 (a), S2 (c), and S3 (e). Note that the evaluation period (horizontal axis) is different for different evaluated sites regarding GPS measurement availability. In (c), the blue ovals demonstrate the periods where the GPS Δr is more significant than the GRACE or model estimates. (Right column) The average land subsidence (negative Δr) computed during the great flood of 2011, July–December 2011 (orange) estimated from the GPS, GRACE, PCR-GLOBWB (PCR-G), and GRACE DA (DA) at three evaluation sites: S1 (b), S2 (d), and S3 (f). For comparison, the average land subsidence computed from all July–December (2002–2018) values is also shown (blue).
Remotesensing 13 00235 g008
Table 1. Characteristics of Global Positioning System (GPS) sites used in this study. The sites with an asterisk (*) are still in operation. The evaluation sites represent the model grid cell location used in the evaluation process. The measurements from GPS sites inside the same PCRaster Global Water Balance (PCR-GLOBWB) model grid cell are averaged. The averaging process results in three (independent) time series of S1–S3 sites (model grid); see, e.g., Figure 3a for their locations on the map.
Table 1. Characteristics of Global Positioning System (GPS) sites used in this study. The sites with an asterisk (*) are still in operation. The evaluation sites represent the model grid cell location used in the evaluation process. The measurements from GPS sites inside the same PCRaster Global Water Balance (PCR-GLOBWB) model grid cell are averaged. The averaging process results in three (independent) time series of S1–S3 sites (model grid); see, e.g., Figure 3a for their locations on the map.
Evaluation SitesGPS SitesData PeriodsLatitude (deg)Longitude (deg)
S1CHGM2010–201218.80398.950
CMUM *2014–201818.76198.932
S2AITB2010–201414.079100.612
ATH22014–201714.082100.613
NIMT2005–201414.043100.714
S3CUSV *2008–201813.736100.534
Table 2. Annual amplitude and statistical values (correlation coefficient (ρ) and root mean square difference (RMSD)) computed from the Gravity Recovery And Climate Experiment (GRACE), PCR-GLOBWB, and GRACE DA (data assimilation)-derived terrestrial water storage variation (ΔTWS) and vertical displacement variation (Δr) time series in the Chao Phraya River and Tonlé Sap Basins. The ρ and RMSD values are computed with respect to the GRACE estimates. The percentage given in the GRACE DA results represents the improvement with respect to the PCR-GLOBWB estimate.
Table 2. Annual amplitude and statistical values (correlation coefficient (ρ) and root mean square difference (RMSD)) computed from the Gravity Recovery And Climate Experiment (GRACE), PCR-GLOBWB, and GRACE DA (data assimilation)-derived terrestrial water storage variation (ΔTWS) and vertical displacement variation (Δr) time series in the Chao Phraya River and Tonlé Sap Basins. The ρ and RMSD values are computed with respect to the GRACE estimates. The percentage given in the GRACE DA results represents the improvement with respect to the PCR-GLOBWB estimate.
ΔTWSΔr
Amp (mm)ρ (-)RMSD (mm)Amp (mm)ρ (-)RMSD (mm)
The Chao Phraya River Basin
GRACE197.35--9.16--
PCR-GLOBWB171.760.9650.937.700.89115.63
GRACE DA187.840.99 (3%)29.76 (42%)8.710.97 (9%)64.41 (44%)
The Tonlé Sap basin
GRACE236.51--8.07--
PCR-GLOBWB139.150.893.996.890.863.33
GRACE DA185.510.93 (4%)2.82 (29%)7.780.93 (8%)2.51 (25%)
Table 3. The annual amplitude and phase estimates from the GPS, GRACE, PCR-GLOBWB, and GRACE DA-derived Δr at three evaluated sites (S1–S3). The average amplitude and phase values are also given in the right column (Avg).
Table 3. The annual amplitude and phase estimates from the GPS, GRACE, PCR-GLOBWB, and GRACE DA-derived Δr at three evaluated sites (S1–S3). The average amplitude and phase values are also given in the right column (Avg).
Amplitude (mm)Phase (month)
S1S2S3AvgS1S2S3Avg
GPS10.267.378.398.674.234.604.574.47
GRACE10.008.638.719.114.154.544.544.41
PCR-GLOBWB9.236.726.967.644.064.194.074.11
GRACE DA9.837.517.568.304.114.224.114.15
Table 4. The ρ, RMSD, and RMS reduction values (with respect to GPS data) estimated from the GRACE, PCR-GLOBWB, and GRACE DA-derived Δr at three evaluated sites (S1–S3). The statistical values associated with the inclusion (top part) and exclusion (bottom part) of atmospheric and oceanic loadings are shown. The average amplitude and phase values are also given in the right column (Avg).
Table 4. The ρ, RMSD, and RMS reduction values (with respect to GPS data) estimated from the GRACE, PCR-GLOBWB, and GRACE DA-derived Δr at three evaluated sites (S1–S3). The statistical values associated with the inclusion (top part) and exclusion (bottom part) of atmospheric and oceanic loadings are shown. The average amplitude and phase values are also given in the right column (Avg).
ρ (-)RMSD (mm)RMS Reduction (%)
S1S2S3AvgS1S2S3AvgS1S2S3Avg
Including atmospheric and oceanic signals
GRACE0.770.780.890.815.704.573.604.6239.3336.7053.4443.16
PCR-GLOBWB0.750.750.840.785.864.774.264.9636.5633.1144.1637.94
GRACE DA0.770.770.860.805.764.644.064.8237.3334.7849.9740.69
Excluding atmospheric and oceanic signals
GRACE0.720.720.850.766.304.984.115.1330.8328.9246.7935.51
PCR-GLOBWB0.720.660.740.716.335.294.935.5229.6222.0832.6328.11
GRACE DA0.740.670.770.736.225.194.675.3632.0224.5634.9130.50
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tangdamrongsub, N.; Šprlák, M. The Assessment of Hydrologic- and Flood-Induced Land Deformation in Data-Sparse Regions Using GRACE/GRACE-FO Data Assimilation. Remote Sens. 2021, 13, 235. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020235

AMA Style

Tangdamrongsub N, Šprlák M. The Assessment of Hydrologic- and Flood-Induced Land Deformation in Data-Sparse Regions Using GRACE/GRACE-FO Data Assimilation. Remote Sensing. 2021; 13(2):235. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020235

Chicago/Turabian Style

Tangdamrongsub, Natthachet, and Michal Šprlák. 2021. "The Assessment of Hydrologic- and Flood-Induced Land Deformation in Data-Sparse Regions Using GRACE/GRACE-FO Data Assimilation" Remote Sensing 13, no. 2: 235. https://0-doi-org.brum.beds.ac.uk/10.3390/rs13020235

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop