激光雷达反演边界层高度方法评估及在北京的应用

边界层高度是影响大气边界层发展和空气污染程度的重要因子,是环境和气候研究的重要参数。本文利用激光雷达对北京地区2011年5月至2012年4月的边界层高度进行探测分析,采用小波协方差方法反演边界层高度,评估了该方法的适用性。得到基于小波协方差方法自动判断边界层高度的最优参数组合,激光雷达与飞机探测结果对比一致性较好;与探空结果相关系数0.88,激光雷达反演的边界层高度略偏高。当激光雷达的垂直分辨率为30 m时,更加适合北京地区的步长和阈值分别为210 m和0.05;当激光雷达的垂直分辨率为15 m时,步长和阈值分别为135 m和0.05。分析期间,不同季节边界层高度日变化有明显的不同,夏季14:00（北京时）左右达到最高,较高的边界层高度可维持3~4 h,平均可达1.30 km;冬季较高边界层高度只能维持2 h左右,平均为1.08 km。有云与无云天气边界层日变化特征以及边界层高度存在显著的差异,云的存在减少了到达地面的直接辐射,抑制了湍流的发展,进一步抑制了边界层的发展;本文也将激光雷达反演边界层高度结果应用于观测时期边界层高度与地面污染的关系研究中,统计得到边界层高度与PM2.5浓度的相关系数为－0.340。     

气溶胶激光雷达和无线电探空观测边界层高度的对比分析

利用北京地区2017年11月至2018年1月连续3个月的激光雷达资料和无线电探空数据,按照清洁天、污染天和多云天3种天气条件,对大气边界层高度的计算方法和结果进行对比分析。结果表明,基于激光雷达消光系数的梯度法、标准差法和小波法都能够较好地提取边界层高度。清洁天标准差法计算的边界层高度高于梯度法和小波法,08:00（当地时间,下同）和20:00由无线电探空得到的清洁天边界层高度平均值分别为1176 m和1224 m。污染天标准差法的计算结果要低于梯度法和小波法,污染天无线电探空得到的边界层高度平均值约为956 m,和清洁天相比降低了两百多米,重污染时最低降低至562 m,逆温层高度和PM2.5浓度具有明显的反相关关系。有云时,梯度法和小波法确定的边界层高度和云高非常接近,标准差法计算的结果略低。总体而言,气溶胶激光雷达计算的边界层高度随着污染等级的提高没有明显的降低趋势,相反在重度污染情况下反而有所增加,这可能是由于污染物的不断堆积导致的。梯度法确定的边界层高度易受到污染物传输过程的影响,略高于逆温层高度。另外,激光雷达确定的边界层高度受到残留层影响时,也会高于逆温层。     

Determination of the Atmospheric Boundary Layer Height from Radiosonde and Lidar Backscatter

The height of the atmospheric boundary layer is derived with the help of two different measuring systems and methods. From radiosoundings the boundary layer height is determined by the parcel method and by temperature and humidity gradients. From lidar backscatter measurements a combination of the averaging variance method and the high-resolution gradient method is used to determine boundary layer heights. In this paper lidar-derived boundary layer heights on a 10 min basis are presented. Datasets from four experiments – two over land and two over the sea – are used to compare boundary layer heights from both methods. Only the daytime boundary layer is investigated because the height of the nighttime stable boundary layer is below the range of the lidar. In many situations the boundary layer heights from both systems coincide within ±200 m. This corresponds to the standard deviation of lidar-derived 10-min values within a 1-h interval and is due to the time and space variability of the boundary layer height. Deviations appear for certain situations and depend on which radiosonde method is applied. The parcel method fails over land surfaces in the afternoon when the boundary layer stabilizes and over the ocean when the boundary layer is slightly stable. An automatic radiosonde gradient method sometimes fails when multiple layers are present, e.g. a residual layer above the growing convective boundary layer. The lidar method has the advantage of continuous tracing and thus avoids confusion with elevated layers. On the other hand, it mostly fails in situations with boundary layer clouds     

Frequency of Boundary-Layer-Top Fluctuations in Convective and Stable Conditions Using Laser Remote Sensing

The planetary boundary-layer (PBL) height is determined with high temporal and altitude resolution from lidar backscatter profiles. Then, the frequencies of daytime thermal updrafts and downdrafts and of nighttime gravity waves are obtained applying a fast Fourier transform on the temporal fluctuation of the PBL height. The principal frequency components of each spectrum are related to the dominant processes occurring at the daytime and nighttime PBL top. Two groups of cases are selected for the study: one group combines daytime cases, measured in weak horizontal wind conditions and dominated by convection. The cases show higher updraft and downdraft frequencies for the shallow, convective boundary layer and lower frequencies for a deep PBL. For cases characterized by strong horizontal winds, the frequencies directly depend on the wind speed. The temporal variation of the PBL height is determined also in the likely presence of lee waves. For nighttime cases, the main frequency components in the spectra do not show a real correlation with the nocturnal PBL height. Altitude fluctuations of the top of the nocturnal boundary layer are observed even though the boundary layer is statically stable. These oscillations are associated with the wind shear effect and with buoyancy waves at the PBL top.     

Sensitivity analysis of planetary boundary layer schemes using the WRF model in Northern Colombia during 2016 dry season

This study conducted meteorological simulations in northern Colombia by analyzing different planetary boundary layer (PBL) schemes available in the numerical Weather Research and Forecasting (WRF) model. The study area included three nested domains with horizontal resolutions of 18 km, 6 km, and 2 km, with 38 vertical levels. The evolution and structure of the PBL were analyzed during the driest months (March, April, and May 2016) and in regions with the highest particulate matter concentrations. Sensitivity analysis of the WRF model was performed with two local and two non-local PBL schemes. The results were validated using observations of the surface air temperature, relative humidity, and surface wind speed collected from three meteorological stations in the area. The PBL heights were experimentally determined using radiosonde data provided by a station located in the center of the study area. Variations in PBL heights were estimated using linear regression methods and minimization of statistical errors for the bulk Richardson number, as well as analysis of vertical temperature and wind profiles. The WRF model reliably reproduced the daily values and diurnal cycles of temperature, relative humidity, and wind speed within the PBL and accounted for the influence of topography and sea breezes. Horizontal heat advection dominates the upwelling of air masses when sea breezes are active. The onshore wind direction starts to change from east to northwest, implying a decay in the land breeze regime. All schemes overestimate the mixing height and tend to underestimate surface air temperature values at night. All show wetter conditions and underestimate wind speed. Although the non-local Yonsei University (YSU) scheme shows the best performance, it also shows the largest sources of errors when determining the behavior of the surface layer during stable conditions. Relative humidity and wind speed estimates provided by the local Mellor‐Yamada‐Nakanishi‐Niino (MYNN) scheme were closer to those recorded at the meteorological stations.     

Evidence of Organized Large Eddies by Ground-Based Doppler Lidar, Sonic Anemometer and Sodar

In an experiment investigating the planetary boundary layer (PBL) wind and temperature fields, and PBL inversion height recorded by various instruments, the results reveal the presence of organized large eddies (OLE) or rolls. The measurements by lidars, anemometers, soundings and sodar gave an overview of the characteristics of the rolls and sources of energy production that maintain them. The experimental results obtained on two consecutive days are compared to model outputs. The agreement is excellent, showing that thermal stratification and wind shear are important factors in the structure and dynamics of OLE. A heterodyne Doppler lidar (HDL) is shown to be a useful tool in the study of OLE.     

春季黄海海雾WRF参数化方案敏感性研究

利用2005—2011年10次春季黄海海雾个例开展WRF模式参数化方案敏感性研究。结果表明:边界层方案对WRF模式雾区模拟结果起决定作用,而微物理方案影响较小,它主要影响海雾浓度和高度。边界层与微物理方案的最佳组合为YSU与Lin方案,最差为Mellor-Yamada与WSM5方案;Mellor-Yamada和QNSE方案模拟的近海面湍流过强,导致边界层过高,不利于海雾的发展与维持;而MYNN与YSU方案刻画的湍流强度与边界层高度合适,有利于海雾发展与维持。MYNN方案虽与YSU方案相当,但在大多数海雾个例中,后者明显优于前者,而在有些个例中却刚好相反。因此对于某一具体海雾个例而言,所用边界层方案仍需在它们之中选择最优者。这些信息可为黄海海雾WRF模式边界层与微物理方案的选择与改进提供参考。     

Evidence of a Mixed-Layer Dynamics Contribution to the Entrainment Process

The internal thermal boundary layer developing over the Mediterranean during a cold-air outbreak associated with a Tramontane event has been studied by means of airborne lidar, in situ sensors, and a modelling approach that consisted of nesting the University of Washington (UW) planetary boundary-layer (PBL) model in an advective zero-order jump model. This approach bypasses some of the deficiencies associated with each model: the absence of the dynamics in the mixed layer for the zero-order jump model and the lack of an inversion at the PBL top for the UW PBL model. Particular attention is given to the parameterization of the entrainment flux at the PBL top. Values of the entrainment closure parameter derived with the model when matching PBL structure observations are much lower than those derived with standard zero-order jump models. They also are in good agreement with values measured in different meteorological situations by other studies. This improvement is a result of the introduction of turbulent kinetic energy production in the mixed layer.     

Marine Boundary Layer And Turbulent Fluxes Over The Baltic Sea: Measurements And Modelling

Two weeks of measurements of the boundary-layer height over a small island (Christiansø) in the Baltic Sea are discussed. The meteorological conditions are characterised by positive heat flux over the sea. The boundary-layer height was simulated with two models, a simple applied high-resolution (2 km × 2 km) model, and the operational numerical weather prediction model HIRLAM (grid resolution of 22.5 km × 22.5 km). For southwesterly winds it was foundthat a relatively large island (Bornholm) lying 20-km upwind of the measuring site influences the boundary-layer height. In this situation the high-resolution simple applied model reproduces the characteristics of the boundary-layer height over the measuring site. Richardson-number based methods using data from simulations with the HIRLAM model fail, most likely because the island and the water fetch to the measuring site are about the size of the grid resolution of the HIRLAM model and therefore poorly resolved. For northerly winds, the water fetch to the measuring site is about 100 km. Both models reproduce the characteristics of the height of the marine boundary layer. This suggests that the HIRLAM model adequately resolves a water fetch of 100 km with respect to predictions of the height of the marine boundary layer.     

A mixed-layer model for regional evaporation

A simple mixed-layer model is developed to describe evaporation into a convective planetary boundary layer (PBL). The model comprises volume budget equations for temperature and humidity, equations to describe transport through the surface layer which is treated as part of the lower boundary, and equations to describe entrainment at the top of the PBL. The ground surface is modelled as a canopy resistance. The model was integrated with canopy resistance, surface-layer resistance and available energy, (R _n – G), input as given functions of time, and the simulated PBL was allowed to grow into an atmosphere with known temperature and humidity profiles.Two variants of the mixed-layer model were tested using data from the KNMI tower site at Cabauw in the Netherlands. These variants differed only in the formulation of entrainment: one used a formulation developed by Driedonks (1982) while the other was a simpler formulation. Simulated evaporation agreed very well with observations irrespective of which entrainment formulation was used, despite discrepancies between simulated and observed PBL height growth which were sometimes quite large for the simpler formulation. Sensitivity analysis of the model confirms that good PBL height-growth predictions are not always a prerequisite for good evaporation predictions.     

Low Tropospheric Layers Over Reunion Island in Lidar-Derived Observations and a High-Resolution Model

In November and December 2008, ground-based mobile lidar (GBML) measurements were carried out on Reunion Island (Indian Ocean, $21^{\circ }07^{\prime }\hbox {S}, 55^{\circ }32^{\prime }\hbox {E}$ , 700 km east of Madagascar) with an ultraviolet (355 nm) aerosol-backscatter lidar. Complex substructures were identified within the planetary boundary layer (PBL). A 500-m-resolution non-hydrostatic model was used to simulate the dynamics of the lower troposphere for two observation periods characteristic of the two main weather regimes in this season: the “trade-wind” regime and the “breeze” regime. The model captured the observed structures with a high degree of realism compared to the GBML. A complete diurnal cycle of the PBL along the south coast of the island during a “trade-wind” day was observed and simulated. The PBL depth was found to be anti-correlated with the wind speed. The model showed that the PBL along the coast behaved as a shallow-water flow in hydraulic theory. As the flow accelerated in response to lateral constriction, conversion of potential into kinetic energy forced the PBL top downwards. This favoured rapid transport of concentrated surface emissions within the contracted surface layer, with a possible impact on air quality. GBML observations were also conducted during the early morning of a “breeze” day on the western slope of the Maïdo mountain (2,200 m), at the top of which a new atmospheric observatory has been in operation since 2012. Both model and GBML revealed two superposed layers. The upper layer, higher than approximately 1,600 m above mean sea level, corresponded to free tropospheric air driven by the trade winds. Below, westerly counterflow advection of humid marine air occurred as a result of wake vortices in the lee of the island. The model suggests that free-tropospheric conditions prevail at the observatory from the second half of the night to mid-morning.     

辽宁绥中地区海陆风的多普勒测风激光雷达观测与特征提取

2019年3月,利用相干多普勒测风激光雷达首次在辽东湾西部绥中地区进行了风廓线测量试验。根据研究区域海岸线走向采用风向的十六分位法定义局地海风和陆风,分析和提取海陆风特征验证了多普勒测风激光雷达在春季季风间断期间观测海陆风的可行性,并计算和分析了大气边界层湍流能量的变化以及回流水平变化等特性。结果表明：1）绥中地区春季存在明显的海陆风环流特征,测风激光雷达观测海陆风出现的时间与地面自动气象站观测的数据较为一致,符合海陆风日的定义。2）海陆风日发生时,水平局地回流指数（RF）较小,1.2 km以下的RF值小于0.5,使得污染物循环累积,较易形成雾霾天气;但是海风时大气边界层的高度可达1 km以上,有利于低层大气污染物向高层扩散,减轻低层大气污染。研究结果为该地边界层参数化方案的设计和污染的防治提供了参考依据。     

Estimate of Boundary-Layer Depth Over Beijing,China, Using Doppler Lidar Data During SURF-2015

Planetary boundary-layer (PBL) structure was investigated using observations from a Doppler lidar and the 325-m Institute of Atmospheric Physics (IAP) meteorological tower in the centre of Beijing during the summer 2015 Study of Urban-impacts on Rainfall and Fog/haze (SURF-2015) field campaign. Using six fair-weather days of lidar and tower data under clear to cloudy skies, we evaluate the ability of the Doppler lidar to probe the urban boundary-layer structure, and then propose a composite method for estimating the diurnal cycle of the PBL depth using the Doppler lidar. For the convective boundary layer (CBL), a threshold method using vertical velocity variance \((\sigma _w^2 >0.1\,\hbox {m}^{2}\hbox {s}^{-2})\) is used, since it provides more reliable CBL depths than a conventional maximum wind-shear method. The nocturnal boundary-layer (NBL) depth is defined as the height at which \(\sigma _w^2\) decreases to 10 % of its near-surface maximum minus a background variance. The PBL depths determined by combining these methods have average values ranging from \(\approx \)270 to \(\approx \)1500 m for the six days, with the greatest maximum depths associated with clear skies. Release of stored and anthropogenic heat contributes to the maintenance of turbulence until late evening, keeping the NBL near-neutral and deeper at night than would be expected over a natural surface. The NBL typically becomes more shallow with time, but grows in the presence of low-level nocturnal jets. While current results are promising, data over a broader range of conditions are needed to fully develop our PBL-depth algorithms.     

Numerical simulation for the effects of PBL and the surface on pollutant concentrations

On the basis of theoretical and experimental results of study of planetary boundary layer (PBL), the physical parameters describing the structure of PBL are calculated by using the data obtained from a meteor-ological tower and the effects of PBL and the surface on pollutant concentrations are numerically simulated with a time-dependent two-dimensional advection and diffusion equation.It is shown that the diurnal variation of PBL results in that of concentration. The height of mixing layer is an important factor to determine the ground-level concentration. As for an elevated point source, the height of mixing layer, growing from lower to higher than the releasing height is a necessary condition for the phenomenon of fumigation. It is also shown that the surface may be considered as a boundary with perfect reflection when Vd ≤ 0.001 m s-1, but has an important effect on concentration and must be carefully dealt with when Vd ≤ 0.01 m s-1.     

Retrieval of the single scattering albedo of Asian dust mixed with pollutants using lidar observations

The vertical distribution of single scattering albedos （SSAs） of Asian dust mixed with pollutants was derived using the multi-wavelength Raman lidar observation system at Gwangju （35.10°N,126.53°E）.Vertical profiles of both backscatter and extinction coefficients for dust and non-dust aerosols were extracted from a mixed Asian dust plume using the depolarization ratio from lidar observations.Vertical profiles of backscatter and extinction coefficients of non-dust particles were input into an inversion algorithm to retrieve the SSAs of non-dust aerosols.Atmospheric aerosol layers at different heights had different light-absorbing characteristics.The SSAs of non-dust particles at each height varied with aerosol type,which was either urban/industrial pollutants from China transported over long distances at high altitude,or regional/local pollutants from the Korean peninsula.Taking advantage of independent profiles of SSAs of non-dust particles,vertical profiles of SSAs of Asian dust mixed with pollutants were estimated for the first time,with a new approach suggested in this study using an empirical determination of the SSA of pure dust.The SSAs of the Asian dust-pollutants mixture within the planetary boundary layer （PBL） were in the range 0.88-0.91,while the values above the PBL were in the range 0.76-0.87,with a very low mean value of 0.76 ± 0.05.The total mixed dust plume SSAs in each aerosol layer were integrated over height for comparison with results from the Aerosol Robotics Network （AERONET） measurements.Values of SSA retrieved from lidar observations of 0.92 ± 0.01 were in good agreement with the results from AERONET measurements.     

The wind in the baroclinic boundary layer with three sublayers incorporating the weak non-linear effect

In considering the weak non-linear effect, and using the small parameter expansion method, the analyt-ical expressions of the wind distribution within PBL (planetary boundary layer) and the vertical velocity at the top of the PBL are obtained when the PBL is divided into three layers and different eddy transfer coefficients K are adopted for the three layers. The conditions of barotropy and neutrality for the PBL are extended to that of baroclinity and non-neutral stratification. An example of a steady circular vortex is used to display the characteristics of the horizontal wind within the PBL and the vertical velocity at the top of the PBL. Some new results have been obtained, indicating that the magnitude of the speed in the lower height calculated by the present model is larger than that by the model in which k is a constant within the whole boundary layer, for example, in the classical Ekman boundary layer model and the model by Wu (1984). The angle between the wind at the top of the PBL and the wind near the surface calculated by the present model is less than that calculated by the single K model. These results are in agreement with the observations.     

Wave-Modified Flux and Plume Dispersion in the Stable Boundary Layer

The effects of a pressure jump and a following internal gravity wave on turbulence and plume diffusion in the stable planetary boundary layer are examined. The pressure jump was accompanied by a sudden increase in turbulence and plume dispersion. The effects of wave perturbations on turbulence statistics are analysed by calculating fluxes and variances with and without the wave signal for averaging times ranging from 1 to 30 min. The wave signals are obtained using a band-pass filter. It is shown that second-order turbulence quantities calculated without first subtracting the wave perturbations from the time are greater than those calculated when the wave signal is separated from the turbulence. Estimates of the vertical dispersion of an elevated tracer plume in the stable boundary layer are made using an elastic backscatter lidar. Plume dispersion observed 25 m downwind of the source increases rapidly with the arrival of the flow disturbances. Measured plume dispersion and plume centreline height correlate with the standard deviation of the vertical velocity but not with the wave signal.     

Development of a Single-Column Model in RegCM4 and Its Preliminary Application for Evaluating PBL Schemes in Simulating the Dry Convection Boundary Layer

A single-column model (SCM) is developed in the regional climate model RegCM4. The evolution of a dry convection boundary layer (DCBL) is used to evaluate this SCM. Moreover, four planetary boundary layer (PBL) schemes, namely the Holtslag-Boville scheme (HB), Yonsei University scheme (YSU), and two University of Washington schemes (UW01, Grenier-Bretherton-McCaa scheme and UW09, Bretherton-Park scheme), are compared by using the SCM approach. A large-eddy simulation (LES) of the DCBL is performed as a benchmark to examine how well a PBL parameterization scheme reproduces the LES results, and several diagnostic outputs are compared to evaluate the schemes. The results show that the SCM is proper constructed. In general, with the DCBL case, the YSU scheme performs best for reproducing the LES results, which include well-mixed features and vertical sensible heat fluxes; the simulated wind speed, turbulent kinetic energy, entrainment flux, and height of the entrainment zone are all underestimated in the UW09; the UW01 has all those biases of the UW09 but larger, and the simulated potential temperature is not well mixed; the HB is the least skillful scheme, by which the PBL height, entrainment flux, height of the entrainment zone, and the vertical gradients within the mixed layer are all overestimated, and a inversion layer near the top of the surface layer is wrongly simulated.Although more cases and further testing are required, these simulations show encouraging results towards the use of this SCM framework for evaluating the simulated physical processes by the RegCM4.     

BJ-RUC系统对北京夏季边界层的预报性能评估

以北京市观象台2010年8月、2011年8月每日3次 (08:00, 14:00, 20:00,北京时,下同) L波段探空秒间隔数据为实况,对BJ-RUC系统 (rapid updated cycle system for the Beijing area) 分析和预报边界层性能进行了初步评估。结果表明:BJ-RUC系统对北京地区夏季白天边界层的细致特征具有较好的预报能力,但也存在明显的系统性误差。08:00边界层偏冷; 14:00和20:00 1 km以下的边界层则显著偏暖, 边界层内明显偏湿。整体上模式对边界层内温度、湿度的预报误差均高于自由大气。该系统对北京地区边界层内早晨 (08:00) 从夜间山风向白天谷风环流过渡、午后 (14:00) 到日落后 (20:00)1500 m以下盛行西南偏南气流的日变化特征具有较强的预报能力。系统预报的14:00边界层顶高度与评估时段内实际对流边界层高度的变化趋势一致。但预报的对流边界层顶偏高,这与BJ-RUC系统采用YSU边界层参数化方案的垂直混合更强有关。