Climate - the "coarse-gridded" state of the coupled ocean - atmosphere system - varies on many time and space scales. The challenge is to relate such variation to specific mechanisms and to produce verifiable quantitative explanations. In this paper, we study the oceanic component of the climate system and, in particular, the different circulation regimes of the mid-latitude win driven ocean on the interannual time scale. These circulations are dominated by two counterrotating, basis scale gyres: subtropical and subpolar. Numerical techniques of bifurcation theory are used to stud the multiplicity and stability of the steady-state solution of a wind-driven, double-gyre, reduced-gravity, shallow water model. Branches of stationary solutions and their linear stability are calculated systematically as parameter are varied. This is one of the first geophysical studies i which such techniques are applied to a dynamical system with tens of thousands of degrees of freedom. Multiple stationary solutions obtain as a result of nonlinear interactions between the two main recirculating cell (cyclonic and anticyclonic) of the large- scale double-gyre flow. These equilibria appear for realistic values of the forcing and dissipation parameters. They undergo Hop bifurcation and transition to aperiodic solutions eventually occurs. The periodic and chaotic behaviour is probably related to an increased number of vorticity cells interaction with each other. A preliminary comparison with observations of the Gulf Stream and Kuroshio Extensions suggests that the intern variability of our simulated mid-latitude ocean is a important factor in the observed interannual variability o these two current systems.

The mathematical structure of a simple climate model is investigated. The model is governed by a system of two nonlinear, autonomous differential equations for the evolution in time of global temperature $T$ and meridional ice-sheet extent $L$. The system's solutions are studied by a combination of qualitative reasoning with explicit calculations, both analytical and numerical. For plausible values of the physical parameters, a branch of periodic solutions obtains, which is both orbitally and structurally stable. The amplitude of the stable periodic solutions in $T$ and $L$ correspond roughly to that obtained from proxy records of Quaternary glaciation cycles. The period of these solutions increases along the branch, until it becomes infinite, while the amplitude of the limiting solution is finite. The limiting solution is a homoclinic orbit formed by the reconnecting separatrix of a saddle. The exchange of stability between the branch of periodic solutions and the steady solution from which it arises is studied by a slight simplification of known methods [20], [21].