Scaling methods of leakage correction in GRACE mass change estimates revisited for the complex hydro-climatic setting of the Indus basin
- 1Hydro-Remote Sensing Applications (H-RSA) Group, Department of Civil Engineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra, India
- 2Institut für Planetare Geodäsie, Technische Universität Dresden, Germany
- 1Hydro-Remote Sensing Applications (H-RSA) Group, Department of Civil Engineering, Indian Institute of Technology Bombay, Mumbai, Maharashtra, India
- 2Institut für Planetare Geodäsie, Technische Universität Dresden, Germany
Abstract. Total Water Storage Change (TWSC) reflects the balance of all water fluxes in a hydrological system. The Gravity Recovery and Climate Experiment/Follow-On (GRACE/GRACE-FO) satellite missions are the only means of observing this state variable, distributed as coefficients of a spherical harmonic (SH) model. The well-known correlated noise in these observations requires filtering, scattering the actual mass changes from their true locations. This effect is known as leakage. This study explores the traditional basin and grid scaling approaches and develops a novel frequency-dependent scaling for leakage correction of GRACE TWSC in a unique, basin-specific assessment for the Indus basin. We harness the characteristics of significant heterogeneity in the Indus basin due to climate and human-induced changes to evaluate the physical nature of these scaling schemes. The most recent, WaterGAP Hydrology Model (WGHM v2.2d), with its two variants, standard (without glacier mass changes) and Integrated (with glacier mass changes), is used to derive scaling factors. For the first time, we explicitly show the effect of inclusion or exclusion of glacier mass changes in the model on the gridded scaling factors. Frequency dependant scaling factors, as a novelty, allow us to compare the differences of scaling seasonal and trend components separately. We employ time series analysis using the Lomb-Scargle periodogram to decompose the time series and visualize the effect of filtering on different time scales. We find that for the Indus basin, where mass changes of different frequencies are localized, frequency-dependent scaling factors for the entire basin perform similar to the gridded scaling factors, while both outperform basin scaling. The property of frequency-dependent scaling to keep the noise level unscaled can be extremely useful for applications requiring a high signal-to-noise ratio of TWSC observations. Apart from these novel developments and insights into the traditional scaling approach, our study encourages the regional scale users to conduct specific assessments for their basin of interest.
- Preprint
(2434 KB) -
Supplement
(71 KB) - BibTeX
- EndNote
Vasaw Tripathi et al.
Status: final response (author comments only)
- RC1: 'Comment on hess-2022-43', A. P. Dimri, 11 Mar 2022
-
RC2: 'Review on hess-2022-43', Henryk Dobslaw, 31 Mar 2022
I read with interest the manuscript
Scaling methods of leakage correction in GRACE mass change estimates revisited for the complex hydro-climatic setting of the Indus basin
written by Vasaw Tripathi and colleagues. The paper adresses the problem of spatial leakage in satellite gravimetry data which arises from the limited ability of sensor data from low-low satellite tracking to accurately resolve steep spatial gradients in surface mass anomalies at scales of a few hundred kilometers and smaller. The resulting systematic error (called spatial leakage) is often mitigated by means of scaling approaches, and the current article is following this avenue of research with a special emphasis on the Indus Basin. The paper is generally well written and complements the existing literature on applications of satellite gravimetry data in this catchment. The study moreover fits nicely into the scope of the journal so that acceptance might be recommended as soon as a number of concerns outlined below are properly addressed.
(1) Rescaling (as performed in this study) is typically applied to allow for regional or even small-scale applications of the GRACE data. For the whole Indus catchment of about 1 million square kilometers, the effect of rescaling should be rather minor (which is also confirmed in the present work by the rather small changes in the time-series when moving from Figure 4 to Figures 14 and 15). It would be thus imperative to discuss in more detail the re-scaled GRACE-results presented in Figure 13, and also to compare them to the CSR mascons pixel by pixel. Note that I do not expect that rescaling will lead to a perfect match with the Mascons, nor do I suggest that the Mascons should be considered as the error-free truth in such an exercise. Instead, it would be important to carve out more explicitly any benefits of applying an elaborate rescaling scheme in lieu of simply downloading and applying the Mascons.
(2) A central result of this study are the scaling factors given in Figure 10. Table 4 gives a nice overview on the interpretation, which should be also reflected in the color scale selected for the figure. Please consider something like bright red or purple for all pixels with zero or even negative scaling, and neutral colors (green or white) for coefficients around one.
(3) Data-driven methods to account for spatial leakage have been proposed in the past by some of the authors (e.g., 10.1002/2017WR021150), and such methods are now routinely applied to approximate spatial leakage in surface mass estimates obtained from GRACE spherical harmonics solutions as disseminated via gravis.gfz-potsdam.de. Since Gaussian filtering is readily available to all the authors of the present paper, it would be quite straightforward to additionally explore the usage of twice-filtered GRACE GSM fields for the assessment of spatial leakage in the Indus Catchment. I believe that such an additional experiment could nicely complement the existing material.
(4) I understand that spatially distributed in situ data is not readily available in a transboundary basin like the Indus catchment. In such a situation, a common approach to demonstrate the applicability of a new method developments would be the usage of simulated sensor data and satellite products, where the true mass variability entering the simulations is known (see, e.g., 10.1007/s10712-015-9338-y). In such a simulation environment, it could be demonstrated to what extent rescaling as proposed here mitigates the adverse effects of spatial leakage.
(5) I am not quite convinced that frequency-dependent scaling (as proposed here) will be "useful for applications requiring a high signal-to-noise ratio". After all, it is any deviation from the seasonal cycle related to either interannual climate variations or hydrometeorological extremes that is of particular interest in many applications of remote sensing data. I do not find arguments in the paper that would help restoring leakage for such signals, so please elaborate your claim a little further.
A number of minor points might be additionally considered during the revision:
line 29: It is surprising to find a PhD thesis cited for such a rather general statement. In case you would like to give credit to the work of this author specifically, please consider citing any of her research papers.
Line 44: Authors should understand that even perfect dealiasing models will not remove all spatially correlated errors. On the contrary, simulation studies (10.1007/s10712-015-9338-y) demonstrate that background model errors in tides and sensor noise of the GRACE accelerometers have an almost similarly large effect on accuracy and spatial resolution of the GRACE monthly solutions available today. This also applies to GRACE-FO, where one of the accelerometer instruments performs worse than expected.
Line 49: Spatial leakage errors do not occur from the truncation of the spherical harmonics expansion but from the limited resolution of the along-track sensor data and the upward continuation of the gravity field from the surface to the orbital height of the satellites. Further, there is some inherent smoothing in the sensor data pre-processing that also reduces spatial gradients in the resulting gravity field estimates. Expanding the cutoff degree will certainly not solve the leakage issues, but render the inversion problem ill-posed.
Line 97: dependent
Line 104: Incomplete sentence. I assume that those glaciers are additionally feeding the Indus via tributary rivers?
Table 1: I am not convinced that "Integrated WGHM" is a very intuitive name for this model experiment. What about WGHM+GGM?
Line 249: In Table 1, eight different model experiments are listed: Please explain more explicitly why just two of them are applied in the following.
- EC1: 'Comment on hess-2022-43', Narendra Das, 19 May 2022
Vasaw Tripathi et al.
Vasaw Tripathi et al.
Viewed
HTML | XML | Total | Supplement | BibTeX | EndNote | |
---|---|---|---|---|---|---|
370 | 93 | 15 | 478 | 22 | 6 | 6 |
- HTML: 370
- PDF: 93
- XML: 15
- Total: 478
- Supplement: 22
- BibTeX: 6
- EndNote: 6
Viewed (geographical distribution)
Country | # | Views | % |
---|
Total: | 0 |
HTML: | 0 |
PDF: | 0 |
XML: | 0 |
- 1