A combined hydrologic and hydraulic modeling approach for testing efficiency of structural flood control measures

The necessity of estimating the degree and spatial extent of positive impacts with regard to protecting communities and properties through potential flood control projects can be considered one of the main reasons for performing flood modeling. This paper presents an overall systematic approach based on the simulation of some extreme event conditions, using a hydrological model to generate the resulting river flows and then using a hydraulic modeling exercise to decide upon floodplain evolution in the case-study area, Bostanli river basin, which has been under the threat of flooding for many years. The potential serviceability of the planned Bostanli Dam in the study area was examined by using the HEC-HMS and HEC-RAS modeling tools, both integrated with GIS functions for spatial operations. The results indicate that the dam construction as planned would have a somewhat positive impact as a potential flood control measure, since it seems to decrease the flood peaks of 68.9 and 158.7 m³/s (that would potentially be generated by 100- and 500-year storm events under current conditions) to 65.5 and 150.7 m³/s (when the dam is in operation), respectively. However, this seems to contribute little to the overall flood mitigation performance in the basin.     

Stability of landslide dams and development of knickpoints

The Wenchuan earthquake triggered many landslides and numerous avalanches and created 100 odd quake lakes. The quake lakes may be removed or preserved. The removal strategy was applied to several large landslide dams, which were dangerous because massive amounts of water pooled up in the quake lakes. The dams could eventually fail under the action of dam outburst flooding, potentially endangering the lives of people in the downstream reaches. This paper studied the stability of landslide dams and the development of knickpoints by field investigations and experiments, and analyzing satellite images. The study concluded that if landslide dams were preserved, they would develop into knickpoints and act as a primary control of riverbed incision and, thus, reduce the potential of new landslide. The stability of landslide dams depends mainly on the development of the step-pool system and stream power of the flood flow. If a landslide dam consists of many boulders, a step-pool system may develop on the spillway channel of the dam, which would maximize the resistance, consume most of the flow energy and consequently protect the dam from incision. The development degree of the step-pool system is represented by a parameter S _p, which was measured with a specially designed instrument. A preservation ratio of landslide dams is defined as the ratio of preserved height after flood scouring to the original height of the dam. For streams with peak flood discharge lower than 30 m³/s, the preservation ratio is linearly proportional to S _p. For rivers with a peak flood discharge higher than 30 m³/s (30–30,000 m³/s), the minimum S _p value for stable channel increases with log p, in which p is the unit stream power. For a landslide dam with a poorly developed step-pool system, S _p is smaller than the minimum value and the outburst flood incises the spillway channel and causes failure of the dam. For preserved landslide dams, sediment deposits in the quake lakes. A landslide dam may develop into a knickpoint if it is stabilized by long-term action of the flow. Large knickpoints can totally change the fluvial processes and river morphology. Uplift of the Qinghai–Tibetan Plateau has caused extensive channel bed incision along almost all rivers. For many rivers, the incision has been partly controlled by knickpoints. Upstream reaches of a knickpoint have a new and unchanging base level. This brings about a transition from degradation to aggradation and from vertical bed evolution to horizontal fluvial process. Multiple and unstable channels are prominent in the reaches, upstream of the knickpoints. If hundreds of landslide dams occurred simultaneously on a reach of a mountain river, the potential energy of bank failure and the slope erosion would be greatly reduced and sediment yield from the watershed may be reduced to nearly zero. The quake lakes may be preserved long term and become beautiful landscapes. Streams with long-term unfilled quake lakes have good aquatic ecology.     

Sulfur isotope studies to quantify sulfate components in water of flooded lignite open pits – Lake Goitsche,Germany

A system of connected lignite mining pits (part of the former Goitsche mining complex, Germany) was flooded with river water between 1999 and 2002. A considerable accumulation of acid associated with oxidized sulfides in sediments was seen as a critical point for the development of the lake water. To characterize the components contributing to the supply of dissolved lake water SO₄ hydro-chemical and isotope investigations with respect to groundwater, pore water in the sulfide bearing sediments, river water and lake water were performed. δ³⁴S of pore water SO₄ that was dominated by oxidized pyrites ranges around −25‰ VCDT and differs strongly from river water SO₄ with about +4.4‰. Thus, interactions between lake water and sediments were particularly pronounced during the first phase of flooding. For this period, a more quantitative estimation of the SO₄ components in the lake water was difficult because of the heterogeneous SO₄ distributions between the different sub-basins of the lake and according to the flooding process itself. Later, a component separation was attempted following mixing of the whole lake, which first occurred in spring 2002. A very heterogeneous groundwater environment with respect to highly variable SO₄ concentrations and δ³⁴S values and changing interaction with the forming lakes proved to be one of the most important limitations in the calculations of the mixing.     

Combined numerical investigation of the Gangda paleolandslide runout and associated dam breach flood propagation in the upper Jinsha River,SE Tibetan Plateau

A large paleolandslide occurred opposite the Gangda village in the upper Jinsha River, SE Tibetan Plateau. Field geological investigations and remote sensing indicated that the Gangda paleolandslide once blocked the Jinsha River. Evidence of river blocking, including landslide dam relics, upstream lacustrine sediments, and downstream outburst sediments, has been well preserved. To understand the river-blocking event including landslide, dam breach, and associated outburst flooding, optically stimulated luminescence (OSL) dating and numerical simulations were performed in this study. OSL dating results showed that the paleolandslide dam was formed at 5.4?±?0.5 ka BP and breached at 3.4?±?0.3 ka BP, indicating that the dam lasted approximately 2000 years. The discrete element method was used to simulate the dynamics of the Gangda rock landslide based on the restored topography, while a fluid–solid coupling model was performed to simulate the landslide dam breaching and flooding. The fluid–solid coupling model can simultaneously reflect the process of landslide-dam collapse and the propagation of outburst flood. The simulated results indicate that the whole landslide process lasted about 60 s with a peak velocity of 38 m/s. It is significant that the simulated morphology of the residual landslide dam and downstream outburst sediments is consistent with the field observations. The combined numerical investigation in this paper provided new insights into the research of landscape evolution and helped to understand the chain disaster of landslide, dam breach, and flooding.

     

Indirect Effects of Salinity and Temperature on Kuwait’s Shrimp Stocks

Discharge of the Shatt Al-Arab is believed to be a dominating component of the northern Arabian Gulf’s ecology and largely responsible for productivity of Kuwait’s fisheries. With major construction of dams on the Tigris and Euphrates Rivers in Turkey, river discharge has been substantially reduced, and flooding essentially eliminated. We attempted to relate river flow and shrimp landings indirectly by correlating 19 years of salinity and temperature data from Kuwait’s waters with corresponding annual shrimp landings. For green tiger prawns (Penaeus semisulcatus), the combination of October salinities and January temperatures provided the best correlation (r = 0.67) with landings in the following shrimping season. For the combined landings of jinga and kiddi shrimps (Metapenaeus affinis and Parapenaeopsis stylifera, respectively), December salinity during season and May temperature prior to season resulted in the best correlation (r = 0.87). Landings of these two species also correlated well with spring and summer temperatures. Under normal conditions, late winter or early spring temperatures prior to fishing season influence recruitment of the green tiger prawn, whereas December salinities during harvest season influence the abundance of jinga–kiddi shrimps. With further reductions in the Shatt Al-Arab discharge and the elimination of flood events, Kuwait’s shrimp landings will most likely decrease over time.     

Fashion and phases of late Pleistocene aggradation and incision in the Alaknanda River Valley, western Himalaya, India

We study the aggradation and incision of the Alaknanda River Valley during the late Pleistocene and Holocene. The morphostratigraphy in the river valley at Deoprayag shows the active riverbed, a cut terrace, and a fill terrace. The sedimentary fabric of the fill terrace comprises four lithofacies representing 1) riverbed accretion, 2) locally derived debris fan, 3) the deposits of waning floods and 4) palaeoflood records. The sedimentation style, coupled with geochemical analysis and Optically Stimulated Luminescence (OSL) dating, indicate that this terrace formed in a drier climate and the river valley aggraded in two phases during 21–18 ka and 13–9 ka. During these periods, sediment supply was relatively higher. Incision began after 10 ka in response to a strengthened monsoon and aided by increase of the tectonic gradient. The cut terrace formed at ~ 5 ka during a phase of stable climate and tectonic quiescence. The palaeoflood records suggest wetter climate 200–300 yr ago when the floods originated in the upper catchment of the Higher Himalaya and in the relatively drier climate ~ 1.2 ka when locally derived sediments from the Lesser Himalaya dominated flood deposits. Maximum and minimum limits of bedrock incision rate at Deoprayag are 2.3 mm/a and 1.4 mm/a.     

Peatland development at the arctic tree line (Québec, Canada) influenced by flooding and permafrost

In this study, we documented the Holocene history of a peat plateau at the arctic tree line in northern Québec using stratigraphic and macrofossil analyses to highlight the effects of geomorphic setting in peatland development. Paludification of the site began about 6800 cal yr BP. From 6390 to 4120 cal yr BP, the peatland experienced a series of flooding events. The location of the peatland in a depression bounded by two small lakes likely explains its sensitivity to runoff. The proximity of a large hill bordering the peatland to the south possibly favored the inflow of mineral-laden water. The onset of permafrost aggradation in several parts of the peatland occurred after 3670 cal yr BP. Uplifting of the peatland surface caused by permafrost stopped the flooding. According to radiocarbon dating of the uppermost peat layers, permafrost distribution progressed from the east to the west of the peatland, indicating differential timing for the initiation of permafrost throughout the peatland. Most of the peatland was affected by permafrost growth during the Little Ice Age. Picea mariana macroremains at 6450 cal yr BP indicate that the species was present during the early stages of peatland development, which occurred soon after the sea regression.     

Formation,failure, and consequences of the Xiaolin landslide dam,triggered by extreme rainfall from Typhoon Morakot,Taiwan

An extreme rainfall event on August 9, 2009, which was close to setting a world record for 48-h accumulated rainfall, induced the Xiaolin deep-seated landslide, which was located in southwestern Taiwan and had volume of 27.6?×?10⁶?m³, and caused the formation of a landslide dam. The landslide dam burst in a very short time, and little information remained afterward. We reconstructed the process of formation and failure of the Xiaolin landslide dam and also inferred the area of the impoundment and topographic changes. A 5?×?5-m digital elevation model, the recorded water stage of the Qishan River, and data from field investigation were used for analysis. The spectral magnitude of the seismic signals induced by the Xiaolin landslide and flooding due to failure of the landslide dam were analyzed to estimate the timing of the dam breach and the peak discharge of the subsequent flood. The Xiaolin landslide dam failure resulted from overtopping. We verified the longevity of the Xiaolin landslide dam at about 2 h relying on seismic signals and water level records. In addition, the inundated area, volume of the impoundment behind the Xiaolin landslide dam, and peak discharge of the flood were estimated at 92.3 ha, 19.5?×?10⁶?m³, and 17?×?10³?m³/s, respectively. The mean velocity of the flood-recession wave front due to the dam blockage was estimated at 28 km/h, and the peak flooding velocity after failure of the dam was estimated at 23 km/h. The Xiaolin landslide provides an invaluable opportunity for understanding the mechanism of deep-seated landslides and flooding processes following a landslide dam failure.     

Stochastic Modelling of the Impact of Flood Protection Measures Along the River Waal in the Netherlands

River flooding is a problem of international interest. In the past few years many countries suffered from severe floods. A large part of the Netherlands is below sea level and river levels. The Dutch flood defences along the river Rhine are designed for water levels with a probability of exceedance of 1/1250 per year. These water levels are computed with a hydrodynamic model using a deterministic bed level and a deterministic design discharge. Traditionally, the safety against flooding in the Netherlands is obtained by building and reinforcing dikes. Recently, a new policy was proposed to cope with increasing design discharges in the Rhine and Meuse rivers. This policy is known as the Room for the River (RfR) policy, in which a reduction of flood levels is achieved by measures creating space for the river, such as dike replacement, side channels and floodplain lowering. As compared with dike reinforcement, these measures may have a stronger impact on flow and sediment transport fields, probably leading to stronger morphological effects. As a result of the latter the flood conveyance capacity may decrease over time. An a priori judgement of safety against flooding on the basis of an increased conveyance capacity of the river can be quite misleading. Therefore, the determination of design water levels using a fixed-bed hydrodynamic model may not be justified and the use of a mobile-bed approach may be more appropriate. This problem is addressed in this paper, using a case study of the river Waal (one of the Rhine branches in the Netherlands). The morphological response of the river Waal to a flood protection measure (floodplain lowering in combination with summer levee removal) is analysed. The effect of this measure is subject to various sources of uncertainty. Monte Carlo simulations are applied to calculate the impact of uncertainties in the river discharge on the bed levels. The impact of the “uncertain” morphological response on design flood level predictions is analysed for three phenomena, viz. the impact of the spatial morphological variation over years, the impact of the seasonal morphological variation and the impact of the morphological variability around bifurcation points. The impact of seasonal morphological variations turns out to be negligible, but the other two phenomena appear to have each an appreciable impact (order of magnitude 0.05–0.1 m) on the computed design water levels. We have to note however, that other sources of uncertainty (e.g. uncertainty in hydraulic roughness predictor), which may be of influence, are not taken into consideration. In fact, the present investigation is limited to the sensitivity of the design water levels to uncertainties in the predicted bed level.     

Environmental Forcing of Phytoplankton in a Mediterranean Estuary (Guadiana Estuary,South-western Iberia): A Decadal Study of Anthropogenic and Climatic Influences

Phytoplankton seasonal and interannual variability in the Guadiana upper estuary was analyzed during 1996–2005, a period that encompassed a climatic controlled reduction in river flow that was superimposed on the construction of a dam. Phytoplankton seasonal patterns revealed an alternation between a persistent light limitation and episodic nutrient limitation. Phytoplankton succession, with early spring diatom blooms and summer–early fall cyanobacterial blooms, was apparently driven by changes in nutrients, water temperature, and turbulence, clearly demonstrating the role of river flow and climate variability. Light intensity in the mixed layer was a prevalent driver of phytoplankton interannual variability, and the increased turbidity caused by the Alqueva dam construction was linked to pronounced decreases in chlorophyll a concentration, particularly at the start and end of the phytoplankton growing period. Decreases in annual maximum and average abundances of diatoms, green algae, and cyanobacteria were also detected. Furthermore, chlorophyll a decreases after dam filling and a decrease in turbidity may point to a shift from light limitation towards a more nutrient-limited mode in the near future.     

云贵高原斜坡地带典型地下水富硫酸盐地区“越层找水”模式及其机理研究

地下水是中国西南云贵高原斜坡地带重要的饮用水源,其中黔中镇宁县北部存在大面积地下水中硫酸盐超标地区,导致了当地出现水质性缺水问题。本次研究通过岩相古地理分析、水化学分析、D、18O、34S同位素测试、岩样测试、岩心观察等技术手段,查明了研究区内地下水的补给源主要为大气降水,地下水中硫酸盐(SO_4~(2-))浓度在30～1100 mg/L,平均值为221.78 mg/L,呈现高度富集SO_4~(2-)的特点;在区域相对隔水层以上,不同类型地下水中SO_4~(2-)浓度大体具有机井上升泉下降泉溶潭地下河出口基岩裂隙泉表层岩溶泉的规律;探讨了地表水、地下水中SO_4~(2-)的来源,表层岩溶泉中硫酸根的来源主要为大气降水,高硫酸盐样品中硫酸根的来源主要为石膏溶解。在此基础上,结合钻井资料,掌握了研究区膏岩层分布及含水层结构特征,通过建立"越层找水"模式,采取下层低硫酸盐含水层,获取合格饮用水源,可有效解决当地水质性缺水问题。     

A Coupled Riverine-Marine Fractionation Model for Dissolved Rare Earths and Yttrium

Fractionation of yttrium (Y) and the rare earth elements (REEs) begins in riverine systems and continues in estuaries and the ocean. Models of yttrium and rare earth (YREE) distributions in seawater must therefore consider the fractionation of these elements in both marine and riverine systems. In this work we develop a coupled riverine/marine fractionation model for dissolved rare earths and yttrium, and apply this model to calculations of marine YREE fractionation for a simple two-box (riverine/marine) geochemical system. Shale-normalized YREE concentrations in seawater can be expressed in terms of fractionation factors ( _ij ) appropriate to riverine environments ( ) and seawater ( ):

The reconstruction of a glacial lake outburst flood using HEC-RAS and its significance for future hazard assessments: an example from Lake 513 in the Cordillera Blanca,Peru

In April 2010, an ice/rockfall into Lake 513 triggered a glacial lake outburst flood (GLOF) along the Chucchun River in the Cordillera Blanca of Peru. This paper reconstructs the hydrological characteristics of this as yet undocumented event using a 1D flood model prepared with HEC-RAS. The principle model inputs were obtained during detailed field surveys of surface characteristics and topography within the river and across the adjacent floodplain; a total of 120 cross-sections were surveyed. These inputs were refined further by eyewitness accounts and additional geomorphological observations. The flood modelling has enabled us to constrain the extent of the water surface and its elevation at each cross-section in addition to defining the peak discharge (580 m³ s^?1). These modelling results show good agreement with other information about the flood including: flood marks and minimum flood levels; the lake displacement wave height; the extent of the flooded area; and the travel time from Lake 513 to the confluence with the Santa River. This demonstrates that the model offers a reliable reconstruction of the basic hydrological characteristics of the GLOF. It provides important information about the flood intensity and significantly improves our ability to model future flood scenarios along both the studied river and within neighbouring catchments. The flood hazard, defined by the flood depth during peak discharge, shows that the majority of the damaged infrastructure (houses, bridges, and a drinking water treatment plant) was only subjected to low or medium flood intensities (defined by a maximum water depth of less than 2 m). These low flood intensities help to explain why the flooding caused comparatively minor damage despite the significant public attention it attracted.     

基于GIS的大连市金普新区洪水淹没分析

采用25 m分辨率的DEM并收集多年洪水位、潮位数据,以保证源数据的精度和可靠性.结合研究区地貌特征和水系分布,基于ArcGIS的空间分析工具采用种子蔓延法进行了有源淹没分析计算,得到淹没区范围和水深分布图.结果表明：淹没区主要分布在水库和河流下游的冲洪积平原和海积平原,淹没面积共318.08 km2,占全区陆域面积的18.6%.洪水淹没深度为0~7.6 m.以水库为种子点的淹没区域较大,淹没水深自种子点到海岸线逐渐增大;以河流为种子点的淹没区域呈带状分布,淹没水深自河流中线向两侧逐渐减小.     

新疆天山南麓柯坪河水文特性与洪水分析

柯坪河是天山南麓典型的雨水和泉水补给河流,河水主要为地下泉水出露补给,流量比较稳定,但暴雨的出现经常引起洪水发生.柯坪河仅有不足一年的水文观测,但河水主要是泉水补给,水量非常稳定,根据野外调查和推算,分析了柯坪河年内水文变化特征.依据洪水调查和文献查证,结合气象观测资料,分析了近百年来的典型洪水事件,建立了不同几率洪水发生的洪峰流量和洪水过程.最近几十年来的气温升高和降水增加,也使暴雨洪水的强度加强,并且频次增加.应加强对极端气候事件引发的洪水的监测和应对,在农业和水库安全运营上应注意和加强建立应对气候变化的措施.     

Quantitative dam break analysis on a reservoir earth dam

Mathematical simulations on dam break or failure using BOSS DAMBRK hydrodynamic flood routing dam break model were carried out to determine the extent of flooding downstream, flood travel times, flood water velocities and impacts on downstream affected residences, properties and environmental sensitive areas due to floodwaters released by failure of the dam structure. Computer simulations for one of the worse case scenarios on dam failure using BOSS DAMBRK software accounted for dam failure, storage effects, floodplains, over bank flow and flood wave attenuation. The simulated results reviewed a maximum flow velocity of 2.40 m/s with a discharge of approximately 242 mз /s occurred at 1.00 km downstream. The maximum discharge increased from 244 m³/s (flow velocity = 1.74 m/s occurred at 8^th. km) to 263 m³/s (flow velocity = 1.37 m/s occurred at 12^th. km); about a 39% drop in flow velocity over a distance of 4.00 km downstream. If the entire dam gives way instantly, some spots stretching from 0.00 km (at dam site) to approximately 3.40 km downstream of the dam may be categorized as “danger zone”, while downstream hazard and economic loss beyond 3.40 km downstream can be classified as “low” or “minimal” zones.     

The Basin of Mexico and its metropolitan area: water abstraction and related environmental problems

The Basin of Mexico is a closed basin of lacustrine character, with an average elevation of 2200 m above sea level. The watershed covers a vast extension in five states. Mexico City and its metropolitan area are located within this basin. The aquifer system is the main source of water supply for more than 20 million people. Water consumption is about 60 m³/s. The aquifer supplies about 43 m³/s from around 1000 wells at 70–200 m depth. Pumping policies have generated subsidence and degradation of the ground water quality in the Basin of Mexico The lacustrian clay layers play an important role in the local hydrogeology, protecting the aquifer from pollution, but the transition and piedmont areas are highly vulnerable to surface pollutants.     

Pairing geotechnics and fluvial hydraulics for the prediction of the hazard zones of an exceptional flooding

The direct consequences of exceptional floods are usually considered to be limited to the maximum flooding zone created downstream. However, considering the magnitude of the flows, the morphology of the flooded zone could undergo deep changes. To predict the hazard zone on a river undergoing exceptional flooding, numerical simulations are widely used. In this article, the simulation of the evolution of river reaches resulting from such catastrophic events is performed by coupling the hydraulic and sediment transport numerical model GSTARS with a developed slope stability model based on the Bishop’s simplified method. This is a novel methodology for the delimitation of hazard zones along riverbanks by taking into consideration not only the flood risks but also the possible induced landslides. Indeed, each section of the river reach is subject to changes caused by the river hydraulics and the associated erosion or sediment deposition and also undergoes profile changes caused by possible landslides. The initial hydraulic and geotechnical characteristics are first defined and then used to test the stability of several slopes of representative sections of the river reaches before the dam break. Validation tests are performed on specific reaches of the Outaouais River (Quebec) undergoing a dam break flood.     

Landslide dam failure and flood hydraulics. Part I: experimental investigation

Landslide dam failure can trigger catastrophic flooding in the downstream. However, field observation of such flooding is rarely available, while laboratory experimental studies are sparse. The mechanism of landslide dam failure and the flood has so far remained insufficiently understood. Here, we present an experimental investigation of landslide dam failure and the flood. A total of 28 runs of experiments are carried out in a flume of 80 m × 1.2 m × 0.8 m, with differing inflow discharge, dam composition, dam geometry, and initial breach dimension. An array of twelve automatic water-level probes is deployed to measure the stage hydrographs along the flume, and the video recording of the dam failure processes facilitates an estimation of the widening of initial breach. Under the present experimental conditions with dams composed of homogeneous materials, landslide dam failure is primarily caused by erosion of overtopping flow, and lateral mass collapse is also considerable during the cause of breach widening. Cohesive clay may act to mitigate the seepage through the dam and thus its subsidence and appreciably modulate the dam failure process and the flood. However, the impacts of clay may be readily overwhelmed by a large inflow discharge and initial breach. Gravels in the dam may appreciably depress the rate of the dam failure process and thus modify the flood. The present work provides new experimental data set for testing mathematical models of the flood flow due to landslide dam failure.     

易贡滑坡堰塞湖溃坝洪水分析

滑坡堰塞坝体主要由块石、碎石土等松散材料组成,随着上游水位的不断上升,极易失稳,一旦决口将对给下游人民的生命财产安全造成极大的威胁。因此,研究堰塞坝溃坝问题具有重要的学术意义和应用价值。2000 年 4 月 9 日,西藏林芝地区波密县易贡藏布河扎木弄沟发生大规模山体滑坡堵塞易贡藏布江,形成坝高60m,长约2500m,库容可达288108m3,体积约28108～30108m3的滑坡堰塞湖, 2000年6月10日堰塞坝溃决。本文以易贡堰塞湖溃坝为例,从连续性方程及Navier Stokes方程出发,结合标准型湍流模型,并采用VOF方法进行自由面处理,基于流体计算软件Fluent模拟分析了溃坝洪水在下游弯曲河道的演进过程及不同位置的流速变化。数值模拟结果与实测资料记录基本一致,表明该模型能够模拟溃坝洪水在地形复杂弯曲河道中的演进过程。