四川森林土壤有機(jī)碳儲(chǔ)量估算方法選擇

張炎周，梁玉喜，代立東，王 瀚，易 輝

(1.四川省林業(yè)調(diào)查規(guī)劃院，四川 成都 610081;2.四川省森林資源和荒漠化監(jiān)測(cè)中心，四川 成都 610081;3.四川農(nóng)業(yè)大學(xué)林學(xué)院，四川 雅安 625014)

1 問(wèn)題的提出

1.1 土壤有機(jī)碳及其分布的基本性質(zhì)

陸地生態(tài)系統(tǒng)碳循環(huán)作為重要組成部分，在全球碳收支平衡中占有主導(dǎo)地位［1］，直接影響全球溫度的變化［2］。陸地生態(tài)系統(tǒng)碳平衡及碳貯存與分布的研究日益成為全球變化研究的焦點(diǎn)。準(zhǔn)確地評(píng)估陸地生態(tài)系統(tǒng)碳循環(huán)不僅是準(zhǔn)確地估算未來(lái)大氣CO2濃度、預(yù)測(cè)氣候變化及其對(duì)陸地生態(tài)系統(tǒng)影響的關(guān)鍵，也是履行《聯(lián)合國(guó)氣候變化框架公約》與《京都議定書(shū)》等國(guó)際公約及制定應(yīng)對(duì)策略的關(guān)鍵點(diǎn)。因此，研究土壤有機(jī)碳含量可以反射影響因子的變化，折射特定時(shí)空的自然現(xiàn)象和人類(lèi)活動(dòng)，對(duì)全球生態(tài)環(huán)境研究具有指導(dǎo)意義［3］。

土壤有機(jī)碳(Soil Orgnic Carbon，SOC)儲(chǔ)量是進(jìn)入土壤的生物殘?bào)w等有機(jī)物質(zhì)的輸入與以土壤微生物分解作用為主的有機(jī)物質(zhì)的損失之間的平衡［4］，現(xiàn)有土壤有機(jī)碳的含量是土壤有機(jī)碳分解速率、作物殘余物數(shù)量與組成、植物根系及其他返還至土壤中有機(jī)物的函數(shù)，土壤有機(jī)碳的庫(kù)存量與進(jìn)入土壤的植物凋落物和地上生物量呈線(xiàn)性正相關(guān)關(guān)系［5～7］。土壤有機(jī)碳主要分布于土壤上層1m深度以?xún)?nèi)，全球土壤上層1m內(nèi)的有機(jī)碳含量為1 220 Gt，相當(dāng)于總現(xiàn)存生物量(自然植被和作物)的1.5倍［8］。然而，全球土壤有機(jī)碳儲(chǔ)量在過(guò)去100 a一直呈下降趨勢(shì)［9］。

1.2 森林土壤有機(jī)碳在全球氣候變化中的重要意義

全球土壤有機(jī)碳量(1 500 Pg～2 000 Pg，1 Pg=1015g)約為陸地生物碳量(620 Pg)的2.4倍，為大氣圈CO2含量的兩倍［10］;土壤碳庫(kù)0.1%的變化將導(dǎo)致大氣圈CO2的濃度發(fā)生百萬(wàn)分之一的變化，全球土壤有機(jī)碳10%的變化，其數(shù)量相當(dāng)于人類(lèi)活動(dòng)30 a排放的CO2量［11］。因此土壤有機(jī)碳的動(dòng)態(tài)及其控制過(guò)程的研究不僅是土地資源可持續(xù)利用的重要基礎(chǔ)，而且可以對(duì)土壤碳循環(huán)與全球氣候變化的相互作用研究具有重要意義［12］。

森林在全球碳平衡中起著重要的作用，全球陸地生態(tài)系統(tǒng)中大約75%的碳量?jī)?chǔ)存在土壤中，而森林土壤又保存了地下土壤碳的40%左右［13］，森林土壤碳占全球土壤碳的73%［14］。森林土壤碳含量大約是森林生物量的2倍～3倍。因此對(duì)土壤碳庫(kù)的研究與植被相比是同樣重要的。但目前土壤碳循環(huán)的研究是陸地碳循壞研究中最不充分的部分。這就說(shuō)明對(duì)森林土壤碳排放進(jìn)行深入研究的必要性，同時(shí)也說(shuō)明對(duì)土壤碳庫(kù)做進(jìn)一步研究的重要性［15］。森林土壤碳主要以有機(jī)碳形式存在。由于土壤有機(jī)碳儲(chǔ)量的巨大庫(kù)容，其較小幅度的變化就可能影響到碳向大氣的排放，以溫室效應(yīng)影響全球氣候變化。正確估算區(qū)域土壤有機(jī)碳含量及其變化是全球氣候變化研究以及全球環(huán)境變化研究的重要內(nèi)容。了解土壤碳循環(huán)是研究全國(guó)陸地生態(tài)系統(tǒng)碳循環(huán)的基礎(chǔ)，精確計(jì)算土壤有機(jī)碳儲(chǔ)量的方法是進(jìn)行土壤碳儲(chǔ)量計(jì)算的基礎(chǔ)，確定土壤有機(jī)碳的儲(chǔ)量、空間分布，對(duì)土壤碳循環(huán)，全球變暖現(xiàn)象的研究具有重要意義。

鑒于森林土壤有機(jī)碳在全球氣候變化研究中的重要性，IPCC也把土壤碳庫(kù)/源作為重要的研究?jī)?nèi)容之一［16］。學(xué)術(shù)界也從土壤碳循環(huán)、碳庫(kù)、碳儲(chǔ)量等不同角度進(jìn)行了土壤碳儲(chǔ)量的研究，并對(duì)各方面的研究進(jìn)展進(jìn)行了總結(jié)和評(píng)述［3，15，17～21］，因此本文對(duì)研究進(jìn)程不再贅述，僅對(duì)研究方法進(jìn)行回顧，并力圖從國(guó)際上公認(rèn)為比較合理的方法中選擇出適合四川森林土壤碳研究實(shí)際的理想方法。

2 土壤有機(jī)碳儲(chǔ)量研究的基本方法

研究方法對(duì)于科學(xué)和準(zhǔn)確地估計(jì)土壤有機(jī)碳儲(chǔ)量，減少不確定性的重要意義不言而喻。在較大區(qū)域內(nèi)，例如全國(guó)或全省范圍內(nèi)，尋求統(tǒng)一而準(zhǔn)確的計(jì)量方法不僅是土壤有機(jī)碳儲(chǔ)量研究的必要手段，也是編制IPCC國(guó)家溫室氣體清單的必要步驟，而且對(duì)于解決全球碳循環(huán)研究中的“失匯(missing sink)”問(wèn)題具有重要意義。世界各國(guó)不同研究者對(duì)全球土壤有機(jī)碳庫(kù)存量的估算所用的方法并無(wú)本質(zhì)的區(qū)別，但由于所用資料來(lái)源與土壤分類(lèi)方式的不同，土壤有機(jī)碳庫(kù)存量的估計(jì)值有較大的差異［22～41］。綜合國(guó)內(nèi)外資料，目前國(guó)內(nèi)外有關(guān)土壤有機(jī)碳儲(chǔ)量研究通常采用生命地帶類(lèi)型法、森林類(lèi)型法、土組法、氣候參數(shù)法、碳擬合法、模型法、相關(guān)關(guān)系估算法、統(tǒng)計(jì)估算法、土壤類(lèi)型法等。各類(lèi)研究方法及其各自特點(diǎn)如下［15］:

2.1 按土壤類(lèi)型的研究方法

土壤類(lèi)型法實(shí)際上是土壤分類(lèi)學(xué)方法，通過(guò)土壤剖面數(shù)據(jù)計(jì)算分類(lèi)單元的土壤有機(jī)碳儲(chǔ)量，根據(jù)各種分類(lèi)層次聚合土壤剖面數(shù)據(jù)，再按照區(qū)域或國(guó)家尺度土壤圖上的面積得到土壤有機(jī)碳蓄積總量。Batjes［41］將世界土壤圖劃分為0.5經(jīng)度×0.5緯度的基本面積單元，每個(gè)單元需要土種分布、土壤深度、土壤容重、有機(jī)碳及礫石含量等數(shù)據(jù)，用來(lái)計(jì)算面積單元的平均碳密度。設(shè)j代表地球表面面積網(wǎng)格單元，i代表土層單元，則各個(gè)面積單元j中的平均有機(jī)碳密度Tjd為:

其中ρi為第i層土壤容重，Pi為第i層土壤有機(jī)碳平均儲(chǔ)量，Di為第i層土壤厚度，Si為大于2 mm的平均礫石含量。然后可以推算出全球區(qū)域面積的土壤有機(jī)碳總量:

其中Aj是網(wǎng)格單元j的面積，Tjd是j單元平均有機(jī)碳密度，n為世界土壤圖面積網(wǎng)格單元總數(shù)(259200個(gè))。

Betjes方法需要具備較完整的全球各類(lèi)土壤理化性質(zhì)數(shù)據(jù)，若這項(xiàng)條件能滿(mǎn)足(實(shí)際上難以做到)則統(tǒng)計(jì)結(jié)果相對(duì)較為準(zhǔn)確可靠。

2.2 生命帶研究方法

生命帶法是按生命地帶土壤有機(jī)碳密度與該類(lèi)型分布面積計(jì)算土壤有機(jī)碳蓄積量。Post［42］使用了可反映全球各主要生命帶的2696個(gè)土壤剖面，其中大多數(shù)來(lái)自美國(guó)土壤保持局的數(shù)據(jù)庫(kù)，其余為其本人所發(fā)表的研究結(jié)果。計(jì)算時(shí)對(duì)于沒(méi)有實(shí)測(cè)容重?cái)?shù)據(jù)的土層，其容重根據(jù)土壤有機(jī)碳的密度與深度關(guān)系來(lái)擬合求出:

其中BD為土壤容重，b0、b1、b2為不同植被類(lèi)型下的已知土壤容重和碳密度所確定的常數(shù)，D是從土表到土層中心的深度，Cf為有機(jī)碳質(zhì)量分?jǐn)?shù)。于是單位面積(1 m2)土層的平均碳密度(C)可由下式計(jì)算:

其中δ2mm為直徑大于2 mm的礫石分?jǐn)?shù)，V為土層體積。用碳密度乘以各個(gè)生命帶所對(duì)應(yīng)的土地面積并累加，可得全球土壤有機(jī)碳總儲(chǔ)量(1 m土層深度)。

使用該方法能較為容易地了解不同生命地帶類(lèi)型的土壤有機(jī)碳庫(kù)蓄積總量，而且各類(lèi)型還可以包含多種土壤類(lèi)型，分布范圍更加廣泛，更能反映氣候因素及植被分布對(duì)土壤有機(jī)碳蓄積的影響。但Post方法中全球植被類(lèi)型與面積難以精確統(tǒng)計(jì)，植被與土壤類(lèi)型并不一一對(duì)應(yīng)，加之土地利用方式在人為影響下不斷變化，這樣統(tǒng)計(jì)中不確定因素增多，計(jì)算誤差也會(huì)較大。不過(guò)，在缺乏土壤剖面資料的情況下推算所得結(jié)果仍具有一定意義。

2.3 GIS估算土壤有機(jī)碳儲(chǔ)量

首先用地理信息系統(tǒng)軟件ARC/INFO將一定比例土壤圖數(shù)字化，建立以土屬為單位的空間數(shù)據(jù)庫(kù)，然后計(jì)算各土壤土屬每個(gè)土層的有機(jī)質(zhì)質(zhì)量分?jǐn)?shù):選取該土屬內(nèi)所有土種的典型土壤剖面，按照土壤發(fā)生層分別采集土壤有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)、土層厚度和容重等數(shù)據(jù)計(jì)算出每個(gè)土層的土壤有機(jī)質(zhì)平均質(zhì)量分?jǐn)?shù)和土層平均深度及其平均容重等，并建立土壤有機(jī)質(zhì)的屬性數(shù)據(jù)庫(kù)，利用ARC/INFO的空間分析功能計(jì)算出各類(lèi)土壤的有機(jī)碳儲(chǔ)量。

各類(lèi)土壤的總碳量:

其中i為土壤類(lèi)型，Ci為第i種土壤類(lèi)型的有機(jī)碳儲(chǔ)量(t)，0.58為碳儲(chǔ)量由有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)乘以Remmelen換算系數(shù)，Si為第i種土壤類(lèi)型的面積，Hj為第i種土壤的j層的土屬平均厚度，Qj為第i種土壤j層的土屬平均有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)，Wj為第i種土壤j層的土屬平均容重。

G1S估算方法可以對(duì)土壤圖進(jìn)行較為精確的類(lèi)型劃分，在此基礎(chǔ)上應(yīng)用ARC/INFO的空間分析功能與上面的公式可估算出比一般方法較準(zhǔn)確的土壤有機(jī)碳儲(chǔ)量，并可繪制其空間分布特征圖。

2.4 模型估算土壤有機(jī)碳儲(chǔ)量

國(guó)際上已經(jīng)開(kāi)發(fā)了多種土壤碳循環(huán)的模型［43～47］。模型的類(lèi)型既有相關(guān)關(guān)系模型和機(jī)理過(guò)程模型也有基于實(shí)測(cè)數(shù)據(jù)和遙感數(shù)據(jù)的模型。盡管統(tǒng)計(jì)分析是土壤碳庫(kù)評(píng)價(jià)中最小化空間變異性的可行方法，但模型卻可以將剖面數(shù)據(jù)外推到相似的土壤和生態(tài)區(qū)域，解決尺度轉(zhuǎn)換的問(wèn)題。模型方法最大的限制性因素是缺乏大量相關(guān)和連續(xù)觀(guān)測(cè)的數(shù)據(jù)，使模型的參數(shù)化和初始化更加困難。隨著實(shí)驗(yàn)方法的改善人們可以通過(guò)積累大量土壤碳動(dòng)力學(xué)的信息改善土壤碳模型，以提高管理土壤有機(jī)碳庫(kù)的能力。

2.5 相關(guān)關(guān)系估算法

相關(guān)關(guān)系估算法主要是通過(guò)分析土壤有機(jī)碳蓄積量與采樣點(diǎn)的各種環(huán)境變量、氣候變量和土壤屬性之間的相關(guān)關(guān)系，建立一定的數(shù)學(xué)統(tǒng)計(jì)關(guān)系，從而實(shí)現(xiàn)在有限數(shù)據(jù)基礎(chǔ)上計(jì)算土壤有機(jī)碳蓄積量的目的［4，48～50］。這種方法要求建立的相關(guān)性較高，可以通過(guò)測(cè)采樣點(diǎn)的一些環(huán)境因子來(lái)得到土壤有機(jī)碳蓄積具有方便、省力和簡(jiǎn)單等優(yōu)點(diǎn)。

建立土壤有機(jī)碳含量與降水、溫度、土壤厚度、土壤質(zhì)地、海拔高度和容重之間的相關(guān)關(guān)系是普遍采用的一種方式。然而它們的相關(guān)關(guān)系并非普遍適用，在不同的地方主要控制因素是不同的，各種相關(guān)性表現(xiàn)不一，因此所確定的統(tǒng)計(jì)關(guān)系需要得到檢驗(yàn)和驗(yàn)證，才能在本區(qū)域上應(yīng)用，這是在實(shí)際應(yīng)用中應(yīng)注意的問(wèn)題。

2.6 統(tǒng)計(jì)估算法

用該方法計(jì)算土壤碳庫(kù)的公式如下:

其中，1為第j個(gè)土種的加權(quán)平均有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)(g·kg)，ci為統(tǒng)計(jì)剖面第i層土壤的有機(jī)質(zhì)質(zhì)量分?jǐn)?shù)(g·kg－1)，H為第i層土壤厚度(cm)，c為碳庫(kù)(kg)，ρ為土壤容重，取平均 1.4 t·m－3，2000/3為換算成平方米的系數(shù)，S為第j個(gè)土種的面積(hm2)［51］。

統(tǒng)計(jì)估算法適用于國(guó)家尺度或區(qū)域尺度的土壤碳庫(kù)計(jì)算，是一種較早應(yīng)用并比較成熟的計(jì)算方法。用以上計(jì)算公式可以推出無(wú)機(jī)碳的估算公式，說(shuō)明它的應(yīng)用與適應(yīng)范圍是比較廣的。

2.7 土壤有機(jī)碳儲(chǔ)量估算方法誤差分析方法

土壤有機(jī)質(zhì)是受各種外在和內(nèi)在因素影響而不斷變化的［52］，依靠大量土壤普查獲取的數(shù)據(jù)只是代表當(dāng)時(shí)土壤的狀況，是靜態(tài)數(shù)值，并不能反映土壤的動(dòng)態(tài)變化。土壤有機(jī)碳蓄積量的估算也只能是逐漸逼近真實(shí)值，不可能用一個(gè)值來(lái)代表土壤碳蓄積量。因此分析土壤碳蓄積量的誤差范圍，以及土壤采樣帶來(lái)的隨機(jī)誤差和系統(tǒng)性誤差，是準(zhǔn)確反映土壤有機(jī)碳實(shí)際存儲(chǔ)量的一個(gè)方法。大部分研究者在估算區(qū)域或國(guó)家、全球尺度土壤有機(jī)碳時(shí)分析了不同土壤類(lèi)型或生態(tài)系統(tǒng)類(lèi)型土壤碳庫(kù)的變異系數(shù)［41，53］和誤差范圍［41，54］，認(rèn)為不確定性范圍在 20% ～50%。概括起來(lái)，計(jì)算土壤碳密度誤差項(xiàng)時(shí)有4種統(tǒng)計(jì)方法［41，49，54］:①土壤碳密度 ± 土壤碳密度的標(biāo)準(zhǔn)差(standard error);②土壤碳密度±土壤碳密度平均值的標(biāo)準(zhǔn)差(standard error of means);③土壤碳密度±1/2的土壤碳密度標(biāo)準(zhǔn)差(1/2standard deviation);④首先使用t檢驗(yàn)計(jì)算每個(gè)土壤亞類(lèi)的碳密度精度范圍La(limit of accuracey):

式中，t是顯著性概率0.05水平下的分布值(置信度為95%，可根據(jù)實(shí)際情況設(shè)置不同的置信度);根據(jù)樣本數(shù)查找數(shù)學(xué)手冊(cè)中的t分布表可得t值，Ds是土壤碳密度的標(biāo)準(zhǔn)差，n為自由度(樣本數(shù))。然后乘以亞類(lèi)的面積得到各個(gè)土壤亞類(lèi)碳蓄積量的誤差范圍，加上所有土壤亞類(lèi)的誤差得到國(guó)家尺度土壤有機(jī)碳蓄積量估計(jì)的總誤差范圍。

除了自然原因無(wú)法抗拒和改變之外，土壤分類(lèi)、土壤觀(guān)測(cè)和實(shí)驗(yàn)、數(shù)據(jù)收集、土壤采樣和計(jì)算方法等人為因素也是產(chǎn)生土壤有機(jī)碳蓄積量估算誤差的重要來(lái)源。此估算方法只是基于此而提出的。該估算方法的研究克服了因區(qū)域差異導(dǎo)致精確估算全球土壤碳庫(kù)量的困難，通過(guò)改善區(qū)域估計(jì)而獲得全球土壤有機(jī)碳儲(chǔ)量的估算方法。

2.8 土壤有機(jī)碳儲(chǔ)量估算方法的運(yùn)用實(shí)踐

在眾多估算方法中以Batjes［41］所用的按土壤類(lèi)型的研究方法和Post［42］的按生命帶方法的研究最有代表性。其中土壤類(lèi)型法由于原理簡(jiǎn)單，數(shù)據(jù)較易獲取，是目前國(guó)內(nèi)外土壤碳儲(chǔ)量估算的常用方法［55～57］。土壤類(lèi)型法即土壤分類(lèi)學(xué)方法，通過(guò)兩種方法［58］估算土壤有機(jī)碳儲(chǔ)量:一是利用區(qū)域(甚或世界)土壤圖和土壤分類(lèi)系統(tǒng)來(lái)估算土壤有機(jī)碳儲(chǔ)量;二是利用每個(gè)土壤剖面代表一個(gè)分類(lèi)學(xué)上的單元，土壤分類(lèi)單元的土壤有機(jī)碳含量通過(guò)土壤剖面數(shù)據(jù)計(jì)算，根據(jù)土壤類(lèi)型、土壤亞類(lèi)、土種、土組等分類(lèi)層次聚合土壤剖面數(shù)據(jù)，再按照區(qū)域或國(guó)家尺度土壤圖上的面積得到土壤有機(jī)碳蓄積總量。前者可以形成統(tǒng)一的土壤有機(jī)碳儲(chǔ)量估算體系，便于匯總和對(duì)比;然而忽略了大量區(qū)域或國(guó)家尺度的土壤類(lèi)型細(xì)節(jié)和土壤多樣性，造成土壤類(lèi)型分布圖和土壤屬性信息不同程度上的缺失。后者優(yōu)越性在于提供有關(guān)土壤類(lèi)型的碳儲(chǔ)量及包含土地利用、植被和氣候等等與土壤有機(jī)碳含量變化有關(guān)的因子信息，便于分辨空間格局，了解有機(jī)碳儲(chǔ)量高低原因;局限性在于土壤剖面實(shí)測(cè)數(shù)據(jù)的缺乏和土壤類(lèi)型的空間變異性導(dǎo)致土壤有機(jī)碳儲(chǔ)量估算的不確定性。

在這些方法中，我國(guó)于東升［59］、王義祥［60］用土壤類(lèi)型推算法進(jìn)行了研究。甘海華［61］、邱建軍［62］運(yùn)用模型也作了這方面的研究;童成立［63］等比較了有機(jī)碳計(jì)算機(jī)模擬模型(SCNC)模型和英國(guó)洛桑模型(ROTHC－26.3)，結(jié)果顯示SCNC的接近真值的效果。趙永存［64］等認(rèn)為回歸克里格預(yù)測(cè)土壤有機(jī)碳的空間分布效果最好，能更好地反映碳密度與地形的關(guān)系以及局部變異。方精云、王紹強(qiáng)［65］的研究工作都以同類(lèi)型土壤碳密度的面積加權(quán)平均值作為各類(lèi)型土壤有機(jī)碳密度，再利用土壤類(lèi)型圖統(tǒng)計(jì)出的各類(lèi)型土壤面積來(lái)估算土壤有機(jī)碳總儲(chǔ)量。潘根興［66］先計(jì)算出各土種剖面的土壤有機(jī)碳密度，然后利用各土種的面積統(tǒng)計(jì)資料來(lái)估算土壤有機(jī)碳儲(chǔ)量。李克讓?zhuān)?7］利用0.5°經(jīng)緯網(wǎng)格分辨率的氣候、土壤和植被數(shù)據(jù)驅(qū)動(dòng)的生物地球化學(xué)模型來(lái)進(jìn)行估算。倪?。?7］基于1:400萬(wàn)土壤植被圖以及其它資料，利用BIMEO3模型對(duì)全國(guó)土壤有機(jī)碳儲(chǔ)量進(jìn)行了估算。

2.9 土壤有機(jī)碳儲(chǔ)量估算結(jié)果的不確定性及其產(chǎn)生的主要原因

當(dāng)前土壤有機(jī)碳儲(chǔ)量的研究仍是陸地碳循環(huán)研究中最不充分的部分，土壤有機(jī)碳儲(chǔ)量研究的不確定性極強(qiáng)，對(duì)土壤碳庫(kù)的估計(jì)誤差很大，成為編制IPCC-LULUCF國(guó)家溫室氣體的難點(diǎn)。例如我國(guó)學(xué)者估算的全國(guó)土壤有機(jī)碳儲(chǔ)量結(jié)果從50 Pg到180 Pg，相差3倍～4倍。產(chǎn)生不確定性的主要原因，歸納如下:

(1)土壤本身的復(fù)雜性和不確定性［21］。土壤是一個(gè)不均勻的二維結(jié)構(gòu)體，在空間上呈現(xiàn)復(fù)雜的鑲嵌性，且與氣候以及陸地植被和生物發(fā)生復(fù)雜的相互作用，因此在研究土壤碳循環(huán)，特別是在區(qū)域尺度上的研究仍面臨著大量需要解決的問(wèn)題，甚至在某些測(cè)定上如土壤呼吸以及土壤內(nèi)植物細(xì)根的周轉(zhuǎn)速率等，至今仍沒(méi)有統(tǒng)一和準(zhǔn)確的觀(guān)測(cè)方法。土壤碳氮含量、質(zhì)地、容重等理化性質(zhì)存在很大的空間差異，氣候、母巖、植被和土地利用對(duì)土壤碳庫(kù)容量的綜合影響也很難確定［65］。

(2)土壤分類(lèi)方法和研究方法的差異［15］。土壤分類(lèi)系統(tǒng)的不統(tǒng)一，采樣方法的差異，以及選用不同的土壤碳蓄積量計(jì)算方法和參數(shù)估計(jì)方法使目前的土壤碳蓄積量的估算存在極大的不一致，土壤實(shí)測(cè)數(shù)據(jù)不充分和缺乏連續(xù)、可靠、完整、統(tǒng)一的土壤剖面數(shù)據(jù)也削弱了碳密度量測(cè)的可行性。各地區(qū)土壤厚度和面積統(tǒng)計(jì)資料來(lái)源不同也是蓄積估算不確定性存在的重要原因之一［68］。不同尺度上的影響因子及主要控制因子也存在很大差異，所以由此得到的土壤碳蓄積機(jī)理過(guò)程模擬及其潛在分解、固定和儲(chǔ)存能力分析都會(huì)有所不同。影響土壤碳儲(chǔ)量估算精確度的因素一般有以下幾個(gè)方面:土壤碳密度的空間變異性，容重受土壤含水量、松緊度等因素影響，也具有空間變異性，對(duì)不同土壤或植被類(lèi)型的面積計(jì)算不準(zhǔn)確，植被與土地利用的變化，統(tǒng)計(jì)樣本偏小，缺乏土壤有機(jī)碳的濃度、土壤容重及礫石數(shù)量等資料等等，這都是造成土壤碳儲(chǔ)量估算差異的重要原因。

(3)研究尺度大小與精確性互相矛盾。國(guó)內(nèi)外對(duì)土壤碳庫(kù)的研究多在宏觀(guān)尺度上進(jìn)行，中小尺度的研究較少。多數(shù)研究者或者研究全球和國(guó)家級(jí)陸地生態(tài)系統(tǒng)中土壤碳庫(kù)儲(chǔ)量，或者探討某個(gè)地帶、大區(qū)、大地形類(lèi)型單元、大流域的土壤碳庫(kù)儲(chǔ)量，研究中小尺度地域單元內(nèi)土壤碳庫(kù)的成果尚不多見(jiàn)。研究尺度大，樣本數(shù)量必然不足，必然影響估算結(jié)果的精確性。

(4)統(tǒng)計(jì)方法單一。土壤碳儲(chǔ)量研究一般按植被類(lèi)型、土壤類(lèi)型、生命帶或模型法來(lái)作統(tǒng)計(jì)，不同研究者所用的各種統(tǒng)計(jì)方法并無(wú)本質(zhì)上的差別，通常都是用各種類(lèi)型的平均碳密度乘以相應(yīng)的土地面積并累加獲得土壤碳庫(kù)儲(chǔ)量，實(shí)際上，即便在同一植被類(lèi)型、土壤類(lèi)型或生命帶內(nèi)，土壤碳密度也可能存在較大差異。用平均碳密度值代替實(shí)際值，過(guò)于簡(jiǎn)單。土壤碳作為土壤的有機(jī)組成部分，其含量在空間上是連續(xù)漸變的，理論上可突破類(lèi)型界限，通過(guò)制作土壤碳含量和容重等的空間分布等值線(xiàn)圖，利用地理信息系統(tǒng)強(qiáng)大的分析功能獲得土壤碳庫(kù)儲(chǔ)量。

3 四川森林土壤有機(jī)碳儲(chǔ)量計(jì)量研究的基礎(chǔ)條件和進(jìn)展

四川地處亞熱帶濕潤(rùn)季風(fēng)氣候區(qū)，地域廣闊，地貌結(jié)構(gòu)復(fù)雜，地勢(shì)起伏劇烈，形成復(fù)雜的氣候、植被的區(qū)域性和垂直地帶性特征，在水平地帶性、垂直地帶性分布規(guī)律和錯(cuò)綜復(fù)雜的生物氣候因素、地質(zhì)地貌因素的影響下，形成多種多樣的土壤類(lèi)型。朱鵬飛等根據(jù)基于土壤發(fā)生學(xué)的土壤分類(lèi)學(xué)原理，將四川省的森林土壤分成13個(gè)土綱、34個(gè)土類(lèi)、51個(gè)亞類(lèi)，構(gòu)成四川森林土壤的完整的分類(lèi)系統(tǒng)［69］。為了研究的方便，四川的學(xué)者以土壤生物氣候特征的區(qū)域性差異為指標(biāo)，以山地土壤垂直帶譜及其組合為基礎(chǔ)，在全省共區(qū)分出3個(gè)土壤地區(qū)和土壤地帶，即四川盆地及其周?chē)降貪駶?rùn)亞熱帶森林土壤地區(qū)(黃壤地帶)、川西南山地季節(jié)性干濕交替亞熱帶森林土壤地區(qū)(紅壤地帶)和川西北高山高原半濕潤(rùn)半干旱森林和草甸草原土壤地區(qū)(水平—垂直復(fù)合地帶，山地燥褐土、山地褐土、山地棕壤和高山草甸草原土)。立足自上世紀(jì)50年代以來(lái)的森林土壤調(diào)查與研究，共觀(guān)察記載4 000多個(gè)完整的土壤剖面，完成150多個(gè)土壤剖面理化性質(zhì)分析，輯錄了各區(qū)主要土類(lèi)(山地黃壤、山地紅壤、山地褐土、山地黃棕壤、山地棕壤、山地暗棕壤、山地灰化土、紫色土)的理化性質(zhì)及其植被分布狀況［69～70］。

在森林土壤有機(jī)碳儲(chǔ)量及其空間分布研究領(lǐng)域，以黃從德等的研究為開(kāi)端。黃從德用生命帶類(lèi)型法，根據(jù)四川12類(lèi)森林、約928個(gè)土壤剖面的理化性質(zhì)實(shí)測(cè)數(shù)據(jù)，結(jié)合全省近年來(lái)的森林資源二類(lèi)調(diào)查資料，對(duì)四川森林生態(tài)系統(tǒng)碳儲(chǔ)量進(jìn)行了宏觀(guān)估算和統(tǒng)計(jì)分析，并用GIS制作出森林土壤有機(jī)碳密度的空間分布圖［71～72］。此外，黃從德、楊萬(wàn)勤等對(duì)四川亞高山暗針葉林、中低山人工林土壤有機(jī)碳含量進(jìn)行了研究，所采用的方法依然是生命帶類(lèi)型法或者土壤類(lèi)型法［73～78］。這些研究成果均在不同側(cè)面反映了四川森林土壤有機(jī)碳儲(chǔ)量的某些特征，但是尚未形成系統(tǒng)的、完整的森林土壤有機(jī)碳儲(chǔ)量估算體系。

森林土壤有機(jī)碳儲(chǔ)量研究中存在問(wèn)題及其產(chǎn)生的原因大致有如下方面:

(1)土壤有機(jī)碳儲(chǔ)量研究缺乏區(qū)域范圍內(nèi)完整準(zhǔn)確的土壤剖面實(shí)測(cè)數(shù)據(jù)。因幅員廣闊，地形復(fù)雜多變，地貌類(lèi)型豐富，只是一定區(qū)域內(nèi)的地質(zhì)條件豐富多樣，土質(zhì)較為復(fù)雜，從而進(jìn)行完整的實(shí)測(cè)坡面數(shù)據(jù)較為困難，難以形成一套系統(tǒng)的實(shí)測(cè)數(shù)據(jù)［79］。

(2)土壤碳氮含量、質(zhì)地、容重等土壤性質(zhì)，以及氣候、地形、植被和土地利用的綜合影響存在相當(dāng)大的時(shí)空變異性，導(dǎo)致各個(gè)時(shí)期土壤有機(jī)質(zhì)研究工作者的數(shù)據(jù)資料均缺乏足夠的準(zhǔn)確性，在很大程度上影響了研究成果的科學(xué)性。

(3)土壤采樣方法的設(shè)計(jì)以及土壤碳蓄積量的計(jì)算方法存在較大差異。由于沒(méi)有統(tǒng)一的、系統(tǒng)的研究方法，導(dǎo)致土壤采樣方法的多樣性;同時(shí)，碳儲(chǔ)量的計(jì)算方法各異，計(jì)算結(jié)果不一致。有機(jī)碳儲(chǔ)量的差異性對(duì)研究結(jié)果的準(zhǔn)確性造成影響。

4 四川開(kāi)展土壤有機(jī)碳儲(chǔ)量計(jì)量研究的方法選擇

目前，為履行《聯(lián)合國(guó)氣候變化框架協(xié)議》和《京都議定書(shū)》，我國(guó)即將著手編制2005年國(guó)家溫室氣體清單。根據(jù)《2006年IPCC國(guó)家溫室氣體清單指南》、《國(guó)家溫室氣體清單優(yōu)良作法指南和不確定性管理》(GPG2000)和《土地利用、土地利用變化和林業(yè)優(yōu)良作法指南》(GPG－LULUCF)等文件，與林業(yè)有關(guān)的計(jì)量指標(biāo)主要體現(xiàn)在土地利用、土地利用變化與林業(yè)(land use，land use change and forestry.LULUCF)部分，土地利用變化與林業(yè)清單編制是國(guó)家溫室氣體清單編制的重要方面，編制LULUCF溫室氣體清單的主要任務(wù)，包含:估算2005年森林和其它木質(zhì)生物質(zhì)碳貯量變化;森林轉(zhuǎn)化溫室氣體排放;農(nóng)田和草地土壤碳變化;完成2005年土地利用變化和林業(yè)溫室氣體排放清單報(bào)告。土壤有機(jī)碳是編制LULUCF清單中極為重要的方面。

不僅四川土壤類(lèi)型復(fù)雜多樣，植被類(lèi)型更是極為豐富，同樣的土類(lèi)上往往有多種甚或十余種群系［69～70，80］，因此根據(jù)四川實(shí)際，森林土壤有機(jī)碳儲(chǔ)量計(jì)量的方法，以土壤類(lèi)型法和生命帶法相結(jié)合的辦法為宜，基本步驟如下:

4.1 土壤有機(jī)碳密度估算

土壤有機(jī)碳密度是指單位面積(1 m2)一定厚度的土層中有機(jī)碳儲(chǔ)量。目前的土壤碳密度計(jì)算以100 cm和20 cm的兩個(gè)不同土壤深度為基礎(chǔ)，己成為評(píng)價(jià)和衡量土壤有機(jī)碳儲(chǔ)量的重要指標(biāo)。某些土壤因土體很薄而不能推算至1 m的土層，則按實(shí)際深度計(jì)算。某一土層i的有機(jī)碳密度SOCi(kg·m－2)計(jì)算公式為［26，41，81］:

某一土體的剖面由k層組成，該剖面的有機(jī)碳密度SOCt的計(jì)算公式為:

式中，Ci為土壤有機(jī)碳含量(%)，Di為容重(g·m－3)，Ei為土層厚度(cm)，Gi為大于2 mm 的石礫所占的體積百分比(%);獲得容重和石礫含量數(shù)據(jù)是計(jì)算中的難點(diǎn)。

Ci由下式求得:Ci=SOMi×0.58

式中，SOMi為土壤有機(jī)質(zhì)含量，0.58為Bemmelen換算系數(shù)［82］。

4.2 土壤有機(jī)質(zhì)含量(SOMi)測(cè)定方法

土壤有機(jī)質(zhì)含量的測(cè)定方法主要有儀器法、干燒法、濕氧化法、重鉻酸鉀氧化法和微波密閉消解法等［83～85］。目前比較精確的測(cè)定方法是儀器法，即用High TOC II分析儀測(cè)量土壤中的有機(jī)碳，該方法具有操作簡(jiǎn)單、樣品用量少、分析結(jié)果誤差小等特點(diǎn)，儀器法的測(cè)量精確度高于重鉻酸鉀氧化法和微波消化法［3］。

［1］ Esw Aren H E，Van Den Berg E V，Reich P.Organic carbon in soils of the world.Soil Sci Soc Am J，1993，57:192 ～194.

［2］ 王紹強(qiáng)，朱松麗.中國(guó)土壤有機(jī)碳庫(kù)及空間分布特征分析［J］.地理學(xué)報(bào)，2000，55(5):533 ～544.

［3 岳曼，常慶瑞，王飛，霍艾迪.土壤有機(jī)碳儲(chǔ)量研究進(jìn)展［J］.土壤通報(bào)，2008，39(5):1173 ～1178.

［4］ Post W M，Izaurralde R C，Mann L K，et al.Montoring and verif-ying changes of organic carbon in soil［J］.Climatic Change，2001，51:73～99.

［5］ Larson W E，Pierce F J.Dowdy R H.The threat of soil erosion to long － term crop production［J］.Science，1983，219:458 ～ 465.

［6］ Rasmussen P E，Allmaras R R，Roager Jr N C.Crop residue influences on soil carbon and nitrogen in a wheat-fallow system［J］.Soil Science Society of American Journal，1980，44:596 ～600.

［7］ Follett R P，Gupta S C，Hunt P G.Conservation practices:relation to the mananement of plant nutrients for crop production［C］/Follett R F，Stewart J W B，Cole C V.Soil Fertility and Organic Matter and Critical Components of Production Systems.SSS A Special Publication No.19.American Society of Anronomy.Madison.WI，1987:19～49.

［8］ Sombroek W G，Nachternaele P O，Hebel A.Amounts，dynamics and sequesterinn of carbon in tropical and subtropical soils［J］.AMBIO，1993，22:417 ～425.

［9］ Rodriguez-Murillo J C.Organic carbon content under different types of land use and soil in peninsular Spain［J］.Biol.Fertil.Soils，2001，33:53 ～61.

［10］ Lai R.World soils and the greenhouse effect［J］.Global Change Newsletter，1999，37:4 ～5.

［11］ Eswaren H E，Van Den Berg E V，Reich P.Organic carbon in soils of the world.Soil Sci Soc Am J，1993，57:192 ～194.

［12］ Kirsch Baum Muf.Will changes in soil organic carbon act as a positive or negative feedback on global warming［J］.Biogeochemistry，2000，48:21 ～51.

［13］ Paul K I，Polglase P J，Nyakuengama J G，Khanna P K.Change in soil carbon following afforestation.Forest Ecology and Management，2002，168:241 ～257.

［14］ Sedjo R A.The carbon cycle and global forest ecosystem［J］.Water Air Soil Pollu，1993，70:295 ～307.

［15］ 周?chē)?guó)模，劉恩斌，佘光輝.森林土壤碳庫(kù)研究方法進(jìn)展.浙江林學(xué)院學(xué)報(bào)，2006，23(2):207 ～216.

［16］ IPCC.Climate Chang 2007:The Scientific Basis.Contribution of Working Group I to the Third Assessment Report of IPCC.2007.

［17］ 劉留輝，邢世和，高承芳，崔紀(jì)超，魏多落，胡爭(zhēng)輝，盧超.國(guó)內(nèi)外土壤碳儲(chǔ)量研究進(jìn)展和存在問(wèn)題及展望［J］.土壤通報(bào)，2009，40(3):697 ～701.

［18］ 陳慶強(qiáng)，沈承德，易惟熙，彭少麟，李志安.土壤碳循環(huán)研究進(jìn)展［J］.地球科學(xué)進(jìn)展，1998，13(6):555 ～563.

［19］ 周廣勝，王玉輝.陸地生態(tài)系統(tǒng)類(lèi)型轉(zhuǎn)變與碳循環(huán)［J］.植物生態(tài)學(xué)報(bào)，2002，26(2):250 ～254.

［20］ 朱連奇，朱小立，李秀霞.土壤有機(jī)碳研究進(jìn)展.河南大學(xué)學(xué)報(bào)，2006，36(3):72 ～75.

［21］ 汪業(yè)勖，趙士洞，牛棟.陸地土壤碳循環(huán)的研究動(dòng)態(tài)［J］.生態(tài)學(xué)雜志，1999，18(5):29 ～35.

［22］ Bollv B.E T Degens，P duvIgneasd，et al.Soope，1979，13:1 ～56.

［23］ Schlesinger W H.Carbon storage in the calishe of arid soils:A case study from Arizona［J］.Soil Sci，1982，298:247 ～ 255.

［24］ Bohn，Hinrich.Estimate of organic carbon in world soils［J］.Soil Sci Soc Am J，1976，40:468 ～470.

［25］ Bohn Hl.Estimate of organic carbon in world soils［J］.Soil Sci Soc Am J，1982，46:1118 ～1119.

［26］ Post W M，et al.Soil carbon pools and world life zones［J］.Nature，1982，298:156 ～159.

［27］ Rozhkov V A，et al.Soil carbon estimates and soil carbon map for Russia［J］.Working paper of IIASA，Laxenburg，Austria，1996.

［28］ 汪業(yè)勖，趙士洞.陸地碳循環(huán)研究中的模型方法［J］.應(yīng)用生態(tài)學(xué)雜志，1998，(9)6:658 ～664.

［29］ 楊昕，王明星.陸面碳循環(huán)研究中若干問(wèn)題的評(píng)述［J］.地球科學(xué)進(jìn)展，2001，16(3):427 ～435.

［30］ 金峰，楊浩，蔡祖聰，等.土壤有機(jī)碳密度及儲(chǔ)量的統(tǒng)計(jì)研究［J］.土壤學(xué)報(bào)，2001，38(4):522 ～528.

［31］ 金峰，楊浩，趙其國(guó).土壤有機(jī)碳儲(chǔ)量及影響因素研究進(jìn)展［J］.土壤，2000(1):11 ～17.

［32］ 張東輝，施明恒，金峰，等.土壤有機(jī)碳轉(zhuǎn)化與遷移研究概述［J］.土壤，2000，(6):305 ～309.

［33］ 王淑平，周廣勝，呂育財(cái)，等.中國(guó)東北樣帶(NECT)土壤碳、氮、磷的梯度分布及其與氣候因子的關(guān)系［J］.植物生態(tài)學(xué)報(bào)，2002，26(5):513 ～517.

［34］ 王艷芬，陳佐忠.人類(lèi)活動(dòng)對(duì)錫林郭勒地區(qū)主要草原土壤有機(jī)碳分布的影響［J］.植物生態(tài)學(xué)報(bào)，1998，22(60:545 ～551.

［35］ 潘根興.中國(guó)土壤有機(jī)碳和無(wú)機(jī)碳庫(kù)量研究［J］.科技通報(bào)，1999，15(5):330 ～332.

［36］ 方精云，劉國(guó)華，徐嵩齡.中國(guó)陸地生態(tài)系統(tǒng)的碳循環(huán)及其全球意義［C］//王庚辰，溫玉璞.溫室氣體濃度和排放檢測(cè)及相關(guān)過(guò)程.北京:中國(guó)環(huán)境科學(xué)出版社，1996:129～139.

［37］ 李克讓?zhuān)踅B強(qiáng)，曹明奎.中國(guó)植被和土壤碳儲(chǔ)量［J］.中國(guó)科學(xué)，2003，33(1):72 ～80;.

［38］ Jian N.Carbon storage in terrestrial ecosystems of China［J］.Climatic Change，2001，49:339 ～358.

［39］ Li Z，Zhao Q G.Organic carbon content and distribution in soil under different land uses in tropical and subtropical China［J］.Plant and Soil，2001，231:175 ～185.

［40］ 周玉榮，于振良，趙士洞.我國(guó)主要森林生態(tài)系統(tǒng)碳儲(chǔ)量和碳平衡［J］.植物生態(tài)學(xué)報(bào)，2000，24(5):518 ～522.

［41］ Batjes N H.Total carbon and nitrogen in the soils of the world［J］.Eur Soil Sci，1996，47:l5l～ 163.

［42］ Post W M，Mann L K.Changes in soil organic carbon and nitrogen as a result of cultivation［M］//BOUWMAN A F.Soils and the Greenhouse Effect.England:John Wiley ＆ Sons，1990:410 ～416.

［43］ Jenkinson D S，Adams D E，Wild A.Model estimates of CO2emissions from soil in response to global warming［J］.Nuture，1991，361:304 ～306.

［44］ King A W.Emanuel R W.Wullschleger D S，et al.A search of the missinn carbon sink:a model of terrestrial biospheric response to land use change and atmospheric CO2［J］.Tellus，1995，471:510～519.

［45］ Parton W J，Schimel D S，Cole C V，et al.Analysis of factors controlling soil organic matter levels in the Great Plains grasslands［J］.Soil Sci Soc.Am，J，1987，51:1173 ～1179.

［46］ Parton W J，Rasmussen P E.Longterm effects of crop management in wheat/fallow(II)century model simulations［J］.Soil Sci Soc Am.J，1994，58:530 ～536.

［47］ Vemap M.Vegetation ecosystem modelling and analysis project:comparing biogeography and biogeochemistry models in continental－scale study of terrestrial ecosystem responses to climate change and CO2，doubling［J］.Global Biogeochem cydes，1995，4:407～437.

［48］ Kern J S.Spatial patterns of soil organic carbon in the continuous United States［J］.Soil Sci Soc Arm.J，1994，68:439 ～ 455.

［49］ Franzmeier D P，Lemme G D，Miles R J.Organic carbon in soils of north central United States［J］.Soil Sci Soc Arm.J，1985，49:702～708.

［50］ Burke I C，Yonkr C M，Parton W J，et.al.Texture，climate，and cultivation effects on soil ogranic matter content in U.S.grassland soils［J］.Soil Sci Soc Arm.J，1989，53:800 ～805.

［51］ 楊麗霞.利用森林生態(tài)系統(tǒng)循環(huán)綜合模型模擬土壤有機(jī)碳動(dòng)態(tài)變化［D］.南京:南京農(nóng)業(yè)大學(xué)，2004.

［52］ 方晰，田大倫，項(xiàng)文化，等.杉木人工林土壤有機(jī)碳的垂直分布特征［J］.浙江林學(xué)院學(xué)報(bào)，2004，21(4):418 ～423.

［53］ Batjes N H，Dijkshoorn J A.Carbon and nitrogen stocks in the soils of the Amazon region［J］.Geoderma，1999，89:273 ～286.

［54］ Stone E L，Harris W G，Brown R B，et al.Carbon storage in Florida Snodosols［J］.Soil Sci Soc Am J，1993，67:179 ～182.

［55］ 王紹強(qiáng)，周成虎，劉紀(jì)遠(yuǎn)，等.東北地區(qū)陸地碳循環(huán)平衡模擬分析［J］.地理學(xué)報(bào)，20010，56(4):390 ～400.

［56］ 王紹強(qiáng)，劉紀(jì)遠(yuǎn).土壤碳蓄積變化的影響因素研究現(xiàn)狀［J］.地球科學(xué)進(jìn)展，2002，17(4):528 ～532.

［57］ 王紹強(qiáng)，周成虎，李克讓?zhuān)?中國(guó)土壤有機(jī)碳庫(kù)及空間分布特征分析［J］.地理學(xué)報(bào)，2000，55(5):533 ～544.

［58］ 邵月紅，潘劍君，許信旺，等.淺談土壤有機(jī)碳密度及儲(chǔ)量的估算方法［J］.土壤通報(bào)，2006，37(5):1007 ～1011.

［59］ 于東升，史學(xué)正，孫維俠等.基于1∶100萬(wàn)土壤數(shù)據(jù)庫(kù)的中國(guó)土壤有機(jī)碳密度及儲(chǔ)量研究［J］.應(yīng)用生態(tài)學(xué)報(bào)，2005，16(12):2279～2283.

［60］ 王義祥，翁伯琦.福建省土壤有機(jī)碳密度和儲(chǔ)量的估算［J］.附件農(nóng)業(yè)科學(xué)，2005，20(1):42 ～45.

［61］ 甘海華，吳順輝，范秀月.廣東土壤有機(jī)碳儲(chǔ)量及空間分布特征［J］.應(yīng)用生態(tài)學(xué)報(bào)，2003，14(9):1499 ～1502.

［62］ 邱建軍，王立剛，唐華俊等.東北三省耕地土壤有機(jī)碳儲(chǔ)量變化的模擬研究［J］.中國(guó)農(nóng)業(yè)科學(xué)，2004，37(8):1166 ～1171.

［63］ 童成立，吳金水，向萬(wàn)勝等.長(zhǎng)江中游稻田土壤有機(jī)碳計(jì)算機(jī)模擬［J］.長(zhǎng)江流域資源與環(huán)境，2002，11(3):229 ～233.

［64］ 趙永存，史學(xué)正，于東升等.不同方法預(yù)測(cè)河北省土壤有機(jī)碳密度空間分布特征的研究［J］.土壤學(xué)報(bào)，2005，42(3):379～385.

［65］ 王紹強(qiáng)，劉紀(jì)遠(yuǎn)，于貴瑞.中國(guó)陸地土壤有機(jī)碳蓄積量估算誤差分析［J］.應(yīng)用生態(tài)學(xué)報(bào)，2003，14(6):797 ～802.

［66］ 潘根興.中國(guó)土壤有機(jī)碳和無(wú)機(jī)碳庫(kù)量研究［J］.科技通報(bào)，1999，15(5):330 ～332.

［67］ 倪健.中國(guó)陸地生態(tài)系統(tǒng)碳儲(chǔ)量研究［J］.氣候變化，2001，49:339～358.

［68］ 王紹強(qiáng)，周成虎，羅承文.中國(guó)陸地自然植被碳量空間分布特征探討［J］.地理科學(xué)進(jìn)展，1999，18(3):238 ～245.

［69］ 朱鵬飛，李德融.四川森林土壤［M］.成都:四川科學(xué)技術(shù)出版社，1989.

［70］ 本書(shū)編著委員會(huì).四川森林［M］.北京:中國(guó)林業(yè)出版社，1990.

［71］ 黃從德.四川森林生態(tài)系統(tǒng)碳儲(chǔ)量及其空間分異特征［D］.四川農(nóng)業(yè)大學(xué)，2008.

［72］ 黃從德，張健，楊萬(wàn)勤，張國(guó)慶，王永軍.四川森林土壤有機(jī)碳儲(chǔ)量的空間分布特征［J］.生態(tài)學(xué)報(bào)，2009，29(3):1217～1225.

［73］ 黃從德，張健，楊萬(wàn)勤，張國(guó)慶.川西低山丘陵區(qū)主要人工林土壤有機(jī)碳密度研究［J］.浙江林業(yè)科技，2009，29(3):5 ～8.

［74］ 鄧仁菊，楊萬(wàn)勤，張健，胡建利，馮瑞芳，簡(jiǎn)毅，林靜.川西亞高山森林土壤有機(jī)層碳氮磷儲(chǔ)量特征［J］.應(yīng)用與環(huán)境生物學(xué)報(bào)，2007，13(4):492 ～496.

［75］ 姜發(fā)艷，孫輝，林波，劉慶.川西亞高山云杉人工林恢復(fù)過(guò)程中表層土壤碳動(dòng)態(tài)變化［J］.應(yīng)用生態(tài)學(xué)報(bào)，2009，(20)11:2581～2587.

［76］ 李江，黃從德，張國(guó)慶.川西退耕還林地苦竹林碳密度、碳貯量及其空間分布［J］.浙江林業(yè)科技，2006，26(4):1 ～5.

［77］ 鮮駿仁，張遠(yuǎn)彬，王開(kāi)運(yùn)，胡庭興，楊華.川西亞高山5種森林生態(tài)系統(tǒng)的碳格局［J］.植物生態(tài)學(xué)報(bào)，2009，33(2):283～290.

［78］ 黃從德，張健，楊萬(wàn)勤，張國(guó)慶.四川人工林生態(tài)系統(tǒng)碳儲(chǔ)量特征［J］.應(yīng)用生態(tài)學(xué)報(bào)，2008，19(8):1644 ～1650.

［79］ 朱連奇，朱小立，李秀霞.土壤有機(jī)碳研究進(jìn)展［J］.河南大學(xué)學(xué)報(bào)，2000，36(3):72 ～75.

［80］ 本書(shū)編輯委員會(huì).四川植被［M］.成都:四川人民出版社，1980.

［81］ Schwartz D，Namrim.Mapping the total organic carbon in the soils of the Congo［J］.Global and Planetary Change，2002，33:77 ～93.

［82］ K Zhang，H Dang，S Tan，X Cheng，Q Zhang.Change in soil organic carbon following the‘grain-for-green’programme in china.Land degradation ＆ development，2010，21:13 ～23.

［83］ 劉光崧，蔣能慧，張連第，等.土壤理化分析與剖面描述［M］.北京:中國(guó)標(biāo)準(zhǔn)出版社，1996.

［84］ 高歧，范彩玲，黃曉書(shū).土壤中總糖量的微波測(cè)定法［J］.土壤通報(bào)，1995，26(4):190 ～191.

［85］ 梁重山，黨志，劉從強(qiáng).土壤沉積物樣品中有機(jī)碳含量的快速測(cè)定［J］.土壤學(xué)報(bào)，2002，39(1):135 ～139.