杨崧 徐连连 2
1 中山大学大气科学学院/南方海洋科学与工程广东省实验室(珠海),珠海 519082
2 广东省气候变化与自然灾害研究重点实验室,珠海 519082
在全球变暖背景下,极端降水强度和频率的变化越来越显著,由此引发了洪涝、干旱和泥石流等自然灾害,严重威胁着人类的生命财产安全(严中伟和杨赤, 2000; 江志红等, 2007; 杨金虎等, 2008;Sillmann et al., 2013; 任福民等, 2014; 孙颖等, 2015;徐影等, 2017; Song et al., 2020; Zhou et al., 2021; 江洁等, 2022; Sun et al., 2022; Xu et al., 2023)。众多国内外研究指出,不同区域的极端降水对气候变暖的响应不同,其中泛南海地区是北半球极端降水变化最显著的区域之一,该区域每年由极端降水造成的经济损失超过总经济损失的40%,严重影响 了 当 地 居 民 的 生 产 和 生 活(Ding et al., 2004;Scoccimarro et al., 2013; Dong et al., 2021; Ge et al.,2021; Xu et al., 2021; Zhang et al., 2021)。
泛南海地区属于亚洲季风系统的关键区域,是亚洲季风与厄尔尼诺—南方涛动(ENSO)相互作用、太平洋与印度洋相互作用以及南北半球相互作用的交汇中心,该区域降水释放的潜热加热可以调节局地哈德莱环流直接影响周边地区的降水,也能够调节南亚高压的位置触发孟加拉湾夏季风的爆发,还可以激发罗斯贝波列影响北美和欧亚大陆东北部的冬季气温(何金海等, 2006; Liu et al., 2013; Li et al., 2016; Chen and Zhai, 2017; Jiang et al., 2017; He and Yang, 2018; 舒 建 川 等, 2019; 杨 崧 等, 2019;Dong and Wang, 2022)。与平均降水相比,泛南海地区的极端降水能够给局地带来更加严重的气象灾害,对周边地区甚至全球的天气气候造成严重影响(Scoccimarro et al., 2013; Xu et al., 2021, 2022)。因此,泛南海地区极端降水的时空变化特征及变异机理,一直是国内外学者关注的焦点。
由于技术条件有限等多种原因,针对泛南海的长期连续观测资料极为匮乏,难以满足南海环境安全保障、生态环境保护以及气候变化研究等需求。近年来在国家863 计划、国家海洋局专项以及国家自然基金等项目支持下,各类气象观测站点逐渐加密,使得泛南海地区的观测数据日益准确,极大地促进了该区域极端气候的研究。例如1994 年开展的“南海季风试验”,中国科学家联合印度尼西亚、马来西亚、菲律宾等多个国家的科学家在泛南海区域开展了大型的地面、探空等加密观测以及海洋浮标、科学考察船等特殊观测项目,所获得的数据在国内外相关研究中得到了广泛的应用。利用两套观测数据(APHRODITE 和CPC-UNI),两套卫星数据(TRMM 和GPCP1DD)以及三套再分析数据(ERA-Interim、ERA5 和JRA55),Kim et al.(2019)分析了1998~2007 年泛南海地区极端降水的空间分布特征,指出APHRODITE 和CPC-UNI 具有较好的一致性。APHRODITE 是基于亚洲地区12000多个观测站点插值成的格点数据,水平分辨率为0.25°×0.25°,时间长度为1951~2015 年(Yatagai et al., 2012)。CPC-UNI 是基于全球30000 多个站点插值成的格点数据,水平分辨率为0.5°×0.5°,时间长度为1979 年至今(Xie et al., 2007)。APHRODITE和CPC-UNI 经常作为参照数据评估模式数据在泛南 海 地 区 的 适 用 性(Ge et al., 2021; Utsumi and Kim, 2022)。
为了方便个人、地区和国家能够以完全相同的方式研究极端气候,气候变化与监测指标专家组(ETCCDI)定义了一系列极端降水指数(Karl et al., 1999)。基于不同时段和不同分辨率的日降水数据,研究者对泛南海区域极端降水的时空变化特征展开了大量的研究,结果表明泛南海区域极端降水的变化具有较大的区域性差异(Endo et al., 2009;Villafuerte and Matsumoto, 2015; 江洁等, 2022; Xu et al., 2022; 杨雯婷等, 2022)。例如,基于APHRODITE的日降水数据,Villafuerte and Matsumoto(2015)发现1951~2007 年东南亚地区年和季节最大降水量在中南半岛和菲律宾中东部地区呈现显著增加的趋势,在海洋性大陆大部分区域呈现减小的趋势。为了更好地刻画极端气候在区域尺度上的变化特征,Iturbide et al.(2020)将全球陆地区域分为49 个子区域,其中泛南海区域包含东南亚、南亚、青藏高原和东亚四个子区域(图1a)。截至目前,少有研究基于不同极端降水指数,全面分析泛南海不同区域极端降水的时空变化特征。
图1 1951~2014 年泛南海地区年平均(a、e)R95p(单位:mm)、(b、f)R20(单位:d)、(c、g)RX1day(单位:mm)和(d、h)RX5day(单位:mm)的空间分布型(左列)以及在四个子区域的平均值(右列)。(a)中的SEA、SAS、TIB 和EAS 分别表示东南亚、南亚、青藏高原以及东亚区域;(e-h)中浅蓝、深蓝、粉色和黄色直条分别表示春季、夏季、秋季和冬季平均值。Fig.1 Spatial patterns of annual mean (ANN) (a) R95p (units: mm), (b) R20 (units: d), (c) RX1day (units: mm), and (d) RX5day (units: mm) during 1951-2014 over the South China Sea and surrounding areas (SCSSA); SEA, SAS, TIB, and EAS in Fig.1a represent Southeast Asia, South Asia, the Tibet Plateau, and East Asia, respectively.Seasonal mean (e) R95p, (f) R20, (g) RX1day, and (h) RX5day averaged over four sub-regions of the SCSSA during 1951-2014; light blue, dark blue, pink, and yellow bars depict spring (MAM), summer (JJA), autumn (SON), and winter (DJF),respectively.
耦合模式比较计划第五和第六阶段(CMIP5和CMIP6)的多模式产品被广泛地应用于预估未来气候的变化。与CMIP5 相比,CMIP6 气候模式对中国区域极端事件的模拟能力有了整体的提升,但与观测仍然存在一定的偏差(Chen et al., 2017;Zhu et al., 2020; Yang et al., 2021)。美国国家航空航天局对CMIP6 气候模式的日产品进行误差订正并降尺度到0.25°×0.25°的水平分辨率,生成了高精度NEX-GDDP-CMIP6 气候数据集(Wood et al.,2002, 2004; Maurer and Hidalgo, 2008; Thrasher et al., 2012)。基于高精度NEX-GDDP-CMIP6 数据得到的泛南海地区极端降水在未来会呈现何种时空分布特征?该问题的认识尚不清晰。
东南亚是泛南海地区对流最活跃的区域,该地区对流活动能够显著影响我国南部的天气气候(何金海等, 2006; Liu et al., 2013; 杨崧等, 2019)。以往研究主要分析热带海温异常对东南亚平均降水年际变率的影响,对极端降水形成机制的研究较少。研究表明,ENSO 期间东南亚更容易发生极端降水,且当ENSO 和印度洋偶极子(IOD)同时发生时,东南亚极端降水强度更大(Villafuerte and Matsumoto,2015; Setiawan et al., 2017; Da Silva and Matthews,2021)。但热带海温异常如何影响东南亚极端降水,这个问题尚未被清晰地揭示。
综上所述,本文将主要解决以下两个科学问题:(1)历史和未来时期,泛南海地区不同极端降水指数会呈现何种时空分布特征?依据CMIP6 时段的划分,本文的历史时段为1951~2014 年,未来时段为2015~2100 年。(2)热带海温异常通过何种途径影响东南亚极端降水?这些问题的研究是深入认识泛南海地区极端降水变异机理的关键所在,对该区域极端降水的预测预警具有十分重要的指示意义。
本研究挑选ETCCDI 的四个极端降水指数(表1),分别是日降水的最大值(RX1day)、连续5 日降水的最大值(RX5day)、极端强降水天数(R20)和非常湿润天(R95p)。RX1day 常用于表征极端强降水的强度。RX5day 一般用于持续性强降水的特征研究。R20 是基于绝对阈值定义的指数,一般用于研究极端强降水的发生频率。R95p 是基于相对阈值定义的指数,表示极端强降水的累计降雨量。本研究利用APHRODITE 的日降水数据,分析1951~2014 年泛南海地区四个极端降水指数的时空变化特征。
表1 ETCCDI 中四个极端降水指数的名称、简写、定义以及单位Table 1 Names, abbreviations, definitions, and units of four extreme precipitation indices obtained from ETCCDI (Expert Team on Climate Change Detection and Indices)
图1 给出了1951~2014 年泛南海地区年平均R95p、R20、RX1day 和RX5day 的空间分布特征以及四个指数在各个子区域的季节平均值。对整个区域而言,年平均R95p、R20、RX1day 和RX5day的区域平均值分别为226.2 mm、7.7 d、37.0 mm和80.5 mm。年平均R95p 和R20、RX1day 以及RX5day 的空间分布相似,相关系数分别是0.94、0.93 和0.93,均通过99%信度水平的显著性检验,较大值分布在东南亚、中国东南部以及青藏高原南坡地区(图1a-d)。由此可见,泛南海地区极端强降水发生的区域,持续性降水的强度也高,极端降水的发生频次也大。
东南亚地区R95p 的季节性差异较小,其他区域R95p 在夏季最大,秋季和春季次之,冬季最小(图1e)。值得注意的是,夏季平均的R95p 在南亚和东亚的值大于东南亚,另外三个季节东南亚的R95p 最大。这是由于亚洲季风和澳大利亚季风系统的影响,泛南海地区的对流活动和降水在春季迅速发展,降水的大值中心从赤道向北移动,在夏季到达南亚和中国南部地区,秋季雨带逐渐向南移动,在冬季回到赤道附近(Zhang et al., 2016; Yang et al., 2019)。东南亚地区四个季节降水超过20 mm的天数大约为4 d。南亚、青藏高原和东亚地区夏季降水超过20 mm 的天数大约7 d、3 d 和5 d,这三个区域冬季降水超过20 mm 的天数几乎为零(图1f)。由此可见,东南亚一年四季都极易发生强降水,南亚和东亚是泛南海夏季高频强降水的发生地区。RX1day 和RX5day 的季节平均特征与R95p 相似(图1g-h),即东南亚强降水和持续性强降水事件的季节性差异较小,其他区域这两种极端事件在夏季最大,冬季最小。
1951~2014 年泛南海地区年平均R95p、R20、RX1day 以及RX5day 趋势的空间分布不均一,且大部分区域的趋势并不显著(图2a-d)。对整个区域而言,年平均R95p 在1951~2014 年呈显著增加的趋势(0.75 mm a-1,p<0.01)。春季、夏季和冬季平均的R95p 均呈显著增加的趋势,而秋季的增加趋势不显著(图2e)。与R95p 相似,除秋季平均以外,年平均和其他季节平均的RX1day在1951~2014 年呈显著增加的趋势(图2f)。年、春季和冬季平均的R20 和RX5day 在1951~2014年呈显著增加的趋势,这两个指数在夏季和秋季的趋势变化不显著(图2g 和h)。
综上所述,1951~2014 年泛南海地区年平均和季节平均的四个极端降水指数的空间分布相似,较大值分布在东南亚、中国东南部以及青藏高原南坡地区。东南亚四个极端降水指数的季节性差异较小,其他区域极端降水指数在夏季最大,秋季和春季次之,冬季最小。泛南海地区1951~2014 年四个极端降水指数趋势的空间分布不均一,且大部分区域趋势的变化不显著。对整个区域而言,年、春季和冬季平均的四个极端降水指数均呈显著增加的趋势,秋季的增加趋势不显著。不同指数在夏季的趋势变化存在差异,R95p 和RX1day 呈显著增加的趋势,而R20 和RX5day 趋势变化不显著。
CMIP5 和CMIP6 气候模式产品是预估未来气候变化的重要手段,但其对区域尺度极端降水的模拟与观测存在偏差(Chen et al., 2017; Zhu et al.,2020; Yang et al., 2021)。为了减小模式的偏差、提供更高水平分辨率的气候数据集,研究者们研发了多种统计和动力降尺度方法。与动力降尺度方法相比,统计降尺度方法需要的计算资源较少,已被广泛地应用于多个研究领域(Wood et al., 2002,2004; Maurer and Hidalgo, 2008; Thrasher et al.,2012; 魏培培等, 2019; Xu and Wang, 2019; Xu et al.,2023)。美国国家航空航天局利用误差订正—空间降尺度(BCSD)方法,将CMIP5 气候模式历史时期和未来情景下的日数据进行误差订正并降尺度到0.25°×0.25°的水平分辨率,生成了高精度NEXGDDP-CMIP5 数据集。与CMIP5 原始值相比,NEXGDDP-CMIP5 对区域尺度极端温度和极端降水的模拟与观测更加接近(Bao and Wen, 2017; Chen et al., 2017; Xu and Wang, 2019)。
BCSD 是基于分位数映射原理发展的统计降尺度方法,它的优势是能够将模式数据的平均值和标准差调整到和观测数据十分接近。例如Xu and Wang(2019)指出,NEX-GDDP-CMIP5 对中国区域日最高温度的平均值和标准差的模拟与观测基本一致(图3)。NEX-GDDP-CMIP6 和NEX-GDDPCMIP5 的生成原理相似,我们进一步基于NEXGDDP-CMIP6 数据集对东南亚的极端降水进行了预估研究。
基于20 个NEX-GDDP-CMIP6 气候模式在SSP1-2.6 和SSP5-8.5 情景下的日降水数据,我们分析了2015~2100 年泛南海地区R95p、R20、RX1day和RX5day 的时空分布特征(图4)。与历史时期相似,SSP1-2.6 情景下年平均的四个极端降水指数较大值分布在东南亚、中国东南部以及青藏高原南坡地区(图3a-d)。SSP5-8.5 情景下极端降水的空间分布也与历史时期相似(图略),且SSP5-8.5 情景下极端降水的值大于SSP1-2.6 情景。全面探索历史时期泛南海地区极端强降水和持续性强降水的时空变化特征并揭示其中的物理机制,对于该区域未来的防灾减灾和气候预警工作至关重要。除了R20 以外,R95p、RX1day 和RX5day 的信噪比在泛南海区域大于1。信噪比定于为多模式的集合平均值与多模式标准差的比,比值大于1 表明模式对泛南海地区极端降水的模拟具有较好的一致性。对整个区域而言,2015~2100 年年平均的四个极端降水指数均呈显著增加的趋势,且SSP5-8.5 情景下的增加趋势大于SSP1-2.6 情景(图4e-h)。因此,泛南海地区极端强降水、持续性强降水以及强降水的发生频次在未来均增加,且温室气体浓度越高,增加趋势越显著。
图4 SSP1-2.6 情景下2015~2100 年泛南海地区年平均(a、e)R95p(单位:mm)、(b、f)R20(单位:d)、(c、g)RX1day(单位:mm)和(d、h)RX5day(单位:mm)的空间分布(左列)以及随时间的变化(右列)。(a-d)中打点区域表示信噪比大于1,信噪比定义为多模式集合平均与多模式标准差的比值;(e-h)蓝线(红线)和蓝色(红色)阴影分别表示SSP1-2.6(SSP5-8.5)情景下20 个模式的平均值和标准差。Fig.4 Spatial patterns of (a) R95p, (b) R20, (c) RX1day, and (d) RX5day during 2015-2100 under the SSP1-2.6 scenario over the SCSSA.The stippling in (a-d) denotes that the signal-to-noise ratio is above 1, the signal-to-noise ratio is defined as the ratio of the NEX-GDDP-CMIP6 multimodel ensemble mean to the corresponding inter-model standard deviation.Time series of (e) R95p, (f) R20, (g) RX1day, and (h) RX5day averaged over the SCSSA during 2015-2100 under the SSP1-2.6 (blue lines) and SSP5-8.5 scenarios (red lines) and their uncertainties (shadings).
图5 给出了2016~2035 年、2046~2065 年和2080~2099 年泛南海四个子区域年平均极端降水指数的百分比变化,参考时段为1995~2014 年。百分比变化定义为研究时段与参考时段的均值之差除以参考时段的均值。在SSP1-2.6 情景下的未来三个时段R95p 在四个子区域变化不大,SSP5-8.5情景下南亚和青藏高原的R95p 增加,且随着时间的推移,增加幅度越大。南亚和青藏高原的R20在未来两种情景的三个时段增加,且SSP5-8.5 情景下的增加幅度大于SSP1-2.6 情景。东南亚的R20 在未来三个时段呈现微弱的增加,而东亚的R20 则轻微减小。在SSP1-2.6 和SSP5-8.5 情景下,RX1day 和RX5day 在东南亚和东亚减小,南亚和青藏高原增加。由此可见,泛南海地区各个指数的变化具有显著的区域性差异,但这四个指数均在南亚和青藏高原增加,在东南亚和东亚减小。大量研究表明,ENSO、局地海温异常、南亚高压的建立时间以及青藏高原热源均能显著影响南亚地区的环流和降水异常(Mao and Wu, 2007; Jiang and Li,2011; Wang et al.2019; Yu et al., 2021)。北极海冰、北大西洋海温异常、ENSO 以及海洋性大陆的对流活动等是影响青藏高原降水异常的重要因子(Yu et al., 2021, 2023)。但南亚和青藏高原的极端降水变化幅度大于东南亚和东亚地区的物理机制仍有待进一步研究。
图5 SSP1-2.6 和SSP5-8.5 情景下泛南海四个子区域年平均(a)R95p、(b)R20、(c)RX1day 和(d)RX5day 的百分比变化。黑色、蓝色和红色盒须图分别表示2016~2035 年、2046~2065 年和2080~2099 年相比于1995~2014 年的变化。百分比变化定义为研究时段与参考时段的均值之差除以参考时段的均值。盒须图展示的是20 个模式变化的10th 分位数、25th 分位数、平均值、75th 分位数和90th 分位数。Fig.5 Box-and-whisker plots of percentage changes in (a) R95p, (b) R20, (c) RX1day, and (d) RX5day during 2016-2035 (black boxes), 2046-2065(blue boxes), and during 2080-2099 (red boxes) with respective to 1995-2014 under the SSP1-2.6 and SSP5-8.5 scenarios averaged over the four subregions.The box-and-whisker plot shows the 10th, 25th, 50th, 75th, and 90th percentiles of the 20 models from the NEX-GDDP-CMIP6 dataset.
综上所述,虽然东南亚是1951~2014 年年平均的各个极端降水指数的大值区,但该区域四个极端降水指数在未来的变化没有其他区域明显。南亚和青藏高原是未来四个极端降水指数变化的大值区,且SSP5-8.5 情景下的变化幅度大于SSP1-2.6 情景。
泛南海地区极端降水的变化具有较大的区域性差异,尤其在东南亚区域,该区域降水释放的潜热加热可以调节局地哈德莱环流直接影响我国华南和西南等地区的降水(何金海等, 2006; Li et al., 2016;Jiang et al., 2017; 舒建川等, 2019)。由于复杂的地形分布、显著的海陆热力差异以及降水的多时空尺度变化特征,模式通常对东南亚地区降水的模拟性能较差(Yang et al., 2019)。深入理解该区域极端降水的变化特征及机理,可为改善模式对该区域降水的模拟能力提供理论支撑,也可为提高我国的短期气候预测水平提供有利的参考依据。
不同区域降水的变化不同,它们的气候效应也存在差别(Chang et al., 2005; Jin et al., 2017; Jiang and Li, 2018)。但以往研究大多将东南亚视为整体,研究区域平均降水的气候效应。简单的区域平均可能低估东南亚降水的年际变率信号,进而低估ENSO 对东亚夏季风的作用(McGregor, 1992;Aldrian and Susanto, 2003)。主模态的提取不仅可以涵盖东南亚所有区域,还能够弥补以往研究多侧重于区域平均而很少涉及降水空间差异性的不足。东南亚降水的空间分布十分复杂,常用的主分量分析方法(EOF/REOF)提取的主模态解释方差偏低(Aldrian and Susanto, 2003; Xu et al., 2019)。最近几十年,K-means、Ward’s method 和自组织映射(SOM)等聚类分析方法被广泛地应用于极端天气气候事件的研究领域。SOM 的优势是:当气候系统存在噪音或者分布十分复杂时,该方法仍能够精确客观地提取气候变量的主模态(Bao and Wallace, 2015; Li et al., 2015; Guo et al., 2019)。
与所有动态聚类方法类似,SOM 方法需要事先指定分类数目。太多的分类数会导致类别之间的物理差异不够显著,分类数太少又不能够完全体现变量的分布特征。伪发现率(FDR)是常用的确定SOM 最优分类数的场显著性检验方法(Johnson et al., 2008; Guo et al., 2019)。显著性水平p选取为0.05。当极端降水被分为K类时,利用FDR 方法分别对组极端降水合成场做显著性检验,其中K的取值从2 开始增加,当K=L时,如果出现场显著性检验的p值大于0.05,即表明有1 对极端降水合成场的差异不显著,由此将SOM 的最优分类数确定为L-1。
APHRODITE 数据的时间长度是1951~2015 年,为了分析最近几十年东南亚极端降水的形成机理,本文利用CPC-UNI 日降水数据。依据相对阈值方法分别挑选东南亚1979~2019 年春季、夏季、秋季和冬季的极端降水日。采用SOM 和FDR 方法对四个季节极端降水日的降水量进行聚类分析(表2)。结果表明,东南亚春季极端降水主要分为北湿南干和全区一致偏湿两种空间分布型,这两种空间分布型出现的频率为38.3%和61.7%。东南亚夏季极端降水主要呈现北湿南干(C1)、全区一致偏湿(C2)和北干南湿(C3)的空间分布特征(图6),且三种主模态的发生频次呈上升的趋势,其中C2 和C3 呈显著增加的趋势,但C1 的上升趋势不显著。值得注意的是,C1 和C3 型极端降水的空间分布及其对应的大尺度环流特征相反,但两种类型极端降水对应的时间序列相关系数不显著,说明C1 和C3 是两种不同类型的极端降水。与夏季极端降水类似,东南亚秋季极端降水主要分为北湿南干、全区一致偏湿和北干南湿三种空间分布型。东南亚冬季一致偏湿和北干南湿分布型发生的频率为75.1%和24.9%。
表2 1979~2019 年春季、夏季、秋季和冬季极端降水的SOM 分型及其频率Table 2 Spatial patterns of SOM modes for spring, summer, autumn, and winter extreme precipitaion over Southeast Asia during 1979-2019 and their frequencies.
图6 1979~2019 年东南亚夏季(a、d)北湿南干(C1)、(b、e)全区一致偏湿(C2)和(c、f)北干南湿(C3)型极端降水的空间分布(第一行)及三种分型对应的主模态发生频次时间序列(第二行)。第一行图右上角的百分数表示每种分型的发生频率,打点区域表示通过了95%信度水平的显著性检验。引自Xu et al.(2022)。Fig.6 Spatial patterns of 1×3 SOM nodes (C1, C2, and C3) for summertime extreme precipitation (units: mm d-1) over Southeast Asia (top panels)and their corresponding time series (bottom panels) during 1979-2019.The percentages in the parentheses of top panels represent the frequency of each pattern.Stippling indicates the anomalies that are statistically significant at the 95% confidence level based on the Student’s t-test.From Xu et al.(2022).
利用回归分析、合成分析、偏相关和偏回归等多种统计分析方法,结合多种诊断分析方法,我们进一步分析了1979~2019 年东南亚夏季不同类型极端降水的大尺度环流特征及机理(图略,详见Xu et al., 2022)。结果表明,C1 型极端降水主要与印度洋的冷海温异常高度相关。印度洋的冷海温异常能造成东南亚南部低层风场辐散和下沉运动以及北部低层风场辐合和上升运动,进一步导致东南亚极端降水呈现北湿南干的空间分布。C2 型极端降水主要与热带北大西洋暖海温异常紧密相关。春季热带北大西洋的暖海温异常能够持续到夏季,通过Matsuno-Gill 响应引起热带东北太平洋和中太平洋的北风异常和冷海温异常。冷海温异常进一步通过Matsuno-Gill 响应造成热带西北太平洋的东风异常和上升运动,导致东南亚极端降水呈现一致偏湿的空间分布型。C3 型极端降水主要由热带太平洋和大西洋的海温异常引起。热带西太平洋、东太平洋和大西洋海温异常的暖—冷—暖空间分布,能够造成东南亚南部低层风场辐合和上升运动,北部低层风场辐散和下沉运动,使得东南亚极端降水呈现北干南湿的空间分布。
此外,我们利用SOM 方法提取了1979~2019年东南亚夏季总降水的主模态(图略),并分析了三种类型极端降水和非极端降水对总降水的贡献(图7)。结果表明,1979~2019 年东南亚夏季总降水可以分为全区一致偏干、海洋性大陆主体偏干、北湿南干、全区一致偏湿和北干南湿五种空间分布型。第1、2 类总降水主要由非极端降水决定,第3~5 类总降水分别由C1、C2 和C3 极端降水和非极端降水共同影响。
本文利用观测数据分析了1951~2014 年泛南海地区年平均和季节平均的四个极端降水指数(R95p、R20、RX1day 和RX5day)的时空变化特征。基于高精度统计降尺度数据(NEX-GDDPCMIP6),预估了该区域各个极端降水指数在SSP1-2.6 和SSP5-8.5 情景下2015-2100 年的变化特征。为了探讨不同区域极端降水对气候变化的响应程度,我们将泛南海地区分为四个子区域,分析了各个子区域年平均极端降水在2016~2035 年、2046~2065 年以及2080~2099 年相比于1995~2014 年的百分比变化。此外,东南亚位于热带西太平洋暖池区,是北半球天气气候变化最敏感的地区之一,同时也是亚洲季风区对流最活跃的地区之一。由于受多种海洋和大气系统的影响,东南亚天气气候复杂多变,极端降水频发。以该区域为例,分析了极端降水的大尺度环流特征及变化机理。本文的主要结论如下:(1)1951~2014 年泛南海地区年平均和季节平均的四个极端降水指数的时空变化特征相似。极端强降水和持续性强降水均分布在东南亚、中国东南部以及青藏高原南坡地区,这些区域同时也是高频强降水发生的地区。东南亚各个极端降水指数的季节性差异较小,其他区域四个极端指数在夏季最大,秋季和春季次之,冬季最小。1951~2014 年泛南海地区各个极端降水指数的趋势空间分布不均一,且大部分区域趋势不显著。(2)和历史时期类似,SSP1-2.6 和SSP5-8.5情景下2015~2100 年泛南海地区极端降水大值区主要分布在东南亚、中国东南部以及青藏高原南坡地区,且SSP5-8.5 情景下极端降水的值大于SSP1-2.6 情景。对整个区域而言,2015~2100 年该区域四个极端降水指数均呈显著增加的趋势,且SSP5-8.5 情景下的趋势比SSP1-2.6 情景显著。虽然东南亚是1951~2014 年极端降水的大值区,但该区域极端降水在未来三个时段(2016~2035 年、2046~2065 年和2080~2099 年)的变化没有其他区域明显,南亚和青藏高原是泛南海地区未来极端降水变化的大值区。(3)1979~2019 年东南亚夏季极端降水主要分为北湿南干、全区一致偏湿和北干南湿的空间分布型。北湿南干型主要受印度洋冷海温异常控制,全区一致偏湿型与热带北大西洋暖海温异常密切相关,热带太平洋和大西洋海温异常是北干南湿型极端降水的主要影响因子。
利用SOM 和FDR 方法,本研究提取了东南亚四个季节极端降水空间分布的主模态,但目前仅分析了该区域夏季极端降水主模态的大尺度环流特征及机理,未来需加强其他三个季节极端降水变异机理的研究。例如,春季是季节转换中一个非常重要的季节,春季气候变异过程对夏季气候有着明显的影响(Li et al., 2016; Yang et al., 2019)。在全球变暖背景下,泛南海地区春季气候状况可能呈现新的变化规律,这些变化会加剧我国南方旱涝预测的难度和不确定性,在今后的工作中需要多关注泛南海地区春季极端降水的变化特征及机理研究。
泛南海不同区域极端降水对气候变暖的响应不同,但本研究仅仅关注了东南亚极端降水的变化机理。未来工作需要加强其他区域极端降水的大尺度环流特征及机理的研究。此外,泛南海地区次季节—季节气候预测不仅对当地人民生活至关重要,也对其周边地区的气候预测有重要影响。研究指出,气候模式对泛南海地区环流和平均降水的次季节—季节变化有一定的预报能力(Zhang et al., 2016),但对该区域极端降水可预测性研究存在不足,亟需开展相关研究。