Solar wind protons can charge exchange with the extensive hydrogen corona of Mars, resulting in a significant flux of energetic neutral atoms (ENAs). As these solar wind hydrogen ENAs precipitate into the upper atmosphere, they can experience electron attachment or detachment, resulting in populations of H− and H+, respectively, with upstream velocity. We seek to characterize the behavior of H− in the ionosphere of Mars through a combination of in situ data analysis and mathematical models. Observations indicate that measurable H− precipitation in the ionosphere of Mars is rare, occurring during only 1.8% of available observations. These events occur primarily during high energy solar wind conditions near perihelion. We also compare H− fluxes to those of H+ and find that H− fluxes are ∼4.5 times less than H+, indicating preferential conversion of hydrogen ENAs to H+. We develop a simple model describing the evolution of the charged and neutral fraction of ENAs and H− ions versus altitude. We find that 0.29 - 0.78% of ENAs are converted to H− for solar wind energies 1 - 3 keV. We also predict that the effects of photodetachment on the H-H− system are non-negligible.

We compare observations of hydrogen (H) and protons associated with Martian proton aurora activity, co-evaluating remote sensing and in situ measurements during these events. Following the currently understood relationship between penetrating protons and H energetic neutral atoms (ENAs) in the formation of proton aurora, we observe an expected correlation between the H Lyman-alpha (Ly-α) emission enhancement (used herein as a proxy for H-ENAs) and penetrating proton flux. However, we observe a notable spread in the trend between these two datasets. We find that this spread is contemporaneous with one of two major impacting events: high dust activity or extreme solar activity. Proton aurora events exhibiting a relative excess in penetrating proton flux compared to Ly-α enhancement tend to correspond with periods of high dust activity. Conversely, proton aurora events exhibiting a relative deficit of penetrating proton flux compared to Ly-α enhancement are qualitatively associated with periods of extreme solar activity. Moreover, we find that the largest proton aurora events occur during concurrent dust storm and solar events, primarily due to the compounding intensified increase in H column density above the bow shock. Finally, we present a simplified empirical estimate for Ly-α emission enhancement during proton aurora events based on observed penetrating proton flux and a knowledge of local dust/solar activity at the time, providing a straightforward method for predicting auroral activity when direct observations are not available. The results of this study advance our understanding of the interconnected relationship between H and protons during Martian proton aurora activity.