“…The existing approaches may lead to numerical instabilities due to the presence of negative eddy viscosity or diffusivity. For example, the tracer diffusivity tensor produced by mesoscale eddies contains one negative eigenvalue (Bachman et al., 2020; Haigh & Berloff, 2021, 2022; Haigh et al., 2021; Kamenkovich et al., 2021; Lu et al., 2022; Ryzhov & Berloff, 2022). Parameterization of mesoscale momentum fluxes by the eddy viscosity operator is also challenging: the crucial physical process of kinetic energy backscatter can be captured only with negative eddy viscosity (Bachman, 2019; Jansen & Held, 2014; Jansen et al., 2019; Juricke et al., 2020).…”