A Data-Driven Approach for Parameterizing Submesoscale Vertical Buoyancy Fluxes in the Ocean Mixed Layer

TL;DR

This paper introduces a data-driven CNN-based method to parameterize submesoscale vertical buoyancy fluxes in ocean models, improving accuracy over traditional physics-based approaches by leveraging high-resolution simulation data.

Contribution

The study presents a novel CNN-based parameterization for submesoscale buoyancy fluxes, capturing complex interactions and dependencies missed by existing physics-based models.

Findings

CNN outperforms baseline parameterizations in offline skill across regions and seasons.

Prediction accuracy decreases during active submesoscales in winter and spring.

Strong dependency on mixed layer depth and large scale strain field was observed.

Abstract

Parameterizations of O(1-10)km submesoscale flows in General Circulation Models (GCMs) represent the effects of unresolved vertical buoyancy fluxes in the ocean mixed layer. These submesoscale flows interact non-linearly with mesoscale and boundary layer turbulence, and it is challenging to account for all the relevant processes in physics-based parameterizations. In this work, we present a data-driven approach for the submesoscale parameterization, that relies on a Convolutional Neural Network (CNN) trained to predict mixed layer vertical buoyancy fluxes as a function of relevant large-scale variables. The data used for training is given from 12 regions sampled from the global high-resolution MITgcm-LLC4320 simulation. When compared with the baseline of a submesoscale physics-based parameterization, the CNN demonstrates high offline skill across all regions, seasons, and filter scales tested in this study. During seasons when submesoscales are most active, which generally corresponds to winter and spring months, we find that the CNN prediction skill tends to be lower than in summer months. The CNN exhibits strong dependency on the mixed layer depth and on the large scale strain field, a variable closely related to frontogenesis, which is currently missing from the submesoscale parameterizations in GCMs.