Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Effect of saltating sediment on flow resistance and bed roughness in overland flow

The acceleration of saltating grains by overland flow causes momentum to be transferred from the flow to the grains, thereby increasing flow resistance and bed roughness. To assess the impact of saltating sediment on overland flow hydraulics, velocity profiles in transitional and turbulent flows on a fixed sand-covered bed were measured using hot-film anemometry. Five discharges were studied. At each discharge, three flows were measured: one free of sediment, one with a relatively low sediment load, and one with a relatively high sediment load. In these flows from 83 to 90 per cent of the sediment was travelling by saltation. As a result, in the sediment-laden flows the near-bed velocities were smaller and the velocity profiles steeper than those in the equivalent sediment-free flows. Sediment loads ranged up to 87·0 per cent of transport capacity and accounted for as much as 20·8 per cent of flow resistance (measured by the friction factor) and 89·7 per cent of bed roughness (measured by the ratio of the roughness length to median grain diameter). It is concluded that saltating sediment has a considerable impact on overland flow hydraulics, at least on fixed granular beds. Saltation is likely to have a relatively smaller effect on overland flow on natural hillslopes and agricultural fields where form and wave resistance dominate. Still, saltation is generally of greater significance in overland flow than in river flow, and for this reason its effect on overland flow hydraulics is deserving of further study. © 1998 John Wiley & Sons, Ltd.     

Characteristics of concentrated flow hydraulics for rangeland ecosystems: implications for hydrologic modeling

Concentrated flow is often the dominant source of water erosion following disturbance on rangelands. Because of the lack of studies that explain the hydraulics of concentrated flow on rangelands, cropland‐based equations have typically been used for rangeland hydrology and erosion modeling, leading to less accurate predictions due to different soil and vegetation cover characteristics. This study investigates the hydraulics of concentrated flow using unconfined field experimental data over diverse rangeland landscapes within the Great Basin Region, United States. The results imply that the overall hydraulics of concentrated flow on rangelands differ significantly from those of cropland rills. Concentrated flow hydraulics on rangelands are largely controlled by the amount of cover or bare soil and hillslope angle. New predictive equations for concentrated flow velocity (R² = 0·47), hydraulic friction (R² = 0·52), and width (R² = 0·4) representing a diverse set of rangeland environments were developed. The resulting equations are applicable across a wide span of ecological sites, soils, slopes, and vegetation and ground cover conditions and can be used by physically‐based rangeland hydrology and erosion models to estimate rangeland concentrated flow hydraulic parameters. Published in 2011. This article is a US Government work and is in the public domain in the USA.     

Assessment of hydrodynamic simulation results for eco‐hydraulic and eco‐hydrological applications: a spatial semivariance approach

Output from a three‐dimensional numerical flow model (SSIIM) is used in conjunction with high‐resolution topographic and velocity data to assess such models for eco‐hydraulic applications in river channel design and habitat appraisal. A new methodology for the comparison between field measurement and model output is detailed. This involves a comparison between conventional goodness‐of‐fit approaches applied to a spatially structured (riffle and pool) sample of model and field data, and a ‘relaxation’ method based upon the spatial semivariance of model/field departures. Conventional assessment indicates that the model predicts point‐by‐point velocity characteristics on a 0·45 m mesh to within ±0·1 m s⁻¹ over 80% of the channel area at low flow, and 50% of the area at high in‐bank flow. When a relative criterion of model fit is used, however, the model appears to perform less well: 60–70% of channel area has predicted velocities that depart from observed velocities by more than 10%. Regression analysis of observed and predicted velocities gives more cause for optimism, but all of these conventional indicators of goodness of fit neglect important spatial characteristics of model performance. Spatial semivariance is a means of supplementing model appraisal in this respect. In particular, using the relaxation approach, results are greatly improved: at a high in‐bank flow, the model results match field measurements to within 0·1 m s⁻¹ for more than 95% of the total channel area, provided that model and field comparisons are allowed within a radius of approximately 1 m from the original point of measurement. It is suggested that this revised form of model assessment is of particular relevance to eco‐hydraulic applications, where some degree of spatial and temporal dynamism (or uncertainty) is a characteristic. The approach may also be generalized to other environmental science modelling applications where the spatial attributes of model fits are of interest. Copyright © 2005 John Wiley & Sons, Ltd.     

The linkage between velocity patterns and sediment entrainment in a forced‐pool and riffle unit

A field‐based project was initiated in order to characterize velocities and sediment entrainment in a forced‐pool and riffle sequence. Three‐dimensional velocities and turbulence intensities were measured with an acoustic Doppler velocimeter at 222 different points at three similar flows that averaged approximately 4·35 m³ s⁻¹ within a large pool–riffle unit on North Saint Vrain Creek, Colorado. Sediment‐sorting patterns were observed with the introduction of 500 tracer particles painted according to initial seeding location. Tracer particles moved sporadically during a 113 day period in response to the annual snowmelt peak flow, which reached a maximum level of 14·8 m³ s⁻¹. Velocity data indicate high instantaneous velocities and turbulence levels in the centre of pools. Patterns of sediment deposition support the notion that stream competence is higher in the pool than the downstream riffle. Flow convergence around a large channel constriction appears to play a major role in multiple processes that include helical flow development and sediment routing, and backwater development with low velocities and turbulence levels above the constriction that may locally limit sediment supply. Jet flow, flow separation, vortex scour and turbulence generation enhance scour in the centre of pools. Ultimately, multiple processes appear to play some role in maintenance of this forced pool and the associated riffle. Copyright © 2008 John Wiley & Sons, Ltd.     

Impacts of sediment load on Manning coefficient in supercritical shallow flow on steep slopes

The Manning equation is one of the most widely used formulae for calculating the velocity of shallow overland flow in hydrological and erosion models. Precise estimation of the Manning's friction coefficient (n) is critical to determining overland flow and soil erosion processes. Few studies have been conducted to quantify the effects of sediment load on Manning's n on steep slopes. This study was conducted to investigate the potential effects of sediment load on Manning's n in a flume with a fixed bed, under wide ranges of hydraulics and sediment loads. Slope gradient varied from 8·7 to 34·2%, unit flow rate from 0·66 to 5·26 × 10^?3 m² s^?1, and sediment load from 0 to 6·95 kg m^?1 s^?1. The Reynolds number ranged from 350 to 5899. Results showed that Manning's n varied in both sediment‐free and sediment‐laden flows ranging from 0·012 to 0·055. The apparent Manning's coefficients of sediment‐laden flow were much greater than those of sediment‐free flow. The mean Manning coefficient of sediment‐laden flow was 51·27% greater than the mean value of sediment‐free flow. For sediment‐laden flow, Manning's n could be estimated with a power function of unit flow discharge and sediment content. Further studies are needed to quantify the potential effects of sediment load on the Manning's n on erodible beds and in fields. Copyright © 2010 John Wiley & Sons, Ltd.     

Hydraulic modelling of fish habitat in urban rivers during high flows

In urban rivers, flow regime and channel morphology are the drivers of physical habitat quality for aquatic species. Peak discharges are increased at high flows as a result of impermeable catchments and channel engineering for flood protection schemes. Hazardous conditions and flashy hydrographs mean that measurement of velocities at high flows is a difficult task. This research uses a three‐dimensional computational fluid dynamics (3D‐CFD) model to simulate hydraulic patterns in two urban river channels. A 3D‐CFD code, called SSIIM, was used to simulate hydraulic conditions in two engineered river reaches of the River Tame, Birmingham, UK. These two sites represent channels with different levels of engineering. Models were calibrated and tested using field measurements. Results show that modelled water surface levels and velocity profiles are well simulated. Calibrated roughness heights are compared with those derived from field measurement of sediment size. Numerical experiments are used to assess the relationship between grid resolution in the vertical dimension and the form of the modelled velocity profiles. Biologists have used laboratory experiments to determine maximum sustainable swimming speeds (MSSS) of fish, often in order to assess what level of a particular pollutant may be tolerable. In this work, simulations of high‐flow hydraulic patterns are used to compare velocity patterns with fish MSSS. Results show that when the water levels rise to fill the first channel of the two‐stage channels at the sites, which occurred 16 times in 2000, MSSS are surpassed in the majority of available habitat, suggesting that excessive velocities at high flows are one factor that limits fish habitat. A comparison between the two reaches shows that there is less available habitat in the more modified reach. Conclusions suggest that an approach that integrates water quality issues and physical channel characteristics must be taken in river rehabilitation schemes, as improvements to water quality alone may not be sufficient to improve habitat quality to the desired level. Copyright © 2002 John Wiley & Sons, Ltd.     

Two‐dimensional and three‐dimensional computational models in hydrodynamic and morphodynamic reconstructions of a river bend: sensitivity and functionality

This study assesses hydrodynamic and morphodynamic model sensitivity and functionality in a curved channel. The sensitivity of a depth‐averaged model to user‐defined parameters (grain size, roughness, transverse bed slope effect, transport relations and secondary flow) is tested. According to the sensitivity analysis, grain size, transverse bed slope effect and sediment transport relations are critical to simulated meander bend morphodynamics. The parametrization of grain size has the most remarkable effect: field‐based grain size parametrization is necessary in a successful morphodynamic reconstruction of a meander bend. The roughness parametrization method affects the distribution of flow velocities and therefore also morphodynamics. The combined effect of various parameters needs further research. Two‐dimensional (2D) and three‐dimensional (3D) reconstructions of a natural meander bend during a flood event are assessed against field measurements of acoustic Doppler current profiler and multi‐temporal mobile laser scanning data. The depth‐averaged velocities are simulated satisfactorily (differences from acoustic Doppler current profiler velocities 5–14%) in both 2D and 3D simulations, but the advantage of the 3D hydrodynamic model is unquestionable because of its ability to model vertical and near‐bed flows. The measured and modelled near‐bed flow, however, differed notably from each other's, the reason of which was left open for future research. It was challenging to model flow direction beyond the apex. The 3D flow features, which also affected the distribution of the bed shear stress, seem not to have much effect on the predicted morphodynamics: the 2D and 3D morphodynamic reconstructions over the point bar resembled each other closely. Although common features between the modelled and measured morphological changes were also found, some specific changes that occurred were not evident in the simulation results. Our results show that short‐term, sub‐bend scale morphodynamic processes of a natural meander bend are challenging to model, which implies that they are affected by factors that have been neglected in the simulations. The modelling of short‐term morphodynamics in natural curved channel is a challenge that requires further study. Copyright © 2014 John Wiley & Sons, Ltd.     

Representing the bed roughness of coarse‐grained streams in computational fluid dynamics

Computational fluid dynamics (CFD) applications are increasingly utilized for modelling complex flow patterns in natural streams and rivers. Although CFD has been successfully implemented to model many complex flow situations in natural stream settings, adequately characterizing the effects of gravel and cobble beds on flow hydraulics in CFD is a difficult challenge due to the scale of roughness lengths and the inadequacy of traditional roughness representations to characterize flow profiles in situations with large roughness elements. An alternative method of representing gravel and cobble beds is presented. Appropriate drag forces associated with different grain sizes are computed and included in the momentum equations to account for the influence of a hydraulically rough bed. Comparisons with field measurements reveal reasonable agreement between measured and modelled profiles of spatially averaged velocity and turbulent kinetic energy, and model fidelity to the non‐logarithmic behaviour of the velocity profiles. The novel method of representing coarse beds expands the utility of CFD for investigating physical processes in natural channels with large bed roughness. Copyright © 2005 John Wiley & Sons, Ltd.     

Potentials and limitations of modelling vertical distributions of root water uptake of an Austrian pine forest on a sandy soil

Root water uptake patterns are often studied with simulation models of the unsaturated soil water flow, as they are difficult to measure directly. Calibration of these models is not straightforward and causes uncertainties in simulated uptake distributions. In this paper we study how uncertainties in the calibration of the SWIF model affect uncertainty intervals in simulated uptake patterns of an Austrian pine stand (Pinus nigra var. nigra) on a sandy soil. After calibrating and validating SWIF with a large data set of more than 125 000 measured soil water contents over a three year period, uncertainty ranges in simulated soil water dynamics and root water uptake distributions were estimated with a Monte Carlo analysis. In general, uncertainties in root uptake patterns were small (typically <2 10⁻⁴ m³ m⁻³ day⁻¹) and were higher for trees with a shallow rooting system (0·8 m) than for trees with a deep rooting system (2·5 m). Uncertainties arose mainly from uncertainties in simulated soil water fluxes and from variations in the reduction of uptake during periods of drought. Uncertainties in soil water contents were far higher (typically 0·01 m³ m⁻³) than uncertainties in uptake, illustrating that uncertainties in uptake parameters and those in the distribution of water uptake hardly affect the modelling of soil water dynamics. Root water uptake models should therefore be validated against measured uptake distributions, which can be determined on sandy soils during dry periods with a high water use when soil fluxes are negligible to uptake. Copyright © 2000 John Wiley & Sons, Ltd.     

Flow modelling in gravel‐bed rivers: rethinking the bottom boundary condition

A key problem in computational fluid dynamics (CFD) modelling of gravel‐bed rivers is the representation of multi‐scale roughness, which spans the range from grain size, through bedforms, to channel topography. These different elements of roughness do not clearly map onto a model mesh and use of simple grain‐scale roughness parameters may create numerical problems. This paper presents CFD simulations for three cases: a plane bed of fine gravel, a plane bed of fine gravel including large, widely‐spaced pebble clusters, and a plane gravel bed with smaller, more frequent, protruding elements. The plane bed of fine gravel is modelled using the conventional wall function approach. The plane bed of fine gravel including large, widely‐spaced pebble clusters is modelled using the wall function coupled with an explicit high‐resolution topographic representation of the pebble clusters. In these cases, the three‐dimensional Reynolds‐averaged continuity and Navier–Stokes equations are solved using the standard k ? ε turbulence model, and model performance is assessed by comparing predicted results with experimental data. For gravel‐bed rivers in the field, it is generally impractical to map the bed topography in sufficient detail to enable the use of an explicit high‐resolution topography. Accordingly, an alternative model based on double‐averaging is developed. Here, the flow calculations are performed by solving the three‐dimensional double‐averaged continuity and Navier‐Stokes equations with the spatially‐averaged 〈k ? ε〉 turbulence model. For the plane bed of fine gravel including large, widely‐spaced pebble clusters, the model performance is assessed by comparing the spatially‐averaged velocity with the experimental data. The case of a plane gravel bed with smaller, more frequent, protruding elements is represented by a series of idealized hypothetical cases. Here, the spatially‐averaged velocity and eddy viscosity are used to investigate the applicability of the model, compared with using the explicit high‐resolution topography. The results show the ability of the model to capture the spatially‐averaged flow field and, thus, illustrate its potential for representing flow processes in natural gravel‐bed rivers. Finally, practical data requirements for implementing such a model for a field example are given. Copyright © 2011 John Wiley & Sons, Ltd.     

Shear velocity estimates in rough‐bed open‐channel flow

Shear velocity u_* is an important parameter in geophysical flows, in particular with respect to sediment transport dynamics. In this study, we investigate the feasibility of applying five standard methods [the logarithmic mean velocity profile, the Reynolds stress profile, the turbulent kinetic energy (TKE) profile, the wall similarity and spectral methods] that were initially developed to estimate shear velocity in smooth bed flow to turbulent flow over a loose bed of coarse gravel (D₅₀ = 1·5 cm) under sub‐threshold conditions. The analysis is based on quasi‐instantaneous three‐dimensional (3D) full depth velocity profiles with high spatial and temporal resolution that were measured with an Acoustic Doppler Velocity Profiler (ADVP) in an open channel. The results of the analysis confirm the importance of detailed velocity profile measurements for the determination of shear velocity in rough‐bed flows. Results from all methods fall into a range of ± 20% variability and no systematic trend between methods was observed. Local and temporal variation in the loose bed roughness may contribute to the variability of the logarithmic profile method results. Estimates obtained from the TKE and Reynolds stress methods reasonably agree. Most results from the wall similarity method are within 10% of those obtained by the TKE and Reynolds stress methods. The spectral method was difficult to use since the spectral energy of the vertical velocity component strongly increased with distance from the bed in the inner layer. This made the choice of the reference level problematic. Mean shear stress for all experiments follows a quadratic relationship with the mean velocity in the flow. The wall similarity method appears to be a promising tool for estimating shear velocity under rough‐bed flow conditions and in field studies where other methods may be difficult to apply. This method allows for the determination of u_* from a single point measurement at one level in the intermediate range (0·3 < h < 0·6). Copyright © 2013 John Wiley & Sons, Ltd.     

A simple concept for calibrating runoff thresholds in quasi‐distributed variable source area watershed models

Most semi‐distributed watershed water quality models divide the watershed into hydrologic response units (HRU) with no flow among them. This is problematic when watersheds are delineated to include variable source areas (VSAs) because it is the lateral flows from upslope areas to downslope areas that generate VSAs. Although hydrologic modellers have often successfully calibrated these types of models, there can still be considerable uncertainty in model results. In this paper, a topographic‐index‐based method is described and tested to distribute effective soil water holding capacity among HRUs, which can be subsequently adjusted using the watershed baseflow coefficient. The method is tested using a version of the Soil and Water Assessment Tool (SWAT) model that simulates VSA runoff and is applied to two watersheds: a New York State (NYS) watershed, and one in the head waters of the Blue Nile Basin (BNB) in Ethiopia. Daily streamflow predicted using effective soil water storage capacities based only on the topographic index were reassuringly accurate in both the NYS watershed (daily Nash Sutcliffe (E) = 0·73) and in the BNB (E = 0·70). Using the baseflow coefficient to adjust the effective soil water storage capacity only slightly improved streamflow predictions in NYS (E = 0·75) but substantially improved the BNB predictions (E = 0·80). By comparison, the standard SWAT model, which uses the traditional look‐up tables to determine a runoff curve number, performed considerably less accurately in un‐calibrated form (E = 0·51 for NYS and E = 0·45 for BNB), but improved substantially when explicitly calibrated to streamflow measurements (E = 0·76 for NYS and E = 0·67 for the BNB). The calibration method presented here provides a parsimonious, systematic approach to using established models in VSA watersheds that reduces the ambiguity inherent in model calibration. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of sediment load on hydraulics of overland flow on steep slopes

Eroded sediment may have significant effects on the hydraulics of overland flow, but few studies have been performed to quantify these effects on steep slopes. This study investigated the potential effects of sediment load on Reynolds number, Froude number, flow depth, mean velocity, Darcy–Weisbach friction coefficient, shear stress, stream power, and unit stream power of overland flow in a sand‐glued hydraulic flume under a wide range of hydraulic conditions and sediment loads. Slope gradients were varied from 8·7 to 34·2%, unit flow rates from 0·66 to 5·26×10^?3 m² s^?1, and sediment loads from 0 to 6·95 kg m^?1 s^?1. Both Reynolds number (Re) and Froude number (Fr) decreased as sediment load increased, implying a decrease in flow turbulence. This inverse relationship should be considered in modeling soil erosion processes. Flow depth increased as sediment load increased with a mean value of 1·227 mm, caused by an increase in volume of sediment‐laden flow (contribution 62·4%) and a decrease in mean flow velocity (contribution 37·6%). The mean flow velocity decreased by up to 0·071 m s^?1 as sediment load increased. The Darcy–Weisbach friction coefficient (f) increased with sediment load, showing that the total energy consumption increased with sediment load. The effects of sediment load on f depended on flow discharge: as flow discharge increased, the influence of sediment load on f decreased due to increased flow depth and reduced relative roughness. Flow shear stress and stream power increased with sediment load, on average, by 80·5% and 60·2%, respectively; however, unit stream power decreased by an average of 11·1% as sediment load increased. Further studies are needed to extend and apply the insights obtained under these controlled conditions to real‐world overland flow conditions. Copyright © 2010 John Wiley & Sons, Ltd.     

Overland runoff erosion dynamics on steep slopes with forages under field simulated rainfall and inflow

Although numerous studies have acknowledged that vegetation can reduce erosion, few process-based studies have examined how vegetation cover affect runoff hydraulics and erosion processes. We present field observations of overland flow hydraulics using rainfall simulations in a typical semiarid area in China. Field plots (5 × 2 m²) were constructed on a loess hillslope (25°), including bare soil plot as control and three plots with planted forage species as treatments—Astragalus adsurgens, Medicago sativa and Cosmos bipinnatus. Both simulated rainfall and simulated rainfall + inflow were applied. Forages reduced soil loss by 55–85% and decreased overland flow rate by 12–37%. Forages significantly increased flow hydraulic resistance expressed by Darcy–Weisbach friction factor by 188–202% and expressed by Manning's friction factor by 66–75%; and decreased overland flow velocity by 28–30%. The upslope inflow significantly increased overland flow velocity by 67% and stream power by 449%, resulting in increased sediment yield rate by 108%. Erosion rate exhibited a significant linear relationship with stream power. M. sativa exhibited the best in reducing soil loss which probably resulted from its role in reducing stream power. Forages on the downslope performed better at reducing sediment yield than upslope due to decreased rill formation and stream power. The findings contribute to an improved understanding of using vegetation to control water and soil loss and land degradation in semiarid environments.     

Estimating flash flood discharge in an ungauged mountain catchment with 2D hydraulic models and dendrogeomorphic palaeostage indicators

There is still wide uncertainty about past flash‐flood processes in mountain regions owing to the lack of systematic databases on former events. This paper presents a methodology to reconstruct peak discharge of flash floods and illustrates a case in an ungauged catchment in the Spanish Central System. The use of dendrogeomorphic evidence (i.e. scars on trees) together with the combined use of a two‐dimensional (2D) numerical hydraulic model and a terrestrial laser scan (TLS) has allowed estimation of peak discharge of a recent flash flood. The size and height distribution of scars observed in the field have been used to define three hypothetical scenarios (S_min or minimum scenario; S_med or medium scenario; and S_max or maximum scenario), thus illustrating the uncertainty involved in peak‐discharge estimation of flash floods in ungauged torrents. All scars analysed with dendrogeomorphic techniques stem from a large flash flood which took place on 17 December 1997. On the basis of the scenarios, peak discharge is estimated to 79 ± 14 m³ s^?1. The average deviation obtained between flood stage and expected scar height was ? 0·09 ± 0·53 m. From the data, it becomes obvious that the geomorphic position of trees is the main factor controlling deviation rate. In this sense, scars with minimum deviation were located on trees growing in exposed locations, especially on unruffled bedrock where the model predicts higher specific kinetic energy. The approach used in this study demonstrates the potential of tree‐ring analysis in palaeohydrology and for flood‐risk assessment in catchments with vulnerable goods and infrastructure. Copyright © 2010 John Wiley & Sons, Ltd.     

Modulation of the flow structure by progressive bedforms in the Kinoshita meandering channel

An in‐house fully three‐dimensional general‐purpose finite element model is applied to solve the hydrodynamic structure in a periodic Kinoshita‐generated meandering channel. The numerical model solves the incompressible Reynolds‐averaged Navier–Stokes equations for mass and momentum, while solving the k ? ε equations for turbulence. The free surface is described by the rigid‐lid approximation (using measured water surface data) for flat (smooth‐bed) and self‐formed (rough‐bed) conditions. The model results are compared against experimental measurements in the ‘Kinoshita channel’, where three‐dimensional flow velocities and turbulence parameters were measured. This validation was carried out for the upstream‐valley meander bend orientation under smooth (flat bed) conditions. After validation, several simulations were carried out to predict the hydrodynamics in conditions where either it was not possible to perform measurements (e.g. applicability of the laboratory acoustic instruments) and to extrapolate the model to other planform configurations. For the flat smooth‐bed case, a symmetric (no skewness) planform configuration was modeled and compared to the upstream‐skewed case. For the self‐formed rough‐bed case, prediction of the hydrodynamics during the progression of bedforms was performed. It appears that the presence of bedforms on a bend has the following effects: (i) the natural secondary flow of the bend is disrupted by the presence of the bedforms, thus depending on the location of the dune, secondary flows might differ completely from the traditional orientation; (ii) an increment on both the bed and bank shear stresses is induced, having as much as 50% more fluvial erosion, and thus a potential increment on the migration rate of the bend. Implications on sediment transport and bend morphodynamics are also discussed in the paper. Copyright © 2013 John Wiley & Sons, Ltd.     

Concentrated flow erosion rates reduced through biological geotextiles

Soil erosion by concentrated flow can cause serious environmental damage. Erosion‐control geotextiles have considerable potential for reducing concentrated flow erosion. However, limited data are available on the erosion‐reducing potential of geotextiles. In this study, the effectiveness of three biological geotextiles in reducing soil losses during concentrated flow is investigated. Hereto, runoff was simulated in a concentrated flow flume, filled with an erodible sandy loam on three slope gradients (13·5, 27·0 and 41·5%). Treatments included three biological geotextiles (borassus, buriti and bamboo) and one bare soil surface. Darcy–Weisbach friction coefficients ranged from 0·01 to 2·84. The highest values are observed for borassus covered soil surfaces, followed by buriti, bamboo and bare soil, respectively. The friction coefficients are linearly correlated with geotextile thickness. For the specific experimental conditions of this study, borassus geotextiles reduced soil detachment rate on average to 56%, buriti geotextiles to 59% and bamboo geotextiles to 66% of the soil detachment rate for bare soil surfaces. Total flow shear stress was the hydraulic parameter best predicting soil detachment rate for bare and geotextile covered surfaces (R² = 0·75–0·84, p <0·001, n = 12–15). The highest resistance against soil detachment was observed for the borassus covered soil surfaces, followed by buriti, bamboo and bare soil surfaces, respectively. Overall, biological geotextiles are less effective in controlling concentrated flow erosion compared with interrill erosion. Copyright © 2009 John Wiley & Sons, Ltd.