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 ]. the orbital variations, on the other [13, Chap. 15]. The next step would be to verify the existence of causal relations between the orbital forcing and the climatic system's behavior.Equilibrium models of the climatic system based on radiation balance, so-called energy-balance models ([5, pp. 461-480], [22]) are now well understood. Insolation intensity appears in them as a bifurcation parameter, and response to its changes can be assumed to be quasi-static. Changes in insolation of the magnitude given by orbital variations, when imposed upon energy-balance models, fail to produce temperature changes larger than fractions of a degree Kelvin [22,.