Effect of the baroclinic tide on impurity diffusion in the ocean

The paper studies the effect of the baroclinic tidal wave on the diffusion of a patch of impurities in a continuously stratified fluid. The turbulent diffusion equation is solved numerically, with wave currents being considered. The diffusion characteristics depend upon the parameters of a baroclinic wave and upon the location of the patch. Translated by Vladimir A. Puchkin.     

苏北辐射沙洲海域水文特征与COD潮扩散计算分析

采用平面二维水动力及污染物输运模型.针对研究海域的水动力特性.应用潮流场作为背景场驱动控制污染物浓度的对流扩散方程,利用有限元方法对苏北辐射沙洲海域大、中、小潮水动力特征及如东洋口港附近污染物输运进行了计算分析。计算结果与实测结果比较表明:计算所得潮流场合理.辐射沙洲区部分潮流脊间水道的潮流运动具有明显的定向往复流形式,潮流强,潮差大.且形成了辐射状潮流流炀;排污口排放的污染物的扩散影响范围不大,不利潮型为小湖,大潮和中潮情况下,污染面积有所减小。     

One-Dimensional Unsteady Analytical Solution of Salinity Intrusion in Estuaries

Based on the one-dimensional salinity transport equation with constant diffusion coefficient, and separated water flow velocity into runoff and tidal current with the single-frequency in an idealized estuary, the simplest unsteady analytical so- lution of salinity intrusion is deduced and the estimation formula of diffusion coefficient is obtained in this paper. The unsteady solution indicates that salinity process in estuaries results from the interaction of runoff and tidal current, and its amplitude is in direct proportion to the product of the velocity of runoff water and the amplitude of tidal flow velocity and in inverse proportion to the diffusion coefficient and the tidal angular frequency, and its phase lag tidal flow with 7/2 which reveals the basic features of the maximum salinity appearing after flood slack and the minimum salinity appearing before ebb slack under the effect of runoff （the advance or lag time is relative to the magnitude of runoff and tidal flow）. According to the measured flow velocity and salinity data, the salinity diffusion coefficient could be estimated. Finally, with the field data of observing sites on the deepwater navigation channel of the Yangtze Estuary, the diffusion coefficient is calculated and a comparative analysis of simulated and measured of salinity process is made. The results show that the solution can comprehensively reflects the basic characteristics and processes of salinity intrusion under the interaction of runoff and tidal flow in estuaries. The solution is not only suitable for theoretical research, but also convenient for estimating reasonable physical parameters and giving the initial condition in the salinity intrusion numerical simulation.     

Undular tidal bore dynamics in the Daly Estuary, Northern Australia

Measurements in the macro-tidal Daly Estuary show that the presence of an undular tidal bore contributed negligibly to the dissipation of tidal energy. No recirculation bubble was observed between a trough and the following wave crest in the lee waves following the undular bore. This differs to stationary undular bores in laboratory experiments at larger Froude numbers where a recirculation bubble exists. Secondary motions and the turbulence generated by the undular bore had no measurable influence on the sediment transport. This situation contrasts with the intense sediment resuspension observed in breaking tidal bores. The tidally averaged sediment budget in the Daly Estuary was controlled by the asymmetry of tidal currents. The undular bore may widen the river by breaking along the banks that it undercuts, leading to bank slippage. A patch of river-wide macro-turbulence of 3-min duration occurred about 20 min after the passage of the bore during accelerating tidal currents.     

水质扩散系数在伶仃洋水域水体交换中的影响分析

选取珠江河口伶仃洋水域为例,建立了贴体曲线坐标系下的二维珠江河口水流数学模型,并采用纯隐格式的混合有限分析法对数学模型进行了离散和求解。在对模型验证的基础上,针对珠江河口伶仃洋水域中的水质扩散系数选取问题进行了分析和探讨。研究表明,主潮汐通道水体的置换率基本不受扩散系数取值的影响,扩散系数的影响主要体现在远离伶仃洋主潮汐通道和示踪剂浓度分界线的水域,并初步推荐了该水域示踪剂扩散系数的取值范围。     

Critical Latitude in Tidal Dynamics Using the Kara Sea as an Example

It is well known that, within the linear nonviscous equations of tidal dynamics, the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components unlimitedly increase when approaching the critical latitude. It is also known that the linear equations of tidal dynamics with a constant and specified vertical eddy viscosity indicate the occurrence of significant tidal velocity shears in the near-bottom layer, which are responsible for increasing the baroclinic tidal energy dissipation, the turbulent kinetic energy, and the thickness of the bottom boundary layer. The first circumstance—the growth of the amplitudes of oscillations of the barotropic and baroclinic tidal velocity components—is due to the elimination in the original equations of small terms, which are small everywhere except for the critical latitude zone. The second circumstance—the occurrence of significant tidal velocity shears—is due to the fact that internal tidal waves, which induce the dissipation of the baroclinic tidal energy and the diapycnal diffusion, are either not taken into account or described inadequately. It is suggested that diapycnal diffusion can lead to the degeneration (complete or partial) of tidal velocity shears, with all the ensuing consequences. The aforesaid is confirmed by simulation results obtained using the QUODDY-4 high-resolution three-dimensional finite-element hydrostatic model along the 66.25° E section, which passes in the Kara Sea across the critical latitude.     

Turbulent momentum and salt transport in the mixing zone of the Elbe estuary

Measurements in the mixing zone of the Elbe estuary were performed during three consecutive tidal cycles with three types of instruments—a moored tripod with velocity and temperature/conductivity/light attenuation sensors, a profiling sonde with similar sensors lowered from an anchored vessel, and instrumented moorings. Acoustic-travel-time sensors were used for velocity measurements.Spectral analysis of 12·8 min pieces of the obtained time series gives results that are consistent with isotropic turbulence for part of the frequency space. Temporal changes of turbulent kinetic energy are correlated with tidal current velocity. A retardation is found between changes in tidal current and turbulent energy. Not all shear stress terms are in similar phase with tidal flow. Mean gradients, Reynolds stress terms, and turbulent salt flux terms are combined to determine eddy viscosity and eddy diffusion coefficients.     

基于FVCOM 的辽东湾水动力特性数值模拟

辽东湾属于半封闭性海湾, 水动力过程具有一定的代表性。本文基于无结构化三角形网格的有限体积海岸海洋模型FVCOM（Finite-Volume Coastal Ocean Model）, 构建了辽东湾及其邻近海域的三维水动力数值模型, 并利用实测数据对6 个潮位验潮站、4 个潮流验潮站的大小潮时刻的潮位、流速、流向进行了对比验证, 该模型能够准确地模拟辽东湾的潮汐、潮流等水动力场情况, 可进一步为研究辽东湾温盐、泥沙、水质、污染物扩散等提供研究基础。     

象山港电厂温排水的实测和数值模拟研究

为了研究象山港电厂温排水的热污染问题,对其附近海域的潮流和水温扩散进行了实测研究,给出了该海域潮流及水温扩散的特征。同时利用近海潮流模式(ECOMSED)建立了象山港潮流三维数值模型,对研究海域的潮流和水温扩散进行了模拟,得到了大小潮涨、落憩时刻的温升场的扩散范围,对比分析发现数模结果与实测结果符合较好。     

A BFG model for calculation of tidal current and diffusion of pollutants in nearshore areas

ＡＢＦＧｍｏｄｅｌｆｏｒｃａｌｃｕｌａｔｉｏｎｏｆｔｉｄａｌｃｕｒｒｅｎｔａｎｄｄｉｆｆｕｓｉｏｎｏｆｐｏｌｌｕｔａｎｔｓｉｎｎｅａｒｓｈｏｒｅａｒｅａｓ￥ＳｈｉＦｅｎｇｙａｎ；ａｎｄＺｈｅｎｇ；Ｌｉａｎｙｕａｎ（ＳｔａｔｅＫｅｙＬａｂｏｒａｔｏｒｙ...     

浅海水域水平扩散系数的推算研究

本文根据拉格朗日流与欧拉流自相关系数的相似关系,用渤中8号平台1987年1月和8月的实测资料,分别推算了水平扩散系数。推算结果表明:研究长时间水平扩散时,把潮流看作湍流,所求的扩散系数实际上是潮流扩散系数,其值在10~6cm~2/s范围内;研究短时间水平扩散时,把潮流看作平流,用富里叶变换的长周期合成值从实测值中扣除,所求的扩散系数消除了长周期潮流部分,可认为是湍流扩散,其值在10~4cm~2/s范围内。     

琼州海峡的一个三维泥沙扩散模式

通过建立一个三维陆架海数值模式来预报近海潮流和泥沙浓度的对流扩散,并应用于琼州海峡跨海工程的前期数值预测。预测结果表明,海区潮流的涨落、岸线形状、海底地形及源强位置与泥沙在近海的输运扩散关系密切,模拟情形基本合理,表明该模式适合应用于河口、近岸等涉海海洋工程海域的流场及泥沙扩散的三维数值预测,为相关的水动力环境要素和水质影响评价提供一定的参考价值。     

现场海域水平扩散系数的探讨 Ⅱ．黄河三角洲埕北海域的计算

根据Ｔａｙｌｏｒ的湍流扩散理论和Ｈａｙ、Ｐａｓｑｕｉｌｌ提出的拉格朗日流和欧拉流的相关函数的相似性，用黄河三角洲埕北海域的实测资料，按四种不同方案计算出该海域的水平扩散系数和谱密度分布。指出：该海区潮流扩散系数主要受潮流的影响，其分布特征与潮流的分布基本一致，所求的水平扩散系数主要是由周期性的潮运动给予潮流扰动所致。方案４是通过富氏变换滤波消去了周期性潮流扰动，本方案计算的水平扩散系数是由周期为６ｈ以下的小尺度涡提供能量。     

A three-dimensional semi-implicit unstructured grid finite volume ocean model

A new three-dimensional semi-implicit finite-volume ocean model has been developed for simulating the coastal ocean circulation, which is based on the staggered C -unstructured non-orthogonal grid in the horizontal direction and z -level grid in the vertical direction. The three-dimensional model is discretized by the semi-implicit finite-volume method, in that the free-surface and the vertical diffusion are semi-implicit, thereby removing stability limitations associated with the surface gravity wave and vertical diffusion terms. The remaining terms in the momentum equations are discretized explicitly by an integral method. The partial cell method is used for resolving topography, which enables the model to better represent irregular topography. The model has been tested against analytical cases for wind and tidal oscillation circulation, and is applied to simulating the tidal flow in the Bohai Sea. The results are in good agreement both with the analytical solutions and measurement results.     

Effects of crab–halophytic plant interactions on creek growth in a S.W. Atlantic salt marsh: A Cellular Automata model

The Bahía Blanca Estuary (38° 50′ S, and 62° 30′ W) presents salt marshes where interactions between the local flora (Sarcocornia perennis) and fauna (Chasmagnathus granulatus) generate some kind of salt pans that alter the normal water circulation and condition its flow and course towards tidal creeks. The crab–vegetation dynamics in the salt marsh presents variations that cannot be quantified in a reasonable period of time. The interaction between S. perennis plant and C. granulatus crab is based on simple laws, but its result is a complex biological mechanism that causes an erosive process on the salt marsh and favors the formation of tidal creeks. To study it, a Cellular Automata model is proposed, based on the laws deduced from the observation of these phenomena in the field, and then verified with measurable data within macroscale time units. Therefore, the objective of this article is to model how the interaction between C. granulatus and S. perennis modifies the landscape of the salt marsh and influences the path of tidal creeks. The model copies the basic laws that rule the problem based on purely biological factors.The Cellular Automata model proved capable of reproducing the effects of the interaction between plants and crabs in the salt marsh. A study of the water drainage of the basins showed that this interaction does indeed modify the development of tidal creeks. Model dynamics would likewise follow different laws, which would provide a different formula for the probability of patch dilation. The patch shape can be obtained changing the pattern that dilates.     

青岛—石臼所沿岸潮流特征分析

本文对青岛—石臼所沿岸潮流特征进行了分析。实际观测资料表明,这一带潮流具有典型的往复流性质,并有较明显的不对称性和较大的垂直变化。文中指出涨落潮流不对称性是由于余流和分潮流影响的结果。讨论了不对称性随余流和分潮流的变化和潮流垂直梯度的特征。不对称和垂直梯度对泥沙运动和污染扩散的影响在文中也作了阐述。     

Comparison of two three-dimensional hydrodynamic modeling systems for coastal tidal motion

A comparison of two three-dimensional numerical modeling systems for tidal elevations and velocities in the coastal waters is presented. The two modeling systems are: (1) the Princeton Ocean Model (POM) and (2) the MIKE 3 flow model. The model performance results for Singapore's coastal waters show that the predicted tidal elevations from the two hydrodynamic modeling systems are almost identical and are in very good agreement with field measurement data. The simulated tidal current velocities match well with field measurement data at the selected stations, but it seems that the POM provides the slightly better simulation, compared to the MIKE 3 flow model. The depth profiles of the velocities obtained from the two modeling systems may be greatly different at some time, due to the vertical diffusion coefficient calculated from different turbulent sub-models in the two modeling systems. The POM generally predicts larger peak tidal velocities. The maximum speed differences for the model results from the two modeling systems occur in the top and differ from time to time and from location to location, reaching up to 20%.     

黄河口门位置及入海径流变化对莱州湾盐度分布的影响

2020年黄河丰水期入海径流量是往年平均值的2倍以上,必然会引起河口水动力和盐度分布的动态变化。本作者基于有限体积海岸海洋模型(Finite Volume Community Ocean Model, FVCOM),模拟2020年黄河冲淡水在丰水期和枯水期的扩散情况,研究黄河口以及莱州湾海域的盐度分布状况变化,以及径流量变化和口门变迁对黄河冲淡水扩散的影响,模型结果与观测结果吻合较好。模拟结果表明,黄河口西北侧潮流沿岸线方向,随着涨落潮呈西北-东南向往复;黄河口以南包括莱州湾的潮流均随着涨落潮呈东北-西南方向往复。高流速区域主要集中在黄河口和莱州湾北部,在0.5m/s以上。在余流作用下,大量的黄河冲淡水会涌入莱州湾,丰水期时27psu等盐线包络面积占到整个莱州湾的1/4左右。径流量和风的变化主要影响羽流的扩散面积,而口门的变迁会改变其扩散方向。黄河冲淡水经北向口门入海主要影响莱州湾区域,经东向口门入海更多地会向北扩散。通过对2020年黄河口及莱州湾海域盐度分布的分析,为黄河入海径流管理及莱州湾渔业资源保护提供科学参考。     

河控型河口盐度层化对悬沙的捕集机制

河控型河口盐度混合和层化是控制悬沙输移扩散的重要动力机制。以珠江磨刀门河口为研究对象,基于2017年洪季三船同步大、小潮水文泥沙观测数据,分析河控型河口水体盐度层化结构的时空变化对悬沙分布的影响机制。结果表明:受径潮动力耦合时空变化影响,河口盐度垂向分布表现出时空差异,即受径流主导的M1站(挂锭角),河口盐度在涨落潮周期内垂向混合均匀,受径潮控制的M2站(口门)在整个潮周期内盐度层化结构明显,口门外侧的M3站,潮动力作用较强,盐度垂向分布随涨落潮变化而变化;悬沙空间分布与盐度分布关系密切,盐度混合均匀利于悬沙垂向均匀分布,而盐度层化则使悬沙倾向于滞留在底层水体中,且在盐度层结界面之下出现高悬沙浓度,悬沙浓度垂向分布曲线呈L字型或抛线型,纵向上表现为高浓度悬沙团抑制在盐水楔前端,盐度层化对悬沙的捕集效应明显。通过对比水体标准化分层系数与水流垂向扩散强度系数发现,两者呈现负相关关系,即标准化分层系数愈大,垂向扩散强度愈小,表明水体层化抑制悬沙垂向扩散强度,而且水体层化程度越高,悬沙垂向扩散抑制程度越大,进而促进了河口水体盐度层化对悬沙捕集作用。本研究有助于揭示河口细颗粒泥沙运动机制及河口拦门沙演变机制,并为磨刀门河口拦门沙治理提供科学依据。     

LAGCARTW: A random walk particle advection-diffusion model

LAGCARTW (Lagrange Cartesian suite for Windows) provides a simple method for simulating and understanding diffusion using computer animations of random walk particles. The model is useful for teaching and demonstration purposes, but may also be used in an investigative mode—the suite optionally produces output files of particle positions which can be used in subsequent analysis. Current velocities (defined as a series of tidal harmonics) are defined over the model domain, and particles can be discharged from a series of outfalls, reflected or absorbed by boundaries, and allowed to decay. The program also allows a variable diffusion coefficient in the vertical plane. The input file can be altered quite easily via a separate interactive program. A series of demonstration files, used to illustrate examples of diffusion problems, are provided.