Ion Escape Dynamics in the Martian Upper Ionosphere: A Self-Similar Plasma Expansion Approach

Abstract This study investigates the dynamics of positive ion escape in the Martian upper ionosphere at altitudes of 850–1200 km using a self-similar hydrodynamic plasma expansion model. The model incorporates Kappa-distributed suprathermal electrons and two oxygen (O) positive ion species, O+ and O₂^+, with plasma parameters based on MAVEN mission observations. Numerical solutions of the coupled nonlinear fluid equations reveal that the O₂^+-e temperature ratio and the suprathermal electron κ −index strongly govern the escape profiles—density, velocity, and electrostatic potential, whereas variations in O+ temperature and density ratios exert only minor influence. The estimated total escape rates are approximately 7. 4 × 1023 ions s−1 for O₂^+ and 5. 7 × 1022 ions s−1 for O+, corresponding to mass loss rates of 0. 039 kg s−1 and 0. 0015 kg s−1, respectively. These values provide quantitative constraints on Martian atmospheric escape processes and align well with the observational estimates. Therefore, our model offers valuable insights into the mechanisms governing ionospheric plasma expansion, thereby contributing to a deeper understanding of the long-term evolution of planetary atmospheres and the potential habitability of Mars and similar terrestrial planets.