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. Extreme flooding impacts millions of people that live within the Amazon floodplain. Global hydrological models (GHMs) are frequently used to assess and inform the management of flood risk, but knowledge on the skill of available models is required to inform their use and development. This paper presents an intercomparison of eight different GHMs freely available from collaborators of the Global Flood Partnership (GFP) for simulating floods in the Amazon basin. To gain insight into the strengths and shortcomings of each model, we assess their ability to reproduce daily and annual peak river flows against gauged observations at 75 hydrological stations over a 19-year period (1997–2015). As well as highlighting regional variability in the accuracy of simulated streamflow, these results indicate that (a) the meteorological input is the dominant control on the accuracy of both daily and annual maximum river flows, and (b) groundwater and routing calibration of Lisflood based on daily river flows has no impact on the ability to simulate flood peaks for the chosen river basin. These findings have important relevance for applications of large-scale hydrological models, including analysis of the impact of climate variability, assessment of the influence of long-term changes such as land-use and anthropogenic climate change, the assessment of flood likelihood, and for flood forecasting systems.

Abstract. We here present GLOFRIM, a globally applicable computational framework for integrated hydrologicalhydrodynamic modelling. GLOFRIM facilitates spatially explicit coupling of hydrodynamic and hydrologic models and caters for an ensemble of models to be coupled. It currently encompasses the global hydrological model PCR-GLOBWB as well as the hydrodynamic models Delft3D Flexible Mesh (DFM; solving the full shallow-water equations and allowing for spatially flexible meshing) and LISFLOOD-FP (LFP; solving the local inertia equations and running on regular grids). The main advantages of the framework are its open and free access, its global applicability, its versatility, and its extensibility with other hydrological or hydrodynamic models. Before applying GLOFRIM to an actual test case, we benchmarked both DFM and LFP for a synthetic test case. Results show that for sub-critical flow conditions, discharge response to the same input signal is near-identical for both models, which agrees with previous studies. We subsequently applied the framework to the Amazon River basin to not only test the framework thoroughly, but also to perform a first-ever benchmark of flexible and regular grids on a large-scale. Both DFM and LFP produce comparable results in terms of simulated discharge with LFP exhibiting slightly higher accuracy as expressed by a Kling-Gupta efficiency of 0.82 compared to 0.76 for DFM. However, benchmarking inundation extent between DFM and LFP over the entire study area, a critical success index of 0.46 was obtained, indicating that the models disagree as often as they agree. Differences between models in both simulated discharge and inundation extent are to a large extent attributable to the gridding techniques employed. In fact, the results show that both the numerical scheme of the inundation model and the gridding technique can contribute to deviations in simulated inundation extent as we control for model forcing and boundary conditions. This study shows that the presented computational framework is robust and widely applicable. GLOFRIM is designed as open access and easily extendable, and thus we hope that other large-scale hydrological and hydrodynamic models will be added. Eventually, more locally relevant processes would be captured and more robust model inter-comparison, benchmarking, and ensemble simulations of flood hazard on a large scale would be allowed for.

Geosci. Model Dev. Discuss., https://doi

Abstract. Large-scale flood events often show spatial correlation in neighbouring basins, and thus can affect adjacent basins simultaneously, as well as result in superposition of different flood peaks. Such flood events therefore need to be addressed with large-scale modelling approaches to capture these processes. Many approaches currently in place are based on either a hydrologic or a hydrodynamic model. However, the resulting lack of interaction between hydrology and hydrodynamics, for instance, by implementing groundwater infiltration on inundated floodplains, can hamper modelled inundation and discharge results where such interactions are important. In this study, the global hydrologic model PCR-GLOBWB at 30 arcmin spatial resolution was one-directionally and spatially coupled with the hydrodynamic model Delft 3D Flexible Mesh (FM) for the Amazon River basin at a grid-by-grid basis and at a daily time step. The use of a flexible unstructured mesh allows for fine-scale representation of channels and floodplains, while preserving a coarser spatial resolution for less floodprone areas, thus not unnecessarily increasing computational costs. In addition, we assessed the difference between a 1-D channel/2-D floodplain and a 2-D schematization in Delft 3D FM. Validating modelled discharge results shows that coupling PCR-GLOBWB to a hydrodynamic routing scheme generally increases model performance compared to using a hydrodynamic or hydrologic model only for all validation parameters applied. Closer examination shows that the 1-D/2-D schematization outperforms 2-D for r 2 and root mean square error (RMSE) whilst having a lower KlingGupta efficiency (KGE). We also found that spatial coupling has the significant advantage of a better representation of inundation at smaller streams throughout the model domain. A validation of simulated inundation extent revealed that only those set-ups incorporating 1-D channels are capable of representing inundations for reaches below the spatial resolution of the 2-D mesh. Implementing 1-D channels is therefore particularly of advantage for large-scale inundation models, as they are often built upon remotely sensed surface elevation data which often enclose a strong vertical bias, hampering downstream connectivity. Since only a one-directional coupling approach was tested, and therefore important feedback processes are not incorporated, simulated discharge and inundation extent for both coupled set-ups is generally overpredicted. Hence, it will be the subsequent step to extend it to a two-directional coupling scheme to obtain a closed feedback loop between hydrologic and hydrodynamic processes. The current findings demonstrating the potential of one-directionally and spatially coupled models to obtain improved discharge estimates form an important step towards a large-scale inundation model with a full dynamic coupling between hydrology and hydrodynamics.

In recent years, a range of global flood models (GFMs) were developed, each utilizing different process descriptions as well as validation data sets and methods. To quantify the magnitude of these differences, studies assessed the performance of GFMs only on the continental and catchment level. Since the default model set-ups resulted in locally marked deviations, there is a clear need for further and especially more standardized research to not only maintain credibility, but also support the application of GFM products by end-users. Consequently, here we conceptually outline the basic requirements and challenges of a Global Flood Model Validation Framework for more standardized model validation and benchmarking. With the proposed framework we hope to encourage the much needed debate, research developments in this direction, and involvement of science with end-users. By means of the framework, it is possible to streamline the data sets used for input and validation as well as the validation approach itself. By subjecting GFMs to more thorough and standardized methods, we think their quality as well as acceptance will increase as a result, especially amongst endusers of their outputs. Otherwise GFMs may only serve a purely scientific purpose of continued 'siloed' model improvement but without practical use. Furthermore, we want to invite GFM developers to make their models more integratable which would allow for representation of more physical processes and even more detailed comparison on a model component basis. We think this is pivotal to not only improve the accuracy of model input data sets, but to focus on the core of each model, the process descriptions. Only if we know more about why GFMs deviate, are we able to improve them accordingly and develop a next generation of models, not only providing first-order estimates of flood extent but supporting the global disaster risk reduction community with more accurate and actionable information.

Fluvial flood events are a major threat to people and infrastructure. Typically, flood hazard is driven by hydrologic or river routing and floodplain flow processes. Since they are often simulated by different models, coupling these models may be a viable way to increase the integration of different physical drivers of simulated inundation estimates. To facilitate coupling different models and integrating across flood hazard processes, we here present GLOFRIM 2.0, a globally applicable framework for integrated hydrologichydrodynamic modelling. We then tested the hypothesis that smart model coupling can advance inundation modelling in the Amazon and Ganges basins. By means of GLOFRIM, we coupled the global hydrologic model PCR-GLOBWB with the hydrodynamic models CaMa-Flood and LISFLOOD-FP. Results show that replacing the kinematic wave approximation of the hydrologic model with the local inertia equation of CaMa-Flood greatly enhances accuracy of peak discharge simulations as expressed by an increase in the Nash-Sutcliffe efficiency (NSE) from 0.48 to 0.71. Flood maps obtained with LISFLOOD-FP improved representation of observed flood extent (critical success index C = 0.46), compared to downscaled products of PCR-GLOBWB and CaMa-Flood (C = 0.30 and C = 0.25, respectively). Results confirm that model coupling can indeed be a viable way forward towards more integrated flood simulations. However, results also sug-gest that the accuracy of coupled models still largely depends on the model forcing. Hence, further efforts must be undertaken to improve the magnitude and timing of simulated runoff. In addition, flood risk is, particularly in delta areas, driven by coastal processes. A more holistic representation of flood processes in delta areas, for example by incorporating a tide and surge model, must therefore be a next development step of GLOFRIM, making even more physically robust estimates possible for adequate flood risk management practices.

No abstract