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.     

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.

     

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.     

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.     

Modeling the longitudinal profiles of streamwise velocity in an open channel with a model patch of vegetation

This paper proposes a model for predicting the longitudinal profiles of streamwise velocities in an open channel with a model patch of vegetation. The governing equation was derived from the momentum equation and flow continuity equation. The model can estimate the longitudinal profiles of velocities both inside and outside a vegetation patch. Laboratory experiments indicate that the longitudinal profiles of velocities inside a patch and in the adjacent bare channel have the same adjustment distance in the longitudinal direction, but the profiles have different trends because the vegetation drag drives the flow from the patch to the adjacent bare channel. The model considers different dimensionless parameters in two flow adjustment regions upstream of and inside the patch. Sixteen sets of experimental data from different sources are used to verify the model. The model is capable of modeling the longitudinal profiles of velocities inside and outside patches of cylinders or cylinder-like plants. Compared to a previous model, the current model improves the modeling accuracy of longitudinal profiles of velocities.

     

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.     

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.     

Flow regimes in gaps within stands of flexible vegetation: laboratory flume simulations

Analyses of results from laboratory flume experiments are presented in which flow within gaps in canopies of flexible, submerged aquatic vegetation simulations is investigated. The aims of the work are (a) to identify the different flow regimes that may be found within such gaps, using Morris’ classical definitions of skimming flow, wake interference flow and isolated roughness flow as a template, (b) to determine the parameter space in which those flow regimes are most consistently delineated, and (c) to provide quantitative measurements of the loci of each flow regime within that parameter space for these experiments. The sedimentary and biological implications of each flow regime are also discussed. The results show that five flow regimes may be identified, expanding on Morris’ original set of three. The five are: (i) skimming flow; (ii) recirculation flow; (iii) boundary layer recovery; (iv) canopy through-flow; and (v) isolated roughness flow, the last being assumed to occur in some cases though it is not directly observed in these experiments. A Reynolds number based on the canopy overflow speed and the gap depth, and the gap aspect ratio are found to be the key parameters that determine these flow regimes, though a Froude number is found to be important for determining bed shear stress, and the length of leaves overhanging the gap from the upstream canopy is found to be important in determining the location of flow recirculation cells within the gap.     

高寒山地甘肃臭草斑块特征与土壤水分的关系

选取石羊河上游高寒退化草地甘肃臭草种群,分析斑块特征、土壤水分以及二者之间的相互关系,以期了解植被对土壤水分变化的适应性和反作用特征.结果表明:甘肃臭草斑块在形成、扩散和稳定阶段种群高度、密度、盖度由中心向边缘逐渐减小,衰退阶段由中心向边缘逐渐增大,甘肃臭草斑块内土壤水分由中心向边缘逐渐增大,斑块内土壤水分大于斑块外;斑块内土壤水分与植被主要生物学特征之间由负相关向正相关过渡.甘肃臭草具有明显的可塑性,通过对干旱区环境的长期适应,改变了小范围内土壤水分分布结构,进而形成繁茂生长的单一优势种群斑块.表3参16     

Should heterogeneity be the basis for conservation? Grassland bird response to fire and grazing.

In tallgrass prairie, disturbances such as grazing and fire can generate patchiness across the landscape, contributing to a shifting mosaic that presumably enhances biodiversity. Grassland birds evolved within the context of this shifting mosaic, with some species restricted to one or two patch types created under spatially and temporally distinct disturbance regimes. Thus, management-driven reductions in heterogeneity may be partly responsible for declines in numbers of grassland birds. We experimentally altered spatial heterogeneity of vegetation structure within a tallgrass prairie by varying the spatial and temporal extent of fire and by allowing grazing animals to move freely among burned and unburned patches (patch treatment). We contrasted this disturbance regime with traditional agricultural management of the region that promotes homogeneity (traditional treatment). We monitored grassland bird abundance during the breeding seasons of 2001-2003 to determine the influence of altered spatial heterogeneity on the grassland bird community. Focal disturbances of patch burning and grazing that shifted through the landscape over several years resulted in a more heterogeneous pattern of vegetation than uniform application of fire and grazing. Greater spatial heterogeneity in vegetation provided greater variability in the grassland bird community. Some bird species occurred in greatest abundance within focally disturbed patches, while others occurred in relatively undisturbed patches in our patch treatment. Henslow's Sparrow, a declining species, occurred only within the patch treatment. Upland Sandpiper and some other species were more abundant on recently disturbed patches within the same treatment. The patch burn treatment created the entire gradient of vegetation structure required to maintain a suite of grassland bird species that differ in habitat preferences. Our study demonstrated that increasing spatial and temporal heterogeneity of disturbance in grasslands increases variability in vegetation structure that results in greater variability at higher trophic levels. Thus, management that creates a shifting mosaic using spatially and temporally discrete disturbances in grasslands can be a useful tool in conservation. In the case of North American tallgrass prairie, discrete fires that capitalize on preferential grazing behavior of large ungulates promote a shifting mosaic of habitat types that maintain biodiversity and agricultural productivity.     

Influence of planform geometry and momentum ratio on thermal mixing at a stream confluence with a concordant bed

The effects of planform geometry and momentum flux ratio on thermal mixing at a stream confluence with concordant bed morphology are investigated based on numerical simulations that can capture the dynamics of large-scale turbulence. In two simulations, the bathymetry and asymmetrical planform geometry are obtained from field experiments and the momentum flux ratio is set at values of one and four. These two conditions provide the basis for studying differences in thermal mixing processes at this confluence when the wake mode and the Kelvin–Helmholtz mode dominate the development of coherent structures within the mixing interface (MI). The effects of channel curvature and angle between the two incoming streams on thermal mixing processes are investigated based on simulations conducted with modified planform geometries. Two additional simulations are conducted for the case where the upstream channels are parallel but not aligned with the downstream channel and for the zero-curvature case where the upstream channels are parallel and aligned with the downstream channel. The simulations highlight the influence of large-scale coherent structures within the MI and of streamwise-oriented vortical (SOV) cells on thermal mixing processes within the confluence hydrodynamics zone. Simulation results demonstrate the critical role played by the SOV cells in promoting large-scale thermal mixing for cases when such cells form in the immediate vicinity of the MI and in modifying the shape of the thermal MI within cross sections of the downstream channel—predictions consistent with empirical measurements of thermal mixing at the confluence. The set of numerical simulations reveal that the degree of thermal mixing occurring within the confluence hydrodynamic zone varies dramatically with planform geometry and incoming flow conditions. In some cases thermal mixing at the downstream end of the confluence hydrodynamic zone is limited to the MI and its immediate vicinity, whereas in others substantial thermal mixing has occurred over most of the cross-sectional area of the flow. Overall, the simulations highlight the flow conditions and the controls of these conditions that influence mixing within the immediate vicinity of a confluence.     

The three-dimensionality of obstructed shear flows

Obstructed shear flows (i.e. those over permeable media) are common in the environment. An archetypal example, flow over a submerged vegetation canopy, is investigated here. Like any flow through complex geometry, canopy flows are characterised by strong spatial gradients. The focus of this experimental study is the three-dimensionality of aquatic canopy flow, in particular that of the coherent interfacial vortices that govern mixing into and out of the canopy. It is shown here that the vortices have a finite lateral scale that is comparable to their vertical scale; both are of the order of the drag length scale of the canopy, (C _D a)⁻¹, where a is the frontal area density and C _D is a bulk drag coefficient. The finite lateral extent of the vortices generates strong lateral hydrodynamic gradients, both instantaneously and in the long-term. The instantaneous gradients, which can contribute greatly to the dispersion of dissolved and particulate species, are far more pronounced. Finally, the potential for canopies to generate differential roughness secondary circulation is examined. In the consideration of vertical scalar transport, this circulation can be of the same order as turbulent diffusion.     

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.     

Using network centrality measures to manage landscape connectivity

We use a graph-theoretical landscape modeling approach to investigate how to identify central patches in the landscape as well as how these central patches influence (1) organism movement within the local neighborhood and (2) the dispersal of organisms beyond the local neighborhood. Organism movements were theoretically estimated based on the spatial configuration of the habitat patches in the studied landscape. We find that centrality depends on the way the graph-theoretical model of habitat patches is constructed, although even the simplest network representation, not taking strength and directionality of potential organisms flows into account, still provides a coarse-grained assessment of the most important patches according to their contribution to landscape connectivity. Moreover, we identify (at least) two general classes of centrality. One accounts for the local flow of organisms in the neighborhood of a patch, and the other accounts for the ability to maintain connectivity beyond the scale of the local neighborhood. Finally, we study how habitat patches with high scores on different network centrality measures are distributed in a fragmented agricultural landscape in Madagascar. Results show that patches with high degree and betweenness centrality are widely spread, while patches with high subgraph and closeness centrality are clumped together in dense clusters. This finding may enable multispecies analyses of single-species network models.     

Statistical detection of multiscale landscape patterns

Detection of discontinuities in landscape patterns is a crucial problem both in ecology and in environmental sciences since they may indicate substantial scale changes in generating and maintaining processes of landscape patches. This paper presents a statistical procedure for detecting distinct scales of pattern for irregular patch mosaics using fractal analysis. The method suggested is based on a piecewise regression model given by fitting different regression lines to different ranges of patches ordered according to patch size (area). Proper shift-points, where discontinuities occur, are then identified by means of an iterative procedure. Further statistical tests are applied in order to verify the statistical significance of the best models selected. Compared to the method proposed by Krummel et al. (1987), the procedure described here is not influenced by subjective choices of initial parameters. The procedure was applied to landscape pattern analysis of irregular patch mosaics (CORINE biotopes) of a watershed within the Map of the Italian Nature Project. Results for three different CORINE patch types are herein presented revealing different scaling properties with special pattern organizations linked to ecological traits of vegetation communities and human disturbance.     

Inferring Metapopulation Dynamics from Patch-Level Incidence of Florida Scrub Plants

Spatial structure and dynamics of multiple populations may explain species distribution patterns in patchy communities with heterogeneous disturbance regimes, especially when species have poor dispersal. The endemic-rich Florida (U.S.A.) rosemary scrub occupies about 4% of the west portion of Archbold Biological Station and occurs scattered within a matrix of less xeric vegetation. Longer fire-return times and higher frequency of open patches in rosemary scrub provide favorable habitat for many plant species. Occupancy of 123 species of vascular plants and ground lichens in 89 patches was determined by repeated site surveys. About two-thirds of the species occurring at more than 14 patches had a significant logistic regression of presence on time-since-fire, patch size, patch isolation, or their interactions. Species with presence related to the interaction between patch isolation and patch size were primarily herbs and small shrubs specializing in rosemary scrub. These results suggest the importance of spatial characteristics of the landscape for population turnover of these species. An incidence-based metapopulation model was used to predict extinction and colonization probabilities of those species with presence in rosemary scrub patches related to the studied spatial variables. This is the first attempt to apply incidence-based metapopulation models to plants. The results showed stronger effects of patch size and patch isolation on extinction probabilities of herbs than on those of woody species. Because of their effect on spatial heterogeneity and habitat availability, fire suppression and habitat destruction may decrease persistence probabilities for these rosemary scrub specialists, many of which are endangered species.     

Exchange flow between open water and floating vegetation

This study describes the exchange flow between a region with open water and a region with a partial-depth porous obstruction, which represents the thermally-driven exchange that occurs between open water and floating vegetation. The partial-depth porous obstruction represents the root layer, which does not penetrate to the bed. Initially, a vertical wall separates the two regions, with fluid of higher density in the obstructed region and fluid of lower density in the open region. This density difference represents the influence of differential solar heating due to shading by the vegetation. For a range of root density and root depths, the velocity distribution is measured in the lab using PIV. When the vertical wall is removed, the less dense water flows into the obstructed region at the surface. This surface flow bifurcates into two layers, one flowing directly through the root layer and one flowing beneath the root layer. A flow directed out of the vegetated region occurs at the bed. A model is developed that predicts the flow rates within each layer based on energy considerations. The experiments and model together suggest that at time- and length-scales relevant to the field, the flow structure for any root layer porosity approaches that of a fully blocked layer, for which the exchange flow occurs only beneath the root layer.     

Evaluating the dimensionality and significance of “active periods” in turbulent environmental flows defined using Lipshitz/Hölder regularity

Active periods within perturbed boundary-layer flows are considered in terms of the local roughness of measured velocity time series and defined in terms of Hölder/Lipshitz exponents. Such events are associated with the passage of energetic, coherent flow structures and are responsible for exerting high turbulent stresses because of the rapid changes in velocity that occur at such times. A method is proposed for assessing the effective dimensionality of such active periods, as well as their significance to the flow field, for a particular choice of flow metric. The method is applied to the turbulent flow through a confluence flow geometry, with velocity samples acquired close to the bed of the channel in a zone of complex mixing. The dimensionality of the active periods is consistent with the observed patterns of sediment entrainment from the bed, with the significance of the active periods decaying away from the erosional zone.     

Does the canopy mixing layer model apply to highly flexible aquatic vegetation? Insights from numerical modelling

Vegetation is a characteristic feature of shallow aquatic flows such as rivers, lakes and coastal waters. Flow through and above aquatic vegetation canopies is commonly described using a canopy mixing layer analogy which provides a canonical framework for assessing key hydraulic characteristics such as velocity profiles, large-scale coherent turbulent structures and mixing and transport processes for solutes and sediments. This theory is well developed for the case of semi-rigid terrestrial vegetation and has more recently been applied to the case of aquatic vegetation. However, aquatic vegetation often displays key differences in morphology and biomechanics to terrestrial vegetation due to the different environment it inhabits. Here we investigate the effect of plant morphology and biomechanical properties on flow–vegetation interactions through the application of a coupled LES-biomechanical model. We present results from two simulations of aquatic vegetated flows: one assuming a semi-rigid canopy and the other a highly flexible canopy and provide a comparison of the associated flow regimes. Our results show that while both cases display canopy mixing layers, there are also clear differences in the shear layer characteristics and turbulent processes between the two, suggesting that the semi-rigid approximation may not provide a complete representation of flow–vegetation interactions.     

Connectivity or demography: Defining sources and sinks in coral reef fish metapopulations

The identity of an individual patch as a source or a sink within a metapopulation is a function of its ability to produce individuals and to disperse them to other patches. In marine systems patch identity is very often defined by dispersal ability alone—upstream patches are sources—while issues of variable habitat quality (which affects local production) are ignored. This can have important ramifications for the science of marine reserve siting. This study develops a spatially explicit source–sink metapopulation model for reef fish and uses it to evaluate the relative importance of connectivity versus demography and how this depends upon the level of local larval retention and the strength of density-dependent recruitment. Elasticity analyses indicated that patch contribution (source or sink) was more sensitive to demographic parameters (particularly survival) than connectivity and this effect was conserved even under strong levels of density-dependence and was generally strengthened as local retention increased. Variability in the relationship between parameter elasticity and local retention was shown to be dependent upon the magnitude of connectivity for an individual patch relative to a critical connectivity value. The proportion of larvae lost due to transport processes was an important parameter which directly affected the magnitude of this critical connectivity value. Patches with connectivity values less than the critical value contributed to the metapopulation largely via production (i.e., local demographics most important). As local retention increased, so did the importance of demographic parameters in these patches. Patches with connectivity values greater than the critical value contributed largely via dispersal of larvae and thus the importance of local demographics decreased as local retention increased.