馬 進(jìn)

一、引言與文獻(xiàn)綜述

隨著經(jīng)濟(jì)全球化發(fā)展,中國農(nóng)業(yè)部門也逐漸加入全球價(jià)值鏈分工體系之中,吸引外商直接投資(FDI)成為參與分工的重要途徑。FDI在助力東道國經(jīng)濟(jì)發(fā)展的同時(shí),也會(huì)帶來一些負(fù)面影響。長期以來,粗放的農(nóng)業(yè)發(fā)展方式使得我國農(nóng)業(yè)面源污染問題日益突出。農(nóng)業(yè)FDI的引入可能會(huì)加劇這一問題。立足中國式現(xiàn)代化背景推進(jìn)鄉(xiāng)村全面振興,在利用FDI推動(dòng)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展的同時(shí),加強(qiáng)農(nóng)業(yè)面源污染治理,是實(shí)現(xiàn)農(nóng)業(yè)綠色化轉(zhuǎn)型和推動(dòng)農(nóng)業(yè)可持續(xù)發(fā)展的重要途徑。此外,由于農(nóng)業(yè)面源污染具有通過水體傳播擴(kuò)散的特點(diǎn),鄰近地區(qū)間可能存在空間溢出效應(yīng)。那么,農(nóng)業(yè)FDI流入會(huì)給我國農(nóng)業(yè)環(huán)境帶來怎樣的影響?是否會(huì)加劇農(nóng)業(yè)面源污染?兩者之間是否存在空間關(guān)聯(lián)性?這是本文研究的主要問題。

FDI是否引致并加劇了環(huán)境污染一直備受關(guān)注且存在較大爭議,形成了“污染天堂”(1)Cole M.A.,Elliott R.J.R.,and Fredriksson P.G.,“Endogenous Pollution Havens: Does FDI Influence Environmental Regulations?” in Scandinavian Journal of Economics,Vol.108,No.1(2006),pp.157-178;Bakhsh K.,Rose S.,and Ali M.F.,et al.,“Economic Growth,CO2 Emissions,Renewable Waste and FDI Relation in Pakistan: New Evidences from 3SLS,” in Journal of Environmental Management,Vol.196(2017),pp.627-632.“污染光環(huán)”(2)宋德勇、易艷春:《外商直接投資與中國碳排放》,《中國人口·資源與環(huán)境》2011年第1期;聶飛、劉海云:《FDI、環(huán)境污染與經(jīng)濟(jì)增長的相關(guān)性研究——基于動(dòng)態(tài)聯(lián)立方程模型的實(shí)證檢驗(yàn)》,《國際貿(mào)易問題》2015年第2期。、綜合環(huán)境效應(yīng)論等不同理論。綜合環(huán)境效應(yīng)理論認(rèn)為FDI對(duì)環(huán)境的影響是由多種效應(yīng)共同作用的結(jié)果,其影響方向不確定(3)周杰琦、汪同三:《外商直接投資、經(jīng)濟(jì)增長和CO2排放——基于中國省際數(shù)據(jù)的實(shí)證研究》,《北京理工大學(xué)學(xué)報(bào)(社會(huì)科學(xué)版)》2014年第3期。。根據(jù)該理論,農(nóng)業(yè)FDI主要通過規(guī)模效應(yīng)、結(jié)構(gòu)效應(yīng)和技術(shù)效應(yīng)三種路徑影響農(nóng)業(yè)生產(chǎn)環(huán)境。從規(guī)模效應(yīng)來看,大部分學(xué)者認(rèn)為農(nóng)業(yè)FDI對(duì)農(nóng)業(yè)經(jīng)濟(jì)增長存在正向促進(jìn)效應(yīng)(4)劉乃郗、韓一軍、王萍萍:《FDI是否提高了中國農(nóng)業(yè)企業(yè)全要素生產(chǎn)率?——來自99801家農(nóng)業(yè)企業(yè)面板數(shù)據(jù)的證據(jù)》,《中國農(nóng)村經(jīng)濟(jì)》2018年第4期。,而農(nóng)業(yè)經(jīng)濟(jì)增長對(duì)于環(huán)境污染的影響方向并不確定(5)馬乃毅、羅珺仁:《中國西部地區(qū)農(nóng)業(yè)經(jīng)濟(jì)增長與環(huán)境污染相關(guān)性研究》,《浙江農(nóng)業(yè)學(xué)報(bào)》2015年第9期。。從結(jié)構(gòu)效應(yīng)來看,農(nóng)業(yè)FDI可以影響農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)和種植業(yè)結(jié)構(gòu)的調(diào)整(6)李泳:《國際直接投資與中國農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)》,《中國農(nóng)村經(jīng)濟(jì)》2006年第5期。,而結(jié)構(gòu)調(diào)整又可以通過集約化經(jīng)營和對(duì)產(chǎn)污單元的調(diào)整影響農(nóng)業(yè)污染排放(7)楊軍、李建琴:《福建省農(nóng)業(yè)經(jīng)濟(jì)增長、農(nóng)業(yè)結(jié)構(gòu)與面源污染關(guān)系研究》,《中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文)》2020年第8期。。從技術(shù)效應(yīng)來看,農(nóng)業(yè)FDI具有技術(shù)溢出效應(yīng)。韋開蕾發(fā)現(xiàn)FDI可以通過競爭示范效應(yīng)、人力資本效應(yīng)以及資本規(guī)模效應(yīng)實(shí)現(xiàn)技術(shù)外溢(8)韋開蕾:《基于FDI溢出效應(yīng)的農(nóng)業(yè)生產(chǎn)技術(shù)效率的地區(qū)差異》,《社會(huì)科學(xué)家》2015年第10期。;但文雁兵發(fā)現(xiàn)在全國層面上,外資和跨國公司逐漸涉足中國農(nóng)業(yè)的研發(fā)和生產(chǎn)環(huán)節(jié),導(dǎo)致農(nóng)業(yè)FDI產(chǎn)生了負(fù)向技術(shù)溢出,抑制了農(nóng)業(yè)技術(shù)的提高(9)文雁兵:《我國農(nóng)業(yè)科技自主創(chuàng)新能力研究——基于產(chǎn)業(yè)關(guān)聯(lián)效應(yīng)和FDI技術(shù)溢出視角》,《科學(xué)學(xué)研究》2015年第7期。。通常認(rèn)為,農(nóng)業(yè)技術(shù)創(chuàng)新可以減少污染排放,提升農(nóng)業(yè)發(fā)展質(zhì)量(10)李紅莉、張俊飚、羅斯炫等:《農(nóng)業(yè)技術(shù)創(chuàng)新對(duì)農(nóng)業(yè)發(fā)展質(zhì)量的影響及作用機(jī)制——基于空間視角的經(jīng)驗(yàn)分析》,《研究與發(fā)展管理》2021年第2期。。

近年來,部分學(xué)者對(duì)農(nóng)業(yè)面源污染的動(dòng)態(tài)演變和空間差異進(jìn)行了研究(11)魯慶堯、王樹進(jìn):《我國農(nóng)業(yè)面源污染的空間相關(guān)性及影響因素研究》,《經(jīng)濟(jì)問題》2015年第12期。,并在空間維度下從經(jīng)濟(jì)增長、產(chǎn)業(yè)結(jié)構(gòu)、技術(shù)進(jìn)步等方面分析了不同經(jīng)濟(jì)因素對(duì)農(nóng)業(yè)面源污染的影響(12)梁偉健、江華、廖文玉等:《農(nóng)業(yè)面源污染與農(nóng)業(yè)經(jīng)濟(jì)增長的空間互動(dòng)效應(yīng)》,《江淮論壇》2018年第3期;閆桂權(quán)、何玉成、張曉恒:《綠色技術(shù)進(jìn)步、農(nóng)業(yè)經(jīng)濟(jì)增長與污染空間溢出——來自中國農(nóng)業(yè)水資源利用的證據(jù)》,《長江流域資源與環(huán)境》2019年第12期。。也有部分學(xué)者開始利用空間計(jì)量模型探討FDI對(duì)環(huán)境的影響(13)徐春華、劉力:《FDI、政府消費(fèi)與CO2排放——基于36國貿(mào)易空間權(quán)重矩陣的空間杜賓模型分析》,《國際經(jīng)貿(mào)探索》2016年第1期;董會(huì)忠、韓沅剛:《開放與綠色理念下如何提升工業(yè)生態(tài)效率?——基于 “污染天堂” 假說的驗(yàn)證》,《商業(yè)研究》2020年第12期。,但是涉及農(nóng)業(yè)環(huán)境的研究較少。

現(xiàn)有文獻(xiàn)對(duì)本文研究具有積極的借鑒意義,但存在以下不足:一是FDI環(huán)境效應(yīng)的研究對(duì)象多集中于工業(yè)部門,農(nóng)業(yè)FDI的環(huán)境效應(yīng)特別是農(nóng)業(yè)FDI對(duì)農(nóng)業(yè)面源污染的研究較少;二是在衡量農(nóng)業(yè)環(huán)境時(shí)往往采用化肥等單一指標(biāo),農(nóng)業(yè)面源污染核算方法和相關(guān)系數(shù)比較陳舊;三是大多關(guān)注農(nóng)業(yè)FDI數(shù)量的影響,如采用外商投資額衡量農(nóng)業(yè)FDI水平,較少探究其質(zhì)量的影響;四是通?？紤]污染物對(duì)本地環(huán)境造成的影響,對(duì)空間溢出效應(yīng)的研究尚不充分。

本文在已有研究基礎(chǔ)上,采用單元調(diào)查法和輸出系數(shù)法核算了2005—2020年我國部分省域農(nóng)業(yè)面源污染物排放量,構(gòu)建FDI數(shù)量和質(zhì)量雙重指標(biāo),將空間自相關(guān)性納入計(jì)量模型,采用空間杜賓模型實(shí)證分析農(nóng)業(yè)FDI對(duì)我國農(nóng)業(yè)面源污染的直接影響和空間效應(yīng),對(duì)促進(jìn)農(nóng)業(yè)綠色可持續(xù)發(fā)展提出政策建議。

二、理論分析

下面從農(nóng)業(yè)FDI對(duì)本地農(nóng)業(yè)面源污染的直接影響和鄰地FDI流入對(duì)本地農(nóng)業(yè)面源污染的空間溢出效應(yīng)兩方面進(jìn)行理論分析。

(一)農(nóng)業(yè)FDI對(duì)本地農(nóng)業(yè)面源污染的影響

農(nóng)業(yè)FDI對(duì)農(nóng)業(yè)面源污染的影響可以從農(nóng)業(yè)FDI數(shù)量和質(zhì)量兩方面來分析。FDI數(shù)量是指外資利用規(guī)模;FDI質(zhì)量的內(nèi)涵極為豐富,本文所涉FDI質(zhì)量主要指生產(chǎn)溢出效應(yīng),包括對(duì)經(jīng)濟(jì)的貢獻(xiàn)和技術(shù)外溢。

1.農(nóng)業(yè)FDI數(shù)量對(duì)農(nóng)業(yè)面源污染的影響

農(nóng)業(yè)FDI可以彌補(bǔ)東道國資金相對(duì)不足的缺陷,帶動(dòng)農(nóng)業(yè)經(jīng)濟(jì)發(fā)展,進(jìn)而影響農(nóng)業(yè)面源污染排放。FDI會(huì)導(dǎo)致一國不同生產(chǎn)部門產(chǎn)品的相對(duì)價(jià)格發(fā)生變動(dòng),使其比較優(yōu)勢產(chǎn)業(yè)結(jié)構(gòu)發(fā)生改變,從而影響一國的環(huán)境狀況。FDI的流入還會(huì)帶來更為清潔的生產(chǎn)技術(shù)。通常認(rèn)為,FDI能通過技術(shù)引進(jìn)、技術(shù)溢出、競爭示范和人才流動(dòng)等方式對(duì)農(nóng)業(yè)生產(chǎn)技術(shù)和污染治理技術(shù)產(chǎn)生正外部性,有效減少農(nóng)業(yè)面源污染排放。據(jù)此,本文提出如下假設(shè):

假設(shè)1:農(nóng)業(yè)FDI流入數(shù)量提升會(huì)影響本地農(nóng)業(yè)面源污染排放水平。

2.農(nóng)業(yè)FDI質(zhì)量對(duì)農(nóng)業(yè)面源污染的影響

農(nóng)業(yè)FDI質(zhì)量越高,其對(duì)農(nóng)業(yè)經(jīng)濟(jì)的貢獻(xiàn)率越高,對(duì)經(jīng)濟(jì)發(fā)展的拉動(dòng)作用越強(qiáng)。因此,在同等數(shù)量下,質(zhì)量越高的FDI對(duì)農(nóng)業(yè)面源污染的影響越大。一方面,FDI質(zhì)量的提升可能會(huì)帶來農(nóng)業(yè)生產(chǎn)規(guī)模的擴(kuò)大,進(jìn)而導(dǎo)致農(nóng)業(yè)環(huán)境的惡化。另一方面,當(dāng)農(nóng)業(yè)FDI質(zhì)量越高時(shí),其流入帶來的清潔技術(shù)含量可能也會(huì)越高,產(chǎn)生的技術(shù)溢出效應(yīng)也會(huì)越強(qiáng),將減少農(nóng)業(yè)污染的排放。據(jù)此,本文提出如下假設(shè):

假設(shè)2:農(nóng)業(yè)FDI流入質(zhì)量提升會(huì)影響本地農(nóng)業(yè)面源污染排放水平。

(二)鄰地農(nóng)業(yè)FDI流入對(duì)本地農(nóng)業(yè)面源污染的影響

相鄰地區(qū)農(nóng)業(yè)FDI流入對(duì)本地農(nóng)業(yè)面源污染的影響途徑大致可以歸納為以下三個(gè)方面:一是鄰地農(nóng)業(yè)FDI通過規(guī)模效應(yīng)對(duì)其農(nóng)業(yè)面源污染的影響會(huì)加劇本地農(nóng)業(yè)面源污染。農(nóng)業(yè)面源污染以水為載體,可溶解的或固態(tài)的農(nóng)業(yè)污染物在降水沖刷或灌溉的作用下,通過地表徑流、農(nóng)田排水和地下滲漏匯入水體,不僅會(huì)造成當(dāng)?shù)厮h(huán)境污染,也極易對(duì)具有相似水系和地形地貌特征的周邊環(huán)境造成破壞。二是相鄰地區(qū)通常具備相似的農(nóng)業(yè)生產(chǎn)條件,鄰地農(nóng)業(yè)FDI對(duì)其生產(chǎn)結(jié)構(gòu)的影響可能會(huì)對(duì)本地起到示范作用,從而影響本地的生產(chǎn)行為。三是鄰地FDI流入帶來的清潔技術(shù)提升,可能會(huì)通過示范作用激勵(lì)本地農(nóng)業(yè)生產(chǎn)者生產(chǎn)更為清潔且符合市場預(yù)期的農(nóng)產(chǎn)品,或者憑借競爭效應(yīng)鼓勵(lì)農(nóng)業(yè)生產(chǎn)者改變生產(chǎn)方式,探索在同等生產(chǎn)水平下減少污染副產(chǎn)品產(chǎn)出的方法。同時(shí),鄰近地區(qū)無論在經(jīng)濟(jì)技術(shù)水平還是在生態(tài)環(huán)境資源方面都較為相似,且由于交通日益便利,進(jìn)一步增加了技術(shù)溢出的可能性。據(jù)此,本文提出如下假設(shè):

假設(shè)3:鄰地農(nóng)業(yè)FDI流入會(huì)影響本地農(nóng)業(yè)面源污染排放水平。

三、研究設(shè)計(jì)

(一)模型構(gòu)建

在可拓展STIRPAT模型的基礎(chǔ)上,借鑒其他學(xué)者(14)高明、吳雪萍、郭施宏:《城市化進(jìn)程、環(huán)境規(guī)制與大氣污染——基于STIRPAT模型的實(shí)證分析》,《工業(yè)技術(shù)經(jīng)濟(jì)》2016年第9期;吳義根、馮開文、李谷成:《人口增長、結(jié)構(gòu)調(diào)整與農(nóng)業(yè)面源污染——基于空間面板STIRPAT模型的實(shí)證研究》,《農(nóng)業(yè)技術(shù)經(jīng)濟(jì)》2017年第3期。的研究思路,引入環(huán)境規(guī)制強(qiáng)度、農(nóng)業(yè)結(jié)構(gòu)等相關(guān)指標(biāo),構(gòu)建實(shí)證分析模型如下:

lnEIit=β0+β1lnFDIit+β2lnSit+β3lnPit+β4lnTit+β5lnRit+β6lnStrit+μit

(1)

其中,i和t分別表示地區(qū)和時(shí)間。EI為被解釋變量,表征農(nóng)業(yè)面源污染水平。FDI為核心解釋變量,表征農(nóng)業(yè)FDI水平。S表示農(nóng)業(yè)經(jīng)濟(jì)規(guī)模,P表示農(nóng)村人口密度,T表示農(nóng)業(yè)技術(shù)水平,R表示環(huán)境規(guī)制強(qiáng)度,Str表示農(nóng)業(yè)結(jié)構(gòu)。

(二)變量說明

1.被解釋變量

本文選擇能夠表征農(nóng)業(yè)面源污染在土地上富集程度的單位面積農(nóng)業(yè)面源污染等標(biāo)排放總量,即農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度(EI)作為被解釋變量。根據(jù)我國農(nóng)業(yè)面源污染的來源和特征,將產(chǎn)污單元?jiǎng)澐譃榛?、畜禽養(yǎng)殖、農(nóng)田固體廢棄物、水產(chǎn)養(yǎng)殖和農(nóng)村生活五大類,對(duì)全國農(nóng)業(yè)面源污染不同污染源的化學(xué)需氧量(COD)、總氮(TN) 和總磷(TP)排放量進(jìn)行估算。其中,化肥、水產(chǎn)養(yǎng)殖、農(nóng)村生活污染的核算采用2021年生態(tài)環(huán)境部發(fā)布的《排放源統(tǒng)計(jì)調(diào)查產(chǎn)排污核算方法和系數(shù)手冊(cè)》中的輸出系數(shù)方法;畜禽養(yǎng)殖和農(nóng)田固體廢棄物污染的核算采用單元調(diào)查法(15)陳敏鵬、陳吉寧、賴斯蕓:《中國農(nóng)業(yè)和農(nóng)村污染的清單分析與空間特征識(shí)別》,《中國環(huán)境科學(xué)》2006年第6期。。

由于各類污染物的環(huán)境標(biāo)準(zhǔn)和污染危害程度不同,無法在同一尺度上對(duì)其污染貢獻(xiàn)度進(jìn)行比較,因此引入等標(biāo)排放量指標(biāo),其計(jì)算公式為:

Pi=Qi/Coi

(2)

其中,Pi為農(nóng)業(yè)面源污染物的等標(biāo)排放量,Qi為污染物的年排放量,Coi為該污染物的環(huán)境質(zhì)量標(biāo)準(zhǔn)。根據(jù)《地表水環(huán)境質(zhì)量標(biāo)準(zhǔn)》(GB3838-2002),COD、TN、TP的III類標(biāo)準(zhǔn)限值分別為20mg·L-1、1.0mg·L-1、0.2mg·L-1。

農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度EI的計(jì)算公式為:

EI=E/AL

(3)

其中,E為COD、TN和TP的等標(biāo)排放量總和,AL為研究區(qū)域的農(nóng)用地面積。

2.核心解釋變量

核心解釋變量為農(nóng)業(yè)FDI數(shù)量(FDIA)和質(zhì)量(FDIP)。農(nóng)業(yè)FDI數(shù)量采用各省份農(nóng)林牧副漁業(yè)實(shí)際利用外商投資額衡量。對(duì)于FDI質(zhì)量內(nèi)涵及測度的系統(tǒng)研究較少,尚未達(dá)成一致。本文參考Wang 和Luo(16)Wang X.and Luo Y.,“Has Technological Innovation Capability Addressed Environmental Pollution from the Dual Perspective of FDI Quantity and Quality? Evidence from China,” in Journal of Cleaner Production,Vol.258 (2020),pp.1-13.的方法,構(gòu)建外資業(yè)績指數(shù)衡量FDI質(zhì)量差異。因?yàn)橥赓Y業(yè)績指數(shù)可以在一定程度上反映一國外資的技術(shù)含量、技術(shù)管理實(shí)力和外商技術(shù)轉(zhuǎn)移的意愿。具體計(jì)算公式如下:

(4)

其中,FDIAi表示i地區(qū)的農(nóng)林牧副漁業(yè)實(shí)際利用外商投資額,FDIA表示全國農(nóng)林牧副漁業(yè)實(shí)際利用外商投資額,GDPi表示i地區(qū)生產(chǎn)總值,GDP表示國內(nèi)生產(chǎn)總值。

3.控制變量

根據(jù)已有研究,選擇控制變量如下:農(nóng)業(yè)經(jīng)濟(jì)規(guī)模(S),以各地農(nóng)林牧副漁業(yè)總產(chǎn)值與其農(nóng)用地面積的比值表示;農(nóng)村人口密度(P),以各地年末常住人口數(shù)與其農(nóng)用地面積比值表示;農(nóng)業(yè)技術(shù)水平(T),以各地農(nóng)業(yè)機(jī)械總動(dòng)力與其農(nóng)用地面積比值表示;環(huán)境規(guī)制強(qiáng)度(R),以各地環(huán)境污染治理投資總額占其GDP的比重表示;農(nóng)業(yè)結(jié)構(gòu)(Str),以各地牧業(yè)總產(chǎn)值與其農(nóng)林牧副漁業(yè)總產(chǎn)值的比值表示。

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

因數(shù)據(jù)限制,本文選取了2005—2020年全國27個(gè)省份(不包括吉林、四川、青海、西藏以及香港、澳門和臺(tái)灣)的平衡面板數(shù)據(jù)作為研究樣本。農(nóng)業(yè)面源污染的各單元調(diào)查指標(biāo)數(shù)據(jù)來自《中國統(tǒng)計(jì)年鑒》;相關(guān)系數(shù)參考大量已有文獻(xiàn)資料得到;農(nóng)業(yè)FDI數(shù)據(jù)來自各省份統(tǒng)計(jì)年鑒,實(shí)際利用外商投資額以人民幣為計(jì)價(jià)單位,使用中美匯率年平均價(jià)進(jìn)行換算,并以2005年為基期采用居民消費(fèi)價(jià)格指數(shù)進(jìn)行平減處理,其他控制變量也進(jìn)行了同樣的平減處理;環(huán)境污染治理投資總額數(shù)據(jù)來自《中國環(huán)境統(tǒng)計(jì)年鑒》;農(nóng)用地面積數(shù)據(jù)來自《中國農(nóng)村統(tǒng)計(jì)年鑒》;其他數(shù)據(jù)來自國家統(tǒng)計(jì)局的地區(qū)分省份年度數(shù)據(jù)。為消除異方差和個(gè)別極端值的影響,對(duì)各變量進(jìn)行了對(duì)數(shù)處理。

四、空間自相關(guān)性檢驗(yàn)和計(jì)量模型選擇

(一)空間自相關(guān)檢驗(yàn)

1.全局莫蘭指數(shù)分析

為判斷我國省域農(nóng)業(yè)面源污染的空間相關(guān)性,利用0—1空間鄰接權(quán)重矩陣,對(duì)2005—2020年我國27個(gè)省份的農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度進(jìn)行空間自相關(guān)檢驗(yàn)(17)囿于篇幅,相關(guān)結(jié)果未列示,留存?zhèn)渌?。。由農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度的全局莫蘭指數(shù)結(jié)果可知,2005—2020年27個(gè)省份的農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度全部通過了1%的顯著性檢驗(yàn),呈現(xiàn)出顯著的空間正相關(guān)性,說明農(nóng)業(yè)面源污染在省域間存在空間溢出效應(yīng)。

2.局部莫蘭指數(shù)分析

基于我國27個(gè)省份的鄰接關(guān)系矩陣,繪制2005年和2020年農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度的莫蘭散點(diǎn)圖,如圖1和圖2所示。由其直線斜率可知,我國27個(gè)省份的農(nóng)業(yè)面源污染具有空間正自相關(guān)效應(yīng);從四個(gè)象限的分布情況可以看出,我國大部分地區(qū)呈現(xiàn)高—高聚集或低—低聚集模式。

圖1 2005年農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度莫蘭散點(diǎn)圖

圖2 2020年農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度莫蘭散點(diǎn)圖

根據(jù)散點(diǎn)圖中各省份的分布象限,得到我國農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度的聚集模式。上海、江蘇、山東、河南、安徽、河北和海南等地屬于高—高聚集模式,貴州、重慶、山西、陜西、內(nèi)蒙古、甘肅、寧夏、新疆和黑龍江等地屬于低—低聚集模式。湖北、湖南、廣東、北京、天津等地區(qū)的聚集模式發(fā)生轉(zhuǎn)變。其中,湖北由低—高聚集模式轉(zhuǎn)變?yōu)楦摺呔奂Ｊ?說明湖北的農(nóng)業(yè)面源污染受到鄰近地區(qū)的影響增加;湖南、廣東則由高—低聚集模式轉(zhuǎn)變?yōu)楦摺呔奂Ｊ?說明近年來以上地區(qū)的農(nóng)業(yè)面源污染在一定程度上影響了周邊地區(qū),導(dǎo)致周邊地區(qū)的面源污染加重;北京由高—高聚集模式轉(zhuǎn)變?yōu)榈汀呔奂Ｊ?說明近年來北京市對(duì)面源污染的治理初顯成效,與周邊地區(qū)相比農(nóng)業(yè)面源污染減少;天津由高—高聚集模式轉(zhuǎn)變?yōu)楦摺途奂Ｊ?與周邊地區(qū)相比農(nóng)業(yè)面源污染增加。

(二)空間計(jì)量模型選擇

首先進(jìn)行LM檢驗(yàn)、LR檢驗(yàn)、Wald檢驗(yàn)以及Hausman檢驗(yàn),根據(jù)檢驗(yàn)結(jié)果,本文選擇固定效應(yīng)的空間杜賓模型進(jìn)行分析。結(jié)合(1)式的設(shè)定,具體空間計(jì)量模型如下:

農(nóng)業(yè)FDI數(shù)量:

lnEIit=ρWlnEIit+β1lnFDIAit+β2lnSit+β3lnPit+β4lnTit+β5lnRit+β6lnStrit

+θ1WlnFDIAit+θ2WlnSit+θ3WlnPit+θ4WlnTit+θ5WlnRit+θ6WlnStrit+εit

(5)

農(nóng)業(yè)FDI質(zhì)量:

lnEIit=ρWlnEIit+β1lnFDIPit+β2lnSit+β3lnPit+β4lnTit+β6lnRit+β6lnStrit

+θ1WlnFDIPit+θ2WlnSit+θ3WlnPit+θ4WlnTit+θ5WlnRit+θ6WlnStrit+εit

(6)

其中,W為空間權(quán)重矩陣,本文采用0—1鄰接空間矩陣,ρ為空間自回歸系數(shù),β、θ為待估計(jì)參數(shù)。設(shè)Y為被解釋變量,X為解釋變量,ρWY為空間滯后項(xiàng),反映鄰接區(qū)域被解釋變量對(duì)該區(qū)域被解釋變量的影響;βX反映該區(qū)域解釋變量對(duì)被解釋變量的影響;θW反映鄰接區(qū)域解釋變量對(duì)該區(qū)域被解釋變量的空間影響。

五、實(shí)證結(jié)果分析

(一)基準(zhǔn)回歸結(jié)果

我國農(nóng)業(yè)FDI數(shù)量和質(zhì)量對(duì)農(nóng)業(yè)面源污染影響的空間杜賓模型實(shí)證結(jié)果如表1和表2所示。由回歸結(jié)果可知,農(nóng)業(yè)FDI數(shù)量和質(zhì)量的時(shí)間固定效應(yīng)模型的擬合優(yōu)度及相關(guān)解釋變量的顯著性均優(yōu)于空間固定效應(yīng)模型和雙向固定效應(yīng)模型。從經(jīng)濟(jì)意義角度考慮,農(nóng)業(yè)面源污染等標(biāo)排放強(qiáng)度受時(shí)間、政策等變化影響比受地理因素影響更為明顯(18)李新、尚杰:《基于空間效應(yīng)視角的農(nóng)業(yè)經(jīng)濟(jì)增長與種植業(yè)面源污染排放的實(shí)證分析》,《中國生態(tài)農(nóng)業(yè)學(xué)報(bào)(中英文)》2022年第9期。。此外,作為本文的核心解釋變量,農(nóng)業(yè)FDI的流入規(guī)模主要受市場準(zhǔn)入標(biāo)準(zhǔn)和外商投資壁壘等影響。因此,本文選擇時(shí)間固定的空間杜賓模型進(jìn)行具體分析。

表1 FDI數(shù)量對(duì)農(nóng)業(yè)面源污染影響的空間杜賓模型回歸結(jié)果

表2 FDI質(zhì)量對(duì)農(nóng)業(yè)面源污染影響的空間杜賓模型回歸結(jié)果

在時(shí)間固定效應(yīng)模型中,農(nóng)業(yè)FDI數(shù)量(lnFDIA)和質(zhì)量(lnFDIP)的系數(shù)均顯著為正,說明兩者上升都會(huì)增加我國農(nóng)業(yè)面源污染排放,驗(yàn)證了本文理論假設(shè)1和假設(shè)2?？臻g系數(shù)WlnFDIA和WlnFDIP均顯著為正,說明相鄰省份農(nóng)業(yè)FDI數(shù)量和質(zhì)量的提升都會(huì)對(duì)本地農(nóng)業(yè)面源污染產(chǎn)生正向的空間傳導(dǎo)作用,即鄰地農(nóng)業(yè)FDI的流入會(huì)惡化本地的農(nóng)業(yè)環(huán)境狀況,驗(yàn)證了假設(shè)3。

根據(jù)實(shí)證結(jié)果推測,農(nóng)業(yè)FDI數(shù)量通過規(guī)模效應(yīng)加劇了我國農(nóng)業(yè)環(huán)境壓力,而由此產(chǎn)生的負(fù)面影響超過了技術(shù)效應(yīng)帶來的正面影響。這可能是由于農(nóng)業(yè)本身的弱質(zhì)性,導(dǎo)致其短期帶來經(jīng)濟(jì)效益的能力較低,吸引FDI的技術(shù)含量也相對(duì)較低,難以發(fā)揮清潔技術(shù)抑制污染排放的能力。此外,推測農(nóng)業(yè)FDI數(shù)量通過結(jié)構(gòu)效應(yīng)對(duì)我國農(nóng)業(yè)環(huán)境產(chǎn)生負(fù)面影響。這是因?yàn)檗r(nóng)業(yè)領(lǐng)域內(nèi)經(jīng)濟(jì)效益較高的部門為畜禽養(yǎng)殖和農(nóng)產(chǎn)品初加工等污染密集型行業(yè),在高利潤回報(bào)率的驅(qū)使下,農(nóng)業(yè)FDI更多地進(jìn)入以上領(lǐng)域,致使養(yǎng)殖和種植需求增加,從而加劇農(nóng)業(yè)面源污染。同時(shí),農(nóng)業(yè)FDI質(zhì)量的提升仍未帶來環(huán)境的改善,說明現(xiàn)階段農(nóng)業(yè)FDI質(zhì)量提升引致的清潔技術(shù)含量不足,也未抵消規(guī)模效應(yīng)所帶來的負(fù)面影響。

從農(nóng)業(yè)FDI對(duì)農(nóng)業(yè)面源污染的空間傳導(dǎo)效應(yīng)來看,鄰地農(nóng)業(yè)FDI流入會(huì)對(duì)本地農(nóng)業(yè)環(huán)境產(chǎn)生負(fù)面影響。鄰地農(nóng)業(yè)FDI流入可能會(huì)通過跨省貿(mào)易造成對(duì)本地農(nóng)產(chǎn)品需求的增加,從而增加本地農(nóng)業(yè)面源污染排放。此外,鄰地農(nóng)業(yè)FDI通過規(guī)模效應(yīng)增加其面源污染排放,該污染通過徑流、淋溶等方式進(jìn)入本地,造成本地農(nóng)業(yè)環(huán)境惡化。

(二)空間效應(yīng)分解

為更準(zhǔn)確地反映農(nóng)業(yè)FDI與農(nóng)業(yè)面源污染的關(guān)系,用偏微分方法對(duì)時(shí)間固定效應(yīng)模型的空間溢出效應(yīng)進(jìn)行分解,結(jié)果如表3和表4所示。

表3 農(nóng)業(yè)FDI數(shù)量對(duì)農(nóng)業(yè)面源污染影響的空間杜賓模型效應(yīng)分解

表3和表4中,直接效應(yīng)是指該地區(qū)內(nèi)各因素變化對(duì)本地區(qū)農(nóng)業(yè)面源污染排放的影響。間接效應(yīng)是指鄰接地區(qū)各因素變化對(duì)本地區(qū)農(nóng)業(yè)面源污染排放的影響。總效應(yīng)為直接效應(yīng)和間接效應(yīng)二者之和,反映的是我國農(nóng)業(yè)面源污染受自身驅(qū)動(dòng)和鄰地交互作用的總影響。

1.農(nóng)業(yè)FDI數(shù)量與FDI質(zhì)量的空間效應(yīng)分解結(jié)果

由表3和表4可知,農(nóng)業(yè)FDI數(shù)量的總效應(yīng)為0.0716,其分解的直接效應(yīng)和間接效應(yīng)分別為0.0220和0.0496,區(qū)域內(nèi)的空間溢出效應(yīng)為直接效應(yīng)的兩倍以上;FDI質(zhì)量的總效應(yīng)為0.0678,小于FDI數(shù)量的總效應(yīng),其分解的直接效應(yīng)和間接效應(yīng)分別為0.0374和0.0304。上述變量的直接效應(yīng)、間接效應(yīng)以及總效應(yīng)均顯著為正,說明農(nóng)業(yè)FDI是影響農(nóng)業(yè)面源污染的重要因素之一。

2.其他控制變量的空間效應(yīng)分解結(jié)果

農(nóng)業(yè)經(jīng)濟(jì)規(guī)模對(duì)農(nóng)業(yè)面源污染的直接效應(yīng)、間接效應(yīng)和總效應(yīng)均顯著為正,說明農(nóng)業(yè)經(jīng)濟(jì)規(guī)模尚未跨越環(huán)境庫茲涅茨曲線拐點(diǎn),其擴(kuò)張會(huì)惡化農(nóng)業(yè)環(huán)境。鄰接地區(qū)可通過經(jīng)濟(jì)輻射影響本地農(nóng)業(yè)面源污染排放,這是由于農(nóng)業(yè)經(jīng)濟(jì)發(fā)展所帶來的化肥、農(nóng)藥等生產(chǎn)要素投入增加和畜禽養(yǎng)殖規(guī)模擴(kuò)大。

農(nóng)村人口密度對(duì)農(nóng)業(yè)面源污染的直接效應(yīng)和總效應(yīng)均顯著為正。這是因?yàn)檗r(nóng)村人口密度越大的省份,其生活污水和生活垃圾所產(chǎn)生的污染強(qiáng)度會(huì)越高。間接效應(yīng)為負(fù)且不顯著,說明農(nóng)村人口密度的空間溢出效應(yīng)不明顯。

農(nóng)業(yè)技術(shù)水平對(duì)農(nóng)業(yè)面源污染的直接效應(yīng)不顯著,說明本地簡單農(nóng)業(yè)機(jī)械的數(shù)量積累和普及率提高不能帶來農(nóng)業(yè)環(huán)境的顯著改善,主要原因是本地綠色機(jī)械化程度不足。間接效應(yīng)和總效應(yīng)顯著為負(fù),是因?yàn)檗r(nóng)業(yè)技術(shù)水平高的地區(qū)能起到示范作用,促進(jìn)周邊區(qū)域?qū)崿F(xiàn)農(nóng)業(yè)資源有效利用,形成省域間綠色農(nóng)業(yè)的良性競爭。從總體看,農(nóng)業(yè)技術(shù)水平的提高對(duì)農(nóng)業(yè)面源污染具有顯著的抑制作用。

環(huán)境規(guī)制強(qiáng)度對(duì)農(nóng)業(yè)面源污染的直接效應(yīng)和總效應(yīng)顯著為負(fù),說明本地環(huán)境標(biāo)準(zhǔn)越高,對(duì)環(huán)境污染治理投資的比例越大,農(nóng)業(yè)面源污染排放越少。這是因?yàn)檗r(nóng)戶在環(huán)境偏好較強(qiáng)的地區(qū)會(huì)更傾向于選擇綠色生產(chǎn)方式,進(jìn)而使農(nóng)業(yè)環(huán)境得到改善。間接效應(yīng)為正且不顯著,可能是隨著環(huán)境管制的加強(qiáng),為了趨利避害,農(nóng)業(yè)產(chǎn)業(yè)會(huì)向周邊環(huán)境規(guī)制水平較低的區(qū)域進(jìn)行轉(zhuǎn)移,進(jìn)而加重鄰地污染。但是環(huán)境規(guī)制的擴(kuò)散效果不明顯,說明省域間存在著較強(qiáng)的行政壁壘。

農(nóng)業(yè)結(jié)構(gòu)對(duì)農(nóng)業(yè)面源污染影響的直接效應(yīng)、間接效應(yīng)和總效應(yīng)均顯著為正。隨著肉蛋奶等產(chǎn)品需求的擴(kuò)大,我國牧業(yè)產(chǎn)值及占農(nóng)林牧漁業(yè)總產(chǎn)值比重不斷上升。通過對(duì)歷年我國農(nóng)業(yè)面源污染排放量的核算,發(fā)現(xiàn)畜禽養(yǎng)殖為我國最大污染源,畜禽養(yǎng)殖比重的提高會(huì)增加農(nóng)業(yè)面源污染。間接效應(yīng)顯著為正,說明在經(jīng)濟(jì)收益的刺激下,鄰接地區(qū)農(nóng)業(yè)結(jié)構(gòu)調(diào)整可能會(huì)推動(dòng)本地區(qū)農(nóng)業(yè)結(jié)構(gòu)的同向變化,從而加劇本地農(nóng)業(yè)面源污染。

(三)穩(wěn)健性檢驗(yàn)(19)囿于篇幅,相關(guān)檢驗(yàn)結(jié)果未列示,留存?zhèn)渌鳌?/h3>

本文采取以下兩種方法,檢驗(yàn)回歸結(jié)果的穩(wěn)健性:第一,更換空間權(quán)重矩陣。根據(jù)各省份省會(huì)的經(jīng)緯坐標(biāo)之間的距離,構(gòu)建反距離權(quán)重矩陣,將該矩陣替換鄰接空間權(quán)重矩陣代入原模型中。第二,縮短時(shí)間窗口。通過對(duì)我國農(nóng)業(yè)面源污染治理政策演進(jìn)特征進(jìn)行分析發(fā)現(xiàn),2011年開始進(jìn)入農(nóng)業(yè)面源污染治理政策的深化發(fā)展階段,農(nóng)業(yè)面源污染治理政策實(shí)施數(shù)量快速增長。(20)華春林、張玖弘、金書泰:《基于文本量化的中國農(nóng)業(yè)面源污染治理政策演進(jìn)特征分析》,《中國農(nóng)業(yè)科學(xué)》2022年第7期?；诖?采用2011—2020年的數(shù)據(jù)進(jìn)行穩(wěn)健性檢驗(yàn)。更換權(quán)重矩陣和縮短時(shí)間窗口后,農(nóng)業(yè)FDI數(shù)量和質(zhì)量的回歸系數(shù)和空間系數(shù)仍均顯著為正。綜合來看,穩(wěn)健性檢驗(yàn)結(jié)果與前文基準(zhǔn)回歸結(jié)果基本一致,說明本文的估計(jì)結(jié)果是穩(wěn)健的。

六、主要結(jié)論與政策建議

本文主要研究結(jié)論如下:第一,我國農(nóng)業(yè)面源污染具有顯著的空間正自相關(guān)性,農(nóng)業(yè)環(huán)境受到本地和鄰地污染排放的雙重影響。農(nóng)業(yè)面源污染排放在我國地理分布上形成一定的聚集區(qū)域,以高—高聚集和低—低聚集為主要特征。第二,現(xiàn)階段農(nóng)業(yè)FDI的流入會(huì)在一定程度上加劇我國環(huán)境污染。FDI數(shù)量和質(zhì)量的提高均會(huì)在顯著增加本地農(nóng)業(yè)面源污染排放的同時(shí),增加鄰地污染排放。FDI數(shù)量對(duì)面源污染的負(fù)面影響比FDI質(zhì)量的影響更大。第三,經(jīng)濟(jì)規(guī)模、人口密度、技術(shù)水平、環(huán)境規(guī)制和農(nóng)業(yè)結(jié)構(gòu)等因素均會(huì)影響農(nóng)業(yè)面源污染。

結(jié)合上述研究,得到以下三點(diǎn)政策啟示:第一,優(yōu)化外資結(jié)構(gòu),提高引資門檻。在提升對(duì)外開放水平和積極引進(jìn)外資的同時(shí),應(yīng)根據(jù)本地區(qū)的環(huán)境承載能力和經(jīng)濟(jì)發(fā)展水平設(shè)置一定的引資門檻。限制高污染農(nóng)業(yè)FDI,鼓勵(lì)高質(zhì)量和清潔技術(shù)密集型FDI的進(jìn)入,倒逼農(nóng)業(yè)生產(chǎn)者通過技術(shù)創(chuàng)新減少污染成本,實(shí)現(xiàn)由“污染天堂”向“污染光環(huán)”效應(yīng)的轉(zhuǎn)變。第二,加強(qiáng)各地區(qū)間的引資合作和治污政策合作。我國農(nóng)業(yè)面源污染具有顯著的空間相關(guān)性。因此,在制定相關(guān)政策時(shí),不僅要關(guān)注各地區(qū)內(nèi)部因素的影響,也要考慮周邊地區(qū)因素的交互影響。盡量降低政策的負(fù)外部性,避免省份間農(nóng)業(yè)面源污染的此消彼長?？紤]到環(huán)境規(guī)制強(qiáng)度的空間溢出效應(yīng)不顯著,應(yīng)著重保障各地農(nóng)業(yè)政策體系的順暢銜接,打破阻礙省域之間融合發(fā)展的行政壁壘,實(shí)現(xiàn)農(nóng)業(yè)面源污染的聯(lián)防聯(lián)控。第三,加快綠色生產(chǎn)方式轉(zhuǎn)型。我國農(nóng)業(yè)尚處在高增長、高污染的發(fā)展階段,技術(shù)水平提升對(duì)農(nóng)業(yè)環(huán)境帶來的正面效應(yīng)不顯著。因此,應(yīng)鼓勵(lì)技術(shù)創(chuàng)新,切實(shí)提高農(nóng)業(yè)生產(chǎn)技術(shù)中的綠色含量,增加農(nóng)產(chǎn)品科技投入,真正發(fā)揮綠色技術(shù)和技術(shù)溢出效應(yīng)對(duì)我國環(huán)境質(zhì)量的改善作用。積極學(xué)習(xí)運(yùn)用“千村示范、萬村整治”工程經(jīng)驗(yàn),對(duì)省域產(chǎn)業(yè)結(jié)構(gòu)進(jìn)行合理布局,積極引導(dǎo)地區(qū)間產(chǎn)業(yè)示范性帶來的空間效應(yīng),盡快實(shí)現(xiàn)農(nóng)業(yè)的集約型內(nèi)涵式發(fā)展。