大洋深层海洋在全球气候变化和生物地球化学循环中扮演着重要的角色。西太平洋紧邻欧亚大陆和我国, 是全球地形最为粗糙的海域之一, 粗糙地形引起的巨大潮能耗散为深海动力过程提供了能量, 产生了丰富的深层多尺度动力过程。此外, 西太平洋上层有丰富的环流系统和暖池, 通过上、下层的潜在动量能量交换会影响深海动力过程, 因此西太平洋是开展深海动力过程研究的理想区域(汪嘉宁等, 2022)。季节内振荡是全球大洋中普遍存在的一个变异特征(Stammer, 1997), 其在上层的产生机制多种多样, 如中尺度涡旋等海洋过程(Zhang et al, 2013; Wang et al, 2016)和马登-朱利安振荡(Madden-Julian Oscillation, MJO)等大气强迫过程(Matthews et al, 2007); 而在中深层海洋, 水体内部的中尺度涡(Shu et al, 2022)和地形罗斯贝波(topographic Rossby waves, TRW)被认为是引起季节内振荡的主要原因。科学家在中国南海(Shu et al, 2016; Quan et al, 2021; Zheng et al, 2021)、墨西哥湾(Oey et al, 2002; Hamilton, 2009)和加拿大海盆(Zhao et al, 2018)等地点陆续观测到了TRW的存在, 研究其变异特征和动力机制。但这些海域的水深基本上均小于3 000 m, 而在开阔大洋数千米的深层由于直接观测比较少, 其季节内振荡的动力机制和与上层海洋的联系仍不明确。

雅浦-马里亚纳海沟连接区(Yap-Mariana Junction, YMJ)位于西太平洋, 在2 500 m以深形成了3个深水通道, 其中东深水通道连接马里亚纳海沟和东马里亚纳海盆, 北深水通道连接西马里亚纳海盆, 南深水通道连接雅浦海沟和西卡罗琳海盆(图 1)。YMJ是深层经向翻转环流上、下分支携带南极上、下绕极深层水入侵西太平洋的首要区域, 对西太深层水的更新和变化具有极其重要的影响, 是西太平洋深海研究的热点区域(Kawabe et al, 2003; Siedler et al, 2004; Wang et al, 2021; Zhou et al, 2022)。

基于雅浦-马里亚纳海沟连接区北深水通道的3套潜标观测数据, Ma等(2019)发现该通道流速和等温线振幅在季节内尺度上随着深度增加而增强, 其中4 200 m处流速振幅达到45 cm/s, 1.5 ℃等温线振幅可达600 m, 该现象被形象地称为“深海风暴”。研究表明该处季节内深层强化现象是由TRW引起的, 观测流速或等温线振幅的垂向分布符合TRW的波解方程, 拟合解释方差达到90%以上, 且季节内振荡具有自上而下位相基本一致的特征。文章还借助高分辨模型数据讨论了TRW的两种能量来源, 一是上层的强涡旋通过位涡守恒激发, 二是低频的深层背景流通过正压和斜压不稳定向季节内变异提供能量。上述工作是深层季节内变异研究的重要进展, 但受制于文章篇幅和数据的长度仍有一些问题有待解决, 比如TRW导致的深层季节内变异是否存在于YMJ其他区域, 是否还有上述两个机制以外的其他过程可以激发TRW, 以及在季节内时间尺度上全水体动力关联的物理图景是什么等。

本文将利用高分辨率的再分析模式数据集, 在较长时间尺度上对YMJ东深水通道(图 1)深层流的季节内变异特征、机制和全水体的动力关联进行进一步分析, 阐明该通道流速和温度的季节内变异机制对全面认知南极水的入侵具有重要科学意义。

1 数据与方法 1.1 GLORYS再分析数据

本文主要使用全球海洋再分析产品GLORYS12V1 (以下简称GLORYS) 2017~2019年逐日和逐月的海表面高度(sea surface height, SSH)、温度、盐度和流速数据。该产品由哥白尼海洋环境监测服务中心(Copernicus Marine Environment Monitoring Service)提供, 产品来源于欧洲海洋模式NEMO (Nucleus for European Modeling of the Ocean), 其经纬度水平分辨率均为1/12°, 垂直方向分为不等间距的50层, 层厚随着深度增加而增大, 在表层10 m层厚约1~2 m, 在100 m层厚约20 m, 在1 000 m层厚约200 m, 在3 000~5 000 m层厚约350~450 m。GLORYS同化了卫星高度计、海表面温度以及温度与盐度剖面等观测数据。

Ma等(2019)评估了GLORYS模式数据对YMJ区域深层温度、流速季节内变异的模拟能力, 发现其可以很好地再现观测的深层季节内变异和TRW过程。Wang等(2021)评估和证实了GLORYS模式数据对YMJ区域深层经向翻转环流季节性路径和流量的模拟能力。以上研究表明本文使用GLORYS数据研究YMJ区域的深层季节内变异机制是可行的。

1.2 计算方法 1.2.1 季节内流速和水平动能计算方法

季节内流速的东西(u′)和南北(v′)分量是通过对原始东西(u)和南北(v)流速进行20~90 d的带通滤波得到的, 滤波使用巴特沃斯滤波器(butterworth filter)。季节内水平动能大小由公式计算得到。

1.2.2 TRW的波解方程和频散关系

根据Ma等(2019)和Shu等(2022)的介绍, TRW水平流速振幅(U)的波解方程为

TRW垂直流速振幅(W_rms)的波解方程为

其中, A₀是常数, z表示深度, λ表示垂向捕获深度的倒数。

在自然坐标系中, x轴平行于等深线, y轴表示地形梯度方向, (k, l)分别为沿(x, y)轴方向波数, 表示波数矢量。TRW频散关系可以近似为, 其中ω表示TRW的频率, H表示水深, 表示地形梯度, N表示布伦特-维萨拉频率, f表示科氏参数。对于短波和水深大于捕获深度的情况, 当时, , 因此频散关系可简化为

其中, θ表示波数矢量K与y轴的顺时针方向夹角, 具体求解方法见Meinen等(1993)。根据公式(3), ω > 0需要满足沿等深线方向的k为负, 因此波数矢量K只能指向y轴向上的第二或者第三象限, 从而θ的范围为0°~180° (Hamilton, 2009)。

TRW水平波长L可由公式(4)计算

1.2.3 涡度计算方法

相对涡度ζ的计算公式为

其中, x和y表示东西和南北方向坐标, 计算选择相距4个模式格点(约为33.3 km)的u和v进行。

位涡守恒公式为

其中, h表示水层厚度。

1.2.4 涡能方程计算方法

我们将使用涡能方程来评估深层平均流不稳定性对深层季节内能量的可能贡献。一般认为, 压力梯度项、开尔文-亥姆霍兹不稳定项量值较小, 季节内总能量的时间变化可以与正压和斜压不稳定项之和相平衡(Oey, 2008; Ma et al, 2019)。季节内总能量包括季节内水平动能和有效位能两部分, 总能量时间导数的表达式为, 正压不稳定项为, 斜压不稳定项为, 其中, ρ为密度, g为重力加速度, 上划线表示120 d的滑动平均项, 异常项(')表示相对于120 d滑动平均的异常值, ρ₀为平均密度。计算采用深层3 200~4 400 m深度范围内和水平以YMJ-E站位为中心1°×1°的范围内区域内平均值。正压和斜压不稳定项的正值表示能量从平均流向季节内尺度转化。

2 结果和讨论 2.1 深层的季节内能量分布

图 2给出了YMJ及其周边海域2017~2019年4 405 m深度上季节内水平动能平均值和标准差值的分布。结果显示YMJ东深水通道位于一个季节内动能平均值的大值区(图 2a), 且具有较大变化(图 2b), 比YMJ其他通道的量值大4~5倍, 这也说明对东通道开展研究的意义。

图 3给出了2017~2019年YMJ东深水通道YMJ-E站位(140.43°E, 10.83°N, 下同)处全水深季节内动能的时间-深度变化。季节内能量在垂向上呈现两个大值区, 分别是上层1 000 m以浅和深层2 500 m以深, 中层1 500~2 500 m处仅在2017年2~4月短暂出现大值。

2.2 TRW及其引起的深层季节内振荡

图 4为YMJ-E站位全水深季节内水平动能的保能谱图, 结果显示季节内能量的显著周期在40~90 d, 深层季节内最大能量对应的周期约为68 d。季节内能量在表层至2 500 m呈现减小趋势, 而后在2 500 m至海底呈现增强趋势, 4 500 m处季节内能量值与表层和次表层相当, 这些定量特征与图 3给出的定性结果是一致的。

东深水通道为东西走向, 受此地形约束通道内深层纬向流速远大于经向流速, 下面我们仅对其纬向流速进行分析。我们对2017~2019年纬向水平流速进行经验正交函数(empirical orthogonal function, EOF)分解, 结果显示其模态一和模态二的解释方差分别为61.2%和24.5%, 第一模态可以较好地解释水平流速的主要变化, 其特征向量代表纬向水平流速的振幅。图 5a给出了2017~2019年纬向水平流速EOF第一模态的平均振幅(黑点), 可以看出振幅在2 800~4 400 m呈现随深度增加而增加的特征, 基本符合TRW的特征。我们使用TRW水平流速的波解方程[公式(1)]对振幅的垂直分布U(z)进行拟合(图 5a黑线), 拟合后的方程解为, 拟合的解释方差为93.4%, 捕获深度为420.1 m, TRW的波解方程可以很好地刻画YMJ-E站位处的水平流速振幅的垂向分布。

我们还利用现场温度等值线的振幅变化反演出垂向速度, 反演时假设垂直对流项主导深层温度的变化, 即忽略水平平流项以及垂直湍流扩散项的作用, 温度变化的控制方程简化为, 其中, T表示现场温度, t表示时间, w表示垂向流速(Ma et al, 2019)。图 6a给出了2017~2019年现场温度的时间-深度变化, 可以发现等温线的振荡随深度的增加而增强, 比如, 在1.45 ℃等温线最大垂向位移可以达到800 m, 而较浅的1.75 ℃等温线最大位移仅200 m。然后进一步利用垂直流速的波解方程[公式(2)]检验其振幅(用垂向流速w的均方根w_rms表示)是否也符合TRW的特征。垂直流速振幅的拟合方程为w, 拟合解释方差为83.6%, 捕获深度为564.6 m (图 6b)。上述分析表明TRW是引起YMJ-E站位深层流速和温度季节内振荡和其强度随深度增加而增强的主要原因。

我们进一步利用TRW的频率关系[公式(3)], 使用N和地形参数等计算TRW频率和周期, 验证其与季节内动能保能谱得到的能量最大值周期是否一致。N的计算基于GLORYS模型2017~2019年2 800~ 4 400 m平均温盐数据, 结果约为5.6×10⁻⁴ s⁻¹, 地形坡度取自经过50 km²中位滤波平滑的ETOPO1地形数据, 其值约为1.1×10⁻²。按照Meinen等(1993)方法, 基于捕获深度和地形参数等迭代求解出波数矢量和地形梯度之间的夹角θ约为11°。将上述数值带入公式(3), TRW的周期约为62 d。而季节内动能保能谱显示该通道处季节内能量的峰值周期为68 d, 两个数值基本接近, 进一步说明了TRW是引起深层季节内能量强化的主要原因。

基于公式(4)计算了YMJ-E站位处TRW的水平波长, 捕获深度取图 5a与6b中捕获深度的平均值492 m, 水平波长约为78 km。我们进一步向后追踪了YMJ-E站位60 d周期TRW的传播路径(Ma et al, 2019), 发现它在YMJ-E站位东南75 km处终止(图 7), 与水平波长数值基本相当。这表明YMJ-E站位附近的TRW主要是局地生成的, 并非由远端传播而来, 下面我们来探寻TRW局地的能量来源。

2.3 TRW的能量来源

图 8给出了YMJ-E站位不同深度上的季节内流速矢量的时间变化图, 结合图 3水平季节内动能的时间深度变化, 我们可以看出虽然深层季节内流速在全部时间段展现出TRW的特征, 但不同时刻水平流速和动能的垂直分布均不相同。我们选取了3个代表性事件: 事件一的时间段为2019年7月16日至2019年9月8日, 此时间段的季节内流速和动能在次表层存在大值区, 后随深度增加逐渐减小, 到2 866 m后又逐渐增大; 事件二的时间段为2017年2月2日至2017年3月24日, 此时间段的季节内流速和动能在中层存在大值区, 后随着深度增加量值逐渐减小, 到3 221 m后深层又逐渐增大; 事件三的时间段为2019年9月25日至2019年11月15日, 此时间段季节内流速和动能最大值位于深层, 上层至3 597 m层的量值均较小。三个事件中季节内能量垂直分布显著差异说明其深层季节内能量来源可能存在着差异, 下面我们将对三个事件进行单独分析来进行探究。

2.3.1 事件一

我们首先使用TRW波解方程分别对事件一时间段的水平(图 5b)和垂直流速(图 6c)振幅进行拟合。水平流速振幅的拟合方程解为U(z)=1.6×10^-5cosh(0.0023z), 解释方差为95.4%, 捕获深度为443.4 m。垂直流速振幅的拟合方程解为w_rms(z)=1.1×10^-6sinh(0.0017z), 解释方差为88.6%, 捕获深度为489.3 m。这说明TRW是事件一深层季节内振荡现象的主要原因。

探究事件一TRW的激发源, 首先检查此时间段内不同水层的涡旋情况。图 9给出了事件一期间代表性流速异常场的水平分布, 结果显示事件一时间段内确有一个次表层气旋涡经过, 该涡旋分布在541~ 1 941 m, 其中涡心在1 452~1 941 m相较于541~ 1 062 m向西北方向倾斜。图 10a和10b分别给出了2017~2019年次表层541 m与深层2533 m相对涡度[公式(5)]的时间变化和两者滑动相关系数的时间变化。滑动相关系数计算以每个时间点为中心, 涵盖前后30 d共60 d的数据来进行(下同)。结果显示在事件一时间段内, 次表层和深层的相对涡度均为正, 两者的相关系数高达0.86。两者的关系可由位涡守恒[公式(6)]进行解释, 次表层的气旋涡引发上升流, 如图 6a所示引起等温线的上升, 从而造成深层层厚h的增加, 在同一地点f不变的情况下, 根据位涡守恒深层的相对涡度也应增加, 两者的变化应相一致。上述分析说明该事件中的深层TRW是由次表层涡旋通过位涡守恒激发。另外需要指出的是在2017~2019年时间段内, 表层没有强涡旋经过YMJ-E站位并能够激发深层TRW过程。

2.3.2 事件二

首先使用TRW波解方程分别对事件二时间段的水平(图 5c)和垂直流速(图 6d)振幅进行拟合。水平流速振幅的拟合方程解为U(z)=6.8×10^-5cosh(0.002z), 解释方差为94.8%, 捕获深度为502.1 m。垂直流速振幅的拟合方程解为w_rms(z)=1.8×10^-7sinh(0.0024z), 解释方差为81.3%, 捕获深度为413.5 m。这说明TRW也是事件二深层季节内振荡现象的主要原因。

探究事件二TRW的激发源, 同样检查此时间段内不同水层的涡旋情况。图 11给出了事件二期间代表性流速异常场的水平分布, 结果显示事件二时间段内YMJ-E站点附近上层并无涡旋产生, 直到中层1 245 m有一个反气旋涡经过YMJ-E站点, 该涡旋一直延续到1 941 m。图 10a和10c分别给出了2017~2019年中层1 684 m与深层2 533 m相对涡度量值的时间变化和二者滑动相关系数的时间变化, 结果显示事件二时间段内, 中层和深层的相对涡度均为负, 两者的相关系数高达0.68。两者的关系可由位涡守恒[公式(6)]进行解释, 中层反气旋涡引发下降流, 如图 6a所示引起等温线的下降, 从而造成深层层厚h的减小, 在同一地点f不变的情况下, 根据位涡守恒深层的相对涡度也应减小, 两者的变化应相一致。上述分析说明该事件中的TRW是由中层涡旋通过位涡守恒激发。

2.3.3 事件三

我们首先使用TRW波解方程分别对事件三时间段的水平(图 5d)和垂直流速(图 6e)振幅进行拟合。水平流速振幅的拟合方程解为, 解释方差为97.2%, 捕获深度为400.6 m。垂直流速振幅的拟合方程解为, 解释方差为87.9%, 捕获深度为411.0 m。这说明TRW也是事件三深层季节内振荡现象的主要原因。

事件三时间段内上层和其他水层流速异常水平分布图均显示大洋内部没有涡旋经过(图略), 这说明涡旋不是TRW的激发源。下面我们使用1.2.4节涡能方程来计算季节内总能量时间导数和深层平均流正压、斜压不稳定性的时间变化(图 12)。结果显示正压不稳定绝对值一般大于斜压不稳定的绝对值, 说明正压不稳定的贡献占优; 事件三时间段内的季节内总能量的时间导数、正压和斜压不稳定项均为正值, 前两项在该时间段内数值较大且变化趋势较为一致, 这说明事件三中平均流主要通过正压不稳定过程向季节内能量转化, 为TRW提供了能量。此外, 2018年9~11月和2017年9~11月期间的情况与事件三相类似, 平均流也会通过正压不稳定过程向TRW提供能量。在事件一中, 正压和斜压不稳定值均为正值, 无法解释季节内总能量的减小。在事件二中, 正压不稳定值为负值, 斜压不稳定值很小, 无法解释季节内总能量的增加。因此, 在事件一和二中, 平均流的正斜压不稳定过程很可能不是TRW的主要能量来源, 这与前面的结论相一致。

3 结论

雅浦-马里亚纳海沟连接区(YMJ)东深水通道是太平洋深层经向翻转环流携带南极绕极水进入西太平洋的首要通道, 该通道流速和温度的多时间尺度变异对西太平洋深层动力状况、深层水更新等具有重要意义。本文基于全球高分辨率同化模式数据, 发现了该通道深层流和温度存在着40~90 d周期的季节内变异, 通过建立全水体的动力关联揭示了深层季节内变异的动力机制。地形罗斯贝波(TRW)是引起YMJ东通道季节内变异的主要动力机制, TRW波解方程可以解释大部分水平和垂直流速振幅垂直分布。TRW引起深层流振幅随深度增加而增加的结构, 在深层4 405 m处季节内流速振幅达23 cm/s, 1.45 ℃等温线振幅可达800 m。地形罗斯贝波在不同时间具有不同的激发机制, 首先, 当次表层有强涡旋经过时, 次表层强涡旋可以通过位涡守恒激发深层的TRW; 第二, 当大洋中层有强涡旋经过时, 中层强涡旋也可通过位涡守恒激发深层的TRW; 第三, 深层平均流可以主要通过正压不稳定向深层季节内TRW提供能量。