Geostrophic turbulent eddies play a crucial role in the oceans, mixing properties such as heat, salt, and geochemical tracers. A useful reduced model for geostrophic turbulence is barotropic (2D) turbulence. The focus of this study is on 2D β -plane turbulence with quadratic drag, which, although arguably the most realistic barotropic model of ocean turbulence, has remained unexplored thus far. We first review and test classical scaling arguments for the eddy diffusivity in three regimes: the strong friction limit, the weak friction/strong β limit, and a transition regime. We then develop a generalized theory by parameterizing the nonlinear eddy–eddy interactions as a stochastic process, which leads to an analytical solution for the eddy diffusivity spectrum, whose integral yields a “bulk” diffusivity. The theory successfully predicts the smooth transition of diffusivity across the three regimes, and echoes with the recent argument that eddy phase propagation relative to the mean flow suppresses the eddy diffusivity. Moreover, the generalized theory reduces to the classical scaling arguments in both the strong friction and strong β limits, which has not been clear from the previous work on diffusivity suppression by flow-relative phase propagation.