内潮耗散与自吸-负荷潮对南海潮波影响的数值研究

利用非结构三角形网格的FVCOM海洋数值模式,在其传统二维潮波方程中加入参数化的内潮耗散项和自吸-负荷潮项,计算了南海及其周边海域的M_2、S_2、K_1和O_1分潮的分布。与实测值的比较表明,引入这两项对模拟准确度的提高有明显效果。根据模式结果本文计算分析了研究海域的潮能输入和耗散。能量输入计算表明,能通量是潮能输入的最主要构成部分,通过吕宋海峡断面进入南海的M_2和K_1分潮能通量分别为38和29GW;半日周期的自吸-负荷潮能量输入以负值居多,而全日周期的自吸-负荷潮能量输入以正值居多,因而自吸-负荷潮减弱了南海的半日潮,并加强了南海的全日潮。引潮力的作用也减弱了半日潮而加强了全日潮,但其作用要小于自吸-负荷潮。潮能耗散的分析显示底摩擦耗散在沿岸浅水区域起主导作用,内潮耗散则主要发生在深水区域。内潮耗散的最大值出现在吕宋海峡,且位于南海之外的海峡东部的耗散量大于位于南海之内的海峡西部的耗散量。对M_2和K_1分潮吕宋海峡的内潮耗散总值分别达到16和23GW。     

Numerical study of tidal dynamics in the South China Sea with adjoint method

We adopt a parameterized internal tide dissipation term to the two-dimensional (2-D) shallow water equations, and develop the corresponding adjoint model to investigate tidal dynamics in the South China Sea (SCS). The harmonic constants derived from 63 tidal gauge stations and 24 TOPEX/Poseidon (T/P) satellite altimeter crossover points are assimilated into the adjoint model to minimize the deviations of the simulated results and observations by optimizing the bottom friction coefficient and the internal tide dissipation coefficient. Tidal constituents M₂, S₂, K₁ and O₁ are simulated simultaneously. The numerical results (assimilating only tidal gauge data) agree well with T/P data showing that the model results are reliable. The co-tidal charts of M₂, S₂, K₁ and O₁ are obtained, which reflect the characteristics of tides in the SCS. The tidal energy flux is analyzed based on numerical results. The strongest tidal energy flux appears in the Luzon Strait (LS) for both semi-diurnal and diurnal tidal constituents. The analysis of tidal energy dissipation indicates that the bottom friction dissipation occurs mainly in shallow water area, meanwhile the internal tide dissipation is mainly concentrated in the LS and the deep basin of the SCS. The tidal energetics in the LS is examined showing that the tidal energy input closely balances the tidal energy dissipation.     

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

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

Tidal energy balance in the Seto Inland Sea

The tidal volume transport in the Seto Inland Sea is calculated. The cross-section where the volume transport of the M₂ tide is zero, is located around the western part of Bisan Strait. The tidal energy dissipation of the M₂ tide by friction is 6.30×10¹⁶ ergs s^–1 in the Seto Inland Sea. The quality factorQ for the M₂ tide is 20.2. The total energy dissipation of the M₂, S₂, K₁ and O₁ tides is 7.99×10¹⁶ ergs s^–1.     

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

基于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%左右。根据计算结果分析了印尼海域的潮汐特征及潮能传播规律,结果显示:爪哇海以外的印尼海域主要为不规则半日潮区;全日潮潮能主要由太平洋传入印尼海域,而半日潮潮能则是从印度洋传入印尼海域。     

北部湾潮波数值研究

利用普林斯顿海洋模式(POM08)建立了北部湾及其临近海区潮汐潮流数值模式,模拟了K1,O1,M2和S2这4个主要分潮,分析了模拟的潮汐和潮流分布特征,从潮波能量的角度讨论了琼州海峡对北部湾潮波系统的影响,并给出北部湾潮能的耗散情况。研究表明,北部湾是典型的全日潮海区,K1和O1分潮在南部湾口形成半个旋转潮波系统,无潮点位于越南顺安附近岸边。琼州海峡中的欧拉潮汐余流为西向流,潮余流造成的水通量约为0.034×106m3/s;余流出海峡西口后,先折向北,然后转向南流出湾外。研究海区中两个强潮流区分别位于琼州海峡和海南岛的西侧,同时这也是两个潮能的高耗散区。北部湾的潮能自南部湾口由外海传入,通过西口涌入琼州海峡,到达海峡东口时日潮波的能量已基本耗散殆尽,在海峡内耗散的4个分潮的潮能约为3.33 GW,相当于北部湾潮能耗散量的35%左右。数值试验表明,琼州海峡作为潮能耗散的重要海区,其存在对于北部湾潮波系统的形成具有较大影响。计算了底边界潮能耗散,结果表明在北部湾和琼州海峡,底边界耗散的潮能分别占该海区总耗散的83%和80%。     

珠江河口潮能收支研究

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.     

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.     

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.     

渤海潮波系统的数值模拟

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

Observed and modelled tidal currents in the New Zealand region

Tidal currents derived from current meter measurements are compared with the output from a barotropic tidal model of the New Zealand region. For the semi‐diurnal constituents there was very good agreement for the M₂ tide and good agreement for the S₂ tide. For the diurnal constituents (K_l, O_l) it was found that as the amplitude of the constituents decreased so did both the model/observation agreement and the accuracy of the observed tidal ellipse parameters. Consequently it was not possible to decide whether differences arose through shortcomings in the model or in the data. However, the overall performance of the model as a prognostic tool for ocean tidal current simulation appears to be good.     

The tidal regime of Shark Bay, Western Australia

A non-linear hydrodynamic model is used to describe the tidal dynamics of Shark Bay, Western Australia. The model is forced by tidal elevations generated by M₂, S₂, K₁ and O₁ constituent data at the open boundaries. The absence of suitable boundary data required a ‘calibration’ of the boundary condition against the known constituent data from within the model domain. The model provides a good match to the available field data, and allows the surface-level and current response to be resolved over the entire domain. Due to a near quarter-wave resonance of the semi-diurnal tide along the eastern Hopeless Reach, which increases the semi-diurnal tide by a factor of 2, the tidal characteristics on each of the Reaches are different: on the eastern Hopeless Reach the tides are mainly semi-diurnal while on the western Freycinet Reach the tides are mainly diurnal. The tidal range is also higher along Hopeless Reach. Tidal harmonics, generated by non-linearity, are important in the shallow regions. The tidal wave is shown to propagate as a progressive wave into the Bay. Substantial phase-lag, attenuation and dissipation occur over the Faure Sill, a major shallow region of the eastern reach of the Bay. Non-linear generation of the M₄ and MS₄ tides is also significant in this region. Depth-averaged residual currents are presented, which show a tidally generated circulation that is enhanced in regions of complex topography. Estimates of tidal dissipation indicate that although the total dissipation is small on a global scale, the areal average is comparable with the Gulf of Carpentaria and approximately one-quarter of the value estimated for the Patagonian Shelf.     

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.     

Three-dimensional modeling of tidal circulation and mixing over the Java Sea

A combination of a three-dimensional hydrodynamic model and in-situ measurements provides the structures of barotropic tides, tidal circulation and their relationship with turbulent mixing in the Java Sea, which allow us to understand the impact of the tides on material distribution. The model retains high horizontal and vertical resolutions and is forced by the boundary conditions taken from a global model. The measurements are composed of the sea level at coastal stations and currents at moorings embedded in Seawatch buoys, in addition to hydrographic data. The simulated tidal elevations are in good agreement with the data for the K₁ and M₂ constituents. The K₁ tide clearly shows the lowest mode resonance in the Java Sea with intensification around the nodal point in the central region. The M₂ tide is secondary and propagates westward from the eastern open boundary, along with a counterclockwise amphidromic point in the western part. The K₁ tide produces a major component of tidal energy, which flows westward and dissipates through the node region near the Karimata Strait. Meanwhile, the M₂ tide dissipates in the entire Java Sea. However, the residual currents are mainly induced by the M₂ tide, which flows westward following the M₂ tidal wave propagation. The tidal mixing is mainly caused by K₁ tide which peaks at the central region and is consistent with the uniform temperature and salinity along the vertical dimension. This mixing is expected to play an important role in the vertical exchange of nutrients and control of biological productivity.     

A study of non-linear tidal propagation in shallow inlet/estuarine systems Part I: Observations

The offshore tide becomes strongly distorted as it propagates into shallow estuarine systems. Observations of sea surface elevation and horizontal currents over periods ranging from three days to one year, at nine stations within Nauset inlet/estuary, document the non-linear interaction of the off-shore equilibrium tidal constituents. Despite strong frictional attenuation within the estuary, the overtides and compound tides of M₂, S₂ and N₂, in particular, reach significant amplitude, resulting in strong tidal distortion. High frequency forced constituents in sea surface are phase-locked, consistently leading the forcing tides by 60–70°, resulting in a persistent distortion where falling tide is longer than rising tide. Forced constituents in currents are more nearly in phase with equilibrium constituents, producing flood currents which are shorter but more intense than ebb currents. A compound fortnightly tide, MS_f, modulates the mean water level such that lowest tides occur during neap phase instead of spring phase. This fortnightly tide can be contaminated by storm surge, changing the phase characteristics of this constituent. Implications of the overtides, compound tides, and lower frequency tides on near-bed, suspended and dissolved material transport are profound.     

Tidal Currents in the Tsushima Straits Estimated from ADCP Data by Ferryboat

Tidal currents in the Tsushima Straits have been analyzed using measurements obtained since February 1997 by an acoustic Doppler current profiler (ADCP) mounted on the ferryboat Camellia. Tidal current constituents (M ₂, S ₂, K ₁, O ₁) are dominant among the ten tidal current constituents (Q ₁, O ₁, P ₁, K ₁, N ₂, M ₂, S ₂, K ₂, MSf, Mf), and generally 1.4–2.1 times stronger at the western channel of the straits than those at the eastern channel. The ratio between amplitude of M ₂, S ₂, K ₁ and O ₁ averaged along the ferryboat track is 1:0.45:0.59:0.51. The major axis directions of tidal current ellipses are generally SW to NE, exceptionally in the vicinity of the Tsushima Islands. Approaching the Tsushima Islands from the Korean Peninsula side, the major axis gradually rotates clockwise. At the western channel, the M ₂ and K ₁ constituents change the rotation direction of current vectors from clockwise to counterclockwise at about 90–130 m depth. The contributions of the tidal currents to the mean kinetic energy and the mean eddy kinetic energy along the ferryboat track are, on average, 0.56 and 0.71, respectively. This suggests that tidal current activities are generally more dominant than the mean current activities and much more dominant than eddy activities. The only region where the eddy activities are comparable to the tidal current activities is located on the east side of the Tsushima Islands. This revised version was published online in July 2006 with corrections to the Cover Date.     

南麂岛附近海域潮汐和潮流的特征

以2008年冬季在浙江近海南麂岛附近投放的4个底锚系观测的水位和流速资料为依据,分析了潮汐和潮流特征。水位谱分析结果显示半日分潮最显著,全日分潮其次;近岸的浅水分潮比离岸大。水位调和分析结果表明:潮汐类型均为正规半日潮,近岸处的平均潮差大于3m,最大可能潮差大于6m,潮汐呈现出显著的低潮日不等和回归潮特征。流速谱分析结果显示半日分潮流最强,全日分潮流其次,且比半日分潮流小得多;近岸浅水分潮流比远离岸显著。流速调和分析结果表明:潮流类型均为正规半日潮流,靠近岸的两个站浅水分潮流较显著;最显著的半日分潮流是M2分潮流,其最大流速介于0.32～0.48m/s之间,全日分潮流均很弱,最大流速小于0.06m/s。M2分潮流均为逆时针旋转,椭圆率越靠近海底越大;最大分潮流流速分布为中上层最大、表层略小、底层最小;最大分潮流流速方向的垂向变化很小,底层比表层略为偏左;最大分潮流流速到达时间随深度的加深而提前,底层比中上层约提前30min。潮流椭圆的垂向分布显示这里的半日分潮流以正压潮流为主;日分潮流则表现出很强的斜压性。     

胶州湾潮汐潮流动边界数值模拟

基于普林斯顿海洋模式,通过干湿网格判别法引入潮汐潮流的漫滩过程,考虑M₂,S₂,K₁,O₁,M₄和MS₄六个主要分潮,建立了胶州湾潮汐潮流数值模拟和预报模型,研究了该海域潮汐潮流特征,并讨论了漫滩对潮流模拟的影响。与实测资料的对比验证表明,该模式能够对胶州湾的潮汐和潮流做出较为合理的预测。给出了胶州湾潮汐、潮流、余流等分布特征,模拟的潮流场以及余流场涡旋等现象与观测符合良好;计算了潮波能通量,从能量角度探讨了潮波的传播特性;对潮位与潮流场演变规律,以及潮能通量的分析表明,胶州湾内的潮波以驻波为主。通过数值试验发现,漫滩过程的引入对胶州湾潮流速度的模拟至关重要,不考虑漫滩过程的模式会夸大或者低估潮流流速。对于滩涂面积广阔的海域来说,潮流数值模式中考虑漫滩的影响是必要的。     

Ocean Tide Models Developed by Assimilating TOPEX/POSEIDON Altimeter Data into Hydrodynamical Model: A Global Model and a Regional Model around Japan

A global ocean tide model (NAO.99b model) representing major 16 constituents with a spatial resolution of 0.5° has been estimated by assimilating about 5 years of TOPEX/POSEIDON altimeter data into barotropic hydrodynamical model. The new solution is characterized by reduced errors in shallow waters compared to the other two models recently developed; CSR4.0 model (improved version of Eanes and Bettadpur, 1994) and GOT99.2b model (Ray, 1999), which are demonstrated in comparison with tide gauge data and collinear residual reduction test. This property mainly benefits from fine-scale along-track tidal analysis of TOPEX/POSEIDON data. A high-resolution (1/12°) regional ocean tide model around Japan (NAO.99Jb model) by assimilating both TOPEX/POSEIDON data and 219 coastal tide gauge data is also developed. A comparison with 80 independent coastal tide gauge data shows the better performance of NAO.99Jb model in the coastal region compared with the other global models. Tidal dissipation around Japan has been investigated for M₂ and K₁ constituents by using NAO.99Jb model. The result suggests that the tidal energy is mainly dissipated by bottom friction in localized area in shallow seas; the M₂ ocean tidal energy is mainly dissipated in the Yellow Sea and the East China Sea at the mean rate of 155 GW, while the K₁ energy is mainly dissipated in the Sea of Okhotsk at the mean rate of 89 GW. TOPEX/POSEIDON data, however, detects broadly distributed surface manifestation of M₂ internal tide, which observationally suggests that the tidal energy is also dissipated by the energy conversion into baroclinic tide.     

江苏海岸中部近岸冬季潮汐和潮流特征

南黄海西侧的江苏海岸近岸区域,素以地形复杂、潮流强劲、悬沙输运剧烈著称,但是较长期的同步潮位和潮流观测数据仍然缺乏,尤其是在近岸(20 km)浅水(20 m)区域。2014年1月在大丰港附近开展了连续潮位和潮流观测,获得的数据揭示了一系列特征。此地潮汐潮流为正规半日潮,浅水分潮显著。平均潮差为3.05 m,最显著的两个分潮为M2和S2分潮,振幅分别为1.45 m和0.52 m。潮流最显著的半日分潮M2分潮和最显著的浅水分潮M4分潮在沿岸方向上振幅分别为0.84m/s和0.12m/s,在跨岸方向上振幅分别为0.24 m/s和0.01 m/s,沿岸方向占绝对优势。潮波的沿岸传播介于前进波和驻波之间,驻波的特征稍强。M2分潮潮流椭圆最大流(长轴)方向为南偏东7.4°。存在冬季沿岸向北的余流,垂向平均值的大小为2.2 cm/s。以上潮汐潮流特征为该区域海洋物质输运研究提供了基础资料。