Numerical Modeling of the Hyperbolic Mild-Slope Equation in Curvilinear Coordinates

The mild-slope equation is familiar to coastal engineers as it can effectively describe wave propagation in nearshore regions. However,its computational method in Cartesian coordinates often renders the model inaccurate in areas with irregular shorelines,such as estuaries and harbors. Based on the hyperbolic mild-slope equation in Cartesian coordinates,the numerical model in orthogonal curvilinear coordinates is developed. The transformed model is discretized by the finite difference method and solved by th...     

Reducing horizontal diffusion errors in σ-coordinate coastal ocean models with a second-order Lagrangian-interpolation finite-difference scheme

When a steep bottom slope exists, it is well known that conventional methods for calculating horizontal diffusion in sigma-coordinate coastal ocean models causes spurious transport (e.g. salinity, temperature, and sediments) and currents. In this study, a second-order accurate finite-difference algorithm and program have been developed to reduce the spurious numerical diffusion errors. In the proposed algorithm, the finite differencing is performed in the x-z coordinate system to approximate the horizontal gradient. Each variable in the finite differential formation is calculated in the sigma-coordinate grid cells using a second-order Lagrangian interpolation polynomial. In conjunction with a stepwise bottom boundary condition, numerical experiments show that the proposed finite-difference scheme considerably reduces numerical errors compared to conventional approaches when dealing with horizontal diffusion over steep topography, which often occurs in coastal oceans and navigation channels.     

Reducing horizontal diffusion errors in σ-coordinate coastal ocean models with a second-order Lagrangian-interpolation finite-difference scheme

When a steep bottom slope exists, it is well known that conventional methods for calculating horizontal diffusion in sigma-coordinate coastal ocean models causes spurious transport (e.g. salinity, temperature, and sediments) and currents. In this study, a second-order accurate finite-difference algorithm and program have been developed to reduce the spurious numerical diffusion errors. In the proposed algorithm, the finite differencing is performed in the x–z coordinate system to approximate the horizontal gradient. Each variable in the finite differential formation is calculated in the sigma-coordinate grid cells using a second-order Lagrangian interpolation polynomial. In conjunction with a stepwise bottom boundary condition, numerical experiments show that the proposed finite-difference scheme considerably reduces numerical errors compared to conventional approaches when dealing with horizontal diffusion over steep topography, which often occurs in coastal oceans and navigation channels.     

Nitrogen and phosphorus species in the coastal and shelf waters of Southeastern North Carolina, Mid-Atlantic U.S. coast

The diversity of small-scale wetlands, high salinity tidal creeks, salt marshes, estuaries, and a wide and shallow shelf with the Gulf Stream close to the break makes the coastal zone of south-eastern North Carolina (U.S.) a natural laboratory for the study of the cycling of nitrogen (N) and phosphorus (P) in coastal and shelf waters. We assessed the summer concentrations, forms, and ratios for each N (total dissolved N, nitrate + nitrite, ammonium and dissolved organic N) and P (total dissolved P, o-phosphate and dissolved organic P) pool as these nutrients travel from tidal creeks, salt marshes and two large estuaries to Long and Onslow Bays. Additionally, we measured ancillary physical (temperature, salinity and turbidity) and chemical (dissolved oxygen, chlorophyll a and pH) water properties. Highest concentrations of all individual N and P compounds were found in the upper parts of each tributary and were attributed to loads from agricultural and urban sources to the coastal watersheds, continuing downstream to receiving estuaries. In all areas, dissolved organic N and P species were predominant constituents of the total dissolved N and P pools (64–97% and 56–93%, respectively). The lower parts of estuaries and surface shelf waters were characterized by oceanic surface values, indicating removal of N and P downstream in all tributaries. The different watershed and hydrological characteristics also determined the different speciation of N and P pools in each estuary. Despite a high level of anthropogenic pressure on the uppermost coastal waters, there is self-regulation in this coastal ecosystem with respect to human perturbations; i.e. significant amounts of the N and P load are retained within estuarine and nearshore waters without reaching the shelf.     

Effect of River Discharge on Bay of Bengal Circulation

The seasonal circulation and mixed layer depths in Bay of Bengal is modeled using the three-dimensional Princeton Ocean Model (POM). Along the coastal boundaries a higher resolution is accomplished using the curvilinear orthogonal grid. Model uses a free-surface and terrain following sigma coordinates. The initial climatological salinity and temperature fields for the model are derived from the World Ocean Atlas-2001(WOA01). The Model is forced with wind stress derived from COADS wind climatology. Bilinear interpolation is used to obtain the initial fields and wind stress to the required model specification. Using the seasonal fields and wind stress the model is integrated for simulating Bay of Bengal circulation. The numerical simulations on climatological scale for monsoon months were conducted to study the evolution of dynamics. The simulations bring out not only the typical characteristic features of fresh water plume along the coast but also intensification of the flow over the monsoon period. The increase in the fresh water flow found to affect only the western parts of the BoB. The opposing currents due to monsoon winds and southward flowing fresh water discharge (FWD) were also delineated. The model results show that the wind stress induced turbulence process is subdued in the presence of strong vertical salinity stratification due to the influence of FWD. The simulated mixed layer depths are in agreement with the reported analytical energy required for mixing values.     

An Application of Finite Volume WENO Scheme to Numerical Modeling of Tidal Current

1 .IntroductionThe development and application of hydrodynamic models to predict circulation in estuarine ,coastal ,and shelf waters has been extremely active over the past three decades (Spauldinget al .,1992 ,1994) . However , most of these models are b…     

海水中腐殖质对溶解态铁的络合及其影响因素

腐殖质(humic substances,简称HS)是地表普遍存在的天然有机物,对海洋中重要的微量营养元素-铁(Fe)的分布及生物地球化学循环具有重要的影响作用。本文对腐殖质的来源、分布及对海水中溶解态铁的迁移转化的影响做了总结,特别论述了其在河口及沿岸水域的行为。大量研究表明河口、沿岸及开放海水中溶解态铁分布的变化可以用腐殖质的浓度及其铁结合能力的变化来解释。腐殖质的络合作用不仅能够阻止溶解态铁(DFe)在河口、沿岸等水域被去除,而且能够通过洋流将DFe迁移至外海及大洋区域,此外还能增加铁的溶解度及对海水中浮游植物的生物可利用性,并且促进铁的氧化还原循环。研究还发现两者之间的络合强度受到盐度、pH等理化因素的影响。盐度是影响HS与DFe配合能力的重要影响因素,盐度增加,导致HS中可以与Fe配合的位点数量降低,配合总量呈现指数降低,而pH的增加可以增加HS与DFe的配合量。另外HS还能影响海水中DFe的氧化还原,并以此影响浮游植物对DFe的吸收利用。因此腐殖质对溶解态铁的有机络合作用是影响其海洋生物地球化学循环的一个重要参数,对进一步研究海水中腐殖质的浓度和分布具有重要的意义。     

杭州湾、长江口余流及其物质输运作用的模拟研究Ⅰ.杭州湾、长江口三维联合模型

根据杭州湾、长江口流场和物质输运的特点,引进正交曲线网格版本ＥＣＯＭ模型,并对斜压梯度力和物质对流扩散的计算作了改进,建立了一个以杭州湾和长江口为整体的三维联合模型,用于潮流、余流和物质输运的计算和研究．     

Benthic indicators and index approaches in the three main estuaries along the French Atlantic coast (Seine, Loire and Gironde)

The three main estuaries in the French Atlantic coast – the Seine, Loire and Gironde, all with high tidal regimes – are interfaces between the continental and the coastal ecosystems. The Seine and Gironde are highly contaminated, whereas the Loire remains in a more natural state. Both the Seine and Gironde have suffered from harbour construction, and as a result, their biological units are extremely compartmentalized. Benthic species and communities have adapted to tolerate temporal physical and chemical changes ( e.g. salinity, substrata, depth, and levels of fine particles and oxygen) and human activities ( e.g. dredging, shipping traffic, and habitat reduction). Although numerous bio-indicators and indices are used to define the ecological quality status (EcoQS) of coastal waters, very few of them were developed specifically for environments with a mosaic of conditions and salinity levels, including freshwater. The main problem appears to be that all the indices for determining anthropogenic stress examine the abundances of stress-tolerant species, which may also be able to tolerate natural stressors such as those occurring in estuaries. This paper takes a look at the development status of the benthic indicators and index approaches used in the three main estuaries along the French Atlantic coast. In addition, it examines the adaptation of the different benthic indicators to the taxonomic sufficiency principle, and the adaptation of the Benthic Opportunistic Polychaetes Amphipods index (BOPA) and the Benthic Opportunistic Annelida Amphipod index (BO2A) for use in the freshwater zones of transitional waters ( i.e. up to the upper limit of the tidal range). Several perspectives are discussed in light of the diverse disturbances and the heterogeneity of such stressed zones, such as the use of multimetric and scoring approaches.     

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary I. A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

ＩＮＴＲＯＤＵＣＴＩＯＮＲｅｓｉｄｕａｌｃｕｒｒｅｎｔａｎｄｉｔｓｉｍｐａｃｔｏｎｍａｓｓｔｒａｎｓｐｏｒｔａｒｅｉｍｐｏｒｔａｎｔｔｏｔｈｅｓｔｕｄｙｏｆｃｏａｓｔａｌｅｎ ｖｉｒｏｎｍｅｎｔ.Ａｌｔｈｏｕｇｈｌｏｔｓｏｆｒｅｓｅａｒｃｈｅｓｈａｖｅｂｅｅｎｄｏｎｅｏｎｔｈｅｒｅｓｉｄｕａｌｃｕｒｒｅｎｔａｎｄｍａｓｓｔｒａｎｓ ｐｏｒｔｉｎｔｈｅＨａｎｇｚｈｏｕＢａｙａｎｄｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ (Ｃａｏ ,1 989;ＣａｏａｎｄＦａｎｇ ,1 986 ;Ｃｈｅｎｅｔａｌ.,1 992 ;ＨｕａｎｄＨ…     

Optimization of the indophenol blue method for the automated determination of ammonia in estuarine waters

Existing automated methods for the determination of ammonia in natural waters suffer from serious ‘salt error’ in estuaries because of changes in pH, ionic strength and optical properties with salinity. A modified automated indophenol blue method is described which minimizes the ‘salt error’ to less than 8% over the entire salinity of estuaries.     

长江口悬沙沉速室内试验研究

在充分认识传统沉降筒缺陷的基础上提出了"大型可温控自动搅拌沉降试验筒"。通过室内系列试验发现:(1)含沙量对长江口细颗粒沉降速度影响最大;(2)温度上升,沉速增加,但不同阶段影响程度有所不同;(3)含沙量越高,盐度对沉速的影响越小,含沙量相同情况下,长江口北槽悬沙枯季水温下盐度对沉速的影响在1.8~5.7倍左右;洪季水温下盐度对沉速的影响在1.5~2.2倍左右。(4)枯季最佳絮凝盐度在7左右,最佳絮凝含沙量为7 kg/m3;洪季最佳絮凝盐度在10~12左右,最佳絮凝含沙量为4.5kg/m3。本研究成果可望加深我们对细颗粒泥沙动力过程相关机理的认识,同时可为相关港口、航道的淤积机理分析,数学模型、物理模型研究工作提供一定技术参考。     

Numerical study on residual current and its impact on mass transport in the Hangzhou Bay and the Changjiang Estuary Ⅰ.A3-D joint model of the Hangzhou Bay and the Changjiang Estuary

A 3-D numerical model is set up in a large domain covering the Hang zhou Bay and the Changjiang Estuary based on the ECOM model in orthogonal curvilinear coordinates.The numerical schemes for baroclinic pressure gradient (BPG) terms and convective terms are improved in the paper according to the characteristics of velocity field and mass transport in the area.The model is validated by the simulations of residual current and salinity transport in the Hangzhou Bay and the Changjiang Estuary.     

A fully nonlinear Boussinesq model in generalized curvilinear coordinates

Based on the fully nonlinear Boussinesq equations in Cartesian coordinates, the equations in generalized coordinates are derived to adapt computations to irregularly shaped shorelines, such as harbors, bays and tidal inlets, and to make computations more efficient in large near-shore regions. Contravariant components of velocity vectors are employed in the derivation instead of the normal components in curvilinear coordinates or original components in Cartesian coordinates, which greatly simplifies the equations in generalized curvilinear coordinates. A high-order finite difference scheme with staggered grids in the image domain is adopted in the numerical model. The model is applied to five examples involving curvilinear coordinate systems. The results of these cases are in good agreement with analytical results, experimental data, and the results from the uniform grid model, which shows that the model has good accuracy and efficiency in dealing with the computations of nonlinear surface gravity waves in domains with complicated geometries.     

Hybrid N-order Lagrangian Interpolation Eulerian?Lagrangian Method for Salinity Calculation

The Eulerian?Lagrangian method (ELM) has been used by many ocean models as the solution of the advection equation, but the numerical error caused by interpolation imposes restriction on its accuracy. In the present study, hybrid N-order Lagrangian interpolation ELM (LiELM) is put forward in which the N-order Lagrangian interpolation is used at first, then the lower order Lagrangian interpolation is applied in the points where the interpolation results are abnormally higher or lower. The calculation results of a step-shaped salinity advection model are analyzed, which show that higher order (N=3?8) LiELM can reduce the mean numerical error of salinity calculation, but the numerical oscillation error is still significant. Even number order LiELM makes larger numerical oscillation error than its adjacent odd number order LiELM. Hybrid N-order LiELM can remove numerical oscillation, and it significantly reduces the mean numerical error when N is even and the current is in fixed direction, while it makes less effect on mean numerical error when N is odd or the current direction changes periodically. Hybrid odd number order LiELM makes less mean numerical error than its adjacent even number order LiELM when the current is in the fixed direction, while the mean numerical error decreases as N increases when the current direction changes periodically, so odd number of N may be better for application. Among various types of Hybrid N-order LiELM, the scheme reducing N-order directly to 1st-order may be the optimal for synthetic selection of accuracy and computational efficiency.     

河口海岸三维双σ坐标斜压水流数值模型研究

在海底地形陡变、垂向密度分层明显的水域,三维σ坐标模式中会出现一种"伪"水平斜压梯度力,并会引起"伪"密度流,以至于影响模拟的精度。垂向上引入双σ坐标变换,建立河口海岸水域三维斜压水流数值模型。数值试验结果表明,在海底地形陡变水域,双σ坐标模式可以减小水平斜压梯度力处理引起的误差。     

秦皇岛河流-海岸水动力和水质耦合模拟分析

基于MIKE FLOOD建立秦皇岛入海河流及其近岸海域一维、二维耦合河流-海岸水动力和水质模型,研究了秦皇岛海域和入海河流的水动力和污染物分布规律。结果表明:(1)海域潮流为半日潮,其中石河口南北支入海口流速变化过程差别较大,但水位变化过程基本一致;(2)整体上河流入海口化学需氧量(COD)与流速呈正相关关系;(3)COD在近岸海域的扩散方向与涨落潮潮流方向一致,各入海河口海域落潮流期间的COD高于涨潮流。如上结论为赤潮爆发机理和应急处置研究提供科学依据和技术支持。     

秋季浙闽沿岸溶解态Cu分布特征及其影响因素

铜(Cu)是海洋浮游植物生长不可或缺的痕量金属之一,对海洋初级生产力起着关键作用。河流是海洋中Cu的重要来源,河口及边缘海对河流输入的Cu起着重要的改造作用,但目前对Cu在浙闽沿岸的生物地球化学行为尚不明确。本研究使用自动固相萃取-电感耦合等离子体联用技术对2021年11月浙闽沿岸及其邻近水域表层水的溶解态铜(dCu)浓度进行分析。结果显示,该区域dCu浓度范围为3.38～26.28 nmol·L^–1,平均浓度为11.66±5.83nmol·L^–1。在研究区域内, dCu的空间分布呈北高南低,近岸高、远岸低的特征。此外, dCu在浙闽沿岸表现出较高的保守性,其与盐度呈显著负相关关系,表明dCu在一定程度上可用于指示人为影响。相关性分析表明,浙闽沿岸及其邻近水域表层dCu的生物地球化学行为和分布可能受到陆源输入、水团输运混合、化学絮凝与吸附等过程的影响。本研究结果有助于进一步理解Cu在海洋中的生物地球化学过程,为探究该区域的生态环境变化提供科学依据。     

辽东湾北部近海大型水母数量分布和温度盐度特征

通过2009-2011年5月下旬-7月下旬辽东湾北部近海10m等深线内的大型水母调查数据,分析了辽东湾北部近海近三年中大型水母资源状况,并探讨了辽东湾大型水母的生态类型。结果显示:辽东湾北部近海大型水母种类主要有海蜇(Rhopilema esculentum)、沙蜇(Nemopilema nomurai)、白色霞水母(Cyanea nozakii)、海月水母(Aurelia sp.1),海蜇和沙蜇是优势种。海蜇幼水母阶段主要集中分布在5m等深线以内的近岸河口水域,随着个体增大有略向深水或密度较小的水域扩散的趋势,仍主要分布在5m等深线两侧水域,属于高温低盐种类。6月份调查海区中发现大量的沙蜇幼水母,随着沙蜇个体增大,7月份调查海区中沙蜇数量大幅度减少。辽东湾海月水母在南部海域出现较多,2010、2011年在北部近海部分海域出现。白色霞水母近几年来辽东湾出现较少,栖息在盐度较高的水域。辽东湾各种大型水母中,沙蜇的生长速度最快。辽东湾海蜇幼水母、沙蜇幼水母的海区出现时间要晚于黄、东海。     

Numerical Tidal Current Modeling Using Fine Boundary Fitted Grids for Tidal Flats

Fine grids with small spacing in boundary-fitted coordinates are normally used to treat thecomputation of fluid dynamics for estuarine areas and tidal flats.However,the adoption of Cartesian com-ponents of velocity vectors in this kind of non-orthogonal coordinates will definitely result in a diffculty insolving implicitly the transformed momentum equations,and also complicate the wet-dry point judgementused for flood areas.To solve this problem,equations in terms of generalized contravariant velocityvectors in curvilinear coordinates are derived in the present study,by virtue of which,anAlternative-Direction-Implicit numerical scheme in non-orthogonal grids would then be easily obtained,and wet-dry point judgement would as well be largely simplified.A comparison is made between the explic-it scheme and implicit scheme,showing that the present model is accurate and numerically stable for com-putations of fluid dynamics for estuarine areas and tidal flats.