Search

Abstract. This study evaluates the impact of increasing resolution on Arctic Ocean simulations using five pairs of matched low- and high-resolution models within the OMIP-2 (Ocean Model Intercomparison Project phase 2) framework. The primary objective is to assess whether a higher resolution can mitigate typical biases in low-resolution models and improve the representation of key climate-relevant variables. We reveal that increasing the horizontal resolution contributes to a reduction in biases in mean temperature and salinity and improves the simulation of the Atlantic water layer and its decadal warming events. A higher resolution also leads to better agreement with observed surface mixed-layer depth, cold halocline base depth and Arctic gateway transports in the Fram and Davis straits. However, the simulation of the mean state and temporal changes in Arctic freshwater content does not show improvement with increased resolution. Not all models achieve improvements for all analyzed ocean variables when spatial resolution is increased so it is crucial to recognize that model numerics and parameterizations also play an important role in faithful simulations. Overall, a higher resolution shows promise in improving the simulation of key Arctic Ocean features and processes, but efforts in model development are required to achieve more accurate representations across all climate-relevant variables.

Advances in coastal modeling and computation provide the opportunity to examine non-hydrostatic and compressible fluid effects at very small scales, but the cost of these new capabilities and the accuracy of these models versus trusted non-hydrostatic codes has yet to be determined. Here the Coastal and Regional Ocean COmmunity model (CROCO, v1.2) and the National Center for Atmospheric Research large-eddy simulation (NCAR-LES) model are compared, with a focus on their simulation accuracy and computational efficiency. These models differ significantly in numerics and capabilities, so they are run on common classic problems of surface-forced, boundary-layer turbulence. In terms of accuracy, we compare turbulence statistics, including the effect of the explicit subgrid-scale (SGS) parameterization, the effect of the second (dilatational) viscosity, and the sensitivity to the speed of sound, which is used as part of the CROCO compressible turbulence formulation. To gauge how far CROCO is from the NCAR-LES, we first compare the NCAR-LES with two other non-hydrostatic Boussinesq approximation LES codes (PALM and Oceananigans), defining the notion and magnitude of accuracy for the LES and CROCO comparison. To judge efficiency of CROCO, strong and weak scaling simulation sets vary different problem sizes and workloads per processor, respectively. Additionally, the effects of 2D decomposition of CROCO and NCAR-LES and supercomputer settings are tested. In summary, the accuracy comparison between CROCO and the NCAR-LES is similar to the NCAR-LES compared to other LES codes. However, the additional capabilities of CROCO (e.g., nesting, non-uniform grid, and realism of ocean configuration in general) and its weakly compressible formulation come with roughly an order of magnitude of additional costs, despite efforts to reduce them by adjusting the second viscosity and speed of sound as far as accuracy allows. However, a new variant of the non-hydrostatic CROCO formulation is cu

Mesoscale and submesoscale processes have crucial impacts on ocean biogeochemistry, importantly enhancing the primary production in nutrient-deficient ocean regions. Yet, the intricate biophysical interplay still holds mysteries. Using targeted high-resolution in situ observations in the South China Sea, we reveal that isopycnal submesoscale stirring serves as the primary driver of vertical nutrient transport to sustain the dome-shaped subsurface chlorophyll maximum (SCM) within a long-lived cyclonic mesoscale eddy. Density surface doming at the eddy core increased light exposure for phytoplankton production, while along-isopycnal submesoscale stirring disrupted the mesoscale coherence and drove significant vertical exchange of tracers. These physical processes play a crucial role in maintaining the elevated phytoplankton biomass in the eddy core. Our findings shed light on the universal mechanism of how mesoscale and submesoscale coupling enhances primary production in ocean cyclonic eddies, highlighting the pivotal role of submesoscale stirring in structuring marine ecosystems.

We study the impact of the Ekman currents and Stokes drift on the surface horizontal transport and mixing properties of the Mediterranean Sea. Finite Size Lyapunov Exponents (FSLE) at the ocean surface are computed over the whole basin using 25 years of satellite altimetry derived geostrophic currents, 10-m wind velocity and wave fields. We find that the transport pathways unveiled by the geostrophic Lagrangian Coherent Structures (LCS) are significantly modified by the ageostrophic currents (i.e., Ekman and Stokes induced velocities), often leading to a decrease of the retention capacity of the eddies. An exhaustive assessment of the regional and temporal variability of the FSLE shows an increase of the horizontal mixing processes, due to the ageostrophic component, up to 37% in regions such as the Gulf of Lion or the Aegean Sea, during the seasons where wind and waves are intense and persistent. Positive trends in the total FSLE (up to 1.2% of the value of FSLE per year in some regions) suggest that Mediterranean Sea has experienced a significant increase in the intensity of the mixing processes over the last decades. Ageostrophic features are considered to play a role in determining the properties of the relative dispersion. Through the analysis of the Lagrangian Anisotropy Index (LAI) using virtual and real pair of drifters, we observe that the particle dispersion is mainly dominated by the zonal flow, and that the ageostrophic currents induce meridional dispersion, particularly in regions where wind and wave are intensified.

© The Author(s), 2022. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Marques, G. M., Loose, N., Yankovsky, E., Steinberg, J. M., Chang, C.-Y., Bhamidipati, N., Adcroft, A., Fox-Kemper, B., Griffies, S. M., Hallberg, R. W., Jansen, M. F., Khatri, H., & Zanna, L. NeverWorld2: an idealized model hierarchy to investigate ocean mesoscale eddies across resolutions. Geoscientific Model Development, 15(17), (2022): 6567–6579, https://doi.org/10.5194/gmd-15-6567-2022., We describe an idealized primitive-equation model for studying mesoscale turbulence and leverage a hierarchy of grid resolutions to make eddy-resolving calculations on the finest grids more affordable. The model has intermediate complexity, incorporating basin-scale geometry with idealized Atlantic and Southern oceans and with non-uniform ocean depth to allow for mesoscale eddy interactions with topography. The model is perfectly adiabatic and spans the Equator and thus fills a gap between quasi-geostrophic models, which cannot span two hemispheres, and idealized general circulation models, which generally include diabatic processes and buoyancy forcing. We show that the model solution is approaching convergence in mean kinetic energy for the ocean mesoscale processes of interest and has a rich range of dynamics with circulation features that emerge only due to resolving mesoscale turbulence., This research has been supported by the US Department of Commerce (grant no. NA18OAR4320123), the Division of Ocean Sciences (grant nos. 1912420, 1912332, 1912357, 1912163, and 1912302), the Division of Atmospheric and Geospace Sciences (grant no. 1852977), and the Climate Program Office (grant nos. NA19OAR4310364, NA19OAR4310365, and NA19OAR4310366).

The ocean mixed layer is the interface between the ocean interior and the atmosphere or sea ice, and plays a key role in climate variability. It is thus critical that numerical models used in climate studies are capable of a good representation of the mixed layer, especially its depth. Here we evaluate the mixed layer depth (MLD) in six pairs of non-eddying (1° resolution) and eddy-rich (up to 1/16°) models from the Ocean Model Intercomparison Project (OMIP), forced by a common atmospheric state. For model validation, we use an updated MLD dataset computed from observations using the OMIP protocol (a constant density threshold). In winter, low resolution models exhibit large biases in the deep water formation regions. These biases are reduced in eddy-rich models but not uniformly across models and regions. The improvement is most noticeable in the mode water formation regions of the northern hemisphere. Results in the Southern Ocean are more contrasted, with biases of either sign remaining at high resolution. In eddy-rich models, mesoscale eddies control the spatial variability of MLD in winter. Contrary to a hypothesis that the deepening of the mixed layer in anticyclones would make the MLD larger globally, eddy-rich models tend to have a shallower mixed layer at most latitudes than coarser models do. In addition, our study highlights the sensitivity of the MLD computation to the choice of a reference level and the spatio-temporal sampling, which motivates new recommendations for MLD computation in future model intercomparison projects.