李奎平王海员杨 洋于卫东李俐俐

(1.自然资源部 第一海洋研究所,山东 青岛266061；2.青岛海洋科学与技术试点国家实验室 区域海洋动力学与数值模拟功能实验室,山东 青岛266237；3.中山大学 大气科学学院,广东 珠海519082)

大气季节内振荡(intra-seasonal oscillation,ISO)是热带地区的一类行星尺度(纬向1～3波)对流扰动系统,它具有明显的季节特性[1-2]。冬季时,ISO事件主要活跃于赤道区域,在印度洋和西太平洋暖池区表现出显著的向东传播。夏季时,ISO事件的向东传播特征减弱,而在南亚季风区存在显著的向北传播。ISO事件是衔接天气尺度和季节尺度的关键桥梁,它对热带气旋、季风系统、印度洋赤道急流、印度洋偶极子和ENSO事件等都具有重要影响[3-7]。

ISO事件包含了复杂的热力和动力过程,在很大程度上表现为一种水汽驱动模态[8-9]。同时,海表湍流热通量和大气辐射通量的反馈机制也是ISO事件产生与维持的关键过程[10-11]。考虑到这些因素都是大气湿静力能(Moist Static Energy,MSE)的控制项,开展MSE诊断研究成为揭示ISO事件内在机理的重要手段。已有一些研究开展了ISO事件向东传播阶段的MSE分析[12-14],发现MSE的季节内变化是ISO事件对流演化的重要前提,它在对流活动之前不断累积(充电),并在对流活动后期逐渐释放(放电)。对各MSE控制项的诊断表明,MSE异常基本上与辐射和湍流加热存在同时正相关,与垂向平流存在反相关,而水平平流则在东侧领先于MSE异常,成为引导ISO向东传播的重要条件。

以上研究大多关注了冬季ISO事件沿赤道的向东传播过程。实际上,夏季ISO事件具有更加复杂的传播特征,它除了沿赤道向东传播,还在亚洲季风区存在东北方向传播以及独立的北向传播[15]。ISO湿/干位相在亚洲季风区的交替向北传播正是夏季风活跃/中断的重要诱因[16]。已有研究指出[17],夏季期间ISO事件的向北传播阶段也存在MSE的充电和放电过程,并且强调了水汽水平平流的重要作用,但是并没有明确控制水汽水平平流的关键物理过程。本文选取孟加拉湾这一典型代表区域,利用大气再分析资料,通过进行MSE和水汽诊断分析,定量探讨夏季ISO事件向北传播阶段的MSE演变特征及关键物理控制过程。

1 资料和方法

1.1 资料

本文用到的资料包括:

1)美国海洋大气局(National Oceanic and Atmospheric Administration,NOAA)提供的大气顶对外长波辐射通量(Outgoing Longwave Radiation,OLR)[18],水平分辨率2.5°×2.5°,时间分辨率为1 d,时间跨度为1979—2016年。OLR是用于表征大气对流的参量,OLR越低表示对流越强。

2)欧洲中期天气预报中心(European Centre for Medium-Range Weather Forecasts,ECMWF)提供的ERA-Interim大气再分析资料[19],包括200～1 000 hPa(共23层)的三维风速、气温、比湿和位势高度,大气顶和海表的长波与短波辐射通量,海表感热与潜热通量。资料水平分辨率1.5°×1.5°,时间分辨率为1 d,时间跨度为1979—2016年。

1.2 方法

大气MSE的表达式为m=Lq+c p T+gz,式中,m表示MSE；L为潜热常数,取2.5×106J·kg-1；q为大气比湿；cp为定压比热,取1 004 J·K-1·kg-1；T为气温；g为重力加速度；z为高度。大气柱垂向积分的MSE变化表示为[20]

式中,〈〉表示200～1 000 hPa的质量积分；V为水平速度；w为垂直速度；p为压力；QT为湍流热通量,即海表感热与潜热之和；Qs为辐射热通量,即大气吸收的长波和短波辐射之和；-〈V·∇m〉称为水平平流项；称为垂向平流项。

为定量揭示控制水汽水平平流的关键物理过程,采用时间滤波将原始物理量A进行如下形式分解:A=Al+Am+As。其中Al为低频变率,通过70 d低通滤波获取,用以表征季风背景；As为高频变率,通过20 d高通滤波获取,主要包括天气尺度过程和准双周振荡等；Am为季节内变化,通过将原始值减掉高频和低频变率后获取,代表了周期20～70 d的ISO信号。据此,水汽水平平流项表达为3种尺度过程的相互作用:

图1 依据参考指数定义的ISO位相示意图Fig.1 Definition of ISO phases based on ISO index

本研究关注孟加拉湾夏季(5—9月)的ISO事件,通过选取典型ISO事件进行合成分析,揭示MSE的演变特征与机理。考虑到每次ISO事件的周期存在差异,采用了按照位相进行合成的方法。首先将孟加拉湾中部(85°～95°E,5°～15°N)OLR进行20～70 d带通滤波,然后将5—9月的时间序列进行标准化处理构建了ISO参考指数。当参考指数低于-1时表示强事件发生,选定该事件作为分析样本,共计104个。针对选定的ISO事件,以ISO参考指数达到极小值时为0°,紧邻的前后两个极大值分别定义为-180°和+180°,指数跨越零点的时刻分别定义为-90°和+90°(图1),其他每10°间隔的位相通过内部插值获得。统计显示,104个ISO事件的平均周期为38 d,因此每10°间隔大致代表1 d。最后,将其它物理量进行20～70 d带通滤波得到其季节内变化信号,并参照上述定义的ISO位相对104个典型事件进行合成。

2 结果与分析

2.1 湿静力能的季节内变化

首先揭示夏季时MSE的季节内变化特征。图2为OLR表征的ISO对流和MSE异常的位相-空间演化,仅填充了MSE异常超过99%显著性水平(t-检验)的区域,并以星号表示了OLR异常超过99%显著性水平的区域。

图2 孟加拉湾夏季ISO事件的周期演化Fig.2 Lifecycle of summertime ISO events in the Bay of Bengal

由图2可见,-180°位相,孟加拉湾被强烈的对流抑制过程控制,MSE呈现显著负异常。进入-150°位相,在赤道西印度洋出现大面积MSE正异常,ISO对流开始形成。随后,MSE正异常中心与ISO对流中心都不断向东移动并逐渐增强,而孟加拉湾的MSE负异常和ISO对流抑制过程则逐渐消褪。-30°位相,MSE正异常中心与ISO对流中心到达赤道东印度洋,随后开始向北转向进入孟加拉湾；与此同时,赤道西印度洋开始出现MSE负异常。0°位相,孟加拉湾MSE正异常与ISO对流达到最强,此后二者在向北移动过程中逐渐减弱。赤道附近的MSE负异常逐渐增强,在+30°位相出现了ISO对流抑制过程。此后,MSE负异常和ISO对流抑制过程基本上重复上述时空演变特征。总体来看,MSE异常与ISO对流信号在空间和时间上都体现出较为一致的演变特征,即ISO对流区域(对流抑制区域)伴随着MSE正(负)异常,并且二者在时间上也呈现出比较同步的增强和减弱。

夏季ISO事件在孟加拉湾主要表现为向北传播,图3为孟加拉湾(85°～95°E)OLR异常和MSE异常的经向-位相演化,仅填充了MSE异常超过99%显著性水平(t-检验)的区域,并以星号表示了OLR异常超过99%显著性水平的区域。图3清晰显示了ISO在孟加拉湾的向北传播阶段伴随的MSE累积和释放过程。在前期的对流抑制位相,MSE呈现显著负异常；伴随着MSE的不断累积,ISO对流逐渐增强,并且都存在显著的向北传播特征；0°位相时,MSE异常和ISO对流都达到极大值；随后,MSE逐渐释放,ISO对流也不断减弱并被对流抑制过程取代。在ISO事件周期演变中,OLR的变化幅度接近50 W·m-2,MSE的振幅大致在16×106J·m-2。在北移过程中,ISO活跃中心大致位于5°～15°N,以下我们将针对此孟加拉湾中部区域(85°～95°E,5°～15°N)开展 MSE诊断分析。

图3 孟加拉湾断面OLR和MSE的经向-位相演化Fig.3 Meridional-phase evolution of OLR and moist static energy along the Bay of Bengal

2.2 湿静力能诊断分析

以上观测事实表明,ISO事件中MSE的累积(释放)是大气对流发展(衰减)的重要条件。MSE的变化受到平流输送、湍流加热和辐射加热等多种物理过程的影响,那么是何种过程主导了ISO北传阶段MSE的季节内变化?我们依据式(1)对MSE的各控制项进行了分析(图4)。从MSE时间变率(∂m/∂t)来看,它在-180°位相时转为正值,表示MSE开始累积(充电),从OLR反映的大气对流来看,它也开始逐渐增强；∂m/∂t在-90°位相时达到最大,表示此时的MSE累积速率最快；随后∂m/∂t逐渐减小,并在0°位相时转为负值,表示MSE充电结束,大气对流达到了最强；此后,∂m/∂t保持负值,MSE逐渐释放(放电),大气对流也逐渐减弱；∂m/∂t在+90°位相时达到最小,表示此时的MSE释放速率最快；+180°位相时,∂m/∂t开始趋近零点,MSE基本结束放电过程,大气对流达到最弱。从MSE各控制项来看,水平平流项(-〈V·∇m〉)的演化领先于∂m/∂t,它早在-180°位相时已体现出显著的充电效应,并维持高值至-90°附近,随后它快速降低并在-40°转为负值,开始导致 MSE放电。辐射热通量(QS)和湍流热通量(Q T)的演化显著滞后于∂m/∂t,它们分别在-80°和-60°位相转为正值,呈现充电效应,并在0°位相附近达到极大值；而在MSE的初期放电阶段(0～+90°),它们依然保持着充电效果。垂向平流项(-〈w·∂m/∂p〉)基本上与辐射/湍流热通量呈现反相关,但是它对MSE变化的贡献相对较小。在MSE季节内演变中,水平平流项、垂向平流项、湍流热通量和辐射热通量的振幅分别为60,11,31和38 W·m-2,这4个分项之和(黑色虚线)与∂m/∂t在时间演化和量值方面都具有较好的一致性,说明了式(1)在本研究中具备适应性。

图4 孟加拉湾中部MSE控制项的位相演化Fig.4 Phase evolution of moist static energy budget terms in the central Bay of Bengal

为清晰对比ISO事件中MSE累积和释放的控制机理,我们将-180°～-10°和+10°～+180°分别定义为充电阶段和放电阶段,并计算各MSE控制项在这两个阶段的贡献大小(图5)。比较发现,在MSE充电阶段,仅有水平平流项体现为显著的正贡献,垂向平流项的贡献接近为零,湍流热通量和辐射热通量则体现为负贡献。类似地,在MSE放电阶段,水平平流项也具有最大贡献,垂向平流项作用较小,湍流热通量和辐射热通量表现为反作用。以上结果说明在孟加拉湾夏季ISO事件中,MSE的充电和放电过程是由水平平流项所主导。事实上,MSE的变化是由水汽、气温和位势3个要素共同决定的,很多研究指出水汽的变化占据主导地位[12-14,21],从表1给出的各组成要素对MSE水平平流的贡献比例可以看出,夏季MSE季节内变化同样是由水汽所决定,这与之前的研究相一致[17]。

图5 孟加拉湾中部MSE各控制项在充电阶段和放电阶段的对比Fig.5 Comparisons of moist static energy budget terms between recharge and discharge phases in the central Bay of Bengal

表1 充电和放电阶段各组成要素对MSE水平平流的贡献Table 1 Percentage contributions of moisture,air temperature and geopotential to horizontal moist static energy advection term in recharge and discharge phases

大气垂向积分的结果表明水汽水平平流(-V·∇q)是孟加拉湾夏季ISO事件中MSE充电和放电的关键控制因素,那么各气压高度上的物理控制过程是否一致?我们进一步分析了水汽水平平流及其2个分量的高度-位相演化(图6),可以看出它们都在MSE的充电(放电)阶段呈现显著的正(负)异常,并且水平平流项在对流层中高层(400～600 hPa)和底层(700～900 hPa)存在2个显著变化中心。通过对比纬向平流(-u∂q/∂x)和经向平流(-v∂q/∂y)两个分量的垂向分布,可以看出对流层中高层的水汽变化中心主要是由纬向平流所致,低层的水汽变化中心则主要由经向平流贡献。

图6 孟加拉湾中部水汽水平平流的垂向-位相演化Fig.6 Vertical-phase evolution of horizontal moisture advection in the central Bay of Bengal

水汽水平平流在对流层低层和中高层存在显著差异,为了明确水汽平流效应的潜在物理控制过程,我们采用时间滤波将比湿、纬向风和经向风分解为低频变率、季节内变率和高频变率,并按照式(2)计算了中高层(400～600 hPa)纬向平流和低层(700～900 hPa)经向平流的各控制项,它们在MSE充电和放电阶段的对比如图7所示。在对流层中高层(图7a),水汽纬向平流在MSE充电和放电阶段的显著差异主要取决于纬向风季节内变化与低频水汽纬向梯度的作用项(-um∂ql/∂x)；在对流层低层(图7b),水汽经向平流在MSE充电和放电阶段的显著差异则主要来源于低频经向风与水汽季节内变率经向梯度的作用项(-vl∂qm/∂y)。

针对对流层中高层和低层的两个关键水汽平流控制项,图8用于阐述它们的潜在物理机制。其中图8a～8c分别为400～600 hPa的低频水汽分布(q l,g·kg-1)、充电阶段和放电阶段的纬向风季节内异常(um,m·s-1),图8d～8f分别为700～900 hPa的低频经向风(vl,m·s-1)、充电阶段和放电阶段的水汽季节内异常(qm,g·kg-1)。夏季时孟加拉湾水汽呈现西低东高的分布(∂ql/∂x＞0,图8a),水汽纬向平流在MSE充电和放电阶段显著对比是由纬向风异常的反号所致。在MSE充电阶段,对流抑制过程导致孟加拉湾被东风异常控制(um＜0,图8b),它对低频水汽的输送作用使孟加拉湾水汽不断增强；在MSE放电阶段,活跃对流导致孟加拉湾被西风异常控制(um＞0,图8c),它对低频水汽的输送作用则会使孟加拉湾水汽不断减少。与水汽纬向平流类似,水汽经向平流在MSE充电和放电阶段的显著对比是由于水汽季节内异常的经向梯度反号所致。夏季时孟加拉湾低层被南风控制(vl＞0,图8d),在MSE充电阶段,孟加拉湾南部(北部)为水汽正(负)异常,呈现南高北低的水汽分布(∂qm/∂y＜0,图8e),它与低频背景风场的相互作用使得孟加拉湾水汽不断增强；而MSE放电阶段水汽异常的南北分布恰好相反(∂qm/∂y＞0,图8f),其与低频背景风场的相互作用使得孟加拉湾水汽不断减少。上述2项水汽平流控制过程最终显著调制了MSE的充电和放电,影响了ISO事件的生命演化过程。

图7 孟加拉湾中部水汽纬向平流和经向平流各控制项在充电和放电阶段的对比Fig.7 Comparisons of zonal and meridional moisture advection terms between recharge and discharge phases in the central Bay of Bengal

图8 纬向平流和经向平流关键控制要素的分布Fig.8 Distribution of critical variables that dominate the zonal and meridional moisture advection

3 结 语

夏季期间ISO事件活跃于孟加拉湾,并存在显著的向北传播。本文采用1979—2016年ERA-Interim大气再分析资料,通过对典型ISO事件进行位相合成,分析了夏季孟加拉湾MSE季节内变化的演变特征和控制机理。结论如下:1)ISO事件在孟加拉湾的向北传播与MSE的演化密切联系,ISO对流的增强(减弱)伴随着MSE的显著充(放)电过程；2)孟加拉湾MSE在充电阶段和放电阶段的显著对比主要是由水平平流项所控制,其中水汽要素占据主导地位；3)水汽水平平流在对流层中高层和低层的控制机理存在显著不同。中高层主要是由纬向平流过程所致,取决于ISO扰动纬向风对低频水汽的输送；低层主要源自经向平流过程,取决于低频经向风对ISO扰动水汽的输送。

除孟加拉湾外,夏季ISO事件在阿拉伯海和中国南海也存在显著的向北传播,当前结果在其他2个海区的适应性需要进一步的研究予以确认。本文强调了不同高度背景风场和水汽分布对MSE季节内演化的重要调控作用,这对理解当前模式中的ISO模拟偏差具有重要指示意义,即提高模式对季风背景的模拟能力或是缩小ISO模拟偏差的重要前提。值得说明的是,水汽水平平流的演化显著领先于ISO对流中心(图6),这是诱导ISO对流向北传播的关键因素,Jiang等[22]强调了对流层低层水汽平流过程的重要性,这里我们同时指出了中高层水汽平流的贡献及其不同的物理控制过程,从而增进了对ISO北传机制的理解。