Modeling residual circulation and stratification in Oujiang River estuary

A 3D,time-dependent,baroclinic,hydrodynamic and salinity model was implemented and applied to the Oujiang River estuarine system in the East China Sea.The model was driven by the forcing of tidal elevations along the open boundaries and freshwater inflows from the Oujiang River.The bottom friction coefficient and vertical eddy viscosity were adjusted to complete model calibration and verification in simulations.It is demonstrated that the model is capable of reproducing observed temporal variability in the water surface elevation and longitudinal velocity,presenting skill coefficient higher than 0.82.This model was then used to investigate the influence of freshwater discharge on residual current and salinity intrusion under different freshwater inflow conditions in the Oujiang River estuary.The model results reveal that the river channel presents a two-layer structure with flood currents near the bottom and ebb currents at the top layer in the region of seawater influenced on north shore under high river flow condition.The river discharge is a major factor affecting the salinity stratification in the estuarine system.The water exchange is mainly driven by the tidal forcing at the estuary mouth,except under high river flow conditions when the freshwater extends its influence from the river’s head to its mouth.     

A three-dimensional hydrodynamic and salinity transport model of estuarine circulation with an application to a macrotidal estuary

A three-dimensional semi-implicit finite volume numerical model has been developed and applied to study tidal circulation and salinity stratification in the region of Oujiang River Estuary, China. The model employs horizontally unstructured grids and boundary-fitted coordinate system in the vertical direction. Governing equations consisting of continuity, momentum, and transport equations are all solved in the integral form of the equations, which provides a better representation of the conservative laws for mass, momentum, and transport in the coastal region with complex geometry and bottom bathymetry. The model performance was firstly quantified with skill assessment statistics on the choice of different parameters and validated with observed tidal elevation, current velocity, direction and salinity data over a spring–neap tidal cycle collected in 2006. Numerical results show that the model with wetting–drying capability successfully simulated the tidal currents and salinity fields with a reasonable accuracy and indicate that the Oujiang River Estuary is a macrotidal estuary with strong tidal mixing. In addition, the model results also show that the Oujiang River Estuary is a well-mixed estuary during spring tide. Then, the numerical simulations were performed to compare the hydrodynamic process and salinity distribution before and after a river training, which was conducted by blocking the south branch of the Oujiang River mouth. The results reveal that with the only north access to the sea, the influence of the blocking project on the flood discharge capacity is limited and the incremental velocity is beneficial to the navigation channel maintenance, although it will cause some scour to the embankment. Furthermore, the redistribution of tidal prism passing in or out the north branch makes a little severe salinity intrusion during high tide or low tide. However, the salinity intrusion is still within acceptable range, although it can cause some adverse effect on water intaking of production and life. The variations of salinity levels in Yueqing Bay situated at the north of the river mouth are not obvious, so the blocking project will not bring damage to local aquiculture. However, significant changes of salinity happen inside or outside of the south branch, so enough attention need to be paid to the changes of environment caused by the salinity variation after the blocking project. Overall, by weighing advantages and disadvantages of the blocking project, it is feasible and the model can be considered as a tool for managing and studying estuarine circulation.     

垂向二维潮流数值模型及其在长江口北槽的应用

应用变网格有限元方法，采用任意四边形等参单元，离散横向积分的Navier—Stokes方程，建立河口潮流垂向二维数学模型。应用此模型，对长江口北槽水域两个观测站的潮流水位、流速垂线分布和整个北槽潮流流速纵向分布进行了模拟。潮流水位、流速垂线分布的模拟值分别与观测站的实测值可以较为吻合，证明了本文模型的实用有效。模拟得到的涨急、落急时刻整个北槽潮流流速纵向分布给出了定性的结果。     

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

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

黄河河口演变(Ⅰ) --(一)河口水文特征

本文利用20世纪后半个世纪的完整实测资料，研究了河口径流量与输沙量的变化，自50年代至90年代呈明显递减。对黄河口的水文特征：水情、沙情、断流、离子流量、滨海潮汐、潮流、余流、温、盐度分布、风暴潮、拦门沙与盐水楔等作出定性和定量的论述。对断流的危害、成因及缓解对策作出分析。河口淤积、延伸、改道对黄河下游的影响也作出时空上的定量分析。对1855年以来近代黄河三角洲海岸线变化，作出了定量研究。     

河口水流结构分解及斜压M4分潮流探讨

河口水流是河口生态环境、河道演变、物质输运等物理过程的根本动力。由于径流、潮波、地形以及气象等因素的影响,河口水流呈现复杂的三维结构。其中既包括淡水注入形成的余流,也包括周期性的潮流、风生流、斜压流及河口非线性作用导致的流动等。为探究河口水流的组成及其潮内变化,基于瓯江口实测资料,利用主成分分析（Principal Component Analysis, PCA）法对河口水流进行分解,探讨了PCA法对河口水流的分解性能及斜压分潮流的高频特征。研究认为,PCA法在河口水流结构研究中既可采用原始数据操作亦可用标准化的数据进行计算。PCA法可分解出斜压成分（河口重力环流型结构）,但不能将正压成分（径流和潮流）分开,径流和潮流二者综合作用的结果体现在主成分的得分之中。主成分的取舍应根据水流结构和累计解释方差综合判断,不宜仅依据累计解释方差。河口斜压流动具有明显的高频特征,近似呈1/4日分潮的周期。     

A Three-dimensional Tidal Model in Boundary-fitted Curvilinear Grids

In hydrodynamic models of marine and estuarine currents the use of boundary-fitted curvilinear grids not only makes the model grids fit to the coastline and bathymetry well, but also makes the kinetic boundary conditions simple and more accurate. Because of these advantages, a three-dimensional tidal model with boundary-fitted curvilinear grids has been developed to simulate both tide and current in estuarine and shelf water. The basic idea is to use a set of coupled σ-stretched and elliptic transformations to map the physical space into a corresponding transformed space such that all boundaries are coincident with co-ordinate lines and the transformed grids are rectangular. The hydrodynamic equations in the transformed space are solved in a rectangular mesh of the transformed grid system. The application of the model in the simulation of the M₂,S₂ , K₁and O₁tidal waves in the Bohai Sea, China, shows that the numerical results are in good agreement with the observations. The method of boundary-fitted curvilinear grids is effective for improving simulation accuracy of current in the estuarine and shallow seas, especially at the coastal regions where the current was usually impractical by a uniform Cartesian grid system with the shoreline and bathymetry represented by numerous stair-steps.     

A finite-element,multi-scale model of the Scheldt tributaries,river, estuary and ROFI

We report on the development and validation of a coupled two- and one-dimensional finite-element model for the Scheldt tributaries, river, estuary and region of fresh water influence (ROFI). The hydrodynamic equations are solved on a single, unstructured, multi-scale mesh stretching from the shelf break to the Scheldt tributaries. The tide is forced on the shelf break and propagates upstream in the riverine network. Upstream boundaries lie on sluices or outside of the region of tidal dominance where daily averaged discharges are imposed. Two-dimensional, depth-averaged shallow water equations are solved by means of the discontinuous Galerkin (DG) method over the marine and estuarine parts of the computational domain. In the rivers, however, one-dimensional equations are dealt with using the DG method with the addition of a technique to cope with confluence points. Model parameters are carefully calibrated, leading to the simulation of wind- and tide-forced flows that are in excellent agreement with available data. The diffusivity in the transport equation is calibrated using time series of salinity at various locations in the estuary. Finally, the Lagrangian residual transport in the estuary and the adjacent coastal zone is investigated. This work is a major step towards an integrated model for studying the dynamics of waterborne contaminants and the water renewal timescales in the Scheldt land-sea continuum.     

Prediction of salinity intrusion in the sheltered estuary of Terengganu River in Malaysia using 1-D empirical intrusion model

Generally one dimensional(1-D) empirical salinity intrusion model is limited to natural alluvial estuary. However,this study attempts to investigate its ability to model a sheltered alluvial estuary of the Terengganu River in Malaysia. The constructed breakwater at the mouth of the river shelters the estuary from direct influence of the open sea. The salinity density along the estuary was collected during the wet and dry seasons for scenarios before and after the constructed breakwater. Moreover, the freshwater discharges, tidal elevations and bathymetry data were also measured as model inputs. A good fit was demonstrated between simulated and observed variables,namely salinity distribution and intrusion length for both scenarios. Thus, the results show that 1-D empirical salinity model can be utilized for sheltered estuarine condition at the Terengganu Estuary, but with an appropriate determination of an initial point. Furthermore, it was observed that the salinity intrusion in the study area is largely dependent on the freshwater discharge rather than tidal elevation fluctuations. The scale of the salinity intrusion length in the study area is proportional to the river discharge of the –1/2 power. It was appeared that the two lines of the 1-D empirical salinity model and discharge power based equation fitted well to each other, with the average predicted minimum freshwater discharge of 150 m3/s is going to be required to maintain acceptable salinity levels during high water slack(HWS) near the water intake station, which is located at 10.63 km from river mouth.     

A three-dimensional hydrodynamic model of estuarine circulation with an application to Southampton Water, UK

A three-dimensional hydrodynamic model has been developed to simulate water mass circulation in estuarine systems. This model is based on the primitive equation in Cartesian coordinates with a terrain-following structure, coupled with a Mellor–Yamada 2.5 turbulence scheme. A fractional-step method is applied and the subset of equations is solved with finite volume and finite element methods. A dry–wet process simulates the presence of the tidal flat at low water. River inputs are introduced using a point-source method. The model was applied to a partially mixed, macrotidal, temperate estuary: Southampton Water, UK. The model is validated by comparisons with sea surface elevation, ADCP measurements and salinity data collected in 2001. The mean spring range 2(M₂ + S₂) and the mean neap range 2(M₂ − S₂) are modelled with an error relative to observation of 12 and 16%, respectively. The unique tidal regime of the system with the presence of the ‘young flood stand’ corresponding to the slackening conditions occurring at mid flood and ‘double high water’ corresponding to an extension of the slackening conditions at high tide is accurately reproduced in the model. The dynamics of the modelled mean surface and bottom velocity closely match the ADCP measurements during neap tides (rms of the difference is 0.09 and 0.01 m s⁻¹ at the bottom and at the surface, respectively), whereas at spring the difference is greater (rms of the difference is 0.25 and 0.20 m s⁻¹ at bottom and surface, respectively). The spatial and temporal variation of the degree of stratification as indicated by salinity distributions compares well with observations.     

Three-dimensional tidal current numerical model of the Oujiang Estuary

The characteristics of three-dimensional (3-D) tidal current in the Oujiang Estuary are investigated according to in situ observations. The Oujiang Estuary has features of irregular coastline, complex topography, many islands, moveable boundary, and submerged dyke, therefore, a 3-D numerical model on an unstructured triangular grid has been developed. The σ coordinate transformation, the moveable boundary and submerged dyke treatment techniques were employed in the model so it is suitable for the tidal simu...     

黄河口二维潮波泥沙有限元数学模型及应用(Ⅰ)——模型及其验证

分析了建立黄河河口泥沙有限元模型的必要性,并结合黄河河口实际情况,建立了黄河河口二维潮流和泥沙数学模型。从数学模型模拟计算得到的潮位、流速和冲淤验证结果可以看出,由模型模拟的潮汐和潮流的等海洋动力特性以及泥沙冲淤和海底变形都与实测资料一致,表明建立的河口泥沙二维有限元数学模型可进行黄河河口的潮汐、潮流和泥沙运动及海底冲淤演变的模拟与预测。     

闽江口的盐、淡水混合

采用盐、淡水混合比的概念研究闽江口的混合状态，并确定了分层、部分混合、高度混合３种类型的分类指标。研究表明，闽江口主要水道的盐、淡水混合以部分混合为主，其次是分层，高度混合出现的机率较少；盐水异重流发育良好，混合状态的变化与潮流速的变化不相一致；盐水入侵时混合比相对较小，后退时相对较大。闽江口是一个径、潮量比值较小的强潮河口，按照一般规律，其盐、淡水混合应该出现高度混合型，但在闽江口的主汊道中发生     

北部湾大风江与南流江河口红树林空间分布格局研究

生长在潮间带的红树植物在河口植物群落构成、海岸防风消浪中具有重要价值。本文基于本地种桐花树胚胎浸泡下沉实验与北部湾南流江和大风江河口段水体盐度、沿线潮间带植物群落结构与地貌分析,探讨红树林在河口空间分布及影响因素。结果主要表明:南流江河口和大风江河口红树林自海向陆基本展现"红树林纯林(桐花树、秋茄、无瓣海桑种类混生)→红树植物与半红树植物(黄槿、苦朗等)混生→红树植物、半红树植物与非红树植物混生→红树植物镶嵌→稀疏红树林小苗"的分布格局,但大风江河口向陆界限主要以红树、红树幼苗及半红树混生为主。此外,红树被浸淹时长是控制河口红树空间分布结构的主要因素。潮水上溯时长影响红树向陆生长的极限位置,宜林滩地是红树发育生长的必要条件。     

A box model of Puget Sound

A classical two-layer box model has been used to calculate volume transports and vertical exchange coefficients for the Main Basin of Puget Sound. High river flow (January–February) and low flow (August–September) calculations, using salinity and runoff observations, show that basically two estuarine types exist within the basin under both flow conditions. Admiralty Inlet, the north entrance to the Main Basin, is similar to a partially mixed estuary with vigorous tidal mixing, so that horizontal and vertical salinity gradients are similar in winter and late summer. Within the deep main basin, two layer transports are proportional to salinity stratification and the vertical exchanges are low. Calculated summer transports are about a factor of two smaller than winter transports for the Main Basin. Model transports agree quite well with daily net transports estimated from current meters. Flushing times calculated by the model also agree with volume replacement times calculated using current meter data and by methods using oxygen deficits in the lower layer.     

Yangtze Brackish water plume—circulation and diffusion

A system of 3-D linearized momentum equations, the continuity equation and a simplified equation of state for sea water are analytically solved. The solution comprises three parts. The first part is a wind-driven current, the second a jet-like current caused by Yangtze River outflow and the last a density flow. The computed current velocity is obtained by solving an advection-diffusion equation of salinity by a finite-difference method. The results show that the Yangtze brackish water plume spreads as a large low salinity tongue. At the surface, the turning of the tongue axis is mainly controlled by wind stress, Yangtze run-off and bottom topography. Near the bottom, the axis changes little all the year around. At about 10–15 m depth, there is a sharp halocline in summer, but at deeper layers the vertical distribution is almost always uniform. The computed salinity distribution agrees fairly well with the observed one.     

珠江伶仃河口湾及邻近内陆架的纵向环流与物质输运分析

根据实测资料分析了珠江伶仃河口湾与邻近内陆架在不同径流影响下的水体混合空间状态和季节变化特征。主要由河口湾表、中层冲淡水和内陆架底层上溯的高盐补偿流构成了河口湾和内陆架之间的净环流,在高径流量和西南大风的情况下,环流下移向内陆架扩展增强。计算表明,7月纵断面上的净环流输运是盐分纵向净通量的控制因素,输运方向指向上游,而其余季节则以向海的净平流输运为主。悬沙净通量主要受净平流及潮抽吸输运控制,潮抽吸输运强度与大小潮有较密切关系。     

Coastal currents adjacent to the Caeté Estuary,Pará Region,North Brazil

We examined and compared tidal currents and water column structure between a near-shore station (12 km from the coast) and an offshore station (32 km from the coast) adjacent to the Caeté River, Pará Region, Brazil. Although the coastal system of Pará is largely influenced by local tides and wind, we found substantial differences in the dominant forcing agents between stations. Water column dynamics at the near-shore station were largely affected by local tidal processes, while differences between surface and bottom layer flows also indicated the importance of gravitational circulation at this station and a substantial influence of the adjacent Caeté River discharge. In comparison, at the offshore station, water column structure was largely influenced by a semi-diurnal tidal flow, an along-coastal current flow (mainly associated with the North Brazil Current) and the dynamics of local wind flow. The near-shore station at low tide showed a high level of stratification; at high tide such stratification was reduced. In comparison, stratification was only apparent within the upper 6 m at the offshore station, the rest of the water column was relatively well-mixed. The stratification within Station 1 at low tide was a result of the bi-directional movement of water discharged from the Caeté River, with lower salinity surface water and high salinity bottom water resulting in an estuarine-like circulation environment. The spatial variability and lack of correlation in current flow and water column structure between the near-shore and offshore stations suggest that a flow field resulting from differences in local circulation, tidal variability and wind persistence separate areas. We argue that this separation may indicate that the offshore station is located in a transition region between the Caeté River waters and the local coastal area.     

A coastal ocean model of semi-implicit finite volume unstructured grid

A two-dimensional coastal ocean model based on unstructured C-grid is built, in which the momentum equation is discretized on the faces of each cell, and the continuity equation is discretized on the cell. The model is discretized by semi-implicit finite volume method, in that the free surface is semi-implicit and the bottom friction is implicit, thereby removing stability limitations associated with the surface gravity wave and friction. The remaining terms in the momentum equations are discretized explicitly by integral finite volume method and second-order Adams-Bashforth method. Tidal flow in the polar quadrant with known analytic solution is employed to test the proposed model. Finally, the performance of the present model to simulate tidal flow in a geometrically complex domain is examined by simulation of tidal currents in the Pearl River Estuary.     

长江口三维潮流数值计算及动力分析

采用三维动力学数学模型对长江口潮流场进行了计算，对长江口潮流流态在一个潮周期内逐时进行了分析，并对长江口南支盐水倒灌从动力学角度进行了初步探讨。计算结果表明：文中模型较好地重演长江口潮流场，可以用于大型复杂河口工程实际中的潮流场的计算和分析。