鲁云龙 冯俊乔 , 贾 凡 , 王庆业 ,

(1. 中国科学院海洋环流与波动重点实验室 青岛 266071; 2. 中国科学院海洋研究所 青岛 266071; 3. 中国科学院大学 北京 100049; 4. 中国科学院海洋大科学研究中心 青岛 266071; 5. 青岛海洋科学与技术试点国家实验室 海洋动力过程与气候功能实验室 青岛 266237)

厄尔尼诺-南方涛动(El Niño and Southern Oscillation, ENSO)是热带太平洋海气作用最强的年际信号, 主要的暖冷ENSO 事件分别对应厄尔尼诺和拉尼娜。进入21 世纪之前, 关于ENSO 传统的较成熟理论主要着眼于热带西太平洋对ENSO 发生发展的影响。一般而言, 厄尔尼诺事件存在两个先决条件: 赤道西太平洋的暖水堆积(Wyrtki, 1985)和西风异常(Luther et al, 1983; McPhaden et al, 1992)。当条件满足时, 异常的暖水可以很容易的向东传播, 造成中东太平洋的异常变暖。这一变暖会通过大气遥相关对全球的天气和气候造成显著的影响。但是, 进入21 世纪后, ENSO 展现出了一系列不同于以往的特征, 最为显著的就是暖海温异常中心位于赤道中太平洋地区的中部型厄尔尼诺事件频繁发生(Ashok et al, 2007), 这对ENSO 传统理论提出了挑战。

除了热带太平洋海气变率对ENSO 产生重要影响外, 在过去十几年中, 许多研究表明, 热带以外的太平洋海气变率可能为ENSO 事件的激发提供一个重要的确定因素。Vimont 等(2001)最先提出季节性足迹机制, 认为北太平洋涛动作为北半球冬季大气变率的第二模态, 其南半支通过调节东北信风的强度, 从而改变潜热通量, 在副热带东北太平洋留下海表面温度异常(sea surface temperature anomalies, SSTA)足迹, 同时与信风相互作用, 在风-蒸发-海表面温度反馈(wind-evaporation-sea surface temperature (SST) feedback, WES, Xie et al, 1994)作用下向西南传播, 在赤道太平洋引起类ENSO 变率。Chiang 等(2004)提出北太平洋SST 经向梯度和热带辐合带(intertropical convergence zone, ITCZ)的变化存在一定联系, 在副热带北太平洋存在一个与大西洋经向模类似的模态, 称之为北太平洋经向模(North Pacific meridional mode, NPMM)。Chang 等(2007)研究表明, 大多数厄尔尼诺事件之前都有明显的北太平洋经向模正位相(东北太平洋暖异常)活动, 北太平洋经向模是北太平洋大气变率影响ENSO 的重要通道。还有一些研究表明, 北太平洋经向模能通过海气过程对中部型厄尔尼诺事件产生重要影响(Yu et al, 2010; Su et al, 2014a)。Su等(2014b)指出2012 年秋季厄尔尼诺事件增暖的突然中断与北太平洋经向模的活动有关, 副热带东北太平洋受北太平洋经向模影响产生SST 冷异常, 其所导致的东风异常和低层辐散抑制了赤道中太平洋的对流活动, 使得2012—2013 年厄尔尼诺夭折。在年代际尺度上, Di Lorenzo 等(2015)的数值实验结果表明太平洋低频变率与NPMM-ENSO 的联合作用相关。

Zhang 等(2014)利用海气热力学耦合模式研究了南太平洋经向模(South Pacific meridional mode, SPMM)暖异常海温信号在风-蒸发-海表面温度反馈下传到热带太平洋的过程, 并提出云-SST 反馈机制在其中的作用。2014 年初, 许多专家预测年底会出现强度堪比1997—1998 年厄尔尼诺的ENSO 暖事件, 然而, 到了夏季暖异常增长停滞; 有研究发现, 在该事件过程中, 与南太平洋经向模有关的副热带东南太平洋冷SSTA 与赤道西太平洋的持续增暖相配合, 减弱了赤道太平洋纬向温度梯度, 增强了赤道中东太平洋的东南信风, 抑制了赤道中东太平洋的海气相互作用, 使得厄尔尼诺发展受到阻碍(Min et al, 2015)。Zheng 等(2018)发现南太平洋经向模东北极子区域的异常海温信号向赤道传播过程中, 风-蒸发-海表面温度反馈不是必须的, 信风的变化也可引起赤道太平洋次表层变化, 从而引起ENSO 事件。

在ENSO 预报方面, Larson 等(2014)肯定了北太平洋经向模作为ENSO 预报因子的有效性, 但也发现其预报技巧很低, 并且Larson 等(2018)提出南太平洋经向模的另一个角色, 不是作为ENSO 的触发因子, 而是作为ENSO 振幅和不确定性的热力学驱动因子。有研究表明, 南北半球的经向模态与不同类型ENSO之间具有重要联系, 南太平洋经向模有利于东部型厄尔尼诺, 北太平洋经向模更有利于中部型厄尔尼诺的发生(Zhang et al, 2014; Lin et al, 2015)。Min 等(2017)也指出南太平洋经向模主要有利于东太平洋SSTA 发展, 北太平洋经向模主要有利于赤道中太平洋SSTA 发展, 并据此可以参考这两个经向模态对两类ENSO 进行预测。

综上可见, 近些年来, 已有学者对北太平洋经向模和南太平洋经向模同ENSO 的关系进行了探索, 但同时对比分析两半球副热带太平洋经向模和热带的关系的研究比较少, 且尚无一致结论。本文将系统分析南北太平洋两个副热带经向模和ENSO 的关系, 深入理解副热带和热带海气相互作用的过程和机制, 并试图建立利用副热带变量预测ENSO 的统计模型。

1 数据与方法

本文所使用的月平均海表面温度资料来自英国气象局Hadley 中心的海冰和海表面温度数据集(HadISST)。月平均海表面高度(sea surface height, SSH)资料来自海洋再分析系统ORAS4(ocean reanalysis system 4), 因为海表面高度涨落是衡量温跃层深度扰动的一种很好的方法, 所以可以用它来代表温跃层的变化(Rebert et al, 1985)。使用的大气资料来自美国国家环境预报中心发布的再分析数据(NCEP-NCAR Reanalysis 1, Kalnay et al, 1996), 其中包括海表面月平均风场和热通量。本文分析的变量期间为1948—2018 年。

首先对所有数据进行线性去趋势处理, 去倾后的各变量减去其季节循环得到变量异常值。对异常序列进行3 个月的滑动平均。

使用的主要分析方法有最大协方差分析(maximum covariance analysis, MCA)、超前滞后相关分析、回归分析等。对相关分析的结果检验采用了Student’s t-test 方法。

2 结果分析

2.1 副热带南、北太平洋经向模

为了获得北太平洋经向模的时空分布, 对10m 风异常和SSTA 的协方差矩阵进行最大协方差分析, 由于热带地区更强的信号会掩盖副热带的特征, 所以在分析之前, 线性去除了冷舌指数(6°S—6°N, 180°—90°W平均SSTA, CTI, Deser et al, 1990), 从而将ENSO 信号剔除。在北太平洋(0°—35°N, 175°E—95°W)对余下数据进行最大协方差分析, 第一模态的方差贡献为72.5%, 相应的时间序列分别标记为 NPC(north principal components)-Wind 和NPC-SST。为了探索NPC-Wind 和NPC-SST 的关系, 计算了两个系数的超前滞后相关(图2a), 结果显示当风超前SST 一个月时, 相关系数最大, 为0.75, 这表明风驱动海洋。

图1 副热带北、南太平洋经向模态(a, b)及其时间序列(c, d) Fig.1 North and south Pacific meridional modes and the time series

如图1a 所示, 分别将10m 风、SSTA、潜热通量和海表面气压异常(sea level pressure anomalies, SLPA) 回归到标准化的NPC-Wind 上, 得到北太平洋经向模的空间结构。北太平洋经向模态的物理解释与已经建立的大西洋经向模态类似。图1a 显示北太平洋经向模的海温正异常中心位于(150°—125°W, 12°—25°N), 负异常中心位于(180°—160°W, 25°—32°N), 定义该结构为北太平洋经向模正位相。除SSTA 呈现偶极子结构外, 海表面气压异常也表现为南北偶极子结构, 其南极(中心位于160°—140°W, 25°—40°N)通过调节副热带高压系统, 使副热带东南太平洋出现异常西南风, 这就抑制了背景东北信风, 对应负潜热通量异常, 潜热释放减少, 有利于正SSTA 发展。

在副热带南太平洋(12°—38°S, 180°—70°W), 对10m 风异常和SSTA 进行相同的最大协方差分析, 第一模态的方差贡献为44.4%, 相应的时间序列标记为SPC-Wind 和SPC-SST。超前滞后相关分析显示(图2d), 与北太平洋经向模相同, 当风超前SST 一个月时相关系数达到最大0.54。同样的, 图1b 显示了10 m异常风场、SSTA、潜热通量异常和海表面气压异常相对于标准化的SPC-Wind 的回归场。海温正异常中心位于(100°—70°W, 20°—35°S), 负异常中心位于(140°—120°W, 32°—45°S), 定义为南太平洋经向模的正位相结构。在中高纬度存在一个海表面气压异常的偶极子结构, 类似于北太平洋经向模的情形, 副热带高压系统的变化调节背景信风, 从而改变潜热通量, 影响海表面温度。

为了研究两个经向模的时间特征, 分别计算了SST 和风指数的季节性变率。由图2b、c 可知北太平洋经向模在北半球冬季(12 月至来年3 月, DJFM)大气变率最强, SST 在冬春季(1—5 月, JFMAM)变率最强。而南太平洋经向模有所不同, 图2e、f 显示, 南太平洋经向模在南半球冬季(5—8 月, MJJA)大气变率最强, 而SST 最大变率发生在南半球秋夏(11 月至来年3 月, NDJFM)。南太平洋经向模风和SST 季节性的较大差异很可能是受副热带东南太平洋混合层深度(mixed layer depth, MLD)的季节性的影响。副热带东南太平洋混合层在南半球冬季较深夏季较浅(Ding et al, 2014), 浅的混合层可能会放大SST 对风的响应, 相反, 更深的混合层会抑制SST 对风的响应。

图2 海表面温度和风指数的超前滞后相关以及风和海表面温度指数的季节性 Fig.2 The lead-lag correlations between PC-SST and PC-Wind, and the seasonality of the PC-SST and PC-Wind

2.2 副热带南北太平洋经向模与ENSO 的关系

为了探索北太平洋经向模、南太平洋经向模与ENSO 的可能关系, 计算了NPC-Wind、SPC-Wind 与Niño3.4 指数的滞后相关, 由于ENSO 的非线性可能会对SSTA 产生影响, 所以以下的分析都是基于风指数(Chiang et al, 2004), 并且风指数包含内部大气变率的信息, 且与 ENSO 密切相关的海洋开尔文波(Kelvin)的产生有关(Thomas et al, 2016)。

图3a 中可以看到, 当Niño3.4 超前NPC-Wind 时, 二者基本没有相关关系, 说明ENSO 对北太平洋经向模事件发生的贡献不大; 但是, 当 Niño3.4 滞后NPC-Wind 8—11 个月时, 二者达到99%的置信度检验, 滞后10 个月时, 相关系数最大, 为0.45。这说明, 北太平洋经向模能够激发热带太平洋ENSO 事件, 北太平洋经向模正位相事件成熟10 个月后, EI Niño 事件可能会达到成熟位相。从图 3b 可以看到, 当Niño3.4 超前SPC-Wind 时, 二者为负相关, 但未达到95%的置信度检验; 当Niño3.4 滞后SPC-Wind 时, 二者为显著正相关, 滞后6 个月左右达到最大, 为0.34, 超过95%置信度检验。这意味着, 如果在南太平洋发生南太平洋经向模正位相事件, 6 个月后, 热带太平洋可能会发生 EI Niño 事件。NPC-Wind、SPC-Wind 和ENSO 的关系表明, 副热带南、北太平洋经向模对ENSO 的激发有着重要的影响, 两个经向模态可能是诱发ENSO 的原因之一。下面我们将利用回归分析讨论副热带和热带联系的可能机理。

2.3 北太平洋经向模与ENSO 联系的可能机制

图4、5 展示了SSTA、海表面风场、海气界面热通量(潜热通量和短波辐射通量, 由于感热和长波辐射通量相对贡献较小, 从而没有展示)、海表面高度异常(sea surface height anomalies, SSHA) 相对于NPC-Wind 2—4 月(FMA)平均指数的滞后回归场。

图3 北(a)、南(b)太平洋经向模态时间序列与Niño3.4 指数的滞后相关 Fig.3 The lead-lag correlation between NPC-Wind (a), SPC-Wind (b), and Niño3.4

图4 3 个月平均的SSTA、10 m 异常风场相对于NPC-Wind 2—4 月(FMA)指数的回归 Fig.4 The lag regressions of 3-month averaged SSTA and surface wind anomalies at 10 m on February-March-April (FMA)-averaged NPC-Wind

由图4a 可以看出, 在1—3 月(JFM)期间, 副热带北太平洋存在大尺度的气旋式环流, 北太平洋经向模暖极子区所在的副热带北太平洋的东南海区首先出现西南风异常, 削弱背景场的东北信风, 使得潜热释放 减少, 而冷极子区信风加强潜热释放增加(图5a1), 可以看出, 这时大气变率最强, 对应的潜热通量的变化最大, 这一过程一直持续到北半球春季, 北太平洋经向模的海温异常达到最大。另外, 在北太平洋经向模冷极子区的西北黑潮延伸体区, 海表面高度出现显著负异常, 表明次表层冷水可通过海水上升运动夹卷至混合层, 促进SSTA 冷异常发展(图5a3)。

北太平洋经向模成熟后, SSTA 暖异常伴随着异常西南风向西南扩展, 这一过程称为风-蒸发-海表面温度反馈。在 5—7 月(MJJ)以及 7—9 月(JAS)期间, 赤道中北太平洋西风异常加强; 赤道西太平洋的新几内亚沿岸区域负海表面温度异常加强, 对比各个通量项, 可以发现这主要是由潜热通量变化引起的。随着时间发展, 赤道中太平洋暖 SSTA 不断增强, 并向东发展, 在赤道太平洋出现Matsuno-Gill 型大气响应(Matsuno, 1966; Gill, 1980), 西太平洋赤道两侧出现气旋式环流异常, 呈现大气罗斯贝波(Rossby)调整,赤道西太平洋盛行异常西风, 并激发东传的海洋暖开尔文波(图5c3, d3), 有利于赤道中东太平洋地区温跃层深度增加和SSTA 增暖。同时, 赤道中太平洋地区短波辐射加强, 有利于SSTA 增长(图5c2, d2)。

在 9—11 月(SON)和 11—1 月(NDJ)期间, 赤道西太平洋西风异常不断向东推进, 在赤道东太平洋, 在SST-温跃层反馈机制作用下, 次表层过程影响 SSTA变化; 同时从海表面高度异常图中可以看到赤道海洋开尔文波持续东传, 最终导致厄尔尼诺发展并达到成熟。在这一期间从短波辐射变化导致的热通量回归场可以看到, 向下短波辐射增强, 促进暖SSTA 发展。

图5 3 个月平均的潜热通量、短波辐射和海表面高度异常相对于NPC-Wind 2-4 月(FMA)指数的回归Fig.5 The lag regressions of 3-month averaged latent heat flux, short wave radiation heat flux, and sea surface height anomalies on February-March-April (FMA)-averaged NPC-Wind

2.4 南太平洋经向模与ENSO 联系的可能机制

类似于北太平洋经向模, 图 6、7 展示了各个变量场相对于SPC-Wind 2—4 月(FMA)平均指数的滞后回归场。

图 6 可以看出, 在 1—3 月(JFM)期间, 副热带东南太平洋有大尺度的气旋式环流异常, 正SSTA 上空存在异常西北风, 削弱了背景场的东南信风, 潜热释放减少。相较于暖极子, 冷极子范围较小, 因此暖极子可能在联系南太平洋经向模与 ENSO 中起更重要的作用; 澳大利亚东南海域背景西风的加强使得潜热释放增加, 产生了冷的 SSTA。这一过程一直持续到春季, 同时东南太平洋气旋式环流异常加强, 暖异常信号在风-蒸发-海表面温度反馈下向西北延伸。这一热力学耦合过程也发生在澳大利亚东岸, 负 SSTA向西北延伸。赤道西风异常不断增强。

5—7 月(MJJ)和 7—9 月(JAS)期间, 从海表高度异常场的回归图中可以看到(图7c3, d3), 沿着赤道的异常西风分量激发向东传播的开尔文波, 促进赤道中东太平洋暖SSTA 发展。同时, 东南太平洋暖异常信号沿着南美洲西海岸继续向西北扩展, 加强了中东赤道太平洋SSTA 发展。在9—11 月(SON)和11—1月(NDJ)期间, 由于开尔文波的到来, 变暖进一步加强, 厄尔尼诺模态不断发展。从短波辐射通量的回归场可以发现, 赤道中太平洋地区短波辐射加强, 有利于厄尔尼诺模态的成熟(图 7e2, f2)。

2.5 回归模型

图 3 超前滞后相关分析结果表明, 在副热带太平洋出现北太平洋经向模正位相事件或者南太平洋经向模正位相事件后, 热带太平洋可能会发生厄尔尼诺事件。该结果为我们提前预测ENSO 事件提供了一个新的途径。据此, 我们利用月平均北太平洋经向模和南太平洋经向模的风场指数(二者关系不大, 相关系数约为0.2), 基于二元线性回归构建了如下预测模型:

其中, Y(t)为预测的 t 时刻 Niño3.4, X1(t-t1)为 t-t1时刻NPC-Wind, X2(t-t2)为 t-t2时刻 SPC-Wind。为了获得提前时间较长且较准确的预测效果, 我们取不同的t1和t2进行组合, 通过比较发现, 当t1=10, t2=6, 即, 当NPC-Wind 和 SPC-Wind 分别提前 10 和 6 个月时, 效果最好。利用最小二乘法计算了研究时段内的多元回归系数, 获得拟合公式为:

图6 3 个月平均的SSTA、10 m 异常风场相对于SPC-Wind 2—4 月(FMA)指数的回归Fig.6 The lag regressions of 3-month averaged SSTA and surface wind anomalies at 10 m on February-March-April (FMA)-averaged SPC-Wind

图7 3 个月平均的潜热通量、短波辐射和海表面高度异常相对于SPC-Wind 2-4 月(FMA)指数的回归Fig.7 The lag regressions of 3-month averaged latent heat flux, short wave radiation heat flux, and sea surface height anomalies on February-March-April (FMA)-averaged SPC-Wind

1950—2018 年, 根据该模型拟合的 Niño3.4 和观测的 Niño3.4 之间的相关系数达 0.57(图 8), 远超过99%的置信度检验, 并强于单独使用NPC-Wind 或者SPC-Wind 预测的情形。有学者也曾提出联合南北太平洋海气变率预测ENSO 有望提升ENSO 预报技巧。比如, You 等(2018)利用CMIP5 模式和再分析资料基于 2—4 月平均的表征南太平洋经向模和北太平洋经向模的风场指数以及赤道暖水体积指数共同构建了预测整个太平洋海盆同年11—1 月平均SSTA 的统计模型; Ding 等(2014)基于由SSTA 表征的南太平洋四极子结构和北太平洋偶极子结构在1—3 月期间的平均指数, 构建了预测同年 11—1 月平均的Niño3.4 的经验模型。本文建立的二元回归模型基于受ENSO 非线性影响较小的风场指数, 且采用北太平洋经向模超前10 个月和南太平洋经向模超前6 个月的指数预测月平均Niño3.4, 考虑了南北两个模态与ENSO 关系的区别; 对月平均Niño3.4 的预测, 不仅能体现ENSO 成熟期的状态, 而且还能了解其时空演变过程。

图8 观测的Niño3.4 指数与模型预测指数Fig.8 Observed and model predicted Niño3.4 index

3 总结与讨论

本文利用观测和再分析资料, 分析了副热带南、北太平洋经向模的基本特征及其与ENSO 的关系。两个经向模态都与副热带高压系统的变化有着重要的联系; 副高引起的信风强度的变化, 可以通过调节海气界面潜热通量, 从而影响海表面温度。相关分析结果表明, 两个经向模态与ENSO 均具有密切联系。当Niño3.4 滞后北太平洋经向模 8—11 个月时, 二者相关系数超过 99%置信度检验, 滞后 10 个月时, 相关最大; Niño3.4 滞后南太平洋经向模5—10 个月时, 二者显著正相关, 滞后6 个月时, 相关最大。在副热带太平洋发生北太平洋经向模正位相事件(北太平洋东南暖, 西北冷)大约 1 年后或者南太平洋经向模正位相事件(南太平洋东北暖, 西南冷)大约半年后, 热带太平洋可能会发生厄尔尼诺事件。据此, 根据表征南太平洋经向模和北太平洋经向模的风场指数, 建立了预测ENSO 的回归模型, 效果良好。

利用回归分析探讨了副热带南北太平洋经向模与ENSO 联系的可能机制。结果发现, 海气热力学过程在联系二者中起主要作用, 同时海洋动力学过程也具有一定贡献。副热带北太平洋暖异常信号在风-蒸发-海表面温度反馈下向西南传播到赤道中太平洋, 在Matsuno-Gill 型响应下, 海气耦合过程使得厄尔尼诺事件快速发展并走向成熟; 南太平洋经向模事件发生后, 东南太平洋暖异常信号沿着南美西海岸向赤道传播, 另外, 澳大利亚东岸冷异常向西北延伸, 赤道太平洋异常西风不断增强, 有助于ENSO 的发展。

当前流行的 CMIP5 模式对热带太平洋背景条件比较敏感, 对热带外信号的捕捉存在很大不确定性(You et al, 2018)。此外, 在全球变暖的背景下, 目前的气候预测模式亦有很多不确定性(Wang et al,2019)。因此, 为了提高对ENSO 的预测, 不仅要重视热带海洋动力和热力过程, 更要重视热带外过程, 特别是最近才开始受到广泛关注的南太平洋。