The fifth-generation Pennsylvania State University–National Center for Atmospheric Research Mesoscale Model is used to simulate Hurricane Bob (1991) using grids nested to high resolution (4 km). Tests are conducted to determine the sensitivity of the simulation to the available planetary boundary layer parameterizations, including the bulk aerodynamic, Blackadar, Medium-Range Forecast (MRF) model, and Burk–Thompson boundary layer schemes. Significant sensitivity is seen, with minimum central pressures varying by up to 16 mb and maximum winds by 15 m s−1. The Burk–Thompson and bulk aerodynamic boundary layer schemes produced the strongest storms while the MRF scheme produced the weakest storm. Simulated horizontal precipitation structures varied substantially between the different PBL schemes, suggesting that accurate forecasts of precipitation in hurricanes can be just as sensitive to the formulation of the PBL as they are to the cloud microphysical parameterizations. Each PBL scheme is different in its formulation of the vertical mixing within the PBL and the surface fluxes, with the exception of the MRF and Blackadar schemes, which share essentially the same surface flux parameterization. Detailed analyses of the PBL schemes describe the key differences in the surface fluxes and how they impact storm intensity. In order to isolate the effects of vertical mixing and surfaces fluxes, simulations were conducted in which each of the surface flux schemes was used in conjunction with the same vertical mixing scheme, and vice versa. These experiments indicate that simulated intensity is largely determined by the surface fluxes rather than by the vertical mixing, with the exception of the MRF PBL case, in which excessively deep vertical mixing acts to dry the lower PBL and reduce hurricane intensity. Simulations that vary only the surface fluxes suggest that the intensity of the simulated hurricane increases with increasing values of the ratio of the exchange coefficients for enthalpy and momentum, Ck/CD. However, even for identical values of Ck/CD, the simulated intensity varies depending on the wind speed dependence of the surface roughness parameter z0.