胡金鹏，金 明，史培军
(1.北京师范大学 应急管理部-教育部减灾与应急管理研究院，北京 100875；
2.北京师范大学 地理科学学部，北京 100875；
3.兰州大学 资源环境学院，甘肃 兰州 730000；
4.青海师范大学 高原科学与可持续发展研究院，青海 西宁 810008)
随着全球增温,人类活动的加剧，干旱的频率、强度和影响范围进一步增加[1-4]，给水资源、农业、植被、环境和人类生活带来了巨大的威胁[5]。我国是一个干旱灾害频发的国家[6-8]，因干旱造成的粮食减产非常严重[9-10]，因此定量评估干旱风险损失对于保护国家粮食安全，开展灾害风险防范[11]至关重要。
近年来，对农作物干旱损失的研究更多偏向于干旱指标与统计数据间的相关性分析[12-15]，这些研究不能完全反映农业干旱损失形成的物理机制，同时由于大部分地区均有灌溉的习惯，通过干旱指标来表征区域范围内的干旱损失并不一定能够反映实际的农业干旱灾情[16]。将作物生长模型应用于农业脆弱性和风险评估研究是一种新兴趋势[17-20]，这种方法可以详细描述灾害形成因素的强度与相应的承灾体损失之间的定量关系。本文选取雨养农业区为研究对象，其农业生产对于气候变化的响应比较敏感，受人工灌溉影响较小，通过采用作物模型与干旱指标相结合的方法构建脆弱性曲线，可以更好的定量研究干旱对农作物产量带来的潜在影响[21]。
甘肃省的陇东地区属于典型的雨养农业区，冬小麦作为陇东地区的主要农作物，由于其生育期较长，受干旱这种缓发性灾害影响较大。DSSAT作为应用最广泛、最有效的作物生长模拟软件之一，可以在日尺度上模拟作物生长过程，因此，本研究采用DSSAT作物模型与SPEI干旱指标相结合的方式构建干旱脆弱性曲线，探究陇东干旱对冬小麦产量的潜在影响，进而实行补偿灌溉，为雨养农业区的抗旱减灾、风险防范提供思路。
1.1 研究区概况
陇东地区地处甘肃省东部,属黄河中游的内陆地区，是半湿润气候向干旱半干旱气候的过渡地带，属于典型的雨养农业区，对气候变化比较敏感,同时也是生态环境比较脆弱的地区[22-23],近年来干旱发生频繁，对区域的农业生产影响显著[24-25]。本文的研究范围为行政区划中所包含的陇东地区15个县区。图1为研究区地形和气象站点分布。
图1 研究区域及站点分布示意图(基于自然资源部标准地图服务网站审图号为GS(2019)1822号的标准地图制作，底图无修改)
1.2 数据来源
本研究所用的数据主要包括：①气象数据由中国气象数据网(http://data.cma.cn)提供，涵盖的数据类型包括1961—2020年逐日和逐月的降水量、日照时间、平均气温、平均最高气温、平均最低气温等。②模型所用土壤数据来自中国土壤数据库(http://vdb3.soil.csdb.cn)和世界土壤数据库(Harmonized World Soil Database)，数据类型包括土壤各层的机械组成成分、有机碳含量、阳离子交换量、全氮含量等。③模型所用田间管理数据由国家气象信息中心提供，时间序列为2010—2012 年，包括作物品种、播种量和播种日期、收获产量和收获日期、生育期各阶段特征、施肥量和施肥日期以及耕种等信息。本研究选用的冬小麦品种为环县站点2010—2012 年田间监测的“环冬1号”。
1.3 研究方法
1.3.1 标准化降水蒸散发指数(SPEI)
SPEI指数是一种新理想化的干旱趋势分析指标，通过计算降水量与潜在蒸散量之间的关系来反应气候的干湿状况。SPEI算法详见文献[26-27]，根据《气象干旱等级：GB/T 20481-2017》并参考前人研究经验，本文基于SPEI值将干旱等级划分为6个等级(表1)。
表1 SPEI指数干旱等级划分
1.3.2 干旱识别
干旱事件的过程主要包括干旱的严重程度和持续时间。目前，游程理论是识别干旱事件最常用的方法，以游程理论为工具，从计算的SPEI指标序列中识别出干旱历时和干旱强度两个特征变量[28-29]。具体步骤如下：
(1)确定三个阈值X0，X1，X2。根据研究经验，选取X0=0，X1=-0.3，X2=-0.5。
(2)将干旱指数小于X1的事件定义为潜在干旱事件，如果潜在干旱事件之间相互连续，则将连续的事件合并为一次事件。如图2识别出5次潜在干旱事件a，b，C1，C2，d。
图2 干旱识别及干旱特征变量的确定
(3)在定义的潜在干旱事件基础上，对于历时为一个月的潜在干旱事件，若其SPEI≥X2则该潜在干旱事件将被忽略，若其SPEI 1.3.3 DSSAT模型校验分析方法 检验作物模型的适用性时，采用标准均方根误差(ENRMS)评价模型的校验效果[30]。当标准均方根误差ENRMS的值小于10%时，可以认为模型模拟的效果很好； (1) (2) 式中：Mi代表实际观测值，Si代表模拟值，M代表实际观测平均值，n为样本数。 1.3.4 联合分布函数Copula 由于干旱特征变量之间并不独立，对于联合分布的建模会变得十分困难，而Copula函数可以将变量的随机性和耦合性很好的分离开来，故本研究采用Copula函数对干旱特征之间的相关性进行建模。常用的Copula联合分布函数类型有3种: 二次型Copula函数、椭圆型Copula函数、阿基米德型Copula函数，其中含有一个参数的阿基米德型Copula函数计算简单，在干旱特征研究方面应用较多[31]，因此本研究选用阿基米德型Copula函数作为干旱特征变量的联合分布函数(表2)。 表2 Copulas函数及参数取值范围 本文采用KS检验(Kolmogorov-Simirnov)对选用的联合分布函数进行拟合相似性评价，采用最短欧氏距离(d2)最小准则对Copula函数进行拟合优度比较。经过检验对比选用Frank Copula函数作为本次研究干旱特征的联合分布函数。 1.3.5 联合重现期 对于灾害事件来说，其重现期T是指事件发生的平均长度，即超过概率的倒数。根据干旱具有多年持续性的特点，传统的基于单变量的重现期计算公式为: (3) (4) 联合分布的二维重现期与单变量重现期不同，其计算方式如下[32-33]： (5) (6) 式中:TD为干旱历时的重现期； 1.3.6 干旱损失计算方法 通过使用DSSAT作物生长模型计算作物在理想气象条件下与干旱条件下的产量差值作为干旱减产的计算依据[34]。计算公式为： (7) 式中：L为干旱产量损失率； 本文将陇东地区1961—2020年气象数据的多年平均数据作为理想背景，同时通过灌溉措施使其满足冬小麦需水要求； 本研究将“环冬1号”小麦品种进行本地化参数校正，把2010年作为小麦品种参数的校正年份。品种参数校验采用DSSAT模型中的最大似然估计模块(GLUE)进行，把模拟作物冬小麦的产量和生育期(播种至开花时间和播种至成熟时间)与实测结果采用均方根误差(ENRMS)进行差异性分析，ENRMS越小表示品种参数越符合实际，选用差异性最小的作为实验品种参数方案，然后将2011、2012年的田间观测资料作为验证资料，表3为模型对生育期(播种至开花时间、播种至成熟时间)和产量的模拟值与观测值的ENRMS，由表 3 可知校正与验证误差基本都满足在 10% 左右，说明模型模拟精度良好，可以使用该品种参数进行研究。 表3 模型检验情况 根据相关研究[35]与国家气象科学数据中心整理的《中国农作物生长发育和农田土壤湿度旬值数据集》进行统计可知，陇东地区冬小麦生长发育期为9月下旬至次年6月下旬，因此本次研究以种植年9月至次年6月作为冬小麦的生长发育期。为了定量分析区域冬小麦生育期的干旱特征，基于游程理论，计算识别了1961—2020年冬小麦生育期的干旱特征变量(干旱历时、干旱强度)(图3)，结果表明1961—2020年陇东地区冬小麦生育期内干旱历时波动性较大，但变化的趋向性不明显，最高历时可达7个月，最低为1个月，平均干旱历时为4个月，占据整个生育期的三分之一； 本文采用干旱历时和干旱强度双变量联合分布的方式来反映干旱胁迫对冬小麦产量的潜在影响，根据识别出的干旱特征(干旱历时、干旱强度)，利用Copula联合分布函数计算出干旱事件的联合重现期，而联合重现期的大小，表明干旱历时和干旱强度耦合达到的年遇型大小，能够反映干旱影响的程度，联合重现期越高，意味着所经历的干旱历时越长或者干旱强度越强。从图4的结果可以发现，大部分干旱事件的联合重现期较低，多为3年一遇的干旱事件，这表明近60年陇东地区冬小麦生育期所经历的干旱事件大多表现为短历时或者低强度的特征，长历时和高强度同时发生的干旱事件较少，在1990—2010年期间，干旱事件的联合重现期增大，最高能达到15年一遇，该时期出现的长历时和高强度的干旱事件相对较多。 图3 干旱历时和干旱强度变化情况 图4 干旱的联合重现期 表4 模拟灌溉方案 为了探索灌溉在降低冬小麦干旱损失风险中的作用，输入不同的灌溉方案来模拟不同的灌溉行为。在不考虑冬小麦生育期出现死亡的情况下，将作物需水量作为模型灌溉情景设计中的主要指标，对于降水量高于作物需水量的时期不给予灌溉，处于降水量低于作物需水量时期给予灌溉，灌溉水资源的总量取决于对应时期降水量与作物需水量的差值。针对研究区域，计算理想背景下多年月平均降水量，根据FAO推荐的作物系数法与Penman公式[36]计算冬小麦在生育期内的需水量,结果如图5所示。 图5 冬小麦月平均需水量与月平均降雨量的变化情况 从总量上来看，生育期内总降水量为371 mm，总需水量为521 mm，灌溉水资源总量为两者之间的差值150 mm； 根据分析，本次实验的灌溉时期为种植年12月至次年5月。模拟不同的灌溉水平，需要通过控制冬小麦生育阶段不同灌溉定额来实现。根据计算出的灌水量，进行了五种补水方案(四种干旱背景灌溉处理 T0-T3 和一种理想背景灌溉处理 CK)。在每个模拟灌溉方案中，在需要灌溉的月份中，每两周灌溉一次水，每次灌溉量为该月总设计灌溉定额的二分之一。对于理想背景下，以冬小麦满足需水量为标准进行灌溉； 农作物的产量变化可以直接反映干旱带来的影响，更高的作物产量损失意味着更强的干旱影响。干旱脆弱性曲线(图6)可以反应不同灌溉水平下干旱频率与相应冬小麦损失之间的关系，本文利用SPEI指数与Copula函数识别计算出不同干旱事件的重现期，通过结合作物生长模型计算各干旱重现期下的平均损失情况。根据分析发现干旱重现期与干旱损失之间符合对数函数关系，函数拟合的相关系数(R2)高于0.65，通过95%显著性检验(p<0.05)，表明干旱频率与干旱损失之间的脆弱性曲线能够定量描述不同灌溉水平下干旱对冬小麦产量的潜在影响。具体函数关系式为： L=a-b×ln(T+c)。 (8) 式中：L表示损失率； 图6 不同灌溉水平下的干旱脆弱性曲线 结果显示，在正常无灌溉措施下，冬小麦多年平均损失率为31%； 表5 不同灌溉水平下干旱年遇型损失(%) 本研究通过将作物模型与农业干旱损失风险曲线结合的方式定量描述陇东地区冬小麦 干旱频率、干旱损失和灌溉水平之间的关系，定量评估干旱对冬小麦造成的潜在影响。通过研究得到以下结论： (1)近60年陇东地区冬小麦生育期所经历的干旱事件大多表现为短历时、低强度的特征，平均干旱历时为4个月，月平均强度达到中度干旱水平，而长历时和高强度同时发生的干旱事件相对较少，多集中在1990—2010年期间。 (2)陇东地区冬小麦从种植年的12月到次年5月一直处于水分亏缺状态，生育期内总降水量为371 mm，总需水量为521 mm，水分亏缺总量达到150 mm。从冬小麦干旱损失模拟的结果来看，在不考虑人工灌溉的情景下，冬小麦多年平均干旱损失接近三分之一，严重时可达到减产七成以上。2年一遇、5年一遇、10年一遇的旱灾损失率分别达到28%,36%,42%。 (3)从灌溉水生产效益上来看，随着灌溉水平的提高，模拟灌溉的水生产效益也在不断提高，冬小麦多年平均损失率从自然无灌溉条件下的31%，逐渐下降到26%、19%、9%，在满足相等灌水增量的前提下，满足75%需水量灌溉水平所带来的水生产效益最大，损失下降速率最快，平均1 mm水资源可以降低0.2%的干旱损失。 以上结论是在非极端平稳环境下进行的研究，干旱频率—产量损失曲线图属于一种内插线性拟合，能准确反应已识别干旱特征区间内的损失情况，由于受样本量的限制，对于发生频率低、严重程度高、持续时间长的极端干旱事件，这种方式不能准确反映出内插范围外的实际情况。后续可以使用随机模拟的方法对干旱事件进行仿真,扩大样本量，结合脆弱性模型,得到更全面的干旱损失概率曲线,完成干旱风险评估。在实际情况中，冬小麦的产量除了受气候影响外，还会受当地灌溉水平与种植品种的不同而不同，本文只进行了统一品种，统一灌溉水平下的模拟研究，要想得到更为准确的模拟产量，还需要考虑实地情况，逐年设置更为准确的情景进行模拟。除了将干旱损失降低作为主要目标之外，还需考虑灌溉水生产效益和水利用效率，进一步优化灌溉方案，提高区域的水资源利用效率。在全球逐渐变暖的背景下，可以考虑将未来气候变化情景与干旱脆弱性曲线相结合，定量评估干旱风险，为区域的综合风险防范提供指导。
ENRMS的结果在10%～20%之间时，可以认为模型模拟效果是良好的；
当ENRMS的结果在20%～30%之间时，认为模型模拟效果一般；
而当ENRMS>30%时，则认为模型模拟效果差，评价指标的计算公式为：
Ts为干旱强度的重现期；
E(L)为干旱事件发生间隔的期望值,等于干旱事件序列的长度除以干旱事件发生的次数。FD(d)为干旱历时D的边缘分布；
FS(s)为干旱强度S的边缘分布函数；
F(d,s)为干旱历时与干旱强度的联合分布函数。
YN为理想背景下的产量；
YD为干旱背景下的产量。
将实际气象数据作为干旱背景，其中非干旱历时的月份采用多年平均气象数据替换。2.1 DSSAT作物模型调试与验证
2.2 干旱特征分析
而干旱强度表现出一定的上升趋势，逐月SPEI值达到-1平均水平，属于中度干旱，最高月尺度SPEI值可达-1.54，属于重度干旱。
2.3 灌溉方案
从时间变化上来看，种植年9—11月期间和次年6月期间，区域的降水量比需水量高，无需采用灌溉措施，在12月到次年5月期间，降水量要比需水量低，冬小麦表现出缺水的状态，需要采用补灌措施进行弥补，具体各月灌溉总量分别为7 mm、5 mm、7 mm、61 mm、83 mm、94 mm。
干旱背景下,则分别取冬小麦对应平均干旱气候环境下需水量的不同比例作为灌溉量,为定量分析灌溉效益，灌溉量分别为 0%、25%、50%、75%比例的需水量(表4)。2.4 冬小麦干旱脆弱性曲线
T表示干旱重现期；
参数b表示干旱损失的变化速率。
当灌溉水资源总量为37.5 mm，满足25%需水量条件时，冬小麦多年平均损失率下降到26%，平均每1 mm灌溉量可以降低0.13%的干旱损失；
当灌溉水资源总量为75 mm，满足50%需水量条件时，冬小麦多年平均损失率为19%，平均每1 mm灌溉量可以降低0.16%的干旱损失；
当灌溉水资源总量为112.5 mm，满足75%需水量条件时，冬小麦多年平均损失率为9%，平均每1 mm灌溉量可以降低0.2%的干旱损失。不同灌溉水平下干旱频率与相应冬小麦损失之间的对应关系如表5。
3.1 结论
3.2 讨论