Studies of ocean surface transport often invoke the “Eulerian-mean hypothesis”: that wave-agnostic general circulation models neglecting explicit surface waves effects simulate the Eulerian-mean ocean velocity time-averaged over surface wave oscillations. Acceptance of the Eulerian-mean hypothesis motivates reconstructing the total, Lagrangian-mean surface velocity by adding Stokes drift to model output. Here, we show that the Eulerian-mean hypothesis is inconsistent, because wave-agnostic models cannot accurately simulate the Eulerian-mean velocity if Stokes drift is significant compared to the Eulerian-mean or Lagrangian-mean velocity. We conclude that Stokes drift should not be added to ocean general circulation model output. We additionally show the viability of the alternative “Lagrangian-mean hypothesis” using a theoretical argument and by comparing a wave-agnostic global ocean simulation with an explicitly wave-averaged simulation. We find that our wave-agnostic model accurately simulates the Lagrangian-mean velocity even though the Stokes drift is significant.

A document by Simone Silvestri. Click on the document to view its contents.

Current eddy-permitting and eddy-resolving ocean models require dissipation to prevent a spurious accumulation of enstrophy at the grid scale. We introduce a new numerical scheme for momentum advection in large-scale ocean models that involves upwinding through a weighted essentially non-oscillatory (WENO) reconstruction. The new scheme provides implicit dissipation and thereby avoids the need for an additional explicit dissipation that may require calibration of unknown parameters. This approach uses the rotational, "vector invariant" formulation of the momentum advection operator that is widely employed by global general circulation models. A novel formulation of the WENO "smoothness indicators" is key for avoiding excessive numerical dissipation of kinetic energy and enstrophy at grid-resolved scales. We test the new advection scheme against a standard approach that combines explicit dissipation with a dispersive discretization of the rotational advection operator in two scenarios: (i) two-dimensional turbulence and (ii) three-dimensional baroclinic equilibration. In both cases, the solutions are stable, free from dispersive artifacts, and achieve increased "effective" resolution compared to other approaches commonly used in ocean models.

We describe CATKE, a parameterization for ocean microturbulence with scales between 1 and 100 meters. CATKE is a one-equation model that predicts diffusive turbulent vertical fluxes a prognostic turbulent kinetic energy (TKE) and a diagnostic mixing length that features a dynamic model for convective adjustment (CA). With its convective mixing length, CATKE predicts not just the depth range where microturbulence acts but also the timescale over which mixing occurs, an important aspect of turbulent convection not captured by convective adjustment schemes. As a result, CATKE can describe the competition between convection and other processes such as baroclinic restractification or biogeochemical production-destruction. We estimate CATKE’s free parameters with a posteriori calibration to eighteen large eddy simulations of the ocean surface boundary layer, and validate CATKE against twelve additional large eddy simulations with stronger and weaker forcing than used during calibration. We find that a CATKE-parameterized single column model accurately predicts the depth structure of buoyancy and momentum at vertical resolutions between 2 and 16 meters and with time steps of 10-20 minutes. We propose directions for future model development, and future efforts to recalibrate CATKE’s parameters against more comprehensive and realistic datasets.

We describe CATKE, a parameterization for fluxes associated with small-scale or “microscale’ ocean turbulent mixing on scales between 1 and 100 meters. CATKE uses a downgradient formulation that depends on a prognostic turbulent kinetic energy (TKE) variable and a diagnostic mixing length scale that includes a dynamic convective adjustment (CA) component. With its dynamic convective mixing length, CATKE predicts not just the depth spanned by convective plumes but also the characteristic convective mixing timescale, an important aspect of turbulent convection not captured by simpler static convective adjustment schemes. As a result, CATKE can describe the competition between convection and other processes such as shear-driven mixing and baroclinic restratification. To calibrate CATKE, we use Ensemble Kalman Inversion to minimize the error between 21~large eddy simulations (LES) and predictions of the LES data by CATKE-parameterized single column simulations at three different vertical resolutions. We find that CATKE makes accurate predictions of both idealized and realistic LES compared to microscale turbulence parameterizations commonly used in climate models.

Parameterizations of unresolved turbulent processes often compromise the fidelity of large-scale ocean models. In this work, we argue for a Bayesian approach to the refinement and evaluation of turbulence parameterizations. Using an ensemble of large eddy simulations of turbulent penetrative convection in the surface boundary layer, we demonstrate the method by estimating the uncertainty of parameters in the convective limit of the popular ‘K-Profile Parameterization’. We uncover structural deficiencies and propose an alternative scaling that overcomes them.

Between 5 and 25% of the total momentum transferred between the atmosphere and ocean is transmitted via the growth of long surface gravity waves called ‘swell’. In this paper, we use large eddy simulations to show that swell-transmitted momentum excites near-inertial waves and drives turbulent mixing that deepens a rotating, stratified, turbulent ocean surface boundary layer. We find that swell-transmitted currents are less effective at producing turbulence and mixing the boundary layer than currents driven by an effective surface stress. Overall, however, the differences between swell-driven and surface-stress-driven boundary layers are relatively minor. In consequence, our results corroborate assumptions made in Earth system models that neglect the vertical structure of swell-transmitted momentum fluxes and instead parameterize all air-sea momentum transfer processes with an effective surface stress.

We explore how neural differential equations (NDEs) may be trained on highly resolved fluid-dynamical models of unresolved scales providing an ideal framework for data-driven parameterizations in climate models. NDEs overcome some of the limitations of traditional neural networks (NNs) in fluid dynamical applications in that they can readily incorporate conservation laws and boundary conditions and are stable when integrated over time. We advocate a method that employs a ‘residual’ approach, in which the NN is used to improve upon an existing parameterization through the representation of residual fluxes which are not captured by the base parameterization. This reduces the amount of training required and providing a method for capturing up-gradient and nonlocal fluxes. As an illustrative example, we consider the parameterization of free convection of the oceanic boundary layer triggered by buoyancy loss at the surface. We demonstrate that a simple parameterization of the process — convective adjustment — can be improved upon by training a NDE against highly resolved explicit models, to capture entrainment fluxes at the base of the well-mixed layer, fluxes that convective adjustment itself cannot represent. The augmented parameterization outperforms existing commonly used parameterizations such as the K-Profile Parameterization (KPP). We showcase that the NDE performs well independent of the time-stepper and that an online training approach using differentiable simulation via the Julia scientific machine learning software stack improves accuracy by an order-of-magnitude. We conclude that NDEs provide an exciting route forward to the development of representations of sub-grid-scale processes for climate science, opening up myriad new opportunities.

We argue that ocean general circulation models and observations based on Ekman or geostrophic balance provide estimates of the Lagrangian-mean ocean velocity field averaged over surface waves — the total time-averaged velocity that advects oceanic tracers, particles, and water parcels. This interpretation contradicts an assumption often made in ocean transport studies that numerical models and observations based on dynamical balances estimate the Eulerian-mean velocity — the velocity time-averaged at a fixed position and only _part_ of the total ocean velocity. Our argument uses the similarity between the wave-averaged Lagrangian-mean momentum equations appropriate at large oceanic scales, and the momentum equations solved by “wave-agnostic” general circulation models that neglect surface wave effects. We further our case by comparing a realistic, global, “wave-agnostic” general circulation ocean model to a wave-averaged Lagrangian-mean general circulation ocean model at eddy-permitting 1/4-degree resolution, and find that the wave-agnostic velocity field is almost identical to the wave-averaged Lagrangian-mean velocity.