林彤，高建崗，王亞華

（1.南京師范大學(xué)地理科學(xué)學(xué)院，南京 210023；2.重慶市勘測(cè)院，重慶 401121）

我國(guó)歷來(lái)重視耕地的利用與保護(hù)，“十分珍惜、合理利用土地，切實(shí)保護(hù)耕地”是我國(guó)的基本國(guó)策[1]。耕地質(zhì)量是自然、經(jīng)濟(jì)、區(qū)位等多種因素綜合作用的結(jié)果，是保障國(guó)家農(nóng)業(yè)生產(chǎn)安全和社會(huì)經(jīng)濟(jì)和諧的重要基石[2]。2021 年中央“一號(hào)文件”提出要落實(shí)最嚴(yán)格的耕地保護(hù)制度，進(jìn)一步明確耕地質(zhì)量保護(hù)的重要性。但隨著城市化的快速推進(jìn)，環(huán)境惡化、用地失衡、耕地質(zhì)量下降等問(wèn)題也日益凸顯，耕地空間分異顯著等問(wèn)題亟待解決。因此，開(kāi)展耕地質(zhì)量空間分布特征分析，研究耕地質(zhì)量空間分異規(guī)律，可為耕地質(zhì)量差異化保護(hù)、國(guó)土空間規(guī)劃編制提供參考。

國(guó)內(nèi)外學(xué)者主要對(duì)耕地質(zhì)量指標(biāo)體系[3?4]、評(píng)價(jià)[5?6]、尺度[7?8]、時(shí)空分異[9?10]、潛力提升[11?13]、時(shí)空演變[14?15]等方面進(jìn)行了深入探索。關(guān)瑜等[16]通過(guò)搭建耕地三維評(píng)價(jià)體系框架，利用熱點(diǎn)分析等方法分析耕地質(zhì)量三維空間關(guān)聯(lián)效應(yīng)，提出耕地利用分區(qū)措施。李子杰等[17]基于太湖縣耕地等級(jí)評(píng)價(jià)結(jié)果，通過(guò)時(shí)間序列方法和空間自相關(guān)解譯了耕地質(zhì)量空間格局。楊君等[18]通過(guò)布設(shè)岳陽(yáng)樓區(qū)采樣點(diǎn)，基于地塊微觀尺度，采用LISA 和聚類(lèi)分析方法探究岳陽(yáng)樓區(qū)耕地質(zhì)量動(dòng)態(tài)演變分析及積聚性特征，并制定耕地空間保護(hù)分區(qū)措施。李志芳等[19]采用空間轉(zhuǎn)換和LISA 分析探究耕地質(zhì)量集聚規(guī)律，并制定相應(yīng)的分區(qū)保護(hù)措施。劉慧芳等[20]從晉中市村級(jí)微觀尺度出發(fā)，采用空間自相關(guān)、地理探測(cè)器等分析耕地質(zhì)量要素聚類(lèi)特征及其影響機(jī)制。史云揚(yáng)等[21]通過(guò)構(gòu)建“三資”多維耕地質(zhì)量評(píng)價(jià)體系，研討了耕地質(zhì)量屬性?xún)?nèi)涵，運(yùn)用局部空間自相關(guān)方法分析了尋甸縣“三資”維度下空間分異格局。楊慧等[22]通過(guò)多維尺度研究陜西省耕地質(zhì)量等別結(jié)構(gòu)特征，并分析了耕地資源本底要素之間的關(guān)系。綜上，國(guó)內(nèi)外關(guān)于耕地質(zhì)量的評(píng)價(jià)方法與指標(biāo)體系建立的相關(guān)研究成果較為豐富，并且學(xué)者們從空間分異角度探究了研究對(duì)象的異質(zhì)性。傳統(tǒng)的研究分析方法很難定量反映耕地質(zhì)量的空間差異性及演變規(guī)律，文獻(xiàn)中從空間聚類(lèi)角度分析耕地質(zhì)量異質(zhì)性的研究較少，且先前的研究往往忽略了耕地質(zhì)量?jī)?nèi)部之間的差異，難以深入探討耕地質(zhì)量空間分異機(jī)制。

鑒于此，本研究基于市縣級(jí)尺度視角，以蘇北地區(qū)所屬各市縣為研究區(qū)域，將2020 年度耕地質(zhì)量三類(lèi)等別指數(shù)作為空間變量，探討蘇北地區(qū)市縣尺度下耕地質(zhì)量空間布局及趨勢(shì)分布，并采用空間聚類(lèi)法判別蘇北地區(qū)分區(qū)保護(hù)優(yōu)先級(jí)，劃定基于空間集聚與異質(zhì)性的分類(lèi)保護(hù)區(qū)，以期為國(guó)土空間規(guī)劃編制、耕地差異化管理和質(zhì)量提升提供理論依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

蘇北地區(qū)位于江蘇省北部（北緯31°09′～31°41′、東經(jīng)119°08′～119°36′），包含徐州、淮安、鹽城、宿遷和連云港5 個(gè)地級(jí)市（圖1），下轄37 個(gè)縣（市、區(qū)），總面積5.48 萬(wàn)km2，是我國(guó)沿海經(jīng)濟(jì)帶的重要組成部分。蘇北地區(qū)地勢(shì)以平原為主，屬溫帶季風(fēng)氣候，光溫水熱資源充足，河流湖泊眾多。蘇北地區(qū)作為江蘇省重要的糧食主產(chǎn)區(qū)域，是中國(guó)重要的“糧倉(cāng)”和“魚(yú)米之鄉(xiāng)”。

圖1 蘇北地區(qū)行政區(qū)劃示意圖Figure 1 Administrative divisions of northern Jiangsu

1.2 數(shù)據(jù)來(lái)源

數(shù)據(jù)來(lái)源主要包含蘇北地區(qū)2020 年度耕地質(zhì)量等別更新數(shù)據(jù)庫(kù)、土地利用現(xiàn)狀調(diào)查數(shù)據(jù)庫(kù)、國(guó)土空間規(guī)劃數(shù)據(jù)庫(kù)、“三調(diào)”成果數(shù)據(jù)庫(kù)，均來(lái)自江蘇省自然資源廳；遙感影像與高程數(shù)據(jù)（0.3 m 分辨率）來(lái)源于地理空間數(shù)據(jù)云網(wǎng)站。

本研究首先匹配校正遙感影像與耕地質(zhì)量等別數(shù)據(jù)庫(kù)坐標(biāo)系，兩者均采用2000 國(guó)家大地坐標(biāo)系；其次，通過(guò)數(shù)據(jù)清洗、掛接、驗(yàn)證檢查等，提取耕地質(zhì)量分析單元數(shù)據(jù)，采用面積加權(quán)法轉(zhuǎn)換耕地質(zhì)量地類(lèi)賦值，通過(guò)趨勢(shì)面分析研究區(qū)域內(nèi)耕地質(zhì)量空間分布。

1.3 研究方法

1.3.1 面積加權(quán)平均法

以蘇北地區(qū)市縣尺度為研究單元，綜合考慮耕地資源的自然稟賦、治理水平、經(jīng)濟(jì)條件、要素功能等特性，采用面積加權(quán)平均法對(duì)市縣區(qū)域耕地質(zhì)量分類(lèi)別指數(shù)進(jìn)行加權(quán)空間平衡轉(zhuǎn)換。具體公式如下：

式中：Ri為第i個(gè)縣域耕地質(zhì)量指數(shù)值；Rij和Sij分別為第i個(gè)縣域中第j塊圖斑的耕地質(zhì)量指數(shù)和圖斑面積；i和j分別為行政村和圖斑序號(hào)；n為圖斑總個(gè)數(shù)。

1.3.2 空間自相關(guān)

空間自相關(guān)描述的是研究區(qū)域內(nèi)測(cè)算單元屬性與周?chē)鷾y(cè)算單元內(nèi)空間要素的依賴(lài)程度，采用Moran′sI指數(shù)來(lái)判斷區(qū)域內(nèi)要素整體性的聚散程度[23]，可表達(dá)為：

式中：n代表研究對(duì)象的個(gè)數(shù)；表示空間單元屬性值的平均值；xi和xj分別表示空間單元i和j的屬性值；Wij為空間鄰接矩陣，表示空間單元i和j之間的鄰接關(guān)系，當(dāng)i和j相鄰時(shí)，Wij為1，不相鄰時(shí)則為0。

局部空間自相關(guān)用于分析空間某要素與相鄰要素的相關(guān)程度，表征其在內(nèi)部空間及鄰域范圍內(nèi)的聚集狀態(tài)，為后期的耕地利用分區(qū)工作提供依據(jù)。局部空間自相關(guān)分析用局域莫蘭指數(shù)Local Moran′sI表示，其計(jì)算公式如下：

式中：Ii為正值表示該空間單元與鄰近單元的屬性值近似，Ii為負(fù)值表示該空間單元與鄰近單元的屬性值差異較大。

2 結(jié)果與分析

2.1 耕地質(zhì)量空間分布特征

基于聚類(lèi)分析與趨勢(shì)性分析研究，得到蘇北區(qū)域耕地質(zhì)量市縣分布特征，如圖2 所示。蘇北地區(qū)市縣尺度下耕地自然等別較為集中，主要以中等地為主，低等地有少量分布。根據(jù)蘇北地區(qū)農(nóng)用地分等參數(shù)，蘇北地區(qū)標(biāo)準(zhǔn)耕作制度分區(qū)屬于江淮平原區(qū)，在熟制上以一年兩熟為主，將水稻作為蘇北地區(qū)的基準(zhǔn)作物，指定作物為水稻、小麥、玉米。其中，蘇北地區(qū)指定作物水稻、小麥、玉米光溫生產(chǎn)潛力指數(shù)分別為1 594、869、2 175，指定作物水稻、小麥、玉米產(chǎn)量比系數(shù)分別為1、1.3、0.94。蘇北全域耕地自然等別在5～7等之間，以6 等地為主。高等別耕地集中分布于蘇北地區(qū)的西北部和東北部，即海州區(qū)、沛縣、銅山縣、灌云縣等地，耕地條件良好，產(chǎn)量較高。低等別耕地分布于南部區(qū)域，主要位于淮安市和鹽城市的縣級(jí)市?？傮w上，蘇北地區(qū)耕地自然等呈現(xiàn)由北向南逐漸遞減的變化趨勢(shì)，在東西方向呈現(xiàn)倒“U”形的變化趨勢(shì)，先逐漸上升后逐漸下降。

圖2 蘇北地區(qū)耕地質(zhì)量等別空間分布Figure 2 Spatial distribution of cultivated land quality grades in northern Jiangsu

蘇北地區(qū)耕地利用等指數(shù)在市縣行政尺度上主要為5～8 等，以6 等地為主。土地利用系數(shù)用來(lái)衡量補(bǔ)充地塊所處區(qū)域自然生產(chǎn)潛力平均發(fā)揮水平及潛力的利用率，反映的是耕地在自然狀態(tài)下的平均生產(chǎn)力。運(yùn)用ArcGIS 軟件，參考修正因素的獲取和量化方法對(duì)等值區(qū)間進(jìn)行修正，按照分等單元計(jì)算各單元的因素分值，重新確定權(quán)重，通過(guò)土地利用系數(shù)修正方法對(duì)原有土地利用系數(shù)進(jìn)行修正，得到各評(píng)價(jià)單元修正因素的分值。低利用等指數(shù)零星分散在泗洪縣、淮陰區(qū)、射陽(yáng)縣等地。高等別指數(shù)分布于銅山縣、灌云縣、東?？h和淮安區(qū)等地，這些區(qū)域基礎(chǔ)設(shè)施、農(nóng)地整治技術(shù)、利用水平較強(qiáng)?？傮w上，蘇北地區(qū)耕地利用等呈現(xiàn)由北向南逐級(jí)遞減的趨勢(shì)。

蘇北地區(qū)耕地經(jīng)濟(jì)等指數(shù)在市縣行政尺度上主要為6～9 等，以7 等和8 等地為主。高等別主要分布在東北部的東?？h和灌云縣、南部的大豐市和淮安區(qū)；低等別耕地分布在泗洪縣、淮陰區(qū)、東臺(tái)市和射陽(yáng)縣等地。總體上，蘇北地區(qū)耕地經(jīng)濟(jì)等在東西方向上呈現(xiàn)出倒“U”形的趨勢(shì)，先逐漸上升后逐漸下降，南北方向上無(wú)明顯變化趨勢(shì)。

2.2 耕地質(zhì)量全局空間自相關(guān)分析

通過(guò)Rook 鄰接原則構(gòu)造空間權(quán)重，采用Open GeoDa 分析蘇北地區(qū)耕地質(zhì)量要素單元全局空間自相關(guān)，分析市縣行政尺度下3 類(lèi)指數(shù)的空間關(guān)聯(lián)程度，結(jié)果見(jiàn)表1。

表1 蘇北地區(qū)耕地質(zhì)量指數(shù)Moran′s I值Table 1 Moran′s I value of cultivated land quality index in northern Jiangsu

由表1 可知，蘇北地區(qū)耕地國(guó)家自然等、國(guó)家利用等、國(guó)家經(jīng)濟(jì)等的Moran′sI值分別為0.357、0.389、0.349，均大于0.2，且顯著性檢驗(yàn)P值均為0.001，明顯小于顯著性水平0.05。從全局Moran′sI值來(lái)看，國(guó)家利用等>國(guó)家自然等>國(guó)家經(jīng)濟(jì)等，表明蘇北地區(qū)耕地質(zhì)量在市縣尺度呈現(xiàn)較強(qiáng)的聚集趨勢(shì)。

2.3 耕地質(zhì)量局部空間自相關(guān)分析

為了進(jìn)一步分析耕地質(zhì)量集聚程度，采用Open GeoDa軟件對(duì)蘇北地區(qū)3類(lèi)耕地等別進(jìn)行局部空間自相關(guān)分析，表2 為市縣行政尺度下3 類(lèi)耕地質(zhì)量指數(shù)局部空間自相關(guān)類(lèi)型與市縣數(shù)量統(tǒng)計(jì)，圖3 為蘇北地區(qū)耕地質(zhì)量指數(shù)局部空間集聚圖。

表2 蘇北地區(qū)耕地質(zhì)量指數(shù)局部空間自相關(guān)類(lèi)型與市縣數(shù)量統(tǒng)計(jì)Table 2 Local spatial autocorrelation types and quantitative statistics of cultivated land quality index in northern Jiangsu

圖3 蘇北地區(qū)耕地質(zhì)量指數(shù)局部空間自相關(guān)（LISA）集聚圖Figure 3 Local spatial autocorrelation（LISA）clustering map of cultivated land quality index in northern Jiangsu

（1）耕地質(zhì)量國(guó)家自然等

根據(jù)表2與圖3a可知，以市縣為研究單元，HH型和LL 型分別包含5、2 個(gè)市縣，占蘇北地區(qū)市縣的18.9%，LH 型和HL 型分別包含3、1 個(gè)市縣，占蘇北地區(qū)市縣的10.8%，NN型包含26個(gè)市縣區(qū)域，占蘇北地區(qū)總數(shù)的70.3%。自然等指數(shù)LISA 分析中正相關(guān)類(lèi)型高于負(fù)相關(guān)類(lèi)型，且正相關(guān)類(lèi)型聚合特性較強(qiáng)。負(fù)相關(guān)類(lèi)型呈現(xiàn)零散分布。

其中，HH 型呈組團(tuán)狀集聚分布，主要分布在贛榆區(qū)、東?？h、海州區(qū)、灌云縣、連云區(qū)等地，HH 型區(qū)域地勢(shì)較為平緩，土壤肥力較高，自然條件較好，道路通達(dá)度高。LL 型主要以片狀形式分布，包含洪澤區(qū)、清江浦區(qū)等區(qū)域。LH和HL型呈零星分布，主要包含豐縣、沭陽(yáng)縣、響水縣、淮陰區(qū)等地區(qū)。

（2）耕地質(zhì)量國(guó)家利用等

根據(jù)表2 與圖3b 可知，以市縣為研究單元，HH型和LL 型分別包含6、3 個(gè)市縣，占蘇北地區(qū)市縣的24.3%，LH 型和HL 型分別各包含2 個(gè)市縣，占蘇北地區(qū)市縣的10.8%，NN型包含24個(gè)市縣區(qū)域，占蘇北地區(qū)總數(shù)的64.9%。利用等指數(shù)LISA 分析中正相關(guān)類(lèi)型高于負(fù)相關(guān)類(lèi)型，且正相關(guān)類(lèi)型聚合特性較強(qiáng)。負(fù)相關(guān)類(lèi)型呈現(xiàn)零散分布。

其中，HH 型呈組團(tuán)狀集聚分布，包括豐縣、東?？h、海州區(qū)、連云區(qū)、灌云縣、贛榆區(qū)等地區(qū)，HH 型區(qū)域耕地利用條件良好，機(jī)械化程度較高，生產(chǎn)便利，便于耕作和管理。LL 型則聚集在泉山區(qū)、鼓樓區(qū)、云龍區(qū)等地區(qū)，LL 型區(qū)域地形起伏大，耕地破碎化程度較高，利用條件較差，不利于農(nóng)業(yè)生產(chǎn)和發(fā)展。LH 型零散分布于沭陽(yáng)縣和響水縣。HL型分布在銅山區(qū)和淮陰區(qū)，管理水平有待進(jìn)一步提高。

（3）耕地質(zhì)量國(guó)家經(jīng)濟(jì)等

根據(jù)表2與圖3c可知，以市縣為研究單元，HH型和LL 型分別包含4、3 個(gè)市縣，占蘇北地區(qū)市縣的18.9%，LH 型和HL 型分別包含2、1 個(gè)市縣，占蘇北地區(qū)市縣的8.1%，NN 型包含27 個(gè)市縣區(qū)域，占蘇北地區(qū)總數(shù)的73.0%。

其中，HH 型通過(guò)組團(tuán)狀集聚分布，主要分布在東海縣、海州區(qū)、連云區(qū)、灌云縣等地，HH 型區(qū)域積極調(diào)整農(nóng)業(yè)生產(chǎn)方式，道路通達(dá)度高，耕地投入產(chǎn)出較高，發(fā)展農(nóng)業(yè)規(guī)模化建設(shè)，經(jīng)濟(jì)效益也隨之提高。LL 型主要以片狀形式分布，主要分布在鼓樓區(qū)、泉山區(qū)和云龍區(qū)等區(qū)域。LH和HL型呈零星分布，主要包含銅山區(qū)、贛榆區(qū)和沭陽(yáng)縣。

今年71歲的蘆廷元，是一名退離休黨員干部，四十三年前的現(xiàn)在，他正是一名風(fēng)華正茂剛走出學(xué)校門(mén)的青年。一個(gè)偶然的機(jī)會(huì)，他跟隨一個(gè)回老家探親的親戚來(lái)新疆游玩，那時(shí)正是國(guó)家三年自然災(zāi)害的最后年頭，看到兵團(tuán)的農(nóng)場(chǎng)一派生機(jī)盎然的景象，平整的條田，筆直的林帶，四通八達(dá)的道路，機(jī)械化生產(chǎn)的模式，給他留下了深刻印象，他當(dāng)時(shí)就產(chǎn)生了留疆的念頭，而后被分配到農(nóng)四團(tuán)八連工作。

2.4 耕地保護(hù)分區(qū)

根據(jù)空間極化理論，聚集分布和離散分布在空間極化過(guò)程中呈現(xiàn)出明顯特征。根據(jù)縣級(jí)耕地質(zhì)量等別指數(shù)局部空間自相關(guān)LISA 分析，構(gòu)建耕地質(zhì)量保護(hù)分區(qū)規(guī)則體系（表3）。根據(jù)LISA 多維度集聚類(lèi)型組合特征，對(duì)蘇北地區(qū)進(jìn)行耕地保護(hù)分區(qū)，劃分為核心保護(hù)區(qū)、質(zhì)量提升區(qū)、綜合改良區(qū)、結(jié)構(gòu)調(diào)整區(qū)4 類(lèi)保護(hù)分區(qū)（圖4），并提出對(duì)應(yīng)的分區(qū)保護(hù)措施。

表3 基于局部空間自相關(guān)類(lèi)型組合的蘇北耕地綜合質(zhì)量分區(qū)Table 3 Partition of comprehensive quality of cultivated land in northern Jiangsu based on the combination of local spatial autocorrelation types

圖4 蘇北地區(qū)耕地保護(hù)分區(qū)Figure 4 Cultivated land protection zone in northern Jiangsu

HH 型為高質(zhì)量集聚度范圍，應(yīng)劃為核心保護(hù)區(qū)，強(qiáng)化其對(duì)周邊耕地的輻射效應(yīng)。LL 型為低質(zhì)量集聚度范圍，應(yīng)劃為綜合改良區(qū)，改善農(nóng)業(yè)生產(chǎn)條件，提高生產(chǎn)功能。HL 型和LH 型高等地與低等地交織環(huán)繞，應(yīng)劃為質(zhì)量提升區(qū)，通過(guò)重點(diǎn)改良舉措，避免低質(zhì)量耕地進(jìn)行非農(nóng)建設(shè)，使其轉(zhuǎn)化為高質(zhì)量耕地。NN型無(wú)明顯集聚性分布特征，應(yīng)劃為結(jié)構(gòu)調(diào)整區(qū)，因地制宜，加強(qiáng)土地整理，采取綜合性調(diào)整措施，優(yōu)化耕地空間布局，提高糧食產(chǎn)量。

（1）核心保護(hù)區(qū)

基于局部空間自相關(guān)類(lèi)型分類(lèi)規(guī)則，蘇北地區(qū)核心保護(hù)區(qū)以6、7 等地為主，主要分布在贛榆區(qū)、灌云縣及周邊區(qū)域。該區(qū)域高質(zhì)量耕地聚集分布，自然稟賦較高，地勢(shì)較為平坦、臨近水源，土壤肥沃，耕地破碎化程度較低。同時(shí)，該區(qū)域擁有完善的道路體系，基礎(chǔ)設(shè)施齊全，耕地利用條件較好，應(yīng)繼續(xù)加強(qiáng)耕地保護(hù)，優(yōu)化種植模式，建設(shè)耕地質(zhì)量動(dòng)態(tài)監(jiān)測(cè)系統(tǒng)，強(qiáng)化其對(duì)周邊耕地的輻射效應(yīng)。

（2）質(zhì)量提升區(qū)

基于局部空間自相關(guān)類(lèi)型分類(lèi)規(guī)則，蘇北地區(qū)質(zhì)量提升區(qū)以5、6 等地為主，主要分布在豐縣、銅山區(qū)、沭陽(yáng)縣、響水縣和淮陰區(qū)。該區(qū)域高等地與低等地交織環(huán)繞、或高等地覆蓋低等地，應(yīng)重點(diǎn)保護(hù)高質(zhì)量耕地；同時(shí)對(duì)低質(zhì)量地塊進(jìn)行土地整治，完善土地整治優(yōu)化措施，提升低質(zhì)量耕地產(chǎn)出率，在生態(tài)資源較好的區(qū)域發(fā)展現(xiàn)代生態(tài)農(nóng)業(yè)；通過(guò)重點(diǎn)改良舉措，避免低質(zhì)量耕地開(kāi)展非農(nóng)建設(shè)，使其轉(zhuǎn)化為高質(zhì)量耕地。

（3）綜合改良區(qū)

基于局部空間自相關(guān)類(lèi)型分類(lèi)規(guī)則，蘇北地區(qū)綜合改良區(qū)以4、5 等地為主，主要分布在邳州市?盱眙縣沿線(xiàn)、灌南縣?東臺(tái)市沿線(xiàn)。綜合改良區(qū)需結(jié)合農(nóng)業(yè)發(fā)展布局，對(duì)低等耕地地塊進(jìn)行土地整治，完善農(nóng)業(yè)基礎(chǔ)配套設(shè)施，改善農(nóng)業(yè)生產(chǎn)條件，提高生產(chǎn)功能；突出區(qū)域傳統(tǒng)產(chǎn)業(yè)與特色產(chǎn)業(yè)，增加農(nóng)產(chǎn)品生產(chǎn)與供應(yīng)，滿(mǎn)足多元化發(fā)展需求，實(shí)現(xiàn)增產(chǎn)增收，提高耕地質(zhì)量。盱眙縣等丘陵山區(qū)可發(fā)展薄殼山核桃等特色種植業(yè)；漣水縣可適當(dāng)發(fā)展蘆筍種植業(yè)。

（4）結(jié)構(gòu)調(diào)整區(qū)

結(jié)構(gòu)調(diào)整區(qū)主要分布在鼓樓區(qū)、泉山區(qū)、云龍區(qū)、洪澤區(qū)、清江浦區(qū)。蘇北地區(qū)結(jié)構(gòu)調(diào)整區(qū)耕地質(zhì)量LISA 關(guān)系無(wú)明顯集聚性分布特征，空間格局表現(xiàn)為隨機(jī)分布。結(jié)構(gòu)調(diào)整區(qū)耕地自然資源投入產(chǎn)出相對(duì)較差，應(yīng)根據(jù)蘇北地區(qū)結(jié)構(gòu)調(diào)整區(qū)的現(xiàn)狀發(fā)展情況，因地制宜，加強(qiáng)土地整理，采取結(jié)構(gòu)調(diào)整措施，優(yōu)化耕地空間布局，提高糧食產(chǎn)量。同時(shí)，對(duì)于泉山區(qū)等生態(tài)農(nóng)業(yè)示范區(qū)的耕地，可加強(qiáng)生態(tài)農(nóng)業(yè)品牌建設(shè)，發(fā)展優(yōu)勢(shì)生態(tài)農(nóng)業(yè)產(chǎn)業(yè)，同時(shí)逐步拓展生態(tài)旅游模式。

3 討論

本研究以蘇北地區(qū)市縣尺度的耕地質(zhì)量等別指數(shù)為空間變量，借助ArcGIS、Open Geoda 等軟件，運(yùn)用趨勢(shì)性分析方法解析了蘇北地區(qū)耕地質(zhì)量空間布局及趨勢(shì)分布。通過(guò)空間LISA 自相關(guān)分析，揭示了蘇北地區(qū)耕地質(zhì)量的空間分異特征，劃定蘇北地區(qū)基于空間集聚與異質(zhì)性的分類(lèi)保護(hù)區(qū)，制定符合實(shí)際的分區(qū)保護(hù)模式。

本研究仍存在不足之處，將從以下方面進(jìn)行改進(jìn)：

（1）深入挖掘耕地質(zhì)量各指標(biāo)之間的空間分異規(guī)律，以進(jìn)一步探析空間分異影響機(jī)制。

（2）深入探究多尺度耕地質(zhì)量空間自相關(guān)性，基于市縣、鄉(xiāng)鎮(zhèn)、村級(jí)尺度探究耕地質(zhì)量空間分異特征，基于宏觀和微觀尺度研究耕地質(zhì)量空間集聚規(guī)律。

4 結(jié)論

（1）蘇北地區(qū)市縣尺度下耕地自然等別分布較為集中，以中等地為主。自然等指數(shù)呈現(xiàn)由北向南逐漸遞減的變化趨勢(shì)，在東西方向呈現(xiàn)倒“U”形的變化趨勢(shì)，先逐漸上升后逐漸下降。利用等指數(shù)呈現(xiàn)由北向南逐級(jí)遞減趨勢(shì)。經(jīng)濟(jì)等指數(shù)在東西方向上呈現(xiàn)倒“U”形的變化趨勢(shì)，先逐漸上升后逐漸下降，南北方向上無(wú)明顯變化趨勢(shì)。

（2）蘇北地區(qū)耕地國(guó)家自然等、國(guó)家利用等、國(guó)家經(jīng)濟(jì)等的Moran′sI值分別為0.357、0.389、0.349，表明蘇北地區(qū)耕地質(zhì)量在市縣尺度呈現(xiàn)較強(qiáng)的聚集趨勢(shì)。受投入產(chǎn)出等經(jīng)濟(jì)條件的影響，人類(lèi)活動(dòng)對(duì)蘇北地區(qū)耕地經(jīng)濟(jì)質(zhì)量影響空間差異明顯。蘇北地勢(shì)以平原為主，土地利用程度較高，其耕地利用質(zhì)量和自然質(zhì)量空間分異明顯，且空間關(guān)聯(lián)度較強(qiáng)。

（3）蘇北地區(qū)耕地質(zhì)量指數(shù)LISA 分析中正相關(guān)類(lèi)型高于負(fù)相關(guān)類(lèi)型，且正相關(guān)類(lèi)型聚集性較強(qiáng)，負(fù)相關(guān)類(lèi)型呈現(xiàn)零散分布。高?高型通過(guò)組團(tuán)狀集聚分布，低?低型主要以片狀、條帶狀形式分布，低?高和高?低型呈零星分布。根據(jù)LISA 多維度集聚類(lèi)型組合特征，蘇北地區(qū)耕地可劃分為核心保護(hù)區(qū)、質(zhì)量提升區(qū)、綜合改良區(qū)、結(jié)構(gòu)調(diào)整區(qū)4 類(lèi)保護(hù)分區(qū)，并實(shí)施對(duì)應(yīng)的分區(qū)保護(hù)措施。