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.     

Simulating Australian Urban Climate in a Mesoscale Atmospheric Numerical Model

We develop an urban canopy scheme coupled to a mesoscale atmospheric numerical model and evaluate the simulated climate of an Australian city. The urban canopy scheme is based on the Town Energy Budget approach, but is modified to efficiently represent the predominately suburban component of Australian cities in regional climate simulations. Energy conservation is improved by adding a simple model of air-conditioning to prevent the urban parametrization acting as an energy sink during the Australian summer. In-canyon vegetation for suburban areas is represented by a big-leaf model, but with a largely reduced set of prognostic variables compared to previous approaches. Although we have used a recirculation/venting based parametrization of in-canyon turbulent heat fluxes that employs two canyon wall energy budgets, we avoid using a fixed canyon orientation by averaging the canyon fluxes after integrating over 180° of possible canyon orientations. The urban canopy scheme is evaluated by simulating the climate for Melbourne, Australia after coupling it to The Air Pollution Model. The combined system was found to predict a realistic climatology of air temperatures and winds when compared with observations from Environmental Protection Authority monitoring stations. The model also produced a plausible partitioning of the urban energy budget when compared to urban flux-tower studies. Overall, the urban canyon parametrization appears to have reasonable potential for studying present and predicting changes in future Australian urban climates in regional climate simulations.     

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

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

Observations of Coherent Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer

Turbulence structures of high Reynolds number flow in the near-neutral atmospheric boundary layer (ABL) are investigated based on observations at Shionomisaki and Shigaraki, Japan. A Doppler sodar measured the vertical profiles of winds in the ABL. Using the integral wavelet transform for the time series of surface wind data, the pattern of a descending high-speed structure with large vertical extent (from the surface to more than 200-m level) is depicted from the Doppler sodar data. Essentially this structure is a specific type of coherent structure that has been previously shown in experiments on turbulent boundary-layer flows. Large-scale high-speed structures in the ABL are extracted using a long time scale (240 s) for the wavelet transform. The non-dimensional interval of time between structures is evaluated as 3.0–6.2 in most cases. These structures make a large contribution to downward momentum transfer in the surface layer. Quadrant analyses of the turbulent motion measured by the sonic anemometer (20-m height) suggest that the sweep motion (high-speed downward motion) plays a substantial role in the downward momentum transfer. In general, the contribution of sweep motions to the momentum flux is nearly equal to that of ejection motions (low-speed upward motions). This contribution of sweep motions is related to the large-scale high-speed structures.     

Impacts of Mixing Processes in Nocturnal Atmospheric Boundary Layer on Urban Ozone Concentrations

A number of open questions remain regarding the role of low-level jets (LLJs) and nocturnal mixing processes in the buildup of tropospheric ozone. The prevalence of southerly winds and LLJs in the U.S. Southern Great Plains during summer makes this region an ideal site for investigating the structure of the nocturnal boundary layer and its impacts on urban air quality. Ozone $(\mathrm{O}_{3})$ and nitrogen oxide concentrations measured at regulatory monitoring sites in the Oklahoma City (OKC) area and simulations with the Weather Research and Forecasting with Chemistry (WRF/Chem) model were analyzed to show how the nocturnal LLJ moderates boundary-layer mixing processes and air quality. Datasets collected during the Joint Urban 2003 campaign, which took place in July 2003 in OKC, provided detailed information about nocturnal boundary-layer structure and dynamics. In general, ${\mathrm{O}_{3}}$ time series show the expected behavior that urban ${\mathrm{O}_{3}}$ concentrations decrease at night due to nitrogen oxide titration reactions, but elevated ${\mathrm{O}_{3}}$ concentrations and secondary ${\mathrm{O}_{3}}$ peaks are also seen quite frequently after sunset. LLJs developed on most nights during the study period and were associated with strong vertical wind shear, which affected the boundary-layer stability and structure. Near-surface ${\mathrm{O}_{3}}$ concentrations are higher during less stable nights when active mixing persists throughout the night. The WRF/Chem model results agree well with the observations and further demonstrate the role of LLJs in moderating nocturnal mixing processes and air quality. The highest nocturnal ${\mathrm{O}_{3}}$ concentrations are linked to a strong LLJ that promotes both nocturnal long-range transport and persistent downward mixing of ${\mathrm{O}_{3}}$ from the residual layer to the surface.     

Systematic Vertical Variation of Mesoscale Fluxes in the Nocturnal Boundary Layer

The vertical mesoscale flux in the nocturnal boundary layer is generally considered to be difficult to estimate because of the small mesoscale vertical velocities and the large random variation of the mesoscale fluxes. However, the mesoscale vertical flux of heat, computed from FLOSSII data, varies quasi-systematically with height, stability and time scale. Such systematic variation requires correction for sonic misalignment and averaging over a large quantity of data. The relation of the mesoscale heat flux to the vertical structure of the nocturnal boundary layer is examined. For the most common conditions, the vertical convergence of the mesoscale heat flux acts to reduce the nocturnal cooling rate. Important uncertainties are discussed as well as the need for improved observations.     

From Near-Neutral to Strongly Stratified: Adequately Modelling the Clear-Sky Nocturnal Boundary Layer at Cabauw

The performance of an atmospheric single-column model (SCM) is studied systematically for stably-stratified conditions. To this end, 11 years (2005–2015) of daily SCM simulations were compared to observations from the Cabauw observatory, The Netherlands. Each individual clear-sky night was classified in terms of the ambient geostrophic wind speed with a \(1\hbox { m} \hbox { s}^{-1}\) bin-width. Nights with overcast conditions were filtered out by selecting only those nights with an average net radiation of less than \(-\,30\hbox { W }\hbox {m}^{-2}\). A similar procedure was applied to the observational dataset. A comparison of observed and modelled ensemble-averaged profiles of wind speed and potential temperature and time series of turbulent fluxes showed that the model represents the dynamics of the nocturnal boundary layer (NBL) at Cabauw very well for a broad range of mechanical forcing conditions. No obvious difference in model performance was found between near-neutral and strongly-stratified conditions. Furthermore, observed NBL regime transitions are represented in a natural way. The reference model version performs much better than a model version that applies excessive vertical mixing as is done in several (global) operational models. Model sensitivity runs showed that for weak-wind conditions the inversion strength depends much more on details of the land-atmosphere coupling than on the turbulent mixing. The presented results indicate that in principle the physical parametrizations of large-scale atmospheric models are sufficiently equipped for modelling stably-stratified conditions for a wide range of forcing conditions.     

Spectral Structure of Small-Scale Turbulent and Mesoscale Fluxes in the Atmospheric Boundary Layer over a Thermally Inhomogeneous Land Surface

Spectral analysis was performed on aircraft observations of a convective boundary layer (CBL) that developed over a thermally inhomogeneous, well-marked mesoscale land surface. The observations, part of the GAME-Siberia experiment, were recorded between April and June 2000 over the Lena River near Yakutsk City. A special integral parameter termed the ‘reduced depth of the CBL’ was used to scale the height of the mixed layer with variable depth. Analysis of wavelet cospectra and spectra facilitated the separation of fluxes and other variables into small-scale turbulent fluctuations (with scales less than the reduced depth of the CBL, approximately 2 km) and mesoscale fluctuations (up to 20 km). This separation approach allows for independent exploration of the scales. Analyses showed that vertical distributions obeyed different laws for small-scale fluxes and mesoscale fluxes (of sensible heat, water vapour, momentum and carbon dioxide) and for other variables (wind speed and air temperature fluctuations, coherence and degree of anisotropy). Vertical profiles of small-scale turbulent fluxes showed a strong decay that differed from generally accepted similarity models for the CBL. Vertical profiles of mesoscale fluxes and other variables clearly showed sharp inflections at the same relative (with respect to the reduced depth of the CBL) height of approximately 0.55 in the CBL. Conventional similarity models for sensible heat fluxes describe both small-scale turbulent and mesoscale flows. The present results suggest that mesoscale motions that reach up to the relative level of 0.55 could be initiated by thermal surface heterogeneity. Entrainment between the upper part of the CBL and the free atmosphere may cause mesoscale motions in that region of the CBL.     

An Aerosol Model of the Marine and Coastal Atmospheric Surface Layer

We present a microphysical model for the surface layer marine and coastal atmospheric aerosols that is based on long-term observations of size distributions for 0.01–100 µm particles. The fundamental feature of the model is a parameterization of amplitudes and widths for aerosol modes of the aerosol size distribution function (ASDF) as functions of fetch and wind speed. The shape of ASDF and its dependence on meteorological parameters, height above sea level (H), fetch (X), wind speed (U) and relative humidity (RH), are investigated. At present, the model covers the ranges H = 0–25 m, U?=?3–18 km s^?1, X?≤?120 km and RH?=?40–98%.

The latest version of the Marine Aerosol Extinction Profiles model (MaexPro) is described and applied to the computation and analysis of the spectral profiles of aerosol extinction coefficients α(λ) in the wavelength band λ?=?0.2–12 µm. MaexPro is based on the aforementioned aerosol model assuming spherically shaped aerosol particles and the well-known Mie theory.

The spectral profiles of α(λ) calculated by MaexPro are in good agreement with observational data and the numerical results obtained from the Navy Aerosol Model (NAM) and the Advanced Navy Aerosol Model (ANAM). Moreover, MaexPro was found to be an accurate and reliable tool for investigating the optical properties of atmospheric aerosols.     

Horizontal and Vertical Turbulent Fluxes Forced by a Gravity Wave Event in the Nocturnal Atmospheric Surface Layer Over the Amazon Forest

A nocturnal gravity wave was detected over a south-western Amazon forest during the Large-Scale Biosphere–Atmosphere experiment in Amazonia (LBA) in the course of the dry-to-wet season campaign on October 2002. The atmospheric surface layer was stably stratified and had low turbulence activity, based on friction velocity values. However, the passage of the wave, an event with a period of about 180–300 s, caused negative turbulent fluxes of carbon dioxide (CO₂) and positive sensible heat fluxes, as measured by the eddy-covariance system at 60 m (≈30 m above the tree tops). The evolution of vertical profiles of air temperature, specific humidity and wind speed during the wave movement revealed that cold and drier air occupied the sub-canopy space while high wind speeds were measured above the vegetation. The analysis of wind speed and scalars high frequency data was performed using the wavelet technique, which enables the decomposition of signals in several frequencies allowed by the data sampling conditions. The results showed that the time series of vertical velocity and air temperature were −90° out of phase during the passage of the wave, implying no direct vertical transport of heat. Similarly, the time series of vertical velocity and CO₂ concentration were 90° out of phase. The wave was not directly associated with vertical fluxes of this variable but the mixing induced by its passage resulted in significant exchanges in smaller scales as measured by the eddy-covariance system. The phase differences between horizontal velocity and both air temperature and CO₂ concentration were, respectively, zero and 180°, implying phase and anti-phase relationships. As a result, the wave contributed to positive horizontal fluxes of heat and negative horizontal fluxes of carbon dioxide. Such results have to be considered in nocturnal boundary-layer surface-atmosphere exchange schemes for modelling purposes.     

夜间近地面稳定边界层湍流间歇与增温

在夜间晴空条件下，近地面大气湍流表现出很强的间歇性，这种间歇现象导致夜间气温短时的急剧下降，随后大幅度增温。近地面大幅度增温表明此时存在着很大的湍流热通量散度，常通量层的概念这时不存在。从各高度层温度和风速变化的曲线上分析，我们发现湍流大多在距离地面较高一点的高度上发生发展，然后向下层传递，尽管上层的湍流可能是由于下层的某一触发机制向上传递而引起的。湍流偶尔也出现自下向上传递的过程，但这一过程较少发生。湍流的这种上下传递说明稳定边界层大气经常处于非平衡状态，在运用相似理论研究稳定边界层大气结构时要特别注意这一情形。     

Intermittency and the Exchange of Scalars in the Nocturnal Surface Layer

The determination of nocturnal surface fluxes in low wind conditions is a major problem for micrometeorological studies. The eddy correlation technique, extensively used in field measurements, becomes inappropriate if not enough turbulent activity exists. At the same time, the phenomenon of turbulence intermittency is responsible for the existence of localized events of short duration within which a large fraction of the total nighttime scalar exchange occurs. The scalar flux within a certain intermittent event varies considerably depending on the window used for the flux calculation. In many cases, events with very different time durations occur in the same night, and therefore, the proper determination of the surface flux would require averaging within data windows of different sizes for each event. In this work, the surface exchanges of temperature, moisture and carbon dioxide are analysed at a micrometeorological tower at southern Brazil. Intermittent turbulence is a common occurrence at the location. The analysis shows that the fluxes vary with turbulence intensity and the estimation technique. A variable-window size method for flux estimation is suggested and shown to cause an increase in the magnitude of the nocturnal surface fluxes     

A Statistical Model for the Prediction of Wind-Speed Probabilities in the Atmospheric Surface Layer

Wind fields in the atmospheric surface layer (ASL) are highly three-dimensional and characterized by strong spatial and temporal variability. For various applications such as wind-comfort assessments and structural design, an understanding of potentially hazardous wind extremes is important. Statistical models are designed to facilitate conclusions about the occurrence probability of wind speeds based on the knowledge of low-order flow statistics. Being particularly interested in the upper tail regions we show that the statistical behaviour of near-surface wind speeds is adequately represented by the Beta distribution. By using the properties of the Beta probability density function in combination with a model for estimating extreme values based on readily available turbulence statistics, it is demonstrated that this novel modelling approach reliably predicts the upper margins of encountered wind speeds. The model’s basic parameter is derived from three substantially different calibrating datasets of flow in the ASL originating from boundary-layer wind-tunnel measurements and direct numerical simulation. Evaluating the model based on independent field observations of near-surface wind speeds shows a high level of agreement between the statistically modelled horizontal wind speeds and measurements. The results show that, based on knowledge of only a few simple flow statistics (mean wind speed, wind-speed fluctuations and integral time scales), the occurrence probability of velocity magnitudes at arbitrary flow locations in the ASL can be estimated with a high degree of confidence.     

Time Scales in the Unstable Atmospheric Surface Layer

Calculation of eddy covariances in the atmospheric surface layer (ASL) requires separating the instantaneous signal into mean and fluctuating components. Since the ASL is not statistically stationary, an inherent ambiguity exists in defining the mean quantities. The present study compares four methods of calculating physically relevant time scales in the unstable ASL that may be used to remove the unsteady mean components of instantaneous time signals, in order to yield local turbulent fluxes that appear to be statistically stationary. The four mean-removal time scales are: (t _c ) based on the location of the maximum in the ogive of the heat flux cospectra, () the location of the zero crossing in the multiresolution decomposition of the heat flux, (t ^*) the ratio of the mixed-layer depth over the convective velocity, and () the convergence time of the vertical velocity and temperature variances. The four time scales are evaluated using high quality, three-dimensional sonic anemometry data acquired at the Surface Layer Turbulence and Environmental Science Test (SLTEST) facility located on the salt flats of Utah’s western desert. Results indicate that and , with t _c achieving values about 2–3 times greater than t ^*. The sensitivity of the eddy covariances to the mean-removal time scale (given a fixed 4-h averaging period during midday) is also demonstrated.     

Large-Scale Motions in the Marine Atmospheric Surface Layer

Multi-level turbulent wind data from the Risø Air-Sea Experiments (RASEX) were used to examine the structure of large-scale motions in the marine atmospheric surface layer. The quadrant technique was used to identify flux events (ejections/sweeps). Ejections, which appear to occur in groups, are seen to occur first at the upper level, moving successively to lower levels with small time delays. A strong correlation between events at different heights suggests that they may all be part of a single large structure. Cross-correlation between velocity signals was used to estimate orientation of the structure using Taylor's hypothesis. The inclination of this structure is shallow ( 15°) near the surface and increases with height. Spatial representations of the fluctuating wind vectors show a structure that is strikingly similar to conceptual models of transverse vortices and shear layers seen in laboratory flows and direct numerical simulation (DNS) of low Reynolds number flows. Spatial visualization of velocity fluctuations during other time periods and conditions clearly shows the existence of shear layers, transverse vortices, plumes, and downdrafts of various sizes and strengths. A quantitative analysis shows an increase in the frequency of shear related events with increasing wind speed.     

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.     

Stability Dependence and Diurnal Change of Large-Scale Turbulence Structures in the Near-Neutral Atmospheric Boundary Layer Observed from a Meteorological Tower

Based on the analysis of observations from a 213-m tall meteorological tower at Tsukuba, Japan, we have investigated the favourable conditions for the predominant existence of large-scale turbulence structures in the near-neutral atmospheric boundary layer (ABL). From the wavelet variance spectrum for the streamwise velocity component ( $u$ ) measured by a sonic anemometer-thermometer at the highest level (200 m), large-scale structures (time-scale range of 100–300 s) predominantly exist under slightly unstable and close to neutral conditions. The emergence of large-scale structures also can be related to the diurnal evolution of the ABL. The large-scale structures play an important role in the overall flow structure of the lower boundary layer. For example, $u$ velocity components at the 200-m and 50-m levels show relatively high correlation with the existence of large-scale structures. Under slightly unstable (near-neutral) conditions, a low-speed region in advance of the high-speed structure shows a positive deviation of temperature and appears as the plume structure that is forced by buoyancy in the heated lower layer. In spite of the difference in buoyancy effects between the near-neutral and unstable cases, large-scale structures are frequently observed in both cases and the same vertical correlation of $u$ components is indicated. However, the vertical wind shear is smaller in the unstable cases. On the other hand, in near-neutral cases, the transport efficiency of momentum at the higher level and the flux contribution of sweep motions are larger than those in the unstable cases.     

A Monte Carlo Model Of The Nocturnal Surface Temperatures In Urban Canyons

A model for the urban canyon is formulated for meteorologicalconditions of weak winds at night time. Thermal radiation, conductivity and convection are simulated by means of the Monte Carlo method. These are the main physical processesof energy transfer that give rise to the characteristic temperaturedistribution in these systems. The model has been satisfactory tested under ideal conditions for which analytical solutions exist.The predictions of the model under morerealistic conditions accurately reproduce the observationalresults. A strong temperature gradient across streets, with the canyon corners up to 4 °C warmer than the canyon centre, is found for the deepest canyons. This theoretical predictionhas been successfully verified with measurementstaken in a number of streets of the city of Granada in Spain.     

Impact of Land Surface Heterogeneity on Mesoscale Atmospheric Dispersion

Prior numerical modelling studies show that atmospheric dispersion is sensitive to surface heterogeneities, but past studies do not consider the impact of a realistic distribution of surface heterogeneities on mesoscale atmospheric dispersion. While these focussed on dispersion in the convective boundary layer, the present work also considers dispersion in the nocturnal boundary layer and above. Using a Lagrangian particle dispersion model (LPDM) coupled to the Eulerian Regional Atmospheric Modeling System (RAMS), the impact of topographic, vegetation, and soil moisture heterogeneities on daytime and nighttime atmospheric dispersion is examined. In addition, the sensitivity to the use of Moderate Resolution Imaging Spectroradiometer (MODIS)-derived spatial distributions of vegetation characteristics on atmospheric dispersion is also studied. The impact of vegetation and terrain heterogeneities on atmospheric dispersion is strongly modulated by soil moisture, with the nature of dispersion switching from non-Gaussian to near-Gaussian behaviour for wetter soils (fraction of saturation soil moisture content exceeding 40%). For drier soil moisture conditions, vegetation heterogeneity produces differential heating and the formation of mesoscale circulation patterns that are primarily responsible for non-Gaussian dispersion patterns. Nighttime dispersion is very sensitive to topographic, vegetation, soil moisture, and soil type heterogeneity and is distinctly non-Gaussian for heterogeneous land-surface conditions. Sensitivity studies show that soil type and vegetation heterogeneities have the most dramatic impact on atmospheric dispersion. To provide more skilful dispersion calculations, we recommend the utilisation of satellite-derived vegetation characteristics coupled with data assimilation techniques that constrain soil-vegetation-atmosphere transfer (SVAT) models to generate realistic spatial distributions of surface energy fluxes.