This paper reviews the evolution of planetary boundary layer (PBL) parameterization schemes that have been used in the operational version of the Hurricane Weather Research and Forecasting (HWRF) model since 2011. Idealized simulations are then used to evaluate the effects of different PBL schemes on hurricane structure and intensity. The original Global Forecast System (GFS) PBL scheme in the 2011 version of HWRF produces the weakest storm, while a modified GFS scheme using a wind-speed dependent parameterization of vertical eddy diffusivity (K_m) produces the strongest storm. The subsequent version of the hybrid eddy diffusivity and mass flux scheme (EDMF) used in HWRF also produces a strong storm, similar to the version using the wind-speed dependent K_m. Both the intensity change rate and maximum intensity of the simulated storms vary with different PBL schemes, mainly due to differences in the parameterization of K_m. The smaller the K_m in the PBL scheme, the faster a storm tends to intensify. Differences in hurricane PBL height, convergence, inflow angle, warm-core structure, distribution of deep convection, and agradient force in these simulations are also examined. Compared to dropsonde and Doppler radar composites, improvements in the kinematic structure are found in simulations using the wind-speed dependent K_m and modified EDMF schemes relative to those with earlier versions of the PBL schemes in HWRF. However, the upper boundary layer in all simulations is much cooler and drier than that in dropsonde observations. This model deficiency needs to be considered and corrected in future model physics upgrades.

The Developmental Testbed Center used the Hurricane Weather Research and Forecasting (HWRF) system to test the sensitivity of tropical cyclone track and intensity forecasts to different convective schemes. A control configuration that employed the HWRF Simplified Arakawa Scheme (SAS) was compared with the Kain‐Fritsch and Tiedtke schemes, as well as with a newer implementation of the SAS. A comprehensive test for Atlantic and Eastern North Pacific storms shows that the SAS scheme produces the best track forecasts. Even though the convective parameterization was absent in the inner 3 km nest, the intensity forecasts are sensitive to the choice of cumulus scheme on the outer grids. The impact of convective‐scale heating on the environmental flow accumulates in time since the hurricane vortex is cycled in the HWRF model initialization. This study shows that, for a given forecast, the sensitivity to cumulus parameterization combines the influence of physics and initial conditions.

One of the well-known challenges of fog forecasting is the high spatio-temporal variability of fog. An ensemble forecast aims to capture this variability by representing the uncertainty in the initial/lateral boundary conditions (ICs/BCs) and model physics. The present study highlights a new operational Ensemble Forecast System (EFS) developed by the Indian Institute of Tropical Meteorology (IITM), Pune, to predict the fog over the Indo-Gangetic Plain (IGP) region using the visibility (Vis) diagnostic algorithm. The EFS framework comprises the WRF model with a 4 km horizontal resolution, initialized by 21 ICs/BCs. The advantages of probabilistic fog forecasting have been demonstrated by comparing control (CNTL) and ensemble-based fog forecasts. The forecast is verified using fog observations from the Indira Gandhi International (IGI) airport during the winter months of 2020–2021 and 2021–2022. The results show that with a probability threshold of 50%, the ensemble forecasts perform better than the CNTL forecasts. The skill scores of EFS are relatively promising, with a Hit Rate of 0.95 and a Critical Success Index of 0.55; additionally, the False Alarm Rate and Missing Rate are low, with values of 0.43 and 0.04, respectively. The EFS could correctly predict more fog events (37 out of 39) compared with the CNTL forecast (31 out of 39) and shows the potential skill. Furthermore, EFS has a substantially reduced error in predicting fog onset and dissipation (mean onset and dissipation error of 1 h each) compared to the CNTL forecasts.

The intensity and structure of simulated tropical cyclones (TCs) are known to be sensitive to the planetary boundary layer (PBL) parameterization in numerical weather prediction models. In this paper, we use an idealized version of the Hurricane Weather Research and Forecast system (HWRF) with constant sea-surface temperature (SST) to examine how the configuration of the PBL scheme used in the operational HWRF affects TC intensity change (including rapid intensification) and structure. The configuration changes explored in this study include disabling non-local vertical mixing, changing the coefficients in the stability functions for momentum and heat, and directly modifying the Prandtl number (Pr), which controls the ratio of momentum to heat and moisture exchange in the PBL. Relative to the control simulation, disabling non-local mixing produced a ~15% larger storm that intensified more gradually, while changing the coefficient values used in the stability functions had little effect. Varying Pr within the PBL had the greatest impact, with the largest Pr (~1.6 versus ~0.8) associated with more rapid intensification (~38 versus 29 m s⁻¹ per day) but a 5–10 m s⁻¹ weaker intensity after the initial period of strengthening. This seemingly paradoxical result is likely due to a decrease in the radius of maximum wind (~15 versus 20 km), but smaller enthalpy fluxes, in simulated storms with larger Pr. These results underscore the importance of measuring the vertical eddy diffusivities of momentum, heat, and moisture under high-wind, open-ocean conditions to reduce uncertainty in Pr in the TC PBL.

The Developmental Testbed Center (DTC) tested two convective parameterization schemes in the Hurricane Weather Research and Forecasting (HWRF) Model and compared them in terms of performance of forecasting tropical cyclones (TCs). Several TC forecasts were conducted with the scale-aware Simplified Arakawa Schubert (SAS) and Grell–Freitas (GF) convective schemes over the Atlantic basin. For this sample of over 100 cases, the storm track and intensity forecasts were superior for the GF scheme compared to SAS. A case study showed improved storm structure for GF when compared with radar observations. The GF run had increased inflow in the boundary layer, which resulted in higher angular momentum. An angular momentum budget analysis shows that the difference in the contribution of the eddy transport to the total angular momentum tendency is small between the two forecasts. The main difference is in the mean transport term, especially in the boundary layer. The temperature tendencies indicate higher contribution from the microphysics and cumulus heating above the boundary layer in the GF run. A temperature budget analysis indicated that both the temperature advection and diabatic heating were the dominant terms and they were larger near the storm center in the GF run than in the SAS run. The above results support the superior performance of the GF scheme for TC intensity forecast.