Impact of global ocean model resolution on sea-level variability with emphasis on interannual time scales

Abstract. Four global ocean/sea-ice simulations driven by the same realistic 47-year daily atmospheric forcing were performed by the DRAKKAR group at 2°, 1°, ½°, and ¼° resolutions. Simulated mean sea-surface heights (MSSH) and sea-level anomalies (SLA) are collocated over the period 1993–2004 onto the AVISO dataset. MSSH fields are compared with an inverse estimate. SLA datasets are filtered and compared over various time and space scales with AVISO regarding three characteristics: SLA standard deviations, spatial correlations between SLA variability maps, and temporal correlations between observed and simulated band-passed filtered local SLA timeseries. Beyond the 2°−1° transition whose benefits are moderate, further increases in resolution and associated changes in subgrid scale parameterizations simultaneously induce (i) strong increases in SLA standard deviations, (ii) strong improvements in the spatial distribution of SLA variability, and (iii) slight decreases in temporal correlations between observed and simulation SLA timeseries. These 3 effects are not only clear on mesoscale (14–180 days) and quasi-annual (5–18 months) fluctuations, but also on the slower (interannual), large-scale variability ultimately involved in ocean-atmosphere coupled processes. Most SLA characteristics are monotonically affected by successive resolution increases, but irregularly and with a strong dependance on frequency and latitude. Benefits of enhanced resolution are greatest in the 1°−½° and ½°−¼° transitions, in the 14–180 day range, and within eddy-active mid- and high-latitude regions. In the real ocean, most eddy-active areas are characterized by a strong SLA variability at all timescales considered here; this localized, broad-banded temporal variability is only captured at ¼° resolution.