Search

This study investigates the parameterization of the geostrophic drag law (GDL) for conventionally neutral atmospheric boundary layers (CNBLs). Utilizing large eddy simulations, we confirm that in CNBLs capped by a potential temperature inversion, the boundary-layer height scales as u ∗ / Nf , where u ∗ represents the friction velocity, N the free-atmosphere Brunt–Väisälä frequency, and f the Coriolis parameter. Additionally, we confirm that the wind gradients normalized by the Brunt–Väisälä frequency have universal profiles above the surface layer. Leveraging these physical insights, we derived analytical expressions for the GDL coefficients A and B, correcting the earlier form of Zilitinkevich and Esau (Q J R Meteorol Soc 131:1863–1892, 2005). These expressions for A and B have been validated numerically, ensuring their accuracy in representing the geostrophic drag coefficient u ∗ / G (G is the geostrophic wind speed) and the cross-isobaric angle. This work extends the range for which the GDL has been validated up to u ∗ / G = [ 0.019 , 0.047 ] . This further supports the application of GDL to CNBLs over a broader range of u ∗ / G , which is useful for meteorological applications such as wind energy. [ABSTRACT FROM AUTHOR]

The impact of clouds on the mixed layer (ML) is critical for understanding the evolution of boundary layer humidity and temperature over the course of a day. We found that accounting for moistening of the cloud layer (CL) by humidity originating in the ML dramatically alters the interaction between the ML and the CL in a one‐dimensional cloud‐topped boundary layer model: Chemistry Land‐surface Atmosphere Soil Slab (CLASS) (Vilà‐Guerau de Arellano et al., 2015, https://doi.org/10.1017/CBO9781316117422). We demonstrate that enabling CLASS to moisten the lower CL improves the prediction of humidity (and the flux of humidity) both above and below the ML top (h). To account for this moistening, we propose a length scale, L, above h, over which mixing of mass fluxes into the environment occurs. The mass fluxes are assumed to decrease linearly from h to a height L meters above h, analogous to a convective plume detraining into the environment at a height‐independent rate. Accounting for this process is accomplished by modifying the differential equations representing the growth of the jumps (sharp changes in humidity and temperature) at h. From analysis of a large number of diurnal Large Eddy simulations (covering approximately 11,000 different early morning initial conditions), we provide a regression model for parameterizing L from early morning weather variables. With the regression‐based estimate of L, the modified model (CLASS‐L) accounts for moistening the lower CL, and as a result, yields improved humidity dynamics, humidity flux profiles, and cloud growth under a broad range of conditions. Plain Language Summary: This study improved the representation of clouds above the mixed layer (ML) top in a simplified cloud‐topped atmospheric boundary layer model, Chemistry Land‐surface Atmosphere Soil Slab (CLASS) (Vilà‐Guerau de Arellano et al., 2015, https://doi.org/10.1017/CBO9781316117422). The CLASS model does not moisten the lower cloud layer (CL) when moisture originating from the ML reaches the CL. To address this, we propose a mixing length scale, L, above the ML top over which moisture originating from the ML is mixed into the environment, and name this model CLASS‐L. A power‐law regression is used to estimate L with regression coefficients based on a large set of daytime diurnal simulations (covering approximately 11,000 different initial conditions) using CLASS‐L and state‐of‐the‐art, three‐dimensional Large Eddy Simulation. With this simple enhancement to CLASS, CLASS‐L dramatically improved humidity dynamics, cloud growth, and humidity flux profiles both above and within the ML under a very broad range of initial conditions. Key Points: We enhance land‐atmosphere‐cloud interactions in a simplified model (Vilà‐Guerau de Arellano et al., 2015, https://doi.org/10.1017/CBO9781316117422)We propose an empirically determined length scale above the mixed layer over which mixing of humidity mass flux with the environment occursThe modifications (Chemistry Land‐surface Atmosphere Soil Slab‐L) dramatically alter the state and flux profiles in cloudy conditions, in agreement with Large Eddy Simulations [ABSTRACT FROM AUTHOR]

An idealized large-eddy simulation of a tropical marine cloud population was performed. At any time, it contained hundreds of clouds, and updraft width in shallow convection emerging from a subcloud layer appeared to be an important indicator of whether specific convective elements deepened. In an environment with 80%–90% relative humidity below the 0°C level, updrafts that penetrated the 0°C level were larger at and above cloud base, which occurred at the lifting condensation level near 600 m. Parcels rising in these updrafts appeared to emerge from boundary layer eddies that averaged ∼200 m wider than those in clouds that only reached 1.5–3 km height. The deeply ascending parcels (growers) possessed statistically similar values of effective buoyancy below the level of free convection (LFC) as parcels that began to ascend in a cloud but stopped before reaching 3000 m (nongrowers). The growers also experienced less dilution above the LFC. Nongrowers were characterized by negative effective buoyancy and rapid deceleration above the LFC, while growers continued to accelerate well above the LFC. Growers occurred in areas with a greater magnitude of background convergence (or weaker divergence) in the subcloud layer, especially between 300 m and cloud base, but whether the convergence actually led to eddy widening is unclear. Significance Statement: Cumulonimbus clouds are responsible for many extreme weather phenomena and are important contributors to Earth's energy balance. However, the processes leading to the growth of individual clouds are not completely understood nor well-represented in weather prediction models. We find that the clouds containing updrafts that start out wider at early stages of their life cycles grow taller, possibly because they are protected more from drier air outside the cloud than narrow clouds. In addition, this work shows how the initial width of clouds might be related to convergence in the lowest part of the atmosphere, at heights where clouds initially develop. However, meteorologists must be careful not to overinterpret these results because numerical simulations inherently include assumptions that may not reflect reality. This reinforces the need to also observe processes occurring at the scales of individual clouds. [ABSTRACT FROM AUTHOR]

To address critical gaps identified by the National Academies of Sciences, Engineering, and Medicine in the current Earth system observation strategy, the 2017–27 Decadal Survey for Earth Science and Applications from Space recommended incubating concepts for future targeted observables including the atmospheric planetary boundary layer (PBL). A subsequent NASA PBL Incubation Study Team Report identified measurement requirements and activities for advancing the maturity of the technologies applicable to the PBL targeted observables and their associated science and applications priorities. While the PBL is the critical layer where humans live and surface energy, moisture, and mass exchanges drive the Earth system, it is also the farthest and most inaccessible layer for spaceborne instruments. Here we document a PBL retrieval observing system simulation experiment (OSSE) framework suitable for assessing existing and new measurement techniques and determining their accuracy and improvements needed for addressing the elevated Decadal Survey requirements. In particular, the benefits of large-eddy simulation (LES) are emphasized as a key source of high-resolution synthetic observations for key PBL regimes: from the tropics, through subtropics and midlatitudes, to subpolar and polar regions. The potential of LES-based PBL retrieval OSSEs is explored using six instrument simulators: Global Navigation Satellite System–Radio Occultation, differential absorption radar, visible to shortwave infrared spectrometer, infrared sounder, Multi-angle Imaging SpectroRadiometer, and microwave sounder. The crucial role of LES in PBL retrieval OSSEs and some perspectives for instrument developments are discussed. [ABSTRACT FROM AUTHOR]

Land–atmosphere interactions play a critical role in both the atmospheric water and energy cycles. Changes in soil moisture and vegetation alter the partitioning of surface water and energy fluxes, influencing diurnal evolution of the planetary boundary layer (PBL). The mixing-diagram framework has proven useful in understanding the evolution of the heat and moisture budget within the convective boundary layer (CBL). We demonstrate that observations from the Department of Energy Atmospheric Radiation Measurement (ARM) Southern Great Plains site provide all of the needed inputs needed for the mixing-diagram framework, allowing us to quantify the impact from the surface fluxes, advection, radiative heating, encroachment, and entrainment on the evolution of the CBL. Profiles of temperature and humidity retrieved from the ground-based infrared spectrometer [Atmospheric Emitted Radiance Interferometer (AERI)] are a critical component in this analysis. Large-eddy simulation results demonstrate that mean mixed-layer values derived are shown to be critical to close the energy and moisture budgets. A novel approach demonstrated here is the use of network of AERIs and Doppler lidars to quantify the advective fluxes of heat and moisture. The framework enables the estimation of the entrainment fluxes as a residual, providing a way to observe the entrainment fluxes without using multiple lidar systems. The high temporal resolution of the AERI observations enables the morning, midday, and afternoon evolution of the CBL to be quantified. This work provides a new way to use observations in this framework to evaluate weather and climate models. Significance Statement: The energy and moisture budget of the planetary boundary layer (PBL) is influenced by multiple sources, and accurately representing this evolution in numerical models is critical for weather forecasts and climate predictions. The mixing-diagram approach, driven by profiling observations as illustrated here, provides a powerful way to quantify the contributions from each of these sources. In particular, the energy and moisture mixed into the PBL from above the PBL can be determined accurately from ground-based remote sensors using this approach. [ABSTRACT FROM AUTHOR]

In large-eddy simulations (LES), it is crucial to ensure that discretization errors do not contaminate the subgrid effect of the turbulence model in a wavelength range larger than the effective resolution. Recently, we showed that a seventh- or eighth-order accuracy is required for advection terms in planetary boundary layer simulations when using conventional gridpoint methods. However, a significant amount of communication between parallel computers is necessary to achieve high-order accuracy in gridpoint methods, and this can degrade computational efficiency. The discontinuous Galerkin method (DGM) is a promising approach for overcoming these limitations. Therefore, this study focuses on the numerical criteria of the DGM at LES from the viewpoint of numerical diffusion and dispersion. We extend our earlier study to the DGM framework and clarify the necessary order of the polynomial (p). We find that p = 4 is required based on the numerical criteria at the grid spacing of O(10) m with sufficiently scale-selective modal filters. The examination of temporal accuracy suggests that the fourth-order is sufficient when a fully explicit temporal scheme is used. In addition, we investigate the effect of hyperupwinding that is usually met when the Rusanov flux is employed in the low Mach number flows. It suggests that the choice of numerical flux has little effect on simulation results when the high-order DGM is used. Furthermore, we perform a series of LES in the planetary boundary layer and confirm that the indication obtained from the criteria holds for an actual LES. [ABSTRACT FROM AUTHOR]

Unsteadiness and horizontal heterogeneities frequently characterize atmospheric motions, especially within convective storms, which are frequently studied using large-eddy simulations (LES). The models of near-surface turbulence employed by atmospheric LES, however, predominantly assume statistically steady and horizontally homogeneous conditions (known as the equilibrium approach). The primary objective of this work is to investigate the potential consequences of such unrealistic assumptions in simulations of tornadoes. Cloud Model 1 (CM1) LES runs are performed using three approaches to model near-surface turbulence: the "semi-slip" boundary condition (which is the most commonly used equilibrium approach), a recently proposed nonequilibrium approach that accounts for some of the effects of turbulence memory, and a nonequilibrium approach based on thin boundary layer equations (TBLE) originally proposed by the engineering community for smooth-wall boundary layer applications. To be adopted for atmospheric applications, the TBLE approach is modified to account for the surface roughness. The implementation of TBLE into CM1 is evaluated using LES results of an idealized, neutral atmospheric boundary layer. LES runs are then performed for an idealized tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. The semi-slip boundary condition, by design, always yields a surface shear stress opposite the horizontal wind at the lowest LES grid level. The nonequilibrium approaches of modeling near-surface turbulence allow for a range of surface-shear-stress directions and enhance the resolved turbulence and wind gusts. The TBLE approach even occasionally permits kinetic energy backscatter from unresolved to resolved scales. Significance Statement: The traditional approach of modeling the near-surface turbulence is not suitable for a tornado characterized by rapid evolution, strongly curved air parcel trajectories, and substantial horizontal heterogeneities. To understand the influence of statistically unsteady and horizontally heterogeneous near-surface conditions on tornadoes, this work adopts a fairly sophisticated approach from the engineering community and implements it into a widely used atmospheric model with necessary modifications. Compared to the traditional approach, the newly implemented approach produces more turbulent near-surface winds, more flexible surface-drag directions, and stronger wind gusts. These findings suggest a simulated tornado is very sensitive to the modeling approach of near-surface turbulence. [ABSTRACT FROM AUTHOR]

For decades, stable boundary layer (SBL) turbulence has proven challenging to measure, parameterize, simulate, and interpret. Uncrewed aircraft systems (UAS) are becoming a reliable method to sample the atmospheric boundary layer, offering new perspectives for understanding the SBL. Moreover, continual computational advances have enabled the use of large-eddy simulations (LES) to simulate the atmosphere at ever-smaller scales. LES is therefore a powerful tool in establishing a baseline framework to understand the extent to which vertical profiles from UAS can represent larger-scale SBL flows. To quantify the representativeness of observations from UAS profiles and eddy-covariance observations within the SBL, we performed a random error analysis using a suite of six large-eddy simulations for a wide range of stabilities. We combine these random error estimates with emulated observations of a UAS and eddy-covariance systems to better inform future observational studies. For each experiment, we estimate relative random errors using the so-called relaxed filtering method for first- and second-order moments as functions of height and averaging time. We show that the random errors can be on the same order of magnitude as other instrument-based errors due to bias or dynamic response. Unlike instrument errors, however, random errors decrease with averaging time. For these reasons, we recommend coupling UAS observations with other ground-based instruments as well as dynamically adjusting the UAS vertical ascent rate to account for how errors change with height and stability. Significance Statement: Weather-sensing uncrewed aircraft systems are rapidly being realized as effective tools to collect valuable observations within the atmospheric boundary layer. To fully capitalize on this novel observational technique, it is necessary to develop an understanding of how well their observations can represent the surrounding atmosphere across various spatial and temporal scales. In this study we quantify the representativeness of atmospheric observations in the stable boundary layer by evaluating the random errors for parameters such as temperature, wind speed, and fluxes as estimated from a suite of large-eddy simulations. Our results can better inform future studies utilizing uncrewed aircraft systems by highlighting how random errors in their observations relate to vertical ascent rate, atmospheric stability, and measurement height. [ABSTRACT FROM AUTHOR]

Large-eddy simulations (LES) of neutrally and stably stratified turbulent flows over urban-type surfaces with relatively low plan area ratios are presented. Numerical experiments were performed for different shapes of streamlined objects and at different static stability. A new method for setting up a numerical experiment aimed at studying the heat and momentum transfer within the roughness layer and investigating the thermal and dynamic interaction between the turbulent flow and the surface as a whole has been proposed. This method enables us to obtain an equilibrium state for values of parameters determining the characteristics of the external turbulent flow chosen beforehand. A strong dependence of the thermal roughness length on stratification was found. We also discuss the physical mechanisms that lead to the maintenance of turbulence above the canopy when the ground surface is strongly cooled. Significance Statement: Using LES, we identify a mechanism that contributes to the maintenance of turbulence in the atmospheric boundary layer under the condition of strong surface cooling. Although these results are obtained for an urban canopy, we believe that the qualitative conclusions should be general for a wide type of surfaces with large-scale roughness elements. We hope that the new results will be useful for improving surface flux schemes in NWP and climate atmospheric models that suffer from attenuated mixing in a very stable boundary layer and the effect of "surface decoupling." The found effect gives a physically justified alternative way to parameterize the air–surface exchange under strong stability compared to the often ad hoc modification of the MOST universal functions. [ABSTRACT FROM AUTHOR]

In the numerical gray zone of the convective boundary layer (CBL), the horizontal resolution is comparable to the size of organized convective circulation. As turbulence becomes partially resolved, gridscale variations of the subgrid-scale (SGS) turbulent fluxes become significant compared to the mean. Previously, such variations have often been ignored in scale-adaptive planetary boundary layer (PBL) schemes developed for the gray zone. This study investigates these variations with respect to height and resolution based on large-eddy simulations. It is found that SGS fluxes exhibit maximum variability at the center of the gray zone, where the resolved and the SGS mean fluxes are approximately equal. A simple analytical model is used to associate such characteristic variations to the nonlinear interactions of the dominant energy-containing mode of CBL turbulence. Examination of the horizontal distribution of the SGS fluxes reveals their preferential location over the updraft edges surrounding the core. A priori analysis further suggests the ability of a scale-similarity closure to reproduce the unique spatial patterns of the SGS fluxes at gray zone resolutions. Four scale-adaptive PBL schemes are evaluated focusing on their representations of the modeled SGS flux variability. Their shared shortcomings as a result of their gradient diffusion–based formulation are exposed. This study suggests that a mixed model consisting of a scale-adaptive PBL scheme to represent the mean, and a scale-similarity component to account for gridscale variability to be advantageous for the gray zone. Significance Statement: As Gresho and Lee's famous quote on numerical schemes goes, "Don't suppress the wiggles. They're telling you something!" In the numerical gray zone, where the grid spacing is comparable to the characteristic length scale of the flow, turbulence is partially resolved. The "wiggles" (or gridscale variability) reflecting the resolved heterogeneity of turbulent fluxes is an outstanding feature of the gray zone. However, they are often overlooked as error bars to the mean. This study uncovers the significance of these error bars, and characterizes and explains their variations with height and grid spacing. A suitable model that captures such variations is investigated to help build better planetary boundary layer schemes. [ABSTRACT FROM AUTHOR]

Imagery and numerical modeling show an abundance of submesoscale oceanic eddies in the upper ocean. Large-eddy simulation (LES) is used to elucidate eddy impacts on the atmospheric boundary layer (ABL) forced by winds, convection, and an eddy with varying radius; the maximum azimuthal eddy speed is 1 m s−1. Simulations span the unstable regime −1/L = [0, ∞], where L is the Monin–Obukhov (M–O) stability parameter. A linearized Ekman model and the LES couple ABL winds to an eddy through rough-wall M–O boundary conditions. The eddy currents cause a surface stress anomaly that induces Ekman pumping in a dipole horizontal pattern. The dipole is understood as a consequence of surface winds aligned or opposing surface currents. In free convection a vigorous updraft is found above the eddy center and persists over the ABL depth. Heterogeneity in surface temperature flux is responsible for the full ABL impact. With winds and convection, current stress coupling generates a dipole in surface temperature flux even with constant sea surface temperature. Wind, pressure, and temperature anomalies are sensitive to an eddy under light winds. The eddy impact on ABL secondary circulations is on the order of the convective velocity scale w * but grows with increasing current speed, decreasing wind, or increasing convection. Flow past an isolated eddy develops a coherent ABL "wake" and secondary circulations for at least five eddy radii downwind. Kinetic energy exchanges by wind work indicate an eddy-killing effect on the oceanic eddy current, but only a spatial rearrangement of the atmospheric wind work. [ABSTRACT FROM AUTHOR]

This study investigates simulated fair-weather shallow cumulus-topped boundary layer (SCTBL) on kilometer- and subkilometer-scale horizontal resolutions, also known as the numerical gray zone of boundary layer turbulence. Based on a priori analysis of a simulated classic SCTBL with large-eddy simulation, its gray zone scale is determined. The dominant length scale of the cloud layer (CL) is found to be the effective cloud diameter, while that of the underlying mixed layer (ML) is the size of organized convection. The two scales are linked by a simple geometric argument based on vertically coherent updrafts, and are quantified through spectral analysis. Comparison to a simulated dry convective boundary layer (CBL) further reveals that the ML gray zone scale does not differentiate between clear and cloudy conditions with the same bulk stability. A posteriori simulations are then performed over a range of resolutions to evaluate the performance of a recently developed scale-adaptive planetary boundary layer (PBL) scheme. Simulation results suggest indifferences to the scale-adaptive capability. Detailed analyses of flux partition reveal that, in the absence of a shallow cumulus scheme, overly energetic resolved fluxes develop in the CL at gray zone and coarse resolutions, and are responsible for overpredicted resolved convection in the ML. These results suggest that modifications are needed for scale-adaptive PBL schemes under shallow cumulus-topped conditions. Significance Statement: Shallow cumulus (ShCu) clouds play an important role in the dynamical and radiative processes of the atmospheric boundary layer. As the grid resolution of modern numerical weather prediction models approach kilometer and subkilometer scales, also known as the gray zone, accurate modeling of ShCu clouds becomes challenging due to difficulties in their parameterization. This study identifies the spatial scale that sets the gray zone of ShCu clouds, providing the key to building better parameterizations. Performance of existing parameterizations developed for clear-sky conditions is evaluated for cloudy conditions, exposing deficiencies and motivating further development. [ABSTRACT FROM AUTHOR]

Accurately representing boundary layer turbulent processes in numerical models is critical to improve tropical cyclone forecasts. A new turbulence kinetic energy (TKE)-based moist eddy-diffusivity mass-flux (EDMF-TKE) planetary boundary layer scheme has been implemented in NOAA's Hurricane Analysis and Forecast System (HAFS). This study evaluates EDMF-TKE in hurricane conditions based on a recently developed framework using large-eddy simulation (LES). Single-column modeling tests indicate that EDMF-TKE produces much greater TKE values below 500-m height than LES benchmark runs in different high-wind conditions. To improve these results, two parameters in the TKE scheme were modified to ensure a match between the PBL and surface-layer parameterizations. Additional improvements were made by reducing the maximum allowable mixing length to 40 m based on LES and observations, by adopting a different definition of boundary layer height, and by reducing nonlocal mass fluxes in high-wind conditions. With these modifications, the profiles of TKE, eddy viscosity, and winds compare much better with LES results. Three-dimensional idealized simulations and an ensemble of HAFS forecasts of Hurricane Michael (2018) consistently show that the modified EDMF-TKE tends to produce a stronger vortex with a smaller radius of maximum wind than the original EDMF-TKE, while the radius of gale-force wind is unaffected. The modified EDMF-TKE code produces smaller eddy viscosity within the boundary layer compared to the original code, which contributes to stronger inflow, especially within the annulus of 1–3 times the radius of maximum wind. The modified EDMF-TKE shows promise to improve forecast skill of rapid intensification in sheared environments. [ABSTRACT FROM AUTHOR]

Boundary layer turbulence in coastal regions differs from that in deep ocean because of bottom interactions. In this paper, we focus on the merging of surface and bottom boundary layers in a finite-depth coastal ocean by numerically solving the wave-averaged equations using a large-eddy simulation method. The ocean fluid is driven by combined effects of wind stress, surface wave, and a steady current in the presence of stable vertical stratification. The resulting flow consists of two overlapping boundary layers, i.e., surface and bottom boundary layers, separated by an interior stratification. The overlapping boundary layers evolve through three phases, i.e., a rapid deepening, an oscillatory equilibrium and a prompt merger, separated by two transitions. Before the merger, internal waves are observed in the stratified layer, and they are excited mainly by Langmuir turbulence in the surface boundary layer. These waves induce a clear modulation on the bottom-generated turbulence, facilitating the interaction between the surface and bottom boundary layers. After the merger, the Langmuir circulations originally confined to the surface layer are found to grow in size and extend down to the sea bottom (even though the surface waves do not feel the bottom), reminiscent of the well-organized Langmuir supercells. These full-depth Langmuir circulations promote the vertical mixing and enhance the bottom shear, leading to a significant enhancement of turbulence levels in the vertical column. [ABSTRACT FROM AUTHOR]

We explore the decoupling physics of a stratocumulus-topped boundary layer (STBL) moving over cooler water, a situation mimicking warm-air advection (WADV). We simulate an initially well-mixed STBL over a doubly periodic domain with the sea surface temperature decreasing linearly over time using the System for Atmospheric Modeling large-eddy model. Due to the surface cooling, the STBL becomes increasingly stably stratified, manifested as a near-surface temperature inversion topped by a well-mixed cloud-containing layer. Unlike the stably stratified STBL in cold-air advection (CADV) that is characterized by cumulus coupling, the stratocumulus deck in the WADV is unambiguously decoupled from the sea surface, manifested as weakly negative buoyancy flux throughout the subcloud layer. Without the influxes of buoyancy from the surface, the convective circulation in the well-mixed cloud-containing layer is driven by cloud-top radiative cooling. In such a regime, the downdrafts propel the circulation, in contrast to that in CADV regime for which the cumulus updrafts play a more determinant role. Such a contrast in convection regime explains the difference in many aspects of the STBLs including the entrainment rate, cloud homogeneity, vertical exchanges of heat and moisture, and lifetime of the stratocumulus deck, with the last being subject to a more thorough investigation in Part II. Finally, we investigate under what conditions a secondary stratus near the surface (or fog) can form in the WADV. We found that weaker subsidence favors the formation of fog whereas a more rapid surface cooling rate does not. Significant Statement: The low-lying blanketlike clouds, called stratocumulus (Sc), reflect much incoming sunlight, substantially modulating Earth's temperature. While much is known about how the Sc evolves when it moves over warmer water, few studies examine the opposite situation of Sc moving over colder water. We used a high-resolution numerical model to simulate such a case. When moving over cold water, the Sc becomes unambiguously decoupled from the water surface, distinctive from its warm counterpart in which the Sc interacts with the water surface via intermittent cauliflower-like clouds called cumulus clouds. Such decoupling influences many aspects of the Sc–sea surface system, which combine to alter the ability of the Sc to reflect sunlight, thereby influencing the climate. This work laid the foundation for future work that quantifies the contribution of such a decoupled Sc regime to Earth's radiative budget and climate change. [ABSTRACT FROM AUTHOR]

This work explores the influence of weighted essentially nonoscillatory (WENO) schemes on Cloud Model 1 (CM1) large-eddy simulations (LES) of a quasi-steady, horizontally homogeneous, fully developed, neutral atmospheric boundary layer (ABL). An advantage of applying WENO schemes to scalar advection in compressible models is the elimination of acoustic waves and associated oscillations of domain-total vertical velocity. Applying WENO schemes to momentum advection in addition to scalar advection yields no further advantage but has an adverse effect on resolved turbulence within LES. As a tool designed to reduce numerically generated spurious oscillations, WENO schemes also suppress physically realistic instability development in turbulence-resolving simulations. Thus, applying WENO schemes to momentum advection reduces vortex stretching, suppresses the energy cascade, reduces shear-production of resolved Reynolds stress, and eventually amplifies the differences between the surface-layer mean wind profiles in the LES and the mean wind profiles expected in accordance with the filtered law of the wall (LOTW). The role of WENO schemes in adversely influencing surface-layer turbulence has inspired a concept of anti-WENO (AWENO) schemes to enhance instability development in regions where energy-containing turbulent motions are inadequately resolved by LES grids. The success in reproducing the filtered LOTW via AWENO schemes suggests that improving advection schemes is a critical component toward faithfully simulating near-surface turbulence and dealing with other "terra incognita" problems. Significance Statement: Turbulent motions are produced through instability development. Advection schemes designed to avoid generating spurious oscillations in flow fields involving sharp gradients may also suppress the development of physically realistic instabilities. This work explores the influence of advection schemes on numerical simulations that resolve turbulent motions. Recommendations are made concerning the use of advection schemes in simulating atmospheric turbulence, which is almost always accompanied with thermodynamic processes involving sharp temperature and moisture gradients. In addition, an advection-scheme-based concept is proposed to reproduce turbulent motions that are inadequately resolved by simulation grids. [ABSTRACT FROM AUTHOR]

Boundary layer turbulent processes affect tropical cyclone (TC) structure and intensity change. However, uncertainties in the parameterization of the planetary boundary layer (PBL) under high-wind conditions remain challenging, mostly due to limited observations. This study presents and evaluates a framework of numerical simulation that can be used for a small-domain [O(5)-km] large-eddy simulation (LES) and single-column modeling (SCM) to study the TC boundary layer. The framework builds upon a previous study that uses a few input parameters to represent the TC vortex and adds a simple nudging term for temperature and moisture to account for the complex thermodynamic processes in TCs. The reference thermodynamic profiles at different wind speeds are retrieved from a composite analysis of dropsonde observations of mature hurricanes. Results from LES show that most of the turbulence kinetic energy and vertical momentum flux is associated with resolved processes when horizontal grid spacing is O(10) m. Comparison to observations of turbulence variables such as momentum flux, effective eddy viscosity, and turbulence length scale show that LES produces reasonable results but highlight areas where further observations are necessary. LES results also demonstrate that compared to a classic Ekman-type boundary layer, the TC boundary layer is shallower, develops steady conditions much quicker, and exhibits stronger wind speed near the surface. The utility of this framework is further highlighted by evaluating a first-order PBL parameterization, suggesting that an asymptotic turbulence length scale of 40 m produces a good match to LES results. [ABSTRACT FROM AUTHOR]

Turbulent flow in a weakly convective marine atmospheric boundary layer (MABL) driven by geostrophic winds Vg = 10 m s−1 and heterogeneous sea surface temperature (SST) is examined using fine-mesh large-eddy simulation (LES). The imposed SST heterogeneity is a single-sided warm or cold front with jumps Δθ = (2, −1.5) K varying over a horizontal x distance of 1 km characteristic of an upper-ocean mesoscale or submesoscale front. The geostrophic winds are oriented parallel to the SST isotherms (i.e., the winds are alongfront). Previously, Sullivan et al. examined a similar flow configuration but with geostrophic winds oriented perpendicular to the imposed SST isotherms (i.e., the winds were across-front). Results with alongfront and across-front winds differ in important ways. With alongfront winds, the ageostrophic surface wind is weak, about 5 times smaller than the geostrophic wind, and horizontal pressure gradients couple the SST front and the atmosphere in the momentum budget. With across-front winds, horizontal pressure gradients are weak and mean horizontal advection primarily balances vertical flux divergence. Alongfront winds generate persistent secondary circulations (SC) that modify the surface fluxes as well as turbulent fluxes in the MABL interior depending on the sign of Δθ. Warm and cold filaments develop opposing pairs of SC with a central upwelling or downwelling region between the cells. Cold filaments reduce the entrainment near the boundary layer top that can potentially impact cloud initiation. The surface-wind–SST-isotherm orientation is an important component of atmosphere–ocean coupling. The results also show frontogenetic tendencies in the MABL. [ABSTRACT FROM AUTHOR]

Recently, large-eddy simulation (LES) has been increasingly employed in meteorological simulations because it is a promising method for turbulent parameterization. However, it is still difficult to affirm that the numerical accuracy required for a dynamical core is fully understood. In this study, we derived two theoretical criteria for the order of accuracy of the advection term in a typical situation of the atmospheric boundary layer, and demonstrate their validity by numerical experiments. In the targeted grid-spacing of O(10) m, we determined the required order of accuracy as follows: Based on the criterion of the numerical diffusion error, the upwind scheme must have at least seventh-order accuracy. The fourth-order central scheme is barely acceptable with fourth-order explicit diffusion, provided that its coefficient is one or two orders of magnitude smaller than the implicit diffusion coefficient of the third-order upwind scheme. Based on the criterion of numerical dispersion error, at minimum, the seventh or eighth order is required. The dispersion error was indirect for the energy spectra, although we expect it may affect the local turbulence mechanism. We also investigated the effects of temporal discretization for compressible models, and found that relatively lower-order time schemes are available up to the O(10) m grid spacing if the time step is sufficiently small due to sound wave limitations. The importance of the derived criteria is that the required order of accuracy increases as the grid spacing decreases. This suggests that considerable care should be taken regarding the numerical error problem for future high-resolution LES. [ABSTRACT FROM AUTHOR]

Abstract Numerous studies have shown that wind turbine wakes within a large wind farm bring about changes to both the dynamics and thermodynamics of the atmospheric boundary layers (ABL). Previously, we investigated the relative humidity budget within a wind farm via field measurements in the near‐wake region and large eddy simulations (LES). The effect of the compounding wakes within a large wind farm on the relative humidity was also investigated by LES. In this study, we investigate how the areas of relative humidity variation, that was observed in the near‐wake, develop downstream in the shadow region of a large wind farm. To this end, LES of a wind farm consisting of 8x6 wind turbines with periodic boundary condition in the lateral direction (inferring an infinitely wide farm) interacting with a stable ABL is carried out. Two wind farm layouts, aligned and staggered, are considered in the analysis and the results from both configurations are compared to each other. It is observed that a decrease of relative humidity underneath the hub height and an increase above the hub height build up within the wind farm, and are maintained in the downstream of the farm for long distances. The staggered farm layout is more effective in keeping a more elongated region of low relative humidity underneath the hub, when compared to the aligned layout.

The wind direction in atmospheric boundary layers changes continuously due to meso-scale weather phenomena. Developing accurate simulations of these changes is essential for understanding their effect on the performance of large wind farms. Our study introduces a new technique to model dynamic wind direction changes obtained from meso-scale simulations or field measurements in micro-scale large eddy simulations. We propose a method in which the simulation domain is treated as a non-inertial rotating reference frame. The primary benefit of our approach is that it is straightforward to implement and reproduces desired wind direction changes excellently. We verified our approach in neutral atmospheric boundary layers and show that the observed boundary-layer characteristics for dynamic wind directions agree very well with those observed for constant mean wind directions when the wind direction is changed slowly such that the flow is quasi-stationary. Further, we show that atmospheric measurements of the wind direction can be reproduced by our method. To underline the importance of the method, we conclude with a representative scenario, which shows that dynamic wind direction changes can affect the performance of large wind farms. • New technique to model dynamic wind direction changes in large eddy simulations. • Approach is validated for simulations of neutral atmospheric boundary layers. • Excellent agreement between imposed and simulated time-varying wind directions. • Performance of large wind farms strongly affected by changing wind directions. [ABSTRACT FROM AUTHOR]

This study explores the response of flow around isolated cuboid buildings to variations in the incoming turbulence arising from changes in atmospheric boundary layer (ABL) stability using a building-resolving large-eddy simulation (LES) technique with explicit representation of building effects through an immersed body force method. An extensive suite of LES for a neutral ABL with different model resolution and advection scheme configurations reveals that at least 6, 12, and 24 grid points per building side are required in order to resolve building-induced vortex shedding, mean-flow features, and turbulence statistics, respectively, with an advection scheme of a minimum of third order. Using model resolutions that meet this requirement, 21 building-resolving simulations are performed under varying atmospheric stability conditions, from weakly stable to convective ABLs, and for different building sizes (H), resulting in LABL/H ~ 0.1-10, where LABL is the integral length scale of the incoming ABL turbulence. The building-induced flow features observed in the canonical neutral ABL simulation, e.g., the upstream horseshoe vortex and the downstream arch vortex, gradually weaken with increasing surface-driven convective instability due to the enhancement of background turbulent mixing. As a result, two local turbulence kinetic energy peaks on the lateral side of the building in nonconvective cases are merged into a single peak in strong convective cases. By considering the ABL turbulence scale and building size altogether, it is shown that the building impact decreases with increasing LABLH, as coherent turbulent structures in the ABL become more dominant over a building-induced flow response for LABL/H > 1. [ABSTRACT FROM AUTHOR]

The diurnally evolving trapped lee wave over a small-scale two-dimensional steep mountain is investigated in large-eddy simulations based on a fully compressible and nonhydrostatic model [Icosahedral Nonhydrostatic (ICON)] with triangular grids of 50-m-edge length. An idealized atmospheric profile derived from a realistic case is designed to account for influences from the stagnant layer near the surface, the stability of the atmospheric boundary layer (ABL) and the upper-level jet. First, simulations were done to bridge from the linear regime to the nonlinear regime by increasing the mountain height, which showed that larger-amplitude lee waves with longer wavelength can be produced in the nonlinear regime than in the linear regime. Second, the effects of the stagnant layer near the surface and the ABL stability were explored, which showed that the stagnant layer or the stable ABL can play a similar wave-absorbing role in the nonlinear regime as in linear theories or simulations. Third, the role of the upper-level jet was explored, indicating that a stronger (weaker) upper-level jet can help to produce longer (shorter) lee waves. The stable ABL with a stagnant layer can more (less) efficiently absorb the longer (shorter) lee waves due to the stronger (weaker) jet, so that the wave response is more sensitive to the wave-absorption layer when an upper-level jet is present. Finally, the momentum budget was analyzed to explore the interaction between the upper and lower levels of the troposphere, which showed that the momentum flux due to the upward-propagating waves and trapped waves varies with the upper-level jet strength and low-level stagnancy and ABL stability. [ABSTRACT FROM AUTHOR]

Eddy diffusivity models are a common method to parameterize turbulent fluxes in the atmospheric sciences community. However, their inability to handle convective boundary layers leads to the addition of a nondiffusive flux component (usually called nonlocal) alongside the original diffusive term (usually called local). Both components are often modeled for convective conditions based on the shape of the eddy diffusivity profile for neutral conditions. This assumption of shape is traditionally employed due to the difficulty of estimating both components based on numerically simulated turbulent fluxes without any a priori assumptions. In this manuscript we propose a novel method to avoid this issue and estimate both components from numerical simulations without having to assume any a priori shape or scaling for either. Our approach is based on optimizing results from a modeling perspective and taking as much advantage as possible from the diffusive term, thus maximizing the eddy diffusivity. We use our method to diagnostically investigate four different large-eddy simulations spanning different stability regimes, which reveal that nondiffusive fluxes are important even when trying to minimize them. Furthermore, the calculated profiles for both diffusive and nondiffusive fluxes suggest that their shapes change with stability, which is an effect that is not included in most models currently in use. Finally, we use our results to discuss modeling approaches and identify opportunities for improving current models. [ABSTRACT FROM AUTHOR]

A three-dimensional, prognostic Atmospheric Boundary Layer Environment–Lattice Boltzmann Model (ABLE-LBM) using the multiple-relaxation-time lattice Boltzmann method was developed for large-eddy simulation of urban boundary layer atmospheric flows. In this article we describe the details of the ABLE-LBM for urban flow, its implementation of complex boundaries, and the subgrid turbulence parameterizations. As a first validation of this newly developed model, the simulation results were evaluated with two wind-tunnel datasets that were collected using particle image velocimetry and Irwin probes, respectively. The ABLE-LBM simulations use the same building layout and Reynolds numbers used in the laboratory wind tunnels. The ABLE-LBM simulations compare favorably to both laboratory studies in terms of the mean wind fields. The turbulent fluxes simulated by the model in the observational planes also agreed reasonably well with the laboratory results. The model produced urban canyon flows and vortices on the lee side and over the building tops that are similar to those of the laboratory studies in strength and location. This validation study using laboratory data indicates that our new ABLE-LBM is a viable approach for modeling atmospheric turbulent flows in urban environments. A numerical implementation using a graphics processing unit shows that real-time simulations are achieved for these two validation cases. [ABSTRACT FROM AUTHOR]

This is the Part II of a two-part study that seeks a theoretical understanding of an empirical relationship for shallow cumulus clouds: subcloud updraft velocity covaries linearly with the cloud-base height. This work focuses on continental cumulus clouds that are more strongly forced by surface fluxes and more deviated from equilibrium than those over oceans (Part I). We use a simple analytical model for shallow cumulus that is well tested against a high-resolution (25 m in the horizontal) large-eddy simulation model. Consistent with a conventional idea, we find that surface Bowen ratio is the key variable that regulates the covariability of both parameters: under the same solar insolation, a drier surface allows for stronger buoyancy flux, triggering stronger convection that deepens the subcloud layer. We find that the slope of the Bowen-ratio-regulated relationship between the two parameters (defined as λ) is dependent on both the local time and the stability of the lower free atmosphere. The value of λ decreases with time exponentially from sunrise to early afternoon and linearly from early afternoon to sunset. The value of λ is larger in a more stable atmosphere. In addition, continental λ in the early afternoon more than doubles the oceanic λ. Validation of the theoretical results against ground observations over the Southern Great Plains shows a reasonable agreement. Physical mechanisms underlying the findings are explained from the perspective of different time scales at which updrafts and cloud-base height respond to a surface flux forcing. [ABSTRACT FROM AUTHOR]

Numerous studies have shown that wind turbine wakes within a large wind farm bring about changes to both the dynamics and thermodynamics of the atmospheric boundary layers (ABL). Previously, we investigated the relative humidity budget within a wind farm via field measurements in the near‐wake region and large eddy simulations (LES). The effect of the compounding wakes within a large wind farm on the relative humidity was also investigated by LES. In this study, we investigate how the areas of relative humidity variation, that was observed in the near‐wake, develop downstream in the shadow region of a large wind farm. To this end, LES of a wind farm consisting of 8x6 wind turbines with periodic boundary condition in the lateral direction (inferring an infinitely wide farm) interacting with a stable ABL is carried out. Two wind farm layouts, aligned and staggered, are considered in the analysis and the results from both configurations are compared to each other. It is observed that a decrease of relative humidity underneath the hub height and an increase above the hub height build up within the wind farm, and are maintained in the downstream of the farm for long distances. The staggered farm layout is more effective in keeping a more elongated region of low relative humidity underneath the hub, when compared to the aligned layout. [ABSTRACT FROM AUTHOR]

A realizable k–ε turbulence model of incompressible fluid is extended for the stable atmosphere after taking account of the buoyancy damping of gravity waves. The new model is consistent with the Monin–Obukhov similarity theory on the stable atmospheric boundary layer (ABL) over a horizontally uniform surface. The model is incorporated into an ABL model to simulate mean flow against observations. Its ABL-model output is compared with the Leipzig dataset, showing the turbulence model works well for a stable ABL. Specifically, the ABL model properly replicates 1) the mixing length, turbulent viscosity, and mean wind; 2) a significant decrease of the mixing length with height in the upper ABL and thus a reasonable altitude of the ABL top; and 3) a sensitivity of the mixing length and turbulent viscosity to atmospheric stability. [ABSTRACT FROM AUTHOR]

The Weather Research and Forecasting (WRF) Model is used in large-eddy simulation (LES) mode to investigate a lake-breeze case occurring on 12 June 2012 over the Lake Taihu region of China. Observational data from 15 locations, wind profiler radar, and the Moderate Resolution Imaging Spectroradiometer (MODIS) are used to evaluate the WRF nested-LES performance in simulating lake breezes. Results indicate that the simulated temporal and spatial variations of the lake breeze by WRF nested LES are consistent with observations. The simulations with high-resolution grid spacing and the LES scheme have a high correlation coefficient and low mean bias when evaluated against 2-m temperature, 10-m wind, and horizontal and vertical lake-breeze circulations. The atmospheric boundary layer (ABL) remains stable over the lake throughout the lake-breeze event, and the stability becomes even stronger as the lake breeze reaches its mature stage. The improved ABL simulation with LES at a grid spacing of 150 m indicates that the non-LES planetary boundary layer parameterization scheme does not adequately represent subgrid-scale turbulent motions. Running WRF fully coupled to a lake model improves lake-surface temperature and consequently the lake-breeze simulations. Allowing for additional model spinup results in a positive impact on lake-surface temperature prediction but is a heavy computational burden. Refinement of a water-property parameter used in the Community Land Model, version 4.5, within WRF and constraining the lake-surface temperature with observational data would further improve lake-breeze representation. [ABSTRACT FROM AUTHOR]

As model grid resolutions move from the mesoscale to the microscale, turbulent structures represented in atmospheric boundary layer simulations change dramatically. At intermediate resolutions, the so-called gray zone, turbulent motions are not resolved accurately, posing a challenge to numerical simulations. The representation of turbulence is also highly sensitive to the choice of closure model. Here, we examine explicit filtering and reconstruction in the gray zone as a technique to better represent atmospheric turbulence. The convective boundary layer is simulated using the Weather Research and Forecasting (WRF) Model with horizontal resolutions ranging from 25 m to 1 km. Four large-eddy simulation (LES) turbulence models are considered: the Smagorinsky model, the TKE-1.5 model, and two versions of the dynamic reconstruction model (DRM). The models are evaluated by their ability to produce consistent mean potential temperature profiles, heat and momentum fluxes, velocity fields, and turbulent kinetic energy spectra as the grids become coarser. The DRM, a mixed model that uses an explicit filtering and reconstruction technique to account for resolvable subfilter-scale (RSFS) stresses, performs very well at resolutions of 500 m and 1 km without any special tuning, whereas the Smagorinsky and TKE-1.5 models produce heavily grid-dependent results. [ABSTRACT FROM AUTHOR]

A canopy model framework is implemented in the Weather Research and Forecasting Model to improve the accuracy of large-eddy simulations (LES) of the atmospheric boundary layer (ABL). The model includes two options that depend on the scale of surface roughness elements. A resolved canopy model, typically used to model flow through vegetation canopies, is employed when roughness elements are resolved by the vertical LES grid. In the case of unresolved roughness, a modified "pseudocanopy model" is developed to distribute drag over a shallow layer above the surface. Both canopy model options are validated against idealized test cases in neutral stability conditions and are shown to improve surface layer velocity profiles relative to simulations employing Monin–Obukhov similarity theory (MOST), which is commonly used as a surface boundary condition in ABL models. Use of the canopy model framework also leads to increased levels of resolved turbulence kinetic energy and turbulent stresses. Because LES of the ABL has a well-known difficulty recovering the expected logarithmic velocity profile (log law) in the surface layer, particular focus is placed on using the pseudocanopy model to alleviate this issue over a range of model configurations. Tests with varying surface roughness values, LES closures, and grid aspect ratios confirm that the pseudocanopy model generally improves log-law agreement relative to simulations that employ a standard MOST boundary condition. The canopy model framework thus represents a low-cost, easy-to-implement method for improving LES of the ABL. [ABSTRACT FROM AUTHOR]

The vertical velocity probability distribution function (PDF) is analyzed throughout the depth of the lower atmosphere, including the subcloud and cloud layers, in four large-eddy simulation (LES) cases of shallow cumulus and stratocumulus. Double-Gaussian PDF closures are examined to test their ability to represent a wide range of turbulence statistics, from stratocumulus cloud layers characterized by Gaussian turbulence to shallow cumulus cloud layers displaying strongly non-Gaussian turbulence statistics. While the majority of the model closures are found to perform well in the former case, the latter presents a considerable challenge. A new model closure is suggested that accounts for high skewness and kurtosis seen in shallow cumulus cloud layers. The well-established parabolic relationship between skewness and kurtosis is examined, with results in agreement with previous studies for the subcloud layer. In cumulus cloud layers, however, a modified relationship is necessary to improve performance. The new closure significantly improves the estimation of the vertical velocity PDF for shallow cumulus cloud layers, in addition to performing well for stratocumulus. In particular, the long updraft tail representing the bulk of cloudy points is much better represented and higher-order moments diagnosed from the PDF are also greatly improved. However, some deficiencies remain owing to fundamental limitations of representing highly non-Gaussian turbulence statistics with a double-Gaussian PDF. [ABSTRACT FROM AUTHOR]

This paper investigates the effects of baroclinic pressure gradients on mean flow and turbulence in the diabatic atmospheric boundary layer (ABL). Large-eddy simulations are conducted where the direction of the baroclinicity, its strength, and the surface buoyancy flux are systematically varied to examine their interacting effects. The thermal wind vector, which represents the vertical change in the geostrophic wind vector resulting from horizontal temperature gradients, significantly influences the velocity profiles, the Ekman turning, and the strength and location of the low-level jet (LLJ). For instance, cold advection and positive (negative) geostrophic shear increased (decreased) friction velocity and changed the LLJ elevation. Given the baroclinicity strength and direction under neutral conditions, a simple reduced model is proposed and validated here to predict the general trends of baroclinic mean winds. The baroclinic effects on turbulence intensity and structure are even more intricate, with turbulent kinetic energy (TKE) profiles displaying an increase of TKE magnitude with height for some cases. When a fixed destabilizing surface heat flux is added, a positive geostrophic shear favors streamwise aligned roll-type structures, which are typical of neutral ABLs. Conversely, a negative geostrophic shear promotes cell-type structures, which are typical of strongly unstable ABLs. Furthermore, baroclinicity increases shear in the outer ABL and tends to make the outer flow more neutral by decreasing the Richardson flux number. These findings are consequential for meteorological measurements and the wind-energy industry, among others: baroclinicity alters the mean wind profiles, the TKE, coherent structures, and the stability of the ABL, and its effects need to be considered. [ABSTRACT FROM AUTHOR]

A new turbulence scheme with two prognostic energies is presented. The scheme is an extension of a turbulence kinetic energy (TKE) scheme following the ideas of Zilitinkevich et al. but valid for the whole stability range and including the influence of moisture. The second turbulence prognostic energy is used only for a modification of the stability parameter. Thus, the scheme is downgradient, and the turbulent fluxes are proportional to the local gradients of the diffused variables. However, the stability parameter and consequently the turbulent exchange coefficients are not strictly local anymore and have a prognostic character. The authors believe that these characteristics enable the scheme to model both turbulence and clouds in the planetary boundary layer. The two-energy scheme was tested in three idealized single-column model (SCM) simulations, two in the convective boundary layer and one in the stable boundary layer. Overall, the scheme performs better than the standard TKE schemes. Compared to the TKE schemes, the two-energy scheme shows a more continuous behavior in time and space and mixes deeper in accordance with the LES results. A drawback of the scheme is that the modeled thermals tend to be too intense and too infrequent. This is due to the particular cutoff formulation of the chosen length-scale parameterization. Long-term three-dimensional global simulations show that the turbulence scheme behaves reasonable well in a full atmospheric model. In agreement with the SCM simulations, the scheme tends to overestimate cloud cover, especially at low levels. [ABSTRACT FROM AUTHOR]

Computational fluid dynamics predictions of urban flow are subject to several sources of uncertainty, such as the definition of the inflow boundary conditions or the turbulence model. Compared to Reynolds-averaged Navier-Stokes (RANS) simulations, large eddy simulations (LES) can reduce turbulence model uncertainty by resolving the turbulence down to scales in the inertial subrange, but the presence of other uncertainties will not be reduced. The objective of this study is to present an initial investigation of the relative importance of these different types of uncertainties by comparing urban flow predictions obtained using RANS and LES to field measurements. The simulations are designed to reproduce measurements performed during the Joint Urban 2003 field experiments. The time-averaged velocity measured at an upstream wind sensor is used to define the inflow boundary condition, and the results are compared to time-averaged measurements at 34 locations in the downtown area. For the turbulence kinetic energy, the LES is found to be more accurate than the RANS in 80% of the available high-frequency measurement locations. For the mean velocity field, this number reduces to 50% of all stations. Comparison of the LES results with a previous inflow uncertainty quantification study for RANS shows that locations where the LES is less accurate than the RANS correspond to locations where the RANS solution is highly sensitive to the inflow boundary conditions. This suggests that inflow uncertainties can be a dominant factor, and that their effect on LES results should be quantified to guarantee predictive capabilities. [ABSTRACT FROM AUTHOR]

The characteristics of wind profiles in a neutral atmospheric boundary layer and their dependence on the geostrophic wind speed Ug, Coriolis parameter f, and surface roughness length z0 are examined utilizing large-eddy simulations. These simulations produce a constant momentum flux layer and a log-law layer above the surface characterized by a logarithmic increase of wind speed with height. The von Kármán constant derived from the mean wind profile is around 0.4 over a wide range of control parameters. The depths of the simulated boundary layer, constant-flux layer, and surface log-law layer tend to increase with the wind speed and decrease with an increasing Coriolis parameter. Immediately above the surface log-law layer, a second log-law layer has been identified from these simulations. The depth of this upper log-law layer is comparable to its counterpart in the surface layer, and the wind speed can be scaled as , as opposed to just in the surface log-law layer, implying that in addition to surface processes, the upper log-law layer is also influenced by Earth’s rotation and large-scale conditions. Here is the friction velocity at the surface, and h is the boundary layer depth. An analytical model is proposed to assist in the interpretation of the log laws in a typical Ekman boundary layer. The physics and implications of the upper log-law layer are discussed. [ABSTRACT FROM AUTHOR]

The stratocumulus cloud–capped boundary layer under a sharp inversion is a challenging regime for large-eddy simulation (LES). Here, data from the first research flight of the Second Dynamics and Chemistry of the Marine Stratocumulus field study are used to evaluate the effect of different LES turbulence closures. Six different turbulence models, including traditional TKE and Smagorinsky models and more advanced models that employ explicit filtering and reconstruction, are tested. The traditional models produce unrealistically thin clouds and a decoupled boundary layer as compared with other more advanced models. Traditional models rely on specified subfilter-scale (SFS) Prandtl and Schmidt numbers to obtain SFS eddy diffusivity from eddy viscosity, whereas dynamic models can compute SFS eddy diffusivity independently through dynamic procedures. The effective SFS Prandtl number in dynamic models is found to be ~0.5 below the cloud and ~10 inside the cloud layer, implying minimized mixing in the cloud. In contrast, the SFS Prandtl number in traditional models is about 1/3 throughout the entire boundary layer, suggesting spuriously strong mixing in the cloud. The SFS Schmidt number in the dynamic models also changes independently from the SFS Prandtl number, whereas in traditional models they are identical, meaning that the efficiency of the turbulent mixing of water content is forced to be the same as that of heat. Since it is very difficult to know in advance the SFS Prandtl and Schmidt numbers in a given flow, dynamic models may provide a more realistic representation of scalar mixing in LES. [ABSTRACT FROM AUTHOR]

Large-eddy simulation (LES) has been an essential tool in the development of theory and parameterizations for clouds, but when applied to stratocumulus clouds under sharp temperature inversions, many LES models produce an unrealistically thin cloud layer and a decoupled boundary layer structure. Here, explicit filtering and reconstruction are used for simulation of stratocumulus clouds observed during the first research flight (RF01) of the Second Dynamics and Chemistry of the Marine Stratocumulus field study (DYCOMS II). The dynamic reconstruction model (DRM) is used within an explicit filtering and reconstruction framework, partitioning subfilter-scale motions into resolvable subfilter scales (RSFSs) and unresolvable subgrid scales (SGSs). The former are reconstructed, and the latter are modeled. Differing from traditional turbulence models, the reconstructed RSFS stress/fluxes can produce backscatter of turbulence kinetic energy (TKE) and, importantly, turbulence potential energy (TPE). The modeled backscatter reduces entrainment at the cloud top and, meanwhile, strengthens resolved turbulence through preserving TKE and TPE, resulting in a realistic boundary layer with an adequate amount of cloud water and vertically coupled turbulent eddies. Additional simulations are performed in the terra incognita, when the grid spacing of a simulation becomes comparable to the size of the most energetic eddies. With 20-m vertical and 1-km horizontal grid spacings, simulations using DRM provide a reasonable representation of bulk properties of the stratocumulus-capped boundary layer. [ABSTRACT FROM AUTHOR]

The role of surface friction on shallow nonprecipitating convection is investigated using a series of large-eddy simulations with varying surface friction velocity and with a cloud identification algorithm. As surface friction intensifies, convective rolls dominate over convective cells and secondary overturning circulation becomes stronger in the subcloud layer, thus transporting more moisture upward and more heat downward between the subcloud and cloud layers. Identifying individual clouds, using the identification algorithm based on a three-dimensional topological analysis, reveals that intensified surface friction increases the number of clouds and the degree of tilting in the downstream direction. Highly intensified surface friction increases wind shear across the cloud base and induces cloud tilting, which leads to a vertically parabolic profile of liquid water mixing ratio instead of the classical two-layer structure (conditionally unstable and trade inversion layers). Furthermore, cloud tilting induces more cloud cover andmore cloud mass fluxmuch above the cloud base (e.g., 0.8 < z < 1.2 km), but less cloud cover and less cloud mass flux in the upper cloud layer (e.g., z > 1.2 km) because of increased lateral entrainment rate. Similarly, profiles of directly measured entrainment and detrainment rates show that detrainment in the lower cloud layer becomes smaller with stronger surface friction. [ABSTRACT FROM AUTHOR]

This paper presents large eddy simulations (LES) of model-scaled neutrally stratified atmospheric boundary layer (ABL) flows for structural engineering applications and examines their statistical properties. A one-k-equation eddy model is used for the subgrid-scale (SGS) motions, and a wall shear model is applied on the ground. The mean streamwise velocity profile is approximately logarithmic, yet near the ground a mismatch persists because of the limited accuracy of the SGS model. The second moments of the turbulence represent well the underlying physics. The spectra of the velocity components at a point are consistent with commonly accepted expressions. The spatial coherences decay exponentially as functions of reduced frequencies. The results suggest that, except for the mean velocity near the ground, the ABL turbulence statistics can be well represented by large eddy simulations with simple SGS and wall models. [ABSTRACT FROM AUTHOR]

Compared to the representation of vertical turbulent mixing through various planetary boundary layer (PBL) schemes, the treatment of horizontal turbulent mixing in the boundary layer has received much less attention. In mesoscale and convective-scale models, subgrid-scale horizontal turbulent mixing has traditionally been associated with mesoscale circulations or eddies. Its parameterization most often adopts the gradient-diffusion model, where the horizontal mixing coefficients are usually set constant, or through the 2D Smagorinsky formulation, or in some cases based on the 1.5-order turbulence kinetic energy (TKE) closure. For horizontal turbulent mixing associated with boundary layer eddies, the traditional schemes are shown to perform poorly. This work investigates the characteristic turbulence velocity and length scales based on analysis of a well-resolved, wide-domain large-eddy simulation of a convective boundary layer (CBL). To improve the representation of horizontal turbulent mixing by CBL eddies, a class of schemes is proposed with different levels of sophistication. The first two schemes can be used together with first-order PBL schemes, while the third uses TKE to define its characteristic velocity scale and can be used together with TKE-based higher-order PBL schemes. The proposed parameterizations are tested a posteriori in idealized simulations of turbulent dispersion of a passive scalar. Comparisons show improved horizontal dispersion by the proposed schemes and further demonstrate the weakness of the existing schemes. [ABSTRACT FROM AUTHOR]

While it has long been a practice to place spires near the inlet of a wind tunnel to quickly develop a turbulent boundary layer with similarities to an atmospheric boundary layer, this has not been the case for creating turbulent boundary layer inflow in large eddy simulations (LESs) of turbulent flows. We carry out LES with the curvilinear immersed boundary method to simulate the flow in a wind tunnel with a series of spires in order to investigate the feasibility of numerically developing inflow conditions from a precursory spire LES and assessing the similarities of the turbulence statistics to those of an atmospheric boundary layer. The simulated mean velocity field demonstrates that a turbulent boundary layer with height equal to the spire height develops very quickly, within five spire heights downstream. The major attribute of using spires for precursory simulations is the spatially evolving coherent structures that form downstream of the spires offering a range of length scales at both the vertical and streamwise directions allowing multiple turbulent inflow conditions to be extracted from a single simulation. While the distribution of length scales far from the spires resembles an atmospheric boundary layer, some turbulence statistics have some significant differences. [ABSTRACT FROM AUTHOR]

The ability to make forecasts depends on atmospheric predictability and the growth of errors. It has recently been shown that the predictability of urban boundary layers differs in important respects from that of the free atmosphere on the mesoscale and larger; in particular, nonlinearity may play a less prominent role in the error evolution. This paper investigates the applicability of linear theory to the error evolution in turbulent street-canyon flow. Using large-eddy simulation, streamwise aspect ratios between 0.15 and 1.50, and identical-twin experiments, it is shown that the growth rate of the error kinetic energy can be estimated from Eulerian averages and that linear theory provides insight into the spatial structure of the error field after saturation. The results should be applicable to cities with deep and closely spaced canyons. Implications for data assimilation and modeling are discussed. [ABSTRACT FROM AUTHOR]

The atmospheric boundary layer (ABL) inflow has significant effects on behaviors of floating offshore wind turbine (FOWT), especially for large-size wind turbine. In this study, numerical investigation of aerodynamics and wakes of a semi-submersible FOWT under ABL inflow is performed. The quasi-equilibrium ABL wind field is generated by large eddy simulations (LES) with sufficient simulation duration. The FOWT wakes are modeled by incorporating LES and actuator line model (ALM) in the computational fluid dynamics (CFD) framework, and the FOWT dynamic responses are simulated by FAST code. A two-way coupling procedure is employed, in which the wind velocity around wind turbine sampled in CFD framework and the wind turbine's body forces and positions solved by FAST code are delivered to each other. The simulation case of a bottom-fixed wind turbine is performed to provide some comparable data. It is revealed that the power variation of FOWT is dominated by atmospheric turbulence, more than platform motions. The slightly enhanced out-of-plane shear force and bending moment of FOWT are caused by platform motions. Owing to the entrance of ambient atmospheric flow into wind turbine wakes, significant deflection of wakes is visualized. In addition, the wake center of FOWT is far away from hub height level due to pitch motion of platform, which is a potential benefit factor for downstream wind turbines. The spatiotemporal characteristics of turbulence intensity in FOWT wakes are complex, and the discrepancies of wakes between floating and bottom-fixed scenarios are not significant. • More realistic ABL inflow is generated by LES, rather than synthetic turbulence generation models. • Coupling of ABL model and dynamics of FOWT is employed to investigate aerodynamics and wakes of FOWT subjected to LES ABL. • Case of bottom-fixed wind turbine is simulated for reference and comparison. • Power variation of FOWT is dominated by atmospheric turbulence, more than platform motions. [ABSTRACT FROM AUTHOR]

The impact of atmospheric boundary layer (ABL) interactions with large-scale stably stratified flow over an isolated, two-dimensional hill is investigated using turbulence-resolving large-eddy simulations. The onset of internal gravity wave breaking and leeside flow response regimes of trapped lee waves and nonlinear breakdown (or hydraulic-jump-like state) as they depend on the classical inverse Froude number, Fr−1 = Nh/ U g, is explored in detail. Here, N is the Brunt-Väisälä frequency, h is the hill height, and U g is the geostrophic wind. The results here demonstrate that the presence of a turbulent ABL influences mountain wave (MW) development in critical aspects, such as dissipation of trapped lee waves and amplified stagnation zone turbulence through Kelvin-Helmholtz instability. It is shown that the nature of interactions between the large-scale flow and the ABL is better characterized by a proposed inverse compensated Froude number, = N( h − z i)/ U g, where z i is the ABL height. In addition, it is found that the onset of the nonlinear-breakdown regime, ≈ 1.0, is initiated when the vertical wavelength becomes comparable to the sufficiently energetic scales of turbulence in the stagnation zone and ABL, yielding an abrupt change in leeside flow response. Finally, energy spectra are presented in the context of MW flows, supporting the existence of a clear transition in leeside flow response, and illustrating two distinct energy distribution states for the trapped-lee-wave and the nonlinear-breakdown regimes. [ABSTRACT FROM AUTHOR]

Mesoscale models are considered to be the state of the art in modeling mountain-wave flows. Herein, the authors investigate the role and accuracy of planetary boundary layer (PBL) parameterizations in handling the interaction between large-scale mountain waves and the atmospheric boundary layer. To that end, recent large-eddy simulation (LES) results of mountain waves over a symmetric two-dimensional bell-shaped hill are used and compared to four commonly used PBL schemes. It is found that one-dimensional PBL parameterizations produce reasonable agreement with the LES results in terms of vertical wavelength, amplitude of velocity, and turbulent kinetic energy distribution in the downhill shooting-flow region. However, the assumption of horizontal homogeneity in PBL parameterizations does not hold in the context of these complex flow configurations. This inappropriate modeling assumption results in a vertical wavelength shift, producing errors of approximately 10 m s−1 at downstream locations because of the presence of a coherent trapped lee wave that does not mix with the atmospheric boundary layer. In contrast, horizontally integrated momentum flux derived from these PBL schemes displays a realistic pattern. Therefore, results from mesoscale models using ensembles of one-dimensional PBL schemes can still potentially be used to parameterize drag effects in general circulation models. Nonetheless, three-dimensional PBL schemes must be developed in order for mesoscale models to accurately represent complex terrain and other types of flows where one-dimensional PBL assumptions are violated. [ABSTRACT FROM AUTHOR]

In many engineering and meteorological applications, atmospheric turbulence within the planetary boundary layer is described in terms of its representative parameters. One such parameter is the structure-function (or structure) parameter that is used to characterize the intensity of turbulent fluctuations of atmospheric flow variables. Structure parameters are derivatives of structure functions, but are used more frequently than the latter ones for practical needs as they do not explicitly include dependence on the separation distance. The structure parameter of potential temperature, which is the subject of this study, describes the spatial variability of the temperature fluctuations. It is broadly represented in theories and models of electromagnetic and acoustic wave propagation in the atmosphere, and forms the basis for the scintillometer measurement concept. The authors consider three methods to compute the potential temperature structure function and structure parameter: the direct method, the true spectral method, and the conventional spectral method. Each method is tested on high-resolution potential temperature datasets generated from large-eddy simulations of a variety of convective boundary layer flow cases reproduced by two representative numerical codes. Results indicate that the popular conventional spectral method routinely exaggerates the potential temperature structure-function parameter, likely due to the unrealistic assumptions underlying the method. The direct method and true spectral method are recommended as the more suitable approaches. [ABSTRACT FROM AUTHOR]

Large-eddy simulations of an idealized diurnal urban heat island are performed using the Weather Research and Forecasting Model. The surface energy balance over an inhomogeneous terrain is solved considering the anthropogenic heat contribution and the differences of thermal and mechanical properties between urban and rural surfaces. Several cases are simulated together with a reference case, considering different values of the control parameters: albedo, thermal inertia, roughness length, anthropogenic heat emission, and geostrophic wind intensity. Spatial distributions of second-moment statistics, including the turbulent kinetic energy (TKE) budget, are analyzed to characterize the structure of the planetary boundary layer (PBL). The effect of each control parameter value on the turbulent properties of the PBL is investigated with respect to the reference case. For all of the analyzed cases, the primary source of TKE is the buoyancy in the lower half of the PBL, the shear in the upper half, and the turbulent transport term at the top. The vertical advection of TKE is significant in the upper half of the PBL. The control parameters significantly influence the shape of the profiles of the transport and shear terms in the TKE budget. Bulk properties of the PBL via proper scaling are compared with literature data. A log-linear relationship between the aspect ratio of the heat island and the Froude number is confirmed. For the first time, the effect of relevant surface control parameters and the geostrophic wind intensity on the bulk and turbulent properties of the PBL is systematically investigated at high resolution. [ABSTRACT FROM AUTHOR]