侯 蕾，彭文啟，劉培斌，陳渠昌，渠曉東，董 飛

（1.中國水利水電科學研究院 水環(huán)境研究所，北京 100038；2.流域水循環(huán)模擬與調控國家重點實驗室，北京 100038；3.北京市水利規(guī)劃設計研究院，北京 100048）

1 研究背景

土地利用變化［1］是人類活動與自然環(huán)境相互作用最直接的表現(xiàn)形式，是人類社會經(jīng)濟活動行為與自然生態(tài)過程交互和鏈接的紐帶［2］。長期以來，國內外學者對土地利用變化從不同角度進行了研究。丑慶川等［3］整理闡述了土地利用變化的分析方法；在流域尺度，近年來定量化分析土地利用變化的水文效應研究為流域水資源管理和規(guī)劃提供了科學依據(jù)［4-7］；于興修等［8］闡述了土地利用變化的環(huán)境效應的研究進展與動向，并重點對土地利用變化的土壤和水環(huán)境效應研究趨勢進行了分析；隨著城市化進程的發(fā)展，土地利用變化對環(huán)境及生態(tài)系統(tǒng)的影響越來越顯著，目前國內外關于土地利用變化的生態(tài)環(huán)境效應的研究主要集中于土地利用變化對大氣成分和區(qū)域氣候變化［9］、生物多樣性［10］、地下水［11-12］、地表徑流［13］和水環(huán)境［14-15］等生態(tài)環(huán)境單一因素的研究；近年來國內的學者針對土地利用變化對區(qū)域整體生態(tài)環(huán)境的效應進行了研究，研究區(qū)域從生態(tài)環(huán)境脆弱區(qū)或經(jīng)濟落后地區(qū)發(fā)展到經(jīng)濟發(fā)達地區(qū)［16-18］。

近幾十年來，永定河流域的生態(tài)環(huán)境問題越來越顯著，尤其20世紀80年代后，隨著經(jīng)濟社會的不斷發(fā)展，該區(qū)域人口增加，取用水量增加，地下水資源量減少，下墊面條件（土地利用/覆被）劇烈變化，永定河流域部分河段出現(xiàn)了斷流，河道徑流量嚴重下降，生態(tài)系統(tǒng)退化。目前，對于永定河流域的土地利用變化的研究已經(jīng)開展，張玉虎等［19］、王海平等［20］采用傳統(tǒng)的分析方法分別對永定河典型山區(qū)-門頭溝區(qū)、永定河（北京段）的土地利用變化的格局以及過程進行了定量研究；張琳等［21］基于1980年、2000年的遙感影像分析了海河流域各支流的土地利用/覆被變化，揭示了海河流域土地利用/覆被變化的規(guī)律。針對以上研究，目前多數(shù)研究集中在分析永定河典型區(qū)域的土地利用變化，而對于永定河上游流域土地利用變化的評價及其生態(tài)環(huán)境的效應研究較少。本研究基于土地利用數(shù)據(jù)，應用ArcGIS平臺，對永定河上游流域土地利用時空動態(tài)變化及其生態(tài)環(huán)境效應進行了分析，旨在進一步認識永定河上游流域下墊面條件的變化情況，為更合理的改善生態(tài)環(huán)境提供科學支撐。

2 研究區(qū)概況

永定河流域是海河流域的七大水系之一，該區(qū)處于干旱和濕潤氣候的過渡地帶，屬中緯度大陸性季風氣候，四季分明，寒暑懸殊。流域多年平均降水不足450 mm，且徑流量年內分配很不均勻，80%左右集中在汛期。該研究區(qū)（圖1）為永定河上游流域（官廳水庫以上），地跨山西、內蒙古、河北等省、自治區(qū)、直轄市。地理坐標為111°58′—116°22′E ，38°50′—41°16′N ，研究區(qū)面積約為 4.3 萬km2，山區(qū)面積占80%。地勢自西北向東南傾斜，西北為山區(qū)，海拔高程一般為500～1000 m，山區(qū)內山巒起伏，地形復雜，河流坡陡流急，且水土流失嚴重。

圖1 永定河上游流域位置

3 數(shù)據(jù)來源及研究方法

3.1 數(shù)據(jù)來源及處理永定河上游流域20世紀70年代末期（1980年）、80年代末期（1990年）、90年代末期（2000年）以及2010年的4期土地利用/覆被數(shù)據(jù)由中國科學院地理科學與資源研究所中國科學院資源環(huán)境科學數(shù)據(jù)中心提供，分辨率為30m×30m。該數(shù)據(jù)是根據(jù)劉紀遠［22］提出的中國科學院土地資源分類系統(tǒng)，利用永定河上游流域20世紀70年代末期（1980年）、80年代末期（1990年）、90年代末期（2000年）以及2010年的Landsat-MSS/TM/ETM遙感圖像，通過人機交互解譯生成，并利用野外調查數(shù)據(jù)進行驗證和糾正，精度可靠。該分類體系主要有兩級（表1）：將研究區(qū)的土地利用類型進行處理，主要對耕地、林地、草地、水域、城鄉(xiāng)、工礦、居民用地和未利用土地6個類型（一級分類系統(tǒng)）進行分析；土地利用二級分類系統(tǒng)有較高分辨率，且體現(xiàn)了生態(tài)差異性，根據(jù)研究區(qū)域特征，參考國內相關研究，對二級分類體系下的土地利用類型所具有的生態(tài)環(huán)境質量進行賦值［23］。

表1 土地利用分類系統(tǒng)及其生態(tài)環(huán)境質量指數(shù)賦值

3.2 研究方法土地資源數(shù)量變化可以用土地利用動態(tài)度來表示。

（1）單一土地利用類型動態(tài)度。單一土地利用類型動態(tài)度［3，24］表達的是某研究區(qū)一定時間范圍內某種土地利用類型的數(shù)量變化情況，其表達式為：

式中：K為研究時段內某一土地利用類型動態(tài)度；Ua、Ub分別為研究期初及研究期末某一種土地利用類型的數(shù)量；T為研究時長，當T為時段設定為年時，K的值就是該研究某種土地利用類型年變化率。

（2）土地利用綜合動態(tài)度。研究區(qū)的綜合土地利用動態(tài)度［24］可表示為：

式中：LUi為監(jiān)測起始時間第i類土地利用類型面積；ΔLUi-j為監(jiān)測時段內第i類土地利用類型轉為非i類土地利用類型面積的絕對值；T為監(jiān)測時段長度。當T的時段設定為年時，LC的值就是該研究區(qū)土地利用年變化率。

土地利用動態(tài)度定量地描述了土地利用的變化速度，對預算未來土地利用變化趨勢有積極作用。

（3）土地利用類型的轉移矩陣。土地利用轉移矩陣是表征土地利用類型之間在一定時段內相互轉移的一種有效的工具，從而總體上把握土地利用的轉移態(tài)勢。在ArcGIS10.2平臺下，對1980年、1990年、2000年以及2010年4期的土地利用數(shù)據(jù)兩兩疊加分析，對所得的相關數(shù)據(jù)進行統(tǒng)計分析，分別獲得1980—1990年、1990—2000年以及2000—2010年以分析流域土地利用時空格局動態(tài)。

單一土地利用轉出率［20］主要反映某一土地利用類型在某一時期內轉化為其他地類的幅度，計算公式為：

式中：Ti為地類i在t0～tk時期內的土地利用轉出率；Tij為地類在t0～tk時期內地類i轉化為j的面積；Lt0為地類i在t0時刻的面積；n為研究區(qū)土地利用類型數(shù)量；單一土地利用轉入率主要反映某一土地利用類型在某一時期內由其他地類轉化而來的幅度，Mi為地類在t0到tk時期內的土地利用轉入率；Mij為地類在t0到tk時期內由地類j轉化為地類i的面積；Lik為地類i在tk時刻的面積。

（4）區(qū)域生態(tài)環(huán)境質量指數(shù)及土地利用變化類型生態(tài)貢獻率［23］。綜合考慮區(qū)域各土地利用所具有的生態(tài)質量及面積比例，定量表征某一區(qū)域內生態(tài)環(huán)境質量的總體狀況，其表達式為：

式中：EVt為區(qū)域生態(tài)環(huán)境質量指數(shù)；LUi和Ci為該區(qū)域內t時期第i種土地利用類型的面積和生態(tài)環(huán)境質量指數(shù)（見表1）；TA為該區(qū)域總面積；N為區(qū)域內土地利用類型數(shù)量。

土地利用變化類型生態(tài)貢獻率指某一種土地利用類型變化所導致的區(qū)域生態(tài)質量的改變，其表達式為：

式中：LEI為土地利用變化類型生態(tài)貢獻率；LEt、LEt+1分別為某種土地利用變化類型所反映的變化初期和末期土地利用類型所具有生態(tài)質量指數(shù)；LA為該變化類型的面積；TA為區(qū)域總面積。

4 結果分析

4.1 永定河上游流域土地利用變化規(guī)律

4.1.1 永定河上游流域土地利用面積變化 永定河上游流域土地利用類型在數(shù)量組成、時間變化及空間構成上存在一定的差異，4期土地利用類型的面積所占比例由大到小排列，為耕地＞草地＞林地＞城鄉(xiāng)、工礦、居民用地＞水域＞未利用土地，在數(shù)量上（圖2），其中以耕地為主，其面積所占比例達到45%左右；草地與林地次之，分別所占比例為29%、20%左右；水域與城鄉(xiāng)、工礦、居民用地面積較少，分別所占比例為1.5%、2.6%左右；植被覆蓋度較高，其中沙地、鹽堿地、沼澤地、裸土地、裸巖石質地等未利用土地的面積極少，所占比例為0.5%左右。

從年變化（如圖3）上分析，30年（1980—2010年）間，永定河上游流域土地利用類型整體上的變化較為緩慢：林地以及城鄉(xiāng)、工礦、居民用地增加，其中林地增加面積約為1 066 km2，城鄉(xiāng)、工礦、居民用地增加面積約為1 254 km2；耕地、草地、水域、未利用土地均減少，耕地減少約為1 261 km2，草地減少約779 km2，水域減少約158 km2，未利用土地減少約121 km2。增加與減少的土地利用面積基本保持一致。

1980—1990年，研究內草地以及城鄉(xiāng)、工礦、居民用地、未利用土地呈增加趨勢，耕地、林地、水域呈減少趨勢；1990—2000年，研究區(qū)內水域以及城鄉(xiāng)、工礦、居民用地呈增加趨勢，耕地、林地、草地以及未利用土地呈減少趨勢；2000—2010年，研究區(qū)內林地與城鄉(xiāng)、工礦、居民用地呈增加趨勢，耕地、草地、水域、未利用土地呈減少趨勢。

4.1.2 永定河上游流域土地利用變化強度 永定河上游流域各土地利用類型各時段動態(tài)度（如表2）起伏不大，變化率的絕對值均不超過10%。從單一土地利用動態(tài)角度，1980—1990年間，草地、未利用土地以及城鄉(xiāng)、工礦、居民用地的動態(tài)度大于0，年增加率分別為0.12%、0.32%、1.06%，耕地、林地、水域的動態(tài)度小于0，年減少率分別為0.04%、0.22%、0.42%。單一土地利用動態(tài)度的絕對值反映了土地利用類型面積變化速度的大小，因此，1980—1990年間，流域內城鄉(xiāng)、工礦、居民用地的變化速度最大，其次為水域、未利用土地，耕地的變化速度最??；1990—2000年間，流域內除了水域面積以及城鄉(xiāng)、工礦、居民用地增加外，其余土地利用類型面積均減少，城鄉(xiāng)、工礦、居民用地變化速度最大，其次為未利用土地，林地、草地的變化速度最小；2000—2010年間，城鄉(xiāng)、工礦、居民用地持續(xù)增加，變化速度最大，另外林地也呈增長趨勢，未利用土地、耕地、草地、水域面積縮減，其中未利用土地減少速度較大。

圖2 永定河上游流域各土地利用類型面積分布圖

圖3 永定河上游流域1980—2010年土地利用面積變化圖

表2 永定河上游流域土地利用動態(tài)變化

整體上，1980—2010年間，隨著城市化的發(fā)展，人類活動對自然干擾日益嚴重，流域內土地利用變化速度最大的是城鄉(xiāng)、工礦、居民用地，呈增加趨勢，其次為未利用土地，呈減少趨勢。

從綜合土地利用動態(tài)度（表2）的角度，1980—1990、1990—2000、2000—2010年的動態(tài)度分別為0.06%、0.19%、0.6%。相對于2000—2010年，1980—1990、1990—2000年綜合土地利用動態(tài)度較小，說明該兩個時段各土地利用類型變化較??；2000—2010年間綜合土地利用動態(tài)度較大，每年每100 km2，土地利用面積變化為0.6 km2，其變化幅度相對較大，土地利用類型變化較大。

本研究對永定河上游流域1980—1990、1990—2000、2000—2010年3個時段土地利用轉換特征進行了分析。

1980—1990年間，永定河上游流域內耕地主要轉化為城鄉(xiāng)、工礦、居民用地，轉化面積為100.61 km2；林地主要轉化為草地，轉化面積為225.29 km2；草地、水域、未利用土地以及城鄉(xiāng)、工礦、居民用地的轉化均較少。

1990—2000年間，各土地利用類型間的轉化面積較少，土地利用變化甚微。

2000—2010年間，耕地轉化為林地，城鄉(xiāng)、工礦、居民用地和草地，轉化面積依次為887.03、829.96和522.82 km2；林地轉化為草地，耕地和城鄉(xiāng)、工礦、居民用地，轉化面積依次為274.79、157.84和91.06 km2；草地轉化為林地，耕地和城鄉(xiāng)、工礦、居民用地，轉化面積依次為854.59、715.06和216.23 km2；水域主要轉化為耕地，轉化面積為132.53 km2；城鄉(xiāng)、工礦、居民用地主要轉化為耕地，轉化面積為124.11 km2；未利用土地主要轉化為耕地、草地，轉化面積為93.12、37.15 km2（見表3）。

綜上，不同時段，隨著自然、經(jīng)濟條件以及人為影響的變化，土地利用格局也伴隨著變化。1980—1990年、1990—2000年這兩個時段永定河上游流域每種土地利用類型的轉出率或轉入率基本上均小于10%（見表4），土地利用類型總體上變化較小，相對而言，2000—2010年，永定河上游流域每種土地利用類型的轉出率或轉入率均大于10%（見表4），土地利用變化強度較大。

表3 2000—2010年研究區(qū)土地利用轉移矩陣 （單位：km2）

表4 研究區(qū)1980—2010年土地利用變化情況

4.2 永定河上游流域生態(tài)環(huán)境變化與驅動因素分析2000—2010年土地利用變化強度較大，文中基于2000年和2010年土地利用變化分析永定河上游流域的生態(tài)環(huán)境變化。根據(jù)公式（3），計算得到研究區(qū)2000年和2010年區(qū)域生態(tài)環(huán)境質量指數(shù)，結果分別為0.4230、0.4290，永定河上游流域生態(tài)環(huán)境質量指數(shù)年均上升0.07%。從計算結果可知，2000—2010年間研究區(qū)總體生態(tài)環(huán)境質量略有好轉，呈上升趨勢，總體上維持著區(qū)域生態(tài)環(huán)境的動態(tài)穩(wěn)定。

2000—2010年致使研究區(qū)域生態(tài)環(huán)境質量改善的主要土地利用變化的驅動因子主要是退耕還林、還草以及林、草地內部的轉化（詳見表5），其中退耕還林、還草的轉化占生態(tài)貢獻率的54.19%，因永定河上游流域是水土流失的重點區(qū)域，退耕還林、還草均勻散落在永定河上游流域，主要分布在河北省張家口、尚義，山西省大同、左云、朔州等地區(qū)（詳見圖4）；林、草地內部的轉化占生態(tài)貢獻率的37.82%。

2000—2010年期間導致研究區(qū)域生態(tài)環(huán)境質量惡化的主要土地利用變化的驅動因子是高覆蓋草地向旱地的轉化、城鄉(xiāng)居民建設用地規(guī)模擴大、毀林造田、圍水造田和耕地內部轉化（詳見表5），其中城鄉(xiāng)居民建設用地擴大主要集中在張家口市市轄區(qū)、大同市市轄區(qū)、朔州市市轄區(qū)及各個縣轄區(qū)（詳見圖4），毀林造田主要集中在各支流的上游區(qū)域，圍水造田沿桑干河、洋河及各大支流呈帶狀分布；其中高覆蓋草地向旱地的轉化占生態(tài)貢獻率的34.23%，建設占地占生態(tài)貢獻率的26.03%，城鄉(xiāng)居民建設用地規(guī)模擴大主要表現(xiàn)是耕地、林地、草地轉化為城鄉(xiāng)居民建設用地。

5 結論

基于1980年，1990年，2000年和2010年4期土地利用數(shù)據(jù)，系統(tǒng)的分析了永定河上游流域土地利用類型的格局變化。

從土地利用的數(shù)量分析，30年間永定河上游流域土地利用的面積所占比例由大到小排列為耕地＞草地＞林地＞城鄉(xiāng)、工礦、居民用地＞水域＞未利用土地。

表5 導致研究區(qū)域生態(tài)環(huán)境改善和惡化的主要土地利用類型變化及貢獻率

圖4 導致研究區(qū)域生態(tài)環(huán)境改善和惡化的主要影響因子空間分布

從土地利用動態(tài)度以及土地利用類型轉移矩陣分析，1980—1990、1990—2000年，土地利用綜合動態(tài)度分別為0.06%、0.19%，每種土地利用類型的轉出率或轉入率基本上均小于10%，土地利用變化相對較小，耕地主要轉化為林地、草地（退耕還林、還草），同時為治理水土流失修建梯田，部分耕地侵占了林地和草地，城鄉(xiāng)建設用地有增加趨勢；2000—2010年，該時段土地利用綜合動態(tài)度較大，達到了6%，每種土地利用類型的轉出率或轉入率均大于10%，由于經(jīng)濟社會以及城鎮(zhèn)化進程的發(fā)展，其中城鄉(xiāng)、工礦、居民用地不斷擴張，建設占地不斷增加，耕地向林、草地轉化明顯，未利用土地呈減少趨勢，土地利用類型面積變化較大。

從生態(tài)環(huán)境變化與驅動因素分析，2000—2010年間永定河上游流域生態(tài)環(huán)境總體上維持著區(qū)域生態(tài)環(huán)境的動態(tài)穩(wěn)定，生態(tài)環(huán)境改善和惡化的兩種趨勢并存，但生態(tài)環(huán)境質量指數(shù)從0.4230上升到0.4290，整體上生態(tài)環(huán)境質量略有好轉，呈上升趨勢。研究結果表明導致生態(tài)環(huán)境改善的主要因素為退耕還林、還草以及林、草地內部的轉化，導致生態(tài)環(huán)境改善的主要因素為高覆蓋草地向旱地的轉化、城鄉(xiāng)居民建設用地規(guī)模擴大、毀林造田、圍水造田和耕地內部轉化等。因此，在土地利用過程中，要限制或取消導致生態(tài)環(huán)境惡化的人類活動，加強生態(tài)環(huán)境改善力度，充分利用生態(tài)系統(tǒng)的自我修復功能，達到恢復或改善生態(tài)環(huán)境的目的。

[1 ]LEEMANS R.Changes in Land use and land cover：A global perspective［J］.Trends in Ecology&Evolution，1995，10（6）：258-259.

[2]劉紀遠，匡文慧，張增祥，等.20世紀80年代末以來中國土地利用變化的基本特征與空間格局［J］.地理學報，2014（1）：3-14.

[3]丑慶川，胡濤，何詩雨，等 .土地利用/覆被變化分析方法研究［J］.安徽農(nóng)業(yè)科學，2015，43（10）：301-304.

[4]王艷君，呂宏軍，施雅風，等.城市化流域的土地利用變化對水文過程的影響——以秦淮河流域為例［J］.自然資源學報，2009，24（1）：30-35.

[5]歐春平，夏軍，王中根，等.土地利用/覆被變化對SWAT模型水循環(huán)模擬結果的影響研究——以海河流域為例［J］.水力發(fā)電學報，2009（4）：124-129.

[6]董國強，楊志勇，于贏東.下墊面變化對流域產(chǎn)匯流影響研究進展［J］.南水北調與水利科技，2013，11（3）：111-117，126.

[7]WOLDESENBET T A，ELAGIB N A，RIBBE L，et al.Hydrological responses to land use/cover changes in the source region of the Upper Blue Nile Basin，Ethiopia［J］ .Science of the Total Environment，2016，575：724-741.

[8]于興修，楊桂山，王瑤 .土地利用/覆被變化的環(huán)境效應研究進展與動向［J］.地理科學，2004（5）：627-633.

[9]陳溪，王子彥，匡文慧.土地利用對氣候變化影響研究進展與圖譜分析［J］.地理科學進展，2011，30（7）：930-937.

［10]ALBERTI M.The effects of urban patterns on ecosystem function［J］.International Regional Science Review，2005，28（2）：168-192.

［11]張喜鳳，張?zhí)m慧，顧娟，等.敦煌地下水時空變異特征及其對土地利用/覆被變化的響應分析［J］.蘭州大學學報：自然科學版，2014，50（3）：311-317.

［12]秦俊桃.武威盆地不同土地利用類型下的地下水位時空動態(tài)分析［J］.中國水利水電科學研究院學報，2014，12（3）：287-290.

［13]袁宇志，張正棟，蒙金華，等.基于SWAT模型的流溪河流域土地利用與氣候變化對徑流的影響［J］.應用生態(tài)學報，2015，26（4）：989-998.

［14]BAJOCCO S，DE ANGELIS A，PERINI L，et al.The impact of land use/land cover changes on land degradation dynamics：a Mediterranean case study［J］.Environmental management，2012，49（5）：980-989.

［15]K?NDLER M，BLECHINGER K，SEIDLER C，et al.Impact of land use on water quality in the upper Nisa catch?ment in the Czech Republic and in Germany［J］.Science of the Total Environment，2017，586：1316-1325.

［16]張芳怡，邢元志，濮勵杰，等.蘇州市土地利用變化的生態(tài)環(huán)境效應研究［J］.水土保持研究，2009，16（5）：98-103.

［17]張揚，劉艷芳，顧漸萍，等.武漢市土地利用覆被變化與生態(tài)環(huán)境效應研究［J］.地理科學，2011，31（10）：1280-1285.

［18]羅慧娟，黃詩峰，臧文斌，等.5.12汶川地震重災區(qū)土地利用變化遙感監(jiān)測與分析［J］.中國水利水電科學研究院學報，2011，9（3）：209-214.

［19]張玉虎，賈海峰，于長青.永定河流域典型區(qū)土地利用/覆被變化［J］.山地學報，2009，27（5）：564-572.

［20]王海平，馮仲科，孫赫，等.北京市永定河流域土地利用/覆蓋變化研究［J］.林業(yè)調查規(guī)劃，2011，36（3）：15-20.

［21]張琳.海河流域下墊面變化情況及趨勢分析［D］.天津：天津大學，2009.

［22]劉紀遠.中國資源環(huán)境遙感宏觀調查與動態(tài)研究［M］.北京：中國科學技術出版社，1996.

［23]李曉文，方創(chuàng)琳，黃金川，等.西北干旱區(qū)城市土地利用變化及其區(qū)域生態(tài)環(huán)境效應-以甘肅河西地區(qū)為例［J］.第四紀研究，2003，23（3）：280-289.

［24]王秀蘭，包玉海.土地利用動態(tài)變化研究方法探討［J］.地理科學進展，1999，18（1）：83-89.