华南地处低纬，濒临南海，受低纬度热带天气系统和中高纬度天气系统的交替影响，天气气候复杂多变。华南虽然年降水比较充沛，属于湿润气候区^[1]，但时空分布不均，兼之太阳辐射强、气温高、土壤蒸发和作物蒸腾强烈，区域性、季节性干旱十分突出，华南沿海地区是我国五大气象干旱中心之一^{[2, 3, 4]}，也是气候变化敏感区和脆弱区之一。在全球变化背景下，极端干旱气候事件呈频发趋势，全球水资源和水循环研究已经成为全球变化科学的新的战略重点^[5]。

文献中对华南干旱的研究可分为阶段降水评估和逐日干旱监测评估两类。降水评估类研究包括对华南区域代表站点年度^[6]、季度^{[7, 8]}、汛期(前汛期)^{[9, 10]}等不同时间尺度降水的分析评定，研究方法包括极差法^[6]、距平法^{[8, 9]}、标准化降水指数^{[7, 10]}等。逐日干旱监测评估研究近年日益得到重视，广东干旱逐日监测评估中发展了逐日土壤水分模拟方法^{[11, 12, 13, 14]}，近年随着综合气象干旱指标CI^[15]及其改进指标^[16]在全国气象干旱监测评估中广泛应用^{[3, 4, 17, 18]}，王春林等^[19]从气象干旱定义“降水持续偏少导致的水分亏缺现象”出发，考虑干旱累积效应^[2]，提出基于标准化前期降水指数(SAPI)的逐日气象干旱指标，克服了基于“等权累加”建立的标准化降水指数(SPI)和综合干旱指数(CI)等指标由于前期降水移出计算窗口而导致的“不合理旱情加剧”问题^[19]。鉴于SAPI本质上是标准化变量，不能充分刻画干旱频率季节性、区域性差异，王春林等^[20]基于SAPI和常年平均相对湿润度指数(M)构建了逐日干旱指数(DI)，并基于DI编制出版了《广东省气象干旱图集》^[20]。DI干旱指标通过广东省气象局业务准入评审，成为广东省气象部门现行气象干旱指标标准。本文采用DI分析华南区域近50年干旱时空特征及其气候变化趋势，期望对进一步开展气候变化影响评估、水资源利用及应对气候变化研究具有参考意义。

本文中的华南包括广东和广西两省。气象资料包括广东86个、广西88个地面气象站从建站至2010年的逐日降水和可能蒸散，其中可能蒸散通过逐日气温、日照时数、相对湿度、风速等资料采用FAO Penman-Monteith公式^[15] 计算。资料来自广东省气象信息中心。

单站逐日气象干旱指数DI表征逐日气象干旱强度，可以刻画逐日气象干旱发生、发展和结束过程。DI定义为：

点击浏览公式

式(1)中SAPI_i是第i日前期降水指数API的标准化变量SAPI^[19]，SAPI计算方法参见《气象干旱等级》国家标准^[15]附录C标准化降水指数SPI，历史样本资料为近30年(1981—2010年)逐日API。API计算公式为：

点击浏览公式

式中,API_i为第i日API，P_i为当日降水量(mm)，API_i-1为前一日的API，k为衰减系数，取经验值0.955。每个站从建站开始逐日滚动计算API，初始API设为0。建站开始后的头4个月API受边界效应影响舍弃不用。

式(1)中 M_i为第i日常年平均相对湿润度指数，M_i根据30a(1981—2010年)逐日降水和可能蒸散量计算^[15]：

点击浏览公式

式中,P_i为第i日30a(1981—2010年)平均降水量(mm)，PE _i为第i日30a(1981—2010年)平均可能蒸散量(mm)，采用FAO Penman-Monteith方法计算^[15]。为了增强M_i年变化曲线平滑性，对P_i、PE_i做2次30d滑动平均处理。M_i理论范围为-1至∞，为了避免因为PE_i接近于0时候导致M_i趋于∞,当M_i>0时采用双曲正切函数式(4),约束M_i变化范围为-1— 1:

点击浏览公式

M_i表征某地常年平均干湿程度及其年变化动态。华南代表性站点30a(1981—2010年)平均逐日相对湿润度指数年变化曲线表明(图 1)，华南各地干湿季节性变化显著，2—3月自北向南渐次进入湿润季节(M_i>0)，10月以后自北向南渐次进入相对干季(M_i ＜0)。区域平均湿润天数(M_i>0)达64%。

图 1 华南代表站点30a(1981—2010年)平均逐日相对湿润度指数M Fig. 1 Daily relative moisture index averaged over 30 years from 1981 to 2010 in representative stations of South China

图选项

月干旱指数MI和年干旱指数YI表征月、年尺度气象干旱强度，其定义为某站月/年内小于0的DI之和除以月/年总天数，即：

点击浏览公式

式中,n为月内总天数，DI_i为逐日干旱指数。YI计算公式同式(5)，n为年内总天数，YI等同于年内各月MI的算术平均。MI、YI用于监测评估月、年尺度气象干旱程度。月气象干旱等级划分标准见表 1。统计表明各等级旱月频率与旱日频率比较接近，因此本文仅讨论各等级旱日频率，旱月频率不重复介绍。

华南平均轻旱、中旱、重旱和特旱日频率分别为12.3%、8.1%、4.2%和1.4%，总旱日频率为26.0%(表 2)，低于基于SPI指标的理论频率(轻旱、中旱、重旱和特旱日频率分别为15%、9.2%、4.4%和2.3%)^[15],这与华南地区总体上比较湿润，年内湿润时段居多有关。

各等级旱日频率具有和降水相反的年变化特征(图 2)，汛期(4—9月)旱日频率低于非汛期(10月至次年3月)。年内最旱的3个月为11、12月和次年1月，这3个月平均降水都少于50 mm，平均月干旱指数MI都小于-0.5，旱日频率都大于40%。值得注意的是，1月份旱日频率最高(49%)、MI最小(-0.63)，比降水最少的12月(30.8 mm)滞后1个月，表明DI指标能够表征气象干旱相对枯水时段的滞后效应。

图 2 近30年(1981—2010年)华南逐月平均旱日频率和降水量 Fig. 2 Monthly drought days and rain of South China(1981—2010 mean)

图选项

华南各等级旱日频率分布总体比较一致，均呈西高东低分布(图 3)。各等级旱日频率高于年平均的区域,即轻旱日频率>12%、中旱日频率>8%、重旱日频率>4%、特旱日频率>1.5的区域,主要分布在广西除桂林以外的大部分地区、广东雷州半岛及沿海地区。年干旱指数YI分布图反映了华南气象干旱综合情况，YI＜-0.3区域和上述各等级旱日频率较高区域相对应。YI＜-0.4区域包括广西西南部的百色、崇左、南宁地区和广东的雷州半岛地区，为华南最旱区域。

比较年干旱指数YI与年降水量分布细节发现，YI2.1 旱日频率年变化特征-0.3的偏旱区域和年降水量少于1600 mm的少雨区域并不完全一致。华南有5个年降水量大于2000mm的降水集中区，即广西的防城-钦州、桂林和广东的阳江、清远、汕尾，值得注意的是这5个降水集中区中，广东汕尾、广西防城-钦州2个区域YI较低，汕尾附近重旱频率较高(图 3)，表明这2个区域尽管年降水较多，但季节分配不均，具有旱涝并重特点。

图 3 近30年(1981—2010年)平均华南轻旱，中旱,重旱,极旱日频率年干旱指数(YI)及降水分布 Fig. 3 Spatial distributions of drought days(%) at light，middle,severe and extreme level,Year drought index and rain in South China averaged over 1981—2010

图选项

近50年(1961—2010年)华南最旱的5a(YI≤-0.5)为：1963、1991、2004、2009、1977，次旱的8年依次为(YI≤-0.4)：1996、2005、1966、1989、2007、1971、1992、1980(图 4)。

图 4 华南1961—2010年年干旱指数(YI) Fig. 4 Yearly drought index(YI)of South China(1961—2010)

图选项

趋势分析表明，近50年华南年平均降水量呈增加趋势，增加强度为14.1mm/10a,轻旱、特旱日呈减少趋势，中旱、重旱呈增加趋势。尽管上述变化趋势都没有达到显著性水平(表 3)，但具有季节性、区域性差异。

从季节性差异看，1—9月降水和MI以增加趋势为主，各等级旱日以减少趋势为主，其中7月份MI增加趋势、各等级旱日减少趋势均达到0.1以上显著水平(表 3，图 5)；而10—12月降水和MI以减小趋势为主，各等级旱日以增加趋势为主，其中10月份降水减少趋势达到0.05显著水平，11月份MI减小趋势及中旱、重旱、总旱日增加趋势达到0.01显著水平(表 3，图 5)。

图 5 近50年(1961—2010年)华南7月、11月月干旱指数(MI)变化趋势 Fig. 5 Linear trends of monthly drought index(MI)of July and Nevomber(1961—2010)of South China (1961—2010)

图选项

从区域分布看，年总旱日趋于增加、减少的站点数各占60%、40%，有11%的站点达到0.1以上显著性检验，不同等级旱日数变化趋势空间分布不同(图 6)。

图 6 华南1961—2010年年干旱指数(YI) Fig. 6 Yearly drought index(YI)of South China(1961—2010)

图选项

轻旱日趋于增加、减少的站点数各占49%、51%，其中增加趋势显著的有5个站(均在广西)，减少趋势显著的有5个站(大部分在广东)(图 6)。

中旱日趋于增加、减少的站点数各占68%、32%，其中增加趋势显著的有17个站(广西中东部)，减少趋势显著的有1个站(茂名)(图 6)。对各站11月份中旱日趋势分析表明，有92%站点呈增加趋势，其中58个站通过0.1信度检验(图略)。

重旱日趋于增加、减少的站点数各占86%、14%，其中增加趋势显著的有30个站(分布在广西中东部和广东的北部)，减少趋势显著的没有(图 6)。对各站11月份重旱日趋势分析表明，有91%站点呈增加趋势，其中59个站通过0.1信度检验(图略)。

特旱日趋于增加、减少的站点数各占36%、64%，其中增加趋势显著的有7个站(集中在广西的崇左、南宁、桂林及广东的东部地区)，减少趋势显著的有10个站(分布在广东东部和北部地区)(图 6)。

综上，各等级旱日显著增加的站点大多集中在广西，旱日显著减少的站点主要集中在广东，表明广西旱于广东的格局可能进一步加剧。

(1)华南近30年(1981—2010年)基于DI指标的总旱日频率平均为26.0%，其中轻旱、中旱、重旱和特旱日分别为12.3%、8.1%、4.2%和1.4%。各等级旱日频率低于基于SPI指标的理论频率(总旱日频率30.8%)^{[14, 15]},表明华南气候总体上比较湿润。

(2)华南各等级旱日频率，具有非汛期(10月至次年3月)高于汛期(4月至9月)、广西高于广东的特征，表明DI指标由于引进平均相对湿润度指数，克服了SAPI、SPI等标准化降水指数局限，能够表征气象干旱季节性、区域性特点。DI指标和SAPI^[19]、改进的综合干旱指数CInew^[16]、土壤水分模拟的指标^{[11, 13, 14]}一样，能够表征气象干旱相对枯水时段的滞后效应。

(3)近50年(1961—2010年) 华南最旱的5a(YI≤-0.5)依次为：1963、1991、2004、2009、1977，与基于年降水量极差法^[6]、季度和汛期降水距平^{[7, 8, 9]}分析的干旱年基本一致。

华南平均年干旱指数YI没有显著的年际变化趋势，降水量有增加趋势，轻旱、特旱日呈减少趋势，中旱、重旱呈增加趋势，但未达到0.1信度水平。

(4)气候变化趋势分析表明，从季节差异看，1至9月降水和MI以增加趋势为主，各等级旱日数减少趋势为主，其中7月份MI增加趋势及各等级旱日减少趋势均达到0.1显著水平；而10至12月降水和MI以减小趋势为主，各等级旱日数以增加趋势为主，其中11月份MI减小趋势及中旱、重旱、总旱日增加趋势均达到0.05显著水平。总体趋势和黄晚华等^[7]基于SPI对华南的研究结果一致。

(5)从区域分布看，年总旱日趋于增加、减少的站点数各占60%、40%，有11%的站点达到0.1以上显著性检验。各等级旱日显著增加的站点大多集中在广西，而旱日显著减少的站点主要集中在广东，表明广西干旱总体上重于广东的格局可能进一步加剧。

气象干旱相对农业干旱、水文干旱和社会经济干旱具有前导性、基础性特点，气象干旱时空特征与所采用的指标体系密切相关，如何客观评价气象干旱指标适用性，学术界尚没有统一观点。进一步研究将结合土壤水分观测资料，分析DI指标敏感性特征和适用性特征，并结合农作物生育期、干旱灾情资料，开展旱涝灾害影响评估及风险区划研究。