基于遙感的海東市樂都區(qū)撂荒耕地空間分布特征研究

摘 要 為識別研究區(qū)撂荒地并了解其影響因素，幫助制定和完善相關(guān)政策，減少耕地負(fù)面影響，進(jìn)一步保障糧食安全，選取青海省東部重要農(nóng)業(yè)產(chǎn)業(yè)基地、具區(qū)域代表性（平原、山地均具備）、能體現(xiàn)區(qū)域特色的地域研究單元海東市樂都區(qū)，以遙感數(shù)據(jù)、青海省地理空間數(shù)據(jù)為信息源，基于農(nóng)作物的生長周期以季度變化構(gòu)建NDVI時間序列數(shù)據(jù)，通過比較不同時間段內(nèi)的NDVI數(shù)值差異來實(shí)現(xiàn)研究區(qū)撂荒地識別研究，并通過核密度分析、空間自相關(guān)分析和冷熱點(diǎn)分析等方法探討研究區(qū)撂荒地的空間分布特征。結(jié)果顯示，研究區(qū)耕地撂荒率近些年來相對穩(wěn)定，撂荒地主要集中在研究區(qū)山間道路兩側(cè)，這些區(qū)域坡耕地坡度較大，從而導(dǎo)致了大量碎化耕地被撂荒。政府可通過制定一系列農(nóng)業(yè)政策，如糧食補(bǔ)貼，促進(jìn)糧食種植，從而保障糧食安全。

關(guān)鍵詞 撂荒地；NDVI時間序列；核密度分析；空間自相關(guān)；糧食安全；青海省海東市樂都區(qū)

中圖分類號：F323.211 文獻(xiàn)標(biāo)志碼：A DOI：10.19415/j.cnki.1673-890x.2024.01.015

中國是農(nóng)業(yè)生產(chǎn)大國[1]。作為人類生存的物質(zhì)基礎(chǔ)，農(nóng)田在糧食生產(chǎn)和生物多樣性保護(hù)方面發(fā)揮著至關(guān)重要的作用。改革開放以來，中央實(shí)施了一系列措施補(bǔ)充建成用地占用的耕地，以保持總面積的動態(tài)平衡，如耕地戰(zhàn)略儲備制度、耕地征用補(bǔ)償平衡等[2]。但在過去十年中，耕地仍然發(fā)生了重大變化，有條件優(yōu)良的耕地被轉(zhuǎn)化為了建設(shè)用地，有退耕還林還草，有坡度高于2°的森林被清理用于耕種，也有各類其他用地轉(zhuǎn)化為耕地。在耕地再分配過程中，新增耕地逐漸向海拔較高、坡度較陡的地區(qū)轉(zhuǎn)移，坡耕地越來越普遍[3]。因此無論是開發(fā)還是開墾，高海拔、坡度較高的耕地補(bǔ)充，已成為耕地資源的重要組成部分，對保障糧食安全具有至關(guān)重要的作用。然而，現(xiàn)有的研究主要集中在平原地帶優(yōu)質(zhì)耕地的研究，對山區(qū)農(nóng)田是否得到長期有效利用關(guān)注不夠[4]。

青海省海東市樂都區(qū)位于西北地區(qū)農(nóng)牧交錯地帶，是青海省東部重要的農(nóng)業(yè)產(chǎn)業(yè)基地，但該區(qū)域大部分是坡耕地，與原始耕地相比，在機(jī)械化、大規(guī)模運(yùn)營和農(nóng)業(yè)支持措施方面處于不利地位，限制了糧食產(chǎn)量，且隨著農(nóng)業(yè)生產(chǎn)材料和勞動力成本的增加，坡耕地的糧食生產(chǎn)越來越具有高投入和低產(chǎn)出的特征[5]。當(dāng)產(chǎn)量不符合農(nóng)民的期望時，條件較差的耕地往往會被撂荒。耕地撂荒破壞了耕地利用的連續(xù)性，是對土地資源的嚴(yán)重浪費(fèi)，對糧食生產(chǎn)和生態(tài)平衡產(chǎn)生了巨大影響。因此本文通過識別研究區(qū)耕地的撂荒情況，分析青海省東部農(nóng)牧交錯地帶撂荒耕地的空間分布格局及其變化，量化耕地的連續(xù)性，以期為該區(qū)域合理管理耕地利用、保障未來糧食安全提供依據(jù)。

1" 研究區(qū)概況

樂都區(qū)是海東市的縣轄區(qū)，位于青海省東部湟水河中下游，全區(qū)轄7鎮(zhèn)12鄉(xiāng)，總面積3 050 km2（圖 1，102°09＇～102°47＇E，36°16＇～36°46＇N）。農(nóng)作物分布在這一地區(qū)，主要類型為小麥、馬鈴薯、油菜等[6]。該區(qū)地貌以山地為主，地勢錯綜復(fù)雜，氣候不穩(wěn)定，云層干擾頻繁，耕地主要分布于湟水河兩岸和南北山地，海拔相差最大可達(dá)2 644 m。此外該區(qū)域經(jīng)濟(jì)不發(fā)達(dá)，人均GDP遠(yuǎn)低于全國平均水平，高比例的農(nóng)民工不再從事農(nóng)業(yè)工作，這可能導(dǎo)致疏于管理的耕地被撂荒。

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

2.1" 數(shù)據(jù)來源

本研究所用遙感數(shù)據(jù)為Sentinel-2遙感影像（分辨率為10 m），來源于GEE平臺（https：//earthengine.google.com/），Google Earth Engine（GEE）集成了大規(guī)模遙感影像數(shù)據(jù)和現(xiàn)有的分類產(chǎn)品，以及大量的函數(shù)和現(xiàn)有算法，可應(yīng)用于大規(guī)模土地監(jiān)測[7]。本文所用NDVI數(shù)據(jù)基于Sentinel-2影像，通過波段計(jì)算獲取得到月最大值NDVI合成數(shù)據(jù)，時間序列NDVI數(shù)據(jù)將用于撂荒地識別。研究區(qū)邊界矢量數(shù)據(jù)來源于國家基礎(chǔ)地理信息中心（www.webmap.cn），樂都區(qū)第三次全國國土調(diào)查數(shù)據(jù)來源于青海省自然資源廳（https：//zrzyt.qinghai.gov.cn/index）。

2.2" 研究方法

2.2.1" NDVI變化監(jiān)測

NDVI（歸一化植被指數(shù)）是一種廣泛用于遙感影像分析的指標(biāo)，可用于估算植被覆蓋度和生長狀況[8]。該指標(biāo)基于植被對不同波長的反射率的變化來計(jì)算，可以通過計(jì)算近紅外波段和紅色波段之間的比率來獲得。NDVI的值范圍從-1～1，值越高表示植被覆蓋度越高。NDVI的變化監(jiān)測通常用于研究植被的生長和變化趨勢，例如在不同季節(jié)或多年的時間段內(nèi)監(jiān)測植被覆蓋度的變化。

本文基于農(nóng)作物的生長周期以季度變化構(gòu)建NDVI時間序列數(shù)據(jù)，通過比較不同時間段內(nèi)的NDVI圖像來實(shí)現(xiàn)，以確定植被的生長或退化情況。具體內(nèi)容為基于樂都區(qū)2019年三調(diào)耕地圖斑，隨機(jī)選取三調(diào)屬性表中耕種樣本和未耕種樣本各50個，通過對比研究區(qū)耕種樣本和非耕種樣本春夏NDVI值變化的差異（表 1），確定兩種地類的最佳分割點(diǎn)。然后連續(xù)調(diào)整分割點(diǎn)直到所有的樣本歸類到對應(yīng)的樣本集中，最終確定分割閾值為0.33。

本文中，我們將耕地一年以內(nèi)未種植的情況叫作輪休，持續(xù)時間超過一年未種植的情況我們作做撂荒。通過運(yùn)用ArcGIS軟件，對研究區(qū)的NDVI數(shù)據(jù)進(jìn)行耕地區(qū)域的掩膜處理，得到耕地地塊的NDVI數(shù)值，對該NDVI數(shù)值進(jìn)行春夏和夏秋的差值運(yùn)算，并進(jìn)行重分類。把NDVI差值大于0.33的判定為本年度耕種地塊，小于0.33的判定為本年度未耕種地塊，將相鄰年份未耕種地塊進(jìn)行疊加，重疊的地塊判定為撂荒地，沒有重疊的地塊判定為當(dāng)年未耕種，具體撂荒地提取流程圖如圖 2。

2.2.2" 核密度分析

核密度分析是一種空間數(shù)據(jù)分析方法，用于確定數(shù)據(jù)分布在空間上的密度模式[9]。在核密度分析中，數(shù)據(jù)集被視為一組離散的事件，這些事件被認(rèn)為在空間上隨機(jī)分布。該分析可以確定在給定區(qū)域內(nèi)存在哪些地方有更高的密度，并可用于探究某一現(xiàn)象的分布模式和熱點(diǎn)區(qū)域的空間分布。本研究通過 ArcGIS 軟件Spatial Analyst Tools 中“Density”下的“Kernel Density”分析工具對研究區(qū)撂荒耕地點(diǎn)狀數(shù)據(jù)按搜索半徑為 3 km 進(jìn)行核密度分析。然后再在 ArcGIS 里利用自然裂點(diǎn)法將核密度分析結(jié)果重分類成 5 個級別，按照核密度從小到大依次為低密度區(qū)、中低密度區(qū)、中密度區(qū)、中高密度區(qū)和高密度區(qū)。

2.2.3" 空間自相關(guān)分析

空間自相關(guān)分析是一種用于研究地理空間數(shù)據(jù)中空間相關(guān)性的統(tǒng)計(jì)方法[10]，其基本思想是檢測地理空間數(shù)據(jù)集中的觀測值之間的相似性或相關(guān)性，通常使用空間權(quán)重矩陣來量化觀測值之間的空間鄰近性，然后計(jì)算各個觀測值與其鄰近觀測值之間的相關(guān)性指標(biāo)。本研究采用全局自相關(guān)分析和冷熱點(diǎn)分析來研究區(qū)域撂荒耕地空間集聚性，二者分別用 Moran's [I] 和 Getis-Ord [Gi?]指數(shù)進(jìn)行分析。

1）全局自相關(guān)分析。本研究用全局自相關(guān)分析了解研究區(qū)撂荒地集中的空間模式和集聚程度，判斷研究區(qū)撂荒地是否集聚。全局Moran's [I]的取值范圍為-1～1，其中正值表示正相關(guān)性且分布集聚，負(fù)值表示負(fù)相關(guān)性且呈分散型分布，而接近于0的值表示不存在空間自相關(guān)。定義如下：

式中，[zi]、[zj]分別為撂荒耕地圖斑[i]、[j]的面積與其平均值的偏差；[wij]為撂荒耕地圖斑[i]和[j]的空間權(quán)重；[S0]為所有空間權(quán)重集合；[n]為撂荒耕地圖斑總數(shù)；[Z]為標(biāo)準(zhǔn)化值；[E[I]]和[V[I]] 分別為期望值和方差。

2）冷熱點(diǎn)分析。冷熱點(diǎn)分析能夠識別和定量評估撂荒耕地空間集聚的范圍和集聚程度。進(jìn)行GIS冷熱點(diǎn)分析的常用方法之一是基于空間統(tǒng)計(jì)指標(biāo)的Getis-Ord [Gi?]統(tǒng)計(jì)量。[Gi?]為正時，表示出現(xiàn)熱點(diǎn)集聚，[Gi?]越大說明高值集聚程度越高；當(dāng) [Gi?]為負(fù)時，表示出現(xiàn)冷點(diǎn)集聚，[Gi?]越小說明低值聚類緊密。其計(jì)算公式為：

式中，[xi]、[xj]分別為撂荒耕地圖斑[i]和[j]的面積。

3" 結(jié)果與分析

3.1" 撂荒耕地分布及變化

研究區(qū)2019—2022年撂荒地空間分布如圖 3所示，從整體分布來看，撂荒地主要集中在東部地區(qū)和湟水河向南北輻射的山間道路，這是因?yàn)椴糠制赂仄露容^大且遠(yuǎn)離居民點(diǎn)，難以利用機(jī)械進(jìn)行耕種，導(dǎo)致撂荒耕地分布較為集中；部分年份處于湟水河兩岸的城鎮(zhèn)地區(qū)撂荒地也較為集中，這是因?yàn)椴糠帜攴莩鞘袛U(kuò)張城鎮(zhèn)建設(shè)用地占用耕地導(dǎo)致耕地處于常年撂荒的現(xiàn)象。研究區(qū)2019—2022年份的撂荒地面積和撂荒率如表 2所示，整體來看近些年份研究區(qū)耕地撂荒率較為穩(wěn)定且有持續(xù)變好的現(xiàn)象，這也與政府部門對耕地的持續(xù)重視和相關(guān)土地政策實(shí)施有關(guān)。

3.2" 撂荒耕地核密度分析

對研究區(qū)2020年撂荒地塊進(jìn)行核密度分析，結(jié)果如圖 4所示。由圖可知，研究區(qū)撂荒耕地分布具有集聚性，撂荒耕地核密度最大值3.04點(diǎn)/km2，最小值為0，平均值為0.27點(diǎn)/km2。研究區(qū)中密度區(qū)域集中在共和鄉(xiāng)、碾伯鎮(zhèn)、蒲臺鄉(xiāng)和洪水鎮(zhèn)，高密度區(qū)域集中在洪水鎮(zhèn)。通過實(shí)地和Google Earth Pro軟件中影像對這些區(qū)域的了解，碾伯鎮(zhèn)處于湟水河兩岸的城鎮(zhèn)地區(qū)，城鎮(zhèn)化的建設(shè)是導(dǎo)致該地區(qū)撂荒的最大因素。共和鄉(xiāng)、蒲臺鄉(xiāng)和洪水鎮(zhèn)的耕地大部分處于湟水河向南北山地輻射的山區(qū)道路兩側(cè)，這些地區(qū)的耕地的坡度更大，離居民點(diǎn)相對較遠(yuǎn)且地塊較為碎片化，農(nóng)耕機(jī)具很難得到利用，這也是該地區(qū)撂荒地較為集中的原因。

3.3" 撂荒耕地空間自相關(guān)分析

1）全局自相關(guān)分析。本文基于核密度值做研究區(qū)撂荒地空間自相關(guān)分析，運(yùn)用Geoda軟件對研究區(qū)撂荒耕地進(jìn)行全局自相關(guān)分析。結(jié)果顯示，Moran′s [I] 為 0.519 0，Z-value值為5.403 1。表明研究區(qū)撂荒地核密度在空間上呈集聚分布的態(tài)勢。

2）冷熱點(diǎn)分析。本文基于核密度值，利用 ArcGIS 的“Hot Spot Anal?ysis（Getis-Ord [Gi?]）”工具對研究區(qū)撂荒耕地進(jìn)行冷熱點(diǎn)分析，結(jié)果顯示，研究區(qū)撂荒耕地[Gi?]指數(shù)為-1.546 18～2.773 27，根據(jù) [Gi?]指數(shù)值的大小，利用自然裂點(diǎn)法將冷熱點(diǎn)分析結(jié)果重分類成 5 級（圖 5）。研究區(qū)低集聚區(qū)和較低集聚區(qū)屬于冷點(diǎn)區(qū)，高集聚區(qū)和較高集聚區(qū)屬于熱點(diǎn)區(qū)，其中低集聚區(qū)分布在下營藏族鄉(xiāng)、城臺鄉(xiāng)和峰堆鄉(xiāng)，較低集聚區(qū)分布在瞿曇鎮(zhèn)、李家鄉(xiāng)、馬營鄉(xiāng)和中嶺鄉(xiāng)，較高集聚區(qū)分布在達(dá)拉土族鄉(xiāng)、壽樂鎮(zhèn)、雨潤鎮(zhèn)、碾伯鎮(zhèn)和中壩藏族鄉(xiāng)，高集聚區(qū)分布在高廟鎮(zhèn)、洪水鎮(zhèn)和蒲臺鄉(xiāng)。從空間分布來看，熱點(diǎn)區(qū)域主要分布在湟水河兩岸的城鎮(zhèn)地區(qū)和研究區(qū)東南部，前者是由于城鎮(zhèn)化建設(shè)導(dǎo)致農(nóng)用地被占用，后者是由于該地區(qū)地形破碎復(fù)雜，溝壑縱橫，耕作條件較差且遠(yuǎn)離居民點(diǎn)，從而導(dǎo)致耕地被撂荒；冷點(diǎn)區(qū)域主要集中在南北的鄉(xiāng)鎮(zhèn)，這是由于研究區(qū)南北山地海拔較高，但坡耕地的坡度更為平緩，也使得該地區(qū)耕地被撂荒的情況較輕。

4" 結(jié)論

以全國第三次土地調(diào)查中的耕地土地利用圖對遙感圖像進(jìn)行掩膜，去除耕地范圍外的變化地物干擾，僅在耕地范圍內(nèi)提取撂荒地，并使用Sentinel系列（分辨率為10 m）數(shù)據(jù)提取的NDVI數(shù)據(jù)進(jìn)行撂荒地識別提取，相較于二調(diào)數(shù)據(jù)和遙感解譯提取耕地信息結(jié)果準(zhǔn)確性更高。結(jié)果顯示：研究區(qū)耕地整體撂荒率近些年相對穩(wěn)定，雖有起伏但沒有大的變化，從空間分布圖和撂荒地核密度分析及空間自相關(guān)分析結(jié)果來看，研究區(qū)耕地主要分布在南北山地和沿湟水河兩岸，而撂荒地主要集中在湟水河向南北輻射的山間道路兩側(cè)，這與研究區(qū)的地形地勢有很大的關(guān)系。湟水河兩岸土地相對平整，部分年份因?yàn)槌鞘薪ㄔO(shè)占用耕地導(dǎo)致撂荒地較為集中；在通向南北山地的山間地帶，由于坡耕地坡度較大導(dǎo)致很多耕地常年被撂荒，這其中有耕作條件的限制，也有當(dāng)?shù)鼐用癯Ｄ暌蕾嚨霓r(nóng)作物收入較低轉(zhuǎn)至其他方式就業(yè)而導(dǎo)致撂荒的因素。

耕地撂荒的原因主要受有人類活動和自然環(huán)境的影響。高海拔且耕地地面的不平整導(dǎo)致了山區(qū)耕地的碎片化，從而抑制了規(guī)?；N植，導(dǎo)致很多耕地常年疏于管理被撂荒。而山區(qū)經(jīng)濟(jì)的發(fā)展滯緩，使得城鄉(xiāng)遷移，農(nóng)業(yè)活動勞動力短缺，進(jìn)一步導(dǎo)致了耕地撂荒現(xiàn)象的加劇。鑒于撂荒耕地在研究區(qū)的分布，政府應(yīng)采取有針對性的措施，改善耕地利用。政府應(yīng)進(jìn)一步增加對農(nóng)村農(nóng)業(yè)的補(bǔ)貼，縮小收入差距；應(yīng)建立相應(yīng)的耕地出讓制度以減少人為因素造成的耕地碎片化，加強(qiáng)集聚分散耕地，提高農(nóng)業(yè)投入產(chǎn)出比，從而保障國家糧食安全。

參考文獻(xiàn)：

[1] 王學(xué)，李秀彬，宋恒飛.我國耕地撂荒問題分析及其對策研究[J].中國土地，2023（2）：15-17.

[2] 陳美球.耕地經(jīng)營者的耕地保護(hù)責(zé)任與權(quán)益探析[J].中國土地，2023（2）：12-14.

[3] 范戰(zhàn)平，趙啟航.耕地占補(bǔ)平衡制度：歷程·問題·建議[J].哈爾濱師范大學(xué)社會科學(xué)學(xué)報(bào)，2023，14（1）：66-69.

[4] 吳峰，李沛鴻，熊凡，等.基于遙感的南方丘陵山區(qū)耕地撂荒研究[J].北京測繪，2022，36（11）：1519-1523.

[5] 龍明順，趙宇鸞，張東麗.山區(qū)耕地細(xì)碎化對農(nóng)戶耕地撂荒的影響[J].農(nóng)業(yè)工程學(xué)報(bào)，2022，38（21）：231-239.

[6] 馬扶林， 周秉榮， 許正福. 基于GIS的青海省主要農(nóng)作物種植區(qū)劃[J]. 現(xiàn)代農(nóng)業(yè)科技， 2016 （12） ： 207-208，212.

[7] 王小娜， 田金炎， 李小娟， 等. 谷歌地球引擎云平臺對遙感發(fā)展的改變[J]. 遙感學(xué)報(bào)， 2022， 26 （2） ：299-309.

[8] 宋憲強(qiáng)， 梁釗雄， 周紅藝等. 基于決策樹與時序NDVI變化檢測的耕地撂荒遙感監(jiān)測——以四川省涼山州普格縣為例[J]. 山地學(xué)報(bào)， 2021， 39 （6） ： 912-921.

[9] 陳曦，彭鳳姣，蔡勇，等.基于ArcGIS核密度分析法的湖北省國家濕地公園空間分布特征及影響因素[J].綠色科技，2021，23（2）：1-3.

[10] 許曉婷，王思恩，焦俏，等.陜西關(guān)中平原撂荒耕地空間分布格局[J].湖北農(nóng)業(yè)科學(xué)，2023，62（3）：75-79.

收稿日期：2023-08-05

基金項(xiàng)目：青海省重點(diǎn)研發(fā)與轉(zhuǎn)化計(jì)劃項(xiàng)目（2022-QY-225）；青海民族大學(xué)創(chuàng)新項(xiàng)目（39M2023014）。

作者簡介：葉鵬帥（1998—），碩士，主要從事耕地資源監(jiān)測與評價研究。E-mail：1796220951@qq.com。

*為通信作者：E-mail：1989990003@qhu.edu.cn。