Ocean component of the first operational version of Hurricane Analysis and Forecast System: Evaluation of HYbrid Coordinate Ocean Model and hurricane feedback forecasts

The first operational version of the coupled Hurricane Analysis and Forecast System (HAFSv1) launched in 2023 consists of the HYbrid Coordinate Ocean Model (HYCOM) and finite-volume cubed-sphere (FV3) dynamic atmosphere model. This system is a product of efforts involving improvements and updates over a 4-year period (2019-2022) through extensive collaborations between the Environmental Modeling Center at the US National Centers for Environmental Prediction (NCEP) and NOAA Atlantic Oceanography and Meteorology Laboratory. To provide two sets of numerical guidance, the initial operational capability of HAFSv1 was configured to two systems-HFSA and HFSB. In this study, we present in-depth analysis of the forecast skills of the upper ocean that was co-evolved by the HFSA and HFSB. We chose hurricane Laura (2020) as an example to demonstrate the interactions between the storm and oceanic mesoscale features. Comparisons performed with the available in situ observations from gliders as well as Argos and National Data Buoy Center moorings show that the HYCOM simulations have better agreement for weak winds than high winds (greater than Category 2). The skill metrics indicate that the model sea-surface temperature (SST) and mixed layer depth (MLD) have a relatively low correlation. The SST, MLD, mixed layer temperature (MLT), and ocean heat content (OHC) are negatively biased. For high winds, SST and MLT are more negative, while MLD is closer to the observations with improvements of about 8%-19%. The OHC discrepancy is proportional to predicted wind intensity. Contrarily, the mixed layer salinity (MLS) uncertainties are smaller and positive for higher winds, probably owing to the higher MLD. The less-negative bias of MLD for high winds implies that the wind-force mixing is less effective owing to the higher MLD and high buoyancy stability (approx. 1.5-1.7 times) than the observations. The heat budget analysis suggests that the maximum heat loss by hurricane Laura was O(< 3 degrees C per day). The main contributor here is advection, followed by entrainment, which act against or with each other depending on the storm quadrant. We also found relatively large unaccountable heat residuals for the in-storm period, and the residuals notably led the heat tendency, meaning that further improvements of the subscale simulations are warranted. In summary, HYCOM simulations showed no systematic differences forced by either HFSA or HFSB.