蒲金芳 王亞楠 張蚌蚌 盧龍輝 張鳳榮 王 數(shù)*

(1.中國農(nóng)業(yè)大學 土地科學與技術學院，北京 100193;2.自然資源部農(nóng)用地質(zhì)量與監(jiān)控重點實驗室，北京 100193;3.西北農(nóng)林科技大學 經(jīng)濟管理學院,陜西 楊凌 712100)

隨著我國城市化、工業(yè)化和信息化進程的不斷加快，我國的城鎮(zhèn)化率從1978年的17.92%上升至2018年的59.58%。城鄉(xiāng)發(fā)展逐漸形成了“以鄉(xiāng)養(yǎng)城”的模式。人口和資源不斷由鄉(xiāng)村涌向城市，城鄉(xiāng)發(fā)展差距不斷擴大，城鄉(xiāng)關系失調(diào)現(xiàn)象顯著，導致鄉(xiāng)村出現(xiàn)耕地撂荒、環(huán)境污染、村莊“空心化”等問題[1-3]。如何科學解決鄉(xiāng)村發(fā)展過程中所出現(xiàn)的一系列問題已成為當前首要任務。鄉(xiāng)村發(fā)展是指鄉(xiāng)村地域經(jīng)濟、社會、聚落等特征不斷演化的過程[4]。土地作為鄉(xiāng)村發(fā)展的基本要素，是人類進行各項社會活動的載體[5]，鄉(xiāng)村發(fā)展過程中存在的問題可以在土地利用中得以反映[6]，土地利用與鄉(xiāng)村發(fā)展密切相關。由此，從土地利用與鄉(xiāng)村發(fā)展關系的角度解決鄉(xiāng)村問題成為鄉(xiāng)村振興的重要途徑。土地利用多功能是指不同土地利用類型相互聯(lián)系與作用而形成的多種功能,體現(xiàn)了土地利用滿足人類多種需求的程度和能力，是研究土地利用的重要視角。提升土地利用多功能是鄉(xiāng)村發(fā)展的有效路徑[7]，而鄉(xiāng)村發(fā)展也驅動著土地利用多功能演變，二者相互影響，共同促進[6]，實現(xiàn)土地利用多功能與鄉(xiāng)村發(fā)展協(xié)調(diào)是鄉(xiāng)村振興的重要內(nèi)容。因此，探索土地利用多功能與鄉(xiāng)村發(fā)展之間的關系，有助于鄉(xiāng)村人地關系協(xié)調(diào)發(fā)展，對于解決鄉(xiāng)村發(fā)展問題與促進鄉(xiāng)村可持續(xù)發(fā)展具有重要意義。

在土地利用多功能的研究中，多將土地利用多功能劃分為社會功能、經(jīng)濟功能和生態(tài)(環(huán)境)功能[8]，并開展時空分異研究，從自然、社會、經(jīng)濟和環(huán)境等多方面探析土地利用多功能的影響因素[9-10]，提出存在問題及制定政策建議等[11]。在鄉(xiāng)村發(fā)展的研究中主要包括定性和定量的研究：從定性的角度[12-13]，探析我國鄉(xiāng)村發(fā)展的歷史脈絡，同時與西方國家鄉(xiāng)村發(fā)展歷程對比，分析我國鄉(xiāng)村發(fā)展存在的問題；從定量的角度[14-16]，基于不同理論基礎構建鄉(xiāng)村發(fā)展測度指標體系，分析空間格局及影響因素，劃分鄉(xiāng)村發(fā)展類型，提出不同的鄉(xiāng)村發(fā)展模式。目前，針對土地利用與鄉(xiāng)村發(fā)展關系的研究主要是探討耕地[17]、林地[18]和農(nóng)村居民點用地[19]等用地類型變化對鄉(xiāng)村發(fā)展的影響，以及農(nóng)村土地制度對鄉(xiāng)村發(fā)展進程的影響[20]，土地利用與鄉(xiāng)村發(fā)展的交互作用研究主要集中在城鄉(xiāng)關系失調(diào)導致鄉(xiāng)村衰落[1]、鄉(xiāng)村土地市場對鄉(xiāng)村發(fā)展進程影響[21-22]和農(nóng)地流轉與農(nóng)村發(fā)展轉型的關系[23]，針對土地利用轉型與鄉(xiāng)村發(fā)展轉型也開展了系統(tǒng)地研究[6,24-25]。但是，對于土地利用與鄉(xiāng)村發(fā)展研究還存在以下不足：一是土地利用系統(tǒng)與鄉(xiāng)村發(fā)展系統(tǒng)之間的交互機理有待進一步厘清；二是基于土地功能視角探討土地利用與鄉(xiāng)村發(fā)展二者之間的耦合關系的研究較少。

江西省萬年縣是稻作文化之源，鄱陽湖地區(qū)綠色有機食品基地，也是承接東部發(fā)達地區(qū)產(chǎn)業(yè)轉移的重要區(qū)域。自2015年縣政府加快推進機械電子、紡織新材料、食品藥品三大主導產(chǎn)業(yè)集群集約高效發(fā)展，并針對村莊逐步開展“秀美鄉(xiāng)村”基礎設施建設項目以來，區(qū)域發(fā)展呈現(xiàn)多產(chǎn)融合的發(fā)展態(tài)勢，鄉(xiāng)村社會經(jīng)濟發(fā)展較快。萬年縣屬于典型的山地丘陵區(qū)，自然資源稟賦差異較大，其土地利用多功能與鄉(xiāng)村發(fā)展不協(xié)調(diào)的問題顯著。因此，本研究擬以江西省萬年縣為研究區(qū)域，2010年和2018年為研究時點，村域為研究單元，在解析土地利用多功能與鄉(xiāng)村發(fā)展交互機理的基礎上，構建土地利用多功能與鄉(xiāng)村發(fā)展評價指標體系，測度土地利用多功能與鄉(xiāng)村發(fā)展水平，借助耦合協(xié)調(diào)度模型，分析萬年縣鄉(xiāng)村土地利用多功能與鄉(xiāng)村發(fā)展的耦合協(xié)調(diào)關系，根據(jù)耦合作用關系劃分鄉(xiāng)村類型，以期為村莊規(guī)劃和鄉(xiāng)村振興戰(zhàn)略實施提供基礎支撐。

1 研究方法

1.1 土地利用多功能與鄉(xiāng)村發(fā)展交互機理

在鄉(xiāng)村發(fā)展轉型過程中，土地利用與鄉(xiāng)村發(fā)展關系密切。土地的不同功能的組合形式屬于土地利用的隱形形態(tài)[26]，是土地利用提供商品和服務及滿足人類社會需求的程度[27]，與鄉(xiāng)村發(fā)展共同刻畫鄉(xiāng)村人地關系地域格局。土地利用多功能主要體現(xiàn)在社會功能、經(jīng)濟功能和生態(tài)功能3個方面[11，28-29]，其區(qū)域土地系統(tǒng)的支撐和保障能力強，食物生產(chǎn)和資源維持水平高，提供生態(tài)環(huán)境質(zhì)量好，越容易得到政府政策幫扶和具備產(chǎn)業(yè)投資條件，從而推動鄉(xiāng)村高水平發(fā)展。而鄉(xiāng)村發(fā)展主要體現(xiàn)在人-地-業(yè)等方面[16，30]，區(qū)域資源稟賦豐富、產(chǎn)業(yè)發(fā)展?jié)摿Υ?、人居環(huán)境和交通區(qū)位優(yōu)越必然驅使人類提高土地利用效率和調(diào)整生產(chǎn)要素結構，提升土地利用多功能水平。同時，二者都受到自然過程和人類活動的作用，而自然過程和人類活動也因土地利用多功能和鄉(xiāng)村發(fā)展的變化而得到反饋。

農(nóng)民是鄉(xiāng)村發(fā)展的主體，改善民生是鄉(xiāng)村發(fā)展的主要目的，提高農(nóng)民福祉是鄉(xiāng)村發(fā)展的深刻內(nèi)涵[31]。而土地利用系統(tǒng)的核心是土地利用，土地利用多功能作為土地利用的全面表達，能較好揭示土地利用的結構和狀態(tài)[6]。因此，土地利用多功能與鄉(xiāng)村發(fā)展的關鍵是土地利用與農(nóng)民福祉的關系，深刻地反映著鄉(xiāng)村人地系統(tǒng)狀態(tài)。如圖1所示，合理的土地利用可以協(xié)調(diào)土地利用系統(tǒng)要素，提高土地效益，從而增進農(nóng)民福祉；而農(nóng)民福祉的提高會激發(fā)農(nóng)民積極性，增加對土地要素的投入，優(yōu)化土地利用結構，提高土地利用效率，從而實現(xiàn)土地可持續(xù)利用。

圖1 土地利用多功能與鄉(xiāng)村發(fā)展交互機理示意圖Fig.1 Diagram of the interactive relationship between land use multi-functions and rural development

土地利用多功能與鄉(xiāng)村發(fā)展交互作用是一個動態(tài)協(xié)調(diào)的過程，二者必然存在正向推動和負向制約的作用，有土地利用多功能與鄉(xiāng)村發(fā)展均為正向推動作用、土地利用多功能正向推動鄉(xiāng)村發(fā)展負向制約或鄉(xiāng)村發(fā)展正向推動土地利用多功能負向抑制、土地利用多功能與鄉(xiāng)村發(fā)展均為負向抑制作用3種情況。因此，將土地利用多功能與鄉(xiāng)村發(fā)展關系分為3類，包括相互促進型，促進與阻礙并存型和相互阻礙型。土地利用多功能與鄉(xiāng)村發(fā)展相互促進是指區(qū)域土地利用系統(tǒng)與鄉(xiāng)村發(fā)展系統(tǒng)具有較好協(xié)調(diào)性，在區(qū)域發(fā)展過程中系統(tǒng)之間產(chǎn)生協(xié)同作用，二者形成相互促進的關系，鄉(xiāng)村人地系統(tǒng)呈正向演化；土地利用多功能與鄉(xiāng)村發(fā)展促進與阻礙并存是指區(qū)域土地利用系統(tǒng)與鄉(xiāng)村發(fā)展系統(tǒng)之間處于失衡狀態(tài)，由土地利用系統(tǒng)和鄉(xiāng)村發(fā)展系統(tǒng)構成的鄉(xiāng)村人地系統(tǒng)處于劇烈變化期，二者交互作用中同時存在推動和制約關系，鄉(xiāng)村人地系統(tǒng)演化方向呈現(xiàn)多向性；土地利用多功能與鄉(xiāng)村發(fā)展相互阻礙是指區(qū)域土地利用系統(tǒng)與鄉(xiāng)村發(fā)展系統(tǒng)中各要素處于無序狀態(tài)，要素之間耦合協(xié)調(diào)水平較低，二者之間形成相互制約的關系。

1.2 土地利用多功能與鄉(xiāng)村發(fā)展評價模型

1.2.1指標標準化與權重確定

由于各個評價指標量綱、數(shù)量級及指標性質(zhì)存在差異，因此采用極差標準化對所有指標數(shù)據(jù)進行標準化處理，其計算公式如下：

正向指標：

(1)

負向指標：

(2)

為了避免主觀權重賦值法的不確定性，客觀合理確定指標權重。因此，采用熵權法確定指標權重，具體計算過程參見文獻[32]，具體結果見表1和表2。

1.2.2土地利用多功能評價指標體系

在已有研究基礎上[9-10，32]，結合研究區(qū)特點，從社會功能、經(jīng)濟功能和生態(tài)功能3個維度，選取9個指標，構建土地利用多功能評價指標體系(表1)。1)社會功能選取人均宅基地面積、人口密度和人均耕地面積3個指標，其中人均宅基地面積和人口密度表征空間承載能力，從人口承載資源條件和規(guī)模方面反映土地系統(tǒng)支撐功能，其值越大，表示土地系統(tǒng)人口承載能力越強；人均耕地面積表征土地生產(chǎn)條件所能提供的基本生活保障功能，其值越大，表示土地系統(tǒng)資源保障能力越強；2)經(jīng)濟功能選取農(nóng)林牧漁業(yè)總產(chǎn)值、土地墾殖率和集體經(jīng)營性建設用地面積3個指標，其中農(nóng)林牧漁業(yè)總產(chǎn)值和土地墾殖率分別從生產(chǎn)能力與資源條件方面綜合表征農(nóng)業(yè)生產(chǎn)能力，其值越大，表示農(nóng)業(yè)經(jīng)濟發(fā)展水平越高；而農(nóng)村集體經(jīng)營性建設用地面積可從規(guī)模和狀態(tài)方面表征鄉(xiāng)村經(jīng)濟發(fā)展狀況[33]，其值越大，表示非農(nóng)經(jīng)濟發(fā)展?jié)摿υ酱螅?)生態(tài)功能選用森林覆蓋率、水網(wǎng)密度指數(shù)和生態(tài)系統(tǒng)服務價值3個指標，森林覆蓋率越高、水網(wǎng)密度指數(shù)越大，表示鄉(xiāng)村維持區(qū)域生態(tài)安全能力和生態(tài)保育功能越強；生態(tài)系統(tǒng)服務價值表征生態(tài)服務功能，其值越大，表示鄉(xiāng)村生態(tài)系統(tǒng)服務價值越高，生態(tài)服務功能越強。

表1 土地利用多功能評價指標體系Table 1 Measurement index system of land use multi-functions

1.2.3土地利用多功能指數(shù)測算

土地利用多功能指數(shù)用于測度研究區(qū)土地利用多功能水平大小，其值越高代表土地利用多功能水平越高，反之亦同。其計算公式如下：

(3)

式中：LUMi為第i個鄉(xiāng)村土地利用多功能測度值；Lij為第i個鄉(xiāng)村的第j個指標的標準化值；Wij為第i個鄉(xiāng)村的第j個指標的權重值；n為指標個數(shù)。

1.2.4鄉(xiāng)村發(fā)展水平評價指標體系

參考已有研究[14]和[16]，結合指標科學性和數(shù)據(jù)的可得性，分別從資源稟賦、產(chǎn)業(yè)發(fā)展、人居環(huán)境和交通區(qū)位4個維度，選取10個指標，測度鄉(xiāng)村發(fā)展水平(表2)。具體考量如下：1)資源稟賦選取主要農(nóng)業(yè)資源面積和主要農(nóng)作物年產(chǎn)量2個指標，農(nóng)業(yè)資源面積從區(qū)域資源本底方面表征村域農(nóng)業(yè)資源稟賦；主要農(nóng)作物年產(chǎn)量是反映不同鄉(xiāng)村區(qū)域農(nóng)業(yè)生產(chǎn)水平，其值越大，表示鄉(xiāng)村資源稟賦狀況越好；2)產(chǎn)業(yè)發(fā)展維度，鑒于數(shù)據(jù)的可獲取性，結合實地調(diào)查結果，選取土地開發(fā)強度、工礦用地面積和鄉(xiāng)村從業(yè)人員占比3個指標能夠反映產(chǎn)業(yè)發(fā)展狀況，其中工礦用地面積反映工礦產(chǎn)業(yè)規(guī)模，其值越大，表示村域的產(chǎn)業(yè)發(fā)展?jié)摿υ酱?；鄉(xiāng)村從業(yè)人員比例和土地開發(fā)強度從側面表征鄉(xiāng)村產(chǎn)業(yè)發(fā)展狀況，鄉(xiāng)村從業(yè)人員比例越高，土地開發(fā)強度越大，表示鄉(xiāng)村產(chǎn)業(yè)發(fā)展情況越好；3)人居環(huán)境選取自來水入戶率和聚落聚集度2個指標，自來水入戶率表征區(qū)域居住保障狀況，居住保障狀況越好，鄉(xiāng)村發(fā)展水平越高；聚落聚集度反映城市化過程中鄉(xiāng)村聚落空間的擴展與演化，集聚度越高，鄉(xiāng)村發(fā)展能力越弱[16]；4)交通區(qū)位選取公路密度、地形起伏度和距縣政府駐地距離3個指標，公路密度和地形起伏度表征區(qū)域交通便利程度，公路密度越大、地形起伏度越小，交通便利程度越高，越有助于鄉(xiāng)村發(fā)展；距縣政府駐地距離表征獲得外部資源的能力，其區(qū)位優(yōu)勢越好，獲得發(fā)展機會可能性越大。

表2 鄉(xiāng)村發(fā)展水平評價指標體系Table 2 Measurement index system of rural development level

1.2.5鄉(xiāng)村發(fā)展水平指數(shù)測算

鄉(xiāng)村發(fā)展指數(shù)用于測度研究區(qū)鄉(xiāng)村發(fā)展水平情況，其值越高代表鄉(xiāng)村發(fā)展狀況越好，反之亦同，其計算公式如下：

(4)

式中：RDi為第i個鄉(xiāng)村發(fā)展水平的測度值；Rij為第i個鄉(xiāng)村的第j個指標的標準化值；Wij為第i個鄉(xiāng)村的第j個指標的權重值；n為指標個數(shù)。

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

1.3.1耦合度模型

耦合度為物理學概念，表示2個或2個以上的系統(tǒng)或運動方式之間相互影響、相互作用的關系，現(xiàn)已被廣泛應用于地理學研究中[36-37]。借助容量耦合概念構建耦合模型，以準確識別土地利用多功能與鄉(xiāng)村發(fā)展的耦合關系，把握鄉(xiāng)村規(guī)律。耦合度模型如下：

(5)

式中：C為耦合度，C∈[0,1]，值越大系統(tǒng)間相互作用越強；L(x)為土地利用多功能綜合值；R(y)為鄉(xiāng)村發(fā)展綜合值。在耦合度模型計算中，鑒于土地利用多功能與鄉(xiāng)村發(fā)展水平測度結果分布區(qū)間存在差異，故利用式(1)對測度結果進行標準化處理，對于值為0的數(shù)據(jù)采用除其以外的最小值代替。參照蔡進等研究[36]，結合研究區(qū)實際情況，將耦合度劃分為4個等級，即低度耦合(0≤C≤0.3)、中度耦合(0.3

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

耦合協(xié)調(diào)度模型指相互影響系統(tǒng)之間的良性關系的耦合程度的大小，體現(xiàn)系統(tǒng)之間耦合狀態(tài)的程度。而土地利用多功能與鄉(xiāng)村發(fā)展之間存在相互影響、相互作用的程度正好反映了二者之間的協(xié)調(diào)程度。由此，構建耦合協(xié)調(diào)度模型如下：

T=α×L(x)+β×R(y)

(6)

(7)

式中：T為土地利用多功能系統(tǒng)和鄉(xiāng)村發(fā)展系統(tǒng)間協(xié)調(diào)指數(shù)；α和β為待定指數(shù)且滿足α+β=1，R為耦合協(xié)調(diào)度。本研究認為土地利用多功能和鄉(xiāng)村發(fā)展同等重要，故α=β=0.5。參照劉愿理等研究[37]，并結合區(qū)域實際情況，將耦合協(xié)調(diào)度劃分為7個等級，并結合兩者系統(tǒng)得分差值劃分為同步和失衡兩類，分別用S和L代表，組合耦合協(xié)調(diào)度和系統(tǒng)得分差值類型，得到耦合協(xié)調(diào)發(fā)展類型14類。并在此基礎上，梳理土地利用多功能與鄉(xiāng)村發(fā)展作用關系，包括相互促進、相互阻礙和促進與阻礙并存3種類型(表3)。

表3 土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)類型劃分Table 3 Division of coordination types between land use multi-functions and rural development

1.4 空間自相關

本研究采用全局莫蘭指數(shù)(Global-Moran’s I)和局部莫蘭指數(shù)(Local-Moran’s I)來分析全局空間自相關和局部空間自相關，用以表征土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)度的空間集聚特征[36]，公式如下：

(8)

(9)

2 研究區(qū)概況與數(shù)據(jù)來源

2.1 研究區(qū)概況

萬年縣位于江西省東北部、鄱陽湖東南岸，地理位置為116°46′41″～117°15′16″ E、28°30′00″～28°54′08″ N，隸屬于上饒市，其北連鄱陽縣、樂平市，東接戈陽縣、貴溪市，南毗余江縣，西鄰余干縣，總面積1 140.76 km2。區(qū)域地貌類型復雜，根據(jù)地貌成因類型及其形態(tài)特征，可劃分低山、丘陵和平原，土地利用類型以耕地、林地和水域為主。2010年全縣三次產(chǎn)業(yè)結構16.9∶56.1∶27，農(nóng)村居民人均可支配收入4 967元；2018年全縣三次產(chǎn)業(yè)結構10.6：57.3：32.1，農(nóng)村居民人均可支配收入13 659元，區(qū)域社會經(jīng)濟發(fā)展較快。全縣共轄12個鄉(xiāng)鎮(zhèn)，考慮到國營農(nóng)場、墾殖場、水庫管理局、城鎮(zhèn)核心區(qū)等地域的特殊性及數(shù)據(jù)缺乏的現(xiàn)實狀況，最終得到128個村級評價單元。

圖2 江西省萬年縣區(qū)位圖Fig.2 Location of Wannian County in Jiangxi Province

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

本研究選擇2010年和2018年為研究時點，以村域為研究對象，所使用的數(shù)據(jù)主要包括： 1)地理空間數(shù)據(jù)。主要包括2010年和2018年萬年縣土地利用現(xiàn)狀圖、三權宗地數(shù)據(jù)、行政村邊界等，來源于萬年縣自然資源局；30 m分辨率DEM數(shù)據(jù)，來源于地理空間數(shù)據(jù)云(http: ∥www.gscloud.cn/)平臺；宅基地面積和集體經(jīng)營性建設用地面積是在土地用調(diào)查數(shù)據(jù)和三權宗地數(shù)據(jù)的基礎上，通過外業(yè)核查方式獲?。?)社會經(jīng)濟數(shù)據(jù)。主要包括統(tǒng)計數(shù)據(jù)和調(diào)研數(shù)據(jù)，統(tǒng)計數(shù)據(jù)來源于《萬年縣統(tǒng)計年鑒(2011)》和《萬年縣統(tǒng)計年鑒(2019)》；調(diào)研數(shù)據(jù)是課題組與12個鎮(zhèn)政府、當?shù)卮逦瘯途用襁M行座談，收集資料和填寫調(diào)查問卷整理得到，包括研究時點村級的鄉(xiāng)村人口、農(nóng)作物產(chǎn)量、水產(chǎn)品產(chǎn)量、鄉(xiāng)村從業(yè)人員、通自來水戶數(shù)等。

3 結果與分析

3.1 土地利用多功能與鄉(xiāng)村發(fā)展測度結果

3.1.1土地利用多功能測度結果

在時間上，2010年萬年縣土地利用多功能指數(shù)為0.19～0.32，均值為0.23；至2018年，土地利用多功能指數(shù)為0.21～0.43，均值為0.30，平均增長26.60%。其中社會功能增長7.37%，經(jīng)濟功能增長115.28%，生態(tài)功能下降7.38%。經(jīng)濟功能增長幅度較大，主要是區(qū)域進行農(nóng)用地流轉，近80%耕地流轉給種糧大戶進行規(guī)?；?jīng)營，提高土地產(chǎn)出，同時產(chǎn)業(yè)用地呈現(xiàn)擴張態(tài)勢，促使土地經(jīng)濟功能提升較快，而生態(tài)功能下降主要是建設用地快速擴張，導致生態(tài)用地劇烈減少，生態(tài)功能下降。在空間上，2010年土地利用多功能呈現(xiàn)多核心的格局特征，主要以較低水平和中等水平鄉(xiāng)村為主，所占比例為57.03%，高水平鄉(xiāng)村主要分布于梓埠鎮(zhèn)、湖云鄉(xiāng)、陳營鎮(zhèn)和裴梅鎮(zhèn)；至2018年，土地利用多功能整體水平提高，空間集聚態(tài)勢明顯，土地利用多功能仍以較低水平和中等水平鄉(xiāng)村為主，所占比例為59.37%，高水平鄉(xiāng)村主要分布于城區(qū)周邊鄉(xiāng)鎮(zhèn)(圖3)。

圖3 研究區(qū)2010—2018年土地利用多功能空間格局Fig.3 Spatial pattern of land use muti-functions in the study area from 2010 to 2018

3.1.2鄉(xiāng)村發(fā)展測度結果

在時間上，2010年萬年縣鄉(xiāng)村發(fā)展指數(shù)為0.15～0.47，均值為0.30；至2018年，鄉(xiāng)村發(fā)展指數(shù)為0.18～0.64，均值為0.42，平均增長40.60%。其中資源稟賦提升110.06%，產(chǎn)業(yè)發(fā)展下降13.31%，人居環(huán)境提升70.24%，交通區(qū)位提升24.01%。產(chǎn)業(yè)發(fā)展下降主要是政府征收土地，發(fā)展雷竹、光伏產(chǎn)業(yè)和開展礦業(yè)整治，區(qū)域吸納就業(yè)能力下降，導致鄉(xiāng)村人口外流，造成鄉(xiāng)村產(chǎn)業(yè)水平有所下降。在空間上，2010年鄉(xiāng)村發(fā)展呈現(xiàn)以城區(qū)周邊鄉(xiāng)鎮(zhèn)為中心向四周遞減的格局特征。鄉(xiāng)村發(fā)展水平主要以中等水平和較高水平為主，所占比例為60.16%。高水平鄉(xiāng)村位于陳營鎮(zhèn)，低水平鄉(xiāng)村主要位于西北部

和東部邊緣區(qū)，包括梓埠鎮(zhèn)、湖云鄉(xiāng)、大源鎮(zhèn)和裴梅鎮(zhèn)；至2018年，鄉(xiāng)村發(fā)展水平呈現(xiàn)以城區(qū)為中心的圈層分布特征突出(圖4)。2010—2018年間，鄉(xiāng)村發(fā)展水平空間格局集聚態(tài)勢明顯，較高水平鄉(xiāng)村逐漸向城區(qū)周邊靠攏。這主要是在區(qū)域發(fā)展過程中，城區(qū)城鎮(zhèn)化輻射作用逐步增強，城鎮(zhèn)周邊鄉(xiāng)村發(fā)展外部環(huán)境較好，區(qū)域形成以點帶面的發(fā)展格局。

3.2 土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)關系

3.2.1耦合度測度結果

2010—2018年萬年縣土地利用多功能與鄉(xiāng)村發(fā)展耦合度增長較快，耦合類型主要以中度耦合為主。2010—2018年，耦合度由0.43提升至0.45，增長4.68%。根據(jù)耦合度類型劃分，2010年耦合度類型以中度耦合為主，總計110個鄉(xiāng)村，所占比例為85.94%；2018年耦合度類型仍以中度耦合為主，總計92個鄉(xiāng)村，所占比例為71.87%。但相較于2010年，低度耦合鄉(xiāng)村數(shù)量由4個上升至7個；高度耦合鄉(xiāng)村數(shù)量由14個上升至29個；無極度耦合鄉(xiāng)村，需進一步加強二者的耦合關系。在空間上，2010—2018年土地利用多功能與鄉(xiāng)村發(fā)展耦合度呈現(xiàn)明顯的擴散趨勢，高度耦合鄉(xiāng)村呈現(xiàn)向城區(qū)周邊聚集為主，并向四周擴散的格局特征；低度耦合鄉(xiāng)村呈現(xiàn)零散化特征(圖5)。表明區(qū)域整體耦合度提高，但仍存在局部落后的情況。

圖5 研究區(qū)2010—2018年土地利用多功能與鄉(xiāng)村發(fā)展耦合度空間格局Fig.5 Spatial pattern of coupling level between land use muti-functions and rural development in the study area from 2010 to 2018

3.2.2耦合協(xié)調(diào)度測度結果

2010—2018年萬年縣土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)度由0.43提升至0.45，增長3.88%。2010年耦合協(xié)調(diào)類型以勉強協(xié)調(diào)為主，總計72個鄉(xiāng)村，所占比例為56.25%；2018年耦合協(xié)調(diào)類型仍以勉強協(xié)調(diào)為主，總計65個鄉(xiāng)村，所占比例為50.78%。相較于2010年，失調(diào)類型鄉(xiāng)村主要以輕度失調(diào)為主，下降31.16%；協(xié)調(diào)類型鄉(xiāng)村主要以勉強協(xié)調(diào)為主，而協(xié)調(diào)度上升幅度最大的鄉(xiāng)村類型為輕度協(xié)調(diào)，上升86.66%。在空間上， 2010—2018年失調(diào)類型鄉(xiāng)村呈現(xiàn)收縮的空間格局特征，其收縮區(qū)域主要位于湖云鄉(xiāng)、大源鎮(zhèn)、青云鄉(xiāng)和裴梅鎮(zhèn)(圖6)。2010—2018年協(xié)調(diào)類型鄉(xiāng)村主要以中度協(xié)調(diào)和輕度協(xié)調(diào)擴張為主，中度協(xié)調(diào)類型鄉(xiāng)村主要分布于陳營鎮(zhèn)、上坊鄉(xiāng)和裴梅鎮(zhèn)，呈現(xiàn)多核心分布；輕度協(xié)調(diào)鄉(xiāng)村呈現(xiàn)局地集聚、整體分散的空間格局特征。結合土地利用多功能與鄉(xiāng)村發(fā)展的耦合度和協(xié)調(diào)度的空間分布，二者具有較強的空間關聯(lián)性。

圖6 研究區(qū)2010—2018年土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)度空間格局Fig.6 Spatial pattern of coupling coordination degree between land use muti-functions and rural development in the study area from 2010 to 2018

3.2.3空間關聯(lián)格局分析

2010年耦合度和協(xié)調(diào)度全局莫蘭指數(shù)分別為0.266 1和0.160 5；2018年耦合度和協(xié)調(diào)度全局莫蘭指數(shù)分別為0.272 5和0.224 8，表明2010—2018年萬年縣土地利用多功能與鄉(xiāng)村發(fā)展的耦合度和協(xié)調(diào)度具有顯著的空間正相關關系，且土地利用多功能和鄉(xiāng)村發(fā)展耦合水平和協(xié)調(diào)度差異在縮小，空間集聚效應增強。從耦合度看，高-高集聚區(qū)主由陳營鎮(zhèn)轉向上坊鄉(xiāng)；低-低集聚區(qū)主要位于萬年縣東部的湖云鄉(xiāng)、梓埠鎮(zhèn)、大源鎮(zhèn)和裴梅鎮(zhèn)，呈收縮態(tài)勢；高-低集聚和低-高集聚區(qū)轉變?yōu)榱阈欠植肌１砻魅f年縣耦合度有顯著的集聚效應，耦合水平局部空間差異減小。從協(xié)調(diào)度看，耦合協(xié)調(diào)度集聚空間分布特征與耦合度相似(圖7)。但土地利用多功能與鄉(xiāng)村發(fā)展耦合協(xié)調(diào)水平在局部空間差異性減小顯著，低-低集聚區(qū)在梓埠鎮(zhèn)有所增加，區(qū)域協(xié)調(diào)水平有待提高。

3.2.4耦合關系類型

參照表3的分類標準，將鄉(xiāng)村類型分為土地利用多功能與鄉(xiāng)村發(fā)展相互阻礙型、土地利用多功能與鄉(xiāng)村發(fā)展相互促進型和土地利用多功能與鄉(xiāng)村發(fā)展促進和阻礙并存型3大類12小類，并繪制空間分布圖(圖8)。

圖8 2010—2018年鄉(xiāng)村耦合協(xié)調(diào)類型圖Fig.8 Diagram of the rural coupling coordination types from 2010 to 2018

1)土地利用多功能與鄉(xiāng)村發(fā)展相互促進型。2010—2018年萬年縣相互促進型鄉(xiāng)村數(shù)量呈現(xiàn)減少趨勢，由71個下降至15個，空間上呈現(xiàn)西縮東擴的格局特征，其中土地利用多功能和鄉(xiāng)村發(fā)展水平均以中等水平或較高水平為主，“雙高”型鄉(xiāng)村較少。表明區(qū)域土地利用系統(tǒng)和鄉(xiāng)村發(fā)展系統(tǒng)處于失衡狀態(tài)，二者耦合協(xié)調(diào)水平有待進一步加強。

2)土地利用多功能與鄉(xiāng)村發(fā)展促進與阻礙并存型。萬年縣促進與阻礙并存型鄉(xiāng)村數(shù)量呈現(xiàn)增加趨勢，由24個增至101個，成為2018年土地利用多功能與鄉(xiāng)村發(fā)展耦合關系的主導類型，且主要以土地利用多功能滯后和勉強協(xié)調(diào)為主。表明萬年縣“秀美鄉(xiāng)村”建設、村莊人居環(huán)境改造等工程的效果顯著，區(qū)域鄉(xiāng)村發(fā)展狀況得到明顯提升。

3)土地利用多功能與鄉(xiāng)村發(fā)展相互阻礙型。萬年縣相互阻礙型鄉(xiāng)村數(shù)量呈現(xiàn)急劇減少的趨勢，由33個下降至12個，耦合協(xié)調(diào)類型以輕度失調(diào)為主，在空間上呈現(xiàn)邊緣化趨勢，主要分布于萬年縣西北部。其主要是西北部鄱陽湖平原的梓埠鎮(zhèn)以傳統(tǒng)農(nóng)業(yè)為主，交通條件較差，致使土地利用多功能和鄉(xiāng)村發(fā)展水平較低，二者之間形成相互制約作用。

4 結論與建議

4.1 結論

1)2010—2018年萬年縣鄉(xiāng)村發(fā)展水平增速高于土地利用多功能，空間格局特征各異。2010—2018年萬年縣土地利用多功能由0.23提升至0.30，增長26.60%，鄉(xiāng)村發(fā)展水平由0.30提升至0.42，增長40.60%；空間上，土地利用多功能呈現(xiàn)多核心分布格局特征，并呈現(xiàn)局地集聚的趨勢，鄉(xiāng)村發(fā)展呈現(xiàn)以城區(qū)為中心集聚態(tài)勢，并向四周逐漸遞減的格局特征。

2)萬年縣土地利用多功能與鄉(xiāng)村發(fā)展耦合關系明顯，耦合度和協(xié)調(diào)度隨時間推移發(fā)生演變。2010—2018年萬年縣土地利用多功能與鄉(xiāng)村發(fā)展耦合度和協(xié)調(diào)度分別以4.68%和3.88%速度增長，耦合類型以中度耦合為主，耦合協(xié)調(diào)類型以勉強協(xié)調(diào)為主，二者在空間上具有一定關聯(lián)性。

3)2010—2018年土地利用多功能與鄉(xiāng)村發(fā)展水平耦合協(xié)調(diào)發(fā)展主要類型以相互阻礙類型轉變?yōu)榇龠M與阻礙并存型，相互促進型和相互阻礙型均呈下降趨勢，這與土地利用多功能增速低于鄉(xiāng)村發(fā)展水平有關。

4.2 建議

土地利用多功能與鄉(xiāng)村發(fā)展耦合關系研究是揭示鄉(xiāng)村地域人地關系的重要視角，也為鄉(xiāng)村發(fā)展分類施策提供新的路徑。從二者耦合的角度準確把握鄉(xiāng)村發(fā)展問題，摸清鄉(xiāng)村發(fā)展規(guī)律，劃分鄉(xiāng)村類型，提出鄉(xiāng)村發(fā)展策略，已成為化解鄉(xiāng)村人地矛盾、促進鄉(xiāng)村發(fā)展的新思路。在村鎮(zhèn)國土空間規(guī)劃中，因地制宜劃分鄉(xiāng)村類型，解決鄉(xiāng)村問題，推動鄉(xiāng)村可持續(xù)發(fā)展具有重要指導作用。

1)土地利用多功能和鄉(xiāng)村發(fā)展相互促進的鄉(xiāng)村，具有較好的自然資源稟賦和社會經(jīng)濟基礎，需要保持土地利用多功能與鄉(xiāng)村發(fā)展水平同步增長，形成正向循環(huán)。結合區(qū)域機械電子、紡織產(chǎn)業(yè)、食品藥品等已有產(chǎn)業(yè)基礎，優(yōu)化產(chǎn)業(yè)用地布局，提高區(qū)域吸納就業(yè)的能力，擴大本地就業(yè)比例，在區(qū)位優(yōu)勢和政策支持的雙重驅動下，增強鄉(xiāng)村發(fā)展內(nèi)生動力，提升鄉(xiāng)村發(fā)展韌性。

2)土地利用多功能和鄉(xiāng)村發(fā)展促進和阻礙作用并存的鄉(xiāng)村，說明土地利用多功能與鄉(xiāng)村發(fā)展處于失衡狀態(tài)。對于土地利用多功能滯后鄉(xiāng)村，應以土地綜合整治為抓手，合理規(guī)劃村莊布局模式，加強農(nóng)田水利工程、水土保持工程、路網(wǎng)工程等農(nóng)業(yè)基礎設施建設，改善農(nóng)業(yè)生產(chǎn)條件，提高農(nóng)業(yè)機械化水平和生產(chǎn)技術，提升土地利用效率，增強區(qū)域土地利用多功能水平；對于鄉(xiāng)村發(fā)展滯后鄉(xiāng)村，應充分利用本底資源，加強人居環(huán)境改造，改善鄉(xiāng)村生產(chǎn)、生活條件，增加外出就業(yè)渠道，引導人員外出務工，提高農(nóng)戶生計能力，從而提高土地利用與鄉(xiāng)村發(fā)展協(xié)調(diào)度。

3)土地利用多功能和鄉(xiāng)村發(fā)展相互阻礙的鄉(xiāng)村，通過構建政府、企業(yè)、農(nóng)民等多主體聯(lián)動機制，制定鄉(xiāng)村幫扶政策，提高政府和企業(yè)等對基礎設建設的投入，發(fā)展水產(chǎn)養(yǎng)殖、油料作物和生豬養(yǎng)殖等特色農(nóng)業(yè)項目，從而提高鄉(xiāng)村自身發(fā)展能力，打破土地利用多功能和鄉(xiāng)村發(fā)展相互阻礙的格局，實現(xiàn)區(qū)域高質(zhì)量發(fā)展。

由于數(shù)據(jù)收集限制，本研究在土地利用多功能和鄉(xiāng)村發(fā)展指標選取的全面性和系統(tǒng)性有待進一步完善。土地利用多功能與鄉(xiāng)村發(fā)展交互作用關系是一個動態(tài)的過程，受到自然過程和人類活動的影響，時空驅動機制是耦合關系變化的重要組成部分，下一步需要加強對耦合協(xié)調(diào)關系變化的驅動機制研究。