共查询到19条相似文献，搜索用时 140 毫秒
1.
坐标旋转是涡动协方差(EC)通量观测数据后处理当中必须进行的步骤,本文利用塔克拉玛干沙漠腹地塔中西面标校站自然沙丘上3 m高度和塔中本站32 m高度2009年4—10月湍流通量数据,对比分析了不同坐标旋转方法在塔中通量观测数据中的应用。结果表明,塔克拉玛干沙漠腹地塔中湍流各通量在旋转订正前后差别不显著,不订正造成的通量计算误差1%。塔中32 m和塔中西站3 m通量计算误差,与两个站点的采集仪器架设高度差别较大有关。即使对于相对平坦的下垫面,也必须要进行倾斜校正,以修正地势、仪器振动、安装倾斜等产生的虚假信息。本研究对塔中涡动资料进行旋转时采用了二次坐标旋转方法。此结果为塔克拉玛干沙漠腹地EC通量观测的风速仪倾斜校正提供参考。 相似文献
2.
不同平均时间对LOPEX10资料涡动相关湍流通量计算结果影响的探讨 总被引：1，自引：0，他引：1
通过处理涡动相关系统观测的近地层湍流脉动量可以获取地—气间感热和潜热通量,然而选择不同平均时间对通量计算的结果有较大影响。采用黄土高原陆面过程野外观测试验(LOPEX10)期间获得的涡动相关系统观测资料,分析了不同平均时间对湍流通量计算的影响,并采用雷诺平均和分解方法推导了平均时间引起的通量差值的数学表达式(Flux Compensation,FC)。结果表明:(1)FC公式可以说明采用不同平均时间数据之间的关系,也可以直接计算低频涡旋对湍流通量的贡献。FC公式计算的结果与直接计算的不同平均时间通量计算之差的相关系数在0.95以上,并可以确定计算湍流通量的最佳平均时间。(2)通过采用Ogive函数确定了计算LOPEX10期间通量的最佳平均时间长度为30min,印证了利用雷诺平均和分解方法计算湍流通量补偿的准确性。(3)通过进一步的数学变换,证明了平均时间对湍流通量计算的影响直接与湍流低频变化相关,FC公式可以用来确定涡动相关观测数据的最佳平均时间,并且在获得较高时间分辨率的湍流通量数据的同时,补偿因平均时间过短而遗漏的低频信息。 相似文献
3.
利用中国科学院那曲高寒气候环境观测研究站2013年9月至2014年8月自动气象站(AWS)和涡动相关系统(EC)的观测资料,基于空气动力学法、地表能量平衡组合法、总体输送法以及涡动相关法,计算了高寒草地下垫面的湍流通量,并对不同方法计算结果间的一致性和差异性进行了分析。结果表明,不同方法计算的湍流通量特征具有明显的差异。地表能量平衡组合法满足能量平衡关系,但在早晨和傍晚层结转换期间,计算的湍流通量出现异常不稳定值;空气动力学法计算的湍流通量在整个观测期与涡动相关法计算的湍流通量相关性最高,但在大气稳定度参数接近0,计算结果不稳定;总体输送法计算的通量数据在地气温差为负值时发散明显,但该方法原理简单,适合在只有常规观测项目的业务气象站或在气象观测项目不全的野外台站使用。空气动力学法和地表能量平衡组合法与涡动相关法的湍流通量的平均偏差相对较小,而总体输送法平均偏差相对较大。研究方法和结果除了为这些方法的使用提供参考外,也为建立长时间通量序列提供了一个较为合理的依据,有助于深入了解高原地气相互作用。 相似文献
4.
涡动相关仪获得的湍流数据质量不仅受到传感器构造及当时环境因子的影响,同时可能受到理论假设带来的误差的影响.对2018年1月1日-2月28日冬季增城丘陵灌木林典型下垫面近地层湍流观测数据进行了包括野点剔除、倾斜订正、感热通量的超声虚温订正、WPL修正以及自动增益控制(AGC)修正等后处理及质量控制.结果表明:常规野点剔除方法对于通量计算结果影响较小,AGC对降水和高湿天气十分敏感,能有效的提高能量闭合率;湍流发展充分性检验与平稳性检验显示,可用于基本研究的潜热通量、C02通量、摩擦速度、感热通量分别占样本总量的90.5％,86.4％,83.1％,87.7％,表明通过本研究质控方法处理后的湍流通量观测数据质量可以得到保证,可用于后续的分析. 相似文献
5.
近地层湍流通量计算对于中尺度数值模式有重要意义,湍流通量的参数化是当前大气边界层研究的重要课题之一.选择青藏高原东缘大理观象台边界层通量观测系统,离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA相似理论迭代方案B),并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析.在大理观象台观测场不同植被随季节交替的状况下,根据边界层铁塔4层高度风速拟合,发现近地层空气动力学粗糙度随季节变化特征明显.将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算.结果表明:稳定度是影响近地层参数化方案精度的重要因素,在不稳定条件下方案B低估了动量通量,方案A优于方案B,而在稳定条件下方案A低估了动量通量,方案B优于方案A,两种方案总体来看误差不大.对于大理边界层通量观测场农田植被交替的环境条件,不同季节下垫面植被类型的差异,以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显着影响.方案B考虑了空气动力学粗糙度z0和热量粗糙度z0h的差异,不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A.针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象,使用了Zeng等1998年提出的用辐射地表温度订正裸土下垫面感热能量方法后,计算结果也有明显改善. 相似文献
6.
近地层湍流通量计算对于中尺度数值模式有重要意义, 湍流通量的参数化是当前大气边界层研究的重要课题之一。选择青藏高原东缘大理观象台边界层通量观测系统, 离线测试了WRF区域模式中的两种常用的近地层参数化方案(MM5相似理论非迭代方案A和ETA 相似理论迭代方案B), 并将参数化方案计算结果与边界层铁塔涡动相关法的观测值进行对比分析。在大理观象台观测场不同植被随季节交替的状况下, 根据边界层铁塔4层高度风速拟合, 发现近地层空气动力学粗糙度随季节变化特征明显。将拟合的空气动力学粗糙度输入模式参数化方案进行通量计算。结果表明:稳定度是影响近地层参数化方案精度的重要因素, 在不稳定条件下方案B低估了动量通量, 方案A优于方案B, 而在稳定条件下方案A低估了动量通量, 方案B优于方案A, 两种方案总体来看误差不大。对于大理边界层通量观测场地农田植被交替的环境条件, 不同季节下垫面植被类型的差异, 以及植被的稀疏对近地层参数化方案湍流通量计算结果的精度有显著影响。方案B考虑了空气动力学粗糙度z0和热量粗糙度z0h的差异, 不稳定条件下感热通量计算结果在裸土或稀少植被条件下明显优于方案A。针对方案B不稳定条件下感热通量计算结果在裸土下垫面仍出现高估的现象, 在使用了(Zeng, et al, 1998)提出的对于使用辐射地表温度在裸土下垫面时的订正方法后, 计算结果也有明显改善。 相似文献
7.
本文叙述的一种最优内插与逐步订正结合的高空风场客观分析方法,是以最优内插的高度场求得的地转风作预备场,使用最优内插权重函数经2—3次订正而成的。 通过对三个个例的三个层次(850、500和300毫巴)上测风资料进行风场客观分析试验,得到了资料拟合误差(全场平均)均在1米/秒左右的满意结果。分析出来的流线图不仅光滑,而且根据分析结果计算的涡度场和散度场与实际天气形势配合一致,并且其量级也是合理的。 相似文献
8.
湍流通量计算方法和误差的比较研究 总被引：4，自引：1，他引：3
近地层的风、温、湿梯度以及辐射和土壤热通量资料在过去的几十年间被广泛地用于计算陆-气感热、潜热通量.然而,不同的计算方法之间存在一些差异,随意的选用可能会造成计算结果存在较大的误差.文中利用绿洲系统能量与水分循环过程观测试验的2005年6月甘肃金塔绿洲中部观测资料,在分析了观测期间近地层微气象特征的基础上,对鲍恩比法、空气动力学法(包含两种不同计算方案)和变分法计算的绿洲农田下垫面的湍流通量进行了比较,并通过敏感性实验和理论分析探讨了不同方法间计算结果的差异和误差的来源.结果表明:不同方法计算的湍流通量特征存在明显差异.鲍恩比法从理论上满足能量平衡关系,但在昼夜交替和夜晚时出现了计算不稳定现象,因此在这些时段不适用.变分法避免了鲍恩比法计算夜间湍流通量的不稳定性,而且二者计算的白天湍流通量有较好的相关性;同时变分法由于加入了辐射和土壤热通量信息,明显改善了空气动力学法白天的能量平衡状况;变分法计算得到的潜热通量对辐射和微气象条件有较为合理的响应;对数据误差的敏感性试验表明,鲍恩比法和空气动力学法应用于绿洲下垫面上对观测误差较为敏感,而变分法相对比较稳定.进一步的分析表明,鲍恩比法在-2<β<0的范围内由于方法本身局限而适用性较差;而空气动力学方法计算结果的差异和误差来源于所应用的相似性函数.因而与鲍恩比法和空气动力学法计算的湍流通量值相比较,变分法的计算结果更加合理和稳定,在绿洲下垫面有较好的适用性.该文研究方法和结果可为这些方法的使用提供参考. 相似文献
9.
在y-z平面上,用积云动力方程组和涡度方程联合的二维模式进行数值试验,得到以下结果: 1.切变环境中,对流运动引起水平涡度的机制是辐合作用和倾斜作用。 2.水平涡度的强度和切变强度成正比。 3.倾斜项形成对流单体南北两侧一对反对称的涡度分布,辐合项使该涡度的水平分布更符合实际观测到的结果。 4.大气密度随高减小是涡度最强值能稳定在对流运动中下部的原因。 5.涡度的最强发展阶段的时间尺度和对流运动发展的时间尺度相当,前者落后约3分钟。 计算结果和观测相比是符合的。 相似文献
10.
理论分析和数值试验结果表明,地表参量的不均匀性对网格区地表湍流通量的计算有重要影响。近中性大气层结条件下,次网格粗糙度长度的变差系数、网格平均粗糙度及参考高度的选取是影响网格湍流通量计算的主要因素,其中次网格粗糙度长度的变差系数对计算偏差起主要决定作用。实际计算表明,某些特定地区(如植被气候过渡带)粗糙度的地表非均匀性引起的计算相对误差可达40%以上,选取特定的参考高度能改善高网格湍流通量计算的效果。非中性大气层结条件下,由地表粗糙度不均匀性所致的平均风速、位温梯度以及近地层大气稳定度的次网格分布都对网格湍流通量(感热通量)计算产生影响,比较而言,相对误差大小对大气稳定度的次网格分布最为敏感。所以,在目前的数值模式中有必要进一步对湍流通量计算过程中由于地表不均匀性产生的计算偏差加以考虑。 相似文献
11.
LIU Li WANG Tijian SUN Zhenhai WANG Qingeng ZHUANG Bingliang HAN Yong LI Shu 《大气科学进展》2012,29(6):1264-1278
Turbulence characteristics of an atmospheric surface layer over a coastal mountain area were investigated under different coordinate frames.Performances of three methods of coordinate rotation:double rotation(DR),triple rotation(TR),and classic planar-fit rotation(PF) were examined in terms of correction of eddy covariance flux.Using the commonly used DR and TR methods,unreasonable rotation angles are encountered at low wind speeds and cause significant run-to-run errors of some turbulence characteristics.The PF method rotates the coordinate system to an ensemble-averaged plane,and shows large tilt error due to an inaccurate fit plane over variable terrain slopes.In this paper,we propose another coordinate rotation scheme.The observational data were separated into two groups according to wind direction.The PF method was adapted to find an ensemble-averaged streamline plane for each group of hourly runs with wind speed exceeding 1.0 m s 1.Then,the coordinate systems were rotated to their respective bestfit planes for all available hourly observations.We call this the PF10 method.The implications of tilt corrections for the turbulence characteristics are discussed with a focus on integral turbulence characteristics,the spectra of wind-velocity components,and sensible heat and momentum fluxes under various atmospheric stabilities.Our results show that the adapted application of PF provides greatly improved estimates of integral turbulence characteristics in complex terrain and maintains data quality.The comparisons of the sensible heat fluxes for four coordinate rotation methods to fluxes before correction indicate that the PF10 scheme is the best to preserve consistency between fluxes. 相似文献
12.
Sonic Anemometer Tilt Correction Algorithms 总被引：72，自引：38，他引：72
The sensitivity of sonic anemometer-derived stress estimates to the tilt of the anemometer is investigated. The largest stress errors are shown to occur for unstable stratification (z/L<0) and deep convective boundary layers. Three methods for determining the tilt angles relative to a mean streamline coordinate system and for computing the tilt-corrected stresses are then compared. The most commonly used method, involving a double rotation of the anemometers' axes, is shown to result in significant run-to-run stress errors due to the sampling uncertainty of the mean vertical velocity. An alternative method, requiring a triple rotation of the anemometer axes, is shown to result in even greater run-to-run stress errors due to the combined sampling errors of the mean vertical velocity and the cross-wind stress. For measurements over the sea where the cross-stream stress is important, the double rotation method is shown to overestimate the surface stress, due to the uncorrected lateral tilt component. A third method, using a planar fit technique, isshown to reduce the run-to-run stress errors due to sampling effects, and provides an unbiased estimate of the lateral stress. 相似文献
13.
青藏高原东部草甸下垫面涡旋相关观测数据的质量控制及评价研究 总被引：2，自引：1，他引：1
目前利用涡旋相关仪观测地表通量的方法已被广泛采用, 但由于涡旋相关法是建立在方程假设简化的基础上, 这对观测数据的质量提出了一定的要求, 所以对观测数据的处理及通量结果的质量状况分析就显得尤为重要。本文以理塘站2006年7月中到8月中的数据为例, 对其进行质量控制, 并将质量控制后的通量计算结果与原始值进行比较。质量控制后感热通量和潜热通量总和有所增长, 动量通量值则受质量控制影响较大。质量评价的结果为:动量通量、感热通量、潜热通量中高质量数据所占比例为68.2%、60.6%、63.3%。 相似文献
14.
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). 相似文献
15.
Renzo Richiardone Elena Giampiccolo Silvia Ferrarese Massimiliano Manfrin 《Boundary-Layer Meteorology》2008,128(2):277-302
We describe the coordinate transformations that can be used to convert the velocity components measured by a set of sonic
anemometers with time-dependent tilt fluctuations into a single, time-independent coordinate system. By applying the planar
fit method (PFM) to each anemometer dataset, it is possible, for planar flows, to locate the flow plane at each measurement
point and compare its orientation with the topography. Installation on a ship is also considered. An application of this method
to intercomparison data has led to the detection of an instrument error due to a misalignment between the assembly of the
sonic transducers and the anemometer pedestal. If this error occurs, pedestal levelling does not guarantee that measurements
are unbiased. A correction method is proposed and the results of two experiments are shown. Flow planarity at different levels
and flow distortion caused by the mast are highlighted. The influence of the error on the evaluation of the Reynolds stresses
using PFM or the double rotation method and the triple rotation method is discussed and the tilt corrected stresses calculated
using the three methods compared. 相似文献
16.
Hyun-Mi Oh Kyung-Ja Ha Ki-Young Heo Kyung-Eak Kim Sang-Jong Park Jae-Seol Shim Larry Mahrt 《Asia-Pacific Journal of Atmospheric Sciences》2010,46(4):513-523
This study presents an evaluation of the atmospheric factors influencing the post-processing for fast-response data of horizontal momentum, vertical wind component, temperature, and water vapor to measure turbulent fluxes. They are observed at the Ieodo ocean research station over the Yellow Sea during the period of October 2004 to February 2008. The post process methods employed here are composed of quality control and tilt correction for turbulent flux measurement. The present result of quality control on the fast-response data shows that total removal ratio of the data generally depends on the factors such as a wind speed, relative humidity, significant wave height, visibility, and stability parameter (z/L). Especially, the removal ratio of water vapor data is significantly increased on light wind and strong stability conditions. The results show that the total removal ratio of water vapor data increases when wind speed is less than 3 m s?1 and wave height is less than 1 m. The total removal ratio of water vapor data also increases with the value of the stability parameter. Three different algorithms of tilt correction methods (double rotation, triple rotation, and planar fit) are applied to correct the tilt of the sonic anemometer used in the observation. Friction velocities in near neutral state are greater than friction velocity in other states. Drag coefficients are categorized in terms of stabilities and seasons. 相似文献
17.
Renmin Yuan Minseok Kang Sung-Bin Park Jinkyu Hong Dongho Lee Joon Kim 《Meteorology and Atmospheric Physics》2011,110(3-4):123-133
An expanded planar-fit (PF) method over complex terrain is presented and applied to coordinate rotation of the eddy-covariance (EC) flux and vertical velocity estimation. Theoretical analysis indicates that PF coefficients depend on wind direction, and an expression of vertical velocity is deduced. We applied the theory using 1?year of observations from the KoFlux site in the Gwangneung Forest in Korea and investigated the influence of wind direction on the PF method. Then, we performed an expansion of the PF method to consider dependence of PF coefficients on wind direction and applied the PF method to every sector. The results show that the PF coefficients and tilt angles over complex terrain vary with wind direction. Two hundred 30-min data sets are sufficient to derive stable PF coefficients over hilly terrain for each sector. The relative difference in eddy-covariance flux between the general planar fit (GPF) and sector planar fit (SPF) is less than 10% for the scalar flux and about 18% for friction velocity. Vertical velocity and vertical advection (VA) terms were also calculated and compared using SPF and GPF methods, and a normal distribution and diurnal trend of real vertical velocity on clear days are presented. 相似文献
18.
Lukas Siebicke Martina Hunner Thomas Foken 《Theoretical and Applied Climatology》2012,109(1-2):109-131
Observations of vegetation–atmosphere exchange of carbon dioxide (CO2) by the eddy covariance (EC) technique are limited by difficult conditions such as nighttime and heterogeneous terrain. Thus, advective flux components are included into the net ecosystem exchange (NEE) budget. However, advection measurements are experimentally challenging and do not always help to solve the night flux problem of the EC technique. This study investigates alternative methods for the observation of horizontal advection, in particular horizontal concentration gradients, as well as different approaches to coordinate rotation and vertical advection. Continuous high-frequency measurements of the horizontal CO2 concentration field are employed and compared to the often used discontinuous sequential sampling. Significant differences were found in the case of 30-min mean concentration values between the conventional discontinuous sampling approach and the complete observation of the time series by continuous sampling. Estimates of vertical advection rely on accurate estimates of vertical wind velocity ( $\emph{w}$ ). Therefore, different approaches to the planar fit coordinate rotation have been investigated. Sector-wise rotation was able to eliminate directional dependencies of mean $\emph{w}$ . Furthermore, the effect of the data set length used for rotation (window length) was investigated and was found to have significant impact on estimates of vertical advection, with larger window lengths yielding about 50% larger vertical advection. A sequential planar fit with controlled window length is proposed to give reproducible results. The different approaches to the measurement and calculation of horizontal and vertical advection presented are applied to data obtained during the exchange processes in mountainous region experiment at the FLUXNET site Waldstein–Weidenbrunnen (DE-Bay). Estimates of NEE including advection are compared to NEE from turbulent and storage flux alone without advection. NEE including vertical advection with sector-wise planar fit rotation and controlled window length and including horizontal advection from continuous gradient measurements, which were comprehensively bias corrected by a new approach, did compare well with the expected night flux error, with meteorological drivers of the fluxes and with soil chamber measurements. Unrealistically large and noisy values of horizontal advection from the conventional discontinuous sampling approach, which lead to unrealistic values of NEE, could be eliminated by the alternative approaches presented. We therefore suggest the further testing of those approaches at other sites in order to improve the accuracy of advection measurements and, subsequently, estimates of NEE. 相似文献
19.
The bandpass eddy covariance method has been used to measure the turbulent flux of scalar quantities using a slow-responsescalar sensor. The method issimilar in principle to the traditional eddy correlation method but includes the estimation of high-frequency components of the flux on the basis of cospectral similarity in the atmospheric surface layer. In order to investigate the performance of the method, measurements of the water vapour flux over a forest with the bandpass eddy covariance method and the direct eddy correlation method were compared. The flux obtained by the bandpass eddy covariance method agreed with that by the eddy correlation method within ±20% for most cases, in spite of a rather slow sensor-response of the adopted hygrometer. This result supports its relevance to a long-term continuous operation, since a stable, low-maintenance,general-purpose sensor canbe utilized for scalar quantities. Oneweak point of the method isits difficulty in principle to measure the correct flux when the magnitude of the sensible heat flux is very small, because the method uses the sensible heat flux as a standard reference for the prediction of undetectable high-frequency components of the scalar flux. An advanced method is then presented to increase its robustness. In the new method, output signals from a slow-response sensor are corrected using empirical frequency-responsefunctions for the sensor,thereby extending the width of the bandpass frequency region where components of the flux are directly measured (not predicted). The advanced method produced correct fluxes for all cases including the cases of small sensible heat flux. The advanced bandpass eddy covariance method is thus appropriate for along-term measurement of the scalar fluxes. 相似文献