Preprints
https://doi.org/10.5194/gmd-2023-39
https://doi.org/10.5194/gmd-2023-39
Submitted as: model evaluation paper
 | 
15 Mar 2023
Submitted as: model evaluation paper |  | 15 Mar 2023
Status: a revised version of this preprint was accepted for the journal GMD and is expected to appear here in due course.

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

Katharina Martha Gallmeier, J. Xavier Prochaska, Peter Cornillon, Dimitris Menemenlis, and Madolyn Kelm

Abstract. We have assembled 2,851,702 nearly cloud-free cutout images (sized 144×144 km2) 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 LCC4320 global-ocean simulation for a one-year period starting on November 17, 2011 but otherwise matched in geography and day-of-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-to-80-km scales (i.e., submesoscale-to-mesoscale). A principal finding is that the LLC4320 simulation reproduces well, over a large fraction of the ocean, the observed distribution of SST patterns, 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 the vicinity of the point at which western boundary currents separate from the continental margin, (2) in the Antarctic Circumpolar Current (ACC), especially in the eastern half of the Indian Ocean, and (3) in an equatorial band equatorward of 15°. It is clear that (1) is a result of premature separation in the simulated western boundary currents. The model output in (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 (3), 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/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.

Katharina Martha Gallmeier et al.

Status: closed

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • CEC1: 'Comment on gmd-2023-39', Astrid Kerkweg, 20 Apr 2023
    • AC1: 'Reply on CEC1', Katharina Gallmeier, 25 Apr 2023
  • RC1: 'Review on gmd-2023-39', Takaya Uchida, 24 Apr 2023
    • AC2: 'Reply on RC1', Katharina Gallmeier, 17 May 2023
    • AC4: 'Reply on RC1', Katharina Gallmeier, 18 May 2023
  • RC2: 'Comment on gmd-2023-39', Anonymous Referee #2, 26 Apr 2023
    • AC3: 'Reply on RC2', Katharina Gallmeier, 17 May 2023

Status: closed

Comment types: AC – author | RC – referee | CC – community | EC – editor | CEC – chief editor | : Report abuse
  • CEC1: 'Comment on gmd-2023-39', Astrid Kerkweg, 20 Apr 2023
    • AC1: 'Reply on CEC1', Katharina Gallmeier, 25 Apr 2023
  • RC1: 'Review on gmd-2023-39', Takaya Uchida, 24 Apr 2023
    • AC2: 'Reply on RC1', Katharina Gallmeier, 17 May 2023
    • AC4: 'Reply on RC1', Katharina Gallmeier, 18 May 2023
  • RC2: 'Comment on gmd-2023-39', Anonymous Referee #2, 26 Apr 2023
    • AC3: 'Reply on RC2', Katharina Gallmeier, 17 May 2023

Katharina Martha Gallmeier et al.

Model code and software

Python Code J. Xavier Prochaska, Katharina Gallmeier, and Madolyn Kelm https://doi.org/10.5281/zenodo.7545904

Katharina Martha Gallmeier et al.

Viewed

Total article views: 716 (including HTML, PDF, and XML)
HTML PDF XML Total BibTeX EndNote
505 190 21 716 6 3
  • HTML: 505
  • PDF: 190
  • XML: 21
  • Total: 716
  • BibTeX: 6
  • EndNote: 3
Views and downloads (calculated since 15 Mar 2023)
Cumulative views and downloads (calculated since 15 Mar 2023)

Viewed (geographical distribution)

Total article views: 695 (including HTML, PDF, and XML) Thereof 695 with geography defined and 0 with unknown origin.
Country # Views %
  • 1
1
 
 
 
 
Latest update: 15 Sep 2023
Download
Short summary
This paper introduces an approach to evaluate the performance of numerical models of ocean circulation. We compare the structure of satellite-derived sea surface temperature anomalies (SSTa) determined by a machine learning algorithm at 10-80 km scales to those output by a high resolution MITgcm run. The simulation over much of the ocean reproduces the observed distribution of SSTa patterns well. This general agreement, alongside few notable exceptions, highlights the potential of this approach.