文献摘录_干旱_0505

1. 合肥市干旱识别及基于Copula的特征值重现期分析_周平（2020）

三阈值法干旱识别及特征值提取见图1，步骤如下。

步骤1 逐时段干旱状态初步识别

的时段为干旱，识别出a、b、c、d共4次干旱。

步骤2 干旱合并与剔除

对历时1个时段的干旱(如a、d)，若其(如a)，则最终被识别为1次干旱，否则，剔除该时段(如d)。对于间隔短于L个时段的相邻干旱，且间隔期，合并为1次干旱(如L=3，b、c可合并为一次干旱b+c)；间隔期为R>R1的时段，则为独立的干旱(如c、d)。

步骤3 干旱特征值提取

记干旱烈度为干旱起始至结束过程中R≤R0的累积水分亏缺量，峰值强度I为过程中水分亏缺最大值。、分别取《气象干旱等级》中轻旱所对应的降水距平百分率的下、上限值，即R0=-40%、R2=-60%；R1为干旱解除阈值，据降水量超过蒸发量和土壤耕作层有效含水率的时段内雨养农业不干旱，取对应R1的月降水量为月水面蒸发量与耕作层有效含水率W之和，即:

W=(θf-θw)h

式中，θf、θw分别为田间持水量、凋萎含水率(体积含水率)，%；h为耕作层厚度，mm。

3 合肥市干旱识别及特征值重现期分析

3.1 干旱识别结果

采用1953~2019年合肥站逐月降水和水面蒸发资料识别干旱。其中1953~2019年合肥站逐月降水取自中国气象数据网(http://data.cma.cn),1963~2019年月水面蒸发量取自合肥水文水资源局设在董铺的蒸发量观测资料,设站前1953~1962年的水面蒸发量以1963~2019年多年平均的月水面蒸发量代替。根据安徽省土壤含水率田间持水量数据和太湖地区典型水稻土土壤水分特征,取合肥市土壤耕作层有效含水率为22mm(太湖地区水稻土耕作层土壤有效体积含水率为15.7%,厚度为14cm[8])。合肥市1953~2019年各月平均蒸发量及R0、R2所对应的降水量见表1。

对于S<10mm干旱,在其过程中遭遇1场中雨量级的降水,所识别的干旱即不存在,故剔除S<10mm的干旱。干旱识别结果受小干旱合并的间隔参数L影响,L=1、3时识别的干旱分别为83、80场。比较识别的干旱和实际旱情,发现L=3识别的干旱更符合实情,80场干旱的S、I的分布函数的拟合见图2。

         游程理论方法通常用于识别气象水文事件，是分析气象水文灾害的时间序列分析方法，本文基于月尺度标准化降水蒸散指数(SPEI)时序性特征，并结合干旱游程理论[81]，提取锡林郭勒盟春季、夏季、秋季和植被生长季干旱事件的特征因子(如图3-14:干旱历时与干旱强度)。

        其中，干旱历时表示一次干旱事件中负游程的起止时间长度，根据锡林郭勒盟干旱指数等级判断,本文干旱阈值R、R1和R2选取为0、-0.5和-1。干旱烈度是在一次干旱过程中，指标值与干旱阈值R1之差的累积量，并无量纲值[82]。

2 锡林郭勒盟多时间尺度干旱及其对植被物候的影响_吕达（2022）

3.3.3基于游程理论方法的干旱特征分析

游程理论方法通常用于识别气象水文事件,是分析气象水文灾害的时间序列分析方法,本文基于月尺度标准化降水蒸散指数(SPEI)时序性特征,并结合干旱游程理论[81],提取锡林郭勒盟春季、夏季、秋季和植被生长季干旱事件的特征因子(如图3-14:干旱历时与干旱强度)。其中,干旱历时表示一次干旱事件中负游程的起止时间长度,根据锡林郭勒盟干旱指数等级判断,本文干旱阈值R、R1和R2选取为0、-0.5和-1。干旱烈度是在一次干旱过程中,指标值与干旱阈值R1之差的累积量,并无量纲值[82]。

        本文干旱识别流程如下:

(1)选取干旱历时与干旱烈度为干旱事件特征因子描述对象，单次干旱小于R1为一次干旱事件，初步识别此月为干旱月份。

(2)如图3-14，第二次干旱事件中,干旱指标超过R2时记为一次干旱过程。但当R2<SPEI<R1时，且干旱历时不足一个月，则将此月记为非干旱月份。

(3)合并干旱过程,如图3-14，第二次干旱事件中,当两次干旱事件间隔小于或等于一个月，且此月R1<SPEI<R0，则将两次干旱事件合并为一次干旱过程，干旱历时为干旱时间长度之和，干旱烈度为各个月份干旱烈度的累加；否则为两个不同的干旱进程。

        为研究锡林郭勒盟植被从春季生长开始期至秋季植被生长结束期之间的干旱事件过程，采用干旱游程理论对锡林郭勒盟春季、夏季、秋季及生长季干旱历时分析，如图3-15所示。

        从图3-15a可以看出，受副热带高压的影响，经过冬季的降温降雪，锡林郭勒盟春季升温并伴随着降水，锡林郭勒盟春季干旱历时在1~2.5个月之间，其干旱持续时间较短[83]。

        春季干旱历时区域平均为1.57个月，尤其以东乌珠穆沁旗及西乌珠穆沁旗以北干旱历时最长，干旱历时高值区多集中在2.33个月。

        除西珠穆沁旗中部和苏尼特右旗南部春季干旱历时较短外，其余各地干旱历时在1.4~2.33个月之间，占总区域面积的87.29%。

        由图3-15a和图3-15b可以看出，1982~2015年锡林郭勒盟春季干旱历时与夏季干旱历时空间分布差异较大，结果表明，锡林郭勒盟夏季干旱平均历时达2.05个月，历时在1.8~2.5个月之间的像元占区域面积的89.54%。

        在阿巴嘎旗、锡林浩特市及苏尼特右旗和苏尼特左旗以西等地夏季干旱历时较长，其中研究区以西荒漠草原地带夏季干旱历时最长，干旱历时超过2.5个月，占区域总面积的5.53%。

        区域东南部太卜寺旗、正蓝旗、苏尼特右旗以南等地夏季干旱历时最短，夏季干旱历时持续1.5~1.8个月的区域占总面积的8.41%。

        由图3-15c可以看出，1982~2015年锡林郭勒盟秋季干旱历时较长，平均历时达2.33个月，历时在1.53~3.63个月之间的像元占区域面积的95%。

        秋季干旱历时呈现研究区以北草甸草原区较短，南部较干旱的沙地草原及荒漠草原地区历时则较长。

        锡林郭勒盟秋季大风日数较多，降水量偏少，在苏尼特右旗、镶黄旗、正镶白旗、正蓝旗及太卜寺旗等地干旱历时较长，主要集中在2.42~3.53个月，占区域总面积的48.34%。区域东北部东乌珠穆沁旗等地秋季干旱历时持续相对较短，占区域面积的21.8%。

 

        由于锡林郭勒盟植被平均生长季长度为春季的四月至秋季的十月，为探究1982~2015年锡林郭勒盟植被生长季干旱历时，选取总计六个月的SPEI值进行提取统计分析，得到锡林郭勒盟34年植被生长季平均干旱历时(图3-15d)。由图3-15d可以看出，34年锡林郭勒盟植被生长季干旱历时在3.05~5.13个月之间，区域平均干旱历时为4.04个月，各地区间干旱历时长度有所不同。历时较高值零星分布在区域大部分地区，其中以苏尼特右旗及苏尼特左旗以西等地区生长季干旱历时最长，达到5个月以上。此区域多以荒漠草原为主，植被覆盖率较低，并伴有大面积的草地退化，干旱持续时间相对较长。区域南部正蓝旗、多伦县等草地生态较好的地区，生长季干旱历时相对较短，多集中在3~3.5个月之间，占区域面积的4%，且以多伦县的历时最短，仅为3.07个月。

        图3-16为1982~2015年锡林郭勒盟春季干旱烈度(a)、夏季干旱烈度(b)、秋季干旱烈度(c)、生长季干旱烈度(d)。

        从图3-16a可以看出，春季干旱烈度范围为7.54~13.82，烈度较大的地区多分布于东乌珠穆沁旗及西乌珠穆沁旗北部，其中较大春季干旱烈度在11~13.82之间像元空间占比为29.06%。

        春季干旱烈度较小的地区分布在研究区东南部，其烈度范围集中在7.54~9之间，占空间面积的19.11%。整体上，南北春季干旱烈度差异化明显，在春季干旱历时较大的区域，春季干旱烈度也较大。

        由图3-16b可以看出，锡林郭勒盟夏季干旱烈度由北向南呈条带状分布，夏季干旱烈度在11.77~15.44之间，烈度较大的地区多分布在区域以南荒漠草原地区，尤其在二连浩特等地夏季干旱烈度达到最高值(15.44)，同时干旱历时也最长。

        夏季干旱烈度较弱地区分布在阿巴嘎旗、东乌珠穆沁旗以北及锡林浩特市等地。

        可以发现，尽管此区域夏季干旱历时相对较长，但旱情较为平缓，干旱强度较小，烈度表现不明显。

        锡林郭勒盟秋季干旱烈度在区域西南正蓝旗、太卜寺旗、多伦县等地较大(图3-16c)，秋季干旱烈度范围在11.48~15.09之间，此区域在浑善达克沙地以南，夏季干旱历时较长，干旱烈度表现较为强烈，这一结果与图雅等[84]在锡林郭勒盟研究旱灾强度，较大旱情的区域空间结果相近。

        由图3-16d可以看出，34年锡林郭勒盟植被生长季干旱烈度在22.62~28.42之间，区域平均干旱烈度为4.04个月，空间上南北生长季干旱烈度较大，研究区中部生长季干旱烈度较小。

        从烈度高值区的分布来看，锡林郭勒盟植被生长季干旱烈度在25~28.42之间，面积占比为59.94%，烈度较高值多集中在区域以南苏尼特左旗、苏尼特右旗、镶黄旗、太卜寺旗等地，尤其在二连浩特市等地，生长季干旱烈度达到最高值。此区域生长季干旱历时均相对较长。

参考文献：

[81]肖名忠,张强,陈晓宏.基于多变量概率分析的珠江流域干旱特征研究[J].地理学报,2012,67(01):83-92.

[82]李天水,王顺,庄文化,刘铁刚.游程理论和Copula函数在二维干旱变量联合分布中的应用[J].干旱区资源与环境,2016,30(06):77-82.

[83]佟斯琴,刘桂香,武娜.1961-2010年锡林郭勒盟气温和降水时空变化特征[J].水土保持通报,2016,36(05):340-345.

3 基于多变量概率分析的珠江流域干旱特征研究_肖名忠（2012）

这篇文献里面有{干旱烈度和干旱强度计算方法}

3.2游程理论_肖名忠(2012)

        干旱一般可以用干旱历时(drought duration)、干旱严重程度(drought severity)和干旱强度(drought intensity)三个特征来表示[9],而这可以用游程理论来得到。一般而言,在一个有限取值的序列中,满足一定条件的同一符号的一个连串称之为一个“游程”,一个游程中同一符号出现的次数称之为游程的长度。有研究认为长历时的SPI小于0的事件也会导致严重的水资源供给及与干旱有关的问题,取0作为游程的截断水平[10],本文采用了与其相同的截断水平(图2)。

        当SPI≤0时定义为发生了一次干旱事件,游程长度为干旱历时D(单位为月),在游程中累积的SPI值与临界值的缺乏值为干旱严重程度(即为干旱事件中累积的干旱缺水状况)S,即,其中干旱严重程度与干旱历时的比值为干旱强度(即为干旱事件中平均的干旱缺水状况)。

        干旱历时、干旱严重程度和干旱强度,任一个可以由另外两个得到,因而可以只用其中两个特征来表示一次干旱,本文选取的是直接得到的干旱历时和干旱严重程度。

        一般来说,当将干旱历时看做离散型分布时,干旱历时的分布可通过几何分布来拟合[11],而如果将干旱历时分布看做连续型分布时,常用指数分布来拟合[10],本文用指数分布来拟合干旱历时。干旱严重程度一般用伽马分布函数来拟合[10-11],为使结果更可靠,本文还同时采用在水文领域应用比较广泛的三变量广义极值分布、三变量广义帕累托分布、皮尔逊三型分布、三变量对数正态分布、五变量Wakeby分布来拟合。各分布函数的参数估计采用线性矩估计[12],最后使用Kolmogorov-Smirnov(K-S)[13]方法进行拟合优度检验,确定最适分布函数。

参考文献：

[9]Mishra Ashok K,Singh Vijay P.A review of drought concepts.Journal of Hydrology,2010,391:202-216.

[10]Shiau J T.Fitting drought duration and severity with two-dimensional Copulas.Water Resources Management,2006,20:795-815.

[11]Mathier L,Perreault L,Bobe B.The use of geometric and gamma-related distributions for frequency analysis of waterdeficit.Stochastic Hydrology and Hydraulics,1992,6(4):239-254.

4 游程理论和Copula函数在二维干旱变量联合分布中的应用*_李天水（2016）

1.2研究方法

1.2.1干旱识别

采用游程理论从降水距平百分率指数(Pa)序列中分离出干旱事件的历时(D)和烈度(S)。Pa值计算方法和对应的干旱等级，根据《气象干旱等级》[14]确定。干旱历时为干旱事件从开始到结束持续的时间，干旱烈度为该次干旱过程降水量的总缺乏值(图1所示干旱事件的阴影面积)。根据已有的研究[6-8]，设定3个干旱指标阈值R0、R1和R2，干旱识别分为3步:

1)初步识别。当Pa值小于R1，则初判此月为干旱。

2)剔除非干旱。对于历时只有1个月的干旱，当其Pa值大于R2，则认为此月非干旱，将其剔除。

3)合并干旱。当两次相邻干旱过程仅间隔1个时段，且该时段的Pa值小于R0，则合并这两次相邻干旱过程为1次干旱，干旱历时为两次干旱历时之和加1，干旱烈度为两次干旱烈度之和;否则为两次独立干旱过程。因此，图1中共有3次干旱过程，即b、c+d、e。

由于游程理论的两类截取水平可以互相转化，故文中采用第一类截取水平识别干旱事件，并同时设定3种不同的截取水平。截取水平A为R1取-30%，截取水平B为R1取-25%，截取水平C为R1取-20%;R0皆取为0，即多年同期平均降水量的Pa值;R2皆取为-40%，即根据《气象干旱等级》发生轻旱时的Pa值。