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

—1—干旱监测指标方法及等级标准干旱监测从不同层面具有多种监测方法和指标，本附件主要介绍8种干旱监测指标及其计算方法，其中最后两种方法为气象卫星遥感干旱监测处理方法和指标。目前在卫星遥感干旱监测业务中，使用的数据源为NOAA-16卫星AVHRR数据，生成卫星遥感干旱监测产品的空间分辨率为1.1公里，监测周期为每旬一次。1、降水量（P）和降水量距平百分率(Pa)1.1原理和计算方法降水量距平百分率（Pa）是指某时段的降水量与常年同期降水量相比的百分率：%100PPPPa……………………(1)其中P为某时段降水量，P为多年平均同期降水量，本标准中取1971～2000年30年气候平均值。niiPnP11……………………(2)其中iP为时段i的降水量，n为样本数，30n。1.2等级划分由于我国各地各季节的降水量变率差异较大，故利用降水量距平百分率划分干旱等级对不同地区和不同时间尺度也有较大差别，表1为适合我国半干旱、半湿润地区的干旱等级标准。—2—表1单站降水量距平百分率划分的干旱等级等级类型降水量距平百分率(Pa)（%）（月尺度）（季尺度）（年尺度）1无旱-50＜Pa-25≤Pa-15≤Pa2轻旱-70＜Pa≤-50-50≤Pa＜-25-30≤Pa＜-153中旱-85＜Pa≤-70-70＜Pa≤-50-40＜Pa≤-304重旱-95＜Pa≤-85-80＜Pa≤-70-45＜Pa≤-405特旱Pa≤-95Pa≤-80Pa≤-452、标准化降水指数（SPI或Z）2.1原理和计算方法标准化降水指数（简称SPI）是先求出降水量分布概率，然后进行正态标准化而得，其计算步骤为：1）假设某时段降水量为随机变量x，则其分布的概率密度函数为：/1)(1)(xexxf，0x……………………(3)01)(dxexx……………………(4)其中：0，0分别为尺度和形状参数，和可用极大似然估计方法求得：AA43/411ˆ……………………(5)ˆ/ˆx……………………(6)其中niixnxA1lg1lg……………………(7)—3—式中ix为降水量资料样本，x为降水量多年平均值。确定概率密度函数中的参数后，对于某一年的降水量0x，可求出随机变量x小于0x事件的概率为：00)()(dxxfxxP……………………(8)利用数值积分可以计算用（3）式代入（8）式后的事件概率近似估计值。2）降水量为0时的事件概率由下式估计：nmxP/)0(……………………(9)式中m为降水量为0的样本数，n为总样本数。3）对分布概率进行正态标准化处理，即将（8）、（9）式求得的概率值代入标准化正态分布函数，即：02/0221)(dxexxPZ……………………(10)对（10）式进行近似求解可得：0.1))(()(123012tdtdtdctctctSZ……………………(11)其中21lnPt，P为（8）式或（9）式求得的概率，并当5.0P时，PP0.1，1S;当5.0P时，1S。515517.20c，802853.01c，010328.02c，432788.11d，189269.02d，001308.03d。由（11）式求得的Z值也就是此标准化降水指数SPI。2.2等级划分—4—由于标准化降水指标就是根据降水累积频率分布来划分干旱等级的，它反映了不同时间和地区的降水气候特点。其干旱等级划分标准具有气候意义，不同时段不同地区都适宜。表2标准化降水指数SPI的干旱等级等级类型SPI值出现频率1无旱-0.5SPI68%2轻旱-1.0SPI≤-0.515%3中旱-1.5SPI≤-1.010%4重旱-2.0SPI≤-1.55%5特旱SPI≤-2.02%3、相对湿润度指数(iM)3.1原理和计算方法相对湿润度指数是某时段降水量与同一时段长有植被地段的最大可能蒸发量相比的百分率，其计算公式：EEPMi……………………(12)式中：P为某时段的降水量，E为某时段的可能蒸散量，用FAOPenman-Monteith或Thornthwaite方法计算，计算方法见参考文献[AllenRichardG.,LuisS.Pereira,DirkRaesandMartinSmith,1998.Cropevapotranspiration-Guidelinesforcomputingcropwaterrequirements-FAOIrrigationanddrainagepaper56,FAO,Rome.]和[马柱国，符淙斌，2001.中国—5—北方地表湿润状况的年际变化趋势，气象学报，59(6)，737-746.]。3.2等级划分相对湿润度指数反映了实际降水供给的水量与最大水分需要量的平衡，故利用相对湿润度指数划分干旱等级不同地区和不同时间尺度也有较大差别，表3为适合我国半干旱、半湿润地区月尺度的干旱等级标准。表3相对湿润度指数Mi的干旱等级等级类型相对湿润度指数Mi1无旱-0.50Mi2轻旱-0.70Mi≤-0.503中旱-0.85Mi≤-0.704重旱-0.95Mi≤-0.855特旱Mi≤-0.954、综合干旱指数Ci4.１原理和计算方法气象干旱综合指数iC是以标准化降水指数、相对湿润指数和降水量为基础建立的一种综合指数：339iCZMZ……………………(13)当0iC，并10P≥0E时（干旱缓和），则iiCC5.0；当200yP毫米（常年干旱气候区，不做干旱监测），0iC。通常50EE，当5E＜5毫米时，则0E=5毫米。式中：3Z、9Z为近30和90天标准化降水指数SPI，由(11)—6—式求得；3M为近30天相对湿润度指数，由（12）式得；5E为近5天的可能蒸散量，用桑斯维特方法（ThornthwaiteMethod）计算，计算方法见参考文献[马柱国，符淙斌，2001.中国北方地表湿润状况的年际变化趋势，气象学报，59(6)，737-746.]。10P为近10天降水量，yP为常年年降水量；α、γ、β为权重系数，分别取0.4、0.8、0.4。通过（13）式，利用逐日平均气温、降水量滚动计算每天综合干旱指数iC进行逐日实时干旱监测。4.２等级划分气象干旱综合指数iC主要是用于实时干旱监测、评估，它能较好地反映短时间尺度的农业干旱情况。表4综合干旱指数iC的干旱等级等级类型iC值干旱对生态环境影响程度1无旱-0.6iC降水正常或较常年偏多，地表湿润，无旱象。2轻旱-1.2iC≤-0.6降水较常年偏少，地表空气干燥，土壤出现水分不足，对农作物有轻微影响。3中旱-1.8iC≤-1.2降水持续较常年偏少，土壤表面干燥，土壤出现水分较严重不足，地表植物叶片白天有萎蔫现象，对农作物和生态环境造成一定影响。4重旱-2.4iC≤-1.8土壤出现水分持续严重不足，土壤出现较厚的干土层，地表植物萎蔫、叶片干枯，果实脱落；对农作物和生态环境造成较严重影响，工业生产、人畜饮水产生一定影响。—7—5特旱iC≤-2.4土壤出现水分长时间持续严重不足，地表植物干枯、死亡；对农作物和生态环境造成严重影响、工业生产、人畜饮水产生较大影响。5、土壤墒情干旱指数５.1原理和计算方法（1）土壤重量含水率计算公式为：%100ddwmmmW……………………(14)式中W为土壤重量含水量，mw为湿土重量，md为干土重量。（2）土壤田间持水量测定和计算方法多采用田间小区灌水法：选择4平方米的小区（2m×2m），除草平整后，做土埂围好；对小区进行灌水，灌水量的计算公式如下：100)(2hswaQ……………………(15)式中Q为灌水量（m3）；α为假设所测土层中的平均田间持水量（%），一般沙土取20%，壤土25%，粘土取27%；ω为灌水前的土壤湿度（%）；ρ为所测深度的土壤容重（m3/m3），一般取1.5；s为小区面积（m2）；h为测定的深度（m）；2为小区需水量的保证系数。在土壤排除重力水后，测定土壤湿度，即田间持水量。土壤排除重力水的时间因土质而异，一般沙性土需1～2天，壤性土需2～3天，粘性土需3～4天。在测定土壤湿度时，每天取样一次，每次取4个重复的平均值，当同一层次前后两次测定的土壤湿度差值＜2.0%时，则第2次的测定值即为该层的田间持水量。（3）土壤相对湿度计算方法：—8—%100cfwR……………………(16)w为土壤重量含水率（%），cf田间持水量(用重量含水率表示)。5.2等级划分表5土壤相对湿度的干旱等级等级类型20厘米深度土壤相对湿度对农作物影响程度1无旱正常R＞60%地表湿润，无旱象2轻旱60%≥R＞50%地表蒸发量较小，近地表空气干燥3中旱50%≥R＞40%土壤表面干燥，地表植物叶片白天有萎蔫现象4重旱40%≥R＞30%土壤出现较厚的干土层，地表植物萎蔫、叶片干枯，果实脱落5特旱R≤30%基本无土壤蒸发，地表植物干枯、死亡6．Palmer干旱指数6.1Palmer干旱指数计算方法Palmer干旱指数可以简称为PDSI(PalmerDroughtSeverityIndex)。它的计算可分为以下七个步骤：一是统计水文帐，由长期气象资料序列计算出月水分平衡各分量的实际值、可能值及平均值，包括蒸散量、潜在蒸散量、径—9—流量、潜在径流量、补水量、潜在补水量、失水量和潜在失水量；二是计算各气候常数和系数，包括蒸散系数、补水系数、径流系数、失水系数和气候特征值；三是计算出水分平衡各分量的气候适宜值，包括气候适宜蒸散量、气候适宜补水量、气候适宜径流量、气候适宜失水量和气候适宜降水量；四是计算水分盈亏值d和水分异常指数Z；五是计算Palmer干旱指数；六是干期（或湿期）结束的度的计算，确立当前期的Palmer指数值。6.2Palmer干旱指数的计算公式PDSI的原理是水分平衡方程，即在“当前情况下达到气候上适宜”（ClimaticallyAppropriateForExistingConditions）的情况下，降水量等于蒸散量与径流量之和再加上（或减去）土壤水分的交换量，表示如下：LRROETP……………………(17)ET为蒸散量，PET为可能蒸散量，R为土壤水分供给量，PR为土壤可能水分供给量，RO为径流，PRO为可能径流，L为土壤水分损失量，PL为土壤可能水分损失量。其中PETET……………………(18)PRR……………………(19)PRORO……………………(20)—10—PLL……………………(21)、、、分别为蒸散系数、土壤水供给系数、径流系数和土壤水损失系数，每地每月分别有四个相应的常系数值，计算如下：)()(PETET……………………(22))()(PRR……………………(23))()(PRORO……………………(24))()(PLL……………………(25)各量上面的横线代表其多年平均值。Palmer指数假定土壤为上下两层模式，除非上层土壤中的水分全部丧失，下层土壤才开始失去水分，且下层土壤的水分不可能全部失去。可能蒸散量PET用Thornthwaite方法计算（计算PET的方法在PDSI程序中是一个独立的部分，可以用其他方法代替）。在计算蒸散量、径流量、土壤水分交换量的可能值与实际值时，需要遵循一系列的规则和假定。另外，土壤有效持水量AWC（AvailableWaterHoldingCapacity）也作为初始输入量。在计算PDSI过程中，实际值与正常值相比的水分距平d表示为实际降水量与气候适宜下降水量的差：PPd……………………(26)水分距平d求出后，又将其与指定地点给定月份的气候权重—11—系数K相乘，得出Palmer-Z指数：dKZ……………………(27)Z被定义为水分异常指数，表示给定地点给定月份，实际气候干湿状况与其多年平均水分状态的偏离程度。其中气候权重系数K的计算公式如下：'12