吉 瑋，王光寶，李 麗，查國(guó)震，張紅月

（1. 江蘇速度遙感大數(shù)據(jù)研究院有限公司，江蘇 南京 210042； 2. 閩江學(xué)院 海洋學(xué)院，福建 福州 350108；3. 江蘇海洋大學(xué) 海洋技術(shù)與測(cè)繪學(xué)院，江蘇 連云港 222005）

生態(tài)系統(tǒng)服務(wù)是指人類(lèi)從生態(tài)系統(tǒng)中所獲取的所有惠益[1]。開(kāi)展生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估，可以為優(yōu)化生態(tài)系統(tǒng)結(jié)構(gòu)提供決策信息，對(duì)促進(jìn)區(qū)域經(jīng)濟(jì)和生態(tài)系統(tǒng)的可持續(xù)發(fā)展具有重要作用。而土地利用變化是制約人類(lèi)社會(huì)可持續(xù)發(fā)展的主要因素之一，土地利用變化能通過(guò)改變生態(tài)系統(tǒng)的功能和結(jié)構(gòu)，對(duì)生態(tài)系統(tǒng)維持及其服務(wù)功能起到?jīng)Q定性作用[2-3]。隨著社會(huì)的不斷進(jìn)步，公眾對(duì)生態(tài)環(huán)境的期望不斷提高，以及“綠水青山就是金山銀山”等科學(xué)論斷的提出，基于土地利用的生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估研究成為國(guó)際生態(tài)學(xué)和生態(tài)經(jīng)濟(jì)學(xué)的研究熱點(diǎn)和前沿領(lǐng)域。

目前，基于土地利用數(shù)據(jù)開(kāi)展生態(tài)系統(tǒng)服務(wù)價(jià)值的評(píng)估主要有3 種：第一種是靜態(tài)評(píng)估，如Gren 等[4]對(duì)歐洲多瑙河流域進(jìn)行了生態(tài)服務(wù)價(jià)值的評(píng)估；蔡中華等[5]使用2010 年中國(guó)生態(tài)系統(tǒng)的統(tǒng)計(jì)數(shù)據(jù)，對(duì)我國(guó)2010 年生態(tài)系統(tǒng)服務(wù)價(jià)值進(jìn)行了再計(jì)算；丁小迪等[6]利用Google Earth 和資源三號(hào)衛(wèi)星影像數(shù)據(jù)針對(duì)山東省濱海濕地的不同生態(tài)服務(wù)功能類(lèi)型，采取市場(chǎng)價(jià)值法、生產(chǎn)價(jià)值法、生態(tài)價(jià)值法等進(jìn)行了初步價(jià)值評(píng)估。第二種是動(dòng)態(tài)評(píng)估，如Mendoza-González 等[7]選取墨西哥灣中部3 個(gè)研究區(qū)，分析了1995-2006 年研究區(qū)土地利用變化及其生態(tài)系統(tǒng)服務(wù)價(jià)值的變化；岳書(shū)平等[8]選取東北樣帶典型區(qū)，分析了1976-2000 年土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的影響；湯峰等[9]運(yùn)用模型法定量評(píng)價(jià)了2009-2016 年高陽(yáng)縣土地利用變化和生態(tài)系統(tǒng)服務(wù)價(jià)值變化及土地碳排放風(fēng)險(xiǎn)；帕茹克·吾斯曼江等[10]分析了1993-2017 年尉犁縣土地利用/覆被變化導(dǎo)致的生態(tài)系統(tǒng)服務(wù)價(jià)值變化。第三種是情景分析，如Nelson 等[11]利用InVEST 模型，分析了1990-2050 年美國(guó)俄勒岡州威拉米特盆地土地利用變化對(duì)生態(tài)系統(tǒng)服務(wù)、生物多樣性保護(hù)和商品生產(chǎn)水平變化的影響。竇攀烽等[12]利用InVEST 模型的情景分析的方法設(shè)計(jì)了6 種情景，分析了2000-2015 年氣候變化和土地利用/覆被變化對(duì)區(qū)域產(chǎn)水量的影響。

黑河市作為黑龍江省的北部門(mén)戶(hù)，雖然工業(yè)基礎(chǔ)薄弱，但是具有生態(tài)環(huán)境優(yōu)美、人文歷史獨(dú)特的相對(duì)優(yōu)勢(shì)。而擁有“北方游獵第一鄉(xiāng)”美譽(yù)的愛(ài)輝區(qū)，近年來(lái)依托“兩山理論”構(gòu)建全域旅游空間體系，努力將生態(tài)優(yōu)勢(shì)轉(zhuǎn)化為經(jīng)濟(jì)優(yōu)勢(shì)[13]，因此對(duì)其開(kāi)展生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估具有重要的指導(dǎo)意義。本文基于2000 年、2007 年、2013 年三期的Landsat 遙感影像，利用監(jiān)督分類(lèi)方法，對(duì)黑河市愛(ài)輝區(qū)進(jìn)行土地利用分類(lèi)；結(jié)合土地利用分類(lèi)結(jié)果，采用中國(guó)陸地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法，參考《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》[14]對(duì)愛(ài)輝區(qū)開(kāi)展了生態(tài)系統(tǒng)服務(wù)價(jià)值動(dòng)態(tài)評(píng)估，以期通過(guò)空間信息技術(shù)對(duì)愛(ài)輝區(qū)區(qū)域規(guī)劃、生態(tài)系統(tǒng)結(jié)構(gòu)優(yōu)化及其可持續(xù)發(fā)展提供數(shù)據(jù)參考和科學(xué)依據(jù)。

1 研究區(qū)概況與數(shù)據(jù)來(lái)源

愛(ài)輝區(qū)隸屬于“胡煥庸線(xiàn)”的北起點(diǎn)黑龍江省黑河市，全區(qū)總面積14 446 km2，地理位置在49°24′38″N ～50°59′11″N，125°49′32″E ～127°36′10″E 之間（圖1），東鄰黑龍江與俄羅斯阿穆?tīng)栔菹嗤?，西接嫩江縣，南聯(lián)孫吳縣，北鄰呼瑪縣[15]。愛(ài)輝區(qū)地處寒溫帶大陸性氣候，春季高溫多風(fēng)，夏季炎熱多雨，秋季降溫極快，冬季寒冷干燥。全年降雨量大約在500 ～600 mm，日照時(shí)數(shù)在2 600 h，無(wú)霜期大約為105 d，年平均氣溫為-0.5 ℃，全年溫差極大，最低氣溫只有-40℃，最高氣溫可達(dá)38℃。愛(ài)輝區(qū)是黑河市政府所在地，同時(shí)也是全國(guó)首批沿邊開(kāi)放城市、國(guó)家級(jí)生態(tài)示范區(qū)和大小興安嶺生態(tài)功能限制開(kāi)發(fā)區(qū)。

圖1 研究區(qū)

本研究中使用的數(shù)據(jù)包括遙感影像數(shù)據(jù)和統(tǒng)計(jì)數(shù)據(jù)兩類(lèi)。其中遙感影像數(shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云平臺(tái)（http://www.gscloud.cn/），由于愛(ài)輝區(qū)可用于土地利用分類(lèi)的高質(zhì)量遙感影像（云量低于1%）有限，本研究從中選取了2000-06-13、2007-09-05 兩期的Landsat TM5 數(shù)據(jù)和2013-08-30 的Landsat8 OLI 數(shù)據(jù)。研究中的統(tǒng)計(jì)數(shù)據(jù)主要來(lái)源于文獻(xiàn)資料。

2 研究方法

2.1 土地利用分類(lèi)與驗(yàn)證

土地利用分類(lèi)方法主要是通過(guò)ENVI 軟件的監(jiān)督分類(lèi)結(jié)合人工目視解譯來(lái)進(jìn)行分類(lèi)和結(jié)果修正。本研究利用ENVI 軟件分別對(duì)2000 年、2007 年和2013 年的Landsat 影像數(shù)據(jù)進(jìn)行圖像預(yù)處理，包括輻射定標(biāo)、大氣校正、校正后地類(lèi)光譜檢查、數(shù)據(jù)鑲嵌和數(shù)據(jù)裁剪。由于Landsat 影像數(shù)據(jù)的幾何精度已經(jīng)能夠滿(mǎn)足本研究的要求，所以不再對(duì)其進(jìn)行幾何校正。對(duì)預(yù)處理后的遙感影像建立分類(lèi)特征庫(kù)，利用監(jiān)督分類(lèi)方法中分類(lèi)效果較好的最大似然分類(lèi)法進(jìn)行研究區(qū)的土地利用分類(lèi)，并對(duì)分類(lèi)結(jié)果進(jìn)行聚類(lèi)處理，結(jié)合人工目視解譯對(duì)分類(lèi)結(jié)果進(jìn)行修正。依據(jù)原國(guó)土資源部土地利用現(xiàn)狀分類(lèi)標(biāo)準(zhǔn)，根據(jù)愛(ài)輝區(qū)土地利用的特征，將土地利用類(lèi)型分為林地、農(nóng)田、草地、水體和其他地類(lèi)共5 種土地利用類(lèi)型。采用總體分類(lèi)精度開(kāi)展分類(lèi)結(jié)果的精度驗(yàn)證，可以通過(guò)ENVI 軟件中的混淆矩陣工具計(jì)算得出。

2.2 生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估方法

本研究采用中國(guó)陸地生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)量方法，參考《森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范》，基于愛(ài)輝區(qū)土地利用分類(lèi)結(jié)果，分別對(duì)林地、農(nóng)田、草地及水體等對(duì)生態(tài)系統(tǒng)服務(wù)貢獻(xiàn)較大的地類(lèi)開(kāi)展了價(jià)值測(cè)算。其他地類(lèi)中包含的建設(shè)用地、裸地等對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值的貢獻(xiàn)較少，本研究對(duì)此不作考慮。

2.2.1 林地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估

1）涵養(yǎng)水源。涵養(yǎng)水源主要包括對(duì)于降雨等水源進(jìn)行保留，防止水分流失，對(duì)其進(jìn)行保存，主要功能是增加可利用的水資源并優(yōu)化水源的質(zhì)量，包含調(diào)節(jié)水量和凈化水質(zhì)兩個(gè)部分。調(diào)節(jié)水量：

式中，UT代表林地每年調(diào)節(jié)水量?jī)r(jià)值（單位：元/a）；Cs代表水庫(kù)建設(shè)單位庫(kù)容投資，依據(jù)牛安逸等[16]的研究，取值為6.11 元/m3；A 代表林地面積（單位：hm2）；X 代表降水量，愛(ài)輝區(qū)的平均降水量為550 mm；E 為林分蒸散量（單位：mm/a），依據(jù)邴龍飛等[17]人研究，用東北地區(qū)平均林地蒸散量代替愛(ài)輝區(qū)的林地蒸散量，取值為450.43 mm；Y 代表地表徑流量，式（2）由孫力等[18]研究地表徑流與降雨量的關(guān)系以及個(gè)人通過(guò)平均計(jì)算得出，算得地表徑流量Y 等于10.51 mm。

凈化水質(zhì)：

式中，US代表林地每年凈化水質(zhì)的價(jià)值（單位：元/a）；K 代表凈化水的費(fèi)用取值2.09 元/t[19]；A,X,E,Y 代表含義和取值與前文相同。

2） 保育土壤。林地對(duì)于土壤的生態(tài)作用主要有兩個(gè)方面，一方面是林地的根部對(duì)于土壤有保護(hù)，防止土壤流失的作用，即固土作用；另一方面則是保持土壤的肥力，即保肥作用。

固土作用：

式中，UG代表林地每年固土價(jià)值（單位：元/a）；A 代含義和取值與前文相同；CT代表挖取運(yùn)送單位體積土壤所需的費(fèi)用，取值為12.6 元/m3[16]；X2代表無(wú)林地土壤侵蝕模數(shù)，取值為77.6 t/hm2[16]；X1代表有林地土壤侵蝕模數(shù)，由于愛(ài)輝區(qū)的林地主要以針闊葉交混林為主[20]，故取值為3 t/hm2；θ 代表林地土壤容重，取值1.257 t/m3[19]。

保肥作用：

式中，UB代表林地每年保肥價(jià)值（單位：元/a）；N,P,K分別代表林地土壤中氮磷鉀的占比；愛(ài)輝區(qū)林地土壤中氮磷鉀含量取值分別為0.073%、0.073%、1.283%[14]；C1代表化肥品種磷酸二銨價(jià)格，取值為2 400 元/t[14]；β1、β2代表化肥品種磷酸二銨含氮磷的比例，取值分別為18%、46%；C2代表化肥品種氯化鉀價(jià)格，取值為2 200 元/t[14]；β3代表化肥氯化鉀中含鉀的比例，取值為47.02%；M 代表林地土壤有機(jī)質(zhì)含量，愛(ài)輝區(qū)林地土壤中有機(jī)質(zhì)的含量是1.723%[19]；C3代表有機(jī)質(zhì)的價(jià)值，取值為320 元/t[14]；A、X1、X2代表含義和取值與前文相同。

3） 固碳釋氧。林地中的植被、土壤動(dòng)物、微生物通過(guò)自身作用，一方面吸收CO2中的碳物質(zhì)，即固碳作用；另一方面則是將CO2中的氧原子以氧氣的形式釋放出來(lái)，即釋氧作用。

固碳作用：

式中，UC代表林地每年固碳價(jià)值（單位：元/a）；CC代表固碳價(jià)格，取值為1 200 元/t[14]；RC代表CO2碳的百分比，取值為27.27%[14]；BN表示林地凈生產(chǎn)力，根據(jù)王琪[21]的研究，計(jì)算得到愛(ài)輝區(qū)的林地凈生產(chǎn)力為3.292 t/(hm2·a)；FT代表林地單位面積土壤年固碳值，取值為2.29 t/(hm2·a)[19]；A 代表含義和取值與前文相同。

釋氧作用：

式中，UY代表林地每年釋氧價(jià)值（單位：元/a）；CY代表氧氣價(jià)格，取值為1 000 元/t[14]，A,BN代表含義和取值與前文相同。

4） 累積營(yíng)養(yǎng)物質(zhì)。林地植物在生化過(guò)程中，通過(guò)吸收大氣、水分、土壤中的氮磷鉀等營(yíng)養(yǎng)物質(zhì)存儲(chǔ)在自身體內(nèi)，這一方面能夠累積營(yíng)養(yǎng)物質(zhì)，另一方面能有效防止水體富營(yíng)養(yǎng)化和污染。計(jì)算公式為：

式中，UYY代表林地每年?duì)I養(yǎng)物質(zhì)累積價(jià)值（單位：元/a）；NYY,PYY,KYY，分別代表林地氮磷鉀的含量；愛(ài)輝區(qū)林地氮磷鉀的含量為0.736%,0.418%,0.033%[19]；A,BN,C1,C2,β1,β2,β3代表含義和取值與前文相同。

5） 凈化大氣環(huán)境。林地生態(tài)系統(tǒng)通過(guò)吸收、過(guò)濾等降低有害物質(zhì)（如氮氧化物、氟化物等）的含量，并可以滯塵，提供有益于生態(tài)的負(fù)離子等功能。

提供負(fù)離子：

式中，UFL代表林地每年提供負(fù)離子價(jià)值（單位：元/a）；H 代表樹(shù)林平均高度，針闊交混林平均高度10 m[19]；KFL代表產(chǎn)生負(fù)離子的費(fèi)用，取值為5.818 5 元/1 018 個(gè)；QFL代表林地中負(fù)離子含量，取值為3 900 個(gè)/cm3；L 代表負(fù)離子平均壽命，取值為10 min[16]；A 代表含義和取值與前文相同。

吸收氟化物：

式中，UF代表林地每年吸收氟化物價(jià)值（單位：元/a）；KF代表氟化物治理費(fèi)用，取值為0.69 元/kg[14]；QF代表林地單位面積年吸收氟化物量，取值為2.5 kg/(hm2·a)[19]；A 代表含義和取值與前文相同。

吸收氮氧化物：

式中，UDY代表林地每年吸收氮氧化物價(jià)值（單位：元/a）；KDY代表氮氧化物治理費(fèi)用，取值為0.63 元/kg[14]；QDY代表林地單位面積年吸收氮氧化物量，取值為6 kg/(hm2·a)[19]；A 代表含義和取值與前文相同。

滯塵作用：

式中，UZC代表林地每年滯塵價(jià)值（單位：元/a）；KZC代表滯塵費(fèi)用，取值為0.15 元/kg[14]；QZC代表林地單位面積年滯塵量，取值為22 kg/(hm2·a)[19]；A 代表含義和取值與前文相同。

6） 生物多樣性保護(hù)。林地的存在為生物提供生存和繁衍生息提供保障，確保生物物種的多樣性。計(jì)算公式為：

式中，UW代表林地每年生物保護(hù)性?xún)r(jià)值（單位：元/a）；SW代表林地單位面積生物多樣性保護(hù)價(jià)值，根據(jù)謝高地等[22]的研究，取值為2 884.6 元/hm2；A 代表含義和取值與前文相同。

7） 森林游憩。林地生態(tài)系統(tǒng)在自身繁衍的同時(shí)，也為人類(lèi)提供文化娛樂(lè)的場(chǎng)所，即游憩價(jià)值。計(jì)算公式為：

式中，UYQ代表林地每年提供的森林游憩價(jià)值（單位：元/a）；SYQ代表林地單位面積森林游憩價(jià)值，取值為1 132.6 元/hm2[22]；A 代表含義和取值與前文相同。

2.2.2 耕地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估

耕地生態(tài)系統(tǒng)服務(wù)價(jià)值的計(jì)算公式如下：

式中，UGD代表耕地每年提供的生態(tài)系統(tǒng)服務(wù)價(jià)值（單位：元/a）；SGD代表耕地單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值，取值為6 114.3 元/hm2[22]；G 代表耕地面積（單位：hm2）。

2.2.3 草地生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估

草地生態(tài)系統(tǒng)服務(wù)價(jià)值的計(jì)算公式如下：

式中，UCD代表草地每年提供的生態(tài)系統(tǒng)服務(wù)價(jià)值（單位：元/a）；SCD代表草地單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值，取值為6 406.5 元/hm2[22]；C 代表草地面積（單位：hm2）。

2.2.4 水體生態(tài)系統(tǒng)服務(wù)價(jià)值評(píng)估

水體生態(tài)系統(tǒng)服務(wù)價(jià)值的計(jì)算公式如下：

式中，UST代表水體每年提供的生態(tài)系統(tǒng)服務(wù)價(jià)值（單位：元/a）；SST代表水體單位面積的生態(tài)系統(tǒng)服務(wù)價(jià)值，取值為4 0 676.4 元/hm2[22]；S 代表水體面積（單位：hm2）。

3 結(jié)果分析與討論

3.1 土地利用及變化分析

利用ENVI 軟件的精度驗(yàn)證計(jì)算功能，得到三期影像的土地利用總體分類(lèi)精度分別為88.38%、89.21%和87.9%，符合精度要求。黑河市愛(ài)輝區(qū)三期土地利用遙感分類(lèi)結(jié)果如圖2 所示。

圖2 愛(ài)輝區(qū)2000、2007、2013 年土地利用分類(lèi)結(jié)果圖

總體上看，愛(ài)輝區(qū)的土地利用類(lèi)型中林地所占的比例最大，占總地類(lèi)面積的一半以上；其次是草地，占總地類(lèi)的1/3 左右；耕地約占總地類(lèi)面積的10%；水體、其他地類(lèi)所占比例較少。林地是生態(tài)系統(tǒng)服務(wù)的主要部分，這符合愛(ài)輝區(qū)作為國(guó)家級(jí)生態(tài)示范區(qū)的標(biāo)準(zhǔn)。

1）耕地面積總體呈下降趨勢(shì)，但下降幅度有所差異。2000-2007 年耕地面積下降8.48%，即15 412.14 hm2；2007-2013 年耕地面積下降18.17%，即30 235.05 hm2。鑒于三期影像數(shù)據(jù)時(shí)間間隔相差不多但下降幅度卻明顯不同的情況，在查閱相關(guān)資料后可見(jiàn)，2000-2007 年間，耕地減少的原因一部分在于農(nóng)村勞動(dòng)力的減少而遭到荒棄，另一部分原因則是被建設(shè)用地占用。而2007-2013 年間，除上述原因外，“退耕還林”等生態(tài)措施同時(shí)影響著耕地面積的變化。

2）林地面積總體變化呈現(xiàn)先減后增趨勢(shì)。2000-2007 年林地面積下降2.11%，即16 243.83 hm2；2007-2013 年愛(ài)輝區(qū)林地面積增加5.11%，即38 492.19 hm2。2000-2007 年，為追求經(jīng)濟(jì)的快速發(fā)展，忽視了林地對(duì)生態(tài)的保護(hù)作用，大量地砍伐木材，雖然林地的面積只減少2.11%，但由于基數(shù)大，故而林木遭到大面積破壞。而2013 年的林地面積卻有大幅增加，與黑河市相繼開(kāi)展的“三年造林會(huì)戰(zhàn)”和“五年綠化”密切相關(guān)[23]。

3）草地面積變化呈現(xiàn)先增后減趨勢(shì)。2000-2007 年草地面積上升6.04%，即27 947.88 hm2；而2007-2013 年草地面積同比下降7.19%。從耕地和林地的變化原因可以分析得出，2000-2007 年草地面積增加一是由于大部分被荒棄的耕地由于缺乏耕作，逐漸被草地取代；二是由于被砍伐的林地直接接觸光照，從而生長(zhǎng)出大量的草體。2007-2013 年草地面積大量減少的原因同樣有兩個(gè)，一是“造林會(huì)戰(zhàn)”使得大量草地變?yōu)榱值?；二是草地作為后備土地?chǔ)蓄資源被開(kāi)發(fā)利用掉。

4）水體的變化呈現(xiàn)先減少后大量增加的現(xiàn)象，2000 年比2007 年水體面積減少7.83%，即812.79 hm2；而2013 年比2007 水體面積增加53.19%，即5 091.03 hm2。在查閱文獻(xiàn)資料后，認(rèn)為造成這一現(xiàn)象的原因主要是：①愛(ài)輝區(qū)屬于寒溫帶大陸性氣候，在我國(guó)極其少見(jiàn)，這種氣候會(huì)引起年降水量極不穩(wěn)定。②根據(jù)相關(guān)資料顯示[24-26]，2007 年黑龍江省出現(xiàn)歷史性干旱，高溫少雨，降水量比常年降低一到三成，2007-06-11～2007-08-10，平均降雨量?jī)H為143 mm，公別拉河蓄水量持續(xù)降低。③2013 年夏季，黑河市出現(xiàn)持續(xù)性大暴雨，發(fā)生了自1984 年以來(lái)最大的洪澇災(zāi)害，黑龍江流域以及各條河流水位大幅度上漲[27-29]，因此2013 年愛(ài)輝區(qū)的水體面積出現(xiàn)明顯增加。

5）其他地類(lèi)從2000-2013 年呈逐年增多趨勢(shì)，2000-2007 年愛(ài)輝區(qū)其他地類(lèi)面積增加57.82%，即4 520.88 hm2；2007-2013 年愛(ài)輝區(qū)其他地類(lèi)面積增加178.04%，即21 970.08 hm2。主要原因是由于隨著城鎮(zhèn)化進(jìn)程的加快，大量的土地被開(kāi)發(fā)成建設(shè)用地。

3.2 生態(tài)系統(tǒng)服務(wù)價(jià)值分析

愛(ài)輝區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值總量呈現(xiàn)先下降后上升的趨勢(shì)（表1）。2000-2007 年，價(jià)值總量下降1.86%，即10.6 億元；2007-2013 年，價(jià)值總量上升4.33%，即24.28 億元。13 a 間生態(tài)系統(tǒng)服務(wù)價(jià)值總體上升了2.40%，即13.68 億元。愛(ài)輝區(qū)各地類(lèi)單位面積生態(tài)服務(wù)價(jià)值表為林地單位面積生態(tài)服務(wù)價(jià)值最高，為6.84 萬(wàn)元/hm2；水體次之，為4.07 萬(wàn)元/hm2；草地為0.64 萬(wàn)元/hm2；耕地的單位面積生態(tài)服務(wù)價(jià)值最低，為0.61 萬(wàn)元/hm2。

表1 愛(ài)輝區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值表（2000、2007、2013 年）

1）耕地生態(tài)系統(tǒng)服務(wù)價(jià)值分析。耕地的生態(tài)系統(tǒng)服務(wù)價(jià)值從整體上看是連續(xù)降低的，下降率分別是8.54%和18.19%，生態(tài)服務(wù)價(jià)值在13 a 間下降了2.8 億元。由于耕地單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值相同，所以其生態(tài)系統(tǒng)服務(wù)價(jià)值變化與土地利用類(lèi)型變化是一致的。

2）林地生態(tài)系統(tǒng)服務(wù)價(jià)值分析。2000-2007 年，林地生態(tài)系統(tǒng)服務(wù)價(jià)值下降率為2.11%，即11.11 億元；2007-2013 年，上升率為5.11%，即26.33 億元。從表1 可以看出，林地的生態(tài)系統(tǒng)服務(wù)價(jià)值直接影響到愛(ài)輝區(qū)生態(tài)系統(tǒng)服務(wù)總價(jià)值。這是因?yàn)榱值氐纳鷳B(tài)系統(tǒng)服務(wù)價(jià)值在所有地類(lèi)中所占比重最大，占總價(jià)值的91.87%以上，林地的面積不僅在所有地類(lèi)中占比最大，單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值也最高。林地單位面積生態(tài)系統(tǒng)服務(wù)價(jià)值高主要體現(xiàn)在其功能類(lèi)別上（如圖3 所示）。其中保育土壤單位面積生態(tài)服務(wù)價(jià)值最高，為4.79 萬(wàn)元/hm2；其后依次是固碳釋氧、涵養(yǎng)水源、生物多樣性保護(hù)、森林游憩、累積營(yíng)養(yǎng)物質(zhì)；最低的是凈化大氣環(huán)境，單位面積生態(tài)服務(wù)價(jià)值0.01 萬(wàn)元/hm2。根據(jù)林地生態(tài)系統(tǒng)服務(wù)價(jià)值表所示（表2），2000-2013 年，林地的生態(tài)功能類(lèi)別和指標(biāo)類(lèi)型的生態(tài)服務(wù)價(jià)值變化和林地的土地利用類(lèi)型變化也呈一致。

圖3 愛(ài)輝區(qū)林地各項(xiàng)生態(tài)系統(tǒng)服務(wù)價(jià)值

表2 林地各生態(tài)功能類(lèi)別價(jià)值（2000 年、2007 年、2013 年）

3）草地生態(tài)系統(tǒng)服務(wù)價(jià)值分析。草地的生態(tài)服務(wù)價(jià)值與草地土地利用類(lèi)型變化趨勢(shì)一樣，先增加后減少。2000-2007 年，草地生態(tài)系統(tǒng)服務(wù)價(jià)值由29.67 億元上升為31.46 億元，上升率為6.03%；2007-2013 年，生態(tài)系統(tǒng)服務(wù)價(jià)值下降至29.19 億元，下降率為7.22%。草地生態(tài)系統(tǒng)服務(wù)價(jià)值總體下降0.48 億元，下降率為1.62%?？傮w上草地生態(tài)系統(tǒng)服務(wù)價(jià)值呈下降趨勢(shì)。

4）水體生態(tài)系統(tǒng)服務(wù)價(jià)值分析。水體的生態(tài)服務(wù)價(jià)值變化趨勢(shì)為先減少后增加。2000-2007 年，水體生態(tài)系統(tǒng)服務(wù)價(jià)值由4.22 億元下降為3.89 億元，下降率為7.82%；2007-2013 年，生態(tài)系統(tǒng)服務(wù)價(jià)值上升至5.96 億元，上升率為53.21%；13 a 間總體上升41.23%。由土地利用變化分析可以得知，2007 年發(fā)生的干旱與2013 年發(fā)生的洪澇災(zāi)害都對(duì)水體生態(tài)系統(tǒng)服務(wù)價(jià)值造成了明顯的影響。

4 結(jié) 語(yǔ)

本文通過(guò)對(duì)黑河市愛(ài)輝區(qū)2000 年、2007 年、2013 年三期遙感影像進(jìn)行土地利用分類(lèi)，采用中國(guó)陸地生態(tài)系統(tǒng)服務(wù)價(jià)值測(cè)量方法及森林生態(tài)系統(tǒng)服務(wù)功能評(píng)估規(guī)范核算了該區(qū)域的生態(tài)系統(tǒng)服務(wù)價(jià)值，可以得出以下結(jié)論：①愛(ài)輝區(qū)2000 年、2007 年和2013 年的生態(tài)系統(tǒng)服務(wù)價(jià)值分別為571.03 億元、560.43 億元、584.71 億元。生態(tài)系統(tǒng)服務(wù)總價(jià)值呈先減少后增加的趨勢(shì)，總體增加13.68 億元。②一個(gè)地區(qū)生態(tài)系統(tǒng)服務(wù)價(jià)值與土地利用類(lèi)型有著密切的關(guān)系。在對(duì)生態(tài)系統(tǒng)服務(wù)價(jià)值貢獻(xiàn)較少的其他地類(lèi)面積連年增加、草地和耕地面積減少的情況下，愛(ài)輝區(qū)的生態(tài)系統(tǒng)服務(wù)價(jià)值的增長(zhǎng)得益于林地面積和水體面積的增加。③從片面地追求城鎮(zhèn)化速度，忽視林地生態(tài)系統(tǒng)服務(wù)價(jià)值到2011 年被命名為“國(guó)家級(jí)生態(tài)示范區(qū)”，愛(ài)輝區(qū)形成了人與自然和諧相處的新型現(xiàn)代化生態(tài)城鎮(zhèn)格局。