An evaluation of the LLC4320 global-ocean simulation based on the submesoscale structure of modeled sea surface temperature fields

Abstract

Abstract. We have assembled 2 851 702 nearly cloud-free cutout images (sized 144 km × 144 km) of sea surface temperature (SST) data from the entire 2012–2020 Level-2 Visible Infrared Imaging Radiometer Suite (VIIRS) dataset to perform a quantitative comparison to the ocean model output from the MIT General Circulation Model (MITgcm). Specifically, we evaluate outputs from the LLC4320 (LLC, latitude–longitude–polar cap) 148∘ global-ocean simulation for a 1-year period starting on 17 November 2011 but otherwise matched in geography and the day of the year to the VIIRS observations. In lieu of simple (e.g., mean, standard deviation) or complex (e.g., power spectrum) statistics, we analyze the cutouts of SST anomalies with an unsupervised probabilistic autoencoder (PAE) trained to learn the distribution of structures in SST anomaly (SSTa) on ∼ 10–80 km scales (i.e., submesoscale to mesoscale). A principal finding is that the LLC4320 simulation reproduces, over a large fraction of the ocean, the observed distribution of SSTa patterns well, both globally and regionally. Globally, the medians of the structure distributions match to within 2σ for 65 % of the ocean, despite a modest, latitude-dependent offset. Regionally, the model outputs reproduce mesoscale variations in SSTa patterns revealed by the PAE in the VIIRS data, including subtle features imprinted by variations in bathymetry. We also identify significant differences in the distribution of SSTa patterns in several regions: (1) in an equatorial band equatorward of 15∘; (2) in the Antarctic Circumpolar Current (ACC), especially in the eastern half of the Indian Ocean; and (3) in the vicinity of the point at which western boundary currents separate from the continental margin. It is clear that region 3 is a result of premature separation in the simulated western boundary currents. The model output in region 2, the southern Indian Ocean, tends to predict more structure than observed, perhaps arising from a misrepresentation of the mixed layer or of energy dissipation and stirring in the simulation. The differences in region 1, the equatorial band, are also likely due to model errors, perhaps arising from the shortness of the simulation or from the lack of high-frequency and/or wavenumber atmospheric forcing. Although we do not yet know the exact causes for these model–data SSTa differences, we expect that this type of comparison will help guide future developments of high-resolution global-ocean simulations.