基于游程理论的农业干旱指数在辽宁省春玉米旱灾损失评估中的应用*

谭方颖, 何 亮, 吕厚荃, 宋迎波, 程 路, 侯英雨

基于游程理论的农业干旱指数在辽宁省春玉米旱灾损失评估中的应用*

谭方颖, 何 亮, 吕厚荃, 宋迎波, 程 路, 侯英雨**

(国家气象中心 北京 100081)

为准确识别农业干旱事件, 评估干旱对农作物产量的影响, 本文以辽宁省为例, 基于52个气象站1961—2015年逐日气象数据, 计算了春玉米全生育期水分亏缺距平指数, 利用游程理论构建了一种新的春玉米全生育期干旱指数, 再结合历年产量损失率构建了区、县级尺度的干旱产量损失评估模型, 并对不同干旱风险下春玉米产量损失进行了评估, 以明确辽宁省春玉米干旱重点防范区域。研究结果表明, 辽宁省春玉米干旱指数呈由西向东递减的经向分布特征, 辽西地区更易发生严重的干旱事件, 且春玉米产量稳定性也最差。春玉米主产区内, 80%以上的区、县春玉米产量损失率与干旱指数呈显著线性正相关, 辽中部分区县和辽东大部两者相关性不显著。相同干旱指数下, 辽西产量损失率最大, 并以朝阳地区为最; 辽南和辽北次之; 辽东湾和中部部分平原地区总体较小; 非水分限制区辽东的凤城市和东港市, 降水偏少反而更有利于春玉米产量提高。辽宁省春玉米干旱重点防范区域主要分布在辽西的朝阳地区, 以及千山山脉以西和以北的辽北、辽南地区, 提高和稳定春玉米产量需增加耐旱品种种植、发展节水灌溉、提高水分利用效率; 降水充沛或灌溉条件较好、产量稳定性较高的辽东大部和辽中部分区县, 可通过提高种植和管理水平、加强其他气象灾害防御等增加春玉米产量。

农业干旱指数; 产量损失评估; 游程理论; 春玉米; 辽宁省

旱灾是世界上也是我国最严重的自然灾害之一。据中国国家统计局数据统计, 1978—2015年我国年平均干旱受灾面积是水灾、风雹、冷冻灾害受灾面积的2~7倍, 对农业的影响远超其他灾害。多年来众多学者做了大量关于干旱的研究[1-6], 包括构建干旱指数、分析干旱特征和风险、评估灾害影响等。与洪涝、冰冻等灾害不同, 干旱发生过程缓慢、反复, 因此从发生到结束往往难以准确识别, 而准确地识别干旱是构建干旱指数、对干旱进行定量分析的前提。农业干旱受到多种因素的影响, 发生机制复杂, 评价方法众多, 目前比较常用的干旱评价指数有Palmer干旱指数、Z指数、作物湿度指数、综合气象干旱指数、遥感指数等[7-13]。以上方法各有优缺点, 其中气象干旱指数无法表达农作物在不同生长发育阶段需水量与抵御干旱能力之间的差异; 遥感指数通常反映的是表层土壤的水分状况, 其监测结果也易受天气因素影响。此外, 多数作物干旱指数在空间上的可比性较差, 科学地构建干旱指数是研究干旱的关键。干旱的影响是非结构性的[14], 具有累积性、多样性和区域性, 而农业干旱的影响研究最终要归结于干旱所造成作物产量下降的程度和范围, 其评估结果除与温度、降水、地形等自然条件有关外, 同时还受到经济发展、管理水平等因素影响, 因此精细化的灾损评估是农业干旱影响研究的方向和落脚点。

游程理论对干旱持续时间以及干旱强度进行了详细的界定, 被广泛应用于水文和气象干旱事件的识别以及干旱监测和评估[15-16], 但其在农业干旱及其损失评估方面的应用研究尚少有发现。作物水分亏缺距平指数[17]既综合考虑了气象、作物、土壤等因素的影响, 又能反映出降水量和作物需水量的匹配关系, 且在时间和空间上也具有较好的可比性。辽宁省地形复杂, 水资源分布不均, 干旱是影响农业生产的首要灾害, 且有研究表明, 降水的空间分布可直接导致灾损程度区域间的差异[8]。综上, 本文以作物水分亏缺距平指数为基础, 利用游程理论从干旱指数序列中提取农业干旱事件, 基于两者构建新的干旱指数, 结合产量损失率构建区、县尺度的春玉米干旱产量损失评估模型, 评估不同风险等级干旱产量损失, 以期探索一种农业干旱评价指数和精细化旱灾损失评估方法, 为辽宁省春玉米干旱风险管理、旱灾保险以及区域防灾减灾工作等实施提供科学依据。

1 研究区域概况与研究方法

1.1 研究区域

考虑到不同地区地形、土壤、降水量等存在差异, 根据产量水平相近、降水特征相似、比邻成片的原则, 参考《中国气象地理区划手册》[18], 本文将辽宁省分成5个区域进行研究, 分别是辽西、辽北、辽中、辽东和辽南(图1)。

图1 研究区域及分区示意图

1.2 数据及来源

研究资料为辽宁省52个气象站1961—2015年基本观测资料, 包括平均气温(℃)、日最高气温(℃)、日最低气温(℃)、平均相对湿度(%)、日降水量(mm)、平均风速(m∙s-1)、平均气压(hPa)、日照时数(h), 来自国家气象信息中心。产量资料为辽宁省51个春玉米主产区、县1981—2015年春玉米种植面积(hm2)、单位面积产量(简称单产, t∙hm-2), 来自中国农业科学院农业信息研究所。

1.3 研究方法

1.3.1 水分亏缺距平指数(CWDIa)的构建

本研究采用水分亏缺距平指数(CWDIa)作为农业干旱指数计算基础, 计算公式如下:

CWDI=×CWDI+×CWDI-1+×CWDI-2+×

CWDI-3+×CWDI-4(2)

式中:P为降水量(mm);I为灌溉量(mm); ETcj为实际蒸散量(mm), ETcj=c×ET0计算, ET0为作物可能蒸散量[采用联合国粮农组织(FAO 1998)推荐的Penman-Monteith公式[19]计算],c为作物所处发育阶段的作物系数或多种作物的平均作物系数, 采用FAO数值或国内临近地区通过试验确定的数值,c值确定方法见国标《GBT 32136—2015农业干旱等级》[17]。

1.3.2 基于游程理论和CWDIa的干旱识别以及干旱指数计算

基于春玉米生育期内逐日CWDIa, 根据国标《GBT 32136—2015 农业干旱等级》, CWDIa大于阈值0.4(CWDI0, 水分亏缺距平触发干旱的阈值)则发生干旱。如图2,s为春玉米发育期内一个干旱事件的开始,a为该干旱事件的结束, 大于CWDI0的面积为该干旱事件的干旱强度。逐年总的干旱强度(total)为春玉米发育期内所有干旱事件强度之和, 并将total作为该年春玉米干旱指数。

图2 游程理论示意图

CWDIa: 作物水分亏缺距平指数; CWDI0: 水分亏缺距平触发干旱的阈值;s: 干旱事件开始时间;a: 干旱事件结束时间;: 干旱事件强度;: 干旱事件历时。CWDIa: crop water deficit abnormal index; CWDI0: drought threshold value of crop water deficit abnormal index;s: the beginning time of a drought event;a: the ending time of a drought event;drought intensity of a drought event;: drought duration of a drought event.

1.3.3 极值概率模型

本文选择国内外广泛应用于极端气象事件的统计分布模型——广义极值分布[20]来拟合干旱强度指数序列, 函数如下所示:

式中:=(-)/,、、分别为广义极值分布(generalized extreme value, GEV)函数的形状参数、尺度参数和位置参数, 其中尺度参数>0。GEV的分布范围取决于形状参数, 不同的形状参数产生不同的极值分布类型, 即极值Ⅰ型、Ⅱ型和Ⅲ型。当=0, 为极值Ⅰ型, 也就是耿贝尔分布(Gumbel分布), 即:

式中:=(-)/,、分别为Gumbel分布函数的尺度参数和位置参数。当>0, 为极值Ⅱ型, 也就是Fréchet分布; 当<0时, 为极值Ⅲ型, 即为逆威尔布分布。

不同重现期下干旱强度估算公式如下:

式中:x为不同重现期的干旱强度(return level);为重现期, 如10年、50年。

参数估计方法采用极大似然法(maximum likelihood method), 同时利用Kolmogorov-Smirnov (K-S)检验法比较数据分布形式与已知分布是否一致, 即对模拟拟合优度检验。

1.3.4 产量损失率计算

利用辽宁省春玉米主产区、县1981—2015年产量序列, 采用“趋势-波动”分解模型, 将产量()时间序列分离成技术以及其他因素影响的趋势项(t)、气候项(W)等。由于干旱是影响辽宁省春玉米产量波动的主要因素[21], 因此, 认为剔除趋势和技术产量部分t后, 即为干旱造成的产量损失W。根据各区县春玉米产量变化特点,t采用线性拟合法、多项式拟合法、差值法等进行拟合, 通过<0.05信度检验的拟合方法和数据给予采用, (Yt)即为春玉米干旱产量损失。为消除各区、县间产量水平差异, 采用相对气象产量表示产量损失率。计算方法如下:

=(-t)/t´100% (7)

2 结果与分析

2.1 春玉米干旱指数分布特征

根据1.3.1和1.3.2节方法, 提取出1961—2015年辽宁省各气象站点干旱事件, 并计算春玉米全生育期干旱指数55年的平均值及最大值, 将结果进行反距离加权插值, 结果见图3。

图3a显示, 辽宁省各站点春玉米干旱指数55年平均值由西部向东部递减, 指数最大的区域在辽西, 为10.7~12.3, 最大值出现在义县; 指数次大的区域为辽南, 为6.4~10.2, 最大值出现在旅顺口; 辽北和辽中干旱指数为5.1~10.1, 总体低于辽西和辽南, 高于辽东; 辽东干旱指数为1.8~5.2, 最小值出现在新宾。可见, 1961—2015年春玉米生育期内, 辽西干旱指数较大, 农业干旱总体较重。

图3b显示, 55年间春玉米干旱指数最大值出现在辽西的建昌(1961年), 辽北的法库(2001年)次之, 干旱指数最大值分别为54.9和48.9; 辽西和辽南（除庄河市和大石桥市）干旱指数最大值均在33.2以上; 此外, 辽中的盘山(1965年)、大洼(2014年)等地也曾发生过较为严重的农业干旱, 干旱指数最大值分别为45.0和44.1; 辽东的各区、县干旱指数最大值仍为最小, 为10.8~28.4。

图3 1961—2015年辽宁省春玉米干旱指数平均值(a)和最大值(b)的空间分析

2.2 春玉米产量稳定性分析

考虑到变异系数能够表达序列数据变化的离散程度, 选取产量序列完整的2000—2015年辽宁省各区、县春玉米单产资料, 并计算其变异系数(图4), 用以分析春玉米产量年际间的差异变化和稳定性。从图4可知, 辽西春玉米单产变异系数最大, 除绥中县、兴城市、凌海市、北镇市、黑山县、彰武县外, 其余区、县单产变异系数均在0.30以上, 为单产年际间差异最大、稳定性最差的区域; 辽北除康平县、法库县和调兵山市单产变异系数在0.30~0.45外, 大部为0.15~0.30, 单产年际间差异仅低于辽西; 辽中、辽东和辽南大部变异系数在0.30以下, 其中, 辽东大部<0.15, 为单产年际间差异最小、最稳定的区域。

图4 2000—2015年辽宁省春玉米单产变异系数空间分布

2.3 春玉米干旱产量损失评估模型

通过2.1和2.2节的分析可知, 辽宁省春玉米干旱指数大小与单产年际间差异的分布具有较明显的相似性, 但又非完全相同, 因此, 要具体阐明各地干旱产量损失风险, 需要构建具有针对性的模型进行精细化评估。

2.3.1 模型构建

通过分析春玉米产量损失率和干旱指数的特点, 综合考虑地理位置、地形特点、抗旱能力、社会经济等因素, 构建辽宁省春玉米各主产区、县干旱产量损失评估模型, 见表1。其中, 对于未发生干旱(干旱指数total=0)、产量有损失(<0)的年份, 如1985年、1986年、1992年、1994年、1995年、2010年、2012年等, 某些区、县却是明显减产年, 经分析该年春玉米生长发育期间发生了较为严重的洪涝或低温灾害[22-24]; 或者某年发生了较为严重的干旱, 但却没有导致产量相应减少, 可能是该地区采取了有效的灌溉、人工增雨等抗旱措施, 干旱损失率受到较大的人为影响。以上两种情况下的年份不参与旱灾损失评估模型的构建。

2.3.2 模型参数分析

根据2.3.1节构建的旱灾产量损失模型, 在干旱指数变化相同的情况下, 参数越大, 产量损失率越大, 说明一定程度上是春玉米对干旱敏感性及承受能力的综合体现。为分析参数的分布特征, 对于辽宁省产量缺失年份较多、未构建损失模型的区、县, 本文采取地理位置相邻和行政区域隶属原则, 对参数进行赋值, 其空间分布见图5。从图5可见,在辽宁省的地域分布特征较为明显。总体表现为, 其绝对值辽西最大, 其中朝阳地区又为辽西之最, 而凌海市、北镇市和黑山县相对较小; 辽南和辽北次大, 辽中较小, 辽东除抚顺外, 大部分区、县产量对干旱指数变化无规律性反应, 凤城市、东港市产量变化对干旱表现为正向响应。以干旱指数增加5.0为例, 辽西大部产量减少10.0%~15.0%, 凌源市、喀左县、建昌县、朝阳县、北票市等地减少15.0%~ 20.5%, 北镇市、凌海市、黑山县等地减少5.0%~ 10.0%; 辽南和辽北大部分区、县产量减少10.0%~ 15.0%; 辽中大部减少10.0%以下; 辽东大部以及辽中的台安县、盘山县等地春玉米产量与干旱指数无显著相关性, 而凤城市、东港市降水偏少的年份春玉米产量反而略有增加。

2.3.3 模型检验

表1显示, 辽宁省大部春玉米旱灾产量损失率与干旱指数呈线性关系, 春玉米各主产区、县逐年干旱指数与产量损失率相关系数达显著水平的区域占80%以上, 相关性较好, 可用于玉米干旱产量损失评估。其中辽西、辽北、辽南、辽中大部以及辽东的抚顺县干旱指数与产量损失率呈正相关, 模型均能通过<0.05的显著性检验; 辽中部分(台安县、盘山县)以及辽东大部区县干旱指数与产量损失率相关性不显著, 而辽东的凤城市和东港市干旱指数与产量损失率呈显著负相关(<0.05)。

2.4 不同干旱风险水平的春玉米产量损失

根据1.3.3节方法, 计算各地不同风险水平的干旱指数, 结合2.3.1节干旱产量损失评估模型, 得到不同干旱风险水平下(10年一遇、50年一遇)辽宁省春玉米产量损失分布(图6)。

图6显示, 在10年、50年一遇干旱风险下, 辽宁省春玉米干旱减产率总体均呈现由西向东逐渐递减的特征。10年一遇干旱风险下, 春玉米产量损失最大的区域为辽西中西部(主要为朝阳市及所辖区县), 损失率在50%以上, 辽西其余大部地区和辽南的瓦房店市、辽北的昌图县和法库县产量损失率为30%~50%; 辽西的锦州市、辽北大部、辽中和辽南大部以及辽东的抚顺市产量损失率为10%~30%; 辽中的沈阳市和辽阳市及其所辖区县、大洼县产量损失率在10%以下; 辽东的凤城市和东港市、辽北的西丰县在10年一遇干旱风险下春玉米产量无损失。

表1 辽宁省春玉米主产区、县干旱产量损失评估模型

为产量损失率,为干旱指数; “/”表示模型未通过相关显著性检验。“” is the yield loss rate, “” is the drought index. “/” means no significant correlation.

图5 辽宁省春玉米产量损失模型(y=ax+b)参数a的空间分布

参数的意义见表1。Meaning the parameterwas shown in the table 1.

50年一遇干旱风险下, 辽西西部、辽北大部、辽南的瓦房店以及辽中的海城市、辽东的抚顺市春玉米产量损失率为100%; 辽西其他区、县产量损失率在50%以上; 辽南大部春玉米产量损失率为20%~50%; 辽中产量损失率相差较大, 北部普遍高于南部; 辽东的凤城市和东港市春玉米在50年一遇干旱风险下产量无损失, 而辽北的西丰县产量损失率在20%以下。

3 结论与讨论

本研究以作物水分亏缺距平指数为基础, 根据基于游程理论的干旱识别方法, 构建了春玉米全生育期干旱指数, 并分析了辽宁省春玉米干旱指数及产量稳定性分布特征, 建立了基于干旱指数和产量损失率的干旱产量损失评估模型。本文构建的农业干旱指数, 理论性和可操作性较强, 弥补了由于缺少历史干旱发生程度、持续时间和损失等定量化信息导致在评价春玉米干旱上的局限性。

图6 辽宁省春玉米10年(a)、50年(b)一遇干旱减产率

1961—2015年辽宁省春玉米干旱指数空间分布显示, 干旱指数平均值总体呈由西部向东部递减的经向分布, 最大值区域为辽西, 其次为辽南, 辽东最小; 55年间干旱指数的最大值出现在辽西的建昌(1961年)。总体来看, 辽西地区更易发生严重的农业干旱事件, 与大多数学者的研究结论一致[25-26], 说明构建的干旱指数能够反映干旱的实际情况。

与干旱指数分布相似, 辽宁省春玉米产量最不稳定的区域为辽西, 最稳定的区域为辽东; 辽北、辽中和辽南产量的稳定性总体差于辽东、好于辽西, 个别区、县年际间产量波动较大, 说明辽宁省春玉米产量的波动与干旱程度关系密切[27]。基于此构建的干旱产量损失评估模型也显示, 辽宁省80%以上春玉米主产区、县产量损失率与干旱指数显著相关, 干旱对春玉米产量的影响在空间上差异明显, 对辽西地区影响最大, 其中又以朝阳地区为最, 其次为辽北和辽南, 以上大部区域又为春玉米高产区, 因此也是辽宁省春玉米干旱重点防范区域, 需增加耐旱品种种植, 加强农田水利与灌溉工程, 发展农业节水灌溉技术, 提高区域整体抗旱能力, 并及时监测预报农业干旱, 适时开展人工增雨, 做好干旱的防灾减灾工作;对辽东湾附近以及中部平原地区影响较小, 辽中的台安县和盘山县以及辽东大部区、县春玉米产量损失与干旱指数无显著相关性, 主要由于台安县和盘山县境内有10余条水系, 灌溉条件极好, 干旱年均可得到不同程度灌溉; 辽东具有特殊的地形及气候特征, 尤其是凤城市和东港市, 降水偏少使春玉米产量不减反增, 这可能与该区域自然降水完全能够满足春玉米生长发育的需要, 且是洪涝、风雹灾害高发区有关[28-29], 可通过提高种植和管理水平、加强自然灾害防御的手段增加产量。

在10年、50年一遇干旱风险下, 辽宁省春玉米干旱减产率总体呈现由西向东逐渐递减的特征。50年一遇干旱风险下, 辽西和辽北除凌海市和西丰县外、辽南的瓦房店和盖州市以及辽中的海城市和辽中县、辽东的抚顺市等地, 产量损失率均在50%以上、甚至绝收。此空间分布特征主要与气候、地形条件等综合因素有关, 与许多学者的研究结论大体一致[28], 但在具体细节上存在一定差别。除评价指标不同外, 本文直接以区、县为研究单位, 而大多数研究是以气象站点为单位再插值到整个区域, 两种研究方法和精度的不同也是导致研究结果存在差异的重要原因。

本文在构建干旱产量损失评估模型时, 一是由于产量数据的限制, 仅针对2015年及以前的历史产量序列较完整的51个春玉米主产区、县进行了建模; 二是虽然已排除典型洪涝、低温冷害等影响较大的年份, 也考虑了重旱年份灌溉的影响, 但未考虑春玉米生长发育期间发生时间短、影响较轻的其他自然灾害、病虫害等对产量的负面影响。今后如能得到完整的历史产量资料以及分类灾损和各区县灌溉数据, 将使春玉米干旱产量损失评估的结果更加精确。

[1] 曹永强, 张兰霞, 张岳军, 等. 基于CI指数的辽宁省气象干旱特征分析[J]. 资源科学, 2012, 34(2): 265–272 CAO Y Q, ZHANG L X, ZHANG Y J, et al. Analysis of drought characteristics of Liaoning Province based on the CI index[J]. Resources Science, 2012, 34(2): 265–272

[2] 姚玉璧, 李耀辉, 石界, 等. 基于GIS的石羊河流域干旱灾害风险评估与区划[J]. 干旱地区农业研究, 2014, 32(2): 21–28 YAO Y B, LI Y H, SHI J, et al. Assessment and division of drought hazard risk in Shiyang River basin based on GIS[J]. Agricultural Research in the Arid Areas, 2014, 32(2): 21–28

[3] 李叶蓓, 陶洪斌, 王若男, 等. 干旱对玉米穗发育及产量的影响[J]. 中国生态农业学报, 2015, 23(4): 383–391 LI Y B, TAO H B, WANG R N, et al. Effect of drought on ear development and yield of maize[J]. Chinese Journal of Eco-Agriculture, 2015, 23(4): 383–391

[4] 胡雪琼, 吉文娟, 张茂松, 等. 云南省冬小麦干旱灾损风险区划[J]. 大气科学学报, 2011, 34(3): 356–362 HU X Q, JI W J, ZHANG M S, et al. Drought risk regionalization for winter wheat in Yunnan Province[J]. Transactions of Atmospheric Sciences, 2011, 34(3): 356–362

[5] ZHANG J Q. Risk assessment of drought disaster in the maize-growing region of Songliao Plain, China[J]. Agriculture, Ecosystems & Environment, 2004, 102(2): 133–153

[6] 李明, 胡炜霞, 张莲芝, 等. 基于SPEI的东北地区气象干旱风险分析[J]. 干旱区资源与环境, 2018, 32(7): 134–139 LI M, HU W X, ZHANG L Z, et al. Risk analysis of meteorological drought in Northeast China based on standardized precipitation evapotranspiration index[J]. Journal of Arid Land Resources and Environment, 2018, 32(7): 134–139

[7] 沈国强, 郑海峰, 雷振锋. 基于SPEI指数的1961—2014年东北地区气象干旱时空特征研究[J]. 生态学报, 2017, 37(17): 5882–5893 SHEN G Q, ZHENG H F, LEI Z F. Spatiotemporal analysis of meteorological drought (1961-2014) in Northeast China using a standardized precipitation evapotranspiration index[J]. Acta Ecologica Sinica, 2017, 37(17): 5882–5893

[8] 刘荣花, 朱自玺, 方文松, 等. 华北平原冬小麦干旱区划初探[J]. 自然灾害学报, 2003, 12(1): 140–144 LIU R H , ZHU Z X, FANG W S, et al. Preliminary study on drought zonation of winter wheat in North China Plain[J]. Journal of Natural Disasters, 2003, 12(1): 140–144

[9] 董朝阳, 刘志娟, 杨晓光. 北方地区不同等级干旱对春玉米产量影响[J]. 农业工程学报, 2015, 31(11): 157–164 DONG C Y, LIU Z J, YANG X G. Effects of different grade drought on grain yield of spring maize in northern China[J]. Transactions of the CSAE, 2015, 31(11): 157–164

[10] NUTINI F, BOSCHETTI M, CANDIANI G, et al. Evaporative fraction as an indicator of moisture condition and water stress status in semi-arid rangeland ecosystems[J]. Remote Sensing, 2014, 6(7): 6300–6323

[11] 宋艳玲, 王建林, 田靳峰, 等. 气象干旱指数在东北春玉米干旱监测中的改进[J]. 应用气象学报, 2019, 30(1): 25–34 SONG Y L, WANG J L, TIAN J F, et al. The spring maize drought index in Northeast China based on meteorological drought index[J]. Journal of Applied Meteorological Science, 2019, 30(1): 25–34

[12] 卫捷, 陶诗言, 张庆云. Palmer干旱指数在华北干旱分析中的应用[J]. 地理学报, 2003, 58(S1): 91–99 WEI J, TAO S Y, ZHANG Q Y. Analysis of drought in northern China based on the Palmer severity drought index[J]. Acta Geographica Sinica, 2003, 58(S1): 91–99

[13] 张调风, 张勃, 刘秀丽, 等. 基于CI指数的甘肃省黄土高原地区气象干旱的变化趋势分析[J]. 冰川冻土, 2012, 34(5): 1076–1083 ZHANG T F, ZHANG B, LIU X L, et al. Trend analysis of the variation of meteorological drought in Loess Plateau of Gansu Province based on comprehensive meteorological drought index[J]. Journal of Glaciology and Geocryology, 2012, 34(5): 1076–1083

[14] MISHRA A K, SINGH V P. A review of drought concepts[J]. Journal of Hydrology, 2010, 391: 202–216

[15] 王晓峰, 张园, 冯晓明, 等. 基于游程理论和Copula函数的干旱特征分析及应用[J]. 农业工程学报, 2017, 33(10): 206–214 WANG X F, ZHANG Y, FENG X M, et al. Analysis and application of drought characteristics based on run theory and Copula function[J]. Transactions of the CSAE, 2017, 33(10): 206–214

[16] 左冬冬, 侯威, 颜鹏程, 等. 基于游程理论和两变量联合分布的中国西南地区干旱特征研究[J]. 物理学报, 2014, 63(23): 53–64 ZUO D D, HOU W, YAN P C, et al. Research on drought in Southwest China based on the theory of run and two-dimensional joint distribution theory[J]. Acta Physica Sinica, 2014, 63(23): 53–64

[17] 吕厚荃, 张玉书, 李茂松, 等. 中华人民共和国国家标准(GB/T 32136—2015). 农业干旱等级[S]. 北京: 中国标准出版社, 2015 LYU H Q, ZHANG Y S, LI M S, et al. National Standards of People’s Republic of China (GB/T 32136–2015). Grade of Agricultural Drought[S]. Beijing: China Standard Press, 2015

[18] 郭进修, 王秀荣, 王维国. 中国气象地理区划手册[M]. 北京: 气象出版社, 2006 GUO J X, WANG X R, WANG W G. China Meteorological Geography Zoning Manual[M]. Beijing: China Meteorological Press, 2006

[19] ALLEN R G, PEREIRA L S, RAES D, et al Crop Evapotranspiration: Guidelines for Computing Crop Water Requirements[M]. Rome: FAO Irrigation and Drainage Paper 56, 1998: 15–86

[20] COLES S. Basics of statistical modeling[M]//An Introduction to Statistical Modeling of Extreme Values. London: Springer London, 2001: 18–44

[21] 单琨, 刘布春, 刘园, 等. 基于自然灾害系统理论的辽宁省玉米干旱风险分析[J]. 农业工程学报, 2012, 28(8): 186–194 SHAN K, LIU B C, LIU Y, et al. Analysis on drought risk of maize based on natural disaster system theory in Liaoning Province[J]. Transactions of the CSAE, 2012, 28(8): 186–194

[22] 江和文, 郭婷婷, 包颖, 等. 辽宁省近50年旱涝灾害的时空特征分析[J]. 水土保持研究, 2012, 19(2): 29–33 JIANG H W, GUO T T, BAO Y, et al. Spatiotemporal distribution of drought and flood disasters over the last 50 years of Liaoning Province[J]. Research of Soil and Water Conservation, 2012, 19(2): 29–33

[23] 陈鹏狮, 张玉书, 纪瑞鹏, 等. 1961—2010年辽宁省玉米低温冷害时空分布特征[J]. 自然灾害学报, 2017, 26(4): 164–171CHEN P S, ZHANG Y S, JI R P, et al. Spatio-temporal distribution characteristics of maize chilling damage in Liaoning Province during 1961-2010[J]. Journal of Natural Disasters, 2017, 26(4): 164–171

[24] 杨若子, 周广胜. 1961—2013年东北三省玉米低温冷害强度的时空分布特征[J]. 生态学报, 2016, 36(14): 4386–4394 YANG R Z, ZHOU G S. Spatio-temporal distribution of maize chilling damage intensity in the three provinces of Northeast China during 1961 to 2013[J]. Acta Ecologica Sinica, 2016, 36(14): 4386–4394

[25] 张淑杰, 张玉书, 李广霞, 等. 辽宁省旱灾的分布特征及其成因分析[J]. 中国农学通报, 2013, 29(5): 199–203 ZHANG S J, ZHANG Y S, LI G X, et al. Distribution characteristics and main causes of drought in Liaoning Province[J]. Chinese Agricultural Science Bulletin, 2013, 29(5): 199–203

[26] 杨晓晨, 明博, 陶洪斌, 等. 中国东北春玉米区干旱时空分布特征及其对产量的影响[J]. 中国生态农业学报, 2015, 23(6): 758–767 YANG X C, MING B, TAO H B, et al. Spatial distribution characteristics and impact on spring maize yield of drought in Northeast China[J]. Chinese Journal of Eco-Agriculture, 2015, 23(6): 758–767

[27] 张淑杰, 周广胜, 张玉书, 等. 辽宁地区玉米耗水量与产量的关系[J]. 干旱地区农业研究, 2018, 36(1): 244–250 ZHANG S J, ZHOU G S, ZHANG Y S, et al. The relationship between yield and water consumption of maize in Liaoning Province[J]. Agricultural Research in the Arid Areas, 2018, 36(1): 244–250

[28] 高晓容, 王春乙, 张继权, 等. 东北地区玉米主要气象灾害风险评价与区划[J]. 中国农业科学, 2014, 47(24): 4805–4820 GAO X R, WANG C Y, ZHANG J Q, et al. Risk assessment and zoning of the main meteorological disasters for maize in Northeast China[J]. Scientia Agricultura Sinica, 2014, 47(24): 4805–4820

[29] 魏新光, 王铁良, 李波, 等. 辽宁省玉米地水分盈亏时空分布特征及灌溉模式分区研究[J]. 农业工程学报, 2018, 34(23): 119–126 WEI X G, WANG T L, LI B, et al. Temporal and spatial distribution characteristics of maize water surplus deficit and irrigation mode partition in Liaoning Province[J]. Transactions of the CSAE, 2018, 34(23): 119–126

Application of agricultural drought index based on Run Theory for the assessment of yield loss in spring maize owing to drought in Liaoning Province*

TAN Fangying, HE Liang, LYU Houquan, SONG Yingbo, CHENG Lu, HOU Yingyu**

(National Meteorological Center, Beijing 100081, China)

Drought is one of the main factors restricting agricultural development and food security in China. Identifying agricultural drought accurately and assessing the impact of drought on crop yield in detail is the core of drought research. The results of this research can provide a theoretical basis for drought disaster prevention and mitigation. The water deficit anomaly index of spring maize across the entire growth period was calculated based on daily weather data from 52 meteorological stations in Liaoning Province from 1961 to 2015, and a new drought index for spring maize over the entire growth period was constructed by using Run Theory. Drought yield loss assessment models at the county or city scales were further constructed based on the drought index and yield loss rate from 1981 to 2015. Finally, key prevention areas for spring maize drought in Liaoning Province were identified. The drought index for spring maize in Liaoning Province exhibited a decreasing distribution from west to east. The likelihood of severe drought events was higher and the interannual variability of spring maize yield was largest in western Liaoning. The rate of yield loss of spring maize had a significant positive linear correlation with the drought index in more than 80% of the counties in the spring maize main producing area, except for most parts of eastern Liaoning and some counties in central Liaoning. Moreover, significant spatial differences were observed in the yield loss of spring maize. With the same drought index, the region with the largest yield loss rate of Liaoning Province was Chaoyang, followed by the southern and northern parts of Liaoning; the Liaodong Bay and the central plain had the lowest yield loss rate. However, in Fengcheng City and Donggang City, where spring maize growth was not restricted by precipitation, the yield of corn was higher in the years with less precipitation. The results of the present study indicated that the government should pay more attention to drought-sensitive zones such as the Chaoyang District, and northern and southern Liaoning west, and north of the Qianshan Mountains. To increase and stabilize spring maize yield in the above areas, it is necessary to increase drought-tolerant varieties, develop water-saving irrigation, and improve water use efficiency.Most counties in eastern Liaoning and some counties in central Liaoning with abundant rainfall or good irrigation conditions can increase spring maize production by improving planting and management methods. Additional attention should be paid to the prevention of damage caused by floods and cold, especially in Fengcheng and Donggang.

Agricultural drought index; Yield loss assessment; Run Theory; Spring maize; Liaoning Province

S166

10.13930/j.cnki.cjea.190675

* 国家重点研发计划项目(2018YFC1507802, 2017YFC1502402)、公益性行业(气象)科研专项(GYHY201306038, GYHY201506001)和国家自然科学基金项目(41705095)资助

侯英雨, 主要从事农业气象相关研究。E-mail: yyhou@cma.gov.cn

谭方颖, 主要从事农业气象相关研究。E-mail: tanfangying0803@163.com

2019-09-17

2019-11-07

* This study was supported by the National Key Research and Development Project of China (2018YFC1507802, 2017YFC1502402), the Special Scientific Research Fund of Meteorology in the Public Welfare Profession of China (GYHY201306038, GYHY201506001) and the National Natural Science Foundation of China (41705095).

, E-mail: yyhou@cma.gov.cn

Nov. 7, 2019

Sep. 17, 2019;

谭方颖, 何亮, 吕厚荃, 宋迎波, 程路, 侯英雨. 基于游程理论的农业干旱指数在辽宁省春玉米旱灾损失评估中的应用[J]. 中国生态农业学报(中英文), 2020, 28(2): 191-199

TAN F Y, HE L, LYU H Q, SONG Y B, CHENG L, HOU Y Y. Application of agricultural drought index based on Run Theory for the assessment of yield loss in spring maize owing to drought in Liaoning Province[J]. Chinese Journal of Eco-Agriculture, 2020, 28(2): 191-199