Variability and geologic associations of volcanic activity on Io in 2001-2016

We investigate the spatial distribution of volcanic activity on Io using nearest neighbor, mean pairwise spacing, and mean latitude statistics with various classification schemes. The analysis confirms previous findings in that the heat dissipation appears to be primarily concentrated in the asthenosphere resulting in a high time-averaged surface heat flux at low latitudes. Our observations show significant spatial deviations do exist from the asthenosphere heat dissipation model while also suggesting a deeper source of magma ascent to be present as well, supporting conclusions from previous analyses of primarily spacecraft data (Veeder et al., 2012; Hamilton, 2013; Davies et al., 2015). From a temporal perspective, there are signs of significant variations in the distribution of global heat flux, as volcanoes undetected, and probably dormant, during the Galileo encounters subsequently erupted and remained active during our observations. We also use the on 3.8-µm radiant intensity timelines of individual hot spots, along with the distribution of extensive lava fields in relation to detected activity, as a means to investigate possible connections between hot spots and short timescale, spatio-temporal variations in the global heat flux distribution. We conclude that while the global heat flux distribution remains relatively constant over decadal timescales, there is evidence that significant deviations do occur potentially as a result of mountain forming processes or triggering mechanisms between eruptions.