Sub-field of view surface thermal modeling of Cassini CIRS observations of Rhea during south polar winter

Rhea's exosphere is thought to originate from sources of carbon, water ice and other volatiles that arrived at Rhea by bombardment. Its seasonal variability is directly driven by polar surface temperatures allowing surface adsorption, and the persistence of source volatiles require that seasonal temperatures remain sufficiently cold to retain them. Cassini CIRS detected temperatures in Rhea's winter polar region of 23 K, amongst the coldest measured in the Solar System, but with a relatively large footprint we seek to add value to these observations by modeling sub field-of-view (FOV) temperature distribution and examining how the coldest scenes evolve on a seasonal basis. A simple, rough surface 1-dimensional thermal model is developed using a digital elevation map as input to a thermal model. We compared averaged rough and flat modeled FOVs to CIRS temperatures for a set of case study observations in the south polar region in winter darkness and found they both agree within expected CIRS error in all cases. We develop an asymmetric estimate of CIRS FOV temperature uncertainty, which is particularly important for very cold spectra to accurately represent the sensitivity of the instrument. This approach provides a conservative upper temperature limit to spectral fits whilst highlighting there is often little information to constrain the lower bound of temperature uncertainty in these cases. The distribution of modeled sub-FOV facet temperatures is explored for the full range of azimuth angles and slope gradients. More extensive and cooler temperatures were found in the rough model, complemented by fewer but warmer areas than the flat model. We find temperature contrasts of individual model facets of up to 15 K warmer and 11 K colder within some CIRS FOV when the rough model was compared flat, in a case study of scenes in winter darkness. This is due to the persistence of the seasonal thermal wave beneath the surface. We tested model sensitivity to thermal input parameters (thermal inertia and bolometric Bond albedo). These values are challenging to constrain due to limited observations and measurement noise and are expected to vary with subtle changes in surface characteristics. We found that within 20° of the pole, temperatures do not rise significantly above 80 K all year, implying that a variety of simple organic molecules, linear amides and other carbon-containing compounds would remain stable on a quasi-permanent basis, potentially until photolytic or