Search

Abstract Mesoscale‐to‐microscale coupling is an important tool for conducting turbulence‐resolving multiscale simulations of realistic atmospheric flows, which are crucial for applications ranging from wind energy to wildfire spread studies. Different techniques are used to facilitate the development of realistic turbulence in the large‐eddy simulation (LES) domain while minimizing computational cost. Here, we explore the impact of a simple and computationally efficient Stochastic Cell Perturbation method using momentum perturbation (SCPM‐M) to accelerate turbulence generation in boundary‐coupled LES simulations using the Weather Research and Forecasting model. We simulate a convective boundary layer (CBL) to characterize the production and dissipation of turbulent kinetic energy (TKE) and the variation of TKE budget terms. Furthermore, we evaluate the impact of applying momentum perturbations of three magnitudes below, up to, and above the CBL on the TKE budget terms. Momentum perturbations greatly reduce the fetch associated with turbulence generation. When applied to half the vertical extent of the boundary layer, momentum perturbations produce an adequate amount of turbulence. However, when applied above the CBL, additional structures are generated at the top of the CBL, near the inversion layer. The magnitudes of the TKE budgets produced by SCPM‐M when applied at varying heights and with different perturbation amplitudes are always higher near the surface and inversion layer than those produced by No‐SCPM, as are their contributions to the TKE. This study provides a better understanding of how SCPM‐M reduces computational costs and how different budget terms contribute to TKE in a boundary‐coupled LES simulation.

The intensity and frequency of wildfires in California (CA) have increased in recent years, causing significant damage to human health and property. In October 2007, a number of small fire events, collectively referred to as the Witch Creek Fire or Witch Fire started in Southern CA and intensified under strong Santa Ana winds. As a test of current mesoscale modeling capabilities, we use the Weather Research and Forecasting (WRF) model to simulate the 2007 wildfire event in terms of meteorological conditions. The main objectives of the present study are to investigate the impact of horizontal grid resolution and planetary boundary layer (PBL) scheme on the model simulation of meteorological conditions associated with a Mega fire. We evaluate the predictive capability of the WRF model to simulate key meteorological and fire-weather forecast parameters such as wind, moisture, and temperature. Results of this study suggest that more accurate predictions of temperature and wind speed relevant for better prediction of wildfire spread can be achieved by downscaling regional numerical weather prediction products to 1 km resolution. Furthermore, accurate prediction of near-surface conditions depends on the choice of the planetary boundary layer parameterization. The MYNN parameterization yields more accurate prediction as compared to the YSU parameterization. WRF simulations at 1 km resolution result in better predictions of temperature and wind speed than relative humidity during the 2007 Witch Fire. In summary, the MYNN PBL parameterization scheme with finer grid resolution simulations improves the prediction of near-surface meteorological conditions during a wildfire event.

Abstract 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.

Abstract Improving the calculation of land‐atmosphere fluxes of heat and water vapor in mountain terrain requires better resolution of thermally driven diurnal winds (i.e., valley, slope winds) due to differential heating by terrain and radiative fluxes. In this study, the Weather Research and Forecasting model is used to simulate flow in large‐eddy simulation (LES) mode over the complex terrain of the Fortress Mountain and Marmot Creek research basins, Kananaskis Valley, Canadian Rockies, Alberta in mid‐summer. The model was used to examine the temporal and spatial evolution of local winds and near‐surface boundary layer processes with variability in topography and elevation. Numerically resolving complex terrain wind flow effects require smaller grid cell size. However, the use of terrain‐following coordinates in most numerical weather prediction models results in large numerical errors when flow over steep terrain is simulated. These errors propagate through the domain and can result in numerical instability. To avoid this issue when simulating flow over steep terrain a local smoothing approach was used, where smoothing is applied only where slope exceeds some predetermined threshold. LES results from local smoothing were compared with a mesoscale model and LES with global smoothing. Simulations are evaluated using sounding data and meteorological stations. The differences in flow patterns and reversals in two mountain basins suggest that valley geometry and volume is relevant to the break up of inversion layers, removal of cold‐air pools, and strength of thermally driven winds.

Abstract This study evaluates the efficacy of current satellite observing systems to detect methane point sources from typical oil and gas production (O&G) facilities using a novel very high-resolution methane concentration dataset generated using a microscale model. Transport and dispersion of typical methane emissions from seven well pads were simulated and the column enhancements for pseudo satellite pixel sizes of 3, 1, and 0.05 km were examined every second of the 2-h simulations (7200 realizations). The detectability of plumes increased with a pixel resolution, but two orders of magnitude change in emission rates at the surface results only in about 0.4%, 1.6%, and 47.8% enhancement in the pseudo-satellite retrieved methane column at 3, 1, and 0.05 km, respectively. Average methane emission rates estimated by employing the integrated mass enhancement (IME) method to column enhancements at 0.05 km showed an underestimation of the mean emissions by 0.2–6.4%. We show that IME derived satellite-based inversions of methane emissions work well for large persistent emission sources (e.g., super emitters), however, the method is ill-suited to resolve short-term emission fluctuations (

Abstract Smoke from wildfires or burning biomass directly affects air quality and weather through modulating cloud microphysics and radiation. A simple wildfire emission coupling of black carbon (BC) and organic carbon (OC) with microphysics was implemented using the Weather Research and Forecasting model's fire module. A set of large‐eddy simulations inspired by unique surface and upper atmospheric observations from the 2021 Santa Coloma de Queralt Fire (Spain) were conducted to investigate the influence of background conditions and interactions between atmospheric and fire processes such as fire smoke, ambient moisture, and latent heat release on the formation and evolution of pyroconvective clouds. While the microphysical impact of BC and OC emissions on the dynamics of fire behavior is minimal on short time scales (

In this study, we focus on the effects of fuel bed representation and fire heat and smoke distribution in a coupled fire-atmosphere simulation platform for two landscape-scale fires: the 2018 Camp Fire and the 2021 Caldor Fire. The fuel bed representation in the coupled fire-atmosphere simulation platform WRF-Fire currently includes only surface fuels. Thus, we enhance the model by adding canopy fuel characteristics and heat release, for which a method to calculate the heat generated from canopy fuel consumption is developed and implemented in WRF-Fire. Furthermore, the current WRF-Fire heat and smoke distribution in the atmosphere is replaced with a heat-conserving Truncated Gaussian (TG) function and its effects are evaluated. The simulated fire perimeters of case studies are validated against semi-continuous, high-resolution fire perimeters derived from NEXRAD radar observations. Furthermore, simulated plumes of the two fire cases are compared to NEXRAD radar reflectivity observations, followed by buoyancy analysis using simulated temperature and vertical velocity fields. The results show that while the improved fuel bed and the TG heat release scheme have small effects on the simulated fire perimeters of the wind-driven Camp Fire, they affect the propagation direction of the plume-driven Caldor Fire, leading to better-matching fire perimeters with the observations. However, the improved fuel bed representation, together with the TG heat smoke release scheme, leads to a more realistic plume structure in comparison to the observations in both fires. The buoyancy analysis also depicts more realistic fire-induced temperature anomalies and atmospheric circulation when the fuel bed is improved.

Accurate estimation of fuels is essential for wildland fire simulations as well as decision-making related to land management. Numerous research efforts have leveraged remote sensing and machine learning for classifying land cover and mapping forest vegetation species. In most cases that focused on surface fuel mapping, the spatial scale of interest was smaller than a few hundred square kilometers; thus, many small-scale site-specific models had to be created to cover the landscape at the national scale. The present work aims to develop a large-scale surface fuel identification model using a custom deep learning framework that can ingest multimodal data. Specifically, we use deep learning to extract information from multispectral signatures, high-resolution imagery, and biophysical climate and terrain data in a way that facilitates their end-to-end training on labeled data. A multi-layer neural network is used with spectral and biophysical data, and a convolutional neural network backbone is used to extract the visual features from high-resolution imagery. A Monte Carlo dropout mechanism was also devised to create a stochastic ensemble of models that can capture classification uncertainties while boosting the prediction performance. To train the system as a proof-of-concept, fuel pseudo-labels were created by a random geospatial sampling of existing fuel maps across California. Application results on independent test sets showed promising fuel identification performance with an overall accuracy ranging from 55% to 75%, depending on the level of granularity of the included fuel types. As expected, including the rare—and possibly less consequential—fuel types reduced the accuracy. On the other hand, the addition of high-resolution imagery improved classification performance at all levels.

The National Center for Atmospheric Research (NCAR) recently updated the comprehensive wind power forecasting system in collaboration with Xcel Energy addressing users’ needs and requirements by enhancing and expanding integration between numerical weather prediction and machine-learning methods. While the original system was designed with the primary focus on day-ahead power prediction in support of power trading, the enhanced system provides short-term forecasting for unit commitment and economic dispatch, uncertainty quantification in wind speed prediction with probabilistic forecasting, and prediction of extreme events such as icing. Furthermore, the empirical power conversion machine-learning algorithms now use a quantile approach to data quality control that has improved the accuracy of the methods. Forecast uncertainty is quantified using an analog ensemble approach. Two methods of providing short-range ramp forecasts are blended: the variational doppler radar analysis system and an observation-based expert system. Extreme events, specifically changes in wind power due to high winds and icing, are now forecasted by combining numerical weather prediction and a fuzzy logic artificial intelligence system. These systems and their recent advances are described and assessed.

Monitoring the fuel moisture content (FMC) of 10 h dead vegetation is crucial for managing and mitigating the impact of wildland fires. The combination of in situ FMC observations, numerical weather prediction (NWP) models, and satellite retrievals has facilitated the development of machine learning (ML) models to estimate 10 h dead FMC retrievals over the contiguous US (CONUS). In this study, ML models were trained using variables from the National Water Model, the High-Resolution Rapid Refresh (HRRR) NWP model, and static surface properties, along with surface reflectances and land surface temperature (LST) retrievals from the Visible Infrared Imaging Radiometer Suite (VIIRS) instrument on the Suomi-NPP satellite system. Extensive hyper-parameter optimization resulted in skillful FMC models compared to a daily climatography RMSE (+44%) and an hourly climatography RMSE (+24%). Notably, VIIRS retrievals played a significant role as predictors for estimating 10 h dead FMC, demonstrating their importance as a group due to their high band correlation. Conversely, individual predictors within the HRRR group exhibited relatively high importance according to explainability techniques. Removing both HRRR and VIIRS retrievals as model inputs led to a significant decline in performance, particularly with worse RMSE values when excluding VIIRS retrievals. The importance of the VIIRS predictor group reinforces the dynamic relationship between 10 h dead fuel, the atmosphere, and soil moisture. These findings underscore the significance of selecting appropriate data sources when utilizing ML models for FMC prediction. VIIRS retrievals, in combination with selected HRRR variables, emerge as critical components in achieving skillful FMC estimates.

Abstract Over the last few years, the concept of incorporating aerial vehicles into the urban environment for diverse purposes has attracted ample interest and investment. These purposes cover a broad spectrum of applications, from larger vehicles designed for passenger transport, to package delivery and inspection/surveillance missions performed by small unmanned drones. While these Advanced Air Mobility (AAM) operations have the potential to alleviate bottlenecks arising from saturated surface transportation networks, there are a number of challenges that need to be addressed to make these operations safe and viable. One challenge is predicting weather effects within the urban environment with the required level of spatiotemporal fidelity, which current operational weather models fail to provide due to the use of coarse grid spacings (a few kilometers) constrained by the predictive performance limitations of traditional computer architectures. Herein, we demonstrate how FastEddy®, a microscale model that exploits the accelerated nature of graphics processing units for high‐performance computing, can be used to understand and predict urban weather impacts from seasonal, day‐to‐day, diurnal, and sub‐hourly scales. To that end, we efficiently perform more than 50 telescoped simulations of microscale urban effects at street‐scale (5 m grid spacing) driven by realistic weather over a 20 km2 region centered at the downtown area of Dallas, Texas. Our analyses demonstrate that urban‐weather interactions at the street‐scale are complex and tightly connected, which is of utmost relevance to AAM operations. These demonstrations reveal the capability of such models to provide real‐time weather hazard avoidance products tailored to capture microscale urban effects.

Wildland fire behavior models are often initiated using the detection information listed in incident reports. This information carries an unknown amount of uncertainty, though it is often the most readily available ignition data. To determine the extent to which the use of detection information affects wildland fire forecasts, this research examines the range of burned area values and propagation directions resulting from different initiation point locations and times. We examined the forecasts for ten Colorado 2018 wildland fire case studies, each initiated from a set of 17 different point locations, and three different starting times (a total of 520 case study simulations). The results show that the range of forecast burned area and propagation direction values is strongly affected by the location of the initiation location, and to a lesser degree by the time of initiation.

Disturbance events can happen at a temporal scale much faster than wildland fire fuel data updates. When used as input for wildland fire behavior models, outdated fuel datasets can contribute to misleading forecasts, which have implications for operational firefighting, mitigation, and wildland fire research. Remote sensing and machine learning methods can provide a solution for on-demand fuel estimation. Here, we show a proof of concept using C-band synthetic aperture radar and multispectral imagery, land cover classes, and tree mortality surveys to train a random forest classifier to estimate wildland fire fuel data in the East Troublesome Fire (Colorado) domain. The algorithm classified over 80% of the test dataset correctly, and the resulting wildland fire fuel data was used to simulate the East Troublesome Fire using the coupled atmosphere—wildland fire behavior model, WRF-Fire. The simulation using the modified fuel inputs, where 43% of original fuels are replaced with fuels representing dead trees, improved the burn area forecast by 38%. This study demonstrates the need for up-to-date fuel maps available in real time to provide accurate prediction of wildland fire spread, and outlines the methodology based on high-resolution satellite observations and machine learning that can accomplish this task.

Abstract As a first step toward achieving full physics urban weather simulation capabilities within the resident‐GPU large‐eddy simulation (LES) FastEddy® model, we have implemented and verified/validated a method for explicit representation of building effects. Herein, we extend the immersed body force method (IBFM) from Chan and Leach (2007, https://doi.org/10.1175/2006JAMC1321.1) to (i) be scale independent and (ii) control building surface temperatures. Through a specific drag‐like term in the momentum equations, the IBFM is able to enforce essentially zero velocities within the buildings, in turn resulting in a no‐slip boundary condition at the building walls. In addition, we propose similar forcing terms in the energy and mass conservation equations that allow an accurate prescription of the building temperature. The extended IBFM is computationally efficient and has the potential to be coupled to building energy models. The IBFM exhibits excellent agreement with laboratory experiments of an array of staggered cubes at a grid spacing of Δ=1 mm, demonstrating the applicability of the method beyond the atmospheric scale. In addition, the IBFM is validated at atmospheric scale through simulations of downtown Oklahoma City ( Δ=2 m) using data collected during the Joint Urban 2003 (JU03) field campaign. Our LES IBFM results for mean wind speed, turbulence kinetic energy, and SF6 transport and dispersion compare well to observations and produce turbulence spectra that are in good agreement with sonic anemometer data. Statistical performance metrics for the JU03 simulations are within the range of other LES models in the literature.

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 ( $\gt$ 40 000 km ^2 ) 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.

A real case long-term nested large eddy simulation (LES) of 25-day duration is performed using the WRF-LES modelling system, with a maximum horizontal grid resolution of 111 m, in order to explore the ability of the model to reproduce the turbulence magnitudes within the first tens of metres of the boundary layer. Sonic anemometer measurements from a 60-m tower installed during the Boundary Layer Late Afternoon and Sunset Turbulence (BLLAST) field campaign are used for verification, which is focused on the turbulent magnitudes in order to assess the success and limitations in resolving turbulent flow characteristics. The mesoscale and LES simulations reproduce the wind speed and direction fairly well, but only LES is able to reproduce the energy of eddies with lifetimes shorter than a few hours. The turbulent kinetic energy in LES simulation is generally underestimated during the daytime, mainly due to a vertical velocity standard deviation that is too low. The turbulent heat flux is misrepresented in the model, probably due to the inaccuracy of the sub-grid scheme.

A modern renewable energy forecasting system blends physical models with artificial intelligence to aid in system operation and grid integration. This paper describes such a system being developed for the Shagaya Renewable Energy Park, which is being developed by the State of Kuwait. The park contains wind turbines, photovoltaic panels, and concentrated solar renewable energy technologies with storage capabilities. The fully operational Kuwait Renewable Energy Prediction System (KREPS) employs artificial intelligence (AI) in multiple portions of the forecasting structure and processes, both for short-range forecasting (i.e., the next six hours) as well as for forecasts several days out. These AI methods work synergistically with the dynamical/physical models employed. This paper briefly describes the methodology used for each of the AI methods, how they are blended, and provides a preliminary assessment of their relative value to the prediction system. Each operational AI component adds value to the system. KREPS is an example of a fully integrated state-of-the-science forecasting system for renewable energy.