Quantifying the relative impact of hydrological and hydraulic modelling parameterizations on uncertainty of inundation maps

ABSTRACT

Flood risk management strongly relies on inundation models for river basin zoning in flood-prone and risk-free areas. Floodplain zoning is significantly affected by the diverse and concurrent uncertainties that characterize the modelling chain used for producing inundation maps. In order to quantify the relative impact of the uncertainties linked to a lumped hydrological (rainfall–runoff) model and a FLO-2D hydraulic model, a Monte Carlo procedure is proposed in this work. The hydrological uncertainty is associated with the design rainfall estimation method, while the hydraulic model uncertainty is associated with roughness parameterization. This uncertainty analysis is tested on the case study of the Marta coastal catchment in Italy, by comparing the different frequency, extent and depth of inundation simulations associated with varying rainfall forcing and/or hydraulic model roughness realizations. The results suggest a significant predominance of the hydrological uncertainty with respect to the hydraulic one on the overall uncertainty associated with the simulated inundation maps.     

Application of Generalized Likelihood Uncertainty Estimation (GLUE) at different temporal scales to reduce the uncertainty level in modelled river flows

ABSTRACT

In this study, the distributed catchment-scale model, DiCaSM, was applied on five catchments across the UK. Given its importance, river flow was selected to study the uncertainty in streamflow prediction using the Generalized Likelihood Uncertainty Estimation (GLUE) methodology at different timescales (daily, monthly, seasonal and annual). The uncertainty analysis showed that the observed river flows were within the predicted bounds/envelope of 5% and 95% percentiles. These predicted river flow bounds contained most of the observed river flows, as expressed by the high containment ratio, CR. In addition to CR, other uncertainty indices – bandwidth B, relative bandwidth RB, degrees of asymmetry S and T, deviation amplitude D, relative deviation amplitude RD and the R factor – also indicated that the predicted river flows have acceptable uncertainty levels. The results show lower uncertainty in predicted river flows when increasing the timescale from daily to monthly to seasonal, with the lowest uncertainty associated with annual flows.     

Water balance components estimation under scenarios of land cover change in the Vea catchment,West Africa

ABSTRACT

The need for a detailed investigation of the Vea catchment water balance components cannot be overemphasized due to its accelerated land-cover dynamics and the associated impacts on the hydrological processes. This study assessed the possible consequences of land-use change scenarios (i.e. business as usual, BAU, and afforestation for the year 2025) compared to the 2016 baseline on the Vea catchment’s water balance components using the Soil and Water Assessment Tool (SWAT) model. The data used include daily climate and discharge, soil and land use/land cover maps. The results indicate that the mean annual water yield may increase by 9.1% under the BAU scenario but decrease by 2.7% under the afforestation scenario; actual evapotranspiration would decrease under BAU but increase under afforestation; and groundwater recharge may increase under both scenarios but would be more pronounced under the afforestation scenario. These outcomes highlight the significance of land-cover dynamics in water resource management and planning at the catchment.     

Detecting patterns of climate change in long-term forecasts of marine environmental parameters

ABSTRACT

Forecasting environmental parameters in the distant future requires complex modelling and large computational resources. Due to the sensitivity and complexity of forecast models, long-term parameter forecasts (e.g. up to 2100) are uncommon and only produced by a few organisations, in heterogeneous formats and based on different assumptions of greenhouse gases emissions. However, data mining techniques can be used to coerce the data to a uniform time and spatial representation, which facilitates their use in many applications. In this paper, streams of big data coming from AquaMaps and NASA collections of 126 long-term forecasts of nine types of environmental parameters are processed through a cloud computing platform in order to (i) standardise and harmonise the data representations, (ii) produce intermediate scenarios and new informative parameters, and (iii) align all sets on a common time and spatial resolution. Time series cross-correlation applied to these aligned datasets reveals patterns of climate change and similarities between parameter trends in 10 marine areas. Our results highlight that (i) the Mediterranean Sea may have a standalone ‘response’ to climate change with respect to other areas, (ii) the Poles are most representative of global forecasted change, and (iii) the trends are generally alarming for most oceans.     

On the Integrability of Limited-Area Numerical Weather Prediction Model ALADIN over Extended Time Periods

This note summarizes results of the first integration of regional numerical weather prediction model ALADIN in a climate mode. The ALADIN model, developed in an international cooperation led by Météo France, is operationally used for weather prediction. The grid step of the model is 12 km; the integration domain covers a major part of Europe. A one-month-long run has been performed with this model on observed boundary conditions (represented by assimilations by the global model ARPEGE). It is demonstrated that no excessive error is generated and accumulated in the model during the integration; hence the model is integrable for extended time periods and may serve a basis for a development towards a regional climate model.     

Stochastic estimation of facies using ground penetrating radar data

Explicitly defining large-scale heterogeneity is a necessary step of groundwater model calibration if accurate estimates of flow and transport are to be made. In this work, neural networks are used to estimate radar facies probabilities from ground penetrating radar (GPR) images, yielding stochastic facies-based models that honour the large-scale architecture of the subsurface. For synthetic GPR images, a neural network was able to correctly identify radar facies with an accuracy of approximately 90%. Manual interpretation of a set of 450 MHz GPR field data from the Borden aquifer resulted in the identification of four radar facies. Of these, a neural network was able to identify two facies with an accuracy of near 80% and one with an accuracy of 44%. The neural network was not able to identify the fourth facies, likely due to the choice of defining facies characteristics. Sequential indicator simulation was used to generate facies realizations conditioned to the radar facies probabilities. Numerical simulations indicate that significant improvements in the prediction of solute transport are possible when GPR is used to constrain the facies model compared to using well data alone, especially when data are sparse.This work was supported by funding to R. Knight under Grant No. DE-FG07–00ER15118-A000, Environmental Management Science Program, Office of Science and Technology, Office of Environment Management, United States Department of Energy (DOE). However, any opinions, findings, conclusions, or recommendations expressed herein are those of the authors and do not necessarily reflect the views of DOE. Further support was provided by a Stanford Graduate Fellowship to S. Moysey. The authors would also like to thank James Irving for his assistance with processing of the radar data.     

Groundwater flow and solute transport at the Mourquong saline-water disposal basin, Murray Basin, southeastern Australia

Saline groundwater and drainage effluent from irrigation are commonly stored in some 200 natural and artificial saline-water disposal basins throughout the Murray-Darling Basin of Australia. Their impact on underlying aquifers and the River Murray, one of Australia's major water supplies, is of serious concern. In one such scheme, saline groundwater is pumped into Lake Mourquong, a natural groundwater discharge complex. The disposal basin is hydrodynamically restricted by low-permeability lacustrine clays, but there are vulnerable areas in the southeast where the clay is apparently missing. The extent of vertical and lateral leakage of basin brines and the processes controlling their migration are examined using (1) analyses of chloride and stable isotopes of water (²H/¹H and ¹⁸O/¹⁶O) to infer mixing between regional groundwater and lake water, and (2) the variable-density groundwater flow and solute-transport code SUTRA. Hydrochemical results indicate that evaporated disposal water has moved at least 100 m in an easterly direction and that there is negligible movement of brines in a southerly direction towards the River Murray. The model is used to consider various management scenarios. Salt-load movement to the River Murray was highest in a "worst-case" scenario with irrigation employed between the basin and the River Murray. Present-day operating conditions lead to little, if any, direct movement of brine from the basin into the river. Electronic Publication     

Direct simulation of solute recycling in irrigated areas

Solute recycling from irrigation can be described as the process that occurs when the salt load that is extracted from irrigation wells and distributed on the fields is returned to the groundwater below irrigated surfaces by deep percolation. Unless the salt load leaves the system by means of drains or surface runoff, transfer to the groundwater will take place, sooner or later. This can lead to solute accumulation and thus to groundwater degradation, particularly in areas where extraction rates exceed infiltration rates (semi-arid and arid regions). Thus, considerable errors can occur in a predictive solute mass budget if the recycling process is not accounted for in the calculation. A method is proposed which allows direct simulation of solute recycling. The transient solute response at an extraction well is shown to be a superposition of solute mass flux contributions from n recycling cycles and is described as a function of the travel time distribution between a recycling point and a well. This leads to an expression for a transient ‘recycling source’ term in the advection–dispersion equation, which generates the effect of solute recycling. At long times, the ‘recycling source’ is a function of the local capture probability of the irrigation well and the solute mass flux captured by the well from the boundaries. The predicted concentration distribution at steady state reflects the maximum spatial concentration distribution in response to solute recycling and can thus be considered as the solute recycling potential or vulnerability of the entire domain for a given hydraulic setting and exploitation scheme. Simulation of the solute recycling potential is computationally undemanding and can therefore, for instance, be used for optimisation purposes. Also, the proposed method allows transient simulation of solute recycling with any standard flow and transport code.     

Pre‐ and post‐test mathematical modelling of a plan‐asymmetric reinforced concrete frame building

Pre‐ and post‐test analyses of the structural response of a three‐storey asymmetric reinforced concrete frame building were performed, aimed at supporting test preparation and performance as well as studying mathematical modelling. The building was designed for gravity loads only. Full‐scale pseudo‐dynamic tests were performed in the ELSA laboratory in Ispra. In the paper the results of initial parametric studies, of the blind pre‐test predictions, and of the post‐test analysis are summarized. In all studies a simple mathematical model, with one‐component member models with concentrated plasticity was employed. The pre‐test analyses were performed using the CANNY program. After the test results became available, the mathematical model was improved using an approach based on a displacement‐controlled analysis. Basically, the same mathematical model was used as in pre‐test analyses, except that the values of some of the parameters were changed. The OpenSees program was employed. Fair agreement between the test and numerical results was obtained. The results prove that relatively simple mathematical models are able to adequately simulate the detailed seismic response of reinforced concrete frame structures to a known ground motion, provided that the input parameters are properly determined. Copyright © 2006 John Wiley & Sons, Ltd.