We present a numerical study of the shape taken by a spherical elastic surface when the volume it encloses is decreased. For the range of 2D parameters where such a surface may model a thin shell of an isotropic elastic material, the mode of deformation that develops a single depression is investigated in detail. It occurs via buckling from sphere toward an axisymmetric dimple, followed by a second buckling where the depression loses its axisymmetry through folding along portions of meridians. For the thinnest shells, a direct transition from the spherical conformation to the folded one can be observed. We could exhibit unifying master curves for the relative volume variation at which first and second buckling occur, and clarify the role of the Poisson's ratio. In the folded conformation, the number of folds and inner pressure are investigated, allowing us to infer shell features from mere observation and/or knowledge of external constraints.