夜间城市大气边界层和气溶胶的相互作用

本文利用能量闭合的二维非线性、非定常模式,结合地面热量平衡方程,研究了夜间城市边界层和气溶胶的相互作用问题。结果表明：气溶胶在夜间对大气低层起保温作用,对大气上层起冷却作用;使大气低层稳定度减小,上层稳定度增加;此外,气溶胶还能够削弱贴地逆温强度。在正常城市气溶胶污染情况下,气溶胶对城市热岛强度影响不大,但可使城市热岛环流稍有增加。夜间城市边界层对气溶胶的反馈作用使大气下层气溶胶浓度减小,上层气溶胶浓度增加。上述部分结论得到了在天津取得的城市热岛观测资料的直接验证。     

Radiative Processes in the Stable Boundary Layer: Part II. The Development of the Nocturnal Boundary Layer

The interaction between radiation and turbulence in the stable boundary layer over land is explored using an idealized model, with a focus on the surface layer after the evening transition. It is shown that finer vertical resolution is required in transitional boundary layers than in developed ones. In very light winds radiative cooling determines the temperature profile, even if similarity functions without a critical Richardson number are used; standard surface similarity theory applied over thick layers then yields poor forecasts of near-surface air temperatures. These points are illustrated with field data. Simulations of the developing nocturnal boundary layer are used to explore the wider role of radiation. Comparatively, radiation is less significant within the developed stable boundary layer than during the transition; although, as previous studies have found, it remains important towards the top of the stable layer and in the residual layer. Near the ground, reducing the surface emissivity below one is found to yield modest relative radiative warming rather than intense cooling, which reduces the potential importance of radiation in the developed surface layer. The profile of the radiative heating rate may be strongly dependent on other processes, leading to quite varied behaviour.     

城市大气边界层的数值研究

本文利用能量闭合的二维非线性、非定常模式,结合地面热量平衡方程,研究了城市热岛和城市热岛环流以及其它相关的城市边界层参量的发展演变规律;同时探讨了在城市污染面源的情况下,污染物的扩散、输送问题.模拟结果得到了天津市城市热岛观测资料的证实.     

Large-Eddy Simulation of air Pollution Dispersion in the Nocturnal Cloud-Topped Atmospheric Boundary Layer

Effects of stratocumulus clouds on the dispersion of contaminants are studied in the nocturnal atmospheric boundary layer. The study is based on a large-eddy simulation (LES) model with a bulk parametrization of clouds. Computations include Lagrangian calculations of atmospheric dispersion of a passive tracer released from point sources at various heights above the ground. The results obtained show that the vertical diffusion is non-Gaussian and depends on the location of a source in the boundary layer.     

北京城区大气边界层高度的反演研究

采用北京市国家大气探测试验基地2017年5月1日至2019年8月31日的微波辐射计观测结果,通过位温气块法反演了该地区的大气边界层高度,统计其日变化和月际变化特征,并与相应的探空反演结果进行对比.结果表明:日间大气边界层高度的变化特征与日照时长对应关系很好.除热力作用外,边界层高度受地面风速的影响较大,近地面的较强大风...     

Observations of Atmospheric Solitary Waves in the Urban Boundary Layer

The simultaneous operation of a three-axis Doppler sodar system in the centralurban area of Rome and two similar systems in the suburban area, forming atriangle about 20 km on each side, provided evidence of solitary-type wavesin the urban boundary layer. Three events, each lasting from a few minutes toabout 30 min, and ranging in depth from the minimum range of the sodar (39 m) to over 500 m, are reported here. Two events were recognizable onall three sodar records while the third event could be observed at the urbanlocation only. Time-height acoustic echo intensity records showed no-echoregions within the wave indicating transport of trapped recirculating air.This is typical of large amplitude solitary waves. The time series plots ofsodar-derived vertical wind velocity revealed a maximum peak-to-peakvariation of about 5 m s^-1 during periods of wave-associated disturbance.The vertical velocity is found to increase with height up to the top of the closedcirculation within the wave and decreases further above. The normalisedamplitude-wavelength relationship for the two events indicates that theobserved waves are close to a strongly nonlinear regime.     

宜昌冬季两次降水过程大气边界层的观测分析

利用宜昌2007年12月10-25日的加密观测资料,分析了两次低值系统经过宜昌时大气边界层的温湿风廓线结构及其日变化特征。结果表明：位温廓线具有明显的日变化特征,对流边界层在白天出现和发展,其高度可达600m,而稳定边界层在夜间出现和发展,其高度可达300m,降水会抑制对流边界层和稳定边界层的发展;湿度廓线结构及其日变化与对流边界层的发展有关,总体上湿度随高度减小,贴近地面的薄层湿度随高度减小较快,而混合层内湿度随高度变化较小,出现降水时,近地层的湿度有明显增加,大气边界层内湿度随高度快速平稳减小;风速廓线结构比较复杂,总体上风速随高度增大,在大气边界层低层有时会出现一个风速极大值,风速廓线没有明显的日变化特征,大气边界层内风向变化较大,但以偏东风为主。     

Estimation Of Atmospheric Water Vapour Flux Profiles In The Nocturnal Unstable Urban Boundary Layer With Doppler Sodar And Raman Lidar

The vertical wind profiles determined by Doppler sodar and the water vapourmixing ratio profiles obtained by Raman lidar are used to estimate the atmosphericwater vapour flux profiles in the nocturnal urban boundary layer under unstableconditions. The experiment was conducted for several nights in the central areaof Rome under a variety of moisture conditions and different urban boundary-layerflow regimes. Despite some scatter in the profiles, the latent heat flux is found tobe positive throughout the depth of the nocturnal urban boundary-layer. Thelayer-averaged flux shows a variation between -4 to +40 W m^-2, whileindividual values of flux in excess of +150 W m^-2 pertain to a case offree convection during cold air advection caused by the sea breeze. The qualityof flux estimates is found to be highly limited by the low sampling rates employedin the experiment resulting in errors to the order of 60%. Therefore, the results mustbe viewed as estimates rather than precise measurements. The skewness profiles ofthe turbulent fluctuations of vertical velocity and water vapour mixing ratio are alsopositive.     

北京冬季大风过程大气边界层特征分析

利用北京中国科学院大气物理研究所325 m气象观测塔的气象梯度资料和湍流资料,分析了2014年11月29日至12月5日北京两次大风过程中气象要素和湍流输送特征的变化。第一次大风过程的强度和持续时间均高于第二次大风过程。强烈的风速垂直切变主要集中在距地面100 m高度范围内,最强风速垂直切变达到0.31 s~(-1)。大风过程中,阵风系数呈现随高度减小的趋势,越接近地面,阵风系数愈大。阵风强度的变化与阵风系数相似,100 m以下高度时,阵风强度随高度增大而减小。大风过程自上而下改变边界层结构,平均动能、湍流动能和摩擦速度最先从上层(280 m)发生变化且迅速增加。近地层由于风速垂直梯度的显著差异,近地层垂直方向的湍流强度最大。大风时各功率谱在低频区(0.01 s~(-1))达到峰值,大风过后各高度的能量都有所下降。     

南京地区大气颗粒物影响近地面臭氧的个例研究

通过对2008年4月2～7日南京地区地面气象观测数据以及两个站点空气质量(O₃、NO_x、PM10)监测资料的分析, 发现O₃和PM10之间存在一定程度的反相关。利用一个光化学箱模式对该个例中大气颗粒物影响近地面臭氧的过程进行模拟, 结果发现大气颗粒物浓度的升高使得气溶胶光学厚度增加20％～40%, 导致NO₂和O₃近地面光解率下降20%～30%, OH和HO₂自由基浓度分别减少20%～50%, 造成O₃净生成率下降30%～40%。研究表明, 颗粒物对光化学过程的抑制造成了大气氧化能力的降低, 是近地面臭氧浓度减少的可能原因。     

Large-Eddy Simulation Of A Nocturnal Stratocumulus-Topped Marine Atmospheric Boundary Layer: An Uncertainty Analysis

A large-eddy simulation (LES) model has been used to study a nocturnalstratocumulus-topped marine atmospheric boundary layer. The main objectivesof our study have been first to investigate the statistical significance of LES-derived data products. Second, to test the sensitivity of our LES results with respect to the representation of subgrid-scale mixing and microphysical processes, and third to evaluate and to quantify the parametric uncertainty arising from the incomplete knowledge of the environmental parameters that are required to specify the initial and boundary conditions of a particular case study. Model simulations were compared with observations obtained in solid stratocumulus during the third flight of the first 'Lagrangian' experiment of the Atlantic Stratocumulus Transition Experiment (ASTEX). Based on these simulations the following conclusions could be drawn. Resolution(50 × 50 × 25 m³) and domain size (3.2 × 3.2 × 1.5 km³) of the LES calculations were adequate from a numerical point of view to represent the essential features of the stratocumulus-topped boundary layer. However, the ensemble runs performed in our study to investigate the statistical significance of LES-derived data products demonstrate that the area-time averaging procedure for the second-order moments produces only a low degree of statistical reliability in the model results. This illustratesthe necessity of having LES model results that are not only of adequate resolution but also of sufficiently large domain. The impact of different subgrid schemes was small, but the primary effects of drizzle were found to influence the boundary-layer structure in a climatologically significant way. The parametric uncertainty analysis revealed that the largest contribution to the variance of the LES-derived data products is due to theuncertainties in the cloud-top jump of total water mixing ratio and the net radiative forcing. The differences between the model and measurements for most of the simulated quantities were within the modelling uncertainties, but the calculated precipitation rate was found to differ significantly from that derived in the observations.     

夜间边界层的一维数值模拟

本文用一维数值模式,模拟了平坦地形情况下小风、晴空夜间边界层内风、温、湿的演变,计算结果与观测事实有较好的一致性。对有高云或少量中、低云存在时的计算结果分析表明,模式对风速的模拟仍然较好,对温度的模拟比实测值可偏低0.3—1℃,对相对湿度的模拟比实测值可偏高5—10％在符合模式使用范围的情况下,本模式可作为夜间边界层预报的一种方法参考使用。     

Directional Shear in the Nocturnal Atmospheric Surface Layer

Stratified nocturnal flow above and within a small valley of approximately 12-m depth and a few hundred metres width is examined as a case study, based on a network of 20 sonic anemometers and a central 20-m tower with eight levels of sonic anemometers. Several regimes of stratified flow over gentle topography are conceptually defined for organizing the data analysis and comparing with the existing literature. In our case study, a marginal cold pool forms within the shallow valley in the early evening but yields to larger ambient wind speeds after a few hours, corresponding to stratified terrain-following flow where the flow outside the valley descends to the valley floor. The terrain-following flow lasts about 10 h and then undergoes transition to an intermittent marginal cold pool towards the end of the night when the larger-scale flow collapses. During this 10-h period, the stratified terrain-following flow is characterized by a three-layer structure, consisting of a thin surface boundary layer of a few metres depth on the valley floor, a deeper boundary layer corresponding to the larger-scale flow, and an intermediate transition layer with significant wind-directional shear and possible advection of lee turbulence that is generated even for the gentle topography of our study. The flow in the valley is often modulated by oscillations with a typical period of 10 min. Cold events with smaller turbulent intensity and duration of tens of minutes move through the observational domain throughout the terrain-following period. One of these events is examined in detail.     

Pollutant Plume Dispersion in the Atmospheric Boundary Layer over Idealized Urban Roughness

The Gaussian model of plume dispersion is commonly used for pollutant concentration estimates. However, its major parameters, dispersion coefficients, barely account for terrain configuration and surface roughness. Large-scale roughness elements (e.g. buildings in urban areas) can substantially modify the ground features together with the pollutant transport in the atmospheric boundary layer over urban roughness (also known as the urban boundary layer, UBL). This study is thus conceived to investigate how urban roughness affects the flow structure and vertical dispersion coefficient in the UBL. Large-eddy simulation (LES) is carried out to examine the plume dispersion from a ground-level pollutant (area) source over idealized street canyons for cross flows in neutral stratification. A range of building-height-to-street-width (aspect) ratios, covering the regimes of skimming flow, wake interference, and isolated roughness, is employed to control the surface roughness. Apart from the widely used aerodynamic resistance or roughness function, the friction factor is another suitable parameter that measures the drag imposed by urban roughness quantitatively. Previous results from laboratory experiments and mathematical modelling also support the aforementioned approach for both two- and three-dimensional roughness elements. Comparing the UBL plume behaviour, the LES results show that the pollutant dispersion strongly depends on the friction factor. Empirical studies reveal that the vertical dispersion coefficient increases with increasing friction factor in the skimming flow regime (lower resistance) but is more uniform in the regimes of wake interference and isolated roughness (higher resistance). Hence, it is proposed that the friction factor and flow regimes could be adopted concurrently for pollutant concentration estimate in the UBL over urban street canyons of different roughness.     

Bulk Mixing and Decoupling of the Nocturnal Stable Boundary Layer Characterized Using a Ubiquitous Natural Tracer

Vertical mixing of the nocturnal stable boundary layer (SBL) over a complex land surface is investigated for a range of stabilities, using a decoupling index ( $0 < D_{rb} < 1$ ) based on the 2–50 m bulk gradient of the ubiquitous natural trace gas radon-222. The relationship between $D_{rb}$ and the bulk Richardson number ( $R_{ib}$ ) exhibits three broad regions: (1) a well-mixed region ( $D_{rb} \approx 0.05$ ) in weakly stable conditions ( $R_{ib} < 0.03$ ); (2) a steeply increasing region ( $0.05 < D_{rb} < 0.9$ ) for “transitional” stabilities ( $0.03 < R_{ib} < 1$ ); and (3) a decoupled region ( $D_{rb} \approx 0.9$ –1.0) in very stable conditions ( $R_{ib} > 1$ ). $D_{rb}$ exhibits a large variability within individual $R_{ib}$ bins, however, due to a range of competing processes influencing bulk mixing under different conditions. To explore these processes in $R_{ib}$ – $D_{rb}$ space, we perform a bivariate analysis of the bulk thermodynamic gradients, various indicators of external influences, and key turbulence quantities at 10 and 50 m. Strong and consistent patterns are found, and five distinct regions in $R_{ib}$ – $D_{rb}$ space are identified and associated with archetypal stable boundary-layer regimes. Results demonstrate that the introduction of a scalar decoupling index yields valuable information about turbulent mixing in the SBL that cannot be gained directly from a single bulk thermodynamic stability parameter. A significant part of the high variability observed in turbulence statistics during very stable conditions is attributable to changes in the degree of decoupling of the SBL from the residual layer above. When examined in $R_{ib}$ – $D_{rb}$ space, it is seen that very different turbulence regimes can occur for the same value of $R_{ib}$ , depending on the particular combination of values for the bulk temperature gradient and wind shear, together with external factors. Extremely low turbulent variances and fluxes are found at 50 m height when $R_{ib} > 1$ and $D_{rb} \approx 1$ (fully decoupled). These “quiescent” cases tend to occur when geostrophic forcing is very weak and subsidence is present, but are not associated with the largest bulk temperature gradients. Humidity and net radiation data indicate the presence of low cloud, patchy fog or dew, any of which may aid decoupling in these cases by preventing temperature gradients from increasing sufficiently to favour gravity wave activity. The largest temperature gradients in our dataset are actually associated with smaller values of the decoupling index ( $D_{rb} < 0.7$ ), indicating the presence of mixing. Strong evidence is seen from enhanced turbulence levels, fluxes and submeso activity at 50 m, as well as high temperature variances and heat flux intermittencies at 10 m, suggesting this region of the $R_{ib}$ – $D_{rb}$ distribution can be identified as a top-down mixing regime. This may indicate an important role for gravity waves and other wave-like phenomena in providing the energy required for sporadic mixing at this complex terrain site.     

Different Turbulent Regimes and Vertical Turbulence Structures of the Urban Nocturnal Stable Boundary Layer

Turbulence in the nocturnal boundary layer(NBL) is still not well characterized, especially over complex underlying surfaces. Herein, gradient tower data and eddy covariance data collected by the Beijing 325-m tower were used to better understand the differentiating characteristics of turbulence regimes and vertical turbulence structure of urban the NBL. As for heights above the urban canopy layer(UCL), the relationship between turbulence velocity scale(V_TKE) and wind speed(V) was con...     

Observations of Wind-Direction Variability in the Nocturnal Boundary Layer

Large sudden wind-direction shifts and submeso variability under nocturnal conditions are examined using a micrometeorological network of stations in north-western Victoria, Australia. The network was located in an area with mostly homogeneous and flat terrain. We have investigated the main characteristics of the horizontal propagation of events causing the wind-direction shift and not addressed in previous studies. The submeso motions at the study site exhibit behaviour typical of flat terrain, such as the lower relative mesovelocity scale and smaller cross-wind variances than that for complex terrain. The distribution of wind-direction shifts shows that there is a small but persistent preference for counter-clockwise rotation, occurring for 55% of the time. Large wind-direction shifts tend to be associated with a sharp decrease in air temperature (74% of the time), which is associated with rising motion of cold air, followed by an increase in turbulent mixing. The horizontal propagation of events was analyzed using the cross-correlation function method. There is no preferred mean wind direction associated with the events nor is there any relationship between the mean wind and propagation directions. The latter indicates that the events are most likely not local flow perturbations advected by the mean flow but are rather features of generally unknown origin. This needs to be taken into account when developing parametrizations of the stable boundary layer in numerical models.     

The Morning Transition of the Nocturnal Boundary Layer

Six years of observations from a surface instrument site have been analysed to determine timings and factors influencing developmental changes in the near-surface wind and turbulent heat fluxes during the morning heating of the atmospheric boundary layer. A simple relationship has been found between near-surface wind speed and screen temperature, together with a predictive equation for the morning transition air temperature. Profile measurements from a probe mounted on a tethered balloon have beenused to supplement the surface data and study the processes underlying these surface relationships. The results have confirmed earlier work and have shown that both before and immediately after morning transition, almost all heating in the surface layer is due to turbulent diffusion from above. In order to explain the mechanisms involved in the relationships, a simple finite difference model has been run and validated against the profile data. The model predictions are compared with observations during both the morning and evening and the differences related to the different temperature profiles. Numerical forecasting rules for the surface wind speed and transition temperature are derived from the results.     

大气边界层气象学研究综述

文中回顾了大气边界层气象学的发展历史，总结了目前大气边界层气象学的主要进展，并指出国内外在未来大气边界层气象学研究方面面临的一些主要科学问题，以及对未来大气边界层气象学的发展方向提出若干建议，同时还指出了大气边界层气象学在思想上和方法上应该注意的一些相关问题。     

Impacts of Aerosol Shortwave Radiation Absorption on the Dynamics of an Idealized Convective Atmospheric Boundary Layer

We investigated the impact of aerosol heat absorption on convective atmospheric boundary-layer (CBL) dynamics. Numerical experiments using a large-eddy simulation model enabled us to study the changes in the structure of a dry and shearless CBL in depth-equilibrium for different vertical profiles of aerosol heating rates. Our results indicated that aerosol heat absorption decreased the depth of the CBL due to a combination of factors: (i) surface shadowing, reducing the sensible heat flux at the surface and, (ii) the development of a deeper inversion layer, stabilizing the upper CBL depending on the vertical aerosol distribution. Steady-state analytical solutions for CBL depth and potential temperature jump, derived using zero-order mixed-layer theory, agreed well with the large-eddy simulations. An analysis of the entrainment zone heat budget showed that, although the entrainment flux was controlled by the reduction in surface flux, the entrainment zone became deeper and less stably stratified. Therefore, the vertical profile of the aerosol heating rate promoted changes in both the structure and evolution of the CBL. More specifically, when absorbing aerosols were present only at the top of the CBL, we found that stratification at lower levels was the mechanism responsible for a reduction in the vertical velocity and a steeper decay of the turbulent kinetic energy throughout the CBL. The increase in the depth of the inversion layer also modified the potential temperature variance. When aerosols were present we observed that the potential temperature variance became significant already around $0.7z_i$ (where $z_i$ is the CBL height) but less intense at the entrainment zone due to the smoother potential temperature vertical gradient.