For a wide range of values of the incoming solar radiation, the Earth features at least two attracting states, which correspond to competing climates. The warm climate is analogous to the present one; the snowball climate features global glaciation and conditions that can hardly support life forms. Paleoclimatic evidences suggest that in the past our planet flipped between these two states. The main physical mechanism responsible for such an instability is the ice-albedo feedback. In a previous work, we defined the Melancholia states that sit between the two climates. Such states are embedded in the boundaries between the two basins of attraction and feature extensive glaciation down to relatively low latitudes. Here, we explore the global stability properties of the system by introducing random perturbations as modulations to the intensity of the incoming solar radiation. We observe noise-induced transitions between the competing basins of attraction. In the weak noise limit, large deviation laws define the invariant measure, the statistics of escape times, and typical escape paths called instantons. By constructing the instantons empirically, we show that the Melancholia states are the gateways for the noise-induced transitions. In the region of multistability, in the zero-noise limit, the measure is supported only on one of the competing attractors. For low (high) values of the solar irradiance, the limit measure is the snowball (warm) climate. The changeover between the two regimes corresponds to a first-order phase transition in the system. The framework we propose seems of general relevance for the study of complex multistable systems. Finally, we put forward a new method for constructing Melancholia states from direct numerical simulations, black which provides a possible alternative with respect to the edge-tracking algorithm.

The Earth is well-known to be, in the current astronomical configuration, in a regime where two asymptotic states can be realised. The warm state we live in is in competition with the ice-covered snowball state. The bistability exists as a result of the positive ice-albedo feedback. In a previous investigation performed on a intermediate complexity climate model we have identified the unstable climate states (Melancholia states) separating the co-existing climates, and studied their dynamical and geometrical properties. The Melancholia states are ice-covered up to the mid-latitudes and attract trajectories initialised on the basins boundary. In this paper, we study how stochastically perturbing the parameter controlling the intensity of the incoming solar radiation impacts the stability of the climate. We detect transitions between the warm and the snowball state and analyse in detail the properties of the noise-induced escapes from the corresponding basins of attraction. We determine the most probable paths for the transitions and find evidence that the Melancholia states act as gateways, similarly to saddle points in an energy landscape.