亞熱帶毛竹林生態系統能量通量及平衡分析

孫 成, 江 洪,2,*, 陳 健, 劉玉莉, 牛曉棟, 陳曉峰, 方成圓

1 浙江農林大學浙江省森林生態系統碳循環與固碳減排重點實驗室, 杭州 311300 2 南京大學國際地球系統科學研究所, 南京 210093

亞熱帶毛竹林生態系統能量通量及平衡分析

孫 成1, 江 洪1,2,*, 陳 健1, 劉玉莉1, 牛曉棟1, 陳曉峰1, 方成圓1

1 浙江農林大學浙江省森林生態系統碳循環與固碳減排重點實驗室, 杭州 311300 2 南京大學國際地球系統科學研究所, 南京 210093

利用開路渦度相關系統和常規氣象觀測儀器，對亞熱帶(浙江省)毛竹林生態系統2011年的凈輻射、顯熱通量、潛熱通量、土壤熱通量以及氣溫、地溫、降雨量等氣象要素進行了連續觀測，定量分析了毛竹林生態系統能量通量的變化和各能量分量的分配特征，并計算了能量閉合度以及波文比。結果表明：毛竹林全年凈輻射為2628.00 MJ/m2，顯熱通量為576.80 MJ/m2，潛熱通量為1666.77 MJ/m2，土壤熱通量為-7.52 MJ/m2，土壤為熱源，各能量分量季節變化明顯，日變化基本呈單峰型曲線變化。顯熱通量占凈輻射的22.0%，潛熱通量占63.4%，毛竹林生態系統潛熱通量為能量散失的主要形式。波文比逐月變化規律不明顯，波動較大，在0.07—1.77之間變化，能量平衡比率法得出毛竹林年能量閉合度為0.85，月平均閉合度為0.84，能量閉合度高于線性回歸法計算結果，但仍有15%的能量不閉合。

毛竹林; 能量通量; 能量閉合度; 湍流能量; 有效能量; 波文比

森林作為地球上最大的陸地生態系統，在全球水循環和能量再分配中都發揮著重要的作用。而森林生態系統能量通量能夠直接影響地面溫度、水分輸送及植被生長發育與生態系統生產力[1]。目前，由于渦動相關法具有不干擾生態系統、時間分辨率高等優點，已逐漸成為生態系統地-氣間能量和物質通量觀測研究的標準方法，并利用渦動相關技術對各種生態系統的水熱過程進行了大量研究[2- 4]。根據渦度相關技術觀測的基本假設，理論上能量平衡閉合度可以作為觀測系統性能和數據質量評價的一個有效途徑。21世紀初在美國佛羅里達州開展的EBEX-2000實驗就是專門針對地表能量通量平衡問題進行的觀測研究[5- 6]。就目前來看，土壤-植被-大氣之間的能量儲存未充分計入是造成地表能量不平衡的主要原因之一[7- 9]。因此，研究毛竹林的能量過程，分析各能量分量的分配特征，有助于揭示生態系統水分和能量平衡過程，對提高毛竹林生態系統生產力有一定的實踐意義。但目前利用渦度相關技術對毛竹林生態系統能量通量的研究尚未見文獻報道。

目前亞熱帶地區通量觀測系統有江西千煙洲和大崗山觀測站，云南哀牢山觀測站，廣東鼎湖山觀測站和湖南會同等觀測站。幾個站點南北跨度較大，主要針對常綠闊葉林和針葉林進行觀測，亞熱帶植被類型以常綠闊葉林為主，混雜有針葉林、混交林、落葉林和竹類。毛竹(Phyllostachysedulis)稈高達13—20 m，胸徑達12—18 cm，節間短，壁厚，筍期3—5月，連年出筍能力強，是我國南方的重要植被類型和經濟林，也是我國竹林中分布最廣、面積最大的一種竹林，因此開展毛竹林的生態系統能量通量研究具有重要意義。

本文對毛竹林的通量觀測正是對亞熱帶通量觀測很好的補充和創新研究。本文利用浙江省安吉縣山川鄉渦度相關系統觀測塔2011年的氣象和能量通量觀測數據，對毛竹林能量通量進行研究分析。旨在為毛竹林生產和區域模型模擬提供基礎數據，并客觀評價渦度相關觀測數據的質量以及準確評估毛竹林生態系統與大氣間的物質和能量交換提供依據。

1 材料與方法

1.1 研究區概況

研究區位于浙江省湖州市安吉縣山川鄉，地理位置30°28′34.5″N，119°40′25.7″E，屬亞熱帶季風氣候。研究區竹林面積2155 hm2，毛竹林面積1693 hm2，占竹林面積78.6%，分布面積很廣。研究區年平均氣溫16.6 ℃，1月溫度最低，平均氣溫-0.4—5.5 ℃，7月溫度最高，平均氣溫24.4—30.8 ℃，年降水量761—1780 mm，年平均相對濕度均在70%以上。觀測塔站點海拔380 m，毛竹群落平均高度15.8 m，平均胸徑14 cm，以4—6年生竹為主，林下灌木草本很少。

毛竹和一般樹木相比，具有生長快，成材時間短(3—6a)，產量高，繁殖更新能力強等特點，具有易于人工種植，一次栽造成功，即可年年抽筍成竹等優點。研究區毛竹林均采用粗放式經營，只在每年筍期3—5月人工挖筍，11月左右鉤稍，減少對毛竹林生態系統的人為干擾。

1.2 觀測儀器

觀測林地建有高40 m的微氣象觀測塔，開路渦度相關系統的探頭安裝在距地面38 m高度上，由三維超聲風速儀(CSAT3, Campbell Inc., 美國) 和開路CO2/H2O分析儀(Li- 7500, LiCor Inc., 美國) 組成，原始采樣頻率為10 Hz,所有數據利用數據采集器(CR1000, Campbell Inc., 美國) 進行存儲，同時在線計算并存儲30 min的CO2通量(FC)、摩擦風速(Ustar)、潛熱通量(LE)和顯熱通量(H)等參數。

氣象數據的觀測采用40 m高的CO2通量觀測系統中的氣象觀測設施，包括3層風速傳感器(010C，metone，美國)、3層大氣溫度和濕度傳感器(HMP45C，Vaisala，Helsinki，芬蘭)，安裝高度分別為1、7、38 m；一臺凈輻射儀傳感器(CNR4，Kipp&Zonen，荷蘭)安裝于38 m，用于采集下行的長/短波輻射、上行的長/短波輻射、凈輻射的數據；2個SI- 111紅外溫度觀測儀置于2 m和23 m，分別用于采集地表和冠層溫度；土壤溫度傳感器(109，Campbell，美國)，分別距土表5、50、100 cm深度處，水平插入各一個探針。林內分布4個自動雨量筒，林外空處分布1個，降雨量來自林外空地雨量數據。土壤熱通量采用土壤熱通量板(HFP01，Hukseflux，荷蘭)進行測定，選取具有很好代表性的點，挖開土壤剖面，分別距土表3、5 cm深度，水平插入，并將土壤回填。本站點土壤熱通量采用5 cm深度觀測數據，同時，3 cm深度的數據用于插補和校正。全部觀測數據保存到數據采集器的TF卡上，同時配置無線傳輸模塊，利用遠程計算機加載虛擬串口通過GPRS連接數據采集器，實時監測傳感器狀態。

1.3 數據處理

在實際觀測中由于受到降水、凝水、昆蟲以及隨機電信號異常等的影響，需要對通量數據進行有條件的剔除。開路渦動相關系統據此計算的顯熱通量(H，MJ/m2)、潛熱通量(LE，MJ/m2)、CO2通量必須進行剔除處理。常規氣象儀器受天氣影響較小，觀測的太陽凈輻射(Rn，MJ/m2)、土壤熱通量(G，MJ/m2)可不做剔除處理。顯熱和潛熱通量根據三維超聲風速儀觀測的風速、風向與虛溫計算得到，CO2通量數據質量與顯熱、潛熱通量數據具有一致性，因此本文以CO2通量數據質量為標準來判斷相應時刻下的能量通量數據。分別經過二次坐標旋轉校正，水汽校正(WPL)，再根據閾值剔除異常通量值[10]，對于各能量分量缺失值要進行插補，缺失1 d內數據采用線性插補法，缺失超過1 d的數據，利用該能量分量與該月凈輻射回歸關系進行插補。在實際分析時，經過數據質量控制和數據插補后的數據都處理成半小時平均資料，并進一步計算。

1.4 能量平衡研究方法

生態系統觀測中，當能量向下進入生態系統時，Rn取正值；反之，Rn取負值。地表能量平衡方程表達式為[11]：

Rn-G-S-Q=H+LE

(1)

式中,Rn為凈輻射；G為土壤熱通量；S為植被冠層熱儲存量；Q為附加能量項的總和，H為顯熱通量，LE為潛熱通量。由于S和Q項小而常被忽略，此時草原能量平衡方程可表示為：

Rn-G=H+LE

(2)

式中,(Rn-G)簡稱為有效能量，(H+LE)簡稱為湍流能量。當有效能量與湍流能量相等時，稱為能量平衡閉合，否則稱為能量平衡不閉合。

目前，國際上常用的能量閉合評價方法有4種，即最小二乘法(OLS)、壓軸回歸法(RMA)、能量平衡比率法(EBR)和能量平衡殘差頻率分布圖法。本文采用能量平衡比率法(EBR)和線性回歸(OLR)分析方法來分析毛竹林能量平衡狀況。EBR是指渦度相關系統直接測定的湍流能量與有效能量的比值。OLR方法是根據最小二乘法原理求出回歸斜率和截距，分析能量平和閉合程度，在理想狀況下有效能量和湍流能量的回歸直線的斜率為1，并通過原點，但通常二者線性關系的截距不能通過原點。

2 結果與分析

2.1 氣象因子的變化

由圖1可以看出，毛竹林全年降雨量整體呈單峰型曲線變化，年降雨量為1543.1 mm，且雨水主要集中在夏季(1062.2 mm)，占全年降雨量的68.8%。最大降雨量出現在6月份為453 mm，最小出現在2月份為16.3 mm。空氣相對濕度均在70%以上，整體上隨降雨量的增加相對濕度會逐漸升高，其中6月份，受梅雨季節影響，降雨量和相對濕度達均達到全年最大值。

全年月平均氣溫、地溫和土壤5 cm處溫度整體呈明顯的單峰型曲線變化，三者變化趨勢基本保持一致。氣溫、地溫和土壤5 cm處溫度最高值均出現在7月份，最低值均出現在1月份。其中2月份的地表溫度變化較緩慢，是由于冬季冰雪覆蓋，溫度回升較氣溫和土壤溫度慢一些，到3月份變化基本同步。

圖1 毛竹林降雨量(P)、空氣相對濕度(RH)、氣溫(Ta)、地溫(Ts)和土壤5cm處溫度(Ts- 5)的月變化

2.2 能量通量日變化特征

陸地生態系統通過感熱和潛熱的形式與大氣進行熱量和水汽交換。不同類型生態系統的群落類型和下墊面不同，造成蒸發散和熱傳導能力的差異，因此生態系統獲得凈輻射能量后，能量在系統內的分配變化特點各異。為分析各季節毛竹林地表各主要能量分量的日變化特征，將有典型季節性代表的1月份、4月份、7月份、10月份的半小時時刻下的能量通量數據作月平均處理，以表征該月的能量通量日變化進程(圖2)。可見亞熱帶毛竹林能量平衡各分量呈顯著的日變化，季節差異也較明顯，主要分量呈明顯的單峰型變化趨勢。

由圖2顯示，4個月全天的能量分量變化均以凈輻射為基礎，呈單峰型曲線變化。凈輻射日最大值出現時刻差異不大，在11:00—12:30。冬季，7:00—7:30凈輻射通量變為正值，即能量開始進入毛竹林成為收入項，11:00達到日最大值258.49 W/m2，到16:00—16:30轉變成負值，即能量成為支出項。一天中凈輻射大于零的平均通量為158.74 W/m2，全天凈輻射平均通量為39.73 W/m2。春季，在6:00—6:30凈輻射通量變為正值，11:30達到日最大值497.92 W/m2，至17:00—17:30轉變成負值。一天中凈輻射大于零的平均通量為300.78 W/m2，全天凈輻射平均通量為114.60 W/m2。夏季，6:00—6:30凈輻射通量變為正值，日最大值為504.41 W/m2，至18:00—18:30轉變為負值。一天中凈輻射大于零的平均通量為291.48 W/m2，全天凈輻射平均通量為131.27 W/m2。秋季，6:30—7:00凈輻射通量變為正值，日最大值為386.95 W/m2，至16:30—15:00轉變為負值。一天中凈輻射大于零的平均通量為218.72 W/m2，全天凈輻射平均通量為68.96 W/m2。

顯熱通量和潛熱通量均與凈輻射具有類似的日變化特征，但其日變化曲線均不如凈輻射曲線平滑，而潛熱通量波動較大可能與間歇性湍流和夜間湍流的低估有關[12]。從顯熱通量來看，冬季由于地面地溫凍結，加之降雨稀少，顯熱通量值較大，最大值為126.96 W/m2。春季，土壤解凍和冰雪融化，雖降雨量多于冬季，但由于凈輻射顯著增加，導致顯熱通量值也增加，達到214.04 W/m2。夏季是亞熱帶的主要降雨時段，空氣相對濕度較高，能量大部分用來水汽傳輸，顯熱通量較小，為96.35 W/m2。秋季，降水量減少，顯熱通量增加，達到162.86 W/m2。各季節的潛熱通量日變化趨勢與顯熱通量非常相似，但在數值上冬季最弱，最大值為38.27 W/m2，春季最大值為154.25 W/m2，夏季最強，達到297.12 W/m2，秋季最大值為143.63 W/m2。

土壤深層向土壤表層或大氣中釋放熱量時土壤熱通量記作負值，土壤為熱源；相反，土壤表層或大氣向土壤深層傳遞熱量時土壤熱通量記作正值，土壤為熱匯。圖2顯示，土壤熱通量日變化有明顯的差異。冬季變化趨勢不明顯，變化范圍為-7.74—-6.76 W/m2，春季凈輻射增加，加之毛竹出筍期呼吸作用強烈，土壤熱通量變化較大，在-7.23—56.68 W/m2，夏季受江南梅雨季節影響，變化范圍為-1.35—43.27 W/m2，秋季為-11.28—7.84 W/m2。土壤熱通量的峰值出現時間較凈輻射延遲約0.5—1 h，因不同季節土壤理化性質不同，土壤熱導率不同，影響土壤吸熱散熱在延遲時間上的差異[13,14]，在熱源/匯上也有差異[15]。

圖2 不同季節毛竹林能量平衡分量的日變化

2.3 能量通量的季節變化

圖3 毛竹林能量分量的月積累

以年為時間尺度，毛竹林的凈輻射、顯熱通量、潛熱通量和土壤熱通量近似呈單峰型變化。但是受中小尺度天氣變化的影響，能量平衡各分量呈現一定程度的鋸齒狀波動，特別是夏季雨期差異較大。由圖3可以看出毛竹林凈輻射的最大值出現在7月(343.67 MJ/m2)，最低值出現在1月(1.6.31 MJ/m2)，全年凈輻射總量為2628.00 MJ/m2。其中6月份凈輻射有明顯的降低，為222.76 MJ/m2，這是由于6月份為江南梅雨季節，大量陰雨天氣造成的，而6月份的降雨量也達到全年最大值453 mm(圖1)。

顯熱通量的變化與凈輻射基本同步，但在6月份沒有出現降低的現象。最大值出現在5月(96.06 MJ/m2)，最低值出現在12月(5.05 MJ/m2)，全年總量為576.80 MJ/m2，占凈輻射總量的22.0%。潛熱通量為蒸散耗熱，1—7月份隨水熱同期逐漸升高，7月份達到峰值284.59 MJ/m2，之后，潛熱通量逐漸下降，全年總量為1666.77 MJ/m2，占凈輻射總量的63.4%。與黃土高原半干旱草原相比[16]，毛竹林的凈輻射高于半干旱草原(2269.23 MJ/m2)，其能量分配方向為潛熱通量大于顯熱通量。土壤熱通量全年變化不明顯，而且在數值上比其他能量分量小1—2個數量級。土壤熱通量在7月份達到正向最大值(18.71 MJ/m2)，1月份達到負向最大值(19.33 MJ/m2)，全年總量為-7.52 MJ/m2，土壤為熱源，僅占凈輻射總量的0.3%。

2.4 能量通量分配特征及波文比

具體分析能量分量一方面可以了解能量閉合狀況，另一方面可依據凈輻射的主要消耗項，分析和判斷下墊面的干濕狀況。將毛竹林2011年能量通量按雨季(6—9月份)和非雨季(1—5月份，10—12月份)分別分析凈輻射與顯熱通量、潛熱通量的關系。

圖4 毛竹林波文比月變化

毛竹林雨季顯熱通量為203.89 MJ/m2，凈輻射為1126.48 MJ/m2，占凈輻射的18.1%，潛熱通量為825.89 MJ/m2，占凈輻射的73.3%；非雨季顯熱通量為372.90 MJ/m2，凈輻射為1501.52 MJ/m2，占凈輻射的24.8%，潛熱通量為840.88 MJ/m2，占凈輻射的56.0%。毛竹林雨季降水量占全年總量的75.5%，地表相對濕潤，空氣相對濕度接近80%，潛熱通量在能量分配中所占比例較大，潛熱通量約為顯熱通量的4.05倍。非雨季由于降水量較少(377.3 mm)，地表相對干燥，空氣相對濕度為75%，潛熱通量約為顯熱通量的2.25倍。全年顯熱通量占凈輻射的22.0%，潛熱通量占凈輻射的63.4%，總體來看，毛竹林潛熱通量在能量分配中起主導作用，潛熱通量約為顯熱通量的3倍，這與華北農田[17]結果相似。

波文比能夠表征大氣-地表能量交換特征，多用于能量平衡計算[18]。鑒于本文中土壤熱通量占凈輻射的比重僅為0.3%，因此凈輻射能量主要分配給顯熱通量和潛熱通量，而顯熱通量與潛熱通量之比即為波文比。其大小決定表明了能量在生態系統中的分配。由圖4可以看出，全年波文比的變化波動較大，1月最大為1.77，7月最小為0.07；1月份顯熱通量大于潛熱通量，其余月份均小于潛熱通量，月平均波文比0.47，年波文比0.35，全年能量分配潛熱通量大于顯熱通量。黃土高原半干旱草原[16]，年波文比在0.5—3.5，年均值接近1，鼎湖山針闊混交林[19]波文比在0.4—3.0之間，長江灘地抑螺防病林[20]的年波文比在0.5—10之間變化。波文比同時受到日出時間、入射凈輻射量、降雨量以及毛竹的生長過程等的共同影響，且有一定的年際效應。

2.5 能量平衡分析

圖5 毛竹林能量的月積累與能量閉合度變化

根據2011年全年數據，利用能量平衡比率法對渦度相關測得的湍流能量(H+LE)與有效能量(Rn-G)進行閉合度分析(圖5)。EBR有明顯的月變化，最大值出現在6月份為1.19，最小值出現在1月份為0.52，除6月份為能量閉合過度現象以外，其余月份均為能量不閉合，月平均閉合度為0.84。當用湍流通量的年總量和有效能量的比值表示能量閉合狀況時，年EBR值為0.85。這表明能量不閉合度為15%，采用渦度相關法測定的湍流通量僅為常規氣象觀測有效能量的85%。其中，雨季的閉合度達到0.96，非雨季為0.78，雨季的能量平衡狀況較好。

理論上能量閉合是成立的，但不論下墊面性質如何，出現能量不閉合的現象都比較嚴重。此現象也是近20年來困擾地氣相互作用實驗研究的主要難點之一，也是生態學家和微氣象學家所關注的焦點問題。而導致能量不閉合的原因有很多種，地表到土壤熱通量板之間的土壤儲存熱會對地表能量閉合度產生重要影響[9,17]，土壤熱通量的觀測距離地表有一定深度，土壤熱通量的相位會隨土壤深度的加深而延遲，導致能量閉合度降低[7,21]，還有在復雜環境下，湍流不充分和損失[22]等。針對毛竹林的能量平衡狀況，還需要通過長期連續的觀測，進一步比較分析。

在半小時尺度上，表1給出了全年逐月線性回歸系數。截距變化為-5.58—16.03 W/m2，斜率變化為0.47—0.68，相關系數變化為0.68—0.89；年平均截距3.50 W/m2，平均斜率0.55，平均相關系數0.81。這一結果低于黃土高原干旱草地的平均截距(17.22 W/m2)、平均斜率(0.69)和相關系數(0.95)[16]，也低于同一地區雷竹林的平均截距(4.71 W/m2)、斜率(0.59)和相關系數(0.84)[23]。國際通量站點的能量閉合度是：斜率在0.55—0.99，相關系數在0.64—0.96[24]，中國通量網中8個站點[25]的截距、斜率和相關系數平均值分別為28(10—79.9 W/m2)、0.67(0.49—0.81)和0.82(0.52—0.94)。毛竹林回歸系數中斜率處于低水平，相關系數處于平均水平，需要進一步提高剔除數據的質量，來提高能量閉合度。

表1 2011年逐月線性回歸系數Table 1 Ordinary linear regression(OLR)coefficients for energy balance

3 結論與討論

亞熱帶毛竹林全年凈輻射為2628.00 MJ/m2，顯熱通量、潛熱通量和土壤熱通量分別為576.80 MJ/m2、1666.77 MJ/m2和-7.52 MJ/m2，土壤為熱源。在能量分配上，毛竹林主要以潛熱通量形式加熱大氣，通過蒸發散耗能，占凈輻射的63.4%，顯熱通量占22.0%，土壤熱通量僅占0.3%。各能量分量日變化基本呈單峰型曲線變化，但顯熱通量和潛熱通量都不如凈輻射平滑，這主要是湍流傳輸的間歇性所致。各能量分量的季節變化明顯，最大值出現在夏季，潛熱通量出現在春季，最小值均出現在冬季。

利用能量平衡比率法計算毛竹林月EBR值在0.52—1.19之間變化，1月份最小，6月份最大；月平均閉合度為0.84，年閉合度為0.85，表征毛竹林能量仍有15%的不閉合現象。而6月份出現過閉合現象，可能因為本月大量陰雨天氣造成數據被剔除較多，加之數據插補及計算方法等綜合原因造成。如果考慮顯熱通量和潛熱通量頻率損失項(≈8%)、冠層和空氣熱量儲存項(≈3%)，閉合度會得到進一步改善，全年能量閉合程度將會達到96%，但能量還是無法完全閉合，本文土壤熱通量的計算沒有考慮地表儲存熱量等，毛竹林能量閉合仍需要進一步觀測研究。對2011年毛竹林能量數據進行逐月線性回歸，平均截距、斜率處于觀測低水平，相關系數處于平均水平。波文比的月變化波動較大，在0.07—1.77之間，最大值為1月1.77，最小值為7月為0.07，月平均為0.47，年波文比為0.35，而全年顯熱通量約為潛熱通量的1/3，兩者結果基本一致。

渦動相關法是目前國際公認的測量通量的主流方法，能夠在小時、日、月、年等不同時間尺度上反映能量變化的過程，揭示大氣-植被-土壤的能量分配特征，以及環境因子的時空變化，但由于受到地形、觀測儀器、物理過程等因素的干擾，其測量結果存在一定的不確定性。該系統在復雜環境的森林生態系統中應用時，其精度還有待于進一步檢驗。因此，毛竹林的能量通量和平衡狀況的研究還有待進一步深入。

[1] Huenneke L F, Anderson J P, Remmenga M, Schlesinger W H. Desertification alters patterns of aboveground net primary production in Chihuahuan ecosystems. Global Change Biology, 2002, 8(3): 247- 264.

[2] Kosugi Y, Katsuyama M. Evapotranspiration over a Japanese cypress forest Ⅱ Comparison of the eddy covariance and water budget methods. Journal of Hydrology, 2007, 334(3/4): 305- 311.

[3] Matsumoto K, Ohta T, Nakai T, Kuwada T, Daikoku K, Shin′ichilida, Yabuki H, Kononov A V, Van der Molen M K, Kodama Y, Maximov T C, Dolman A J, Hattori S. Energy consumption and evapotranspiration at several boreal and temperate forests in the Far East. Agricultural and Forest Meteorology, 2008, 148(12): 1978- 1989.

[4] Oishi A C, Oren R, Stoy P C. Estimating components of forest evapotranspiration: A footprint approach for scaling sap flux measurements. Agricultural and Forest Meteorology, 2008, 148(11): 1719- 1732.

[5] Mauder M, Oncley S P, Vogt R, Weidinger T, Ribeiro L, Bernhofer C, Foken T, Kohsiek W, De Bruin H A R, Liu H P. The energy balance experiment EBEX- 2000 Part Ⅱ: Intercomparison of eddy-covariance sensors and post-field data processing methods. Boundary-Layer Meteorology, 2007, 123(1): 29- 54.

[6] Oncley S P, Foken T, Vogt R, Kohsiek W, De Bruin H A R, Bernhofer C, Christian A, Gorsel E V, Granta D, Feigenwinter C, Lehner I, Liebethal C, Liu H P, Mauder M, Pitacco A, Luis R, Weidinger T. The energy balance experiment EBEX- 2000. Part I: overview and energy balance. Boundary-Layer Meteorology, 2007, 123(1): 1- 28.

[7] Gao Z Q, Fan X G, Bian L G. An analytical solution to one-dimensional thermal conduction-convection in soil. Soil Science, 2003, 168(2): 99- 107.

[8] Heusinkveld B G, Jacobs A F G, Holtslag A A M, Berkowicz S M. Surface energy balance closure in an arid region: role of soil heat flux. Agricultural and Forest Meteorology, 2004, 122(1/2): 21- 37.

[9] Yang K, Wang J M. A temperature prediction-correction method for estimating surface soil heat flux from soil temperature and moisture data. Science in China Series D: Earth Sciences, 2008, 51(5): 721- 729.

[10] Falge E, Baldocchi D, Olson R, Anthoni P, Aubinet M, Bernhofer C, Burba G, Ceulemans R, Celment R, Dolman H, Granier A, Gross P, Grunwald T, Hollinger D, Jensen N O, Katul G, Keronen P, Kowalski A, Lai C T, Law B E, Meyers T, Moncrieff J, Moors E, Munger J W, Pilegaard K, Rannik U, Rebmann C, Suker A, Tenhunen J, Tu K, Verma S, Vesala T, Wilson K, Wofsy S. Gap filling strategies for long term energy flux data sets. Agricultural and Forest Meteorology, 2001, 107(1): 71- 77.

[11] Baldocchi D, Falge E, Gu L, Olson R, Hollinger D, Running S, Anthoni P, Bernhofer C, Davis K, Vans R E, Fuentes J, Goldstein A, Katul G, Law B, Lee X H, Malhi Y, Meyers T, Munger W, Oechel W, Paw K T, Pilegaard K, Schmid H P, Valentini R, Verma S, Vesala T, Wilson K, Wofsy S. FLUXNET: A new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide, water vapor, and energy flux densities. Bulletin of the American Meteorological Society, 2001, 82(11): 2415- 2434.

[12] Zhang J H, Ding Z H, Han S J, Zou C J, Zhou Y M. Turbulence regime near the forest floor of a mixed broad leaved/Korean pine forest in Changbai Mountains. Journal of Forestry Research, 2002, 13(2): 119- 122.

[13] 王旭, 周國逸, 張德強, 王春林. 南亞熱帶針闊混交林土壤熱通量研究. 生態環境, 2005, 14(2): 260- 265.

[14] 尹光彩, 王旭, 周國逸, 張德強. 鼎湖山針闊混交林土壤熱狀況研究. 華南農業大學學報, 2006, 27(3): 16- 20.

[15] 王美蓮, 崔學明, 韓鵬, 周梅. 大興安嶺原始林區土壤熱通量變化特征的初探. 內蒙古農業大學學報, 2010, 31(4): 139- 142.

[16] 岳平, 張強, 楊金虎, 李宏宇, 孫旭映, 楊啟國, 張建忠. 黃土高原半干旱草地地表能量通量及閉合率. 生態學報, 2011, 31(22): 6866- 6876.

[17] 郭曉峰, 康凌, 蔡旭暉, 朱彤. 華南農田下墊面地氣交換和能量收支的觀測研究. 大氣科學, 2006, 30(3): 453- 46.3

[18] Flerchinger G N, Cooley K R. A ten-year water balance of a mountainous semi-arid watershed. Journal of Hydrology, 2000, 237(1/2): 86- 99.

[19] 王春林, 周國逸, 王旭, 周傳艷, 于貴瑞. 鼎湖山針闊葉混交林生態系統能量平衡分析. 熱帶氣象學報, 2007, 23(6): 643- 651

[20] 王昭艷. 長江灘地抑螺防病林生態系統能量平衡與水汽通量研究. 北京: 中國林業科學研究院, 2008

[21] 岳平, 張強, 牛生杰, 成華, 王西育. 半干旱草原下墊面能量平衡特征及土壤熱通量對能量閉合率的影響. 氣象學報, 2012, 70(1): 136- 143.

[22] Lee X H. On micrometeorological observations of surface-air exchange over tall vegetation. Agricultural and Forest Meteorology, 1998, 91(1/2): 39- 49.

[23] 陳云飛, 江洪, 周國模, 孫成, 陳健. 高效經營雷竹林生態系統能量通量過程及閉合度. 應用生態學報, 2013, 24(4): 1063- 1069.

[24] Wilson K, Goldstein A, Falge E, Aubinet M, Baldocchi D, Berbigier P, Bernhofer C, Ceulemans R, Dolman H, Field C, Grelle A, Ibrom A, Law B E, Kowalski A, Meyers T, Moncrieff J, Monson R, Oechel W, Tenhunen J, Valentini R, Shashi V. Energy balance closure at FLUXNET sites. Agricultural and Forest Meteorology, 2002, 113(1/4): 223- 243.

[25] 李正泉, 于貴瑞, 溫學發, 張雷明, 任傳友, 伏玉玲. 中國通量觀測網絡(China FLUX)能量平衡閉合狀況的評價. 中國科學D輯,地球科學 2004, 34(增刊Ⅱ): 46- 56.

Energy flux and balance analysis ofPhyllostachysedulisforest ecosystem in subtropical China

SUN Cheng1, JIANG Hong1,2,*, CHEN Jian1, LIU Yuli1, NIU Xiaodong1, CHEN Xiaofeng1, FANG Chengyuan1

1ZhejiangProvincialKeyLaboratoryofCarbonCyclinginForestEcosystemsandCarbonSequestration,ZhejiangAgricultureandForestryUniversity,Hangzhou311300,China2InternationalInstituteforEarthSystemScience,NanjingUniversity,Nanjing210093,China

The atmospheric boundary layer is a passage between the atmosphere and the land surface, through the atmospheric momentum, energy and substances are transported upwards and downwards in two layers on earth surface.In the lower atmospheric planetary boundary layer turbulent transport is important connection for substance and energy exchange between the atmosphere and the earth surface, currently it is one focus of earth system science. By using open-path eddy covariance system and micro-climate instruments, the net radiation, sensitive heat flux, latent heat flux, soil heat flux, air temperature, soil temperature, and precipitation were conducted in Mao bamboo forest ecosystem of subtropical in China since 2011. But the energy flux and analysis of Mao bamboo forest ecosystem reported in the literature at home and abroad is very few, in subtropical regions to establish long-term continuous observation of the flux tower, enhancing the research is highly requirement. The diurnal and monthly variation of energy balance as well as the distribution pattern of each component were analyzed, and energy closure and Bowen ratio also were calculated. As an important index to evaluate the reliability of eddy covariance measurements, energy balance analysis has been widely accepted by the community. Using two methods of OLR(Ordinary least spuares) and EBR(Energy balance ratio), energy imbalance characteristics of the Mao bamboo forest ecosystem were systematically analyzed compared in the paper. The results showed that: the yearly net radiation of the Mao bamboo forest ecosystem was 2628.00 MJ/m2, and the sensitive flux heat, latent heat flux, and soil heat flux were 576.80 MJ/m2, 1666.77 MJ/m2, and -7.52 MJ/m2, respectively. The results suggest that soil of Mao bamboo forest was a heat source. The energy balance components had obvious seasonal pattern, diurnal variations of energy balance components existed on expected single peaks. The sensible and latent heat fluxes almost have 22.0% and 63.4% of the net radiation, respectively, indicating that latent heat flux was the main form of energy loss in this forest ecosystem. Monthly variation of the Bowen ratio is slight, and fluctuated from 0.07 to 1.77, the yearly enrgy balance closure of the Mao bamboo forest ecosystem was 0.85, and the mean monthly averge was 0.84, yet energy imclosure of 15% still existed for the energy balance of Mao bamboo forest ecosystem.

mao bamboo forest; energy flux; energy balance closure; turbulent flux; available energy; bowen ratio

國家自然科學重大基金項目(61190114); 國家“973”重點基礎研究發展規劃基金項目(2011CB302705, 2010CB950702, 2010CB428503); 國家自然科學基金項目(41171324); 科技部重大國際合作項目(20073819)

2013- 08- 27;

2014- 07- 02

10.5846/stxb201308272161

*通訊作者Corresponding author.E-mail: jianghong_china@hotmail.com

孫成, 江洪, 陳健, 劉玉莉, 牛曉棟, 陳曉峰, 方成圓.亞熱帶毛竹林生態系統能量通量及平衡分析.生態學報,2015,35(12):4128- 4136.

Sun C, Jiang H, Chen J, Liu Y L, Niu X D, Chen X F, Fang C Y.Energy flux and balance analysis ofPhyllostachysedulisforest ecosystem in subtropical China.Acta Ecologica Sinica,2015,35(12):4128- 4136.