王松年，王云琦,2＊，王 凱，馮印成，王杰帥

（1.北京林業(yè)大學(xué)水土保持學(xué)院，重慶三峽庫(kù)區(qū)森林生態(tài)系統(tǒng)教育部野外科學(xué)觀測(cè)研究站，北京 100083；2.北京林業(yè)大學(xué)水土保持學(xué)院，重慶縉云山三峽庫(kù)區(qū)森林生態(tài)系統(tǒng)國(guó)家定位觀測(cè)研究站，北京 100083；3.蘇州高新區(qū)（虎丘區(qū)）林業(yè)站，江蘇 蘇州 215000）

自工業(yè)革命以來(lái)，大氣中CO2等溫室氣體濃度的不斷增加引起氣候變化和全球增溫，從而使人們開(kāi)始關(guān)注陸地生態(tài)系統(tǒng)的碳循環(huán)過(guò)程[1]。森林生態(tài)系統(tǒng)是陸地生態(tài)系統(tǒng)的主體，其占比陸地生態(tài)系統(tǒng)的31%[2-3]，其作為陸地生態(tài)系統(tǒng)的重要碳匯[4]，對(duì)實(shí)現(xiàn)碳中和起到極為關(guān)鍵的作用。因此，對(duì)于明晰森林生態(tài)系統(tǒng)的碳交換量的變化顯得尤為重要。

大氣湍流是生態(tài)系統(tǒng)和大氣進(jìn)行氣體和能量交換的主要媒介，渦協(xié)方差相關(guān)法的誕生就是為了記錄這一過(guò)程的物質(zhì)和能量的量[5]，是目前測(cè)定大氣與陸地生態(tài)系統(tǒng)之間氣體和能量交換的最直接的方法[6-7]。在實(shí)際應(yīng)用過(guò)程中，渦度相關(guān)技術(shù)的監(jiān)測(cè)受到多種環(huán)境因素的影響，如生態(tài)系統(tǒng)類(lèi)型、觀測(cè)高度、下墊面地形、大氣因素等都會(huì)影響其監(jiān)測(cè)結(jié)果的準(zhǔn)確性[5,8]。需要通過(guò)對(duì)數(shù)據(jù)處理后的能量閉合和源區(qū)分析對(duì)其的適用性和準(zhǔn)確性進(jìn)行驗(yàn)證。許多學(xué)者都利用渦度相關(guān)技術(shù)，對(duì)復(fù)雜的森林生態(tài)系統(tǒng)與大氣之間的CO2交換進(jìn)行過(guò)探究[9-10]，同時(shí)針對(duì)不同類(lèi)型的生態(tài)系統(tǒng)和站點(diǎn)進(jìn)行比較分析。目前大部分的研究表明，森林生態(tài)系統(tǒng)為碳匯，只有一小部分的研究得出的結(jié)論為碳源且包括幼齡林在內(nèi)[11-12]。盡管森林生態(tài)系統(tǒng)的碳匯作用得到了極大多數(shù)學(xué)者的廣泛認(rèn)可，但由于環(huán)境因子的影響，使得全球森林生態(tài)系統(tǒng)與大氣間的CO2交換量在不斷變化，所以，對(duì)于森林生態(tài)系統(tǒng)與大氣之間的碳交換監(jiān)測(cè)依舊重要[13]。

重慶縉云山地處三峽庫(kù)區(qū)的尾端，是重要的生態(tài)屏障，其中針闊葉混交林是該地區(qū)的典型林分之一，也是我國(guó)西南山地森林生態(tài)系統(tǒng)的典型代表。王倩、王杰帥等都曾利用二次坐標(biāo)旋轉(zhuǎn)對(duì)其碳通量進(jìn)行分析[10,14]。本研究以2019.11—2020.10縉云山針闊葉混交林通量監(jiān)測(cè)數(shù)據(jù)為基礎(chǔ)，坐標(biāo)旋轉(zhuǎn)處理上采用平面擬合法，通過(guò)湍流分析，能量閉合和源區(qū)分析，來(lái)對(duì)渦度相關(guān)對(duì)其的適用性進(jìn)行探究，同時(shí)對(duì)該地區(qū)碳通量的時(shí)間變化進(jìn)行分析。以期為預(yù)測(cè)未來(lái)碳通量變化下的森林動(dòng)態(tài)和經(jīng)營(yíng)提供科學(xué)依據(jù)。

1 研究區(qū)概況

研究區(qū)所在地為重慶市北碚區(qū)縉云山國(guó)家級(jí)自然保護(hù)區(qū)國(guó)家定位觀測(cè)研究站內(nèi)（106°17'～106°24' E、29°41'～29°52' N），該保護(hù)區(qū)的整體面積約7 600 hm2，海拔200～952 m，地處三峽庫(kù)區(qū)的尾端，是三峽庫(kù)區(qū)的重要生態(tài)屏障?？N云山氣候?yàn)榈湫偷膩啛釒Ъ撅L(fēng)性濕潤(rùn)氣候，降雨充沛，年平均降水1 611.8 mm，主要集中在生長(zhǎng)季（4—10月）；年平均蒸發(fā)量 777.1 mm；該地區(qū)受局部小氣候影響潮濕多霧，平均全年有霧日達(dá)89.8 d。研究區(qū)為針闊葉混交林森林生態(tài)系統(tǒng)，主要的優(yōu)勢(shì)樹(shù)種為四川山礬（Symplocos setchuensis Brand）、廣東山胡椒（Lindera kwangtungensis Allen）、楊桐（Adinandra bockiana Pritzel ex Diels）、馬尾松（Pinus massoniana Lamb）和杉木（Cunninghamia lanceolata Hook）。

2 材料與方法

2.1 通量塔的儀器配置

重慶縉云山通量塔位于觀測(cè)站的針闊葉混交林內(nèi)，塔高35 m，配備的觀測(cè)系統(tǒng)包括兩個(gè)部分，即開(kāi)路式渦度相關(guān)監(jiān)測(cè)系統(tǒng)和微氣象觀測(cè)系統(tǒng)，具體的儀器布設(shè)詳見(jiàn)表1。

表1 通量塔儀器布設(shè)Table 1 Instrument layout of flux tower

渦度相關(guān)法可以直接測(cè)定植物群落與大氣之間的氣體交換通量，其原理是通過(guò)測(cè)定大氣中湍流運(yùn)動(dòng)產(chǎn)生的風(fēng)速脈動(dòng)和物理量脈動(dòng)的協(xié)方差來(lái)測(cè)定物質(zhì)或能量的通量。

大氣與植物群落之間的CO2湍流通量的公式為：

Fc（mg·m?2·s?1）為 CO2的湍流通量，ρc′為CO2密度脈動(dòng)，ω′為垂直風(fēng)速的脈動(dòng)量。

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

原始數(shù)據(jù)處理主要在EddyPro-7.0.6中進(jìn)行，具體方法如下：

（1）采用Vicker提出的處理辦法對(duì)原始數(shù)據(jù)進(jìn)行質(zhì)量檢測(cè)，并將野點(diǎn)及異常值進(jìn)行剔除。

（2） 采 用Mauder和Foken[15]提 出的0～2級(jí)分級(jí)法作為本研究數(shù)據(jù)的質(zhì)量分級(jí)標(biāo)準(zhǔn)。

（3）坐標(biāo)校正，考慮到研究區(qū)位于山區(qū)，下墊面不平穩(wěn)，且生態(tài)系統(tǒng)結(jié)構(gòu)較為復(fù)雜，所以采用平面擬合法[16]對(duì)其進(jìn)行校正。

（4）通量校正，主要包括：超聲虛溫修正、WPL校正、頻率響應(yīng)校正、夜間通量校正。

（5）數(shù)據(jù)插補(bǔ)，采用由Max Planck Institute for Biogeochemistry開(kāi)發(fā)的REddyProcWeb在線數(shù)據(jù)處理工具對(duì)不連續(xù)的通量數(shù)據(jù)進(jìn)行插補(bǔ)。

后期的數(shù)據(jù)處理和制圖分別在Excel、MATLAB、Tovi和Origin中進(jìn)行。

2.3 能量閉合處理

根據(jù)能量守恒定律，輸入生態(tài)系統(tǒng)的能量將會(huì)以某種形式輸出，因此能量的閉合與否是評(píng)判通量數(shù)據(jù)是否可用的衡量標(biāo)準(zhǔn)之一。森林生態(tài)系統(tǒng)中的能量閉合（EBR）計(jì)算公式為：

Rn為太陽(yáng)凈輻射，G為土壤熱通量，H為顯熱通量，LE為潛熱通量，F(xiàn)s為冠層儲(chǔ)熱，Q為其他能量輸出項(xiàng)。

由于冠層儲(chǔ)熱和其他能量輸出項(xiàng)所占比例很少，且部分?jǐn)?shù)據(jù)在測(cè)量中存在較大難度，所以大量的學(xué)者在研究能量閉合時(shí)都將其忽略不計(jì)，本研究中只利用（Rn-G）和（H+LE）來(lái)進(jìn)行能量閉合的計(jì)算。

3 結(jié)果與分析

3.1 湍流數(shù)據(jù)質(zhì)量評(píng)價(jià)

質(zhì)量評(píng)價(jià)是0和1的為高質(zhì)量通量數(shù)據(jù)，圖中顯示為藍(lán)色；質(zhì)量評(píng)價(jià)是2的數(shù)據(jù)為中低質(zhì)量數(shù)據(jù)，圖中顯示為黃色（圖1）。圖中2019.12.06—2019.12.19數(shù)據(jù)缺失是由于通量塔斷電導(dǎo)致。質(zhì)量評(píng)價(jià)為0和1的高質(zhì)量數(shù)據(jù)占比為77%，質(zhì)量評(píng)價(jià)為2的中低質(zhì)量數(shù)據(jù)占比16%，缺失數(shù)據(jù)為7%。結(jié)果表明，77%的數(shù)據(jù)趨于理想的穩(wěn)態(tài)條件，并且很好的符合湍流觀測(cè)的發(fā)展特點(diǎn)。

圖1 顯熱通量質(zhì)量等級(jí)分級(jí)Fig.1 Quality classification of sensible heat flux

由圖2可見(jiàn)，質(zhì)量評(píng)價(jià)為0和1的高質(zhì)量數(shù)據(jù)占比為61%，質(zhì)量評(píng)價(jià)為2的中低質(zhì)量數(shù)據(jù)占比20%，缺失數(shù)據(jù)為19%。結(jié)果表明，61%的數(shù)據(jù)趨于理想的穩(wěn)態(tài)條件。并且很好的符合湍流觀測(cè)的發(fā)展特點(diǎn)。

圖2 潛熱通量質(zhì)量等級(jí)分級(jí)Fig.2 Quality classification of latent heat flux

碳通量的數(shù)據(jù)質(zhì)量容易受到降水的影響，樣地所在地全年多雨（4—10月為雨季），且非雨季時(shí)，由于受山地小氣候影響，該山地長(zhǎng)時(shí)期水霧籠罩。從圖3可以看出，CO2湍流的高質(zhì)量數(shù)據(jù)占比為62%，中低質(zhì)量數(shù)據(jù)為23%，缺失數(shù)據(jù)為15%。符合湍流觀測(cè)的發(fā)展特點(diǎn)。

圖3 碳通量質(zhì)量等級(jí)分級(jí)Fig.3 Quality classification of carbon flux

3.2 能量閉合與通量足跡分析

3.2.1 能量閉合分析 分別對(duì)25 m和35 m處的通量觀測(cè)系統(tǒng)所觀測(cè)到的生長(zhǎng)季通量數(shù)據(jù)進(jìn)行了能量閉合分析（圖4）。該觀測(cè)系統(tǒng)下層能量閉合率為0.73，上層能量閉合率為0.82，均在國(guó)際認(rèn)可的能量閉合范圍之內(nèi)（0.7～0.9），表明能量閉合情況良好。

圖4 能量閉合分析Fig.4 Diagram of energy closure analysis

3.2.2 通量足跡分析 本研究中的通量足跡如圖5所示。左圖為全年白天通量足跡模擬，右圖為全年夜間通量足跡模擬。有色區(qū)域即為通量貢獻(xiàn)區(qū)，顏色越亮代表通量觀測(cè)的數(shù)據(jù)來(lái)源越密集，閉合圈線代表不同累積通量貢獻(xiàn)率的范圍。因?yàn)橥克^測(cè)的通量數(shù)據(jù)一般來(lái)源于上風(fēng)向，因此貢獻(xiàn)區(qū)形狀為橢圓，長(zhǎng)軸在東北和西南方向，同時(shí)可以看出高亮區(qū)域大部分集中在東北側(cè)，即全年大部分月份主要通量觀測(cè)區(qū)域來(lái)源于東北方。全年風(fēng)速風(fēng)向變化如風(fēng)玫瑰圖（圖6）所示，主風(fēng)向來(lái)自東北方向，這與通量的高貢獻(xiàn)區(qū)相吻合，說(shuō)明通量數(shù)據(jù)來(lái)源可信且質(zhì)量良好。

圖5 全年通量足跡（左為白天，右為夜間）Fig.5 Annual flux footprint （day on the left and night on the right）

圖6 風(fēng)速玫瑰圖Fig.6 Wind speed rose

通量的峰值年貢獻(xiàn)區(qū)的年平均距離為82.25 m，10%、30%、50%、70%和90%累計(jì)通量貢獻(xiàn)區(qū)年平均距離分別為28.23、70.31、107.22、150.38 m和225.29 m（表2）。由表2可以看出9月份為貢獻(xiàn)區(qū)距離最大月，峰值距離為91.12 m，90%累計(jì)貢獻(xiàn)距離為249.59 m，最小月份為4月，峰值距離為83.43 m，90%累計(jì)貢獻(xiàn)距離為228.55 m。貢獻(xiàn)區(qū)距離在正常范圍內(nèi)，且不同月份間相差不大。

表2 通量貢獻(xiàn)區(qū)距離Table 2 Distance of flux contribution area

3.3 碳通量的變化分析

3.3.1 各月份碳通量的日變化特征 利用MATLAB對(duì)全年半小時(shí)通量數(shù)據(jù)進(jìn)行平均，得到各月份同一時(shí)刻下半小時(shí)的平均值，得到全年各月同一時(shí)刻的Fc的日變化，由圖7可以看出，各月份的Fc日變化大體相同，均成U字形，F(xiàn)c為負(fù)值時(shí)，說(shuō)明生態(tài)系統(tǒng)在吸收CO2，當(dāng)Fc為正值時(shí)，說(shuō)明生態(tài)系統(tǒng)在釋放CO2，也就是說(shuō)在全年各個(gè)月份當(dāng)中，縉云山針闊葉混交林生態(tài)系統(tǒng)在白天都為碳匯，在夜間則為碳源。主要原因是受到輻射的影響，從而使生態(tài)系統(tǒng)在光合與呼吸作用之間相互切換。其中碳匯的波動(dòng)范圍為：?0.97～?0.01 mg·m?2·s?1，最小值出現(xiàn)在8月的12:00，最大值出現(xiàn)在5月的19:00。碳源多集中在夜間，其波動(dòng)區(qū)間是：0.04～0.32 mg·m?2·s?1，最大值出現(xiàn)在 12月的0:00點(diǎn)，最小值出現(xiàn)在1月份的6:30。

圖7 各月份碳通量日變化Fig.7 Daily variation of carbon flux in different months

研究期間各月份碳匯峰值分布范圍為11:30～15:00，6月最早達(dá)到碳匯峰值，1月最晚達(dá)到峰值；碳源峰值主要分布范圍為 20:00—次日 3:00。除了11月，其他各月均出現(xiàn)兩個(gè)碳匯峰值，即在中午段Fc出現(xiàn)增大的現(xiàn)象。碳匯能力最大的是8月份，日均 Fc 為?0.24 mg·m?2·s?1，最小的是12月份，F(xiàn)C 為 0.03 mg·m?2·s?1。每個(gè)月份碳匯碳源的轉(zhuǎn)化時(shí)刻也不盡相同，最早由碳源轉(zhuǎn)為碳匯的是7月和8月均為6:00，最晚為12月時(shí)間是7:30；最早由碳匯轉(zhuǎn)變?yōu)樘荚吹氖?月（18:00），最晚為5、6、7、8和9月（19:30）。其中，日固碳時(shí)間最長(zhǎng)的月份為7月，日固碳時(shí)間為13.5 h，最短為12月，日固碳9.5 h。

3.3.2 各季節(jié)碳通量的日變化特征 按照月份的不同，劃分為四季，分別是春季（3—5月）、夏季（6—8月）、秋季（9—11月）、冬季（12—次年2月），分別求同一時(shí)刻下的Fc的平均值。得到碳通量各季節(jié)的日變化（圖8）。各季節(jié)的Fc大致呈現(xiàn)U字形的日變化，在季節(jié)上也會(huì)出現(xiàn)午間的雙峰變化。夏季峰值最早出現(xiàn)，發(fā)生時(shí)間為11:00，春季最不明顯。不同季節(jié)的碳匯峰值排列如下：夏季（?0.83 mg·m?2·s?1）>春季（?0.65 mg·m2·s?1）>秋季（?0.62 mg·m?2·s?1）>冬季（?0.38 mg·m?2·s?1）；日碳匯時(shí)長(zhǎng)夏季（13.5 h）>春季（12 h）=秋季（12 h） > 冬季（10.5 h），可見(jiàn)縉云山針闊葉混交林生態(tài)系統(tǒng)冬季的碳匯能力要小于其他季節(jié)。

圖8 各季節(jié)碳通量的日變化Fig.8 Seasonal variation of carbon flux

3.3.3 碳通量全年月累積變化特征 圖9為凈生態(tài)系統(tǒng)交換量（NEE）的全年月累積變化。全年的月平均NEE為?73.95 g·m?2，年累積NEE 值為-887.40 g·m?2，可知縉云山針闊葉混交林生態(tài)系統(tǒng)為碳匯。不同月份NEE的累積值也不盡相同。最高值出現(xiàn)在 12月份，為 22.21 g·m?2·mon?1；最低值出現(xiàn)在8月份，為?172.27 g·m?2·mon?1。全年所有月份中，只有12月和1月為正值，其他季節(jié)均為負(fù)值，說(shuō)明該生態(tài)系統(tǒng)冬季處于碳源過(guò)程，而其它季節(jié)處于碳匯過(guò)程。碳匯月份之間也存在差異，最大固碳月份為8月（?172.27 g·m?2·mon?1）；最小固碳月份為11月（?3.71 g·m?2·mon?1）。碳匯月份從 2月開(kāi)始，碳匯能力逐漸增強(qiáng)，到8月份到達(dá)峰值；從9月開(kāi)始逐漸減弱，一直到11月降為最低。從植物的生長(zhǎng)季角度來(lái)看，生長(zhǎng)季（4月—10月）與非生長(zhǎng)季之間的固碳能力也存在著較大差異。生長(zhǎng)季的Fc累積變化量為：?826.2 g·m?2，非生長(zhǎng)季則只有：?61.2 g·m?2。這是因?yàn)樵谏L(zhǎng)季節(jié)，植物需要消耗大量的養(yǎng)分，所以，需要進(jìn)行強(qiáng)烈的光合作用來(lái)提供，從而需要吸收大量的CO2。

圖9 碳通量月累積變化Fig.9 monthly cumulative variation of carbon flux

3.3.4 碳通量變化的影響因子 通過(guò)F檢驗(yàn)以及碳通量與主要環(huán)境因子的擬合關(guān)系表明，光合有效輻射和風(fēng)速是影響碳通量變化的主要因子，R2分別為0.75和0.43（圖10）。光輻射作為光合作用的必要條件，其對(duì)植物光合至關(guān)重要。隨著光輻射強(qiáng)度的增加，碳通量的值逐漸減小，植物吸收CO2的能力逐漸增強(qiáng)。風(fēng)速作為影響渦度相關(guān)測(cè)量的主要?dú)庀笠蜃?，在本研究中，?duì)碳通量的影響也很大，風(fēng)速越小碳通量越平穩(wěn)?？諝庀鄬?duì)濕度主要是通過(guò)影響植物的氣孔導(dǎo)度從而間接影響碳通量的變化，空氣濕度過(guò)高時(shí)植物的氣孔會(huì)有不同程度的閉合?？諝鉁囟葘?duì)碳通量變化的影響不顯著。

圖10 碳通量與主要環(huán)境因子的關(guān)系Fig.10 The relationship between carbon flux and major environmental factors

4 討論

4.1 渦度相關(guān)在縉云山針闊葉混交林的適用性

評(píng)判渦度相關(guān)在一個(gè)地區(qū)的適應(yīng)性，主要包括：湍流質(zhì)量的評(píng)價(jià)、能量閉合分子以及通量足跡模擬。

本研究中H、LE、Fc的高質(zhì)量湍流數(shù)據(jù)占比分別為77%、61%、62%。均超過(guò)60%。僅有少量的數(shù)據(jù)需要進(jìn)行剔除插補(bǔ)。湍流的平穩(wěn)狀態(tài)，表明了針闊葉混交林生態(tài)系統(tǒng)近地面湍流發(fā)展平穩(wěn)。王杰帥[14]在之前的研究中利用二次坐標(biāo)旋轉(zhuǎn)所得到的湍流質(zhì)量評(píng)價(jià)分別為：77.5%、65.7%、68.5%，另外比較國(guó)內(nèi)已有研究，吳志祥[17]研究了海南橡膠林生態(tài)系統(tǒng)的湍流數(shù)據(jù)質(zhì)量評(píng)價(jià)分別為：63.4%、52.6%和59.1%與本研究結(jié)果相似。說(shuō)明本研究中湍流趨于理想的穩(wěn)態(tài)條件并且很好的符合湍流觀測(cè)的發(fā)展特點(diǎn)。

關(guān)于能量閉合，張璇[18]、王杰帥[14]等均在縉云山針闊葉混交林生態(tài)系統(tǒng)進(jìn)行過(guò)研究，能量閉合率分別為：0.84、0.78。在國(guó)內(nèi)外相關(guān)的森林生態(tài)系統(tǒng)研究中，吳亞楠[19]關(guān)于小浪底山地人工生態(tài)系統(tǒng)的通量研究中，能量閉合率為0.77；李桐[20]在對(duì)九寨溝針闊葉混交林生態(tài)系統(tǒng)的通量研究中得出能量閉合率為0.71；Y.Kosugi等[21]在研究中，能量閉合率為0.74。均與本研究得到類(lèi)似結(jié)果，說(shuō)明縉云山針闊葉混交林生態(tài)系統(tǒng)能量閉合良好，渦度相關(guān)法在該研究區(qū)具有良好的適用性。

最后是通量足跡模擬，渦度相關(guān)法在測(cè)量時(shí)，需要生態(tài)系統(tǒng)下墊面具有足夠的代表性[22]，由于觀測(cè)方法受到大氣層和下墊面等環(huán)境因素的限制，通量足跡的函數(shù)模擬的貢獻(xiàn)峰值區(qū)主要來(lái)自上風(fēng)向[8,23]。本研究所得到通量足跡模擬的峰值貢獻(xiàn)區(qū)方向與研究期內(nèi)的主風(fēng)向一致，都是東北方向，說(shuō)明通量數(shù)據(jù)信號(hào)來(lái)源正確。本研究中通量足跡的90%累計(jì)貢獻(xiàn)區(qū)年平均距離為225.29 m，相比于其他平緩的林地[24]和農(nóng)田[25-26]生態(tài)系統(tǒng)來(lái)說(shuō)偏小，這是因?yàn)樯降厍鹆甑禺?dāng)中的森林生態(tài)系統(tǒng)的下墊面通常是非均一的，且生態(tài)系統(tǒng)類(lèi)型較為復(fù)雜[18]，從而導(dǎo)致下墊面附近的大氣在摩擦作用下產(chǎn)生環(huán)流，另外不同的觀測(cè)高度也會(huì)使得通量貢獻(xiàn)區(qū)距離不同，通常情況下觀測(cè)高度越高，通量貢獻(xiàn)區(qū)的距離就越大[27]，這也是本研究中通量貢獻(xiàn)區(qū)偏小的原因。

4.2 縉云山針闊葉混交林碳通量的變化

本研究期內(nèi)碳通量的日變化中碳通量的波動(dòng)范圍是?0.97～0.38 mg·m?2·s?1，最小值出現(xiàn)在 8月的12:00，最大值出現(xiàn)在5月的23:00，且波動(dòng)范圍在植物的生長(zhǎng)季更為廣泛，這是由于植物在生長(zhǎng)季需要吸收大量的CO2來(lái)進(jìn)行光合作用，從而合成有機(jī)物供應(yīng)自身的生長(zhǎng)，同時(shí)夜間呼吸也更為旺盛。2016年在對(duì)浙江天目山常綠闊葉落葉混交林生態(tài)系統(tǒng)的檢測(cè)中，所得結(jié)果為：?1.00～0.30 mg·m?2·s?1[28]；2019 年對(duì)廣東鼎湖山的針闊葉混交林的監(jiān)測(cè)結(jié)果中碳通量的變化范圍是?1.18～0.37 mg·m?2·s?1[29]，與本研究結(jié)果相似。說(shuō)明氣候條件和生態(tài)系統(tǒng)類(lèi)型相似的情況下，碳通量的變化也相近。而陳曉峰[30]在對(duì)安吉毛竹林的監(jiān)測(cè)中，所得結(jié)果為?0.52～0.20 mg·m?2·s?1，碳匯能力弱于本研究結(jié)果。說(shuō)明不同類(lèi)型的森林生態(tài)系統(tǒng)，碳匯能力之間存在著較大差異。

本研究中碳通量成雙峰變化，這主要是由于植物的“光合午休”造成的。因?yàn)橹形绻庹諒?qiáng)烈，致使氣溫升高、空氣干燥，植物葉片為減少蒸騰，防止水分蒸發(fā)，被迫關(guān)閉氣孔來(lái)減少水分散失，同時(shí)降低與光合作用有關(guān)的酶活性，這種現(xiàn)象被稱(chēng)為“光合午休”[31]。引起植物光合午休的主要原因是光輻射過(guò)高，過(guò)高的光輻射會(huì)引起葉綠體的回避反應(yīng)，從而使植物的光合作用受到抑制[32]，此外在高溫的脅迫下會(huì)使植物的氣孔關(guān)閉[33]，從而導(dǎo)致其光合能力下降。徐勇峰[34]和牛曉棟[35]在對(duì)洪澤湖濕地生態(tài)系統(tǒng)和天目山常綠闊葉混交林的研究中也得到了相近的結(jié)果，說(shuō)明植物在相應(yīng)的環(huán)境脅迫下會(huì)出現(xiàn)光合午休的現(xiàn)象。光合午休的現(xiàn)象在夏季最為明顯，這是因?yàn)橄募咎?yáng)高度角大，太陽(yáng)輻射強(qiáng)度大，能夠更快地使植物葉片因高溫干旱而關(guān)閉氣孔。春季之所以不明顯，一方面是由于春季太陽(yáng)輻射強(qiáng)度較小，另一方面是因?yàn)槎镜牧值卣羯l(fā)作用不強(qiáng)烈，從而使土壤含水量充足以便于植物在春季吸收，用于蒸騰和生長(zhǎng)。另外，由于本研究中所觀測(cè)的生態(tài)系統(tǒng)為針闊葉混交林生態(tài)系統(tǒng)，植物種類(lèi)并非單一，由于不同植物的光飽和點(diǎn)不同，其吸收CO2所達(dá)到峰值的時(shí)間也不近相同[36-37]。

5 結(jié)論

（1）研究期內(nèi)，湍流因子中的H、LE和Fc的湍流數(shù)據(jù)質(zhì)量評(píng)價(jià)較好，高質(zhì)量數(shù)據(jù)占比分別77%、61%和62%，趨于理想的穩(wěn)態(tài)條件；在生長(zhǎng)季內(nèi)，上、下層能量閉合率良好分別為：0.82和0.73，均在國(guó)際山地森林生態(tài)系統(tǒng)能量閉合認(rèn)可范圍內(nèi)（0.7～0.9）。

（2）研究期內(nèi)，通量足跡模擬高貢獻(xiàn)區(qū)域與研究期內(nèi)主風(fēng)向（東北）一致，表明通量監(jiān)測(cè)的數(shù)據(jù)來(lái)源可信且質(zhì)量良好。

（3）研究期內(nèi)，該生態(tài)系統(tǒng)處于碳匯過(guò)程，年固碳量為887.40 g·m?2，月均日變化和各季節(jié)日變化大多為雙峰變化。光合有效輻射和風(fēng)速是影響碳通量變化的主要因子，R2分別為0.75和0.43。空氣溫度對(duì)碳通量變化的影響不顯著。