干旱监测指标方法及等级标准

1、干旱监测指标方法及等级标准 干旱监测从不同层面具有多种监测方法和指标，本附件主要介绍8种干旱监测指标及其计算方法，其中最后两种方法为气象卫星遥感干旱监测处理方法和指标。目前在卫星遥感干旱监测业务中，使用的数据源为NOAA-16卫星AVHRR数据，生成卫星遥感干旱监测产品的空间分辨率为1.1公里，监测周期为每旬一次。1、 降水量（）和降水量距平百分率()1.1 原理和计算方法 降水量距平百分率（）是指某时段的降水量与常年同期降水量相比的百分率： (1)其中为某时段降水量，为多年平均同期降水量，本标准中取19712000年30年气候平均值。 (2)其中为时段的降水量，为样本数，。1.2 等级划分由

2、于我国各地各季节的降水量变率差异较大，故利用降水量距平百分率划分干旱等级对不同地区和不同时间尺度也有较大差别，表1为适合我国半干旱、半湿润地区的干旱等级标准。表1 单站降水量距平百分率划分的干旱等级等级类型降水量距平百分率(Pa)（%）（月尺度）（季尺度）（年尺度）1无旱-50Pa-25Pa-15Pa2轻旱-70Pa-50-50Pa-25-30Pa-153中旱-85Pa-70-70Pa-50-40Pa-304重旱-95Pa-85-80Pa-70-45Pa-405特旱Pa-95Pa-80Pa-452、标准化降水指数（SPI或Z）2.1 原理和计算方法标准化降水指数（简称SPI）是先求出降水量分布

3、概率，然后进行正态标准化而得，其计算步骤为： 1）假设某时段降水量为随机变量x，则其分布的概率密度函数为：， (3) (4)其中：，分别为尺度和形状参数，和可用极大似然估计方法求得： (5) (6)其中 (7)式中为降水量资料样本，为降水量多年平均值。确定概率密度函数中的参数后，对于某一年的降水量，可求出随机变量小于事件的概率为： (8)利用数值积分可以计算用（3）式代入（8）式后的事件概率近似估计值。 2）降水量为0时的事件概率由下式估计： (9)式中m为降水量为0的样本数，n为总样本数。 3）对分布概率进行正态标准化处理，即将（8）、（9）式求得的概率值代入标准化正态分布函数，即： (10

4、)对（10）式进行近似求解可得： (11)其中 ，为（8）式或（9）式求得的概率，并当时，;当时，。，。由（11）式求得的值也就是此标准化降水指数SPI。2.2 等级划分由于标准化降水指标就是根据降水累积频率分布来划分干旱等级的，它反映了不同时间和地区的降水气候特点。其干旱等级划分标准具有气候意义，不同时段不同地区都适宜。表2 标准化降水指数SPI的干旱等级等级类型SPI值出现频率1无旱-0.5<SPI 68%2轻旱-1.0<SPI-0.515%3中旱-1.5<SPI-1.010%4重旱-2.0<SPI-1.55%5特旱 SPI-2.02%3、 相对湿润度指数()3.1

5、 原理和计算方法相对湿润度指数是某时段降水量与同一时段长有植被地段的最大可能蒸发量相比的百分率，其计算公式： (12)式中： 为某时段的降水量，为某时段的可能蒸散量，用FAO Penman-Monteith或Thornthwaite 方法计算，计算方法见参考文献Allen Richard G., Luis S. Pereira, Dirk Raes and Martin Smith,1998. Crop evapotranspiration - Guidelines for computing crop water requirements - FAO Irrigation and drain

6、age paper 56, FAO ,Rome.和马柱国，符淙斌，2001. 中国北方地表湿润状况的年际变化趋势，气象学报，59(6)，737-746. 。3.2 等级划分相对湿润度指数反映了实际降水供给的水量与最大水分需要量的平衡，故利用相对湿润度指数划分干旱等级不同地区和不同时间尺度也有较大差别，表3为适合我国半干旱、半湿润地区月尺度的干旱等级标准。表3 相对湿润度指数Mi的干旱等级等级类型相对湿润度指数Mi1无旱-0.50< Mi2轻旱-0.70< Mi -0.503中旱-0.85< Mi -0.704重旱-0.95< Mi -0.855特旱 Mi -0.954、

7、综合干旱指数Ci 4.原理和计算方法气象干旱综合指数是以标准化降水指数、相对湿润指数和降水量为基础建立的一种综合指数： (13)当，并时（干旱缓和），则；当毫米（常年干旱气候区，不做干旱监测），。通常，当5毫米时，则=5毫米。式中： 、为近30和90天标准化降水指数SPI，由(11)式求得；为近30天相对湿润度指数，由（12）式得；为近5天的可能蒸散量，用桑斯维特方法（Thornthwaite Method）计算，计算方法见参考文献马柱国，符淙斌，2001. 中国北方地表湿润状况的年际变化趋势，气象学报，59(6)，737-746.。为近10天降水量，为常年年降水量；、为权重系数，分别取0.4

8、、0.8、0.4。通过（13）式，利用逐日平均气温、降水量滚动计算每天综合干旱指数进行逐日实时干旱监测。4. 等级划分气象干旱综合指数主要是用于实时干旱监测、评估，它能较好地反映短时间尺度的农业干旱情况。表4 综合干旱指数的干旱等级等级类型值干旱对生态环境影响程度1无旱-0.6<降水正常或较常年偏多，地表湿润，无旱象。2轻旱-1.2< -0.6降水较常年偏少，地表空气干燥，土壤出现水分不足，对农作物有轻微影响。3中旱-1.8< -1.2降水持续较常年偏少，土壤表面干燥，土壤出现水分较严重不足，地表植物叶片白天有萎蔫现象，对农作物和生态环境造成一定影响。4重旱-2.4<

9、-1.8土壤出现水分持续严重不足，土壤出现较厚的干土层，地表植物萎蔫、叶片干枯，果实脱落；对农作物和生态环境造成较严重影响，工业生产、人畜饮水产生一定影响。5特旱 -2.4土壤出现水分长时间持续严重不足，地表植物干枯、死亡；对农作物和生态环境造成严重影响、工业生产、人畜饮水产生较大影响。5、 土壤墒情干旱指数.1 原理和计算方法（1）土壤重量含水率计算公式为： (14)式中W为土壤重量含水量，mw为湿土重量，md为干土重量。（2）土壤田间持水量测定和计算方法多采用田间小区灌水法：选择4平方米的小区（2m×2m），除草平整后，做土埂围好；对小区进行灌水，灌水量的计算公式如下： (15)

10、式中Q为灌水量（m3）；为假设所测土层中的平均田间持水量（%），一般沙土取20%，壤土25%，粘土取27%；为灌水前的土壤湿度（%）；为所测深度的土壤容重（m3/m3），一般取1.5；s为小区面积（m2）；h为测定的深度（m）；2为小区需水量的保证系数。在土壤排除重力水后，测定土壤湿度，即田间持水量。土壤排除重力水的时间因土质而异，一般沙性土需12天，壤性土需23天，粘性土需34天。在测定土壤湿度时，每天取样一次，每次取4个重复的平均值，当同一层次前后两次测定的土壤湿度差值2.0%时，则第2次的测定值即为该层的田间持水量。（3）土壤相对湿度计算方法： (16)为土壤重量含水率（%），田间持水量

11、(用重量含水率表示)。5.2 等级划分表5 土壤相对湿度的干旱等级等级类型20厘米深度土壤相对湿度对农作物影响程度1无旱正常R60%地表湿润，无旱象2轻旱60%R50%地表蒸发量较小，近地表空气干燥3中旱50%R40%土壤表面干燥，地表植物叶片白天有萎蔫现象4重旱40%R30%土壤出现较厚的干土层，地表植物萎蔫、叶片干枯，果实脱落5特旱R30%基本无土壤蒸发，地表植物干枯、死亡6Palmer干旱指数6.1 Palmer干旱指数计算方法Palmer干旱指数可以简称为PDSI(Palmer Drought Severity Index)。它的计算可分为以下七个步骤：一是统计水文帐，由长期气象资料序

12、列计算出月水分平衡各分量的实际值、可能值及平均值，包括蒸散量、潜在蒸散量、径流量、潜在径流量、补水量、潜在补水量、失水量和潜在失水量； 二是计算各气候常数和系数，包括蒸散系数、补水系数、径流系数、失水系数和气候特征值；三是计算出水分平衡各分量的气候适宜值，包括气候适宜蒸散量、气候适宜补水量、气候适宜径流量、气候适宜失水量和气候适宜降水量；四是计算水分盈亏值d和水分异常指数Z ；五是计算Palmer干旱指数；六是干期（或湿期）结束的度的计算，确立当前期的Palmer指数值。6.2 Palmer干旱指数的计算公式PDSI的原理是水分平衡方程，即在“当前情况下达到气候上适宜”（Climaticall

13、y Appropriate For Existing Conditions）的情况下，降水量等于蒸散量与径流量之和再加上（或减去）土壤水分的交换量，表示如下： (17)为蒸散量，为可能蒸散量，为土壤水分供给量，为土壤可能水分供给量，为径流，为可能径流，为土壤水分损失量，为土壤可能水分损失量。其中 (18) (19) (20) (21)、分别为蒸散系数、土壤水供给系数、径流系数和土壤水损失系数，每地每月分别有四个相应的常系数值，计算如下： (22) (23) (24) (25)各量上面的横线代表其多年平均值。Palmer指数假定土壤为上下两层模式，除非上层土壤中的水分全部丧失，下层土壤才开始失去

14、水分，且下层土壤的水分不可能全部失去。可能蒸散量PET用Thornthwaite方法计算（计算PET的方法在PDSI程序中是一个独立的部分，可以用其他方法代替）。在计算蒸散量、径流量、土壤水分交换量的可能值与实际值时，需要遵循一系列的规则和假定。另外，土壤有效持水量AWC（Available Water Holding Capacity）也作为初始输入量。在计算PDSI过程中，实际值与正常值相比的水分距平d表示为实际降水量与气候适宜下降水量的差： (26)水分距平d求出后，又将其与指定地点给定月份的气候权重系数K相乘，得出Palmer-Z指数： (27)Z被定义为水分异常指数，表示给定地点给定

15、月份，实际气候干湿状况与其多年平均水分状态的偏离程度。其中气候权重系数K的计算公式如下： (28) (29) (30)其中i为第i个月，K随着地点和月份而变化。K值确定后，根据公式（31）来决定月PDSI值，计算公式如下： (31)公式（31）中的0.897和为持续因子，它们实际上影响到了PDSI对水分盈缺的灵敏度，起始月份的PDSI计算公式为 (32)在得到上述的统计经验计算公式之后，利用当前干（湿）期结束概率（）来决定实际的PDSI值。它的意义是计算需要得到（失去）多少水分才能使当前干（湿）期回到正常状态（即PDSI值在-0.5和0.5之间），设需要得到（失去）的水分为ZE，则有 (33)

16、 (34)由此可以得到ZE (35)这里认为大于-0.15的Z值对于结束干旱有效，定义有效增湿量， (36)同理，定义有效增干量， (37)使当前干（湿）期结束的概率可以看作是得到（失去）的水分与确切地结束当前干（湿）期所需要得到（失去）的水分的百分比率。然而，的计算不能简单地用或与ZE相比，因为在一段比常年湿的趋势中可能会出现一个干月，由此带来一个负的，所以这里又设计了一个有效增湿（干）累积量V，计算如下：若使当前干期结束 (38)若使当前湿期结束 (39)则当前干（湿）期结束的概率由下式得到 (40)在一个干（湿）期里，当达到100%时，则开始一个回算（Backtracking）计算过程，

17、以此来决定PDSI值。这里要注意到，在前面计算PDSI值时，实际上设立了三个指数X1、X2和X3来分别统计，三个指数的值都是按照式（31）来计算，X1为初始湿期的PDSI值，X2为初始干期的PDSI值，X3为当前确立干（湿）期的PDSI值，实际的PDSI值是根据一系列的规则从X1、X2和X3之中挑选出来的，而这些规则就是通过来决定的Palmer, W.C., 1965. Meteorological Drought. Research Paper No.45, US Weather Bureau （NOAA Library and Information Services Division,

18、Washington, D.C. 20852）.、Karl, T.R., 1983. Some spatial characteristics of drought duration in the United States. J. Climate Appl. Meteor., 22, 1356-1366.和Karl, T. R., 1986. The Sensitivity of the Palmer Drought Severity Index and Palmers Z-Index to their Calibration Coefficients Including Potential

19、 Evapotranspiration. J. Climate Appl. Meteor., 25, 77-86.。6.2 Palmer指数干旱等级及划分命名表6 Palmer指数干旱等级划分命名级别干旱等级指数范围名称危害程度12345无旱轻旱中旱重旱特旱无危害轻微危害中等危害严重危害特重危害X-0.99-1.00X-1.99-2.00X-2.99-3.00X-3.99X-4.007土壤湿度遥感模型（热惯量法）7.1 原理和公式该模型基于土壤的热惯量，它是度量土壤热惰性的物理量，反映土壤阻止热变化的能力。土壤热惯量随土壤密度、热传导率、热容量的增加而增加，而土壤密度、热传导率、热容量等特性的

20、变化在一定条件下取决于土壤含水量的变化，所以土壤热惯量与土壤含水量之间呈正相关关系。此外，土壤湿度控制着土壤表层温度日较差，土壤日较差与土壤含水量之间呈负相关关系，土壤日较差可以利用卫星遥感数据获得。因此，对于裸土和低植被覆盖区域，可利用气象卫星数据和实测土壤墒情资料，运用热惯量模型反演土壤表层湿度。用统计方法建立土壤湿度遥感模型，主要有线形模型和幂函数模型，在业务应用中为了简化计算直接使用日校差，拟合公式为： (41) (42)其中为土壤湿度，为拟合系数，为NOAA16卫星四通道（10.511.5µm）白天和夜间亮温差。 对NOAA-16卫星资料进行定标、定位、几何纠正等处理后，计算每日