One of the most important spatial scales in the ocean circulation is the internal Rossby radius of deformation L D ; it ranges from 50 to 100 km at mid-latitudes to a few km in the polar regions (Hallberg, 2013). At this scale, perturbations are amplified on mean flows through mixed barotropic/baroclinic instability, giving rise to ocean eddies. Interactions between these eddies and the mean flow can lead to up-gradient momentum transport affecting the strength and separation of ocean western boundary currents such as the Kuroshio and Agulhas (Chassignet et al., 2020).Most climate models, in particularly those used in CMIP5 and CMIP6, do not resolve ocean processes at the scale L D as the spatial grid size used is too large; typically 1° (Eyring et al., 2016). The main reason is computational costs, as doubling the horizontal resolution increases these costs roughly by a factor 10. Effects of subgrid-scale processes are hence parameterized in these models. For example, the effect of ocean eddies on tracer transport is represented by the Gent- McWilliams (Gent et al., 1995) scheme, but such a scheme cannot capture, for example, the up-gradient momentum transport. Hence, western boundary flows are too weak and diffuse, and do not separate at the correct location (Chassignet et al., 2020).Over the last few years, first simulations have been performed with global climate models, where the ocean model component has a resolution of 0.1°, which is smaller than L D for many locations on the globe (Chang Abstract We introduce a "symbiotic" ocean modeling strategy that leverages data-driven and machine learning methods to allow high-and low-resolution dynamical models to mutually benefit from each other. In this work we mainly focus on how a low-resolution model can be enhanced within a symbiotic model configuration. The broader aim is to enhance the representation of unresolved processes in low-resolution models, while simultaneously improving the efficiency of high-resolution models. To achieve this, we use a grid-switching approach together with hybrid modeling techniques that combine linear regression-based methods with nonlinear echo state networks. The approach is applied to both the Kuramoto-Sivashinsky equation and a single-layer quasi-geostrophic ocean model, and shown to simulate short-term and long-term behavior better than either purely data-based methods or low-resolution models. By maintaining key flow characteristics, the hybrid modeling techniques are also able to provide higher quality initial conditions for high-resolution models, thereby improving their efficiency.Plain Language Summary Models of the ocean vary in complexity. Some are very detailed and manage to show oceanic vortices, whereas others are very efficient but coarse, and unable to compute such vortices. The idea in this paper is to let these different model types work together and benefit from each other, as if in a symbiosis. With knowledge of differences between the detailed and coarse model we can use machine learning techniques to improv...