張永強(qiáng) 牛嘉宸

東北農(nóng)業(yè)大學(xué)，黑龍江 哈爾濱 150030

一、問題的提出

我國糧食供求長期處于“緊平衡”狀態(tài)，人口規(guī)模加劇糧食供給壓力，土地資源和水資源等要素約束趨緊，化肥成為保障糧食供給安全的重要手段?；实倪^量投入不僅對環(huán)境產(chǎn)生負(fù)外部性[1]，也阻礙了我國農(nóng)業(yè)可持續(xù)發(fā)展進(jìn)程[2]。自2014年以來，中央一號(hào)文件連續(xù)8年對農(nóng)業(yè)投入品過量使用引起農(nóng)村生態(tài)環(huán)境突出問題予以關(guān)注，2015年國家出臺(tái)《到2020年化肥使用量零增長行動(dòng)方案》，強(qiáng)調(diào)要加大農(nóng)業(yè)面源污染防治力度，大力推進(jìn)化肥減量提效，減少不合理投入，實(shí)現(xiàn)化肥農(nóng)藥使用量零增長、負(fù)增長。

農(nóng)戶是農(nóng)業(yè)生產(chǎn)的基礎(chǔ)單元和直接主體，其生產(chǎn)決策行為是影響農(nóng)業(yè)生產(chǎn)環(huán)境質(zhì)量的重要因素，也是我國農(nóng)村環(huán)境問題治理的源頭。青壯年農(nóng)戶“非農(nóng)化”流失使得農(nóng)業(yè)勞動(dòng)力在要素稟賦和知識(shí)技能等方面出現(xiàn)明顯的結(jié)構(gòu)性短板。李旻、趙連閣等研究表明，勞動(dòng)力老齡化、女性化在體力、文化程度和有效勞動(dòng)供給時(shí)間等方面存在的劣勢會(huì)對農(nóng)業(yè)生產(chǎn)效率產(chǎn)生約束效應(yīng)[3-4]。鄭宏運(yùn)等指出，理性農(nóng)戶為實(shí)現(xiàn)生產(chǎn)收益最大化，通常會(huì)選擇增施化肥以彌補(bǔ)勞動(dòng)力結(jié)構(gòu)性缺失對農(nóng)業(yè)生產(chǎn)的不利影響，但化肥等生產(chǎn)要素的長期錯(cuò)配會(huì)降低農(nóng)業(yè)生態(tài)效率，從而對生產(chǎn)要素投入規(guī)模產(chǎn)生負(fù)反饋效應(yīng)，并進(jìn)一步使農(nóng)戶追加化肥等要素投入形成惡性循環(huán)[5]。當(dāng)化學(xué)投入品用量超過環(huán)境承載力時(shí)，農(nóng)戶最終會(huì)在“理性”目標(biāo)下引起非理性行為而造成環(huán)境污染[6]。胡雪枝、林本喜等研究則認(rèn)為生產(chǎn)技術(shù)和組織方式的改進(jìn)降低了農(nóng)業(yè)勞動(dòng)強(qiáng)度和對農(nóng)戶的體力需求，使農(nóng)業(yè)生產(chǎn)逐漸成為標(biāo)準(zhǔn)化活動(dòng)，因而農(nóng)業(yè)勞動(dòng)力結(jié)構(gòu)轉(zhuǎn)變不會(huì)對農(nóng)業(yè)生產(chǎn)造成顯著影響[7-8]。田紅宇等研究發(fā)現(xiàn)，土地流轉(zhuǎn)和農(nóng)業(yè)社會(huì)化服務(wù)的調(diào)節(jié)作用可以緩解勞動(dòng)力質(zhì)量下滑對農(nóng)業(yè)生產(chǎn)技術(shù)效率的不利影響[9]，使農(nóng)戶減輕對化肥等要素的依賴。那么，農(nóng)業(yè)勞動(dòng)力的結(jié)構(gòu)失衡是否會(huì)引發(fā)農(nóng)戶化肥施用行為的改變并對化肥面源污染強(qiáng)度產(chǎn)生影響？本文分析我國化肥面源污染強(qiáng)度的時(shí)空格局和動(dòng)態(tài)演進(jìn)趨勢，構(gòu)建STIRPAT模型，探究勞動(dòng)力結(jié)構(gòu)轉(zhuǎn)變與化肥面源污染強(qiáng)度的關(guān)系，以檢驗(yàn)城鎮(zhèn)化和勞動(dòng)力非農(nóng)轉(zhuǎn)移背景下勞動(dòng)力結(jié)構(gòu)變化在農(nóng)業(yè)生產(chǎn)中引起的環(huán)境效應(yīng)。

二、數(shù)據(jù)來源及處理

本文選擇2004—2018年全國31個(gè)?。ㄗ灾螀^(qū)、直轄市）的面板數(shù)據(jù)作為研究樣本。研究所需的化肥施用折純量、糧食作物播種面積、農(nóng)作物總播種面積、農(nóng)業(yè)機(jī)械總動(dòng)力等數(shù)據(jù)來自于《中國農(nóng)村統(tǒng)計(jì)年鑒》。農(nóng)業(yè)總產(chǎn)值、農(nóng)林牧漁業(yè)總產(chǎn)值、農(nóng)村居民人均可支配收入及各地區(qū)年末人口數(shù)等數(shù)據(jù)均來自《中國統(tǒng)計(jì)年鑒》。鄉(xiāng)村從業(yè)人口與農(nóng)業(yè)從業(yè)人口數(shù)據(jù)來自各?。ㄗ灾螀^(qū)、直轄市）的統(tǒng)計(jì)年鑒。測算老齡化、女性化和人力資本水平等核心解釋變量采用的數(shù)據(jù)均來自《中國人口與就業(yè)統(tǒng)計(jì)年鑒》。

勞動(dòng)力在農(nóng)業(yè)生產(chǎn)要素中處于支配性地位，其結(jié)構(gòu)特征與化肥施用水平密切相關(guān)。勞動(dòng)力結(jié)構(gòu)的概念體系豐富，包含年齡、性別、人力資本水平、勞動(dòng)部門、職業(yè)和經(jīng)濟(jì)形式等多個(gè)外延。鑒于本文研究的勞動(dòng)力主體是留守在農(nóng)村地區(qū)參與農(nóng)作物種植的農(nóng)戶，其勞動(dòng)部門、職業(yè)和經(jīng)濟(jì)形式相對統(tǒng)一，因此參考魏娟的研究方法，從年齡、性別和人力資本水平三方面考察勞動(dòng)力結(jié)構(gòu)特征[10]。其中，老齡化率（old）用于表示年齡結(jié)構(gòu)特征，為保證數(shù)據(jù)可獲得性，本文采用各地區(qū)鄉(xiāng)村65歲及以上人口數(shù)量占15歲及以上人口數(shù)量的百分比表示；女性化率（female）用于表示性別結(jié)構(gòu)特征，采用各地區(qū)鄉(xiāng)村15歲及以上人口中女性人口的比例表示；人力資本水平（edu），采用各地鄉(xiāng)村6歲及以上人口的平均受教育年限表示，其中小學(xué)、初中、高中、大專及以上水平的受教育年限分別按6、9、12、15.5年計(jì)算。

對于控制變量的選取，本文參考吳義根和欒健等的研究方法，選取人口密度（density）表示人口規(guī)模（P），人均農(nóng)業(yè)產(chǎn)值（agdp）表示地區(qū)經(jīng)濟(jì)水平（A），農(nóng)業(yè)種植結(jié)構(gòu)（ps）及化肥施用技術(shù)水平（tech）表示技術(shù)水平（T）作為控制變量。此外，將非農(nóng)轉(zhuǎn)移率（nlr）、農(nóng)民收入水平（income）、測土配方施肥政策（cp）和化肥零增長行動(dòng)政策（fp）等可能影響農(nóng)戶化肥施用決策的因素納入STIRPAT模型的考察范圍[11-12]。其中，人口密度（density）采用各地區(qū)單位耕地供給的人口數(shù)表示（即各地區(qū)年末人口數(shù)與各地區(qū)農(nóng)作物總播種面積的比值）。人均農(nóng)業(yè)產(chǎn)值（agdp），采用農(nóng)業(yè)GDP與各地區(qū)年末人口數(shù)的比值表示，其反映了農(nóng)業(yè)經(jīng)濟(jì)發(fā)達(dá)程度，對當(dāng)?shù)剞r(nóng)業(yè)生產(chǎn)經(jīng)營方式、管理水平、種植結(jié)構(gòu)產(chǎn)生直接影響。農(nóng)業(yè)種植結(jié)構(gòu)（ps），采用各地區(qū)糧食作物播種面積占農(nóng)作物總播種面積的比例表示。相比于糧食作物，經(jīng)濟(jì)作物的化肥利用效率較低，經(jīng)濟(jì)作物種植規(guī)模的擴(kuò)大會(huì)增加化肥投入數(shù)量，加重化肥面源污染?；适┯眉夹g(shù)水平（tech），采用各地區(qū)單位化肥折純量的有效產(chǎn)出表示，即各地區(qū)農(nóng)林牧漁總產(chǎn)值與化肥施用折純量之比。化肥施用技術(shù)水平的提高有利于提升化肥利用效率，使農(nóng)戶在投入更少化肥的前提下獲得更多產(chǎn)出。非農(nóng)轉(zhuǎn)移率（nlr），用各地區(qū)鄉(xiāng)村從業(yè)人口與農(nóng)業(yè)從業(yè)人口的差值占鄉(xiāng)村從業(yè)人口數(shù)量的比重表示。農(nóng)村勞動(dòng)力非農(nóng)轉(zhuǎn)移比例的擴(kuò)大，一方面提高了勞動(dòng)力要素的稀缺性，增強(qiáng)了化肥對勞動(dòng)力的替代效應(yīng)，引起化肥施用量增加，加劇化肥面源污染；另一方面促進(jìn)農(nóng)村土地流轉(zhuǎn)，通過提高土地經(jīng)營規(guī)?；潭?、化肥利用效率，減少化肥不合理投入。農(nóng)民收入水平（income），用各地區(qū)農(nóng)村居民人均可支配收入表示。收入上漲一方面會(huì)提高農(nóng)民對化肥等生產(chǎn)要素的購買能力，并通過追加化肥施用量以獲得更多農(nóng)業(yè)產(chǎn)出和經(jīng)營收益；另一方面會(huì)促使農(nóng)戶選擇質(zhì)量更高、環(huán)境危害更少的化學(xué)肥料，提高肥料利用效率的同時(shí)減輕面源污染。減污政策主要包含測土配方施肥政策（cp）和化肥零增長行動(dòng)政策（fp），以虛擬變量形式表示。其中，2005年中央一號(hào)文件提出要大力推廣測土配方施肥，因此各地區(qū)測土配方施肥政策（cp）的變量值在2005年之前取0，2005年以后取1。2015年我國制定了《到2020年化肥使用量零增長行動(dòng)方案》，減少化肥不合理投入，因此化肥零增長行動(dòng)政策（fp）的變量值在2015年之前取0，2015年以后取1。機(jī)械化的效率優(yōu)勢有利于彌補(bǔ)勞動(dòng)力體力和勞動(dòng)時(shí)間的不足，減輕農(nóng)戶對化肥的依賴。因此，本文引入機(jī)械化水平（power）作為虛擬變量，取各年度農(nóng)業(yè)機(jī)械總動(dòng)力中位數(shù)為衡量標(biāo)準(zhǔn)，各省份當(dāng)年農(nóng)業(yè)機(jī)械總動(dòng)力若大于中位數(shù)取1，反之取0，用于考察機(jī)械化在勞動(dòng)力結(jié)構(gòu)變化對化肥面源污染影響中的調(diào)節(jié)效應(yīng)。此外，考慮通貨膨脹可能對研究結(jié)果準(zhǔn)確性產(chǎn)生干擾，研究中涉及的經(jīng)濟(jì)指標(biāo)數(shù)據(jù)均換算為2004年的不變價(jià)格。

（一）化肥面源污染強(qiáng)度測算

化肥污染物主要通過土壤酸化、空氣揮發(fā)、水體污染等方式對生態(tài)環(huán)境產(chǎn)生負(fù)面影響。參考現(xiàn)有研究成果，本文重點(diǎn)考察化肥施用對水體的污染，采用單元調(diào)查評估法對全國31個(gè)省份的化肥面源污染強(qiáng)度進(jìn)行測算。由于化肥施用過程中，鉀肥不會(huì)直接產(chǎn)生面源污染，故本文選取化肥總氮（TN）和總磷（TP）污染物排放量作為核定化肥面源污染程度的指標(biāo)。污染物排放總量和強(qiáng)度的計(jì)算公式如下：

式中，E為氮、磷污染物的總排放量；EUi為不同種類化肥的施用量；ρi表示化肥的氮、磷產(chǎn)污系數(shù)，各種類化肥產(chǎn)污系數(shù)見表1。ηi為化肥流失系數(shù)，本文參照賴斯蕓等的地區(qū)劃分確定氮肥和磷肥流失系數(shù)，各地區(qū)化肥流失情況見表2[13-14]。EI表示化肥面源污染排放強(qiáng)度；AL表示農(nóng)作物總播種面積。

表1 不同種類化肥產(chǎn)污系數(shù)

表2 化肥流失情況分區(qū)統(tǒng)計(jì)表

（二）STIRPAT模型構(gòu)建

為研究經(jīng)濟(jì)發(fā)展對環(huán)境產(chǎn)生的壓力，Paul R.Ehrlich等在1971年提出IPAT等式，其基本形式為：

I=PAT

等式中，I為環(huán)境負(fù)荷，表示經(jīng)濟(jì)社會(huì)發(fā)展對環(huán)境產(chǎn)生的影響；P、A、T分別表示人口規(guī)模、經(jīng)濟(jì)水平和技術(shù)水平。但等式中其他條件維持恒定不變時(shí)，各影響因素對因變量的影響效應(yīng)是成等比例變化的，與一般經(jīng)濟(jì)運(yùn)行規(guī)律不相符。因此，Dietz等在2003年對IPAT恒等式進(jìn)行了完善，得到STIRPAT模型，其一般形式為：

其中，a為STIRPAT模型的系數(shù)，b、c、d分別為P（人口規(guī)模）、A（經(jīng)濟(jì)水平）和T（技術(shù)水平）等變量對環(huán)境影響的彈性，e為模型的隨機(jī)擾動(dòng)項(xiàng)。STIRPAT模型由于允許在已有的P（人口規(guī)模）、A（經(jīng)濟(jì)水平）和T（技術(shù)水平）等自變量基礎(chǔ)上引入其他可能對環(huán)境沖擊造成影響的驅(qū)動(dòng)因子作為解釋變量而得到廣泛應(yīng)用。可以看到，當(dāng)自變量的環(huán)境彈性、模型系數(shù)和隨機(jī)擾動(dòng)項(xiàng)取值均為1時(shí)，STIRPAT模型就轉(zhuǎn)化為原始的IPAT模型。本文基于勞動(dòng)力結(jié)構(gòu)變化對化肥面源污染的影響機(jī)理，在STIRPAT模型一般形式的基礎(chǔ)上引入其他可能影響化肥面源污染強(qiáng)度的解釋變量，在等式兩邊進(jìn)行對數(shù)化處理，將模型轉(zhuǎn)化為線性形式，構(gòu)建計(jì)量模型1：

式中，EIij表示第 i省第 j年的化肥面源污染強(qiáng)度，densityij、agdpij、psij、techij、nlrij、incomeij、oldij、femaleij、eduij、cpij、fpij分別表示第i省第j年的人口密度、人均農(nóng)業(yè)GDP、種植結(jié)構(gòu)、化肥施用技術(shù)水平、非農(nóng)轉(zhuǎn)移率、農(nóng)民收入水平、老齡化率、女性化率、人力資本水平、測土配方施肥政策和化肥零增長行動(dòng)政策情況。α0至α11為各變量的待估計(jì)系數(shù)，εij為隨機(jī)擾動(dòng)項(xiàng)。

為驗(yàn)證機(jī)械投入與勞動(dòng)力結(jié)構(gòu)變化對化肥面源污染的協(xié)同作用，在模型1的基礎(chǔ)上，利用機(jī)械總動(dòng)力（power）的虛擬變量構(gòu)建機(jī)械化水平及其與老齡化、女性化水平的交互項(xiàng)，得到計(jì)量模型2：

（三）化肥面源污染強(qiáng)度測算結(jié)果

本文基于單元分析的綜合調(diào)查評價(jià)方法和清單分析法，測算2004—2018年我國31個(gè)省份（自治區(qū)、直轄市）的化肥面源污染強(qiáng)度，結(jié)果見表3。

表3 2004—2018年各?。ㄗ灾螀^(qū)、直轄市）化肥面源污染強(qiáng)度

（四）化肥面源污染時(shí)空分布特征

本文根據(jù)化肥面源污染程度將各省（區(qū)、市）依次劃分為極重度污染（EI＞78 kg·hm-2）、重度污染（54 kg·hm-2＜EI≤78kg·hm-2）、中度污染（30 kg·hm-2＜EI≤54 kg·hm-2）和輕度污染（EI＜kg·hm-2）四類區(qū)域?；拭嬖次廴局囟纫陨现饕性诒本?、天津、江蘇、浙江、福建、廣東等地區(qū)，其中北京、天津、江蘇等地存在極重度污染情況。北京2017年化肥面源污染強(qiáng)度達(dá)到107.04 kg·hm-2，為統(tǒng)計(jì)期內(nèi)最高值。這些地區(qū)大多相互接壤，存在明顯的空間集聚特征及空間相關(guān)性。化肥面源污染中度地區(qū)主要集中在吉林、遼寧、河北、寧夏、陜西、山東、湖南、湖北、云南。其他地區(qū)整體保持在輕度污染狀態(tài)。2018年黑龍江化肥面源污染強(qiáng)度6.50 kg·hm-2，為統(tǒng)計(jì)期內(nèi)最低值。從時(shí)間演變趨勢看，2004—2008年我國化肥面源污染的重心逐漸由東北向西南地區(qū)轉(zhuǎn)移，2008—2016年保持不變。此后，隨著東北地區(qū)化肥面源污染強(qiáng)度逐漸減弱，化肥面源污染的重心再次向西南地區(qū)轉(zhuǎn)移。北京及東南沿海一帶仍是我國化肥面源污染嚴(yán)重地區(qū)，主要原因在于這些地區(qū)多為經(jīng)濟(jì)發(fā)達(dá)地區(qū)，經(jīng)濟(jì)的高速發(fā)展使得本地大量勞動(dòng)力向非農(nóng)產(chǎn)業(yè)轉(zhuǎn)移，糧食等農(nóng)產(chǎn)品供給能力不足。同時(shí)，外地勞動(dòng)力跨省流入現(xiàn)象明顯，糧食等農(nóng)產(chǎn)品作為剛需商品，供給壓力增大，化肥過量施用成為短時(shí)間內(nèi)提升農(nóng)產(chǎn)品供給能力的有效途徑。

為進(jìn)一步驗(yàn)證化肥面源污染存在空間依賴性，本文引入Moran’s I檢驗(yàn)用于表征鄰近省份化肥面源污染的空間相關(guān)性及差異性，檢驗(yàn)結(jié)果見表4。根據(jù)檢驗(yàn)結(jié)果可知，各年度化肥面源污染強(qiáng)度具有高度的空間相關(guān)性并通過顯著性檢驗(yàn)。各年度化肥面源污染強(qiáng)度的Moran指數(shù)均為正值，說明化肥面源污染存在空間集聚特征，而非隨機(jī)分布，鄰近地域的化肥面源污染存在關(guān)聯(lián)特征，這也更好地解釋了化肥面源污染重度以上地區(qū)呈現(xiàn)出相互毗鄰的原因。

表4 Moran’s I檢驗(yàn)結(jié)果

三、農(nóng)業(yè)勞動(dòng)力結(jié)構(gòu)變化與化肥面源污染關(guān)系的實(shí)證分析

（一）基于全國層面的實(shí)證分析

本文采用Stata軟件對各變量組成的面板數(shù)據(jù)進(jìn)行Hausman檢驗(yàn)，根據(jù)檢驗(yàn)結(jié)果，利用固定效應(yīng)模型對面板數(shù)據(jù)進(jìn)行回歸分析，回歸結(jié)果見表5。

表5 勞動(dòng)力結(jié)構(gòu)對化肥面源污染的影響回歸結(jié)果

分析結(jié)果發(fā)現(xiàn)，其一，老齡化率對化肥面源污染強(qiáng)度存在顯著的負(fù)向影響，且在1%的水平上通過顯著性檢驗(yàn)，說明隨著農(nóng)業(yè)勞動(dòng)力老齡化程度的提高，化肥面源污染強(qiáng)度逐漸降低。雖然老年勞動(dòng)力在體力上存在弱質(zhì)性，農(nóng)業(yè)生產(chǎn)受限，但隨著年齡的增加，老年勞動(dòng)力農(nóng)耕經(jīng)驗(yàn)積累效應(yīng)增強(qiáng)，有助于提升其生產(chǎn)效率，彌補(bǔ)體力上的不足，從而減少化學(xué)投入品的使用。其二，女性化率對化肥面源污染強(qiáng)度存在顯著的正向影響，且在10%的水平上通過顯著性檢驗(yàn)，說明農(nóng)業(yè)勞動(dòng)力中女性群體比例的擴(kuò)大，會(huì)加重化肥面源污染強(qiáng)度。隨著青壯年勞動(dòng)力非農(nóng)轉(zhuǎn)移的規(guī)模不斷擴(kuò)大，留守婦女逐漸成為農(nóng)村家庭從事農(nóng)業(yè)生產(chǎn)的重要力量，需要兼顧農(nóng)業(yè)生產(chǎn)和家庭責(zé)任?；瘜W(xué)投入品可以提高農(nóng)業(yè)生產(chǎn)效率，在保證產(chǎn)量的前提下縮短勞動(dòng)時(shí)間，為女性農(nóng)戶照顧老人、教育子女等提供了時(shí)間保障。其三，人力資本水平對化肥面源污染強(qiáng)度存在顯著的負(fù)向影響，且在5%的水平上通過顯著性檢驗(yàn)。本文對人力資本水平的測度主要考察農(nóng)業(yè)勞動(dòng)力的平均受教育年限，隨著受教育水平的提升，農(nóng)業(yè)勞動(dòng)力對現(xiàn)代化農(nóng)業(yè)生產(chǎn)技術(shù)和機(jī)械化設(shè)備的接受程度普遍提高，農(nóng)業(yè)新技術(shù)和機(jī)械設(shè)備的投用，可以從整體上促進(jìn)農(nóng)業(yè)生產(chǎn)效率提升，緩解農(nóng)戶對化學(xué)投入品的依賴，從而減輕化肥面源污染。

為考察機(jī)械化和勞動(dòng)力結(jié)構(gòu)對化肥面源污染影響的協(xié)同作用，本文在模型中引入機(jī)械化水平與老齡化、女性化指標(biāo)的交互項(xiàng)再次進(jìn)行回歸分析，結(jié)果見表6?；貧w分析結(jié)果顯示，在機(jī)械化生產(chǎn)的協(xié)同作用下，老齡化、女性化和人力資本水平指標(biāo)分別在1%、5%和10%的置信水平上顯著，且回歸系數(shù)正負(fù)效應(yīng)均未發(fā)生變化，但女性化率的系數(shù)數(shù)值遠(yuǎn)高于老齡化和人力資本水平，說明女性化率對化肥面源污染的促進(jìn)作用遠(yuǎn)大于老齡化和人力資本水平的抑制作用。因此，總體上，勞動(dòng)力結(jié)構(gòu)變化本身不利于減輕化肥面源污染。分析機(jī)械化與老齡化、女性化水平的交互項(xiàng)可以發(fā)現(xiàn)，老齡化與機(jī)械化交互項(xiàng)對化肥面源污染影響的回歸系數(shù)為0.054并在10%的水平上顯著，女性化與機(jī)械化交互項(xiàng)的回歸系數(shù)為-0.149并在10%的水平上通過顯著性檢驗(yàn)，說明農(nóng)業(yè)機(jī)械化投入增強(qiáng)了老齡化造成化肥面源污染的可能性，但對減輕勞動(dòng)力女性化引起的化肥面源污染具有顯著的積極調(diào)節(jié)效應(yīng)?？赡艿脑蚴抢淆g農(nóng)戶所經(jīng)營的土地通常表現(xiàn)出細(xì)碎分散的特征，機(jī)械化在農(nóng)業(yè)生產(chǎn)中的效率優(yōu)勢無法充分發(fā)揮，而且其對現(xiàn)代農(nóng)業(yè)技術(shù)和設(shè)備的學(xué)習(xí)和接受能力較低，機(jī)械投入以簡易型機(jī)械為主，操作熟練程度不高，機(jī)械化作業(yè)較為粗糙，無效率損耗較多，追加化肥投入可以用最小的成本彌補(bǔ)和降低這部分損失，但同時(shí)也加劇了化肥面源污染。而農(nóng)業(yè)機(jī)械的投入可以緩解女性農(nóng)村勞動(dòng)力因體力限制、家務(wù)負(fù)擔(dān)等原因加大化肥投入而帶來的化肥面源污染。整體而言，機(jī)械化投入對女性化引起化肥面源污染的負(fù)向調(diào)節(jié)效應(yīng)遠(yuǎn)高于對老齡化引起化肥面源污染的正向調(diào)節(jié)效應(yīng)，機(jī)械化的推廣與應(yīng)用有利于防治化肥面源污染。

表6 全國層面勞動(dòng)力結(jié)構(gòu)對化肥面源污染影響的回歸結(jié)果

（二）基于地區(qū)層面的實(shí)證分析

為進(jìn)一步分析中國不同地區(qū)勞動(dòng)力結(jié)構(gòu)對化肥面源污染的影響效應(yīng)，按照國家統(tǒng)計(jì)局的分類依據(jù)，將我國劃分為東、中、西和東北四個(gè)地區(qū)經(jīng)濟(jì)帶分別進(jìn)行回歸檢驗(yàn)，考察在機(jī)械化協(xié)同作用下，勞動(dòng)力結(jié)構(gòu)對化肥面源污染的影響效應(yīng)?；貧w結(jié)果見表7。

表7 地區(qū)層面勞動(dòng)力結(jié)構(gòu)對化肥面源污染影響的回歸結(jié)果

分析回歸結(jié)果可以發(fā)現(xiàn)，勞動(dòng)力結(jié)構(gòu)因素對化肥面源污染的影響效應(yīng)存在地區(qū)差異。中部、西部和東北地區(qū)老齡化率對化肥面源污染影響的回歸系數(shù)分別為-0.667、-0.081和-0.378且通過顯著性檢驗(yàn)，說明在我國中部、西部和東北地區(qū)，老齡化率對化肥面源污染的影響存在負(fù)向效應(yīng)，其作用程度在中部地區(qū)最強(qiáng)，東北地區(qū)次之，西部地區(qū)較弱。老齡化水平的提升對緩解化肥面源污染存在促進(jìn)作用，老齡農(nóng)戶憑借長期積累的農(nóng)業(yè)生產(chǎn)經(jīng)驗(yàn)可以熟練掌握化肥的施用時(shí)間、施用方式和施用數(shù)量，以更少的化肥投入獲得理想的農(nóng)業(yè)產(chǎn)出，從而減輕化肥面源污染。

女性化率在中部地區(qū)對化肥面源污染影響的回歸系數(shù)為1.977，并在1%的水平上通過顯著性檢驗(yàn)，說明農(nóng)業(yè)勞動(dòng)力女性比例的提升不利于中部地區(qū)化肥面源污染的治理。而女性化率在東北地區(qū)對化肥面源污染影響的回歸系數(shù)為-1.950，并在1%的水平上通過顯著性檢驗(yàn)，說明農(nóng)業(yè)勞動(dòng)力女性比例的提升有利于東北地區(qū)化肥面源污染的緩解。

人力資本水平對于化肥面源污染的作用表現(xiàn)出區(qū)域異質(zhì)性，東部地區(qū)人力資本水平對化肥面源污染影響的回歸系數(shù)為-0.219并在10%的水平上通過顯著性檢驗(yàn)，說明人力資本水平的提升有利于減輕其化肥面源污染。而中部、西部和東北地區(qū)人力資本水平對化肥面源污染影響的回歸系數(shù)為0.345、0.129和0.223并分別在5%、10%和10%的水平上通過顯著性檢驗(yàn)，說明在這些地區(qū)人力資本水平的提升明顯加劇了化肥面源污染。張淑輝等研究表明，人力資本水平與當(dāng)?shù)亟?jīng)濟(jì)水平存在密切聯(lián)系，不同等級的人力資本水平對綠色農(nóng)業(yè)生產(chǎn)的影響作用在各地區(qū)間存在差異[15]。東部地區(qū)經(jīng)濟(jì)發(fā)達(dá)，農(nóng)村人口有機(jī)會(huì)獲得質(zhì)量更好、時(shí)間更長的教育機(jī)會(huì)，為該地區(qū)高等人力資本的形成和積累創(chuàng)造了有利條件。高等人力資本水平的勞動(dòng)力在推動(dòng)?xùn)|部地區(qū)綠色農(nóng)業(yè)生產(chǎn)中發(fā)揮主要作用，其對現(xiàn)代農(nóng)業(yè)技術(shù)和化學(xué)投入品的應(yīng)用更為合理，有助于減輕化肥面源污染。而在中部和東北地區(qū)，初等和中等人力資本在推動(dòng)綠色農(nóng)業(yè)生產(chǎn)中作用顯著，義務(wù)教育階段的培養(yǎng)目標(biāo)與農(nóng)業(yè)生產(chǎn)知識(shí)和技能需求的脫節(jié)導(dǎo)致初等人力資本對綠色農(nóng)業(yè)生產(chǎn)的抑制效應(yīng)高于中等人力資本的促進(jìn)效應(yīng)。西部地區(qū)受整體經(jīng)濟(jì)水平的影響，具有較高人力資本水平的勞動(dòng)力傾向于轉(zhuǎn)移至勞動(dòng)報(bào)酬更高的非農(nóng)產(chǎn)業(yè)就業(yè)，農(nóng)村留守勞動(dòng)力整體教育水平較低，制約農(nóng)戶合理采用農(nóng)業(yè)技術(shù)成果和科學(xué)使用化肥農(nóng)藥，加大了化肥面源污染排放。

老齡化與機(jī)械化水平交互項(xiàng)在中部地區(qū)和東北地區(qū)的回歸系數(shù)分別為0.557和0.272并在1%的水平上通過了顯著性檢驗(yàn)，表明上述地區(qū)在機(jī)械化的調(diào)節(jié)作用下老齡化對化肥面源污染排放產(chǎn)生了推動(dòng)效應(yīng)，加重了化肥面源污染。女性化與機(jī)械化水平交互項(xiàng)在中部地區(qū)和東北地區(qū)的回歸系數(shù)分別為-1.604和-0.963并在1%的水平上通過了顯著性檢驗(yàn)，說明在上述地區(qū)農(nóng)業(yè)機(jī)械的投入可以緩解女性農(nóng)村勞動(dòng)力弱質(zhì)性引起的化肥面源污染?？傮w而言，勞動(dòng)力結(jié)構(gòu)變化在中部地區(qū)和東北地區(qū)對化肥面源污染防治的影響效應(yīng)最強(qiáng)，在東部地區(qū)的影響效應(yīng)最弱，農(nóng)業(yè)機(jī)械化投入在我國中部和東北地區(qū)具有明顯的調(diào)節(jié)作用。

四、結(jié)論與政策建議

（一）結(jié)論

本文利用2004—2018年31個(gè)省份的面板數(shù)據(jù)，分析了我國化肥面源污染的特征，通過構(gòu)建STIRPAT模型，從農(nóng)業(yè)勞動(dòng)力的性別結(jié)構(gòu)、年齡結(jié)構(gòu)和人力資本水平等方面考察近年來農(nóng)業(yè)勞動(dòng)力結(jié)構(gòu)轉(zhuǎn)變對化肥面源污染的影響效應(yīng)。研究結(jié)果表明：

老齡化水平和人力資本水平對化肥面源污染強(qiáng)度存在顯著的負(fù)向影響，且老齡化水平比人力資本水平的影響效應(yīng)更高。農(nóng)戶隨著年齡和受教育水平的提升，可以在一定程度上緩解化肥面源污染程度。女性化率對化肥面源污染強(qiáng)度存在顯著的正向影響，且影響效應(yīng)明顯高于老齡化水平和人力資本水平。農(nóng)業(yè)勞動(dòng)力女性比例擴(kuò)大，加劇化肥面源污染程度。農(nóng)業(yè)現(xiàn)代化及機(jī)械化生產(chǎn)的普及，有利于減輕女性勞動(dòng)力因彌補(bǔ)自身弱質(zhì)性追加化肥投入帶來的化肥面源污染，但老齡農(nóng)戶因彌補(bǔ)土地經(jīng)營規(guī)模細(xì)碎化、機(jī)械操作不熟練造成機(jī)械生產(chǎn)效率低、損耗高的不足，進(jìn)而增加化肥投入，存在加劇化肥面源污染的可能。勞動(dòng)力結(jié)構(gòu)變化對化肥面源污染的影響存在地區(qū)差異，勞動(dòng)力結(jié)構(gòu)變化對化肥面源污染防治的影響效應(yīng)在中部地區(qū)最強(qiáng)，西部地區(qū)最弱，農(nóng)業(yè)機(jī)械化投入在我國中部和東北地區(qū)調(diào)節(jié)作用明顯。

（二）政策建議

完善農(nóng)村公共服務(wù)和社會(huì)保障體系建設(shè)，減輕農(nóng)村老年和婦女群體生活負(fù)擔(dān)。從研究結(jié)果看，老齡農(nóng)戶由于非農(nóng)就業(yè)機(jī)會(huì)選擇減少，可以更加專注于農(nóng)業(yè)生產(chǎn)，相對農(nóng)業(yè)勞動(dòng)時(shí)間延長，化肥利用效率較高；女性勞動(dòng)力務(wù)農(nóng)時(shí)還要贍養(yǎng)老人、撫育子女等，通常需要通過追加化肥投入，彌補(bǔ)農(nóng)業(yè)勞動(dòng)時(shí)間不足。因此，政府要進(jìn)一步完善農(nóng)村地區(qū)公共服務(wù)和社會(huì)保障體系，提高鄉(xiāng)鎮(zhèn)衛(wèi)生院醫(yī)療服務(wù)質(zhì)量，為農(nóng)村老年居民從事農(nóng)業(yè)生產(chǎn)提供健康保障；建立農(nóng)村留守老人關(guān)愛服務(wù)體系，穩(wěn)定老年群體留鄉(xiāng)從事農(nóng)業(yè)生產(chǎn)意愿；引入社會(huì)資本，加強(qiáng)農(nóng)村托幼機(jī)構(gòu)建設(shè)并提高服務(wù)質(zhì)量，提升農(nóng)村中小學(xué)基礎(chǔ)設(shè)施及后勤服務(wù)水平，拓展義務(wù)教育學(xué)校課中、課后服務(wù)內(nèi)容，減輕女性勞動(dòng)力的家庭負(fù)擔(dān)。

創(chuàng)新優(yōu)化農(nóng)技推廣體系，針對老年和女性農(nóng)戶不同特點(diǎn)，提供差異化農(nóng)民技術(shù)培訓(xùn)服務(wù)。農(nóng)事經(jīng)驗(yàn)的積累及化肥施用技術(shù)水平的提升有助于減輕化肥面源污染。因此，在農(nóng)業(yè)勞動(dòng)力呈現(xiàn)老齡化和女性化的背景下，要面向老年和女性農(nóng)戶提供有針對性、差異化的農(nóng)業(yè)技術(shù)推廣服務(wù)。可針對老年勞動(dòng)力偏好，采用田間學(xué)校等開放型推廣模式，技術(shù)人員在田間一對一、手把手開展技術(shù)培訓(xùn)服務(wù)，鼓勵(lì)其積極參與農(nóng)技觀摩活動(dòng)。加強(qiáng)農(nóng)業(yè)技術(shù)應(yīng)用與推廣，為女性農(nóng)戶提供農(nóng)機(jī)社會(huì)化服務(wù)專項(xiàng)補(bǔ)貼，培育更多女性農(nóng)業(yè)生產(chǎn)精英。

加快培育新型職業(yè)農(nóng)民，拓展農(nóng)村義務(wù)教育體系，多層面提升農(nóng)村人力資本水平[16]。從短期看，提升農(nóng)戶農(nóng)業(yè)職業(yè)素養(yǎng)，有助于提升農(nóng)業(yè)生產(chǎn)效率，促進(jìn)化肥施用減量增效，因此要以提升農(nóng)業(yè)生產(chǎn)效率為導(dǎo)向，政府、市場和社會(huì)三方合力，構(gòu)建新型職業(yè)農(nóng)民培訓(xùn)體系，全方位提升農(nóng)戶職業(yè)素質(zhì)，重視對農(nóng)業(yè)生產(chǎn)者，特別是老齡農(nóng)戶開展農(nóng)業(yè)機(jī)械裝備的應(yīng)用技能培訓(xùn)，提升高層次新型職業(yè)農(nóng)民儲(chǔ)備。擴(kuò)大西部和東北地區(qū)新型職業(yè)農(nóng)民培訓(xùn)覆蓋規(guī)模和培訓(xùn)力度，為有志于長期從事農(nóng)業(yè)生產(chǎn)的農(nóng)戶提供進(jìn)修等多樣化培訓(xùn)形式，提升新型職業(yè)農(nóng)民人力資本水平。借助現(xiàn)代信息技術(shù)促進(jìn)優(yōu)質(zhì)教育資源在城鄉(xiāng)之間合理流動(dòng)，提高農(nóng)村義務(wù)教育學(xué)生的可及性。鼓勵(lì)我國中部地區(qū)的農(nóng)村中小學(xué)，在基礎(chǔ)教育階段根據(jù)地區(qū)農(nóng)業(yè)發(fā)展特點(diǎn)，將農(nóng)業(yè)生產(chǎn)基礎(chǔ)知識(shí)引入課堂，打造具有農(nóng)業(yè)特色的鄉(xiāng)土校本課程，全面提升初、中等人力資本對農(nóng)業(yè)生產(chǎn)知識(shí)的認(rèn)知水平，增強(qiáng)農(nóng)村義務(wù)教育對人力資本的“造血”功能。