摘要 采用喀什地區(qū)2010、2015和2020年Landsat-TM影像數(shù)據(jù)，結(jié)合GIS技術(shù)及植被覆蓋度計算方法，分析2010—2020年喀什地區(qū)土地利用變化與植被覆蓋度變化特征。結(jié)果表明，2010—2020年喀什地區(qū)土地利用格局主要表現(xiàn)為耕地和林地面積的持續(xù)增加以及草地和未利用地面積的不斷減少；分析土地類型之間的相互轉(zhuǎn)化發(fā)現(xiàn)，區(qū)域內(nèi)土地轉(zhuǎn)化以耕地、林地和草地的變化尤為突出，草地面積明顯減少的主要原因是向林地、耕地和未利用地發(fā)生轉(zhuǎn)化，林地面積的增加主要源于未利用地、耕地和草地的轉(zhuǎn)化；2010—2020年喀什地區(qū)植被覆蓋度總體變化幅度較小且主要以中低覆蓋度為主。

關(guān)鍵詞 土地利用；植被覆蓋；變化特征；喀什地區(qū)

中圖分類號 F301.2 文獻標識碼 A 文章編號 0517-6611（2023）09-0053-05

doi：10.3969/j.issn.0517-6611.2023.09.014

Abstract Using Landsat TM image data from Kashgar region in 2010， 2015， and 2020， combined with GIS technology and vegetation coverage calculation methods， the characteristics of land use change and vegetation coverage change in Kashgar region from 2010 to 2020 were analyzed.The results showed that the land use pattern in Kashgar region from 2010 to 2020 was mainly manifested by the continuous increase in the area of cultivated land and forest land， as well as the continuous decrease in the area of grassland and unused land.Analyzing the mutual transformation between land types， it was found that the transformation of cultivated land， forest land and grassland was particularly prominent in the region. The main reason for the significant decrease in grassland area was the transformation to forest land， cultivated land， and unused land. The increase in forest land area mainly stems from the transformation of unused land， cultivated land and grassland.From 2010 to 2020， the overall change in vegetation coverage in Kashgar region was relatively small and mainly characterized by medium to low coverage.

Key words Land use;Vegetation coverage;Change characteristics;Kashgar region

基金項目 陜西省土地工程建設(shè)集團內(nèi)部科研項目（DJNY-YB-2023-11）。

作者簡介 馬鋼（1987—），男，河南南陽人，工程師，在讀博士，從事土地工程、生態(tài)學、草學領(lǐng)域研究。

區(qū)域內(nèi)的土地利用變化會對生態(tài)環(huán)境、氣候和生物多樣性等產(chǎn)生影響，在人與自然環(huán)境的交互關(guān)系中有著重要作用［1］，是全球共同關(guān)注的熱點問題［2］。隨著溫室效應(yīng)加劇，全球變暖已成為人類所面臨的最為嚴峻的挑戰(zhàn)。社會經(jīng)濟的快速發(fā)展，人們對土地的過度開發(fā)利用，導(dǎo)致賴以生存的生態(tài)環(huán)境日益惡化，水土資源短缺、生物量大量縮減、土地荒漠化、土壤鹽漬化及生物多樣性降低等一系列生態(tài)問題頻發(fā)［3］。

生態(tài)環(huán)境是人類社會賴以生存和發(fā)展的基礎(chǔ)，同時不合理的人類活動也會反作用于生態(tài)環(huán)境，影響生態(tài)環(huán)境的組成、結(jié)構(gòu)及其功能，嚴重時則會導(dǎo)致生態(tài)環(huán)境質(zhì)量發(fā)生系統(tǒng)退化或病變［4］。人類活動作用于生態(tài)環(huán)境，最直接的方式是通過改變土地利用方式影響生態(tài)環(huán)境質(zhì)量狀況。土地利用的變化會引起區(qū)域內(nèi)物質(zhì)能量循環(huán)發(fā)生巨大的變化，給區(qū)域生態(tài)環(huán)境帶來重大影響，并產(chǎn)生明顯的生態(tài)環(huán)境效應(yīng)［5］。自1995年推出土地利用/土地覆被變化（LUCC）計劃以來，土地利用變化的生態(tài)環(huán)境影響就已經(jīng)成為國內(nèi)外眾多學者關(guān)注的熱點問題［6］，其研究主要包含土地利用及其驅(qū)動力［7］、空間格局［8］、土地利用變化分析［9］等。近年來，隨著RS和GIS技術(shù)的迅速發(fā)展及研究的不斷深入，許多學者用多種遙感指數(shù)對不同區(qū)域及各種類型的生態(tài)環(huán)境進行研究［10-12］，研究尺度大到洲際、小到縣域，為多尺度生態(tài)環(huán)境動態(tài)監(jiān)測提供了有效手段［13-15］，使生態(tài)環(huán)境質(zhì)量評價理論和方法逐步得到完善。植被對全球氣候變化具有十分重要的意義，屬于國際地圈生物圈計劃（IGBP）的核心內(nèi)容，植被是陸地生態(tài)系統(tǒng)的重要組成部分，對于水土保持、生態(tài)環(huán)境和土壤的改善具有積極作用，植被覆蓋指數(shù)（NDVI）是相關(guān)研究中最常用和最重要的指標，開展關(guān)于NDVI的研究對全球變化和生態(tài)環(huán)境的可持續(xù)發(fā)展十分重要［16-17］。

綠洲在干旱區(qū)尤為重要，同樣也最為復(fù)雜和敏感，人類在干旱區(qū)開展的生產(chǎn)生活等一系列活動主要集中在綠洲內(nèi)。綠洲區(qū)域內(nèi)的土地利用對當?shù)厣鐣?jīng)濟的發(fā)展具有積極意義［18］。近年來已經(jīng)開展了大量與土地利用相關(guān)的研究，杜宏茹［19］對我國綠洲的土地利用空間變化及其驅(qū)動因素進行了研究；陳瓊等［20］從土地利用現(xiàn)狀、變化及效應(yīng)等方面對黃河流域河源區(qū)進行了全面分析；孔冬冬等［21］探究了青藏高原植被的NDVI空間格局變化及其物候變化；李婭麗等［22］基于NDVI的MODIS 13Q1數(shù)據(jù)產(chǎn)品，對區(qū)域內(nèi)的地表溫度與植被覆蓋之間的關(guān)系進行研究。該研究通過2010、2015、2020年3期遙感影像提取植被覆蓋數(shù)據(jù)，對喀什地區(qū)土地利用格局的演變和區(qū)域植被覆蓋度進行分析，為干旱區(qū)綠洲水土資源的可持續(xù)利用和生態(tài)環(huán)境的發(fā)展建設(shè)提供技術(shù)支撐和科學依據(jù)。

1 資料與方法

1.1 研究區(qū)概況

研究區(qū)位于新疆維吾爾自治區(qū)西部（35°20′～40°18′N、73°20′～79°57′E），地處我國中緯度大陸腹地、西北內(nèi)陸干旱區(qū)，東臨塔克拉瑪干大沙漠，東北與阿克蘇地區(qū)柯坪縣、阿瓦提縣相連（圖1）。該區(qū)面積達16.2萬km2，邊境線總長888 km，南北長535 km，東西寬750 km［23］。

該區(qū)氣候主要為高原干旱、半干旱大陸性氣候，高寒、干旱，日照時間長，太陽輻射強；雨熱同季，降水量小，蒸發(fā)量大；日溫差大，年溫差大；氣候地理分布差異大，氣候特征垂直變化較為明顯［24］。植被資源較為貧乏，主要以草本植物圓葉鹽爪爪（Kalidium schrenkianum）和絹蒿（Artemisiaka schgaricum）居多，而白刺（Nitraria tangutorum Bobr）、小蓬（Nanophyton erinaceum）、狗尾草（Setaria viridis）、鹽生草（Halogeton glomeratus）分布較少，灌木植物主要以檉柳（Tamarix chinensis Lour）和鹽穗木（Halostachys belangeriana）為主，植被特征整體較為單一，甚至存在單個樣地無樣本現(xiàn)象［25］。

1.2 研究方法

1.2.1 土地利用分析。根據(jù)全國土地利用/覆蓋分類系統(tǒng)，參考Tan等［26］的研究方法對新疆維吾爾自治區(qū)喀什地區(qū)土地利用進行遙感解譯，結(jié)合流域Landsat-TM影像數(shù)據(jù)的空間分辨率及流域特征，確定區(qū)分出6個一級分類和19個小二級分類。一級分類主要分為耕地、林地、草地、水域、建設(shè)用地和未利用土地，利用ArcGIS 10.7中的空間分析工具對影像進行轉(zhuǎn)移矩陣分析。

1.2.2 植被覆蓋度分析。該研究中植被覆蓋度指數(shù)的提取所用到的數(shù)據(jù)為MOD13Q1產(chǎn)品中的EVI數(shù)據(jù)，空間分辨率為250 m×250 m，時間序列為2010年1月—2020年12月。對獲得的EVI遙感數(shù)據(jù)進行數(shù)據(jù)格式轉(zhuǎn)換、鑲嵌、投影轉(zhuǎn)換及研究區(qū)提取等預(yù)處理，為降低噪音信息對數(shù)據(jù)影響，還對其進行Savitzky-Golay濾波和MVC合成處理，獲得代表植被生長最好狀況的年EVI數(shù)據(jù)［27］。

像元二分模型是目前植被覆蓋度反演的有效方法，其計算公式如下：

Fc=EVI-EVIsoilEVI-EVIveg×100%（1）

式中，F(xiàn)c為植被覆蓋度；EVIsoil為研究區(qū)純裸土像元EVI值；EVIveg為純植被像元EVI值。該研究分別取研究區(qū)EVI圖像直方圖的5%處和95%處EVI值代表EVIsoil值和EVIveg值；根據(jù)Fc對植被覆蓋度進行劃分：0≤Fc≤25%，低覆蓋區(qū)；25%＜Fc≤50%，中低覆蓋區(qū)；50%＜Fc≤75%，中覆蓋區(qū)；Fc＞75%，高覆蓋區(qū)。

1.3 數(shù)據(jù)分析

該研究采用ENVI 5.3軟件對遙感數(shù)據(jù)進行格式轉(zhuǎn)換、鑲嵌、投影轉(zhuǎn)換及研究區(qū)提取等處理，對質(zhì)量較差數(shù)據(jù)以及缺失的像元進行平滑等預(yù)處理并結(jié)合研究區(qū)特有的生態(tài)環(huán)境特點［28］，利用ArcGIS 10.7軟件對土地利用數(shù)據(jù)和植被覆蓋度數(shù)據(jù)進行重分類分析［29］；利用Microsoft excel 2019對數(shù)據(jù)進行統(tǒng)計與處理［30］，利用Origin 2017分析其變化趨勢并繪圖［31］。

2 結(jié)果與分析

2.1 土地利用變化時空分布格局

從圖2可以看出，2010—2020年喀什地區(qū)土地利用發(fā)生顯著變化。結(jié)合10年來土地利用面積變化（表1）可以發(fā)現(xiàn)，喀什地區(qū)的草地和未利用地面積占很大比例，表現(xiàn)主要的特征為耕地和林地面積的持續(xù)增加以及草地和未利用地面積的不斷減少。

由表1可知，2010—2015年水域和建設(shè)用地的面積變化并不明顯，耕地面積增加了0.11萬km2，相比于其他土地類型，草地和未利用地的面積有所減少。2015—2020年，林地的面積有了明顯的增加，至2020年增加了22.97%；草地和未利用地的面積卻持續(xù)減少。綜合這2個時期的土地利用面積變化分析，2010—2020年新增耕地面積為0.14萬 km2；林地面積新增0.40萬km2，增長率高達28.17%，主要集中于喀什地區(qū)的中部區(qū)域；未利用地面積減少了3.67%；草地面積減少了0.32萬km2，與2010年相比減少了6.96%；水域和建設(shè)用地在此期間面積發(fā)生緩慢增加。

2.2 土地類型之間的相互轉(zhuǎn)化

從2010—2020年的土地類型轉(zhuǎn)移矩陣（表2～3）可以看出，2010—2020年研究區(qū)域的土地利用格局發(fā)生顯著變化，其中耕地、林地和草地的變化尤為突出。2010—2015年，耕地與其他土地類型之間的轉(zhuǎn)化十分明顯，耕地增加的來源主要是林地、草地，轉(zhuǎn)化比例分別為11.41%、3.57%，耕地轉(zhuǎn)出方向也主要是林地和草地；林地增加的來源主要是耕地、草地、未利用地，轉(zhuǎn)化比例分別為11.06%、14.10%、27.80%，林地轉(zhuǎn)出方向主要是耕地、草地、未利用地；草地增加的來源最主要是未利用地，轉(zhuǎn)化比例為6.96%，但同時轉(zhuǎn)出方向也主要是未利用地，且總體上草地是減少趨勢；水域增加的主要來源是草地和未利用地，轉(zhuǎn)出方向也主要是草地和未利用地，整體上水域僅增加了0.84%；建設(shè)用地增加了36.71%；未利用地主要在林地、草地、水域3個土地類型間轉(zhuǎn)移，整體是減少趨勢。2015—2020年，耕地增加的主要來源還是林地和草地，主要轉(zhuǎn)出方向也是林地和草地；與2010—2015年相同，林地增加的來源主要還是耕地、草地和未利用地，轉(zhuǎn)出方向也一致，整體上轉(zhuǎn)入明顯大于轉(zhuǎn)出，導(dǎo)致林地明顯增加；草地增加和轉(zhuǎn)出方向均主要為未利用地，整體是減少趨勢；水域增加主要來自草地和未利用地，也主要向草地和未利用地轉(zhuǎn)化，基本保持動態(tài)平衡；建設(shè)用地增加了18.52%，與2010—2015年相比，建設(shè)用地增加幅度減少；整體上，未利用地不斷被利用為林地，減少了2.57%。

2.3 植被覆蓋度變化 由圖3和圖4可知，2010—2020年，喀什地區(qū)植被覆蓋度總體變化幅度較小且主要以中低覆蓋度為主，低覆蓋度區(qū)域面積占比表現(xiàn)為2015年最大，中低覆蓋度區(qū)域面積占比表現(xiàn)為2020年最大，高覆蓋度和中覆蓋度區(qū)域面積占比均表現(xiàn)為2010年最大，且高覆蓋度和中覆蓋度區(qū)域主要分布于人工綠洲區(qū)域。說明2010—2020年喀什地區(qū)生態(tài)環(huán)境狀況相對較為穩(wěn)定，但有所下降。

3 結(jié)論

該研究采用喀什地區(qū)2010、2015和2020年Landsat-TM影像數(shù)據(jù)，結(jié)合GIS技術(shù)及植被覆蓋度指數(shù)計算方法，分析2010—2020年喀什地區(qū)土地利用變化與植被覆蓋度變化特征。結(jié)果表明，2010—2020年喀什地區(qū)土地利用格局主要表現(xiàn)為耕地和林地面積的持續(xù)增加以及草地和未利用地面積的不斷減少；分析土地類型之間的相互轉(zhuǎn)化發(fā)現(xiàn)，區(qū)域內(nèi)土地轉(zhuǎn)化以耕地、林地和草地的變化尤為突出，草地面積明顯減少的主要原因是向林地、耕地和未利用地發(fā)生轉(zhuǎn)化，林地面積的增加主要源于未利用地、耕地和草地的轉(zhuǎn)化；2010—2020年喀什地區(qū)植被覆蓋度總體變化幅度較小且主要以中低覆蓋度為主。由此可見，近年來喀什地區(qū)生態(tài)環(huán)境狀況有所下降。

參考文獻

［1］ 吳建國，呂佳佳.土地利用變化對生物多樣性的影響［J］.生態(tài)環(huán)境，2008，17（3）：1276-1281.

［2］ 丘君，陳利頂，傅伯杰.土地利用/覆被變化對生物多樣性的影響［C］//中國地理學會自然地理專業(yè)委員會.土地覆被變化及其環(huán)境效應(yīng)學術(shù)研討會論文集.北京：星球地圖出版社，2002：55-63.

［3］ 曹建廷，李原園.虛擬水及其對解決我國水資源短缺問題的啟示［J］.科技導(dǎo)報，2004（3）：15-16，14.

［4］ 馬世駿，王如松.社會-經(jīng)濟-自然復(fù)合生態(tài)系統(tǒng)［J］.生態(tài)學報，1984，4（1）：1-9.

［5］ 鐘海燕.鄱陽湖區(qū)土地利用變化及其生態(tài)環(huán)境效應(yīng)研究［D］.南京：南京農(nóng)業(yè)大學，2011.

［6］ TURNER B L II，SKOLE D，SANDERSON S，et al.Landuse and landcover change：Science/research planning［J］.Global change report（Sweden），1995，43：669-679.

［7］ 董建紅，張志斌，笪曉軍，等.“三生”空間視角下土地利用轉(zhuǎn)型的生態(tài)環(huán)境效應(yīng)及驅(qū)動力：以甘肅省為例［J］.生態(tài)學報，2021，41（15）：5919-5928.

［8］ 王少華.鄭州沿黃旅游區(qū)土地利用變化及其生態(tài)環(huán)境效應(yīng)評價研究［D］.開封：河南大學，2016.

［9］ 伏洋，肖建設(shè)，校瑞香，等.基于RS和GIS的西寧市LUCC分析及模擬預(yù)測［J］.農(nóng)業(yè)工程學報，2009，25（7）：211-218.

［10］ 高志強，劉紀遠，莊大方.基于遙感和GIS的中國土地利用/土地覆蓋的現(xiàn)狀研究［J］.遙感學報，1999，3（2）：134-138.

［11］ 何春陽，陳晉，陳云浩，等.土地利用/覆蓋變化混合動態(tài)監(jiān)測方法研究［J］.自然資源學報，2001，16（3）：255-262.

［12］ 史培軍，宮鵬，李曉兵，等.土地利用/覆蓋變化研究的方法與實踐［M］.北京：科學出版社，2000.

［13］ 李秀彬.全球環(huán)境變化研究的核心領(lǐng)域：土地利用/土地覆被變化的國際研究動向［J］.地理學報，1996，51（6）：553-558.

［14］ 劉盛和，吳傳鈞，沈洪泉.基于GIS的北京城市土地利用擴展模式［J］.地理學報，2000，55（4）：407-416.

［15］ 賀秋華.江蘇濱海土地利用/覆蓋變化及其生態(tài)環(huán)境效應(yīng)研究［D］.南京：南京師范大學，2011.

［16］ 高彥春，王長耀.水文循環(huán)的生物圈方面（BAHC計劃）研究進展［J］.地理科學進展，2000，19（2）：97-103.

［17］ 王曉利.中國沿海極端氣候變化及其對NDVI的影響特征研究［D］.煙臺：中國科學院煙臺海岸帶研究所，2017.

［18］ 石培基，周堂，張寧.綠洲土地利用的景觀生態(tài)學分析：以武威市古浪縣綠洲為例［J］.干旱區(qū)資源與環(huán)境，2011，25（8）：53-58.

［19］ 杜宏茹.生態(tài)城市建設(shè)進程中土地利用變化趨勢［J］.國土與自然資源研究，2004（1）：9-10.

［20］ 陳瓊，張鐿鋰，劉峰貴，等.黃河流域河源區(qū)土地利用變化及其影響研究綜述［J］.資源科學，2020，42（3）：446-459.

［21］ 孔冬冬，張強，黃文琳，等.1982—2013年青藏高原植被物候變化及氣象因素影響［J］.地理學報，2017，72（1）：39-52.

［22］ 李婭麗，汪小欽，陳蕓芝，等.福建省地表溫度與植被覆蓋度的相關(guān)性分析［J］.地球信息科學學報，2019，21（3）：445-454.

［23］ 阿布都克日木·阿巴司，秦榕，伊力達爾江·吐生，等.喀什地區(qū)1961—2010年氣候變化特征［J］.沙漠與綠洲氣象，2012，6（6）：34-40.

［24］ 熱汗古麗·吾買爾，滿蘇爾·沙比提，陸吐布拉·依明.喀什地區(qū)近10年地下水資源時空動態(tài)變化分析［J］.干旱區(qū)資源與環(huán)境，2011，25（7）：63-68.

［25］ 王文杰，張永福．基于GIS的新疆喀什地區(qū)耕地變化及驅(qū)動力研究［J］.湖北農(nóng)業(yè)科學，2016，55（17）：4560-4564．

［26］ TAN K L，WANG X F，GAO H J.Analysis of ecological effects of comprehensive treatment in the Tarim River Basin using remote sensing data［J］.Mining science and technology（China），2011，21（4）：519-524.

［27］ 徐長春，陳亞寧，李衛(wèi)紅，等.塔里木河流域近50年氣候變化及其水文過程響應(yīng)［J］.科學通報，2006，51（S1）：21-30.

［28］ 李孟倩，蘇世宗.基于Landsat8-OLI遙感影像水稻種植區(qū)提取方法研究［J］.華北理工大學學報（自然科學版），2019，41（3）：26-33，59.

［29］ 靳士科.遙感技術(shù)在核電陸域生態(tài)環(huán)境狀況評價中的應(yīng)用研究［J］.環(huán)境科學與管理，2020，45（12）：163-166.

［30］ 謝勇堅.應(yīng)用Microsoft Excel對分析數(shù)據(jù)進行處理［J］.西南給排水，2001，23（4）：5-12.

［31］ 單文坡，盧海霞，鄭輝.Excel和Origin在環(huán)境監(jiān)測數(shù)據(jù)處理中的應(yīng)用［J］.石家莊職業(yè)技術(shù)學院學報，2008，20（2）：58-60，67.