Influence of surface roughness of dune bedforms on flow and turbulence characteristics

The current paper investigates the flow and turbulence characteristics over dune bedforms by means of laboratory experiments, where spatially dense and temporally high frequency velocity measurements were done. Although similar studies are available in the literature, the focus and novelty of the current study is to assess the influence of surface roughness of the dune bedforms on the nearbed flow. For direct comparison, two different surface roughness heights over idealized, fixed-shaped, high-angled dune bedforms were tested; one with a hydraulically-smooth surface, and the other with a fully-rough surface. Spatial variation of time-averaged flow as well as turbulence statistics were examined, which was complemented by streamline plots and spectral analyses. The results are interpreted from sediment entrainment and sediment transport points of view. The results show that increased dune surface roughness reduces the nearbed flow velocity, but increases the flow velocities at upper regions. The upward directed flow near the dune crests becomes stronger in the case of smooth surface, while the re-attachment point moves further downstream compared to the rough wall case. It is concluded that the roughness of the dune surface affects the nearbed flow and turbulence characteristics qualitatively and quantitatively, which is shown to have direct consequences on sediment entrainment characteristics.     

Influence of Zostera marina canopies on unidirectional flow,hydraulic roughness and sediment movement

Seagrasses develop extensive or patchy underwater meadows in coastal areas around the world, forming complex, highly productive ecosystems. Seagrass canopies exert strong effects on water flow inside and around them, thereby affecting flow structure, sediment transport and benthic ecology. The influence of Zostera marina canopies on flow velocity, turbulence, hydraulic roughness and sediment movement was evaluated through laboratory experiments in 2 flumes and using live Z. marina and a mobile sand bed. Profiles of instantaneous velocities were measured and sediment movement was identified upstream, within and downstream of patches of different sizes and shoot density and at different free-stream velocities. Flow structure was characterised by time-averaged velocity, turbulence intensity and Turbulent Kinetic Energy (TKE). When velocity data were available above the canopy, they were fitted to the Law of the Wall and shear velocities and roughness lengths were calculated. When a seagrass canopy was present, three layers were distinguishable in the water column: (1) within canopy represented by low velocities and high turbulence; (2) transition zone around the height of the canopy, where velocities increased, turbulence decreased and TKE was high; and (3) above canopy where velocities were equal or higher than free-stream velocities and turbulence and TKE were lower than below. Shoot density and patch-width influenced this partitioning of the flow when the canopy was long enough (based on flume experiments, at least more than 1 m-long). The enhanced TKE observed at the canopy/water interface suggests that large-scale turbulence is generated at the canopy surface. These oscillations, likely to be related to the canopy undulations, are then broken down within the canopy and high-frequency turbulence takes place near the bed. This turbulence ‘cascade’ through the canopy may have an important impact on biogeochemical processes. The velocity above the canopy generally followed a logarithmic profile. Roughness lengths were higher above the canopy than over bare sand and increased with increasing distance from the leading edge of the canopy; however, they were still small (<1 cm) compared to other studies in the literature. Within and downstream of the canopy, sediment movement was observed at velocities below the threshold of motion. It was likely caused by the increased turbulence at those positions. This has large implications for sediment transport in coastal zones where seagrass beds develop.     

Observations of bedforms on a dissipative macrotidal beach

Field measurements of wave ripples and megaripples were made with a Sand Ripple Profiler in the surf and shoaling zones of a sandy macrotidal dissipative beach at Perranporth, UK in depths 1–6 m and significant wave heights up to 2.2 m. A frequency domain partitioning approach allowed quantification of height (η), length (λ) and migration rate of ripples and megaripples. Wave ripples with heights up to 2 cm and wavelengths ～20 cm developed in low orbital velocity conditions (u _m?<?0.65 m/s) with mobility number ψ?<?25. Wave ripple heights decreased with increasing orbital velocity and were flattened when mean currents were >0.1 m/s. Wave ripples were superimposed on top of megaripples (η?=?10 cm, λ?=?1 m) and contributed up to 35 % of the total bed roughness. Large megaripples with heights up to 30 cm and lengths 1–1.8 m developed when the orbital velocity was 0.5–0.8 m/s, corresponding to mobility numbers 25–50. Megaripple heights and wavelengths increased with orbital velocity but reduced when mean current strengths were >0.15 m/s. Wave ripple and megaripple migrations were generally onshore directed in the shoaling and surf zones. Onshore ripple migration rates increased with onshore-directed (+ve) incident wave skewness. The onshore migration rate reduced as offshore-directed mean flows (undertow) increased in strength and reached zero when the offshore-directed mean flow was >0.15 m/s. The migration pattern was therefore linked to cross-shore position relative to the surf zone, controlled by competition between onshore-directed velocity skewness and offshore-directed mean flow.     

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.     

Alteration of bottom roughness by benthic organisms in a sandy coastal environment

This work deals with the potential influence of benthic communities on the sediment dynamics of the coastal zone, and specifically with the modifications to bottom roughness caused by communities and their effects on wave propagation across the coastal profile. Time-series of video observations of the sea bottom on the Ebro delta coast were analysed in order to estimate the bottom roughness associated with physical and biological morphological components and the sediment reworking rates caused by epifaunal organisms. Biological roughness was mainly caused by ophiuroids and tanatocenosis of bivalve and gastropod shells, which changed their abundance during the study period. The total biological roughness (K_bio) ranged between 0.27 and 0.81 cm and represented a significant part (<20%) of the total form drag roughness. Flattening of ripples caused by bioturbation was observed under low-energy conditions. Surface sediment perturbation and bioturbation rates were also estimated. Based on these observations a wave propagation model was applied in order to carry out a sensitivity analysis of the significance of biological roughness on wave dissipation in the study area under different wave conditions.     

Application of third generation shallow water wave models in a tidal environment

Wave modeling was performed in the German Bight of the North Sea during November 2002, using the spectral wave models, namely the K-model and Simulating WAves Nearshore (SWAN), both developed for applications in environments of shallow water depths. These models mainly differ with respect to their dissipation source term expressions and in exclusion or inclusion of nonlinear wave–wave interactions. The K-model uses nonlinear dissipation and bottom dissipation, and neglects quadruplet wave–wave interaction whereas, SWAN includes, besides bottom dissipation, dissipation by white-capping and depth induced wave breaking and triad wave–wave interaction. The boundary spectra were extracted from the WAM model results of a North Sea hindcast of the HIPOCAS project, wind fields, tidal current and water level variations from the results of models used in the Belawatt project. The purpose of this study was to test the performance of both wave models to see whether they were able to predict near-shore wave conditions accurately. The runs were performed with and without tidal current and level variations to determine their effect on the waves. Comparisons of model results with buoy measurements show that taking into account tides and currents improve the spectral shape especially in areas of high current speeds. Whereas SWAN performed better in terms of spectral shape, especially in case of two peaked spectra, the K-model showed better results in terms of integrated parameters.Responsible Editor: Hans Burchard     

Impact of tidal mixing with different scales of bottom roughness on the general circulation

The current study deals with a parameterization of diapycnal diffusivity in an ocean model. The parameterization estimates the diapycnal diffusivity depending on the location of tidal-related energy dissipation over rough topography. The scheme requires a bottom roughness map that can be chosen depending on the scales of topographic features. Here, we implement the parameterization on an ocean general circulation model, and we examine the sensitivity of the modeled circulations to different spatial scales of the modeled bottom roughness. We compare three simulations that include the tidal mixing scheme using bottom roughness calculated at three different ranges of spatial scales, with the largest scale varying up to 200?km. Three main results are discussed. First, the dependence of the topographic spectra with depth, characterized by an increase in spectral energy over short length scales in the deep ocean, influences the vertical profile of the diffusivity. Second, the changes in diffusivities lead to different equilibrium solutions in the Atlantic meridional overturning circulation and bottom circulation. In particular, the lower cell of the Atlantic overturning and the bottom water transport in the Pacific Ocean are stronger for stronger diffusivities at the corresponding basins and depths, and the strongest when using the small-scale roughness map. Third, a comparison of the density fields of the three simulations with the density field of World Ocean Atlas dataset, from which the models are initialized, shows that among the simulations with three different roughness maps, the one using small-scale bottom roughness map has the smallest density bias.     

Impact of tidal mixing with different scales of bottom roughness on the general circulation

The current study deals with a parameterization of diapycnal diffusivity in an ocean model. The parameterization estimates the diapycnal diffusivity depending on the location of tidal-related energy dissipation over rough topography. The scheme requires a bottom roughness map that can be chosen depending on the scales of topographic features. Here, we implement the parameterization on an ocean general circulation model, and we examine the sensitivity of the modeled circulations to different spatial scales of the modeled bottom roughness. We compare three simulations that include the tidal mixing scheme using bottom roughness calculated at three different ranges of spatial scales, with the largest scale varying up to 200 km. Three main results are discussed. First, the dependence of the topographic spectra with depth, characterized by an increase in spectral energy over short length scales in the deep ocean, influences the vertical profile of the diffusivity. Second, the changes in diffusivities lead to different equilibrium solutions in the Atlantic meridional overturning circulation and bottom circulation. In particular, the lower cell of the Atlantic overturning and the bottom water transport in the Pacific Ocean are stronger for stronger diffusivities at the corresponding basins and depths, and the strongest when using the small-scale roughness map. Third, a comparison of the density fields of the three simulations with the density field of World Ocean Atlas dataset, from which the models are initialized, shows that among the simulations with three different roughness maps, the one using small-scale bottom roughness map has the smallest density bias.

     

Three‐dimensional flow structure around small‐scale bedforms in a simulated gravel‐bed environment

Pebble clusters are common small‐scale morphological features in gravel‐bed rivers, occupying as much as 10 per cent of the bed surface. Important links exist between the presence of pebble clusters and the development of flow structures. These links are poorly understood at the three‐dimensional level. Particularly neglected has been the effect of clusters on the lateral flow characteristics. A laboratory study was conducted using a hydraulic flume, within which simulated pebble clusters were superimposed onto a plane bed of gravel material. High‐resolution three‐dimensional flow data were collected above the bed at two different flow depths using an acoustic Doppler velocimeter. The results present evidence of the importance of lateral flow in the development of turbulent flow structure. Narrow regions of high lateral and downstream turbulence intensity exist to both sides of clusters and in a three‐dimensional separation zone in their lee. This may indicate the presence of horseshoe‐type vortical structures analogous to those identified in less hydraulically rough environments. However, it is likely that these structures are more complicated given the mutual interference of the surrounding medium. The lateral flow was also identified as a key component in the upwelling identified by other authors in the lee of pebble clusters. The results of the vertical flow analysis confirm the hypothesis that six regions with distinct vertical flow characteristics exist above clusters: flow acceleration up the stoss‐side of the cluster; recirculation behind the cluster in the wake region; vortex shedding from the pebble crest and shear layer; flow reattachment downstream of the cluster; upwelling of flow downstream of the point of reattachment; and recovery of flow. Copyright © 2001 John Wiley & Sons, Ltd.     

Lateral sediment transport by bedforms in a meander bend

Bedforms and their migration are among the most important features determining bed load transport on point bars in sand-bedded rivers. A case study has been carried out in the Danish river Gels Å in a bend with equilibrium point bar topography. Bedforms (dunes) were repeatedly surveyed to determine the bedform migration on the point bar surface. Bedforms occur in two categories in the bend: primary dunes and secondary dunes migrating transversely in the troughs between primary dunes. The study shows that the migration of these bedforms can account for the entire lateral transport known to concentrate the transport in the convex half of the cross section at the bend exit. The bedform transport is induced by both dune categories, partly by turning of the primary dunes towards the convex bank and partly by the migration of secondary dunes. In certain areas, the two categories operate together and thus make the lateral transport particularly large; it is shown that almost 60 per cent of the entire lateral transport in the bend takes place within less than 30 per cent of the streamwise extent of the point bar.     

Transport of nonsorbing solutes in a streambed with periodic bedforms

Previous studies of hyporheic zone focused largely on the net mass transfer of solutes between stream and streambed. Solute transport within the bed has attracted less attention. In this study, we combined flume experiments and numerical simulations to examine solute transport processes in a streambed with periodic bedforms. Solute originating from the stream was subjected to advective transport driven by pore water circulation due to current–bedform interactions as well as hydrodynamic dispersion in the porous bed. The experimental and numerical results showed that advection played a dominant role at the early stage of solute transport, which took place in the hyporheic zone. Downward solute transfer to the deep ambient flow zone was controlled by transverse dispersion at the later stage when the elapsed time exceeded the advective transport characteristic time t_c (= L/u_c with L being the bedform length and u_c the characteristic pore water velocity). The advection-based pumping exchange model was found to predict reasonably well solute transfer between the overlying water and streambed at the early stage but its performance deteriorated at the later stage. With dispersion neglected, the pumping exchange model underestimated the long-term rate and total mass of solute transfer from the overlying water to the bed. Therefore both advective and dispersive transport components are essential for quantification of hyporheic exchange processes.     

Inverse modeling of coastal aquifers using tidal response and hydraulic tests

Remediation of contaminated aquifers demands a reliable characterization of hydraulic connectivity patterns. Hydraulic diffusivity is possibly the best indicator of connectivity. It can be derived using the tidal response method (TRM), which is based on fitting observations to a closed-form solution. Unfortunately, the conventional TRM assumes homogeneity. The objective of this study was to overcome this limitation and use tidal response to identify preferential flowpaths. Additionally, the procedure requires joint inversion with hydraulic test data. These provide further information on connectivity and are needed to resolve diffusivity into transmissivity and storage coefficient. Spatial variability is characterized using the regularized pilot points method. Actual application may be complicated by the need to filter tidal effects from the response to pumping and by the need to deal with different types of data, which we have addressed using maximum likelihood methods. Application to a contaminated artificial coastal fill leads to flowpaths that are consistent with the materials used during construction and to solute transport predictions that compare well with observations. We conclude that tidal response can be used to identify connectivity patterns. As such, it should be useful when designing measures to control sea water intrusion.     

Influence of geometrical variations on morphodynamic equilibria in short tidal basins

Ocean Dynamics - The existence of cross-sectionally averaged morphodynamic equilibria of tidal inlets is investigated, using a cross-sectionally averaged model, and their sensitivity to variations...     

Flow resistance and hydraulic geometry in bedrock rivers with multiple roughness length scales

Many models of incision by bedrock rivers predict water depth and shear stress from discharge; conversely, palaeoflood discharge is sometimes reconstructed from flow depth markers in rock gorges. In both cases, assumptions are made about flow resistance. The depth–discharge relation in a bedrock river must depend on at least two roughness length scales (exposed rock and sediment cover) and possibly a third (sidewalls). A conceptually attractive way to model the depth–discharge relation in such situations is to partition the total shear stress and friction factor, but it is not obvious how to quantify the friction factor for rough walls in a way that can be used in incision process models. We show that a single flow resistance calculation using a spatially averaged roughness length scale closely approximates the partitioning of stress between sediment and rock, and between bed and walls, in idealized scenarios. Both approaches give closer fits to the measured depth–discharge relations in two small bedrock reaches than can be achieved using a fixed value of Manning's n or the Chézy friction factor. Sidewalls that are substantially rougher or smoother than the bed have a significant effect on the partitioning of shear stress between bed and sidewalls. More research is needed on how best to estimate roughness length scales from observable or measurable channel characteristics. © 2019 John Wiley & Sons, Ltd.     

Effect of seepage on flow and bedforms dynamics

This study, using an experimental approach, focuses on the effect of downward seepage on a threshold alluvial channel morphology and corresponding turbulent flow characteristics. In all the experiments, we observed that the streamwise time‐averaged velocities and Reynolds shear stresses were increased under the influence of downward seepage. Scales of eddy length and eddy turnover time were significantly increased with the application of downward seepage, leading to sediment transport and initiation of bedforms along the channel length. As the amount of seepage discharge increased, eddy length and turnover time were further increased, causing the development of larger bedforms. It was revealed that the geometry of bedforms was linked with the size of eddies. In this work, statistics of bedform dynamics are presented in terms of multi‐scalar bedforms in the presence of seepage. These multi‐scalar ubiquitous bedforms cast a potential impact on flow turbulence as well as stream bed morphology in channels. We used wavelet to analyse temporally lagged spatial bed elevation profiles that were obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuations at different length scales. A spatial cross‐correlation analysis, based on the wavelet coefficients, was performed on these bed elevation datasets to observe the effect of downward seepage on the dynamic behaviour of bedforms at different length scales. It was found that celerity of bedforms reduced with increase in seepage percentage. Bedform celerity was best approximated by a probability density function such as Rayleigh distribution under varying downward seepage. Further, statistical analysis of physical parameters of bedforms ascertained that the reduction in bedform celerity was a result of increased bedform size. Copyright © 2017 John Wiley & Sons, Ltd.     

Influence of variable salinity conditions in a tidal creek on riparian groundwater flow and salinity dynamics

While recent studies have revealed that tidal fluctuations in an estuary significantly affect groundwater flows and salt transport in the riparian zone, only seawater salinity in the estuary has been considered. A numerical study is conducted to investigate the influence of estuarine salinity variations on the groundwater flow and salt dynamics in the adjacent aquifer to extend our understanding of these complex and dynamic systems. Tidal salinity fluctuations (synchronous with estuary stage) were found to alter the magnitude and distribution of groundwater discharge to the estuary, which subsequently impacted on groundwater salinity patterns and residence times, especially in the riparian zone. The effects of salinity fluctuations were not fully captured by adopting a constant, time-averaged estuarine salinity. The modelling analysis also included an assessment of the impact of a seasonal freshwater flush in the estuary, similar to that expected in tropical climates (e.g. mean estuary level during flood significantly greater than average), on adjacent groundwater flow and salinity conditions. The three-month freshwater flushing event temporarily disrupted the salt distribution and re-circulation patterns predicted to occur under conditions of constant salinity and tidal water level fluctuations in the estuary. The results indicate that the salinity variations in tidal estuaries impact significantly on estuary–aquifer interaction and need to be accounted for to properly assess salinity and flow dynamics and groundwater residence times of riparian zones.     

Influence of varying tidal prism on hydrodynamics and sedimentary processes in a hypertidal salt marsh creek

Recent initiatives directing tidal power development in the Bay of Fundy have raised questions about far‐field environmental impacts related to energy extraction. It is understood that commercial scale tidal power installations in the Minas Passage will result in an overall decrease in tidal amplitude in the Minas Basin. Corresponding changes in sedimentation patterns may or may not be within the natural range of variability, and it is hypothesized that intertidal sedimentation rates will demonstrate a non‐linear response to modification of the tidal energy regime. This research considers acoustic Doppler velocimeter (ADV) and optical backscatter sensor (OBS) data from a sheltered tidal creek in the Minas Basin, for analysis of tidal characteristics in a hypertidal creek environment over spring and neap tidal cycles. Sediment deposition in the creek was also measured. Results show a first‐order control of topography on flow magnitude in the tidal creek, which impacts net sediment deposition through resuspension and removal of newly introduced material. This study demonstrates that tides which peak around the bankfull level show reduced early ebb stage turbulence and flow velocity and encourage an extended depositional period. The dynamics of marshfull tides may be responsible for the maximum sediment deposition in tidal creeks, providing large amounts of material that is eventually distributed to and deposited on marsh surfaces. Copyright © 2012 John Wiley & Sons, Ltd.     

Influence of bedding orientation on the anisotropy of hydraulic conductivity in a well-sorted fluvial sediment

The paper describes a permeameter test method for determination of the hydraulic conductivity(AT) along multi-directions in fluvial sediments with cross beddings.Unlike existing in-situ permeameter methods that determine hydraulic conductivity for submerged streambeds,our method was intended to measure hydraulic conductivity of exposed streambeds or fluvial sediments.The method was applied to the Wei River,Shaanxi Province,Central China for characterization of the anisotropy of K in a well-sorted fluvial sediment.The results illustrated that even in well-sorted sediments,cross-bedding and sediment fabrication(or texture) can lead to varied K values along different measurement directions.The K value was the largest along the dip direction(or the major direction) that is parallel to the orientation of cross bedding and the smallest in the direction perpendicular to the bedding(or the minor direction). The K value in a given direction between the major and minor direction often fell in the range bounded by the K values in the major and minor directions.The anisotropy ratio of K(the ratio of K value between the major and minor directions) in two trenches for this well-sorted fluvial sediment was up to 1.14 to 1.23,respectively.Our results also demonstrated that even for well-sorted sediments,the K values between two sampling points only about 10 cm apart can differ.It is clear that the K distribution strongly correlates to the bedding orientation.