A grand challenge from the wind energy industry is to provide reliable forecasts on mountain winds several hours in advance at microscale (∼100 m) resolution. This requires better microscale wind-energy physics included in forecasting tools, for which field observations are imperative. While mesoscale (∼1 km) measurements abound, microscale processes are not monitored in practice nor do plentiful measurements exist at this scale. After a decade of preparation, a group of European and U.S. collaborators conducted a field campaign during 1 May–15 June 2017 in Vale Cobrão in central Portugal to delve into microscale processes in complex terrain. This valley is nestled within a parallel double ridge near the town of Perdigão with dominant wind climatology normal to the ridges, offering a nominally simple yet natural setting for fundamental studies. The dense instrument ensemble deployed covered a ∼4 km × 4 km swath horizontally and ∼10 km vertically, with measurement resolutions of tens of meters and seconds. Meteorological data were collected continuously, capturing multiscale flow interactions from synoptic to microscales, diurnal variability, thermal circulation, turbine wake and acoustics, waves, and turbulence. Particularly noteworthy are the extensiveness of the instrument array, space–time scales covered, use of leading-edge multiple-lidar technology alongside conventional tower and remote sensors, fruitful cross-Atlantic partnership, and adaptive management of the campaign. Preliminary data analysis uncovered interesting new phenomena. All data are being archived for public use.

Accurately representing flow across the mesoscale to the microscale is a persistent roadblock for completing realistic microscale simulations. The science challenges that must be addressed to coupling at these scales include the following: 1) What is necessary to capture the variability of the mesoscale flow, and how do we avoid generating spurious rolls within the terra incognita between the scales? 2) Which methods effectively couple the mesoscale to the microscale and capture the correct nonstationary features at the microscale? 3) What are the best methods to initialize turbulence at the microscale? 4) What is the best way to handle the surface-layer parameterizations consistently at the mesoscale and the microscale? 5) How do we assess the impact of improvements in each of these aspects and quantify the uncertainty in the simulations? The U.S. Department of Energy Mesoscale-to-Microscale-Coupling project seeks to develop, verify, and validate physical models and modeling techniques that bridge the most important atmospheric scales determining wind plant performance and reliability, which impacts many meteorological applications. The approach begins with choosing case days that are interesting for wind energy for which there are observational data for validation. The team has focused on modeling nonstationary conditions for both flat and complex terrain. This paper describes the approaches taken to answer the science challenges, culminating in recommendations for best approaches for coupled modeling.

Multiscale atmospheric simulations can be computationally prohibitive, as they require large domains and fine spatiotemporal resolutions. Grid‐nesting can alleviate this by bridging mesoscales and microscales, but one turbulence scheme must run at resolutions within a range of scales known as the terra incognita (TI). TI grid‐cell sizes can violate both mesoscale and microscale subgrid‐scale parametrization assumptions, resulting in unrealistic flow structures. Herein we assess the impact of unrealistic lateral boundary conditions from parent mesoscale simulations at TI resolutions on nested large eddy simulations (LES), to determine whether parent domains bias the nested LES. We present a series of idealized nested mesoscale‐to‐LES runs of a dry convective boundary layer (CBL) with different parent resolutions in the TI. We compare the nested LES with a stand‐alone LES with periodic boundary conditions. The nested LES domains develop ∼20% smaller convective structures, while potential temperature profiles are nearly identical for both the mesoscales and LES simulations. The horizontal wind speed and surface wind shear in the nested simulations closely resemble the reference LES. Heat fluxes are overestimated by up to ∼0.01 K m s−1 in the top half of the PBL for all nested simulations. Overestimates of turbulent kinetic energy (TKE) and Reynolds stress in the nested domains are proportional to the parent domain's grid‐cell size, and are almost eliminated for the simulation with the finest parent grid‐cell size. Based on these results, we recommend that LES of the CBL be forced by mesoscale simulations with the finest practical resolution.

Coupling between mesoscale and large eddy simulation (LES) is critically important for many atmospheric model applications, from predictions of wind energy to fire propagation. The grid-nesting technique enables bridging between vastly different scales without incurring prohibitive computational costs. However, the transition from coarser to finer resolutions often requires a large number of grid points from inflow boundaries for the development of fine-scale turbulence features in the LES domain. Recently, the cell perturbation method (CPM) was developed to reduce the turbulence development region with high computational efficiency. Herein, we explore a new method based on the CPM that uses force perturbations in both the horizontal and vertical directions (Force Cell Perturbation Method) instead of the potential temperature perturbations in the original CPM, as an attempt to further explore the performance of the random perturbation techniques. This approach is tested for a neutral and a convective atmospheric boundary layer under idealized conditions. Overall, similar performance is found between the optimal configurations of the CPM and the Force Cell Perturbation Method pointing to the robustness of this family of perturbation methods in accelerating turbulence generation in nested domains. Vertical force perturbations performed better than horizontal force perturbations for both atmospheric stability conditions. The CPM performed best under convective stability conditions. The combination of the force and potential temperature perturbations is found to provide no additional performance improvement over the stand-alone application of the individual methods.