海洋内波发展演变数值试验

利用二维非静力数值模式,对海洋跃层中剪切背景流下内波做了线性和非线性模式的数值试验,结果表明:在本文所取的初始场和环境参数下,海洋内波密度扰动大值中心始终出现在密度跃层附近,被跃层所俘获,这与实际观测相一致;线性情况下,流函数扰动中心与密度扰动中心有很好地配合,密度扰动正、负中心附近分别存在流函数扰动的正、负单圈环流中心,且单圈环流垂直贯穿整个水体。线性模式中内波扰动强度在调整适应期后一直呈指数增长发展,剪切背景流是不稳定的;而非线性情况则呈准线性增长发展,以后增长停止而进入稳定期,且后者的增长也比前者要慢;非线性效应抑制了内波增长而使之趋于稳定,具有维稳作用。随着积分时间增加,线性模型中内波波包波形基本不变。波包的负、正振幅大体相同,整个波包宽度也基本不变,仅振幅有所增长;内波波包传播速度与其振幅无关,且传播速度基本不变。非线性情况下内波波包的形状则发生改变,其负、正振幅之比越来越大,波峰前的陡降也越来越强;并且内波传播速度与其振幅和波形有关,通常非线性效应越强,内波传播越慢。从本文中海洋内波的波长、周期和传播速度来看,该内波属于高频内波(短周期内波),性质为内重力惯性波。     

浅海内波及声场起伏数值研究

为探讨声学方法监测(反演)海洋内波的可行性,分析ASIAEX2001亚洲海海洋声学实验定点温度链实验数据,结果表明：实验期间存在M2内潮波、周期10～30min非线性内波现象,并观测到双温跃层同相传播“双非线性内波”。分别对垂直和平行于内波传播方向传播声场进行了数值计算分析。垂直于内波传播方向的声传播计算结果表明：简正波波数时间演化函数的频谱与内波引起的声速剖面变化频谱基本一致;简正波系数(模值)时间演化函数的频谱与内波引起的声源位置处的声速演化函数的频谱基本一致;高阶简正波波数差随时间的变化要比低阶简正波波数差小约1个量级。平行于内波传播方向的声传播计算结果表明：内潮波主要导致相邻简正波间耦合;非线性内波包能够导致跨号甚至垮多号简正波间耦合。讨论了2种海洋内波声学监测方法。     

基于MODIS遥感影像的安达曼海内波特征参数分布及生成周期研究

本文基于2017年634幅MODIS影像分析了安达曼海3个典型区域的内波空间分布特征，定量统计了波峰线长度、波包面积等特征参数，利用射线追踪法探讨了内波的潜在激发源并推算了内波的生成周期。研究表明，安达曼海北部海域的内波空间尺度较小，前导波波峰线的平均长度约为107 km，平均波包面积约为1 860 km2，内波的传播方向主要为东向以及西南向。安达曼海中部海域内波前导波波峰线的平均长度约为133 km，平均波包面积约为3 503 km2，超过70%的内波沿东偏北方向传播。苏门答腊岛北部海域内波前导波波峰线的平均长度约为131 km，平均波包面积约为2 997 km2，内波的传播方向主要为东向、东北向及东南向。安达曼海共有7个潜在内波激发源，内波的生成时间间隔介于11.5～13 h，具有明显的半日周期特征。     

近海潮成内波波包结构的分析

在文献[1]中,采用波包的观点给出了确定内波包的方向和群速的方法。本文基于1984年7月在L₃站所作的87.5小时的观测资料,研究了内波包的结构。结果表明,沿单一方向传播的潮成内波包是不连续的,存在四个波包,其间间以三个噪声区。四个波包的平均传播方向和速度相应为-90°和50cm/s。波包的频率特征与等温线及上、下层平均流的变差流的频谱特性相一致。通过模式对比证明,我们的二层模式与孤立子理论的二层模式的物理意义完全一致,即两者都反映出,在沿波传播的方向上,潮成内波的切变效应与其波幅成正比,不过我们的模式更为适用。     

东海西部陆架海区潮内波特征的初步分析

我国对浅海内波现象的调查和研究已有二十多年的历史.六十年代初,中国科学院海洋研究所在东海进行过多次单船定点连续观测和三船同步定点连续观测,并用调和分析的方法对观测资料进行过分析和研究,计算过该海区潮周期内波的各要素.近年来,中国科学院声学研究所在黄海也进行过多次短期定点连续观测,并报道过与声场有关的高频段内波谱分析的结果,但没有涉及潮周期内波.     

海洋内波破碎问题的研究

从理论、观测、数值实验和实验室实验四个方面对国内外近20年来关于海洋内波破碎问题的研究成果进行了分析总结.数值实验和实验室实验表明:中高频内波破碎时,初始的不稳定是二维的,当最终有横向对流卷团形成时,能量开始大量耗散,这时不稳定发展成为三维的;从初始的二维不稳定到对流卷团的产生这一过程,到底是一个剪切不稳定过程,还是一个对流不稳定过程,或者是对流不稳定和剪切不稳定共同存在的一个过程,取决于海水的层化、地形、背景剪切流和内波的自身性质.现场曾观测到内孤立波破碎时存在的剪切不稳定过程,数值研究模拟出了内孤立波破碎时存在的对流不稳定过程.现有的海洋内波破碎判据主要是关于中高频海洋内波的.理论分析侧重于确定线性或弱非线性内波的破碎机制和破碎条件.     

内波远距离传播演化的实验研究

为了研究内波远距离传播过程中的演化规律,本文采用图像测速法(PIV)分别对内波近场和远场的速度场进行测量。实验中同时采用两台CCD相机对实验区域进行拍摄,根据实验结果对内波能量和垂向模态结构进行计算分析。实验结果表明,在近场区域生成的内波主要表现为内波射线结构。内波射线在经过海表面反射后,其能量在空间上出现非对称结构,能量在加强区域较减弱区增加约15%。在远场,内波射线结构不再清晰,内波主要表现为低模态内波结构。内波射线在反射时能量衰减显著,损失约为50%;低模态内波可以离开内波生成源地远距离传播,传播过程中能量损失较小,在远场传播过程中(第一模态内波半波长的距离)能量损失约20%。低模态内波的传播相速度介于垂向第一模态和第二模态相速度之间。     

基于多源遥感数据的日本海内波特征研究

日本海特殊的地理位置和复杂的地形使得该海域内波表征极为复杂，遥感是大范围观测内波的有效手段，已被广泛应用于内波的探测研究。本文利用MODIS、GF-1和ENVISAT ASAR遥感影像，开展了日本海内波特征研究。通过提取内波波峰线，生成了日本海内波空间分布图；获取了内波的波峰线长度和传播速度，并基于非线性薛定谔方程反演了内波振幅。研究结果表明，日本海内波分布范围宽广，不仅大陆架沿海区内波分布密集，深海盆地也探测到了大量内波；日本海北部45°N附近海域有少量内波出现，利用高分影像探测到朝鲜陆架浅海区有大量小尺度内波，大和海盆、大和隆起的西南部海域没有发现内波。日本海内波波峰线长达100多千米，深海区的传播速度大于1 m/s；浅海区内波振幅约10 m左右，深海区可达60 m以上。     

失稳背景流下对称型海洋内波的生成演变

本文采用无海底地形但考虑海洋跃层和剪切背景流的二维非静力准不可压缩方程组的数值模式，开展失稳垂向剪切背景流下线性和非线性对称型海洋内波生成演变的数值实验，并对结果进行分析、比较和讨论。研究结果表明，线性内波强度随积分时间始终呈指数增长，并有内波的对称不稳定；而非线性内波强度则在发展期呈准线性增长，最终进入稳定期。线性增长比非线性增长要快得多，非线性效应具有维稳作用。对该线性和非线性对称型内波，在跃层附近位密度扰动均有大值中心，即其为跃层所俘获，这与实际观测相一致；流函数与位密度扰动两者均有很好配合，位密度扰动的正、负中心分别相应于流函数的上升、下沉运动，表明有从海底向上的斜对流发生，且以跃层为顶盖。对线性内波来说，随积分时间增加，其波形大体不变，其正、负振幅也大体相同，并有符号相反原地增长的两个倾斜环流圈，而在它们之间则有较强倾斜上升流。非线性内波波形随积分时间改变，倾斜环流圈数目也在增加，最终形成负环流强于正环流的结果，并导致流函数、位密度扰动水平梯度剧增，其可视为间断。     

内波吸引子的数值模拟

密度稳定层结的流体中产生的内波沿着由内波固有频率、流体浮力频率等因素所确定的特征线(或内波射线)传播。边界上的反射不改变内波的频率,从而也不改变反射后的内波特征线与重力方向所成的夹角。侧边界倾斜的封闭容器内,内波能量沿特征线传播的过程中经侧壁、表面和底面的反射可能会集中在一个封闭的轨道上,形成内波吸引子。该现象已经得到水槽试验、线性理论和数值试验的验证。本文利用非线性非静压的环流模式MITgcm,模拟了二维封闭区域中(1,1)-吸引子和(2,1)-吸引子的形成过程,并讨论初值条件对它们的影响。稳定的(1,1)-吸引子其极限环两侧流速出现很强的剪切流。当减小地形的坡度时,由于线性因素的增加,吸引子的结构不变,但吸引子厚度在相空间中的收缩速度加快。对于(2,1)-吸引子,由于轨道所成的两个环中间的节点耗散了部分能量,吸引子的收敛速度较慢。节点处,流体速度始终为0,但存在强烈混合,流体浮力频率呈现振幅较大的周期变化。     

深海地形调查中声速剖面质量控制方法

通过对深海多波束地形调查中实测声速剖面数据质量的研究,探讨声速剖面仪检定以及实际应用中的质量控制措施。利用抛弃式温盐深仪和同海区国际剖面浮标的剖面数据,对声速剖面仪实测数据进行比对,分析声速剖面数据的固定偏差或线性偏差,完成声速剖面数据修正,提高了实测声速剖面数据的应用水平,为大洋专项调查中声速剖面数据的质量控制,提供了一种有效参考方法。     

Solitary wave propagation influenced by submerged breakwater

The form of Boussinesq equation derived by Nwogu （1993） using velocity at an arbitrary distance and surface elevation as variables is used to simulate wave surface elevation changes. In the numerical experiment, water depth was divided into five layers with six layer interfaces to simulate velocity at each layer interface. Besides, a physical experiment was carried out to validate numerical model and study solitary wave propagation.“Water column collapsing”method （WCCM） was used to generate solitary wave. A series of wave gauges around an impervious breakwater were set-up in the flume to measure the solitary wave shoaling, run-up, and breaking processes. The results show that the measured data and simulated data are in good agreement. Moreover, simulated and measured surface elevations were analyzed by the wavelet transform method. It shows that different wave frequencies stratified in the wavelet amplitude spectrum. Finally, horizontal and vertical velocities of each layer interface were analyzed in the process of solitary wave propagation through submerged breakwater.     

Mooring observations of internal solitary waves in the deep basin west of Luzon Strait

A total of 137 internal solitary waves （ISWs） are captured during a field experiment conducted in the deep basin west of the Luzon Strait （LS） from March to August, 2010. Mooring observations reveal that a fully developed ISW owns a maximum westward velocity of more than 1.8 m/s and an amplitude of about 200 m. The ISWs in the South China Sea （SCS） are most active in July, which may be due to the strong stratification in summer. Most of the ISW episodes are detected around and after the 1st or 15th lunar day, indicating that the ISW in the SCS is triggered by astronomic tides. Half part of the observed ISWs were detected around 19：00 local time, which can be explained by the fact that type-a ISWs emerged in the evening at roughly the same time each day. The propagation direction of the ISWs and the astronomic tides in the LS show that the area south of the Batan Island is probably the main source region of the type-a ISWs, while the area south of Itbayat Island and south of the Batan Island is likely the main source region of the type-b ISWs observed at the mooring. Moreover, for the resonance of semidiurnal internal tides emitting from the double ridges in the LS, the underwater ridge south of the Itbayat Island and south of the Batan Island is believed to favor the generation of the energetic ISWs.     

Time-domain analysis of second-harmonic generation of primary Lamb wave propagation in an elastic plate

Within the second-order perturbation approximation, this paper investigates the physical process of generation of the time-domain second harmonic by a primary Lamb wave waveform in an elastic plate. The present work is performed based on the preconditions that the phase velocity matching is satisfied and that the transfer of energy from the primary Lamb wave to the double frequency Lamb wave is not zero. It investigates the influences of the difference between the group velocities of the primary Lamb wave and the double frequency Lamb wave, the propagation distance and the duration of the primary Lamb wave waveform on the envelope shape of the time-domain second harmonic. It finds that the maximum magnitude of the envelope of the second-harmonic waveform can grow within some propagation distance even if the condition of group velocity matching is not satisfied. Our analyses also indicate that the maximum magnitude of the envelope of the second-harmonic waveform is kept constant beyond a specific propagation distance. Furthermore, it concludes that the integration amplitude of the time-domain second-harmonic waveform always grows with propagation distance within the second-order perturbation. The present research yields new physical insight not previously available into the effect of generation of the time-domain second harmonic by propagation of a primary Lamb wave waveform.     

The Relationship between Propagation Time and Sound Velocity Profile for Positioning Seafloor Reference Points

This paper focuses on the relationship between the propagation time and the integral of the sound velocity profile (SVP) with respect to depth (SVP area for short) for positioning seafloor reference points. We proved a linear relationship between the propagation time and the SVP area and defined its mathematical expression (ST law). We showed in three simulations of possible variations in SVPs (random errors, internal waves, or a combination of both) that this ST law is verified. Using this new law in a simulated experiment, the ranging residuals due to sound velocity change through time are significantly improved. With this method, monitoring the SVP while acquiring travel times to a submerged transponder can significantly reduce the positioning errors of the transponder.     

Artificial Upwelling of Deep Seawater Using the Perpetual Salt Fountain for Cultivation of Ocean Desert

Deep seawater in the ocean contains a great deal of nutrients. Stommel et al. have proposed the notion of a “perpetual salt fountain” (Stommel et al., 1956). They noted the possibility of a permanent upwelling of deep seawater with no additional external energy source. If we can cause deep seawater to upwell extensively, we can achieve an ocean farm. We have succeeded in measuring the upwelling velocity by an experiment in the Mariana Trench area using a special measurement system. A 0.3 m diameter, 280 m long soft pipe made of PVC sheet was used in the experiment. The measured data, a verification experiment, and numerical simulation results, gave an estimate of upwelling velocity of 212 m/day. This revised version was published online in July 2006 with corrections to the Cover Date.     

Propagation Mechanisms of Incident Tsunami Wave in Jiangsu Coastal Area,Caused by Eastern Japan Earthquake on March 11,2011

At 13:46 on March 11, 2011(Beijing time), an earthquake of Mw=9.0 occurred in Japan. By comparing the tsunami data from Guanhekou marine station with other tsunami wave observation gathered from southeast coastal area of China, it was evident that, only in Guanhekou, the position of the maximum wave height appeared in the middle part rather than in the front of the tsunami wave train. A numerical model of tsunami propagation based on 2-D nonlinear shallow water equations was built to study the impact range and main causes of the special tsunami waveform discovered in Jiangsu coastal area. The results showed that nearly three-quarters of the Jiangsu coastal area, mainly comprised the part north of the radial sand ridges, reached its maximum tsunami wave height in the middle part of the wave train. The main cause of the special waveform was the special underwater topography condition of the Yellow Sea and the East China Sea area, which influenced the tsunami propagation and waveform significantly. Although land boundary reflection brought an effect on the position of the maximum wave height to a certain extent, as the limits of the incident waveform and distances between the observation points and shore, it was not the dominant influence factor of the special waveform. Coriolis force's impact on the tsunami waves was so weak that it was not the main cause for the special phenomenon in Jiangsu coastal area. The study reminds us that the most destructive wave might not appear in the first one in tsunami wave train.     

粤西阳江海域春夏季沿岸激流现象及机理分析

通过分析2018–2019年粤西阳江沿岸流海域多站点周年观测水文气象实测资料,发现2019年春季和夏季阳江沿岸流海域存在激流现象。研究结果表明:（1）2019年5月5日凌晨6时,观测站点2 m水深处流速达到164.7 cm/s,9 m水深处流速达到127.6 cm/s。2019年8月1日凌晨4时至5时,阳江沙扒海域2 m水深处流速达到161.8 cm/s,9 m水深处流速达到156.6 cm/s。（2）粤西沿岸流阳江20～30 m水深海域春夏季突发性强流具有典型的激流特征。激流在涨急时刻发生在海洋表层,持续2～4 h。（3）在西南风与东北风转换期间,粤西沿岸海域容易形成海水幅聚带,近岸海域海平面上升,外海海域海平面下降,强劲的自岸向外水平压强梯度力导致近岸海水加强向西运动,从而产生激流。     

Comparison of experimental irregular water wave elevation and kinematic data with new hybrid wave model predictions

Surface water wave elevations and kinematics from four unidirectional irregular wave trains, with a Pierson and Moskowitz or JONSWAP random wave spectrum, were measured in the laboratory using resistance wave probes and a laser Doppler anemometer. The wave elevation data, velocity time series, extreme (largest) wave horizontal velocity profiles and extreme wave acceleration fields are compared with the predictions of a new wave kinematics model, named the hybrid wave model. Irregular waves are commonly viewed as the summation of many linear wave components of different frequencies, but more accurate predictions of downstream surface elevations (wave evolution) and wave kinematics are attained by considering the non-linear interactions among wave components. The hybrid wave model incorporates these non-linear wave component interactions, and its wave evolution predictions and kinematics estimates are compared with laboratory measurements in this study. Linear random wave theory, Wheeler stretching and linear extrapolation wave kinematic prediction techniques are also compared. Comparisons between measurements and hybrid wave model estimates demonstrate its improved capability to predict velocity and acceleration fields and wave evolution in two-dimensional irregular waves.