李玉恒，黃惠倩，郭桐冰，杜國明

（1.中國科學院地理科學與資源研究所區(qū)域可持續(xù)發(fā)展分析與模擬重點實驗室，北京 100101；2.中國科學院大學資源與環(huán)境學院，北京100049；3.東北農(nóng)業(yè)大學公共管理與法學院，黑龍江 哈爾濱 150030）

1 引言

“零饑餓”是聯(lián)合國2030年可持續(xù)發(fā)展議程17項目標之一，聚焦于在全球范圍內(nèi)消除一切形式的饑餓和營養(yǎng)不良，保證所有人——尤其是兒童及弱勢群體享有充足且營養(yǎng)的食物。糧食安全關(guān)系到人類永續(xù)發(fā)展和前途命運，是保障世界和平與發(fā)展、構(gòu)建人類命運共同體的重要基礎(chǔ)。作為糧食生產(chǎn)的載體，耕地保護與可持續(xù)利用已成為保障糧食安全、支撐經(jīng)濟社會平穩(wěn)健康發(fā)展的關(guān)鍵途徑。作為世界最大的發(fā)展中國家，我國人多地少、人均耕地資源稀缺的基本國情要求實行最嚴格的耕地保護制度。自1981年國務(wù)院明確提出耕地保護的概念以來，我國構(gòu)建了耕地資源數(shù)量、質(zhì)量、生態(tài)“三位一體”的保護格局，穩(wěn)步提升了糧食生產(chǎn)能力，實現(xiàn)了連續(xù)多年的豐產(chǎn)豐收，支撐了工業(yè)化、城鎮(zhèn)化及農(nóng)業(yè)現(xiàn)代化發(fā)展。當前，我國面臨著耕地“非糧化”與“非食物化”、邊際化、細碎化、質(zhì)量退化、自然災(zāi)害破壞等諸多挑戰(zhàn)[1-3]。

黑土是指擁有黑色或暗黑色腐殖質(zhì)表土層的土壤，是地球上珍貴的土壤資源。黑土耕地是重要的農(nóng)業(yè)生產(chǎn)要素，是一種性狀好、肥力高、適宜農(nóng)耕的優(yōu)質(zhì)土地①《國家黑土地保護工程實施方案（2021—2025年）》。，被稱作“耕地中的大熊貓”。東北地區(qū)是世界主要黑土帶之一，也是我國糧食主產(chǎn)區(qū)和重要的商品糧基地，糧食產(chǎn)量占全國的1/4，商品量占全國的1/4，調(diào)出量占全國的1/3，是國家糧食安全的“壓艙石”[4]。東北黑土區(qū)范圍涉及遼寧、吉林、黑龍江以及內(nèi)蒙古東部的部分地區(qū)，總面積55.6萬km2。然而，長期大規(guī)模、高強度的開墾導致東北黑土地土壤肥力普遍下降，黑土地資源量在減少、質(zhì)在退化的趨勢愈發(fā)明顯，威脅國家糧食安全與區(qū)域生態(tài)安全，亟需科學推進黑土耕地的保護與有序利用[5]。

當前，我國啟動了東北黑土地保護利用工程，制定了《東北黑土地保護性耕作行動計劃（2020—2025年）》 《國家黑土地保護工程實施方案（2021—2025年）》。耕地保護是一項系統(tǒng)工程，耕地數(shù)量、質(zhì)量與生態(tài)面臨著源于自身屬性和外部環(huán)境變化的多重壓力。已有研究顯示，耕地屬性、土壤類型、農(nóng)戶耕作行為等內(nèi)在因素與土地制度與政策、建設(shè)用地擴張、農(nóng)業(yè)污染、氣候變化等外界因素相互交織，共同影響著耕地保護方式與利用成效[6-13]?？茖W構(gòu)建耕地保護體系，提升耕地應(yīng)對外界擾動沖擊和不確定風險的韌性，對于保障耕地數(shù)量、品質(zhì)與糧食生產(chǎn)能力具有重要的現(xiàn)實意義。本文以黑龍江省拜泉縣為例，將韌性理論引入到耕地保護研究，基于村級調(diào)查數(shù)據(jù)，綜合評價了多重壓力脅迫下東北黑土區(qū)耕地韌性，揭示了耕地韌性空間分布特征、關(guān)鍵短板因素，提出了旨在提升黑土區(qū)耕地韌性的路徑措施與政策建議。

2 耕地韌性研究框架

韌性是系統(tǒng)的基本屬性，由于未來發(fā)展的不確定性和不可預(yù)測性，抵御外界風險挑戰(zhàn)沖擊的韌性對于維持系統(tǒng)穩(wěn)定性具有至關(guān)重要的作用。1973年，生態(tài)學家HOLLING將韌性這一概念率先應(yīng)用于生態(tài)學領(lǐng)域，指生態(tài)系統(tǒng)受到外界壓力沖擊后，最大程度吸收擾動沖擊并恢復(fù)到新的均衡狀態(tài)的能力[14]。隨著學科的發(fā)展，韌性研究逐漸應(yīng)用到經(jīng)濟學、社會學、災(zāi)害學等領(lǐng)域，成為解析社會—生態(tài)系統(tǒng)可持續(xù)性的重要工具之一[15-18]。

韌性理論指出系統(tǒng)內(nèi)核要素的類型、數(shù)量、質(zhì)量及其耦合方式構(gòu)成了韌性基礎(chǔ)，在與外部環(huán)境交互作用過程中，系統(tǒng)展現(xiàn)出對壓力沖擊的抵御能力、適應(yīng)能力與轉(zhuǎn)型能力[19]。耕地系統(tǒng)是由相互影響、相互作用的若干要素按照一定結(jié)構(gòu)組成的、具有特定功能的有機整體。在與外界環(huán)境系統(tǒng)交互作用過程中，耕地具有通過自身調(diào)節(jié)，抵御擾動沖擊而保持系統(tǒng)穩(wěn)定、功能不變的韌性，而耕地非農(nóng)化、土質(zhì)退化與污損化則是擾動沖擊作用于耕地系統(tǒng)，超出了韌性調(diào)控閾值而產(chǎn)生的直接后果。

韌性是實現(xiàn)耕地系統(tǒng)功能穩(wěn)定、資源可持續(xù)利用的重要支撐。在面對不同壓力沖擊時，耕地系統(tǒng)的響應(yīng)方式與程度不盡相同，從而展現(xiàn)出差異的韌性。相較于耕地質(zhì)量評價、耕地生態(tài)與資源安全，耕地韌性體現(xiàn)了耕地系統(tǒng)對外界環(huán)境變化沖擊與風險挑戰(zhàn)的應(yīng)對能力。當前，韌性研究多集中于經(jīng)濟、生態(tài)、城鄉(xiāng)、災(zāi)害等學科領(lǐng)域，涉及耕地韌性的研究相對較少。孟麗君等[20]基于表象韌性和潛在韌性，從穩(wěn)定增長能力、高效利用能力、抵抗干擾能力、資源支撐能力4個方面建立了指標體系，評價了河北曲周縣1973—2014年耕地利用系統(tǒng)韌性。李窚琪等[21]從耕地系統(tǒng)的功能性、穩(wěn)定性、冗余度和適應(yīng)性4個維度構(gòu)建了指標體系，評價了昆明市2001—2017年耕地系統(tǒng)韌性。已有研究主要依賴于統(tǒng)計數(shù)據(jù)，且缺乏對耕地面臨壓力的判斷，不利于揭示不同壓力下的耕地韌性。此外，村域是耕地保護與農(nóng)業(yè)生產(chǎn)的基礎(chǔ)單元，以村域為單元開展耕地韌性研究，有利于揭示不同村域耕地立地條件、農(nóng)田基礎(chǔ)設(shè)施建設(shè)水平、農(nóng)戶生計方式等因素對耕地韌性的影響。

基于此，本文從韌性發(fā)生學視角，構(gòu)建了基于“壓力—狀態(tài)—響應(yīng)”的耕地韌性解釋模型。如圖1所示，耕地內(nèi)核系統(tǒng)包括地理區(qū)位、立地條件、土壤類型、土壤養(yǎng)分等因素，形成了耕地的“狀態(tài)”，包括耕地適宜性、耕地質(zhì)量、耕作方式、生產(chǎn)力等現(xiàn)狀特征。制度設(shè)計與宏觀政策、非農(nóng)用地需求與污染、農(nóng)戶生計決策、氣候變化與災(zāi)害等因素交織，構(gòu)成了對耕地系統(tǒng)的“壓力”?！绊憫?yīng)”指在壓力作用時及作用后，為維系耕地功能不變、系統(tǒng)穩(wěn)定所采取的應(yīng)對措施如土地整治、土壤修復(fù)、耕地流轉(zhuǎn)、農(nóng)技投入等。因而，耕地韌性的發(fā)生過程可分解為：擾動源，即外界風險挑戰(zhàn)作用于耕地系統(tǒng)并產(chǎn)生相關(guān)影響；擾動中，即耕地系統(tǒng)抵御壓力沖擊的能力及在壓力影響下的系統(tǒng)狀態(tài)；擾動后，即耕地系統(tǒng)從擾動沖擊中恢復(fù)正常狀態(tài)的舉措，以及在沖擊后尋求發(fā)展轉(zhuǎn)型的能力。

圖1 基于“壓力—狀態(tài)—響應(yīng)”的耕地韌性模型Fig.1 Cultivated land resilience model based on “pressure-state-response”

3 研究方法與數(shù)據(jù)

3.1 拜泉縣概況

拜泉縣隸屬黑龍江省齊齊哈爾市，介于北緯47°17′～47°54′、東經(jīng)125°29′～126°31′之間。該縣地處黑龍江省中西部、小興安嶺與松嫩平原的交匯地帶，地勢東高西低，具有典型漫川漫崗的地貌特征，海拔最高處365 m，最低點92 m，境內(nèi)主要河流為雙陽河、通肯河。拜泉縣轄7鎮(zhèn)9鄉(xiāng)，186個行政村，總?cè)丝?5萬人，其中農(nóng)村常住人口17萬人①www.baiquan.gov.cn.，國土總面積3 599.15 km2，其中耕地面積2 057.13 km2②拜泉縣第三次國土調(diào)查。。作為東北黑土區(qū)的典型縣域，拜泉縣土壤類型以黑土為主，黑鈣土、草甸土為輔，土層厚度1 m以上。氣候?qū)僦袦貛Т箨懶约撅L氣候，雨熱同期，年均降雨量490 mm，年均溫度1.2℃，年均日照2 730 h，無霜期122天①。

3.2 數(shù)據(jù)來源及處理方式

本文所用到的經(jīng)濟、人口、農(nóng)業(yè)數(shù)據(jù)均來源于2020年筆者及團隊對拜泉縣186個行政村的實地調(diào)研與問卷調(diào)查。經(jīng)核查有32份問卷由于相關(guān)數(shù)據(jù)失實，判為無效問卷，最終獲得154份有效問卷。坡度、高程等地形數(shù)據(jù)來源于地理空間數(shù)據(jù)云提供的30 m分辨率數(shù)字高程模型（DEM）。村域?qū)用娴母孛娣e、有排澇設(shè)施農(nóng)田比例等地理空間數(shù)據(jù)來自于拜泉縣第三次國土調(diào)查和村莊地籍數(shù)據(jù)。

在村莊人口問卷中，存在村內(nèi)純務(wù)農(nóng)人數(shù)大于村內(nèi)勞動力人數(shù)的情況，根據(jù)研究區(qū)60周歲以上老人務(wù)農(nóng)較為普遍的現(xiàn)實情況，在純務(wù)農(nóng)人數(shù)的指標統(tǒng)計中沒有統(tǒng)計有勞動力人口中的純務(wù)農(nóng)人口，而是實際生活中仍然從事農(nóng)業(yè)勞作的人口。在農(nóng)業(yè)生產(chǎn)問卷中，由于各村種植作物情況不同，在計算種植每畝耕地的投入成本及糧食產(chǎn)量指標時，先根據(jù)每村種植各類作物面積的占比賦權(quán)重，再基于所得權(quán)重計算得出每畝耕地的實際情況；有排澇設(shè)施農(nóng)田比例是根據(jù)第三次國土調(diào)查數(shù)據(jù)中提取的各村水田、水澆地面積占總耕地面積的比值計算得出。結(jié)合本文研究目的，耕地平均坡度僅提取集體所有耕地，對于權(quán)屬類型為國有的耕地進行了裁剪處理，再通過ArcGIS 10.6的柵格處理、區(qū)域分析、表面分析等工具從DEM柵格數(shù)據(jù)中提取。

在這個科技發(fā)達的時代，許多新的設(shè)施已經(jīng)投入到電網(wǎng)的使用中去。大城市的整個電力系統(tǒng)已經(jīng)煥然一新，這些設(shè)施的使用也很好的經(jīng)濟的控制了供電方面的浪費。所以，盡最大可能改善農(nóng)村配網(wǎng)的設(shè)施可以有效發(fā)展農(nóng)村電力問題，還可以為農(nóng)村節(jié)省一筆不必要的浪費。

3.3 研究方法

（1）評價指標體系。本文從韌性發(fā)生學視角構(gòu)建了基于“壓力—狀態(tài)—響應(yīng)”的耕地韌性評價模型，以村域為單元對拜泉縣耕地韌性進行評價。耕地保護聚焦于耕地數(shù)量、質(zhì)量和生態(tài)三位一體，具有層次性和系統(tǒng)性。耕地系統(tǒng)穩(wěn)定性與可持續(xù)性面臨著多方面的壓力。首先，地形條件通過對水熱條件的再分布影響耕地生產(chǎn)能力，并間接影響著耕地建設(shè)、開發(fā)和改良等的經(jīng)濟效益[23]。東北黑土區(qū)坡耕地面積較大，以漫崗長坡耕地為主，匯水面積大，易形成水蝕，直接影響到耕地質(zhì)量與生產(chǎn)力。其次，隨著農(nóng)戶生計方式的非農(nóng)化轉(zhuǎn)型，農(nóng)戶不再以地為生的特征凸顯，導致耕地投入減少、土地粗放及低效利用等問題出現(xiàn)，對耕地的有效利用與穩(wěn)定產(chǎn)出帶來巨大挑戰(zhàn)。

根據(jù)耕地內(nèi)核系統(tǒng)構(gòu)成要素及拜泉縣耕地現(xiàn)狀特征，黑土耕地的“狀態(tài)”包括戶均地塊數(shù)量、每畝耕地糧食產(chǎn)量、每畝耕地純收入。這些狀態(tài)要素決定了在面臨立地條件、農(nóng)戶生計轉(zhuǎn)型壓力沖擊時，為實現(xiàn)耕地保護目標所采取的響應(yīng)措施，包括單位面積耕地的勞動力投入、機械投入、農(nóng)資投入及農(nóng)田減災(zāi)防災(zāi)設(shè)施建設(shè)?；谝陨戏治觯疚臉?gòu)建了黑土區(qū)耕地韌性的評價指標體系，包括3個維度9項指標（表1）。

表1 黑土區(qū)耕地韌性評價指標體系Tab.1 Evaluation index system of cultivated land resilience of the black soil region

（2）熵值法。本文采用熵值法確定耕地韌性指標體系各指標所占權(quán)重，測算多重壓力脅迫下拜泉縣的耕地韌性，具體測算方法如下：

首先，假設(shè)某一研究區(qū)域有n個村莊，h項測評指標，建立原始指標矩陣為：X= {xij}n×h（1≤i≤n，1≤j≤h），xij為第i個村莊的第j項指標值。對指標數(shù)據(jù)進行無量綱化處理，利用離差標準化法將指標體系中各項原始指標轉(zhuǎn)換為無量綱化指標。本文中涉及壓力指標以及耕地集中連片程度為負向指標，其余數(shù)據(jù)為正向指標，標準化處理結(jié)果Zij的計算如式（1）、式（2）所示，各評估指標經(jīng)標準化處理后的值介于[ 0，1]。

其次，確定各評價指標權(quán)重。在對指標進行如式（1）與式（2）的歸一化處理基礎(chǔ)上，計算各指標熵值Ej，如式（3）所示。

最后，利用式（4）計算熵值冗余度，在式（5）的支持下得到權(quán)重Wj，接下來按式（6）進行耕地韌性測算，得到各村莊的耕地韌性水平Ri。

（3）多重壓力脅迫模型。為探究黑土區(qū)耕地在不同壓力脅迫下的韌性情況，本文分別對立地條件、農(nóng)戶生計轉(zhuǎn)型兩種壓力脅迫下的耕地韌性進行測算，并比較分析了單個壓力以及兩種壓力共同作用下的耕地韌性變化情況。

在本文中，兩種壓力指標Pij均為負向指標，首先利用離差標準化法對壓力指標進行無量綱化處理，處理結(jié)果均為正數(shù)。其次依次對不同壓力脅迫下的黑土區(qū)耕地韌性進行測算，基于式（5）得出各指標權(quán)重Wj，根據(jù)指標體系表1對指標數(shù)據(jù)進行標準化處理得到壓力指標Pij，狀態(tài)指標Sij，響應(yīng)指標RPij，利用式（7）進行m（1≤m≤3）種壓力條件下的耕地韌性測算，得到不同情形下的耕地韌性指數(shù)Lim。經(jīng)過多重壓力脅迫模型計算得出的韌性指數(shù)越小，則表示該類型壓力的負效應(yīng)越明顯。

（4）冷熱點分析技術(shù)。冷熱點分析屬于空間聚類方法，通過判斷高值或低值要素在空間上發(fā)生聚類的位置識別要素空間聚集的分布規(guī)律。為進一步揭示拜泉縣耕地韌性的空間集聚性，本文利用ArcGIS 10.6平臺的Getis-OderGi*功能對拜泉縣不同鄉(xiāng)村耕地韌性數(shù)值進行分析，探測變化量的冷熱點區(qū)域。冷熱點分析具體模型如式（8）、式（9）所示。

式（8）—式（9）中：xj是每個統(tǒng)計網(wǎng)格的屬性值；wij是i和j的空間權(quán)重，在探測中采取空間搜索半徑法，以100 m為相關(guān)距離構(gòu)建出空間權(quán)重矩陣；指空間權(quán)重平均值；S為空間權(quán)重標準差。

4 研究結(jié)果

4.1 多重壓力脅迫下的耕地韌性

測算結(jié)果顯示，拜泉縣耕地韌性為0.421，總體處于中等偏下水平，體現(xiàn)了立地條件與農(nóng)戶生計轉(zhuǎn)型雙重壓力脅迫下，拜泉縣耕地韌性不足的現(xiàn)實境況，也印證了長期高強度利用下東北黑土區(qū)耕地質(zhì)量退化的現(xiàn)狀特征，表現(xiàn)為黑土層在“變薄、變瘦、變硬”[24-25]。

表2展示了多重壓力脅迫下的拜泉縣耕地韌性水平。由模型1和模型2可知，立地條件壓力、農(nóng)戶生計轉(zhuǎn)型壓力脅迫下的耕地韌性分別為0.409和0.418，表明耕地韌性受立地條件的影響相對較大。拜泉縣地勢丘陵起伏、漫川漫崗，水土流失比較嚴重。坡耕地因坡緩、坡長，匯水面積大，容易引發(fā)水蝕，坡耕地水土流失面積占全縣耕地面積的70%，破壞了黑土土壤結(jié)構(gòu)，直接影響到耕地質(zhì)量與生產(chǎn)力。此外，坡耕地的侵蝕溝日益增多，耕地被切割得支離破碎，嚴重影響了機耕作業(yè)，勞動生產(chǎn)率下降。

表2 耕地韌性測算結(jié)果Tab.2 Calculation results of cultivated land resilience

研究進一步測算了雙重壓力組合情形下的耕地韌性。如模型3結(jié)果所示，在立地條件不變的情形下，疊加農(nóng)戶生計轉(zhuǎn)型壓力，拜泉縣耕地韌性由0.409增長到0.421，表明農(nóng)戶生計轉(zhuǎn)型的壓力脅迫并未對耕地韌性產(chǎn)生顯著負效應(yīng)。調(diào)查數(shù)據(jù)顯示，拜泉縣農(nóng)村居民年人均純收入為8 651元，低于齊齊哈爾市（16 689元）和黑龍江省平均水平（14 982元）。然而拜泉縣人均13畝耕地，畝均純收入達325元，人均務(wù)農(nóng)收入占到了農(nóng)戶人均純收入的一半，因而農(nóng)戶愛惜耕地，能夠有效管護、利用耕地，及時完成農(nóng)業(yè)生產(chǎn)各項投入，從而保證了耕地質(zhì)量與生產(chǎn)力。研究結(jié)果表明，在立地條件不變的情形下，農(nóng)戶對耕地的投入與管護有助于提升耕地質(zhì)量與生產(chǎn)力，從而提高了耕地韌性。

4.2 耕地韌性空間分布

圖2展示了拜泉縣耕地韌性的空間分布情況，耕地韌性較高的鄉(xiāng)村主要集中于拜泉縣西部的永勤鄉(xiāng)、富強鎮(zhèn)、愛農(nóng)鄉(xiāng)、豐產(chǎn)鄉(xiāng)、時中鄉(xiāng)，北部的長春鎮(zhèn)以及東部的三道鎮(zhèn)，該區(qū)域耕地坡度較低，地勢較緩，水蝕壓力較小，適宜農(nóng)作物的規(guī)?；N植，耕地韌性平均值較高，為0.442。然而，拜泉縣中部的拜泉鎮(zhèn)、興華鄉(xiāng)、興國鄉(xiāng)與新生鄉(xiāng)的耕地坡度較高，平均坡度大于10°，受水土流失影響，上述鄉(xiāng)鎮(zhèn)的耕地韌性平均值為0.342，低于全縣平均水平。

圖2 拜泉縣耕地韌性空間分布圖Fig.2 Spatial distribution of cultivated land resilience in Baiquan County

研究進一步分析了以鄉(xiāng)鎮(zhèn)為單元的耕地韌性空間分布情況。如圖3所示，拜泉鎮(zhèn)、興華鄉(xiāng)與興國鄉(xiāng)耕地韌性水平相對較低，測算結(jié)果分別0.313、0.355和0.298，均顯著低于拜泉縣耕地平均韌性水平。上述地區(qū)多位于漫川漫崗地帶，耕作條件較差，加之臨近縣城，農(nóng)戶生計非農(nóng)化特征明顯，抑制了耕地投入。調(diào)查數(shù)據(jù)顯示，上述鄉(xiāng)鎮(zhèn)的每畝機械投入與化肥投入僅為31.8元和37.9元，分別低于全縣57.7元和65.4元的平均水平。這些鄉(xiāng)鎮(zhèn)的糧食平均產(chǎn)量為177.3 kg，畝均純收入為132.1元，分別低于全縣畝均產(chǎn)出270.3 kg和畝均純收入325.2元的水平。

圖3 拜泉縣各鄉(xiāng)鎮(zhèn)耕地韌性空間分布Fig.3 Spatial distribution of cultivated land resilience at township level in Baiquan County

永勤鄉(xiāng)、富強鎮(zhèn)、三道鎮(zhèn)耕地韌性較高，測算結(jié)果分別為0.479、0.462和0.466，均高于全縣的耕地韌性水平。該地區(qū)地勢平坦，有利于機械化大規(guī)模耕作，農(nóng)戶的耕地投入較大。調(diào)查顯示上述鄉(xiāng)鎮(zhèn)每畝耕地的機械投入與化肥投入分別為60.0元和78.9元，均高于全縣平均水平。糧食畝均產(chǎn)量與畝均純收入分別為334.3 kg和405.1元，顯著高于全縣平均水平。

4.3 耕地韌性空間集聚性

為進一步揭示拜泉縣耕地韌性的空間集聚性，本文基于拜泉縣各村耕地韌性數(shù)據(jù)，運用ArcGIS平臺揭示了拜泉縣耕地韌性高值與低值的空間聚集情況，通過運用空間“冷點—熱點”模型，得到拜泉縣耕地韌性高低值的冷熱點分布情況（圖4）。研究發(fā)現(xiàn)，在90%的置信水平下，耕地韌性冷點區(qū)集中于拜泉鎮(zhèn)及其周圍的興華鄉(xiāng)與興國鄉(xiāng)，而熱點區(qū)集中于拜泉縣西部的愛農(nóng)鄉(xiāng)、永勤鄉(xiāng)、豐產(chǎn)鄉(xiāng)，北部的長春鎮(zhèn)以及東部的三道鎮(zhèn)。由于縣城所在地拜泉鎮(zhèn)及周邊地區(qū)耕地破碎化程度較高，不利于機械化規(guī)模化土地經(jīng)營，耕地生產(chǎn)效益低下，拜泉鎮(zhèn)耕地畝均純收入僅為65.3元，遠低于全縣325.2元的平均水平。縣城周邊居民多從事非農(nóng)產(chǎn)業(yè)，對耕地投入相對較少，調(diào)查數(shù)據(jù)顯示，拜泉鎮(zhèn)、興華鄉(xiāng)與興國鄉(xiāng)每畝耕地的機械投入與化肥投入平均值僅為36.5元與37.9元，顯著低于全縣57.7元與65.4元的平均水平，導致耕地產(chǎn)出降低。

圖4 拜泉縣耕地韌性冷熱點分析圖Fig.4 Analysis of cold and hot spots of cultivated land resilience in Baiquan County

耕地韌性熱點區(qū)集中于愛農(nóng)鄉(xiāng)、永勤鄉(xiāng)、長春鎮(zhèn)、三道鎮(zhèn)，上述鄉(xiāng)鎮(zhèn)距離縣城較遠，耕地廣闊平坦，易于開展規(guī)模化經(jīng)營，農(nóng)戶對耕地投入較高，每畝機械投入與化肥投入平均值達60.9元與66.1元，高于全縣平均水平，土地生產(chǎn)效益較好，上述鄉(xiāng)鎮(zhèn)畝均純收入平均值為400.7元，高于全縣325.2元的平均水平。

4.4 耕地韌性短板因素

本文采用雷達圖展示了拜泉縣耕地韌性的構(gòu)成要素。如圖5所示，拜泉縣耕地韌性的主要支撐要素為每畝耕地的糧食產(chǎn)量、化肥與機械投入成本、每畝耕地純收入以及耕地集中連片程度。調(diào)查數(shù)據(jù)顯示，拜泉縣每畝耕地糧食產(chǎn)量為270.3 kg，高于全國173.1 kg的平均水平。每畝耕地化肥投入65.4元、機械投入57.7元，均高于全市平均水平，有助于提升耕地質(zhì)量，保障糧食生產(chǎn)水平。當前，拜泉縣人均13畝耕地，每畝耕地純收入325.2元，人均務(wù)農(nóng)純收入達4 227元，占到了農(nóng)戶人均純收入（8 651元）的一半，有利于保障對耕地的有效投入。經(jīng)調(diào)查測算，拜泉縣耕地流轉(zhuǎn)率為52.9%，種糧大戶、農(nóng)業(yè)合作社、涉農(nóng)企業(yè)等成為耕地經(jīng)營的主力軍，通過土地整治，提高了耕地平整度與集中連片水平，有助于實施農(nóng)業(yè)生產(chǎn)經(jīng)營規(guī)?；⑷套鳂I(yè)機械化、社會服務(wù)專業(yè)化，顯著提升了耕地產(chǎn)出效益與抵御市場波動、自然災(zāi)害沖擊的能力。

圖5 拜泉縣耕地韌性構(gòu)成雷達圖Fig.5 Radar map of cultivated land resilience composition in Baiquan County

然而，有排澇設(shè)施的農(nóng)田比例以及每畝耕地的純務(wù)農(nóng)人數(shù)為拜泉縣耕地韌性的短板因素，揭示了立地條件與農(nóng)戶生計轉(zhuǎn)型雙重壓力脅迫下，拜泉縣農(nóng)業(yè)勞動力投入不足、耕地抗災(zāi)害能力不強的現(xiàn)狀，降低了耕地韌性。調(diào)查數(shù)據(jù)顯示，拜泉縣鄉(xiāng)村有排澇設(shè)施的農(nóng)田比例僅為4%，缺乏抵御洪澇災(zāi)害的能力，增大了農(nóng)作物受洪澇災(zāi)害影響而導致產(chǎn)量下跌的風險。此外，拜泉縣農(nóng)村人口老齡化程度較高，60歲以上老齡人口占全縣常住人口數(shù)34%，大量青壯年勞動力外出務(wù)工導致農(nóng)業(yè)勞動力以老齡化人口為主，習慣采用傳統(tǒng)耕作方式，耕地保護意識薄弱，接受先進農(nóng)業(yè)生產(chǎn)知識與黑土地保護技術(shù)的積極性與主動性普遍較差，不利于黑土地保護性耕作，弱化了耕地韌性。

5 結(jié)論與討論

黑土耕地是重要的農(nóng)業(yè)生產(chǎn)要素，加強黑土地保護，關(guān)系到我國糧食安全、土地資源永續(xù)利用和農(nóng)業(yè)可持續(xù)發(fā)展，具有極其重要的戰(zhàn)略意義。本文以黑龍江省拜泉縣為例，基于覆蓋全縣16個鄉(xiāng)鎮(zhèn)、154個行政村的調(diào)查數(shù)據(jù)，從韌性發(fā)生學視角構(gòu)建了“壓力—狀態(tài)—響應(yīng)”模型，綜合度量了黑土耕地韌性。研究結(jié)果顯示，拜泉縣耕地韌性為0.421，處于中等偏下水平，揭示了黑土區(qū)耕地韌性不足的現(xiàn)狀特征。立地條件壓力對耕地韌性具有顯著負作用，耕地韌性較高的村莊集中于遠離縣城且地形坡度較低的地區(qū)，靠近縣城且地形坡度較大的村莊耕地韌性相對較低。農(nóng)田排澇設(shè)施短缺、耕地勞動力投入不足是拜泉縣耕地韌性的關(guān)鍵短板因素。

耕地保護是一項系統(tǒng)工程，提升耕地韌性需要“內(nèi)外兼修”。首先，坡耕地水蝕（面蝕和溝蝕）導致的水土流失依然是東北黑土區(qū)耕地保護面臨的重要挑戰(zhàn)之一，不僅降低了水土資源質(zhì)量，而且造成了生態(tài)環(huán)境惡化，亟需開展坡耕地區(qū)域排水技術(shù)，構(gòu)建完善的排水系統(tǒng)[26]。東北黑土區(qū)地形具有坡緩、坡長的特點，一般坡度在10°以下，坡長為500～2 000 m，黑土土壤疏松、抗蝕能力弱，長期不合理的開發(fā)利用以及降雨集中，導致了黑土耕地水土流失日趨嚴重。據(jù)2018年水土保持動態(tài)監(jiān)測資料顯示，東北黑土區(qū)水土流失面積為22.16萬km2，占土地總面積的20.48%，其中水蝕面積14.14萬km2，占水土流失總面積的63.8%[27]?；诖耍叫栳槍谕粮氐牡匦?、地貌等條件，通過工程、農(nóng)藝、農(nóng)機及生物措施，加強坡耕地綜合防護與治理，有效控制水土和養(yǎng)分流失、耕地質(zhì)量退化。

其次，落實黑土區(qū)耕地保護，保障糧食生產(chǎn)安全，還應(yīng)進一步強化農(nóng)業(yè)設(shè)施建設(shè)與涉農(nóng)投入。著力加強以農(nóng)田水利為重點的農(nóng)業(yè)基礎(chǔ)設(shè)施建設(shè)與管理維護，在低洼易澇地區(qū)，修建條田化排水、截水排澇設(shè)施，有效減輕積水對農(nóng)作物播種和生長的不利影響，增強東北黑土區(qū)糧食生產(chǎn)抗逆減災(zāi)能力。鼓勵和支持分散農(nóng)戶通過土地轉(zhuǎn)讓、入股、合作社代耕等方式，將土地向農(nóng)業(yè)合作社、家庭農(nóng)場、種糧大戶等新型農(nóng)業(yè)經(jīng)營主體流轉(zhuǎn)，積極開展土地整治工程，提升耕地集中連片水平，推進黑土區(qū)農(nóng)業(yè)規(guī)?；?、機械化與智能化發(fā)展，切實提高耕地的投入產(chǎn)出效益與農(nóng)業(yè)抵御市場風險的能力。

韌性理論為科學認知耕地系統(tǒng)提供了全新視角。針對外界環(huán)境的復(fù)雜性、動態(tài)性以及未來發(fā)展的不確定性和不可預(yù)測性，有效實施耕地保護應(yīng)強調(diào)基于韌性的治理理念。一方面要注重不同地域耕地內(nèi)核系統(tǒng)與外界環(huán)境系統(tǒng)的交互作用方式與過程，診斷識別影響耕地質(zhì)量與生產(chǎn)力的關(guān)鍵短板因素并制定應(yīng)對措施；另一方面要培育增強對發(fā)展環(huán)境變化和外界不確定性的研判能力，及時調(diào)整耕地保護策略，實現(xiàn)耕地保護轉(zhuǎn)型，從而有效規(guī)避風險挑戰(zhàn)或?qū)p失降到最低。