生態(tài)系統(tǒng)服務功能是目前生態(tài)學研究熱點之一，通過直觀的貨幣形式來表述生態(tài)系統(tǒng)的價值，使人們對生態(tài)系統(tǒng)的重要性有了更新的認識。土壤侵蝕是近年來全世界面臨的一個嚴峻問題，因此，生態(tài)系統(tǒng)對保持土壤的功能越來越得到各國專家和學者的重視，在這方面已經作了比較深入的研究［１－２］。

西藏地處祖國西南邊陲，由于海拔高、氣候惡劣等自然因素，環(huán)境極為脆弱，水土流失在東南部表現極為明顯，這給下游地區(qū)造成較大的危害，為此，國家曾啟動了長江中上游天然林保護工程。西藏工布自然保護區(qū)位于西藏林芝地區(qū)境內，是迄今為止國家批復投資最大的林業(yè)自然保護區(qū)建設工程［３］，它不僅保護了獨特的原始森林和自然景觀，而且也保護了該區(qū)域的環(huán)境。就保持土壤功能而言，目前在該區(qū)域所做工作較少，筆者采用國家林業(yè)局頒布的行業(yè)標準《森林生態(tài)系統(tǒng)服務功能評估規(guī)范》方法，利用部分實測數據和相似區(qū)域的公開發(fā)表數據，對工布自然保護區(qū)生態(tài)系統(tǒng)土壤保持功能及其經濟價值進行研究，有助于深入認識工布保護區(qū)保護土壤的重要性，加強對生態(tài)系統(tǒng)的保護，保證其生態(tài)服務功能的正常發(fā)揮。

１研究區(qū)域與方法

１．１研究區(qū)概況

西藏工布自然保護區(qū)位于西藏自治區(qū)林芝地區(qū)境內，地理范圍為東經９３°１３′～９４°５０′，北緯２８°４０′～３０°２０′，處于南迦巴瓦峰西、喜馬拉雅山東段與念青唐古拉山脈之間，東起雅魯藏布大峽谷國家級自然保護區(qū)西界，西至工布江達縣城，南接隆子縣、墨脫縣縣界，北靠嘉黎縣縣界［３］，總面積約２１　５５８．１６ｋｍ２，涉及到林芝地區(qū)的工布江達、米林、朗縣、林芝共４個縣１９個鄉(xiāng)鎮(zhèn)。從植被類型來看，工布自然保護區(qū)包含三類：東部色季拉山為典型山地森林垂直帶譜；北部尼洋河流域受溫度的影響，形成從半濕潤森林到半干旱高山草甸的各種植被類型；南部雅魯藏布江流域受水分的作用，形成從濕潤森林到荒漠的各種植被類型［３］。從經濟社會來看，林芝地區(qū)行署所在地八一鎮(zhèn)處于本自然保護區(qū)中，是林芝地區(qū)政治、經濟、文化、交通中心，人類活動相對頻繁，同時，保護區(qū)中的尼洋河和雅魯藏布江河谷地帶，也是城鎮(zhèn)、村莊、農業(yè)等土地利用較為集中地帶，是今后林芝地區(qū)社會經濟發(fā)展的主要區(qū)域［３］。

１．２研究方法

１．２．１數據來源

云冷杉數據來源于西藏林芝高山暗針葉林國家重點野外站臺（以下簡稱林芝站）多年觀測的數據，其他林分數據來源于相鄰地區(qū)（主要為川西和云南）公開發(fā)表的論文數據，相關的經濟指標數據來源于《西藏統(tǒng)計年鑒》。

１．２．２測量方法

本文計算方法均采用２００８年國家林業(yè)局發(fā)布的行業(yè)標準《森林生態(tài)系統(tǒng)服務功能評估規(guī)范》中所確定的評估公式。

１．２．３測量項目

１）年固土量Ｇ固土＝Ａ（Ｘ２?Ｘ１）式中，Ｇ固土為森林年固土量（ｔ／年），Ａ為森林面積（ｈｍ２），Ｘ２為無林地土壤侵蝕模數［ｔ／（ｈｍ２•年）］，Ｘ１為林地土壤侵蝕模數［ｔ／（ｈｍ２•年）］。２）固土年價值Ｕ固土＝Ｃ土•Ｇ固土／Ｐ式中，Ｕ固土為固土年價值（元／年），Ｃ土為挖取和運輸單位體積土方所需的費用（元／ｍ３）（參考西藏工程預計，本文取值５．７１４元／ｍ３），Ｇ固土為年固土量（ｔ／年），Ｐ為土壤容重（ｔ／ｍ３）（本文取值為１．２ｔ／ｍ３）。３）年保肥量ＧＮ＝ＡＮ（Ｘ２?Ｘ１）＝ＮＧ固土ＧＰ＝ＡＰ（Ｘ２?Ｘ１）＝ＰＧ固土ＧＫ＝ＡＫ（Ｘ２?Ｘ１）＝ＫＧ固土ＧＭ＝ＡＭ（Ｘ２?Ｘ１）＝ＭＧ固土式中，ＧＮ、ＧＰ、ＧＫ和ＧＭ分別為每年氮、磷、鉀、有機質的流失量（ｔ／年），Ｎ、Ｐ、Ｋ和Ｍ分別為土壤中氮、磷、鉀、有機質的含量（％）。４）保肥年價值Ｕ肥＝ＧＮＣ１／Ｒ１＋ＧＰＣ１／Ｒ２＋ＧＫＣ２／Ｒ３＋ＧＭＣ３式中，Ｕ肥為保肥年價值（元／年），Ｒ１為磷酸二銨化肥中含氮量（１４％），Ｒ２為磷酸二銨化肥中含磷量（１５．０１％），Ｒ３為氯化鉀化肥中含鉀量（５０％），Ｃ１、Ｃ２、Ｃ３分別為磷酸二銨、氯化鉀、有機質的市場價格（元／ｔ）。

２結果與分析

２．１年固土量和固土年價值

２．１．１森林年固土量

林芝站實測無林地年地表徑流量為２２．６ｍｍ，年平均含沙量為１２．１７ｍｇ／ｍＬ，即可得到無林地土壤侵蝕模數為２７．５ｔ／（ｈｍ２•年）［４－５］。工布自然保護區(qū)中森林主要分布區(qū)域坡度相近，年降水量也相似，所以無林地土壤侵蝕模數均采用此數據。根據林芝站的數據，冷杉林年地表徑流量為０．８６ｍｍ，地表徑流中泥沙含量為６．９４ｍｇ／ｍＬ，計算出冷杉林的侵蝕模數為０．６０ｔ／（ｈｍ２•年）［６－８］。云杉數據采用工布自然保護區(qū)南伊溝實測數據，該林分年地表徑流量為０．９７ｍｍ，地表徑流中泥沙含量為８．６４ｍｇ／ｍＬ，計算出云杉林的侵蝕模數＝０．８４ｔ／（ｈｍ２•年）［６－８］。其他林分沒有實測數據，據相關文獻規(guī)定［９－１３］，針葉林的土壤侵蝕模數是云冷杉的平均值，即７１７．４６ｋｇ／（ｈｍ２•年）；而闊葉樹為針葉樹的１．５倍，即１?。埃罚叮保梗耄纾ǎ瑁恚?bull;年）?！　谋恚焙蛨D１可以看出，工布自然保護區(qū)森林年固土量為１?。矗玻常慈fｔ。其中，冷杉林年固土量最大，為８４３．６萬ｔ，占總量的５９．２７％；其次是高山松（１４８．７萬ｔ）和云杉林（１１８．０萬ｔ），分別占總量的１０．４５％和８．２９％。單位面積森林年固土量為２６．８０ｔ／（ｈｍ２•年）。

２．１．２非森林系統(tǒng)對于其他生態(tài)系統(tǒng)年固土量

仍然按照公式進行計算，其中，各生態(tài)系統(tǒng)的侵蝕模數數據來源于相關文獻［１４－１６］，從計算結果（表２）可以看出，工布自然保護區(qū)除森林系統(tǒng)之外，其他生態(tài)系統(tǒng)年固土總量為５０５．９２萬ｔ。其中，灌木林固土量最大２２７．４０萬ｔ，占總量的４４．９５％；其次是草地，三類草地共占總量的３３．９６％（圖２）。從表２還可以看出，不同的生態(tài)系統(tǒng)單位面積年固土量差異較大。其中，固土能力最大的是其他類型的土地，１ｈｍ２每年固土量高達２８．４８ｔ，這里的其他類型土地主要是指農村居民點，它的潛在侵蝕模數較大，但其上面多有建筑物，所以固土能力最強；其次是水域，由于河流的切割作用，一般水域周邊坡度較大，所以潛在侵蝕模數大，但由于被水覆蓋，所以幾乎不可能產生侵蝕現象；固土能力最弱的是無林地和未成林地。