原文文，張勁松**，孟 平，同小娟，潘慶梅，何方杰，李劍俠

開路與閉路渦度相關(guān)系統(tǒng)對森林生態(tài)系統(tǒng)甲烷通量觀測的比較*

原文文1，張勁松1**，孟 平1，同小娟2，潘慶梅1，何方杰1，李劍俠3

（1.中國林業(yè)科學(xué)研究院林業(yè)研究所，北京 100091；2.北京林業(yè)大學(xué)，北京 100083；3.國有濟(jì)源市南山林場，河南 454650）

作為渦度相關(guān)技術(shù)觀測的兩種主要技術(shù)手段，開路（OPEC）和閉路（CPEC）兩種渦度相關(guān)系統(tǒng)在觀測森林生態(tài)系統(tǒng)CH4通量過程中存在較大的不確定性，本研究利用OPEC和CPEC兩種渦度相關(guān)系統(tǒng)，對黃河小浪底人工混交林CH4通量進(jìn)行連續(xù)觀測，選取生長旺季2016年7月24日?8月5日連續(xù)14d數(shù)據(jù)，對比兩種觀測系統(tǒng)的功率譜和協(xié)譜，估算閉路渦度相關(guān)系統(tǒng)的延遲時間，并分析其在連續(xù)晴天和連續(xù)雨天的通量觀測結(jié)果。結(jié)果表明：CPEC系統(tǒng)的功率譜和協(xié)譜在所有頻率上與OPEC系統(tǒng)基本一致，在慣性副區(qū)功率譜符合?5/3相似規(guī)律，協(xié)譜符合?4/3相似規(guī)律；以O(shè)PEC系統(tǒng)為“準(zhǔn)”標(biāo)準(zhǔn)，CPEC系統(tǒng)觀測CH4通量的延遲時間合適流速范圍內(nèi)的5個不同流速（40、37.5、35.5、33.5、31.5L·min?1）分別為4.6、7.7、5.3、10.8和14.3s，平均延遲時間為8～9s；與OPEC觀測系統(tǒng)測定的CH4通量相比，CPEC系統(tǒng)觀測結(jié)果晴天偏低12%；雨天高出32%。OPEC觀測系統(tǒng)適用于晴天CH4通量觀測。經(jīng)校正，消除延遲影響后的CPEC觀測系統(tǒng)可用于測定雨天CH4通量，以彌補(bǔ)OPEC觀測系統(tǒng)缺測的值。兩種系統(tǒng)并行觀測、相互彌補(bǔ)，可望獲得更完整、更高質(zhì)量的CH4通量數(shù)據(jù)。

渦度相關(guān)技術(shù)；延遲時間；譜分析；甲烷通量；森林生態(tài)系統(tǒng)

甲烷（CH4）是一種痕量溫室氣體，其增溫潛勢是CO2的28倍[1]，引起各國政府和科學(xué)界的廣泛重視。由于CH4源/匯復(fù)雜，全球CH4收支評估依然存在較大不確定性[2?5]，準(zhǔn)確測量CH4通量成為相關(guān)研究關(guān)注的重點(diǎn)問題。目前渦度相關(guān)方法觀測CO2、H2O以及熱量通量的研究已經(jīng)逐步成熟[6?8]，也已運(yùn)用到CH4通量的觀測中[9]。隨著高頻響應(yīng)的CH4濃度分析儀的出現(xiàn)，基于渦度相關(guān)法的CH4通量長期連續(xù)觀測研究在國際上日漸增加[10?12]。

根據(jù)氣體濃度傳感器的性能不同，基于渦度相關(guān)法觀測通量的分析儀分為開路式（OPEC）和閉路式（CPEC）兩大類。截至2017年底，全球不同地區(qū)已建立了900多個通量觀測站點(diǎn)（http://fluxnet. fluxdata.org/sites/site-summary），目前在國際通量觀測網(wǎng)的觀測中，OPEC和CPEC兩種觀測系統(tǒng)成為通量的主要觀測手段，AsiaFlux通量網(wǎng)主要以CPEC為主，其中有幾個觀測站點(diǎn)用OPEC與CPEC兩種觀測系統(tǒng)并行觀測，Ameri-Flux通量網(wǎng)是兩種觀測系統(tǒng)并行存在，Euro-Flux通量網(wǎng)以CPEC為主，Ko-Flux通量網(wǎng)主要以O(shè)PEC為主。兩種渦度相關(guān)系統(tǒng)并行觀測也可以互相彌補(bǔ)數(shù)據(jù)缺失問題，但是，系統(tǒng)本身的差異，可能導(dǎo)致觀測結(jié)果的不同，同時，兩種不同觀測系統(tǒng)觀測結(jié)果的不確定性研究也一直都是一個重要的研究課題。目前，基于渦度相關(guān)法CO2通量觀測及其數(shù)據(jù)處理過程中不確定性分析已逐漸成熟[13?16]，但是其不確定性是否等同于CH4通量，也成為相關(guān)研究人員廣泛關(guān)注的熱點(diǎn)問題之一[17?18]。目前，大部分研究主要集中在CH4通量的時空變化特征分析和各系統(tǒng)本身的不確定性分析，國內(nèi)目前僅見宋霞等關(guān)于兩種系統(tǒng)觀測CO2通量的對比研究報道[19]，國外部分研究大多集中于濕地、水稻田等生態(tài)系統(tǒng)[20?22]，而對于兩種不同觀測系統(tǒng)觀測森林生態(tài)系統(tǒng)CH4通量不確定性分析鮮少報道。因此，對比分析OPEC與CPEC兩種觀測系統(tǒng)的長期CH4通量觀測結(jié)果，可以為中國森林生態(tài)系統(tǒng)的CH4通量觀測提供適宜的觀測手段。

OPEC和CPEC觀測系統(tǒng)在實(shí)際應(yīng)用中的優(yōu)缺點(diǎn)主要表現(xiàn)在對觀測環(huán)境的適應(yīng)性、設(shè)備維護(hù)和觀測結(jié)果等方面的差異[17]。一般認(rèn)為OPEC以高頻率響應(yīng)為主要優(yōu)勢[23]，功耗低，基本測算原理符合大氣湍流運(yùn)動特征；對信號強(qiáng)度要求較高，在雨天或空氣質(zhì)量較差的條件下，數(shù)據(jù)無法使用，若希望得到高質(zhì)量的通量數(shù)據(jù)，需要定期清潔傳感器鏡頭[24?26]。而CPEC恰好彌補(bǔ)了OPEC在實(shí)際應(yīng)用中所存在的缺陷，CPEC相對較穩(wěn)定，不容易受到外界環(huán)境的干擾，并且可以在弱湍流通量條件下全天候?qū)ι稚鷳B(tài)系統(tǒng)進(jìn)行長期穩(wěn)定的CH4通量觀測[27?30]。為了滿足高頻率數(shù)據(jù)采集的要求，CPEC的傳感器需要用高流量泵對采樣氣體抽樣，而野外電力限制以及站點(diǎn)配置有時難以保證其對供電功率和供電質(zhì)量的要求，這在一定程度上限制了CPEC系統(tǒng)在野外的應(yīng)用[30]。同時由于抽氣管進(jìn)口和氣體分析儀之間有一定的距離，導(dǎo)致CPEC系統(tǒng)的紅外分析儀測定氣體濃度的時間滯后于超聲風(fēng)速儀測定風(fēng)速的時間，即時間延遲問題。因此，精確估測CPEC系統(tǒng)紅外分析儀的延遲時間是準(zhǔn)確計算通量的前提[31]，CO2的延遲時間在7～8s[19]，而關(guān)于森林生態(tài)系統(tǒng)CH4通量觀測的延遲時間鮮少報道，因此，對兩個系統(tǒng)進(jìn)行對比分析，估算CPEC系統(tǒng)的延遲時間，可為準(zhǔn)確觀測森林生態(tài)系統(tǒng)CH4通量提供理論基礎(chǔ)，并為準(zhǔn)確計算全球源/匯提供科學(xué)基礎(chǔ)。

第八次全國森林資源清查結(jié)果顯示，中國現(xiàn)有森林面積2.08億hm2，其中人工林面積為0.69億hm2，占森林總面積的33.2%。研究表明，中國人工林對全國森林碳匯的貢獻(xiàn)率超過80%[32]。全球尺度上由于森林生態(tài)系統(tǒng)CH4估算不完整或估算偏差造成的丟失源大約為10Tg·a?1[33]。目前大部分的CH4通量研究多集中在濕地、水稻田等生態(tài)系統(tǒng)[34?36]，而關(guān)于森林生態(tài)系統(tǒng)CH4通量的研究鮮少報道。華北山區(qū)地處暖溫帶氣候區(qū)，一直是中國林業(yè)生態(tài)工程建設(shè)的重點(diǎn)區(qū)域，迫切需要準(zhǔn)備估算該森林CH4通量，闡明其在全球氣候變暖中的作用。因此，需要長期準(zhǔn)確測算森林生態(tài)系統(tǒng)CH4通量，了解源匯轉(zhuǎn)換格局，為評估與預(yù)測CH4收支狀況提供數(shù)據(jù)基礎(chǔ)。為此，本研究以華北南部丘陵山地栓皮櫟?側(cè)柏?刺槐人工混交林生態(tài)系統(tǒng)為例，采用OPEC與CPEC兩種觀測系統(tǒng)長期連續(xù)觀測CH4通量，通過實(shí)時觀測CH4通量數(shù)據(jù)的譜分析比較OPEC與CPEC兩種觀測系統(tǒng)的性能，評價兩種觀測系統(tǒng)觀測CH4通量數(shù)據(jù)質(zhì)量及觀測結(jié)果的可靠性，確定CPEC觀測CH4通量相對于OPEC系統(tǒng)的延遲時間；并通過半小時CH4平均通量結(jié)果的比較，分析兩種觀測系統(tǒng)觀測結(jié)果的差異，以期為中國不同森林生態(tài)系統(tǒng)的CH4通量觀測選取適宜的觀測手段提供科學(xué)依據(jù)，同時為兩種渦度相關(guān)系統(tǒng)并行觀測以互相彌補(bǔ)數(shù)據(jù)缺失問題提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究在國家林業(yè)局黃河小浪底森林生態(tài)定位研究站進(jìn)行，該站位于河南省濟(jì)源市（35°01'N，112°28'E，海拔410m），屬暖溫帶亞濕潤季風(fēng)氣候。年平均氣溫12.4℃，年平均降水量641.7mm。受季風(fēng)氣候影響，降水季節(jié)分配不均勻，5?9月平均降水量438.0mm，占全年的68.3%。生長季風(fēng)向以東北偏東、西南為主。主要樹種為栓皮櫟（）、側(cè)柏（）、刺槐（），林齡分別為45、43和41a，平均株高分別為10.5、8.2和9.3m。生長季植被葉面積指數(shù)（LAI）最大為6.3。通量觀測塔周圍1.8 km2范圍內(nèi)平均坡度為14°，林木覆蓋率約為96%，栓皮櫟、側(cè)柏、刺槐所占比例分別為80%、8%和12%[37]。

1.2 渦度相關(guān)觀測系統(tǒng)

渦度相關(guān)系統(tǒng)包括開路式渦度相關(guān)系統(tǒng)（OPEC）和閉路式渦度相關(guān)系統(tǒng)（CPEC）。

開路式渦度相關(guān)系統(tǒng)（OPEC）主要由CSAT3型三維超聲風(fēng)速儀（Campbel Sci.Inc.，USA）、快速響應(yīng)的LI-7500型紅外CO2/H2O分析儀（Li-Cor.Inc.，USA）和LI-7700型紅外CH4/H2O分析儀（Li-Cor.Inc.，USA）組成，安裝高度距地面36m。上述儀器均與CR3000型數(shù)據(jù)采集器（Campbell Sci.Inc.，USA）相連，晝夜連續(xù)自動采集，原始數(shù)據(jù)采樣頻率為10 Hz，每30 min輸出1組平均值。

閉路式渦度相關(guān)系統(tǒng)（CPEC）與OPEC系統(tǒng)最大區(qū)別是需要用220 V大容量真空泵把氣體抽入分析儀，氣泵流量為5m3·h?1，進(jìn)氣口在開路分析儀附近。進(jìn)氣管道長40m，內(nèi)徑為5mm，抽氣管道共接有4個過濾芯，分別是：（i）進(jìn)氣口處有一直徑為100μm的過濾芯，（ii）管道末端進(jìn)口處有一直徑5μm和一直徑2μm的梅花形過濾芯，（iii）在分析儀進(jìn)氣口處有一直徑為2μm的金屬過濾芯，主要過濾空氣中顆粒物等雜質(zhì)，保持分析儀內(nèi)光腔清潔（衰蕩時間μ > 6μm·s?1）。

1.3 冠層湍流譜

湍流是由各種渦旋疊加而成，廣泛存在于大氣中，在大氣邊界層中觀測到的渦旋是時間和空間上的廣闊結(jié)構(gòu)。譜分析反映了這些渦旋對湍流總動能的貢獻(xiàn)。由于湍流具有普遍統(tǒng)計規(guī)律的假設(shè)[24]，其慣性副區(qū)內(nèi)的風(fēng)速u分量譜的表達(dá)式為

使用Welch方法求解三維方向風(fēng)速（u、v、w），CH4濃度的功率譜Sx(f)和協(xié)譜Cwx(f)，對功率譜和協(xié)譜都分別利用變量x的方差以及垂直風(fēng)速w和x的協(xié)方差進(jìn)行歸一化處理，橫坐標(biāo)為習(xí)慣的無量綱頻率。再通過對數(shù)坐標(biāo)下的功率譜和協(xié)譜點(diǎn)陣曲線的分析，比較點(diǎn)陣曲線斜率是否滿足?5/3和?4/3[40]，以確定不同變量的功率譜在慣性副區(qū)內(nèi)的斜率，就可以把握渦度相關(guān)系統(tǒng)儀器對高低頻信號的響應(yīng)能力。

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

利用Matlab2016b對原始觀測數(shù)據(jù)（10Hz）進(jìn)行湍流譜計算，根據(jù)計算結(jié)果采用Origin2018軟件進(jìn)行繪制湍流譜。利用Eddypro軟件對采樣頻率為10 Hz的原始湍流數(shù)據(jù)進(jìn)行再處理。數(shù)據(jù)處理過程包括異常值（含野點(diǎn)）去除、延遲時間校正、坐標(biāo)旋轉(zhuǎn)、頻率損失訂正、感熱通量超聲虛位溫修正以及WPL修正等必要的修正[41?45]。同時，剔除了夜間摩擦速度＜0.35m·s?1的通量數(shù)據(jù)[29]。當(dāng)儀器出現(xiàn)故障、人為影響以及雨天或清晨有露水時，所得的異常值應(yīng)當(dāng)剔除。對以上經(jīng)過初步篩選的數(shù)據(jù)計算出各月平均值和方差，與平均值相差超過3倍方差的數(shù)值也要剔除掉。小于2h的缺失數(shù)據(jù)用線性內(nèi)插法插補(bǔ)，大于2h的缺失數(shù)據(jù)采用平均日變化法（MDV）進(jìn)行插補(bǔ)[46]。

2 結(jié)果與分析

2.1 兩種觀測系統(tǒng)CH4通量觀測數(shù)據(jù)的頻譜分析

選擇典型晴天7月24日大氣層結(jié)條件下半小時原始數(shù)據(jù)，對OPEC和CPEC系統(tǒng)進(jìn)行功率譜和協(xié)譜分析。從圖1功率譜a、b可以看出，在對數(shù)坐標(biāo)條件下，OPEC和CPEC系統(tǒng)的CH4的功率譜圖低頻區(qū)有明顯峰值，在高頻區(qū)有較清楚的慣性副區(qū)，且在慣性副區(qū)內(nèi)斜線斜率基本符合?5/3定律[47]。功率譜高頻端出現(xiàn)“上翹”現(xiàn)象，說明兩種系統(tǒng)的儀器均具有高頻噪音，可能與天氣、地形等狀況有關(guān)。且上翹曲線下面積很小，對湍流通量測定沒有明顯的影響。以O(shè)PEC系統(tǒng)觀測為“準(zhǔn)”標(biāo)準(zhǔn)，在所有頻率范圍內(nèi)CPEC與OPEC系統(tǒng)分析儀的頻率響應(yīng)基本一致，但當(dāng)頻率＞1時，CPEC系統(tǒng)觀測的CH4氣體功率譜出現(xiàn)衰減現(xiàn)象，可能是因?yàn)镃PEC采樣管道具有高頻過濾的作用，這與CPEC系統(tǒng)觀測CO2氣體功率譜研究結(jié)果相似[48?49]。CPEC與OPEC系統(tǒng)觀測的CH4協(xié)譜對比如圖1c、d所示，由圖可以看出，CPEC系統(tǒng)LGR-CH4氣體快速分析儀與OPEC系統(tǒng)Li-7700 CH4/H2O紅外分析儀所觀測的CH4協(xié)譜在慣性副區(qū)內(nèi)與?4/3斜線基本一致[47]，由此可見，CPEC與OPEC系統(tǒng)分析儀的頻率響應(yīng)基本一致，CPEC系統(tǒng)的LGR-CH4氣體快速分析儀對CH4高頻信號的響應(yīng)能力是能夠滿足通量觀測要求的。

（a）CPEC系統(tǒng)功率譜Spectra of CPEC; （b）OPEC系統(tǒng)功率譜Spectra of OPEC; （c）CPEC系統(tǒng)垂直風(fēng)速和CH4的協(xié)譜 Cospectrum of CPEC;（d）OPEC系統(tǒng)垂直風(fēng)速和CH4的協(xié)譜 Cospectrum of OPEC

注：觀測時間為2016年7月24日，晴天。

Note: The observed time is in 9：00?16：00, July 24th in 2016, sunny day.

2.2 CPEC系統(tǒng)紅外分析儀的延遲時間

時間延遲是影響CPEC系統(tǒng)準(zhǔn)確測量CH4通量的重要因素，一般通過垂直風(fēng)速與CH4濃度最大協(xié)方差來決定[23，41，50?53]，通常在湍流運(yùn)動劇烈時，所計算的延遲時間較可靠。延遲時間的長短主要受氣流速度、管道長度和管道內(nèi)徑影響。因此，本研究在2016年7月24日?8月1日分別設(shè)置5個不同采樣流速（流速40L·min?1，14：00?14：30；流速37.5L·min?1，12：00?12：30；流速35.5L·min?1，10：30?11：00；流速33.5L·min?1，13：00?13：30；流速31.5L·min?1，11：30?12：00）的半小時原始數(shù)據(jù)，以O(shè)PEC系統(tǒng)為準(zhǔn)標(biāo)準(zhǔn)[9]，計算垂直風(fēng)速與CH4濃度最大協(xié)方差來確定延遲時間（圖2）。當(dāng)流速為40L·min?1時，CPEC系統(tǒng)的垂直風(fēng)速與CH4濃度最大協(xié)方差值出現(xiàn)的時間在0.13s，此時OPEC系統(tǒng)最大值出現(xiàn)在4.74s，此流速下計算的延遲時間為4.6s；同上相應(yīng)的流速設(shè)置下其延遲時間分別為7.7、5.3、10.8和14.3s，平均延遲時間在8～9s。因此，以O(shè)PEC系統(tǒng)的觀測結(jié)果為準(zhǔn)標(biāo)準(zhǔn)，CPEC系統(tǒng)測定CH4通量的延遲平均時間在8～9s，比CO2的延遲時間（7～8s）略長[19]，主要是因?yàn)镃H4是一種痕量溫室氣體，森林生態(tài)系統(tǒng)中含量很少，并且，流速越小滯后時間越長。

圖2 CPEC系統(tǒng)在不同流速下相對于OPEC系統(tǒng)最大協(xié)方差決定延遲時間

2.3 CPEC與OPEC兩種觀測系統(tǒng)所測CH4通量的比較

基于白天半小時通量數(shù)據(jù)，分別對典型晴天（2016年7月22?25日）和雨天（2016年8月2?5日）CPEC與OPEC系統(tǒng)CH4通量觀測結(jié)果進(jìn)行對比分析，結(jié)果見圖3。由圖3a可見，連續(xù)晴天，兩套系統(tǒng)觀測CH4通量模擬直線斜率為0.87（R2=0.8336，P<0.01），觀測結(jié)果差異不大，與OPEC系統(tǒng)相比，CPEC系統(tǒng)所測CH4通量偏低。但是連續(xù)雨天（圖3b），OPEC系統(tǒng)觀測結(jié)果低于CPEC系統(tǒng)（R2=0.4024，P<0.01）。

對比分析兩套系統(tǒng)的觀測結(jié)果可以看出，雖然CPEC系統(tǒng)外接抽氣管道的衰減作用造成觀測結(jié)果較低于OPEC系統(tǒng)，但是均可滿足通量觀測的需求，觀測結(jié)果可靠。降雨對OPEC系統(tǒng)觀測信號造成影響，導(dǎo)致其觀測通量值低于正常值。因此，連續(xù)雨天OPEC觀測結(jié)果不可靠，但可以通過CPEC觀測結(jié)果進(jìn)行訂正或直接插補(bǔ)。

（a）7月22?25日（連續(xù)晴天）Continuous sunny day (July 22nd to 25th)；（b）8月2?5日（連續(xù)雨天）Continuous rainy day (Aug. 2nd to 5th)

3 結(jié)論與討論

3.1 討論

在評價CPEC系統(tǒng)觀測結(jié)果時，通常假定OPEC系統(tǒng)的觀測通量結(jié)果為比較標(biāo)準(zhǔn)[9]。為了確定CPEC系統(tǒng)觀測的可靠性，對CPEC系統(tǒng)和OPEC系統(tǒng)的功率譜和頻譜進(jìn)行對比分析，結(jié)果表明，在對數(shù)坐標(biāo)下，CPEC系統(tǒng)和OPEC系統(tǒng)在所有頻率范圍內(nèi)頻率響應(yīng)基本一致。低頻區(qū)有明顯峰值，在高頻區(qū)有較清楚的慣性副區(qū)且在慣性副區(qū)功率譜滿足?5/3定律，協(xié)譜滿足?4/3規(guī)律，在高頻端均出現(xiàn)“上翹”現(xiàn)象。這兩套系統(tǒng)儀器均存在高頻噪音，可能與天氣、地形等狀況有關(guān)，上翹曲線下面積很小，且高頻噪音與垂直風(fēng)速無相關(guān)性，因此，對湍流通量測定沒有明顯的影響[44]。當(dāng)頻率大于1Hz時，CPEC系統(tǒng)的功率譜和頻譜均小于OPEC系統(tǒng)，這可能是因?yàn)镃PEC系統(tǒng)采樣管道具有高頻過濾的作用，導(dǎo)致氣體濃度在管道內(nèi)發(fā)生衰減。這與兩套系統(tǒng)觀測CO2通量的譜分析結(jié)果一致[19]。同時，由于本研究觀測區(qū)位于華北低丘山地，下墊面地形起伏且植被非均勻，因此湍流譜分布并不完全等同于理想條件下譜分布。因此，CPEC系統(tǒng)的LGR-CH4氣體快速分析儀對CH4高頻信號的響應(yīng)能力是能夠滿足通量觀測要求的，觀測結(jié)果可靠。

用CPEC系統(tǒng)連續(xù)觀測時，由于抽氣管進(jìn)口和氣體分析儀之間有一定距離，因此CPEC系統(tǒng)的分析儀測定氣體濃度的時間滯后于超聲風(fēng)速儀測定風(fēng)速的時間，精確估測時間延遲成為CPEC系統(tǒng)準(zhǔn)確觀測CH4通量的重要前提。而延遲時間的長短主要受氣流速度、管道長度、管道內(nèi)徑的影響。以O(shè)PEC系統(tǒng)觀測結(jié)果為比較標(biāo)準(zhǔn)，本研究在合適采樣流速范圍內(nèi)選取5個不同流速（40、37.5、35.5、33.5、31.5L·min?1）計算的延遲時間分別為4.6、7.7、5.3、10.8和14.3s，平均延遲時間為8～9s，比CO2的延遲時間（7～8s）略長，主要是因?yàn)镃H4是一種痕量溫室氣體，在森林生態(tài)系統(tǒng)中含量很低。并且，流速越小，延遲時間越長。因此，在后續(xù)研究中，只要在CH4通量計算過程中給予時間延遲的訂正，CPEC系統(tǒng)的時間延遲不會造成明顯的觀測誤差。

兩種渦度相關(guān)系統(tǒng)并行觀測也可以互相彌補(bǔ)數(shù)據(jù)缺失問題，但是，系統(tǒng)本身的差異，可能導(dǎo)致觀測結(jié)果的不同。因此，本研究對比分析兩種渦度相關(guān)系統(tǒng)在連續(xù)晴天和雨天觀測的通量結(jié)果。晴天條件下，兩種系統(tǒng)的觀測結(jié)果基本一致，CPEC系統(tǒng)對CH4通量的觀測結(jié)果是OPEC系統(tǒng)的88%，這可能是CPEC系統(tǒng)抽樣管道的衰減，其觀測值低于OPEC系統(tǒng)觀測結(jié)果，但是均可滿足通量觀測的需求，觀測結(jié)果可靠。但是在連續(xù)降雨條件下，兩種系統(tǒng)的觀測結(jié)果不一致，CPEC系統(tǒng)的觀測結(jié)果高于OPEC系統(tǒng)。這是由于連續(xù)降雨條件下，OPEC系統(tǒng)LI-7700分析儀傳感器鏡面清潔度降低，導(dǎo)致捕捉到的CH4信號值逐漸降低，觀測的通量數(shù)據(jù)低于真實(shí)值?？梢姡琌PEC系統(tǒng)受連續(xù)雨天等惡劣天氣條件影響較大。因此，在今后的通量觀測試驗(yàn)中，兩種系統(tǒng)可以作為并行觀測相互彌補(bǔ)的觀測技術(shù)手段。以O(shè)PEC觀測的通量數(shù)據(jù)為主，連續(xù)雨天采用CPEC系統(tǒng)觀測數(shù)據(jù)進(jìn)行訂正或者直接插補(bǔ)。

3.2 結(jié)論

（1）對比CPEC和OPEC系統(tǒng)觀測CH4通量的譜分析，實(shí)時數(shù)據(jù)譜分析結(jié)果表明兩種系統(tǒng)觀測結(jié)果一致，功率譜和頻譜在慣性副區(qū)分別滿足斜率為?5/3和?4/3的規(guī)律，在高頻段雖有“上翹”現(xiàn)象，但不影響通量觀測結(jié)果。

（2）對合適流速范圍內(nèi)的5個不同流速（40、37.5、35.5、33.5、31.5L·min?1）計算CPEC系統(tǒng)觀測CH4通量的延遲時間分別為4.6、7.7、5.3、10.8和14.3s，平均延遲時間為8～9s。

（3）分析不同天氣條件下兩種系統(tǒng)觀測CH4通量結(jié)果表明，與OPEC觀測系統(tǒng)測定的CH4通量相比，CPEC系統(tǒng)觀測結(jié)果晴天偏低12%；雨天高出32%。OPEC觀測系統(tǒng)適用于晴天CH4通量觀測。經(jīng)校正，消除延遲影響后的CPEC觀測系統(tǒng)可用于測定雨天CH4通量，以彌補(bǔ)OPEC觀測系統(tǒng)缺測的值。兩種系統(tǒng)并行觀測、相互彌補(bǔ)，可望獲得更完整、更高質(zhì)量的CH4通量數(shù)據(jù)。

[1] IPCC.Climate change 2013:the physical science basic. Contribution of working group I to the fifth assessment report of the intergovernmental panel on climate change[J]. Cambridge:Cambridge University,2013.

[2] Bridgham S D,Cadillo-Quiroz H,Keller J K,et al.Methane emissions from wetlands:biogeochemical,microbial,and modeling perspectives from local to global scales[J].Global Change Biology,2013,19(5):1325-1346.

[3] Kai F M,Tyler S C,Randerson J T,et al.Reduced methane growth rate explained by decreased Northern Hemisphere microbial sources[J].Nature,2011,476(7359):194-197.

[4] Bousquet P,Ciais P,Miller J B,et al.Contribution of anthropogenic and natural sources to atmospheric methane variability[J]. Nature,2006,443(7110):439-443.

[5] Solomon E A,Kastner M,Macdonald I R,et al.Considerable methane fluxes to the atmosphere from hydrocarbon seeps in the Gulf of Mexico[J].Nature Geoscience,2009,2(8):561-565.

[6] 蘇榮瑞,劉凱文,耿一風(fēng),等.江漢平原稻-油連作系統(tǒng)冠層CO2通量變化特征[J].中國農(nóng)業(yè)氣象,2012,33(3):362-367.

Su R R,Liu K W,Geng Y F,et al.CO2flux variation over canopy rice-rape succession system in Jianghan Plain[J]. Chinese Journal of Agrometeorology,2012,33(3):362-367.(in Chinese)

[7] 賈志軍,姬興杰.三江平原稻田蒸散量模擬研究[J].中國農(nóng)業(yè)氣象,2014,35(4):380-388.

Jia Z J,Ji X J.Simulation study on evapotranspiration of rice paddy in the Sanjiang Plain[J].Chinese Journal of Agrometeorology, 2014,35(4):380-388.(in Chinese)

[8] 劉渡,李俊,同小娟,等.華北平原冬小麥/夏玉米輪作田能量閉合狀況分析[J].中國農(nóng)業(yè)氣象,2012,33(4):493-499.

Liu D,Li J,Tong X J,et al.Analysis of the energy balance closure in a winter wheat/summer maize double cropping system in the North China Plain[J].Chinese Journal of Agrometeorology,2012,33(4):493-499.(in Chinese)

[9] Baldocchi D,Detto M,Sonnentag O,et al.The challenges of measuring methane fluxes and concentrations over a peatland pasture[J].Agricultural and Forest Meteorology,2012,153: 180-187.

[10] Olson D M,Griffis T J,Noormets A,et al. Interannual, seasonal,and retrospective analysis of the methane and carbon dioxide budgets of a temperate peatland[J].Journal of Geophysical Research:Biogeosciences,2013,118(1):226-238.

[11] Iwata H,Kosugi Y,Ono K,et al.Cross-Validation of open-path and closed-path eddy-covariance techniques for observing methane fluxes[J].Boundary-Layer Meteorology,2014,151 (1):95-118.

[12] Yu X,Song C,Sun L,et al.Growing season methane emissions from a permafrost peatland of northeast China: observations using open-path eddy covariance method[J]. Atmospheric Environment,2017,153:135-149.

[13] Masahito U,Kota Y,Kentaro T.A cool-temperate young larch plantation as a net methane source:a 4-year continuous hyperbolic relaxed eddy accumulation and chamber measurements[J].Atmospheric Environment,2018,184:110-120.

[14] Wang G X,Hirata R,Hirano T,et al.Micrometeorological measurement of methane flux above a tropical peat swamp forest[J].Agricultural and Forest Meteorology,2018,256-257: 353-361.

[15] 胡誠,張彌,肖薇,等.通量及其不確定性對農(nóng)業(yè)區(qū)高塔CO2濃度模擬的影響[J].中國農(nóng)業(yè)氣象,2017,38(8):469-480.

Hu C,Zhang M,Xiao W,et al.Effect of flux and its uncertainty on tall tower CO2concentration simulation in the agricultural domain[J].Chinese Journal of Agrometeorology, 2017,38(8):469-480.(in Chinese)

[16] 王春林,周國逸,王旭,等.復(fù)雜地形條件下渦度相關(guān)法通量測定修正方法分析[J].中國農(nóng)業(yè)氣象,2007,28(3):233-240.

Wang C L,Zhou G Y,Wang X,et al.Analysis of correction method on eddy flux measurement over complex terrain[J].Chinese Journal of Agrometeorology,2007,28(3): 233-240.(in Chinese)

[17] Kroon P S,Vesala T,Grace J.Flux measurements of CH4and N2O exchanges[J].Agricultural and Forest Meteorology, 2010,150(6):740-747.

[18] Eugster W,Zeyer K,Zeeman M,et al.Methodical study of nitrous oxide eddy covariance measurements using quantum cascade laser spectrometery over a Swiss forest[J]. Biogeosciences,2007,4(5):927-939.

[19] 宋霞,于貴瑞,劉允芬,等.開路與閉路渦度相關(guān)系統(tǒng)通量觀測比較研究[J].中國科學(xué):地球科學(xué),2004,34(S2): 67-76.

Song X,Yu G R,Liu Y F,et al.Comparison of flux measurement by open-path and close-path eddy covariance systems[J].Science in China:Earth Science,2004,34(S2): 67-76.(in Chinese)

[20] Ayaka S,Masayuki I,Takashi H,et al.Ecosystem-scale methane flux in tropical peat swamp forest in Indonesia [J].Global Change Biology,2018:1-14.

[21] Taylor M A,Celis G,Ledman J D,et al.Methane efflux measured by eddy covariance in Alaskan upland tundra undergoing permafrost degradation[J].Journal of Geophysical Research:Biogeosciences,2018,123(7):2695-2710.

[22] Desjardins R L,Worth D E,Pattey E,et al.The challenge of reconciling bottom-up agricultural methane emissions inventories with top-down measurements[J].Agricultural and Forest Meteorology,2018,48:48-59.

[23] Leuning R,Moncrieff J.Eddy-covariance CO2flux measurements using open and close-path CO2analyzers: corrections for analyzer water vapor sensitivity and damping of fluctuations in air sampling tubes[J].Boundary-Layer Meteorology,1990,53:63-76.

[24] Gharavi M,Buckley S G.Diode laser absorption spectroscopy measurement of lines strengths and pressure broadening coefficients of the methane 2×3 band at elevated temperatures[J].Mol Spectrosc,2005,229:78-88.

[25] Li-Cor Inc.Li-7700 open path CH4analyzer instruction manual[J].2010,LI-COR Inc.,Lincoln:1-8.

[26] Detto M,Katul G G.Simplified expressions for adjusting higher-order turbulent statistics obtained from open path gas analyzers[J].Boundary-Layer Meteorology,2007,122(1):205-216.

[27] Desai A R,Xu K,Tian H,et al.Landscape-level terrestrial methane flux observed from a very tall tower[J].Agricultural and Forest Meteorology,2015,201:61-75.

[28] Iwata H,Harazono Y,Ueyema M,et al.Methane exchange in a poorly-drained black spruce forest over permafrost observed using the eddy covariance technique[J].Agricultural and Forest Meteorology,2015,214-215:157-168.

[29] Asakawa T,Kanno N,Tonokura K.Diode laser detection of greenhouse gases in the near-Infrared region by wavelength modulation spectroscopy:pressure dependence of the detection sensitivity[J].Sensors,2010,10:4686-4699.

[30] Detto M,Verfaillie J,Anderson F,et al.Comparing laser-base open and close path gas analyzers to measure methane fluxes using the eddy covariance method[J].Agricultural and Forest Meteorology,2011,151:1312-1324.

[31] Kroon P S,Hensen A,Hjj J,et al.Suitability of quantum cascade laser spectrometry for CH4and N2O eddy covariance measurements[J].Biogeosciences Discussions, 2007,4(2):1137-1165.

[32] Fang J Y,Liu G H,Zhu B,et al.Carbon budgets of three temperate forest ecosystems in Dongling Mt[J].Science in China,2007,50(1):92-101.

[33] Megonigal J P,Guenther A B.Methane emissions from upland forest soils and vegetation[J].Tree Physiology, 2008,28:491-498.

[34] 何方杰,韓輝邦,馬學(xué)謙,等.隆寶灘沼澤濕地不同區(qū)域的甲烷通量特征及影響因素[J].生態(tài)環(huán)境學(xué)報,2019,8(4): 803-811.

He F J,Han H B,Ma X Q,et al.Characteristics and influence factors of CH4flux in different areas of Longbaotan marsh wetland[J].Ecology and Environmental Sciences,2019,8(4): 803-811.(in Chinese)

[35] 奉小明,王凱,鄭循化,等.亞熱帶地區(qū)蔬菜地甲烷凈交換通量研究[J].中國生態(tài)農(nóng)業(yè)學(xué)報,2018,26(8):1091-1099.

Feng X M,Wang K,Zheng X H,et al.Net methane flux exchange in subtropical vegetable fields[J].Chinese Journal of Eco-agriculture,2018,26(8):1091-1099.(in Chinese)

[36] 張強(qiáng),蔣國慶,孫睿,等.張掖濕地甲烷通量動態(tài)特征及其影響因子[J].生態(tài)學(xué)報,2017,37(17):5681-5690.

Zhang Q,Jiang G Q,Sun R,et al.CH4flux variations and main factors from a reed wetland oasis-desert area in Zhangye,China[J].Ecologica Scinica,2017,37(17):5681-5690. (in Chinese)

[37] Tong X,Ping M,Zhang J,et al.Ecosystem carbon exchange over a warm-temperate mixed plantation in the Lithoid hilly area of the north China[J].Atmospheric Environment,2012, 49(3):257-267.

[38] H?gstr?m Ulf.Analysis of turbulence structure in the surface layer with a modified similarity formulation for near neutral conditions[J].Journal of the Atmospheric Sciences,1990,47 (16):1949-1972.

[39] Garratt J R,Dessler A J,Houghton J T,et al.The atmospheric boundary layer[M].Cambridge:Cambridge University Press, 1992:316.

[40] Kolmogorov A N.The local structure of turbulence in an incompresible fluid at very high reynolds number[J].Soviet Physics Uspekhi,1941,30:301-305.

[41] Vickers D,Mahrt L.Quality control and flux sampling problems for tower and aircraft data[J].Atmos.Oceanic Technol,1997,14:512-526.

[42] Wilczak J M,Oncley S P,Stage S A.Sonic anemometer tilt correction algorithms[J].Boundary Layer Meteorology,2001, 99:127-150.

[43] Ibrom A,Dellwik E,Larsen S E,et al.On the use of the Webb-Pearman-Leuning theory for closed-path eddy correlation measurements[J].Tellus B.,2007,59:937-946.

[44] Ibrom A,Dellwik E,Flvbjerg H,et al.Strong low-pass filtering effects on water vapor flux measurements with closed-path eddy correlation system[J].Agricultural and Forest Meteorology,2007,147:140-156.

[45] Mauder M,Foken T.Documentation and instruction manual of the eddy covariacne software package TK2[J]. Arbeitsergebn,Unversity of Bayreuth,Dept.of Micrometeorology, 2004,26:1-60.

[46] Long K D,Flanagan L B,Cai T.Diurnal and seasonal variation in methane emissions in a northern Canadian peatland measured by eddy covariance[J].Global Change Biology,2010,16:2420-2435.

[47] Falge E,Baldocchi D,Olson R,et al.Gap filling strategies for defensible annual sums of net ecosystem exchange [J].Agricultural and Forest Meteorology,2001,107 (8):43-69.

[48] Aubinet M,Joly L,Loustau D,et al.Dimensioning IRGA gas sampling system:laboratory and field experiments[J]. Atmospheric Measurement Techniques,2014,9(3):1361-1367.

[49] Leuning R,King K M.Comparison of eddy-covariance measurements of CO2,fluxes by open- and closed-path CO2analysers[J].Boundary-Layer Meteorology,1992,59(3):297-311.

[50] Yasuda Y,Watanabe T.Comparative measurements of CO2flux over a forest using closed-path and open-path CO2analysers[J].Boundary-Layer Meteorology,2001,100:191-208.

[51] Peltola O,Mammarella I,Haapanala S,et al.Field intercomparison of four methane gas analyzers suitable for eddy covariance flux measurements[J].Biogeosciences, 2013,10(6):3749-3765.

[52] Leuning R A Y,Judd M J.The relative merits of open-and closed-path analyses for measurement of eddy fluxes [J].Global Change Biology,2010,2(3):241-253.

[53] Lee X,Black T A,Novak M D.Comparison of flux measurements with open- and closed-path gas analyzers above an agricultural field and a forest floor[J]. Boundary- Layer Meteorology,1994,67(1-2):195-202.

Comparison of CH4Flux Measurement by Open- and Close- Path Eddy Covariance System

YUAN Wen-wen1，ZHANG Jin-song1，MENG Ping1，TONG Xiao-juan2，PAN Qing-mei1，HE Fang-jie1，LI Jian-xia3

(1. Research Institute of Forestry Chinese Academy of Forestry, Beijing 100091, China; 2. Beijing Forestry University, Beijing 100083; 3. Nanshan National Forest Farm, Henan 454650)

s: The eddy covariance technique provides a useful tool to directly measure CH4exchange between the vegetation and the atmosphere. The open-path (OPEC) and close-path (CPEC) systems have uncertainties in measuring methane flux in forests. Therefore, it is necessary to compare the measured CH4flux obtained from the OPEC and CPEC systems. In this study, CH4flux was measured using the OPEC and CPEC systems in a mixed plantation in the Xiaolangdi area during the growing season of 2016 (July 24th to August 5th). The (co)spectra were analyzed and the time lag was estimated. Moreover, CH4flux from the OPEC and CPEC systems was discussed. The results showed that the (co)spectra showed a good agreement with the ?5/3 for spectra and ?4/3 for co-spectra. Under different sampling flow rate (40, 37.5, 35.5, 33.5 and 31.5L·min?1) conditions, the time lag was 4.6, 7.7, 5.3, 10.8 and14.3s, respectively, and the average lag time was 8?9 s. CH4flux derived from the CPEC system was 12% lower than that measured by the OPEC system in sunny days, while CH4flux obtained by the CPEC system was 32% higher than that measured by the OPEC in rainy days. The OPEC system was more suitable measuring CH4flux in comparison with the CPEC system in sunny days. The lag time should be corrected for the CPEC system in rainy days since it can be used for compensate for missing values in the OPEC system. Higher quality data of CH4flux can be obtained by using the CPEC and OPEC systems together.

Eddy covariance；Lag time；Spectrum analysis；Methane flux；Forest ecosystem

10.3969/j.issn.1000-6362.2019.11.001

原文文,張勁松,孟平,等.開路與閉路渦度相關(guān)系統(tǒng)對森林生態(tài)系統(tǒng)甲烷通量觀測的比較[J].中國農(nóng)業(yè)氣象,2019,40(11):669-677

2019?03?06

。E-mail：zhangjs@caf.ac.cn

中央級公益性科研院所基本科研業(yè)務(wù)費(fèi)專項資金（CAFYBB2017ZX002-1）

原文文（1989?），女，研究方向?yàn)榱謽I(yè)氣象。E-mail：yw7198205@163.com