Improving estimates of water resources in a semi-arid region by assimilating GRACE data into the PCR-GLOBWB hydrological model

Abstract. An accurate estimation of water resources dynamics is crucial for proper management of both agriculture and the local ecology, particularly in semi-arid regions. Imperfections in model physics, uncertainties in model land parameters and meteorological data, as well as the human impact on land changes often limit the accuracy of hydrological models in estimating water storages. To mitigate this problem, this study investigated the assimilation of terrestrial water storage variation (TWSV) estimates derived from the Gravity Recovery And Climate Experiment (GRACE) data using an ensemble Kalman filter (EnKF) approach. The region considered was the Hexi Corridor in northern China. The hydrological model used for the analysis was PCR-GLOBWB, driven by satellite-based forcing data from April 2002 to December 2010. The impact of the GRACE data assimilation (DA) scheme was evaluated in terms of the TWSV, as well as the variation of individual hydrological storage estimates. The capability of GRACE DA to adjust the storage level was apparent not only for the entire TWSV but also for the groundwater component. In this study, spatially correlated errors in GRACE data were taken into account, utilizing the full error variance–covariance matrices provided as a part of the GRACE data product. The benefits of this approach were demonstrated by comparing the EnKF results obtained with and without taking into account error correlations. The results were validated against in situ groundwater data from five well sites. On average, the experiments showed that GRACE DA improved the accuracy of groundwater storage estimates by as much as 25 %. The inclusion of error correlations provided an equal or greater improvement in the estimates. In contrast, a validation against in situ streamflow data from two river gauges showed no significant benefits of GRACE DA. This is likely due to the limited spatial and temporal resolution of GRACE observations. Finally, results of the GRACE DA study were used to assess the status of water resources over the Hexi Corridor over the considered 9-year time interval. Areally averaged values revealed that TWS, soil moisture, and groundwater storages over the region decreased with an average rate of approximately 0.2, 0.1, and 0.1 cm yr−1 in terms of equivalent water heights, respectively. A particularly rapid decline in TWS (approximately −0.4 cm yr−1) was seen over the Shiyang River basin located in the southeastern part of Hexi Corridor. The reduction mostly occurred in the groundwater layer. An investigation of the relationship between water resources and agricultural activities suggested that groundwater consumption required to maintain crop yield in the growing season for this specific basin was likely the cause of the groundwater depletion.

[1]  Hubert H. G. Savenije,et al.  A Comparison of Global and Regional GRACE Models for Land Hydrology , 2008 .

[2]  Jeffrey P. Walker,et al.  THE GLOBAL LAND DATA ASSIMILATION SYSTEM , 2004 .

[3]  J. Zhu,et al.  Nonlinear effects of locally heterogeneous hydraulic conductivity fields on regional stream-aquifer exchanges , 2015 .

[4]  M. Watkins,et al.  GRACE Measurements of Mass Variability in the Earth System , 2004, Science.

[5]  Xu-sheng Wang,et al.  Satellite-based estimates of groundwater depletion in the Badain Jaran Desert, China , 2015, Scientific Reports.

[6]  Zizhan Zhang,et al.  Long-Term Groundwater Variations in Northwest India From Satellite Gravity Measurements , 2014 .

[7]  Geir Evensen,et al.  The Ensemble Kalman Filter: theoretical formulation and practical implementation , 2003 .

[8]  I. Shiklomanov,et al.  Assessment of water resources and water availability in the world , 1997 .

[9]  C. Jekeli.  Alternative methods to smooth the Earth's gravity field , 1981 .

[10]  Zong-Liang Yang,et al.  Retrieving snow mass from GRACE terrestrial water storage change with a land surface model , 2007 .

[11]  G. Wang,et al.  Environmental degradation in the Hexi Corridor region of China over the last 50 years and comprehensive mitigation and rehabilitation strategies , 2003 .

[12]  M. Bierkens,et al.  Global modeling of withdrawal, allocation and consumptive use of surface water and groundwater resources , 2013 .

[13]  Robin Wardlaw,et al.  Application of Multi-Criterion Decision Making Analysis to Integrated Water Resources Management , 2013, Water Resources Management.

[14]  J. Kusche,et al.  A systematic impact assessment of GRACE error correlation on data assimilation in hydrological models , 2016, Journal of Geodesy.

[15]  S. Swenson,et al.  Post‐processing removal of correlated errors in GRACE data , 2006 .

[16]  Zong-Liang Yang,et al.  Multisensor snow data assimilation at the continental scale: The value of Gravity Recovery and Climate Experiment terrestrial water storage information , 2010 .

[17]  J. Kusche,et al.  Calibration/Data Assimilation Approach for Integrating GRACE Data into the WaterGAP Global Hydrology Model (WGHM) Using an Ensemble Kalman Filter: First Results , 2013, Surveys in Geophysics.

[18]  B. Scanlon,et al.  Uncertainty in evapotranspiration from land surface modeling, remote sensing, and GRACE satellites , 2014 .

[19]  A. Huete,et al.  Overview of the radiometric and biophysical performance of the MODIS vegetation indices , 2002 .

[20]  T. McMahon,et al.  Updated world map of the Köppen-Geiger climate classification , 2007 .

[21]  Brian C. Gunter,et al.  Data assimilation of GRACE terrestrial water storage estimates into a regional hydrological model of the Rhine River basin , 2014 .

[22]  M. Rodell,et al.  Assimilation of gridded terrestrial water storage observations from GRACE into a land surface model , 2016 .

[23]  P. Döll,et al.  Global‐scale assessment of groundwater depletion and related groundwater abstractions: Combining hydrological modeling with information from well observations and GRACE satellites , 2014 .

[24]  Grzegorz Michalak,et al.  GFZ RL05: An Improved Time-Series of Monthly GRACE Gravity Field Solutions , 2014 .

[25]  M. Bierkens,et al.  Global monthly water stress: 1. Water balance and water availability , 2011 .

[26]  Y. Hong,et al.  The TRMM Multisatellite Precipitation Analysis (TMPA): Quasi-Global, Multiyear, Combined-Sensor Precipitation Estimates at Fine Scales , 2007 .

[27]  F. Bryan,et al.  Time variability of the Earth's gravity field: Hydrological and oceanic effects and their possible detection using GRACE , 1998 .

[28]  S. M. de Jong,et al.  Calibrating a large‐extent high‐resolution coupled groundwater‐land surface model using soil moisture and discharge data , 2014 .

[29]  M. Cheng,et al.  Variations in the Earth's oblateness during the past 28 years , 2004 .

[30]  George J. Huffman,et al.  Estimates of Root-Mean-Square Random Error for Finite Samples of Estimated Precipitation , 1997 .

[31]  Frank Flechtner,et al.  Status of the GRACE Follow-On Mission , 2013 .

[32]  Shaozhong Kang,et al.  China's food security is threatened by the unsustainable use of water resources in North and Northwest China , 2014 .

[33]  M. Rodell,et al.  Assimilation of terrestrial water storage from GRACE in a snow‐dominated basin , 2011 .

[34]  H. Malano,et al.  Seasonal Water Allocation: Dealing with Hydrologic Variability in the Context of a Water Rights System , 2013 .

[35]  Jasmeet Judge,et al.  A particle batch smoother for soil moisture estimation using soil temperature observations , 2015 .

[36]  R. Kalin,et al.  Multi-objective optimization for sustainable groundwater resource management in a semiarid catchment , 2001 .

[37]  S. Swenson,et al.  Accuracy of GRACE mass estimates , 2006 .

[38]  T. D. Mitchell,et al.  An improved method of constructing a database of monthly climate observations and associated high‐resolution grids , 2005 .

[39]  J. Shao,et al.  The Role of Ground Water in Arid/Semiarid Ecosystems, Northwest China , 2005, Ground water.

[40]  D. Chambers,et al.  Estimating Geocenter Variations from a Combination of GRACE and Ocean Model Output , 2008 .

[41]  Niko E. C. Verhoest,et al.  Satellite-Scale Snow Water Equivalent Assimilation into a High-Resolution Land Surface Model , 2010 .

[42]  M. Rodell,et al.  Assimilation of GRACE Terrestrial Water Storage Data into a Land Surface Model: Results for the Mississippi River Basin , 2008 .

[43]  S. M. de Jong,et al.  Large-scale groundwater modeling using global datasets: a test case for the Rhine-Meuse basin , 2011 .

[44]  J. Thepaut,et al.  The ERA‐Interim reanalysis: configuration and performance of the data assimilation system , 2011 .

[45]  A. Weerts,et al.  Particle filtering and ensemble Kalman filtering for state updating with hydrological conceptual rainfall‐runoff models , 2006 .

[46]  H. Steinfeld,et al.  Livestock's long shadow: environmental issues and options. , 2006 .

[47]  S. Bettadpur.  Insights into the Earth System mass variability from CSR-RL05 GRACE gravity fields , 2012 .

[48]  T. Carlson,et al.  On the relation between NDVI, fractional vegetation cover, and leaf area index , 1997 .

[49]  E. Wood,et al.  Development of a 50-Year High-Resolution Global Dataset of Meteorological Forcings for Land Surface Modeling , 2006 .

[50]  R. Koster,et al.  Assimilation of GRACE terrestrial water storage into a land surface model: Evaluation and potential value for drought monitoring in western and central Europe , 2012 .

[51]  Randal D. Koster,et al.  Bias reduction in short records of satellite soil moisture , 2004 .

[52]  R. Houborg,et al.  Drought indicators based on model‐assimilated Gravity Recovery and Climate Experiment (GRACE) terrestrial water storage observations , 2012 .

[53]  M. Bierkens,et al.  The benefits of using remotely sensed soil moisture in parameter identification of large‐scale hydrological models , 2013 .

[54]  C. Kummerow,et al.  The Tropical Rainfall Measuring Mission (TRMM) Sensor Package , 1998 .

[55]  J. Whitaker,et al.  Distance-dependent filtering of background error covariance estimates in an ensemble Kalman filter , 2001 .

[56]  P. Shi,et al.  Daily precipitation changes in the semi-arid region over northern China , 2004 .

[57]  X. S. Wang,et al.  The characteristics of ground-water resources and their changes under the impacts of human activity in the arid Northwest China—a case study of the Shiyang River Basin , 2005 .

[58]  Niels Drost,et al.  PCR-GLOBWB 2: a 5 arc-minute global hydrologicaland water resources model , 2017 .

[59]  Brian C. Gunter,et al.  Assessing total water storage and identifying flood events over Tonlé Sap basin in Cambodia using GRACE and MODIS satellite observations combined with hydrological models , 2016 .

[60]  Z. Huo,et al.  The Response of Water-Land Environment to Human Activities in Arid Minqin Oasis, Northwest China , 2007 .