大興安嶺天保工程區(qū)生態(tài)系統(tǒng)服務(wù)變化研究

　　摘要：利用Landsat和MODIS遙感數(shù)據(jù)和InVEST模型，評(píng)估了1990―2015年大興安嶺天保工程區(qū)實(shí)施前后生態(tài)系統(tǒng)土壤風(fēng)蝕量、生境質(zhì)量和碳儲(chǔ)量的時(shí)空變化，并通過(guò)地理探測(cè)器方法對(duì)生態(tài)系統(tǒng)服務(wù)時(shí)空變化分布與趨勢(shì)進(jìn)行驅(qū)動(dòng)因子探測(cè)和交互作用探測(cè)分析。結(jié)果表明：①?gòu)目臻g分布來(lái)看，額爾古納市南部、鄂溫克族自治旗西北部、科爾沁右翼前旗的土壤風(fēng)蝕量大，生境質(zhì)量差，碳儲(chǔ)量低。②2000年天保工程實(shí)施之后，森林面積顯著增加，土壤風(fēng)蝕量減少1.14%，生境質(zhì)量增加0.49%，碳儲(chǔ)量減少幅度為0.12%，低于2000年以前的0.64%。③土壤風(fēng)蝕量主要受植被類(lèi)型、土壤類(lèi)型的影響，生態(tài)系統(tǒng)類(lèi)型、植被類(lèi)型對(duì)生境質(zhì)量、碳儲(chǔ)量空間異質(zhì)性的影響顯著大于其他因素，海拔、坡度2種地形因子對(duì)土壤風(fēng)蝕量、碳儲(chǔ)量的解釋能力很小。溫度、降水2種氣候因子對(duì)土壤風(fēng)蝕量、生境質(zhì)量和碳儲(chǔ)量的解釋能力不固定，但仍是不可或缺的因素。

　　關(guān)鍵詞：天然林保護(hù)工程;生態(tài)系統(tǒng)服務(wù);遙感數(shù)據(jù);地理探測(cè)器;大興安嶺地區(qū)

　　森林生態(tài)系統(tǒng)是保證健康環(huán)境的重要部分，森林提供的生態(tài)系統(tǒng)服務(wù)是維持社會(huì)經(jīng)濟(jì)發(fā)展和健康生態(tài)系統(tǒng)的基礎(chǔ)[1]。由于社會(huì)經(jīng)濟(jì)發(fā)展與人口增長(zhǎng)，土地退化、土壤侵蝕和毀林造成的生態(tài)環(huán)境惡化問(wèn)題一直困擾著中國(guó)生態(tài)安全[2]。為了減少這些問(wèn)題的影響，國(guó)務(wù)院2000年實(shí)施了世界上最大和最全面的生態(tài)保護(hù)工程—天然林保護(hù)工程(天保工程)，天保工程森林面積占中國(guó)天然林面積的53%，主要目標(biāo)是解決天然林的恢復(fù)發(fā)展問(wèn)題，保護(hù)環(huán)境與促進(jìn)社會(huì)和諧發(fā)展[3]。眾多學(xué)者對(duì)國(guó)內(nèi)天保工程生態(tài)系統(tǒng)服務(wù)進(jìn)行了評(píng)析。范琳等[4]通過(guò)建立森林生態(tài)系統(tǒng)綜合評(píng)價(jià)指標(biāo)體系，對(duì)山西省天保工程進(jìn)行了效益評(píng)價(jià)。肖中琪等[5]結(jié)合遙感與隨機(jī)森林，分析了新疆天山西部天保工程生態(tài)系統(tǒng)類(lèi)型、森林生物量的時(shí)空變化。

　　黃繼紅等[6]選取6項(xiàng)服務(wù)、9項(xiàng)指標(biāo)對(duì)陜西省寧陜縣天保工程的顯著成效進(jìn)行了分析。傳統(tǒng)的天保工程生態(tài)系統(tǒng)服務(wù)研究側(cè)重于及價(jià)值量及時(shí)空變化研究，但是很少有人對(duì)其中的原因進(jìn)行定量分析。定量分析生態(tài)系統(tǒng)服務(wù)對(duì)生態(tài)保護(hù)與修復(fù)、區(qū)域可持續(xù)發(fā)展具有科學(xué)指導(dǎo)意義。大興安嶺天保工程區(qū)自然條件優(yōu)越，是中國(guó)重要的生態(tài)屏障，是中國(guó)木材的重要產(chǎn)區(qū)和戰(zhàn)略資源儲(chǔ)備基地;是黑龍江水系的源頭，穩(wěn)定的森林生態(tài)系統(tǒng)為中下游地區(qū)提供了寶貴的工農(nóng)業(yè)生產(chǎn)和生活用水。以往對(duì)大興安嶺天保工程區(qū)研究多集中在生態(tài)系統(tǒng)類(lèi)型變化研究[7]，生態(tài)系統(tǒng)服務(wù)時(shí)空變化格局分析[8]，或者單一生態(tài)系統(tǒng)服務(wù)時(shí)空變化與原因分析[9]。

　　目前缺乏對(duì)天保工程實(shí)施前后生態(tài)系統(tǒng)服務(wù)的全面評(píng)估與原因分析。本研究利用長(zhǎng)期定位觀(guān)測(cè)、遙感數(shù)據(jù)與InVEST模型，對(duì)比大興安嶺天保工程區(qū)生態(tài)系統(tǒng)類(lèi)型變化，從土壤風(fēng)蝕量、生境質(zhì)量、碳儲(chǔ)量3個(gè)方面評(píng)估1990―2015年天然林保護(hù)工程實(shí)施前后生態(tài)系統(tǒng)服務(wù)演變情況，探討生態(tài)系統(tǒng)服務(wù)時(shí)空變化的影響因素。以期為今后工程區(qū)規(guī)劃和行動(dòng)、自然資源合理開(kāi)發(fā)利用、區(qū)域生態(tài)文明建設(shè)提供重要的科學(xué)基礎(chǔ)。

　　1研究數(shù)據(jù)與方法

　　1.1 研究區(qū)概況

　　大興安嶺天保工程區(qū)位于118°14′~127°44′E、45°13′~53°56′N，行政范圍包括內(nèi)蒙古自治區(qū)興安盟、呼倫貝爾市和黑龍江省大興安嶺地區(qū)，面積約為26萬(wàn)km2。屬于溫帶大陸性季風(fēng)氣候，冬季漫長(zhǎng)寒冷，夏季短暫多雨，太陽(yáng)輻射資源豐富，全年均溫−2.5~6.8℃，年均降水量310~750mm，該區(qū)跨寒溫帶針葉林、溫帶針闊葉混交林和溫帶草原三大植被區(qū)域，植被類(lèi)型以森林為主，其中主要分布針葉林、闊葉林以及針闊混交林，是中國(guó)的重點(diǎn)林區(qū)。

　　1.2 數(shù)據(jù)源及處理

　　本文涉及的數(shù)據(jù)包括：①1990、2000、2015年的生態(tài)系統(tǒng)類(lèi)型數(shù)據(jù)來(lái)源于國(guó)家地球系統(tǒng)科學(xué)數(shù)據(jù)中心

　　，三期土地覆被數(shù)據(jù)以L(fǎng)andsat為數(shù)據(jù)源，采用面向?qū)ο蠓椒ǖ玫剑瑢⒀芯繀^(qū)劃分為森林、草地、農(nóng)田、濕地、城鎮(zhèn)及裸地6個(gè)一級(jí)類(lèi)，落葉針葉林、針闊混交林、草叢、湖泊、河流等40個(gè)二級(jí)類(lèi)，空間分辨率為30m，利用ArcGIS10.6裁剪，得到大興安嶺地區(qū)土地覆被數(shù)據(jù)。②30m分辨率數(shù)字高程模型(DEM)數(shù)據(jù)來(lái)源于美國(guó)國(guó)家航空航天局(https://search.asf.alaska.edu/#/)，利用ArcGIS10.6提取高程、坡度、坡向等地形因子。③土壤數(shù)據(jù)：土壤質(zhì)地?cái)?shù)據(jù)來(lái)源于地理空間數(shù)據(jù)云，土壤深度數(shù)據(jù)來(lái)自于全國(guó)第二次土壤普查數(shù)據(jù)集。

　�、躆ODISNDVI數(shù)據(jù)，來(lái)自MOD13A3的植被指數(shù)數(shù)據(jù)集，時(shí)間分辨率為月，空間分辨率為250m×250m，來(lái)源自NASA/EOSLPDAAC數(shù)據(jù)分發(fā)中心的MODIS產(chǎn)品，湖泊分布、河流分布、居民地和道路分布數(shù)據(jù)均采用ArcGIS10.6從生態(tài)系統(tǒng)類(lèi)型數(shù)據(jù)中提取。⑤氣象數(shù)據(jù)：日最高均溫、日最低均溫來(lái)自于中國(guó)氣象數(shù)據(jù)網(wǎng)，采用克里格插值的方法生成1km柵格數(shù)據(jù)，年平均降水、年平均氣溫、植被類(lèi)型、土壤類(lèi)型空間分布柵格數(shù)據(jù)來(lái)自于中科院資源環(huán)境數(shù)據(jù)中心，用于影響因素分析。參考相關(guān)文獻(xiàn)[9~14]，并結(jié)合大興安嶺天保工程區(qū)實(shí)際情況，選取生態(tài)系統(tǒng)類(lèi)型、降水量、溫度、植被類(lèi)型、海拔、坡度、土壤類(lèi)型作為影響因素。

　　1.3 生態(tài)系統(tǒng)服務(wù)變化評(píng)估方法

　　1.3.1土壤風(fēng)蝕量土壤風(fēng)蝕量采用修正風(fēng)蝕方程(RWEQ)進(jìn)行估算。在充分考慮氣候條件、植被狀況、地表粗糙度、土壤可蝕性、土壤結(jié)皮等要素情況下，RWEQ模型可較好定量評(píng)估土壤侵蝕模數(shù)[10]。

　　1.3.2生境質(zhì)量InVEST模型生物多樣性模塊采用生境質(zhì)量表征生物多樣性維持。對(duì)不同因子賦權(quán)重以確定不同因子對(duì)生境質(zhì)量的影響程度。采用因子層次分析法對(duì)天保工程區(qū)進(jìn)行生境質(zhì)量評(píng)價(jià)。

　　2結(jié)果分析

　　2.1 生態(tài)系統(tǒng)類(lèi)型變化特征

　　大興安嶺天保工程區(qū)生態(tài)系統(tǒng)類(lèi)型包括森林、草地、濕地、農(nóng)田、城鎮(zhèn)和其它，以森林為主，森林比另外5種生態(tài)系統(tǒng)的分割度低，聚合程度高，連通性好，主要分布在大興安嶺中心地區(qū)。1990―2015年，森林面積呈現(xiàn)先減少后增加的趨勢(shì)。1990―2000年，森林面積減少了0.29%，2000年森林面積為173632.71km2，主要是向農(nóng)田轉(zhuǎn)化，轉(zhuǎn)出面積為895.50km2，占總轉(zhuǎn)出面積的41.21%，面積減少的生態(tài)系統(tǒng)類(lèi)型還有草地、濕地。

　　2015年，森林面積達(dá)175185.50km2，增加了0.89%，2000―2015年，森林主要是由濕地、草地轉(zhuǎn)化而來(lái)，分別占總轉(zhuǎn)入面積的44.47%、30.13%，面積減少的生態(tài)系統(tǒng)類(lèi)型主要有草地、濕地。大興安嶺天保工程實(shí)施之后，森林生態(tài)系統(tǒng)得到很大程度的恢復(fù)與改善，森林總面積不斷增加。1990―2015年，農(nóng)田、城鎮(zhèn)、其它生態(tài)系統(tǒng)面積一直呈增加趨勢(shì)。

　　2.2 生態(tài)系統(tǒng)服務(wù)時(shí)空分布特征

　　從空間分布來(lái)看，3個(gè)時(shí)期研究區(qū)生態(tài)系統(tǒng)的土壤風(fēng)蝕量、生境質(zhì)量、碳儲(chǔ)量表現(xiàn)出一致性與差異性，額爾古納市南部、鄂溫克族自治旗西北部、科爾沁右翼前旗的土壤風(fēng)蝕量大，生境質(zhì)量差，碳儲(chǔ)量低。在時(shí)空分布上，土壤風(fēng)蝕量呈現(xiàn)先增加后減少的趨勢(shì)，1990、2000、2015年的土壤風(fēng)蝕量均值分別為912.49、1084.08、1071.71t/km2。1990―2000年，土壤風(fēng)蝕量增加了171.59t/km2，占比18.80%，增加區(qū)域主要分布在呼瑪縣東部、額爾古納市、牙克石市以及鄂溫克族自治旗北部。

　　2000―2015年，土壤風(fēng)蝕量減少了12.37t/km2，占比1.14%，減少區(qū)域主要分布在額爾古納市南部、鄂溫克自治旗、科爾沁右翼前旗中部及東南部。生境質(zhì)量呈現(xiàn)先減少后增加的趨勢(shì)，1990、2000、2015年生境質(zhì)量平均值分別為60.00、59.03、59.32。1990―2000年，研究區(qū)四周生境質(zhì)量有變好趨勢(shì)，中部呈現(xiàn)變差趨勢(shì)，總體變差，占比1.62%，2000―2015年，生境質(zhì)量整體有變好趨勢(shì)，增加0.49%。1990―2015年研究區(qū)內(nèi)大部分面積總碳儲(chǔ)量保持不變，從1990年的6222.74Tg到2000年的6190.44Tg再到2015年的6182.86Tg，大興安嶺天保工程區(qū)平均碳儲(chǔ)量則呈現(xiàn)出一直下降趨勢(shì)。但從下降幅度0.64%到0.12%，2000―2015年碳儲(chǔ)量的下降幅度比1990―2010年的緩慢。

　　2.3 生態(tài)系統(tǒng)服務(wù)時(shí)空分布影響因素

　　除海拔、坡度與土壤風(fēng)蝕量，坡度與碳儲(chǔ)量不顯著相關(guān)外，其余因素與生態(tài)系統(tǒng)服務(wù)均顯著相關(guān)(P<0.01)。植被類(lèi)型、土壤類(lèi)型對(duì)土壤風(fēng)蝕量的解釋能力明顯大于其他因子，生態(tài)系統(tǒng)類(lèi)型、植被類(lèi)型對(duì)生境質(zhì)量、碳儲(chǔ)量的解釋能力明顯大于其他因子。海拔、坡度2種地形因子對(duì)土壤風(fēng)蝕量、碳儲(chǔ)量的解釋能力很小。由于年際間氣候處于動(dòng)態(tài)變化中，溫度、降水2種氣候因子對(duì)土壤風(fēng)蝕量、生境質(zhì)量和碳儲(chǔ)量的解釋能力不固定，但仍是不可或缺的因素。

　　在因子間的交互作用方面，各因子表現(xiàn)出較強(qiáng)的交互作用，對(duì)土壤風(fēng)蝕量解釋能力最大的是植被類(lèi)型協(xié)同土壤類(lèi)型，各年度解釋能力達(dá)到0.26~0.29;除生態(tài)系統(tǒng)類(lèi)型與植被類(lèi)型、生態(tài)系統(tǒng)類(lèi)型與土壤類(lèi)型外，其余因素兩兩交互作用對(duì)生境質(zhì)量各年度的解釋能力均超過(guò)0.99;生態(tài)系統(tǒng)類(lèi)型協(xié)同降水量、溫度、植被類(lèi)型、土壤類(lèi)型均能對(duì)碳儲(chǔ)量達(dá)到很好的解釋作用，各年度解釋力達(dá)到0.89~0.93。

　　3討論

　　自2000年大興安嶺天保工程實(shí)施以來(lái)，生態(tài)系統(tǒng)類(lèi)型發(fā)生較大變化，森林面積顯著提高，生態(tài)系統(tǒng)服務(wù)能力得到改善。1990―2015年，土壤風(fēng)蝕量呈現(xiàn)先增加后減少的趨勢(shì)，相反，生境質(zhì)量先減少后增加，碳儲(chǔ)量則一直減少，但2000年后的減少幅度小于2000年前。說(shuō)明大興安嶺天保工程的實(shí)施對(duì)生態(tài)系統(tǒng)有顯著成效。但是對(duì)比1990年，2015年的土壤風(fēng)蝕量、生境質(zhì)量、碳儲(chǔ)量分別增加159.22t/km2，減少0.68，減少39.88Tg，恢復(fù)能力較慢，沒(méi)有恢復(fù)到1990年的水平，這與大興安嶺地區(qū)氣候寒冷，生態(tài)系統(tǒng)服務(wù)恢復(fù)能力較其他溫暖濕潤(rùn)區(qū)較低有關(guān)。

　　土壤論文范例：綠肥不同還田方式對(duì)土壤溫室氣體排放的影響

　　4結(jié)論

　　1)從空間分布來(lái)看，額爾古納市南部、鄂溫克族自治旗西北部、科爾沁右翼前旗的土壤風(fēng)蝕量大，生境質(zhì)量差，碳儲(chǔ)量低。

　　2)2000年天保工程實(shí)施之后，森林面積顯著增加，土壤風(fēng)蝕量減少1.14%，生境質(zhì)量增加0.49%，碳儲(chǔ)量減少幅度為0.12%，低于2000年前的0.64%。說(shuō)明天保工程的實(shí)施有助于生態(tài)系統(tǒng)服務(wù)的恢復(fù)與改善。

　　3)根據(jù)地理探測(cè)器統(tǒng)計(jì)結(jié)果q值分析，土壤風(fēng)蝕量主要受植被類(lèi)型、土壤類(lèi)型的影響，生態(tài)系統(tǒng)類(lèi)型、植被類(lèi)型對(duì)生境質(zhì)量、碳儲(chǔ)量空間異質(zhì)性的影響顯著大于其他因素，是時(shí)空變化的重要因素。海拔、坡度2種地形因子對(duì)土壤風(fēng)蝕量、碳儲(chǔ)量的解釋能力很小。溫度、降水2種氣候因子對(duì)土壤風(fēng)蝕量、生境質(zhì)量和碳儲(chǔ)量的解釋能力不固定，但仍是不可或缺的因素。

　　參考文獻(xiàn)(References)：

　　CaoSX,LiuYJ,SuWetal.ThenetecosystemservicesvalueinmainlandChina[J].ScienceChinaEarthSciences,2018,61:595-603.

　　作者：鄭樹(shù)峰1,2，王麗萍2，臧淑英1

轉(zhuǎn)載請(qǐng)注明來(lái)自發(fā)表學(xué)術(shù)論文網(wǎng)：http:///nylw/27615.html