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).

This document summarizes the physics schemes used in two configurations of the first version of the operational Hurricane Analysis and Forecast System (HAFSv1) at NOAA NCEP. The physics package in HAFSv1 is the same as that used in NCEP global forecast system (GFS) version 16 except for an additional microphysics scheme and modifications to sea surface roughness lengths, boundary layer scheme, and the entrainment rate in the deep convection scheme. Those modifications are specifically designed for improving the simulation of tropical cyclones (TCs). The two configurations of HAFSv1 mainly differ in the adopted microphysics schemes and TC-specific modifications in addition to model initialization. Experiments are made to highlight the impacts of TC-specific modifications and different microphysics schemes on HAFSv1 performance. Challenges and developmental plans of physics schemes for future versions of operational HAFS are discussed.

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.

A modification to the mixing length formulation in a planetary boundary layer (PBL) scheme is introduced to improve the intensity forecast of tropical cyclones (TCs) in a basin-scale Hurricane Analysis and Forecast System (HAFS) for the real-time experiment in 2021. The 2020 basin-scale HAFS with the physics suite of the NCEP operational Global Forecast System performs well in terms of the reduced root-mean-square (RMS) errors in track and intensity except for the mean intensity bias, compared with NCEP operational hurricane models. To address the large intensity bias issue, the vertical mixing length near the surface used in the PBL scheme is increased to follow the similarity theory, consistent with that used in the surface layer scheme. Test results show that the RMS error and bias in intensity are further reduced without the degradation of the track forecast. An idealized one-dimensional TC PBL model is used to understand the model response to the modification, indicating that the radial wind is strengthened to dynamically balance the enhanced downward momentum mixing. This is also exhibited in the case study of a three-dimensional HAFS simulation, with the improved vertical distribution of the simulated wind speed in the eyewall area. Given the improvement, the modification has been implemented in one of the configurations of the first version of the operational HAFS at NCEP. Finally, the adjustment of the parameterization of diffusion and mixing in TC simulations is discussed. Significance Statement: A modification to the mixing length formulation in a PBL scheme is described, which improves the intensity forecast of tropical cyclones simulated in the Hurricane Analysis and Forecast System (HAFS). Retrospective tests indicate that the modification can reduce the root-mean-square error and bias of the simulated TC intensity by 5%–10% and 50%, respectively. This modification has been implemented in one of the operational configurations of HAFS, version 1, at NCEP, improving the hurricane model guidance.

A new physically based horizontal mixing-length formulation is introduced and evaluated in the Hurricane Weather and Research Forecasting (HWRF) Model. Recent studies have shown that the structure and intensity of tropical cyclones (TCs) simulated by numerical models are sensitive to horizontal mixing length in the parameterization of horizontal diffusion. Currently, many numerical models including the operational HWRF Model formulate the horizontal mixing length as a fixed fraction of grid spacing or a constant value, which is not realistic. To improve the representation of the horizontal diffusion process, the new formulation relates the horizontal mixing length to local wind and its horizontal gradients. The resulting horizontal mixing length and diffusivity are much closer to those derived from field measurements. To understand the impact of different mixing-length formulations, we analyze the evolutions of an idealized TC simulated by the HWRF Model with the new formulation and with the current formulation (i.e., constant values) of horizontal mixing length. In two real-case tests, the HWRF Model with the new formulation produces the intensity and track forecasts of Hurricanes Harvey (2017) and Lane (2018) that are much closer to observations. Retrospective runs of hundreds of forecast cycles of multiple hurricanes show that the mean errors in intensity and track simulated by HWRF with the new formulation can be reduced approximately by 10%.