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.

Groundwater is the world's largest accessible source of fresh water. It plays a vital role in satisfying basic needs for drinking water, agriculture and industrial activities. During times of drought groundwater sustains baseflow to rivers and wetlands, thereby supporting ecosystems. Most global-scale hydrological models (GHMs) do not include a groundwater flow component, mainly due to lack of geohydrological data at the global scale. For the simulation of lateral flow and groundwater head dynamics, a realistic physical representation of the groundwater system is needed, especially for GHMs that run at finer resolutions. In this study we present a global-scale groundwater model (run at 6' resolution) using MODFLOW to construct an equilibrium water table at its natural state as the result of long-term climatic forcing. The used aquifer schematization and properties are based on available global data sets of lithology and transmissivities combined with the estimated thickness of an upper, unconfined aquifer. This model is forced with outputs from the land-surface PCRaster Global Water Balance (PCR-GLOBWB) model, specifically net recharge and surface water levels. A sensitivity analysis, in which the model was run with various parameter settings, showed that variation in saturated conductivity has the largest impact on the groundwater levels simulated. Validation with observed groundwater heads showed that groundwater heads are reasonably well simulated for many regions of the world, especially for sediment basins (R² = 0.95). The simulated regional-scale groundwater patterns and flow paths demonstrate the relevance of lateral groundwater flow in GHMs. Inter-basin groundwater flows can be a significant part of a basin's water budget and help to sustain river baseflows, especially during droughts. Also, water availability of larger aquifer systems can be positively affected by additional recharge from inter-basin groundwater flows.

A considerable number of river basins around the world lack sufficient ground observations of hydro-meteorological 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 situ-observed 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 Multi-Source 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.

[1] We explore the possibility of using remotely sensed soil moisture data and in situ discharge observations to calibrate a large-extent hydrological model. The model used is PCR-GLOBWB-MOD, which is a physically based and fully coupled groundwater-land surface model operating at a daily basis and having a resolution of 30 arc sec (about 1 km at the equator). As a test bed, we use the combined Rhine-Meuse basin (total area: about 200,000 km²), where there are 4250 point-scale observed groundwater head time series that are used to verify the model results. Calibration is performed by simulating 3045 model runs with varying parameter values affecting groundwater head dynamics. The simulation results of all runs are evaluated against the remotely sensed soil moisture time series of SWI (Soil Water Index) and field discharge data. The former is derived from European Remote Sensing scatterometers and provides estimates of the first meter profile soil moisture content at 30 arc min resolution (50 km at the equator). From the evaluation of these runs, we then introduce a stepwise calibration approach that considers stream discharge first, then soil moisture, and finally verify the resulting simulation to groundwater head observations. Our results indicate that the remotely sensed soil moisture data can be used for the calibration of upper soil hydraulic conductivities determining simulated groundwater recharge of the model. However, discharge data should be included to obtain full calibration of the coupled model, specifically to constrain aquifer transmissivities and runoff-infiltration partitioning processes. The stepwise approach introduced in this study, using both discharge and soil moisture data, can calibrate both discharge and soil moisture, as well as predicting groundwater head dynamics with acceptable accuracy. As our approach to parameterize and calibrate the model uses globally available data sets only, it opens up the possibility to set up large-extent coupled groundwater-land surface models in other basins or even globally.

Land surface and hydrologic models (LSMs/HMs) are used at diverse spatial resolutions ranging from catchment-scale (1–10 km) to global-scale (over 50 km) applications. Applying the same model structure at different spatial scales requires that the model estimates similar fluxes independent of the chosen resolution, i.e., fulfills a flux-matching condition across scales. An analysis of state-of-the-art LSMs and HMs reveals that most do not have consistent hydrologic parameter fields. Multiple experiments with the mHM, Noah-MP, PCR-GLOBWB, and WaterGAP models demonstrate the pitfalls of deficient parameterization practices currently used in most operational models, which are insufficient to satisfy the flux-matching condition. These examples demonstrate that J. Dooge's 1982 statement on the unsolved problem of parameterization in these models remains true. Based on a review of existing parameter regionalization techniques, we postulate that the multiscale parameter regionalization (MPR) technique offers a practical and robust method that provides consistent (seamless) parameter and flux fields across scales. Herein, we develop a general model protocol to describe how MPR can be applied to a particular model and present an example application using the PCR-GLOBWB model. Finally, we discuss potential advantages and limitations of MPR in obtaining the seamless prediction of hydrological fluxes and states across spatial scales.