農(nóng)業(yè)機械投入對化肥減量化的作用

1、農(nóng)業(yè)機械投入對化肥減量化的作用 基于全國232個地級市數(shù)據(jù)，運用有限混合模型FMM分析了農(nóng)業(yè)機械投入對化肥減量化的作用效果。結果說明：總樣本由農(nóng)業(yè)機械化低水平樣本和高水平樣本組成；機械強度與化肥強度之間存在非線性關系；在機械化低水平階段，機械強度與化肥強度呈正相關關系；在機械化高水平階段，機械強度與化肥強度呈負相關關系；要通過增加農(nóng)業(yè)機械投入促進化肥減量化，關鍵在于轉變農(nóng)戶施肥需求。建議加快研發(fā)推廣闊型農(nóng)業(yè)機械，推廣秸稈還田、化肥機械深施等技術，改善土壤理化性質(zhì)，提高化肥利用率，實現(xiàn)化肥減量化?；蕼p量化；農(nóng)業(yè)機械化；施肥方式；作用效果；有限混合模型近年來，化肥過度施用造成了資源浪費、環(huán)境污染

2、等一系列問題，制約了農(nóng)業(yè)可持續(xù)開展。研究農(nóng)業(yè)化肥減量化的有效途徑具有重要現(xiàn)實意義。為了加快實現(xiàn)化肥減量化目標，農(nóng)業(yè)農(nóng)村部于2021年2月發(fā)布了?到2021年化肥使用量零增長行動方案?，為化肥減量化工作提供了政策條件。同時，秸稈還田、化肥機械深施、種肥同播等技術的大型農(nóng)業(yè)機械相繼出現(xiàn)，為農(nóng)業(yè)化肥施用方式的轉變提供了技術條件。由此，增加農(nóng)業(yè)機械投入是否能促進傳統(tǒng)化肥施用方式轉向精細化施肥進而實現(xiàn)化肥減量化？針對該問題的思考，有助于為實現(xiàn)化肥減量化提供有益參考。1文獻綜述目前，對于我國化肥減量化實現(xiàn)路徑的研究仍較少，已有研究大多數(shù)為定性研究，從化肥施用對象、施用總量等角度分析總結我國化肥施用現(xiàn)狀，根

3、據(jù)當前不同作物的化肥施用強度提出減量化的作物對象優(yōu)先順序，也有學者關注中國化肥減量化的目標值，基于農(nóng)產(chǎn)品供應和水資源壓力的雙重約束設定中國化肥減量化總體目標。對于農(nóng)戶化肥施用行為影響因素的研究，主要集中在農(nóng)戶的認知風險、個人特征、家庭特征層面，大多直接采用傳統(tǒng)線性回歸分析進行影響因素識別，或從結構驅(qū)動效應和效率驅(qū)動效應的角度對化肥施用強度的變化進行分解；有學者研究了水稻種植這一單一作物生產(chǎn)類型中的農(nóng)業(yè)機械化技術對稻農(nóng)化肥減量化施用行為的影響，認為生產(chǎn)方式的轉變是實現(xiàn)化肥減量化的關鍵途徑。關于機械化開展帶來的系列影響，從已有成果來看，主要集中在機械化與農(nóng)業(yè)總產(chǎn)出、糧食產(chǎn)出、生產(chǎn)效率、糧食播種面積

4、、農(nóng)業(yè)生產(chǎn)中的勞動力要素替代等方面，認為農(nóng)業(yè)機械化開展能夠增加農(nóng)業(yè)產(chǎn)出，對糧食生產(chǎn)存在正的空間溢出效應，對農(nóng)業(yè)經(jīng)濟開展方式轉變具有正向影響，但對糧食播種面積沒有顯著影響。少量研究從定性角度關注了機械化與生態(tài)環(huán)境的關系，認為農(nóng)業(yè)機械化開展與農(nóng)業(yè)資源利用效率提高、農(nóng)村環(huán)境改善和農(nóng)業(yè)可持續(xù)開展密切相關，應大力研發(fā)節(jié)能、低碳農(nóng)業(yè)機械，開展精準耕整、播種、施肥等技術，但仍未明確涉及對化肥減量化影響的作用效果?？傮w而言，關于化肥減量化實現(xiàn)路徑的研究仍處于起步階段，存在較大空白，鮮有針對農(nóng)業(yè)機械投入對于化肥減量化效果的討論，因此，及時展開這方面的研究，除了為農(nóng)業(yè)生產(chǎn)安排和政策制定提供參考，對于補充新的研究案

5、例也具有一定意義。為分析農(nóng)業(yè)機械投入能否促進化肥減量化，本文假設一個地區(qū)的農(nóng)業(yè)機械化水平開展具有階段性，并且不同的階段存在一定概率混合，在不同階段農(nóng)戶的化肥施用方式和施用需求不同，農(nóng)業(yè)機械投入與化肥施用存在多階段演變關系。在此根底上，以2021年232個地級市數(shù)據(jù)為樣本，運用有限混合模型FMM識別我國農(nóng)業(yè)機械化開展水平的異質(zhì)性構成，分析農(nóng)業(yè)機械投入與化肥減量化兩者間關系的階段性變化，為加快實現(xiàn)我國化肥減量化目標提供經(jīng)驗證據(jù)和政策建議。2模型與方法化肥施用屬于農(nóng)業(yè)生產(chǎn)中的投入行為，經(jīng)濟學上可用追求利潤最大化的農(nóng)戶生產(chǎn)模型來描述農(nóng)戶選擇行為：在保證同等產(chǎn)出的情況下投入最少生產(chǎn)要素。然而農(nóng)戶的生產(chǎn)投

6、入受到給定的技術條件限制，不同的機械水平帶來耕作模式變化，也導致化肥投入方式發(fā)生變化。由于農(nóng)業(yè)機械化開展水平的階段性，本文認為農(nóng)業(yè)機械投入與化肥施用兩者之間存在多階段演變關系。農(nóng)戶層面，農(nóng)戶施肥方式在不同農(nóng)業(yè)機械化水平階段有所不同。在農(nóng)業(yè)機械化低水平階段，農(nóng)戶使用小型簡易機械進行表土疏松、作物收割，化肥施用方式以撒施、表施為主，屬于傳統(tǒng)施肥方式，養(yǎng)分損失較大，化肥利用率較低。在此階段，隨著機械投入的增加，農(nóng)戶為獲得最大利潤，提高化肥強度，直至化肥投入的邊際收益為零。隨著農(nóng)業(yè)機械化開展，進入農(nóng)業(yè)機械大型化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，農(nóng)戶可采用秸稈還田、土壤深耕、種肥同播、根施肥料等技術，能夠改

7、進土壤營養(yǎng)特性，減少肥料流失，提高化肥利用率。因此，在農(nóng)機大型化高水平階段，隨著機械投入的增加，施肥方式轉變，為獲得利潤最大化，農(nóng)戶所需要投入的化肥更少，化肥強度隨機械強度提高而降低。然而農(nóng)業(yè)生產(chǎn)與機械化開展水平往往情況多樣，較為復雜。不同農(nóng)戶的農(nóng)業(yè)機械投入不同，其施肥方式不同，其開展轉型的階段并不能清晰劃分。因此，可將一個地市的化肥施用情況理解為不同機械化水平施肥行為的混合：農(nóng)業(yè)機械化低水平的傳統(tǒng)施肥行為和農(nóng)業(yè)機械化高水平的施肥行為按照一定概率比例的混合，且每一類行為都有各自的分布。如果一個地市機械化低水平傳統(tǒng)施肥行為的概率比例大于機械化高水平施肥行為的概率比例，那么該市處于農(nóng)業(yè)機械化低水平

8、階段；反之，處于高水平階段。同理，一個地市層面的樣本，是由處于不同機械化水平階段的地市按照一定概率比例的混合，各階段地市服從各自分布，地市之間具有異質(zhì)性。3實證分析本文收集了2021年232個地市數(shù)據(jù)，通過有限混合模型，定量考察在不同農(nóng)業(yè)機械化水平階段，農(nóng)業(yè)機械投入與化肥減量化之間的關系，以期揭示農(nóng)業(yè)機械投入水平與化肥減量化之間的關系。本文被解釋變量選取化肥強度，即單位農(nóng)作物播種面積對應的化肥折純使用量。關鍵解釋變量選取機械強度，即單位農(nóng)作物播種面積對應的農(nóng)業(yè)機械總動力投入。根據(jù)農(nóng)戶生產(chǎn)模型，化肥施用和機械投入還受到種植經(jīng)濟結構、城鎮(zhèn)化程度和人力資本的影響，因此本文將以上變量納入控制變量。對于

9、種植經(jīng)濟結構，經(jīng)濟作物單位耕種面積帶來的經(jīng)濟效益普遍高于糧食作物，選取單位耕種面積對應的種植業(yè)產(chǎn)值作為種植經(jīng)濟結構的代理變量。對于城鎮(zhèn)化程度，選取常住城鎮(zhèn)人口占常住總人口比例作為城鎮(zhèn)化率。對于人力資本，采用每萬人在校學生數(shù)代理人力資本。根據(jù)2021年2月開始實施?到2021年化肥使用量零增長行動方案?，本文以此作為化肥減量化工作標志性時間節(jié)點。本文通過各省統(tǒng)計年鑒收集了2021年全國232個地級市數(shù)據(jù)，涵蓋糧食主產(chǎn)省份中的內(nèi)蒙古、吉林、安徽、山東、江蘇、河北、河南、湖南、黑龍江和非主產(chǎn)省份中的廣東、浙江、海南、廣西、-、甘肅、貴州、陜西。表1顯示了232個地級市5個指標的描述性統(tǒng)計結果，圖1顯

10、示了總樣本中機械強度的概率密度分布。對于機械強度，總樣本分布并不完全服從正態(tài)分布，考慮存在至少兩個分布的疊加。為進行對照，先對包含232個地市的總樣本進行傳統(tǒng)OLS線性回歸，系數(shù)估計結果見表2，殘差分布見圖2。傳統(tǒng)OLS線性回歸顯示，在其他條件保持不變的情況下，隨著農(nóng)業(yè)機械投入增加，化肥強度下降，但不具有統(tǒng)計顯著性，且殘差分布不符合正態(tài)分布，總樣本分布存在起碼兩個分布的疊加。為此，有必要采用FMM作進一步分析。總樣本可能由兩個及以上子分布混合組成，但當子分布數(shù)目較多時，會增加有限混合模型估計的難度，且本文數(shù)據(jù)中，當子分布數(shù)目大于兩個時，估計結果并不顯著。結合機械化水平存在兩階段特征，本文選擇子

11、分布數(shù)目為兩個。FMM回歸結果顯示見表2，總樣本分布由子分布一占0.779，子分布二占0.221的概率比例混合組成，且子分布一的機械強度均值7.652小于子分布二9.618，故子分布一對應農(nóng)業(yè)機械化低水平階段，子分布二對應農(nóng)業(yè)機械化高水平階段。在農(nóng)業(yè)機械化低水平階段，農(nóng)戶使用低水平農(nóng)業(yè)機械，增加低水平機械投入未能轉變農(nóng)戶的傳統(tǒng)施肥方式，化肥強度與機械強度呈正相關關系。從邊際效應看，在保持其他因素不變的情況下，農(nóng)業(yè)機械強度每提高1個單位kW/hm2，化肥強度提高10.070個單位kg/hm2，并在1%的水平下具有統(tǒng)計顯著性。在農(nóng)業(yè)機械化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，帶來秸稈還田、化肥機械深施技術

12、的改進，有助于改進土壤肥力，轉變農(nóng)戶的施肥需求，化肥強度與機械強度呈負相關關系。從邊際效應看，在保持其他因素不變的情況下，農(nóng)業(yè)機械強度每提高1個單位kW/hm2，化肥強度降低9.306個單位kg/hm2。實證說明，在農(nóng)業(yè)機械化開展過程中，機械強度與化肥強度之間存在非線性兩階段演變關系，增加農(nóng)業(yè)機械投入能否促進化肥減量化，關鍵在于能否轉變農(nóng)戶的施肥需求。在機械化低水平階段，機械投入以簡易型機械為主，可實現(xiàn)的技術層次較低，如簡單的表土翻耕、收割，并未轉變農(nóng)戶施肥方式，機械強度與化肥強度呈正相關關系。在機械化高水平階段，大型農(nóng)業(yè)機械出現(xiàn)，秸稈還田、化肥深施、精準根施等技術得到應用，化肥利用率提高，化

13、肥施用需求下降，機械強度與化肥強度呈負相關關系。為減少農(nóng)業(yè)種植經(jīng)濟結構不同導致化肥強度存在差異而對實證結果造成影響，本文對糧食主產(chǎn)省份的123個地市進行FMM分析見表3。結果說明，在機械化低水平階段，機械強度與化肥強度之間呈正相關關系，從邊際效應看，保持其他因素不變，農(nóng)業(yè)機械強度每增加1個單位kW/hm2，化肥強度增加17.253個單位kg/hm2；在機械化高水平階段，機械強度與化肥強度呈負相關關系，從邊際效應看，保持其他因素不變，農(nóng)業(yè)機械強度增加1個單位kW/hm2，化肥強度減少9.364個單位kg/hm2，與理論設想一致。4研究結論本文假設機械化開展水平具有階段性，由于機械化水平開展會導致

14、農(nóng)戶化肥施用方式發(fā)生轉變，然而由于實際情況的復雜性，一個地區(qū)的農(nóng)業(yè)生產(chǎn)是農(nóng)業(yè)機械化高水平階段生產(chǎn)行為和低水平階段生產(chǎn)行為的概率混合；在機械化低水平階段，農(nóng)戶沿用傳統(tǒng)施肥方式；在機械化高水平階段，農(nóng)戶逐漸轉向精細化施肥。增加農(nóng)業(yè)機械投入對于化肥減量化的作用是非線性的兩階段演變過程。通過有限混合模型對232個地市組成的總樣本進行分類，分別探討在不同農(nóng)業(yè)機械化水平下農(nóng)業(yè)機械投入對化肥減量化的作用效果，結果說明：總樣本由77.9%的農(nóng)業(yè)機械化低水平階段樣本和22.1%的機械化高水平階段樣本組成；在機械化低水平階段，機械強度與化肥強度呈正相關關系，保持其他因素不變，機械強度提高1個單位kW/hm2，化肥強度提高10.070個單位kg/hm2；在機械化高水平階段，機械強度與化肥強度呈負相關關系，保持其他因素不變，農(nóng)業(yè)機