摘 要 為了厘清中國(guó)城市耕地?zé)o機(jī)肥的氮素投入狀況和經(jīng)濟(jì)產(chǎn)出效益,采用區(qū)域種植業(yè)經(jīng)濟(jì)產(chǎn)出與區(qū)域無(wú)機(jī)肥氮素總投入量的比值作為無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率指示因子,研究了2000—2019年間全國(guó)、各城市、九大農(nóng)業(yè)區(qū)、東中西部、糧食產(chǎn)區(qū)的無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率?;诤嗣芏惹€、空間自相關(guān)分析等方法分析耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率的時(shí)空差異特征。結(jié)果表明:1)總體來(lái)看各城市和區(qū)域耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率存在長(zhǎng)期向好傾向;2)從區(qū)域尺度上看,東中西地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)“東部高西部低”分布格局;3)從空間分異特征來(lái)看,在全局范圍內(nèi)呈顯著的空間正相關(guān)性,且存在集聚分布態(tài)勢(shì)和空間依賴性,空間集聚程度表現(xiàn)出極化效應(yīng)有所減緩。
關(guān)鍵詞 無(wú)機(jī)肥;氮素經(jīng)濟(jì)產(chǎn)出效率;動(dòng)態(tài)演變;時(shí)空分異
中圖分類號(hào):S963.9 文獻(xiàn)標(biāo)志碼:A DOI:10.19415/j.cnki.1673-890x.2022.11.031
氮元素是維持生命的關(guān)鍵元素。改革開(kāi)放以來(lái),中國(guó)農(nóng)業(yè)成功利用占世界7%的耕地養(yǎng)活占世界22%的人口,取得了舉世矚目的成就,無(wú)機(jī)肥的高投入在其中發(fā)揮了重要的作用[1]。但無(wú)機(jī)肥的過(guò)量施用或不當(dāng)施用會(huì)引起耕地土壤酸化、耕層變淺[2],造成局部空氣污染、產(chǎn)生溫室效應(yīng),造成水體富營(yíng)養(yǎng)化和地下水污染等農(nóng)業(yè)面源污染問(wèn)題[3]。農(nóng)業(yè)發(fā)展對(duì)生態(tài)環(huán)境的脅迫效應(yīng)十分劇烈[4],進(jìn)一步加強(qiáng)耕地?zé)o機(jī)肥氮素管理對(duì)提高區(qū)域種植業(yè)生產(chǎn)率和改善耕地質(zhì)量具有積極的影響。耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率反映了種植業(yè)無(wú)機(jī)肥氮素投入的經(jīng)濟(jì)效益,探討中國(guó)城市耕地?zé)o機(jī)肥的氮素投入狀況,測(cè)算耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率,有助于了解各城市耕地利用發(fā)展的高質(zhì)量成效,為中國(guó)城市合理施用無(wú)機(jī)肥和有效控制無(wú)機(jī)肥面源污染提供科學(xué)依據(jù)。
1? 材料與方法
1.1? 研究區(qū)域概況
為科學(xué)反映各地級(jí)及以上城市在全國(guó)范圍內(nèi)的空間分異特征,本文將直轄市納入研究樣本,涉及我國(guó)333個(gè)地級(jí)行政區(qū)(包括293個(gè)地級(jí)市、7個(gè)地區(qū)、30個(gè)自治州、3個(gè)盟),以及4個(gè)直轄市(北京、上海、重慶、天津)。
研究區(qū)域按3種類別進(jìn)行劃分,1)按照地理分區(qū),以省級(jí)行政區(qū)為標(biāo)準(zhǔn)分為東中西地區(qū);2)按照糧食生產(chǎn)情況劃分,分為糧食主產(chǎn)區(qū)和非糧食主產(chǎn)區(qū);3)根據(jù)《中國(guó)綜合農(nóng)業(yè)區(qū)劃》,將全國(guó)劃分為9個(gè)農(nóng)業(yè)區(qū)。
1.2? 數(shù)據(jù)來(lái)源
本研究所需數(shù)據(jù)包括全國(guó)各城市的氮肥施用折純量、復(fù)合肥施用折純量、耕地面積、種植業(yè)總產(chǎn)值,時(shí)間尺度為2000—2019年,數(shù)據(jù)來(lái)源于我國(guó)的統(tǒng)計(jì)年鑒、農(nóng)村統(tǒng)計(jì)年鑒、農(nóng)業(yè)統(tǒng)計(jì)年鑒、城市統(tǒng)計(jì)年鑒、城市建設(shè)統(tǒng)計(jì)年鑒以及各城市(州)統(tǒng)計(jì)年鑒、各城市(州)國(guó)民經(jīng)濟(jì)與社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)。
1.3? 測(cè)算方法
1.3.1? 無(wú)機(jī)肥氮素投入量測(cè)算方法
我國(guó)農(nóng)業(yè)生產(chǎn)上含氮元素的無(wú)機(jī)肥肥料主要包括氮肥和復(fù)合肥。氮肥施用折純量和復(fù)合肥施用折純的含氮量構(gòu)成無(wú)機(jī)肥氮元素的輸入源,本研究中根據(jù)我國(guó)常用的不同配方復(fù)合肥中氮含量折算,折純的復(fù)合肥中氮素含量為32.2%[5]。
[Nfer=NF+CF×rN]? ?(1)
式(1)中,[NF]和[CF]分別表示氮肥和復(fù)合肥的施用折純量;[rN]表示復(fù)合肥中氮素的含量,無(wú)機(jī)肥氮素投入涉及的數(shù)據(jù)主要為地市級(jí)行政單位氮、磷、鉀、復(fù)合肥施用量及無(wú)機(jī)肥施用總量(折純量),耕地面積等,數(shù)據(jù)來(lái)源為我國(guó)歷年地級(jí)行政區(qū)及省份的統(tǒng)計(jì)年鑒。部分缺失數(shù)據(jù)用省級(jí)數(shù)據(jù)按耕地面積插值替代[6]。
1.3.2? 無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率測(cè)算方法
本文參考?xì)W盟國(guó)際上提出的氮素利用率(NUE),確定采用區(qū)域種植業(yè)經(jīng)濟(jì)產(chǎn)出與區(qū)域無(wú)機(jī)肥氮素總投入量的比值作為無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率(Nitrogen economic output efficiency)指示因子來(lái)評(píng)估研究中國(guó)城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率大小。
[NEOE=AOVNinput×100%]? (2)
式(2)中,NEOE表示區(qū)域無(wú)機(jī)肥區(qū)域氮素經(jīng)濟(jì)產(chǎn)出效率,AOV表示區(qū)域種植業(yè)的經(jīng)濟(jì)產(chǎn)出(億元),數(shù)據(jù)來(lái)源于各地市級(jí)行政單元統(tǒng)計(jì)年鑒中農(nóng)林牧漁業(yè)總產(chǎn)值分類中的農(nóng)業(yè)總產(chǎn)值;[Ninput]表示農(nóng)作物種植過(guò)程中無(wú)機(jī)肥氮素總投入量(萬(wàn)t)。
2? 結(jié)果與分析
2.1? 無(wú)機(jī)肥氮素投入量時(shí)間演變特征
中國(guó)耕地?zé)o機(jī)肥氮素投入量從2000年的2 662萬(wàn)t增長(zhǎng)到2019年的2 906萬(wàn)t,年均上升率為0.46%,氮肥施用折純量氮投入從2000年的2 358萬(wàn)t下降到2019年的2 171萬(wàn)t,年均上升率為-0.4%,在全國(guó)城市耕地?zé)o機(jī)肥氮素投入量中貢獻(xiàn)率從2000年88.59%的降低到2019年的74.71%;復(fù)合肥施用折純量從2000年的862萬(wàn)t上升到2019年的2 128萬(wàn)t,復(fù)合肥施用折純量氮素投入量從2000年的304萬(wàn)t增長(zhǎng)到2019年的735萬(wàn)t,增長(zhǎng)率為142%,年均上升率為7.09%,在全國(guó)城市耕地?zé)o機(jī)肥氮素投入量中貢獻(xiàn)率從2000年11.41%的增長(zhǎng)到2019年的25.29%。
2.2? 無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率時(shí)間演變特征
2.2.1? 中國(guó)城市無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率時(shí)間趨勢(shì)
研究期間中國(guó)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率總體上呈現(xiàn)持續(xù)上漲趨勢(shì),均值為12.59,最高值為2019年的21.79,最低值為2000年的5.43,從2000年的5.43增長(zhǎng)為2019年的21.79,年均增長(zhǎng)率為15.05%。中國(guó)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率的變動(dòng)主要分為3個(gè)階段。第1階段為2000—2003年,該階段出效率分布在[5,6],這段時(shí)期的經(jīng)濟(jì)產(chǎn)出水平呈現(xiàn)緩慢上升趨勢(shì)。相對(duì)整個(gè)時(shí)期來(lái)說(shuō)較低;第2階段為2006—2015年,該階段效率分布在[6,18],這段時(shí)期的氮素產(chǎn)出水平呈現(xiàn)勻速上升趨勢(shì),相對(duì)整個(gè)研究時(shí)期來(lái)說(shuō)屬于中效水平;第3階段為2017—2019年,該階段效率分布在[18,22],這段時(shí)期增長(zhǎng)幅度較大,可能與政府重視糧食安全,推進(jìn)農(nóng)業(yè)可持續(xù)和綠色高質(zhì)量發(fā)展有關(guān)。隨著經(jīng)濟(jì)的快速發(fā)展,農(nóng)業(yè)科技進(jìn)步、農(nóng)業(yè)機(jī)械化水平得到提升,2017—2019年的經(jīng)濟(jì)效益較高。從總體上來(lái)看,我國(guó)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率總體上呈現(xiàn)持續(xù)上升趨勢(shì),亟需繼續(xù)提高無(wú)機(jī)肥的氮素產(chǎn)出效率,以實(shí)現(xiàn)經(jīng)濟(jì)效益和生態(tài)環(huán)境保護(hù)的可持續(xù)協(xié)調(diào)發(fā)展。
為了進(jìn)一步研究我國(guó)城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率隨時(shí)間變化趨勢(shì),本文根據(jù)2000、2005、2010、2015、2019年城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率核密度曲線展現(xiàn)出其動(dòng)態(tài)變化。從峰值分布可以看出,我國(guó)城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)明顯的單峰分布態(tài)勢(shì),說(shuō)明2000—2019年相較于處于低效的城市數(shù)量而言,較高效的數(shù)量有明顯的增加;從時(shí)間維度上看,5個(gè)年份的低效率峰值正逐年向右移動(dòng),說(shuō)明在研究時(shí)期內(nèi),低效率城市的效率值得到大幅度的提升。就波峰陡峭程度而言,也可以看出低效率峰值逐年降低,高效率峰值則逐年升高??傮w上也說(shuō)明在研究時(shí)期內(nèi),我國(guó)城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率兩極分化現(xiàn)象正不斷減弱,其中低效城市效率提高的速度較快,當(dāng)?shù)卣孕杓訌?qiáng)區(qū)域間的協(xié)同合作和協(xié)調(diào)發(fā)展。
2.2.2? 不同區(qū)域耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率時(shí)間變化
從東中西地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率來(lái)看,2000—2019年三大地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率總體上呈現(xiàn)上升的時(shí)間變化趨勢(shì)。考察期內(nèi),東中西地區(qū)效率均值分別為13.64、12.08、9.64,總體上呈現(xiàn)“東部高西部低”分布格局。東部地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率從2000年的6.6增長(zhǎng)為2019年的23.89,除2004年有所下降外,整體上逐年上升,年均增長(zhǎng)率為13.09%。中部地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率從2000年的4.21增長(zhǎng)為2019年的18.04,年均增長(zhǎng)率為16.45%;西部地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率從2000年的5.75增長(zhǎng)到2019年的24.78,年均增長(zhǎng)率為16.55%。2000—2016年,中西部地區(qū)變動(dòng)趨勢(shì)相對(duì)一致,但是從2017年開(kāi)始逐漸拉開(kāi)差距,西部地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率上升幅度較大,但中部增長(zhǎng)幅度較小。
進(jìn)一步分析糧食主產(chǎn)區(qū)和非糧食主產(chǎn)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率,其變化趨勢(shì)與全國(guó)及東中西三大區(qū)域變化趨勢(shì)一致,呈現(xiàn)“波動(dòng)式”上升演進(jìn)趨勢(shì)??疾炱趦?nèi)糧食主產(chǎn)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率平均值大于非糧食主產(chǎn)區(qū),均值分別為11.54、11.71。糧食主產(chǎn)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率從2000年的5.32增長(zhǎng)為2019年的20.64,年均增長(zhǎng)率為14.39%;非糧食主產(chǎn)區(qū)的耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率從2000年的5.66增長(zhǎng)為2019年的23.99,年均增長(zhǎng)率為16.19%。
2.2.3? 分農(nóng)業(yè)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率時(shí)間變化
從分農(nóng)業(yè)區(qū)層面看,中國(guó)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率差異較為明顯。總體上2000—2019年中國(guó)各農(nóng)業(yè)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)波動(dòng)性上升趨勢(shì)。青藏高原區(qū)2019年效率值達(dá)到70.65且2000—2019年耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率平均值位列第一,表明青藏高原區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率處于生產(chǎn)前沿面,用相對(duì)較少的無(wú)機(jī)肥氮素投入量獲得經(jīng)濟(jì)收益。華南區(qū)2019年耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率值達(dá)到31.85且2000—2019年效率平均值位列第二,表明華南用相對(duì)較少的無(wú)機(jī)肥氮素投入量獲得較高的經(jīng)濟(jì)收益。耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率平均值小于10的省份有2個(gè),依次為北方干旱半干旱區(qū)(9.48)、黃土高原區(qū)(6.49),表明相對(duì)來(lái)講相同的無(wú)機(jī)肥氮素投入量獲得較少的經(jīng)濟(jì)收益,可能與區(qū)域土壤、氣候、用水條件、種植業(yè)發(fā)展方式有一定的關(guān)系。2000—2019年四川盆地、青藏高原區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率增長(zhǎng)較快,年均增長(zhǎng)率均達(dá)到20%以上;其他農(nóng)業(yè)區(qū)年均增長(zhǎng)率均達(dá)到10%以上。總體上九大農(nóng)業(yè)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)波動(dòng)性上升趨勢(shì),亟需提高耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率,大力實(shí)施糧食產(chǎn)能提升工程,集中支持適宜區(qū)域、重點(diǎn)品種,因地制宜發(fā)展現(xiàn)代化農(nóng)業(yè)。
2.3? 無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率空間變化特征
2.3.1? 全局空間自相關(guān)分析
本文采用Geoda軟件對(duì)Moran’s I散點(diǎn)圖、LISA聚集圖等進(jìn)行分析,從而發(fā)現(xiàn)城市無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率在空間分布上存在的空間聚集特征。
全局空間正相關(guān)顯著。2000—2019年各城市耕地氮素經(jīng)濟(jì)產(chǎn)出效率的 Moran’s I指數(shù)全部為正,而且Z值均大于0,P值均為0.001,基本上都通過(guò)了1%水平下顯著性檢驗(yàn),在全局范圍內(nèi)呈顯著的空間正相關(guān)性。地區(qū)間不平衡性尚存,高值集聚區(qū)增多。結(jié)合2000—2019年各城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率指數(shù)散點(diǎn)圖,發(fā)現(xiàn)其分布由集中趨于分散,這也在一定程度上表明各城市無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率水平極化現(xiàn)狀嚴(yán)重,有向均衡化發(fā)展態(tài)勢(shì)靠近的趨勢(shì)。
2.3.2? 局部空間自相關(guān)分析
局部空間正相關(guān)不顯著。根據(jù)LISA顯著性圖可知,2000—2019年典型年份300左右以上城市的耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率不顯著,P值小于0.01的城市介于20~40個(gè),這些城市主要分布在海南省、黃土高原及青藏高原等地區(qū),在城市全局范圍內(nèi)呈顯著的空間正相關(guān)性,且存在集聚分布態(tài)勢(shì)和空間依賴性??臻g聚集類型不顯著。2000—2019年,270個(gè)以上的地級(jí)行政區(qū)耕地氮素經(jīng)濟(jì)產(chǎn)出效率聚集類型不顯著,2000—2019年H-H型的地級(jí)行政區(qū)數(shù)量10~15個(gè),主要集中在東北平原區(qū)北部地級(jí)市、青藏高原區(qū)南部地區(qū)。
3? 結(jié)論與討論
1)研究期間,全國(guó)各城市耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率存在長(zhǎng)期向好傾向,從2000年的5.43增長(zhǎng)為2019年的21.79,均值為12.59,年均增長(zhǎng)率為15.05%,總體上呈現(xiàn)持續(xù)上升趨勢(shì),低效率城市效率提高的速度較快,兩極分化現(xiàn)象正不斷減弱。亟需繼續(xù)提高耕地利用的氮素效率,在測(cè)土配方的基礎(chǔ)上,平衡綜合使用無(wú)機(jī)肥和有機(jī)肥;調(diào)整要素配置,合理控制無(wú)機(jī)肥消費(fèi)總量,確定生產(chǎn)資源配置的最佳方案;更好地推廣和教育,提高農(nóng)民對(duì)土壤肥力和無(wú)機(jī)肥施用的環(huán)保認(rèn)知,實(shí)現(xiàn)經(jīng)濟(jì)效益和生態(tài)環(huán)境保護(hù)的可持續(xù)協(xié)調(diào)發(fā)展。
2)2000—2019年?yáng)|中西地區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率總體上呈現(xiàn)上升的時(shí)間變化趨勢(shì)。糧食主產(chǎn)區(qū)、非糧食主產(chǎn)區(qū)和九大農(nóng)業(yè)區(qū)耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)“波動(dòng)式”上升演進(jìn)趨勢(shì)??疾炱趦?nèi),耕地?zé)o機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率呈現(xiàn)“東高中低”分布格局,糧食主產(chǎn)區(qū)效率均值比非糧食主產(chǎn)區(qū)高,青藏高原區(qū)和華南區(qū)的平均值較其他農(nóng)業(yè)區(qū)高。各區(qū)域應(yīng)根據(jù)當(dāng)?shù)刭Y源優(yōu)勢(shì)和農(nóng)業(yè)特色,因地制宜發(fā)展高產(chǎn)高效農(nóng)作物。通過(guò)建立和發(fā)展特色農(nóng)業(yè)基地,開(kāi)發(fā)優(yōu)勢(shì)農(nóng)產(chǎn)品,積極推進(jìn)特色農(nóng)業(yè)產(chǎn)業(yè)化和規(guī)模經(jīng)營(yíng)發(fā)展。同時(shí),積極推進(jìn)農(nóng)業(yè)與二三產(chǎn)業(yè)融合發(fā)展,不斷延伸農(nóng)業(yè)產(chǎn)業(yè)鏈和價(jià)值鏈,促進(jìn)高技術(shù)產(chǎn)業(yè)發(fā)展。利用互聯(lián)網(wǎng)平臺(tái),帶動(dòng)特色農(nóng)業(yè)產(chǎn)業(yè)的發(fā)展。
3)從空間分異特征來(lái)看,整體上呈現(xiàn)出東部高西部低的特征。從全局自相關(guān)看,在全局范圍內(nèi)呈顯著的空間正相關(guān)性,且存在集聚分布態(tài)勢(shì)和空間依賴性,各城市無(wú)機(jī)肥氮素經(jīng)濟(jì)產(chǎn)出效率水平逐步向均衡化發(fā)展態(tài)勢(shì)靠近,高值集聚現(xiàn)象不斷增強(qiáng)。從局部自相關(guān)看,2000—2019年典型年份300左右以上城市的耕地氮素經(jīng)濟(jì)產(chǎn)出效率不顯著。東部地區(qū)和中部地區(qū)大部分城市的生產(chǎn)率聚集類型不顯著。亟需推進(jìn)種植業(yè)結(jié)構(gòu)優(yōu)化升級(jí)和綠色高質(zhì)量轉(zhuǎn)型,推廣耕地利用可再生的綠色清潔能源,加強(qiáng)區(qū)域間協(xié)同發(fā)展,建立跨區(qū)域無(wú)機(jī)肥減量與環(huán)境保護(hù)合作機(jī)制。
參考文獻(xiàn):
[1]? 朱兆良,孫波,楊林章,等.我國(guó)農(nóng)業(yè)面源污染的控制政策和措施[J].科技導(dǎo)報(bào),2005(4):47-51.
[2]? 張福鎖.我國(guó)農(nóng)田土壤酸化現(xiàn)狀及影響[J].民主與科學(xué),2016(6):26-27.
[3]? 李海鵬.中國(guó)農(nóng)業(yè)面源污染的經(jīng)濟(jì)分析與政策研究[D].武漢:華中農(nóng)業(yè)大學(xué),2007.
[4]? 李海鵬,羅麗,張雄,等.中國(guó)農(nóng)業(yè)能源效率動(dòng)態(tài)演變及其影響因素[J].中國(guó)人口·資源與環(huán)境,2020,30(12):105-115.
[5]? 高偉,郭懷成,后???中國(guó)大陸市域人類活動(dòng)凈氮輸入量(NANI)評(píng)估[J].北京大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,50(5):951-959.
[6]? 田沛佩,盧宏瑋,李丹,等.2008~2017年中國(guó)地級(jí)市化肥施用碳足跡的時(shí)空演變格局[J].中國(guó)環(huán)境科學(xué),2021,41(2):967-973.
收稿日期:2022-04-03
作者簡(jiǎn)介:韓莉(1998—),女,河南信陽(yáng)人,在讀碩士,主要研究方向?yàn)橘Y源經(jīng)濟(jì)與管理。E-mail: 1076458661@qq.com。