The intense plume activity at the South Pole of Enceladus together with the recent detection of libration hints at an internal water ocean underneath the outer ice shell. However, the interpretation of gravity, shape, and libration data leads to contradicting results regarding the depth of ocean/ice interface and the total volume of the ocean. Here we develop an interior structure model consisting of a rocky core, an internal ocean, and an ice shell, which satisfies simultaneously the gravity, shape, and libration data. We show that the data can be reconciled by considering isostatic compensation including the effect of a few hundred meter thick elastic lithosphere. Our model predicts that the core radius is 180–185 km, the ocean density is at least 1030 kg/m3, and the ice shell is 18–22 km thick on average. The ice thicknesses are reduced at poles decreasing to less than 5 km in the south polar region.