Search NASA⌕ Search

SEARCH · Search NASA

Results for “atmospheric boundary layer”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Assessing the sensitivity of land-atmosphere coupling strength to boundary and surface layer parameters in the WRF model over Amazon

Modeling tools can be used to diagnose regional land-atmosphere (L-A) coupling strength in the absence of sufficient observations, but subject to uncertainties associated with parameters in model physical parameterizations. Different sensitivity analysis (SA) approaches may lead to different conclusions about the underlying sensitivities. In this study, we quantify simulation uncertainties related to parameter perturbations, and use different approaches to conduct parameter SA on the WRF model pertaining to L-A coupling strength for simulations over the Amazon region. A total of twenty parameters from the Yonsei University (YSU) planetary boundary layer (PBL) and the revised MM5 surface layer (SL) schemes were selected in this analysis. Three different SA methods, the Morris One-at-A-Time (MOAT) method, the Multivariate Adaptive Regression Splines (MARS) method, and the Sobol’ method, were employed to analyze seven WRF-simulated variables and five L-A coupling metrics. Results show that 1) parameter perturbations cause large simulation uncertainties which are comparable to those in the observations; 2) three different SA methods give consistent L-A coupling strength outcomes; 3) six out of the twenty parameters contribute 80%–95% of the total variance in the metrics analyzed, and first-order effects dominate over interaction effects; 4) the twelve variables/metrics of interest show similar sensitivity patterns to the selected parameters, which is consistent across all the methods used. Physical mechanisms for how the sensitive parameters act in determining the L-A coupling strength and associated variables also are illustrated. Our results will help quantifying L-A coupling strength and establishing a basis for parameter calibration over the Amazon region.

Wang, Chen↗
DOE OSTI2022

A Vertically Resolved MSE Framework Highlights the Role of the Boundary Layer in Convective Self-Aggregation

Abstract Convective self-aggregation refers to a phenomenon in which random convection can self-organize into large-scale clusters over an ocean surface with uniform temperature in cloud-resolving models. Previous literature studies convective aggregation primarily by analyzing vertically integrated (VI) moist static energy (MSE) variance. That is the global MSE variance, including both the local MSE variance at a given altitude and the covariance of MSE anomalies between different altitudes. Here we present a vertically resolved (VR) MSE framework that focuses on the local MSE variance to study convective self-aggregation. Using a cloud-resolving simulation, we show that the development of self-aggregation is associated with an increase of local MSE variance, and that the diabatic and adiabatic generation of the MSE variance is mainly dominated by the boundary layer (BL; the lowest 2 km). The results agree with recent numerical simulation results and the available potential energy analyses showing that the BL plays a key role in the development of self-aggregation. Additionally, we find that the lower free troposphere (2–4 km) also generates significant MSE variance in the first 15 days. We further present a detailed comparison between the global and local MSE variance frameworks in their mathematical formulation and diagnostic results, highlighting their differences.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Doppler-Lidar Evaluation of HRRR-Model Skill at Simulating Summertime Wind Regimes in the Columbia River Basin during WFIP2

Complex-terrain locations often have repeatable near-surface wind patterns, such as synoptic gap flows and local thermally forced flows. An example is the Columbia River Valley in east-central Oregon-Washington, a significant wind-energy-generation region and the site of the Second Wind-Forecast Improvement Project (WFIP2). Data from three Doppler lidars deployed during WFIP2 define and characterize summertime wind regimes and their large-scale contexts, and provide insight into NWP model errors by examining differences in the ability of a model [NOAA’s High-Resolution Rapid-Refresh (HRRR-version1)] to forecast wind-speed profiles for different regimes. Seven regimes were identified based on daily time series of the lidar-measured rotor-layer winds, which then suggested two broad categories. First, in three regimes the primary dynamic forcing was the large-scale pressure gradient. Second, in two regimes the dominant forcing was the diurnal heating-cooling cycle (regional sea-breeze-type dynamics), including the marine intrusion previously described, which generates strong nocturnal winds over the region. The other two included a hybrid regime and a non-conforming regime. For the large-scale pressure-gradient regimes, HRRR had wind-speed biases of ~1 m s -1 and RMSEs of 2-3 m s -1 . Errors were much larger for the thermally forced regimes, owing to the premature demise of the strong nocturnal flow in HRRR. Thus, the more dominant the role of surface heating in generating the flow, the larger the errors. Major errors could result from surface heating of the atmosphere, boundary-layer responses to that heating, and associated terrain interactions. Finally, measurement/modeling research programs should be aimed at determining which modeled processes produce the largest errors, so those processes can be improved and errors reduced.

17 WIND ENERGY↗
DOE OSTI2024

Characterizing model uncertainties in simulated coast-to-offshore wind over the northeast U.S. using multi-platform measurements from the TCAP field campaign

Numerical weather prediction (NWP) models, such as the Weather Research and Forecasting (WRF) model, are widely used to provide estimates of the offshore wind energy resource owing to their large spatial coverage compared to available observations. Nevertheless, spatiotemporal distribution of model biases is highly dependent on factors including model configuration, location, and the interplay of multi-scale physical processes. Here, in this study, we focus on the characterization of model uncertainties in simulated coast-to-offshore winds over the northeast U.S., by varying sea surface temperature (SST) forcings, surface layer (SL) and planetary boundary layer (PBL) parameterizations, as well as identifying biases that may be directly passed from initial and boundary conditions. Multiple measurements, including aircraft data collected during the U.S. Department of Energy's Two-Column Aerosol Project (TCAP) experiment, are used to constrain the model results and facilitate quantitative comparisons. Our analysis indicates while SST forcing has notable impacts on simulated air temperature and moisture within PBL, the modeled winds are in general more sensitive to the choices of SL and PBL physics than to SST. The model’s forcing data not only controls the vertical dependence of wind speed errors, but also alters regional variability in wind speed’s spatial correlation. Bias comparisons between ERA5 reanalysis and ensemble simulations revealed significant similarity, particularly in wind speed biases during winter, underscoring their dependency on initial and boundary conditions. Coastal and offshore near-surface wind speed biases tend to exhibit much higher similarity in winter than in summer due to the presence of much stronger and more persistent synoptic wind conditions. This study highlights the importance of accurate atmospheric forcing and parameterization choices in improving wind forecasts and suggests the potential for extrapolating coastal wind biases to offshore locations, aiding wind energy forecasting and informing the Wind Forecast Improvement Project-3 (WFIP3).

17 WIND ENERGY↗
DOE OSTI2023

Using observational mean-flow data to drive large-eddy simulations of a diurnal cycle at the SWiFT site

Reproducing realistic date- and site-specific unsteady wind conditions in large-eddy simulations is becoming increasingly useful in wind energy. How to run a large-eddy simulation to match observed conditions, however, remains an open research question. One approach that has received considerable attention is mesoscale-to-microscale coupling, in which information about the mesoscale weather, most commonly acquired from a mesoscale numerical weather model, is passed on to a microscale model. In this paper, we demonstrate how the recently developed profile-assimilation technique, a form of mesoscale-to-microscale coupling, can be used to drive large-eddy simulations solely based on observed mean-flow profiles at a single location, bypassing the need for auxiliary mesoscale simulations. The new approach is evaluated for a diurnal cycle at the Scaled Wind Farm Technology site. Observed mean-flow profiles from the ground up to a height of 2 km are reconstructed by aggregating measurements from multiple instruments, and gaps in the data are infilled with natural neighbor interpolation. We perform nine simulations using various forcing approaches to deal with data limitations. The results show that it is indeed possible to drive microscale large-eddy simulation with observations using the profile-assimilation technique, notwithstanding large gaps in virtual potential temperature measurements. However, profile assimilation with vertical smoothing of the error between the desired and actual profiles is required. Without that smoothing, the microscale simulations develop unrealistically high turbulence levels under many situations. Finally, we show that simulated mesoscale data can account for missing observations, although care is needed as both data sources are not necessarily compatible.

17 WIND ENERGY↗
DOE OSTI2025

Exploring the impact of surface topography on Rayleigh-Bénard dry convection in the Pi cloud chamber using OpenFOAM: In cylindrical and rectangular geometries

The Pi convection-cloud chamber can generate steady-state turbulence in both rectangular and cylindrical shapes via Rayleigh-Bénard convection (RBC) by maintaining warm bottom and cold top surfaces. Although most experiments in the Pi chamber were conducted in cylindrical shapes, all previous Pi chamber simulations were conducted in a rectangular shape due to the limitations of those models to discretize a cylindrical domain when using the finite difference method therein. Here, we use OpenFOAM, an open-source finite-volume-based Computational Fluid Dynamics (CFD) software package, to conduct Large-Eddy Simulation (LES) of dry RBC in the Pi chamber at high Rayleigh numbers (10 8 to 10 9 ). Results show that large-scale circulation (LSC) direction varies in the chamber with a constant side wall temperature. Imposing a slight temperature imbalance at the side wall ranging from 0.1 to 0.7 degrees can lock the LSC, aligning better with Pi chamber observations, particularly at higher Rayleigh numbers. In addition, we examine the impact of surface topography on LSC and heat transfer in RBC systems within cylindrical and rectangular shapes under varying conditions. Results show that roughing top/bottom surfaces by adding bars of a few tens millimeters height can strengthen thermal plumes and enhance temperature fluctuations in the chamber. Furthermore, we observe that different bar height configurations lead to notable changes in LSC orientation and thermal stratification, highlighting the complex interactions between surface features and convection patterns. This finding highlights how surface topography and chamber geometry affect Rayleigh-Bénard convection, improving understanding of turbulent heat transfer and atmospheric boundary-layer processes. Direct Numerical Simulations (DNS) are also conducted to validate LES results. In conclusion, while LES effectively captures qualitative behaviors seen in DNS, it tends to underestimate velocity variances near walls, illustrating a trade-off between computational efficiency and accuracy.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Effect of the integral length scales of turbulent inflows on wind turbine loads

As wind turbines become larger, the fluctuations in the inflow become increasingly influential in the turbine structural loading. These fluctuations are characterized by the integral length scale, which represents the average size of the largest energy-containing turbulent eddies. Current design standards neglect the varying integral length scales that characterize inflows of wind turbines in operation. Using large-eddy simulations, we generate turbulent inflows of varying integral length scales and quantify the loads of the IEA 15-MW reference wind turbine. Results illustrate that the impact of turbulence on rotor and tower loads is up to 10 times greater than the impact of the mean shear profile. Increasing integral length scales from 0.3x to 0.5x the rotor diameter reduces blade root flapwise moments and rotor and tower loads. Increasing integral length scales from 0.5x to 1.7x the rotor diameter increases the rotor aerodynamic thrust force and the blade root flapwise shear loads and decreases the tower base torsional moment and the tilting and yawing rotor aerodynamic moments. Additionally, turbulence intensity has a greater impact on wind turbine loads than integral length scales. In conclusion, findings indicate that design standards should consider varying integral length scales for accurate wind turbine loading characterization in turbulent inflow conditions.

17 WIND ENERGY↗
DOE OSTI2025

Unsteady Land-Sea Breeze Circulations in the Presence of a Synoptic Pressure Forcing

Unsteady land-sea breezes (LSBs) that result from time-varying surface temperature contrasts Δθ(t) are explored in the presence of a constant synoptic pressure forcing, M g , oriented from sea to land (α = 0°) or land to sea (α = 180°). Large eddy simulations reveal the development of four distinctive regimes, depending on the joint interaction between M g , α, and Δθ(t) in modulating the fine-scale dynamics. Time lags, computed as the shifts that maximize correlation coefficients of the velocity between the unsteady and the corresponding steady scenarios at Δθ = Δθ max , are found to be significant and to extend 2 hr longer for α = 0° compared to α = 180°. These diurnal dynamics result in nonequilibrium conditions that are significantly affected by the flow history, and that behave differently over the two patches for the different α’s. Turbulence is found to be out of equilibrium with the mean flow, and the mean itself is found to be out of equilibrium with the thermal forcing. The sea surface heat flux is consistently more sensitive than its land counterpart to the time-varying external forcing Δθ(t), and more so for synoptic forcing from land to sea (α = 180°). Hence, although the land reaches equilibrium faster, the sea patch is found to exert a stronger control on the turbulence-mean flow equilibrium response. Finally, the vertical velocity profile at the shore and shore-normal velocity transects at the first grid level are shown to encode the multiscale regimes of the LSBs evolution and can thus be used to identify these regimes using k-means clustering.

58 GEOSCIENCES↗
DOE OSTI2026

Pinned Clouds over Industrial Sources of Heat during TRACER

An analysis of stereo-camera data from the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign in Houston, Texas, reveals the existence of pinned clouds sitting atop industrial heat sources. These are not plumes of vapor and condensate emanating from stacks: The pinned shallow cumuli have cloud bases at about the lifting condensation of near-surface air. On many mornings, the pinned clouds are the only shallow cumuli in the field of view of the stereo cameras, persisting by themselves for at least an hour. On those mornings, the lower atmosphere straddles the boundary between stable and unstable, and the waste heat from industrial facilities is able to pin moist convection overhead. When solar heating of the surface becomes sufficient later in the morning, those mornings transition to having widespread shallow moist convection. Occurring within the field of view of the stereo cameras and other TRACER instruments, these pinned clouds represent steady plume-like moist convection triggered by known heat sources in a well-characterized atmosphere, making them a rich target for further study. Significance Statement This article introduces steady-state cumulus clouds observed during the Tracking Aerosol Convection Interactions Experiment (TRACER) campaign that were pinned to industrial heat sources. These pinned clouds were present in the early morning only when conditions were primed for moist convection, but for which there was not yet enough solar heating to trigger widespread convection. They persist in a statistical steady state and therefore provide a natural setting to explore the mechanisms controlling convective dynamics, which play a central role in Earth’s atmospheric circulation and radiative balance.

Atmosphere↗
DOE OSTI2020

Assimilating Scanning Radar Data into High-Resolution Models

This report documents the findings of a project aiming to improve the initial conditions of km-scale simulations of deep convective storms by assimilating cloud-scale weather radar observations. More accurate numerical analyses will increase the effectiveness of research efforts using LES cloud models as tool to better understand land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development. Work reported on herein includes: i) assessment of radar data sets for use in data assimilation (‘DA’) experiments, ii) quality control of the radar data set and format conversion to one acceptable to DA schemes utilized by the Weather Research and Forecasting (WRF) model, and iii) examination of the sensitivity of the numerical representation of cloud-scale wind and microphysical features to a variety of tunable DA parameters.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Assimilating Scanning Radar Data into Cloud-Scale Models

The aim of this project is to improve the initial conditions of large eddy simulations (LES) of deep convective storms by assimilating scanning radar reflectivity and radial velocity data as well as local environmental observations. Accomplishing this goal often results in more accurate numerical forecasts of clouds and precipitation than if the background initial meteorological conditions are solely relied upon for convection initiation and intensification, as is typical of most operational and research models. Conducted in a storm-scale LES framework, these more accurate numerical analyses can serve as tool to better understand and parameterize land-atmosphere coupling, boundary layer turbulence, and cloud processes, each used for model parameterization development of convection. Our work examined the sensitivity of the numerical representation of cloud-scale wind and microphysical structures to: i) the assimilation of radiosonde observations characterizing the near-storm inflow, ii) the assimilation of radar data of varying resolution for model simulations with grid spacings of 1-km and 0.25-km, and iii) to the assimilation of three-dimensional plan-projection-indicator (PPI) and range-height indicator (RHI) radar volumes. Few studies have explored deep convective storm analyses and forecasts produced by assimilation of radar observations to cloud-scale LES. Furthermore, to our knowledge, no storm-scale DA study has explored the value of assimilating high-resolution sector RHI volumes to determine if they can improve the retrieval of the vertical kinematic and microphysical storm structure aloft over the use of traditionally utilized PPI scans.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Reconciling Heliostat and Solar Panel Wind Loads To Develop a Reduced-Order Model for Heliostat Wind Loads

This work has been carried out as part of the HelioCon Field Deployment subtask, with the aim to distill published literature in papers and standards on heliostat and solar panel wind loads to develop a reduced order model for heliostat wind loads. It is intended to be understood and applied by structural engineers and heliostat designers as inputs to determine wind loads on heliostats due to the underlying derivations and factors that influence the complex structural dynamics and wind-induced behavior of heliostats.

14 SOLAR ENERGY↗
DOE OSTI2025

Parts per Billion Concentrations of Volatile Organic Gases collected with a PTRMS at the EPCAPE Mount Soledad site from 05-20-2023 to 06-15-2023

This dataset contains the parts per billion concentrations of various volatile organic gases. Data was collected using a Proton Transfer Reaction Mass Spectrometer (PTRMS). Concentrations for all compounds were calculated using a reaction constant (k-value) of 2.0E-9 cm^3/s. On June 14, a calibration was preformed for the gases methanol, acetone, isoprene, butanone, benzene and toluene. Calibrations curves were used to correct the concentrations of these gases. Negative concentration values were replaced with zeros. Other gases were not calibrated for. Data has undergone one minute averaging. Time data is in coordinated universal time (UTC). Time has been recorded in seconds since 1/1/1904.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Particle composition measurements for ultrafine particles collected at the EPCAPE Mount Soledad site from 04-27-2023 to 06-13-2023 using a Thermal Desorption Chemical Ionization Mass Spectrometer

The dataset contains particle composition data for both the positive and negative reagent ion modes of the Thermal Desorption Chemical Ionization Mass Spectrometer (TDCIMS). The dataset is split into two directories: one for particles with diameters of 30 nm and the other for particles with diameters of less than 100 nm. The positive reagent ion mode uses H3O+ as the reagent ion, and ionization usually occurs through hydrogen addition. The negative mode uses O2- as the reagent ion. Negative mode ionization generally occurs through hydrogen abstraction, but O2- addition is also possible. Ion concentrations were normalized to total ion counts, and unknown ions were then removed from the data. The name of each ion fraction time series includes the mass to charge ratio and the chemical formula for the ion. Time is recorded in seconds since 1/1/1904. The time zone is UTC.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Modeling Wind in Agrivoltaics and its Impact on Eddy Covariance Flux Measurements

Silicon Ranch is conducting a research study at their 138MW Bancroft Station agrivoltaics site, the largest agrivoltaics research array in the US, in which they are using an eddy covariance flux tower to measure the carbon budget. They are interested in finding out if the presence of the panels changes the effectiveness of these measurements. To investigate this topic, NLR performed numerical simulations of wind through the agrivoltaic array to identify under what environmental conditions (wind speed, wind direction, temperatures) and panel tilt angles the mean wind speeds at the measurement heights are significantly altered by the presence of the solar panels. NLR adapted and used the PVade simulation tool [1] to model wind in Silicon Ranch's agrivoltaic array and found that at the height of the flux tower (6m), the wind speed measured is not significantly altered by the presence of solar panels. The percent difference is less than 4% and is generally greatest during daytime conditions at higher wind speeds and perpendicular wind direction. Closer to the ground, the flow is significantly altered (wind speeds are reduced) by the presence of solar panels and the percent difference increases in strong winds. The work confirmed that eddy covariance flux towers can be used within the context of a solar array field.

14 SOLAR ENERGY↗
DOE OSTI2021

An Overview of Atmospheric Features Over the Western North Atlantic Ocean and North American East Coast—Part 2: Circulation, Boundary Layer, and Clouds

The Western North Atlantic Ocean (WNAO) is a complex land-ocean-atmosphere system that experiences a wide range of atmospheric phenomena, which in turn drive unique aerosol transport pathways, cloud morphologies, and boundary layer variability. This work, Part 2 of a 2-part paper series, provides an overview of the atmospheric circulation, boundary layer variability, three-dimensional cloud structure, and precipitation over the WNAO; the companion paper (Part 1) focused on chemical characterization of aerosol, gas, and wet deposition. Seasonal changes in the circulation and sea surface temperature explain a clear transition in cloud morphologies from small shallow cumulus clouds, convective clouds, and tropical storms in summer, to stratus/stratocumulus and multi-layer cloud systems associated with winter storms. Synoptic variability in cloud fields is estimated using satellite-based weather states, and the role of postfrontal conditions (cold-air outbreaks) in the development of stratiform clouds is further analyzed. Precipitation is persistent over the ocean, with a regional peak over the Gulf Stream path, where offshore sea surface temperature gradients are large and surface fluxes reach a regional peak. Satellite data show a clear annual cycle in cloud droplet number concentration along the coast in winter, and minima in summer, suggesting a marked annual cycle in aerosol-cloud interactions. The performance of four climate models are assessed against satellite cloud retrievals. The models qualitatively reproduce the annual cycle in cloud cover and liquid water path, but with large discrepancies across models, especially in the extra-tropics. The paper concludes with a summary of outstanding issues and recommendations for future work.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Developing and testing capabilities for simulating cases with heterogeneous land/water surfaces in a novel atmospheric large eddy simulation code

Large eddy simulations (LES) are the primary computational tool used to simulate high Reynolds number three-dimensional turbulent flows. In the context of earth system sciences, particularly atmospheric science, LES are uniquely able to resolve the scales of atmospheric motion that are key for building process-level understanding of boundary layer turbulence, atmosphere-surface interaction, clouds, and cloud-aerosol-chemistry interaction, and are a core limited-area modeling capability. Increasing demands are being placed on LES code bases as growing high performance computing resources allow LES to address a wider range of scientific problems. In addition, LES are emerging as a source of high-quality machine learning training data. These demands necessitate an agile and extensible code base that allows the model to quickly adapt to emergent needs. However, LES have largely relied on legacy Fortran code bases that lack flexibility. A new, Python-based LES capability called Predicting INteractions of Aerosol and Clouds in Large Eddy Simulation (PINACLES) has been developed as part of the Department of Energy’s Earth System Model Development (ESMD) program area’s Enabling Aerosol-cloud interactions at Global convection-permitting scalES (EAGLES) project. PINACLES was developed from the ground up with a philosophy of maximizing scientific throughput, by attempting to optimize for both model throughput and software extensibility. The initial development of PINACLES delivered a state-of-the-art idealized LES capability solving the non-hydrostatic anelastic equations of motion with doubly periodic boundary conditions and idealized homogenous surface boundary conditions. Here we provide a final report on the outcomes of a fiscal year 2021 Seed Laboratory Directed Research Project that extended PINACLES in two key ways. First, PINACLES was coupled to a state-of-the-art land surface model enabling it to simulate spatially inhomogeneous land-atmosphere interactions that are known to control key atmospheric processes. Second, the dynamical core of PINACLES was modified to permit non-periodic boundary conditions. This model enhancement enables simulation of realistic cases with boundary conditions prescribed from atmospheric reanalysis and enables nested simulations conducted on a hierarchy of computational domains with increasing resolution. Together, these extensions to PINACLES make it a formidable modeling capability and expand its potential application to diverse components of DOE’s atmospheric science portfolio.

42 ENGINEERING↗
DOE OSTI2022

A Case Study of the Weather Research and Forecasting Model Applied to the Joint Urban 2003 Tracer Field Experiment. Part III: Boundary-Layer Parametrizations

Numerical weather prediction is often used to supply the mean wind and turbulence fields for atmospheric transport and dispersion plume models as they provide dense geographic coverage in comparison to typically sparse monitoring networks. Here, the Weather Research and Forecasting (WRF) model 4.0 was run over the month-long period of the Joint Urban 2003 field campaign conducted in Oklahoma City. We compare three different simulations in their ability to reproduce the observations, each using a different boundary-layer parametrization. Specifically, we examine the Mellor–Yamada–Janjic (MYJ), Yonsei University (YSU), and Mellor–Yamada–Nakanishi–Niino (MYNN) boundary-layer parametrizations. All three predict the wind speed well during the day but overpredict it at night. The MYNN parametrization is better than MYJ at predicting the daytime turbulence in the surface layer, but both underpredict the nocturnal turbulence. Additionally, the MYJ parametrization is best at predicting the reciprocal Obukhov length, while MYNN and YSU both significantly overpredict thermal stability. Reconstructing the reciprocal Obukhov length from other simulated parameters produces more accurate values for both parametrizations. All three models overpredict the boundary-layer height, particularly under convective conditions. The MYJ parametrization overestimates boundary-layer height the most, while YSU and MYNN have comparable performance with MYNN having an advantage in predicting the stable boundary-layer height. Several days were found where the WRF simulations predict significant deviations from the prevailing diurnal pattern in wind direction, which are not found in the observations.

54 ENVIRONMENTAL SCIENCES↗