川滇农牧交错带极端降水的时空演变及突变分析

周 丽，张 彬，吴 彬

1. 内江师范学院 地理与资源科学学院，四川 内江 641112；2. 内江市东兴区气象局，四川 内江 641100

在世界许多地区，极端降水是最重要的气候变量之一[1]，降水对自然和社会环境有着实质性的影响． 政府间气候变化专门委员会(IPCC)认为，由于气候变化和人类活动，极端降水事件的数量和强度将继续增加[2-3]． 例如，2000-2019年，国际灾害数据库(EM-DAT：The Emergency Events Database)报告全球共发生3 174起洪水灾害，远远超过1960-1999年报告的1 807起洪水灾害． 显然，在地形复杂的地区，极端降水强度的增加会导致更多的地貌灾害，如山洪、水土流失、山体滑坡、泥石流等． 因此，对极端降水事件的空间分布和趋势进行评价是十分重要和迫切的． 极端降水事件引起了科学界越来越多的关注． 在过去几十年中，极端降水事件的增加趋势已在气候较温暖的大多数地区得到了检测[4]． 然而，全球各地极端降水的空间格局有很大的差异[5]． 以往研究表明，极端降水对气候变暖具有高敏感性，尤其是在山地和高原地区[6]，这种趋势在极端降水事件较为常见的中国西南地区最为明显[7]．

川滇农牧交错带位于中国西南部，东经98°08′-104°58′，北纬26°28′-34°19′，包括川西36个县和滇西北7个县[8]，除少数县的部分地区外，43个县大部分属于青藏高原． 川滇农牧交错带地形复杂，以高原和高山峡谷为主要地貌类型，包括3个地貌单元：松潘高原、四川西部高原和横断山脉(图1)． 南部的横断山脉表现为南北向高山与纵谷交错分布[9]，大多数河流由北向南流动． 横断山脉阻挡了东西向的东亚季风，为南北向的南亚季风提供了通道[10]． 川西高原和松潘高原分别位于中部和北部地区． 川滇农牧交错带是青藏高原到云贵高原和四川盆地的主要过渡地带，是中国滑坡、山洪、泥石流的主要发生地区，极端降水在这些地貌灾害中起着重要作用[11-12]．

目前针对西南地区极端降水的时空变化特点，国内外学者做了许多工作[13-15]，研究区包括四川省[16-18]、云南省[19]． 但以上研究区域范围较大，对极端强降水反应敏感的川滇农牧交错带研究甚少，而对中国北方农牧交错带的研究较多[20-21]． 因此，本文拟研究近45年来川滇农牧交错带极端降水的时空演变及趋势，对理解研究区与极端降水相关的灾害和制定适当的适应和缓解战略至关重要．

审图号：GS(2016)2884号，下同．图1 川滇农牧交错带各气象站点分布

1 数据来源与方法

1.1 数据来源

川滇农牧交错带DEM数据来源于地理空间数据云 (http：//www.gscloud.cn/)，气候数据是由中国气象科学数据共享服务网(http：//cdc.cma.gov.cn)提供的逐日降水数据序列． 选取1971-2015年川滇农牧交错带20个分布均匀、无缺测情况的站点(图1)． 采用气候变化检测和指数专家小组定义的极端降水指数中的4个指标(表1)，对川滇农牧交错带极端降水事件近45年来的时空变化特征进行分析．

表1 本研究所采用的极端降水指数

1.2 研究方法

采用泰森多边形法[22]和一元线性回归方程计算各极端降水系列的气候倾向率[23]，并利用Mann-Kendall[24]和滑动T检验[25]对极端降水指数时间序列的突变特性进行分析；运用反距离权重空间插值法(IDW)对川滇农牧交错带极端降水进行空间分析．

2 结果与分析

2.1 极端降水的时间变化特征

通过泰森多边形法赋予研究区各个站点不同的权重系数，分别计算川滇农牧交错带全区域以及松潘高原、川西高原和横断山区的极端降水指数． 图2-3分别显示了4个极端降水指数的年际变化序列．

2.1.1 强降水量和强降水日数

由图2a可以发现，川滇农牧交错带全区域R95pTOT有微弱上升趋势，经检验达到0.1显著性水平． 松潘高原、川西高原和横断山区R95pTOT均有较弱的上升趋势，其中川西高原上升趋势相对最大，经检验达到0.1显著性水平，说明川滇农牧交错带强降水量在增加，更易发生山洪、泥石流等灾害． 全区域R95pTOT多年均值达到144.56 mm，松潘高原、川西高原和横断山区R95pTOT多年均值分别为143.69，130.27，182.45 mm，横断山区R95pTOT比松潘高原和川西高原高27.00%和40.06%． 因此，大于95%分位数的强降水占该区域降水的比例在19.23%～21.69%，其中横断山区和松潘高原在年降水总量呈下降趋势的前提下，强降水量仍表现出上升趋势． 由图2b可见，川滇农牧交错带全区域(经检验达到0.1显著性水平)、松潘高原、川西高原以及横断山区的R95pD均表现出波动上升趋势． 其中，松潘高原R95pD多年均值和气候倾向率均为最大． 横断山区R95pD多年均值和气候倾向率均为最小，同时，横断山区R95pTOT为最大，揭示出横断山区是研究区强降水量最大且在时间上最集中的区域，是川滇农牧交错带最易发生山洪、泥石流等灾害的区域．

图2 强降水量(a)和强降水日数(b)的年际变化

2.1.2 一日最大降水量和五日最大降水量

由图3可见，川滇农牧交错带全区域、松潘高原、川西高原和横断山区RX1day均呈上升趋势，且川西高原上升幅度最大，经检验达到0.05显著性水平；同时，横断山区RX1day多年均值最大，高于川西高原43.25%，加之气候倾向率为正，导致该区域一日最大降水量仍为全区域最高． 全区域、松潘高原和川西高原RX5day均呈波动上升趋势，且川西高原上升幅度最大，而横断山区则呈波动下降趋势． 因此川滇农牧交错带短历时的极端降水呈上升趋势，且横断山区在时间上的集中性比其他区域表现得更为明显．

图3 一日最大降水量(a)和五日最大降水量(b) 的年际变化

2.2 极端降水的空间变化特征

运用反距离权重空间插值法将各站点极端降水指数的气候倾向率进行可视化，如图4-5所示．

2.2.1 强降水量和强降水日数

总体来看R95pTOT变化的空间分布(图4a)，川滇农牧交错带中东部升高，西北、西南两侧降低；川西高原表现为以上升趋势为主，高值区位于川西高原东部的小金和康定附近；松潘高原仍以上升趋势为主，横断山区西南部和东南部主要呈下降趋势． R95pD与R95pTOT有较相似的空间分布(图4b)，上升的高值区分布在川西高原东部的小金和康定以及横断山区东部的木里，经检验均达到0.05或0.1显著性水平． 下降大值区分布在横断山区西部的贡山． 川滇农牧交错带东部地区R95pTOT与R95pD均呈现显著上升趋势，说明该地区受极端降水的影响在加大．

图4 川滇农牧交错带强降水量(a)和强降水日数(b)变化的空间分布

2.2.2 一日最大降水量和五日最大降水量

由图5可见，RX5day较RX1day的下降区有明显扩大，上升区则有明显缩小． 全区域大部分站点的RX1day呈上升趋势，松潘高原和川西高原呈上升趋势的站点分别占各区域的66.67%和72.73%，仅横断山区以下降趋势为主，约66.67%的站点表现为下降趋势． RX5day 空间变化差异明显，松潘高原和横断山区以下降趋势为主，约占各所在区域站点总数的66.67%和83.33%，而川西高原主要表现为上升趋势，约占川西高原站点总数的81.82%． 整体来说，RX1day，RX5day 以上升趋势为主，表明川滇农牧交错带极端强降水的量和强度均呈上升趋势，且极端强降水事件的发生在较短历时的情况下增加更为显著．

图5 川滇农牧交错带一日最大降水量(a)和五日最大降水量(b)变化的空间分布

2.3 极端降水的突变分析

极端降水的显著变化趋势对川滇农牧交错带洪灾风险有可能造成重要影响，因此，通过Mann-Kendall检验，选取有显著变化趋势的木里、理塘两个典型站点，对其极端降水情况做突变检验． 图6-7为近45年川滇农牧交错带4个极端降水指标典型站点的M-K检验，UF和UB分别为按顺序和逆序时间序列计算的统计量序列[13]． 若UF(UB)大于0，表示序列呈上升趋势，反之则呈下降趋势． 图中红色虚线为临界线，临界值为±1.96，即达到0.05的显著性水平，绝对值大于1.96的UF(UB)表示下降或上升趋势显著． 当UF和UB在临界线之间出现唯一交点，则交点对应的时刻为突变开始时间；当UF和UB出现多个交点，利用滑动T检验法做进一步显著性检验，能通过检验的交点对应的年份为序列发生突变的年份[26]；当UF或UB超过临界线，超过临界线的范围则为出现突变的时间区域．

2.3.1 强降水量和强降水日数

由图6a可知，UF与UB在1984-1985年和1991年出现3次交点，均落在临界线内，且UF值大于0，说明木里站R95pTOT可能的突变发生在1984-1985和1991年左右，且R95pTOT呈上升趋势，因此有必要对突变年份进行进一步确定，对UF和UB曲线的3个交点进行滑动T检验，仅1984年的交点通过显著性检验，说明1984年木里站R95pTOT出现了突然上升的过程． 图6b中，UF与UB在1976年存在1个交点，且交点落在临界线之间，说明木里站R95pD在1976年前后发生了1次突变，且在1976年之后UF值大于0，表明在1976-2015年有持续增加趋势，1984年开始增加趋势通过0.05的显著性水平检验，说明木里站R95pD在1976年发生了1次明显的突变后呈增长趋势，并在1984年后呈明显的持续增长． 从图6可知，木里站R95pTOT突变比 R95pD 突变发生晚，且R95pD上升趋势较R95pTOT显著．

图6 川滇农牧交错带典型站点强降水量(a)和强降水日数(b)的Mann-Kendall检验

2.3.2 一日最大降水量和五日最大降水量

由图7a可知，木里站RX1day于1971-1986年期间处于上升—下降持续波动状态，1986年开始出现稳定上升(UF>0)，且在1997年先后出现显著上升． 理塘站RX5day(图7(b))在1979年发生突变后，持续上升(UF>0)，并于1989年出现显著上升．

图7 川滇农牧交错带典型站点一日最大降水量(a)和五日最大降水量(b)的Mann-Kendall检验

从川滇农牧交错带两个典型站点4个极端降水指标的M-K突变检测可知，木里站和理塘站发生突变的时间较接近，整体来看，突变的时间大致发生在20世纪70-80年代．

3 结论与讨论

(1) 川滇农牧交错带全区域的平均极端降水指数中，R95pTOT，R95pD，RX1day，RX5day均呈上升趋势；松潘高原R95pD为各区域最高，而横断山区除R95pD外，其他各极端降水指数均大于松潘高原和川西高原；同时，川西高原各指数年际变化比松潘高原和横断山区更剧烈．

(2) R95pTOT与R95pD的空间分布相似，上升的大值区均集中在川西高原东部和横断山区东部，下降大值区均分布在横断山区西南部． 同时，RX1day，RX5day以上升趋势为主，且随着持续时间的延长，极端降水量的上升区有明显的缩小，下降区有明显的扩大，说明极端强降水事件的发生在较短历时的情况下增加更为显著．

(3) 近45年来川滇农牧交错带典型站点极端降水指数突变的时间大致发生在20世纪70-80年代． R95pD，RX5day发生突变的时间最早，分别发生在1976年和1979年左右；其次为R95pTOT，RX1day，突变依次发生1984年和1986年． 需要强调的是，以上4个指标均发生突变上升趋势． 同时，两个典型站点发生突变的时间较接近．

研究过程中发现，川滇农牧交错带极端降水强度有增大的趋势，且较短历时强降水事件增加更明显． 因此，川滇农牧交错带在未来极端降水指数保持上升趋势的情况下，对显著变化的站点须做深入研究，将极端降水强度指标与坡度、坡向、植被覆盖等因子相结合，综合分析该区域地貌灾害的致灾因子及等级． 另外，该区域由于地形复杂，降水具有典型的局地性特征，应在强降水事件高发区域增加雨量观测站的空间密度，使其观测数据更具代表性，以便更好地因地制宜实施防灾减灾措施，为当地农业生产服务．