This thesis considers the heavy ion dynamics due to ion-cyclotron resonance energization processes that take place in the turbulent region of the Earth’s topside, high latitude ionosphere. We simulate the impact of this transverse heating process upon energies and velocity distribution functions of outflowing oxygen ions (O+) in the approximate altitude range of 800 km to 15,000 km. To do so most effectively, we use a single particle tracing model that precisely reproduces the small-scale wave-particle interaction of broadband extremely low frequency (BBELF) waves with the ions’ cyclotron motions, leading to the upward acceleration of ions in type-II ion outflows and ion conics. Instead of employing the guiding center approximation, as is u...