High-Resolution Rapid Refresh Research Articles

Probabilistic Forecasts of Atmospheric River events using the HRRR Ensemble

The nine-member High-Resolution Rapid Refresh Ensemble (HRRRE) is evaluated for its ability to forecast five Atmospheric River (AR) events that impacted California in February–March 2019. Two sets of retrospective HRRRE simulations are conducted, a control with the standard set of perturbations (initial and boundary conditions, stochastic parameters, and physics tendency), and an experiment with initial and boundary perturbations only. Reliability plots suggest the HRRRE control represents the observed Stage IV precipitation frequency well at 6-h to 24-h lead times, and rank histograms suggest the ensemble is slightly underdispersive. The HRRRE overpredicts precipitation frequency at the higher (25 mm) threshold. These results suggest the HRRRE is a useful tool to quantify probabilistic forecasts of AR events in this region. Removing stochastic physics perturbations did not substantially impact probabilistic forecasts, suggesting most of the ensemble spread is from initial and boundary condition perturbations. Spatially, ensemble precipitation coefficient of variance is lower (less forecast uncertainty) over the Sierra Nevada range than other regions, suggesting that these ensemble perturbations have a smaller impact on precipitation processes occurring over the Sierra Nevada range. More work should be conducted to understand the impacts of other model perturbations, such as microphysics, on ensemble performance, and to improve Stage IV accuracy with frozen precipitation in mountainous regions.

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.

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.

Analysis and Prediction of Summer Rainfall over Southwestern Utah

Abstract While many studies have examined intense rainfall and flash flooding during the North American Monsoon (NAM) in Arizona, Nevada, and New Mexico, less attention has focused on the NAMS’s extension into southwestern Utah. This study relates flash flood reports and Multi-Radar Multi-Sensor (MRMS) precipitation across southwestern Utah to atmospheric moisture content and instability analyses and forecasts from the High-Resolution Rapid Refresh (HRRR) model during the 2021–23 monsoon seasons. MRMS quantitative precipitation estimates over southwestern Utah during summer depend largely on the areal coverage from the KICX WSR-88D radar near Cedar City, UT. Those estimates are generally consistent with the limited number of precipitation gauge reports in the region except at extended distances from the radar. A strong relationship is evident between days with widespread precipitation and afternoons with above average precipitable water (PWAT) and convective available potential energy (CAPE) estimated from HRRR analyses across the region. Time-lagged ensembles of HRRR forecasts (initialization times from 03–06 UTC) that are 13–18 h prior to the afternoon period when convection is initiating (18–21 UTC) are useful for situational awareness of widespread precipitation events after adjusting for underprediction of afternoon CAPE. Improved skill is possible using random forest classification relying only on PWAT and CAPE to predict days experiencing excessive (upper quartile) precipitation. Such HRRR predictions may be useful for forecasters at the Salt Lake City National Weather Service Forecast Office to assist issuing flash flood potential statements for visitors to national parks and other recreational areas in the region.

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.

Evaluation of Wintertime Precipitation Estimates and Forecasts in the Mountains of Colorado

Abstract Wintertime precipitation poses many observational and forecasting challenges, especially in the complex topography of the western United States where radar beam blockage and difficulty siting in situ observations yields more sparse observations than in the eastern United States. Uncertainty in western U.S. winter precipitation is known to be high, so much so that some studies have found model simulated precipitation to produce similar or better large-scale estimates of annual precipitation than gridded observational products during climatologically anomalous years. This study evaluates high-resolution gridded precipitation estimates from Multi-Radar Multi-Sensor (MRMS) and Stage IV as well as forecasts from NOAA’s High-Resolution Rapid Refresh (HRRR) model in the Colorado Rocky Mountains. Gridded precipitation estimates and forecasts are compared with in situ SNOTEL measurements for two seasons of wintertime precipitation. The influence of forecast length, lead time, and model elevation on seasonal precipitation predictions from the HRRR are investigated. Additional comparisons are made with the relatively dense network of observations deployed in Colorado’s East River Watershed during the Study of Precipitation, the Lower Atmosphere and Surface for Hydrometeorology (SPLASH) campaign. Gridded products and forecasts are found to underestimate cold-season precipitation by 25%–65% relative to in situ and aircraft measurements, with longer forecast periods and lead times (6–24 h) having smaller biases (25%–30%) than shorter forecast periods and lead times (55%–65%). The assessment of multiple years of observations indicates that these biases are related more to the data and methods used to create the gridded products and forecasts than to precipitation characteristics. Significance Statement In the mountainous western United States, it is very challenging to both observe and forecast wintertime precipitation, yet snowfall plays an important role in providing the region’s annual water supply. This study aims to increase our understanding of the biases in observations and forecasts of snowfall in the Colorado Rocky Mountains, which can in turn impact forecasts of water availability for the ensuing warm season. In this study we find high-resolution gridded precipitation estimates and forecasts to underestimate cold-season precipitation when compared with in situ observing stations, with longer-range forecasts (e.g., daily) being the least biased. These findings were consistent over two years of study and have broad implications for the hydrologic modeling and water management communities.

NYSolarCast: A solar power forecasting system for New York State

Solar energy generation capacity will need to be greatly increased to meet aggressive clean energy targets by New York State (NYS), which are 70% renewable energy (RE) generation by 2030, and 100% by 2040. Because solar energy is a variable, weather-driven resource, accurate forecasts of solar energy generation both in nowcast (intra-day) and day-ahead time horizons are necessary for electric utilities and independent system operators, as they maintain grid stability and maximize RE use.Meeting this need for NYS, a gridded, open-source solar power forecasting system called NYSolarCast was developed (https://github.com/ncar/NYSolarCast_delivery). NYSolarCast makes 15-min resolution predictions of global horizontal irradiance (GHI) on a 3-km grid covering all of NYS, which are then used to estimate solar power generation both for select utility-scale PV plants (15-min resolution) and for zone-aggregated distributed PV (1-h resolution). The statewide GHI forecasts are made by applying the StatCast statistical forecasting model, which is trained on over two years of GHI observations from pyranometers at all 126 NYS Mesonet stations and gridded numerical weather prediction forecasts from both the Weather Research and Forecasting model tuned for solar applications (WRF-Solar®) and NOAA’s operational High-Resolution Rapid Refresh (HRRR) model. Forecasts are made at each NYS Mesonet site and then blended outward into the rest of the grid. This paper gives an overview of NYSolarCast performance for intra-day (0–6-h) GHI and power forecasts during a one-year period.

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.

Classification of Warm-Season Precipitation in High-Resolution Rapid Refresh (HRRR) Model Forecasts over the Contiguous United States

Abstract This study uses the convective adjustment time scale to identify the climatological frequency of equilibrium and nonequilibrium convection in different parts of the contiguous United States (CONUS) as modeled by the operational convection-allowing High-Resolution Rapid Refresh (HRRR) forecast system. We find a qualitatively different climatology in the northern and southern domains separated by the 40°N parallel. The convective adjustment time scale picks up the fact that convection over the northern domains is governed by synoptic flow (leading to equilibrium), while locally forced, nonequilibrium convection dominates over the southern domains. Using a machine learning algorithm, we demonstrate that the convective adjustment time-scale diagnostic provides a sensible classification that agrees with the underlying dynamics of equilibrium and nonequilibrium convection. Furthermore, the convective adjustment time scale can indicate the model quantitative precipitation forecast (QPF) quality, as it correctly reflects the higher QPF skill for precipitation under strong synoptic forcing. This diagnostic based on the strength of forcing for convection will be employed in future studies across different parts of CONUS to objectively distinguish different weather situations and explore the potential connection to warm-season precipitation predictability. Significance Statement An objective classification metric that can delineate a wide range of forecasts into distinct scenarios can serve as a valuable tool. This study represents a pioneering effort in utilizing the convective adjustment time scale to identify the climatological frequency of warm-season precipitation under varying levels of synoptic forcing in different parts of the contiguous United States (CONUS). The results demonstrate that the convective adjustment time scale is a robust metric for categorizing precipitation events and establishing a direct link to their predictability. Overall, this study provides a valuable framework for future studies focused on the CONUS domain, offering guidance on how to employ the convective adjustment time scale to classify weather regimes and explore the influence of environmental conditions on predictability of convection.

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.

Methane point source quantification using MethaneAIR: a new airborne imaging spectrometer

Abstract. The MethaneSAT satellite instrument and its aircraft precursor, MethaneAIR, are imaging spectrometers designed to measure methane concentrations with wide spatial coverage, fine spatial resolution, and high precision compared to currently deployed remote sensing instruments. At 12 960 m cruise altitude above ground (13 850 m above sea level), MethaneAIR datasets have a 4.5 km swath gridded to 10 m × 10 m pixels with 17–20 ppb standard deviation on a flat scene. MethaneAIR was deployed in the summer of 2021 in the Permian Basin to test the accuracy of the retrieved methane concentrations and emission rates using the algorithms developed for MethaneSAT. We report here point source emissions obtained during a single-blind volume-controlled release experiment, using two methods. (1) The modified integrated mass enhancement (mIME) method estimates emission rates using the total mass enhancement of methane in an observed plume combined with winds obtained from Weather Research Forecast driven by High-Resolution Rapid Refresh meteorological data in Large Eddy Simulations mode (WRF-LES-HRRR). WRF-LES-HRRR simulates winds in stochastic eddy-scale (100–1000 m) variability, which is particularly important for low-wind conditions and informing the error budget. The mIME can estimate emission rates of plumes of any size that are detectable by MethaneAIR. (2) The divergence integral (DI) method applies Gauss's theorem to estimate the flux divergence fields through a series of closed surfaces enclosing the sources. The set of boxes grows from the upwind side of the plume through the core of each plume and downwind. No selection of inflow concentration, as used in the mIME, is required. The DI approach can efficiently determine fluxes from large sources and clusters of sources but cannot resolve small point emissions. These methods account for the effects of eddy-scale variation in different ways: the DI averages across many eddies, whereas the mIME re-samples many eddies from the LES simulation. The DI directly uses HRRR winds, while mIME uses WRF-LES-HRRR wind products. Emissions estimates from both the mIME and DI methods agreed closely with the single-blind volume-controlled experiments (N = 21). The York regression between the estimated emissions and the released emissions has a slope of 0.96 [0.84, 1.08], R = 0.83 and N = 21, with 30 % mean percentage error for the whole dataset, which indicates that MethaneAIR can quantify point sources emitting more than 200 kg h−1 for the mIME and 500 kg h−1 for the DI method. The two methods also agreed on methane emission estimates from various uncontrolled sources in the Permian Basin. The experiment thus demonstrates the powerful potential of the MethaneAIR instrument and suggests that the quantification method should be transferable to MethaneSAT if it meets the design specifications.

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.

Analysis of Connected Vehicle Data to Quantify National Mobility Impacts of Winter Storms for Decision Makers and Media Reports

Traditional techniques of monitoring roadway mobility during winter weather have relied on embedded road sensors, roadside cameras, radio reports from public safety staff, or public incident reports. However, widely available connected vehicle (CV) data provides government agencies and media with a unique opportunity to monitor the mobility impact of inclement weather events in near real-time. This study presents such a use case that analyzed over 500 billion CV records characterizing the spatial and temporal impact of a winter storm that moved across the country from 21 to 26 December 2022. The analysis covered 97,000 directional miles of interstate roadway and processed over 503 billion CV records. At the storm’s peak on 22 December at 5:26 PM Eastern Time, nearly 4800 directional miles of interstate roadway were operating under 45 mph, a widely accepted indicator of degraded interstate conditions. The study presents a methodological approach to systematically assess the mobility impact of this winter event on interstate roadways at a national and regional level. The paper then looks at a case study on Interstate 70, a 4350 directional mile route passing through ten states. Statewide comparison showed Ohio was most impacted, with 9% of mile-hours operating below 45 mph on 23 December. High-Resolution Rapid Refresh weather data provided by the National Oceanic and Atmospheric Administration was integrated into the analysis to provide a visualization of the storm’s temporal path and severity. We believe the proposed metrics and visualizations are effective tools for communicating the severity and geographic impact of extreme weather events to broad non-technical audiences.

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.

Impact of Biomass Burning Organic Aerosol Volatility on Smoke Concentrations Downwind of Fires.

Biomass burning particulate matter (BBPM) affects regional air quality and global climate, with impacts expected to continue to grow over the coming years. We show that studies of North American fires have a systematic altitude dependence in measured BBPM normalized excess mixing ratio (NEMR; ΔPM/ΔCO), with airborne and high-altitude studies showing a factor of 2 higher NEMR than ground-based measurements. We report direct airborne measurements of BBPM volatility that partially explain the difference in the BBPM NEMR observed across platforms. We find that when heated to 40-45 °C in an airborne thermal denuder, 19% of lofted smoke PM1 evaporates. Thermal denuder measurements are consistent with evaporation observed when a single smoke plume was sampled across a range of temperatures as the plume descended from 4 to 2 km altitude. We also demonstrate that chemical aging of smoke and differences in PM emission factors can not fully explain the platform-dependent differences. When the measured PM volatility is applied to output from the High Resolution Rapid Refresh Smoke regional model, we predict a lower PM NEMR at the surface compared to the lofted smoke measured by aircraft. These results emphasize the significant role that gas-particle partitioning plays in determining the air quality impacts of wildfire smoke.

Using Object-Based Verification to Assess Improvements in Forecasts of Convective Storms between Operational HRRR Versions 3 and 4

Abstract The object-based verification procedure described in a recent paper by Duda and Turner was expanded herein to compare forecasts of composite reflectivity and 6-h precipitation objects between the two most recent operational versions of the High-Resolution Rapid Refresh (HRRR) model, versions 3 and 4, over an expanded set of warm season cases in 2019 and 2020. In addition to analyzing all objects, a reduced set of forecast–observation object pairs was constructed by taking the best forecast match to a given observation object for the purposes of bias-reduction and unequivocal object comparison. Despite the apparent signal of improved scalar metrics such as the object-based threat score in HRRRv4 compared to HRRRv3, no statistically significant differences were found between the models. Nonetheless, many object attribute comparisons revealed indications of improved forecast performance in HRRRv4 compared to HRRRv3. For example, HRRRv4 had a reduced overforecasting bias for medium- and large-sized reflectivity objects, and all objects during the afternoon. HRRRv4 also better replicated the distribution of object complexity and aspect ratio. Results for 6-h precipitation also suggested superior performance in HRRRv4 over HRRRv3. However, HRRRv4 was worse with centroid displacement errors and more severely overforecast objects with a high maximum precipitation amount. Overall, this exercise revealed multiple forecast deficiencies in the HRRR, which enables developers to direct development efforts on detailed and specific endeavors to improve model forecasts. Significance Statement This work builds upon the authors’ prior work in assessing model forecast quality using an alternative verification method—object-based verification. In this paper we verified two versions of the same model (one an upgrade from the other) that were making forecasts covering the same time window, using the object-based verification method. We found that the updated model was not statistically significantly better, although there were indications it performed better in certain aspects such as capturing the change in the number of storms during the daytime. We were able to identify specific problem areas in the models, which helps us direct model developers in their efforts to further improve the model.

Diagnosing Factors Leading to an Incorrect Supercell Thunderstorm Forecast

Abstract On 28 April 2019, hourly forecasts from the operational High-Resolution Rapid Refresh (HRRR) model consistently predicted an isolated supercell storm late in the day near Dodge City, Kansas, that subsequently was not observed. Two convection-allowing model (CAM) ensemble runs are created to explore the reasons for this forecast error and implications for severe weather forecasting. The 40-member CAM ensembles are run using the HRRR configuration of the WRF-ARW Model at 3-km horizontal grid spacing. The Gridpoint Statistical Interpolation (GSI)-based ensemble Kalman filter is used to assimilate observations every 15 min from 1500 to 1900 UTC, with resulting ensemble forecasts run out to 0000 UTC. One ensemble only assimilates conventional observations, and its forecasts strongly resemble the operational HRRR with all ensemble members predicting a supercell storm near Dodge City. In the second ensemble, conventional observations plus observations of WSR-88D radar clear-air radial velocities, WSR-88D diagnosed convective boundary layer height, and GOES-16 all-sky infrared brightness temperatures are assimilated to improve forecasts of the preconvective environment, and its forecasts have half of the members predicting supercells. Results further show that the magnitude of the low-level meridional water vapor flux in the moist tongue largely separates members with and without supercells, with water vapor flux differences of 12% leading to these different outcomes. Additional experiments that assimilate only radar or satellite observations show that both are important to predictions of the meridional water vapor flux. This analysis suggests that mesoscale environmental uncertainty remains a challenge that is difficult to overcome. Significance Statement Forecasts from operational numerical models are the foundation of weather forecasting. There are times when these models make forecasts that do not come true, such as 28 April 2019 when successive forecasts from the operational High-Resolution Rapid Refresh (HRRR) model predicted a supercell storm late in the day near Dodge City, Kansas, that subsequently was not observed. Reasons for this forecast error are explored using numerical experiments. Results suggest that relatively small changes to the prestorm environment led to large differences in the evolution of storms on this day. This result emphasizes the challenges to operational severe weather forecasting and the continued need for improved use of all available observations to better define the atmospheric state given to forecast models.

Evaluation of Ensemble Snowfall Forecasts Using Operationally Used Snow-to-Liquid Ratios in Five Winter Storms

Abstract The prediction of snow accumulation remains a forecasting challenge. While the adoption of ensemble numerical weather prediction has enabled the development of probabilistic guidance, the challenges associated with snow accumulation, particularly snow-to-liquid ratio (SLR), still remain when building snow-accumulation tools. In operations, SLR is generally assumed to either fit a simple mathematical relationship or conform to a historic average. In this paper, the impacts of the choice of SLR on ensemble snow forecasts are tested. Ensemble forecasts from the nine-member High-Resolution Rapid Refresh Ensemble (HRRRE) were used to create 24-h snowfall forecasts for five snowfall events associated with winter cyclones. These snowfall forecasts were derived from model liquid precipitation forecasts using five SLR relationships. These forecasts were evaluated against daily new snowfall observations from the Community Collaborative Rain Hail and Snow network. The results of this analysis show that the forecast error associated with individual members is similar to the error associated with choice of SLR. The SLR with the lowest forecast error showed regional agreement across nearby observations. This suggests that, while there is no one SLR that works best everywhere, it may be possible to improve ensemble snow forecasts if regions where SLRs perform best can be determined ahead of time. The implications of these findings for future ensemble snowfall tools will be discussed. Significance Statement Snowfall prediction remains a challenge. Computer models are used to address the inherent uncertainty in forecasts. This uncertainty includes aspects like the location and rate of snowfall. Meteorologists run multiple similar computer models to understand the range of possible weather outcomes. One aspect of uncertainty is the snow-to-liquid ratio, or the ratio of snow depth to the amount of liquid water it melts into. This study tests how common predictions of snow-to-liquid ratio impact snowfall forecasts. The results show that snow-to-liquid ratio choices are as impactful as the models’ differing snow rate or snow location forecasts, and that no particular snow-to-liquid ratio is most accurate. These results underscore the importance of better snow-to-liquid ratio prediction to improve snowfall forecasts.