Mixed mode oscillations in a conceptual climate model

Much work has been done on relaxation oscillations and other simple oscillators in conceptual climate models. However, the oscillatory patterns in climate data are often more complicated than what can be described by such mechanisms. This paper examines complex oscillatory behavior in climate data through the lens of mixed-mode oscillations. As a case study, a conceptual climate model with governing equations for global mean temperature, atmospheric carbon, and oceanic carbon is analyzed. The nondimensionalized model is a fast/slow system with one fast variable (corresponding to ice volume) and two slow variables (corresponding to the two carbon stores). Geometric singular perturbation theory is used to demonstrate the existence of a folded node singularity. A parameter regime is found in which (singular) trajectories that pass through the folded node are returned to the singular funnel in the limiting case where $\epsilon = 0$. In this parameter regime, the model has a stable periodic orbit of type $1^s$ for some $s>0$. To our knowledge, it is the first conceptual climate model demonstrated to have the capability to produce an MMO pattern.