scholarly journals Unraveling the time-dependent relevance of input model uncertainties for a lumped hydrologic model of a pre-alpine karst system

Author(s):  
Daniel Bittner ◽  
Beatrice Richieri ◽  
Gabriele Chiogna

AbstractUncertainties in hydrologic model outputs can arise for many reasons such as structural, parametric and input uncertainty. Identification of the sources of uncertainties and the quantification of their impacts on model results are important to appropriately reproduce hydrodynamic processes in karst aquifers and to support decision-making. The present study investigates the time-dependent relevance of model input uncertainties, defined as the conceptual uncertainties affecting the representation and parameterization of processes relevant for groundwater recharge, i.e. interception, evapotranspiration and snow dynamic, on the lumped karst model LuKARS. A total of nine different models are applied, three to compute interception (DVWK, Gash and Liu), three to compute evapotranspiration (Thornthwaite, Hamon and Oudin) and three to compute snow processes (Martinec, Girons Lopez and Magnusson). All the input model combinations are tested for the case study of the Kerschbaum spring in Austria. The model parameters are kept constant for all combinations. While parametric uncertainties computed for the same model in previous studies do not show pronounced temporal variations, the results of the present work show that input uncertainties are seasonally varying. Moreover, the input uncertainties of evapotranspiration and snowmelt are higher than the interception uncertainties. The results show that the importance of a specific process for groundwater recharge can be estimated from the respective input uncertainties. These findings have practical implications as they can guide researchers to obtain relevant field data to improve the representation of different processes in lumped parameter models and to support model calibration.

2005 ◽  
Vol 5 ◽  
pp. 31-35 ◽  
Author(s):  
J. Götzinger ◽  
A. Bárdossy

Abstract. Water balance models provide significant input to integrated models that are used to simulate river basin processes. However, one of the primary problems involves the coupling and simultaneous calibration of rainfall-runoff and groundwater models. This problem manifests itself through circular arguments - the hydrologic model is modified to calculate highly discretized groundwater recharge rates as input to the groundwater model which provides modeled base flow for the flood-routing module of the rainfall-runoff model. A possibility to overcome this problem using a modified version of the HBV Model is presented in this paper. Regionalisation and optimization methods lead to objective and efficient calibration despite large numbers of parameters. The representation of model parameters by transfer functions of catchment characteristics enables consistent parameter estimation. By establishing such relationships, models are calibrated for the parameters of the transfer functions instead of the model parameters themselves. Simulated annealing, using weighted Nash-Sutcliffe-coefficients of variable temporal aggregation, assists in efficient parameterisations. The simulations are compared to observed discharge and groundwater recharge modeled by the State Institute for Environmental Protection Baden-Württemberg using the model TRAIN-GWN.


Author(s):  
Ashley E. Van Beusekom ◽  
Lauren E. Hay ◽  
Andrew R. Bennett ◽  
Young-Don Choi ◽  
Martyn P. Clark ◽  
...  

Abstract Surface meteorological analyses are an essential input (termed ‘forcing’) for hydrologic modeling. This study investigated the sensitivity of different hydrologic model configurations to temporal variations of seven forcing variables (precipitation rate, air temperature, longwave radiation, specific humidity, shortwave radiation, wind speed, and air pressure). Specifically, the effects of temporally aggregating hourly forcings to hourly daily-average forcings were examined. The analysis was based on 14 hydrological outputs from the Structure for Unifying Multiple Modeling Alternatives (SUMMA) model for the 671 Catchment Attributes and MEteorology for Large-sample Studies (CAMELS) basins across the contiguous United States (CONUS). Results demonstrated that the hydrologic model sensitivity to temporally aggregating the forcing inputs varies across model output variables and model locations. We used Latin Hypercube sampling to sample model parameters from eight combinations of three influential model physics choices (three model decisions with two options for each decision, i.e., eight model configurations). Results showed that the choice of model physics can change the relative influence of forcing on model outputs and the forcing importance may not be dependent on the parameter space. This allows for model output sensitivity to forcing aggregation to be tested prior to parameter calibration. More generally, this work provides a comprehensive analysis of the dependence of modeled outcomes on input forcing behavior, providing insight into the regional variability of forcing variable dominance on modeled outputs across CONUS.


Author(s):  
Christopher J. Arthurs ◽  
Nan Xiao ◽  
Philippe Moireau ◽  
Tobias Schaeffter ◽  
C. Alberto Figueroa

AbstractA major challenge in constructing three dimensional patient specific hemodynamic models is the calibration of model parameters to match patient data on flow, pressure, wall motion, etc. acquired in the clinic. Current workflows are manual and time-consuming. This work presents a flexible computational framework for model parameter estimation in cardiovascular flows that relies on the following fundamental contributions. (i) A Reduced-Order Unscented Kalman Filter (ROUKF) model for data assimilation for wall material and simple lumped parameter network (LPN) boundary condition model parameters. (ii) A constrained least squares augmentation (ROUKF-CLS) for more complex LPNs. (iii) A “Netlist” implementation, supporting easy filtering of parameters in such complex LPNs. The ROUKF algorithm is demonstrated using non-invasive patient-specific data on anatomy, flow and pressure from a healthy volunteer. The ROUKF-CLS algorithm is demonstrated using synthetic data on a coronary LPN. The methods described in this paper have been implemented as part of the CRIMSON hemodynamics software package.


Author(s):  
Pal Kostka ◽  
Zsolt Techy ◽  
James J. Sienicki

Hydrogen combustion may represent a threat to containment integrity in a VVER-440/213 plant owing to the combination of high pressure and high temperature. A study has been carried out using the GASFLOW 2.1 three-dimensional CFD code to evaluate the hydrogen distribution in the containment during a beyond design basis accident. The VVER-440/213 containment input model consists of two 3D blocks connected via one-dimensional (1D) ducts. One 3D block contains the reactor building and the accident localization tower with the suppression pools. Another 3D block models the air traps. 1D ducts represent the check valves connecting the accident localization tower with the air traps. The VVER pressure suppression system, called “bubbler condenser,” was modeled as a distributed heat sink with water thermodynamic properties. This model accounts for the energy balance. However, it is not currently possible to model dynamic phenomena associated with the water pools (e.g., vent clearing, level change). The GASFLOW 2.1 calculation gave detailed results for the spatial distribution of thermal-hydraulic parameters and gas concentrations. The range and trend of the parameters are reasonable and valuable. There are particularly interesting circulation patterns around the steam generators, in the bubbler tower and other primary system compartments. In case of the bubbler tower, concentration and temperature contour plots show an inhomogeneous distribution along the height and width, changing during the accident. Hydrogen concentrations also vary within primary system compartments displaying lower as well as higher (up to 13–20% and higher) values in some nodes. Prediction of such concentration distributions was not previously possible with lumped parameter codes. GASFLOW 2.1 calculations were compared with CONTAIN 1.2 (lumped parameter code) results. Apart from the qualitatively similar trends, there are, for the time being, quantitative differences between the results concerning, for example, pressure histories, or the total amount of steam available in the containment. The results confirm the importance of detailed modeling of the containment, as well as of the bubbler condenser and sump water pools. The study showed that modeling of hydrogen distribution in the VVER-440/213 containment was possible using the GASFLOW 2.1 code with reasonable results and remarkable physical insights.


2010 ◽  
Vol 11 (3) ◽  
pp. 781-796 ◽  
Author(s):  
Jonathan J. Gourley ◽  
Scott E. Giangrande ◽  
Yang Hong ◽  
Zachary L. Flamig ◽  
Terry Schuur ◽  
...  

Abstract Rainfall estimated from the polarimetric prototype of the Weather Surveillance Radar-1988 Doppler [WSR-88D (KOUN)] was evaluated using a dense Micronet rain gauge network for nine events on the Ft. Cobb research watershed in Oklahoma. The operation of KOUN and its upgrade to dual polarization was completed by the National Severe Storms Laboratory. Storm events included an extreme rainfall case from Tropical Storm Erin that had a 100-yr return interval. Comparisons with collocated Micronet rain gauge measurements indicated all six rainfall algorithms that used polarimetric observations had lower root-mean-squared errors and higher Pearson correlation coefficients than the conventional algorithm that used reflectivity factor alone when considering all events combined. The reflectivity based relation R(Z) was the least biased with an event-combined normalized bias of −9%. The bias for R(Z), however, was found to vary significantly from case to case and as a function of rainfall intensity. This variability was attributed to different drop size distributions (DSDs) and the presence of hail. The synthetic polarimetric algorithm R(syn) had a large normalized bias of −31%, but this bias was found to be stationary. To evaluate whether polarimetric radar observations improve discharge simulation, recent advances in Markov Chain Monte Carlo simulation using the Hydrology Laboratory Research Distributed Hydrologic Model (HL-RDHM) were used. This Bayesian approach infers the posterior probability density function of model parameters and output predictions, which allows us to quantify HL-RDHM uncertainty. Hydrologic simulations were compared to observed streamflow and also to simulations forced by rain gauge inputs. The hydrologic evaluation indicated that all polarimetric rainfall estimators outperformed the conventional R(Z) algorithm, but only after their long-term biases were identified and corrected.


2012 ◽  
Vol 16 (9) ◽  
pp. 3083-3099 ◽  
Author(s):  
H. Xie ◽  
L. Longuevergne ◽  
C. Ringler ◽  
B. R. Scanlon

Abstract. Irrigation development is rapidly expanding in mostly rainfed Sub-Saharan Africa. This expansion underscores the need for a more comprehensive understanding of water resources beyond surface water. Gravity Recovery and Climate Experiment (GRACE) satellites provide valuable information on spatio-temporal variability in water storage. The objective of this study was to calibrate and evaluate a semi-distributed regional-scale hydrologic model based on the Soil and Water Assessment Tool (SWAT) code for basins in Sub-Saharan Africa using seven-year (July 2002–April 2009) 10-day GRACE data and multi-site river discharge data. The analysis was conducted in a multi-criteria framework. In spite of the uncertainty arising from the tradeoff in optimising model parameters with respect to two non-commensurable criteria defined for two fluxes, SWAT was found to perform well in simulating total water storage variability in most areas of Sub-Saharan Africa, which have semi-arid and sub-humid climates, and that among various water storages represented in SWAT, water storage variations in soil, vadose zone and groundwater are dominant. The study also showed that the simulated total water storage variations tend to have less agreement with GRACE data in arid and equatorial humid regions, and model-based partitioning of total water storage variations into different water storage compartments may be highly uncertain. Thus, future work will be needed for model enhancement in these areas with inferior model fit and for uncertainty reduction in component-wise estimation of water storage variations.


2017 ◽  
Author(s):  
Miao Jing ◽  
Falk Heße ◽  
Wenqing Wang ◽  
Thomas Fischer ◽  
Marc Walther ◽  
...  

Abstract. Most of the current large scale hydrological models do not contain a physically-based groundwater flow component. The main difficulties in large-scale groundwater modeling include the efficient representation of unsaturated zone flow, the characterization of dynamic groundwater-surface water interaction and the numerical stability while preserving complex physical processes and high resolution. To address these problems, we propose a highly-scalable coupled hydrologic and groundwater model (mHM#OGS) based on the integration of two open-source modeling codes: the mesoscale hydrologic Model (mHM) and the finite element simulator OpenGeoSys (OGS). mHM#OGS is coupled using a boundary condition-based coupling scheme that dynamically links the surface and subsurface parts. Nested time stepping allows smaller time steps for typically faster surface runoff routing in mHM and larger time steps for slower subsurface flow in OGS. mHM#OGS features the coupling interface which can transfer the groundwater recharge and river baseflow rate between mHM and OpenGeoSys. Verification of the coupled model was conducted using the time-series of observed streamflow and groundwater levels. Moreover, we force the transient model using groundwater recharge in two scenarios: (1) spatially variable recharge based on the mHM simulations, and (2) spatially homogeneous groundwater recharge. The modeling result in first scenario has a slightly higher correlation with groundwater head time-series, which further validates the plausibility of spatial groundwater recharge distribution calculated by mHM in the mesocale. The statistical analysis of model predictions shows a promising prediction ability of the model. The offline coupling method implemented here can reproduce reasonable groundwater head time series while keep a desired level of detail in the subsurface model structure with little surplus in computational cost. Our exemplary calculations show that the coupled model mHM#OGS can be a valuable tool to assess the effects of variability in land surface heterogeneity, meteorological, topographical forces and geological zonation on the groundwater flow dynamics.


2016 ◽  
Vol 9 (3) ◽  
Author(s):  
Ata Joodavi ◽  
Mohammad Zare ◽  
Ezzat Raeisi ◽  
Mohammad Bagher Ahmadi

Author(s):  
Mehmet Cüneyd Demirel ◽  
Julian Koch ◽  
Gorka Mendiguren ◽  
Simon Stisen

Hydrologic models are conventionally constrained and evaluated using point measurements of streamflow, which represents an aggregated catchment measure. As a consequence of this single objective focus, model parametrization and model parameter sensitivity are typically not reflecting other aspects of catchment behavior. Specifically for distributed models, the spatial pattern aspect is often overlooked. Our paper examines the utility of multiple performance measures in a spatial sensitivity analysis framework to determine the key parameters governing the spatial variability of predicted actual evapotranspiration (AET). Latin hypercube one-at-a-time (LHS-OAT) sampling strategy with multiple initial parameter sets was applied using the mesoscale hydrologic model (mHM) and a total of 17 model parameters were identified as sensitive. The results indicate different parameter sensitivities for different performance measures focusing on temporal hydrograph dynamics and spatial variability of actual evapotranspiration. While spatial patterns were found to be sensitive to vegetation parameters, streamflow dynamics were sensitive to pedo-transfer function (PTF) parameters. Above all, our results show that behavioral model definition based only on streamflow metrics in the generalized likelihood uncertainty estimation (GLUE) type methods require reformulation by incorporating spatial patterns into the definition of threshold values to reveal robust hydrologic behavior in the analysis.


2021 ◽  
Vol 25 (11) ◽  
pp. 5805-5837
Author(s):  
Oscar M. Baez-Villanueva ◽  
Mauricio Zambrano-Bigiarini ◽  
Pablo A. Mendoza ◽  
Ian McNamara ◽  
Hylke E. Beck ◽  
...  

Abstract. Over the past decades, novel parameter regionalisation techniques have been developed to predict streamflow in data-scarce regions. In this paper, we examined how the choice of gridded daily precipitation (P) products affects the relative performance of three well-known parameter regionalisation techniques (spatial proximity, feature similarity, and parameter regression) over 100 near-natural catchments with diverse hydrological regimes across Chile. We set up and calibrated a conceptual semi-distributed HBV-like hydrological model (TUWmodel) for each catchment, using four P products (CR2MET, RF-MEP, ERA5, and MSWEPv2.8). We assessed the ability of these regionalisation techniques to transfer the parameters of a rainfall-runoff model, implementing a leave-one-out cross-validation procedure for each P product. Despite differences in the spatio-temporal distribution of P, all products provided good performance during calibration (median Kling–Gupta efficiencies (KGE′s) > 0.77), two independent verification periods (median KGE′s >0.70 and 0.61, for near-normal and dry conditions, respectively), and regionalisation (median KGE′s for the best method ranging from 0.56 to 0.63). We show how model calibration is able to compensate, to some extent, differences between P forcings by adjusting model parameters and thus the water balance components. Overall, feature similarity provided the best results, followed by spatial proximity, while parameter regression resulted in the worst performance, reinforcing the importance of transferring complete model parameter sets to ungauged catchments. Our results suggest that (i) merging P products and ground-based measurements does not necessarily translate into an improved hydrologic model performance; (ii) the spatial resolution of P products does not substantially affect the regionalisation performance; (iii) a P product that provides the best individual model performance during calibration and verification does not necessarily yield the best performance in terms of parameter regionalisation; and (iv) the model parameters and the performance of regionalisation methods are affected by the hydrological regime, with the best results for spatial proximity and feature similarity obtained for rain-dominated catchments with a minor snowmelt component.


Sign in / Sign up

Export Citation Format

Share Document