Relaxed Eddy Accumulation Simulations of Aerosol Number Fluxes and Potential Proxy Scalars

Direct eddy-covariance measurements of aerosol number fluxes obtained during the 2007 CHATS field experiment in Dixon, California, USA are compared with relaxed eddy accumulation simulations using temperature and water vapour concentration as proxy scalars. After a brief discussion of the limited time response of the aerosol measurement, the applicability of temperature and water vapour concentration as proxy scalars for aerosol number concentration is investigated by evaluating scalar and spectral correlation coefficients as simple measures of scalar similarity. In addition, the proportionality factor b, which compensates for the use of a constant sampling flow in relaxed eddy accumulation, is derived from the time series of aerosol number, temperature and water vapour, and its variability is analyzed. The reduction of the b factor due to application of a deadband, i.e. the rejection of data when the vertical wind speed is close to zero, is evaluated for all three studied scalars, and compared with published functional relationships. In this study, using temperature or water vapour as proxy scalars for aerosol number shows no advantage over the use of a constant b factor. Thus, it is suggested to apply a deadband H _REA  = w′/σ _w  = 0.6 to 0.8 (where w′ is the vertical velocity fluctuation and σ _w is its standard deviation), to use a theoretical b factor based on a parameterization that includes a stability dependence, and to calculate the deadband effect according to a derived relation for aerosol relaxed eddy accumulation.     

海表动量和热量通量的数值模拟

本文分析了用IAP两层大气环流模式模拟的海表动量和热量通量,并将其同Han等和Esbensen等的气候资料比较.模拟的热量通量与观测估计值有类似的水平分布和季节变化,但在中低纬地区有偏多的热量由海洋向大气输送,尤其1月北半球中纬大气需要从海洋获得过分多的热量;模拟的海表动量通量和气候估计值也类似,但模拟的北半球冬季中高纬西风动量通量中心位置偏东,赤道中西太平洋和大西洋的东风动量通量偏弱,南半球环绕南极的西风带模拟得过分弱(尤其在7月份).本文还检验了基本变量的日变化和日际变化对计算海表动量和湍流扩散热量通量的影响,结果表明用月平均,日平均和每小时的基本量计算的动量和湍流扩散热量通量依次增大,尤其在中高纬地区更明显.     

黑河中上游不同下垫面动量总体输送系数和地表粗糙度对比分析

利用廓线法计算了黑河中上游地区盈科农田站、冰沟稀疏草地站、阿柔牧场站和大冬树山垭口积雪观测站的总体输送系数和地表粗糙度。结果表明,地表粗糙度与植被覆盖度和高度以及下垫面的性质有关,夏季地表粗糙度大小是农田站最大,其次是牧场站和稀疏草地站,高寒草甸站最小。下垫面状况还影响动量总体输送系数对稳定度的依赖程度,地表粗糙度大的地区强于地表粗糙度小的地区。     

Numerical Studies with a Regional Atmospheric Climate Model Based on Changes in the Roughness Length for Momentum and Heat Over Antarctica

A regional atmospheric climate model is used toexamine the effect of changes in the roughnesslengths of momentum (_z0m) and heat (_z0h)on the structure of the lower atmosphere and on thesurface energy fluxes over Antarctica. Fourexperiments were carried out in which _z0mand/or _z0h were altered with respect to acontrol experiment. The changes consisted of (1) alowering of _z0m from a field aggregated froma vegetation map with an orographic correction basedon the European Centre for Medium-Range WeatherForecasts _z0m field, to a constant value of10^-3 m; and (2) a lowering of _z0h from a valueequal to _z0m to a constant value of 10^-3 mor a value dependent on the wind speed via a surfacerenewal model. A reduction of _z0m results in theexpected increase in near-surface wind speed. It alsoresults in an increase in the depth of the layer in whichsouth-easterly near-surface winds prevail, and in adecrease in the strength of the large-scale flow overthe continent, in particular in summer. In theescarpment region a decrease of _z0m is foundto result in too high wind speeds. Surface temperatureson average decrease while atmospheric temperaturesincrease, resulting in an increase of near-surfacestatic stability. Changes in roughness lengths donot significantly change the temperature profiles.The surface fluxes, on average found reduced, aremodelled best by using the _z0h based on thesurface renewal method.     

珠海凤凰山森林下垫面干季和湿季气象要素的对比分析与动量和感热交换系数的参数化研究

利用珠海凤凰山陆气相互作用观测塔站2014年11月至2016年5月的观测数据,对比分析了干湿季森林下垫面能量通量和气象要素的变化特征,分析了在不同稳定度下3个风向范围（315°～45°、45°～135°和135°～225°）的动量和感热交换系数随冠层表面风速的变化特征,并对动量和感热交换系数进行了参数化研究。结果表明:干季感热和潜热通量值相当,湿季潜热远大于感热。干季和湿季的夜晚都出现负感热现象,感热从大气向森林输送。相对湿度的变化幅度大,与该地气象状况密切相关,相对湿度的垂直梯度夜晚较大,白天较小。干季的气温垂直梯度比湿季的明显。风速在冬季变化平缓,夏季变化剧烈,低层风速随高度变化梯度明显,高层较紊乱。各高度风向差异不大。中性和近中性状态下,在风向为315°～45°、45°～135°和135°～225°时,动量交换系数C_dn分别为0.05、0.0055和0.022,感热交换系数C_hn分别为0.0055、0.003和0.004。在稳定和不稳定状态下,动量交换系数C_d、感热交换系数C_h随冠层表面风速v明显发生变化,稳定条件下,C_d、C_h随v的增大而增大;不稳定条件下,C_d、C_h随v的增大而减小。分不同风向对森林冠层C_d、C_h在稳定和不稳定条件下与v的关系进行了拟合,得到了参数化公式。     

Statistics of Surface-Layer Turbulence Over Terrain with Metre-Scale Heterogeneity

Refuge has patchy vegetation in sandy soil. During midday and at night, the surface sources and sinks for heat and moisture may thus be different. Although the Sevilleta is broad and level, its metre-scale heterogeneity could therefore violate an assumption on which Monin-Obukhov similarity theory (MOST) relies. To test the applicability of MOST in such a setting, we measured the standard deviations of vertical (_w) and longitudinal velocity (_u), temperature (_t), and humidity (_q), the temperature-humidity covariance (¯tq), and the temperature skewness (S_t). Dividing the former five quantities by the appropriate flux scales (u_*, _*, and q_*) yielded the nondimensional statistics _w/u_*, _u/u_*, _t/|t_*|, _q/|q_*|, and ¯tq/t_*q_*. _w/u_*, _t/|t_*|, and S_t have magnitudes and variations with stability similar to those reported in the literature and, thus, seem to obey MOST. Though _u/u_* is often presumed not to obey MOST, our _u/u_* data also agree with MOST scaling arguments. While _q/|q_*| has the same dependence on stability as _t/|t_*|, its magnitude is 28% larger. When we ignore ¯tq/t_*q_* values measured during sunrise and sunset transitions – when MOST is not expected to apply – this statistic has essentially the same magnitude and stability dependence as (_t/t_*)². In a flow that truly obeys MOST, (_t/t_*)², (_q/q_*)², and ¯tq/t_*q_* should all have the same functional form. That (_q/q_*)² differs from the other two suggests that the Sevilleta has an interesting surface not compatible with MOST. The sources of humidity reflect the patchiness while, despite the patchiness, the sources of heat seem uniformly distributed.     

Determination of Roughness Lengths for Heat and Momentum Over Boreal Forests

The roughness length for momentum (z_0m), zero-plane displacementheight (d), and roughness length for heat (z_0h) are importantparameters used to estimate land-atmosphere energy exchange. Although many different approaches have been developed to parameterizemomentum and heat transfer, existing parameterizations generally utilizehighly simplified representations of vegetation structure. Further, a mismatch exists between the treatments used for momentum and heat exchange and those used for radiative energy exchanges. In this paper, parameterizations are developed to estimate z_0m, d, and z_0h for forested regimes using information related to tree crown density and structure. The parameterizations provide realistic representationfor the vertical distribution of foliage within canopies, and include explicit treatment for the effects of the canopy roughness sublayer and leaf drag on momentum exchange. The proposed parameterizationsare able to realistically account for site-to-site differences in roughness lengths that arise from canopy structural properties.Comparisons between model predictions and field measurements show good agreement, suggesting that the proposed parameterizations capture the most important factors influencing turbulent exchange of momentumand heat over forests.     

Optimized Computation of Transfer Coefficients in Surface Layer with Different Momentum and Heat Roughness Lengths

A new method for computing the surface transfer coefficients is proposed, based on state-of-the-art empirical flux-profile relationships. The influence of the roughness length ratio is first demonstrated with the classical iterative calculation method. Then a non-iterative algorithm is developed, taking into account the difference between momentum and heat roughness lengths.The new method is validated by comparison with the reference iterative computation. The large gain-in computer processing time (CPU) time gain for the calculation of surface fluxes in Eulerian grid models is finally assessed.     

Examination of a Correction Procedure for the Flow Attenuation in Orthogonal Sonic Anemometers

In this study we examine the calibration of wind speed measured by a sonic anemometer, with an orthogonal probe configuration, with regard to the approximation of the flow angle, and the wind-speed dependence of the flow attenuation attributed to the transducer shadow. The flow angle should be calculated by the iterative method when the attenuation is relatively high. For a probe manufactured by Kaijo Co. TR-61C, the wind-speed dependence of the transducer shadow effect is formulated from the results of the wind-tunnel experiment. Assuming the equation is applicable to field observations, significant errors possibly remain especially when the wind speed is low, and /or the angle between the flow vector and the sonic path is small, if the wind-speed dependence in measurement errors is neglected.     

On the Parameterisation of the Effective Roughness Length for Momentum Transfer over Heterogeneous Terrain

We develop a parameterisation for the effective roughness length of terrain that consists of a repeating sequence of patches, in which each patch is composed of strips of two roughness types. A numerical model with second-order closure in the turbulent stress is developed and used to show that: (i) the normalised Reynolds stress develops as a self-similar profile; (ii) the mixing-length parameterisation is a good first-order approximation to the Reynolds stress. These findings are used to characterise the blending layer, where the stress adjusts smoothly from its local surface value to its effective value aloft. Previous studies have assumed that this adjustment occurs abruptly at a single level, often called the blending height. The blending layer is shown to be characterised by height scales that arise naturally in linear models of surface layer flow over roughness changes, and calculations with the numerical model show that these height scales remain appropriate in the nonlinear regime. This concept of the blending layer allows the development of a new parameterisation of the effective roughness length, which gives values for the effective roughness length that are shown to compare well with both atmospheric measurements and values determined from the second-order model.     

Analysis and numerical simulation of the meteorological observations at a tropical coastal site Chennai in India

Meteorological measurements were carried out at North Chennai semi rural area during pre-monsoon period as a part of an air quality study program. Analysis of the data showed the effects of coastal terrain namely the land-sea breeze circulation, temperature cooling during the sea breeze, difference in onset times at these sites etc. Sea breeze onset was observed with a sharp turning of the wind from westerly to south easterly associated with rise in wind speed. Advection speed of the front was about 2.0 m s^− 1. A simple mesoscale meteorological model (MAM-I) developed at Kalpakkam for coastal atmospheric dispersion estimation was used to simulate the observed characteristics. All the major features observed could be simulated by the model while significant difference was noticed in sea breeze frontal movement. MAM results were also inter-compared with MM5. There were no significant differences in the estimate of mean parameters by both the models. It is concluded that the simple model, which takes less run time in a desktop PC, is adequate enough for practical application of providing wind field for plume dispersion models at coastal sites.     

风杯风速仪过高效应的订正研究

该文基于对风杯风速仪动力方程的数值计算和风杯风速仪与超声风速仪平行对比的观测结果,发现专用于近地面层精细风廓线测量的轻型低阈值风杯风速仪测风的过高效应引起的u-error约为1%,因风速平均方法不同（标量平均和矢量平均）引起的DP-error是风杯风速仪测风过高的主要原因,在小风情况下更为明显。DP-error经修正后可获得合理的风速值。     

Power Spectra and Cospectra for Wind and Scalars in a Disturbed Surface Layer at the Base of an Advective Inversion

This paper reports power spectra and cospectra of windspeed and several scalars measured at two heights nearthe base of an advective inversion. The inversion hadformed over a paddy field downwind of an extensive dryregion. Winds over the paddy field were variable instrength and direction, as a result of convectivemotions in the atmospheric boundary layer passing overfrom the dry region upwind. Fetch over the rice waslarge enough that advective effects on the transportprocesses were small at the upper level and negligibleat the lower level. Results from the lower level areinterpreted in terms of a horizontally homogeneous,but disturbed, surface layer.Power spectra of longitudinal and lateral velocitywere substantially enhanced at low frequencies. Theresulting vertical motions added only a small amountto the spectrum of vertical velocity but this stronglyaffected scalar power spectra and cospectra. Thesewere all substantially enhanced over a range of lowfrequencies. We also found that differences in lowerboundary conditions cause differences among scalarspectra at low frequencies.Our analysis shows that the spectra and cospectra havethree components, characterized by different scalingregimes. We call these the ILS (inner-layer scaling),OLS (outer-layer scaling) and CS (combined scaling)components. Of these, the CS component had notpreviously been identified. We identify CS componentsof spectra by their independence of height andfrequency. Spectra with these characteristics had beenpredicted by Kader and Yaglom for a layer of theatmosphere where spectral matching between ILS and OLSwas proposed. However, we find that the velocity andscalar scales used by Kader and Yaglom do not fit ourresults and that their concept of a matching layer isincompatible with our application. An alternativebasis for this behaviour and alternative scales areproposed.We compare our decomposition of spectra into ILS, CSand OLS components with an extended form of Townsend'shypothesis, in which wind and scalar fluctuations aredivided into active and inactive components. Wefind the schemes are compatible if we identify all OLSspectral components as inactive, and all CS and ILScomponents as active.By extending the implications of our results toordinary unstable daytime conditions,we predict that classical Monin–Obukhovsimilarity theory should be modified. We find that theheight of the convective boundary layer is animportant parameter when describing transportprocesses near the ground, and that the scalar scalein the ILS part of the spectrum, which includes theinertial subrange, is proportional to observationheight times the local mean scalar gradient, and notthe Monin–Obukhov scalar scale parameter. The formerdepends on two stability parameters: the Monin–Obukhovstability parameter and the ratio of the inner-layerand outer-layer velocity scales. The outer-layer scalecan reflect disturbances by topographically-inducededdying as well as by convective motions.     

Stability Dependence of Height Scales and Effective Roughness Lengths of Momentum and Heat Transfer Over Roughness Changes

It is widely accepted that the correct formulation of an effective roughness length, defined as the area average of the roughness length in heterogeneous terrain, relies upon the appropriate de-termination of a height scale. At this height a meteorological quantity is approximately in equilibrium with local surface conditions and independent of horizontal position. This research note determines explicitly the different height scales from the perturbation solutions of flow velocity and temperature, as well as the fluxes of momentum and heat, in a stratified boundary layer. These solutions are derived from the asymptotic approximation theory and shown to capture major characteristics of momentum and heat transfer over heterogeneous terrain with changes of the underlying roughness lengths. The effective roughness lengths can then be computed by use of these height scales. The dependence of height scales and effective roughness lengths upon stratification is also discussed briefly.     

敦煌戈壁不同近地层大气稳定度下的感热交换系数

首先利用空气动力学法,计算分析了敦煌戈壁2000年9月—2001年9月的感热交换系数Ch的特征,计算得到Ch的年平均值为0.00203±0.00045。敦煌戈壁的Ch在近中性和中性条件下(-0.01<ζ<0.01,ζ为Monin-Obukhov大气稳定度参数)为常数,Chn≈0.0020;在稳定条件下(ζ≥0.01),Ch比较小,随地表风速的增大而增大,趋于Chn,它们之间存在的关系是:Ch=0.0012 0.0003ln(V);在不稳定条件下(ζ≤-0.01),Ch比较大,随地表风速的增大而减小,也趋于Chn,它们之间存在的关系是:Ch=0.0026-0.0002ln(V)。本文还利用地表温度、一层气温和一层风速近似计算了总体Richardson数Rb,判定大气的层结稳定度,再根据该稳定度来选择感热交换系数。     

Analysis of Model Results for the Turning of the Wind and Related Momentum Fluxes in the Stable Boundary Layer

The turning of wind with height and the related cross-isobaric (ageostrophic) flow in the thermally stable stratified boundary layer is analysed from a variety of model results acquired in the first Global Energy and Water Cycle Experiment (GEWEX) Atmospheric Boundary Layer Study (GABLS1). From the governing equations in this particular simple case it becomes clear that the cross-isobaric flow is solely determined by the surface turbulent stress in the direction of the geostrophic wind for the quasi-steady state conditions under consideration. Most models indeed seem to approach this relationship but for very different absolute values. Because turbulence closures used in operational models typically tend to give too deep a boundary layer, the integrated total cross-isobaric mass flux is up to three times that given by research numerical models and large-eddy simulation. In addition, the angle between the surface and the geostrophic wind is typically too low, which has important implications for the representation of the larger-scale flow. It appears that some models provide inconsistent results for the surface angle and the momentum flux profile, and when the results from these models are removed from the analysis, the remaining ten models do show a unique relationship between the boundary-layer depth and the surface angle, consistent with the theory given. The present results also imply that it is beneficial to locate the first model level rather close to the surface for a proper representation of the turning of wind with height in the stable boundary layer.     

城市边界层动量和保守物通量的特征

利用2005年1-5月北京325 m气象塔47 m高度的湍流脉动资料(风速、温度、水汽和CO2),对城市边界层冠层内的湍流运动统计特征(相似关系、高阶矩、通量和谱等)进行了分析。其中,谱分析的结果表明,城市冠层内稳定度对湍流谱的影响比较小,而水平风速的影响比较大。因此,速度和温度的相似关系在夜间稳定条件下也成立。但是,由于水汽和CO2还受其他因素的影响,相似关系并不适用。更高阶矩的研究表明它们的陡峭度与偏斜度之间存在平方关系。而水汽和CO2之间也存在差异,它们的通量日变化特征明显不同,CO2通量的日变化更能体现人类活动的影响。同时,感热通量、潜热通量和CO2通量存在季节变化,尤其是潜热通量季节差异很大。     

A Note on the Flux-Variance Similarity Relationships for Heat and Water Vapour in the Unstable Atmospheric Surface Layer

Atmospheric surface layer (ASL) experiments over the past 10 years demonstrate that the flux-variance similarity functions for water vapour are consistently larger in magnitude than their temperature counterpart. In addition, latent heat flux calculations using the flux-variance method do not compare as favorably to eddy- correlation measurements when compared to their sensible heat counterpart. These two findings, in concert with measured heat to water vapour transport efficiencies in excess of unity, are commonly used as evidence of dissimilarity between heat and water vapour transport in the unstable atmospheric surface layer. In this note, it is demonstrated that even if near equality in flux-profile similarity functions for heat and water vapour is satisfied, the flux-variance similarity functions for water vapour are larger in magnitude than temperature for a planar, homogeneous, unstably-stratified turbulent boundary-layer flow.     

On the Temporal Distribution of Bursts and Events in the Stable Marine Atmospheric Surface Layer

A methodology for tracing the bursting phenomenon occurring in the turbulent boundary-layer is applied to stable marine atmospheric boundary-layer (MABL) data, collected within the framework of the Coupled Boundary Layer Air-Sea Transfer, Low wind component experiment at Nantucket Island, 2003, Massachusetts, USA. For the definition of the events contributing to the bursting phenomenon, the quadrant analysis of the momentum flux with varying threshold value is utilized. The present work aims at the better understanding of the time scales of the microstructures in the MABL, by providing the number of the events and the groups and their mean duration and the time difference of all quadrants of the quadrant analysis, for different thresholds and under different meteorological conditions. The interrelation of the different time scales as well as the correlation between the quadrants are investigated and discussed. Also, it is demonstrated that the procedure of grouping events leads to bursts with a mean duration that is independent of the threshold value, for the most important quadrants of the momentum transfer (2 and 4). It is found that the combination of near-neutral conditions with low wind speed permits the development of larger duration events within the stable MABL. Finally, it is shown that the groups of events maintain adequate time and flux coherence only for quadrants 2 and 4.