Field studies of velocity,salinity and suspended solids concentration in a shallow tidal channel near tidal flap gates

The design and operation of mathematical models of solute mixing and sediment transport in estuaries rely heavily on the provision of good-quality field data. We present some observations of salinity, suspended sediment concentration and velocity at one of the tidal limits of a semi-enclosed tidal lagoon in Southern England (Pagham Harbour, West Sussex, UK) where the natural processes of tidal incursion and solute mixing have been heavily modified as a result of the construction of sea walls dating back to the 18th Century. These observations, made immediately downstream of two parallel tidal flap gates by conductivity-temperature-depth (CTD) profiler, and also using bed-mounted sensor frames to measure velocity at 2 fixed depths, have yielded a set of results covering 11 tidal cycles over the period 2002–04. It is clear from the results obtained that over a typical tidal cycle, the greatest vertical salinity gradients occur in the 1–2 h immediately after the onset of the flood tide, and that subsequently, energetic mixing acts to rapidly break down this stratification. Under moderate-to-high fresh water flows (>0.5 m³/s), the break-down in vertical salinity gradient is more gradual, while under low fresh water flows (<0.2 m³/s), the vertical salinity gradient is generally less pronounced. Estimates of Richardson number during the early flood-tide period reveal values that vary rapidly between <1 and about 20, with lower values occurring after around 1.5–2 h after low water. Observations of suspended sediment concentration vary widely even for similar tidal and fresh water flow conditions, revealing the possible influence of wind speed, the storage effects of the water in the lagoon downstream of the observation site, and the complexity of the hydrodynamics downstream of tidal flap gates. The data also show that most of the sediment transport is landward, and occurs during flood tides, with estimated total tidal landward flood tide flux of fine sediment of the order of 50–120 kg under low fresh water flow conditions. These observations, which reinforce the results presented in Warner et al. (2004) and elsewhere, can help to provide information about the appropriate techniques for managing sediments and pollutants, including nutrients from sewage effluent waters, in estuaries where hydraulic flap gates are used to control the entry of fresh water over the tidal cycle.     

基于FVCOM的泉州湾海域三维潮汐与潮流数值模拟

基于FVCOM海洋数值模式,采用非结构的三角形网格和有限体积法,建立了泉州湾海域高分辨率(26 m)的三维潮汐、潮流数值模型。模拟结果同2个验潮站和3个连续测流站的观测资料符合良好,较好地反映了泉州湾内潮汐、潮流运动的变化状况和分布特征,给出了M₂、S₂、K₁、O₁ 4个主要分潮的同潮图、表层潮流椭圆分布,以及模拟区域内最大可能潮差、表层最大可能潮流流速和潮余流分布。分析表明,4个分潮的最大潮汐振幅和迟角差分别为219 cm和19°,85 cm和25°,26 cm和12°,26 cm和9°;石湖港以东海域的潮波为逆时针旋转的驻波,以西海域为前进波;最大可能潮差由湾口的8.0m向湾内增加至8.8 m。湾内潮流类型为规则半日潮流,落潮最大流速大于涨潮最大流速,北乌礁水道为强流区,表层最大可能潮流流速为2.4 m/s;湾口潮流运动以逆时针方向的旋转流形式为主,湾内的潮流运动以往复流形式为主,长轴走向主要沿着水道方向,与等深线和海岸线平行;四个分潮流表层最大流速分别为1.4 m/s,0.58 m/s,0.12 m/s,0.10 m/s。余流流速大小与潮流强弱有密切的联系,表、中、底层最大余流流速分别为26 cm/s,20 cm/s,16 cm/s,三者在水平方向基本呈北进南出的分布形态。     

基于调和分析法与ANFIS系统的综合潮汐预报模型

港口沿岸地区以及河流入海口等地区的精确潮汐预报对于各种海洋工程作业有着非常重要的意义。潮汐水位的变化受到众多复杂因素的影响,而且这些复杂的因素往往有着较强的实变性和非线性。为了进一步提高沿岸港口码头等水域的潮汐水位的预测精度,本文提出了一种基于调和分析模型与自适应神经模糊推理系统相结合的模块化潮汐水位预测模型;并采用相关分析确定整个预测模型的输入维数;模块化将潮汐分解为两部分：由天体引潮力形成的天文潮部分和由各种天气以及环境因素引起非天文潮部分。其中调和分析法用于天文潮部分的预测,ANFIS用于预测具有较强非线性的非文潮部分。模块化综合了两种方法的优势,即调和分析法能够实现长期、稳定的天文潮预报,ANFIS能够以较高的精度实现潮汐非线性拟合与预测。模型使用ANFIS模型和调和分析模型分别对潮汐的非天文潮和天文潮部分进行仿真预测,然后将两部分的预测结果综合形成最终的潮汐预测值。此外,本文选用三种不同的模糊规则生成方法（grid partition （GP）,fuzzy c-means （FCM） and sub-clustering （SC））生成完整的ANFIS系统,并使用实测数据进行验证用以选取最优的ANFIS预测模型。最后将最优的ANFIS模型与调和分析模型相结合进行潮汐水位的最终预报。仿真实验选用Fort Pulaski潮汐观测站的实测潮汐值数据进行预报的仿真实验,仿真结果验证了该模型的可行性与有效性并取得了良好的效果,具有较高的预报精度。     

Tidal Propagation in the Mandovi–Zuari Estuarine Network, West Coast of India: Impact of Freshwater Influx

The Mandovi–Zuari estuarine network on the west coast of India consists of shallow strongly converging channels, that receive large seasonal influx of fresh water due to the monsoons. The main channels, the Mandovi and Zuari estuaries, connect the network to the Arabian Sea. Observations show that tidal amplitude in the channels remains unchanged over large distances (|mS40 km) from the mouths of the main channels and then decays rapidly over approximately 10 km near the head. To understand the dynamics behind this behaviour, a numerical model for tidal propagation has been used that simulates the observed tidal elevations well. Momentum balance in the model is predominantly between pressure gradient and friction. In the region of undamped propagation, the model behaviour is consistent with the theory that geometric amplification balances frictional decay leaving the tide unchanged. This balance breaks down near the upstream end, where channels are narrowest, and mean velocity associated with freshwater influx is sufficiently large to prevent upstream propagation of tide. This leads to rapid decay in tidal amplitude. The model also shows that the mean water-level rises in the upstream direction, in the region of the decay.     

Modelling tidal currents on the coast of Portugal

A tide circulation model of the Atlantic coast of the Iberian Peninsula has been constructed. This regional numerical model covers the whole continental shelf. The finite element computational grid is made of some 16,300 triangular elements with sizes ranging between 13 km (on the offshore boundary) and 1 km (near the coast), with local refinements on the continental shelf and in the area of Figueira da Foz. This site was selected as experimental site for the study of waves and currents in the frame of the MAST2/WAVEMOD research project.Boundary conditions along the three oceanic limits of this widely open domain are obtained from the North-Atlantic component of a World Ocean tidal numerical model known as FES94 [Le Provost, C., Genco, M.L., Lyard, F., Vincent, P. and Canceil, P. (1994) Spectroscopy of the world ocean tides from a finite element hydrodynamic model, Journal of Geophysical Research, Topex-Poseidon Special Issue]. A new radiation-like boundary condition has been introduced in the modelling system used (Telemac-2D), which solves the Shallow Water Equations (SWE), in order to interface the two models and to allow for the tidal wave to leave the northern limit without reflection. Model calibration has been performed on the dominating M₂ constituent. The introduction of the astral static potential generating the tide in the SWE improved this regional model.A long duration run (1 month) has been performed, the model being forced by the eight major tide constituents. Harmonic analysis of results has been performed on 17 tide constituents, due to non-linear interactions of constituents on the continental shelf. Comparisons with the FES94 model on one hand, and with a set of coastal tide gauges on the other hand, are good. A database of tidal harmonics is now available for forecasting sea levels and currents in this area.This work has shown that diurnal shelf trapped tide waves exist in some places along the Portuguese continental shelf, which induce diurnal dominant tidal currents in these places (North of Figueira da Foz).     

海平面变化背景下三大河口咸潮入侵特征及变化浅析

基于2009—2019年实测资料分析研究我国三大河口咸潮入侵的特征及变化规律。研究结果表明:2009年以来,珠江口每年都会受到咸潮入侵的威胁,2019年咸潮入侵程度为近8年来最为严重的一年,全禄水厂超标时间超过210 h;钱塘江口除2014年外均发生较强咸潮入侵,其中2019年咸潮入侵程度为近9年最为严重的一年,南星水厂超标时间超过550 h;长江口咸潮入侵程度整体呈下降趋势,2015年后咸潮入侵次数明显减少,平均每年3次,主要原因为长江口北支倒灌影响,2019年咸潮入侵程度为近5年最为严重的一年,青草沙水库超标时间超过470 h。对河口咸潮入侵影响最大的因素为径流,其次是潮汐和海平面变化;受潮周期影响,咸潮入侵具有日变化和半月变化周期;受径流影响,咸潮入侵具有明显的季节变化和年际变化规律,我国三大河口咸潮入侵机制共性和差异性并存。海平面上升和风暴潮增水会加剧磨刀门咸潮入侵灾害,在高海平面期及风暴潮多发期需格外注意咸潮入侵灾害,做好预警应对。     

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

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

江苏如东西太阳沙及烂沙洋海域潮流泥沙数值模拟

基于不规则三角形网格有限差分法并考虑波浪及其破碎作用,建立了平面二维潮流场和泥沙场数学模型.该模型对有望建设成深水码头和深水航道的江苏如东西太阳沙和烂沙洋海域的潮流场和泥沙场进行了细化数值模拟.数值模拟结果表明：（1）本海区潮流基本上是顺深槽流动的往复流,潮流流速大,烂沙洋北水道和西太阳沙附近大潮涨落潮最大流速分别在2 m/s和1 m/s以上;（2）本海区的潮平均水体含沙量在0.5 kg/m^3以下,落潮含沙量大于涨潮含沙量;（3）小浪对水体含沙量影响很小,大浪作用下水体含沙量明显增加.     

不同周期分潮簇对复杂河网潮位分布的影响*

在不规则半日潮为主的河口, 半日和全日周期的分潮簇在很大程度上影响了潮位的振幅。但当潮波向三角洲河网内部传播时, 由于不同周期分潮簇的振幅沿程变化速率存在差异, 特别是半月周期分潮簇振幅作用得到增强, 使得不同周期分潮簇对河网内潮位分布的影响具有明显的区域性特征。通过珠江三角洲河网一维水动力模型得到研究区域内高频且长周期的潮位资料, 基于小波分解重构出主要周期分潮簇的振幅, 采用谱聚类图像分割算法, 通过像素分割计算确定综合相似度, 从图像学等角度论证了不同周期分潮簇与复杂河网潮位分布的空间关联性。研究发现, 半日和全日周期的分潮簇影响了复杂河网口门潮位振幅的量值, 但半月周期的分潮簇却影响了复杂河网整体潮位分布的形态特征。     

Estuarine salt conservation: a Lagrangian approach

Understanding tidal and mean flow transport mechanisms that maintain the estuarine salt balance against seaward transport by river flow is one of the fundamental questions of estuarine oceanography. Previous salt transport studies have failed to resolve this problem for two principal methodological reasons, in addition to the inherent variability of estuaries. First, salt transport expansions used to represent the salt balance have included a large number of terms that could not be related to any theory of estuarine circulation and whose physical meaning was thus unclear. Second, it has proven difficult to measure small mean velocities in the presence of much larger tidal variations. A new analysis method that expresses Lagrangian scalar property transport in terms of time and space variations of Eulerian variables is much simpler than expansions previously used and yields new insights into estuarine transport processes. This approach (like previous expansions) is valid for narrow systems in which cross-channel transport processes are weak and the ratio of tidal amplitude to mean depth is small. A surface boundary condition is used to eliminate one class of troublesome terms. The equivalence of the tidal-cycle average transport of tidal variations in salinity and the Stokes transport of mean salinity is then employed to show that the vertically integrated tidal flow plays no role in the overall salt balance. That is, seaward transport of salt caused by the river flow can only be balanced by landward transport resulting from correlations between shear and stratification in the mean, tidal and overtide flows. In a hypothetical inviscid estuary, tide and overtide flows would be vertically uniform, and outward fluvial salt transport could only be balanced by shear and stratification in the mean gravitational circulation. In a more realistic example with strong friction, the gravitational circulation would be severely damped, and inward transport could only be accomplished by correlations of shear and stratification in the tide and overtide flows.     

潮致泥沙全沙净输运解析模式

泥沙运动作为水流和底床相互作用的纽带,是河流、河口及海岸工程研究的重要内容。在潮波作用明显的河口、海岸地区,周期性的动力作用下的泥沙运动具有往复和可逆的特征,因此研究这类水域的泥沙的净输运更具有实际的意义。基于泥沙输运和流速呈指数关系假设,建立潮流环境下的泥沙全沙净输运的解析解公式,并对该公式的计算结果和数值计算以及数学模型的结果进行了检验和验证,结果表明本研究提出的公式能较好地反应潮流环境下的泥沙净输运。由此,基于本公式采用潮流分潮调和常数可计算得到全沙净输运,并可以分析各分潮流及其相互作用与泥沙净输运的关系。研究结果显示,在受径流影响较大的半日潮河口,S2、MS4、M2三潮相互作用对全沙净输运的贡献显著高于通常的潮流不对称作用,即M2、M4的相互作用。此外,河口区域径流导致的余流对泥沙净输运的贡献不可忽略,特别是在洪季,大径流量条件下往往导致余流较大,其对泥沙净输运的贡献所占比例较大。     

Strong tidal currents observed near the bottom in the Suruga Trough,central Japan

Current measurements carried out at the depth of 4 m above the sea bottom near the northern edge of the Suruga Trough in the early fall of 1985 indicated the existence of strong semidiurnal tidal currents, which were considered to be associated with internal tides. In order to examine the spatial structure of the bottom intensified tidal flow, more detailed current observations were carried out at three or four depths at two stations along the main axis of the Suruga Trough during about 70 days from August to October 1988. We obtained the following results: (1) the variations of the current velocity caused by the semidiurnal and diurnal internal tides are evident in all of the records, and the orientation of the major axis of each tidal ellipse nearly coincides with that of the main axis of the trough; (2) the semidiurnal internal tide is dominant over the diurnal internal tide at 4 m above the sea bottom at both stations; (3) at the northern station the semidiurnal internal tide is dominant over the diurnal internal tide, whereas they are nearly equal at the southern station except at 4 m above the sea bottom; (4) the biharmonic internal tides with 1/3 day and 1/4 day periods, are found near the sea bottom and the major axis of the tidal ellipse is perpendicular to the orientation of the main axis of the Suruga Trough.     

一维水动力模型在河口潮流推算中的应用

针对河口地区潮流特点,建立了潮流量推算一维水动力数学模型。该模型将Sa in t-V enan t方程组进行求解,并根据涨、落潮流的水动力特性,分别进行涨、落潮糙率率定,准确计算出逐时段河道各断面的水位和流量。应用该模型对宁波三江口地区的潮流量进行推算,结果表明,计算的潮流量与该地区实测的潮流量吻合得较好。     

Concurrent simulation of the eddying general circulation and tides in a global ocean model

This paper presents a five-year global simulation of HYCOM, the HYbrid Coordinate Ocean Model, that simultaneously resolves the eddying general circulation, barotropic tides, and baroclinic tides with 32 layers in the vertical direction and 1/12.5° (equatorial) horizontal grid spacing. A parameterized topographic wave drag is inserted into the model and tuned so that the surface tidal elevations are of comparable accuracy to those in optimally tuned forward tide models used in previous studies. The model captures 93% of the open-ocean sea-surface height variance of the eight largest tidal constituents, as recorded by a standard set of 102 pelagic tide gauges spread around the World Ocean. In order to minimize the impact of the wave drag on non-tidal motions, the model utilizes a running 25-h average to approximately separate tidal and non-tidal components of the near-bottom flow. In contrast to earlier high-resolution global baroclinic tide simulations, which utilized tidal forcing only, the simulation presented here has a horizontally non-uniform stratification, supported by the wind- and buoyancy forcing. The horizontally varying stratification affects the baroclinic tides in high latitudes to first order. The magnitude of the internal tide perturbations to sea surface elevation amplitude and phase in a large box surrounding Hawai’i is quite similar to that observed in satellite altimeter data, although the exact locations of peaks and troughs in the modeled perturbations differ from those in the observed perturbations.     

A study of non-linear tidal propagation in shallow inlet/estuarine systems Part II: Theory

The generation of tidal asymmetries is clarified via numerical integration of the one-dimensional equations for channel geometries characteristic of shallow estuaries. Channels without tidal flats develop a time asymmetry characterized by a longer falling than rising tide. This behavior is enhanced by strong friction and large channel cross-sectional area variability over a tidal cycle. Resulting tidal currents have a shorter, intense flood and a longer, weak ebb (flood-dominant). Addition of tidal flats to the channels can produce a longer rising tide and stronger ebb currents (ebb-dominant), if the area of tidal flats is large enough to overcome the effects of time-variable channel geometry. Weaker friction with flats can also produce this asymmetry.Despite the physical complexity of these systems, essential features of estuarine tidal response can be recovered from one-dimensional models. Shallow estuaries are shown to have a system response leading to stable, uniform senses of tidal asymmetry (either flood- or ebb-dominated, due to phase-locking of forced tidal constituents), with down-channel development in magnitude of asymmetry. These concepts are illustrated by modeling idealized representations of tidal channels at Nauset Inlet, MA, and Wachapreague Inlet, VA, which have flood- and ebb-dominance, respectively.     

High-resolution simulations of a macrotidal estuary using SUNTANS

The parallel, finite-volume, unstructured-grid SUNTANS model has been employed to study the interaction of the tides with complex bathymetry in the macrotidal Snohomish River estuary. The unstructured grid resolves the large-scale, O(10 km) tidal dynamics of the estuary while employing 8 m grid-resolution at a specific region of interest in the vicinity of a confluence of two channels and extensive intertidal mudflats to understand detailed local intratidal flow processes. After calibrating tidal forcing parameters to enforce a match between free surface and depth-averaged velocities at several locations throughout the domain, we analyze the complex dynamics of the confluence and show that the exposure of the intertidal mudflats during low tide induces a complex flow reversal. When coupled with the longitudinal salinity gradient, this flow reversal results in a highly variable salinity field, which has profound implications for local mixing, stratification and the occurrence of fine-scale flow structures. This complex flow is then used as a testbed from which to describe several challenges associated with high resolution modeling of macrotidal estuaries, including specification of high resolution bathymetry, specification of the bottom stress, computation of the nonhydrostatic pressure, accurate advection of momentum, and the influence of the freshwater inflow. The results indicate that with high resolution comes the added difficulty of requiring more accurate specification of boundary conditions. In particular, the bottom bathymetry plays the most important role in achieving accurate predictions when high resolution is employed.     

2008年黄海绿潮路径的数值模拟

利用FVCOM模式拉格朗日粒子跟踪模块模拟了2008年5月到7月黄海绿潮漂移路径,根据遥感图像选取粒子跟踪的初始位置和初始时刻,模式中加入M2,S2,K1,01四个主要分潮和Quickscat每日风场资料作为驱动,模拟的粒子运动的主要路径和到达青岛近岸的时间与遥感图像对比都比较吻合。该结论进一步验证了前人关于青岛绿潮来...     

Influence of the nodal tide on the morphological response of estuaries

The nodal tide causes a periodic variation in the mean tidal range, with a period of 18.6 years and amplitude of about 4% of the mean range. As the cross-sectional area of the channel in an estuary has been observed to be related to the tidal prism, such a variation of the tidal range can significantly influence the morphological development of tide-dominated estuaries. Field observations from the Humber Estuary in England confirm that there is indeed an 18.6 year cycle in the total water volume of the estuary. In this estuary, with a relatively small morphological time scale, the response to the nodal tide variation appears to dominate the long-term large-scale sediment budget changes of the estuary. The semi-empirical model ESTMORF was able to reproduce these observations, in terms of the observed amplitude and phase lag (relative to the variation of the tidal amplitude) of the morphological response. It also showed good agreement with the along estuary variation in both the magnitude and phase of the response. To examine the character of the response in more detail, a simplified analytical solution has been derived from the governing non-linear equations. This has been investigated in terms of the sensitivity of the model to the various model parameters and by comparison with the historic data. The results demonstrate how the system response is strongly dependent on the morphological time scale of the estuary. Furthermore, the response is not uniform but varies along the length of the estuary, being most marked at the mouth and decaying up estuary. When examining changes in geomorphological features, such as the intertidal area, an appreciation of the role of the nodal tide is shown to be essential if historic data are to be properly interpreted.     

河口回水动力的演变过程及其影响机制

河口回水区长度的时空演变对防洪、供水、航运等水资源高效开发利用具有重要指导意义。本文在经典河流回水理论的基础上考虑潮汐动力的影响,聚焦潮汐动力引起的回水效应问题,基于一维水动力解析模型,重新定义河口回水区上游界限（回水界）,并以长江河口为例,探讨河口回水动力演变过程及其影响机制。结果表明:回水界距天生港的距离（即回水区长度）与上游流量、外海边界振幅分别具有显著的负相关和正相关关系,且基本为线性关系;回水界对径潮动力的响应比潮区界更为敏感,能有效表征河口感潮河段径潮动力格局演变;长江河口回水界位置具有明显的季节性差异,其春季和秋季回水界位于芜湖附近（春季和秋季分别距天生港419 km和367 km）,冬季回水界位于感潮河段以上区域,夏季潮汐动力引起的回水效应基本可忽略;长江河口潮流因子和径潮相互作用因子控制的余水位梯度的季节性变化,是长江河口回水界位置变动的主导因素。     

ECOMSED模式在杭州湾海域流场模拟中的应用

针对杭州湾独特的喇叭型强潮河口湾的特点,基于Blumberg等(1996)的ECOMSED模式,引入动边界技术,建立杭州湾三维动边界的潮流模型.模型以正交曲线坐标下三维非线性水动力方程为基本方程,应用Mellor和Yamada的2.5阶湍流闭合模型计算紊动黏滞系数,嵌入Grant和Madsen的底边界层模型考虑波浪对底部应力的作用,采用干湿网格法模拟潮流漫滩过程;综合考虑径流,风应力,密度流和M₂,S₂,K₁,O₁四个主要分潮和M₄,S₄,MS₄三个浅水分潮的作用,从而提高杭州湾潮流模拟的精度.通过验潮站调和常数和多次海流连续观测资料的验证,表明该文建立的模型可以更好的用于杭州湾流场的预报模拟.