疏勒河流域中游地区玉米性状及单产的影响因子研究

刘焕才， 史书琦， 李 曼， 张艳芳， 韩 丽

（1.山西师范大学地理科学学院，山西 太原 030031；2.陕西师范大学地理科学与旅游学院，陕西 西安 710062）

粮食安全是国家安全的基础，是经济发展、社会稳定的基石［1］。玉米作为我国主要粮食作物之一，2020 年国内玉米产量达到2.66×109t，约占全国粮食产量的38.93%。随着粮食需求量增长、耕地面积的刚性减少，提高玉米单产成为保障我国粮食安全的重要途径之一［2］。相关学者针对玉米单产的影响因素研究主要涉及自然因素和人为因素两方面［3-4］。在自然因素方面，气候变化背景下，干旱、强降水和高温等因素对玉米性状、单产影响显著［5］，且表现出明显地域差异，例如：李祎君［6］利用通径分析法探讨气象灾害对东北平原玉米单产及品质的影响，提出干旱已成为威胁玉米气候产量的首要灾害；江铭诺［7］等利用世界粮食作物研究模型（WOFOST）、奇异值分解（SVD）等方法分析华北平原不同区域夏玉米潜在产量，发现该区域夏玉米潜在产量减少与太阳总辐射下降有关；陈上［8］基于美国的农业技术转移决策支持系统（DSSAT）和机器学习算法预测黄土高原玉米产量，认为降水的年内分配不均和频发的干旱是限制玉米产量的主要因素；Chen等［9］基于DSSATCERES-Maize模型和极端气候指数探讨长江流域雨养农业区玉米单产对气候变化的响应，指出温度与玉米单产的相关性最高。在人为因素方面，合理的灌溉、施肥和播期等因素是玉米增产的重要途径，例如：Meng 等［10］通过探讨净光合速率、气孔导度以及蒸腾速率变化与灌水量的关系，发现适宜的灌水量有利于提高玉米干物质积累量和单产；吴晓丽等通过研究不同施肥模式下玉米各器官对氮、锌的吸收和累积，得出合理的施肥量和施肥技术可提高肥料利用率，进而促进玉米生长发育的结论［11-13］；薛庆禹等［14］通过设置不同播期下夏玉米的试验，发现夏玉米生育期随播期推迟而缩短。

目前，众多学者普遍选用各类作物模型研究作物生长发育过程及产量［15］，例如：澳大利亚的农业生产系统模拟器（APSIM）、美国的农业技术转移决策支持系统（DSSAT）、联合国粮食及农业组织发布的作物-水生产力模型（AquaCrop）、荷兰瓦赫宁根大学与国际水稻研究所联合开发的作物生长模拟系统（ORYZA）。其中，DSSAT模型作为全球应用最广泛的作物模型之一，能够在大量数据基础上综合考虑作物、环境和管理措施三者间的相互关系［16］，对特定区域作物品种的遗传参数进行本地化调试。在我国，多地已经利用该模型开展了相关研究工作，并取得良好适用效果，如：韩智博等［17］利用DSSAT 模型开展未来气候变化对黑河流域玉米产量的影响研究，并提出玉米的适应性措施；王兴鹏等［18］基于DSSAT 模型开展南疆地区膜下滴灌种植模式下的棉花灌溉定额研究，为播期及生育期灌溉管理提供借鉴。

西北干旱区土地资源丰富，光热充足，是我国粮食生产的战略后备区，但以干旱缺水为代表的众多因素长期制约着当地农业生产［19］。地处疏勒河流域中游地区的甘肃省玉门市作为我国西北干旱区的典型农业灌溉区和重要玉米产区，目前的研究多集中于气候变化对作物产量的影响，而将自然和人为因素对玉米性状及气候产量影响定量化的综合研究较少［20-22］。自然因素主要通过最高温、最低温、太阳辐射和CO2等影响因子对粮食作物的生长发育和结实率产生重要影响［5-9］；另一方面，自然因素对粮食生长发育的作用会受到农户田间管理行为的影响［23-24］，即人为因素，例如：价格支持政策、包膜控释肥、喷滴灌、调整农时、补种补苗和覆盖棚膜等，可见自然因素和人为因素会共同作用于玉米性状和产量。因此，基于数据可获取性，本文选取最高温、最低温、太阳辐射、降水、风速等自然因素和化肥施用量、有效灌溉面积等人为因素，利用DSSAT-CERES-Maize 模型定量评估其对玉米性状及气候单产的影响，以此探究当地玉米的适宜生长条件，为疏勒河流域乃至西北干旱区保障粮食安全提供科学支撑。

1 材料与方法

1.1 研究区概况

甘肃省玉门市（39°40′～41°00′N，96°15′～98°30′E）位于疏勒河流域中游、河西走廊西端。该地光照充足，热量丰沛，但降水偏少，属大陆性中温带干旱气候［25］。玉米作为当地主要粮食作物之一，播种面积约占粮食作物总播种面积的26.76%，研究区玉米物候期尺度下的气象因子均值为：最高温29.01 ℃，最低温13.52 ℃，风速3.37 m·s-1，降水量56.85 mm。受区域降水少、年内年际变化大、蒸发大等因素影响，玉米生长用水高度依赖灌溉，疏勒河是其灌溉用水的主要水源。

1.2 数据来源

本文选取的数据类型有气象数据、土壤属性数据和田间管理数据。气象数据来源于全国温室共享平台，包括1990—2020 年玉门逐日气象数据（最高温、最低温、风速、降水量和日照时数）。土壤属性数据来源于世界土壤数据库（Harmonized World Soil Database）的中国土壤数据集（http://vdb.soil.csdb.cn/），包括玉门地区的土壤质地组成、有机质含量、容重和土层全氮含量等。田间管理数据涉及生长阶段、施肥和灌溉等要素的技术信息，主要通过实地调研获取（表1）。

表1 田间管理要素技术信息Tab.1 Technical information of field management elements

1.3 研究方法

1.3.1 太阳辐射太阳辐射对玉米的光合作用有重要作用，是作物生长模型的关键输入数据，其演变与日序时次的日照时长存在较高一致性［26］，因此本文基于日照时数利用埃斯屈郎经验公式［27］计算太阳辐射。计算公式如下：

式中：δ为太阳倾角（rad）；J为日序；dr为日地相对距离；Ф为纬度（rad）；Ws为日落时角（rad）；Rmax为地球外辐射（MJ·m-2）；Rs为太阳辐射（MJ·m-2）；n为实际日照时数（h）；N为最大可能日照时数（h）；n/N为相对日照；as、bs为回归常数（as=0.25、bs=0.5），当阴天（n=0）时，到达地球表面的地球外辐射的透过系数为as；当晴天（n=N）时，到达地球表面的地球外辐射透过系数为as+bs。

1.3.2 非参数Mann-Kendall检验法目前，长期趋势分析包括增长和下降趋势的确定，研究表明，在分析趋势时，由于变量总体分布无法精确，非参数统计方法是效果较好的常用方法之一［28］。非参数Mann-Kendall检验法［29］（M-K检验）是计算序列数据长期趋势的非参数统计方法，该方法无需数据服从标准正态分布［30］，可有效消除数据相关性和异常值对检验结果的影响［31］，以适用范围广、人为性少、定量化程度高而著称［32］。具体步骤如下：

设X1、X2、X3、…、Xn为时间序列变量，n为时间序列长度，检验统计量（UFK）公式为：

式中：S为检验的统计量；VAR(S)为统计量S的方差；Xi、Xj分别为i、j年相应的测量值，i＞j，UF1=0。按时间序列X的逆序Xn，Xn-1，Xn-2，…,X1，重复上述计算，得到UBK，即：UBK=-UFK（K=n,n-1, …, 1），UB1=0。若UFK＞0，表示序列呈上升趋势；若UFK＜0，表示呈下降趋势；当UFK或UBK≥U0.05=±1.96 时，表示通过了0.05 置信度显著性检验；当UFK或UBK≥U0.01=±2.56时，表示通过了0.01置信度显著性检验；UFK=UBK的交叉点为突变点。

1.3.3 玉米气候单产的分离方法玉米单产与农田种植技术、土壤肥力和病虫害等因素密切相关，为分析气象因子与单产的关系，本文利用5 a滑动平均法［32］分离趋势单产与气候单产。对于样本容量为n的单产序列，5 a滑动平均序列表示为：

设立科创板，是将资本市场进一步细分，可促进A股市场分化，主板市场上不去，海外市场又太拥挤，一些海外的中概股公司要回归国内资本市场，也没有完全对应的承接。科创板将提供一个新的资本市场入口。同时有利于克服相当一部分高科技上市公司（即俗称“独角兽”）估值偏低的弊端。实际上，高科技企业才是推动股市繁荣的根本力量。单独设立科创板，将充分肯定其客观价值，推动上市公司脱虚向实，回归资本市场的本质属性，客观反映经济，扩大直接融资，刺激科技创新，利好实体经济。

式中：Y(w)为玉米气候单产（kg·hm-2）；Y为实际单产（kg·hm-2）；Y(t)为趋势单产（kg·hm-2），Yj(t)为第j年趋势产量滑动平均值（kg·hm-2）；k为滑动平均时间步长，取k=5；xj-k+1为滑动步长年间第（j-k+1）年的实际单产（kg·hm-2）。

1.3.4 敏感度分析本文将影响因子数据进行无量纲化处理后，分析玉米性状、气候单产对其敏感度［33］，公式如下：

式中：Pheno 为玉米性状或气候单产；Tmax为最高温（℃）；Tmin为最低温（℃）；SRAD为太阳辐射（MJ·m-2）；PRCP 为降水量（mm）；WIND 为风速（m·s-1）；Fer 为施肥量（t）；Irr 为有效灌溉面积（hm2）；β0为截距项；β1、β2、β3、β4、β5、β6和β7分别为气象单产或性状对最高温、最低温、太阳辐射、降水量、风速、施肥量和有效灌溉面积的敏感度。若β＞0，说明影响因子对玉米气象单产或性状表现出正向敏感度，为促进效应；若β＜0，则说明影响因子对玉米气象单产或性状表现出负向敏感度，为抑制效应；且|β|越大表示气象玉米单产和性状对影响因子越敏感。基于本文中玉米气象单产及性状对影响因子的敏感度变化范围为-0.85～0.85，为方便比较，将敏感度进行A～-A 分级排序，共计17 级，相邻两级间相差0.1，由大到小依次为：A＞B＞C＞D＞E＞F＞G＞H＞I＞-H＞-G＞-F＞-E＞-D＞-C＞-B＞-A，其中I级敏感度的变化范围为-0.05～0.05。

1.3.5 CERES-Maize模型CERES-Maize 是DSSAT模型中的子模型，通过以日为时间步长模拟玉米在不同种植模式下性状、单产动态形成过程。模型输入数据包括气象数据、土壤数据、作物品种遗传参数以及农田管理实践数据等，其中受数据可获得性影响，太阳辐射数据需根据日照时数估算获得；玉米品种特性遗传参数通过最大似然估计模块GLUE模拟获得，包括：P1（从出苗到幼年结束所需的积温）、P2（光周期敏感参数）、P5（灌浆特性参数）、G2（单株最大穗粒数）、G3（最大灌浆速率参数）和PHINT（完成一片叶生长所需参数）。本文中，模型的输出数据涉及包括玉米单产在内的17 种玉米性状，根据各类性状与影响因子的相关程度对相关性状进行分类（表2），其中，收获时单产和籽粒氮是反映玉米生产系统质量的核心性状［34］；影响玉米物质积累的各生长阶段物积速率和时长是玉米籽粒氮、单产的反映性状［35-39］。

表2 影响因子及其相关玉米性状Tab.2 Influencing factors and related maize traits

作物参数调试运转过程中，本文利用“试错法”对玉门市的玉米参数值进行3 次调试［40-41］，每次调试采用30000 次运转，使其模拟单产等无限趋近于实际单产。使用归一化均方根差（nRMSE）来度量模拟值与实测值的相对差异程度，并用一致性指数（D）检验模拟值和实测值的吻合度［42］，公式如下：

式中：Pi和Oi分别为模拟值和观测值；Oˉ为观测平均值；i为序列值；n为样本数。若nRMSE 值低于10%时，则模拟值与观测值的相对误差小，模拟精度高；若nRMSE 值介于10%～20%，则模拟效果较好；若nRMSE 值介于20%～30%，则模拟效果一般；若nRMSE 值高于30%，则表明模拟效果差。D值越接近1，表明模拟值与观测值的一致性越高，模型的模拟效果越好，反之亦然。

通过1990—2020 年玉米实际单产与模拟单产的对比分析发现（图1），nRMSE 值为0.22%，D值为0.99，说明其品种参数对玉米单产的模拟精度高，可用于研究区玉米生产潜力模拟研究。

图1 单产的模拟值与实测值拟合结果Fig.1 Fitting results of simulated and measured yield per unit area

2 结果与分析

2.1 玉米性状及单产对影响因子的敏感度

图2 1990—2020年玉门市玉米各生长阶段气象因子的Mann-Kendall检验Fig.2 Mann-Kendall of meteorological factors at each growth stage of maize in Yumen City from 1990 to 2020

2.1.2 水肥时间变化特征玉门市玉米化肥施用量变化趋势可分为2 个阶段（图3），第一阶段为1990—2012年，玉米化肥施用量随着年份增加呈现出极显著增加趋势，用量增加了259.00%；第二阶段为2013—2020 年，化肥施用量呈现出减少趋势，用量减少了29.79%。玉米有效灌溉面积总体上呈现出增加趋势，其中1990—1991、2005—2009、2017—2020年增速较大，分别为118.65 hm2·a-1、57.46 hm2·a-1、42.99 hm2·a-1，其余年份增速则较为缓慢。

图3 1990—2020年玉门市玉米化肥施用量和有效灌溉面积Fig.3 Fertilizer application amount and effective irrigated area of maize in Yumen City from 1990 to 2020

2.1.3 玉米单产时间变化特征将玉米的实际单产分离为趋势单产和气候单产（图4），1990—2020 年玉米实际单产多年平均值为10078.75 kg·hm-2，总体呈小幅增加趋势。趋势单产呈现出增加趋势，在20世纪90 年代，当地政府快速推广优质品种、制定农业生产条件管理办法，旱地农田建设有一定成效，因此单产增速较快；2000—2005年由于水资源利用效率较低，玉米单产呈现出减少趋势；进入2006 年后，随着玉米制种品种更新换代和农田水利建设有一定成效［43］，单产变化较小，进入稳定状态。气候单产多年来平均值为190.35 kg·hm-2，呈现出不显著减少趋势，但波动较大，变化范围为-2071.64～1158.21 kg·hm-2，其中1992 年和2019 年为丰产年，2002年和2020年是较为明显的低产年。

图4 玉米实际单产、趋势单产以及气候单产年际变化Fig.4 Interannual variations of actual yield per unit area,trend yield per unit area and climate yield per unit area

2.1.4 玉米性状及气候单产对影响因子的敏感度玉米不同性状及气候单产对影响因子响应程度存在差异，降水表现尤为明显，籽粒氮和气候单产等对其有较高敏感度，但收获时器官占比和全生育期蒸腾量对其敏感度都较低；最高温是对性状影响程度最大的影响因子，其中粒数、籽粒氮和气候单产都表现为较高正向敏感度（表3）。

表3 物候期玉米性状及气候单产对影响因子的敏感度Tab.3 Sensitivity of maize traits and climatic yield per unit area to meteorological influencing factors at phenological period

分别选取敏感度排序靠前的3个影响因子作为玉米性状及气候单产的关键影响因子，分析在玉米各生长阶段，气候单产及性状对关键影响因子敏感度可知（表4）。仅气候单产而言，对最高温和降水表现出强正向敏感度，对最低温表现出强负向敏感度；玉米性状则对穗期花粒期最低温、化肥施用量和有效灌溉面积表现为强正向敏感度，穗期、花粒期最高温表现为强负向敏感度。在不同生长阶段，同一性状对相同影响因子响应程度也存在差异，其中穗期气象因子敏感度最高，因此穗期是对玉米生长发育造成影响的关键时期。

表4 各生长阶段玉米性状及气候单产对影响因子敏感度Tab.4 Sensitivity of maize traits and climatic yield per unit area to influencing factors at each growth stage

2.2 玉米性状及气候单产对影响因子的敏感性试验

为进一步定量化探究最高温、最低温、太阳辐射、风速和水肥等影响因子对玉米性状及气候单产的影响程度，选择气候条件、玉米单产最接近多年平均值的2017 年为典型年份开展相关敏感性模拟试验，影响因子的参数设置参考当地实际情况（表5）。

表5 影响因子参数设置Tab.5 Parameter setting of influencing factors

2.2.1 玉米性状及气候单产对温度因子的敏感性低温和高温均会影响玉米生长发育，当温度变化＜-1.28 ℃时，低温将会导致花粒期物积速率降低，苗穗期时长延长，根、果重占比降低，茎重占比增加，最终引起气候单产降低；当温度变化＞1.12 ℃时，高温延长了苗穗期时长，玉米果重占比降低，根、茎、叶重占比增加，使得籽粒氮和气候单产表现出降低趋势（图5）。因此，玉门市玉米的适宜温度变化范围为-1.28～1.12 ℃，即：最高温变化范围为14.80～38.72 ℃，最低温变化范围为-0.38～22.32 ℃，在此范围内，温度升高将加快玉米苗穗期物积速率，缩短苗穗期时长，增加果重占比，使得玉米气候单产和籽粒氮增加。

图5 温度变化时玉米性状及气候单产变化率Fig.5 Change rate of maize traits and climatic yield per unit area under temperature change

2.2.2 玉米性状及气候单产对太阳辐射因子的敏感性当太阳辐射变化＞-0.15 MJ·m-2时，随太阳辐射的增强，苗期物积速率加快、时长延长，穗期物积速率无明显表现但时长缩短，花粒期物积速率和时长均表现为增加趋势（图6），在此基础上，在太阳辐射变化＞0.30 MJ·m-2时，随太阳辐射的增强，粒重/叶面积、粒数/叶面积、籽粒氮和气候单产均表现为增加趋势。因此，玉门市玉米的适宜太阳辐射变化幅度为：0.30～0.50 MJ·m-2，即：太阳辐射变化范围为3.93～25.15 MJ·m-2，在此范围内，太阳辐射增强，花前干物质积累量增加，花后玉米生长发育运输速率加快，使得玉米籽粒氮和气候单产上升。

图6 太阳辐射变化时玉米性状及气候单产变化率Fig.6 Change rate of maize traits and climatic yield per unit area under the change of solar radiation

当太阳辐射变化＜-0.15 MJ·m-2时，受弱光条件影响，玉米穗期物积速率降低、时长增加，导致花粒期物积速率和时长、粒重/叶面积、粒数/叶面积、籽粒氮、气候单产均表现为减少趋势。

2.2.3 玉米性状及气候单产对风速因子的敏感性风速过大或过小均不利于玉米生长发育。当风速变化≤-0.39 m·s-1时，随风速增大，尽管苗穗期物积速率呈现上升趋势，但花粒期物积速率和时长均表现为减少趋势，进而引起玉米蒸腾量和籽粒氮减少（图7）；当风速变化≥0.30 m·s-1时，随风速增大，苗穗期时长缩短，花粒期时长延长，籽粒氮上升，但全生育期蒸腾量过度下降，物质积累量较少，气候单产减少。

当风速变化介于-0.39～0.30 m·s-1时，随风速增大，苗穗期物积速率表现为波动趋势，时长为减少趋势；蒸腾量由增加转变为减少趋势；花粒期物积速率和时长、气候单产、籽粒氮均表现为增加趋势。因此，研究区玉米的适宜风速变化为：-0.39～0.30 m·s-1，即：风速变化范围为0～3.81 m·s-1，随风速增大，花粒期物积速率上升，时长延长，籽粒氮和气候单产上升。

2.2.4 玉米性状及气候单产对水肥因子的敏感性玉米的气候单产、籽粒氮、籽粒数和籽粒重均随水肥施入量的增加呈现出“中高周低”的分布特征（图8），按各性状数值大小，将水肥处理分为5级：（1）W1F1；（2）W0F1、W1F0、W1F2、W2F1；（3）W0F2、W2F0；（4）W2F2；（5）W0F0。当水肥投入量处于W1F1（灌溉15 mm，施肥10 kg·hm-2）时，玉米各性状数值达到最大，即：中等程度的灌溉和施肥有助于增产。但无水、无肥、高水或高肥均会引起玉米不同程度减产，其中，当灌溉量处于W1或施肥量处于F1时，减少幅度最小；当灌溉施肥量处于W0F0（灌溉0 mm，施肥0 kg·hm-2）时，减少幅度最大。尽管合理灌溉施肥条件对玉米性状有明显的促进作用，但不同性状对水肥的响应程度不同，由大到小排序为气候单产＞籽粒数＞籽粒氮＞籽粒重。

图8 水肥变化时玉米性状及气候单产变化Fig.8 Changes of maize traits and climatic yield per unit area under the change of water and fertilizer

3 讨论

最高温、最低温、降水等自然因素和化肥施用量、有效灌溉面积等人为因素在玉米生长发育中起着关键作用［44］，但由于地理环境差异，不同地区玉米对主要影响因子的响应形式存在差异［45］。本文指出地处西北地区的疏勒河中游地区玉米性状和气候单产对苗期最高温表现为正向敏感度，这是由于苗期最高温并未超过幼苗适宜生长的上限温度35～40 ℃［46］，因此，最高温的适当升高可以提高该地区玉米幼苗发芽率、根重占比和物质积累速率［47］；相反，近31 a 来，当地苗期最低温的下限温度范围为-1.5～7.5 ℃，尽管亦呈现一定程度的增加趋势，但总体仍低于玉米苗期最低温的适宜下限温度6～10 ℃［46］，不利于玉米增产。由于玉米属于喜温作物，穗期、花粒期最低温适当的增加可以加快物质积累速率，从而提高玉米单产［48］；相反，在降水稀少的西北内陆地区，穗期、花粒期极端高温会加快水分蒸发速率，缩短生育期长度，使玉米减产［49］。除此之外，降水、有效灌溉面积和化肥施用量在全生育期对玉米性状和气候单产均起到正向促进作用，其中玉米是生育期内需水较多的旱作作物［50］，特别是对于西北干旱区而言，有效降水的增多可以大大促进干物质积累、植株吐丝和授粉，降低水胁迫产生的可能性［51］，同时合理灌溉也可以缓解降水不足对玉米单产的负面影响［52］，促进玉米生长发育。2013年以来玉门市化肥施用量显著降低，但由于农田水利建设的有效推进，可提高化肥利用效率，进而释放玉米高产潜力［53］。

由此可见，该地区极端恶劣的自然条件使得当地玉米性状和气候单产对众多自然因素和人为因素都存在极高敏感度，表现出一定独特性。因此，为保障疏勒河中游地区的粮食安全，根据本文研究，该地区需提高应对极端高温或极端降水等事件的能力，同时建议高效利用疏勒河流域的热量资源，通过地表覆膜、调整播期、土地秋翻等人工保温措施，规避低温胁迫，降低生产风险；通过推行节水灌溉、秸秆还田等农田管理措施，提高水资源利用效率，改善土壤肥力，提高玉米生产潜力，发展与气候变化相适应的优质农业。

4 结论

本文在研究我国疏勒河中游地区玉米种植对气象因子与水肥因子的敏感度基础上，利用DSSATCERES-Maize模型定量评估主要影响因子对玉米性状和气候单产的影响，得到了如下结论：

（1）1990—2020年地处疏勒河中游地区的玉门市最高温、最低温和降水量变化呈上升趋势，太阳辐射和风速呈下降趋势，气候暖湿化特征显著。玉米化肥施用量1990—2012 年呈现出增加趋势，2013—2020年呈现出减少的变化趋势，有效灌溉面积则呈现出持续增加趋势。

（2）玉米实际单产总体呈现出小幅增加趋势，其中趋势单产受政府快速推广优质品种和农田水利建设等政策的影响，呈现出增加-减少-稳定的变化趋势；玉米气候单产总体上则呈现出减少趋势，但波动较大，敏感度分析结果表明，其对最高温和降水表现出强正向敏感度，对最低温表现出强负向敏感度；玉米性状的敏感度分析结果表明，其对化肥施用量、有效灌溉面积和穗期、花粒期的最低温表现为强正向敏感度，对穗期、花粒期最高温表现为强负向敏感度。

（3）通过主要影响因子对玉米性状、气候单产影响模拟试验研究发现，最高温、最低温、风速和水肥因子与性状、气候单产存在“倒U型”关系，即：随着上述影响因子数值的增大，其对玉米性状、气候单产的影响由促进效应转变为抑制效应；太阳辐射与性状、气候单产存在“U 型”关系，即：随着太阳辐射增强，对玉米性状和气候单产的影响由抑制效应转变为促进效应。疏勒河中游地区玉米植株生长发育的适宜条件为：最高温（14.80～38.72 ℃）、最低温（-0.38～22.32 ℃）、太阳辐射（3.93～25.15 MJ·m-2）、风速（0～3.81 m·s-1）以及实际降水基础上的灌溉量（15 mm）和施肥量（10 kg·hm-2）。