Modeling Io's sublimation-driven atmosphere: Gas dynamics and radiation emission

Io's sublimation-driven atmosphere is modeled using the direct simulation Monte Carlo method. These rarefied gas dynamics simulations improve upon earlier models by using a three-dimensional domain encompassing the entire planet computed in parallel. The effects of plasma impact heating, planetary rotation, and inhomogeneous surface frost are investigated. Circumplanetary flow is predicted to develop from the warm subsolar region toward the colder night-side. The non-equilibrium thermal structure of the atmosphere, including vibrational and rotational temperatures, is also presented. Io's rotation leads to an asymmetric surface temperature distribution which is found to strengthen circumplanetary flow near the dusk terminator. Plasma heating is found to significantly inflate the atmosphere on both day- and night-sides. The plasma energy flux also causes high temperatures at high altitudes but permits relatively cooler temperatures at low altitudes near the dense subsolar point due to plasma energy depletion. To validate the atmospheric model, a radiative transfer model was developed utilizing the backward Monte Carlo method. The model allows the calculation of the atmospheric radiation from emitting/absorbing and scattering gas using an arbitrary scattering law and an arbitrary surface reflectivity. The model calculates the spectra in the v ₂ vibrational band of SO ₂ which are then compared to the observational data.