An official website of the United States government
Abstract 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.
more »
« less
This content will become publicly available on August 16, 2026
Evaluating WRF Multiscale Wind Simulations in Complex Terrain: Insights From the Perdigão Field Campaign
Abstract Accurate microscale flow simulations are essential for assessing wind characteristics in complex terrain. This study evaluates a large ensemble of multiscale simulations, including large‐eddy simulations (LES), using the Weather Research and Forecasting model (WRF) over Perdigão, Portugal, driven by boundary conditions from multiple global data sets. Simulations are compared with data from the Perdigão field campaign, including radiosonde and flux tower measurements. Results show that LES, using high‐resolution topography and land use data, better replicate flow features and dynamics, providing valuable insights for wind resource quantification and mapping. We identify that model performance varies spatially. The RMSE of wind speed at 10 m at ridge towers is 5.65 m , while at valley towers, it is lower (2.28 m ), and variation across runs is greater for higher wind speeds. Temporally, surface winds show substantial variability throughout the day, posing greater modeling challenges during nighttime and synoptic transitions. This variability is not observed at 100 m, where topographic effects are less dominant and RMSE remains consistent across runs. Simulations driven by hourly boundary conditions perform best. However, drawing general conclusions about optimal turbulence modeling in the gray zone remains challenging due to microscale meteorology in complex terrain. Wind field characteristics are sensitive to turbulence scheme choice, particularly in the boundary layer, while above it, wind behavior is mainly influenced by boundary conditions. These results help identify key factors driving model variability and biases, which may guide future model developments to enhance wind flow simulation accuracy and reliability.
more »
« less
- Award ID(s):
- 2236504
- PAR ID:
- 10653225
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Atmospheres
- Volume:
- 130
- Issue:
- 15
- ISSN:
- 2169-897X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present a new ensemble of 36 numerical experiments aimed at comprehensively gauging the sensitivity of nested large-eddy simulations (LES) driven by large-scale dynamics. Specifically, we explore 36 multiscale configurations of the Weather Research and Forecasting (WRF) Model to simulate the boundary layer flow over the complex topography at the Perdigão field site, with five nested domains discretized at horizontal resolutions ranging from 11.25 km to 30 m. Each ensemble member has a unique combination of the following input factors: (i) large-scale initial and boundary conditions, (ii) subgrid turbulence modeling in thegray zoneof turbulence, (iii) subgrid-scale (SGS) models in LES, and (iv) topography and land-cover datasets. We probe their relative importance for LES calculations of velocity, temperature, and moisture fields. Variance decomposition analysis unravels large sensitivities to topography and land-use datasets and very weak sensitivity to the LES SGS model. Discrepancies within ensemble members can be as large as 2.5 m s−1for the time-averaged near-surface wind speed on the ridge and as large as 10 m s−1without time averaging. At specific time points, a large fraction of this sensitivity can be explained by the different turbulence models in the gray zone domains. We implement a horizontal momentum and moisture budget routine in WRF to further elucidate the mechanisms behind the observed sensitivity, paving the way for an increased understanding of the tangible effects of the gray zone of turbulence problem. Significance StatementSeveral science and engineering applications, including wind turbine siting and operations, weather prediction, and downscaling of climate projections, call for high-resolution numerical simulations of the lowest part of the atmosphere. Recent studies have highlighted that such high-resolution simulations, coupled with large-scale models, are challenging and require several important assumptions. With a new set of numerical experiments, we evaluate and compare the significance of different assumptions and outstanding challenges in multiscale modeling (i.e., coupling large-scale models and high-resolution atmospheric simulations). The ultimate goal of this analysis is to put each individual assumption into the wider perspective of a realistic problem and quantify its relative importance compared to other important modeling choices.more » « less
-
The large-scale deployment of wind energy underscores the critical need for accurate resource characterization to reduce uncertainty in power estimates and to enable the installation of wind farms in increasingly complex terrains. Accurate wind resource assessment in peri-urban and moderately complex terrains remains a significant challenge due to spatial heterogeneity in surface terrain features and atmospheric thermal stability. This study investigates the influence of surface complexity and atmospheric stratification on vertical wind profiles at a utility-scale wind turbine site in Cedar Rapids, Iowa. One year of multi-level wind data from a 106-meter-tall meteorological tower were analyzed to quantify variations in the wind shear exponent α, wind direction veer, and horizontal turbulence intensity (TI) across open-field and complex-surface wind sectors and four thermal stability classes, defined by the bulk Richardson number Rib. The results show that the wind shear exponent α increases systematically with atmospheric stability. Over the open-field terrain, α ranges from 0.11 in unstable conditions to 0.45 in strongly stable conditions, compared to 0.17 and 0.40 over the complex surface. A pronounced diurnal variation in α was observed, particularly during the summer months. Wind veer was greatest and exceeded 30° under strongly stable conditions over open terrain. Elevated TI values peaked at 32 m in height due to flow separation and wake turbulence from nearby vegetation and sloping terrain. These findings highlight the importance of incorporating terrain-induced and thermally driven variability into wind resource assessments to improve power prediction and turbine siting in complex heterogeneous terrain environments.more » « less
-
Abstract Understanding the influence of roughness and terrain slope on stable boundary layer turbulence is challenging. This is investigated using observations collected from October to November of 2018 during the Stable Atmospheric Variability ANd Transport (SAVANT) field campaign conducted in a shallow sloping Midwestern field. We analyze the turbulence velocity scale and its variation with the mean wind speed using observations up to 10–20 m on four meteorological towers located along a shallow gully. The roughness length for momentum over this complex terrain varied with wind direction from 0.0049 m to a maximum of 0.12 m for winds coming through deciduous trees present in the field. The variation of the turbulence velocity with wind speed shows a transition from a weak wind regime to a stronger wind regime, as reported by past studies. This transition is not observed for winds coming from the tree area, where turbulence is enhanced even for weak wind speeds. For weak stratification and stronger winds, the turbulent velocity scale increased with an increase in roughness while the terrain slope is seen to have a weak influence. The sizes of the dominant turbulent eddies seen from the vertical velocity power spectra are observed to be larger for winds coming through the tree area. The turbulence enhancement by the trees is found to be strong within a fetch distance of 7 times the tree height and not observable at 16 times of the tree height.more » « less
-
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.more » « less
