本文所用的福州站点(26.08°N,119.28°E)月均降水稳定同位素数据来自全球大气降水同位素网络(Global Network of Isotope in Precipitation,GNIP),数据起止时间为1985—1992年。同位素δ值均以维也纳标准海洋水(V-SMOW)的千分差表示:

δ¹⁸O(or δD)=[(R_水样/R_V-SMOW)-1]*1000‰

式中:R_水样、R_V-SMOW分别为样品与V-SMOW中的¹⁸O/¹⁶O及D/H比值。

福州站点的气象数据来自中国气象数据网,OLR数据源自国家卫星气象中心提供的风云二E地球同步轨道气象卫星,可见光红外自旋扫描辐射仪所处理的月均向外长波福射(OLR)数据,空间分辨率为1°*1°。卫星定位于105°E,能更好的反映我国及周边区域的水热状况及ENSO背景下水汽源区的演变过程。

本文用Niño 3.4区海洋尼诺指数(Oceanic Niño index,ONI),即海表温度指数(Sea Surface Temperature,SST)距平的3个月滑动平均值来定义ENSO事件,参考南方涛动指数(Southern Oscillation Index,SOI)计算海气强度指数。依美国气候预测中心(CPC)标准,ONI指数连续5个月≥0.5℃,则定义为厄尔尼诺事件(El Niño),ONI指数连续5个月≤-0.5℃,则定义为拉尼娜事件(La Niña)。依照ONI与SOI累计值的标准差计算海气综合强度指数(Oceanic Niño & Southern Oscillation Composite Index,OSI):

OSI_n =Std(∑ONI)-Std(∑SOI)

将海气综合强度指数(OSI)作为划分ENSO事件强度的依据,参照陈虹颖等^[17]的分类方法,按表1标准将ENSO事件分为极强、强、中等、弱、极弱5个等级,字母E代表各阶段厄尔尼诺事件,字母L代表各阶段拉尼娜事件,1985—1992年共经历了2次厄尔尼诺事件和2次拉尼娜事件,划分结果见表2。

表1 ENSO事件的海气综合强度指标

表2 厄尔尼诺事件和拉尼娜事件统计

对福州站的δ¹⁸ O、降水、气温数据分别按厄尔尼诺年、正常年、拉尼娜年计算采样时段内多年月平均值,结果如图1所示。

总体上,福州站δ¹⁸ O值表现为夏半年较低冬半年较高的趋势,由1月至12月表现为先波动下降再急剧上升的过程。由各事件年来看,δ¹⁸ O值总体表现出厄尔尼诺年>正常年>拉尼娜年的特征,特别是在夏季,这一特征更为突出,证明降水中δ¹⁸ O值的水平敏感的响应着ENSO事件对季风气候的影响。由图1可见,厄尔尼诺年、正常年、拉尼娜年的气温差异并不显著,但各事件年降水的差异则十分明显,与各事件年δ¹⁸ O值的变化有着较好的对应。依据各事件年降水量数据,计算雨季4—10月各事件年降水量所占比例(图2),可见在前期(4—7月)厄尔尼诺年降水量所占比例稳定在40%左右,而在后期(8—10月),特别是9—10月则是拉尼娜年降水比例显著增多,因而厄尔尼诺年降水显著集中于雨季前期,而拉尼娜年降水在雨季后期占主导地位。由此可见,ENSO事件对大气降水δ¹⁸ O的影响主要通过影响降水的年内分布,特别是雨季降水的时间分布状况而实现。

图1 福州大气降水δ18 O、降水量、温度的多年平均变化

图2 雨季各事件年平均降水量对比

大气降水中δD和δ¹⁸ O之间的线性关系被称作大气降水线方程。因降水中的同位素分馏状况及比率变化受降水区自然环境要素的影响,大气降水线方程代表着稳定同位素与地理特征之间的独特关系,深刻的反映着不同地域独特的降水特征^[18]。依据78组月均数据,得到福州站厄尔尼诺年、正常年、拉尼娜年大气降水线方程为:

厄尔尼诺年:δD=8.03δ¹⁸O+12.87(R²=0.89)

正常年:δD=8.21δ¹⁸O+10.38(R²=0.92)

拉尼娜年:δD=8.68δ¹⁸O+16.80(R²=0.98)

图3为福州δD-δ¹⁸ O间的线性关系,可见厄尔尼诺年的斜率小于正常年,拉尼娜年的斜率与截距均大于正常年,表明ENSO事件对降水中稳定同位素的分馏状况有着规律性的影响,造成这一差异的原因可能是由于ENSO事件改变了福州降水的年内分配状况。厄尔尼诺年降水集中于雨季前期,是旱季向雨季转换的季节,蒸发旺盛,雨水在降落过程中受到了二次蒸发,局地再循环水汽在降水中所占比例也较高,从而使厄尔尼诺年的斜率与截距降低。拉尼娜年降水集中于雨季后期,水汽蒸发趋缓,局地再循环水汽在降水中的比例也随之下降,从而推高了拉尼娜年的斜率与截距。但这一结论与上海^[12]站的大气降水线方程不同,其厄尔尼诺年的斜率与截距在数值上大于非厄尔尼诺年。由此可见,ENSO事件对大气环流的影响极为复杂,其对大气降水线方程的影响因地区而异,这一规律还有待于进一步研究。与Yurtsever等^[19]获得的全球大气降水线方程(δD=8.17δ¹⁸O+10.56)相比,正常年斜率与截距与之极为相近,非正常年斜率与截距则发生轻微的反向变化; 与Liu等^[20]提出的中国降水线方程(δD=7.48δ¹⁸O+1.01)、柳鉴容等^[21]所报道的中国东部季风区降水线方程(δD=7.46δ¹⁸O+0.90)相比,斜率与截距均较高,与珠江流域大气降水线方程(δD=8.08δ¹⁸O+11.00)相近^[18]。综合来看,与全球、珠江流域大气降水线方程相近,而较中国各地大气降水线方程显著较高的斜率与截距,体现出暖湿性的气候特征。

图3 福州大气降水线方程

在我国季风区,降水同位素组成受环境因素的制约,在不同季节形成不同的氢氧稳定同位素特征^[22],其中温度与降水量对同位素分馏比率的影响最为显著。因全年尺度下的环境效应并不显著,依据福州站月均数据,分别计算各事件年旱雨两季的降水量效应与温度效应(表3)。

表3 福州大气降水δ18 O与气温、降水的关系

福州站δ¹⁸O的环境效应在不同事件年有着不同的表现,总体上雨季以降水量效应为主,旱季以反温度效应为主。在雨季,各事件年δ¹⁸O与降水量的关系都通过0.05显著性水平检验,证明雨季降水量效应占主导地位,正常年R²_δ-P高达0.89,拉尼娜年略有减少,而厄尔尼诺年则减弱至0.30,δ¹⁸O与温度的关系全部未通过显著性检验。这应是由于ENSO事件导致水汽源区湿度出现异常,从而影响到水汽中的稳定同位素分馏状况所致。在旱季,厄尔尼诺年与正常年δ¹⁸O与温度的关系通过显著性检验,且厄尔尼诺年R²_δ-T为0.69,拉尼娜年反温度效应则并不存在,各事件年的降水量效应也都并不存在。总体而言,ENSO事件对福州环境效应的影响主要体现在:厄尔尼诺年减弱雨季降水量效应而增强旱季反温度效应,拉尼娜年则削弱旱季反温度效应,这一规律应与水汽源区的年际差异有关。

在影响降水稳定同位素特征的诸多因素中,水汽源区变化首先影响着稳定同位素的分馏状况,水汽团运移过程与降水区环境特征都是在其基础上,对降水中的稳定同位素特征进行再塑造^[23]。就水汽源区变化的影响而言,ENSO事件无疑是最为显著的因子,对水汽源区的季节性演变、区域乃至全球的降水、气温等气候要素的异常有着密切关系^[24],并直接影响着水循环中稳定同位素的分馏状况与特征。由于降水同位素的水平及特征在一个较长的时期内都处于较为稳定的状态,综合考虑降水同位素与气象卫星数据的起止时间,因ENSO事件的多为跨年度发生,且转换时间一般为夏半年,选定2015—2016年、2014—2015年、2010—2011年作为厄尔尼诺年、正常年与拉尼娜年的分析时间,以7月、10月、1月、4月OLR场图代表水汽源区的年内变化状况。需要说明的是,2015—2016年厄尔尼诺事件为极强,其OLR场表现较拉尼娜年更为极端,但也更加凸显其变化规律,结果如附图1所示。

同月份在不同ENSO事件背景下,水汽源区的状况有着显著差异,其OLR值总体上呈现出厄尔尼诺年>正常年>拉尼娜年的特征,这与图2中各站点在不同年份背景下表现出的δ¹⁸ O值水平差异相一致,证明ENSO事件背景下水汽源区的变化对δ¹⁸ O值的年际变化有着很大的影响。

具体来看,7月代表的夏季是气温最高、降水较多的季节,此时厄尔尼诺年水汽源区的OLR场表现为高值区域,无明显的水汽中心,而正常年与拉尼娜年则建立了明显的ITCZ。正常年孟加拉湾与西太平洋水汽源连为一体,西太平洋OLR低值区程度更深,水汽辐合更强,而拉尼娜年则相对较弱,孟加拉湾低值区辐合更强,表明正常年与拉尼娜年水汽源区特征差异显著。这一时期各站点正常年与拉尼娜年降水量总体大于厄尔尼诺年,由沿海向内陆这一特点更趋明显,图1中的δ¹⁸ O值同样表现为正常年>拉尼娜年>厄尔尼诺年,与OLR场特征一致,这一时期也是OLR与δ¹⁸ O值的差异表现最为显著的时期。

10月代表的秋季是气温与降水逐步降低的阶段,此时水汽源区已经南移,副高已经控制华南地区,厄尔尼诺年与正常年ITCZ显著减弱,而拉尼娜年则十分明显。这一时期降水显著减少,δ¹⁸ O值则出现了拉尼娜年与厄尔尼诺年先上升再下降的现象,特别是在11月演变为拉尼娜年>正常年>厄尔尼诺年,与OLR场一致。此时拉尼娜年辐合带范围更广、更偏北且程度更深,海温相对较高且水汽运移路线相对于正常年与厄尔尼诺年更短,稳定同位素在水汽团运移过程中贫化程度更低,气温降低的同时δ¹⁸ O值下降并不多。

1月代表的冬季是是气温与降水最低的季节,10°N地区由副高控制,ITCZ退至赤道地区,厄尔尼诺年辐合中心极为微弱,而正常年与拉尼娜年则较为明显。这一时期拉尼娜年辐合带保持着范围更广、更偏北且程度更深的特点,但这一时期福州主要受北方大陆气团影响,水汽源区的变化与δ¹⁸ O的变化并不同步。

4月代表的春季是气温与降水逐步升高的阶段,此时副高减弱,ITCZ北移,厄尔尼诺年辐合中心逐渐加强,正常年与拉尼娜年较为明显,正常年ITCZ加深,厄尔尼诺年ITCZ至5月方才建立,降水也随之增多。这一时期是厄尔尼诺年ITCZ最强的时期,随着孟加拉湾水汽源区的发展,厄尔尼诺年δ¹⁸ O值逐渐上升,这一时期δ¹⁸ O值水平全部为厄尔尼诺年>正常年>拉尼娜年的状态。

基于以上对OLR场各事件年内变化规律的分析,为进一步确证水汽源区的位置及其在不同ENSO背景下的总体变化状况,依据各事件年12个月份的月均OLR场,得到各事件年的年均OLR场,并分别计算厄尔尼诺年/拉尼娜年与正常年的差值,结果见附图2。

从ENSO事件的发展过程来看,与正常年相比,厄尔尼诺与拉尼娜年ITCZ的变化规律呈现出相反的趋势。厄尔尼诺年ITCZ在4月开始建立并加深,至7月达到最盛,而10月则较弱; 正常年在7月和4月左右是ITCZ程度最深的阶段; 拉尼娜年ITCZ则在10月达到最盛,1月仍保持着较大的范围。因此,厄尔尼诺年ITCZ主要建立并加深在春夏季,且总体偏南,规模较小; 而拉尼娜年ITCZ则在秋季达到一年中的峰值,总体偏北,范围较大,不同事件年ITCZ的变化特征与福州站δ¹⁸ O的环境效应规律相吻合。厄尔尼诺年对降水和δ¹⁸ O值的影响主要体现在雨季前期,而拉尼娜年主要体现在雨季的后期,这与2.1与2.3中的分析结果一致。

厄尔尼诺年与正常年的OLR差值主要体现在西太平洋,印度洋水汽源区差值与范围相对较小,主要位于赤道附近地区,其OLR均值较正常年偏高10～15 W/m²。拉尼娜年OLR差值同样体现在西太平洋水汽源区,且范围更广,强度更大,其OLR均值较正常年偏低20 W/m²左右,而印度洋水源区的波动范围与强度相较于厄尔尼诺年都明显减弱。拉尼娜年水源区位置要更靠北,接近10°N左右,西太平洋水汽源区的波动范围与强度都远大于印度洋水源区,因此印度洋水源区较为稳定,西太平洋水汽源区位置与对流强度的变化是ENSO事件年中影响研究区降水稳定同位素变化的主要因子。