陳云飛，江 洪，2,*，周國模，楊 爽，陳 健

(1.浙江農(nóng)林大學浙江省森林生態(tài)系統(tǒng)碳循環(huán)與固碳減排重點實驗室，浙江 311300;2.南京大學國際地球系統(tǒng)科學研究所，南京 210093)

溫室氣體排放使得地球生物化學循環(huán)碳收支失衡，森林生態(tài)系統(tǒng)在吸收CO2維持生物圈碳氧平衡過程中作用巨大，可計量的森林碳匯日益成為國際研究的熱點。研究生物圈的碳循環(huán)機理和全球陸地生態(tài)系統(tǒng)碳收支及其對環(huán)境變化的響應(yīng)已成為當今一系列大型國際計劃共同關(guān)注的核心問題，科學界急切的要揭示“碳失匯”問題[1]，這部分未知的碳匯有待研究。CO2和H2O通過植被的光合作用和呼吸作用在土壤-植被-大氣圈空間層次上時刻進行著交換，以渦度相關(guān)技術(shù)為主體對土壤-植被-大氣間的CO2/H2O和能量通量以及生態(tài)系統(tǒng)碳水循環(huán)的關(guān)鍵過程進行長期連續(xù)的觀測，所獲取的觀測數(shù)據(jù)將被用來量化和對比分析研究區(qū)域內(nèi)的生態(tài)系統(tǒng)碳收支與平衡特征及其對環(huán)境變化的響應(yīng)[2]。目前渦度相關(guān)技術(shù)已在全球范圍內(nèi)廣泛的應(yīng)用于陸地生態(tài)系統(tǒng)的碳通量吸收與排放測定中，并取得了很好的成效，該方法已成為國際通量觀測網(wǎng)(FLUXNET)的標準方法。中國陸地生態(tài)通量觀測網(wǎng)絡(luò)(ChinaFLUX)也已經(jīng)利用該技術(shù)開展了廣泛的觀測[3-4]。

亞熱帶季風氣候分布在北緯25°—35°亞熱帶大陸東岸地區(qū)，分布地區(qū)包括日本、韓國，美國、巴西的部分地區(qū)，中國的秦嶺-淮河以南的區(qū)域，并以中國面積最大。在亞熱帶通量觀測的系統(tǒng)還很有限，只有江西千煙州人工針葉林觀測站站以及哀牢山地常綠闊葉林觀測系統(tǒng)，兩站點南北跨度大僅對兩大植被類型進行了觀測，而亞熱帶植被類型以常綠闊葉林為主，混雜有落葉林、針葉林、混交林和竹類，亞熱帶季風區(qū)域植被類型豐富，同時多丘陵，森林覆蓋率高，有巨大的碳匯潛力和重要的碳源匯角色。竹林的通量觀測正是對亞熱帶通量觀測的補充和深入研究。

雷竹(Phyllostachysviolascens)作為經(jīng)濟竹類在長江以南均有種植，且種植面積日益擴大。竹林采伐期短，更新代謝快，高效經(jīng)營，收獲竹筍經(jīng)濟效益顯著，同時應(yīng)對氣候變化，增加森林碳匯有著積極的意義。獨特的管理方式和生理生態(tài)過程使其兼具森林、農(nóng)田兩大系統(tǒng)的特點，雷竹林生態(tài)系統(tǒng)較一般森林，作為常綠樹種全年光合能力強，采伐期短更新代謝快，經(jīng)濟林型種植面積廣大。很多研究表明，青幼年森林的碳吸收能力要強于衰老森林。北美黑云杉的采伐對碳源匯的影響[5]，EOBS站點為采伐前林分年齡105a微弱碳匯(6±4)Cg·m-2·a-1， HBS75站點青幼年林型林分年齡33a為強的碳匯(143±35)Cg·m-2·a-1。因此，從生物學角度雷竹始終為青幼年的林齡，應(yīng)具有較強碳匯潛力，但竹林的通量過程與固碳能力國內(nèi)外文獻報道很少，有待于深入研究。在亞熱帶地區(qū)建立長期連續(xù)觀測的通量塔，對揭示竹林碳源匯過程以及影響碳平衡機制有深刻意義，對區(qū)域尺度下特色生態(tài)系統(tǒng)的碳平衡研究有啟發(fā)指導(dǎo)作用。

本文以浙江省臨安太湖源雷竹林生態(tài)系統(tǒng)觀測塔的一年通量數(shù)據(jù)為依據(jù)，結(jié)合雷竹高效經(jīng)營覆蓋耕作科學的管理措施[6]。分析雷竹林系統(tǒng)CO2通量變化過程和碳收支量，及其人為干擾和環(huán)境因子對CO2通量變化的影響。

1 材料與方法

1.1 觀測站點

試驗區(qū)位于浙江省臨安市太湖源鎮(zhèn)，地理位置N 30°18′169″ E119°34′104″屬亞熱帶季風氣候，溫暖濕潤，雨熱同期，氣候特征春季多雨，夏濕熱有梅雨期，秋氣爽，冬干冷。全年降水量1600mm，年平均氣溫16℃，全年日照時數(shù)1900h。試驗地土壤以紅壤為主，海拔185m，坡向北偏東35°，坡度2—3°，雷竹林建園歷史10 a。觀測塔下雷竹林群落平均高度4.5m，平均胸徑4cm，以2、3年生竹為主，總蓋度80％，立竹密度每畝1500株，林下灌木草本很少，有竹葉糠皮覆蓋，為人工經(jīng)營的筍用經(jīng)濟林。雷竹是禾本科竹亞科剛竹屬竹種，由于早春打雷即出筍而得名，實為早竹之變種。自然狀態(tài)下雷竹，稈高7—11m，徑達4—6cm。節(jié)間較短而均勻，幼竿密被白粉生長季以1a為1個周期，連年出筍能力強，在出筍期與筍芽分化期要求有充足降水[7]。

人工高效經(jīng)營的雷竹主要措施是充足肥料供給、充足水分灌溉、鉤梢剪枝、間伐老竹、覆蓋增溫等。其中冬季采用雙層法覆蓋增溫是關(guān)鍵技術(shù)，用礱糠、稻草雙層覆蓋，其厚度30cm(下面鋪稻草10—15cm，上面覆蓋礱糠15—20cm)。自然狀態(tài)下，2月下旬到3月上旬土壤溫度10℃以上，竹出筍最低溫度需9—10℃，最適溫度20—24℃[8]，有效積溫200℃。4月份覆蓋的稻草腐爛入土，礱糠及時移除。覆蓋增溫后，11月底竹林覆蓋稻草礱糠保持地溫10℃以上，使雷竹提前發(fā)筍12月底開始萌動。1月中下旬大量出土，恰趕春節(jié)前后上市，價格較高。雷竹林每年都要收獲竹筍留養(yǎng)新竹，伐走等量的4a老竹，相當于間伐的管理措施，間伐量約占竹林的20％—30％，始終保持一定的立竹密度和1—3a的竹齡，這樣竹林生態(tài)系統(tǒng)處于發(fā)育的青幼時期，其碳匯能力較森林和農(nóng)田系統(tǒng)有明顯優(yōu)勢。

1.2 觀測儀器

觀測竹林建有高20m的微氣象觀測塔，開路渦度相關(guān)系統(tǒng)的探頭安裝在17 m高度上，由三維超聲風溫儀(CAST3， Campbell Inc.， USA) 和開路CO2PH2O 分析儀(Li-7500， LiCor Inc.， USA) 組成，原始采樣頻率為10Hz，數(shù)據(jù)傳輸給數(shù)據(jù)采集器(CR1000， Campbell Inc.， USA) 進行存儲，同時根據(jù)渦度相關(guān)原理在線計算并存儲30min 的CO2通量(FC)、潛熱通量(LE)和顯熱通量(HS)等統(tǒng)計量。常規(guī)氣象觀測系統(tǒng)，包括3層風速(010C， met one， USA) 、3層大氣溫度和濕度(HMP45C，Vaisala， Helsinki， Finland)，安裝高度分別為1m，5m， 17m，2個SI-111紅外溫度分別置于1.5m和5m，用于采集地表和冠層溫度，凈輻射儀(CNR4，Kipp&Zonen)傳感器安裝高度17m，用于采集上行/下行的長波/短波輻射、凈輻射的數(shù)據(jù)，此外還有土壤熱通量(HFP01，Hukseflux)觀測深度3cm，5cm;土壤含水量(CS616，Campbell,USA)觀測深度5,50,100cm；土壤溫度(109，Campbell,USA)觀測深度5,50,100cm。常規(guī)氣象觀測系統(tǒng)數(shù)據(jù)采樣頻率為0.5 Hz，通過數(shù)據(jù)采集器(CR1000，Campbell Inc,USA) 每30 min自動記錄平均風速、溫度、氣壓、凈輻射等常規(guī)氣象信息。常規(guī)氣象數(shù)據(jù)用于插補通量觀測數(shù)據(jù)的缺失。7層CO2/H2O廓線觀測系統(tǒng)(LI-840，Li CorInc,USA)觀測CO2/H2O從土壤到植被再到大氣的擴散過程，為閉路系統(tǒng)，采樣頻率為0.2Hz，同時輸出30min平均值。全部觀測數(shù)據(jù)保存到數(shù)據(jù)采集器的TF卡上，同時配置無線傳輸模塊利用遠程計算機加載虛擬串口通過手機信號連接數(shù)據(jù)采集器，可以實時監(jiān)測傳感器狀態(tài)，30min數(shù)據(jù)可以實時傳輸。

1.3 研究方法

開路渦度系統(tǒng)數(shù)據(jù)采集器可在線計算虛溫訂正和空氣密度變化訂正。但是未考慮地形和儀器傾斜影響，當下墊面有傾斜度時，由于地球引力作用，順著山坡走向大氣會發(fā)生匯流、漏流現(xiàn)象，此時平均垂直風速并不為零。很多研究都是通過旋轉(zhuǎn)風向坐標軸來計算通量[9-11]。在中尺度大氣環(huán)流時，可以旋轉(zhuǎn)坐標軸迫使平均垂直風速為零[12]。根據(jù)風向、儀器底座、主風向地形坡度等建立一新的坐標軸參考系統(tǒng)，經(jīng)過坐標旋轉(zhuǎn)后，凈生態(tài)系統(tǒng)交換量即為渦度相關(guān)通量、儲存項、水平和垂直方向的平流3項的和。常見的有二次坐標旋轉(zhuǎn)、三次坐標旋轉(zhuǎn)和平面擬合，本文采用二次坐標旋轉(zhuǎn)(DR)，使坐標系x軸與平均水平風方向平行，從而使平均側(cè)向風速度和平均垂直風速度為零。

本文的觀測方法采用的是渦度相關(guān)的微氣象觀測方法，通過測定大氣中湍流運動產(chǎn)生的風速脈動和氣體濃度脈動，計算二者協(xié)方差求解通量值[13]。

凈生態(tài)系統(tǒng)碳交換量(NEE)主要是指生態(tài)系統(tǒng)中植物光合作用、植被冠層空氣的碳儲存和生態(tài)系統(tǒng)呼吸消耗的碳排放引起的生態(tài)系統(tǒng)碳儲量變化[14]

微氣象學研究以大氣為對象，當CO2從大氣進入到生態(tài)系統(tǒng)時，定義NEE符號為負；當CO2從生態(tài)系統(tǒng)排放到大氣中時，定義NEE符號為正。生態(tài)學家研究的重點在生態(tài)系統(tǒng)的變化，對于通量的符號規(guī)定為:氣體由大氣圈進入生態(tài)系統(tǒng)的通量值符號為負，反之為正[15]。

凈生態(tài)系統(tǒng)生產(chǎn)力(NEP)，定義的符號剛好和NEE符號相反，生態(tài)系統(tǒng)總交換量(GEE)與生態(tài)系統(tǒng)總初級生產(chǎn)力(GPP)符號也是相反的，這樣陸地和大氣之間的氣體交換過程中的關(guān)系，可用下列方程描述：

GEE=NEE-RE

(2)

NEE=Fc+Fs

(3)

RE=REnight+REday

(4)

式中,Fc為大氣和生態(tài)系統(tǒng)冠層的碳通量，即渦度探頭觀測值，F(xiàn)s為冠層內(nèi)的碳儲存通量(雷竹林高度不足8m年際變化可忽略)；式中RE生態(tài)系統(tǒng)呼吸包括植物自養(yǎng)呼吸以及土壤微生物分解土壤有機質(zhì)和凋落物呼吸通量，分為白天與夜晚計算。夜間生態(tài)系統(tǒng)完全為CO2排放狀態(tài)，NEE數(shù)值上就等于生態(tài)系統(tǒng)呼吸值RE。白天的生態(tài)系統(tǒng)呼吸是通過擬合半小時時長的5cm土壤溫度數(shù)據(jù)與夜間NEE數(shù)據(jù)，建立溫度與生態(tài)系統(tǒng)呼吸的關(guān)系式，推算REday。根據(jù)公式(2)，(4)。分解NEE可以分解成GEE和RE，進而計算出GEE。

1.4 數(shù)據(jù)校正與插補

基于渦度原理開路系統(tǒng)采集的10Hz原始數(shù)據(jù)，因為用紅外氣體儀觀測的CO2氣體濃度是相當于干空氣的質(zhì)量混合比，大氣的溫度、壓力、濕度發(fā)生變化會引起CO2質(zhì)量濃度的變化，需要根據(jù)理想氣體狀態(tài)方程校正氣體密度為摩爾質(zhì)量比。需要進行水汽校正，即WPL校正。根據(jù)垂直平均風速為零假設(shè)，做坐標軸旋轉(zhuǎn)校正。以上校正應(yīng)用EdiResoftwar設(shè)定參數(shù)，加載模塊，計算成為30min步長數(shù)據(jù)。EdiRe software是由愛丁堡大學編寫，專門處理渦度觀測數(shù)據(jù)，并做相關(guān)校正。

在實際觀測中由于受到降水、凝水、昆蟲以及隨機電信號異常等的影響，需要對通量數(shù)據(jù)進行質(zhì)量控制，結(jié)合渦度相關(guān)法通量觀測原理和ChinaFLUX推薦篩選標準，目前通用的渦度數(shù)據(jù)處理，有平均日變化法(MDV)，查表法(LookUp-table)，非線性回歸法(NLR)。本文中應(yīng)用非線性的經(jīng)驗方程，光合有效輻射數(shù)據(jù)偏差較大，插補數(shù)據(jù)的結(jié)果不理想相關(guān)性很小，因此采用平均日變化法。平均日變化法是對缺失數(shù)據(jù)用相鄰幾天同時刻數(shù)據(jù)的平均值進行查補。該方法首先確定平均時段的長度，研究發(fā)現(xiàn)白天取14d、夜間取7d的平均時間長度時偏差最小[16]。

渦度相關(guān)開路系統(tǒng)數(shù)據(jù)易受降水、電信號、儀器故障影響造成數(shù)據(jù)缺失或者不合理，一般來說文獻報道的國際通量網(wǎng)絡(luò)各個站點數(shù)據(jù)缺失和不合理比例在17％—50％[17]。美洲通量網(wǎng)(AmeriFlux)的站點白天缺失和不合理數(shù)據(jù)在20％—35％[18]。本站點白天數(shù)據(jù)有效率在80％左右；夜間湍流不充分CO2過多的沉積在林冠下部不能被探頭檢測到，篩選風摩擦系數(shù)較大的數(shù)據(jù)，有效數(shù)據(jù)在48％左右。表1說明除8月份數(shù)據(jù)有效率偏低外，其余各月數(shù)據(jù)有效率要高于通量觀測網(wǎng)絡(luò)系統(tǒng)的平均水平。

表1 各月有效數(shù)據(jù)統(tǒng)計

2 結(jié)果與分析

2.1 常規(guī)氣象因子的季節(jié)變化

觀測數(shù)據(jù)從2010年10月到2011年9月，構(gòu)成全年數(shù)據(jù)。日平均最高溫2011年7月8日30.21℃，最低溫2011年1月17日-5.50℃。月平均氣象因子見表2。全年最高氣溫7月，最低氣溫1月；光合有效輻射變化與氣溫相同；降雨量最大6月，最小12月，降雨與該地區(qū)往期持平，降雨量變化較氣溫變化年際差異大。

表2 各月氣象因子變化

2.2 凈生態(tài)系統(tǒng)碳交換量

從圖1月每時刻平均的日進程可以看出，NEE呈“U”形變化，表現(xiàn)為在白天通量負值，夜間通量正值，雷竹林生態(tài)系統(tǒng)月尺度上的NEE符號變化時刻和最大的碳排放、碳吸收時刻都有明顯的差異。各月通量符號變化，由正值轉(zhuǎn)為負值在日出1—1.5h，春冬季與夏秋季相比延遲1h；由負值轉(zhuǎn)為正值時刻較集中，發(fā)生在17:00—18:00。9—11月NEE由正轉(zhuǎn)負發(fā)生在6:30—7:00，NEE最大負值出現(xiàn)在11:00前后。在12月—翌年2月NEE由正轉(zhuǎn)負發(fā)生在7:30—8:00，NEE最大負值出現(xiàn)在12:00前后。3—8月的通量符號變化時刻和最大負值時刻介于秋季與冬季之間，NEE由正轉(zhuǎn)負發(fā)生在7:00—7:30，NEE最大負值出現(xiàn)在11:00—11:30。

圖1 CO2通量各月平均日變化

從全年NEE 變化看，年碳吸收總量126.303gC·m-2·a-1，總體是碳匯。這要小于周國模應(yīng)用生物量調(diào)查法對毛竹固碳能力509gC·m-2·a-1，杉木為348gC·m-2·a-1的估算[19]。李靜潔等研究江西紅壤丘陵地區(qū)雙季稻碳收支[20]，在不同施肥措施下估算早稻田兩季稻的碳匯能力為311—629gC·m-2·a-1，雷竹年固碳能力要小于水稻，同時也小于北方農(nóng)田系統(tǒng)198—318gC·m-2·a-1 [21]年吸收能力。具體季節(jié)分析12—2月即冬季覆蓋月份NEE為正，雷竹林表現(xiàn)為碳源，春季、夏季、秋季均為碳匯。全年碳通量，吸收量表現(xiàn)為雙峰變化曲線，第1峰出現(xiàn)在6月(-21.39gCm-2·月-1)；第二峰出現(xiàn)在11月(-35.89gC·m-2·月-1)。最高排放量1月(16.061gC·m-2·月-1)。

6月全月凈積累總量要小于11月，這可能是由于亞熱帶6月份開始進入梅雨期，影響雷竹的光合生產(chǎn)能力。同緯度氣候下千煙洲站點6月出現(xiàn)的凈積累量也小于11月的凈積累量[22]。7—8月受高溫影響有明顯的“午休”現(xiàn)象，CO2通量一般在10:0—11:00達到峰值,之后降低，09—11月適宜的水熱條件使碳吸收持續(xù)增加可為冬季竹筍萌發(fā)積蓄有機物。7—8月份雖然有較強的初級生產(chǎn)力，但由于呼吸過程排放的碳也比較大，碳凈吸收低于5—6月和9—11月。這在很多類型的生態(tài)系統(tǒng)普遍存在，在竹類系統(tǒng)也得到了驗證。12月進入冬季開始覆蓋增溫，人為經(jīng)營措施干擾NEE開始減小。

從圖2可以看出RE呈單峰變化，在夏季較高冬季較低。圖5中顯示夜間NEE與土壤5cm溫度的散點圖，擬合曲線顯示夜間NEE與土壤5cm溫度成指數(shù)關(guān)系變化，相關(guān)系數(shù)0.3。7月土壤溫度最高月份，RE達到峰值151.05 gC·m-2·月-1，在氣溫較低的春冬季呼吸較低。從季節(jié)上分析，分為冬季(12—2月)，春季(3—5月)，夏季(6—8月)，秋季(9—11月)。各季節(jié)RE占全年總量的比例分別為16.7％，21.2％，37.0％，25.1％；各季節(jié)NEE占全年比例分別為-21.3％，26.6％，25.8％，68.9％；各季節(jié)GEE占全年比例為12.8％，21.8％，35.8％，29.6％。全年RE貢獻最大為夏季，夏季生態(tài)系統(tǒng)呼吸是冬季的2.3倍；全年NEE貢獻最大為秋季超過全年其余月份之和，而冬季為負貢獻降低了全年碳匯總量，春夏季NEE貢獻為年平均水平；全年GEE最大貢獻為秋季，最小為夏季，變化與RE一致。

全年GEE季節(jié)變化特征明顯，峰值出現(xiàn)在6—7月，可能與溫度、水分影響有關(guān)；之后開始下降，但秋季生長明顯大于春季，冬季GPP最小。根據(jù)雷竹的生長特性，竹筍出土吐葉到成竹經(jīng)30—40 d，5月大量新竹開始光合生產(chǎn)，6月下旬到7月梅雨期，為雷竹生長提供充足的水分供給。冬季的防雪鉤梢和溫度影響，可能是生產(chǎn)力開始下降的主要影響因素。覆蓋月份為12—3月，出筍是1—3月，竹筍在適宜的水熱條件下爆發(fā)性生長，竹筍出土呼吸加大，同時覆蓋后土壤和覆蓋物分解呼吸大大超出覆蓋前的土壤呼吸，有研究發(fā)現(xiàn)覆蓋后是覆蓋前的5—10倍[23]，造成冬季CO2排放不降。土壤呼吸、覆蓋物呼吸以及植物呼吸都土壤溫度密切相關(guān)，1月溫度降低，生態(tài)系統(tǒng)呼吸隨之降低。6月GEE與RE同時上升出現(xiàn)峰值，新竹不斷增強的光合作用吸收了一定CO2，同時生態(tài)系統(tǒng)呼吸達到了較大值，造成NEE不為最高值。9月GEE與RE同時開始下降時，NEE開始增加，并在11月達到全年的最高峰值。在春季NEE隨GEE與RE的升高而增加，在秋季NEE隨GEE與RE的下降而增加，出現(xiàn)兩個峰值。由此可見人工高效經(jīng)營下，溫度是主要的影響因子，光照次之。雷竹喜濕怕澇，在受水分脅迫時人工會及時補水排澇，水分影響這里先不考慮。

圖2 全年各月累積凈生態(tài)系統(tǒng)CO2交換量(NEE)、生態(tài)系統(tǒng)呼吸(RE)和總生態(tài)系統(tǒng)CO2交換量(GEE)變化過程

圖3 固碳效率變化

2.3 竹林固碳效率

在一定的時間尺度上，生態(tài)系統(tǒng)碳交換通量與碳固定和排放速率根據(jù)選定的研究界面和對象能夠相互直接轉(zhuǎn)化[24]，在數(shù)值上GEE可等同為生態(tài)系統(tǒng)固碳速率，RE即為生態(tài)系統(tǒng)碳排放速率，NEE即生態(tài)系統(tǒng)凈固碳速率，量綱不變?nèi)詾間C·m-2·a-1。相對于植物生理研究，定義有碳利用效率(CUE)，就是植物凈積累的碳占總光合生產(chǎn)中固定碳的比例[25]，即有CUE=NPP/GEP[26]，該指標表示植物碳同化能力的強弱，但是對于整個生態(tài)系統(tǒng)的碳收支研究缺失土壤呼吸和植物凋落物等異養(yǎng)呼吸碳排放計算，不能代表整個生態(tài)系統(tǒng)碳收支水平。如果將大氣-植被-土壤看作一個連續(xù)的整體，研究在生態(tài)系統(tǒng)氣體交換過程中生態(tài)系統(tǒng)的碳收支，將生態(tài)系統(tǒng)凈固碳速率與生態(tài)系統(tǒng)固碳速率比值，可以暫時定義為生態(tài)系統(tǒng)碳固定效率(ECSE)，反映生態(tài)系統(tǒng)氣體交換中固定下來碳的比例，即ECSE=NEE/GEE。從圖3中看出，生態(tài)系統(tǒng)的碳利用效率從1月份開始為負值表明，是極強的碳釋放，冬季雷竹雖然有一定的光合能力但生態(tài)系統(tǒng)呼吸大大超過光合生產(chǎn)的凈積累量；3月開始，變?yōu)檎缔D(zhuǎn)，為固碳并保持在10—16％；夏季7—8月固碳能力下降，碳排放的比例增加，進入秋季達到峰值，有生態(tài)系統(tǒng)氣體交換過程有33％被固定下來，全年固碳效率為11％。

2.4 高效經(jīng)營覆蓋措施對CO2通量的影響

圖4 土壤5cm溫度與空氣1m溫度的比較

冬季覆蓋技術(shù)是雷竹高效經(jīng)營的關(guān)鍵技術(shù)，覆蓋的主要作用是增加土壤溫度。夏季土壤表層溫度高于深層溫度，在冬季則相反，土壤溫度隨深度增加而升高。在冬季增加覆蓋措施后5cm土壤溫度要高于1m處氣溫5—10℃。從圖4中看出在12—2月，地表5cm土壤溫度明顯高出空氣溫度，平均溫度都維持在10℃以上，這就加大了生態(tài)系統(tǒng)的呼吸，同時冬季光合有效輻射為全年最低，形成冬季月份碳源。

圖5 土壤5cm溫度變化與夜間NEE的關(guān)系

全年土壤溫度變化范圍9—26℃，波動范圍不大，夜間NEE可以代表生態(tài)系統(tǒng)呼吸即包括植物暗呼吸、土壤微生物呼吸以及凋落物分解的CO2排放，最高呼吸在0.6—0.7mg·m-2·s-1。全年的土壤5cm溫度與經(jīng)摩擦風速篩選后夜間NEE關(guān)系并不明顯，因為點數(shù)太多，將各個溫度下的夜間NEE做平均處理，得到圖5。從圖5中，可以看出生態(tài)系統(tǒng)與土壤5cm成指數(shù)關(guān)系，但相關(guān)系數(shù)只有0.3，成指數(shù)關(guān)系這與之前有關(guān)文獻報道的研究結(jié)果較一致[16,27]。ChinaFLUX站點也有相似的報道，陳述悅研究華北麥田發(fā)現(xiàn)5cm地溫同土壤呼吸相關(guān)性最好[28]，劉允芬等研究千煙洲人工針葉林土壤溫度與夜間NEE數(shù)值呈指數(shù)關(guān)系相關(guān)系數(shù)0.38[29]。相關(guān)系數(shù)不大的原因可能有兩個，首先冬季的覆蓋增溫措施干擾了土壤自然的呼吸，使溫度變化范圍減小，生物菌肥的使用加大了呼吸；再有溫度升高使植物自養(yǎng)呼吸增加，生態(tài)系統(tǒng)呼吸也隨之增加。

假定雷竹光合能力不變，數(shù)據(jù)顯示冬季無覆蓋時近地1m空氣溫度與土壤5cm溫度相差不大，用1m處空氣溫度代替無覆蓋的5cm土壤溫度，根據(jù)與夜間NEE的指數(shù)方程模擬無覆蓋的生態(tài)系統(tǒng)呼吸。以溫度為影響因子比較無覆蓋與有覆蓋的碳排放，模擬結(jié)果說明12—2月單位面積上，覆蓋后分別多增加21.56g，26.52g，20.53g的排放(表3)。因此若雷竹林無覆蓋可少排放CO2，相應(yīng)冬季月份可由碳源轉(zhuǎn)匯。

冬季覆蓋增加的碳排放部分來自覆蓋的有機物分解。從表4中可以計算不包括地下部分生長計算碳收支平衡，其中雷竹每株干重計以3kg計算，每年每畝伐300株，地上部分碳密度參照周國模研究毛竹碳密度[19]以500g/kg計算，竹筍收獲以1500kg，含水90％計算，總年收獲量為525kg碳，即787.50Cg·m-2·a-1。竹筍收獲占14％，竹材收獲占86％，竹材的碳可以作為碳庫多年儲存，竹筍的碳運出系統(tǒng)后又將排放到大氣中，則雷竹林系統(tǒng)為十分微弱的碳源。雷竹林NEE較水稻田、北方農(nóng)田都小，但收獲量缺明顯高出其他農(nóng)田系統(tǒng)，可見高效經(jīng)營措施對收獲量產(chǎn)出的影響。

覆蓋物投入的有機碳含量以及礱糠收回，稻草則腐爛入土的比例需要進一步實驗，對竹林系統(tǒng)碳庫影響都還不甚清楚。目前有關(guān)覆蓋物對雷竹林碳儲量影響、雷竹土壤有機質(zhì)轉(zhuǎn)化研究的文獻還很少見，初步估計，每年每公頃投入覆蓋的稻草及礱糠等有機物料共100 t左右，大量的有機質(zhì)輸入對土壤有機碳庫的組成分解及轉(zhuǎn)化應(yīng)產(chǎn)生重要影響[30-32]。有文獻研究15a種植雷竹有機碳轉(zhuǎn)化變化，覆蓋經(jīng)營后極大地提高了土壤有機質(zhì)含量，每年每畝表層有機質(zhì)增加達到616 kg，占投入物料的23.8％。該值略低于南方紅壤區(qū)作物秸稈的平均腐殖化系數(shù)(0.22—0.34，平均 0.28)，可能是由于礱糠 C/N 比較高，不易腐殖化[33]。每畝1.6t有機碳投入，以23.8％轉(zhuǎn)化系數(shù)計算，將有380kg有機碳進入系統(tǒng)，即增加了碳排放571Cg·m-2·a-1，這些有機碳在系統(tǒng)中如何分配的有待深入研究。

表3 模擬覆蓋與無覆蓋碳排放變化

表4 高效經(jīng)營中投入產(chǎn)出(以666.7m2計)

3 結(jié)論與展望

本文分析了太湖源人工經(jīng)營的雷竹林碳通量各月的平均日變化趨勢，季節(jié)碳通量的變化趨勢，同時探究了冬季覆蓋增溫過程中對碳通量變化以及碳平衡的影響。研究表明：

人工經(jīng)營的雷竹林生態(tài)系統(tǒng)全年各月碳通量平均日變化趨勢明顯；2010年10月到2011年9月碳收支為-126.303gC·m-2·a-1，RE為1108.845 gC·m-2·a-1，GEE為-1235.15gC·m-2·a-1。整體為碳匯，冬季覆蓋的12—2月份為碳源，其余月份為碳匯。通量過程具有明顯的季節(jié)性，雷竹林生態(tài)系統(tǒng)月尺度上的NEE符號變化時刻和最大的碳排放、碳吸收時刻都有明顯的差異，通量符號由正值轉(zhuǎn)為負值在日出1h—1.5h，春冬季與夏秋季相比延遲1h；由負值轉(zhuǎn)為正值時刻較集中，發(fā)生在17:00—18:00；最大碳匯時刻秋季在11:00，冬季在12:00，春夏介于11:00—12:00之間。溫度是影響碳通量的主要因子，5cm處土壤溫度與生態(tài)系統(tǒng)呼吸存在明顯的指數(shù)關(guān)系。在春季NEE隨NEE與RE的同時升高而增加，在秋季NEE隨NEE與RE的同時下降而增加，出現(xiàn)兩個峰值并有碳匯總量秋季大于春季。秋季竹林NEE貢獻達68.9％，同時固碳效率最高為33％。冬季覆蓋通過投入有機物料，提高地溫，同時增大了生態(tài)系統(tǒng)呼吸，增加了碳排放。

竹林具有很強的森林碳匯潛力，合理的耕作措施，科學采伐老竹，保持生態(tài)系統(tǒng)活力，可為陸地生態(tài)系統(tǒng)維持碳氧平衡應(yīng)對氣候變化做出積極的貢獻。另一方面，對于竹林生態(tài)系統(tǒng)碳循環(huán)研究，本文僅對雷竹林的碳源匯做了結(jié)論，對與影響碳吸收、排放的各個分量未做計算。同時本站點僅對雷竹林一年的碳通量過程研究進行初步觀測研究，時間較短，因碳通量過程有明顯的年際變化易受氣候變化影響，要通過長期的連續(xù)觀測和相應(yīng)的輔助實驗，來做深入研究。

[1]Tans P P,Fung I Y,Takahashi T.Observation constraints on the global atmospheric CO2budget.Science,1990,247(4949):1431-1438.

[2]Berbigier P,Bonnefond J M,Mellmann P.CO2and water vapour fluxes for 2 years above Euroflux forest site.Agricultural and Forest Meteorology,2001,108(3):183-197.

[3]Yu G R,SunX M.Flux Measurement and Research of Terrestrial Ecosystem in China.Beijing:Science Press,2008.

[4]Yu G R,Sun X M.Principles of Flux Measurement in Terrestrial Ecosystems.Beijing:Higher Education Press,2006:4-4.

[5]Payeur-Poirier J L,Coursolle C,Hank A,Giasson M A.CO2fluxes of a boreal black spruce chronosequence in eastern North America.Agricultural and Forest Meteorology,2012,153:94-105

[6]FangW,He J C,Lu K X,Chen J H.Cultivation techniques of early shooting and high yielding for Lei Bamboo Sprout.Journalof Zhejiang Forestry College,1994,11(2):121-128.

[7]Wang Z T,FangW,He J C.High Yield Cultivation Techniques of Lei Bamboo shoots.Beijing:China Forestry Press,1993.

[8]Xu L L,Hu Y B,Hu Y C,Sun C M,Yu M C,Liang K,Lang J B.Lei bamboo production development and covered technology.Shanghai Agricultural Science and Technology,2008,(2):85-86.

[9]Baldocchi D D,Hicks B B,Meyers T P.Measuring biosphere-atmosphere exchanges of biologically related gases with micrometeorological methods.Ecology,1988,69(5):1331-1340.

[10]McMillen R T.An eddy correlation technique with extended applicability to non-simple terrain.Boundary-Layer Meteorology,1988,43(3):231-245.

[11]FokenTh,Wichura B.Tools for quality assessment of surface-based flux measurements.Agricultural and Forest Meteorology,1995,78(1/2):83-105.

[12]Paw U K T,Baldocchic D D,Meyers T P,Wilson K B.Correction of eddy-covariance measurements incorporating both advective effects and density fluxes.Boundary-Layer Meteorology,2000,97(3):487-511.

[13]Reynolds O.On the dynamical theory of incompressible viscous fluids and the determination of the criterion.Philosophical Transactions of the Royal Society of LondonA,1895,186:123-164.

[14]Garratt J R.Limitations of the eddy-correlation technique for the determination of turbulent fluxes near the surface.Boundary-Layer Meteorology,1975,8(3/4):255-259.

[15]Song X,Liu Y F,Xu XF,Yu GR,Wen X F.Comparison study on carbon dioxide,water and heat fluxes of the forest ecosystem in red earth hilly zone over winter and spring.Resources Science,2004,26(3):96-104.

[16]FalgeE,Baldocchi D,Olson R,Anthoni P,Aubinet M,Bernhofer C,Burba G,Ceulemans R,Clement R,Dolman H,Granier A,Gross P,Grünwald T,Hollinger D,Jensen N O,Katul G,Keronen P,Kowalski A,Ta Lai C,Law B E,Meyers T,Moncrieff J,Moors E,Munger J W,Pilegaard K,Rannik ü,Rebmann C,Suyker A,Tenhunen J,Tu K,Verma S,Vesala T,Wilson K,Wofs S.Gap filling strategies for defensible annual sums of net ecosystem exchange.Agricultural and Forest Meteorology,2001,107:43-69.

[17]Katul G,Lai C T,Sch?fer K,Vidakovic B,Ellsworth D,Oren R.Multiscale analysis of vegetation surfacefluxes:from seconds to years.Advances in Water Resources,2001,24(9/10):1119-1132.

[18]Wilson K B,Baldocchi D D.Comparing independent estimates of carbon dioxide exchange over 5 years at a deciduous forest in the southeastern United States.Journal of Geophysical Research,2001,106(D24):34167-34178.

[19]Zhou G M.Carbon Storage,Fixation and Distribution in Mao Bamboo(Phyllostachyspubescens) Stands Ecosystem.Zhejiang:Zhejiang University,2006.

[20]Li J J,Pan G X,Li LQ,Zhang X H.Estimation of net carbon balance and benefits of rice-rice cropping farm of a red earth paddy under long term fertilization experiment from Jiangxi,China.Journal of Agro-Environment Science,2009,28(12):2520-2525.

[21]LiJ,YuQ,Sun X M,Tong X J,Ren C Y,Wang J,Liu E M,Zhu Z L,Lin G R.The north China plain farmland ecological system and its environmental regulation mechanism carbon exchange.Science in ChinaSeries D:Earth Sciences,2006,36(Supplement):210-223.

[22]Liu Y F,Song X X,Sun X M,Wen X F,Chen Y R.CO2flux seasonal changes and environmental factors of influence form Qianyanzhou of Coniferous.Science in ChinaSeries D:Earth Sciences,2006,36(Supplement):109-117.

[23]Zhou G M,Jiang P K,Xu Q F.Bamboo Forest Ecosystem Carbon Fixation and Transformation.Beijing:Science Press,2010.

[24]Yu G R,Wang Q F,Liu Y C,Liu Y H.Conceptual framework of carbon sequestration rate and potential increment of carbon sink of regional terrestrial ecosystem and scientific basis for quantitative carbon Authentification.Progress in Geography,2011,30(7):772-774.

[25]Cabral O M R,Gash J H C,Rocha H R,Marsden C,Ligo M A V,Freitas H C,Tatsch J D,Gomes E.Fluxes of CO2above a plantation of Eucalyptus in southeast Brazil.Agricultural and Forest Meteorology,2011,151(1):49-59.

[26]ChoudhuryB J.Modeling radiation-and carbon-use efficiencies of maize,sorghum,and rice.Agricultural and Forest Meteorology,2001,106(4):317-330.

[24]LloydJ,Taylor J A.On the temperature dependence of soil respiration.Functional Ecology,1994,8(3):315-323.

[27]LloydJ,Taylor J A.On the temperature dependence of soil respiration.Functional Ecology,1994,8(3):315-323.

[28]ChenSS,LiJ,Lu P L,Wang Y H,Yu Q.Soil respiration characteristics in winter wheat field in North China Plain.Chinese Journal of Applied Ecology,2004,15(9):1552-1560.

[29]Liu Y F,Yu G R,Wen X F,Wang Y H,Song X,Li J,Sun X M,Yang F T,Chen Y R,Liu Q L.The season of CO2flux variation characteristics form Subtropical plantation ecological system.Science in ChinaSeries D:Earth Sciences,2006,36(Supplement):91-102.

[30]Jiang P K,Zhou GM,Xu Q F.Effect of intensive cultivation on the carbon pool of soil inPhyllostachyspraecoxstands.Scientia Silvae Sinicae,2002,(6):6-11

[31]Jiang P K,Xu Q F.Effects of fertilization on soil active organic carbon underPhyllostachyspraecoxstand.Chinese Journal of Applied Ecology,2005,(2):253-256.

[32]Cai R R,HuangF,Sun D,Qin H,Yang F,Zhuang S Y,Zhou G M,Cao Z H.Temporal and spatial variation of soil organic matters inPhyllostachyspraecoxstands with intensive cultivation management.Journal of Zhejiang Forestry College,2007,(4):450-455.

[33]Li G D,Liu G Q,Zhuang S R,Gui R Y.Changes of organic matter in soils planted Lei Bamboo with different years.Chinese Journal of Soil Science,2010,41(4):847-848.

參考文獻:

[3]于貴瑞,孫曉敏,等.中國陸地生態(tài)系統(tǒng)碳通量觀測技術(shù)及時空變化特征.北京:中國科學出版社,2008.

[4]于貴瑞,孫曉敏.陸地生態(tài)系統(tǒng)通量觀測的原理與方法.北京:高等教育出版社,2006:4.

[6]方偉,何鈞潮,盧可學,陳健華.雷竹早產(chǎn)高效栽培技術(shù).浙江林學院學報,1994,11(2):121-128.

[7]汪祖潭,方偉,何鈞潮.雷竹筍用林高產(chǎn)栽培技術(shù).北京:中國林業(yè)出版社,1993.

[8]徐繆暢,胡永斌,胡余楚,孫春明,俞茂昌,梁侃,朗進寶.雷竹生產(chǎn)發(fā)展與竹園覆蓋技術(shù).上海農(nóng)業(yè)科技,2008,(2):85-86.

[15]宋霞,劉允芬,徐小鋒,于貴瑞,溫學發(fā).紅壤丘陵區(qū)人工林冬春時段碳、水、熱通量的觀測與分析.資源科學,2004,26(3):96-104.

[19]周國模.毛竹林生態(tài)系統(tǒng)中碳儲量、固定及其分配與分布的研究[D].浙江:浙江大學,2006.

[20]李潔靜,潘根興,李戀卿,張旭輝.紅壤丘陵雙季稻稻田農(nóng)田生態(tài)系統(tǒng)不同施肥下碳匯效應(yīng)及收益評估.農(nóng)業(yè)環(huán)境科學學報,2009,28(12):2520-2525.

[21]李俊,于強,孫曉敏,同小娟,任傳友,王靖,劉恩民,朱治林,于貴瑞.華北平原農(nóng)田生態(tài)系統(tǒng)碳交換及其環(huán)境調(diào)控機制.中國科學D輯:地球科學,2006,36(增刊):210-223.

[22]劉允芬,宋霞,孫曉敏,溫學法,陳永潤.千煙洲人工針葉林CO2通量季節(jié)變化及其環(huán)境因子的影響.中國科學D輯:地球科學,2004,34(增刊):109-117.

[23]周國模,姜培坤,徐秋芳.竹林生態(tài)系統(tǒng)中碳的固定與轉(zhuǎn)化.北京:科學出版社,2010.

[24]于貴瑞,王秋鳳,劉迎春,劉穎慧.區(qū)域尺度陸地生態(tài)系統(tǒng)固碳速率和增匯潛力概念框架及其定量認證科學基礎(chǔ).地理科學進展,2011,30(7):772-774.

[28]陳述說,李俊,陸佩玲,王迎紅,于強.華北平原麥田土壤呼吸特征.應(yīng)用生態(tài)學報,2004,15(9):1552-1560.

[29]劉允芬,于貴瑞,溫學法,王迎紅,宋霞,李菊,孫曉敏,楊風亭,陳永瑞,劉琪瓃.千煙洲中亞熱帶人工林生態(tài)系統(tǒng)CO2通量的季節(jié)變異特征.中國科學D輯:地球科學,2006,36(增刊):91-102.

[30]姜培坤,周國模,徐秋芳.雷竹高效栽培措施對土壤碳庫的影響.林業(yè)科學,2002,(6):6-11.

[31]姜培坤,徐秋芳.施肥對雷竹林土壤活性有機碳的影響.應(yīng)用生態(tài)學報,2005,(2):253-256.

[32]蔡榮榮,黃芳,孫達,秦華,楊芳,莊舜堯,周國模,曹志洪.集約經(jīng)營雷竹林土壤有機質(zhì)的時空變化.浙江林學院學報,2007,(4):450-455.

[33]李國棟,劉國群,莊舜堯,桂仁意.不同種植年限下雷竹林土壤的有機質(zhì)轉(zhuǎn)化.土壤通報,2010,41(4):847-848.