天山中部云杉天然林水源涵養功能定量評估
——以烏魯木齊河流域為例

丁程鋒,張繪芳,李 霞,李偉濤，高亞琪,*

1 新疆林業科學院現代林業研究所，烏魯木齊 830000 2 新疆農業大學草業與環境科學學院，烏魯木齊 830052 3 滁州學院地理信息與旅游學院, 滁州 239000

?

天山中部云杉天然林水源涵養功能定量評估
——以烏魯木齊河流域為例

丁程鋒1,2,張繪芳1,李 霞2,李偉濤3，高亞琪1,*

1 新疆林業科學院現代林業研究所，烏魯木齊 830000 2 新疆農業大學草業與環境科學學院，烏魯木齊 830052 3 滁州學院地理信息與旅游學院, 滁州 239000

為定量分析天山中部流域尺度云杉天然林水源涵養功能,以烏魯木齊河流域為研究區,使用InVEST模型,研究云杉林及其它主要地類水源涵養量的大小,并對云杉林在不同分布面積、不同地形因子條件下的水源涵養量變化進行定量分析。結果表明：①InVEST模型可較好的確定研究區流域尺度水源涵養量。基于InVEST模型的模擬,烏魯木齊河流域云杉林區水源涵養總量為4.93×106m3,占研究區水源涵養總量(2.41×107m3)的20.46%,林區平均水源涵養深度為54.25 mm；②云杉林區水源涵養量的大小在海拔、坡度、坡向上的變化與云杉林空間分布格局一致,每公頃水源涵養量隨海拔升高先增大后降低、隨坡度增大而降低；陰坡水源涵養能力最強；③研究區各地類中以林地水源涵養量最大(水源涵養量544.78 m3/hm2),隨著云杉林覆蓋率的不斷增加,水源涵養總量、平均水源涵養深度及單位面積涵養量均呈增加趨勢；④研究區水源涵養量的貢獻率與森林面積密切相關,云杉林面積每增加1%,研究區與云杉林區水源涵養總量分別增加0.437×106m3、0.522×106m3；加強對研究區云杉林的保護與撫育管理,才能使森林發揮持續穩定的水源涵養生態服務功能。

水源涵養功能；InVEST模型；云杉林；天山中部；烏魯木齊河流域

水源涵養是森林的重要生態功能,以往對森林水源涵養功能的研究大多基于林分尺度或坡面尺度[1],量化流域尺度森林水源涵養功能,是水資源管理的需求也是多年來森林水源涵養功能研究的瓶頸。為正確認識森林的水源涵養功能,眾多學者通過固定樣地調查、長期野外監測方式,采用水量平衡法[2]、蓄水能力法[3]、綜合評價法[4]、林冠截留法[5]等對林分尺度、坡面尺度下的云杉林水源涵養功能進行研究,但對流域尺度云杉林水源涵養功能的定量評估卻鮮見報道。位于烏魯木齊河流域中山帶的云杉天然林,擔負著涵養水源、調蓄山區降水、保護冰川安全和維護區域生態平衡的重大生態責任[6]。定量分析天山中部流域尺度云杉林變化條件下的水源涵養功能,對推進干旱區森林生態功能的定量評估及水資源的科學管理都具有重要意義。

隨著GIS技術與生態系統服務模型的發展,流域尺度水源涵養功能定量評估已成為可能[7]。InVEST(Integrated Valuation of Ecosystem Services and Tradeoffs)模型是生態系統服務模型中應用最成熟的模型之一,該模型不僅可將評價結果以地圖的形式直觀表達出來,并可以進行多尺度、多情景的分析[8]。目前已被廣泛應用于不同區域生態系統服務功能[9-11]、水質凈化[12]、土壤侵蝕與保持[13]、生境質量評價[14]等方面的研究。

自2010年起,InVEST模型逐漸被應用于森林水源涵養評價,傅斌等[15]利用InVEST模型中的水源涵養模型對四川都江堰市水源涵養功能進行了分析,并采用綜合指數法對研究區水源涵養的重要性進行了評價；余新曉等[16]以北京山區為研究區,在森林資源二類調查數據基礎上,對森林景觀類型進行了劃分,并利用InVEST模型對不同森林景觀類型的水源涵養功能進行了評價；王紀偉等[17]以漢江上游為研究區,利用InVEST模型對該區域各森林生態系統的水源涵養量進行了計算,并分析了該區域水源涵養功能的空間變化趨勢；白楊等[7]以白洋淀流域為研究區,通過建立森林變化情景,利用InVEST模型對不同情景下的水源涵養功能進行了評價；李明陽等[18]基于實測數據和遙感數據,利用InVEST模型,對武夷山國家自然保護區內森林土壤涵養水源的能力進行了分析。上述研究為使用InVEST模型進行新疆天山中部云杉天然林森林水源涵養功能定量評估提供了思路。本文在確定研究區云杉林分布范圍變化規律的基礎上,基于流域多年平均降水與蒸散發數據,利用InVEST模型,綜合考慮氣候、土壤、地形等因素,定量評估流域尺度云杉林的水源涵養功能,為水資源有效利用及更大尺度水源涵養功能預測與評估提供依據。

1 研究區概況

研究區位于天山中部86°45′—87°56′ E,43°00′—44°07′ N間的烏魯木齊河流域(圖1)。山區(西白楊溝口以上)流域集水面積1070 km2,其中冰川面積約38 km2,流域山體高大,地形陡峻,平均海拔3006 m,平均坡度48.5°,年徑流量2.37億m3。流域南起烏魯木齊河上游天山山脈的依連哈比尕山分水嶺,北至古爾班通古特沙漠南緣的東道海子,西接頭屯河流域,東到烏拉泊和柴窩堡洼地之間的分水嶺,由南至北各支流泉溝匯入東道海子[19]。流域屬典型大陸性氣候,多年平均氣溫2.1 ℃,多年平均降水449.2 mm,降水年內分配不均,主要集中在夏季,5—9月占全年降水量比重最大,約為70%。植被分布以中山帶雪嶺云杉(PiceaschrenkianaFisch. et Mey.)為主,研究區內云杉林面積約90.85 km2(圖2)。流域內土壤垂直地帶變化明顯,依次為高山草甸土、灰褐色森林土、黑鈣土、山地栗鈣土。

圖1 研究區示意圖Fig.1 The location of the study area

圖2 云杉林空間分布圖Fig.2 The spatial distribution map of spruce forest

2 數據與方法

2.1 數據來源與處理

本研究所需主要數據有：多年平均降水量、多年潛在蒸散發、土壤深度、土壤可利用有效水、土壤飽和導水率、土地利用類型圖、流域及子流域邊界等,研究所需其數據來源及處理方法見表1。

表1 數據來源及處理方法

2.2 方法

2.2.1 產水量模擬

InVEST模型的產水量模塊是基于水量平衡法開發而來,將每個柵格上的降水量減去實際蒸散發后的水量,以柵格為單元定量評價不同地塊的產水能力[16]。計算公式如下：

(1)

式中,Yxj為森林類型j中單元格x上的年產水量(mm)；AETxj為森林類型j中單元格x上的年實際蒸散量(mm)；Px為單元格x上的年降水量(mm)。

(2)

式中,Rxj為土地利用類型j上柵格單元x的干燥指數,無量綱,表示潛在蒸發量與降雨量的比值。

(3)

式中,k為作物系數,即蒸散系數,是作物蒸散量與潛在蒸散量的比值；ET0為潛在蒸散發量。

ωx為修正植被年可利用水量與降水量的比值,無量綱：

(4)

式中,Z為zhang系數,是表征多年平均降水特征用的一個常數,是模型的關鍵參數,默認值是9.433；AWCx為可利用水。

2.2.2 水源涵養計算

用InVEST模型中的產水量模塊計算年產水量之后,根據DEM計算地形指數D,考慮不同土壤的滲透性、不同土地利用類型地表徑流流速系數的影響,最后計算水源涵養量[16],計算公式如下：

WR= min (1, 249/V) × min (1, 0.9×D/3) × min (1,Ksoil/300) ×Y

(5)

式中,WR為多年平均水源涵養量(mm)；D為地形指數,無量綱；Ksoil為土壤飽和導水率(cm/d)；V為流速系數,由相關文獻獲得[15]；Y為產水量。

(6)

式中,Watershed pixel count為集水區柵格數量,無量綱；Soil depth為土壤深度(mm)；Percent slope為百分比坡度。

2.2.3 地形因子區間劃分

為定量評估云杉林區水源涵養與地形因子關系,對海拔、坡向和坡度進行區間劃分,并利用ArcGIS空間統計工具統計各區間范圍的水源涵養量。其中,海拔劃分以1600 m為起始基準,每增加200 m為1級,共劃分為7級：1600—1800、1800—2000、2000—2200、2200—2400、2400—2600、2600—2800、2800—3000 m；坡向以正北方向為0°,按45°角,順時針旋轉劃分6個方向：陰坡(0°—45°,315°—360°)、東北坡(45°—90°)、東南坡(90°—135°)、陽坡(135°—225°)、西南坡(225°—270°)、西北坡(270°—315°)；坡度劃分為6個等級：平坡(<5°)、緩坡(5°—15°)、斜坡(15°—25°)、陡坡(25°—35°)、急坡(35°—45°)、險坡(>45°)。

2.2.4 情景模擬

研究區云杉林以中齡林為主,幼齡林分布較少[26],在高海拔地帶,由于氣候條件等多方面的原因,云杉林分布稀疏,生長不良；在低海拔地帶,因降水較少,氣溫較高,云杉林生長也不好,唯中海拔地帶降水充沛,氣候溫和,適宜云杉林生長[27]。為研究云杉林分布變化對水源涵養的影響,根據張毓濤、劉貴峰等[26,28]研究成果(云杉林主要分布在海拔1800—2700 m、坡度10°—50°的陰坡半陰坡區域)和天山云杉的生物學特征[27]設定云杉林變化的5種情景,分別利用InVEST模型對不同情景下的產水量進行模擬,根據2.2.2的方法計算各情景下的水源涵養量(各情景模擬過程中,其它數據不變,只改變土地利用類型數據),最后運用線性趨勢分析法分析不同情景間的水源涵養量差異。設定的云杉分布5種情景如下(圖3為云杉林現狀分布(情景1)、云杉林覆蓋最大(情景3)、無云杉林(情景5)的情景圖)：

情景1 現有云杉林實際分布,此時,云杉林覆蓋率為8.17%；

情景2 海拔1800—2700 m、坡度10°—50°區域內陰坡全部納為林地,此時,云杉林覆蓋率12.26%；

情景4 現有半陰坡區域的林地轉換為草地,此時,云杉林覆蓋率5.73%；

情景5 現有云杉林全部轉換為草地,此時,云杉林覆蓋率0%。

圖3 不同情景土地利用類型圖Fig.3 The land use map of different situations

3 結果與分析

3.1 研究區水源涵養量總體分析

3.1.1 水源涵養量空間分布規律

根據2.2.1的方法,得到研究區產水量,在此基礎之上,根據2.2.2的方法計算得到研究區水源涵養量空間分布圖(圖4)。由圖4得出,研究區多年平均水源涵養深度在0—118.54 mm之間,水源涵養總量為2.41×107m3。

水源涵養功能總體呈西高東低、北高南低的趨勢,水源涵養量較高的區域主要有兩個：①云杉林分布區,由于森林覆蓋,林下土壤層較厚,有利于土壤孔隙度發育,且土壤表層枯落物較多,形成良好的土壤結構及通風狀況,其土壤下滲、持水能力較強[15,17],有利于對水分的截留；②研究區中西部區域,因該區多年平均降水量較大、且草地廣泛分布所致。

2.3.1 色譜條件 色譜柱：辛烷基鍵合硅膠為填充劑（150×4.6mm，5μm）；流動相：0.01mol·L-1磷酸氫二鈉溶液（用磷酸調節 pH 值至 7.6）-乙腈（75∶25）；檢測波長：302nm；柱溫：30℃；流速：1.0mL·min-1。

水源涵養功能較差的區域主要有：①研究區南部高海拔區,該區為裸巖和永久冰川,幾乎無土壤、植被分布,降水截留能力較差；但該區域有少量水源涵養功能較高的地方,疊加土地利用類型圖可知,散布水源涵養較高的地方屬于草地分布區；②流域北部出山口,該區是人類活動集中區,建設用地和耕地的水源涵養能力較差[15],區域降水少、蒸發量大,故水源涵養量較低。

3.1.2 不同土地利用類型的水源涵養量分析

利用ArcGIS空間統計分析工具,統計不同土地利用類型每公頃水源涵養量,得到表2。由表2可知,不同土地利用類型水源涵養總量依次為：草地(1394.43×104m3)>裸地(508.08×104m3)>林地(494.72×104m3)>冰雪(35.47×104m3)>水體(1.98×104m3)>其它(0.39×104m3),分別占研究區水源涵養總量的57.26%、20.86%、20.32%、1.46%、0.08%、0.02%,水源涵養總量的大小與各種地類面積有很大關系[29],但不同土地類型的水源涵養總量的大小與其水源涵養能力的大小并無緊密關系[30]。

圖4 研究區水源涵養分布圖Fig.4 Water conservation in the study area

圖5 云杉林區水源涵養分布圖Fig.5 Water conservation in forest region

不同土地利用類型的水源涵養能力大小由每公頃水源涵養量來衡量,由表2知,各土地利用類型每公頃水源涵養量大小順序為：林地(544.78 m3/hm2)>草地(275.29 m3/hm2)>裸地(118.85 m3/hm2)>其它(88.23 m3/hm2)>冰雪(50.68 m3/hm2)>水體(23.41 m3/hm2),林地水源涵養能力最強,這主要是由于林地土壤孔隙度較大,且對土壤結構改善作用明顯。

表2 不同土地利用類型水源涵養量

3.2 云杉林區水源涵養量分析

3.2.1 云杉林區水源涵養量分析

在確定研究區水源涵養空間分布圖基礎上,以云杉林林區為邊界,裁剪得到林區水源涵養分布圖(圖5)。由圖5看出,云杉林林區多年平均水源涵養深度在2.4—118.54 mm之間；研究區林地覆蓋率8.17%,水源涵養總量為4.93×106m3,占研究區水源涵養總量的20.46%,林地平均水源涵養深度54.25 mm,云杉林區水源涵養功能呈現出北高南低、西高東低的趨勢。

3.2.2 云杉林區水源涵養與地形因子關系

在研究區,地形微環境影響水、土和熱量的分布,由于水熱環境的不同使森林分布格局、地表狀況及土壤狀況發生變化[17]。根據2.2.3的方法統計不同地形因子不同區間范圍內的水源涵養量,分別得到表3—表5。

表3 不同海拔水源涵養量

表4 不同坡度水源涵養量

表5 不同坡向水源涵養量

由表3可知,不同海拔云杉林水源涵養量的貢獻率分別為：0.05%、5.07%、20.19%、32.58%、31.11%、10.80%、0.21%,以海拔2200—2400 m處貢獻率最大,這與云杉林在不同海拔區間分布的面積比有很大關系；就每公頃水源涵養量來說,以海拔2000—2200 m處每公頃水源涵養能力最強,達593.47 m3/hm2,其次是2200—2400 m,為573.21 m3/hm2,海拔1600—1800 m與2800—3000 m處每公頃水源涵養能力較弱,分別為329.66、368.08 m3/hm2,這與不同海拔降水及森林分布有關。

由表4可知,不同坡度(平坡、緩坡、斜坡、陡坡、急坡、險坡)條件下,云杉林水源涵養量的貢獻率分別為：0.56%、6.06%、18.95%、35.31%、29.66%、9.46%,這與云杉林在不同坡度分布的面積比有很大關系；就單位面積水源涵養量來說,平坡(652.12 m3/hm2)>緩坡(637.86 m3/hm2)>斜坡(591.93 m3/hm2)>陡坡(549.62 m3/hm2)>急坡(519.52 m3/hm2)>險坡(474.61 m3/hm2),在重力作用下,土壤水分沿坡向下運動并在平緩的地方聚集[31],因此,林區單位面積水源涵養能力隨著坡度的增大而逐漸減弱。

由表5知,不同坡向(陰坡、西北、東北、東南、西南、陽坡)條件下,云杉林水源涵養量貢獻率分別為：63.20%、16.99%、12.29%、1.66%、4.55%、1.30%,這與云杉林在不同坡向分布面積比有很大關系。越接近陰坡,森林分布狀況越好,降水較陽坡多,且每次降水后,陽坡蒸散發強烈,經一段時間后,陽坡土壤含水量較陰坡要小[32],因此,陽坡水源涵養總量小于陰坡；就單位面積水源涵養量來說,越接近陰坡,其水源涵養能力越強,陰坡每公頃水源涵養量達567.90 m3/hm2,陽坡為399.94 m3/hm2,因陰坡有森林分布,土層厚度大,森林可有效增加土壤孔隙度、減小容重[33],陽坡無森林分布,土壤較薄、土層緊實,蓄積水分的能力較陰坡弱,因此,單位面積水源涵養量較陰坡低。

3.3 不同云杉林面積設定情景下水源涵養特征分析

3.3.1 研究區水源涵養變化

根據2.2.4設定的不同情景,分別統計各情景下水源涵養總量及平均水源涵養深度,得到表6。

表6 不同情景的水源涵養量

由表6知,不同情景水源涵養量順序為：情景3>情景2>情景1>情景4>情景5,隨著云杉林覆蓋率的增大(0—197.06 km2),研究區平均水源涵養深度由17.95 mm增加到27.10 mm,水源涵養總量由1.98×107m3增加到2.99×107m3。隨著研究區云杉林覆蓋率增大,其水源涵養功能增大。

以情景1云杉林現實覆蓋率及水源涵養量為基準,采用線性趨勢分析法,分析不同情景與情景1的研究區水源涵養總量變化差異,得到圖6及如下方程：

ΔWR= 0.0437ΔF- 0.0505R2= 0.9951P<0.01

式中,ΔF為云杉林面積變化率；ΔWR為研究區水源涵養總量變化。計算表明云杉林面積變化,對研究區水源涵養總量變化影響極顯著,云杉林覆蓋率每增加1%,研究區水源涵養總量增加0.437×106m3。

圖6 云杉林變化率與研究區水源涵養關系圖 Fig.6 Relationship between the change rate of spruce forest and water conservation in the study area

圖7 云杉林變化率與林區水源涵養關系圖 Fig.7 Relationship between the change rate of spruce forest and water conservation in forest area

3.3.2 云杉林面積變化及其水源涵養量的響應特征

運用ArcGIS空間統計分析功能,統計不同情景下林區水源涵養量,見表7。由表7知,不同情景林區水源涵養功能大小順序為：情景3>情景2>情景1>情景4>情景5；在未來林地面積增加的兩種情景下(情景2、情景3),林區平均水源涵養深度從54.25 mm分別上升到59.74、62.63 mm,林區水源涵養總量從4.93×106m3分別上升到8.13×106、12.28×106m3,其單位面積涵養量從542.65 m3/hm2分別上升到596.74、623.16 m3/hm2；當未來林地面積減少(情景4)直至消失(情景5)時,林區平均水源涵養深度從54.25 mm下降到53.59 mm,水源涵養總量減少了1.52×106m3。

表7 不同情景的林區水源涵養量

-表示林地消失,不計算水源涵養量

以情景1云杉林現實覆蓋率及水源涵養量為基準,采用線性趨勢分析法,分析不同情景與情景1的林區水源涵養量變化差異,得到圖7及如下方程：

ΔWRf= 0.5215ΔF- 0.2852R2= 0.9995P<0.05

式中,ΔF為云杉林面積變化率；ΔWRf林區水源涵養量變化。計算表明云杉林面積變化,對林區水源涵養量有顯著的影響,云杉林覆蓋率每增加1%,林區水源涵養量增加0.522×106m3。

綜上所述,隨著云杉林面積的不斷增加,林區水源涵養功能不斷增強,平均水源涵養深度、水源涵養總量及單位面積涵養量均呈現增加趨勢,云杉林具有顯著的水源涵養能力。

4 討論

(1)有關學者研究表明：常綠針葉林最大水源涵養深度為321.94 mm,平均水源涵養深度為83.29 mm[17],這與本文研究結果有一定差異,可能與森林類型、局地氣候、局地土壤條件等因子不同有關。本文運用InVEST模型對流域尺度云杉天然林水源涵養功能進行了定量評估,云杉林區水源涵養功能呈現出北高南低、西高東低的趨勢,這是否與云杉林的生長分布及其微環境有關？還有待進一步研究。此外,由于森林水源涵養功能的復雜性,未來還需對以下問題進行研究：①如何充分利用現有林分尺度、坡面尺度的觀測研究成果,實現森林水源涵養功能的尺度推移？②在使用InVEST模型過程中,如何解決同一土壤類型的物理屬性(土壤飽和導水率、可利用有效水)在空間分布上存在的差異？怎樣設計更大密度的野外觀測及調查方案,才能有效模擬和彌補這一不足？

(2)縱觀各學者研究成果,對不同土地利用類型的水源涵養能力觀點基本一致,森林、草地等植被覆蓋率較高的區域水源涵養能力較強,裸地、建設用地等硬質地面及植被稀疏區水源涵養能力較低[15,30],但流域水體水源涵養量最低,這與本研究結果一致。其原因是河道兩側的河床土壤粒徑較大,河道水分消耗于垂直方向的滲透與水平方向的側向流,水分易流失,且在冬季河道兩側水分向河道匯集,而在夏季則是河流水分向河道兩側流失,河流作為水分傳輸的載體,對水分的截留能力相對較弱。

(3)吳丹等[3]、劉璐璐等[29]、孫清琳等[34]分別對江西省、四川省、浙江省的森林水源涵養與地形因子的關系進行了研究,結果均表明水源涵養量的貢獻率與森林面積分布有密切關系,每公頃水源涵養量隨坡度增大而降低；這與本文研究結果一致。但對每公頃水源涵養量與海拔的關系卻不盡相同,這可能是不同區域、不同海拔森林分布的狀況不同所致。

5 結論

(1)InVEST模型可較好的確定烏魯木齊河流域的水源涵養量。InVEST模型模擬表明,研究區水源涵養總量2.41×107m3,其中,占研究區總面積8.17%的云杉林區,水源涵養總量為4.93×106m3,達研究區水源涵養總量的20.46%,林區平均水源涵養深度為54.25 mm。

(2)云杉林區水源涵養量的貢獻率與森林面積分布密切相關,其水源涵養量貢獻率在海拔、坡度、坡向上的分布與云杉林空間分布格局一致；云杉林每公頃水源涵養量隨海拔升高先增大后降低、隨坡度增大而降低,越接近陰坡,水源涵養能力越強。

(3)研究區各地類中以林地水源涵養能力最強(每公頃水源涵養量為544.78 m3/hm2),隨著云杉林覆蓋率的不斷增加,水源涵養功能不斷增強,平均水源涵養深度、水源涵養總量及單位面積涵養量均呈增加趨勢。

(4)隨著森林面積增加,研究區水源涵養總量增加。云杉林面積每增加1%,研究區與云杉林區水源涵養總量分別增加0.437×106m3、0.522×106m3。因此,對研究區進行有計劃的撫育造林,逐步擴大森林面積,可有效增加森林水源涵養生態服務功能。

[1] 張志強, 余新曉, 趙玉濤, 秦永勝. 森林對水文過程影響研究進展. 應用生態學報, 2003, 14(1): 113- 116.

[2] 張彪, 李文華, 謝高地, 肖玉. 北京市森林生態系統的水源涵養功能. 生態學報, 2008, 28(11): 5619- 5624.

[3] 吳丹, 邵全琴, 劉紀遠. 江西泰和縣森林生態系統水源涵養功能評估. 地理科學進展, 2012, 31(3): 330- 336.

[4] 莫菲, 李敘勇, 賀淑霞, 王曉學. 東靈山林區不同森林植被水源涵養功能評價. 生態學報, 2011, 31(17): 5009- 5016.

[5] 賀淑霞, 李敘勇, 莫菲, 周彬, 高廣磊. 中國東部森林樣帶典型森林水源涵養功能. 生態學報, 2011, 31(12): 3285- 3295.

[6] 趙傳燕, 別強, 彭煥華. 祁連山北坡青海云杉林生境特征分析. 地理學報, 2010, 65(1): 113- 121.

[7] 白楊, 鄭華, 莊長偉, 歐陽志云, 徐衛華. 白洋淀流域生態系統服務評估及其調控. 生態學報, 2013, 33(3): 711- 717.

[8] 馬良, 金陶陶, 文一惠, 吳秀芹, 劉桂環. InVEST模型研究進展. 生態經濟, 2015, 31(10): 126- 131, 179- 179.

[9] Nelson E, Mendoza G, Regetz J, Polasky S, Tallis H, Cameron D, Chan K M, Daily G C, Goldstein J, Kareiva P M, Lonsdorf E, Naidoo R, Ricketts T H, Shaw M. Modeling multiple ecosystem services, biodiversity conservation, commodity production, and tradeoffs at landscape scales. Frontiers in Ecology and the Environment, 2009, 7(1): 4- 11.

[10] Fisher B, Turner R K, Burgess N D, Swetnam R D, Green J, Green R E, Kajembe G, Kulindwa K, Lewis S L, Marchant R, Marshall A R, Madoffe S, Munishi P K T, Morse-Jones S, Mwakalila S, Paavola J, Naidoo R, Ricketts T, Rouget M, Willcock S, White S, Balmford A. Measuring, modeling and mapping ecosystem services in the Eastern Arc Mountains of Tanzania. Progress in Physical Geography, 2011, 35(5): 595- 611.

[11] Shaw M R, Pendleton L, Cameron D R, Morris B, Bachelet D, Klausmeyer K, Mackenzie J, Conklin D R, Bratman G N, Lenihan J, Haunreiter E, Daly C, Roehrdanz P R. The impact of climate change on California′s ecosystem services. Climatic Change, 2011, 109(Supplement 1): 465- 484.

[12] Goldstein J H, Caldarone G, Duarte T K, Ennaanay D, Hannahs N, Mendoza G, Polasky S, Wolny S, Daily G C. Integrating ecosystem-service tradeoffs into land-use decisions. Proceedings of the National Academy of Sciences of the United States of America, 2012, 109(19): 7565- 7570.

[13] 周彬, 余新曉, 陳麗華, 張振明, 呂錫芝, 范敏銳. 基于InVEST模型的北京山區土壤侵蝕模擬. 水土保持研究, 2010, 17(6): 9- 13, 19- 19.

[14] 包玉斌, 劉康, 李婷, 胡勝. 基于InVEST模型的土地利用變化對生境的影響——以陜西省黃河濕地自然保護區為例. 干旱區研究, 2015, 32(3): 622- 629.

[15] 傅斌, 徐佩, 王玉寬, 彭怡, 任靜. 都江堰市水源涵養功能空間格局. 生態學報, 2013, 33(3): 789- 797.

[16] 余新曉, 周彬, 呂錫芝, 楊之歌. 基于InVEST模型的北京山區森林水源涵養功能評估. 林業科學, 2012, 48(10): 1- 5.

[17] 王紀偉, 劉康, 甕耐義. 基于In-VEST模型的漢江上游森林生態系統水源涵養服務功能研究. 水土保持通報, 2014, 34(5): 213- 217.

[18] 李明陽, 吳軍, 時宇, 余超. 武夷山國家級自然保護區森林土壤涵水量估測及空間分析. 西南林業大學學報, 2014, 34(1): 1- 7.

[19] 姚海燕, 張民. 烏魯木齊河流域水管理信息系統. 水文, 2006, 26(2): 78- 80, 96- 96.

[20] 張東啟, 周尚哲. 天山烏魯木齊河源1號冰川對降水影響的定量分析. 冰川凍土, 2000, 22(3): 243- 249.

[21] 張寅生, 康爾泗, 劉朝海. 天山烏魯木齊河流域山區氣候特征分析. 冰川凍土, 1994, 16(4): 333- 341.

[22] 王紅娟. 基于GIS和RS的分布式融雪徑流模型在干旱區的應用研究——以烏魯木齊河為例[D]. 烏魯木齊: 新疆大學, 2004: 42- 43.

[23] 新疆維吾爾自治區農業廳. 新疆土種志. 烏魯木齊: 新疆科技衛生出版社, 1993: 6- 28.

[24] 新疆維吾爾自治區農業廳, 土壤普查辦公室. 新疆土壤. 北京: 科學出版社, 1996: 75- 116.

[25] Feinstein A R, Cicchetti D V. High agreement but low Kappa: I. the problems of two paradoxes. Journal of Clinical Epidemiology, 1990, 43(6): 543- 549.

[26] 張毓濤, 李吉玫, 常順利, 李翔, 蘆建江. 天山中部天山云杉種群空間分布格局及其與地形因子的關系. 應用生態學報, 2011, 22(11): 2799- 2806.

[27] 陸平, 嚴賡雪. 新疆森林. 烏魯木齊: 新疆人民出版社, 1989: 46- 507.

[28] 劉貴峰, 臧潤國, 郭仲軍, 巴哈爾古麗·阿尤甫, 張新平, 成克武, 白志強. 不同經度天山云杉群落物種豐富度隨海拔梯度變化. 應用生態學報, 2008, 19(7): 1407- 1413.

[29] 劉璐璐, 邵全琴, 劉紀遠, 楊存建. 瓊江河流域森林生態系統水源涵養能力估算. 生態環境學報, 2013, 22(3): 451- 457.

[30] 張海波, 張明陽, 王克林, 秦建新, 符靜. 南方丘陵山地帶水源涵養功能變化特征. 農業現代化研究, 2014, 35(3): 345- 348.

[31] 王佑民. 中國林地枯落物持水保土作用研究概況. 水土保持學報, 2000, 14(4): 108- 113.

[32] 胡振華, 王治國. 晉西黃土殘塬區坡面的日蒸散模型. 中國水土保持科學, 2003, 1(1): 95- 98.

[33] 李海軍, 張毓濤, 張新平, 韓占濤, 李吉玫. 天山中部不同林齡天然云杉林地表土壤入滲性能研究. 干旱區資源與環境, 2011, 25(5): 197- 203.

[34] 孫清琳, 李延森, 郭繼凱, 吳秀芹. 云和縣森林生態系統水源涵養功能評估. 北京大學學報: 自然科學版, 2015, 51(5): 888- 896.

Quantitative assessment of water conservation function of the natural spruce forest in the central Tianshan Mountains: a case study of the Urumqi River Basin

DING Chengfeng1,2, ZHANG Huifang1, LI Xia2, LI Weitao3, GAO Yaqi1,*

1InstituteofModernForestry,XinjiangAcademyofForestry,Urumqi830000,China2CollegeofPrataculturalandEnvironmentalScience,XinjiangAgriculturalUniversity,Urumqi830052,China3SchoolofGeographyInformationandTourism,ChuzhouUniversity,Chuzhou239000,China

To quantitatively investigate the water conservation function of the natural spruce forest at the watershed scale in the central Tianshan Mountains, the Urumqi River Basin was chosen as a study area. The InVEST model was used to investigate the capacity for water conservation of the spruce forest and other land use types, as well as quantitatively analyze water conservation changes of the spruce forest under the conditions of different distribution areas and terrain factors. The results showed that: (1) The InVEST model could adequately assess water conservation at the watershed scale. The InVEST model simulation showed that the total water conservation of the spruce forest in the Urumqi River Basin was 4.93 × 106m3, which accounted for 20.46% of the total water conservation in the study area (2.41 × 107m3). In addition, the average water conservation depth of forest area was 54.25 mm. (2) The contribution rate of water conservation in the slope, aspect, and elevation was consistent with the spatial distribution pattern of the spruce forest. At first, per hectare water conservation increased with altitude and then decreased, and it decreased with an increase in slope. The water conservation ability in the shade slope was greatest. (3) The spruce forest showed the greatest capacity, which reached 544.78 m3/hm2. The total, average depth and unit area of water conservation exhibited an increasing trend with increasing spruce forest coverage. (4) The contribution rate of water conservation in the spruce forest was closely related to forest area. Total water conservation in study area and spruce area increased to 0.437 × 106m3and 0.522 × 106m3with a 1% increase of spruce forest area, respectively. Therefore, the forest plays a vital role in maintaining sustainable and stable water conservation and ecological service functions, which necessitate effective measures to strengthen the protection and management of spruce forests.

water conservation function; InVEST model; the spruce forest; the Central Tianshan Mountains; Urumqi River Basin

新疆公益性科研院所基本科研業務經費資助項目(XMBM000001953);安徽省滁州學院培育項目(2014PY02)

2016- 04- 13; 網絡出版日期:2017- 02- 22

10.5846/stxb201604130673

*通訊作者Corresponding author.E-mail: gyq611003@163.com

丁程鋒，張繪芳，李霞，李偉濤，高亞琪.天山中部云杉天然林水源涵養功能定量評估——以烏魯木齊河流域為例.生態學報,2017,37(11):3733- 3743.

Ding C F, Zhang H F, Li X, Li W T, Gao Y Q.Quantitative assessment of water conservation function of the natural spruce forest in the central Tianshan Mountains: a case study of the Urumqi River Basin.Acta Ecologica Sinica,2017,37(11):3733- 3743.