西藏青稞主要种植区干旱时空分布及致灾危险性评估.pdf

第46卷第7期2023年7月ARIDLANDGEOGRAPHYVol.46No.7Jul.2023doi：10.12118/j.issn.10006060.2022.535西藏青稞主要种植区干旱时空分布及致灾危险性评估史继清1，甘臣龙2，周刊社1，袁雷1，张东东1（1.西藏自治区气候中心，西藏 拉萨850000；2.墨竹工卡县气象局，西藏 拉萨850000）摘要：西藏青稞主要种植区是西藏地区主要粮食作物的产区，在气候变化和人类活动加剧的背景下，该区的农业深受干旱影响。根据 气象干旱等级（GB/T 20481-2017）国家标准，计算西藏青稞主要种植区25个气象站点全生育期的逐日气象干旱综合指数（MCI），统计19812020年各站点的气象干旱过程，并进一步分析了西藏青稞主要种植区各生育期气象干旱的时空变化特征及致灾危险性等级分布。结果表明：（1）近40 a来西藏青稞主要种植区呈现气温显著上升（P-0.5轻旱-1.0MCI-0.5中旱-1.5MCI-1.0重旱-2.0MCI-1.5特旱MCI-2.011007期史继清等：西藏青稞主要种植区干旱时空分布及致灾危险性评估态，该方法既可以取空间不同站点进行分解，也可以对同一站点的不同时间、不同高度的多种要素进行综和分析，现广泛应用于气象学领域28；同时，EOF分解的展开、收敛速度均很快，变量场的信息能够很容易地集中在前几个模态上；累计方差贡献率达到一定程度的前几个模态具有一定物理意义，可有效反映干旱强度时空场的变化特征29，基本思路是将包含g个空间点h个年份的观测值构成的要素场（Xgh）分解为相互正交的空间函数矩阵（Vgp）和时间函数矩阵（Tph）的乘积，即:Xgh=p=1g(VgpTph),p=1,2,g（6）其中，Vgp=v11v12v1pvg1vg2vgp（7）Tph=t11t12t1pth1th2thp（8）2结果与分析2.1 西藏青稞主要种植区气温、降水变化趋势由图2可知，19812020年西藏青稞主要种植区多年平均降水量为490.92 mm，且以7.71 mm(10a)-1的速率呈不显著增加趋势，其中，1983年和1998年分别是近40 a该研究区年平均降水量最少（350.11 mm）和最多（611.20 mm）的年份。该研究区多年平均气温为6.52，并以0.34(10a)-1的速率呈显著的上升趋势（P0.01），其中1997年和2009年分别是过去40 a 该研究区年平均气温最低（5.46）和最高（7.58）的年份。由图3可知，近40 a西藏青稞主要种植区春季多年平均降水量为 84.99 mm，占年平均降水量的17.31%，并以4.83 mm(10a)-1的速率呈显著增加趋势（P0.05）；夏季多年平均降水量为296.86 mm，占年平均降水量的60.47%，并以6.71 mm(10a)-1的速率呈不显著增加趋势；秋、冬季多年平均降水量分别为94.69 mm和14.39 mm，分别以-3.22 mm(10a)-1和-0.42 mm(10a)-1的速率呈不显著的减少趋势。春、夏、秋、冬季多年平均气温分别为 6.58、14.04、7.05 和-1.63，并分别以0.27(10a)-1、0.26(10a)-1、0.38(10a)-1、0.46(10a)-1的速率呈现显著的上升趋势（P0.01）。2.2 干旱覆盖范围的年际变化特征19812020年西藏青稞主要种植区各生育期干旱覆盖范围的变化趋势如图4所示。播种-分蘖期和抽穗-成熟期，干旱覆盖范围均呈不明显的下降趋势，多年均值分别为 42.7%和31.4%。播种-分蘖期，各年份均有一定范围的干旱发生；从各年代来看，19912000年平均干旱覆盖范围最广，为51.6%，20012010年最小，为32.4%。抽穗-成熟期，1983年干旱覆盖范围最广，为84.0%，2000、2004 年以及 20162017 年干旱覆盖范围为0，表明这些年份研究区无干旱发生；各年代中，19811990年干旱发生范围最大，为45.3%，20012010年最小，为24.0%。图2 19812020年西藏青稞主要种植区年平均降水量与年平均气温的变化Fig.2 Changes of annual mean precipitation and annual mean temperaturein highland barley growing areas of Tibet from 1981 to 2020110146卷图4 19812020年西藏青稞主要种植区各生育期干旱覆盖范围的变化Fig.4 Change of drought coverage in different growth stages in highland barley growing areas of Tibet from 1981 to 2020图3 19812020年西藏青稞主要种植区季节平均降水量与平均气温的变化Fig.3 Changes of seasonal mean precipitation and mean temperaturein highland barley growing areas of Tibet from 1981 to 202011027期史继清等：西藏青稞主要种植区干旱时空分布及致灾危险性评估分蘖-抽穗期和全生育期，干旱覆盖范围均呈现显著（P0.05）的下降趋势，多年均值分别为51.3%和36.8%。分蘖-抽穗期共有6 a覆盖范围超过 80.0%，其中，1982 年和 1986 年覆盖范围均为88.0%，发生范围最广，2000年无干旱发生；各年代中，19811990 年干旱最普遍，平均覆盖范围为72.0%，20012010 年 平 均 覆 盖 范 围 最 少，为32.4%。全生育期，1992年干旱覆盖范围最广，为92.0%，2004年和2008年无干旱发生，其余年份干旱覆盖范围不一；各年代中，19811990年干旱覆盖范围最广，为 54.2%，20012010 年最小，仅为19.6%。2.3 干旱持续日数的年际变化特征由图5可知，播种-分蘖期、分蘖-抽穗期、抽穗-成熟期和全生育期的干旱持续日数均呈现不显著的下降趋势，多年均值分别为24 d、26 d、25 d和60 d。播种-分蘖期，2000年干旱持续日数最长，为35 d，2004年最短，为14 d；各年代干旱持续日数介于2325 d之间。分蘖-抽穗期，干旱持续日数最长的是1988年的39 d，最短的是2000年的0 d，无干旱发生；各年代中，19811990年干旱平均持续日数最长，约30 d，20012010年最短，约20 d，其余年代均为28 d。抽穗-成熟期，1983年干旱持续日数最长，为49 d，2000、2004年和20162017年均无干旱发生；各年代中，19811990年干旱平均持续日数最长，约31 d，20012010年最短，约21 d。全生育期，干旱持续日数最长的是1984年的136 d，最短的是 2000 年和 2008 年的 0 d；各年代中，19811990年干旱平均持续日数最长，约73 d，20012010年最短，约44 d，其余年代在55 d左右。2.4 干旱频率的空间分布特征由图6可以看出，整体上，分蘖-抽穗期干旱频率最高，播种-分蘖期次之，抽穗-成熟期最低。播种-分蘖期，林芝市、错那站的干旱发生率普遍偏低，均在30.0%以下；研究区东部、西部干旱发生率多在40.0%60.0%之间，占总站点数的68.0%。分蘖-抽穗期，介于50.0%70.0%之间的站点数最多，占比达64.0%，主要分布在研究区中西部和昌都市中南部；干旱发生率在30.0%以下和70.0%以上的站点均为1个，分别出现在波密站和江孜站。抽穗-成熟期，干旱发生率均在50.0%以下（除隆子站外），低值区主要位于昌都市大部、林芝市中部地区。全生图5 19812020年西藏青稞主要种植区各生育期的干旱持续日数的变化趋势Fig.5 Change trends of drought duration days in different growth stages of highland barleygrowing areas in Tibet from 1981 to 2020110346卷育期，干旱发生率高值区集中在日喀则市中东部、昌都市西部以及山南市北部。2.5 干旱强度的时空特征利用EOF对西藏青稞主要种植区各生育期干旱强度的标准值进行分解（表3、图7），分析结果表明，播种-分蘖期、分蘖-抽穗期、抽穗-成熟期和全生育期前5个特征向量特征值的累计贡献率分别达到71.2%、72.1%、83.5%和79.3%，其中播种-分蘖期和全生育期仅有第一个特征根的误差范围不重叠，分蘖-抽穗期和抽穗-成熟期则有前2个特征根的误差范围不重叠。通过North显著性检验，播种-分蘖期（第一载荷量）、分蘖-抽穗期（第一、二载荷量）、抽穗-成熟期（第一、二载荷量）和全生育期（第一载荷量）的累计贡献率分别接近23.1%、51.8%、68.7%和35.8%，可以很好地解释近40 a来西藏青稞主要种植区各生育期干旱强度的主要时空分布特征。播种-分蘖期（图7a），第一模态（EOF1）的干旱强度负值区集中在雅鲁藏布江一线大部、林芝站、米林站，其余为正值区。从EOF1时间序列（图8a）可以看出，19891998年和20072020年干旱强度时间系数在正负间交替变化，说明这2个时间段研究区干旱强度的强弱变化明显；19982006年均为正值，表明该时间段研究区干旱强度减弱，结合图2发现，19982005年降水量偏多，因此对干旱强度产生一定的影响。分蘖-抽穗期、抽穗-成熟期和全生育期（图7b、d、f）的EOF1在整个研究区大体上为正值（除抽穗-成熟期的聂拉木站和全生育期的丁青站以外），表明近40 a来研究区在这3个生育期的干旱强度变化趋势空间上具有高度的一致性，即干旱强度在研究区普遍增强或减弱。高值中心主要位于雅鲁藏布江一线的山南市北部、拉萨市南部和日喀则市东部，是干旱强度变化反应最敏感的3个区域。从图8bd可以得出，分蘖-抽穗期、抽穗-成熟期和全生育期的时间系数分别在19992004年、19952004年和19962004年为正值，表明这3个时间段研究区干旱强度减弱，与相同时段的降水量偏多有一定的关系；其余时间段时间系数正负交替变化。分蘖-抽穗期的第二模态（EOF2）（图7c）的干旱强度整体表现为研究区中部正值中心与东部、西部边缘负值中心呈相反的空间分布特征；山南市、日喀则市东部、拉萨市南部、林芝市局部边缘地区图6 19812020年西藏青稞主要种植区各生育期干旱频率的空间分布Fig.6 Spatial distributions of drought frequency in different growth stages of highland barleygrowing areas in Tibet from 1981 to 202011047期史继清等：西藏青稞主要种植区干旱时空分布及致灾危险性评估以及丁青站是正值区，其余为负值区。结合图8b发现，20002006年和20172020年主要表现为研究区中部干旱强度减弱，而东部、西部边缘增强；20072012年呈现相反的变化趋势。抽穗-成熟期的第二模态（EOF2）（图7e）的干旱强度表现为“西南负、东北正”的反向分布型，即西南东北差异型；负值区主要集中在雅鲁藏布江一线南段、沿喜马拉雅山一带以及左贡站、芒康站、察隅站，其他属于正值区。由图8c发现，19812020年的时间系数在正负间振荡，表明研究区干旱强度年际间的变化趋势显著。2.6 干旱致灾危险性等级的空间分布利用ArcGIS 10.4软件，分别将播种-分蘖期、分蘖-抽穗期、抽穗-成熟期和全生育期的干旱致灾危险性划分为低、较低、中等、较高、高等5个等级（图9）。干旱致灾危险性在播种-分蘖期，大体上沿洛隆八宿左贡一线，西部地区呈现由东向西逐渐加重的趋势；分蘖-抽穗期旱灾危险性均偏高；抽穗-成熟期，整体上呈现研究区西高东低的分布特征；全生育期旱灾危险性处于中等或偏低。播种-分蘖期，干旱致灾危险性中高等级集中分布在日喀则市中部，其他地区危险性等级偏低，其中林芝市大部和山南市南部最低。分蘖-抽穗期，雅鲁藏布江一线大部、昌都市西南部以及隆子站的干旱致灾危险性最高；林芝市中部和聂拉木站、错那站处于中等危险性，其他地区属于较高等级。抽穗-成熟期，除雅鲁藏布江一线大部和错那站、隆子站、八宿站外，其他地区旱灾的危险性多属于中等级别以下，其中聂拉木站和芒康站危险性最低。全生育期，干旱致灾危险性等级均处于中等以下，林芝站和米林站最低；中等危险性区域主要位于雅鲁藏布江一线，其他区域危险性较低。3讨 论（1）结合前面的干旱过程覆盖范围，分析各种空间模态干旱强度随时间的变化特征。播种-分蘖期，EOF1在1990、1992、1997年和2015年这4 a均为干旱强年。分蘖-抽穗期，EOF1 在 1982、19861987、2000年和20082009年其时间系数绝对值较高，结合图4b，发现1982、19861987年和2009年为明显的干旱强年，其中2009年的时间系数绝对值最高，当年对应的干旱覆盖范围也大（84%），这与袁雷等22对西藏西部地区雨季得出的“干旱强度在19811990年和20012010年偏高，在2009年达到峰值”的结轮较一致；EOF2 在 1993、1995 年和2018年为干旱强年。抽穗-成熟期，结合图4c，发现EOF1在1983、1986、1992年和2009年这4 a为突出的干旱强年，其中1983年的时间系数绝对值最高，同时干旱影响范围为最大值（84%）；EOF2在19821983、1986年和1992年为典型的干旱强年；近40 a来，抽穗-成熟期2个模态的干旱强年与袁雷等22对西藏沿江一线雨季得出的“干旱强度在19811990年达到最大”的结轮较一致。全生育期，结合图4d，发现 EOF1 在 19821983、19861987、1992 年和2009年为显著的干旱强年。（2）文中增加了20142020年时间序列，得出表3 EOF前5个载荷量方差贡献率及累计方差贡献率Tab.3 Variance contribution rate and cumulative variance contribution rate of the first five loads of EOF生育期播种-分蘖期抽穗-成熟期EOF第一载荷量第二载荷量第三载荷量第四载荷量第五载荷量第一载荷量第二载荷量第三载荷量第四载荷量第五载荷量方差贡献率/%23.118.813.68.66.861.57.25.44.84.6累计方差贡献率/%23.141.955.864.471.261.568.774.178.983.5生育期分蘖-抽穗期全生育期EOF第一载荷量第二载荷量第三载荷量第四载荷量第五载荷量第一载荷量第二载荷量第三载荷量第四载荷量第五载荷量方差贡献率/%37.913.98.16.65.635.824.47.37.14.7累计方差贡献率/%37.951.859.966.572.135.860.267.474.679.3注：EOF为经验正交函数法。110546卷19811990年干旱覆盖范围最广、持续日数最长、发生频率最大，20012010年均最小（短），与熊俊楠9等得出“旱灾在19912005年发生较为频繁，2002年后每年发生频次均低于多年平均值”的结论不尽相同，分析发现时间序列的增加并未改变干旱频率的低发时间段，出现高发频率干旱年份不一致的原因可能是本文研究仅涉及全年的部分月份（49月），虽然13月和1012月的干旱灾害发生相对较少，但也会发生。不同生育期干旱发生频率的高值区主要集中在日喀则市东部、拉萨市南部、昌都市西南部和山南市中部地区，低值区位于林芝市中部，这与熊俊楠等9得出的区域范围较为吻合，但与袁雷等22对西藏雨季干旱频率的结论有所差别，原因可能是时间序列、研究区域、干旱指标的选取不同造成。（3）综合4个生育期干旱致灾危险性等级，西藏青稞主要种植区干旱致灾危险性较高的区域主要分布在日喀则市中东部、拉萨市南部和山南市北部，该区域多位于雅鲁藏布江一线，降水相对偏少，蒸散较大，更易发生干旱；较低区域主要集中在林芝市，该区域属于热带湿润和半湿润气候，降水丰富、蒸发量相对较小，干旱事件相对少发。熊俊楠注：EOF1、EOF2分别为经验正交函数法（EOF）的第一模态、第二模态。下同。图7 19812020年西藏青稞主要种植区各生育期干旱强度EOF不同模态空间分布特征Fig.7 Spatial distribution characteristics of different modes of drought intensity EOF in different growth stagesof highland barley growing areas in Tibet from 1981 to 202011067期史继清等：西藏青稞主要种植区干旱时空分布及致灾危险性评估图9 19812020年西藏青稞主要种植区各生育期干旱致灾危险性等级的空间分布Fig.9 Spatial distributions of drought hazard grades in different growth stages of highlandbarley growing areas in Tibet from 1981 to 2020图8 19812020年西藏青稞主要种植区各生育期干旱强度EOF不同模态对应的时间系数Fig.8 Time coefficient of EOF of different modes in different growth stages of highlandbarley growing areas in Tibet from 1981 to 2020110746卷等9基于历史气象灾害事件、史继清等11利用自然水分亏缺率指数对西藏地区气象干旱空间分布特征所得出的结论与本文结论较一致；而金建新等21的研究结果认为“干旱指数从西向东逐渐递增，从西部地区的特干旱气候区依次向东部的湿润半湿润及湿润气候区过渡”，与本文结论有所差异，究其原因在于选取的指标不一致、两者的计算方法存在差异等。研究区分蘖-抽穗期致灾危险性总体较高，可能由于分蘖-抽穗期是青稞需水关键期，且此时期降水量较小、蒸发量大、气温回升快，最终导致了该生育期干旱最为严重，进而证实了熊俊楠等9得出“旱灾多发生在上半年，从1月开始逐渐递增，在6月达到峰值”的研究结论。4结 论（1）19812020年西藏青稞主要种植区的气温呈显著上升趋势，降水呈不显著增加趋势，季节性气温和降水变化趋势差异明显。多年来秋冬季降水呈减少趋势，春夏季则呈现不同速率的增加趋势；四季的年际气温变化均呈显著上升趋势，其中冬季气温上升速率最快，达到0.46(10a)-1。（2）19812020年西藏青稞主要种植区不同生育期干旱覆盖范围和持续日数总体呈现不同程度的下降趋势，干旱强度均有所增强（除播种-分蘖期外）；分蘖-抽穗期干旱频率最高，播种-分蘖期次之，抽穗-成熟期最低。（3）通过EOF分析的西藏青稞主要种植区干旱强度分布特点：播种-分蘖期呈现中部与中西部相反；分蘖-抽穗期呈现全区一致、中部与中西部相反；抽穗-成熟期呈现全区一致（除聂拉木站外）、西南与东北相反；全生育期呈现全区一致的多种模态。中部与中西部相反的典型年份为2015年（播种-分蘖期，定日站、拉孜站严重干旱）和1995年（分蘖-抽穗期，昌都市大部、林芝市东部严重干旱），西南与东北相反典型年份为1986年（日喀则市东部严重干旱）。（4）西藏青稞主要种植区干旱致灾危险性等级在播种-分蘖期和分蘖-抽穗期时，整体上呈现中部低两边高的格局特征；处于抽穗-成熟期和全生育期，整体上呈现西高东低的的分布特征。参考文献（References）1姜大膀,王晓欣.对IPCC第六次评估报告中有关干旱变化的解读J.大气科学学报,2021,44(5):650-653.Jiang Dabang,Wang Xiaoxin.A brief interpretation of drought change from IPCCSixth Assessment ReportJ.Transactions of Atmospheric Sciences,2021,44(5):650-653.2王晓利,张春艳,侯西勇.19612017年环渤海地区气象干旱时空特征及致灾危险性评估J.生态学报,2019,39(13):4647-4659.Wang Xiaoli,Zhang Chunyan,Hou Xiyong.Spatial-temporal characteristics and hazard risks of meteorological drought inCircum-Bohai-Sea region from 1961 to 2017J.Acta Ecologica Sinica,2019,39(13):4647-4659.3Freitas A A,Drumond A,Carvalho V S B,et al.Drought assessment in So Francisco River Basin,Brazil:Characterizationthrough SPI and associated anomalous climate patternsJ.Atmosphere,2022,13(1):41,doi:10.3390/atmos13010041.4Yang B,Kong L,Lai C G,et al.A framework on analyzing long-term drought changes and its influential factors based on the PDSIJ.Atmosphere,2022,13(7):1151,doi:103390/atmos13071151.5Suroso,Nadhilah D,Ardiansyah,et al.Drought detection in JavaIsland based on standardized precipitation and evapotranspirationindex(SPEI)J.Journal of Water and Climate Change,2021,12(6):2734-2752.6费龙,邓国荣,张洪岩,等.基于降水Z指数的朝鲜降水及旱涝时空特征J.自然资源学报,2020,35(12):3051-3063.FeiLong,Deng Guorong,Zhang Hongyan,et al.Spatiotemporal patterns of precipitation and drought and flood using Z-index in Democratic People s Republic of KoreaJ.Journal of Natural Resources,2020,35(12):3051-3063.7周丽,谢舒蕾,吴彬.基于CI和强度分析方法的四川冬春季干旱事件变化特征J.自然灾害学报,2020,29(3):36-44.ZhouLi,Xie Shulei,Wu Bin.Variation characteristics of the winter andspring drought events in Sichuan based on CI and intensity analysisJ.Journal of Natural Disasters,2020,29(3):36-44.8曲学斌,杨钦宇,王慧清,等.基于MCI的内蒙古气象干旱强度特征分析J.气象与环境科学,2019,42(4):47-54.Qu Xuebin,Yang Qinyu,Wang Huiqing,et al.Characteristics of meteorological drought intensity in Inner Mongolia based on MCIJ.Meteorological and Environmental Sciences,2019,42(4):47-54.9熊俊楠,刘志奇,范春捆,等.19832013年西藏自治区气象灾害时空分布特征与变化趋势J.冰川冻土,2017,39(6):1221-1231.Xiong Junnan,Liu Zhiqi,Fan Chunkun,et al.Temporaland spatial distribution characteristics and changing trend of meteorologic disaster in Tibet Autonomous Region from 1983 through2013J.Journal of Glaciology and Geocryology,2017,39(6):1221-1231.11087期史继清等：西藏青稞主要种植区干旱时空分布及致灾危险性评估10 史继清,豆永丽,杨霏云,等.西藏地区潜在蒸散量时空格局特征及影响因素研究J.干旱区研究,2021,38(3):724-732.ShiJiqing,Dou Yongli,Yang Feiyun,et al.Temporal and spatial pattern characteristics of potential evapotranspiration in Tibet and itsinfluencing factorsJ.Arid Zone Research,2021,38(3):724-732.11 史继清,杨霏云,边多,等.基于干旱灾害风险综合评估指数的西藏主要农区青稞干旱时空格局J.中国农学通报,2021,37(2):80-87.Shi Jiqing,Yang Feiyun,Bian Duo,et al.Temporaland spatial patterns of barley in Tibet major agricultural areasbased on the comprehensive assessment index of drought disasterriskJ.Chinese Agricultural Science Bulletin,2021,37(2):80-87.12 高佳佳,徐薇,边央,等.基于FAO Penman-Monteith方法的西藏春青稞干旱趋势研究J.中国农学通报,2021,37(11):79-86.Gao Jiajia,Xu Wei,Bian Yang,et al.The drought tendency ofspring highland barley on Tibet Plateau:A research via the FAOPenman-Monteith methodJ.Chinese Agricultural Science Bulletin,2021,37(11):79-86.13 宋善允,王鹏祥.西藏气候M.北京:气象出版社,2013:11-13,385-387.Song Shanyun,Wang Pengxiang.Tibet climateM.Beijing:China Meteorology Press,2013:11-13,385-387.14 蒙强,刘静霞,李玉庆,等.土壤水分下限调控灌溉对春青稞生长、产量和水分利用的影响J.干旱地区农业研究,2020,38(1):173-182.Meng Qiang,Liu Jingxia,Li Yuqing,et al.Effectsof regulated irrigation of lower limits of soil moisture on growth,yield and water use efficiency of spring highland barleyJ.Agricultural Research in the Arid Areas,2020,38(1):173-182.15 时学双,李法虎,闫宝莹,等.不同生育期水分亏缺对春青稞水分利用和产量的影响J.农业机械学报,2015,46(10):144-151,265.Shi Xueshuang,Li Fahu,Yan Baoying,et al.Effects ofwater deficit at different growth stages on water use and yield ofspring highland barleyJ.Journal of Agricultural Machinery,2015,46(10):144-151,265.16 郑盐源.日土县气象灾害对青稞种植的影响及其防御对策J.现代农业科技,2017(23):172-176.Zheng Yanyuan.Effects ofmeteorological disasters on highland barley cultivation in RutogCounty and its countermeasuresJ.Modern Agricultural Scienceand Technology,2017(23):172-176.17 邓国卫,孙俊,赖江,等.四川水稻不同生育期干旱与产量灰色关联分析J.干旱气象,2022,40(5):814-822.Deng Guowei,Sun Jun,Lai Jiang,et al.Grey correlation analysis of drought andyield at different growth stages of rice in Sichuan ProvinceJ.Journal of Arid Meteorology,2022,40(5):814-822.18 黄岩,李晶,王莹,等.不同生育期干旱对玉米生长及产量的影响模拟J.农业灾害研究,2019,9(6):47-49,92.Huang Yan,LiJing,Wang Ying,et al.Impact simulation of drought on maizegrowth and yield in different growth stagesJ.Study on Agricultural Disaster,2019,9(6):47-49,92.19 曹永强,王怡涵,冯兴兴,等.河北省夏玉米不同生育期干旱时空分析J.华北水利水电大学学报(自然科学版),2020,41(4):1-9.Cao Yongqiang,Wang Yihan,Feng Xingxing,et al.Spatio-temporal analysis of drought in different growth stages of summermaize in Hebei ProvinceJ.Journal of North China University ofWater Resources and Electric Power(Natural Science Edition),2020,41(4):1-9.20 吴乾慧.19602015年北方冬麦区冬小麦生育期农业干旱演变特征及风险区划D.兰州:西北师范大学,2018.Wu Qianhui.Agricultural drought evolution and risk regionalization of winter wheat in north China during growth period from 1960 to 2015D.Lanzhou:Northwest Normal University,2018.21 金建新,张娜,桂林国.西藏地区干旱指标的时空演变J.水土保持研究,2019,26(5):377-380.Jin Jianxin,Zhang Na,Gui Linguo.Temporal and spatial variations of the drought index to the Tibetan PlateauJ.Research of Soil and Water Conservation,2019,26(5):377-380.22 袁雷,刘依兰,马鹏飞.基于标准化降水指数的19812013年西藏干旱时空特征分析J.中国农学通报,2015,31(25):228-234.Yuan Lei,Liu Yilan,Ma Pengfei.Temporal and spatial patterns of drought based on standard precipitation index(SPI)in Tibet during 19812013J.Chinese Agricultural Science Bulletin,2015,31(25):228-234.23 孙晓光,王腾,李白萍,等.基于标准化降水指数的西藏昌都市旱涝时空特征J.中国农学通报,2017,33(8):112-118.SunXiaoguang,Wang Teng,Li Baiping,et al.Spatial and temporalcharacteristics of drought-flood in Chamdo based on standard precipitation index(SPI)J.Chinese Agricultural Science Bulletin,2017,33(8):112-118.24 王圆圆,扎西央宗.利用条件植被指数评价西藏植被对气象干旱的响应J.应用气象学报,2016,27(4):435-444.Wang Yuanyuan,Zhaxi Yangzong.Assessing vegetation response to meteorological drought in Tibet Autonomous Region using vegetation indexJ.Journal of Applied Meteorological Science,2016,27(4):435-444.25 杜军,袁雷,周刊社,等.西藏自治区县级青稞种植气候适宜性区划M.北京:气象出版社,2017.Du Jun,Yuan Lei,Zhou Kanshe,et al.Climate suitability zoning for highland barley cultivation at County TibetM.Beijing:China Meteorology Press,2017.26张强,邹旭恺,肖风劲,等.气象干旱等级M.北京:气象出版社,2008:33-53.Zhang Qiang,Zou Xukai,Xiao Fengjin,et al.Meteorological drought gradeM.Beijing:China Meteorology Press,2008:33-53.27Jung J.Spatial analysis of the vulnerability to meteorological hazards in KoreaJ.Journal of Climate Research,2018,13(3):211-229.28 史继清,边多,杨霏云,等.西藏地区潜在蒸散量变化特征及灰色模型预测初探J.干旱区地理,2021,44(6):1570-1579.Shi110946卷Jiqing,Bian Duo,Yang Feiyun,et al.Variation characteristics ofpotential