Towards the construction of representative regional hydro(geo)logical numerical models: Modelling the upper Danube basin as a starting point

Introduction: Pressure on groundwater resources is increasing rapidly by population growth and climate change effects. Thus, it is urgent to quantify their availability and determine their dynamics at a global scale to assess the impacts of climate change or anthropogenically induced pressure, and to support water management strategies. In this context, regional hydrogeological numerical models become essential to simulate the behavior of groundwater resources. However, the construction of global hydrogeological models faces a lot of challenges that affect their accuracy.Methods: In this work, using the German portion of the Upper Danube Basin (∼43,000 km2) we outline common challenges encountered in parameterizing a regional-scale groundwater model, and provide an innovative approach to efficiently tackle such challenges. The hydrogeological model of the Danube consists of the groundwater finite element code OpenGeoSys forced by the groundwater recharge of the surface hydrological model mHM.Results: The main novelties of the suggested approach are 1) the use of spectral analyses of the river baseflow and a steady state calibration taking as reference the topography to constraint the hydraulic parameters and facilitate the calibration process, and 2) the calibration of the hydraulic parameters for a transient state model by considering parameters derived from the piezometric head evolution.Discussion/conclusion: The results show that the proposed methodology is useful to build a reliable large-scale groundwater model. Finally, the suggested approach is compared with the standard one used by other authors for the construction of global models. The comparison shows that the proposed approach allows for obtaining more reliable results, especially in mountainous areas

Multimodel assessment of flood characteristics in four large river basins at global warming of 1.5, 2.0 and 3.0 K above the pre-industrial level

This study assesses the flood characteristics (timing, magnitude and frequency) in the pre-industrial and historical periods, and analyzes climate change impacts on floods at the warming levels of 1.5, 2.0 and 3.0 K above the pre-industrial level in four large river basins as required by the Paris agreement. Three well-established hydrological models (HMs) were forced with bias-corrected outputs from four global climate models (GCMs) for the pre-industrial, historical and future periods until 2100. The long pre-industrial and historical periods were subdivided into multiple 31-year subperiods to investigate the natural variability. The mean flood characteristics in the pre-industrial period were derived from the large ensemble based on all GCMs, HMs and 31-year subperiods, and compared to the ensemble means in the historical and future periods. In general, the variance of simulated flood characteristics is quite large in the pre-industrial and historical periods. Mostly GCMs and HMs contribute to the variance, especially for flood timing and magnitude, while the selection of 31-year subperiods is an important source of variance for flood frequency. The comparison between the ensemble means shows that there are already some changes in flood characteristics between the pre-industrial and historical periods. There is a clear shift towards earlier flooding for the Rhine (1.5 K scenario) and Upper Mississippi (3.0 K scenario). The flood magnitudes show a substantial increase in the Rhine and Upper Yellow only under the 3.0 K scenario. The floods are projected to occur more frequently in the Rhine under the 1.5 and 2.0 K scenarios, and less frequently in the Upper Mississippi under all scenarios

On the choice of calibration metrics for “high-flow” estimation using hydrologic models

Calibration is an essential step for improving the accuracy of simulations generated using hydrologic models. A key modeling decision is selecting the performance metric to be optimized It has been common to use squared error performance metrics, or normalized variants such as Nash-Sutcliffe efficiency (NSE), based on the idea that their squared-error nature will emphasize the estimates of high flows. However, we conclude that NSE-based model calibrations actually result in poor reproduction of high-flow events, such as the annual peak flows that are used for flood frequency estimation. Using three different types of performance metrics, we calibrate two hydrological models at a daily step, the Variable Infiltration Capacity (VIC) model and the mesoscale Hydrologic Model (mHM), and evaluate their ability to simulate high-flow events for 492 basins throughout the contiguous United States. The metrics investigated are (1) NSE, (2) Kling-Gupta efficiency (KGE) and its variants, and (3) annual peak flow bias (APFB), where the latter is an application-specific metric that focuses on annual peak flows. As expected, the APFB metric produces the best annual peak flow estimates; however, performance on other high-flow-related metrics is poor. In contrast, the use of NSE results in annual peak flow estimates that are more than 20 % worse, primarily due to the tendency of NSE to underestimate observed flow variability. On the other hand, the use of KGE results in annual peak flow estimates that are better than from NSE, owing to improved flow time series metrics (mean and variance), with only a slight degradation in performance with respect to other related metrics, particularly when a non-standard weighting of the components of KGE is used. Stochastically generated ensemble simulations based on model residuals show the ability to improve the high-flow metrics, regardless of the deterministic performances. However, we emphasize that improving the fidelity of streamflow dynamics from deterministically calibrated models is still important, as it may improve high-flow metrics (for the right reasons). Overall, this work highlights the need for a deeper understanding of performance metric behavior and design in relation to the desired goals of model calibration.Open access journalThis item from the UA Faculty Publications collection is made available by the University of Arizona with support from the University of Arizona Libraries. If you have questions, please contact us at [email protected]

Shifts in flood generation processes exacerbate regional flood anomalies in Europe

Anomalies in the frequency of river floods, i.e., flood-rich or -poor periods, cause biases in flood risk estimates and thus make climate adaptation measures less efficient. While observations have recently confirmed the presence of flood anomalies in Europe, their exact causes are not clear. Here we analyse streamflow and climate observations during 1960-2010 to show that shifts in flood generation processes contribute more to the occurrence of regional flood anomalies than changes in extreme rainfall. A shift from rain on dry soil to rain on wet soil events by 5% increased the frequency of flood-rich periods in the Atlantic region, and an opposite shift in the Mediterranean region increased the frequency of flood-poor periods, but will likely make singular extreme floods occur more often. Flood anomalies driven by changing flood generation processes in Europe may further intensify in a warming climate and should be considered in flood estimation and management.publishedVersio

Improving global hydrological simulations through bias-correction and multi-model blending

There is an immediate need to develop accurate and reliable global hydrological forecasts in light of the future vulnerability to hydrological hazards and water scarcity under a changing climate. As a part of the World Meteorological Organization's (WMO) Global Hydrological Status and Outlook System (HydroSOS) initiative, we investigated different approaches to blending multi-model simulations for developing holistic operational global forecasts. The ULYSSES (mULti-model hYdrological SeaSonal prEdictionS system) dataset, to be published as “Global seasonal forecasts and reforecasts of river discharge and related hydrological variables ensemble from four state-of-the-art land surface and hydrological models” is used in this study. The first step for improving these forecasts is to investigate ways to improve the model simulations, as global models are not calibrated for local conditions. The analysis was performed over 119 different catchments worldwide for the baseline period of 1981–2019 for three variables: evapotranspiration, surface soil moisture and streamflow. This study evaluated blending approaches with a performance metric based (weighted) averaging of the multi-model simulations, using the catchment's Kling-Gupta Efficiency (KGE) for the variable to define the weight. Hydrological model simulations were also bias-corrected to improve the multi-model blending output. Weighted blending in conjunction with bias-correction provided the best improvement in performance for the catchments investigated. Applying modelled weights during blending original simulations improved performance over ungauged catchments. The results indicate that there is potential to successfully and easily implement the bias-corrected weighted blending approach to improve operational forecasts globally. This work can be used to improve water resources management and hydrological hazard mitigation, especially in data-sparse regions

High-resolution drought simulations and comparison to soil moisture observations in Germany

Germany\u27s 2018–2020 consecutive drought events resulted in multiple sectors – including agriculture, forestry, water management, energy production, and transport – being impacted. High-resolution information systems are key to preparedness for such extreme drought events. This study evaluates the new setup of the one-kilometer German drought monitor (GDM), which is based on daily soil moisture (SM) simulations from the mesoscale hydrological model (mHM). The simulated SM is compared against a set of diverse observations from single profile measurements, spatially distributed sensor networks, cosmic-ray neutron stations, and lysimeters at 40 sites in Germany. Our results show that the agreement of simulated and observed SM dynamics in the upper soil (0–25 cm) are especially high in the vegetative active period (0.84 median correlation R) and lower in winter (0.59 median R). The lower agreement in winter results from methodological uncertainties in both simulations and observations. Moderate but significant improvements between the coarser 4 km resolution setup and the ≈ 1.2 km resolution GDM in the agreement to observed SM dynamics is observed in autumn (+0.07 median R) and winter (+0.12 median R). Both model setups display similar correlations to observations in the dry anomaly spectrum, with higher overall agreement of simulations to observations with a larger spatial footprint. The higher resolution of the second GDM version allows for a more detailed representation of the spatial variability of SM, which is particularly beneficial for local risk assessments. Furthermore, the results underline that nationwide drought information systems depend both on appropriate simulations of the water cycle and a broad, high-quality, observational soil moisture database

Great Lakes Runoff Intercomparison Project Phase 3: Lake Erie (GRIP-E)

Hydrologic model intercomparison studies help to evaluate the agility of models to simulate variables such as streamflow, evaporation, and soil moisture. This study is the third in a sequence of the Great Lakes Runoff Intercomparison Projects. The densely populated Lake Erie watershed studied here is an important international lake that has experienced recent flooding and shoreline erosion alongside excessive nutrient loads that have contributed to lake eutrophication. Understanding the sources and pathways of flows is critical to solve the complex issues facing this watershed. Seventeen hydrologic and land-surface models of different complexity are set up over this domain using the same meteorological forcings, and their simulated streamflows at 46 calibration and seven independent validation stations are compared. Results show that: (1) the good performance of Machine Learning models during calibration decreases significantly in validation due to the limited amount of training data; (2) models calibrated at individual stations perform equally well in validation; and (3) most distributed models calibrated over the entire domain have problems in simulating urban areas but outperform the other models in validation

Understanding each other's models: a standard representation of global water models to support improvement, intercomparison, and communication

Global water models (GWMs) simulate the terrestrial water cycle, on the global scale, and are used to assess the impacts of climate change on freshwater systems. GWMs are developed within different modeling frameworks and consider different underlying hydrological processes, leading to varied model structures. Furthermore, the equations used to describe various processes take different forms and are generally accessible only from within the individual model codes. These factors have hindered a holistic and detailed understanding of how different models operate, yet such an understanding is crucial for explaining the results of model evaluation studies, understanding inter-model differences in their simulations, and identifying areas for future model development. This study provides a comprehensive overview of how state-of-the-art GWMs are designed. We analyze water storage compartments, water flows, and human water use sectors included in 16 GWMs that provide simulations for the Inter-Sectoral Impact Model Intercomparison Project phase 2b (ISIMIP2b). We develop a standard writing style for the model equations to further enhance model improvement, intercomparison, and communication. In this study, WaterGAP2 used the highest number of water storage compartments, 11, and CWatM used 10 compartments. Seven models used six compartments, while three models (JULES-W1, Mac-PDM.20, and VIC) used the lowest number, three compartments. WaterGAP2 simulates five human water use sectors, while four models (CLM4.5, CLM5.0, LPJmL, and MPI-HM) simulate only water used by humans for the irrigation sector. We conclude that even though hydrologic processes are often based on similar equations, in the end, these equations have been adjusted or have used different values for specific parameters or specific variables. Our results highlight that the predictive uncertainty of GWMs can be reduced through improvements of the existing hydrologic processes, implementation of new processes in the models, and high-quality input data

Rhine flood stories: Spatio‐temporal analysis of historic and projected flood genesis in the Rhine River basin

AbstractThe genesis of floods in large river basins often is complex. Streamflow originating from precipitation and snowmelt and different tributaries can superimpose and cause high water levels, threatening cities and communities along the riverbanks. For better understanding the mechanisms (origin and composition) of flood events in large and complex basins, we capture and share the story behind major historic and projected streamflow peaks in the Rhine River basin. Our analysis is based on hydrological simulations with the mesoscale Hydrological Model forced with both meteorological observations and an ensemble of climate projections. The spatio‐temporal analysis of the flood events includes the assessment and mapping of antecedent liquid precipitation, snow cover changes, generated and routed runoff, areal extents of events, and the above‐average runoff from major sub‐basins up to 10 days before a streamflow peak. We introduce and assess the analytical setup by presenting the flood genesis of the two well‐known Rhine floods that occurred in January 1995 and May 1999. We share our extensive collection of event‐based Rhine River flood genesis, which can be used in‐ and outside the scientific community to explore the complexity and diversity of historic and projected flood genesis in the Rhine basin. An interactive web‐based viewer provides easy access to all major historic and projected streamflow peaks at four locations along the Rhine. The comparison of peak flow genesis depending on different warming levels elucidates the role of changes in snow cover and precipitation characteristics in the (pre‐)Alps for flood hazards along the entire channel of the Rhine. Furthermore, our results suggest a positive correlation between flood magnitudes and areal extents of an event. Further hydro‐climatological research is required to improve the understanding of the climatic impact on the Rhine and beyond.The genesis of riverine floods in large river basins often is complex. Streamflow originating from precipitation and snowmelt and different tributaries can superimpose and cause high water levels threatening cities and communities along the riverbanks. In this study, we capture and share the story behind major historic and projected streamflow peaks in the large and complex basin of the Rhine River. https://doi.org/10.5281/zenodo.3239055https://github.com/ERottler/rhine-flood-genesishttp://natriskchange.ad.umwelt.uni-potsdam.de:3838/rhine-flood-genesishttps://b2share.eudat.eu/records/72d7a4f5d38043d1a137228b39c7ecc

Making the most out of a hydrological model data set: sensitivity analyses to open the model black-box

In this work, we investigate methods for gaining greater insight from hydrological model runs conducted for uncertainty quantification and model differentiation. We frame the sensitivity analysis questions in terms of the main purposes of sensitivity analysis: parameter prioritization, trend identification, and interaction quantification. For parameter prioritization, we consider variance-based sensitivity measures, sensitivity indices based on the L1-norm, the Kuiper metric, and the sensitivity indices of the DELSA methods. For trend identification, we investigate insights derived from graphing the one-way ANOVA sensitivity functions, the recently introduced CUSUNORO plots, and derivative scatterplots. For interaction quantification, we consider information delivered by variance-based sensitivity indices. We rely on the so-called given-data principle, in which results from a set of model runs are used to perform a defined set of analyses. One avoids using specific designs for each insight, thus controlling the computational burden. The methodology is applied to a hydrological model of a river in Belgium simulated using the well-established Framework for Understanding Structural Errors (FUSE) on five alternative configurations. The findings show that the integration of the chosen methods provides insights unavailable in most other analyses