莫桑比克海峡及其邻近海区正压潮流数值模拟与能量收支分析

莫桑比克海峡及其邻近海区是全球海洋潮流和潮能耗散最强的海区之一。文章利用高分辨率通用环流模式对该海区的正压潮流进行模拟, 并对该海区潮能通量和潮能耗散特征进行分析。结果表明, 莫桑比克海峡及其邻近海区的潮波主要是半日分潮占主导地位, 全日分潮可忽略不计, M₂分潮形成1个左旋潮波系统和1个右旋潮波系统, S₂分潮形成1个左旋潮波系统。莫桑比克海峡和马达加斯加岛南部等绝大数区域的M₂和S₂半日潮流是逆时针旋转, 在马达加斯加岛顶部等局部区域是顺时针旋转, 而且在海峡通道等复杂地形处潮流流速量级较大。潮能通量矢量主要来自东边界, 大部分潮能通量沿马达加斯岛北部传入莫桑比克海峡区域, 其中经过马达加斯加岛北部和进入莫桑比克海峡的M₂ (S₂)分潮的潮能通量分别为156.86GW (40.53GW)和148.07GW (36.05GW), S₂分潮潮能通量的量级大约为M₂分潮的1/5~1/4。底摩擦耗散主要发生莫桑比克海峡和马达加斯加岛南北部, 其中莫桑比克海峡M₂ (S₂)分潮的底摩擦耗散为1.762GW (0.460GW), 占其底部总耗散的43.74% (39.72%)。     

黄海、渤海TOPEX/Poseidon高度计资料潮汐伴随同化

首先将大约10a的TOPEX/Poseidon(T/P)高度计资料沿星下轨迹点做潮汐调和分析,提取得到各分潮的调和常数,利用伴随同化方法,同化到二维非线性潮汐数值模式中,模拟了黄海、渤海区域M₂,S₂,O₁,K₁等4个潮汐分潮,并根据计算结果给出了各分潮的同潮图.将计算值与观测值的进行偏差统计,结果表明计算值与验潮站资料符合良好.研究过程中做了两类试验:一类试验是针对不同的参数进行优化,一类试验是针对不同的资料进行同化.第一类试验表明:将开边界条件和底摩擦系数同时作为模型优化的控制参数,其结果明显优于单独优化开边界条件;第二类试验表明:同时同化高度计资料与验潮站资料,比单独同化其中任一种资料,对模式计算结果都有较好的改进.研究结果表明,采用伴随同化方法,利用T/P高度计资料和验潮站资料作为同化数据能有效改进模拟结果,用来反演黄海、渤海的潮波系统是可行的.     

基于近30 a多源卫星高度计数据的东中国海潮汐信息提取

我国HY-2B卫星已成功运行3 a多,本文首次将 HY-2B测高数据用于计算潮汐。将HY-2B与相同时间段的Jason-3在东中国海分别进行潮汐信息提取,验证了其结果的一致性。建立了基于 10 颗国内外卫星高度计(TOPEX/Poseidon、Jason-1、Jason-2、Jason-3、ENVISAT、ERS-1、ERS-2、Sentinel-3A、Sentinel-3B、HY-2B)数据的时间序列,得到东中国海近30 a时间序列的较高空间分辨率网格化海面高度。利用该数据提取了东中国海12个分潮(O₁、K₁、Q₁、P₁、M₂、S₂、N₂、K₂、S_A、S_Sa、M_m、M_f)的潮汐调和常数,并将 4 个主要分潮M₂、S₂、K₁、O₁的调和常数...     

太平洋潮波特征比较

本文对TOPEX/Poseidon高度计资料直接分析得到4个主要分潮(M2、K1、S2和O1)的调和常数,将其与全球大洋潮波模式TPXO6.2的模拟结果以及太平洋中48个验潮站观测资料的分析结果进行了系统比较,得出高度计资料直接分析结果与潮波模式模拟结果总体比较一致.模拟出的无潮点的位置和高度计资料直接分析结果有差异,且K1和O1分潮差异较大.与站点结果比较表明TPXO6.2模式模拟结果与验潮站观测结果的振幅绝对偏差小于2cm的站点的百分比达到85%,迟角绝对偏差小于40°的站点的百分比达到70%以上,比高度计资料直接分析结果精确;矢量均方根误差比较表明,太平洋中部结果比整个太平洋结果准确,太平洋矢量均方根误差的值和其他研究者、其他模式的结果近似.     

TOPEX/POSEIDON高度计资料用于中国近海M2分潮的同化反演

以整个中国海为研究区域,把开边界取在TOPEX/POSEIDON(T/P)高度计资料数据点,利用趋近插值法将高度计资料同化到二维非线性潮汐数值模式中,反演了中国近海M2分潮的分布.计算结果与60个验潮站资料比较,符合良好振幅平均绝对偏差为6.7cm,相角偏差为5.5°,均方根偏差为12.6cm.根据计算结果给出东中国海和南海M2分潮的同潮图,与相关文献比较,符合良好.研究结果表明,以本方法利用T/P资料来反演中国海潮汐是可行的.     

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

利用球坐标系中的二维非线性潮波方程组,数值计算了渤、黄、东海全海区的全日及半日潮汐潮流。沿岸81个潮位站的计算与实测值的比较表明,M₂分潮振幅差平均为7.2cm,相角差为6.4°,m₁分潮振幅差平均为2.6cm,相角差为7.4°,计算与实测符合良好。潮流的比较结果表明,计算与实测的符合程度也是比较好的。文中给出的同潮图同Fang(1986)给出的实测与数值的综合结果基本一致。本计算还证实或首次给出了若干圆流点。如对M₂分潮流,证实了在北黄海山东北部近海及南黄海北部各存在一对圆流点,并在浙江北部近海新发现一对圆流点;对m₁分潮流在苏北浅滩外侧发现一个圆流点,另外在东海东北部(济州岛东南)新给出两个圆流点,东海东南部的弱流区存在三个圆流点,此外,文中还分别讨论了M₂及m₁分潮能通量的传播和消耗情况,并指出从太平洋经吐噶喇海峡及冲绳至宫古岛之间的水道传入东海的m₁分潮,在遇到陆坡的阻挡后,其中有相当部分潮能被反射回太平洋。     

南海潮汐潮流的数值模拟

本文用二维球坐标数值模式计算了南海m₁[=(K₁+O₁)/2]和M₂分潮的分布.计算范围从2°N到25°N,99°E到121°30'E,坐标的经向纬向、格距均为1°/4.计算结果与92个实测站进行比较符合良好,m₁分潮振幅的平均误差为4cm,迟角为7°.M₂分潮振幅的平均误差为9cm,迟角为12°.根据计算结果给出南海m₁和M₂分潮的潮汐、潮流、潮余流和潮能通量分布图.     

渤黄东海潮能通量与潮能耗散

利用同化高度计资料和沿岸验潮站资料对潮汐数值模式进行同化,根据同化后的数值模式结果,对渤黄东海中的潮能通量和潮能耗散进行了研究.M2分潮从太平洋进入渤黄东海的潮能为122.499GW,占4个主要分潮进入总量的79%.黄海是半日分潮潮能耗散的主要海区.全日分潮则主要耗散在东海.全日分潮在遇到陆坡的阻挡以后有一部分潮能沿着冲绳海槽向西南传播,并有一部分潮能反射回太平洋,其中O1分潮通过C3断面反射回太平洋的潮能,约占其传入东海潮能的44%.     

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

利用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个分潮的潮汐、潮流、潮余流和潮能通量密度分布图.     

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

把利用正交潮响应方法对 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资料同化到潮汐数值模式中去对模式进行校正是有效的     

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.     

吕宋海峡M2内潮的三维数值模拟

A three-dimensional isopycnic-coordinate internal tidal model is employed to investigate the generation,propagation, vertical structure and energy conversion of M2 internal tides in the Luzon Strait(LS) with mooring observations. Simulated results, especially the tidal current amplitudes, agree well with observations,demonstrating the reasonability and accuracy of the model. Results indicate that M2 internal tides mainly propagate into three directions horizontally, i.e., eastward towards the western Pacific Ocean, westward towards the Dongsha Island and southwestward towards the South China Sea Basin. In the horizontal direction, tidal current amplitudes decrease as distance increases away from the LS; in the vertical direction, they show an obvious decreasing tendency with depth. Between the double ridges of the LS, a clockwise gyre of M2 baroclinic energy flux appears, which is caused by reflections of M2 internal tides at supercritical topographies, and resonance of M2 internal tides happens along 19.5° and 21.5°N due to the heights and separation distance of the double ridges. The total energy conversion in the LS is about 14.20 GW.     

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

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

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

吕宋海峡由于剧烈变化的地形成为内潮产生的源地,内潮是海洋混合的重要原因。为了认知南海的内潮能通量分布,对南海的内潮有更好的理解,本文利用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。内潮主要产生在地形变化剧烈的地方,海域的大部分地区内潮能量很小。在吕宋海峡中部,全日分潮能通量要小于南部地区,而半日分潮则有较大值。     

吕宋海峡及周边海域湍流混合的时空分布特征研究

利用1992—2002年的温盐深数据与2012—2016年的Argo数据,基于细尺度参数化方法研究了吕宋海峡及周边海域(12°—30°N,115°—129°E)湍流混合的时空分布特征,并分析了地形粗糙度、内潮以及风输入的近惯性能通量对湍流混合的影响。结果表明,吕宋海峡和东海陆坡处具有强混合的特征,扩散率高达4×10~(-3) m~2/s,主要是由内潮产生导致的,其中吕宋海峡主要是M2、K1和O1内潮的贡献,而东海陆坡处主要是M_2内潮的贡献;南海北部也呈现较强的混合,且陆坡处的混合比海盆高1—2个量级;南海中央海盆和离岸的菲律宾海混合较弱,扩散率为O (10-5 m2/s)。此外,在研究区域内,湍流混合的年际变化和季节变化均不明显,且混合扩散率与风输入的近惯性能通量未表现出明显的季节相关。     

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.     

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.     

吐噶喇海峡内潮的能量收支和大小潮变化

吐噶喇海峡是西北太平洋重要的内潮产生区域,该区域内产生的内潮对于东海陆架和西北太平洋的混合和物质输运有十分重要的作用。水平分辨率为3km的JCOPE-T(JapanCoastalOcean PredictabilityExperiment—Tides)水动力学模式的结果表明,吐噶喇海峡的内潮主要产生在地形变化剧烈的海山和海岛附近,其引起的等密面起伏振幅可达30m。吐噶喇海峡的内潮在垂直于等深线方向分为两支向外传播:一支向西北方向传播,进入东海陆架后迅速减小;另一支向东南方向传播,进入西北太平洋。吐噶喇海峡潮能丰富,其在约半个月内的平均输入的净正压潮能通量为13.92GW,其中约有3.73GW转化为内潮能量。生成的内潮能量有77.2%在当地耗散,传出的内潮能通量为0.84GW,主要通过西北和东南两个边界传出。该区域潮能通量有显著的大小潮变化,大潮期间输入的正压潮净能通量和产生的内潮能通量均约为小潮期间的2倍,但其主要产生区域基本不变,且内潮能量耗散比率均在产生的内潮通量的76%—79%。另外,内潮能通量的传播方向也没有发生变化,仍主要通过西北和东南两个边界传出。因此,大小潮的变化仅影响吐噶喇海峡处产生的内潮能量的大小,不影响其产生区域、传播方向和耗散比率。     

Mode-1 internal tides in the western North Atlantic Ocean

Mode-1 internal tides were observed the western North Atlantic using an ocean acoustic tomography array deployed in 1991–1992 centered on 25°N, 66°W. The pentagonal array, 700-km across, acted as an antenna for mode-1 internal-tides. Coherent internal-tide waves with O(1 m) displacements were observed traveling in several directions. Although the internal tides of the region were relatively quiescent, they were essentially phase locked over the 200–300 day data record lengths. Both semidiurnal and diurnal internal waves were detected, with wavenumbers consistent with those calculated from hydrographic data. The M₂ internal-tide energy flux was estimated to be about 70 W m⁻¹, suggesting that mode-1 waves radiate 0.2 GW of energy, with large uncertainty, from the Caribbean island chain at this frequency. A global tidal model (TPXO 5) suggested that 1–2 GW is lost from the M₂ barotropic tide over this region, but the precise value was uncertain because the complicated topography makes the calculation problematic. In any case, significant conversion of barotropic to baroclinic tidal energy does not occur in the western North Atlantic basin. It is apparent, however, that mode-1 internal tides have very weak decay and retain their coherence over great distances, so that ocean basins may be filled up with such waves. Observed diurnal amplitudes were an order of magnitude larger than expected. The amplitude and phase variations of the K₁ and O₁ constituents observed over the tomography array were consistent with the theoretical solutions for standing internal waves near their turning latitude. The energy densities of the resonant diurnal internal waves were roughly twice those of the barotropic tide at those frequencies.     

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.