Abstract. We present PCR-GLOBWB 2, a global hydrology and water resources model. Compared to previous versions of PCR-GLOBWB, this version fully integrates water use. Sector-specific water demand, groundwater and surface water withdrawal, water consumption, and return flows are dynamically calculated at every time step and interact directly with the simulated hydrology. PCR-GLOBWB 2 has been fully rewritten in Python and PCRaster Python and has a modular structure, allowing easier replacement, maintenance, and development of model components. PCR-GLOBWB 2 has been implemented at 5 arcmin resolution, but a version parameterized at 30 arcmin resolution is also available. Both versions are available as open-source codes on https://github.com/UU-Hydro/PCR-GLOBWB_model (Sutanudjaja et al., 2017a). PCR-GLOBWB 2 has its own routines for groundwater dynamics and surface water routing. These relatively simple routines can alternatively be replaced by dynamically coupling PCR-GLOBWB 2 to a global two-layer groundwater model and 1-D–2-D hydrodynamic models. Here, we describe the main components of the model, compare results of the 30 and 5 arcmin versions, and evaluate their model performance using Global Runoff Data Centre discharge data. Results show that model performance of the 5 arcmin version is notably better than that of the 30 arcmin version. Furthermore, we compare simulated time series of total water storage (TWS) of the 5 arcmin model with those observed with GRACE, showing similar negative trends in areas of prevalent groundwater depletion. Also, we find that simulated total water withdrawal matches reasonably well with reported water withdrawal from AQUASTAT, while water withdrawal by source and sector provide mixed results.

Abstract. A considerable number of river basins around the world lack sufficient ground observations of hydrometeorological data for effective water resources assessment and management. Several approaches can be developed to increase the quality and availability of data in these poorly gauged or ungauged river basins; among them, the use of Earth observations products has recently become promising. Earth observations of various environmental variables can be used potentially to increase knowledge about the hydrological processes in the basin and to improve streamflow model estimates, via assimilation or calibration. The present study aims to calibrate the large-scale hydrological model PCRaster GLOBal Water Balance (PCR-GLOBWB) using satellite-based products of evapotranspiration and soil moisture for the Moroccan Oum er Rbia River basin. Daily simulations at a spatial resolution of 5 × 5 arcmin are performed with varying parameters values for the 32-year period 1979-2010. Five different calibration scenarios are inter-compared: (i) reference scenario using the hydrological model with the standard parameterization, (ii) calibration using in situobserved discharge time series, (iii) calibration using the Global Land Evaporation Amsterdam Model (GLEAM) actual evapotranspiration time series, (iv) calibration using ESA Climate Change Initiative (CCI) surface soil moisture time series and (v) step-wise calibration using GLEAM actual evapotranspiration and ESA CCI surface soil moisture time series. The impact on discharge estimates of precipitation in comparison with model parameters calibration is investigated using three global precipitation products, including ERA-Interim (EI), WATCH Forcing methodology applied to ERA-Interim reanalysis data (WFDEI) and MultiSource Weighted-Ensemble Precipitation data by merging gauge, satellite and reanalysis data (MSWEP).Results show that GLEAM evapotranspiration and ESA CCI soil moisture may be used for model calibration resulting in reasonable discharge estimates (NSE values from 0.5 to 0.75), although better model performance is achieved when the model is calibrated with in situ streamflow observations. Independent calibration based on only evapotranspiration or soil moisture observations improves model predictions to a lesser extent. Precipitation input affects discharge estimates more than calibrating model parameters. The use of WFDEI precipitation leads to the lowest model performances. Apart from the in situ discharge calibration scenario, the highest discharge improvement is obtained when EI and MSWEP precipitation products are used in combination with a step-wise calibration approach based on evapotranspiration and soil moisture observations. This study opens up the possibility of using globally available Earth observations and reanalysis products of precipitation, evapotranspiration and soil moisture in large-scale hydrological models to estimate discharge at a river basin scale.

Abstract. The coarse spatial resolution of global hydrological models (typically  >  0.25°) limits their ability to resolve key water balance processes for many river basins and thus compromises their suitability for water resources management, especially when compared to locally tuned river models. A possible solution to the problem may be to drive the coarse-resolution models with locally available high-spatial-resolution meteorological data as well as to assimilate ground-based and remotely sensed observations of key water cycle variables. While this would improve the resolution of the global model, the impact of prediction accuracy remains largely an open question. In this study, we investigate the impact of assimilating streamflow and satellite soil moisture observations on the accuracy of global hydrological model estimations, when driven by either coarse- or high-resolution meteorological observations in the Murrumbidgee River basin in Australia. To this end, a 0.08° resolution version of the PCR-GLOBWB global hydrological model is forced with downscaled global meteorological data (downscaled from 0.5° to 0.08° resolution) obtained from the WATCH Forcing Data methodology applied to ERA-Interim (WFDEI) and a local high-resolution, gauging-station-based gridded data set (0.05°). Downscaled satellite-derived soil moisture (downscaled from  ∼  0.5° to 0.08° resolution) from the remote observation system AMSR-E and streamflow observations collected from 23 gauging stations are assimilated using an ensemble Kalman filter. Several scenarios are analysed to explore the added value of data assimilation considering both local and global meteorological data. Results show that the assimilation of soil moisture observations results in the largest improvement of the model estimates of streamflow. The joint assimilation of both streamflow and downscaled soil moisture observations leads to further improvement in streamflow simulations (20 % reduction in RMSE). Furthermore, results show that the added contribution of data assimilation, for both soil moisture and streamflow, is more pronounced when the global meteorological data are used to force the models. This is caused by the higher uncertainty and coarser resolution of the global forcing. We conclude that it is possible to improve PCR-GLOBWB simulations forced by coarse-resolution meteorological data with assimilation of downscaled spaceborne soil moisture and streamflow observations. These improved model results are close to the ones from a local model forced with local meteorological data. These findings are important in light of the efforts that are currently made to move to global hyper-resolution modelling and can help to advance this research.