羅 維，易海杰，李紅舉，賈文濤，馮兆忠

1 中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室，北京 100085 2 國土資源部土地整理中心，北京 100035

洋河流域土地利用時(shí)空變異及其對生態(tài)服務(wù)功能價(jià)值的影響

羅 維1,*，易海杰1，李紅舉2，賈文濤2，馮兆忠1

1 中國科學(xué)院生態(tài)環(huán)境研究中心城市與區(qū)域生態(tài)國家重點(diǎn)實(shí)驗(yàn)室，北京 100085 2 國土資源部土地整理中心，北京 100035

洋河流域位于我國北方典型的農(nóng)牧交錯(cuò)地帶,是首都重要水源涵養(yǎng)地和生態(tài)保護(hù)屏障,也是2022年冬奧會(huì)舉辦地之一。對流域土地利用變化及其生態(tài)服務(wù)功能價(jià)值影響的研究可為首都水源保護(hù)以及冬奧會(huì)場地環(huán)境質(zhì)量的改善提供重要理論依據(jù)和決策支持。利用ArcGIS技術(shù)分析了洋河流域1990—2013年間的土地遙感數(shù)據(jù),估算了這些變化對流域生態(tài)服務(wù)價(jià)值的影響。結(jié)果表明：1990—2013年洋河流域林地和未利用地面積呈增加,而水域、草地則呈減少的趨勢。由單一土地類型利用動(dòng)態(tài)絕對值說明,除林地以外2000—2013年洋河流域其他土地類型受人類或自然環(huán)境變化的影響明顯高于1990—2000年。土地轉(zhuǎn)移主要發(fā)生在流域西部的興和縣、陽高縣、天鎮(zhèn)縣和下游的張家口市、宣化縣和懷來縣。近年來流域土地面積變化幅度和變化范圍均增加。1990—2013年間洋河流域土地利用的年總生態(tài)服務(wù)價(jià)值經(jīng)歷先增加后減小、總體減少的過程,但林地和未利用地的生態(tài)服務(wù)價(jià)值卻一直增加。由于林地面積增加從而導(dǎo)致1990—2013年流域的原材料生產(chǎn)、氣候調(diào)節(jié)、氣體調(diào)節(jié)、生物多樣性保護(hù)和提供美學(xué)景生態(tài)功能占年總價(jià)值的比值增加。林地、草地和耕地是洋河流域生態(tài)服務(wù)的主要貢獻(xiàn)者,也是生態(tài)服務(wù)價(jià)值變化的敏感因子。

農(nóng)牧交錯(cuò)帶；水源地；生態(tài)服務(wù)功能；價(jià)值評(píng)估；土地利用變化

Abstract： The Yanghe Watershed is located in a typical farming pastoral zone of North China. It is an important water-source conservation site and an ecological protection barrier for Beijing. Furthermore, the 2022 Winter Olympic Games will be held in Zhangjiakou City, which is a key region of the Yanghe Watershed. A study of land use changes and their ecosystem service-function values of the watershed can provide important theories for effective decision-making regarding water resource protection and environmental quality improvement in fields of the areas where the games will be held. ArcGIS was used to analyze remote sensing data of land-use changes that occurred within the Yanghe Watershed from 1990 to 2013, and ecosystem service-function values of the lands were estimated. The results showed that forest and unused lands within the watershed increased, whereas grassland and aquatic areas decreased during the study period. Based on absolute values of the changing rate of a single land-use type, apart from forestland, the other land-use types were more greatly affected by human or natural changes in 2000—2013 than in 1990—2000. Land transfers mainly occurred in Xihe, Yanggao, Tianzhen, Xuanhua, and Huailia Counties, and the land change degree and distribution increased in recent years. The annual value of total ecosystem service functions of the lands increased during 1990—2000 and decreased during 2000—2013; however, overall it decreased from 1990 to 2013. Service function values of grassland, aquatic areas, and farmland decreased, whereas that of forest and unused land increased during 1990—2013. Percentages of values of raw material production, climate regulation, air conditioning, biodiversity conservation, and esthetic landscape increased owing to increases of forestland within the watershed during the study period. Forest-, grass-, and farmlands were not only the main contributors but also sensitive factors to ecosystem service function values of the lands within the Yanghe Watershed from 1990 to 2013.

KeyWords: farming pastoral zone; water source site; ecosystem service function; value estimation; land use change

土地利用變化改變土地原有的生態(tài)服務(wù)價(jià)值和功能[1-2],深入研究土地利用時(shí)空格局變化及其對生態(tài)服務(wù)價(jià)值及功能的影響可揭示土地利用對人類活動(dòng)程度的響應(yīng),對于理解人與自然的關(guān)系具有重要的意義。我國經(jīng)濟(jì)的快速發(fā)展和人口的迅速增長[3-4],已使土地利用結(jié)構(gòu)和生態(tài)服務(wù)價(jià)值發(fā)生顯著變化[5-6],倍受國內(nèi)外學(xué)者的高度關(guān)注[7- 10]。我國北方農(nóng)牧交錯(cuò)帶是一條生態(tài)脆弱的過渡帶和生態(tài)安全屏障帶,也是以生態(tài)保護(hù)為主體功能的限制性開發(fā)區(qū),具有獨(dú)特的地理分布格局和土地利用變化過程,已成為研究的熱點(diǎn)[1,11]。流域是跨區(qū)聯(lián)系密切、生態(tài)關(guān)聯(lián)性很高的有機(jī)整體,以流域?yàn)橥暾芯繂卧芯课覈狈睫r(nóng)牧交錯(cuò)帶土地利用及其對生態(tài)服務(wù)功能價(jià)值的影響將更具系統(tǒng)性和挑戰(zhàn)性。近年來我國北方農(nóng)牧交錯(cuò)帶內(nèi)大多數(shù)流域的土地利用和生態(tài)服務(wù)價(jià)值已發(fā)生劇烈變化,利用土地利用變化速度指數(shù)、生態(tài)服務(wù)價(jià)值指數(shù)、生態(tài)價(jià)值貢獻(xiàn)率指數(shù)和敏感性指數(shù)揭示帶內(nèi)典型流域的長時(shí)間序列土地利用變化及其生態(tài)服務(wù)功能價(jià)值的影響對于指導(dǎo)我國北方農(nóng)牧交錯(cuò)帶流域土地資源開發(fā)利用和生態(tài)環(huán)境保護(hù)將具有重要的理論意義和現(xiàn)實(shí)價(jià)值。

洋河流域是我國北方農(nóng)牧交錯(cuò)帶的典型流域,是首都北京重要的水源保護(hù)地和生態(tài)保護(hù)屏障,對京津冀大都市圈的生態(tài)環(huán)境保護(hù)具有很重要的地位和作用[12]。目前已有學(xué)者對洋河流域部分河段的土地利用變化進(jìn)行了初步研究[6, 13- 16],但是對整個(gè)流域土地利用時(shí)空變異及其對生態(tài)服務(wù)功能價(jià)值影響的研究仍相當(dāng)匱乏,這不利于指導(dǎo)整個(gè)區(qū)域的協(xié)調(diào)發(fā)展。另外,隨著洋河流域的張家口市與北京市聯(lián)合舉辦2022年冬奧會(huì)籌備工作的逐步展開,國家將對北京-張家口市的原有土地利用格局做出重大調(diào)整；充分理解洋河流域土地時(shí)空變異及其對生態(tài)服務(wù)價(jià)值和功能的影響可為冬奧會(huì)場地生態(tài)環(huán)境的合理規(guī)劃建設(shè)提供重要參考。本文以遙感解譯數(shù)據(jù)為基礎(chǔ),對整個(gè)洋河流域長時(shí)間序列(1990—2013年)土地利用格局和過程的變化及其生態(tài)系統(tǒng)服務(wù)功能價(jià)值的影響進(jìn)行深入地研究,以期為我國北方農(nóng)牧交錯(cuò)帶流域土地資源的合理利用與管理、為2022年北京-張家口冬奧會(huì)場地的生態(tài)環(huán)境保護(hù)、甚至整個(gè)京津冀大都市圈的可持續(xù)發(fā)展提供重要的理論依據(jù)和決策支持。

1 研究區(qū)域概況

洋河流域地處我國北方西部草原牧區(qū)和東部農(nóng)耕區(qū)之間的農(nóng)牧交錯(cuò)地帶,其經(jīng)緯度113°29′32″E—115°43′40″E,39°59′42″N—41°15′3″N之間,流域西起內(nèi)蒙古自治區(qū)和山西省,由西北向東南進(jìn)流入河北省境內(nèi),在張家口市的懷來縣與桑干河一起匯入永定河,最終注入北京官廳水庫,屬于永定河支流。洋河干流全長約262 km,流域總面積1.62×104km2,沿途經(jīng)過內(nèi)蒙古自治區(qū)烏蘭察布市的興和縣(1個(gè)縣)、山西省大同市陽高縣和天鎮(zhèn)縣(2個(gè)縣)與河北省張家口市(1個(gè)縣)及其管轄8個(gè)縣(圖1)。

圖1 洋河流域地理位置Fig.1 Map of the Yanghe Watershed, China

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

2.1 洋河流域土地利用數(shù)據(jù)來源及分析

1990、2000和2013年洋河流域土地利用數(shù)據(jù)采用中國科學(xué)院資源環(huán)境數(shù)據(jù)中心提供的1∶10萬土地利用數(shù)字化地圖,并結(jié)合相應(yīng)時(shí)段美國陸地衛(wèi)星Landsat TM數(shù)字影像進(jìn)行檢查和校正。參照中國科學(xué)院資源環(huán)境數(shù)據(jù)中心土地分類系統(tǒng),結(jié)合洋河流域土地利用類型特點(diǎn),將研究區(qū)土地類型分為6個(gè)一級(jí)類型：耕地、林地、草地、水域、城鎮(zhèn)用地、工礦用地和未利用地。為保證獲取數(shù)據(jù)的解譯質(zhì)量和一致性,利用GPS在野外對洋河流域土地利用解譯結(jié)果進(jìn)行了選點(diǎn)驗(yàn)證,統(tǒng)計(jì)各土地利用類型的解譯精度均在85%以上。根據(jù)數(shù)字地圖和高程(DEM)數(shù)據(jù),利用ArcGIS 10.2對1990、2000和2013年洋河流域土地?cái)?shù)據(jù)進(jìn)行分析。

2.2 洋河流域不同土地的生態(tài)服務(wù)功能價(jià)值

2.2.1 生態(tài)服務(wù)價(jià)值變化

根據(jù)生態(tài)服務(wù)價(jià)值功能的估算方法[11],其計(jì)算公式如下：

(1)

(2)

(3)

式中,ESVK表示土地類型K的生態(tài)服務(wù)價(jià)值(百元),本研究中6種土地利用類型分別為：耕地K=1、林地K=2、草地K=3、水域K=4、未利用地K=5、城鎮(zhèn)工礦用地K=6。由于城鎮(zhèn)工礦用地的生態(tài)服務(wù)價(jià)值系數(shù)為零[17,18],因而未參與本文生態(tài)服務(wù)價(jià)值的計(jì)算。ESVF代表生態(tài)服務(wù)功能F的服務(wù)價(jià)值(百元),ESV代表流域生態(tài)服務(wù)總價(jià)值(百元),AK表示K類土地的面積(km2),VVCKF表示K類土地的F種生態(tài)服務(wù)功能的價(jià)值系數(shù)(百元 km-2a-1),本研究涉及9種生態(tài)服務(wù)功能,分別是食物生產(chǎn)F=1,原材料生產(chǎn)F=2,氣體調(diào)節(jié)F=3,氣候調(diào)節(jié)F=4,水源涵養(yǎng)F=5,廢物處理F=6,土壤形成與保護(hù)F=7,生物多樣性保護(hù)F=8,提供美學(xué)景觀F=9。利用平均加權(quán)法求得洋河流域生態(tài)服務(wù)價(jià)值的修正系數(shù)：山西為0.46、內(nèi)蒙為0.44、河北為1.02,因而整個(gè)洋河流域修正系數(shù)0.778[17-18],結(jié)合謝高地等對全國土地生態(tài)功能服務(wù)價(jià)值系數(shù)的研究[17-18],得到洋河流域的生態(tài)服務(wù)價(jià)值系數(shù)如表1所示。

表1 洋河流域生態(tài)系統(tǒng)服務(wù)價(jià)值系數(shù)(VVCKF)

2.2.2 生態(tài)貢獻(xiàn)率

分析1990—2000,2000—2013,1990—2013年3個(gè)時(shí)期洋河流域各類土地的生態(tài)貢獻(xiàn)率,可以確定流域生態(tài)服務(wù)價(jià)值的主要貢獻(xiàn)因子和敏感因子。生態(tài)貢獻(xiàn)率計(jì)算公式如下[19]：

(4)

式中,SKT表示K類土地在時(shí)間段T的生態(tài)服務(wù)貢獻(xiàn)率(%); ΔESVKT表示K類土地在時(shí)間段T內(nèi)的生態(tài)服務(wù)價(jià)值變化量(萬元)。

2.2.3 敏感性

用敏感性指數(shù)(CCS)表征生態(tài)服務(wù)價(jià)值隨時(shí)間變化對價(jià)值指數(shù)變化的依賴程度[20- 23],若CCS>1,說明研究區(qū)域內(nèi)的生態(tài)服務(wù)價(jià)值ESV對CCS富有彈性；如果CCS<1,則不富有彈性,說明價(jià)值系數(shù)適合研究區(qū)域。CCS越大,表明準(zhǔn)確性越關(guān)鍵。本文在洋河流域生態(tài)服務(wù)價(jià)值系數(shù)基礎(chǔ)上,上下調(diào)整50%[20],其計(jì)算公式(5)所示[20]：

(5)

式中，ESVI和ESVJ表示價(jià)值系數(shù)調(diào)整前后的總價(jià)值(萬元)；VVCKI和VVCKJ代表值系數(shù)調(diào)整前后K類土地的價(jià)值系數(shù)。

2.2.4 單一土地利用動(dòng)態(tài)度

單一土地利用動(dòng)態(tài)度是刻畫不同土地利用類型在一定時(shí)間段內(nèi)的變化速度和幅度的指標(biāo),反映人類活動(dòng)對單一土地利用類型的影響,其公式為[24]：

(6)

式中，k值為t1到t2時(shí)段內(nèi)某土地利用類型動(dòng)態(tài)度(%),ua、ub為研究初期t1及研究末期t2某類型土地面積(km2),t1為1990和2000年,t2為2000和2013年。

3 結(jié)果與分析

3.1 洋河流域土地利用的時(shí)空演變

1990、2000和2013年洋河流域的土地面積變化如表2所示。洋河流域土地主要類型以耕地、林地和草地為主,分別占總面積的43%,33%和17%以上,與中國北方農(nóng)牧交錯(cuò)帶其他地區(qū)的主要土地類型基本一致[24-25]。1990—2013年洋河流域的耕地、草地和水域面積減少,林地、未利用地面和城鎮(zhèn)工礦用地面積增加,城鎮(zhèn)工礦用地變化幅度相對較大,這可能與2013年洋河流域核心區(qū)域張家口市的均GDP比1990年增加35.252萬元,人口增長23.62萬[26],耕地被大量建設(shè)成住宅和公共設(shè)施有關(guān)[3-4]；該結(jié)論與其他研究相類似。另外,北方的退耕還林還草政策[27-28]是導(dǎo)致洋河流域耕地減少,林地增加的原因之一。

表2 1990—2013年期間洋河流域各土地利用類型面積和占總面積的百分比

土地的轉(zhuǎn)移情況如圖2所示,土地的轉(zhuǎn)移主要發(fā)生在耕地、林地和草地之間。1990—2000年耕地的轉(zhuǎn)出主要發(fā)生在陽高縣和天鎮(zhèn)縣的上游和洋河的張家口河段,在陽高縣和天鎮(zhèn)縣主要是耕地轉(zhuǎn)化成林地和草地,而在洋河的張家口段主要是轉(zhuǎn)化為城鎮(zhèn)礦工用地；林地和草地的轉(zhuǎn)移主要發(fā)生在陽高縣、天鎮(zhèn)縣和尚義縣和興和縣,在陽高縣、天鎮(zhèn)縣林地主要向耕地轉(zhuǎn)移,其次是草地,在興和縣主要向草地轉(zhuǎn)移；草地主要向耕地轉(zhuǎn)移(圖2)。2000—2013年興和縣發(fā)生大面積的耕地轉(zhuǎn)移,在陽高縣、天鎮(zhèn)縣和興和縣耕地主要轉(zhuǎn)化為林地、草地、其次是城鎮(zhèn)工礦用地,在洋河中下游張家口市、宣化縣和懷來縣耕地主要轉(zhuǎn)向城鎮(zhèn)礦工用地；在小清水河上游林地向草地轉(zhuǎn)移。草地在興和縣主要轉(zhuǎn)向耕地,在清水河的張家口市主要轉(zhuǎn)向城鎮(zhèn)礦工用地(圖2)。不同時(shí)期尚義縣、陽高縣和天鎮(zhèn)縣土地主要在林地、草地和耕地之間轉(zhuǎn)移,在張家口市和宣化縣主要由其他用地向城鎮(zhèn)礦工用地轉(zhuǎn)移(圖2)。不同時(shí)期洋河流域的土地變化主要發(fā)生在興和縣、陽高縣、天鎮(zhèn)縣和張家口市,但是2000—2013年比1990—2000年的土地變化范圍更大(圖2)。

圖2 不同時(shí)期洋河流域土地轉(zhuǎn)移和變化分布,Fig.2 Distribution of land transfer and change at different time within the Yanghe Watershed, China1：耕地,2：林地,3：草地,4：水域,5：城鎮(zhèn)礦工用地,6：未利用地

1990—2013年洋河流域各市縣不同土地利用類型面積的變化如表3所示,張北縣和尚義縣的耕地略有增加；其他縣均在減少,減少最多的是陽高縣和張家口市,這是由于旱地的減少所致[14]。萬全縣、興和縣和陽高縣的林地面積增加,增加最大的是陽高縣。興和縣、天鎮(zhèn)縣和懷安縣的草地增加,其中興和縣增加最大。興和縣、陽高縣和張家口市的水域面積減少,其他縣市增加。未利用地的面積變化主要集中在城鎮(zhèn)礦工用地面積較大的懷來縣、懷安縣、張家口市、天鎮(zhèn)縣和陽高縣。每個(gè)縣市的城鎮(zhèn)礦工用地面積都在增加,增加最快的是張家口市,其次是懷來縣和宣化縣。不同時(shí)期耕地、草地、水域、未利用地和城鎮(zhèn)礦工用地的面積變化幅度為：2000—2013年>1990—2000年,而在城鎮(zhèn)面積較大的張家口市、陽高縣、宣化縣和崇禮縣的林地面積變化幅度為：2000—2013>1990—2000,其他縣市則相反,這主要是由于2000年后我國北方實(shí)施的退耕還林還草所致[15]。

3.2 洋河流域土地利用類型變化的速度

1990—2013年間洋河流域單一土地利用類型變化的速度(動(dòng)態(tài)度)k由大到小分別為：城鎮(zhèn)工礦用地>林地>水域>未利用地>耕地>草地(表4)。2000—2013年間流域城鎮(zhèn)工礦用地的k最大,其次是未利用地。除林地以外,2000—2013年其他單一土地利用動(dòng)態(tài)度k的絕對值均高于1990—2000年,說明該時(shí)段洋河流域大多數(shù)土地類型明顯受人類或自然環(huán)境變化的影響。

表3 1990—2013年洋河流域各縣市不同土地類型面積變化

注：負(fù)號(hào)表示面積變化是朝著減少方向進(jìn)行,絕對值大小表示速度的快慢變化

3.3 生態(tài)敏感性

從1990年開始洋河流域林地的敏感性指數(shù)呈逐年增加的趨勢,耕地的敏感性指數(shù)呈減小趨勢(表5),說明林地和耕地生態(tài)服務(wù)功能價(jià)值變化對洋河流域土地生態(tài)服務(wù)總價(jià)值具有放大和縮小效應(yīng)；草地的敏感性指數(shù)先減小后增加(表5),說明草地的價(jià)值系數(shù)變化對區(qū)域內(nèi)生態(tài)服務(wù)總價(jià)值先具縮小效應(yīng)而后具放大效應(yīng)；水域和未利用地的敏感性指數(shù)基本不變(表5),說明二者的價(jià)值系數(shù)變化對區(qū)域內(nèi)生態(tài)服務(wù)總價(jià)值沒有明顯影響。1990—2013年洋河流域各類土地的敏感度指數(shù)均小于1(表5),說明流域生態(tài)系統(tǒng)服務(wù)價(jià)值對本研究采用的價(jià)值系數(shù)缺乏彈性,本文采用的價(jià)值系數(shù)適合洋河流域的情況,計(jì)算的服務(wù)價(jià)值是可信的。

3.4 洋河流域土地利用生態(tài)服務(wù)價(jià)值變化

1990—2000年洋河流域土地利用年度總體生態(tài)服務(wù)價(jià)值增加2557.6萬元, 增長了0.353%,2000—2013年減少5504.1萬元, 減少了0.758%(表6)。除了林地和未利用地,1990—2013年洋河流域其他幾類土地類型的生態(tài)服務(wù)價(jià)值均減少。

表5 1990—2013年間洋河流域各土地類型生態(tài)服務(wù)價(jià)值的敏感性指數(shù)(Ccs)

表6 1990—2013年洋河流域土地利用生態(tài)服務(wù)價(jià)值ESV及其變化

負(fù)號(hào)表示減小

1990—2013年洋河流域不同縣市土地的生態(tài)服務(wù)功能如表7所示,生態(tài)服務(wù)價(jià)值最大的是崇禮縣,其次是興和縣和宣化縣。1990—2013年生態(tài)服務(wù)價(jià)值在尚義縣、天鎮(zhèn)縣和陽高縣的呈先減小后增大的“V”型變化,在萬全縣、興和縣和涿鹿縣則相反。陽高縣1990—2013年的生態(tài)服務(wù)價(jià)值增加了587.2萬元,其他縣均在減少,其中興和縣、懷來縣、宣化縣和張家口市減量較大,主要是由于草地和水域減少所致[16]。

表7 1990—2013年間洋河流域各縣市生態(tài)服務(wù)價(jià)值ESV及其變化(萬元)

負(fù)號(hào)表示減小

3.5 洋河流域各生態(tài)功能服務(wù)價(jià)值的變化

近20年洋河流域各項(xiàng)功能價(jià)值均減少,減少總量達(dá)2946.4萬元/a,其中廢物處理價(jià)值減小了589.0萬元/年,減幅最大,其次減幅由高到低的順序依次分別為土壤形成與保護(hù)功能價(jià)值、水源涵養(yǎng)功能價(jià)值和生物多樣性保護(hù)價(jià)值(表8),這種減小趨勢與河北省其他地區(qū)各項(xiàng)生態(tài)功能價(jià)值增加的變化規(guī)律不同[10],但與內(nèi)蒙古地區(qū)變化相似[9],可能與研究區(qū)的地理位置和生態(tài)條件等有關(guān)。1990—2013年洋河流域的四類一級(jí)生態(tài)服務(wù)功能中調(diào)節(jié)服務(wù)功能和支持服務(wù)功能年總價(jià)值為372547.1—375379.9 萬元和222950.3—244607.7萬元,分別占所有土地年服務(wù)價(jià)值的51%和30%以上。1990—2013年洋河流域的9種二級(jí)功能服務(wù)中,土壤形成與保護(hù)功能價(jià)值最大,其次分別是生物多樣性保護(hù)、水源涵養(yǎng)、氣候調(diào)節(jié)、氣體調(diào)節(jié)、廢物處理、原材料生產(chǎn)、提供美學(xué)景觀和糧食生產(chǎn)的價(jià)值。

表8 1990—2013年間洋河流域各生態(tài)服務(wù)功能價(jià)值ESV及所占百分比

3.6 洋河流域不同土地生態(tài)服務(wù)價(jià)值變化的貢獻(xiàn)率

1990—2000年和2000—2013年流域不同土地生態(tài)服務(wù)價(jià)值變化的貢獻(xiàn)率最大的是林地和耕地,分別為83.171%和43.295%。1990—2013年耕地(55.215%)、林地(12.355%)和草地(20.957%)的生態(tài)貢獻(xiàn)率最大,三者之和超過88%(表9)。1990—2000和2000—2013年兩個(gè)時(shí)期耕地、林地的貢獻(xiàn)率變化分別是36.576%和43.768%,說明1990—2013年耕地和林地是洋河流域生態(tài)服務(wù)價(jià)值變化主要貢獻(xiàn)因子和的敏感因子。

表9 1990—2013年洋河流域不同土地類型生態(tài)服務(wù)價(jià)值變化的貢獻(xiàn)率SKT

4 結(jié)論

(1)1990—2013年洋河流域的主要土地利用類型是林地、耕地和草地,三者占整個(gè)流域面積的93%以上,在該期間內(nèi)耕地、草地和水域的面積減少,其中耕地減少最大,林地、未利用地和城鎮(zhèn)工礦用地的面積增加,城鎮(zhèn)工礦用地增幅度最大；流域增加的林地主要集中于陽高縣、萬全縣和興和縣,其他區(qū)縣市的林地則減小。除林地外,2000—2013年洋河流域其他單一土地利用類型變化的速度高于1990—2000年。

(2)1990—2013年洋河流域土地的轉(zhuǎn)移主要發(fā)生在林地、耕地、草地和城鎮(zhèn)礦工用地之間。流域西部的陽高縣、天鎮(zhèn)縣、興和縣和尚義縣,耕地主要向草地和林地轉(zhuǎn)移,草地主要向耕地轉(zhuǎn)移,林地主要向草地和耕地轉(zhuǎn)移；張家口市和宣化縣和懷來縣的耕地主要轉(zhuǎn)向城鎮(zhèn)礦工用地。2000—2013年流域土地面積變化幅度和變化范圍皆高于1990—2000年,說明近年來人類活動(dòng)對流域土地面積變化的影響在逐漸擴(kuò)大。

(3)1990—2013年洋河流域土地利用的年總生態(tài)服務(wù)價(jià)值減小,耕地、草地和水域的生態(tài)服務(wù)價(jià)值呈減少趨勢,林地和未利用地的生態(tài)服務(wù)價(jià)值呈增加趨勢。林地生態(tài)服務(wù)價(jià)值增加主要集中于陽高縣、萬全縣和興和縣,因?yàn)槿h境內(nèi)的林地面積增加。林地和耕地對洋河流域生態(tài)服務(wù)價(jià)值變化有放大和縮小效應(yīng),是主要貢獻(xiàn)因子和敏感因子。

(4)1990—2013年流域各項(xiàng)生態(tài)服務(wù)功能價(jià)值呈1990—2000年增加,2000—2013年減小,總體減小的趨勢。洋河流域一級(jí)生態(tài)服務(wù)功能價(jià)值由大到小依次是：調(diào)節(jié)服務(wù)功能>支持服務(wù)功能>供給服務(wù)功能>文化服務(wù)功能；二級(jí)生態(tài)服務(wù)功能價(jià)值的順序：土壤形成與保護(hù)功能>生物多樣性保護(hù)>水源涵養(yǎng)>氣候調(diào)節(jié)>氣體調(diào)節(jié)>廢物處理>原材料生產(chǎn)>提供美學(xué)景觀>糧食生產(chǎn)。

[1] 陳海燕, 邵全琴, 安如. 1980s—2005年內(nèi)蒙古地區(qū)生態(tài)系統(tǒng)服務(wù)功能價(jià)值的時(shí)空變化分析. 地球信息科學(xué)學(xué)報(bào), 2014, 16(4): 582- 591.

[2] 元媛, 靳占忠, 劉宏娟, 譚莉梅, 劉浩杰, 劉慧濤, 張廣錄, 劉金銅. 河北省38°N生態(tài)樣帶生態(tài)系統(tǒng)服務(wù)功能時(shí)空變化. 自然資源學(xué)報(bào), 2011, 26(7): 1166- 1179.

[3] 崔曉臨, 白紅英, 王濤, 郭琳, 張加利. 北方農(nóng)牧交錯(cuò)帶靖邊縣土地利用變化及其驅(qū)動(dòng)因素. 水土保持通報(bào), 2013, 33(1): 34- 38, 43- 43.

[4] 毛留喜, 宇振榮, 程序, 王利文, 劉國彬. 北方農(nóng)牧交錯(cuò)帶人口脅迫與耕地利用的相互關(guān)系. 農(nóng)業(yè)工程學(xué)報(bào), 2000, 16(4): 11- 14.

[5] 劉紀(jì)遠(yuǎn), 張?jiān)鱿? 莊大方, 王一謀, 周萬村, 張樹文, 李仁東, 江南, 吳世新. 20世紀(jì)90年代中國土地利用變化時(shí)空特征及其成因分析. 地理研究, 2003, 22(1): 1- 12.

[6] 馬振剛. 基于遙感技術(shù)的洋河流域土地覆被變化特征分析[D]. 長春: 東北師范大學(xué), 2008.

[7] Behrman K D, Juenger T E, Kiniry J R, Keitt T H. Spatial land use trade-offs for maintenance of biodiversity, biofuel, and agriculture. Landscape Ecology, 2015, 30(10): 1987- 1999.

[8] Michelsen O. Assessment of land use impact on biodiversity-Proposal of a new methodology exemplified with forestry operations in Norway. The International Journal of Life Cycle Assessment, 2008, 13: 22- 31.

[9] Riedener E, Rusterholz H P, Baur B. Land-use abandonment owing to irrigation cessation affects the biodiversity of hay meadows in an arid mountain region. Agriculture, Ecosystems & Environment, 2014, 185: 144- 152.

[10] Trisurat Y, Alkemade R, Verburg P H. Projecting land-use change and its consequences for biodiversity in northern Thailand. Environmental Management, 2010, 45(3): 626- 639.

[11] 馬泉來, 高鳳杰, 張志民, 韓文文, 劉洋, 單培明. 東北農(nóng)林交錯(cuò)區(qū)土地利用景觀及生態(tài)服務(wù)價(jià)值變化. 水土保持通報(bào), 2016, 36(1): 265- 273.

[12] 李文潔. 北京市對張家口市開展生態(tài)補(bǔ)償研究[D]. 保定: 河北大學(xué), 2010.

[13] 李黎黎, 馬振剛, 王寶鈞. 基于遙感和GIS技術(shù)的洋河流域土地利用/土地覆被變化特征分析. 安徽農(nóng)業(yè)科學(xué), 2008, 36(16): 6924- 6926.

[14] 郭洪峰, 許月卿, 田媛. 張家口市土地利用結(jié)構(gòu)特征及其動(dòng)態(tài)變化. 水土保持通報(bào), 2013, 33(3): 259- 264.

[15] 田璐, 許月卿, 孫丕苓. 退耕還林還草工程對土地利用/覆被變化及景觀格局的影響——以張家口市為例. 中國農(nóng)業(yè)大學(xué)學(xué)報(bào), 2015, 20(4): 205- 213.

[16] 常上, 李珊珊, 宋豫秦. 張家口市生態(tài)服務(wù)價(jià)值對土地利用變化的響應(yīng). 北京大學(xué)學(xué)報(bào): 自然科學(xué)版, 2015, 51(6): 1149- 1156.

[17] 謝高地, 甄霖, 魯春霞, 肖玉, 陳操. 一個(gè)基于專家知識(shí)的生態(tài)系統(tǒng)服務(wù)價(jià)值化方法. 自然資源學(xué)報(bào), 2008, 23(5): 911- 919.

[18] 謝高地, 肖玉, 甄霖, 魯春霞. 我國糧食生產(chǎn)的生態(tài)服務(wù)價(jià)值研究. 中國生態(tài)農(nóng)業(yè)學(xué)報(bào), 2005, 13(3): 10- 13.

[19] 封建民, 郭玲霞. 陜西省神木縣土地利用格局和生態(tài)服務(wù)價(jià)值變化. 水土保持通報(bào), 2014, 34(6): 293- 298.

[20] 趙晴, 趙旭陽, 劉征. 石家莊市土地利用變化及其生態(tài)服務(wù)功能響應(yīng). 水土保持通報(bào), 2015, 35(3): 242- 249, 255- 255.

[21] Gascoigne W R, Hoag D, Koontz L, Tangen B A, Shaffer T L, Gleason R A. Valuing ecosystem and economic services across land-use scenarios in the Prairie Pothole Region of the Dakotas, USA. Ecological Economics, 2011, 70(10): 1715- 1725.

[22] 彭文甫, 樊淑云, 周介銘, 趙景峰, 楊存建. 基于遙感與GIS的土地利用變化對生態(tài)服務(wù)價(jià)值的影響. 中國農(nóng)學(xué)通報(bào), 2014, 30(5): 195—202.

[23] 劉桂林, 張落成, 張倩. 長三角地區(qū)土地利用時(shí)空變化對生態(tài)系統(tǒng)服務(wù)價(jià)值的影響. 生態(tài)學(xué)報(bào), 2014, 34(12): 3311- 3319.

[24] 吳琳娜, 楊勝天, 劉曉燕, 羅婭, 周旭, 趙海根. 1976年以來北洛河流域土地利用變化對人類活動(dòng)程度的響應(yīng). 地理學(xué)報(bào), 2014, 69(1): 54- 63.

[25] 張偉科, 封志明. 楊艷昭, 孫小舟, 游珍. 北方農(nóng)牧交錯(cuò)帶土地利用/覆被變化分析——以西遼河流域?yàn)槔? 資源科學(xué), 2010, 32(3): 573- 579.

[26] 張家口市統(tǒng)計(jì)局. 張家口經(jīng)濟(jì)年鑒2014(65年特別版). 北京: 中國統(tǒng)計(jì)出版社, 2014.

[27] 周德成, 趙淑清, 朱超. 退耕還林還草工程對中國北方農(nóng)牧交錯(cuò)區(qū)土地利用/覆被變化的影響——以科爾沁左翼后旗為例. 地理科學(xué), 2012, 32(4): 442- 449.

[28] 汪芳甜, 安萍莉, 蔡璐佳, 黃鑫鑫, 郝晉珉. 基于RS與GIS的內(nèi)蒙古武川縣退耕還林生態(tài)成效監(jiān)測. 農(nóng)業(yè)工程學(xué)報(bào), 2015, 31(11): 269- 277.

TemporalandspatialvariationoflanduseandtheirinfluenceonecosystemservicefunctionvaluesintheYangheWatershed,China

LUO Wei1,*，YI Haijie1，LI Hongju2，JIA Wentao2，F(xiàn)ENG Zhaozhong1

1StatekeylaboratoryofUrbanandRegionalEcology,ResearchCenterforEco-EnvironmentalSciences,ChineseAcademyofSciences,Beijing100085,China2LandConsolidationandRehabilitationCenterMinistryofLandandResources,Beijing100035,China

國家自然科學(xué)基金項(xiàng)目(41571479,41271502,C031001)；國土資源部土地整理中心外協(xié)項(xiàng)目：土地整治工程的生態(tài)效益評(píng)價(jià)研究(2017-20)

2016- 05- 18; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

日期:2017- 03- 27

*通訊作者Corresponding author.E-mail: luow@rcees.ac.cn

10.5846/stxb201605180962

羅維，易海杰，李紅舉，賈文濤，馮兆忠.洋河流域土地利用時(shí)空變異及其對生態(tài)服務(wù)功能價(jià)值的影響.生態(tài)學(xué)報(bào),2017,37(16):5342- 5351.

Luo W, Yi H J, Li H J, Jia W T, Feng Z Z.Temporal and spatial variation of land use and their influence on ecosystem service function values in the Yanghe Watershed, China.Acta Ecologica Sinica,2017,37(16):5342- 5351.