In this study, it is shown (i) that, as a result of the nonlinearity of the seawater equation of state, unresolved scales represent a major source of uncertainties in the computation of the large-scale horizontal density gradient from the large-scale temperature and salinity fields, and (ii) that the effect of these uncertainties can be simulated using random processes to represent unresolved temperature and salinity fluctuations. The results of experiments performed with a low resolution global ocean model show that this parameterization has a considerable effect on the average large-scale circulation of the ocean, especially in the regions of intense mesoscale activity. The large-scale flow is less geostrophic, with more intense associated vertical velocities, and the average geographical position of the main temperature and salinity fronts is more consistent with observations. In particular, the simulations suggest that the stochastic effect of the unresolved temperature and salinity fluctuations on the large-scale density field may be sufficient to explain why the Gulf Stream pathway systematically overshoots in non-stochastic low resolution ocean models.