Theoretical models indicate that water ice is stable in the shallow subsurface (depths of <1-2 m) of Mars at high latitudes. These models have been mainly supported by the observed presence of large concentrations of hydrogen detected by the Gamma Ray Spectrometer suite of instruments on the Mars Odyssey spacecraft. The models and measurements are consistent with a water-ice table that steadily increases in depth with decreasing latitude. More detailed modelling has predicted that the depth at which water ice is stable can be highly variable, owing to local surface heterogeneities such as rocks and slopes, and the thermal inertia of the ground cover. Measurements have, however, been limited to the footprint (several hundred kilometres) of the Gamma Ray Spectrometer suite, preventing the observations from documenting more detailed water-ice distributions. Here I show that by observing the seasonal temperature response of the martian surface with the Thermal Emission Imaging System on the Mars Odyssey spacecraft, it is possible to observe such heterogeneities at subkilometre scale. These observations show significant regional and local water-ice depth variability, and, in some cases, support distributions in the subsurface predicted by atmospheric exchange and vapour diffusion models. The presence of water ice where it follows the depth of stability under current climatic conditions implies an active martian water cycle that responds to orbit-driven climate cycles. Several regions also have apparent deviations from the theoretical stability level, indicating that additional factors influence the ice-table depth. The high-resolution measurements show that the depth to the water-ice table is highly variable within the potential Phoenix spacecraft landing ellipses, and is likely to be variable at scales that may be sampled by the spacecraft.