張 玥，代亞強(qiáng)，陳媛媛，柯新利

?土地保障與生態(tài)安全?

中國耕地多功能耦合協(xié)調(diào)時(shí)空演變及其驅(qū)動(dòng)因素

張 玥，代亞強(qiáng)，陳媛媛，柯新利※

（華中農(nóng)業(yè)大學(xué)公共管理學(xué)院，武漢 430070）

為探究耕地多功能耦合協(xié)調(diào)演變規(guī)律及其驅(qū)動(dòng)因素，該研究依據(jù)2000－2018年中國30個(gè)省級(jí)行政區(qū)統(tǒng)計(jì)數(shù)據(jù)，綜合評價(jià)耕地“生產(chǎn)-生活-生態(tài)”三維功能水平并采用核密度估計(jì)和可視化制圖刻畫耕地多功能時(shí)空特征，進(jìn)而采用耦合協(xié)調(diào)度模型分析耕地多功能耦合協(xié)調(diào)關(guān)系，最終采用地理探測器識(shí)別耕地多功能耦合協(xié)調(diào)演變的驅(qū)動(dòng)因素并解析其驅(qū)動(dòng)機(jī)制。研究結(jié)果表明：1）耕地生產(chǎn)功能水平、生活功能水平和生態(tài)功能水平分別處于[0.030, 0.608]、[0.042, 0.672]和[0.058, 0.897]范圍，生產(chǎn)功能呈現(xiàn)上升-下降-上升波動(dòng)變化，生活功能水平先上升后下降，生態(tài)功能水平則整體小幅提升，各項(xiàng)功能均具有極化特征。耕地生產(chǎn)功能、生活功能和生態(tài)功能空間格局分別呈現(xiàn)“東高西低”“北高南低”和“西高東低”分布特征。2）耕地多功能耦合協(xié)調(diào)水平處于[0.093, 0.554]范圍，整體穩(wěn)步提升但水平仍然較低，多數(shù)省份耕地多功能仍為失調(diào)狀態(tài)。耕地多功能耦合協(xié)調(diào)空間格局逐漸呈現(xiàn)集聚提升的演化特征，東北地區(qū)、京津冀地區(qū)和長三角地區(qū)耕地多功能耦合協(xié)調(diào)水平較高，西北地區(qū)則處于較低水平。3）耕地多功能耦合協(xié)調(diào)度演變受經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平影響顯著，其中，經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)力較強(qiáng)，農(nóng)業(yè)發(fā)展水平和工業(yè)發(fā)展水平次之。各驅(qū)動(dòng)因素交互作用類型均為雙因子增強(qiáng)或非線性增強(qiáng)，因素交互作用正向強(qiáng)化了各驅(qū)動(dòng)因素對耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)力。研究結(jié)果為充實(shí)豐富耕地多功能研究提供了實(shí)證指導(dǎo)，同時(shí)為推進(jìn)耕地資源差異化治理和耕地多功能互促提升提供了決策支撐。

土地利用；耕地；多功能；耦合協(xié)調(diào)度；時(shí)空演變；驅(qū)動(dòng)因素

0 引 言

耕地資源是人類賴以生存和發(fā)展的基礎(chǔ)自然資源和重要空間載體，其功能衍生與需求層次升級(jí)和社會(huì)經(jīng)濟(jì)革新密切相關(guān)[1]。伴隨中國由農(nóng)耕文明過渡為工業(yè)文明進(jìn)而走向生態(tài)文明，人的需求也經(jīng)歷著“生存-享受-發(fā)展”的漸次演進(jìn)，耕地多功能已然成為耕地利用的合理方向。從家庭聯(lián)產(chǎn)承包責(zé)任制改革著力解決人民溫飽問題，到改革開放推進(jìn)工業(yè)化和城鎮(zhèn)化進(jìn)程，再到鄉(xiāng)村振興和生態(tài)文明建設(shè)提升人類福祉、促進(jìn)可持續(xù)發(fā)展，耕地利用承擔(dān)著促進(jìn)經(jīng)濟(jì)發(fā)展、維護(hù)社會(huì)穩(wěn)定和保障生態(tài)安全等多重任務(wù)[2-3]。與此同時(shí)，耕地資源的稀缺性日益凸顯，糧食安全、生態(tài)保護(hù)與城鎮(zhèn)發(fā)展的用地矛盾和空間競爭愈發(fā)激烈，嚴(yán)重催化了耕地資源隱性流失，抑制了耕地生產(chǎn)潛能發(fā)揮，同時(shí)削弱了其社會(huì)貢獻(xiàn)價(jià)值并引發(fā)生態(tài)系統(tǒng)退化[4-5]，有違耕地保護(hù)底線原則和永續(xù)利用科學(xué)理念。當(dāng)前，傳統(tǒng)的、單一的耕地生產(chǎn)功能管理模式已難以適配多元化的耕地利用現(xiàn)實(shí)形勢，同時(shí)因其忽視耕地利用復(fù)合價(jià)值、抑制微觀主體保護(hù)動(dòng)力而未能充分助力耕地保護(hù)實(shí)踐。因此，在中國新發(fā)展階段與新發(fā)展格局背景下，耕地多功能協(xié)調(diào)利用不僅是順應(yīng)社會(huì)經(jīng)濟(jì)發(fā)展、滿足人民美好需求的必然選擇，同時(shí)是平衡多方利益關(guān)系、強(qiáng)化耕地保護(hù)實(shí)效的重要手段。對此，積極開展耕地多功能實(shí)證研究與機(jī)制探討，尤其是探清耕地多功能耦合協(xié)調(diào)規(guī)律表征與現(xiàn)實(shí)癥結(jié)，對提升耕地利用綜合效能、促進(jìn)耕地管理創(chuàng)新轉(zhuǎn)型具有重要的科學(xué)意義和現(xiàn)實(shí)意義。

長期以來，耕地多功能研究已逐漸形成“內(nèi)涵界定-現(xiàn)狀評價(jià)-管理調(diào)控”的研究脈絡(luò)，并在宏觀-中觀-微觀不同研究尺度上加以表達(dá)。在內(nèi)涵界定方面，多功能概念起源于農(nóng)業(yè)領(lǐng)域，并逐漸應(yīng)用至土地利用、景觀管理、生態(tài)系統(tǒng)服務(wù)等領(lǐng)域[6-8]。耕地系統(tǒng)是自然、經(jīng)濟(jì)、社會(huì)和生態(tài)子系統(tǒng)高度耦合的復(fù)合系統(tǒng)，耕地多功能是指在復(fù)雜因素交互作用下耕地提供產(chǎn)品和服務(wù)以滿足人類需求的各項(xiàng)能力[9-10]。依據(jù)耕地多功能內(nèi)涵對其類型加以細(xì)分，耕地多功能不僅包括提供糧食、蔬菜等重要農(nóng)作物的生產(chǎn)能力，同時(shí)包括保障農(nóng)民生活、維持農(nóng)村穩(wěn)定等的生活功能和提供氣候調(diào)節(jié)、水源涵養(yǎng)、環(huán)境凈化等生態(tài)系統(tǒng)服務(wù)的生態(tài)能力，兼具阻隔城鎮(zhèn)空間、承載鄉(xiāng)土文化等其他衍生功能[11-12]。在現(xiàn)狀評價(jià)方面，眾多學(xué)者多基于耕地多功能內(nèi)涵界定及類型劃分構(gòu)建評價(jià)體系，采用各類綜合評價(jià)方法、三角模型、空間收斂性法等從多角度揭示耕地多功能時(shí)空演變特征及規(guī)律[13-16]。在管理調(diào)控方面，基于不同研究單元的耕地多功能評價(jià)結(jié)果，耕地多功能管理可綜合運(yùn)用分區(qū)調(diào)控、重點(diǎn)保護(hù)、典型模式推廣等手段進(jìn)行合理引導(dǎo)和科學(xué)調(diào)控[2,17]?，F(xiàn)有耕地多功能研究已初步構(gòu)建從內(nèi)涵認(rèn)知到管理應(yīng)用的研究體系，為本文提供了重要的理論基礎(chǔ)和實(shí)踐支撐。耕地利用系統(tǒng)是經(jīng)濟(jì)-社會(huì)-生態(tài)子系統(tǒng)構(gòu)成的復(fù)合系統(tǒng)，耕地多功能共存于同一區(qū)域?qū)嶓w且具有此消彼長的內(nèi)在關(guān)聯(lián)，因而耕地多功能研究需充分關(guān)注并著力厘清功能結(jié)構(gòu)性、共生性和協(xié)同性特征表達(dá)及其演化邏輯。然而，現(xiàn)有耕地多功能研究多關(guān)注功能時(shí)空分異特征和權(quán)衡協(xié)同關(guān)系，并多在區(qū)域、省級(jí)、市級(jí)尺度開展實(shí)證研究，因而其未能從全局視角充分研判耕地多功能耦合協(xié)調(diào)總體水平及其階段性規(guī)律，而且研究結(jié)論和實(shí)踐建議具有明顯的地域性、典型性特征。耕地多功能耦合協(xié)調(diào)水平是衡量耕地子功能交互影響程度的重要指標(biāo)，其時(shí)序更替性和空間特定性可反映耕地利用綜合效能的實(shí)現(xiàn)程度。中國不同地區(qū)自然資源稟賦、耕地利用模式、社會(huì)經(jīng)濟(jì)水平等各具特點(diǎn)，耕地利用主導(dǎo)功能及其與其他衍生功能關(guān)系不盡相同，針對國家層面省域尺度的宏觀研究有待開展。同時(shí)，鮮少研究定量識(shí)別耕地多功能耦合協(xié)調(diào)的驅(qū)動(dòng)因素，尤其是未能全面考察耕地利用系統(tǒng)“內(nèi)生-外緣”作用力以深入剖析耕地多功能演進(jìn)的內(nèi)在機(jī)理，難以為耕地多功能優(yōu)化管控和統(tǒng)籌協(xié)調(diào)提供充分的決策支撐。耕地利用系統(tǒng)是受社會(huì)經(jīng)濟(jì)環(huán)境影響的復(fù)雜巨系統(tǒng)，要素間相互關(guān)聯(lián)、彼此制約的關(guān)系推動(dòng)著各子功能演化重塑，影響著耕地多功能發(fā)揮程度和協(xié)調(diào)穩(wěn)態(tài)。因此，診斷耕地多功能耦合協(xié)調(diào)演變“內(nèi)生-外緣”驅(qū)動(dòng)因素，深入解析各因素作用機(jī)制，可為促進(jìn)耕地功能協(xié)同融合提供科學(xué)依據(jù)。

鑒于此，本文構(gòu)建“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價(jià)體系，采用2000－2018年中國30個(gè)省級(jí)行政區(qū)面板數(shù)據(jù)開展耕地多功能實(shí)證探究，刻畫并分析耕地多功能時(shí)序特征和空間格局，并采用地理探測器識(shí)別耕地多功能耦合協(xié)調(diào)度時(shí)空演變的驅(qū)動(dòng)因素，以期為制定耕地優(yōu)化配置方案、實(shí)施耕地資源差異化治理提供科學(xué)依據(jù)和政策建議。

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

1.1 研究方法

本研究構(gòu)建“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價(jià)模型，綜合評價(jià)2000－2018年中國30個(gè)省級(jí)行政區(qū)耕地多功能水平，采用核密度估計(jì)和可視化制圖刻畫耕地多功能時(shí)空演進(jìn)特征，進(jìn)而采用耦合協(xié)調(diào)度模型分析耕地多功能耦合協(xié)調(diào)關(guān)系及其演變規(guī)律，最終采用地理探測器識(shí)別耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)因素（圖1）。

1.1.1 耕地多功能評價(jià)模型

耕地系統(tǒng)是自然-人工的復(fù)合系統(tǒng)，具有經(jīng)濟(jì)、社會(huì)、生態(tài)等多重屬性，其功能因人類利用活動(dòng)而存在。中國耕地利用和保護(hù)的重心在于保障糧食供給、維持社會(huì)穩(wěn)定和維護(hù)生態(tài)安全3個(gè)維度，分別對應(yīng)耕地的生產(chǎn)、生活和生態(tài)3種功能。結(jié)合現(xiàn)有研究成果[2,17-18]，依據(jù)系統(tǒng)性、科學(xué)性和數(shù)據(jù)可獲取性等原則，本文選取12項(xiàng)評價(jià)指標(biāo)綜合評價(jià)耕地多功能，指標(biāo)體系及具體釋義如表1所示。

1）生產(chǎn)功能。生產(chǎn)功能是耕地最基本、最核心的功能，反映耕地生產(chǎn)活動(dòng)提供農(nóng)產(chǎn)品的能力[18]。本文綜合考慮耕地利用的實(shí)際產(chǎn)出與耕作狀況，選取糧食單產(chǎn)水平、耕地地均產(chǎn)值、復(fù)種指數(shù)和土地墾殖率4項(xiàng)指標(biāo)予以衡量。其中，糧食單產(chǎn)水平與耕地地均產(chǎn)值分別表征實(shí)物產(chǎn)量水平和經(jīng)濟(jì)產(chǎn)出水平以反映耕地利用的實(shí)際產(chǎn)出，復(fù)種指數(shù)和土地墾殖率則分別表征農(nóng)業(yè)生產(chǎn)活動(dòng)強(qiáng)度和耕地開發(fā)更新程度以反映耕地利用的耕作狀況。

2）生活功能。生活功能是耕地最主要的衍生功能之一，反映耕地保障糧食安全和農(nóng)民生活的能力[19]。為體現(xiàn)耕地利用滿足城鄉(xiāng)居民糧食需求和保障糧食穩(wěn)定自給的糧食安全保障功能，本文選取人均糧食保障率予以反映，具體參照國際慣例，以每年人均400 kg糧食為基準(zhǔn)進(jìn)行測算[2]。為體現(xiàn)耕地利用承載農(nóng)村剩余勞動(dòng)力和滿足促進(jìn)農(nóng)村社會(huì)經(jīng)濟(jì)發(fā)展需求的農(nóng)民生活保障功能，本文選取家庭農(nóng)業(yè)收入比重、農(nóng)業(yè)從業(yè)人員比重和人均農(nóng)業(yè)機(jī)械化水平3項(xiàng)指標(biāo)予以反映。為與下文農(nóng)業(yè)發(fā)展水平驅(qū)動(dòng)因素中的第一產(chǎn)業(yè)從業(yè)人員數(shù)進(jìn)行區(qū)分，此處采用從業(yè)人員比重值衡量耕地生活功能的勞動(dòng)人員承載能力。

3）生態(tài)功能。生態(tài)功能反映耕地維持生物多樣性、提供生態(tài)系統(tǒng)服務(wù)和維護(hù)農(nóng)田生態(tài)健康的能力，不僅與耕地自身生態(tài)本底有關(guān)，同時(shí)受農(nóng)業(yè)生產(chǎn)方式的影響[20-21]。本文選取農(nóng)田生態(tài)系統(tǒng)多樣性和耕地占生態(tài)用地比例2項(xiàng)指標(biāo)反映耕地的生態(tài)本底條件[4,18,22]，其中，農(nóng)田生態(tài)系統(tǒng)多樣性通過各品種作物播種面積與農(nóng)作物播種面積之比加以測算。同時(shí)，本文選取人均耕地生態(tài)承載力和耕地利用化學(xué)負(fù)荷反映耕地的生態(tài)系統(tǒng)服務(wù)供給和農(nóng)田生態(tài)環(huán)境狀況。

表1 耕地多功能評價(jià)指標(biāo)體系

依據(jù)耕地多功能評價(jià)指標(biāo)體系，采用極差標(biāo)準(zhǔn)化法對原始數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化處理，采用熵值法確定指標(biāo)權(quán)重，最后采用線性加權(quán)法計(jì)算得到各省耕地多功能水平。計(jì)算方法如下：

式中為第個(gè)省份第項(xiàng)評價(jià)指標(biāo)標(biāo)準(zhǔn)化數(shù)值；為第項(xiàng)評價(jià)指標(biāo)權(quán)重；為第個(gè)省份的耕地多功能綜合分值，取值區(qū)間為[0,1]，數(shù)值越大，則耕地多功能水平越高。

1.1.2 核密度估計(jì)

核密度估計(jì)（kernel density estimation，KDE）主要優(yōu)勢在于不依賴數(shù)據(jù)分布的先驗(yàn)知識(shí)，不設(shè)置參數(shù)模型的基本假設(shè)，從數(shù)據(jù)本身出發(fā)用連續(xù)的密度函數(shù)描述變量階段分布特征和時(shí)序演進(jìn)規(guī)律[23-24]。因此，借助Eviews8.0軟件，本文采用核密度估計(jì)方法刻畫耕地多功能階段性、動(dòng)態(tài)性特征，具體計(jì)算方法可參考文獻(xiàn)[23]。

1.1.3 耦合協(xié)調(diào)度模型

耦合協(xié)調(diào)度常用以表征兩個(gè)及以上子系統(tǒng)在發(fā)展演化過程中的協(xié)調(diào)一致程度[25]。本文構(gòu)建耕地多功能耦合協(xié)調(diào)度模型，用以定量評價(jià)耕地生產(chǎn)功能、生活功能和生態(tài)功能的交互作用。在耦合協(xié)調(diào)度模型中，耦合度表征子系統(tǒng)之間的相互作用強(qiáng)弱程度；協(xié)調(diào)度表征子系統(tǒng)相互作用關(guān)系中良性互動(dòng)強(qiáng)弱程度[26-27]；最終，耦合協(xié)調(diào)度表征子系統(tǒng)耦合交互的協(xié)同程度。耦合協(xié)調(diào)度具體計(jì)算方法及類型劃分依據(jù)可參考文獻(xiàn)[28-29]。

1.1.4 地理探測器

耕地多功能耦合協(xié)調(diào)度時(shí)空演變受內(nèi)生驅(qū)動(dòng)和外在驅(qū)動(dòng)的雙重影響，識(shí)別耕地多功能耦合協(xié)調(diào)度時(shí)空演變的驅(qū)動(dòng)因素對實(shí)現(xiàn)耕地充分利用和有效保護(hù)具有重要的參考意義。依據(jù)現(xiàn)有研究[4,18,30-32]，本文從經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平4個(gè)方面考察各因素對耕地多功能耦合協(xié)調(diào)度時(shí)空演變的影響程度（表2）。依據(jù)代表性和數(shù)據(jù)可獲取性原則，在經(jīng)營主體條件方面，選取農(nóng)村居民人均純收入和農(nóng)村居民消費(fèi)水平作為表征指標(biāo)，反映耕地綜合利用理念的認(rèn)知程度和功能實(shí)現(xiàn)程度的保障條件[33]；在農(nóng)業(yè)發(fā)展水平方面，選取第一產(chǎn)業(yè)從業(yè)人員數(shù)和第一產(chǎn)業(yè)增加值占比作為表征指標(biāo)，反映農(nóng)業(yè)發(fā)展集約化、現(xiàn)代化水平。為與上文耕地多功能評價(jià)指標(biāo)體系中的農(nóng)業(yè)從業(yè)人員比重進(jìn)行區(qū)分，此處采用第一產(chǎn)業(yè)從業(yè)人員數(shù)和增加值占比共同反映農(nóng)業(yè)發(fā)展水平驅(qū)動(dòng)作用；在工業(yè)發(fā)展水平方面，選取第二產(chǎn)業(yè)增加值占比和規(guī)模以上工企業(yè)個(gè)數(shù)作為表征指標(biāo)，反映工業(yè)發(fā)展現(xiàn)狀對耕地多功能發(fā)揮的經(jīng)濟(jì)支撐和需求引導(dǎo)；在城鎮(zhèn)化水平方面，選取城鎮(zhèn)人口占比和城鎮(zhèn)道路面積作為表征指標(biāo)，反映城鎮(zhèn)化發(fā)展對耕地資源本底條件及其功能水平的外力改造和強(qiáng)度加載。

表2 耕地多功能耦合協(xié)調(diào)度驅(qū)動(dòng)因素指標(biāo)選取

地理探測器（geographical detector）是探析地理現(xiàn)象空間分異形成機(jī)理的重要方法，其優(yōu)勢在于不設(shè)置線性假設(shè)和條件限制，可客觀探測定性因子和定量因子的驅(qū)動(dòng)作用及不同因子間的交互作用[34]。因此，本文首先運(yùn)用自然斷點(diǎn)法對驅(qū)動(dòng)因素進(jìn)行離散化和類別化處理，進(jìn)而采用地理探測器定量識(shí)別各驅(qū)動(dòng)因素作用程度及因素交互作用。地理探測器具體計(jì)算過程和驅(qū)動(dòng)因素交互作用類型判別依據(jù)可參考文獻(xiàn)[34]。

1.2 研究區(qū)域與數(shù)據(jù)來源

由于西藏、香港、澳門及臺(tái)灣數(shù)據(jù)缺失，本文選取中國30個(gè)省級(jí)行政區(qū)為研究單元，其中包括22個(gè)省、4個(gè)自治區(qū)、4個(gè)直轄市。本文使用的數(shù)據(jù)包括統(tǒng)計(jì)數(shù)據(jù)與行政區(qū)劃數(shù)據(jù)。其中，統(tǒng)計(jì)數(shù)據(jù)均來源于2001－2019年《中國統(tǒng)計(jì)年鑒》《中國農(nóng)村統(tǒng)計(jì)年鑒》《中國環(huán)境統(tǒng)計(jì)年鑒》、中國環(huán)境數(shù)據(jù)庫。因上述指標(biāo)相鄰年份數(shù)據(jù)具有關(guān)聯(lián)性、較少發(fā)生異常變動(dòng)，因此個(gè)別年份缺失數(shù)據(jù)采用趨勢外推或插值法進(jìn)行近似獲取[35]。行政區(qū)劃數(shù)據(jù)來源于中科院地理科學(xué)與資源研究所資源環(huán)境科學(xué)與數(shù)據(jù)中心（http:// www. resdc.cn/）。

2 結(jié)果與分析

2.1 耕地多功能評價(jià)

2.1.1 時(shí)序演變特征

依據(jù)耕地各功能水平評價(jià)結(jié)果，繪制2000－2018年耕地生產(chǎn)、生活和生態(tài)功能水平核密度曲線圖（圖2）。1）從耕地生產(chǎn)功能核密度曲線來看，曲線重心經(jīng)歷“右-左-右”遷移軌跡，曲線始終為單一波峰且高度逐漸下降，表明研究期內(nèi)中國耕地生產(chǎn)功能水平整體呈現(xiàn)上升-下降-上升的波動(dòng)趨勢，各省生產(chǎn)功能處于極化狀態(tài)且差異有所擴(kuò)大；2）從耕地生活功能核密度曲線來看，2000－2015年重心持續(xù)向右遷移，2015－2018年重心輕微向左遷移，曲線保持單一波峰且高度持續(xù)下降，表明研究期內(nèi)中國耕地生活功能水平整體呈現(xiàn)先上升后下降的趨勢，各省生活功能呈現(xiàn)兩極分化狀態(tài)且差異逐漸擴(kuò)大；3）從耕地生態(tài)功能核密度曲線來看，曲線重心輕微向右遷移，曲線單一波峰于2000－2015年持續(xù)下降，后于2015－2018年上升，表明中國耕地生態(tài)功能水平整體呈現(xiàn)小幅提升，省域間耕地生態(tài)功能水平差異呈現(xiàn)先擴(kuò)大后縮小的特征。

圖2 2000－2018年中國耕地多功能時(shí)序特征

2.1.2 空間格局特征

為直觀展現(xiàn)各功能空間分布，繪制耕地生產(chǎn)功能、生活功能和生態(tài)功能空間格局圖。從耕地生產(chǎn)功能來看，研究期內(nèi)中國耕地生產(chǎn)功能水平處于[0.030, 0.608]范圍，高水平區(qū)集中位于黃淮海平原并逐漸向東北平原和長江中下游地區(qū)延伸，低水平區(qū)則分散位于西北地區(qū)和西南地區(qū)，大體呈現(xiàn)“東高西低”的空間非均衡格局（圖3）。具體來看：1）2000年，耕地生產(chǎn)功能高水平區(qū)集中在黃淮海平原，得益于其平坦廣闊的地形和較為完善的農(nóng)業(yè)設(shè)施；2）2005－2015年，耕地生產(chǎn)功能高水平區(qū)集中在黃淮海平原和長江中下游地區(qū)，四川盆地耕地生產(chǎn)功能有所下降而云貴地區(qū)、華南部分地區(qū)則有所提升；3）2018年，耕地生產(chǎn)功能高水平區(qū)集中位于東北平原、黃淮海平原和長江中下游地區(qū)，低水平區(qū)位于西北地區(qū)和西南地區(qū)，耕地生產(chǎn)功能受限于崎嶇地形、水熱條件較差、農(nóng)業(yè)設(shè)施等條件[36-37]?？傮w而言，研究期內(nèi)耕地生產(chǎn)功能空間格局以黃淮海平原、長江中下游地區(qū)等重要農(nóng)業(yè)區(qū)為核心輻射延展，耕地生產(chǎn)功能整體水平有所提升。

從耕地生活功能來看，研究期內(nèi)中國耕地生活功能水平處于[0.042, 0.672]范圍，高水平區(qū)分散位于黃淮海平原、東北平原和新疆地區(qū)，隨后延伸至江漢平原和長江中下游地區(qū)，低水平區(qū)則位于京津地區(qū)、云貴地區(qū)和華南地區(qū)，大體呈現(xiàn)“北高南低”的空間非均衡格局（圖4）。1）2000－2005年，耕地生活功能高水平省份為新疆、黑龍江和江蘇，新疆產(chǎn)業(yè)結(jié)構(gòu)以第一產(chǎn)業(yè)為主，其余兩省為重要的糧食主產(chǎn)區(qū)，耕地資源可為農(nóng)民生活提供較大保障；2）2010－2015年，耕地生活功能高水平區(qū)以上述三省為核心延伸至北方干旱半干旱區(qū)、東北平原、黃淮海地區(qū)和長江中下游北部地區(qū)；3）2018年，耕地生活功能較高水平區(qū)向長江中下游南部地區(qū)拓展，該空間格局與耕地生產(chǎn)功能高水平區(qū)布局具有相似性。研究期內(nèi)，耕地生活功能低水平區(qū)分散位于云貴地區(qū)、華南地區(qū)和西北三省。其中，華南地區(qū)經(jīng)濟(jì)發(fā)達(dá)、農(nóng)民生計(jì)選擇多樣，云貴地區(qū)和西北三省人口外遷嚴(yán)重，耕地生活保障作用較為薄弱。

從耕地生態(tài)功能來看，研究期內(nèi)中國耕地生態(tài)功能水平處于[0.058, 0.897]范圍，高水平區(qū)集中位于東北地區(qū)和西南地區(qū)并逐漸延伸至黃淮海平原和四川盆地及周邊地區(qū)，低水平區(qū)則位于京津冀地區(qū)和長江中部地區(qū)，大體呈現(xiàn)“西高東低”的空間格局（圖5）。具體來看：1）2000－2005年，耕地生態(tài)功能高水平區(qū)集中位于東北地區(qū)和西南地區(qū)，其中，東北地區(qū)耕地資源豐富、規(guī)模經(jīng)營較為成熟，西南地區(qū)化肥、農(nóng)藥等農(nóng)資投放量較少。同時(shí)，內(nèi)蒙古人均耕地資源較多且規(guī)模經(jīng)營范圍較廣，因而耕地生態(tài)功能穩(wěn)定保持高水平。耕地生態(tài)功能低水平區(qū)則位于京津冀地區(qū)、長三角地區(qū)和長江中部地區(qū)，上述地區(qū)耕地資源約束趨緊且利用強(qiáng)度增大，因而其耕地生態(tài)系統(tǒng)愈加脆弱、生態(tài)功能逐漸退化；2）2010－2015年，耕地生態(tài)功能高水平區(qū)向黃淮海平原和四川盆地及周邊地區(qū)拓展，同時(shí)在生態(tài)文明理念倡導(dǎo)與經(jīng)濟(jì)發(fā)展基礎(chǔ)保障的雙重作用下，低水平區(qū)生態(tài)逐漸改善；3）2018年，耕地生態(tài)功能仍然保持“西高東低”的格局特征。

注：西藏、香港、澳門、臺(tái)灣數(shù)據(jù)暫缺。下同。

圖4 2000－2018年中國耕地生活功能空間格局

2.2 耕地多功能耦合協(xié)調(diào)度評價(jià)

從耕地多功能耦合協(xié)調(diào)時(shí)序演變來看，研究期內(nèi)中國耕地多功能耦合協(xié)調(diào)水平處于[0.093, 0.554]范圍，整體呈現(xiàn)穩(wěn)步提升趨勢但水平仍然較低，多數(shù)省份耕地多功能為失調(diào)狀態(tài)，表明中國耕地多功能正在由無序沖突轉(zhuǎn)為有序發(fā)展，但耕地多功能耦合協(xié)調(diào)水平仍存在較大提升空間。從耕地多功能耦合協(xié)調(diào)空間格局來看，研究期內(nèi)中國耕地多功能耦合協(xié)調(diào)度逐漸呈現(xiàn)集聚提升的演化特征（圖6）。具體來看：1）2000年，除內(nèi)蒙古、黑龍江和江蘇耕地多功能為輕度失調(diào)外，中度失調(diào)省份多位于京津冀地區(qū)、東南沿海地區(qū)和中部地區(qū)，嚴(yán)重失調(diào)和極度失調(diào)省份位于西北地區(qū)；2）2005－2015年，吉林耕地多功能提升至勉強(qiáng)協(xié)調(diào)類型，黑龍江和上海提升至瀕臨失調(diào)類型，輕度失調(diào)省份多位于四川盆地和長江中下游地區(qū)；3）2018年，東北地區(qū)耕地多功能耦合協(xié)調(diào)水平較高，西北地區(qū)雖有所提升但耦合協(xié)調(diào)水平仍然較低。由此可見，耕地多功能退耦現(xiàn)象有所改善，多功能逐漸形成協(xié)同提升的良性演化關(guān)系。

圖5 2000－2018年中國耕地生態(tài)功能空間格局

圖6 2000－2018年中國耕地多功能耦合協(xié)調(diào)度空間格局

耕地多功能是耕地利用的內(nèi)在屬性和客觀需求，其實(shí)質(zhì)是耕地系統(tǒng)與經(jīng)濟(jì)、社會(huì)和環(huán)境持續(xù)交換物質(zhì)、能量與信息的過程[9]。耕地多功能是單項(xiàng)功能關(guān)聯(lián)反饋、互動(dòng)融合而形成的有機(jī)整體，并非單項(xiàng)功能的獨(dú)立發(fā)揮與機(jī)械疊加。在耕地資源長期稀缺的前提下，耕地利用需求不斷衍生引致耕地功能競爭消長、損益平衡，宏觀政策制定執(zhí)行與微觀農(nóng)戶決策行動(dòng)均折射為耕地多功能耦合協(xié)調(diào)水平的動(dòng)態(tài)演化和區(qū)域異質(zhì)分布[9]。京津冀地區(qū)、長三角地區(qū)等經(jīng)濟(jì)發(fā)達(dá)區(qū)農(nóng)業(yè)機(jī)械化、專業(yè)化程度較高，同時(shí)，其綠色的市場消費(fèi)導(dǎo)向和先進(jìn)的生態(tài)保護(hù)理念進(jìn)一步促進(jìn)耕地多功能有序分化[11,38]。東北地區(qū)具有優(yōu)越的資源稟賦條件和農(nóng)業(yè)生產(chǎn)基礎(chǔ)，加之其積極推廣采用保護(hù)性耕作技術(shù)和治理修復(fù)措施，因而該地區(qū)耕地多功能水平較高。西北地區(qū)耕地利用與改造難度較大，同時(shí)該地區(qū)技術(shù)推廣較為滯后、市場化程度提升緩慢，相對滯后的耕地利用認(rèn)知程度和改造能力阻礙了耕地多功能優(yōu)化提升，仍需合理引導(dǎo)和積極促進(jìn)。

2.3 耕地多功能耦合協(xié)調(diào)度驅(qū)動(dòng)因素

2.3.1 單因子驅(qū)動(dòng)探測

在揭示中國耕地多功能耦合協(xié)調(diào)度演變規(guī)律的基礎(chǔ)上，采用地理探測器模型探析各驅(qū)動(dòng)因素對耕地多功能耦合協(xié)調(diào)度的驅(qū)動(dòng)力（表3）。各驅(qū)動(dòng)因素均通過1%的顯著性水平檢驗(yàn)，表明經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調(diào)度演變均具有驅(qū)動(dòng)作用，但不同因素的驅(qū)動(dòng)力不盡相同，經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)力較強(qiáng)，農(nóng)業(yè)發(fā)展水平和工業(yè)發(fā)展水平的驅(qū)動(dòng)力次之。

2.3.2 雙因子交互探測

在識(shí)別耕地多功能耦合協(xié)調(diào)度演變驅(qū)動(dòng)因素的基礎(chǔ)上，進(jìn)一步探究各驅(qū)動(dòng)因素交互作用對其影響（圖7）。具體如下：1）研究期內(nèi)，各驅(qū)動(dòng)因素交互作用類型均為雙因子增強(qiáng)或非線性增強(qiáng)，表明雙因素交互作用正向強(qiáng)化了各驅(qū)動(dòng)因素對耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)力。同時(shí)，雙因子增強(qiáng)類型由8組減至2組，非線性增強(qiáng)類型由20組增至26組，表明驅(qū)動(dòng)因素交互作用漸趨增強(qiáng)、合力作用趨勢逐漸顯現(xiàn)；2）農(nóng)村居民人均純收入（1）、農(nóng)村居民消費(fèi)水平（2）、第一產(chǎn)業(yè)從業(yè)人員數(shù)（3）與其他驅(qū)動(dòng)因素交互作用有所增強(qiáng)，第一產(chǎn)業(yè)增加值占比（4）、第二產(chǎn)業(yè)增加值占比（5）與其他驅(qū)動(dòng)因素交互作用有所減弱，進(jìn)一步驗(yàn)證，經(jīng)營主體條件和城鎮(zhèn)化水平對耕地多功能耦合協(xié)調(diào)度演變的驅(qū)動(dòng)力較為主導(dǎo)，兩者與農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平共同作用引致耕地多功能耦合協(xié)調(diào)度呈現(xiàn)上述演變規(guī)律。因此，本文將從經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平四個(gè)方面具體闡述耕地多功能耦合協(xié)調(diào)度演變驅(qū)動(dòng)機(jī)制。

表3 耕地多功能耦合協(xié)調(diào)演變驅(qū)動(dòng)因素探測結(jié)果

圖7 耕地多功能耦合協(xié)調(diào)演變驅(qū)動(dòng)因素交互作用探測結(jié)果

2.3.3 驅(qū)動(dòng)機(jī)制解析

耕地系統(tǒng)是開放的、受人為干預(yù)明顯的復(fù)雜巨系統(tǒng)，其多功能耦合協(xié)調(diào)度演變受社會(huì)經(jīng)濟(jì)諸多要素影響，即內(nèi)生驅(qū)動(dòng)和外生驅(qū)動(dòng)綜合作用（圖8），具體驅(qū)動(dòng)機(jī)制如下：

1）經(jīng)營主體條件反映農(nóng)村居民對耕地綜合利用理念的認(rèn)知程度和對耕地綜合利用開展的保障條件，是影響耕地功能發(fā)揮的基礎(chǔ)因素。一方面，農(nóng)村居民開始適應(yīng)并追求多元化從業(yè)方式，農(nóng)村剩余勞動(dòng)力加快析出并逐步獲取穩(wěn)定就業(yè)機(jī)會(huì)，以致農(nóng)村居民逐漸脫離對耕地資源的生存依賴，有助于激發(fā)耕地流轉(zhuǎn)行為，實(shí)現(xiàn)規(guī)?；?jīng)營，對耕地生產(chǎn)功能和生活功能提升具有促進(jìn)作用[32,39]；另一方面，農(nóng)村居民感知市場需求信息和消費(fèi)升級(jí)導(dǎo)向的敏銳性增強(qiáng)，逐漸關(guān)注耕地生態(tài)功能提升和生態(tài)效益挖掘，進(jìn)一步促進(jìn)了耕地多功能相互融合與良性互饋。

2）農(nóng)業(yè)發(fā)展水平反映農(nóng)業(yè)發(fā)展集約化、現(xiàn)代化水平，是影響耕地功能發(fā)揮的重要因素之一。伴隨著傳統(tǒng)農(nóng)業(yè)發(fā)展理念的逐步改觀，現(xiàn)代農(nóng)業(yè)不僅要求耕地資源發(fā)揮糧食供給作用，同時(shí)提倡耕地利用經(jīng)濟(jì)-社會(huì)-生態(tài)效益的統(tǒng)籌兼顧。農(nóng)業(yè)發(fā)展水平提高為耕地多功能利用創(chuàng)造了良好的基礎(chǔ)保障，促進(jìn)多功能協(xié)調(diào)發(fā)展并正向反饋于農(nóng)業(yè)可持續(xù)發(fā)展。但是，農(nóng)業(yè)發(fā)展水平加速了農(nóng)業(yè)機(jī)械化推廣與水利化改造建設(shè)，雖有助于高效利用耕地，但一定程度上破壞了耕地本底生態(tài)，加重了耕地利用負(fù)擔(dān)，進(jìn)而影響生態(tài)功能的保持和發(fā)揮，以致耕地多功能難以有效協(xié)同。

3）工業(yè)發(fā)展水平反映工業(yè)發(fā)展對耕地多功能顯化的支撐能力和引導(dǎo)作用，是影響耕地功能發(fā)揮的誘導(dǎo)因素。工業(yè)發(fā)展刺激消費(fèi)需求增長和消費(fèi)結(jié)構(gòu)升級(jí)，誘發(fā)耕地利用投入結(jié)構(gòu)、種植結(jié)構(gòu)、經(jīng)營方式等發(fā)生改變[30]，驅(qū)動(dòng)傳統(tǒng)耕地利用方式轉(zhuǎn)型。另一方面，工業(yè)快速發(fā)展會(huì)攫取甚至破壞耕地資源，對耕地資源數(shù)量和質(zhì)量均產(chǎn)生負(fù)面影響。為穩(wěn)定并提高有限耕地資源的經(jīng)濟(jì)產(chǎn)出，大量追加化肥、農(nóng)藥等化學(xué)生產(chǎn)要素是促進(jìn)生產(chǎn)的直接手段之一。但是，高強(qiáng)度、高能耗的耕地利用方式抑制了耕地生態(tài)功能的發(fā)揮，從而影響耕地多功能全面協(xié)調(diào)發(fā)展。

4）城鎮(zhèn)化水平反映城鎮(zhèn)化發(fā)展對耕地利用的改造方向和強(qiáng)度，是影響耕地功能發(fā)揮的關(guān)鍵因素。城鎮(zhèn)化發(fā)展逐漸破除城鄉(xiāng)發(fā)展壁壘，城鎮(zhèn)人口聚集增加和生活水平改善提升對耕地利用提出了更高品質(zhì)的效益需求，推動(dòng)了耕地多功能協(xié)調(diào)發(fā)展[4,30]。但同時(shí)，城市空間擴(kuò)張迫使耕地向城鎮(zhèn)建設(shè)用地轉(zhuǎn)變，對耕地生產(chǎn)本底造成了負(fù)面影響。另外，城鎮(zhèn)化發(fā)展也會(huì)導(dǎo)致農(nóng)業(yè)勞動(dòng)力質(zhì)量降低、農(nóng)業(yè)發(fā)展活力喪失，雖耕地多重效益需求膨脹但多功能利用開展存在一定難度，因此城鎮(zhèn)化發(fā)展會(huì)對耕地多功能培育和發(fā)展產(chǎn)生負(fù)向影響。

圖8 耕地多功能耦合協(xié)調(diào)演變驅(qū)動(dòng)機(jī)制分析框架

3 討 論

耕地是自然-人工充分交互的復(fù)雜巨系統(tǒng)，其多功能發(fā)揮具有明顯的階段分化和地域分異特征。因此，本文從時(shí)空二維視角綜合評價(jià)耕地多功能水平及其演變規(guī)律，同時(shí)進(jìn)一步考察耕地多功能耦合協(xié)調(diào)關(guān)系，通過定量探測驅(qū)動(dòng)因素和定性解構(gòu)驅(qū)動(dòng)機(jī)制相結(jié)合的方法深入探究耕地多功能耦合協(xié)調(diào)演變的作用機(jī)理，嘗試遵循“表征刻畫-動(dòng)因識(shí)別-機(jī)理剖析”的研究邏輯豐富耕地多功能理性認(rèn)識(shí)。耕地多功能為耕地利用提供了合理方向并為耕地價(jià)值顯化提供了具體路徑，耕地多功能管理應(yīng)結(jié)合資源優(yōu)勢與發(fā)展環(huán)境以制定精細(xì)化、差異化的管控方案及保障措施。

在“高效利用+高度保護(hù)”的可持續(xù)理念引導(dǎo)下，耕地多功能管理不僅應(yīng)該針對性滿足區(qū)域社會(huì)經(jīng)濟(jì)增長和居民生活改善對農(nóng)業(yè)和耕地的適度要求，充分利用社會(huì)多元需求和現(xiàn)代化發(fā)展提供的有利機(jī)遇顯化自身價(jià)值，同時(shí)應(yīng)強(qiáng)調(diào)耕地質(zhì)量提升和生態(tài)系統(tǒng)穩(wěn)定以實(shí)現(xiàn)耕地永續(xù)利用。在耕地多功能利用發(fā)展較為活躍的地區(qū)，政府和各經(jīng)營主體應(yīng)充分利用地區(qū)資源優(yōu)勢和市場優(yōu)勢，通過釋放市場消費(fèi)潛力、激發(fā)耕地多功能利用活力。在耕地資源有限但市場需求旺盛的背景下，該類地區(qū)應(yīng)充分利用發(fā)展耕地多功能的先行條件，著力打造生態(tài)綠色農(nóng)業(yè)、現(xiàn)代設(shè)施農(nóng)業(yè)、典型特色農(nóng)業(yè)等，通過推動(dòng)多樣化、高端化、綠色化的農(nóng)業(yè)產(chǎn)業(yè)融合，保障經(jīng)營主體收益、助力耕地多功能價(jià)值顯化。同時(shí)，這類地區(qū)應(yīng)自下而上地推進(jìn)耕地多功能管理制度完善和路徑創(chuàng)新，為耕地可持續(xù)利用提供良好的制度支撐。在耕地多功能耦合協(xié)調(diào)水平較低的地區(qū)，決策者可以通過加強(qiáng)規(guī)劃布局引導(dǎo)、指標(biāo)控制、經(jīng)濟(jì)補(bǔ)貼、監(jiān)測預(yù)警等多種方式實(shí)行耕地多功能管理，在不降低耕地商品性生產(chǎn)功能的前提下，激勵(lì)農(nóng)戶參與耕地衍生功能的培育和挖掘，發(fā)揮農(nóng)戶保護(hù)與利用耕地的主觀能動(dòng)性，因勢利導(dǎo)、因地制宜地實(shí)現(xiàn)耕地利用轉(zhuǎn)型。

4 結(jié) 論

本文以2000－2018年中國30個(gè)省級(jí)行政區(qū)為研究單元，開展“生產(chǎn)-生活-生態(tài)”三維耕地多功能評價(jià)并采用核密度估計(jì)和可視化制圖刻畫耕地多功能時(shí)空演進(jìn)特征，進(jìn)而采用耦合協(xié)調(diào)度模型分析耕地多功能耦合協(xié)調(diào)關(guān)系及其演變規(guī)律，最終采用地理探測器識(shí)別耕地多功能耦合協(xié)調(diào)演變的驅(qū)動(dòng)因素。研究主要結(jié)論如下：

1）從時(shí)序演變特征來看，耕地生產(chǎn)功能水平呈現(xiàn)上升-下降-上升的波動(dòng)變化，生活功能水平先上升后下降，生態(tài)功能水平則整體小幅提升，各功能均具有極化特征。從空間格局特征來看，耕地生產(chǎn)功能、生活功能和生態(tài)功能分別呈現(xiàn)“東高西低”“北高南低”和“西高東低”的空間非均衡格局。

2）從時(shí)序演變來看，中國耕地多功能耦合協(xié)調(diào)水平整體穩(wěn)步提升但水平仍然較低，多數(shù)省份耕地多功能為失調(diào)狀態(tài)，表明中國耕地多功能正在由無序沖突轉(zhuǎn)為有序發(fā)展，但耕地多功能耦合協(xié)調(diào)水平仍需進(jìn)一步提升。從空間格局來看，中國耕地多功能耦合協(xié)調(diào)度逐漸呈現(xiàn)集聚提升的演化特征。

3）耕地多功能耦合協(xié)調(diào)度演變受經(jīng)營主體條件、農(nóng)業(yè)發(fā)展水平、工業(yè)發(fā)展水平和城鎮(zhèn)化水平影響顯著，但不同因素的驅(qū)動(dòng)力不盡相同。各驅(qū)動(dòng)因素交互作用類型均為雙因子增強(qiáng)或非線性增強(qiáng)，且非線性增強(qiáng)類型逐漸增多，表明耕地多功能耦合協(xié)調(diào)演變?yōu)閮?nèi)生驅(qū)動(dòng)和外在驅(qū)動(dòng)雙重影響的結(jié)果。

本文針對長時(shí)序、省級(jí)尺度耕地多功能開展實(shí)證研究，受限于統(tǒng)計(jì)數(shù)據(jù)和調(diào)研資料的可獲取性，并未于更細(xì)尺度、更長時(shí)序揭示耕地多功能時(shí)空規(guī)律及演化機(jī)理，也未能將耕地多功能評價(jià)結(jié)果與農(nóng)戶走訪調(diào)查情況進(jìn)行對比修正。因此，未來研究將采用多源數(shù)據(jù)、多種技術(shù)手段全面考察并深入挖掘耕地多功能演化表征和驅(qū)動(dòng)機(jī)制，同時(shí)借助專家訪談、實(shí)地調(diào)研等多種研究方法輔助修正研究結(jié)果，以期充分揭示中國耕地多功能結(jié)構(gòu)性、趨勢性、系統(tǒng)性的發(fā)展規(guī)律和內(nèi)在機(jī)理，為實(shí)現(xiàn)耕地資源優(yōu)化配置和功能協(xié)同提升提供決策支撐。

[1] 鄒利林，李裕瑞，劉彥隨，等. 基于要素視角的耕地“三生”功能理論建構(gòu)與實(shí)證研究[J]. 地理研究，2021，40(3)：839-855.

ZOU Lilin, LI Yurui, LIU Yansui, et al. Theory building and empirical research of production-living-ecological function of cultivated land based on the elements[J]. Geographical Research, 2021, 40(3): 839-855. (in Chinese with English abstract)

[2] 熊昌盛，張永蕾，王雅娟，等. 中國耕地多功能評價(jià)及分區(qū)管控[J]. 中國土地科學(xué)，2021，35(10)：104-114.

XIONG Changsheng, ZHANG Yonglei, WANG Yajuan, et al. Multi-function evaluation and zoning control of cultivated land in China[J]. Chinese Land Science, 2021, 35(10): 104-114. (in Chinese with English abstract)

[3] 朱從謀，李武艷，杜瑩瑩，等. 浙江省耕地多功能價(jià)值時(shí)空變化與權(quán)衡-協(xié)同關(guān)系[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(14)：263-272.

ZHU Congmou, LI Wuyan, DU Yingying, et al. Spatial-temporal change, trade-off and synergy relationships of cropland multifunctional value in Zhejiang Province, China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(14): 263-272. (in Chinese with English abstract)

[4] 張英男，龍花樓，戈大專，等. 黃淮海平原耕地功能演變的時(shí)空特征及其驅(qū)動(dòng)機(jī)制[J]. 地理學(xué)報(bào)，2018，73(3)：518-534.

ZHANG Yingnan, LONG Hualou, GE Dazhuan, et al. Spatio-temporal characteristics and dynamic mechanism of farmland functions evolution in the Huang-Huai-Hai Plain[J]. Acta Geographica Sinica, 2018, 73(3): 518-534. (in Chinese with English abstract)

[5] 安悅，譚雪蘭，李印齊，等. 洞庭湖地區(qū)耕地功能時(shí)空演變特征及影響因素研究[J]. 地理科學(xué)，2022，42(7)：1272-1282.

AN Yue, TAN Xuelan, LI Yinqi, et al.Spatio-temporal evolution characteristics and influencing factors of cultivated land functions in the Dongting Lake Aera[J]. Scientia Geographica Sinica, 2022, 42(7): 1272-1282. (in Chinese with English abstract)

[6] 宋小青，歐陽竹. 耕地多功能內(nèi)涵及其對耕地保護(hù)的啟示[J]. 地理科學(xué)進(jìn)展，2012，31(7)：859-868.

SONG Xiaoqing, OUYANG Zhu. Connotation of multifunctional cultivated land and its implications for cultivated land protection[J]. Progress in Geography, 2012, 31(7): 859-868. (in Chinese with English abstract)

[7] HELMING K, PEREZ-SOBA M, TABBUSH P. Sustainability Impact Assessment of Land Use Changes[M]. Berlin: Springer, 2008: 375-404.

[8] MANDER U, HELMING K, WIGGERING H. Multifunctional Land Use: Meeting Future Demands for Landscape Goods and Services[M]. Berlin: Springer, 2007: 1-13.

[9] 宋小青，歐陽竹. 中國耕地多功能管理的實(shí)踐路徑探討[J]. 自然資源學(xué)報(bào)，2012，27(4)：540-551.

SONG Xiaoqing, OUYANG Zhu. Route of multifunctional cultivated land management in China[J]. Journal of Natural Resources, 2012, 27(4): 540-551. (in Chinese with English abstract)

[10] 宋小青，李心怡. 區(qū)域耕地利用功能轉(zhuǎn)型的理論解釋與實(shí)證[J]. 地理學(xué)報(bào)，2019，74(5)：992-1010.

SONG Xiaoqing, LI Xinyi. Theoretical explanation and case study of regional cultivated land use function transition[J]. Acta Geographica Sinica, 2019, 74(5): 992-1010. (in Chinese with English abstract)

[11] 范業(yè)婷，金曉斌，項(xiàng)曉敏，等. 蘇南地區(qū)耕地多功能評價(jià)與空間特征分析[J]. 資源科學(xué)，2018，40(5)：980-992.

FAN Yeting, JIN Xiaobin, XIANG Xiaomin, et al. Evaluation and spatial characteristics of arable land multifunction in southern Jiangsu[J]. Resources Science, 2018, 40(5): 980-992. (in Chinese with English abstract)

[12] WANG X J, WANG D Y, WU S Z, et al. Cultivated land multifunctionality in undeveloped peri-urban agriculture areas in China: Implications for sustainable land management[J]. Journal of Environmental Management, 2023, 325: 116500.

[13] ZHANG S Y, HU W Y, LI M R, et al. Multiscale research on spatial supply-demand mismatches and synergic strategies of multifunctional cultivated land[J]. Journal of Environmental Management, 2021, 299: 113605.

[14] JIANG G H, WANG M Z, QU Y B, et al. Towards cultivated land multifunction assessment in China: Applying the “influencing factors-functions-products-demands” integrated framework[J]. Land Use Policy, 2020, 99: 104982.

[15] 馬才學(xué)，金瑩，柯新利，等. 基于全排列多邊形圖示法的湖北省耕地多功能強(qiáng)度與協(xié)調(diào)度典型模式探究[J]. 中國土地科學(xué)，2018，32(4)：51-58.

MA Caixue, JIN Ying, KE Xinli, et al. Study on the typical modes of farmland multi-functional intensity and coordination in Hubei province based on the entire-array-polygon indicator method[J]. Chinese Land Science, 2018, 32(4): 51-58. (in Chinese with English abstract)

[16] 向敬偉，廖曉莉，宋小青，等. 中國耕地多功能的區(qū)域收斂性[J]. 資源科學(xué)，2019，41(11)：1959-1971.

XIANG Jingwei, LIAO Xiaoli, SONG Xiaoqing, et al. Regional convergence of cultivated land multifunctions in China[J]. Resources Science, 2019, 41(11): 1959-1971. (in Chinese with English abstract)

[17] 方瑩，王靜，孔雪松，等. 耕地利用多功能權(quán)衡關(guān)系測度與分區(qū)優(yōu)化——以河南省為例[J]. 中國土地科學(xué)，2018，32(11)：57-64.

FANG Ying, WANG Jing, KONG Xuesong, et al. Trade-off relation measurement and zoning optimization of multi-functionality of cultivated land use: A case study of Henan province[J]. Chinese Land Science, 2018, 32(11): 57-64. (in Chinese with English abstract)

[18] 宋小青，吳志峰，歐陽竹. 1949年以來中國耕地功能變化[J]. 地理學(xué)報(bào)，2014，69(4)：435-447.

SONG Xiaoqing, WU Zhifeng, OUYANG Zhu. Changes of cultivated land function in China since 1949[J]. Acta Geographica Sinica, 2014, 69(4): 435-447. (in Chinese with English abstract)

[19] 胡偉艷，朱慶瑩，張安錄，等. 總量與結(jié)構(gòu)視角耕地多功能對農(nóng)業(yè)經(jīng)濟(jì)增長的影響：以湖北省為例[J]. 中國土地科學(xué)，2018，32(5)：62-70.

HU Weiyan, ZHU Qingying, ZHANG Anlu, et al. Impacts of multifunctional farmland on agricultural economic growth at county-level from the perspective of quantity and structure: Evidence from Hubei province[J]. Chinese Land Science, 2018, 32(5): 62-70. (in Chinese with English abstract)

[20] 許多藝，濮勵(lì)杰，黃思華，等. 江蘇省耕地多功能時(shí)空動(dòng)態(tài)分析及對耕地?cái)?shù)量變化響應(yīng)研究[J]. 長江流域資源與環(huán)境，2022，31(3)：575-587.

XU Duoyi, PU Lijie, HUANG Sihua, et al. Spatial and temporal dynamic analysis of cultivated land multifunction in Jiangsu Province and its response to cultivated land change[J]. Resources and Environment in the Yangtze Basin, 2022, 31(3): 575-587. (in Chinese with English abstract)

[21] 高星，宋昭穎，李晨曦，等. 城鄉(xiāng)梯度下的耕地多功能價(jià)值空間分異特征[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(16)：251-259.

GAO Xing, SONG Zhaoying, LI Chenxi, et al. Spatial differentiation characteristics of cultivated land multifunctional value under urban-rural gradient[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(16): 251-259. (in Chinese with English abstract)

[22] 盧新海，唐一峰，易家林，等. 基于空間計(jì)量模型的耕地利用轉(zhuǎn)型對農(nóng)業(yè)經(jīng)濟(jì)增長影響研究[J]. 中國土地科學(xué)，2019，33(6)：53-61.

LU Xinhai, TANG Yifeng, YI Jialin, et al. Study on the impact of cultivated land use transition on agricultural economic growth based on spatial econometric model[J]. Chinese Land Science, 2019, 33(6): 53-61. (in Chinese with English abstract)

[23] 盧新海，楊喜，陳澤秀. 中國城市土地綠色利用效率測度及其時(shí)空演變特征[J]. 中國人口·資源與環(huán)境，2020，30(8)：83-91.

LU Xinhai, YANG Xi, CHEN Zexiu. Measurement and temporal-spatial evolution characteristics of urban land green use efficiency in China[J]. China Population, Resources and Environment, 2020, 30(8): 83-91. (in Chinese with English abstract)

[24] LU X H, KUANG B, LI J. Regional difference decomposition and policy implications of China’s urban land use efficiency under the environmental restriction[J]. Habitat International, 2018, 77: 32-39.

[25] 徐維祥，李露，周建平，等. 鄉(xiāng)村振興與新型城鎮(zhèn)化耦合協(xié)調(diào)的動(dòng)態(tài)演進(jìn)及其驅(qū)動(dòng)機(jī)制[J]. 自然資源學(xué)報(bào)，2020，35(9)：2044-2062.

XU Weixiang, LI Lu, ZHOU Jianping, et al. The dynamic evolution and its driving mechanism of coordination of rural rejuvenation and new urbanization[J]. Journal of Natural Resources, 2020, 35(9): 2044-2062. (in Chinese with English abstract)

[26] CAI J, LI X P, LIU L J, et al. Coupling and coordinated development of new urbanization and agro-ecological environment in China[J]. Science of the Total Environment, 2021, 776: 145837.

[27] LI W J, WANG Y, XIE S Y, et al. Coupling coordination analysis and spatiotemporal heterogeneity between urbanization and ecosystem health in Chongqing municipality, China[J]. Science of the Total Environment, 2021, 791: 148311.

[28] 朱慶瑩，陳銀蓉，胡偉艷，等. 中國土地集約利用與區(qū)域生態(tài)效率耦合協(xié)調(diào)度時(shí)空格局[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2020，36(4)：234-243.

ZHU Qingying, CHEN Yingrong, HU Weiyan, et al. Spatiotemporal pattern of coupling coordination degree between land intensive use and regional ecological efficiency in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(4): 234-243. (in Chinese with English abstract)

[29] 蓋美，張福祥. 遼寧省區(qū)域碳排放-經(jīng)濟(jì)發(fā)展-環(huán)境保護(hù)耦合協(xié)調(diào)分析[J]. 地理科學(xué)，2018，38(5)：764-772.

GAI Mei, ZHANG Fuxiang. Regional carbon emissions, economic development and environmental protection coupling in Liaoning province[J]. Scientia Geographica Sinica, 2018, 38(5): 764-772. (in Chinese with English abstract)

[30] 張文斌，張志斌，董建紅，等. 多尺度視角下耕地利用功能轉(zhuǎn)型及驅(qū)動(dòng)力分析：以甘肅省為例[J]. 地理科學(xué)，2021，41(5)：900-910.

ZHANG Wenbin, ZHANG Zhibin, DONG Jianhong, et al. Transformation and driving forces of cultivated land utilization function from a multi-scale perspective in Gansu Province[J]. Scientia Geographica Sinica, 2021, 41(5): 900-910. (in Chinese with English abstract)

[31] 衛(wèi)新東，林良國，羅平平，等. 耕地多功能耦合協(xié)調(diào)發(fā)展時(shí)空格局與驅(qū)動(dòng)力分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(4)：260-269.

WEI Xindong, LIN Liangguo, LUO Pingping, et al. Spatiotemporal pattern and driving force analysis of multi-functional coupling coordinated development of cultivated land[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(4): 260-269. (in Chinese with English abstract)

[32] 李怡，方斌，李裕瑞，等. 城鎮(zhèn)化進(jìn)程中耕地多功能權(quán)衡/協(xié)同關(guān)系演變及其驅(qū)動(dòng)機(jī)制[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(8)：244-254.

LI Yi, FANG Bin, LI Yurui, et al. Trade-off and synergy evolution of farmland functions and its dynamic mechanism in the process of urbanization[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 244-254.

[33] 張玥，代亞強(qiáng)，陳媛媛，等. 中國耕地低碳利用效率時(shí)空演變及其驅(qū)動(dòng)因素[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(8)：234-243.

ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, et al. Spatial-temporal evolution and driving factors of low-carbon use efficiency of cultivated land in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(8): 234-243. (in Chinese with English abstract)

[34] 王勁峰，徐成東. 地理探測器：原理與展望[J]. 地理學(xué)報(bào)，2017，72(1)：116-134.

WANG Jinfeng, XU Chengdong. Geodetector: Principle and prospective[J]. Acta Geographica Sinica, 2017, 72(1): 116-134. (in Chinese with English abstract)

[35] 匡兵，范翔宇，盧新海. 中國耕地利用綠色轉(zhuǎn)型效率的時(shí)空分異特征及其影響因素[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2021，37(21)：269-277.

KUANG Bing, FAN Xiangyu, LU Xinhai. Spatial-temporal differentiation characteristics of the efficiency of green transformation of cultivated land use and its affecting factors in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2021, 37(21): 269-277. (in Chinese with English abstract)

[36] 陳正發(fā)，史東梅，何偉，等. 1980-2015年云南坡耕地資源時(shí)空分布及演變特征分析[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2019，35(15)：256-265.

CHEN Zhengfa, SHI Dongmei, HE Wei, et al. Spatio-temporal distribution and evolution characteristics of slope farmland resources in Yunnan from 1980 to 2015[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2019, 35(15): 256-265. (in Chinese with English abstract)

[37] 李升發(fā)，李秀彬. 中國山區(qū)耕地利用邊際化表現(xiàn)及其機(jī)理[J]. 地理學(xué)報(bào)，2018，73(5)：803-817.

LI Shengfa, LI Xiubin. Economic characteristics and the mechanism of farmland marginalization in mountainous areas of China[J]. Acta Geographica Sinica, 2018, 73(5): 803-817. (in Chinese with English abstract)

[38] 付慧，劉艷軍，孫宏日，等. 京津冀地區(qū)耕地利用轉(zhuǎn)型時(shí)空分異及驅(qū)動(dòng)機(jī)制[J]. 地理科學(xué)進(jìn)展，2020，39(12)：1985-1998.

FU Hui, LIU Yanjun, SUN Hongri, et al. Spatiotemporal characteristics and dynamic mechanism of cultivated land use transition in the Beijing-Tianjin-Hebei region[J]. Progress in Geography, 2020, 39(12): 1985-1998. (in Chinese with English abstract)

[39] 秦彥杰，劉欣，趙艷霞，等. 1949年以來河北省耕地利用轉(zhuǎn)型特征及過程[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2022，38(12)：260-269.

QIN Yanjie, LIU Xin, ZHAO Yanxia, et al. Characteristics and process of the cultivated land use transformation in Hebei Province since 1949[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2022, 38(12): 260-269. (in Chinese with English abstract)

Spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development in China

ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, KE Xinli※

(,,430070,)

To explore the spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development, this paper applying the data of 30 provinces in China from 2000 to 2018 evaluated cultivated land production-living-ecological functions and revealed the spatial-temporal evolution of cultivated land multifunction with kernel density estimation and visualization mapping. Furthermore, this paper adopted coupling coordination degree model to analyze the cultivated land multifunctional coupling coordination development. Finally, this paper explored driving factors and driving mechanism of cultivated land multifunctional coupling coordination development by geographical detector. The results showed that: 1) from the perspective of temporal evolution, cultivated land production function level showed a fluctuation evolution trend of “rising-falling-rising”, the living function level showed an evolution trend of “rising-falling”, and the ecological function level increased slightly. All cultivated land functions had polarized characteristics. From the perspective of spatial pattern, the cultivated land production function level changed within [0.030, 0.608], showing a non-equilibrium spatial pattern of “high in the east, low in the west”. The living function level changed within [0.042, 0.672], showing a non-equilibrium spatial pattern of “high in the north, low in the south”. The ecological function level changed within [0.058, 0.897], showing a non-equilibrium spatial pattern of “high in the west, low in the east”. The areas with high level of production and living function were concentrated in the Northeast Plain, the Huang-Huai-Hai Plain and the Yangtze Region. The areas with low level of production function were located in the northwest and southwest regions, and the areas with low level of living function were located in the Beijing-Tianjin, Yunnan-Guizhou and southern regions in China. The area with high level of ecological function areas were located in the northeast and southwest regions, while the low-level areas were located in the Beijing-Tianjin-Hebei region, the Yangtze River Delta region and the central region of the Yangtze River. 2) The cultivated land multifunctional coupling coordination level changed within [0.093, 0.554], and the overall level had steadily improved but was still low. The cultivated land multifunctional coupling coordination levels of most provinces were still out of balance. The spatial pattern of cultivated land multifunctional coupling coordination level showed an evolution trend of agglomeration and promotion. In detail, the cultivated land multifunctional coupling coordination level of the northeast China, Beijing-Tianjin-Hebei region and the Yangtze River Delta were relatively high, while that was low in the northwestern part of China. 3) The evolution of cultivated land multifunctional coupling coordination level was significantly affected by the business entities conditions, agricultural development level, industrial development level and urbanization level. In short, the driving force of the business entities conditions and urbanization level were more significant than that of the agricultural development level and industrial development level. Moreover, the interaction type of two driving factors was two-factor enhancement or nonlinear enhancement, which indicated that the driving force of various factors had obvious divergence and interaction of factors positively strengthened the single-factor driving force, indicated that the divergence and interaction of various driving factors positively strengthened the single-factor driving force. In conclusion, the research results provide the guiding effect for deepening the study on cultivated land multifunction and decision support for promoting the differential governance of cultivated land resources and mutual promotion of cultivated land functions.

land use; cultivated land; multifunction; coupling coordination degree; spatial-temporal evolution; driving factors

2022-09-23

2023-02-10

國家自然科學(xué)基金面上項(xiàng)目（41971240）；國家社科基金后期資助項(xiàng)目（19FGLB071）；中央高?；究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金資助項(xiàng)目（2662022XCZX012）

張玥，博士生，研究方向?yàn)橥恋刭Y源管理。Email：yuez@webmail.hzau.edu.cn

柯新利，博士，教授，博士生導(dǎo)師，研究方向?yàn)橥恋乩脙?yōu)化。Email：kexl@mail.hzau.edu.cn

10.11975/j.issn.1002-6819.202209185

S28; F301.21

A

1002-6819(2023)-07-0244-12

張玥，代亞強(qiáng)，陳媛媛，等. 中國耕地多功能耦合協(xié)調(diào)時(shí)空演變及其驅(qū)動(dòng)因素[J]. 農(nóng)業(yè)工程學(xué)報(bào)，2023，39(7)：244-255. doi：10.11975/j.issn.1002-6819.202209185 http://www.tcsae.org

ZHANG Yue, DAI Yaqiang, CHEN Yuanyuan, et al. Spatial-temporal evolution and driving factors of cultivated land multifunctional coupling coordination development in China[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2023, 39(7): 244-255. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.202209185 http://www.tcsae.org