基于CWDI蒙冀半干旱区近60a谷子干旱时空变化特征 (1).pdf

第 4 期 赵凌暄等：基于 CWDI 蒙冀半干旱区近 60a 谷子干旱时空变化特征 419 中国农业气象（Chinese Journal of Agrometeorology）2024 年doi:10.3969/j.issn.1000-6362.2024.04.008 赵凌暄,王靖,李扬,等.基于 CWDI 蒙冀半干旱区近 60a 谷子干旱时空变化特征J.中国农业气象,2024,45(4):419-430 基于 CWDI 蒙冀半干旱区近 60a 谷子干旱时空变化特征 赵凌暄，王 靖*，李 扬，王孝贤，赵熙玲，陈仁伟，胡 琦，张佳莹，王慧烨，赵庚云（中国农业大学资源与环境学院，北京 100193）摘要：谷子是半干旱区的特色作物，具有耐旱稳产的特点。蒙冀半干旱区是谷子的优势产区，种植面积及总产量均居全国前列，但干旱严重威胁着该区谷子生产。为揭示干旱对谷子生长过程的影响，提高谷子生产应对干旱风险的能力，选取研究区内 27 个气象站点 19612019 年谷子生长季（59 月）逐日气象数据，采用作物水分亏缺指数（CWDI）作为干旱指标，结合干旱强度、干旱频率、干旱站次比分析谷子干旱的时空变化特征。结果表明：（1）近 60a，蒙冀半干旱区谷子干旱强度呈下降趋势。干旱的发生频率和站次比均随干旱等级的升高而降低。（2）在谷子各生育阶段，干旱强度和干旱频率表现为拔节期苗期成熟期灌浆期抽穗期。（3）空间分布上，研究区内谷子干旱强度和干旱频率均表现为由西北向东南逐渐降低，其中中旱、重旱和特旱的发生范围均减小。干旱强度和频率的降低有利于蒙冀半干旱区的谷子生产。关键词：半干旱区；CWDI 干旱指数；谷子；生育阶段；干旱强度与频率 Spatial and Temporal Variation Characteristics of Millet Drought in Semi-arid Region of Inner Mongolia and Hebei Based on CWDI in Recent 60 Years ZHAO Ling-xuan,WANG Jing,LI Yang,WANG Xiao-xian,ZHAO Xi-ling,CHEN Ren-wei,HU Qi,ZHANG Jia-ying,WANG Hui-ye,ZHAO Geng-yun(College of Resources and Environment,China Agricultural University,Beijing 100193,China)Abstract:Millet is a characteristic crop of semi-arid regions,with high drought tolerance and stable yield.The semi-arid region of Inner Mongolia and Hebei is the dominant millet production region,with the highest planting area and total yield of millet in China.Drought is threatening millet production in this region.In order to reveal the impact of drought on the millet growth processes and increase the capacity of millet production to cope with drought,the crop water deficit index(CWDI)was used to investigate millet drought in this study.Based on the daily meteorological data from 27 meteorological stations in the study area from 1961 to 2019,the temporal and spatial variation characteristics of the millet drought were analyzed by combining drought intensity,drought frequency,and drought station ratio.The results showed that:（1）in the past 60 years,the drought intensity of millet in the semi-arid area of Inner Mongolia and Hebei showed a significant downward trend.The frequency and station ratio of millet drought decreased with increased drought level.（2）The drought intensity and frequency in each growth stage of millet was in order of:jointing periodseedling periodmaturity periodgrain filling periodheading period.（3）The drought intensity and frequency of millet decreased gradually from northwest to southeast while the occurrence range of moderate drought,severe drought,and extreme drought decreased with time.The reduced intensity and *收稿日期：20230606 基金项目：国家重点研发计划项目（2021YFD1901104）*通讯作者：王靖，教授，主要从事农业生产系统模拟与气候变化影响评估研究，E-mail： 第一作者联系方式：赵凌暄，E-mail： 中 国 农 业 气 象 第 45 卷 420 frequency of drought is in favor of millet production in semi-arid region of Inner Mongolia and Hebei.Key words:Semi-arid region;Crop water deficit index;Millet;Growth stage;Drought intensity and frequency 干旱半干旱区约占全球陆地面积的 41%1，占中国陆地面积的 50%以上2，其生态环境较为脆弱、对气候变化敏感且适应能力弱3。蒙冀半干旱区位于中国华北地区，是典型的农牧交错区，蒸发量大、降水量少，雨养农业是其主要种植模式4。谷子（Setaria italica）作为典型的抗旱和环境友好型作物，因其耐旱、耐瘠薄的特点，成为干旱半干旱地区的重要特色作物5。谷子还具有耐储藏、营养价值丰富、粮草兼收的优点，在中国北方旱区农业生产中占有重要地位67。据统计，中国谷子播种面积约 100 万 hm2，年产量超过 200 万 t，单产与总产均居世界首位8；其中，内蒙古、河北和山西三省谷子种植总面积和总产量居全国前列，占中国谷子产量的 67.1%9，并形成了河北省张家口市、河北省承德市、内蒙古赤峰市的谷子优势特色产区1011。鉴于气候条件显著影响谷子生产，研究蒙冀半干旱区谷子干旱时空变化特征，有助于开展谷子种植区划并优化布局，应对日益严重的粮食危机，为中国北方干旱地区发展特色农业提供理论支撑，对农业供给侧改革和种植结构优化布局等具有重要意义。干旱胁迫对谷子生长影响研究较多集中于 CO2浓度变化、品种基因型等对谷子抗旱生理特性影响方面，如李君亮等12研究发现轻度干旱胁迫下，大气 CO2浓度的增加提高了叶片抗氧化酶活性，改变了叶片渗透调节物质含量，降低了干旱对谷子的不利影响；于国红等13研究不同基因型谷子品种抗旱性，揭示了谷子对干旱胁迫的调控机制。但极少有针对谷子种植区干旱时空变化特征的分析研究，缺乏干旱影响谷子生长发育过程的深入理解，制约了具有时空针对性的谷子抗旱措施发展。目前用于评价作物干旱程度的指标较多，其中作物水分亏缺指数（CWDI）不仅综合气象、土壤和作物三方面因素考虑了降水量和蒸散量对作物干旱的影响，并且考虑了作物生育期内干旱的累积效应，较为直观地反映作物水分亏缺程度14。张晓芳等1516利用CWDI 干旱指数分别分析了中国北方地区和南方地区春玉米生育期内的干旱时空变化特征；王连喜等1719利用 CWDI 干旱指数分析了河南省、江苏省及关中地区冬小麦生长季及需水关键期内的干旱变化特征；肖楠舒等20则利用 CWDI 对环渤海地区鲜食葡萄干旱风险进行了评估。以上结果表明，作物水分亏缺指数对不同地区、不同粮食作物及经济作物的干旱时空变化均具有较好的普适性。因此，为了揭示干旱对谷子生长发育过程的影响，了解蒙冀半干旱区谷子干旱应对的重点区域和重点时段，本研究应用 CWDI 指数，通过干旱强度、干旱频率和干旱站次比，分析蒙冀半干旱区谷子生长季（59 月）旬尺度干旱的时空演变特征，以期为谷子科学生产和节水灌溉提供理论指导，对蒙冀半干旱区谷子种植布局和防灾增收提供参考。1 材料与方法 1.1 研究区概况 研 究 区 域 位 于 华 北 地 区 蒙 冀 半 干 旱 区（38514531N，1095112042E），包括内蒙古自治区（包头市、呼和浩特市、乌兰察布市、锡林郭勒盟和赤峰市）、河北省（张家口市、承德市和保定市）以及山西省（大同市）（图 1）。该区属温带半湿润半干旱大陆性季风气候，冬季寒冷干燥，夏季炎热多雨、雨热同期，全年多大风天气，且风沙集中在春季21。研究区内年平均温度和年降水量空间分布差异大，年平均温度 213，年降水量 170 图 1 研究区域及气象站点分布示意图 Fig.1 The study area and the distribution of meteorological stations 第 4 期 赵凌暄等：基于 CWDI 蒙冀半干旱区近 60a 谷子干旱时空变化特征 421520mm，由东南向西北递减，年参考蒸散量 10002000mm22。1.2 数据来源 选取研究区内 27 个气象站点 19612019 年谷子生长季（5 月中旬9 月下旬）逐日气象观测资料，包括日平均温度（）、最高温度（）、最低温度（）、风速（ms1）、降水量（mm）、平均空气相对湿度（%）和日照时数（h）。数据来源于中国气象局国家气象信息中心地面气象资料观测数据集，该数据已经过严格的质量控制2324。参考前人研究结果2526，蒙冀半干旱区的谷子生长季（5 月中旬9 月下旬）在旬尺度上划分为 5个生育阶段，即苗期（5 月中旬6 月中旬）、拔节期（6 月下旬7 月中旬）、抽穗期（7 月下旬8 月中旬）、灌浆期（8 月下旬9 月上旬）和成熟期（9 月中旬下旬）。1.3 研究方法研究方法 1.3.1 作物水分亏缺指数（CWDI）作物水分亏缺指数（CWDI）是考虑连续五旬水分亏缺状况27，即当旬及向前推 4 旬的 CWDI 值按一定权重累加。CWDI 值越大，表示作物的水分亏缺越严重，即干旱程度越高15。计算式为 ii 1i 2i 3i 4CWDIaCWDIbCWDIcCWDI dCWDIeCWDI（1）式中，CWDIi、CWDIi1、CWDIi2、CWDIi3和CWDIi4分别为某时段当旬及向前推四旬的水分亏缺指数；a、b、c、d 和 e 分别为五旬 CWDI 的权重系数，参考黄晚华等2829研究结果，分别取值为 0.3、0.25、0.2、0.15 和 0.1。其中第 i 旬作物水分亏缺指数（CWDIi）计算式为 ciiciiciiciiETP100%ETPETCWDI0 ETP（2）式中，ETci为第 i 旬作物需水量（mm）；Pi为第i 旬降水量（mm）。作物需水量计算式为 cc0ETk ET （3）式中，kc为谷子作物系数，参考前人研究结果30，蒙冀半干旱区谷子生长季各月份作物系数如表1所示；ET0为参考作物蒸散量（mm），采用1998年联合国粮食及农业组织（FAO）推荐的Penman-Monteith公式3133 n2sa020.408RG900uee/T273ET 10.34u （4）式中，Rn为地表净辐射（MJm2d1）；G为土壤热通量（MJm2d1）；T为日平均气温（）；u2为2m高处风速（ms1）；es为饱和水汽压（kPa）；ea为实际水汽压（kPa）；为饱和水汽压曲线斜率（kPa1）；为干湿表常数（kPa1）。各项参数具体计算过程参考前人研究34。表 1 蒙冀半干旱区谷子生长季内各月作物系数 Table 1 Monthly crop coefficients during the growing season of millet 5 月 May 6 月 Jun 7 月 Jul 8 月 Aug 9 月 Sep 生长季平均 Mean kc 0.45 0.74 1.08 1.00 0.63 0.83 1.3.2 干旱等级划分 根据国家农业干旱等级现行标准（GB/T32136 2015）35及其他文献结果19,25，将蒙冀半干旱区谷子生长季内干旱划分为5个等级（表2）。表 2 干旱等级划分标准 Table 2 Division standard for millet drought grade 干旱等级 Drought grade CWDI 无旱 None drought CWDI35%轻旱 Slight drought 35%CWDI50%中旱 Moderate drought 50%CWDI65%重旱 Severe drought 65%CWDI80%特旱 Extreme drought 80%CWDI 1.3.3 干旱频率 干旱频率是评价某一等级干旱在谷子生长季内发生的频率高低，对某一站点，干旱频率用干旱发生次数和统计范围内总年数的比值来表示36，即 nF100N%（5）式中，F为干旱发生频率（%）；n为统计时段内某等级干旱发生次数；N为统计范围内总年数，本研究取值59（19612019年）。1.3.4 干旱站次比 干旱站次比反映某一等级干旱发生时，研究区内影响范围大小，即某等级干旱发生站点数与研究 中 国 农 业 气 象 第 45 卷 422 区内所有站点数的比值37，即 mP100M%（6）式中，P为干旱站次比（%）；m为某等级干旱发生站点的数量；M为研究区内气象站总数量（M=27）。2 结果与分析 2.1 区域干旱强度时空分布特征 2.1.1 时间变化 图2为19612019年蒙冀地区谷子生长季区域 图 2 19612019 年研究区域谷子生长季作物水分亏缺指数（CWDI）的年际变化 Fig.2 The change of growing season averaged crop water deficit index（CWDI）of millet from 1961 to 2019 in the study area CWDI的变化趋势。19612019年研究区谷子生长季平均CWDI在31%54%区间变化，最高值发生在1972年，最低值在2003年。从年代平均看，1970s和2000s干旱强度最高，1990s干旱强度最低。整体上，19612019年研究区谷子干旱强度呈不显著下降变化。图3为19612019年蒙冀半干旱区谷子生长季内旬尺度CWDI的变化。由图可知，研究区域内谷子生长季内水分亏缺呈先升高、后降低、再升高的变化规律，苗期后期（6月中旬）干旱强度最高，平均CWDI为57%，抽穗期（8月上旬）干旱强度最低，平均CWDI为39%。苗期（5月中旬6月中旬）谷子平均CWDI为47%，拔节期（6月下旬7月中旬）为49%，抽穗期为40%，灌浆期为41%，成熟期为47%。2.1.2 空间分布 由图4可见，19612019年研究区内谷子干旱强度呈西北向东南逐渐降低的规律。西北部常年处于高强度干旱（包含中旱、重旱和特旱）区域，东南部站点常年处于轻旱，甚至是无旱。随年份变化，蒙冀半干旱区干旱高强度范围有向西北方向收缩趋势，且强度有所降低；低强度范围由东南方向西北扩张；沿东北西南方中间地带逐渐由高强度地区转为低强度地区。随谷子生育进程，干旱高强度区域 图 3 19612019 年研究区域谷子生育阶段内 CWDI 旬尺度变化 Fig.3 Ten-day CWDI at each growth stage of millet from 1961 to 2019 in the study area 注：箱线图中上下两端端点表示 27 个气象站点累计 CWDI 的最大值和最小值；箱体表示 25%75%四分位；表示各生育期内 27 个气象站点的 CWDI 中位值；表示各生育期内 27 个气象站点的 CWDI 均值；短线表示标准误差。Note:E-is the first ten-day of a month；M-is the middle ten-day of a month；L-is the last ten-day of a month.The upper and lower lines of the box diagram represent the maximum and minimum values of accumulative CWDI for 27 meteorological stations.The box represents 25%75%interquartile range.is median of CWDI for 27 meteorological stations at each growth stage.is average value of CWDI for 27 meteorological stations at each growth stage.The short line represents the standard error.第4期 赵凌暄等：基于CWDI蒙冀半干旱区近60a谷子干旱时空变化特征 423 图 4 研究区域不同年代谷子各生育阶段 CWDI 的空间分布 Fig.4 Spatial distribution of CWDI at different periods and growth stages during millet growing season from 1961 to 2019 in the study area 注：图中 a、b、c、d、e、f 分别代表 19611970、19711980、19811990、19912000、20012010 和 20112019 年；1、2、3、4、5 分别代表谷子的苗期、拔节期、抽穗期、灌浆期和成熟期。Note:The figures a,b,c,d,e,and f represent the result during the period of 19611970,19711980,19811990,19912000,20012010 and 20112019,respectively;1,2,3,4,and 5 represent the seedling stage,jointing stage,heading stage,grain filling stage,and maturity stage of millet,respectively.的空间分布均表现为先收缩再扩张的趋势，且干旱强度基本稳定在中旱以上。2.2 区域干旱频率的时空分布特征 2.2.1 时间变化 5个干旱等级发生频率随谷子生育阶段的变化如图5所示，其中总旱发生频率为轻旱、中旱、重旱和特旱发生频率的总和。由图可见，谷子各生育阶段干旱发生频率均随干旱等级的升高而降低。苗期和拔节期干旱频率大致持平且为整个生育期内的最高水平，干旱频率分别为70%和71%；抽穗期和 中 国 农 业 气 象 第45卷 424 图 5 研究区域 19612019 年谷子各生育阶段不同等级干旱发生频率 Fig.5 Drought frequency of each drought grade during the growth stage of millet from 1961 to 2019 in the study area 注：短线表示标准误差。Note:The short line represents the standard error.灌浆期干旱频率相近且为生育期内的最低水平，干旱频率分别为55%和56%，成熟期干旱频率稍有升高但仍低于苗期和拔节期，干旱频率为67%。谷子各生育阶段无旱发生频率由高到低依次为抽穗期、灌浆期、成熟期、苗期和拔节期，发生频率分别为 45%、44%、33%、30%和29%。2.2.2 空间分布 由图6可见，19612019年蒙冀半干旱区干旱频率的空间分布特征与干旱强度一致，均呈西北向东南递减的规律。轻旱和中旱发生频率为西北、东南两侧高而中间低，轻旱发生频率为东南高于西北，中旱发生频率为西北高于东南。重旱和特旱发生频率表现为由西北向东南递减，重旱高频地区的范围明显大于特旱，但重旱的发生频率较特旱低20%。研究区西北部作为高强度干旱区，发生干旱多以等级较高的重旱和特旱为主，几乎不发生轻旱，而中部及东南部则以等级较低的轻旱和中旱为主。研究区西北部无旱发生频率低于30%，局部地区（如满都拉、二连浩特）低于10%，东南部除丰宁、承德、围场三站的无旱频率高于65%外，其余地区发生无旱的频率均介于30%60%。2.3 区域干旱站次比的时间变化特征 由图7可见，19612019年研究区内干旱站次比随年份呈波动性变化，且随干旱等级的升高而降低，由轻旱到特旱的多年平均站次比依次为23%、20%、图 6 研究区域 19612019 年谷子各等级干旱发生频率的空间分布 Fig.6 Spatial distribution of drought frequency in each drought grade of millet from 1961 to 2019 in the study area 第4期 赵凌暄等：基于CWDI蒙冀半干旱区近60a谷子干旱时空变化特征 425 图 7 19612019 年谷子各等级干旱站次比变化趋势 Fig.7 The change of the ratio of stations with different drought grades to total study station from 1961 to 2019 注：*表示 P0.05；y1，y2，y3，y4，y5和 y6分别代表无旱、轻旱、中旱、重旱、特旱和总旱的站次比趋势方程。Note:*represents P0.05;y1,y2,y3,y4,y5,and y6 represent the change of the ratio of stations with none drought,light drought,moderate drought,severe drought,extreme drought,and total drought to total study station,respectively.12%和5%。研究时段内，轻旱的干旱站次比随时间变化表现出显著上升趋势（P0.05），而中旱，重旱、特旱和总旱的站次比随时间变化均表现出显著下降趋势（P0.05），其中总旱的站次比为轻旱、中旱、重旱和特旱站次比的总和。表明该地区中旱、重旱、特旱的发生强度和发生面积均在减小，并且正在向轻旱甚至无旱转变。3 结论与讨论 3.1 结论结论（1）19612019年，蒙冀半干旱区谷子生育期内干旱强度呈先升高、后降低、再升高的变化趋势，总体呈下降趋势。苗期后期干旱强度达到最高，抽穗期干旱强度最低。干旱强度呈西北向东南逐渐降低的规律，随年代际干旱高强度范围向西北方向收缩，低强度范围由东南向西北扩张。（2）蒙冀半干旱区内各等级干旱发生频率随干旱等级升高而降低。谷子苗期和拔节期干旱频率最高，抽穗期和灌浆期干旱频率最低。总旱、重旱、特旱发生频率均为由西北向东南递减，轻旱和中旱的发生频率为西北、东南两侧高而中间地区低。西北部干旱多以重旱和特旱为主，几乎不发生轻旱，中部及东南部以轻旱和中旱为主，几乎不发生重旱和特旱。（3）19612019年，蒙冀半干旱区谷子生育期内干旱站次比随干旱等级升高而降低。研究区内干旱站次比即干旱发生范围显著减小（P0.05），中旱、重旱、特旱发生强度和发生面积均在减小，且向轻旱甚至无旱转变。3.2 讨论讨论 19612019年蒙冀半干旱区谷子生长季内的干旱强度呈显著下降趋势，这与目前中国干旱程度正在不断加深的结论相悖3841。但有研究表明，一方面，尽管全球变暖使大多数地区的气温均呈显著升高趋势，但蒸发皿蒸发量、潜在蒸散和作物需水量反而呈减少趋势42，由于影响环境水分变化的因素是多样的，所以该现象的机制目前尚不完全明确；另一方面，夏季风雨带的持续北移造成了中国北方和西部地区出现持续性多雨的格局43，使中国干旱半干旱区呈现较明显的暖湿化趋势4447。结合这两方面的原因，作物需水量的减少和降水量的增多会导致作物生长季的干旱减弱。另外，华北地区48、东北地区49的小麦50、玉米51、棉花52等作物的需水量也均呈下降趋势，从而导致生长季内干旱程度减弱，这与本研究结果一致。汪士为53基于温度植被干旱指数（TVDI）研究发现近20a内蒙古干旱程度有所缓解；陈敏等54基于气象干旱综合指数（MCI）研究发现19612016年中国北方不同等级干旱均存在“先增多、后减少”的年代际变化特征，与本研究结果基本一致。但也有一些学者的研究结果与之相悖，如张煦庭等55基于标准化降水蒸散指数（SPEI）研究发现，19602015年内蒙古东部地区干旱程度加深；杜波波等56基于标准化降水指数（SPI）研究发现内蒙古锡林郭勒盟在19802018年干旱频次、持续时间均呈增加趋势，说明SPI、SPEI等基于降水量单一气象要素的干旱指数无法考虑干旱的累积效应及其对农作物的影 中 国 农 业 气 象 第45卷 426 响，而CWDI、TVDI等基于综合气象因素的干旱指数能够充分反映干旱的连续性、累积性、区域差异性和对作物的影响程度。本研究表明，谷子生长期内的干旱强度呈先升高、后降低、再升高的趋势。干旱强度在苗期逐渐加剧，因为该阶段谷子需水量开始增加57，而此时降水量处于较低水平，因此干旱程度随着谷子生长而加剧，到拔节期时干旱程度开始降低，但此期降水量低于需水量，干旱程度仍较强。抽穗期时，虽谷子需水量仍然较高，但此时正值雨季，降水量大于谷子需水量，因此干旱得以缓解，干旱强度逐渐降低。生育后期，谷子叶片逐渐干枯脱落，叶片蒸腾减弱，谷子需水量整体趋于下降，但此时雨季已过，降水量水平较拔节期、抽穗期时大幅度降低58，因此干旱强度逐渐增大。结合干旱强度和干旱频率在谷子各生育阶段的分布特点，建议重点在谷子苗期和拔节期时进行适当补灌，以减少干旱对谷子生产的不利影响。另外，研究区内各站点在当地温度条件满足谷子正常生长所需时，可将播期延迟至5月中旬，以使谷子的需水关键期（抽穗期和灌浆期）与降水相匹配。研究区内西北部，如满都拉、二连浩特、阿巴嘎旗等地，因常年处于高强度干旱状态，不适宜谷子种植，应减少或摒弃谷子种植；中部如呼和浩特、赤峰、巴林左旗等地，除播种时进行一次灌溉外，还应在苗期和拔节期分别进行3050mm的补灌，水源充足的地区可在灌浆期再进行3050mm的补灌；在研究区东南部如承德、丰宁、围场等地，除播种时进行一次灌溉外，后期可完全雨养，无需补灌。此外，作物系数反映的是作物需水量多少的生理特性，尽管对同一种作物而言，品种、生育期、土壤环境、气象因子等因素不同，都会导致作物系数有所区别5960。但由于目前针对谷子在不同地区、不同田间管理措施下作物系数的研究较少。本研究中采用同一种作物系数计算研究区内27个站点的作物水分亏缺指数，可通过后续针对谷子作物系数的田间试验结果来改进。另外由于研究区内气象因素的空间差异，导致不同站点谷子的实际生育期存在一定差异，但本研究为了统一比较气象因子变化对谷子干旱的影响，因此仅在旬尺度上划分了研究区域内谷子的生育阶段，后续可根据谷子的实际生育期精细划分研究亚区，分别对其展开研究。参考文献 References 1 Huang J P,Yu H P,Dai A G,et al.Drylands face potential threat under 2C global warming targetJ.Nature Climate Change,2017,7(6):417-422.2 马柱国,符淙斌,杨庆,等.关于我国北方干旱化及其转折性变化J.大气科学,2018,42(4):951-961.Ma Z G,Fu Z B,Yang Q,et al.Drying trend in Northern China and its shift during 1951-2016J.Chinese Journal of Atmospheric Sciences,2018,42(4):951-961.(in Chinese)3 Fraser E D G,Dougill A J,Hubacek K,et al.Assessing vulnerability to climate change in dryland livelihood systems:conceptual challenges and interdisciplinary solutionsJ.Ecology and Society,2011,16(3):3.4 赵尊练,严小良.地中海地区节水高效设施农业及其对我国干旱半干旱地区农业发展的借鉴J.农业工程学报,2003(2):12-17.Zhao Z L,Yan X L.Water-saving and protected agriculture in the Mediterranean area and its reference to arid and semi-arid agriculture of ChinaJ.Transactions of the CSAE,2003(2):12-17.(in Chinese)5 杜虎平.黄土高原小杂粮生产发展与产业开发研究D.咸阳:西北农林科技大学,2005.Du H P.Study on production development and industrialization of minor crops in the Loess Plateau of ChinaD.Xianyang:Northwest Sci-tech University of Agriculture and Forestry,2005.(in Chinese)6 张雪峰.中国谷子产业发展问题研究D.哈尔滨:东北农业大学,2013.Zhang X F.Studies on the issues of millet industry development in ChinaD.Harbin:Northeast Agricultural University,2013.(in Chinese)7 刘猛,张新仕,王桂荣,等.河北省主要杂粮种植成本与收益分析J.河北农业科学,2020,24(1):1-4+25.Liu M,Zhang X S,Wang G R,et al.Analysis on planting cost and benefit of miscellaneous grain crops in Hebei provinceJ.Journal of Hebei Agricultural Sciences,2020,24(1):1-4+25.(in Chinese)8 弓俊.小杂粮产业发展存在的问题与对策J.农家参谋,2022,715(1):43-45.Gong J.Problems and countermeasures in the development of small grains industryJ.Farm Counsellor,2022,715(1):43-45.(in Chinese)第4期 赵凌暄等：基于CWDI蒙冀半干旱区近60a谷子干旱时空变化特征 4279 曲佳佳,张蕙杰,麻吉亮.中国杂粮生产及贸易形势展望J.农业展望,2021,17(5):78-85.Qu J J,Zhang H J,Ma J L.Prospect of production and trade of Chinas miscellaneous grainsJ.Agricultural Outlook,2021,17(5):78-85.(in Chinese)10 刘猛,刘斐,赵宇,等.河北省杂粮生产现状及生产布局展望J.农学学报,2021,11(3):93-100.Liu M,Liu F,Zhao Y,et al.Production status and production layout prospect of coarse grain in Hebei provinceJ.Agricultural Journal,2021,11(3):93-100.(in Chinese)11 王志春,陈素华,王松虎,等.赤峰地区谷子气候适宜度评价模型研究J.东北农业科学,2022,47(4):122-126.Wang Z C,Chen S H,Wang S H,et al.Study on climatic suitability evaluation model of millet in Chifeng districtJ.Journal of Northeast Agricultural Sciences,2022,47(4):122-126.(in Chinese)12 李君亮,王士博,李亚军,等.CO2浓度升高对干旱胁迫下谷子细胞结构和抗逆生理的影响J.应用生态学报,2023,34(5):1281-1289.Li J L,Wang S B,Li Y J,et al.Effects of elevated CO2 concentration on cell structure and stress resistance physiology of Setaria italic under drought stressJ.Chinese Journal of Applied Ecology,2023,34(5):1281-1289.(in Chinese)13 于国红,刘朋程,郝洪波,等.不同基因型谷子对干旱胁迫的调控机制J.植物营养与肥料学报,2022,28(1):157-167.Yu G H,Liu P C,Hao H B,et al.Regulation mechanism of drought resistance in different genotypes of foxtail milletJ.Journal of Plant Nutrition and Fertilizers,2022,28(1):157-167.(in Chinese)14 李崇瑞,游松财,武永峰,等.改进作物水分亏缺指数用于东北地区春玉米干旱灾变监测J.农业工程学报,2019,35(21):175-185.Li C R,You S C,Wu Y F,et al.Improved cr