印度尼西亚近海潮汐潮流的数值模拟

利用FVCOM海洋数值模式计算了印尼近海的M2,S2,K1,O1分潮的分布,计算范围从20°S～20°N,90°～150°E,计算网格分辨率在印尼海域岛屿平均为1/12度,在大陆边界平均为1/5度,在开边界平均为1/2度.计算结果与104个TOPEX/Poseidon卫星高度计交叉点数据和79个验潮站数据进行比较,符合良好;与高度计交叉点比较,M2分潮振幅的均方根差为6 cm,迟角为7°;S2分潮的振幅偏差为3 cm,迟角偏差为8°;K1分潮振幅的偏差为6 cm,迟角偏差为10°;O1分潮振幅偏差为3 cm,迟角偏差为10°.根据计算结果给出了4个分潮的潮汐、潮流、潮余流和潮能通量密度分布图.     

渤海潮波系统的数值模拟

利用二维非线性潮波方程组,讨论了渤黄海主要分潮(全日潮、半日潮及浅水分潮) 数值模拟中的有关问题。数值模拟中同时考虑了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°。实验结果较好地体现了渤黄海潮波系统的特征。     

基于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°,结果显示各分潮模拟结果合理,较真实地反映了渤海海域四个分潮传播情况。     

两种不同水动力模型对渤海潮汐模拟的比较研究

本文在考虑相同物理过程以及网格分辨率基础上,分别利用FVCOM (Finite-volume coastal ocean model)模型和DELFT 3D模型对渤海潮汐运动进行数值模拟,并将两种模型模拟结果与实测数据进行对比。模拟结果显示,两种模型都能较好模拟出渤海潮波系统分布情况。FVCOM模拟的M2分潮振幅和迟角均方根误差分别为8.63cm、7.41°;DELFT 3D模拟的M2分潮振幅和迟角均方根误差分别为8.97cm、8.39°。最后采用FVCOM模型对辽东湾海域进行局部加密,分析表明非结构网格局部加密方案可以进一步提高加密区域的潮汐模拟精度。     

利用T/P 卫星高度计资料调和分析南海潮汐信息

利用j,v模型调和分析1992～2002年共10 a的TOPEX/Poseidon(T/P)海面高度距平资料,提取了南海K1,O1,P1,Q1,M2,S2,N2和K2等8个主要分潮的潮汐调和常数。分析比较了卫星上下行轨道的19个交叉点的振幅和迟角,其中M2,S2,K1和O1的平均向量均方根偏差分别是1.5,1.1,2.5和1.4 cm;将交叉点的调和常数与TPXO7.2模式的结果进行了比较,结果表明M2,S2,K1和O1分潮振幅的绝对平均误差均小于3 cm,迟角的最大绝对平均误差为7.8°。选取了与卫星轨道较近的8个验潮站,对验潮站的实测数据调和常数和本文所得调和常数进行了比较,结果显示K1分潮的向量均方根偏差为4.7 cm,M2分潮的向量均方根偏差为3.7 cm。论文结果表明利用j,v模型调和分析方法对南海海域卫星高度计资料进行潮汐信息提取是可靠的,并可为局部重力场的研究提供海洋潮汐改正数据,有一定的参考价值。     

温州鳌江近海建设工程环境影响潮汐潮流数值模拟(英文)

收集了近年来鳌江口附近海域多个工程不同阶段5个潮位站的3~5年潮位实测数据和部分海流实测资料,通过对鳌江口附近海域的不同年份的水位资料进行潮汐调和常数分析,鳌江近海海域主要为半日潮区,其中M2分潮的振幅在170 cm~193 cm;迟角在260°~280°之间,这些站的2007年、2010年、2011年调和常数分析结果相比,主要的半日分潮M2、S2、N2,全日分潮K1、O1及浅水分潮M4、MS4、及M6等分潮振幅、迟角的最大变化分别在1.8 cm~4.4 cm和3°~7°之间。在初步掌握了鳌江口潮汐潮流特征的基础上,采用无结构的三角形网格和有限体积法的FVCOM海洋数值模型,进行模拟结果验证,计算结果与实测数据符合良好。构建重点年份建设工程合拢产生新的岸线水深的潮汐潮流场,刻画鳌江口建设工程的叠加影响。     

温州鳌江近海建设工程环境影响潮汐潮流数值模拟

收集了近年来鳌江口附近海域多个工程不同阶段5个潮位站的3-5年潮位实测数据和部分海流实测资料,通过对鳌江口附近海域的不同年份的水位资料进行潮汐调和常数分析,鳌江近海海域主要为半日潮区,其中M2分潮的振幅在170~193 cm;迟角在260°~280°之间,这些站的2007年、2010年、2011年调和常数分析结果相比,主要的半日分潮M2、S2、N2,全日分潮K1、O1及浅水分潮M4、MS4、及M6等分潮振幅、迟角的最大变化分别在1.8~4.4 cm和3°~7°之间。在初步掌握了鳌江口潮汐潮流特征的基础上,采用无结构的三角形网格和有限体积法的FVCOM海洋数值模型,进行模拟结果验证,计算结果与实测数据符合良好。构建重点年份建设工程合拢产生新的岸线水深的潮汐潮流场,刻画鳌江口建设工程的叠加影响。     

基于FVCOM的厦门湾及其周边海域三维潮流数值模拟

基于非结构三角形网格、干-湿判别技术和有限体积法的FVCOM（finite volume coastal oceanmodel）海洋数值模式,建立了厦门湾及其周边海域高分辨率（30 m）的三维潮汐、潮流数值模型.模拟结果同该海域2个验潮站和4个连续海流站的观测资料符合良好,较好地反映了厦门湾及其周边海域潮汐、潮流运动的变化状况和分布特征,并给出了M2、S2、K1、O1共4个主要分潮的同潮图、表层潮流椭圆、最大可能潮流流速及表、底层潮余流分布.厦门湾及其周边海域属正规半日潮类型,4个分潮的最大潮汐振幅分别为200、65、36、29 cm,厦门湾内外迟角差分别为20°、25°、18°、10°;镇海角至围头角连线东南侧湾口区为逆时针旋转的驻波,西北侧湾内为前进潮波.湾内潮流属正规半日潮流,湾口区潮流以逆时针方向的旋转流运动为主,湾内各水道为往复潮流,椭圆长轴与水道走向一致,4个分潮流表层最大流速分别为201、51、34、25 cm/s;九龙江口区3条港道内的流速以南港最大;表层余流大于底层余流,二者水平分布形态基本一致,都为北进南出.     

基于FVCOM的温州近海潮汐潮流数值模拟

基于有限体积法海洋数值模型(FVCOM),构建了温州近海潮汐潮流数值模式,模式模拟区域为(120°24′00″~121°19′12″E,27°21′00″~28°24′00″N),模式水平分辨率由近岸河口区的50m,逐渐增加至开边界附近的2km。模式模拟并分析了温州近海的M2,S2,N2,K1,O1五个主要分潮。利用温州近海实测资料对模拟结果进行了验证,模拟与实测符合良好;其中与4个验潮站资料比较,M2,S2,N2,K1,O1五个主要分潮的振幅绝均差和迟角绝均差分别为4.84cm和5.14°,2.19cm和3.35°,5.18cm和4.38°,0.64cm和3.67°,0.59cm和4.61°;与9个海流连续观测站比较,流速绝均差为11.71cm/s,流向绝均差为9.66°。在模拟结果较好地反映温州近海潮汐、潮流运动状况的基础上,本文给出了各模拟分潮的潮汐同潮图和潮流椭圆分布、潮汐和潮流类型分布以及最大可能潮流分布等。     

中国南海北部潮汐主要分潮的变化趋势分析

基于中国南海北部6个长期验潮站逐时水位资料,对该海域O1、K1、M2和S2的4个主要分潮的振幅和迟角变化趋势分别进行分析。结果表明:各验潮站的O1、K1和M2分潮振幅呈现显著的周期性变化,S2分潮振幅年变化量基本在毫米级,具有较高的稳定性。厦门站1954~1997年期间的M2分潮振幅存在18.98年的周期变化,与交点潮18.61年的变化周期相吻合。基于最小二乘原理,通过曲线拟合方式能够对主要分潮振幅在一定时间内进行比较准确的预报。     

Cotidal charts and tidal power input atlases of the global ocean from TOPEX/Poseidon and JASON-1 altimetry

The global distributions of eight principal tidal constituents, M₂ , S₂ , K₁ , O₁ , N₂ , K₂ , P₁ , and Q₁ , are derived using TOPEX/Poseidon and JASON-1（T/P-J） satellite altimeter data for 16 a. The intercomparison of the derived harmonics at 7000 subsatellite track crossover points shows that the root mean square （RMS） values of the tidal height differences of the above eight constituents range from 1.19 cm to 2.67 cm, with an average of about 2 cm. The RMS values of the tidal height differences between T/P-J solutions and the harmonics from ground measurements at 152 tidal gauge stations for the above constituents range from 0.34 cm to 1.08 cm, and the relative deviations range from 0.031 to 0.211. The root sum square of the RMS differences of these eight constituents is 2.12 cm, showing the improvement of the present model over the existing global ocean tidal models. Based on the obtained tidal model the global ocean tidal energetics is studied and the global distribution of the tidal power input density by tide-generating force of each constituent is calculated, showing that the power input source regions of semidiurnal tides are mainly concentrated in the tropical belt between 30 S and 30 N, while the power input source regions of diurnal tides are mainly concentrated off the tropic oceans. The global energy dissipation rates of the M₂ , S₂ , K₁ , O₁ , N₂ , P₁ , K₂ and Q₁ tides are 2.424, 0.401, 0.334, 0.160, 0.113, 0.035, 0.030 and 0.006 TW, respectively. The total global tidal dissipation rate of these eight constituents amounts to 3.5 TW.     

Estimation of baroclinic tide energy available for deep ocean mixing based on three-dimensional global numerical simulations

The global distributions of the major semidiurnal (M₂ and S₂) and diurnal (K₁ and O₁) baroclinic tide energy are investigated using a hydrostatic sigma-coordinate numerical model. A series of numerical simulations using various horizontal grid spacings of 1/15–1/5° shows that generation of energetic baroclinic tides is restricted over representative prominent topographic features. For example, nearly half of the diurnal (K₁ and O₁) baroclinic tide energy is excited along the western boundary of the North Pacific from the Aleutian Islands down to the Indonesian Archipelago. It is also found that the rate of energy conversion from the barotropic to baroclinic tides is very sensitive to the horizontal grid spacing as well as the resolution of the model bottom topography; the conversion rate integrated over the global ocean increases exponentially as the model grid spacing is reduced. Extrapolating the calculated results in the limit of zero grid spacing yields the estimate of the global conversion rate to be 1105 GW (821, 145, 102, 53 GW for M₂, S₂, K₁, and O₁ tidal constituents, respectively). The amount of baroclinic tide energy dissipated in the open ocean below a depth of 1000 m, in particular, is estimated to be 500–600 GW, which is comparable to the mixing energy estimated by Webb and Suginohara (Nature 409:37, 2001) as needed to sustain the global overturning circulation.     

Phase distribution of tidal constituents around New Zealand

The phase distributions of the M₂, S₂, K₁, and O₁ tidal constituents around New Zealand are plotted from existing harmonic analyses of tidal heights. Both semidiurnal constituents exhibit a complete 360° range of phase around New Zealand, with complex areas of rapid phase change through or near the strait separating the two main islands. The K₁ amphidrome and that for O₁, which previously were thought to be centred on New Zealand, are shown to be located cast of New Zealand. The distributions plotted highlight areas where tidal observations are lacking.     

珠江河口潮能收支研究

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.     

Numerical study of baroclinic tides in Luzon Strait

The spatial and temporal variations of baroclinic tides in the Luzon Strait (LS) are investigated using a three-dimensional tide model driven by four principal constituents, O₁, K₁, M₂ and S₂, individually or together with seasonal mean summer or winter stratifications as the initial field. Barotropic tides propagate predominantly westward from the Pacific Ocean, impinge on two prominent north-south running submarine ridges in LS, and generate strong baroclinic tides propagating into both the South China Sea (SCS) and the Pacific Ocean. Strong baroclinic tides, ∼19 GW for diurnal tides and ∼11 GW for semidiurnal tides, are excited on both the east ridge (70%) and the west ridge (30%). The barotropic to baroclinic energy conversion rate reaches 30% for diurnal tides and ∼20% for semidiurnal tides. Diurnal (O₁ and K₁) and semidiurnal (M₂) baroclinic tides have a comparable depth-integrated energy flux 10–20 kW m⁻¹ emanating from the LS into the SCS and the Pacific basin. The spring-neap averaged, meridionally integrated baroclinic tidal energy flux is ∼7 GW into the SCS and ∼6 GW into the Pacific Ocean, representing one of the strongest baroclinic tidal energy flux regimes in the World Ocean. About 18 GW of baroclinic tidal energy, ∼50% of that generated in the LS, is lost locally, which is more than five times that estimated in the vicinity of the Hawaiian ridge. The strong westward-propagating semidiurnal baroclinic tidal energy flux is likely the energy source for the large-amplitude nonlinear internal waves found in the SCS. The baroclinic tidal energy generation, energy fluxes, and energy dissipation rates in the spring tide are about five times those in the neap tide; while there is no significant seasonal variation of energetics, but the propagation speed of baroclinic tide is about 10% faster in summer than in winter. Within the LS, the average turbulence kinetic energy dissipation rate is O(10⁻⁷) W kg^{− 1} and the turbulence diffusivity is O(10⁻³) m²s⁻¹, a factor of 100 greater than those in the typical open ocean. This strong turbulence mixing induced by the baroclinic tidal energy dissipation exists in the main path of the Kuroshio and is important in mixing the Pacific Ocean, Kuroshio, and the SCS waters.     

Characteristics of Tides in the Bay of Bengal

A vertically integrated 2D numerical model was developed for the simulation of major tidal constituents (M₂, S₂, N₂, K₁ and O₁) in the Bay of Bengal. The bathymetry for the model domain was derived from an improved ETOPO5 dataset prepared in our earlier work. The simulated tidal elevations showed good agreement with the hourly tide gauge observations at Paradip, Visakhapatnam, and Chennai. The amplitudes and phases of M₂, S₂, K₁, and O₁ at the coastal stations, obtained from harmonic analysis of simulated tides, were found to agree well with those obtained from Admiralty Tide Tables with the RMS misfit 9.2, 5.6, 2.9 and 3.1 cm, respectively. In the Bay of Bengal, semi-diurnal tides (M₂, S₂, and N₂) attain highest amplitudes (180, 80, 30 cm, respectively) in the Gulf of Martaban while amplitudes of diurnal tides (K₁, O₁) reaches maximum (20, 12 cm, respectively) in the Malacca Strait. The continental shelf in the head bay and along the southern coast of Myanmar is about 200 km wide and the amplitudes of semi-diurnal tides are doubled in these regions while the diurnal tides amplify only marginally, which is consistent with Clarke and Battisti theory. In the north eastern end of the head bay and the Gulf of Martaban, the geometrical configuration of the coastline, in addition to the wide continental shelf, could contribute to the amplification of both semi-diurnal and diurnal constituents. In the Malacca Strait, the amplitudes of both semi-diurnal and diurnal tides are found to increase gradually from the northern end to the 2.5°N and decreases towards southern boundary. The co-tidal and co-range charts of M₂ and S₂ tidal constituents also show the presence of two degenerate amphidromic points in the head bay. A virtual amphidromic point for M₂ is identified in the Malacca Strait.     

烟墩角海域二测站周日水温及潮流变化特征

为加深对山东半岛东部养殖区的海流、海温变化特征,及其与该地区生态特征分布之间关系的认识,应用短期资料的潮流准调和分析方法,计算了烟墩角海域两个测站的O1,K1,M2,S2,M4,MS46个主要分潮的北、东分量潮流调和常数,并给出了各测站在各层的潮流椭圆要素。结果表明:该海区水温呈现层化结构,在1天内呈现两个波动,与N-S向潮流变化关系密切。海区除湾内表层外属日潮性质外,潮流属于正规半日潮流性质,半日分潮流的北分量大于东分量,呈现典型往复流特征,最大流流向多为南北向。分析结果也反映出表层、中层的浅水分潮在观测海流中所占的份额高于底层的浅水分潮。     

泰国湾及邻近海域潮汐潮流的数值模拟

本文基于FVCOM(Finite-Volume Coastal Ocean Model)模式,模拟了泰国湾及其周边海域K1、O1、M2和S2四个主要分潮。采用47个验潮站实测调和常数与模拟结果进行比较,所得4个分潮的均方差分别为4.06cm、3.76cm、8.22cm和4.71cm,符合良好。根据计算结果分析了泰国湾及其周边海域的潮汐、潮流的分布特征和潮波的传播特征。数值试验表明,现有的数字水深资料(ETOPO1,ETOPO5,DBDB-V)的准确度不足以合理地模拟泰国湾潮波。