Evan A. Kalina

To assess the performance and scalability of the Unified Forecast System (UFS) Short-Range Weather (SRW) application, case studies are chosen to cover a wide variety of forecast applications. Here, model forecasts of Hurricane Barry (July 2019) are examined and analyzed. Several versions of the Global Forecast System (GFS) and Rapid Refresh Forecast System (RRFS) physics suites are run in the UFS-SRW at grid spacings of 25 km, 13 km, and 3 km. All model configurations produce significant track errors of up to 350 km at landfall. The track errors are investigated, and several commonalities are seen between model configurations. A westerly bias in the environmental steering flow surrounding the tropical cyclone (TC) is seen across forecasts, and this bias is coincident with a warm sea surface temperature (SST) bias and overactive convection on the eastern side of the forecasted TC. Positive feedback between the surface winds, latent heating, moisture, convection, and TC intensification is initiated by this SST bias. The asymmetric divergent flow induced by the excess convection results in all model TC tracks being diverted to the east as compared to the track derived from reanalysis. The large differences between runs using the same physics packages at different grid spacing suggest a deficiency in the scalability of these packages with respect to hurricane forecasting in vertical wind shear.

The HAFS model is an effort under the NGGPS and UFS initiatives to create the next generation of hurricane prediction and analysis system based on FV3-GFS. It has been validated extensively using traditional verification indicators such as tracker error and biases, intensity error and biases, and the radii of gale, damaging and hurricane strength winds. While satellite images have been used to verify hurricane model forecasts, they have not been used on HAFS. The community radiative transfer model CRTM is used to generate model synthetic satellite images from HAFS model forecast state variables. The 24 forecast snapshots in the mature stage of hurricane Dorian in 2019 are used to generate a composite model synthetic GOES-R infrared brightness image. The composite synthetic image is compared to the corresponding composite image generated from the observed GOES-R data, to evaluate the model forecast TC vortex intensity, size, and asymmetric structure. Results show that the HAFS forecast TC Dorian agrees reasonably well with the observation, but the forecast intensity is weaker, its overall vortex size smaller, and the radii of its eye and maximum winds larger than the observed. The evaluation results can be used to further improve the model. While these results are consistent with those obtained by traditional verification methods, evaluations based on composite satellite images provide an additional benefit with richer information because they have near-real-times spatially and temporally continuous high-resolution data with global coverage. Composite satellite infrared images could be used routinely to supplement traditional verification methods in the HAFS and other hurricane model evaluations. Note since this study only evaluated one hurricane, the above conclusions are only applicable to the model behavior of the mature stage of hurricane Dorian in 2019, and caution is needed to extend these conclusions to expect model biases in predicting other TCs. Nevertheless, the consistency between the evaluation using composite satellite images and the traditional metrics, of hurricane Dorian, shows that this method has the potential to be applied to other storms in future studies.

The High-Resolution Rapid Refresh Ensemble (HRRRE) is a 36-member ensemble analysis system with 9 forecast members that utilizes the Advanced Research version of the Weather Research and Forecasting (ARW-WRF) dynamic core and the physics suite from the operational Rapid Refresh/High-Resolution Rapid Refresh deterministic modeling system. A goal of HRRRE development is a system with sufficient spread among members, comparable in magnitude to the random error in the ensemble mean, to represent the range of possible future atmospheric states. HRRRE member diversity has traditionally been obtained by perturbing the initial and lateral boundary conditions of each member, but recent development has focused on implementing stochastic approaches in HRRRE to generate additional spread. These techniques were tested in retrospective experiments and in the May 2019 Hazardous Weather Testbed Spring Experiment (HWT-SE). Results show a 6%–25% increase in the ensemble spread in 2-m temperature, 2-m mixing ratio, and 10-m wind speed when stochastic parameter perturbations are used in HRRRE (HRRRE-SPP). Case studies from HWT-SE demonstrate that HRRRE-SPP performed similar to or better than the operational High-Resolution Ensemble Forecast system, version 2 (HREFv2), and the nonstochastic HRRRE. However, subjective evaluations provided by HWT-SE forecasters indicated that overall, HRRRE-SPP predicted lower probabilities of severe weather (using updraft helicity as a proxy) compared to HREFv2. A statistical analysis of the performance of HRRRE-SPP and HREFv2 from the 2019 summer convective season supports this claim, but also demonstrates that the two systems have similar reliability for prediction of severe weather using updraft helicity.

The intensity and structure of simulated tropical cyclones (TCs) are known to be sensitive to the planetary boundary layer (PBL) parameterization in numerical weather prediction models. In this paper, we use an idealized version of the Hurricane Weather Research and Forecast system (HWRF) with constant sea-surface temperature (SST) to examine how the configuration of the PBL scheme used in the operational HWRF affects TC intensity change (including rapid intensification) and structure. The configuration changes explored in this study include disabling non-local vertical mixing, changing the coefficients in the stability functions for momentum and heat, and directly modifying the Prandtl number (Pr), which controls the ratio of momentum to heat and moisture exchange in the PBL. Relative to the control simulation, disabling non-local mixing produced a ~15% larger storm that intensified more gradually, while changing the coefficient values used in the stability functions had little effect. Varying Pr within the PBL had the greatest impact, with the largest Pr (~1.6 versus ~0.8) associated with more rapid intensification (~38 versus 29 m s⁻¹ per day) but a 5–10 m s⁻¹ weaker intensity after the initial period of strengthening. This seemingly paradoxical result is likely due to a decrease in the radius of maximum wind (~15 versus 20 km), but smaller enthalpy fluxes, in simulated storms with larger Pr. These results underscore the importance of measuring the vertical eddy diffusivities of momentum, heat, and moisture under high-wind, open-ocean conditions to reduce uncertainty in Pr in the TC PBL.

This paper reviews the evolution of planetary boundary layer (PBL) parameterization schemes that have been used in the operational version of the Hurricane Weather Research and Forecasting (HWRF) model since 2011. Idealized simulations are then used to evaluate the effects of different PBL schemes on hurricane structure and intensity. The original Global Forecast System (GFS) PBL scheme in the 2011 version of HWRF produces the weakest storm, while a modified GFS scheme using a wind-speed dependent parameterization of vertical eddy diffusivity (K_m) produces the strongest storm. The subsequent version of the hybrid eddy diffusivity and mass flux scheme (EDMF) used in HWRF also produces a strong storm, similar to the version using the wind-speed dependent K_m. Both the intensity change rate and maximum intensity of the simulated storms vary with different PBL schemes, mainly due to differences in the parameterization of K_m. The smaller the K_m in the PBL scheme, the faster a storm tends to intensify. Differences in hurricane PBL height, convergence, inflow angle, warm-core structure, distribution of deep convection, and agradient force in these simulations are also examined. Compared to dropsonde and Doppler radar composites, improvements in the kinematic structure are found in simulations using the wind-speed dependent K_m and modified EDMF schemes relative to those with earlier versions of the PBL schemes in HWRF. However, the upper boundary layer in all simulations is much cooler and drier than that in dropsonde observations. This model deficiency needs to be considered and corrected in future model physics upgrades.

Thunderstorms that produce hail accumulations at the surface can impact residents by obstructing roadways, closing airports, and causing localized flooding from hail-clogged drainages. These storms have recently gained an increased interest within the scientific community. However, differences that are observable in real time between these storms and storms that produce nonimpactful hail accumulations have yet to be documented. Similarly, the characteristics within a single storm that are useful to quantify or predict hail accumulations are not fully understood. This study uses lightning and dual-polarization radar data to characterize hail accumulations from three storms that occurred on the same day along the Colorado–Wyoming Front Range. Each storm’s characteristics are verified against radar-derived hail accumulation maps and in situ observations. The storms differed in maximum accumulation, either producing 22 cm, 7 cm, or no accumulation. The magnitude of surface hail accumulations is found to be dependent on a combination of in-cloud hail production, storm translation speed, and hailstone melting. The optimal combination for substantial hail accumulations is enhanced in-cloud hail production, slow storm speed, and limited hailstone melting. However, during periods of similar in-cloud hail production, lesser accumulations are derived when storm speed and/or hailstone melting, identified by radar presentation, is sufficiently large. These results will aid forecasters in identifying when hail accumulations are occurring in real time.

The Developmental Testbed Center (DTC) tested two convective parameterization schemes in the Hurricane Weather Research and Forecasting (HWRF) Model and compared them in terms of performance of forecasting tropical cyclones (TCs). Several TC forecasts were conducted with the scale-aware Simplified Arakawa Schubert (SAS) and Grell–Freitas (GF) convective schemes over the Atlantic basin. For this sample of over 100 cases, the storm track and intensity forecasts were superior for the GF scheme compared to SAS. A case study showed improved storm structure for GF when compared with radar observations. The GF run had increased inflow in the boundary layer, which resulted in higher angular momentum. An angular momentum budget analysis shows that the difference in the contribution of the eddy transport to the total angular momentum tendency is small between the two forecasts. The main difference is in the mean transport term, especially in the boundary layer. The temperature tendencies indicate higher contribution from the microphysics and cumulus heating above the boundary layer in the GF run. A temperature budget analysis indicated that both the temperature advection and diabatic heating were the dominant terms and they were larger near the storm center in the GF run than in the SAS run. The above results support the superior performance of the GF scheme for TC intensity forecast.

Unique data from seven flights of the Coyote small unmanned aircraft system (sUAS) were collected in Hurricanes Maria (2017) and Michael (2018). Using NOAA’s P-3 reconnaissance aircraft as a deployment vehicle, the sUAS collected high-frequency (>1 Hz) measurements in the turbulent boundary layer of hurricane eyewalls, including measurements of wind speed, wind direction, pressure, temperature, moisture, and sea surface temperature, which are valuable for advancing knowledge of hurricane structure and the process of hurricane intensification. This study presents an overview of the sUAS system and preliminary analyses that were enabled by these unique data. Among the most notable results are measurements of turbulence kinetic energy and momentum flux for the first time at low levels (<150 m) in a hurricane eyewall. At higher altitudes and lower wind speeds, where data were collected from previous flights of the NOAA P-3, the Coyote sUAS momentum flux values are encouragingly similar, thus demonstrating the ability of an sUAS to measure important turbulence properties in hurricane boundary layers. Analyses from a large-eddy simulation (LES) are used to place the Coyote measurements into context of the complicated high-wind eyewall region. Thermodynamic data are also used to evaluate the operational HWRF model, showing a cool, dry, and thermodynamically unstable bias near the surface. Preliminary data assimilation experiments also show how sUAS data can be used to improve analyses of storm structure. These results highlight the potential of sUAS operations in hurricanes and suggest opportunities for future work using these promising new observing platforms.

In recent years, hail accumulations from thunderstorms have occurred frequently enough to catch the attention of the National Weather Service, the general public, and news agencies. Despite the extreme nature of these thunderstorms, no mechanism is currently in place to obtain adequate reports, measurements, or forecasts of accumulated hail depth. To better identify and forecast hail accumulations, the Colorado Hail Accumulation from Thunderstorms (CHAT) project was initiated in 2016 with the goals of collecting improved and more frequent hail depth reports on the ground as well as studying characteristics of storms that produce hail accumulations in Colorado. A desired outcome of this research is to identify predictors for hail-producing thunderstorms typically occurring along the Colorado Front Range that might be used as operational nowcast products in the future. During the 2016 convective season, we asked amateur meteorologists to send general information, photos, and videos on hail depth using social media. They submitted over 58 reports in Colorado with information on location, time, depth, and areal coverage of hail accumulations. We have analyzed dual-polarization radar and lightning mapping array data from 32 thunderstorms in Colorado, which produced between 0.5 and 50 cm of hail accumulation on the ground, to identify characteristics unique to storms with hail accumulations. This preliminary analysis shows how enhanced in-cloud hail presence and surface accumulation can be tracked throughout the lifetime of a thunderstorm using dual-polarization radar and lightning data, and how hail accumulation events are associated with large in-cloud ice water content, long hailfall duration, or a combination of these.

Thunderstorms that produce surface hail accumulations, sometimes as large as 60 cm in depth, have significantly affected the residents of the Front Range and High Plains of Colorado and Wyoming by creating hazardous road conditions and endangering lives and property. To date, surface hail accumulation is not part of a routine forecasting or monitoring system. Extensive coordinated hail accumulation reports and operational products designed to identify deep hail accumulating storms in real time are lacking. Kalina et al. used dual-polarization WSR-88D radar observations to calculate hail depth and hail accumulations but never validated the algorithm. This study shows how 20 quality-controlled hail depth reports from the hail depth database built by the Colorado Hail Accumulation from Thunderstorms (CHAT) project are being used to validate the Kalina et al. radar-based hail accumulation algorithm for operational application. The validated algorithm shows increased correlations between radar-derived and reported accumulations for hail depth reports not included in the validation. Furthermore, increases in computational efficiency have allowed the improved algorithm to be used operationally. With an improved hail accumulation algorithm, thunderstorms that produce hail accumulations are more frequently detected than previously reported.