On the sensitivity of hurricane intensity and structure to horizontal tracer advection schemes in FV3

We investigate the sensitivity of hurricane intensity and structure to the horizontal tracer advection in the Geophysical Fluid Dynamics Laboratory (GFDL) Finite-Volume Cubed-Sphere Dynamical Core (FV3). We compare two schemes, a monotonic scheme and a less diffusive positive-definite scheme. The positive-definite scheme leads to significant improvement in the intensity prediction relative to the monotonic scheme in a suite of five-day forecasts that mostly consist of rapidly intensifying hurricanes. Notable storm structural differences are present: the radius of maximum wind (RMW) is smaller and eyewall convection occurs farther inside the RMW when the positive-definite scheme is used. Moreover, we find that the horizontal tracer advection scheme affects the eyewall convection location by affecting the moisture distribution in the inner-core region. This study highlights the importance of dynamical core algorithms in hurricane intensity prediction.

Rapid observations of ocean dynamics and stratification along a steep island coast during Hurricane María

Hurricanes are extreme storms that affect coastal communities, but the linkages between hurricane forcing and ocean dynamics remain poorly understood. Here, we present full water column observations at unprecedented resolution from the southwest Puerto Rico insular shelf and slope during Hurricane María, representing a rare set of high-frequency, subsurface, oceanographic observations collected along an island margin during a hurricane. The shelf geometry and orientation relative to the storm acted to stabilize and strengthen stratification. This maintained elevated sea-surface temperatures (SSTs) throughout the storm and led to an estimated 65% greater potential hurricane intensity contribution at this site before eye passage. Coastal cooling did not occur until 11 hours after the eye passage. Our findings present a new framework for how hurricane interaction with insular island margins may generate baroclinic processes that maintain elevated SSTs, thus potentially providing increased energy for the storm.

Understanding Error Distributions of Hurricane Intensity Forecasts during Rapid Intensity Changes

The characteristics of official National Hurricane Center (NHC) intensity forecast errors are examined for the North Atlantic and east Pacific basins from 1989 to 2018. It is shown how rapid intensification (RI) and rapid weakening (RW) influence yearly NHC forecast errors for forecasts between 12 and 48 h in length. In addition to being the tail of the intensity change distribution, RI and RW are at the tails of the forecast error distribution. Yearly mean absolute forecast errors are positively correlated with the yearly number of RI/RW occurrences and explain roughly 20% of the variance in the Atlantic and 30% in the east Pacific. The higher occurrence of RI events in the east Pacific contributes to larger intensity forecast errors overall but also a better probability of detection and success ratio. Statistically significant improvements to 24-h RI forecast biases have been made in the east Pacific and to 24-h RW biases in the Atlantic. Over-ocean 24-h RW events cause larger mean errors in the east Pacific that have not improved with time. Environmental predictors from the Statistical Hurricane Intensity Prediction Scheme (SHIPS) are used to diagnose what conditions lead to the largest RI and RW forecast errors on average. The forecast error distributions widen for both RI and RW when tropical systems experience low vertical wind shear, warm sea surface temperature, and moderate low-level relative humidity. Consistent with existing literature, the forecast error distributions suggest that improvements to our observational capabilities, understanding, and prediction of inner-core processes is paramount to both RI and RW prediction.

Estimation of Hurricane Intensity from ATMS-Derived Temperature Anomaly using Machine Learning

The warm-core structure is one of the basic characteristics that vary during the different stages of tropical cyclones (TCs). The warm core structure of the TCs during2016-2019 over the Atlantic Ocean was derived based on the observations of the ATMS onboard S-NPP. From linear regression, the mean prediction error (MPE) is 39.04 mph for Vmax and 14.47 hPa for Pmin. The root-mean-square error(RMSE) is 42.70 mph for the maximum sustained wind (Vmax) and 77.69 hPa for the minimum sea-level pressure (Pmin). Several machine learning (ML) techniques are used to develop the Atlantic TC intensity (Vmax and Pmin) estimation models. The support vector machine (SVM) model has the best performance with the MPE of 14.62 mph for Vmaxan 7.66 hPa for Pmin, and the RMSE of 19.91 mph for Vmax and 10.58 hPa for Pmin. Adding latitude and day of year (DOY) can further improve the estimation of Vmax by decreasing MPE to 13.01mph and RME to 17.33 mph using SVM. Best estimation of Pminoccurs when adding the day of year to the training process, as the MPE is 7.23 hPa and RMS is 9.88 hPa. Other TC information, such as longitude and local time, does not help to improve the performance of the hurricane intensity estimation models significantly.

Evaluation of the Four-Dimensional Ensemble-Variational Hybrid Data Assimilation with Self-Consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the Gridpoint Statistical Interpolation (GSI)-based four-dimensional ensemble-variational (GSI-4DEnVar) hybrid data assimilation system for the Hurricane Weather Research and Forecasting model (HWRF) model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-based three-dimensional ensemble-variational (GSI-3DEnVar) hybrid data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) ensemble Kalman filtering data assimilation system. In addition, our proposed initialization framework discards the empirical intensity correction in the vortex initialization package that is employed by the GSI-3DEnVar initialization framework in operational HWRF. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, minimum sea level pressure (MSLP), and maximum surface wind (MSW) forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. With assimilating high-resolution observations within the hurricane inner-core region, GSI-4DEnVar can produce the initial hurricane intensity reasonably well without the empirical vortex intensity correction. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts.

Evaluation of the Four-dimensional Ensemble-Variational Hybrid Data Assimilation with Self-consistent Regional Background Error Covariance for Improved Hurricane Intensity Forecasts

The feasibility of a hurricane initialization framework based on the GSI-4DEnVar data assimilation system for the HWRF model is evaluated in this study. The system considers the temporal evolution of error covariances via the use of four-dimensional ensemble perturbations that are provided by high-resolution, self-consistent HWRF ensemble forecasts. It is different from the configuration of the GSI-3DEnVar data assimilation system, similar to that used in the operational HWRF, which employs background error covariances provided by coarser-resolution global ensembles from the NCEP GFS ensemble Kalman filtering data assimilation system. Data assimilation and numerical simulation experiments for Hurricanes Joaquin (2015), Patricia (2015), and Matthew (2016) are conducted during their intensity changes. The impacts of two initialization frameworks on the HWRF analyses and forecasts are compared. It is found that GSI-4DEnVar leads to a reduction in track, MSLP, and MSW forecast errors in all of the HWRF simulations, compared with the GSI-3DEnVar initialization framework. Further diagnoses with Hurricane Joaquin indicate that GSI-4DEnVar can significantly alleviate the imbalances in the initial conditions and enhance the performance of the data assimilation and subsequent hurricane intensity and precipitation forecasts.