Spatial and Temporal Characteristics of Horizontal Rolls and Cells in the Atmospheric Boundary Layer Based on Radar and in Situ Observations

This paper is the second of a series devoted to the observation and analysis of coherent structures in the cloudy Atmospheric Boundary Layer (ABL) such as horizontal rolls or thermal cells. In the first paper, the TRAC (Turbulence Radar Aircraft Cells) experiment which is the observational support of this investigation based on coupled radar-aircraft measurements, was presented along with an overview of the main results of this campaign held in June 1993 in France. Here the analysis is focused on the spatial characteristics (length-scale, orientation ... ) of the coherent structures, their temporal evolution and vertical distribution deduced from the radar reflectivity fields acquired in clear air at several levels for five different ABLs. For that, an original and efficient image processing method able to extract the major mode of the organisation was developed. These characteristics are examined in relation to the dynamic and thermodynamic state of the ABL using mean and turbulent information as observed by the in-situ aircraft. These experimental results are the basis of a future Large Eddy Simulation modeling of an organised ABL which is in progress and will be the concern of the third paper in the series.     

Atmospheric boundary-layer structure observed during a haze event due to forest-fire smoke

During a haze event in Baltimore, U.S.A. from July 6 to 8, 2002, smoke from forest fires in the Québec region (Canada), degraded air quality and impacted upon local climate, decreasing solar radiation and air temperature. The smoke particles in and above the atmospheric boundary layer (ABL) served as a tracer and provided a unique opportunity to investigate the ABL structure, especially entrainment. Elastic backscatter lidar measurements taken during the haze event distinctly reveal the downward sweeps (or wisps) of smoke-laden air from the free atmosphere into the ABL. Visualisations of mechanisms such as dry convection, the entrainment process, detrainment, coherent entrainment structures, and mixing inside the ABL, are presented. Thermals overshooting at the ABL top are shown to create disturbances in the form of gravity waves in the free atmosphere aloft, as evidenced by a corresponding ripple structure at the bottom of the smoke layer. Lidar data, aerosol ground-based measurements and supporting meteorological data are used to link free atmosphere, mixed-layer and ground-level aerosols. During the peak period of the haze event (July 7, 2002), the correlation between time series of elastic backscatter lidar data within the mixed layer and the scattering coefficient from a nephelometer at ground level was found to be high (R=0.96 for z =324 m, and R=0.89 for z=504 m). Ground-level aerosol concentration was at a maximum about 2 h after the smoke layer intersected with the growing ABL, confirming that the wisps do not initially reach the ground.     

Numerical Study Of The Impact Of Coherent Structures On Vertical Transfers In The Atmospheric Boundary Layer

In two preceding papers, coherent structures of theatmospheric boundary layer (ABL), such as rollvortices or cells, were investigated through radar andaircraft observations collected during the TRAC-93(Turbulence Radar Aircraft Cells) experiment held inFrance in June 1993. The analysis of this experimentaldata set provided information on the spatialcharacteristics of these organisations (length scale,orientation, type ... ), their temporal and verticalevolution, and their relation with the dynamic andthermodynamic conditions of the ABL. For the thirdpaper in this series, a large eddy simulation model is used to examine the impact of thecoherent structures on the ABL vertical fluxes. Theanalysis of the simulated horizontal fields is madewith two-dimensional auto and cross-correlationsapplied on different pertinent ABL variables. Theresults emphasise a directional anisotropy of theseorganised fields throughout the ABL, much morepronounced in the heat flux fields, not only at thelength scale of organisations but also at theturbulence scales. This finding has an importantconsequence for traditional ABL flux measurementsbased on the hypothesis of isotropic and homogeneousturbulence. It can explain part of the underestimationof the surface fluxes often mentioned in theliterature. This approach makes it possible tomodify the concept of diffusion time (in chemicalmodelling) and could also lead to revised ABLparameterisations in Range Scale models.     

大气边界层阵风相干结构的产生条件

壁湍流相干结构的发现是近代湍流研究的重大进展之一,从20世纪50年代开始,在大气边界层湍流中也发现了相干结构——对流云街,并进行了系统的研究。近些年来,人们发现在近地层湍流中也存在相干结构。利用北京325 m气象塔对城市下垫面中大风和小风天气的风速分析,发现较有规律的周期3~6 min的阵风,且有明显的相干结构,而对不同下垫面的阵风研究,均发现存在这种相干结构,这种阵风相干结构对通量输送有不可忽视的作用。本文利用2012年4月甘肃省民勤县巴丹吉林沙漠观测塔的超声风速和平均场风速、温度观测资料,对阵风相干结构的产生条件进行了分析。采用傅立叶变换,将三维超声风速按频率分成基流(周期10分钟以上)、阵风扰动(周期1到10分钟)、湍流脉动(周期小于1分钟)三部分,结合平均场的资料分析发现:阵风相干结构出现在静力中性、不稳定甚至略微稳定的条件下,或者说机械作用主导的大气边界层,阵风区就会出现相干结构,热力作用对其有抑制和干扰的作用。从而,阵风的相干结构和壁面相干结构都出现在中性条件下,是机械湍流的现象,都主导着动量能量的输运。阵风区的相干结构并不等同于对流云街,他们出现在不同的大气稳定度条件下且尺度不同。     

高性能机载云粒子成像仪研制及应用

机载云粒子成像仪是目前直接观测云中粒子谱分布和形状特征的关键设备,在云物理结构探测、遥感反演验证、数值模式云物理过程改进和人工影响天气等领域具有十分重要的作用。因此,针对现有进口仪器性能和应用方面的不足,开发更高性能的云粒子探测仪器十分必要。通过2011年启动的国家重大科学仪器设备开发专项—机载云粒子谱仪与成像仪研制项目,研制出国产高性能机载云粒子测量系统。经过多年的研制和反复测试,解决和显著改进了激光光束均匀化、弱信号探测、多路并行处理及微粒消衍射等关键技术,从而提高了探测器像元的光功率均匀性和一致性,使光斑区域内平均光功率密度提高约3倍,从而改进了粒子成像能力。采用更高分辨率国产光电线阵探测器,测量响应时间明显缩短,对小云粒子的探测能力明显提高。研制的仪器通过一系列基础参数测试、系统优化、环境适应性实验和累计60多架次的飞行测试。对2018年11月5日进行的两架次对比探测结果的初步分析表明,新研制的云粒子成像仪对小粒子浓度的测量精度较进口仪器提高一个量级左右,测量数据起伏更小,稳定性更高,并且可正确获得更清晰可靠的云粒子形态特征。      

Creating Surface Flux Maps from Airborne Measurements: Application to the Mackenzie Area GEWEX Study MAGS 1999

The objective of this study is to produce two-dimensional maps of the sensible and the latent heat fluxes from airborne measurements, based on the analysis of a flight pattern, called grid flights. A footprint model with along-wind and cross-wind components was used to project the measured fluxes onto the surface map. The method was applied to measurements over Arctic tundra during the Mackenzie Area GEWEX (Global Energy and Water Cycle Experiment) Study (MAGS) 1999. The resulting flux estimates were computed by integration of a wavelet transform, and corrected for long wavelength losses using information from 100 km long regional runs that were conducted close to the grid flights. The random flux error was estimated based on the flight length that is represented in each map element, and a map resolution of 3 × 3 km was chosen in order to keep the average random error of the latent heat flux below 25%. The random error of the sensible heat flux was smaller by a factor of 1.4 on average. An analysis of airborne flux measurements at different altitudes showed no significant increase of flux estimates for measurement heights below 90 m. Thus, the fluxes measured at heights between 48 and 64 m were not corrected for vertical flux divergence. The resulting flux maps provide quantitative two-dimensional estimates of the energy exchange between the surface and the atmosphere during the snow melt period in an Arctic environment, which are well-suited for calibration and validation of numerical models.     

Airborne measurements of the impact of ground-based glaciogenic cloud seeding on orographic precipitation

Data from in situ probes and a vertically-pointing mm-wave Doppler radar aboard a research aircraft are used to study the cloud microphysical effect of glaciogenic seeding of cold-season orographic clouds. A previous study (Geerts et al., 2010) has shown that radar reflectivity tends to be higher during seeding periods in a shallow layer above the ground downwind of ground-based silver iodide (AgI) nuclei generators. This finding is based on seven flights, conducted over a mountain in Wyoming (the Unites States), each with a no-seeding period followed by a seeding period. In order to assess this impact, geographically fixed flight tracks were flown over a target mountain, both upwind and downwind of the AgI generators. This paper examines data from the same flights for further evidence of the cloud seeding impact. Composite radar data show that the low-level reflectivity increase is best defined upwind of the mountain crest and downwind of the point where the cloud base intersects the terrain. The main argument that this increase can be attributed to AgI seeding is that it is confined to a shallow layer near the ground where the flow is turbulent. Yet during two flights when clouds were cumuliform and coherent updrafts to flight level were recorded by the radar, the seeding impact was evident in the flight-level updrafts (about 610 m above the mountain peak) as a significant increase in the ice crystal concentration in all size bins. The seeding effect appears short-lived as it is not apparent just downwind of the crest.     

Coherent Structures in the Atmospheric Surface Layer under Stable and Unstable Conditions

Simultaneous measurements of the instantaneous values of absolute temperatureat seven heights within the lower 36 m of the atmospheric boundary layer underdifferent stability conditions were carried out, accompanied by measurements ofthe wind velocity components at two levels and of solar radiation flux at the surface.The data obtained allow one to investigate individual convective cells known ascoherent structures (CS). Outside the CS, i.e., during quiet periods, an instanttemperature profile is in close agreement with the dry-adiabatic lapse rate, butwithin CS the temperature changes much faster with height, and the shape ofthe profile varies significantly.A method was developed to transform temperature records from sensors atseveral heights into an other form, namely, into temporal variations of theheights of isothermal surfaces. Since coherent structures were found to advectwith the mean wind velocity, these temporal height variations may be transformedinto the spatial ones, i.e., into the xoz-plane section of the temperature field.In such a pictorial presentation coherent structures look like asymmetric columnsof heat, penetrating the whole atmospheric surface layer.Coherent structures also exist in the stable stratified surface layer, but they have aninverse asymmetry and occupy only the lower several metres. Wavelike activitydominates in the upper part of the stable surface layer.The characteristic time of surface-layer adjustment to the rapid changes of solarradiation (due to cloud shadows or cloud gaps) was found to be on the order ofone minute. Such a time interval is required for coherent structure to reach the topof surface layer.     

基于旋翼无人机的大气边界层环境气象垂直观测及订正方法的研究

基于旋翼无人机开展大气边界层观测可为气象要素和大气污染物垂直结构的研究提供具有高时空解析能力的新方法,有助于深入理解低层大气物理化学变化机制。本文详述了旋翼无人机在开展大气边界层环境气象垂直观测实验的应用及优势。基于自主研发的旋翼无人机环境气象观测平台,通过开展传感器在无人机上不同的搭载位置,以及旋翼无人机与探空气球、高塔的对比观测实验,明确了旋翼无人机对气象环境观测的影响及合理的搭载方式。研究进一步在湖北重污染天气条件下开展了0—1000 m的大气边界层垂直观测,并研发了基于旋翼无人机姿态数据的大气边界层气象要素及污染物垂直观测的订正方法。结果表明:实验获取了2—10 m垂直分辨率的高质量大气廓线数据,可精细捕捉大气边界层及其逆温层高度和污染物浓度等要素的垂直变化特征。本文旨在为无人机观测的科研应用提供一种技术可行且数据可靠的观测手段。     

‘A Case Study on a Polar Cold Air Outbreak over Fram Strait using a Mesoscale Weather Prediction Model’

The mesoscale weather prediction model ’Lokal-Modell’ (LM) of the Deutscher Wetterdienst is applied to the situation of an Arctic cold air outbreak in the Fram Strait region in April 1998. Observations are available from a flight along 5⁰E carried out during the ARTIST campaign. Initial and time-dependent boundary data for the simulation are taken from a larger scale operational model system. Using the standard configuration of LM, the simulation reproduced the propagation of cold air and the characteristic structure of the atmospheric boundary layer (ABL) in fair agreement with the observations. However, a detailed comparison revealed three basic problems. Firstly, there is evidence that the available data on sea-ice conditions were insufficient approximations to the true state for several reasons. A modification of the sea-ice data towards observations revealed that parts of the discrepancies were due to the original sea-ice data. Secondly, a control run with the model in its standard configuration shows an insufficient warming of the ABL downstream of the ice edge due to underestimation of surface heat fluxes. A simple modification of the approach for the scalar roughness length resulted in the strongest benefit, while comparative studies showed only a slight sensitivity to different types of parametrisation of turbulent mixing or the inclusion of an additional moist convection parametrisation. Thirdly, in all the simulations the deepening of the convective ABL downstream of the ice edge is weaker than observed. This may be partly due to the thermal stratification above the ABL in the analysis data, which is more stable than observed; but it may also be a hint to the fact that processes near the inversion are insufficiently parametrised in mesoscale models with resolutions as used in LM. The simulated cloud layer in the convective ABL is similar to that observed with respect to condensate content, a sharply defined cloud top, a diffuse lower bound, and continuous light precipitation.     

Comparison of Kansas Data with High-resolution Large-eddy Simulation Fields

The surface layer of an atmospheric boundary layer (ABL) is most accessible to field measurements and hence its ensemble-mean structure has been well established. The Kansas field measurements were the first detailed study of this layer, providing numerous benchmark statistical profiles for a wide range of stability states. Large-eddy simulation (LES), in contrast, is most suitable for studying the mixed layer of the ABL where the energy-containing range of the vertical velocity field is well resolved. In the surface layer, typical large-eddy simulations barely resolve the energy-containing vertical-velocity fields and hence do not provide sufficient data for a detailed analysis.We carried out a nested-mesh simulation of a moderately convective ABL (-z_i/L = 8) in which the lower 6% of the boundary layer had an effective grid resolution of 512³. We analyze the LES fields above the 6th vertical grid level (z = 23 m) where the vertical velocity field has a well formed inertial subrange, for a detailed comparison with the Kansas results. Various terms in the budgets of turbulent kinetic energy, temperature variance, Reynolds stress, temperature flux, and some higher order moments are compared. The agreement is generally quite good; however, we do observe certain discrepancies, particularly in the terms involving pressure fluctuations.     

Impact of terrain heterogeneity on near-surface turbulence structure

This study investigates the impact of terrain heterogeneity on local turbulence measurements using 18 months of turbulence data taken on a 30 m tower at the SIRTA mixed land-use observatory under varying stability conditions and fetch configurations. These measurements show that turbulence variables such as the turbulent kinetic energy or momentum fluxes are strongly dependent on the upstream complexity of the terrain (presence of trees or buildings, open field). However, using a detection technique based on wavelet transforms which permits the isolation of the large-scale coherent structures from small-scale background fluctuations, the study shows that, for all stability conditions, whatever the upstream complexity of the terrain, the coherent structures display universal properties which are independent of the terrain nature: the frequency of occurrence, time duration of the coherent structures, the time separation between coherent structures and the relative contribution of the coherent structures to the total fluxes (momentum and heat) appear to be independent of the upstream roughness. This is an important result since coherent structures are known to transport a large portion of the total energy. This study extends to all stability conditions a numerical study by Fesquet et al. [Fesquet, C., Dupont, S., Drobinski, P., Barthlott, C., Dubos, T., 2008. Impact of terrain heterogeneities on coherent structures properties: experimental and numerical approaches. In: 18th Symposium on Boundary Layers and Turbulence. No. 11B.1. Stockholm, Sweden., Fesquet, C., Dupont, S., Drobinski, P., Dubos, T., Barthlott, C., in press. Impact of terrain heterogeneity on coherent structure properties: numerical approach. Bound.-Layer Meteorol.] conducted in neutral conditions which shows that a reason for such behavior is that the production of local active turbulence in an internal boundary layer associated with coherent structure originating from the outer layer and impinging onto the ground is not sensitive to the nature of the terrain.     

Atmospheric Aerosol Characterization Over Naples During 2000–2003 EARLINET Project: Planetary Boundary-Layer Evolution and Layering

The evolution of the planetary boundary layer and the influence of local circulation phenomena over Naples (southern Italy, 40.838° N, 14.183° E, 118 m above sea level) have been studied by systematic lidar measurements of aerosol optical properties and vertical distributions carried out from May 2000 to August 2003, in the course of the EARLINET project. In particular, our data show the development of aerosol layers typically located in the range between 1,000 and 2,300 m, with variable thickness. The optical properties of the observed layers have been determined. In order to analyse the evolution of the planetary boundary layer, detailed observations of complete diurnal cycles have also been performed. The analysis of lidar measurements of vertical profiles of wind speed and wind direction and air mass back-trajectories allowed us to characterize the sea-breeze circulation influence on both the planetary boundary-layer evolution and the observed aerosol vertical distribution.     

A comparison of tropical oceanic heat fluxes determined by airborne eddy correlation and shipboard bulk aerodynamic techniques

Vertical turbulent transports of latent and sensible heat through the tropical marine boundary layer were measured with ship and airborne instruments in the GATE experimental area. The measurements from the two systems are compared for cases of simultaneous sampling in the same locations during undisturbed weather and during the wakes of convective disturbances. The paired average fluxes from the airborne eddy correlation measurements and the shipboard bulk aerodynamic measurements, for undisturbed weather, are related with correlation coefficients in the range of 0.6 to 0.9; the correlations depend primarily on stratifications of the aircraft data according to sampling altitude (15–153 m msl) and flight pattern relative to the mean wind. The agreement between the two types of measurements is best for ship data paired with aircraft samples from the lowest altitudes (15 and 30 m) and from alongwind rather than crosswind flights, as is appropriate since the stationary ships monitor the alongwind characteristics of the turbulence. The use of long (1–3 hr) versus short (10 min) ship samples did not significantly affect the comparability with the aircraft samples during undisturbed periods.The good baseline comparability demonstrated by the undisturbed-weather comparisons was applied to interpret the results from the wakes behind showers. Although these paired measurements were few, not only the sensible heat, but also the latent heat flux was shown to decrease substantially with altitude from 10 m to at least 150 m in the wakes, where the mixing was relatively intense. Variations in the fluxes were much greater at the higher level than nearer the sea surface.     

Simulation of a Stably Stratified Atmospheric Boundary Layer Using One-Dimensional Turbulence

One-dimensional turbulence (ODT) is a single-column simulation in which vertical motions are represented by an unsteady advective process, rather than their customary representation by a diffusive process. No space or time averaging of mesh-resolved motions is invoked. Molecular-transport scales can be resolved in ODT simulations of laboratory-scale flows, but this resolution of these scales is prohibitively expensive in ODT simulations of the atmospheric boundary layer (ABL), except possibly in small subregions of a non-uniform mesh.Here, two methods for ODT simulation of the ABL on uniform meshes are described and applied to the GABLS (GEWEX Atmospheric Boundary Layer Study; GEWEX is the Global Energy and Water Cycle Experiment) stable boundary-layer intercomparison case. One method involves resolution of the roughness scale using a fixed eddy viscosity to represent subgrid motions. The other method, which is implemented at lower spatial resolution, involves a variable eddy viscosity determined by the local mesh-resolved flow, as in multi-dimensional large-eddy simulation (LES). When run at typical LES resolution, it reproduces some of the key high-resolution results, but its fidelity is lower in some important respects. It is concluded that a more elaborate empirically based representation of the subgrid physics, closely analogous to closures currently employed in LES of the ABL, might improve its performance substantially, yielding a cost-effective ABL simulation tool. Prospects for further application of ODT to the ABL, including possible use of ODT as a near-surface subgrid closure framework for general circulation modeling, are assessed.     

Model evaluation of the atmospheric boundary layer and mixed-layer evolution

In the present study, an attempt is made to assess the atmospheric boundary-layer (ABL) depth over an urban area, as derived from different ABL schemes employed by the mesoscale model MM5. Furthermore, the relationship of the mixing height, as depicted by the measurements, to the calculated ABL depth or other features of the ABL structure, is also examined. In particular, the diurnal evolution of ABL depth is examined over the greater Athens area, employing four different ABL schemes plus a modified version, whereby urban features are considered. Measurements for two selected days, when convective conditions prevailed and a strong sea-breeze cell developed, were used for comparison. It was found that the calculated eddy viscosity profile seems to better indicate the mixing height in both cases, where either a deep convective boundary layer develops, or a more confined internal boundary layer is formed. For the urban scheme, the incorporation of both anthropogenic and storage heat release provides promising results for urban applications.     

大气边界层的研究——从均匀到非均匀

主要从理论、试验角度对大气边界层的研究进行了简要的回顾,尤其对近年来倍受关注的非均匀地表条件下的边界层研究进行了较详细的介绍,并对大气边界层的研究进行了展望。20世纪50年代,Monin和Obukhov提出的Monin-Obukhov相似理论,为现代大气边界层物理学奠定了基础。之后进行的多次大型综合性外场试验,不仅揭示和发现了大气边界层中的许多现象和规律,也验证和发展了Monin-Obukhov相似理论。80年代以后,非均匀地表条件下的大气边界层研究越来越受到关注,一些新的现象被揭示,新的问题被提出,同时也提出了一些新的概念和研究方法。未来10 a,除大气湍流机理等需长期坚持研究外,还应进一步开展非均匀复杂下垫面边界层研究、城市边界层研究和局地边界层气象预报研究。     

北京秋季一次降雪前污染天气的激光雷达观测研究

以2009年11月5~8日北京地区发生的一次特殊天气形势下的重污染天气过程为例,研究分析本次污染特点和大气边界层结构特征以及此天气过程的大气温度和相对湿度结构特点。激光雷达是探测大气边界层及气溶胶的一个高效工具,利用ALS300激光雷达系统测量信号,应用Fernald方法反演大气消光系数,根据反演的气溶胶消光系数的最大突变,即最大递减率的高度来确定大气边界层的高度。利用其观测的退偏比分析大气污染物特性。利用微波辐射计数据,确定大气温度和湿度时空特征。研究结果表明:在本次污染天气下,大气具有很强的逆温结构,逆温最大可达近1 K(100 m)-1,500 m以上的大气相对湿度很低,在这种天气特征下的大气边界层高度在400 m左右,非常稳定。污染结束降雪开始前,大气逆温结构消失,大气湿度大幅度增加,接近饱和。根据lidar(light detection and ranging)退偏比的分析,本次污染天气是一次典型的烟尘类颗粒物的污染,污染具有区域性特点。PM2.5(空气动力学当量直径小于等于2.5μm的颗粒物)与AOT(Aerosol Optical Thickness)之间有明显的线性关系,相关系数达到0.72。该lidar系统能够反演出秋季降雪前本次污染天气背景下北京城区上空的大气污染特性和大气边界层高度。     

The Effects of Vegetation Density on Coherent Turbulent Structures within the Canopy Sublayer: A Large-Eddy Simulation Study

Large-eddy simulation has become an important tool for the study of the atmospheric boundary layer. However, since large-eddy simulation does not simulate small scales, which do interact to some degree with large scales, and does not explicitly resolve the viscous sublayer, it is reasonable to ask if these limitations affect significantly the ability of large-eddy simulation to simulate large-scale coherent structures. This issue is investigated here through the analysis of simulated coherent structures with the proper orthogonal decomposition technique. We compare large-eddy simulation of the atmospheric boundary layer with direct numerical simulation of channel flow. Despite the differences of the two flow types it is expected that the atmospheric boundary layer should exhibit similar structures as those in the channel flow, since these large-scale coherent structures arise from the same primary instability generated by the interaction of the mean flow with the wall surface in both flows. It is shown here that several important similarities are present in the two simulations: (i) coherent structures in the spanwise-vertical plane consist of a strong ejection between a pair of counter-rotating vortices; (ii) each vortex in the pair is inclined from the wall in the spanwise direction with a tilt angle of approximately 45°; (iii) the vortex pair curves up in the streamwise direction. Overall, this comparison adds further confidence in the ability of large-eddy simulation to produce large-scale structures even when wall models are used. Truncated reconstruction of instantaneous turbulent fields is carried out, testing the ability of the proper orthogonal decomposition technique to approximate the original turbulent field with only a few of the most important eigenmodes. It is observed that the proper orthogonal decomposition reconstructs the turbulent kinetic energy more efficiently than the vorticity.