退耕前后20年延安市土地利用變化及驅動因素分析

王曉飛，孫智輝，劉志超，曹雪梅

（延安市氣象局，陜西 延安 716000）

土地利用與覆蓋變化（LUCC）是全球氣候、環境變化研究的重要內容之一［1，2］，成為當前學術界研究的熱點。改革開放以來中國經濟快速發展，工業化、城市化進程加快，使得LUCC 發生了顯著變化，特別是在黃土高原，1999 年以后實施了包括天然林保護工程（Natural forest protection program，NFPP）和退耕還林/還草工程（Grain for green program，GGP）等，將對植被變化、水資源安全以及可持續發展產生深遠影響。

延安市地處黃土高原，為中國半干旱半濕潤過渡帶，對氣候變化［3］和人為干擾［4］較為敏感。1999年開始實施了天然林保護工程、退耕還林工程和“三北”防護林工程，至2020 年底延安市實施退耕還林面積71.83 萬hm2，占全市國土總面積的19.4%，占全國計劃內退耕還林面積的2.5%，占全省計劃內退耕還林面積的27.0%。工程涉及28.6 萬農戶、124.8 萬農村人口。延安市完成營造林142.27 萬hm2，林業用地面積增至289.24 萬hm2，占國土面積的78.2%，有林地面積增至171.47 萬hm2，森林覆蓋率達53.07%，成為陜西省“一葉肺”［5］。針對延安市土地利用、植被覆蓋度變化、生態效益及其對濕度、地表溫度、徑流的影響，以及植被對氣候的響應等，國內學者在退耕還林工程實施早期多利用退耕后遙感數據或低分辨率數據進行研究［6-19］。延安市在1999 年全面開展退耕還林工程［20，21］，退耕以后林草的生長還需要很長的過程。因此，本研究將2000 年作為分界線，利用延安市1980、1990、2000 年遙感影像數據反映退耕前土地利用變化，利用2005、2010、2015、2020 年遙感數據反映退耕后土地利用變化，采用1981—2020 年延安市降水、溫度等氣候數據及城鎮化、人均GDP 等社會發展經濟數據，全面分析延安市土地利用變化及其影響因素，以期為該市生態建設提供參考。

1 研究區概況與數據

1.1 研究區概況

延安市轄1 市2 區10縣，國土面積3.7 萬km2，人口226 萬人。延安市年降水量520 mm，年平均氣溫9.7 ℃，年日照時數2 400 h，年蒸發量1 450～1 700 mm，是降水量的3 倍左右。根據延安市地形、氣候特點，可分為南部塬區和北部丘陵溝壑區，北部為半干旱區，南部為半濕潤區。延安市北部土壤侵蝕嚴重，主要實施了退耕還林和“三北”防護林工程，因此，本研究將延安市作為一個整體進行分析，同時對延安市北部單獨進行土地利用分析。延安市北部包括吳起、志丹、安塞、子長、延川、延長和寶塔7 縣（區、市），總面積20 837 km2，總人口134.96 萬人。

1.2 數據

遙感數據來源于中國科學院資源環境科學與數據中心數據平臺（http：//www.resdc.cn），人口數據來源于全國第七次人口普查公告［22］，經濟數據來源于《2020 年延安市國民經濟和社會發展統計公報》［23］。氣象數據由延安市各縣（區、市）氣象站提供。

2 方法

2.1 分類方法

全國土地利用類型遙感監測空間分布數據（http：//www.resdc.cn/Datalist1.aspx？FieldTyepID=1，3）是基于美國陸地衛星Landsat TM 影像，通過人工目視解譯生成。土地利用類型包括耕地、林地、草地、水域、居民地和未利用土地6 個一級類型以及25個二級類型。數據經過嚴格質量控制，其中土地利用的一級類型綜合評價精度達94.3%，二級類型分類精度達91.2%以上［24］。李艷忠等［25］對黃河中游地區進行精度驗證，其中9 種類別各選取了50 個樣本點，驗證整體精度達92.9%，滿足分析要求。延安市位于黃河中游，因此本研究沒有進行精度驗證，對分類結果重新合并處理，共獲取了包括耕地、林地、草地、未利用地、建設用地和水域6 類數據。林地和草地按二級分類結果進行分析，林地包括有林地、灌木林、疏林地、其他林地，草地包括高覆蓋草地、中覆蓋草地、低覆蓋草地。

2.2 土地利用轉移矩陣

國土空間轉移矩陣用來刻畫區域空間類型之間相互轉換的數量關系，反映在人類活動作用下國土空間開發和保護的方向及結構，揭示國土空間格局的變化特征［26］。計算公式如下。

式中，aij為研究初期i類國土空間轉為末期j類國土空間的總面積；n為國土空間利用類型數量。本研究運用ArcGIS 10.5 軟件的空間分析模塊，對2020 年和2000 年土地利用數據進行統計和疊加分析，得到退耕后20 年延安市土地利用類型面積轉移矩陣，通過矩陣分析各時段內各地類之間的數量轉換關系和土地利用動態變化情況。

2 結果與分析

2.1 延安市土地利用現狀

由圖1 可以看出，延安市區域內的土地類型以草地為主，占總面積的46.7%，其次是耕地面積和林地，分別占總面積的25.2%、26.8%。水域、建設用地、未利用土地面積占比為1.3%。林草面積占比達73.5%。延安市北部草地、耕地、林地面積占比分別為51.8%、30.6%、16.7%。表明1980—2020 年耕地面積減少，林地、草地面積增加主要是在國家實行退耕還林政策后發生了明顯變化。

2.2 延安市土地利用變化

2.2.1 耕地變化 由表1、表2 可知，延安市1980—2020 年耕地面積減少了2 094.61 km2，在退耕還林前，1980—2000 年耕地面積緩慢增加，增加了327.24 km2，在2000 年達到峰值。退耕還林后至2020 年耕地面積減少了2 421.84 km2。隨著退耕還林工程的實施，耕地面積減少主要集中在2000—2010年，減少了2 279.11 km2，占總退耕面積的94.1%。延安市北部1980—2020 年耕地面積減少了2 040.51 km2，在退耕還林前，1980—2000 年耕地面積增加了46.00 km2，退耕還林后至2020 年耕地面積減少了2 086.51 km2。延安市北部退耕面積占全市退耕面積的86.15%。

2.2.2 林地變化 延安市林地面積1980 年以來呈增加趨勢，至2020 年增加了923.44 km2，增加時段主要發生在2000—2020年，增加了875.41 km2，占林地總增加面積的94.8%（表1）。延安市北部林地面積1980—2020 年增加了858.11 km2，退耕還林 后至2020 年林地面積增加了825.89 km2，占北部林地總增加面積的96.2%（表2）。表明林地增加面積主要集中在北部，占延安市林地總增加面積的92.9%。

延安市林地二級分類有林地、灌木林、疏林地、其他林地，其面積變化如下，有林地面積保持基本穩定，其他3 類林地面積分別增加了233.16、95.21、608.10 km2，以其他林地面積變化最為明顯，占林地增加面積的65.9%，其次是灌木林，其增加面積占林地增加面積的25.2%。

灌木林和其他林地增加面積占林地增加面積的91.1%，因此只分析灌木林和其他林地面積的變化時段。灌木林面積增加主要集中在2000—2010年，而其他林地面積增加集中在1995—2020年，灌木林面積變化與退耕還林工程有關，而其他林地與農業結構調整有關。

2020 年延安市北部有林地、灌木林、疏林地、其他林地面積占其國土面積比重分別為0.1%、9.0%、3.6%和3.1%，以灌木林面積最大，有林地面積最小。

2.2.3 草地變化 延安市1980—2020 年草地面積增加了13 125.66 km2，但在退耕還林前草地面積減少了236.36 km2，退耕還林后草地面積增加了1 362.02 km2（表1）。草地面積變化主要是由于高、中覆蓋草地面積增加明顯，高覆蓋草地面積由2000年的1 905.46 km2增加到2020 年的2 244.09 km2；中覆蓋草地面積由2000年的10 302.38 km2增加到2020年的11 488.79 km2，分別增加了338.63、1 186.41 km2。同樣，退耕還林后草地面積增加也主要反映在延安市北部，面積增加了1 099.32 km2，占全市退耕還林后草地增加面積的80.7%（表2）。

表2 退耕前后20 年延安市北部土地利用變化情況 （單位：km2）

2020 年延安市北部高覆蓋、中覆蓋、低覆蓋草地面積分別為410.67、6782.36、3 536.71 km2，分別占延安市北部國土面積的19.8%、32.7%、17.1%，草地以中低覆蓋為主。草地面積變化主要是高、中覆蓋草地面積增加明顯，分別增加了194.37、919.83 km2。低覆蓋草地面積2020 年比1980 年略有減少，減少了19.97 km2。高覆蓋草地面積增加主要發生在2000—2010年，而中覆蓋草地面積變化集中在2005—2010 年。

2.2.4 建設用地、未利用地、水域 延安市建設用地面積呈持續增加趨勢，由1980 年的116.16 km2增加到2020 年的280.88 km2，增幅達142%。變化明顯時段主要發生在2005—2010 年和2015—2020年，分別增加了49.28 km2和74.79 km2（表1）。

未利用地面積在2010 年前穩定少變，在2010 年以后有所增加，特別是2015—2020年，面積增加了26.73 km2（表1）。

水域面積在1980—2010 年有少量減少，2010 年后開始增加，2015—2020 年增加了14.34 km2，水域面積正在緩慢增長（表1）。

2.2.5 空間轉移特征 基于GIS 空間分析方法和土地利用分類數據，獲取研究區2000—2020 年土地利用變化格局［26，27］。由表3 可知，2000—2020 年研究區域耕地面積總體上呈減少趨勢，耕地轉出面積為2 772.84 km2，主要轉出類型為草地、林地、建設用地，轉出面積分別為1 953.13、702.50、106.10 km2，另轉出水域面積6.48 km2和未利用地面積4.63 km2。林地面積增加，主要由耕地和草地轉入，面積分別為702.50、342.43 km2。草地面積增加明顯，主要由耕地、林地轉入，轉入面積分別為1 953.13、87.31 km2。建設用地面積也有所增加，主要由耕地、草地和林地轉入。

表3 2000—2020 年延安市土地利用轉移矩陣 （單位：km2）

如圖2 所示，土地轉化為林草地主要發生在延安市北部，草地分布比較廣泛，林地較為集中。轉化為林地的區域發生在吳起縣東北、志丹縣中部和東北、安塞區北部、子長市西部和清澗河兩岸，延川縣中部和延長縣中部，寶塔區東部等地，特別是在子長市、安塞區、志丹縣、吳起縣等地北部邊界，處于400 mm 降水線波動影響范圍內［28，29］，形成了林帶。而在吳起縣西北年降水量少于400 mm，耕地轉為林地面積很少。建設用地面積增加多集中在延安市城區、縣城區和北部重點鄉鎮，南部塬區農村建設用地面積也有增加。

圖2 2000—2020 年延安市土地利用轉置情況

3 驅動因素分析

3.1 氣候

3.1.1 降水 如圖3 所示，1981 年以來，延安市年降水量呈波動增加的趨勢，氣候傾向率為16.138 mm/10 年。延安市1981—2020 年平均降水量為515.5 mm，20 世紀80、90 年代，21 世紀00、10 年代的年降水量分別為525.7、443.3、520.4、572.6 mm，退耕以來的20 年降水量偏多，特別是21 世紀10 年代降水量明顯偏多。

圖3 1981—2020 年延安市降水量變化趨勢

退耕還林前后20 年延安市平均降水量比較（2001—2020 年平均降水量減去1981—2000 年平均降水量）如圖4 所示，退耕后降水量偏多62.6 mm，增加值在延安市南部區域較少，高于平均值的區域主要分布在延安市中北部。子長市、延川縣降水量差值達80 mm 以上，寶塔區、志丹縣降水量差值超過70 mm，安塞區、延長縣降水量差值分別為69.7 mm和65.2 mm，均在平均值以上，吳起縣較少，盡管降水量差值為50 mm，其對林草生長也十分有利。

圖4 延安市退耕后與退耕前降水量差值分布

聶桐等［30］通過對延安市植被覆蓋度時空變化及驅動力研究，發現降水量、溫度是植被覆蓋及其動態變化的重要影響因素，當降水量為547～554 mm，氣溫為10.49～10.64 ℃時，植被覆蓋好，適合植被生長。如圖5 所示，延安市降水量空間分布上有很大差異，南部降水量多，而西北和東部黃河沿岸降水量少。統計發現退耕以來，延安市北部20 年平均降水量為523.3 mm，安塞區、寶塔區降水量達550 mm；特別是21 世紀10 年代，延安市北部降水量達555.7 mm，除吳起縣降水量僅為500 mm外，其他縣區降水量均達550 mm 以上，特別是寶塔區降水量超過600 mm，充沛的降水量為植被生長提供了良好的水分條件。

圖5 2001—2020 年延安市平均降水量分布

3.1.2 溫度 由圖6 可以看出，1981 年以來，年平均溫度在1981—1993 年快速上升，1994—2020 年溫度變化平緩。從年代際變化來看，20 世紀80、90 年代，21 世紀00、10 年代的年平均溫度分別為9.1、9.7、10.1、10.0 ℃。退耕以來延安市處于相對暖期。

圖6 1981—2020 年延安市平均溫度變化趨勢

由圖7 可以看出，在空間分布上，延安市溫度呈西北低、東部黃河沿岸和南部溫度高的特點，吳起縣、志丹縣溫度低于9.0 ℃，對植被生長有一定的影響。

圖7 2001—2020 年延安市平均溫度分布

3.2 延安市人口

根據第七次全國人口普查公報，延安市總人口從1982 年的153.3 萬人增加到2020 年的228.3 萬人，總體呈增加趨勢，但增長率和人口增加數量呈降低趨勢。城市人口增加非常快，從18.9 萬人增加到140.1 萬人，而農村人口在1982—2000 年始終維持在較高水平，明顯多于城鎮人口，1982—2010 年農村人口均多于城鎮人口。2000—2010 年人口數量變化較大，鄉村人口大幅減少，城鎮人口快速增加；2010—2020 年城鄉人口結構實現了逆轉，到2020 年城鎮人口遠多于鄉村人口。

表4 延安市人口變化

3.3 延安市經濟

延安市國民生產總值在1980 年為4.33 億元，到2020 年發展為1 601.5 億元，實現了飛躍式發展。但從年變化來看，在2004 年以后進入高速發展期，2006 年國民生產總值突破500 億元，2011 年突破千億元大關，2018 年超過1 500 億元。城鄉居民人均可支配收入在2005—2015 年有了明顯增加，城鎮居民人均可支配收入由7 471 元增加到33 127元，鄉村人均可支配收入由2 195 元增加到9 198元，并在2016 年突破萬元。

3.4 延安市政策因素

制定地權制度、經營機制、社會經濟政策等是區域政策和土地政策法規影響土地利用結構的直接方式，也是影響土地利用方式的重要因素。20 世紀80年代實行包產到戶，提高了勞動生產率，隨后實施農業產業結構調整，發展林果、草畜、設施農業等高效產業，大力發展二三產業，城市化水平以及社會經濟發展水平得到較快提高。同時，積極開展退耕還林工程和治溝造地土地整治重大工程，從而增加農民收入，進而推動區域內的社會經濟發展。

4 小結

1）延安市土地類型以草地為主，占總面積的46.7%，其次是耕地和林地，分別占總面積的25.2%、26.8%。水域、建設用地、未利用土地占比1.3%。1980—2000 年延安市耕地面積增加了327.24 km2，退耕還林后減少了2 421.84 km2。耕地變化主要集中在延安市北部，占全市退耕面積的86.15%。

2）林地面積增加了923.43 km2，增加時段主要在2000 年以后，增加了875.4 km2，占比達94.8%。林地增加面積主要集中在延安市北部，占林地總增加面積的92.9%。林地二級分類面積變化以其他林地最為明顯，占林地增加面積的65.9%，其次是灌木林，增加面積占比為25.2%。轉化為林地的區域在吳起縣東北部、志丹縣中部和東北部、安塞區北部、子長市西部和清澗河兩岸，延川縣中部和延長縣中部，寶塔區東部等地。受氣候影響，在延安市北部邊界400 mm 降水帶，出現耕地轉化為林地的地帶，生態得到明顯改善。

3）1980—2020 年草地面積增加了1 125.66 km2，1980—2000 年草地面積減少了236.36 km2，退耕還林后草地面積增加了1 362.02 km2。草地面積變化主要是因為高、中覆蓋草地面積增加明顯，分別增加了338.63、1 186.41 km2，但還是以中低覆蓋草地為主。

4）建設用地呈持續增加趨勢，由1980 年的116.16 km2增加到2020 年的280.88 km2，增加幅度達142%。明顯變化時段主要發生在2005—2010 年和2015—2020年，分別增加了49.28 km2和74.79 km2。建設用地增加多集中在延安市城區、縣城區和北部少量重點鄉鎮，南部塬區農村建設用地面積也有增加。延安市北部人口多集中到縣城和延安市城區，鄉鎮和農村居住人口大幅減少。

5）20 世紀90 年代，受干旱氣候和人口增加雙重影響，延安市出現毀草種地現象。實行退耕還林政策后至2010年，耕地面積大幅減少，政策對土地轉換起到決定性作用。隨著蘋果、設施農業發展，農民收入增加，減少了對土地的依賴，經濟水平和城鎮率提高，農民走向城市生活，確保了退耕成果。但在2015年，農村人口的減少造成了未利用土地的大幅增加，耕地沒有得到充分利用。氣候對土地利用變化也有明顯影響，2000 年以來，延安市處在降水多、溫度高的暖濕時段，特別是北部降水量，增加了70～80 mm，是林地、高中覆蓋草地面積增加的主要原因。

延安市北部中低覆蓋草地面積較大，生態比較脆弱，加之氣候還在持續變暖，氣候不確定性增加，嚴重制約著延安市的生態恢復。因此，為實現延安市生態的可持續發展和鄉村振興，還需要持續關注，開展科學研究，提供對策建議。