2D numerical analysis of the influence of near-bank vegetation patches on the bed morphological adjustment

This paper investigates the influence of near-bank vegetation patches on the bed morphological adjustment in open channel flow systems. The 2D depth-averaged hydro-morphological model is adopted for this investigation, which is first validated by laboratory experimental data measured in an open channel with a single near-bank vegetation patch. The validated model is then applied for extensive numerical simulations, with the aim of conducting a systematic analysis of the influence of different geometric controlling parameters on the bed morphological evolution. The controlling parameters taken into account for numerical analysis include the angle of repose value (RAV) of sediment, vegetation density (VD), patch length (PL) and patch width (PW). The numerical results and analysis show that: (1) the RAV of sediment with slope-failure parametrization only influences the shape of the transverse bed topography in the junction region; (2) increase in VD, PL and PW that substantially enhances flow blockage effect encourages the growth of the pool adjacent to the patch in three dimensions; (3) increase in VD, PL and PW produces analogous retrogressive erosion (erosion toward the upstream) in the pool region, presumably due to the increase in flow resistance. Additional numerical experiments suggest that the staggered-order distribution of multiple patches might be an optimal choice for channel restoration and conservation since pools and riffles with larger scales can be produced.

     

Hydrodynamics and turbulence in emergent and sparsely vegetated open channel flow

This present study reports the results of an experimental study characterizing thorough variation of turbulent hydrodynamics and flow distribution in emergent and sparsely vegetated open channel flow. An emergent and rigid sparse vegetation patch with regular spacing between stems along the flow and transverse directions was fixed in the central region of the cross-section of open channel. Experiments were conducted in subcritical flow conditions and velocity measurements were obtained with an acoustic Doppler Velocimetry system. Large variations of the turbulence intensities, Reynolds shear stress, turbulent kinetic energy and vortical motions are found in and around the vegetation patch. At any cross-section through the interior of the vegetation patch, streamwise velocity decreases with increase in streamwise length and the velocity profiles converge from the log-law to a linear profile with increasing slope. Time-averaged lateral and vertical velocities inside the vegetation patch increase with increasing streamwise distance and converge from negative values to positive values. Turbulence intensities interior of the sparse vegetation patch are more than those of without the vegetation patch. Similar to the trend of streamwise velocity profiles inside the vegetation, turbulence intensities and longitudinal-normal Reynolds shear stress profile decreases with streamwise direction. In the interior of the vegetation patch and downstream of the trailing edge, turbulent kinetic energy profiles are exhibiting irregular fluctuations and the maximum values are occurring in the outer layer. Analysis of flow distribution confirms sparse vegetation patch is inducing a serpentine flow pattern in its vicinity. At the leading edge, flow is rushing towards the right hand sidewall, and at the trailing edge, flow is turning to the left hand sidewall. In between the leading and trailing edges, the streamlines are following a zig-zag fashion at varied degree along the streamwise and lateral directions. Immediate upstream of the leading edge and in the interior of the vegetation patch, vortex motion is clearly visible and the vortices are stretched along the width of the channel with streamwise direction.     

Characteristics of flow structure of free-surface flow in a partly obstructed open channel with vegetation patch

Free-surface flows over patchy vegetation are common in aquatic environments. In this study, the hydrodynamics of free-surface flow in a rectangular channel with a bed of rigid vegetation-like cylinders occupying half of the channel bed was investigated and interpreted by means of laboratory experiments and numerical simulations. The channel configurations have low width-to-depth aspect ratio (1.235 and 2.153). Experimental results show that the adjustment length for the flow to be fully developed through the vegetation patch in the present study is shorter than observed for large-aspect-ratio channels in other studies. Outside the lateral edge of the vegetation patch, negative velocity gradient (\(\partial \overline{u}/\partial z < 0\)) and a local velocity maximum are observed in the vertical profile of the longitudinal velocity in the near-bed region, corresponding to the negative Reynolds stress (\(- \overline{{u^{\prime}w^{\prime}}} < 0\)) at the same location. Assuming coherent vortices to be the dominant factor influencing the mean flow field, an improved Spalart–Allmaras turbulence model is developed. The model improvement is based on an enhanced turbulence length scale accounting for coherent vortices due to the effect of the porous vegetation canopy and channel bed. This particular flow characteristic is more profound in the case of high vegetation density due to the stronger momentum exchange of horizontal coherent vortices. Numerical simulations confirmed the local maximum velocity and negative gradient in the velocity profile due to the presence of vegetation and bed friction. This in turn supports the physical interpretation of the flow processes in the partly obstructed channel with vegetation patch. In addition, the vertical profile of the longitudinal velocity can also be explained by the vertical behavior of the horizontal coherent vortices based on a theoretical argument.     

Spatial evolution of coherent motions in finite-length vegetation patch flow

A number of experimental studies on submerged canopy flows have focused on fully-developed flow and turbulent characteristics. In many natural rivers, however, aquatic vegetation occurs in patches of finite length. In such vegetated flows, the shear layer is not formed at the upstream edge of the vegetation patch and coherent motions develop downstream. Therefore, more work is neededz to reveal the development process for large-scale coherent structures within vegetation patches. For this work, we considered the effect of a limited length vegetation patch. Turbulence measurements were intensively conducted in open-channel flows with submerged vegetation using Particle Image Velocimetry (PIV). To examine the transition from boundary-layer flow upstream of the vegetation patch to a mixing-layer-type flow within the patch, velocity profiles were measured at 33 positions in a longitudinal direction. A phenomenological model for the development process in the vegetation flow was developed. The model decomposed the entire flow region into four zones. The four zones are the following: (i) the smooth bed zone, (ii) the diverging flow zone, (iii) the developing zone and (iv) the fully-developed zone. The PIV data also confirmed the efficiency of the mixing-layer analogy and provided insight into the spatial evolution of coherent motions.     

Measurements and modeling of open-channel flows with finite semi-rigid vegetation patches

The hydrodynamics of flows through a finite length semi-rigid vegetation patch (VP) were investigated experimentally and numerically. Detailed measurements have been carried out to determine the spatial variation of velocity and turbulence profiles within the VP. The measurement results show that an intrusion region exists in which the peak Reynolds stress remains near the bed. The velocity profile is invariant within the downstream part of the VP while the Reynolds stress profile requires a longer distance to attain the spatially invariant state. Higher vegetation density leads to a shorter adjustment length of the transition region, and a higher turbulence level within the VP. The vegetation density used in the present study permits the passing through of water and causes the peak Reynolds stress and turbulence kinetic energy each the maximum at the downstream end of the patch. A 3D Reynolds-averaged Navier–Stokes model incorporating the Spalart–Allmaras turbulence closure was employed subsequently to replicate the flow development within the VP. The model reproduced transitional flow characteristics well and the results are in good agreement with the experimental data. Additional numerical experiments show that the adjustment length can be scaled by the water depth, mean velocity and maximum shear stress. Empirical equations of the adjustment lengths for mean velocity and Reynolds stress were derived with coefficients quantified from the numerical simulation results.     

Compound channel flow with a longitudinal transition in hydraulic roughness over the floodplains

Flows in a compound open-channel (two-stage geometry with a main channel and adjacent floodplains) with a longitudinal transition in roughness over the floodplains are experimentally investigated in an 18 m long and 3 m wide flume. Transitions from submerged dense vegetation (meadow) to emergent rigid vegetation (wood) and vice versa are modelled using plastic grass and vertical wooden cylinders. For a given roughness transition, the upstream discharge distribution between main channel and floodplain (called subsections) is also varied, keeping the total flow rate constant. The flows with a roughness transition are compared to flows with a uniformly distributed roughness over the whole length of the flume. Besides the influence of the downstream boundary condition, the longitudinal profiles of water depth are controlled by the upstream discharge distribution. The latter also strongly influences the magnitude of the lateral net mass exchanges between subsections, especially upstream from the roughness transition. Irrespective of flow conditions, the inflection point in the mean velocity profile across the mixing layer is always observed at the interface between subsections. The longitudinal velocity at the main channel/floodplain interface, denoted \(U_{int}\), appeared to be a key parameter for characterising the flows. First, the mean velocity profiles across the mixing layer, normalised using \(U_{int}\), are superimposed irrespective of downstream position, flow depth, floodplain roughness type and lateral mass transfers. However, the profiles of turbulence quantities do not coincide, indicating that the flows are not fully self-similar and that the eddy viscosity assumption is not valid in this case. Second, the depth-averaged turbulent intensities and Reynolds stresses, when scaled by the depth-averaged velocity \(U_{d,int}\) exhibit two plateau values, each related to a roughness type, meadow or wood. Lastly, the same results hold when scaling by \(U_{d,int}\) the depth-averaged lateral flux of momentum due to secondary currents. Turbulence production and magnitude of secondary currents are increased by the presence of emergent rigid elements over the floodplains. The autocorrelation functions show that the length of the coherent structures scales with the mixing layer width for all flow cases. It is suggested that coherent structures tend to a state where the magnitude of velocity fluctuations (of both horizontal vortices and secondary currents) and the spatial extension of the structures are in equilibrium.     

Analytical model for vertical velocity profiles in flows with submerged shrub-like vegetation

An analytical solution for the vertical profiles of the horizontal velocity of channel flow with submerged shrub-like vegetation is investigated in this paper. At first, a shape function is proposed to fit the diameter change of different types of shrub-like vegetation. Using the momentum theorem and the mixing-length turbulence model, an analytical solution for the vertical profile of the horizontal velocity within the vegetation is obtained. The velocity distribution of the whole column is determined in tandem with the logarithmic velocity profile above the vegetation. The solution is compared with experimental data in excellent agreement. The results show that the flow above the vegetation has a logarithmic velocity profile while the flow within the vegetation is impacted greatly by the shape and density of vegetation. The flows within shrub-like vegetations are non-uniform and vary inversely with the shrub diameter.     

Double-average mean flow and local turbulence intensity profiles from PIV measurements for an open channel flow with rigid vegetation

This contribution presents particle image velocimetry measurements for an open channel stationary uniform and fully developed flow of water over a horizontal flat bed of uniform glass beads in presence of a staggered array of vertical cylindrical stems. The main objective was to explore and quantify the influence of the stems-to-flow relative submergence, h _v /h, over the mean flow and local turbulence intensities. A comparison with measurements for the non-vegetated flow over the same granular bed is presented. Results indicate a remarkable influence of h _v /h over the whole flow field. The time-average mean flow presents a strong spatial variation in the layer of the flow occupied by the stems. The local velocity fluctuations are strongly affected by the presence of the stems, with regions in between the stems where they reach peaks that are several times larger than those encountered in the flow in absence of vegetation. The turbulence intensity profiles are noticeably different when compared to those measured in the non-vegetated flow conditions. From previous works it was possible to derive an equation for the mean velocity, U _v , of the flow through the vegetated layer of height h _v . The prediction of this equation is in good agreement with the uniform value for the double-average longitudinal velocity profile in this layer. A final brief discussion about the possible impact of these vegetated-flow features on the sediment transport is presented.     

Experimental investigation of bubbly flow and turbulence in hydraulic jumps

Many environmental problems are linked to multiphase flows encompassing ecological issues, chemical processes and mixing or diffusion, with applications in different engineering fields. The transition from a supercritical flow to a subcritical motion constitutes a hydraulic jump. This flow regime is characterised by strong interactions between turbulence, free surface and air–water mixing. Although a hydraulic jump contributes to some dissipation of the flow kinetic energy, it is also associated with increases of turbulent shear stresses and the development of turbulent eddies with implications in terms of scour, erosion and sediment transport. Despite a number of experimental, theoretical and numerical studies, there is a lack of knowledge concerning the physical mechanisms involved in the diffusion and air–water mixing processes within hydraulic jumps, as well as on the interaction between the free-surface and turbulence. New experimental investigations were undertaken in hydraulic jumps with Froude numbers up to Fr = 8.3. Two-phase flow measurements were performed with phase-detection conductivity probes. Basic results related to the distributions of void fraction, bubble frequency and mean bubble chord length are presented. New developments are discussed for the interfacial bubble velocities and their fluctuations, characterizing the turbulence level and integral time scales of turbulence representing a “lifetime” of the longitudinal bubbly flow structures. The analyses show good agreement with previous studies in terms of the vertical profiles of void fraction, bubble frequency and mean bubble chord length. The dimensionless distributions of interfacial velocities compared favourably with wall-jet equations. Measurements showed high turbulence levels. Turbulence time scales were found to be dependent on the distance downstream of the toe as well as on the distance to the bottom showing the importance of the lower (channel bed) and upper (free surface) boundary conditions on the turbulence structure.     

Vegetation cover influence on the morphology and migration patterns of river mouths

The present study focuses on the effects of vegetation cover changes on the dynamic morphology of seven southeastern Mediterranean river mouths. The methodology used comprised monitoring and mapping by GIS techniques, with data derived from historic aerial photographs, which were applied in the investigation of the morphological spatial and temporal migration patterns of the mouths, and subsequent analysis of the vegetation cover changes influencing them. Vegetation cover adjacent to river mouths influences river mouth morphology through five primary mechanisms: a) bank vegetation; b) dune advancement toward the shoreline; c) changes in the beach??s micro-topography; as well as d) long-term continuous channel migration through permanent vegetation patches; and e) channel switching through permanent vegetation patches. The five mechanisms are part of a system of interactions between channel water flow and fluvial processes; coastal sediment transport and coastal processes; and the evolution of plant communities. In the interplay between these factors they all affect and are being affected by one another. In many river mouths artificial channel diversion is often needed due to uncontrolled channel migration. It is demonstrated that vegetation cover can serve as a mean of ??soft?? channel regulation. Therefore, a better understanding of the five influencing mechanisms may aid in controlling and managing river mouth migration patterns. The study contributes to the knowledge about bank vegetation as a tool of ??soft?? channel regulation and thus can contribute to the improvement of coastal zone management.     

The effects of Pinus tabuliformis on soil detachment under different influencing factors in the Loess Plateau of China

In situ experiments were undertaken to examine the impacts of Pinus tabuliformis on soil detachment under different influencing factors. Experimental sites with slopes containing vegetation cover, plant roots and bare ground were investigated for slope gradients of 8.7%, 17.6%, 26.8%, 36.4% and 46.6% and flow discharges of 0.25, 0.5, 1.0, 1.5 and 2.0?m³?s^?1. Results from our study provide details on the relationship between soil detachment and both slope gradient and flow discharge. Soil detachment rates increased with an increase in slope and discharge, and discharge was identified to have a stronger influence on soil detachment than slope. In contrast, the combination of litter cover and roots played an important role in reducing soil detachment; detachment rate decreased by about 55% relative to bare slopes; and plant roots had a greater impact on detachment reduction. The presence of vegetation cover can decrease rill erodibility by 83% and increase critical shear stress by 224% compared with bare slopes. These results provide valuable information on the importance of woodland in soil detachment control, and may help to improve vegetation construction in seriously eroded regions of the Loess Plateau.     

Large eddy simulation of flow and scalar transport in a vegetated channel

Predicting flow and mass transport in vegetated regions has a broad range of applications in ecology and engineering practice. This paper presents large eddy simulation (LES) of turbulent flow and scalar transport within a fully developed open-channel with submerged vegetation. To properly represent the scalar transport, an additional diffusivity was introduced within the canopy to account for the contribution of stem wakes, which were not resolved by the LES, to turbulent diffusion. The LES produced good agreement with the velocity and concentration fields measured in a flume experiment. The simulation revealed a secondary flow distributed symmetrically about the channel centerline, which differed significantly from the circulation in a bare channel. The secondary circulation accelerated the vertical spread of the plume both within and above the canopy layer. Quadrant analysis was used to identify the form and shape of canopy-scale turbulent structures within and above the vegetation canopy. Within the canopy, sweep events contributed more to momentum transfer than ejection events, whereas the opposite occurred above the canopy. The coherent structures were similar to those observed in terrestrial canopies, but smaller in scale due to the constraint of the water surface.     

A new view of dispersion in well-mixed estuaries

A two-dimensional simulation of Delaware estuary hydrodynamics has been constructed. This simulation has been achieved through a rational estimate of the character of natural turbulence. Non-homogeneous velocities, on the cross-section, are employed in two-dimensional, laterally homogeneous species mass balances. In turn, concentration profiles are interpreted in the form of classical, one-dimensional dispersion coefficients. Variation of dispersion as a function of both freshwater inflow and longitudinal distance was generated. Variation of dispersion in time within a tidal cycle was found to be insignificant while no significant variation from one tidal cycle to the next has been detected.The modeling process involves the solution of tractable equations by implicit numerical methods and is capable of being excited by a wide range of input conditions.A study of the sensitivity of dispersion due to vertical mass diffusion revealed that longitudinal mixing characteristics are inversely proportional to vertical eddy diffusivity and analysis of the numerical results showed the dispersion coefficient is essentially insensitive to variation of longitudinal mass diffusivity. This leads to the conclusion that turbulent diffusivity of mass in the longitudinal direction may be taken as constant for most purposes in the study of a two-dimensional species mass balance model.A field program was carried out near the Delaware Memorial Bridge to collect velocity profiles. Substantial portions of the scheme have been verified (i.e. one- and two-dimensional tidal dynamic models) through the use of these data.     

Measuring Protected‐Area Isolation and Correlations of Isolation with Land‐Use Intensity and Protection Status

Abstract: Protected areas cover over 12% of the terrestrial surface of Earth, and yet many fail to protect species and ecological processes as originally envisioned. Results of recent studies suggest that a critical reason for this failure is an increasing contrast between the protected lands and the surrounding matrix of often highly altered land cover. We measured the isolation of 114 protected areas distributed worldwide by comparing vegetation‐cover heterogeneity inside protected areas with heterogeneity outside the protected areas. We quantified heterogeneity as the contagion of greenness on the basis of NDVI (normalized difference vegetation index) values, for which a higher value of contagion indicates less heterogeneous land cover. We then measured isolation as the difference between mean contagion inside the protected area and mean contagion in 3 buffer areas of increasing distance from the protected‐area border. The isolation of protected areas was significantly positive in 110 of the 114 areas, indicating that vegetation cover was consistently more heterogeneous 10–20 km outside protected areas than inside their borders. Unlike previous researchers, we found that protected areas in which low levels of human activity are allowed were more isolated than areas in which high levels are allowed. Our method is a novel way to assess the isolation of protected areas in different environmental contexts and regions.     

Hydrodynamics and bed stability of open channel flows with submerged foliaged plants

In this work we investigate experimentally and numerically the flow structure around foliaged plants deployed in a channel with gravels on the bed under submerged and barely submerged conditions. Velocity and Reynolds stress were measured by using a NORTEK Vectrino profiler. Visual observation shows that the initial motion of gravels is easier to be triggered under the condition of flow with barely submerged vegetation. This is confirmed by the measured velocity, Reynolds stress and total kinetic energy (TKE) profiles. The velocity exhibits a speed up in the near-bed region, and the associated Reynolds stress and TKE increase there. A 3D numerical model is then verified against the experiments and used to investigate systematically the effect of degree of submergence of foliaged plants on the channel bed shear stress. The results show that the maximum bed shear stress occurs when the foliage is situated slightly below the water surface, which can enhance channel bed instability.     

Energy dissipation, flow resistance and gas-liquid interfacial area in skimming flows on moderate-slope stepped spillways

With the re-evaluation and revision of a number of design floods, several embankment overtopping protection systems have been developed and a common technique is the construction of a stepped spillway on the downstream slope. For such moderate slope stepped channels, detailed air–water flow measurements were performed in a large facility with a focus on the rate of energy dissipation, flow resistance, air–water interfacial areas and re-aeration rates. Past and present experimental results showed a significant aeration of the flow. The median dimensionless residual head was about 3 × d_c for the 21.8° sloping chute and smaller than that for flatter slopes (θ = 3.4° and 15.9°). The flow resistance results yielded an equivalent Darcy friction factor of about 0.25 implying a larger flow resistance for the 21.8° slope angle than for smaller slope angles. The re-aeration rate was deduced from the integration of the mass transfer equation using measured air–water interfacial areas and air–water flow velocities. The results suggested an increasing re-aeration rate with increasing rate of energy dissipation. The stepped invert contributed to intense turbulence production, free-surface aeration and large interfacial areas. The experimental data showed however some distinctive seesaw pattern in the longitudinal distribution of air–water flow properties with a wave length of about two step cavities. While these may be caused by the interactions between successive adjacent step cavities and their interference with the free-surface, the existence of such “instabilities” implies that the traditional concept of normal flow might not exist in skimming flows above moderate-slope stepped spillways.     

Turbulence and waves inside flexible-wall systems designed for biological studies

Simultaneous flow measurements were performed by hot-wire anemometry inside and outside (a) a closed bag, and (b) a flow-through system designed for metabolic experiments of Fucus vesiculosus communities. Since water movement is considered an important parameter in such biological studies, the walls of the systems were flexible in order to establish flow conditions within the enclosed water body comparable to those in the natural environment. In situ experiments in the Baltic Sea at a water depth of 2 m showed that energy spectra inside the systems were comparable to those outside for a variety of flow and wave conditions. Thus, biological data from glexible wall systems, carefully designed to meet specific natural flow requirements, can be taken as reliable input data for natural ecosystem modeling.Contribution No. 147 of the Joint Research Programme 95 Interaction Sea-Sea Bottom, Kiel University, Kiel (FRG), and Hawaii Institute of Geophysics Contribution No. 803.     

Benefits of protected areas for nonbreeding waterbirds adjusting their distributions under climate warming

Climate warming is driving changes in species distributions and community composition. Many species have a so-called climatic debt, that is, shifts in range lag behind shifts in temperature isoclines. Inside protected areas (PAs), community changes in response to climate warming can be facilitated by greater colonization rates by warm-dwelling species, but also mitigated by lowering extirpation rates of cold-dwelling species. An evaluation of the relative importance of colonization-extirpation processes is important to inform conservation strategies that aim for both climate debt reduction and species conservation. We assessed the colonization-extirpation dynamics involved in community changes in response to climate inside and outside PAs. To do so, we used 25 years of occurrence data of nonbreeding waterbirds in the western Palearctic (97 species, 7071 sites, 39 countries, 1993–2017). We used a community temperature index (CTI) framework based on species thermal affinities to investigate species turnover induced by temperature increase. We determined whether thermal community adjustment was associated with colonization by warm-dwelling species or extirpation of cold-dwelling species by modeling change in standard deviation of the CTI (CTI_SD). Using linear mixed-effects models, we investigated whether communities in PAs had lower climatic debt and different patterns of community change than communities outside PAs. For CTI and CTI_SD combined, communities inside PAs had more species, higher colonization, lower extirpation, and lower climatic debt (16%) than communities outside PAs. Thus, our results suggest that PAs facilitate 2 independent processes that shape community dynamics and maintain biodiversity. The community adjustment was, however, not sufficiently fast to keep pace with the large temperature increases in the central and northeastern western Palearctic. Our results underline the potential of combining CTI and CTI_SD metrics to improve understanding of the colonization-extirpation patterns driven by climate warming.     

Comparison of Two Types of Metapopulation Models in Real and Artificial Landscapes

Abstract: Application of metapopulation models is becoming increasingly widespread in the conservation of species in fragmented landscapes. We provide one of the first detailed comparisons of two of the most common modeling techniques, incidence function models and stage-based matrix models, and test their accuracy in predicting patch occupancy for a real metapopulation. We measured patch occupancies and demographic rates for regional populations of the Florida scrub lizard (   Sceloporus woodi ) and compared the observed occupancies with those predicted by each model. Both modeling strategies predicted patch occupancies with good accuracy ( 77–80%) and gave similar results when we compared hypothetical management scenarios involving removal of key habitat patches and degradation of habitat quality. To compare the two modeling approaches over a broader set of conditions, we simulated metapopulation dynamics for 150 artificial landscapes composed of equal-sized patches (2–1024 ha) spaced at equal distances (50–750 m). Differences in predicted patch occupancy were small to moderate (<20%) for about 74% of all simulations, but 22% of the landscapes had differences openface> 50%. Incidence function models and stage-based matrix models differ in their approaches, assumptions, and requirements for empirical data, and our findings provide evidence that the two models can produce different results. We encourage researchers to use both techniques and further examine potential differences in model output. The feasibility of obtaining data for population modeling varies widely among species and limits the modeling approaches appropriate for each species. Understanding different modeling approaches will become increasingly important as conservation programs undertake the challenge of managing for multiple species in a landscape context.