張鵬靜，熊 濤,2*，朱 侯

（1.華中農(nóng)業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院，湖北 武漢 430070；2.華中農(nóng)業(yè)大學(xué) 信息與價(jià)格預(yù)測(cè)運(yùn)籌中心，湖北 武漢 430070；3.中山大學(xué) 信息管理學(xué)院，廣東 廣州 510006）

民為國(guó)基，谷為民命，糧食安全是國(guó)家安全的重要基礎(chǔ)。而保障糧食安全首先要保護(hù)耕地。為此，黨和政府出臺(tái)了一系列措施保護(hù)18億畝耕地紅線不受損，2022年中央一號(hào)文件提出要落實(shí)“長(zhǎng)牙齒”的耕地保護(hù)硬措施，實(shí)行耕地保護(hù)黨政同責(zé)，嚴(yán)守18 億畝耕地紅線；2023 年中央一號(hào)文件進(jìn)一步強(qiáng)調(diào)需加強(qiáng)耕地保護(hù)和用途管控，嚴(yán)格耕地占補(bǔ)平衡管理，實(shí)行部門(mén)聯(lián)合開(kāi)展補(bǔ)充耕地驗(yàn)收評(píng)定和“市縣審核、省級(jí)復(fù)核、社會(huì)監(jiān)督”機(jī)制，確保補(bǔ)充的耕地?cái)?shù)量相等、質(zhì)量相當(dāng)、產(chǎn)能不降。但是，耕地拋荒現(xiàn)象仍層出不窮，中國(guó)農(nóng)村家庭追蹤調(diào)查數(shù)據(jù)發(fā)現(xiàn)，2017 和2019 年約11.8%的小農(nóng)戶有拋荒耕地的行為[1]，遙感大數(shù)據(jù)監(jiān)測(cè)1992-2015年中國(guó)耕地拋荒面積達(dá)到了559170.26km2，占總耕地面積的18.59%[2]。耕地拋荒增加了糧食生產(chǎn)壓力，對(duì)糧食安全造成威脅。

改革開(kāi)放以來(lái)，中國(guó)雖然屢創(chuàng)經(jīng)濟(jì)奇跡，但是經(jīng)濟(jì)粗放型發(fā)展帶來(lái)的環(huán)境污染問(wèn)題已經(jīng)成為影響中國(guó)可持續(xù)發(fā)展的巨大障礙[3]。2022 年339 個(gè)地級(jí)以上城市中，空氣質(zhì)量超標(biāo)占比為37.2%，平均超標(biāo)天數(shù)比例為13.5%；468 個(gè)監(jiān)測(cè)降水的城市（區(qū)、縣）中出現(xiàn)酸雨的城市比例為33.8%；監(jiān)測(cè)地表水的3629個(gè)國(guó)考斷面中，Ⅳ類以下共占比12.1%；1890個(gè)國(guó)家地下水環(huán)境質(zhì)量考核的點(diǎn)位中，Ⅴ類占比22.4%；質(zhì)量為高等地的耕地面積比例僅為31.24%①數(shù)據(jù)來(lái)源于中華人民共和國(guó)生態(tài)環(huán)境部《2022中國(guó)生態(tài)環(huán)境狀況公報(bào)》。。這使我們逐漸意識(shí)到良好的生態(tài)環(huán)境是實(shí)現(xiàn)永續(xù)發(fā)展的基礎(chǔ)。

習(xí)近平總書(shū)記在黨的二十大報(bào)告中強(qiáng)調(diào)：“必須牢固樹(shù)立和踐行綠水青山就是金山銀山的理念，站在人與自然和諧共生的高度謀劃發(fā)展”。這意味著保護(hù)生態(tài)環(huán)境就是保護(hù)生產(chǎn)力，改善生態(tài)環(huán)境就是發(fā)展生產(chǎn)力。相反，如果生態(tài)環(huán)境受損，無(wú)疑會(huì)給社會(huì)經(jīng)濟(jì)發(fā)展帶來(lái)嚴(yán)重的負(fù)面影響，在農(nóng)業(yè)生產(chǎn)方面尤為突出。2014-2022 年臭氧污染導(dǎo)致中國(guó)冬小麥、單季稻、早稻和晚稻的產(chǎn)量分別下降了9.0%～17.0%、6.0%～8.2%、4.8%～9.4%和3.9%～7.0%，經(jīng)濟(jì)損失總計(jì)749 億元/年～1719 億元/年[4]；土壤銅污染會(huì)抑制水稻生長(zhǎng)，使單位面積穗數(shù)和每穗穎花數(shù)下降而使水稻減產(chǎn)[5]。除此之外，環(huán)境污染還會(huì)影響人體健康，2017 年中國(guó)因大氣污染、水污染造成的健康損失分別為9632.92 億元和681.04億元，損失總價(jià)值為當(dāng)年GDP的1.218%[6]，通過(guò)損耗勞動(dòng)力健康資本，降低勞動(dòng)力供給和勞動(dòng)生產(chǎn)率[7]，從而影響農(nóng)業(yè)生產(chǎn)。

對(duì)于農(nóng)業(yè)生產(chǎn)而言，可觀的產(chǎn)出回報(bào)是農(nóng)戶保持生產(chǎn)動(dòng)力的關(guān)鍵，充足的勞動(dòng)力供應(yīng)是進(jìn)行農(nóng)業(yè)生產(chǎn)的基礎(chǔ)。空氣、水和土壤是作物生長(zhǎng)不可或缺的要素，清潔的空氣和干凈的水也是人體健康的保證，空氣污染、水污染和土壤污染作為環(huán)境污染的重要組成部分，通過(guò)干擾作物生長(zhǎng)和人體健康，最終會(huì)影響農(nóng)戶家庭的農(nóng)業(yè)生產(chǎn)決策。但是已有文獻(xiàn)分別探討了環(huán)境污染對(duì)作物單產(chǎn)[8，9]、居民健康[10]以及勞動(dòng)力供給[11]等方面的影響，而在單產(chǎn)降低導(dǎo)致農(nóng)業(yè)生產(chǎn)收益減少以及健康受損制約勞動(dòng)力供給時(shí)，農(nóng)戶是否會(huì)將能力之外以及成本收益較低的耕地拋荒，無(wú)論在理論還是實(shí)證上都很少有研究對(duì)此進(jìn)行探究，本文力圖彌補(bǔ)該研究不足。故本文將耕地拋荒作為研究對(duì)象，利用中國(guó)勞動(dòng)力動(dòng)態(tài)調(diào)查（CLDS）2014、2016和2018年三期面板數(shù)據(jù)和地級(jí)市環(huán)境污染數(shù)據(jù)，在空氣、水和土壤視域下評(píng)估環(huán)境污染對(duì)耕地拋荒的影響效應(yīng)。此外，對(duì)環(huán)境污染與耕地拋荒關(guān)系的探討更是從實(shí)踐上證明“保護(hù)生態(tài)環(huán)境就是保護(hù)生產(chǎn)力”。

相較于已有研究，本文的邊際貢獻(xiàn)主要為:第一，在現(xiàn)有文獻(xiàn)研究環(huán)境污染對(duì)作物單產(chǎn)和人體健康的基礎(chǔ)上，將耕地拋荒納入環(huán)境污染影響評(píng)估的研究范疇，實(shí)證分析環(huán)境污染對(duì)耕地拋荒的影響；第二，在實(shí)證檢驗(yàn)環(huán)境污染與耕地拋荒的關(guān)系之前，構(gòu)建農(nóng)戶耕地利用理論模型，從理論上推導(dǎo)環(huán)境污染對(duì)耕地拋荒的影響。

一、文獻(xiàn)綜述

耕地拋荒可以認(rèn)為是在一系列因素的綜合作用下，農(nóng)戶在一定時(shí)期內(nèi)對(duì)可耕種的耕地停止種植或者減少耕種,從而使耕地處于一種荒蕪或未充分利用的狀態(tài)[12]。空氣、水和土壤是作物和人類賴以生存的基礎(chǔ)，而環(huán)境污染會(huì)干擾作物生長(zhǎng)影響作物單產(chǎn)、損害人體健康制約勞動(dòng)力供給，所以從理論上來(lái)說(shuō)，環(huán)境污染會(huì)改變農(nóng)戶的農(nóng)業(yè)生產(chǎn)決策，對(duì)耕地拋荒造成一定的影響。

作物單產(chǎn)是衡量耕地生產(chǎn)力的重要指標(biāo)，也是農(nóng)戶計(jì)算成本收益率的重要參考因素。作為理性的“經(jīng)濟(jì)人”，當(dāng)作物收益難以覆蓋種植成本時(shí)，拋荒就成為了農(nóng)戶家庭的最優(yōu)生產(chǎn)策略。所以，單產(chǎn)越低的耕地，越有可能被拋荒[13]。但是，以單產(chǎn)為切入點(diǎn)研究耕地拋荒問(wèn)題時(shí)，以往研究主要關(guān)注地塊特征，如坡度高、土壤質(zhì)量差的耕地由于生產(chǎn)力表現(xiàn)不佳，所以被拋荒的概率較高[14]，而很少考慮地塊以外的因素。自然科學(xué)領(lǐng)域的研究成果表明，作物生長(zhǎng)是與自然界交換能量的過(guò)程，環(huán)境污染會(huì)直接作用于作物的生物物理性狀，導(dǎo)致光合作用等其他生理活動(dòng)受損，干擾能量交換，從而造成作物產(chǎn)量下降[8-9]。例如，臭氧使中國(guó)小麥、水稻和玉米的單產(chǎn)分別損失了33%、23%和9%[15]；土壤鉛污染使單株小麥籽粒產(chǎn)量減少7.0%～26.4%[16]；2008年河北省邢臺(tái)市先于村因灌溉污水導(dǎo)致40公頃玉米減產(chǎn)[17]。當(dāng)作物單產(chǎn)因環(huán)境污染影響而降低時(shí)，農(nóng)戶會(huì)重新考慮耕地的種植決策，對(duì)于產(chǎn)量收益低于種植成本的地塊，停止種植、拋荒這些耕地顯然是合理的。所以，環(huán)境污染會(huì)通過(guò)降低作物單產(chǎn)影響農(nóng)戶拋荒耕地。

勞動(dòng)力是進(jìn)行農(nóng)業(yè)生產(chǎn)不可或缺的要素，充足的勞動(dòng)力供應(yīng)是實(shí)現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展的基礎(chǔ)。隨著中國(guó)經(jīng)濟(jì)高速增長(zhǎng)，農(nóng)業(yè)產(chǎn)業(yè)與非農(nóng)產(chǎn)業(yè)的比較收益差距逐漸拉大，為追求高收入，大量農(nóng)村青壯年勞動(dòng)力由農(nóng)村轉(zhuǎn)移進(jìn)城市，從農(nóng)業(yè)產(chǎn)業(yè)轉(zhuǎn)向非農(nóng)產(chǎn)業(yè)，部分農(nóng)村地區(qū)開(kāi)始出現(xiàn)勞動(dòng)力短缺現(xiàn)象[18]。以往研究在探討影響耕地拋荒的因素時(shí)，大多認(rèn)為農(nóng)業(yè)勞動(dòng)力非農(nóng)轉(zhuǎn)移造成農(nóng)業(yè)生產(chǎn)缺乏必需的勞動(dòng)力是導(dǎo)致耕地拋荒的重要原因[19-21]。但與此同時(shí)，也有研究表明健康水平越低的農(nóng)戶家庭越傾向于拋荒耕地[22]，原因在于勞動(dòng)力的健康水平?jīng)Q定了農(nóng)戶農(nóng)業(yè)生產(chǎn)時(shí)的勞動(dòng)力供給量，健康狀態(tài)越差，可提供的勞動(dòng)力供給就越少，農(nóng)戶只能將能力之外的耕地拋荒。仔細(xì)研究可以發(fā)現(xiàn)，勞動(dòng)力供給不足是導(dǎo)致農(nóng)戶耕地拋荒的重要因素，而勞動(dòng)力數(shù)量減少和勞動(dòng)力質(zhì)量降低是體現(xiàn)勞動(dòng)力供給不足的兩個(gè)方面。農(nóng)業(yè)勞動(dòng)力非農(nóng)轉(zhuǎn)移致使從事農(nóng)業(yè)生產(chǎn)的勞動(dòng)力數(shù)量大幅減少，勞動(dòng)力的健康水平限制勞動(dòng)參與率，最終都表現(xiàn)為因勞動(dòng)力供給受限而拋荒耕地。但是大多研究主要集中于勞動(dòng)力數(shù)量方面研究耕地拋荒，通過(guò)勞動(dòng)力質(zhì)量研究耕地拋荒的較少。近些年來(lái)，環(huán)境污染與人體健康的話題備受關(guān)注，作為一種典型的負(fù)外部效應(yīng)，環(huán)境污染的出現(xiàn)顯著增加了人類的健康風(fēng)險(xiǎn)[23-25]，特別是在一些低收入、醫(yī)療條件落后的地區(qū)，環(huán)境污染威脅人類健康的問(wèn)題更加突出[26]。對(duì)于醫(yī)療基礎(chǔ)設(shè)施差和經(jīng)濟(jì)發(fā)展水平低的農(nóng)村地區(qū)來(lái)說(shuō)[27]，農(nóng)村人口的健康水平更易受到環(huán)境污染的侵?jǐn)_。當(dāng)家庭農(nóng)業(yè)勞動(dòng)力健康水平因受環(huán)境污染影響而降低時(shí)，迫于體力限制，農(nóng)戶不得不將部分耕地拋荒。所以，環(huán)境污染會(huì)通過(guò)損害農(nóng)業(yè)勞動(dòng)力健康影響農(nóng)戶拋荒耕地。

環(huán)境污染通過(guò)降低作物單產(chǎn)、損害人體健康影響農(nóng)戶拋荒耕地，但值得注意的是農(nóng)戶在生產(chǎn)特征和個(gè)人特征等方面具有極強(qiáng)的差異，從而表現(xiàn)為環(huán)境污染對(duì)耕地拋荒的影響在不同農(nóng)戶群體中呈現(xiàn)異質(zhì)性。第一，糧食作物和經(jīng)濟(jì)作物的收益差距較大，種植經(jīng)濟(jì)作物的利潤(rùn)率更高。雖然環(huán)境污染對(duì)糧食作物和經(jīng)濟(jì)作物的單產(chǎn)均會(huì)產(chǎn)生顯著的負(fù)向影響，但是經(jīng)濟(jì)作物的價(jià)格更高，收入更多[28]。出現(xiàn)環(huán)境污染后，種植糧食作物會(huì)先跨過(guò)成本收益平衡點(diǎn)，所以，理論上來(lái)說(shuō)主營(yíng)糧食作物的農(nóng)戶家庭的耕地拋荒程度可能更高；第二，老齡化將成為中國(guó)未來(lái)一段時(shí)期的基本國(guó)情，且農(nóng)村地區(qū)的人口老齡化情況遠(yuǎn)比城鎮(zhèn)地區(qū)嚴(yán)重[29]。隨著年齡增長(zhǎng)，人體機(jī)能退化，老年人對(duì)周圍環(huán)境的質(zhì)量變化更敏感，暴露于污染環(huán)境中時(shí)的健康風(fēng)險(xiǎn)更高[30]。所以相對(duì)于年輕人，老年人更易受到環(huán)境污染的影響，其耕地拋荒程度可能更高，在人口老齡化程度不斷提高的農(nóng)村地區(qū)，環(huán)境污染對(duì)耕地拋荒的影響可能會(huì)愈加嚴(yán)峻。

二、理論分析

參考Cui[31]的研究，本文從作物單產(chǎn)和健康風(fēng)險(xiǎn)兩條影響路徑構(gòu)建農(nóng)戶耕地利用理論模型，以剖析環(huán)境污染對(duì)耕地拋荒的影響。假設(shè)農(nóng)戶在作物種植前做出生產(chǎn)決策，選擇是否拋荒耕地以及拋荒哪些耕地以使其利潤(rùn)最大化。令A(yù)為農(nóng)戶拋荒面積，L為耕地總面積，則拋荒程度為a=A L。假定農(nóng)產(chǎn)品銷售價(jià)格是外生給定的，農(nóng)戶是價(jià)格接受者，給定農(nóng)產(chǎn)品價(jià)格為p。作物產(chǎn)量Q同時(shí)取決于拋荒程度a、環(huán)境污染E、勞動(dòng)力投入農(nóng)業(yè)生產(chǎn)時(shí)間T（代替健康風(fēng)險(xiǎn)變量，健康風(fēng)險(xiǎn)越高，投入農(nóng)業(yè)生產(chǎn)的時(shí)間就越少）以及其他決定產(chǎn)量的因素α（包括農(nóng)業(yè)生產(chǎn)資料等）。農(nóng)業(yè)生產(chǎn)成本由邊際成本c和固定成本B組成。

根據(jù)Prishchepov 等[32]和Li等[13]的研究，當(dāng)農(nóng)戶拋荒耕地時(shí)，將優(yōu)先拋荒單產(chǎn)較低的耕地，而后才拋荒單產(chǎn)較高的耕地。因此，假設(shè)作物產(chǎn)量會(huì)隨著拋荒程度的增加而減少，并且減少的速度將遞增，即假設(shè)環(huán)境污染E是一個(gè)對(duì)產(chǎn)量有著負(fù)面影響的單維變量，即環(huán)境污染將降低邊際拋荒耕地產(chǎn)量的減少速度，即當(dāng)且僅當(dāng)種植一單位作物的收入大于等于其種植成本時(shí)，耕地才不會(huì)被拋荒[33]，實(shí)現(xiàn)利潤(rùn)最大化時(shí)滿足環(huán)境污染對(duì)耕地拋荒的影可以通過(guò)對(duì)一階條件兩邊的E進(jìn)行微分來(lái)得出：

由以上假設(shè)可以得出環(huán)境污染對(duì)耕地拋荒具有正向影響：

三、數(shù)據(jù)來(lái)源、變量說(shuō)明與實(shí)證模型

1.數(shù)據(jù)來(lái)源

文章所用數(shù)據(jù)來(lái)源于中山大學(xué)社會(huì)科學(xué)調(diào)查中心2014、2016 和2018 年進(jìn)行的“中國(guó)勞動(dòng)力動(dòng)態(tài)調(diào)查（CLDS）”。該數(shù)據(jù)采用多階段、多層次、與勞動(dòng)力規(guī)模成比例的抽樣方法，涵蓋了29 個(gè)省、自治區(qū)、直轄市（不包括西藏和海南），由社區(qū)模塊、家庭模塊和個(gè)人模塊三部分組成，包括了個(gè)人成長(zhǎng)經(jīng)歷、家庭經(jīng)濟(jì)狀況以及社區(qū)發(fā)展等內(nèi)容，樣本具有較強(qiáng)的代表性和適用性，首先將家庭模塊和社區(qū)模塊合并為一個(gè)數(shù)據(jù)集，然后再合并成面板數(shù)據(jù)。在數(shù)據(jù)分析過(guò)程中，僅保留了農(nóng)村戶口的樣本，去除了關(guān)鍵變量為缺失和異常的觀測(cè)值，最后獲得了覆蓋中國(guó)26個(gè)省的17265戶樣本。

2.變量描述

根據(jù)以往文獻(xiàn)構(gòu)建變量的方法，結(jié)合本文研究?jī)?nèi)容的需要，主要使用了以下變量來(lái)檢驗(yàn)環(huán)境污染與耕地拋荒之間的關(guān)系。

（1）被解釋變量。本文的被解釋變量為拋荒程度。以農(nóng)戶拋荒面積與其耕地面積的比值指代拋荒程度。

（2）核心解釋變量。本文的核心解釋變量為環(huán)境污染，主要從空氣污染、水污染和土壤污染3 個(gè)方面構(gòu)造環(huán)境污染變量。中國(guó)的環(huán)境污染主要來(lái)源于工業(yè)污染[34-35]，工業(yè)三廢（廢氣、廢水和廢渣）的排放是造成空氣污染、水污染和土壤污染的重要原因[36-38]。但由于地級(jí)市層面的工業(yè)固體廢棄物（廢渣）排放量數(shù)據(jù)缺失，本文參照張可等[39]的研究，利用工業(yè)二氧化硫排放量、工業(yè)煙（粉）塵排放量和工業(yè)廢水排放量使用熵值法構(gòu)建環(huán)境污染指數(shù)，數(shù)據(jù)來(lái)源于《中國(guó)城市統(tǒng)計(jì)年鑒》。農(nóng)戶在做耕地利用決策時(shí)一般會(huì)根據(jù)之前的生產(chǎn)經(jīng)營(yíng)狀況做出判斷，所以核心解釋變量使用t-1 期的環(huán)境污染指數(shù)，環(huán)境污染指數(shù)越大代表當(dāng)?shù)丨h(huán)境污染情況越嚴(yán)重。

（3）控制變量。為了排除其他混雜因素對(duì)估計(jì)結(jié)果的影響，參考Chen等[40]的研究，控制了可能影響耕地拋荒的戶主特征、家庭特征和村莊特征。其中戶主特征為戶主年齡；家庭特征為家庭債務(wù)、醫(yī)療支出占比、農(nóng)業(yè)收入、農(nóng)業(yè)收入占比、糧食作物補(bǔ)貼、耕地面積、農(nóng)業(yè)生產(chǎn)機(jī)械化程度和土地確權(quán)等；村莊特征為村莊常住勞動(dòng)力占比、村莊是否提供統(tǒng)一購(gòu)買生產(chǎn)資料服務(wù)、村莊是否提供外出務(wù)工服務(wù)、村莊行政區(qū)域內(nèi)是否有信用合作社、村莊是否處于大中等城市郊區(qū)以及作物生長(zhǎng)期是否需要灌溉等。變量的描述性統(tǒng)計(jì)見(jiàn)表1。

表1 變量定義及描述性統(tǒng)計(jì)

3.計(jì)量模型設(shè)定

本文使用雙向固定效應(yīng)模型估計(jì)環(huán)境污染對(duì)耕地拋荒的影響，計(jì)量模型如下：

式（4）中，下標(biāo)i=1,…,n表示農(nóng)戶個(gè)體，t=2014,2016,2018 表示觀察年份,Abandonment_ratioit表示農(nóng)戶i在第t年的耕地拋荒程度，Envi_pollutionit-1表示農(nóng)戶i所在地級(jí)市第t-1 年的環(huán)境污染指數(shù)，Xit為控制變量，μi為農(nóng)戶個(gè)體固定效應(yīng)，λt為時(shí)間固定效應(yīng)，εit為回歸殘差項(xiàng)，β0為截距項(xiàng)，β1和β2為回歸系數(shù)。

四、實(shí)證結(jié)果

1.基準(zhǔn)回歸

運(yùn)用Stata 軟件對(duì)式（4）進(jìn)行回歸。在基準(zhǔn)回歸之前，對(duì)變量之間可能存在的多重共線性問(wèn)題進(jìn)行了檢驗(yàn)，結(jié)果顯示方差膨脹因子（VIF）的取值范圍為1～2.36，遠(yuǎn)小于10，說(shuō)明變量間不存在嚴(yán)重的多重共線性問(wèn)題。表2 報(bào)告了環(huán)境污染對(duì)耕地拋荒影響的基準(zhǔn)回歸結(jié)果。表2 列（1）沒(méi)有加入時(shí)間固定效應(yīng)，環(huán)境污染指數(shù)變量在5%的水平上顯著，且系數(shù)為正，說(shuō)明環(huán)境污染對(duì)耕地拋荒產(chǎn)生了顯著的正向影響。表2 列（2）加入時(shí)間固定效應(yīng)，環(huán)境污染指數(shù)變量仍在5%的水平上顯著，且系數(shù)為正，表明環(huán)境污染越嚴(yán)重，農(nóng)戶的耕地拋荒程度就會(huì)越高?；鶞?zhǔn)回歸結(jié)果表明，環(huán)境污染會(huì)顯著正向影響農(nóng)戶拋荒耕地，與第2 節(jié)理論分析結(jié)果一致。環(huán)境污染的負(fù)外部性會(huì)對(duì)農(nóng)業(yè)生產(chǎn)造成一定的負(fù)面影響，導(dǎo)致農(nóng)戶家庭改變耕地利用決策，即在勞動(dòng)力承受范圍內(nèi)保證作物種植有利可得的情況下拋荒部分耕地，隨著環(huán)境污染程度的加深，農(nóng)戶會(huì)拋荒更多的耕地。另外，在控制變量中發(fā)現(xiàn)，農(nóng)業(yè)收入和糧食作物補(bǔ)貼變量均在1%的水平上顯著為負(fù)，農(nóng)業(yè)收入越高，獲得糧食作物補(bǔ)貼越多的農(nóng)戶家庭的耕地拋荒程度越低，與實(shí)踐一致。

表2 環(huán)境污染對(duì)耕地拋荒的影響：基準(zhǔn)結(jié)果N=17265

除此之外，為了比較不同類型環(huán)境污染對(duì)耕地拋荒的影響，考慮到變量的單位不同，利用極差標(biāo)準(zhǔn)化方法將工業(yè)廢水排放量、工業(yè)二氧化硫排放量和工業(yè)煙（粉）塵排放量進(jìn)行標(biāo)準(zhǔn)化處理，然后再放入模型回歸，回歸結(jié)果見(jiàn)表3。結(jié)果顯示工業(yè)二氧化硫排放量和工業(yè)煙（粉）塵排放量對(duì)耕地拋荒產(chǎn)生了顯著的正向影響，而工業(yè)廢水排放量對(duì)耕地拋荒的影響系數(shù)為正，但不顯著?？梢钥闯霾煌愋铜h(huán)境污染對(duì)耕地拋荒的影響存在差異，氣體污染物對(duì)耕地拋荒的影響更明顯，可能的原因是氣體的擴(kuò)散速度較快，污染范圍更廣，可以在短時(shí)間內(nèi)干擾作物生長(zhǎng)和人體健康，從而對(duì)耕地拋荒產(chǎn)生了顯著影響。所以在處理環(huán)境污染影響農(nóng)業(yè)生產(chǎn)的問(wèn)題時(shí)，可以優(yōu)先考慮氣體污染物的治理，此類污染物的及時(shí)處理可以在較大程度上減輕環(huán)境污染對(duì)農(nóng)業(yè)生產(chǎn)的負(fù)面影響，從而減少農(nóng)戶拋荒耕地。

表3 不同類型環(huán)境污染對(duì)耕地拋荒的影響N=17265

2.穩(wěn)健性檢驗(yàn)

為檢驗(yàn)基準(zhǔn)回歸結(jié)果的穩(wěn)健性，文章分別使用替換核心解釋變量、被解釋變量、增加核心解釋變量滯后期、提高聚類層次以及使用平衡面板數(shù)據(jù)的方法估計(jì)環(huán)境污染對(duì)耕地拋荒程度的影響。

第一，文章的環(huán)境污染指數(shù)利用的是《中國(guó)城市統(tǒng)計(jì)年鑒》中工業(yè)廢水排放量、工業(yè)二氧化硫排放量和工業(yè)煙（粉）塵排放量三個(gè)指標(biāo)進(jìn)行構(gòu)建的，為了更好地反映農(nóng)村地區(qū)的生態(tài)環(huán)境，更加準(zhǔn)確地估計(jì)環(huán)境污染對(duì)農(nóng)戶耕地拋荒的影響，使用村級(jí)問(wèn)卷中“當(dāng)年是否發(fā)生過(guò)因?yàn)榄h(huán)境污染問(wèn)題與政府或相關(guān)企業(yè)交涉的事件？”問(wèn)題的回答作為環(huán)境污染的代理變量替換環(huán)境污染指數(shù)，如果村莊當(dāng)年的確因?yàn)榄h(huán)境污染問(wèn)題與政府或相關(guān)企業(yè)交涉過(guò)，那么說(shuō)明當(dāng)?shù)卮_實(shí)存在環(huán)境污染。表4列（1）為估計(jì)結(jié)果，可以看出環(huán)境污染變量在10%水平上顯著，且系數(shù)為正，表明環(huán)境污染會(huì)導(dǎo)致農(nóng)戶增加拋荒程度。同時(shí)，考慮到農(nóng)業(yè)生產(chǎn)對(duì)特定環(huán)境污染類型的敏感性，文章還利用村莊問(wèn)卷中對(duì)于當(dāng)?shù)厣鷳B(tài)環(huán)境問(wèn)題更加細(xì)化的回答作為核心解釋變量進(jìn)行估計(jì)，具體為當(dāng)年是否因?yàn)榭諝馕廴?、土壤污染或水污染與政府或相關(guān)企業(yè)交涉，如果村莊曾因上述三種污染中的至少一種與政府或相關(guān)企業(yè)交涉過(guò)，則認(rèn)為當(dāng)?shù)卮嬖诃h(huán)境污染。估計(jì)結(jié)果見(jiàn)表4列（2），可以發(fā)現(xiàn)環(huán)境污染變量仍在10%的水平上顯著，且系數(shù)為正，同樣說(shuō)明環(huán)境污染會(huì)顯著正向影響農(nóng)戶拋荒耕地。

表4 穩(wěn)健性檢驗(yàn)結(jié)果

第二，使用耕地拋荒面積替換拋荒程度進(jìn)行穩(wěn)健性檢驗(yàn)，表4 列（3）為估計(jì)結(jié)果，結(jié)果顯示環(huán)境污染指數(shù)變量在10%的水平上顯著，且系數(shù)為正，說(shuō)明環(huán)境污染越嚴(yán)重農(nóng)戶的耕地拋荒面積越多。

第三，雖然文章使用農(nóng)戶農(nóng)業(yè)生產(chǎn)前一年的環(huán)境污染指數(shù)變量進(jìn)行估計(jì)，但是考慮到環(huán)境污染對(duì)耕地拋荒的影響可能具有一定的時(shí)間延續(xù)性，即耕地拋荒可能不是由t-1 期的環(huán)境污染引起的，而是由t-1 期之前的環(huán)境污染導(dǎo)致的，為此文章進(jìn)一步加入環(huán)境污染指數(shù)變量的滯后項(xiàng)檢驗(yàn)基準(zhǔn)回歸結(jié)果的穩(wěn)健性。根據(jù)AIC 和BIC 兩個(gè)信息準(zhǔn)則確定最佳滯后階數(shù)，當(dāng)滯后2 階時(shí)AIC 值和BIC 值最小，因此將t-2 期環(huán)境污染指數(shù)和t-3 期環(huán)境污染指數(shù)納入式（4）進(jìn)行估計(jì)，估計(jì)結(jié)果見(jiàn)表4列（4）?？梢钥闯鰐-1期環(huán)境污染指數(shù)變量在1%的水平上顯著，且系數(shù)為正，同樣說(shuō)明環(huán)境污染會(huì)顯著正向影響農(nóng)戶拋荒耕地，與基準(zhǔn)回歸結(jié)果一致。

第四，考慮到同一村莊內(nèi)農(nóng)戶之間的隨機(jī)擾動(dòng)項(xiàng)可能相關(guān)，將聚類層次從農(nóng)戶提升至村莊，表4列（5）估計(jì)結(jié)果顯示環(huán)境污染指數(shù)變量在10%的水平上顯著，且系數(shù)為正，同樣說(shuō)明環(huán)境污染越嚴(yán)重，農(nóng)戶的耕地程度就會(huì)越高。

第五，非平衡面板的最大問(wèn)題在于，如果一些樣本不再出現(xiàn)的原因是內(nèi)生的（即與擾動(dòng)項(xiàng)相關(guān)），則會(huì)導(dǎo)致所使用的樣本不再具有代表性，為此，文章僅保留2014 年、2016 年和2018 年均存在的樣本構(gòu)建平衡面板檢驗(yàn)環(huán)境污染對(duì)耕地拋荒的影響，表4 列（6）估計(jì)結(jié)果顯示，環(huán)境污染指數(shù)變量在1%的水平上顯著，且系數(shù)為正，說(shuō)明環(huán)境污染會(huì)顯著正向影響農(nóng)戶拋荒耕地。以上，說(shuō)明本文的基準(zhǔn)回歸結(jié)果具有一定的穩(wěn)健性。

3.工具變量法

在考慮農(nóng)戶個(gè)人、家庭和村莊層面協(xié)變量的基礎(chǔ)上，進(jìn)一步利用工具變量法解決可能因遺漏變量和反向因果關(guān)系造成的內(nèi)生性問(wèn)題。近些年來(lái)，為了提高環(huán)境質(zhì)量，政府開(kāi)始實(shí)施環(huán)境污染治理投資，文章使用《中國(guó)環(huán)境統(tǒng)計(jì)年鑒》2012、2014 和2016 年，即t-2 期的省級(jí)工業(yè)污染源治理投資額作為工具變量，顯然t-2 期工業(yè)污染源治理投資額與t-1 期環(huán)境污染指數(shù)直接相關(guān)，但在第t期農(nóng)戶耕地拋荒決策中外生，能夠滿足工具變量有效性的兩個(gè)條件。

表5 匯報(bào)了內(nèi)生性檢驗(yàn)結(jié)果。Kleibergen-Paap rk LM統(tǒng)計(jì)量的P值小于0.001，拒絕了工具變量識(shí)別不足的假設(shè)；Kleibergen-Paap rk Wald F統(tǒng)計(jì)值為13.787，大于10，拒絕了工具變量為弱工具變量的假設(shè)。第一階段的回歸結(jié)果表明工具變量系數(shù)在1%水平上顯著為負(fù)，說(shuō)明t-2期工業(yè)污染源治理投資額越高，t-1期的環(huán)境污染指數(shù)越低。第二階段估計(jì)結(jié)果顯示，環(huán)境污染指數(shù)變量在1%水平上顯著為正，表明環(huán)境污染越嚴(yán)重，農(nóng)戶的耕地拋荒程度就會(huì)越高。說(shuō)明文章的研究結(jié)論是可靠的。

表5 內(nèi)生性檢驗(yàn)結(jié)果 N=17265

4.機(jī)制檢驗(yàn)

環(huán)境污染不僅會(huì)直接影響作物生長(zhǎng)，致使單產(chǎn)降低，還會(huì)增加勞動(dòng)力健康風(fēng)險(xiǎn)，最終導(dǎo)致農(nóng)戶拋荒耕地?；诖耍恼聵?gòu)造作物單產(chǎn)和勞動(dòng)力健康風(fēng)險(xiǎn)變量檢驗(yàn)環(huán)境污染影響耕地拋荒的微觀機(jī)制。首先，考慮到不同作物的差異，文章選擇村莊第一或第二主產(chǎn)糧食作物是玉米、小麥或水稻的農(nóng)戶家庭為研究對(duì)象，將玉米單產(chǎn)、小麥單產(chǎn)和水稻單產(chǎn)作為作物單產(chǎn)的代理變量。參考前人的研究[22]，根據(jù)家庭農(nóng)業(yè)勞動(dòng)力的自評(píng)健康狀況，將選擇“非常健康”“健康”“一般”的農(nóng)業(yè)勞動(dòng)力認(rèn)定為“健康”，并賦值為1，將選擇“比較不健康”“非常不健康”的農(nóng)業(yè)勞動(dòng)力認(rèn)定為“不健康”，并賦值為0，然后計(jì)算家庭農(nóng)業(yè)勞動(dòng)力的健康比例，家庭農(nóng)業(yè)勞動(dòng)力健康比例=健康家庭農(nóng)業(yè)勞動(dòng)力數(shù)量/家庭農(nóng)業(yè)勞動(dòng)力總數(shù)量。

在具體的估計(jì)過(guò)程中，考慮到環(huán)境污染對(duì)作物單產(chǎn)和農(nóng)業(yè)勞動(dòng)力健康的影響可能存在一定的時(shí)間滯后性，文章將t-2 期環(huán)境污染指數(shù)和t-3 期環(huán)境污染指數(shù)變量同時(shí)加入了模型之中，估計(jì)結(jié)果見(jiàn)表6。列（1）、（3）、（5）、（7）沒(méi)有加入環(huán)境污染指數(shù)變量的滯后項(xiàng)，環(huán)境污染指數(shù)變量分別在1%和10%的水平上顯著，且系數(shù)為負(fù)，表明環(huán)境污染會(huì)顯著降低作物單產(chǎn)和農(nóng)業(yè)勞動(dòng)力健康。列（2）、（4）、（6）、（8）加入t-2 期環(huán)境污染指數(shù)和t-3 期環(huán)境污染指數(shù)，t-1 期環(huán)境污染指數(shù)變量仍對(duì)作物單產(chǎn)和農(nóng)業(yè)勞動(dòng)力健康產(chǎn)生了顯著的負(fù)向影響，無(wú)論是否加入環(huán)境污染指數(shù)變量的滯后項(xiàng)，t-1期環(huán)境污染指數(shù)變量對(duì)作物單產(chǎn)和農(nóng)業(yè)勞動(dòng)力健康均產(chǎn)生了顯著的負(fù)向影響，由此可知，環(huán)境污染會(huì)通過(guò)降低作物單產(chǎn)和增加農(nóng)業(yè)勞動(dòng)力健康風(fēng)險(xiǎn)兩條機(jī)制路徑導(dǎo)致農(nóng)戶拋荒耕地。

表6 機(jī)制分析結(jié)果

5.異質(zhì)性分析

（1）單產(chǎn)雖然是影響作物成本收益率的重要因素，但價(jià)格的高低也決定了作物種植利潤(rùn)的多寡。與經(jīng)濟(jì)作物相比，糧食作物的價(jià)格相對(duì)較低[28]，在環(huán)境污染同時(shí)干擾經(jīng)濟(jì)作物和糧食作物時(shí)，糧食作物的成本收益率會(huì)首先因單產(chǎn)的降低而越過(guò)成本收益平衡點(diǎn)，所以理論上來(lái)說(shuō)種植糧食作物的農(nóng)戶更容易受到環(huán)境污染的影響而拋荒耕地。根據(jù)村莊農(nóng)業(yè)收入主要來(lái)源是否來(lái)自于糧食生產(chǎn)，將農(nóng)戶家庭劃分為糧食作物組和非糧食作物組，分析環(huán)境污染在種植不同類型作物的農(nóng)戶群體中對(duì)耕地拋荒的異質(zhì)性影響。圖1 估計(jì)結(jié)果顯示環(huán)境污染指數(shù)變量?jī)H對(duì)以糧食生產(chǎn)為農(nóng)業(yè)收入主要來(lái)源的農(nóng)戶家庭的耕地拋荒程度產(chǎn)生了顯著的正向影響。這說(shuō)明環(huán)境污染對(duì)糧食作物種植戶的耕地拋荒影響更強(qiáng)烈，環(huán)境污染會(huì)阻礙糧食作物種植的可持續(xù)發(fā)展。所以在治理環(huán)境污染造成的耕地拋荒問(wèn)題時(shí)，應(yīng)重點(diǎn)關(guān)注糧食作物種植比例較大的地區(qū)，如糧食主產(chǎn)區(qū)。

圖1 環(huán)境污染對(duì)耕地拋荒的異質(zhì)性影響

（2）老齡化是世界各國(guó)普遍的發(fā)展趨勢(shì)，中國(guó)的年輕新一代接受高等教育后返回農(nóng)村務(wù)農(nóng)的數(shù)量很少，均以走進(jìn)“城市”為目標(biāo)，家中留守的老人逐漸成為了農(nóng)業(yè)生產(chǎn)經(jīng)營(yíng)的主要力量。樣本中2014 年戶主平均年齡為54.158 歲，2016 年戶主平均年齡為55.479 歲，2018 年戶主平均年齡為57.341歲，農(nóng)業(yè)生產(chǎn)經(jīng)營(yíng)逐漸呈現(xiàn)老齡化趨勢(shì)。隨著年齡的增長(zhǎng)，人體機(jī)能逐漸衰退，老年人在污染環(huán)境中的健康風(fēng)險(xiǎn)更高[28]。所以在環(huán)境污染的情況下，老年農(nóng)戶的耕地拋荒程度可能會(huì)更高。按樣本戶主平均年齡55.63歲進(jìn)行劃分，探討環(huán)境污染在不同年齡組對(duì)耕地拋荒的異質(zhì)性影響。圖1估計(jì)結(jié)果顯示環(huán)境污染指數(shù)在年齡高于均值組中對(duì)耕地拋荒產(chǎn)生了顯著的正向影響，在年齡低于均值組中并不顯著①95%水平的置信區(qū)間下，年齡高于均值組中環(huán)境污染指數(shù)變量系數(shù)的顯著性水平為10%，其置信區(qū)間包含0，為了在圖中清晰地表示年齡高于均值組中環(huán)境污染指數(shù)變量的顯著性，在按年齡分組的回歸模型中將置信區(qū)間改為90%，其他模型保持95%水平的置信區(qū)間不變。，這說(shuō)明老年農(nóng)戶容易受到環(huán)境污染的影響而拋荒耕地。老齡化是中國(guó)未來(lái)一段時(shí)期的基本國(guó)情，農(nóng)村人口的老齡化趨勢(shì)更加嚴(yán)峻，如果不及時(shí)對(duì)環(huán)境污染進(jìn)行治理，將會(huì)導(dǎo)致更多的耕地拋荒，嚴(yán)重阻礙農(nóng)業(yè)可持續(xù)發(fā)展。

五、結(jié)論與政策啟示

耕地作為農(nóng)業(yè)生產(chǎn)的基礎(chǔ)，是保證中國(guó)人民能把“飯碗”牢牢端在自己手上的首要條件。中國(guó)人口占全球總?cè)丝诘?0%，但只擁有全球7%的耕地，耕地資源十分匱乏。而在此情形下，中國(guó)耕地拋荒現(xiàn)象仍層出不窮，如何緩解耕地拋荒趨勢(shì)，保護(hù)好18億畝耕地紅線，對(duì)保障中國(guó)糧食安全具有重大意義。文章使用中國(guó)勞動(dòng)力動(dòng)態(tài)調(diào)查(CLDS) 2014、2016 和2018 年三期農(nóng)戶家庭數(shù)據(jù)，匹配2013、2015 和2017 年地級(jí)市環(huán)境污染指數(shù)（利用《中國(guó)城市統(tǒng)計(jì)年鑒》中工業(yè)廢水排放量、工業(yè)二氧化硫排放量和工業(yè)煙（粉）塵排放量數(shù)據(jù)構(gòu)造），構(gòu)建農(nóng)戶耕地利用理論模型，使用雙向固定效應(yīng)模型實(shí)證檢驗(yàn)了環(huán)境污染對(duì)耕地拋荒的影響和機(jī)制。主要研究結(jié)論為：（1）環(huán)境污染對(duì)耕地拋荒具有顯著的正向影響，替換核心解釋變量、被解釋變量、增加核心解釋變量滯后期、提高聚類層次、使用平衡面板數(shù)據(jù)以及工具變量法后，研究結(jié)論依然成立；（2）機(jī)制分析發(fā)現(xiàn)環(huán)境污染不僅顯著降低了作物單產(chǎn)，同時(shí)也顯著增加了農(nóng)業(yè)勞動(dòng)力的健康風(fēng)險(xiǎn)，繼而導(dǎo)致農(nóng)戶拋荒耕地；（3）異質(zhì)性分析顯示，環(huán)境污染對(duì)糧食作物種植戶、老年農(nóng)戶的耕地拋荒產(chǎn)生了顯著影響，中國(guó)農(nóng)業(yè)生產(chǎn)主要群體更易受到環(huán)境污染的影響。習(xí)近平總書(shū)記指出，綠水青山就是金山銀山，保護(hù)生態(tài)環(huán)境就是保護(hù)生產(chǎn)力，改善生態(tài)環(huán)境就是發(fā)展生產(chǎn)力。本文的研究發(fā)現(xiàn)從實(shí)踐上證明了這一結(jié)論，環(huán)境污染對(duì)農(nóng)業(yè)生產(chǎn)的負(fù)面影響導(dǎo)致農(nóng)戶拋荒耕地，而耕地是實(shí)現(xiàn)農(nóng)業(yè)可持續(xù)發(fā)展的基礎(chǔ)，所以良好的生態(tài)環(huán)境是保證生產(chǎn)力永續(xù)的前提條件，相反，生態(tài)環(huán)境質(zhì)量的降低無(wú)疑會(huì)損害生產(chǎn)力的發(fā)展。

通過(guò)上述研究結(jié)論得出如下政策啟示：

第一，政府部門(mén)在評(píng)估環(huán)境污染造成的糧食生產(chǎn)損失時(shí)，通常僅關(guān)注環(huán)境污染導(dǎo)致的糧食單產(chǎn)下降。本文研究發(fā)現(xiàn)，環(huán)境污染會(huì)導(dǎo)致農(nóng)戶拋荒耕地，使得該地塊的糧食產(chǎn)量為零。因此，相比于環(huán)境污染誘發(fā)的糧食單產(chǎn)降低，環(huán)境污染導(dǎo)致的耕地拋荒對(duì)糧食生產(chǎn)的影響可能更大。學(xué)界和政府部門(mén)在評(píng)估環(huán)境污染對(duì)糧食生產(chǎn)的影響時(shí)，不僅要考慮環(huán)境污染對(duì)糧食單產(chǎn)的影響，還要明確環(huán)境污染對(duì)耕地拋荒的影響。

第二，政府部門(mén)在分配環(huán)境污染稅收收入時(shí)，需量化環(huán)境污染對(duì)農(nóng)業(yè)生產(chǎn)的負(fù)面影響，對(duì)處于環(huán)境污染區(qū)的農(nóng)業(yè)生產(chǎn)主體進(jìn)行相應(yīng)的補(bǔ)償，彌補(bǔ)環(huán)境污染造成的農(nóng)業(yè)生產(chǎn)損失，保證農(nóng)戶獲得穩(wěn)定的農(nóng)業(yè)收入，繼而減少耕地拋荒，穩(wěn)定糧食生產(chǎn)。同時(shí)，杜絕環(huán)境污染的產(chǎn)生，并治理已有的環(huán)境污染是解決環(huán)境污染阻礙農(nóng)業(yè)可持續(xù)發(fā)展問(wèn)題的關(guān)鍵。氣體污染物因其快速傳播的特征對(duì)耕地拋荒的影響更顯著，因此在污染治理實(shí)踐中應(yīng)優(yōu)先考慮氣體污染物。

第三，政府部門(mén)在解決環(huán)境污染問(wèn)題時(shí)，應(yīng)該優(yōu)先關(guān)注老齡化程度較高的地區(qū)和糧食作物種植比例較大的地區(qū)，重點(diǎn)幫扶老齡農(nóng)戶和糧食種植戶，確保此類農(nóng)戶能夠從糧食生產(chǎn)中獲得長(zhǎng)期穩(wěn)定的收益，緩解環(huán)境污染造成的耕地拋荒，最終保障我國(guó)糧食安全。