Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δx) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δx ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation. Significance Statement: The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.

Wind energy applications including wind resource assessment, wind power forecasting, and wind plant optimization require high-resolution mesoscale simulations. High resolution mesoscale simulations are essential for accurate characterization of atmospheric flows over heterogeneous land use and complex terrain. Under such conditions, the assumption of grid-cell homogeneity, used in one-dimensional planetary boundary layer (1D PBL) parameterizations, breaks down. However, in most numerical weather prediction (NWP) models, boundary layer turbulence is parameterized using 1D PBL parameterizations. We have therefore developed a three-dimensional (3D) PBL parameterization to better account for horizontal flow heterogeneities. We have implemented and tested the 3D PBL parameterization in the Weather Research and Forecasting (WRF) numerical weather prediction model. The new parameterization is validated using observations from the Wind Forecast Improvement 2 (WFIP 2) project and compared to 1D PBL results.

The Mesoscale to Microscale Coupling team, part of the U.S. Department of Energy Atmosphere to Electrons (A2e) initiative, has studied various important challenges related to coupling mesoscale models to microscale models for the use case of wind energy development and operation. Several coupling methods and techniques for generating turbulence at the microscale that is subgrid to the mesoscale have been evaluated for a variety of cases. Case studies included flat-terrain, complex-terrain, and offshore environments. Methods were developed to bridge the terra incognita, which scales from about 100 m through the depth of the boundary layer. The team used wind-relevant metrics and archived code, case information, and assessment tools and is making those widely available. Lessons learned and discerned best practices are described in the context of the cases studied for the purpose of enabling further deployment of wind energy.

Wind plant wake impacts can be estimated with a number of simulation methodologies, each with its own fidelity and sensitivity to model inputs. In turbine-free mesoscale simulations, hub-height wind speeds often significantly vary with the choice of a planetary boundary layer (PBL) scheme. However, the sensitivity of wind plant wakes to a PBL scheme has not been explored because, as of the Weather Research and Forecasting model v4.3.3, wake parameterizations were only compatible with one PBL scheme. We couple the Fitch wind farm parameterization with the new NCAR 3DPBL scheme and compare the resulting wakes to those simulated with a widely used PBL scheme. We simulate a wind plant in pseudo-steady states under idealized stable, neutral, and unstable conditions with matching hub-height wind speeds using two PBL schemes: MYNN and the NCAR 3DPBL. For these idealized scenarios, average hub-height wind speed losses within the plant differ between PBL schemes by between −0.20 and 0.22 m s⁻¹, and correspondingly, capacity factors range between 39.5 %–53.8 %. These simulations suggest that PBL schemes represent a meaningful source of modeled wind resource uncertainty; therefore, we recommend incorporating PBL variability into future wind plant planning sensitivity studies as well as wind forecasting studies.

The 2020 wildfire season (May through December) in the United States was exceptionally active, with the National Interagency Fire Center reporting over 10 million acres ( > 40 000 km²) burned. During the September 2020 wildfire events, large concentrations of smoke particulates were emitted into the atmosphere. As a result, smoke was responsible for ∼10%–30% reduction in solar power production during peak hours as recorded by the California Independent System Operator (CAISO) sites. In this study, we focus on a 9 d period in September when wildfire smoke had a profound impact on solar energy production. During the smoke episodes, hour-ahead forecasts utilized by CAISO did not include the effects of smoke and therefore overestimated the expected power production by ∼10%–50%. Here we use multiple observational networks and a numerical weather prediction (NWP) model to show that the wildfire events of 2020 had a significantly detrimental influence on solar energy production due to high aerosol loading. We find that including the contribution of biomass burning particles greatly improves the day-ahead solar energy bias forecast of both global horizontal irradiance and direct normal irradiance by nearly ∼50%. Our results suggest that a more comprehensive treatment of aerosols, including biomass burning aerosols, in NWP models may be an important consideration for energy grid balancing, in addition to solar resource assessment, as solar power reliance increases.

On 30 December 2021, the Marshall Fire devastated the Boulder, Colorado region. The fire initiated in fine fuels in open space just southeast of Boulder and spread rapidly due to the strong, downslope winds that penetrated into the Boulder Foothills. Despite the increasing occurrence of wildland‐urban interface (WUI) disasters, many questions remain about how fires progress through vegetation and the built environment. To help answer these questions for the Marshall Fire, we use a coupled fire‐atmosphere model and Doppler on Wheels (DOW) observations to study the fire's progression as well as examine the physical drivers of its spread. Evaluation of the model using the DOW suggests that the model is able to capture general characteristics of the flow field; however, it does not produce as robust of a hydraulic jump as the one observed. Our results highlight limitations of the model that should be addressed for successful WUI simulations.