南海潮汐的伴随同化数值模拟

把利用正交潮响应方法对 2 4 8个周期超过 6年的南中国海的TOPEX/Poseidon卫星高度计资料进行潮波分析提取的沿轨分潮调和常数同化到二维非线性潮汐数值模式中去 ,优化模型中的开边界条件和底摩擦系数 ,模拟了南海m1 和M2 分潮的潮汐。所用的同化方法是伴随同化。根据计算结果给出了m1 和M2 分潮的同潮图。计算结果与 5 9个验潮站资料的比较结果是 :m1 分潮的振幅和迟角的平均绝对误差分别是 4.8cm和 8.7°;M2 分潮的振幅和迟角的平均绝对误差分别是 4.3cm和 1 1 .0°,表明计算结果与验潮站资料符合良好。研究结果表明 ,利用伴随同化方法把TOPEX/Poseidon资料同化到潮汐数值模式中去对模式进行校正是有效的     

渤海开边界潮汐的伴随法反演

潮汐潮流数值模拟中的一个主要难点在于开边界条件的确定。本文采用伴随法 ,由渤海沿岸 1 9个验潮站的潮汐调和常数来反演渤海海域的开边界条件 ,以实现渤海潮波的数值模拟。计算所得调和常数与实测值之差的绝对平均值 :m1 潮波振幅差为 1 4cm ,迟角差为5 0°;M2 潮波振幅差为 2 4cm ,迟角差为 5 0°。数值模拟结果正确地反映了渤海m1 和M2 潮波的基本特征     

渤海潮波系统的数值模拟

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

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.     

Oceanic tidal angular momentum and Earth's rotation variations

Luni-solar tides affect Earth's rotation in a variety of ways. We give an overview of the physics and focus on the excitation of Earth rotational variations by ocean tides under the conservation of angular momentum. Various models for diurnal and semidiurnal tidal height and tidal current fields have been derived, following a legacy of a number of theoretical tide models, from the Topex/Poseidon (T/P) ocean altimetry data. We review the oceanic tidal angular momenta (OTAM) predicted by these T/P models for the eight major tides (Q₁, O₁, P₁, K₁, N₂, M₂, S₂, K₂), and their excitations on both Earth's rotational speed variation (in terms of length-of-day or UT1) and polar motion (prograde diurnal/semidiurnal components and retrograde semidiurnal components). These small, high-frequency effects have been unambiguously observed in recent years by precise Earth rotation measurements via space geodetic techniques. Here we review the comparison of the very-long-baseline-interferometry (VLBI) data with the T/P OTAM predictions. The agreement is good with discrepancies typically within 1 – 2 microseconds for UT1 and 10 – 30 microarcseconds for polar motion. The eight tides collectively explain the majority of subdaily Earth rotation variance during the intensive VLBI campaign Cont94. This establishes the dominant role of OTAM in exciting the diurnal/semidiurnal polar motion and paves the way for detailed studies of short-period non-OTAM excitations, such as atmospheric and oceanic angular momentum variations, earthquakes, the atmospheric thermal tides, Earth librations, and the response of the mantle lateral inhomogeneities to tidal forcing. These studies await further improvements in tide models and Earth rotation measurements.     

Tidal oscillations in the Caspian Sea

Long-term hourly data from 12 tide gauge stations were used to examine the character of tidal oscillations in the Caspian Sea. Diurnal and semidiurnal tidal peaks are well-defined in sea level spectra in the Middle and South Caspian basins. High-resolution spectral analysis revealed that the diurnal sea level oscillations in the Middle Caspian Basin have a gravitational origin, while those in the South Caspian Basin are mainly caused by radiational effects: the amplitude of diurnal radiational harmonic S₁ is much higher than those of gravitational harmonics О₁, P₁, and K₁. In the North Caspian Basin, there are no gravitational tides and only weak radiational tides are observed. A semidiurnal type of tide is predominant in the Middle and South Caspian basins. Harmonic analysis of the tides for individual annual series with subsequent vector averaging over the entire observational period was applied to estimate the mean amplitudes and phases of major tidal constituents. The amplitude of the M₂ harmonic reaches 5.4 cm in the South Caspian Basin (at Aladga). A maximum tidal range of 21 cm was found at the Aladga station in the southeastern part of the Caspian Sea, whereas the tidal range in the western part of the South Caspian Basin varies from 5 to 10 cm.     

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.     

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

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

利用TOPEX/Poseidon高度计资料提取的太平洋M2内潮能通量分布

基于内波动力学方程,提出利用TOPEX/Poseidon高度计资料提取内潮的方法.利用该方法,结合1992年10月到2002年6月共10a的TOPEX/Poseidon高度计资料和Levitus(1998)资料,给出了整个太平洋M₂内潮能通量的分布,并与观测资料进行比测,两者符合较好.同时也发现沿整个太平洋边界M₂内潮能通量向大洋内部输入的总功率为58.4GW,其中北太平洋对此贡献为30.2GW,南太平洋为28.2GW,可见南、北太平洋的贡献是基本相等的.东太平洋的总量为17.8GW,而西太平洋为40.6GW,两者差异较大(以160°W作为东、西太平洋分界线).     

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.     

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

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

半日潮流作用下悬移质泥沙的运动特征及其影响因素研究

本文通过建立一维水深平均悬沙模型,对典型潮流控制的水道内悬沙运动特征进行研究。模型以泥沙再悬浮、沉降和平流为主要物理过程,动力因素包含M₂、S₂分潮及余流,采用湄洲湾2007年8月潮位、潮流、悬沙、底质同步观测资料进行分析和验证。通过三角傅里叶分析,将悬沙的时间序列分解为12个主要的谐波分量,其中主要分量包括:M₂分潮作用下产生的具有M₂倍潮角速度的1/4日分潮项,M₂与S₂分潮共同作用下且角速度为两分潮角速度之和的1/4日分潮项,及水平悬沙梯度、余流与M₂分潮共同作用下具有M₂分潮角速度的半日潮项。悬沙在时间上的平均值受到余流、悬沙水平梯度、M₂分潮流及悬沙起动条件等因素控制。余流导致了悬沙序列中相邻周期之间的不对称性。反映泥沙特性的参量对悬沙的曲线特征具有重要影响,泥沙沉降速度影响悬沙的相位,并影响其振幅;再悬浮有关的参量仅影响各谐波分量的振幅,但不影响相位。     

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

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

基于z坐标海洋模式的全球正压潮波数值模拟研究

The performance of a z-level ocean model, the Modular Ocean Model Version 4(MOM4), is evaluated in terms of simulating the global tide with different horizontal resolutions commonly used by climate models. The performance using various sets of model topography is evaluated. The results show that the optimum filter radius can improve the simulated co-tidal phase and that better topography quality can lead to smaller rootmean square(RMS) error in simulated tides. Sensitivity experiments are conducted to test the impact of spatial resolutions. It is shown that the model results are sensitive to horizontal resolutions. The calculated absolute mean errors of the co-tidal phase show that simulations with horizontal resolutions of 0.5° and 0.25° have about 35.5% higher performance compared that with 1° model resolution. An internal tide drag parameterization is adopted to reduce large system errors in the tidal amplitude. The RMS error of the best tuned 0.25° model compared with the satellite-altimetry-constrained model TPXO7.2 is 8.5 cm for M_2. The tidal energy fluxes of M_2 and K_1 are calculated and their patterns are in good agreement with those from the TPXO7.2. The correlation coefficients of the tidal energy fluxes can be used as an important index to evaluate a model skill.     

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

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

利用多源高度计资料提取南海低模态内潮能量

吕宋海峡由于剧烈变化的地形成为内潮产生的源地,内潮是海洋混合的重要原因。为了认知南海的内潮能通量分布,对南海的内潮有更好的理解,本文利用21世纪以来发射的多颗高度计卫星:J2、J1T、GFO以及EN,提取了吕宋海峡附近内潮的能通量。研究使用了调和分析和高通滤波等方法来提取第一模态内潮,主要提取K_1,K_2,M_2,N_2,O_1,P_1,Q_1和S_2八个分潮。同时结合WOA数据对能通量进行计算。结果表明,目标区域潮汐以全日分潮为主,所选区域的全日分潮中K_1所占比例最大;半日分潮中M_2分潮最强,而内潮的能通量则是M_2分潮所占最大,在吕宋海峡区域M_2能通量为6.45GW。内潮主要产生在地形变化剧烈的地方,海域的大部分地区内潮能量很小。在吕宋海峡中部,全日分潮能通量要小于南部地区,而半日分潮则有较大值。     

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.     

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.     

Modelling the tides and their impacts on the vertical stratification over the Sofala Bank,Mozambique

The Sofala Bank, a wide shelf located along the central coast of Mozambique, hosts tides with high amplitudes. The Regional Ocean Modelling System (ROMS) was used to analyse the tidal currents on the bank and to investigate their effects on the stratification and generation of tidal fronts. During spring tides, barotropic tidal currents with maximum values ranging from 40 cm s^–1 to 70 cm s^–1 are found on the central bank. The major axis of the tidal ellipses for M₂ and S₂ follow a cross-shelf direction with mainly anticlockwise rotation. Similar to observations, three distinct regimes occur: (i) a warm well-mixed region on the inner shelf where the depths are <30 m; (ii) a wellmixed colder region above the shelf edge; and (iii) a stratified region offshore. The model shows that the tides lead to cooling where two criteria are satisfied: the Simpson and Hunter parameter log₁₀(h/U³) <3.2 and the depth h >30 m. The shelf edge of the bank is important for internal tide generation. Two frontal structures result, one offshore between cooler mixed waters and warmer stratified waters and the other in shallow inshore waters, between cooler mixed waters and solar heated mixed waters.