Abstract. Aerosol-cloud interactions are explored using 1 km simulations of a case study of predominantly closedcell SE Pacific stratocumulus clouds. The simulations include realistic meteorology along with newly implemented cloud microphysics and sub-grid cloud schemes. The model was critically assessed against observations of liquid water path (LWP), broadband fluxes, cloud fraction (f c ), droplet number concentrations (N d ), thermodynamic profiles, and radar reflectivities.Aerosol loading sensitivity tests showed that at low aerosol loadings, changes to aerosol affected shortwave fluxes equally through changes to cloud macrophysical characteristics (LWP, f c ) and cloud albedo changes due solely to N d changes. However, at high aerosol loadings, only the N d albedo change was important. Evidence was also provided to show that a treatment of sub-grid clouds is as important as order of magnitude changes in aerosol loading for the accurate simulation of stratocumulus at this grid resolution.Overall, the control model demonstrated a credible ability to reproduce observations, suggesting that many of the important physical processes for accurately simulating these clouds are represented within the model and giving some confidence in the predictions of the model concerning stratocumulus and the impact of aerosol. For example, the control run was able to reproduce the shape and magnitude of the observed diurnal cycle of domain mean LWP to within ∼ 10 g m −2 for the nighttime, but with an overestimate for the daytime of up to 30 g m −2 . The latter was attributed to the uniform aerosol fields imposed on the model, which meant that the model failed to include the low-N d mode that was observed further offshore, preventing the LWP removal through precipitation that likely occurred in reality. The boundary layer was too low by around 260 m, which was attributed to the driving global model analysis. The shapes and sizes of the observed bands of clouds and open-cell-like regions of low areal cloud cover were qualitatively captured. The daytime f c frequency distribution was reproduced to within f c = 0.04 for f c >∼ 0.7 as was the domain mean nighttime f c (at a single time) to within f c = 0.02. Frequency distributions of shortwave top-of-the-atmosphere (TOA) fluxes from the satellite were well represented by the model, with only a slight underestimate of the mean by 15 %; this was attributed to near-shore aerosol concentrations that were too low for the particular times of the satellite overpasses. TOA long-wave flux distributions were close to those from the satellite with agreement of the mean value to within 0.4 %. From comparisons of N d distributions to those from the satellite, it was found that the N d mode from the model agreed with the higher of the two observed modes to within ∼ 15 %.