摘要 基于2004—2018年三大玉米種植區(qū)的面板數(shù)據(jù)，運用SBM模型方向性距離函數(shù)中的Malmquist-luenberger生產(chǎn)率指數(shù)和以產(chǎn)量為導(dǎo)向的DEA-malmquist指數(shù)，實證研究我國玉米種區(qū)生態(tài)效率的區(qū)域差異，以及在環(huán)境制約下我國玉米不同種區(qū)的真實生產(chǎn)狀況。隨后通過Tobit面板模型，分析不同玉米種區(qū)的生態(tài)效率和生產(chǎn)效率指數(shù)在相同影響因素下的差異。結(jié)果表明，我國玉米種區(qū)整體為粗放式的發(fā)展模式，北部種區(qū)整體生態(tài)效率較高，生產(chǎn)方式較為低碳;中部種區(qū)的平均生態(tài)效率較高，但技術(shù)效率退步明顯;西部種區(qū)的生態(tài)效率增長緩慢，技術(shù)投入不足。受災(zāi)率、生態(tài)治理水平和環(huán)保投入因素對各玉米種區(qū)的生態(tài)效率和傳統(tǒng)M指數(shù)存在趨同的影響作用，但不同種植區(qū)域還是存在較大的影響差異。

關(guān)鍵詞 玉米;SBM模型;地區(qū)差異;影響因素

中圖分類號 S-9;F326.11? 文獻標識碼 A? 文章編號 0517-6611（2022）05-0196-04

doi：10.3969/j.issn.0517-6611.2022.05.050

開放科學(xué)（資源服務(wù)）標識碼（OSID）：

Study on Regional Differences and Influencing Factors of Ecological Efficiency of Maize Planting Areas in China

CHEN? Yan-jun

（Economics and Management School of Yangtze University，Jingzhou，Hubei 434023）

Abstract Based on the panel data of three corn growing areas from 2004 to 2018， this paper used the Malmquist-luenberger productivity index and yield oriented DEA-malmquist index in the directional distance function of SBM model to empirically study the regional differences of ecological efficiency of corn growing areas in China， and the real production situation of different corn growing areas under environmental constraints. Then， through Tobit panel model， the differences of eco efficiency and production efficiency index of different maize seed areas under the same influence factors were analyzed. The results showed that the overall development mode of maize seed area in China was extensive， the overall ecological efficiency of northern seed area was higher， and the production mode was low-carbon;the average ecological efficiency of central seed area was higher， but the technical efficiency declined significantly;the ecological efficiency of western seed area increased slowly， and the technical input was insufficient. The disaster rate， ecological management level and environmental protection input factors have a convergence effect on the ecological efficiency and traditional M index of maize planting areas， but there are still large differences in different planting areas.

Key words Maize;SBM model;Regional difference;Influencing factors

基金項目 國家社科基金項目“以生態(tài)安全和糧食安全并重為導(dǎo)向的糧食補貼政策研究”（17BGL250）。

作者簡介 陳衍?。?994—），男，湖北武漢人，碩士研究生，研究方向：農(nóng)業(yè)經(jīng)濟理論與政策。

收稿日期 2021-06-04

玉米作為我國的三大主糧之一，是飼料用糧和工業(yè)用糧的重要來源。伴隨玉米高速增產(chǎn)，以化肥農(nóng)藥為主的農(nóng)業(yè)消耗品也會隨之大量投入，糧食高產(chǎn)的背后是農(nóng)業(yè)生態(tài)資源的損失。據(jù)國家統(tǒng)計局公布的數(shù)據(jù)顯示，自臨時收儲政策實施以來，我國玉米種植平均化肥折純量至2015年增幅達19.12%，特別是以東北三省為代表的北部種區(qū)，農(nóng)用消耗品施用量占全國總施用量的40.96%，增幅近10%[1]。雖然目前國家對玉米實行價補分離政策，對農(nóng)業(yè)供給側(cè)調(diào)減玉米種植面積，倡導(dǎo)糧食產(chǎn)業(yè)以綠色低碳的生產(chǎn)方式來提升單產(chǎn)和生產(chǎn)效率。但最近對糧食產(chǎn)業(yè)的研究[2]發(fā)現(xiàn)，我國玉米各種區(qū)依舊存在生產(chǎn)經(jīng)營粗放、種植占比過大等問題，與我國糧食產(chǎn)業(yè)要發(fā)展協(xié)調(diào)好資源消耗、高效生產(chǎn)和環(huán)境保護的這一目標仍存在差距。

農(nóng)業(yè)生態(tài)效率是衡量農(nóng)業(yè)生態(tài)化發(fā)展水平的測度標準，其目標是為實現(xiàn)在減少環(huán)境影響下同時創(chuàng)造農(nóng)業(yè)生產(chǎn)活動的總價值。隨著環(huán)境問題在糧食產(chǎn)業(yè)中越來越受到學(xué)者的關(guān)注，于婷等[3]使用SBM-undesirable 模型，在碳排放等環(huán)境條件約束下測度不同糧食主產(chǎn)區(qū)的生態(tài)效率，發(fā)現(xiàn)不同區(qū)域間的農(nóng)業(yè)生態(tài)效率存在較大差異，且呈不斷擴大趨勢;魯慶堯等[4]通過SBM模型分析發(fā)現(xiàn)我國糧食生產(chǎn)的三大區(qū)域生態(tài)效率區(qū)域性差異明顯，呈由西部、中部到東部顯著的遞減趨勢。較多學(xué)者對我國糧食產(chǎn)業(yè)整體以及糧食主產(chǎn)區(qū)的發(fā)展狀況進行測度與評估，但由于不同糧食作物生長和適應(yīng)環(huán)境存在差異，對特定糧食產(chǎn)物的生態(tài)效率研究較為缺乏。玉米在我國各地的種植范圍廣泛，不同種植區(qū)資源差異稟賦。因此研究玉米不同產(chǎn)區(qū)的生態(tài)效率，對推動各玉米種植區(qū)的低碳綠色發(fā)展有一定的引導(dǎo)作用。

基于以往對玉米不同生態(tài)類型區(qū)[5]的研究，筆者將我國玉米產(chǎn)區(qū)分為北部春播區(qū)、中部黃淮海平原種區(qū)和西部種區(qū)進行區(qū)域差異的研究，以玉米生產(chǎn)中碳排放總量作為限制因素，使用SBM模型方向性距離函數(shù)中的Malmquist-luenberger生產(chǎn)率指數(shù)（ML指數(shù)）和傳統(tǒng)以產(chǎn)量為導(dǎo)向的DEA-malmquist指數(shù)（M指數(shù)），分別測度我國各主要種區(qū)2004—2018年的生態(tài)效率和生產(chǎn)效率，對比反映玉米產(chǎn)業(yè)的現(xiàn)狀和差異;通過面板Tobit模型，分析不同產(chǎn)區(qū)的玉米生態(tài)效率和傳統(tǒng)M指數(shù)的影響因素，以實證結(jié)果更清晰地提出針對性建議，從而為實現(xiàn)不同產(chǎn)區(qū)玉米高質(zhì)量發(fā)展與資源環(huán)境可持續(xù)的目標提供參考。

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

1.1 研究方法

依據(jù)Chung等[6]的研究，結(jié)合SBM模型方向性距離函數(shù)構(gòu)建Malmquist-luenberger指數(shù)（ML指數(shù)），同時對玉米生產(chǎn)中的期望產(chǎn)出和非期望產(chǎn)出進行測度分析。從T期到T+1期包含非期望因素的Malmquist-luenberger指數(shù)函數(shù)，表述如下：

MLT+1 T=ST C（xT+1，aT+1，bT+1）

ST C（xT，aT，bT）×

ST+1 C（xT+1，aT+1，bT+1）

ST+1 C（xT，aT，bT）12

=ST+1 C（xT+1，aT+1，bT+1）

ST C（xT，aT，bT）

×

ST C（xT+1，aT+1，bT+1）

ST+1 C（xT+1，aT+1，bT+1）×

ST C（xT，aT，bT）

ST C（xT，aT，bT）12

=MLEFFCH（xT+1，aT+1，bT+1;x，a，b）×

MLTECH（xT+1，aT+1，bT+1;x，a，b）

（1）

式中，MLT+1 T表示從T期到T+1期的玉米生產(chǎn)生態(tài)效率變化情況。當MLT+1 T>1時，表示玉米生態(tài)效率從T到T+1期有所提高，反之則降低;MLT+1 T可以分解為技術(shù)效率變化指數(shù)（MLEFFCHT+1 T）和技術(shù)進步指數(shù)（MLTECHT+1 T）。

1.2 指標選取

參照以往對糧食產(chǎn)業(yè)的碳排放[7]研究，化肥比例在我國玉米生產(chǎn)六大碳源中的增長比例最大，與氮肥施用排放量一起所占總量的59.46%，其次是柴油占12.67%，農(nóng)藥占7.42%，種子占6.90%，這4項投入因素占碳排放總量的86.45%，可以很大程度地反映玉米生產(chǎn)中的非期望產(chǎn)出。碳排放系數(shù)參照美國橡樹嶺國家實驗室以及IPCC[8-11]公布的化肥0.895 6 kg/kg，農(nóng)藥0.592 7 kg/kg，柴油4.934 1 kg/kg，種子3.85 kg/kg來計算碳排放總量。選定的投入產(chǎn)出具體指標見表1。

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

根據(jù)我國玉米種植區(qū)分區(qū)[13]情況，選擇北部春播區(qū)包括河北、山西、內(nèi)蒙古、遼寧、吉林、黑龍江，中部黃淮海平原種區(qū)包括江蘇、安徽、山東、河南、湖北，西部種區(qū)包括重慶、四川、貴州、云南、陜西、甘肅、新疆。以上種植區(qū)2004—2018年玉米投入產(chǎn)出的數(shù)據(jù)均來源于《中國統(tǒng)計年鑒》《中國農(nóng)村統(tǒng)計年鑒》《全國農(nóng)產(chǎn)品成本收益資料匯編》《中國物價年鑒》（部分省份年度價格數(shù)據(jù)缺失以當年全國平均價格替代）。

2 結(jié)果與分析

2.1 玉米種植區(qū)生態(tài)效率測度結(jié)果與分析

根據(jù)獲取的投入產(chǎn)出數(shù)據(jù)，利用公式并借助MAXDEA軟件測算出我國玉米主要種植區(qū)的生態(tài)效率值。為突出種植過程中非期望因素的影響，使用DEAP2.1軟件測算出未考慮環(huán)境要素，以產(chǎn)量為導(dǎo)向的傳統(tǒng)Malmquist指數(shù)，并與生態(tài)效率的結(jié)果進行對比，結(jié)果見表2。

由表2可知，從分區(qū)情況來看，北部和中部種區(qū)的生態(tài)效率演變較為相近，增速相對較明顯;西部種區(qū)的增長相對較緩。不同種區(qū)的生態(tài)效率出現(xiàn)緩增的情況，可能與國家對“鐮刀彎地區(qū)”玉米生產(chǎn)適時調(diào)減種植面積，而相對減少了環(huán)境污染的負面效益，以及推進玉米市場結(jié)構(gòu)轉(zhuǎn)變的政策有關(guān)。進一步對生態(tài)效率進行分解為技術(shù)效率（MLEFFCH）和技術(shù)進步指數(shù)（MLTECH），以揭示驅(qū)動玉米種植區(qū)生態(tài)效率變化的因素。結(jié)果顯示，北部種區(qū)和中部平原種區(qū)的生態(tài)效率增長主要依靠的是前沿技術(shù)的進步，但北部種區(qū)的技術(shù)效率變化不顯著，而中部種區(qū)的技術(shù)效率呈退步趨勢，體現(xiàn)出的是大機械化、粗放式經(jīng)營的增長方式;西部種區(qū)的增長依靠的是技術(shù)效率提高，而前沿技術(shù)指數(shù)呈現(xiàn)衰退的狀態(tài)，表明西部種區(qū)的機械化投入程度不足。同時，參照鄭麗楠等[14]對我國農(nóng)業(yè)生態(tài)效率空間劃分，將玉米種植區(qū)按生態(tài)效率值劃分為相對高效區(qū)（ML指數(shù)大于等于1.03）、相對中效區(qū)（ML指數(shù)大于等于1且小于1.03）、相對低效區(qū)（ML指數(shù)小于1）。從分省市區(qū)的情況來看，根據(jù)各省市生態(tài)效率變化趨勢分析，內(nèi)蒙古、遼寧、河南、云南相對高效;河北、吉林、黑龍江、安徽、山東、重慶、四川、貴州、陜西、甘肅、新疆相對中效;山西、江蘇、湖北、新疆相對低效。

根據(jù)表2中玉米種植區(qū)的生態(tài)效率值與產(chǎn)量導(dǎo)向的Malmquist指數(shù)值，對比分析發(fā)現(xiàn)，未考慮環(huán)境要素的M指數(shù)增長均值為4.23%，加入非期望產(chǎn)出后ML指數(shù)年均增長為1.68%，增長效率明顯降低2.55%，因此在不考慮環(huán)境要素下，對玉米產(chǎn)業(yè)生產(chǎn)效率的測度結(jié)果會帶來一定的偏差;參考Fre等[15]和相關(guān)學(xué)者的研究，證明當投入一定時，若ML指數(shù)>M指數(shù)，則意味著期望產(chǎn)出的增長速度高于非期望產(chǎn)出的增長，那么該生產(chǎn)單位相對低碳，反之則相對高碳不環(huán)保。各玉米種區(qū)的數(shù)據(jù)對比可知，北部春播區(qū)（0.2%）的生產(chǎn)方式相對低碳，而中部平原種區(qū)（-2.39%）、西部山地種區(qū)（-5.03%）的生產(chǎn)方式相對高碳，指數(shù)結(jié)果證明玉米各種植區(qū)的生產(chǎn)方式差異明顯。

2.2 玉米種植區(qū)生態(tài)效率的影響因素分析

2.2.1 指標選取。

為進一步探討對生態(tài)效率的影響因素，選取宏觀環(huán)境、農(nóng)業(yè)發(fā)展和政策支持3個方面進行考量，指標選取見表3。

2.2.2 模型構(gòu)建。

構(gòu)建以生態(tài)效率值為因變量，各影響因素為自變量的面板模型;同時構(gòu)建以產(chǎn)量為導(dǎo)向的傳統(tǒng)Malmquist指數(shù)為因變量，相同影響因素為自變量的面板模型，對比結(jié)果是否具有差異性。為消除數(shù)據(jù)的不穩(wěn)定性，所有變量均取對數(shù)形式。

lnCEE it=β 0+β 1lngdp i，t+β 2lnarea i，t+β 3lnurb i，t+β 4lngovs i，t+β 5lnmach i，t+β 6lndis i，t+β 7lnECI i，t+β 8lnEPM i，t+ε i，t（2）

式中，CEE i，t為選取的各種區(qū)玉米的生態(tài)效率;β 0為常數(shù)項，β為回歸系數(shù)，i=1，2，…，18代表玉米產(chǎn)出省份，t代表時序，ε i，t代表隨機擾動項。

2.2.3 回歸過程與結(jié)果分析。

因變量取值均為效率值具有截斷特征，使用傳統(tǒng)OLS回歸模型可能存在結(jié)果的偏差。故選擇Tobit面板模型進行回歸，使用Stata 16.0軟件運算結(jié)果見表4。

（1）北部種區(qū)是我國玉米種植的重要地區(qū)，其中東北三省的玉米產(chǎn)量占全國玉米總產(chǎn)量的33.1%[16]。通過對北部種區(qū)生態(tài)效率的影響因素分析，生態(tài)治理水平有顯著的正向影響，而機械密度、受災(zāi)率對生態(tài)效率有顯著的負向影響;對比M指數(shù)的回歸結(jié)果，城鎮(zhèn)化水平、機械密度、生態(tài)治理水平對生產(chǎn)效率有顯著的正向影響，產(chǎn)業(yè)結(jié)構(gòu)對生產(chǎn)效率有顯著的負向影響。從結(jié)果對比來看，北部玉米種區(qū)的差異性主要體現(xiàn)在機械密度，機械密度投入的加大，有利于提升玉米的生產(chǎn)效率，因此喬丹等[5]建議加大玉米產(chǎn)業(yè)機械化投入和使用。但從糧食產(chǎn)業(yè)高質(zhì)量發(fā)展角度，認為過度使用機械會增大環(huán)境的負擔，使用不當和效率不高也是對資源的消耗和浪費，不利于北部玉米種區(qū)生態(tài)效率的提升，因此合理布局機械化生產(chǎn)和指導(dǎo)顯得尤為重要。

（2）中部玉米種區(qū)是我國第二大產(chǎn)區(qū)，通過對中部種區(qū)的生態(tài)效率影響因素分析，結(jié)果顯示環(huán)保投入對中部玉米種區(qū)的生態(tài)效率有顯著的正向影響，而負向影響較為明顯的是受災(zāi)率。因中部種區(qū)為黃淮海平原地區(qū)，降雨分布不均勻，容易受到春旱和夏澇等自然災(zāi)害的影響[13]，因此玉米產(chǎn)業(yè)的生態(tài)效率受到自然災(zāi)害的負面影響較大;對比M指數(shù)的結(jié)果，機械密度對M指數(shù)有顯著的正向影響，產(chǎn)業(yè)結(jié)構(gòu)對M指數(shù)有顯著的負面影響;從結(jié)果對比來看，中部玉米種區(qū)的生產(chǎn)效率影響因素分析結(jié)果與北部種區(qū)相似，農(nóng)業(yè)技術(shù)的利用較為充分，土地利用率高，種植活動對傳統(tǒng)農(nóng)資的依賴大，不利于中部種區(qū)生態(tài)效率的提升。因此，加大對中部地區(qū)的環(huán)保投入力度，對促進中部玉米生態(tài)效率有重要的推動作用。

（3）西部玉米種區(qū)多為山地地形，土地較為貧瘠，經(jīng)濟水平落后。該種區(qū)生態(tài)效率的影響因素結(jié)果顯示，環(huán)保投入的正向影響顯著，而城鎮(zhèn)化水平、財政投入、受災(zāi)率對西部種區(qū)生態(tài)效率的負向影響顯著。西部種區(qū)包含我國較多經(jīng)濟不發(fā)達地區(qū)，城鎮(zhèn)化的發(fā)展雖然會加快區(qū)域經(jīng)濟的發(fā)展，但也加速農(nóng)村的勞動力、土地和資金徑直流入城鎮(zhèn)地區(qū)，導(dǎo)致西部農(nóng)村的資源配置失衡，不利于西部玉米種區(qū)的生態(tài)效率提升;另一方面，財政對農(nóng)業(yè)的投入傾斜增大，對西部地區(qū)農(nóng)業(yè)發(fā)展有重要的推動作用，但同時會刺激糧農(nóng)的種植積極性，可能會存在為提高產(chǎn)量而過度使用農(nóng)資，增加環(huán)境污染的風險性，從而降低西部玉米產(chǎn)業(yè)的生態(tài)效率。對比M指數(shù)的結(jié)果發(fā)現(xiàn)，西部種區(qū)的受災(zāi)率增加會同時顯著降低西部玉米種區(qū)的生態(tài)效率和生產(chǎn)效率，農(nóng)業(yè)自然災(zāi)害直接影響的是作物產(chǎn)量，削弱農(nóng)業(yè)生產(chǎn)者的信心;環(huán)保投入的增加會顯著提升西部種區(qū)的生態(tài)、生產(chǎn)效率，基于十四五規(guī)劃要積極推進農(nóng)業(yè)生產(chǎn)更加低碳，因環(huán)保治理投資力度的增大，有利于西部玉米產(chǎn)業(yè)生態(tài)效率的有效增長。

3 結(jié)論與建議

3.1 結(jié)論

利用SBM模型中的ML指數(shù)法，基于面板數(shù)據(jù)測算了我國玉米三大種植區(qū)的生態(tài)效率，并同時測算了以產(chǎn)量為導(dǎo)向的傳統(tǒng)Malmquist值，將兩類指數(shù)的結(jié)果進行對比;運用Tobit面板模型對兩類指數(shù)進行相同影響因素分析。

（1）加入非期望因素測算的生態(tài)效率更能真實地反映各種植區(qū)的生產(chǎn)現(xiàn)狀，且呈現(xiàn)明顯的區(qū)域差異。我國玉米產(chǎn)業(yè)整體還是處于粗放式經(jīng)營的發(fā)展模式，其中生態(tài)效率高效區(qū)數(shù)量較少且分布不均，但有較大的改善空間。北部種區(qū)整體生態(tài)效率較高，增長明顯且生產(chǎn)方式較為低碳，其生態(tài)效率的增長主要依靠于技術(shù)進步;中部種區(qū)的平均生態(tài)效率較高，但技術(shù)效率退步明顯，生產(chǎn)方式為高碳的傳統(tǒng)方式;西部種區(qū)的生態(tài)效率增長緩慢，技術(shù)投入不足，且生產(chǎn)方式是犧牲環(huán)境的高碳方式，但整體提升空間較大。

（2）雖然在受災(zāi)率、生態(tài)治理水平和環(huán)保投入方面對玉米種植區(qū)的生態(tài)效率和生產(chǎn)效率存在趨同的影響作用，但產(chǎn)業(yè)機構(gòu)、城鎮(zhèn)化水平、機械密度和財政投入等因素對不同玉米種植區(qū)的兩類指數(shù)影響存在明顯差異。北部種區(qū)機械密度的投入增加雖然會促進玉米產(chǎn)量導(dǎo)向的生產(chǎn)效率提升，但同時也會降低玉米產(chǎn)業(yè)的生態(tài)效率;中西部種區(qū)對產(chǎn)業(yè)結(jié)構(gòu)的調(diào)整，加大投入財政力度可能會促進種區(qū)的生產(chǎn)效率和產(chǎn)量，但不利于該種區(qū)生態(tài)效率的增長。

3.2 政策建議

由于玉米種區(qū)的個體異質(zhì)性，我國玉米各種植區(qū)應(yīng)根據(jù)自身發(fā)展現(xiàn)狀，并結(jié)合自身優(yōu)勢，根據(jù)我國糧食產(chǎn)業(yè)高質(zhì)量發(fā)展方向，因地制宜地提出發(fā)展策略來提高我國各地玉米產(chǎn)業(yè)的生態(tài)效率。

（1）北部種區(qū)作為玉米重要種植區(qū)，繼續(xù)保持低碳生產(chǎn)方式的同時，要合理布局機械投入結(jié)構(gòu)，配合玉米規(guī)?；N植模式，進一步提高機械化使用效率，提高玉米種植質(zhì)量和產(chǎn)品質(zhì)量，利用種植資源的優(yōu)勢積極研發(fā)玉米新品種，并形成北部種區(qū)優(yōu)質(zhì)玉米的品牌化效應(yīng)。另外，北部種區(qū)要增強生態(tài)治理的能力，提高土地的利用效率，對農(nóng)業(yè)生產(chǎn)建立災(zāi)害預(yù)警機制，提高生產(chǎn)活動中災(zāi)害的防范能力和管理水平。

（2）中部種區(qū)要積極轉(zhuǎn)變生產(chǎn)方式，加快綠色新技術(shù)的引入，提高技術(shù)投入的經(jīng)濟效率和力度，從而轉(zhuǎn)變傳統(tǒng)的大機械化模式。優(yōu)化農(nóng)業(yè)生產(chǎn)資料的布局，提高資源利用效率，降低使用不合理所帶來的資源浪費和環(huán)境污染的可能性。另外中部種區(qū)要加強環(huán)保投入力度，加強玉米生產(chǎn)對環(huán)保技術(shù)的認知能力，利用各種形式的宣傳環(huán)保教育活動，提高種植戶的環(huán)保意識。

（3）西部種區(qū)要加大機械化投入力度，增加現(xiàn)代機械的應(yīng)用和研究，推進以新型農(nóng)資替代傳統(tǒng)碳源物資并提高資源的利用率，推動農(nóng)業(yè)綠色技術(shù)和綠色肥料的使用，提高耕地質(zhì)量，積極轉(zhuǎn)變傳統(tǒng)的高碳生產(chǎn)方式。另外對于西部種區(qū)，城鎮(zhèn)化發(fā)展進程中要平衡農(nóng)業(yè)資源的合理分配，注重財政補貼結(jié)構(gòu)的優(yōu)化與調(diào)整，加大對節(jié)能環(huán)保低耗農(nóng)資研發(fā)傾斜，同時加強政策的驅(qū)動作用。

參考文獻

[1] 蔡榮，汪紫鈺.中國玉米生產(chǎn)的環(huán)境效率及其時空分異：兼論玉米臨時收儲政策的環(huán)境影響[J].農(nóng)林經(jīng)濟管理學(xué)報，2019，18（5）：569-578.

[2] 崔愛民，張久剛，張虎，等.我國玉米生產(chǎn)現(xiàn)狀及發(fā)展變革[J].中國農(nóng)業(yè)科技導(dǎo)報，2020，22（7）：10-19.

[3] 于婷，郝信波.糧食主產(chǎn)區(qū)農(nóng)業(yè)生態(tài)效率時空特征及改善路徑研究[J].生態(tài)經(jīng)濟，2018，34（9）：104-110.

[4] 魯慶堯，王樹進，孟祥海.基于SBM模型的我國糧食生產(chǎn)生態(tài)效率測度與PS收斂檢驗[J].農(nóng)村經(jīng)濟，2020（12）：24-32.

[5] 喬丹，陸遷.不同生態(tài)類型區(qū)玉米生產(chǎn)技術(shù)效率及有偏演進模式[J].華南農(nóng)業(yè)大學(xué)學(xué)報（社會科學(xué)版），2016，15（5）：28-36.

[6] CHUNG Y H，F(xiàn)RE R，GROSSKOPF S.Productivity and undesirable outputs：A directional distance function approach[J].Journal of environmental management，1997，51（3）：229-240.

[7] 王海娜.我國玉米生產(chǎn)碳排放效率研究[D].長春：吉林大學(xué)，2018：44-45.

[8] WEST T O，MARLAND G.A synthesis of carbon sequestration，carbon emissions，and net carbon flux in agriculture：Comparing tillage practices in the United States[J].Agriculture ecosystem and environment，2002，91（1/2/3）：217-232.

[9] 李波，張俊飚，李海鵬.中國農(nóng)業(yè)碳排放時空特征及影響因素分解[J].中國人口·資源與環(huán)境，2011，21（8）：80-86.

[10] 田云，張俊飚，李波.中國農(nóng)業(yè)碳排放研究：測算、時空比較及脫鉤效應(yīng)[J].資源科學(xué)，2012，34（11）：2097-2105.

[11] WEST T O，POST W M.Soil organic carbon sequestration rates by tillage and crop rotation：A global data analysis[J].Soil science society of America journal，2002，66（6）：1930-1946.

[12] 周志花.利用LCA法核算農(nóng)作物生產(chǎn)碳足跡[D].北京：中國農(nóng)業(yè)科學(xué)院，2018.

[13] 王淑萍.我國玉米種植產(chǎn)區(qū)的劃分[J].養(yǎng)殖技術(shù)顧問，2010（6）：74.

[14] 鄭麗楠，洪名勇.中國農(nóng)業(yè)生態(tài)效率的時空特征及驅(qū)動因素[J].江西財經(jīng)大學(xué)學(xué)報，2019（5）：46-56.

[15] FRE R，GROSSKOPF S，PASURKA C A JR.Accounting for air pollution emissions in measures of state manufacturing productivity growth[J].Journal of regional science，2001，41（3）：381-409.

[16] 李明文，王振華，張廣勝.東北玉米種植結(jié)構(gòu)調(diào)整與糧食高質(zhì)量增長：基于全要素生產(chǎn)率視角[J].農(nóng)業(yè)現(xiàn)代化研究，2019，40（5）：745-754.