第一作者:赵新来(1992-),男,安徽阜阳人,在读硕士生,研究方向为草地遥感与地理信息系统。E-mail:zhaoxinlai@163.com
利用甘南州及周边地区33个气象台站1965-2013年的降水和温度数据,比较分析了降水距平指数(Pa)、标准化降水指数(SPI)和标准化降水蒸散指数(SPEI)对高寒草地干旱响应的有效性,并结合历史旱灾统计资料进行了结果验证。结果表明,从Pa指数、SPI指数、SPEI指数与草地NDVI的相关分析可以看出,SPI6和SPEI6与草地NDVI关系较好,Pa、SPI1、SPEI1、SPI3、SPEI3不足是波动不稳定,SPI12、SPEI12对干旱的反应表现出一定的滞后性。SPI6和SPEI6均适于定量研究高寒草地的干旱状况,结合历史统计资料验证,两项干旱指数监测效果相对较好。本研究结果将为有关青藏高原高寒草地的干旱研究提供有益的参考。
Using precipitation and temperature data from 33 meteorological stations around Gannan in 1965-2013, the percentage of precipitation anomalies (Pa), standardized precipitation index (SPI), and standardized precipitation evapotranspiration index (SPEI) were calculated. To select a suitable index to describe dry climatic factors in alpine meadows, we analysed the response of these indexes to drought. In addition, verification was conducted by comparing historical disaster statistics with the drought indexes. The correlation analysis of SPI, SPEI, and Pa with grassland normalized difference vegetation index (NDVI) showed that the reaction of Pa, SPI1, SPEI1, SPI3, and SPEI3 to precipitation was too sensitive; however, SPI12 and SPEI12 were slightly hysteretic. SPI6 and SPEI6 were both well suited to quantitatively evaluating the dry climatic factors in this area. Based on historical statistics, two drought indexes were somewhat better. These results provide effective data in support of further research on drought risk analysis in alpine grasslands.
在全球气候变暖的大背景下, 全球范围内特大干旱、高温等极端天气气候事件发生的频率和强度都呈现增加趋势[1, 2]。其中, 由天然降水异常引起的干旱是人类必须面对的重大环境问题[3]。干旱在气象学上有干旱气候和干旱灾害两种理解。干旱气候是指蒸发量比降水量大得多的一种气候; 干旱灾害是指某一地理范围在某一具体时段内的降水量比多年平均降水量显著减少, 导致该地区的经济活动和人类活动受到较大危害的现象。1997年, 美国气象学会(American Meteorological Society, AMS)将干旱分为4类:气象干旱或气候干旱、农业干旱、水文干旱和社会经济干旱, 这受到了众多研究干旱的科学工作者的普遍认可[4]。之后, 有研究学者又提出了生态干旱, 并认为气象干旱主要反映了降水强度及其概率特征, 而农业干旱、生态干旱、水文干旱和社会经济干旱则分别表征了气象干旱对农业、生态、水资源和社会经济活动的影响程度[5]。
青藏高原东部处于温带大陆性气候、温带季风气候和高原山地气候交汇处, 常年降水偏少, 尤其近50年来降水呈现非周期振荡, 干旱事件频繁发生[6, 7]。极端的气候变化对该区域草地植被的生长也产生了较为明显的影响[8, 9, 10, 11]。2015年11月27日联合国粮食与农业组织(Food and Agriculture Organization of the United Nations, FAO)召开亚洲草地可持续管理专家咨询会议(Regional Consultation on Sustainable Grassland and Pasture Management in Asia), 就干旱区的草地恢复、气候变化以及草地保护等问题进行讨论, 指出干旱区草地面临严峻的生态环境问题。相对于森林等其它生态系统, 草地生态系统较为脆弱, 植被类型以草地为主, 气候的变化对依赖降水涵养的天然草地的影响尤为明显, 同时对草原上的主要经济产业的影响更为明显[12, 13, 14]。目前, 干旱领域的研究主要集中在以平原或山区种植业为主的农业干旱, 统计资料和研究方法已经较为完善, 对草地, 尤其是高寒草地的干旱研究鲜见报道。因此, 选择适宜的干旱指数监测草地干旱成为研究的关键之一。
干旱指数是表征干旱程度的标准, 也是旱情描述的数值表达, 但由于干旱地面监测既要考虑气象因素, 也要考虑水文、土壤、作物、水利灌溉等诸多因素, 因而提出精确而普遍使用的干旱指数有很大困难。为了监测干旱, 研究者先后提出了帕默尔干旱指数(palmer drought severity index, PDSI)、地表水分供应指数(surface water supply index, SWSI) 和标准化降水指数(standardized precipitation index, SPI) 等50多种干旱指数[15]。其中SPI由于具有计算简单、时间尺度灵活和对干旱时空变化指示效果较好等优点而被广泛应用[16, 17, 18]。近年来Vicente-Serrano等[19]提出的标准化降水蒸散指数(standardized precipitation evapotranspiration index, SPEI)考虑了降水与蒸散, 也保留了SPI、PDSI指数对干旱事件的敏感度, 受到广泛的关注。因此, SPI和SPEI常被学者用来研究由气候因素变化所导致的干旱等极端气候现象, 可以较好地反映干旱强度与持续时间, 而且能反映不同时间尺度和区域的变化情况。利用降雨变异指数(Z指数)、帕默尔干旱指数和不同时间尺度的SPI(SPI1、SPI2、SPI3、SPI6、SPI9、SPI12和SPI24)作为干旱的定量度量, 将植被健康指数(vegetation condition index, VCI)作为植被对干旱的响应, 分析它们的相关关系, 最后得出SPI6与VCI的相关性最大(R2=0.54)[20]。利用SPI分析了美国大平原北部地区草地和作物对干旱的响应, 探究了不同时间尺度的SPI(SPI1、SPI2、SPI3、SPI6、SPI9和SPI12)与生长季草地或农作物(5-10月)的NDVI线性相关关系, 结果显示, 草地中SPI3和SPI6的相关性较高, 在9月份相关系数最高(R=0.774)[21]。近年来, 国内关于草地生态系统对干旱现象响应的研究逐渐兴起。有学者利用内蒙古43个气象台站2004-2011年5月-8月降水量和降水距平百分率观测资料结合调查统计资料对内蒙古草地干旱损失做了评价[22]。以内蒙古草地为研究区域研究了近50年来荒漠草原气候干燥度的空间分布及其演变特征[23]。以内蒙古草地为研究区域, 利用SPI6干旱指数探究了不同草地类型草地净初级生产力对干旱现象的响应[24]。利用SPEI指数对青藏高原高寒草地自然保护区的气候年际变化作了研究[25]。关于草地对干旱响应的研究日趋增多, 然而研究的草地类型大多是典型草原, 关于高寒草地对干旱的响应还相对少见。为此, 本研究利用地面气象数据、遥感数据和历史统计资料, 集中探究干旱指数在高寒草地干旱监测中的适宜性。
甘南藏族自治州地处青藏高原与黄土高原过渡地带(33° 06'-35° 34' N, 100° 45'-104° 45' E), 位于东亚季风区、西北干旱区和青藏高原高寒区的交汇处, 总面积约4.5万km2, 高寒草地面积约2.75万km2, 草地类型主要有高寒草甸和山地草甸, 丰富的草地资源不仅对国民经济可持续发展有举足轻重的作用, 对黄河、长江上游甚至中下游的生态环境与经济社会的可持续发展也具有重要作用。
本研究中的数据来源主要有气象数据、遥感卫星数据、历史统计资料以及其它矢量数据等。其中, 气象数据是甘南州及其周边33个气象台站1965-2013年逐月降水和温度资料, 数据来源于中国气象数据网(http://data.cma.cn/)。遥感数据是基于Terra 和Aqua 两颗星观测所得16天合成MOD13Q1 L3级产品数据, 覆盖研究区的影像轨道号为H26V05, 空间分辨率为250 m, 数据来自NASA数据应用网站(http://glovis.usgs.gov/), 文中所用到的数据为2000年5月-9月时段数据。研究中用到的年鉴统计数据主要是《中国气象灾害大典· 甘肃卷》、《中国暴雨洪水及干旱年鉴》等[26, 27]。除此之外, 研究用到的数据还有甘南地区1:50万草地类型图、1:400万行政分区图等空间矢量数据(图1)。
1.3.1 Pa指数 降水距平百分率可以用来反映某时段降水与周期平均状态的偏离程度, 能直接反映降水异常引起的干旱。
Pa=
式中:P为某时段降水量,
![]() | 表1 Pa、SPI、SPEI干旱分级标准[28, 33, 34, 35] Table 1 Classification standard of drought based on Pa、SPI and SPEI[28, 33, 34, 35] |
1.3.2 SPI指数 SPI指数是运用数学方法将降水量的累积频率分布转化为标准正态分布, 具有无量纲、标准化特点, 并适用于不同地区[29, 30, 31]。
设某一时间尺度的降水量为x, 假设降水量序列服从Gamma分布, 则其概率密度函数为:
p(x)=
式中:β 和α 分别是尺度参数和形状参数, 均大于零, 用极大似然估计方法求得, 函数表达式见式(2)和式(3)。
α =
β =
式中:A=ln(
令q=
H(x)=q+(1-q)G(x) (4)
将累积概率分布H(x)转换到标准正态分布, 即得到对应的SPI值:
当0< H(x)≤ 0.5时, 令k=
SPI=-
当0.5< H(x)< 1时, 令k=
SPI=k-
式中:c0=2.515 517, c1=0.802 853, c2=0.010 328; d1=1.432 788, d2=0.189 269, d3=0.001 308。
因为SPI是通过计算给定时间尺度降水的累积概率来计算干旱指数的, 所以利用当前月和前n个月(n=1, 2, 3, …)的降水量来计算时间尺度为(n+1)个月的SPI, 简称SPIn+1。Guttman[32]认为时间尺度为1-24个月最具有实际应用的价值, 并认为1-2个月时间尺度的SPI适用于气象干旱, 1-6个月的时间尺度的SPI适用于农业干旱, 大于6个月时间尺度的SPI较适用于水文干旱的分析应用。SPI指数优点是计算稳定性良好, 时间尺度灵活, 对时空变化的指示效果良好, 能有效地反映旱涝状况[33]。
1.3.3 SPEI指数 SPEI是Vicente-Serrano在SPI的基础上引入潜在蒸散项构建的, 融合了SPI和PDSI的优点, 计算公式如下:
SPEI=W-
式中:W为蒸散降水推导函数的累积概率函数值, 详细步骤可见参考文献[19]。与SPI一样, SPEI具有空间一致性、多时间尺度且计算简单的特点, 同时还综合考虑了降水和温度因子, 适用于气候变暖背景下干旱检测与评估[34, 35]。
1.3.4 归一化指数 归一化植被指数(normalized difference vegetation index, NDVI)是根据植被对光谱反射曲线的特点发展而来的比值型植被指数, 可使植被从水体和土壤中分离出来, 是植被生长状态及植被覆盖度的有效指示因子[36]。在区域干旱程度增加时, 植被生长会因水分不足而受到影响, 其表现在遥感影像上的特征是NDVI相应地减小, 相反地, 当干旱消退时, 植被NDVI也相应地增大。因此, 植被对干旱的响应可以表征在植被NDVI上。
NDVI=
式中:NIR和R分别是近红外波段和红光波段的地面反射率。
根据1:50万的草地类型图, 剔除非草地区域的18个台站点, 选取剩下的位于高寒草地区域的15个其它台站点作为样本。该15个台站点的草地类型相同, 均为高寒草地, 海拔地形相近, 最大程度上消除了其它因素的影响。图2、图3和图4展示了2000年草地生长季(5月-9月)NDVI和Pa、不同尺度SPI、SPEI(1, 3, 6, 9个月)的线性相关关系及显著性检验结果。用到的相关系数是Pearson简单相关系数, 显著性检验采用0.01显著性水平下的双尾检验。在5、8月份, Pa指数和高寒草地NDVI显著相关, 且均通过显著性检验(P< 0.01), 而在7和9月Pa和高寒草地NDVI的相关关系不显著(P> 0.01)(图2)。5月份, SPI1、SPI3和 SPI6与高寒草地NDVI具有显著的相关性, 6月份, SPI3、SPI6和高寒草地NDVI具有显著的相关性; 7月份, SPI3、SPI6、SPI12和高寒草地NDVI具有显著的相关性; 8月份, 4个时间尺度的SPI指数和NDVI均有显著的相关性; 9月份, SPI6、SPI12和高寒草地NDVI具有显著的相关关系(图3)。5月份, SPEI1、SPEI3和SPEI6与高寒草地NDVI具有显著的相关性; 6月份, SPEI1、SPEI3、SPEI6和高寒草地NDVI具有显著相关性; 7月份, SPEI3、SPEI6和高寒草地NDVI具有显著相关性; 8月份, 4个时间尺度的SPI指数和NDVI均有显著相关性; 9月份, SPEI12和高寒草地NDVI具有显著相关性。综合比较, 在5月-9月, 在不同尺度干旱指数中, SPI6和SPEI6与高寒草地NDVI均有相对显著的相关关系。进而表明, 当气候发生变化时, 植被生长状况也相应的变化, 即表现为植被NDVI的变化, 可以用干旱指数和NDVI的相关性高低来衡量干旱指数的适宜性。因此, 从结果可以看出, SPI6和SPEI6两种干旱指数更适宜衡量该地区的干旱程度。
绘制出近50年的时间序列图(图5), 参考表1的分级标准, 同时查阅《中国气象灾害大典》等历史灾害记载资料[26, 27], 列出近年50年甘南地区旱灾统计与干旱指数检测结果(表2)。查询到的历史记录共37条, 其中SPI6检测结果和历史记录相符合的有27条, 检验精度为72.97%; SEPI6检测结果和历史记录相符合的有26 条, 检验精度为70.27%。将研究结果和此检验结果综合分析, 可以得出, SPI和SPEI6干旱指数均适宜衡量该地区的干旱程度, 且检验精度都在70%以上。
![]() | 表2 近50年甘南地区旱灾统计资料和干旱指数检测结果比较 Table 2 Historical statistic of drought and corresponding drought index |
借助气象数据, 以农业干旱风险分析为基础, 对高寒草地干旱气候响应下的干旱指数进行了比较分析。结果发现, Pa指数在评估草地生长季7月和9月时, 与草地NDVI的相关性不明显。在7月份大多样点的Pa< 0, 表明7月份降水相对不足, 而NDVI并未相应地明显减小。7月份正是高寒草地的盛草期, 可能是因为此时相对偏少的降水并不是主要限制因素, 而温度正适合草地生长。在9月份, 过半样点的Pa> 0, 表明9月份降水相对充足, 而此时Pa与NDVI的相关关系并不明显, 可能是因为9月份气温降低, 此时, 温度是限制草地生长的主要因子。
比较分析不同尺度的SPI和SPEI, 其中SPI1、SPEI1与草地NDVI在7月和9月相关性不显著; SPI3、SPEI3与草地NDVI在9月的相关性不显著; SPI12与草地NDVI在5月和6月的相关性不显著; SPEI12与草地NDVI在5月、6月、7月的相关性不显著; SPI6、SPEI6与草地NDVI在5月-9月的显著相关。SPI1和SPEI1与草地NDVI的关系和Pa指数的结果相似, 由此可以看出, Pa、SPI1和SPEI1的不足是波动不稳定。SPI3和SPEI3与草地NDVI的相关关系5月-8月显著, 在9月份不显著, 但略好于Pa与9月份草地NDVI的相关性。与Pa相比, SPI6和SPEI6在5月-9月与草地NDVI的相关性显著且较为稳定。由SPI12和SPEI12与草地NDVI的相关性的结果可以看出, 它们具有明显延后的缺点。
甘南地区的旱灾资料记录有37条, 和SPI6、SPEI6指数的监测结果进行比较验证, 两种指数的精度分别是72.97%和70.27%, 进一步证实了SPI6和SPEI6在甘南地区的监测效果较好。而本研究分析中只考虑降水因素的SPI6指数与考虑降水和温度的SPEI6指数和NDVI的相关性结果相似。多数研究认为, 气温和降水是影响植被周期性NDVI变化的重要因素, 但也有研究表明干旱和半干旱地区降水因子对植被NDVI的影响作用明显大于温度因子[37, 38, 39, 40]。甘南藏族自治州位于青藏高原东部, 气温常年相对较低, 蒸散量较小[41], 因此才会得出SPI6和SPEI6的监测效果相近。
通过草地NDVI对Pa指数、SPI指数、SPEI指数的响应分析, SPI6和SPEI6与草地NDVI显著相关, Pa、SPI1、SPEI1、SPI3、SPEI3的不足是波动不稳定, SPI12、SPEI12对干旱的反应表现出一定的滞后性。SPI6和SPEI6均适于定量研究高寒草地的干旱状况, 结合历史统计资料验证, 两项指数监测效果相对较好。以上研究结果将为日后我国高寒牧区的干旱生态风险评价工作提供一定的数据支持和学术参考。
The authors have declared that no competing interests exist.
[1] |
|
[2] |
|
[3] |
|
[4] |
|
[5] |
|
[6] |
|
[7] |
|
[8] |
|
[9] |
|
[10] |
|
[11] |
|
[12] |
|
[13] |
|
[14] |
|
[15] |
|
[16] |
|
[17] |
|
[18] |
|
[19] |
|
[20] |
|
[21] |
|
[22] |
|
[23] |
|
[24] |
|
[25] |
|
[26] |
|
[27] |
|
[28] |
|
[29] |
|
[30] |
|
[31] |
|
[32] |
|
[33] |
|
[34] |
|
[35] |
|
[36] |
|
[37] |
|
[38] |
|
[39] |
|
[40] |
|
[41] |
|