農(nóng)業(yè)環(huán)境質(zhì)量對(duì)綠色農(nóng)業(yè)發(fā)展影響的空間效應(yīng)分析

林雨仙 鄢 貞 * 胡偉斌 白 麗

（1.浙江大學(xué)中國(guó)農(nóng)村發(fā)展研究院 浙江 杭州 郵編：310058；2.浙江大學(xué)公共管理學(xué)院農(nóng)業(yè)與農(nóng)村發(fā)展智能計(jì)算機(jī)實(shí)驗(yàn)室 浙江 杭州 郵編：310058；3.吉林大學(xué)生物與農(nóng)業(yè)工程學(xué)院 吉林 長(zhǎng)春 郵編：130022）

一、引 言

隨著全球?qū)ι鷳B(tài)環(huán)境與可持續(xù)發(fā)展的日益重視，農(nóng)業(yè)生產(chǎn)中化肥農(nóng)藥施用帶來(lái)農(nóng)業(yè)面源污染，造成生態(tài)環(huán)境風(fēng)險(xiǎn)已成為各界熱議的話題。降低對(duì)化肥農(nóng)藥依賴程度，減肥增效、高產(chǎn)高效成為農(nóng)業(yè)可持續(xù)發(fā)展與生態(tài)安全重要路徑之一[1]。以歐盟、美國(guó)為代表的發(fā)達(dá)國(guó)家和地區(qū)肥藥施用強(qiáng)度變化表現(xiàn)出隨著經(jīng)濟(jì)發(fā)展呈現(xiàn)快速增加，而后逐漸降低的趨勢(shì)[2]。為保障糧食安全，發(fā)展中國(guó)家農(nóng)業(yè)生產(chǎn)中化肥農(nóng)藥施用強(qiáng)度，仍處于發(fā)達(dá)國(guó)家經(jīng)歷的前半段“快速增長(zhǎng)期”。以中國(guó)為例，自改革開放以來(lái)，社會(huì)經(jīng)濟(jì)迅速發(fā)展與人口增加，同時(shí)帶來(lái)糧食安全保障壓力，也加重了農(nóng)業(yè)生產(chǎn)對(duì)生態(tài)環(huán)境的破壞與污染，造成嚴(yán)峻的農(nóng)產(chǎn)品質(zhì)量安全問(wèn)題。

我國(guó)出臺(tái)了一系列綠色農(nóng)業(yè)發(fā)展的指導(dǎo)意見和政策，突出綠色農(nóng)業(yè)發(fā)展在解決農(nóng)業(yè)源頭生態(tài)環(huán)境污染和促進(jìn)肥藥施用減量的重要性，逐漸形成并建立健全我國(guó)“三品一標(biāo)”綠色農(nóng)業(yè)產(chǎn)業(yè)發(fā)展體系。經(jīng)過(guò)30 多年發(fā)展，綠色食品標(biāo)志許可體系推廣應(yīng)用，對(duì)于農(nóng)產(chǎn)品質(zhì)量安全的提升、農(nóng)業(yè)質(zhì)量與效益改善以及農(nóng)民增收都起到了重要作用[3]。從根本上來(lái)看，生態(tài)安全的產(chǎn)地環(huán)境是綠色農(nóng)產(chǎn)品生產(chǎn)和綠色農(nóng)業(yè)可持續(xù)發(fā)展的基礎(chǔ)性資源[4]。由于各地原有農(nóng)業(yè)環(huán)境質(zhì)量狀況、農(nóng)業(yè)發(fā)展水平以及經(jīng)濟(jì)技術(shù)條件存在差異，綠色農(nóng)業(yè)仍存在區(qū)域間發(fā)展不平衡問(wèn)題?，F(xiàn)有研究尚未從微觀生產(chǎn)主體角度探究農(nóng)業(yè)環(huán)境質(zhì)量對(duì)可持續(xù)生產(chǎn)方式的空間作用機(jī)理。因而，農(nóng)業(yè)面源污染產(chǎn)生的農(nóng)業(yè)環(huán)境質(zhì)量約束，是否促進(jìn)農(nóng)業(yè)經(jīng)營(yíng)主體從事農(nóng)資減量、環(huán)境友好、品質(zhì)優(yōu)良綠色食品的標(biāo)準(zhǔn)化生產(chǎn)管理，進(jìn)而推進(jìn)農(nóng)業(yè)產(chǎn)業(yè)轉(zhuǎn)型升級(jí)，且這一轉(zhuǎn)變是某一區(qū)域特殊性還是會(huì)存在波特假說(shuō)下技術(shù)創(chuàng)新帶來(lái)的溢出效應(yīng)，還有待于進(jìn)一步研究和實(shí)證檢驗(yàn)。如何平衡生態(tài)環(huán)境保護(hù)和綠色高質(zhì)量現(xiàn)代農(nóng)業(yè)發(fā)展，以綠色食品為切入點(diǎn)的空間差異性和農(nóng)業(yè)環(huán)境質(zhì)量約束研究尤其重要。

本文選取中國(guó)2011—2020 年的31 個(gè)省份（港澳臺(tái)除外）267 個(gè)城市新增綠色食品數(shù)作為研究對(duì)象，對(duì)各省域生態(tài)環(huán)境質(zhì)量水平與果蔬類綠色食品的時(shí)空分布進(jìn)行識(shí)別，并構(gòu)建空間杜賓模型檢驗(yàn)生態(tài)環(huán)境水平對(duì)綠色食品影響的空間效應(yīng)。本研究可能存在以下兩點(diǎn)創(chuàng)新：其一，首次對(duì)獲得綠色食品標(biāo)志許可的歷年微觀數(shù)據(jù)進(jìn)行采集，嘗試分析全國(guó)大尺度下綠色食品數(shù)量水平的空間集聚情況；其二，對(duì)近年來(lái)中國(guó)農(nóng)業(yè)綠色發(fā)展過(guò)程中表現(xiàn)出的顯著的空間溢出性和示范效應(yīng)展開探討。因而，本文通過(guò)對(duì)中國(guó)各省域綠色食品分布的空間差異性分析、空間效應(yīng)檢驗(yàn)以及其與農(nóng)業(yè)環(huán)境質(zhì)量的響應(yīng)關(guān)系分析，從實(shí)踐上驗(yàn)證 “兩山”理論的科學(xué)指導(dǎo)性，對(duì)中國(guó)農(nóng)業(yè)高質(zhì)量發(fā)展的理論發(fā)展具有一定借鑒意義。

二、綠色農(nóng)業(yè)發(fā)展的理論背景

生態(tài)環(huán)境具有準(zhǔn)公共品和信任品屬性[5,6]，農(nóng)業(yè)生產(chǎn)主體一般不會(huì)積極主動(dòng)進(jìn)行環(huán)境友好型生產(chǎn)。波特假說(shuō)認(rèn)為環(huán)境保護(hù)政策主要通過(guò)促進(jìn)企業(yè)技術(shù)創(chuàng)新或引進(jìn)創(chuàng)新技術(shù)來(lái)影響經(jīng)濟(jì)，短期來(lái)看成本增加，但長(zhǎng)期而言可以提高生產(chǎn)效率，增強(qiáng)競(jìng)爭(zhēng)力，在實(shí)現(xiàn)經(jīng)濟(jì)增長(zhǎng)的同時(shí)保護(hù)生態(tài)環(huán)境[7]。在綠色農(nóng)業(yè)發(fā)展進(jìn)程中，農(nóng)業(yè)面源污染治理帶來(lái)的農(nóng)業(yè)環(huán)境質(zhì)量提升對(duì)高品質(zhì)農(nóng)業(yè)生產(chǎn)轉(zhuǎn)型在長(zhǎng)短期可能存在截然不同的影響。從短期看，政策規(guī)制下，生產(chǎn)者從農(nóng)資密集型農(nóng)業(yè)生產(chǎn)方式向勞動(dòng)力或技術(shù)密集型生產(chǎn)方式轉(zhuǎn)換，可能存在投入成本增加，并且市場(chǎng)“優(yōu)質(zhì)優(yōu)價(jià)”未能有效實(shí)現(xiàn)溢價(jià)收益，難以抵消成本[8]；但從長(zhǎng)期看，生產(chǎn)環(huán)境持續(xù)改善后生產(chǎn)技術(shù)的規(guī)模效應(yīng)發(fā)揮作用，以及生態(tài)農(nóng)產(chǎn)品市場(chǎng)優(yōu)質(zhì)優(yōu)價(jià)激勵(lì)，可以促進(jìn)高質(zhì)量農(nóng)業(yè)有效、長(zhǎng)久發(fā)展。因而我國(guó)一系列環(huán)境政策如“耕地分類管理”“地下水污染防治”“農(nóng)村人居環(huán)境整治”“肥藥雙減”等源頭治理政策，從長(zhǎng)遠(yuǎn)看可以促進(jìn)以綠色食品為代表的高質(zhì)量農(nóng)業(yè)生產(chǎn)方式快速發(fā)展。

綠色農(nóng)業(yè)生產(chǎn)技術(shù)改進(jìn)不僅體現(xiàn)在保障農(nóng)產(chǎn)品質(zhì)量安全方面，同時(shí)通過(guò)規(guī)范化管理對(duì)生態(tài)環(huán)境存在嚴(yán)格要求，因而農(nóng)業(yè)中環(huán)境友好型生產(chǎn)技術(shù)同工業(yè)創(chuàng)新技術(shù)有著相似溢出效應(yīng)[9,10]。根據(jù)地理學(xué)第一定律，由于在地理形態(tài)上，相鄰地區(qū)在自然資源條件、經(jīng)濟(jì)發(fā)展水平和技術(shù)水平等方面有很強(qiáng)的相似性，一個(gè)地區(qū)環(huán)境友好型技術(shù)對(duì)鄰近地區(qū)存在著示范效應(yīng)[11]，并且，示范效應(yīng)可能會(huì)引起本地區(qū)當(dāng)期農(nóng)業(yè)生產(chǎn)調(diào)整，促進(jìn)本區(qū)改善農(nóng)業(yè)生產(chǎn)環(huán)境，對(duì)農(nóng)產(chǎn)品生產(chǎn)過(guò)程進(jìn)行標(biāo)準(zhǔn)化管理，形成同群效應(yīng)或?qū)W習(xí)帶來(lái)的技術(shù)溢出效應(yīng)?，F(xiàn)有文獻(xiàn)認(rèn)為綠色生產(chǎn)行為和生態(tài)環(huán)境具有典型的正外部性[12]，例如退耕還林等政策在非市場(chǎng)條件下具有外溢效應(yīng)[13]。而在多數(shù)農(nóng)戶經(jīng)濟(jì)行為的研究中，沒(méi)有考慮到個(gè)體生產(chǎn)行為的外部性[14]。因而在考慮市場(chǎng)條件下的農(nóng)戶或者農(nóng)業(yè)經(jīng)營(yíng)主體生態(tài)可持續(xù)的生產(chǎn)行為時(shí)，應(yīng)充分討論區(qū)域間外溢效應(yīng)產(chǎn)生的正外部性帶來(lái)生態(tài)價(jià)值增加，正確估計(jì)生態(tài)、綠色可持續(xù)農(nóng)業(yè)的生產(chǎn)管理技術(shù)帶來(lái)潛在的生態(tài)價(jià)值轉(zhuǎn)換作用。因此本文探究生態(tài)環(huán)境約束對(duì)綠色食品的空間溢出效應(yīng)是可行且具有意義的。

三、綠色農(nóng)業(yè)發(fā)展與農(nóng)業(yè)環(huán)境質(zhì)量的時(shí)空演化特征

（一）數(shù)據(jù)來(lái)源

中國(guó)綠色食品，是指產(chǎn)自優(yōu)良生態(tài)環(huán)境、按照綠色食品標(biāo)準(zhǔn)生產(chǎn)、實(shí)行全程質(zhì)量控制并獲得綠色食品標(biāo)志使用權(quán)的安全、優(yōu)質(zhì)食用農(nóng)產(chǎn)品及相關(guān)產(chǎn)品。

當(dāng)前水果蔬菜為代表的新鮮農(nóng)產(chǎn)品產(chǎn)量占比逐年增加，市場(chǎng)發(fā)展蓬勃；同時(shí)，受到農(nóng)業(yè)面源污染的影響，果蔬產(chǎn)品農(nóng)殘超標(biāo)等農(nóng)產(chǎn)品質(zhì)量安全問(wèn)題較為嚴(yán)重[15]，因而本文以獲得綠色食品標(biāo)志使用權(quán)的果蔬類產(chǎn)品為研究對(duì)象。

本文選擇2011—2020 年全國(guó)267 個(gè)地級(jí)市（不包括港澳臺(tái)地區(qū)）的果蔬類綠色食品作為研究對(duì)象，綠色食品的相關(guān)數(shù)據(jù)采集來(lái)自中國(guó)綠色食品網(wǎng)站，并與企查查、天眼查網(wǎng)站信息匹配，利用文本語(yǔ)義識(shí)別獲得綠色食品標(biāo)志使用權(quán)主體的基本信息，共獲得75634 條數(shù)據(jù)，選取新鮮果蔬類34 452 條作為本文研究對(duì)象。由于綠色食品標(biāo)志使用權(quán)有效期為3 年，本文的有效綠色食品數(shù)量年度統(tǒng)計(jì)在綠色食品標(biāo)志使用有效期內(nèi)均計(jì)數(shù)，新增數(shù)量為當(dāng)年獲得標(biāo)志使用權(quán)的產(chǎn)品計(jì)數(shù)。各省份的經(jīng)濟(jì)統(tǒng)計(jì)數(shù)據(jù)與農(nóng)資施用量整理自《中國(guó)統(tǒng)計(jì)年鑒》《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》和布瑞克數(shù)據(jù)庫(kù)等。

（二）綠色農(nóng)產(chǎn)品認(rèn)證時(shí)空演化特征

1.綠色農(nóng)產(chǎn)品認(rèn)證基本情況

2011—2020 年中國(guó)果蔬類綠色食品數(shù)量呈現(xiàn)波動(dòng)上升態(tài)勢(shì)，2020 年果蔬綠色食品8091 個(gè)，占綠色食品總量的49.51%，相較于2011 年增加了269%。由圖1 可以看出，在2015 年之后有較大增幅，2015 年至2020 年累計(jì)新增綠色食品32 162 個(gè)。從增速看，2015 年后增長(zhǎng)率明顯上升，綠色食品數(shù)量迅速增加，后在2017 年增長(zhǎng)率回落并維持在15%水平。從新增綠色食品的生產(chǎn)經(jīng)營(yíng)主體看，2015 年后呈現(xiàn)穩(wěn)步上升的趨勢(shì)，但增速小于綠色食品的產(chǎn)品數(shù)量增速，且在2016 年達(dá)到波峰76.47%后回落至2020 年的22.2%。

圖1 2011—2020 年中國(guó)果蔬類綠色食品情況

2.空間演化趨勢(shì)

中國(guó)各省在2011—2020 年的果蔬類綠色食品數(shù)量均呈上升趨勢(shì)，2015 年后全國(guó)各地區(qū)綠色食品數(shù)量提升較快，中國(guó)綠色食品進(jìn)入一個(gè)快速發(fā)展的新階段。同時(shí)也可以看到，各省之間綠色食品數(shù)量差異在擴(kuò)大，部分西部省份果蔬類綠色食品數(shù)量較低，且在近十年內(nèi)僅存在輕微變化。2011 年中國(guó)果蔬類綠色食品發(fā)展集中在京津冀地區(qū)和山東半島，山東遼寧與京津冀占比52.41%，逐漸往江蘇、浙江等果蔬生產(chǎn)大省擴(kuò)散，且在東部沿海地區(qū)城市群表現(xiàn)出較強(qiáng)的輻射帶動(dòng)效應(yīng)，到2017 年，山東、江蘇、安徽、浙江等四地新增果蔬類綠色食品數(shù)量占到當(dāng)年全國(guó)總新增綠色食品的41.70%。2015 年后，中部、西南地區(qū)增勢(shì)明顯，2015 年中部與西部地區(qū)的新增數(shù)占全國(guó)新增數(shù)的比例是25.19%，到2020 年，中部、西部地區(qū)占比已經(jīng)達(dá)到了33.34%。2017 年后，我國(guó)果蔬類綠色食品已形成了冀江浙滬魯連片發(fā)展、西藏湖南貴州等省異軍突起的發(fā)展態(tài)勢(shì)。綠色食品發(fā)展水平在空間上存在明顯的差異性，東部沿海地區(qū)發(fā)展較早且水平較高，中西部地區(qū)在近年來(lái)迅速發(fā)展，并取得了一定的效果，這一現(xiàn)象與黃漫宇、張長(zhǎng)魯?shù)葘W(xué)者的早期研究結(jié)論較相似[16,17]。

3.空間自相關(guān)性

為檢驗(yàn)變量是否存在空間相關(guān)性，本文用Moran’s I 統(tǒng)計(jì)量進(jìn)行了空間統(tǒng)計(jì)分析，該統(tǒng)計(jì)值反映空間鄰近區(qū)域單元屬性值的相似程度，其取值范圍一般在[-1, 1]之間，越接近1，表示區(qū)域單元間的關(guān)系越密切；越接近-1，說(shuō)明區(qū)域單元間的差異越大或分布越不集中。運(yùn)用Geoda 軟件，基于地理距離空間權(quán)重矩陣，本文計(jì)算了2011—2020 年中國(guó)367 個(gè)城市果蔬類綠色食品新增水平的 Moran’s I 統(tǒng)計(jì)值（表1）?？梢园l(fā)現(xiàn)，我國(guó)綠色食品申請(qǐng)不是隨機(jī)分布的，申請(qǐng)數(shù)量近10 年來(lái)一直具有明顯的空間集聚效應(yīng)，且在2015 年后逐漸增強(qiáng)。

表1 我國(guó)城市果蔬類綠色食品數(shù)量水平和農(nóng)業(yè)環(huán)境質(zhì)量水平Moran’s I 指數(shù)

為進(jìn)一步識(shí)別各省綠色食品數(shù)量水平的局部集聚特征，本文繪制了中國(guó)城市一級(jí)的果蔬類綠色食品數(shù)量的LISA 集聚圖①相關(guān)LISA 集聚圖表在正文中略去，可向作者交流索要。。整體上，中國(guó)各區(qū)域果蔬類綠色食品以“高-高”集聚和“低-低”集聚為主，“高-高”集聚的城市數(shù)逐年增多，主要分布在東部沿海地區(qū)，如京津冀地區(qū)、山東半島、長(zhǎng)三角地區(qū)，說(shuō)明中國(guó)綠色食品發(fā)展的區(qū)域集聚態(tài)勢(shì)加強(qiáng)且向高水平發(fā)展，輻射帶動(dòng)效應(yīng)較強(qiáng)。

從時(shí)間演化可知，2011 年全國(guó)果蔬類綠色食品數(shù)量較低，僅有河北、天津和山東部分城市存在“高-高”集聚，河北和山西的交界地帶以及山東的部分城市表現(xiàn)出“低-高”集聚，甘肅、四川、廣東和云南零散幾個(gè)城市表現(xiàn)出“高-低”集聚。2013—2017 年，“高-高”聚集往江蘇、浙江等部分城市拓展，江蘇和浙江的“低-高”聚集地區(qū)也逐漸增加。黑龍江與內(nèi)蒙古自治區(qū)出現(xiàn)的“低-低”聚集地區(qū)變得不顯著，但可以看到，“高-低”集聚在西部地區(qū)如青海、四川、陜西等的占比增加。2019 年全國(guó)綠色食品數(shù)量增長(zhǎng)明顯，“高-高”聚集繼續(xù)往山西、江浙擴(kuò)張，同時(shí)安徽也出現(xiàn)了“高-高”聚集現(xiàn)象；“低-低”聚集在貴州、青海等地縮減；“低-高”聚集在“高-高”聚集地區(qū)附近；“高-低”聚集在云南有所擴(kuò)張且分布集中。

（三）農(nóng)業(yè)環(huán)境質(zhì)量的時(shí)空演化特征

1.農(nóng)業(yè)環(huán)境基本情況

本文綜合已有文獻(xiàn)的指標(biāo)體系[18-20]，選取化肥施用量和農(nóng)藥施用量作為衡量我國(guó)農(nóng)業(yè)環(huán)境質(zhì)量的代理變量。作為勞動(dòng)力密集型和農(nóng)資集約型的果蔬產(chǎn)業(yè)產(chǎn)量猛增的同時(shí)，對(duì)農(nóng)藥化肥依賴程度非常明顯。從全國(guó)來(lái)看，如表2 所示，我國(guó)農(nóng)資化肥施用量在2011—2020 年間呈現(xiàn)出先增長(zhǎng)后下降的變化趨勢(shì)，農(nóng)藥施用量在2015 年后開始下降，化肥施用量也在2016 年回落，但仍保持在較高水平。

表2 2011—2020 年我國(guó)果蔬生產(chǎn)環(huán)境質(zhì)量情況

2.空間演化特征與相關(guān)性分析

全國(guó)各省份農(nóng)藥施用量在統(tǒng)計(jì)期間的變化不明顯，施用較多省份集中在果蔬高產(chǎn)區(qū)的中東部等省份（廣東、山東、湖北），對(duì)于農(nóng)產(chǎn)品綠色認(rèn)證體系普及與推廣，仍存在較強(qiáng)的環(huán)境約束外部負(fù)效應(yīng)；而化肥施用量則有明顯的變化，2015 年后下降較為明顯的（尤其在2017 年）有吉林、天津、河北、安徽、山西、甘肅、內(nèi)蒙古、遼寧等地。

我國(guó)農(nóng)業(yè)環(huán)境質(zhì)量的全局Moran’s I 值見表1。由表1 可以發(fā)現(xiàn)：Moran’s I 指數(shù)在大部分年份顯著為正，說(shuō)明中國(guó)各市綠色環(huán)境規(guī)制存在顯著的正向空間自相關(guān)，在鄰近地區(qū)表現(xiàn)出趨同性。

四、綠色農(nóng)業(yè)發(fā)展的空間響應(yīng)機(jī)制分析

（一）變量選擇與計(jì)量模型構(gòu)建

1.變量選擇

本文選取每年新增果蔬類綠色食品數(shù)量作為被解釋變量。解釋變量方面，由于綠色食品產(chǎn)地環(huán)境標(biāo)準(zhǔn)分別對(duì)綠色食品產(chǎn)地的空氣質(zhì)量、農(nóng)田灌溉水質(zhì)量和土壤環(huán)境的質(zhì)量的各項(xiàng)指標(biāo)、濃度限值做了明確規(guī)定，特別對(duì)生產(chǎn)過(guò)程中的投入品如農(nóng)藥、肥料等生產(chǎn)資料使用方面的規(guī)定，基于已有文獻(xiàn)[16,21-23]，本文將各省農(nóng)藥/化肥施用量作為衡量綠色農(nóng)業(yè)生產(chǎn)環(huán)境質(zhì)量情況的代理變量。為了剔除內(nèi)生性問(wèn)題，將滯后一期的果蔬農(nóng)藥/化肥施用量作為核心變量。此外，將化肥施用量分解為氮肥、磷肥、鉀肥以及復(fù)合肥，考察化肥施用結(jié)構(gòu)的影響。

控制變量方面，經(jīng)濟(jì)越發(fā)達(dá)地區(qū)，對(duì)生態(tài)環(huán)境感知度越高。本文根據(jù)Ji 等的研究，選取人均GDP表示當(dāng)?shù)亟?jīng)濟(jì)發(fā)展情況[24]；第一產(chǎn)業(yè)GDP 占比表示農(nóng)業(yè)重要性程度；綠色農(nóng)業(yè)政策解釋為當(dāng)?shù)厥欠駥?shí)施過(guò)生態(tài)循環(huán)、肥藥雙減等對(duì)農(nóng)業(yè)生產(chǎn)方式存在一定規(guī)制的政策條例；政府農(nóng)林水事務(wù)支出用以表示政府對(duì)農(nóng)業(yè)、生態(tài)環(huán)境等外部公共轉(zhuǎn)移支付，可能對(duì)當(dāng)?shù)氐纳鷳B(tài)環(huán)境和農(nóng)業(yè)面源污染存在潛在影響；以居民消費(fèi)價(jià)格指數(shù)CPI 的增速來(lái)反映市場(chǎng)銷售前景；城鎮(zhèn)居民可支配收入表示市場(chǎng)消費(fèi)潛力，市場(chǎng)消費(fèi)潛力越大地區(qū)，對(duì)高品質(zhì)農(nóng)產(chǎn)品需求越大，對(duì)綠色食品需求可能存在影響。用果蔬產(chǎn)量來(lái)表示該地區(qū)的果蔬產(chǎn)業(yè)發(fā)展水平，以消除潛在影響綠色食品數(shù)量和農(nóng)藥化肥施用量的其他社會(huì)經(jīng)濟(jì)因素，盡量減少遺漏變量可能造成的高估核心解釋變量的誤差偏差。具體解釋見表3。

2.空間杜賓模型的選擇

本文首先運(yùn)用空間計(jì)量檢驗(yàn)及相關(guān)統(tǒng)計(jì)量進(jìn)行LR 檢驗(yàn)和Wald 檢驗(yàn)等，確定空間杜賓模型，Hausman 檢驗(yàn)結(jié)果為正值且通過(guò)了顯著水平為1%的檢驗(yàn)，選擇固定效應(yīng)模型，并通過(guò)效應(yīng)檢驗(yàn)在個(gè)體（ind）、時(shí)間（time）和雙固定（both）中選擇雙固定模型，受限篇幅，表格未在正文列出。模型構(gòu)建如下：

其中，Y代表綠色食品數(shù)量水平（果蔬類產(chǎn)品數(shù)），X代表解釋變量為化肥、農(nóng)藥施用量和政府綠色農(nóng)業(yè)支持政策、財(cái)政轉(zhuǎn)移支付、經(jīng)濟(jì)發(fā)展水平、市場(chǎng)消費(fèi)潛力、果蔬產(chǎn)量等控制變量，W為空間權(quán)重矩陣，ρ表示空間溢出效應(yīng)（鄰近地區(qū)的綠色食品數(shù)量水平對(duì)本地區(qū)的影響），β代表解釋變量回歸系數(shù)向量，δ代表鄰近地區(qū)的解釋變量對(duì)本地區(qū)綠色食品數(shù)量水平發(fā)展的影響，ε為隨機(jī)擾動(dòng)項(xiàng)。

（二）生態(tài)環(huán)境對(duì)綠色食品發(fā)展影響的空間溢出效應(yīng)

1.空間相關(guān)性檢驗(yàn)結(jié)果

如表4 所示，不管是否考慮化肥的施用結(jié)構(gòu)，綠色果蔬產(chǎn)品數(shù)的空間自相關(guān)系數(shù)rho 值均為正，且通過(guò)了顯著水平為1%的檢驗(yàn)，說(shuō)明綠色果蔬產(chǎn)品數(shù)在全國(guó)范圍內(nèi)確實(shí)存在顯著的空間依賴性。由此可見，地理位置的鄰近程度對(duì)我國(guó)農(nóng)業(yè)綠色食品發(fā)展水平正向作用顯著，驗(yàn)證了相鄰區(qū)域有著相似的農(nóng)業(yè)環(huán)境，農(nóng)業(yè)綠色發(fā)展存在溢出效應(yīng)和示范作用，有利于區(qū)域間資源要素的合理配置和綠色食品協(xié)同發(fā)展。

表4 空間杜賓模型結(jié)果

空間杜賓模型估計(jì)結(jié)果顯示，農(nóng)藥施用量的系數(shù)在5%水平下顯著為負(fù)，這表明環(huán)境規(guī)制政策下農(nóng)藥施用量的減少有效促進(jìn)了綠色食品與農(nóng)業(yè)綠色轉(zhuǎn)型發(fā)展。變量的空間滯后項(xiàng)中，W*AF 系數(shù)顯著為負(fù)，表明化肥施用量的降低對(duì)綠色食品數(shù)量存在正向溢出效應(yīng)。在考慮化肥施用結(jié)構(gòu)的情況下，復(fù)合肥施用的顯著增加是促進(jìn)綠色食品數(shù)量水平提升的主要原因，在考慮空間滯后作用后，氮肥的減量使用對(duì)區(qū)域綠色食品數(shù)量存在顯著的正向溢出效應(yīng)。

2.農(nóng)業(yè)生態(tài)環(huán)境對(duì)綠色認(rèn)證的影響效應(yīng)分解

為了更加全面地分析農(nóng)業(yè)綠色食品發(fā)屬水平的空間響應(yīng)機(jī)理，本文從三個(gè)效應(yīng)方面對(duì)模型系數(shù)進(jìn)行分解，具體結(jié)果見表5。

表5 空間杜賓模型的直接效應(yīng)、間接效應(yīng)與總效應(yīng)

農(nóng)藥施用量的直接效應(yīng)為-33.65，這驗(yàn)證了本文之前提出的環(huán)境規(guī)制效應(yīng)的假說(shuō)，由于減藥的農(nóng)業(yè)生產(chǎn)管理方式促進(jìn)農(nóng)業(yè)生態(tài)環(huán)境改善，更有利于生產(chǎn)經(jīng)營(yíng)主體從事綠色生態(tài)型生產(chǎn)，激勵(lì)綠色食品認(rèn)證。這一是由于國(guó)家與行業(yè)相關(guān)政策和規(guī)定，例如《農(nóng)藥管理?xiàng)l例》中對(duì)禁止施用劇毒、高毒農(nóng)藥等的宣傳，對(duì)農(nóng)業(yè)主體或者農(nóng)戶更為深入的管控；二是由于果蔬中農(nóng)藥殘留問(wèn)題更為嚴(yán)重[25]，因此其減量效果更為顯著；三是新型生物農(nóng)藥施用成本更高[26]，所以農(nóng)戶或者主體的減藥行為和效應(yīng)更明顯。

氮肥施用量的直接效應(yīng)與溢出效應(yīng)分別為-2.40 和-25.18，且其溢出效應(yīng)在1%水平顯著；磷肥和鉀肥施用量的直接效應(yīng)與溢出效應(yīng)基本為正，但未通過(guò)顯著性水平檢驗(yàn)；復(fù)合肥使用量的直接效應(yīng)在5%的顯著水平下為正，這與我國(guó)農(nóng)藥減量施用、優(yōu)化氮磷鉀肥配比的環(huán)境規(guī)制政策相一致，表明在一系列生態(tài)環(huán)境規(guī)制政策下，我國(guó)農(nóng)業(yè)綠色發(fā)展水平得到提升，且具有地區(qū)溢出效應(yīng)。

五、研究結(jié)論

本文以中國(guó)267 個(gè)地級(jí)及以上城市為研究單元，利用2011—2020 年果蔬類綠色食品數(shù)量數(shù)據(jù)，采用空間自相關(guān)分析對(duì)中國(guó)城市綠色食品發(fā)展水平差異的時(shí)空演變特征進(jìn)行了分析，并利用空間計(jì)量模型對(duì)2011—2020 年我國(guó)30 個(gè)省級(jí)綠色食品發(fā)展水平的空間演變進(jìn)行實(shí)證分析，得到以下主要結(jié)論：（1）綠色食品數(shù)量的時(shí)空演化特征明顯，2011—2020 年間中國(guó)果蔬綠色食品數(shù)呈現(xiàn)波動(dòng)上升態(tài)勢(shì)，不同省之間的差異在增大，東部部分省份增長(zhǎng)迅速。（2）中國(guó)各省市綠色食品分布有著較強(qiáng)的空間自相關(guān)，表現(xiàn)出很強(qiáng)的地區(qū)依賴性和差異性，京津冀地區(qū)、山東半島主要表現(xiàn)為“高-高”聚集，并逐漸向江浙地區(qū)擴(kuò)展；西部地區(qū)、西南地區(qū)和東北表現(xiàn)“低-低”聚集，內(nèi)蒙古、黑龍江等地出現(xiàn)“高-低”的空間聚集特征。（3）從綠色農(nóng)產(chǎn)品空間溢出路徑識(shí)別看，生態(tài)環(huán)境對(duì)綠色食品發(fā)展產(chǎn)生顯著的規(guī)制作用，農(nóng)資施用對(duì)周邊地區(qū)綠色食品發(fā)展產(chǎn)生較為明顯的溢出效應(yīng)，且化肥施用結(jié)構(gòu)的優(yōu)化調(diào)整會(huì)進(jìn)一步促進(jìn)農(nóng)業(yè)綠色發(fā)展。

基于上述研究結(jié)論，本文提出生態(tài)環(huán)境約束下的綠色食品發(fā)展的對(duì)策建議。

第一，發(fā)揮核心省域輻射帶動(dòng)作用，協(xié)調(diào)區(qū)域綠色食品發(fā)展。未來(lái)在區(qū)域農(nóng)業(yè)綠色食品發(fā)展和高質(zhì)量發(fā)展的政策制定方面要重視省市間的依賴和聯(lián)系，重點(diǎn)培育我國(guó)中西部地區(qū)核心城市，發(fā)展“低-低”聚集區(qū)域；積極推進(jìn)“高-高”聚集區(qū)域發(fā)展，發(fā)揮其對(duì)鄰近地區(qū)的帶動(dòng)作用，使周邊地區(qū)與當(dāng)?shù)毓餐l(fā)展，進(jìn)而促進(jìn)全國(guó)整體高質(zhì)量綠色農(nóng)業(yè)的協(xié)調(diào)有序發(fā)展。

第二，加大政策支持與推廣力度，改善農(nóng)業(yè)環(huán)境質(zhì)量。由于環(huán)境保護(hù)政策執(zhí)行的不到位以及農(nóng)民環(huán)境保護(hù)意識(shí)的缺乏等，優(yōu)質(zhì)可持續(xù)的農(nóng)業(yè)生產(chǎn)環(huán)境資源缺乏，應(yīng)從產(chǎn)地環(huán)境源頭治理和政策激勵(lì)兩方面入手，加強(qiáng)保護(hù)宣傳活動(dòng)，對(duì)生產(chǎn)者進(jìn)行教育，同時(shí)通過(guò)各種優(yōu)惠政策激勵(lì)各主體的綠色生產(chǎn)行為，促使農(nóng)業(yè)生產(chǎn)環(huán)境的改善，為農(nóng)業(yè)綠色與高質(zhì)量發(fā)展提供良好的自然環(huán)境基礎(chǔ)。

第三，因地制宜促進(jìn)綠色農(nóng)業(yè)高質(zhì)量發(fā)展。依據(jù)不同區(qū)域的自然經(jīng)濟(jì)特點(diǎn)，因地制宜地推行農(nóng)業(yè)綠色生產(chǎn)方式，加大力度促進(jìn)區(qū)域間要素流動(dòng)，形成資源要素的跨域優(yōu)勢(shì)互補(bǔ)和合理分工，開發(fā)生產(chǎn)健康綠色的農(nóng)產(chǎn)品，統(tǒng)一資源環(huán)境保護(hù)與農(nóng)業(yè)經(jīng)濟(jì)發(fā)展，促進(jìn)綠色、生態(tài)農(nóng)業(yè)的示范效應(yīng)發(fā)揮，實(shí)現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展。