潘永泰，劉貴華，劉星，于孟飛，馬琳*

（1西藏大學(xué) 青藏高原生態(tài)與環(huán)境研究中心，拉薩 850000；2西藏大學(xué) 理學(xué)院，拉薩 850000；3中國(guó)科學(xué)院 武漢植物園&生植物與流域生態(tài)重點(diǎn)實(shí)驗(yàn)室，武漢 430074；4武漢大學(xué) 生命科學(xué)學(xué)院，武漢 430000；5中南民族大學(xué) 生命科學(xué)學(xué)院，武漢 430074）

凍土中富含有機(jī)碳，其所含碳的總量幾乎是大氣中碳的兩倍［1］.在充當(dāng)溫室氣體“匯”的同時(shí)，一定條件下凍土也發(fā)揮“源”的作用.研究證實(shí)：在一定溫度范圍內(nèi)，土壤呼吸速率隨溫度的升高而上升［2］，表明氣候變暖將會(huì)導(dǎo)致更多的CO2排放到大氣之中，從而使溫室效應(yīng)進(jìn)一步加強(qiáng).因此，近40年來(lái)，氣候變暖對(duì)土壤呼吸速率的影響受到了廣泛的關(guān)注［3］.

青藏高原地區(qū)人為干擾強(qiáng)度較低，其生態(tài)系統(tǒng)各反饋過(guò)程摻雜的人為因素較少，這為環(huán)境與氣候變化間關(guān)系研究提供了絕佳的素材［4］.在過(guò)去的50年里，青藏高原的溫度以每10年0.2℃的速度上升，約為觀測(cè)到的全球變暖速度的兩倍，導(dǎo)致了顯著的永久凍土融化和冰川退縮現(xiàn)象［5］.當(dāng)前，氣候變暖仍是青藏高原地區(qū)面臨的重要生態(tài)問(wèn)題.然而目前為止，關(guān)于青藏高原河流濕地中土壤溫室氣體釋放速率及其溫度敏感性方面的研究仍較匱乏.

土壤呼吸溫度敏感性通常被描述為Q10，表示溫度每升高10℃土壤呼吸速率的變化率，它在一定程度上決定了陸地生態(tài)系統(tǒng)碳儲(chǔ)量對(duì)全球變暖的響應(yīng)［6］.普遍認(rèn)為，全球變暖會(huì)增加凈初級(jí)生產(chǎn)力和土壤呼吸.如果土壤呼吸的溫度敏感性高于凈初級(jí)生產(chǎn)力，陸地生態(tài)系統(tǒng)將成為全球變暖背景下的碳源，大氣CO2濃度和全球氣溫將進(jìn)一步升高.相反，陸地生態(tài)系統(tǒng)成為碳匯，全球升溫趨勢(shì)放緩.目前，凈初級(jí)生產(chǎn)力對(duì)溫度升高的響應(yīng)機(jī)制相關(guān)的研究較多［7］，而Q10高度不確定，在許多方面研究不充分并有爭(zhēng)議.

Q10受氣候、土壤和生物因素等的影響，但各因子的影響程度在不同研究中存在分歧.一方面，Q10被證明受到溫度、降水、土壤生化條件的影響［3］.溫度和水分是影響Q10的重要因素，較多的研究表明：Q10在低溫時(shí)較高，高溫時(shí)較低，但Q10也可隨溫度上升而增大，或隨溫度變化而保持穩(wěn)定［3］；土壤碳礦化量隨著溫度升高增大；Q10隨水分升高而顯著升高.為了準(zhǔn)確模擬全球土壤碳庫(kù)對(duì)全球變暖的響應(yīng)，有必要了解Q10空間驅(qū)動(dòng)的氣候或土壤因素［8］.在區(qū)域范圍內(nèi)以土壤有機(jī)碳為基礎(chǔ)的微生物呼吸速率可以預(yù)測(cè)Q10的變化［9］；另一方面，在生態(tài)系統(tǒng)尺度上，Q10存在較大的空間變異，去除混淆因素后，Q10之間存在相似性［8］.室內(nèi)培養(yǎng)能只改變溫度而保持其他因子一致，雖然室內(nèi)培養(yǎng)中采集土樣時(shí)會(huì)對(duì)土壤有擾動(dòng)，但它可以屏蔽溫度以外其他因子（根系呼吸、pH、土壤水分、孔隙和養(yǎng)分限制等）對(duì)土壤呼吸的影響，所以室內(nèi)培養(yǎng)被認(rèn)為是對(duì)有機(jī)碳分解的溫度敏感性的最小偏差估計(jì)［10］.

本文以青藏高原拉薩河流域河岸帶濕地為研究對(duì)象，在全流域內(nèi)采集了26個(gè)位點(diǎn)的土壤樣本，采用室內(nèi)連續(xù)變溫培養(yǎng)法，通過(guò)在短時(shí)間內(nèi)連續(xù)改變同一土壤樣品的溫度并測(cè)量每個(gè)溫度下的呼吸速率［11］，計(jì)算了土壤樣的CO2通量和Q10值，探究了氣候、土壤生化條件、植被、人口經(jīng)濟(jì)、地理等因素對(duì)Q10的影響，以期揭示青藏高原河岸帶濕地土壤Q10的驅(qū)動(dòng)因子，為青藏高原濕地土壤碳動(dòng)態(tài)預(yù)測(cè)提供科學(xué)數(shù)據(jù).

1 材料與方法

1.1 樣品采集及數(shù)據(jù)收集

拉薩河發(fā)源于念青唐古拉山南麓嘉黎里彭措拉孔馬溝.流經(jīng)那曲、當(dāng)雄、林周、墨竹工卡、達(dá)孜、城關(guān)、堆龍德慶，至曲水縣，是雅魯藏布江中游一條較大的支流，全長(zhǎng)495 km，流域面積31760 km2；最大流量2830 m3?s-1，最小流量20 m3?s-1，年平區(qū)流量287 m3?s-1；海拔高度由源頭5500 m到河口3580 m，是世界上最高的河流之一，其補(bǔ)給類型為融雪和降雨補(bǔ)給.

根據(jù)全國(guó)地理信息資源目錄服務(wù)系統(tǒng)將拉薩流域分成12個(gè)流域（圖1）.流域內(nèi)徑流的形成經(jīng)過(guò)自然降雨到通過(guò)地表徑流和地下徑流匯集到流域出口斷面中，流交匯區(qū)域可以解釋流域中部分碳氮的來(lái)源，樣點(diǎn)均勻分布在流域的各個(gè)河流域交匯處來(lái)代表各個(gè)子流域的土壤特性.根據(jù)交通便利及采集可實(shí)現(xiàn)性原則，自拉薩河源頭到河口及其主要支流共選擇26個(gè)位點(diǎn)（其中R7河流匯入口改道，采集樣品標(biāo)注為無(wú)水R7和改道后河流R7+）.于2021年5月，每個(gè)采樣位點(diǎn)取表層0～10 cm河岸帶土壤約500 g，用2 mm篩子除去植物殘?bào)w和石塊，裝入自封袋中放入到4℃的車載冰箱，用干冰包裹在72 h內(nèi)帶回實(shí)驗(yàn)室后-20℃下保存.

圖1 拉薩河流域12條高海拔河流采樣點(diǎn)位置（紅圈）Fig.1 Location of samplingsites（red circles）in thetwelvehigh-elevation rivers of the Lhasa River Basin

本文所用的數(shù)據(jù)中，氣象數(shù)據(jù)來(lái)自WorldClim；植被數(shù)據(jù)：凈初級(jí)生產(chǎn)量（NPP）、植被歸一化指數(shù)（NDVI）來(lái)自國(guó)家青藏高原科學(xué)數(shù)據(jù)中心，蒸騰量（ET）來(lái)自“一帶一路”及其毗鄰區(qū)域地表蒸散發(fā)數(shù)據(jù)集（2015）全球變化科學(xué)研究數(shù)據(jù)出版系統(tǒng)；人口經(jīng)濟(jì)數(shù)據(jù)：人口（Pop）、夜晚燈光量（NL）和全民生產(chǎn)總值（GDP）來(lái)自國(guó)家青藏高原科學(xué)數(shù)據(jù)中心；地理數(shù)據(jù)：流域（Watershed，WS）劃分、河流網(wǎng)絡(luò)、道路等數(shù)據(jù)來(lái)自全國(guó)地理信息資源目錄服務(wù)系統(tǒng)，海拔（Alt）來(lái)自地理空間數(shù)據(jù)云30 m，土地利用數(shù)據(jù)（LU）來(lái)自資源環(huán)境科學(xué)與數(shù)據(jù)中心（2015）；坡度（Slope）、道 路 緩 沖 區(qū)（Roadbuff）、河 流 緩 沖 區(qū)（Riverbuff）依據(jù)以上數(shù)據(jù)在Arcgis中運(yùn)行獲得；雙評(píng)價(jià)（系資源環(huán)境承載能力和國(guó)土空間開(kāi)發(fā)適宜性評(píng)價(jià)）數(shù)據(jù)、降雨侵蝕歸一化指數(shù)（RF）、水源涵養(yǎng)量（WC）、水土流失（W&SL）均采用《生態(tài)保護(hù)紅線劃定指南（2017）》在Arcgis中通過(guò)柵格計(jì)算獲得.

1.2 室內(nèi)培養(yǎng)和變溫設(shè)定

室內(nèi)采用連續(xù)變溫培養(yǎng)法模擬生長(zhǎng)季節(jié)的晝夜溫度動(dòng)態(tài)，探討土壤呼吸對(duì)于升溫和降溫的不同響應(yīng).參考文獻(xiàn)［12］稍加改進(jìn)：從每個(gè)土樣稱取20 g鮮土混勻，重復(fù)3次，并將其分成3份加水；放置到250 mL棕色瓶中維持60%的含水量，同時(shí)放置3個(gè)空瓶作為空白培養(yǎng)，20℃培養(yǎng)箱中復(fù)蘇7 d.根據(jù)實(shí)際采樣環(huán)境的土壤溫度，設(shè)定溫度從0℃增至20℃，再?gòu)?0℃逐步下降到0℃，總共培養(yǎng)23 d.在0、5、10、15、20℃下分別密封培養(yǎng)48、30、20、10、6 h，再用注射器等量置換方法抽取頂空10 mL氣體，通過(guò)氣相色譜儀（Agilent 7890A）同時(shí)測(cè)定氣體中CO2濃度，后續(xù)降溫培養(yǎng)時(shí)間則與升溫相同.待每個(gè)設(shè)定溫度密封培養(yǎng)結(jié)束后，培養(yǎng)瓶敞口3 h使土壤微生物適應(yīng)新的溫度，將培養(yǎng)瓶蓋上橡膠塞，立即用氣泵通過(guò)進(jìn)-出針頭給培養(yǎng)瓶置換新鮮空氣（約12～15 min），使所有培養(yǎng)瓶?jī)?nèi)初始空氣狀況一致.采用稱重法保持水分恒定，稱重選擇在每次采氣之后立即稱量以維持原取土培養(yǎng)重量.

1.3 室內(nèi)土壤指標(biāo)測(cè)定

用1 mol?L-1HCl除去樣品中的無(wú)機(jī)碳后，使用元素分析儀（Elementar Vario MICRO，Germany）測(cè)定樣品的有機(jī)碳（organic carbon，OC）和總氮（total nitrogen，TN）.MBC和MBN采用氯仿熏蒸和K2SO4提?。?3］，通過(guò)總有機(jī)碳-氮分析儀（Multi N/C 3100-OC/TN Analyzer，Analytik Jena，Germany）上機(jī)測(cè)定.

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

CO2排放或吸收速率的計(jì)算方法如下［12］：

式中：R為CO2排放速率，μg·g-1·h-1；V為培養(yǎng)瓶頂空體積，mL；ΔC=（C培養(yǎng)后-C培養(yǎng)前）/t，μmol·mol-1·h-1；ρ為標(biāo)準(zhǔn)CO2密度，1.977 kg·m-3；T為培養(yǎng)溫度，℃；W為土樣質(zhì)量（干重），g；c為土壤有機(jī)碳含量，μg·g-1；22.4為1 mol CO2氣體在標(biāo)準(zhǔn)狀況下的體積，L·mol-1.

土壤溫度敏感性Q10在較大尺度上呈現(xiàn)相似性［14］，T代表培養(yǎng)溫度，通過(guò)公式（2）、（3）推導(dǎo)，以一次指數(shù)模型得到的Q10在不同溫度條件下呈現(xiàn)固定值.

為了探究Q10隨溫度的變化的關(guān)系，如公式（4）在10℃（如20～10℃）溫度區(qū)間計(jì)算Q10值.

式中T1和T2分別為低溫和高溫下的CO2排放速率［12］.

1.5 統(tǒng)計(jì)方法

統(tǒng)計(jì)分析利用SPSS13.0來(lái)進(jìn)行，多重比較均采用Duncan′s新復(fù)極差法測(cè)驗(yàn)，顯著水平為P＜0.05.用克里金（Kriging）插值法對(duì)測(cè)量生化指標(biāo)進(jìn)行空間運(yùn)算.王勁峰等提出的地理檢測(cè)器是一種空間方差分析方法［15］，最初設(shè)計(jì)用于測(cè)試空間分層異質(zhì)性，現(xiàn)在廣泛用于探索導(dǎo)致異質(zhì)性的驅(qū)動(dòng)因素的解釋.與傳統(tǒng)的線性模型不同，地理探測(cè)器不受線性假設(shè)和共線多變量的限制.本文使用地理檢測(cè)器中的q-統(tǒng)計(jì)量來(lái)表示每個(gè)因素（X）對(duì)土壤溫度敏感性Q10在降溫階段（20～10℃，Y）的影響.q統(tǒng)計(jì)量表示因子X(jué)對(duì)結(jié)果Y的解釋力，可以理解為X解釋Y的100q%.q統(tǒng)計(jì)量的數(shù)學(xué)公式［15］如下：

式中：h=1，…，L為變量Y或因子X(jué)的分層（Strata），即分類或分區(qū)；Nh和N分別為層h和全區(qū)的單元數(shù)和分別是層h和全區(qū)的Y值的方差；SSW和SST分別為層內(nèi)方差之和（Sum of Squares Within）和全區(qū)總方差（Sum of Squares Total）；q的值域?yàn)椋?，1］，值越大說(shuō)明Y的空間分異性越明顯，如果分層由自變量X生成，則q值越大表示自變量X對(duì)屬性Y的解釋力越強(qiáng)，反之則越弱.

2 結(jié)果

2.1 土壤理化性質(zhì)

拉薩河岸帶土壤中TN和OC的分布具有相似性，在中游區(qū)域，呈現(xiàn)西邊流域高東邊流域低的特點(diǎn)，匯入雅魯藏布江處較低，流域發(fā)源地較高［圖2（a）、（b）］，TN含量為0.16～2.07 mg·g-1，OC含量為1.40～25.26 mg·g-1.西南部流域土壤C/N比低于東北部流域土壤，而北部流域土壤C/N比則高于南部流域土壤［圖2（c）］.土壤微生物量氮（MBN）含量為1.7～20.92μg·g-1，南部流域土壤中MBN高于北部流域土壤［圖2（d）］；土壤微生物量碳（MBC）含量為68.42～154.26μg·g-1，呈現(xiàn)西北部和中部流域MBC高而南部流域低的特點(diǎn)［圖2（e）］.MBC/MBN自北向南逐漸減小［圖2（f）］，MBC/OC則是自西向東逐漸減?。蹐D2（g）］，MBN/TN最高出現(xiàn)在中部流域［圖2（h）］.

圖2 土壤理化特性流域分布Fig.2 Watershed distribution of soil physicochemical properties

2.2 變溫條件下河岸帶土壤CO2排放特征

拉薩河流域CO2排放速率隨著溫度的升高而增加，在20℃達(dá)到峰值，其后逐漸下降（圖3）.下游各位點(diǎn)CO2平均排放速率顯著高于中游和上游各位點(diǎn)的平均排放速率，而中游和上游差異不明顯.

圖3 變溫條件下河岸帶土壤CO2排放特征Fig.3 Characteristicsof CO2 emission fromriparian soil under variable temperature

2.3 CO2排放的溫度敏感性指數(shù)Q10變化規(guī)律

通過(guò)公式（2）計(jì)算Q10，隨著區(qū)間溫度的升高，土壤溫度敏感性逐步增加（圖4）.升溫（U10～20℃）和降溫（D20～10℃）聚為一類，選取降溫（D20～10℃）進(jìn)行后續(xù)分析.在流域內(nèi)匯入口和下游西邊流域（堆龍曲）表現(xiàn)出相似的特征，在下游東邊流域（墨竹工卡流域）和上游流域（麥地藏布流域）的土壤溫度敏感性較低.CO2排放的溫度敏感性指數(shù)Q10（D20～10℃）與OC和TN呈顯著負(fù)相關(guān)，與MBN/OC和MBN/TN呈顯著正相關(guān)（表1），土壤有機(jī)碳（OC）與土壤總氮（TN）呈顯著正相關(guān).

表1 Q10（D20～10℃）與土壤理化性質(zhì)的Pearson相關(guān)系數(shù)Tab.1 Pearson correlation coefficient of Q10（D20-10℃）with physicochemical soil properties

圖4 CO2排放的Q10變化規(guī)律和Q10（D20～10℃）流域分布Fig.4 Variations of Q10 for CO2 emission and watershed distribution of Q10（D20-10℃）

2.4 基于地理探測(cè)器的土壤呼吸溫度敏感性Q10空間分異性分析

2.4.1 各因子對(duì)土壤呼吸溫度敏感性Q10空間分異性的影響及顯著性分析

通過(guò)氣象、植被、地理、雙評(píng)價(jià)、生化、人口經(jīng)濟(jì)六個(gè)方面衡量Q10的影響因素，用因子探測(cè)確定風(fēng)險(xiǎn)因素對(duì)Q10的影響（圖5）.氣象因素中，降雨（80.76%）＞太陽(yáng)輻射（74.68%）＞5月均溫（25.78%）＞5月最高溫（25.59%）＞年均溫（25.58%）＞5月最低溫（25.01%），P＜0.001；降雨和太陽(yáng)輻射對(duì)Q10有強(qiáng)烈影響，且兩者之間有顯著差異；溫度對(duì)于其影響較弱，6個(gè)參數(shù)均與Q10關(guān)系顯著（圖5Ⅰ）.在雙評(píng)價(jià)中，降雨侵蝕歸一化指數(shù)（79.67%）＞水源涵養(yǎng)量（31.97%）＞水土流失敏感性（13.09%）＞地形起伏因子（5.28%），P＜0.01，前兩個(gè)因素主要與降雨有關(guān)，是主要影響因素，且與Q10關(guān)系顯著（圖5Ⅳ）.地理參數(shù)中，流域?qū)10的影響最大（88.95%，P＜0.001），海拔與Q10關(guān)系顯著（19.53%，P＜0.001），其他因素則不顯著（圖5Ⅲ）.生化因素中OC（88.16%）＞MBC/OC（78.72%）＞MBN/TN（77.30%）＞TN（72.37%）＞C/N（35.14%）＞MBN（8.70%）＞MBC（6.80%）＞MBC/MBN（5.89%），均顯著與Q10相關(guān)（P＜0.01，圖5Ⅴ）.流域內(nèi)的人口量對(duì)Q10也有較大影響（74.43%，P＜0.001，圖5Ⅵ，表2）.

表2 土壤呼吸溫度敏感性驅(qū)動(dòng)因子描述的統(tǒng)計(jì)分析Tab.2 Statistical analysisof soil respiratory temperature sensitivity driver description

圖5 氣象、植被、地理、雙評(píng)價(jià)、生化、人口經(jīng)濟(jì)對(duì)拉薩河流域Q10的解釋力Fig.5 Explanatory power of meteorology，vegetation，geography，dual evaluation，biochemistry，and population economy on Q10 in the Lhasa River Basin

2.4.2 顯著因子對(duì)土壤呼吸溫度敏感性空間分異性的交互作用

交互作用表現(xiàn)出兩個(gè)及兩個(gè)以上因素的聯(lián)合影響（表3），首先檢驗(yàn)變量的獨(dú)立性，結(jié)果表明生化因素、人口、太陽(yáng)輻射、降雨、海拔、流域可以用于交互分析.生化因素和流域的交互作用影響最大（0.90～0.96），人口分布與有機(jī)碳交互作用影響較強(qiáng)（0.88～0.94），此外氣象因素中降雨和太陽(yáng)輻射與其他因素的聯(lián)合對(duì)Q10有強(qiáng)烈影響.水土流失與C/N、MBC/MBN、MBC的交互作用對(duì)土壤敏感性Q10的影響較弱.

表3 成對(duì)因素間的交互作用對(duì)Q10的影響Tab.3 Effectsof interactionsbetween pairwisefactorson Q10

3 討論

流域濕地河岸帶土壤在20℃培養(yǎng)下，下游土壤呼吸速率33.6～74.84μg·g-1·h-1，是中游土壤呼吸速率和上游土壤呼吸速率的兩倍（圖3）.溫度對(duì)拉薩河流域土壤CO2排放具有顯著影響，溫度越高土壤CO2排放量越高，在青藏高原高寒草地［9］，也發(fā)現(xiàn)類似的規(guī)律.

目前，對(duì)于Q10的研究仍然存在諸多不確定性.不同生態(tài)類型Q10值不同：Q10在生態(tài)系統(tǒng)碳循環(huán)模型中被假定為一個(gè)常數(shù)，在CLM（Community Land Model）模 型 中 為1.5，在CASA（Ames-Stanford Approach）模型中為2.但Q10不是常數(shù)，在不同生態(tài)系統(tǒng)間差異顯著：在森林生態(tài)系統(tǒng)中，Q10存在較大空間變異（1.09～5.51）［2］；在青藏高原草甸生態(tài)系統(tǒng)中，Q10介于1.4～2.4，與土壤有機(jī)碳和土壤全氮含量無(wú)顯著的相關(guān)關(guān)系，也無(wú)明顯的空間格局［9］.本文以河岸帶濕地為主，溫度和水分是主要影響Q10的因素；不同研究方法的Q10計(jì)算公式不同［14］，文中公式（2）、（3）計(jì)算Q10值在2.59～39.38，總體比森林生態(tài)系統(tǒng)和草甸生態(tài)系統(tǒng)要高.Q10還受到溫度、水分、呼吸底物、土壤微生物、土壤礦物學(xué)性質(zhì)和土壤質(zhì)地的影響，對(duì)于溫度的響應(yīng)還有待進(jìn)一步研究.Q10對(duì)溫度的響應(yīng)不同，對(duì)溫度變化的敏感性以寒帶（2.3～10.3）最高，溫帶（1.1～14.2）次之，熱帶（1.4～4.6）和亞熱帶（1.4～4.6）最低.

文中依據(jù)公式（2）～（4）計(jì)算得出的Q10值隨著培養(yǎng)溫度的升高也逐步升高.首先，Q10受到土壤中底物（有機(jī)碳和總氮）影響.Q10值在青藏高原拉薩河流域凍土中的空間變化格局均與土壤有機(jī)碳和全氮含量的空間分布特征相對(duì)應(yīng)（圖2、4）.隨著培養(yǎng)溫度升高，底物中微生物活性增加，加速底物消耗.土壤呼吸溫度敏感性（Q10）較大的區(qū)域，土壤碳氮含量較低，呈現(xiàn)顯著的負(fù)相關(guān)關(guān)系［16］.在低海拔區(qū)域，環(huán)境基礎(chǔ)溫度較高，隨著培養(yǎng)溫度升高，消耗底物中碳氮，Q10與土壤中碳氮含量表現(xiàn)出正相關(guān)關(guān)系［16］.因?yàn)榈偷行钥赡芤种泼富钚?，從而抑制有機(jī)碳分解的升溫響應(yīng).

文中Q10受土地利用（LU）單一因素的影響較小，這與前人實(shí)地研究Q10值與土壤溫度呈顯著負(fù)相關(guān)的觀點(diǎn)相悖［17］，但按土地利用方式參與降雨侵蝕歸一化指數(shù)（RF）、水源涵養(yǎng)量（WC）、水土流失（W&SL）指標(biāo)的柵格計(jì)算，這些因素對(duì)Q10有顯著驅(qū)動(dòng)影響（圖6），說(shuō)明土地利用變化的是通過(guò)降雨間接影響Q10.降雨和太陽(yáng)輻射對(duì)Q10的影響較大（圖5）.青藏高原拉薩河流域河岸帶濕地，表現(xiàn)出高蒸騰，強(qiáng)太陽(yáng)輻射，低降雨量，使得該區(qū)域土壤中水分成為限制因素.研究表明：土壤中60%～70%的含水量有利于土壤呼吸的作用的進(jìn)行［18］，水分過(guò)低或過(guò)高均會(huì)抑制土壤中CO2的釋放.因此，僅考慮生化因素的影響不能完全解釋土壤溫度敏感性Q10.生化因素和流域的交互作用可解釋90%的Q10，人口分布與有機(jī)碳交互作用可解釋Q10的95%.

在自然條件下，Q10受到多因素交互作用的影響，文中兩因子共同作用時(shí)均會(huì)增加對(duì)研究區(qū)土壤呼吸空間分布的解釋能力，表現(xiàn)為雙因子增強(qiáng)協(xié)同作用.流域和降雨與各因子間的交互作用對(duì)Q10的影響最強(qiáng)，表明自然因素對(duì)Q10的影響起主要作用.研究區(qū)Q10同時(shí)受到氣象、生化、流域、人口等多個(gè)因子相互作用和影響，各要素通過(guò)其自身或與相互耦合脅迫來(lái)影響Q10的空間變化.地理探測(cè)器雖具有獨(dú)特優(yōu)勢(shì)，能反映驅(qū)動(dòng)因素之間交互作用（表2），任何兩種因子對(duì)Q10值的交互作用均大于單一驅(qū)動(dòng)因子的單獨(dú)作用，但它同時(shí)也存在一定缺陷，影響因子的分類或分級(jí)方法定性成分較大，結(jié)果直接影響其運(yùn)算精度［19］，故地理探測(cè)器還不能較好地解釋各影響因子之間的相互作用機(jī)理.

4 結(jié)論

（1）CO2排放速率和Q10空間分布特征：下游CO2平均排放速率顯著高于中游和上游；Q10在流域內(nèi)匯入口和下游堆龍曲表現(xiàn)出與排放速率相似的特征，在下游墨竹工卡流域和上游麥地藏布流域.

（2）Q10驅(qū)動(dòng)因子探測(cè)：降雨、流域、海拔、有機(jī)碳、總氮、人口對(duì)Q10的影響較大.

（3）Q10交互作用探測(cè)：流域與其他顯著因子的交互作用可解釋90%以上的Q10，有機(jī)碳、降雨和太陽(yáng)輻射與其他顯著因子的交互作用可分別解釋80%、83%和79%以上的Q10.