Thermal modeling of cometary nuclei

A model of the sublimation of volatile ices from a cometary nucleus is presented which includes the effects of (1) diurnal heating and cooling, (2) rotation period and pole orientation, (3) the thermal properties of the ice and subsurface layers, and (4) the contributions from coma opacity, scattering and thermal emission where the properties of the coma are derived from the integrated rate of volatile production by the nucleus. In applying the model to the case of the 1986 apparition of Halley's comet, it is found that the generation of a cometary dust coma increases the total energy reaching the Halley nucleus due to the greater geometrical cross-section of the coma as compared with the bare nucleus. The calculated coma opacity of Halley is about 0.2 at 1 AU from the sun and 1.2 at perihelion. Possible consequences of the results obtained for the generation of nongravitational forces, volatile production rates for comets and cometary lifetimes against sublimation are discussed.