The second Wind Forecast Improvement Project (WFIP2) is a multiagency field campaign held in the Columbia Gorge area (October 2015–March 2017). The main goal of the project is to understand and improve the forecast skill of numerical weather prediction (NWP) models in complex terrain, particularly beneficial for the wind energy industry. This region is well known for its excellent wind resource. One of the biggest challenges for wind power production is the accurate forecasting of wind ramp events (large changes of generated power over short periods of time). Poor forecasting of the ramps requires large and sudden adjustments in conventional power generation, ultimately increasing the costs of power. A Ramp Tool and Metric (RT&M) was developed during the first WFIP experiment, held in the U.S. Great Plains (September 2011–August 2012). The RT&M was designed to explicitly measure the skill of NWP models at forecasting wind ramp events. Here we apply the RT&M to 80-m (turbine hub-height) wind speeds measured by 19 sodars and three lidars, and to forecasts from the High-Resolution Rapid Refresh (HRRR), 3-km, and from the High-Resolution Rapid Refresh Nest (HRRRNEST), 750-m horizontal grid spacing, models. The diurnal and seasonal distribution of ramp events are analyzed, finding a noticeable diurnal variability for spring and summer but less for fall and especially winter. Also, winter has fewer ramps compared to the other seasons. The model skill at forecasting ramp events, including the impact of the modification to the model physical parameterizations, was finally investigated.
This work details a master oscillator power amplifier (MOPA) microjoule-class pulsed coherent Doppler lidar system configuration designed to measure line-of-sight wind velocities and backscatter intensity of atmospheric aerosols. The instrument is unique in its form factor. It consists of two physically separated modules connected by a 10 m umbilical cable. One module hosts the transceiver, which is composed of the telescope, transmit/receive (T/R) switch, and high-gain optical amplifier, and is housed in a small box (34.3 cm X 34.3 cm X 17.8 cm). The second module contains the data acquisition system and several electro-optical components. This form factor enables deployments on platforms that are otherwise inaccessible by commercial and research instruments of similar design. In this work, optical, electrical, and data acquisition components and configurations of the lidar are detailed and two example deployments are presented. The first deployment describes measurements of a controlled wildfire burn from a small aircraft to measure vertical plume dynamics and fire inflow conditions during summer in Florida. The second presents measurements of the marine boundary layer height and vertical velocity and variance profiles from the Research Vessel (R/V) Thomas Thompson. The new instrument has enabled greater flexibility in field campaigns where previous instruments would have been too costly or space prohibitive to deploy.
The spatiotemporal scales of turbulence in the boundary layer flow vary with the static stability of the atmosphere. For example, turbulence structures in unstable conditions are often found to be large and coherent, whereas the turbulent structures in stable conditions are relatively small and sporadic. To better understand the behavior of such turbulent structures under different atmospheric stability conditions in the lower atmosphere, herein, we use scanning lidar data collected over several months during the US Department of Energy's second Wind Forecast Improvement Project (WFIP 2) in 2016. The radial velocity data from scanning lidar were gathered along the conical surface (PPI mode) and were analyzed using proper orthogonal decomposition (POD) and spectral energy methods. The results from these analyses show that the turbulent structures in daytime unstable conditions are different from those present during nighttime stable conditions in terms of amount of energy and coherent structures. The POD energy is distributed more evenly in all the higher spatial modes in the daytime unstable case compared to the nighttime stable case. Besides characterizing the turbulence in the flow, the lidar data were also used to compare the simulated flow field generated from the WRF-LES framework (Weather Research and Forecasting - large-eddy simulation) with the scanning lidar data. The results from both simulated and lidar scanned data show similar structures - orienting the streak-like structure along the mean wind. This indicates that the streak-like structures can exist in the real atmosphere during unstable atmospheric condition cases.
Recently, two air quality campaigns were conducted in the southwestern United States to study the impact of transported ozone, stratospheric intrusions, and fire emissions on ground-level ozone concentrations. The California Baseline Ozone Transport Study (CABOTS) took place in May – August 2016 covering the central California coast and San Joaquin Valley, and the Fires, Asian, and Stratospheric Transport Las Vegas Ozone Study (FAST-LVOS) was conducted in the greater Las Vegas, Nevada area in May – June 2017. During these studies, nearly 1000 hours of ozone and aerosol profile data were collected with the NOAA TOPAZ lidar. A Doppler wind lidar and a radar wind profiler provided continuous observations of atmospheric turbulence, horizontal winds, and mixed layer height. These measurements allowed us to directly observe the degree to which ozone transport layers aloft were entrained into the boundary layer and to quantify the resulting impact on surface ozone levels. Mixed layer heights in the San Joaquin Valley during CABOTS were generally below 1 km above ground level (AGL), while boundary layer heights in Las Vegas during FAST-LVOS routinely exceeded 3 km AGL and occasionally reached up to 4.5 km AGL. Consequently, boundary layer entrainment was more often observed during FAST-LVOS, while most elevated ozone layers passed untapped over the San Joaquin Valley during CABOTS.
Urban wind island effect (UWI) is defined as a phenomenon in which boundary layer mean wind speeds in an urban area are noticeably higher than its neighboring rural areas. Unlike urban heat island effect which has been extensively studied, the UWI was only recently observed in a modeling study. Here we study existence of the UWI over Indianapolis, Indiana using wind profile measurements from two Doppler wind lidars (DWL) that were deployed in climatologically upwind and downwind of the city. Under certain atmospheric conditions higher wind speeds and turbulence were observed at the downwind site over the entire urban boundary layer outside the urban canopy layer.
In the present paper, we assessed the potential for local wind, solar PV, and energy storage to provide baseload (constant, uninterrupted) power in every county of the contiguous United States. The amount of available capacity between 2020 and 2050 was determined via a least-cost optimization model that took into account changing costs of constituent technologies and local meteorological conditions. We found that, by 2050, the potential exists for about 6.8 TW of renewable baseload power at an average cost of approximately $50 / MWh, which is competitive with current wholesale market rates for electricity. The optimal technology configurations constructed always resulted in over two hours of emergency energy reserves, with the amount increasing as the price of energy storage falls. We also found that, given current price decline trajectories, the model has a tendency to select more solar capacity than wind over time. A second part of the study performed three million simulations followed by a regression analysis to generate an online map-based tool that allows users to change input costs assumptions and compute the cost of renewable baseload electricity in every contiguous US county.
During the second Wind Forecast Improvement Project (WFIP2; October 2015–March 2017, held in the Columbia River Gorge and Basin area of eastern Washington and Oregon states), several improvements to the parameterizations used in the High Resolution Rapid Refresh (HRRR – 3 km horizontal grid spacing) and the High Resolution Rapid Refresh Nest (HRRRNEST – 750 m horizontal grid spacing) numerical weather prediction (NWP) models were tested during four 6-week reforecast periods (one for each season). For these tests the models were run in control (CNT) and experimental (EXP) configurations, with the EXP configuration including all the improved parameterizations. The impacts of the experimental parameterizations on the forecast of 80 m wind speeds (wind turbine hub height) from the HRRR and HRRRNEST models are assessed, using observations collected by 19 sodars and three profiling lidars for comparison. Improvements due to the experimental physics (EXP vs. CNT runs) and those due to finer horizontal grid spacing (HRRRNEST vs. HRRR) and the combination of the two are compared, using standard bulk statistics such as mean absolute error (MAE) and mean bias error (bias). On average, the HRRR 80 m wind speed MAE is reduced by 3 %–4 % due to the experimental physics. The impact of the finer horizontal grid spacing in the CNT runs also shows a positive improvement of 5 % on MAE, which is particularly large at nighttime and during the morning transition. Lastly, the combined impact of the experimental physics and finer horizontal grid spacing produces larger improvements in the 80 m wind speed MAE, up to 7 %–8 %. The improvements are evaluated as a function of the model's initialization time, forecast horizon, time of the day, season of the year, site elevation, and meteorological phenomena. Causes of model weaknesses are identified. Finally, bias correction methods are applied to the 80 m wind speed model outputs to measure their impact on the improvements due to the removal of the systematic component of the errors.
To improve parameterizations of the turbulence dissipation rate (ϵ) in numerical weather prediction models, the temporal and spatial variability of ϵ must be assessed. In this study, we explore influences on the variability of ϵ at various scales in the Columbia River Gorge during the WFIP2 field experiment between 2015 and 2017. We calculate ϵ from five sonic anemometers all deployed in a ∼4 km2 area as well as from two scanning Doppler lidars and four profiling Doppler lidars, whose locations span a ∼300 km wide region. We retrieve ϵ from the sonic anemometers using the second-order structure function method, from the scanning lidars with the azimuth structure function approach, and from the profiling lidars with a novel technique using the variance of the line-of-sight velocity. The turbulence dissipation rate shows large spatial variability, even at the microscale, especially during nighttime stable conditions. Orographic features have a strong impact on the variability of ϵ, with the correlation between ϵ at different stations being highly influenced by terrain. ϵ shows larger values in sites located downwind of complex orographic structures or in wind farm wakes. A clear diurnal cycle in ϵ is found, with daytime convective conditions determining values over an order of magnitude higher than nighttime stable conditions. ϵ also shows a distinct seasonal cycle, with differences greater than an order of magnitude between average ϵ values in summer and winter.
The Wind Forecast Improvement Project 2 (WFIP2), in the complex terrain of the Columbia Gorge, focuses on a set of weather phenomena that poses particular challenges for wind and wind power forecasting. It further aims to understand and improve the skill of weather forecast models, particularly the National Oceanic and Atmospheric Administration's (NOAA's) High Resolution Rapid Refresh (HRRR), in complex terrain. The project also includes an extensive field campaign in the Columbia River Gorge, during which data from many different instruments were collected. WFIP2 is funded by the U.S. Department of Energy (DOE). Scientists from four U.S. national laboratories, NOAA, Vaisala, and universities contribute to the project with wide and varied interests and skill sets. Therefore, coordinated verification & validation (V&V) efforts across these members allows for the development of a clear picture of model improvements and scientific findings within WFIP2. The WFIP2 V&V team is tasked with providing tools, methods, and guidance to enable repeatable, metrics-based assessment of the Weather Research and Forecasting (WRF) model and associated modeling suites for analysis and forecasting of mesoscale weather phenomena that are important for wind energy in the Columbia River Gorge and other parts of the continental United States. This report summarizes what has been worked on and accomplished within the duration of the WFIP2 project, 4 years, by the V&V team.
Annually and seasonally averaged wind profiles from three Doppler lidars were obtained from sites in the Columbia River basin of east-central Oregon and Washington, a major region of wind-energy production, for the Second Wind Forecast Improvement Project (WFIP2) experiment. The profile data are used to quantify the spatial variability of wind flows in this area of complex terrain, to assess the HRRR-NCEP model's ability to capture spatial and temporal variability of wind profiles, and to evaluate model errors. Annually averaged measured wind speed differences over the 70-km extent of the lidar measurements reached 1 m s(-1) within the wind-turbine rotor layer, and 2 m s(-1) for 200-500 m AGL. Stronger wind speeds in the lowest 500 m occurred at sites higher in elevation, farther from the river, and farther west-closer to the Cascade Mountain barrier. Validating against the lidar data, the HRRR model underestimated strong wind speeds (>12 m s(-1)) and, consequently, their frequency of occurrence, especially at the two lowest-elevation sites, producing annual low biases in rotor-layer wind speed of 0.5 m s(-1). The RMSE between measured and modeled winds at all sites was about 3 m s(-1) and did not degrade significantly with forecast lead time. The nature of the model errors was different for different seasons. Moreover, although the three sites were located in the same basin terrain, the nature of the model errors was different at each site. Thus, if only one of the sites had been instrumented, different conclusions would have been drawn as to the major sources of model error, depending on where the measurements were made.
The Second Wind Forecast Improvement Project (WFIP2) is a U.S. Department of Energy (DOE)- and National Oceanic and Atmospheric Administration (NOAA)-funded program, with private-sector and university partners, which aims to improve the accuracy of numerical weather prediction (NWP) model forecasts of wind speed in complex terrain for wind energy applications. A core component of WFIP2 was an 18-month field campaign that took place in the U.S. Pacific Northwest between October 2015 and March 2017. A large suite of instrumentation was deployed in a series of telescoping arrays, ranging from 500 km across to a densely instrumented 2 km x 2 km area similar in size to a high-resolution NWP model grid cell. Observations from these instruments are being used to improve our understanding of the meteorological phenomena that affect wind energy production in complex terrain and to evaluate and improve model physical parameterization schemes. We present several brief case studies using these observations to describe phenomena that are routinely difficult to forecast, including wintertime cold pools, diurnally driven gap flows, and mountain waves/wakes. Observing system and data product improvements developed during WFIP2 are also described.
Ground-based Doppler-lidar instrumentation provides atmospheric wind data at dramatically improved accuracies and spatial/temporal resolutions. These capabilities have provided new insights into atmospheric flow phenomena, but they also should have a strong role in NWP model improvement. Insight into the nature of model errors can be gained by studying recurrent atmospheric flows, here a regional summertime diurnal sea breeze and subsequent marine-air intrusion into the arid interior of Oregon–Washington, where these winds are an important wind-energy resource. These marine intrusions were sampled by three scanning Doppler lidars in the Columbia River basin as part of the Second Wind Forecast Improvement Project (WFIP2), using data from summer 2016. Lidar time–height cross sections of wind speed identified 8 days when the diurnal flow cycle (peak wind speeds at midnight, afternoon minima) was obvious and strong. The 8-day composite time–height cross sections of lidar wind speeds are used to validate those generated by the operational NCEP–HRRR model. HRRR simulated the diurnal wind cycle, but produced errors in the timing of onset and significant errors due to a premature nighttime demise of the intrusion flow, producing low-bias errors of 6 m s−1. Day-to-day and in the composite, whenever a marine intrusion occurred, HRRR made these same errors. The errors occurred under a range of gradient wind conditions indicating that they resulted from the misrepresentation of physical processes within a limited region around the measurement locations. Because of their generation within a limited geographical area, field measurement programs can be designed to find and address the sources of these NWP errors.