We are interested in climate models introduced by Sellers in 1969 which takes the form of some nonlinear parabolic equation with a degenerate diffusion coefficient. We investigate here some inverse problem issue that consists in recovering the so-called insolation function. We not only solve the uniqueness question but also provide some strong stability result, more precisely unconditional Lipschitz stability.