Evaluation of Groundwater Storage Variations Estimated from GRACE Data Assimilation and State-of-the-Art Land Surface Models in Australia and the North China Plain

The accurate knowledge of the groundwater storage variation (ΔGWS) is essential for reliable water resource assessment, particularly in arid and semi-arid environments (e.g., Australia, the North China Plain (NCP)) where water storage is significantly affected by human activities and spatiotemporal climate variations. The large-scale ΔGWS can be simulated from a land surface model (LSM), but the high model uncertainty is a major drawback that reduces the reliability of the estimates. The evaluation of the model estimate is then very important to assess its accuracy. To improve the model performance, the terrestrial water storage variation derived from the Gravity Recovery And Climate Experiment (GRACE) satellite mission is commonly assimilated into LSMs to enhance the accuracy of the ΔGWS estimate. This study assimilates GRACE data into the PCRaster Global Water Balance (PCR-GLOBWB) model. The GRACE data assimilation (DA) is developed based on the three-dimensional ensemble Kalman smoother (EnKS 3D), which considers the statistical correlation of all extents (spatial, temporal, vertical) in the DA process. The ΔGWS estimates from GRACE DA and four LSM simulations (PCR-GLOBWB, the Community Atmosphere Biosphere Land Exchange (CABLE), the Water Global Assessment and Prognosis Global Hydrology Model (WGHM), and World-Wide Water (W3)) are validated against the in situ groundwater data. The evaluation is conducted in terms of temporal correlation, seasonality, long-term trend, and detection of groundwater depletion. The GRACE DA estimate shows a significant improvement in all measures, notably the correlation coefficients (respect to the in situ data) are always higher than the values obtained from model simulations alone (e.g., ~0.15 greater in Australia, and ~0.1 greater in the NCP). GRACE DA also improves the estimation of groundwater depletion that the models cannot accurately capture due to the incorrect information of the groundwater demand (in, e.g., PCR-GLOBWB, WGHM) or the unavailability of a groundwater consumption routine (in, e.g., CABLE, W3). In addition, this study conducts the inter-comparison between four model simulations and reveals that PCR-GLOBWB and CABLE provide a more accurate ΔGWS estimate in Australia (subject to the calibrated parameter) while PCR-GLOBWB and WGHM are more accurate in the NCP (subject to the inclusion of anthropogenic factors). The analysis can be used to declare the status of the ΔGWS estimate, as well as itemize the possible improvements of the future model development.

[1]  W. Feng,et al.  Evaluation of groundwater depletion in North China using the Gravity Recovery and Climate Experiment (GRACE) data and ground‐based measurements , 2013 .

[2]  P. Döll,et al.  Groundwater use for irrigation - a global inventory , 2010 .

[3]  J. Kusche,et al.  Improving drought simulations within the Murray-Darling Basin by combined calibration/assimilation of GRACE data into the WaterGAP Global Hydrology Model , 2018 .

[4]  Wenji Zhao,et al.  Subregional‐scale groundwater depletion detected by GRACE for both shallow and deep aquifers in North China Plain , 2015 .

[5]  Xubin Zeng,et al.  Impact of Modified Richards Equation on Global Soil Moisture Simulation in the Community Land Model (CLM3.5) , 2009 .

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

[7]  V. Zlotnicki,et al.  Australian water mass variations from GRACE data linked to Indo-Pacific climate variability , 2011 .

[8]  M. Rodell,et al.  Benefits and pitfalls of GRACE data assimilation: A case study of terrestrial water storage depletion in India , 2017, Geophysical research letters.

[9]  O. Batelaan,et al.  Estimation of GRACE water storage components by temporal decomposition. , 2017 .

[10]  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 .

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

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

[13]  I. Hoteit,et al.  A two-update ensemble Kalman filter for land hydrological data assimilation with an uncertain constraint. , 2017 .

[14]  Byron D. Tapley,et al.  Long-term groundwater storage change in Victoria, Australia from satellite gravity and in situ observations , 2016 .

[15]  Y. Hong,et al.  Have GRACE satellites overestimated groundwater depletion in the Northwest India Aquifer? , 2016, Scientific Reports.

[16]  A. Pitman,et al.  Modelling evapotranspiration during precipitation deficits: identifying critical processes in a land surface model , 2015 .

[17]  Bridget R. Scanlon,et al.  Use of flow modeling to assess sustainability of groundwater resources in the North China Plain , 2013 .

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

[19]  J. Kusche,et al.  Separation of large scale water storage patterns over Iran using GRACE, altimetry and hydrological data , 2014 .

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

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

[22]  K. Lambeck,et al.  GRACE estimates of sea surface height anomalies in the Gulf of Carpentaria, Australia , 2008 .

[23]  P. Döll,et al.  A global hydrological model for deriving water availability indicators: model tuning and validation , 2003 .

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

[25]  A. Dijk,et al.  Accounting for spatial correlation errors in the assimilation of GRACE into hydrological models through localization , 2017 .

[26]  Guillaume Ramillien,et al.  Basin‐scale, integrated observations of the early 21st century multiyear drought in southeast Australia , 2009 .

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

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

[29]  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.

[30]  Chuang Xu,et al.  Characterizing Drought and Flood Events over the Yangtze River Basin Using the HUST-Grace2016 Solution and Ancillary Data , 2017, Remote. Sens..

[31]  P. Döll,et al.  Sensitivity of simulated global-scale freshwater fluxes and storages to input data, hydrological model structure, human water use and calibration , 2014 .

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

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

[34]  F. Chiew,et al.  Modelling climate-change impacts on groundwater recharge in the Murray-Darling Basin, Australia , 2010 .

[36]  M. Decker Development and evaluation of a new soil moisture and runoff parameterization for the CABLE LSM including subgrid‐scale processes , 2015 .

[37]  Edwin Sutanudjaja,et al.  Towards a global land subsidence map , 2015 .

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

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

[40]  Ibrahim Hoteit,et al.  Assessing sequential data assimilation techniques for integrating GRACE data into a hydrological model , 2017 .

[41]  J. Famiglietti,et al.  Satellite-based estimates of groundwater depletion in India , 2009, Nature.

[42]  Eric F. Wood,et al.  Global analysis of seasonal streamflow predictability using an ensemble prediction system and observations from 6192 small catchments worldwide , 2013 .

[43]  B. Timbal,et al.  The Millennium Drought in southeast Australia (2001–2009): Natural and human causes and implications for water resources, ecosystems, economy, and society , 2013 .

[44]  Martha C. Anderson,et al.  Assimilation of Gridded GRACE Terrestrial Water Storage Estimates in the North American Land Data Assimilation System , 2016 .

[45]  Impact of climate forcing uncertainty and human water use on global and continental water balance components , 2016 .

[46]  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 .

[47]  J. McCallum,et al.  Impacts of climate change on groundwater in Australia: a sensitivity analysis of recharge , 2010 .

[48]  M. Webber,et al.  The South–North Water Transfer Project: remaking the geography of China , 2017 .

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

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

[51]  Abhijit Mukherjee,et al.  A Data Assimilation Perspective on India’s Terrestrial Water Storage Trends , 2017 .

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

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

[54]  P. Tregoning,et al.  Tropical cyclones and the ecohydrology of Australia's recent continental‐scale drought , 2012 .

[55]  A. Mostafaie,et al.  Large-Scale Total Water Storage and Water Flux Changes over the Arid and Semiarid Parts of the Middle East from GRACE and Reanalysis Products , 2017, Surveys in Geophysics.

[56]  G. Balsamo,et al.  The WFDEI meteorological forcing data set: WATCH Forcing Data methodology applied to ERA‐Interim reanalysis data , 2014 .

[57]  Khandu,et al.  Exploring the influence of precipitation extremes and human water use on total water storage (TWS) changes in the Ganges‐Brahmaputra‐Meghna River Basin , 2016 .

[58]  S. Swenson,et al.  Satellites measure recent rates of groundwater depletion in California's Central Valley , 2011 .

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

[60]  Shin‐Chan Han,et al.  On the use of the GRACE normal equation of inter-satellite tracking data for estimation of soil moisture and groundwater in Australia , 2018 .

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

[62]  Brian C. Gunter,et al.  Improving estimates of water resources in a semi-arid region by assimilating GRACE data into the PCR-GLOBWB hydrological model , 2016 .

[63]  E. Kowalczyk,et al.  The ACCESS coupled model: description, control climate and evaluation , 2013 .

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

[65]  P. Tregoning,et al.  Assessment of GRACE satellites for groundwater estimation in Australia , 2012 .

[66]  Eloise Kendy,et al.  Groundwater depletion: A global problem , 2005 .

[67]  Enhancement of water storage estimates using GRACE data assimilation with particle filter framework , 2017 .

[68]  Ge Sun,et al.  A COMPARISON OF SIX POTENTIAL EVAPOTRANSPIRATION METHODS FOR REGIONAL USE IN THE SOUTHEASTERN UNITED STATES 1 , 2005 .

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

[70]  Petra Döll,et al.  Value of river discharge data for global-scale hydrological modeling , 2007 .