來亞男，嚴(yán)俊霞，李洪建

（山西大學(xué)黃土高原研究所，山西 太原 030006）

土壤呼吸（Rs）是將植物固定的CO2釋放到大氣的主要途徑。據(jù)估計，在全球范圍內(nèi)，Rs每年釋放的CO2量是工業(yè)釋放CO2的10倍[1-2]。由于土壤巨大的呼吸量，因此，對不同生態(tài)系統(tǒng)的Rs研究越來越受到關(guān)注[3-5]。在土壤溫度作用的驅(qū)使下，Rs具有明顯的日變化特點(diǎn)，不同地區(qū)Rs的日變化幅度不同，主要取決于土壤溫度的變化范圍和幅度，這給精確估算日均Rs帶來很大的困難。受儀器和勞動力的限制，目前絕大多數(shù)土壤呼吸測定很難做到24 h的連續(xù)測定，而是以不同時間的測定結(jié)果代表日土壤呼吸值。有研究表明，用白天測定的Rs代表日均Rs可能高估Rs[6]。他們認(rèn)為白天12 h的Rs平均值比24 h的Rs平均值大10%。也有研究報道，可用12 h的平均值代替24 h的平均值[7]；夜間Rs的平均值比白天的小8%[8]。因?yàn)榇蠖鄶?shù)Rs測定都在白天進(jìn)行且不連續(xù)。因此，白天不同時段所測得的Rs的代表程度如何，是許多研究所關(guān)注的問題。

本試驗(yàn)對太原盆地3種土地利用方式的土壤呼吸、土壤溫度進(jìn)行了3次24 h的定位測定，旨在研究和比較不同土地利用方式下不同時間段的Rs日動態(tài)及其差異；分析土壤溫度對Rs日變化的影響，明確估算日均Rs的最佳時間窗口。

1 材料和方法

1.1 試驗(yàn)地概況

試驗(yàn)在山西省農(nóng)業(yè)科學(xué)院農(nóng)業(yè)環(huán)境與資源研究所北營試驗(yàn)基地進(jìn)行。試驗(yàn)所在區(qū)域海拔785m。氣候?yàn)榕瘻貛О敫珊荡箨懶约撅L(fēng)氣候，年平均氣溫9.5℃；最冷月1月平均氣溫-6.6℃，最熱月7月平均氣溫23.5℃。多年平均降水量478 mm，降水年內(nèi)分布極為不均，夏季降水量占全年降水量的60%、秋季23%、春季14%和冬季3%。試驗(yàn)田土壤均為沖積黃土，土層深厚，具有廣泛的代表性。3種土地利用類型分別種植檸條、中藥材和玉米。樣地土壤理化性狀：土壤容重1.38～1.381g/cm3；田間持水量21.7%～22.7%；有機(jī)質(zhì)3.91%～3.99%；全氮0.11%；速效磷6.80～12.292 mg/kg；速效鉀 23.0～277.7 mg/kg。

1.2 土壤呼吸及環(huán)境因子測定

用LI-COR 6400便攜式光合作用儀連接LI-COR-09土壤氣室測定Rs。3種土地利用方式Rs測定均為3個樣點(diǎn)。首先在每個樣地選取一個區(qū)域，面積約15 m2，每次測定都選定區(qū)域進(jìn)行。測定樣點(diǎn)呈三角形分布，每個樣點(diǎn)測定1次，每次3個循環(huán)，9個數(shù)據(jù)，取其平均值作為該次測定值進(jìn)行分析。測定于6月22日、8月4日和10月10日進(jìn)行，6月22日共測定9次，其余2 d為11次。10 cm深度的土壤溫度（T10）由系統(tǒng)土壤溫度探針（6400-12）測定。在測定Rs的同時，將溫度探針垂直插入所測土體附近10 cm深，系統(tǒng)將自動同步記錄T10和Rs數(shù)據(jù)，用來分析Rs與T10的關(guān)系。采用傳統(tǒng)烘干法測定0～10 cm深度的土壤水分。

1.3 數(shù)據(jù)分析

簡單方差分析用于分析3個樣地之間的Rs，T10和土壤水分（Ws，%）均值差異程度。日尺度土壤呼吸和溫度的變異系數(shù)CV（%）用于分析同一樣地內(nèi)各指標(biāo)的日變異程度。用指數(shù)模型（Rs=αeβT；式中，Rs為土壤呼吸，T為實(shí)測土壤溫度；α，β為擬合參數(shù)）分析土壤呼吸（Rs）和土壤溫度（T10）的因子關(guān)系。用偏差（deviation）來衡量各時間段Rs的平均值（Rr）與日平均Rs（Rm）的差異大小，其計算公式為deviation＝（Rr-Rm）/Rm×100%。

數(shù)據(jù)處理用Microsoft Excel 2007和SPSS 22.0軟件，用Sigmaplot 12.5繪制圖形。

2 結(jié)果與分析

2.1 土壤呼吸的日變化趨勢

3次Rs日變化測定時的T10有比較明顯的日變化特點(diǎn)。早晨，土壤溫度較低，隨氣溫升高，土壤溫度亦逐漸升高，土壤溫度的最高值出現(xiàn)在15：00—16：00，最低值出現(xiàn)在 6：00—8：00。Rs同樣具有明顯的日變化，3種土地利用方式Rs的日變化在3次測定時均呈單峰趨勢，且不對稱（圖1）。在6月份，檸條地、藥材地和玉米地的Rs日變化趨勢比較一致，最大值均出現(xiàn)在12：00左右，最低值均出現(xiàn)在6：00左右。藥材地Rs的CV最高，為25.38%，檸條地和玉米地的CV差異不大，在14%左右。在8月份，3種土地利用方式下的Rs最小值出現(xiàn)在晚上或者早晨溫度相對較低的時候，檸條地和藥材地的Rs最大值出現(xiàn)在 12：00—14：00，而玉米地的 Rs在21：00左右達(dá)到最大，主要由于玉米地8月3日降水之前的灌溉，測定時的土壤水分相對其他樣地較高，水分過多抑制Rs，隨后經(jīng)過一天的水分蒸發(fā)，水分降低到抑制Rs的臨界值之下，最終因?yàn)橥砩蠝囟戎饾u下降，Rs也隨之逐漸下降。檸條地和藥材地Rs的CV比玉米地的分別高12.68%和12.10%。在10月份，3種土地利用方式下的土壤呼吸最大值均出現(xiàn)在 14：00—16：00，最小值均出現(xiàn)在 6：00，土壤呼吸速率日變化趨勢比較一致。檸條地Rs的CV最低，為16.35%，藥材地和玉米地Rs的CV差異不大，在23%左右。

2.2 土壤呼吸的日變化差異

3種土地利用類型的Rs日變化趨勢相一致，但Rs值有所不同，3個測定日3種土地利用方式的土壤呼吸差異性檢驗(yàn)結(jié)果如表1所示。6月份藥材地的T10顯著低于檸條地和玉米地的T10，而藥材地的Rs卻顯著高于檸條地和玉米地的Rs。其原因可能是因?yàn)樗幉牡氐腤s顯著高于檸條地和玉米地的Ws，檸條地和玉米地的Rs受到土壤水分干旱脅迫所致。8月份3種土地利用方式下的Rs總體差異不顯著（P＞0.05），玉米地的Rs略高于檸條地和藥材地的Rs。3種土地利用方式下的T10和Ws總體差異顯著（P＜0.05），藥材地的T10顯著高于玉米地的T10，而與檸條地的T10無顯著差異，玉米地的Ws顯著高于藥材地和檸條地的Ws。10月份Rs的日均值為玉米地＞檸條地＞藥材地，其中，檸條地和玉米地的Rs顯著高于藥材地的Rs（P＜0.05）。3種土地利用方式下的T10相互之間無顯著差異。

從3次測定的結(jié)果來看，3種不同土地利用方式下的土壤呼吸呈現(xiàn)先增加后減低的趨勢。3個樣地的Rs日均值最高值出現(xiàn)在8月，而檸條地和玉米地的T10的最高值出現(xiàn)在6月，藥材地則出現(xiàn)在8月，檸條地、藥材地和玉米地的Ws的最高值分別出現(xiàn)在10，6，8月，土壤呼吸的最高值并未出現(xiàn)在土壤溫度或土壤水分最高的月份，而是出現(xiàn)在土壤溫度和土壤水分適中的月份，表明土壤呼吸受土壤溫度和土壤水分的共同影響。

表1 土壤呼吸、土壤溫度及土壤水分的日變化方差分析

2.3 土壤呼吸日變化與土壤溫度（T10）日變化的 關(guān)系

表2 3種土地利用方式的土壤呼吸日變化與土壤溫度日變化的關(guān)系

分別用各次測定的T10為自變量、Rs為因變量擬合關(guān)系模型，結(jié)果列于表2。由表2可知，3個測定日除了8月測定的玉米地的T10與Rs相關(guān)關(guān)系顯著外（P＜0.05），其余測定日3種土地利用類型的T10與Rs相關(guān)關(guān)系都不顯著（P＞0.05）。分析其原因可能與土壤溫度的變化滯后于土壤呼吸的日變化有關(guān)，如土壤呼吸的最大值多出現(xiàn)在12：00左右，而土壤溫度的最大值則出現(xiàn)在14：00—15：00時。當(dāng)考慮了溫度的滯后效應(yīng)后重新進(jìn)行分析時，發(fā)現(xiàn)T10與Rs的關(guān)系方程的決定系數(shù)都有很大的提高，進(jìn)一步表明了土壤溫度的日變化比土壤呼吸的日變化有滯后性。

2.4 不同時間段土壤呼吸值與土壤呼吸日均值的比較

由表3可知，3個測定日檸條地和藥材地的Rs日均值與 8：00—10：00（8 月除外）和 16：00—18：00這2個時間段的平均值最為接近，偏差在±10%左右；3個測定日玉米地全天各時間段的Rs平均值與日均值的偏差變異不大，6月 8：00—10：00和14：00—18：00，8 月和 10 月的 8：00—16：00 各時間段的均值都能用于估算Rs日均值，偏差在±15%以內(nèi)。

表3 白天不同時段土壤呼吸的平均值與24 h土壤呼吸平均值的比值分析 %

3 結(jié)論與討論

在溫帶森林或其他土地利用條件下，受氣溫日變化的影響，土壤呼吸常常表現(xiàn)出明顯的日變化特點(diǎn)[9-10]，即使在熱帶地區(qū)也是如此[11]。YU等[12]研究表明，松嫩平原西部草甸草原的6個樣地的土壤呼吸均具有明顯的日變化特點(diǎn)，最低值出現(xiàn)在21：00—23：00 或黎明前，最高值出現(xiàn)在 11：00—13：00。本研究中，3種土地利用方式的土壤呼吸均有明顯的日變化趨勢，最低值出現(xiàn)在23：00或黎明前后，最高值出現(xiàn)在12：00—16：00。不同土地利用方式不同季節(jié)最高值和最低值出現(xiàn)的時間不盡相同，可能的原因是研究區(qū)不同時間土壤微生物活性和根系呼吸對溫度和水分的敏感性不同所致。

一般情況下，一天之內(nèi)土壤水分及土壤理化性質(zhì)都不會有明顯變化。因此，土壤呼吸的日變化主要受土壤溫度日變化的影響[13]。本研究中，Rs隨土壤溫度的升高增加，但相關(guān)關(guān)系均不顯著；但是去除土壤溫度的滯后效應(yīng)后，相關(guān)系數(shù)得到了極大的提高，表明土壤溫度對Rs的影響具有滯后性，其滯后程度與土壤水分有關(guān)[14]。Rs的最高值并未出現(xiàn)在土壤溫度或土壤水分最高的月份，而是出現(xiàn)在土壤溫度和土壤水分適中的月份，表明土壤呼吸受土壤溫度和土壤水分的共同影響[15]。

日尺度土壤呼吸值是精確估算區(qū)域累計CO2釋放量的一個重要指標(biāo)[16]。許多研究把在單一時間點(diǎn)觀測的土壤呼吸當(dāng)作測定當(dāng)天甚至測定當(dāng)周或2周的平均土壤呼吸，這種方法會高估或低估土壤呼吸[17-18]。有學(xué)者認(rèn)為，9：00—11：00 的土壤呼吸值可以代表日平均值，其誤差在0.9%～1.5%[19]。YU等[12]研究結(jié)果表明，松嫩平原西部草甸草原的6個樣地的 7：00—9：00 和 17：00—19：00 間測定的土壤呼吸平均值能夠代表24 h土壤呼吸的日平均值；DAVIDSON等[16]利用溫帶闊葉混交林連續(xù)2 d觀測的數(shù)據(jù)研究表明，9：00—12：00測定的值能夠代表日平均值[20]。本研究中，3個測定日期3種土地利用類型的 Rs日均值與 8：00—10：00 和 16：00—18：00這2個時間段的平均值最為接近，偏差在±10%左右，用該時段的值可作為日均值計算呼吸通量值。

［1］RAICH J W，POTTER C S.Global patterns of carbon dioxide emissions from soils[J].Global Biogeochemical Cycles，1995，9（1）：23-36.

［2］ SCHLESINGER W H，ANDREWS J A.Soil respiration and the global carbon cycle[J].Biogeochemistry，2000，48（1）：7-20.

［3］WANG C，MA Y，TROGISCH S，et al.Soil respiration is driven by fine root biomass alonga forest chronosequence in subtropical China[J].Journal of Plant Ecology，2017，10（1）：36-46.

［4］PHILLIPS C L，BOND-LAMBERTY B，DESAI A R，et al.Erratum to:The value of soil respiration measurements for interpreting and modeling terrestrial carbon cycling[J].Plant&Soil，2016，1（1）：1-25.

［5］ LIU X，ZHENG J，ZHANG D，et al.Biochar has no effect on soil respiration across Chinese agricultural soils[J].Science of the Total Environment，2016，554/555：259-265.

［6］QI Y，XU M.Separating the effects of moisture and temperature on soil CO2efflux in a coniferous forest in the Sierra Nevada mountains[J].Plant&Soil，2001，237（1）：15-23.

［7］LAW B E，KELLIHER F M，BALDOCCHI D D，et al.Spatial and temporal variation in respiration in a young ponderosa pine forest during a summer drought[J].Agricultural&Forest Meteorology，2001，110（1）：27-43.

［8］ FLANAGAN L B，JOHNSON B G.Interacting effects of temperature，soil moisture and plant biomass production on ecosystemrespiration in a northern temperate grassland[J].Agricultural&Forest Meteorology，2005，130（3/4）：237-253.

［9］XU M，QI Y.Soil-surface CO2efflux and its spatial and temporal variations in a young ponderosa pine plantation in northern California[J].Global Change Biology，2001a，7（6）：667-677.

［10］楊晶，黃建輝，詹學(xué)明，等.農(nóng)牧交錯區(qū)不同植物群落土壤呼吸的日動態(tài)觀測與測定方法比較 [J].植物生態(tài)學(xué)報，2004，28（3）：318-325.

［11］ ELENEIDE D S，PATRICK M，YADVINDER M，et al.Soil CO2effluxin a tropical forest in the central Amazon[J].Global Change Biology，2010，10（5）：601-617.

［12］YU X L，LI X J，XU L S，et al.Diurnal variation of soil CO2efflux and its optimal measuring time-window of temperate meadow steppes in Western Songnen Plain，China[J].Chinese Geographical Science，2016，26（4）：518-526.

［13］李洪建.不同生態(tài)系統(tǒng)土壤呼吸與環(huán)境因子的關(guān)系研究[D].太原：山西大學(xué)，2008.

［14］RIVEROS IREGUI D A，EMANUEL R E，MUTH D J，et al.Diurnal hysteresis between soil CO2and soil temperature is controlled by soil water content[J].Geophysical Research Letters，2007，34（17）：138.

［15］劉惠，趙平.華南丘陵區(qū)典型土地利用類型土壤呼吸日變化[J].生態(tài)環(huán)境學(xué)報，2008，17（1）：249-255.

［16］嚴(yán)俊霞，李洪建，湯億，等.小尺度范圍內(nèi)植被類型對土壤呼吸的影響[J].環(huán)境科學(xué)，2009，30（11）：3121-3129.

［17］荊雪鍇，嚴(yán)俊霞，李洪建.不同土地利用下的土壤呼吸及其與環(huán)境因子的關(guān)系[J].山西農(nóng)業(yè)科學(xué)，2016，44（8）：1151-1157.

［18］嚴(yán)俊霞，秦作棟，張義輝，等.土壤溫度和水分對油松林土壤呼吸的影響[J].生態(tài)學(xué)報，2009，29（12）：6366-6376.

［19］XU M，QI Y.Spatial and seasonal variations of Q10determined by soil respiration measurements at a Sierra Nevadan Forest[J].Global Biogeochemical Cycles，2001，15（3）：199-209.

［20］DAVIDSON E A，SAVAGE K，VERCHOT L V，et al.Minimizing artifacts and biases in chamber-based measurements of soil respiration [J].Agricultural&Forest Meteorology，2002，113（1/4）：21-37.