Most of the ocean’s kinetic energy is contained within the mesoscale
eddy field. Models that do not resolve these eddies tend to parameterize
their impacts such that the parameterized transport of buoyancy and
tracers reduces the large-scale available potential energy and spreads
tracers. However, the parameterizations used in the ocean components of
current generation Earth System Models (ESMs) rely on an assumption of a
flat ocean floor even though observations and high-resolution modelling
show that eddy transport is sensitive to the potential vorticity
gradients associated with a sloping seafloor. We show that buoyancy
transport coefficient diagnosed from idealized eddy-resolving
simulations is indeed reduced over both prograde and retrograde bottom
slopes (topographic wave propagation along or against the mean flow,
respectively) and that the reduction can be skilfully captured by a
mixing length parameterization by introducing the topographic Rhines
scale as a length scale. This modified ‘GM’ parameterization enhances
the strength of thermal wind currents over the slopes in
coarse-resolution, non-eddying, simulations. We find that in realistic
global coarse-resolution simulations the impact of topography is most
pronounced at high latitudes, enhancing the mean flow strength and
reducing temperature and salinity biases. Reducing the buoyancy
transport coefficient further with a mean-flow dependent eddy efficiency
factor, has notable effects also at lower latitudes and leads to
reduction of global mean biases.