An integrated water balance model for assessing water scarcity in a data-sparse interfluve in eastern India

Abstract

The objective of this study is to measure the balance of water demand versus water resource availability in an interfluve of West Bengal, India to support water resource planning, particularly of inter-basin transfers. Surface water availability was modelled using the US Soil Conservation Service curve number (SCS-CN) approach, whilst groundwater availability was modelled based on water-level fluctuations and the rainfall infiltration method. Water use was modelled separately for the agricultural, industrial, and domestic sectors using a predominantly normative approach and water use to availability ratios calculated for different administrative areas within the interfluve. Overall, the approach suggested that the interfluve receives 327 × 10⁶ m³ year^-1 of excess water after satisfying these sectoral demands, but that the eastern part of the study area is in deficit. However, a sensitivity analysis carried on the approach to several assumptions in the model suggested changed circumstances would produce surplus/deficit ranging from ?215 × 10⁶ to 435 × 10⁶ m³ year^-1 . The approach could have potential for localised water balance modelling in other Indian catchments.

Editor D. Koutsoyiannis; Associate editor D. Hughes     

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.     

Geomorphological effects of catastrophic flooding on a small subtropical Estuary

The Mgeni Estuary is situated on the east coast of South Africa. Tidal salinity changes commonly extend 2.5 km upstream from the mouth. The subtropical climatic regime causes seasonal flooding, the geomorphological effects of which are increased by a steep river gradient. A 120-150 year flood event during September 1987, with an estimated peak discharge of 10800 cumecs, resulted in the erosion of 1.86 × 10⁶ m³ of sediment from the estuary. The normal high tide volume of the estuary (0.35 × 10⁶ m³) increased to 1.85 × 10⁶ m³ after the flood. The post-flood intertidal volume (tidal prism) of the estuary was 0.43 × 10⁶ m³ compared to 0.19 × 10⁶ m³ before the event. The nature and extent of erosion in the estuary was controlled by the composition of the estuary banks and bed material. Gravel and sand substrates were preferentially eroded from the estuary bed before the cohesive bank materials were undercut in sections of the estuary. Mud and fine sand were deposited on overbank areas when the flood waned and a lag of coarse sand and gravel was produced on the estuary bed. Approximately 46 per cent of the bedload sediment supplied to the estuary since 1917 was retained in the estuary until the 1987 flood. The remainder was transported through the system and into the Indian Ocean. Results show that catastrophic floods play an important role in sedimentation in small estuaries.     

Soil erosion and sediment interception by check dams in a watershed for an extreme rainstorm on the Loess Plateau,China

The magnitude of soil erosion and sediment load reduction efficiency of check dams under extreme rainstorms is a long-standing concern. The current paper aims to use check dams to deduce the amount of soil erosion under extreme rainstorms in a watershed and to identify the difference in sediment interception efficiency of different types of check dams. Based on the sediment deposition at 12 check dams with 100% sediment interception efficiency and sub-catchment clustering by taking 12 dam-controlled catchments as clustering criteria, the amount of soil erosion resulting from an extreme rainstorm event on July 26, 2017 (named “7·26” extreme rainstorm) was estimated in the Chabagou watershed in the hill and gully region of the Loess Plateau. The differences in the sediment interception efficiency among the check dams in the watershed were analyzed according to field observations at 17 check dams. The results show that the average erosion intensity under the “7–26” extreme rainstorm was approximately 2.03 × 10⁴ t/km², which was 5 times that in the second largest erosive rainfall in 2017 (4.15 × 10³ t/km²) and 11–384 times that for storms in 2018 (0.53 × 10² t/km² - 1.81 × 10³ t/km²). Under the “7–26” extreme rainstorm, the amount of soil erosion in the Chabagou watershed above the Caoping hydrological station was 4.20 × 10⁶ t. The sediment interception efficiency of the check dams with drainage canals (including the destroyed check dams) and with drainage culverts was 6.48 and 39.49%, respectively. The total actual sediment amount trapped by the check dams was 1.11 × 10⁶ t, accounting for 26.36% of the total amount of soil erosion. In contrast, 3.09 × 10⁶ t of sediment were input to the downstream channel, and the sediment deposition in the channel was 2.23 × 10⁶ t, accounting for 53.15% of the total amount of soil erosion. The amount of sediment transport at the hydrological station was 8.60 × 10⁵ t. The Sediment Delivery Ratio (SDR) under the “7·26” extreme rainstorm was 0.21. The results indicated that the amount of soil erosion was huge, and the sediment interception efficiency of the check dams was greatly reduced under extreme rainstorms. It is necessary to strengthen the management and construction technology standards of check dams to improve the sediment interception efficiency and flood safety in the watershed.     

Effect of bed topography on soil aggregates transport by rill flow

After its formation, a rill may remain in the field for months, often receiving lower flow rates than the formative discharge. The objective of this work was to evaluate the rill flow transport capacity of soil aggregates at discharges unable to erode the rill, and to analyse the influence of the rill macro‐roughness on this transport process. A non‐erodible rill was built in which roughness was reproduced in detail. In order to assess only the rill macro‐roughness, a flat channel with a similar micro‐roughness to that in the rill replica was built. Rill and channel experiments were carried out at a slope of 8 and at six discharges (8·3 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the rill, and eight discharges (1·6 × 10^?5 to 5·2 × 10^?4 m³ s^?1) in the channel. Non‐erodible aggregates of three sizes (1–2, 3–5 and 5–10 mm) were released at the inlet of the rill/channel. The number of aggregates received at the outlet was registered. The number and position of the remaining aggregates along the rill/channel were also determined. The rill flow was a major sediment transport mechanism only during the formation of the rill, as during that period the power of the flow was great enough to overcome the influence of the macro‐roughness of the rill bed. At lower discharges the transport capacity in the previously formed rill was significantly less than that in the flat channel under similar slope and discharge. This was determined to be due to local slowing of flow velocities at the exit of rill pools. Copyright © 2006 John Wiley & Sons, Ltd.     

A remote sensing contribution to hydrologic modelling in arid and inaccessible watersheds,Pishin Lora basin,Pakistan

The lack of adequate field measurements often hampers the construction and calibration of rainfall‐runoff models over many of the world's watersheds. We adopted methodologies that rely heavily on readily available remote sensing datasets as viable alternatives for assessing, managing, and modelling of such remote and inadequately gauged regions. The Soil and Water Assessment Tool was selected for continuous (1998–2005) rainfall‐runoff modelling of one such area, the northeast part of the Pishin Lora basin (NEPL). Input to the model included satellite‐based Tropical Rainfall Measuring Mission precipitation data, and modelled runoff was calibrated against satellite‐based observations, the latter included: (i) monthly estimates of the water volumes impounded by the Khushdil Khan (latitude 30°40′N, longitude 67°40′E), and the Kara Lora (latitude 30°34′N, longitude 66°52′E) reservoirs, and (ii) inferred wet versus dry conditions in streams across the NEPL. Calibrations were also conducted against observed flow reported from the Burj Aziz Khan station at the NEPL outlet (latitude 30°20′N; longitude 66°35′E). Model simulations indicate that (i) average annual precipitation (1998–2005), runoff and recharge in the NEPL are 1300 × 10⁶ m³, 148 × 10⁶ m³, and 361 × 10⁶ m³, respectively; (ii) within the NEPL watershed, precipitation and runoff are high for the northeast (precipitation: 194 mm/year; runoff: 38 × 10⁶ m³/year) and northwest (134 mm/year; 26 × 10⁶ m³/year) basins compared to the southern basin (124 mm/year; 8 × 10⁶ m³/year); and (3) construction of delay action dams in the northeast and northwest basins could increase recharge from 361 × 10⁶ m³/year up to 432 × 10⁶ m³/year and achieve sustainable extraction. The adopted methodologies are not a substitute for traditional approaches, but they could provide first‐order estimates for rainfall, runoff, and recharge in the arid and semi‐arid parts of the world that are inaccessible and/or lack adequate coverage with field data. Copyright © 2011 John Wiley & Sons, Ltd.     

The adjustment of the meandering pattern of the lower Jordan River to change in water discharge

The annual discharge of the Lower Jordan River has decreased from about 1250 × 10⁶ m³ y^?1 to about 300 × 10⁶ m³ y^?1 due to water exploitation. The decrease in water discharge was not followed by a similar decrease in the suspended sediment. Evidence from air photographs and maps from 1920s to early 1970s show that the river is adjusting itself by shortening its course. The shortening is not distributed uniformly along the valley length. Local effects obscure the effect of decreasing discharge. Since the early 1970s the channel is becoming longer again due to a change in the dominant flow.     

Debris-flow-dominated sediment transport through a channel network after wildfire

Field studies that investigate sediment transport between debris-flow-producing headwaters and rivers are uncommon, particularly in forested settings, where debris flows are infrequent and opportunities for collecting data are limited. This study quantifies the volume and composition of sediment deposited in the arterial channel network of a 14-km² catchment (Washington Creek) that connects small, burned and debris-flow-producing headwaters (<1 km²) with the Ovens River in SE Australia. We construct a sediment budget by combining new data on deposition with a sediment delivery model for post-fire debris flows. Data on deposits were plotted alongside the slope–area curve to examine links between processes, catchment morphometry and geomorphic process domains. The results show that large deposits are concentrated in the proximity of three major channel junctions, which correspond to breaks in channel slope. Hyperconcentrated flows are more prominent towards the catchment outlet, where the slope–area curve indicates a transition from debris flow to fluvial domains. This shift corresponds to a change in efficiency of the flow, determined from the ratio of median grain size to channel slope. Our sediment budget suggests a total sediment efflux from Washington Creek catchment of 61 × 10³ m³. There are similar contributions from hillslopes (43 ± 14 × 10³ m³), first to third stream order channel (35 ± 12 × 10³ m³) and the arterial fourth to fifth stream order channel (31 ± 17 × 10³ m³) to the total volume of erosion. Deposition (39 ± 17 × 10³ m³) within the arterial channel was higher than erosion (31 ± 17 × 10³ m³), which means a net sediment gain of about 8 × 10³ m³ in the arterial channel. The ratio of total deposition to total erosion was 0.44. For fines <63 μm, this ratio was much smaller (0.11), which means that fines are preferentially exported. This has important implications for suspended sediment and water quality in downstream rivers. © 2019 John Wiley & Sons, Ltd.     

Application of <Superscript>137</Superscript>Cs fingerprinting technique to interpreting sediment production records from reservoir deposits in a small catchment of the Hilly Loess Plateau,China

According to variations of 137Cs and clay contents, 44 flood couplets were identified in a profile of reservoir deposit with a vertical length of 28.12 m in the Yuntaishan Gully. Couplet 27 at the middle of the profile had the highest average 137Cs content of 12.65 Bq kg-1, which indicated the 1963s' deposits, then 137Cs content decreased both downward and upward in the profile. The second top and bottom couplets had average 137Cs contents of 2.15 Bq kg-1 and 0.92 Bq kg-1, respectively. By integrated analysis of reservoir construction and management history, variations of 137Cs contents over the profile, sediment yields of flood couplets and rainfall data during the period of 1958-1970, individual storms related to the flood couplets were identified. 44 floods with a total sediment yield of 2.36×104 m3 occurred and flood events in a year varied between 1 and 10 times during the period of 1960-1970. 7-10 flood events occurred during the wet period of 1961-1964 with very wet autumn, while only 1-2 events during the dry period of 1965-1969. Average annual specific sediment yield was 1.29×104 t km-2 a-1 for the Yuntaishan Gully during the period of 1960-1970, which was slightly higher than 1.11 ×104 t km-2 a-1 for the Upper Yanhe River Basin above the Ganguyi Hydrological Station and slightly lower than 1.40 ×104 t km-2 a-1 for the nearby Zhifang Gully during the same period. Annual specific sediment yields for the Yuntaishan Gully were correlated to the wet season's rainfalls well.     

Impact of anthropogenic drivers on subaqueous topographical change in the Datong to Xuliujing reach of the Yangtze River

Changes of subaqueous topography in shallow offshore water pose safety risks for embankments,navigation,and ports.This study conducted measurements of subaqueous topography between Datong and Xuliujing in the Yangtze River using a Sea Bat 7125 multi-beam echo sounder,and the channel change from 1998 to 2013 was calculated using historical bathymetry data.The study revealed several important results:(1)the overall pattern of changes through the studied stretch of the river was erosion–deposition–erosion.Erosion with a volume 700×10~6m~3occurred in the upper reach,deposition of about 204×10~6m~3occurred in the middle reach,and erosion of about 602×10~6m~3occurred in the lower reach.(2)Dunes are the most common microtopographic feature,accounting for 64.3%of the Datong to Xuliujing reach,followed by erosional topography and flat river topography,accounting for 27.6%and 6.6%,respectively.(3)Human activities have a direct impact on the development of the microtopography.For instance,the mining of sand formed holes on the surface of dunes with lengths of 20–35 m and depths of 3–5 m.We concluded that the overall trend of erosion(net erosion volume of 468×10~6m~3)occurred in the study area mainly because of the decreased sediment discharge following the closure of the Three Gorges Dam.However,other human activities were also impact factors of topographic change.Use of embankments and channel management reduced channel width,restricted river meandering,and exacerbated the erosion phenomenon.     

Finite-Amplitude thermal convection within a self-gravitating fluid sphere

Abstract

Finite-difference calculations have been carried out to determine the structure of finite-amplitude thermal convection within a self-gravitating fluid sphere with uniform heat release. For a fixed-surface boundary condition single-cell convection breaks up into double-cell convection at a Rayleigh number of 3 × 10⁴, at a Rayleigh number of 5 × 10⁵ four-cell convection is observed. With a free-surface boundary condition only single cell convection is obtained up to a Rayleigh number of 5 × 10⁶.     

Palaeoflood hydrology of the Salt river,Arizona

Geomorphic evidence along bedrock-confined reaches of the Salt River in east-central Arizona provides a record of the river's largest flood events. Fine-grained flood slackwater deposits accumulated at channel margin irregularities several metres above the low-flow channel. Discharges associated with flow events responsible for the deposits were estimated by computer flow modelling. These estimates document flood magnitudes in excess of gauged historic streamflows. Relative and radiocarbon dating suggest that a flood record in excess of 600 y is preserved in the slackwater sequences. A prominent flood scar cut into grussy hillslope soils allows the extension of the prehistoric flood record to several thousand years. A maximum discharge estimate of 4600 m³s^?1 affixed to the flood scar represents the largest flood event in the record, and is given a minimum recurrence interval of 1000–2000 y. The 1952 flood is the largest historic flow event experienced along the study reach and is estimated at 2900 m³s^?1. Two palaeoflood events preserved in the slackwater stratigraphy exceed the 1952 event, and are given recurrence intervals of 300 and 600 y. The latter flood event had an estimated discharge of 3200 m³s^?1. It is apparent that discharge estimates affixed to these infrequent, large-magnitude flood events approach a maximum with decreased probabilities (large recurrence intervals). This suggests that a physical limit on discharge may exist within the Salt River drainage basin and is perhaps directly related to drainage basin size.     

Leakage properties of the multi-aquifer system underlying the Greater Onitsha water scheme borehole field,Anambra State,Nigeria

Abstract

Leakage properties and the potential for land subsidence due to groundwater withdrawal from a multi-aquifer water supply system were investigated by applying leaky type curve and one dimensional consolidation models to drawdown data that were obtained during a pumping test experiment in an aquifer-aquitard system. The producing aquifer has transmissivity and storativity values of 5.3 × 10^?3 m² s^?1 and 9.54 × 10^?4 respectively. It is recharged through leakage at a rate of 5.67 × 10^?8 m s^?1, giving a leakage amount of more than 0.007 m³ s^?1. Drainage of the aquifer-aquitard system could result in aquitard compaction of between 50 and 180 mm year^?1 for pumping periods of 6 and 22 h day^?1, respectively. The observed leakage has important implications for land subsidence problems and waste disposal practices in the area.     

SUR LES PROJECTS FUTURS DE L'AISH / ON THE FUTURE PLANS OF IAHS

Abstract

The baseflow characteristics of some of the numerous small basins in southeastern Nigeria have been analysed to estimate the developable groundwater in the basins. It is shown that from 5.62 × 10⁴ to 1.59 × 10⁶ m³ of groundwater can be developed per square kilometre of basin per annum. The relationship between the baseflow characteristics and other attributes of the basins, such as geology and stream density, were studied statistically, leading to the development of empirical equations for predicting the hydrological features of the several ungauged streams in the region. It is shown, for example, that the basin geology (represented as the percentage of sands), the drainage density, the basin area, the baseflow depletion rate and the total groundwater stored in the basin, Q_tp, are related by the equation:

Q_tp = ?1.85 × 10⁹?7.96 × 10⁸ dd+4.18 × 10⁷ gf?2.01 × 10⁶ df+6.25 × 10⁵ wa

where dd is drainage density; gf geological factor; df depletion factor; and wa basin area.     

Yields of suspended sediment and dissolved solids from the Andean basins of Ecuador

Abstract

Water discharge and suspended and dissolved sediment data from three rivers (Napo, Pastaza and Santiago) in the Ecuadorian Amazon basin and a river in the Pacific basin (Esmeraldas) over a 9-year period, are presented. This data set allows us to present: (a) the chemical weathering rates; (b) the erosion rates, calculated from the suspended sediment from the Andean basin; (c) the spatio-temporal variability of the two regions; and (d) the relationship between this variability and the precipitation, topography, lithology and seismic activity of the area. The dissolved solids load from the Esmeraldas basin was 2 × 10⁶ t year^-1, whereas for the Napo, Pastaza and Santiago basins, it was 4, 2 and 3 × 10⁶ t year^-1, respectively. For stations in the Andean piedmont of Ecuador, the relationship between surface sediment and the total sediment concentration was found to be close to one. This is due to minimal stratification of the suspended sediment in the vertical profile, which is attributed to turbulence and high vertical water speeds. However, during the dry season, when the water speed decreases, sediment stratification appears, but this effect can be neglected in the sediment flux calculations due to low concentration rates. The suspended sediment load in the Pacific basin was 6 × 10⁶ t year^-1, and the total for the three Amazon basins was 47 × 10⁶ t year^-1. The difference between these contributions of the suspended sediment load is likely due to the tectonic uplift and the seismic and volcanic dynamics that occur on the Amazon side.

Editor Z.W. Kundzewicz

Citation Armijos, E., Laraque, A., Barba, S., Bourrel, L., Ceron, C., Lagane, C., Magat, P., Moquet, J.-S., Pombosa, R., Sondag, F., Vauchel, P., Vera, A., and Guyot, J.L., 2013. Yields of suspended sediment and dissolved solids from the Andean basins of Ecuador. Hydrological Sciences Journal, 58 (7), 1478–1494.     

Catchment-scale storm velocity: quantification,scale dependence and effect on flood response

Abstract

The concept of “catchment-scale storm velocity” quantifies the rate of storm motion up and down the basin accounting for the interaction between the rainfall space–time variability and the structure of the drainage network. It provides an assessment of the impact of storm motion on flood shape. We evaluate the catchment-scale storm velocity for the 29 August 2003 extreme storm that occurred on the 700 km²-wide Fella River basin in the eastern Italian Alps. The storm was characterized by the high rate of motion of convective cells across the basin. Analysis is carried out for a set of basins that range in area from 8 to 623 km² to: (a) determine velocity magnitudes for different sub-basins; (b) examine the relationship of velocity with basin scale and (c) assess the impact of storm motion on simulated flood response. Two spatially distributed hydrological models of varying degree of complexity in the representation of the runoff generation processes are used to evaluate the effects of the storm velocity on flood modelling and investigate model dependencies of the results. It is shown that catchment-scale storm velocity has a non-linear dependence on basin scale and generally exhibits rather moderate values, in spite of the strong kinematic characteristics of individual storm elements. Consistently with these observations and for both models, hydrological simulations show that storm motion has an almost negligible effect on the flood response modelling.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Nikolopoulos, E.I., Borga, M., Zoccatelli, D., and Anagnostou, E.N., 2014. Catchment-scale storm velocity: quantification, scale dependence and effect on flood response. Hydrological Sciences Journal, 59 (7), 1363–1376. http://dx.doi.org/10.1080/02626667.2014.923889     

Quantification of the hydraulic diffusivity of a bentonite-sand mixture using the water head decrease measured upon sudden flow interruption

This paper presents a new modelling approach to quantify the hydraulic diffusivity of low-permeability unconsolidated porous media under confined saturated-flow conditions in the laboratory. The derived analytical solution for the transient variation of the hydraulic head after flow interruption was applied to experimental data obtained from continuous measurements of the water pressure at two locations in the soil column. Three soil samples made of a mixture of natural bentonite (at different mass fractions) and medium sand were studied during a series of stepwise constant flow rates of water. The numerical results well fit the experimentally measured decrease of the dimensionless hydraulic head. The study shows that the increase of the mass fraction of bentonite in the soil sample from 10 to 30% is accompanied by a strong decrease of the hydraulic diffusivity from 2.4 × 10⁻² to 1.1 × 10⁻³ m² s⁻¹, which is clearly due to the decrease of the hydraulic conductivity of the soil sample. The specific storages obtained for each of the three samples are in the same order of magnitude and seem to decrease with the increase of mass fraction of bentonite. However, they clearly reflect the predominant portion of the compressibility of the porous medium compared with that of water. Compared with reported literature values for clayey soils, the specific storage values in this study are slightly higher, varying within the range of 2 × 10⁻³ to 8.1 × 10⁻³ m⁻¹._. The experimental results also give insight into the limitations of the modelling approach. In the case of low-permeability soils (K < 2 × 10⁻⁶ ms⁻¹) and steady-flow conditions with low Reynolds numbers, for example, Re < 0.003, it is recommended to choose a time duration for flow interruption between subsequent flow rate steps of longer than 5 s. For high-permeability porous media, to increase the precision of the quantified hydraulic diffusivity, it might be useful to select a measuring frequency significantly higher than 1 Hz.     

Influence of storm-related sediment storage on the sediment delivery from tributary catchments in the upper Waipaoa River,New Zealand

Although much is known about overall sediment delivery ratios for catchments as components of sediment production and sediment yield, little is known about the component of temporary sediment storage. Sediment delivery ratios focused on the influence of storm-related sediment storage are measured at Matakonekone and Oil Springs tributaries of the Waipaoa River basin, east coast of New Zealand. The terrace deposits of both tributaries show abundant evidence of storm-related sedimentation, especially sediment delivered from Cyclone Bola, a 50 year return rainfall event which occurred in 1988. The sediment delivery ratio is calculated by dividing the volume of sediment transported from a tributary to the main stream by the volume of sediment generated at erosion sites in the tributary catchment. Because the sediment delivery volume is unknown, it can be calculated as the difference between sediment generation volume and sediment storage volume in the channel reach of the tributary. The volume of sediment generated from erosion sites in each tributary catchment was calculated from measurements made on aerial photographs dating from 1960 (1:44 000) and 1988 (1:27 000). The volume of sediment stored in the tributary can be calculated from measurements of cross-sections located along the tributary channel, which are accompanied by terrace deposits dated by counting annual growth rings of trees on terrace surfaces. Sediment delivery ratios are 0·93 for both Matakonekone catchment and Oil Springs catchment. Results indicate that Oil Springs catchment has contributed more than twice the volume of sediment to the Waipaoa River than the Matakonekone catchment (2·75 × 10⁶ m³ vs 1·22 × 10⁶ m³). Although large volumes of sediment are initially deposited during floods, subsequent smaller flows scour away much of these deposits. The sediment scouring rate from storage is 1·25 × 10⁴ m³ a⁻¹ for Matakonekone stream and 0·83 × 10⁴ m³ a⁻¹ for Oil Springs stream. Matakonekone and Oil Springs channels respond to extreme storms by instantaneously aggrading, then gradually excavating the temporarily stored sediment. Results from Matakonekone and Oil Springs streams suggest a mechanism by which event recurrence interval can strongly influence the magnitude of a geomorphic change. Matakonekone stream with its higher stream power is expected to excavate sediment deposits more rapidly and allow more rapid re-establishment of storage capacity. Copyright © 1999 John Wiley & Sons, Ltd.     

Genetic programming approach for flood routing in natural channels

In recognition of the non‐linear relationship between storage and discharge existing in most river systems, non‐linear forms of the Muskingum model have been proposed, together with methods to calibrate the model parameters. However, most studies have focused only on routing a typical hypothetical flood hydrograph characterized by a single peak. In this study, we demonstrate that the storage–discharge relationship adopted for the non‐linear Muskingum model is not adequate for routing flood hydrographs in natural channels, which are often characterized by multiple peaks. As an alternative, an evolutionary algorithm‐based modelling approach, i.e. genetic programming (GP), is proposed, which is found to route complex flood hydrographs accurately. The proposed method is applied for constructing a routing model for a channel reach along the Walla Walla River, USA. The GP model performs extremely well with a root‐mean‐square error (RMSE) of 0·73 m³ s^?1 as against an RMSE of 3·26 m³ s^?1 for routing the multi‐peaked hydrograph. The advantage of GP lies in the fact that, unlike other models, it establishes the routing relationship in an easy and simple mathematical form. Copyright © 2007 John Wiley & Sons, Ltd.     

Characterizing replenishment water,lake water and groundwater interactions by numerical modelling in arid regions: a case study of Shahu Lake

In order to maintain the scenic and eco-environmental values of a lake, we need to characterize its water interactions. Shahu Lake was used as a case study to show the interactions among replenishment water, lake water and groundwater in an arid region. Shahu Lake is located in the Ningxia Hui Autonomous Region of northwest China and has an area of 13.96 km² and an average depth of 2.2 m. The groundwater modelling software MODFLOW was used. The analysis results show that hydraulic connectivity among replenishment water, lake water and groundwater is the crucial driving factor that affects the water level in Shahu Lake. The lake water level is highly sensitive to the volume of replenishment water. The groundwater is of great importance in balancing the water level in the lake and preventing it from drying up. It was determined that 13.8 × 10⁶ m³/yr is the optimal volume of replenishment water for Shahu Lake in order to maintain the lake level at its normal state and also to make the best use of available water resources on a long-term basis. Understanding of the water interactions can promote effective management of water resources in Shahu Lake.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR D. Hughes