The next generation Hurricane Analysis and Forecast System (HAFS) has been developed recently in the National Oceanic and Atmospheric Administration (NOAA) to accelerate the improvement of tropical cyclone (TC) forecasts within the Unified Forecast System (UFS) framework. The finite-volume cubed sphere (FV3) based convection-allowing HAFS Stand-Alone Regional model (HAFS-SAR) was successfully implemented during Hurricane Forecast Improvement Project (HFIP) real-time experiments for the 2019 Atlantic TC season. HAFS-SAR has a single large 3-km horizontal resolution regional domain covering the North Atlantic basin. A total of 273 cases during the 2019 TC season are systematically evaluated against the best track and compared with three operational forecasting systems: Global Forecast System (GFS), Hurricane Weather Research and Forecasting model (HWRF), and Hurricanes in a Multi-scale Ocean-coupled Non-hydrostatic model (HMON). HAFS-SAR has the best performance in track forecasts among the models presented in this study. The intensity forecasts are improved over GFS, but show less skill compared to HWRF and HMON. The radius of gale force wind is over-predicted in HAFS-SAR, while the hurricane force wind radius has lower error than other models.

This paper describes the vortex initialization (VI) currently used in NCEP operational hurricane models (HWRF and HMON, and possibly HAFS in the future). The VI corrects the background fields for hurricane models: it consists of vortex relocation, and size and intensity corrections. The VI creates an improved background field for the data assimilation and thereby produces an improved analysis for the operational hurricane forecast. The background field after VI can be used as an initial field (as in the HMON model, without data assimilation) or a background field for data assimilation (as in HWRF model).

Realistic wind information is critical for accurate forecasts of landfalling hurricanes. In order to provide more realistic near-surface wind forecasts of hurricanes over coastal regions, high-resolution land–sea masks are considered. As a leading hurricane modeling system, the National Centers for Environmental Prediction (NCEP) Hurricane Weather Research Forecast (HWRF) system has been widely used in both operational and research environments for studying hurricanes in different basins. In this study, high-resolution land–sea mask datasets are introduced to the nested domain of HWRF, for the first time, as an attempt to improve hurricane wind forecasts. Four destructive North Atlantic hurricanes (Harvey and Irma in 2017; and Florence and Michael in 2018), which brought historic wind damage and storm surge along the Eastern Seaboard of the United States and Northeastern Gulf Coast, were selected to demonstrate the methodology of extending the capability to HWRF, due to the introduction of the high-resolution land–sea masks into the nested domains for the first time. A preliminary assessment of the numerical experiments with HWRF shows that the introduction of high-resolution land–sea masks into the nested domains produce significantly more accurate definitions of coastlines, land-use, and soil types. Furthermore, the high-resolution land–sea mask not only improves the quality of simulated wind information along the coast, but also improves the hurricane track, intensity, and storm-size predictions.

The goal of this paper is to introduce a new multi-storm atmosphere/ocean coupling scheme that was implemented and tested in the Basin-Scale Hurricane Weather Research and Forecasting (HWRF-B) model. HWRF-B, an experimental model developed at the National Oceanic and Atmospheric Administration (NOAA) and supported by the Hurricane Forecast Improvement Program, is configured with multiple storm-following nested domains to produce high-resolution predictions for several tropical cyclones (TCs) within the same forecast integration. The new coupling scheme parallelizes atmosphere/ocean interactions for each nested domain in HWRF-B, and it may be applied to any atmosphere/ocean coupled system. TC forecasts from this new hydrodynamical modeling system were produced in the North Atlantic and eastern North Pacific from 2017–2019. The performance of HWRF-B was evaluated, including forecasts of TC track, intensity, structure (e.g., surface wind radii), and intensity change, and simulated sea-surface temperatures were compared with satellite observations. Median forecast skill scores showed significant improvement over the operational HWRF at most forecast lead times for track, intensity, and structure. Sea-surface temperatures cooled by 1–8 °C for the five HWRF-B case studies, demonstrating the utility of the model to study the impact of the ocean on TC intensity forecasting. These results show the value of a multi-storm modeling system and provide confidence that the multi-storm coupling scheme was implemented correctly. Future TC models within NOAA, especially the Unified Forecast System’s Hurricane Analysis and Forecast System, would benefit from the multi-storm coupling scheme whose utility and performance are demonstrated in HWRF-B here.

A statistical analysis is performed on the high-frequency (3⅓ s) output from NOAA’s cloud-permitting, high-resolution operational Hurricane Weather Research and Forecasting (HWRF) Model for all tropical cyclones (TCs) in the North Atlantic Ocean basin over a 3-yr period (2017–19). High-frequency HWRF forecasts of TC track and 10-m maximum wind speed (Vmax) exhibited large fluctuations that were not captured by traditional low-frequency (6 h) model output. Track fluctuations were inversely proportional to Vmax, with average values of 6–8 km. The Vmax fluctuations were as high as 20 kt (10.3 m s⁻¹) in individual forecasts and were a function of maximum intensity, with a standard deviation of 5.5 kt (2.8 m s⁻¹) for category-2 hurricanes and smaller fluctuations for tropical storms and major hurricanes. The radius of Vmax contracted or remained steady when TCs rapidly intensified in high-frequency HWRF forecasts, consistent with observations. Running-mean windows of 3–9 h were applied at synoptic times to smooth the high-frequency HWRF output to investigate its utility to operational forecasting. Smoothed high-frequency HWRF output improved Vmax forecast skill by up to 8% and produced a more realistic distribution of 6-h intensity change when compared with low-frequency, instantaneous output. Furthermore, the high-frequency track forecast output may be useful for investigating characteristics of TC trochoidal motions.

The National Oceanic and Atmospheric Administration’s (NOAA) cloud-permitting high-resolution operational Hurricane Weather and Research Forecasting (HWRF) model includes the sophisticated hybrid grid-point statistical interpolation (GSI) and Ensemble Kalman Filter (EnKF) data assimilation (DA) system, which allows assimilating high-resolution aircraft observations in tropical cyclone (TC) inner core regions. In the operational HWRF DA system, the flow-dependent background error covariance matrix is calculated from the HWRF self-cycled 40-member ensemble. This DA system has proved to provide improved initial TC structure and therefore improved TC track and intensity forecasts. However, the uncertainties from the model physics are not taken into account in the FY2017 version of the HWRF DA system. In order to further improve the HWRF DA system, the stochastic physics perturbations are introduced in the HWRF DA, including the cumulus convection scheme, the planetary boundary layer (PBL) scheme, and model surface physics (drag coefficient), for HWRF-based ensembles. This study shows that both TC initial conditions and TC track and intensity forecast skills are improved by adding stochastic model physics in the HWRF self-cycled DA system. It was found that the improvements in the TC initial conditions and forecasts are the results of ensemble spread increases which realistically represent the model background error covariance matrix in HWRF DA. For all 2016 Atlantic storms, the TC track and intensity forecast skills are improved by about ~3% and 6%, respectively, compared to the control experiment. The case study shows that the stochastic physics in HWRF DA is especially helpful for those TCs that have inner-core high-resolution aircraft observations, such as tail Doppler radar (TDR) data.