Weather Research And Forecasting Ensemble Research Articles

Wind Speed Ramp Rate Predictions Using Wind Farm SCADA Data Assimilation and a WRF Ensemble

This study presents a novel method for improving wind power ramp events forecasts up to six hours ahead by utilizing data assimilation of SCADA measurements with an ensemble of Weather Research and Forecasting (WRF) models estimates. Leveraging data from nine wind farms in France and Belgium, the approach aims to improve WRF model predictions for wind speed and ramp event timing. The methodology employs grid and observational nudging techniques, enhancing model accuracy by incorporating real-time observational data. Key findings demonstrate that nudging significantly reduces Mean Absolute Error (MAE), decreases the Time Distortion Index (TDI), and increases the Probability of Detection (POD) of ramp events. Nudged ensemble members outperform non-nudged counterparts, exhibiting better accuracy in identifying true ramp events and reducing false alarms. MAE, TDI and POD improvements are as high as 3.7%, 8.5% and 37%, respectively. The study also explores the benefits of an ensemble approach, highlighting improved accuracy in predicting ramp rate magnitudes and providing valuable insights for grid stability management. This research contributes to wind power forecasting, showcasing the importance of integrating SCADA data into predictive models.

A Flood Forecasting Framework Coupling a High Resolution WRF Ensemble With an Urban Hydrologic Model

Numerical weather prediction (NWP) models such as the Weather Research and Forecasting (WRF) model are increasingly used over the Indian region to forecast extreme rainfall events. However, studies which explore the application of high-resolution rainfall simulations obtained from the WRF model in urban hydrology are limited. In this paper, the utility of a model coupling framework to predict urban floods is explored through the case study of Bangalore city in India. This framework is used to simulate multiple extreme events that occurred over the city for the monsoons of years 2020 and 2021. To address the uncertainty from the WRF model, a 12-member convection permitting ensemble is used. Model configurations using Kain Fritsch and WSM6 parameterization schemes could simulate the spatial and temporal pattern of the selected event. The city is easily flooded with rainfall events above a threshold of 60 mm/day and to capture the response of the urban catchment, the Personal Computer Storm Water Management Model (PCSWMM) is used in this study. Flood forecasts are created using the outputs from the WRF ensemble and the Global Forecasting System (GFS). The high temporal and spatial resolution of the rainfall forecasts (<4 km at 15-min intervals), has proved critical in reproducing the urban flood event. The flood forecasts created using the WRF ensemble indicate that flooding and water levels are comparable to the observed whereas the GFS underestimates these to a large extent. Thus, the coupled WRF–PCSWMM modelling framework is found effective in forecasting flood events over an Indian city.

Objective identification of thunderstorm gust fronts in numerical weather prediction models for fire weather forecasting

Abrupt changes in wind direction and speed can dramatically impact wildfire development and spread, endangering firefighters. A frequent cause of such wind shifts is outflow from thunderstorms and organised convective systems; thus, their identification and prediction present critical challenges for fire weather forecasters. Here, we develop a methodology and implement it in a software tool that can identify and depict convective outflow boundaries in high-resolution numerical weather prediction (NWP) models to provide guidance for fire weather forecasting. The tool can process model output, objectively identify gust fronts, and graphically display detected gust fronts and similar boundaries in NWP model forecasts. The tool is demonstrated with output from the Weather Research and Forecasting (WRF) model from the operational High-Resolution Rapid Refresh (HRRR) forecasting system and from a WRF ensemble run at the National Center for Atmospheric Research that can provide probabilistic information about model-predicted gust fronts. The tool can identify outflow boundaries in model forecasts of convective events occurring in both simple and complex terrain, both with and without concurrent wildfire activity. With accurate underlying model forecast output, the tool can reliably reveal areas of potential gust front activity and thus provide valuable guidance to incident meteorologists and command personnel.

The making of the New European Wind Atlas – Part 1: Model sensitivity

Abstract. This is the first of two papers that document the creation of the New European Wind Atlas (NEWA). It describes the sensitivity analysis and evaluation procedures that formed the basis for choosing the final setup of the mesoscale model simulations of the wind atlas. The suitable combination of model setup and parameterizations, bound by practical constraints, was found for simulating the climatology of the wind field at turbine-relevant heights with the Weather Research and Forecasting (WRF) model. Initial WRF model sensitivity experiments compared the wind climate generated by using two commonly used planetary boundary layer schemes and were carried out over several regions in Europe. They confirmed that the most significant differences in annual mean wind speed at 100 m a.g.l. (above ground level) mostly coincide with areas of high surface roughness length and not with the location of the domains or maximum wind speed. Then an ensemble of more than 50 simulations with different setups for a single year was carried out for one domain covering northern Europe for which tall mast observations were available. We varied many different parameters across the simulations, e.g. model version, forcing data, various physical parameterizations, and the size of the model domain. These simulations showed that although virtually every parameter change affects the results in some way, significant changes in the wind climate in the boundary layer are mostly due to using different physical parameterizations, especially the planetary boundary layer scheme, the representation of the land surface, and the prescribed surface roughness length. Also, the setup of the simulations, such as the integration length and the domain size, can considerably influence the results. We assessed the degree of similarity between winds simulated by the WRF ensemble members and the observations using a suite of metrics, including the Earth Mover's Distance (EMD), a statistic that measures the distance between two probability distributions. The EMD was used to diagnose the performance of each ensemble member using the full wind speed and direction distribution, which is essential for wind resource assessment. We identified the most realistic ensemble members to determine the most suitable configuration to be used in the final production run, which is fully described and evaluated in the second part of this study (Dörenkämper et al., 2020).

Validation of Snow Multibands in the Comma Head of an Extratropical Cyclone Using a 40-Member Ensemble

Abstract This paper investigates the ability of the Weather Research and Forecasting (WRF) Model in simulating multiple small-scale precipitation bands (multibands) within the extratropical cyclone comma head using four winter storm cases from 2014 to 2017. Using the model output, some physical processes are explored to investigate band prediction. A 40-member WRF ensemble was constructed down to 2-km grid spacing over the Northeast United States using different physics, stochastic physics perturbations, different initial/boundary conditions from the first five perturbed members of the Global Forecast System (GFS) Ensemble Reforecast (GEFSR), and a stochastic kinetic energy backscatter scheme (SKEBS). It was found that 2-km grid spacing is adequate to resolve most snowbands. A feature-based verification is applied to hourly WRF reflectivity fields from each ensemble member and the WSR-88D radar reflectivity at 2-km height above sea level. The Method for Object-Based Diagnostic Evaluation (MODE) tool is used for identifying multibands, which are defined as two or more bands that are 5–20 km in width and that also exhibit a >2:1 aspect ratio. The WRF underpredicts the number of multibands and has a slight eastward position bias. There is no significant difference in frontogenetical forcing, vertical stability, moisture, and vertical shear between the banded versus nonbanded members. Underpredicted band members tend to have slightly stronger frontogenesis than observed, which may be consolidating the bands, but overall there is no clear linkage in ambient condition errors and band errors, thus leaving the source for the band underprediction motivation for future work.

Simulating Climate Change Impacts on Surface Water Resources Within a Lake‐Affected Region Using Regional Climate Projections

AbstractThis study aims to assess the impact of climate change on water resources in a large watershed within the Laurentian Great Lakes region, using the fully integrated surface‐subsurface model HydroGeoSphere. The hydrologic model is forced with an ensemble of high‐resolution climate projections from the Weather Research and Forecasting (WRF) model. The latter has been extended with an interactive lake model (FLake) to capture the effect of the Great Lakes on the regional climate. The WRF ensemble encompasses two different moist physics configurations at resolutions of 90, 30, and 10 km, as well as four different initial and boundary conditions, so as to control for natural climate variability. The integrated hydrologic model is run with a representative seasonal cycle, which effectively controls natural climate variability, while remaining computationally tractable with a large integrated model. However, the range of natural variability is also investigated, as are the impacts of climate model resolution and bias correction. The two WRF configurations show opposite climate change responses in summer precipitation, but similar responses otherwise. The hydrologic simulations generally follow the climate forcing; however, due to the memory of the subsurface, the differences in summer propagate throughout the entire seasonal cycle. This results in a set of dry scenarios with reduced streamflow and water availability year‐round and a set of wet scenarios with increased streamflow for all times excluding the spring peak, which does not increase. Most of the analysis focuses on streamflow, but changes in the seasonal cycle of baseflow and groundwater recharge are also analyzed.

How many RCM ensemble members provide confidence in the impact of land‐use land cover change?

ABSTRACTRegional climate models (RCMs) include both terrestrial and atmospheric compartments and thereby allow studying land–atmosphere feedback, in particular, the impact of land‐use land cover driven by biogeophysical processes on regional climate. In this study, a method is developed to separate the signals from the noise in RCM simulations of the effects of changes in land use, using perturbed initial boundary conditions (PICs). We want to know how many ensemble members are required to identify robust and statistically significant land‐use land cover change (LULCC) effects from RCM LULCC studies. The method is applied to a case study of urbanization and deforestation, for which LULCC scenarios are implemented in the RCM Weather Research and Forecasting (WRF). Based on WRF ensemble simulations with PICs for 2010, the signal‐to‐noise ratio (SNR) is used to identify areas with pronounced effect of an LULCC or, rather, the parametrization of the land‐use classes. While in the urbanization scenarios clear and statistically significant signals are found for air temperature and for both latent‐ and sensible heat (SNR values up to 24), the effects are less pronounced for precipitation, and for deforestation in general (SNR values < 1). For the case study of urbanization and precipitation, the impact of the ensemble size is studied in order to derive robust conclusions about the effects of LULCC on precipitation. We conclude that single RCM realizations of different land‐use representations are not sufficient to derive LULCC‐induced signals, particularly not for precipitation. Small ensemble sizes led to concluding there were significant LULCC‐induced precipitation signals, but these disappeared when the ensemble size was increased. Our regional analysis suggests the need for ensemble sizes well above 10 for precipitation.

An improved multi-step forecasting model based on WRF ensembles and creative fuzzy systems for wind speed

Accurate wind speed forecasting, which strongly influences the safe usage of wind resources, is still a critical issue and a huge challenge. At present, the single-valued deterministic NWP forecast is primarily adopted by wind farms; however, recent techniques cannot meet the actual needs of grid dispatch in many cases. This paper contributes to a new multi-step forecasting method for operational wind forecast, 96-steps of the next day, termed the CS-FS-WRF-E model, which is based on a Weather Research and Forecasting (WRF) ensemble forecast, a novel Fuzzy System, and a Cuckoo Search (CS) algorithm. First, the WRF ensemble, which considers three horizontal resolutions and four initial fields, using a 0.5° horizontal grid-spacing Global Forecast System (GFS) model output, is constructed as the basic forecasting results. Then, a novel fuzzy system, which can extract the features of these ensembles, is built under the concept of membership degrees. With the help of CS optimization, the final model is constructed using this evolutionary algorithm to adjust and correct the results obtained based on physical laws, yielding the best forecasting performance and outperforming individual ensemble members and all of the other models for comparison.

An observation-constrained multi-physics WRF ensemble for simulating European mega heat waves

Abstract. Many climate models have difficulties in properly reproducing climate extremes, such as heat wave conditions. Here we use the Weather Research and Forecasting (WRF) regional climate model with a large combination of different atmospheric physics schemes, in combination with the NOAH land-surface scheme, with the goal of detecting the most sensitive physics and identifying those that appear most suitable for simulating the heat wave events of 2003 in western Europe and 2010 in Russia. In total, 55 out of 216 simulations combining different atmospheric physical schemes have a temperature bias smaller than 1 °C during the heat wave episodes, the majority of simulations showing a cold bias of on average 2–3 °C. Conversely, precipitation is mostly overestimated prior to heat waves, and shortwave radiation is slightly overestimated. Convection is found to be the most sensitive atmospheric physical process impacting simulated heat wave temperature across four different convection schemes in the simulation ensemble. Based on these comparisons, we design a reduced ensemble of five well performing and diverse scheme configurations, which may be used in the future to perform heat wave analysis and to investigate the impact of climate change during summer in Europe.

Extreme precipitation in WRF during the Newcastle East Coast Low of 2007

In the context of regional downscaling, we study the representation of extreme precipitation in the Weather Research and Forecasting (WRF) model, focusing on a major event that occurred on the 8th of June 2007 along the coast of eastern Australia (abbreviated “Newy”). This was one of the strongest extra-tropical low-pressure systems off eastern Australia in the last 30 years and was one of several storms comprising a test bed for the WRF ensemble that underpins the regional climate change projections for eastern Australia (New South Wales/Australian Capital Territory Regional Climate Modelling Project, NARCliM). Newy provides an informative case study for examining precipitation extremes as simulated by WRF set up for regional downscaling. Here, simulations from the NARCliM physics ensemble of Newy available at ∼10 km grid spacing are used. Extremes and spatio-temporal characteristics are examined using land-based daily and hourly precipitation totals, with a particular focus on hourly accumulations. Of the different physics schemes assessed, the cumulus and the boundary layer schemes cause the largest differences. Although the Betts-Miller-Janjic cumulus scheme produces better rainfall totals over the entire storm, the Kain-Fritsch cumulus scheme promotes higher and more realistic hourly extreme precipitation totals. Analysis indicates the Kain-Fritsch runs are correlated with larger resolved grid-scale vertical moisture fluxes, which are produced through the influence of parameterized convection on the larger-scale circulation and the subsequent convergence and ascent of moisture. Results show that WRF qualitatively reproduces spatial precipitation patterns during the storm, albeit with some errors in timing. This case study indicates that whilst regional climate simulations of an extreme event such as Newy in WRF may be well represented at daily scales irrespective of the physics scheme used, the representation at hourly scales is likely to be physics scheme dependent.

Regional climate hindcast simulations within EURO-CORDEX: evaluation of a WRF multi-physics ensemble

Abstract. In the current work we present six hindcast WRF (Weather Research and Forecasting model) simulations for the EURO-CORDEX (European Coordinated Regional Climate Downscaling Experiment) domain with different configurations in microphysics, convection and radiation for the time period 1990–2008. All regional model simulations are forced by the ERA-Interim reanalysis and have the same spatial resolution (0.44°). These simulations are evaluated for surface temperature, precipitation, short- and longwave downward radiation at the surface and total cloud cover. The analysis of the WRF ensemble indicates systematic temperature and precipitation biases, which are linked to different physical mechanisms in the summer and winter seasons. Overestimation of total cloud cover and underestimation of downward shortwave radiation at the surface, mostly linked to the Grell–Devenyi convection and CAM (Community Atmosphere Model) radiation schemes, intensifies the negative bias in summer temperatures over northern Europe (max −2.5 °C). Conversely, a strong positive bias in downward shortwave radiation in summer over central (40–60%) and southern Europe mitigates the systematic cold bias over these regions, signifying a typical case of error compensation. Maximum winter cold biases are over northeastern Europe (−2.8 °C); this location suggests that land–atmosphere rather than cloud–radiation interactions are to blame. Precipitation is overestimated in summer by all model configurations, especially the higher quantiles which are associated with summertime deep cumulus convection. The largest precipitation biases are produced by the Kain–Fritsch convection scheme over the Mediterranean. Precipitation biases in winter are lower than those for summer in all model configurations (15–30%). The results of this study indicate the importance of evaluating not only the basic climatic parameters of interest for climate change applications (temperature and precipitation), but also other components of the energy and water cycle, in order to identify the sources of systematic biases, possible compensatory or masking mechanisms and suggest pathways for model improvement.

Improving the spin-up of regional EnKF for typhoon assimilation and forecasting with Typhoon Sinlaku (2008)

The Running-In-Place (RIP) method is implemented in the framework of the Local Ensemble Transform Kalman Filter (LETKF) coupled with the Weather Research and Forecasting (WRF) model. RIP aims at accelerating the spin-up of the regional LETKF system when the WRF ensemble is initialised from a global analysis, which is obtained at a coarser resolution and lacks features related to the underlying mesoscale evolution. The RIP method is further proposed as an outer-loop scheme to improve the nonlinear evolution of the ensemble when the characteristics of the error statistics change rapidly owing to strong nonlinear dynamics. The impact of using RIP as an outer-loop for the WRF-LETKF system is evaluated for typhoon assimilation and prediction with Typhoon Sinlaku (2008) as a case study. For forecasts beyond one day, the typhoon track prediction is significantly improved after RIP is applied, especially during the spin-up period of the LETKF assimilation when Sinlaku is developing rapidly from a severe tropical storm to a typhoon. The impact of the dropsondes is significantly increased by RIP at early assimilation cycles. Results suggest that these improvements are because of the positive impact on the environmental condition of the typhoon. Results also suggest that using the RIP scheme adaptively allows RIP to be used as an outer-loop for the WRF-LETKF with further improvements.

Application of Object-Based Verification Techniques to Ensemble Precipitation Forecasts

Abstract Both the Method for Object-based Diagnostic Evaluation (MODE) and contiguous rain area (CRA) object-based verification techniques have been used to analyze precipitation forecasts from two sets of ensembles to determine if spread-skill behavior observed using traditional measures can be seen in the object parameters. One set consisted of two eight-member Weather Research and Forecasting (WRF) model ensembles: one having mixed physics and dynamics with unperturbed initial and lateral boundary conditions (Phys) and another using common physics and a dynamic core but with perturbed initial and lateral boundary conditions (IC/LBC). Traditional measures found that spread grows much faster in IC/LBC than in Phys so that after roughly 24 h, better skill and spread are found in IC/LBC. These measures also reflected a strong diurnal signal of precipitation. The other set of ensembles included five members of a 4-km grid-spacing WRF ensemble (ENS4) and five members of a 20-km WRF ensemble (ENS20). Traditional measures suggested that the diurnal signal was better in ENS4 and spread increased more rapidly than in ENS20. Standard deviations (SDs) of four object parameters computed for the first set of ensembles using MODE and CRA showed the trend of enhanced spread growth in IC/LBC compared to Phys that had been observed in traditional measures, with the areal coverage of precipitation exhibiting the greatest growth in spread with time. The two techniques did not produce identical results; although, they did show the same general trends. A diurnal signal could be seen in the SDs of all parameters, especially rain rate, volume, and areal coverage. MODE results also found evidence of a diurnal signal and faster growth of spread in object parameters in ENS4 than in ENS20. Some forecasting approaches based on MODE and CRA output are also demonstrated. Forecasts based on averages of object parameters from each ensemble member were more skillful than forecasts based on MODE or CRA applied to an ensemble mean computed using the probability matching technique for areal coverage and volume, but differences in the two techniques were less pronounced for rain rate and displacement. The use of a probability threshold to define objects was also shown to be a valid forecasting approach with MODE.

Initial Condition Sensitivity of Western Pacific Extratropical Transitions Determined Using Ensemble-Based Sensitivity Analysis

Abstract An ensemble Kalman filter based on the Weather Research and Forecasting (WRF) model is used to generate ensemble analyses and forecasts for the extratropical transition (ET) events associated with Typhoons Tokage (2004) and Nabi (2005). Ensemble sensitivity analysis is then used to evaluate the relationship between forecast errors and initial condition errors at the onset of transition, and to objectively determine the observations having the largest impact on forecasts of these storms. Observations from rawinsondes, surface stations, aircraft, cloud winds, and cyclone best-track position are assimilated every 6 h for a period before, during, and after transition. Ensemble forecasts initialized at the onset of transition exhibit skill similar to the operational Global Forecast System (GFS) forecast and to a WRF forecast initialized from the GFS analysis. WRF ensemble forecasts of Tokage (Nabi) are characterized by relatively large (small) ensemble variance and greater (smaller) sensitivity to the initial conditions. In both cases, the 48-h forecast of cyclone minimum SLP and the RMS forecast error in SLP are most sensitive to the tropical cyclone position and to midlatitude troughs that interact with the tropical cyclone during ET. Diagnostic perturbations added to the initial conditions based on ensemble sensitivity reduce the error in the storm minimum SLP forecast by 50%. Observation impact calculations indicate that assimilating approximately 40 observations in regions of greatest initial condition sensitivity produces a large, statistically significant impact on the 48-h cyclone minimum SLP forecast. For the Tokage forecast, assimilating the single highest impact observation, an upper-tropospheric zonal wind observation from a Mongolian rawinsonde, yields 48-h forecast perturbations in excess of 10 hPa and 60 m in SLP and 500-hPa height, respectively.