Changing Chemistry, Changing Dynamics: Modeling the Response of Jupiter's magnetosphere to Io's volcanic eruptions

I am modeling the Io plasma torus around Jupiter. The planet's powerful gravity flexes the moon, inducing high levels of volcanism. The material ejected from these volcanoes forms a torus of gas along the moon's orbit, and the gas particles are ionized by energetic electrons, forming plasma that responds to Jupiter's titanic magnetosphere. My model simulates the motion of this plasma with respect to the planet's magnetic field, and seeks to explore several mysteries about the Jovian system. Where do the energetic electrons come from? Why does the plasma persistently lag behind the magnetic field lines rotating with the planet? How does changing the size and length of eruptions on Io affect plasma conditions elsewhere in the torus, such as at Europa, the next moon out? In addition to answering questions specific to Jupiter, my model looks to provide more information on the general response of plasma inside a magnetosphere to changing conditions, such as the rate of solar wind entering Earth's atmosphere helping determine the magnitude of aurora in Alaskan skies.