This paper outlines a pathway to understanding eddy diffusivity in stratified flows by considering scalar-velocity cospectra, also called scalar flux spectra (E<sub>w<em>θ</em></sub>). Its main objective is to unify existing scalings of eddy diffusivity into a quantitative recipe. We construct a synthetic model scalar-velocity cospectrum for stratified turbulence applicable to the ocean. Its integral over all wavenumbers, normalized by the mean scalar gradient, equals eddy diffusivity and scales principally with horizontal Froude number Fr<sub><em>h</em></sub> and buoyancy Reynolds number Re<sub>b</sub>. Our model spectrum reproduces published mixing coefficient estimates from oceanic microstructure and direct numerical simulation (DNS), as well as a new analysis of existing eddy covariance data, over the entire accessible range of Fr<sub><em>h</em></sub>, or equivalently, Ozmidov/Thorpe length scale ratios. However, scaling between eddy diffusivities and Fr<sub><em>h</em></sub> for strongly stratified flows, Fr<sub><em>h</em></sub> < 1, appears to differ between most DNS and oceanic measurements. We particularly emphasize the importance of the spectral slope <em>s<sub>b</sub></em> between the Ozmidov and integral length scales, finding best agreement with ocean observations when assuming <em>s<sub>b</sub></em> ≈ −4/3, but <em>s<sub>b</sub></em> ≈ 1/3 for DNS for the stratified wavenumber range. We argue that this difference may be due to the length scale of the forcing. From direct measurements of scalar-velocity cospectra in the ice-ocean boundary layer, we find that additional contributions to eddy diffusivity at low Froude numbers may result from background processes not scaling with Fr<em><sub>h</sub></em> or Re<sub>b</sub>, which warrants further investigation into boundary layer processes.