涡动相关观测与数据处理基础知识系列之一:通量塔的选址与建塔的基本原则
近年来,采用涡动相关(eddy-covariance,EC)方法测量温室气体通量的站点数量在迅速增加,但是要在科学目的、工程标准、安装运行成本和实用性之间做出平衡,寻找到最佳的解决方法,仍是一个具有挑战的工作。从观测结果准确性和精确度来说,选址、建塔等站点设计的环节是重中之重。1、位置选择站点选址的基本原则是,该站点能够尽量观测到全部的研究对象,这涉及到两个问题,一个是方向,一个是架设高度。首先是确定观测区域近几年的主风向,可以参考近几年的气象数据。由于中国大部分地区是季风气候,一般在春夏和秋冬会有两个主风向,这时候要考虑通量仪器的架设方向,实验观测的主要周期等。如果仪器架设方向可以随主风向的改变方便调整,或者实验周期是明确区分了春夏或者秋冬,那么在选址时可以选在观测对象的下风向,这样可以尽可能多的观测到目标对象;如果不能改变通量仪器的架设方向,且是长期定位观测,那尽量将观测地点选址在观测对象的中央位置,或者沿主风向的中点位置,这样可以尽可能的在不改变仪器方向和位置的前提下,观测到尽可能多的研究对象。确定架设高度要满足通量仪器的基本观测条件, 即满足湍流运动的充分交换。一般的架设高度是下垫面冠层高度的1.5到2倍(具体确定观测高度的经验法则见图 1);在相对平坦和均匀的下垫面条件下,观测距离大约是观测有效高度的100倍(风浪区原理),具体范围需要根据footprint源区计算,随着湍流运动强度和下垫面情况会有所改变。图 1 确定观测高度的经验法则通量源区代表性分析(Footprint分析)是检验一个通量站质量的重要手段,可以用来进行实验方案的设计指导,观测数据的质量控制,以及通过特定传感器的源区分布和来自感兴趣下垫面(植被)的通量贡献,从而对观测结果进行分析解释。图 2 Footprint分析2、下垫面的影响2.1植被类型涡动相关法测量温室气体通量要求仪器安装在常通量层内,而常通量层假设要求稳态大气、下垫面与仪器之间没有任何源或者汇、足够长的风浪区和水平均匀的下垫面等基本条件。在涡动相关传感器能监测到的“源区域”内植被类型均匀一致的情况下,其观测到的通量结果是比较有意义的,可以用来解释生态系统的温室气体收支情况。但当涡动相关传感器的“源区域”覆盖到不同植被类型时,情况就会变得复杂起来。一个极端的例子是:某站点周围具有两种不同的森林植被类型,每天周期性地,白天,风从一种植被类型吹向另一种;夜间,则正好相反。那么,该站点观测得到的通量资料的日平均值将毫无意义。这种极端的情况虽然极少出现,但许多站点都会有微妙的风向变化,在数据分析时需要做仔细考虑。此外,光、土壤湿度、土壤结构、叶面积以及物种种类组成的空间异质性会导致温室气体源/汇强度的水平梯度。而其植被类型的变化也会造成表面粗糙度的变化,当风通过不同粗糙度或者不同源/汇强度表面的区域时,就会产生非常明显的平流效应(Raupach & Finnigan, 1997 Baldocchi et al., 2000)。图 3 不同下垫面的地表粗糙度(参考 于贵瑞&孙晓敏,2006)地表植被类型的突然变化会导致气流的变化,如气流在从高大森林向低矮草地移动时,会在森林边缘形成回流区(如图 4所示),导致近地面和上方气流方向不一致,其水平长度尺度(距离)等于冠层高度的2-5倍(Detto et al., 2008)。图 4森林边缘附近湍流结构的概念模型(参考Detto et al., 2008)2.2冠层高度通量足迹Footprint描述了EC系统能够观测到的“源区域”,提供了每个表面元素对测量的垂直通量的相对贡献。Footprint取决于观测高度、表面粗糙度和大气稳定度等。如图 5所示,通常来说,传感器的观测高度越高,就越能观测到更远、更广的区域(Horst & Weil, 1994),也便于捕捉植物冠层上方混合良好的边界层中的通量交换。但是观测高度也不是越高越好,在大气层结稳定的条件下(如夜间),过高的观测高度可能会使观测到的“源区域”超出感兴趣的研究区域。因此应该预先计算并确保来自感兴趣区域的通量贡献至少为90%(Gö ckede et al., 2004),在稳定条件下至少50%的时间以确保适当的数据覆盖不同的风向和不同的天气条件。图 5观测高度与通量足迹基于Munger(2012)等确定塔/测量高度(hm)的原则(如图 1),可能存在准确测量实际观测高度和冠层高度的困难,需要考虑后期调整高度的可能性。观测高度必须用三维超声风速计测量路径的中心来确定,其值取决于感兴趣的生态系统的冠层高度(hc),冠层高度值不需要特别准确:采用主要冠层的平均预期高度是合理的。对于冠层高度在生长季节中快速变化的农田、草地和种植园以及同样具有快速变化特性的冰雪下垫面,塔架设计必须考虑允许通过改变塔架高度(例如伸缩式塔架设计)或通过移动传感器来改变测量高度。随着时间的推移为了确保相同的通量观测源区,可以考虑改变测量高度,遵循的原则是测量高度与冠层高度的0.76倍之间的差值保持在一个确定数值的±10%左右。但这种调整的频率不用特别频繁,最多在植被生长期或在积雪季节每隔一周进行。假设在植被生长期开始时的裸土,其测量高度为2 m,在冠层高度达到1.2 米前,不需要改变测量高度;在植被达到1.2米后(例如增加约0.5-0.8米)开始提高测量高度,然后保持测量高度与冠层高度的0.76倍之间的差值保持在一个确定数值。改变表面高度(由于生长和积雪)以及改变测量高度必须准确记录,因为这必须在后期数据处理中考虑。2.3地形影响EC法测量通量假设了地形水平,这样可以保证地形的坐标系和传感器坐标系方向一致,避免平流、泄流效应的影响。图 6复杂地形对EC观测的影响在复杂的地形条件下,风吹过小山时会引起气流的辐合或辐散运动,产生平流效应(Kaimail & Finnigan, 1994)。存在有局地风场影响的站点,在夜间大气稳定,垂直湍流输送和大气混合作用较弱,CO2的水平和垂直平流效应的影响是很重要的(于贵瑞&孙晓敏,2006)。Mordukhovish & Tsvang(1966)的研究表明,斜坡地形能导致水平异质和通量的辐散。对于设在地势较高的观测塔,在夜间对流比较弱时,通常会因CO2沿斜坡泄流而造成大气传输的通量低估,最后导致生态系统净生产力的估算偏高;对于在地势较低沟谷中的观测塔,其问题更加复杂,如果外部的大气平流/泄流通过观测界面进入生态系统,会高估光合作用吸收CO2的能力;如果外部的大气平流/泄流不能通过观测界面,而是从观测界面下部直接进入生态系统,则会在生态系统中暂时储存,最终输出生态系统,造成对呼吸作用的高估。在大多数情况下,实际地形难以满足地形水平的假设,这就需要进行坐标旋转,以消除平流项的影响。当安装铁塔的斜坡坡度特别大时,可以考虑将原本应水平安装的超声风速计调整为与地面平行。3、塔及塔附属设施的影响3.1塔体本身塔本身对观测的影响可分为塔本身对风场的影响,以及塔的偏转、震荡对测量过程的影响两种。3.1.1 对风场的影响自然气流无论是经过几十米的观测塔,还是遇到几毫米的仪器翼梁或电缆,各种尺度的障碍物都会使流线发散,从而导致用于计算通量的流线分离,称为流体失真,流动失真以难以看见的方式影响测量,其影响只能在塔的设计建造阶段进行最小化。在塔的迎风侧(上游),风速受到影响会有所降低。受流动失真影响的逆风距离与障碍物大小的立方成比例,并随着距离的立方体而减小(Wyngaard, 1981, 1988)。在塔的背风侧(下游),风速也减弱,这种效果随着风速的增加而减小(湍流的更快速重构)并且受到障碍物的长度和宽度的影响。图 7 展示了在高塔的迎风侧观察到的风向上的偏转与加速, 图 8则展示了高塔顶部和底部方向迥异的风向。这是由于在背风侧下方产生的回流区造成的,障碍物(塔)尺寸越大,回流区就越容易发展得更大。在塔基通量观测中,森林生态系统的观测常需要10m以上的高塔作为基础,容易导致回流区的产生,回流也增加了向上流动的倾向,并加强了烟囱效应,这可能会显著影响风的测量和干扰混合比梯度。图 7 在塔的迎风侧观察到风向上偏转和加速(引自Sanuki and Tsuda, 1957)图 8 塔顶部的西风流(离地面10米)和离地面2米处的东风回流(引自Vaucher et al., 2004)在建造塔时,尽量选择塔身纤细、结构较少的铁塔,避免对风场的影响,也要注意控制林窗的大小,避免人为形成回流区域。此外,应该尽量减少树木和树枝的移除,因为它们对风的阻力作用可以减少这些回流区域的形成。选择纤细塔体的同时也要保证塔体足够坚固,以确保安全的维护通道和应对整个观测周期中的极端环境。当塔架底座和结构由于受到外界辐射而加热引起对流循环时,可以观察到烟囱效应。这增强了气流的垂直偏转,从而使更多的空气向上移动。烟囱效应取决于基础和塔的质量和热容量、塔的形状、对树冠的干扰程度(清理/切割塔构造的树木)和站点的净辐射量等。烟囱效应是不可避免的,应尽量减少混凝土基础和塔架结构,塔的的横截面也尽量不超过2 x 3 m (Munger et al., 2012)。塔体结构对经过气流的扭曲变形和烟囱效应应该通过专业的方式或通过建模方法(Griessbaum & Schmidt,2009)进行调查(Serafimovich et al., 2011)。3.1.2 对测量过程的影响塔体本身随风速的运动会导致测量中的系统不确定性;塔的移动应限制在0.02 m s-1(即测量风速的精度),并且不应具有在1到20 Hz之间与风向共同变化的力矩(谐波效应);快速响应加速度设备可用于量化塔运动,逐点校正还需要快速响应测斜仪测量以确定旋转速率以及加速度;由于在塔上工作的人员而导致的塔架运动不会随着风或标量交换而变化,但可能会扰乱风场。3.2塔上横臂在1976年的国际湍流对比实验中,一些报告显示直径0.05 m的水平支撑结构造成的平均上升风速为0.1 m/s (Dyer, 1981),它大到足以使涡动相关测量无效。因此,风速计安装臂的尺寸也要尽量小,只需要提供一个安全稳定的测量平台就可以了。王国华等利用成熟的计算流体软件,对布置多个支撑观测仪器的支架所导致的大气边界层风场失真进行定量仿真。他们发现,当支架间距小于6倍的支架直径D或来流风向角小于30°时支架附近流场受到明显的相互干扰。通过对不同来流风向及支架间距离模拟结果的对比分析,认为使用多支架进行多点联合观测时,支架应沿垂直于观测地点常年来流主风向的展向布置。为避免不同支架相互干扰,支架间的最小距离L应大于9倍的支架截面直径。此外,横臂本身需要足够稳定以支撑仪表,可以通过增加侧臂和拉索的方式,以避免横臂的扭矩和振荡。3.3塔下建筑物3.1.1一节讨论了塔体本身对风速和风向造成扭曲从而影响风场的作用,塔下其他障碍物(如设备房间、供电小屋等)也存在这种作用,如图 9 所示。图 9 从障碍物侧面看的迎风流畸变和背风侧流畸变的概念图(引自Davies and Miller, 1982)回流效应在高大的森林冠层中最为明显,但较矮的草地和作物冠层也必须考虑,特别是在附近存放其他设备的房屋的情况下。因此,应尽可能地减少这种流动变形源,在不可减少的情况下,障碍物应远离观测塔,避免对风场的影响。参考文献1. Raupach M R , Finnigan J J . The influence of topography on meteorological variables and surface-atmosphere interactions[J]. Journal of Hydrology, 1997, 190(3-4):182-213.2. Baldocchi D , Falge E , Wilson K . A spectral analysis of biosphere-atmosphere trace gas flux densities and meteorological variables across hour to multi-year time scales. 2000.3. 于贵瑞, 孙晓敏. 陆地生态系统通量观测的原理与方法[M]. 高等教育出版社, 2006.4. Detto M, Katul G G, Siqueira M, et al. The structure of turbulence near a tall forest edge: The backward‐facing step flow analogy revisited[J]. Ecological Applications, 2008, 18(6): 1420-1435.5. Horst T W, Weil J C. How far is far enough?: The fetch requirements for micrometeorological measurement of surface fluxes[J]. Journal of Atmospheric and Oceanic Technology, 1994, 11(4): 1018-1025.6. Gö ckede M, Rebmann C, Foken T. A combination of quality assessment tools for eddy covariance measurements with footprint modelling for the characterisation of complex sites[J]. Agricultural and Forest Meteorology, 2004, 127(3-4): 175-188.7. Munger J W, Loescher H W, Luo H. Measurement, tower, and site design considerations[M]//Eddy Covariance. Springer, Dordrecht, 2012: 21-58.8. Kaimal J C, Finnigan J J. Atmospheric boundary layer flows: their structure and measurement[M]. Oxford university press, 1994.9. Mordukhovich M I, Tsvang L R. Direct measurement of turbulent flows at two heights in the atmospheric ground layer(Atmospheric turbulence statistical characteristics dependence on stratification and elevation from heat flux and wind friction stress characteristics)[J]. ACADEMY OF SCIENCES, USSR, IZVESTIYA, ATMOSPHERIC AND OCEANIC PHYSICS, 1966, 2: 477-486.10. Wyngaard J C. The effects of probe-induced flow distortion on atmospheric turbulence measurements[J]. Journal of Applied Meteorology and Climatology, 1981, 20(7): 784-794.11. Wyngaard J C. The effects of probe-induced flow distortion on atmospheric turbulence measurements: Extension to scalars[J]. Journal of Atmospheric Sciences, 1988, 45(22): 3400-3412.12. Sanukii M, Tsuda N. What are we measuring on the top of a tower?[J]. Papers in Meteorology and Geophysics, 1957, 8(1): 98-101.13. Vaucher G T, Cionco R, Bustillos M, et al. 7.3 FORECASTING STABILITY TRANSITIONS AND AIR FLOW AROUND AN URBAN BUILDING–PHASE I[J]. 2004.14. Griessbaum F, Schmidt A. Advanced tilt correction from flow distortion effects on turbulent CO2 fluxes in complex environments using large eddy simulation[J]. Quarterly Journal of the Royal Meteorological Society: A journal of the atmospheric sciences, applied meteorology and physical oceanography, 2009, 135(643): 1603-1613.15. Serafimovich A, Thomas C, Foken T. Vertical and horizontal transport of energy and matter by coherent motions in a tall spruce canopy[J]. Boundary-Layer Meteorology, 2011, 140(3): 429-451.16. Dyer A J . Flow distortion by supporting structures[J]. 1981, 20(2):243-251.17. 王国华, 贾淑明, 郑晓静. 观测支架引起的大气边界层风场的失真规律[J]. 兰州大学学报: 自然科学版, 2012, 48(5):71-78.Davies M E, Miller B L. Wind effects on offshore platforms-a summary of wind tunnel studies[R]. National Maritime Inst., Feltham (UK), 1982.为了保障各位老师同学从仪器维护的工作中解放出来,做数据的使用者,把更多的时间和精力用在数据深度分析和科学价值发掘方面,我们特提供以下技术服务:站点长期正式运维基于站点管理、工作流程/规范、设备安全、系统优化、设备/数据预警、站点/设备监控、数据分析、科研成果凝练和挖掘等多方面综合执行。站点短期巡检发现目前设备安装、使用、维护、运行状态等影响数据质量的问题。数据远程综汇系统升级建立系统平台,对站点运行状态和数据质量进行预警、监控等。数据整理分析和深度挖掘通过数据整理、插补和分析,形成数据质量分析报告;同时深入挖掘数据背后的科学信息,可以多方面地支撑文章写作、项目申请、专利以及软件著作权申请等工作。通量观测技术培训(涡动相关系统、闪烁仪系统等)根据用户的实际需求,可以有针对性地培训涡动通量观测和设备运行的基本原理,数据处理的基本流程,通量数据处理软件介绍及实际操作演示,通量、气象设备日常维护以及仪器标定,站点选址等相关内容。提供远程视频和上门现场培训等多种方案。