The multi-layer quasi-geostrophic model of the wind-driven ocean gyres is numerically investigated using a combination of long-time runs (200 years) needed for accurate statistics, spatial resolutions (grid interval of less than one kilometer) needed for accurate representation of mesoscale eddies, and large Reynolds number (Re > 104) needed for more realistic flow regimes. We gradually increased the Reynolds number by lowering the eddy viscosity and analysed the corresponding changes of the large-scale circulation, energetics and eddy fluxes, with the goal to understand how the nonlinear eddy dynamics affects the large-scale ocean circulation, as more and more degrees of freedom become dynamically available. Three- and six-layer configurations of the model are considered in order to understand effects of higher baroclinic modes. A parameter sensitivity study is also carried out to show that the explored flow regime is robust.

As Re increases, most properties of the flow show no signs of approaching an asymptote, and the following tendencies are found. The time-mean flow properties tend to an asymptote in the three-layer model but not in the six-layer one, suggesting that higher baroclinic modes are dynamically more active at larger Re. The eddy kinetic and potential energies grow faster in the six-layer case. The intensity of the eddy forcing (eddy flux divergence) increases with Re. The inter-gyre eddy potential vorticity flux is predominantly northward and up-gradient for all Re studied. A comparison of the three- and six-layer model solutions revealed an inhibitory influence of high baroclinic modes on the penetration length of the eastward jet extension of the western boundary currents and on the strength of the adjacent recirculation zones. In large-Re regimes, the population of eddies is mostly sustained by the eddy generation at the eastern end of the eastward jet rather than in its central section. Finally, by studying the numerical convergence of the solutions, we found the empirical dependency between the eddy viscosity and the required grid resolution: halving the viscosity requires halving the grid spacing.