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

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

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

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

1 研究區(qū)概況與數(shù)據(jù)

1.1 研究區(qū)概況

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

1.2 數(shù)據(jù)

遙感數(shù)據(jù)來源于中國科學(xué)院資源環(huán)境科學(xué)與數(shù)據(jù)中心數(shù)據(jù)平臺(tái)（http：//www.resdc.cn），人口數(shù)據(jù)來源于全國第七次人口普查公告［22］，經(jīng)濟(jì)數(shù)據(jù)來源于《2020 年延安市國民經(jīng)濟(jì)和社會(huì)發(fā)展統(tǒng)計(jì)公報(bào)》［23］。氣象數(shù)據(jù)由延安市各縣（區(qū)、市）氣象站提供。

2 方法

2.1 分類方法

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

2.2 土地利用轉(zhuǎn)移矩陣

國土空間轉(zhuǎn)移矩陣用來刻畫區(qū)域空間類型之間相互轉(zhuǎn)換的數(shù)量關(guān)系，反映在人類活動(dòng)作用下國土空間開發(fā)和保護(hù)的方向及結(jié)構(gòu)，揭示國土空間格局的變化特征［26］。計(jì)算公式如下。

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

2 結(jié)果與分析

2.1 延安市土地利用現(xiàn)狀

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

2.2 延安市土地利用變化

2.2.1 耕地變化 由表1、表2 可知，延安市1980—2020 年耕地面積減少了2 094.61 km2，在退耕還林前，1980—2000 年耕地面積緩慢增加，增加了327.24 km2，在2000 年達(dá)到峰值。退耕還林后至2020 年耕地面積減少了2 421.84 km2。隨著退耕還林工程的實(shí)施，耕地面積減少主要集中在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 年以來呈增加趨勢(shì)，至2020 年增加了923.44 km2，增加時(shí)段主要發(fā)生在2000—2020年，增加了875.41 km2，占林地總增加面積的94.8%（表1）。延安市北部林地面積1980—2020 年增加了858.11 km2，退耕還林 后至2020 年林地面積增加了825.89 km2，占北部林地總增加面積的96.2%（表2）。表明林地增加面積主要集中在北部，占延安市林地總增加面積的92.9%。

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

灌木林和其他林地增加面積占林地增加面積的91.1%，因此只分析灌木林和其他林地面積的變化時(shí)段。灌木林面積增加主要集中在2000—2010年，而其他林地面積增加集中在1995—2020年，灌木林面積變化與退耕還林工程有關(guān)，而其他林地與農(nóng)業(yè)結(jié)構(gòu)調(diào)整有關(guān)。

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。高覆蓋草地面積增加主要發(fā)生在2000—2010年，而中覆蓋草地面積變化集中在2005—2010 年。

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

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

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

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

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

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

圖2 2000—2020 年延安市土地利用轉(zhuǎn)置情況

3 驅(qū)動(dòng)因素分析

3.1 氣候

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

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

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

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

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

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

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

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

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

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

3.2 延安市人口

根據(jù)第七次全國人口普查公報(bào)，延安市總?cè)丝趶?982 年的153.3 萬人增加到2020 年的228.3 萬人，總體呈增加趨勢(shì)，但增長率和人口增加數(shù)量呈降低趨勢(shì)。城市人口增加非?？?，從18.9 萬人增加到140.1 萬人，而農(nóng)村人口在1982—2000 年始終維持在較高水平，明顯多于城鎮(zhèn)人口，1982—2010 年農(nóng)村人口均多于城鎮(zhèn)人口。2000—2010 年人口數(shù)量變化較大，鄉(xiāng)村人口大幅減少，城鎮(zhèn)人口快速增加；2010—2020 年城鄉(xiāng)人口結(jié)構(gòu)實(shí)現(xiàn)了逆轉(zhuǎn)，到2020 年城鎮(zhèn)人口遠(yuǎn)多于鄉(xiāng)村人口。

表4 延安市人口變化

3.3 延安市經(jīng)濟(jì)

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

3.4 延安市政策因素

制定地權(quán)制度、經(jīng)營機(jī)制、社會(huì)經(jīng)濟(jì)政策等是區(qū)域政策和土地政策法規(guī)影響土地利用結(jié)構(gòu)的直接方式，也是影響土地利用方式的重要因素。20 世紀(jì)80年代實(shí)行包產(chǎn)到戶，提高了勞動(dòng)生產(chǎn)率，隨后實(shí)施農(nóng)業(yè)產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，發(fā)展林果、草畜、設(shè)施農(nóng)業(yè)等高效產(chǎn)業(yè)，大力發(fā)展二三產(chǎn)業(yè)，城市化水平以及社會(huì)經(jīng)濟(jì)發(fā)展水平得到較快提高。同時(shí)，積極開展退耕還林工程和治溝造地土地整治重大工程，從而增加農(nóng)民收入，進(jìn)而推動(dòng)區(qū)域內(nèi)的社會(huì)經(jīng)濟(jì)發(fā)展。

4 小結(jié)

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

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

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

4）建設(shè)用地呈持續(xù)增加趨勢(shì)，由1980 年的116.16 km2增加到2020 年的280.88 km2，增加幅度達(dá)142%。明顯變化時(shí)段主要發(fā)生在2005—2010 年和2015—2020年，分別增加了49.28 km2和74.79 km2。建設(shè)用地增加多集中在延安市城區(qū)、縣城區(qū)和北部少量重點(diǎn)鄉(xiāng)鎮(zhèn)，南部塬區(qū)農(nóng)村建設(shè)用地面積也有增加。延安市北部人口多集中到縣城和延安市城區(qū)，鄉(xiāng)鎮(zhèn)和農(nóng)村居住人口大幅減少。

5）20 世紀(jì)90 年代，受干旱氣候和人口增加雙重影響，延安市出現(xiàn)毀草種地現(xiàn)象。實(shí)行退耕還林政策后至2010年，耕地面積大幅減少，政策對(duì)土地轉(zhuǎn)換起到?jīng)Q定性作用。隨著蘋果、設(shè)施農(nóng)業(yè)發(fā)展，農(nóng)民收入增加，減少了對(duì)土地的依賴，經(jīng)濟(jì)水平和城鎮(zhèn)率提高，農(nóng)民走向城市生活，確保了退耕成果。但在2015年，農(nóng)村人口的減少造成了未利用土地的大幅增加，耕地沒有得到充分利用。氣候?qū)ν恋乩米兓灿忻黠@影響，2000 年以來，延安市處在降水多、溫度高的暖濕時(shí)段，特別是北部降水量，增加了70～80 mm，是林地、高中覆蓋草地面積增加的主要原因。

延安市北部中低覆蓋草地面積較大，生態(tài)比較脆弱，加之氣候還在持續(xù)變暖，氣候不確定性增加，嚴(yán)重制約著延安市的生態(tài)恢復(fù)。因此，為實(shí)現(xiàn)延安市生態(tài)的可持續(xù)發(fā)展和鄉(xiāng)村振興，還需要持續(xù)關(guān)注，開展科學(xué)研究，提供對(duì)策建議。