波生运动和波湍相互作用对海洋上混合层影响的数值探究

本文将波生运动和波湍相互作作用的参数化方案嵌入一维垂向混合模式GOTM中,并与不考虑波浪效应以及仅考虑波浪破碎的试验结果进行对比,发现不考虑波浪效应时,海表温度模拟结果偏高,混合层深度模拟结果偏浅,偏差在夏季尤其明显。波浪破碎对湍流的增强作用仅限于上层几米甚至仅限于表层,对整个混合层的温度分布和混合层深度影响不大。波生运动和波湍相互作用则有效增强海洋上层的湍流强度,改善模式高估海表温度而低估混合层深度的问题,温度分布模拟结果降低了上层温度同时增大了次表层温度,与观测更加相符。波生运动和波湍相互作用增大了海洋上层的湍流剪切生成项、湍动能、耗散率和湍流输运系数,两者对上混合层的温度分布、湍流强度和湍流输运作用的改善结果十分相似。波生运动和波湍相互作用的影响在冬季并不明显,此时可能有其他混合机制占主导地位。     

夏季长江口及其邻近海域湍流特征分析

利用2019年7月在长江口科学考察实验研究夏季航段(NORC2019-03-02)中获得的MSS90L湍流剖面仪的直接观测数据,本文计算并分析了该断面的湍动能耗散率ε和垂向湍扩散系数K_Z的分布情况。湍动能耗散率的大小为1.72×10?10～2.95×10?5 W/kg;垂向湍扩散系数的大小为3.24×10?7～4.55×10?2 m2/s。湍动能耗散率和垂向湍扩散系数的分布相似,均为上层最强,底层次之,中层最弱。上层由于风应力的作用,使得湍动能耗散率和垂向湍扩散系数较大;温跃层处层化较强,抑制了湍动能的耗散和垂向上的湍混合。盐度锋面的次级环流会促使低盐水团脱离,锋面引起的垂向环流会加强海洋的湍混合。低盐水团与外界的能量交换较少,湍动能耗散率较弱。长江口海区存在明显的上升流和下降流,它们是由锋面的次级环流产生的;上升流和下降流的存在促进湍动能的耗散与湍混合。     

海浪破碎对海洋上混合层中湍能量收支的影响

海浪破碎产生一向下输入的湍动能通量,在近海表处形成一湍流生成明显增加的次层,加强了海洋上混合层中的湍流垂向混合。为了研究海浪破碎对混合层中湍能量收支的影响,文中分析了海浪破碎对海洋上混合层中湍流生成的影响机制,采用垂向一维湍封闭混合模式,通过改变湍动能方程的上边界条件,引入了海浪破碎产生的湍动能通量,并分别对不同风速下海浪破碎的影响进行了数值研究,分析了混合层中湍能量收支的变化。当考虑海浪破碎影响时,近海表次层中的垂直扩散项和耗散项都有显著的增加,该次层中被耗散的湍动能占整个混合层中耗散的总的湍能量的92.0%,比无海浪破碎影响的结果增加了近1倍;由于平均流场切变减小,混合层中的湍流剪切生成减小了3.5%,形成一种存在于湍动能的耗散和垂直扩散之间的局部平衡关系。在该次层以下,局部平衡关系与壁层定律的结论一致,即湍动能的剪切生成与耗散相平衡。研究结果表明,海浪破碎在海表产生的湍动能通量影响了海洋上混合层中的各项湍能量收支间的局部平衡关系。     

风作用下孤立波在珊瑚礁上传播变形机制数值模拟研究

基于二维不可压缩两相流模型建立了数值风浪水槽,采用SST k-ω雷诺时均湍流模型,研究了风作用下孤立波在珊瑚礁上的传播变形规律。将计算结果与实验数据对比,证明了该两相流模型计算孤立波在珊瑚礁上传播的准确性,并进一步分析了不同风速对珊瑚礁上孤立波传播变形的影响。结果表明：风的作用会使波面发生随机脉动特征。当地波高随风速的增大而增大;当地波高关于风速的变化梯度随入射波高的增大而增大。风的作用会加快孤立波的传播并且使孤立波提前发生破碎;孤立波开始破碎的位置随风速的增大向远离礁坪的方向移动。反射系数随风速的增大而增大;反射系数关于风速的变化梯度随入射波高的增大而减小;透射系数随风速的增大呈增大趋势。平底区波峰剖面同一水深处的水平流速随风速的增大而增大;且一定的风速不改变水平流速沿水深的变化梯度。有风时波面上方的矢量密度和大小均明显高于无风时且与风速呈正相关,并且波峰上方气流不再循环。随着风速的增大,水气交界面附近的正涡量和负湍流剪应力减小,负涡量和正湍流剪应力增大。水体动能、势能和总能达到高值的时间随风速的增大而减少;水体动能、势能和总能随风速的增大而增大,并且风速对水体动能的相对影响大于势能。     

内孤立波破碎所致混合的实验研究

为定量分析内孤立波破碎的混合过程,本文在二维内波水槽中进行了两层流体第一模态内孤立波在斜坡上破碎的实验,运用粒子图像测速技术(PIV)测量内孤立波传播、破碎、反射过程的流场,计算涡度、湍动能和湍耗散率。结果表明不同振幅内波在不同角度斜坡上破碎时各个量的分布特征十分相似,各组实验各要素时间序列中均有两个峰值,分别发生于非线性增强和破碎时刻。得到破碎时湍耗散率与内孤立波振幅的关系为:较小振幅内波的湍耗散率与振幅呈2次关系,无因次振幅增大到0.9湍耗散率趋于不变;与斜坡角度的关系为:对于小振幅内波斜坡角度增大,破碎程度降低,耗散率减小;振幅较大时,存在一个角度使破碎程度最大。破碎引起的湍耗散率的量级在10–7到10–4m2/s3之间,比实测海洋中内孤立波传播界面和内潮遇地形破碎的湍耗散大1个量级。     

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

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

渤海弱层化期湍流混合观测特征分析

利用2017年9月在渤海共享航次中取得的湍流混合直接观测数据,本文研究了渤海海域湍流混合的空间分布特征及有关的影响因素。9月观测海区水体垂向层结较弱,莱州湾受黄河冲淡水影响出现高温低盐结构,位于渤海中央浅滩南北两侧洼地的双中心冷水结构依旧存在。湍流观测结果表明湍动能耗散率在10~(-9)～10~(-5)W/kg之间变化,统计上满足对数正态分布。耗散率强值区出现在辽东湾及渤海湾湾口近岸处,相应的垂向湍扩散系数约为10~(-6)～10~(-2)m~2/s。垂向上,水体表、底层混合较强,进一步研究发现弱层化水体的平均湍动能耗散率〈ε〉与风速和正压潮流速的大小存在正相关关系。另一方面,耗散率ε与浮性频率N近似满足ε=2.0×10~(-8)+3.0×10~(-7)(N~2/N_0~2)~(-5)的拟合函数关系,反映了层化对水体垂向混合的抑制作用。     

不同破碎波气液混合区发展特征及其运动特性研究

了解波浪破碎速度场的分布特性对于波浪破碎物理机制的研究极为重要,同时,对比研究不同类型破碎产生的气液混合区的演化特征有利于波浪白冠覆盖率模型的完善。在实验室水槽,生成了深水临界波、单次崩破波和单次卷破波,采用图像测速技术获取了波浪破碎图像、波面下水体和气液混合区速度场。结果表明,崩破波的水平向速度u和垂直向速度v在波峰前和波峰后的分布极为不对称,其水平向最大速度u_max并不位于峰顶,而是在主导波峰前0.7η_{m ax}处;卷破波的u_max出现在波峰峰顶前端极小的区域内,且该区域与周围区域的速度梯度极大。崩破波和卷破波生成的气液混合区发展特征也存在差异：崩破波的u_max值大、影响区域长、混合区厚度较小、各区域影响时间短;而卷破波的各项特征参数与崩破波形成对比。     

沙质岸坡上波浪传播变形及动力响应特性研究

波浪作用对水库岸坡稳定性有重要影响。为了解波浪在岸坡地形中的传播演变机制和孔隙水压力响应特性,在波浪水槽末端铺设长6 m、坡度1∶16的斜坡沙床进行试验。通过改变入射波浪参数,测量斜坡段各处波面形态,采集斜坡段不同位置处孔隙水压力,分析了波浪在沙质岸坡上浅水变形区域内波面变化特征、波能演变规律以及岸坡土体孔压特征。结果表明：随着入射波浪厄塞尔数的增大,波浪浅水变形更加明显,波形不对称性加剧,各阶谐波之间互相作用更加强烈;水深较大区域,岸坡渗透作用大于浅水变形作用,波高呈现减小趋势;浅水变形剧烈区,浅水变形作用大于岸坡渗透作用,波高呈现增大趋势,最终破碎;孔压随入射波高与波周期的增大而增大,岸坡不同位置处孔压沿深度衰减速率和随波高增长速率均不同;岸坡孔压沿深度衰减速率与入射波周期呈现出正相关关系,与波高并无太大关系。     

变性绕极深层水入侵时南极普里兹湾湍流混合特征研究

基于中国第28、29和31次南极科学考察中的CTD数据,利用Thorpe尺度方法计算了普里兹湾及其附近海域湍动能耗散率,分析了其分布特征,并对当地的水团结构进行研究。结果表明,普里兹湾及其附近海域中,前两个航次观测中次表层湍动能耗散率强度在陆架坡折区域达到最大。在水团分布方面,在第28和29航次中均观测到了变性绕极深层水陆架入侵现象,水团分别向上涌升至海表以下100 m和200 m深度,向南均可达到67.5°S处。普里兹湾陆架坡折区域次表层湍动能耗散率强度分布与当地水团结构存在良好对应关系。研究认为变性绕极深层水入侵陆架,会使该深度水体变得不稳定,发生水体交换现象,最终造成该区域湍流混合强度加强。     

Turbulence and energy dissipations of surf-zone spilling breakers

A laboratory study on the turbulence and wave energy dissipations of spilling breakers in a surf zone is presented. Instantaneous velocity fields of propagating breaking waves on a 1/20 slope were measured using Particle Image Velocimetry (PIV). Due to the large region of the evolving wave breaking generated turbulent flow, seven PIV fields of view (FOVs) were mosaicked to form a continuous flow field in the surf zone. Mean and turbulence quantities were extracted by ensemble averaging 25 repeated instantaneous measurements at each FOV. New results for distribution and evolution of turbulent kinetic energy, mean flow energy, and total energy across the surf zone were obtained from analyzing the data. The turbulence dissipation rate was estimated based on several different approaches. It was found that the vertical distribution of the turbulence dissipation rate decays exponentially from the crest level to the bottom. The resulting energy budget and energy flux were also calculated. The calculated total energy dissipation rate was compared to that based on a bore approximation. It was found that the ratio of turbulence dissipation rate to total energy dissipation rate was about 0.01 in the outer surf zone and increased to about 0.1 after the breaking waves transformed into developed turbulent bores in the inner surf zone.     

Effect of Stokes drift on upper ocean mixing

Stokes drift is the main source of vertical vorticity in the ocean mixed layer. In the ways of Coriolis - Stokes forcing and Langmuir circulations, Stokes drift can substantially affect the whole mixed layer. A modified Mellor-Yamada 2. 5 level turbulence closure model is used to parameterize its effect on upper ocean mixing conventionally. Results show that comparing surface heating with wave breaking, Stokes drift plays the most important role in the entire ocean mixed layer, especially in the subsurface layer. As expected, Stokes drift elevates both the dissipation rate and the turbulence energy in the upper ocean mixing. Also, ilffluence of the surface heating, wave breaking and wind speed on Stokes drift is investigated respectively. Research shows that it is significant and important to assessing the Stokes drift into ocean mixed layer studying. The laboratory observations are supporting numerical experiments quantitatively.     

Modeling of the eddy viscosity by breaking waves

Breaking wave induced nearsurface turbulence has important consequences for many physical and biochemical processes including water column and nutrients mixing,heat and gases exchange across air-sea interface.The energy loss from wave breaking and the bubble plume penetration depth are estimated.As a consequence,the vertical distribution of the turbulent kinetic energy(TKE),the TKE dissipation rate and the eddy viscosity induced by wave breaking are also provided.It is indicated that model results are found to be consistent with the observational evidence that most TKE generated by wave breaking is lost within a depth of a few meters near the sea surface.High turbulence level with intensities of eddy viscosity induced by breaking is nearly four orders larger than υwl(=κuwz),the value predicted for the wall layer scaling close to the surface,where uw is the friction velocity in water,κ with 0.4 is the von Kármán constant,and z is the water depth,and the strength of the eddy viscosity depends both on wind speed and sea state,and decays rapidly through the depth.This leads to the conclusion that the breaking wave induced vertical mixing is mainly limited to the near surface layer,well above the classical values expected from the similarity theory.Deeper down,however,the effects of wave breaking on the vertical mixing become less important.     

Effect of surface wave breaking on the surface boundary layer of temperature in the Yellow Sea in summer

Near-surface enhancement of turbulent mixing and vertical mixing coefficient for temperature owing to the effect of surface wave breaking is investigated using a two-dimensional (2-D) ocean circulation model with a tidal boundary condition in an idealized shelf sea. On the basis of the 2-D simulation, the effect of surface wave breaking on surface boundary layer deepening in the Yellow Sea in summer is studied utilizing a 3-D ocean circulation model. A well-mixed temperature surface layer in the Yellow Sea can be successfully reconstructed when the effect of surface wave breaking is considered. The diagnostic analysis of the turbulent kinetic energy equation shows that turbulent mixing is enhanced greatly in the Yellow Sea in summer by surface wave breaking. In addition, the diagnostic analysis of momentum budget and temperature budget also show that surface wave breaking has an evident contribution to the turbulent mixing in the surface boundary layer. We therefore conclude that surface wave breaking is an important factor in determining the depth of the surface boundary layer of temperature in the Yellow Sea in summer.     

基于精细温度观测的南海东北部陆坡-深海盆底层湍流混合*

南海是存在强湍流混合的边缘海之一, 但前人对南海湍流混合的研究更多关注的是中上层, 对底层则鲜有关注。本文基于高分辨率温度传感器于2019年5月在南海东北部22个站位海底上方0.5m处持续观测4.4d的温度数据, 分析了2216~3200m深度范围内底层海水温度的时间变化特征, 并探讨了地形粗糙度和内潮对底层湍流混合的影响。分析结果表明, 南海东北部各站位底层海水的温度变化量级约为10^-4~10^-3℃; 温度变化趋势与正压潮变化趋势不同, 温度能谱显示多数站位在全日和半日频带区间出现谱峰, 温度变化更多地受斜压潮影响, 全日、半日内潮起主要调制作用。陆坡-深海盆过渡区及深海盆底层的湍动能耗散率量级为10^-10~10^-9m²∙s^-3, 涡扩散系数量级为10^-4~10^-3m²∙s^-1。观测数据未能显示底层湍流混合与地形粗糙度存在明显的相关性。底层湍流混合的空间分布与过去观测到的南海北部深海盆内潮的南北不对称性分布一致。     

大洋细结构的一种可能的机制Ⅲ.对混合过程动能耗散的估计

继第部分之后研究了惯性内波和近惯性内波由f～的作用所致的剪切不稳定引起的破碎机制。物理上,该机制很象存在由风应力所致薄表面涡旋漂流层时表面波的破碎与饱和过程。惯性内波和近惯性内波的破碎产物与小尺度湍流一起形成了混合块,它与Gregg等人(1986)的持久混合观测结果一致。依据Thorpe(1973)实验的结果作者提出了一个估计湍流动能耗散率和消衰时间的方法。结果表明,在剪切不稳定中近惯性内波在湍动耗散中起了关键作用,而惯性内波引起非常弱的湍动耗散。使用内波能量谱的标准总能量密度估计出的近惯性内波的耗散率和消衰时间与PATCHEX测量结果非常一致。文中还讨论了几个与此破碎机制有关的问题。     

Simulation of the ocean surface mixed layer under the wave breaking

A one-dimensional mixed-layer model, including a Mellor-Yamada level 2.5 turbulence closure scheme, was implemented to investigate the dynamical and thermal structures of the ocean surface mixed layer in the northern South China Sea. The turbulent kinetic energy released through wave breaking was incorporated into the model as a source of energy at the ocean surface, and the influence of the breaking waves on the mixed layer was studied. The numerical simulations show that the simulated SST is overestimated in summer without the breaking waves. However, the cooler SST is simulated when the effect of the breaking waves is considered, the corresponding discrepancy with the observed data decreases up to 20% and the MLD calculated averagely deepens 3.8 m. Owing to the wave-enhanced turbulence mixing in the summertime, the stratification at the bottom of the mixed layer was modified and the temperature gradient spread throughout the whole thermocline compared with the concentrated distribution without wave breaking.