(南京農(nóng)業(yè)大學(xué) 經(jīng)濟(jì)管理學(xué)院,江蘇 南京210095)

農(nóng)業(yè)是國民經(jīng)濟(jì)的基礎(chǔ),是國民經(jīng)濟(jì)其他部門賴以存在和發(fā)展的基礎(chǔ)。水稻是我國最重要的糧食作物之一,其種植總面積約占糧食作物種植面積的30%,產(chǎn)量接近糧食總產(chǎn)量的1/2,在農(nóng)業(yè)生產(chǎn)中具有很高的地位。2014年，中央經(jīng)濟(jì)工作會議指出我國經(jīng)濟(jì)發(fā)展已進(jìn)入新常態(tài),正從高速增長轉(zhuǎn)向中高速增長,經(jīng)濟(jì)增長方式從粗放型增長轉(zhuǎn)向集約型增長。我國“十三五”規(guī)劃明確將綠色發(fā)展納入五大發(fā)展理念,強(qiáng)調(diào)加快改善生態(tài)環(huán)境質(zhì)量。2018年“中央一號”文件提出加強(qiáng)農(nóng)業(yè)面源污染防治,開展農(nóng)業(yè)綠色發(fā)展行動,實(shí)現(xiàn)投入品減量化、有機(jī)肥替代化肥,推進(jìn)鄉(xiāng)村綠色發(fā)展,打造人與自然和諧共生發(fā)展新格局。本文以可持續(xù)發(fā)展理論為基礎(chǔ),研究水稻的綠色發(fā)展情況。

隨著我國城鎮(zhèn)化進(jìn)程的加快,農(nóng)村勞動力大量向城市轉(zhuǎn)移,造成了農(nóng)村人力資本不足。目前農(nóng)村老齡化的問題日趨嚴(yán)重,我國社會中的“人口紅利”逐漸消失。此外,國內(nèi)農(nóng)業(yè)生產(chǎn)成本快速攀升,各種農(nóng)業(yè)生產(chǎn)資料、土地資金、農(nóng)業(yè)勞動力價(jià)格不斷上漲。鑒于種種原因,只能依靠提高生產(chǎn)效率來轉(zhuǎn)變生產(chǎn)方式,走產(chǎn)出高效、產(chǎn)品安全、資源節(jié)約、環(huán)境友好的現(xiàn)代農(nóng)業(yè)發(fā)展道路。從20世紀(jì)80年代以來,隨著土地資源的稀缺,從事農(nóng)業(yè)的勞動力數(shù)量不斷下降,導(dǎo)致化肥投入急劇增加,我國迅速成為化肥生產(chǎn)量和消費(fèi)量最多的國家。然而,由于不合理施用化肥所帶來的面源污染問題、農(nóng)產(chǎn)品質(zhì)量安全問題,對農(nóng)業(yè)發(fā)展和農(nóng)民生活帶來了巨大沖擊。隨著農(nóng)業(yè)環(huán)境問題越來越嚴(yán)重,我國逐步認(rèn)識到不合理施用化肥對農(nóng)業(yè)可持續(xù)發(fā)展的影響越來越大。水資源是影響我國農(nóng)業(yè)發(fā)展的重要因素，高效節(jié)水灌溉有助于提高農(nóng)業(yè)生產(chǎn)效益,節(jié)約水資源,保障農(nóng)產(chǎn)品增產(chǎn)增收。然而,高效節(jié)水灌溉或選擇合理的灌溉方式是否有助于減少氮磷流失、提高化肥的利用率、減少化肥使用量,以減輕化肥污染、緩減農(nóng)業(yè)環(huán)境問題仍有待考察。

近年來,眾多學(xué)者對環(huán)境約束下的全要素生產(chǎn)率進(jìn)行了大量研究,并將其定義為綠色全要素生產(chǎn)率。陳詩一將CO2排放量作為投入要素測算了我國工業(yè)38個(gè)行業(yè)的環(huán)境效率,發(fā)現(xiàn)考慮環(huán)境約束的實(shí)際全要素生產(chǎn)率比傳統(tǒng)不考慮環(huán)境因素的估算值低[1];潘丹將我國30個(gè)省份分為中、東、西三個(gè)部分,測度了1998—2009年的農(nóng)業(yè)綠色全要素生產(chǎn)率,發(fā)現(xiàn)我國東部的綠色生產(chǎn)率水平高于中西部,中西部地區(qū)呈現(xiàn)出農(nóng)業(yè)經(jīng)濟(jì)發(fā)展緩慢和資源消耗、環(huán)境破壞的“雙重惡化”情況[2];吳傳清運(yùn)用SBM模型和GML指數(shù)測度了1997—2015年長江經(jīng)濟(jì)帶農(nóng)業(yè)綠色生產(chǎn)率,并對其影響因素作了實(shí)證分析,發(fā)現(xiàn)長江經(jīng)濟(jì)帶農(nóng)業(yè)綠色全要素生產(chǎn)率有所提升,但低于全國水平,機(jī)械化水平、人力資本存量和財(cái)政支出水平對長江經(jīng)濟(jì)帶農(nóng)業(yè)綠色全要素生產(chǎn)率具有顯著正向影響[3];展進(jìn)濤測度了考慮氮磷流失的農(nóng)業(yè)綠色生產(chǎn)率,并檢驗(yàn)了農(nóng)業(yè)綠色全要素生產(chǎn)率與糧食安全保障程度之間的關(guān)系,發(fā)現(xiàn)農(nóng)業(yè)綠色全要素生產(chǎn)率與糧食安全存在雙向因果關(guān)系[4]。

現(xiàn)有文獻(xiàn)雖然取得了一些豐碩的成果,但仍存在著不足,本文在此基礎(chǔ)上做了一些改進(jìn):①現(xiàn)有研究大多對農(nóng)業(yè)綠色全要素生產(chǎn)率、工業(yè)綠色全要素生產(chǎn)率等做了整體測度,極少有學(xué)者研究某一具體農(nóng)作物的綠色全要素生產(chǎn)率,本文將研究視角聚焦于水稻的綠色全要素生產(chǎn)率。②現(xiàn)有研究僅僅將灌溉設(shè)施水平作為其中一個(gè)變量進(jìn)行回歸檢驗(yàn),對灌溉設(shè)施水平與綠色全要素生產(chǎn)率的內(nèi)在關(guān)系并沒有進(jìn)行深入探討,本文著重分析了灌溉設(shè)施水平對水稻綠色全要素生產(chǎn)率的內(nèi)在影響。

鑒于此,本文首先對我國2004—2016年水稻化肥投入造成的環(huán)境成本進(jìn)行測算,基于我國18個(gè)省份的面板數(shù)據(jù)進(jìn)行實(shí)證分析,運(yùn)用ML指數(shù)測算我國不同主產(chǎn)區(qū)水稻綠色全要素生產(chǎn)率的變動與差異情況。其次，運(yùn)用面板OLS回歸模型探討灌溉設(shè)施水平對水稻綠色全要素生產(chǎn)率的影響。最后，根據(jù)研究結(jié)果提出基于綠色視角下促進(jìn)整個(gè)水稻產(chǎn)業(yè)生產(chǎn)效率提升的合理建議,提出必要的約束措施,落實(shí)可持續(xù)發(fā)展觀念。

1 研究方法與數(shù)據(jù)來源

1.1 研究方法

根據(jù)Chung[5]等學(xué)者的研究方法,引入動態(tài)效率Malmquist-Luenberger(ML)指數(shù),見式(1)。與傳統(tǒng)的Malmquist指數(shù)相比,ML指數(shù)可測算考慮環(huán)境污染引起的非期望產(chǎn)出下的全要素生產(chǎn)率。本文將化肥作為環(huán)境污染物,測算化肥環(huán)境成本下的水稻全要素生產(chǎn)率。

(1)

(2)

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

本文選取我國水稻主產(chǎn)區(qū)三大區(qū)域中的18個(gè)省份2004—2016年的統(tǒng)計(jì)數(shù)據(jù)進(jìn)行分析。根據(jù)《全國優(yōu)勢農(nóng)產(chǎn)品區(qū)域布局規(guī)劃(2008—2015年)》中的劃分方法,我國目前水稻主產(chǎn)區(qū)可分為三大區(qū)域:長江流域、東南沿海、東北平原。長江流域包括河南、江蘇、安徽、湖北、湖南、江西、重慶、四川、貴州、云南10個(gè)省市；東南沿海包括浙江、福建、廣東、廣西、海南5個(gè)省區(qū)；東北平原包括吉林、遼寧、黑龍江3個(gè)省。選取的投入—產(chǎn)出變量數(shù)據(jù)來自相關(guān)年份的《全國農(nóng)產(chǎn)品成本收益資料匯編》、《中國統(tǒng)計(jì)年鑒》和《中國農(nóng)村統(tǒng)計(jì)年鑒》等。

1.3 描述性分析

水稻種植的化肥投入施用量變化:首先是化肥施用結(jié)構(gòu)變化。從表1可知,2016年與2004年相比,氮肥總施用量減少了30.5%,磷肥比2004年減少了約3/4,鉀肥的施用量減少了約17.6%,但復(fù)混肥卻比2004年增加了約70.2%,導(dǎo)致化肥總施用量增長了約22.2%。因此,化肥的施用結(jié)構(gòu)雖然發(fā)生了變化,但化肥的總施用量在近年仍以年均1.7%的速度呈現(xiàn)上升的態(tài)勢。

表1 2004—2016年氮肥、磷肥、鉀肥、復(fù)混肥、化肥總施用量(萬t)

年份氮肥磷肥鉀肥復(fù)混肥化肥總和20041031.17202.59151.97391.351777.0820051009.88214.72179.70544.761949.062006960.41172.68172.28630.161935.532007969.36154.55178.16704.112006.182008955.41158.58130.96722.671967.622009963.33137.20119.87792.402012.802010933.71118.80124.06895.872072.442011912.16111.45121.88930.722076.212012866.98100.73120.57988.122076.402013819.9982.13119.231050.262071.612014847.8980.35131.901121.962182.092015822.5263.14129.631228.692243.972016789.9549.36129.221315.252283.78

注:數(shù)據(jù)來源于《全國農(nóng)產(chǎn)品成本收益資料匯編》，表2同。

我國三大主產(chǎn)區(qū)化肥施用強(qiáng)度變化:由表2可知,我國三大主產(chǎn)區(qū)中,東南沿海年均化肥施用強(qiáng)度最高,年均化肥施用強(qiáng)度為3.50kg/hm2,其次是長江流域,為3.06kg/hm2,東北平原最低,年均化肥施用強(qiáng)度為1.40kg/hm2。東南沿海屬于我國經(jīng)濟(jì)發(fā)達(dá)地區(qū),土地要素稀缺,農(nóng)民傾向于增加化肥的要素投入來提高產(chǎn)量；長江流域施用強(qiáng)度近年來逐漸升高,可能是由于長江流域稻區(qū)的病蟲害偏重發(fā)生,且發(fā)生范圍廣、程度重、影響大,因此化肥的施用量在近年日益增加；東北平原自然環(huán)境條件好、水資源充沛、土地資源豐富，化肥施用強(qiáng)度較低。

表2 2004—2016年中國三大主產(chǎn)區(qū)化肥施用量(kg/hm2)

水稻常用的灌溉方式:我國具有復(fù)雜的地貌和特殊的氣候,這些自然條件決定了我國必須走灌溉農(nóng)業(yè)的發(fā)展道路。水稻屬于習(xí)水農(nóng)作物,對水的需求非常高,其用水量約占農(nóng)業(yè)用水總量的70%[6]。發(fā)展水稻節(jié)水栽培,實(shí)現(xiàn)稻作可持續(xù)發(fā)展對我國糧食安全供應(yīng)和社會穩(wěn)定具有十分重要的意義[7]。節(jié)水灌溉技術(shù)多種多樣,目前在水稻生產(chǎn)上常用的灌溉技術(shù)主要有濕潤灌溉、干濕交替灌溉、水稻強(qiáng)化栽培體系、覆蓋旱種等。

1.4 指標(biāo)計(jì)算方法

投入變量:根據(jù)既有文獻(xiàn)和數(shù)據(jù)的可獲得性,本文選擇水稻播種面積、勞動用工量、總物質(zhì)與服務(wù)費(fèi)用三個(gè)投入指標(biāo)。水稻播種總面積比耕地面積更準(zhǔn)確,因?yàn)樗紤]了重復(fù)播種情況,反映土地投入狀況更為準(zhǔn)確[8]。為剔除價(jià)格因素的影響,本文使用各地區(qū)綜合生產(chǎn)資料價(jià)格指數(shù)對每公頃物質(zhì)與服務(wù)費(fèi)用進(jìn)行平減。由于勞動用工量、總物質(zhì)與服務(wù)費(fèi)用數(shù)據(jù)各省沒有單獨(dú)統(tǒng)計(jì),因此本研究采用如下的處理方法:勞動用工量按每公頃用工數(shù)量乘以水稻播種面積計(jì)算,物質(zhì)與服務(wù)費(fèi)用按每公頃物質(zhì)與服務(wù)費(fèi)用乘以水稻播種面積計(jì)算。

產(chǎn)出變量:①期望產(chǎn)出。本文選取水稻總產(chǎn)量為期望產(chǎn)出。②非期望產(chǎn)出。本文選取化肥施用的總環(huán)境成本為非期望產(chǎn)出。環(huán)境污染損失計(jì)算的方法有很多種,經(jīng)過幾十年的研究,形成了一系列從不同角度、不同經(jīng)濟(jì)理論評估環(huán)境損失的方法,如直接市場法、替代市場法、支付意愿法、能值分析法等。本文根據(jù)賴?yán)鸞9]的研究,基于能值分析理論,采取劑量影響法對我國水稻種植化肥施用的環(huán)境成本測算進(jìn)行估算。

對化肥施用的污染進(jìn)行分類,在此基礎(chǔ)上對各種可能存在的污染物的產(chǎn)生途徑、污染方式進(jìn)行歸納。本文將化肥對環(huán)境的影響分為三類,具體有大氣污染、土壤污染和水體污染三類。其中,大氣污染主要有NH3、N2O、NOx,土壤污染主要有硝酸鹽、硫酸鹽、重金屬元素鎘,水體污染主要有硝態(tài)氮、磷酸鹽。

歸納國內(nèi)外有關(guān)農(nóng)田氮磷素運(yùn)移轉(zhuǎn)化比例和污染物產(chǎn)生劑量的研究成果數(shù)據(jù)(表3),對我國化肥施用的營養(yǎng)物質(zhì)流向和污染產(chǎn)生物的影響劑量進(jìn)行分析,計(jì)算公式為:

Dosei=M×Cei×(We/Wf)

(3)

式中,Dosei為污染物i的產(chǎn)生劑量(t);M為氮肥或磷肥的施用折純量(t);Cei為營養(yǎng)元素的流轉(zhuǎn)系數(shù);We為所產(chǎn)生污染物的分子量;Wf為N或P2O5的分子量。

表3 氮肥、磷肥的養(yǎng)分運(yùn)移轉(zhuǎn)化比例(%)

對大氣、土壤和水體等環(huán)境質(zhì)量影響采用傷殘調(diào)整生命年法,估算化肥施用帶來污染物所造成的人類健康影響,計(jì)算方法為:

DALYi=Cdi×Dosei

(4)

式中,DALYi為污染物i造成的損害年累計(jì)數(shù)(a);Cdi為單位污染物引致的生命損害年數(shù)(a/kg);Dosei為污染物i的產(chǎn)生劑量。

依據(jù)國內(nèi)外有關(guān)單位勞動力的能值消費(fèi)數(shù)據(jù),乘以生命損害年累計(jì)數(shù),估算化肥施用對環(huán)境質(zhì)量影響的能值,最后匯總得出化肥環(huán)境影響的總能值。

(5)

式中,U為化肥環(huán)境影響的總能值成本，單位為sej;Emergyi為污染物的能值成本;Cm為單位勞動力的年能值消費(fèi)量;DALYi為某種污染物造成的損害年累計(jì)數(shù)。

依據(jù)我國各年份的能值貨幣比率數(shù)據(jù),折算化肥施用的綜合環(huán)境成本,計(jì)算公式為:

Emdollar=U/Cg

(6)

式中,Emdollar為化肥環(huán)境影響的宏觀經(jīng)濟(jì)價(jià)值，單位為￥;Cg為當(dāng)年的能值/貨幣比率,即一國家或地區(qū)單位時(shí)間內(nèi)使用的能值與GDP之比，單位為sej/￥;U為化肥環(huán)境影響的總能值成本，單位為sej。通過上述公式逐步計(jì)算,結(jié)果見表4。

表4 水稻種植化肥施用環(huán)境污染的宏觀經(jīng)濟(jì)價(jià)值

2004—2016年,我國水稻化肥施用的總環(huán)境成本,即宏觀經(jīng)濟(jì)價(jià)值,以年均5.9%的速度快速增加,表明由于化肥不合理施用造成的環(huán)境污染日益加劇,化肥不合理施用的問題仍存在,因此必須更加重視化肥不合理投入及其帶來的對環(huán)境質(zhì)量的影響。水稻綠色生產(chǎn)率測度投入—產(chǎn)出變量的描述性統(tǒng)計(jì)見表5。

表5 水稻種植化肥施用投入—產(chǎn)出變量的描述性統(tǒng)計(jì)

2 結(jié)果及分析

2.1 環(huán)境因素對全要素生產(chǎn)率測度影響

本文采用Malmquist-Luenberger指數(shù)研究考慮化肥環(huán)境成本下水稻全要素生產(chǎn)率的變動率。比較考慮化肥環(huán)境成本前后水稻全要素生產(chǎn)率值,發(fā)現(xiàn)大多數(shù)地區(qū)化肥環(huán)境成本下全要素生產(chǎn)率值(ML)低于不考慮化肥環(huán)境成本下全要素生產(chǎn)率(M)的值,說明我國水稻發(fā)展是以不合理施用化肥為基礎(chǔ)的,如果忽視環(huán)境作用,M指數(shù)會高估水稻的全要素生產(chǎn)率。

表6 2004—2016年考慮環(huán)境因素前后水稻全要素生產(chǎn)率指數(shù)結(jié)果

2.2 我國水稻綠色生產(chǎn)率指數(shù)結(jié)果分析

對三大主產(chǎn)區(qū)而言(表7),東北平原施肥強(qiáng)度最低,水稻綠色生產(chǎn)率均值明顯高于長江流域和東南沿海,說明東北平原能有效協(xié)調(diào)水稻生產(chǎn)和環(huán)境保護(hù)的關(guān)系,是三大主產(chǎn)區(qū)中水稻綠色生產(chǎn)率唯一達(dá)到有效的產(chǎn)區(qū),貫徹了水稻生產(chǎn)環(huán)境友好型發(fā)展。長江流域由于病蟲害影響較嚴(yán)重,施肥強(qiáng)度較高,因此水稻綠色生產(chǎn)率均值沒有達(dá)到有效。但從數(shù)據(jù)看,近三年來長江流域水稻綠色增長率一直處在生產(chǎn)前沿面上,表明有意識地控制了化肥的施用強(qiáng)度,提高了水稻綠色生產(chǎn)率。東南沿海地區(qū)經(jīng)濟(jì)發(fā)達(dá)、農(nóng)業(yè)生產(chǎn)占社會總生產(chǎn)比重小、土地資源稀缺,因此只能依靠增加化肥、農(nóng)藥等要素來增加水稻產(chǎn)量,未實(shí)現(xiàn)勞動力、土地、資源等要素的有效配置,水稻綠色生產(chǎn)率最小。因此，東南沿海地區(qū)應(yīng)深化技術(shù)改革,提升化肥的利用效率,選育抗病蟲害高效的品種,這才有利于促進(jìn)水稻發(fā)展由粗放加快向集約轉(zhuǎn)型。

表7 2004—2016年中國三大主產(chǎn)區(qū)水稻綠色生產(chǎn)率指數(shù)結(jié)果

(續(xù)表7)

地區(qū)200520062007200820092010201120122013201420152016均值湖南0.8731.0081.0771.1331.2060.8301.0111.0110.9611.1281.0621.0491.024江西0.9111.0121.0280.9980.5262.0681.0341.0341.0751.1251.0800.9531.027長江流域0.9120.9541.0341.0080.9511.0200.9860.9970.9521.0391.0311.0030.990浙江0.9521.0180.9611.0370.9680.9901.0000.9780.9141.0060.9951.0030.985福建0.7610.8661.0000.9660.8130.8880.8100.7880.7840.8781.0000.9590.872廣東0.8130.9800.9370.9201.0080.9441.0460.7930.7731.0690.9680.9340.927廣西0.7690.9801.0780.9730.9840.9260.9080.9910.9731.0280.9181.0060.958海南0.4011.1360.9580.9590.9270.9071.0300.9830.8961.0200.9791.0170.910東南沿海0.7110.9920.9860.9700.9370.9300.9550.9020.8650.9980.9710.9830.930吉林0.9651.0020.9171.0320.7721.0211.0220.7800.9750.9981.0811.0070.959遼寧0.8420.7201.6880.9870.9540.7141.1191.0031.1151.1581.1161.0531.013黑龍江0.9871.1801.0440.9791.2201.1391.1010.9900.9971.0171.0141.0711.059東北平原0.9290.9481.1730.9990.9650.9401.0800.9181.0271.0551.0691.0431.010

2.3 灌溉設(shè)施水平對水稻綠色生產(chǎn)率的影響

灌溉對水稻綠色生產(chǎn)率的影響機(jī)理:考察灌溉設(shè)施水平對水稻綠色生產(chǎn)率的影響,需要從期望產(chǎn)出和非期望產(chǎn)出兩個(gè)方面來考慮。一方面，因灌溉能明顯抑制農(nóng)田排水徑流量和滲漏量[10],同時(shí)改善土壤的通氣狀況,水稻在一定時(shí)期內(nèi)處于水分脅迫狀態(tài),限制了土壤有害物質(zhì)的產(chǎn)生,高效控制了氮磷流失[11],提升了化肥的利用率,因此灌溉能減少化肥的施用量,從而減少化肥投入造成的非期望產(chǎn)出,會促進(jìn)水稻綠色生產(chǎn)率的提升。另一方面,改善灌溉設(shè)施水平能保證水分的有效利用,從而保證期望產(chǎn)出的增加,會提升水稻的綠色生產(chǎn)率,機(jī)理見圖1。

圖1 灌溉對水稻綠色生產(chǎn)率的影響機(jī)理

實(shí)證分析:參考已有文獻(xiàn)和已有數(shù)據(jù)的可獲得性,分析節(jié)水灌溉對水稻綠色全要素生產(chǎn)率的影響,本文將選取農(nóng)業(yè)經(jīng)營規(guī)模、收入依賴程度、機(jī)械總動力、農(nóng)村經(jīng)濟(jì)發(fā)展水平作為控制變量。選取的變量為:①灌溉設(shè)施水平(ai)。用有效灌溉率表示,即有效灌溉面積占耕地面積的比重表示。②農(nóng)業(yè)經(jīng)營規(guī)模(as)。參照李燕的研究[12],用勞均耕地面積表示。③收入依賴程度(ae)。參照吉星星的研究[13],用家庭經(jīng)營收入占家庭純收入比重表示。④農(nóng)業(yè)機(jī)械化(am)。用農(nóng)業(yè)機(jī)械總動力表示。⑤農(nóng)村經(jīng)濟(jì)發(fā)展水平(ad)。參照潘丹的做法[2],用農(nóng)村居民純經(jīng)營性收入對數(shù)表示。本文按照我國三大水稻主產(chǎn)區(qū)長江流域、東南沿海、東北平原分別進(jìn)行檢驗(yàn),并選取以下面板OLS回歸模型。

TFPit=β0+β1aiit+β2asit+β3aeit+β4amit+β5adit+εit

(7)

式中,因變量TFPit為處理后的水稻綠色全要素生產(chǎn)率累積增長指數(shù),本文測算的水稻綠色全要素生產(chǎn)率指數(shù)是以上一年為1的環(huán)比指數(shù),參照李谷成的研究方法,將其轉(zhuǎn)化為2004年為1的累積增長指數(shù),然后進(jìn)行對數(shù)轉(zhuǎn)換作為被解釋變量,轉(zhuǎn)換公式為ln(1+TFP)[14]。模型中,i,t分別為省份和時(shí)間；β0為常數(shù)項(xiàng)；β1、β2、β3、β4、β5分別為各變量的待估計(jì)參數(shù)；εit為隨機(jī)擾動項(xiàng)。由Stata得出各變量的系數(shù)值。采用Stata進(jìn)行實(shí)證分析,模型回歸結(jié)果見表8。

表8 灌溉設(shè)施水平對水稻綠色生產(chǎn)率的影響

注:*、**、***分別表示在10%、5%、1%的水平上顯著,括號里的數(shù)據(jù)為p值。

由表8可知,基于全樣本數(shù)據(jù)分析,灌溉設(shè)施水平、農(nóng)業(yè)經(jīng)營規(guī)模、收入依賴程度、農(nóng)業(yè)機(jī)械化、農(nóng)村經(jīng)濟(jì)發(fā)展水平均通過了顯著性檢驗(yàn)。其中,農(nóng)業(yè)經(jīng)營規(guī)模、收入依賴程度對水稻綠色全要素生產(chǎn)率均具有正向影響,但農(nóng)業(yè)機(jī)械化對水稻綠色全要素生產(chǎn)率具有負(fù)面影響,原因在于我國幅員遼闊、地形復(fù)雜,導(dǎo)致水稻農(nóng)作物生產(chǎn)的綜合機(jī)械化水平偏低。因此，要因地制宜推進(jìn)農(nóng)業(yè)生產(chǎn)機(jī)械化。農(nóng)村經(jīng)濟(jì)發(fā)展水平對水稻綠色全要素生產(chǎn)率具有負(fù)面影響,原因是在經(jīng)濟(jì)發(fā)展的起初階段,為了增加水稻的收入與產(chǎn)出,人們往往依賴于自然資源的消耗和生產(chǎn)要素投入增加,忽視了對資源環(huán)境的保護(hù),從而導(dǎo)致水稻綠色全要素生產(chǎn)率較低,這與潘丹[2]的結(jié)論相符。

從總體上看,主要解釋變量灌溉設(shè)施水平對水稻綠色生產(chǎn)率具有正向影響,灌溉設(shè)施水平每增加1%,水稻綠色生產(chǎn)率將在5%的顯著性水平下增加22.6%,表明節(jié)水灌溉對于水稻綠色增長率的提高具有顯著的促進(jìn)作用。改善灌溉設(shè)施水平可提高灌溉用水效率,減少氮磷流失,緩解化肥的面源污染問題,提高水稻綠色生產(chǎn)率?？傊?提高灌溉設(shè)施水平可緩解由于不合理施用化肥所導(dǎo)致的農(nóng)業(yè)面源污染問題,促進(jìn)水稻綠色生產(chǎn)率的提高。

從各主產(chǎn)區(qū)的模型結(jié)果看,長江流域和東北平原顯著為正,表明長江流域和東北平原灌溉設(shè)施水平每增加1%,水稻綠色生產(chǎn)率分別在1%和10%顯著性水平上分別增加38.3%和92.3%。然而,東南沿海顯著為負(fù),灌溉設(shè)施水平每增加1%,水稻綠色生產(chǎn)率將降低33.4%。原因是：首先，東南沿海地區(qū)地形復(fù)雜,多為丘陵低山等,導(dǎo)致用水灌溉效率低；其次,災(zāi)害繁多對節(jié)水灌溉設(shè)施損害嚴(yán)重,導(dǎo)致節(jié)水灌溉設(shè)施普遍較薄弱,用水灌溉效率低。因此,我國要加大東南沿海的節(jié)水灌溉設(shè)施力度,積極采用灌溉新技術(shù),提升用水灌溉效率,積極抵御自然災(zāi)害對農(nóng)作物的影響。

3 結(jié)論

本文基于2004—2016年我國18個(gè)省份的面板數(shù)據(jù),利用Malmquist-Luenberger指數(shù)測算了我國化肥環(huán)境成本下的水稻綠色生產(chǎn)率。在此基礎(chǔ)上,利用面板OLS回歸模型分析了灌溉設(shè)施水平對于水稻綠色生產(chǎn)率的影響。得出以下結(jié)論:①化肥施用結(jié)構(gòu)雖然發(fā)生了變化,但化肥的總施用量在近幾年卻以年均1.7%的速度呈現(xiàn)上升的態(tài)勢。在我國三大主產(chǎn)區(qū)中,東南沿海地區(qū)年均化肥施用強(qiáng)度最高,其次是長江流域,最低的是東北平原。②東北平原是三大主產(chǎn)區(qū)中水稻綠色生產(chǎn)率唯一達(dá)到有效的產(chǎn)區(qū),長江流域和東北平原均沒有實(shí)現(xiàn)水稻生產(chǎn)和環(huán)境的協(xié)調(diào)發(fā)展。③總體看,灌溉設(shè)施水平有助于提升水稻綠色生產(chǎn)率；從各個(gè)主產(chǎn)區(qū)來看,長江流域和東北平原的灌溉設(shè)施水平有助于提升水稻綠色生產(chǎn)率,東南沿海的灌溉設(shè)施水平對水稻綠色生產(chǎn)率的提升起到了阻礙作用。

4 建議

基于本文對2004—2016年我國18個(gè)三大水稻種植區(qū)域18個(gè)省份水稻綠色生產(chǎn)率的研究結(jié)論，本文從政府視角出發(fā)，提出以下三個(gè)方面的建議。

首先，灌溉種類繁多，政府應(yīng)因地制宜地選擇和推廣農(nóng)業(yè)節(jié)水灌溉技術(shù)。水稻灌溉種類主要分為滴灌、霧噴灌、膜下灌與滲灌4大類，應(yīng)根據(jù)不同的地區(qū)的種植特點(diǎn)采用不同的灌溉方式，而不能都采取最普通的灌溉方式。在本文研究的三大主產(chǎn)區(qū)之中，東南沿海地區(qū)的灌溉設(shè)施水平最為薄弱，尤其需要注重節(jié)水灌溉技術(shù)在東南沿海地區(qū)的發(fā)展與應(yīng)用。當(dāng)?shù)卣畱?yīng)加強(qiáng)完善水利設(shè)施，重視灌溉對于水稻發(fā)展的作用。

其次，加大農(nóng)業(yè)科研投入，增加對灌溉、施肥技術(shù)、水稻品種等的研發(fā)，提升灌溉用水效率和化肥利用率，注重水肥合一，提升水稻的綠色生產(chǎn)率。隨著人們生活水平的提高，吃飯已經(jīng)不再是人們的最主要訴求，健康、綠色、安全逐漸成為人們選擇稻米的標(biāo)準(zhǔn)，因此在提升水稻綠色生產(chǎn)率方面投入更多的科研顯得尤為重要。

第三，加強(qiáng)對資源環(huán)境友好的生產(chǎn)行為補(bǔ)貼。政府可以對秸稈還田、節(jié)水灌溉、施用農(nóng)家肥等環(huán)境友好型生產(chǎn)行為給予一定的補(bǔ)貼，引導(dǎo)生產(chǎn)者改變生產(chǎn)方式，采納有利于農(nóng)業(yè)資源環(huán)境的生產(chǎn)行為。此外，為響應(yīng)國家號召，地方政府也應(yīng)推出農(nóng)作物垃圾分類回收的政策。例如，農(nóng)藥包裝廢棄物的有償回收，不僅可以保護(hù)耕地環(huán)境與水資源，還可以節(jié)省農(nóng)戶的種植成本。