全球气候变暖将会对温度和降雨产生影响，从而直接影响土壤的湿度^[1]，导致农业用水压力与日俱增。而我国农业生产用水占总用水量的80%^[2]，农业水资源短缺和灌溉水利用率低的问题一直存在，因此明确气候变化对农作物的灌溉需水量的影响很有必要。目前，气候变化对农作物灌溉需水量的影响受到许多学者的关注。如，Gunter Wriedt^[3] 运用作物模型(EPIC)及作物分布和种植面积的统计数据，研究了欧洲不同地区的净灌溉需水量，并根据水分运输效率和灌溉管理，模拟出总的灌水量是田间需水量的1.3—2.5倍；C.S. De Silva^[4]研究发现在A2和B2情景下，20世纪50年代斯里兰卡水稻灌溉需水量分别增加了23%和13%；王卫光等^[5]基于统计降尺度模型(SDSM)，研究了气候变化背景下长江中下游水稻水分需求变化。还有不少学者运用作物系数法和GIS来研究气候变化对作物水分的影响，如，J. A. Rodriguez Diaz^[6]对西班牙的瓜达尔基维尔河流域的作物灌溉需水量进行预测，发现在气候变化条件下，20世纪50年代该流域农业灌溉用水显著增加20%；韩冰等^[7]以辽宁营口灌区为例，分析了气候变化对水稻生育 期和灌溉需水量的影响，发现过去60a水稻生育期缩短了12.7 d。由于气温上升导致农田腾发量增加，该区域水稻灌溉需水量都呈上升趋势。同时，国内外学者在气候变化对作物物候的影响上做了大量的研究。Estrella,Xiao 和Wang jing^{[8, 9, 10]}研究发现气候变暖促使冬小麦品种更替，且缩短了冬小麦的生育期。国内学者^{[11, 12, 13, 14]}研究得出中国北方小麦物候期变化主要受温度影响，气候变暖导致播种期推迟，成熟期提前，从而导致生育期缩短。

华北平原是中国冬小麦主产区之一，气候变化必然会影响到冬小麦生育期和其对水分的需求，从而影响到其产量和品质。在华北地区普遍升温的背景下，冬小麦灌溉需水量区域变化研究是一个重要研究内容，影响冬小麦灌溉需水量的驱动因素也需进一步明确。因此，本文以华北区为例，利用1980—2010年4个典型站点(图 1)的冬小麦生育数据和气象资料，研究了冬小麦生育期的变化和灌溉需水量的主要影响因素及其规律，为合理配置灌溉计划和提高水资源的利用效率提供重要的参数，为冬小麦种植制定适应气候变化措施提供科学的依据。

图 1 华北平原典型气象站点分布图 Fig. 1 The geographical location of meteorological stations in the North China Plain

图选项

华北地区主要粮食作物有小麦、玉米、高粱、谷子、薯类等，是我国重要的粮棉产区，在研究区选取4个有代表性站点，分别是北京密云站、天津宝坻站、石家庄栾城站和邢台南宫站，各个站基本情况见表 1。

气象数据来源于中国气象科学共享网。选取数据记录较完整的4个代表性地面气象站1981—2010年逐日气象资料，其中包括平均气温(℃)、最高气温(℃)、最低气温(℃)、日照时数(h)、降水量(mm)、风速(m/s)、平均相对湿度(%)。对于气象数据的缺失，温度(平均、最高、最低温度)缺测值利用五日滑动平均法进行插补，降水量缺测值利用附近站点数据进行线性插补^[15]。

冬小麦生育期数据来源于北京密云站、天津宝坻站、石家庄栾城站、邢台南宫站1981—2010年的观测资料。因冬小麦生育期观测站和地面气象站(北京站、天津站、石家庄站、南宫站)并不完全一致，所以本文采用取就近匹配原则，选取了4个有连续的30a观测序列资料的农业气象站地与之对应。

本文将冬小麦的整个生育期划分为5个生育阶段：播种—出苗期，出苗—拔节期、拔节—抽穗期、抽穗—乳熟期、乳熟—成熟期。根据冬小麦多年发育期资料，采用生育期日数的多年平均值来代表当地一般生育期日数。生育期日数按照“儒略日(Julian Day)”计算，即日期按照每年365d转化为“日序day”，1月1日即为1，以此类推。

冬小麦生育阶段有效降水量即冬小麦某一生育阶段内，降水中实际补充到小麦根层土壤的净水量，代表总降水量中的有效部分。本文利用美国农业部土壤保持局推荐的有效降水量分析方法^{[16, 17, 18]}:

点击浏览公式

式中，Pe为有效降水量(mm/d)，P为总降水量(mm/d)。

冬小麦生育期需水量的确定采用FAO推荐的公式^[19]：

点击浏览公式

式中，ET_C为作物需水量(mm/d)，K_c为作物系数，ET₀为参考作物蒸散量(mm/d)。

参考作物蒸散量(ET₀)的计算采用FAO(1998)推荐的Penman-Monteith公式计算^[19]。

点击浏览公式

式中，ET₀为参考作物日蒸散量(mm/d)；R_n为作物表面净辐射量(MJ·m^-2·d^-1)；G为土壤热通量( MJ·m^-2·d^-1 )；Δ为饱和水汽压与温度关系曲线的斜率(kPa/℃)；γ为湿度计常数(kPa/℃)；T为平均气温(℃)；U₂为在地面以上2m高处的风速(m/s)；e_s为空气饱和水汽压(kPa)；e_a为空气实际水汽压(kPa)。

根据冬小麦生长发育特征，在FAO推荐值^[19]基础上，根据分段单值平均法计算生育中期的作物系数。将生育期划分为初始生长期即从播种到作物覆盖率接近10％，此阶段作物系数为K_cini=0.7(播种—出苗期，出苗—拔节期，其中越冬期为0.4)；快速发育期也是生育中期即从覆盖率10％到充分覆盖，此阶段作物系数为K_cmin=0.15(拔节—抽穗期，抽穗—乳熟期)；成熟期即从叶片开始变黄到生理成熟或收获，此阶段作物系数K_cend=0.4(乳熟期—成熟期)。根据当地气候条件，将生育中期作物系数进行修正，即：

点击浏览公式

式中，U₂为2m高度处的日平均风速(m/s)，RH_min为日最低相对湿度的平均值，h为作物的平均高度(m)。

灌溉需水量可以认为是生育期内总的作物需水量与生育期内总的有效降水量之差^{[18, 20, 21]}，用公式表达为：

点击浏览公式

式中，N为生育期日数(d)，ET_C为作物需水量(mm/d)，P_e为有效降水量(mm/d)，I_s为灌溉需水量(mm)。

作物灌溉需水量、生育期和气象要素的变化趋势采用气候趋势分析方法，气候倾向率由最小二乘法获得，计算公式为：

点击浏览公式

式中，X ^ _i为要素的拟合值；a为回归常数；b为回归系数，b×10称为气候倾向率，表示气候要素每10a的变化量；i为自然数列，代表年份序号；气候要素序列X ^ _i与自然数列i之间的相关系数称为气候趋势系数，其显著性检验本文采用相关系数检验法。

作物需水量和灌溉需水量的年际波动是规划和设计农田水利工程的重要依据，可用逐年作物灌溉需水量的标准差或离散系数来衡量其年际波动程度，但标准差与样本均值有关，当两组样本的均值存在显著差异时，用标准差便不能准确分辨哪组数据的离散程度更高，而离散系数指标既考虑了样本标准差大小又考虑了样本均值大小，便于两组数据离散程度的比较，采用离散系数来衡量作物灌溉需水量和生育期日数的年际波动程度，其计算公式为^[22]：

点击浏览公式

式中，Cv为离散系数；σ为标准差；μ为均值。Cv绝对值越大，表明年际波动越大，反之亦然。

采用趋势分析法对1981—2010年的4个典型站冬小麦生长季内各气象要素的变化进行了分析，总体上，冬小麦播种—出苗期有效降水量变化趋势较显著，其中北京密云站、天津宝坻站和石家庄栾城站呈增加趋势，而北京密云站通过a=0.01的显著性检验(表 2)。冬小麦其他生育阶段有效降水量基本呈增加趋势，但不存在显著性差异，其中北京密云站各生育阶段的有效降水量升幅范围(0.1—6.29mm/10a)高于其他站点。由此可见，华北平原冬小麦生长季内有效降水量年际变化趋势不一致，且年际变化波动较大，与前人研究结果一致^[23]。

已有的研究结果显示，过去50a华北区气温上升幅度较大的地区之一^[24]。冬小麦生育期阶段，除了播种—出苗期温度有小幅下降趋势以外，其他阶段温度都呈上升趋势，且各站点不同生育阶段温度的升幅范围为0.02—0.08℃/10a。4个站点的风速变化较不一致，石家庄栾城站冬小麦生育阶段风速都呈减小趋势，且通过了a=0.01的显著性检验，而其他站只有拔节—抽穗期风速呈减小趋势，其余生育阶段都呈增加趋势，但是趋势不显著。此外，石家庄栾城站和邢台南宫站的日照时数和相对湿度呈下降趋势。

由表 3可以看出，冬小麦的播种期和出苗期以及石家庄栾城站和邢台南宫站冬小麦成熟期均有推迟趋势，且北京密云站推迟幅度较大，为6.12—6.84d/10a，均通过了a=0.01的显著性检验；可以看出高纬度站点的变化趋势较大，且显著性较好。其他生育阶段生育期日数均表现为减小趋势，即生育阶段呈提前趋势，平均每10a提前0.71—3.87d。与前人的研究结果一致^{[25, 26, 27]}。从年际波动来看，整体变化较一致，拔节期离散系数最大，为0.03—0.06；而成熟期离散系数最小，为0.01—0.02。6个生育时期的年际波动程度依次为：拔节期>抽穗期>播种期=出苗期=乳熟期>成熟期。

由表 4可以看出，高纬度站点北京密云站和天津宝坻站的全生育期长度比其他两个站点长。而近30a来，4个站点的全生育期均表现为缩短，且北京密云站和天津宝坻站生育期缩短最为显著，均通过了a=0.01的显著性检验。

有研究发现，冬小麦播种期的变化与平均温度呈正相关，而与日照时数呈负相关^[25]。播种期推迟还可能和前茬玉米推迟收获有关，而拔节期提前的原因是日照时数减少，积温的增加导致冬小麦成熟期推迟^{[25, 26]}。气候变化使我国的冬小麦种植区向北扩展，且适宜种植的品种向减弱冬性方向演化^{[28, 29]}。由于近30年来，农作物的种植制度和品种更替的变化，对作物生育期的影响也是不可忽视的。

对1981—2010年北京密云站和石家庄栾城站冬小麦灌溉需水量变化的气候倾向率进行分析(表 5)，可以得出，在抽穗期之前，北京密云站和天津宝坻站所需的灌溉需水量较多，这一生育阶段灌溉需水量有从南向北增加的趋势；而从抽穗期到成熟期刚好相反。从4个站点的需水强度(日平均灌溉需水量)来看(图 2)，需水强度较大的生育阶段是拔节—抽穗期和抽穗—乳熟期，是其他生育阶段的3—6倍。从年际波动来看，南部站点的年际波动系数较大，但4个站点不同生育阶段的系数大小变化较一致。播种—出苗期的年际波动系数最大，北京密云站、天津宝坻站、石家庄栾城站和邢台南宫站分别为0.53，0.44，0.91和0.61；而最小的是播种—成熟期，分别是0.12，0.12，0.17和0.14。6个生育阶段的年际波动程度依次为：播种—出苗期>乳熟—成熟期>抽穗—乳熟期>拔节—抽穗期>出苗—拔节期>播种—成熟期。

图 2 典型站点不同生育阶段需水强度 Fig. 2 Water demand intensity in different growth stages for meteorological stations

图选项

华北地区冬小麦灌溉需水量在空间上从北向南、自东向西逐渐递减趋势；在时间上，东西部地区灌溉需水量变化趋势相反，东部地区呈逐渐增加趋势，而西部地区呈减小趋势。冬小麦各生育阶段的灌溉需水量随时间的变化不同，站点播种—出苗，拔节—抽穗阶段灌溉需水量均表现为不同程度的减少，但是趋势不显著，仅天津宝坻站通过了a=0.05的显著性检验；而出苗—拔节和乳熟—成熟期则表现为增加趋势。就整个生育期而言，西部站点(北京密云站和石家庄栾城站)有减少趋势，分别减少6.72mm/10a和8.30mm/10a，石家庄栾城站的变化幅度较大；而东部站点(天津宝坻站和邢台南宫站)的趋势正好相反，分别增加了2.60mm/10a和7.08mm/10a。但灌溉需水量变化均未通过显著性检验。

上述结果表明：灌溉需水量及变化趋势，都存在明显的区域差异。本文利用Pearson相关分析法分析灌溉需水量和各气象要素之间的相关关系(表 6)，发现两者之间存在着一定的相关关系，而且显著性较好。4个站点之间同一气象要素的变化趋势存在一定的差异，不同气象要素与灌溉需水量的相关关系也不同，其中有效降水量对灌溉需水量的变化最为明显，它直接参与灌溉需水量的计算，与灌溉需水量关系密不可分。灌溉需水量同有效降水量、相对湿度呈负相关，且相关关系极显著(a=0.01)，与日照时数、平均温度和风速呈显著正相关，生育期长度虽没有直接参与灌溉需水量的计算过程，但在考虑气候变化对灌溉需水量的影响时，考虑生育期长度的变化是很有必要的^[30]。北京密云站和天津宝坻站生育期长度同灌溉需水量存在微弱负相关关系，而石家庄栾城站和邢台南宫站存在显著负相关关系，整体趋势一致。说明考虑生育期长度变化，可以使灌溉需水量计算的结果更加精确。由于华北区的降雨主要集中6—9月，冬小麦生育期推迟，可以有效减少冬小麦后期的灌溉量。

为了深入研究4个典型站点灌溉需水量变化的主要原因，本文选用线性逐步回归的统计方法分析了灌溉需水量与主要气候要素之间的关系。其中，有效降水量和相对湿度的变化对其影响最为突出，两者对它得贡献率都为负的。4个站点线性回归方程的决定系数均达到了0.8以上。石家庄栾城站和邢台南宫站，除了播种—成熟期外，其他生育阶段的决定系数都高达0.9以上。

每个生育阶段灌溉需水量的变化的主要气象要素有不同，拔节—抽穗期最主要影响要素是相对湿度；抽穗—乳熟期，北京站和其他站有所区别，分别是平均温度和相对湿度，平均温度对它得贡献率是 正的；其他生育阶段都是有效降水量。其次是日照 时数、风速和生育期长度。过去的30a，生育期内日照时数减少和风速降低导致作物需水量呈减小趋势。而需水量是灌溉需水量计算中的主要因子之一，必然导致灌溉需水量也会减少。

(1)过去30a，冬小麦播种期和出苗期均有推迟趋势，且北京密云站推迟幅度较大，为6.12—6.84d/10a，可以看出高纬度地区的变化趋势显著。其他生育阶段生育期日数均表现为减小趋势，即生育阶段呈提前趋势。

(2)华北地区冬小麦灌溉需水量在空间上北多南少，有自东向西递减趋势；在时间上，东西部地区灌溉需水量变化趋势相反，东部地区呈逐渐增加趋势，而西部地区呈减小趋势。冬小麦各生育阶段的灌溉需水量随时间趋势不同，播种—出苗、拔节—抽穗、抽穗—乳熟阶段灌溉需水量均表现为不同程度的减少趋势，而出苗—拔节和乳熟—成熟期则表现为增加趋势。就整个生育期而言，西部站点(北京密云站和石家庄栾城站)有减少趋势，分别减少6.72mm/10a和8.3mm/10a；而东部站点(天津宝坻站和邢台南宫站)的趋势正好相反，分别增加2.6mm/10a和7.08mm/10a。在抽穗期之前，北京密云站和天津宝坻站所需的灌溉需水量较多，这一生育阶段灌溉需水量有从南向北增加的趋势；而从抽穗期到成熟期刚好相反。从四个站点的日平均灌溉需水量来看，需水强度较大的生育阶段相同，是拔节—抽穗期和抽穗—乳熟期。6个生育阶段灌溉需水量的年际波动程度依次为：播种—出苗期>乳熟—成熟期>抽穗—乳熟期>拔节—抽穗期>出苗—拔节期>播种—成熟期。

(3)冬小麦灌溉需水量同有效降水量、相对湿度呈负相关，且相关关系极显著，与生育期长度存在微负相关关系；与日照时数、平均温度和风速呈显著正相关。若考虑生育期长度变化，可以使灌溉需水量计算的结果更加精确。每个生育阶段灌溉需水量的变化的主要气象影响要素不同，拔节—抽穗期最主要影响要素是相对湿度；抽穗—乳熟期，北京密云站和石家庄栾城站有所区别，分别是平均温度和相对湿度，平均温度对它得贡献率是正的；其他生育阶段都是有效降水量，其次是日照时数、风速和生育期长度。

影响灌溉需水总量的因素较复杂，前人研究的影响要素有作物播种面积^{[3, 31]}、种植制度、海拔和纬度^[32]等。而本文从气象要素的角度考虑，明确冬小麦不同生育期的灌溉需水量及其影响因素。本文研究得出，有效降水量是影响冬小麦生育期内灌溉需水量的最主要指标。华北区站点近30年冬小麦生育期内的有效降水量表现为增加趋势，与前人研究的近50年华北区域年降水趋势不同^[33]，主要是因为20世纪90年代以来冬小麦关键生育期的降水相对有所增加^{[34, 35]}。同时由表 6看出，冬小麦关键生育期的灌溉需水量主要受相对湿度影响，与梁丽乔、杜加强等^{[36, 37]}研究松嫩平原及黄河上游生长季参考作物蒸散发对相对湿度最为敏感的结果相一致。其原因是这一时段有效降水虽有一定幅度的增加，但并不能完全满足冬小麦的生长需求。由于相对湿度的显著增加造成土壤湿度增加，同时降低田间温度，减小蒸发，对田间节水起到了持续保墒的关键性作用。也有学者研究表明作物水分变化主要受日照和风速^{[38, 39, 40]}影响。主要是由于日照减少使到达地面的能量减弱，导致地面蒸发减少；风速降低会使空气与土壤中的水分交换强度变低，有利于土壤水分的保持。大部分前人研究是基于作物全生育期，而本文细化到作物的每个生育阶段，同时结合气象要素和生育期来研究作物水分的内在变化影响因素，能使结果更加全面准确。由于数据搜集有限，本文结果在其它地区的适用性较弱，有待进一步研究。运用灌溉需水量的结果与农业实际生产、农机措施相结合在节水的前提下提高作物产量将是下一步的研究重点。

在明确冬小麦生育期和灌溉需水量变化的基础上，有利于采用针对性措施，有效应对气候变化和促进小麦生产。如非充分灌溉技术，做到及时在小麦关键需水期补充供水，满足田间持水量下限55—65%的灌溉量^[32]，可以有效提高水分利用效率和产量^[41]；通过合理施肥可以使土壤保持较高的含水量，产生水肥耦合效应^{[42, 43, 44]}；选用晚熟品种延长冬小麦生育期，提高生物量累积，从而达到提高产量的目的。