On the Light-Absorbing Surface Layer of Cometary Nuclei: II. Thermal Modeling

The classical way to treat absorption of solar light in thermophysical modeling of cometary nuclei (and other ice-rich bodies such as jovian satellites) has been to assume complete opaqueness of the surface material. However, as shown by Davidsson and Skorov (2002, Icarus156, 223-248), substantial light penetration can occur in porous ice even if it is very dusty, implying that gradual absorption of energy in a surface layer should be accounted for.We present a thorough comparison between a surface energy absorption model and a layer energy absorption model, for various combinations of heliocentric distances, conductivities, opacities, pore sizes, and rotational periods relevant for cometary nuclei, by fully solving the coupled differential equations of heat transfer and gas diffusion. We find substantial differences between the models in terms of gas production rate, thermal lag angle, surface temperature, and the origin of coma molecules. For example, the surface energy absorption model overestimates the total gas production by a factor of 2-7, underestimates the lag angle by a factor of 2-3, and places the origin of coma molecules at the surface, instead of the near-surface interior.