抛物线型背景流对内孤立波特征和能量影响的数值模拟研究

本文基于南海东北部观测的抛物线型背景流,设计了8种形式的抛物线型背景流,利用IGW模式研究了其对内孤立波(ISW)的特征和能量学的影响。研究结果表明:背景流对波包中ISW数目没有影响,但减小了ISW的相速度;对于下边界在主温跃层附近或在其上的抛物线型背景流,ISW振幅和最大位移深度均增加;随着抛物背景流曲率减小,ISW振幅、斜压与正压能比值减小,同时ISW相速度、正压能、斜压能、KE/APE都增加;如果抛物背景流底部延伸至海底,且曲率减小,则ISW振幅、相速度减小,同时正压能、斜压能、KE/APE增加;在整个深度上的背景流,随着下层曲率减小和上层曲率增大,ISW振幅、相速度、斜压与正压能比值、斜压能、KE/APE均增加。     

Numerical simulation of interaction between internal solitary waves and submerged ridges

Direct numerical simulations are performed to study the transformation of internal solitary waves (ISWs) of depression type propagating over an underwater ridge in a two-layer fluid system. Bottom ridges with relatively smooth vertex are employed to represent sills in natural lakes and oceans. Consistent with previous experiments, three interaction types (weak, moderate and strong) are observed to be based on the energy loss. In addition, the moderate interaction are found to be categorized into transmitted and reflected type according to their distinct transformation process. General flow characteristics for ISW–ridge interaction in the benthic boundary layer and in the pycnocline is monitored and analysed. A modified degree of blocking B_m considering both the nonlinear effect of incident ISWs and the blockage effect of the submerged ridge is proposed. Different ISW–ridge interactions are discovered to be linked with B_m. Maximum wave-induced velocities, wave energy losses, reflected and transmitted wave amplitudes are found to have a self-similar feature with B_m. The maximum energy loss is up to 35% and the maximum wave-induced velocity can reach 1.8 times of the phase speed of the incident ISW. Empirical equations are obtained based on the data fitting to predict some useful physical parameters during ISW–ridge interaction.     

A new method to estimate phase speed and vertical velocity of internal solitary waves in the South China Sea

A new method of estimating the phase speed and vertical velocity of internal solitary waves (ISW) based on dynamic governing equations using Acoustic Doppler Current Profiler (ADCP) data was developed. This method was applied to a representative ISW case, which was captured in the South China Sea in the 2007 summer experiment. The result shows that this ISW had a phase speed of approximately 2.6?m?s^?1. This method provides an excellent alternative in estimating the phase speed of ISWs, especially useful in the absence of vertical stratification needed in solving the Korteweg–de Vries (KdV)-type equation, e.g., long-term mooring observation. Analyses show that the new method is fairly self-consistent, and it can be applied when only a part of the ISW observations is available. The computed vertical velocity of the ISW using the new method has a good agreement with the ADCP observations both in magnitude and pattern.     

Internal solitary waves in the East China Sea

A European Space Agency＇ s ENVISAT advanced synthetic aperture radar （ASAR） image covering Zhejiang coastal water in the East China Sea （ECS） was acquired on 1 August 2007. This image shows that there are about 20 coherent internal solitary wave （ISW） packets propagating southwestward toward Zhejiang coast. These ISW packets are separated by about 10 kin, suggesting that these ISWs are tide-generated waves. Each ISW packet contains 5-15 wave crests. The wavelengths of the wave crests within the ISW packets are about 300 m. The lengths of the leading wave crests are about 50 km. The ISW amplitude is estimated from solving KdV equation in an ideal two-layer ocean model. It is found that the ISW amplitudes is about 8 m. Further analysis of the ASAR image and ocean stratification profiles show that the observed ISWs are depression waves. Analyzing the tidal current finds that these waves are locally generated. The wavelength and amplitude of the ECS ISW are much smaller than their counter- parts in the South China Sea （SCS）. The propagation speed of the ECS ISW is also an order of magnitude smaller than that of the SCS ISW. The observed ISWs in the ECS happened during a spring tide period.     

Monthly variation of some parameters about internal solitary waves in the South China sea

In this paper, by non-dimensional analysis, it is found that finite-depth theory is more appropriate to the study of internal solitary waves (ISWs) in the South China Sea (SCS) than shallow-water theory. The 1-degree grid data of monthly mean temperature and salinity data at standard levels in the SCS are used to solve the linearized vertical eigenvalue problem. The nonlinear parameter and the wave phase speed are computed, then the nonlinear phase speed and the characteristic half-width of ISWs are calculated respectively by two different theories to investigate the difference between these two parameters in the SCS. The nonlinearity is the strongest near the continental slope of the SCS or islands where the bottom topography changes sharply, it is stronger in summer than that in winter; it increases (decreases) as pycnocline depth deepens (shallows), stratification strengthens (weakens) and pycnocline thickness thins (thickens). The nonlinear wave phase speed and the characteristic half-width are the largest in deep sea area, they then reduce peripherally in shallower water. The nonlinear wave phase speed in the SCS changes slightly with time, but the characteristic half-width changes somewhat larger with time. In most of the SCS basin, the nonlinear wave phase speed derived from shallow-water theory is very close to that derived from finite-depth theory, but the characteristic half-width derived from shallow-water theory is about 0.2–0.6 times larger than that derived from finite-depth theory. The ISW induced horizontal current velocity derived from shallow-water theory is larger than that derived from finite-depth theory. Some observed and numerical modeled ISW characteristic half-widths are compared with those derived from shallow-water and finite-depth theories, respectively. It is shown that, the characteristic half-widths derived from finite-depth theory agree better with observational and numerical modeled results than those derived from shallow-water theory in most cases, finite-depth theory is more applicable to the estimation of ISW characteristic half-widths in the northern SCS. It is also suggested that, to derive the precise ISW parameters in further study, the physical non-dimensional ratios which are related with ISW characteristic half-width, amplitude, thermocline and water depths should be calculated, so that an appropriate theory can be chosen for estimation.     

A modeling study of acoustic propagation through movingshallow-water solitary wave packets

Propagation of 400-Hz sound through continental-shelf internal solitary wave packets is shown by numerical simulation to be strongly influenced by coupling of normal modes. Coupling in a packet is controlled by the mode coefficients at the point where sound enters the packet, the dimensions of the waves and packet, and the ambient depth structures of temperature and salinity. In the case of a moving packet, changes of phases of the incident modes with respect to each other dominate over the other factors, altering the coupling over time and thus inducing signal fluctuations. The phasing within a moving packet varies with time scales of minutes, causing coupling and signal fluctuations with comparable time scales. The directionality of energy flux between high-order acoustic modes and (less attenuated) low-order modes determines a gain factor for long-range propagation. A significant finding is that energy flux toward low-order modes through the effect of a packet near a source favoring high-order modes will give net amplification at distant ranges. Conversely, a packet far from a source sends energy into otherwise quiet higher modes. The intermittency of the coupling and of high-mode attenuation via bottom interaction means that signal energy fluctuations and modal diversity fluctuations at a distant receiver are complementary, with energy fluctuations suggesting a source-region packet and mode fluctuations suggesting a receiver-region packet. Simulations entailing 33-km propagation are used in the analyses, imitating the SWARM experiment geometry, allowing comparison with observations     

内孤立波与平顶海山作用的能量耗散实验研究

内波破碎引起的能量耗散和混合是海洋内部的重要物理过程。通过在二维内波水槽进行实验室实验,分析内波与地形的作用,探究内孤立波与平顶海山地形作用时波要素、能量以及湍耗散率的时空变化。本实验利用重力塌陷法在两层流体中制造第一模态内孤立波,通过粒子图像测速技术(particle image velocimetry, PIV)获得内孤立波与地形作用时的流场结构,定量分析整个作用过程。结果表明,地形会改变波形甚至引起破碎,内波与地形作用时,振幅和能量密度会在内孤立波爬坡时迅速增大,在地形前缘产生强烈能量耗散。入射波的能量与塌陷高度呈二次函数关系,透射波能量随地形升高减小,反射波能量随地形升高增大。地形前缘局地湍耗散率极值时间序列在部分实验中呈双峰结构,对应内孤立波界面处剪切加强引起湍流耗散和波后缘翻转破碎。破碎引起的地形前缘区域平均湍耗散率量级在10~(-5)m~2/s~3,局地湍耗散率极值与入射波振幅呈指数关系,所有实验中局地湍耗散率的最大值接近10~(-3) m~2/s~3量级。     

喀拉海峡潮-地相互作用激发非线性内波的数值模拟

Nonlinear internal waves(NIWs) are ubiquitous around the Kara Sea, a part of the Arctic Ocean that is north of Siberia. Three hot spot sources for internal waves, one of which is the Kara Strait, have been identified based on Envisat ASAR. The generation and evolution of the NIWs through the interactions of the tide and topography across the strait is studied based on a nonhydrostatic numerical model. The model captures most wave characteristics shown by satellite data. A typical inter-packets distance on the Barents Sea side is about 25 km in summer, with a phase speed about 0.65 m/s. A northward background current may intensify the accumulation of energy during generation, but it has little influence on the other properties of the generated waves. The single internal solitary wave(ISW) structure is a special phenomenon that follows major wave trains, with a distance about 5–8 km. This wave is generated with the leading wave packets during the same tidal period. When a steady current toward the Kara Sea is included, the basic generation process is similar, but the waves toward the Kara Sea weaken and display an internal bore-like structure with smaller amplitude than in the control experiment. In winter, due to the growth of sea ice, stratification across the Kara Strait is mainly determined by the salinity, with an almost uniform temperature close to freezing. A pycnocline deepens near the middle of the water depth(Barents Sea side), and the NIWs process is not as important as the NIWs process in summer. There is no fission process during the simulation.     

Circulation and melting under the Ross Ice Shelf: estimates from evolving CFC,salinity and temperature fields in the Ross Sea

Temperature, salinity and chlorofluorocarbons (CFCs) 11, 12 and 113 were measured on a line of stations along the front of the Ross Ice Shelf in the austral summers of 1984, 1994 and 2000. Water mass distributions were similar each year but with high variability in the cross-sectional areas. CFC concentrations increased and salinity decreased with time throughout the water column. CFC saturation levels in the shelf and surface waters also increased with time and ranged from 43% to 90%. The undersaturation was due to inflow of low-CFC modified Circumpolar Deep Water, gas exchange limited by sea ice cover and isolation of water from the atmosphere beneath the ice shelf. The residence time of dense shelf waters resulting from sea ice formation is less well constrained by the chemical data than is the strong flow into the Ross Ice Shelf cavity. Shelf waters are transformed over about 3.5 years, by net basal melting of the ice shelf, into fresher Ice Shelf Water (ISW), which emerges as a large plume near the central ice front at temperatures below the sea surface freezing point. We estimate an average ISW production rate of 0.86 Sv and an average net basal melt rate of 60 km³/year for the Ross Ice Shelf exceeding a 300 m draft (75% of the ice cavity) during recent decades from box and stream tube models fit to all of the CFC and salinity data. Model fits to the individual data sets suggest ISW production and net basal melt rate variability due to interannual changes on a shorter time scale than our observations. ISW production based on the CFC budget is better constrained than net basal melting based on thermohaline data, with a heat budget yielding a rate of only 20 km³/yr. Reconciling differences between apparent freshwater and temperature changes under the ice shelf involves considerations of mixing, freezing and the flow of meltwater across the ice shelf grounding line.     

南海北部陆架区内孤立波向岸传播过程研究

南海北部是全球海洋中内孤立波最强和最为活跃的海域。然而,内孤立波在传入陆架区后,其形态发生显著变化,其传播演变过程表现出高度的复杂性。本研究综合卫星图像和数值模式手段研究了内孤立波在向岸传播过程中的空间变化特征。可见光卫星图像研究结果显示,南海北部陆架区存在三种形态的内孤立波,分别为第一模态下凹型内孤立波、第一模态上凸型内孤立波和第二模态内孤立波。受水深和层结变化的控制,它们的分布区域显著不同。基于MITgcm的数值模拟研究表明,上凸型内孤立波由第一模态下凹内孤立波经过极性转换过程发展而来,而第二模态内孤立波由第一模态下凹内孤立波与急剧变浅地形相互作用而产生。     

Assessment of theoretical approaches to derivation of internal solitary wave parameters from multi-satellite images near the Dongsha Atoll of the South China Sea

This study assesses the accuracy and the applicability of the Korteweg-de Vries(KdV) and the nonlinear Schr?dinger(NLS) equation solutions to derivation of dynamic parameters of internal solitary waves(ISWs) from satellite images. Visible band images taken by five satellite sensors with spatial resolutions from 5 m to 250 m near the Dongsha Atoll of the northern South China Sea(NSCS) are used as a baseline. From the baseline, the amplitudes of ISWs occurring from July 10 to 13, 2017 are estimate...     

青岛外海潮致第二模内孤立波的观测研究

根据在青岛附近海域开展的现场观测数据介绍了一种可能的新的第二模内孤立波的产生机制。温度链观测数据显示在2 h内孤立波特征由显著的第一模内孤立波变为了显著的第二模内孤立波。所观测到的第二模内孤立波的波面起伏和它的垂直结构同KdV(Korteweg-de Vries)方程的理论结果吻合良好,对应涨潮引起的内波非线性系数、频散系数以及Ursell数的变化和第二模内孤立波的产生相匹配。这些结果表明,所观测到的第二模内孤立波可能是由于涨潮导致的局地层结快速变化引起的。进一步分析表明,局地层结快速变化可能由潮汐推动黄海冷水团边缘锋面移动引起。     

一个典型南海北部第二模态内孤立波的观测分析

第二模态内孤立波在海洋中极少被观测到。本文基于潜标高时空分辨率观测数据,对南海北部陆架区的一个典型第二模态内孤立波进行了分析。结果表明,该第二模态内孤立波的流核出现在135 m深度处,其最大水平流速为0.66 m/s,传播方向为西偏北58°。沿传播方向的内孤立波流速分布在80～170 m的深度范围内,而与传播方向相反的逆流出现在海表和海底附近。垂向模态分析表明,该第二模态内孤立波水平流速的垂向结构与理论结果吻合良好。能量计算结果显示其动能密度的垂向积分可达14 kJ/m2,而波峰线方向单位长度上的动能估算值为5.98 MJ/m。尽管该第二模态内孤立波的动能比陆架区第一模态内孤立波小1个量级,但其高达0.045 s-1的流速垂向剪切约为典型第一模态内孤立波的2倍,表明其导致的混合可能更强。     

海洋内孤立波预警监测识别技术及其在流花16-2油田群开发中的应用

我国南海内孤立波发生频繁,由于其振幅大、流速强、能量集中等特点,已成为深海油气开发海上安装工程必须考虑的重大风险之一。依托现场守护船,提出一种内波流结构单体式监测、识别技术,采用KdV方程实现内波流预警,并成功应用于流花16-2油田群开发项目海管管线终端(PLET)安装、浮式生产储卸油装置(FPSO)锚系及水下管缆回接、空气潜水等高风险作业的内波流预警安全保障中,同时基于此次内波流预警监测结果对该海域的内波生成源和生成机制进行了探讨。应用期间采用该技术共监测、识别并发布内波预警450次。应用结果表明,约93%的内波预警时刻误差在10 min以内,平均误差为±3.90 min。统计分析结果表明,该海域的内孤立波主要来源于巴士海峡而非局地的潮地相互作用,大致可以分为两类:一类是直接产生于巴士海峡,西传至该海域;另一类是局地产生,主要是由巴士海峡产生的内波西传至陆坡浅化分裂生成或是由巴士海峡产生的内潮西传至陆坡由于非线性变陡机制产生。该内波流监测技术在流花16-2油田群开发项目的成功应用,可为南海其他类似项目提供直接的借鉴。     

A study of internal solitary waves observed on the continental shelf in the northwestern South China Sea

Based on in-situ time series data from the acoustic Doppler current profiler (ADCP) and thermistor chain in Wenchang area, a sequence of internal solitary wave (ISW) packets was observed in September 2005, propagating northwest on the continental shelf of the northwestern South China Sea (SCS). Corresponding to different stratification of the water column and tidal condition, both elevation and depression ISWs were observed at the same mooring location with amplitude of 35 m and 25 m respectively in different days. Regular arrival of the remarkable ISW packets at approximately the diurnal tidal period and the dominance of diurnal internal waves in the study area, strongly suggest that the main energy source of the waves is the diurnal tide. Notice that the wave packets were all riding on the troughs and shoulders of the internal tides, they were probably generated locally from the shelf break by the evolution of the internal tides due to nonlinear and dispersive effects.     

Spatial and temporal variation process of seabed dynamic response induced by the internal solitary wave

Internal solitary wave(ISW) is often accompanied by huge energy transport, which will change the pore water pressure in the seabed. Based on the two-dimensional Biot consolidation theory, the excess pore water pressure in seabed was simulated, and the spatiotemporal distribution characteristics of excess pore water pressure was studied. As the parameters of both ISW and seabed can affect the excess pore water pressure, the distribution of pore water pressure showed both dissipation and phase lag...     

南海东沙海域内孤立波形态研究

本文基于反射地震数据和MODIS遥感数据,对南海东沙海域内孤立波及孤立波群的形态有了系统的认识。内孤立波存在上升型和下降型两种极性波,又根据波形分成“钟形”、“平底形”和“碗形”三种类型,孤立波在波谷处的形状与孤立波振幅有关。单个内孤立波在传播一段时间后,受到各种因素的影响,会发育成波群。东沙海域的不同位置存在两种波群:“有序型”波群和“复杂型”波群,结合它们的位置及波群传播过程,认为这两种波群可能是孤立波从深海向陆架的整个传播过程中的两个阶段,“有序型”波群在被东沙岛阻碍后,受到各种海底地形、东沙环礁、波-波相互作用的影响,转变成“复杂型”波群。     

Experimental study of wave–wave nonlinear interactions using the wavelet-based bicoherence

The wavelet-based bicoherence, which is a new and powerful tool in the analysis of nonlinear phase coupling, is used to study the nonlinear wave–wave interactions of breaking and non-breaking gravity waves propagating over a sill. Two cases of mechanically generated random waves based on Jonswap spectra are used for this purpose. Values of relative depth, k_ph (k_p is the wave number of the spectral peak and h is the water depth) for this study range between 0.38 and 1.22. The variations of wavelet-based total bicoherence for the test cases indicate that the degree of quadratic phase coupling increases in the shoaling region consistent with a wave profile that is pitched shoreward, relative to a vertical axis as seen in the experiments, but decreases in the de-shoaling region. For the non-breaking case, the degree of quadratic phase coupling continues to increase until waves reach the top of the sill. Breaking waves, however, achieve their highest level of quadratic phase coupling immediately before incipient breaking and the degree of phase coupling decreases sharply following breaking. In addition the wavelet-based bicoherence spectra provide evidence of the harmonics' growth which is reflected in the energy spectra. The bicoherence spectra also show that quadratic phase coupling between modes within the peak frequency as well as between modes of the peak frequency and its higher harmonics are dominant in the shoaling region, even though there are relatively high levels of quadratic phase coupling occurring between other frequencies. Furthermore, using the temporal resolution property of the wavelet-based bicoherence, we find that the quadratic wave interactions occur more readily during segments of time with large change of wave amplitude, rather than those segments having large wave amplitudes, but small gradients in amplitude.     

基于大涡模拟耦合模式的小尺度海气相互作用研究

大气和海洋是影响地球气候系统的两个重要因素，它们之间的相互作用是海洋和大气研究的重要课题，海气耦合模式则是研究海气相互作用的重要工具，而海气耦合模式重点考虑的参数是海气通量。针对传统的大尺度海气耦合模式缺少湍流尺度分析的问题，本文使用并行大涡模拟海气耦合模式（The Parallelized Large-Eddy Simulation Model，PALM），在小尺度上探究风速对海气通量及湍流动能收支（Turbulence Kinetic Energy Budget，TKE Budget）的影响，设置了5、10和15m/s三种地转风速度对大气边界层（Atmospheric Boundary Layer，ABL）和海洋混合层（Oceanic Mixed Layer，OML）进行海气耦合模拟。研究表明：海气通量的分布与风速大小密切相关，风速越大，净热通量和浮力通量相对越大，由于温度上升导致海水蒸发加剧，使得大气的淡水通量增大；海洋湍流动能收支各项在近海面处受风速影响较大，且随着深度加深而逐渐减弱。本研究初步展示了小尺度海气耦合模式在海气通量研究中的应用，对进行小尺度海气相互作用研究具有一定的意义。     

An investigation on internal solitary waves in a two-layer fluid: Propagation and reflection from steep slopes

Experimental investigations on internal solitary wave (ISW) propagation and their reflection from a smooth uniform slope were conducted in a two-layered fluid system with a free surface. A 12-meter-long wave flume was in use which incorporated with: (1) a movable vertical gate for generating ISW; (2) six ultrasonic probes for measuring the fluctuation of an ISW; and (3) a steep uniform slope (from one of θ=30°, 50°, 60°, 90°, 120° and 130°) much greater than those ever published in the literature. This paper presents the wave profile properties of the ISW recorded in the flume and their nonlinear features in comparison with the existing Korteweg de Vries (KdV) and modified Korteweg-de Vries (MKdV) theories. Experimental results show that the KdV theory is suitable for most small-amplituded ISWs and MKdV theory is appropriate for the reflected ISWs from various uniform slopes. In addition, both the amplitude-based reflection coefficient and reflected energy approach a constant value asymptotically when plotted against the slope and the characteristic length ratio, respectively. The reflected wave amplitudes calculated from experimental data agree well with those reported elsewhere. The optimum reflection coefficient is found within the limit of 0.85 for wave amplitude, among the test runs from steep normal slope of 30° to inverse angle of 130°, and around 0.75 for the reflected wave energy, produced by an ISW on a vertical wall.