干旱是对人类社会影响最严重的气候灾害之一[1], 它具有出现频率高、持续时间长、波及范围广的特点。干旱的频繁发生和长期持续不但会给国民经济特别是农业生产等带来巨大的损失, 还会造成水资源短缺、荒漠化加剧、沙尘暴频发等诸多不利影响。因此监测、研究、预测干旱变化具有重要的意义。
Dai等[2]利用1870-2002年Palmer指数的计算结果分析表明, 20世纪70年代后期以来, 全球极端干旱面积增加了一倍多, 这种干旱趋势的加剧主要是由于干旱化地区降水减少和气温升高造成的。在我国, 干旱灾害发生频繁, 1949-2006年平均每年受旱面积2122×104 hm2, 约占各种气象灾害受灾面积的60%, 加上近几十年来经济迅速发展、人口增长等原因, 干旱给社会造成的不利影响和对人类生存环境的危害日趋严重。我国北方地区继1997年发生了大范围的干旱后, 1999-2002年又连续4年少雨干旱。频繁发生的干旱以及干旱对社会和环境产生复杂和深远的影响[3]。2006年夏季, 重庆发生百年一遇特大伏旱, 四川遭遇近50年最重干旱[4]。这些严重的干旱事件及其引发的不利影响引起了政府决策部门和公众的广泛关注。特别是对于一些气候脆弱区, 如果干旱化趋势进一步加重, 不仅使可利用的水资源减少, 人类无法继续生存, 而且荒漠化和沙漠化将导致植被的消失, 对周边环境也会产生更加恶劣的负面作用[5]。不少学者对我国的干旱变化进行了大量研究。王志伟等[6]根据1950-2000年我国629个站逐月降水资料, 采用Z指数作为旱涝等级划分标准计算了干旱发生的范围, 指出51年来我国北方主要农业区干旱面积呈扩大趋势, 特别是华北等地干旱面积扩大迅速, 形势严峻, 而且极端年份干旱面积显著扩大。马柱国等[7-8]利用地表湿润指数分析了1951-2000年北方地区极端干湿事件的演变规律, 指出近10年来东北和华北地区极端干旱频率显著增加, 而极端湿润发生的频率相对减少; 与气温变化趋势的比较还发现, 极端干旱的频发区往往对应着增温明显的地区, 而西北地区极端干旱发生频率出现了下降的趋势。王绍武等[9]、施雅风等[10]的研究表明我国西北地区有变湿的趋势。
以上学者虽然都在不同程度上对我国干旱变化特征做了有意义的研究, 但大部分研究都集中在局部地区或者研究中缺少最近几年最新的干旱演变情况分析。本文利用《气象干旱等级》国家标准[11]推荐使用的综合气象干旱指数, 详细分析了1951-2006年全国及各区域干旱的变化特征。
1 资料和方法 1.1 资料本研究中所用到的资料包括我国地面606个台站1951-2006年的逐日降水量和平均气温资料 (图 1), 大部分站点属于国家基准站和基本站, 个别为一般气象站。资料由国家气象信息中心提供, 均经过严格质量检验。台站变迁、仪器变化、观测方法改变等是影响气候资料均一性的重要因素。从我国基本、基准气象台站的历史沿革来看, 约有一半的台站均有过迁移[12]; 在仪器和观测方法上, 气温、降水的观测仪器和观测规范也发生过多次变化, 但多发生在20世纪60年代中期之前[12-13]。为了减少资料的非均一性, 本文在资料选取时规定:月内逐日资料缺测达5 d的, 该月资料定为缺测; 1年中有1个月缺测时, 该年年值为缺测; 资料时间序列应包括气候平均场 (1971-2000年), 时间长度至少为40年。原始资料库中共包含731个台站, 本文剔除了不满足上述条件的125个站点资料, 最终选取了606个站点 (图 1), 剔除站点数占总站点数的17%。选取的606个地面观测台站覆盖了我国大陆绝大部分地区。本文选取的资料集中, 20世纪50年代初期站点数量较少, 仅有150~350个左右, 而且空间分布不均, 西部地区除了新疆有少量台站外, 其他地区记录很少; 50年代中期以后收录的观测站点数迅速增加, 西部的站点也逐渐增多, 1956年以后我国西部所有省区都有了观测站; 60年代以后资料站点数均在600个以上, 各年的站点数量变化不大。经过上述质量控制后的资料集时间序列长, 覆盖率较高, 具有较好的代表性和可靠性。为了分析不同区域的干旱变化特点, 本文分别讨论了东北、华北、西北东部、西北西部、长江中下游、华南、西南和青藏高原8个地区近半个世纪干旱的变化特征 (图 1)。在讨论各区域的干旱长期变化时, 采用了格点化计算面积的方法, 在一定程度上有助于增加时间序列的可靠性[14]。
1.2 方法
本文中采用《气象干旱等级》国家标准GB/T 20481-2006中推荐使用的综合气象干旱指数来统计分析近50多年来我国的干旱时空分布特征[11]。综合气象干旱指数是利用近30 d (相当于月尺度) 和近90 d (相当于季尺度) 降水量标准化降水指数, 以及近30 d相对湿润指数进行综合而得, 该指标既反映短时间尺度 (月) 和长时间尺度 (季) 降水量气候异常情况, 又反映短时间尺度 (影响农作物) 水分亏欠情况。该指标适合实时气象干旱监测和历史同期气象干旱评估。综合气象干旱指数IC的计算见式(1)。
(1) |
式(1)中, Z30, Z90分别为近30 d和近90 d标准化降水指数SPI值, 计算方法见《气象干旱等级》国家标准GB/T 20481-2006中附录C[11]; M30为近30 d相对湿润度指数, 由式(2) 计算; a为近30 d标准化降水系数, 由达轻旱以上级别Z30的平均值除以历史出现的最小Z30值得到, 平均取0.4; b为近90 d标准化降水系数, 由达轻旱以上级别Z90的平均值除以历史出现最小Z90值得到, 平均取0.4; c是近30 d相对湿润系数, 由达轻旱以上级别M30的平均值, 除以历史出现最小M30值得到, 平均取0.8。
相对湿润度指数的计算见式(2):
(2) |
式(2)中, P是某时段的降水量; PE为某时段的可能蒸散量, 本文中用Thornthwaite方法计算得到。Thornthwaite方法求算可能蒸散量是以月平均温度为主要依据, 并考虑纬度因子 (日照长度) 建立的经验公式, 需要输入的因子少, 计算方法简单:
(3) |
式(3)中, Ti是月平均气温, 单位: ℃; H是年热量指数; A为常数。年热量指数H和常数A分别由相应的经验公式计算, 具体方法见文献[11]。
当月平均气温Ti≤0 ℃时, 月热量指数Hi=0, 可能蒸散量PE=0(单位: mm/月)。
另外, 本文采用估算干旱面积的方法来分析近50多年我国的干旱变化特征, 具体计算方法如下。
在计算干旱影响范围时, 先将干旱指数值的站点资料格点化, 转化为对应的2°×2°格点资料, 计算干旱影响范围时, 按照IC≤-0.6为出现干旱的区域。每个IC≤-0.6格点覆盖的面积近似为
(4) |
式(4)中, xgrid和ygrid分别为经度和纬度格距, 110.0×111.0 km2为赤道地区1°×1°网格的近似面积, cosΨi是对纬度的订正, 则第j年n个格点的总面积Sj为
(5) |
对发生过干旱的格点按式(4)、(5) 统计面积, 以此表示干旱影响范围的大小。考虑到干旱资料为地面站点资料, 在沙漠和青藏高原高寒地区站点分布稀少, 会造成部分格点资料无有效值。为了减少因此带来的误差, 将估算的格点面积又转化为其占有效格点面积的百分率。计算第j年有效观测区格点的覆盖面积为SNj, IC≤-0.6的出现干旱的格点的覆盖面积为Sj, 则第j年的干旱面积百分率SSj为
(6) |
此结果即最终的干旱面积百分率。
2 全国及各区域干旱的长期变化特征 2.1 全国干旱的长期变化图 2显示了基于IC指数统计的全国干旱面积百分率历年变化。从年代际变化看, 在近半个多世纪中, 我国干旱较重的时期主要出现在20世纪60年代、70年代后期至80年代前期、80年代中后期以及90年代后期至21世纪初。其中最为严重的干旱出现在2001年。根据Kendall'stau的趋势计算[15]结果分析, 整体而言, 全国干旱面积在近56年中没有显著增加或减少的趋势, 趋势值为-0.04/10a (表 2), 没有通过95%信度检验。
2.2 各区域干旱的长期变化
总体而言, 全国干旱面积在近56年中没有显著增加或减少的长期变化趋势, 但不同地区的变化特点有很大差异。下面分别讨论东北、华北、西北地区东部、西北地区西部、长江中下游、华南、西南和青藏高原近半个世纪的干旱变化特征 (图 3)。
从各区域近56年来干旱面积的变化看, 东北和华北地区有明显增加趋势 (图 3)。东北地区的干旱面积变化趋势达到2.31%/10a, 华北地区达2.50%/10a, 均通过了95%信度检验, 说明东北和华北地区近半个多世纪以来的干旱化趋势十分明显。特别是自20世纪90年代中后期至21世纪初, 东北和华北地区的干旱面积之大是近半个多世纪罕见的, 2001年两个地区的干旱面积均为历年最大, 东北和华北地区分别为58.9%和62.7%。而且在这一时期, 大范围的严重干旱持续了4~5年之久。西北地区东部的干旱面积在近半个世纪没有明显的增加或减少趋势。但在20世纪90年代中后期至21世纪初也出现了连续数年的大范围干旱。继20世纪90年代后期至21世纪初的严重干旱后, 2003年北方地区降水较常年明显偏多, 华北和西北东部地区的干旱面积出现了显著减少, 但最近几年又出现了明显增多, 加之前期连续数年的干旱, 水资源短缺的现象日益严重。西北地区西部、长江中下游、华南和西南地区的干旱面积也没有显著的变化趋势存在, 但存在着明显的年代际变化。其中西北地区西部在20世纪80年代中期以后干旱面积有较明显的减少。青藏高原地区的干旱面积有减少趋势, 趋势值为-2.14%/10a, 尤其是近20多年减少趋势明显 (图 3), 这与20世纪80年代中后期以后, 青藏高原大部分地区年降水量和相对湿度明显的增加趋势相对应[16]。
3 干旱持续时间和发生次数除了干旱程度外, 持续时间和发生频次也是干旱特征的重要组成部分。根据《气象干旱等级国家标准》[11]规定, 统计了全国606个站点1951-2006年历年发生的干旱过程, 以此来分析干旱的持续时间和发生次数的变化特征。图 4反映了1951-2006年, 全国干旱过程的最长持续时间分布。干旱的平均持续时间的空间分布 (图略) 与其最长持续时间的空间分布相似。干旱持续时间长的几个中心分别位于东北地区西部、华北地区中部和云南等地, 最长持续时间一般有4个月以上。另外, 东北地区中东部、华北大部、黄淮、华南地区中南部和西南地区东部部分地区干旱最长持续时间一般也有3个月左右, 可见干旱灾害之严重。统计表明, 对于东北、华北、西北地区东部、华南及四川盆地等地的大部分地区, 在近半个多世纪里, 最长的干旱事件大多发生在1980年以后的二十多年里 (图 5)。
分别计算1951-1978年、1979-2006年两段时间的年平均干旱日数之差 (图 6), 结果表明:东北地区西部和南部、华北、西北地区东部、华南及云南南部等地的后28年的干旱发生日数比前28年明显增多, 即干旱发生更加频繁了。在长江中下游地区和西南地区中南部干旱日数有所减少。全国其余大部分地区干旱日数没有出现明显增多或减少的变化。
4 干旱变化与区域增暖
近50年来, 干旱化趋势明显的华北、东北地区和西北地区东部部分地区增暖也十分显著。近半个世纪以来年平均气温资料分析结果反映, 我国大部分地区都出现了气温升高的现象, 其中北方大部分地区增暖趋势显著, 华北大部、东北大部及西北部分地区的增温幅度达到0.4~0.8 ℃/10a左右, 南方大部分地区没有显著的冷暖变化趋势[17]。近50年来气温上升最多的地区是东北、华北和西北的部分地区 (图 7)。从利用Kendall' stau的趋势检验计算方法[15]得到的结果可以看到 (图 7和图 8), 我国干旱化趋势最显著的地区与增暖幅度最大的地区有很大的一致性, 这说明区域增暖在干旱变化趋势中起着一定的作用。这与Dai等[2]的研究结果相一致, 即20世纪70年代后期以来, 全球变暖的气候变化进程不仅带来了地面气温的升高, 还使得一些地区干旱化趋势加剧。
另外, 干旱的发生和大范围蔓延会抑制植被的生长甚至带来植被的大量减少, 增加下垫面的沙源条件, 从而引起风沙活动加剧、沙尘天气频繁发生, 导致沙漠化进程加剧。研究表明[18], 20世纪90年代后期至21世纪初的几年北方地区连续出现的大范围严重干旱, 导致植被覆盖范围的迅速减小, 这可能是造成这几年北方地区沙尘天气频发的重要原因之一。近半个世纪以来, 我国干旱化明显的地区主要出现在北方地区, 北方地区干旱化趋势的演变会影响到地面植被的生长, 并导致生态环境的恶化[19]。虽然在2003年以后北方部分地区的降水量较之前几年有所增加, 但由于受20世纪末至21世纪初连续数年严重干旱的后续影响, 同时北方大部地区处于半干旱、半湿润气候区, 属于气候脆弱区, 干旱的威胁仍然相当严重。
5 小结综上所述, 1951-2006年, 我国的干旱总体上没有显著增加或减少的趋势, 但在不同地区存在显著的区域差异。东北和华北地区干旱有显著加重趋势, 在这些地区, 20世纪90年代后期至21世纪初的严重干旱在近半个多世纪中是十分罕见的。另外, 近50年来, 东北、华北和西北地区东部的大部分地区持续时间最长的干旱事件多发生在1980年以后的二十多年中, 同时, 干旱发生也更加频繁。另外, 近50年来, 我国北方地区区域增暖显著, 同时干6期邹旭恺等:近半个世纪我国干旱变化的初步研究685旱化趋势也最明显。
[1] | Woodhouse C A, Overpeck J T, 2000 years of drought variability in the central United States. Bull Amer Meteor Soc, 1998, 79: 2693–2714. DOI:10.1175/1520-0477(1998)079<2693:YODVIT>2.0.CO;2 |
[2] | Dai A G, Trenberth K E, Qian T T, A Global data set of Palmer Drought Severity Index for 1870--2002: Relationship with soil moisture and effects of surface warming. J Hydrometeorology, 2004, 5: 1117–1130. DOI:10.1175/JHM-386.1 |
[3] | 张强.干旱及其影响∥全国气候影响评价 (2002).北京:气象出版社, 2003:12-18. |
[4] | 邹旭恺.干旱及其影响∥全国气候影响评价 (2006).北京:气象出版社, 2007:12-18. |
[5] | 谢安, 孙永罡, 白人海. 中国东北近50年干旱发展及对全球气候变暖的响应. 地理学报, 2003, 58, (增刊): 75–82. |
[6] | 王志伟, 翟盘茂. 中国北方近50年干旱变化特征. 地理学报, 2003, 58, (增刊): 61–68. |
[7] | 马柱国, 华丽娟, 任小波. 中国近代北方极端干湿事件的演变规律. 地理学报, 2003, 58, (增刊): 69–74. |
[8] | 马柱国. 华北干旱化趋势及转折性变化与太平洋年代际振荡的关系. 科学通报, 2007, 52, (10): 1199–1206. |
[9] | 王绍武, 龚道溢, 翟盘茂.西部地区的气候变化 (第二章)//秦大河.中国西部地区的环境演变.北京:科学出版社, 2002:31-80. |
[10] | 施雅风, 沈永平, 胡汝骥. 西北气候由暖干向暖湿转型的信号、影响和前景初步探讨. 冰川冻土, 2002, 24, (3): 219–226. |
[11] | 张强, 邹旭恺, 肖风劲, 等.气象干旱等级.GB/T20481-2006,中华人民共和国国家标准.北京:中国标准出版社, 2006:1-17. |
[12] | 李庆祥, 刘小宁, 张洪政, 等. 定点观测气候序列的均一性研究. 气象科技, 2003, 31, (1): 3–10. |
[13] | 吴增祥. 气象台站历史沿革信息及其对观测资料序列均一性影响的初步分析. 应用气象学报, 2005, 16, (4): 461–467. |
[14] | Jones P D, Hulme M, Calculating regional climatic time series for temperature and precipitation: Methods and illustrations. Int J Climate, 1996, 16: 361–377. DOI:10.1002/(ISSN)1097-0088 |
[15] | Wang X L, Swail V R, Changes of extreme wave heights in Northern hemisphere oceans and related atmospheric circulation regimes. J Climate, 2001, 14: 2204–2221. DOI:10.1175/1520-0442(2001)014<2204:COEWHI>2.0.CO;2 |
[16] | 牛涛, 刘洪利, 宋燕, 等. 青藏高原气候由暖干到暖湿时期的年代际变化特征研究. 应用气象学报, 2005, 16, (6): 763–771. |
[17] | 翟盘茂, 章国材. 气候变化与气象灾害. 科技导报, 2004, 7: 11–14. |
[18] | Zou X K, Zhai P M, Relationship between vegetation coverage and spring dust storms over northern China. J Geophys Res, 2004, 109: D03104. DOI:10.1029/2003JD003913 |
[19] | 朱震达. 全球变化与荒漠化. 地学前缘, 1997, 4: 213–219. |