A surface-aware projection basis for quasigeostrophic flow
Recent studies indicate that altimetric observations of the ocean's mesoscale eddy field reflect the combined influence of surface buoyancy and interior potential vorticity anomalies. The former have a surface-trapped structure, while the latter have a more grave form. To assess the relative importance of each contribution to the signal, it is useful to project the observed field onto a set of modes that separates their influence in a natural way. However, the surface-trapped dynamics are not well-represented by standard baroclinic modes; moreover, they are dependent on horizontal scale. Here we derive a modal decomposition that results from the simultaneous diagonalization of the energy and a generalisation of potential enstrophy that includes contributions from the surface buoyancy fields. This approach yields a family of orthonomal bases that depend on two parameters: the standard baroclinic modes are recovered in a limiting case, while other choices provide modes that represent surface and interior dynamics in an efficient way. For constant stratification, these modes consist of symmetric and antisymmetric exponential modes that capture the surface dynamics, and a series of oscillating modes that represent the interior dynamics. Motivated by the ocean, where shears are concentrated near the upper surface, we also consider the special case of a quiescent lower surface. In this case, the interior modes are independent of wavenumber, and there is a single exponential surface mode that replaces the barotropic mode. We demonstrate the use and effectiveness of these modes by projecting the energy in a set of simulations of baroclinic turbulence.