# Detection of visible-wavelength aurora on Mars

Elise W Knutsen, Timothy H Mcconnochie, Mark Lemmon, Chris Donaldson, Raymond Francis, Carey Legett, Shayla B Viet, Lauriane Soret, Daniel Toledo, Victor Apéstigue, et al.

► To cite this version:

Elise W Knutsen, Timothy H Mcconnochie, Mark Lemmon, Chris Donaldson, Raymond Francis, et al. Detection of visible-wavelength aurora on Mars. Science Advances, 2025, 11 (20), pp.eads1563. 10.1126/sciadv.ads1563. hal-05073417

HAL Id: hal-05073417

https://hal.science/hal-05073417v1

Submitted on 19 May 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

PLANETARY SCIENCE

# Detection of visible-wavelength aurora on Mars

Elise W. Knutsen<fnref n="1" />*, Timothy H. McConnochie², Mark Lemmon², Chris Donaldson³,
Raymond Francis⁴, Carey Legett⁵, Shayla B. Viet⁶, Lauriane Soret⁷, Daniel Toledo⁸,
Victor Apéstigue⁸, Olivier Witasse⁹, Franck Montmessin¹⁰, Rebecca Jolitz¹¹,
Nicolas M. Schneider¹¹, Leslie Tamppari⁴, Agnès Cousin¹², Roger C. Wiens¹³, Sylvestre Maurice¹²,
James F. Bell III¹⁴, Olivier Forni¹², Jeremie Lasue¹², Paolo Pilleri¹², Tanguy Bertrand¹⁵,
Priya Patel⁴,¹⁶, Susanne Schröder¹⁷, Shannon Curry¹¹, Christina O. Lee¹⁸, Ali Rahmati¹⁸

Mars hosts various auroral processes despite the planet's tenuous atmosphere and lack of a global magnetic field. To date, all aurora observations have been at ultraviolet wavelengths from orbit. We describe the discovery of green visible-wavelength aurora, originating from the atomic oxygen line at 557.7 nanometers, detected with the SuperCam and Mastcam-Z instruments on the Mars 2020 Perseverance rover. Near-real-time simulations of a Mars-directed coronal mass ejection (CME) provided sufficient lead-time to schedule an observation with the rover. The emission was observed 3 days after the CME eruption, suggesting that the aurora was induced by particles accelerated by the moving shock front. To our knowledge, detection of aurora from a planetary surface other than Earth has never been reported, nor has visible aurora been observed at Mars. This detection demonstrates that auroral forecasting at Mars is possible, and that during events with higher particle precipitation, or under less dusty atmospheric conditions, aurorae will be visible to future astronauts.

Copyright © 2025 The
Authors, some rights
reserved; exclusive
licensee American
Association for the
Advancement of
Science. No claim to
original U.S.
Government Works.
Distributed under a
Creative Commons
Attribution License 4.0
(CC BY).

INTRODUCTION

Planetary aurorae have been observed on every planet with an at-
mosphere in our Solar System. In response to the precipitation of
energetic particles, emissions span x-ray to radio wavelengths and
can be used to constrain the incoming particle energy distribution
(1-3) and as a means of remotely sensing the conditions of the emis-
sion source region (4, 5).

Mars is, as far as we know, the only planet in our Solar System with a hybrid magnetosphere, with elements of both the induced magnetosphere originating from the interplanetary magnetic field (IMF) lines being draped around Mars' conductive ionosphere and regions of crustal magnetization (6, 7) mostly confined to the ancient southern highlands (8). This composite field structure creates a complex magnetic morphology leading to diverse auroral phenomena, where questions of solar particle injection, precipitation, and transport remain open, particularly with respect to magnetotail instabilities and acceleration mechanisms and the role of magnetic topology and magnetic reconnection (9-11).

Several morphologies of aurorae have so far been observed on
Mars: localized discrete and patchy aurora (12, 13), global diffuse
aurora (14), dayside proton aurora (15), and large-scale sinuous aurora

(10, 16). These can be categorized by the precipitating particle ener-
gies, namely, solar wind aurora (proton aurora), suprathermal aurora
(discrete, sinuous aurora and patchy aurora, with electron energies of
tens to hundreds of electron volts (10, 17, 18), and charged solar ener-
getic particle (SEP) aurora (diffuse aurora). SEP and suprathermal
aurora are observed intermittently on the martian nightside. Con-
fined discrete aurorae appear frequently but are short-lived [tens of
minutes to hours (16, 19) and are mostly observed near the crustal
fields (18, 19)], while the patchy aurorae appear in nonmagnetized
regions as well and often in the north where there are no crustal
fields. The recently discovered sinuous aurora appears in 3% of ob-
servations from the Emirates Mars Mission Emirates Mars Ultraviolet
Spectrograph and is believed to indicate a highly dynamic magneto-
tail (10, 16). Diffuse aurorae are strongly correlated with SEP events
(20, 21) and can be forecasted whether the solar event is identified
shortly after eruption.

SEPs are accelerated by two main processes: magnetic reconnec-
tion processes during solar flare activity (impulsive SEP events) and
wave acceleration by coronal mass ejection (CME)-driven shocks
(gradual SEP events) (22, 23). Impulsive events accelerate the parti-
cles near the Sun, while gradual events accelerate particles as the
CME-driven shock front propagates radially through the Solar Sys-
tem (24). SEPs travel along the Parker spiral (25) and can thus be
detected at Mars if the planet is magnetically connected to a local
acceleration site (at the Sun for impulsive events or to the moving
shock source for gradual events). The travel time of impulsive SEPs
to Mars is on the order of a few tens of minutes (for relativistic elec-
trons) to hours (e.g., for ~1-MeV protons) (26, 27). Meanwhile, the
travel time of gradual SEPs depends on when Mars is magnetically
connected to the shock source as the interplanetary CME (ICME)
propagates outward. If the ICME is Mars directed, then an impact
can occur within 3 days after the initial eruption for a fast ICME
(e.g., >1000 km/s) (28). The durations of impulsive and gradual SEP
events vary; impulsive events can last for a few hours, while gradual
events can disturb the planetary environment for several days
(22, 29).

<footnote n="1">¹University of Oslo, CENSSS, Oslo, Norway. ²Space Science Institute, Boulder, CO, USA. ³Malin Space Science Systems, San Diego, CA, USA. ⁴Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA, USA. ⁵Los Alamos National Laboratory, Los Alamos, NM, USA. ⁶Norwegian University of Science and Technology, Trondheim, Norway. ⁷Laboratoire de Physique Atmosphérique et Planétaire, STAR Institute, Université de Liège, Liege, Belgium. ⁸Instituto Nacional De Técnica Aeronspacial (INTA), Madrid, Spain. ⁹ESA/ESTEC, Noordwijk, Netherlands. ¹⁰LATMOS/IPSL, UVSQ Université Paris-Saclay, Sorbonne Université, CNRS, Paris, France. ¹¹Laboratory for Atmospheric and Space Physics, University of Colorado, Boulder, CO, USA. ¹²IRAP-Institut de recherche en astrophysique et planétologie, Toulouse, France. ¹³Purdue University, West Lafayette, IN, USA. ¹⁴School of Earth and Space Exploration, Arizona State University, Tempe, AZ, USA. ¹⁵LESIA, Paris Observatory, Paris, France. ¹⁶Mullard Space Science Laboratory, University College London, England. ¹⁷Deutsches Zentrum für Luft- und Raumfahrt (DLR), Institute of Optical Sensor Systems, Berlin, Germany. ¹⁸Space Sciences Laboratory, University of California, Berkeley, CA, USA.</footnote>

<footnote>*Corresponding author. Email: elisewkn@uio.no</footnote>
Mars' magnetosphere largely deflects and decelerates solar wind particles, but SEPs (and particularly higher-energy (a few mega-electron volts to ~100 MeV) protons (30) penetrate deep into the martian atmosphere [e.g., (29–31)]). The material pile-up in front of the leading shock of ICMEs can disturb the fields around Mars sufficiently to increase coupling between Mars' magnetosphere and atmosphere, allowing for the injection of a large flux of particles, subsequently leading to auroral emissions (32). Using a Monte Carlo Particle Transport In Planetary atmospheres model, Nakamura et al. (33) demonstrated that the auroral peak emission altitude varies depending on the precipitating particle species; proton-induced diffuse aurorae being both lower and brighter than electron-induced aurorae. Theoretical work has shown that 100-keV electrons are capable of inducing auroral emission at 60-km altitude (21, 34), yet models predict that only a small fraction of the SEP electrons (energies in the 10- to 200-keV range) penetrates into the atmosphere (30). It has been hypothesized that an interplanetary disturbance, caused, for example, by ICMEs, might be required to increase the access of external particles to the atmosphere (32).

SEP-driven aurorae have so far been observed to be long-lasting (timescale of days) and cover the entire nightside disk of Mars (14, 21). The correlation between the impact of solar events at Mars and the appearance of SEP aurora is well established (14, 21), but individual events have not previously been targeted in advance for the purpose of auroral detection. Gradual, proton-dominated ICME SEP events are ideal for this purpose; we therefore designed our observation strategy around targeting diffuse SEP aurora driven by ICME-related events.

All previous aurora observations have been made from orbit and have led to the identification of several emission lines in the ultra-violet (UV) domain. The dominant lines are the CO Cameron, CO₂⁺ UV doublet, the Fox-Duffendack-Barker (FDB) bands, and the far-UV oxygen emissions at 130 and 135 nm. The less intense atomic oxygen O(¹S-³P) line at 297.2 nm is also often detected (12, 14, 16, 18). It shares the same upper state as the 557.7-nm green line [O(¹S → ¹D)], which has been observed in photodissociation-induced martian dayglow (35, 36). The 557.7-nm line is predicted to be present in aurorae (18, 37) but has so far never been observed on Mars' nightside and has never been in response to space weather events at Mars. On Venus, it has been found in nightside emissions (38), and the appearance correlated well with the impact of solar wind disturbances (39). Fainter blue and red aurora on Mars produced from the visible end of the FDB (around 400 nm) bands and the O(¹D → ³P) atomic oxygen doublet (630 to 636.4 nm), respectively, were first predicted by Lilensten et al. (37). At that time, only discrete suprathermal aurora had been detected, so they (37) estimated an emission altitude of around 140 km, coinciding with the observed UV suprathermal aurora. The visible diffuse SEP aurora targeted here is expected to exhibit an altitude distribution similar to that of the UV SEP aurora reported in (14, 21).

Despite the large number of visible-wavelength instruments at Mars, few have the necessary combination of high spectral resolution, low-light sensitivity, and scheduling flexibility to reliably detect transient auroral events. Surface assets are more limited in terms of power, data transfer, and geographic coverage compared to orbiters, but the dynamic nature of rover planning and operations allows for the implementation of a more reactive observation strategy that can take advantage of near-real-time propagation analysis of ICME events. Starting in May 2023, we made several attempts to react to

such simulated SEP events and to observe with the Mars 2020 Perseverance rover instruments at times when the likelihood of auroral emission was highest. The fourth such attempt, in March 2024, yielded the first detection of visible-wavelength aurorae. We report on that detection here.

## INSTRUMENTS AND METHODS

### Solar storm selection

Several solar wind parameters can be used to gauge the local space weather conditions at Mars. Sudden increases in ion density and temperature can indicate the presence of ICME-driven SEPs (40), and IMF strength and the solar wind dynamic pressure are considered to be important factors and possible indicators of imminent auroral activity (10, 32, 41, 42). To trigger diffuse aurora, kilo-electron volt-mega-electron volt SEPs are required (14, 33, 34). We focus here on gradual SEP events, as impulsive SEPs arrive at Mars too fast to plan responsive observations and since CMEs are routinely analyzed and simulated by the Moon to Mars (M2M) Space Weather Analysis Office (https://science.gsfc.nasa.gov/674/m2m/) at the NASA Goddard Space Flight Center, supported by the heliospheric modeling carried out by the NASA Community Coordinated Modeling Center (CCMC). This publicly available information is crucial in the planning of these targeted observations.

We use the NASA CCMC Space Weather Database Of Notifications, Knowledge, Information (DONKI) (https://kauai.ccmc.gsfc.nasa.gov/DONKI) to evaluate ongoing Mars-directed ICME events. DONKI provides near-real-time simulation results generated from the solar corona-solar wind model, Wang-Sheeley-Arge (WSA)-ENLIL (43-45), together with the geometric “cone” CME tool (46). The WSA-ENLIL model provides a time-dependent description of the background solar wind and interplanetary magnetic field from the Sun out to 2 astronomical units (AU), while the cone tool provides estimates for the initial CME parameters (angular width, initial speed at 21.5 solar radii, and solar source latitude and longitude). To plan ahead for SEP aurora observation campaigns, we monitor Mars-facing solar active regions and routinely review DONKI for new Mars-directed ICME event simulations. For the particular event reported here, we were notified of the ICME activity through the Mars Space Weather Alert Notification email distributed to several operational Mars mission projects, courtesy of the Mars Atmosphere and Volatile Evolution (MAVEN) mission.

To initiate an aurora rover observation, we evaluate the simulated ICME properties from the CCMC simulation and target ICME events that fulfill certain requirements. These requirements were initially based on statistical solar wind properties at Mars (47) and by studying previous Mars-impacting ICMEs. After two unsuccessful attempts, the criteria were refined to increase the chances of auroral detection. The criteria used for this detection (visualized in Fig. 1) were as follows: an initial CME ejection velocity > 800 km/s estimated from white light coronagraph imaging, modeled values of ion density > 15 ions/cm³, IMF strength > 10 nT, and plasma temperature > 2 × 10⁵ K at Mars.

When an ICME fulfilling the criteria listed above is identified, we attempt to schedule an observation with the Perseverance rover's SuperCam and Mastcam-Z instruments as close in time as possible after the predicted peak time of the solar wind ion density. Note that it is not always possible to schedule an observation for a selected ICME due to a combination of factors: limited rover resources,
<figure><img src="image_3.png" /><figcaption>**Fig. 1.** Simulation of the solar wind plasma and IMF enhancement at Mars during the ICME event. The figure is adapted from simulation results generated by the NASA CCMC using the WSA-ENLIL + cone model. The left figure shows the solar wind ion density and velocity in the top and the IMF magnitude and plasma temperature in the bottom, with horizontal dotted lines marking the corresponding solar event selection criteria as described in the "Solar storm selection" section. The right figure shows a time snapshot of the modeled solar wind ion density (scaled by solar distance squared) in the ecliptic plane from 0.1 to 2 AU at the time of peak ion density at Mars, coinciding with the maximum compression of the solar wind at 1.5 AU by the ICME structure. The full simulation and solar activity details are available at https://kauai.ccmc.gsfc.nasa.gov/DONKI/view/WSA-ENLIL/29602/1.</figcaption></figure>

timing of the event simulations posted to DONKI, timing of when or whether the Mars space weather alert notifications are issued, or the estimated ICME arrival time at Mars might align poorly with the martian day/night and rover planning cycles.

## Observation strategy

The Mars 2020 Perseverance rover (48) landed in Jezero crater on 18 February 2021, at 18.5° north and 77° east longitude. For time-keeping purposes, each martian day is numbered, with mission sol 0 being the day of landing. Of its seven instruments, SuperCam and Mastcam-Z have the best combination of sensitivity and operational flexibility to detect diffuse aurora. The Radiation and Dust Sensor (RDS) (49), which is part of the Mars Environmental Dynamics Analyzer (MEDA) instrument suite (50), covers the wavelengths of UV auroral lines, but since RDS sensitivity was optimized for solar flux, its chances of detecting the weaker UV emission are low. The rover observation strategy includes acquisition of visible spectra with SuperCam and color images with Mastcam-Z, with instrument settings tailored to our assumptions on the appearance of diffuse aurora from the surface.

## Perseverance SuperCam

SuperCam is a multitechnique spectrometer designed for laser-induced breakdown (LIBS) and Raman spectroscopy and visible-infrared reflectance spectroscopy (51, 52). Only the passive mode at wavelengths below 1 µm is relevant to aurora observations. Super-Cam's visible-range spectrometer covers the 535- to 853-nm range, with a spectral resolution of 0.35 nm full width at half maximum (FWHM) in the green spectral range and 0.6 nm FWHM in the "orange" channel (which includes the 630-nm red emission line). SuperCam's 535- to 853-nm range includes an intensifier which can amplify weak optical signals and therefore makes visible-wavelength auroral emissions detectable above the dark noise from charge-coupled device (CCD) thermal electrons. SuperCam covers several

hypothesized auroral lines, such as the forbidden atomic oxygen transitions at 557.7 and 630 nm. However, the red (630 nm) line is estimated to be roughly 20 times fainter than the green (557.7 nm) line (37) and consequently more difficult to detect.

The SuperCam aurora observations consist of sets of 75 spectra each with 0.75-s integration time, for which means and medians but not individual spectra are downlinked. The complete observation for the data presented here are four sets of spectra: a "dark" set, followed by two illuminated sets, followed by a final dark set, with a total duration of 5 min. Dark sets are those which have the intensifier switched off to measure the detector background.

To estimate uncertainties, we perturb the best fit in-band emission line spectrum, which contains 32 data points surrounding the theoretically expected line center, by all possible sets of 32 adjacent residual data points in their original order drawn from the full set of residuals, repeating the line fit with each perturbation and thereby building up a probability distribution of fit results. The full set of residuals contains the 266 data points of residuals generated by the continuum plus emission line fit. Maintaining perturbing data points in their original order and adjacency ensures that the effects of autocorrelations in instrument noise are taken into account in the probability distribution.

## Perseverance Mastcam-Z

Mastcam-Z is a pair of multispectral, stereoscopic imaging cameras capable of providing broadband and narrowband color images and direct solar images using neutral density filters. It has a selectable field of view (FOV) ranging from 6.2° × 4.6° to 25.6° × 19.2° (53, 54). Each camera has a CCD detector with 1648 × 1200 photoactive pixels and Bayer pattern RGB (red, green and blue) microfilters. Bayer patterns have 4-pixel unit cells, each with two green (544 ± 42 nm) pixels, one red (630 ± 43 nm), and one blue (480 ± 46 nm). This arrangement of unit cells is noteworthy for the processing strategy described in the "Sky brightness" section. The Mastcam-Z aurora
observation included four stereo pairs taken above the western horizon using 120-s exposures, spanning roughly 9 min in total duration, immediately after the SuperCam observation. The first and last were acquired with solar filters ($10^5$ and $10^6$ neutral density for the right and left cameras, respectively) to serve as contemporaneous dark measurements, while the middle two were taken with clear filters. Mastcam-Z was at 26-mm zoom, giving a FOV of $25.6^\circ \times 19.2^\circ$ and using the full resolution of $1608 \times 1200$ pixels plus 40 virtual dark columns.

## Mars Express memory errors

To investigate whether SEPs affected Mars close to the time of observation, we used an error detection and correction (EDAC) counter from the Mars Express (MEx) orbiter, which was made available in near real time. EDACs cumulatively count bit flips in spacecraft memories (55). Bit flips can be caused by the impact of charged energetic particles (solar wind particles with energies < 10 keV do not trigger an EDAC count), thus the solar cycle modulation of galactic cosmic rays is discernible in EDAC counters as a sinusoidal curve (56), while SEP events manifest as sudden increases (57, 58). The necessary energy required to increase the EDAC counter depends on the spacecraft structure and material and will therefore differ from spacecraft to spacecraft (59). Sánchez-Cano et al. (57) found that protons with energies > 20 MeV are required to trigger an EDAC increment on MEx, which correspond to the proton energies found in (33) to contribute to the diffuse aurora.

## MAVEN SEP

To further study the energetic particles impacting Mars during the rover observation, the SEP instrument on the MAVEN satellite was used. The MAVEN/SEP instrument measures the energy spectrum and angular distribution of solar energetic electrons and ions and consists of two sensors mounted on the corners of the spacecraft, each with a FOV of by $42^\circ \times 31^\circ$. Each MAVEN/SEP sensor is a double-ended telescope that contains a solid-state detector that measures fluxes of electrons from 20 to 1000 keV and ions (mostly protons) from 20 to 6000 keV (60), thus covering the range of electron energies modeled to reproduce the vertical auroral emission

profiles observed by the MAVEN Imaging Ultraviolet Spectrograph (14, 21).

## RESULTS

### The solar storm

The analysis of the solar activity observations and the WSA-ENLIL + cone simulation results that is relevant to the Mars aurora activity reported here was performed by the NASA M2M Space Weather Analysis Office and the CCMC, respectively, and can be found at the DONKI simulation page and the web links within (https://kauai.ccmc.gsfc.nasa.gov/DONKI/view/WSA-ENLIL/29602/1). Here, we briefly provide some of the analysis details and simulation results for space weather context.

On 15 March 2024, a long-lasting C4.9 flare with an accompanying CME erupted from Active Region 13599. The CME was initially observed from the SOHO LASCO C2 coronagraph at 02:10 UTC, while the flare was observed by GOES to peak at 03:57 UTC. On the basis of initial CME observations, the ICME was parametrized with the cone tool to have an initial velocity of 1121 km/s. The modeled arrival time of the ICME leading edge was estimated to reach Mars on 18 March at 00:51 UTC, $\pm$7 hours, with the ion density peak estimated to occur around 04:00 UTC. The modeled solar wind plasma and magnetic field values substantially exceeded the solar storm selection criteria described in the “Solar storm selection” section, as shown in Fig. 1.

Figure 2 shows the impact of energetic particles at Mars. There was no increase in MEx errors or MAVEN/SEP flux immediately after 15 March, indicating that few SEPs affected Mars as a result of the flare. MEx registered an increase in errors late on 17 March, about 6 hours ahead of the MAVEN/SEP flux increase. The difference in peak times between the EDAC and MAVEN/SEP could be explained by EDACs likely being triggered by higher-energy particles (mega-electron volts), reaching Mars ahead of the bulk of the ICME which consists of less the energetic particles detected by MAVEN (kilo-electron volts). The forecast ion density peak corresponds very well with the MAVEN/SEP flux increase. MAVEN/SEP

<figure><img src="image_4.png" /><figcaption>**Fig. 2.** Energetic solar particles affecting Mars during the 18 March event. Left: MEx memory error time series. The orange dotted line represents the number of errors per day. Right: MAVEN/SEP electron and ion flux timeline. The color bars represent ion (top) and electron (bottom) fluxes. For all figures, the black dashed line indicates the time of the rover observations.</figcaption></figure>
shows a simultaneous broadband increase in SEP electron and ion, similar to the SEP signatures of previous ICMEs affecting Mars (61).

## The rover observations

The rover activity was initiated on 18 March (mission sol 1094), with the instruments switching on at 06:30 UTC, corresponding to 00:15 local mean solar time (LMST). Diffuse aurora is assumed to be spatially uniform, and the intensity recorded at the surface is influenced by viewing elevation angle and aerosol absorption. The observations presented here were made in the Northern Hemisphere early fall at $L_s = 219^\circ$, during Mars' dusty season. Because of recent regional dust storm activity, visible-band vertical-column dust opacity ($\tau$) was 1.1, slightly higher than seasonal normal, as measured by routine direct imaging of the Sun during the day (62). On the basis of preliminary assessments of stellar photometry and aerosol-scattered starlight, the true nighttime opacity may have been notably higher.

SuperCam was pointed at 20° elevation angle, to maximize the received auroral intensity. The aurora is expected to be weak and spatially uniform, so by pointing the spectrometer at a certain elevation angle, we intersect a longer path length through the emitting layer compared to if the instrument was pointed at zenith. The optimal angle depends on the dust loading and the altitude of the emitting layer. Mastcam-Z was pointed at a 15.5° elevation angle to capture the local topography at the bottom edge of the image as a visual contrast to the sky illumination. The center azimuth angle for both instruments was 252° relative to local north.

An illuminated Phobos was in the sky during the observation, which led to an increase in sky brightness. To better characterize the background night sky with no auroral emission, a follow-up observation was conducted on sol 1107 with identical sky pointing, Phobos at a similar position, comparable dust loading, and additional images acquired while Phobos was fully eclipsed.

## Spectral line identification

The SuperCam spectra were acquired over a period of 5 min beginning at 00:32 LMST (06:46 UTC). The data are processed as follows, starting from engineering data records stored in NASA's Planetary Data System (PDS): (i) take the median at each pixel position for each acquired 75 spectra set yielding one spectrum for each set; (ii) subtract the nearest dark spectrum from each illuminated spectrum; (iii) take the mean of the illuminated sets; (iv) perform a 3-pixel wide median filter to reduce signal spikes.

We use the same instrument response function as is used to calibrate the PDS SuperCam calibrated data records. This PDS instrument response function is as reported in (63) but with an additional scaling factor of 1.25 that applied when detector temperatures are below 0°C. This scaling factor is based on prelaunch laboratory testing of the SuperCam body unit (52). We use a wavelength calibration with pixel wavelength assignment from (63), which used a dataset with detector temperatures nearly identical to those of the sol 1094 observation. Note that routine LIBS observations have observed no temporal or temperature-related change in wavelength calibration, consistent with expectations from (52).

We fit a third-order polynomial continuum and then a Gaussian emission line shape to these data, with line center and integrated line intensity as the free parameters. For the width of the emission line, we use the instrumental line width of 0.35 nm FWHM from (52). The data and fit are shown in Fig. 3. With this fit, we obtain a line intensity of $101 \pm 14$ (1σ uncertainty) Rayleighs ($R$) and a line center at $557.843 \pm 0.024$ nm ($1R = 7.96 \times 10^8$ photons m$^{-2}$ sr$^{-1}$ s$^{-1}$).

<figure><img src="image_5.png" /><figcaption>**Fig. 3. SuperCam average spectrum and best-fit model.** The average of 75 × 2 spectra from the sol 1094 aurora detection observation shown in black, with the corresponding best-fit model for the continuum and the 557.7 nm emission in green. Observed spectra from nondetections on sols 790, 900, and 1107 are shown for comparison in gray, with the signal levels offset so that their averages between 559 and 566 nm are matched to that of sol 1094.</figcaption></figure>

The well-known "557.7-nm" O(1S → 1D) line has a vacuum wavelength of 557.889 nm (64), which is within 2σ of what we obtain from the observation. The vacuum wavelength is the reference for the SuperCam wavelength calibration.

To validate the above results, we performed identical data processing procedures for two previous detection attempts and the follow-up observation (sols 790, 900, and 1107), all of which yielded nondetections. The signal levels and uncertainties for all four observations are plotted in Fig. 4 for comparison with the corresponding Mastcam-Z measurements.

## Sky brightness

Mastcam-Z observations started at 00:43 LMST (06:57 UTC). Analysis of the clear filter images began with removing the dark current determined from the before-and-after dark images, where detector self-heating resulted in a cold and warm measurement. Then, all pixels of $\ge 5\sigma$ from the mean were masked and discarded from further consideration; these were largely radiation strikes but also included pixels with variable dark current. The remaining pixels were divided by an instrumental flat field and binned spatially to increase precision at the expense of resolution.

Phobos was fully illuminated during the observation, and scattering of its light by aerosols contributed to 0.6 to 1.9 digital numbers (DN) of the signal. The aerosol-scattered Phobos light was red-enhanced, causing the overall night sky color, at that time, to be very similar to a much less bright daytime martian sky. The Mastcam-Z signal from aerosol-scattered starlight (after excluding Phobos contribution) was found using previous dark sky images (e.g., sol 900) and was estimated to be ~0.5 DN or ~1.2 $\times 10^{-6}$ Wm$^{-2}$ µm$^{-1}$ sr$^{-1}$ [table 7 of (54)] across the visible spectrum. This value is presumably dependent on aerosol opacity, especially at shorter wavelengths where aerosol dust has substantial absorption. The scattered starlight has a similar signal in each color channel (thus gray in color). Any residual did not contribute to the excess green signal (positive or negative) and was therefore not removed but might be of interest for future studies of nighttime atmospheric phenomena.

To create an image of the aurora-illuminated martian night sky with the Phobos contribution removed, a Phobos illumination model was made with the discrete-ordinate-method radiative transfer model (65) using Hapke parameters from (66) and the shape model from (67). The resulting irradiance was verified against a Phobos image from sol 319. The model was calculated and sampled at each pixel and then averaged in the same way as the observations.
<figure><img src="image_6.png" /><figcaption>**Fig. 4.** Evidence of excess green signal in the Martian sky on March 18th. Profiles of excess green signal for all four aurora detection attempts (**left**), and images of the martian sky with aurora on sol 1094 (**middle**) and on a reference sol without aurora (**right**). Both images, which have identical color stretches, aerosol-scattered Phobos-light removal, and 2-pixel Gaussian smoothing, show the martian night sky contrasted against nearby topography seen at the bottom edge. See the Supplementary Materials for further image processing details. In each image, there are three squares which show how different noise-free uniform signals would appear if measured by Mastcam-Z and presented in the color stretch used for the images. The top squares represent the 557.7-nm aurora with no other sky illumination. The middle squares show the observed average signal in the center of each image after Phobos signal removal. The bottom squares show the average signals in the center of the images, with the excess green aurora signal removed for the aurora image and added for the aurora-free image. Thus, the bottom squares illustrate the colors that an aurora-free sky would have had on sol 1094 for the left image and that an aurora-illuminated sky would have had on sol 1108 for the right image. Left: The average excess green signal (calculated as described in the text) in the left image as a function of elevation angle. Mastcam-Z profiles and the model are shown as lines, while the SuperCam radiance measurements are indicated by diamonds. The colors represent different mission sols. Only sol 1094 (solid green) yielded a positive detection. The shaded green and gray areas represent the Mastcam-Z instrumental uncertainty of the best fit and the 95% confidence interval including uncertainties due to corrections for scattered Phobos light, respectively. The orange dashed line shows an auroral line radiative transfer model calculation fit to the sol 1094 SuperCam measurement.</figcaption></figure>

Corrected versions of the images were calculated by subtracting the modeled Phobos light. The two corrected pairs were then averaged, and the result can be seen in the left of Fig. 4. The corrected image includes aurora, scattered starlight, and any residual. The diffuse aurora covers the entire FOV of Mastcam-Z and causes the green-yellow hue of the image. There is no evidence for spatial nonuniformity of the aurora, and the brightness changes only with a viewing angle. With an integration time of 2 min, any temporal variability on shorter timescales is not detectable but cannot be ruled out on the basis of this observation. Deimos appears by coincidence in the upper right of aurora-free image. Because of the long exposure time, two stars create visible trails in the image.

To isolate the excess green signal and display it as a function of elevation angle, the RGB Bayer pattern array was separated into four color channels: two for green and one each for red and blue. A weighted average of the red and blue channels was created, which was multiplied by 1.05 to account for the higher sensitivity of the green channel (54). The amount of excess green light was determined by the difference between the green and the average of the red and blue channels, effectively removing white starlight and other residual white light. The unmasked pixels were then averaged across detector regions. Consistency along brightness profiles from all four images demonstrated that instrumental measurement uncertainty was <0.05 DN. Profiles from the first and second left and right images were averaged, creating the emission profile in the left of Fig. 4. The intensity is higher at larger elevation angles, where the path length through the dusty lower atmosphere from a high-altitude source is shorter.

With Phobos light removed and starlight observed to be spec-
trally uniform, there is no other likely source for an excess green

signal, and so we attribute it to the atomic oxygen O(¹S → ¹D) green 557.7-nm auroral emission. We therefore map the excess green DN signal to Rayleighs using the Mastcam-Z quantum efficiency at this wavelength together with the Mastcam-Z optical characteristics [c.f. figure 20 and table 1 of (54), respectively]. At a 20° elevation angle, the excess green signal is 0.119 ± 0.004 DN (1σ) or 82 ± 2 R with instrumental uncertainties only. The uncertainty contributed by the Phobos light model is about 0.02 DN or 14 R (1σ).

Two other nights with nondetections were processed similarly; sol 900 was moonless, while sol 790 included Phobos light. In these cases, and on sol 1107, the excess green was <0.01 DN as seen in Fig. 4, and the RMS residual in the elevation profiles was <0.02 DN.

MEDA RDS did not detect aurora on 18 March 2024. Assuming the aurora was spatially uniform, that the UV and visible lines had the same relative intensities as the strong event of September 2017 reported in (21), and using the 16.5 transitional probability factor (21), we estimate an integrated UV brightness of $4 \times 10^{-8} \text{ Wm}^{-2}$ in the RDS sensor covering the 250- to 319-nm range for the September 2017 event in (21), which is two to three times below the sensitivity threshold for the UV RDS sensor. The 18 March 2024 event was an order of magnitude weaker than the September 2017 event (see Discussion), so no MEDA RDS detection would be expected for either event.

To put the observed sky brightness in context, we made prelimi-
nary aurora brightness calculations using the multiple-scattering
pseudo-spherical radiative transfer code of (68). We replaced the
1270-nm O₂ (¹Δg) volume emission from (68) with the 557.7-nm
line volume emission with a vertical profile shaped to match the
relative brightness by altitude of the 297.2-nm emission observed
for the September 2017 event (21) and used aerosol column opacity
as observed for mission sol 1094 (i.e., $\tau = 1.1$). The original version of this code, developed by Smith *et al.* (69), was designed for any wavelength from 300 to 3900 nm; therefore, no other changes to the code used by Clancy *et al.* (68) were required. The results proved to be insensitive to the aerosol vertical profile and insensitive to aerosol properties other than total column opacity. Using this code, we calculate the brightness-versus-elevation angle profile in Fig. 4 (orange dashed line) and estimate that the observed 101 R intensity would have been 2.1 kR if viewed from orbit on the limb. Applying the transitional probability factor of 16.5 from (35), we estimate that the $O(^1S \rightarrow ^3P)$ line at 297.2 nm would produce a signal of 130 R on the limb, 10 times fainter than the September 2017 event (21) observed by MAVEN.

The radiative transfer model also indicates that, as expected, the same aurora emission rate would yield substantially brighter aurora as observed from the surface if atmospheric conditions were clear. For example, with $\tau = 0.3$ (near the minimum value observed by Perseverance), the 101 R would have been approximately 200 R at the brightest elevation angle (10° due to the low opacity), which is just below the threshold of what may be perceptible to the human eye [Lilensten *et al.* (37) give a threshold of 250 R]. If the solar event was as intense as the September 2017 event, the green line emission would lead to a sky brightness of about 2 kR in $\tau = 0.3$ clear conditions or 1 kR with conditions as dusty as during the current observation. Both situations would be visible to the human eye as a glow above the horizon.

## DISCUSSION

On 15 March 2024, a minor solar flare was followed by a larger ICME, modeled by NASA CCMC to produce a SEP event at Mars on 18 March. A space weather alert notification was issued, which provided enough lead time to schedule SuperCam and Mastcam-Z auroral observations in response to the expected SEP impact. SuperCam identified the 557.7-nm $O(^1S \rightarrow ^1D)$ green emission line with an intensity of $101 \pm 14$ R. An excess green signal in the Mastcam-Z image was attributed to the same emission.

We discard the possibility of dust or clouds causing the observed emission line, as aerosol scattered light would lead to a broad-band feature. We also find it unlikely that this is a nightglow phenomenon resulting from dayside ionization, as no similar signal was found in the previous detection attempts or in the follow-up observation. In addition, green martian nightglow manifests through the $O_2$ Herzberg II emission (70) and the much weaker Chamberlain system, with closest peaks at 552 and 558 nm, respectively (71). As there is no other likely source of green light, the 557.7-nm emission line is attributed to an auroral event.

The $0.12 \pm 0.02$ DN ($1\sigma$) at 20° elevation angle observed by Mastcam-Z (the same angle as SuperCam's observation) equates to $82 \pm 14$ R. This is fully consistent with the SuperCam measurement, as the cross-calibration is not expected to be perfect between the two instruments. The two measurements were also conducted some minutes apart, so temporal variations may also have affected the difference in retrieved brightness. The apparent elevation angle gradient retrieved from the Mastcam-Z image matches radiative transfer model expectations for a horizontally uniform source at ionospheric altitudes.

Considering that C-class flares are unlikely to produce SEPs (72) and that the MAVEN/SEP and MEx EDAC peak coincided with the

estimated arrival of the ICME, we conclude that the observed aurora was mainly induced by particles accelerated by the ICME-driven shock. To determine the relative contributions to the emission from kilo-electron volt electrons and mega-electron volt protons, further investigation is required. While the brightness of this event was dimmed by dust, events under better viewing conditions or more intense particle precipitation might be above the threshold for human vision and visible to future astronauts.

Studying auroras at multiple wavelengths is a powerful tool, as the ratio between emissions can be used to constrain the incoming particle energy distribution (1, 2), and provides insight into the atmospheric composition and pressure. In particular, the ratio of green to red intensity can be used as an indicator of the incoming particle energy distribution; higher energies favor the $O(^1S \rightarrow ^1D)$ 557.7-nm green line over the $O(^1D \rightarrow ^3P)$ 630- to 636-nm red doublet, as higher-energy particles penetrate deeper into the atmosphere where the emission from the $O(^1D \rightarrow ^3P)$ transition is suppressed by collisional quenching. During this particular event, no other emission lines were detected by SuperCam, but an upper limit for the 630-nm red line was derived to be ~50 R, although its expected brightness is a factor of 20 lower than the 557.7-nm line.

This detection confirms the prediction from UV observations that emissions also occur at visible wavelengths [e.g., (4, 73)] and demonstrates that auroral forecasting is possible. The presence of visible aurora opens a new avenue for the study of space weather events at Mars; visible imagers and spectrometers can be simpler and lower cost than the UV instrumentation used to date for the study of martian aurora, especially if using wide-angle optics with narrow-band filters tuned to specific visible auroral emission lines, such as the Mars Aurora and Dust Camera included in the proposed M-MATISSE mission (74).

The detection of visible-wavelength aurora from the surface also enables the potential future investigation of temporal auroral variability and spatial structures, observations not easily obtained from orbit with current instrumentation. Spatial and/or temporal variations in the emission intensity will reveal the degree of rapid dynamics in the magnetotail. Follow-up studies will include adapted observation strategies to attempt to capture any variability, as well as in-depth analysis of concurrent plasma and SEP measurements from orbit and other surface assets, and comparisons with data acquired during additional space weather events which followed this first detection.

## Supplementary Materials

This PDF file includes:

Supplementary Text

Figs. S1 and S2

References

## REFERENCES AND NOTES

1. M. H. Rees, D. Luckey, Auroral electron energy derived from ratio of spectroscopic emissions 1. Model computations. *J. Geophys. Res.* **79**, 5181–5186 (1974).

2. L. Soret, B. Hubert, J. C. Gérard, S. Jain, K. Chirakkil, R. Lillis, J. Deighan, Quantifying the electron energy of Mars aurorae through the oxygen emission brightness ratio at 130.4 and 135.6 nm. *J. Geophys. Res. Planets* **129**, e2023JE008214 (2024).

3. L. Soret, J. C. Gérard, L. Libert, V. I. Shematovich, D. V. Bisikalo, A. Stiepen, J.-L. Bertaux, SPICAM observations and modeling of Mars aurorae. *Icarus* **264**, 398–406 (2016).

4. F. González-Galindo, J.-C. Gérard, L. Soret, J.-Y. Chaufray, A. Fedorova, M. Holmstrom, F. Lefèvre, M. Á. López-Valverde, F. Montmessin, Airglow and aurora in the martian
atmosphere: Contributions by the Mars Express and ExoMars TGO missions. *Space Sci. Rev.* **220**, 42 (2024).

5. L. J. Paxton, Y. Zhang, H. Kil, R. K. Schaefer, "Exploring the upper atmosphere: Using optical remote sensing" in *Geophysical Monograph Series*, W. Wang, Y. Zhang, L. J. Paxton, Eds. (Wiley, 2021), pp. 487–522.

6. M. H. Acuña, J. E. Connerney, N. F. Ness, R. P. Lin, D. Mitchell, C. W. Carlson, J. McFadden, K. A. Anderson, H. Reme, C. Mazelle, D. Vignes, P. Wasilewski, P. Cloutier, Global distribution of crustal magnetization discovered by the Mars global surveyor MAG/ER experiment. *J. Geophys. Res.* **284**, 790–793 (1999).

7. D. A. Brain, F. Bagenal, M. H. Acuña, J. E. P. Connerney, Martian magnetic morphology: Contributions from the solar wind and crust. *J. Geophys. Res. Space Phys.* **108**, 2002JA009482 (2003).

8. D. L. Mitchell, R. J. Lillis, R. P. Lin, J. E. P. Connerney, M. H. Acuña, A global map of Mars' crustal magnetic field based on electron reflectometry. *J. Geophys. Res.* **112**, E01002 (2007).

9. B. Sánchez-Cano, M. Lester, D. J. Andrews, H. Opgenoorth, R. Lillis, F. Leblanc, C. M. Fowler, X. Fang, O. Vaisberg, M. Mayyasi, M. Holmberg, J. Guo, M. Hamrin, C. Mazelle, K. Peter, M. Pätzold, K. Stergiopoulou, C. Goetz, V. N. Ermakov, S. Shuvalov, J. A. Wild, P. L. Blely, M. Mendillo, C. Bertucci, M. Cartacci, R. Orosei, F. Chu, A. J. Kopf, Z. Girazian, M. T. Roman, Mars' plasma system. Scientific potential of coordinated multipoint missions: "The next generation". *Exp. Astron.* **54**, 641–676 (2022).

10. R. J. Lillis, J. Deighan, K. Chirakkil, S. Jain, M. Fillingim, M. Chaffin, G. Holsclaw, R. Susarla, D. Brain, H. al Matroushi, F. Lootah, H. al Mazmi, Y. Dong, N. Schneider, A. Azari, R. Ramstad, M. Nauth, Y. Ma, J. Halekas, J. Espley, S. Curry, Sinuous aurora on Mars: A link to the tail current sheet? *J. Geophys. Res. Space Phys.* **129**, e2024JA032477 (2024).

11. Y. Harada, J. S. Halekas, J. P. McFadden, J. Espley, G. A. DiBraccio, D. L. Mitchell, C. Mazelle, D. A. Brain, L. Andersson, Y. J. Ma, D. E. Larson, S. Xu, T. Hara, S. Ruhunusiri, R. Livi, M. B. Jakosky, Survey of magnetic reconnection signatures in the martian magnetotail with MAVEN: TAIL RECONNECTION AT MARS. *J. Geophys. Res. Space Phys.* **122**, 5114–5131 (2017).

12. J.-L. Beraux, F. Leblanc, O. Witasse, E. Quemerais, J. Lilensten, S. A. Stern, B. Sandel, O. Korablev, Discovery of an aurora on Mars. *Nature* **435**, 790–794 (2005).

13. K. Chirakkil, R. J. Lillis, J. Deighan, M. S. Chaffin, S. K. Jain, D. A. Brain, M. O. Fillingim, R. Susarla, G. Holsclaw, X. Fang, N. M. Schneider, H. A. Mazmi, H. A. Matroushi, M. Gacesa, N. El-Kork, E. Thiemann, J. S. Halekas, EMM EMUS observations of FUV aurora on Mars: Dependence on magnetic topology, local time, and season. *J. Geophys. Res. Planets* **129**, e2024JE008336 (2024).

14. N. M. Schneider, J. I. Deighan, S. K. Jain, A. Stiepen, A. I. F. Stewart, D. Larson, D. L. Mitchell, C. Mazelle, C. O. Lee, R. J. Lillis, J. S. Evans, D. Brain, M. H. Stevens, W. E. McClintock, M. S. Chaffin, M. Crismani, G. M. Holsclaw, F. Lefèvre, D. Y. Lo, J. T. Clarke, F. Montmessin, B. M. Jakosky, Discovery of diffuse aurora on Mars. *Science* **350**, aad0313 (2015).

15. J. Deighan, S. K. Jain, M. S. Chaffin, X. Fang, J. S. Halekas, J. T. Clarke, N. M. Schneider, A. I. F. Stewart, J. Y. Chaufray, J. S. Evans, M. H. Stevens, M. Mayyasi, A. Stiepen, M. Crismani, W. E. McClintock, G. M. Holsclaw, D. Y. Lo, F. Montmessin, F. Lefèvre, B. M. Jakosky, Discovery of a proton aurora at Mars. *Nat. Astron.* **2**, 802–807 (2018).

16. R. J. Lillis, J. Deighan, D. Brain, M. Fillingim, S. Jain, M. Chaffin, S. England, G. Holsclaw, K. Chirakkil, H. A. Matroushi, F. Lootah, H. A. Mazmi, E. Thiemann, F. Eparvier, N. Schneider, S. Curry, First synoptic images of FUV discrete aurora and discovery of sinuous aurora at Mars by EMUS. *Geophys. Res. Lett.* **49**, e2022GL099820 (2022).

17. D. A. Brain, J. S. Halekas, L. M. Peticolas, R. P. Lin, J. G. Luhmann, D. L. Mitchell, G. T. Delory, S. W. Bougher, M. H. Acuña, H. Rème, On the origin of aurorae on Mars. *Geophys. Res. Lett.* **33**, 2005GL024782 (2006).

18. L. Soret, J.-C. Gérard, N. Schneider, S. Jain, Z. Milby, B. Ritter, B. Hubert, T. Weber, Discrete aurora on Mars: Spectral properties, vertical profiles, and electron energies. *J. Geophys. Res. Space Phys.* **126**, e2021JA029495 (2021).

19. Y. Harada, Y. Fujiwara, R. J. Lillis, J. Deighan, H. Nakagawa, B. Sánchez-Cano, M. Lester, Y. Futaana, M. Holmström, R. A. Frahm, Discrete aurora and the nightside ionosphere of Mars: An EMM-MEX conjunction of FUV imaging, ionospheric radar sounding, and suprathermal electron measurements. *Earth Planets Space* **76**, 64 (2024).

20. B. M. Jakosky, J. M. Grebowsky, J. G. Luhmann, J. Connerney, F. Eparvier, R. Ergun, J. Halekas, D. Larson, P. Mahaffy, J. McFadden, D. F. Mitchell, N. Schneider, R. Zurek, S. Bougher, D. Brain, Y. J. Ma, C. Mazelle, L. Andersson, D. Andrews, D. Baird, D. Baker, J. M. Bell, M. Benna, M. Chaffin, P. Chamberlin, Y.-Y. Chaufray, J. Clarke, G. Collinson, M. Combi, F. Crary, T. Cravens, M. Crismani, S. Curry, D. Curtis, J. Deighan, G. Delory, R. Dewey, G. DiBraccio, C. Dong, Y. Dong, P. Dunn, M. Elrod, S. England, A. Eriksson, J. Espley, S. Evans, X. Fang, M. Fillingim, K. Fortier, C. M. Fowler, J. Fox, H. Gröller, S. Guzewich, T. Hara, Y. Harada, G. Holsclaw, S. K. Jain, R. Jolitz, F. Leblanc, C. O. Lee, Y. Lee, F. Lefèvre, R. Lillis, R. Livi, D. Lo, M. Mayyasi, W. McClintock, T. McEnulty, R. Modolo, F. Montmessin, M. Morooka, A. Nagy, K. Olsen, W. Peterson, A. Rahmati, S. Ruhunusiri, C. T. Russell, S. Sakai, J.-A. Sauvaud, K. Seki, M. Steckiewicz, M. Stevens, A. I. F. Stewart, A. Stiepen, S. Stone, V. Tenishev, E. Thiemann, R. Tolson, D. Toublanc, M. Vogt, T. Weber, P. Withers, T. Woods, R. Yelle, MAVEN observations of the response of Mars to an interplanetary coronal mass ejection. *Science* **350**, aad0210 (2015).

21. N. M. Schneider, S. K. Jain, J. Deighan, C. R. Nasr, D. A. Brain, D. Larson, R. Lillis, A. Rahmati, J. S. Halekas, C. O. Lee, M. S. Chaffin, A. Stiepen, M. Crismani, J. S. Evans, M. H. Stevens, D. Y. Lo, W. E. McClintock, A. I. F. Stewart, R. V. Yelle, J. T. Clarke, G. M. Holsclaw, F. Lefèvre, F. Montmessin, B. M. Jakosky, Global aurora on Mars during the September 2017 space weather event. *Geophys. Res. Lett.* **45**, 7391–7398 (2018).

22. M. Desai, J. Giacalone, Large gradual solar energetic particle events. *Living Rev. Sol. Phys.* **13**, 3 (2016).

23. D. V. Reames, Element abundances in impulsive solar energetic-particle events. *Universe* **9**, 466 (2023).

24. D. V. Reames, Particle acceleration at the Sun and in the heliosphere. *Space Sci. Rev.* **90**, 413–491 (1999).

25. M. J. Owens, R. J. Forsyth, The heliospheric magnetic field. *Living Rev. Sol. Phys* **10**, 5 (2013).

26. Y. Futaana, S. Barabash, M. Yamauchi, S. McKenna-Lawlor, R. Lundin, J. G. Luhmann, D. Brain, E. Carlsson, J.-A. Sauvaud, J. D. Winningham, R. A. Frahm, P. Wurz, M. Holmström, H. Gunell, E. Kallio, W. Baumjohann, H. Lammer, J. R. Sharber, K. C. Hsieh, H. Andersson, A. Grigoriev, K. Brinkfeldt, H. Nilsson, K. Asamura, T. L. Zhang, A. J. Coates, D. R. Linder, D. O. Kataria, C. C. Curtis, B. R. Sandel, A. Fedorov, C. Mazelle, J.-J. Thocaven, M. Grande, H. E. J. Koskinen, T. Sales, W. Schmidt, P. Riihela, J. Kozyra, N. Krupp, J. Woch, M. Fränz, E. Dubinin, S. Orsini, R. Cerulli-Irelli, A. Mura, A. Millio, M. Maggi, E. Roelof, P. Brandt, K. Szego, J. Scherrer, P. Bochsler, Mars Express and Venus Express multi-point observations of geoeffective solar flare events in December 2006. *Planet. Space Sci.* **67**, 873–880 (2008).

27. A. Posner, R. D. Strauss, Warning time analysis from SEP simulations of a two-tier REleAE system applied to Mars exploration. *Space Weather* **18**, e2019SW002354 (2020).

28. C. O. Lee, B. M. Jakosky, J. G. Luhmann, D. A. Brain, M. L. Mays, D. M. Hassler, M. Holmström, D. E. Larson, D. L. Mitchell, C. Mazelle, J. S. Halekas, Observations and impacts of the 10 September 2017 solar events at Mars: An overview and synthesis of the initial results. *Geophys. Res. Lett.* **45**, 8871–8885 (2018).

29. F. Leblanc, J. G. Luhman, R. E. Johnson, E. Chassefiere, Some expected impacts of a solar energetic particle event at Mars. *J. Geophys. Res. Space Phys.* **107**, SIA5-1-SIA5-10 (2002).

30. R. D. Jolitz, C. F. Dong, A. Rahmati, D. A. Brain, C. O. Lee, R. J. Lillis, S. M. Curry, B. M. Jakosky, Test particle model predictions of SEP electron transport and precipitation at Mars. *J. Geophys. Res. Space Phys.* **126**, e2021JA029132 (2021).

31. J. S. Halekas, D. A. Brain, J. G. Luhmann, G. A. DiBraccio, S. Ruhunusiri, Y. Harada, C. M. Fowler, D. L. Mitchell, J. E. P. Connerney, J. R. Espley, C. Mazelle, B. M. Jakosky, Flows, fields, and forces in the Mars-solar wind interaction. *J. Geophys. Res. Space Phys.* **122**, 11,320–11,341 (2017).

32. J. G. Luhmann, C. F. Dong, Y. J. Ma, S. M. Curry, S. Xu, C. O. Lee, T. Hara, J. Halekas, Y. Li, J. R. Gruesbeck, J. Espley, D. A. Brain, C. T. Russell, B. M. Jakosky, Martian magnetic storms. *J. Geophys. Res. Space Phys.* **122**, 6185–6209 (2017).

33. Y. Nakamura, N. Terada, F. Leblanc, A. Rahmati, H. Nakagawa, S. Sakai, S. Hiruba, R. Kataoka, K. Murase, Modeling of diffuse auroral emission at Mars: Contribution of MeV protons. *J. Geophys. Space Phys.* **127**, e2021JA029914 (2022).

34. J.-C. Gérard, L. Soret, V. I. Shematovich, D. V. Bisikalo, S. W. Bougher, The Mars diffuse aurora: A model of ultraviolet and visible emissions. *Icarus* **288**, 284–294 (2017).

35. J.-C. Gérard, S. Aoki, Y. Willame, L. Gkouvelis, C. Depiesse, I. R. Thomas, B. Ristic, A. C. Vandaele, F. Daerden, B. Hubert, J. Mason, M. R. Patel, J. J. López-Moreno, G. Bellucci, M. A. López-Valverde, B. Beeckman, Detection of green line emission in the dayside atmosphere of Mars from NOMAD-TGO observations. *Nat. Astron.* **4**, 1049–1052 (2020).

36. L. Soret, J. C. Gérard, S. Aoki, L. Gkouvelis, I. R. Thomas, B. Ristic, B. Hubert, Y. Willame, C. Depiesse, A. C. Vandaele, M. R. Patel, J. P. Mason, F. Daerden, J. J. López-Moreno, G. Bellucci, The Mars oxygen visible dayglow: A martian year of NOMAD/UVIS observations. *J. Geophys. Res. Planets* **127**, e2022JE007220 (2022).

37. J. Lilensten, D. Bernard, M. Barthélémy, G. Gronoff, C. Simon Wedlund, A. Opitz, Prediction of blue, red and green aurorae at Mars. *Planet. Space Sci.* **115**, 48–56 (2015).

38. T. G. Slanger, P. C. Cosby, D. L. Huestis, T. A. Bida, Discovery of the atomic oxygen green line in the Venus night airglow. *Science* **291**, 463–465 (2001).

39. C. L. Gray, N. J. Chanover, T. G. Slanger, K. Molaverdikhani, The effect of solar flares, coronal mass ejections, and solar wind streams on Venus' 5577Å oxygen green line. *Icarus* **233**, 342–347 (2014).

40. M. E. Cuesta, D. J. McComas, L. Y. Khoo, R. Bandyopadhyay, T. Sharma, M. M. Shen, J. S. Rankin, A. T. Cummings, J. R. Szalay, C. M. S. Cohen, N. A. Schwadron, R. Chhiber, F. Pecora, W. H. Matthaeus, R. A. Leske, M. L. Stevens, Correlation of coronal mass ejection shock temperature with solar energetic particle intensity. *Astrophys. J.* **964**, 114 (2024).

41. Z. Girazian, N. M. Schneider, Z. Milby, X. Fang, J. Halekas, T. Weber, S. K. Jain, J. C. Gérard, L. Soret, J. Deighan, C. O. Lee, Discrete aurora at Mars: Dependence on upstream solar wind conditions. *J. Geophys. Res. Space Phys.* **127**, e2021JA030238 (2022).
42. X. Fang, Y. Ma, N. Schneider, Z. Girazian, J. Luhmann, Z. Milby, S. Jain, Y. Dong, S. Curry, B. Jakosky, Discrete aurora on the nightside of Mars: Occurrence location and probability. *J. Geophys. Res. Space Phys.* **127**, e2021JA029716 (2022).

43. D. Odstrcil, V. J. Pizzo, Three-dimensional propagation of coronal mass ejections (CMEs) in a structured solar wind flow: 1. CME launched within the streamer belt. *J. Geophys. Res. Space Phys.* **104**, 483–492 (1999).

44. D. Odstrcil, Modeling 3-D solar wind structure. *Adv. Space Res.* **32**, 497–506 (2003).

45. C. N. Arge, J. G. Luhmann, D. Odstrcil, C. J. Schrijver, Y. Li, Stream structure and coronal sources of the solar wind during the May 12th, 1997 CME. *J. Atmos. Sol. Terr. Phys.* **66**, 1295–1309 (2004).

46. G. Millward, D. Biesecker, V. Pizzo, C. A. de Koning, An operational software tool for the analysis of coronagraph images: Determining CME parameters for input into the WSA-Enil heliospheric model. *Space Weather* **11**, 57–68 (2013).

47. R. Ramstad, S. Barabash, Y. Futaana, M. Yamauchi, H. Nilsson, M. Holmström, Mars under primordial solar wind conditions: Mars Express observations of the strongest CME detected at Mars under solar cycle #24 and its impact on atmospheric ion escape. *Geophys. Res. Lett.* **44**, 10,805–10,811 (2017).

48. K. A. Farley, K. H. Williford, K. M. Stack, L. Bhartia, A. Chen, M. de la Torre, K. Hand, Y. Goreva, D. C. K. Herd, R. Hueso, Y. Liu, J. N. Maki, G. Martínez, R. C. Moeller, A. Nelessen, C. E. Newman, D. Nunes, A. Ponce, N. Spanovich, P. A. Willis, L. W. Beegle, J. F. Bell III, A. J. Brown, S. E. Hamran, J. A. Hurowitz, S. Maurice, D. A. Paige, J. A. Rodriguez-Manfredi, M. Schulte, C. R. Wiens, Mars 2020 mission overview. *Space Sci. Rev.* **216**, 142 (2020).

49. V. Apestigue, A. Gonzalo, J. J. Jiménez, J. Boland, M. Lemmon, J. R. de Mingo, E. García-Menendez, J. Rivas, J. Azcue, L. Bastide, N. Andrés-Santiuste, J. Martínez-Oter, M. González-Guerrero, A. Martin-Ortega, D. Toledo, F. J. Alvarez-Rios, F. Serrano, B. Martín-Vodopivec, J. Manzano, R. L. Heredero, I. Carrasco, S. Aparicio, Á. Carretero, D. R. M. Donald, L. B. Moore, M. Á. Alcacer, J. A. Fernández-Viguri, I. Martín, M. Yela, M. Álvarez, P. Manzano, J. A. Martin, J. C. Del Hoyo, M. Reina, R. Urqui, J. A. Rodriguez-Manfredi, M. de la Torre Juárez, C. Hernandez, E. Cordoba, R. Leiter, A. Thompson, S. Madsen, M. D. Smith, D. Viúdez-Moreiras, A. Saiz-Lopez, A. Sánchez-Lavega, L. Gomez-Martín, G. M. Martinez, F. J. Gómez-Elvira, I. Arruego, Radiation and dust sensor for Mars environmental dynamic analyzer onboard M2020 rover. *Sensors* **22**, 2907 (2022).

50. J. A. Rodriguez-Manfredi, M. de la Torre Juárez, A. Alonso, V. Apéstigue, I. Arruego, T. Atienza, D. Banfield, J. Boland, M. A. Carrera, L. Castañer, J. Ceballos, H. Chen-Chen, A. Cobos, P. G. Conrad, E. Cordoba, T. del Rio-Gaztelurrutria, A. de Vicente-Retortillo, M. Dominguez-Pumar, S. Espejo, A. G. Fairen, A. Fernández-Palma, F. Tardíz, F. Ferri, E. Fischer, A. García-Manchado, M. García-Villadangos, M. Genzer, S. Giménez, J. Gómez-Elvira, F. Gómez, S. D. Guzewich, A. M. Harri, C. D. Hernández, M. Hieta, R. Hueso, I. Jaakonaho, J. J. Jiménez, V. Jiménez, A. Larman, R. Leiter, A. Lepinette, M. T. Lemmon, G. López, S. N. Madsen, T. Mäkinen, M. Marín, J. Martin-Soler, G. Martínez, A. Molina, L. Mora-Sotomayor, J. F. Moreno-Álvarez, S. Navarro, C. E. Newman, O. Ortega, M. C. Parrondo, V. Peinado, A. Peña, I. Pérez-Grande, S. De los Hoyos, O. Pla-García, J. Polkko, M. Postigo, O. Prieto-Ballesteros, S. C. R. Rafkin, M. Ramos, M. I. Richardson, J. Romeral, C. Romero, K. D. Runyon, A. Saiz-Lopez, A. Sánchez-Lavega, I. Sard, J. T. Schofield, E. Sebastian, M. D. Smith, R. J. Ulrich, L. K. Tamppari, A. D. Thompson, D. Toledo, F. Torrero, J. Torres, R. Urqui, T. Velasco, D. Viúdez-Moreiras, S. Zurita, The MEDA team, A suite of environmental sensors for the Mars 2020 mission. *Space Sci. Rev.* **217**, 48 (2021).

51. S. Maurice, R. C. Wiens, P. Bernardi, P. Cais, S. Robinson, T. Nelson, O. Gasnault, J. M. Reess, M. Deleuze, F. Rull, J. A. Manrique, S. Abbaki, R. B. Anderson, Y. André, S. M. Angel, G. Arana, T. Battault, P. Beck, K. Benzerara, S. Bernard, J. P. Berthias, O. Beyssac, M. Bonafous, B. Bousquet, M. Boutillier, A. Cadu, K. Castro, F. Chapron, B. Chide, K. Clark, E. Clavé, S. Clegg, E. Cloutis, C. Collin, E. C. Cordoba, A. Cousin, J. C. Dameury, W. D’Anna, Y. Daydou, A. Debus, L. Deflores, D. Ehouck, D. Depp, G. de Los Santos, C. Donny, A. Doressoundiram, G. Dromart, B. Dubois, A. Dufour, M. Dupieux, M. Egan, J. Ervin, C. Fabre, A. Fau, W. Fischer, O. Forni, T. Fouchet, J. Frydenvang, S. Gauffre, M. Gauthier, V. Gharakanian, O. Gilard, I. Gontijo, R. Gonzalez, D. Granena, S. Grotzinger, R. Hassen-Khodja, M. Heim, Y. Hello, G. Hervet, O. Humeau, X. Jacob, S. Jacquinod, J. R. Johnson, D. Kouach, G. Lacombe, N. Lanza, L. Lapauw, J. Laserna, J. Lasue, L. le Deit, S. le Mouélic, E. le Comte, Q. M. Lee, C. Legett IV, R. Leville, E. Lewin, C. Leyrat, G. Lopez-Reyes, R. Lorenz, B. Lucero, J. M. Madariaga, S. Madsen, M. Madsen, N. Mangold, F. Manni, J. F. Mariscal, J. Martinez-Frias, K. Mathieu, R. Mathon, K. P. McCabe, T. McConnochie, S. M. McLennan, J. Mekki, N. Melikechi, P. Y. Meslin, Y. Micheau, Y. Michel, J. M. Michel, D. Mimoun, A. Misra, G. Montagnac, C. Montaron, F. Montmessin, J. Moros, V. Mousset, Y. Morizet, N. Murdoch, R. T. Newell, H. Newsom, N. Nguyen Tuong, A. M. Ollila, G. Otttent, L. Oudda, L. Pares, J. Parisot, Y. Parot, R. Pérez, D. Pheav, L. Picot, P. Pilleri, C. Pilorget, P. Pinet, G. Pont, F. Poulet, C. Quantin-Nataf, B. Quertier, D. Rambaud, W. Rapin, P. Romano, L. Roucayrol, C. Royer, M. Ruellan, B. F. Sandoval, V. Sautter, M. J. Schoppers, S. Schröder, H. C. Seran, S. K. Sharma, P. Sobron, A. Soudki, A. Sridhar, D. Standardovsky, S. Storms, N. Striebig, M. Tatat, M. Toplis, I. Torre-Fdez, N. Toulemont, C. Velasco, M. Veneranda, D. Venhaus, C. Virmontois, M. Viso, P. Willis, K. W. Wong, The

SuperCam instrument suite on the Mars 2020 rover: Science objectives and mast-unit description. *Space Sci. Rev.* **217**, 47 (2021).

52. R. C. Wiens, S. Maurice, S. H. Robinson, A. E. Nelson, P. Cais, P. Bernardi, R. T. Newell, S. Clegg, S. K. Sharma, S. Storms, J. Deming, D. Beckman, A. M. Ollila, O. Gasnault, R. B. Anderson, Y. André, S. Michael Angel, G. Arana, E. Auden, P. Beck, J. Becker, K. Benzerara, S. Bernard, O. Beyssac, L. Borges, B. Bousquet, K. Boyd, M. Caffrey, J. Carlson, D. Castro, J. Celis, B. Chide, K. Clark, E. Cloutis, E. C. Cordoba, A. Cousin, M. Dale, L. Deflores, D. Delapp, M. Deleuze, M. Dirmyer, C. Donny, G. Dromart, M. George Duran, M. Egan, J. Ervin, C. Fabre, A. Fau, W. Fischer, O. Forni, T. Fouchet, J. Fresquez, J. Frydenvang, D. Gasway, I. Gontijo, J. Grotzinger, X. Jacob, S. Jacquinod, J. R. Johnson, R. A. Klisiewicz, J. Lake, N. Lanza, J. Laserna, J. Lasue, S. le Mouélic, C. Legett IV, R. Leveille, E. Lewin, G. López-Reyes, R. Lorenz, E. Lorigny, S. P. Love, B. Lucero, J. M. Madariaga, M. Madsen, S. Madsen, N. Mangold, J. A. Manrique, J. P. Martinez, J. Martinez-Frias, K. P. McCabe, T. H. McConnochie, J. M. McGlown, S. M. McLennan, N. Melikechi, P. Y. Meslin, J. M. Michel, D. Mimoun, A. Misra, G. Montagnac, F. Montmessin, V. Mousset, N. Murdoch, H. Newsom, L. A. Ott, Z. R. Ousnamer, L. Pares, Y. Parot, R. Pawluczyk, C. Glen Peterson, P. Pilleri, P. Pinet, G. Pont, F. Poulet, C. Provost, B. Quertier, H. Quinn, W. Rapin, J. M. Reess, A. H. Regan, A. L. Reyes-Newell, P. J. Romano, C. Royer, F. Rull, B. Sandoval, J. H. Sarrao, V. Sautter, M. J. Schoppers, S. Schröder, D. Seit, T. Shepherd, P. Sobron, B. Dubois, V. Sridhar, M. J. Toplis, I. Torre-Fdez, I. A. Trettel, M. Underwood, A. Valdez, J. Valdez, D. Venhaus, P. Willis, The SuperCam instrument suite on the NASA Mars 2020 rover: Body unit and combined system tests. *Space Sci. Rev.* **217**, 4 (2021).

53. J. F. Bell, J. N. Maki, G. L. Mehall, M. A. Ravine, M. A. Caplinger, Z. J. Bailey, S. Brylow, J. A. Schaffner, K. M. Kinch, M. B. Madsen, A. Winhold, A. G. Hayes, P. Corlies, C. Tate, M. Barrington, E. Cisneros, K. Paris, K. Crawford, C. Rojas, L. Mehall, J. B. Proton, N. Cluff, R. G. Deen, B. Betts, E. Cloutis, A. J. Coates, A. Colaprete, K. S. Edgett, B. L. Ehlmann, S. Fagents, J. P. Grotzinger, C. Hardgrove, K. E. Herkenhoff, B. Horgan, R. Jaumann, J. R. Johnson, M. Lemmon, G. Paar, M. Caballo-Perucha, S. Gupta, C. Traxler, F. Preusker, M. S. Rice, M. S. Robinson, N. Schmitz, R. Sullivan, M. J. Wolff, The Mars 2020 Perseverance rover mast camera zoom (Mastcam-Z) multispectral, stereoscopic imaging investigation. *Space Sci. Rev.* **217**, 24 (2021).

54. A. G. Hayes, P. Corlies, C. Tate, M. Barrington, J. F. Bell, J. N. Maki, M. Caplinger, M. Ravine, K. M. Kinch, K. Herkenhoff, B. Horgan, J. Johnson, M. Lemmon, G. Paar, M. S. Rice, E. Jensen, T. M. Kubacki, E. Cloutis, R. Deen, B. L. Ehlmann, E. Lakdawalla, R. Sullivan, A. Winhold, A. Parkinson, Z. Bailey, J. van Beek, P. Caballo-Perucha, E. Cisneros, D. Dixon, C. Donaldson, O. B. Jensen, J. Kuik, K. Lapo, A. Magee, M. Merusi, M. Mollerup, N. Scudder, C. Seeger, E. Stanish, M. Starr, M. Thompson, N. Turenne, K. Winchell, Pre-flight calibration of the Mars 2020 rover Mastcam zoom (Mastcam-Z) multispectral, stereoscopic imager. *Space Sci. Rev.* **217**, 29 (2021).

55. P. P. Shirvani, N. R. Saxena, E. J. McCluskey, Software-implemented EDAC protection against SEUs. *IEEE Trans. Reliab.* **49**, 273–284 (2000).

56. E. W. Knutsen, O. Witasse, B. Sanchez-Cano, M. Lester, R. F. Wimmer-Schweingruber, M. Denis, J. Godfrey, A. Johnstone, Galactic radiation modulation at Mars and beyond measured with EDACs on Mars Express and Rosetta. *Astron. Astrophys.* **650**, A165 (2021).

57. B. Sánchez-Cano, O. Witasse, E. W. Knutsen, D. Meggi, S. Viet, M. Lester, R. F. Wimmer-Schweingruber, M. Pinto, R. Moissi, J. Benkhoff, H. Opgenoorth, U. Auster, J. de Bruijne, P. Collins, G. de Marchi, D. Fischer, Y. Futaana, J. Godfrey, D. Heyner, M. Holmstrom, A. Johnstone, S. Joyce, D. Lakey, S. Martinez, D. Milligan, E. Montagnon, D. Müller, S. A. Livi, T. Prusti, J. Raines, I. Richter, S. Schmid, P. Schmitz, H. Svedham, M. G. G. T. Taylor, E. Tremolizzo, D. Titov, C. Wilson, S. Wood, J. Zender, Solar energetic particle events detected in the housekeeping data of the European Space Agency’s Spacecraft Flotilla in the Solar System. *Space Weather* **21**, e20235W003540 (2023).

58. S. Viet, E. W. Knutsen, F. Montmessin, O. Witasse, B. S. Cano, M. Lester, R. F. Wimmer-Schweingruber, “Detecting and investigating space weather events around Mars with EDAC counters” in *73rd International Astronautical Congress Proceedings* (International Astronautical Federation, 2022).

59. B. Zhang, S. Zhang, G. Shen, C. Tuo, X. Zhang, H. Zhang, L. Quan, C. Tian, D. Hou, P. Zhou, W. Ji, Monitor of the single event upsets and linear energy transfer of space radiation on the Beidou navigation satellites. *Open Astron.* **32**, 20220206 (2023).

60. D. E. Larson, R. J. Lillis, C. O. Lee, P. A. Dunn, K. Hatch, M. Robinson, D. Glaser, J. Chen, D. Curtis, C. Tiu, R. P. Lin, J. G. Luhmann, B. M. Jakosky, The MAVEN solar energetic particle investigation. *Space Sci. Rev.* **195**, 153–172 (2015).

61. D. Zhao, J. Guo, H. Huang, H. Lin, Y. Hong, X. Feng, J. Cui, Y. Wei, Y. Wang, Y. Feng, L. Li, L. Liu, Interplanetary coronal mass ejections from MAVEN orbital observations at Mars. *Astrophys. J.* **923**, 4 (2021).

62. J. F. Bell III, J. N. Maki, S. Alwmark, B. L. Ehlmann, S. A. Fagents, J. P. Grotzinger, S. Gupta, A. Hayes, K. E. Herkenhoff, B. H. N. Horgan, J. R. Johnson, K. B. Kinch, M. T. Lemmon, M. B. Madsen, J. I. Núñez, G. Paar, M. Rice, J. W. Rice Jr, N. Schmitz, R. Sullivan, A. Vaughan, M. J. Wolff, A. Bechtold, T. Bosak, L. E. Duflot, A. G. Fairén, B. Garczynski, R. Jaumann, M. Merusi, C. Million, E. Ravanis, D. L. Shuster, J. Simon, M. S. Clair, C. T. Sate, W. Walter, B. Weiss, A. M. Bailey, T. Bertrand, O. Beyssac, A. J. Brown, P. Caballo-Perucha, M. A. Caplinger, C. M. Caudill, F. Cary, E. Cisneros, E. A. Cloutis, N. Cluff, P. Corlies,
K. Crawford, S. Curtis, R. Deen, D. Dixon, C. Donaldson, M. Barrington, M. Ficht, S. Fleron, M. Hansen, D. Harker, R. Howson, J. Huggett, S. Jacob, E. Jensen, O. B. Jensen, M. Jodhpurkar, J. Joseph, C. Juarez, L. C. Kah, O. Kanine, J. Kristensen, T. Kubacki, K. Lapo, A. Magee, M. Maimone, G. L. Mehall, L. Mehall, J. Mollerup, D. Viúdez-Moreiras, K. Paris, K. E. Powell, F. Preusker, J. Proton, C. Rojas, D. Sallurday, K. Saxton, E. Scheller, C. H. Seeger, M. Starr, N. Stein, N. Turenne, J. Van Beek, A. G. Winhold, R. Yingling, Geological, multispectral, and meteorological imaging results from the Mars 2020 Perseverance rover in Jezero crater. *Sci. Adv.* 8, eabo4856 (2022).

63. C. Legett, R. T. Newell, A. L. Reyes-Newell, A. E. Nelson, P. Bernardi, S. C. Bender, O. Forni, D. M. Venhaus, S. M. Clegg, A. M. Ollila, P. Pilleri, V. Sridhar, S. Maurice, R. C. Wiens, Optical calibration of the SuperCam instrument body unit spectrometers. *Appl. Opt.* **61**, 2967–2974 (2022).

64. A. Kramida, Y. Ralchenko, NIST Atomic Spectra Database, NIST Standard Reference Database 78, 1999.

65. I. Laszlo, K. Stammes, W. J. Wiscombe, S.-C. Tsay, "The discrete ordinate algorithm, DISORT for radiative transfer" in *Light Scattering Reviews, Volume 11*, A. Kokhanovsky, Ed. (Springer, 2016), pp. 3–65.

66. D. P. Simonelli, M. Wisz, A. Switala, D. Adinolfi, J. Veverka, P. C. Thomas, P. Helfenstein, Photometric properties of phobos surface materials from viking images. *Icarus* **131**, 52–77 (1998).

67. C. M. Ernst, R. T. Daly, R. W. Gaskell, O. S. Barnouin, H. Nair, B. A. Hyatt, M. M. al Asad, K. K. W. Hoch, High-resolution shape models of phobos and deimos from stereophotoclinometry. *Earth Planets Space* **75**, 103 (2023).

68. R. T. Clancy, M. D. Smith, F. Lefèvre, T. H. M. Connochie, B. J. Sandor, M. J. Wolff, S. W. Lee, S. L. Murchie, A. D. Toigo, H. Nair, T. Navarro, Vertical profiles of Mars 1.27 µm O₂ dayglow from MRO CRISM limb spectra: Seasonal/global behaviors, comparisons to LMDGCM simulations, and a global definition for Mars water vapor profiles. *Icarus* **293**, 132–156 (2017).

69. M. D. Smith, M. J. Wolff, R. T. Clancy, A. Kleinböhl, S. L. Murchie, Vertical distribution of dust and water ice aerosols from CRISM Limb-geometry observations. *J. Geophys. Res. Planets* **118**, 321–334 (2013).

70. L. Soret, F. González-Galindo, J. C. Gérard, I. R. Thomas, B. Ristic, Y. Willame, A. C. Vandaele, B. Hubert, F. Lefèvre, F. Daerden, M. R. Patel, Ultraviolet NO and visible O₂ nightglow in the Mars southern winter polar region: Statistical study and model comparison. *J. Geophys. Res. Planets* **129**, e2024JE008620 (2024).

71. A. Migliorini, G. Piccioni, J. C. Gérard, L. Soret, T. G. Slanger, R. Politi, M. Snels, P. Drossart, F. Nuccilli, The characteristics of the O₂ Herzberg II and chamberlain bands observed with VIRTIS/Venus Express. *Icarus* **223**, 609–614 (2013).

72. M. Dierckxsens, K. Tziotziou, S. Dalla, I. Patsou, M. S. Marsh, N. B. Crosby, O. Malandraki, G. Tsiropoula, Relationship between solar energetic particles and properties of flares and CMEs: Statistical analysis of solar cycle 23 events. *Sol. Phys.* **290**, 841–874 (2015).

73. S. A. Haider, K. K. Mahajan, S. W. Bougher, N. M. Schneider, J. Deighan, S. K. Jain, J. C. Gérard, Observations and modeling of martian auroras. *Space Sci. Rev.* **218**, 32 (2022).

74. B. Sanchez-Cano, the M-MATISSE team, *The M-MATISSE mission: Mars Magnetosphere Atmosphere Ionosphere and Space weather SciencEAn ESA Medium class (M7) candidate*.

75. M. T. Lemmon, S. D. Guzewich, T. McConnochie, A. de Vicente-Retortillo, G. Martínez, M. D. Smith, J. F. Bell III, D. Wellington, S. Jacob, Large dust aerosol sizes seen during the 2018 martian global dust event by the Curiosity rover. *Geophys. Res. Lett.* **46**, 9448–9456 (2019).

76. M. J. Wolff, M. D. Smith, R. T. Clancy, R. Arvidson, M. Kahre, F. Seelos IV, S. Murchie, H. Savijärvi, Wavelength dependence of dust aerosol single scattering albedo as observed by the compact reconnaissance imaging spectrometer. *J. Geophys. Res.* **114**, E00D04 (2009).

**Acknowledgments:** We acknowledge the CCMC at the NASA Goddard Space Flight Center for the use of the DONKI tool, https://kaui.ccmc.gsfc.nasa.gov/DONKI/, and the MAVEN mission for providing the Mars Space Weather Alert Notification that notified us of the 18 March SEP event. This research was supported in part by NASA's Mars Exploration Program. We acknowledge M. D. Smith for providing the original code for the multiple-scattering pseudo-spherical radiative transfer code described in Discussion and M. J. Wolff for providing the aerosol scattering data sets for that code. **Funding:** E.W.K. is supported by the Research Council of Norway, grant no. 309835, Centre for Space Sensors and Systems (CENSSS), through their SFI Centre for Research-based Innovation program. T.H.M. is funded by the NASA Mars 2020 Participating Scientist Program grant 80NSSC21K0490 and California Institute of Technology Jet Propulsion Laboratory Mars 2020 Participating Scientist Program subcontract 1708435. R.C.W. acknowledges support from the NASA Mars Exploration Program through contract NNH13ZDA0180, subcontracted through LANL. L.S. is supported by the Belgian National Fund for Scientific Research (FNRS). R.J., S.C., C.O.L., and A.R. acknowledge funding through the MAVEN project funded through the NASA Mars Exploration program. For L.T. and R.F., this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration (80NM0018D004).

**Author contributions:** Conceptualization: E.W.K., T.H.M., L.S., S.M., and T.B. Methodology: E.W.K., T.H.M., R.F., M.L., J.L., O.F., A.C., S.M., T.B., and J.F.B. Software: T.H.M., O.F., R.C.W., R.F., A.R., C.O.L., and P.Pi. Validation: E.W.K., T.H.M., M.L., R.F., A.C., A.R., S.M., J.F.B., C.L., and P.Pi. Formal analysis: E.W.K., T.H.M., M.L., O.W., J.L., C.L., and A.R. Investigation: E.W.K., T.H.M., M.L., D.T., V.A., L.T., R.J., R.F., O.W., C.D., J.L., R.C.W., P.Pa., S.M., T.B., and J.F.B. Resources: T.H.M., F.M., R.F., O.W., R.C.W., A.R., and T.B. Data curation: E.W.K., T.H.M., R.F., S.C., R.J., A.R., and T.B. Writing—original draft: E.W.K., T.H.M., M.L., and P.Pa. Writing—review and editing: E.W.K., T.H.M., M.L., S.B.V., N.M.S., S.S., J.F.B., L.S., R.C.W., T.B., S.M., C.O.L., F.M., R.F., P.Pi., J.L., R.C.W., D.T., C.L., A.C., and P.Pa. Visualization: E.W.K., T.H.M., M.L., S.B.V., R.J., and T.B. Supervision: E.W.K., T.H.M., A.C., R.C.W., J.F.B., S.C., A.R., R.C.W., and T.B. Project administration: E.W.K., J.L., R.C.W., A.C., and S.M. Funding acquisition: T.H.M., R.C.W., and J.F.B. **Competing interests:** The authors declare that they have no competing interests.

**Data and materials availability:** All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Perseverance SuperCam and Mastcam-Z and MAVEN/SEP data are available at the NASA PDS. The MEx EDAC data are available at the online repository Zenodo: https://zenodo.org/records/13364069.

Submitted 7 August 2024

Accepted 3 April 2025

Published 14 May 2025

10.1126/sciadv.advs1563