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.

We present the first observations of the dayside coronal oxygen emission in far ultraviolet (FUV) measured by the Emirates Mars Ultraviolet Spectrometer (EMUS) onboard the Emirates Mars Mission (EMM). The high sensitivity of EMUS is providing an opportunity to observe the tenuous oxygen corona in FUV, which is otherwise difficult to observe. Oxygen resonance fluorescence emission at 130.4 nm provides a measurement of the upper atmospheric and exospheric oxygen. More than 500 oxygen corona profiles are constructed using the long-exposure time cross-exospheric mode (OS4) of EMUS observations. These profiles range from ~200 km altitude up to several Mars radii (>6 RM) across all seasons and for two Mars years. Our analysis shows that OI 130.4 nm is highly correlated with solar irradiance (solar photoionizing and 130.4 nm illuminating irradiances) as well as changes in the Sun-Mars distance. The prominent short term periodicity in oxygen corona brightness is consistent with the solar rotation period (quasi-27-days). A comparison between the perihelion seasons of Mars Year (MY) 36 and MY 37 shows interannual variability with enhanced emission intensities during MY 37, due to the rise of Solar Cycle 25. These observations show a highly variable oxygen corona, which has significant implications on constraining the photochemical escape of atomic oxygen from Mars.

We present a comprehensive study of the nightside discrete electron aurora phenomenon on Mars, utilizing observations from EMUS onboard EMM. The oxygen emission at 130.4 nm is by far the brightest FUV auroral emission line observed at Mars. We identify auroral pixels in OI 130.4 nm disk observations, with higher sensitivity than previously possible. Our statistical analysis reveals regional, SZA, local time, and seasonal dependencies of auroral occurrence. Higher occurrence of aurora is observed in regions of open magnetic topology and vertical crustal magnetic fields. Aurora occurs more frequently closer to the terminator and is more likely on the dusk versus dawn sides of the night hemisphere. A pronounced auroral feature appears close to midnight local times in the southern hemisphere, consistent with the “spot” of energetic electron fluxes previously identified in the MGS data. The auroral spot is more frequent after midnight than before. Additionally, some regions on Mars are “aurora voids” where essentially no aurora occurs. The non-crustal field aurora exhibits a seasonal dependence, with major enhancements around Ls 235° (near perihelion) and Ls 30°. This is in line with the seasonal variability in ionospheric TEC observed by Mars Express, which is in turn related to the variability of solar irradiance and thermospheric density. Aurora occurrence also shows an increase with the rise of Solar Cycle 25. These observations not only shed light on where and when Martian aurora occurs, but also add to our understanding of Mars’ magnetic environment and its interaction with the heliospheric environment.

The D/H ratio in water on Mars, Rwater, is 4–6× the Earth ratio, signifying past water loss to space. Recently, measurements have revealed high values of the D/H ratio in hydrogen, Ratomic, in the thermosphere during southern summer. Here, we use a photochemical model to explore the potential drivers of Ratomic, testing three: thermospheric temperatures, excess mesospheric water, and changing insolation. We find that Ratomic can achieve values between 15× the Earth ratio (due to water) and 25× the Earth ratio (due to temperature). The effects arise because H escape is diffusion-limited, while D escape is energy-limited. Our results underscore how Ratomic reflects mesospheric dynamics, and the need for concurrent measurements of mesospheric water, thermospheric temperatures, and Ratomic to understand seasonal changes in the martian water cycle and atmospheric loss.

We evaluate for the first time the large-scale influence of the upstream interplanetary magnetic field (IMF) magnitude and orientation on Martian proton aurora. We identify a correlation between proton aurora activity and IMF magnitude that varies seasonally, with highest correlations around Mars southern summer solstice when proton aurora activity reaches an annual high. We identify preferentially increased proton aurora emission enhancements at IMF cone and clock angles consistent with the shape of the Parker spiral at Mars’s heliocentric distance. An increased proton aurora occurrence rate is observed for near-radial IMF orientations, particularly evident when the hydrogen corona and bow shock are annually decreased. Lastly, our results suggest that dayside magnetic reconnection may influence proton aurora, which are observed to have preferentially higher occurrence under strong ±Bz orientations. These findings enhance our understanding of the upstream IMF’s influence on Martian proton aurora as an important driver for auroral brightness and occurrence.

Proton aurora are the most commonly observed yet least studied type of aurora at Mars. In order to better understand the physics and driving processes of Martian proton aurora, we undertake a multi-model comparison campaign. We compare results from four different proton/hydrogen precipitation models with unique abilities to represent Martian proton aurora: Jolitz model (3-D Monte Carlo), Kallio model (3-D Monte Carlo), Bisikalo/Shematovich et al. model (1-D kinetic Monte Carlo), and Gronoff et al. model (1-D kinetic). This campaign is divided into two steps: an inter-model comparison and a data-model comparison. The inter-model comparison entails modeling five different representative cases using similar constraints in order to better understand the capabilities and limitations of each of the models. Through this step we find that the two primary variables affecting proton aurora are the incident solar wind particle flux and velocity. In the data-model comparison, we assess the robustness of each model based on its ability to reproduce a MAVEN/IUVS proton aurora observation. All models are able to effectively simulate the data. Variations in modeled intensity and peak altitude can be attributed to differences in model capabilities/solving techniques and input assumptions (e.g., cross sections, 3-D versus 1-D solvers, and implementation of the relevant physics and processes). The good match between the observations and multiple models gives a measure of confidence that the appropriate physical processes and their associated parameters have been correctly identified, and provides insight into the key physics that should be incorporated in future models.

Although deuterium (D) on Mars has received substantial attention, the deuterated ionosphere remains relatively unstudied. This means that we also know very little about non-thermal D escape from Mars, since it is primarily driven by excess energy imparted to atoms produced in ion-neutral reactions. Most D escape from Mars is expected to be non-thermal, highlighting a gap in our understanding of water loss from Mars. In this work, we set out to fill this knowledge gap. To accomplish our goals, we use an upgraded 1D photochemical model that fully couples ions and neutrals and does not assume photochemical equilibrium. To our knowledge, such a model has not been applied to Mars previously. We model the atmosphere during solar minimum, mean, and maximum, and find that the deuterated ionosphere behaves similarly to the H-bearing ionosphere, but that non-thermal escape on the order of 8000-9000 cm-2s-1 dominates atomic D loss under all solar conditions. The total fractionation factor, f, is 0.04–0.07, and integrated water loss is 147–158 m GEL. This is still less than geomorphological estimates. Deuterated ions at Mars are likely difficult to measure with current techniques due to low densities and mass degeneracies with more abundant H ions. Future missions wishing to measure the deuterated ionosphere in situ will need to develop innovative techniques to do so.

Mars’ ultraviolet airglow has been used to study its upper atmosphere for over four decades. Identifying variations in emission features has provided information on composition, density and temperature. The Emirates Mars Ultraviolet Spectrometer onboard the Emirates Mars Mission observes Mars’ airglow at Far and Extreme UV wavelengths. Variations in disk emission features are studied, with a focus on O I 1304 Å, CO Fourth Positive Group and C I 1561 Å. All show variations with local time and emission angle as expected. Dawn-dusk asymmetry observed is attributed to local time differences in advection. Variations in the brightness of several dayglow features, including 1304 Å, with irregular shapes are noted in around 25% of the disk images. These display some local time and hemispheric asymmetry in their occurrence rates. Examination of their spatial structure, occurrence, and spectra suggests these are associated with variations in composition and photoelectron flux.

Much of the water that once flowed on the surface of Mars was lost to space long ago, and the total amount lost remains unknown. Clues to the amount lost can be found by studying hydrogen (H) and its isotope deuterium (D), both of which are produced when atmospheric water molecules H$_2$O and HDO dissociate. The freed H and D atoms then escape to space at different rates due to their different masses, leaving an enhanced D/H ratio. The rate of change of D/H is referred to as the fractionation factor $f$. Both the D/H ratio and $f$ are necessary to estimate water loss; thus, if we can constrain the range of $f$, we will be able to estimate water loss more accurately. In this study, we use a 1D photochemical model of the Martian atmosphere to determine how $f$ depends on assumed temperature and water vapor profiles. We find that for most Martian atmospheric conditions, $f$ varies between $10^{-1}$ and $10^{-5}$; for the standard Martian atmosphere, $f=0.002$ for thermal escape processes, and $f\approxeq0.06$ when both thermal and non-thermal escape are considered. Using these results, we estimate that Mars has lost at minimum 66-123 m GEL of water. Our results demonstrate that the value of $f$ is almost completely controlled by the amount of non-thermal escape of D, and that photochemical modeling studies that include fractionation must thus model both neutral and ion processes throughout the atmosphere.

Lower atmosphere variations in the martian water vapour and hydrogen abundance during the Mars Year (MY) 34 C storm from LS=326.1◦-333.5◦ and their associated effect on hydrogen escape are investigated using a multi-spacecraft assimilation of atmospheric retrievals into a Martian global circulation model. Elevation of the hygropause and associated increase in middle atmosphere hydrogen at the peak of the MY 34 C storm led to a hydrogen escape rate of around 1.4×109 cm−2s−1 , meaning the MY 34 C storm enhanced water loss rates on Mars to levels observed during global-scale dust storms. The water loss rate during the MY 34 C storm (a loss of 15% of the total annual water loss during only 5% of the year) was three times stronger than the weak MY 30 C storm assimilation, demonstrating that interannual variations in C storm strength must be considered when calculating the integrated loss of water on Mars.

The Mars Atmosphere and Volatile Evolution (MAVEN) mission instrument suite includes an ultraviolet echelle spectrograph with high-spectral resolution designed to resolve D and H Lyman-α emissions. The high-spectral resolution mode was previously characterized in the lab and in the cruise phase to Mars and had been calibrated using observations and models of interplanetary hydrogen Lyman-α emissions. This work presents improved characterizations of the high-spectral resolution mode using in-orbit observations that allow for more robust detections of the faint D Lyman-α emission line. Additionally, the instrument was re-calibrated using simultaneous and comparable observations made with the Hubble Space Telescope high-spectral resolution instrument. Comparisons to Lyman-α observations made with the low-resolution UV channel on the spectrometer, that had been calibrated with stars, showed consistency in the brightness values for measurements obtained at similar observational conditions. The combined upgrades to the faint-emission fitting and new calibration techniques of the MAVEN echelle channel have resulted in an improved data-reduction pipeline with favorable implications for the science utility of D and H Lyman-α emissions.

Water escape on Mars has recently undergone a paradigm shift with the discovery of unexpected seasonal variations in the population of hydrogen atoms in the exosphere where thermal escape occurs and results in water lost to space. This discovery led to the hypothesis that, contradicting the accepted pathway, atomic hydrogen in the exosphere was not only produced by molecular hydrogen but mostly by high altitude water vapor. Enhanced presence of water at high altitude during southern spring and summer, due to atmospheric warming and intensified transport, favors production of H through photolysis ionized chemistry of water molecules and thus appears to be the main cause of the observed seasonal variability in escaping hydrogen. This hypothesis is supported by the observation of large concentrations of water vapor between 50 km and 150 km during the southern summer solstice and global dust events. Using a simplified yet representative air parcel transport model, we show that in addition to the formation of atomic hydrogen from water photolysis above 80 km, a major fraction of the exospheric hydrogen is formed at altitudes as low as 60 km and is then directly advected to the upper atmosphere. Comparing the injection modes of a variety of events (global dust storm, perihelion periods, regional storm), we conclude that southern spring/summer controls H production and further ascent into the upper atmosphere on the long term with direct implication for water escape.