王云英, 裴薇薇, 郭小偉, 杜巖功*

(1.中國科學(xué)院西北高原生物研究所, 青海 西寧 810008； 2.中國科學(xué)院大學(xué), 北京 100049)

濕地生態(tài)系統(tǒng)土壤和植被碳儲(chǔ)量豐富，碳密度高[1]，高寒濕地是青藏高原高寒生態(tài)系統(tǒng)長(zhǎng)期適應(yīng)高寒氣候環(huán)境的主要植被類型[2]。青藏高原濕地約占中國濕地總面積三分之一[3]，對(duì)氣候變化較為敏感，可作為氣候變化的指示器[4]。同時(shí)，氣候變化對(duì)濕地面積和分布也有重要影響[5]。作為青藏高原特殊生態(tài)系統(tǒng)類型，高寒濕地海拔高、溫度低[6]，面積約13.3×106hm2，在全球碳平衡中發(fā)揮著重要作用，是青藏高原重要水源涵養(yǎng)地[7]。

高寒濕地碳循環(huán)和水循環(huán)之間存在密切耦合關(guān)系[8]，可為環(huán)境補(bǔ)充地下水、為野生動(dòng)物提供棲息地等[9]。較高的水分利用效率表明可在有限時(shí)間內(nèi)獲得更多光合產(chǎn)物[10]。2015—2017年三年期間大九湖泥潭濕地生態(tài)系統(tǒng)水分利用效率日均值大于0，為3.26 gC·kg-1H2O[11]。季節(jié)尺度上，金河灣濕地水分利用效率在冬季呈現(xiàn)最高值，而在夏季呈現(xiàn)最低值，全年均值為4.09 mg·g-1[12]。生態(tài)系統(tǒng)碳呼吸是陸地和大氣[13]、草地生態(tài)系統(tǒng)碳交換[14]過程的主要環(huán)節(jié)，其包括地上部分植物呼吸和土壤呼吸兩部分[15]。土壤呼吸動(dòng)態(tài)變化直接影響著全球碳收支平衡[16]。植物呼吸排放的CO2量約占整個(gè)生態(tài)系統(tǒng)的碳交換的30%～70%[17]。

作為一種調(diào)控因子，氣象因子平均氣溫、年降水量對(duì)濕地面積變化、景觀格局改變[9]、青藏高原環(huán)境變化[18]等有著重要反饋?zhàn)饔?。此外，氣象因子光合有效輻射、若爾蓋濕地生態(tài)系統(tǒng)碳吸收主要受光合有效輻射影響[19]。土壤溫度是影響青藏高原高寒濕地[20]、三角洲濱海濕地[21]、黃河三角洲蘆葦濕地碳通量的主要調(diào)控因子[22]。光量子通量密度是影響鄱陽湖南磯濕地凈生態(tài)系統(tǒng)CO2交換量的主要調(diào)控因子[19]。有研究表明，當(dāng)太陽輻射強(qiáng)度較低時(shí),光合有效輻射的增加會(huì)引起CO2凈交換量增加[23]。增溫促進(jìn)了河源濕地生態(tài)系統(tǒng)呼吸[24]，使得小泊湖濕地生態(tài)系統(tǒng)整個(gè)生長(zhǎng)季呼吸速率提升了24.50%[25]。青海湖高寒濕地生態(tài)系統(tǒng)水分利用效率主要受降水量變化調(diào)控，與溫度之間呈現(xiàn)極顯著正相關(guān)關(guān)系[26]。草甸濕地生態(tài)系統(tǒng)水分利用效率受到凈輻射和空氣溫度影響[8]。土壤溫度由于其自身對(duì)環(huán)境變化較為敏感，進(jìn)而作為青海湖濕地生態(tài)系統(tǒng)CO2通量的主要環(huán)境因子[27]。

目前，大部分研究證明濕地生態(tài)系統(tǒng)總體表現(xiàn)為碳匯[28]。碳通量作為高寒濕地生態(tài)系統(tǒng)研究重點(diǎn)之一，對(duì)于碳通量變化及其調(diào)控機(jī)制研究涉及的很多，但現(xiàn)階段對(duì)于高寒帕米爾苔草濕地生態(tài)系統(tǒng)基于通量變化水分利用效率研究涉及的較少。故本研究借助于2004—2009年間青海海北高寒帕米爾苔草濕地生態(tài)系統(tǒng)碳通量數(shù)據(jù)，計(jì)算出生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力和蒸散量的變化，進(jìn)而計(jì)算出水分利用效率值，探討高寒草甸濕地生態(tài)系統(tǒng)碳通量和水分利用效率變化。并耦合氣象因子進(jìn)一步明晰高寒草甸生態(tài)系統(tǒng)碳通量數(shù)據(jù)和水分利用效率的主要?dú)庀笳{(diào)控因子。

1 材料與方法

1.1 研究區(qū)概況

本項(xiàng)研究是在中國科學(xué)院海北高寒草甸生態(tài)系統(tǒng)定位站(海北站，37°36′51.12″N，101°20′13.92″E，海拔3 235 m)附近的典型高寒帕米爾苔草濕地生態(tài)系統(tǒng)上進(jìn)行。海北站地處青藏高原東北隅祁連山北支冷龍嶺東段南麓的大通河河谷。區(qū)域多年平均氣溫為-1.7℃，年降水量較高，約580 mm，降水主要集中于暖季的5—9月，占年降水量的80%。濕地植被以適冷、濕生、多年生地下芽、具有發(fā)達(dá)通氣組織的草本植物莎草科、毛茛科植物為主,濕地中央以帕米爾薹草(Carexpamirensis)為優(yōu)勢(shì)種,邊緣以西藏嵩草(Kobresiatibetica)為優(yōu)勢(shì)種。

1.2 研究方法與數(shù)據(jù)分析的采集

1.2.1生態(tài)系統(tǒng)呼吸速率和生態(tài)系統(tǒng)CO2碳凈交換量(Net ecosystem exchange,NEE)數(shù)據(jù)和生態(tài)系統(tǒng)呼吸速率(Ecosystem respiration,Re) 站點(diǎn)通量數(shù)據(jù)來自國家生態(tài)系統(tǒng)觀測(cè)研究網(wǎng)絡(luò)科技資源服務(wù)系統(tǒng)(http://rs.cern.ac.cn/)，數(shù)據(jù)時(shí)間尺度為2004—2009年。利用通量數(shù)據(jù)Re、NEE和潛熱通量(Latent heat flux,LE)計(jì)算生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力(Gross primary productivity，GPP)和蒸散量[29]。

蒸散量ET計(jì)算公式如下：

ET=LE/λ
λ=(2500-2.4Ta)×103

其中λ為蒸發(fā)潛熱,單位是J·kg-1,式中Ta是空氣溫度,單位是℃。LE是潛熱通量，單位是W·m-2。

總初級(jí)生產(chǎn)力GPP的計(jì)算公式如下：

GPP=Re-NEE

式中，GPP為植被總初級(jí)生產(chǎn)力，單位為gC·m-2；Re為生態(tài)系統(tǒng)呼吸量，單位為gC·m-2；NEE為生態(tài)系統(tǒng)CO2凈交換量，單位為gC·m-2。

水分利用效率的計(jì)算公式如下：

式中,WUE為水分利用效率,單位為gC·m-2·mm-1；GPP和ET分別為生態(tài)系統(tǒng)總初級(jí)生產(chǎn)力和實(shí)際蒸散量,單位分別為gC·m-2和mm。

1.2.2氣象數(shù)據(jù) 本研究中氣象數(shù)據(jù)均來自于青海海北高寒草地生態(tài)系統(tǒng)國家野外科學(xué)觀測(cè)研究站氣象觀測(cè)站。

1.3 數(shù)據(jù)處理

數(shù)據(jù)整理采用Excel，均值、標(biāo)準(zhǔn)誤、顯著性檢驗(yàn)、結(jié)構(gòu)方程模型的計(jì)算等采用R語言進(jìn)行計(jì)算處理，進(jìn)一步結(jié)合氣象因子結(jié)合結(jié)構(gòu)方程模型探討生態(tài)系統(tǒng)CO2碳凈交換量、生態(tài)系統(tǒng)呼吸速率及水分利用效率變化與氣象因子的關(guān)系。最初通過選取間接影響達(dá)到顯著性水平的氣象因子，最后參與結(jié)構(gòu)方程模型中的為氣象因素：風(fēng)速、相對(duì)濕度、空氣溫度、降水量和凈輻射。比較擬合系數(shù)(CFI)，Tucker-Lewis指數(shù)(TLI)和近似誤差均方根(RMSEA)用來模型評(píng)估。最大似然法，多被結(jié)構(gòu)方程采用的卡方估計(jì)來評(píng)估模型卡方值(χ2)和路徑系數(shù)。方程的擬合是在R中使用lavaan包實(shí)現(xiàn)的。運(yùn)用Origin 2019進(jìn)行圖件的繪制。

2 結(jié)果與分析

2.1 濕地生態(tài)系統(tǒng)碳水通量和水分利用效率年際變化特征

年際尺度上，三者均呈極顯著上升趨勢(shì)(P<0.01，圖1)。NEE年內(nèi)變化呈“V”型波動(dòng)式變化趨勢(shì)，且下降幅度遠(yuǎn)大于上升幅度生態(tài)系統(tǒng)呼吸速率和水分利用效率變化呈倒“V”型變化趨勢(shì)。6年來，NEE,Re和WUE均值分別為0.33 gC·m-2·d-1，1.68 gC·m-2·d-1，0.52 gC·m-2·m-1。

圖1 CO2通量、生態(tài)系統(tǒng)呼吸速率和水分利用效率年際變化Fig.1 Interannual variation of ecosystem net exchange of CO2,ecosystem respiration rate and water use efficiency

2.2 濕地生態(tài)系統(tǒng)碳水通量和水分利用效率季度變化特征

季節(jié)尺度上，連續(xù)6年間非生長(zhǎng)季生態(tài)系統(tǒng)碳排放總量大于生長(zhǎng)季碳吸收總量(圖2a)。生態(tài)系統(tǒng)呼吸速率和水分利用效率均呈現(xiàn)出生長(zhǎng)季值大于非生長(zhǎng)季值(圖2b；圖2c)。整體而言，生長(zhǎng)季NEE,Re和WUE與非生長(zhǎng)季變化呈現(xiàn)出較大差異(圖2d)。

圖2 生態(tài)系統(tǒng)CO2碳凈交換量、生態(tài)系統(tǒng)呼吸速率和水分利用效率季節(jié)變化Fig.2 Seasonal changes of ecosystem net exchange of CO2,Ecosystem respiration and water use efficiency

2.3 高寒濕地生態(tài)系統(tǒng)水分利用效率對(duì)氣象因子的響應(yīng)

結(jié)構(gòu)方程模型結(jié)果表明，在整個(gè)研究期間，降水量、風(fēng)速、凈輻射、空氣溫度和相對(duì)濕度是影響高寒濕地NEE的主要影響因子，且均達(dá)到極顯著性水平(P<0.01)，相對(duì)濕度對(duì)生態(tài)系統(tǒng)NEE的顯著性影響最大，直接影響和間接影響效應(yīng)系數(shù)值分別為-0.21和-0.09，凈輻射影響次之，直接和間接影響的效應(yīng)值分別為-0.16和0.19；相對(duì)濕度和空氣溫度是影響高寒濕地生態(tài)系統(tǒng)呼吸速率的極顯著正效應(yīng)影響因子(P<0.01)。直接和間接效應(yīng)值分別為0.20,0.72和0.22,0.04；降水量和空氣溫度是顯著影響高寒濕地生態(tài)系統(tǒng)水分利用效率變化的主要影響因子(P<0.05)。兩者直接和間接效應(yīng)值分別為0.05,0.08和0.02,0.01。

圖3 環(huán)境因子對(duì)生態(tài)系統(tǒng)CO2凈交換量、生態(tài)系統(tǒng)呼吸速率和水分利用效率影響Fig.3 Structural equation diagram of the influence of environmental factors on ecosystem net exchange of CO2,ecosystem respiration rate and water use efficiency注：實(shí)線表示正相關(guān)，虛線表示負(fù)相關(guān)；“***”和“*”表示在P<0.01和P<0.05達(dá)到顯著性水平，P表示降水量，Ws表示風(fēng)速，Rn表示凈輻射，Ta表示空氣溫度，RH表示相對(duì)濕度Note:The solid line represents the positive correlation,the dashed line represents the negative correlation;“***” and “**” indicate significant levels at P<0.01 and P<0.05,P indicates precipitation,Ws indicates wind speed,and Rn indicates net radiation,Ta means air temperature,RH means relative humidity

3 討論

本研究中，2004—2009連續(xù)6年間生態(tài)系統(tǒng)呼吸速率、生態(tài)系統(tǒng)CO2凈交換量和水分利用效率均呈極顯著上升趨勢(shì)。季節(jié)尺度上，6年間非生長(zhǎng)季生態(tài)系統(tǒng)碳排放總量大于生長(zhǎng)季碳吸收總量。這與對(duì)寒區(qū)濱河濕地對(duì)碳通量的季節(jié)變化研究結(jié)果一致[30]。究其原因可能是因?yàn)樵谥参锷L(zhǎng)季晝夜溫差大，有助于植物光合和光合產(chǎn)物的積累[31]。對(duì)于生態(tài)系統(tǒng)呼吸速率和水分利用效率均呈現(xiàn)出生長(zhǎng)季值大于非生長(zhǎng)季值?？赡苁且?yàn)樵谥脖簧L(zhǎng)季，溫濕度條件適宜時(shí)，其光合作用較強(qiáng)，植被生長(zhǎng)良好，故植物生態(tài)系統(tǒng)呼吸速率和水分利用效率在植被生長(zhǎng)季較高[32]。

植被自身和微生物的活動(dòng)也會(huì)影響生態(tài)系統(tǒng)CO2凈交換量[16]。本研究中，凈輻射和相對(duì)濕度、風(fēng)速對(duì)生態(tài)系統(tǒng)CO2凈交換量存在極顯著負(fù)效應(yīng)影響。凈輻射通過影響植物光合作用來影響高寒濕地夏季生態(tài)系統(tǒng)CO2凈交換量[16]。青海隆寶灘濕地生態(tài)系統(tǒng)CO2凈交換量與土壤濕度之間存在顯著負(fù)相關(guān)關(guān)系，濕地水位降低有利于氧氣的進(jìn)入，加速了土壤有機(jī)物質(zhì)分解，進(jìn)而導(dǎo)致CO2排放增加[4]。

生態(tài)系統(tǒng)呼吸速率是草地生態(tài)系統(tǒng)進(jìn)行碳交換的主要環(huán)節(jié)[11]。本研究中，生態(tài)系統(tǒng)呼吸速率主要受到相對(duì)濕度和氣溫的影響，且均呈正相關(guān)影響。氣溫的變化會(huì)引起植物物候和植被生產(chǎn)力的變化，進(jìn)而伴隨著生態(tài)系統(tǒng)呼吸速率的變化[33]。神農(nóng)架大九湖泥潭濕地土壤溫度通過影響植物和微生物呼吸酶活性、對(duì)氣體和底物擴(kuò)散速率等影響生長(zhǎng)季夜間生態(tài)呼吸速率[34]。

降水量和溫度是調(diào)控水分利用效率變化的主要因子，且與兩者之間均呈正相關(guān)。大九湖泥潭濕地生態(tài)系統(tǒng)夏季水分利用效率變化與土壤溫度之間呈現(xiàn)正相關(guān)關(guān)系，其余季度呈負(fù)相關(guān)關(guān)系，究其原因是因?yàn)橥寥罍囟戎苯佑绊懙乇硭终舭l(fā)和通過影響酶活性和影響氣體的擴(kuò)散速度來影響水分利用效率變化[15]。降水量的變化對(duì)高寒濕地生態(tài)系統(tǒng)水分利用效率的影響存在閾值效應(yīng)，在未達(dá)到閾值時(shí)，兩者之間呈正相關(guān)關(guān)系，超過閾值后，水分利用效率變化基本不變，這是因?yàn)樯鷳B(tài)系統(tǒng)生產(chǎn)力不僅受到降水量的控制[35]，還受到溫度的控制，進(jìn)而影響水分利用效率值[26]。

4 結(jié)論

年際尺度上，生態(tài)系統(tǒng)CO2凈交換量、生態(tài)系統(tǒng)呼吸速率和水分利用效率均呈極顯著上升趨勢(shì)；季節(jié)尺度上，連續(xù)6年間，非生長(zhǎng)季生態(tài)系統(tǒng)碳排放總量大于生長(zhǎng)季碳吸收總量。針對(duì)于生態(tài)系統(tǒng)呼吸速率和水分利用效率均呈現(xiàn)出生長(zhǎng)季值大于非生長(zhǎng)季值。相對(duì)濕度和凈輻射、相對(duì)濕度和空氣溫度、降水量和空氣溫度分別是影響青海海北高寒濕地生態(tài)系統(tǒng)CO2通量、生態(tài)系統(tǒng)呼吸速率和水分利用效率的主要?dú)庀笠蜃印?/p>