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

利用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。湍动能耗散率和垂向湍扩散系数的分布相似,均为上层最强,底层次之,中层最弱。上层由于风应力的作用,使得湍动能耗散率和垂向湍扩散系数较大;温跃层处层化较强,抑制了湍动能的耗散和垂向上的湍混合。盐度锋面的次级环流会促使低盐水团脱离,锋面引起的垂向环流会加强海洋的湍混合。低盐水团与外界的能量交换较少,湍动能耗散率较弱。长江口海区存在明显的上升流和下降流,它们是由锋面的次级环流产生的;上升流和下降流的存在促进湍动能的耗散与湍混合。     

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

利用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)的拟合函数关系,反映了层化对水体垂向混合的抑制作用。     

北冰洋次表层暖水形成机制的研究

分析了北冰洋加拿大海盆上层海水温度结构特征,针对在海冰覆盖区域普遍发生的次表层暖水现象,建立了冰海耦合的一维柱形模式,成功模拟了次表层暖水的垂向结构和形成机制,证明了太阳辐射加热和表层冷却是次表层暖水的成因。通过数值实验定量分析了形成次表层暖水各因子的作用,表明太阳辐射是形成次表层暖水的能量来源;长波辐射、气温和比湿等其他大气因子对次表层暖水强度有较大影响;在厚冰区不能形成次表层暖水,海冰越薄,水道面积越大,则冰下次表层海水温度升高得就越快,薄冰和冰间水道是形成次表层暖水的主要能量通道;通过模式的拟合给出湍扩散系数垂向变化的估值,冰下表层海水垂向湍扩散系数为5.0×10-3m2/s左右,而次表层垂向湍扩散系数却快速降低到1.0×10-6m2/s,表明跃层导致的稳定性增大是形成次表层暖水尖峰的关键因素。     

西北太平洋137°E断面温度场和盐度场的时空特征

本文利用日本气象厅在137°E断面获得的水温和盐度长期观测资料,分析了该断面温度场和盐度场的时空特征.结果表明,137°E断面的温度场和盐度场都存在着明显的季节差异和年际变化.冬季,温度场变化的关键区位于3°～18°N的300m以浅海域,而盐度场变化的关键区则位于18°～34°N的300m以浅海域.夏季,温度场变化的关键区位于3°～16°N的300m以浅海域,而盐度场则有两个关键区,分别位于3°～18°N的200m以浅海域和24°～34°N的300m以浅海域.温度场的年际变化与ENSO循环相联系,而盐度场的年际变化则比较复杂.     

南极普里兹湾海域湍流扩散系数估计

基于Thorpe尺度方法,利用CTD数据,计算了南极普里兹湾海域的Thorpe尺度和湍流扩散系数,分析了观测区域(64°～69°S,66°～80°E)湍流翻转现象的强弱及分布。结果表明,在海底和地形粗糙区存在较大的Thorpe尺度(较强湍流翻转)和湍流扩散系数,湍流扩散系数最大值能达到10^-2m²/s量级,比平坦开阔海洋高2～3个数量级,部分观测站位的湍流扩散系数和湍动能耗散率表现出大-小-大的垂向分布结构,总水深较深的区域尤为明显;深水区域的浮力频率在海表面到500 m层比较大,浅水区域该现象不明显;湍动能耗散率在(67.25°S,73°E)周围和经度为78°E的各站位都表现相对较大,能达到10^-6 w/kg量级,个别站位甚至能达到10^-5 w/kg量级。     

湍能封闭模型在渤海潮流模拟中的应用

将一个三维湍能封闭模型应用于渤海潮流模拟，通过流速分解的初边值方法，计算了渤海全海区潮波系统、潮流结构，结果与观测符合较好。计算同时给出了湍粘性系数和湍动能的时空分布，比较常湍粘性系数和湍的垂直抛物分布及湍能封闭模拟表明，湍粘性系数的不同选取会对流场结构有显著影响。     

夏季黄东海硝酸盐垂向扩散通量的分析

湍流扩散过程导致的硝酸盐垂向输运对海水表层的浮游植物生长和初级生产力的大小有着重要影响。本文基于2018年夏季黄、东海水文环境、硝酸盐浓度和湍动能耗散率的同步、原位数据，分析了海域温度、盐度和硝酸盐的空间分布特征，结果表明营养盐含量丰富的黄海冷水团、长江冲淡水、东海北部底层混合水与黑潮次表层水是影响研究海域硝酸盐分布的主要水团。利用垂向湍扩散硝酸盐通量公式，计算了三个选定断面上的硝酸盐垂向扩散通量，其高值区与湍流扩散系数的高值区的位置基本一致。针对存在明显硝酸盐跃层的站位，计算得到跨硝酸盐跃层的垂向通量F_ND的范围在-9.78—36.60mmol/（m²·d）之间，在黄海冷水团区，夏季温跃层限制了该区营养盐向近表层的湍流垂向扩散；东海北部底层混合水区，湍流垂向扩散向上层补充了大量硝酸盐，促进了跃层之上浮游植物的生长；黑潮次表层水影响海区，夏季中层水体混合较弱，跨跃层的垂向通量也普遍偏低。开展硝酸盐垂向扩散通量的计算与分析，对进一步明确营养盐的输运机制有着重要研究意义。     

McCreary和Fickian湍扩散形式差异的分析

用矩阵方法,求解McCreary和Fickian铅垂向湍扩散形式下,无界海洋条件时在风应力和深层浮力作用下的线性解。得出不同的湍扩散系数A情况下McCreary与Fickian湍扩散形式的流场结构。指出随湍扩散系数的增加,McCreary和Fickian湍扩散形式下u分量结构及其量值差异增大。在Fickian湍扩散形式下,能量穿越密度跃层的传递能力较强。当A=0.1cm~2s~(-3)时,在风应力和深层浮力作用下两种湍扩散形式造成流场的差异很小。提出在小湍扩散系数下,用McCreary湍扩散形式研究热带海洋环流,其结果和Fickian形式应相当一致。     

黄海西部海洋湍流的季节变化特征分析

在2006—2007年开展的"中国近海海洋综合调查与评价"项目中,作者利用自由下降湍流剖面仪MSS60在南黄海海区分别进行了夏、冬、秋季三个航次的微尺度湍流观测,并计算分析了该海区的湍动能耗散系数ε,湍扩散系数κ等。通过与温度、流速分布图对比,结果表明三个季节的湍流混合趋势大体一致。在沿岸浅水区,混合作用比较强烈。而深水区湍流混合的垂直分布明显地表现出三层结构,混合较强的上混合层和底混合层,及相对较弱的中层。风混合和潮混合是黄海湍流混合的主要形式。风的影响主要表现在海洋上层,潮流的影响则表现于底层。     

垂向混合参数化方案对模拟黄海夏季上层温度结构的影响研究

采用POMgcs(Princeton Ocean Model with generalized coordinate system)和MITgcm(MIT General Circulation Model)两个海洋数值模式,研究了M-Y2.0、基于固壁近似假定的M-Y2.5、基于波浪破碎作用的M-Y2.5和KPP 4种垂向混合参数化方案对模拟黄海夏季上层温度结构的影响。结果表明,M-Y2.0和基于固壁近似假定的M-Y2.5方案低估了黄海上层的湍动能,模拟的黄海夏季温度上混合层的效果与实测相比均偏浅,不能够很好地重构黄海夏季温度的垂直结构。而基于波浪破碎作用的M-Y2.5和KPP方案均可以增加海洋上层湍动能的输入量,模拟的黄海夏季温度上混合层的效果与实测较为一致。故推测黄海夏季的上层结构是受波浪混合和流场剪切等物理机制共同调节的,若通过合理的垂向混合参数化方案将这些物理机制的作用加以体现,将会较真实地模拟和重构出黄海夏季海温上层结构。     

Zooplankton community structure in a cyclonic and mode-water eddy in the Sargasso Sea

Mesoscale eddies are important suppliers of nutrients to the surface waters of oligotrophic gyres, but little is known about the biological response, particularly that of higher trophic levels, to these physical perturbations. During the summers of 2004 and 2005, we followed the development of a cyclonic eddy and an anti-cyclonic mode-water eddy in the Sargasso Sea. Zooplankton (>150 μm) were collected across both eddies in 9 discrete depth intervals between 0 and 700 m. Comparison of the abundance of major taxa of mesozooplankton in the upper 150 m at eddy center and outside the eddies (day and night) indicated that the cyclone and mode-water eddy supported similar mesozooplankton communities, with several taxa significantly higher in abundance inside than outside the eddies, when compared with the Bermuda Atlantic Time-series Study site as representative of mean conditions. In both eddies copepod peak abundance occurred in the 50-100 m depth interval, coincident with the chlorophyll a maximum, suggesting elevated food concentration in the eddies may have influenced zooplankton vertical distribution. The two eddies differed in the strength of diel vertical migration of zooplankton, as indicated by the ratio of night:day abundance in the epipelagic zone, which was higher at the center of the mode-water eddy for most taxa. Over the sampling interval of 1-2 months, abundance of the three most common taxa (copepods, chaetognaths, and ostracods) decreased in the cyclone and increased in the mode-water eddy. This further supports previous findings that over the sampling period the cyclone was in a decay phase, while the mode-water eddy was sustaining nutrient fluxes and high phytoplankton concentrations. A more detailed analysis of community structure in the mode-water eddy indicated the 0-700 m integrated abundance of doliolids was significantly higher inside the mode-water eddy than outside. The presence of a mesopelagic (200-700 m) layer of lepadid barnacle cyprids in this eddy highlights the potential of eddies to transport and disperse biota. We conclude that when compared with average ambient conditions (as measured at BATS), mesoscale eddies can influence zooplankton behavior and alter zooplankton community structure which can affect food-web interactions and biogeochemical cycling in the open ocean.     

Three-dimensional structure of an observed cyclonic mesoscale eddy in the Northwest Pacific and its assimilation experiment

Mesoscale eddies play an important role in modulating the ocean circulation. Many previous studies on the threedimensional structure of mesoscale eddies were mainly based on composite analysis, and there are few targeted observations for individual eddies. A cyclonic eddy surveyed during an oceanographic cruise in the Northwest Pacific Ocean is investigated in this study. The three-dimensional structure of this cyclonic eddy is revealed by observations and simulated by the four-dimensional variational data assimilation(4 DVAR) system combined with the Regional Ocean Modeling System. The observation and assimilation results together present the characteristics of the cyclonic eddy. The cold eddy has an obvious dual-core structure of temperature anomaly.One core is at 50–150 m and another is at 300–550 m, which both have the average temperature anomaly of approximately-3.5°C. The salinity anomaly core is between 250 m and 500 m, which is approximately-0.3. The horizontal velocity structure is axis-asymmetric and it is enhanced on the eastern side of the cold eddy. In the assimilation experiment, sea level anomaly, sea surface temperature, and in situ measurements are assimilated into the system, and the results of assimilation are close to the observations. Based on the high-resolution assimilation output results, the study also diagnoses the vertical velocity in the mesoscale eddy, which reaches the maximum of approximately 10 m/d. The larger vertical velocity is found to be distributed in the range of 0.5 to 1 time of the normalized radius of the eddy. The validation of the simulation result shows that the 4 DVAR method is effective to reconstruct the three-dimensional structure of mesoscale eddy and the research is an application to study the mesoscale eddy in the Northwest Pacific by combining observation and assimilation methods.     

基于2005~2019年卫星遥感观测的南大洋印度洋扇区中部涡旋特征分布研究

在南大洋印度洋扇区中部海域,除了地形控制(凯尔盖朗高台),南极绕极流和厄加勒斯回流的汇合流进一步加强了下游的斜压剪切强度,导致涡旋能量显著增强,因此,对该海域涡旋的研究有助于了解该海域的涡旋特征以及地形与涡旋的分布关系。基于2005~2019年卫星遥感数据,对该海域涡旋特征进行统计,并对涡旋产生地分布、跨锋面涡旋的移动状况进行分析,同时结合Argo剖面数据,进一步剖析涡旋内部水文分布特征。结果表明:该海域涡旋生命周期多在20 d以内(64.25%),平均半径多在30~100 km(96.13%);平均半径与平均振幅呈正相关关系(相关系数R=0.55);生命周期越大的涡旋平均传播距离也越大。2014年开始涡旋数量明显增加,主要由短寿命涡旋(<30 d)数量增加所贡献。反之, 21世纪10年代后期年平均涡动能异常呈减小趋势。涡旋产生地随着寿命增长,逐渐从亚南极锋与南极绕极流南部边界之间的锋面区域向亚南极锋以北移动。跨锋面涡旋中,暖涡向高纬,冷涡向低纬移动,大部分具有携带水团移动的能力。由涡旋内部水文特征分析结果可知,不同极性的涡旋能够实现完全不同来源水团的远距离输送,对同一来源水团,气旋涡具有抬升作用,而反气旋涡具有压沉作用。该研究工作有助于提升对南大洋涡旋特征及变动的认识,为进一步的涡旋动力研究提供支撑。     

Seasonal variation of eddy kinetic energy in the South China Sea

Mesoscale eddy activity and its modulation mechanism in the South China Sea (SCS) are investigated with newly reprocessed satellite altimetry observations and hydrographic data.The eddy kinetic energy ...     

On the variability of vertical eddy heat flux in the upper ocean

Ocean eddies produce strong vertical heat flux (VHF) in the upper ocean, exerting profound influences on the climate and ecosystem. Currently, mooring array provides a standard way to estimate the eddy-induced VHF (EVHF) based on the adiabatic potential density equation. Apart from the validity of adiabatic assumption, it remains unclear to what extent the estimated EVHF at a single location within a limited time period is representative of its climatological mean value. In this study, we analyzed the above issue by systematically evaluating the variability of EVHF simulated by a 1-km ocean model configured over the Kuroshio Extension. It is found that the EVHF at a single location exhibits pronounced variability. Even averaged over one year that is comparable to the current maintenance capacity of mooring array, the EVHF still deviates significantly from its climatological mean value. For more than 49% of locations in our computational domain (31°–40°N, 149°–166°E), the discrepancy between the one-year mean EVHF and its climatological mean value at the peaking depth is larger than the climatological mean itself. The mesoscale eddies play a dominant role in the variability of EVHF but contribute little to the climatological mean EVHF; the opposite is true for submesoscale eddies. Our findings indicate that nested mooring array allowing for isolating the effects of submesoscale eddies will be useful to obtain climatological mean EVHF.     

A Study of Eddy Viscosity Coefficient in Numerical Tidal Simulation

Based on the fluid motion equations, the physical meaning of eddy viscosity coefficient and the rationality of the Boussinesq hypothesis are discussed in this paper. The effect of the coefficient on numerical stability is analyzed briefly. A semi-enclosed rectangular sea area, with an orthogonal spur dike, is applied in a 2-D numerical model to study the effect of horizontal eddy viscosity coefficient (AH). The computed result shows that AH has little influence on the tidal level and averaged flow velocity, but has obvious influence on the intensity and the range of return flow around near the spur dike. Correspondingly, a wind-driven current pool and an annular current are applied in a 3-D numerical model respectively to study the effect of vertical eddy viscosity coefficient (Av). The computed result shows that the absolute value of Av is inversely proportional to that of horizontal velocity, and the vertical gradient value of Av determines the vertical distribution of horizontal velocity. The distrib     

不同布设间距和来流速度下方型人工鱼礁上升流效应的数值模拟

合理的人工鱼礁组合可以有效改善投放水域的流场效应,提高投放水域底层与上层水体之间的扰动。通过使用并行大涡模拟模式及被动示踪物模块,并通过调整人工鱼礁布设间距,研究了在不同背景流速条件(0.1、0.5、0.6和1.0 m/s)下,在不同的横向间距(1L、2L、3L)(L表示人工鱼礁的边长)或纵向间距(1L、2L、3L、4L、5L)情况下,方型人工鱼礁对上升流体积、营养盐的抬升和垂向涡黏系数的影响。研究结果表明,在同一布置条件下,单排布置下的三块人工鱼礁形成的上升流体积大小与来流速度成正相关,体积随来流流速增加而增大6.4%～80.5%;在同一流速条件下,上升流体积大小与纵向布置的间距成正比,与横向布置的间距成反比;在横向布置条件下,当来流速度为1.0 m/s、布设间距为1L时,上升流体积参数最佳。总体来说,上升流体积参数、示踪物浓度差和垂向涡黏性系数均显示横向布置优于纵向布置,相较于布设间距,来流速度是影响上升流体积最重要的因素。     

Eddy diffusivity and coherent mesoscale eddy analysis in the Southern Ocean

The spatial distribution of eddy diffusivity, basic characteristics of coherent mesoscale eddies and their relationship are analyzed from numerical model outputs in the Southern Ocean. Mesoscale fluctuation information is obtained by a temporal-spatial filtering method, and the eddy diffusivity is calculated using a linear regression analysis between isoneutral thickness flux and large-scale isoneutral thickness gradient. The eddy diffusivity is on the order of O (103 m2/s) with a significant spatial variation, and it is larger in the area with strong coherent mesoscale eddy activity. The mesoscale eddies are mainly located in the upper ocean layer, with the average intensity no larger than 0.2. The mean radius of the coherent mesoscale cyclonic (anticyclonic) eddy gradually decays from (121.2±10.4) km ((117.8±9.6) km) at 30°S to (43.9±5.3) km ((44.7±4.9) km) at 65°S. Their vertical penetration depths (lifespans) are deeper (longer) between the northern side of the Subpolar Antarctic Front and 48°S. The normalized eddy diffusivity and coherent mesoscale eddy activity show a significant positive correlation, indicating that coherent mesoscale eddy plays an important role in eddy diffusivity.     

Energetics of lateral eddy diffusion/advection：Part Ⅰ. Thermodynamics and energetics of vertical eddy diffusion

Two important nonlinear properties of seawater thermodynamics linked to changes of water density, cabbeling and elasticity(compressibility), are discussed. Eddy diffusion and advection lead to changes in density; as a result, gravitational potential energy of the system is changed. Therefore, cabbeling and elasticity play key roles in the energetics of lateral eddy diffusion and advection. Vertical eddy diffusion is one of the key elements in the mechanical energy balance of the global oceans. Vertical eddy diffusion can be conceptually separated into two steps: stirring and subscale diffusion. Vertical eddy stirring pushes cold/dense water upward and warm/light water downward; thus, gravitational potential energy is increased. During the second steps, water masses from different places mix through subscale diffusion, and water density is increased due to cabbeling. Using WOA01 climatology and assuming the vertical eddy diffusivity is equal to a constant value of 2×103 Pa2/s, the total amount of gravitational potential energy increase due to vertical stirring in the world oceans is estimated at 263 GW. Cabbeling associated with vertical subscale diffusion is a sink of gravitational potential energy, and the total value of energy lost is estimated at 73 GW. Therefore, the net source of gravitational potential energy due to vertical eddy diffusion for the world oceans is estimated at 189 GW.