葉思菁，宋長(zhǎng)青，程 鋒，張蕾娜，程昌秀，張 超，楊建宇，朱德海

中國(guó)耕地健康產(chǎn)能綜合評(píng)價(jià)與試點(diǎn)評(píng)估研究

葉思菁1,2，宋長(zhǎng)青1,2，程 鋒3，張蕾娜3，程昌秀1,2※，張 超4，楊建宇4，朱德海4

（1. 北京師范大學(xué)地表過程與資源生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京，100875；2. 北京師范大學(xué)地理科學(xué)學(xué)部，北京，100875；3. 自然資源部國(guó)土整治中心，北京，100035；4. 自然資源部農(nóng)用地質(zhì)量與監(jiān)控重點(diǎn)實(shí)驗(yàn)室，北京，100083）

評(píng)價(jià)中國(guó)耕地健康產(chǎn)能狀況，分析危及耕地健康的區(qū)域性關(guān)鍵問題，對(duì)探索可持續(xù)的耕地利用路徑，維護(hù)國(guó)家糧食安全與生態(tài)安全具有重要意義。耕地健康產(chǎn)能是耕地自然系統(tǒng)與利用系統(tǒng)綜合作用的結(jié)果，直接計(jì)算耕地健康產(chǎn)能狀況是復(fù)雜困難的。因此，該文嘗試采用一種簡(jiǎn)化的方法評(píng)價(jià)耕地健康產(chǎn)能，即分別評(píng)價(jià)耕地健康與耕地產(chǎn)能狀況，并通過耦合耕地健康與產(chǎn)能，分析耕地利用合理程度與生產(chǎn)能力之間的相互作用特征。其中耕地健康狀況強(qiáng)調(diào)從不同方面評(píng)價(jià)人對(duì)耕地利用的適宜程度并體現(xiàn)區(qū)域耕地不合理利用的主導(dǎo)因素。該文以全國(guó)范圍65個(gè)縣（市）作為研究對(duì)象，應(yīng)用地形特征、土壤性狀、耕作條件、環(huán)境狀況和生物特性5類指標(biāo)整合影響耕地健康的關(guān)鍵因子并評(píng)價(jià)耕地健康狀況，在此基礎(chǔ)上，通過分析指標(biāo)標(biāo)準(zhǔn)差收斂特征與評(píng)價(jià)因子的評(píng)分規(guī)則定義各指標(biāo)的臨界健康與亞健康閾值，將試點(diǎn)縣劃分為健康、亞健康與不健康3類，進(jìn)而分析不同健康狀況的分布特征。同時(shí)，綜合區(qū)域氣候條件、耕地自然條件與耕作技術(shù)水平評(píng)價(jià)耕地產(chǎn)能水平，并分析區(qū)域耕地健康-耕地產(chǎn)能作用特征。結(jié)果表明，試點(diǎn)縣耕地亞健康及不健康類型表現(xiàn)出顯著區(qū)域性特征，不同試點(diǎn)縣限制耕地健康的關(guān)鍵要素存在差異；一方面，耕作條件的亞健康狀況普遍（數(shù)量占比約40%），表現(xiàn)為耕地破碎化嚴(yán)重以及農(nóng)田基礎(chǔ)設(shè)施不完善，也反映出高標(biāo)準(zhǔn)農(nóng)田建設(shè)的迫切需求；另一方面，生物特性形式嚴(yán)峻，表明在氣候條件-土壤結(jié)構(gòu)-耕作技術(shù)等因素相互作用下耕地保障動(dòng)物健康生長(zhǎng)的能力正面臨退化。此外，通過耕地健康與耕地產(chǎn)能耦合分析，該文將試點(diǎn)縣耕地健康產(chǎn)能特征劃分為“健康-高產(chǎn)能”、“健康-低產(chǎn)能”、“亞健康-高產(chǎn)能”、“亞健康-低產(chǎn)能”、“不健康-高產(chǎn)能”、“不健康-低產(chǎn)能”6類，并分析不同類型區(qū)域耕地保護(hù)的重點(diǎn)方向，其中“不健康/亞健康-高產(chǎn)能”區(qū)域分布廣泛，是中國(guó)實(shí)施藏糧于地、藏糧于技戰(zhàn)略的關(guān)鍵地帶，應(yīng)兼顧耕地健康與耕地產(chǎn)能，加強(qiáng)提升耕地健康狀況的方法研究與實(shí)踐。該項(xiàng)研究為耕地健康產(chǎn)能提升提供科學(xué)依據(jù)，為不同區(qū)域耕地資源保護(hù)路徑設(shè)計(jì)提供有效支撐。

土地利用；地理信息系統(tǒng)；耕地；耕地健康；耕地產(chǎn)能；氣候變化

0 引 言

糧食安全問題是關(guān)系人類可持續(xù)發(fā)展基本需求的世界性問題[1]。FAO報(bào)告指出，近年來全球仍有近10億人缺乏飲食保障[2]，并且隨著人口數(shù)量、社會(huì)經(jīng)濟(jì)狀況、氣候條件等因素變化，至2050年全球糧食需求將增長(zhǎng)近70%[3]（也有研究提出糧食需求增長(zhǎng)100%～110%[4]）。2019年8月，IPCC《氣候變化與土地報(bào)告》則指出耕地退化在危及糧食安全同時(shí)也會(huì)加劇氣候變化，可持續(xù)耕地利用對(duì)減少土壤侵蝕，消除饑餓以及應(yīng)對(duì)氣候變化具有重要作用[5]。這為全球科學(xué)家和政策制定者提出了一個(gè)艱巨的挑戰(zhàn)—如何在糧食產(chǎn)量持續(xù)增加的同時(shí)，保護(hù)耕地生態(tài)系統(tǒng)穩(wěn)定，滿足農(nóng)業(yè)和社會(huì)經(jīng)濟(jì)的可持續(xù)發(fā)展[6-7]。而中國(guó)作為世界人口大國(guó)，糧食安全問題長(zhǎng)期以來更是備受全球眾多學(xué)者與研究機(jī)構(gòu)廣泛關(guān)注[8]。事實(shí)證明，近年來即便快速城鎮(zhèn)化與經(jīng)濟(jì)發(fā)展引起優(yōu)質(zhì)耕地資源減少、耕地土壤退化等問題，中國(guó)用全球7%的耕地養(yǎng)活22%的人口，實(shí)現(xiàn)糧食產(chǎn)量的“十三連增”，為實(shí)現(xiàn)聯(lián)合國(guó)千年發(fā)展目標(biāo)做出重要貢獻(xiàn)[9-10]。然而，在糧食產(chǎn)量不斷增長(zhǎng)同時(shí)，化肥農(nóng)藥使用量超標(biāo)、重金屬污染[11]、東北土層變薄[12]、土壤有機(jī)質(zhì)含量下降[13]、地下水位降低[14]等問題也日益加劇，嚴(yán)重制約中國(guó)耕地可持續(xù)利用與生態(tài)安全。因此，評(píng)價(jià)耕地健康產(chǎn)能狀況，分析不同區(qū)域危及耕地健康的關(guān)鍵問題，對(duì)探索可持續(xù)的耕地利用路徑，維護(hù)國(guó)家糧食安全與生態(tài)安全具有重要意義。

從概念與內(nèi)涵角度梳理，國(guó)內(nèi)外與耕地健康產(chǎn)能相關(guān)的研究主要集中在土壤健康、生態(tài)系統(tǒng)健康、耕地產(chǎn)能等方面。在土壤健康方面，1941年，Leopold首先提出土地健康是指土地被人類占用但其功能沒有被破壞的狀態(tài)，土地有機(jī)體具有自我更新能力，并使用土地疾病來描繪土地功能紊亂[15]。Rapport與Sparling等則提出相比于土壤質(zhì)量，土壤健康更側(cè)重于土壤內(nèi)在屬性的狀態(tài)而忽略其對(duì)確切功能的考量（后者正是土壤質(zhì)量評(píng)價(jià)的重點(diǎn)）[16-17]。Hodgson與Kumar等也分別探討了耕地健康的內(nèi)涵[18-19]。在生態(tài)系統(tǒng)健康方面，眾多學(xué)者提出通過活力、組織結(jié)構(gòu)和恢復(fù)力等特征來表達(dá)生態(tài)系統(tǒng)所具有穩(wěn)定性和持續(xù)性[20-22]。在耕地產(chǎn)能方面，眾多學(xué)者認(rèn)為耕地產(chǎn)能即耕地綜合生產(chǎn)能力，是一定時(shí)期、一定經(jīng)濟(jì)與地域條件下形成的耕地生產(chǎn)能力，并研究不同尺度下耕地的理論產(chǎn)能、可實(shí)現(xiàn)產(chǎn)能與實(shí)際產(chǎn)量相互關(guān)系[23-26]。上述研究均已形成相對(duì)完整的概念和內(nèi)涵，并對(duì)耕地健康產(chǎn)能有重要借鑒作用，但未能對(duì)自然-社會(huì)-經(jīng)濟(jì)綜合影響的耕地健康產(chǎn)能進(jìn)行系統(tǒng)性的理論梳理。在上述研究基礎(chǔ)上，近年來中國(guó)學(xué)者開始逐漸關(guān)注耕地健康產(chǎn)能的概念與內(nèi)涵。葉興慶等提出“五毒產(chǎn)能”，歸納了5種不健康的耕地產(chǎn)能挖掘模式[27]。李強(qiáng)等人提出耕地健康內(nèi)涵[28]，并構(gòu)建由耕地質(zhì)量、產(chǎn)能和土壤環(huán)境組成的耕地健康多層次模糊評(píng)價(jià)體系，評(píng)價(jià)研究區(qū)耕地的健康、亞健康和不健康狀態(tài)[29]。本文采用鄖文聚等的觀點(diǎn)[30]，認(rèn)為健康產(chǎn)能是耕地利用過程中各種功能的正常體現(xiàn)，始終具備持續(xù)、穩(wěn)定地進(jìn)行農(nóng)業(yè)生產(chǎn)的能力，具有自我調(diào)節(jié)和對(duì)損害的恢復(fù)能力，在一個(gè)完整的生產(chǎn)周期內(nèi)能夠通過耕地的自凈能力實(shí)現(xiàn)自我更新，且生產(chǎn)的農(nóng)產(chǎn)品不會(huì)對(duì)人體健康造成影響。

從評(píng)價(jià)方法角度梳理，20世紀(jì)60年代以來，美國(guó)、英國(guó)、荷蘭、澳大利亞等國(guó)均開展了土地評(píng)價(jià)方面的研究工作，但多以土地分類和土地潛力分類為主，例如美國(guó)于1961年采用劃分土地利用潛力等級(jí)的方法，從土壤特征出發(fā)進(jìn)行土地潛力評(píng)價(jià)，分為潛力級(jí)、亞潛力級(jí)和潛力單位3級(jí)[31]。1976年，F(xiàn)AO頒布的“土地評(píng)價(jià)綱要”從土地適宜性角度，應(yīng)用綱、類、亞類、單元4級(jí)評(píng)價(jià)土地生產(chǎn)的適宜性程度、限制性因素與管理需求，并被廣泛應(yīng)用于各國(guó)土地評(píng)價(jià)項(xiàng)目[32-35]。至80年代末，隨著人口增長(zhǎng)、土地退化和環(huán)境問題的日益加劇，土地可持續(xù)利用問題逐漸成為該領(lǐng)域研究重點(diǎn)，而“土地評(píng)價(jià)綱要”由于缺乏空間、過程和動(dòng)態(tài)分析，未能指出土地利用持續(xù)性問題，逐漸無法滿足應(yīng)用需求[36-38]。1993年，F(xiàn)AO制定“可持續(xù)土地利用管理大綱”，從生產(chǎn)、穩(wěn)定、可行、可承受、保護(hù)性進(jìn)行土地可持續(xù)評(píng)價(jià)[39]。2001年起，康奈爾土壤健康團(tuán)隊(duì)綜合考慮數(shù)據(jù)獲取的難易程度與農(nóng)戶的可接受度，確定了13項(xiàng)指標(biāo)（包括物理指標(biāo)5項(xiàng)，生物指標(biāo)4項(xiàng)，化學(xué)指標(biāo)4項(xiàng)）綜合評(píng)定土壤健康狀況[40]。Vincent等人應(yīng)用縮減回歸（RR），主成分回歸（PCR），偏最小二乘回（PLSR），等方法分析了10種土壤理化屬性指標(biāo)的相互作用特征，并提出土壤有機(jī)碳，土壤容重，碳氮比與導(dǎo)電率是影響土壤質(zhì)量的最主要變量[41]。然而，上述研究更多關(guān)注“土壤在生態(tài)系統(tǒng)和土地利用范圍內(nèi)發(fā)揮作用以維持生物生產(chǎn)能力、維持環(huán)境質(zhì)量、促進(jìn)植物和動(dòng)物健康的能力”[42]，缺乏對(duì)影響耕地可持續(xù)利用的各種核心要素的綜合分析。2016年，美國(guó)自然保護(hù)協(xié)會(huì)針對(duì)土壤退化嚴(yán)重問題制定土壤健康路線圖[43]，從科學(xué)，經(jīng)濟(jì)與政策3方面提出建議，以期推動(dòng)美國(guó)農(nóng)田管理模式轉(zhuǎn)變，為耕地保護(hù)與利用提供指導(dǎo)。該研究將增強(qiáng)氣候-土壤-耕作系統(tǒng)綜合認(rèn)知分析作為未來重要發(fā)展方向。此外，由于中國(guó)幅員遼闊，區(qū)域氣候-土壤-基礎(chǔ)設(shè)施-自然資源-社會(huì)經(jīng)濟(jì)條件差異化特征復(fù)雜，也提升了耕地評(píng)價(jià)的難度。2008年，農(nóng)業(yè)部以指導(dǎo)農(nóng)業(yè)生產(chǎn)、提升耕地產(chǎn)量為目標(biāo)頒布《耕地地力調(diào)查與質(zhì)量評(píng)價(jià)技術(shù)規(guī)程》[44]，從耕地地力和土壤環(huán)境質(zhì)量2個(gè)方面評(píng)價(jià)耕地質(zhì)量，取得良好的應(yīng)用成果，但一定程度上未能體現(xiàn)氣候條件與基礎(chǔ)設(shè)施條件對(duì)農(nóng)業(yè)生產(chǎn)的影響。原國(guó)土資源部也先后頒布《農(nóng)用地質(zhì)量分等規(guī)程》[45]與《土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范》[46]，以耕地資源保護(hù)為目標(biāo)系統(tǒng)全面地調(diào)查農(nóng)用地產(chǎn)能狀況、土壤養(yǎng)分元素、土壤重金屬環(huán)境、灌溉水環(huán)境、大氣干濕沉降物等指標(biāo)，但存在相關(guān)成果缺乏統(tǒng)一梳理整合，無法反映土壤動(dòng)物活性等問題。

綜上所述，耕地健康產(chǎn)能是耕地自然系統(tǒng)與利用系統(tǒng)綜合作用的結(jié)果，其中耕地產(chǎn)能較大程度上是由自然系統(tǒng)主導(dǎo)的（如光照、溫度、降水、土壤理化性質(zhì)等），而耕地健康問題則主要由于人類的不合理利用。并且由于耕地健康產(chǎn)能的內(nèi)涵豐富，影響耕地健康產(chǎn)能的因素眾多且作用過程復(fù)雜，直接計(jì)算耕地健康產(chǎn)能狀況是困難的。因此，本文嘗試采用一種簡(jiǎn)化的方法評(píng)價(jià)耕地健康產(chǎn)能，即分別評(píng)價(jià)耕地健康與耕地產(chǎn)能狀況，并通過耦合耕地健康與產(chǎn)能分析耕地利用合理程度與生產(chǎn)能力之間的相互作用特征。其中耕地健康狀況強(qiáng)調(diào)從不同方面評(píng)價(jià)人對(duì)耕地利用的適宜程度并體現(xiàn)區(qū)域耕地不合理利用的主導(dǎo)因素。在此基礎(chǔ)上，本文根據(jù)標(biāo)準(zhǔn)耕作制度[47]差異在全國(guó)范圍選取65個(gè)縣（市）作為研究區(qū)域，從地形特征、土壤性狀、耕作條件、環(huán)境狀況和生物特性5個(gè)維度梳理影響耕地健康的關(guān)鍵因子并評(píng)價(jià)耕地健康狀況，綜合區(qū)域氣候條件、耕地自然條件與耕作技術(shù)水平評(píng)價(jià)耕地產(chǎn)能水平，并分析區(qū)域耕地健康-耕地產(chǎn)能作用特征，以期為不同區(qū)域耕地資源保護(hù)路徑設(shè)計(jì)提供支持。

1 方 法

1.1 耕地健康評(píng)價(jià)方法

耕地健康格局由包含土壤結(jié)構(gòu)、耕作技術(shù)、土壤養(yǎng)分含量、污染狀況、生物群落狀態(tài)等在內(nèi)的多因素間復(fù)雜相互作用形成，并且這種復(fù)雜相互作用可能由于區(qū)域水熱條件與自然-社會(huì)-經(jīng)濟(jì)格局差異而表現(xiàn)出不同的特征。因此，耕地健康保護(hù)與提升需要厘清影響耕地可持續(xù)利用的關(guān)鍵因素。本文耕地健康評(píng)價(jià)在溫良友[48-49]等耕地質(zhì)量評(píng)價(jià)體系基礎(chǔ)上面向不同區(qū)域分異特征梳理修正評(píng)價(jià)指標(biāo)（表1），應(yīng)用地形特征系數(shù)、土壤性狀系數(shù)、耕作條件系數(shù)、環(huán)境狀況系數(shù)和生物特性系數(shù)5類指標(biāo)表達(dá)人對(duì)耕地利用的合理程度，評(píng)價(jià)單元為耕地地塊圖斑。其中，地形特征旨在從地形形態(tài)對(duì)水肥供給與運(yùn)移的影響角度評(píng)價(jià)耕地地理位置的合理度；土壤性狀綜合體現(xiàn)耕地土壤的物理和化學(xué)特性對(duì)作物生長(zhǎng)的適宜性；耕作條件體現(xiàn)人對(duì)耕地進(jìn)行的基礎(chǔ)設(shè)施建設(shè)和改造狀況；環(huán)境狀況表達(dá)耕地面臨的污染程度；生物特性表達(dá)耕地利用對(duì)動(dòng)物健康生長(zhǎng)的影響，應(yīng)用物種指示方法[50]，通過單位面積土壤蚯蚓數(shù)量計(jì)算。

表1 耕地健康評(píng)價(jià)指標(biāo)體系（除環(huán)境狀況）與賦分規(guī)則

在實(shí)際評(píng)價(jià)過程中，本研究首先依據(jù)區(qū)域水熱條件，以“不損傷耕地利用可持續(xù)性前提下最大限度發(fā)揮耕地生產(chǎn)潛力”為標(biāo)準(zhǔn)制定不同縣區(qū)的適宜種植制度[47]。隨后應(yīng)用德爾菲法針對(duì)不同二級(jí)分區(qū)遴選區(qū)域備選因子并設(shè)計(jì)各評(píng)價(jià)因子權(quán)重與評(píng)價(jià)指標(biāo)權(quán)重；同時(shí)我們預(yù)先結(jié)合經(jīng)驗(yàn)知識(shí)定義各評(píng)價(jià)因子的評(píng)分規(guī)則（表1），以將不同量綱的評(píng)價(jià)因子觀測(cè)/調(diào)查數(shù)據(jù)轉(zhuǎn)換到同一標(biāo)準(zhǔn)[51-52]。其中地形特征，土壤性狀與耕作條件通過評(píng)價(jià)因子加權(quán)平均計(jì)算；環(huán)境狀況系數(shù)計(jì)算在應(yīng)用單項(xiàng)污染指數(shù)法分別評(píng)價(jià)土壤與灌溉水中某種元素的污染程度基礎(chǔ)上采用“1+土+水”累加模型突出高濃度污染元素對(duì)耕地環(huán)境的影響。如式（1）所示

式中土和水分別為土壤重金屬污染與灌溉水環(huán)境影響系數(shù)；表示《土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)GB15618—2008》[53]中砷（As）、鉻（Cr）、鎘（Cd）、銅（Cu）、汞（Hg）、鉛（Pb）、鎳（Ni）和鋅（Zn）8個(gè)重金屬元素指標(biāo)的實(shí)測(cè)含量，mg/kg；則為GB15618中規(guī)定的相應(yīng)指標(biāo)的標(biāo)準(zhǔn)值，mg/kg；P表示對(duì)應(yīng)分值（表2），當(dāng)小于或等于1時(shí)，評(píng)價(jià)指標(biāo)符合標(biāo)準(zhǔn)，等級(jí)為1級(jí)，P為0；當(dāng)大于1時(shí)，指標(biāo)含量高于GB15618標(biāo)準(zhǔn)值，P隨著增大而減??；取值為P最小值,且當(dāng)大于5時(shí)耕地退出使用，環(huán)境狀況系數(shù)為0。同理，與分別表示《農(nóng)田灌溉水質(zhì)標(biāo)準(zhǔn)GB5084—2005》[54]中水質(zhì)評(píng)價(jià)指標(biāo)（例如總砷、六價(jià)鉻、總鎘、總汞、總鉛等）的實(shí)測(cè)含量與標(biāo)準(zhǔn)值，mg/L，當(dāng)水j大于1時(shí)，指標(biāo)含量高于標(biāo)準(zhǔn)值，等級(jí)為2級(jí)，P為-0.1。生物特性則計(jì)算“土壤蚯蚓數(shù)量對(duì)應(yīng)分值/100”。

表2 耕地環(huán)境狀況評(píng)價(jià)指標(biāo)體系與賦分規(guī)則

1.2 耕地產(chǎn)能評(píng)價(jià)方法

耕地產(chǎn)能是指一定時(shí)期內(nèi)一個(gè)地區(qū)在一定經(jīng)濟(jì)技術(shù)條件下，由各生產(chǎn)要素特別是長(zhǎng)期要素（自然條件，技術(shù)條件）綜合投入所形成的、可以穩(wěn)定實(shí)現(xiàn)一定產(chǎn)量的糧食產(chǎn)出能力[55]。產(chǎn)量是指實(shí)際種植的糧食作物，在應(yīng)對(duì)作物生長(zhǎng)期各種自然、經(jīng)濟(jì)、社會(huì)、技術(shù)條件變化下，在耕地面積上實(shí)現(xiàn)的實(shí)際糧食產(chǎn)出數(shù)量。從邏輯關(guān)系上看，耕地產(chǎn)能是本，糧食產(chǎn)量是相。耕地產(chǎn)能是糧食產(chǎn)量的前提和基礎(chǔ)，只有耕地產(chǎn)能得到根本提升和良好保育，糧食產(chǎn)量才有可能提高。產(chǎn)量是產(chǎn)能的有效釋放和挖掘，產(chǎn)量接近產(chǎn)能或高于產(chǎn)能，表現(xiàn)為耕地過渡挖掘利用；產(chǎn)量遠(yuǎn)低于產(chǎn)能，則表現(xiàn)為耕地利用不足，如圖1示。張晉科等通過對(duì)比耕地生產(chǎn)能力與糧食現(xiàn)實(shí)產(chǎn)量分析全國(guó)耕地利用程度，提出利用程度高的區(qū)域主要分布在東南沿海省份、東北省份以及自然條件差的青海、西藏兩省，而利用低的省份主要分布在中國(guó)北部以及南部的云貴高原區(qū)[56]。

圖1 耕地產(chǎn)能與作物產(chǎn)量關(guān)系

本文耕地產(chǎn)能評(píng)價(jià)通過綜合區(qū)域氣候條件系數(shù)，耕地自然條件系數(shù)與技術(shù)水平系數(shù)表達(dá)天地相互作用形成的耕地生產(chǎn)力水平和產(chǎn)出效果。在實(shí)際評(píng)價(jià)過程中，本研究一方面以種植制度為參考，將物候觀測(cè)數(shù)據(jù)與氣象數(shù)據(jù)導(dǎo)入AEZ模型[57-58]，計(jì)算不同縣級(jí)行政區(qū)相關(guān)作物對(duì)應(yīng)的光溫/氣候生產(chǎn)潛力與產(chǎn)量比系數(shù)。其中根據(jù)耕地灌溉條件區(qū)分灌溉農(nóng)業(yè)與雨養(yǎng)農(nóng)業(yè)，相應(yīng)計(jì)算光溫生產(chǎn)潛力或氣候生產(chǎn)潛力，并應(yīng)用試驗(yàn)田產(chǎn)量數(shù)據(jù)修正；用于將不同作物的單產(chǎn)轉(zhuǎn)換到同一種作物。另一方面應(yīng)用德爾菲法遴選區(qū)域備選因子并設(shè)計(jì)各評(píng)價(jià)因子權(quán)重（表3），并預(yù)先定義各評(píng)價(jià)因子的評(píng)分規(guī)則（表1），以將不同量綱的評(píng)價(jià)因子觀測(cè)/調(diào)查數(shù)據(jù)轉(zhuǎn)換到同一標(biāo)準(zhǔn)。耕地產(chǎn)能指數(shù)計(jì)算公式如下

表3 耕地產(chǎn)能評(píng)價(jià)指標(biāo)體系

2 試驗(yàn)數(shù)據(jù)

3 試驗(yàn)與分析

3.1 耕地健康格局分析

本文嘗試從5種評(píng)價(jià)指標(biāo)維度，分析限制耕地健康的關(guān)鍵因素。圖2a～2e分別列出了縣級(jí)地形特征系數(shù)、土壤性狀系數(shù)、耕作條件系數(shù)、環(huán)境狀況系數(shù)與生物特性系數(shù)平均值與其相應(yīng)標(biāo)準(zhǔn)差的相關(guān)性特征。如圖所示，一方面，各系數(shù)平均值與標(biāo)準(zhǔn)差無顯著的線性相關(guān)特征，2均小于0.48；另一方面，隨著平均值增大，各系數(shù)標(biāo)準(zhǔn)差均趨于減小，這表明對(duì)于某種系數(shù)平均水平較高的試點(diǎn)縣，其縣內(nèi)不同耕地地塊對(duì)應(yīng)的該系數(shù)的離散程度低，表現(xiàn)出更高的一致性。在此基礎(chǔ)上，本文將這種標(biāo)準(zhǔn)差收斂特征與各評(píng)價(jià)因子的評(píng)分規(guī)則相結(jié)合，定義各系數(shù)的臨界健康閾值與臨界亞健康閾值，如圖2f所示。以耕作條件系數(shù)為例，其平均值大于75的樣點(diǎn)，相應(yīng)標(biāo)準(zhǔn)差顯著收斂，樣點(diǎn)分布表現(xiàn)出更高的聚集特征（即高頻特征）；并且隨著平均值減小，樣點(diǎn)頻率逐漸降低，高頻數(shù)據(jù)與低頻數(shù)據(jù)在平均值約為50處產(chǎn)生間斷（這種間斷可以視為自然分類法的基礎(chǔ)[59]）；此外，相關(guān)評(píng)分規(guī)則中應(yīng)用分值80（或90）表征評(píng)價(jià)因子較好狀態(tài)，分值60（或70）表征中等狀態(tài)，分值40（或50）表征較差狀態(tài)；綜上，本文定義耕作條件系數(shù)的臨界健康閾值為75，臨界亞健康閾值為50。

a. 地形特征系數(shù)相關(guān)性分析 a. Correlation analysis of topographic feature coefficientb. 土壤性狀系數(shù)相關(guān)性分析 b. Correlation analysis of soil property coefficientc. 耕作條件系數(shù)相關(guān)性分析 c. Correlation analysis of farming condition coefficient

d. 環(huán)境狀況系數(shù)相關(guān)性分析 d. Correlation analysis of environmental status coefficiente. 生物特性系數(shù)相關(guān)性分析 e. Correlation analysis of biological property coefficientf. 臨界健康-臨界亞健康閥值設(shè)計(jì) f. Critical threshold of health and sub-health status

如圖3所示，本文依據(jù)臨界健康閾值與臨界亞健康閾值將試點(diǎn)縣劃分為健康（A類），亞健康（B類）與不健康（C類）3類，對(duì)于5種系數(shù)均大于（或等于）臨界健康閾值的試點(diǎn)縣歸于A類；對(duì)于存在一個(gè)或多個(gè)系數(shù)小于臨界亞健康閾值的試點(diǎn)縣歸于C類，其它則歸于B類。對(duì)于A類試點(diǎn)縣，影響耕地可持續(xù)利用的關(guān)鍵因素均處于良好水平，耕地安全在一定時(shí)期內(nèi)得以保障；對(duì)于B類試點(diǎn)縣，部分因素由于耕地不合理利用而逐漸退化并成為耕地健康隱患，應(yīng)側(cè)重厘清制約耕地利用的短板與驅(qū)動(dòng)因子，及其預(yù)防與提升策略；對(duì)于C類試點(diǎn)縣，耕地健康面臨嚴(yán)重威脅，應(yīng)以治理為主，通過土地整治工程與減小耕地利用強(qiáng)度（例如休耕）等方法改善耕地健康狀況，并加強(qiáng)耕地健康評(píng)價(jià)指標(biāo)監(jiān)測(cè)。在空間分布特征方面，A類與B類試點(diǎn)縣多分布在中國(guó)東部地區(qū)，例如長(zhǎng)江中下游區(qū)，江南區(qū)，東北區(qū)以及黃淮海區(qū)等；其中A類試點(diǎn)縣集中分布在東部平原地區(qū)，如黃淮平原，江淮平原，沿江平原，松嫩平原等；C類試點(diǎn)縣橫跨中國(guó)北部地區(qū)，在云貴高原與四川盆地交界地帶以及南部沿海地區(qū)也存在樣例，并且相關(guān)地形顯著為山地或丘陵；而B類試點(diǎn)縣作為A類與C類間的過渡類型，分布最為廣泛，在山地區(qū)與平原區(qū)均有樣例。

圖4a～4c分別展示3類試點(diǎn)縣對(duì)應(yīng)的評(píng)價(jià)指標(biāo)分布特征，圖4d則統(tǒng)計(jì)不同評(píng)價(jià)指標(biāo)健康、亞健康以及不健康狀態(tài)對(duì)應(yīng)的試點(diǎn)縣數(shù)量。如圖所示，不同試點(diǎn)縣限制耕地健康的關(guān)鍵要素表現(xiàn)出顯著區(qū)域性特征，耕地健康保護(hù)路徑研究需要面向區(qū)域?qū)嶋H“短板”問題。對(duì)于不同類型試點(diǎn)縣，生物特性系數(shù)分布離散程度均顯著高于其它評(píng)價(jià)指標(biāo)；大部分B類試點(diǎn)縣的“短板”在于生物特性系數(shù)與耕作條件系數(shù)，其它評(píng)價(jià)指標(biāo)基本處于健康狀態(tài)；而大部分C類試點(diǎn)縣的生物特性系數(shù)低于亞健康臨界閾值，其它評(píng)價(jià)指標(biāo)低于亞健康臨界閾值的現(xiàn)象很少。由圖4d可見，本研究中存在地形特征、土壤性狀與環(huán)境狀況處于亞健康或不健康狀態(tài)的試點(diǎn)縣樣例，但這種現(xiàn)象并不普遍。比較而言，一方面，耕作條件的亞健康狀況更為普遍，表現(xiàn)為耕地破碎化嚴(yán)重以及農(nóng)田基礎(chǔ)設(shè)施不完善，也反映出高標(biāo)準(zhǔn)農(nóng)田建設(shè)的迫切需求。另一方面，生物特性形式嚴(yán)峻，表明耕地保障動(dòng)物健康生長(zhǎng)的能力正面臨退化，這可能與近40a來耕地利用強(qiáng)度持續(xù)加大，農(nóng)藥、化肥、除草劑大量使用，農(nóng)機(jī)裝備趨于中重型，缺乏深耕細(xì)作及土壤有機(jī)質(zhì)含量下降有關(guān)，也與蚯蚓數(shù)量區(qū)域分異特征以及缺乏長(zhǎng)期觀測(cè)實(shí)踐經(jīng)驗(yàn)有一定聯(lián)系。此外，通過土地平整、規(guī)?；?jīng)營(yíng)、加強(qiáng)基礎(chǔ)設(shè)施建設(shè)等能夠在較長(zhǎng)期間內(nèi)有效提升并維持耕作條件水平，而生物特性狀況面臨成因復(fù)雜、觀測(cè)體系不完善、改良困難、變化性強(qiáng)等難題，亟需進(jìn)一步完善評(píng)價(jià)指標(biāo)與觀測(cè)評(píng)價(jià)方法，加強(qiáng)多因素綜合作用過程認(rèn)知研究。

圖3 試點(diǎn)縣健康-亞健康-不健康類型空間分布特征

a. A類試點(diǎn)縣系數(shù)分布特征 a. Coefficient characteristics of counties belong to class Ab. B類試點(diǎn)縣系數(shù)分布特征 b. Coefficient characteristics of counties belong to class Bc. C類試點(diǎn)縣系數(shù)分布特征 c. Coefficient characteristics of counties belong to class Cd. 不同健康狀態(tài)對(duì)應(yīng)試點(diǎn)縣數(shù)量 d. Pilot county count of different health status

3.2 耕地產(chǎn)能格局分析

圖5 試點(diǎn)縣耕地產(chǎn)能格局

3.3 耕地健康-產(chǎn)能耦合分析

耕地健康與耕地產(chǎn)能各有側(cè)重又緊密聯(lián)系：一方面，耕地健康側(cè)重體現(xiàn)區(qū)域耕地不合理利用的主導(dǎo)因素，耕地產(chǎn)能則側(cè)重刻畫天地綜合作用形成的耕地生產(chǎn)力水平；另一方面，耕地健康是影響耕地產(chǎn)能的重要因素，不健康的耕地往往較大程度上限制作物生長(zhǎng)過程對(duì)區(qū)域水-熱-肥條件的利用，而對(duì)耕地產(chǎn)能的過度挖掘也易于引發(fā)耕地健康問題。因此分析耕地健康與耕地產(chǎn)能的關(guān)聯(lián)特征，對(duì)認(rèn)知區(qū)域耕地保護(hù)的重點(diǎn)方向具有重要參考作用。

本文通過建立耕地健康類型與耕地產(chǎn)能指數(shù)的組合關(guān)系，將試點(diǎn)縣耕地健康產(chǎn)能特征劃分為“健康-高產(chǎn)能”、“健康-低產(chǎn)能”、“亞健康-高產(chǎn)能”、“亞健康-低產(chǎn)能”、“不健康-高產(chǎn)能”、“不健康-低產(chǎn)能”6類，并統(tǒng)計(jì)不同類型試點(diǎn)縣對(duì)應(yīng)歸一化耕地平均產(chǎn)能以及地形特征、土壤性狀、耕作條件、環(huán)境狀況、生物特性平均值，如圖6所示。對(duì)比統(tǒng)計(jì)結(jié)果，健康（-低產(chǎn)能/高產(chǎn)能）類耕地相較于同等產(chǎn)能水平的不健康、亞健康耕地顯著表現(xiàn)出更高的平均產(chǎn)能，反映出耕地健康對(duì)耕地產(chǎn)能的正向促進(jìn)作用；健康/亞健康-高/低產(chǎn)能4種類型試點(diǎn)縣對(duì)于生物特性指標(biāo)差異大，而對(duì)其他4項(xiàng)指標(biāo)差異較小，反映出生物特性指標(biāo)對(duì)耕地健康-產(chǎn)能耦合關(guān)系的指示性作用；相比不健康-低產(chǎn)能類型，不健康-高產(chǎn)能類型的土壤性狀與生物特性水平更低，反映出不健康狀態(tài)下產(chǎn)能高的耕地面臨更大的退化風(fēng)險(xiǎn)。此外，不健康-低產(chǎn)能類型的地形特征，耕作條件與環(huán)境狀況均顯著較低，這可能表現(xiàn)出該類耕地?zé)o序開墾、投入不足、粗放利用的問題。“健康-高產(chǎn)能”地區(qū)在長(zhǎng)江中下游及黃淮海南部地區(qū)形成聚集，是中國(guó)糧食核心產(chǎn)區(qū)的重要體現(xiàn)，同時(shí)也是重點(diǎn)防護(hù)區(qū)域，需要加強(qiáng)觀測(cè)以防止耕地過度利用引發(fā)土壤、環(huán)境、資源、生態(tài)等方面問題。“亞健康-高產(chǎn)能”地區(qū)具有良好的糧食產(chǎn)出能力，對(duì)于中國(guó)糧食安全有重要影響，是中國(guó)實(shí)施藏糧于地、藏糧于技戰(zhàn)略的關(guān)鍵地帶；但不合理的利用已經(jīng)危及耕地可持續(xù)性，表現(xiàn)出耕地健康-產(chǎn)能不平衡特征，因此，對(duì)于這類區(qū)域耕地保護(hù)的重點(diǎn)方向應(yīng)兼顧耕地健康與耕地產(chǎn)能，加強(qiáng)提升耕地健康狀況的方法研究與實(shí)踐，同時(shí)以不惡化健康狀況為目標(biāo)研究合理的耕地利用強(qiáng)度?！安唤】?低產(chǎn)能”類型主要分布在中國(guó)北部地區(qū)，對(duì)于這類區(qū)域應(yīng)重點(diǎn)關(guān)注耕地健康修復(fù)，例如通過休耕、輪作等方式修復(fù)土壤性狀問題嚴(yán)重區(qū)域或生態(tài)嚴(yán)重退化地區(qū)。

注：耕地健康產(chǎn)能類型分別為：不健康-低產(chǎn)能類型、不健康-高產(chǎn)能類型、亞健康-低產(chǎn)能類型、亞健康-高產(chǎn)能類型、健康-低產(chǎn)能類型、健康-高產(chǎn)能類型。

4 討 論

耕地健康產(chǎn)能是耕地自然系統(tǒng)與利用系統(tǒng)綜合作用的結(jié)果，既包含影響耕地健康格局的自然環(huán)境要素和社會(huì)經(jīng)濟(jì)要素，也包含人類利用導(dǎo)致的不穩(wěn)定要素。本文在整合已有評(píng)價(jià)指標(biāo)體系與評(píng)價(jià)方法成果基礎(chǔ)上，計(jì)算全國(guó)65個(gè)試點(diǎn)縣影響耕地健康的關(guān)鍵要素格局，并通過耦合耕地健康與耕地產(chǎn)能格局分析兩者間相互作用特征，以期為不同區(qū)域耕地資源保護(hù)路徑設(shè)計(jì)提供信息支持，但仍存在不足：

1）生物特性狀況是氣候條件-土壤結(jié)構(gòu)-土壤養(yǎng)分-耕作技術(shù)等多因素作用的結(jié)果，面臨成因復(fù)雜、改良困難、變化性強(qiáng)等難題。蚯蚓活動(dòng)對(duì)提高和保持土壤肥力，土壤結(jié)構(gòu)和團(tuán)聚體穩(wěn)定性具有重要作用[60-61]。實(shí)驗(yàn)表明，草地中蚯蚓的缺失增加了土壤容重和抗剪強(qiáng)度，進(jìn)而大大降低了土壤有機(jī)質(zhì)、土壤水分和入滲率[62]。泥漿，氨基無機(jī)肥料，酸化肥料，有機(jī)磷農(nóng)藥以及高強(qiáng)度耕作均可能對(duì)蚯蚓的繁殖或生長(zhǎng)產(chǎn)生不利影響[63-69]。同時(shí)，土壤微生物群落是創(chuàng)造和保持良好土壤結(jié)構(gòu)的關(guān)鍵，例如通過改善通風(fēng)和腐殖質(zhì)影響顆粒團(tuán)聚體的形成。絲狀真菌和放線菌尤其能吸附土壤顆粒，使其形成團(tuán)聚體。此外，細(xì)胞外代謝物，如多糖、脂質(zhì)和蛋白質(zhì)，具有穩(wěn)定聚集物的功能[70-73]。但由于計(jì)算微生物物種多樣性的難度，以及解釋土壤呼吸測(cè)試的難度而應(yīng)用較少[74-76]。此外，研究表明土壤有機(jī)碳是土壤生物群落動(dòng)力學(xué)的最佳代表，在農(nóng)田的肥力、土壤水分有效性和團(tuán)聚體穩(wěn)定性方面發(fā)揮著關(guān)鍵作用[77]。本文從數(shù)據(jù)易獲取性角度考慮，選擇土壤蚯蚓數(shù)量表征該項(xiàng)指標(biāo)，存在評(píng)價(jià)指標(biāo)單一、觀測(cè)經(jīng)驗(yàn)不足、觀測(cè)體系不完善、觀測(cè)數(shù)據(jù)積累較少等不足，因此在未來的研究中亟需進(jìn)一步加強(qiáng)多因素作用過程認(rèn)知研究，研究微生物群落、有益微生物種群數(shù)等指標(biāo)作用特征，完善評(píng)價(jià)指標(biāo)與觀測(cè)調(diào)查技術(shù)。

2）耕地系統(tǒng)是一個(gè)由多種要素組成、組織結(jié)構(gòu)復(fù)雜、演化方向多樣，具有區(qū)域差異特征、區(qū)域綜合特征的復(fù)雜巨系統(tǒng)。以資源可持續(xù)利用為目標(biāo)的耕地健康產(chǎn)能評(píng)價(jià)研究已經(jīng)進(jìn)入系統(tǒng)階段。傳統(tǒng)的多要素加權(quán)平均方法沿襲自然科學(xué)研究“還原論”經(jīng)典的思維邏輯，因在很大程度上忽略了要素的內(nèi)在關(guān)聯(lián)和相互作用關(guān)系，不能表達(dá)耕地整體特征而逐漸無法滿足應(yīng)用需求[78-79]。因此，耕地健康產(chǎn)能評(píng)價(jià)指標(biāo)綜合方法研究需要在已有數(shù)據(jù)成果基礎(chǔ)上，促發(fā)從模式思維向數(shù)據(jù)思維轉(zhuǎn)變，通過研究評(píng)價(jià)因子空間分異特征及因子間相互作用過程認(rèn)知耕地系統(tǒng)變化規(guī)律，創(chuàng)建耦合的土地計(jì)算方法體系。

3）在耕地健康評(píng)價(jià)指標(biāo)設(shè)計(jì)方面，本文從地形特征、土壤性狀、耕作條件、環(huán)境狀況和生物特性5個(gè)方面計(jì)算表達(dá)人對(duì)耕地利用的適宜程度，這與國(guó)外土壤健康、土壤質(zhì)量評(píng)價(jià)研究側(cè)重土壤物理、化學(xué)與生物性質(zhì)[80]是一致的。此外考慮到耕地侵占、破碎化嚴(yán)重、技術(shù)水平不足等中國(guó)特色問題，本文進(jìn)一步集成地形特征與耕作條件方面的評(píng)價(jià)。本文試點(diǎn)研究表明耕作條件的亞健康狀況較為普遍，反映出高標(biāo)準(zhǔn)農(nóng)田建設(shè)的迫切需求。耕地健康內(nèi)涵豐富，涉及因素眾多，例如耕作方式（輪作/保護(hù)性耕作/免耕/覆蓋作物）、利用強(qiáng)度、水資源承載力、農(nóng)作物安全、大氣環(huán)境等均是影響耕地健康的重要因素。后續(xù)研究中我們將完善評(píng)價(jià)指標(biāo)體系，更全面表達(dá)耕地健康狀況。

4）在評(píng)價(jià)因子數(shù)據(jù)獲取方面，本研究以縣級(jí)行政區(qū)為單元開展評(píng)價(jià)研究，相關(guān)數(shù)據(jù)多來源于已有資料成果（例如縣級(jí)土壤普查成果、縣級(jí)地球化學(xué)調(diào)查成果、縣級(jí)耕地質(zhì)量等級(jí)評(píng)價(jià)成果、縣級(jí)水利普查數(shù)據(jù)等），存在數(shù)據(jù)年份不一致的現(xiàn)象；少量數(shù)據(jù)來源于實(shí)地采樣檢測(cè)或外業(yè)調(diào)查（例如土壤理化性質(zhì)檢測(cè)，灌溉保障程度外業(yè)調(diào)查等），存在數(shù)據(jù)采集人工成本高、數(shù)據(jù)更新周期長(zhǎng)等問題；此外，部分評(píng)價(jià)因子（例如災(zāi)害防治水平和農(nóng)藝管理水平）的判定標(biāo)準(zhǔn)以定性描述為主，數(shù)據(jù)采集難度大，并受外業(yè)人員主觀認(rèn)識(shí)、知識(shí)水平及對(duì)當(dāng)?shù)厍闆r了解程度影響較大，不利于調(diào)查成果的橫向?qū)Ρ?。因此，需要進(jìn)一步研究基于遙感與地面?zhèn)鞲芯W(wǎng)絡(luò)的評(píng)價(jià)因子數(shù)據(jù)采集獲取方法與技術(shù)[81-84]，發(fā)展?jié)M足客觀快速觀測(cè)需求的評(píng)價(jià)因子指標(biāo)體系，加強(qiáng)數(shù)據(jù)獲取準(zhǔn)確性、時(shí)效性與統(tǒng)一可比性。

5）閾值確定是國(guó)內(nèi)外生態(tài)系統(tǒng)健康、土壤健康、土壤質(zhì)量等方面研究廣泛面臨的難題。相關(guān)方法包括基于專業(yè)知識(shí)的德爾菲法，基于指標(biāo)空間格局的多項(xiàng)式擬合方法[85]以及基于長(zhǎng)期觀測(cè)數(shù)據(jù)的機(jī)理模型方法[86]，其中機(jī)理模型方法最為客觀，但由于機(jī)理過程復(fù)雜在實(shí)際應(yīng)用中往往針對(duì)單個(gè)評(píng)價(jià)因子，在多因子綜合研究中難以應(yīng)用。本文對(duì)于耕地健康-亞健康臨界閾值確定嘗試采用將前兩種方法相結(jié)合，以專家知識(shí)為主導(dǎo)并應(yīng)用指標(biāo)標(biāo)準(zhǔn)差收斂特征修正，目的是將地塊尺度的因子內(nèi)在屬性評(píng)分與縣域尺度指標(biāo)一致性特征結(jié)合，但由于缺乏機(jī)理模型而存在可解釋性不足的問題。因此，在未來亟需進(jìn)一步加強(qiáng)評(píng)價(jià)指標(biāo)狀態(tài)-壓力-響應(yīng)過程研究，發(fā)展耕地健康評(píng)價(jià)閾值確定方法。

5 結(jié) 論

本文以全國(guó)范圍65個(gè)縣（市）作為研究區(qū)域，應(yīng)用地形特征、土壤性狀、耕作條件、環(huán)境狀況和生物特性5類指標(biāo)整合影響耕地健康的關(guān)鍵因素并評(píng)價(jià)耕地健康狀況，在此基礎(chǔ)上通過分析指標(biāo)標(biāo)準(zhǔn)差收斂特征與評(píng)價(jià)因子的評(píng)分規(guī)則，定義各指標(biāo)的臨界健康閾值與臨界亞健康閾值，將試點(diǎn)縣劃分為健康、亞健康與不健康3類，進(jìn)而分析不同健康狀況的分布特征。同時(shí)，綜合區(qū)域氣候條件、耕地自然條件與耕作技術(shù)水平評(píng)價(jià)耕地產(chǎn)能水平，并分析研究區(qū)域耕地健康-耕地產(chǎn)能作用特征。研究表明，試點(diǎn)縣耕地亞健康及不健康類型表現(xiàn)出顯著區(qū)域性特征，不同試點(diǎn)縣限制耕地健康的關(guān)鍵要素存在差異；一方面，耕作條件的亞健康狀況更為普遍，表現(xiàn)為耕地破碎化嚴(yán)重以及農(nóng)田基礎(chǔ)設(shè)施不完善，也反映出高標(biāo)準(zhǔn)農(nóng)田建設(shè)的迫切需求；另一方面，生物特性形式嚴(yán)峻，表明在氣候條件-土壤結(jié)構(gòu)-耕作方式等因素相互作用下耕地保障動(dòng)物健康生長(zhǎng)的能力正面臨退化。此外，通過耕地健康與耕地產(chǎn)能耦合分析，本文將試點(diǎn)縣耕地健康產(chǎn)能特征劃分為“健康-高產(chǎn)能”、“健康-低產(chǎn)能”、“亞健康-高產(chǎn)能”、“亞健康-低產(chǎn)能”、“不健康-高產(chǎn)能”、“不健康-低產(chǎn)能”6類，并分析不同類型區(qū)域耕地保護(hù)的重點(diǎn)方向。

致 謝

感謝張鳳榮，孔祥斌，吳克寧，高培超，沈石對(duì)本文指標(biāo)體系的討論與建議。感謝北京師范大學(xué)地理科學(xué)學(xué)部地理數(shù)據(jù)與應(yīng)用分析中心地學(xué)高性能計(jì)算平臺(tái)支持[https://gda.bnu.edu.cn/]。

[1]Rosegrant M W, Cline S A. Global food security: Challenges and policies[J]. Science, 2003(302): 1917－1919.

[2]Food and Agriculture Organization (FAO). The State of Food Insecurity in the World[R]. Rome: Food and Agriculture Organization, 2012.

[3]Alexandratos N, Bruinsma J. World Agriculture Towards 2030/2050: The 2012 Revision[R]. Rome: Food and Agriculture Organization, 2012.

[4]Tilman D, Balzer C, Hill J, et al. Global food demand and the sustainable intensification of agriculture[J]. Proceedings of the National Academy of Sciences of the United States of America, 2011, 108(50): 20260–20264.

[5]IPCC. Climate Change and Land[R]. (2019-06-12) [2019-07-28] https://www.ipcc.ch/ report/srccl/.

[6]Valujeva K, O’Sullivan L, Gutzler C, et al. The challenge of managing soil functions at multiple scales: An optimisation study of the synergistic and antagonistic trade-offs between soil functions in Ireland[J]. Land Use Policy, 2016,58: 335－347.

[7]Coyle C, Creamer R E, Schulte R P O, et al. A functional land management conceptual framework under soil drainage and land use scenarios[J]. Environmental Science & Policy, 2016, 56: 39－48.

[8]Brown L. Who will Feed China? Wake-up Call for a Small Planet[M].New York: WW Norton & Company, 1995.

[9]Un, United Nations. The Millennium Development Goals Report 2015 [J]. Agenda, 2015, 23(81): 155-157.

[10]Deng X Z, Huang J K, Rozelle S, et al. Impact of urbanization on cultivated land changes in China[J]. Land Use Policy, 2015(5): 1－7.

[11]劉霈珈，吳克寧，羅明，等. 農(nóng)用地土壤重金屬超標(biāo)評(píng)價(jià)與安全利用分區(qū)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2016，32(23)：254－262. Liu Peijia, Wu Kening, Luo Ming, et al. Evaluation of agricultural land soil heavy metal elements exceed standards and safe utilization zones[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(23): 254－262. (in Chinese with English abstract)

[12]范昊明，蔡強(qiáng)國(guó)，陳光，等. 世界三大黑土區(qū)水土流失與防治比較分析[J]. 自然資源學(xué)報(bào)，2005, 20(3)：387－393. Fan Haoming, Cai Qiangguo, Chen Guang, et al. Comparative study of the soil erosion and control in the three major black soil regions in the world[J]. Journal of Natural Resources, 2005, 20(3): 387－393. (in Chinese with English abstract)

[13]趙業(yè)婷，常慶瑞，李志鵬，等. 1983—2009年西安市郊區(qū)耕地土壤有機(jī)質(zhì)空間特征與變化[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2013，29(2)：254－262. Zhao Yeting, Chang Qingrui, Li Zhipeng, et al. Spatial characteristics and changes of soil organic matter for cultivated land in suburban area of Xi’an from 1983 to 2009[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(2): 254－262. (in Chinese with English abstract)

[14]Leegood R C, Evans J R, Furbank R T. Food security requires genetic advances to increase farm yield[J]. Nature, 2010, 464(7290): 831.

[15]LeoPold A. Wilderness as land laboratory[J]. Living Wldemess, 1941, 6(2): 3.

[16]Rapport D J, McCullum J, Miller M H. Soil health: Its relationship to ecosystem health[C]. //Biological Indicators of Soil Health. Wallingford: CAB International, 1997: 29－47.

[17]Sparling G P. Soil microbial biomass, activity and nutrient cycling as indicators of soil health[C]//Biological Indicators of Soil Health. Wallingford: CAB International, 1997: 97－119.

[18]Hodgson J. Grazing Management: Science into Practice[M]. New York: Longman Scientific & Technical, 1990.

[19]Kumar S, Garg A K, Aulakh M S. Effect of conservation agriculture practices on physical, chemical and biological attributes of soil health under soybean–rapeseed rotation[J]. Agricultural Research, 2016, 5(2): 145－161.

[20]Rapport D J. What constitutes ecosystem health?[J]. Perspectives in Biology & Medicine, 1989, 33(1): 120－132.

[21]Rapport D J, Buckingham D, Cairns J, et al. Ecosystem health: The concept, the ISEH and the important tasks ahead[J]. Ecosystem Health, 1999, 5(2): 82－90.

[22]Costanza R, Norton B G, Haskell B J. Ecosystem Health: New Goals for Environmental Management[M]. Washington DC: Island Press, 1992.

[23]伍育鵬，鄖文聚，鄒如. 耕地產(chǎn)能核算模型的研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24(2)：108－113. Wu Yupeng, Yun Wenju, Zou Ru. Model for calculation of cultivated land arable[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(2): 108－113. (in Chinese with English abstract)

[24]伍育鵬，鄖文聚，鄒如. 基于耕地質(zhì)量評(píng)價(jià)的耕地產(chǎn)能核算方法及分布[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24(1)：85－89. Wu Yupeng, Yun Wenju, Zou Ru. Productivity calculation methods and distribution regulation based on cultivated land quality appraisal[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(1): 85－89. (in Chinese with English abstract)

[25]王洪波，鄖文聚，吳次芳，等. 基于農(nóng)用地分等的耕地產(chǎn)能監(jiān)測(cè)體系研究[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2008，24(4)：122－126. Wang Hongbo, Yun Wenju, Wu Cifang, et al. Farmland productivity monitoring system based on farmland gradation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2008, 24(4): 122－126. (in Chinese with English abstract)

[26]張紅富，周生路，吳紹華，等. 省域尺度耕地產(chǎn)能空間分異規(guī)律及其影響因子[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2010，26(8)：308－314. Zhang Hongfu, Zhou Shenglu, Wu Shaohua, et al. Provincial scale spatial variation of cultivated land production capacity and its impact factors[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2010, 26(8): 308－314. (in Chinese with English abstract)

[27]葉興慶. “健康產(chǎn)能”的增長(zhǎng)如何跑贏“有毒產(chǎn)能”的退出[J]. 中國(guó)發(fā)展觀察，2014(2)：10－11.

[28]李強(qiáng)，嚴(yán)金明，趙燁. 耕地健康也需診斷[J]. 中國(guó)土地，2011(11)：37－39.

[29]李強(qiáng)，彭文英，王建強(qiáng)，等. 鄉(xiāng)鎮(zhèn)企業(yè)發(fā)達(dá)區(qū)耕地健康評(píng)價(jià)與驅(qū)動(dòng)機(jī)理研究[J]. 自然資源學(xué)報(bào)，2015，30(9)：1499－1510. Li Qiang, Peng Wenying, Wang Jianqiang, et al. Health assessment and driving mechanism analysis of cultivated land in the township enterprises developed region[J]. Journal of Natural Resources, 2015, 30(9): 1499－1510. (in Chinese with English abstract)

[30]鄖文聚，梁夢(mèng)茵，湯懷志. 提升耕地質(zhì)量重在健康產(chǎn)能建設(shè)[J]. 中國(guó)土地，2015(3)：22－23.

[31]Klingebiel A A, Montgomery P H. Land Capability Classification[M]. US Department of Agriculture Handbook. 1961: 210.

[32]Food and Agriculture Organization (FAO). A Framework for Land Evaluation[R]. FAO Soil Bulletin, 32, Rome. 1976.

[33]Food and Agriculture Organization (FAO). Report on the Agro-Ecological Zones Project[R]. Rome: World Soil Resources Reports, 1978.

[34]Food and Agriculture Organization (FAO). Guidelines: Land Evaluation for Rainfed Agriculture[R]. Rome: FAO Soil Bulletin, 1983.

[35]Food and Agriculture Organization (FAO). Guidelines: Land Evaluation for irrigated Agriculture[R]. Rome: FAO Soil Bulletin, 1985.

[36]傅伯杰，陳利頂，馬誠(chéng). 土地可持續(xù)利用評(píng)價(jià)的指標(biāo)體系與方法[J]. 自然資源學(xué)報(bào)，1997，12(2)：113－118. Fu Bojie, Chen Liding, Ma Cheng. The index system and method of land sustainable use evaluation[J]. Journal of Natural Resources, 1997, 12(2): 113－118. (in Chinese with English abstract)

[37]傅伯杰，劉世梁，馬克明. 生態(tài)系統(tǒng)綜合評(píng)價(jià)的內(nèi)容與方法[J]. 生態(tài)學(xué)報(bào)，2001，21(11)：1885－1892. Fu Bojie, Liu Shiliang, Ma Keming. The contents and methods of integrated ecosystem assessment[J]. Acta Ecologica Sinica, 2001, 21(11): 1885－1892. (in Chinese with English abstract)

[38]Acton D F. A program to assess and monitor soil quality in Canada[J]. Research Branch, Agriculture Canada, 1994: 5－8.

[39]Food and Agriculture Organization (FAO). FESLM: An international Framework for evaluating sustainable land management[R]. Rome: World Soil Resources Report No. 73, 1993.

[40]盛豐. 康奈爾土壤健康評(píng)價(jià)系統(tǒng)及其應(yīng)用[J]. 土壤通報(bào)，2014，45(6)：1289－1296. Sheng Feng. Introduction and application of cornell soil health assessment[J]. Chinese Journal of Soil Science, 2014, 45(6): 1289－1296.

[41]Vincent D P O, Lal R. A standardized soil quality index for diverse field conditions[J]. Science of the Total Environment, 2016, 541: 424－434.

[42]Doran J W, Zeiss M R. Soil health and sustainability: managing the biotic component of soil quality[J]. Applied Soil Ecology, 2000(15): 3－11.

[43]The Nature Conservancy. Rethink Soil: A Roadmap to U. S. Soil Health[EB/OL]. (2017-03-23)[2019-07-28] www.nature. org/soil.

[44]耕地地力調(diào)查與質(zhì)量評(píng)價(jià)技術(shù)規(guī)程：NY/T1634—2008[S]. 北京：中華人民共和國(guó)農(nóng)業(yè)部，2008.

[45]農(nóng)用地質(zhì)量分等規(guī)程：GBT28407—2012[S]. 北京：中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì)，2012.

[46]土地質(zhì)量地球化學(xué)評(píng)價(jià)規(guī)范：DZ/T 0295—2016[S]. 北京：中國(guó)地質(zhì)調(diào)查局，2016.

[47]安萍莉，張鳳榮，陳阜. 農(nóng)用地分等定級(jí)中標(biāo)準(zhǔn)耕作制度的確定[J]. 地理與地理信息科學(xué)，2002，18(2)：45－48. An Pingli, Zhang Fengrong, Chen Fu. Determination of the standard farming system for farmland grading[J]. Geography and Territorial Research, 2002, 18(2): 45－48 (in Chinese with English abstract).

[48]溫良友，孔祥斌，張蚌蚌，等. 基于可持續(xù)發(fā)展需求的耕地質(zhì)量評(píng)價(jià)體系構(gòu)建與應(yīng)用[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(10)：234－242Wen Liangyou, Kong Xiangbin, Zhang Bangbang, et al. Construction and application of arable land quality evaluation system based on sustainable development demand[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(10): 234－242. (in Chinese with English abstract)

[49]鄖文聚，沈重陽(yáng). 中國(guó)農(nóng)用地質(zhì)量發(fā)展研究報(bào)告（2018）[M]. 北京: 中國(guó)農(nóng)業(yè)大學(xué)出版社，2018：19－32.

[50]劉焱序，彭建，汪安，等. 生態(tài)系統(tǒng)健康研究進(jìn)展[J]. 生態(tài)學(xué)報(bào)，2015，35(18)：5920－5930. Liu Yanxu, Peng Jian, Wang An, et al. New research progress and trends in ecosystem health[J]. Acta Ecologica Sinica, 2015, 35(18): 5920－5930. (in Chinese with English abstract)

[51]Yao X C, Mokbel M F, Ye S J, et al. Spatial coding-based approach for partitioning big spatial data in Hadoop[J]. Computers & Geosciences, 2017, 106: 60－67.

[52]Ye S J, Zhu D H, Yao X C, et al. Development of a highly flexible mobile GIS-based system for collecting arable land quality data[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2014, 5(14): 4432－4441.

[53]土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)：GB15618—2008[S]. 北京：中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì)，2008.

[54]農(nóng)田灌溉水質(zhì)標(biāo)準(zhǔn)：GB5084—2005[S]. 北京：中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì)，2005.

[55]鄖文聚，張蕾娜. 如何加強(qiáng)耕地質(zhì)量保護(hù)與建設(shè)管理[J]. 中國(guó)土地，2012(1)：33－34.

[56]張晉科，張鳳榮，張琳，等. 中國(guó)耕地的糧食生產(chǎn)能力與糧食產(chǎn)量對(duì)比研究[J]. 中國(guó)農(nóng)業(yè)科學(xué)，2006，39(11)：2278－2285. Zhang Jinke, Zhang Fengrong, Zhang Lin, et al. Comparison between the potential grain productivity and the actual grain yield of cultivated lands in mainland China[J]. Scientia Agricultura Sinica, 2006, 39(11): 2278－2285. (in Chinese with English abstract)

[57]Albersen P J, Fischer G, Keyzer M A, et al. Estimation of Agricultural Production Relations in the LUC Model for China[M]. Laxenburg, Austria: International Institute for Applied Systems Analysis, 2002.

[58]Fischer G, Van Velthuizen H, Shah M, et al. Global Agro-ecological Assessment for Agriculture in the 21st Century: Methodology and Results[M]. Laxenburg, Austria: International Institute for Applied Systems Analysis, 2002.

[59]Jiang B. Head/Tail Breaks: A New Classification Scheme for Data with a Heavy-Tailed Distribution[J]. The Professional Geographer, 2012, 65(3): 482－494.

[60]Edwards C A, Lofty J R. Biology of Earthworms[M] Second ed. London: Chapman and Hall, 1977: 309.

[61]Marinissen J C Y. Earthworm populations and stability of soil-structure in a silt loam soil of a reclaimed polder in the Netherlands[J]. Agriculture Ecosystems & Environment, 1994, 51: 75－87.

[62]Clements R O, Murray P J. Sturdy the impact of 20 years’ absence of earthworms and three levels of N fertilizer on a grassland soil environment[J]. Agriculture Ecosystems & Environment, 1991, 36: 75－85.

[63]Curry J P. Some effects of animal manures on earthworms in grassland[J]. Pedobiologia, 1976, 16: 425－438.

[64]Lofs-Holmin A. Influence of agricultural practices on earthworms ()[J]. Acta Agriculturae Scandinavica, 1983, 33: 225－234.

[65]Hansen S, Engelstad F. Earthworm populations in a cool and wet district as affected by tractor traffic and fertilisation Appl[J]. Applied Soil Ecology, 1999, 13: 237－250.

[66]Edwards C A, Bohlen P J. The effects of toxic chemicals on earthworms[J]. Reviews of Environmental Contamination and Toxicology, 1992, 125: 23－99.

[67]Ma W C, Brussaard L, Ridder de J A. Long-term effects of nitrogenous fertilizers on grassland earthworms (): their relation to soil acidification[J]. Agriculture Ecosystems & Environment, 1990(30): 71－80.

[68]Davis BNK. Earthworms their Ecology and Relationships with Soils and Land Use[J]. Environmental Pollution, 1985, 42(1): 94.

[69]Hugh R, Reidun P, Ragnar E, et al. Soil structure, organic matter and earthworm activity in a comparison of cropping systems with contrasting tillage, rotations, fertilizer levels and manure use[J]. Agriculture Ecosystems & Environment, 2008, 124(3): 275－284.

[70]Stenberg Bo. Monitoring soil quality of arable land: Microbiological indicators[J]. Acta Agriculturae Scandinavica, Section B-Soil & Plant Science, 1999, 49(1): 1－24.

[71]Lynch J M. Interactions between biological processes, cultivation and soil structure[J]. Plant and Soil, 1984, 76(1): 307－318.

[72]Lynch J M, Bragg E. Microorganisms and Soil Aggregate Stability//Stewart B A(eds). Advances in Soil Science. vol 2.[M]. New York: Springer, 1985.

[73]Gupta V V S R, Germida J J. Distribution of microbial biomass and its activity in different soil aggregate size classes as affected by cultivation[J]. Soil Biology and Biochemistry, 1988, 20(6): 777－786.

[74]Arshad M A, Martin S. Identifying critical limits for soil quality indicators in agro-ecosystems[J]. Agriculture Ecosystems and Environment, 2002, 88(2): 153－160.

[75]Vincent D P O, Lal R. Using meta-analyses to assess pedo-variability under different land uses and soil management in central Ohio, USA[J]. Geoderma, 2014, (232/233/234): 56－68.

[76]Vincent D P O, Lal R. A standardized soil quality index for diverse field conditions[J]. Science of The Total Environment, 2016, 541: 424－434.

[77]Mcbratney A, Field D J, Koch A. The dimensions of soil security[J]. Geoderma, 2014, 213: 203－213.

[78]宋長(zhǎng)青，程昌秀，史培軍. 新時(shí)代地理復(fù)雜性的內(nèi)涵[J]. 地理學(xué)報(bào)，2018，73(7)：1204－1213. Song Changqing, Cheng Changxiu, Shi Peijun. Geography complexity: New connotations of geography in the new era[J]. Acta Geographica Sinica, 2018, 73(7): 1204－1213. (in Chinese with English abstract)

[79]程昌秀，史培軍，宋長(zhǎng)青，高劍波. 地理大數(shù)據(jù)為地理復(fù)雜性研究提供新機(jī)遇[J]. 地理學(xué)報(bào)，2018，73(8)：1397－1406. Cheng Changxiu, Shi Peijun, Song Changqing, Gao Jianbo. Geographic big-data: A new opportunity for geography complexity study[J]. Acta Geographica Sinica, 2018, 73(8): 1397－1406. (in Chinese with English abstract)

[80]Drobnik T, Greiner L, Keller A, et al. Soil quality indicators-From soil functions to ecosystem services[J]. Ecological Indicators, 2018, 94(1): 151－169.

[81]葉思菁. 大數(shù)據(jù)環(huán)境下遙感圖譜應(yīng)用方法研究—以作物干旱監(jiān)測(cè)為例[J]. 測(cè)繪學(xué)報(bào)，2018，47(6)：892. Ye Sijing. Research on application of remote sensing Tupu: Take monitoring of meteorological disaster for example[J]. Acta Geodaetica et Cartographica Sinica, 2018, 47(6): 892. (in Chinese with English abstract)

[82]葉思菁，張超，王媛，等. GF-1遙感大數(shù)據(jù)自動(dòng)化正射校正系統(tǒng)設(shè)計(jì)與實(shí)現(xiàn)[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2017(33)：266－273. Ye Sijing, Zhang Chao, Wang Yuan, et al. Design and implementation of automatic orthorectification system based on GF-1 big data[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2017(33): 266－273. (in Chinese with English abstract)

[83]Ye Sijing, Yan Tailai, Yue Yanli et al. Developing a reversible rapid coordinate transformation model for the cylindrical projection[J]. Computers & Geosciences, 2016, 89(4): 44－56.

[84]Ye Sijing, Cheng Changxiu, Song Changqing, et al. RDCRMG: A raster dataset clean & reconstitution multi-grid architecture for remote sensing monitoring of vegetation dryness[J]. Remote Sensing, 2018(10): 1－24.

[85]Juhos K, Czigány S, Madarász B, et al. Interpretation of soil quality indicators for land suitability assessment-A multivariate approach for Central European arable soils[J]. Ecological Indicators, 2019(99): 261－272.

[86]Clivot H, Mouny J, Duparque A, et al. Modeling soil organic carbon evolution in long-term arable experiments with AMG model[J]. Environmental Modelling & Software, 2019, 118: 99－113.

Cultivated land health-productivity comprehensive evaluation and its pilot evaluation in China

Ye Sijing1,2, Song Changqing1,2, Cheng Feng3, Zhang Leina3, Cheng Changxiu1,2※, Zhang Chao4, Yang Jianyu4, Zhu Dehai4

(1.100875;2.100875; 3.100035; 4.100083)

Evaluation of the healthy productivity of cultivated land in China, analysis of key problems of different regions that endangers cultivated land health is of great significance to explore sustainable cultivated land utilization path and safeguard national food security and ecological security. The evaluation of cultivated land healthy productivity should be executed base on integration of existing evaluation indicators system and evaluation method achievements and clarification of key factor states and driving processes that affecting cultivated land health from the perspective of cultivated land components. And interaction characteristics between cultivated land health and cultivated land productivity should be analyzed by coupling them with each other. In this paper, we selected 65 counties (cities) nationwide as research regions with considering the difference of climate conditions in different regions of China, and used 5 categories of indicators including topographic characteristics, soil properties, tillage conditions, environmental conditions and biological characteristics to present the cultivated land health status. And on that basis, we defined critical health threshold and critical sub-health threshold of different indicators by analyzing counties’ indicator standard deviation convergence characteristics and scoring rules. We divided all pilot counties into three statuses: Health, sub-health and unhealth, and analyzed spatial distribution features of counties with different health status. Furthermore, regional climate conditions indicator, cultivated land natural conditions indicator and tillage technique level indicator had been synthetically integrated to evaluate the productivity of cultivated land, and regional cultivated land health-productivity interaction characteristic had been analyzed. Experimental result showed that counties whose cultivated land in sub-healthy or unhealthy status expressed significant regional features, and the key factors of limiting cultivated land health were different among pilot counties. On the one hand, the sub-health status cases caused by inferior farming conditions indicator were more common (accounting for nearly 40% of pilot counties) than other indicators, which was not only manifested by serious fragmentation of cultivated land and imperfect farmland infrastructure, but also reflected the urgent requirement for promoting high-standard farmland construction. On the other hand, the form of biological characteristics was severe, indicating that under the interaction of climatic conditions, soil structure and tillage technology, the ability of cultivated land to guarantee the healthy growth of animals was facing degradation. Furthermore, through coupling analysis of cultivated land health and cultivated land productivity, we divided the cultivated land healthy productivity characteristics of pilot counties into six categories: “health-high productivity”, “health-low productivity”, “sub-health-high productivity”, “sub-health-low productivity”, “unhealth-high productivity” and “unhealth-low productivity”, and respectively analyzed the regional farmland protection directions for each category. “Sub-health/unhealth-high productivity” type was widely distributed, which was the core region for implementation of “storing grain on land and storing grain in technology” strategy in China. Therefore, it is necessary to give consideration to the health and productivity of cultivated land in this region, and strengthen the research and practice on methods to improve the health of cultivated land. This study provides scientific basis for improving cultivated land healthy productivity and provides effective support for designing protection paths of cultivated land resources in different regions.

land use; GIS; cultivated land; cultivated land health; cultivated land productivity; climate change

葉思菁，宋長(zhǎng)青，程 鋒，張蕾娜，程昌秀，張 超，楊建宇，朱德海. 中國(guó)耕地健康產(chǎn)能綜合評(píng)價(jià)與試點(diǎn)評(píng)估研究[J].農(nóng)業(yè)工程學(xué)報(bào)，2019，35(22)：66－78.doi：10.11975/j.issn.1002-6819.2019.22.008 http://www.tcsae.org

Ye Sijing, Song Changqing, Cheng Feng, Zhang Leina, Cheng Changxiu, Zhang Chao, Yang Jianyu, Zhu Dehai. Cultivated land health-productivity comprehensive evaluation and its pilot evaluation in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(22): 66－78. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2019.22.008 http://www.tcsae.org

2019-08-12

2019-09-26

“美麗中國(guó)”科學(xué)院先導(dǎo)科技專項(xiàng)（XDA23100303）；國(guó)家自然科學(xué)基金資助項(xiàng)目（41801300）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金（2018NTST03）。

葉思菁，博士，講師，主要從事耕地安全評(píng)價(jià)與耕地可持續(xù)利用路徑研究。Email：yesj@bnu.edu.cn

程昌秀，博士，博士生導(dǎo)師，主要從事空間信息技術(shù)及其在土地資源領(lǐng)域應(yīng)用研究。Email：chengcx@bnu.edu.cn

10.11975/j.issn.1002-6819.2019.22.008

F323.21

A

1002-6819(2019)-22-0066-13