In the present paper, the geometrical nonlinear behaviour of the sandwich shallow spherical shells made of the functionally graded carbon nanotube-reinforced composites is assessed. In order to solve the proposed problem, the equilibrium, compatibility and constitutive equations are extended based on the first-order shear deformation theory. The obtained system of nonlinear differential equations is solved by employing an analytical approach based on the modified interactions method, which results a closed-form relationship between external load and deformation of the shell. Firstly, the obtained closed-form formula is validated by comparing the corresponding results with those obtained by other methods such as finite element method. After the validation procedure, the mentioned closed-form solution is employed to perform a comprehensive parametric study incorporating the mechanical and geometrical properties, such as the effects of CNTs agglomerating, depth and boundary conditions of the shell. The obtained results in the parametric studies show that all of the parameters can affect the buckling load, stability region and intensity of the instability (difference of upper and lower limit loads). [ABSTRACT FROM AUTHOR]