賈丙瑞，周廣勝，2，* ，蔣延玲，王 宇，殷曉潔，胡天宇

(1.中國科學(xué)院植物研究所植被與環(huán)境變化國家重點(diǎn)實(shí)驗(yàn)室，北京 100093;2.中國氣象科學(xué)研究院，北京 100081)

北方林又稱泰加林，分布于45°—70°N之間的寒溫帶地區(qū)，橫貫歐亞和北美大陸的北部。北方林向南延伸到我國大興安嶺北部，即我國的寒溫針葉林。作為我國寒溫針葉林的主要植被類型，興安落葉松林面積占全國森林面積的30%，被認(rèn)為在碳匯中占據(jù)重要份額［1］。北方林是地球上第二大陸地生物群區(qū)，來自地面植被觀測(cè)、大氣CO2和O2濃度監(jiān)測(cè)、衛(wèi)星遙感信息應(yīng)用、生態(tài)和大氣模型模擬等方面的研究均表明，北半球中高緯度的森林生態(tài)系統(tǒng)是一個(gè)巨大的碳匯，固定了全球碳循環(huán)中大部分“未知碳匯”［2-3］。

彭少麟等［4］研究表明，土壤碳大量積累并不是來自大的輸入量，而是其它因素(如低溫)對(duì)分解者的限制起主導(dǎo)作用。北方林生態(tài)系統(tǒng)地處寒溫帶氣候區(qū)，因此，北方林生態(tài)系統(tǒng)擁有全球最大的土壤碳庫，約占全球土壤碳庫總量的23.4%［5］。已有研究表明，溫度是寒溫針葉林生態(tài)系統(tǒng)土壤呼吸作用的主要調(diào)控因子［6-11］。王慶豐等［12］對(duì)移栽自不同緯度的興安落葉松研究表明，土壤呼吸作用對(duì)溫度的敏感性隨緯度升高而增加;其他研究也發(fā)現(xiàn)，寒冷氣候區(qū)土壤呼吸作用對(duì)溫度升高更加敏感［13-15］。因此，北方林生態(tài)系統(tǒng)因其巨大的土壤碳庫和對(duì)溫度變化的響應(yīng)極為敏感而備受關(guān)注。我國大興安嶺林區(qū)近50年來氣溫呈顯著增加趨勢(shì)［16-18］。如果全球持續(xù)變暖，寒溫針葉林生態(tài)系統(tǒng)潛在的巨大碳庫將可能成為大氣CO2的重要來源。因此，寒溫針葉林在氣候變化與全球碳收支中具有十分重要的作用［19］。

土壤呼吸作用作為土壤碳主要輸出途徑和重要的大氣CO2源，對(duì)其精確測(cè)定已成為全球變化研究中的關(guān)鍵問題之一。目前，土壤呼吸作用野外原位測(cè)定通常依據(jù)觀測(cè)樣方大小不同隨機(jī)選取3—5個(gè)觀測(cè)樣點(diǎn)［7，12，20-24］，少數(shù)研究設(shè)置更多的觀測(cè)樣點(diǎn)數(shù)［9，25］。已有研究表明，寒溫針葉林生態(tài)系統(tǒng)土壤呼吸作用空間變動(dòng)性很大［8，10，20，26-28］。土壤呼吸作用的空間變化增加了其測(cè)定的不確定性，為準(zhǔn)確估算土壤呼吸作用，必須解決土壤呼吸作用小尺度上的空間變異性［29］。同時(shí)，土壤呼吸作用時(shí)空動(dòng)態(tài)變化特征及其驅(qū)動(dòng)因子研究對(duì)于預(yù)測(cè)生態(tài)系統(tǒng)對(duì)氣候變化的響應(yīng)至關(guān)重要［8，22］。本研究對(duì)大興安嶺原始寒溫針葉林生態(tài)系統(tǒng)土壤呼吸作用的日、季動(dòng)態(tài)特征及其影響因子進(jìn)行了分析，并比較同一區(qū)域不同觀測(cè)樣方個(gè)數(shù)土壤呼吸作用的變異性及其取樣代表性，為深入研究寒溫針葉林生態(tài)系統(tǒng)碳動(dòng)態(tài)特征和準(zhǔn)確估算土壤碳排放量提供科學(xué)依據(jù)與基礎(chǔ)資料。

1 材料和方法

1.1 區(qū)域概況

本研究在中國科學(xué)院植物研究所中國北方林森林生態(tài)系統(tǒng)定位研究站(51°47'N，123°01'E，海拔773 m)進(jìn)行。該站位于我國最大的寒溫針葉林保護(hù)區(qū)——黑龍江呼中國家級(jí)自然保護(hù)區(qū)內(nèi)，屬寒溫帶大陸性季風(fēng)氣候。年均氣溫－4.4℃，在－2.7— －5.3℃之間波動(dòng)。降水量為481.6 mm，主要集中在6—8月，無霜期80—100 d。棕色針葉林土是該地區(qū)最具代表性的土壤類型，地貌類型為大興安嶺北部石質(zhì)中低山山地，坡度平緩，一般坡度在15°以下。地帶性植被類型為寒溫針葉林，以興安落葉松(Larix gmelinii)為單優(yōu)勢(shì)樹種，并混有少量白樺(Betula platyphylla)，平均樹高15 m，郁閉度0.52，林下灌木主要有杜香(Ledum palustre)、越桔(Vaccinium vitis-idaea)等。

1.2 研究方法

采用便攜式光合作用測(cè)量系統(tǒng)Li-6400(LI-COR Inc.，USA)配合土壤氣室(Li-6400-09)測(cè)定土壤呼吸作用，土壤呼吸作用測(cè)定前一天將10 cm直徑PVC環(huán)插入土壤表面下2 cm，保持不動(dòng)，并剪掉植被地上部分，共設(shè)置9個(gè)重復(fù)，每個(gè)重復(fù)連續(xù)測(cè)定3次，取其平均值。2006年和2007年6月至10月，每月進(jìn)行1次土壤呼吸作用測(cè)定，其中2006年7月5—6日和8月10—11日、2007年7月22—23日進(jìn)行了土壤呼吸作用全天日動(dòng)態(tài)測(cè)定，其他為白天(6:00—18:00)日動(dòng)態(tài)測(cè)定，每隔2 h測(cè)定1次。為了解土壤呼吸作用的空間異質(zhì)性及其取樣代表性，選擇平坦的樣地布設(shè)40個(gè)觀測(cè)樣點(diǎn)，分為5列8行，列和行間距均為2 m，2009年7月4日分別于3:00、7:00和11:00測(cè)定土壤呼吸作用。測(cè)定完成后收集40個(gè)觀測(cè)點(diǎn)的地表凋落物，并挖取每個(gè)樣點(diǎn)0—30 cm土壤，用0.5 mm孔徑篩沖洗分別揀出土層中的粗根(＞2 mm)和細(xì)根(≤2 mm)，帶回實(shí)驗(yàn)室烘干(75℃)稱重。氣象觀測(cè)塔距觀測(cè)樣地約100 m，每半小時(shí)自動(dòng)記錄1次，包括1.5 m氣溫(HMP45C，Vaisala，F(xiàn)inland);表層、5 cm 和 10 cm 土壤溫度(107 L，Campbell Scientific Inc.，USA)等。

1.3 統(tǒng)計(jì)分析

運(yùn)用SPSS統(tǒng)計(jì)軟件雙變量相關(guān)分析中的Pearson相關(guān)分析，進(jìn)行土壤呼吸作用與不同層次溫度在日、季尺度的相關(guān)性比較。為求取土壤呼吸作用、根系和凋落物40個(gè)觀測(cè)樣點(diǎn)不同取樣個(gè)數(shù)的變異性，分別以1個(gè)樣點(diǎn)觀測(cè)值以及相鄰2個(gè)、3個(gè)，…，10個(gè)樣點(diǎn)平均值計(jì)算變異系數(shù)，用于描述樣地空間異質(zhì)性造成的取樣變異性，變異系數(shù)=標(biāo)準(zhǔn)差/平均值×100。

2 結(jié)果和分析

2.1 氣溫與地溫的日、季動(dòng)態(tài)變化

寒溫針葉林生態(tài)系統(tǒng)6月至10月1.5 m氣溫日動(dòng)態(tài)呈單峰型曲線(圖1)，最高值出現(xiàn)在午后14:00左右，最低值通常出現(xiàn)在 4:30 左右;6—10 月月均溫分別為 14.47、15.81、16.15、5.29、0.19 ℃，6 月、7 月和 8 月氣溫逐漸升高，9月和10月氣溫則大幅降低且夜間氣溫降到0℃以下。圖1給出了寒溫針葉林生態(tài)系統(tǒng)6—10月5 cm土壤溫度日動(dòng)態(tài)變化，其中6—9月呈單峰型曲線，10月變幅很小，基本在0℃上下波動(dòng)。6—9月5 cm土壤溫度最高值出現(xiàn)在16:00左右，滯后于氣溫2 h;最低值出現(xiàn)在6:00左右，滯后于氣溫約1.5h。6—8月1.5 m月平均氣溫比5 cm土壤溫度高約5℃，9月和10月兩者月均溫則接近。

圖1 6—10月1.5 m氣溫和5 cm土溫月平均日動(dòng)態(tài)Fig.1 Mean diurnal dynamics of air temperature at 1.5 m height and soil temperature at 5 cm depth from June to October.Dots are monthly means with standard error bars

2.2 土壤呼吸作用的日、季動(dòng)態(tài)及其影響因子

圖2給出了2006年7月5—6日和8月10—11日、2007年7月22—23日土壤呼吸作用日動(dòng)態(tài)變化，一天內(nèi)土壤呼吸作用均呈單峰型曲線變化，高峰值出現(xiàn)在16:00左右，最低值出現(xiàn)在6:00左右，日變動(dòng)幅度分別為 2.88、5.31 μmol CO2m－2s－1和 4.42 μmol CO2m－2s－1，與 5 cm 土壤溫度日動(dòng)態(tài)變化趨勢(shì)相近，滯后于氣溫日動(dòng)態(tài)變化(圖1)。10:00和22:00的土壤呼吸作用最接近一天的平均值。因此，可以選擇10:00左右測(cè)定的土壤呼吸作用平均值作為日平均值。姚玉剛等［30］對(duì)熱帶雨林研究也表明，9:00—11:00土壤呼吸作用接近日平均值。

圖2 2006年和2007年土壤呼吸作用日動(dòng)態(tài)變化Fig.2 Diurnal dynamics of soil respiration rate on Jul.5—6，2006，Aug.10—11，2006 and Jul.22—23，2007

圖3給出了2006年和2007年生長季節(jié)6—10月土壤呼吸作用動(dòng)態(tài)變化，6月份隨著天氣轉(zhuǎn)暖和植被生長土壤呼吸作用逐漸增大，8月份氣溫最高時(shí)土壤呼吸作用也達(dá)到最大值，2006年和2007年分別為8.19μmol CO2m－2s－1和6.89 μmol CO2m－2s－1，隨后 9 月份氣溫迅速下降土壤呼吸作用也降低，到10月份氣溫波動(dòng)于0℃左右，寒溫針葉林生態(tài)系統(tǒng)只維持微弱的土壤呼吸作用，約為 0.61 μmol CO2m－2s－1。2006 年和 2007 年 6—10月土壤呼吸作用季節(jié)動(dòng)態(tài)變化幅度分別為7.38 μmol CO2m－2s－1和 6.48 μmol CO2m－2s－1。興安落葉松人工林土壤呼吸作用與本研究結(jié)果相近，6—10月分別為 4.30、4.88、5.88、3.47、2.76 μmol CO2m－2s－1，最大值出現(xiàn)在8月份［31］。小興安嶺興安落葉松沼澤林［24］和寒溫帶落葉林［21］也表現(xiàn)出相似的季節(jié)變化特征。

圖3 2006年和2007年土壤呼吸作用季節(jié)動(dòng)態(tài)變化Fig.3 Seasonal dynamics of soil respiration rate in 2006 and 2007

土壤呼吸作用平均值與氣溫和3個(gè)層次土壤溫度相關(guān)性分析表明(表1)，在日尺度上，土壤呼吸作用與3個(gè)層次土壤溫度的相關(guān)性均達(dá)到顯著和極顯著水平，Pearson相關(guān)系數(shù)(R)在0.737以上;土壤呼吸作用與氣溫的相關(guān)系數(shù)為0.502—0.868，土壤呼吸作用與土壤溫度的相關(guān)性好于氣溫。在季節(jié)尺度上，2006年土壤呼吸作用與土壤溫度和氣溫的相關(guān)性均達(dá)到極顯著水平，其中與3個(gè)層次土壤溫度的相關(guān)系數(shù)在0.846以上，而與氣溫的相關(guān)系數(shù)為0.687;2007年土壤呼吸作用與3個(gè)層次土壤溫度相關(guān)性均達(dá)到極顯著水平且相關(guān)系數(shù)在0.615以上，與氣溫的相關(guān)性達(dá)到顯著水平，相關(guān)系數(shù)為0.513。因此，不論在日尺度上還是季節(jié)尺度上，寒溫針葉林生態(tài)系統(tǒng)土壤呼吸作用主要受土壤溫度調(diào)控，尤其是5 cm土壤溫度。對(duì)興安落葉松人工林［32］和興安落葉松沼澤林［24］研究均表明，土壤溫度比氣溫更能準(zhǔn)確地反映土壤呼吸作用的動(dòng)態(tài)變化。

表1 土壤呼吸作用與不同層次溫度的Pearson相關(guān)系數(shù)(R)Table 1 Pearson correlation coefficients(R)of soil respiration rate with temperature at different layers

2.3 土壤呼吸作用的空間分布特征

基于3:00、7:00和11:00以2 m間距(5列8行)測(cè)定的40個(gè)觀測(cè)樣點(diǎn)土壤呼吸作用，分別以1個(gè)樣點(diǎn)觀測(cè)值以及相鄰2個(gè)、3個(gè)、…、10個(gè)樣點(diǎn)平均值計(jì)算的土壤呼吸作用變異系數(shù)如圖4。1個(gè)觀測(cè)樣點(diǎn)變異系數(shù)即40個(gè)樣方在2 m間距(5列8行)范圍內(nèi)的空間變異性，3:00、7:00和11:00分別為35.5%、27.6%和23.0%，說明該研究地點(diǎn)以10 cm直徑圓形樣方測(cè)定土壤呼吸作用時(shí)，樣方間存在較大的空間變異性。3:00、7:00 和11:00 變動(dòng)幅度分別為5.97、4.96 μmol CO2m－2s－1和5.84 μmol CO2m－2s－1。隨著觀測(cè)樣點(diǎn)個(gè)數(shù)增加，土壤呼吸作用變異系數(shù)呈逐漸遞減趨勢(shì)，3:00由35.5%減少到19.8%，7:00由27.6%減少到11.1%，11:00由23.0%減少到9.1%。3個(gè)觀測(cè)樣點(diǎn)時(shí)土壤呼吸作用變異系數(shù)迅速降低，3:00、7:00和11:00變異系數(shù)分別為27.0%、16.8%和 13.0%;6 個(gè)觀測(cè)樣點(diǎn)時(shí)，3:00、7:00 和 11:00 變異系數(shù)分別為 22.8%、11.1%和 11.8%，基本接近7—10個(gè)觀測(cè)樣點(diǎn)時(shí)的最小變異系數(shù)。因此，以10 cm直徑圓形樣方測(cè)定時(shí)，6個(gè)觀測(cè)樣點(diǎn)基本能夠反映寒溫針葉林的土壤呼吸作用變化情況。

對(duì)于40個(gè)樣方的根系和凋落物來說，與土壤呼吸作用表現(xiàn)出相似的空間變異性(圖4)。其中粗根的空間變異性最大，達(dá)到103.8%;總根系和凋落物次之，分別為60.3%和47.1%;細(xì)根空間變異性最小，為32.1%。隨著觀測(cè)樣點(diǎn)個(gè)數(shù)增加，變異系數(shù)均逐漸降低，且降低幅度按粗根、總根、凋落物和細(xì)根順序依次遞減，3個(gè)觀測(cè)樣點(diǎn)時(shí)變異系數(shù)分別為 53.1%、33.4%、34.1%和 18.2%，6個(gè)觀測(cè)樣點(diǎn)時(shí)變異系數(shù)分別為 37.8%、27.3%、25.7%和16.7%，與7—10個(gè)觀測(cè)樣點(diǎn)時(shí)的變異系數(shù)基本持平。

3:00、7:00和11:00土壤呼吸作用與細(xì)根、粗根、總根及凋落物的相關(guān)性分析表明(圖5)，3個(gè)時(shí)段土壤呼吸作用均與細(xì)根的相關(guān)性最好，相關(guān)系數(shù)(R2)分別為0.351、0.396和0.322，說明該區(qū)域土壤呼吸作用的空間變異性主要由于細(xì)根差異引起，粗根和凋落物的影響并不大。

3 討論

本研究土壤呼吸作用最大值出現(xiàn)在夏季生長旺盛期8月份氣溫最高時(shí)，2006年和2007年分別為8.19 μmol CO2m－2s－1和 6.89 μmol CO2m－2s－1。興安落葉松林 4 大類常見微生物(細(xì)菌、放線菌、真菌、自生固氮菌)最大數(shù)量也同樣出現(xiàn)在夏季［33］。10月份氣溫波動(dòng)于0℃左右，寒溫針葉林生態(tài)系統(tǒng)只維持微弱的土壤呼吸作用，約為 0.61 μmol CO2m－2s－1，接近于解凍前興安落葉松林的土壤呼吸作用平均值(0.5 μmol CO2m－2s－1)［12］。分布于其他緯度的北方林生態(tài)系統(tǒng)土壤呼吸作用變動(dòng)范圍為 0.4—6.9 μmol CO2m－2s－1(48°13'N)［10］、0.6—9.3 μmol CO2m－2s－1(53°42'N)［6］、0.1—7 μmol CO2m－2s－1(53°59'N)［9］，本研究結(jié)果基本接近這些范圍。

圖4 土壤呼吸作用變異系數(shù)、根系和凋落物變異系數(shù)與重復(fù)次數(shù)的關(guān)系Fig.4 Relationship of sample number with coefficients of variation for soil respiration rate，root and litter

圖5 土壤呼吸作用與細(xì)根、粗根、總根及凋落物的關(guān)系Fig.5 Relationships of soil respiration rate with fine root，coarse root，total root and litter

本研究區(qū)域40個(gè)樣方土壤呼吸作用的空間變異幅度為4.96—5.97μmol CO2m－2s－1，大于或接近土壤呼吸作用日動(dòng)態(tài)變化幅度(2.88—5.31 μmol CO2m－2s－1)，而小于土壤呼吸作用季節(jié)動(dòng)態(tài)變化幅度(6.48—7.38 μmol CO2m－2s－1)，與Buchmann［20］研究結(jié)果相近。以10 cm直徑圓形樣方測(cè)定土壤呼吸作用時(shí)，40個(gè)樣點(diǎn)在3:00、7:00和11:00 3個(gè)時(shí)間段的變異系數(shù)分別為35.5%、27.6%和23.0%，說明寒溫針葉林生境下用10 cm直徑圓形樣方取樣時(shí)各樣方間的空間變異性較大，通過增加觀測(cè)樣點(diǎn)個(gè)數(shù)可使土壤呼吸作用變異系數(shù)迅速降低，6個(gè)觀測(cè)樣點(diǎn)基本能夠反映寒溫針葉林的土壤呼吸作用變化情況。Fang等［34］佛羅里達(dá)濕地松人工林研究表明，在25×25 m樣地內(nèi)12個(gè)樣方(直徑7.1cm)的土壤呼吸作用變異系數(shù)高達(dá)55%。Buchmann［20］對(duì)挪威云杉林同樣研究表明，夏季生長旺盛期土壤呼吸作用最大時(shí)，同一研究地點(diǎn)不同觀測(cè)樣方(直徑10 cm)間土壤呼吸作用的變異系數(shù)最高，達(dá)40%。Yim等［27］對(duì)落葉松人工林研究表明，在30 m×30 m樣地內(nèi)50個(gè)樣方(直徑12.6 cm)兩個(gè)觀測(cè)日期的土壤呼吸作用變異系數(shù)平均值為28%。Rodeghiero和Cescatti［28］對(duì)兩個(gè)寒溫針葉林生態(tài)系統(tǒng)研究表明，挪威云杉林和銀冷杉林觀測(cè)樣方(直徑10 cm)間土壤呼吸作用變異系數(shù)平均值分別為41.5%和28.9%。另外，土壤呼吸作用空間變異性與其觀測(cè)樣方的大小有關(guān)［27，35］。

對(duì)于相對(duì)均一的生態(tài)系統(tǒng)來說，溫度和水分通常是影響土壤呼吸作用時(shí)間動(dòng)態(tài)變化的主導(dǎo)因素，而根系、凋落物、土壤有機(jī)質(zhì)和孔隙度等反映了土壤呼吸作用的空間變化特征［27，34］。土壤呼吸作用主要由根系呼吸和微生物呼吸組成，森林生態(tài)系統(tǒng)根系呼吸所占比例平均為45%—50%［36］。細(xì)根分布不均是該區(qū)域土壤呼吸作用空間變動(dòng)的主要因素之一，能解釋土壤呼吸作用32.2%—39.6%變化。究其原因，根系呼吸速率隨直徑增加而呈明顯下降趨勢(shì)［37-38］，細(xì)根呼吸速率大于粗根。Cook等［39］研究也表明，隨著與樹干距離的增加土壤呼吸作用呈指數(shù)形式遞減，主要與細(xì)根密度的減少有關(guān)。

［1］Gower ST，Richards J H.Larches:deciduous conifers in an evergreen world.Bioscience，1990，40(11):818-826.

［2］Fang J Y，Piao SL，Zhao SQ.The carbon sink:the role of the middle and high latitudes terrestrial ecosystems in the northern hemisphere.Acta Phytoecologica Sinica，2001，25(5):594-602.

［3］Yu G R，Li H T，Wang SQ.Global Change，Carbon Cycle and Storage in Terrestrial Ecosystem.Beijing:China Meteorological Press，2003.

［4］Peng SL，Li Y L，Ren H，Zhao P.Progress in research on soil respiration under the global change.Advances in Earth Sciences，2002，17(5):705-713.

［5］IPCC.Climate change 2001:The Scientific Basis.Intergovernmental Panel on Climate Change.2001.http://www.grida.no/climate/ipcc_tar/wg1/index.htm.

［6］Russell C A，Voroney R P.Carbon dioxide efflux from the floor of a boreal aspen forest.I.Relationship to environmental variables and estimates of C respired.Canadian Journal of Soil Science，1998，78(2):301-310.

［7］Morén A S，Lindroth A.CO2exchange at the floor of a boreal forest.Agricultural and Forest Meteorology，2000，101(1):1-14.

［8］Rayment M B，Jarvis P G.Temporal and spatial variation of soil CO2efflux in a Canadian boreal forest.Soil Biology and Biochemistry，2000，32(1):35-45.

［9］Swanson R V，F(xiàn)lanagan L B.Environmental regulation of carbon dioxide exchange at the forest floor in a boreal black spruce ecosystem.Agricultural and Forest Meteorology，2001，108(3):165-181.

［10］Khomik M，Arain M A，McCaughey J H.Temporal and spatial variability of soil respiration in a boreal mixedwood forest.Agricultural and Forest Meteorology，2006，140(1/4):244-256.

［11］Karhu K，F(xiàn)ritze H，Tuomi M，Vanhala P，Spetz P，Kitunen V，Liski J.Temperature sensitivity of organic matter decomposition in two boreal forest soil profiles.Soil Biology and Biochemistry，2010，42(1):72-82.

［12］Wang Q F，Wang C K，Tan L H.Vernal soil respiration of Larix gmelinii Rupr.forests transplanted from a latitudinal transect.Acta Ecologica Sinica，2008，28(5):1883-1892.

［13］Lloyd J，Taylor JA.On the temperature dependence of soil respiration.Functional Ecology，1994，8(3):315-323.

［14］Peng SS，Piao SL，Wang T，Sun J Y，Shen Z H.Temperature sensitivity of soil respiration in different ecosystems in China.Soil Biology and Biochemistry，2009，41(5):1008-1014.

［15］Liao Y，Yang Z F，Xia X Q，Jiang H C.Research on temperature sensitivity of soil respiration and different active organic carbon fractions of Qinghai-Tibet Plateau permafrost.Earth Science Frontiers，2011，18(6):85-93.

［16］Zhou M，Yu X X，F(xiàn)eng L.Analysis on the driving force of the permafrost degeneration in Daxinganling.Journal of Arid Land Resources and Environment，2002，16(4):44-47.

［17］Sun F H，Yang X Q，Lu S，Yang S Y.The contrast analysis on the average and extremum temperature trend in northeast China.Scientia Meteorologica Sinica，2006，26(2):157-163.

［18］Wang Y，Zhou G S，Jia B R，Na JH.Climate change characteristics in Chinese boreal forest region from 1954 to 2005.Chinese Journal of Applied Ecology，2008，19(5):942-948.

［19］Jia B R，Zhou G S.Advance in the studies of the response of boreal forest to climate change.Advances in Earth Sciences，2009，24(6):668-674.

［20］Buchmann N.Biotic and abiotic factors controlling soil respiration rates in Picea abies stands.Soil Biology and Biochemistry，2000，32(11/12):1625-1635.

［21］Lee M S，Nakane K，Nakatsubo T，Koizumi H.Seasonal changes in the contribution of root respiration to total soil respiration in a cool-temperate deciduous forest.Plant and Soil，2003，255(1):311-318.

［22］Li H J，Yan J X，Yue X F，Wang M B.Significance of soil temperature and moisture for soil respiration in a Chinese mountain area.Agricultural and Forest Meteorology，2008，148(3):490-503.

［23］Yan JX，Qin Z D，Zhang Y H，Li H J.Effect of soil temperature and moisture on soil CO2efflux in a Pinus tabulaeformis forest.Acta Ecologica Sinica，2009，29(12):6366-6376.

［24］Mu C C，Cheng W，Sun X X，Wu Y X.Seasonal variation of emission fluxes of CO2，N2Oand CH4from Larix gemlinii swamps soils in Xiaoxing'an Mountains of China.Scientia Silvae Sinicae，2010，46(7):7-15.

［25］Zhang J B，Shangguan T L，Meng Z Q.Changes in soil carbon flux and carbon stock over a rotation of poplar plantations in northwest China.Ecological Research，2011，26(1):153-161.

［26］Thierron V，Laudelout H.Contribution of root respiration to total CO2efflux from the soil of a deciduous forest.Canadian Journal of Forest Research，1996，26(7):1142-1148.

［27］Yim M H，Joo SJ，Shutou K，Nakane K.Spatial variability of soil respiration in a larch plantation:estimation of the number of sampling points required.Forest Ecology and Management，2003，175(1/3):585-588.

［28］Rodeghiero M，Cescatti A.Spatial variability and optimal sampling strategy of soil respiration.Forest Ecology and Management，2008，255(1):106-112.

［29］Han G X，Zhou G S.Review of spatial and temporal variations of soil respiration and driving mechanisms.Chinese Journal of Plant Ecology，2009，33(1):197-205.

［30］Yao Y G，Zhang Y P，Yu G R，Sha L Q，Deng Y，Tan Z H.Representative time selection analysis on daily average value of soil respiration in a tropical rain forest.Journal of Nanjing Forestry University:Natural Science Edition，2011，35(4):74-78.

［31］Wang X W.Study on Key Process of Carbon Cycle in Larix Gmelinii Plantation［D］.Harbin:Northeast Forestry University，2006:12-43.

［32］Wang C K，Yang JY.Rhizospheric and heterotrophic components of soil respiration in six Chinese temperate forests.Global Change Biology，2007，13(1):123-131.

［33］Feng B P，Gao R H，Zhang Q L，Zhao Q M.Study on dynamic variation of soil microorganism under different management modes in Larix Gmelini forest.Journal of Inner Mongolia Agricultural University:Natural Science Edition，2009，30(4):74-79.

［34］Fang C，Moncrieff JB，Gholz H L，Clark K L.Soil CO2efflux and its spatial variation in a Florida slash pine plantation.Plant and Soil，1998，205(2):135-146.

［35］Davidson E A，Savage K，Verchot L V，Navarro R.Minimizing artifacts and biases in chamber-based measurements of soil respiration.Agricultural and Forest Meteorology，2002，113(1/4):21-37.

［36］Hanson P J，Edwards N T，Garten CT，Andrews JA.Separating root and soil microbial contributions to soil respiration:a review of methods and observations.Biogeochemistry，2000，48(1):115-146.

［37］Chapman SB.Some interrelationships between soil and root respiration in lowland Calluna heathland in southern England.The Journal of Ecology，1979，67(1):1-20.

［38］Jia SX，Zhao Y L，Ding G Q，Sun Y，Xu Y，Wang Z Q.Relationship among fine-root morphology，anatomy，tissue nitrogen concentration and respiration in different branch root orders in Larix gmelinii and Fraxinus mandshurica.Chinese Bulletin of Botany，2010，45(2):174-181.

［39］Cook F J，Thomas S M，Kelliher F M，Whitehead D.A model of one-dimensional steady-state carbon dioxide diffusion from soil.Ecological Modelling，1998，109(2):155-164.

參考文獻(xiàn):

［2］方精云，樸世龍，趙淑清.CO2失匯與北半球中高緯度陸地生態(tài)系統(tǒng)的碳匯.植物生態(tài)學(xué)報(bào)，2001，25(5):594-602.

［3］于貴瑞，李海濤，王紹強(qiáng).全球變化與陸地生態(tài)系統(tǒng)碳循環(huán)和碳蓄積.北京:氣象出版社，2003.

［4］彭少麟，李躍林，任海，趙平.全球變化條件下的土壤呼吸效應(yīng).地球科學(xué)進(jìn)展，2002，17(5):705-713.

［12］王慶豐，王傳寬，譚立何.移栽自不同緯度的落葉松 (Larix gmelinii Rupr.)林的春季土壤呼吸.生態(tài)學(xué)報(bào)，2008，28(5):1883-1892.

［15］廖艷，楊忠芳，夏學(xué)齊，蔣宏忱.青藏高原凍土土壤呼吸溫度敏感性和不同活性有機(jī)碳組分研究.地學(xué)前緣，2011，18(6):85-93.

［16］周梅，余新曉，馮林.大興安嶺林區(qū)多年凍土退化的驅(qū)動(dòng)力分析.干旱區(qū)資源與環(huán)境，2002，16(4):44-47.

［17］孫鳳華，楊修群，路爽，楊素英.東北地區(qū)平均、最高、最低氣溫時(shí)空變化特征及對(duì)比分析.氣象科學(xué)，2006，26(2):157-163.

［18］王宇，周廣勝，賈丙瑞，那濟(jì)海.1954—2005年中國北方針葉林分布區(qū)的氣候變化特征.應(yīng)用生態(tài)學(xué)報(bào)，2008，19(5):942-948.

［19］賈丙瑞，周廣勝.北方針葉林對(duì)氣候變化響應(yīng)的研究進(jìn)展.地球科學(xué)進(jìn)展，2009，24(6):668-674.

［23］嚴(yán)俊霞，秦作棟，張義輝，李洪建.土壤溫度和水分對(duì)油松林土壤呼吸的影響.生態(tài)學(xué)報(bào)，2009，29(12):6366-6376.

［24］牟長城，程偉，孫曉新，吳云霞.小興安嶺落葉松沼澤林土壤CO2，N2O和CH4的排放規(guī)律.林業(yè)科學(xué)，2010，46(7):7-15.

［29］韓廣軒，周廣勝.土壤呼吸作用時(shí)空動(dòng)態(tài)變化及其影響機(jī)制研究與展望.植物生態(tài)學(xué)報(bào)，2009，33(1):197-205.

［30］姚玉剛，張一平，于貴瑞，沙麗清，鄧云，譚正洪.熱帶雨林土壤呼吸測(cè)定代表性時(shí)段研究.南京林業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版，2011，35(4):74-78.

［31］王秀偉.興安落葉松人工林碳循環(huán)關(guān)鍵過程的研究［D］.哈爾濱:東北林業(yè)大學(xué)，2006:12-43.

［33］馮保平，高潤宏，張秋良，趙全民.不同經(jīng)營方式下興安落葉松林土壤微生物年季動(dòng)態(tài)研究.內(nèi)蒙古農(nóng)業(yè)大學(xué)學(xué)報(bào):自然科學(xué)版，2009，30(4):74-79.

［38］賈淑霞，趙妍麗，丁國泉，孫玥，許旸，王政權(quán).落葉松和水曲柳不同根序細(xì)根形態(tài)結(jié)構(gòu)、組織氮濃度與根呼吸的關(guān)系.植物學(xué)報(bào)，2010，45(2):174-181.