田 静 张 凡 李任之 高 超 吴 昊邢 飞 高建华 汪亚平,①

(1.南京大学地理与海洋科学学院 江苏南京 210023; 2.华东师范大学河口海岸学国家重点实验室 上海 200241)

湍流通常被认为是一种高能旋转的涡流运动状态, 其引起的混合速率远高于单独的分子扩散过程(Thorpe, 2007)。在河口系统中, 湍流混合是控制物质(如沉积物、污染物质、营养盐等)及能量(如动量、热量等)交换的重要物理过程, 对河口沉积、地貌演化、生物化学过程等有着重要的影响(Bakhoday-Paskyabiet al, 2018)。湍流混合相关的特征参数包括摩阻流速、拖曳系数、湍动能、湍动能剪切生成、湍动能耗散率、浮力通量等。其中, 摩阻流速和拖曳系数是描述底边界层流体特性的重要参数(Soulsbyet al, 1981), 其准确性是影响河口水动力数值模型精度的主要因素(Davieset al, 1994; Guoet al, 1998)。而湍动能剪切生成、湍动能耗散率、浮力通量等则是湍流能量平衡的主控因子(Trowbridgeet al, 1999; Rippethet al, 2001)。近年来关于湍流能量平衡的研究显示, 河口海岸浅水区中的湍流能量大多处于非局地平衡状态, 即在忽略湍动能时空变化的前提下, 湍动能剪切生成无法与湍动能耗散率、浮力通量相平衡。Talke 等(2013)在斯诺霍米什河口的观测中发现, 近底部摩擦是湍动能的源, 该处产生的湍动能向上运输, 导致近底部湍动能剪切生成远超过湍动能耗散率与浮力通量之和。随着距底距离的增加, 湍动能的垂向扩散使得湍动能剪切生成小于湍动能耗散率与浮力通量之和。Li等(2022)利用安装在底层的声学多普勒点式流速仪(acoustic doppler velocimeters, ADV)对长江口进行湍流观测, 发现底层水体在涨潮期间出现了湍流局地平衡, 落潮期间则呈现出湍动能耗散率远大于湍动能剪切生成的现象, 并把这种现象归因于湍动能对流项的影响。林姚坤等(2018)在长江口的观测显示,受水体层化的影响, 水体的能量关系仅在水体底部混合层内符合湍动能局地平衡。Wang 等(2015)通过对梅里马克河口的多年观测发现, 羽流虽然会造成河口水体的分层, 但也会产生强剪切并形成强的湍流, 且强剪切效应对湍流发育的影响更加明显。因此,河口海岸地区的湍流能量受垂向扩散、水平对流、水体层化等多种因素的影响, 从而呈现出复杂的特征。湍流结构和能量转换过程的研究对认识河口物质混合、构建海洋模型有重要意义。

河口湍流研究主要依赖于高分辨率高频现场观测。声学多普勒点式流速仪ADV 由于具有高频的突出优点, 被广泛应用于河口湍流观测中(Voulgariset al, 1998)。但由于其只能进行单点观测, 因此很难获取完整水柱的高质量湍流剖面资料。传统声学多普勒海流剖面仪ADCP 虽然可以获取高频流速剖面数据, 但由于受采样频率限制, 其声学信号噪音较大,致使湍流信号被噪音覆盖。此外, ADCP 的垂向流速是通过多个波束合成得到的, 其水平流速对湍流信号的影响难以被完全摒除(Staceyet al, 1999)。同时,在水深较浅的河口海岸区域, 微结构剖面仪等主流湍流剖面观测仪器的使用受限, 无法进行高密度的长期观测(刘志宇, 2009)11。Nortek 公司新近推出的新型5 波束Signature 1000kHz AD2CP 具有宽带、高频、低噪等特点, 且其第5 波束可直接观测垂向流速, 由此可以直接计算雷诺应力、湍动能等湍流特征参数,显著优化了四波束方差法(Guerraet al, 2017)。

长江是世界第三大河, 其河口地区的水动力结构受到河流冲淡水、风力驱动的陆架环流、潮流、波浪以及人类活动等多种过程的共同控制(沈焕庭等,1999; 时钟, 2001), 湍流效应显著(熊龙兵, 2014)。受限于观测技术、仪器精度以及水深等多重因素, 长江口近海区域的湍流剖面资料十分匮乏, 人们对长江口水体的湍流结构知之甚少。为此, 本文尝试在长江口近海海域投放座底式三角架, 采用AD2CP 和ADV 进行同步观测, 获取长江口湍流剖面结构特征, 为深入开展长江口物质混合研究和模型开发等提供关键资料。

1 材料与方法

1.1 现场观测

为减少径流的影响, 本文选择枯季在长江口不同水深区域的两个站位进行现场观测(图1a 与图1b)。T1 站位(122°00.690′E, 31°06.553′N)位于长江口南槽拦门沙海域, 平均水深8.2 m, 观测时间段为2019 年1 月12 日18:00 至 1 月 16 日 22:00; T2 站位(122°57.110′E,30°56.379′N)位于长江水下三角洲前缘, 平均水深44.7 m,观测时间为2019年12月10日14:00至12月16日20:00。

图1 研究区域(a)、观测站位(b)、现场仪器安装示意图(c 和d)以及AD2CP 俯视图(e)Fig.1 The study area (a), the location of observation stations (b), schematic diagram of field instrument installation (c, d) and top view of the AD2CP (e)

观测期间, ADV 和光学后向散射浊度计(optical backscatter point sensor, OBS)固定于座底式三脚架上,ADV 采样方式为Burst 模式, 每10 min 以16 Hz 的采样频率测量高频流速512 s; OBS 采样频率为1 Hz, 采样周期为5 min, 连续测20 s (图1c)。新型5 波束声学多普勒流速剖面仪(Nortek Signature 1 000 kHz AD2CP)被固定于另外一个相近的座底式四脚架上进行仰视观测(图1d), 采样方式为 Burst 模式, 每20 min 以8 Hz 的采样频率测量高频流速512 s (表1)。同时, 在距离海底三脚架约200 m 处布置WatchKeeper Buoy 型波浪浮标(加拿大AXYS 公司生产)以直接获得每小时的波浪数据。此外, 将一台温盐深仪Seabird19+CTD 和一台OBS 安装于船载观测架中,利用绞车每1～2 h 进行垂向剖面观测, 以获取海水的温度、盐度、浊度剖面数据, 同时每小时进行底、中、表三层现场水样采集。

表1 观测仪器及参数设置Tab.1 Instruments and their set up in the field observations

1.2 流速数据处理

AD2CP 与ADV 测量时易受到环境因素的干扰(如大颗粒悬浮物、水体气泡、边界回声等因素), 现场观测数据中会存在不同程度的噪声信息, 因此需要通过相空间阈值法对原始Beam坐标系下的高频流速数据进行预处理, 以获取精确湍流脉动数据(鲁远征等, 2012)。

随后, 将高频流速数据转换变为地球坐标系(eastnorth-up, ENU)。然后进行经验正交分解以获得主流向、次流向的高频流速数据, 包含平均流速、波浪脉动流速以及湍流脉动流速, 即

其中,u、v、w分别为主流向、次流向以及垂向瞬时流速(单位: m/s);u、v、w为潮流引起的平均流速(单位: m/s);u～、v～、w～为波浪运动引起的波浪脉动流速(单位: m/s);u′、v′、w′为湍 流 运 动引起的湍流脉动流速(单位: m/s)。AD2CP 与ADV 获取的数据分别采用改进的成对差分法(Feddersenet al, 2007)与相位法(Brickeret al, 2007)对波浪脉动流速与湍流脉动流速进行分离(图2)。

图2 剔除波浪前后的湍流频率能谱密度 φ( f)举例Fig.2 Examples of wave-turbulence spectral decomposition

由于AD2CP 固定在海底, 其倾角值较小, 可采用Dewey 等(2007)的优化四波束方差法计算雷诺应力:

其中,θ为AD2CP 波束的倾角, 为25°;φ2=φˆ3, 代表AD2CP 的横摇;φ3=-φˆ2, 代表AD2CP 的纵摇;ui′2为沿波束方向湍流脉动流速的方差,i= 1, 2,…,5分别代表AD2CP 的第1、3、4、2、5 波束。在海洋研究中, 一般只考虑水平动量的垂直通量(刘志宇, 2009)5,因此下文中提到的雷诺应力为式(7)、式(8)之和。

1.3 悬沙数据处理

采用现场水样的悬沙浓度来标定OBS 的浊度数据Tur(单位: NTU), 结果见图3a 和图3b。

图3 悬沙浓度与各站位仪器标定曲线Fig.3 Suspended sediment concentration calibration curves

ADV 和AD2CP 除了能够测量高频流速, 其记录的回声强度还能反映悬浮沉积物后向散射强度。ADV的回声强度与悬浮沉积物浓度之间的关系如下(Voulgariset al, 2004)1664:

其中,IADV为ADV 测得的回声强度;C为悬沙浓度(单位: mg/L);a、b为回声强度与悬沙浓度之间通过线性回归拟合得到的斜率和截距。基于AD2CP 的声学剖面信号反演悬沙浓度的方法与ADV 类似, 但需要对沿声波波束路径上的信号衰减进行修正(汪亚平等,1999)。声波在水体中的衰减机制主要包括声波传输损失、散射损失和水体吸收损失(Urick, 1983; Thorneet al, 1991), 因此引入“范围归一化”回声强度(E)的概念(Lohrmann, 2001):

其中,H为仪器探头距水面距离(单位: m);H0为测量水体距水面距离(单位: m);Lθ为波束倾角, 对于AD2CP 倾斜波束(即第1～4 波束)而言,θL=25 °; 对于垂向波束(即第 5 波束)而言,θL=0 °。ADV 与AD2CP 利用回声强度反演悬沙浓度的结果见图3c～3f。

1.4 水体层化指标

河口水体温度、盐度或者悬浮沉积物不均匀混合会产生垂向密度梯度从而引发层化效应(Simpsonet al, 1990)。通常用无量纲的Richardson 数Ri 来表征水体的层化程度, 具体计算为势能与动能的比值(Turneret al, 1974):

其中,sρ为沉积物密度, 其值取2 650 kg/m3。

1.5 湍流特征参数

采用雷诺应力法确定摩阻流速, 表达式为:

假定每个采样间隔期间水体温度、盐度基本不变,水体密度变化仅由悬沙浓度变化引起(由 ADV、AD2CP 反演获得高频悬沙浓度), 可计算出浮力通量B(湍流垂向混合转化的势能, 单位: m2/s3):

其中,ρ′为由悬沙浓度变化引起的水体高频密度波动(单位: kg/m³)。

2 结果

2.1 水动力时空特征

观测期间, T1 站位(平均水深8.2 m)有效波高为0.15～0.45 m, 处于弱波浪环境(图4a)。受局部地形影响, T1 站位表现为往复流的特征, 涨、落潮平均流向分别为305°、125°, 平均历时分别约为5.6、6.9 h, 最大流速分别为1.74、1.26 m/s (图4b 和图4c)。在涨、落潮转流时刻, 底层水体的潮流转向比上层水体提前约1～2 h。水体温度在7.9～8.7 °C 之间, 由于冬季海面气温低, 水体呈现一定程度的逆温现象(图4d)。受径流影响, 水体盐度为2.2～17.8 (图4e), 最高值出现于落急时刻的中下层水体, 最低值出现于落憩、初涨时刻的表层水体。涨潮时, 高盐度水体由底层入侵,导致水层盐度升高的时刻由底及表产生了显著相位滞后。水体密度受温度、盐度以及悬沙浓度的影响,处于1 001.5～1 014.0 kg/m³范围内(图4f)。水体悬沙浓度为41～781 mg/L, 其中近底层水体悬沙浓度较高(图4g)。

图4 T1 (a～j)和T2 (k～t)站位的有效波高、主流向流速、次流向流速、温度、盐度、密度、悬沙浓度、浮力频率N2、流速剪切S2 以及理查森数Ri 剖面图Fig.4 The vertical profiles of significant wave height, velocity along the major-axis, velocity along the minor-axis, temperature,salinity, water density, suspended sediment concentration, buoyancy frequency N2, the square vertical shear S2, and Richardson number Ri at T1 (a～j) station and T2 (k～t) station, respectively

观测期间, T2 站位(平均水深44.7 m)有效波高为0.29～2.34 m, 11 日15:00 至12 日16:00 间波浪显著增强, 最大有效波高达2.34 m (图4k)。潮流为旋转流,涨、落潮平均流向分别为344°、164°, 平均历时分别约为6.0、6.5 h。T2 站位距入海口较远, 受径流影响作用较弱, 涨、落潮最大流速分别为0.92、0.93 m/s(图4l 和图4m。由于受AD2CP 仪器自身观测范围的限制, 有效数据范围仅能覆盖T2 站位中下层水体),温度为14.2～18.2 °C (图4n), 盐度为30.6～33.2 (图4o),密度为1 022.4～1 024 kg/m³ (图4p), 悬沙浓度相对较低, 为37～277 mg/L (图4q)。

浮力频率N为度量水体垂向自由振动频率的物理量, 用于表征层化水体的稳定性。T1 站位N2处于1.3×10-7～4.5×10-2s-2范围, 其低值大多位于底层水体, 而高值多出现于中上层水体, 水体层化、稳定性较强(图4h)。水体流速剪切平方S2呈现较为显著潮周期变化, 且由底及表呈现相位滞后现象(图4i)。由于受到海底摩擦效应的影响, 底层水体流速剪切平方S2一直处于较高水平。根据层状剪切流动的线性理论,当Ri>0.25 时, 水体流速剪切作用强, 促使湍流发育;当Ri<0.25 时, 水体层化结构稳定性较强, 抑制湍流发育(Zhanget al, 2018)。如图4j 所示, 水体除了底层混合相对均匀外, 基本处于强层化状态(图中黑线为Ri=0.25 的分界线, 下同)。T2 站位N2为2.2×10-6～2.4×10-3s-2(图4r),S2为1.2×10-4～6.1×10-3s-2(图4s);距底10 m 以下水体混合均匀, 以上水体层化较强(图4t)。

摩阻流速u*可反映近底边界层水体切应力的大小, 是影响湍流混合的主要因素。T1 站位摩阻流速在潮周期内出现两次峰、谷值, ADV 与AD2CP 所测摩阻流速平均值分别为1.3×10-2m/s 和8.3×10-3m/s, 其变化趋势与底部流速相位基本一致(图5a 和图5b),与之相关的拖曳系数Cd平均值分别为2.9×10-3与1.4×10-3(图5c)。T2 站位ADV1、ADV2 与AD2CP所测摩阻流速平均值分别为 1.6×10-2、1.2×10-2、8.0×10-3m/s (图5e), 与之相应的拖曳系数平均值为3.2×10-3、1.9×10-3和5.0×10-4(图5f)。数据分析结果显示, ADV 和AD2CP 计算的结果较为吻合, 且T1 站位拖曳系数与前人在南槽的观测结果数量级基本一致(汪亚平等, 2006; 兰庭飞等, 2019)。

2.2 湍流剖面特征

T1 站位由ADV 与AD2CP 计算的近底部雷诺应力范围分别为1.3×10-6～1.0×10-3m2/s2和7.5×10-7～1.0×10-3m2/s2(图6a)。总体而言, 雷诺应力的大小主要受水平流速的控制, 两种仪器观测得到的雷诺应力差异并不大。雷诺应力剖面时间序列呈现出一定的周期变化, 在 1.5×10-7～3.8×10-3m2/s2范围内变动,大多数时刻处于10-5～10-3量级之间(图6b), 其高值主要出现于海底与海表附近, 低值则主要出现于中部水体。这可能是由于受到海底摩擦效应的影响, 底部湍流能够充分发育, 增强了水平动量的垂向通量;而退潮时河流上层的羽流与水体剪切作用强烈, 加之海表附近可能受波浪能量输入的影响, 雷诺应力相应增大, 故而出现了底、表两个雷诺应力的高值区。

图6 T1 (a～e)和T2 (f～j)站位近底部雷诺应力SR 时间序列和所观测水体剖面的雷诺应力SR、湍动能剪切生成P、湍动能耗散率ε、浮力通量B 剖面图Fig.6 The time series of Reynolds stress SR near the bottom and vertical profiles of Reynolds stress SR, turbulent kinetic energy production P, turbulent kinetic energy dissipation ε, and buoyancy flux B at T1 (a～e) station and T2 (f～j) station, respectively

T2 站位雷诺应力范围为3.2×10-7～2.0 ×10-3m2/s2,主要集中在10-5～10-3量级间, 高值区域主要分布在水体底层, 部分时刻上层也出现高值, 而低值区域主要集中于中层水体(图6f)。同时, 在强波浪时刻, 整层水体雷诺应力均有增大, 且相对于强波浪存在约4 h 的相位滞后。

T1、T2 站位的湍动能剪切生成P集中在10-8～10-3m2/s3(图6c 和图6h), 高值区主要分布在底、表层水体。湍流剪切生成受雷诺应力与水平流速的垂向梯度共同影响, 湍流剪切生成与雷诺应力的相位不尽相同。T1 站位湍动能耗散率ε集中在 10-7～10-3m2/s3, 同样是底、表层水体较高(图6d); T2 站位湍动能耗散率现象与T1 相似(图6i)。T1 站位浮力通量B为10-10～10-5m2/s3(图6e), 其高值集中于底层水体, 低值集中于中层水体; T2 站位浮力通量B为10-11～10-6m2/s3(图6j), 分布相对较均匀, 在强波浪时刻整个水体出现浮力通量较高的现象。整体而言,两站位由悬沙浓度变化引起的浮力通量较同层位的湍动能剪切生成、湍动能耗散率小1～2 个数量级, 可忽略不计。

2.3 剪切和层化对湍流混合的影响

为了探究剪切和层化对湍流结构变化的影响,本文将多个潮周期相位平均为一个完整的潮周期(Voulgariset al, 2004)1672。因T1 站位观测数据覆盖了整个水层, 故对其进行重点分析(图7a～7e)。T1 站位底层水体为强剪切、弱分层, 湍流混合作用强烈。在相对时间-12.5～ -9.5 h 与0～3 h 的落憩至初涨时间段内, 整个水层具有较高的剪切不稳定性, Ri 虽大于经典临界值0.25, 但小于1, 水体混合相对较为均匀。存在强剪切的分层水体可发生剪切不稳定。在这一过程中, 湍流持续混合, 平均流的动能转化为湍动能耗散与浮力通量。在相对时间-3～ -1 h、8.5～11.5 h的涨憩至落急时间段内, 流速垂向剪切较弱, 而水体层化较强, Ri 在整个水体内均处于高值, 这抑制了湍流的发育, 限制了湍动能的产生以及湍流混合,因而出现了低湍动能剪切生成和低湍动能耗散率的现象。熊龙兵等(2014)也观测到类似的现象。T2 站位近底层水体混合均匀, 湍流充分发育; 中上层水体Ri 量级大多处于100～102范围内, 水体较弱的剪切以及较强的层化状态一定程度上抑制了湍流的发育(图7f～7j)。

图7 相位平均后T1 (a～e)和T2 (f～j)站位的浮力频率N2、流速剪切S2、理查森数Ri、湍动能剪切生成P 以及湍动能耗散率ε 剖面图Fig.7 The phase-averaged vertical profiles of buoyancy frequency N2, square vertical shear S2, Richardson number Ri, turbulent kinetic energy production P, and turbulent kinetic energy dissipation ε at T1 (a～e) and T2 (f～j) stations, respectively

3 讨论

其中, 左式代表ETK时间变化率;T为能量通量的散度,包括脉动压力做功项、对流项以及黏性输运项。假设湍动能方程中的T、B可忽略不计, 则存在着P与ε间的局地平衡。

将多个潮周期的湍流垂向结构数据在涨、落潮阶段分别进行平均化处理, 用以进行局地平衡特征分析, 并计算参数来表征湍动能剪切生成P与湍动能耗散率ε的差异情况(图8)。结果表明, 随着距底高度增加, T1 站位的湍动能剪切生成与湍动能耗散率均呈现出先减小后增大的趋势(图8a, 8b), 且在底层与上层水体处二者差异最大(图8c)。涨潮阶段, 距底0.3 m 处的ADV 测得湍动能剪切生成明显大于湍动能耗散率(分别为3.4×10-4m2/s3与8.2×10-5m2/s3)。AD2CP 在距底3 m 处测得的湍动能剪切生成与湍动能耗散率差异不大, 处于同一量级(分别为 5.8×10-6m2/s3和5.4×10-6m2/s3), 达到局地平衡。而距底3 m 以上水体的湍动能耗散率均大于湍动能剪切生成。落潮阶段的垂向结构与涨潮阶段较为相似, 距底0.3 m 处的湍动能剪切生成大于湍动能耗散率, 距底2 m 处达到局地平衡(图8a)。但落潮阶段较涨潮阶段湍动能剪切生成与湍动能耗散率分别小了49.0%与47.6% (图8b)。湍动能剪切生成与湍动能耗散率在涨、落潮周期内的不对称分布可能为潮汐应变所致。潮汐应变使得涨潮时水体层化减弱而落潮时水体层化加强, 从而导致涨潮时湍流作用强于落潮(熊龙兵等,2014)。此外, T1站位落潮流弱于涨潮流, 也可能是湍动能剪切生成与湍动能耗散率涨、落潮呈不对称分布原因之一。

图8 T1 和T2 站位湍流垂向结构Fig.8VerticalstructureofturbulenceatT1andT2stations

在近岸水体中, 湍动能的源主要有两部分, 一部分来自于海底底部摩擦, 另一部分则是海表波浪的能量输入(Grassoet al, 2012; Talkeet al, 2013)。T1 站位底摩擦产生的湍动能向上运输, 导致近底部的湍动能剪切生成大于湍动能耗散率与浮力通量之和。由于湍动能在水体中的垂向扩散, 湍动能剪切生成与湍动能耗散率均随着距底高度的增加而减小, 涨潮时距底3 m 处、落潮时距底2 m 处达到局地平衡, 上层水体则表现出湍动能耗散率大于湍动能剪切生成的现象。另外, 由于落潮时河流羽流与水体的剪切作用增强, 并且海表处常常伴随波浪的能量输入, 中上层水体的湍动能剪切生成与湍动能耗散率又逐渐增大。在本文中, 浮力通量是基于温盐恒定、密度变化仅由悬沙浓度变化引起的假定条件计算得到, 比湍动能剪切生成与湍动能耗散率小1～2 数量级, 可忽略不计。但刘兴泉等(2004)、王鹏皓(2020)表明, 长江口区域内温盐结构可能发生活跃变化, 因此, 在水体层化结构稳定的T1 站位, 无法排除浮力通量对湍动能非局地平衡的影响。综上所述, T1 站位的湍动能非局地平衡可能为温盐等斜压作用引起的浮力通量与对流共同作用所致。

T2 站位的湍流垂向结构与T1 站位相似, 湍动能剪切生成与湍动能耗散率, 二者差异在上层水体处最大(图8f), 落潮阶段的湍动能剪切生成与湍动能耗散率比涨潮阶段分别小了39.1%与24.7% (图8d和图8e)。T2 站位水体基本处于混合均匀或弱层化状态(浮力通量可忽略不计), 观测期间有强波浪事件(有效波高大于2 m), 在强波浪作用下, 波浪破碎通过影响脉动压力做功与黏性输运, 对湍流场造成显著影响(Scullyet al, 2016)。因此T2 站位可能是由对流、脉动压力做功与黏性输运共同作用引起的。

4 结论

本文在长江口及其邻近海域不同水深的区域(8.2 m, 44.7 m)采用新型5 波束Signature 1 000 kHz AD2CP 和ADV 为主要仪器进行座底式观测, 得出以下主要结论:

(1) 基于ADV 和AD2CP 获得的近底部边界层参数(摩阻流速、拖曳系数以及雷诺应力)结果基本一致,两种仪器的组合观测能在一定程度上提高湍流剖面观测的时空解析度。

(2) 长江口南槽受径流作用影响显著, 在水体混合较为均匀的落憩至初涨时刻, Ri 的量级范围约为10-6～100, 强剪切作用促进了湍流混合; 在水体层化结构稳定性较强的涨憩至落急时刻, Ri 约为100～104,湍流混合被抑制。水下三角洲前缘海域由于距入海口较远, 其水体混合程度强于近岸海域, Ri 多处于10-6～102。

(3) 研究区湍动能主要有底层摩擦与海面波浪两种能量输入源。底摩擦作用致使底层出现P>ε的湍流非局地平衡现象, 过剩的湍动能通过拟序结构向上迁移, 致使底层上方水体的湍流特征参数逐渐减弱且P的衰减程度强于ε。由于波浪能量的输入以及落潮时河口上层羽流所产生的强剪切, 海表处P、ε增大且出现P<ε的现象, 水体呈现向下消耗的湍流耗散结构。受近岸径流作用的影响, 长江口南槽湍流非局地平衡现象的主因可能为湍动能方程中温盐变化引起的浮力通量以及对流效应。而长江水下三角洲前缘海域湍流非局地平衡现象可能是对流和强波浪作用影响的脉动压力做功、黏性输运因素导致的。

由于现场观测的时空局限性, 在未来研究中, 将进一步开展空间尺度上湍流能量通量的散度研究、完善不同季节长江口海域的观测, 这对于充分了解湍流混合、认知河口物质输运、构建高精度模型具有重要意义。

致谢 南京大学唐杰平、陈德志、卢婷、兰庭飞与华东师范大学盛辉、汤碧璇、常洋参加了现场观测及样品处理分析工作, 同济大学涂俊彪、孟令鹏在数据处理方面给予诸多帮助, 谨致谢忱。