1960—2019年珠江流域多尺度旱涝特征研究.docxVIP

1960—2019年珠江流域多尺度旱涝特征研究 摘要：珠江流域旱涝灾害对生态环境和社会经济造成巨大影响，认识其旱涝时空演变特征对及时制定减灾策略具有重要意义。基于1960—2019年43个气象站点逐月降水资料和标准化降水指数(SPI)，采用ANUSPLIN气象插值方法、Sen-Median趋势分析、Mann-Kendall检验、Hurst指数和经验正交函数分析法(EOF)，研究了珠江流域降水量及多尺度旱涝时空变化特征。结果表明：1)年均降水量总体呈不显著上升趋势( 随着全球气候变暖的加剧，我国干旱问题日益凸显，干旱灾害发生频率高、持续时间长和影响范围广由于干旱成因复杂，易受人类活动影响，常用干旱指数来描述干旱现象。普遍使用的干旱指数有Palmer指数、相对湿润指数、降水距平百分率、标准化降水指数(Standardized Precipitation Index，SPI)、标准化降水蒸散指数(Standardized PrecipitationEvapotranspiration Index,SPEI)，其中SPI计算仅需降水数据，能够较好的反映气象干旱强度和持续时间，适用于不同地区、不同尺度的气象干旱情况分析，因而得到广泛应用上述研究成果对揭示珠江流域旱涝时空演变特征具有重要意义。然而，已有的研究时段相对较短，而且从不同时空尺度研究旱涝演变特征及变化趋势的研究尚不多见。因此，本文基于珠江流域1960—2019年的气象资料，运用最成熟的SPI作为干旱指标，结合ANUSPLIN插值方法、Sen-Median趋势分析、Mann-Kendall检验、Hurst指数及经验正交函数(Empirical Orthogonal Function,EOF)分析法，研究近60年珠江流域年际及季节旱涝时空变化特征及未来趋势，探究珠江流域旱涝的主要空间模态及对应的时间系数变化，以期为流域旱涝灾害的监测及预警提供依据。1 研究区概况珠江流域位于东经102°14′～115°53′、北纬21°31′～26°49′，流经滇、黔、桂、粤、湘、赣等省(区)，流域面积约44万km2 数据及方法2.1 数据来源及处理收集珠江流域1960—2019年43个气象站点逐月降水观测资料，资料来源于中国气象数据共享服务网()，气象站点空间分布如图1所示。为了保证数据的完整性、一致性，数据处理过程中对极少部分缺测数据采用同一站点的相邻月份进行线性插补。基于站点的月降水量，在Excel中统计为年降水量。研究区90 mSRTM(Shuttle Radar Topography Mission) DEM数据来源于地理空间数据云平台(/)，对DEM数据进行镶嵌拼接，然后重采样为500 m，最后利用研究区边界进行裁剪，得到珠江流域地形数据。2.2 研究方法2.2.1 标准化降水指数  SPI最早由McKee等式中：当式中：2.2.2 ANUSPLIN插值法  ANUSPLIN是澳大利亚科学家Hutchinson基于薄盘样条理论编写的针对气象数据曲面拟合的专用插值软件2.2.3 趋势分析方法  Sen-Median趋势分析及Mann-Kendall检验方法已被广泛用于长时间序列气候变化和水文分析的趋势检验2.2.4 Hurst指数  重标极差(Rescaled Range Analysis Method，R/S)分析法由英国水文学家Hurst在研究尼罗河水文问题时提出，后Mandelbrot对该方法进行了理论方面的补充和完善2.2.5 经验正交函数分析法  经验正交函数分析法(EOF)3 结果与分析3.1 降水量时空变化特征3.1.1 时间变化特征  从图2可以看出，近60 a珠江流域年均降水量总体呈不显著上升趋势(3.1.2 空间变化特征  从图3(a)可以看出，珠江流域近60 a年均降水量空间分布具有显著的差异性，总体呈东北高、西南低的空间分布格局。受地形和大气环流的影响，年均降水量空间分布总体呈现随着距离水汽源地越远及地形的抬升而减少的趋势。降水量大于2 000 mm的高值区主要分布在柳江流域、北江流域、桂江流域等区域；降水量小于1 000 mm的低值区集中分布在南盘江流域、右江流域，北盘江流域也有零星分布。基于Sen-Median +Mann-Kendall趋势分析方法，从像元尺度研究近60 a珠江流域年降水量变化趋势的空间分布格局。其中，3.2 SPI时空变化特征3.2.1 时间变化特征  基于珠江流域气象站点的年际SPI，计算1960—2019年年均SPI并绘制图4曲线。从图4可以看出，1960—2019年珠江流域年际S