趙子龍，張宏升＊，康 凌

1 北京大學物理學院大氣與海洋科學系，氣候與海氣實驗室，北京 100871

2 北京大學環(huán)境科學與工程學院，北京 100871

1 引 言

碳是地球上儲量最豐富的元素之一，廣泛分布于大氣、海洋、地殼沉積巖和生物體中［1］.碳循環(huán)是各碳庫之間碳的轉化和流動，包括在物理、化學和生物過程及其相互作用驅動下，各種形態(tài)的碳在各個子系統(tǒng)內部的遷移轉化過程，以及發(fā)生在子系統(tǒng)之間的交換過程［2］.陸地碳循環(huán)是全球碳循環(huán)研究中最重要的內容之一，在全球碳收支中占主導地位［3］.自然界自發(fā)進行的碳循環(huán)，在較長的時間尺度上始終處于一種均衡的狀態(tài)，然而自工業(yè)化進程以來，伴隨著全球人口的迅猛增長以及工農業(yè)的快速發(fā)展，大氣溫室氣體濃度迅速增加，其中CO2是產生溫室效應導致全球氣候變暖的最主要原因之一［4］.中國瓦里關大氣本底基準觀象臺的觀測結果表明：從1994年至2005年，大氣中CO2日平均及月平均濃度均呈逐年增加的趨勢，濃度年平均值約增長了16.48mmol·mol－1，相當于4.6%［5］.大氣中 CO2濃度的增長，將會影響全球氣候變化，從而給人類社會帶來深遠影響［6］.而研究陸地碳循環(huán)機制及其對全球變化的響應，是預測大氣CO2含量及氣候變化的基礎.近幾十年來的研究表明，北半球中高緯度的陸地生態(tài)系統(tǒng)是一個巨大的碳匯［7－13］，溫度和水分是決定陸地生態(tài)系統(tǒng)碳收支的關鍵環(huán)境因子.

大氣中的碳和水熱循環(huán)離不開大氣湍流的影響.大量的觀測實驗使人們對平坦均一下墊面的溫濕度湍流結構有了較為系統(tǒng)的認識，得到了許多有益的結論［14－16］：溫度和濕度的歸一化標準差與穩(wěn)定度參數z／L滿足－1／3次方律關系；溫度和濕度能譜在不穩(wěn)定條件下，高頻段滿足－2／3次方律，中頻段有一平緩區(qū)域；溫度和濕度分別與水平、垂直方向風速之間的互譜符合MOS理論，在慣性副區(qū)符合－4／3次方律；Roth等［17］認為低頻段溫度能譜集中在一個相對狹窄的范圍內，基本不隨z／L變化.

我國西北戈壁地區(qū)地域遼闊，大氣中CO2及水熱通量等有著獨特的演變規(guī)律［18］，對全球碳循環(huán)的貢獻不可忽略.以往的實驗研究主要集中在下墊面為森林、草原和農田等植被較為豐富的地區(qū)［19－23］，對干燥的戈壁下墊面的溫度、濕度及CO2濃度的湍流特征的報道并不多見，尤其是CO2濃度的研究在文獻中較少涉及.

本文利用大氣邊界層實驗觀測資料，采用渦動相關法分析研究了我國西北戈壁地區(qū)低層大氣碳收支和水熱輸送特征、湍流統(tǒng)計特征和湍流譜結構，試圖得到戈壁下墊面地區(qū)的碳及水熱循環(huán)規(guī)律.

2 理論依據

根據近地面層相似性理論，不穩(wěn)定條件下，溫度歸一化標準差σT／T＊與大氣穩(wěn)定度參數z／L滿足－1／3次指數關系：

其中，T＊＝－，α 為常數.Wyngaard等［24］利用美國中部Kansas平坦草地下墊面的觀測資料得到系數α＝0.95.

類似地，不穩(wěn)定條件下，比濕和CO2濃度等標量的歸一化標準差σq／q＊與大氣穩(wěn)定度參數z／L也滿足－1／3次方律的關系，即

其中，q＊＝－，γ 為常數.Hogstorm 等［25］給出γ＝1.04.

溫度譜的慣性副區(qū)的形式可表達為［26］：

其中，β＝0.8為普適常數，ε為湍流動能耗散率，N為湍流熱能耗散率，k為波數.利用Taylor假設將波數轉化為頻率，引入無因次耗散率：

和無因次頻率f＝nz／u－（下同）進行歸一化，得到

歸一化后的湍流溫度能譜密度在慣性副區(qū)與無因次頻率f的－2／3次方成正比.Kaimal等［27］在Kansas實驗給出中性層結條件下溫度能譜密度的擬合曲線為：濕度和CO2濃度能譜密度存在類似關系.

物理量之間的互譜常常用來分析不同尺度的渦旋對相應通量的貢獻情況.Wyngaard等［28］提出互譜在慣性副區(qū)符合－4／3次方律，這一結論已經得到了眾多實驗的支持.Kaimal等［27］在Kansas實驗中給出了中性層結條件下風速互譜近似滿足：

垂直風速和溫度互譜的曲線滿足［27］：

3 資料獲取與處理

本文利用我國西北玉門鎮(zhèn)附近戈壁地區(qū)的大氣邊界層的觀測數據.該實驗站位于40°09′27.3″N，97°17′58.6″E，海拔高度為1200m；實驗站東側為荒涼的戈壁灘，西距小型居民生活區(qū)約3km，南距祁連山約30km，北距北山山地約20km，地形自南往北緩慢傾斜，東北方向有極小的坡度下降，坡度約為1°，間有自然沖溝和小起伏，局部地形是平坦開闊的戈壁地區(qū).全年主導風向為東風，西風次之.

觀測時間為2006年6月27日—7月27日，觀測平臺為102m氣象鐵塔.本文涉及的實驗觀測項目主要包括H2O和CO2的大氣湍流及4層的風速、溫度和濕度梯度.在30m高度架設Campbell公司生產的CSAT3型超聲風溫儀和Li－Cor公司生產的Li－7500型開路式 CO2／H2O 分析儀.在10m，30m，60m和100m高度架設美國MetOne公司生產的034B型風杯風速計，同時安裝有溫度和濕度傳感器.

湍流觀測自動、連續(xù)，采樣頻率為10Hz，分析時使用每小時的前半小時結果，并對湍流資料進行了嚴格的篩選，即：去除了風向與超聲探頭迎風軸朝向的夾角大于60°、風向與水平面夾角大于5°、風速小于2m·s－1、摩擦速度小于0.01m·s－1、熱通量小于5W·m－2以及其它非定常的數據組.塔層風速、溫度和濕度梯度觀測，為10min的平均值，取30min滑動平均.考慮不同天氣狀況對太陽輻射有較大影響，從而增加了影響碳及水熱通量變化的因素，為避免天氣狀況對分析碳收支及水熱通量規(guī)律可能引起的不確定性，本文只選取晴好少云天氣（6月27日到7月7日）的觀測數據進行分析，并剔除有降水的觀測時次.

4 結果與分析

4.1 平均量特征

圖1給出2006年6月27日—7月7日戈壁地區(qū)近地層感熱通量隨時間變化特征.受太陽輻射的影響，戈壁地區(qū)白天熱量自下向上輸送，感熱通量在正午前后達到極大值；夜間，熱量由上向下輸送，在日出前后達到極小值.30m高度感熱通量的最大值為463.8W·m－2，出現在7月3日13時，最小值為－53.2W·m－2，出現在6月30日00時.

圖2給出了2006年6月27日—7月7日戈壁地區(qū)近地層潛熱通量隨時間的變化特征.盡管戈壁地區(qū)空氣干燥，水汽含量低，潛熱通量的獲取受到儀器精度的限制，但還是可以從圖2看出，30m高度處的潛熱通量在午后向地面?zhèn)鬏?這種現象產生的原因是在戈壁干旱地區(qū)，地表干燥，空氣濕度普遍較小，當空氣的水平運動帶來水汽時，會向地面進行水汽輸送，形成逆濕，從而造成潛熱通量的負值.

圖1 2006年6月27日—7月7日30m高度戈壁地區(qū)感熱通量隨時間變化Fig.1 Temporal variations of the sensible heat flux at the 30mheight over the Gobi Desert from June 27to July 7，2006

圖3給出了2006年6月27日—7月7日戈壁地區(qū)近地層CO2通量的時間序列.CO2通量平均值為－0.199mg·m－2·s－1，整體表現是碳匯；絕對值的最大值出現在6月29日16時，數值為－1.10mg·m－2·s－1，與劉冉等［29］2004年在鹽生荒漠下墊面觀測到CO2通量最大值為－1.12mg·m－2·s－1的結果相近，是平均值的5.5倍.白天正午前后有兩個向下輸送的峰值，出現的時間比CO2濃度峰值約滯后4～5小時.

圖4給出了2006年6月27日到2006年7月7日戈壁地區(qū)近地層CO2通量的平均日變化.可見，白天CO2通量向下極大值為－0.74mg·m－2·s－1，與劉輝志等［22］在吉林通榆觀測到的退化草原下墊面在生長季的結果為－0.8mg·m－2·s－1相近；夜間向上，其數值小于白天.這與生長期的森林、草地等［30－32］植被區(qū)域下墊面的觀測結果大致相同，說明夏季戈壁下墊面的土壤具有一定的碳吸收功能［33－34］.

圖4 2006年6月27日—7月7日30m高度戈壁地區(qū)CO2通量平均日變化Fig.4 Diurnal variation of the CO2flux at the 30m height over the Gobi Desert averaged from June 27to July 7，2006

4.2 溫度、濕度和CO2的湍流統(tǒng)計特征

圖5a給出了2006年6月27日－7月7日戈壁下墊面近地層溫度歸一化標準差σT／T＊隨穩(wěn)定度參數z／L的變化特征.可見，不穩(wěn)定層結下溫度歸一化標準差與穩(wěn)定度參數z／L呈－1／3次方律關系，與Wyngaard等［24］的結果相符合.擬合公式為：

σT／T＊＝1.00（－z／L）－1／3，z／L＜0 （10）在穩(wěn)定層結下，數據的離散性較大.

圖5b和圖5c分別給出了2006年6月27日—7月7日戈壁下墊面近地層比濕歸一化標準差σq／q＊和CO2濃度歸一化標準差σc／c＊隨穩(wěn)定度參數z／L的變化特征，直線是 Hogstorm等［25］的結果.可見，比濕歸一化標準差和CO2濃度歸一化標準差隨z／L變化的關系較為離散，究其原因，一是由于戈壁地區(qū)水汽含量低，濕度觀測相對誤差較大，戈壁地區(qū)CO2的變化較小，其觀測精度不足；另一種可能是由于比濕與CO2濃度是被動標量［35］，其歸一化標準差的變化規(guī)律與溫度等主動標量不同［35－36］，對于被動標量，McBean等［35］建議使用經驗公式Lq和Lc替代L來計算穩(wěn)定度參數，但目前還沒有公認的經驗公式，故本文仍使用了傳統(tǒng)的穩(wěn)定度參數z／L.

4.3 能譜特征

圖6a，6b，6c分別給出了2006年6月27日—7月7日近地層不同穩(wěn)定度條件下的溫度、濕度和CO2濃度的歸一化能譜密度隨無因次頻率的變化，圖中實線為Kaimal等［27］給出的Kansas實驗溫度能譜密度中性層結下的擬合曲線.可以看到，溫度、濕度和CO2濃度的能譜在慣性副區(qū)均滿足－2／3次方律關系，溫度和CO2濃度的特征相似，而濕度則有較大的差別.

濕度能譜密度的峰值頻率略低，說明戈壁地區(qū)水汽輸送貢獻最大的渦旋尺度略大.在低頻段，不同穩(wěn)定度的情況下，濕度能譜的譜線較為分散，而溫度和CO2濃度能譜的譜線則集中在很窄的范圍內，難以區(qū)分，說明大氣穩(wěn)定度對戈壁地區(qū)溫度和CO2濃度的影響小于濕度，低頻段的溫度能譜和CO2濃度能譜還受其它尺度的控制［16］.

圖7a，7b，7c分別給出2006年6月27日—7月7日近地層戈壁地區(qū)不同穩(wěn)定度下水平方向風速與溫度、水汽和CO2濃度的歸一化互譜特征，圖中實線為Kaimal等［27］給出的Kansas實驗中性層結下wθ互譜擬合曲線.

圖5 2006年6月27日－7月7日30m高度戈壁下墊面溫度、濕度、CO2濃度歸一化標準差σT／T＊ 、σq／q＊ 、σc／c＊ 隨穩(wěn)定度參數z／L的變化特征Fig.5 Distribution of the normalized standard deviation of T（σT／T），q（σq／q＊ ）and c（σc／c＊ ）against z／Lat the 30mheight over the Gobi Desert from June 27to July 7，2006

圖7a所示的uθ互譜的峰值頻率與Kansas實驗［27］中性層結條件下的wθ互譜的峰值頻率相比數值略大，不同穩(wěn)定度下所有譜線的峰值頻率非常接近，約為0.2.低頻段uθ互譜所有譜線都集中在較窄的范圍內，隨著無因次頻率的降低，譜線的下降比uc互譜更快，說明湍渦尺度大小對水平方向熱量輸送的影響比對水平方向CO2通量輸送的影響更大.uθ互譜和uc互譜非常的相似，但與uq互譜區(qū)別較大.

uq互譜在低頻區(qū)域不同穩(wěn)定度情況下譜線分布在一個比較寬的范圍內，峰值頻率的分布也比較分散，從0.1到0.2之間，說明在不同穩(wěn)定度的條件下，影響水平方向水汽通量輸送的渦旋尺度有較大的差別.慣性副區(qū)的區(qū)間相比uc互譜較小.

由圖7c可以看到，不同穩(wěn)定度參數下，uc互譜在慣性副區(qū)的曲線滿足－4／3次方律，不同穩(wěn)定度參數情況下的譜線分布較為集中，峰值頻率、慣性副區(qū)起始位置、低頻區(qū)下降特征等都與Kaimal的平坦草原下墊面wθ的結果非常接近.

圖8a，8b，8c分別給出了2006年6月27日—7月7日戈壁地區(qū)近地層不同穩(wěn)定度下垂直方向風速與溫度、水汽和CO2濃度的歸一化互譜特征，圖中實線為Kaimal等［27］給出的Kansas實驗中性層結下互譜擬合曲線.由圖7和圖8可見，各個互譜在慣性副區(qū)都符合－4／3次方律的關系.

由圖8a，與Kansas實驗代表的平坦草原下墊面的互譜特征相比，無論是慣性副區(qū)的區(qū)間范圍還是峰值頻率、峰值大小、低頻段下降特征等兩者都很相似，說明平坦戈壁下墊面和平坦草原下墊面在垂直方向感熱輸送上存在很強的相似性.分析圖8a和圖7a，可以看出垂直和水平方向湍流交換對熱量輸送方面的影響相近.

對比圖8b所示的wq互譜特征與圖7b所示的uq互譜特征，在慣性副區(qū)區(qū)間大小、區(qū)間位置和峰值頻率大小等特征上差別不大，但wq互譜低頻區(qū)和峰值大小有所降低，說明在夏季戈壁地區(qū)，水汽通量輸送主要受水平方向湍流交換的影響.wq互譜的峰值頻率分布也比較分散，低頻區(qū)不同穩(wěn)定度下譜線的分布比uq互譜要集中.

由圖8c和圖7c可見，在慣性副區(qū)區(qū)間大小和位置、低頻段下降趨勢、峰值頻率大小等特征上wc互譜和uc互譜具有很強的相似性，但與uc互譜相比，wc互譜峰值大小略低，說明戈壁地區(qū)CO2通量輸送受水平方向湍流的影響較大.

5 結 論

本文利用渦動相關法分析了2006年夏季6月27日—7月27日在我國西北玉門鎮(zhèn)附近戈壁地區(qū)碳和水熱循環(huán)特征，以及大氣湍流特征，主要結果如下：

（1）受太陽輻射、下墊面特征和干旱氣候的影響，戈壁地區(qū)近地層感熱通量較大，日變化規(guī)律明顯，在正午前后達到極大值.觀測期間感熱通量的極大值通常出現于13時.

圖8 2006年6月27日—7月7日30m高度戈壁地區(qū)不同穩(wěn)定度下wθ、wq、wc歸一化互譜特征（f＝nz／ˉu）Fig.8 Normalized wθ，wq，wc cospectra under different stabilities at the 30mheight over the Gobi Desert from June 27to July 7，2006

（2）夏季戈壁地區(qū)近地層潛熱通量相對較小，在白天大部分時間內僅為感熱通量的10%左右，這與戈壁地區(qū)干旱的氣候條件有關.然而實驗觀測顯示午后潛熱通量時常為負值，說明此時潛熱從空中向地面輸送，主要是受平流帶來的水汽影響.

（3）夏季戈壁地區(qū)近地層CO2濃度的平均值為557mg·m－3.CO2通量白天向下傳輸，最大值為1.10mg·m－2·s－1，夜間向上傳輸，其值很小，均不超過0.05mg·m－2·s－1.平均值為－0.199mg·m－2·s－1，整體表現為碳匯.CO2通量的日變化規(guī)律，與生長期的森林、草地等有植被區(qū)域下墊面的觀測結果大致相同.

（4）通過對夏季戈壁地區(qū)近地層的溫度、比濕、CO2濃度的湍流統(tǒng)計特征分析可知，在慣性副區(qū)，溫度歸一化標準差σT／T＊與穩(wěn)定度參數z／L滿足－1／3次方律的關系；而比濕歸一化標準差σq／q＊和CO2濃度歸一化標準差σc／c＊與穩(wěn)定度參數z／L的關系表現離散，可能的原因之一是比濕和CO2濃度這樣的被動標量的規(guī)律需要用與溫度T這樣的主動標量不同的穩(wěn)定度參數來進行描述，然而具體的參數化公式還有待研究，同時也不排除測量儀器的精度有限.

（5）夏季戈壁地區(qū)近地層的溫度、濕度、CO2濃度的能譜和uθ、uq、uc、wθ，wq，wc的互譜特征顯示：在慣性副區(qū)，各能譜均滿足－2／3次方律；各互譜均滿足－4／3次方律；溫度能譜與CO2濃度能譜相似，與Kaimal等［27］給出的Kansas實驗溫度能譜相近，而濕度能譜則表現出峰值頻率略低和低頻段譜線較為分散的特征；相對于uq互譜來說，uc互譜與uθ互譜相似，主要表現為不同穩(wěn)定度情況下譜線的峰值頻率較為接近，低頻段譜線都集中在較窄的范圍內；wc互譜與wθ互譜相似，在慣性副區(qū)區(qū)間范圍、峰值頻率和峰值大小等特征上都與平坦草原下墊面的互譜wθ相近；對比uc互譜與wc互譜、uq互譜與wq互譜、uθ互譜與wθ互譜，可知水汽和CO2的輸送主要受水平方向湍流的影響，而水平和垂直方向湍流都對熱量輸送起重要作用.

（References）

［1］翁金桃.碳酸鹽巖在全球碳循環(huán)過程中的作用.地球科學進展，1995，10（2）：154－158.Weng J T.The effect of carbonate rocks on global carbon cycle.Advan.Earth Sci.（in Chinese），1995，10（2）：154－158.

［2］陳泮勤，黃耀，于貴瑞.地球系統(tǒng)碳循環(huán).北京：科學出版社，2008.Chen P Q，Huang Y，Yu G R.Carbon Cycle of Earth System（in Chinese）.Beijing：Science Press，2008.

［3］Watson R T，Verardo D J.Land Use，Land－Use Change and Forestry.United Kingdom：Cambridge University Press，2000.

［4］Bolin B，Degens E T.The Global Biogeochemical Carbon Cycles.New York：John Wiley ＆Sons，1979.

［5］趙玉成，溫玉璞，德力格爾等.青海瓦里關大氣CO2本底濃度的變化特征.中國環(huán)境科學，2006，26（1）：1－5.Zhao Y C，Wen Y P，Deliger，et al.Change characteristics of atmospheric CO2background concentration in Waliguan Qinghai.China Environ.Sci.（in Chinese），2006，26（1）：1－5.

［6］Dolman A J，van der Werf G R，van der Molen M K，et al.A carbon cycle science update since IPCC AR－4.Ambio，2010，39（5－6）：402－412.

［7］Schimel D S，House J I，Hibbard K A，et al.Recent patterns and mechanisms of carbon exchange by terrestrial ecosystems.Nature，2001，414（6860）：169－172.

［8］Pacala S W，Hurtt G C，Baker D，et al.Consistent land－and atmosphere－based US carbon sink estimates.Science，2001，292（5525）：2316－2320.

［9］Dixon R K，Solomon A M，Brown S，et al.Carbon pools and flux of global forest ecosystems.Science，1994，263（5144）：185－190.

［10］Battle M，Bender M L，Tans P P，et al.Global carbon sinks and their variability inferred from atmospheric O2andδ13C.Science，2000，287（5462）：2467－2470.

［11］Fan S，Gloor M，Mahlman J，et al.A large terrestrial carbon sink in North America implied by atmospheric and oceanic carbon dioxide data and models.Science，1998，282（5388）：442－446.

［12］Schimel D，Melillo J，Tian H Q，et al.Contribution of increasing CO2and climate to carbon storage by ecosystems in the United States.Science，2000，287（5460）：2004－2006.

［13］李玉強.沙漠化對沙地生態(tài)系統(tǒng)碳氮儲量與平衡及能量的影響［博士論文］.北京：中國科學院研究生院，2006.Li Y Q.Impacts of desertification on carbon and nitrogen storages，carbon and nitrogen balances，and biomass energy in Horqin sand land ecosystems［Ph.D.thesis］（in Chinese）.Beijing：Graduate University of Chinese Academy of Sciences，2006.

［14］Haugen D A.Workshop on micrometeorology.Bull.Amer.Meteor.Soc.，1973，54（2）：108－110.

［15］Panofsky H A，Dutton J A.Atmospheric Turbulence：Models and Methods for Engineering Applications.New York：A Wiley－Interscience Pub.，1984.

［16］劉明星，張宏升，宋星灼等.不同下墊面溫度和濕度湍流譜特征研究.北京大學學報 （自然科學版），2008，44（3）：391－398.Liu M X，Zhang H S，Song X Z，et al.Spectral characteristics of atmospheric turbulence over various surface conditions.Acta Scientiarum Naturalium Universitatis Pekinensis（in Chinese），2008，44（3）：391－398.

［17］Roth M，Oke T R.Turbulent transfer relationships over an urban surface.1.spectral characteristics.Quart.J.Roy.Meteor Soc.，1993，119（513）：1071－1104.

［18］何玉斐.戈壁下墊面塔層湍流輸送與碳通量特征研究［碩士論文］.北京：北京大學物理學院大氣與海洋科學系，2008.He Y F.Study on characteristics of turbulent transfer and carbon dioxide fluxes in tower layer over Gobi Desert area［Master′s thesis］（in Chinese）.Beijing：Department of Atmospheric and Oceanic Sciences，School of Physics of Peking University，2008.

［19］王文杰，于景華，毛子軍等.森林生態(tài)系統(tǒng)CO2通量的研究方法及研究進展.生態(tài)學雜志，2003，22（5）：102－107.Wang W J，Yu J H，Mao Z J，et al.Study method and research advance in forest CO2flux.Chinese J.Ecol.（in Chinese），2003，22（5）：102－107.

［20］李婧，劉樹華，茅宇豪等.不同生態(tài)系統(tǒng)CO2通量和濃度特征分析研究.地球物理學報，2006，49（5）：1298－1307.Li Q，Liu S H，Mao Y H，et al.Characteristics of CO2flux and concentration in different ecosystems.Chinese J.Geophys.（in Chinese），2006，49（5）：1298－1307.

［21］陸龍驊，程彥杰，卞林根等.長江三角洲典型稻作區(qū)近地層二氧化碳等湍流通量的觀測研究.地球物理學報，2003，46（6）：751－759.Lu L H，Cheng Y J，Bian L G，et al.A study on the turbulence fluxes of the surface layer CO2，sensitive and latent exchange over the typical rice field，Changjiang delta.Chinese J.Geophys.（in Chinese），2003，46（6）：751－759.

［22］劉輝志，涂鋼，董文杰等.半干旱地區(qū)地氣界面水汽和二氧化碳通量的日變化及季節(jié)變化.大氣科學，2006，30（1）：108－118.Liu H Z，Tu G，Dong W J，et al.Seasonal and diurnal variations of the exchange of water vapor and CO2between the land surface and atmosphere in the semi－arid area.Chinese J.Atmos.Sci.（in Chinese），2006，30（1）：108－118.

［23］姜紀峰，延曉冬，黃耀等.半干旱區(qū)農田和草地與大氣間二氧化碳和水熱通量的模擬研究.氣候與環(huán)境研究，2006，11（3）：413－424.Jiang J F，Yan X D，Huang Y，et al.Simulation of CO2and sensible／latent heat fluxes exchange between land surface and atmosphere over cropland and grassland in semi－arid region.Climatic and Environmental Research （in Chinese），2006，11（3）：413－424.

［24］Wyngaard J C，Cote O R，Izumi Y.Local free convection，similarity，and the budgets of shear stress and heat flux.J.Atmos Sci.，1971，28（7）：1171－1182.

［25］Hogstorm U， Smedman－Hogstorm A S. Turbulence mechanisms at an agricultural site.Bound－Layer Meteor.，1974，7（3）：373－389.

［26］Corrsin S. On the spectrum of isotropic temperature fluctuations in an isotropic turbulence.J.Appl.Phys.，1951，22（4）：469－473.

［27］Kaimal J C， Wyngaard J C，Izumi Y，et al.Spectral characteristics of surface－layer turbulence.Quart.J.Roy.Meteor.Soc.，1972，98（417）：563－589.

［28］Wyngaard J C，Coté O R.Cospectral similarity in the atmospheric surface layer.Quart.J.Roy.Meteor.Soc.，1972，98（417）：590－603.

［29］劉冉，王勤學，唐立松等.鹽生荒漠地表水熱與二氧化碳通量的季節(jié)變化及驅動因素.生態(tài)學報，2009，29（1）：67－75.Liu R，Wang Q X，Tang L S，et al.Seasonal variation in water，heat and CO2fluxes and its driving forces over a saline desert.Acta Ecol.Sin.（in Chinese），2009，29（1）：67－75.

［30］高志球，卞林根，陸龍驊等.水稻不同生長期稻田能量收支、CO2通量模擬研究.應用氣象學報，2004，15（2）：129－140.Gao Z Q，Bian L G，Lu L H，et al.Modeling of energy and CO2fluxes during different growth periods over rice field by using SIB2.J.Appl.Meteor.Sci.（in Chinese），2004，15（2）：129－140.

［31］卞林根，高志球，陸龍驊等.長江下游農業(yè)生態(tài)區(qū)CO2通量的觀測試驗.應用氣象學報，2005，16（6）：828－834.Bian L G，Gao Z Q，Lu L H，et al.Measurements of CO2fluxes over two different underlying surfaces in an agricultural ecosystem over lower basins of the Yangtze.J.Appl.Meteor.Sci.（in Chinese），2005，16（6）：828－834.

［32］徐玲玲，張憲洲，石培禮等.青藏高原高寒草甸生態(tài)系統(tǒng)凈二氧化碳交換量特征.生態(tài)學報，2005，25（8）：1948－1952.Xu L L，Zhang X Z，Shi P L，et al.Net ecosystem carbon dioxide exchange of alpine meadow in the Tibetan plateau from august to october.Acta Ecol.Sin.（in Chinese），2005，25（8）：1948－1952.

［33］趙哈林，李玉強，周瑞蓮.沙漠化對科爾沁沙質草地土壤呼吸速率及碳平衡的影響.土壤學報，2009，46（5）：809－816.Zhao H L，Li Y Q，Zhou R L.Effects of desertification on soil respiration rate and carbon balance in Horqin sandy grassland.Acta Pedol.Sin.（in Chinese），2009，46（5）：809－816.

［34］Pan G，Guo T.Pedogenic carbonate of aridic soil in China and its significance in carbon sequestration in terrestrial systems.∥Lal R，Kimble J，Eswaran H eds.Global Climate Change and Pedogenic Carbonates.Florida：Lewis Publishers Inc，2000.

［35］McBean G A.The variations of the statistics of wind，temperature and humidity fluctuations with stability.Bound－Layer Meteor.，1971，1（4）：438－457.

［36］Roth M.Turbulent transfer relationships over an urban surface.II：Integral statistics.Quart.J.Roy.Meteor.Soc.，1993，119（513）：1105－1120.