High-Resolution Rapid Refresh Model Research Articles

Forecasting High Wind Events in the HRRR Model over Wyoming and Colorado. Part I: Evaluation of Wind Speeds and Gusts

Abstract Strong wind events cause significant societal damage ranging from loss of property and disruption of commerce to loss of life. Over portions of the United States, the strongest winds occur in the cold season and may be driven by interactions with the terrain (downslope winds, gap flow, and mountain wave activity). In the first part of this two-part series, we evaluate the High-Resolution Rapid Refresh (HRRR) model wind speed and gust forecasts for the 2016–22 winter months over Wyoming and Colorado, an area prone to downslope windstorms and gap flows due to its complex topography. The HRRR model exhibits a positive bias for low wind speeds/gusts and a large negative bias for strong wind speeds/gusts. In general, the model misses many strong wind events, but when it does predict strong winds, there is a high false alarm probability. An analysis of proxies for surface winds is conducted. Specifically, 700- and 850-mb (1 mb = 1 hPa) geopotential height gradients are found to be good proxies for strong wind speeds and gusts at two wind-prone locations in Wyoming. Given the good agreement between low-level height gradients and surface wind speeds yet a strong negative bias for strong wind speeds and gusts, there is a potential shortcoming in the boundary layer physics in the HRRR model with regard to predicting strong winds over complex terrain, which is the focus of the second part of this two-part study. Last, the sites with the largest strong wind speed bias are found to mostly sit on the leeward side of high mountains, suggesting that the HRRR model performs poorly in the prediction of downslope windstorms. Significance Statement We investigate the performance of the High-Resolution Rapid Refresh (HRRR) model with respect to strong wintertime wind speeds and gusts over the complex terrain of Wyoming and Colorado. We show that the overall performance of the HRRR model is low with regard to strong wind speed and wind gust forecasts across the investigated winter seasons, with a large negative bias in predicted strong wind speeds and gusts and a small positive bias for weak wind speeds and gusts. The largest biases are found to be on the leeward side of high mountains, indicating poor prediction of downslope winds. This study also utilizes National Weather Service forecasting metrics to understand their performance with respect to strong wind forecasts, and we find that they provide skill in forecasting these events.

Forecasting High Wind Events in the HRRR Model over Wyoming and Colorado. Part II: Sensitivity of Surface Wind Speeds to Model Resolution and Physics

Abstract Strong wind events can cause severe economic loss and societal impacts. Peak winds over Wyoming and Colorado occur during the wintertime months, and in the first part of this two-part series, it was shown that the High-Resolution Rapid Refresh (HRRR) model displays large negative biases with respect to strong wind events. In this part of the study, we address two questions: 1) does increasing the horizontal resolution improve the representation of strong wind events over this region, and 2) are the biases in HRRR-forecasted winds related to the selected planetary boundary layer (PBL), surface layer (SL), and/or land surface model (LSM) parameterizations? We conduct Weather Research and Forecasting (WRF) Model simulations to address these two main questions. Increasing the horizontal resolution leads to a small improvement in the simulation of strong wind speeds over the complex terrain of Wyoming and Colorado. In general, changes in the PBL, SL, and LSM parameterizations show much larger changes in simulated wind speeds compared to increasing the model resolution alone. Specifically, changing from the Mellor–Yamada–Nakanishi–Niino scheme to the Mellor–Yamada–Janjić PBL and SL schemes results in nearly no change in the r2 values, but there is a decrease in the magnitude of the strong wind speed bias (from −12.52 to −10.16 m s−1). We attribute these differences to differences in the diagnosis of surface wind speeds and mixing in the boundary layer. Further analysis is conducted to determine the value of 1-km forecasts of strong winds compared with wind speed diagnostics commonly used by the National Weather Service. Significance Statement Motivated by prior studies showing low skill in the prediction of strong winds using state-of-the-art weather forecast models, in this study, we aim to investigate two questions: 1) does increasing the horizontal resolution improve the prediction of strong wind events over the complex terrain of Wyoming and Colorado, and 2) are the biases in High-Resolution Rapid Refresh (HRRR) forecasted winds related to the selected planetary boundary layer, surface layer, and/or land surface model parameterizations? We find that increasing the horizontal resolution provides a slight improvement in the prediction of strong winds. Further, considerable improvement in the prediction of strong winds is found for varying boundary layer, surface layer, and land surface parameterizations.

Analysis of bias correction of HRRR model outputs for offshore wind power ramp events

This study examines the impact of bias correction on the prediction of offshore wind power ramp events (WPREs) using NOAA's High-Resolution Rapid Refresh (HRRR) model. The analysis focuses on New England's offshore region, using Light Detection and Ranging (LiDAR) wind observations and the power curve derived from the nearby Block Island wind farm as a reference. The study characterizes bias in terms of amplitude and phase, finding a 30-min systematic lag error and a marginal underestimation of high wind speeds. Perfect bias correction treatments including both the amplitude and the phase are applied to HRRR wind speed, and their impact on WPREs is assessed using categorical performance measures. A large discrepancy in the number of cases with HRRR having much fewer cases than the reference is found, which is attributed to the model's temporal and spatial resolution. The results show that the amplitude bias correction has a marginal impact on reducing frequency bias, while the dynamic time warping (DTW) used in the phase bias correction shows promising results but limited performance. The study concludes that phase bias correction can be a feasible option to reduce bias in WPREs in offshore power systems.

An Evaluation of NOAA Modeled and In Situ Soil Moisture Values and Variability across the Continental United States

Abstract Estimates of soil moisture from two National Oceanic and Atmospheric Administration (NOAA) models are compared to in situ observations. The estimates are from a high-resolution atmospheric model with a land surface model [High-Resolution Rapid Refresh (HRRR) model] and a hydrologic model from the NOAA Climate Prediction Center (CPC). Both models produce wetter soils in dry regions and drier soils in wet regions, as compared to the in situ observations. These soil moisture differences occur at most soil depths but are larger at the deeper depths below the surface (100 cm). Comparisons of soil moisture variability are also assessed as a function of soil moisture regime. Both models have lower standard deviations as compared to the in situ observations for all soil moisture regimes. The HRRR model’s soil moisture is better correlated with in situ observations for drier soils as compared to wetter soils—a trend that was not present in the CPC model comparisons. In terms of seasonality, soil moisture comparisons vary depending on the metric, time of year, and soil moisture regime. Therefore, consideration of both the seasonality and soil moisture regime is needed to accurately determine model biases. These NOAA soil moisture estimates are used for a variety of forecasting and societal applications, and understanding their differences provides important context for their applications and can lead to model improvements. Significance Statement Soil moisture is an essential variable coupling the land surface to the atmosphere. Accurate estimates of soil moisture are important for forecasting near-surface temperature and moisture, predicting where clouds will form, and assessing drought and fire risks. There are multiple estimates of soil moisture available, and in this study, we compare soil moisture estimates from two different National Oceanic and Atmospheric Administration (NOAA) models to in situ observations. These comparisons include both soil moisture amount and variability and are conducted at several soil depths, in different soil moisture regimes, and for different seasons and years. This comprehensive assessment allows for an accurate assessment of biases within these models that would be missed when conducting analyses more broadly.

Verification of the Global Forecast System, North American Mesoscale Forecast System, and High-Resolution Rapid Refresh Model Near-Surface Forecasts by Use of the New York State Mesonet

Abstract Global Forecast System (GFS), North American Mesoscale Forecast System (NAM), and High-Resolution Rapid Refresh (HRRR) 2-m temperature, 10-m wind speed, and precipitation accumulation forecasts initialized at 1200 UTC are verified against New York State Mesonet (NYSM) observations from 1 January 2018 through 31 December 2021. NYSM observations at 126 site locations are used to calculate standard error statistics (e.g., forecast error, root-mean-square error) for temperature and wind speed and contingency table statistics for precipitation across forecast hours, meteorological seasons, and regions. The majority of the focus is placed on the first 18 forecast hours to allow for comparison among all three models. A daily NYSM station-mean temperature error analysis identified a slight cold bias at temperatures below 25°C in the GFS, a cool-to-warm bias as forecast temperatures warm in the HRRR, and a warm bias at temperatures above 30°C in each model. Differences arise when considering temperature biases with respect to lead times and seasons. Wind speeds are overforecast at all ranges in each season, and forecast wind speeds ≥ 18 m s−1 are rarely observed. Performance diagrams indicate overall good forecast performance at precipitation thresholds of 0.1–1.5 mm, but with a high frequency bias in the GFS and NAM. This paper provides an overview of deterministic forecast performance across New York State, with the aim of sharing common biases associated with temperature, wind speed, and precipitation with operational forecasters and is the first step in developing a real-time model forecast uncertainty prediction tool.

Methods for Validating HRRR Simulated Cloud Properties for Different Weather Phenomena Using Satellite and Radar Observations

Abstract In this study, we evaluate the ability of the High-Resolution Rapid Refresh (HRRR) model to forecast cloud characteristics through comparison of observed and simulated satellite brightness temperatures (BTs) and radar reflectivity during different weather phenomena in December 2021: the Mayfield, Kentucky, tornado on 11 December, a heavy snow event in Minnesota from 10 to 11 December, and the Midwest derecho on 15 December. This is done to illustrate the importance of examining model accuracy across a range of weather phenomena. Observation and forecast objects were created using the Method for Object-Based Diagnostic Evaluation (MODE). HRRR accurately depicted the spatial displacements between observation cloud (defined using BTs) and radar reflectivity objects, namely, the centers of cloud objects are to the east of the radar objects for the tornado and derecho events, and generally west of the radar objects for the snow event. However, HRRR had higher (less intense) simulated BTs and higher (more intense) radar reflectivity than the observations for the tornado event. Simulated radar reflectivity is higher and BTs are lower than the observations during the middle of the snow event. Also, simulated radar reflectivity is higher and BTs are lower than the observations during the derecho event. Of the three weather events, the HRRR forecasts are most accurate for the snow event, based on the object-based threat score, followed by the derecho and tornado events. The tornado event has lower accuracy because matches between paired simulated and observation objects are worse than for the snow event, with less similarity in size forecast objects and greater distance between paired object centers. Significance Statement The purpose of this study is to assess the accuracy of forecast cloud and radar objects, defined using simulated satellite brightness temperatures and radar reflectivity, from the High-Resolution Rapid Refresh (HRRR) model. This assessment was conducted for a tornado, snow, and derecho event from December 2021. Results from these three events indicate that the HRRR model accurately represents the observed displacement between the center of cloud and radar objects for the tornado and derecho events, and is the most accurate overall for the snow event.

Case study of a bore wind-ramp event from lidar measurements and HRRR simulations over ARM Southern Great Plains

The rapid change of wind speed and direction on 21 August 2017 is studied using Doppler lidar measurements at five sites of the Atmospheric Radiation Measurement (ARM) Southern Great Plains (SGP) facility in north-central Oklahoma. The Doppler lidar data were investigated along with meteorological variables such as temperature, humidity, and turbulence available from the large suite of instrumentation deployed at the SGP Central Facility (C1) during the Land-Atmosphere Feedback Experiment in August 2017. Lidar measurements at five sites, separated by 55–70 km, allowed us to document the development and evolution of the wind flow over the SGP area, examine synoptic conditions to understand the mechanism that leads to the ramp event, and estimate the ability of the High-Resolution Rapid Refresh model to reproduce this event. The flow feature in question is an atmospheric bore, a small-scale phenomenon that is challenging to represent in models, that was generated by a thunderstorm outflow northwest of the ARM SGP area. The small-scale nature of bores, its impact on power generation, and the modeling challenges associated with representing bores are discussed in this paper. The results also provide information about model errors between sites of different surface and vegetation types.

Precipitation Vertical Structure Characterization: A Feature-Based Approach

Abstract The three-dimensional (3D) structure of precipitation systems is highly dependent on hydrometeor formation processes and microphysics. This study aims to characterize distinct vertical profiles of precipitation regimes by relying on the availability of a high-quality, spatially dense radar network and its capability to observe the 3D structure of the storms. A deep-learning-based framework, coupled with unsupervised clustering methods, is developed to identify types of precipitation structures irrespective of their physical properties. A 6-month period of 3D reflectivity profiles from the Multi-Radar Multi-Sensor (MRMS) network is used to identify different regimes and investigate their properties with respect to the underlying environmental conditions. Dominant features retrieved from radar reflectivity profiles using convolutional neural-network-based autoencoders are employed to identify similar-looking vertical structures using coupled k-means and agglomerative clustering algorithms. The k-means method identifies distinct groups, while the agglomerative clustering visualizes intercluster relationships. The framework identifies 18 clusters that can be broadly combined into five groups of varied echo-top heights. The 18 clusters demonstrate variability with respect to structural features and precipitation rate/type, implying that profiles in each group belong to a physically different precipitation regime. An independent analysis of the regime properties is conducted by matching the MRMS reflectivity profiles with environmental parameters derived from the High-Resolution Rapid Refresh model forecasts. The distribution of the environmental variables confirms cluster-specific feature properties, confirming the physics-based regime separation across the clusters and their dependence on the vertical structure. The identified precipitation regimes can assist in developing physics-guided retrievals and studying precipitation regimes. Significance Statement This study proposes a systematic model to identify precipitation profiles of distinct vertical structures and evaluate their dependence on environmental conditions. The model was developed using ground-based radar observations; however, there is potential to extend this model to reflectivity profiles from both ground- and satellite-based sensors. In addition, the identified precipitation regime clusters could be a proxy for the vertical structure of precipitation systems and assist in determining the structural variability within traditional precipitation type classification (e.g., convective versus stratiform). Moreover, identifying the precipitation regimes could also be used to improve satellite-based precipitation retrievals. Finally, a better understanding of precipitation structure would also help improve the initialization of climate models.

A Kinematic Modeling Study of the Reorganization of Snowfall between Cloud-Top Generating Cells and Low-Level Snowbands in Midlatitude Winter Storms

Abstract This paper explores whether particles within uniformly spaced generating cells falling at terminal velocity within observed 2D wind fields and idealized deformation flow beneath cloud top can be reorganized consistent with the presence of single and multibanded structures present on WSR-88D radars. In the first experiment, two-dimensional wind fields, calculated along cross sections normal to the long axis of snowbands observed during three northeast U.S. winter storms, were taken from the initialization of the High-Resolution Rapid Refresh model. This experiment demonstrated that the greater the residence time of the particles in each of the three storms, the greater particle reorganization occurred. For experiments with longer residence times, increases in particle concentrations were nearly or directly collocated with reflectivity bands. For experiments with shorter residence times, particle reorganization still conformed to the band features but with less concentration enhancement. This experiment demonstrates that the combination of long particle residence time and net convergent cross-sectional flow through the cloud depth is sufficient to reorganize particles into locations consistent with precipitation bands. Increased concentrations of ice particles can then contribute, along with any dynamic forcing, to the low-level reflectivity bands seen on WSR-88D radars. In a second experiment, the impact of flow deformation on the reorganization of falling ice particles was investigated using an idealized kinematic model with stretching deformation flow of different depths and magnitudes. These experiments showed that deformation flow provides for little particle reorganization given typical deformation layer depths and magnitudes within the comma head of such storms. Significance Statement Past research with vertically pointing and scanning radars presents two different perspectives regarding snowfall organization in winter storms. Vertically pointing radars often observe cloud-top generating cells with precipitation fallstreaks descending into a broad stratiform echo at lower altitudes. In contrast, scanning radars often observe snowfall organized in quasi-linear bands. This work attempts to provide a connection between these two perspectives by examining how two-dimensional convergent and deformation flow occurring in winter storms can contribute to the reorganization of snowfall between cloud top and the ground.

Competing Forecast Verification: Using the Power-Divergence Statistic for Testing the Frequency of “Better”

Abstract When testing hypotheses about which of two competing models is better, say A and B, the difference is often not significant. An alternative, complementary approach, is to measure how often model A is better than model B regardless of how slight or large the difference. The hypothesis concerns whether or not the percentage of time that model A is better than model B is larger than 50%. One generalized test statistic that can be used is the power-divergence test, which encompasses many familiar goodness-of-fit test statistics, such as the loglikelihood-ratio and Pearson X2 tests. Theoretical results justify using the distribution for the entire family of test statistics, where k is the number of categories. However, these results assume that the underlying data are independent and identically distributed, which is often violated. Empirical results demonstrate that the reduction to two categories (i.e., model A is better than model B versus model B is better than A) results in a test that is reasonably robust to even severe departures from temporal independence, as well as contemporaneous correlation. The test is demonstrated on two different example verification sets: 6-h forecasts of eddy dissipation rate (m2/3 s−1) from two versions of the Graphical Turbulence Guidance model and for 12-h forecasts of 2-m temperature (°C) and 10-m wind speed (m s−1) from two versions of the High-Resolution Rapid Refresh model. The novelty of this paper is in demonstrating the utility of the power-divergence statistic in the face of temporally dependent data, as well as the emphasis on testing for the “frequency-of-better” alongside more traditional measures.

Verification of Rapid Refresh and High-Resolution Rapid Refresh Model Variables in Tornadic Tropical Cyclones

Abstract Tropical cyclone tornadoes (TCTORs) are a hazard to life and property during landfalling tropical cyclones (TCs). The threat is often spread over a wide area within the TC envelope and must be continually evaluated as the TC moves inland and dissipates. To anticipate the risk of TCTORs, forecasters may use high-resolution, rapidly updating model analyses and short-range forecasts such as the Rapid Refresh (RAP) and High-Resolution Rapid Refresh (HRRR), and an ingredients-based approach similar to that used for forecasting continental midlatitude tornadoes. Though RAP and HRRR errors have been identified in typical midlatitude convective environments, this study evaluates the performance of the RAP and the HRRR within the TC envelope, with particular attention given to sounding-derived parameters previously identified as useful for TCTOR forecasting. A sample of 1730 observed upper-air soundings is sourced from 13 TCs that made landfall along the U.S. coastline between 2017 and 2019. The observed soundings are paired with their corresponding model gridpoint soundings from the RAP analysis, RAP 12-h forecast, and HRRR 12-h forecast. Model errors are calculated for both the raw sounding variables of temperature, dewpoint, and wind speed, as well as for the selected sounding-derived parameters. Results show a moist bias that worsens with height across all model runs. There are also statistically significant underpredictions in stability-related parameters such as convective available potential energy (CAPE) and kinematic parameters such as vertical wind shear.

Evaluation of a cloudy cold-air pool in the Columbia River basin in different versions of the High-Resolution Rapid Refresh (HRRR) model

Abstract. The accurate forecast of persistent orographic cold-air pools in numerical weather prediction models is essential for the optimal integration of wind energy into the electrical grid during these events. Model development efforts during the second Wind Forecast Improvement Project (WFIP2) aimed to address the challenges related to this. We evaluated three versions of the National Oceanic and Atmospheric Administration (NOAA) High-Resolution Rapid Refresh model with two different horizontal grid spacings against in situ and remote sensing observations to investigate how developments in physical parameterizations and numerical methods targeted during WFIP2 impacted the simulation of a persistent cold-air pool in the Columbia River basin. Differences amongst model versions were most apparent in simulated temperature and low-level cloud fields during the persistent phase of the cold-air pool. The model developments led to an enhanced low-level cloud cover, resulting in better agreement with the observations. This removed a diurnal cycle in the near-surface temperature bias at stations throughout the basin by reducing a cold bias during the night and a warm bias during the day. However, low-level clouds did not clear sufficiently during daytime in the newest model version, which leaves room for further model developments. The model developments also led to a better representation of the decay of the cold-air pool by slowing down its erosion.

Application of geostationary satellite and high-resolution meteorology data in estimating hourly PM2.5 levels during the Camp Fire episode in California

BACKGROUND: Wildland fire smoke contains large amounts of PM2.5 that can traverse tens to hundreds of kilometers, resulting in significant deterioration of air quality and excess mortality and morbidity in downwind regions. Estimating PM2.5 levels while considering the impact of wildfire smoke has been challenging due to the lack of ground monitoring coverage near the smoke plumes. AIM: We aim to estimate total PM2.5 concentration during the Camp Fire episode, the deadliest wildland fire in California history. METHODS: Our random forest (RF) model combines calibrated low-cost sensor data (PurpleAir) with regulatory monitor measurements (Air Quality System, AQS) to bolster ground observations, Geostationary Operational Environmental Satellite-16 (GOES-16)’s high temporal resolution to achieve hourly predictions, and oversampling techniques (Synthetic Minority Oversampling Technique, SMOTE) to reduce model underestimation at high PM2.5 levels. In addition, meteorological fields at 3 km resolution from the High-Resolution Rapid Refresh model and land use variables were also included in the model. RESULTS: Our AQS-only model achieved an out of bag (OOB) R2 (RMSE) of 0.84 (12.00 μg/m3) and spatial and temporal cross-validation (CV) R2 (RMSE) of 0.74 (16.28 μg/m3) and 0.73 (16.58 μg/m3), respectively. Our AQS + Weighted PurpleAir Model achieved OOB R2 (RMSE) of 0.86 (9.52 μg/m3) and spatial and temporal CV R2 (RMSE) of 0.75 (14.93 μg/m3) and 0.79 (11.89 μg/m3), respectively. Our AQS + Weighted PurpleAir + SMOTE Model achieved OOB R2 (RMSE) of 0.92 (10.44 μg/m3) and spatial and temporal CV R2 (RMSE) of 0.84 (12.36 μg/m3) and 0.85 (14.88 μg/m3), respectively. CONCLUSIONS: Hourly predictions from our model may aid in epidemiological investigations of intense and acute exposure to PM2.5 during the Camp Fire episode. KEYWORDS: PM2.5, GOES16, Wildland fire, Remote sensing, AOD, SMOTE, Weighted Random Forest

Radar Reflectivity–Based Model Initialization Using Specified Latent Heating (Radar-LHI) within a Diabatic Digital Filter or Pre-Forecast Integration

Abstract A technique for model initialization using three-dimensional radar reflectivity data has been developed and applied within the NOAA 13-km Rapid Refresh (RAP) and 3-km High-Resolution Rapid Refresh (HRRR) regional forecast systems. This technique enabled the first assimilation of radar reflectivity data for operational NOAA forecast models, critical especially for more accurate short-range prediction of convective storms. For the RAP, the technique uses a diabatic digital filter initialization (DFI) procedure originally deployed to control initial inertial gravity wave noise. Within the forward-model integration portion of diabatic DFI, temperature tendencies obtained from the model cloud/precipitation processes are replaced by specified latent heating–based temperature tendencies derived from the three-dimensional radar reflectivity data, where available. To further refine initial conditions for the convection-allowing HRRR model, a similar procedure is used in the HRRR, but without DFI. Both of these procedures, together called the “Radar-LHI” (latent heating initialization) technique, have been essential for initialization of ongoing precipitation systems, especially convective systems, within all NOAA operational versions of the 13-km RAP and 3-km HRRR models extending through the latest implementation upgrade at NCEP in 2020. Application of the latent heat–derived temperature tendency induces a vertical circulation with low-level convergence and upper-level divergence in precipitation systems. Retrospective tests of the Radar-LHI technique show significant improvement in short-range (0–6 h) precipitation system forecasts, as revealed by reflectivity verification scores. Results presented document the impact on HRRR reflectivity forecasts of the radar reflectivity initialization technique applied to the RAP alone, HRRR alone, and both the RAP and HRRR. Significance Statement The large forecast uncertainty of convective situations, even at short lead times, coupled with the hazardous weather they produce, makes convective storm prediction one of the most significant short-range forecast challenges confronting the operational numerical weather prediction community. Prediction of heavy precipitation events also requires accurate initialization of precipitation systems. An innovative assimilation technique using radar reflectivity data to initialize NOAA operational weather prediction models is described. This technique, which uses latent heating specified from radar reflectivity (and can accommodate lightning data and other convection/precipitation indicators), was first implemented in 2009 at NOAA/NCEP and continues to be used in 2022 in the NCEP-operational RAP and HRRR models, making it a backbone of the NOAA rapidly updated numerical weather prediction capability.

The High-Resolution Rapid Refresh (HRRR): An Hourly Updating Convection-Allowing Forecast Model. Part II: Forecast Performance

Abstract The High-Resolution Rapid Refresh (HRRR) is a convection-allowing implementation of the Advanced Research version of the Weather Research and Forecast (WRF-ARW) Model that covers the conterminous United States and Alaska and runs hourly (for CONUS; every 3 h for Alaska) in real time at the National Centers for Environmental Prediction. The high-resolution forecasts support a variety of user applications including aviation, renewable energy, and prediction of many forms of severe weather. In this second of two articles, forecast performance is documented for a wide variety of forecast variables and across HRRR versions. HRRR performance varies across geographical domain, season, and time of day depending on both prevalence of particular meteorological phenomena and the availability of both conventional and nonconventional observations. Station-based verification of surface weather forecasts (2-m temperature and dewpoint temperature, 10-m winds, visibility, and cloud ceiling) highlights the ability of the HRRR to represent daily planetary boundary layer evolution and the development of convective and stratiform cloud systems, while gridded verification of simulated composite radar reflectivity and quantitative precipitation forecasts reveals HRRR predictive skill for summer and winter precipitation systems. Significant improvements in performance for specific forecast problems are documented for the upgrade versions of the HRRR (HRRRv2, v3, and v4) implemented in 2016, 2018, and 2020, respectively. Development of the HRRR model data assimilation and physics paves the way for future progress with operational convective-scale modeling. Significance Statement NOAA’s operational hourly updating convection-allowing model, the High-Resolution Rapid Refresh (HRRR), is a key tool for short-range weather forecasting and situational awareness. Improvements in assimilation of weather observations, as well as in physics parameterizations, has led to improvements in simulated radar reflectivity and quantitative precipitation forecasts since the initial implementation of HRRR in September 2014. Other targeted development has focused on improved representation of the diurnal cycle of the planetary boundary layer, resulting in improved near-surface temperature and humidity forecasts. Additional physics and data assimilation changes have led to improved treatment of the development and erosion of low-level clouds, including subgrid-scale clouds. The final version of HRRR features storm-scale ensemble data assimilation and explicit prediction of wildfire smoke plumes.

An Evaluation of Surface Wind and Gust Forecasts from the High-Resolution Rapid Refresh Model

Abstract We utilized high temporal resolution, near-surface observations of sustained winds and gusts from two networks, the primarily airport-based Automated Surface Observing System (ASOS) and the New York State Mesonet (NYSM), to evaluate forecasts from the operational High-Resolution Rapid Refresh (HRRR) model, versions 3 and 4. Consistent with past studies, we showed the model has a high degree of skill in reproducing the diurnal variation of network-averaged wind speed of ASOS stations, but also revealed several areas where improvements could be made. Forecasts were found to be underdispersive, deficient in both temporal and spatial variability, with significant errors occurring during local nighttime hours in all regions and in forested environments for all hours of the day. This explained why the model overpredicted the network-averaged wind in the NYSM because much of that network’s stations are in forested areas. A simple gust parameterization was shown not only to have skill in predicting gusts in both networks but also to mitigate systemic biases found in the sustained wind forecasts. Significance Statement Many users depend on forecasts from operational models and need to know their strengths, weaknesses, and limitations. We examined generally high-quality near-surface observations of sustained winds and gusts from the nationwide Automated Surface Observing System (ASOS) and the New York State Mesonet (NYSM) and used them to evaluate forecasts from the previous (version 3) and current (version 4) operational High-Resolution Rapid Refresh (HRRR) model for a selected month. Evidence indicated that the wind forecasts are excellent yet imperfect and areas for further improvement remain. In particular, we showed there is a high degree of skill in representing the diurnal variation of sustained wind at ASOS stations but insufficient spatial and temporal forecast variability and overprediction at night everywhere, in forested areas at all times of day, and at NYSM sites in particular, which are more likely to be sited in the forest. Gusts are subgrid even at the fine grid spacing of the HRRR (3 km) and thus must be parameterized. Our simple gust algorithm corrected for some of these systemic biases, resulting in very good predictions of the maximum hourly gust.

Estimating the economic impacts of improved wind speed forecasts in the United States electricity sector

Each year the U.S. government makes significant investments in improving weather forecast models. In this paper, we use a multidisciplinary approach to examine how utilities can benefit from improved wind-speed forecasts to more efficiently use wind-generated electricity and subsequently increase economic activity. Specifically, we examine how improvements to the National Oceanic and Atmospheric Administration's high-resolution rapid refresh model (HRRR) wind forecasts can provide (1) cost savings for utilities and (2) increase in real household income. To do so, we compare 12-h-ahead wind forecasts with real-time observations for two HRRR model transitions (i.e., when one model is operational, the other is being tested). We compare estimates of actual and predicted wind power under the publicly available and developmental models, with reduced forecast errors allowing for better utility decision-making and lower production costs. We then translate potential cost savings into electricity price changes, which are entered as exogenous shocks to eight regional computable general equilibrium models constructed for the U.S. Overall, we find that households would have seen a potential $60 million increase in real income for our sample (13% of all contiguous U.S. land-based turbine capacity), which had the updated HRRR models been in place during the two transition periods; applying our estimated savings for the sample of turbines to the entire array of turbines shows a potential real household income increase in approximately $384 million during these time frames.

Using Cloud Computing to Analyze Model Output Archived in Zarr Format

Abstract Numerical weather prediction centers rely on the Gridded Binary Second Edition (GRIB2) file format to efficiently compress and disseminate model output as two-dimensional grids. User processing time and storage requirements are high if many GRIB2 files with size O(100 MB, where B = bytes) need to be accessed routinely. We illustrate one approach to overcome such bottlenecks by reformatting GRIB2 model output from the High-Resolution Rapid Refresh (HRRR) model of the National Centers for Environmental Prediction to a cloud-optimized storage type, Zarr. Archives of the original HRRR GRIB2 files and the resulting Zarr stores on Amazon Web Services (AWS) Simple Storage Service (S3) are available publicly through the Amazon Sustainability Data Initiative. Every hour, the HRRR model produces 18- or 48-hourly GRIB2 surface forecast files of size O(100 MB). To simplify access to the grids in the surface files, we reorganize the HRRR model output for each variable and vertical level into Zarr stores of size O(1 MB), with chunks O(10 kB) containing all forecast lead times for 150 × 150 gridpoint subdomains. Open-source libraries provide efficient access to the compressed Zarr stores using cloud or local computing resources. The HRRR-Zarr approach is illustrated for common applications of sensible weather parameters, including real-time alerts for high-impact situations and retrospective access to output from hundreds to thousands of model runs. For example, time series of surface pressure forecast grids can be accessed using AWS cloud computing resources approximately 40 times as fast from the HRRR-Zarr store as from the HRRR-GRIB2 archive. Significance Statement The rapid evolution of computing power and data storage have enabled numerical weather prediction forecasts to be generated faster and with more detail than ever before. The increased temporal and spatial resolution of forecast model output can force end users with finite memory and storage capabilities to make pragmatic decisions about which data to retrieve, archive, and process for their applications. We illustrate an approach to alleviate this access bottleneck for common weather analysis and forecasting applications by using the Amazon Web Services (AWS) Simple Storage Service (S3) to store output from the High-Resolution Rapid Refresh (HRRR) model in Zarr format. Zarr is a relatively new data storage format that is flexible, compressible, and designed to be accessed with open-source software either using cloud or local computing resources. The HRRR-Zarr dataset is publicly available as part of the AWS Sustainability Data Initiative.

Application of geostationary satellite and high-resolution meteorology data in estimating hourly PM2.5 levels during the Camp Fire episode in California

Wildland fire smoke contains large amounts of PM2.5 that can traverse tens to hundreds of kilometers, resulting in significant deterioration of air quality and excess mortality and morbidity in downwind regions. Estimating PM2.5 levels while considering the impact of wildfire smoke has been challenging due to the lack of ground monitoring coverage near the smoke plumes. We aim to estimate total PM2.5 concentration during the Camp Fire episode, the deadliest wildland fire in California history. Our random forest (RF) model combines calibrated low-cost sensor data (PurpleAir) with regulatory monitor measurements (Air Quality System, AQS) to bolster ground observations, Geostationary Operational Environmental Satellite-16 (GOES-16)’s high temporal resolution to achieve hourly predictions, and oversampling techniques (Synthetic Minority Oversampling Technique, SMOTE) to reduce model underestimation at high PM2.5 levels. In addition, meteorological fields at 3 km resolution from the High-Resolution Rapid Refresh model and land use variables were also included in the model. Our AQS-only model achieved an out of bag (OOB) R2 (RMSE) of 0.84 (12.00 μg/m3) and spatial and temporal cross-validation (CV) R2 (RMSE) of 0.74 (16.28 μg/m3) and 0.73 (16.58 μg/m3), respectively. Our AQS + Weighted PurpleAir Model achieved OOB R2 (RMSE) of 0.86 (9.52 μg/m3) and spatial and temporal CV R2 (RMSE) of 0.75 (14.93 μg/m3) and 0.79 (11.89 μg/m3), respectively. Our AQS + Weighted PurpleAir + SMOTE Model achieved OOB R2 (RMSE) of 0.92 (10.44 μg/m3) and spatial and temporal CV R2 (RMSE) of 0.84 (12.36 μg/m3) and 0.85 (14.88 μg/m3), respectively. Hourly predictions from our model may aid in epidemiological investigations of intense and acute exposure to PM2.5 during the Camp Fire episode.

Forecast Skill for California Heavy Precipitation Periods from the High-Resolution Rapid Refresh Model and the Coupled Ocean–Atmosphere Mesoscale Prediction System

Abstract Precipitation forecasts from the High-Resolution Rapid Refresh model (HRRR) of the National Centers for Environmental Prediction (NCEP) and the Navy’s Coupled Ocean–Atmosphere Mesoscale Prediction System (COAMPS) are examined during heavy precipitation periods in California. Precipitation forecast discrepancies between the two models are examined during a recent heavy winter precipitation episode in California from 6 to 8 December 2019. The skill of initial 12-h precipitation forecasts is examined objectively from 1 December 2018 to 28 February 2019 from the HRRR, COAMPS, and NCEP’s North American Mesoscale Forecast System (NAM-3km). The HRRR exhibited lower seasonal biases and higher skill based on several metrics applied to a sample of 48 12-h periods during California’s second wettest winter season during the past 20 years. Overall, the NAM-3km and COAMPS exhibited a large wet bias over the interior mountain regions while the HRRR model indicated a dry bias along the northern coastal region. All models tended to underestimate precipitation along the coastal mountains of Northern California. To highlight the regional and localized nature of forecast skill, the fraction skill score (FSS) metric is applied across ranges of spatial scales and precipitation values. For the domain as a whole, the HRRR had higher precipitation forecast skill compared to the other two models, particularly within radial distances of 20–30 km and moderate (10–50 mm) precipitation totals. FSS computed locally highlights the HRRR’s overall higher skill as well as enhanced skill in the southern half of the state.