毛中明，楊豪毅，劉陳思雨，馬 珠

基于LMDI方法的恩施州化肥施用驅(qū)動(dòng)因素分析

毛中明，楊豪毅，劉陳思雨，馬 珠

（中南民族大學(xué)經(jīng)濟(jì)學(xué)院，武漢 430074）

近年來(lái)由化肥過(guò)量施用造成的農(nóng)業(yè)面源污染問(wèn)題日益突出，加大農(nóng)業(yè)污染治理力度，探究化肥施用強(qiáng)度的變化特征及其驅(qū)動(dòng)因素，可以為減肥增效、推進(jìn)農(nóng)業(yè)綠色發(fā)展提供科學(xué)依據(jù)。研究以湖北省恩施州為實(shí)證對(duì)象，在研究其時(shí)空變化特征的基礎(chǔ)上，通過(guò)綜合考慮化肥施用效率、區(qū)域施肥結(jié)構(gòu)、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積等因素，構(gòu)建對(duì)數(shù)平均迪氏指數(shù)法（Logarithmic Mean Divisia Index，LMDI）指數(shù)分解模型與面板回歸模型。結(jié)果表明：1）近10 a來(lái)，恩施州化肥施用量呈現(xiàn)先增長(zhǎng)后平穩(wěn)再下降趨勢(shì)，氮、磷、鉀比例由2007年的1∶0.30∶0.15增長(zhǎng)到2018年的1∶0.39∶0.26，逐漸向世界發(fā)達(dá)國(guó)家水平靠攏。截至2018年，恩施州化肥施用共有1個(gè)低強(qiáng)度區(qū)、2個(gè)中強(qiáng)度區(qū)、4個(gè)高強(qiáng)度區(qū)、1個(gè)超高強(qiáng)度區(qū)。2）2016年恩施州化肥施用量比上年減少7 008.32 t，2017年減少3 341.65 t，2018年減少15 059.81 t。其中化肥施用效率提高使化肥施用量在近3 a分別減少11 159.44、14 051.95和8 644.21 t，復(fù)種指數(shù)使化肥施用量在2017年減少5 229.76 t，2018年減少30 889.98 t，2016年效應(yīng)數(shù)值不大。農(nóng)業(yè)技術(shù)發(fā)展使化肥施用量增加了4 089.71、11 136.97和12 215.74 t，區(qū)域結(jié)構(gòu)的驅(qū)動(dòng)效應(yīng)不大。3）經(jīng)濟(jì)發(fā)展、農(nóng)業(yè)結(jié)構(gòu)調(diào)整、金融規(guī)模擴(kuò)大對(duì)化肥施用量起正向作用，人口數(shù)量增加對(duì)化肥施用量起負(fù)向作用。建議轉(zhuǎn)變農(nóng)業(yè)技術(shù)發(fā)展方向由增產(chǎn)至高效增產(chǎn)，合理調(diào)整有效耕地面積，進(jìn)一步提升化肥施用效率，培育生態(tài)安全綠色發(fā)展模式，促進(jìn)化肥減量增效。

農(nóng)業(yè)；化肥；時(shí)空特征；LMDI模型；驅(qū)動(dòng)因素；恩施州

0 引 言

在增產(chǎn)導(dǎo)向型農(nóng)業(yè)政策推動(dòng)下，中國(guó)農(nóng)業(yè)糧食產(chǎn)量從1980年32 055.53萬(wàn)t增長(zhǎng)到2019年66 384.21萬(wàn)t，實(shí)現(xiàn)了快速增長(zhǎng)[1]。與此同時(shí)，存在另外一個(gè)易被忽視又極為重要的事實(shí)：中國(guó)化肥施用量從1980年1 269.40萬(wàn)t增長(zhǎng)到2019年5 403.59萬(wàn)t，化肥施用強(qiáng)度由1980年的82 kg/hm2增長(zhǎng)到2019年的293 kg/hm2[2]，大量化肥施用帶來(lái)水體污染[3]、土壤污染等一系列農(nóng)業(yè)面源污染[4]問(wèn)題。因此，客觀分析農(nóng)業(yè)化肥施用強(qiáng)度影響因素及其作用機(jī)制，找準(zhǔn)降低化肥施用強(qiáng)度有效途徑，為“十四五”期間涉及農(nóng)業(yè)綠色發(fā)展有關(guān)決策提出針對(duì)性政策建議[5]，成為當(dāng)前農(nóng)業(yè)發(fā)展的研究重點(diǎn)之一。

為了控制化肥過(guò)度施用所帶來(lái)的面源污染，國(guó)內(nèi)外在化肥安全用量及其影響因素方面做了大量研究。Green等[6]通過(guò)Logit模型確定耕作制度、作物品種、信貸獲取、非農(nóng)勞動(dòng)力收入是影響肥料采用的主要因素，因此影響這些變量的政策對(duì)肥料用量的影響最大。Endale[7]通過(guò)消費(fèi)模型得出，戶主受教育程度是影響肥料施用最重要的變量。畜牧業(yè)所有權(quán)，擁有土地的多少，信貸數(shù)量以及接受過(guò)義務(wù)教育的家庭成員數(shù)量是對(duì)肥料施用產(chǎn)生積極影響的其他因素。同時(shí)確定了解決肥料使用及其消耗問(wèn)題的優(yōu)先干預(yù)措施領(lǐng)域，供應(yīng)方干預(yù)的最優(yōu)先領(lǐng)域是化肥價(jià)格。幾乎50％的農(nóng)民認(rèn)為價(jià)格是他們施用化肥最大的限制。這就需要考慮其他替代手段，例如作物特定的肥料部分補(bǔ)貼和現(xiàn)金轉(zhuǎn)移。非洲的肥料施用率僅為亞洲、歐洲和美洲的肥料使用率的10%～20%，具有較強(qiáng)的借鑒意義，Burke等[8]通過(guò)對(duì)“非洲綠色革命”的研究與實(shí)地考察，得出發(fā)展農(nóng)作物多樣化，提高化肥施用效率，可以有效減少化肥施用。國(guó)內(nèi)學(xué)者主要從種植結(jié)構(gòu)[9]、播種面積[10-11]、施肥強(qiáng)度[12]等方面進(jìn)行分析。欒江等[13]利用指標(biāo)分解方法分析種植結(jié)構(gòu)調(diào)整、施肥強(qiáng)度增加、播種面積變化對(duì)化肥施用量變動(dòng)的影響，但僅使用3個(gè)指標(biāo)衡量化肥施用量變動(dòng)還不夠全面。趙明正等[14]從全國(guó)整體、分區(qū)域、分品種3個(gè)維度出發(fā)分析種植結(jié)構(gòu)、區(qū)域結(jié)構(gòu)、播種面積對(duì)于化肥施用的影響，指出種植結(jié)構(gòu)變化是導(dǎo)致中國(guó)化肥施用量下降的主要驅(qū)動(dòng)因素。上述研究中所選取指標(biāo)較為寬泛，且多是基于田間試驗(yàn)、農(nóng)戶調(diào)查等方法，而綜合考慮各因素與農(nóng)業(yè)化肥施用間聯(lián)系方向與作用程度大小的研究甚少。地處鄂西的恩施土家族苗族自治州是一個(gè)典型的農(nóng)耕區(qū)，農(nóng)業(yè)在全州國(guó)民經(jīng)濟(jì)中的地位極其重要，其化肥施用在湖北一直處于高位狀態(tài)。因此本文通過(guò)構(gòu)建化肥施用驅(qū)動(dòng)因素深層次分解的方法，對(duì)恩施州2007—2018年化肥施用的時(shí)空變化特征及驅(qū)動(dòng)因素、深層次影響因素進(jìn)行研究。一方面，完善化肥施用影響評(píng)價(jià)方法，另一方面為恩施州實(shí)行化肥施用總量控制和農(nóng)業(yè)綠色發(fā)展提供政策依據(jù)。

1 數(shù)據(jù)與方法

1.1 研究區(qū)概況

恩施州位于湖北省西南部（109°29′10″E，30°15′12″N），地處湘、鄂、渝三?。ㄊ校┙粎R處，地形錯(cuò)綜復(fù)雜，地勢(shì)高低懸殊，呈現(xiàn)出極其明顯的氣候垂直地域差異。2018年恩施州土地總面積2 406 026 hm2，其中農(nóng)業(yè)用地面積367 233 hm2。全州巖溶發(fā)育強(qiáng)烈，暗河伏流多，地表水資源總量為124.92億m3，地下水資源總量44.13億m3，為當(dāng)?shù)丶Z食安全與農(nóng)業(yè)灌溉提供良好的保障。恩施州大氣環(huán)境良好，生態(tài)價(jià)值突出，2018年恩施州空氣質(zhì)量?jī)?yōu)良天數(shù)為352 d，縣市城區(qū)可吸入顆粒物（PM10）濃度均值為34g/m3，細(xì)顆粒物（PM2.5）濃度均值為23g/m3。近年來(lái)，恩施州大量建設(shè)生態(tài)農(nóng)業(yè)產(chǎn)業(yè)園、無(wú)公害綠色農(nóng)產(chǎn)品生產(chǎn)基地，建立并完善生態(tài)安全格局，為該地及周邊地區(qū)提供良好的自然環(huán)境。

恩施州包含恩施市、利川市、建始縣、巴東縣、宣恩縣、咸豐縣、來(lái)鳳縣、鶴峰縣8個(gè)縣市。2005以來(lái)恩施州經(jīng)濟(jì)發(fā)展趨勢(shì)具有階段性的特征，2005—2011年恩施州經(jīng)濟(jì)高速發(fā)展，2012—2018年恩施州經(jīng)濟(jì)發(fā)展速度逐漸放緩。GDP增長(zhǎng)速度與湖北省GDP增長(zhǎng)速度趨勢(shì)大體相同，但2009年、2012—2015年恩施州GDP增長(zhǎng)速度高于湖北省GDP增長(zhǎng)速度。2005以來(lái)恩施州產(chǎn)業(yè)結(jié)構(gòu)調(diào)整明顯，第一產(chǎn)業(yè)整體呈現(xiàn)下降趨勢(shì)，第二、三產(chǎn)業(yè)大體呈現(xiàn)穩(wěn)步上升趨勢(shì)，在2011年第二產(chǎn)業(yè)占比超過(guò)第一產(chǎn)業(yè)占比。

近年來(lái)（2005—2018年），農(nóng)林牧漁業(yè)占恩施州各縣市生產(chǎn)總值比例的平均水平呈現(xiàn)下降趨勢(shì)，均值由2005年40.88%下降到2018年的17.65%（圖1）。利川市、宣恩縣、建始縣、咸豐縣農(nóng)林牧漁業(yè)占比超過(guò)平均水平，其余四縣市低于平均水平。

1.2 LMDI指數(shù)分解模型

采用對(duì)數(shù)平均迪氏指數(shù)法（Logarithmic Mean Divisia Index，LMDI）[15]方法探討化肥施用的驅(qū)動(dòng)因素，該方法具有完全分解、無(wú)殘差、無(wú)零值等優(yōu)點(diǎn)。根據(jù)LMDI方法構(gòu)建恩施州模型（式（1））和分縣域模型（式（2））：

式中C表示年恩施州化肥施用總量，kg；表示所屬恩施州的8個(gè)縣域（恩施市、利川市、建始縣、巴東縣、宣恩縣、咸豐縣、來(lái)鳳縣、鶴峰縣）；C表示縣化肥施用量，kg；C表示年縣的化肥施用量，kg；Q表示年縣的農(nóng)作物產(chǎn)量，t；S表示年縣的農(nóng)作物種植面積，hm2；T表示年縣的耕地面積，hm2；Q表示年恩施州的農(nóng)作物產(chǎn)量，t；S表示年恩施州的農(nóng)作物種植面積，hm2；T表示年恩施州的耕地面積，hm2。故表示化肥施用量與農(nóng)作物產(chǎn)量的比值（kg/t），由LMDI指數(shù)分解模型性質(zhì)，可表征為化肥施用效率；表示各個(gè)縣農(nóng)作物產(chǎn)量與恩施州農(nóng)作物產(chǎn)量的比值，即區(qū)域結(jié)構(gòu)；表示農(nóng)作物產(chǎn)量與農(nóng)作物種植面積的比值（t/hm2），即農(nóng)業(yè)技術(shù)發(fā)展；表示農(nóng)作物種植面積與耕地面積的比值，即復(fù)種指數(shù)；是耕地面積，hm2。

由式（1）可知，恩施州化肥施用量受恩施州化肥施用效率、區(qū)域結(jié)構(gòu)、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積5個(gè)因素的影響。一段時(shí)期后恩施州化肥施用量的差異可以分解為

式中ΔC、ΔC、ΔC、ΔC、ΔC分別為化肥施用效率、區(qū)域結(jié)構(gòu)、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積對(duì)化肥施用量的影響。

對(duì)于某個(gè)縣，該縣的化肥施用量?jī)H受該縣化肥施用效率、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積4個(gè)因素影響。一段時(shí)期后縣化肥施用量的差異可以分解為

根據(jù)LMDI分解技術(shù)的基本原理，測(cè)度逐年各驅(qū)動(dòng)因素效應(yīng)，有

其中：

1.3 影響因素面板回歸模型

由式（4）可知，恩施州各縣的化肥施用量受化肥施用效率、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積4個(gè)因素所驅(qū)動(dòng)。進(jìn)一步分析，這些驅(qū)動(dòng)因素背后有更深層次的影響因素，如圖2所示。

1）人口數(shù)量，人口的增多會(huì)導(dǎo)致城鎮(zhèn)化速度加快，進(jìn)而會(huì)造成人多地少的矛盾，增大復(fù)種指數(shù)是解決這種矛盾的最有效的方法，所以人口數(shù)量在一定程度上影響了復(fù)種指數(shù)[16-17]。2）經(jīng)濟(jì)發(fā)展，一方面影響著農(nóng)業(yè)勞動(dòng)者的機(jī)會(huì)成本，導(dǎo)致農(nóng)業(yè)勞動(dòng)者轉(zhuǎn)向其他行業(yè)，使部分耕地成為荒地，另一方面，經(jīng)濟(jì)發(fā)展的同時(shí)人們會(huì)更加注重生活質(zhì)量，環(huán)境忍受度也會(huì)隨之改變，從而減少化肥施用量[18]。3）農(nóng)業(yè)結(jié)構(gòu)，湖北省多項(xiàng)促進(jìn)農(nóng)民增收的政策，帶動(dòng)了農(nóng)林牧漁業(yè)、發(fā)展，促進(jìn)了農(nóng)業(yè)結(jié)構(gòu)調(diào)整。農(nóng)業(yè)結(jié)構(gòu)的變化會(huì)帶動(dòng)耕地面積的變化，林牧漁業(yè)的發(fā)展使部分耕地轉(zhuǎn)為園林、養(yǎng)殖地等，實(shí)現(xiàn)多樣化經(jīng)營(yíng)，從而耕地面積減少[19]。4）農(nóng)民收入水平，一方面，農(nóng)民收入水平的提高，使得農(nóng)民有更多的資金購(gòu)買化肥，增加了化肥施用量，另一方面，由于農(nóng)民收入水平的提高，更多的農(nóng)民開(kāi)始購(gòu)買新品種化肥，提高化肥施用效率，又減少了化肥施用量[20]。5）金融規(guī)模，金融規(guī)模的擴(kuò)大促進(jìn)內(nèi)源融資、外源融資等多種融資方式的產(chǎn)生，這些融資方式可以降低融資約束和融資成本，更好地滿足了企業(yè)對(duì)資金的需求，大量的資金注入對(duì)農(nóng)業(yè)技術(shù)創(chuàng)新以及產(chǎn)品創(chuàng)新提供基礎(chǔ)，進(jìn)而減少化肥施用量[21-22]。根據(jù)上述分析，設(shè)定面板回歸模型如下：

式中是化肥施用量，t；是人口數(shù)量（萬(wàn)人），用來(lái)考察人口變化的影響；是人均GDP（萬(wàn)元），用來(lái)考察經(jīng)濟(jì)發(fā)展的影響；是農(nóng)業(yè)產(chǎn)值（萬(wàn)元）占農(nóng)林牧漁總產(chǎn)值（萬(wàn)元）之比，用來(lái)考察農(nóng)業(yè)結(jié)構(gòu)的影響；是農(nóng)村居民人均可支配收入（元），用來(lái)考察農(nóng)民收入水平的影響；是金融機(jī)構(gòu)年末存貸款總額（萬(wàn)元）與當(dāng)?shù)谿DP（萬(wàn)元）的比值，用來(lái)考察金融規(guī)模的影響；1、1～5為系數(shù)，ε誤差。

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

本文所采用的數(shù)據(jù)中，農(nóng)作物產(chǎn)量來(lái)源于2015－2018年的《全國(guó)農(nóng)產(chǎn)品成本收益資料匯編》；種植面積、耕地面積來(lái)源于2015－2018年的《中國(guó)農(nóng)村統(tǒng)計(jì)年鑒》；化肥施用量、人口數(shù)量、GDP、農(nóng)業(yè)產(chǎn)值、金融機(jī)構(gòu)年末貸款額來(lái)源于2009－2018年的《恩施州統(tǒng)計(jì)年鑒》。由于已有資料并沒(méi)有給出復(fù)合肥中氮、磷、鉀的比例，根據(jù)市場(chǎng)調(diào)查情況，并參考有關(guān)文獻(xiàn)[23-24]和聯(lián)合國(guó)糧農(nóng)組織對(duì)中國(guó)復(fù)合肥中養(yǎng)分比例的計(jì)算方法，本文化肥施用量數(shù)據(jù)中的復(fù)合肥按氮磷鉀含量比例1∶1∶1處理。

2 結(jié)果與分析

2.1 恩施州化肥施用時(shí)空分異

2.1.1 化肥施用時(shí)間變化特征

2007—2018年恩施州化肥施用總量呈現(xiàn)先增長(zhǎng)再下降趨勢(shì)（圖3），由2007年的184 653 t增加到2012年的291 554 t達(dá)到最高點(diǎn)，其中氮肥和復(fù)合肥上升趨勢(shì)明顯，磷肥、鉀肥稍有上升；經(jīng)過(guò)3a平穩(wěn)調(diào)整，在2015年開(kāi)始實(shí)現(xiàn)連續(xù)3 a負(fù)增長(zhǎng)，2018年施肥總量為258 534 t。說(shuō)明2015年農(nóng)業(yè)部印發(fā)的《到2020年化肥使用量零增長(zhǎng)行動(dòng)方案》以及提出的相關(guān)政策措施取得顯著減肥成效。

另外，氮（N）、磷（P2O5）、鉀（K2O）化肥比例由2007年的1∶0.30∶0.15增長(zhǎng)到2018年的1∶0.39∶0.26，高于中國(guó)平均氮磷鉀化肥比例1∶0.37∶0.24[25]并逐漸向世界發(fā)達(dá)國(guó)家水平1:0.50:0.50靠攏[26]。

2.1.2 化肥施用空間變化特征

根據(jù)恩施州各縣化肥施用量與播種面積的比值，分別計(jì)算各縣2007年以來(lái)的平均化肥施用強(qiáng)度。以國(guó)家生態(tài)鄉(xiāng)鎮(zhèn)建設(shè)要求[27]化肥施用強(qiáng)度不超過(guò)250 kg/hm2的標(biāo)準(zhǔn)，以及張福鎖等[28]對(duì)化肥施用強(qiáng)度等級(jí)以100 kg間隔劃分法，按照化肥施用強(qiáng)度<250、250～350、351～450、>450 kg/hm2的變化范圍，把恩施州各縣化肥投入分為低施肥強(qiáng)度區(qū)、中施肥強(qiáng)度區(qū)、高施肥強(qiáng)度區(qū)、超高施肥強(qiáng)度區(qū)4個(gè)類型，2007－2018年化肥施用強(qiáng)度類型的空間分布如圖4所示。

1）由圖4可知，化肥施用強(qiáng)度區(qū)域差異越來(lái)越大，由2007年的1個(gè)低強(qiáng)度、6個(gè)中強(qiáng)度、1個(gè)高強(qiáng)度3種類型增長(zhǎng)至2018年1個(gè)低強(qiáng)度、2個(gè)中強(qiáng)度、4個(gè)高強(qiáng)度、1個(gè)超高強(qiáng)度4種類型。2007年僅有恩施市為低施肥強(qiáng)度區(qū)，巴東縣為高施肥強(qiáng)度區(qū)，其余各縣均為中施肥強(qiáng)度區(qū)；2012年恩施市上升為中強(qiáng)度施肥區(qū)，利川市和咸豐縣從中強(qiáng)度上升至高強(qiáng)度區(qū)，而鶴峰縣的化肥施用強(qiáng)度超過(guò)了450 kg/hm2，成為恩施州首個(gè)超高強(qiáng)度施肥區(qū)；2015年恩施州的化肥施用達(dá)到頂峰，超高強(qiáng)度施肥區(qū)域有咸豐縣、建始縣、鶴峰縣3個(gè)，但是來(lái)鳳縣化肥施用有所下降成為了唯一的低強(qiáng)度施肥區(qū)，其余地區(qū)均為中強(qiáng)度；2018年由于零增長(zhǎng)政策的實(shí)施顯現(xiàn)成效[29]，各地施肥強(qiáng)度較2015年均有不同程度的下降，但是鶴峰縣施肥強(qiáng)度依然高達(dá)555.20 kg/hm2，是施肥強(qiáng)度最低縣來(lái)鳳縣207.44 kg/hm2的2.676倍，為全州唯一超高強(qiáng)度施肥區(qū)。咸豐縣、建始縣降由超高強(qiáng)度至高強(qiáng)度，其余地區(qū)沒(méi)有變化。

2）高強(qiáng)度化肥施用區(qū)與低強(qiáng)度化肥施用區(qū)的差距逐漸趨于平穩(wěn)。2007年，巴東縣和恩施市分別是全州唯一高強(qiáng)度施肥區(qū)和低強(qiáng)度施肥區(qū)；從2012年開(kāi)始，恩施州最高最低強(qiáng)度施肥區(qū)域分別是鶴峰縣和來(lái)鳳縣。2007、2012、2015、2018這4 a的化肥施用最高區(qū)域與最低區(qū)域的強(qiáng)度比值分別為1.593、2.563、2.565、2.676，呈先上升后平穩(wěn)的趨勢(shì)。說(shuō)明所屬恩施州的各縣在農(nóng)業(yè)發(fā)展方面的步伐較為一致。

2.2 恩施州化肥施用驅(qū)動(dòng)因素分析

2.2.1 LMDI全州模型分解

由表1可知，2016年恩施州化肥施用量比上年減少7 008.32 t，2017年恩施州化肥施用量比上年減少3 341.65 t，2018年恩施州化肥施用量比上年減少15 059.81 t。其中化肥施用效率使化肥施用量在近3a分別減少11 159.44、14 051.95和8 644.21 t，說(shuō)明恩施州化肥施用效率在逐年提高。復(fù)種指數(shù)使化肥施用量在2017年減少3 745.02 t，2018年減少32 321.67 t，說(shuō)明農(nóng)用耕地得到了更合理的利用，2016年效應(yīng)數(shù)值不大。農(nóng)業(yè)技術(shù)發(fā)展使化肥施用量增加了3 396.37、9 587.52和13 767.38 t，說(shuō)明恩施州近些年大力發(fā)展農(nóng)業(yè)技術(shù)導(dǎo)致化肥施用量水漲船高。區(qū)域結(jié)構(gòu)的驅(qū)動(dòng)效應(yīng)不大。

表1 全州化肥施用量驅(qū)動(dòng)因素分解結(jié)果

2.2.2 LMDI分縣模型分解

由表2可知，2015—2016年除鶴峰縣外其他縣的化肥施用量均為下降，其中下降超1 000 t的有恩施市、利川市、建始縣和咸豐縣。2016—2017年，巴東縣和咸豐縣的化肥施用量為上升，且增幅都超過(guò)1 000 t，其余縣化肥施用量均為下降，下降超過(guò)1 000 t的有恩施市和鶴峰縣。2017—2018年宣恩縣化肥施用量為上升，僅40 t，其余各縣均為下降，下降超過(guò)1 000 t的有恩施市、建始縣、巴東縣和鶴峰縣，下降最多的是建始縣為8 283 t。

表2 分縣（市）化肥施用量驅(qū)動(dòng)因素分解結(jié)果

由表2和圖5可知，化肥施用效率變化使2015－2016年全州化肥施用量下降了11 159.44 t，2016－2017年全州化肥施用量下降了14 051.95 t，2017－2018年全州化肥施用量下降了8 644.21 t。如圖5所示，該效應(yīng)在2015－2016年使化肥施用量下降超過(guò)2 000 t的是恩施市和咸豐縣；2016－2017年使化肥施用量下降超過(guò)2 000 t的是恩施市、利川市和鶴峰縣；2017－2018年使化肥施用量下降超過(guò)2 000 t的是建始縣。恩施市作為全州經(jīng)濟(jì)發(fā)展最好的地區(qū)，其經(jīng)濟(jì)水平最高，獲得務(wù)農(nóng)相關(guān)技術(shù)培訓(xùn)最多，這些都是提高化肥施用效率的重要因素[30]。其余地區(qū)則是政府實(shí)施增加農(nóng)民收入、鼓勵(lì)土地流轉(zhuǎn)政策等環(huán)境使然，提高了化肥施用效率[31]。該效應(yīng)在2016－2017年沒(méi)有使哪個(gè)地區(qū)化肥施用量大量增長(zhǎng)，2017－2018年使利川市化肥施用量增加超過(guò)2 000 t，這是由于利川市在2017年減少了復(fù)合肥的施用，氮磷鉀比例低于之前年份導(dǎo)致化肥結(jié)構(gòu)不合理從而降低了化肥施用效率。

表2表明，農(nóng)業(yè)技術(shù)發(fā)展在2015－2016年使全州化肥施用量增加了4 089.71 t，2016－2017年使全州化肥施用量增加了11 136.97 t，2017－2018年使全州化肥施用量增加了12 215.74 t。如圖6所示，該效應(yīng)在2015－2016年使化肥施用量下降超過(guò)2 000 t的是建始縣；2016－2018年沒(méi)有使化肥施用量下降超過(guò)2 000 t的地區(qū)。由于建始縣起步晚進(jìn)展慢，經(jīng)過(guò)努力在近年才形成有一定市場(chǎng)、規(guī)模的生產(chǎn)基地、專業(yè)合作社，同時(shí)政府按照“一產(chǎn)一業(yè)”、“一村一品”、“一人一專”的要求，做好了對(duì)農(nóng)技推廣人員的培訓(xùn)工作，結(jié)合基層農(nóng)技推廣體系改革與建設(shè)項(xiàng)目實(shí)施方案，強(qiáng)化技術(shù)指導(dǎo)、輻射帶動(dòng)農(nóng)戶增收，實(shí)現(xiàn)減肥增效[32-33]。該效應(yīng)在2015－2016年使化肥施用量增加超過(guò)2 000 t的是恩施市和鶴峰縣，2016－2017年使化肥施用量增加超過(guò)2 000 t的是利川市、咸豐縣和鶴峰縣，2017－2018年使化肥施用量增加超過(guò)2 000 t的是恩施市、利川市、宣恩縣和建始縣。由于上述各縣的技術(shù)提升推動(dòng)糧食需求，導(dǎo)致化肥施用量也不斷增加[14]。

2015－2016年，復(fù)種指數(shù)的降低使得全州化肥施用量減少了384.78 t，2016－2017年使全州化肥施用量減少了5 229.76 t，2017－2018年使全州化肥施用量減少了30 889.98 t。如圖7所示，該效應(yīng)在2015－2016年使化肥施用量下降超過(guò)2 000 t的是鶴峰縣；2016－2017年使化肥施用量下降超過(guò)2 000 t的是建始縣；2017－2018年使化肥施用量下降超過(guò)2 000 t的是恩施市、利川市、建始縣和宣恩縣。其中建始縣和鶴峰縣多處山地，機(jī)耕率和耕地的有效灌溉率很難提升，導(dǎo)致播種面積減少?gòu)亩档蛷?fù)種指數(shù)，恩施市、利川市和宣恩縣人口相對(duì)密集，糧食需求較大，需提升耕地集約度從而改變復(fù)種指數(shù)[35]。該效應(yīng)在2015－2016年使化肥施用量增加超過(guò)2 000 t的是咸豐縣，2016－2017年使化肥施用量增加超過(guò)2 000 t的是咸豐縣，2017－2018年沒(méi)有使化肥施用量增加超過(guò)2 000 t的地區(qū)。咸豐縣把玉米高產(chǎn)高效模式作為高產(chǎn)創(chuàng)建的主要技術(shù)進(jìn)行廣泛推廣，擴(kuò)大種植面積，使復(fù)種指數(shù)增加[36]。

耕地面積在2015－2016年使全州化肥施用量增加了445.51 t，2016－2017年使全州化肥施用量增加了4 799.75 t，2017－2018年使全州化肥施用量增加了12 225.46 t。如圖8所示，該效應(yīng)在2015－2016年使化肥施用量下降超過(guò)2 000 t的是恩施市、利川市、巴東縣、宣恩縣、來(lái)鳳縣和鶴峰縣；2016－2018年沒(méi)有使任何地區(qū)的化肥施用量下降超過(guò)2 000 t。該效應(yīng)在2016－2017年使化肥施用量增加超過(guò)2 000 t的是建始縣，2017－2018年使化肥施用量增加超過(guò)2 000 t的是恩施市、建始縣。說(shuō)明上述地區(qū)的耕地面積較之前有較大增長(zhǎng)。

2.3 恩施州化肥施用量影響因素面板數(shù)據(jù)分析

2.3.1 變量測(cè)度

本文使用2009－2018年間恩施州所屬各縣市的面板數(shù)據(jù)作為研究樣本，對(duì)恩施州化肥施用影響因素進(jìn)行分析。表3為各變量的統(tǒng)計(jì)性描述。

2.3.2模型估計(jì)結(jié)果

采用固定效應(yīng)模型對(duì)面板數(shù)據(jù)進(jìn)行回歸。恩施州自2008年以來(lái)實(shí)行“東部雜糧西部米”的發(fā)展目標(biāo)和思路，由于不同農(nóng)作物所需要化肥施用量不同，故本文把東部（巴東、建始、恩施、鶴峰）和西部（利川、咸豐、宣恩、來(lái)鳳）分別進(jìn)行回歸?；貧w結(jié)果如表4。

對(duì)恩施州整體回歸，經(jīng)濟(jì)發(fā)展在5%的水平上顯著為正，說(shuō)明隨著經(jīng)濟(jì)的發(fā)展，恩施州化肥施用量在不斷上升。農(nóng)業(yè)結(jié)構(gòu)在10%的水平顯著且系數(shù)為正，說(shuō)明農(nóng)業(yè)產(chǎn)值提高會(huì)促進(jìn)化肥的施用量。人口數(shù)量在1%水平上顯著且系數(shù)為負(fù)，說(shuō)明隨著人口的增長(zhǎng)和城市化的推進(jìn)，耕地面積和復(fù)種指數(shù)發(fā)生相應(yīng)的變化從而導(dǎo)致化肥施用的減少。金融規(guī)模在10%水平上顯著為正，說(shuō)明金融資金注入較完善的地區(qū)化肥施用量相對(duì)較高。整體模型中，農(nóng)民可支配收入沒(méi)有通過(guò)顯著性檢驗(yàn)。

表3 各影響因素的描述性評(píng)價(jià)

表4 整體模型和分組模型估計(jì)結(jié)果

注：*、**、***分別表示結(jié)果在10%、5%、1%的水平上顯著

Note: *, ** and *** indicate that the results are significant at the level of 10%, 5% and 1%, respectively.

對(duì)東部地區(qū)數(shù)據(jù)進(jìn)行回歸，其中農(nóng)民可支配收入在5%的水平上顯著且系數(shù)為負(fù)，說(shuō)明高收入的農(nóng)民有能力和意愿購(gòu)買更高效的化肥。金融規(guī)模在1%水平顯著且系數(shù)為正，說(shuō)明農(nóng)業(yè)金融大大促進(jìn)農(nóng)業(yè)技術(shù)提升。對(duì)西部地區(qū)數(shù)據(jù)進(jìn)行回歸，結(jié)果顯示，經(jīng)濟(jì)發(fā)展不顯著，且系數(shù)為負(fù)，與東部地區(qū)以及整體相反，可見(jiàn)西部地區(qū)的經(jīng)濟(jì)發(fā)展重心不在農(nóng)業(yè)，化肥施用與此關(guān)聯(lián)不大。農(nóng)業(yè)結(jié)構(gòu)在10%的水平顯著為正，說(shuō)明西部地區(qū)農(nóng)業(yè)產(chǎn)值占比越來(lái)越大。人口數(shù)量在5%的水平顯著且系數(shù)為負(fù)，說(shuō)明西部地區(qū)人口外流導(dǎo)致化肥對(duì)勞動(dòng)的替代從而增加化肥施用。金融規(guī)模在10%的水平顯著為正，說(shuō)明西部地區(qū)的農(nóng)業(yè)技術(shù)發(fā)展朝向仍為增產(chǎn)階段，較少考慮減肥問(wèn)題。綜上，東西部地區(qū)的回歸檢驗(yàn)結(jié)果較為一致，相當(dāng)于通過(guò)穩(wěn)健性檢驗(yàn)。

3 討 論

本文采用驅(qū)動(dòng)因素測(cè)度模型來(lái)測(cè)量化肥施用效率、區(qū)域施肥結(jié)構(gòu)、農(nóng)業(yè)技術(shù)發(fā)展、復(fù)種指數(shù)、耕地面積對(duì)化肥施用的作用方向與大小。解釋變量選取不夠全面，影響農(nóng)戶施肥行為的原因較復(fù)雜，自變量種類多，比如農(nóng)戶的施肥習(xí)慣或慣性、土地規(guī)模經(jīng)營(yíng)的比例等因素也會(huì)對(duì)化肥施用造成一定的影響，農(nóng)村勞動(dòng)力素質(zhì)水平相對(duì)較低，其固有的傳統(tǒng)施肥習(xí)慣難以改變，不利于科學(xué)施肥技術(shù)的推進(jìn)。而一定的條件下擴(kuò)大經(jīng)營(yíng)規(guī)模，推進(jìn)適度規(guī)模經(jīng)營(yíng)，有助于提高化肥施用效率。LMDI方法能夠防止模型外變量的干擾，只分解所用的解釋變量，這在一定程度上彌補(bǔ)了上述不足。如果過(guò)細(xì)考慮驅(qū)動(dòng)化肥施用的所有因素，會(huì)導(dǎo)致數(shù)據(jù)獲得性變低以及操作難度增大，模型結(jié)果缺乏針對(duì)性。本文得出化肥施用效率、農(nóng)業(yè)技術(shù)發(fā)展和復(fù)種指數(shù)為驅(qū)動(dòng)化肥施用量變化的主要因素在一定程度上較真實(shí)地反映出恩施州化肥施用環(huán)境。這不僅能夠給當(dāng)?shù)剞r(nóng)戶行為、政府政策制定提供借鑒意義，同時(shí)對(duì)各地市的化肥減量控制都有一定參考價(jià)值。

化肥施用量還存在大量深層次影響因素。比如區(qū)域空間差異，恩施州作為第一產(chǎn)業(yè)占比較高地區(qū)，其化肥施用量在湖北一直排名首位，但是州下各縣市又由于主要種植農(nóng)作物與地形等因素導(dǎo)致不同地區(qū)化肥施用同樣存在差異。比如本文時(shí)空差異部分提到，2018年鶴峰縣施肥強(qiáng)度有555.20 kg/hm2，來(lái)鳳縣僅為207.44 kg/hm2。同樣地，經(jīng)濟(jì)發(fā)展、金融規(guī)模等并不只與近年來(lái)發(fā)展勢(shì)頭迅猛的第二、三產(chǎn)業(yè)有關(guān)聯(lián)，它們?cè)谀撤N意義上與化肥施用量的多少也息息相關(guān)，有資料顯示，經(jīng)濟(jì)較發(fā)達(dá)的海南、北京與經(jīng)濟(jì)欠發(fā)達(dá)的青海、西藏、貴州等相比，其化肥施用量差異巨大[36]，本文正是采用實(shí)證分析法證明這些因素對(duì)化肥施用影響顯著?；适┯眠^(guò)量帶來(lái)大量污染，倡導(dǎo)減肥、綠色發(fā)展農(nóng)業(yè)刻不容緩。應(yīng)該從宏微觀方面綜合考慮，微觀層面上，根據(jù)不同地區(qū)的情況，因地制宜，宏觀層面上從深層次方面挖掘原因，有針對(duì)性地實(shí)施減肥增效，促使糧食生產(chǎn)和環(huán)境保護(hù)雙贏發(fā)展。

4 結(jié) 論

1）恩施州化肥施用在2007－2018年呈先上升后下降趨勢(shì)。隨著農(nóng)業(yè)經(jīng)濟(jì)的發(fā)展，化肥施用強(qiáng)度區(qū)域差異越來(lái)越大，2007年僅有巴東縣1個(gè)高強(qiáng)度區(qū)，2018年有巴東縣、建始縣、宣恩縣、咸豐縣4個(gè)高強(qiáng)度區(qū)和1個(gè)鶴峰縣超高強(qiáng)度區(qū)。2007、2012、2015、2018年這4 a的化肥施用最高區(qū)域與最低區(qū)域的強(qiáng)度比值呈先上升后平穩(wěn)的趨勢(shì)。

2）就恩施州整體而言，2016年恩施州化肥施用量比上年減少7 008.32 t，主要驅(qū)動(dòng)因素為化肥施用效率和區(qū)域結(jié)構(gòu)；2017年恩施州化肥施用量比上年減少3 341.65 t，主要驅(qū)動(dòng)因素為化肥施用效率和復(fù)種指數(shù)；2018年恩施州化肥施用量比上年減少15 059.81 t，主要驅(qū)動(dòng)因素為化肥施用效率和復(fù)種指數(shù)。就不同縣域而言，2015－2016年下降超1 000 t的有恩施市、利川市、建始縣和咸豐縣；2016－2017年下降超1 000 t的有恩施市、鶴峰縣；2017－2017年下降超1 000 t的有恩施市、建始縣、巴東縣和鶴峰縣。

3）整體上恩施州化肥施用量隨著經(jīng)濟(jì)的發(fā)展而增加，農(nóng)業(yè)結(jié)構(gòu)調(diào)整使農(nóng)業(yè)占比提高繼而增加化肥施用，金融規(guī)模擴(kuò)大對(duì)化肥施用同樣起正向作用，但隨著人口的增多，化肥施用量顯著減少。對(duì)東部地區(qū)而言，農(nóng)民可支配收入增多促使農(nóng)民購(gòu)買高效化肥從而減少化肥的施用，農(nóng)村金融規(guī)模的擴(kuò)大依然使化肥施用量增加，其余因素影響效應(yīng)不夠明顯。對(duì)西部地區(qū)而言，不同于恩施州整體與東部地區(qū)，經(jīng)濟(jì)發(fā)展雖然不顯著，但其對(duì)化肥施用產(chǎn)生負(fù)向影響，這是由于相對(duì)農(nóng)業(yè)來(lái)說(shuō)西部地區(qū)更注重其他產(chǎn)業(yè)的發(fā)展，農(nóng)業(yè)結(jié)構(gòu)與金融規(guī)模的變化對(duì)化肥施用產(chǎn)生正向影響，人口的增加會(huì)導(dǎo)致化肥施用量減少。

[1] 馬紅坤，曹原，毛世平. 歐盟共同農(nóng)業(yè)政策的綠色轉(zhuǎn)型軌跡及其對(duì)我國(guó)政策改革的鏡鑒[J]. 農(nóng)村經(jīng)濟(jì)，2019(3)：135-144.

Ma Hongkun, Cao Yuan, Mao Shiping. The green transformation track of EU common agricultural policy and its reference to China's policy reform[J]. Rural Economy, 2019(3): 135-144. (in Chinese with English abstract)

[2] 劉琦，趙明正. 農(nóng)業(yè)現(xiàn)代化進(jìn)程中農(nóng)業(yè)要素使用強(qiáng)度變化規(guī)律研究：基于全球29個(gè)主要農(nóng)業(yè)國(guó)家的國(guó)際經(jīng)驗(yàn)[J]. 農(nóng)業(yè)經(jīng)濟(jì)問(wèn)題，2018(3)：23-32.

Liu Qi, Zhao Mingzheng. Theoretical analysis and empirical study on international experience of agricultural modernization: Based on the panel data of 29 major agricultural countries[J]. Issues in Agricultural Economy, 2018(3): 23-32. (in Chinese with English abstract)

[3] 陳慧，高麗萍，廖慶喜，等. 肥料減量深施對(duì)土壤N2O排放和冬油菜產(chǎn)量的影響[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(21)：80-87.

Chen Hui, Gao Liping, Liao Qingxi, et al. Effects of reduced and deep fertilizer on soil N2O emission and yield of winter rapeseed[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(21): 80-87. (in Chinese with English abstract)

[4] 展曉瑩，張愛(ài)平，張晴雯. 農(nóng)業(yè)綠色高質(zhì)量發(fā)展期面源污染治理的思考與實(shí)踐[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(20)：1-7.

Zhan Xiaoying, Zhang Aiping, Zhang Qingwen. Controlling agricultural non-point source pollution: Thinking andpractice in the era of agricultural green high quality development[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(20): 1-7. (in Chinese with English abstract)

[5] 金書(shū)秦，牛坤玉，韓冬梅. 農(nóng)業(yè)綠色發(fā)展路徑及其“十四五”取向[J]. 改革，2020(2)：30-39.

Jin Shuqin, Niu Kunyu, Han Dongmei. The Path of agricultural green development and its orientation in the 14th five-year plan period[J]. Reform, 2020(2): 30-39. (in Chinese with English abstract)

[6] Green D A G, Ng'ong ‘ola D H. Factors affecting fertilizer adoption in less developed countries: An application of multivariate logistic analysis in Mala?i[J]. Journal of Agricultural Economics, 1993, 44(1): 99-109.

[7] Endale K. Fertilizer consumption and agricultural productivity in Ethiopia[R/OL]. Working Papers 003, Ethiopian Development Research Institute, 2011. [2021-01-10 ]https://ideas.repec.org/p/etd/wpaper/003.html.

[8] Burke W J, Jayne T S, Black J R. Factors explaining the low and variable profitability of fertilizer application to maize in Zambia[J]. Agricultural Economics, 2017, 48(1): 115-126.

[9] 劉莉，劉靜. 基于種植結(jié)構(gòu)調(diào)整視角的化肥減施對(duì)策研究[J]. 中國(guó)農(nóng)業(yè)資源與區(qū)劃，2019，40(1)：17-25.

Liu Li, Liu Jing. Study on the path of chemical fertilizer reduction from the perspective of planting structure adjustment[J]. Chinese Journal of Agricultural Resources and Regional Planning, 2019, 40(1): 17-25. (in Chinese with English abstract)

[10] 鄭微微，徐雪高. 江蘇省化肥施用強(qiáng)度變化驅(qū)動(dòng)因子分解及其影響因素分析[J]. 華中農(nóng)業(yè)大學(xué)學(xué)報(bào)：社會(huì)科學(xué)版，2017(4)：55-62，147.

Zheng Weiwei, Xu Xuegao. Analysis on the driving factor decomposition and its influencing factors of fertilizer application intensity change in Jiangsu province[J]. Journal of Huazhong Agricultural University: Social Sciences Edition, 2017(4): 55-62, 147. (in Chinese with English abstract)

[11] 龔琦，王雅鵬. 我國(guó)農(nóng)用化肥施用的影響因素：基于省際面板數(shù)據(jù)的實(shí)證分析[J]. 生態(tài)經(jīng)濟(jì)，2011(2)：33-38，43.

Gong Qi, Wang Yapeng. Influential factors of chemical fertilizer employing in China: Empirical study based on provincial panel data[J]. Ecological Economy, 2011(2): 33-38, 43. (in Chinese with English abstract)

[12] 王珊珊，張廣勝，李秋丹，等. 我國(guó)糧食主產(chǎn)區(qū)化肥施用量增長(zhǎng)的驅(qū)動(dòng)因素分解[J]. 農(nóng)業(yè)現(xiàn)代化研究，2017，38(4)：658-665.

Wang Shanshan, Zhang Guangsheng, Li Qiudan, et al. Driving factor decomposition of fertilizer application growth in China’s main grain producing areas[J]. Research of Agricultural Modernization, 2017, 38(4): 658-665. (in Chinese with English abstract)

[13] 欒江，仇煥廣，井月，等. 我國(guó)化肥施用量持續(xù)增長(zhǎng)的原因分解及趨勢(shì)預(yù)測(cè)[J]. 自然資源學(xué)報(bào)，2013，28(11)：1869-1878.

Luan Jiang, Qiu Huangunag, Jing Yue, et al. Decomposition of factors contributed to the increase of China’s chemical fertilizer use and projections for future fertilizer use in China[J]. Journal of Natural Resources, 2013, 28(11): 1869-1878. (in Chinese with English abstract)

[14] 趙明正，趙翠萍，李天祥，等.“零增長(zhǎng)”行動(dòng)背景下中國(guó)化肥使用量下降的驅(qū)動(dòng)因素研究：基于LMDI分解和面板回歸分析[J]. 農(nóng)業(yè)技術(shù)經(jīng)濟(jì)，2019(12)：118-130.

Zhao Mingzheng, Zhao Cuiping, Li Tianxiang, et al. Driving forces analysis for China’s recent fertilizer use reduction under the “Zero Growth ”action plan: Estimation based on the LMDI approach and panel regression analysis[J]. Journal of Agrotechnical Economics, 2019(12): 118-130. (in Chinese with English abstract)

[15] Ang B W. The LMDI approach to decomposition analysis: A practical guide[J]. Energy Policy, 2005, 33: 867-871.

[16] 曹銳，錢(qián)海婷. 新型城鎮(zhèn)化發(fā)展視角下農(nóng)村人力資源現(xiàn)狀與發(fā)展對(duì)策研究[J]. 農(nóng)業(yè)經(jīng)濟(jì)，2020(3)：65-67.

Cao Rui, Qian Haiting. Research on the current situation and development countermeasures of rural human resources from the perspective of new-type urbanization development[J]. Agricultural Economy, 2020(3): 65-67. (in Chinese with English abstract)

[17] 申健，常慶瑞，李粉玲，等. 2000－2013年關(guān)中地區(qū)耕地復(fù)種指數(shù)遙感動(dòng)態(tài)監(jiān)測(cè)[J]. 農(nóng)業(yè)機(jī)械學(xué)報(bào)，2016，47(8)：280-287.

Shen Jian, Chang Qingrui, Li Fenling, et al. Dynamic monitoring of cropping index in Guanzhong area using remote sensing in 2000-2013[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(8): 280-287. (in Chinese with English abstract)

[18] 姚成勝，黃琳，呂晞，等. 基于能值理論的中國(guó)耕地利用集約度時(shí)空變化分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2014，30(8)：1-12.

Yao Chengsheng, Huang Lin, Lyu Xi, et al. Temporal and spatial change of cultivated land use intensity in China based on emergy theory[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2014, 30(8): 1-12. (in Chinese with English abstract)

[19] 孫亮，商蕾. 巴西農(nóng)業(yè)發(fā)展及其現(xiàn)代化[J]. 世界農(nóng)業(yè)，2014(5)：160-162.

Sun Liang, Shang Lei. Agricultural development and modernization in Brazil[J]. World Agriculture, 2014(5): 160-162. (in Chinese with English abstract)

[20] 尹峻，黨敬淇，孫小霞. 價(jià)格與替代效應(yīng)對(duì)農(nóng)業(yè)生產(chǎn)中化肥施用的影響：以玉米和圓白菜為例[J]. 中國(guó)農(nóng)業(yè)大學(xué)學(xué)報(bào)，2020，25(5)：209-220.

Yin Jun, Dang Jingqi, Sun Xiaoxia. Effect of price and substitution on fertilizer use in agricultural production:An empirical analysis of corn and cabbage[J]. Journal of China Agricultural University, 2020, 25(5): 209-220. (in Chinese with English abstract)

[21] 李志平. 促進(jìn)農(nóng)村地區(qū)技術(shù)進(jìn)步的模塊化金融創(chuàng)新研究：以湖北為例[J]. 科技管理研究，2010，30(14)：100-103.

Li Zhiping. Module financial innovation for technological advances: A case study of Hubei rural area[J]. Science and Technology Management Research, 2010, 30(14): 100-103. (in Chinese with English abstract)

[22] 張惠茹. 農(nóng)村金融市場(chǎng)靈活性缺口：理論與實(shí)證分析[J]. 農(nóng)村經(jīng)濟(jì)，2012(7)：46-50.

Zhang Huiru. Ruralfinancial market flexibility gap: Theoretical and empirical analysis[J]. Rural Economy, 2012(7): 46-50. (in Chinese with English abstract)

[23] 宋大平，左強(qiáng)，劉本生，等.農(nóng)業(yè)面源污染中氮排放時(shí)空變化及其健康風(fēng)險(xiǎn)評(píng)價(jià)研究：以淮河流域?yàn)槔齕J]. 農(nóng)業(yè)環(huán)境科學(xué)學(xué)報(bào)，2018，37(6)：1219-1231.

Song Daping, Zuo Qiang, Liu Bensheng, et al. Estimation of spatio-temporal variability and health risks of nitrogen emissions from agricultural non-pointsource pollution：A case study of the Huaihe River Basin, China[J]. Journal of Agro-Environment Science, 2018, 37(6): 1219-1231. (in Chinese with English abstract)

[24] 趙建勛，程燚. 安徽省肥料施用現(xiàn)狀與對(duì)策[J]. 安徽農(nóng)學(xué)通報(bào)(上半月刊)，2011，17(1)：102-103，131.

Zhao Jianxun, Cheng Yi. Current situation and countermeasure of fertilizer application in Anhui province[J]. Anhui Agricultural Science Bulletin, 2011, 17(1): 102-103,131. (in Chinese with English abstract)

[25] 宋秀杰，程大軍，張?chǎng)?，? 北京種植業(yè)結(jié)構(gòu)調(diào)整及化肥面源污染控制[A]. 中國(guó)環(huán)境科學(xué)學(xué)會(huì). 2010中國(guó)環(huán)境科學(xué)學(xué)會(huì)學(xué)術(shù)年會(huì)論文集(第一卷) [C]. 中國(guó)環(huán)境科學(xué)學(xué)會(huì)，2010.

[26] Motesharezadeh B, Etesami H, Bagheri-Novair S, et al. Fertilizer consumption trend in developing countries vs. developed countries[J/OL]. Environmental Monitoring and Assessment, 2017, 189(3): 103.

[27] 環(huán)境保護(hù)部. 關(guān)于印發(fā)《國(guó)家級(jí)生態(tài)鄉(xiāng)鎮(zhèn)申報(bào)及管理規(guī)定（試行）》的通知(環(huán)發(fā) [2010]75) [EB/OL]. (2010-06-23) [2019-02-20].http://www.jhepb.gov.cn/html/news/20150818/15303.shtml.

[28] 張福鎖，王激清，張衛(wèi)峰，等. 中國(guó)主要糧食作物肥料利用率現(xiàn)狀與提高途徑[J]. 土壤學(xué)報(bào)，2008，45(5)：915-924.

Zhang Fusuo, Wang Jiqing, Zhang Weifeng, et al. Current situation and improving ways of fertilizer utilization rate of main grain crops in China[J]. Acta Pedologica Sinica, 2008, 45(5): 915-924. (in Chinese with English abstract)

[29] 農(nóng)業(yè)部.農(nóng)業(yè)部關(guān)于印發(fā)《到2020年化肥使用量零增長(zhǎng)行動(dòng)方案》的通知[J]. 中華人民共和國(guó)農(nóng)業(yè)部公報(bào)，2015(3)：19-27.

[30] 王萍萍，韓一軍，張益. 中國(guó)農(nóng)業(yè)化肥施用技術(shù)效率演變特征及影響因素[J]. 資源科學(xué)，2020，42(9)：1764-1776.

Wang Pingping, Han Yijun, Zhang Yi. Characteristics of change and influencing factors of the technical efficiency of chemical fertilizer use for agricultural production in China[J]. Resources Science, 2020, 42(9): 1764-1776. (in Chinese with English abstract)

[31] 劉小玄，鄭京海. 國(guó)有企業(yè)效率的決定因素：1985-1994[J]. 經(jīng)濟(jì)研究，1998(1)：39-48.

Liu Xiaoxuan, Zheng Jinghai. Determinants of the efficiency of state owned enterprises: 1985-1994[J]. Economic Research Journal, 1998(1): 39-48. (in Chinese with English abstract)

[32] 張典香. 建始縣農(nóng)業(yè)技術(shù)推廣體系建設(shè)的基本思路與對(duì)策[J]. 中國(guó)農(nóng)業(yè)信息，2015(5)：124.

Zhang Dianxiang. The basic idea and countermeasure of construction of agricultural technology extension system in Jianshi county[J]. China Agricultural Informatics, 2015(5): 124. (in Chinese with English abstract)

[33] 江方波. 建始縣農(nóng)業(yè)特色優(yōu)質(zhì)資源產(chǎn)業(yè)化發(fā)展思考[J]. 現(xiàn)代農(nóng)村科技，2015(19)：5-6.

Jiang Fangbo. Reflection on industrialization development of agricultural characteristic high quality resources in Jianshi county[J]. Modern Rural Science and Technology, 2015(19): 5-6. (in Chinese with English abstract)

[34] 曹銀貴，周偉，王靜，等. 基于主成分分析與層次分析的三峽庫(kù)區(qū)耕地集約利用對(duì)比[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(4)：291-296.

Cao Yingui, Zhou Wei, Wang Jing, et al. Comparative on regional cultivated land intensive use based on principal component analysis and analytic hierarchy process in Three Gorges Reservoir Area[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(4): 291-296. (in Chinese with English abstract)

[35] 洪凱，郭發(fā)吉，覃信立，等. 咸豐縣玉米高產(chǎn)高效模式探究[J]. 現(xiàn)代農(nóng)業(yè)科技，2014(2)：85-88.

Hong Kai, Guo Faji, Qin Xinli, et al. Study on high yield and high efficiency model of maize in Xianfeng county[J]. Modern Agricultural Science and Technology, 2014(2): 85-88. (inChinese with English abstract)

[36] 劉欽普.中國(guó)化肥施用強(qiáng)度及環(huán)境安全閾值時(shí)空變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017，33(6)：214-221.

Liu Qinpu. Spatio-temporal changes of fertilization intensity and environmental safety threshold in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017, 33(6): 214-221. (in Chinese with English abstract)

Analysis of driving factors of chemical fertilizer application by LMDI method at Enshi of Hubei Province in China

Mao Zhongming, Yang Haoyi, Liu Chensiyu, Ma Zhu

(,,430074,)

Excessive application of chemical fertilizer has caused the ever-increasing agricultural non-point source pollution in recent years. Great efforts were made to control agricultural pollution, and thereby explore the changing characteristics and driving factors for the application intensity of chemical fertilizer. The amount of chemical fertilizer was widely expected to be reduced, further to enhance the efficiency of chemical fertilizer application and green development of agriculture. Taking the Enshi Prefecture of Hubei Province in central China as an empirical case, this study aims to determine the spatio-temporal variation characteristics and the driving force of chemical fertilizer application. Logarithmic Mean Divisia Index (EMDI) exponential decomposition and panel regression models were established to comprehensively consider the efficiency of fertilizer application, regional fertilization structure, agricultural technology development, multiple planting index, and agricultural acreage. The results showed as follows: 1) The fertilizer application in the study area increased firstly, then steadily, and finally decreased over the last ten years. The ratio of nitrogen, phosphorus and potassium increased from 1:0.30:0.15 in 2007 to 1:0.39:0.26 in 2018, gradually approaching the level of developed countries in the world. There were one low-intensity zone, two medium-intensity zones, four high-intensity zones, and one ultra-high intensity zone by 2018 for chemical fertilizer application in the study area. 2) Chemical fertilizer application in the study area decreased by 7 008.32 t in 2016, 3 341.65 t in 2017, and 15 059.81 t in 2018, compared with the previous year. The efficiency increase of fertilizer application reduced the amount of fertilizer application by 11 159.44, 14 051.95 and 8 644.21 t in the last three years. The multiple cropping index reduced the amount of fertilizer application by 5 229.76 t in 2017, and 30 889.98 t in 2018. There was no significant effect value in 2016. The development of agricultural technology increased the amount of fertilizer application by 4 089.71, 11 136.97 and 12 215.74 t in 2016, 2017 and 2018. But there was no significant driving effect of regional structure. 3) There were positive effects of three factors to reduce the application amount of chemical fertilizer, including the economic development, adjustment of agricultural structure, and expansion of financial scale. But the population increase presented negative effects in this case. It is suggested to change the direction of agricultural technology development, from merely increasing production to multi production with high efficiency. The effective arable land area should also be reasonably adjusted, in order to further improve the efficiency of chemical fertilizer application. A model of ecological security and green development should be fostered to promote the reduction of chemical fertilizers and the increase of efficiency in modern agriculture.

agriculture; fertilizers; spatio-temporal characteristics; LMDI model; driving factors; Enshi prefecture

毛中明，楊豪毅，劉陳思雨，等. 基于LMDI方法的恩施州化肥施用驅(qū)動(dòng)因素分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(13)：111-120.

10.11975/j.issn.1002-6819.2021.13.013 http://www.tcsae.org

Mao Zhongming, Yang Haoyi, Liu Chensiyu, et al. Analysis of driving factors of chemical fertilizer application by LMDI method at Enshi of Hubei Province in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(13): 111-120. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2021.13.013 http://www.tcsae.org

2021-02-24

2021-05-10

湖北省社會(huì)科學(xué)基金項(xiàng)目（HBSK2019ZXYB036）；中國(guó)科學(xué)院--國(guó)家民委農(nóng)業(yè)信息技術(shù)研究與開(kāi)發(fā)聯(lián)合實(shí)驗(yàn)室2019年度招標(biāo)課題（PJW060011906）

毛中明，博士，教授，研究方向?yàn)楫a(chǎn)業(yè)組織理論。Email：782673444@qq.com

10.11975/j.issn.1002-6819.2021.13.013

F327

A

1002-6819(2021)-13-0111-10