A low-order quasigeostrophic double-gyre ocean model is subjected to an aperiodic forcing that mimics time dependence dominated by interdecadal variability. This model is used as a prototype of an unstable and nonlinear dynamical system with time-dependent forcing to explore basic features of climate change in the presence of natural variability. The study relies on the theoretical framework of nonautonomous dynamical systems and of their pullback attractors (PBAs), that is, of the time-dependent invariant sets attracting all trajectories initialized in the remote past. The existence of a global PBA is rigorously demonstrated for this weakly dissipative nonlinear model. Ensemble simulations are carried out and the convergence to PBAs is assessed by computing the probability density function (PDF) of localization of the trajectories. A sensitivity analysis with respect to forcing amplitude shows that the PBAs experience large modifications if the underlying autonomous system is dominated by small-amplitude limit cycles, while less dramatic changes occur in a regime characterized by large-amplitude relaxation oscillations. The dependence of the attracting sets on the choice of the ensemble of initial states is then analyzed. Two types of basins of attraction coexist for certain parameter ranges; they contain chaotic and nonchaotic trajectories, respectively. The statistics of the former does not depend on the initial states whereas the trajectories in the latter converge to small portions of the global PBA. This complex scenario requires separate PDFs for chaotic and nonchaotic trajectories. General implications for climate predictability are finally discussed.

# Ocean & coupled ocean

Exploring the pullback attractors of a low-order quasigeostrophic ocean model: The deterministic case. Journal of Climate. 2016;29 (11) :4185-4202.Abstract

.
Low-frequency variability of the large-scale ocean circulation: a dynamical systems approach. Reviews of Geophysics. 2005;43.Abstract

.
Spatio-temporal variability in a mid-latitude ocean basin subject to periodic wind forcing. Atmosphere-ocean. 2007;45 (4) :227–250.Abstract

.
Data assimilation in meteorology and oceanography. Advances in Geophysics. 1991;33 :141–266.

.
Dynamical origin of low-frequency variability in a highly nonlinear midlatitude coupled model. Journal of Climate. 2006;19 (24).Abstract

.
Topics in Geophysical Fluid Dynamics: Atmospheric Dynamics, Dynamo Theory and Climate Dynamics. Springer-Verlag, New York/Berlin; 1987 pp. 485.

.
A highly nonlinear coupled mode of decadal variability in a mid-latitude ocean–atmosphere model. Dynamics of Atmospheres and Oceans. 2007;43 (3) :123–150.

.
Multiple equilibria and oscillatory modes in a mid-latitude ocean-forced atmospheric model. Nonlinear Processes in Geophysics. 2012;19 (5) :479–499.

.
Trends, interdecadal and interannual oscillations in global sea-surface temperatures. Climate Dynamics. 1998;14 (7) :545–569.Abstract

.
El Niño Southern Oscillation and the annual cycle: Subharmonic frequency-locking and aperiodicity. Physica D. 1996;98 :442–465.

.
Baroclinic and barotropic aspects of the wind-driven ocean circulation. Physica D: Nonlinear Phenomena. 2002;167 (1) :1–35.Abstract

.
Midlatitude ocean-atmosphere interaction in an idealized coupled model. Climate Dynamics. 2002;19 (8) :693–711.

.
Atmospheric circulations induced by a midlatitude SST front: A GCM study. Journal of Climate. 2012;25 (6) :1847–1853.Abstract

.
Low-Frequency Variability in Shallow-Water Models of the Wind-Driven Ocean Circulation. Part I: Steady-State Solution. Journal of Physical Oceanography. 2003;33 (4).Abstract

.
A mechanistic model of mid-latitude decadal climate variability. Physica D. 2008;237 :584–599.

.