Flux difference splitting for open channel flows

A finite difference scheme based on flux difference splitting is presented for the solution of the one-dimensional shallow-water equations in open channels, together with an extension to two-dimensional flows. A linearized problem, analogous to that of Riemann for gas dynamics, is defined and a scheme, based on numerical characteristic decomposition, is presented for obtaining approximate solutions to the linearized problem. The method of upwind differencing is used for the resulting scalar problems, together with a flux limiter for obtaining a second-order scheme which avoids non-physical, spurious oscillations. The scheme is applied to a one-dimensional dam-break problem, and to a problem of flow in a river whose geometry induces a region of supercritical flow. The scheme is also applied to a two-dimensional dam-break problem. The numerical results are compared with the exact solution, or other numerical results, where available.

A numerical scheme for two-dimensional, open channel flows with non-rectangular geometries

An algorithm based on flux difference splitting is presented for the solution of two-dimensional, open channel flows. A transformation maps a non-rectangular, physical domain into a rectangular one. The governing equations are then the shallow water equations, including terms of slope and friction, in a generalized coordinate system. A regular mesh on a rectangular computational domain can then be employed. The resulting scheme has good jump capturing properties and the advantage of using boundary/body-fitted meshes. The scheme is applied to a problem of flow in a river whose geometry induces a region of supercritical flow.

Analytic Benchmark Solutions for Open-Channel Flows

An important test of the quality of a computational model is its ability to reproduce standard test cases or benchmarks. For steady open–channel flow based on the Saint Venant equations some benchmarks exist for simple geometries from the work of Bresse, Bakhmeteff and Chow but these are tabulated in the form of standard integrals. This paper provides benchmark solutions for a wider range of cases, which may have a nonprismatic cross section, nonuniform bed slope, and transitions between subcritical and supercritical flow. This makes it possible to assess the underlying quality of computational algorithms in more difficult cases, including those with hydraulic jumps. Several new test cases are given in detail and the performance of a commercial steady flow package is evaluated against two of them. The test cases may also be used as benchmarks for both steady flow models and unsteady flow models in the steady limit.

Flow controls on lowland river macrophytes: A review

We review the current status of knowledge regarding the role that flow parameters play in controlling the macrophyte communities of temperate lowland rivers. We consider both direct and indirect effects and the interaction with other factors known to control macrophyte communities. Knowledge gaps are identified and implications for the management of river systems considered. The main factors and processes controlling the status of macrophytes in lowland rivers are velocity (hence also discharge), light, substrate, competition, nutrient status and river management practices. We suggest that whilst the characteristics of any particular macrophyte community reflect the integral effects of a combination of the factors, fundamental importance can be attributed to the role of discharge and velocity in controlling instream macrophyte colonisation, establishment and persistence. Velocity and discharge also appear to control the relative influence of some of the other controlling factors. Despite the apparent importance of velocity in determining the status of macrophyte communities in lowland rivers, relatively little is understood about the nature of the processes controlling this relationship. Quantitative knowledge is particularly lacking. Consequently, the ability to predict macrophyte abundance and distribution in rivers is still limited. This is further complicated by the likely existence of feedback effects between the growth of macrophytes and velocity. Demand for water resources increases the pressure on lowland aquatic ecosystems. Despite growing recognition of the need to allocate water for the needs of instream biota, the inability to assess the flow requirements of macrophyte communities limits the scope to achieve this. This increases the likelihood of overexploitation of the water resource as other users, whose demands are quantifiable, are prioritised. (C) 2008 Elsevier B.V. All rights reserved.

Integration of remote sensing and GIS for modelling flash floods in Wadi Hudain catchment, Egypt

This paper presents a new approach to modelling flash floods in dryland catchments by integrating remote sensing and digital elevation model (DEM) data in a geographical information system (GIS). The spectral reflectance of channels affected by recent flash floods exhibit a marked increase, due to the deposition of fine sediments in these channels as the flood recedes. This allows the parts of a catchment that have been affected by a recent flood event to be discriminated from unaffected parts, using a time series of Landsat images. Using images of the Wadi Hudain catchment in southern Egypt, the hillslope areas contributing flow were inferred for different flood events. The SRTM3 DEM was used to derive flow direction, flow length, active channel cross-sectional areas and slope. The Manning Equation was used to estimate the channel flow velocities, and hence the time-area zones of the catchment. A channel reach that was active during a 1985 runoff event, that does not receive any tributary flow, was used to estimate a transmission loss rate of 7·5 mm h−1, given the maximum peak discharge estimate. Runoff patterns resulting from different flood events are quite variable; however the southern part of the catchment appears to have experienced more floods during the period of study (1984–2000), perhaps because the bedrock hillslopes in this area are more effective at runoff production than other parts of the catchment which are underlain by unconsolidated Quaternary sands and gravels. Due to high transmission loss, runoff generated within the upper reaches is rarely delivered to the alluvial fan and Shalateen city situated at the catchment outlet. The synthetic GIS-based time area zones, on their own, cannot be relied on to model the hydrographs reliably; physical parameters, such as rainfall intensity, distribution, and transmission loss, must also be considered.

Holocene settlement shifts and paleoenvironments on the Central Iranian Plateau: investigating linked systems

For thousands of years, humans have inhabited locations that are highly vulnerable to the impacts of climate change, earthquakes, and floods. In order to investigate the extent to which Holocene environmental changes may have impacted on cultural evolution, we present new geologic, geomorphic, and chronologic data from the Qazvin Plain in northwest Iran that provides a backdrop of natural environmental changes for the simultaneous cultural dynamics observed on the Central Iranian Plateau. Well-resolved archaeological data from the neighbouring settlements of Zagheh (7170—6300 yr BP), Ghabristan (6215—4950 yr BP) and Sagzabad (4050—2350 yr BP) indicate that Holocene occupation of the Hajiarab alluvial fan was interrupted by a 900 year settlement hiatus. Multiproxy climate data from nearby lakes in northwest Iran suggest a transition from arid early-Holocene conditions to more humid middle-Holocene conditions from c. 7550 to 6750 yr BP, coinciding with the settlement of Zagheh, and a peak in aridity at c. 4550 yr BP during the settlement hiatus. Palaeoseismic investigations indicate that large active fault systems in close proximity to the tell sites incurred a series of large (MW ~7.1) earthquakes with return periods of ~500—1000 years during human occupation of the tells. Mapping and optically stimulated luminescence (OSL) chronology of the alluvial sequences reveals changes in depositional style from coarse-grained unconfined sheet flow deposits to proximal channel flow and distally prograding alluvial deposits sometime after c. 8830 yr BP, possibly reflecting an increase in moisture following the early-Holocene arid phase. The coincidence of major climate changes, earthquake activity, and varying sedimentation styles with changing patterns of human occupation on the Hajiarab fan indicate links between environmental and anthropogenic systems. However, temporal coincidence does not necessitate a fundamental causative dependency.

Water quality trends in the Windrush catchment: Nitrogen speciation and sediment interactions

For the predominantly agricultural River Windrush catchment, spatial variations in concentrations of nitrogen species and suspended sediment were strongly related to geology and land use. Temporal patterns of NO3- and NO2- concentrations during the three year study were highly correlated with seasonal variations in baseflow. Suspended sediment concentrations were mainly controlled by storm discharge. Variations in total ammonium concentrations reflected both flow controls. Suspended sediment effects total ammonium and organic nitrogen transport to the aquatic system, and in-stream cycling processes. Organic nitrogen did not display consistent seasonal variations, but concentrations occasionally exceeding those of NO3-. Overall, NO3- and organic nitrogen were the most important at 60% and -40%, of total nitrogen load, respectively. Future assessments of agriculture impact on river water quality should consider the total nitrogen load, and not solely that of NO3-.

Physico-chemical controls on phosphorus cycling in two lowland streams. Part 1 - the water column

This paper investigates phosphorus (P) transport and transformation dynamics in two contrasting sub-catchments of the River Kennel, England. Samples were collected daily under baseflow and hourly under stormflow conditions using autosamplers for 2 years and analysed for a range of determinands (full P fractionation, suspended sediment (SS), cations, pH, alkalinity, temperature and oxygen). Concentrations of SRP, SUP, PP and SS were higher in the flashy River Enborne (means of 0.186, 0.071, 0.101 and 34 mg l(-1), respectively) than the groundwater-fed River Lambourn (0.079, 0.057, 0.028 and 9 mg l(-1), respectively). A seasonal trend in the daily P dataset was evident, with lower concentrations during intermediate flows and the spring (caused by a dilution effect and macrophyte uptake) than during baseflow conditions. However, in the hourly P dataset, highest concentrations were observed during storm events in the autumn and winter (reflecting higher scour with increased capacity to entrain particles). Storm events were more significant in contributing to the total P load in the River Enborne than the River Lambourn, especially during August to October, when dry antecedent conditions were observed in the catchment. Re-suspension of P-rich sediment that accumulated within the channel during summer low flows might account for these observations. It is suggested that a P-calcite co-precipitation mechanism was operating during summer in the River Lambourn, while adsorption by metal oxyhydroxide groups was an important mechanism controlling P fractionation in the River Enborne. The influence of flow conditions and channel storage/release mechanisms on P dynamics in these two lowland rivers is assessed. (C) 2004 Elsevier B.V. All rights reserved.

Salt marsh morphodynamics: An investigation of tidal flows and marsh channel equilibrium

Although numerous field studies have evaluated flow and transport processes in salt marsh channels, the overall role of channels in delivering and removing material from salt marsh platforms is still poorly characterised. In this paper, we consider this issue based on a numerical hydrodynamic model for a prototype marsh system and on a field survey of the cross-sectional geometry of a marsh channel network. Results of the numerical simulations indicate that the channel transfers approximately three times the volume of water that would be estimated from mass balance considerations alone. Marsh platform roughness exerts a significant influence on the partitioning of discharge between the channel and the marsh platform edge, alters flow patterns on the marsh platform due to its effects on channel-to-platform transfer and also controls the timing of peak discharge relative to marsh-edge overtopping. Although peak channel discharges and velocities are associated with the flood tide and marsh inundation, a larger volume of water is transferred by the channel during ebb flows, a portion of which transfer takes place after the tidal height is below the marsh platform. Detailed surveys of the marsh channels crossing a series of transects at Upper Stiffkey Marsh, north Norfolk, England, show that the total channel cross-sectional area increases linearly with catchment area in the inner part of the marsh, which is consistent with the increase in shoreward tidal prism removed by the channels. Toward the marsh edge, however, a deficit in the total cross-sectional area develops, suggesting that discharge partitioning between the marsh channels and the marsh platform edge may also be expressed in the morphology of marsh channel systems.

Extraction of tidal channel networks from airborne scanning laser altimetry and aerial photography

The study of the morphodynamics of tidal channel networks is important because of their role in tidal propagation and the evolution of salt-marshes and tidal flats. Channel dimensions range from tens of metres wide and metres deep near the low water mark to only 20-30cm wide and 20cm deep for the smallest channels on the marshes. The conventional method of measuring the networks is cumbersome, involving manual digitising of aerial photographs. This paper describes a semi-automatic knowledge-based network extraction method that is being implemented to work using airborne scanning laser altimetry (and later aerial photography). The channels exhibit a width variation of several orders of magnitude, making an approach based on multi-scale line detection difficult. The processing therefore uses multi-scale edge detection to detect channel edges, then associates adjacent anti-parallel edges together to form channels using a distance-with-destination transform. Breaks in the networks are repaired by extending channel ends in the direction of their ends to join with nearby channels, using domain knowledge that flow paths should proceed downhill and that any network fragment should be joined to a nearby fragment so as to connect eventually to the open sea.