<jats:p>The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice, and plays a key role in climate variability. Numerical models used in climate studies should therefore have a good representation of the mixed layer, especially its depth (MLD). Here we use simulations from the Ocean Model Intercomparison Project (OMIP), which have been forced by a common atmospheric state, to assess the realism of the simulated MLDs. For model validation, an updated MLD dataset has been computed from observations using the fixed density threshold recommended by the OMIP protocol. We evaluate the influence of horizontal resolution by using six pairs of simulations, non-eddying (typically 1&amp;#176; resolution) and eddy-rich (1/10&amp;#176; to 1/16&amp;#176; resolution). In winter, low resolution models exhibit large biases in the deep water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode water formation regions of the northern hemisphere, where the eddy-rich models produce a more robust MLD and deep biases are reduced. The Southern Ocean offers a more contrasted view, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability of MLD in winter. Contrary to an hypothesis that the deepening of the MLD in anticyclones would make the MLD deeper globally, eddy-rich models tend to have a shallower MLD in the zonal mean. In summer, a deep MLD bias is found in all the non-eddying models north of the equator; this bias is greatly reduced at high resolution. In addition, our study highlights the sensitivity of the MLD computation to choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects.</jats:p>