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…     

σ-Coordinates hydrodynamic numerical model for coastal and ocean three-dimensional circulation

A three-dimensional finite-difference hydrodynamic model has been developed using σ-coordinate for the vertical dimension. An explicit scheme for temporal integration and a staggered grid for spatial discretization have been adopted. The model has been tested against analytical or literature cases for wind and tide induced circulation. Results are in good agreement both with analytical solutions under idealised conditions and with results from the model of Shankar et al. (1996).     

长江口中肋骨条藻赤潮的分布与特点

对1990年6月9～25日在东海赤潮多发区海域调查期间所观察到的5次中肋骨条藻赤潮的特点以及与环境要素变化的关系作了初步分析。结果表明:(1)中肋骨条藻是多发区海域夏季的敏感赤潮生物,风速减小、气压下降以及良好的海况条件有利于赤潮的发生;(2)赤潮发生时发生区边界存在明显的界面变化特征,这种空间界面的环境要素变化与赤潮发生的时间变化特点是一致的。     

内潮动力机制的数值研究

结合南沙群岛海域出现的内潮现象 ,利用一个两层数值模式 ,通过数值试验来探讨潮汐性质、水体的层结程度、跃层所在的深度、风应力等因子对陆架坡折处内潮特性的影响 ,结果表明 :内潮的波动周期只与海区的潮汐性质有关 ;风应力、非线性效应以及正压潮源地的深度等都可改变内潮的波动振幅 ;内潮的波动振幅随海区的水体层结程度及温跃层所在的深度不同而改变     

广西近海潮汐和海流的观测分析与数值研究——Ⅰ.观测与分析

1992 - 1 996年 ,对广西近海水深浅于 2 0 m的整个海域 ,利用 2 8个锚碇浮标和 2 3个周日连续观测站 ,对水位和海流进行了观测。根据观测结果 ,对该海域的潮汐性质 ,潮流特征和余流分布进行了分析 ,并对潮汐与潮流性质差异的动因和影响余流变化的主要因子进行初步探讨。结果表明 :广西近海的潮汐性质属于正规日潮。潮流性质大部分属于不正规半日潮或不正规日潮。余流主要由潮余流和风海流组成。潮流大小潮变化和风的变化是导致余流变化的主要原因。潮流绕海角运动 ,流速增强 ,并能诱导经向离岸流     

Application of a coupled discontinuous–continuous Galerkin finite element shallow water model to coastal ocean dynamics

A coupled discontinuous–continuous Galerkin (DG–CG) shallow water model is compared to a continuous Galerkin generalized wave-continuity equation (GWCE) based model for the coastal ocean, whereby local mass imbalance typical of GWCE-based solutions is eliminated using the coupled DG–CG approach. Two mass imbalance indicators for the GWCE-based model are presented and analyzed. The indicators motivate discussion on the suitability of using a GWCE-based model versus the locally conservative coupled DG–CG model. Both realistic and idealized test problems for tide, wind, and wave-driven circulation form the basis of the study. For the problems studied, coupled DG–CG solutions retain the robustness of well-documented solutions from GWCE-based models and also capture the dynamics driven by small-scale, highly advective processes which are problematic for GWCE-based models. Issues associated with the coupled DG–CG model are explored, including increased cost due to increased degrees of freedom, the necessary application of slope limiters, as well as the actual coupling process.     

Wave Hindcast for the Neighbouring Seas of Korea Based on Loosely Coupled Wave-Tide-Surge Model

A hindcast simulation of 75 typhoons and winter monsoons which affected the coastal areas of Korean Peninsula is performed by use of a third generation ocean wave prediction model, WAM-cycle 4 model, loosely coupled with a com-bined tide and surge model. Typhoon wind fields are derived from the planetary marine boundary layer model for effective neutral winds embedding the vortical storm wind from the parameterized Rankin vortex type model in the limited areas of the overall modeled region. The hindcasted results illustrate that significant wave heights (SWH) considering the wave-tide-surge coupled process are significantly different from the results via the decoupled case especially in the region of the estuaries of the Changjiang Estuary, The Hangzhou Bay, and the southwestern tip of Korean Peninsula. This extensive model simulation is the first attempt to investigate the strong wave-tide-surge interaction for the shallow depth area along the coasts of the Yellow Sea and the East China Sea Continental     

Detecting red tides in the eastern Seto inland sea with satellite ocean color imagery

A multi-spectral classification scheme is proposed to identify water with red tide(s) using satellite ocean color imagery obtained by the Sea-viewing Wide Field-of-view Sensor (SeaWiFS). The study area was the eastern Seto Inland Sea in Japan, where serious red tides frequently occur. “Background Ocean Colors” (BOCs hereafter), or colors of water around a red tide or those of the water before/after a red tide, are calculated as the monthly climatological average of normalized water-leaving radiances (nLw) with 0.01 degree spatial resolution with SeaWiFS imagery. Criteria for detecting red-tide pixels are established from analyses of characteristics of the nLws (in the 443, 490, 510, and 555 nm bands) anomalies from BOCs and the nLw spectra together with the red-tide records in Osaka Bay. The proposed scheme can efficiently indicate the presence or absence of red tides for independent match-ups with 83% accuracy. Additional validations of specific events indicate that the algorithm performed well in the study area. These results suggest that the scheme is appropriate to detect red tides in the optically complex coastal water of the eastern Seto Inland Sea.     

Multilevel finite-difference model for three-dimensional hydrodynamic circulation

A three-dimensional finite-difference multilevel hydrodynamic model is developed using an explicit scheme on a staggered grid. The model has been tested against four cases, namely (i) wind-induced circulation (ii) density-driven circulation (iii) seiche oscillation in a closed basin and (iv) tide-induced circulation in a open channel. The results obtained in the present study compare well with those obtained from the corresponding analytical solutions under idealised conditions for the above four cases. The model was also tested against the case of circulation induced by wind and Coriolis force and the results obtained are compared with the results of Davies and Owen (1979).     

江苏近海风暴潮增水数值模拟研究

为研究江苏近海海域风暴潮的特性以及为该海域风暴潮增水变化机理及后报做铺垫,本文基于FVCOM(Finite Volume Coast and Ocean Model)海洋模式和Jelesnianski圆形台风风场模型,建立了江苏近海风暴潮数值模型,并对江苏近海的天文潮以及1109号台风和1210号台风引起的风暴潮进行模拟。结合验潮站水位观测,研究了连云港站和吕泗站的天文潮和风暴潮增水过程。我们将风暴潮与天文潮非线性作用下的风暴潮增水和纯风暴潮增水过程进行对比,讨论了天文潮与1109号和1210号台风风暴潮之间的非线性作用引起的增水特征。结果均表明,在天文潮高潮时,天文潮和风暴潮之间的非线性作用可以抑制增水,在天文潮低潮时,天文潮和风暴潮之间的非线性作用有利于增水。除了气象因子以及天文潮和风暴潮之间的非线性作用外,该海区的地理环境也对台风风暴潮增水产生影响。因此对江苏近海的海岸线变化和浅滩地形变化进行敏感性试验,结果表明,本文所设计的海岸线变化对该海域的风暴潮增水影响较小,江苏沿海岸线的向外推移使得江苏海域风暴潮的增水略微上涨,而本文所设计的地形的变化对风暴潮增水影响较大。     

渤海湾水位、流场的一个数值预报模式

本文利用普林斯顿海洋模式(POM),建立了渤、黄海风增水预报模式和渤海湾水位海流预报模式;利用正交曲线网格提高重点区域的分辨率;采用美国海洋大气局(NOAA)全球预报风场和气压场作为模式表面强迫场,将计算域边界上的天文潮预报值与风增水模式预报的余水位相叠加构建模式的边界条件,在正压条件下,模拟了渤海湾2002年的水位流场过程。结果表明,模式能够较好地再现计算域内天文潮和综合水位的预报,域内10个潮位站模式与实测分析的m1和M2分潮的振幅与迟角差均不超过5.1cm和6.3°,15个潮流站模式与实测分析的m1和M2分潮流的振幅与迟角差均不超过7.5cm/s和15.8°,模式预报的水位值与塘沽站实测值非常接近,预报精度较单纯的天文潮预报有明显提高。     

Open boundary conditions for regional tidal simulations

The performance of regional tide model simulations is examined in relation to the choice of open boundary conditions. Three barotropic open boundary conditions, clamped elevation, clamped normal velocity, and Flather, give similar results when the prescribed values are exact; however, Flather is much less sensitive to errors in the prescribed values. Of particular concern, it was found that with a phase error between the two boundaries, both the clamped conditions resulted in magnitude errors in the unclamped variable (although the simulation remained stable).A modified flow relaxation scheme for the depth-varying prognostic variables is presented. This implementation allows the transmission of a range of vertical modes while retaining realistic topography at the boundary. It was found to be an excellent internal tide boundary condition in tests comparing simulations of different domain length encompassing a ridge and sloping bottom, and in a comparison to an analytical solution. Mass is conserved without any artificial volume constraint.     

风暴潮模拟中潮位对风拖曳力系数的影响研究

近年来,应用数值模型模拟台风引起的风暴潮运动越来越普遍,模型中对于风拖曳力系数的确定,一般都从相对风速出发,可引用的公式也较多,但这些公式很少考虑潮位变化对此系数的影响.在强潮河口、海岸海域,潮位变幅大,最高潮位甚至可达风速参考高度(10m)的近一半,如长江口和杭州湾.在数值模拟中不考虑风暴潮和天文潮共同引起的潮位变化,会造成风应力高潮时被低估、低潮时被高估的现象,从而影响风暴潮模拟的精度.为此本文对现有的风拖曳力系数加以改进,提出了考虑潮位影响的风拖曳力系数表达式,并应用于长江口、杭州湾9711号台风风暴潮的模拟中,增水模拟结果得到了明显改善,可进一步推广应用于强潮河口、海岸的风暴潮增水模拟中.     

The accuracy of surface elevations in forward global barotropic and baroclinic tide models

This paper examines the accuracy of surface elevations in a forward global numerical model of 10 tidal constituents. Both one-layer and two-layer simulations are performed. As far as the authors are aware, the two-layer simulations and the simulations in a companion paper (Deep-Sea Research II, 51 (2004) 3043) represent the first published global numerical solutions for baroclinic tides. Self-consistent forward solutions for the global tide are achieved with a convergent iteration procedure for the self-attraction and loading term. Energies are too large, and elevation accuracies are poor, unless substantial abyssal drag is present. Reasonably accurate tidal elevations can be obtained with a spatially uniform bulk drag c_d or horizontal viscosity K_H, but only if these are inordinately large. More plausible schemes concentrate drag over rough topography. The topographic drag scheme used here is based on an exact analytical solution for arbitrary small-amplitude terrain, and supplemented by dimensional analysis to account for drag due to flow-splitting and low-level turbulence as well as that due to breaking of radiating waves. The scheme is augmented by a multiplicative factor tuned to minimize elevation discrepancies with respect to the TOPEX/POSEIDON (T/P)-constrained GOT99.2 model. The multiplicative factor may account for undersampled small spatial scales in bathymetric datasets. An optimally tuned multi-constituent one-layer simulation has an RMS elevation discrepancy of 9.54 cm with respect to GOT99.2, in waters deeper than 1000 m and over latitudes covered by T/P (66N to 66S). The surface elevation discrepancy decreases to 8.90 cm (92 percent of the height variance captured) in the optimally tuned two-layer solution. The improvement in accuracy is not due to the direct surface elevation signature of internal tides, which is of small amplitude, but to a shift in the barotropic tide induced by baroclinicity. Elevations are also more accurate in the two-layer model when pelagic tide gauges are used as the benchmark, and when the T/P-constrained TPXO6.2 model is used as a benchmark in deep waters south of 66S. For Antarctic diurnal tides, the improvement in forward model elevation accuracy with baroclinicity is substantial. The optimal multiplicative factor in the two-layer case is nearly the same as in the one-layer case, against initial expectations that the explicit resolution of low-mode conversion would allow less parameterized drag. In the optimally tuned two-layer M₂ solution, local values of the ratio of temporally averaged squared upper layer speed to squared lower layer speed often exceed 10.     

基于海气耦合模型的渤海两种类型风暴潮数值模拟与对比分析

渤海一年四季都易受到由温带风暴和热带气旋所致风暴潮的影响。为了缓解风暴潮灾害对海岸地区人员生命财产的影响,十分有必要了解大型风暴潮的发生过程和机制。目前大部分研究主要局限于单一的温带风暴潮或台风风暴潮。本文利用所构建的海气耦合数值模型研究了发生于渤海的两种类型的风暴潮,对发生在渤海的2次典型强风暴潮过程进行了模拟。由WRF模型模拟得到的风场强度和最低海平面气压与实测数据吻合较好,由ROMS模型模拟得到的风暴潮期间水位变化过程与潮位站观测结果也吻合较好。对两种类型风暴潮期间的风场结钩、海面风应力、海洋表面平均流场以及水位分布进行了分析对比,并将耦合模型结果与非耦合模型结果进行了对比。研究表明,渤海两种类型风暴潮期间的风场结钩、海面风应力、海洋表面平均流场以及水位分布等均存在巨大差异。渤海风暴潮的强度主要由海洋表面的驱动力所决定,但同时也受海岸地形地貌的影响。     

Implementation of method of lines to predict water levels due to a storm along the coastal region of Bangladesh

In this study, the method of lines (MOL) has been applied to solve two-dimensional vertically integrated shallow water equations in Cartesian coordinates for the prediction of water levels due to a storm surge along the coast of Bangladesh. In doing so, the partial derivatives with respect to the space variables were discretized by the finite difference (central) method to obtain a system of ordinary differential equations (ODEs) with time as independent variable. The classical fourth-order Runge–Kutta method was used to solve the obtained system of the ODEs. We used a nested finite difference scheme, where a high resolution fine grid model (FGM) capable of incorporating all major islands along the coastal region of Bangladesh was nested into a coarse grid model (CGM) covering up to 15°N latitude of the Bay of Bengal. The boundaries of the coast and islands were approximated through proper stair step. Appropriate tidal condition over the model domain was generated by forcing the sea level to be oscillatory with the constituent M ₂ along the southern open boundary of the CGM omitting wind stress. Along the northeast corner of the FGM, the Meghna River discharge was taken into account. The developed model was applied to estimate water levels along the coast of Bangladesh due to the interaction of tide and surge associated with the April 1991 storm. We also computed our results employing the standard finite difference method (FDM). Simulated results show the MOL performs well in comparison with the FDM with regard to CPU time and stability, and ensures conformity with observations.     

海上中距离跨海高程传递方法比较

远离岸线的中距离(30~50 km)海域缺少测绘基础资料,高程传递尚未形成相应的作业体系。为满足海上工程建设的精度要求,采用同步验潮平均海面传递法、GNSS/水准传递法以及DTU18全球平均海面模型传递法,在山东省半封闭海域与开阔海域进行了中距离海上高程传递实验。实验结果表明,15天同步验潮平均海面传递法与GNSS/水准法的精度相当(厘米级),与DTU18全球平均海面模型传递法进行直接传递的差异在15 cm内,可为不同测绘基础条件下海上风电高程基准传递提供作业参考。     

龙口市人工岛周边海域波浪场和水动力场的数值模拟

为了科学的分析离岸人工岛群对周边海岸动力的影响,本文以龙口市人工岛为例,运用MIKE21数值模拟软件建立了龙口离岸人工岛及附近海域的水动力模型和波浪模型,根据波浪的生成、成长和传播原理,针对最有可能形成较大波浪的W、WNW、NW三个方向,采用频率为0.5%(重现期为200a)的高潮位叠加频率为0.5%(重现期为200a)的W、WNW和NW向风作用下的波浪场和水动力场进行数值模拟研究,给出了高潮位条件下人工岛周边不同波向对应的有效波高和最大波高值,以及最大流速和流向数据,可为人工岛群的地坪标高确立和防浪建筑物的布设提供科学依据。     

基于MIKE SA模型的厦门西港海域溢油影响的数值模拟研究

基于MIKE SA溢油模块,以燃料油为油种,建立了厦门西港海域溢油模型,模拟静风、主导风向（东北东风）和不利风向（西南风）3种风场条件下,一个潮周期内涨急、高潮、落急和低潮4个时段发生10 t溢油后油膜的漂移路径和影响范围.结果显示,发生在厦门西港海域的溢油在海面的漂移过程主要受潮流和风的影响,其中潮流起着主导作用.不同风向条件下,24 h内油膜的影响范围不同,静风条件下溢油浓度超一类（或二类,≥0.05 mg/dm3）、超三类（≥0.30 mg/dm3）和超四类（≥0.50 mg/dm3）的总影响面积分别为31.33、19.63和11.74 km2;主导风向条件下溢油浓度超一类（或二类）、超三类和超四类的总影响面积分别为99.62、69.01和8.99 km2;不利风向溢油浓度超一类（或二类）、超三类和超四类的总影响面积分别为8.38、5.05和2.10 km2.该预测结果可给出溢油事故发生后的影响范围、影响程度和影响敏感目标的时间,可为溢油事故应急决策的制定及溢油损害评估提供科学决策和支持,提升厦门海域环境风险管理应急能力建设.     

使用间隔气象数据的余水位分析

余水位是潮汐预报中的"报不准"部分。在现有条件下,许多验潮站在对水位数据进行连续采集的同时,对气象数据还使用间隔时间采集的方式(若干小时采集一次)。使用间隔气象数据,按风向不同对气象数据、余水位数据进行分组,对风与余水位关系进行定性分析,对风-压与余水位关系进行拟合分析,并结合实例,对方法进行说明和检视。