楊 紅，劉合滿，曹麗花，連玉珍

(1. 西藏農(nóng)牧學(xué)院，資源與環(huán)境學(xué)院，西藏 林芝 860000；2.信陽農(nóng)林學(xué)院，河南 信陽 464000)

【研究意義】近百年來，以氣候變暖為主要標(biāo)志的全球氣候變化成為影響陸地生態(tài)系統(tǒng)和環(huán)境健康的重要因素，已引起國內(nèi)外學(xué)者的廣泛重視。大氣CO2濃度自19世紀(jì)以來，呈快速增加趨勢，大氣CO2濃度由1832年的284.3 μl/L增加到2014年的397.7 μl/L，僅用了182年，其年增加率高達(dá)0.6 μl/L，且自1988年開始，大氣CO2水平已超過IPCC所認(rèn)定的350 μl/L的安全水平。廣大國際學(xué)者認(rèn)為，全球氣候變暖及異常溫度事件的發(fā)生主要是因?yàn)榇髿釩O2濃度升高引起的，而全球溫度的增加與大氣CO2濃度的增加是同步的。陸地生態(tài)系統(tǒng)作為非常重要的碳庫，其在碳循環(huán)過程中的碳排放(CO2emission)是大氣CO2濃度增加的一個(gè)重要碳源，其微小的變化將對大氣CO2濃度產(chǎn)生明顯的擾動效應(yīng)，從而影響氣候變化。其中土壤碳庫作為陸地生態(tài)系統(tǒng)最重要的有機(jī)碳庫(Organic carbon pool)，據(jù)估計(jì)，全球1 m深土壤碳庫容量高達(dá)1500 Pg C(1 Pg = 1012g)[1]，分別是大氣和生物碳庫的2和3倍[2]，土壤碳庫微小的改變將可能引起大氣CO2濃度的顯著變化，故土壤碳循環(huán)及對環(huán)境變化的響應(yīng)關(guān)系將成為國內(nèi)外學(xué)者研究的焦點(diǎn)之一。森林下層枯落物在地上地下的養(yǎng)分分配、固持養(yǎng)分、涵養(yǎng)水源、降低水汽蒸發(fā)、保護(hù)土壤結(jié)構(gòu)穩(wěn)定性及維持土壤生物食物來源等方面起著重要作用。因此，森林下層枯落物的變化將極大地改變土壤理化形狀及生物群落結(jié)構(gòu)，從而在一定程度上影響土壤碳排放。【前人研究進(jìn)展】劉毓敏等[3]研究表明，相對于添加枯落層而言，去除枯落層導(dǎo)致土壤總碳含量和土壤總氮含量分別減少11.5 %和15.7 %，土壤有機(jī)碳含量下降12.76 %，微生物數(shù)量下降了9.14 %。覃志偉等[4]研究表明，相對于去除枯落層而言，添加枯落層顯著增加了34.02 %的土壤呼吸速率，年累積碳排放量增加了33.54 %。目前，國內(nèi)外對添加和移除枯落層對土壤碳排放影響的研究工作不少，但主要集中在溫?zé)釒^(qū)域，對青藏高原高寒氣候條件下，特別是過熟森林下層添加和去除枯落層對土壤呼吸影響的研究略顯不足，尚需進(jìn)一步研究。【本研究切入點(diǎn)】青藏高原平均海拔為4000 m以上，加之其自然歷史發(fā)育年輕，地形復(fù)雜、面積廣闊以及多樣的氣候類型，使青藏高原成為全球氣候變化的敏感區(qū)，在亞洲乃至全球氣候變化中均具有非常重要的作用。因此，在全球氣候變化的大背景下，青藏高原土壤碳庫將積極響應(yīng)于全球氣候變化，對青藏高原高寒森林生態(tài)系統(tǒng)表層碳庫穩(wěn)定性產(chǎn)生重要影響，從而加速大氣溫室氣體濃度升高，增加全球氣候變化壓力。【擬解決的關(guān)鍵問題】故基于西藏高原高寒氣候條件，森林生態(tài)系統(tǒng)枯落物分解速率較慢，以及枯落層土壤是下層土壤有機(jī)碳(Soil organic carbon, SOC)重要來源的科學(xué)背景下，為了定量評價(jià)不同森林土壤碳排放及枯落層碳排放對整個(gè)土壤碳排放的貢獻(xiàn)，本研究選擇藏東南色季拉山2種主要類型森林為研究對象，研究移除和保留枯落層后的土壤碳排放變化特征。旨在為科學(xué)闡明西藏高原不同類型森林枯落層對土壤碳排放的貢獻(xiàn)及未來可能的變化提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

試驗(yàn)地位于藏東南色季拉山，針葉林(N：29.654037°，E：94.373345°，A：3129 m)主要以高山松為建群種，林下多荀子、薔薇等小灌木和野草莓、禾草等草本植物；闊葉林(N：29.653723°，E：94.368423°，A：3121 m)以高山櫟為建群種，林下多蕨類等植物。

1.2 樣地選擇及研究方法

樣地選擇：在同一海拔高度選擇坡度、朝向一致的針葉林(針葉)和闊葉林(闊葉)2種典型林型，分別在每個(gè)林型選擇有代表性的1個(gè)樣點(diǎn)。

保留和去除枯落層：在每個(gè)樣點(diǎn)設(shè)置2個(gè)處理，即保留和去除枯落層。保留枯落層：本樣點(diǎn)針葉林枯落層厚度為12.0 cm，闊葉林枯落層厚度為6.5 cm，為保持針闊葉林枯落層厚度一致，在闊葉林枯落層上添加5.5 cm厚度的闊葉林枯落物；去除枯落層：在距離添加枯落層處理約25 cm處，去除地表枯落層。在每個(gè)處理上設(shè)置2個(gè)直徑20 cm，高10 cm的土壤環(huán)，其中土壤環(huán)插土壤深度約4 cm，土壤環(huán)中的植物齊地面剪除，為減小人為干擾，提前24 h在各處理上設(shè)置土壤環(huán)。

土壤呼吸(Soil respiration,Rs)測定：使用LI-8100A(Li-COR, Lincoln, NE, USA)開路式土壤碳通量測量系統(tǒng)測定針葉林和闊葉林的Rs速率。針葉林和闊葉林各用一臺LI-8100開路式土壤碳通量測量系統(tǒng)進(jìn)行Rs測定，LI-8100測定時(shí)間設(shè)置為每3 min采集一次數(shù)據(jù)，每小時(shí)對保留和去除枯落層處理交叉各采集10次數(shù)據(jù)，每個(gè)樣點(diǎn)采集3 d數(shù)據(jù)，每1 d從8：00開始采集，至20：00結(jié)束。

土壤溫濕度測定：使用EM50土壤水分?jǐn)?shù)據(jù)采集器+ECH2O土壤濕度傳感器(EM50, USA)進(jìn)行土壤溫濕度數(shù)據(jù)采集。ECH2O土壤濕度傳感器設(shè)置：在針葉和闊葉林樣點(diǎn)保留和去除枯落層處理上各使用一個(gè)ECH2O土壤濕度傳感器進(jìn)行土壤溫濕度的采集。

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

采用Excel 2007進(jìn)行數(shù)據(jù)處理，不同處理間土壤呼吸、土壤溫度和土壤體積含水量之間的差異分析采用單因素方差分析(One-way ANOVA)法進(jìn)行。作圖及相關(guān)性分析采用Origin 9.0(Originlab公司，美國)進(jìn)行。

2 結(jié)果與分析

2.1 不同林型土壤基本性質(zhì)

由表1可知，從土壤層次來看，不同林型SOC、易氧化有機(jī)碳(Readily oxidizable organic carbon, ROC)及全氮(Total nitrogen, TN)含量均表現(xiàn)為：0～5 cm>5～10 cm，土壤pH值表現(xiàn)為：0～5 cm<5～10 cm。0～5 cm層次針葉林SOC、ROC即TN含量分別較5～10 cm層次高出了34.54 %、30.90 %和24.67 %，5～10 cm層次土壤pH值較0～5 cm層次高出了9.54 %。

表1 不同林型表層土壤基本性質(zhì)

從不同植被類型來看，0～5 cm和5～10 cm層次SOC、ROC及TN含量均表現(xiàn)為：針葉林>闊葉林，土壤pH值表現(xiàn)為：針葉林<闊葉林。將兩層次SOC、ROC及TN含量進(jìn)行平均后比較，針葉林SOC、ROC及TN平均含量分別較闊葉林高出了13.73 %、16.09 %和45.66 %，針葉林土壤pH值較高闊葉林低了6.75 %?？梢姡哚樔~林土壤酸化較闊葉林嚴(yán)重，這主要與2種林型不同樹種的分泌物類型、化合物組成等相關(guān)。

2.2 Rs速率變化特征

由圖1-a可知，以闊葉林為例，同一林型，同一處理，在不同測定日期所測Rs速率存在差異，在9月28日、9月29日和10月7日，闊葉林保留枯落層處理的Rs速率變化范圍分別為：(5.67±0.60)～(6.95±0.53)、(4.72±0.53)～(6.23±0.39)和(3.63±0.40)～(5.81±0.22)μmol·m-2·s-1，其日變異系數(shù)分別為6.39 %、7.25 %和11.91 %；闊葉林去除枯落層處理的Rs速率變化范圍分別為：(1.73±0.10)～(1.98±0.10)、(1.24±0.46)～(1.90±0.12)和(1.40±0.30)～(1.82±0.20)μmol·m-2·s-1，其日變異系數(shù)分別為3.91 %、9.64 %和7.73 %。不同測定日期Rs速率最大值和最小值出現(xiàn)的時(shí)間有所差異。在9月28日、9月29日和10月7日，闊葉林保留枯落層Rs速率最大值分別出現(xiàn)在18：00、14：00和17：00，9月28日最小值分別出現(xiàn)在13：00，9月29日和10月7日最小值均出現(xiàn)在20：00。可見，不同測定日期的天氣情況是決定Rs速率的關(guān)鍵因素。針葉林2種處理的變化特征與闊葉林的變化規(guī)律基本一致。

圖1 Rs速率日變化特征Fig.1 Diurnal variation characteristics of soil respiration rate

2.3 枯落層Rs貢獻(xiàn)率

將本研究中，供試3 d的Rs速率按日平均，可比較2種林型Rs速率差異，同時(shí)可評價(jià)兩種林型地表枯落層對整個(gè)Rs的貢獻(xiàn)，如圖2-a可知，不同測定日期日均Rs速率均表現(xiàn)為：針葉林>闊葉林。9月28日、9月29日和10月7日，針葉林日均Rs速率分別為(6.86±0.52)、(6.40±0.63)和(5.90±0.44)μmol·m-2·s-1，闊葉林分別為(6.38±0.58)、(5.40±0.54)和(4.61±0.39)μmol·m-2·s-1，針葉林Rs速率較闊葉林分別高出了7.05 %、15.68 %和21.83 %。這與羅璐等[5]在神農(nóng)架海拔梯度上對4種典型森林Rs的研究結(jié)果不一致，她的研究表明，闊葉林Rs速率顯著高于針葉林。但Hibbard等[6]對闊葉林林和針葉林林年平均Rs速率研究發(fā)現(xiàn)，闊葉林Rs速率為2.40 μmol·m-2·s-1，針葉林Rs速率為2.42 μmol·m-2·s-1，可見本研究結(jié)果與此一致。存在這種差異的原因可能與研究區(qū)域的土壤質(zhì)量、氣候類型、生物條件及森林的成熟程度等有關(guān)。臧逸飛等[7]研究表明，Rs與SOC、TN呈極顯著正相關(guān)。而本研究中針葉林SOC、ROC及TN含量均高于闊葉林，故土壤呼吸速率亦表現(xiàn)為：針葉林>闊葉林。

圖2 不同林型枯落層對Rs的貢獻(xiàn)Fig.2 Contribution of different types of litter to soil respiration

不同林型Rs速率之間差異明顯[8]，Witkamp等[9]研究表明，Rs速率與枯落物的質(zhì)和量及分解速率之間具有顯著的相關(guān)關(guān)系。本研究結(jié)果表明，在供試3 d的枯落層日均Rs速率均表現(xiàn)為：闊葉林>針葉林，9月28日、29日和10月7日，針葉林枯落層日均Rs速率分別為：(3.48±0.37)、(2.99±0.43)和(2.95±0.36)μmol·m-2·s-1，分別占到總Rs速率的50.59 %、46.76 %和49.93 %；闊葉林枯落層日均Rs速率分別為：(4.51±0.37)、(3.62±0.44)和(2.98±0.28) μmol·m-2·s-1，分別占到總Rs速率的70.53 %、67.00 %和64.25 %。且在9月28日，闊葉林日均Rs速率最高，較針葉林高出了22.78 %，其次為9月29日，闊葉林日均Rs速率較針葉林高出了17.32 %，10月7日闊葉林日均Rs速率最低，較針葉林高出了0.83 %。

綜上可知，供試3 d枯落層日均Rs速率均表現(xiàn)為：闊葉林>針葉林，且對總Rs的貢獻(xiàn)率也表現(xiàn)為：闊葉林>針葉林。這主要可能與研究區(qū)域地上地下的生物量的多少、種類、枯落物的質(zhì)和量及生物代謝過程相關(guān)[10-11]。Landsberg等[12]發(fā)現(xiàn)，闊葉林的枯落物較針葉林的枯落物更容易分解。故本研究闊葉林枯落層Rs貢獻(xiàn)率大于針葉林。

2.4 土壤溫、濕度日變化特征

如圖3-a、b所示，在供試的3 d內(nèi)，8：00-20：00時(shí)間段，闊葉林枯落層和去除枯落層處理的土壤溫度分別分布在10.80～15.30和8.30～15.40 ℃之間。供試3 d內(nèi)，枯落層和去除枯落層處理最高溫度均出現(xiàn)在15：00-16：00之間。由于本研究所測時(shí)間段為8：00-20：00，故1 d內(nèi)溫度變化幅度較小，9月28日、9月29日和10月7日枯落層土壤溫度日變化變異系數(shù)分別為5.70 %、3.46 %和9.28 %，去除枯落層處理土壤溫度日變化變異系數(shù)分別為9.07 %、11.43 %和15.14 %；可見，同1 d內(nèi)去除枯落層處理日變異性明顯大于枯落層，說明去除枯落層后的土壤更容易受到外界環(huán)境條件變化的擾動，故枯落層具有抵抗外界干擾，防止土壤溫度異常變化的作用。土壤升溫階段主要為8：00-14：00，供試3 d內(nèi)，9月28日、29日和10月7日枯落層處理平均升溫速率達(dá)到0.34、0.21和0.39 ℃/h，然后在15：00-20：00期間快速降溫，平均降溫速率分別為0.27、0.17和0.50 ℃/h；去除枯落層處理的升溫和降溫時(shí)間段與枯落層處理一致，即平均升溫速率達(dá)到0.43、0.36和0.47 ℃/h，平均降溫速率分別為0.47、0.58和1.03 ℃/h。

圖3 土壤溫、體積含水量日變化特征(2017年)Fig.3 Characteristics of daily change of soil temperature and volume water content(2017)

如圖3-c、d所示，針葉林土壤溫度日變化特征與闊葉林的變化規(guī)律一致，9月28日、29日和10月7日針葉林枯落層平均升溫速率達(dá)到0.46、0.35和0.59 ℃/h，然后在15：00-20：00期間快速降溫，平均降溫速率分別為0.33、0.29和0.18 ℃/h；去除枯落層處理平均升溫速率達(dá)到0.12、0.11和0.16 ℃/h，平均降溫速率分別為0.03、0.06和0.11 ℃/h。

如圖3-e、f所示，測試3 d內(nèi)，枯落層和去除枯落層土壤體積含水率分別在0.013～0.015和0.13～0.18 m3·m-3之間，去除枯落層處理的土壤體積含水率明顯高于枯落層。土壤體積含水率日變化幅度極小，枯落層日變異系數(shù)分別為4.68 %、2.22 %和1.96 %，去除枯落層處理日變異系數(shù)分別為6.92 %、3.31 %和8.64 %。可見，同1 d內(nèi)去除枯落層日變異性明顯大于枯落層處理。相對于土壤溫度日變異性而言，土壤體積含水率日變異性較小。供試3 d內(nèi)，枯落層和去除枯落層處理每1 d土壤體積含水率在最大和最小值出現(xiàn)時(shí)間上存在差異。

由圖3-g、h可知，針葉林土壤體積含水率日變化特征與闊葉林基本一致。測試3 d內(nèi)，枯落層和去除枯落層處理土壤體積含水率分別在0.049～0.071和0.11～0.12 m3·m-3之間，土壤體積含水率日變化幅度極小，枯落層日變異系數(shù)分別為2.51 %、2.44 %和2.70 %，去除枯落層處理日變異系數(shù)分別為1.08 %、0.75 %和2.19 %。

2.5 Rs對土壤溫、體積含水量的響應(yīng)

土壤溫度在影響土壤酶及土壤微生物活性的同時(shí)，通過影響呼吸底物供應(yīng)對Rs產(chǎn)生影響[13]。宋啟亮等[14]對大興安嶺5種類型低質(zhì)林Rs速率日變化的研究發(fā)現(xiàn)，5種類型低質(zhì)林樣地Rs速率與土壤溫度的關(guān)系都適合指數(shù)模型(R2在0.73～0.82范圍) ，邱睿等[15]研究亦表明，Rs速率與土壤溫度之間呈極顯著指數(shù)相關(guān)關(guān)系。將本研究供試3 d的Rs速率及土壤溫度相關(guān)數(shù)據(jù)進(jìn)行擬合，可知，Rs速率與土壤溫度之間呈指數(shù)相關(guān)關(guān)系，這與前人研究結(jié)果一致。由相關(guān)系數(shù)可知，供試3 d內(nèi)，10月7日Rs速率與溫度之間的擬合程度最高，相關(guān)系數(shù)R值分布在0.442～0.912，其余2 d內(nèi)，除了9月28日闊葉林去除枯落層處理的相關(guān)系數(shù)R值為0.728外，其它的相關(guān)系數(shù)R值均小于0.555。表明Rs速率及溫度受天氣變化影響明顯，不同日期所測Rs速率及溫度差異較大。

土壤濕度是影響Rs速率的另一個(gè)重要因素。曹興等[16]研究表明，Rs速率與土壤濕度為極顯著的正相關(guān)關(guān)系。在本研究中，Rs速率與土壤體積含水量之間的關(guān)系可用指數(shù)函數(shù)進(jìn)行擬合，除10月7日，闊葉林去除枯落層處理Rs速率與體積含水量的相關(guān)性達(dá)顯著水平外(P<0.05，n=13)，其它均未達(dá)顯著水平(P>0.05，n=13)。這主要可能與本研究區(qū)域森林郁閉度較大，土壤體積含水量較高，加之森林蓄水持水能力較強(qiáng)，故土壤體積含水量對Rs的影響并不顯著。

表2 Rs與土壤溫度的關(guān)系

表3 Rs速率與土壤體積含水量的關(guān)系

表4 土壤溫度、體積含水量與土壤呼吸的復(fù)合關(guān)系

2.6 土壤溫度、體積含水量與土壤呼吸的復(fù)合關(guān)系

土壤呼吸速率的變化并非受某一單一因子的影響，而是受多種因素的共同作用，其中，土壤溫度和土壤體積含水量的影響尤為重要[17]，且關(guān)系復(fù)雜，諸多研究者在做回歸分析時(shí)選擇的關(guān)系擬合模型差異較大，筆者采用雙變量回歸模型y=at+bw+k(t為土壤溫度，w為土壤體積含水量)對土壤呼吸速率與土壤溫度和體積含水量的關(guān)系進(jìn)行回歸分析，結(jié)果表明，闊葉林添加枯落層處理，土壤溫度和土壤體積含水量兩者擬合的雙變量模型可共同解釋土壤呼吸速率日變化的19.00 %～58.21 %，闊葉林去除枯落層處理可共同解釋土壤呼吸速率日變化的8.05 %～81.76 %；針葉林添加枯落層處理，可共同解釋土壤呼吸速率日變化的2.87 %～76.89 %，去除枯落層處理可共同解釋土壤呼吸速率日變化的17.14 %～87.33 %。由此可知，不同測定日期，土壤呼吸與土壤溫度和體積含水量之間的擬合程度差異較大，且移除枯落層處理更容易受到土壤溫度和體積含水量變化的影響。同時(shí)，相對于土壤體積含水量的擬合程度而言，土壤呼吸更容易受到土壤溫度的影響。

3 討 論

從不同土壤層次來看，2種林型SOC、ROC及TN含量均表現(xiàn)為：0～5 cm>5～10 cm，土壤pH值表現(xiàn)為：0～5 cm<5～10 cm；相對于整個(gè)土體而言，由于較低的pH可增加土壤養(yǎng)分的有效性[18]，加之0～5 cm深度有是本研究中SOC、ROC及TN含量最高的層次，故0～5 cm土壤深度是C、N養(yǎng)分有效性最高的層次，亦是土壤CO2排放的重要來源。從不同林型來看，0～5和5～10 cm層次SOC、ROC及TN含量均表現(xiàn)為：針葉林>闊葉林，土壤pH值表現(xiàn)為：針葉林<闊葉林。本研究結(jié)果與彭艷等[19]對不同林型土壤養(yǎng)分的研究結(jié)果一致，其原因是針葉林枯落物較闊葉林枯落物不易分解，這種較緩慢的分解更有利于針葉林土壤對C、N等養(yǎng)分的固持[20]。

供試3 d內(nèi)，Rs速率以9月28日最高，10月7日最低，且不同測定日期Rs速率最大值和最小值出現(xiàn)的時(shí)間有所差異。表明土壤呼吸日變化受天氣變化的影響較大。不同測定日期日均總Rs速率均表現(xiàn)為：針葉林>闊葉林，枯落層呼吸速率均表現(xiàn)為：闊葉林>針葉林；就總土壤呼吸而言，在本研究中，針葉林土壤SOC、ROC等呼吸底物的含量明顯大于闊葉林，而土壤呼吸速率與底物之間呈正相關(guān)關(guān)系[21]，因此，針葉林土壤呼吸高于闊葉林。供試3 d內(nèi)，針葉枯落層呼吸對總Rs速率的貢獻(xiàn)率分別為：50.59 %、46.76 %和49.93 %；闊葉林枯落層呼吸對總Rs速率的貢獻(xiàn)率分別為：70.53 %、67.00 %和64.25 %。就枯落層呼吸而言，由于針葉中木質(zhì)素含量高于闊葉林，碳氮比遠(yuǎn)大于闊葉林[22]，因此，針葉較闊葉不易分解，故枯落層呼吸表現(xiàn)為闊葉林大于針葉林，闊葉林枯落層呼吸對總Rs速率的貢獻(xiàn)率亦高于針葉林。

在供試的3 d內(nèi)，8：00-20：00時(shí)間段，闊葉林枯落層和去除枯落層處理的土壤溫度分別分布在10.80～15.30和8.30～15.40 ℃之間。且由枯落層和去除枯落的日變異系數(shù)可知，去除枯落層后的土壤更容易受到外界環(huán)境條件變化的擾動?？萋鋵雍腿コ萋鋵犹幚硗寥荔w積含水率分別在0.013～0.015和0.13～0.18 m3·m-3之間，去除枯落層處理的土壤體積含水率明顯高于枯落層。針葉林土壤溫、濕度日變化特征與闊葉林的變化規(guī)律一致。Rs速率與土壤溫、體積含水量之間呈指數(shù)相關(guān)。Rs速率與土壤體積含水量之間的關(guān)系可用指數(shù)函數(shù)進(jìn)行擬合。由雙變量回歸模型可知，不同測定日期，Rs與土壤溫度和體積含水量之間的擬合程度差異較大，且移除枯落層處理更容易受到土壤溫度和體積含水量變化的影響。同時(shí)，相對于土壤體積含水量的擬合程度而言，Rs更容易受到土壤溫度的影響。

4 結(jié) 論

本研究以藏東南色季拉山針-闊葉天然林為研究對象，針對高寒針-闊葉天然林枯落層呼吸對總Rs的貢獻(xiàn)做了相關(guān)研究。結(jié)果表明，2種林型SOC、ROC及TN含量具有明顯的表聚效應(yīng)(0～5 cm>5～10 cm)，土壤pH值表現(xiàn)為：0～5 cm<5～10 cm，整體而言，針葉林表層SOC、ROC及TN含量均高于闊葉林。不同測定日期總Rs速率之間存在差異，但供試3 d內(nèi)均表現(xiàn)為：針葉林>闊葉林，枯落層呼吸速率均表現(xiàn)為：闊葉林>針葉林，枯落層呼吸對總Rs速率的貢獻(xiàn)率表現(xiàn)為：針葉林<闊葉林；8：00-20：00時(shí)間段，闊葉林枯落層和去除枯落層處理的土壤溫度分別分布在10.80～15.30和8.30～15.40 ℃之間。且由枯落層和去除枯落的日變異系數(shù)可知，去除枯落層后的土壤更容易受到外界環(huán)境條件變化的擾動?？萋鋵雍腿コ萋鋵犹幚硗寥荔w積含水率分別在0.013～0.015和0.13～0.18 m3·m-3之間，去除枯落層處理的土壤體積含水率明顯高于枯落層；通過擬合方程可知，供試3 d的Rs速率與土壤溫、濕度之間均呈指數(shù)相關(guān)關(guān)系，由雙變量回歸模型可知，不同測定日期，Rs與土壤溫度和體積含水量之間的擬合程度差異較大，且移除枯落層處理更容易受到土壤溫度和體積含水量變化的影響。同時(shí)，相對于土壤體積含水量的擬合程度而言，Rs更容易受到土壤溫度的影響。