Boundary-Layer Meteorology最新文献

Taylor’s Frozen Turbulence Hypothesis (TH) is a critical assumption in turbulent theory and practice which allows time series of point measurements of turbulent variables to be translated to the spatial domain via the mean wind. Using a 3D array of fibre-optic distributed temperature sensing in the atmospheric surface layer over an idealized desert site we present a systematic investigation of the applicability of Taylor’s Hypothesis to atmospheric surface layer flows over a variety of conditions: unstable, near-neutral, and stable atmospheric stabilities; and multiple measurement heights between the surface and 3 m above ground level. Both spatially integrated and spatially scale-dependent eddy velocities are investigated by means of time-lagged streamwise two-point correlations and compared to the mean Eulerian wind. We find that eddies travel slower than predicted by TH at small spatial separations, as predicted by TH at separations typically between 5 and 16 m, and faster than predicted by TH at larger spatial separations. In unstable atmospheric conditions the spatial separation at which eddy velocity is larger than Eulerian velocity decreases with height.

Wildland fire–atmosphere interaction generates complex turbulence patterns, organized across multiple scales, which inform fire-spread behaviour, firebrand transport, and smoke dispersion. Here, we utilize wavelet-based techniques to explore the characteristic temporal scales associated with coherent patterns in the measured temperature and the turbulent fluxes during a prescribed wind-driven (heading) surface fire beneath a forest canopy. We use temperature and velocity measurements from tower-mounted sonic anemometers at multiple heights. Patterns in the wavelet-based energy density of the measured temperature plotted on a time–frequency plane indicate the presence of fire-modulated ramp–cliff structures in the low-to-mid-frequency band (0.01–0.33 Hz), with mean ramp durations approximately 20% shorter and ramp slopes that are an order of magnitude higher compared to no-fire conditions. We then investigate heat- and momentum-flux events near the canopy top through a cross-wavelet coherence analysis. Briefly before the fire-front arrives at the tower base, momentum-flux events are relatively suppressed and turbulent fluxes are chiefly thermally-driven near the canopy top, owing to the tilting of the flame in the direction of the wind. Fire-induced heat-flux events comprising warm updrafts and cool downdrafts are coherent down to periods of a second, whereas ambient heat-flux events operate mainly at higher periods (above 17 s). Later, when the strongest temperature fluctuations are recorded near the surface, fire-induced heat-flux events occur intermittently at shorter scales and cool sweeps start being seen for periods ranging from 8 to 35 s near the canopy top, suggesting a diminishing influence of the flame and increasing background atmospheric variability thereat. The improved understanding of the characteristic time scales associated with fire-induced turbulence features, as the fire-front evolves, will help develop more reliable fire behaviour and scalar transport models.

This study investigates flow variability at different scales and its effects on the dispersion of a passive scalar in a regular street network by means of direct numerical simulations (DNS), and compared to wind tunnel (WT) measurements. Specific scientific questions addressed include: (i) sources of variability in the flow at street-network scale, (ii) the effects of such variability on both puff and continuous localised releases, (iii) additional sources of uncertainty related to experimental setups and their consequences. The street network modelled here consists of an array of rectangular buildings arranged uniformly and with periodic horizontal boundary conditions. The flow is driven by a body force at an angle of 45 degrees relative to the streets in the network. Sources of passive scalars were located near ground level at three different types of locations: a short street, an intersection between streets and a long street. Flow variability is documented at different scales: small-scale intra-street variations linked with local flow topology; inter-street flow structure differences; street-network scale variability; and larger-scale spatial variations associated with above-canopy structures. Flow statistics and the dispersion behaviour of both continuous and short-duration (puff) releases of a passive scalar in the street network are analysed and compared with the results of wind-tunnel measurements. Results agree well with the experimental data for a source location in an intersection, especially for flow statistics and mean concentration profiles for continuous releases. Larger differences arise in the comparisons of puff releases. These differences are quantified by computing several puff parameters including time of arrival, travel time, rise and decay times. Reasons for the differences are discussed in relation to the underlying flow variability identified, differences between the DNS and WT setup and uncertainties in the experimental setup. Implications for the propagation of short-duration releases in real urban areas are discussed in the light of our findings. In particular, it is highlighted that in modelling singular events such as accidental releases, characterising uncertainties is more meaningful and useful than computing ensemble averages.

A single-column turbulence model for stratified atmospheric boundary layer (ABL), which solves the transport equations of turbulence probability density function (PDF) using a Lagrangian stochastic modeling (LSM) approach, is proposed in this study. This study adopts previously developed stochastic differential equations (SDEs) for particle velocity and temperature and extends the LSM to simulate inhomogeneous turbulence. The proposed LSM is tested for its ability to fully simulate statistics of inhomogeneous stratified turbulence. In the model, particles evolve by SDEs, and turbulence statistics are calculated by averaging the properties of particles. The model provides a full representation of turbulence PDF and simulates turbulent transport without any modeling assumption. The model performance is evaluated against large-eddy simulation (LES) results in the simulations of convective and stable ABL cases. For the convective ABL, LSM realistically simulates the entrainment process with the temperature and heat flux profiles that closely match with LES. The joint PDF simulated by LSM reproduces a curved and highly skewed shape, and some distinct features, like the asymmetric distribution of vertical velocity and the separation of the PDF in the entrainment zone, are simulated. LSM also reproduces the entrainment enhancement by wind shear in the simulation of sheared convective ABL. The LSM simulation of stable ABL predicts realistic turbulence intensity and mean field profiles, where Gaussian-like PDFs are simulated both in LSM and LES.

Abstract

Recent studies have highlighted the importance of accurate meteorological conditions for urban transport and dispersion calculations. In this work, we present a novel scheme to compute the meteorological input in the Quick Urban & Industrial Complex (QUIC) diagnostic urban wind solver to improve the characterization of upstream wind veer and shear in the Atmospheric Boundary Layer (ABL). The new formulation is based on a coupled set of Ordinary Differential Equations (ODEs) derived from the Reynolds Averaged Navier–Stokes (RANS) equations, and is fast to compute. Building upon recent progress in modeling the idealized ABL, we include effects from surface roughness, turbulent stress, Coriolis force, buoyancy and baroclinicity. We verify the performance of the new scheme with canonical Large Eddy Simulation (LES) tests with the GPU-accelerated FastEddy solver in neutral, stable, unstable and baroclinic conditions with different surface roughness. Furthermore, we evaluate QUIC calculations with and without the new inflow scheme with real data from the Urban Threat Dispersion (UTD) field experiment, which includes Lidar-based wind measurements as well as concentration observations from multiple outdoor releases of a non-reactive tracer in downtown New York City. Compared to previous inflow capabilities that were limited to a constant wind direction with height, we show that the new scheme can model wind veer in the ABL and enhance the prediction of the surface cross-isobaric angle, improving evaluation statistics of simulated concentrations paired in time and space with UTD measurements.

It has long been known that under the right circumstances, inertial particles (such as sand, dust, pollen, or water droplets) settling through the atmospheric boundary layer can experience a net enhancement in their average settling velocity due to their inertia. Since this enhancement arises due to their interactions with the surrounding turbulence it must be modelled at coarse scales. Models for the enhanced settling velocity (or deposition) of the dispersed phase that find practical use in mesoscale weather models are often ad hoc or are built on phenomenological closure assumptions, meaning that the general deposition rate of particles is a key uncertainty in these models. Instead of taking a phenomenological approach, exact phase-space methods can be used to model the physical mechanisms responsible for the enhanced settling, and these individual mechanisms can be estimated or modelled to build a more general parameterization of the enhanced settling of inertial particles. In this work, we use direct numerical simulations (DNS) and phase-space methods as tools to evaluate the efficacy of phenomenological modeling approaches for the enhanced settling velocity of inertial particles for particles with varying friction Stokes numbers and settling velocity parameters. We use the DNS data to estimate profiles of a drift–diffusion based parameterization of the fluid velocity sampled by the particles, which is key for determining the settling velocity behaviour of particles with low to moderate Stokes number. We find that by increasing the settling velocity parameter at moderate friction Stokes number, the magnitude of preferential sweeping is modified, and this behaviour is explained by the drift component of the aforementioned parameterization. These profiles indicate that that when eddy-diffusivity-like closures are used to represent turbulent transport, empirical corrections used in phenomenological models may be potentially compensating for their incompleteness. Finally, we discuss opportunities for reinterpreting phenomenological approaches for use in coarse-scale weather models in terms of the exact phase-space approach.

We introduce an analytical model that describes the vertical structure of Ekman boundary layer flows coupled to the Monin-Obukhov Similarity Theory (MOST) surface layer representation, which is valid for conventionally neutral (CNBL) and stable (SBL) atmospheric conditions. The model is based on a self-similar profile of horizontal stress for both CNBL and SBL flows that merges the classic 3/2 power law profile with a MOST-consistent stress profile in the surface layer. The velocity profiles are then obtained from the Ekman momentum balance equation. The same stress model is used to derive a new self-consistent Geostrophic Drag Law (GDL). We determine the ABL height (h) using an equilibrium boundary layer height model and parameterize the surface heat flux for quasi-steady SBL flows as a function of a prescribed surface temperature cooling rate. The ABL height and GDL equations can then be solved together to obtain the friction velocity ((u_*)) and the cross-isobaric angle ((alpha _0)) as a function of known input parameters such as the Geostrophic wind speed and surface roughness ((z_0)). We show that the model predictions agree well with simulation data from the literature and newly generated Large Eddy Simulations (LES). These results indicate that the proposed model provides an efficient and relatively accurate self-consistent approach for predicting the mean wind velocity distribution in CNBL and SBL flows.

We examine measurements in the very stable boundary layer using tower data and a network of flux stations in the Shallow Cold Pool experiment. Submeso motions in the very stable boundary layer significantly modulate the turbulent heat fluxes and also directly contribute to the submeso vertical heat flux. Time series include well-defined submeso structures such as microfronts, wave-like motions, and meandering but also include complex structures that are difficult to isolate. These structures significantly influence the time and height variation of the turbulent heat flux. From the 19 flux stations distributed across the shallow valley, we find that the surface heat flux with low wind speeds varies significantly on a horizontal scale of 100 m, or less, related partly to the modest topography. For this dataset, the turbulent surface heat fluxes for low wind speeds are closely related to submeso variations of the wind speed but not significantly related to variations of the stratification.

Abstract

Turbulent exchanges between sea ice and the atmosphere are known to influence the melting rate of sea ice, the development of atmospheric circulation anomalies and, potentially, teleconnections between polar and non-polar regions. Large model errors remain in the parametrization of turbulent heat fluxes over sea ice in climate models, resulting in significant uncertainties in projections of future climate. Fluxes are typically calculated using bulk formulae, based on Monin-Obukhov similarity theory, which have shown particular limitations in polar regions. Parametrizations developed specifically for polar conditions (e.g. representing form drag from ridges or melt ponds on sea ice) rely on sparse observations and thus may not be universally applicable. In this study, new data-driven parametrizations have been developed for surface turbulent fluxes of momentum, sensible heat and latent heat in the Arctic. Machine learning has already been used outside the polar regions to provide accurate and computationally inexpensive estimates of surface turbulent fluxes. To investigate the feasibility of this approach in the Arctic, we have fitted neural-network models to a reference dataset (SHEBA). Predictive performance has been tested using data from other observational campaigns. For momentum and sensible heat, performance of the neural networks is found to be comparable to, and in some cases substantially better than, that of a state-of-the-art bulk formulation. These results offer an efficient alternative to the traditional bulk approach in cases where the latter fails, and can serve to inform further physically based developments.

Abstract

Large-eddy simulations (LES) above forests and cities typically constrain the simulation domain to the first 10–20% of the Atmospheric Boundary Layer (ABL), aiming to represent the finer details of the roughness elements and sublayer. These simulations are also commonly driven by a constant pressure gradient term in the streamwise direction and zero stress at the top, resulting in an unrealistic fast decay of the total stress profile. In this study, we investigate five LES setups, including pressure and/or top-shear driven flows with and without the Coriolis force, with the aim of identifying which option best represents turbulence profiles in the atmospheric surface layer (ASL). We show that flows driven solely by pressure not only result in a fast-decaying stress profile, but also in lower velocity variances and higher velocity skewnesses. Top-shear driven flows, on the other hand, better replicate ASL statistics. Overall, we recommend, and provide setup guidance for, simulation designs that include both a large scale pressure forcing and a non-zero stress and scalar flux at the top of the domain, and that also represent the Coriolis force. Such setups retain all the forces used in typical full ABL cases and result in the best match of the profiles of various statistical moments.