Impulsive waves caused by subaerial landslides

This paper presents the experimental results of impulsive waves caused by subaerial landslides. A wide range of effective parameters are considered and studied by performing 120 laboratory tests. Considered slide masses are both rigid and deformable. The effects of bed slope angle, water depth, slide impact velocity, geometry, shape and deformation on impulse wave characteristics have been inspected. The impulse wave features such as amplitude, period and also energy conversation are studied. The effects of slide Froude number and deformation on energy conversation from slide into wave are also investigated. Based on laboratory measured data an empirical equation for impulse wave amplitude and period have been presented and successfully verified using available data of previous laboratory works.     

Curvature-induced secondary flow in submarine channels

The curvature-driven secondary flow in sinuous submarine channels has been a subject of considerable interest and controversy. Here, results from numerical model studies involving saline flow in laboratory-scale channels are presented. A 3D finite volume model of density and turbidity currents is used and simulations are run with different inflow discharges and channel-axis slopes. The simulation results show strong influence of bend wave length, channel gradient, confinement and cross sectional shape on the structure of secondary flow in submarine channels. Major findings are: (i) reversal of secondary flow in submarine channels is strongly associated with a tight bend characterized by a smaller wave length to width ratio or larger wave number, (ii) for the same inflow condition and planform characteristics, a trapezoidal channel cross section is more favorable to secondary flow reversal than a rectangular cross section, (iii) lateral convection resulting from the interaction between in-channel and overbank flows leads to the reversal of secondary flow in an unconfined channel at a lower channel slope than in a confined channel with the same dimensions, (iv) flow discharge has only nominal effect on the secondary flow in submarine channels.     

Modeling the buoyant flow of heated water discharged from surface and submerged side outfalls in shallow and deep water with a cross flow

In this study, a three-dimensional model was used to numerically study the buoyant flow, along with its mixing characteristics, of heated water discharged from the surface and submerged side outfalls in shallow and deep water with a cross flow. Hydraulic experimental data were used to evaluate the applicability of the model. The simulation results agree well with the experimental results, particularly, the jet trajectories, the dimensions of the recirculating zone, and the distribution of the dimensionless excess temperature. The level of accuracy of the simulation results of the present study is nearly identical to that of the results conducted by McGuirk and Rodi (1978). If the heated water is discharged into shallow water where the momentum flux ratio and the discharge densimetric Froude number are high, the submerged discharge method is better than the surface discharge method in terms of the scale of the recirculating zone and the minimum dilution. In deep water, where the momentum flux ratio and discharge densimetric Froude number are low, however, the submerged discharge method had few advantages. In shallow water, the discharge jet is deflected by the ambient cross flow, while forcing the ambient flow to bend towards the far bank for the full depth. For a submerged discharge in shallow water, the recirculating zone is the largest in the lowest layer but becomes smaller in the upper layer. As the water depth increases, the ambient flow goes over the jet and diminishes the blocking effect, thereby decreasing the bending of the jet.     

Formation of tidal starting-jet vortices through idealized barotropic inlets with finite length

This paper presents a surface particle image velocimetry study to investigate the dynamics of shallow starting-jet dipoles formed by tidal flow through inlets and their interaction with vorticity formed at the inlet channel lateral boundaries. Vortical structure in the flow field is identified using a local swirl strength criterion evaluated from the two-dimensional flow field. The starting jet dipole vortices and vortices formed as the lateral boundary layers are expelled during flow reversal are characterized by their trajectory, size, and circulation. Using these quantities, a model is developed to predict the size and strength of the expelled lateral boundary layer vortices based on the inlet velocity, channel length, and width of the lateral boundary layer. The expelled boundary layer vortices are found to disrupt the formation of the primary tidal jet dipole through two mechanisms. First, because the boundary layer vortices themselves form a dipole with each half of the starting-jet dipole, the starting-jet vortices are pulled apart and advected away from the inlet mouth early in the tidal cycle, resulting in a reduction in the spin-up time and the amount of vorticity input during starting-jet vortex formation. Second, the advection of each dipole away from the inlet disconnects each starting-jet vortex from the starting jet; hence, the vortices are not fed by fluid in the jet or energized by shear in the jet boundary layers. These influences of the lateral boundary layer on the starting-jet vortices’ formation and propagation are found to be a function of the channel length L, maximum velocity U, and tidal period T, resulting in a predictive value to characterize their trajectory, strength, and evolution.     

Initial mixing of inclined dense jet in perpendicular crossflow

A comprehensive experimental investigation for an inclined ( $60^{\circ }$ to vertical) dense jet in perpendicular crossflow—with a three-dimensional trajectory—is reported. The detailed tracer concentration field in the vertical cross-section of the bent-over jet is measured by the laser-induced fluorescence technique for a wide range of jet densimetric Froude number $Fr$ and ambient to jet velocity ratios $U_r$ . The jet trajectory and dilution determined from a large number of cross-sectional scalar fields are interpreted by the Lagrangian model over the entire range of jet-dominated to crossflow-dominated regimes. The mixing during the ascent phase of the dense jet resembles that of an advected jet or line puff and changes to a negatively buoyant thermal on descent. It is found that the mixing behavior is governed by a crossflow Froude number $\mathbf{F} = U_r Fr$ . For $\mathbf{F} < 0.8$ , the mixing is jet-dominated and governed by shear entrainment; significant detrainment occurs and the maximum height of rise $Z_{max}$ is under-predicted as in the case of a dense jet in stagnant fluid. While the jet trajectory in the horizontal momentum plane is well-predicted, the measurements indicate a greater rise and slower descent. For $\mathbf{F} \ge 0.8$ the dense jet becomes significantly bent-over during its ascent phase; the jet mixing is dominated by vortex entrainment. For $\mathbf{F} \ge 2$ , the detrainment ceases to have any effect on the jet behavior. The jet trajectory in both the horizontal momentum and buoyancy planes are well predicted by the model. Despite the under-prediction of terminal rise, the jet dilution at a large number of cross-sections covering the ascent and descent of the dense jet are well-predicted. Both the terminal rise and the initial dilution for the inclined jet in perpendicular crossflow are smaller than those of a corresponding vertical jet. Both the maximum terminal rise $Z_{max}$ and horizontal lateral penetration $Y_{max}$ follow a $\mathbf{F}^{-1/2}$ dependence in the crossflow-dominated regime. The initial dilution at terminal rise follows a $S \sim \mathbf{F}^{1/3}$ dependence.     

Instabilities in non-rotating and rotating shallow shear flows

Numerical simulations for the wave radiation effect on the linear and nonlinear instabilities of rotating and non-rotating shallow flows are conducted using shallow-water equations. At a low convective Froude number, the results of the instabilities is a string of eddies. The coalescence between the neighbouring eddies decides the transverse mixing of the shallow shear flow. At a higher convective Froude number, the development of the shear flow is characterized by wave radiation and the production of shocklets. The radiation of waves in the non-rotating shallow flow is a phenomenon analogous to the radiation of sound in gas dynamics. In the rotating flow on the other hand, the shallow-flow instabilities are intensified due to rotational interference within a window of instability over a narrow range of Rossby numbers.     

The generation of internal waves by tidal flow over continental shelf/slope topography

Experiments were performed to investigate the generation of internal waves induced by a barotropic tidal flow over continental shelf/slope topography in a two-layer stratified fluid. The interaction of the barotropic tide with the continental slope resulted in the formation of both linear and nonlinear waveforms, ranging from a linear wave of depression to a highly nonlinear internal bolus. The type of wave response was strongly dependent on the tidal forcing intensity and the position of the density interface relative to the shelf depth. Based on the values of the internal Froude number and the layer depth ratio, we delineate four distinct generation regimes, each with a distinct wave response.     

Effects of artificial eutrophication on the metabolism of the Japanese oyster Crassostrea gigas

The morphology of the primary tube feet in 15 species of comatulid (unstalked) crinoids from coral reefs in the Palau Islands and Lizard Island, Great Barrier Reef, was investigated using close-up underwater photographs of the tube feet taken in the natural habitat. Measurements of length of the tube feet and their spacing along the pinnule were taken from these photographs. Tube feet of species of the family Comasteridae have a mean length of 0.75 mm and a mean spacing of 6 tube feet mm^-1. Tube feet of non-comasterid species occurring in the same environments have a mean length of 0.55 mm and a mean spacing of about 8 tube feet mm^-1. The relationship between spacing and length of the tube feet is highty significant for both the Palau and Lizard Island samples (P<0.01). Species having longer and more widely spaced tube feet live partly concealed within the infrastructure of the reef and hold the arms and pinnules in a multidirectional posture. Species having shorter and more closely spaced tube feet perch on top of reef pinnacles or alcyonarians and form planar filtration fans normal to unidirectional currents or wave oscillations. Longer and more widely spaced tube feet in species dwelling within the reef infrastructure provide more efficient filtration in the slow, meandering flow prevailing there. Closer spacing of the tube feet in species exposed to near-mainstream flow provides a more efficient filtration mechanism at higher flow velocities. Reduced length of the tube feet in these species may be a consequence of closer spacing of the pinnules. Differentiation of these co-occurring species in spacing and length of the tube feet implies differentiation in food particles captured. This may in some cases constitute resource partitioning. Most species which overlap in living habits are significantly different in spacing and length of the tube feet.     

Analysis of the velocity field in a large rectangular channel with lateral shockwave

In this work the authors describe the main characteristics of the velocity field of hydraulic jumps in a very large channel where lateral shockwaves occur. Experiments were carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). Extensive flow velocity measurements were investigated in order to have a clearer understanding of both hydraulic jump development and lateral shockwave formation in a very large channel. Eight experiments were performed in a 4m wide rectangular channel; the experiments differed in the inlet Froude number F ₀ and the jump type. Seven tests were carried out with undular jumps and one with a roller jump. The flow velocity and the flow free surface measurements were taken using a two-dimensional Acoustic Doppler Velocimeter (ADV) and an ultrasonic profiler, respectively. The experimental results can be summarized as follow: (i) the formation of well developed lateral shockwaves similar to those of oblique jumps were observed; (ii) the comparison of the experimental and theoretical data shows that the classic shockwave theory is sufficiently confirmed in the analyzed range of Reynolds number, taking into account the experimental errors and the difference between the theoretical and experimental assumptions; (iii) the transversal flow velocity profiles in the recirculating zone show a good agreement with the numerical simulations presented in literature in the case of a separated turbulent boundary layer over a flat plate. This conclusion enables us to confirm the hypothesis that the lateral shockwaves in the channel are the result of a boundary layer which, as observed, forms on the channel sidewalls.     

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.     

Froude scaling limitations in modeling of turbidity currents

The scaling problem associated with the modeling of turbidity currents has been recognized but is yet to be explored systematically. This paper presents an analysis of the dimensionless governing equations of turbidity currents to investigate the scale effect. Three types of flow conditions are considered: (i) conservative density current; (ii) purely depositional turbidity current; and (iii) mixed erosional/depositional turbidity current. Two controlling dimensionless numbers, the Froude number and the Reynolds number, appear in the non-dimensional governing equations. When densimetric Froude similarity is satisfied, the analysis shows that the results would be scale-invariant for conservative density current under the rough turbulent condition. In the case of purely depositional flows, truly scale-invariant results cannot be obtained, as the Reynolds-mediated scale effects appear in the bottom boundary conditions of the flow velocity and sediment fall velocity. However, the scale effect would be relatively modest. The Reynolds effect becomes more significant for erosional or mixed erosional/depositional turbidity currents as Reynolds-mediated scale effects also appear in the sediment entrainment relation. Numerical simulations have been conducted at three different scales by considering densimetric Froude scaling alone as well as combined densimetric Froude and Reynolds similarity. Simulation results confirm that although the scaling of densimetric Froude number alone can produce scale-invariable results for conservative density currents, variations occur in the case of turbidity currents. The results become scale invariant when densimetric Froude and Reynolds similarities are satisfied simultaneously.     

Growth of the staghorn coral Acropora formosa at Houtman Abrolhos, Western Australia

The growth (extension rate, number of radial branches, skeletal mass, branch diameter) of the␣staghorn coral Acropora formosa (Dana, 1846) was examined at four sites on the Beacon Island platform at Houtman Abrolhos, in subtropical Western Australia (28°S). Sites were at depths of 7 to 11 m, with variable exposure to weather and swell conditions. Two sites on the western reef slope were partly exposed to the oceanic swell, and two sites in the lagoon were largely protected from wave action. Linear extension rate between 1994 and 1995 varied significantly between sites, with greater linear extension at the more protected lagoonal sites. However, accumulation of skeletal mass per branch and number of newly initiated radial branches did not vary significantly between the sites. Carbonate was deposited in similar amounts, but either as porous, rapidly extending branches, or as denser branches which extended more slowly. Branch extension rate over 11.5 mo ranged from a mean of 50.3 mm (range=13 to 93 mm) at a reef slope site to a mean of 76.0 mm (range=31 to 115 mm) at a sheltered lagoonal site. Mean extension rates were almost twice that previously reported for this species in Houtman Abrolhos (37 to 43 mm yr⁻¹) from a shallower site where environmental conditions were apparently sub-optimal. Growth was within the range reported for A. formosa from tropical sites, which is consistent with the relatively high calcification and reef-accretion rates recorded for Houtman Abrolhos in geological and metabolic studies. The role of reduced coral growth-rate in limiting coral reef formation at high latitudes remains equivocal. Received: 19 November 1997 / Accepted: 5 May 1998     

Calculation of Entrainment in Density Jumps

The density jump in a two-layer channel flow of miscible fluids, in which one of the layers is infinitely deep and at rest, is analyzed using the momentum flux and mass flux conservation equations. The analysis yields simple equations relating the heights upstream and downstream of the jump with the upstream Froude number and the rate of entrainment into the moving layer, as well as a relation between the maximum possible entrainment and the upstream Froude number. The analysis also shows that when the flow down-stream of the jump is controlled by an obstruction or a contraction, the entrainment rate into the jump depends solely on the dimensionless obstruction height contraction ratio and the upstream Froude number.     

Experimental study of recirculating flows generated by lateral shock waves in very large channels

Due to the lack of data on hydraulic-jump dynamics in very large channels, the present paper describes the main characteristics of the velocity field and turbulence in a large rectangular channel with a width of 4 m. Although a hydraulic jump is always treated as a wave that is transversal to the channel wall, in the case of this study it has a trapezoidal front shape, first starting from a point at the sidewalls and then developing downstream in an oblique manner, finally giving rise to a trapezoidal shape. The oblique wave front may be regarded as a lateral shockwave that arises from a perturbation at a certain point of the lateral wall and travels obliquely toward the centreline of the channel. The experimental work was carried out at the Coastal Engineering Laboratory of the Water Engineering and Chemistry Department of the Technical University of Bari (Italy). In addition to the hydraulic jump formation, a large recirculating flow zone starts to develop from the separating point of the lateral shock wave and a separate boundary layer occurs. Intensive measurements of the streamwise and spanwise flow velocity components along one-half width of the channel were taken using a bidimensional Acoustic Doppler Velocimeter (ADV). The water surface elevation was obtained by means of an ultrasonic profiler. Velocity vectors, transversal velocity profiles, turbulence intensities and Reynolds shear stresses were all investigated. The experimental results of the separated boundary layer were compared with numerical predictions and related work presented in literature and showed good agreement. The transversal velocity profiles indicated the presence of adverse pressure gradient zones and the law of the wall appears to govern the region around the separated boundary layer.     

Integral Model for Turbulent Buoyant Jets in Unbounded Stratified Flows Part 2: Plane Jet Dynamics Resulting from Multiport Diffuser Jets

An integral model for the plane buoyant jet dynamics resulting from the interaction of multiple buoyant jet effluxes spaced along a diffuser line is considered as an extension of the round jet formulation that was proposed in Part I. The receiving fluid is given by an unbounded ambient environment with uniform density or stable density stratification and under stagnant or steady sheared current conditions. Applications for this situation are primarily for submerged multiport diffusers for discharges of liquid effluents into ambient water bodies, but also for multiple cooling tower plumes and building air-conditioning. The CorJet model formulation describes the conservation of mass, momentum, buoyancy and scalar quantities in the turbulent jet flow in the plane jet geometry. It employs an entrainment closure approach that distinguishes between the separate contributions of transverse shear and of internal instability mechanisms, and contains a quadratic law turbulent pressure force mechanism. But the model formulation also includes several significant three-dimensional effects that distinguish actual diffuser installations in the water environment. These relate to local merging processes from the individual multiple jets, to overall finite length effects affecting the plume geometry, and to bottom proximity effects given by a “leakage factor” that measures the combined affect of port height and spacing in allowing the ambient flow to pass through the diffuser line in order to provide sufficient entrainment flow for the mixing downstream from the diffuser. The model is validated in several stages: First, comparison with experimental data for the asymptotic, self-similar stages of plane buoyant jet flows, i.e. the plane pure jet, the pure plume, the pure wake, the advected line puff, and the advected line thermal, support the choice of the turbulent closure coefficients contained in the entrainment formulation. Second, comparison with data for many types of non-equilibrium flows with a plane geometry support the proposed functional form of the entrainment relationship, and also the role of the pressure force in the jet deflection dynamics. Third, the observed behavior of the merging process from different types of multiport diffuser discharges in both stagnant and flowing ambient conditions and with stratification appears well predicted with the CorJet formulation. Fourth, a number of spatial limits of applicability, relating to terminal layer formation in stratification or transition to passive diffusion in a turbulent ambient shear flow, have been proposed. In sum, the CorJet integral model appears to provide a mechanistically sound, accurate and reliable representation of complex buoyant jet mixing processes, provided the condition of an unbounded receiving fluid is satisfied.     

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.     

Large-scale turbulent structure of uniform shallow free-surface flows

The hydrodynamics of super- and sub-critical shallow uniform free-surface flows are assessed using laboratory experiments aimed at identifying and quantifying flow structure at scales larger than the flow depth. In particular, we provide information on probability distributions of horizontal velocity components, their correlation functions, velocity spectra, and structure functions for the near-water-surface flow region. The data suggest that for the high Froude number flows the structure of the near-surface layer resembles that of two-dimensional turbulence with an inverse energy cascade. In contrast, although large-scale velocity fluctuations were also present in low Froude number flow its behaviour was different, with a direct energy cascade. Based on our results and some published data we suggest a physical explanation for the observed behaviours. The experiments support Jirka’s [Jirka GH (2001) J Hydraul Res 39(6):567–573] hypothesis that secondary instabilities of the base flow may generate large-scale two-dimensional eddies, even in the absence of transverse gradients in the time-averaged flow properties.     

Characterization of storm wave asymmetries with functional data analysis

Functional data analysis (FDA) is a set of tools developed to perform statistical analysis on data having a functional form. In our case we consider the one-dimensional wave surface profiles registered during a North-Sea storm as functional data. The data is split into 20 min intervals within which an individual wave is defined as the profile between two consecutive downcrossings. After registration of these individual waves to the interval [0, 1], the mean wave profile for the entire 20 min interval is obtained along with the first two derivatives of this mean profile. We analyze the shape of these mean waves and their derivatives and show how they change as a function of the significant wave height, which is a measure of the severity of the sea for the corresponding time interval. We also look at the evolution of the energy, as represented by the phase diagram, as a function of significant wave height. The results show the asymmetry in vertical and horizontal scales for real data. Comparison with a Gaussian wave simulation model calculated from the actual wave spectra shows important differences in symmetry and shape of the average wave and seem to indicate that the greatest difference in the distribution of energy during the wave cycle lies in the second and fourth quarters of the wave period. FDA can be applied to derive information on the individual and average wave profiles and also provide an understanding of the variation in energy throughout the wave phase.     

Inorganic skeleton of the demosponge Halichondria panicea. Seasonality in spicule production in the Baltic Sea

Seasonality in spicule production of the demosponge Halichondria panicea (Pallas, 1766) was studied in the Western Baltic from July 1992 to July 1993. Spicule production increased in parallel with increasing temperatures in spring, being most pronounced during summer after release of larvae. Large proportions of up to 30% immature spicules occurred in the population during skeleton production peaks, and this resulted in reduction of mean spicule dimensions by ≃20%. Effects on spicule width were more obvious than on spicule length. In winter, when wave forces were largest, higher percentages of mineral skeleton were observed: 50% of dry weight compared to 40% in summer. Enhanced spicule densities in the tissues were not attained by marked spicule production but through a reduction in organic content by shrinking. In spring, growth of body tissue was faster than the build-up of skeleton, so that the relative amount of the latter was reduced. The indirect, growth-related developments of the skeleton were inversely proportional to temperature and directly to salinity, silicate concentrations and wave energy. Spicule size and the number of immature spicules were correlated with temperature and only slightly with salinity. No correlation with oxygen concentration or pH could be detected. Received: 3 April 1997 / Accepted: 8 August 1997     

An ecological risk assessment model for a pulsed contaminant emission into a wetland channel flow

As a continuation of the modelling on ecological degradation and hydraulic dispersion of pollutant emission into an idealized two-dimensional free-surface wetland flow (Zeng, L., Chen, G.Q., 2009b. Ecological degradation and hydraulic dispersion of contaminant in wetland. Ecol. Model., doi:10.1016/j.ecolmodel.2009.10.024), an ecological risk assessment model for the typical case of a pulsed contaminant emission into a realistic three-dimensional wetland channel flow is presented in this paper for the fate of cross-sectional mean concentration under environmental dispersion. An environmental dispersion model for the mean concentration is devised as an extension of Taylor’s classical analysis on dispersion in fluid flows. The velocity distribution and the environmental dispersivity in the fully developed steady flow through the wetland is found and illustrated with limiting cases covering various known solutions for the porous media flow between parallel plates, flow in a shallow wetland, sweeping flow in a densely vegetated wetland, and single phase flow in a channel. Obtained by Aris’s method of moments, the environmental dispersivity is shown characterized with multi-scale asymptotic time variations with stem dominated stage, transitional stage, and width-depth-stem dominated stage. Based on the solution for the evolution of contaminant cloud in the wetland channel flow, critical length and duration of the contaminant cloud with concentration beyond given environmental standard level are concretely illustrated for typical pollutant constituents in wastewater emission. Under the same emission intensity and environmental standard, the duration of contaminant cloud in the wetland channel is revealed shorter than that in a free surface wetland, due to the lateral effect.