Tristan J. Shepherd

A 15-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of 29–30 June 2012 that impacted the northeastern urban corridor of the USA is presented. This event generated 1100 reports of damaging winds, generated significant wind gusts over an extensive area of up to 500 000 km², caused several fatalities, and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global-warming experiments to examine the sensitivity of historical, societally important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic contexts. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as is manifest in high composite reflectivity, cREF), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBCs), start date, use of nudging, compiler choice, damping, and number of vertical levels). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g., convective available potential energy and vertical wind shear) and critical manifestations of deep convection, e.g., vertical velocities, cold-pool generation, and how those properties relate to the correct characterization of the associated atmospheric hazards (wind gusts and hail). Use of a double-moment, seven-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model-simulated wind gusts and convective structure to the observations of this event. All ensemble members, however, fail to capture the intensity of the event in terms of the spatial extent of convection and the production of high near-surface wind gusts. We further show very high sensitivity to the LBCs employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim compared to ERA5. Excess convective available potential energy (CAPE) in all ensemble members after the derecho passage leads to excess production of convective cells, wind gusts, cREF > 40 dBZ, and precipitation during a frontal passage on the subsequent day. This event proved very challenging to forecast in real time and to reproduce in the 15-member hindcast simulation ensemble presented here. Future work could examine if simulations with other initial and lateral boundary conditions can achieve greater fidelity.

Projected changes to the El Niño Southern Oscillation (ENSO) climate mode have been explored using global Earth system models (ESMs). Regional expressions of such changes have yet to be fully advanced and may require the use of regional downscaling. Here, we employ regional climate modeling (RCM) using the Weather Research and Forecasting (WRF) model at convection-permitting resolution and nested in output from the HadGEM2 ESM. We quantify ENSO teleconnections to temperature and precipitation anomalies in historical and future climate scenarios over eastern North America. Two paired simulations are run, a strong El Niño (positive ENSO phase) and a weak La Niña (negative ENSO phase), for the historical and future years. The HadGEM2 direct output and HadGEM2-WRF simulation output are compared to the anomalies derived from the NOAA ENSO Climate Normals dataset. The near-surface temperature and precipitation differences by ENSO phase, as represented by the HadGEM2-WRF historical simulations, show a poor degree of association with the NOAA ENSO Climate Normals, in part because of the large biases in the HadGEM2 model. Downscaling with the WRF model does improve the agreement with the observations, and large discrepancies remain. The model chain HadGEM2-WRF reverses the sign of the ENSO phase response over eastern North America under simulations of the future climate with high greenhouse gas forcing, but due to the poor agreement with the NOAA ENSO Climate Normals it is difficult to assign confidence to this prediction.

An 11-member ensemble of convection-permitting regional simulations of the fast-moving and destructive derecho of June 29 – 30, 2012 that impacted the northeastern urban corridor of the US is presented. This event generated 1100 reports of damaging winds, significant wind gusts over an extensive area of up to 500,000 km², caused several fatalities and resulted in widespread loss of electrical power. Extreme events such as this are increasingly being used within pseudo-global warming experiments that seek to examine the sensitivity of historical, societally-important events to global climate non-stationarity and how they may evolve as a result of changing thermodynamic and dynamic context. As such it is important to examine the fidelity with which such events are described in hindcast experiments. The regional simulations presented herein are performed using the Weather Research and Forecasting (WRF) model. The resulting ensemble is used to explore simulation fidelity relative to observations for wind gust magnitudes, spatial scales of convection (as manifest in high composite reflectivity), and both rainfall and hail production as a function of model configuration (microphysics parameterization, lateral boundary conditions (LBC), start date, and use of nudging). We also examine the degree to which each ensemble member differs with respect to key mesoscale drivers of convective systems (e.g. convective available potential energy and vertical wind shear) and critical manifestations of deep convection; e.g. vertical velocities, cold pool generation, and how those properties relate to correct characterization of the associated atmospheric hazards (wind gusts and hail). Here, we show that the use of a double-moment, 7-class scheme with number concentrations for all species (including hail and graupel) results in the greatest fidelity of model simulated wind gusts and convective structure against the observations of this event. We further show very high sensitivity to the LBC employed and specifically that simulation fidelity is higher for simulations nested within ERA-Interim than ERA5.

Output from 6 months of high-resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low-level jets (LLJs) over Iowa for winter and spring in the contemporary climate. Low-level jets affect rotor plane aerodynamic loading, turbine structural loading and turbine performance, and thus accurate characterization and identification are pertinent. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20–530 m a.g.l. (above ground level) indicate the presence of an LLJ in at least one of the 14 700 4 km×4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJs are most frequently associated with stable stratification and low turbulent kinetic energy (TKE) and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 m s⁻¹, and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying, with a mean duration of 3.5 h. The highest frequency occurs in the topographically complex northwest of the state in winter and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the (i) LLJ definition and (ii) vertical resolution at which the WRF output is sampled is examined. LLJ definitions commonly used in the literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with definition. Use of different definitions identifies both different frequencies of LLJs and different LLJ events. Further, when the model output is down-sampled to lower vertical resolution, the mean jet core wind speed height decreases, but spatial distributions of regions of high frequency and duration are conserved. Implementation of a polynomial interpolation to extrapolate down-sampled output to full-resolution results in reduced sensitivity of LLJ characteristics to down-sampling.

Output from high resolution simulations with the Weather Research and Forecasting (WRF) model are analyzed to characterize local low level jets (LLJ) over Iowa. Analyses using a detection algorithm wherein the wind speed above and below the jet maximum must be below 80 % of the jet wind speed within a vertical window of approximately 20 m–530 m a.g.l. indicate the presence of a LLJ in at least one of the 14700 4 km by 4 km grid cells over Iowa on 98 % of nights. Nocturnal LLJ are most frequently associated with stable stratification and low TKE and hence are more frequent during the winter months. The spatiotemporal mean LLJ maximum (jet core) wind speed is 9.55 ms⁻¹ and the mean height is 182 m. Locations of high LLJ frequency and duration across the state are seasonally varying with a mean duration of 3.5 hours. LLJ are most frequent in the topographically complex northwest of the state in winter, and in the flatter northeast of the state in spring. Sensitivity of LLJ characteristics to the: i) LLJ definition and ii) vertical resolution at which the WRF output is sampled are examined. LLJ definitions commonly used in LLJ literature are considered in the first sensitivity analysis. These sensitivity analyses indicate that LLJ characteristics are highly variable with LLJ definition. Further, when the model output is down-sampled to lower vertical resolution, the maximum LLJ wind speed and mean height decrease, but spatial distributions of regions of high frequency and duration are conserved.

Continued growth of wind turbine physical dimensions is examined in terms of the implications for wind speed, power and shear across the rotor plane. High-resolution simulations with the Weather Research and Forecasting model are used to generate statistics of wind speed profiles for scenarios of current and future wind turbines. The nine-month simulations, focused on the eastern Central Plains, show that the power scales broadly as expected with the increase in rotor diameter (D) and wind speeds at hub-height (H). Increasing wind turbine dimensions from current values (approximately H = 100 m, D = 100 m) to those of the new International Energy Agency reference wind turbine (H = 150 m, D = 240 m), the power across the rotor plane increases 7.1 times. The mean domain-wide wind shear exponent (α) decreases from 0.21 (H = 100 m, D = 100 m) to 0.19 for the largest wind turbine scenario considered (H = 168 m, D = 248 m) and the frequency of extreme positive shear (α > 0.2) declines from 48% to 38% of 10-min periods. Thus, deployment of larger wind turbines potentially yields considerable net benefits for both the wind resource and reductions in fatigue loading related to vertical shear.

The Weather Research and Forecasting (WRF) Model has been extensively used for wind energy applications, and current releases include a scheme that can be applied to examine the effects of wind turbine arrays on the atmospheric flow and electricity generation from wind turbines. Herein we present a high-resolution simulation using two different wind farm parameterizations: 1) the “Fitch” parameterization that is included in WRF releases and 2) the recently developed Explicit Wake Parameterization (EWP) scheme. We compare the schemes using a single yearlong simulation for a domain centered on the highest density of current turbine deployments in the contiguous United States (Iowa). Pairwise analyses are applied to diagnose the downstream wake effects and impact of wind turbine arrays on near-surface climate conditions. On average, use of the EWP scheme results in small-magnitude wake effects within wind farm arrays and faster recovery of full WT array wakes. This in turn leads to smaller impacts on near-surface climate variables and reduced array–array interactions, which at a systemwide scale lead to summertime capacity factors (i.e., the electrical power produced relative to nameplate installed capacity) that are 2%–3% higher than those from the more commonly applied Fitch parameterization. It is currently not possible to make recommendations with regard to which wind farm parameterization exhibits higher fidelity or to draw inferences with regard to whether the relative performance may vary with prevailing climate conditions and/or wind turbine deployment configuration. However, the sensitivities documented herein to the wind farm parameterization are of sufficient magnitude to potentially influence wind turbine array siting decisions. Thus, our research findings imply high value in undertaking combined long-term high-fidelity observational studies in support of model validation and verification.

Impacts from current and future wind turbine (WT) deployments necessary to achieve 20% electricity from wind are analyzed using high resolution numerical simulations over the eastern USA. Theoretical scenarios for future deployments are based on repowering (i.e. replacing with higher capacity WTs) thus avoiding competition for land. Simulations for the contemporary climate and current WT deployments exhibit good agreement with observed electricity generation efficiency (gross capacity factors (CF) from simulations = 45–48%, while net CF for WT installed in 2016 = 42.5%). Under the scenario of quadrupled installed capacity there is a small decrease in system-wide efficiency as indicated by annual mean CF. This difference is approximately equal to that from the two simulation years and may reflect saturation of the wind resource in some areas. WT modify the local near-surface climate in the grid cells where they are deployed. The simulated impact on near-surface climate properties at both the regional and local scales does not increase with increasing WT installed capacity. Climate impacts from WT are modest compared to regional changes induced by historical changes in land cover and to the global temperature perturbation induced by use of coal to generate an equivalent amount of electricity.

High-resolution simulations are conducted with the Weather Research and Forecasting Model to evaluate the sensitivity of wake effects and power production from two wind farm parameterizations [the commonly used Fitch scheme and the more recently developed Explicit Wake Parameterization (EWP)] to the resolution at which the model is applied. The simulations are conducted for a 9-month period for a domain encompassing much of the U.S. Midwest. The two horizontal resolutions considered are 4 km × 4 km and 2 km × 2 km grid cells, and the two vertical discretizations employ either 41 or 57 vertical layers (with the latter having double the number in the lowest 1 km). Higher wind speeds are observed close to the wind turbine hub height when a larger number of vertical layers are employed (12 in the lowest 200 m vs 6), which contributes to higher power production from both wind farm schemes. Differences in gross capacity factors for wind turbine power production from the two wind farm parameterizations and with resolution are most strongly manifest under stable conditions (i.e., at night). The spatial extent of wind farm wakes when defined as the area affected by velocity deficits near to wind turbine hub heights in excess of 2% of the simulation without wind turbines is considerably larger in simulations with the Fitch scheme. This spatial extent is generally reduced by increasing the horizontal resolution and/or increasing the number of vertical levels. These results have important applications to projections of expected annual energy production from new wind turbine arrays constructed in the wind shadow from existing wind farms.