渤海垂直湍流混合强度季节变化的数值模拟

渤海为极浅陆架海 ,其中湍流耗散作用显著。将三维斜压陆架海模式 HAMSOM应用于渤海 ,以渤海周边台站每天 4次的常规气象资料作为风和热驱动 ,渤海海峡开边界以 5个主要分潮调和常数计算水位强迫 ,计算了渤海 1982年水文要素和流场变化 ,并用模式以湍的局地平衡理论封闭计算出垂直湍流粘性的时空分布。结果表明 :渤海湍流混合冬强夏弱 ,变化幅度较大 ( 10～ 2 0 0 cm2 / s) ,这是风搅拌和潮混合的湍流输入在密度层化调整下的结果 ;风的作用在冬季强于潮的作用 ,而底层则由潮混合控制呈现半月周期 ;渤海湍粘性系数的空间分布十分复杂 ,这是在渤海地形和岸形轮廓限制下 ,由一定大气条件驱动的流场和密度场导致的湍流混合强度不同所致     

Physical processes associated with high primary production in Saanich Inlet, British Columbia

Saanich Inlet, British Columbia, has long been known for the presence, in most years, of anoxic bottom water. One factor contributing to this anoxia is a high level of primary production, which occurs as a major spring bloom followed by sporadic ‘mini-blooms’ throughout the summer and early fall. The process(es) by which new production is refueled after nutrient exhaustion caused by the spring bloom are not well known, since Saanich is an inverse estuary and vertical mixing driven by winds and tides is low. This study presents new observational evidence that strongly suggests that the dominant mechanism of nutrient resupply during the summer months is intermittent advective exchange, driven by pressure gradients set up by strong tidal mixing in passages outside Saanich Inlet itself. A simple box model is formulated to illustrate this mechanism. When driven by annual freshwater forcing and deepwater renewal functions characteristic of the region and measured tides for 1975, the model predicts resupply of nitrate during most of the periods observed in 1975 observations (Deep-Sea Res. 24 (1977) 775). This ‘action-at-a distance’ nutrient resupply mechanism, involving strong but localized turbulent mixing and subsequent distribution of the products of mixing over large-horizontal distances by pressure-gradient-driven flow, is likely important in other coastal regions where the estuarine circulation is weak.     

Modélisation sur vingt ans (1976–1995) de la production phytoplanctonique en baie de Seine

This paper presents an ecological modelling of the Bay of Seine (Eastern Channel) over the last twenty years, chosen as a typical case of eutrophication in a river plume. In the physical sub-model, the Bay is divided into 42 boxes and water fluxes between them are calculated automatically using Ifremer's “Elise” software. A two-layer, vertical thermohaline model is then linked with the horizontal circulation scheme in order to take vertical stratification into account. The biological submodel deals with two chemical elements, nitrogen and silicon, and splits phytoplankton into diatoms and flagellates. Results from this ecological model point out the spatial concordance of highest phytoplanktonic concentrations with the river plume spreading in the bay. Contrary to diatoms, flagellate production appears to be mainly confined to the eastern bay, due to the vertical haline stratification in front of the river mouth. As far as the whole bay is concerned, the interannual fluctuations of diatom production are related to the level of spring insolation, whereas silicon inputs regulate diatom production in the river plume. The flagellate summer production in the plume is enhanced by high water temperature and high N/Si ratios, which appear during dry years with low discharge regimes. Finally, interannual increase of flagellate production could be related to gradual increase of nitrogen loadings, contrasting with silicon loadings, which remained stable for twenty years.     

Simulations of Annual Cycle of Phytoplankton Production and the Utilization of Nitrogen in the Yellow Sea

A nutrient dynamic model coupled with a 3D physical model has been developed to study the annual cycle of phytoplankton production in the Yellow Sea. The biological model involves interactions between inorganic nitrogen (nitrate and ammonium), phosphate and phytoplankton biomass. The model successfully reproduces the main features of phytoplankton-nutrient variation and dynamics of production. 1. The well-mixed coastal water is characterized by high primary production, as well as high new production. 2. In summer, the convergence of tidal front is an important hydrodynamic process, which contributes to high biomass at frontal areas. 3. The evolution of phytoplankton blooms and thermocline in the central region demonstrate that mixing is a dominant factor to the production in the Yellow Sea. In this simulation, nitrate- and ammonium-based productions are estimated regionally and temporally. The northern Yellow Sea is one of the highly ranked regions in the Yellow Sea for the capability of fixing carbon and nitrogen. The annual averaged f-ratio of 0.37 indicates that regenerated production prevails over the Yellow Sea. The result also shows that phosphate is the major nutrient, limiting phytoplankton growth throughout the year and it can be an indicator to predict the bloom magnitude. Finally, the relative roles of external nutrient sources have been evaluated, and benthic fluxes might play a significant role in compensating 54.6% of new nitrogen for new production consumption.     

江苏近海潮流湍流切应力研究

潮流的流速分布和湍流切应力的分布密切相关。为了研究潮流湍流切应力的特征和变化规律,本文从潮流运动方程推导出随潮流呈周期性变化的湍流切应力随深度的分布表达式。从江苏近海现场实测流速剖面中计算得到一个潮周期内的湍流切应力的变化过程。通过将实测值与本文理论解的比较发现,计算值能很好地反映出实测值的变化特征,尤其是能正确地反映潮流在加速和减速过程中,湍流切应力偏离线性分布呈现上凹和下凹的现象,以及只有在加减速转换的时候切应力才呈现线性分布的特点。最后,通过分析湍流切应力振幅及相位沿水深的变化情况,认为反映潮流周期、涡粘性和水深之间关系的参数Ri是决定湍流切应力弯曲程度及分布形态的重要参数。对于近海潮流,水深是影响湍流切应力偏离线性程度的主要因素。     

Quantifying cross-shelf and vertical nutrient flux in the Coastal Gulf of Alaska with a spatially nested,coupled biophysical model

The Coastal Gulf of Alaska (CGOA) is productive, with large populations of fish, seabirds, and marine mammals; yet it is subject to downwelling-favorable coastal winds. Downwelling regions in other parts of the world are typically much less productive than their upwelling counterparts. Alternate sources of nutrients to feed primary production in the topographically complex CGOA are poorly known and difficult to quantify. Here we diagnose the output from a spatially nested, coupled hydrodynamic and lower trophic level model of the CGOA, to quantify both horizontal and vertical nutrient fluxes into the euphotic zone. Our nested model includes both nitrogen and iron limitation of phytoplankton production, and is driven by a fine-scale atmospheric model that resolves the effects of local orography on the coastal winds. Results indicate significant “rivers” of cross-shelf nitrogen flux due to horizontal advection, as well as “fountains” of vertical transport over shallow banks due to tidal mixing. Using these results, we constructed a provisional budget of nutrient transport among subregions of the CGOA. Contrary to expectations, this budget reveals substantial upwelling of nutrients over major portions of the shelf, driven by local wind-stress curl. These effects are large enough to overwhelm the smaller downwelling flux at the coast throughout the growing season. Vertical mixing by winds and tides, and horizontal flux from the deep basin, are other substantial contributors to nutrients above the 15-m horizon. These findings help to explain the productivity of this coastal ecosystem.     

An ecosystem model with iron limitation of primary production in the equatorial Pacific at 140°W

We construct a one-dimensional ecosystem model (nitrate, ammonium, phytoplankton, zooplnakton and detritus) with simple physics and biology in order to focus on the structural relations and intrinsic properties of the food web that characterizes the biological regime in the central equatorial Pacific at 140°W. When possible, data collected during the EgPac and other cruises were used to calibrate model parameters for two simulations that differ in the limiting nutrient, i.e. nitrogen or iron. Both simulations show annual results in good agreement with the data, but phytoplankton biomass and primary production show a more pronounced annual variability when iron is used as the limiting nutrient. This more realistically reproduces the variability of biological production and illustrates the greater coupling between vertical physical processes and biological production when the limiting nutrient is iron rather than nitrogen. The iron simulation also illustrates how iron supply controls primary production variability, how grazingbalances primary production and controls phytoplankton biomass, and how both iron supply and grazingcontrol primary production. These results suggest that it is not possible to capture primary production variability in the central equatorial Pacific with biological models using nitrogen as the limiting nutrient. Other indirect results of this modeling study were: (1) partitioning of export production between dissolved and particulate matter is almost equal, suggesting that the importance of DOC export may have been previously overestimated; (2) lateral export of live biomass has to be taken into account in order to balance the nitrogen budget on the equator at 140°W; and (3) preferential uptake of ammonium (i.e. nitrate uptake inhibition by ammonium) associated with high regeneration of nitrogen (low f ratio as a consequence of the food web structure imposed by iron limitation) largely accounts for the surface build-up of upwelled nitrate.     

Mathematical modeling of turbulent jets of deep-water sewage discharge into coastal basins

The basic properties of the dynamic model of a turbulent jet formed by a deep-water sewage discharge into the stratified environment of coastal regions are considered. The model developed was used to estimate the parameters of a floating-up jet of deep wastewater discharge from Sand Island into the basin of Mamala Bay (Hawaii) depending on the season and discharge operation mode. The estimates of the float-up depths of the jet and the initial dilution of the jet were estimated on the basis of model calculations using experimental data on the vertical profiles of the water temperature and salinity under the actual conditions of stratification in the study region. It is shown that the further propagation of the wastewater jet depends on tidal events and internal waves generated by tides. The appearance of turbulent jets at the sea surface was recorded. The model estimates of the parameters of the wastewater discharge were compared with the results of experimental measurements. Good agreement was found, which indicates that the physical mechanisms of the propagation of turbulent jets in a stratified medium are adequately described by the model.     

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

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

Combined tidal velocity and duration asymmetries as a determinant of water transport and residual flow in Paranaguá Bay estuary

Net transports of water, salt and suspended particulate matter (SPM) for a cross-section in front of the Paranaguá Harbour (Paranaguá Bay, Brazil) are presented for eight distinct tidal cycles. Data include measurements over single spring and neap tidal cycles, during both wet and dry seasons. The main drive forces of circulation and SPM dynamics are identified. Advective transport dominated under moderate to high vertical salinity stratification and weak currents, while mixing processes dominated under well-mixed conditions generated by high currents and low freshwater input. Under partially mixed conditions, both advective and mixing processes were important. The tide-induced residual circulation dictated the magnitude and direction of residual currents and net transports of water and salt, but not of SPM transport. The SPM dynamics was intrinsically related to cyclical processes of erosion, resuspension and deposition driven by tidal currents. The turbulent mixing intensity conditioned the vertical mixing of SPM. Resuspension and vertical mixing were conspicuous in spring cycles, while the horizontal advection preponderated in the neaps. Lags between maximum currents and SPM peak concentrations occurred, with more pronounced hysteresis during ebb periods.     

湍流局地平衡假设的新推论--齐次湍流动能输运方程封闭模型与应用(Ⅱ)

在湍流局地平衡假设下 ,建立了齐次湍能输运方程封闭模型 (HKE) ,并在平板边界层的两种经典流动中加以检验 ,给出 HKE封闭下的流速、湍流动能和湍流混合系数剖面的形式解。结果表明 ,HKE可以避免在流速剪切为零时的无湍流混合问题 ,其解与 L aufer湍流实验吻合 ,因而HKE模型比混合长理论有更合理的内涵。文中还给出 HKE封闭的浅海动力学模型 ,以湍应力和水位梯度力的平衡为运动的基本受力平衡 ,进行了模型的量阶分析和运动分析 :当阻尼频率和运动频率同量阶时 ,惯性运动不可忽略 ;在潮振荡占优的浅海中 ,对流非线性相对于惯性运动为小量 ;当阻尼频率足够大时 ,科氏力项相对于湍应力也可能为小量。     

Physical and biological modeling in the Gulf Stream region: Part II. Physical and biological processes

Mesoscale physical and biological processes are examined at the Gulf Stream front by means of a 4-D simulation including physical and biological data assimilation. The data assimilated are from Leg 1 of the Fall BIOSYNOP cruise, 21 Sept.–8 Oct. 1988, and GULFCAST data for the same period. Focus is on the vertical velocities at the front, the vertical and horizontal transports of nutrients and plankton, and the impact of these transports on phytoplankton biomass, production and organic particle export. It was found that while jet meandering enhances new production at the front, primary production and phytoplankton concentration at the front are not significantly enhanced over those of Slope water. Winds during this period also have little impact on productivity at the front, due to their high temporal variability. Ring–stream interactions, however, significantly increase the net vertical and meridional transports of nutrients and plankton and can lead to phytoplankton patchiness at the front. This emphasizes the importance of submesoscale events between interacting mesoscale physical features in the transport of nutrients and plankton, and in explaining the observations. The enhanced phytoplankton concentrations observed during BIOSYNOP are found to be primarily due to advection (convergence) rather than in situ biological growth.     

三维海洋溢油预测模型的建立(英文)

提出一个海洋溢油三维物理归宿和输运的动力学综合模型。该模型包含了一系列的数学公式来描述对流、湍扩散、表面扩展、铅直扩散、乳化和蒸发过程。每一公式的建立是独立的并且与相关过程、环境和其它参数相联系。该模型需要输入流场作为输运的媒介 ,这可从感兴趣区域的三维潮和风驱动的流体动力学模型获得。模型用来预测和后报溢油在海洋环境中的归宿和输移 ,可为溢油应急反应和环境影响评价服务。     

A dispersion equation and dispersion coefficient for material transport in Inland Seas

To analyse material transport in inland seas, a horizontal two-dimensional dispersion equation is derived, and the dispersion coefficient due to the combined effect of vertical turbulent mixing and vertical shear of both a steady current and a tidal current is studied. In the present study, the assumption that velocity is uniform in horizontal planes is not necessary, and velocity has a free vertical profile; thus the dispersion coefficient formulated is general, and is represented by a tensor of the second order. The properties of the dispersion coefficient in the horizontal two-dimensional dispersion model are also investigated, and it is shown that the time-averaged dispersion coefficient due to the tidal current over a tidal period is approximately half that due to the steady current, if the velocity amplitude and the vertical profile of the tidal current are the same as those of the steady current (a similar result was presented byBowden (1965) for horizontal one-dimensional models). Finally, the dispersion coefficient in Hiuchi-Nada (Hiuchi Sound) in the central part of the Seto Inland Sea is evaluated by using the model. The values of the dispersion coefficient in that region range from 10³ cm² s^–1 to 10⁵ cm² s^–1 when vertical turbulent diffusivity is taken to be 50 cm² s^–1.     

Similitude of the hydraulic model experiment for tidal mixing

The treatment of diffusion due to the tidal current in the near-shore ocean, and the similitude in the hydraulic model experiment are studied.In broad and shallow tidal bays and in the coastal seawaters near irregular boundaries, horizontal eddy-currents induced geometrically and turbulence caused by their cascading have predominant effect on dispersion of river and waste waters. These turbulent diffusion processes are similarly reproduced in the Froude models of turbulent resume, by adding the similitude for the self-similar structure of the spectral density of turbulence, or the eddy diffusivity. The similitude means to take the scale ratios of the time and the vertical length as the two-thirds power of the scale ratio of the horizontal length.Similitudes are also derived for the system of the gravitational circulation, the stratification and the salt-mass transport in partially and well-mixed estuaries. Generally, the vertical eddy diffusivity must be exaggerated by a half power of the model distortion externally by some methods of agitation. When the tidal bay is broad and very shallow, the Rayleigh number and the Hansen number are small, and the effect of density current and stratification on the flushing is small. Instead the effect of local eddies, geometrically induced tidal residual circulations become predominant. In this special case, there is no need to satisfy the similitude for density difference and vertical shear effects on dispersion.     

Three-dimensional Layer-integrated Modelling of Estuarine Flows with Flooding and Drying

Details are given of the refinement and application of a thee-dimensional (3-D) layer-integrated numerical model of tidal circulation, with the aim of simulating severe tidal conditions for practical engineering applications. The mode splitting strategy has been used in the model. A set of depth-integrated 2-D equations are first solved to give the pressure gradient, and the layer-integrated 3-D equations are then solved to obtain the vertical distributions of the flow velocities. Attention has been given to maintaining consistency of the physical quantities derived from the 2-D and 3-D equations. A TWO=layer mixing length turbulence model for the vertical shear stress distribution has been included in the model. Emphasis has been focused on applying the model to a real estuary, which is geometrically complicated and has large tidal ranges giving rise to extensive flooding and drying. The model has been applied to three examples, including: wind-driven flow in a rectangular lake, tidal circulation in a model rectangular harbour, and tidal circulation in a large estuary. Favourable results have been obtained for both the simple and complex flow beds.     

Combined tidal and wind driven flows and bedload transport over a flat bottom

The combined tidal and wind driven flow and resulting sediment transport in the ocean over a flat bottom at intermediate water depth has been investigated, using a simple one dimensional two-equation turbulence closure model. This model has been verified against field measurements of a tidal flow in the Celtic Sea. The tidal velocity ellipses and the time series of the horizontal velocity components at given elevations above the bottom are well predicted through the water column although there are some deviations between the predicted and measured velocities near the bottom due to the uncertainty of the bottom roughness. For the combined tidal and wind driven flows the velocity profiles, turbulent kinetic energy profiles and surface particle trajectories are predicted for weak and strong winds. Furthermore, the bottom shear stress and the resulting bedload transport have been predicted; the parts of the particle trajectories in the close vicinity of the bottom where the bedload transport exists are displayed. Finally, the direction and magnitude of the surface drift, the depth-averaged mean velocity and the mean bedload transport are given, and the effect of the bottom roughness on the sea surface drift is investigated.     

Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

Turbulent mixing in a stratified estuarine tidal channel: Hikapu Reach,Pelorus Sound,New Zealand

Channel constrictions within an estuary can influence overall estuary-sea exchange of salt or suspended/dissolved material. The exchange is modulated by turbulent mixing through its effect on density stratification. Here we quantify turbulent mixing in Hikapu Reach, an estuarine channel in the Marlborough Sounds, New Zealand. The focus is on a period of relatively low freshwater input but where density stratification still persists throughout the tidal cycle, although the strength of stratification and its vertical structure vary substantially. The density stratification increases through the ebb tide, and decreases through the flood tide. During the spring tides observed here, ebb tidal flow speeds reached 0.7?m?s^?1 and the buoyancy frequency squared was in the range 10^?5 to 10^?3?s^?2. Turbulence parameters were estimated using both shear microstructure and velocimeter-derived inertial dissipation which compared favourably. The rate of dissipation of turbulent kinetic energy reached 1?×?10^?6?m²?s^?3 late in the ebb tide, and estimates of the gradient Richardson number (the ratio of stability to shear) fell as low as 0.1 (i.e. unstable) although the results show that bottom-boundary driven turbulence can dominate for periods. The implication, based on scaling, is that the mixing within the channel does not homogenise the water column within a tidal cycle. Scaling, developed to characterise the tidal advection relative to the channel length, shows how riverine-driven buoyancy fluxes can pass through the tidal channel section and the stratification can remain partially intact.     

Three-dimensional modeling of tidal circulation and mixing over the Java Sea

A combination of a three-dimensional hydrodynamic model and in-situ measurements provides the structures of barotropic tides, tidal circulation and their relationship with turbulent mixing in the Java Sea, which allow us to understand the impact of the tides on material distribution. The model retains high horizontal and vertical resolutions and is forced by the boundary conditions taken from a global model. The measurements are composed of the sea level at coastal stations and currents at moorings embedded in Seawatch buoys, in addition to hydrographic data. The simulated tidal elevations are in good agreement with the data for the K₁ and M₂ constituents. The K₁ tide clearly shows the lowest mode resonance in the Java Sea with intensification around the nodal point in the central region. The M₂ tide is secondary and propagates westward from the eastern open boundary, along with a counterclockwise amphidromic point in the western part. The K₁ tide produces a major component of tidal energy, which flows westward and dissipates through the node region near the Karimata Strait. Meanwhile, the M₂ tide dissipates in the entire Java Sea. However, the residual currents are mainly induced by the M₂ tide, which flows westward following the M₂ tidal wave propagation. The tidal mixing is mainly caused by K₁ tide which peaks at the central region and is consistent with the uniform temperature and salinity along the vertical dimension. This mixing is expected to play an important role in the vertical exchange of nutrients and control of biological productivity.