生態(tài)系統(tǒng)服務(wù)功能是目前生態(tài)學(xué)研究熱點(diǎn)之一，通過(guò)直觀的貨幣形式來(lái)表述生態(tài)系統(tǒng)的價(jià)值，使人們對(duì)生態(tài)系統(tǒng)的重要性有了更新的認(rèn)識(shí)。土壤侵蝕是近年來(lái)全世界面臨的一個(gè)嚴(yán)峻問(wèn)題，因此，生態(tài)系統(tǒng)對(duì)保持土壤的功能越來(lái)越得到各國(guó)專(zhuān)家和學(xué)者的重視，在這方面已經(jīng)作了比較深入的研究［１－２］。

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

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

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

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

１．２研究方法

１．２．１數(shù)據(jù)來(lái)源

云冷杉數(shù)據(jù)來(lái)源于西藏林芝高山暗針葉林國(guó)家重點(diǎn)野外站臺(tái)（以下簡(jiǎn)稱(chēng)林芝站）多年觀測(cè)的數(shù)據(jù)，其他林分?jǐn)?shù)據(jù)來(lái)源于相鄰地區(qū)（主要為川西和云南）公開(kāi)發(fā)表的論文數(shù)據(jù)，相關(guān)的經(jīng)濟(jì)指標(biāo)數(shù)據(jù)來(lái)源于《西藏統(tǒng)計(jì)年鑒》。

１．２．２測(cè)量方法

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

１．２．３測(cè)量項(xiàng)目

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

２結(jié)果與分析

２．１年固土量和固土年價(jià)值

２．１．１森林年固土量

林芝站實(shí)測(cè)無(wú)林地年地表徑流量為２２．６ｍｍ，年平均含沙量為１２．１７ｍｇ／ｍＬ，即可得到無(wú)林地土壤侵蝕模數(shù)為２７．５ｔ／（ｈｍ２•年）［４－５］。工布自然保護(hù)區(qū)中森林主要分布區(qū)域坡度相近，年降水量也相似，所以無(wú)林地土壤侵蝕模數(shù)均采用此數(shù)據(jù)。根據(jù)林芝站的數(shù)據(jù)，冷杉林年地表徑流量為０．８６ｍｍ，地表徑流中泥沙含量為６．９４ｍｇ／ｍＬ，計(jì)算出冷杉林的侵蝕模數(shù)為０．６０ｔ／（ｈｍ２•年）［６－８］。云杉數(shù)據(jù)采用工布自然保護(hù)區(qū)南伊溝實(shí)測(cè)數(shù)據(jù)，該林分年地表徑流量為０．９７ｍｍ，地表徑流中泥沙含量為８．６４ｍｇ／ｍＬ，計(jì)算出云杉林的侵蝕模數(shù)＝０．８４ｔ／（ｈｍ２•年）［６－８］。其他林分沒(méi)有實(shí)測(cè)數(shù)據(jù)，據(jù)相關(guān)文獻(xiàn)規(guī)定［９－１３］，針葉林的土壤侵蝕模數(shù)是云冷杉的平均值，即７１７．４６ｋｇ／（ｈｍ２•年）；而闊葉樹(shù)為針葉樹(shù)的１．５倍，即１　０７６．１９ｋｇ／（ｈｍ２•年）。　　從表１和圖１可以看出，工布自然保護(hù)區(qū)森林年固土量為１　４２３．４萬(wàn)ｔ。其中，冷杉林年固土量最大，為８４３．６萬(wàn)ｔ，占總量的５９．２７％；其次是高山松（１４８．７萬(wàn)ｔ）和云杉林（１１８．０萬(wàn)ｔ），分別占總量的１０．４５％和８．２９％。單位面積森林年固土量為２６．８０ｔ／（ｈｍ２•年）。

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

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