Assuming our Solar System as typical, exomoons may outnumber exoplanets. If their habitability fraction is similar, they would thus constitute the largest portion of habitable real estate in the Universe. Icy moons in our Solar System, such as Europa and Enceladus, have already been shown to possess liquid water, a prerequisite for life on Earth. We intend to investigate under what circumstances small, icy moons may sustain subsurface oceans and thus be "subsurface habitable". We pay specific attention to tidal heating. We made use of a phenomenological approach to tidal heating. We computed the orbit averaged flux from both stellar and planetary (both thermal and reflected stellar) illumination. We then calculated subsurface temperatures depending on illumination and thermal conduction to the surface through the ice shell and an insulating layer of regolith. We adopted a conduction only model, ignoring volcanism and ice shell convection as an outlet for internal heat. In doing so, we determined at which depth, if any, ice melts and a subsurface ocean forms. We find an analytical expression between the moon's physical and orbital characteristics and the melting depth. Since this expression directly relates icy moon observables to the melting depth, it allows us to swiftly put an upper limit on the melting depth for any given moon. We reproduce the existence of Enceladus' subsurface ocean; we also find that the two largest moons of Uranus (Titania & Oberon) could well sustain them. Our model predicts that Rhea does not have liquid water. Habitable exomoon environments may be found across an exoplanetary system, largely irrespective of the distance to the host star. Small, icy subsurface habitable moons may exist anywhere beyond the snow line. This may, in future observations, expand the search area for extraterrestrial habitable environments beyond the circumstellar habitable zone.

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist.