Volatile distribution inversion and rotation analysis of comet 103P/Hartley 2 using nongravitational effects

The rotational states of comet 103P/Hartley 2 (103P) during its perihelion passage are strongly modulated by nongravitational forces. By fitting its attitude-variation measurements, this study aims to infer the spatial distributions of H_2O and CO_2 ices on the nucleus and to assess how these inhomogeneities drive the nucleus's rotational behavior. We adopted a dust-mantle thermal model of the nucleus that incorporates H_2O and CO_2 ices as the principal volatiles. A coupled attitude-thermal simulation framework was developed, in which a neural-network surrogate greatly accelerates the iterative computations. The nucleus was partitioned into three morphological regions---large end, waist, and small end---and a localized activity model described the uneven sublimation across these zones. Our model successfully reproduces both the observed rotational variations and the measured volatile production rates of 103P during its 2010 encounter. At the observation instant, in addition to the pronounced CO_2 sublimation at the sunlit small end, the dust mantle covering the large end insulated the underlying CO_2, allowing it to continue sublimating on the shadowed side. We find that the H_2O ice-front depth is on the order of millimeters, whereas the CO_2 ice front lies between centimeters and decimeters, with the large end exhibiting stronger activity. The fitted CO_2 production rate closely matches the observations, while direct H_2O sublimation contributes only a minor fraction. The total H_2O production rate of 103P may primarily originate from H_2O ice ejected by CO_2 sublimation at the nucleus's ends. Moreover, CO_2 sublimation at the large end produces the majority of the recoil torque, driving the nucleus's enhanced long-axis excitation following the encounter.