Velocity profiles in a salt marsh canopy

Flow velocity profiles, measured in aSpartina anglica canopy in a laboratory flume, change with the location of measurement and plant stem density. The shear velocity above the canopy is larger than that within the canopy. The reduction ofu _* within the canopy will favor the deposition of cohesive sediment. The reducedu _* and flow turbulence within the canopy can enhance particle flocculation and settling velocity. The canopy exerts a strong influence on the concentration, settling velocity of the flocs, and deposition rate of the suspended sediment through effects on bed shear stress and turbulence of flow within the canopy.     

Observations of turbulence under weakly and highly stratified conditions in the Ariake Sea

Observations of turbulence, stratification, and mean current were made using a microstructure profiler and an acoustic Doppler current profiler (ADCP) during four cruises at a central location in the Ariake Sea, under weakly and strongly stratified conditions. Continuous measurements of the dissipation rate of turbulent kinetic energy (TKE), ε, were made. These revealed that frictional bed turbulence with quarterdiurnal variation in the bottom boundary layer (BBL) was one of the most energetic sources of vertical mixing in the sea. Thickness of the BBL was strongly confined by the stable stratification. We investigate a relationship between the BBL height h and the Ozmidov scale. We present a systematic argument that describes the vertical structure and characteristic scales of velocity and turbulence inside the frictional BBL, where the stratification persisted. Considerable deviation of observed vertical shear from the law of the wall indicated a modification of turbulent scales by the stratification. Shear stress calculated from the velocity data using vertical integration of the equation of motion was found to decrease approximately linearly with height. The TKE production rate P, estimated using the shear stress, was highly correlated with the dissipation rate. The buoyancy contribution to TKE balance in the BBL was quantified in terms of the flux Richardson number R _f as R _f?=?0.12.     

Hydrodynamics and associated sediment transport over coastal wetlands in Quanzhou Bay,China

Coastal salt marshes represent an important coastal wetland system.In order to understand the differences between boundary layer parameters of vegetated and unvegetated areas,as well as the mechanisms of sediment transport,several electromagnetic current meters (AEM HR,products of Alec Electronics Co.Ltd.) were deployed in coastal wetlands in Quanzhou Bay,China,to measure current velocity.During the low tide phase,the surficial sediment was collected at 10 m intervals.In situ measurements show that the current velocities on the bare flat were much higher than those in the Spartina alterniflora marsh.Current velocity also varied with distance from marsh edge and plant canopy height and diameter.Around 63% of the velocity profiles in the tidal creek can be described by a logarithmic equation.Over the bare flat and Spartina alterniflora marsh,a logarithmic profile almost occurs during the flood tide phase.Sediment analysis shows that mean grain size was 6.7 Φ along the marsh edge,and surface sediments were transported from bare flat to marsh;the tidal creeks may change this sediment transport pattern.The hydrodynamics at early flood tide and late ebb tide phases determined the net transport direction within the study area.     

Numerical study of turbulence and wave damping induced by vegetation canopies

Vegetation canopies control mean and turbulent flow structure as well as surface wave processes in coastal regions. A non-hydrostatic RANS model based on NHWAVE (Ma et al., 2012) is developed to study turbulent mixing, surface wave attenuation and nearshore circulation induced by vegetation. A nonlinear k − ϵ model accounting for vegetation-induced turbulence production is implemented to study turbulent flow within the vegetation field. The model is calibrated and validated using experimental data from vegetated open channel flow, as well as nonbreaking and breaking random wave propagation in vegetation fields. It is found that the drag-related coefficients in the k − ϵ model C_fk and C_fϵ can greatly affect turbulent flow structure, but seldom change the wave attenuation rate. The bulk drag coefficient C_D is the major parameter controlling surface wave damping by vegetation canopies. Using the empirical formula of Mendez and Losada (2004), the present model provides accurate predictions of vegetation-induced wave energy dissipation. Wave propagation through a finite patch of vegetation in the surf zone is investigated as well. It is found that the presence of a finite patch of vegetation may generate strong pressure-driven nearshore currents, with an onshore mean flow in the unvegetated zone and an offshore return flow in the vegetated zone.     

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

Validating a turbulence closure against estuarine microstructure measurements

A simple k–ε turbulence closure is introduced which has no stability functions but instead a Richardson number-dependent turbulent Prandtl number. Its free parameters are determined in a comparison with microstructure observations from a stratified and sheared tidal estuary and laboratory measurements. The closure is able to simulate observed turbulent dissipation rates (ε) and turbulent length scales (l_th) in regions of strong mean shear and small gradient Richardson number (R_g) to within factors of 2–3. It fails in regions of small shear and large R_g, presumably because of the dominance of internal wave-driven mixing. Additional simulations with a k–ε closure with stability functions taken from Canuto et al. [Canuto, V.M., Howard, A., Cheng, Y., Dubovikov, M.S., 2001. Ocean turbulence I: one-point closure model. Momentum and heat vertical diffusivities. J. Phys. Oceanogr. 31, 1413–1426] and with the closure of Baumert and Peters [Baumert, H., Peters, H., 2004. Turbulence closure, steady state, and collapse into waves. J. Phys. Oceanogr. 34, 505–512] show poor performance. Establishing a valid 1:1 comparison of simulated and observed ε and l_th requires nudging the model velocity and density toward observed values because free model integrations quickly diverge from the observations. Steady state gradient Richardson numbers are constrained to a range of 0.18–0.25, while flux Richardson numbers are constrained to the range of 0.1–0.22. The closure output is rather insensitive to such parameter variations. The simulations are sensitive, however, to the treatment of the observed velocity and density used to nudge the model. Good closure performance requires averaging the measured tidal flow over about an hour, a time scale for which conventional numerical models of estuarine circulations should be able to match observed shears. In the closure simulations the TKE balance stays close to a production–dissipation balance. The time rate of change and vertical diffusion of TKE are small, of the same order of magnitude, and vary in magnitude relative to each other systematically across the water column.     

生态空心块体水沙动力效应研究

生态空心块体是一种兼有调控水流与生态修复功能的工程结构物,在海堤、丁坝等近海工程中被广泛应用.通过室内水槽试验,设计多种来水来沙工况,研究生态空心块体与柔性植被共同影响下的水沙动力特性,重点探究生态空心块体内的水沙动力效应.结果 表明:生态空心块体组成的坝体结构缓流效果显著,坝体内外侧水流流速分布差异明显,坝体内部的水...     

黄河三角洲互花米草整治对潮滩湿地防灾功效和沉积动力的影响

互花米草作为黄河三角洲外来入侵物种, 2010年开始爆发式蔓延, 2020年面积最大时可超6 000 hm²,对潮间带生物多样性造成了严重影响。自2020年起,黄河三角洲自然保护区开始大规模的互花米草整治与盐地碱蓬生态修复工作。互花米草具有较强的消浪、缓流、促淤、固滩能力,短期内、大范围的植被去除,可能导致波流和沉积动力环境的突变。基于现场观测数据和考虑植被作用的浪、流、沙耦合数值模型,研究分析了互花米草整治和盐地碱蓬修复对黄河三角洲湿地动力过程的影响,重点关注风暴作用下黄河三角洲潮滩湿地的防灾减灾能力和沉积效应变化。研究结果表明,互花米草可有效抵御极端风暴潮,最大风暴潮位衰减率可达15 cm/km。互花米草整治将导致湿地防灾减灾能力显著下降,堤前风暴潮峰值可抬升10～30 cm,同时潮间带的泥沙捕集能力大幅下降。受限于植被密度、高度,盐地碱蓬完全修复后湿地衰减风暴潮、捕集泥沙的能力显著低于原互花米草湿地。研究建议,互花米草整治工程应考虑短期内植被去除导致的近岸风暴潮抬升、侵蚀加剧等潜在风险。     

A Linked Physical and Biological Framework to Assess Biogeochemical Dynamics in a Shallow Estuarine Ecosystem

The littoral zone of Chesapeake Bay contains a mosaic of shallow vegetated and nonvegetated habitats with biotic components that are sensitive to changes in biological and physical driving factors. Static and dynamic modelling frameworks provide an integrative way to study complex hydrodynamic and biogeochemical processes in linked estuarine habitats. In this study we describe a spatial simulation model developed and calibrated relative to a specific littoral zone, estuarine ecosystem. The model consisted of four distinct habitats that contained phytoplankton, sediment microalgae, Zostera marina (eelgrass), and Spartina alterniflora. There was tidal exchange of phytoplankton, particulate and dissolved organic carbon and dissolved inorganic nitrogen between the littoral zone ecosystem and the offshore channel. Physical exchange and biogeochemical transformations within the habitats determined water column concentrations in each habitat. Predicted subtidal water column concentrations and Z. marina and S. alterniflora biomass were within the variability of validation data and the predicted annual rates of net primary production were similar to measured rates. Phytoplankton accounted for 17%, sediment microalgae 46%, the Z. marina community 24% and S. alterniflora 13% of the annual littoral zone primary production. The linked habitat model provided insights into producer, habitat and ecosystem carbon and nitrogen properties that might not have been evident with stand-alone models. Although it was an intra-ecosystem sink for particulate carbon, the seagrass habitat was a DOC source and responsible for over 30% of the littoral zone carbon and nitrogen primary production. The model predicted that the Goodwin Islands littoral zone was a sink of channel derived POC, but a source of DOC to the surrounding estuary. The framework created in this study of estuarine ecosystem dynamics is applicable to many different aquatic systems over a range of spatial and temporal scales.     

Large-scale turbulence under a solitary wave

The structure of large-scale turbulence under a broken solitary wave on a 1 in 50 plane slope was studied. Three-component velocity measurements were taken at different heights above a smooth bed in the middle surf zone using an acoustic Doppler velocimeter. The measured data showed that turbulent velocity components were well correlated in the middle part of the water column. The velocity correlations could be produced by an oblique vortex similar to the obliquely descending eddy observed previously by other investigators. The vertical distributions of the relative values of the components of the Reynolds stress tensor showed that the structure of turbulence evolved continuously between the free surface and the bottom. The evolution was related to transition from two-dimensional to three-dimensional flow structures and the effect of the solid bottom on flow structures. Time histories of measured turbulent kinetic energy and turbulence stresses showed episodic turbulent events near the free surface but more sporadic turbulence in the lower layer. Large or intense turbulent events were found to have short duration and time lag relative to the wave crest point. These events also maintained good correlations between the turbulence velocity components close to the bottom.Instantaneous turbulent velocity fields were measured near the bottom at the same cross-shore location by using a stereoscopic particle image velocimetry system. These measurements showed that the near-bed flow field was characterized by large-scale, coherent flow structures that were the sources of most of the turbulent kinetic energy and turbulence stresses. The types of organized flow structures observed included vortices and downbursts of turbulence descending directly from above, lateral spreading of turbulent fluid along the bed, and formation of vortices in shear layers between fluid streams. A common feature of the organized flow structures near the bed was the large turbulence velocities in the longitudinal and transverse directions, which reflected the influence of a solid bottom on the breaking-wave-generated turbulence arriving at the bed.     

长江口湍流剖面的观测与分析

小尺度湍流过程对河口物质输运与能量交换至关重要。受传统观测方法的限制, 河口浅水区域的剖面观测资料至今较为匮乏, 进而限制了湍流过程的研究。为此, 采用新型5波束声学多普勒流速剖面仪(Nortek Signature 1 000 kHz AD2CP)在长江口开展湍流剖面观测, 获取高频、低噪的高质量湍流剖面数据, 并与声学多普勒点式流速仪(acoustic doppler velocimeters, ADV)同步观测的数据进行对比。结果表明, 通过AD2CP与ADV获得的近底部边界层摩阻流速u_*、拖曳系数C_d、雷诺应力S_R等特征参数基本一致, 底摩擦与波浪能量为河口区域湍动能的主要输入源。湍流垂向结构存在显著的非局地平衡, 即温盐等斜压作用引起的浮力通量、对流项以及强波浪作用影响的脉动压力做功、黏性输运等因素可能为长江口湍流非局地平衡的主因。     

Modeling of the Turbulence in the Water Column under Breaking Wind Waves

Past studies have shown that there is a wave-enhanced, near-surface mixed-layer in which the dissipation rate is greater than that derived from the “law of the wall”. In this study, turbulence in water columns under wind breaking waves is investigated numerically and analytically. Improved estimations of dissipation rate are parameterized as surface source of turbulent kinetic energy (TKE) for a more accurate modelling of vertical profile of velocity and TKE in the water column. The simulation results have been compared with the experimental results obtained by Cheung and Street (1988) and Kitaigorodskii et al. (1983), with good agreement. The results show that the numerical full model can well simulate the near-surface wave-enhanced layer and suggest that the vertical diffusive coefficients are highly empirical and related to the TKE diffusion, the shear production and the dissipation. Analytical solutions of TKE are also derived for near surface layer and in deep water respectively. Near the surface layer, the dissipation rate is assumed to be balanced by the TKE diffusion to obtain the analytical solution; however, the balance between the dissipation and the shear production is applied at the deep layer. The analytical results in various layers are compared with that of the full numerical model, which confirms that the wave-enhanced layer near the surface is a diffusion-dominated region. The influence of the wave energy factor is also examined, which increases the surface TKE flux with the wave development. Under this region, the water behavior transits to satisfy the classic law of the wall. Below the transition depth, the shear production dominantly balances the dissipation. This revised version was published online in August 2006 with corrections to the Cover Date.     

A numerical study of swash flows generated by bores

The dynamic processes of bore propagation over a uniform slope are studied numerically using a 2-D Reynolds Averaged Navier–Stokes (RANS) solver, coupled to a non-linear k − ε turbulence closure and a volume of fluid (VOF) method. The dam-break mechanism is used to generate bores in a constant depth region. Present numerical results for the ensemble-averaged flow field are compared with existing experimental data as well as theoretical and numerical results based on non-linear shallow water (NSW) equations. Reasonable agreement between the present numerical solutions and experimental data is observed. Using the numerical results, small-scale bore behaviors and flow features, such as the bore collapse process near the still-water shoreline, the ‘mini-collapse’ during the runup phase and the ‘back-wash bore’ in the down-rush phase, are described. In the case of a strong bore, the evolution of the averaged turbulence kinetic energy (TKE) over the swash zone consists of two phases: in the region near the still-water shoreline, the production and the dissipation of TKE are roughly in balance; in the region farther landwards of the still-water shoreline, the TKE decay rate is very close to that of homogeneous grid turbulence. On the other hand, in the case of a weak bore, the bore collapse generated turbulence is confined near the bottom boundary layer and the TKE decays at a much slower rate.     

Laboratory measurements of large-scale near-bed turbulent flow structures under spilling regular waves

The instantaneous turbulent velocity field created by the breaking of spilling regular waves on a plane slope was measured in a plane running parallel to the slope using particle image velocimetry. The measurement plane was located at a height of about 1 mm above the bed. The measurement area encompassed the region where the large eddies generated at incipient wave breaking impinged on the bottom inside the surf zone. A total of 30 trials were conducted under identical experimental conditions. In each trial, six consecutive wave cycles were recorded. The measured velocity fields were separated into a mean flow and a turbulence component by ensemble averaging. The instantaneous turbulent velocity fields were analyzed to determine the occurrence frequency, location, geometry and evolution of the large eddies, and their contributions to instantaneous shear stresses, turbulent kinetic energy and turbulence energy fluxes. The motion of single glass spheres along the bed was also investigated. The two-phase flow measurements showed that the velocity and displacement of large solid particles on a smooth bed were significantly affected by the magnitude and direction of turbulence velocities. Overall, this study has examined the kinematic and dynamic properties of large eddies impinging on the bed and the interaction of these large-scale turbulent flow structures with the mean flow. The study has also highlighted the important role of large eddies in sediment transport.     

Spatial and temporal variations in sediment grain size in tidal wetlands, Yangtze Delta: On the role of physical and biotic controls

To examine the spatial and temporal variability of sediment grain size in exposed tidal wetlands with ample sediment supply, we sampled sediments and measured hydrodynamics, accretion/erosion rates, and vegetation characteristics in the Yangtze Delta. Sediment grain size exhibited a landward/upward decreasing trend. This trend is mainly attributed to attenuation of hydrodynamics. A 630-day series of daily surface sediment sampling at a fixed site on an unvegetated intertidal flat revealed significant seasonal and storm-cyclic changes in grain size. This temporal variability was related to alternating accretion/erosion events, with erosion associated with coarser grain size. Such temporal dynamics were not present in vegetation, where sediment remained fine grained throughout the year. In the marsh, vegetation cover enables the trapping of fine-grained sediments in the following ways: (a) adherence of suspended sediments onto plants; (b) deposition of suspended sediments stimulated by attenuation of hydrodynamics through plant obstruction; and (c) prevention of resuspension of fine-grained deposits due to the protection of the plant canopy. The influence of vegetation on sediment grain size was clearly seen when comparing sediment trapped by different vegetation types and seasonal patterns of trapped sediment on different vegetation canopy densities. The relatively high plant biomass of the recently introduced Spartina alterniflora enhanced the trapping effect, whereas plant degradation due to buffalo grazing reduced the trapping effect. We conclude that for exposed tidal wetlands with ample sediment supply such as the Yangtze Delta, the spatial and temporal variability of sediment grain size is governed predominantly by physical controls on the unvegetated flat and predominantly by biophysical interaction of hydrodynamics and vegetation in the salt marsh, rather than by sediment supply.     

A pair of Langmuir cells in a laboratory tank (I) wind-only experiment

Three velocity components of subsurface flow, observed in a rectangular tank under the action of a constant wind speed, are measured systematically at mesh points distributed uniformly over a vertical cross-section of the tank. Measurements are carried out for two cases: 1) reference wind speedU _r =7.5 m/s and fetchF=10 m; and 2)U _r =10 m/s andF=25 m. A pair of Langmuir cells is observed for both cases; downwelling zones are found along both of the sidewalls and an upwelling zone in the centre of the tank. Near the water surface, the vertical momentum flux is dominated by the Reynolds stress resulting from small-scale turbulence, while over the entire cross-section except near the surface, the Reynolds stress due to the Langmuir cells dominates the vertical momentum flux. As the result of the occurrence of this Langmuir cells, the vertical momentum flux, which consists of both mean advection and small-scale turbulence, is markedly inhomogeneous in the spanwise direction; for example, the largest vertical flux of the order of the wind stress is observed in the downwelling zone near one sidewall, while at the centre of the tank, the vertical momentum flux occupies only 30% of the wind stress. This indicates that a pair of Langmuir cells plays more important role than small-scale turbulence in the mixing process in a greater part of the wind-wave tank.Address after April 1, 1992: Department of Civil Engineering, Hiroshima Institute of Technology, Miyake 2-1-1, Saeki-ku, Hiroshima 731-51, Japan.     

Numerical simulation of solitary and randomwave propagation through vegetation based on VOFmethod

A vertical two-dimensional numerical model has been applied to solving the Reynolds Averaged Navier- Stokes （RANS} equations in the simulation of current and wave propagation through vegetated and non- vegetated waters. The k-e model is used for turbulence closure of RANS equations. The effect of vegeta- tion is simulated by adding the drag force of vegetation in the flow momentum equations and turbulence model. To solve the modified N-S equations, the finite difference method is used with the staggered grid system to solver equations. The Youngs＇ fractional volume of fluid （VOF） is applied tracking the free sur- face with second-order accuracy. The model has been tested by simulating dam break wave, pure current with vegetation, solitary wave runup on vegetated and non-vegetated channel, regular and random waves over a vegetated field. The model reasonably well reproduces these experimental observations, the model- ing approach presented herein should be useful in simulating nearshore processes in coastal domains with vegetation effects.     

Direct numerical simulations of flow over two-dimensional and three-dimensional ripples and implication to sediment transport: Steady flow

A well resolved and highly accurate direct numerical simulation (DNS) solver has been developed to understand the implication of hydrodynamics to sediment transport. In the first part of the study we focus on steady flow over two-dimensional and three-dimensional ripples at two Reynolds numbers Re_τ = 180 and 400 (defined by channel half-height and wall-friction velocity) in a channel geometry. The DNS scheme is based on a fourth-order vertical velocity and second-order vertical vorticity formulation, which resolves the difficulties in pressure boundary condition encountered when solving the Navier–Stokes equations. The complex boundary introduced due to the ripples has been imposed in the Cartesian domain using an elegant immersed boundary method. Detailed hydrodynamic analysis has revealed turbulence statistics (in particular, the higher order) and henceforth, the flow structures are sensitive — whether the ripples are two-dimensional or three-dimensional. The importance of fluctuating component of the bottom stress in addition to its mean component; and its significance to sediment transport and ripple migration speed have been investigated.     

The structure and dynamics of the Bottom Boundary Layer in shallow sea areas without tidal influence: an experimental approach

This report describes extensive investigations of the near bottom layer of the Western Baltic (Mecklenburg Bight, Darss Sill and Arkona Basin) which were conducted over a 5 year period to determine the typical structure, vertical thickness, vertical turbulence structure, and spatial and temporal variability of this water mass with regard to the area's particular hydrographic conditions. Series of vertical profiles were obtained using the microstructure profiler MSS86, which is capable of measuring high resolution profiles of temperature, conductivity, current shear, light attenuation and pressure down to the seafloor. The near bottom current structure was simultaneously measured with conventional current metres at fixed depths. A typical vertical density structure of the near bottom layer was found. At all investigation sites the Bottom Boundary Layer was separated from the overlying water mass by a well pronounced thermohaline pycnocline. A homogeneous water layer was situated above the bottom with a mean thickness of 2.2 m and typical variation between 0.5 and 3.5 m. The thickness of both the homogeneous layer and of the near bottom layer vary considerably. It is suggested that horizontal advection is responsible for these fluctuations in thickness. The variation in thickness of the Homogeneous Layer is independent of the local mean current velocity, wind speed and energy dissipation rate. Over periods of about 2 days the thickness of the Homogeneous Layer is determined by the average wind speed. The Bottom Boundary Layer shows its own characteristic dynamic, which is largely decoupled from that of the remaining water body. A logarithmic layer was generally not resolved by the current measurements. From dissipation rate measurements, the wall layer was determined to be 0.9 m thick. There was no significant correlation between the dissipation rate and the local wind speed, or between the dissipation rate and local mean current u₁₀₀. This means that any simple parameterisation relating u₁₀₀ or friction velocity to the locally produced turbulence and consequently to the resuspension of sediment is probably not applicable to shallow sea areas with properties like the Western Baltic. The investigation of sediment concentration in the BBL illustrates the importance of local effects combined with advection. The sediment stratified layer covers only the bottom most 50 cm.     

Modeling of the eddy viscosity by breaking waves

Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing,heat and gases exchange across air-sea interface.The energy loss from wave breaking and the bubble plume penetration depth are estimated.As a consequence,the vertical distribution of the turbulent kinetic energy(TKE),the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided.It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface.High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than υwl(=κuwz),the value predicted for the wall layer scaling close to the surface,where uw is the friction velocity in water,κ with 0.4 is the von Kármán constant,and z is the water depth,and the strength of the eddy viscosity depends both on wind speed and sea state,and decays rapidly through the depth.This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer,well above the classical values expected from the similarity theory.Deeper down,however,the effects of wave breaking on the vertical mixing become less important.