This paper investigates the sensitivity of sea breeze (SB) simulations to combinations of boundary-layer turbulence and land-surface process parameterizations implemented in the MM5 mesoscale meteorological mode for an observed SB case over the Swedish west coast. Various combinations from four different planetary boundary layer (PBL) schemes [Blackadar, Gayno-Seaman (GS), Eta, MRF], and two land surface model (LSM) schemes (SLAB, Noah) with different complexity are designed to simulate a typical SB case over the Swedish west coast. The simulations are conducted using two-way interactively nested grids. Simulated 10-m winds are compared against observed near-surface wind data from the GOTE2001 campaign to examine the diurnal cycle of wind direction and speed for SB timing. The SB (vertical) circulation is also compared in the different experiments. The results show that the different combinations of PBL and LSM parameterization schemes result in different SB timing and vertical circulation characteristics. All experiments predict a delayed SB. The vertical component of the SB circulation varies in the experiments, among which the GS PBL scheme produces the strongest SB circulation. Evident differences between the SLAB and Noah LSMs are also found, especially in maximum of updraft and downdraft velocities of the SB vertical circulation. The results have significant implications for convective initiation, air quality studies and other environmental problems in coastal areas.