城市地表湍流观测数据的后处理及质量控制

涡动相关仪基于其观测原理,要求观测仪器所在地表平坦均匀,并要求湍流平稳及充分发展等.为了准确合理的运用观测数据,必须对原始数据进行严格的后处理及质量控制,城市复杂非均匀下垫面上的湍流观测数据进行后处理及质量控制更为重要.本文对2010年7月25日-8月5日架设在南京市委党校33 m高度上的涡动相关观测数据进行了包括野点的剔除、坐标旋转修正、WPL修正以及夜里静稳条件下的修正等后处理及质量控制.结果表明:由于观测期间天气状况较好,本次观测的数据质量相对较好,从而野点值剔除,对地表感热和潜热通量测量结果影响甚微.使用二次坐标旋转及平面拟合方法的结果对通量影响较为一致.而且,为了潜热和CO2通量数据的准确性,空气密度效应订正是必须的.本文观测数据基于夏季城市下垫面,所有时次湍流呈现不稳定层结,经过WPL订正,CO2通量增幅为6.5％,潜热通量增幅为12％;观测期间风速不大,主要风频方向上观测点附近无明显高大建筑,且观测点位置较高,σu,v,w/u*与-z/L基本上满足1/3次方律.经湍流平稳性IST检验,感热及动量通量质量较好.     

珠海凤凰山陆气相互作用观测塔通量数据的质量控制与评价

针对广东省珠海凤凰山森林生态系统2015年的湍流原始数据,分析了各修正过程对通量的影响。结果表明:野点较多,需要考虑去除;频率响应修正对通量结果影响不大;倾斜修正的影响在3%左右;超声虚温修正对感热通量的影响变化幅度达30%;WPL修正对碳通量结果影响可达25%左右。经过对观测数据的质量评价,感热通量约有85%,潜热通量和碳通量约有80%的数据可用于通量分析。感热通量在1年中的变化幅度不大;潜热通量变化显著,尤其在季风后迅速增大;碳通量变化与森林下垫面密切相关,夏季变化幅度大,冬季变化较小。      

基于通量方差法对青藏高原东侧高寒草甸地表通量的估算

利用位于青藏高原东侧理塘大气综合观测站2008年观测资料,分析了高寒草甸下垫面上地表通量的时间变化特征,确定了温度、水汽和CO₂的归一化标准差在不稳定情况下随稳定度变化的通量方差关系,应用通量方差法对感热、潜热和CO₂通量进行了计算,并与涡旋相关系统的观测结果进行了比较。结果表明：地表通量月平均日变化呈较为规则的日循环特征,季节变化特征也很明显,雨季（5-9月）潜热大于感热,干季则以感热为主,CO₂通量以6-9月值最大。在不稳定条件下,温度、水汽和CO₂的归一化标准差随稳定度的变化均满足-1/3规律,其通量方差相似性常数分别为1.2,1.4和0.9。通量方差法估算出的通量值与涡旋相关观测得到的通量值有较好的一致性,但感热通量的效果优于潜热通量和CO₂通量。该方法高估了感热通量尤其是潜热通量,而低估了CO₂通量。采用直接观测的感热通量值计算潜热通量和CO₂通量可改善计算结果。     

南疆沙漠腹地大气边界层湍流通量特征的观测研究

利用新疆塔中站2006年4月、8月的三维风速。温度和水汽脉动资料,运用涡旋相关法计算得到了春、夏季塔中10m高度的动量、感热和潜热通量。结果表明,塔中地区地表热量输送以感热输送为主。春季每天的最大感热通量变化范围为120—320W·m^-2,月平均值为220W·m^-2;夏季最大感热通量的变化范围为140—340W·m^-2,月平均值为230W·m^-2。感热通量值在夜间为负,白天为正,符号的改变出现在日出、日落前后。夏季潜热通量最大值一般为20—60W·m^-2,平均值为27W·m^-2,潜热通量比感热通量小一个量级。春季动量通量的平均值为-0．063W·m^-2,夏季动量通量的平均值为-0．091W·m^-2。日变化规律比较明显,日出后,动量向下传输增大,在09-10时（地方时）出现一个最大值,随后动量向下传输并开始减小。     

南沙海域大气湍流通量输送特征分析

利用１９９４年９月在南沙群岛的渚碧礁（１０°５２′Ｎ，１１４°０４′Ｅ）观测的湍流资料，采用涡旋相关法计算，给出了南沙海域的曳力系数，动量通量，感热通量和潜热通量，并与其它海气相互作用实验及ＨＥＩＦＥ的结果作了比较，结果表明，该海域的曳力系数ＣＤ＝（１．５４±０．２４）×１０＾－３。大气向海洋输送动量通量与水平风速挟方成正比，该海区无论白在还是晚上都是将其贮存的热量以潜热或感热的形式输送到大气中。     

北京325米气象塔上CO2梯度观测数据质量控制与评价

近些年涡度相关系统在城市通量研究中得到了广泛应用,因城市下垫面的特殊性和复杂性,以及系统观测原理和仪器精度存在着局限性,使得结果存在10%以上偏差,故必须对原始观测数据进行质量控制和评价。本文对北京325 m气象塔七层高度上的CO₂(二氧化碳)、能量等通量进行了长期观测,研究了涡度相关技术在城市环境通量观测中的适用性,设计出了一套适合于复杂城市下垫面上的涡度相关系统资料前处理和质量控制方案,并对比了不同修正方案效果。分析结果表明,城市下垫面环境中计算CO₂通量的最优周期为30 min;二次坐标旋转法优于平面拟合法;频率响应修正后的感热通量提高5.21%,潜热通量和CO₂通量均提高9.42%。空气密度脉动修正Liu法(Liu,2005)优于WPL法(Webb et al.,1980);湍流谱在惯性副区满足-2/3次方定律,协方差谱满足-4/3次方定律。经过这套前处理和质量控制方案,原始数据中有79%能够用于基础研究,该质量控制与评价体系可为复杂城市下垫面通量研究提供参考。     

高温和台风影响下的苏州城市粗糙子层湍流通量特征

利用苏州地区2011年12月20日—2012年8月13日的湍流观测资料对不同季节、高温、台风强天气过程下的湍流特征进行分析。结果表明:城市地区不同季节动量通量、感热通量、潜热通量日变化明显,各通量的夏季平均值、最大值均高于冬春季,夏季感热通量日最大值为160.2 W·m^-2,感热在城市地表能量平衡中的作用大于潜热,各季节潜热通量平均值仅为感热通量的40%~45%。降水量和植被覆盖度影响地表能量平衡,尤其影响地表热量在感热和潜热之间的分配。在高温天气过程中,感热通量增加明显,其峰值约是夏季平均的1.93倍。由于水汽较少,潜热通量明显减少,约为夏季日平均值的60%。速度三分量谱中u谱与w谱在低频区存在两个峰值,说明在城市复杂下垫面里,湍流激发机制中存在低频过程的影响。在台风天气过程中,动量通量大且变化快,感热输送弱,潜热输送波动大。速度谱w基本不符合"-5/3"次律,惯性子区最小且向高频移动,这和台风内部的复杂上升下沉气流有关。     

2002年南海夏季风爆发期间南海北部海气通量分析与比较

基于2002年南海夏季风期间西沙海气通量观测试验, 采用梯度廓线法计算了观测期间的动量通量、潜热通量、感热通量及其相关要素, 并进行了分析与对比.结果表明: 南海夏季风爆发前后, 南海北部海气通量及其相关要素等都存在显著差异, 不同天气形势下也有较大变化; 南海北部海域在全球海气相互作用中的影响可能不如西太平洋显著; 平均而言, 三次试验的南海北部海气通量变化较小, 感热通量与潜热通量之间存在比较稳定的比例关系, 鲍恩比(Bowen Ratio)维持在大约0.05.     

青藏高原典型下垫面地表能量通量的模型估算与验证

青藏高原地区地表能量通量的估算与验证对高原及其周边地区能量和水循环研究具有重要意义,地表能量平衡系统SEBS(Surface Energy Balance System)模型为研究高原非均匀地表区域地表能量通量提供了一种行之有效的方法。基于中国科学院那曲高寒气候环境观测研究站(简称那曲站)、中国科学院纳木错多圈层综合观测研究站(简称纳木错站)和中国科学院珠穆朗玛大气与环境综合观测研究站(简称珠峰站) 2008年辐射资料、大气边界层塔站观测资料,结合MODIS卫星数据,利用SEBS模型估算地表能量通量,并用站点地表能量通量观测资料进行验证。结果表明,模型估算的感热通量和土壤热通量与站点实测值具有较好的一致性,且感热通量和土壤热通量的估算精度明显优于潜热通量;感热通量的估算精度最高,那曲站、纳木错站和珠峰站的均方根误差分别为54. 98,37. 37和27. 10 W·m~(-2);而模型估算的潜热通量验证结果偏差较大和站点实测数据存在"能量不闭合"问题相关。鉴于在地表能量通量观测中广泛存在"能量不闭合"的问题,利用波文比校正方法校正站点实测潜热通量。研究表明波文比校正方法可以明显改善地表通量观测数据"能量不闭合"的问题,那曲站、纳木错站和珠峰站的能量闭合率分别提高了19. 4%,21. 4%和19. 1%;与原始站点实测潜热通量相比,校正后的潜热通量与SEBS模型估算结果一致性较好,3个站点潜热通量的均方根误差分别减少了6. 78,33. 48和29. 30 W·m~(-2)。     

海-气界面通量计算模式的敏感性试验研究

利用<2000年南海海气通量观测>外场观测资料对海气通量计算模式进行了敏感性试验.试验分析结果表明:感热通量对观测要素的观测误差比较敏感,因而其计算结果的误差较大;相反,潜热通量和动量通量对观测要素的观测误差相对不太敏感,因此,它们的计算结果比较可靠.试验结果也证实了表面通量计算的合理性,其结论可用于对实测资料进行误差检测.     

Surface turbulent flux measurements over the Loess Plateau for a semi-arid climate change study

In order to provide high quality data for climate change studies, the data quality of turbulent flux measurements at the station of SACOL (Semi-Arid Climate & Environment Observatory of Lanzhou University), which is located on a semi-arid grassland over the Loess Plateau in China, has been analyzed in detail. The effects of different procedures of the flux corrections on CO₂, momentum, and latent and sensible heat fluxes were assessed. The result showed that coordinate rotation has a great influence on the momentum flux but little on scalar fluxes. For coordinate rotation using the planar fit method, different regression planes should be determined for different wind direction sectors due to the heterogeneous nature of the ground surface. Sonic temperature correction decreased the sensible heat flux by about 9%, while WPL correction (correction for density fluctuations) increased the latent heat flux by about 10%. WPL correction is also particularly important for CO₂ fluxes. Other procedures of flux corrections, such as the time delay correction and frequency response correction, do not significantly influence the turbulent fluxes. Furthermore, quality tests on stationarity and turbulence development conditions were discussed. Parameterizations of integral turbulent characteristics (ITC) were tested and a specific parameterization scheme was provided for SACOL. The ITC test on turbulence development conditions was suggested to be applied only for the vertical velocity. The combined results of the quality tests showed that about 62%–65% of the total data were of high quality for the latent heat flux and CO₂ flux, and as much as about 76% for the sensible heat flux. For the momentum flux, however, only about 35% of the data were of high quality.     

EVALUATION OF CORRECTIONS ON TURBULENT FLUXES OBTAINED BY EDDY COVARIANCE METHOD IN HIGH WINDS

This paper investigates processing of fast-response data and corrections of turbulent fluxes obtained by using eddy covariance method based on data collected at an offshore observation tower during three Cold-intrusion(CI)events in the South China Sea in 2010. This study presents the data processing procedure in detail and compares frictional velocities(u*), sensible heat fluxes(H) and latent heat fluxes(LE) yielded by using different averaging periods and different coordinate rotation methods; evaluates the sonic temperature correction for sensible heat flux and the Webb correction for latent heat flux as a function of 10 m wind speed(u10) during the CIs. The results show(1) that the different averaging periods of 30 min and 10 min cause biases of u*(H, LE) within 5%(15%, 62%). The values of u*(H,LE) averaged from 30 mins are mostly larger than those averaged from 10 mins. We suggest that the averaging period of 10 min is not sufficiently long to capture all scale eddies and recommend 30 min averaging period in calculating turbulent fluxes using eddy covariance method during CIs;(2) that the values of u*(H, LE) obtained from double rotation(DR2) and those obtained from planar fit rotation(PF) have good agreements and correlation coefficients between them are larger than 0.99. Because PF method requires unchanged environment and it is easier to apply DR2 method, we suggest DR2 coordinate rotation method in processing fast-response data; and(3) that the median values of frictional velocity(sensible heat flux and latent heat flux) binned according to 2 m s~(-1) intervals of u_(10) increase(decrease,increase) by less than 9%(4%, 10%) by Coriolis corrections(sonic temperature corrections, Webb corrections), which decreases(decreases, increases) with increasing u10 when u10 are 5-17 m s~(-1).     

高山草甸下垫面夏季近地层能量输送及微气象特征

利用青藏高原东坡理塘站2007年6～8月的观测资料, 分析了高原东坡草甸下垫面夏季近地层气象要素和湍流通量日变化特征, 并用涡动相关法估算地面的曳力系数。结果表明: 水平风速、 动量通量、 摩擦速度等均在下午最大, 早晨最小。二氧化碳浓度表现为早晚高、 中午低的日变化特征, 比湿的最大值出现在早晨。地表辐射、 热量平衡各分量最高值出现在中午, 最低值出现在早晨。地表反照率表现出早晚高中午低的 “U” 型分布, 日平均值为0.164。夏季地面热源强度在白天午后表现为强的热源, 在夜里表现为弱的冷、 热源交替出现。夏季近地层地气热量交换中, 感热输送作用小, 潜热输送占主要地位。     

沙漠地区沙尘天气近地层湍流输送特征分析

利用1996－04－20－50－10在内蒙古自治区吉兰泰气象站以及1997－04－05－27在内蒙古自治区阿拉善右旗气象站观测到的风，温，湿廓线资料，应用空气动力学方法，计算了每天各观测时次的流动量通量，湍流感热通量和湍流潜热通量，并且比较分析了荒漠下垫面不同沙尘天气状况下的湍流动量通量，感热通量和潜热通量的日变化。结果表明，在沙漠地区春季无降水天气状况下，晴空条件下，白天热力湍流基于动力湍流，在近地层最重要的物质交换是湍流感热交换，浮尘天气条件下，近地层湍流输送较晴空天气条件下弱；扬沙天气条件下，近地层湍流感热交换和湍流动力量流感热通量的峰值大于交换，并且湍流动量通量的峰值较晴空峰值大将近一个量级，湍流感热通量的峰值大于浮尘天气峰值，但小于晴空条件下峰值，沙尘暴天气条件下，湍流动量通量和感热通量一样都是最重要的湍流交换，湍流交换强于扬沙天气，强的感热对沙漠地局地性沙尘暴的产生和加强起着重要作用。     

2014年夏季青藏高原狮泉河与林芝降水低频振荡及陆—气过程日变化特征

根据“第三次青藏高原大气科学试验” 2014年7、8月青藏高原西南部狮泉河站、东南部林芝站的3 m涡动相关系统原始数据和10 Hz湍流资料以及中国气象局台站观测资料、JRA-55（Japanese 55-year Reanalysis）逐日再分析资料、GPCP（Global Precipitation Climatology Project）全球降水逐日观测资料,分别讨论了这两个站在10~20天低频振荡的天气背景下其干、湿位相近地层气象要素的日变化特征以及湍流变化特征。结果表明:两站高低空环流场、水汽通量场、热源的10~20天低频分量在其干、湿位相期间的配置相反。低频地表感热和潜热的不同变化对降水的影响分别在高原西部和东部有不同表现。狮泉河站的低频振荡在纬向上自西向东传播,而林芝站的低频振荡在纬向上自东向西传播,结果表明这两个站分别存在两种不同起源的低频振荡。两站干、湿位相的近地面气象要素以及湍流通量具有明显的日变化特征,通常温度极大值出现在午后14时（北京时,下同）,但狮泉河站干、湿位相的温度极大值均出现在夜间20时;由波文比可知,狮泉河站湿位相全天以潜热为主导,干位相期间,06时之前以潜热为主,06时之后以感热为主;林芝站干、湿位相均为08时之前以感热为主,08时之后以潜热为主。两站湍流平均动能与平均风速正相关,垂直动量表现为向下传输,热量和水汽表现为向上传输。     

Quality analysis applied on eddy covariance measurements at complex forest sites using footprint modelling

Summary Measuring turbulent fluxes with the eddy covariance method has become a widely accepted and powerful tool for the determination of long term data sets for the exchange of momentum, sensible and latent heat, and trace gases such as CO₂ between the atmosphere and the underlying surface. Several flux networks developed continuous measurements above complex terrain, e.g. AmeriFlux and EUROFLUX, with a strong focus on the net exchange of CO₂ between the atmosphere and the underlying surface. Under many conditions basic assumptions for the eddy covariance method in its simplified form, such as stationarity of the flow, homogeneity of the surface and fully developed turbulence of the flow field, are not fulfilled. To deal with non-ideal conditions which are common at many FLUXNET sites, quality tests have been developed to check if these basic theoretical assumptions are valid.In the framework of the CARBOEUROFLUX project, we combined quality tests described by Foken and Wichura (1996) with the analytical footprint model of Schmid (1997). The aim was to identify suitable wind sectors and meteorological conditions for flux measurements. These tools were used on data of 18 participating sites. Quality tests were applied on the fluxes of momentum, sensible and latent heat, and on the CO₂-flux, respectively. The influence of the topography on the vertical wind component was also checked. At many sites the land use around the flux towers is not homogeneous or the fetch may not be large enough. So the relative contribution of the land use type intended to be measured was also investigated. Thus the developed tool allows comparative investigations of the measured turbulent fluxes at different sites if using the same technique and algorithms for the determination of the fluxes as well as analyses of potential problems caused by influences of the surrounding land use patterns.     

青藏高原不同地区夏季近地层能量输送与微气象特征比较分析

利用2009年7月在青藏高原理塘、林芝、海北、拉萨获得的气象观测资料,对比分析了这些地区近地层气象要素、辐射收支及湍流通量日变化特征。结果表明:无论是高原东部、中部还是北部,无论是高原台地还是高山峡谷区,7月份近地层各气象要素、湍流通量、辐射收支都有明显的日变化。各地区的地表辐射、感热、潜热等最高值都出现在中午,最低值出现在早晨。地表反照率日变化均呈早晚高中午低的“U”型分布。地面热源强度在白天均为热源,正午为强热源,在夜间表现为弱的冷、热源交替出现。7月份近地层地气热量交换中,感热输送作用小,潜热输送占主导地位。动量通量和摩擦速度均在风速较大的下午较大,风速小的早晨小。      

On the Significance of the Webb Correction to Fluxes

For establishing correct mass or energy balances at the Earth's surface, detailed and correct measurements of air constituent fluxes are needed. Flux measurements obtained from the eddy covariance technique have to pass several corrections of different relevance in order to give correct flux data. One of these corrections, the Webb correction, is analysed herein from latent heat flux and CO₂ flux data recorded during two field experiments. The significance of this correction for the latent heat flux data varies with the air humidity and the Bowen ratio. The correction changes the latent heat flux values only a little, but significantly (by 2 to 3%). For other air constituents (like CO₂), the Webb correction is much more important (20 to 30% of the flux).     

Evaluation of the profile and the resistance method for estimation of surface fluxes of momentum,sensible and latent heat

Summary Hourly lysimetric and micrometeorological data taken over a grass surface at the Meteorological Research Unit, Cardington U.K. have been analysed. A temperature difference and measurements of wind speed at only one height, combined with an independently estimated effective roughness length allowed sensible heat and momentum fluxes determination by the profile method on an hourly basis. The estimates are compared with direct measurements of sensible heat and friction velocity obtained by the eddy correlation method. The sensible and latent heat fluxes are also modelled by the resistance method. Equations based on the Monin—Obukhov similarity theory are used to account for stability effects through various forms of parameterization Aerodynamic and surface resistances, necessary for the Penman—Monteith equation are calculated from routinely measured meteorological data. The profile method for estimation of sensible heat flux and friction velocity is found to work excellently on the discussed daytime experimental data which correspond mainly to near neutral or slightly unstable conditions.Surface latent and sensible heat fluxes can also be described very well by the resistance method. A slightly better estimate of the sensible heat flux is achieved when stability corrections are taken into account. On the contrary Penman-Monteith equation for estimating latent heat flux is insensitive to adjustments for atmospheric stability.The comparison of the various methods leads to the establishment of empirical relationships which correlate various quantities such as soil heat flux, resistances, evapotranspiration etc. to routinely measured meteorological data.With 8 Figures