Proper orthogonal decomposition reduced-order model of the global oceans

Abstract A reduced-order model (ROM) of the global oceans is developed by projecting the hydrostatic Boussinesq equations of motion onto a proper orthogonal decomposition (POD) basis. Three-dimensional POD modes are calculated from the ocean fields of an ensemble climate reanalysis dataset. The coefficients in the POD ROM are calculated using a regression approach. The performance of various POD ROM configurations are assessed. Each configuration is derived from an alternate sea-water equation of state, linking the density and temperature fields. POD ROM variants incorporating an equation of state in which density is a quadratic function of temperature, are able to reproduce the statistics of the large-scale structures at a fraction of the computational cost required to numerically simulate this flow. Due to the speed and efficiency of calculation, such reduced-order models of the global geophysical system will enable researchers and policy makers to assess the physical risk for a broader range of potential future climate scenarios.