The main aim of the paper is to investigate the transitions of the thermohaline circulation in a spherical shell in a parameter regime which only allows transitions to multiple equilibria. We find that the first transition is either continuous (Type-I) or drastic (Type-II) depending on the sign of the transition number. The transition number depends on the system parameters and l(c), which is the common degree of spherical harmonics of the first critical eigenmodes, and it can be written as a sum of terms describing the nonlinear interactions of various modes with the critical modes. We obtain the exact formulas of this transition number for l(c) = 1 and l(c) = 2 cases. Numerically, we find that the main contribution to the transition number is due to nonlinear interactions with modes having zero wave number and the contribution from the nonlinear interactions with higher frequency modes is negligible. In our numerical experiments we encountered both types of transition for Le < 1 but only continuous transition for Le > 1. In the continuous transition scenario, we rigorously prove that an attractor in the phase space bifurcates which is homeomorphic to the 2l(c) dimensional sphere and consists entirely of degenerate steady state solutions.