河套灌区春小麦生产水足迹影响因子敏感性及贡献率分析

孙世坤，刘文艳，刘 静，王玉宝，陈帝伊，吴普特

（1西北农林科技大学旱区农业水土工程教育部重点实验室，陕西杨凌 712100；2西北农林科技大学旱区节水农业研究院，陕西杨凌 712100；3河海大学水文水资源学院，南京 210098）

河套灌区春小麦生产水足迹影响因子敏感性及贡献率分析

孙世坤1，2，刘文艳1，刘 静3，王玉宝1，2，陈帝伊1，吴普特1，2

（1西北农林科技大学旱区农业水土工程教育部重点实验室，陕西杨凌 712100；2西北农林科技大学旱区节水农业研究院，陕西杨凌 712100；3河海大学水文水资源学院，南京 210098）

【目的】农业水资源高效利用是保障国家粮食和水资源安全的重要途径，作物用水评价是农业用水管理的主要研究课题之一。水足迹为农业用水评价提供了新的指标体系，研究旨在通过对作物生产水足迹影响因素进行定量评价，进而进行水足迹调控，以便实现农业水资源高效利用。【方法】以河套灌区为研究区域，基于水足迹概念体系，利用改进的水足迹量化方法对河套灌区春小麦生产水足迹进行量化并分析其在研究时段内的演变特征；利用单因素轮换（One-At-A-Time，OAT）敏感性分析方法和贡献率分析方法探究气候、农业生产投入因子和水资源利用效率对春小麦生产水足迹变化的驱动力。【结果】春小麦生产水足迹在研究时段内呈显著下降趋势，从1981年的4.71 m3·kg-1，下降到2010年的1.52 m3·kg-1；同时年际变化呈现出较为明显的阶段性特征，可划分为波动性下降期（1981—1987年）、快速下降期（1988—1995年）、缓慢下降期（1996—2010年），该变化规律与灌区农业生产和灌溉水平的发展特征基本一致；从水足迹的蓝、绿水构成来看，河套灌区春小麦生产水足迹中蓝水足迹比例超过90%，而绿水足迹比例不足10%，这与灌区农业生产用水特征相一致。敏感性分析显示，水足迹对各个因子的敏感性差别十分显著，日照时数、相对湿度、降水量、灌溉水利用系数和单位面积化肥用量在±20%波动的情况下，春小麦生产水足迹的波动范围分别为±30%、±24%、±2%、±63%和±4%，春小麦生产水足迹对灌溉水利用系数、日照时数和相对湿度的敏感性较高，对其余影响因子的敏感性较低。贡献率分析显示，研究时段内相对湿度的减少和降水量的增加促使春小麦生产水足迹的增加。而日照时数的下降、化肥使用量以及灌区用水效率的提高促使了春小麦生产水足迹的降低。定量分析结果显示，化肥和灌溉水利用系数对春小麦生产水足迹变化的贡献率分别为-36.89%和-39.42%，而气候因子的综合贡献率仅为 2.80%，对研究时段内春小麦生产水足迹下降贡献率最大的是灌溉水利用系数，其次为单位面积化肥用量，而相对湿度、日照时数和降水量三者的贡献率相近，贡献率最小的是日照时数，这主要是因为日照时数在研究时段内的变率较小。【结论】气候、农业生产资料投入和水资源利用效率是影响作物生产水足迹的主要因素，就河套灌区而言，农业生产资料投入和水资源利用效率的提高是促使灌区春小麦生产水足迹下降的主要因素，而气候因子在研究时段内对春小麦生产水足迹的影响较小。

水足迹；敏感性分析；贡献率分析；河套灌区；春小麦

0　引言

【研究意义】水资源是重要的自然资源和经济资源，随着人口增长与各行业用水增加，水资源供需矛盾不断加剧。作为人口和农业大国，中国面临着严峻的水资源和粮食安全压力。在气候变化的背景下，全球气温、降水和水资源分布格局都将发生变化，从而对区域水资源状况和农业生产带来不确定性［1］。近年来，随着中国经济与社会的快速发展，水资源短缺形势愈来愈严峻。农业是中国主要的水资源消耗行业，其用水占全国用水总量的 60%左右［2］，而灌溉用水占农业用水的近90%，灌区水资源利用效率的提高对于提高中国农业水资源利用效率，缓解农业用水压力具有重要意义。同时，灌溉农业对保障中国粮食安全发挥着重要作用，占全国耕地面积49%的灌区生产出约占全国总量 75%的粮食［3］。因此，合理评价灌区作物生产过程中水资源利用情况可为提高灌区农业水资源管理和调控水平提供理论依据，实现农业水资源的持续、高效利用。水足迹作为一个较新的水资源评价指标体系，从水资源消耗数量、类型和效率等几个层面评价了人类生产、消费活动对水资源的多重影响［4］，有别于传统的用水量评价指标，它是一个体现消耗的水量、水源类型以及污染量和污染类型的多层面指标。同时，作物生产水足迹的量化是进行农产品虚拟水贸易分析的基础，开展水足迹与虚拟水的研究有助于理解淡水资源的全球属性及消费和贸易对水资源使用的影响。对作物生产水足迹影响因子敏感性进行研究有助于明晰作物生产水足迹的主要影响因素，并为水足迹调控和农业水资源管理提供理论依据。【前人研究进展】目前，国际上对水足迹的研究大致集中在3个方面：一是水足迹理论框架及量化方法的研究；二是水足迹的案例研究；三是基于水足迹的水资源评价和管理研究。在水足迹理论及量化方法研究方面：HOEKSTRA等［4］探讨了水足迹在水资源管理中的应用前景，并初步构建了水足迹量化与可持续评价框架。MEKONNEN等［5］构建了基于栅格尺度的动态水分平衡模型，并在此基础上对全球作物的生产水足迹进行了量化。在水足迹案例研究方面：龙爱华等［6］、马静等［7］、HOEKSTRA等［8］、CAZCARRO等［9］分别以省（区）、国家和全球为空间分析单元，研究了区域、国家和全球的消费水足迹以及行业水足迹（旅游业）。在基于水足迹的水资源评价方面：ZHANG等［10］指出水足迹为明确人类消费活动对水资源的影响提供了新的途径，通过进行水足迹评价，可以帮助管理机构实现水资源的可持续利用。ERCIN等［11］分析了未来人口、经济发展及生产和消费结构变化情景下，全球消费和生产水足迹的变化情况，为政府和相关机构制定水资源管理策略提供参考［12］。上述研究对于水足迹理论和应用的发展提供了有益的借鉴，但目前对于农业生产水足迹的量化方法还比较单一，针对不同尺度和评价目标的量化体系尚未建立，以往对作物生产水足迹的量化主要通过计算田间尺度作物蒸散量求得，没有考虑灌溉水在输配水及田间灌水过程中的损失，因而无法反映区域尺度农业生产过程中的水资源利用状况。同时，有关水足迹的研究大多还停留在对水足迹的量化和评价方面，对作物生产水足迹影响因子敏感性及贡献率的研究还较少。ZHUO等［13］评价了降水量、作物系数、生育期等对作物生产水足迹的影响，开始了作物生产水足迹影响因素的初步探讨。【本研究切入点】农业生产是一个复杂的自然、社会交融的生产过程，气候等自然因素对作物生产水足迹产生重要的影响，但“绿色革命”以来农业生产投入对作物单产水平的提高发挥了重要作用，而作物单产是影响作物生产水足迹的重要因素，同时灌溉水利用效率水平也将直接影响作物生产水足迹。因此，综合评价气候、农业生产因子和灌溉水平对作物生产水足迹的影响对于全面评价作物生产水足迹的影响因素，进而提出水足迹调控标准和措施具有重要意义。【拟解决的关键问题】以河套灌区为研究区域，基于水足迹概念体系，利用改进的水足迹量化方法对河套灌区春小麦生产水足迹进行量化并分析其在研究时段内的演变特征；利用单因素轮换（One-At-A-Time，OAT）敏感性分析方法和贡献率分析方法探究气候、农业生产投入因子和灌溉水平对春小麦生产水足迹变化的驱动力，以期为灌区水足迹调控和水资源管理提供依据。

1　数据与方法

1.1研究区概况

河套灌区位于内蒙古巴彦淖尔市（图1），灌区地处干旱、半干旱荒漠草原地带，大陆性气候特征明显，冬季严寒少雪，夏季高温干热，年平均气温6—8℃，灌区降雨量少，年平均降水量为130—215 mm，年蒸发量为2 100—2 300 mm［14］。灌区有效灌溉面积57.4万hm2，灌区主要种植作物为春小麦、玉米和向日葵等。引黄水是灌区主要的水源，河套灌区多年平均引黄水量为50亿m3左右。但由于黄河年径流量的下降，同时随着黄河统一调度，河套灌区配水量日益减少，灌区水资源供需矛盾不断加剧［15］。

图1　河套灌区示意图Fig.1 Location of Hetao irrigation district

1.2数据及来源

研究所需数据包括：气象数据、农业统计资料和农业用水资料。研究所需的气象数据包括河套灌区1981—2010年月平均最高气温、月平均最低气温、相对湿度、风速、日照时数、降水量等参数，数据来自中国气象科学数据共享网和巴彦淖尔市气象局［16］。研究区作物产量、播种面积、单产、生产资料投入等农业生产数据来自灌区管理年报和《内蒙古统计年鉴》［17］。灌区历年引、排水、灌溉定额、灌溉水利用系数等数据来源于河套灌区管理年报、《巴彦淖尔市水资源公报》等资料［18］。

1.3作物生产水足迹计算方法

作物生产水足迹可表述为生产作物单位产量（通常为经济产量）所消耗的广义水资源量，它由蓝水足迹（作物生产所消耗的地表水或地下水，即灌溉水）、绿水足迹（作物生产过程消耗的有效降水）和灰水足迹（用于稀释作物生产过程排放的污染物，使其达到环境标准所需的水资源量）组成。其中蓝水足迹与绿水足迹为作物消耗性用水，而灰水足迹为改善环境性用水。由于河套灌区所面临的主要环境问题为土壤盐渍化，而目前国际上关于农业灰水足迹的污染源多界定为化肥、农药等造成的水体污染，因此，本研究主要关注春小麦生产过程中的消耗性用水（即蓝水足迹与绿水足迹），暂不考虑灰水足迹［4，12］。目前国际上对作物生产水足迹的量化主要通过计算田间尺度作物蒸散量求得，没有考虑灌溉水在输配水及田间灌水过程中的损失，因而无法反映区域尺度农业生产过程中的水资源利用状况。针对该方法的不足，本研究在已有研究的基础上，根据作物生产水足迹定义，提出了改进的作物生产蓝水足迹计算方法，即将灌溉水损失量中的渠道蒸发、深层渗漏后随地下水流出研究区域外的灌溉水量纳入蓝水足迹量化之中，从而反映区域作物生产过程中的水资源利用特征，主要计算过程如下［15］：

式中，WFgreen为春小麦绿水足迹（m3·kg-1）；WFblue为春小麦蓝水足迹（m3·kg-1）；Wgreen为春小麦单位面积消耗的绿水资源量（m3·hm-2）；Wblue为春小麦单位面积消耗的蓝水资源量（m3·hm-2）；10为单位换算系数；ETc为春小麦生育期蒸发蒸腾量（mm）；Pe为春小麦生育期有效降水量（mm）；IRG为作物单位面积消耗的灌溉水量（m3·hm-2）；Y为作物单位面积产量（kg·hm-2）。

Pe根据 CROPWAT模型中内嵌的有效降水量计算公式进行量化，该公式为美国农业部土壤保持局（USDA SCS）提出的方法［19］：

式中，Pdec为旬降水量（mm）；Pe（dec）为旬有效降水量（mm）。根据作物生育期对Pe（dec）进行求和即可得到作物整个生育期的有效降水量，沈成［20］、AUGUSTIN［21］等人的研究表明，该方法在河套灌区具有较好的适用性。

式中，WA为灌区种植业灌溉水总消耗量（m3）；α为春小麦灌溉水消耗量占灌区灌溉水消耗总量的比例；AW为春小麦的播种面积（hm2）。

WA根据灌区水量平衡进行计算，计算过程见下式：

式中，WIn为灌区引水量（m3）；WP为灌区降水补给量（m3）；WG为地下水补给量（m3）；Wout为灌区排水量（m3）；ΔW为灌区内部蓄水的变化量（m3）；WPe为灌区种植业消耗的有效降水量（m3）；WI为工业生产耗水量（m3）；WL为城乡生活耗水量（m3）；WE为生态环境耗水量（m3）。

α根据下式计算［15］：

1.4水足迹敏感性分析

因子敏感性是表征一个或几个相关影响因子发生变化时，对因变量影响程度的定量参数［22-24］。论文利用OAT（One-at-a-Time）敏感性分析法分析春小麦生产水足迹对其影响因素的敏感性。OAT分析是一种筛选式的试验方法，综合了局部和全局的敏感性分析方法，通过改变准则因子的相关数据来评估因变量对自变量的敏感性［25-27］。

1.5因子贡献率分析

根据水足迹基本理论，作物生产水足迹是由作物耗水量与单位面积产量共同决定，作物单产受当地热量、水分、辐射等农业气候资源，以及化肥、农药、机械化水平等生产资料因子的影响，用函数表示其关系为［15］：

式中，Ci、Si分别为影响作物单位面积产量的气象因子及生产资料投入因子。

作物生产过程的水资源消耗主要与区域气候条件及灌溉设施和技术水平有关，作物耗水量可用下式表示［15］：

式中，CWU为作物单位面积水资源消耗量（m3·hm-2）；Si为作物生产过程中与水资源消耗相关的生产变量，如灌溉水利用系数。

因此，作物生产水足迹可用下式表示［15］：

对等式两边分别取对数进行线性化可得：

对上式进行多元线性回归得到以下回归方程：

式中，e为回归方程常数项，其他参数与上同。

根据弹性系数的定义，利用Eviews 8.0软件对上式进行多元线性回归，回归系数即为各因子的弹性系数（一定时期内相互联系的两个指标变化速度的比率），农业生产资料和气象因子变化对水足迹变化的贡献率可按下式计算［15，28-29］：

式中，δi为第 i项生产资料对作物生产水足迹变化的贡献率；αi为第i项生产资料的弹性系数；ΔSi为第i项生产资料投入的变化量；ΔWF为作物生产水足迹的变化量；γi为第i项气象因子对作物生产水足迹变化的贡献率；bi为第i项气象因子的弹性系数；ΔCi为第i项气象因子的变化量。

2　结果

2.1春小麦生产水足迹年际变化分析

春小麦生产过程的水资源消耗有两种类型：蓝水（灌溉水）和绿水（有效降水）。对河套灌区春小麦生产过程中的水资源消耗分析显示，在研究时段内，春小麦生育期水资源消耗主要依靠灌溉水供给，而由于降水稀少，春小麦消耗的降水较少。从年代变化来看，由于区域灌溉系统运行水平的提高，单位面积灌溉水消耗量呈波动下降的趋势，而降水消耗量受降水波动性的影响，其在研究时段内也呈现出较为明显的波动趋势（图2）。

根据灌区气象、农业生产和灌区用水数据对河套灌区1981—2010年的春小麦生产水足迹进行量化。利用趋势分析和M-K检验［30-31］分析了春小麦生产水足迹的时间演变特征。图3中生产水足迹在研究时段内呈显著下降趋势，从1981年的4.71 m3·kg-1，下降到2010年的1.52 m3·kg-1，年均下降率为0.11 m3·kg-1。M-K趋势检验值为-6.71，说明下降趋势达到极显著水平（P＜0.01）。同时，春小麦生产水足迹年际变化呈现出较为明显的阶段性特征：1981—1987年间主要呈现出波动性下降的趋势；1988—1995年间呈现出快速下降的变化趋势；1995年之后呈现出较为平缓的下降趋势。从水足迹的蓝、绿水构成来看，河套灌区春小麦生产水足迹中蓝水足迹比例超过90%，而绿水足迹比例不足10%，灌区春小麦生产主要依靠蓝水资源。从二者在研究时段内的变化趋势来看，春小麦生产蓝水足迹下降趋势较春小麦生产绿水足迹更为显著。

图2　区域尺度春小麦单位面积耗水量年际变化Fig.2 Interannual variation of water consumption of spring wheat at regional scale

图3　春小麦生产水足迹年际变化Fig.3 Interannual variation of water footprint of spring wheat

春小麦生产水足迹是由春小麦生产过程中的水资源消耗量和春小麦单位面积产量共同决定，通过分析研究时段二者的变化规律可以发现（图 4），春小麦生产水足迹的快速下降阶段与春小麦单产大幅度提高（伴随着化肥、农药使用量的增加）的阶段基本吻合，进入1995年之后，春小麦单产增速放缓，其水足迹的下降主要是在单位面积耗水量下降的带动下实现的。

2.2影响因子敏感性分析

鉴于河套灌区农业生产主要依靠蓝水资源，春小麦蓝水足迹占春小麦生产水足迹总量的比例在90%以上，春小麦蓝水足迹与春小麦总水足迹的影响因子敏感性具有大体一致的特性。因此，本研究着重分析春小麦总水足迹的因子敏感性及贡献率。利用 OAT方法对春小麦生产水足迹的影响因素进行敏感性分析，结果显示（图 5）：对于气象因子而言，春小麦生产水足迹与日照时数和降水量呈正相关关系，与相对湿度呈现负相关关系，即春小麦生产水足迹随着日照时数和降水量的增加而增加，随着相对湿度的增加而减小。对于农业生产因子而言，春小麦生产水足迹与化肥使用量和灌溉水利用效率呈负相关关系，即春小麦生产水足迹随着化肥使用量的增加、灌溉水利用系数的提高而降低。

图4　春小麦单位面积产量及耗水量年际变化Fig.4 Interannual variation of yield and water consumption per unit area of spring wheat

气象因子中，水足迹对日照时数的敏感性最大，日照时数在±20%的波动的情况下，春小麦生产水足迹的变化范围约为±30%；相对湿度的敏感性次之，在±20%的波动的情况下，春小麦生产水足迹的变化范围约为±24%；而水足迹对降水的敏感性较低。农业生产水平因子中，水足迹对灌溉水利用率的敏感性较高，灌溉水利用系数在±20%变动的幅度下，春小麦生产水足迹的变动幅度达到±60%以上；但是春小麦生产水足迹对化肥使用量的敏感性较低，单位面积化肥用量在±20%的波动的情况下，春小麦生产水足迹的变化范围约为±4%左右。由此可见，水足迹对各个因子的敏感性差别十分显著，春小麦生产水足迹对灌溉水利用系数、日照时数和相对湿度的敏感性较高，对其余影响因子的敏感性较低。

2.3影响因子贡献率分析

作物生产水足迹影响因子在研究时段内的变化对春小麦生产水足迹产生一定的影响，为了评价影响因子变化对春小麦生产水足迹变化的贡献程度，采用贡献率分析方法定量分析了各个影响因子对春小麦生产水足迹变化的贡献程度。影响因子变化率分析结果显示：相对湿度和日照时数在研究时段内的变率分别为-3.04%和-0.35%，而降水量、化肥使用量和灌溉水利用系数在研究时段内的变率分别为31.40%、406.70% 和26.22%，这说明灌区相对湿度和日照时数在研究时段内呈现出下降趋势，而降水量、化肥使用量和灌溉水利用系数在研究时段内呈增加态势。

上述影响因子在研究时段内的变化将会对春小麦生产水足迹产生影响，各影响因子贡献率分析结果表明（表1），相对湿度和降水量的贡献率为正值，日照时数、单位面积化肥用量及灌溉水利用系数的贡献率为负值。这说明，研究时段内相对湿度的减少和降水量的增加促使春小麦生产水足迹的增加。而日照时数的下降、化肥使用量以及灌区用水效率的提高促使了春小麦生产水足迹的降低。定量分析结果表明，对春小麦生产水足迹值贡献率最大的是灌溉水利用系数，其次为单位面积化肥用量，而相对湿度、日照时数和降水量三者的贡献率相近，贡献率最小的是日照时数，这主要是因为日照时数在研究时段内的变率最小。

图5　春小麦生产水足迹影响因子敏感性Fig.5 Sensitivity of influencing factors of wheat water footprint

进一步分析因子贡献率差异的原因可以发现，灌溉水利用系数的弹性系数和研究时段内的变率均较大，所以其贡献率也最大。相对湿度和日照时数的弹性系数虽然较大，但是其研究时段内的变率较小，所以对水足迹的贡献也较小。此外，虽然相对湿度的敏感系数为负值，但是由于年际减少的变化趋势，其对春小麦生产水足迹的贡献率为正。

表1　春小麦生产水足迹影响因子贡献率分析Table 1 Contribution rate analysis of influencing factors of wheat water footprint

3　讨论

气象因子对作物生产水足迹的影响主要通过其对作物单位面积产量和耗水量的作用实现。关于气象因子对作物耗水量的影响，已有学者进行了系统的研究。陈素英等［32］研究指出相对湿度的变化将会影响作物蒸散和气孔导度，湿度越大将会导致大气蒸发能力下降，作物需水量下降。同时湿度增加将会导致作物气孔导度减小，影响作物光合作用，降低作物生物量［32］。日照时数影响着作物生育期的太阳辐射量，从而对作物生物量的积累产生影响，另一方面，根据水气扩散理论，太阳辐射将为作物蒸散提供能量，因此日照时数的改变将会影响太阳辐射从而影响作物产量和作物耗水量。赵璐等［27］对川中丘陵区近 60年参考作物蒸发蒸腾量变化成因进行了研究，结果显示日照时数和风速减小对参考作物蒸发蒸腾量变化的贡献率分别为-5.403%和-2.054%；相对湿度减小和温度增加对参考作物蒸发蒸腾量变化的贡献率分别为 2.085%和 0.562%，研究指出日照时数和风速减小是引起该区域参考作物蒸发蒸腾量下降的主要原因。

气候因子对作物生产水足迹的影响途径如下：随着空气相对湿度的增加大气蒸发能力下降，作物耗水量也随之下降，从而促使作物生产水足迹下降。日照时数决定着作物生育期的太阳辐射量，日照时数的改变将会影响太阳辐射从而影响作物产量和作物耗水量，在二者共同作用下对作物生产水足迹产生影响。降水作为绿水的来源，只有转化为土壤有效水并被作物利用的那部分降水构成了作物绿水足迹，由于河套灌区地处西北内陆干旱区，夏季高温干旱、冬季严寒少雪，降水稀少，农业生产主要依靠灌溉，因此该区域作物生产水足迹受降水的影响较小。农业生产资料投入及灌溉水平对作物生产水足迹影响途径如下：施肥能提高土壤肥力，是提高作物单位面积产量的有效措施，因此化肥使用量的增加将会促使作物生产水足迹的下降。而灌溉水利用系数则直接决定了蓝水资源的利用效率，灌溉水利用系数的提高将会使作物生育期的水资源消耗量下降，从而促使作物生产水足迹的下降。

河套灌区春小麦生产水足迹、单位面积产量和耗水量的年际变化结果显示，21世纪90年代之后，由于作物单位面积产量年均增加幅度放缓，因而通过提高作物单产水平以降低作物生产水足迹的空间有限，同时LIU等［33］和ZHANG等［34］的研究表明由于化肥和农药的“边际效益”递减，二者使用所带来的单产水平提高幅度越来越小。这一点也从研究中水足迹对化肥的敏感性较低这一结论得到佐证，而且使用过多的化肥、农药将会带来农田生态环境的恶化。因此，河套灌区作物生产水足迹调控的重点和发展方向是提高农业水资源利用率和利用效率，降低作物生产的水资源消耗量。春小麦生产水足迹影响因子敏感性分析结果表明，对于河套灌区而言，春小麦生产水足迹对灌溉水利用效率敏感性较高，因此通过农田水利设施和雨水积蓄利用设施建设提高灌溉水和降水利用系数，同时合理控制化肥和农药用量，在维持高产的同时，减小化肥、农药对区域土壤、水环境的影响，是实现区域农业水资源高效、可持续利用的主要措施。

本研究利用敏感性和贡献率分析对水足迹影响因素的敏感性和贡献率进行了研究，并得到了一些具有参考价值的初步结论，该结论可为灌区实施水足迹调控，实现水资源定量化管理提供有效的途径和方向。但本研究还存在以下局限性：数据是水足迹评价研究的基础，水足迹量化及评价所涉及的数据面广、量大，对数据准确性、完整性和有效性都有较高的要求，由于数据统计口径和误差，导致基础数据本身存在一定的不确定性，同时受限于基础数据的时间序列较短，可能导致敏感性和贡献率分析结果存在一定的误差。另外，本研究利用传统的水量平衡进行水足迹计算，未来可借助开发适用于灌区的分布式水文模型对水足迹进行量化，可为分析农业水文循环过程和机理提供借鉴。

4　结论

（1）在研究时段内，河套灌区春小麦生育期水资源消耗主要依靠灌溉水供给。从年代变化来看，由于区域灌溉系统运行水平的提高，单位面积灌溉水消耗量呈波动下降的趋势，而降水消耗量受降水波动性的影响，其在研究时段内也呈现出较为明显的波动趋势。春小麦生产水足迹在研究时段内呈显著下降趋势，年均下降率为0.11 m3·kg-1，同时其年际变化呈现出较为明显的阶段性特征，可划分为波动性下降期（1981—1987年）、快速下降期（1988—1995年）、缓慢下降期（1996—2010年），这个变化规律与灌区农业生产和灌溉水平的发展特征基本一致。

（2）河套灌区春小麦生产水足迹与日照时数和降水量呈现正相关关系，与相对湿度、化肥使用量和灌溉水利用系数呈负相关关系；敏感性定量分析显示，春小麦生产水足迹对灌溉水利用系数最敏感，其次是日照时数、相对湿度、化肥使用量和降水量。贡献率分析结果显示农业生产资料（表现为化肥的投入）和水资源利用效率（表现为灌溉水利用系数）的提高是促使河套灌区春小麦生产水足迹在研究时段内下降的主要因素，而气候因素由于在研究时段内变率较小，其对春小麦生产水足迹的影响较小。

（3）作物生产水足迹对化肥的敏感性较低表明：由于“边际效益”递减，化肥和农药的使用所带来的单产水平提高幅度越来越小，而且过多的使用化肥、农药将会带来农田生态环境的恶化。根据因子敏感性和贡献率分析结果，河套灌区作物生产水足迹调控的方向是提高农业水资源利用率和利用效率，降低作物生产的水资源消耗量。

References

［1］ IPCC. Climate Change 2014: Impacts，Adaptation，and Vulnerability Contribution of Working Group II to the Fifth Assessment Report. Cambridge: Cambridge University Press，2014.

［2］ 中华人民共和国水利部. 2010中国水资源公报. 北京: 中国水利水电出版社，2011. Ministry of Water Resources of China. China Water Resources Bulletin 2010. Beijing: China Water & Power Press，2011. （in Chinese）

［3］ 齐学斌，黄仲冬，乔冬梅，张现超，李平，Andersen M N. 灌区水资源合理配置研究进展. 水科学进展，2015，26（2）: 1-9. QI X B，HUANG Z D，QIAO D M，ZHANG X C，LI P，ANDERSEN M N. Research advances on the reasonable water resources allocation in irrigation district. Advances in Water Science，2015，26（2）: 1-9. （in Chinese）

［4］ HOEKSTRA A Y，CHAPAGAIN A K，ALDAYA M M，MEKONNEN M M. The Water Footprint Assessment Manual: Setting the Global Standard. London: Earthscan，2011.

［5］ MEKONNEN M M，HOEKSTRA A Y. The green，blue and grey water footprint of crops and derived crop products. Hydrology and Earth System Sciences，2011，15（5）: 1577-1600.

［6］ 龙爱华，张志强，徐中民，苏志勇. 甘肃省水资源足迹与消费模式分析. 水科学进展，2005，16（3）: 418-425. LONG A H，ZHANG Z Q，XU Z M，SU Z Y. Analysis of Chinese cultivated land’s spatial-temporal changes and causes in recent 30 years. Advances in Water Science，2005，16（3）: 418-425. （in Chinese）

［7］ 马静，汪党献，来海亮，王茵. 中国区域水足迹的估算. 资源科学，2005，27（5）: 96-100. MA J，WANG D X，LAI H L，WANG Y. Water footprint-an application in water resources research. Resources Science，2005，27（5）: 96-100. （in Chinese）

［8］ HOEKSTRA A Y，MEKONNEN M M. The water footprint of humanity. PNAS，2012，109（9）: 3232-3237.

［9］ CAZCARRO I，HOEKSTRA A Y，SÁNCHEZ CHÓLIZ J. The water footprint of tourism in Spain. Tourism Management，2014，40（6）: 90-101.

［10］ ZHANG G P，HOEKSTRA A Y，MATHEWS R E. Water footprint assessment （WFA） for better water governance and sustainabledevelopment. Water Resources and Industry，2013（1/2）: 1-6.

［11］ ERCIN A E，HOEKSTRA A Y. Water footprint scenarios for 2050: A global analysis. Environment International，2014，64（1）: 71-82.

［12］ 孙世坤，王玉宝，吴普特，赵西宁. 小麦生产水足迹区域差异及归因分析. 农业工程学报，2015，31（13）: 142-148. SUN S K，WANG Y B，WU P T，ZHAO X N. Spatial variability and attribution analysis of water footprint of wheat in China. Transactions of the Chinese Society of Agricultural Engineering，2015，31（13）: 142-148. （in Chinese）

［13］ ZHUO L，MEKONNEN M M，HOEKSTRA A Y. Sensitivity and uncertainty in crop water footprint accounting: A case study for the Yellow River basin. Hydrology and Earth System Sciences，2014，18（6）: 2219-2234.

［14］ 白岗栓，张蕊，耿桂俊，任志宏，张沛琪，史吉刚. 河套灌区农业节水技术集成研究. 水土保持通报，2010，31（1）: 149-154. BAI G S，ZHANG R，GENG G J，REN Z H，ZHANG P Q，SHI Z G. Integrating agricultural water-saving technologies in Hetao irrigation district. Bulletin of Soil and Water Conservation，2010，31（1）: 149-154. （in Chinese）

［15］ SUN S K，WU P T，WANG Y B，ZHAO X N，LIU J，ZHANG X H. The impacts of interannual climate variability and agricultural inputs on water footprint of crop production in an irrigation district of China. Science of the Total Environment，2013，444（8）: 498-507.

［16］ 国家气象局信息中心. 中国气象科学数据共享服务网［DB/OL］. ［2014-06-02］. http://cdc.cma.gov.cn/home.do，2010. National Meteorological Information Center. China Meteorological Data Sharing Service System［DB/OL］. ［2014-06-02］. http://cdc.cma. gov.cn/home.do，2010. （in Chinese）

［17］ 内蒙古自治区统计局. 内蒙古统计年鉴 1981-2010. 北京: 中国统计出版社，1982-2011. Inner Mongolia Statistical Bureau. Inner Mongolia Statistical Yearbook 1981-2010. Beijing: China Statistics Press，1982-2011. （in Chinese）

［18］ 巴彦淖尔市水务局. 巴彦淖尔市水资源公报［DB/OL］. ［2014-03-02］. http://www.htgq.gov.cn/news/show.asp？id=4379. Water Authority in Bayannur. Water Resources Bulletin in Bayannur ［DB/OL］. ［2014-03-02］. http://www.htgq.gov.cn/news/show.asp？id =4379. （in Chinese）

［19］ DOLL P，SIEBERT S. Global modeling of irrigation water requirements. Water Resources Research，2002，38（4）: 1037-1048.

［20］ 沈成. 基于 CROPWAT模型河套灌区套种模式非充分灌溉制度模拟与区域决策［D］. 内蒙古: 内蒙古农业大学，2012. SHEN C. Insufficient irrigation system simulation and regional strategy under interplanting mode based on CROPWAT model in Hetao irrigation area［D］. Inner Mongolia: Inner Mongolia Agricultural University，2012. （in Chinese）

［21］ AUGUSTIN L K. 基于GIS和CROPWAT的内蒙古河套灌区作物的最优灌溉制度研究［D］.北京: 中国农业大学，2015. AUGUSTIN LK. Optimal irrigation scheduling for crops: Based on GIS and CROPWAT application in Hetao district of Inner Mongolia Autonomous Region，China［D］. Beijing: China Agricultural University，2015. （in Chinese）

［22］ 黄清华，张万昌. SWAT模型参数敏感性分析及应用. 干旱区地理，2010，33（1）: 8-15. HUANG Q H，ZHANG W C. Application and parameters sensitivity and analysis of SWAT model. Arid Land Geography，2010，33（1）: 8-15. （in Chinese）

［23］ 田雨，雷晓辉，蒋云钟，李薇. 水文模型参数敏感性分析方法研究评述. 水文，2010，30（4）: 9-12. TIAN Y，LEI X H，JIANG Y Z，LI W. Comment on parameter sensitivity analysis of hydrological model. Journal of China Hydrology，2010，30（4）: 9-12. （in Chinese）

［24］ HAMBY D M. A review of techniques for parameter sensitivity analysis of environmental models. Environmental Monitoring and Assessment，1994，32（2）: 135-154.

［25］ SUN X Y，NEWHAM L. Three complementary methods for sensitivity analysis of a water quality model. Environmental Modelling & Software，2012，37（1）: 19-29.

［26］ 曹雯，申双和，段春锋. 西北地区生长季参考作物蒸散变化成因的定量分析. 地理学报，2011，66（3）: 407-415. CAO W，SHEN S H，DUAN C F. Quantification of the causes for reference crop evapotranspiration changes in growing season in Northwest China. Acta Geographica Sinica，2011，66（3）: 407-415. （in Chinese）

［27］ 赵璐，梁川，崔宁博，刘利文. 川中丘陵区参考作物蒸发蒸腾量近60年变化成因研究. 水利学报，2013，44（2）: 183-190. ZHAO L，LIANG C，CUI N B，LIU L W. Attribution analyses of ET0change in hilly area of central Sichuan in recent 60 years. Journal of Hydraulic Engineering，2013，44（2）: 183-190. （in Chinese）

［28］ 刘忠，黄峰，李保国. 2003-2011年中国粮食增产的贡献因素分析.农业工程学报，2013，29（23）:1-8. LIU Z，HUANG F，LI B G. Investigating contribution factors to China’s grain output increase in period of 2003 to 2011. Transactions of the Chinese Society of Agricultural Engineering，2013，29（23）: 1-8 （in Chinese）

［29］ 王丹. 气候变化对中国粮食安全的影响与对策研究［D］. 武汉: 华中农业大学，2009. WANG D. Impact of climate change on Chinese grain security and countermeasures［D］. Wuhan: Huazhong Agricultural University，2009. （in Chinese）

［30］ HAMED K H. Trend detection in hydrologic data: The Mann-Kendall trend test under the scaling hypothesis. Journal of Hydrology，2008，349（3/4）: 350-363.

［31］ 张淑杰，张玉书，隋东，蔡福，武晋雯，纪瑞鹏，陈鹏狮，刘庆婺.东北地区参考蒸散量的变化特征及其成因分析. 自然资源学报，2010，25（10）:1750-1761. ZHANG S J，ZHANG Y S，SUI D，CAI F，WU J W，JI R P，CHEN P S，LIU Q W. Changes in reference evapotranspiration and its causes in Northeast China. Journal of Natural Resources，2010，25（10）: 1750-1761. （in Chinese）

［32］ 陈素英，张喜英，邵立威，孙宏勇. 农业技术和气候变化对农作物产量和蒸散量的影响. 中国生态农业学报，2011，19（5）: 1039-1047. CHEN S Y，ZHANG X Y，SHAO L W，SUN H Y. Effects of climate change and agricultural technology improvement on evapotranspiration and crop yield. Chinese Journal of Eco-Agriculture，2011，19（5）: 1039-1047. （in Chinese）

［33］ LIU X J，ZHANG Y，HAN W X，TANG A H，SHEN J L，CUI Z L，VITOUSEK P，ERISMAN J W，GOULDDING K，CHRISTIE P，FANGMEIER A，ZHANG F S. Enhanced nitrogen deposition over China. Nature，2013，494（7306）: 459-462.

［34］ ZHANG F S，CHEN X P，VITOUSEK P. An experiment for the world. Nature，2013，497（7447）: 33-35.

（责任编辑 杨鑫浩）

Sensitivity and Contribution Rate Analysis of the Influencing Factors of Spring Wheat Water Footprint in Hetao Irrigation District

SUN Shi-kun1，2，LIU Wen-yan1，LIU jing3，WANG Yu-bao1，2，CHEN Di-yi1，WU Pu-te1，2
（1Key Laboratory of Agricultural Soil and Water Engineering in Arid and Semiarid Areas of Ministry of Education，Northwest A&F University，Yangling 712100，Shaanxi；2Institute of Water Saving Agriculture in Arid regions of China，Northwest A&F University，Yangling 712100，Shaanxi；3College of Hydrology and Water Resources，Hohai University，Nanjing 210098）

【Objective】The efficient utilization of agricultural water resources is a key measure to guarantee national food and water security. Crop water use efficiency evaluation is one of the main research fields of agricultural water management. Water footprint provides a new index for agricultural water use evaluation，and the quantitative evaluation on the influencing factors of crop water footprint will be helpful to the implementation of water footprint control and improvement of the agricultural water use efficiency. 【Method】Based on the concept of water footprint，the water footprint of spring wheat in Hetao irrigation district was quantified by using an improved calculation method and the temporal variation was analyzed. Sensitivity and contribution rate analysis were used to quantify the relationship between crop water footprint and its influencing factors. 【Result】The results show that water footprint of wheat declined significantly during the study period. It decreased from 4.71 m3·kg-1in 1981 to 1.52 m3·kg-1in 2010. The variation of water footprint of wheat displayed an obvious stage characteristic. It can be divided into three stages: fluctuate declining stage （1981-1987），rapid declining stage （1988-1995） and slow declining stage （1996-2010）. And this variation characteristic was consistent with the variation of agricultural production and irrigation level in Hetao irrigation district. The blue water footprint accounted for the larger proportion （more than 90%），while for the share of green water footprint it was relatively small. Therefore，the production of wheat in the Hetao irrigation district mainly depended on blue water （irrigation water）. Sensitivity analysis shows that the difference of sensitivity between the influencing factors was significant. The variations of water footprint of spring wheat was ±30%，±24%，±2%，±63% and ±4% when sunshine hours，relative humidity，precipitation，irrigation water use coefficient and fertilizer rate per unit area varied at ±20%. Irrigation water use coefficient is the most sensitivity factor of wheat water footprint，following by sunshine hours，relative humidity，fertilizer usage and precipitation. Contribution analysis results show that the decline of relative humidity and the increase of precipitation led to the increase of wheat water footprint. On the contrary，the decrease of sunshine hours combined with the increase of fertilizer usage and irrigation water use efficiency led to the decrease of wheat water footprint. The contribution rates of fertilizer and irrigation water use coefficient were -36.89% and -39.42%，respectively，while the contribution rate of the total climatic factors was 2.80%. The increase of utilization coefficient of irrigation water had the largest contribution rate to the decrease of wheat water footprint during the study period，followed by fertilizer usage. The relative humidity，sunshine hours and precipitation had similar contribution rate to the variation of wheat water footprint. 【Conclusion】 The main kinds of influencing factors of crop water footprint are climate，agricultural production inputs and water use efficiency. As for Hetao irrigation district，the improvement of agricultural production and water use efficiency are the major driving forces that cause the variation of wheat water footprint in Hetao irrigation district，while the climate factors have little effect on wheat water footprint. The results of this study could provide reference for the water footprint control.

water footprint； sensitivity analysis； contribution rate analysis； Hetao irrigation district； spring wheat

2015-12-14；接受日期：2016-05-27

国家自然科学基金（51409218）、中央高校基本科研业务费专项资金（2014YB050）、陕西省自然科学基础研究计划（2016JQ5092）

联系方式：孙世坤，E-mail：ssksci@163.com。通信作者王玉宝，E-mail：wyb0406@sina.com