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.     

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.     

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.     

A σ-coordinate transport model coupled with rotational boussinesq-type equations

This paper describes a σ-coordinate scalar transport model coupled with a Boussinesq-type hydrodynamic model. The Boussinesq model has the ability to calculate both three-dimensional velocity distributions and the water surface motion. To capture ‘dispersion’ processes in open channel flow, horizontal vorticity effects induced by a bottom shear stress are included in the Boussinesq model. Thus, a reasonable representation of vertical flow structure can be captured in shallow and wavy flow fields. To solve the coupled Boussinesq and scalar transport system, a finite-volume method, based on a Godunov-type scheme with the HLL Riemann solver, is employed. Basic advection and advection–diffusion numerical tests in a non-rectangular domain were carried out and the computed results show good agreement with analytic solutions. With quantitative comparisons of dispersion experiments in an open channel, it is verified that the proposed coupled model is appropriate for both near and far field scalar transport predictions. From numerical simulations in the surf zone, physically reasonable results showing expected vertical variation are obtained.     

Spatial decoupling of facilitation and competition at the origin of gapped vegetation patterns

Spatially periodic vegetation patterns, forming gaps, bands, labyrinths, or spots, are characteristic of arid and semiarid landscapes. Self-organization models can explain this variety of structures within a unified conceptual framework. All these models are based on the interplay of positive and negative effects of plants on soil water, but they can be divided according to whether they assume the interactions to be mediated by water redistribution through runoff/diffusion or by plants' organs. We carried out a multi-proxy approach of the processes operating in a gapped pattern in southwest Niger dominated by a shrub species. Soil moisture within the root layer was monitored in time and space over one month of the rainy season. Soil water recharge displayed no spatial variation with respect to vegetation cover, but the stock half-life under cover was twice that of bare areas. A kernel of facilitation by the aboveground parts of shrubs was parameterized, and soil water half-life was significantly correlated to the cumulated facilitative effects of shrubs. The kernel range was found to be smaller than the canopy radius (81%). This effect of plants on soil water dynamics, probably through a reduction of evaporation by shading, is shown to be a better explanatory variable than potentially relevant soil and topography parameters. The root systems of five individuals of Combretum micranthum G. Don were excavated. Root density data were used as a proxy to parameterize a kernel function of interplant competition. The range of this kernel was larger than the canopy radius (125%). The facilitation-to-competition range ratio, reflecting the above-to-belowground ratio of plant lateral extent, was smaller than 1 (0.64), a result supporting models assuming that patterning may emerge from an adaptation of plant morphology to aridity and shallow soils by means of an extended lateral root system. Moreover, observed soil water gradients had directions opposite to those assumed by alternative mathematical models based on underground water diffusion. This study contributes to the growing awareness that combined facilitative and competitive plant interactions can induce landscape-scale patterns and shape the two-way feedback loops between environment and vegetation.     

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.     

Stokes drift through corals

Motivated by shallow ocean waves propagating over coral reefs, we investigate the drift velocities due to surface wave motion in an effectively inviscid fluid that overlies a saturated porous bed of finite depth. Previous work in this area either neglects the large-scale flow between layers (Phillips in Flow and reactions in permeable rocks, Cambridge University Press, Cambridge, 1991) or only considers the drift above the porous layer (Monismith in Ann Rev Fluid Mech 39:37–55, 2007). Overcoming these limitations, we propose a model where flow is described by a velocity potential above the porous layer and by Darcy’s law in the porous bed, with derived matching conditions at the interface between the two layers. Both a horizontal and a novel vertical drift effect arise from the damping of the porous bed, which requires the use of a complex wavenumber k. This is in contrast to the purely horizontal second-order drift first derived by Stokes (Trans Camb Philos Soc 8:441–455, 1847) when working with solely a pure fluid layer. Our work provides a physical model for coral reefs in shallow seas, where fluid drift both above and within the reef is vitally important for maintaining a healthy reef ecosystem (Koehl et al. In: Proceedings of the 8th International Coral Reef Symposium, vol 2, pp 1087–1092, 1997; Monismith in Ann Rev Fluid Mech 39:37–55, 2007). We compare our model with field measurements by Koehl and Hadfield (J Mar Syst 49:75–88, 2004) and also explain the vertical drift effects as documented by Koehl et al. (Mar Ecol Prog Ser 335:1–18, 2007), who measured the exchange between a coral reef layer and the (relatively shallow) sea above.

     

Experimental investigation of channel flow through idealized isolated tree-like vegetation

Environmental Fluid Mechanics - Riparian and floodplain tree-like emergent vegetation alter significantly the flow field and lead to complicated three-dimensional flow patterns, characterized by...     

Simulation of exchange flow between open water and floating vegetation using a modified RNG <Emphasis Type="Italic">k</Emphasis>-<Emphasis Type="Italic">ε</Emphasis> turbulence model

The paper focuses on the numerical simulation of the exchange flow between open water and floating vegetation, which plays an important role in maintaining the ecological balance by transporting nutrient matter. The simulation was conducted using a new solver developed upon OpenFOAM. A modified RNG k-ε turbulence model, which is expected to model both the high- and low-Reynolds number flows correctly, was used to determine the eddy viscosity. Several particular terms were added into the momentum equations and turbulence model equations to model the effects of vegetation and buoyancy. Among these terms, the term for the effect of vegetation in the ε-equation was re-modelled. The model was validated by properly predicting the profiles of mean velocity and turbulent kinetic energy for flows through suspended canopies. The density flow between open and vegetated water was simulated with the same conditions as those of the experiment conducted by Zhang and Nepf. The predicted results agreed well with the experimental data and provided more detailed information of such exchange flow. The convection between the root layer and the layer beneath the roots, which was not observed in the experiment, was observed in the numerical simulation.     

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 solution of two-dimensional advection–dispersion equation with spatio-temporal coefficients for point sources in an infinite medium using Green’s function method

In the present study analytical solutions of a two-dimensional advection–dispersion equation (ADE) with spatially and temporally dependent longitudinal and lateral components of the dispersion coefficient and velocity are obtained using Green’s Function Method (GFM). These solutions describe solute transport in infinite horizontal groundwater flow, assimilating the spatio-temporal dependence of transport properties, dependence of dispersion coefficient on velocity, and the particulate heterogeneity of the aquifer. The solution is obtained in the general form of temporal dependence and the source term, from which solutions for instantaneous and continuous point sources are derived. The spatial dependence of groundwater velocity is considered non-homogeneous linear, whereas the dispersion coefficient is considered proportional to the square of spatial dependence of velocity. An asymptotically increasing temporal function is considered to illustrate the proposed solutions. The solutions are validated with the existing solutions derived from the proposed solutions in three special cases. The effect of spatially/temporally dependent heterogeneity on the solute transport is also demonstrated. To use the GFM, the ADE with spatio-temporally dependent coefficients is reduced to a dispersion equation with constant coefficients in terms of new position variables introduced through properly developed coordinate transformation equations. Also, a new time variable is introduced through a known transformation.     

Modelling the spatial distribution of montane and subalpine forests in the central Alps using digital elevation models

The analysis of complex interactions between spatial distribution patterns of site factors and vegetation types is crucial for understanding high mountain ecosystems, especially in the view of a changing climate. Therefore, in the present study, a GIS and remote sensing-based approach is followed to produce a vegetation map for a study area in the Western Alps (Switzerland). Two major forest alliances are chosen for analysis: subalpine coniferous forest Vaccinio-Piceion/Larici-Pinetum cembrae and montane oak forest Quercion pubescenti-petraeae. As spatial information on site factors is commonly lacking in mountain areas, the use of a digital elevation model (DEM) is a potential substitute for use in vegetation analyses: it highly correlates with temperature, moisture, geomorphological processes and disturbance factors. Thus, it is important to analyse the capabilities of a DEM for indicating habitat conditions in a landscape characterised by high topodiversity and a patchwork of microclimatic habitats.For the purpose of identifying the potential of landform parameters for the indication of forest habitat structures in the present study, 24 primary and secondary landform parameters have been derived, indicating temperature and moisture distribution, exposure towards wind, snow, etc. Quantitative analyses were performed using statistical means such as contingency correlation coefficients and principal components analysis. The results formed the basis for the development of parallel-epiped-vegetation models (PED) used to simulate the spatial distribution patterns of the subalpine coniferous and the montane oak forest. It can be shown that topographic variables derived from a DEM at a spatial resolution of 25 m are very useful for indicating habitats of large forest types. Additionally potential forest sites in the cultural landscape, removed by human logging, can be reconstructed.Inaccuracies within the simulation results can partly be attributed to the insufficient parameterisation of geomorphologic activity and to poor spatial resolution of the DEM as compared to the vegetation data. Although the lack of information on the human dimension leads to some uncertainties in the interpretation of spatial patterns of vegetation, the exclusive use of topographic variables in vegetation models for the indication of forest habitats is very promising.     

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.     

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.     

Modification of a stream temperature model with Beer's law and application to GaoShan Creek in Taiwan

A physics-based stream temperature model [Tung, C.P., Lee, T.Y., Yang, Y.C., 2006. Modelling climate-change impacts on stream temperature of Formosan Landlocked Salmon habitat. Hydrol. Process. 20, 1629–1649] was improved by incorporating shading effects caused by both cliff terrain and riverbank dense vegetation to simulate hourly stream temperature variations in 1 day. Daily maximal stream temperature is a critical factor to the habit distribution of the Formosan Landlocked Salmon, an important and endangered species. Currently, it only can be found in ChiChiaWan Creek and GaoShan Creek in Taiwan. The former stream temperature model only considers the shading effects of cliff terrain and works well for ChiChiaWan Creek, but overestimates stream temperatures of GaoShan Creek having dense riverbank vegetative covers. The model was modified with the Beer's law and a parameterization scheme to describe the diminishing of the incident solar radiation to take vegetative shading effects into account. Simulation results of GaoShan Creek show the success of this improvement. The shading effects induced by both terrain and vegetation can significantly affect stream temperature distributions. Simulation experiments were conducted to indicate shading effects are varied in different watersheds and seasons.     

Sensitivity of modelled gross primary productivity to topographic effects on surface radiation: A case study in the Cotter River Catchment, Australia

Gross primary productivity (GPP) is a critical response variable for many environmental problems, including terrestrial carbon accounting and the calculation of catchment water balances. Various approaches for modelling GPP have been developed and applied at continental and landscapes scales, but little attention has been given to the sensitivity of GPP to the spatial scale of its driving variables. A key driving variable is surface radiation (R_s) which is influenced by both meso-scale factors (latitude, time of year, cloudiness) and the topographic variables of slope, aspect and horizon shading. We compared the sensitivity of modelled GPP to two different sources of surface radiation (R_s): (1) the ANUCLIM method which only captures the meso-scaled factors; and (2) the SRAD method which incorporates the topographic effects GPP was calculated using the radiation use efficiency (RUE) model (Roderick et al., 2001) to discern general patterns of vegetation productivity at a sub-catchment (i.e. sub-water shed) scale. The radiation use efficiency approach uses the normalized difference vegetation index (NDVI) derived from satellite data (MODIS TERRA), along with estimates of solar radiation at the top of the atmosphere (R_o) and canopy (R_s). In this approach, R_o and R_s capture the influence of diffuse irradiance in canopy photosynthesis and vegetation productivity respectively. This research showed that R_s calculated using the SRAD program provides important discrimination of GPP regimes at a sub-catchment scale, as the result of minimum and maximum daily radiation varying between shaded and exposed surfaces. However, mean daily radiation at a whole-of-catchment scale did not differ between the two sources as the differences in the minimum and maximum daily values tend to cancel each other out. Applications of GPP models therefore need to consider whether topographic factors are important and select the appropriate source of R_s values. GPP models should also reflect understanding of radiation use efficiency. However, further research is required especially with respect to the influence of water stress on plant response.     

A mathematical model for type II profile of concentration distribution in turbulent flows

This paper presents a mathematical model to investigate type II profile of suspension concentration distribution (i.e., the concentration profile where the maximum concentration appears at some distance above the bed surface) in a steady, uniform turbulent flow through open-channels. Starting from the mass and momentum conservation equations of two-phase flow, a theoretical model has been derived. The distribution equation is derived considering the effects of fluid lift force, drag force, particle inertia, particle–particle interactions, particle velocity fluctuations and drift diffusion. The equation is solved numerically and is compared with available experimental data as well as with other models existing in the literature. Good agreement between the observed value and computed result, and minimum error in comparison to other models indicate that the present model can be applied in predicting particle concentration distribution for type II profile for a wide range of flow conditions. The proposed model is also able to show the transition from type I profile to type II profile.     

Turbulent secondary flows in wall turbulence: vortex forcing,scaling arguments,and similarity solution

Spanwise surface heterogeneity beneath high-Reynolds number, fully-rough wall turbulence is known to induce a mean secondary flow in the form of counter-rotating streamwise vortices—this arrangement is prevalent, for example, in open-channel flows relevant to hydraulic engineering. These counter-rotating vortices flank regions of predominant excess(deficit) in mean streamwise velocity and downwelling(upwelling) in mean vertical velocity. The secondary flows have been definitively attributed to the lower surface conditions, and are now known to be a manifestation of Prandtl’s secondary flow of the second kind—driven and sustained by spatial heterogeneity of components of the turbulent (Reynolds averaged) stress tensor (Anderson et al. J Fluid Mech 768:316–347, 2015). The spacing between adjacent surface heterogeneities serves as a control on the spatial extent of the counter-rotating cells, while their intensity is controlled by the spanwise gradient in imposed drag (where larger gradients associated with more dramatic transitions in roughness induce stronger cells). In this work, we have performed an order of magnitude analysis of the mean (Reynolds averaged) transport equation for streamwise vorticity, which has revealed the scaling dependence of streamwise circulation intensity upon characteristics of the problem. The scaling arguments are supported by a recent numerical parametric study on the effect of spacing. Then, we demonstrate that mean streamwise velocity can be predicted a priori via a similarity solution to the mean streamwise vorticity transport equation. A vortex forcing term has been used to represent the effects of spanwise topographic heterogeneity within the flow. Efficacy of the vortex forcing term was established with a series of large-eddy simulation cases wherein vortex forcing model parameters were altered to capture different values of spanwise spacing, all of which demonstrate that the model can impose the effects of spanwise topographic heterogeneity (absent the need to actually model roughness elements); these results also justify use of the vortex forcing model in the similarity solution.     

植被-侵蚀动力学模型参数的确定及在黄土高原的应用

植被-侵蚀动力学是研究流域植被与侵蚀在人类活动影响下演变规律的一门新的边缘学科。由植被-侵蚀动力学模型可绘制植被-侵蚀状态图,以预测停止人为干扰后植被及侵蚀变化趋势,评价水土保持措施成效,提出小流域优化治理方略。植被-侵蚀动力学模型4个参数a、c、b、f是流域重要特征参数,其值的确定是进行植被-侵蚀动力学模拟及绘制植被-侵蚀状态图的前提和基础。它们依赖于流域气候、土壤和地形地貌。本文以黄土高原典型小流域为例,采用试算法确定了它们的取值,通过相关分析发现参数a与P、T,参数c、b、f与D50、S之间存在明显的相关关系,并得到了相关关系式。将其用于吕二沟小流域植被-侵蚀动力学模拟,结果与实际情况符合较好。绘制了各参数的取值分布图。最后以燕儿沟小流域为例,将以上相关关系式用于植被-侵蚀状态图的绘制,及小流域自然状况、植被-侵蚀状态及水土保持效果的评价。     

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.