Vertical eddy diffusivity (VED) in the planetary boundary layer (PBL) has a significant impact on forecasts of tropical cyclone (TC) structure and intensity. VED uncertainties in PBL parameterizations can be partly attributed to the model’s inability to represent roll vortices (RV). In this study, RV effects on turbulent fluxes derived from a large eddy simulation (LES) by Li et al. (Geophys. Res. Lett., 2021, 48, e2020GL090703) are added to the VED parameterization of the PBL scheme within the operational Hurricane Weather Research and Forecasting (HWRF) model. RV contribution to VED is parameterized through a coefficient and varies with the RV intensity and velocity scale. A modification over land has also been implemented. This modified VED parameterization is compared with the original wind-speed-dependent VED scheme in HWRF. Retrospective HWRF forecasts of Hurricanes Florence (2018) and Laura (2020) are analyzed to evaluate the impacts of the modified VED scheme on landfalling hurricane forecasts. Results show that the modified PBL scheme with the RV effect leads to an improvement in 10-m maximum wind speed forecasts of 14%–31%, with a neutral to positive improvement for track forecasts. Improved wind structure and precipitation forecasts against observations are also noted with the modified PBL scheme. Further diagnoses indicate that the revised PBL scheme enhances moist entropy in the boundary layer over land, leading to improved TC intensity prediction compared to the original scheme.