基于國家耕地質(zhì)量等別指數(shù)和空間自相關(guān)的蘇北地區(qū)耕地保護(hù)分區(qū)

林彤，高建崗，王亞華

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

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

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

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

1 材料與方法

1.1 研究區(qū)概況

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

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

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

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

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

1.3 研究方法

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

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

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

1.3.2 空間自相關(guān)

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

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

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

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

2 結(jié)果與分析

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

基于聚類分析與趨勢性分析研究，得到蘇北區(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ū)域，主要位于淮安市和鹽城市的縣級市。總體上，蘇北地區(qū)耕地自然等呈現(xiàn)由北向南逐漸遞減的變化趨勢，在東西方向呈現(xiàn)倒“U”形的變化趨勢，先逐漸上升后逐漸下降。

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

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

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

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

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

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

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

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

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

表2 蘇北地區(qū)耕地質(zhì)量指數(shù)局部空間自相關(guān)類型與市縣數(shù)量統(tǒng)計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ì)量國家自然等

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

基于局部空間自相關(guān)類型分類規(guī)則，蘇北地區(qū)綜合改良區(qū)以4、5 等地為主，主要分布在邳州市?盱眙縣沿線、灌南縣?東臺市沿線。綜合改良區(qū)需結(jié)合農(nóng)業(yè)發(fā)展布局，對低等耕地地塊進(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)，滿足多元化發(fā)展需求，實現(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)系無明顯集聚性分布特征，空間格局表現(xiàn)為隨機(jī)分布。結(jié)構(gòu)調(diào)整區(qū)耕地自然資源投入產(chǎn)出相對較差，應(yīng)根據(jù)蘇北地區(qū)結(jié)構(gòu)調(diào)整區(qū)的現(xiàn)狀發(fā)展情況，因地制宜，加強(qiáng)土地整理，采取結(jié)構(gòu)調(diào)整措施，優(yōu)化耕地空間布局，提高糧食產(chǎn)量。同時，對于泉山區(qū)等生態(tài)農(nóng)業(yè)示范區(qū)的耕地，可加強(qiáng)生態(tài)農(nóng)業(yè)品牌建設(shè)，發(fā)展優(yōu)勢生態(tài)農(nóng)業(yè)產(chǎn)業(yè)，同時逐步拓展生態(tài)旅游模式。

3 討論

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

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

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

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

4 結(jié)論

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

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

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