As a long-term energy source, tidal heating in subsurface oceans of icy satellites can influence their thermal, rotational, and orbital evolution, and the sustainability of oceans. We present a new theoretical treatment for tidal heating in thin subsurface oceans with overlying incompressible elastic shells of arbitrary thickness. The stabilizing effect of an overlying shell damps ocean tides, reducing tidal heating. This effect is more pronounced on Enceladus than on Europa because the effective rigidity on a small body like Enceladus is larger. For the range of likely shell and ocean thicknesses of Enceladus and Europa, the thin shell approximation of Beuthe (2016) is generally accurate to less than about 4%.The time-averaged surface distribution of ocean tidal heating is distinct from that due to dissipation in the solid shell, with higher dissipation near the equator and poles for eccentricity and obliquity forcing respectively. This can lead to unique horizontal shell thickness variations if the shell is conductive. The surface displacement driven by eccentricity and obliquity forcing can have a phase lag relative to the forcing tidal potential due to the delayed ocean response. For Europa and Enceladus, eccentricity forcing generally produces greater tidal amplitudes due to the large eccentricity values relative to the obliquity values. Despite the small obliquity values, obliquity forcing generally produces larger phase lags due to the generation of Rossby-Haurwitz waves. If Europa's shell and ocean are respectively 10 and 100 km thick, the tide amplitude and phase lag are 26.5 m and $<1$ degree for eccentricity forcing, and $<2.5$ m and $<18$ degrees for obliquity forcing. Measurement of the obliquity phase lag (e.g. by Europa Clipper) would provide a probe of ocean thickness

Laser altimetry is a powerful tool for addressing the major objectives of planetary physics and geodesy, and have been applied in planetary explorations of the Moon, Mars, Mercury, and the asteroids Eros, and Itokawa. The JUpiter Icy Moons Explorer (JUICE), led by European Space Agency (ESA), has started development to explore the emergence of habitable worlds around gas giants. The Ganymede Laser Altimeter (GALA) will be the first laser altimeter for icy bodies, and will measure the shape and topography of the large icy moons of Jupiter, (globally for Ganymede, and using flyby ground-tracks for Europa and Callisto). Such information is crucial for understanding the formation of surface features and can tremendously improve our understanding of the icy tectonics. In addition, the GALA will infer the presence or absence of a subsurface ocean by measuring the tidal and rotational responses. Furthermore, it also improves the accuracy of gravity field measurements reflecting the interior structure, collaborating with the radio science experiment. In addition to range measurements, the signal strength and the waveform of the laser pulses reflected from the moon's surface contain information about surface reflectance at the laser wavelength and small scale roughness. Therefore we can infer the degrees of chemical and physical alterations, e.g., erosion, space weathering, compaction and deposition of exogenous materials, through GALA measurements without being affected by illumination conditions. JUICE spacecraft carries ten science payloads including GALA. They work closely together in a synergistic way with GALA being one of the key instruments for understanding the evolution of the icy satellites Ganymede, Europa, and Callisto.