Ensembles of numerical weather prediction (NWP) model predictions are used for a variety of forecasting applications. Such ensembles quantify the uncertainty of the prediction because the spread in the ensemble predictions is correlated to forecast uncertainty. For atmospheric transport and dispersion and wind energy applications in particular, the NWP ensemble spread should accurately represent uncertainty in the low-level mean wind. To adequately sample the probability density function (PDF) of the forecast atmospheric state, it is necessary to account for several sources of uncertainty. Limited computational resources constrain the size of ensembles, so choices must be made about which members to include. No known objective methodology exists to guide users in choosing which combinations of physics parameterizations to include in an NWP ensemble, however. This study presents such a methodology. The authors build an NWP ensemble using the Advanced Research Weather Research and Forecasting Model (ARW-WRF). This 24-member ensemble varies physics parameterizations for 18 randomly selected 48-h forecast periods in boreal summer 2009. Verification focuses on 2-m temperature and 10-m wind components at forecast lead times from 12 to 48 h. Various statistical guidance methods are employed for down-selection, calibration, and verification of the ensemble forecasts. The ensemble down-selection is accomplished with principal component analysis. The ensemble PDF is then statistically dressed, or calibrated, using Bayesian model averaging. The postprocessing techniques presented here result in a recommended down-selected ensemble that is about half the size of the original ensemble yet produces similar forecast performance, and still includes critical diversity in several types of physics schemes.

As integration of solar power into the national electric grid rapidly increases, it becomes imperative to improve forecasting of this highly variable renewable resource. Thus, a team of researchers from the public, private, and academic sectors partnered to develop and assess a new solar power forecasting system, Sun4Cast. The partnership focused on improving decision-making for utilities and independent system operators, ultimately resulting in improved grid stability and cost savings for consumers. The project followed a value chain approach to determine key research and technology needs to reach desired results. Sun4Cast integrates various forecasting technologies across a spectrum of temporal and spatial scales to predict surface solar irradiance. Anchoring the system is WRF-Solar, a version of the Weather Research and Forecasting (WRF) numerical weather prediction (NWP) model optimized for solar irradiance prediction. Forecasts from multiple NWP models are blended via the Dynamic Integrated Forecast (DICast) System, which forms the basis of the system beyond about 6 h. For short-range (0–6 h) forecasts, Sun4Cast leverages several observation-based nowcasting technologies. These technologies are blended via the Nowcasting Expert System Integrator (NESI). The NESI and DICast systems are subsequently blended to produce short- to midterm irradiance forecasts for solar array locations. The irradiance forecasts are translated into power with uncertainties quantified using an analog ensemble approach and are provided to the industry partners for real-time decision-making. The Sun4Cast system ran operationally throughout 2015 and results were assessed. This paper analyzes the collaborative design process, discusses the project results, and provides recommendations for best-practice solar forecasting.

No abstract

The Sun4Cast solar power forecasting system, designed to predict solar irradiance and power generation at solar farms, is composed of several component models operating on both the nowcasting (0–6 h) and day-ahead forecast horizons. The different nowcasting models include a statistical forecasting model (StatCast), two satellite-based forecasting models [the Cooperative Institute for Research in the Atmosphere Nowcast (CIRACast) and the Multisensor Advection-Diffusion Nowcast (MADCast)], and a numerical weather prediction model (WRF-Solar). It is important to better understand and assess the strengths and weaknesses of these short-range models to facilitate further improvements. To that end, each of these models, including four WRF-Solar configurations, was evaluated for four case days in April 2014. For each model, the 15-min average predicted global horizontal irradiance (GHI) was compared with GHI observations from a network of seven pyranometers operated by the Sacramento Municipal Utility District (SMUD) in California. Each case day represents a canonical sky-cover regime for the SMUD region and thus represents different modeling challenges. The analysis found that each of the nowcasting models perform better or worse for particular lead times and weather situations. StatCast performs best in clear skies and for 0–1-h forecasts; CIRACast and MADCast perform reasonably well when cloud fields are not rapidly growing or dissipating; and WRF-Solar, when configured with a high-spatial-resolution aerosol climatology and a shallow cumulus parameterization, generally performs well in all situations. Further research is needed to develop an optimal dynamic blending technique that provides a single best forecast to energy utility operators.

Solar power can provide substantial power supply to the grid; however, it is also a highly variable energy source due to changes in weather conditions, i.e. clouds, that can cause rapid changes in solar power output. Independent systems operators (ISOs) and regional transmission organizations (RTOs) monitor the demand load and direct power generation from utilities, define operating limits and create contingency plans to balance the load with the available power generation resources. ISOs, RTOs, and utilities will require solar irradiance forecasts to effectively and efficiently balance the energy grid as the penetration of solar power increases. This study presents a cloud regime-dependent short-range solar irradiance forecasting system to provide 15-min average clearness index forecasts for 15-min, 60-min, 120-min and 180-min lead-times. A k-means algorithm identifies the cloud regime based on surface weather observations and irradiance observations. Then, Artificial Neural Networks (ANNs) are trained to predict the clearness index. This regime-dependent system makes a more accurate deterministic forecast than a global ANN or clearness index persistence and produces more accurate predictions of expected irradiance variability than assuming climatological average variability.

Flows in the atmospheric boundary layer are turbulent, characterized by a large Reynolds number, the existence of a roughness sublayer and the absence of a well-defined viscous layer. Exchanges with the surface are therefore dominated by turbulent fluxes. In numerical models for atmospheric flows, turbulent fluxes must be specified at the surface; however, surface fluxes are not known a priori and therefore must be parametrized. Atmospheric flow models, including global circulation, limited area models, and large-eddy simulation, employ Monin–Obukhov similarity theory (MOST) to parametrize surface fluxes. The MOST approach is a semi-empirical formulation that accounts for atmospheric stability effects through universal stability functions. The stability functions are determined based on limited observations using simple regression as a function of the non-dimensional stability parameter representing a ratio of distance from the surface and the Obukhov length scale (Obukhov in Trudy Inst Theor Geofiz AN SSSR 1:95–115, 1946), $$z/L$$ z / L . However, simple regression cannot capture the relationship between governing parameters and surface-layer structure under the wide range of conditions to which MOST is commonly applied. We therefore develop, train, and test two machine-learning models, an artificial neural network (ANN) and random forest (RF), to estimate surface fluxes of momentum, sensible heat, and moisture based on surface and near-surface observations. To train and test these machine-learning algorithms, we use several years of observations from the Cabauw mast in the Netherlands and from the National Oceanic and Atmospheric Administration’s Field Research Division tower in Idaho. The RF and ANN models outperform MOST. Even when we train the RF and ANN on one set of data and apply them to the second set, they provide more accurate estimates of all of the fluxes compared to MOST. Estimates of sensible heat and moisture fluxes are significantly improved, and model interpretability techniques highlight the logical physical relationships we expect in surface-layer processes.

No abstract

Wildland fire decision support systems require accurate predictions of wildland fire spread. Fuel moisture content (FMC) is one of the important parameters controlling the rate of spread of wildland fire. However, dead FMC measurements are provided by a relatively sparse network of remote automatic weather stations (RAWS), while live FMC is relatively infrequently measured manually. We developed a high resolution, gridded, real-time FMC data sets that did not previously exist for assimilation into operational wildland fire prediction systems based on ML. We used surface observations of live and dead FMC to train machine learning models to estimate FMC based on satellite observations. Moderate Resolution Imaging Spectrometer Terra and Aqua reflectances are used to predict the live and dead FMC measured by the Wildland Fire Assessment System and RAWS). We evaluate multiple machine learning methods including multiple linear regression, random forests (RFs), gradient boosted regression and artificial neural networks. The models are trained to learn the relationships between the satellite reflectances, surface weather and soil moisture observations and FMC. After training on data corresponding to the temporally and spatially nearest grid points to the irregularly spaced surface FMC observations, the machine learning models could be applied to all grid cells for a gridded product over the Conterminous United States (CONUS). The results show generally that the rule-based approaches have the lowest errors likely due to the sharp decision boundaries among the predictors, and the RF approach that utilizes bagging to avoid over-fitting has the lowest error on the test dataset. The errors are typically between 25%−33% the typical variability of the FMC data, which indicate the skill of the RF in estimating the FMC based on satellite data and surface characteristics. The FMC gridded product based on the RF runs operationally daily over CONUS and can be assimilated into WRF-Fire for more accurate wildland fire spread predictions.