基于水足迹的作物生产生态效率评价——以陕西省为例

惠 婞,陈晓楠,宋健峰

西北农林科技大学经济管理学院, 杨凌 712100

中国是一个农业大国,并非农业强国,在农产品刚性需求增长的情况下,农业生产面临着土地、气候和水等自然资源的制约,同时农业生产中化肥、农药的使用对环境造成损害。然而传统的农业生产函数以新古典经济学C-D生产函数或其变形为主要形式,解释变量包括土地、劳动力、农业机械总动力等[1],往往不考虑自然资源因素,特别是在陕西省相比土地资源更为稀缺的水资源,同时也不考虑过度施用化肥、农药造成的环境影响。

随着水资源的日益短缺和环境污染的加剧,在自然资源研究领域,对用水大户—农业生产中水资源消耗和农业点、面源污染的研究日益深入,水足迹概念和理论的提出将这一研究推向了新的高度。在虚拟水理论的基础上,借鉴生态足迹的概念,Hoekstra和Hung[2]以及Hoekstra和Chapagain[3]提出并不断完善了水足迹理论。水足迹(Water Footprint,简称WF)指生产产品所需要的水资源总量,该部分水资源如果来源于暂时留在未饱和土壤中,未形成径流的那部分降雨资源,称为绿水；如果来源于江、河、湖泊(灌溉用水),称为蓝水；如果生产造成污染,理论上使用淡水资源稀释污染物使其达到水质标准需要的水称为灰水[4]。作物生产水足迹概念综合考虑了作物生产中的主要资源投入和环境影响,表达了保障可持续发展的作物生产用水总量。

同时,结合资源和环境投入的生态效率作为衡量可持续发展的重要工具,受到广泛关注并取得丰富的成果。生态效率的概念最早由Schaltegger和Sturm[5]提出,后经世界可持续发展工商理事会(World Business Council for Sustainable Development,WBSCD)推广。经济合作与发展组织(OECD)将生态效率定义为“利用生态资源满足人类需求的效率”,使用“产品或服务价值与环境影响的比值”表达。该定义被普遍接受,用以说明和评估生态效率。

目前对于生态效率的评价主要有两种方法。 一种是OECD提出的比值法,另一种是基于边界的效率评价方法,数据包络分析(Data Envelopment Analysis,DEA)。初期的生态效率评价主要针对具体的产品或生产过程,比值法被广泛使用。比值法的关键在于衡量经济价值和环境影响,生产的环境影响往往有很多方面,生命周期评价(Life Cycle Assessment,LCA)通常被用来汇总这些环境影响。DEA是一种非参数的前沿方法,适用于多投入多产出的生产率评价[6]。Kuosmanen和Kortelainen[7]首次在生态效率评价中使用DEA方法。当涉及多个经济产出和环境影响时,相比单一指标比值法,DEA可以给出更综合的评价结果,目前已成为评价生态效率的主流方法[8]。

具体到有关作物生产生态效率的实证研究,Reith和Guidry[9]使用比值法评价美国南部一个可持续农业综合体(MSAC)的生态效率。Thanawong等[10]使用比值法比较不同灌溉条件下泰国东北部Lam Sieo Yai流域稻米生产的生态效率。Mülle等[11]通过评价生态效率研究新西兰Plenty湾猕猴桃的可持续生产。Gómez-Limón等[12]以及Picazo-Tadeo等[13]提出使用定向距离函数和DEA技术评估生态效率,并应用于西班牙橄榄种植农场。Hoang和Alauddin[14]使用DEA方法分析经合组织国家的农业生态效率。Beltrán-Esteve等[15]综合LCA和DEA方法评价西班牙柑橘农场的生态效率。Rosano等[16]使用随机边界的双曲距离函数研究亚马逊农业生态效率。Han等[17]使用元边界分析混合模型分析中国省际农业生态效率。Coluccia等[18]使用DEA方法评价意大利农业生态效率。

国内学者吴小庆等[19]、朱玉林等[20]和程翠云等[21]使用比值法分别评价盆栽水稻、湖南农业生态系统和我国2003—2010年的农业生态效率。许朗等[22]运用DEA方法测算出中国13个粮食主产区2000—2012年的农业生态效率,并运用Malmquist指数方法进行动态分析。洪开荣等[23]利用网络DEA模型对我国2005—2013年30个省市农业生态系统效率进行测度。在DEA模型基础上,一些研究将经济增长或农业总产值作为期望产出,将环境负效应作为非期望产出,采用SBM模型对农业生态效率进行评价[24- 26]。除了测算农业生态效率外,很多学者使用两阶段方法进一步分析影响农业生态效率的因素。王宝义和张卫国[27]采用DEA-Tobit等模型方法,研究了中国农业生态效率的省际差异和影响因素。赵哲等[28]使用超效率DEA方法对呼伦贝尔地区草牧业进行生态效率评价和影响因素分析。

与DEA并列的另一种基于边界的效率评价方法,随机前沿方法(Stochastic Frontier Approach,SFA)在各领域中被广泛应用,但却很少被用于生态效率评价。可能的原因在于SFA需要对函数形式进行正确的设定,同时其不适用于多产出模型[6]。Song等[29]将水足迹引入生产函数,可以克服了这一缺点,使SFA方法的使用成为可能。郑德凤等[30]在使用DEA方法评价甘肃省农业生态效率时进行了类似的处理,但其将作物生产灰水足迹作为非期望产出,同时未区分蓝水和绿水。与DEA方法一样,SFA方法相比比值法能提供综合的效率评价结果,同时避免比值法中因需要汇总环境影响存在的客观赋权问题。SFA方法作为一种参数分析方法,其确定的生产函数形式,相比DEA方法,可以排除随机误差的影响,效率估计和影响因素分析可以一步完成,还可以进一步分析各投入要素的效率和产出弹性,提供更深入的分析。本文将该方法应用于我国西北半干旱区的陕西省的作物生产,使用陕西省各区域1985—2018年面板数据,将作物生产的绿水、蓝水和灰水足迹(分别代表作物生产使用的降水资源、灌溉水资源和环境影响)引入生产函数,使用SFA方法评价作物生产生态效率(考虑资源投入和环境影响的作物生产技术效率)。该研究除了可为WF-SFA方法在农业生态效率评价上的应用提供新的证据外,对于作物生产水足迹,尤其是综合考虑化肥和农药对灰水足迹进行了更准确的测算,同时,使用随机前沿生态函数的参数估计进一步研究了各水足迹投入的生产效率,有效区分对降水资源的利用效率,灌溉水资源投入的效率和作物生产的环境影响,为生态效率改进提供针对性方案。

1 研究方法和数据来源

1.1 研究方法

1.1.1作物生产水足迹测算

作物生产总水足迹等于绿水、蓝水和灰水足迹之和。

WF=GWF+BWF+EWF

(1)

式中,WF为作物生产总水足迹,m3；GWF、BWF、EWF分别为作物生产绿水足迹、蓝水足迹和灰水足迹,m3。

作物生产绿水足迹和蓝水足迹由不同农作物单位面积消耗的绿水资源量和蓝水资源量按照作物种植面积加权加总得到。

(2)

式中,GWi为i作物单位面积消耗的绿水资源量,m3/hm2；BWi为i作物单位面积消耗的蓝水资源量,m3/hm2；Ai为i作物种植面积。本研究计算陕西省各区域各农田作物(小麦、稻谷、玉米、薯类、大豆、棉花、油菜、花生、向日葵、麻类、糖料、烤烟、其他小杂粮、蔬菜和瓜类)和经济林果(茶、桑和果)生产水足迹。

参考吴普特等[31]和孙世坤等[32]的计算方法,作物单位面积消耗的绿水和蓝水资源量按下式测算：

(3)

式中,ETc为作物蒸发蒸腾量,mm；Pe为作物生育期有效降水量,mm；10为单位转化系数,将单位由水深转化为单位面积水量。

ETc计算公式如下：

ETc=Kc×ET0

(4)

式中,Kc为作物系数；ET0为参考作物蒸发蒸腾量,mm。ETc使用FAO 推荐的CropWat 8.0 软件计算得到,该软件根据当地气象资料,利用Penman-Monteith公式计算参考作物蒸发蒸腾量。

作物生育期有效降水量由生育期各旬有效降水量累加得到,旬有效降水量Pe(dec)按照美国农业部土壤保持局推荐方法计算：

(5)

式中,P为作物生育期旬降水量,mm。

农业生产中施用的化肥、农药等除被植物利用的部分外,还有很大一部分通过各种形式散失到土壤和水环境中,造成污染。稀释作物生产各污染物所需的最大淡水量为作物生产灰水足迹[33],即：

EWF=max{EWi}

(6)

式中：EWi代表第i种污染物灰水足迹,m3。

各污染物灰水足迹采用Hoekstra等[34]的方法计算：

(7)

式中,Pi为第i种污染物的使用量,包括氮肥、农药等,kg；αi为第i种污染物淋溶率,%；Cmaxi为第i种污染物环境最大容许浓度,kg/m3；Cnati为第i种污染物的自然本底浓度,kg/m3。

农业生产施用的化肥中主要含有氮、 磷、钾,其中比例最大的是氮肥。考虑灰水足迹测算符合“短板原理”[33],本研究选择对灰水足迹影响较大的2个污染物,氮肥和石油类农药,测算作物生产灰水足迹。参考Mekonnen和Hoekstra[4]的研究,氮肥的溶淋率选择10%,采用US-EPA的标准,氮的环境最大容许浓度为10mg/L,假设自然水体当中氮的浓度为0。参考曹连海等[35]的研究,石油类农药的溶淋率为75%,环境最大容许浓度为10mg/L,假设自然水体当中石油类农药的浓度为0。

1.1.2作物生产生态效率评价方法

SFA由Aigner等[36]以及Meeusen和Broeck[37]提出,随后,Reifschneider和Stevenson[38]、Kumbhakar等[39]以及Battese和Coelli[40-41]对随机前沿生产函数形式和方法进行了改进,为其后的实证研究建立了模型方法的数理基础。参考Battese和Coelli[41]模型形式,本文建立如下随机前沿生产函数模型：

Yit=f(Xit)exp(vit-uit)=f(Kit,Lit,GWFit,BWFit,EWFit,T)exp(vit-uit)
uit=Zitδ+ωit

(8)

式中,Yit为第i个区域第t年的农业增加值；Kit为资本投入,用农业机械总动力代表；Lit为劳动投入,用农业从业人口代表；T为时间变量。变量的选取参考潘丹和应瑞瑶[24],郑德凤等[30]和侯孟阳和姚顺波[26]等人的研究,除了传统农业生产中资本、劳动的投入外,该模型将农业生产中的资源投入及环境影响考虑在内,因此将作物生产水足迹作为投入指标。

(9)

式中,TEit为作物生态生产的技术效率即生态效率。当uit=0时,TEit=1,表示没有效率损失,如果uit>0,0

本文选择采用形式最为灵活,可近似反映任何生产技术的超越对数(Translog)生产函数作为随机前沿模型的具体形式,对数形式的函数式如下：

(10)

具体的技术欠效率测度模型定义为：

uit=δ0+δ1FSit+δ2PGit+δ3ISit+δ4GSit+δ5T

(11)

式中,FS为第一产业产值占比；PG为人均生产总值；IS为水田和水浇地面积占耕地面积比例；GS为粮食作物面积占总种植面积比例。欠效率测度模型变量的选择参考程翠云等[21]和洪开荣等[23]等人的研究,同时增加考虑了水资源要素。

式(10)和(11)中的未知参数用最大似然估计法联合估计得到。令：

(12)

1.2 数据来源及处理

本文选取1985—2018年陕西省各区域面板数据为研究对象,通过加总各农田作物和经济林果生产水足迹得到陕西省各区域作物生产水足迹,并估算作物生产生态效率。农业增加值、农业机械总动力、农业从业人口、各农作物播种面积、各经济林果种植面积、氮肥施用折纯量、地区生产总值、第一产业产值、人均生产总值、水田、水浇地、耕地面积等数据主要来自《陕西统计年鉴》(1986—2019年),部分数据来自陕西各市历年统计年鉴,缺失值采用插值法和比例法填补。农业增加值和人均生产总值按照1985年不变价格折算。各气象站的月最低气温、月最高气温、平均风速、平均相对湿度、日照时数、降水量等数据来源于中国气象科学数据共享服务网(http://cdc.cma.gov.cn/home.do)。

2 结果分析

2.1 陕西省作物生产水足迹测算

1985—2018年陕西省作物生产水足迹变化情况如图1所示。研究期内陕西省作物生产总水足迹呈上升趋势,从170.66亿m3增长到381.92亿m3,增长了2.24倍。同期,陕西省的作物总种植面积也缓慢增长,从4856. 27千hm2增长到5361.65千hm2,增长1.10倍,说明陕西省作物生产单位面积用水强度(水足迹角度)不断增强。

水足迹衡量的是作物生产中所需要的水资源总量而非实际用水量,由图1可知陕西省作物生产水足迹的增长主要在于灰水足迹部分,说明陕西省作物生产的环境影响呈现快速增长趋势。1985年,陕西省绿、蓝、灰水足迹分别为103.47亿m3、34.16亿m3和33.03亿m3,分别占比60.63%、20.02%和19.35%。而在2018年,绿、蓝、灰水足迹分别为165.95亿m3、62.05亿m3和153.92亿m3,分别占比43.45%、16.25%和40.30%。绿、蓝、灰水足迹年均增长率分别为1.83%、2.47%和11.09%。

图1 陕西省作物生产水足迹变动趋势Fig.1 Water footprint of crop production in Shaanxi Province

2.2 陕西省作物生产生态效率评价

2.2.1模型检验

使用软件Frontier Version 4.1对式(10)和(11)进行估计以测算陕西省作物生产生态效率。模型中各变量的描述性统计见表1,模型估计结果见表2。

表1 变量的描述性统计

由表2可知,γ估计值为0.7927,显著大于0,接近于1,并在1%的显著水平通过了t检验,说明有79.27%的随机误差来自技术欠效率的影响,该模型是整体有效的。同时进行了一系列的检验以确定随机前沿函数是否适合以及哪种具体的函数形式更优,结果见表3。表3显示所有备择假设均被拒绝,可知本文的随机前沿生产函数形式是合理有效的。

表2 陕西省农业随机前沿生产函数模型估计结果

表3 模型检验结果

2.2.2陕西省作物生产生态效率测算

1985—2018年陕西省各区域作物生产生态效率平均值为0.7477(见图2),年际间的波动具体可分为两个阶段：1985—2003年,作物生产生态效率波动下跌,历史最高值在1985年,为0.8729,在2003年达到最低值0.6324；此后进入波动增长期,2018年效率值为0.8077,表明目前陕西省作物生产实际产出离前沿产出边界还有不小距离,作物生产生态效率有一定的提升空间。

图2 陕西省作物生产生态效率变动趋势(1985—2018年) Fig.2 Crop production eco-efficiency of Shaanxi Province from 1985 to 2018

目前中国作物生产生态效率损失的原因并不在于产出不足,而在于资源消耗过多和环境污染物排放过多[24]。陕西省作物生产生态效率整体下跌的最主要原因是灰水投入的大幅增长(见图1),陕西省耕地面积在2003年之前持续减少(从1985年的3627.07千hm2减少至2003年的2795.82千hm2),耕地的减少限制了农业产值的增长,但同期的化肥、农药施用量却持续增长。尤其在2003之前数年,国内正处于“以地生财”的征地高峰,西部大开发促使陕西省耕地面积大幅减少[42],单位面积化肥施用折纯量持续上升(从最低点1996年的8594.35kg/hm2增长至2003年的9592.95kg/hm2),同时,2003年高温、洪涝等气候灾害造成作物大幅度减产[19],也是陕西省作物生产生态生态效率出现极低值的重要原因。之后,国家严格对耕地保护,加大了对农业资源环境的管制力度,陕西省2003年之后几年耕地面积仅微幅下降,并在2008年至今一直处于增长通道,使得作物产量相当程度得以保证的同时,化肥农药的使用效率也不断提高,单位面积化肥施用折纯量在2003年高点后持续降低,2016年,陕西省出台化肥农药使用减量化意见,在一定程度上促进了生态农业的发展。

2.3 陕西省作物生产生态效率影响因素分析

从表2的估计结果可以看出,陕西省作物生产生态效率损失影响因素中,第一产业产值占比、人均生产总值和水浇地占比与效率损失显著负相关,时间变量与效率损失显著正相关。表明,第一产业产值占比、人均生产总值和水浇地占比的提高可以有效的改进作物生产生态效率。过度追求经济效益会使得人们采取粗放的经营模式,而随着经济发展,人们对环境的重视程度增强,不再片面追求农业的经济效益,同时也更有能力控制和治理环境污染,这会促使农业向生态化的方向转变,从而改进作物生产生态效率,这一结果与王宝义和张卫国[28]的研究一致,同时这一格局在农业主产区更易形成。水浇地占比代表了区域的灌溉技术强弱,较大的水浇地占比一方面反映了该区域降水资源相对匮乏,更加重视节约用水,同时也表明该区域有相应的技术保障农业用水,进而促进了农业产值增长。另外,粮食作物面积占比没有通过显著性检验,对作物生产生态效率没有影响。

3 讨论

3.1 陕西省作物生产水足迹区域差异与原因分析

1985—2018年陕西省各区域年均作物生产水足迹如图3所示。渭南市的作物生产总水足迹最大,为59.93.86亿m3,其次是咸阳市(51.07亿m3)和榆林市(35.26亿m3)。杨凌示范区的作物生产水足迹最小,为0.65亿m3,其次是铜川市(5.81亿m3)和商洛市(13.08亿m3)。陕西省水足迹区域差异明显,绿水占比从大到小的区域排序依次为安康市、延安市、榆林市、商洛市、铜川市、汉中市、宝鸡市、咸阳市、西安市、渭南市和杨凌示范区。蓝水占比从大到小的区域排序依次为杨凌示范区、西安市、渭南市、咸阳市、宝鸡市、榆林市、汉中市、铜川市、商洛市、安康市和延安市。灰水占比从大到小的区域排序为渭南市、杨凌示范区、西安市、咸阳市、宝鸡市、铜川市、商洛市、汉中市、榆林市、延安市和安康市。

图3 陕西省各区域年均作物生产水足迹Fig.3 Regional annual average water footprint related to crop production of Shaanxi Province

区域水足迹的大小与农作物种植面积有密切关系,陕西省单位面积水足迹最大的区域是杨凌示范区(7939.01m3/hm2),其次是渭南市(7333.65m3/hm2)和咸阳市(6962.68m3/hm2)。单位面积水足迹最小的三个区域分别是榆林市(3857.08m3/hm2),延安市(4251.01m3/hm2)和商洛市(4322.03m3/hm2)。单位面积水足迹与气候、种植结构和农业污染相关。其中,气候因素中光照和降雨量对作物的生长发育和产量产生直接影响,日照时数和风速等通过影响相对湿度和作物冠层水汽扩散速率间接影响作物生长,从而影响作物单位面积水足迹；不同作物的虚拟水含量不同,作物种植结构调整将导致水足迹的变化；此外,农业污染物的数量变化导致稀释所需的淡水量变化。

陕西省各区域种植结构、降雨量与单位面积化肥、农药施用量如图4所示。其中,低耗水作物、中耗水作物和高耗水作物按照前文计算得到的陕西省各作物需水量进行分类,低耗水作物主要包括玉米、薯类和小杂粮,中耗水作物主要包括小麦、大豆和油料,高耗水作物主要包括经济林果、蔬菜、瓜类、稻谷和棉花。整体而言陕西省的降水量从北到南递增,陕北地区的延安市和榆林市,以及陕南地区的安康市和商洛市灌溉条件限制,是其蓝水占比低,绿水占比高的主要原因。区域比较可以发现,陕北的榆林市降雨资源最少,种植结构中低耗水作物面积占比最大,匹配度较高,而安康市和商洛市虽拥有较丰富的降水资源,但其作物种植集中在中、低耗水却低产值的粮食作物,存在绿水资源利用不足的问题。关中地区渭南市、杨凌示范区、西安市和咸阳市较高的单位面积化肥和农药投入是其灰水占比高的主要原因。其中比较特殊的是,与其他区域不同,渭南市化肥农药的大量使用并没有伴随农业产值的增长,同时渭南市较好的灌溉条件并没有被充分利用以生产高耗水高产值经济作物,导致资源与种植结构不匹配的效率低下的问题。

图4 陕西省各区域种植结构、降雨量与单位面积化肥、农药施用量Fig.4 Regional planting structure, rainfall, and application amount of fertilizer and pesticide per unit area of Shaanxi Province

3.2 陕西省作物生产生态效率区域差异与变化分析

如图5所示,1985—2018年间,作物生产生态效率平均值最高的区域是西安市,效率值为0.9109,其次为咸阳市(0.8934)和杨凌示范区(0.8901)。作物生产生态效率最低的区域是铜川市,效率值为0.4644,其次为商洛市(0.6220)和安康市(0.6843)。为进一步分析陕西省作物生产生态效率区域差异与年际变化,选择作物生产生态效率变化的关键节点年份(1985,1993,1998,2003,2010和2018年),陕西各区域关键节点年份作物生产生态效率分布见图6。

图5 陕西省各区域作物生产生态效率箱图(1985—2018年)Fig.5 Box diagram of eco-efficiency of crop production in various regions of Shaanxi Province (1985—2018)

图6 陕西省关键节点年份作物生产生态效率分布Fig.6 Distribution of crop production eco-efficiency in key node years of Shaanxi Province

如前文所述,从图6中也可看出陕西省整体作物生产生态效率,以2003年为分界点,经历了先下降后上升的U型曲线,同时也存在一定的区域差异。其中陕北地区的延安市和榆林市、关中地区西部的宝鸡市、陕南地区西部的汉中市呈现类似的 U型曲线,而关中地区中部的咸阳市、西安市和杨凌示范区一直维持较高的生态效率,陕南地区东部的安康市和商洛市,以及关中地区东部的渭南市生态效率持续下降至2003年后长时间处于低效率区间,尤其是商洛市至今仍无恢复迹象,最值得关注的是铜川市,从1985年至今,其作物生产生态效率处于持续低效率水平。

比较各投入要素效率可以发现,整体而言,铜川市各投入效率均属于低效率区域,尤其是农业机械动力、绿水足迹和灰水足迹。其中,铜川市单位农业机械动力增加值为0.049元/kW,约为效率最高的咸阳市的50%,单位绿水和灰水增加值分别为0.420元/m3和0.601元/m3,约为效率最高的杨凌示范区的20%和38%。渭南市除绿水效率较高外,其他要素效率均属于低效率区域,尤其是农业机械动力、蓝水足迹和灰水足迹。其中,单位农业机械动力增加值为0.044元/kW,单位蓝水增加值为0.866元/m3,仅约为效率最高的延安市的22%,单位灰水增加值为0.513元/m3。商洛市的低效率主要表现在劳动力要素上,单位劳动力增加值为0.067元/人,约为效率最高的杨凌示范区的37%。安康市的低效率主要表现在绿水足迹上,单位绿水增加值为0.342元/人。

具体分析各低生态效率区域影响因素发现,铜川市作为一个工业市,其农业人口占比和一产非常低,同时其农业耕地和灌溉条件差,水浇地和水田占比排名倒数第一,是其农业生产落后,进而农业生态效率水平低下的根本原因。渭南市灌溉条件较好,但经济整体落后,低效率(产值角度)的灌溉和化肥农药施用导致较低的生态效率。商洛市和安康市比较类似,农业人口占比和一产占比均排在陕西各区域前列,然而灌溉条件的匮乏和落后的经济水平以及农业生产方式严重阻碍作物生产生态效率的提升。

3.3 陕西省作物生产各投入的产出弹性分析

基于陕西省生态农业随机前沿模型模拟结果,进一步分析了各投入的产出弹性,结果见图7。从投入的产出弹性总和值的变化可以看出,陕西省作物生产经历了从规模报酬递增(1985—2000年),到规模报酬递减(2001—2018年)2个阶段,与耕地面积的变化基本对应。2018年陕西省作物生产各投入产出弹性总和值为0.904,说明考虑资源投入和环境影响时,在目前技术水平和产业组织下,作物生产规模扩大带来的产出增长在缩小。

图7 陕西省作物生产各投入的产出弹性Fig.7 The elasticities of output with respect to the inputs of crop production in Shaanxi Province

各弹性中,资本的产出弹性连年增长,目前是拉动产出增长最有力的要素。劳动投入的产出弹性不断降低,在-1.1以下,农业劳动力投入与产出增长负相关,原因在于,一方面,单纯使用劳动力数量无法准确衡量作物生产中的劳动投入,另一方面,农业的升级发展并不依赖劳动投入的数量,专业的高素质农业劳动力是农业发展的关键。目前,陕西省存在大量冗余的农业劳动力。三种水足迹投入中,蓝、绿水产出弹性均处于正值区域,绿水的产出弹性波动下降,蓝水的产出弹性增长至1995年后长期稳定,说明目前陕西省作物生产无论在自然水资源还是灌溉水资源上都还存在缺口。灰水的产出弹性也处于正值区域,但呈现下降趋势,说明目前化肥、农药等投入对产出仍有很大贡献,但其作用不断弱化,造成的环境影响不容小觑。

具体到陕西省各区域三种水足迹的产出弹性分布如图8所示。榆林市、渭南市和西安市的绿水产出弹性由正(1985—2003年)转为负(2004—2018年),表明这些区域作物种植对绿水资源已经过了最佳匹配点,增加消耗绿水作物的种植对产值增长无益。其中,西安市和渭南市拥有较高的蓝水产出弹性,在这些区域增加灌溉投入可以有效地提高农业产出,而榆林市拥有较高的灰水产出弹性,说明榆林市的化肥农药使用还处于高效率区间,可考虑加大使用。铜川市和商洛市都有较高的绿水产出弹性,如果这两个区域目前的作物生产没有充分地利用降雨资源,种植结构和产值结构均需要调整。延安市、安康市和榆林市的蓝水产出弹性相比较小,在这些区域增加灌溉投入的效果较弱。铜川市和杨凌示范区的灰水产出弹性由正(1985—2003年)变为负(2004—2018年),说明肥料和农药的投入已经对产出增长无益,同时造成环境污染,需要严格控制。

图8 陕西省各区域作物生产水足迹的产出弹性Fig.8 The regional output elasticities of water footprints of crop production in Shaanxi Province

4 结论

本文将绿水、蓝水和灰水足迹作为投入要素引入作物生产函数,分别代表降水资源、灌溉水资源和作物生产的环境影响,通过建立随机前沿作物生产函数模型,测算1985—2018年陕西省各区域的作物生产生态效率,并联合估计第一产业产值占比、人均生产总值、水浇地占比和粮食作物面积占比等因素对生态效率的影响。得到以下研究结论：

(1) 1985—2018年间陕西省作物生产总水足迹呈上升趋势,作物生产单位面积用水强度不断增强。总水足迹中,绿水、蓝水、灰水均有所增长,灰水的增长最为显著。单位面积水足迹及蓝、绿、灰水占比与气候、种植结构,灌溉条件和农业污染相关。

(2) 研究期内,陕西省作物生产生态效率以2003年为分界点呈现先降后升区域,2018年效率值为0.8077,表明目前陕西省作物生产实际产出离前沿产出边界还有距离,作物生产生态效率有一定的提升空间。陕西省作物生产生态效率的U型波动对应政府对耕地资源保护和农业污染管控政策变化。

(3) 第一产业产值占比、人均生产总值和水浇地占比的提高可以有效的改进作物生产生态效率,粮食作物面积占比对作物生产生态效率没有显著影响。作物生产生态效率低下的区域集中在生产落后的农业区和农业条件差的工业区,农业耕地质量差,灌溉条件匮乏,低效率(产值角度)的灌溉和化肥农药施用以及落后农业生产方式是导致作物生产生态效率低下的主要原因。

(4) 研究期内,陕西省作物生产经历了从规模报酬递增,到规模报酬递减2个阶段,与耕地面积的变化基本对应。2018年陕西省作物生产各投入产出弹性总和值为0.904,说明考虑资源投入和环境影响时,在目前技术水平和产业组织下,作物生产规模扩大带来的产出增长在缩小。短期内,考虑各投入要素的产出弹性调整种植结构和化肥农药施用量是提高作物生产生态效率的有效途径,而促进农业产业升级和技术革新是推动作物生产生态效率提升的根本原因。