三峽庫區(qū)典型馬尾松林生態(tài)系統(tǒng)碳分配格局研究

薛沛沛，李 彬，王軼浩，馮大蘭，朱小龍*，盧 斌

(1.重慶市林業(yè)科學(xué)研究院，三峽庫區(qū)森林保護與恢復(fù)重慶市市級重點實驗室，重慶 400036;2.北京林業(yè)大學(xué)，北京 100083)

21世紀(jì)，全球變化問題已成為聯(lián)合國氣候變化峰會和APEC等世界各國首腦會議的主題。全球變化引起了研究者對陸地生態(tài)系統(tǒng)中碳平衡、碳儲量的關(guān)注。森林作為陸地生態(tài)系統(tǒng)的主體，也是地球上最大的碳庫，貯存了陸地生態(tài)系統(tǒng)總碳庫的56%［1，2］，森林生態(tài)系統(tǒng)地上部分和地下部分碳庫分別占陸地生態(tài)系統(tǒng)碳庫的80%和40%［3］。

隨著全球變暖趨勢的日益明顯，森林固碳問題越來越受到各國學(xué)者的重視［4］，并對此做了大量的研究［5－9］，關(guān)于馬尾松林分碳儲量和碳分配格局的報道也不少［10，11］。然而，三峽庫區(qū)作為馬尾松最主要的適生區(qū)之一，對千萬畝馬尾松林的碳分配格局研究報道很少。為此，本文運用野外實地采樣和室內(nèi)分析相結(jié)合的方法，并結(jié)合文獻資料，對三峽庫區(qū)馬尾松林碳儲量及分配格局進行了探討，以期為三峽庫區(qū)馬尾松林碳分配的進一步研究和正確評價馬尾松林的生態(tài)效益提供科學(xué)依據(jù)。

1 材料與方法

本研究采取基于文獻的研究數(shù)據(jù)資料和實測方法相結(jié)合，對三峽庫區(qū)馬尾松林的碳分配進行研究。

1.1 文獻資料

張治軍等對重慶鐵山坪馬尾松林生態(tài)系統(tǒng)碳貯量及其分配特征的研究［12］。

1.2 實測方法

1.2.1 研究區(qū)概況

研究區(qū)位于重慶市萬州區(qū)新田林場，長江以南，方斗山脈以北，東臨武陵鎮(zhèn)、高峰鎮(zhèn)，西接石柱縣，南與百安辦接壤，交通較為方便。地處北緯30°36'19.8″，東經(jīng) 108°23'43.3″。海拔1 080.3 m。立地條件較好，黃壤，中性至微酸性，中山。全場林業(yè)用地面積3 262.3 hm2，活立木總蓄積為43.5萬m3。森林植被以馬尾松、杉木為主，灌木有青岡、巨木、柃木、山礬、青冬、杜鵑等，蓋度為20%，地被物有茅草、蕨類，蓋度為25%。

1.2.2 標(biāo)準(zhǔn)地選擇

2011年1月～2月，在40 a～45 a生的馬尾松林中，設(shè)置4個具有代表性的面積為20 m×20 m的標(biāo)準(zhǔn)地。標(biāo)準(zhǔn)地基本情況見表1。

表1 新田林場馬尾松林喬木群落的一些參數(shù)Table1 Parameters of arbor layer of Pinus massoniana forest in Xintian State Forest Farm

1.2.3 樣品采集

1.2.3.1 喬木層生物量調(diào)查 選擇標(biāo)準(zhǔn)木進行樹干解析。(1)樹干調(diào)查:樹干解析方法參照測樹學(xué)［13］，并現(xiàn)場稱鮮質(zhì)量。在樹木地徑處以及全高的0.25、0.5和0.75處分別截取3 cm～5 cm厚的圓盤，裝入自封袋做好標(biāo)記，帶回實驗室。(2)枝葉調(diào)查:以枝基徑 ＜1 cm、1 cm ～2 cm、2 cm ～4 cm、＞4 cm為標(biāo)準(zhǔn)進行分級，稱取鮮質(zhì)量;每級選取3個標(biāo)準(zhǔn)枝稱取帶葉枝總鮮質(zhì)量;摘凈葉，分別稱取枝葉鮮質(zhì)量;按等級分別裝入自封袋做好標(biāo)記，帶回實驗室。(3)根系調(diào)查:分不同方向(樹干基部的坡上、坡下、左、右)、層次(0 ～20 cm、20 cm ～40 cm、40 cm～60 cm、60 cm ～80 cm、80 cm ～100 cm 等層次，直至無根系分布)挖取樹木根系，用水沖洗，風(fēng)干后分d＜1 cm、1 cm ～2 cm、2 cm ～5 cm、5 cm ～10 cm、＞10 cm 5類，分層分類稱鮮質(zhì)量后，隨機抽取1 kg樣品，裝入自封袋做好標(biāo)記，帶回實驗室。

1.2.3.2 林下植被生物量調(diào)查 在每個標(biāo)準(zhǔn)地內(nèi)，按對角線法設(shè)3 m×3 m的樣方4個，采用“樣方收獲法”［14，15］測定灌木層生物量。由于林地地表幾乎沒有草本層，因此本研究在碳儲量測算中草本層忽略不計。按對角線法設(shè)置1 m×1 m的樣方4個，凋落物按未分解、半分解測定其鮮質(zhì)量。

1.2.3.3 土壤有機碳調(diào)查 每一標(biāo)準(zhǔn)地內(nèi)按“S”形布點方法，挖取3個土壤剖面，分0～20 cm、20 cm～40 cm、40 cm ～60 cm、60 cm ～100 cm 等4個層次取樣，每層取樣不少于1 kg。同時用100 cm3環(huán)刀測定土壤容重。土樣風(fēng)干，去除礫石等雜物后過2 mm篩備用。

1.2.4 數(shù)據(jù)分析與處理

1.2.4.1 植物樣品 帶回的植物樣品，分別在80℃下烘干至恒重，稱重。根據(jù)模型W=a(D2H)b和每木檢尺數(shù)據(jù)，求算林分喬木樹干、枝、葉和根生物量。按照方運霆等對廣東鼎湖山馬尾松純林(林齡約為50 a以上)各器官與林下植被及凋落物碳含量的測定結(jié)果，對本地區(qū)植被的碳貯量進行估算［16］，喬木層各器官的碳含量分別為:樹干61.38%，樹枝53.33%，樹葉44.63%，樹根59.13%。木質(zhì)物殘體層以及林下植被層碳儲量則采用常見的0.45的碳轉(zhuǎn)化率［17］。然后根據(jù)標(biāo)準(zhǔn)地密度計算每公頃馬尾松林的碳儲量。

1.2.4.2 土壤樣品 以重鉻酸鉀法—外加熱法測定有機碳含量［18］;再根據(jù)土壤有機碳含量、土壤密度和土層深度計算土壤有機碳儲量。

土壤有機碳儲量計算公式為:

式中:S表示土壤有機碳儲量，i表示土壤層次，Ci表示土壤有機碳含量，di表示土壤密度，Di表示土壤深度。本研究土壤密度為1.28 g·cm－3～1.55 g·cm－3，因此，本研究土壤各層次的有機碳貯量用密度的平均值1.41 g·cm－3進行換算。

2 結(jié)果與分析

2.1 馬尾松林生物量及其分配特征

生物量是植物群落最重要的數(shù)量特征之一，直接反映了生態(tài)系統(tǒng)生產(chǎn)者的物質(zhì)生產(chǎn)量，也是生態(tài)系統(tǒng)生產(chǎn)力的重要體現(xiàn)［19］。

測定結(jié)果表明，新田林場40 a～45 a生馬尾松林的總生物量平均為145.48 t·hm－2(見表2)。資料顯示，鐵山坪林場45 a～51 a生馬尾松林總生物量平均為142.06 t·hm－2。兩個林場的馬尾松生物量表現(xiàn)出如下特征:(1)喬木層生物量占馬尾松林生態(tài)系統(tǒng)生物量的絕對優(yōu)勢(新田林場占85.68%，鐵山坪林場占89.91%)，灌木層次之(新田林場占11.42%，鐵山坪林場占5.61%)，凋落物層分別占2.90%(未分解層占64.27%，半分解層占35.73%)和2.98%(未分解層占 62.03%，半分解層占37.97%)。由于實測試驗是在年初進行的，所以未進行草本層的采集和調(diào)查，因此本文不進行討論。可見，在森林生態(tài)系統(tǒng)生長空間尺度上，生物量分布呈現(xiàn)由上到下遞減的規(guī)律［20］。

表2 新田林場馬尾松林分不同層次生物量分配情況Table 2 Biomass distribution of stand in different layers of Pinus massoniana forest in Xintian State Forest Farm

調(diào)查發(fā)現(xiàn)，新田林場馬尾松林的總生物量差異不大，標(biāo)準(zhǔn)地4的生物量最大，為153.80 t·hm－2，其它3個樣方的生物量排序依次為標(biāo)準(zhǔn)地1＞標(biāo)準(zhǔn)地2＞標(biāo)準(zhǔn)地3。但鐵山坪林場各標(biāo)準(zhǔn)地的總生物量相差較大，生物量最大樣方1(161.11 t·hm－2)是生物量最小樣方(115.73 t·hm－2)的1.39倍，這可能與樹木密度、樹木分布或立地條件等有關(guān)。

喬木層是人工生態(tài)系統(tǒng)中最大的物質(zhì)生產(chǎn)者，它的生物量特點直接反映了林分有機物質(zhì)的生產(chǎn)與積累水平［21］。喬木層生物量中，樹干生物量占整株樹木生物總量的60%以上，馬尾松各器官生物量占各自總生物量的比例依次為干＞枝＞根＞葉(見表3)。由此可見，在喬木層生物量中，樹干的比重最高，這與張尚炬［22］對沿海馬尾松臺灣相思次生林生態(tài)系統(tǒng)生物量及碳貯量的研究結(jié)果一致。

表3 新田林場馬尾松林喬木層的生物量及其器官分配Table 3 The biomass and organ distribution of arbor layer in Pinus massoniana forest in Xintian State Forest Farm

2.2 馬尾松林各組分碳含量及其分配特征

新田林場馬尾松林各組分碳含量見表4。表4顯示，馬尾松林的碳分配以喬木層占的比例最大，所占比例88.79%;灌木層碳儲量為7.37 t·hm－2，凋落物層的碳儲量所占比例為2.28%。碳儲量在三峽庫區(qū)馬尾松林的分配同生物量的分配一致，表現(xiàn)出從上到下依次遞減的趨勢。

表4 新田林場馬尾松林林各組分碳含量Table4 Carbon storage of of different parts in Pinus massoniana forest in Xintian State Forest Farm

根據(jù)馬尾松的生物量和碳含量，得出各器官的碳儲量。由表5可以看出，新田林場馬尾松林喬木層的碳儲量均值為73.29 t·hm－2，比鄧華平研究的35 a生的馬尾松林喬木層碳儲量的115.52 t·hm－2減少了42.33 t·hm－2。資料表明，鐵山坪馬尾松林喬木層的碳儲量均值為76.06 t·hm－2。40 a～45 a生的馬尾松林碳儲量與45 a～51 a生馬尾松林相差2.77 t·hm－2。由此說明馬尾松林在45 a之后生長比較緩慢。同時也可以看出三峽庫區(qū)的馬尾松林較其它地區(qū)生長較慢，生物量積累較慢，這可能與未進行人工經(jīng)營有關(guān)，從而導(dǎo)致碳儲量偏低。這與夏鑫［23］對馬尾松老齡林的研究具有高度的一致性。馬尾松林喬木層各器官碳儲量的排序為:樹干＞樹枝＞樹根＞樹葉。

表5 新田林場馬尾松林喬木層各器官的碳儲量Table 5 Carbon storage of different organ of arbor layer in Pinus massoniana forest in Xintian State Forest Farm

2.3 不同層次土壤有機碳儲量及其分配特征

馬尾松林各樣方土壤碳儲量隨著土層深度的增加而逐漸減小。由表6可看出，土壤碳儲量以表層土(0 ～20 cm)最大，均值為 68.61 t·hm－2，占土壤層總碳儲量的55.45%。20 cm～40 cm土壤碳儲量25.61 t·hm－2，占總碳儲量的20.69%;40 cm ～60 cm土壤碳儲量15.62 t·hm－2，占總碳儲量的12.63%;60 cm ～100 cm 土壤碳儲量 13.90 t·hm－2，占總碳儲量的11.23%。可見，0～40 cm土層，在森林有機碳儲量的研究中，占有舉足輕重的地位。鐵山坪0～62 cm土壤層總有機碳貯量為114.96 t·hm－2，腐殖質(zhì)層有機碳含量最高，平均為9.68%，而2 cm～12 cm土層有機碳含量平均為2.17%，42 cm～62 cm土層有機碳含量僅為0.76%?？梢?，三峽庫區(qū)的土壤有機碳貯量基本上是由表層向深層逐漸降低。

表6 新田林場不同層次土壤有機碳儲量Table 6 Carbon storage of different soil layers in Xintian State Forest Farm

2.4 馬尾松林生態(tài)系統(tǒng)碳儲量的空間分布

新田林場馬尾松林生態(tài)系統(tǒng)總碳儲量平均為206.28 t·hm－2，其中，植被層占 40.01%，土壤層占59.99%。植被層中，喬木層占35.53%，灌木層占3.57%，凋落物層占1.05%。鐵山坪林場馬尾松林生態(tài)系統(tǒng)的總有機碳貯量為197.78 t·hm－2，其中土壤有機碳貯量占58.13%?？梢?，森林土壤是森林生態(tài)系統(tǒng)以及陸地生態(tài)系統(tǒng)最主要的碳庫之一。此結(jié)果與韓娟娟［24］對子午嶺遼東櫟類群落碳儲量的研究結(jié)果相似。

因此，研究表明，三峽庫區(qū)馬尾松林碳儲量的分布以喬木層和土壤碳庫為主，但與其它研究相比，三峽庫區(qū)馬尾松林的碳儲量比較低，而且40 a～51 a生的馬尾松林已達成熟期。為了增加森林的固碳能力，應(yīng)加強對森林的管理、經(jīng)營以及改造。

3 結(jié)論與討論

3.1 調(diào)查方法

目前，研究森林生態(tài)系統(tǒng)植冠層生物量和碳儲量的方法主要有樣地清查法、模型模擬法、遙感估算法等。樣地清查法適用于小尺度森林生態(tài)系統(tǒng)的研究，實測所有樹木的生物量和碳儲量頗為困難。模型模擬法主要用于大尺度森林生態(tài)系統(tǒng)碳循環(huán)研究，但由于一些生態(tài)學(xué)過程的特征參數(shù)不易獲得或難以把握，以及可靠的觀測數(shù)據(jù)的可獲得性標(biāo)準(zhǔn)不一致，模型化很難［25］。遙感估算法能滿足不同的研究在時間和空間尺度上的要求，適合于區(qū)域尺度上的研究。通過比較分析，本研究采用樣地清查法。

不同的研究人員用不同方法對森林生態(tài)系統(tǒng)的生物量和碳儲量進行研究，所得的結(jié)果也存在差異［26］。如設(shè)置樣方面積的大小、樣方選擇的代表性等。本研究所應(yīng)用的標(biāo)準(zhǔn)地20 m×20 m、灌木樣方3 m×3 m、凋落物樣方1 m×1 m，是參照大量文獻［27～29］的基礎(chǔ)上，根據(jù)實際情況而確定的。

3.2 植被碳含量

過去的區(qū)域和國家尺度的森林生態(tài)系統(tǒng)碳儲量的估算中，國內(nèi)外研究者大多采用0.5來作為所有森林類型的平均含碳率［30～32］，也有采用0.45作為平均含碳率的［33］。極少根據(jù)樹木不同組分的不同碳含量進行計算。如劉應(yīng)芳等［34］對蜀南竹海風(fēng)景區(qū)毛竹林生態(tài)系統(tǒng)碳儲量及其空間分配特征進行研究時指出，為了更準(zhǔn)確地計算毛竹林碳儲量，應(yīng)根據(jù)毛竹不同器官采用不同的含碳量作為轉(zhuǎn)換系數(shù)。前人對鼎湖山保護區(qū)馬尾松林的碳素含量進行了測定，因此，本研究中，喬木層各器官的碳含量分別按樹干 61.38%、樹枝 53.33%、樹葉 44.63、樹根59.13%進行換算。

群落的總生物量主要取決于喬木層，喬木層的生物量主要以樹干為主。三峽庫區(qū)馬尾松林分的生物量與表現(xiàn)出相同的特征，該研究結(jié)果與方運霆等［35］對馬尾松和荷木混交林、馬煒等［36］對長白落葉松、周國模等［37］對毛竹林的研究以及張駿［38］對中國中亞熱帶東部森林生態(tài)系統(tǒng)的研究結(jié)果一致。

3.3 土壤碳含量

土壤碳儲量所用公式不同，獲得的結(jié)果不盡相同。本研究所采用的土壤碳儲量計算方法為大部分土壤碳儲量研究中常用的方法。研究結(jié)果表明土壤碳儲量隨土層深度的增加而逐漸降低，這與大多數(shù)研究的結(jié)果呈現(xiàn)出相似的規(guī)律，如楊曉梅等［20］對不同林地土壤有機碳儲量及垂直分布特征的研究、李正才等［39］對北亞熱帶土地利用變化對土壤有機碳垂直分布特征及儲量的影響的研究。

3.4 展望

本研究實測的4個標(biāo)準(zhǔn)地中，在植被碳儲量上，標(biāo)準(zhǔn)地4占據(jù)相對優(yōu)勢，是由于標(biāo)準(zhǔn)地4的土壤碳含量高的原因。對植被層來說，標(biāo)準(zhǔn)地1的碳儲量最大，說明新田林場密度為1 500株·hm－2的馬尾松林表現(xiàn)相對較好，在以后的研究中，可以對產(chǎn)生該現(xiàn)象的具體原因做詳細研究。但與其它地區(qū)對馬尾松林的研究相比，三峽庫區(qū)馬尾松林生態(tài)系統(tǒng)的碳儲量相對偏低，這可能與未對其進行經(jīng)營管理有關(guān)，也可能與該研究所采用的方法有關(guān)。本研究所選擇的新田林場和鐵山坪林場，雖然具有較好的代表性，但由于三峽庫區(qū)包括范圍廣，地勢起伏，立地條件相差較大，因此，在今后的研究中，應(yīng)該選取在三峽庫區(qū)多處設(shè)置樣地，應(yīng)進行不同取樣方法、不同數(shù)據(jù)處理和計算方法、不同森林經(jīng)營方法等方面的比較和研究。

［1］ 吳慶標(biāo)，王效科，段曉男，等.中國森林生態(tài)系統(tǒng)植被固碳現(xiàn)狀和潛力［J］.生態(tài)學(xué)報，2008，28(2):517 ～524.

［2］ Houghton R A.Aboveground forest biomass and the global carbon balance［J］.Global Change Biology，2005，11:945 ～958.

［3］ Dixon R K.Carbon pools and flux of global forest ecosystems［J］.Science，1994，263:185 ～190.

［4］ Matamala R，M A.Gonzalez-Meler，Jastrow J D，et al.Impacts of fine root turnover on forest NPP and soil C sequestration potential［J］.Science，2003，302:1385 ～1387.

［5］ Christine M.Johnson，Ima C G，Vieira，Daniel J.Zarin，et al.Carbon and nutrient storage in primary and secondary forests in eastern Amazonia［J］.Forest Ecology and Management，2001，147:245 ～252.

［6］ Catharina J E.Schulp，Gert-Jan Nabuurs，Peter H.Verburg，et al.Effect of tree species on carbon stocks in forest floor and mineral soil and implications for soil carbon inventories［J］.Forest Ecology and Management，2008，256:482 ～490.

［7］ 程堂仁，馮菁，馬欽彥，等.甘肅小隴山森林植被碳庫及其分配特征［J］.生態(tài)學(xué)報，2008，28(1):33 ～44.

［8］ 張萍.北京森林碳儲量研究［D］.北京:北京林業(yè)大學(xué)，2009.

［9］ 趙敏，周廣勝.中國森林生態(tài)系統(tǒng)的植物碳貯量及其影響因子分析［J］.地理科學(xué)，2004，24(1):50 ～54.

［10］ 鄧華平，耿賡，王正超.豫南35年生馬尾松林生態(tài)系統(tǒng)碳庫特征及其分配［J］.中南林業(yè)科技大學(xué)學(xué)報，2010，30(6):5～9.

［11］ 張國慶，黃從德，郭恒，等.不同密度馬尾松人工林生態(tài)系統(tǒng)碳儲量空間分布格局［J］.浙江林業(yè)科技，2007，27(6):10～14.

［12］ 張治軍，張小全，王彥輝，等.重慶鐵山坪馬尾松林生態(tài)系統(tǒng)碳貯量及其分配特征［J］.林業(yè)科學(xué)，2009，45(5):49 ～53.

［13］ 孟憲宇.測樹學(xué)［M］.北京:中國林業(yè)出版社，2006.

［14］ 丁貴杰，王鵬臣，嚴(yán)仁發(fā).馬尾松紙漿商品用材林生物量變化規(guī)律和模型研究［J］.林業(yè)科學(xué)，1998，34(1):33 ～41.

［15］ 馮宗煒，陳楚瑩，張家武，等.湖南會同地區(qū)馬尾松林生物量的測定［J］.林業(yè)科學(xué)，1982，18(2):127 ～134.

［16］ 方運霆，莫江明.鼎湖山馬尾松林生態(tài)系統(tǒng)碳素分配和貯量的研究［J］.廣西植物，2002，22(4):305 ～310.

［17］ IPCC.IPCC Guidelines for National Greenhouse Gas Inventories，Prepared by the National Greenhouse Gas Inventories Program∥Eggleston H.S.，Buendia L.，Miwa K.，Ngara T.and Tanabe K.eds.Hayama:IPCC/IGES，Japan，2006.

［18］ 中國科學(xué)院南京土壤研究所.土壤理化分析［M］.上海:上海科學(xué)技術(shù)出版社，1978.

［19］ 張峰，上官鐵梁，李素珍.關(guān)于灌木生物量建模方法的改進［J］.生態(tài)學(xué)雜志，1993，12(6):67 ～69.

［20］ 楊曉梅，程積民，孟蕾.黃土高原天然柴松林碳儲量與碳密度特征［J］.中國水土保持科學(xué)，2010，8(2):41 ～45.

［21］ 朱宇林，溫遠光，曹福亮，等.短周期尾巨按連栽林分生產(chǎn)力的研究［J］.西北農(nóng)業(yè)大學(xué)學(xué)報，2006，28(1):90 ～94.

［22］ 張尚炬.沿海馬尾松臺灣相思次生林生態(tài)系統(tǒng)生物量及碳貯量的研究［D］.福建:福建農(nóng)林大學(xué)，2008.

［23］ 夏鑫.馬尾松人工林的老齡林生態(tài)系統(tǒng)生物量和碳貯量研究［D］.福建:福建農(nóng)林大學(xué)，2008.

［24］ 韓娟娟，程積民，萬惠娥，等.子午嶺遼東櫟群落碳儲量研究［J］.西北林學(xué)院學(xué)報，2010，25(5):18 ～23.

［25］ IPCC.IPCC meeting on current scientific understanding of the p rocesses affecting terrestrial carbon stocks and human influences upon them//Schimel D，Manning M，IPCC Expert Meeting Report.IPCC Secretariat，Geneva，Switzerland，2003:34.

［26］ Kurts Pregitzer，Eugenies Uuskirchen.Carbon cycling and storage in world forests:biome patterns related to forest age［J］.Global Change Biology，2004，(10):2052 ～2077.

［27］ 肖春波，王海，范凱峰，等.崇明島不同年齡水杉人工林生態(tài)系統(tǒng)碳儲量的特點及估測［J］.上海交通大學(xué)學(xué)報(農(nóng)業(yè)科學(xué)版)，2010，28(1):30 ～34.

［28］ 吳鵬飛，朱波，劉世榮，等.不同林齡榿－柏混交林生態(tài)系統(tǒng)的碳儲量及其分配［J］.應(yīng)用生態(tài)學(xué)報，2008，19(7):1419～1424.

［29］ 鮑春生，白艷，青梅，等.興安落葉松天然林生物生產(chǎn)力及碳儲量研究［J］.內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報，2010，31(2):77 ～82.

［30］ Brown S，Lugo A E.The storage and production of organic matter in tropical forests and their role in the global carbon cycle［J］.Biotropica，1982，14:161 ～179.

［31］ 方精云.中國森林生產(chǎn)力及其對全球氣候變化的響應(yīng)［J］.植物生態(tài)學(xué)報，2000，24(5):513 ～517.

［32］ 馬欽彥，陳遐林，王娟，等.華北主要森林類型建群種的含碳率分析［J］.北京林業(yè)大學(xué)學(xué)報，2002，24(5～6):96～100.

［33］ 趙士洞，汪業(yè)勖，于振良，等.中國森林生態(tài)系統(tǒng)碳循環(huán)研究［J］.中國生態(tài)學(xué)會通訊，2000(特刊):50 ～52.

［34］ 劉應(yīng)芳，黃從德，陳其兵.蜀南竹海風(fēng)景區(qū)毛竹林生態(tài)系統(tǒng)碳儲量及其空間分配特征［J］.四川農(nóng)業(yè)大學(xué)學(xué)報，28(2):136～140.

［35］ 方運霆，莫江明，黃忠良，等.鼎湖山馬尾松、荷木混交林生態(tài)系統(tǒng)碳素積累和分配特征［J］.熱帶亞熱帶植物學(xué)報，2003，11(1):47～52.

［36］ 馬煒，孫玉軍，郭孝玉，等.不同林齡長白落葉松人工林碳儲量［J］.生態(tài)學(xué)報，2010，30(17):4659 ～4667.

［37］ 周國模，吳家森，姜培坤.不同管理模式對毛竹林碳貯量的影響［J］.北京林業(yè)大學(xué)學(xué)報，2006，28(6):51 ～54.

［38］ 張駿.中國中亞熱帶東部森林生態(tài)系統(tǒng)生產(chǎn)力和碳儲量研究［D］.浙江:浙江大學(xué)，2008.

［39］ 李正才，徐德應(yīng)，傅懋毅，等.北亞熱帶土地利用變化對土壤有機碳垂直分布特征及儲量的影響［J］.林業(yè)科學(xué)研究，2007，20(6):744 ～749.