Effect of input data in hydraulic modeling for flood warning systems

The main objective of this study was to quantify the error associated with input data, including various resolutions of elevation datasets and Manning’s roughness for travel time computation and floodplain mapping. This was accomplished on the test bed, the Grand River (Ohio, USA) using the HEC-RAS model. LiDAR data integrated with survey data provided conservative predictions, whereas coarser elevation datasets provided a positive difference in the travel time (11.03–15.01%) and inundation area (32.56–44.52%). The minimum differences in travel time and inundation area were 0.50–4.33% and 3.55–7.16%, respectively, when the result from LiDAR integrated with survey data was compared with a 10-m DEM integrated with survey data. The results suggest that a 10-m DEM in the channel and LiDAR data in the floodplain combined with survey data would be appropriate for a flood warning system. Additionally, Manning’s roughness of the channel section was found to be more sensitive than that of the floodplain. The decrease in inundation area was highest (8.97%) for the lower value of Manning’s roughness.     

Numerical modeling of flow in riverine basins using an improved dynamic roughness coefficient

Various parameters such as bed and bank materials, shape and irregularity of the section, vegetation, river meanders, plan of the river path etc. affect the flow hydraulic resistance. In open channel hydraulics the effects of all these parameters are generally considered as the roughness coefficient. The Manning’s equation is one of the most practical equations to flow resistance analysis, in which the surface roughness is defined by Manning coefficient. Since many parameters are effective on the value of this coefficient, in this research study it was tried to define the roughness coefficient somehow that it be able to dynamically change with different river and hydraulic conditions. The collected data in Karun River (Iran) for two periods were used as the case study. It is shown that the accuracy of model predictions for water surface elevations were improved more than 13% in error estimation in comparison with the corresponding results obtained for a constant roughness coefficient. The roughness coefficient (n) for Karun River was also estimated using the empirical method proposed by Cowan for two different dry and wet periods. These values were then successfully compared with the average corresponding roughness coefficients calculated by the numerical model for those periods.     

Factors affecting sediment trapping in vegetated filter strips: simulation study using VFSMOD

Soil and water conservation practices have been promoted for a long time, in order to sustain agricultural activities and prevent environmental pollution. Vegetated filter strips (VFS) have been used to reduce sediment pollution into water bodies at or near the pollutant source. However, factors effecting VFS performance under natural conditions have not been well understood owing to the physical, time and financial limitations of field experiments. The use of well‐validated simulation models to understand the performance of VFS and factors affecting sediment deposition is highly justified. The objective of this research is to investigate sediment trapping in VFS and to study various factors affecting VFS performance using the simulation model VFSMOD, which was developed by researchers at University of North Carolina. Recently, VFSMOD has been validated successfully by using 21 filters with varying length, slope and vegetated cover. A wide range of five parameters was selected for the simulations, namely filter length, filter slope, manning roughness coefficient, soil type and characteristics of incoming sediment from adjacent fields. Computer simulations revealed that the length of filter is the most significant factor affecting sediment trapping in VFS. The relative increase in trapping efficiencies was not linearly related to an increase in filter length. Inflow sediment class also has a major influence on sediment trapping in VFS. The trapping efficiency of clay sediments in a 15 m length VFS was 47% compared with 92% for silt from incoming sediment. Manning roughness coefficient had a moderate effect on sediment trapping and was more significant in short filters. Land slope and soil type of VFS had a minor influence on the performance of VFS. Copyright © 2001 John Wiley & Sons, Ltd.     

Modeling river dune development and dune transition to upper stage plane bed

Large asymmetric bedforms known as dunes commonly dominate the bed of sand rivers. Due to the turbulence generation over their stoss and lee sides, dunes are of central importance in predicting hydraulic roughness and water levels. During floods in steep alluvial rivers, dunes are observed to grow rapidly as flow strength increases, undergoing an unstable transition regime, after which they are washed out in what is called upper stage plane bed. This transition of dunes to upper stage plane bed is associated with high transport of bed sediment in suspension and large decrease in bedform roughness. In the present study, we aim to improve the prediction of dune development and dune transition to upper stage plane bed by introducing the transport of suspended sediment in an existing dune evolution model. In addition, flume experiments are carried out to investigate dune development under bed load and suspended load dominated transport regimes, and to get insight in the time scales related to the transition of dunes to upper stage plane bed. Simulations with the extended model including the transport of suspended sediment show significant improvement in the prediction of equilibrium dune parameters (e.g. dune height, dune length, dune steepness, dune migration rate, dune lee side slope) both under bed load dominant and suspended load dominant transport regimes. The chosen modeling approach also allows us to model the transition of dunes to upper stage plane bed which was not possible with the original dune evolution model. The extended model predicts change in the dune shapes as was observed in the flume experiments with decreasing dune heights and dune lee slopes. Furthermore, the time scale of dune transition to upper stage plane bed was quite well predicted by the extended model. Copyright © 2015 John Wiley & Sons, Ltd.     

Evaluation of the parameters affecting the roughness coefficient of sewer pipes with rigid and loose boundary conditions via kernel based approaches

One of the important issues in water transport and sewer systems is determining the flow resistance and roughness coefficient.An accurate estimation of the roughness coefficient is a substantial issue in the design and operation of hydraulic structures such as sewer pipes,the calculation of water depth and flow velocity,and the accurate characterization of energy losses.The current study,applies two kernel based approaches [Support Vector Machine(SVM) and Gaussian Process Regression(GPR)] to develop roughness coefficient models for sewer pipes.In the modeling process,two types of sewer bed conditions were considered:loose bed and rigid bed.In order to develop the models,different input combinations were considered under three scenarios(Scenario 1:based on hydraulic characteristics,Scenarios2 and 3:based on hydraulic and sediment characteristics with and without considering sediment concentration as input).The results proved the capability of the kernel based approaches in prediction of the roughness coefficient and it was found that for prediction of this parameter in sewer pipes Scenario 3 performed better than Scenarios 1 and 2.Also,the sensitivity analysis results showed that Dgr(Dimensionless particle number) for a rigid bed and w_b/y(ratio of deposited bed width,w_b,to flow depth,y) for a loose bed had the most significant impact on the modeling process.     

The effect of asymmetric dune roughness on tidal asymmetry in the Weser estuary

The bed of estuaries is often characterized by ripples and dunes of varying size. Whereas smaller bedforms adapt their morphological shape to the oscillating tidal currents, large compound dunes (here: asymmetric tidal dunes) remain stable for periods longer than a tidal cycle. Bedforms constitute a form roughness, that is, hydraulic flow resistance, which has a large-scale effect on tidal asymmetry and, hence, on hydrodynamics, sediment transport, and morphodynamics of estuaries and coastal seas. Flow separation behind the dune crest and recirculation on the steep downstream side result in turbulence and energy loss. Since the energy dissipation can be related to the dune lee slope angle, asymmetric dune shapes induce variable flow resistance during ebb and flood phases. Here, a noncalibrated numerical model has been applied to analyze the large-scale effect of symmetric and asymmetric dune shapes on estuarine tidal asymmetry evaluated by residual bed load sediment transport at the Weser estuary, Germany. Scenario simulations were performed with parameterized bed roughness of symmetric and asymmetric dune shapes and without dune roughness. The spatiotemporal interaction of distinct dune shapes with the main drivers of estuarine sediment and morphodynamics, that is, river discharge and tidal energy, is shown to be complex but substantial. The contrasting effects of flood- and ebb-oriented asymmetric dunes on residual bed load transport rates and directions are estimated to be of a similar importance as the controls of seasonal changes of discharge on these net sediment fluxes at the Lower Weser estuary. This corroborates the need to consider dune-induced directional bed roughness in numerical models of estuarine and tidal environments.     

Parameterization of near-bed processes under collinear wave and current flows from a two-phase sheet flow model

Sediment transport models require appropriate representation of near-bed processes. We aim here to explore the parameterizations of bed shear stress, bed load transport rate and near-bed sediment erosion rate under the sheet flow regime. To that end, we employ a one-dimensional two-phase sheet flow model which is able to resolve the intrawave boundary layer and sediment dynamics at a length scale on the order of the sediment grain. We have conducted 79 numerical simulations to cover a range of collinear wave and current conditions and sediment diameters in the range 210–460 μm$\mathrm{\mu }\mathrm{m}$. The numerical results confirm that the intrawave bed shear stress leads the free stream velocity, and we assess an explicit expression relating the phase lead to the maximum velocity, wave period and bed roughness. The numerical sheet flow model is also used to provide estimates for the bed load transport rate and to inspect the near-bed sediment erosion. A common bed load transport rate formulation and two typical reference concentration approaches are assessed. A dependence of the bed load transport rate on the sediment grain diameter is observed and parameterized. Finally, the intrawave near-bed vertical sediment flux is further investigated and related to the time derivative of the bed shear stress.     

Experimental study on the effects of artificial bed roughness on turbidity currents over abrupt bed slope change

A series of experimental observations are presented in the current study to discuss the effects of artificial bed roughness on the turbidity current flowing in a rectangular channel with an abrupt change in bed slope.For this purpose,two different types of elements,sinusoidal and trapezoidal,with various heights and arrangements are considered as artificial bed roughness.A Vectrino velocity meter was used to measure the velocity and sediment concentration profiles.The effects of inlet sediment concentration on front velocity,body velocity,unit discharge,sediment concentration,and suspended load transport rate also were investigated.Accurate equations were developed for estimation of the velocity of a turbidity current over smooth and rough beds.The unexpected experimental results showed that unlike the effect of roughness height,a change in the roughness arrangement has no significant influence on the velocity of a turbidity current.Also,the effect of bed roughness on the front velocity of a denser current is more significant.     

The settling dynamics of flocculating mud and sand mixtures: part 2—numerical modelling

Estuarine and coastal sediment transport is characterised by the transport of both sand-sized particles (of diameter greater than 63?μm) and muddy fine-grained sediments (silt, diameter less than 63?μm; clay, diameter less than 2?μm). These fractions are traditionally considered as non-cohesive and cohesive, respectively, because of the negligible physico-chemical attraction that occurs between sand grains. However, the flocculation of sediment particles is not only caused by physico-chemical attraction. Cohesivity of sediment is also caused by biology, in particular the sticky extra-cellular polymeric substances secreted by diatoms, and the effect of biology in binding sediment particles can be much larger than that of physico-chemical attraction. As demonstrated by Manning (2008) and further expanded in part 1 of this paper (Manning et al., submitted), the greater binding effect of biology allows sand particles to flocculate with mud. In many estuaries, both the sand and fine sediment fractions are transported in significant quantities. Many of the more common sediment transport modelling suites now have the capability to combine mud and sand transport. However, in all of these modelling approaches, the modelling of mixed sediment transport has still essentially separated the modelling of sand and mud fractions assuming that these different fractions do not interact except at the bed. However, the use of in situ video techniques has greatly enhanced the accuracy and reliability of settling velocity measurements and has led to a re-appraisal of this widely held assumption. Measurements of settling velocity in mixed sands presented by Manning et al. (2009) have shown strong evidence for the flocculation of mixed sediments, whilst the greater understanding of the role of biology in flocculation has identified mechanisms by which this mud-sand flocculation can occur. In the first part of this paper (Manning et al., submitted), the development of an empirical flocculation model is described which represents the interaction between sand and mud particles in the flocculation process. Measurements of the settling velocity of varying mud-sand mixtures are described, and empirical algorithms governing the variation of settling velocity with turbulence, suspended sediment concentration and mud-sand content are derived. The second part of this paper continues the theme of examination of the effects of mud-sand interaction on flocculation. A 1DV mixed transport model is developed and used to reproduce the vertical transport of mixed sediment fractions. The 1DV model is used to reproduce the measured settling velocities in the laboratory experiments described in the part 1 paper and also to reproduce measurements of concentration of mixed sediments in the Outer Thames. In both modelling exercises, the model is run using the algorithms developed in part 1 and repeated using an assumption of no interaction between mud and sand in the flocculation process. The results of the modelling show a significant improvement in the ability of the 1DV to reproduce the observed sediment behaviour when the empirical equations are used. This represents further strong evidence of the interaction between sand and mud in the flocculation process.     

Multiple time scales of alluvial rivers carrying suspended sediment and their implications for mathematical modeling

Flow, sediment transport and bed deformation in alluvial rivers normally exhibit multiple time scales. Enhanced knowledge of the time scales can facilitate better approaches to the understanding of the fluvial processes. Yet prior studies of the time scales are based upon the concept of sediment transport capacity at low concentrations, which however is not generally applicable. This paper presents new formulations of the time scales of fluvial flow, suspended sediment transport and bed deformation, under the framework of shallow water hydrodynamics, non-capacity sediment transport and the theory of characteristics for the hyperbolic governing equations. The time scale of bed deformation in relation to that of flow depth is demonstrated to delimit the applicability region of mathematical river models, and the time scale of suspended sediment transport relative to that of the pertinent flow information is analyzed to address if the concept of sediment transport capacity is applicable. For shallow flows with high sediment concentrations, bed deformation may considerably affect the flow and a fully coupled model is normally required. In contrast, for deep flows at low sediment concentrations, a decoupled model is mostly justified. This pilot study of the time scales delivers a new theoretical basis, on which the interaction between flow, suspended sediment transport and bed deformation can be potentially better characterized.     

Quantification of uncertainty in design water levels due to uncertain bed form roughness in the Dutch river Waal

Hydrodynamic river models are applied to design and evaluate measures for purposes such as safety against flooding. The modelling of river processes involves numerous uncertainties, resulting in uncertain model results. Knowledge of the type and magnitude of these uncertainties is crucial for a meaningful interpretation of the model results. Uncertainty in the hydraulic roughness due to bed forms is one of the main contributors to the uncertainty in the modelled water levels. The aim of this study was to quantify the uncertainty in the bed form roughness under design conditions and quantify the effect on the design water levels in the Dutch river Waal. Five roughness models that predict bed form roughness based on measured bed form and flow characteristics were extrapolated to design conditions. The results show that the 95% confidence interval of the predicted Nikuradse roughness values under design conditions ranges from 0.32 to 1.03 m. This uncertainty was propagated through the two‐dimensional hydrodynamic model, WAQUA, by means of a Monte Carlo simulation for an idealized schematization of the Dutch river Waal. The uncertain bed form roughness results in an uncertainty in the design water levels, with a 95% confidence interval of 0.53 m, which is significant for Dutch river management practice. The uncertainty in the bed form roughness was mainly caused by a lack of knowledge about the physical process of bed form evolution that causes roughness. An improved estimation of bed form roughness can significantly reduce the uncertainty in the design water levels. Copyright © 2012 John Wiley & Sons, Ltd.     

Capturing the essential spatial variability in distributed hydrological modelling: Hydraulic roughness

The spatial variability of each parameter affecting storm runoff must be accounted for in distributed modelling. The objective of the work reported here is to assess the effects of using distributed versus lumped hydraulic roughness coefficients in the modelling of direct surface runoff. A spatially variable data set composed of Manning roughness coefficients is used to model direct surface runoff. To assess the information content (as measured by entropy) of spatially variable data and its significance in distributed modelling, various degrees of smoothing are applied. The error resulting from smoothing the hydraulic roughness coefficients is determined by modelling overland flow using a finite element solution. The Manning roughness coefficients were taken from field measurements of the Manning roughness coefficient at 0.6 m on a 14 m hillslope. These values were then used in a numerical simulation of outflow hydrographs to investigate the dependence of error on spatial variability. Our study focuses on the characteristics of spatial data used in distributed hydrological modelling. The field sites have fractal dimensions of ≈? 1.4, which is close to a Brownian variation. The sampling interval that captures the essential spatial variability of the Manning roughness coefficient does not seem to matter due to its Brownian variation in the field sites. Hence due to the nearly uniform random distribution, measurements at 0.6 m intervals are not necessary and larger intervals would yield results that are just as acceptable provided the mean value together with a uniformly random distribution is maintained for any size of finite element or sampling resolution. Because detailed measurements of hydraulic roughness are not practically available for deterministic catchment modelling, it is important to know that larger sampling resolutions may be used than 0.6 m.     

A probabilistic model for sediment entrainment:The role of bed irregularity

A generalized probabilistic model is developed in this study to predict sediment entrainment under the incipient motion, rolling, and pickup modes. A novelty of the proposed model is that it incorporates in its formulation the probability density function of the bed shear stress, instead of the near-bed velocity fluctuations, to account for the effects of both flow turbulence and bed surface irregularity on sediment entrainment. The proposed model incorporates in its formulation the collective effects of three para-meters describing bed surface irregularity, namely the relative roughness, the volumetric fraction and relative position of sediment particles within the active layer. Another key feature of the model is that it provides a criterion for estimating the lift and drag coefficients jointly based on the recognition that lift and drag forces acting on sediment particles are interdependent and vary with particle protrusion and packing density. The model was validated using laboratory data of both fine and coarse sediment and was compared with previously published models. The study results show that all the examined models perform adequately for the fine sediment data, where the sediment particles have more uniform gra-dation and relative roughness is not a factor. The proposed model was particularly suited for the coarse sediment data, where the increased bed irregularity was captured by the new parameters introduced in the model formulation. As a result, the proposed model yielded smaller prediction errors and physically acceptable values for the lift coefficient compared to the other models in case of the coarse sediment data.     

Measurements and modelling of beach groundwater flow in the swash-zone: a review

This paper reviews research on beach groundwater dynamics and identifies research questions which will need to be answered before swash zone sediment transport and beach profile evolution can be successfully modelled. Beach groundwater hydrodynamics are a result of combined forcing from the tide and waves at a range of frequencies, and a large number of observations exist which describe the shape and elevation of the beach watertable in response to tidal forcing at diurnal, semi-diurnal and spring-neap tidal frequencies. Models of beach watertable response to tidal forcing have been successfully validated; however, models of watertable response to wave forcing are less well developed and require verification. Improved predictions of swash zone sediment transport and beach profile evolution cannot be achieved unless the complex fluid and sediment interactions between the surface flow and the beach groundwater are better understood, particularly the sensitivity of sediment transport processes to flow perpendicular to the permeable bed.The presence of a capillary fringe, particularly when it lies just below the sand surface, has influences on beach groundwater dynamics. The presence of a capillary fringe can have a significant effect on the exchange of water between the ocean and the coastal aquifer, particularly in terms of the storage capacity of the aquifer. Field and laboratory observations have also shown that natural groundwater waves usually propagate faster and decay more slowly in aquifers with a capillary fringe, and observations which suggest that horizontal flows may also occur in the capillary zone have been reported. The effects of infiltration and exfiltration are generally invoked to explain why beaches with a low watertable tend to accrete and beaches with a high watertable tend to erode. However, the relative importance of processes such as infiltration losses in the swash, changes in the effective weight of the sediment, and modified shear stress due to boundary layer thinning, are not yet clear. Experimental work on the influence of seepage flows within sediment beds provides conflicting results concerning the effect on bed stability. Both modelling and experimental work indicates that the hydraulic conductivity of the beach is a critical parameter. However, hydraulic conductivity varies both spatially and temporally on beaches, particularly on gravel and mixed sand and gravel beaches. Another important, but poorly understood, consideration in beach groundwater studies is the role of air encapsulation during the wetting of beach sand.     

General velocity formula of boundary layer above mobile sediment bed induced by asymmetric waves

Sediment movement in the wave boundary layer above a mobile sediment bed is complex.A velocity formula for the boundary layer is proposed for sheet flow induced by asymmetric waves above a mobile sediment bed.The formula consists of a free stream velocity and a defect function which contains a phase-lead,boundary layer thickness and mobile sediment bed.Phase-lag of sediment movement is considered in the formula for the mobile sediment bed.The formula needs six dependent variables about asymmetric wave and sediment characteristics.Asymmetry effects on parameters(orbital amplitude,roughness height,bed shear stress,and boundary layer thickness)are properly considered such that the formula can yield velocity differences among onshore,offshore,acceleration,and deceleration stages.The formula estimates the net boundary layer velocity resulting from the mobile sediment bed and asymmetric boundary layer thickness.In addition,a non-constant phase-lead also contributes to the net boundary layer velocity in asymmetric oscillatory sheet flow.Results of the formula are as good as that of a two-phase numerical model.Sheet flow transport induced by asymmetric waves,and the offshore net sediment transport rate with a large phase-lag under velocity-skewed waves,can be adequately estimated by the formula with a power sediment concentration function.     

Patterns of sediment slug translation and dispersion following typhoon‐induced disturbance,Oyabu Creek,Kyushu, Japan

A six‐year monitoring programme characterized the migration/dispersion patterns of sediment slugs generated following typhoon‐induced disturbances in 1993 and 1997 along a single‐thread gravel‐bed stream, Oyabu Creek, on Kyushu Island, Japan. This laterally con?ned creek comprises rif?e–pool sequences with intervening bedrock outcrops. The passage of sediment pulses associated with sediment slug processes re?ected, and was controlled by, the rif?e–pool structures which provided channel bed roughness, the volume of sediment stored along valley ?oors, and the distribution of bedrock outcrops. Changes to bed material size following major sediment inputs during the disturbance events also exerted an in?uence on subsequent sediment slug processes. The sequence of rainfall events, together with changes to channel bed structure, induced different phases in the sediment slug processes. The capacity of a reach to store or trap sediment, as recorded by the longitudinal structure of the channel, varied during these differing phases. Copyright © 2004 John Wiley & Sons, Ltd.     

Modeling discharge and sediment concentrations after landscape interventions in a humid monsoon climate: The Anjeni watershed in the highlands of Ethiopia

Increasing population and intensification of agriculture increase erosion rates and often result in severe land degradation and sedimentation of reservoirs. Finding effective management practices to counteract the increasing sediment load is becoming increasingly urgent especially in the Ethiopian highlands where the construction of the hydroelectric Grand Renaissance Dam on the Blue Nile is underway. In this paper, we examine the results of 9 years of a watershed experiment in which discharge and sediment losses were observed in the 113 ha Anjeni watershed of the Blue Nile Basin. The study period encompasses conditions before, during, and after the installation of graded Fanya‐Juu (“throw uphill” bunds) soil and water conservation practices (SWCP), which had the ultimate goal of creating terraces. We use a saturation‐excess runoff model named the parameter‐efficient distributed model as a mathematical construct to relate rainfall with discharge and sediment losses at the outlet. The parameter‐efficient distributed model is based on landscape units in which the excess rainfall becomes direct runoff or infiltrates based on topographic position or hardpan characteristics. Deviations in this rainfall–discharge–sediment loss relationship are ascribed to the changes in infiltration characteristics caused by SWCPs on the hillslopes. With this technique, we found that in the Anjeni basin, the Fanya‐Juu SWCPs are only effective in increasing the infiltration and thereby reducing the direct runoff and sediment concentrations in the first 5 years. At the end of the 9‐year observation period, the direct runoff and sediment concentrations were barely reduced compared to the levels before SWCP were installed. In addition, we found that the model structure based on landscape units was able to represent the varying runoff and erosion processes during the 9 years well by varying mainly the portion of degraded land (and thereby representing the effectiveness of the Fanya‐Juu to reduce runoff by increasing infiltration).     

Analysis of surface roughness in relation to soil loss and runoff at high rainfall intensities

The decay of roughness is an important factor governing surface processes such as infiltration and soil erosion. Thus the decay of surface roughness under different surface conditions was investigated and related to quantitative amounts of soil loss, runoff and sediment concentration in a laboratory experiment. Rainfall with an intensity of 128 mm/h was applied to a bare or mulched surfaces of a sandy loam soil with known surface roughness at specified time intervals. The decay of roughness as expressed by roughness ratio, in this experiment, was better predicted when related to an exponential function of the square root of cumulative kinetic energy of rainfall rather than with the cumulative rainfall. The roughness decay equations in literature did not predict breakdown under mulched surfaces accurately. Thus the exponent parameters of the roughness decay equations were adjusted to reflect the reduced decay occurring under mulched surfaces. In a bare soil, regression equations expressing the dependent variables as a function of initial roughness index were significant, but with low coefficients of determination, being 0·39 for soil loss, 0·12 for runoff and 0·36 for sediment concentration. In addition to initial roughness index, cumulative kinetic energy of rainfall was further included in the regressions. This led to an increase in coefficients of determination, which was 0·81 for soil loss, 0·74 for runoff and 0·49 for sediment concentration. The coefficients of determination (0·87 for soil loss, 0·85 for runoff and 0·51 for sediment concentration) were further increased when the final roughness index was included in addition to initial roughness index and cumulative kinetic energy in the regressions. This work shows that soil loss and runoff could be predicted from bare soil surface provided the initial roughness and the energy of rainfall is known. However, field verifications of these relationships are needed under different tillage tools and under natural rainfall. Copyright © 2002 John Wiley & Sons, Ltd.     

The impact of submersed aquatic vegetation on the development of river mouth bars

This work is inspired by the sudden resurgence of the submersed aquatic vegetation (SAV) bed in the Chesapeake Bay (USA). Because the SAV bed occurs at the mouth of the Bay's main tributary (Susquehanna River), it plays a significant role in modulating sediment and nutrient inputs from the Susquehanna to the Bay. Previous model studies on the impact of submersed aquatic vegetation on the development of river mouth bars lacked a complete mechanistic understanding. This study takes advantage of new advances in 3D computational models that include explicit physical-sedimentological feedbacks to obtain this understanding. Specifically, we used Delft3D, a state-of-the-art hydrodynamic model that provides fine-scale computations of three-dimensional flow velocity and bed shear stress, which can be linked to sediment deposition and erosion. Vegetation is modeled using a parameterization of hydraulic roughness that depends on vegetation height, stem density, diameter, and drag coefficient. We evaluate the hydrodynamics, bed shear stresses, and sediment dynamics for different vegetation scenarios under conditions of low and high river discharge. Model runs vary the vegetation height, density, river discharge, and suspended-sediment concentration. Numerical results from the idealized model show that dense SAV on river mouth bars substantially diverts river discharge into adjacent channels and promotes sediment deposition at ridge margins, as well as upstream bar migration. Increasing vegetation height and density forms sandier bars closer to the river mouth and alteration of the bar shape. Thus, this study highlights the important role of SAV in shaping estuarine geomorphology, which is especially relevant for coastal management. © 2019 John Wiley & Sons, Ltd.     

Effects of rock fragment cover on hydrological processes under rainfall simulation in a semi‐arid region of China

Rock fragment cover has long been an important agricultural crop production technique on the Loess Plateau, China. Although this approach plays an important role in controlling hydrological processes and preventing soil erosion, inconsistent results have been recovered in this field. In this study, we investigated the effects of rock fragment cover on infiltration, run‐off, soil erosion, and hydraulic parameters using rainfall simulation in the field in a semi‐arid region of China. Two field plots encompassing 6 rock fragment coverages (0%, 10%, 20%, 25%, 30%, and 40%), as well as 2 rock fragment positions and sizes were exposed to rainfall at a particular intensity (60 mm h^?1). The results of this study showed that increasing the rock fragment coverage with rock fragments resting on the soil surface increased infiltration but decreased run‐off generation and sediment yield. A contrasting result was found, however, when rock fragments were partially embedded into the soil surface; in this case, a positive relationship between rock fragment coverage and run‐off rate as well as a nonmonotonic relationship with respect to soil loss rate was recovered. The size of rock fragments also exerted a positive effect on run‐off generation and sediment yield but had a negative effect on infiltration. At the same time, both mean flow velocity and Froude number decreased with increasing rock fragment coverage regardless of rock fragment position and size, whereas both Manning roughness and Darcy–Weisbach friction factor were positively correlated. Results show that stream power is the most sensitive hydraulic parameter affecting soil loss. Combined with variance analysis, we concluded that the order of significance of rock fragment cover variables was position followed by coverage and then size. We also quantitatively incorporated the effects of rock fragment cover on soil loss via the C and K factors in the Revised Universal Soil Loss Equation. Overall, this study will enable the development of more accurate modelling approaches and lead to a better understanding of hydrological processes under rock fragment cover conditions.