黄土高原北部地区榆林市旱涝急转特征及其对作物产量的影响*

宋惠敏，张雪艳，马 欣

(1.中国农业科学院 农业环境与可持续发展研究所，北京 100081；2.中国科学院地理科学与资源研究所 陆地表层格局与模拟院重点实验室，北京 100101)

旱涝急转事件发生规律及其灾损评估具有重要意义。国内外关于旱灾或涝灾发生规律的研究很多，对旱涝急转事件性规律的研究较少[1-2]。近年，以2011年春夏季长江中下游地区严重的旱涝急转事件为契机，对长江中下游地区旱涝急转特征[3]和成因[4]进行了探索，提出大气环流异常为主的归因结果。在淮河流域，程智等[5]揭示了淮河流域旱涝急转的气候特征，唐明等[6]和张效武等[7]提出了安徽省旱涝急转发生规律和淮北地区旱涝急转的应对措施建议。孙小婷等[8]还对我国西南地区的旱涝急转事件长周期规律进行了分析。在旱涝急转事件的评估方法方面，主要有基于气候因素中降水量设计的长周期旱涝急转指数(LDFAI)[9]、标准化降水指数(SPI)[5]、旬降水距平百分率(Pa)[10]、标准降水蒸发指数(SPEI)[11]，或采用降水距平和一组降水阈值[12]等对区域旱涝急转事件进行分析。其中，LDFAI的绝对值可以反映旱涝急转事件的强度，Pa能反映作物生育期内的旱涝情况[13-14]。然而，当前对旱涝急转事件发生规律的研究主要集中在长流流域、淮河流域和西南地区，对我国北方地区的旱涝急转事件规律关注较少。

水旱灾害是农业生产最主要的挑战之一，每年因干旱和洪涝灾害造成的农业减产在10%以上，在极端的条件下可以接近20%[15]。在我国北方地区，5—6月是玉米、大豆和薯类的播种期、出苗期和三叶期，作物需水量较少，7—8月是作物的开花期和成熟期，作物需水量较高，旱涝急转事件的发生会对作物产量产生一定影响[16-17]。研究表明，旱涝急转对南方水稻产量具有正/负补偿效应[18]，水稻生育后期重旱下的旱涝急转会导致作物严重减产[19-20]。但对于北方地区，特别是生态脆弱地区的农业生产，旱涝急转给作物产量带来的影响仍未进行深入的研究。

黄土高原北部的榆林市位于陕西省北部，水土流失严重、农业用地破碎、多种作物混种，是我国北方地区典型的农业生态环境脆弱区[21]，夏季降水的异常会给当地农业生产带来较大影响[22]，严重的旱涝急转灾害会使农业生产遭受巨大损失。历史上榆林市曾多次发生旱涝急转事件[23-25]，但发生规律尚不清楚[14,26-28]、对农业的影响尚不明确[28-29]。鉴于此，本研究选择榆林市作为北方农业生态环境脆弱区的代表，开展旱涝急转事件发生规律和对农业影响评估的研究，以期为北方地区生态脆弱区应对旱涝急转事件提供参考。

1 材料与方法

1.1 研究区概况

榆林市位于陕西省北部(36°57′～39°35′N，107°28′～111°15′E)，是典型的温带季风气候，多年平均降水量约为429 mm，受季风影响，降水主要集中在夏季，常以暴雨形式出现[30]。榆林市主要作物有玉米、大豆、薯类等，2018年榆林市玉米总产量占榆林市粮食总产量的56.52%，大豆在榆林市粮食作物总产量中排名第三，榆林市马铃薯占全陕西省马铃薯播种面积的50%以上[31]，这三类作物在榆林市的农业发展中有着至关重要的地位。

图1 研究区及气象站点分布示意图

1.2 资料来源与处理

榆林市降水量数据来自中国气象数据共享网(http://data.cma.cn/)，研究区内站点包括府谷、神木、佳县、米脂、子洲、绥德、吴堡、清涧、榆林、横山、靖边和定边。对原始数据进行质量整理和补差，将1981—2018年作物生育期(5—8月)逐日降水量整理成逐月降水量，对1981—2018年的主要作物玉米、大豆和薯类的产量进行标准化处理。

榆林市干旱、洪涝灾害历史数据来源于1981—2018年的《神木县志》[32]《子洲县志》[24]《府谷县志》[33]《佳县志》[34]《米脂县志》[25]《靖边县志》[35]《吴堡县志》[36]《清涧县志》[23]《定边县志》[37]《绥德县志》[38]《中国气象灾害大典·陕西卷》[39]。

1.3 方法

本研究中长周期变化的时间尺度为2个月，“旱转涝”事件通常指5—6月发生较长时间连续无降水的干旱事件，7—8月突发强降水或连阴雨，发生水涝事件；反之，则为“涝转旱”[9]。榆林市典型作物生育期主要在5—8月，故对每个气象站点5—6月和7—8月的降水量标准化后计算长周期旱涝急转指数(LDFAI)：

LDFAI=(R78-R56)×(|R56|+|R78|)×1.8-|R56+R78|。

(1)

式中：R56为5—6月标准化降水量；R78为7—8月标准化降水量；(|R56-R78|)是旱涝急转强度项；(|R56|+|R78|)是旱涝强度项；1.8-|R56+R78|为权重系数，其作用是增加旱涝急转事件所占权重，降低全旱或全涝事件权重。采用LDFAI值的绝对值(ILDFAI)表达旱涝急转事件的强度，ILDFAI值越大，说明旱涝急转事件越严重[40]。根据榆林市气候及农业特征，使用降水距平百分率(Pa)辅助判断旱涝情况，定义Pa>25%为偏涝，定义Pa<-25%为偏旱[29]。

计算榆林市12个气象站点1981—2018年逐年作物生育期LDFAI，采用线性倾向法、M-K突变检验分析时间序列变化，使用ARCGIS软件运用反距离权重插值法(IDW)对旱涝急转事件进行空间演变分析，基于非旱涝急转事件，使用SPSS对不同旱涝急转类型的作物产量进行独立性T检验。分析过程中将各年代依次表示为1980年代(1981—1990)、1990年代(1991—2000)、2000年代(2001—2010)、2010年代(2011—2018)。

2 结果与分析

2.1 榆林市作物生育期旱涝急转事件时间变化及趋势分析

1981—2018年，榆林市LDFAI在时间上以0.148/10年的速率变化(图2a)，ILDFAI以0.076/10年的速率变化(图2b)。LDFAI在年代上变化幅度逐渐减小，在1上下波动；ILDFAI在年代上具有一定的阶段性，1980年代平均强度相对较弱，1990年代平均强度相对最高，从1990年代往后逐年代降低。图2c和2d可以看出，LDFAI和ILDFAI的正序列曲线UF在1981—1994年间的年份均为负值，这表明LDFAI和ILDFAI在1981—1994年间波动下降，其中1983—1987年超过-1.96(α=0.05的显著性水平信度线)，表明1983—1987年后榆林市“旱转涝”事件发生频率显著减少，“涝转旱”事件显著增多，旱涝急转发生强度显著降低；1994年以后在显著性水平信度线的正值区间内平稳波动，表明1994年以后“旱转涝”事件发生频率趋于增多，旱涝急转强度升高。UFk和UBk序列在显著水平信度线内交叉点的数量可确定突变[41]，榆林市38年间LDFAI和ILDFAI的UF、UB正反序列曲线在显著性水平信度线(α=0.05)之间有多个交点，表明旱涝急转事件及强度未发生突变。

图2 1981—2018年榆林市作物生育期内长周期旱涝急转指数(LDFAI)及强度(ILDFAI)的变化及M-K检验结果

表1 1981—2018年榆林市作物生育期内旱涝急转事件特征

2.2 榆林市典型年份作物生育期旱涝及LDFAI特征

根据表1中的LDFAI结果，结合《中国气象大典·陕西卷》和榆林市各县县志中查到的在5—8月有旱灾、又发生涝灾的年份，选出旱转涝、涝转旱发生典型年份如下：①对于旱转涝年，1994—1995年最为典型。1994年前半年大旱，到7月26日和8月5日出现两次历史罕见大暴雨，定边县经济损失达1亿元以上[37]；1994年秋至1995年7月，榆林市连旱300 d，1995年8月底至9月上旬，大部分地区出现了十多天的连阴雨天气，数万公顷马铃薯就地腐烂[39]。②对于涝转旱年，其中1991年最为严重，5月下旬榆林市部分地区突发暴雨，7月底榆林市大部分沟溪断流，万余公顷耕地绝收[36]。综上，旱涝急转年的5—6月Pa和7—8月Pa之差的绝对值在50%以上，典型旱转涝年份5—8月降水量偏多，但5—6月偏少30%以上，7—8月偏多30%以上，LDFAI均大于1.8；典型涝转旱年份除1983年外其余年份5—8月降水量偏少但接近正常，5—6月偏多30%以上，旱情严重的年份7—8月偏少40%以上，LDFAI绝对值小于0.7。可见，LDFAI和降水距平的组合可以较好地识别出榆林市长周期旱涝急转事件的典型年份。

2.3 作物生育期旱涝急转事件空间分布特征

从图3a和图3b可得，榆林市1981—2018年多年平均旱涝急转指数及强度在空间上分布一致，均呈现为“东高西低”的特征。多年平均高值区分布在以吴堡、神木为代表的东部地区，低值区分布在以定边、靖边、横山为代表的西部山区。图3c和图3d可得，榆林市旱转涝事件与涝转旱事件发生频率相当，且分布规律相对一致，整体呈现“南高北低”的特点。其中，旱转涝事件发生频率最高的是子洲，发生频次为10次，府谷和佳县发生频次相对较少，为5次；涝转旱事件子洲发生频次较高，为9次，府谷发生频次最少，为5次。总体来看，子洲旱涝急转事件总发生频次为19次，发生频率全区最高，超过平均值的35.71%，府谷发生频次为10次，低于全区发生平均次数的28.57%，发生频率全区最低。

2.4 榆林市作物生育期内旱涝急转事件对作物产量的影响

在表2中，按旱涝急转事件划分标准，将1981—2018年榆林市12个区县的旱涝急转事件进行统计发现，12个区县累计发生“旱转涝”事件74次，“涝转旱”事件76次。将作物总产量和单产按照旱转涝(D-F)、涝转旱(F-D)、非旱涝急转分为三类，基于非旱涝急转比较不同类型旱涝急转事件下作物产量的变化。遭遇旱转涝事件，作物产量变化不明显(P>0.05)；遭遇涝转旱事件，粮食单产、大豆单产、薯类单产分别显著减少17.85%(P<0.05)、16.63%(P<0.05)、19.46%(P<0.01)。

图3 1981—2018年榆林市作物生育期平均LDFAI、ILDFAI及旱涝急转事件发生频次的空间分布示意图

表2 榆林市作物生育期内旱涝急转事件对作物产量的影响

3 讨论

通过榆林市5—8月降水变化发现，作物生育期旱涝急转发生频率为34.37%。其成因机制异常复杂，南部丘陵沟壑区复杂的地形必然会加强对流天气的发生和发展程度[40]。榆林市5—8月的水汽主要来源于印度洋和太平洋，深受西太平洋副热带高压和南边界的水汽输送影响[40,42]。在涝转旱年涝期(5—6月)，榆林市位于副高西侧和蒙古低压底部，偏弱的西风带使中高纬度冷空气南下，西南暖湿气流充足北上，降水激增；旱期(7—8月)，副高虽然偏强但在西太平洋上，西风带偏强阻止了中高纬度冷空气南下，使榆林市降水偏少。旱转涝年机制则相反。

由表2可知，粮食单产、大豆单产、薯类单产分别显著减少17.85%(P<0.05)、16.63%(P<0.05)、19.46%(P<0.01)。遭遇涝转旱事件，前期降水过多导致植株芽期的分化与根系养分的积累都无法达到正常水平，生育期长的大豆播种时间越晚其产量损失越严重[16]，宜选育短生育期品种大豆保证作物产量。马铃薯结薯期对温度较敏感，且需水量占全生育期耗水量的50%以上，干热的天气会严重抑制块茎生长[17]，应对马铃薯进行阴坡种植，控温保水，并选择早熟品种以避免气象灾害威胁。大豆、马铃薯和玉米前期幼苗期以地下根系生长为主，需水量较少[17,43]，前期轻微的干旱在一定程度上能缓解、补偿后期湿涝对大豆的胁迫[44]，后期结实期需水量较高，因此大豆、马铃薯和玉米受旱转涝事件产量变化不显著(P>0.05)，为避免前期旱情对玉米、大豆的影响，宜选择抗旱品种进行覆膜湿播，生育期后期及时补充灌溉。为保障农民收入，需进一步丰富政策性农业气象灾害保险体系。现有政策性农业气象灾害保险仅针对单一作物设保，间种或套种的其他作物暂未纳入其中[45]，因此可考虑将玉米、马铃薯、大豆的混种比例纳入气象灾害保险产品设计参考范围。

4 结论

经过对1981—2018年间榆林市作物生育期旱涝急转事件进行分析，发现榆林市旱涝急转事件的发生频率为34.37%。旱涝急转强度具有阶段性特征，1980年代平均强度相对较弱，1990年代平均强度相对最高，从1990年代往后逐年代降低。榆林市多年平均旱涝急转指数及强度在空间上分布一致，均呈现为“东高西低”的特征，多年平均高值区分布在东部地区，多年平均低值区分布在西部山区。涝转旱事件对作物产量影响显著，发生涝转旱事件，粮食单产、大豆单产、薯类单产分别显著减少17.85%(P<0.05)、16.63%(P<0.05)、19.46%(P<0.01)。根据榆林市农业生产特点提出旱涝急转农业综合应对措施建议，包括加强作物灌溉管理、作物品种优化和丰富农业气象灾害保险体系，以期为北方地区生态脆弱区应对旱涝急转事件提供参考。