边界层方案对华北低层O3垂直分布模拟的影响

利用WRF-Chem模式,采用3种边界层参数化方案 (YSU, MYJ和ACM2),针对1个晴空、静稳日 (2013年8月26日20:00—27日20:00(北京时)) 进行模拟,着重分析不同边界层参数化方案对夜间残留层形成及日出前后O₃浓度垂直分布形式的模拟效果,并与固城站地面及垂直同步观测资料进行对比。结果表明:3种边界层参数化方案均能够模拟出温度及风速的区域分布形式以及风温垂直结构的变化特征;相比之下,MYJ方案模拟的夜间边界层高度较YSU方案和ACM2方案明显偏高,该对比结果可能是导致近地面污染物浓度模拟差异的重要原因;在夜间稳定层结至日出后稳定状态打破的边界层结构演变过程中,采用YSU方案和ACM2方案模拟的温度和风速垂直扩线形式与观测结果更为接近;同样采用非局地闭合的YSU方案和同时考虑局地和非局地闭合的ACM2方案,对于边界层高度内O₃浓度垂直分布形式的模拟效果具有明显优势。     

WRF模式中QNSE方案的湍流长度尺度系数的调整试验研究

边界层参数化方案中湍流混合对数值模拟起着重要的作用,湍流混合作用的恰当描述对于温度、湿度、风场以及降水的准确模拟至关重要。我国长江中下游流域人口密集,暴雨灾害频发,很有必要寻找一种适合该地区降水模拟的边界层参数化方案。本文运用WRF（Weather Research and Forecasting）中尺度数值模式,以QNSE（Quasi-Normal Scale Elimination）边界层参数化方案为基础,将其中湍流混合长度尺度系数调整为可变参数。本文将Noh et al.（2003）提出的Prandtl公式与Janji?提出的修正湍流长度尺度系数的方法相结合,通过考虑非局地项的强迫、地表稳定度与边界层高度对湍流长度尺度系数的影响,强调大气的动力结构特征与热力结构特征对湍流混合的共同影响,从而提高QNSE边界层参数化方案在不同地理环境下的模拟能力。文章通过进行长江中下游地区的典型暴雨试验,将调整参数后的QNSE方案与原方案进行比较,重点分析调整参数后的方案与原方案对关键基本气象要素场、边界层结构特征以及降水的模拟能力,并将模拟结果与观测结果进行对比,结果表明:调整参数后的方案一定程度上改进了地表温度、边界层结构以及降水的模拟效果。进一步研究表明,调整参数后的方案主要通过改变边界层混合缓解水汽混合比、位温模拟方面的误差。     

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

边界层高度是影响大气边界层发展和空气污染程度的重要因子,是环境和气候研究的重要参数。本文利用激光雷达对北京地区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。     

Data Assimilation Strategies in the Planetary Boundary Layer

We investigate the effect of the assimilation of surface and boundary-layer mass-field observations on the planetary boundary layer (PBL) within a one-dimensional (1D) version of the non-hydrostatic Fifth-Generation Pennsylvania State University/National Center for Atmospheric Research Mesoscale Model (MM5). We focus on the vertical extent and effects of mass-field nudging within the PBL based on surface observations, and the added value of assimilating column mass observations within the PBL. Model experiments for dynamic initialization and dynamic analysis are conducted and composited for 29 May, 6 June, and 7 June 2002 during the International H₂O Project (IHOP) over the Southern Great Plains, U.S.A. Advantages are found when the data assimilation uses the innovation (the difference between the modelled value and the observed value) calculated by comparing the surface mass-field observation to the model value at the 2-m observation height rather than at the lowest model level. It is shown that this innovation can be applied throughout the model-diagnosed PBL via nudging during free-convective conditions because of the well-mixed nature of the PBL. However, in stable conditions, due to decreased vertical mixing the surface innovation may be best applied only in a shallow layer adjacent to the surface. Surface air-temperature innovations were also applied to the top soil-layer temperature to minimize disruption to the surface energy balance. In combination with the surface observations, the use of within-PBL mass-field data assimilation improves the simulated PBL structure.     

An analytical study of diurnal wind-structure variations in the boundary layer and the low-level nocturnal jet

An analytical framework is proposed for studying variations in the diurnal wind structure in the planetary boundary layer (PBL) and the evolution of the low-level nocturnal jet. A time-dependent eddy-diffusivity coefficient corresponding to solar input is proposed, and an appropriate coordinate transformation ensures that mixing height varies continuously with ground heat-flux changes. The solution exhibits the receding character of the daytime PBL as evening approaches, thereby dividing the PBL into two regimes — the one just above the ground, representing the nocturnal boundary layer, and the region above it. It is assumed that inertial oscillations (IO) are triggered in the upper layer at about the time of sunset when the reversal in the direction of ground heat flux is felt in the upper layer. Two approaches are adopted to determine the characteristic features of IO and the evolution of the nocturnal low-level jet. The first one is based on the physical principle that release of horizontal momentum due to deviation from the geostrophic wind gives rise to the IO. The solution captures all the characteristic features of the IO, such as phase shift and decreasing amplitude of the IO with increasing height. According to this analysis the IO is triggered at a level as soon as the top of the receding boundary layer leaves that level. The solution is discontinuous with respect to the vertical coordinate. In the second approach we solve an initial-value problem to determine the solution in the upper layer, assuming that at about the time of sunset there is a rapid collapse of the daytime PBL to the steady, nocturnal boundary layer. The assumption is based on the mixing-height profiles prepared from climatological data collected at Delhi. The solution for the nocturnal boundary-layer regime is then obtained as a boundary-value problem. The solutions so obtained are continuous throughout the domain of interest and exhibit the characteristic features of an IO. The analysis leads to the conditions under which a low-level nocturnal jet is produced and provides quantitative estimates of the parameters, such as length of night, latitude, mixing height at sunset and nocturnal mixing height, that are conducive to the generation of a jet. The nocturnal wind profile produced by this approach compares well both with typical atmospheric data observed at Delhi and with output from a mesoscale numerical model. There is still some uncertainty related to the time of initiation of the IO as a function of latitude.     

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

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

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边界层参数化方案的垂直混合更强有关。     

The sensitivity of k-ɛ model computations of the neutral planetary boundary layer to baroclinicity

While the importance of baroclinicity in determining the structure of the planetary boundary layer (PBL) is well recognized, the actual effect of baroclinicity on the structure is not well understood. Results based on simulations obtained using the turbulent kinetic energy-dissipation rate of turbulent kinetic energy closure model of the turbulent flow in a neutral baroclinic PBL provide additional insight into the role of baroclinicity. The baroclinic PBL is characterized by significant shear production of turbulent kinetic energy throughout the complete boundary-layer depth. The turbulent mixing length is bounded by the presence of a stable temperature inversion layer indicating that the depth of the baroclinic PBL is determined by the inversion height. Significant turbulent shear stresses exist throughout the baroclinic PBL and the air is relatively well-mixed except in the surface layer.     

Aircraft study of aerosol vertical distributions over Beijing and their optical properties

A set of 152 vertical profiles of aerosol number concentration and size distribution with diameter ranging from 0.12 to 3.0 μm observed by the airborne optical spectrometer probe in Beijing, China, between February 2005 and September 2006 is analysed and discussed. The statistic of aerosol number concentration ( N _a) reveals a high aerosol number density in this region with average surface level number concentration ( N ₀) of about 6600 cm⁻³ (0.12–3.0 μm). The average vertical profile of N _a approximately satisfies an exponential decline function with a scale height of 1419 m. The N _a profiles are influenced by the structures of planetary boundary layer (PBL) significantly and two typical types of N _a profile under different conditions of PBL are presented and parametrized in this study. The observations of aerosol size distribution show that, in most cases the aerosol size distributions are not very sensitive to altitude, with effective radii ranging from 0.16 to 0.28 μm. Comparison between aircraft-derived aerosol optical depth (AOD) and Moderate Resolution Imaging Spectroradiometer-derived AOD shows good agreement. The Mie model calculations suggest that the surface level number concentration, the PBL height and the structure of PBL can influence the AOD significantly.     

北京等地对流层低层晴天气溶胶微物理特征

1982年4—7月对北京等地对流层低层晴天气溶胶数浓度、尺度谱分布进行了飞机观测。结果表明:来自西、西北和北部气团的气溶胶数浓度相当低（平均为10-1—100 /厘米3）。尺度谱分布用dn/dD=C/DA拟合，相关显著水平α≤0.01，幂指数A值多小于3（占总数73%以上），比Junge给出的A=4值小些。鉴于直径D=0.5—2微米粒子在大气中的重要性，也按dn/dD=C′/DB进行拟合，多数测值α≤0.01，少数α≤0.05。幂指数A、B和D=0.5—1微米粒子相对数浓度f1在水平、垂直方向上是相当均匀稳定的。受城市影响时，混合层内气溶胶数浓度明显增大，尺度谱分布斜率变陡（A、B增大），特别是和“气-粒”转化有关的D=1微米左右粒子数浓度增加非常显著，比不受城市影响时，增加了1—2个量级，最大可增加3个量级。可见，城市影响是造成混合层气溶胶数浓度增大、尺度谱分布变陡、局地大气污染加重、低层能见度变坏的重要原因。     

基于微波辐射计和气溶胶激光雷达的边界层高度研究及应用

分别基于微波辐射计温湿度廓线资料的气块法、位温法和比湿法,地面气象资料的罗氏法及气溶胶激光雷达数据的梯度法,计算得出广州地区大气边界层高度,对比分析5种边界层高度结果及其与气象条件、空气质量之间的关系,结合典型大气污染过程分析边界层高度对PM_2.5、O₃浓度的影响。结果显示：(1)利用位温法、气块法、罗氏法、比湿法和梯度法计算得出广州地区平均边界层高度分别为2 207 m、1 239 m、901 m、717 m和660 m,位温法显著高估了广州地区的边界层高度;(2)利用气块法得出的混合层高度日变化能够较好地表征白天大气边界层演变特征,利用气块法和比湿法得出的白天混合层高度与近地面O₃浓度有显著的正相关关系,相关系数在0.5以上,在O₃污染防治中,应同时考虑边界层内垂直输送的影响;(3)利用梯度法得出的边界层高度在污染天气时与PM_2.5浓度的相关性较好,能较好地表现出大气污染情况,在PM_2.5污染天气过程分析中具有较好的应用价值。     

Impact of Turbulent Mixing in the Stratocumulus-Topped Boundary Layer on Numerical Weather Prediction

The impact of enhanced turbulent mixing induced by radiative cooling at the top of the stratocumulus-topped boundary layer (STBL) on numerical weather prediction is examined. An additional term involving top-down turbulent mixing via in-cloud radiative cooling is applied to the Yonsei University (YSU) planetary boundary layer (PBL) parameterization scheme using a top-down diffusivity profile and cloud-top entrainment. The modified scheme is evaluated in an advection fog case over the Yellow Sea of Korea using the Weather Research and Forecasting (WRF) model and in global medium-range forecasts using the Global/Regional Integrated Model system (GRIMs). In the fog case simulation, consideration of the additional top-down mixing parameterization in the YSU PBL simulates less formation and more rapid dispersion of the fog. As a result, the modified scheme simulates a drier and warmer boundary layer and a moister and cooler layer above the PBL. The modified algorithm also improves surface temperature prediction over the Yellow Sea accompanying early dissipation of the fog. In the global medium-range forecast experiment, the modified scheme simulates overall enhanced PBL mixing over the STBL in the tropics and subtropical ocean, showing drier and warmer regions near the surface and moister and cooler regions above the PBL, resulting in prediction of reduced low level cloud amount and increased downward shortwave radiation at the surface. The modified scheme appears to improve systematic bias in temperature and humidity in the lower troposphere compared to the control simulation.     

热带西太平洋边界层的多层结构分析

本文主要利用第一次“中美热带西太平洋海气相互作用联合调查”所获得的海洋水文、气象及大气探测资料,分析了热带西太平洋暖水区边界层大气层结特征,指出热带西太平洋边界层大气具有多层结构和多层送温现象,边界层内各层厚度及气温、位温、混合比随高度的变化趋势如下:超绝热层位于地面以上50m左右,气温、位温、混合比迅速减小,系绝热不稳定层;第一层为混合层,平均厚度为351m,位温、混合比变化不明显;混合层上又是逆温层大多为30—100m厚,位温迅速增大,混合比迅速减小,系绝对稳定层;700—1000m为浅对流云活跃层,     

The role of radiation in an inversion-capped planetary boundary layer

In an inversion-capped planetary boundary layer (PBL), the structure of the turbulent fluxes as well as the height of the inversion are determined by the interaction of turbulent mixing in the PBL, large-scale subsidence above the PBL and radiational cooling. Here the sensitivity of the inversion height and of the turbulent fluxes due to radiational processes is investigated with the aid of a three-layered model for a well mixed PBL. For an example of the Trade-Wind region, the inversion height (i.e., the difference between surface pressure and pressure at the inversion level) varies between 46 and 257 mb and the surface flux of moist static energy between 417 and 99 W m^-2, if the (mean) radiative net flux divergence for both the inversion and the well-mixed layer is changed over a reasonable range of values. None of the parameterization schemes existing in the literature is able to describe these radiational effects in an appropriate way. This is due to the fact that these parameterizations are either not or not flexibly enough linked to the thermodynamical model parameter. Therefore the demand for an adequate parameterization of the radiational influence in a well-mixed PBL under a subsidence inversion is obvious.     

大规模风电场建成后对风能资源影响的研究

考虑了千万千瓦级风电基地建成后风电机群对近地面层风速的影响,采用Frandsen研究了大规模风电场内部风速损失时所使用的方法,在内边界层已经充分发展成新边界层的区域内,对轮毂高度65m处风速Uh进行了计算。结果表明,风电场建成后研究区内,风速Uh与未建场时的65m风速U0相比变小,存在风速损失,该风速损失随着U0的增大而减小,与风电机的推力系数CT性质有关;大规模风电场建成后,Uh在3～20m.s-1范围内的平均风功率密度与未建场时U0在此范围内的平均风功率密度相比损失约为58.45%,这与建场地区建场前65m处风速值大小以及各风速值出现的概率有关。     

Statistical Variability of Dispersion in the Convective Boundary Layer: Ensembles of Simulations and Observations

A Lagrangian particle dispersion model (LPDM) driven by velocity fields from large-eddy simulations (LESs) is used to determine the mean and variability of plume dispersion in a highly convective planetary boundary layer (PBL). The total velocity of a “particle” is divided into resolved and unresolved or random (subfilter scale, SFS) velocities with the resolved component obtained from the LES and the SFS velocity from a Lagrangian stochastic model. This LPDM-LES model is used to obtain an ensemble of dispersion realizations for calculating the mean, root-mean-square (r.m.s.) deviation, and fluctuating fields of dispersion quantities. An ensemble of 30 realizations is generated for each of three source heights: surface, near-surface, and elevated. We compare the LPDM calculations with convection tank experiments and field observations to assess the realism of the results. The overall conclusion is that the LPDM-LES model produces a realistic range of dispersion realizations and statistical variability (i.e., r.m.s. deviations) that match observations in this highly convective PBL, while also matching the ensemble-mean properties. This is true for the plume height or trajectory, vertical dispersion, and the surface values of the crosswind-integrated concentration (CWIC), and their dependence on downstream distance. One exception is the crosswind dispersion for an elevated source, which is underestimated by the model. Other analyses that highlight important LPDM results include: (1) the plume meander and CWIC fluctuation intensity at the surface, (2) the applicability of a similarity theory for plume height from a surface source to only the very strong updraft plumes—not the mean height, and (3) the appropriate variation with distance of the mean surface CWIC and the lower bound of the CWIC realizations for a surface source.     

北京地区夏季边界层结构日变化的高分辨模拟对比

使用WRF中尺度数值模式, 分别选用两种不同的边界层参数化方案 (MYJ, YSU) 和3种陆面参数化方案 (SLAB, Noah, RUC), 对2004年7月1日08:00—7月4日20:00 (北京时) 北京地区夏季边界层结构进行1 km的高分辨模拟。对比分析了近地面层风场、温度场以及边界层的日变化特征, 结果发现:WRF模式基本模拟出了北京夏季边界层的日变化特征; 在边界层方案中, MYJ方案描述的边界层结构较YSU方案合理; Noah陆面模式较好地反映了城市的热岛效应; 无降水时, 风速及边界层高度对于陆面过程不敏感, 而降水发生后, 陆面过程对于边界层结构的影响增大; 各方案模拟的城区风速明显偏大, 这是因为没有充分考虑城市建筑物的阻力作用。     

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.     

塔克拉玛干沙漠腹地大气边界层参数化方案的模拟评估

沙尘起沙、沉降、传输均受到沙漠地区大气边界层条件的制约。沙漠地区观测资料匮乏,限制大气边界层模拟效果的检验和评估。利用WRFV3.7.1中尺度数值模式中5种边界层参数化方案(ACM2、BL、MYJ、MYNN2.5、YSU),模拟2014年4月塔克拉玛干沙漠大气边界层特征,并与塔中80 m塔及风廓线雷达晴朗天气下的观测资料对比分析。结果表明:5种方案均能模拟出近地面气温及地表温度,边界层高度,感热、潜热、地表热通量的变化趋势,但未能模拟出边界层风速的日变化趋势,温风湿廓线能较好的反映晴日沙漠地区边界层结构的变化特征,但未模拟出风速随高度变化趋势。沙漠地区下垫面干燥,热容量低,晴天极易形成对流不稳定边界层,非局地湍流参数化方案,ACM2方案是沙漠地区大气边界层模拟较为合理的选择。     

长江流域一次暴雨过程的模拟对行星边界层参数化的敏感性研究

用MM5对长江流域的一次暴雨进行模拟考察其对行星边界层参数化的敏感性。不同的边界层参数化表现在不同的地表通量和垂直混合设计上，本文分析了MM5中4个主要方案地表通量和垂直混合参数化方案的不同以及它们对降水强度、落区和时间的影响。研究还发现，地表通量对暴雨模拟结果的影响比垂直混合方案要大。