Multidecadal variability of the overturning circulation in presence of eddy turbulence

At low-resolution, idealized ocean circulation models forced by prescribed differential surface heat fluxes show spontaneous multidecadal variability depending critically on eddy diffusivity coefficients. The existence of this critical threshold in the range of observational estimates legitimates some doubt on the relevance of such intrinsic oscillations in the real ocean. Through a series of numerical simulations with increasing resolution up to eddy-resolving ones (10 km) and various diapycnal diffusivities, this multidecadal variability proves a generic ubiquitous feature, at least in model versions with a flat bottom. The mean circulation largely changes in the process of refining the horizontal grid (along with the associated implicit viscosity and diffusivity), and the spatial structure of the variability is largely modified, but there is no clear influence of the resolution on the main oscillation period. The interdecadal variability appears even more robust to low vertical diffusivity and overturning when mesoscale eddies are resolved. The mechanism previously proposed for these oscillations, involving westward-propagating baroclinically unstable Rossby waves in the subpolar region and its feedback on the mean circulation, appears unaffected by mesoscale turbulence and is simply displaced following the polar front