Decadal large-scale low-frequency variability of the ocean circulation due to its nonlinear dynamics remains a big challenge for theoretical understanding and practical ocean modeling. This paper presents a novel fully data driven approach that addresses this challenge. Proposed is non-Markovian low-order methodology with stochastic closure and use of mode decomposition by multichannel Singular Spectrum Analysis. The multilayer stochastic linear model is obtained from the coarse-grained eddy-resolving ocean model solution, and with high accuracy it reproduces the main statistical properties of the decadal variability. The proposed methodology does not depend on the governing fluid dynamics equations and geometry of the problem, and it can be extended to other ocean models and ultimately to the real data.
BackgroundMidlatitude atmosphere and ocean possess not only significant interannual variability but also several large-scale variability modes on decadal and interdecadal timescales: North Atlantic Oscillation, Atlantic Multidecadal Oscillation, and Pacific Decadal Oscillation. Oceanic evidence of these modes is found in the sea surface temperature [Hansen and Bezdek, 1996], hydrography [Qiu and Joyce, 1992], and altimetry [Qiu and Chen, 2005], and it is mostly expressed in structural changes of the eastward jet extensions and the adjacent recirculation zones of the main western boundary currents, such as the Gulf Stream and Kuroshio. Physical origins of these large-scale low-frequency variability (LFV) modes remain unclear, and it is not even known to what extent these origins are intrinsic atmospheric, intrinsic oceanic, or coupled oceanic-atmospheric. One of the reasons, why coupled ocean-atmosphere comprehensive general circulation models (GCMs) cannot distinguish between these scenarios, is inability of their oceanic model components to simulate mesoscale eddies and their impact on the large-scale currents accurately and for a long time. On the other hand, seasonally forced, eddy-resolving ocean-only GCMs exhibit significant intrinsic LFV variability of mesoscale activity [Penduff et al., 2011], intergyre heat transport [Hall et al., 2004], sea level height [Cabanes et al., 2006], western boundary currents [Taguchi et al., 2007], and meridional overturning circulation [Biastoch et al., 2008]. This suggests that not only significant part of the LFV is intrinsic to the ocean but also oceanic mesoscale eddies could play significant role in driving it.Present theoretical understanding of the intrinsic LFV of the ocean starts from the null hypothesis that the ocean integrates high-frequency atmospheric forcing and responds in terms of the red variability spectrum [Hasselman, 1976]. However, over the last 20 years it became evident that in the presence of stochastic, seasonal, or constant atmospheric forcing, the most significant intrinsic LFV of the ocean is not red but operates on interannual-to-interdecadal time scales. One idea is that this variability can be understood in terms of early bifurcations...