魏文俊 尤文忠 張慧東 魏忠平 趙 剛

(遼寧省林業(yè)科學(xué)研究院，沈陽(yáng)，110032)

土壤有機(jī)碳是陸地生態(tài)系統(tǒng)中最大的碳庫(kù)[1－2]，約為植被碳儲(chǔ)量的 2.5～3 倍[3]，為大氣碳庫(kù)的2～3倍[4]。土壤和凋落物每年以CO2形式進(jìn)入大氣的碳通量是化石燃燒排放 CO2的 10倍[5－6]。森林生態(tài)系統(tǒng)作為陸地生態(tài)系統(tǒng)的主體，其土壤有機(jī)碳主要分布于1 m深度的土層內(nèi)，對(duì)氣候變化和人類活動(dòng)干擾較敏感[7]。氣候變化導(dǎo)致的森林土壤碳庫(kù)的微小變化可能造成全球CO2通量的顯著變化[8－9]，因此，森林土壤碳庫(kù)在全球碳循環(huán)中發(fā)揮著重要的作用[10－11]。

由于受到植被類型、氣候條件、地形變化、土壤屬性、發(fā)育母質(zhì)以及土地利用方式等自然和人為因素的綜合影響，森林土壤有機(jī)碳儲(chǔ)量的分布存在著顯著的區(qū)域不均勻性[12－15]。遼西地區(qū)地處科爾沁沙地東南邊緣，是我國(guó)東北地區(qū)生態(tài)環(huán)境最脆弱地區(qū)[16]，大部分區(qū)域?qū)侔敫珊荡箨懶约撅L(fēng)氣候[17]。位于其邊緣的醫(yī)巫閭山系陰山山脈的余脈，是內(nèi)蒙高原到遼河平原的三大屏障之一[18]，屬于半濕潤(rùn)大陸性季風(fēng)氣候。區(qū)內(nèi)分布的大面積天然油松林是東亞地區(qū)所特有[19]，同時(shí)該區(qū)也屬于我國(guó)油松(Pinus tabulaeformis)天然分布的最北端。醫(yī)巫閭山陽(yáng)坡由于歷史原因，植被被嚴(yán)重破壞，強(qiáng)烈地土壤侵蝕形成了大面積的裸露基巖;陰坡土層深厚植被完好[18]。在生態(tài)系統(tǒng)的物質(zhì)循環(huán)中，土壤碳素和氮素被緊密地聯(lián)系在一起。以醫(yī)巫閭山國(guó)家級(jí)自然保護(hù)區(qū)陰坡的天然油松混交林為研究對(duì)象，分析不同植被組成的天然油松混交林土壤碳氮的分布格局，揭示土壤屬性和植被覆蓋對(duì)土壤碳氮分布規(guī)律的綜合影響，為生態(tài)脆弱邊緣區(qū)對(duì)全球氣候變化的響應(yīng)提供科學(xué)依據(jù)。

1 研究地概況

醫(yī)巫閭山國(guó)家級(jí)自然保護(hù)區(qū)位于遼寧省北寧市和義縣交界處，地理位置 N41°26'～ 41°46'，E121°31'～121°46'。地處陰山山系松嶺山脈，是遼西低山丘陵的一部分，海拔200～800 m，最高峰望海寺峰高867 m。該區(qū)屬于暖溫帶半濕潤(rùn)季風(fēng)氣候，年平均氣溫8.2℃，平均降水量572.5 mm，平均蒸發(fā)量1965 mm，無(wú)霜期160～180 d。母巖多為花崗巖、片麻巖。土壤主要為棕色森林土。醫(yī)巫閭山植被區(qū)劃屬冀遼山地丘陵油松、遼東櫟(Quereus liatungensis)、槲櫟(Q.aliena)林區(qū)，大面積天然油松林為東亞地區(qū)特有[16]。醫(yī)巫閭山天然油松林屬溫性針葉林。林內(nèi)植物種類組成較簡(jiǎn)單，根據(jù)伴生的林下植物種類組成可分為:苔草(Carex nanella)－油松林、胡枝子(Lespedeza hedysaroides)－油松林，迎紅杜鵑(Rhododendron latam)－油松林，荊條(Vitex negundo)－油松林、華北繡線菊(Spiraea fritschiana)－油松林5種類型。

2 研究方法

2.1 樣地布設(shè)與植被調(diào)查

選擇荊條－油松林(YⅠ)、胡枝子－油松林(YⅡ)、華北繡線菊－油松林(YⅢ)3種類型作為研究對(duì)象。每種森林類型隨機(jī)設(shè)置3塊20 m×20 m標(biāo)準(zhǔn)樣地，記錄樣地海拔高度、坡位、坡向、坡度等，調(diào)查喬木的胸徑、樹高、郁閉度等;每個(gè)標(biāo)準(zhǔn)樣地內(nèi)上、中、下坡位置分別設(shè)置3個(gè)2 m×2 m的灌木樣方，調(diào)查灌木的種類、地徑、高度、蓋度等;每個(gè)灌木樣方內(nèi)設(shè)置1 m×1 m的草本樣方，調(diào)查草本的種類、平均高度、蓋度等;標(biāo)準(zhǔn)樣地的基本情況(見表1)。

2.2 土壤調(diào)查與樣品采集

在每個(gè)標(biāo)準(zhǔn)地內(nèi)上、中、下坡分別設(shè)置3個(gè)土壤剖面，按照基本發(fā)生層劃分土壤剖面層次，記錄各個(gè)層次的土壤深度、顏色、結(jié)構(gòu)、質(zhì)地、緊實(shí)度、根量、石礫質(zhì)量分?jǐn)?shù)等;按照[0，10]、(10，20]、(20，40]、(40，60]、(60，100]cm 五個(gè)層次，分別用環(huán)刀、鋁盒采集土樣，用布袋采用四分法采集各層的混合土樣，帶回實(shí)驗(yàn)室，測(cè)定土壤自然含水量、密度、pH值、有機(jī)碳質(zhì)量分?jǐn)?shù)和全氮質(zhì)量分?jǐn)?shù)等。

表1 天然油松林標(biāo)準(zhǔn)樣地概況

利用烘干法測(cè)定土壤自然含水量，先稱取鋁盒土樣鮮重，再將采集的土樣在105℃下恒溫烘干至恒質(zhì)量，計(jì)算出土壤自然含水量;利用環(huán)刀法測(cè)定土壤密度，稱取環(huán)刀土樣鮮質(zhì)量，并根據(jù)烘干法測(cè)定得到的土壤自然含水量，進(jìn)而計(jì)算出土壤密度;用于測(cè)pH值的布袋土樣在室溫下自然風(fēng)干，揀出其中的動(dòng)植物殘?bào)w、根系和石塊等雜物，過2 mm篩，然后利用DENVER酸度計(jì)(UB－7)進(jìn)行測(cè)定;用于測(cè)定土壤有機(jī)碳和全氮的布袋土樣同樣在室溫下自然風(fēng)干，揀出植物殘?bào)w、根系和石塊等雜物，研磨后，過0.149 mm 篩，然后利用KCrO7－濃 H2SO4外加熱法測(cè)定土壤有機(jī)碳質(zhì)量分?jǐn)?shù)，利用FOSS全自動(dòng)凱式定氮儀(Kjeltec 8400)測(cè)定土壤全氮質(zhì)量分?jǐn)?shù)。

2.3 數(shù)據(jù)處理

采用SPSS13.0軟件和Excel軟件進(jìn)行統(tǒng)計(jì)分析。對(duì)不同類型天然油松林和不同土壤剖面深度土壤有機(jī)碳質(zhì)量分?jǐn)?shù)和全氮進(jìn)行顯著性檢驗(yàn)，并用LSD法進(jìn)行多重比較。對(duì)影響表層[0，10]cm土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的土壤屬性(土壤密度、全氮質(zhì)量分?jǐn)?shù)、碳氮比、pH值)和植被覆蓋因子(枯落物有機(jī)碳質(zhì)量分?jǐn)?shù)、全氮質(zhì)量分?jǐn)?shù)、碳氮比)進(jìn)行偏相關(guān)分析。在對(duì)土壤有機(jī)碳質(zhì)量分?jǐn)?shù)進(jìn)行正態(tài)分布檢驗(yàn)的前提下，對(duì)影響表層[0，10]cm土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的因子進(jìn)行逐步回歸分析，進(jìn)而進(jìn)行通徑分析[20]。

3 結(jié)果與分析

3.1 遼西天然油松林土壤有機(jī)碳分布特征

遼西3種類型天然油松林土壤平均有機(jī)碳質(zhì)量分?jǐn)?shù)均隨著剖面深度的增加而降低(表2)，平均有機(jī)碳質(zhì)量分?jǐn)?shù)降低的幅度不盡相同，[0，40]cm土層降低最快，而40 cm以下土層降低趨緩。不同土壤剖面層次內(nèi)3種類型天然油松林土壤有機(jī)碳質(zhì)量分?jǐn)?shù)分布規(guī)律也不同，[0，40]cm土層平均土壤有機(jī)碳質(zhì)量分?jǐn)?shù)從大到小依次為 YⅡ、YⅢ、YⅠ，(40，60]cm 土層從大到小依次為 YⅠ、YⅢ、YⅡ，60 cm 以下土層從大到小依次為 YⅠ、YⅡ、YⅢ。[0，20]cm 土層平均土壤有機(jī)碳質(zhì)量分?jǐn)?shù)YⅡ和YⅢ均顯著大于YⅠ，尤其在[0，10]cm 土層，YⅡ平均土壤含碳量約為YⅠ的2.3倍，YⅢ約為 YⅠ的2倍。

3種類型天然油松林土壤有機(jī)碳變異系數(shù)表現(xiàn)為:[0，20]cm 表層土壤從大到小依次為 YⅠ、YⅢ、YⅡ;20 cm以下土壤從大到小依次為YⅠ、YⅡ、YⅢ(見表2)。在[0，100]cm 土壤剖面內(nèi)，YⅠ有機(jī)碳垂直分布格局的差異性顯著高于YⅡ和YⅢ。

3種類型天然油松林土壤有機(jī)碳方差分析結(jié)果表明。對(duì)不同類型天然油松林而言，YⅠ與YⅡ、YⅢ之間的差異性顯著，而YⅡ與YⅢ之間差異性不顯著;對(duì)不同剖面深度而言，40 cm以下土層剖面深度之間間土壤有機(jī)碳差異不顯著，其它土壤剖面深度之間均呈現(xiàn)顯著的差異(表2)。

表2 遼西3種類型天然油松林土壤有機(jī)碳質(zhì)量分?jǐn)?shù)

3.2 遼西天然油松林土壤全氮分布特征

由表3可知。3種類型天然油松林平均土壤全氮質(zhì)量分?jǐn)?shù)隨著剖面深度的增加而降低，與平均土壤有機(jī)碳質(zhì)量分?jǐn)?shù)變化相似，平均土壤全氮在[0，40]cm土層降低最快，在40 cm以下土層降低趨緩。不同土壤剖面層次內(nèi)，3種類型天然油松林土壤全氮分布規(guī)律也不盡相同，[0，10]cm土層深度，土壤全氮質(zhì)量分?jǐn)?shù)從大到小依次為 YⅡ、YⅢ、YⅠ，(10，60]cm土層深度全氮質(zhì)量分?jǐn)?shù)從大到小依次為YⅢ、YⅡ、YⅠ，60 cm 以下土層全氮質(zhì)量分?jǐn)?shù)從大到小依次為 YⅡ、YⅠ、YⅢ。表層土壤全氮質(zhì)量分?jǐn)?shù)YⅡ、YⅢ均顯著大于 YⅠ，尤其在[0，10]cm 土層內(nèi)，YⅡ約為 YⅠ的1.96倍，YⅢ約為 YⅠ的1.9倍。

在[0，100]cm土壤剖面深度內(nèi)，3種類型天然油松林土壤全氮質(zhì)量分?jǐn)?shù)變異系數(shù)總體上表現(xiàn)為YⅢ＞YⅡ＞YⅠ。與土壤有機(jī)碳的垂直空間異質(zhì)性不盡相同。

土壤全氮質(zhì)量分?jǐn)?shù)方差分析結(jié)果顯示。對(duì)于不同類型天然油松林而言，土壤全氮質(zhì)量分?jǐn)?shù)YⅠ與YⅡ、YⅢ之間的差異性顯著，而YⅡ與YⅢ之間的差異性不顯著，與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)變化規(guī)律一致;對(duì)于不同剖面深度而言，20 cm以下土壤剖面深度之間，土壤全氮差異不顯著，其它剖面深度之間均呈現(xiàn)顯著的差異。

表3 遼西3種類型天然油松林土壤全氮質(zhì)量分?jǐn)?shù)

3.3 遼西3種類型天然油松林表層土壤屬性和植被覆蓋因子與土壤碳質(zhì)量分?jǐn)?shù)的關(guān)聯(lián)度

對(duì)影響遼西地區(qū)3種類型天然油松林表層([0，10]cm)土壤有機(jī)碳質(zhì)量分?jǐn)?shù)(y)的土壤屬性(土壤密度(x1)、土壤全氮質(zhì)量分?jǐn)?shù)(x2)、土壤碳氮比(x3)、土壤pH值(x4))和植被覆蓋因子(枯落物有機(jī)碳質(zhì)量分?jǐn)?shù)(x5)、枯落物全氮質(zhì)量分?jǐn)?shù)(x6)、枯落物碳氮比(x7))進(jìn)行偏相關(guān)分析(見表4)。對(duì)于土壤屬性而言，土壤全氮、碳氮比與有機(jī)碳質(zhì)量分?jǐn)?shù)呈極顯著的正相關(guān)關(guān)系，而土壤密度、pH值呈顯著的負(fù)相關(guān)關(guān)系;對(duì)植被覆蓋因子而言，枯落物全氮質(zhì)量分?jǐn)?shù)、碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)呈顯著的正相關(guān)關(guān)系，枯落物有機(jī)碳質(zhì)量分?jǐn)?shù)則呈顯著的負(fù)相關(guān)系。

表4 土壤有機(jī)碳質(zhì)量分?jǐn)?shù)與其影響因子偏相關(guān)分析

對(duì)表層[0，10]cm土壤有機(jī)碳質(zhì)量分?jǐn)?shù)進(jìn)行正態(tài)性檢驗(yàn)，由于樣本數(shù)n=18屬于小樣本，Shapiro－Wilk檢驗(yàn)方法結(jié)果適用于小樣本檢驗(yàn)，進(jìn)行Shapiro－Wilk 檢驗(yàn)得到統(tǒng)計(jì)量為 0.924，顯著性水平 ɑ=0.424＞0.05，所以表層[0，10]cm 土壤有機(jī)碳質(zhì)量分?jǐn)?shù)服從正態(tài)分布，可以對(duì)其進(jìn)行逐步回歸分析;通過逐步回歸分析方法，從土壤屬性和植被覆蓋因子中逐步選擇和剔除影響因子(見表5)，建立最優(yōu)回歸方程為:y=14.686x2+1.83x3－26.428，R2=0.999，剩余因子 e=0.0316，該值較小，說明土壤全氮質(zhì)量分?jǐn)?shù)和碳氮比兩個(gè)土壤屬性因子對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的影響較大;在此基礎(chǔ)上，進(jìn)一步運(yùn)用通徑分析方法分析土壤全氮質(zhì)量分?jǐn)?shù)和碳氮比對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的直接影響作用和間接影響作用。通徑分析結(jié)果(見表6)。土壤全氮比土壤碳氮比對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的直接影響大。土壤全氮、土壤碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)之間的簡(jiǎn)單相關(guān)系數(shù)分別為0.992、0.793。說明土壤全氮通過土壤碳氮比間接的影響表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)。

表5 土壤有機(jī)碳質(zhì)量分?jǐn)?shù)逐步回歸系數(shù)

表6 土壤有機(jī)碳質(zhì)量分?jǐn)?shù)直接通徑系數(shù)和間接通徑系數(shù)的分解

4 結(jié)論與討論

森林土壤是一個(gè)多相、多界面的復(fù)雜系統(tǒng)。土壤碳氮的變化與植被覆蓋類型、氣候變化、土壤自身屬性、凋落物分解、土壤呼吸以及人類活動(dòng)干擾等眾多相互聯(lián)系和相互影響的生物地球化學(xué)過程密切關(guān)聯(lián)，因此，森林土壤有機(jī)碳分布格局存在著很強(qiáng)的不均勻性。在生態(tài)系統(tǒng)的物質(zhì)循環(huán)中，土壤碳素和氮素被緊密地聯(lián)系在一起。

在[0，100]cm土壤剖面深度內(nèi)，3種類型天然油松林土壤有機(jī)碳質(zhì)量分?jǐn)?shù)和全氮隨著剖面深度的增加而降低。土壤有機(jī)碳、全氮質(zhì)量分?jǐn)?shù)在表層均存在一定厚度富集層，深層質(zhì)量分?jǐn)?shù)有一定下降，質(zhì)量分?jǐn)?shù)低而穩(wěn)定[21－22]，但不隨深度增加而線性降低[22－25]。土壤有機(jī)碳與全氮均存在極顯著線性相關(guān)關(guān)系[26]，其空間分異規(guī)律也表現(xiàn)出一定的相似性[27]，土壤有機(jī)碳空間垂直分異較全氮更大，與白軍紅等[28]研究結(jié)果正好相反。有機(jī)碳和全氮的積累主要存在于地表[29]。遼西3種類型天然油松林土壤有機(jī)碳和全氮主要積累于表層20 cm土層深度內(nèi)，其地表枯落物層厚度可達(dá)6 cm以上，并且油松自身側(cè)根發(fā)達(dá)且分布較淺，根系分布密集層為10～20 cm[30]，可見植被的根系分泌物和凋落物的分解是影響土壤有機(jī)碳循環(huán)的重要驅(qū)動(dòng)力[31]。深層土壤有機(jī)物的分解礦化速率較表層土壤顯著降低[32－33]，其土壤有機(jī)碳氮垂直分布與根系分布及其分泌物密切相關(guān)[34－35]，導(dǎo)致3種類型天然油松林深層土壤有機(jī)碳和全氮差異不顯著，受到母質(zhì)層屬性的影響較大，而其它土層之間的差異性均顯著，有機(jī)碳氮剖面中的分布模式與根系在土層中的垂直分布規(guī)律基本一致[31]。因此，與深層土壤相比，表層土壤對(duì)氣候變化、植被覆蓋狀況和人為干擾等外界因素的影響更敏感，受到影響后，表層土壤的變化更劇烈，而深層土壤則比較穩(wěn)定。

3種類型天然油松林土壤碳氮分布的影響潛力分析。對(duì)土壤屬性而言，土壤全氮、碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)成極顯著的正相關(guān)關(guān)系，而土壤容重、pH值成顯著的負(fù)相關(guān)關(guān)系。對(duì)植被覆蓋而言，地表枯落物全氮、碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)成顯著的正相關(guān)關(guān)系，與枯落物有機(jī)碳質(zhì)量分?jǐn)?shù)成顯著的負(fù)相關(guān)關(guān)系;土壤全氮和碳氮比兩個(gè)土壤屬性對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的影響較大，土壤全氮比土壤碳氮比對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的直接作用更大，土壤全氮通過土壤碳氮比對(duì)土壤有機(jī)碳質(zhì)量分?jǐn)?shù)的間接作用更大。本文研究證明了土壤碳氮空間分異規(guī)律具有相似性。土壤碳氮比是土壤氮素礦化能力的主要標(biāo)志[27]，其比值低有利于微生物在有機(jī)質(zhì)分解過程中的養(yǎng)分釋放和土壤中的有效氮增加，而植物組織中的碳氮比決定了進(jìn)入土壤中的枯落物的分解速率和分解量[36]。土壤pH值通過影響微生物的活動(dòng)而顯著影響著土壤對(duì)碳氮的固定和累積能力，它是影響土壤有機(jī)碳及全氮的空間分布的環(huán)境因子之一，表層土壤pH值與有機(jī)碳質(zhì)量分?jǐn)?shù)呈現(xiàn)顯著負(fù)相關(guān)性[28]。

YⅠ與YⅡ、YⅢ之間各層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)和全氮質(zhì)量分?jǐn)?shù)均存在顯著的差異，而YⅡ與YⅢ之間均沒有差異。[0，20]cm表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)YⅡ、YⅢ明顯大于YⅠ;20 cm土層以下，三者平均土壤有機(jī)碳質(zhì)量分?jǐn)?shù)差異縮小，且逐漸成為YⅠ大于YⅡ、YⅢ。說明植被覆蓋對(duì)表層土壤碳氮具有不可忽視的影響，也證明了表層土壤有機(jī)碳主要來(lái)源于地表枯落物層分解釋放和植被根系分泌物補(bǔ)給[31]。[0，20]cm表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)YⅡ和YⅢ明顯大于YⅠ，可見胡枝子和華北繡線菊兩種灌木對(duì)提高醫(yī)巫閭山天然油松林土壤有機(jī)碳質(zhì)量分?jǐn)?shù)和全氮的作用優(yōu)于荊條，其中胡枝子屬于固氮植物，與劉潔等[37]得到固氮植物對(duì)提高水蝕風(fēng)蝕交錯(cuò)區(qū)土壤全氮質(zhì)量分?jǐn)?shù)的研究結(jié)果一致。

3種天然油松林土壤碳氮分布均隨著剖面深度的增加而降低。而同一土壤剖面層次內(nèi)，三者的平均土壤碳氮質(zhì)量分?jǐn)?shù)變化規(guī)律不盡相同。[0，20]cm土層，平均土壤碳氮質(zhì)量分?jǐn)?shù)YⅡ、YⅢ均顯著大于YⅠ;(20，100]cm土層，三者平均土壤碳氮質(zhì)量分?jǐn)?shù)差異縮小。土壤有機(jī)碳空間異質(zhì)性YⅠ顯著高于YⅡ、YⅢ，而土壤全氮空間分布格局差異性YⅢ高于YⅠ、YⅡ。

不同植被類型，土壤碳、氮質(zhì)量分?jǐn)?shù)均表現(xiàn)為YⅠ與YⅡ、YⅢ之間的差異性顯著，而YⅡ與YⅢ之間差異性不顯著。不同剖面深度，土壤有機(jī)碳質(zhì)量分?jǐn)?shù)，(40，100]cm土壤深度間差異不顯著外，其它剖面深度之間均呈現(xiàn)顯著的差異;土壤全氮質(zhì)量分?jǐn)?shù)，[0，20]cm土壤深度間呈現(xiàn)顯著的差異外，其它土壤剖面深度之間差異均不顯著。

表層土壤全氮質(zhì)量分?jǐn)?shù)、碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)呈極顯著的正相關(guān)關(guān)系，而土壤密度、pH值呈顯著的負(fù)相關(guān)關(guān)系;地表枯落物全氮質(zhì)量分?jǐn)?shù)、碳氮比與土壤有機(jī)碳質(zhì)量分?jǐn)?shù)呈顯著的正相關(guān)關(guān)系，枯落物有機(jī)碳質(zhì)量分?jǐn)?shù)則呈顯著的負(fù)相關(guān)關(guān)系。

土壤全氮質(zhì)量分?jǐn)?shù)和碳氮比對(duì)表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)影響較大，表層土壤有機(jī)碳質(zhì)量分?jǐn)?shù)主要受土壤全氮的影響。

[1]Post W M，Peng T H，Emanuel W R，et al.The global carbon cycle[J].American Scientist，1990，78:310－326.

[2]Batjes N H.Total Carbon and Nitrogen in the soils of the world[J].European Journal of Soil Science，2014，65(1，SI):10－21.

[3]Post W M，Emanuel W R，Zinke P，et al.Soil carbon pools and world life zones[J].Nature，1982，298:156－159.

[4]Davidson E A，Trumbore S E，Amundson R.Soil warming and organic Carbon content[J].Nature，2000，408:789－790.

[5]IPCC，(intergovernmental panel on climate change).scientific assessment of climate change.World Meteorological Organization/U-nited Nations environmental program[M].Cambridge:cambridge university press，1990.

[6]James W，Raich Christopher S，Potter.Global patterns of carbon dioxide emissions from soils[J].Global Biogeochemical Cycles，1995，9(1):23－36.

[7]Fang C，Smith P，Moncrieff J B，et al.Similar response of labile and resistant soil organic matter pools to changes in temperature[J].Nature，2005，433:57－59.

[8]Rustad L E，Huntington T G，Boone R D.Controls on soil respiration:Implications for climate change[J].Biogeochemistry，2000，48(1):1－6.

[9]Janzen H H.Carbon cycling in earth systems-a soil science perspective[J].Agriculture Ecosystems ＆ Environment，2004，104(3):399－417.

[10]Zhou C，Zhou Q，Wang S.Estimating and analyzing the spatial distribution of soil organic Carbon in China[J].Ambio，2003，32(1):6－12.

[11]Wang S H，Tian H Q，Liu J Y，et al.Pattern and change of soil organic carbon storage in China:1960s－ 1980s[J].Tellus B，2003，55(2):416－427.

[12]鄭敬剛，吳國(guó)璽，何明珠，等.阿拉善荒漠區(qū)植物多樣性與土壤理化性質(zhì)相關(guān)性研究[J].干旱區(qū)資源與環(huán)境，2009，23(3):151－155.

[13]尉海東，馬祥慶，劉愛琴，等.森林生態(tài)系統(tǒng)碳循環(huán)研究進(jìn)展[J].中國(guó)生態(tài)農(nóng)業(yè)學(xué)報(bào)，2007，15(2):188－192.

[14]李春莉，趙萌莉，韓國(guó)棟，等.不同放牧壓力下荒漠草原土壤有機(jī)碳特征及其與植被之間關(guān)系的研究[J].干旱區(qū)資源與環(huán)境，2008，22(5):134－138.

[15]鄭淑華，趙萌莉，韓國(guó)棟，等.不同放牧壓力下典型草原土壤物理性質(zhì)與植被關(guān)系的研究[J].干旱區(qū)資源與環(huán)境，2005，19(z1):199－203.

[16]張寶東，孫占祥，白偉，等.基于灰色關(guān)聯(lián)度分析的遼西地區(qū)典型種植模式評(píng)價(jià)[J].生態(tài)學(xué)雜志，2010，29(12):2452－2456.

[17]王江山.遼西北地區(qū)氣候與生態(tài)環(huán)境的相互影響[J].安徽農(nóng)業(yè)科學(xué)，2009，37(31):15376－15378.

[18]于晶，曹同，郭水良，等.醫(yī)巫閭山自然保護(hù)區(qū)苔蘚植物區(qū)系成分與地理分布特征研究[J].植物研究，2001，21(1):38－41.

[19]曹偉，劉云龍，朱彩霞，等.遼寧醫(yī)巫閭山種子植物區(qū)系研究[J].植物研究，2004，24(2):152－157.

[20]明道緒.通徑分析[M].雅安:四川農(nóng)業(yè)大學(xué)出版社，1990.

[21]孫志高，劉景雙，于君寶.三江平原小葉章濕地土壤中堿解氮和全氮含量的季節(jié)變化特征[J].干旱區(qū)資源與環(huán)境，2009，23(8):145－149.

[22]黨亞愛.黃土高原南北主要類型土壤有機(jī)碳氮庫(kù)分布特征研究[D].楊凌:西北農(nóng)林科技大學(xué)，2008.

[23]李玉強(qiáng)，趙哈林，移小勇，等.沙漠化過程中科爾沁沙地植物－土壤系統(tǒng)碳氮儲(chǔ)量動(dòng)態(tài)[J].環(huán)境科學(xué)，2006，27(4):635－640.

[24]郝玉光，劉芳，張偉華，等.烏蘭布和沙區(qū)人工綠洲土壤養(yǎng)分動(dòng)態(tài)研究[J].干旱區(qū)資源與環(huán)境，2009，23(7):176－181.

[25]邢雅玲.賽罕烏拉自然保護(hù)區(qū)不同植被類型下土壤碳氮特征分析[D].呼和浩特:內(nèi)蒙古農(nóng)業(yè)大學(xué)，2011.

[26]馬玉紅，郭勝利，楊雨林，等.植被類型對(duì)黃土丘陵區(qū)流域土壤有機(jī)碳氮的影響[J].自然資源學(xué)報(bào)，2007，22(1):97－105.

[27]胡啟武，歐陽(yáng)華，劉賢德.祁連山北坡垂直帶土壤碳氮分布特征[J].山地學(xué)報(bào)，2006，24(6):654－661.

[28]白軍紅，鄧偉，朱顏明，等.霍林河流域濕地土壤碳氮空間分布特征及生態(tài)效應(yīng)[J].應(yīng)用生態(tài)學(xué)報(bào)，2003，14(9):1494－1498.

[29]賈曉紅，李新榮，李元壽.干旱沙區(qū)植被恢復(fù)中土壤碳氮變化規(guī)律[J].植物生態(tài)學(xué)報(bào)，2007，31(1):66－74.

[30]劉春江，郭旭臨，徐振泉.北京西山地區(qū)人工油松栓皮櫟混交林根系研究初報(bào)[J].北京林學(xué)院學(xué)報(bào)，1985(1):77－84.

[31]孫文義，郭勝利.天然次生林與人工林對(duì)黃土丘陵溝壑區(qū)深層土壤有機(jī)碳氮的影響[J].生態(tài)學(xué)報(bào)，2010，30(10):2611－2620.

[32]Franzluebbers A J，Hons F M，Zuberer D A.Seasonal changes in soil microbial biomass and mineralizable C and N in wheat management systems[J].Soil Biology and Biochemistry，1994，26(11):1469－1475.

[33]Van Dam D，Van Breemen N，Veldkamp E.Soil organic Carbon dynamics:variability with depth in forested and deforested soils under pasture in Costa Rica[J].Biogeochemistry，1997，39(3):343－375.

[34]程云環(huán)，韓有志，王慶成，等.落葉松人工林細(xì)根動(dòng)態(tài)與土壤資源有效性關(guān)系研究[J].植物生態(tài)學(xué)報(bào)，2005，29(3):403－410.

[35]武春華，陳云明，王國(guó)梁，等.黃土丘陵區(qū)典型植物群落根系垂直分布與環(huán)境因子的關(guān)系[J].中國(guó)水土保持科學(xué)，2008，6(3):65－70，87.

[36]高志勤，傅懋毅.不同毛竹林土壤碳氮養(yǎng)分的季節(jié)變化特征[J].浙江林學(xué)院學(xué)報(bào)，2006，23(3):248－254.

[37]劉潔，劉小林，張興昌，等.水蝕風(fēng)蝕交錯(cuò)區(qū)不同類型植被對(duì)土壤碳氮影響的研究[J].水土保持研究，2008，15(6):53－56.