^1,2Ery C. Hughes,²Katherine de Kleer,²John Eiler,³Francis Nimmo,⁴Kathleen Mandt,⁵Amy E. Hofmann
Journal of Geophysical Research (Planets)(in Press) Open Access Link to Article [https://doi.org/10.1029/2023JE008086]
¹Te Pū Ao, GNS Science, National Isotope Centre and Avalon, Lower Hutt, Aotearoa New Zealand
²Division of Geological and Planetary Science, Caltech, Pasadena, CA, USA
³Earth & Planetary Sciences Department, University of California Santa Cruz, Santa Cruz, CA, USA
⁴NASA Goddard Space Flight Center, Greenbelt, MD, USA
⁵Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA
Published by arrangement with John Wiley & Sons

Stable isotope fractionation of sulfur offers a window into Io’s tidal heating history, which is difficult to constrain because Io’s dynamic atmosphere and high resurfacing rates leave it with a young surface. We constructed a numerical model to describe the fluxes in Io’s sulfur cycle using literature constraints on rates and isotopic fractionations of relevant processes. Combining our numerical model with measurements of the 34S/32S ratio in Io’s atmosphere, we constrain the rates for the processes that move sulfur between reservoirs and model the evolution of sulfur isotopes over time. Gravitational stratification of SO2 in the upper atmosphere, leading to a decrease in 34S/32S with increasing altitude, is the main cause of sulfur isotopic fractionation associated with loss to space. Efficient recycling of the atmospheric escape residue into the interior is required to explain the 34S/32S enrichment magnitude measured in the modern atmosphere. We hypothesize this recycling occurs by SO2 surface frost burial and SO2 reaction with crustal rocks, which founder into the mantle and/or mix with mantle-derived magmas as they ascend. Therefore, we predict that magmatic SO2 plumes vented from the mantle to the atmosphere will have lower 34S/32S than the ambient atmosphere, yet are still significantly enriched compared to solar-system average sulfur. Observations of atmospheric variations in 34S/32S with time and/or location could reveal the average mantle melting rate and hence whether the current tidal heating rate is anomalous compared to Io’s long-term average. Our modeling suggests that tides have heated Io for >1.6 Gyr if Io today is representative of past Io.