基于POM模式与blending同化法建立中国近海潮汐模型

利用POM海洋数值模式建立了中国近海（2°N-41°N,99°E～132°E）分辨率为5′×5′的潮汐模型,模式采用blending同化法同化了由10年TOPEX／Poseidon测高数据反演的潮汐参数与沿岸52个验潮站观测。精度分析表明建立的潮汐模型的8分潮RSS为12．5cm。     

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

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

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了4个主要分潮(M2,S2,K1,O1)和两个浅水分潮(M4,MS4)。分析表明,在渤黄海潮波系统数值模拟中,稳定后选取14 d的数值模拟结果进行调和分析能够取得最佳(最合理)的调和分析结果。计算出调和常数的模拟值与实测值之差的绝对平均值:M2分潮的振幅差为4cm,迟角差为3.3°,S2分潮的振幅差为2cm,迟角差为4.2°,K1 分潮的振幅差为1cm,迟角差为3.7°,O1分潮的振幅差为2 cm,迟角差为5.5°。实验结果较好地体现了渤黄海潮波系统的特征。     

Ocean Tides Observed from A GPS Receiver on Floating Sea Ice Near Chinese Zhongshan Station,Antarctica

Due to limit of coverage in TOPEX/Poseidon (T/P) satellite and sparseness of in-situ tide gauges around Antarctica, the accuracy of global ocean tide models in Antarctic seas is relatively poorer than in low- and mid-latitude regions. To better understand ocean tides in Prydz Bay, east Antarctica, a GPS receiver was deployed on floating sea ice to measure tide-induced ice motion in multiple campaigns. Four online Precise Point Positioning (PPP) services are used to process the GPS data in the kinematic PPP mode, and UTide software is used to separate the major tidal constituents. Comparison between results from different processing methods (relative processing solutions from Track, kinematic PPP solutions from online services) and with bottom pressure gauge (BPG) shows that, high-accuracy tidal information can be obtained from GPS observations on floating sea ice, the root-sum-square (RSS) for the eight major constituents (O1, K1, P1, Q1, M2, S2, N2, K2) is below 4 cm. We have also studied the impacts of data span and filter edge effects at daily boundaries on the accuracy of tide estimates, and found that to obtain reliable tide estimates and neglect the filter edge effects, continuous observation longer than 30 days is necessary. Our study suggests that GPS provides an independent method to estimate tides in Prydz Bay, and can be an alternative to tidal gauges, which are costly and hard to maintain in Antarctica.     

A Kalman Filter Approach to Realize the Lowest Astronomical Tide Surface

In this paper , we present a novel Kalman filter approach to combine a hydrodynamic model-derived lowest astronomical tide (LAT) surface with tide gauge record-derived LAT values. In the approach, tidal water levels are assimilated into the model. As such, the combination is guided by the model physics. When validating the obtained “Kalman-filtered LAT realization” at all tide gauges, we obtained an overall root-mean-square (RMS) difference of 15.1 cm. At the tide gauges not used in the data assimilation, the RMS is 17.9 cm. We found that the assimilation reduces the overall RMS difference by ～ 31% and ～ 22%, respectively. In the Dutch North Sea and Wadden Sea, the RMS differences are 6.6 and 14.8 cm (all tide gauges), respectively. Furthermore, we address the problem of LAT realization in intertidal waters where LAT is not defined. We propose to replace LAT by pseudo-LAT, which we suggest to realize similarly as LAT except that all water level boundary conditions and assimilated tidal water levels have to be enlarged by a constant value that is removed afterward. Using this approach, we obtained a smooth reference surface for the Dutch Wadden Sea that fits LAT at the North Sea boundary within a few centimeters.     

Tidal current fluctuations in the Soya Current

Long-term current measurements were carried out near the Soya Strait in the Okhotsk Sea during a period from February 1980 to September 1982. The data were divided into five segments, each being 150 days long, and the tidal ellipse parameters of major axis, minor axis, orientation, and phase for the four major constituents (M2, S2, K1 and O1 tides) were calculated at each segment. The major axis of the mean tidal ellipse averaged over five segments was 29.9 cm sec^–1 for O1 tide, 28.3 cm sec^–1 for K1 tide, 10.4 cm sec^–1 for M2 tide, and 3.7 cm sec^–1 for S2 tide. The phase and orientation of the tidal ellipse were much stable. But, the root mean square deviations of the major axis reached 20% of the mean values for all four constituents. The tidal currents estimated from the sea level records at Wakkanai and Esashi along the Hokkaido coast in the Okhotsk Sea show that their amplitudes and phases are in good agreement with the observed ones for all four constituents.     

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

海道测量中几种水位改正方法的比较与分析

针对多波束测深易出现的因水位改正不完善导致的相邻测深条带间的拼接断层,分别采用天文潮预报、基于余水位配置的海洋潮汐推算以及基于日平均海面订正的海洋潮汐推算等方法进行水深测量水位改正。结果表明,后两种方法均适用于多波束水深测量水位改正。     

永暑礁海区潮汐海流特征分析

将永暑礁水域1992年逐时潮位和1993年5月-2002年5月各航次的海流及2002年海面逐时风等观测资料进行统计和分析,讨论该岛礁水域潮汐和海流变化特征并对综合极值海流进行初步估算。结果表明,该岛礁水域全日分潮大于半日分潮,潮汐性质属不正规日潮,最大可能潮差较大。春季该岛礁水域实测海流变化较大;最大流速75cm／s。潮流性质为不正规日潮流,随水深的增加日潮流的性质有所减弱;主要日分潮流多作顺时针方向旋转,主要半日潮流则作反时针方向旋转。春季余流较稳定,速度在10cm／s左右,流向多为偏S-SW向。初步探讨了表层不同方位上综合极值海流的分布,不同方位上综合极值海流相差较大,NE、SSW方位上综合极值海流占份量较显著,显示出季风海流的特征。     

The Adriatic Sea M2 and K1 tides by 3D model and data assimilation

In the present work we explore the impact of assimilating local tide-gauge and altimetric data on the quality of predicting the major Adriatic tides (M₂ and K₁). To that end we compute optimal tidal open boundary conditions for a 3D high-resolution finite-element model by using an incremental assimilation formalism. The essence of the method is the use of two dynamical models where the solution in the complex 3D high-resolution model is sought via assimilation of prediction errors into the simpler 2D model with explicit inverse. In the central numerical experiment, harmonic constants from 12 tide gauges are assimilated and the results are analysed at 31 locations, hence 19 independent ones. The data assimilation contributes to the reduction of maximum amplitude error from 5.6 to 0.5 cm for M₂ and from 3.9 to 0.1 cm for K₁. The assimilation procedure is repeated by assimilating suitably processed Topex/Poseidon altimeter data, again validating the outcome at 31 tide gauge locations. The result was very similar to the gauge-data assimilation outcome. The model output is also validated with the current data, not used in the assimilation. At two locations and at three depths the model was able to reproduce the major and the minor semi-axes of tidal ellipses, as well as their orientations very well.     

Data Assimilation Modeling of the Barotropic Tides in the Korea/Tsushima Strait

During 1999–2000, 13 bottom mounted acoustic Doppler current profilers (ADCPs) and 12 wave/tide gauges were deployed along two lines across the Korea/Tsushima Strait, providing long-term measurements of currents and bottom pressure. Tidally analyzed velocity and pressure data from the moorings are used in conjunction with other moored ADCPs, coastal tide gauge measurements, and altimeter measurements in a linear barotropic data assimilation model. The model fits the vertically averaged data to the linear shallow water equations in a least-squares sense by only adjusting the incoming gravity waves along the boundaries. Model predictions are made for the O₁, P₁, K₁, μ₂, N₂, M₂, S₂, and K₂ tides. An extensive analysis of the accuracy of the M₂ surface-height predictions suggests that for broad regions near the mooring lines and in the Jeju Strait the amplitude prediction errors are less than 0.5 cm. Elsewhere, the analysis suggests that errors range from 1 to 4 cm with the exception of small regions where the tides are not well determined by the dataset. The errors in the model predictions are primarily caused by bias error in the model’s physics, numerics, and/or parameterization as opposed to random errors in the observational data. In the model predictions, the highest ranges in sea level height occur for tidal constituents M₂, S₂, K₁, O₁, and N₂, with the highest magnitudes of tidal velocities occurring for M₂, K₁, S₂, and O₁. The tides exhibit a complex structure in which diurnal constituents have higher currents relative to their sea level height ranges than semi-diurnal constituents.     

基于光学遥感的海岛潮间带和湿地信息提取--以东沙岛(礁)为例

基于高分辨率遥感图像,采用相似三角形原理,结合海岛多年潮汐数据进行了海岛潮间带的确定;利用克里格（Kriging）插值方法,对水深图中海岛浅水区域的水深点数据进行单元网格化处理,在遥感图像的辅助下,进行了海岛湿地范围的界定。在此基础上,确定了东沙岛潮间带和湿地的范围,提取了各自的面积,并对结果进行了分析。     

Accuracy assessment of global ocean tide models in the South China Sea using satellite altimeter and tide gauge data

In this study, to meet the need for accurate tidal prediction, the accuracy of global ocean tide models was assessed in the South China Sea (0°–26°N, 99°–121°E). Seven tide models, namely, DTU10, EOT11a, FES2014, GOT4.8, HAMTIDE12, OSU12 and TPXO8, were considered. The accuracy of eight major tidal constituents (i.e., Q₁, O₁, P₁, K₁, N₂, M₂, S₂ and K₂) were assessed for the shallow water and coastal areas based on the tidal constants derived from multi-mission satellite altimetry (TOPEX and Jason series) and tide gauge observations. The root mean square values of each constituent between satellite-derived tidal constants and tide models were found in the range of 0.72–1.90 cm in the deep ocean (depth>200 m) and 1.18–5.63 cm in shallow water area (depth<200 m). Large inter-model discrepancies were noted in the Strait of Malacca and the Taiwan Strait, which could be attributable to the complicated hydrodynamic systems and the paucity of high-quality satellite altimetry data. In coastal regions, an accuracy performance was investigated using tidal results from 37 tide gauge stations. The root sum square values were in the range of 9.35–19.11 cm, with the FES2014 model exhibiting slightly superior performance.     

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

利用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个全日分潮的无潮点,证实了宗谷海峡全日潮无潮点的存在,首次模拟得到津轻海峡的全日潮无潮点;还给出了整个计算海域内最大可能潮差和潮汐余水位的分布特征。     

基于FVCOM的钦州湾三维潮流数值模拟

基于采用不规则三角网格和有限体积方法的FVCOM模式,建立钦州湾三维潮流数值模型来重现钦州湾的潮位和潮流变化状况.根据模拟结果计算得到了较以往更为精细的钦州湾K1、O1、M2、S2分潮的同潮图,潮汐最大振幅分别为112、96、50和15cm;最大可能流速的分布基本与等深线一致,龙门港附近最大可能潮流流速可达200cm·s-1;钦州湾的外湾口海域开阔,一般为旋转流,近岸海区及水道、河口等多为往复流,K1和M2分潮流椭圆长轴的分布与地形密切相关,旋转方向均为顺时针,流速极值出现在龙门港区.由潮余流场的分布特点可以看出,钦州湾西槽和中槽的西侧区域是其主要的出水通道,东槽和中槽的东侧区域则是主要的进水通道.     

Three-dimensional numerical simulation for tide and tidal current in the Beibu Gulf

ＩＮＴＲＯＤＵＣＴＩＯＮＩｎｅａｒｌｙ 1 96 0’ｓ,ｔｈｅｔｉｄｅａｎｄｔｉｄａｌｃｕｒｒｅｎｔｉｎｔｈｅＢｅｉｂｕＧｕｌｆｗｅｒｅｏｂｓｅｒｖｅｄａｎｄａｎａｌｙｓｅｄｂｙＣｈｉｎａｉｎｃｏｏｐｅｒａｔｉｏｎｗｉｔｈＶｉｅｔｎａｍ1) .ＴｈｅｓｙｓｔｅｍａｔｉｃｓｔｕｄｉｅｓｏｆｔｉｄｅａｎｄｔｉｄａｌｃｕｒｒｅｎｔｉｎｔｈｅＢｅｉｂｕＧｕｌｆｗｅｒｅｆｉｒｓｔｃａｒｒｉｅｄｏｕｔｂｙＦａｎｇ (1 986 ) .Ｔｈｅｈｉｓｔｏｒｙｏｆｎｕｍｅｒｉｃａｌｓｔｕｄｙｏｆｔｉｄｅａｎｄｔｉｄａｌｃｕｒｒｅｎｔ…     

Design, Optimization and Numerical Modelling of A Novel Floating Pendulum Wave Energy Converter with Tide Adaptation

A novel floating pendulum wave energy converter (WEC) with the ability of tide adaptation is designed and presented in this paper. Aiming to a high efficiency, the buoy''s hydrodynamic shape is optimized by enumeration and comparison. Furthermore, in order to keep the buoy''s well-designed leading edge always facing the incoming wave straightly, a novel transmission mechanism is then adopted, which is called the tidal adaptation mechanism in this paper. Time domain numerical models of a floating pendulum WEC with or without tide adaptation mechanism are built to compare their performance on various water levels. When comparing these two WECs in terms of their average output based on the linear passive control strategy, the output power of WEC with the tide adaptation mechanism is much steadier with the change of the water level and always larger than that without the tide adaptation mechanism.     

Effect of Sea Level Variation on Calculation of Design Water Level

The long-term variation and seasonal variation of sea level have a notable effect on the calculation of engineering water level. Such an effect is first analyzed in this paper. The maximal amplitude of inter-annual anomaly of monthly mean sea level along the China coast is larger than 60 cm. Both the storm surge disaster and cold wave disaster are seasonal disasters in various regions, so the water level corresponding to the 1% of the cumulative frequency in the cumulative frequency curve of hourly water level data for different seasons in various sea areas is different from design water level, for example, the difference between them reaches maximum in June, July and August for northern sea area, and maximum in September, October and November for Southern China Sea. The hourly water level data of 19 gauge stations along the China coast are analyzed. Firstly, the annual mean sea level for every station is obtained; secondly, linear chan ging rates of annual mean sea level are obtained with the stochasti     

胶州湾的潮汐与潮流

胶州湾是我国沿海的一个重要港湾,长期以来对其进行了较多的调查研究工作。大港验潮站自1926年起开始观测(1936-1946年曾中断),目前已积累了44年的长期水位资料。山角底验潮站自1967年开始观测,也积累了16年的水位资料。我们还系统地搜集了我所1958-1959年在薛家岛、东洋嘴、团岛、麦岛、大公岛分别进行的四个月、两个月及半个月的观测资料,以及山东海洋学院1975年在红岛船厂(阴岛)、黄岛客运码头分别进行的三个月及两个月的水位观测资料和团岛湾、后岔湾、大石头等地的潮汐调和常数。此外,我们还搜集了胶州湾及其附近海区共二十三个测站的海流连续观测资料,其中十九个站的资料是由我所在1957-1959年观测的,三个站是1982年“全国海岸带调查”中测得的,另一个站是华东水利学院提供的资料。这些测站,除六个站为一昼夜连续观测外,其余十七个站均为二昼夜以上连续观测(有一个站为十五昼夜连续观测)。观测站位见图1。 我们把搜集到的这些潮汐、潮流资料进行了调和分析,得到了各主要分潮的调和常数及潮汐、潮流特征值。本文以此为基础,结合实测资料,分析了胶州湾潮汐、潮流的基本特征。     

深圳海域潮汐海啸波耦合数值研究

以COMCOT海啸模式和TPXO7.1全球潮汐模式为基础,采用三层嵌套网格,建立了南海海啸与潮汐耦合计算模型,分析深圳海域海啸和潮汐相互作用。潮汐计算结果与实测数据吻合较好,高、低潮位平均误差小于15 cm,20 cm;在潮汐验证的基础上,以马尼拉海沟潜在地震海啸源为案例,进行8.0,9.0级地震海啸与潮汐耦合情景模拟计算,计算结果表明,9级地震海啸在深圳海域外海波高为140~150 cm,如先行波为正波发生在高潮时将产生异常高潮位,负波发生在低潮时将产生异常低潮位,线性叠加计算结果偏大,在25.0 cm之内,到达时间差异小于6 min。