Applications of the Mild-Slope Equation to Tidal Computations in the Taiwan Strait

Employing harmonic analysis of tidal data in the Taiwan Strait, the cross-strait tidal characteristics are completely illustrated. Based on the two dimensional mild-slope equation which can be reduced to the shallow-water wave equation, a finite element model (Tsay et al., 1989) is applied to investigate the characteristics of tides in the Taiwan Strait. The co-range and equi-phase charts of major tidal constituents, such as M₂, S₂, N₂, and K₁, are reproduced. Anomalous amplification of semidiurnal tides in the Taiwan Strait is verified. With rotation effects neglected and by applying a non-reflective condition on the open boundaries, the numerical results of phase-lag and co-range distributions show very good agreement with observed data for semidiurnal tides in the Taiwan Strait. Due to crude representation of the topography at two ends along the China coast, computed tidal distributions deviate from the observations. However, both computed amplitudes and phase-lags compare very well with observed data along the central half of the China coast.     

浙江近海潮汐潮流的数值模拟

用三维陆架海模式(HAMSOM)对浙江近海的潮汐、潮流进行了数值模拟,并采用网格嵌套和动边界技术对原模式作了改进,以提高计算的精度,改进后的模式在浙江近海的应用中被证明是成功的.沿岸50个潮位站计算与实测值的比较表明,加入动边界以后的小区域细网格计算较之粗网格以及未加动边界以前精度普遍提高,比较的均方差结果为:M₂分潮振幅差4.6cm,相角差7.14°;S₂分潮振幅差5.0cm,相角差5.4°;K₁分潮振幅差2.25cm,相角差5.76°;O₁分潮振幅差1.56cm,相角差5.5°,可见计算与实测符合良好.另外,选取了105个实测潮流点,比较了表层M₂和K₁分潮流调和常数分量Ucosξ,Usinξ,Vcosη,Vsinη的实测值与计算值的偏差,结果表明计算与实测的符合程度较好.在此基础上,给出了各主要分潮的潮位同潮图、潮流同潮图、潮汐性质、潮流性质、潮流椭圆和潮流的运动形式等,发现4个主要分潮M₂,S₂,K₁,O₁在本区内均未出现无潮点;M₂分潮流在29°18'N,122°46'E处有一个圆流点.此外还得到了一些有意义的结论,都与实测情况符合良好,从而对整个浙江沿海区域的潮汐潮流特性有了一个全面认识.     

Effect of Sea Level Variation on Tidal Characteristic Values for the East China Sea

1 .IntroductionTheglobalairtemperatureroseabout 0 .5～ 0 .6°Coverthepast 2 0thcentury ,andtheglobalmeansealevelincreasedbyabout2 0cmduringtheperiod .Theregionalmeansealevelriseswiththerisingglobalmeansealevel.Zuoetal.( 1 997)indicatedthatthemeanrisingrateofabsolutemeansealevelalongtheChinacoastontheassumptionofunifiedisostaticdatumis 2mm a .Woodworth( 1 999)analyzedsealevelspanning 1 76 8tothepresentinLiverpool,andobtainedaseculartrendforheperiodupto 1 880of0 .39± 0 .1 7mm a ,andatrendfort…     

Wavelet and artificial neural network analyses of tide forecasting and supplement of tides around Taiwan and South China Sea

In multi-resolution analysis (MRA) by wavelet function Daubechies (db), we decompose the signal in two parts, the low and high-frequency contents. We remove the high-frequency content and reconstruct a new “de-noise” signal by using inverse wavelet transform. The calculation of tidal constituent phase-lags was made to determine the input and output data patterns used in building network structure of Artificial Neuron-Network (ANN) model. The “de-noise” signal was, then, used as the input data to improve the forecasting accuracy of the ANN model. The wavelet spectrum, conventional energy spectrum (fast Fourier transform, FFT), and harmonic analysis were used to analyze the characteristics of tidal data.Using only a very short-period data as a training data set in Artificial Neuron-Network Back-Propagate (ANN-BP) model, the developed ANN+Wavelet model can accurately predict or supply the missing tide data for a long period (1–5 years). The results also show that the concept of tidal constituent phase-lags can improve ANN model of tidal forecasting and data supplement. The addition of the wavelet analysis to ANN method can prominently improve the prediction quality.     

A barotropic model of the linear semidiurnal tide over a continental slope

Semidiurnal tides, and especially the lunar tide M₂, are dominant dynamics in the Bay of Biscay. Strong tidal currents are associated with the presence of a significant continental slope. By combining Newton's gravitation laws and Euler's equations, Laplace's equations contain the astronomical forcing responsible for the observed semidiurnal tides. In shallow waters, this direct forcing is often neglected. We study here its influence on the tidal dynamics over the continental slope through the development of a simple model describing the barotropic semidiurnal dynamics on a transect perpendicular to the slope. This new model results from the combination of two different models, i.e. the one developed by Rosenfeld and Beardsley (1987), which takes into account the tide-generating force, and that of Battisti and Clarke (1982), which neglects it. A first model is developed by neglecting the direct astronomical forcing in equations: it consists in solving a second-order homogeneous propagation equation for the barotropic semidiurnal tide and needs only coastal conditions as well as the knowledge of the along-slope wave number of the solution. For a mean slope typical of the South Brittany area, this non-forced model provides results in accordance with those of Battisti and Clarke and Le Cann (1990): in particular, in the upper part of the slope, it shows a polarization inversion of tidal ellipses characteristic of the tidal dynamics observed in this area. Then, the direct astronomical forcing is kept in equations. The simple model developed without this forcing is fitted in order to solve the resulting forced propagation equation for the barotropic tide. The solution of this second model is the sum of a forced wave responding to the direct astronomical forcing and of a free wave generated at the coastal boundary. Under the same boundary conditions, the results obtained with the influence of the tide-generating force are then compared with those obtained without it. This comparison allows one to apprehend the importance of the direct astronomical forcing on tidal dynamics across the slope: in particular, the main difference appears in deep waters where this forcing induces a phase-lag between the plain and the shelf for the sea-surface slope.     

Tide Forecasting of Tides Around Taiwan by Artificial Neural Network Method and Wavelet Analysis

In multiresolution analysis(MRA)by wavelet function Daubechies(db),we decompose the signal to two parts,the low and high frequency content.The high-frequency content of the data is removed first and a new "de-noise" signal is reconstructed by using inverse wavelet transform.The wavelet spectrum and harmonic analysis were used to analyze the characteristics of tidal data before constructing the input and output structure of ANN model.That is,the concept of tidal constituent phase-lags was introduced and the new "de-noise" signal was used as the input data set of ANN and the forecasting accuracy of ANN model is significantly improved.     

时差法水位改正的精度评估方法研究

为实现对时差法水位改正结果的精度评估,根据协方差传播律,推导了时差法水位改正的误差方程,讨论了评估潮时差确定精度的方法,并通过数据试验对结果进行了验证。试验结果表明,最小二乘拟合法比相关系数法更能满足时差法水位改正精度评估对潮时差的要求;根据《海道测量规范》要求,试验海域潮时差中误差的限差为6.2min,验潮站网潮时差闭合差中误差的限差为1.0min;水位改正方差反映出了水位改正结果的合理性,可用以评估水位改正精度。     

A validation concept for cohesive sediment transport model and application on Lianyungang Harbor,China

Morphology evolution and hydro-sedimentological interactions in muddy coastal environments are long term processes. These processes are closely related to suspended sediment transport driven by waves and tidal currents. In the traditional calibration/verification methodology for cohesive sediment transport models, time-series data of suspended sediment concentration (SSC) measured during one or several tidal periods are commonly used as major validation references. But the disadvantage of this approach includes that it cannot filter the noises caused by the stochastic nature of short term hydrodynamics induced by waves and the varying properties of bottom sediments; besides a phase-lag phenomenon is often observed between sediment transport and hydrodynamics in the short term. On the contrary, a stable relationship between sedimentation and hydrodynamics is usually found in the long term. For a specific weather, the SSC values often agree well with local wave heights for muddy coasts. Therefore, in this paper a conceptual quantity defined as “representative SSC” was brought forward, and a model validation concept, including calibration and verification approaches, was proposed, in which calibration is performed against yearly-representative SSC values and the long term transport trend, and verification by using representative SSC values corresponding to different wave conditions. A numerical simulation was set up, and a real-life engineering application, Lianyungang Harbor, China, was executed to elaborate the proposed validation concept. Finally, the characteristics of SSC distribution around Lianyungang Harbor were discussed.     

Tidal data inversion: interpolation and inference

Initial efforts in applying inverse methods to studies of ocean tides have focused on making the best use of a small number of observations to map tidal fields in a large area. As such, inversion can be viewed as an objective analysis scheme which uses a dynamically appropriate spatial covariance, derived from the shallow water equations, to interpolate and smooth a sparse data set. Data from recent altimetry missions are not sparsely distributed relative to tidal wavelengths in the open ocean, apparently reducing the need for complicated dynamically based interpolation schemes. Altimetric data sets are also quite large, making application of rigorous inversion methods to global tidal modeling a challenging computational problem. We describe here a new iterative solution scheme which allows us for the first time to fit the full set of TOPEX/Poseidon cross-over differences. The resulting solution (TPXO.3) fits validation tide gauges significantly better than previous inverse solutions. TPXO.3 also reduces residual cross-over variances relative to other recent inverse and empirical solutions, particularly in shallow water where improvements are dramatic. With the new solution approach very significant improvements in global tidal models should be possible in shallow areas and in the vicinity of complex bathymetry, where high-accuracy tidal modeling remains a challenging problem. With the recent improvements in the definition of tidal elevations in the open ocean it should now also be possible to resolve some long unanswered questions about tidal energetics and dynamics. Inverse methods provide a natural framework for addressing these issues, and making inferences about tidal dynamics. In particular, by bringing data and dynamics together in a single solution, we can rigorously test the consistency of the two. We present results of global and local inversions which suggest that over elongated bathymetric features oriented perpendicular to tidal flows, energy dissipation in the open ocean is significantly enhanced, presumably due to conversion of barotropic tidal motions into baroclinic modes. For the M₂ tide our preliminary results suggest that perhaps as much as 0.5 TW of energy is dissipated in this manner. However, due to the simplified linear dynamics and limited spatial resolution used for our inversion, there are significant uncertainties associated with these results. A more careful application of inverse methods to make more rigorous inferences about tidal energetics, including use of more reasonable prior dynamics, and the highest possible spatial resolution, should allow for closure of the tidal energy budget within the next few years.     

渤海、黄海、东海M2潮汐潮流的三维数值模拟

利用建立的一种新的半隐半显三维数值格式,将渤海、黄海、东海作为一个整体,采用球面坐标系下的三维潮波方程组,考虑了引潮力的作用,数值模拟了渤海、黄海、东海的M₂分潮的潮汐与潮流,结果较好地体现了渤海、黄海、东海M₂分潮的特征.通过比较65个验潮站的实测值与计算值,所得计算结果的振幅差平均为6.4cm,相角差为6.1°,计算与实测符合良好.本文给出的问潮图与Fang于1986年给出的实测占数值综合结果基本一致.对选取的47个测流站,比较了各层潮流调和常数Ucosζ、Usinζ、Vcosη、Vsinη的计算值与实测值的偏差,偏差绝对值的平均在2.6～4.9cm/s之间.并比较分析了潮流的垂直结构,所得结果与实测符合较好.首次揭示出回流点的水平位置不随深度变化这一特性.最后给出了M₂分潮的潮能消耗.     

北印度洋半日潮波的数值模拟研究

基于FVCOM(Finite Volume Coast and Ocean Model)模型,建立北印度洋海域(31^°~102^°E,16^°S~31^°N)的M₂和S₂分潮潮波数值模式,研究北印度洋半日潮潮汐、潮流分布特征。对底摩擦系数进行数值试验,利用代价函数梯度下降法,得到分潮调和常数向量均方根偏差(RMSE)的变化曲线,逼近并确定最优的底摩擦系数。将采用该系数的模拟结果与TOPEX/Poseidon卫星高度计交叉点的调和常数数据、国际海道测量组织(IHO)及部分文献中的验潮站数据进行比较与验证,一致性较好。其中对比卫星数据的振幅偏差为2~4 cm、迟角偏差为7^°~8^°,与验潮站数据的振幅偏差为3~6 cm、迟角偏差为8^°~9^°。根据模拟结果,分析了北印度洋海域M₂和S₂分潮潮波传播特征和潮流椭圆的空间分布特征等。M₂分潮潮波在阿拉伯海南部有1个无潮点,在波斯湾内有2个无潮点,最大振幅超过80 cm;潮流在西北印度洋和孟加拉湾中部大多为顺时针旋转,其余海域大多为逆时针旋转;流速在阿拉伯海东北部、安达曼海、波斯湾和孟加拉湾北部较大,最大流速为160 cm/s,其他海域较小。S₂分潮的潮波传播特征、无潮点的位置和潮流椭圆的空间分布特征等都与M₂分潮类似,但潮波振幅和潮流流速等都相对M₂分潮较小。研究完善了北印度洋海域2个主要半日分潮M₂和S₂的整体特征。     

Silting in the Lower Courses of Tidal Sluices in China

Serious sediment deposition often occurs after the construction of tidal sluices in snrall or medium-sized tidal muddy estuaries, so desilting or dredging is needed to meet the demands of flood diseharge, saltwater retaining, and navigation in those areas. In this paper, the problem of sediment deposition induced by construction of tidal sluices is analyzed.Different problems of silting near tidal sluices for different types of estuaries are summed up. at the same time, corresponding methods are given to solve these probleras, and a few successful examples are also given. The idea of comprehensive regulation and utilization of estuaries is put forward, and some proposals for solution of sediment deposition in this kind of estuaries are made.     

西北太平洋的一种潮汐数值同化模型

利用FVCOM海洋数值模式,在球坐标系统下考虑非线性效应和天体引潮力的影响,基于非结构的三角形网格建立了包括中国近海、日本海、鄂霍次科海和部分西北太平洋海域的高分辨率海洋潮汐数值模型,并采用趋近法同化84个沿岸验潮站的观测资料。模拟结果与175个验潮站的实测结果拟合良好,M₂,S₂,K₁,O₁四个主要分潮振幅和迟角的绝对平均误差分别为4.0 cm和5.6°,2.4 cm和7.5°,2.6 cm和6.3°,1.5 cm和5.0°。依据调和分析结果给出了4个主要分潮的同潮图分布,得到8个半日分潮和5个全日分潮的无潮点,证实了宗谷海峡全日潮无潮点的存在,首次模拟得到津轻海峡的全日潮无潮点;还给出了整个计算海域内最大可能潮差和潮汐余水位的分布特征。     

印度尼西亚海域潮波的数值研究

基于ROMS模式构建了模拟区域为(15.52°S-7.13°N,110.39°~134.15°E)水平分辨率为2′的潮波数值模式,分别模拟了印尼海域M2、S2、K1、O1四个主要分潮。模拟结果与29个卫星高度计交叠点上的调和常数进行比较,符合较好。M2分潮的振幅均方根差为3.4cm,迟角均方根差为5.9°;S2分潮的振幅均方根差为1.7cm,迟角均方根差为6.3°;K1分潮振幅均方根差为1.1cm,迟角均方根差为5.8°;O1分潮振幅均方根差为1.2cm,迟角均方根差为4.4°。M2、S2、K1、O1分潮向量均方根差分别为3.8cm、2.4cm、1.9cm和1.3cm,模拟结果的相对偏差在10%左右。根据计算结果分析了印尼海域的潮汐特征及潮能传播规律,结果显示:爪哇海以外的印尼海域主要为不规则半日潮区;全日潮潮能主要由太平洋传入印尼海域,而半日潮潮能则是从印度洋传入印尼海域。     

涠洲岛附近海域溢油数值模拟

采用潮流控制方程和湍封闭方程构建完整的三维浅海潮波定解方程组,建立涠洲岛附近海域的三维潮流数值模型来计算涠洲岛附近海域的潮流和潮位变化情况,采用大海域计算得出的调和常数值作为开边界的输入值,从潮位验证、潮流验证和流场的变化规律等多方面情况看,计算结果与实际监测符合良好,具有较高的精确度,模拟计算结果较真实地反映了涠洲岛附近海域的潮位变化和潮流运动状况.在潮流模拟验证正确的情况下,建立溢油预测模型,采用欧拉-拉格朗日追踪方法,进行油膜中心轨迹的预测,并预测了油膜漂移的扩展面积、扫海面积和残留量随时间的变化值.     

基于FVCOM 的渤海潮波数值模拟

基于有限体积法海洋数值模型(FVCOM),对渤海当前水深岸线状况下的潮汐潮流进行了数值计算。模式采用不规则三角形网格,较好地提高了黄河口处网格分辨率,模拟了渤海海域K1,O1,M2和S2四个主要分潮。利用渤海沿岸19个验潮站的资料对模拟结果进行了验证,K1分潮振幅绝均差2.39 cm,迟角绝均差4.36°,O1分潮振幅绝均差1.40 cm,迟角绝均差4.29°,M2分潮振幅绝均差为3.55 cm,迟角绝均差为5.69°,S2分潮振幅绝均差1.72 cm,迟角绝均差8.86°,结果显示各分潮模拟结果合理,较真实地反映了渤海海域四个分潮传播情况。     

珠江河口潮能收支研究

Tidal energy budget in the Zhujiang(Pearl River) Estuary(ZE) is evaluated by employing high-resolution baroclinic regional ocean modeling system(ROMS). The results obtained via applying the least square method on the model elevations are compared against the tidal harmonic constants at 18 tide stations along the ZE and its adjacent coast. The mean absolute errors between the simulation and the observation of M_2, S_2, K_1 and O_1 are 4.6, 2.8, 3.2 and 2.8 cm in amplitudes and 9.8°, 15.0°, 4.6° and 4.6° in phase-lags, respectively. The comparisons between the simulated and observed sea level heights at 11 tide gauge stations also suggest good model performance. The total tidal energy flux incoming the ZE is estimated to be 343.49 MW in the dry season and larger than 336.18 MW in the wet season, which should due to higher mean sea level height and heavier density in the dry season. M_2, K_1, S_2, O_1 and N_2, the top five barotropic tidal energy flux contributors for the ZE,import 242.23(236.79), 52.97(52.08), 24.49(23.96), 16.22(15.91) and 7.10(6.97) MW energy flux into the ZE in dry(wet) season, successively and respectively. The enhanced turbulent mixing induced by eddies around isolated islands and sharp headlands dominated by bottom friction, interaction between tidal currents and sill topography or constricted narrow waterways together account for the five energy dissipation hotspots, which add up to about 38% of the total energy dissipation inside the ZE.     

基于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,三者在水平方向基本呈北进南出的分布形态。     

A method of accounting for tidal changes in regional climates of a water basin under conditions of an ice-free Barents Sea

To account for tidal variations in the regional climate of a water basin, we propose adding up the vertical eddy diffusivity, determined by wind and thermohaline forcings, and the diapycnal diffusivity, determined from the solution to the problem of the internal tidal wave (ITW) dynamics. This approach agrees with the approximation of “weak interaction” between turbulence of various origins. Then, the hydrothermodynamics equations are integrated with and without regard for ITW-induced diapycnal diffusion until a quasistationary solution is reached. Next we compare these solutions, found by using the 3D finite-element hydrostatic model QUODDY-4. This comparison shows that the contribution of tides to the formation of the Barents Sea climate in summer is not negligible with respect to certain hydrological characteristics. We present the fields of the dynamic topography of a free surface, surface current velocities, and seawater temperature and salinity at the depth of the pycnocline in the sea to illustrate the occurrence of tidal effects.