Rainfall/runoff simulation with 2D full shallow water equations: Sensitivity analysis and calibration of infiltration parameters

2016 ◽  
Vol 536 ◽  
pp. 496-513 ◽  
Author(s):  
Javier Fernández-Pato ◽  
Daniel Caviedes-Voullième ◽  
Pilar García-Navarro
2018 ◽  
Vol 20 (4) ◽  
pp. 898-916 ◽  
Author(s):  
J. Fernández-Pato ◽  
J. L. Gracia ◽  
P. García-Navarro

Abstract In this work, a distributed two-dimensional (2D) shallow water (SW) flow model is combined with a fractional-order version of the Green-Ampt (FOGA) infiltration law to improve rainfall/runoff simulation in real catchments. The surface water model is based on a robust finite volume method on triangular grids that can handle flow over dry bed and multiple wet/dry fronts. When supplied with adequate infiltration laws, this model can provide useful information in surface hydrology. The classical Green-Ampt law is generalized by using a Caputo fractional derivative of order less than or equal to 1 in Darcy's law. The novelty of this combination is that, on the one hand, the distributed SW simulation provides a detailed surface water distribution and, on the other hand, the FOGA model offers the possibility to model infiltration rates not monotonically decreasing. In order to obtain the best results, a non-uniform order of the fractional derivative depending on the cumulative infiltration and the existence of available surface water is proposed for realistic cases. This allows significant improvement of previous published numerical results in the literature for several storm events in catchments where the infiltration process occurs.


2020 ◽  
Author(s):  
Pelagiya Belyakova ◽  
Ekaterina Vasil'eva ◽  
Andrey Aleksyuk ◽  
Vitaly Belikov ◽  
Boris Gartsman ◽  
...  

<p>In the Russian part of Western Caucasus heavy rainfall episodes frequently occur, leading to flash floods that often cause fatalities and severe damage. As soon as climate change is expected to increase the risk of flash floods it is necessary to improve flood forecasting and flood risk mapping as well as other precautionary measures. For this scope the better knowledge of catchment response on heavy precipitation is needed using rainfall-runoff simulation and further hydrodynamic modelling of inundation of urbanized areas.</p><p>There is a number of models used for flash flood simulation. In this study we used an available unit hydrograph model KW-GIUH [1] and a hydrodynamic model STREAM 2D CUDA [2]. KW-GIUH model only schematically describes overland flow over the catchment, nonlinear character of response is introduced via kinematic-wave approximation of the travel time. STREAM 2D CUDA is based on numerical solution of shallow water equations in a two-dimensional formulation according to the original algorithm using the exact solution of the Riemann problem [2], due to which the calculation is performed for the entire catchment without special allocation of the channel network. Models were tested on several flash flood events on the river Adagum (6-7 July 2012, catastrophic flood in the Krymsk town) and the Zapadny Dagomys river (25 June 2015, 24-25 October 2018, Sochi).</p><p>Comparison of simulation results was done as the same input data set was used. Input data included DEM HydroSHEDS, measured hourly precipitation and runoff volumes observed on gauges and estimated after high-water marks. Also 10-min water levels from a regional automated flood monitoring system of the Krasnodar Territory were applied. Simulated runoff volumes and peak timing were analyzed. For the Zapadny Dagomys river a forecasting calculation was done using precipitation forecast from COSMO-Ru. For the Adagum river STREAM 2D CUDA allowed to conduct an experiment to assess possible effect from potential reservoir-traps in the tributaries. The results of the rainfall-runoff simulation by the KW-GIUH model can be used as inflow to the boundary of the area for hydrodynamic modeling using STREAM 2D CUDA, also for operational use. Scenario calculations with changing hydraulic conditions at the catchment can be simulated using the STREAM 2D CUDA model itself.</p><p>The flood simulation was supported by the Russian Science Foundation under grant №17-77-30006. Data processing from an automated flood monitoring system in the Krasnodar Territory is funded by Russian Foundation for Basic Research and the Krasnodar Territory, grant № 19-45-233007.</p><p>References:</p><ol><li>Lee K.T., Cheng N.K., Gartsman B.I., Bugayets A.N. (2009): A current version of the model of a unit hydrograph and its use in Taiwan and Russia, Geography and Natural Resources, Volume 30, issue 1, pp. 79–85. https://doi.org/10.1016/j.gnr.2009.03.015</li> <li>Aleksyuk A.I., Belikov V.V. (2017): Simulation of shallow water flows with shoaling areas and bottom discontinuities, Computational Mathematics and Mathematical Physics, Volume 57, issue 2, pp. 318–339. https://doi.org/10.1134/S0965542517020026</li> </ol>


2015 ◽  
Vol 74 (11) ◽  
pp. 7295-7305 ◽  
Author(s):  
Asier Lacasta ◽  
Mario Morales-Hernández ◽  
Javier Murillo ◽  
Pilar García-Navarro

Author(s):  
Xiao-Hua Zhu ◽  
Xiao-Hua Zhu ◽  
Ze-Nan Zhu ◽  
Ze-Nan Zhu ◽  
Xinyu Guo ◽  
...  

A coastal acoustic tomography (CAT) experiment for mapping the tidal currents in the Zhitouyang Bay was successfully carried out with seven acoustic stations during July 12 to 13, 2009. The horizontal distributions of tidal current in the tomography domain are calculated by the inverse analysis in which the travel time differences for sound traveling reciprocally are used as data. Spatial mean amplitude ratios M2 : M4 : M6 are 1.00 : 0.15 : 0.11. The shallow-water equations are used to analyze the generation mechanisms of M4 and M6. In the deep area, velocity amplitudes of M4 measured by CAT agree well with those of M4 predicted by the advection terms in the shallow water equations, indicating that M4 in the deep area where water depths are larger than 60 m is predominantly generated by the advection terms. M6 measured by CAT and M6 predicted by the nonlinear quadratic bottom friction terms agree well in the area where water depths are less than 20 m, indicating that friction mechanisms are predominant for generating M6 in the shallow area. Dynamic analysis of the residual currents using the tidally averaged momentum equation shows that spatial mean values of the horizontal pressure gradient due to residual sea level and of the advection of residual currents together contribute about 75% of the spatial mean values of the advection by the tidal currents, indicating that residual currents in this bay are induced mainly by the nonlinear effects of tidal currents.


Sign in / Sign up

Export Citation Format

Share Document