Statistical Inversion of Surface Ocean Kinematics from Sea Surface Temperature Observations

Abstract

AbstractThe sea surface temperature (SST) record provides a unique view of the surface ocean at high spatiotemporal resolution and holds useful information on the kinematics underlying the SST variability. To access this information, we develop a new local matrix inversion method that allows us to quantify the evolution of a given SST perturbation with a response function and to estimate velocity, diffusivity, and decay fields associated with it. The matrix inversion makes use of the stochastic climate model paradigm—we assume that SST variations are governed by a linear transport operator and a forcing that has a relatively short autocorrelation time scale compared to that of SST. We show that under these assumptions, the transport operator can be inverted from the covariance matrices of the underlying SST data. The accuracy of the results depends on the length of the time series, and in general the inverted properties depend on the spatial and time resolution of the SST data. Future studies could use the methodology to explore the interannual variability of SST anomalies; to estimate the scale dependency of ocean mixing; and to estimate anomaly propagation, both at the surface and in the interior. The methodology can be easily used with any gridded observations or model output with adequate time and spatial resolution, and it is not restricted to SST. The inversion code is written in Python and distributed as a MicroInverse package through GitHub and the Python Package Index.