黃土高原復雜地形上高質(zhì)量湍流通量數(shù)據(jù)獲取方法

苑廣輝， 張鐳， 李遙， 梁捷寧

半干旱氣候變化教育部重點實驗室, 蘭州大學大氣科學學院, 蘭州　730000

?

黃土高原復雜地形上高質(zhì)量湍流通量數(shù)據(jù)獲取方法

苑廣輝， 張鐳*， 李遙， 梁捷寧

半干旱氣候變化教育部重點實驗室, 蘭州大學大氣科學學院, 蘭州730000

摘要利用蘭州大學半干旱氣候與環(huán)境觀測站(Semi-Arid Climate and Environment Observatory of Lanzhou University，SACOL)湍流觀測資料，分析了二次坐標旋轉(zhuǎn)(double rotation，DR)、平面擬合(planar fit，PF)和分風區(qū)平面擬合(fetch planar fit，F(xiàn)PF)在復雜地形上的適用性，總結出一套適用于SACOL的總體湍流特征參數(shù)化方案．經(jīng)過超聲虛溫訂正、坐標旋轉(zhuǎn)、空氣密度脈動訂正以及平穩(wěn)性檢驗、總體湍流特征檢驗、總體質(zhì)量分級處理，摩擦速度(u*)、感熱通量、潛熱通量、CO2通量高質(zhì)量數(shù)據(jù)所占比例分別為45%～62%、66%～68%、62%～65%、52%～54%．采用DR得到的高質(zhì)量數(shù)據(jù)比例與采用PF相比，u*提高了17%，后三種通量略降低2%～3%．PF和FPF兩種結果的差別主要體現(xiàn)在u*上，只考慮主導風向數(shù)據(jù)DR得到的u*質(zhì)量仍最好．綜合兼顧數(shù)據(jù)質(zhì)量和計算工作量，在復雜地形上處理湍流觀測資料的最優(yōu)坐標旋轉(zhuǎn)方法是DR．

關鍵詞黃土高原； 渦動相關； 通量訂正； 質(zhì)量控制； 分風區(qū)平面擬合

1引言

在近地層，湍流是大氣運動的主要形式，引起各物理屬性在地-氣之間的交換輸送，地表熱量、水汽、CO2等物質(zhì)的傳輸過程受到廣泛關注(Jiang et al.，2013；劉樹華等，2005a, 2005b；徐自為等，2008)，截至2014年4月全球通量觀測站網(wǎng)絡(FLUXNET)已有683個站點(http:∥fluxnet.ornl.gov/)．渦動相關法(eddy covariance，EC)可以用來長期連續(xù)觀測CO2、水汽和熱通量，廣泛應用于地氣間通量交換的測量，已成為FLUXNET測量通量的主要技術手段(Baldocchi et al.，2001)．EC建立在一定的假設之上，如湍流平穩(wěn)、下墊面水平均勻、湍流充分發(fā)展、存在常通量層等(Foken and Wichura，1996)，而在黃土高原復雜下墊面上，這些假設條件通常難以得到滿足，地形起伏對EC測量通量有重要影響(王介民等，2007)．如何準確地計算近地層通量一直受到地球科學界的重視(丁一匯，1997；姜海梅等，2013；劉樹華等，2005c, 2009)，已有不少學者對湍流資料處理和質(zhì)量控制進行研究：Mauder等(2006)通過對LITFASS-2003實驗中14個測點的數(shù)據(jù)進行總體質(zhì)量檢驗，發(fā)現(xiàn)80%的潛熱通量為高質(zhì)量數(shù)據(jù)；G?ckede等(2008)將拉格朗日隨機足跡模型與質(zhì)量評價相結合應用到CarboEurope的25個森林測站來研究通量測量的空間代表性、儀器以及坐標旋轉(zhuǎn)方法對湍流通量的影響；王少影等(2009)討論了在綠洲和戈壁測站不同坐標旋轉(zhuǎn)方法對摩擦速度、感熱、潛熱通量和湍流平穩(wěn)性檢驗以及總體湍流特征檢驗的影響；朱治林等(2004)討論了在非平坦、非均勻下墊面上，儀器安裝不垂直在EC計算通量時可能產(chǎn)生的誤差及其校正方法，對比分析了不同坐標旋轉(zhuǎn)方法的校正效果和適用條件；在平坦均勻下墊面上，諶志剛等(2008)的分析發(fā)現(xiàn)，平面擬合方法優(yōu)于三次坐標旋轉(zhuǎn)方法；姜海梅等(2012)運用EBEX-2000實驗的湍流、凈輻射和土壤觀測資料，運用渦動相關方法分析了非均勻灌溉引起的熱內(nèi)邊界層發(fā)展條件下近地層感熱、潛熱通量特征，并對有無灌溉兩種條件下的能量閉合度進行了對比分析．黃土高原半干旱地區(qū)面積廣闊，其陸氣相互作用不僅對西北干旱氣候的形成以及東亞季風環(huán)流有著不可忽視的影響，同時對全球氣候和大氣環(huán)流的變化也可能產(chǎn)生比較重要的作用(楊文治和邵明安，2000)．這一區(qū)域的陸氣相互作用問題已經(jīng)成為我國迫切需要研究的重要基礎性科學問題之一．在黃土高原復雜下墊面上，地形起伏，溝壑縱橫，植被稀疏，湍流更為復雜，而關于湍流資料處理和質(zhì)量控制的研究較少，獲取高質(zhì)量通量數(shù)據(jù)對深入認識該地區(qū)地氣交換特征非常必要．

本文主要討論在黃土高原復雜地形上各種湍流資料處理方法對EC觀測結果的影響，應用超聲虛溫訂正、坐標旋轉(zhuǎn)、空氣密度脈動訂正以及平穩(wěn)性檢驗、總體湍流特征檢驗、總體質(zhì)量分級處理方法，對蘭州大學半干旱氣候與環(huán)境觀測站(Semi-Arid Climate and Environment Observatory of Lanzhou University，SACOL)2008年11月1日至21日的10 Hz原始湍流觀測數(shù)據(jù)進行了質(zhì)量控制，有效地減少了因地形等原因造成的湍流觀測數(shù)據(jù)附加誤差，使其湍流觀測數(shù)據(jù)質(zhì)量得到提高．根據(jù)SACOL風向特點使用分風區(qū)平面擬合方法，著重關注二次坐標旋轉(zhuǎn)、平面擬合及分風區(qū)平面擬合的適用性，進而分析SACOL所代表的黃土高原復雜地形上近地層湍流特征．

2觀測資料

SACOL建立于2005年，位于海拔高度為1965.8 m的萃英山頂(35.95°N，104.13°E)，距蘭州市中心約48 km，位于中國黃土高原半干旱區(qū)，下墊面為典型的黃土高原殘塬地貌，溝壑縱橫，土壤為第四紀黃土風蝕形成的灰鈣土．擁有國際先進的觀測儀器，是繼中國科學院吉林通榆站之后，第二個由我國自主建設的半干旱區(qū)長期觀測站．已被批準加入國際協(xié)同觀測計劃項目(The Coordinated Enhan-ced Observing Period)，并作為此計劃的全球協(xié)同加強觀測站之一(Huang et al.，2008；梁捷寧等，2014)．

通量觀測場地較為平坦，東西方向約200 m，南北方向約1000 m．地表為長芒草、冷蒿、賴草等短小植被覆蓋，冬季地表植被高約0.10 m(Huang et al.，2008；梁捷寧等，2014)．

渦動相關系統(tǒng)分別利用三維超聲風溫儀(CSAT3，Campbell)測量u，v，w三維風速和超聲虛溫Ts，開路紅外CO2/H2O氣體分析儀(LI7500，LI-COR)測量CO2和H2O密度，感應器采樣頻率為10 Hz，觀測高度為3 m，數(shù)據(jù)采集器(CR5000，Campbell)．

2008年夏季對觀測儀器進行校準，此后數(shù)據(jù)可信度較高．冬季下墊面植被低矮稀疏，對湍流通量影響較小，故選取2008年11月1日至21日的10 Hz原始觀測數(shù)據(jù)，資料完好率高，由于11月3日數(shù)據(jù)中有大量異常值出現(xiàn)，所以剔除11月3日全天數(shù)據(jù)．

3分析方法

3.1渦動相關方法

渦動相關方法由澳大利亞微氣象學家Swinbank于1951年提出(Swinbank，1951)，某物理量X的垂直湍流通量定義為

(1)

(2)

(3)

(4)

(5)

3.2.1野點剔除和插補

3.2.2超聲虛溫訂正

超聲風速計測得的溫度為超聲虛溫Ts，受濕度影響．計算溫度和感熱通量時應考慮濕度訂正．Schotanus等(1983)提出超聲虛溫的訂正方法，Aubinet等(2011)對其進行簡化，

(6)

(7)

3.2.3坐標旋轉(zhuǎn)

坐標旋轉(zhuǎn)方法主要有二次坐標旋轉(zhuǎn)(double rotation，DR)、三次坐標旋轉(zhuǎn)(triple rotation，TR)和平面擬合(planar fit，PF)，三次坐標旋轉(zhuǎn)在計算應力時誤差較大，已不推薦使用(Kaimal and Finnigan，1994)．

(1)二次坐標旋轉(zhuǎn)

(8)

u1=u0cosγ+v0sinγ,

(9)

v1=-u0sinγ+v0cosγ,

(10)

w1=w0,

(11)

(12)

u2=u1cosα+w1sinα,

(13)

w2=-u1sinα+w1cosα,

(14)

v2=v1,

(15)

(2)平面擬合

PF(Wilczak et al.，2001)是根據(jù)一個較長的時段，本文選取2008年1日至21日(3日除外)20天的數(shù)據(jù)，確定一個與地面平行的平均風場，并將各時次的u，v，w旋轉(zhuǎn)到該平面上．利用u，v，w的30 min平均值計算方程組，

(16)

(17)

(19)

(20)

3.2.4WPL訂正

溫度和濕度擾動會引起關注的微量氣體濃度變化，影響湍流通量測量，對水汽通量(Fv)和CO2通量(Fc)需進行WPL訂正以消除密度效應的影響．

(21)

(22)

其中μ=1.6，為干空氣與水汽分子量之比，σ=ρv/ρa，為水汽和干空氣密度之比，ρa為干空氣密度(Webb et al.，1980)．WPL訂正在水汽通量上是加一個感熱通量訂正項，在CO2通量上是加一個水汽通量和一個感熱通量訂正．

4結果與討論

4.1超聲虛溫訂正和WPL訂正

表1　超聲虛溫訂正結果(百分比為)

圖1為WPL訂正前后潛熱通量和碳通量對比．由圖1a, 1c可知，潛熱通量日變化為白天高，最高可達120 W·m-2，夜間基本低于5 W·m-2．WPL訂正使?jié)摕嵬靠傮w增加了7.4%，白天增加5.6%，夜間降低3.4%，這與感熱輸送方向的晝夜差異有關．白天，向上輸送的感熱加熱空氣，使觀測高度處水汽密度減小，導致EC對水汽通量的觀測值偏低，夜間則相反．由圖1b, 1d可看出，向下的CO2通量在正午前后達到最高，向上的最高值出現(xiàn)在0時左右．白天，植物光合作用吸收大氣中的CO2，形成向下的CO2通量，并于正午達到最強；夜間，植物呼吸作用釋放出CO2，使得CO2通量向上傳遞，形成一峰一谷的日變化趨勢．總體上，WPL訂正使CO2通量降低了72.2%，白天向下的CO2通量降低64.7%，夜間向上的CO2通量降低23.4%．

4.2DR和PF得到的通量

圖2A為原始、DR和PF得到的觀測資料經(jīng)過超聲虛溫訂正和WPL訂正后的u*、Fc、LvE和H的對比．不同坐標旋轉(zhuǎn)方法對Fc、LvE和H的結果影響很??；坐標旋轉(zhuǎn)剔除了儀器傾斜引起的側(cè)向應力的影響，DR和PF分別使u*減小6%和3%．圖2B比較了兩種坐標旋轉(zhuǎn)方法得到的通量值，可以看出，兩種方法的H和LvE相差不多，PF得到的Fc較DR偏低，而PF得到的u*比DR高6.18%，這表明兩種坐標旋轉(zhuǎn)方法得到的三維風速有較大差異．

4.3分風區(qū)平面擬合

4.3.1分風區(qū)平面擬合

利用20天的風向做風向玫瑰圖(圖3)，可知SACOL主導風向為東南風和西北風．沿不同風向，地形有較大差異，為此分別選取風向范圍為100°～150°和280°～330°的觀測資料做平面擬合，稱為分風區(qū)平面擬合法(fetch planar fit，F(xiàn)PF)．

4.3.2三種坐標旋轉(zhuǎn)的比較

圖1　潛熱通量(a, c)和碳通量(b, d)的WPL訂正(2008年11月1—6日，3日除外)Fig.1　WPL correction of latent heat flux (a, c) and CO2 flux (b, d) (1—6 Nov 2008 except 3rd)

圖2　DR、PF得到的各通量對比(2008年11月1—21日，3日除外)Fig.2　Comparison of DR and PF fluxes during 1—21 Nov 2008 (except 3rd)

4.4質(zhì)量控制

4.4.1湍流平穩(wěn)性檢驗

根據(jù)Foken和Wichura(1996)提出的方法，將30 min長度的時間窗區(qū)分成M個(M=6)時長為5 min的子窗區(qū)，每個子窗區(qū)有N個(N=3000)數(shù)據(jù)點，計算每個子窗區(qū)的協(xié)方差，

(23)

圖3　SACOL 2008年11月1—21日風向頻數(shù)分布Fig.3　Frequency distribution of wind directions during 1—21 Nov 2008

6個子窗區(qū)的平均協(xié)方差為：

(24)

(25)定義湍流平穩(wěn)性檢驗指數(shù)：

(26)

表2　湍流平穩(wěn)性和ITC檢驗質(zhì)量劃分

圖4　三種坐標旋轉(zhuǎn)對比(a) 未訂正、PF訂正、FPF訂正垂直速度比較； (b) DR、(c)PF、(d)FPF對u*訂正結果.Fig.4　Comparison of three coordinate rotations(a) Vertical velocity without correction, with PF correction and PDF correction； (b) u* with DR correction； (c) u* with PF correction； (d) u* with FPF correction.

4.4.2總體湍流特征檢驗(ITC)

由于u*<0.1 m·s-1時EC對湍流通量觀測存在較大誤差(Zuo et al.，2009)，在考察無量綱標準差與大氣穩(wěn)定度關系時，剔除了u*<0.1 m·s-1時的資料．由圖6，無論是采用哪種坐標旋轉(zhuǎn)方法，垂直速度、溫度的無量綱標準差對穩(wěn)定度的擬合情況較好．張宏升等(2004)認為σw/u*和σT/T*主要受近地層局地特征尺度影響，尺度相對較小，而σu/u*并不完全取決于近地層局地特征尺度，在強不穩(wěn)定條件下，其影響尺度為混合層尺度．Wyngaard和Coté(1971)、Kaimal等(1982)的研究也表明σw/u*能較好地滿足相似性關系且與地形無關．這是由于垂直方向以小尺度高頻湍渦為主，尺度小的湍渦對地形變化的適應較快，地形起伏變化及下墊面物理特性差異對垂直風速的統(tǒng)計量影響較小，即黃土高原復雜下墊面不穩(wěn)定條件下的σw/u*與平坦下墊面接近且較好地滿足相似性關系(Moraes，2000；Al-Jiboori，2001)．水平方向的風速脈動主要由尺度較大的準水平湍流產(chǎn)生，一般為幾百米甚至更大，對地形的適應較慢，觀測到的氣流一般會“記憶”著上風方向的地形特點，產(chǎn)生較大的方差(趙鳴等，1991)．根據(jù)兩種坐標旋轉(zhuǎn)方法得到的相似性關系得出SOCAL的無量綱參數(shù)和穩(wěn)定度的關系，即c1，c2的值(表3)．

根據(jù)Foken和擬合出SACOL的c1，c2分別計算模擬出的歸一化標準差，得出ITC指數(shù)，進行質(zhì)量分類如圖7．采用DR進行坐標旋轉(zhuǎn)時，SACOL參數(shù)化方案得到u，w，T相對于Foken參數(shù)化方案得到的質(zhì)量明顯提高，即高質(zhì)量頻率分布增大，低質(zhì)量頻率分布減?。捎肞F進行坐標旋轉(zhuǎn)時，SACOL參數(shù)化方案得到的w和T質(zhì)量有明顯提高，但u的質(zhì)量無明顯變化．SACOL參數(shù)化方案對兩種坐標旋轉(zhuǎn)方法的適用性較好．不管是采用哪種坐標旋轉(zhuǎn)方法和參數(shù)化方案，w的質(zhì)量都要比u和T高，這說明湍流方差相似理論對w的適用性最好．

4.4.3總體質(zhì)量

將資料進行湍流平穩(wěn)性檢驗和ITC檢驗是為了篩選出高質(zhì)量數(shù)據(jù)，以用于進一步的研究．為了方便使用，Lee等(2005)提出了一套總體質(zhì)量的劃分方法，見表4．質(zhì)量等級為1～3的為高質(zhì)量數(shù)據(jù)，可用于基本研究，例如參數(shù)化方案的發(fā)展；質(zhì)量等級為4～6的為中等質(zhì)量數(shù)據(jù)，可用于長期觀測資料處理；質(zhì)量等級為7～9的為低質(zhì)量數(shù)據(jù)，應舍棄，必要時對缺失數(shù)據(jù)做插補．

圖6　不穩(wěn)定條件下T、u、w的無量綱標準差與大氣穩(wěn)定度z/L在(a，b，c)DR和(d，e，f)PF兩種坐標旋轉(zhuǎn)方法下的相似性關系Fig. 6　Standard deviations of temperature normalized by T* , and horizontal (u) , vertical velocity (w) normalized by u* as a function of stability (z/L) under unstable stratifications of DR (a,b,c) and PF (d,e,f)

表5給出了分別采用DR、PF和FPF進行坐標變換時，所得湍流通量總體質(zhì)量分布情況．66%～68%的感熱通量、62%～65%的潛熱通量和52%～54%的CO2通量為高質(zhì)量數(shù)據(jù)，坐標旋轉(zhuǎn)對這三者的總體質(zhì)量影響不大，PF相比DR得到的高質(zhì)量數(shù)據(jù)所占比例提高2%～3%，但對于摩擦速度，PF得到的高質(zhì)量u*只占45%，這與Zuo等(2009)的研究結果一致：復雜地形導致PF得到的u*質(zhì)量不高，需進一步研究以提高u*質(zhì)量；而DR得到的高質(zhì)量u*占到62%，相較于PF，DR能夠?qū)⒏哔|(zhì)量數(shù)據(jù)所占比例提高17%．為了與FPF得到的u*數(shù)據(jù)質(zhì)量進行比較，選取主導風向資料進行分析．由于使用DR和PF兩種坐標旋轉(zhuǎn)方法時，利用全部風向數(shù)據(jù)確定旋轉(zhuǎn)角度，即便有非主導風向的影響，但主導風向在決定旋轉(zhuǎn)角度時仍占主導地位，因此當剔除非主導風向時，PF和DR得到的u*質(zhì)量也得到提高．FPF得到的高質(zhì)量u*占61%，比PF有明顯提高，但DR得到的u*質(zhì)量仍然是最高的．Wilczak等(2001)指出：平均垂直速度的抽樣誤差可能導致DR的傾角估計誤差，增加縱向應力的隨機噪聲，從而導致應力計算的不確定性，PF使用的數(shù)據(jù)量增加，有效減小了隨機誤差．朱治林等(2004)研究表明：在復雜的地形條件下，PF方法是不合適的．綜合兼顧數(shù)據(jù)質(zhì)量和計算工作量，在復雜地形上處理湍流觀測資料的最優(yōu)坐標旋轉(zhuǎn)方法是DR．

表3　不同參數(shù)化方案各穩(wěn)定度對應的c1，c2

表4　總體湍流質(zhì)量分級

圖7　不同的參數(shù)化方案下兩種坐標旋轉(zhuǎn)的ITC檢驗質(zhì)量等級頻率分布(a, b)Foken參數(shù)化方案；(c, d) SACOL參數(shù)化方案. (a, c) DR; (b, d) PF.Fig.7　Frequency distributions of quality level for ITC test using different parameterization schemes(a, b) Foken parameterization scheme; (c, d) SACOL parameterization scheme.

質(zhì)量等級H(%)LvE(%)Fc(%)u*(%)主導風區(qū)u*(%)DRPFDRPFDRPFDRPFFPFPFDR1～34～67～966.69.124.468.18.223.762.222.615.265.220.913.952.527.719.954.126.519.561.819.918.345.026.828.261.120.818.156.219.724.274.513.511.9

5結論

利用SACOL的湍流通量觀測資料，比較分析了DR、PF、FPF坐標旋轉(zhuǎn)方法在黃土高原復雜地形上的適用性．應用超聲虛溫訂正、坐標旋轉(zhuǎn)、空氣密度脈動訂正以及平穩(wěn)性檢驗、總體湍流特征檢驗、總體質(zhì)量分級處理對原始觀測數(shù)據(jù)進行處理和質(zhì)量控制，有效地減少了因地形等原因造成的湍流觀測數(shù)據(jù)附加誤差，使湍流觀測數(shù)據(jù)質(zhì)量得到提高．

(1) DR和PF得到的感熱通量經(jīng)過超聲虛溫訂正分別減小了7.3%和5.9%，且穩(wěn)定度會影響超聲虛溫訂正效果．WPL訂正使?jié)摕嵬吭黾恿?.4%、CO2通量減小72.2%．

(2) 經(jīng)過DR和PF訂正，u*分別減小6%和3%．只考慮主導風向，DR對u*的訂正效果與風向無關，隨著湍流交換強度增大(u*>0.3 m·s-1)，PF和FPF逐漸呈現(xiàn)出SE風區(qū)的訂正值高于NW風區(qū)的訂正值，PF在SE風區(qū)使u*增大了9.23%，在NW風區(qū)使u*減小3.86%；而FPF在SE風區(qū)使u*增大了10.09%，在NW風區(qū)使u*減小1.18%．

(3) 在ITC檢驗中得到的SACOL參數(shù)化方案對兩種坐標旋轉(zhuǎn)方法的適用性都比較好．

(4) 采用DR得到的高質(zhì)量數(shù)據(jù)比例與采用PF相比，u*提高了17%，感熱通量、潛熱通量、CO2通量略降低2%～3%．PF和FPF兩種結果的差別主要體現(xiàn)在u*上，只考慮主導風向數(shù)據(jù)DR得到的u*質(zhì)量仍最好．綜合兼顧數(shù)據(jù)質(zhì)量和計算工作量，在復雜地形上處理湍流觀測資料的最優(yōu)坐標旋轉(zhuǎn)方法是DR．

致謝感謝蘭州大學半干旱氣候與環(huán)境觀測站(SACOL)為本文提供渦動相關等數(shù)據(jù)資料．感謝審稿專家的評議意見．

References

Al-Jiboori M H, Xu Y M, Qian Y F. 2001. Turbulence characteristics over complex terrain in west China.Boundary-LayerMeteorology, 101(1): 109-126.

Aubinet M, Vesala T, Papale D. 2011. Eddy Covariance: A Practical Guide to Measurement and Data Analysis. New York: Springer.

Baldocchi D, Falge L, Gu L H, et al. 2001. FLUXNET: A new tool to study the temporal and spatial variability of ecosystem-scale carbon dioxide, water vapor, and energy flux densities.Bull.Amer.Meteor.Soc., 82(11): 2415-2434.

Bao J. 2012. Study on land-atmosphere interaction in surface layer over semi-arid Loess Plateau of China [Ph. D. thesis] (in Chinese). Lanzhou: Lanzhou University.

Chen Z G, Bian L G, Lu L H, et al. 2008. Contrast and application of tilt correction methods for eddy covariance measurement.Meteor.Sci.Technol. (in Chinese), 36(3): 355-359.

Ding Y H. 1997. On some aspects of estimates of the surface fluxes.QuarterlyJournalofAppliedMeteorology(in Chinese), 8(S): 29-35.

Foken T, Wichura B. 1996. Tools for quality assessment of surface-based flux measurements.AgriculturalandForestMeteorology, 78(1-2): 83-105.Foken T. 2008. The energy balance closure problem: An overview.Ecol.Appl., 18(6): 1351-1367.

G?ckede M, Foken T, Aubinet M, et al. 2008. Quality control of CarboEurope flux data-Part 1: Coupling footprint analyses with flux data quality assessment to evaluate sites in forest ecosystems.Biogeosciences, 5(2): 433-450.

Huang J P, Zhang W, Zuo J Q, et al. 2008. An overview of the semi-arid climate and environment research observatory over the loess plateau.Adv.Atmos.Sci., 25(6): 906-921.

Jiang H M, Liu S H, Liu H P. 2012. A study on energy budget characteristics over a heterogeneously irrigated cotton field.ChineseJ.Geophys. (in Chinese), 55(2): 428-440, doi: 10.6038/j.issn.0001-5733.2012.02.007.

Jiang H M, Liu S H, Zhang L, et al. 2013. A study of turbulent heat flux corrections and energy balance closure problem on the surface layer in EBEX-2000.ActaScientiarumNaturaliumUniversitatisPekinensis(in Chinese), 49(3): 443-451.

Jiang H M, Liu S H, Liu H P. 2013. Spectra and cospectra of turbulence in an internal boundary layer over a heterogeneously irrigated cotton field.ActaMeteor.Sin., 27(2): 233-248.

Kaimal J C, Finnigan J J. 1994. Atmospheric Boundary Layer Flows: Their Structure and Measurement. Oxford: Oxford University Press.

Kaimal J C, Eversole R A, Lenschow D H, et al. 1982. Spectral characteristics of the convective boundary layer over uneven terrain.J.Atmos.Sci., 39(5): 1098-1114. Lee X H, Massman W, Law B. 2005. Handbook of Micrometeorology: A Guide for Surface Flux Measurement and Analysis. Netherlands: Springer.Liang J N, Zhang L, Tian P F, et al. 2014. Impact of low-level jets on turbulent in nocturnal boundary layer over complex terrain of the Loess Plateau.ChineseJ.Geophys. (in Chinese), 57(5): 1387-1398, doi: 10.6038/cjg20140504.

Liu H P, Peters G, Foken T. 2001. New equations for sonic temperature variance and buoyancy heat flux with an omnidirectional sonic anemometer.Boundary-LayerMeteorology, 100(3): 459-468. Liu S H, Li J, Liu H P, et al. 2005a. Characteristics of macroturbulence variables in EBEX-2000.Chin.J.Atmos.Sci. (in Chinese), 29(4): 503-509.

Liu S H, Liu H P, Li J, et al. 2005b. Characteristics of turbulence dissipation rates, characteristic length scales and structure parameters in EBEX-2000.Chin.J.Atmos.Sci. (in Chinese), 29(3): 475-481.

Liu S H, Hu Y, Hu F, et al. 2005c. Numerical simulation of land-atmosphere interaction and oasis effect over oasis-desert.ChineseJ.Geophys. (in Chinese), 48(5): 1019-1027, doi: 10.3321/j.issn:0001-5733.2005.05.007.

Liu S H, Mao Y H, Hu F, et al. 2009. A comparative study of computing methods of turbulent fluxes on different underling surfaces.ChineseJ.Geophys. (in Chinese), 52(3): 616-629.

Mauder M, Liebethal C, G?ckede M, et al. 2006. Processing and quality control of flux data during LITFASS-2003.Boundary-LayerMeteorology, 122(1): 67-88.

Moraes O L L. 2000. Turbulence characteristics in the surface boundary layer over the South American Pampa.Boundary-LayerMeteorology, 96(3): 317-335.

Schotanus P, Nieuwstadt F T M, De Bruin H A R. 1983. Temperature measurement with a sonic anemometer and its application to heat and moisture fluxes.Boundary-LayerMeteorology, 26(1): 81-93.

Swinbank W C. 1951. The measurement of vertical transfer of heat and water vapor by eddies in the lower atmosphere.J.Meteor., 8(3): 135-145.

Vickers D, Mahrt L. 1997. Quality control and flux sampling problems for tower and aircraft data.J.Atmos.Oceanic.Technol., 14(3): 512-526.Wang J M, Wang W Z, Ao Y H, et al. 2007. Turbulence flux measurements under complicated conditions.Adv.EarthSci. (in Chinese), 22(8): 791-797, doi: 10.3321/j.issn:1001-8166.2007.08.004.

Wang S Y, Zhang Y, Lü S H, et al. 2009. The preliminary study on turbulence data quality control of Jinta oasis.PlateauMeteorology(in Chinese), 28(6): 1260-1273.

Webb E K, Pearman G I, Leuning R. 1980. Correction of flux measurements for density effects due to heat and water vapour transfer.Quart.J.Roy.Meteor.Soc., 106(447): 85-100.

Wilczak J M, Oncley S P, Stage S A. 2001. Sonic anemometer tilt correction algorithms.Boundary-LayerMeteorology, 99(1): 127-150.

Wyngaard J C, Coté O R. 1971. The budgets of turbulent kinetic energy and temperature variance in the atmospheric surface layer.J.Atmos.Sci., 28(2): 190-201.

Xu Z W, Liu S M, Gong L J, et al. 2008. A study on the data processing and quality assessment of the eddy covariance system.Adv.EarthSci. (in Chinese), 23(4): 357-370, doi: 10.3321/j.issn:1001-8166.2008.04.005.

Yang W Z, Shao M A. 2000. The Research of the Loess Plateau Soil Moisture (in Chinese). Beijing: Science Press.

Zhang H S, Li F Y, Chen J Y. 2004. Statistical characteristics of atmospheric turbulence in different underlying surface conditions.PlateauMeteorology(in Chinese), 23(5): 598-604.Zhao M, Miao M Q, Wang Y C. 1991. Boundary-Layer Meteorology (in Chinese). Beijing: China Meteorological Press.Zhu Z L, Sun X M, Yuan G F, et al. 2004. Correcting method of eddy covariance fluxes over non-flat surfaces and its application in China FLUX.Sci.ChinaSer.D:EarthSci., 48(S1): 42-50. Zuo J Q, Huang J P, Wang J M, et al. 2009. Surface turbulent flux measurements over the loess plateau for a semi-arid climate change study.Adv.Atmos.Sci., 26(4): 679-691.

附中文參考文獻

鮑婧. 2012. 黃土高原半干旱區(qū)近地層陸氣相互作用研究[博士論文]. 蘭州: 蘭州大學. 諶志剛, 卞林根, 陸龍驊等. 2008. 渦度相關儀傾斜訂正方法的比較及應用. 氣象科技, 36(3): 355-359.

丁一匯. 1997. 地表通量的計算問題. 應用氣象學報, 8(S): 29-35.

姜海梅, 劉樹華, 劉和平. 2012. 非均勻灌溉棉田能量平衡特征研究. 地球物理學報, 55(2): 428-440, doi: 10.6038/j.issn.0001-5733.2012.02.007.

姜海梅, 劉樹華, 張磊等. 2013. EBEX-2000湍流熱通量訂正和地表能量平衡閉合問題研究. 北京大學學報(自然科學版), 49(3): 443-451.

梁捷寧, 張鐳, 田鵬飛等. 2014. 黃土高原復雜地形上邊界層低空急流對近地層湍流的影響. 地球物理學報, 57(5): 1387-1398, doi: 10.6038/cjg20140504.

劉樹華, 李潔, 劉和平等. 2005a. 在EBEX-2000實驗資料中的湍流宏觀量特征. 大氣科學, 29(4): 503-509.

劉樹華, 劉和平, 李潔等. 2005b. 在EBEX-2000實驗資料中湍流耗散率、長度尺度和結構參數(shù)特征. 大氣科學, 29(3): 475-481.

劉樹華, 胡予, 胡非等. 2005c. 沙漠-綠洲-陸-氣相互作用和綠洲效應的數(shù)值模擬. 地球物理學報, 48(5): 1019-1027, doi: 10.3321/j.issn:0001-5733.2005.05.007.

劉樹華, 茅宇豪, 胡非等. 2009. 不同下墊面湍流通量計算方法的

比較研究. 地球物理學報, 52(3): 616-629.

王介民, 王維真, 奧銀煥等. 2007. 復雜條件下湍流通量的觀測與分析. 地球科學進展, 22(8): 791-797, doi: 10.3321/j.issn:1001-8166.2007.08.004.

王少影, 張宇, 呂世華等. 2009. 金塔綠洲湍流資料的質(zhì)量控制研究. 高原氣象, 28(6): 1260-1273.

徐自為, 劉紹民, 宮麗娟等. 2008. 渦動相關儀觀測數(shù)據(jù)的處理與質(zhì)量評價研究. 地球科學進展, 23(4): 357-370, doi: 10.3321/j.issn:1001-8166.2008.04.005. 楊文治, 邵明安. 2000. 黃土高原土壤水分研究. 北京: 科學出版社. 張宏升, 李富余, 陳家宜. 2004. 不同下墊面湍流統(tǒng)計特征研究. 高原氣象, 23(5): 598-604. 趙鳴, 苗曼倩, 王彥昌. 1991. 邊界層氣象學教程. 北京: 氣象出版社. 朱治林, 孫曉敏, 袁國富等. 2004. 非平坦下墊面渦度相關通量的校正方法及其在ChinaFLUX中的應用. 中國科學D輯: 地球科學, 34(S2): 37-45.

(本文編輯何燕)

doi:10.6038/cjg20160604 中圖分類號P401, P404

收稿日期2015-04-28，2016-03-03收修定稿

基金項目國家自然科學基金項目(41475008)和國家重大科學研究計劃項目(2012CB955302)資助.

作者簡介苑廣輝，女，1991年生，主要從事大氣邊界層和大氣湍流研究.E-mail：yuangh09@lzu.edu.cn

*通訊作者張鐳，教授，博士生導師，主要從事大氣物理與大氣環(huán)境研究.E-mail：zhanglei@lzu.edu.cn

Method of acquiring high-quality surface turbulent fluxes over the Loess Plateau

YUAN Guang-Hui, ZHANG Lei*, LI Yao, LIANG Jie-Ning

KeyLaboratoryforSemi-AridClimateChangeoftheMinistryofEducation,CollegeofAtmosphericSciences,LanzhouUniversity,Lanzhou730000,China

AbstractThe eddy covariance (EC) technique measuring the turbulent exchanges of heat, moisture, CO2 and momentum between surface and atmosphere has been used widely. However, the use of EC is based on some assumptions which do not exist in practice. So the results are not accurate without necessary corrections. Using data collected at the Semi-Arid Climate and Environment Observatory of Lanzhou University (SACOL), this work studies a method of acquiring high quality surface turbulent flux. In this approach, eddy covariance data processing includes eliminating despikes, coordinate rotations, sonic temperature correction, and WPL correction (correction for density fluctuations). Double rotation and planar fit are used to do coordinate rotations. The fluctuations of sonic temperature include the effect of humidity on the speed of sound, and should be converted into actual temperature. WPL correction was required to latent heat flux and CO2 flux for the density effects due to heat and water vapor transfer on turbulent flux measurements. Quality control is performed by the stationarity test, integral turbulent characteristics (ITC) test and overall quality is controlled to permit selecting high quality data. The results show that sonic temperature correction decreases sensible heat flux by about 7.3% when using DR, but 5.9% when using PF. The stability can influence the results of sonic temperature correction. WPL increases latent heat flux by 7.4% and decreases CO2 flux by 72.2%, respectively. Three coordinate rotations have great influence on momentum but little on scalar fluxes. The values of u* obtained from DR and PF are decreased by 6% and 3%, respectively. Only considering the dominant wind direction, the wind direction has no relationship with the correction of u* by DR which reduces u* by 5%. As the turbulent exchange strengthened (u*> 0.3 m·s-1), PF and FPF gradually exhibit that the corrected values in southeast are bigger than the northwest area, PF in the southeast wind area increases u* by 9.23%, yet decreases 3.86% in the northwest area and FPF in southeast wind area increases u* by 10.09%, yet decreases 1.18% in the northwest area. The differences of the steady state tests between DR and PF are mainly in u*. A parameterization scheme of DR and PF for SACOL provided in the ITC test works well. The overall quality shows that about 45%～62% of the total data is higher for u*, 66%～68% for sensible heat flux, 62%～65% for latent flux and 52%～54% for CO2 flux. The proportion of the high quality of u* obtained by DR is 17% higher than PF, while that of the latter three kinds of fluxes obtained by PF is 2%～3% higher than DR. The difference between PF and FPF is mainly in u*. Comparing the three coordinate rotations in the dominant wind direction, DR still obtains the best quality of u*. The use of DR is recommended in the complex terrain for reducing calculation and improving the data quality.

KeywordsLoess Plateau; Eddy covariance; Flux correction; Quality control; Fetch planar fit

苑廣輝， 張鐳， 李遙等. 2016. 黃土高原復雜地形上高質(zhì)量湍流通量數(shù)據(jù)獲取方法. 地球物理學報，59(6):1971-1982,doi:10.6038/cjg20160604.

Yuan G H, Zhang L, Li Y,et al. 2016. Method of acquiring high-quality surface turbulent fluxes over the Loess Plateau.ChineseJ.Geophys. (in Chinese),59(6):1971-1982,doi:10.6038/cjg20160604.