重構土體土壤呼吸及其水熱影響因子日變化特征

張 劍，雷 娜

(1.長安大學環境科學與工程學院，陜西西安710064;2.陜西省土地工程建設集團有限責任公司，陜西西安710075;3.西北農林科技大學水土保持研究所，陜西楊凌712100)

在氣候變暖日益加劇的大背景下，土壤呼吸通過參與生態系統碳循環影響氣候變化，因此，國內外眾多學者開展了針對不同類型生態系統、不同土地利用類型［1-3］、不同地域、退化土壤等土壤呼吸特征的研究［4-5］。學者們的研究一般都是從最小時間研究尺度——土壤呼吸日變化開始的。土壤呼吸日變化規律的研究可為進一步掌握土壤呼吸季節變化及其影響因子、準確核算區域碳排放情況提供數據支撐。眾多學者研究結果均顯示溫度和水分是土壤呼吸日變化的重要環境因子［6-7］，有許多研究表明土壤呼吸速率與大氣溫度或土壤不同深度溫度之間具有顯著的相關關系，一般可用線性函數、指數函數等形式來描述［8-9］。

我國耕地資源緊缺，目前主要通過土地整治工程進行補充，隨著耕地后備資源的減少，土地整治已從簡單的土地平整、配套設施修建發展到土體重構。在陜西省寶雞市、華陰市等地存在大量的荒石灘地，已有工程通過添加不同成土材料和改良材料在荒石灘地進行土體重構，重構土體將成為今后土地整治的主要方向以及補充耕地資源的重要手段，與此同時土體重構改變了區域下墊面條件，導致新形成的土體土壤呼吸產生變化，進而對區域氣候產生影響，但現階段研究重構土體土壤呼吸還不多見。本研究以添加蛭石、頁巖、沙和砒砂巖的4種重構土體為研究對象，每月份選取1個典型日，使用土壤碳通量測定系統對土壤呼吸進行觀測，同時對土壤水分、溫度等水熱影響因子進行測定，以期探明4種重構土體土壤呼吸和水熱影響因子日變化特征及其之間的關系，為重構土體土壤呼吸的研究提供科學參考，同時為準確評估區域CO2排放及制定合理的CO2減排措施提供數據支撐。

1 材料與方法

1.1 試驗樣地概況

試驗樣地位于陜西省寶雞市眉縣湯峪鎮上王村(107°53＇50″E，34°8＇33＇N)。該試驗樣地是荒石灘地整治工程的示范區域，總面積8.00 hm2，實現新增耕地6.80 hm2，新增產值11.20萬元。選取砒砂巖、沙、泥質頁巖、蛭石4種材料，將選取的4種材料粉碎過10 mm篩，消毒滅菌處理后，分別與構建土源混合(構建土體所使用的土源為當地常見土壤類型塿土)，形成蛭石、頁巖、沙和砒砂巖與土壤混合層(50 cm)，進而形成礫石+蛭石+塿土型、礫石+頁巖+塿土型、礫石+沙+塿土型和礫石+砒砂巖+塿土型4種重構土體(以下簡稱蛭石、頁巖、沙、砒砂巖重構土體)，容重分別為 1.55，1.36，1.52 g·cm-3和1.38 g·cm-3，進行長期定位試驗。其中蛭石、頁巖、砒砂巖、沙用量均為1×10-3m3·m-2，添加比例均為2%(體積比)。4個試驗樣地規格均為20 m×30 m。4個試驗樣地內均埋設3個內徑為10 cm×20 cm的土壤呼吸環，保證呼吸環高出地面2 cm。

1.2 研究方法

2017年11月至2018年10月，每月選取1個典型日對4個試驗樣地全部土壤呼吸環進行測量，每個典型日的測量時間為 6∶00、8∶00、10∶00、12 ∶00、14 ∶00、16 ∶00、18 ∶00、20 ∶00。土壤呼吸測定使用土壤碳通量測量系統(LI-8100)，測定土壤碳通量、土壤5 m處溫度和10 cm處含水量，每個土壤呼吸環重復測定3次，測量時長為4 min，降水量等主要通過試驗小區3 m外安裝的HOBO氣象站自動獲取。

1.3 數據分析方法

土壤呼吸與土壤溫度之間的關系用指數模型擬合，土壤呼吸與含水量之間的關系采用二次曲線模型擬合，土壤呼吸與土壤溫度和土壤體積含水量之間的關系采用冪-指數模型擬合。

式中，RS為土壤呼吸速率(μmol·m-2·s-1)，T為土壤溫度(℃)，W為土壤體積含水量(%);a、b、c分別為擬合參數。

2 結果與分析

2.1 重構土體土壤呼吸日變化

4種重構土體土壤呼吸在觀測月份的日變化均呈現出單峰曲線，觀測月份的最低值出現在6∶00或20 ∶00，最高值出現在10 ∶00、12 ∶00或14 ∶00。這主要是由4種重構材料蛭石、頁巖、沙、砒砂巖的物質特性決定的。在觀測的12個月中，添加蛭石、頁巖、沙、砒砂巖的4種重構土體日變化幅度的范圍分別為 0.37～3.42、0.52～4.34、0.71～4.62、0.55～4.03μmol·m-2·s-1。觀測12個月內4種重構土體日變化幅度的排序為添加沙＞頁巖＞砒砂巖＞蛭石，可見添加沙的重構土體土壤呼吸日變化幅度最大，頁巖次之，砒砂巖較小，蛭石日變化幅度最小(圖1)。

1月份，氣溫達到一年中的最低值，作物和微生物生理活動基本停止，4種重構土體土壤呼吸值在觀測時段內達到最低，日變化幅度也最小。添加蛭石、頁巖、沙、砒砂巖的重構土體土壤呼吸最小值分別為 0.04、0.06、0.04 μmol·m-2·s-1和 0 μmol·m-2·s-1，最大值均出現在12∶00，其日變化幅度分別為 0.37、0.59、0.71 μmol·m-2·s-1和 0.55 μmol·m-2·s-1(圖 1c)。

8月份，屬于作物的生長旺季，氣溫高且雨量充沛，4種重構土體土壤呼吸值達到試驗觀測期間最大值，觀測時段最大值均出現在12∶00，最小值出現在6∶00和20∶00。添加蛭石、頁巖、沙、砒砂巖的重構土體土壤呼吸值最高分別為9.88、10.47、14.94 μmol·m-2·s-1和 12.19 μmol·m-2·s-1，日變化幅度分別為 3.42、4.34、4.62 μmol·m-2·s-1和4.03 μmol·m-2·s-1(圖 1j)。

2.2 重構土體土壤呼吸與水熱影響因子日變化的關系

2.2.1 土壤呼吸與土壤溫度日變化關系 指數模型能夠表征添加蛭石、頁巖、沙、砒砂巖4種重構土體土壤呼吸與土壤5 cm處溫度之間的日變化關系。4種重構土體土壤呼吸速率與土壤溫度存在顯著的正指數相關關系。觀測時段內添加蛭石的重構土體土壤溫度能夠解釋土壤呼吸日變化的41.0%～96.0%，添加頁巖的重構土體土壤溫度能夠解釋土壤呼吸日變化的53.0%～98.0%，添加沙的重構土體土壤溫度能夠解釋土壤呼吸日變化的76.0%～98.0%，添加砒砂巖的重構土體土壤溫度能夠解釋土壤呼吸日變化的66.0%～96.0%(表1)。其中添加沙重構土體土壤溫度對土壤呼吸日變化的解釋能力最高，添加蛭石的重構土體土壤溫度對土壤呼吸日變化的解釋能力最低。

2.2.2 土壤呼吸與土壤體積含水量日變化關系

二次曲線模型能夠解釋4種重構土體土壤呼吸與土壤體積含水量之間日變化的關系，二者之間呈顯著負相關，但是解釋能力顯著低于溫度的解釋能力。觀測時段內添加蛭石的重構土體土壤體積含水量單因子能夠解釋土壤呼吸日變化的23.0%～46.0%，添加頁巖的重構土體土壤體積含水量單因子能夠解釋土壤呼吸日變化的28.0%～53.0%，添加沙的重構土體土壤體積含水量單因子能夠解釋土壤呼吸日變化的26.0%～44.0%，添加砒砂巖的重構土體土壤體積含水量單因子能夠解釋土壤呼吸日變化的24.0%～50.0%(表2)。添加頁巖重構土體土壤體積含水量單因子對土壤呼吸日變化的解釋能力最高，添加蛭石重構土體土壤體積含水量單因子對土壤呼吸日變化的解釋能力最低。

2.2.3 土壤呼吸與土壤溫度、體積含水量日變化關系 冪-指數模型能夠表征4種重構土體土壤呼吸與土壤溫度、體積含水量之間的日變化關系，土壤呼吸與水熱因子之間呈顯著正相關。

觀測時段內添加蛭石的重構土體土壤溫度和體積含水量雙因子能夠解釋土壤呼吸日變化的56.0%～94.0%，添加頁巖的重構土體土壤溫度和體積含水量雙因子能夠解釋土壤呼吸日變化的67.0%～53.0%，添加沙的重構土體土壤溫度和體積含水量雙因子能夠解釋土壤呼吸日變化的83.0%～98.0%，添加砒砂巖的重構土體土壤溫度和土壤體積含水量雙因子能夠解釋土壤呼吸日變化的56.0%～99.0%(表3)。其中添加沙的重構土體土壤溫度和土壤體積含水量雙因子對土壤呼吸日變化的解釋能力最高，添加蛭石的重構土體土壤溫度和土壤體積含水量雙因子對土壤呼吸日變化的解釋能力最低。

圖1 測定期內各月份不同重構土體土壤呼吸日變化Fig.1 Diurnal variation of soil respiration in different reconstructed soils in each month

表1 重構土體土壤呼吸速率(RS)與土壤溫度(Ts)日變化之間的關系(RS=a e bT)Table 1 Relationship betweendaily soil respiration rate and temperature for reconstituted soil mass

表2 重構土體呼吸速率(RS)與土壤體積含水量(WS)日變化之間的關系(RS=aw2+bw+c)Table 2 Relationship between daily soil respiration rate and volumetric water content for reconstituted soil mass

表3 重構土體呼吸速率(RS)與土壤溫度(Ts)、體積含水量(WS)日變化的關系(RS=a e bTwc)Table 3 Relationship between daily soil respiration rate and temperature and volumetric water content for reconstituted soil mass

3 討論

3.1 不同重構材料對土壤呼吸的影響

國內外已研究了高達60種礦物材料應用于土壤質量提升，并形成系列改良產品如礦物肥料、生長劑、營養物質載體等［10］。現階段眾多學者研究認為砒砂巖能夠改良土壤性狀，增加土壤的持水性和保水性，吸附土壤重金屬污染，進而提高作物產量［11］。學者們研究表明，土壤中適當摻沙可增加土壤的通透性、滲水性、保濕性，使土壤中空氣含量增加［12］。康倍銘等［13］研究了添加一定比例的頁巖對土壤水穩性團聚體、孔隙度、土壤有機質、有效含水量和水分蒸發等理化性質的影響，結果表明改良后的土壤孔隙度、有效含水量和水分蒸發變化不明顯，但土壤有機質、水穩性團聚體顯著提高。蛭石作為土壤改良劑被學者們廣泛認可，研究表明蛭石能夠改善土壤的理化性質及生物營養成分，并在物理性質方面促使土壤形成良好的結構，提高蓄水保墑能力。

不同重構材料的性質不同，會導致土壤呼吸表現出一定的差異性。砒砂巖、頁巖黏粒含量高，是成土材料，蛭石作為改良材料，3種材料均可增加土壤的保水性和持水性。水分占據了土壤的空隙，導致土壤中CO2含量相對較少，土壤呼吸不強烈，有利于環境保護。同時沙作為土壤改良材料，能夠使土壤中空氣含量增加從而增加土壤呼吸強度，不利于環境保護。

3.2 重構土體土壤呼吸與水熱影響因子之間的關系

溫度是影響土壤呼吸速率的重要因子。學者們認為指數模型能夠解釋土壤呼吸速率的絕大部分變異［14-16］，土壤溫度可通過影響植物根系以及土壤中微生物的活動間接影響土壤呼吸。本研究發現4種重構土體呼吸日變化與土壤溫度呈顯著正相關，與大部分學者得到的結論一致［17-19］。土壤體積含水量對土壤呼吸的影響比較復雜，學者們在不同地域、不同研究時段、不同生態系統等得出的結論往往差異較大［20-21］。土壤體積含水量除能夠通過影響根系的生長以及土壤中微生物的活動影響土壤呼吸速率外，還可以通過影響土壤中CO2的含量以及傳輸過程影響土壤呼吸速率。有研究指出土壤體積含水量與土壤之間的關系可用對數模型、二次曲線和線性關系模擬，二者呈顯著負相關［22-23］;也有研究指出，當土壤發生水分脅迫時，水分可成為影響土壤呼吸的主導因素［24］。本研究發現二次曲線的模擬效果好于線性模型，但體積含水量與土壤呼吸的相關系數較小，即對土壤呼吸的影響較小。多數研究均表明土壤呼吸日變化一般受到土壤溫度和土壤水分的共同影響，因而采用二者復合來表征土壤呼吸的日變化情況［25-26］。有學者對中國溫帶森林土壤呼吸的6種生態系統進行研究，發現雙因素模型的解釋能力一般高于50.0%，解釋能力通常情況下高于單因素模型［27］。本研究中得到冪-指數模型即土壤溫度、體積含水量對土壤呼吸日變化的解釋能力大于56.0%，其和指數模型的解釋能力遠高于二次曲線模型，但大部分情況下，冪-指數模型的解釋能力高于指數模型，這與前人得到的結論類似。

4 結 論

1)4種重構土體土壤呼吸在觀測月份的日變化均呈現單峰曲線，呼吸日變化模式基本相同，觀測月份的最低值出現在6∶00和20∶00，最高值出現在10∶00、12∶00和14∶00。8月份，4種重構土體土壤呼吸值達到試驗觀測期間最大值，12月份4種重構土體土壤呼吸值在觀測時段內達到最低，日變化幅度也最小，觀測12個月份內4種重構土體日變化幅度的排序為添加沙＞頁巖＞砒砂巖＞蛭石，添加蛭石的重構土體更有利于生態環境的保護。

2)4種重構土體土壤溫度(指數模型)、土壤體積含水量(二次曲線模型)單因素以及土壤溫度和含水量雙因素(冪-指數模型)均能夠解釋土壤呼吸的日變化，但是解釋能力不同。雙因素模型解釋能力一般高于50.0%，單因素模型中溫度的解釋能力顯著高于水分，水分的解釋能力不強，一般低于50.0%，表征重構土體水熱影響因子與呼吸之間的關系最優的是雙因素模型。