黄土高原复杂地形上高质量湍流通量数据获取方法

苑广辉， 张镭， 李遥， 梁捷宁

半干旱气候变化教育部重点实验室, 兰州大学大气科学学院, 兰州　730000



黄土高原复杂地形上高质量湍流通量数据获取方法

苑广辉， 张镭*， 李遥， 梁捷宁

半干旱气候变化教育部重点实验室, 兰州大学大气科学学院, 兰州730000

摘要利用兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University，SACOL)湍流观测资料，分析了二次坐标旋转(double rotation，DR)、平面拟合(planar fit，PF)和分风区平面拟合(fetch planar fit，FPF)在复杂地形上的适用性，总结出一套适用于SACOL的总体湍流特征参数化方案．经过超声虚温订正、坐标旋转、空气密度脉动订正以及平稳性检验、总体湍流特征检验、总体质量分级处理，摩擦速度(u*)、感热通量、潜热通量、CO2通量高质量数据所占比例分别为45%～62%、66%～68%、62%～65%、52%～54%．采用DR得到的高质量数据比例与采用PF相比，u*提高了17%，后三种通量略降低2%～3%．PF和FPF两种结果的差别主要体现在u*上，只考虑主导风向数据DR得到的u*质量仍最好．综合兼顾数据质量和计算工作量，在复杂地形上处理湍流观测资料的最优坐标旋转方法是DR．

关键词黄土高原； 涡动相关； 通量订正； 质量控制； 分风区平面拟合

1引言

在近地层，湍流是大气运动的主要形式，引起各物理属性在地-气之间的交换输送，地表热量、水汽、CO2等物质的传输过程受到广泛关注(Jiang et al.，2013；刘树华等，2005a, 2005b；徐自为等，2008)，截至2014年4月全球通量观测站网络(FLUXNET)已有683个站点(http:∥fluxnet.ornl.gov/)．涡动相关法(eddy covariance，EC)可以用来长期连续观测CO2、水汽和热通量，广泛应用于地气间通量交换的测量，已成为FLUXNET测量通量的主要技术手段(Baldocchi et al.，2001)．EC建立在一定的假设之上，如湍流平稳、下垫面水平均匀、湍流充分发展、存在常通量层等(Foken and Wichura，1996)，而在黄土高原复杂下垫面上，这些假设条件通常难以得到满足，地形起伏对EC测量通量有重要影响(王介民等，2007)．如何准确地计算近地层通量一直受到地球科学界的重视(丁一汇，1997；姜海梅等，2013；刘树华等，2005c, 2009)，已有不少学者对湍流资料处理和质量控制进行研究：Mauder等(2006)通过对LITFASS-2003实验中14个测点的数据进行总体质量检验，发现80%的潜热通量为高质量数据；Göckede等(2008)将拉格朗日随机足迹模型与质量评价相结合应用到CarboEurope的25个森林测站来研究通量测量的空间代表性、仪器以及坐标旋转方法对湍流通量的影响；王少影等(2009)讨论了在绿洲和戈壁测站不同坐标旋转方法对摩擦速度、感热、潜热通量和湍流平稳性检验以及总体湍流特征检验的影响；朱治林等(2004)讨论了在非平坦、非均匀下垫面上，仪器安装不垂直在EC计算通量时可能产生的误差及其校正方法，对比分析了不同坐标旋转方法的校正效果和适用条件；在平坦均匀下垫面上，谌志刚等(2008)的分析发现，平面拟合方法优于三次坐标旋转方法；姜海梅等(2012)运用EBEX-2000实验的湍流、净辐射和土壤观测资料，运用涡动相关方法分析了非均匀灌溉引起的热内边界层发展条件下近地层感热、潜热通量特征，并对有无灌溉两种条件下的能量闭合度进行了对比分析．黄土高原半干旱地区面积广阔，其陆气相互作用不仅对西北干旱气候的形成以及东亚季风环流有着不可忽视的影响，同时对全球气候和大气环流的变化也可能产生比较重要的作用(杨文治和邵明安，2000)．这一区域的陆气相互作用问题已经成为我国迫切需要研究的重要基础性科学问题之一．在黄土高原复杂下垫面上，地形起伏，沟壑纵横，植被稀疏，湍流更为复杂，而关于湍流资料处理和质量控制的研究较少，获取高质量通量数据对深入认识该地区地气交换特征非常必要．

本文主要讨论在黄土高原复杂地形上各种湍流资料处理方法对EC观测结果的影响，应用超声虚温订正、坐标旋转、空气密度脉动订正以及平稳性检验、总体湍流特征检验、总体质量分级处理方法，对兰州大学半干旱气候与环境观测站(Semi-Arid Climate and Environment Observatory of Lanzhou University，SACOL)2008年11月1日至21日的10 Hz原始湍流观测数据进行了质量控制，有效地减少了因地形等原因造成的湍流观测数据附加误差，使其湍流观测数据质量得到提高．根据SACOL风向特点使用分风区平面拟合方法，着重关注二次坐标旋转、平面拟合及分风区平面拟合的适用性，进而分析SACOL所代表的黄土高原复杂地形上近地层湍流特征．

2观测资料

SACOL建立于2005年，位于海拔高度为1965.8 m的萃英山顶(35.95°N，104.13°E)，距兰州市中心约48 km，位于中国黄土高原半干旱区，下垫面为典型的黄土高原残塬地貌，沟壑纵横，土壤为第四纪黄土风蚀形成的灰钙土．拥有国际先进的观测仪器，是继中国科学院吉林通榆站之后，第二个由我国自主建设的半干旱区长期观测站．已被批准加入国际协同观测计划项目(The Coordinated Enhan-ced Observing Period)，并作为此计划的全球协同加强观测站之一(Huang et al.，2008；梁捷宁等，2014)．

通量观测场地较为平坦，东西方向约200 m，南北方向约1000 m．地表为长芒草、冷蒿、赖草等短小植被覆盖，冬季地表植被高约0.10 m(Huang et al.，2008；梁捷宁等，2014)．

涡动相关系统分别利用三维超声风温仪(CSAT3，Campbell)测量u，v，w三维风速和超声虚温Ts，开路红外CO2/H2O气体分析仪(LI7500，LI-COR)测量CO2和H2O密度，感应器采样频率为10 Hz，观测高度为3 m，数据采集器(CR5000，Campbell)．

2008年夏季对观测仪器进行校准，此后数据可信度较高．冬季下垫面植被低矮稀疏，对湍流通量影响较小，故选取2008年11月1日至21日的10 Hz原始观测数据，资料完好率高，由于11月3日数据中有大量异常值出现，所以剔除11月3日全天数据．

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坐标旋转

坐标旋转方法主要有二次坐标旋转(double rotation，DR)、三次坐标旋转(triple rotation，TR)和平面拟合(planar fit，PF)，三次坐标旋转在计算应力时误差较大，已不推荐使用(Kaimal and Finnigan，1994)．

(1)二次坐标旋转

(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)是根据一个较长的时段，本文选取2008年1日至21日(3日除外)20天的数据，确定一个与地面平行的平均风场，并将各时次的u，v，w旋转到该平面上．利用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订正使潜热通量总体增加了7.4%，白天增加5.6%，夜间降低3.4%，这与感热输送方向的昼夜差异有关．白天，向上输送的感热加热空气，使观测高度处水汽密度减小，导致EC对水汽通量的观测值偏低，夜间则相反．由图1b, 1d可看出，向下的CO2通量在正午前后达到最高，向上的最高值出现在0时左右．白天，植物光合作用吸收大气中的CO2，形成向下的CO2通量，并于正午达到最强；夜间，植物呼吸作用释放出CO2，使得CO2通量向上传递，形成一峰一谷的日变化趋势．总体上，WPL订正使CO2通量降低了72.2%，白天向下的CO2通量降低64.7%，夜间向上的CO2通量降低23.4%．

4.2DR和PF得到的通量

图2A为原始、DR和PF得到的观测资料经过超声虚温订正和WPL订正后的u*、Fc、LvE和H的对比．不同坐标旋转方法对Fc、LvE和H的结果影响很小；坐标旋转剔除了仪器倾斜引起的侧向应力的影响，DR和PF分别使u*减小6%和3%．图2B比较了两种坐标旋转方法得到的通量值，可以看出，两种方法的H和LvE相差不多，PF得到的Fc较DR偏低，而PF得到的u*比DR高6.18%，这表明两种坐标旋转方法得到的三维风速有较大差异．

4.3分风区平面拟合

4.3.1分风区平面拟合

利用20天的风向做风向玫瑰图(图3)，可知SACOL主导风向为东南风和西北风．沿不同风向，地形有较大差异，为此分别选取风向范围为100°～150°和280°～330°的观测资料做平面拟合，称为分风区平面拟合法(fetch planar fit，FPF)．

4.3.2三种坐标旋转的比较

图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质量控制

4.4.1湍流平稳性检验

根据Foken和Wichura(1996)提出的方法，将30 min长度的时间窗区分成M个(M=6)时长为5 min的子窗区，每个子窗区有N个(N=3000)数据点，计算每个子窗区的协方差，

(23)

图3　SACOL 2008年11月1—21日风向频数分布Fig.3　Frequency distribution of wind directions during 1—21 Nov 2008

6个子窗区的平均协方差为：

(24)

(25)定义湍流平稳性检验指数：

(26)

表2　湍流平稳性和ITC检验质量划分

图4　三种坐标旋转对比(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)，在考察无量纲标准差与大气稳定度关系时，剔除了u*<0.1 m·s-1时的资料．由图6，无论是采用哪种坐标旋转方法，垂直速度、温度的无量纲标准差对稳定度的拟合情况较好．张宏升等(2004)认为σw/u*和σT/T*主要受近地层局地特征尺度影响，尺度相对较小，而σu/u*并不完全取决于近地层局地特征尺度，在强不稳定条件下，其影响尺度为混合层尺度．Wyngaard和Coté(1971)、Kaimal等(1982)的研究也表明σw/u*能较好地满足相似性关系且与地形无关．这是由于垂直方向以小尺度高频湍涡为主，尺度小的湍涡对地形变化的适应较快，地形起伏变化及下垫面物理特性差异对垂直风速的统计量影响较小，即黄土高原复杂下垫面不稳定条件下的σw/u*与平坦下垫面接近且较好地满足相似性关系(Moraes，2000；Al-Jiboori，2001)．水平方向的风速脉动主要由尺度较大的准水平湍流产生，一般为几百米甚至更大，对地形的适应较慢，观测到的气流一般会“记忆”着上风方向的地形特点，产生较大的方差(赵鸣等，1991)．根据两种坐标旋转方法得到的相似性关系得出SOCAL的无量纲参数和稳定度的关系，即c1，c2的值(表3)．

根据Foken和拟合出SACOL的c1，c2分别计算模拟出的归一化标准差，得出ITC指数，进行质量分类如图7．采用DR进行坐标旋转时，SACOL参数化方案得到u，w，T相对于Foken参数化方案得到的质量明显提高，即高质量频率分布增大，低质量频率分布减小．采用PF进行坐标旋转时，SACOL参数化方案得到的w和T质量有明显提高，但u的质量无明显变化．SACOL参数化方案对两种坐标旋转方法的适用性较好．不管是采用哪种坐标旋转方法和参数化方案，w的质量都要比u和T高，这说明湍流方差相似理论对w的适用性最好．

4.4.3总体质量

将资料进行湍流平稳性检验和ITC检验是为了筛选出高质量数据，以用于进一步的研究．为了方便使用，Lee等(2005)提出了一套总体质量的划分方法，见表4．质量等级为1～3的为高质量数据，可用于基本研究，例如参数化方案的发展；质量等级为4～6的为中等质量数据，可用于长期观测资料处理；质量等级为7～9的为低质量数据，应舍弃，必要时对缺失数据做插补．

图6　不稳定条件下T、u、w的无量纲标准差与大气稳定度z/L在(a，b，c)DR和(d，e，f)PF两种坐标旋转方法下的相似性关系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进行坐标变换时，所得湍流通量总体质量分布情况．66%～68%的感热通量、62%～65%的潜热通量和52%～54%的CO2通量为高质量数据，坐标旋转对这三者的总体质量影响不大，PF相比DR得到的高质量数据所占比例提高2%～3%，但对于摩擦速度，PF得到的高质量u*只占45%，这与Zuo等(2009)的研究结果一致：复杂地形导致PF得到的u*质量不高，需进一步研究以提高u*质量；而DR得到的高质量u*占到62%，相较于PF，DR能够将高质量数据所占比例提高17%．为了与FPF得到的u*数据质量进行比较，选取主导风向资料进行分析．由于使用DR和PF两种坐标旋转方法时，利用全部风向数据确定旋转角度，即便有非主导风向的影响，但主导风向在决定旋转角度时仍占主导地位，因此当剔除非主导风向时，PF和DR得到的u*质量也得到提高．FPF得到的高质量u*占61%，比PF有明显提高，但DR得到的u*质量仍然是最高的．Wilczak等(2001)指出：平均垂直速度的抽样误差可能导致DR的倾角估计误差，增加纵向应力的随机噪声，从而导致应力计算的不确定性，PF使用的数据量增加，有效减小了随机误差．朱治林等(2004)研究表明：在复杂的地形条件下，PF方法是不合适的．综合兼顾数据质量和计算工作量，在复杂地形上处理湍流观测资料的最优坐标旋转方法是DR．

表3　不同参数化方案各稳定度对应的c1，c2

表4　总体湍流质量分级

图7　不同的参数化方案下两种坐标旋转的ITC检验质量等级频率分布(a, b)Foken参数化方案；(c, d) SACOL参数化方案. (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.

质量等级H(%)LvE(%)Fc(%)u*(%)主导风区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坐标旋转方法在黄土高原复杂地形上的适用性．应用超声虚温订正、坐标旋转、空气密度脉动订正以及平稳性检验、总体湍流特征检验、总体质量分级处理对原始观测数据进行处理和质量控制，有效地减少了因地形等原因造成的湍流观测数据附加误差，使湍流观测数据质量得到提高．

(1) DR和PF得到的感热通量经过超声虚温订正分别减小了7.3%和5.9%，且稳定度会影响超声虚温订正效果．WPL订正使潜热通量增加了7.4%、CO2通量减小72.2%．

(2) 经过DR和PF订正，u*分别减小6%和3%．只考虑主导风向，DR对u*的订正效果与风向无关，随着湍流交换强度增大(u*>0.3 m·s-1)，PF和FPF逐渐呈现出SE风区的订正值高于NW风区的订正值，PF在SE风区使u*增大了9.23%，在NW风区使u*减小3.86%；而FPF在SE风区使u*增大了10.09%，在NW风区使u*减小1.18%．

(3) 在ITC检验中得到的SACOL参数化方案对两种坐标旋转方法的适用性都比较好．

(4) 采用DR得到的高质量数据比例与采用PF相比，u*提高了17%，感热通量、潜热通量、CO2通量略降低2%～3%．PF和FPF两种结果的差别主要体现在u*上，只考虑主导风向数据DR得到的u*质量仍最好．综合兼顾数据质量和计算工作量，在复杂地形上处理湍流观测资料的最优坐标旋转方法是DR．

致谢感谢兰州大学半干旱气候与环境观测站(SACOL)为本文提供涡动相关等数据资料．感谢审稿专家的评议意见．

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. 黄土高原半干旱区近地层陆气相互作用研究[博士论文]. 兰州: 兰州大学. 谌志刚, 卞林根, 陆龙骅等. 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实验资料中湍流耗散率、长度尺度和结构参数特征. 大气科学, 29(3): 475-481.

刘树华, 胡予, 胡非等. 2005c. 沙漠-绿洲-陆-气相互作用和绿洲效应的数值模拟. 地球物理学报, 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. 金塔绿洲湍流资料的质量控制研究. 高原气象, 28(6): 1260-1273.

徐自为, 刘绍民, 宫丽娟等. 2008. 涡动相关仪观测数据的处理与质量评价研究. 地球科学进展, 23(4): 357-370, doi: 10.3321/j.issn:1001-8166.2008.04.005. 杨文治, 邵明安. 2000. 黄土高原土壤水分研究. 北京: 科学出版社. 张宏升, 李富余, 陈家宜. 2004. 不同下垫面湍流统计特征研究. 高原气象, 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

*通讯作者张镭，教授，博士生导师，主要从事大气物理与大气环境研究.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. 黄土高原复杂地形上高质量湍流通量数据获取方法. 地球物理学报，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.