中國城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響*——基于中國13個(gè)糧食主產(chǎn)區(qū)2009—2018年面板數(shù)據(jù)

尚 杰, 吉雪強(qiáng), 陳璽名

中國城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響*——基于中國13個(gè)糧食主產(chǎn)區(qū)2009—2018年面板數(shù)據(jù)

尚 杰, 吉雪強(qiáng), 陳璽名

(東北林業(yè)大學(xué)經(jīng)濟(jì)管理學(xué)院 哈爾濱 150040)

城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率會(huì)產(chǎn)生重要影響。為更全面和準(zhǔn)確分析城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響, 本文基于中國13個(gè)糧食主產(chǎn)區(qū)2009—2018年面板數(shù)據(jù), 利用3階段非期望產(chǎn)出SBM-DEA模型、熵值法、Tobit回歸等多種模型和方法, 在理論分析基礎(chǔ)上從多個(gè)角度就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行研究。結(jié)果表明: 1)外生環(huán)境影響和隨機(jī)干擾對農(nóng)業(yè)生態(tài)效率測度會(huì)產(chǎn)生顯著影響, 并進(jìn)一步作用于城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響。2)中國13個(gè)糧食主產(chǎn)區(qū)2009—2018年間農(nóng)業(yè)生態(tài)效率呈現(xiàn)持續(xù)增長趨勢, 剔除外生環(huán)境影響和隨機(jī)干擾后, 年度效率均值由2009年的0.53提升至2018年的0.80; 但不同地區(qū)間存在較大差距, 2018年河北等地區(qū)農(nóng)業(yè)生態(tài)效率達(dá)到前沿面, 但是內(nèi)蒙古農(nóng)業(yè)生態(tài)效率只有0.45。3)城鎮(zhèn)化發(fā)展整體上能夠促進(jìn)農(nóng)業(yè)生態(tài)效率的提高, 城鎮(zhèn)化綜合指數(shù)在1%的水平上顯著正向影響農(nóng)業(yè)生態(tài)效率, 影響系數(shù)為0.60。4)城鎮(zhèn)化各指標(biāo)對農(nóng)業(yè)生態(tài)效率的影響存在差異, 從各指標(biāo)情況來看, 人口城鎮(zhèn)化率負(fù)向影響農(nóng)業(yè)生態(tài)效率, 但不顯著; 城鎮(zhèn)居民人均可支配收入、城鎮(zhèn)經(jīng)濟(jì)密度正向影響農(nóng)業(yè)生態(tài)效率; 二三產(chǎn)業(yè)產(chǎn)值占GDP比重、建成區(qū)面積比重、人均建成區(qū)面積等指標(biāo)負(fù)向影響農(nóng)業(yè)生態(tài)效率。為進(jìn)一步發(fā)揮城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的促進(jìn)作用, 本文根據(jù)研究結(jié)果從人口、經(jīng)濟(jì)、土地城鎮(zhèn)化3個(gè)方面提出了政策建議。

城鎮(zhèn)化; 農(nóng)業(yè)生態(tài)效率; 糧食主產(chǎn)區(qū); DEA; Tobit

城鎮(zhèn)化問題是國家發(fā)展中的重要問題, 《2019年國務(wù)院政府工作報(bào)告》《2019年新型城鎮(zhèn)化建設(shè)重點(diǎn)任務(wù)》等多份政府文件都體現(xiàn)了對城鎮(zhèn)化的關(guān)注, 強(qiáng)調(diào)城鄉(xiāng)協(xié)調(diào)發(fā)展的新型城鎮(zhèn)化是中國城鎮(zhèn)化的主要方向。城鎮(zhèn)化離不開鄉(xiāng)村的支持, 在中國城鎮(zhèn)化速度加快的同時(shí), 農(nóng)村空心化、農(nóng)村生態(tài)環(huán)境惡化等問題逐漸突顯[1-2]。農(nóng)業(yè)生態(tài)環(huán)境改善是鄉(xiāng)村振興的重要內(nèi)容, 更是生態(tài)文明理念貫徹落實(shí)的重要實(shí)踐。如何推動(dòng)新型城鎮(zhèn)化與農(nóng)業(yè)生態(tài)環(huán)境改善相協(xié)調(diào)是中國農(nóng)村發(fā)展過程中的重要問題。通過研究城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響, 可以更有針對性地為農(nóng)業(yè)生態(tài)環(huán)境改善提出建議, 從而為城鄉(xiāng)生態(tài)環(huán)境協(xié)調(diào)發(fā)展提供理論支撐。

從現(xiàn)有研究來看, 城鎮(zhèn)化是世界范圍的研究熱點(diǎn), 而農(nóng)業(yè)生態(tài)效率作為農(nóng)業(yè)生態(tài)環(huán)境狀況的重要反映得到了諸多學(xué)者的認(rèn)可。城鎮(zhèn)化對農(nóng)業(yè)農(nóng)村發(fā)展的重要影響受到了學(xué)者們的關(guān)注, Westlund[3]、Deng等[4]和Zhang等[5]就城鎮(zhèn)化對農(nóng)村農(nóng)業(yè)發(fā)展的負(fù)面影響進(jìn)行了探討; Andersen[6]、Zhou等[7]就城鎮(zhèn)化對農(nóng)業(yè)生產(chǎn)和農(nóng)民生活水平的促進(jìn)作用進(jìn)行了分析。農(nóng)業(yè)生態(tài)效率是農(nóng)業(yè)可持續(xù)發(fā)展能力的重要指標(biāo), 諸多學(xué)者利用農(nóng)業(yè)生態(tài)效率這一指標(biāo)對農(nóng)業(yè)生態(tài)環(huán)境問題進(jìn)行了探討。在農(nóng)業(yè)生態(tài)效率的測度方面, 大部分學(xué)者選擇了無需預(yù)設(shè)函數(shù)關(guān)系、能夠減少研究主觀爭議性的數(shù)據(jù)包絡(luò)分析模型, 現(xiàn)有研究中學(xué)者們分別使用偏好錐的DEA模型[8]、SBM-DEA模型[9]、SBM-Undesirable擴(kuò)展模型[10]、超效率DEA模型[11]等DEA模型對農(nóng)業(yè)生態(tài)效率進(jìn)行了測度, 取得了一定成果, 但是這些DEA模型卻未能考慮外生環(huán)境影響和隨機(jī)干擾對農(nóng)業(yè)生態(tài)效率測度所產(chǎn)生的影響。

在農(nóng)業(yè)生態(tài)效率影響因素研究方面, 大多數(shù)研究立足于農(nóng)業(yè)生產(chǎn)環(huán)境或者農(nóng)村發(fā)展?fàn)顩r, 較少探討對農(nóng)業(yè)農(nóng)村發(fā)展有著重要影響的城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響。侯孟陽等[12]分析了農(nóng)村勞動(dòng)力轉(zhuǎn)移對農(nóng)業(yè)生態(tài)效率的影響; 王寶義等[13]探討了人均農(nóng)業(yè)增加值、農(nóng)業(yè)規(guī)?；健⑥r(nóng)業(yè)受災(zāi)率、農(nóng)業(yè)機(jī)械密度、農(nóng)民家庭經(jīng)營收入比等因素與農(nóng)業(yè)生態(tài)效率之間的關(guān)系; 胡平波等[11]研究了農(nóng)旅融合對農(nóng)業(yè)生態(tài)效率的影響; 曹俊文等[14]分析了農(nóng)業(yè)經(jīng)濟(jì)水平、農(nóng)業(yè)公共投資、農(nóng)業(yè)勞動(dòng)力教育程度、農(nóng)業(yè)機(jī)械化、政府規(guī)制和技術(shù)進(jìn)步等因素對農(nóng)業(yè)生態(tài)效率的影響。

就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率產(chǎn)生的影響進(jìn)行研究,可以為中國農(nóng)業(yè)生態(tài)環(huán)境改善提供有效建議。然而從現(xiàn)有研究來看, 農(nóng)業(yè)生態(tài)效率影響因素的研究有所局限, 較少有研究就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行深入探討。即使少部分研究涉及到城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響, 也多從單一人口城鎮(zhèn)化角度進(jìn)行分析, 忽略了經(jīng)濟(jì)城鎮(zhèn)化及土地城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響[15]。通過對現(xiàn)有研究的梳理, 筆者發(fā)現(xiàn)以下特征: 1)在農(nóng)業(yè)生態(tài)效率測度上大多數(shù)研究采取了DEA方法, 但是少有研究考慮外生環(huán)境影響和隨機(jī)干擾; 2)專門就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響的研究較為稀少, 尤其缺乏全方位深入探討城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響的研究, 且現(xiàn)有研究涉及的農(nóng)業(yè)生態(tài)效率測度較多受到外生環(huán)境影響, 難以真實(shí)反映城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響情況。因此, 本文在現(xiàn)有研究基礎(chǔ)上, 首先對城鎮(zhèn)化、農(nóng)業(yè)生態(tài)效率、城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行初步理論分析, 其次對相關(guān)模型和方法、變量設(shè)定及數(shù)據(jù)來源進(jìn)行介紹, 然后結(jié)合2009—2018年中國13個(gè)糧食主產(chǎn)區(qū)實(shí)際數(shù)據(jù)進(jìn)行農(nóng)業(yè)生態(tài)效率、城鎮(zhèn)化綜合指數(shù)測度, 并就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行回歸分析, 最后就相關(guān)結(jié)論進(jìn)行總結(jié)并提出建議。相對現(xiàn)有研究本文具有以下創(chuàng)新: 1)使用3階段非期望產(chǎn)出SBM-DEA模型測度農(nóng)業(yè)生態(tài)效率, 在測度過程中減少外生環(huán)境影響和隨機(jī)干擾對農(nóng)業(yè)生態(tài)效率測度的影響, 提高研究準(zhǔn)確性; 2)從人口、經(jīng)濟(jì)、土地城鎮(zhèn)化3個(gè)方面展開城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響研究, 較現(xiàn)有研究更為深入且更為全面, 由于本文所探討的農(nóng)業(yè)生態(tài)效率是剔除外生環(huán)境影響及隨機(jī)干擾后的效率值, 其結(jié)果更能真實(shí)反映城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的作用。

1 理論分析

1.1 城鎮(zhèn)化及城鎮(zhèn)化指標(biāo)

城鎮(zhèn)化, 也可以稱為城市化, 大多數(shù)學(xué)者從人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化、土地城鎮(zhèn)化3個(gè)方面進(jìn)行城鎮(zhèn)化研究[16-22]。本文同樣從人口、經(jīng)濟(jì)、土地城鎮(zhèn)化對城鎮(zhèn)化進(jìn)行研究, 以便于更具體地研究城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響。本文在參考現(xiàn)有研究基礎(chǔ)上以人口城鎮(zhèn)化率、城鎮(zhèn)居民人均可支配收入作為人口城鎮(zhèn)化狀況的替代指標(biāo)來反映農(nóng)村人口流入城市情況以及城鎮(zhèn)居民生活水平, 以城鎮(zhèn)經(jīng)濟(jì)密度、第二三產(chǎn)業(yè)占GDP比重作為經(jīng)濟(jì)城鎮(zhèn)化的替代指標(biāo)反映城市產(chǎn)業(yè)集中情況以及地區(qū)產(chǎn)業(yè)結(jié)構(gòu), 以建成區(qū)面積比重、人均建成區(qū)面積作為土地城鎮(zhèn)化的替代指標(biāo)反映城市土地的擴(kuò)張情況, 具體指標(biāo)設(shè)置如表1所示。并且利用熵值法對6個(gè)指標(biāo)進(jìn)行計(jì)算以得到城鎮(zhèn)化綜合指數(shù), 從而更全面地反映城鎮(zhèn)化整體發(fā)展情況。

表1 城鎮(zhèn)化綜合指數(shù)評價(jià)指標(biāo)體系

1.2 農(nóng)業(yè)生態(tài)效率及其指標(biāo)體系

農(nóng)業(yè)生態(tài)效率是指在一定的農(nóng)業(yè)投入要素組合下, 以盡可能小的資源消耗和環(huán)境污染, 得到盡可能多的農(nóng)業(yè)產(chǎn)出[23]。考慮到農(nóng)業(yè)勞動(dòng)力是農(nóng)業(yè)生產(chǎn)的主要力量, 土地資源是農(nóng)業(yè)活動(dòng)開展的重要基礎(chǔ), 而化肥、農(nóng)藥、農(nóng)膜、農(nóng)業(yè)機(jī)械動(dòng)力是農(nóng)業(yè)實(shí)踐中重要生產(chǎn)資料, 且農(nóng)業(yè)生產(chǎn)離不開對水資源的利用, 為此本文選取勞動(dòng)力、化肥、農(nóng)藥、農(nóng)膜、農(nóng)業(yè)機(jī)械動(dòng)力、灌溉6項(xiàng)指標(biāo)作為農(nóng)業(yè)生態(tài)效率的投入指標(biāo)?？紤]農(nóng)業(yè)總產(chǎn)值是農(nóng)業(yè)生產(chǎn)成果的重要反映, 而農(nóng)業(yè)生產(chǎn)過程碳排放和化肥、農(nóng)藥、農(nóng)膜等生產(chǎn)資料流失是造成農(nóng)業(yè)生態(tài)環(huán)境惡化的重要因素, 在結(jié)合王寶義等[10,13]、侯孟陽等[12]學(xué)者研究的基礎(chǔ)上, 選擇農(nóng)業(yè)總產(chǎn)值作為農(nóng)業(yè)生態(tài)效率期望產(chǎn)出指標(biāo), 選擇碳排放和污染排放作為農(nóng)業(yè)生態(tài)效率非期望產(chǎn)出指標(biāo), 構(gòu)建考慮非期望產(chǎn)出的農(nóng)業(yè)生態(tài)效率投入-產(chǎn)出指標(biāo)體系, 見表2。由于農(nóng)業(yè)生態(tài)環(huán)境建設(shè)會(huì)受到地區(qū)宏觀經(jīng)濟(jì)水平、政府政策支持情況、地區(qū)科技發(fā)展情況等宏觀外部環(huán)境的干擾, 而Fried認(rèn)為利用3階段DEA模型可以很好地解決效率值測度中外部環(huán)境干擾的問題[24]。為此本文使用3階段非期望產(chǎn)出SBM-DEA模型進(jìn)行農(nóng)業(yè)生態(tài)效率測度, 以提升農(nóng)業(yè)生態(tài)效率測度準(zhǔn)確性。

表2 農(nóng)業(yè)生態(tài)效率投入-產(chǎn)出指標(biāo)體系

1.3 城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響理論分析

城鎮(zhèn)化對農(nóng)業(yè)農(nóng)村發(fā)展會(huì)產(chǎn)生重要的影響, 同樣也會(huì)對農(nóng)業(yè)生態(tài)效率產(chǎn)生重要影響, 下面結(jié)合現(xiàn)有研究成果與社會(huì)實(shí)際就城鎮(zhèn)化各方面對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行初步理論分析。

1)人口城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率產(chǎn)生的影響。隨著人口城鎮(zhèn)化水平的提高, 大量勞動(dòng)力離開農(nóng)村, 留在農(nóng)村的勞動(dòng)力不得不使用較多的農(nóng)業(yè)機(jī)械以彌補(bǔ)勞動(dòng)力的減少[25], 而農(nóng)業(yè)機(jī)械使用數(shù)量的增多將造成農(nóng)業(yè)生態(tài)環(huán)境惡化; 而另一方面城鎮(zhèn)化對于促進(jìn)農(nóng)民環(huán)保行為改善具有積極效應(yīng)[26], 這種環(huán)保觀念有可能通過入城農(nóng)民與留守農(nóng)民之間的聯(lián)系而影響留守農(nóng)民在農(nóng)業(yè)生產(chǎn)中的行為選擇, 從而減少農(nóng)業(yè)生產(chǎn)中化學(xué)制品的應(yīng)用。當(dāng)城鎮(zhèn)居民人均可支配收入提升時(shí)會(huì)推動(dòng)居民消費(fèi)水平升級, 而居民消費(fèi)水平升級可能增加對綠色農(nóng)產(chǎn)品的需求, 綠色消費(fèi)增加可以為綠色農(nóng)業(yè)產(chǎn)業(yè)化發(fā)展提供巨大的市場空間[27], 從而使得農(nóng)民生產(chǎn)時(shí)重視農(nóng)業(yè)生態(tài)環(huán)境保護(hù), 減少化肥、農(nóng)藥等化學(xué)制品的使用。

2)經(jīng)濟(jì)城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率產(chǎn)生的影響。鄉(xiāng)鎮(zhèn)企業(yè)向城鎮(zhèn)的聚集將出現(xiàn)大量財(cái)富向城鎮(zhèn)的直接轉(zhuǎn)移[28], 而城鎮(zhèn)經(jīng)濟(jì)聚集也將進(jìn)一步促使農(nóng)村生產(chǎn)要素加速向城市流動(dòng), 農(nóng)業(yè)生產(chǎn)要素流出會(huì)影響到農(nóng)業(yè)生產(chǎn), 如城市周邊農(nóng)業(yè)用地的變更對農(nóng)業(yè)生產(chǎn)會(huì)產(chǎn)生負(fù)外溢效應(yīng)[29]。二三產(chǎn)業(yè)的發(fā)展對農(nóng)業(yè)生態(tài)環(huán)境的影響也是多方面的, 工業(yè)發(fā)展能夠豐富農(nóng)業(yè)生產(chǎn)資料, 在一定程度上會(huì)降低化肥等產(chǎn)品使用成本, 使得農(nóng)戶在農(nóng)業(yè)生產(chǎn)中使用更多化肥、農(nóng)藥, 從而對農(nóng)業(yè)生態(tài)環(huán)境造成不良影響; 但是工業(yè)發(fā)展也能提高農(nóng)業(yè)機(jī)械動(dòng)力投入, 從而提升農(nóng)業(yè)生產(chǎn)效率, 但是農(nóng)業(yè)機(jī)械投入的增多又會(huì)增加柴油等能源消耗[13], 從而加大農(nóng)業(yè)生產(chǎn)碳排放; 此外信息產(chǎn)業(yè)的發(fā)展可以推動(dòng)農(nóng)業(yè)信息化建設(shè)以提高農(nóng)業(yè)生產(chǎn)能力, 農(nóng)業(yè)信息化水平可以提高農(nóng)業(yè)綠色全要素生產(chǎn)率[30]。

3)土地城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率產(chǎn)生的影響。土地城鎮(zhèn)化對耕地可持續(xù)利用會(huì)產(chǎn)生嚴(yán)重的負(fù)面影響[31], 城市建成區(qū)面積的擴(kuò)大將減少農(nóng)業(yè)用地, 可能對農(nóng)業(yè)產(chǎn)出產(chǎn)生不良影響; 但是建成區(qū)面積的增多能夠?yàn)榫G色農(nóng)產(chǎn)品的銷售創(chuàng)造更好的交通和儲(chǔ)存條件, 從而推動(dòng)綠色農(nóng)業(yè)的發(fā)展。人均建成區(qū)面積作為個(gè)人所能享受到城市服務(wù)的重要反映, 當(dāng)人均建成區(qū)面積越大時(shí), 個(gè)人越有可能享受到更多的城市服務(wù), 更易受到城市公益環(huán)保宣傳的影響形成綠色消費(fèi)習(xí)慣, 使得其有更多意愿購買綠色農(nóng)產(chǎn)品。

從以上分析中可以看出城鎮(zhèn)化發(fā)展對農(nóng)業(yè)生態(tài)效率的影響是復(fù)雜的、多方面的, 僅依靠經(jīng)驗(yàn)分析很難把握城鎮(zhèn)化與農(nóng)業(yè)生態(tài)效率之間的關(guān)系, 為了更準(zhǔn)確了解城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響, 需要進(jìn)行進(jìn)一步的量化分析。

2 方法介紹、變量選擇與數(shù)據(jù)來源

2.1 方法介紹

本文利用SBM-DEA模型對農(nóng)業(yè)生態(tài)效率進(jìn)行定量測度; 利用熵值法得到城鎮(zhèn)化綜合指數(shù); 利用Tobit回歸方法就城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響進(jìn)行分析。

Tone在提出SBM模型之后, 又提出非期望產(chǎn)出SBM模型[32-33], 如下式所示:

在非期望產(chǎn)出SBM模型基礎(chǔ)上, 構(gòu)建3階段非期望產(chǎn)出SBM-DEA模型。

第1階段, 采用非期望產(chǎn)出SBM模型, 計(jì)算初始效率和松弛;

第2階段, 類似Fried等[24]構(gòu)建SFA模型找出環(huán)境變量并根據(jù)結(jié)果調(diào)整投入產(chǎn)出量?；貧w方程如下:

第3階段: 使用調(diào)整后投入產(chǎn)出數(shù)據(jù)運(yùn)用非期望SBM-DEA模型進(jìn)行效率測算。

2.2 變量選擇

1)被解釋變量: 以農(nóng)業(yè)生態(tài)效率作為被解釋變量, 其投入產(chǎn)出指標(biāo)體系的構(gòu)建如表2所示。

在農(nóng)業(yè)生態(tài)效率的投入產(chǎn)出指標(biāo)體系中包括7項(xiàng)投入指標(biāo)、3項(xiàng)產(chǎn)出指標(biāo), 各項(xiàng)指標(biāo)所使用替代變量如表2所示。需要特別說明的是勞動(dòng)力投入, 本文用農(nóng)業(yè)從業(yè)人員來表征, 結(jié)合前人的做法[13], 用下式計(jì)算農(nóng)業(yè)從業(yè)人員: 農(nóng)業(yè)從業(yè)人員=農(nóng)林牧漁業(yè)從業(yè)人員(以第一產(chǎn)業(yè)從業(yè)人員作為替代指標(biāo))′農(nóng)業(yè)總產(chǎn)值/農(nóng)林牧漁業(yè)總產(chǎn)值, 參考李寧等[36]、黃瑪蘭等[37]的做法將第一產(chǎn)業(yè)從業(yè)人員作為農(nóng)林牧漁業(yè)從業(yè)人員的替代指標(biāo)計(jì)算農(nóng)業(yè)勞動(dòng)力投入。產(chǎn)出指標(biāo)之中的農(nóng)業(yè)期望產(chǎn)出指標(biāo)用農(nóng)業(yè)總產(chǎn)值(可比價(jià)格)來表征。在對非期望產(chǎn)出中的農(nóng)業(yè)面源污染量化計(jì)算中, 參考吳小慶等[8]的做法采用化肥污染量、農(nóng)藥污染量及農(nóng)膜殘留量估算農(nóng)業(yè)面源污染水平?；饰廴玖扛鶕?jù)歷年化肥施用量與化肥利用率數(shù)據(jù)計(jì)算得到, 即化肥污染量=化肥施用量×化肥流失率, 農(nóng)藥污染量和農(nóng)膜殘留率的算法類似; 其中根據(jù)賴斯蕓等[38]、吳小慶等[8]學(xué)者的做法設(shè)置化肥流失率為65%、農(nóng)藥污染率50%、地膜殘留率為10%。本文選取化肥、農(nóng)藥、農(nóng)膜、農(nóng)用柴油、農(nóng)業(yè)灌溉、農(nóng)業(yè)耕作6項(xiàng)指標(biāo)估算農(nóng)業(yè)碳排放量, 采取的方法是相應(yīng)指標(biāo)乘以排放系數(shù)。參考已有研究[13], 設(shè)置排放源排放系數(shù): 化肥0.90、農(nóng)藥4.93、農(nóng)膜5.18、柴油0.59、農(nóng)業(yè)灌溉20.48、農(nóng)業(yè)耕作312.60。

在利用3階段非期望產(chǎn)出SBM-DEA對農(nóng)業(yè)生態(tài)效率進(jìn)行測度時(shí)需要對環(huán)境變量進(jìn)行選擇。結(jié)合中國農(nóng)業(yè)發(fā)展特點(diǎn), 選擇宏觀經(jīng)濟(jì)發(fā)展情況、政府政策支持、技術(shù)支撐3個(gè)方面作為環(huán)境因素。在宏觀經(jīng)濟(jì)發(fā)展方面, 選擇能反映農(nóng)業(yè)繁榮程度的地區(qū)第一產(chǎn)業(yè)產(chǎn)值作為代理變量。在政府對產(chǎn)業(yè)政策支持方面, 選擇地方財(cái)政支農(nóng)支出作為政府對農(nóng)業(yè)發(fā)展政策方面的代理變量。在產(chǎn)業(yè)技術(shù)支撐方面, 考慮到農(nóng)業(yè)建設(shè)受到農(nóng)業(yè)、生物、信息等諸多技術(shù)影響, 選擇研究與發(fā)展內(nèi)部經(jīng)費(fèi)作為技術(shù)投入代理變量。

2)解釋變量: 本文選擇兩組變量作為城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率回歸分析的解釋變量, 以實(shí)現(xiàn)從城鎮(zhèn)化整體到各指標(biāo)對農(nóng)業(yè)生態(tài)效率影響的研究。第一組為城鎮(zhèn)化綜合指數(shù), 該指數(shù)是基于熵值法對表1中人口城鎮(zhèn)化率、城鎮(zhèn)居民人均可支配收入、城鎮(zhèn)經(jīng)濟(jì)密度、第二三產(chǎn)業(yè)占GDP的比重、建成區(qū)面積比重、人均建成區(qū)面積等6項(xiàng)城鎮(zhèn)化指標(biāo)進(jìn)行綜合運(yùn)算后獲得, 以此來反映城鎮(zhèn)化的整體情況。第二組則為城鎮(zhèn)化6項(xiàng)指標(biāo), 即上述人口城鎮(zhèn)化率等6項(xiàng)指標(biāo)。

2.3 數(shù)據(jù)來源

本文以中國13個(gè)糧食主產(chǎn)區(qū)(河北、內(nèi)蒙古、遼寧、吉林、黑龍江、江蘇、安徽、江西、山東、河南、湖北、湖南、四川)為研究對象, 主要是考慮到北京等直轄市農(nóng)業(yè)功能退化, 而西藏等地則以牧業(yè)為主, 為此選擇我國主要糧食產(chǎn)區(qū)進(jìn)行研究, 可以更好地反映城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響。本文基礎(chǔ)數(shù)據(jù)來源于《中國統(tǒng)計(jì)年鑒》《中國農(nóng)村統(tǒng)計(jì)年鑒》《中國人口和就業(yè)統(tǒng)計(jì)年鑒》《中國城市建設(shè)統(tǒng)計(jì)年鑒》《中國科技統(tǒng)計(jì)年鑒》。

3 結(jié)果與分析

3.1 農(nóng)業(yè)生態(tài)效率測度

1)第1階段農(nóng)業(yè)生態(tài)效率測度。本文以DEA．SOLVER PRO5.0為農(nóng)業(yè)生態(tài)效率測度平臺(tái)。Battese等[39]認(rèn)為, 規(guī)模報(bào)酬可變假設(shè)適用于微觀的研究對象。根據(jù)本文研究對象較為宏觀的特點(diǎn), 選擇規(guī)模報(bào)酬不變的非期望產(chǎn)出SBM-DEA模型進(jìn)行農(nóng)業(yè)生態(tài)效率測度, 按照默認(rèn)設(shè)置, 設(shè)定期望產(chǎn)出和非期望產(chǎn)出總權(quán)重均為1, 測度結(jié)果如表3所示。從第1階段測度的農(nóng)業(yè)生態(tài)效率情況來看, 我國主要糧食產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率值不斷提升, 反映出我國農(nóng)業(yè)生態(tài)環(huán)境改善, 這與現(xiàn)有大部分研究成果所測度的結(jié)果較為相似, 顯示第1階段測度結(jié)果較為準(zhǔn)確。在2018年多數(shù)糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率達(dá)到前沿, 顯示出較高的效率水平。

2)第2階段指標(biāo)數(shù)據(jù)調(diào)整。本階段以第1階段測度獲得的關(guān)于農(nóng)業(yè)生態(tài)效率各投入產(chǎn)出指標(biāo)冗余變量為被解釋變量, 以2009—2018年各糧食主產(chǎn)區(qū)第一產(chǎn)業(yè)產(chǎn)值、財(cái)政支農(nóng)支出、研究經(jīng)費(fèi)支出3個(gè)環(huán)境變量為解釋變量, 利用Coelli等[40]給出的隨機(jī)前沿分析軟件Frontier 4.1, 計(jì)算環(huán)境變量對農(nóng)業(yè)生態(tài)效率的各投入產(chǎn)出指標(biāo)變量冗余的影響, 得表4結(jié)果。

從表4可知, 所選主要環(huán)境變量大部分通過5%水平顯著性檢驗(yàn), 同時(shí), 大部分分析結(jié)果的值在顯著的情況下趨近于1, 表明環(huán)境因素和隨機(jī)誤差對農(nóng)業(yè)生態(tài)效率冗余有顯著影響。從表中系數(shù)可知第一產(chǎn)業(yè)產(chǎn)值提升, 化肥投入冗余、農(nóng)藥投入冗余、農(nóng)業(yè)機(jī)械動(dòng)力投入冗余、碳排放產(chǎn)出冗余、面源污染產(chǎn)出冗余顯著增加; 財(cái)政支農(nóng)支出增加時(shí)農(nóng)藥投入冗余與農(nóng)膜投入冗余顯著增加; 科技研究經(jīng)費(fèi)增多時(shí)灌溉投入冗余顯著增加。

3)第3階段農(nóng)業(yè)生態(tài)效率測度。將經(jīng)過調(diào)整的投入產(chǎn)出數(shù)據(jù)再次進(jìn)行測度, 得第3階段農(nóng)業(yè)生態(tài)效率, 見表3。第3階段的農(nóng)業(yè)生態(tài)效率與第1階段未剔除外生環(huán)境影響和隨機(jī)干擾的農(nóng)業(yè)生態(tài)效率相比有較大變動(dòng), 這表明外生環(huán)境影響和隨機(jī)干擾對農(nóng)業(yè)生態(tài)效率測度有明顯影響, 如果不剔除外生環(huán)境影響和隨機(jī)干擾, 農(nóng)業(yè)生態(tài)效率值將受到宏觀環(huán)境等因素干擾, 無法呈現(xiàn)真實(shí)情況, 最終影響研究的準(zhǔn)確性, 因此在農(nóng)業(yè)生態(tài)效率測度中只有剔除外生環(huán)境影響和隨機(jī)干擾, 才能真實(shí)反映農(nóng)業(yè)生態(tài)效率情況, 從而更準(zhǔn)確地為城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響的研究提供有效數(shù)據(jù)支撐。

表3 2009—2018年第1、3階段中國主要糧食產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率

整體上看, 2009—2018年間第1階段與第3階段測度的農(nóng)業(yè)生態(tài)效率整體趨勢表現(xiàn)為上升, 反映出中國糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)環(huán)境的改善, 這與大多數(shù)學(xué)者研究結(jié)果是一致的。大部分糧食主產(chǎn)區(qū)第3階段測度的農(nóng)業(yè)生態(tài)效率高于第1階段, 表明這些糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率受到外部壞境制約。其中, 內(nèi)蒙古、遼寧、吉林、江西4地區(qū)第3階段農(nóng)業(yè)生態(tài)效率低于第1階段, 表明其原來較高的農(nóng)業(yè)生態(tài)效率是得益于外部環(huán)境, 而非效率值真實(shí)體現(xiàn)。就不同地區(qū)效率來看, 地區(qū)之間存在較大差距, 2018年河北、黑龍江、江蘇、山東、河南、湖北、四川等糧食主產(chǎn)區(qū)第3階段效率達(dá)到效率前沿, 而內(nèi)蒙古農(nóng)業(yè)生態(tài)效率只有0.45。在以下分析中將使用剔除了外生環(huán)境影響和隨機(jī)干擾的第3階段農(nóng)業(yè)生態(tài)效率數(shù)據(jù)作為研究基礎(chǔ)。

3.2 城鎮(zhèn)化綜合指數(shù)測度

使用Excel軟件, 利用熵值法結(jié)合人口城鎮(zhèn)化率、城鎮(zhèn)居民人均可支配收入、城鎮(zhèn)經(jīng)濟(jì)密度、第二三產(chǎn)業(yè)占GDP的比重、建成區(qū)面積比重、人均建成區(qū)面積等城鎮(zhèn)化指標(biāo)進(jìn)行城鎮(zhèn)化綜合指數(shù)測度以全面反映城鎮(zhèn)化狀況。首先利用極差法對2009—2018年糧食主產(chǎn)區(qū)城鎮(zhèn)化指標(biāo)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化, 而后利用Excel軟件結(jié)合熵值法對城鎮(zhèn)化6個(gè)指標(biāo)的權(quán)重進(jìn)行測度, 分別為0.11、0.17、0.27、0.07、0.21、0.16。利用各指標(biāo)權(quán)重對標(biāo)準(zhǔn)化后指標(biāo)數(shù)據(jù)進(jìn)行調(diào)整, 得到城鎮(zhèn)化綜合指數(shù)(表5)。2009—2018年中國糧食主產(chǎn)區(qū)城鎮(zhèn)化綜合指數(shù)均值由2009年0.24提升至2018年0.51, 反映出中國城鎮(zhèn)化水平不斷提升, 且整體提升速度較快, 這與中國實(shí)際發(fā)展相符合。從不同地區(qū)情況來看, 城鎮(zhèn)化綜合指數(shù)最高為江蘇, 其2018年城鎮(zhèn)化綜合指數(shù)為0.90, 遠(yuǎn)超過其他省市。

3.3 城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響分析

利用Stata 15軟件, 結(jié)合上文第3階段測度獲得的2009—2018年各糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率與上文測度獲得的城鎮(zhèn)化綜合指數(shù)進(jìn)行Tobit回歸。為了更好地了解城鎮(zhèn)化從整體到各指標(biāo)對農(nóng)業(yè)生態(tài)效率的影響, 以第3階段農(nóng)業(yè)生態(tài)效率為被解釋變量, 進(jìn)行兩次Tobit回歸, 第1次回歸解釋變量為城鎮(zhèn)化綜合指數(shù), 第2次回歸解釋變量為城鎮(zhèn)化6項(xiàng)指標(biāo)。為了驗(yàn)證估計(jì)結(jié)果的穩(wěn)定性, 在Tobit回歸基礎(chǔ)上進(jìn)行Ols估計(jì), 被解釋變量與解釋變量設(shè)置與Tobit回歸一致。模型估計(jì)結(jié)果見表6, 可知Tobit估計(jì)和Ols估計(jì)整體保持一致, 可見估計(jì)結(jié)果具有穩(wěn)健性。

表4 第2階段環(huán)境變量對農(nóng)業(yè)生態(tài)效率冗余變量影響分析結(jié)果

*、**和***表示在10%、5%和1%水平顯著;0為模型截距項(xiàng),1是第一產(chǎn)業(yè)產(chǎn)值參數(shù)估計(jì)值,2是財(cái)政支農(nóng)支出參數(shù)估計(jì)值,3是研究經(jīng)費(fèi)支出參數(shù)估計(jì)值,2是模型中由隨機(jī)噪聲和無效效應(yīng)符合作用造成的偏差,是無效影響因素對個(gè)體效率差異的解釋程度。*, ** and *** represent significant difference at levels of 10%, 5%, 1%, respectively.0 is the model intercept term;1 is the first industry output value parameter estimate;2 is the financial expenditure parameter estimate;3 is the research expenditure parameter estimate.2is the deviation caused by the coincidence of random noise and invalid effect in the model, andis the explanation degree of invalid influence factors to individual efficiency difference.

從表6可知, 城鎮(zhèn)化綜合指數(shù)在1%水平顯著正向影響農(nóng)業(yè)生態(tài)效率, 即隨著城鎮(zhèn)化綜合指數(shù)的提升, 糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率將有所提高。雖然從各指標(biāo)情況來看, 城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響較為復(fù)雜, 但是城鎮(zhèn)化發(fā)展整體上有利于促進(jìn)農(nóng)業(yè)生態(tài)效率的提升。其中城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的正向影響主要得益于城鎮(zhèn)居民人均可支配收入及城鎮(zhèn)經(jīng)濟(jì)密度對農(nóng)業(yè)生態(tài)效率的正向影響。

人口城鎮(zhèn)化率負(fù)向影響農(nóng)業(yè)生態(tài)效率, 但不顯著。這與現(xiàn)有人口城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響的研究結(jié)果有所不同, 主要是由于現(xiàn)有研究在農(nóng)業(yè)生態(tài)效率測度時(shí)未能剔除外生環(huán)境影響及隨機(jī)因素干擾, 農(nóng)業(yè)生態(tài)效率測度準(zhǔn)確性不足, 難以有效反映農(nóng)業(yè)生態(tài)效率真實(shí)情況。在未剔除外生環(huán)境影響及隨機(jī)干擾前, 人口城鎮(zhèn)化率提升促進(jìn)了各地區(qū)宏觀環(huán)境的改善, 如推動(dòng)地區(qū)整體經(jīng)濟(jì)增長, 故人口城鎮(zhèn)化率對農(nóng)業(yè)生態(tài)效率產(chǎn)生了正向影響。在剔除外生環(huán)境影響和隨機(jī)干擾后, 人口城鎮(zhèn)化率對農(nóng)業(yè)生態(tài)效率的真實(shí)影響顯現(xiàn), 由于農(nóng)村勞動(dòng)力向城市流動(dòng), 使得農(nóng)業(yè)勞動(dòng)力減少, 為了減輕勞動(dòng)力不足對農(nóng)業(yè)生產(chǎn)的影響, 農(nóng)業(yè)機(jī)械將獲得更多使用, 而機(jī)械設(shè)備使用所排放的廢棄物將加劇農(nóng)業(yè)生態(tài)環(huán)境污染。

表5 2009—2018年中國糧食主產(chǎn)區(qū)城鎮(zhèn)化綜合指數(shù)

城鎮(zhèn)居民人均可支配收入在1%水平顯著正向影響農(nóng)業(yè)生態(tài)效率, 即隨著城鎮(zhèn)居民人均可支配收入的提高, 農(nóng)業(yè)生態(tài)效率將隨之增長。其原因可能是當(dāng)市民收入提升時(shí), 消費(fèi)水平也將隨之提高, 對農(nóng)產(chǎn)品質(zhì)量的要求也將更高, 有機(jī)食品、綠色食品等生態(tài)農(nóng)產(chǎn)品將更加獲得市場青睞, 在市場利益引導(dǎo)下農(nóng)業(yè)生產(chǎn)部門將隨之調(diào)整, 更加重視農(nóng)業(yè)生產(chǎn)綠色化以迎合市場綠色消費(fèi)需求。

城鎮(zhèn)經(jīng)濟(jì)密度在1%的水平顯著正向影響農(nóng)業(yè)生態(tài)效率, 即隨著二三產(chǎn)業(yè)向城鎮(zhèn)集中, 農(nóng)業(yè)生態(tài)效率將有所提升, 其原因可能是二三產(chǎn)業(yè)聚集于城市而減少了對農(nóng)業(yè)用地的侵占, 從而提升了農(nóng)業(yè)產(chǎn)出, 此外由于二三產(chǎn)業(yè)集中城鎮(zhèn), 其生產(chǎn)過程中所排放的廢棄物對農(nóng)業(yè)生態(tài)環(huán)境的污染減輕。

二三產(chǎn)業(yè)產(chǎn)值占GDP比重在1%的水平上顯著負(fù)向影響農(nóng)業(yè)生態(tài)效率, 即隨著二三產(chǎn)業(yè)比重提升, 農(nóng)業(yè)生態(tài)效率將有所下降, 其原因可能是二三產(chǎn)業(yè)比重提升加劇了環(huán)境污染, 顯示出中國二三產(chǎn)業(yè)清潔生產(chǎn)能力仍然有待提升。

建成區(qū)城鎮(zhèn)面積在5%的水平顯著負(fù)向影響農(nóng)業(yè)生態(tài)效率, 人均建成區(qū)面積在1%水平顯著負(fù)向影響農(nóng)業(yè)生態(tài)效率, 即隨著建成區(qū)城鎮(zhèn)面積和人均建成區(qū)面積擴(kuò)大, 農(nóng)業(yè)生態(tài)效率將下降。其原因可能是城市用地?cái)U(kuò)大將侵占農(nóng)業(yè)用地, 而隨著農(nóng)業(yè)用地減少, 農(nóng)業(yè)生產(chǎn)受到影響, 一方面農(nóng)業(yè)用地不足直接造成農(nóng)業(yè)產(chǎn)出下降, 另一方面也可能造成更多化學(xué)制品使用, 從而加劇農(nóng)業(yè)生態(tài)環(huán)境污染。

4 結(jié)論與建議

4.1 研究結(jié)論

本文在剔除外生環(huán)境影響和隨機(jī)干擾基礎(chǔ)上從多個(gè)角度探討城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響, 主要研究結(jié)論如下:

1)外生環(huán)境影響和隨機(jī)干擾對于農(nóng)業(yè)生態(tài)效率測度會(huì)產(chǎn)生顯著影響, 并進(jìn)一步影響城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的作用。第1階段農(nóng)業(yè)生態(tài)效率與第3階段農(nóng)業(yè)生態(tài)效率存在較大差距, 剔除外生環(huán)境影響和隨機(jī)干擾后糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率出現(xiàn)明顯變化。此外, 在與現(xiàn)有研究對比中, 部分城鎮(zhèn)化指標(biāo)對剔除外生環(huán)境和隨機(jī)干擾后農(nóng)業(yè)生態(tài)效率影響方向與現(xiàn)有研究相比出現(xiàn)變化。

2)中國糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率持續(xù)增長, 但是地區(qū)之間存在差距。根據(jù)剔除了外生環(huán)境影響和隨機(jī)干擾的第3階段農(nóng)業(yè)生態(tài)效率可以看出, 中國糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率呈現(xiàn)持續(xù)增長趨勢, 但是不同地區(qū)間農(nóng)業(yè)生態(tài)效率存在差距, 部分糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率達(dá)到效率前沿, 但仍有部分地區(qū)農(nóng)業(yè)生態(tài)效率處于較低水平。

表6 城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率影響的估計(jì)結(jié)果

*、**、***表示在10%、5%、1%水平顯著。*, **, *** represent significant effects at levels of 10%, 5%, 1%, respectively.

3)城鎮(zhèn)化整體上能夠促進(jìn)農(nóng)業(yè)生態(tài)效率的提升, 具體城鎮(zhèn)化指標(biāo)對農(nóng)業(yè)生態(tài)效率的影響復(fù)雜, 因此探討城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響不能僅僅從單一方面進(jìn)行考慮, 因當(dāng)多角度進(jìn)行分析。根據(jù)分析可知, 城鎮(zhèn)化綜合指數(shù)顯著正向影響農(nóng)業(yè)生態(tài)效率; 從各指標(biāo)來看, 人口城鎮(zhèn)化率負(fù)向影響農(nóng)業(yè)生態(tài)效率, 但不顯著; 城鎮(zhèn)居民人均可支配收入、城鎮(zhèn)經(jīng)濟(jì)密度正向影響農(nóng)業(yè)生態(tài)效率; 二三產(chǎn)業(yè)產(chǎn)值占GDP比重、建成區(qū)面積比重、人均建成區(qū)面積負(fù)向影響農(nóng)業(yè)生態(tài)效率。

4.2 政策建議

根據(jù)研究可知城鎮(zhèn)化整體上能夠提升農(nóng)業(yè)生態(tài)效率, 為此應(yīng)當(dāng)進(jìn)一步推進(jìn)城鎮(zhèn)化, 并不斷提升城鎮(zhèn)化質(zhì)量, 下面從人口城鎮(zhèn)化、經(jīng)濟(jì)城鎮(zhèn)化、土地城鎮(zhèn)化3個(gè)方面提出建議。

1)人口城鎮(zhèn)化穩(wěn)定推進(jìn), 切實(shí)提升居民收入水平。人口城鎮(zhèn)化是城鎮(zhèn)化基礎(chǔ), 雖然人口城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率會(huì)產(chǎn)生一定負(fù)向影響, 但是過剩的農(nóng)業(yè)人口進(jìn)入城市可以減少人類活動(dòng)對農(nóng)業(yè)生態(tài)環(huán)境的直接污染, 可為城市發(fā)展提供動(dòng)力, 而城市發(fā)展又能為農(nóng)業(yè)生態(tài)環(huán)境改善提供更多資金。為此, 需要逐步引導(dǎo)農(nóng)村過剩勞動(dòng)力向城市轉(zhuǎn)移, 在此過程中應(yīng)當(dāng)重視農(nóng)業(yè)勞動(dòng)力流失對農(nóng)業(yè)發(fā)展的影響, 積極推動(dòng)先進(jìn)耕種技術(shù)推廣及綠色肥料應(yīng)用, 降低農(nóng)業(yè)勞動(dòng)力減少對農(nóng)業(yè)生產(chǎn)的負(fù)面影響。此外, 從研究中可以看出居民可支配收入的提升能夠有效提高農(nóng)業(yè)生態(tài)效率, 為此應(yīng)通過穩(wěn)定就業(yè)、降低稅負(fù)、拓寬居民投資渠道、加強(qiáng)社會(huì)保障等途徑提升居民可支配收入, 提高居民生活水平。

2)經(jīng)濟(jì)城鎮(zhèn)化清潔發(fā)展, 有效降低產(chǎn)業(yè)聚集污染。從研究中可以看出城鎮(zhèn)經(jīng)濟(jì)密度的提高能夠促進(jìn)農(nóng)業(yè)生態(tài)效率提升, 為此應(yīng)當(dāng)在提升城鎮(zhèn)產(chǎn)業(yè)清潔生產(chǎn)水平和加強(qiáng)產(chǎn)區(qū)環(huán)保體系建設(shè)基礎(chǔ)上促進(jìn)二三產(chǎn)業(yè)向城鎮(zhèn)地區(qū)集聚, 減少二三產(chǎn)業(yè)生產(chǎn)對農(nóng)業(yè)生態(tài)環(huán)境的影響。促進(jìn)二三產(chǎn)業(yè)向城鎮(zhèn)集聚應(yīng)當(dāng)充分考慮產(chǎn)業(yè)園區(qū)環(huán)境承載力等因素, 避免產(chǎn)業(yè)于單一區(qū)域的過度聚集, 綜合考量產(chǎn)業(yè)特點(diǎn)和城鎮(zhèn)空間分布, 通過提升產(chǎn)業(yè)園區(qū)清潔生產(chǎn)技術(shù)、強(qiáng)化廢棄物回收利用、優(yōu)化產(chǎn)業(yè)園區(qū)位置選擇等方式減輕經(jīng)濟(jì)城鎮(zhèn)化帶來的城市環(huán)境問題。

3)土地城鎮(zhèn)化適度實(shí)現(xiàn), 推動(dòng)大小城市與小城鎮(zhèn)協(xié)調(diào)發(fā)展。從研究中可以看出土地城鎮(zhèn)化會(huì)制約農(nóng)業(yè)生態(tài)效率提升, 為此應(yīng)當(dāng)推動(dòng)土地城鎮(zhèn)化適度發(fā)展, 在充分考慮地區(qū)發(fā)展實(shí)際和產(chǎn)業(yè)發(fā)展特征基礎(chǔ)上, 因地制宜推進(jìn)衛(wèi)星城、小城鎮(zhèn)建設(shè), 避免對城市規(guī)模的片面追求, 避免形成畸大畸小城鎮(zhèn)格局, 以實(shí)現(xiàn)地區(qū)城鎮(zhèn)協(xié)調(diào)發(fā)展, 減少土地城鎮(zhèn)化過程對農(nóng)業(yè)用地的過度擠壓。

[1] 項(xiàng)繼權(quán), 周長友. “新三農(nóng)”問題的演變與政策選擇[J]. 中國農(nóng)村經(jīng)濟(jì), 2017, (10): 13–25 XIANG J Q, ZHOU C Y. Evolution and policy choices of new three rural issues[J]. China Rural Economy, 2017, (10): 13–25

[2] 張貴友. 鄉(xiāng)村振興背景下“空心村”治理對策研究——基于安徽省的調(diào)查[J]. 江淮論壇, 2019, (5): 37–42 ZHANG G Y. Study on the countermeasures of “Hollow Village” governance under the background of vitalization of rural areas — Investigation based on Anhui Province[J]. Jianghuai Tribune, 2019, (5): 37–42

[3] WESTLUND H. Urban futures in planning, policy and regional science: Are we entering a post-urban world?[J]. Built Environment, 2014, 40(4): 447–457

[4] DENG X Z, HUANG J K, ROZELLE S, et al. Impact of urbanization on cultivated land changes in China[J]. Land Use Policy, 2015, 45: 1–7

[5] ZHANG R J, JIANG G H, ZHANG Q. Does urbanization always lead to rural hollowing? Assessing the spatio-temporal variations in this relationship at the county level in China 2000–2015[J]. Journal of Cleaner Production, 2019, 220: 9–22

[6] ANDERSEN H T. Planning on the edge: The context for planning at the rural-urban fringe by Nick Gallent; Johan Andersson; Marco Bianconi[J]. The Town Planning Review, 2007, 78(4): 539–540

[7] ZHOU Y, LIU Y S, WU W X. Disaster risk: Strengthen China’s flood control[J]. Nature, 2016, 536(7617): 396

[8] 吳小慶, 王亞平, 何麗梅, 等. 基于AHP和DEA模型的農(nóng)業(yè)生態(tài)效率評價(jià)——以無錫市為例[J]. 長江流域資源與環(huán)境, 2012, 21(6): 714—719 WU X Q, WANG Y P, HE L M, et al. Agricultural eco-efficiency evaluation based on AHP and DEA model — A case of Wuxi City[J]. Resources and Environment in the Yangtza Basin, 2012, 21(6): 714—719

[9] 張楊, 陳娟娟. 農(nóng)業(yè)生態(tài)效率的國際比較及中國的定位研究[J]. 中國軟科學(xué), 2019, (10): 165–172 ZHANG Y, CHEN J J. International comparison of agroecological efficiency and China’s positioning study[J]. China Soft Science, 2019, (10): 165–172

[10] 王寶義, 張衛(wèi)國. 中國農(nóng)業(yè)生態(tài)效率測度及時(shí)空差異研究[J]. 中國人口?資源與環(huán)境, 2016, 26(6): 11–19 WANG B Y, ZHANG W G. A research of agricultural eco-efficiency measure in China and space-time differences[J]. China Population, Resources and Environment, 2016, 26(6): 11–19

[11] 胡平波, 鐘漪萍. 政府支持下的農(nóng)旅融合促進(jìn)農(nóng)業(yè)生態(tài)效率提升機(jī)理與實(shí)證分析——以全國休閑農(nóng)業(yè)與鄉(xiāng)村旅游示范縣為例[J]. 中國農(nóng)村經(jīng)濟(jì), 2019, (12): 85–104 HU P B, ZHANG Y P. The mechanism of improving agricultural eco-efficiency by the integration of agriculture and tourism supported by the government: Taking the national leisure agriculture and rural tourism demonstration counties as an example[J]. China’s Rural Economy, 2019, (12): 85–104

[12] 侯孟陽, 姚順波. 中國農(nóng)村勞動(dòng)力轉(zhuǎn)移對農(nóng)業(yè)生態(tài)效率影響的空間溢出效應(yīng)與門檻特征[J]. 資源科學(xué), 2018, 40(12): 2475–2486 HOU M Y, YAO S B. spatial spillover effects and threshold characteristics of rural labor transfer on agricultural eco-efficiency in China[J]. Resources Science, 2018, 40(12): 2475–2486

[13] 王寶義, 張衛(wèi)國. 中國農(nóng)業(yè)生態(tài)效率的省際差異和影響因素——基于1996~2015年31個(gè)省份的面板數(shù)據(jù)分析[J]. 中國農(nóng)村經(jīng)濟(jì), 2018, (1): 46–62 WANG B Y, ZHANG W G. Cross-provincial differences in determinants of agricultural eco-efficiency in China: An analysis based on panel data from 31 provinces in 1996–2015[J]. Chinese Rural Economy, 2018, (1): 46–62

[14] 曹俊文, 曾康. 低碳視角下長江經(jīng)濟(jì)帶農(nóng)業(yè)生態(tài)效率及影響因素研究[J]. 生態(tài)經(jīng)濟(jì), 2019, 35(8): 115–119 CAO J W, ZENG K. Study on agricultural eco-efficiency and its influencing factors in the Yangtze River economic belt from the perspective of low carbon[J]. Ecological Economy, 2019, 35(8): 115–119

[15] 鄭麗楠, 洪名勇. 中國農(nóng)業(yè)生態(tài)效率的時(shí)空特征及驅(qū)動(dòng)因素[J]. 江西財(cái)經(jīng)大學(xué)學(xué)報(bào), 2019, (5): 46–56 ZHENG L N, HONG M Y. Spatial and temporal characteristics and driving factors of agricultural ecological efficiency in China[J]. Journal of Jiangxi University of Finance and Economics, 2019, (5): 46–56

[16] 薛俊菲, 陳雯, 張蕾. 中國市域綜合城市化水平測度與空間格局研究[J]. 經(jīng)濟(jì)地理, 2010, 30(12): 2005–2011 XUE J F, CHEN W, ZHANG L. Comprehensive measurement and spatial pattern of urbanization level of China’s cities[J]. Economic Geography, 2010, 30(12): 2005–2011

[17] 范虹玨, 劉祖云. 中國城鎮(zhèn)化空間發(fā)展態(tài)勢研究——基于人口、土地、經(jīng)濟(jì)城鎮(zhèn)化協(xié)調(diào)發(fā)展的視角[J]. 內(nèi)蒙古社會(huì)科學(xué): 漢文版, 2014, 35(1): 95–100 FAN H Y, LIU Z Y. A study on the spatial development of urbanization in China-based on the perspective of coordinated development of population, land and economic urbanization[J]. Inner Mongolia Social Sciences: Han Version, 2014, 35(1): 95–100

[18] 劉旭曄. 城鎮(zhèn)化對耕地可持續(xù)利用的影響研究——基于湖北省縣級面板數(shù)據(jù)分析[J]. 中南財(cái)經(jīng)政法大學(xué)學(xué)報(bào), 2015, (2): 20–27LIU X Y. Study of the influence of urbanization on sustainable utilization of cultivated land: Based on the analysis of the county level data in Hubei Province[J]. Journal of Zhongnan University of Economics and Law, 2015, (2): 20–27

[19] 賀三維, 邵璽. 京津冀地區(qū)人口-土地-經(jīng)濟(jì)城鎮(zhèn)化空間集聚及耦合協(xié)調(diào)發(fā)展研究[J]. 經(jīng)濟(jì)地理, 2018, (1): 95–102 HE S W, SHAO X. Spatial clustering and coupling coordination of population-land-economic urbanization in Beijing-Tianjin-Hebei Region[J]. Economic Geography, 2018, (1): 95–102

[20] 李新剛, 孫鈺. 人口、土地、經(jīng)濟(jì)城鎮(zhèn)化動(dòng)態(tài)影響效應(yīng)研究[J]. 經(jīng)濟(jì)問題, 2018, (3): 106–111 LI X G, SUN Y. The dynamic effect of population, land, economic urbanization[J]. On Economic Problems, 2018, (3): 106–111

[21] 盧新海, 柯楠, 匡兵, 等. 中部地區(qū)土地城鎮(zhèn)化水平差異的時(shí)空特征及影響因素[J]. 經(jīng)濟(jì)地理, 2019, 39(4): 192–198 LU X H, KE N, KUANG B, et al. Spatial-temporal features and influencing factors of difference in land urbanization level of central China[J]. Economic Geography, 2019, 39(4): 192–198

[22] 蔡之兵, 鄧仲良. 中國空間城鎮(zhèn)化水平測度及其影響因素分析——基于省級地理加權(quán)回歸模型的證據(jù)[J]. 鄭州大學(xué)學(xué)報(bào): 哲學(xué)社會(huì)科學(xué)版, 2020, 53(266): 48–54 CAI Z B, DENG Z L. An analysis of China’s spatial urbanization level and its influencing factors: Based on geographically weighted regression model of Chinese provinces[J]. Journal of Zhengzhou University: Philosophy and Social Sciences Edition, 2020, 53(266): 48–54

[23] STIGSON B. Eco-efficiency: Creating More Value with Less Impact[R]. North Yorkshire, UK: WBCSD, 2000: 5–36

[24] FRIED H O, LOVELL C A K, SCHMIDT S S, et al. Accounting for environmental effects and statistical noise in data envelopment analysis[J]. Journal of Productivity Analysis, 2002, 17(1): 157–174

[25] 張艦, 亞伯拉罕·艾賓斯坦, 瑪格麗特·麥克米倫, 等. 農(nóng)村勞動(dòng)力轉(zhuǎn)移、化肥過度使用與環(huán)境污染[J]. 經(jīng)濟(jì)社會(huì)體制比較, 2017, (3): 149–160 ZHANG J, EBENSTEIN A, MCMILLAN M, et al. Migration, excessive fertilizer use and environmental consequences[J]. Comparative Economic and Social Systems, 2017, (3): 149–160

[26] 何興邦, 周葵. 城鎮(zhèn)化對進(jìn)城農(nóng)民環(huán)保行為影響研究[J]. 中國人口?資源與環(huán)境, 2016, 26(8): 70–78 HE X B, ZHOU K. Effects of urbanization on the environmental behaviors of peasants in cities[J]. China Population, Resources and Environment, 2016, 26(8): 70–78

[27] 劉洪波. 以綠色消費(fèi)助推社會(huì)主義新農(nóng)村建設(shè)[J]. 行政論壇, 2009, 16(5): 87–89 LIU H B. Promote the construction of new socialist countryside with green consumption[J]. Administrative Tribune, 2009, 16(5): 87–89

[28] 韓紀(jì)江, 孔祥智. 城鎮(zhèn)化進(jìn)程對農(nóng)村經(jīng)濟(jì)的負(fù)面效應(yīng)淺議[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2001, (7): 26–30 HAN J J, KONG X Z. Discussions on the negative impacts of urbanization on rural economy[J]. Issues in Agricultural Economy, 2001, (7): 26–30

[29] 課題組. 農(nóng)地城市流轉(zhuǎn)對周邊區(qū)域農(nóng)業(yè)生產(chǎn)的負(fù)外溢效應(yīng)研究[J]. 農(nóng)業(yè)經(jīng)濟(jì)問題, 2010, (12): 98–101 Discussion Group. Study on the negative spillover effects on local agricultural production arising from the rural-urban land conversion process[J]. Issues in Agricultural Economy, 2010, (12): 98–101

[30] 高楊, 牛子恒. 農(nóng)業(yè)信息化、空間溢出效應(yīng)與農(nóng)業(yè)綠色全要素生產(chǎn)率——基于SBM-ML指數(shù)法和空間杜賓模型[J]. 統(tǒng)計(jì)與信息論壇, 2018, 33(10): 66–75 GAO Y, NIU Z H. Agricultural informatization, spatial spillover effect and agricultural green total factor productivity: Based on the method of SBM-ML index and spatial Durbin model[J]. Statistics & Information Forum, 2018, 33(10): 66–75

[31] 劉旭曄, 張開華, 王穎. 城鎮(zhèn)化快速發(fā)展背景下耕地可持續(xù)利用動(dòng)態(tài)評價(jià)——基于湖北省2003~2012年縣級層面數(shù)據(jù)分析[J]. 農(nóng)村經(jīng)濟(jì)與科技, 2015, 26(5): 15–19 LIU X Y, ZHANG K H, WANG Y. Dynamic evaluation of sustainable use of cultivated land under the background of rapid development of urbanization — Based on the county level data analysis of Hubei Province from 2003 to 2012[J]. Rural Economy and Science-Technology, 2015, 26(5): 15–19

[32] TONE K. A slacks-based measure of efficiency in data envelopment analysis[J]. European Journal of Operational Research, 2001, 130(3): 498–509

[33] TONE K. Dealing with Undesirable Outputs in DEA: A Slacks-based Measure (SBM) Approach[M]. Toronto: North American Productivity Workshop, 2004.

[34] 羅登躍. 三階段DEA模型管理無效率估計(jì)注記[J]. 統(tǒng)計(jì)研究, 2012, 29(4): 104–107 LUO D Y. A note on estimating managerial inefficiency of three-stage DEA model[J]. Statistical Research, 2012, 29(4): 104–107

[35] 陳巍巍, 張雷, 馬鐵虎, 等. 關(guān)于三階段DEA模型的幾點(diǎn)研究[J]. 系統(tǒng)工程, 2014, 32(9): 144–149 CHEN W W, ZHANG L, MA T H, et al. Research on three-stage DEA model[J]. Systems Engineering, 2014, 32(9): 144–149

[36] 李寧, 何文劍, 仇童偉, 等. 農(nóng)地產(chǎn)權(quán)結(jié)構(gòu)、生產(chǎn)要素效率與農(nóng)業(yè)績效[J]. 管理世界, 2017, (3): 44–62 LI N, HE W J, QIU T W, et al. Farmland property rights structure, production factor efficiency and agricultural performance[J]. Management World, 2017, (3): 44–62

[37] 黃瑪蘭, 李曉云, 游良志. 農(nóng)業(yè)機(jī)械與農(nóng)業(yè)勞動(dòng)力投入對糧食產(chǎn)出的影響及其替代彈性[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào): 社會(huì)科學(xué)版, 2018, (2): 37–45 HUANG M L, LI X Y, YOU L Z. The influence of agricultural machinery and agricultural labor input on grain output and its substitution elasticity[J]. Journal of Huazhong Agricultural University: Social Sciences Edition, 2018, (2): 37–45

[38] 賴斯蕓, 杜鵬飛, 陳吉寧. 基于單元分析的非點(diǎn)源污染調(diào)查評估方法[J]. 清華大學(xué)學(xué)報(bào): 自然科學(xué)版, 2004, 44(9): 1184–1187 LAI S Y, DU P F, CHEN J N. Evaluation of non-point source pollution based on unit analysis[J]. Journal of Tsinghua University: Science and Technology, 2004, 44(9): 1184–1187

[39] BATTESE G E, COELLI T J. Frontier production functions, technical efficiency and panel data: With application to paddy farmers in India[J]. The Journal of Productivity Analysis, 1992, 3(1/2): 153–169

[40] COELLI T J. RAO D S P. Total factor productivity growth in agriculture: A malmquist index analysis of 93 countries, 1980–2000[J]. Agricultural Economics, 2005, 32(S1): 115–134

Study on the impact of China’s urbanization on agricultural ecological efficiency: Based on panel data of 13 major grain-producing regions in China from 2009 to 2018*

SHANG Jie, JI Xueqiang, CHEN Ximing

(School of Economics and Management, Northeast Forestry University, Harbin 150040, China)

Urbanization is an important factor to impact agricultural ecological efficiency. However, most of the existing studies have discussed the impact of urbanization on agricultural ecological efficiency from a single perspective. Further, these studies are often affected by the impact of exogenous environment or other random factors during the measurement of agricultural ecological efficiency. To analyze the impact of urbanization on agricultural ecological efficiency more comprehensively and accurately, the present study used the three-stage SBM-DEA model, which provides simultaneously evaluation of both the undesirable and desirable outputs, entropy value method, and Tobit model, based on panel data obtained from 13 major grain-producing regions in China from 2009 to 2018. Various models and methods were used to study the impact of urbanization on agricultural ecological efficiency from multiple perspectives based on theoretical analysis. Our results showed that: 1) the impact of exogenous environment and random interference significantly impacted the agricultural ecological efficiency. After excluding these impacts, the agricultural ecological efficiency of the 13 major grain-producing regions in China changed significantly. Additionally, unlike the existing studies, the impact of urbanization on agricultural ecological efficiency too changed, after the effects of exogenous environment and random interference were removed. 2) The agricultural ecological efficiency of China’s 13 major grain-producing regions showed a continuous growth trend from 2009 to 2018. The average annual efficiency, after excluding the impact of exogenous environment and random interference, increased from 0.53 in 2009 to 0.80 in 2018. However, there was a big gap in the agricultural ecological efficiency among different regions. For example, in 2018, the agricultural ecological efficiency of Hebei and other regions reached the frontier (1.00), but that of Inner Mongolia was only 0.45. 3) Urbanization as a whole improved the agricultural ecological efficiency. The comprehensive index of urbanization significantly affected agricultural ecological efficiency (at a level of 1%, impact coefficient: 0.60). 4) Various indicators of urbanization differently impacted the agricultural ecological efficiency. The population urbanization rate showed an insignificant negative impact on the agricultural ecological efficiency. The per capita disposable income of urban residents and urban economic density significantly positively affected the agricultural ecological efficiency. The output value of the secondary and tertiary industries, proportion of built-up area and per capita built-up area had significantly negative effect on the agricultural ecological efficiency. To further elucidate the role of urbanization in promoting agricultural eco-efficiency, this article proposed suggestions to modify the relevant policies, based on the research results from three aspects of population, economy, and land urbanization.

Urbanization; Agro-ecological efficiency; Main grain production areas; DEA; Tobit

F323.22

10.13930/j.cnki.cjea.200135

尚杰, 吉雪強(qiáng), 陳璽名. 中國城鎮(zhèn)化對農(nóng)業(yè)生態(tài)效率的影響——基于中國13個(gè)糧食主產(chǎn)區(qū)2009—2018年面板數(shù)據(jù)[J]. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文), 2020, 28(8): 1265-1276

SHANG J, JI X Q, CHEN X M. Study on the impact of China’s urbanization on agricultural ecological efficiency: Based on panel data of 13 major grain-producing regions in China from 2009 to 2018[J]. Chinese Journal of Eco-Agriculture, 2020, 28(8): 1265-1276

* 國家自然科學(xué)基金項(xiàng)目(71573036)資助

尚杰, 主要從事農(nóng)業(yè)經(jīng)濟(jì)、環(huán)境科學(xué)與資源利用研究。E-mail: 13576274061@163.com

2020-02-27

2020-05-26

* This work was supported by the National Natural Science Foundation of China (71573036).

, SHANG Jie, E-mail: 13576274061@163.com

Feb. 27, 2020;

May 26, 2020