An assessment of large-scale flood modelling based on LiDAR data

Large-scale flood modelling approaches designed for regional to continental scales usually rely on relatively simple assumptions to represent the potentially highly complex river bathymetry at the watershed scale based on digital elevation models (DEMs) with a resolution in the range of 25–30 m. Here, high-resolution (1 m) LiDAR DEMs are employed to present a novel large-scale methodology using a more realistic estimation of bathymetry based on hydrogeomorphological GIS tools to extract water surface slope. The large-scale 1D/2D flood model LISFLOOD-FP is applied to validate the simulated flood levels using detailed water level data in four different watersheds in Quebec (Canada), including continuous profiles over extensive distances measured with the HydroBall technology. A GIS-automated procedure allows to obtain the average width required to run LISFLOOD-FP. The GIS-automated procedure to estimate bathymetry from LiDAR water surface data uses a hydraulic inverse problem based on discharge at the time of acquisition of LiDAR data. A tiling approach, allowing several small independent hydraulic simulations to cover an entire watershed, greatly improves processing time to simulate large watersheds with a 10-m resampled LiDAR DEM. Results show significant improvements to large-scale flood modelling at the watershed scale with standard deviation in the range of 0.30 m and an average fit of around 90%. The main advantage of the proposed approach is to avoid the need to collect expensive bathymetry data to efficiently and accurately simulate flood levels over extensive areas.     

Influence of diversion angle on water and sediment flow into diversion channel

Accumulation of the sediment in the stream of the diversion channels adversely affects its operational systems. Diversion channels are often constructed perpendicular to the main river. In this study, the water flow and sediment transport in the diversion channel with different angles were investigated in an attempt to maximize water discharge and minimize sediment discharge. A physical model with movable bed was used to simulate water and sediment flow with five diversion angles (θ) between (30°–90°). Moreover, three bed width ratios (B_r) (the relation between diversion to main channel bed width) between 30% and 50% and five total discharges between (7.25 L/s to 12.25 L/s) were considered for each case of (θ). The results showed, up to 10%, increasing in proportion discharge ratios for 30 and 45 diversion angles compared with 90° diversion angle. The results also showed that the lowest diversion sediment concentration was provided by the (θ) of 30°. Across all scenarios, the average proportion concentration reduction was 64%, compared with 90° diversion angle. Closer observation of the diversion system mechanism confirmed that decreased (θ) result in decreased sediment concentrations in the diversion channels. In conclusion, the diversion channel water and sediment discharge could be effectively managed by changing the (θ) to 30° or 45° instead of 90°.     

Comparison of two modeling approaches for groundwater–surface water interactions

An assessment of interactions between groundwater and surface water was carried out by applying two different modeling approaches to a small‐scale study area in the municipality of Havelock, Quebec. The first approach involved a commonly used sequential procedure that consists in determining the daily recharge rate using a quasi 2D infiltration model (HELP), applied in the next step as an imposed flux to a 3D finite‐element groundwater flow model. The flow model was calibrated under steady‐state and transient conditions against measured water levels. The second approach was based on a recently developed physically based, 3D fully coupled groundwater–surface water flow model (CATHY) applied to the entire flow domain in an integrated manner. Implementation, calibration, and results of the simulations for both approaches are presented and discussed. For equal annual precipitation (1038 mm/y) and evapotranspiration (556 mm/y), the second approach computed a recharge rate of 233 mm/y (8.9% higher than the first approach) and a net upward flow from the fractured aquifer (the first approach predicted a net downward flow to the rock). The simulated annual discharge was similar for the two approaches (9.6% difference). Both approaches were found to be useful in understanding the interactions between groundwater and surface water, although limitations are apparent in the sequential procedure's inability to account for surface–subsurface feedbacks, for instance near stream reaches where groundwater discharge is prevalent. The decoupled, two‐model approach provides disaggregated surface, vadose, and aquifer flows, and a simple aperçu at the different components of total discharge. The fully coupled model accounts for continuous water exchanges between the land surface, subsurface, and stream channel in a more complex manner, and produces a better match against observed data. Copyright © 2012 John Wiley & Sons, Ltd.     

Aquifer Vulnerability Assessment Based on Sequence Stratigraphic and 39Ar Transport Modeling

A large‐scale groundwater flow and transport model is developed for a deep‐seated (100 to 300 m below ground surface) sedimentary aquifer system. The model is based on a three‐dimensional (3D) hydrostratigraphic model, building on a sequence stratigraphic approach. The flow model is calibrated against observations of hydraulic head and stream discharge while the credibility of the transport model is evaluated against measurements of ³⁹Ar from deep wells using alternative parameterizations of dispersivity and effective porosity. The directly simulated 3D mean age distributions and vertical fluxes are used to visualize the two‐dimensional (2D)/3D age and flux distribution along transects and at the top plane of individual aquifers. The simulation results are used to assess the vulnerability of the aquifer system that generally has been assumed to be protected by thick overlaying clayey units and therefore proposed as future reservoirs for drinking water supply. The results indicate that on a regional scale these deep‐seated aquifers are not as protected from modern surface water contamination as expected because significant leakage to the deeper aquifers occurs. The complex distribution of local and intermediate groundwater flow systems controlled by the distribution of the river network as well as the topographical variation (Tóth 1963) provides the possibility for modern water to be found in even the deepest aquifers.     

Two-dimensional numerical modeling of dam-break flows over natural terrain using a central explicit scheme

A two-dimensional (2D) numerical model has been developed to solve shallow water equations for simulation of dam-break flows. The spatial derivatives are discretized using a well-balanced explicit central upwind conservative scheme. The scheme is Riemann solver free and guarantees the positivity of the flow depth over complex topography if the Courant number is kept less than 0.25. The time integration is performed by Euler’s scheme. The model is verified against analytical results for water surface elevation and discharge for three benchmark test cases. A good agreement between analytical solutions and computed results is observed. The property of well-balancing in still water over an uneven bottom is also confirmed. The model is then validated by simulating a laboratory experiment in which a dam break flow propagates over a triangular obstacle. The model performance was found to be satisfactory. A dam break laboratory experimental test case on a frictionless horizontal bottom is also simulated for 2D validation of the model, and good agreement between simulation and the experimental data is observed. The suitability of the proposed model for real life applications is demonstrated by simulating the Malpasset dam-break event, which occurred in 1959 in France. The computed arrival time of the flood wave front and the maximum flow depths at various observation points matched well with the measurements on a 1/400 scale physical model. The overall performance indicates that this model can be applied for simulation of dam-break waves in real life cases.     

DEFORMATION OF RIVER BED AND THE CHARACTERISTICS OF SEDIMENT TRANSPORT DURING HYPER-CONCENTRATED FLOOD IN THE YELLOW RIVER

LINTRODUCTIONTheYellowforeriswellknownasaheavilysilt-caacingriverintheworld.Haaer-concentratedfloodsoftenoccurinitsmasterstemandaswellasthetriblltaries.ThecharacteristicsofdeformationandsedimellttranSportdifferfromreachtOreach.Duetohighsedimelltconcelltration,StrongfluvialactionandthenatUreofunSteadysedimenttranSPOrt,problemsandabnormalphenomenonareoflencreated.ThelaterbringinimpacttOnoodcontrolOfthelowerreach.Thus,itishelpfultoenhancetheunderstandingofthemotionlawsOfhaer-concentrat…     

The topographic data source of digital terrain models as a key element in the accuracy of hydraulic flood modelling

The effects of the topographic data source and resolution on the hydraulic modelling of floods were analysed. Seven digital terrain models (DTMs) were generated from three different altimetric sources: a global positioning system (GPS) survey and bathymetry; high‐resolution laser altimetry data LiDAR (light detection and ranging); and vectorial cartography (1:5000). Hydraulic results were obtained, using the HEC‐RAS one‐dimensional model, for all seven DTMs. The importance of the DTM's accuracy on the hydraulic modelling results was analysed within three different hydraulic contexts: (1) the discharge and water surface elevation results from the hydraulic model; (2) the delineation of the flooded area; and (3) the relative sensitivity of the hydraulic model to changes in the Manning's n roughness coefficient. The contour‐based DTM was the least accurate with a root mean square error (RMSE) of 4·5 m in the determination of the water level and a variation of up to 50 per cent in the estimation of the inundated area of the floodplain. The GPS‐based DTM produced more realistic water surface elevation results and variations of up to 8 per cent in terms of the flooded area. The laser‐based model's RMSE for water level was 0·3 m, with the flooded area varying by less than 1 per cent. The LiDAR data also showed the greatest sensitivity to changes in the Manning's roughness coefficient. An analysis of the effect of mesh resolution indicated an influence on the delineation of the flooded area with variations of up to 7·3 per cent. In addition to determining the accuracy of the hydraulic modelling results produced from each DTM, an analysis of the time–cost ratio of each topographic data source illustrates that airborne laser scanning is a cost‐effective means of developing a DTM of sufficient accuracy, especially over large areas. Copyright © 2006 John Wiley & Sons, Ltd.     

Study of the flood control scheduling scheme for the Three Gorges Reservoir in a catastrophic flood

The Three Gorges Project is the world's largest water conservancy project. According to the design standards for the 1,000‐year flood, flood diversion areas in the Jingjiang reach of the Yangtze River must be utilized to ensure the safety of the Jingjiang area and the city of Wuhan. However, once these areas are used, the economic and life loss in these areas may be very great. Therefore, it is vital to reduce this loss by developing a scheme that reduces the use of the flood diversion areas through flood regulation by the Three Gorges Reservoir (TGR), under the premise of ensuring the safety of the Three Gorges Dam. For a 1,000‐year flood on the basis of a highly destructive flood in 1954, this paper evaluates scheduling schemes in which flood diversion areas are or are not used. The schemes are simulated based on 2.5‐m resolution reservoir topography and an optimized model of dynamic capacity flood regulation. The simulation results show the following. (a) In accord with the normal flood‐control regulation discharge, the maximum water level above the dam should be not more than 175 m, which ensures the safety of the dam and reservoir area. However, it is necessary to utilize the flood diversion areas within the Jingjiang area, and flood discharge can reach 2.81 billion m³. (b) In the case of relying on the TGR to impound floodwaters independently rather than using the flood diversion areas, the maximum water level above the dam reaches 177.35 m, which is less than the flood check level of 180.4 m to ensure the safety of the Three Gorges Dam. The average increase of the TGR water level in the Chongqing area is not more than 0.11 m, which indicates no significant effect on the upstream reservoir area. Comparing the various scheduling schemes, when the flood diversion areas are not used, it is believed that the TGR can execute safe flood control for a 1,000‐year flood, thereby greatly reducing flood damage.     

Developing channel and floodplain dimensions with limited data: a case study in the Amazon Basin

This research builds on the concept of hydraulic geometry and presents a methodology for estimating bankfull discharge and the hydraulic geometry coefficients and exponents for a station using limited data; only stage‐discharge and Landsat imagery. The approach is implemented using 82 streamflow gauging locations in the Amazon Basin. Using the estimated values for the hydraulic geometry relations, bankfull discharge, discharge data above bankfull and upstream drainage area at each site, relationships for estimating channel and floodplain characteristics as a function of drainage area are developed. Specifically, this research provides relationships for estimating bankfull discharge, bankfull depth, bankfull width, and floodplain width as a function of upstream drainage area in the Amazon Basin intended for providing reasonable cross‐section estimates for large scale hydraulic routing models. The derived relationships are also combined with a high resolution drainage network to develop relationships for estimating cumulative upstream channel lengths and surface areas as a function of the specified minimum channel width ranging from 2 m to 1 km (i.e. threshold drainage areas ranging from 1 to 431,000 km²). At the finest resolution (i.e. all channels greater than 2 m or a threshold area of 1 km²), the Amazon Basin contains approximately 4.4 million kilometers of channels with a combined surface area of 59,700 km². The intended use of these relationships is for partitioning total floodable area (channels versus lakes and floodplain lakes) obtained from remote sensing for biogeochemical applications (e.g. quantifying CO₂ evasion in the Amazon Basin). Copyright © 2011 John Wiley & Sons, Ltd.     

Surveyed and modelled one‐year morphodynamics in the braided lower Tana River

Field measurements and morphodynamic simulations were carried out along a 5‐km reach of the sandy, braided, lower Tana River in order to detect temporal and spatial variations in river bed modifications and to determine the relative importance of different magnitude discharges on river bed and braid channel evolution during a time span of one year, i.e. 2008–2009. Fulfilling these aims required testing the morphodynamic model's capability to simulate changes in the braided reach. We performed the simulations using a 2‐D morphodynamic model and different transport equations. The survey showed that more deposition than erosion occurred during 2008–2009. Continuous bed‐load transport and bed elevation changes of ±1 m, and a 70–188‐m downstream migration of the thalweg occurred. Simulation results indicated that, during low water periods, modifications occurred in both the main channel and in other braid channels. Thus, unlike some gravel‐bed rivers, the sandy lower Tana River does not behave like a single‐thread channel at low discharge. However, at higher discharge, i.e. exceeding 497 m³/s, the river channel resembled a single‐thread channel when channel banks confined the flow. Although the spring discharge peaks caused more rapid modifications than slower flows, the cumulative volumetric changes of the low water period were greater. The importance of low water period flows for channel modifications is emphasized. Although the 2‐D model requires further improvements, the results were nevertheless promising for the future use of this approach in braided rivers. Copyright © 2013 John Wiley & Sons, Ltd.     

ON THE SHRINKAGE OF RIVER CHANNEL

ＯＮＴＨＥＳＨＲＩＮＫＡＧＥＯＦＲＩＶＥＲＣＨＡＮＮＥＬＣＨＥＮＤｏｎｇ１,ＣＡＯＷｅｎｈｏｎｇ２ａｎｄＺＨＡＮＧＱｉｓｈｕｎ３ＡＢＳＴＲＡＣＴＡｌｏｎｇｗｉｔｈｔｈｅｒａｐｉｄｄｅｖｅｌｏｐｍｅｎｔｏｆｓｏｃｉａｌｅｃｏｎｏｍｙ,ｄｅｖｅｌｏｐｍｅ...     

PUBLICATIONS OF THE IAHS

Abstract

Gauging stations where the stage—discharge relationship is affected by hysteresis due to unsteady flow represent a challenge in hydrometry. In such situations, the standard hydrometric practice of fitting a single-valued rating curve to the available stage—discharge measurements is inappropriate. As a solution to this problem, this study provides a method based on the Jones formula and nonlinear regression, which requires no further data beyond the available stage—discharge measurements, given that either the stages before and after each measurement are known along with the duration of each measurement, or a stage hydrograph is available. The regression model based on the Jones formula rating curve is developed by applying the monoclinal rising wave approximation and the generalized friction law for uniform flow, along with simplifying assumptions about the hydraulic and geometric properties of the river channel in conjunction with the gauging station. Methods for obtaining the nonlinear least-squares rating-curve estimates, while factoring in approximated uncertainty, are discussed. The broad practical applicability and appropriateness of the method are demonstrated by applying the model to: (a) an accurate, comprehensive and detailed database from a hydropower-generated highly dynamic flow in the Chattahoochee River, Georgia, USA; and (b) data from gauging stations in two large rivers in the USA affected by hysteresis. It is also shown that the model is especially suitable for post-modelling hydraulic and statistical validation and assessment.     

基于改善水质的浅水湖泊引调水模式的评价指标

以浅水湖泊引清调水的目的为基础,结合生态水力学要求,兼考虑经济性因子,尝试性地建立了浅水湖泊调水方案评价指标体系.评价指标包括湖泊水质改善效果指标、水力条件指标和经济效益指标.以玄武湖调水模式为例,采用二维水质水量模型对玄武湖不同的引水规模、引水方式、引水口和出水口流量分配的引调水方案进行模拟计算,统计出各个评价指标的结果,综合分析评价出相对最优的引调水方案.该评价方法可推广到其它类似的浅水湖泊引清调水模式的评价中.     

Delineating and quantifying ground water discharge zones using streambed temperatures

Streambed temperature mapping, hydraulic testing using minipiezometers, and geochemical analyses of interstitial water of the streambed were used to delineate the pattern of ground water discharge in a sandy streambed and to develop a flux-based conceptual model for ground water/surface water interactions. A new and simple empirical method was used to relate fluxes obtained from minipiezometer data to streambed temperatures. The relationship allowed flux to be calculated at locations where only streambed temperature measurements were made. Slug testing and potentiomanometer measurements at 34 piezometers indicated ground water discharge ranged from 0.03 to 446 L/m2/day (and possibly as high as 7060 L/m2/day) along a 60 m long by 11 to 14 m wide reach of river. Complex but similar plan-view patterns of flux were calculated for both summer and winter using hundreds of streambed temperatures measured on a 1 by 2 m grid. The reach was dominated by ground water discharge and 5% to 7% of the area accounted for approximately 20% to 24% of the total discharge. < 12% of the total area consisted of recharge zones or no-discharge zones. A conceptual model for ground water/surface water interactions consisting of five different behaviors was developed based on the magnitude and direction of flux across the surface of the streambed. The behaviors include short-circuit discharge (e.g., high-flow springs), high discharge (e.g., preferential flowpaths), low to moderate discharge, no discharge (e.g., horizontal hyporheic or ground water flow), and recharge. Geological variations at depth played a key role in determining which type of flow behavior occurred in the streambed.     

NUMERICAL SIMULATION OF FLOOD DETENTION OPERATION IN MULTIPLY CONNECTED RIVER CHANNEL NETWORK

IINTRODUCTIONAccuratepredictionoftheconsequencesofaflooddetelltion-diversionoperationincomplexriverchannelnetworkisofgreatimportancetothedecisionmakingprocess.InthisregardtWotypesofmodelavailableforsimulatingsuchoperationsarefarfromsatisfactory,duetotheintricatenatUreofprocessesinvolvedintheseoperations.Thefirstmodeltype,i.e.,thewidelyusedhydrologicalmodels,dependsonsystematicwaterstageanddischargemeasurementsatgaugestations,andtheirpredictionsarelimitedtovariationsimwaterstageanddischar…     

Groundwater–surface water interactions in a large semi‐arid floodplain: implications for salinity management

Flow regulation and water diversion for irrigation have considerably impacted the exchange of surface water between the Murray River and its floodplains. However, the way in which river regulation has impacted groundwater–surface water interactions is not completely understood, especially in regards to the salinization and accompanying vegetation dieback currently occurring in many of the floodplains. Groundwater–surface water interactions were studied over a 2 year period in the riparian area of a large floodplain (Hattah–Kulkyne, Victoria) using a combination of piezometric surface monitoring and environmental tracers (Cl⁻, δ²H, and δ¹⁸O). Despite being located in a local and regional groundwater discharge zone, the Murray River is a losing stream under low flow conditions at Hattah–Kulkyne. The discharge zone for local groundwater, regional groundwater and bank recharge is in the floodplain within ∼1 km of the river and is probably driven by high rates of transpiration by the riparian Eucalyptus camaldulensis woodland. Environmental tracers data suggest that the origin of groundwater is principally bank recharge in the riparian zone and a combination of diffuse rainfall recharge and localized floodwater recharge elsewhere in the floodplain. Although the Murray River was losing under low flows, bank discharge occurred during some flood recession periods. The way in which the water table responded to changes in river level was a function of the type of stream bank present, with point bars providing a better connection to the alluvial aquifer than the more common clay‐lined banks. Understanding the spatial variability in the hydraulic connection with the river channel and in vertical recharge following inundations will be critical to design effective salinity remediation strategies for large semi‐arid floodplains. Copyright © 2005 John Wiley & Sons, Ltd.     

Effects of river discharge on hyporheic exchange flows in salmon spawning areas of a large gravel‐bed river

The flow magnitude and timing from hydroelectric dams in the Snake River Basin of the Pacific north‐western US is managed in part for the benefit of salmon. The objective of this research was to evaluate the effects of Hells Canyon Dam discharge operations on hydrologic exchange flows between the river and riverbed in Snake River fall Chinook salmon spawning areas. Interactions between river water and pore water within the upper 1 m of the riverbed were quantified through the use of self‐contained temperature and water level data loggers suspended inside of piezometers. The data were recorded at 20 min intervals over a period of 200 days when the mean daily discharge was 218–605 m³ s^?1, with hourly stage changes as large as 1·9 m. Differences in head pressure between the river and riverbed were small, often within ± 2 cm. Measured temperature gradients in the riverbed indicated significant interactions between the surface and subsurface water. At the majority of sites, neither hydraulic nor temperature gradients were significantly affected by either short‐ or long‐term changes in discharge operations from Hells Canyon Dam. Only 2 of 14 study sites exhibited acute flux reversals between the river and riverbed resulting from short‐term, large magnitude changes in discharge. The findings suggest that local scale measurements may not be wholly explanatory of the hydrological exchange between the river and riverbed. The processes controlling surface water exchange at the study sites are likely to be bedform‐induced advective pumping, turbulence at the riverbed surface, and large‐scale hydraulic gradients along the longitudinal profile of the riverbed. By incorporating the knowledge of hydrological exchange processes into water management planning, regional agencies will be better prepared to manage the limited water resources among competing priorities that include salmon recovery, flood control, irrigation supply, hydropower production, and recreation. Copyright © 2007 John Wiley & Sons, Ltd.     

洞庭湖城陵矶水道水力几何形态的研究

根据１９５１－１９８８年洞庭湖城陵矶站的水文测验资料,运用Ｌ．Ｂ．Ｌｅｏｐｏｌｄ河床力几何形态原理,建立洞庭湖出口－城陵矶水道河相关系式,研究该水道水力几何形态的特点及变化。研究表明,与河流水道相比,洞庭湖出口水道河宽指数ｂ随流量的变化较小,而水深指数ｆ及流速指数ｍ随流量的变化较大,河床横面具有窄深的特点。     

Assessment of river flow with significant lateral inflow through reverse routing modeling

The discharge hydrograph estimation in rivers based on reverse routing modeling and using only water level data at two gauged sections is here extended to the most general case of significant lateral flow contribution, without needing to deploy rainfall–runoff procedures. The proposed methodology solves the Saint‐Venant equations in diffusive form also involving the lateral contribution using a “head‐driven” modeling approach where lateral inflow is assumed to be function of the water level at the tributary junction. The procedure allows to assess the discharge hydrograph at ends of a selected river reach with significant lateral inflow, starting from the stage recorded there and without needing rainfall data. Specifically, the MAST 1D hydraulic model is applied to solve the diffusive wave equation using the observed stage hydrograph at the upstream section as upstream boundary condition. The other required data are (a) the observed stage hydrograph at the downstream section, as benchmark for the parameter calibration, and (b) the bathymetry of the river reach, from the upstream section to a short distance after the downstream gauged section. The method is validated with different flood events observed in two river reaches with a significant intermediate basin, where reliable rating curves were available, selected along the Tiber River, in central Italy, and the Alzette River, in Luxembourg. Very good performance indices are found for the computed discharge hydrographs at both the channel ends and along the tributaries. The mean Nash‐Sutcliffe value (NS_q) at the channel ends of two rivers is found equal to 0.99 and 0.86 for the upstream and downstream sites, respectively. The procedure is also validated on a longer stretch of the Tiber River including three tributaries for which appreciable results are obtained in terms of NS_q for the computed discharge hydrographs at both the channel ends for three investigated flood events.     

Late Holocene flood magnitudes in the Lower Rhine river valley and upper delta resolved by a two-dimensional hydraulic modelling approach

Palaeoflood hydraulic modelling is essential for quantifying ‘millennial flood’ events not covered in the instrumental record. Palaeoflood modelling research has largely focused on one-dimensional analysis for geomorphologically stable fluvial settings because two-dimensional analysis for dynamic alluvial settings is time consuming and requires a detailed representation of the past landscape. In this study, we make the step to spatially continuous palaeoflood modelling for a large and dynamic lowland area. We applied advanced hydraulic model simulations (1D–2D coupled set-up in HEC-RAS with 950 channel sections and 108 × 10³ floodplain grid cells) to quantify the extent and magnitude of past floods in the Lower Rhine river valley and upper delta. As input, we used a high-resolution terrain reconstruction (palaeo-DEM) of the area in early mediaeval times, complemented with hydraulic roughness values. After conducting a series of model runs with increasing discharge magnitudes at the upstream boundary, we compared the simulated flood water levels with an inventory of exceeded and non-exceeded elevations extracted from various geological, archaeological and historical sources. This comparison demonstrated a Lower Rhine millennial flood magnitude of approximately 14,000 m³/s for the Late Holocene period before late mediaeval times. This value exceeds the largest measured discharges in the instrumental record, but not the design discharges currently accounted for in flood risk management.