Effects of waves on the initiation of headland-associated sandbanks

Linear sandbanks appear in the lee of coastal headlands where the hydrodynamics are dominated by strong tidal currents and the seabed is characterized by an abundance of sands. They may develop as symmetrical sandbanks on either sides of the headland or as an unique banner bank. The present study numerically investigates the combined effects of waves and tide on the initial development of headland-associated sandbanks. A morphological model based on the coupling of the wave propagation module SWAN (Simulating WAves Nearshore) with the three-dimensional circulation module COHERENS (COupled Hydrodynamical-Ecological model for REgioNal and Shelf seas) is applied to an idealized Gaussian shaped headland for waves conditions varying in heights and directions at the offshore boundary. The coupling considers the effects of the interactions between the wave and current bottom boundary layers, namely the enhanced levels of turbulence near the bottom and the increase of the total bottom shear stress. Waves substantially modify the initial development of sandbanks formed by suspension narrowing their width and reorienting them along the side of the headland. They weakly impact the morphogenesis of sandbanks by bedload favoring on a short-time scale the growth of symmetric circular-shaped features and a central depositional spit prolonging the headland tip. Waves of transverse directions toward the tip of the headland contribute to the initiation by suspension of a well-developed feature in the headland side of low energy limiting the seabed evolution in the exposed area.     

Water and sediment movement around a coastal headland: Portland Bill,southern UK

Numerical simulations of tidal flow and sand transport around a coastal headland (Portland Bill, southern UK) were undertaken to investigate patterns of sand transport during the development of tidally induced transient eddies. Results obtained from a 2-D finite-element hydrodynamic model (TELEMAC-2D) were combined with a sediment transport model (SEDTRANS), to simulate the sand transport processes around the headland. Simulation of the tidal flow around Portland Bill has shown the formation and evolution of tidally induced transient eddies, around the headland. During the evolution of these transient eddies, no current-induced bedload (transport) eddy is formed for either side of the headland. Net bedload sand transport direction, around a coastal headland, is the result of instantaneous gradients in bedload transport rates, during flood and ebb flows, rather than the average (residual) flow. Thus, the use of residual (water) circulation to describe patterns of sediment movement as bedload is not an appropriatedapproach. In the case study presented here, the distinct characteristics of the coastal and seabed morphology around the Isle of Portland (i.e. headland shape and the bathymetry) indicate that these parameters can be influencing tidal (flow) and sediment dispersion around the headland. Such an interpretation has broader implications and applications to headland-associated sandbanks elsewhere.Responsible Editor: Hans Burchard     

Post‐fire thin debris flows: sediment transport and numerical modelling

The creation of a hydrophobic layer in the soil during ?res in semi‐arid environments inhibits the in?ltration of rainfall. This leads to increased rates of runoff and associated sediment transport. When the hydrophobic layer is deposited beneath the soil surface, a perched water table develops which may cause thin (1–2 cm) hillslope failures that are distinguishable from features caused by rilling and sheet?ow. Evidence for these failures was observed after a ?re near Santa Barbara, California. The amount of sediment eroded from some hillslopes was substantial, with 290 kg of sediment per metre width of hillslope delivered to the valley ?oor. The mechanics of these failures are examined with a numerical model that incorporates a stability analysis with subsurface ?ow routing along a typical hillslope pro?le. The model correctly predicts the location of the failures as well as the rainfall amount necessary to trigger them. Copyright © 2003 John Wiley & Sons, Ltd.     

Secondary flows and the pool-riffle unit: A case study of the processes of meander development

The major bedforms of gravel bed rivers, whether braided, meandering or straight, may be usefully resolved into pool-riffle units, each comprising a single scour pool together with an associated depositional shoal downstream. At low flow, the latter may be characterized by a single emergent bar-head and submerged riffle, or by a variety of remnant braid bars, depending upon the type of channel. Identification of pool-riffle units and observation of associated flow structures on a small meandering stream in northwest England has demonstrated important links between bedforms, flow patterns and channel change. Each unit appears to be associated with a systematic pattern of secondary flows, which are able to modify the bedforms and initiate meander development. Feedback links between plan morphology, flow patterns and erosive and depositional forces within these units ensure that each stage of meander growth has a characteristic style of channel change. Consequently, meanders tend to evolve by regular cycles of increasing curvature and complexity.     

Zonal shear flows with a free surface: Hamiltonian formulation and linear and nonlinear stability

General theoretical results via a Hamiltonian formulation are developed for zonal shear flows with the inclusion of the vortex stretching effect of the deformed free surface (equivalent barotropic model). These results include a generalization of the Flierl–Stern–Whitehead zero angular momentum theorem for localized nonlinear structures (whether or not on a β-plane), and sufficient conditions for linear and nonlinear stability in the Liapunov sense–the latter are given as estimates in terms of an L ²-type perturbation norm which are global in time and are derived via bounds on the equilibrium potential vorticity gradient.     

Numerical modelling of the suspended sediment transport in Torres Strait

One-dimensional vertical and three-dimensional fine-resolution numerical models of sediment transport have been developed and applied to the Torres Strait region of northern Australia. The one-dimensional model, driven by measured waves and currents, was calibrated against measured suspended sediment concentrations using a sequential data assimilation algorithm. The algorithm produced a good match between model and data, but this was achieved only by allowing some temporal variability in parameter values, suggesting that there were underlying uncertainties in the model structure and forcing data. Implications of the assimilation results to the accuracy of the numerical modelling are discussed and the need for observational programmes having an extensive spatial and temporal coverage is highlighted. The three-dimensional sediment model, driven by modelled waves and currents, simulates sediment transport over the shelf during the monsoon and trade-wind seasons covering 1997–2000. The model predicts strong seasonal variability of the sediment transport on the shelf attributed to seasonally varying hydrodynamics, and illustrates significant inter-annual variability of the sediment fluxes driven by extreme events. The developed model provides a platform for testing scientific hypothesis. With additional calibration, including uncertainty analysis, it can also be used in a management context.     

Aeolian processes across transverse dunes. II: modelling the sediment transport and profile development

This paper discusses a model which simulates dune development resulting from aeolian saltation transport. The model was developed for application to coastal foredunes, but is also applicable to sandy deserts with transverse dunes. Sediment transport is calculated using published deterministic and empirical relationships, describing the influence of meteorological conditions, topography, sediment characteristics and vegetation. A so-called adaptation length is incorporated to calculate the development of transport equilibrium along the profile. Changes in topography are derived from the predicted transport, using the continuity equation. Vegetation height is incorporated in the model as a dynamic variable. Vegetation can be buried during transport events, which results in important changes in the sediment transport rates. The sediment transport model is dynamically linked to a second-order closure air flow model, which predicts friction velocities over the profile, influenced by topography and surface roughness. Modelling results are shown for (a) the growth and migration of bare, initially sine-shaped dunes, and (b) dune building on a partly vegetated and initially flat surface. Results show that the bare symmetrical dunes change into asymmetric shapes with a slipface on the lee side. This result could only be achieved in combination with the secondorder closure model for the calculation of air flow. The simulations with the partly vegetated surfaces reveal that the resulting dune morphology strongly depends on the value of the adaptation length parameter and on the vegetation height. The latter result implies that the dynamical interaction between aeolian activity and vegetation (reaction to burial, growth rates) is highly relevant in dune geomorphology and deserves much attention in future studies. Copyright © 1999 John Wiley & Sons, Ltd.     

Significance of particle interaction to the modelling of cohesive sediment transport in rivers

Aggregation processes of fine sediments have rarely been integrated in numerical simulations of cohesive sediment transport in riverine systems. These processes, however, can significantly alter the hydrodynamic characteristics of suspended particulate matter (SPM), modifying the particle settling velocity, which is one of the most important parameters in modelling suspended sediment dynamics. The present paper presents data from field measurements and an approach to integrate particle aggregation in a hydrodynamic sediment transport model. The aggregation term used represents the interaction of multiple sediment classes (fractions) with corresponding multiple deposition behaviour. The k–ε–turbulence model was used to calculate the coefficient of vertical turbulent mixing needed for the two‐dimensional vertical‐plane simulations. The model has been applied to transport and deposition of tracer particles and natural SPM in a lake‐outlet lowland river (Spree River, Germany). The results of simulations were evaluated by comparison with field data obtained for two levels of river discharge. Experimental data for both discharge levels showed that under the prevailing uniform hydraulic conditions along the river reach, the settling velocity distribution did not change significantly downstream, whereas the amount of SPM declined. It was also shown that higher flow velocities (higher fluid shear) resulted in higher proportions of fast settling SPM fractions. We conclude that in accordance with the respective prevailing turbulence structures, typical aggregation mechanisms occur that continuously generate similar distribution patterns, including particles that settle toward the river bed and thus mainly contribute to the observed decline in the total SPM concentration. In order to determine time‐scales of aggregation and related mass fluxes between the settling velocity fractions, results of model simulations were fitted to experimental data for total SPM concentration and of settling velocity frequency distributions. The comparison with simulations for the case of non‐interacting fractions clearly demonstrated the practical significance of particle interaction for a more realistic modelling of cohesive sediment and contaminant transport. Copyright © 2004 John Wiley & Sons, Ltd.     

Calculation of primary and secondary flow and boundary shear stresses in a meandering channel

Turbulent flow in a meandering channel is computed with two Computational Fluid Dynamics (CFD) codes solving the Navier–Stokes equations by employing different turbulence closure approaches. The first CFD code solves the steady Reynolds-Averaged Navier–Stokes equations (RANS) using an isotropic turbulence closure. The second code is based on the concept of Large Eddy Simulation (LES). LES resolves the large-scale turbulence structures in the flow and is known to outperform RANS models in flows in which large-scale structures dominate the statistics. The results obtained from the two codes are compared with experimental data from a physical model study. Both, LES and RANS simulation, predict the primary helical flow pattern in the meander as well as the occurrence of an outer-bank secondary cell. Computed primary as well as secondary flow velocities are in reasonably good agreement with experimental data. Evidence is given that the outer-bank secondary cell in a meander bend is the residual of the main secondary cell of the previous bend. However, the RANS code, regardless of the turbulence model employed, overpredicts the size and strength of the outer-bank secondary cell. Furthermore, only LES is able to uphold the outer-bank second secondary cell beyond the bend apex until the exit of the bend as turbulence anisotropy contributes to its persistence. The presence of multiple secondary cells has important consequences for the distribution of shear stresses along the wetted perimeter of the channel, and thereby the sediment transport in meandering channels. Consequently, even though LES is expected to compute the bed-shear stresses along the wetted perimeter of the channel with a higher degree of accuracy than the RANS model, comparisons between LES and RANS computed wall shear stresses agree well. These findings are useful for practitioners who need to rely on RANS model predictions of the flow in meandering channels at field scale.     

Application of lattice gas techniques to the study of sediment erosion and transport caused by laminar sheetflow

In order to gain a deeper understanding of the dynamics of erosion and sediment transport on hillslopes, it seems important to clarify the role of some basic mechanisms involved in these processes. While there is evidence that this cannot be done using the theoretical framework of river hydraulics, the use of numerical analysis could be of considerable help. The nature of the problem requires a technique capable of solving Navier–Stokes equations at low Reynolds number, with geometrically complex boundaries and solid particles moving inside the flow field. All these requirements make a novel method, known as lattice gas automaton LGA, a natural candidate for the study of the hydrodynamics of sheetflows. However, due to the recent introduction of this technique, there is a lack of a clear definition of its operational limits. Considering the case of a viscous sheetflow on an erodible rough boundary, we argue that by using LGA the stream Reynolds number can be increased only at the expense of a reduction of the boundary shearing stress. Accordingly, LGA cannot profitably be used to study the beginning of sediment motion and transport. On the other hand, a further evolution of LGA, known as the lattice Boltzmann method, seems highly promising for the numerical study of the erosion processes that eventually lead to drainage network evolution along hillslopes. © 1997 by John Wiley & Sons, Ltd.     

Assessment of the impact of spatio‐temporal attributes of wetlands on stream flows using a hydrological modelling framework: a theoretical case study of a watershed under temperate climatic conditions

Wetlands play a significant role on the hydrological cycle, reducing flood peaks through water storage functions and sustaining low flows through slow water release ability. However, their impacts on water resources availability and flood control are mainly driven by wetland type (e.g. isolated wetland—IW—and riparian wetland—RW) and location within a watershed. Consequently, assessing the qualitative and quantitative impact of wetlands on hydrological regimes has become a relevant issue for scientists as well as stakeholders and decision‐makers. In this study, the distributed hydrological model, HYDROTEL, was used to investigate the role and impact of the geographic distribution of isolated and RWs on stream flows of the Becancour River watershed of the St Lawrence Lowlands, Quebec, Canada. The model was set up and calibrated using available datasets (i.e. DEM, soil, wetland distribution, climate, land cover, and hydrometeorological data for the 1969–2010 period). Different wetland theoretical location tests (WTLT) were simulated. Results were used to determine whether stream flow parameters, related to peak flows and low flows, were related to: (i) geographic location of wetlands, (ii) typology of wetlands, and (iii) seasonality. The contribution of a particular wetland was assessed using intrinsic characteristics (e.g. surface area, typology) and extrinsic factors (e.g. location in the watershed landscape and seasonality). Through these investigations, the results suggest, to some extent, that both IWs and RWs impact landscape hydrology. The more IWs are located in the upper part of the watershed, the greater their effect on both on high flow damping and low flow support seems to be. The more RWs are connected to a main stream, the greater their effect is. Our modelling results indicate that local landscape conditions may influence the wetland effect; promoting or limiting their efficiency, and thus their impacts on stream flows depend on a combined effect of wetland and landscape attributes. Copyright © 2015 John Wiley & Sons, Ltd.     

Physically based modelling of sediment generation and transport under a large rainfall simulator

A series of large rainfall simulator experiments was conducted in 2002 and 2003 on a small plot located in an experimental catchment in the North Island of New Zealand. These experiments measured both runoff and sediment transport under carefully controlled conditions. A physically based hydrological modelling system (SHETRAN) was then applied to reproduce the observed hydrographs and sedigraphs. SHETRAN uses physically based equations to represent flow and sediment transport, and two erodibility coefficients to model detachment of soil particles by raindrop erosion and overland flow erosion. The rate of raindrop erosion also depended on the amount of bare ground under the simulator; this was estimated before each experiment. These erodibility coefficients were calibrated systematically for summer and winter experiments separately, and lower values were obtained for the summer experiments. Earlier studies using small rainfall simulators in the vicinity of the plot also found the soil to be less erodible in summer and autumn. Limited validation of model parameters was carried out using results from a series of autumn experiments. The modelled suspended sediment load was also sensitive to parameters controlling the generation of runoff from the rainfall simulator plot; therefore, we found that accurate runoff predictions were important for the sediment predictions, especially from the experiments where the pasture cover was good and overland flow erosion was the dominant mechanism. The rainfall simulator experiments showed that the mass of suspended sediment increased post‐grazing, and according to the model this was due to raindrop detachment. The results indicated that grazing cattle or sheep on steeply sloping hill‐country paddocks should be carefully managed, especially in winter, to limit the transport of suspended sediment into watercourses. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling of hydrodynamics and cohesive sediment transport in Tanshui River estuarine system,Taiwan

A laterally averaged two-dimensional numerical model is used to simulate hydrodynamics and cohesive sediment transport in the Tanshui River estuarine system. The model handles tributaries as well as the main stem of the estuarine system. Observed time series of salinity data and tidally averaged salinity distributions have been compared with model results to calibrate the turbulent diffusion coefficients. The overall model verification is achieved with comparisons of residual currents and salinity distribution. The model reproduces the prototype water surface elevation, currents and salinity distributions. Comparisons of the suspended cohesive sediment concentrations calculated by the numerical model and the field data at various stations show good agreement. The validated model is applied to investigate the tidally averaged salinity distributions, residual circulation and suspended sediment concentration under low flow conditions in the Tanshui River estuarine system. The model results show that the limit of salt intrusion in the mainstem estuary is located at Hsin-Hai bridge in Tahan Stream, 26 km from the River mouth under Q75 flow. The null point is located at the head of salt intrusion, using 1 ppt isohaline as an indicator. The tidally averaged sediment concentration distribution exhibits a local maximum around the null point.     

Barchan ripples under unidirectional water flows in the laboratory: formation and planar morphology

The purpose of the present study is to investigate the formation and planar shape of barchan ripples generated by laboratory unidirectional water currents and to compare the morphology with that of barchans formed by laboratory waves and flows in natural environments. A thin veneer of sand as a sediment source was placed initially at the upstream part of a recirculating flume and the change in bed configurations by the flow was recorded by a video camera. Two types of formative process were observed: one was that barchan ripples grew from barchanoid bedforms and the other was that they developed from a small aggregate of sand particles. The barchan ripples in the present test had the same characteristics in the height–width relation and in the migration speed as previously reported from the desert environment and wave‐flume studies. An examination of the planar shape of the barchans led to the result that the barchan ripples had larger values of body‐length/width, compared with those of barchans formed by water waves in the laboratory and by airflow in natural deserts. The horn‐length/width of the barchan ripples was smaller than that of barchan dunes in deserts but larger than that of barchans in oscillatory laboratory flows. Copyright © 2005 John Wiley & Sons, Ltd.     

Predictability of river flow and suspended sediment transport in the Mississippi River basin: a non‐linear deterministic approach

As the Mississippi River plays a major role in fulfilling various water demands in North America, accurate prediction of river flow and sediment transport in the basin is crucial for undertaking both short‐term emergency measures and long‐term management efforts. To this effect, the present study investigates the predictability of river flow and suspended sediment transport in the basin. As most of the existing approaches that link water discharge, suspended sediment concentration and suspended sediment load possess certain limitations (absence of consensus on linkages), this study employs an approach that presents predictions of a variable based on history of the variable alone. The approach, based on non‐linear determinism, involves: (1) reconstruction of single‐dimensional series in multi‐dimensional phase‐space for representing the underlying dynamics; and (2) use of the local approximation technique for prediction. For implementation, river flow and suspended sediment transport variables observed at the St. Louis (Missouri) station are studied. Specifically, daily water discharge, suspended sediment concentration and suspended sediment load data are analysed for their predictability and range, by making predictions from one day to ten days ahead. The results lead to the following conclusions: (1) extremely good one‐day ahead predictions are possible for all the series; (2) prediction accuracy decreases with increasing lead time for all the series, but the decrease is much more significant for suspended sediment concentration and suspended sediment load; and (3) the number of mechanisms dominantly governing the dynamics is three for each of the series. Copyright © 2005 John Wiley & Sons, Ltd.     

The contrasted evolution of high and low flows and precipitation indices in the Duero basin (Spain)

Abstract

Trends in high and low flows are valuable indicators of hydrological change because they highlight changes in various parts of the frequency distribution of streamflow series. This enables improved assessment of water availability in regions with high seasonal and inter-annual variability. There has been a substantial reduction in water resources in the Duero basin (Iberian Peninsula, Spain) and other areas of the Mediterranean region during the last 50 years, and this is likely to continue because of climate change. In this study, we investigated the evolution and trends in high and low flows in the Spanish part of the Duero basin, and in equivalent or closely-related precipitation indices for the period 1961–2005. The results showed a general trend of decrease in the frequency and magnitude of high flows throughout most of the basin. Moreover, the number of days with low flows significantly increased over this period. No clear relationship was evident between the evolution of high/low flows and changes in the distribution frequencies of the precipitation series. In contrast to what was expected, the number of days with heavy precipitation and the mean annual precipitation did not show significant trends across the basin, and the number of days without rainfall decreased slightly. The divergence between precipitation and runoff evolution was more accentuated in spring and summer. In the absence of trends in precipitation, it is possible that reforestation processes in the region, and increasing temperatures in recent decades, could be related to the decreasing frequency of high flows and the increasing frequency of low flows.

Editor Z.W. Kundzewicz; Associate editor S. Grimaldi

Citation Morán-Tejeda, E., López-Moreno, J.I., Vicente-Serrano, S.M., Lorenzo-Lacruz, J. and Ceballos-Barbancho, A., 2012. The contrasted evolution of high and low flows and precipitation indices in the Duero basin (Spain). Hydrological Sciences Journal, 57 (4), 591–611.     

Theoretical/numerical model for the transport of non-uniform suspended sediment in open channels

In this paper, we address the transport of multi-disperse suspended sediment mixtures in open channels, via the use of the two-fluid model. To that end, we extend previously developed frameworks for the dilute and non-dilute transport of suspended sediment. Within the scope of the Reynolds-averaged Navier-Stokes (RANS) equations, these modeling frameworks comprise mass and momentum equations for both phases (water and sediment). Here, we calculate the distribution of total volumetric concentration of sediment using two approaches: (1) by considering the mixture as represented by a single size; we call this approach Partial two-fluid model for uniform sediments (PTFMU); and (2) by combining the volumetric concentration of the sediment corresponding to several particle size classes; we call this approach Partial two-fluid model for non-uniform sediments (PTFMNU). In the second approach, we propose a methodology for the computation of the overall velocity of the disperse phase as a function of the velocities of each size class. k-ε type closures to account for the turbulence in the carrier phase (water) are applied. We also consider the coupling between the two phases through the drag force. Velocities of the carrier and disperse phases, and concentrations for each sediment class size are numerically solved by integrating the differential equations over control volumes. In order to validate our models, we compare numerical results to experimental data of Einstein and Chien [H.A. Einstein, N. Chien, Effects of heavy sediment concentration near the bed on velocity and sediment distribution, MRD sediment series report, University of California, Berkley, 1955] and Taggart et al. [W.C. Taggart, C.A. Yermoli, S. Montes, A. Ippen, Effects of sediment size and gradation on concentration profiles for turbulent flow, Massachusetts Institute of Technology, 1972]. Results of mean velocity of the carrier phase are in close agreement with the experimental data. For the prediction of sediment concentrations, we observe that there is a difference in the results using the two approaches mentioned above. We additionally obtain values of the Schmidt number needed to improve the agreement between predictions of the distribution of suspended sediment and the experimental data, and discuss the effect of sediment size and increasing sediment concentration on the values of the Schmidt number.