Abstract
Dust storms on Mars play a role in transporting water from its lower to upper atmosphere, seasonally enhancing hydrogen escape. However, it remains unclear how water is diurnally transported during a dust storm and how its elements, hydrogen and oxygen, are subsequently influenced in the upper atmosphere. Here, we use multi-spacecraft and space telescope observations obtained during a major dust storm in Mars Year 33 to show that hydrogen abundance in the upper atmosphere gradually increases because of water supply above an altitude of 60âkm, while oxygen abundance temporarily decreases via water ice absorption, catalytic loss, or downward transportation. Additionally, atmospheric waves modulate dust and water transportations, causing alternate oscillations of hydrogen and oxygen abundances in the upper atmosphere. If dust- and wave-driven couplings of the Martian lower and upper atmospheres are common in dust storms, with increasing escape of hydrogen, oxygen will less efficiently escape from the upper atmosphere, leading to a more oxidized atmosphere. These findings provide insights regarding Marsâ water loss history and its redox state, which are crucial for understanding the Martian habitable environment.
Similar content being viewed by others
Introduction
It is believed that, over billions of years, Mars has lost water through atmospheric escape. A classic Martian water loss mechanism is that H2 molecules that are chemically produced from water vapor in the lower atmosphere diffuse into the upper atmosphere, and their subsequent dissociation results in hydrogen escape1. As the time scale of H2 diffusion is 1200 years2, it was believed that Marsâ hydrogen escape rates have been steady within a shorter time scale. However, many observations have demonstrated that hydrogen escape rates from the Martian upper atmosphere vary seasonally due to seasonal change in vertical water vapor distribution3,4,5,6,7. Additionally, recent observations revealed that major dust storms, which expand regionally or globally on Mars, directly transfer water vapor from the lower to the upper atmosphere and rapidly change the vertical water distribution, increasing hydrogen escape5,6,7,8,9,10,11,12,13. Because major regional dust storms occur seasonally and global dust storms occasionally on Mars, monitoring Marsâ hydrogen upper atmosphere and dust storms is crucial for understanding its water escape mechanisms. Recent studies have reported that the amount of oxygen atoms likely decreases during global/regional dust storms14,15,16. This is of interest because oxygen atoms, which are mainly produced by CO2 photodissociation in the upper atmosphere, may be affected by the change in the vertical water distribution during major dust storms. This implies that some complexities exist when dust storms expand from the lower to the upper atmosphere. On this basis, comprehensive observations of hydrogen and oxygen in Marsâ upper atmosphere, and dust and water vapor in Marsâ lower atmosphere are needed.
In this paper, we analyze multi-spacecraft and space telescope observation data obtained during one of the major regional dust storms in Mars Year 33 (June 2015âMay 2017) and show that the global abundance of hydrogen atoms in the upper atmosphere gradually increases while that of oxygen atoms decreases, responding to rapid variations of lower atmospheric conditions such as air temperature, dust and water ice opacities, and water vapor mixing ratio during the dust storm. We also show a presence of alternately and periodically varying hydrogen and oxygen abundances in the upper atmosphere during the dust storm, which is likely imposed by atmospheric waves propagating from the lower atmosphere.
Results and Discussion
A regional dust storm event in 2016
In early September 2016, a regional dust storm rapidly expanded on Mars, transporting a large amount of water vapor to altitudes greater than 60 km12,17,18. During this period, we observed the upper atmosphere using the Hisaki space telescope, which carried the Extreme Ultraviolet Spectroscope for Exospheric Dynamics instrument19,20. We analyzed the HI Ly-β, OI 1304âà , and OI 1356âà emissions to examine the variability of the total amount of hydrogen and oxygen atoms in Marsâ upper atmosphere. We also obtained various measurements between the Mars surface and upper atmosphere using several instruments and spacecraft, including the Mars Climate Sounder21 (MCS) on the Mars Reconnaissance Orbiter, the Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Mars22 (SPICAM) on the Mars Express (MEX), the Extreme Ultraviolet Monitor23 (EUVM) and the Neutral Gas and Ion Mass Spectrometer24 (NGIMS) on the Mars Atmosphere and Volatile Evolution (MAVEN), and the Rover Environmental Monitoring Station25 (REMS) on the Mars Curiosity rover. These data analyses are described in detail in the Methods section.
On September 3, the regional dust storm started to expand (LS~216°), during which dust opacity and air temperature increased rapidly (~3â5 days) in the 20â60âkm altitude range (Figs. 1c, b, g, and f). Simultaneously, air pressure started to increase at all altitudes, although less rapidly than dust opacity and air temperature (Fig. 1nâp). Meanwhile, the H2O ice opacity began to decrease at the 20â40âkm altitude after storm onset, whereas it swiftly increased at altitudes of 40â60âkm (Fig. 1k, j). These observations indicate that because of the atmospheric warming caused by the lifted dust, the entire atmosphere expanded and the H2O ice cloud layer, usually present at altitudes of 10â20âkm, was also raised to a higher altitude. Near the surface (0â20âkm altitude), only small variations were observed for dust opacity, indicating an already dusty lower atmosphere. This was consistent with the fact that Mars was approaching the southern summer, during which the atmosphere is seasonally dusty, as shown in Supplementary Fig. 1.
Comprehensive observations of the middle and upper atmosphere
Figure 2 shows the time series of middle and upper atmospheric observational data obtained via multiple spacecraft during the dust storm period. Figure 2jâm show the time series of the MCS air pressure, H2O ice opacity, air temperature, and dust opacity measurements averaged at an altitude range of 60â80âkm, which are also shown in Fig. 1m, i, e, and a, respectively. In this range, according to the air pressure data, the atmosphere started to expand on the day of storm onset and increased as rapidly as at the lower altitudes (Fig. 2j). The H2O ice opacity also rapidly surged for ~5 days after the storm onset (Fig. 2k), suggesting that the H2O ice cloud layer was lifted higher than 60âkm. Meanwhile, the storm onset in dust opacity lagged a day in this altitude range, and its increasing speed, according to the peak opacity data (September 9 at 40â60âkm and September 18 at 60â80âkm), was also substantially delayed (~9 days). Figure 2n reveals the volume mixing ratio of water vapor at the four different altitude ranges observed by MEX/SPICAM. These measurements were assumed to be representative of the entire atmosphere, as they exhibited similar increasing tendencies between the northern and southern hemispheres at 20â60âkm (Supplementary Fig. 2). Overall, the volume mixing ratio gradually increased with a slow increase in dust opacity. Notably, this ratio exceeded 100 ppm at an altitude of 60â90âkm, which is a preferable condition for enhancing hydrogen escape from the upper atmosphere9,26. The rapid increase of H2O ice opacity and the slow increase of dust opacity in this altitude range may indicate that the lifted dust acted as cloud condensation nuclei. The high-altitude water vapor condensed on dust grains, and therefore, as the H2O ice opacity increased the dust opacity decreased during the formation of H2O ice clouds. The time lag in dust opacity may indicate that H2O ice clouds were being formed in this altitude range (60â80âkm).
Figure 2a shows the daily solar ultraviolet irradiances of 305âà , 1025âà , and 1305âà of the MAVEN EUVM l3 data. These data have ~5% variation during the period as a result of solar rotation and active regions on the solar surface. Figure 2g exhibits the time series of hydrogen (HI Ly-β) airglow in the upper atmosphere, as observed by the Hisaki telescope. Overall, Ly-β brightness started to increase gradually ~3 days after the storm onset and varied periodically over the period, revealing a behavior similar to the slow increase in dust opacity and water vapor mixing ratio at the 60â80âkm altitude. Note that this variation pattern was different from that of the solar 1025âà irradiance (Fig. 2a), suggesting a surge in hydrogen abundance in the upper atmosphere owing to high-altitude water vapor rather than an increased emission rate âg-factor.â The disk brightness of the airglow approximately doubled over ~20 days, indicating that the column density of hydrogen atoms at least doubled, considering that Ly-β is an optically thick emission. The optical depth of Ly-β reaches near unity at 200âkm altitude based on the H profile using the Chamberlain model27 and exospheric parameters in ref. 28. Thus, our observation most likely reflects the hydrogen column density above 200âkm altitude. The time scales of the hydrogen upper atmospheric response (~3 days) and its increase (~20 days) are comparable to those of vertical diffusion of hydrogen atoms dissociated from the high-altitude water at 60âkm to the exosphere as shown by a photochemical model26. The time scale for the hydrogen increase also corresponded to a seasonal scale of ÎLS~15°. This value was close to the time lag between the Mars perihelion (LSâ=â251°) at which airborne dust absorbs the solar flux the most (and thus dust storms would be the most common) and the nominal seasonal peak of hydrogen escape rates (LSâ~â270°)7,13,29. If the slow vertical diffusion of hydrogen atoms from an altitude of 60âkm is common for every dust storm, the slow increase in hydrogen abundance in the upper atmosphere may explain the time lag between the Mars perihelion and the seasonal hydrogen escape peak every Mars year. Another explanation for the source of hydrogen atoms is that water vapor was transported to the altitude of the ionosphere. Although no mixing ratio was retrieved above 90âkm with the SPICAM observations during the dust storm period, it is possible that water vapor was transported to the ionosphere at 100â150âkm as recently observed by MAVEN10 and the Trace Gas Orbiter4 and that the water-origin hydrogen atoms increased in the upper atmosphere through collisions with ionospheric plasma10. Because most of the hydrogen atoms produced at this altitude diffuse upward to escape Mars and its time scale is short, hydrogen airglow would vary with the water mixing ratio. Our observations show that the increasing trends of the hydrogen airglow brightness above 200âkm and water vapor mixing ratio up to 90âkm are similar to each other (Fig. 2g, n). If the water vapor was transported to the altitude of the ionosphere and exhibited a similar variation to that of the lower altitude, hydrogen atoms may originate from the 100â150âkm altitude.
The variation trend of the OI 1304 and OI 1356 brightness was similar to each other but different from that of the HI Ly-β airglow brightness (Fig 2h, i). The two oxygen emissions exhibited similar variations, likely because they mainly reflect variations of O abundance variations under the relatively small variations (~5%) of the solar 1304 à and 305 à fluxes (Fig. 2a) that affect the resonant scattering and photoelectron impact excitations of O atoms, respectively. OI 1304 is mainly excited by resonant scattering and is thus a direct indicator of the O column density, but it is an optically thick emission at Mars and its brightness is not linearly related to O column density13,30. On the other hand, because OI 1356 is an optically thin emission30, its variations mainly reflect the total abundance of O atoms in the upper atmosphere. According to the rolling average data (green line in Fig. 2), the OI 1356 brightness rapidly declined by a factor of ~3 in ~6 days after the storm onset and then recovered to the original level after another ~10 days. This suggested that the oxygen abundance in the upper atmosphere temporarily decreased by a factor of three after the dust storm. OI 1304 brightness exhibited a similar variation, but its time scale, including both decline and recovery, was shorter than that of OI 1356. This may result from OI 1356 only being excited by photoelectron impact with the peak intensity at ~130âkm, while OI 1304 was excited by the photoelectron impact at ~130âkm as well as by the resonant scattering of O atoms in a large altitude range by the solar 1304 flux30. Figure 2b shows NGIMS observations of O number density variation between 150â300âkm. Variations of the O number density at 180â190âkm, 190â200âkm, and 200â210âkm altitudes are also shown in Fig. 2câe. Although the in-situ NGIMS measurements were not directly comparable with the entire disk observations of the Hisaki telescope, the O density variations in the three altitudes overall showed a similar trend to the O airglow variations. Although periodic variations were present, the O number densities decreased by a factor of ~1.5 between the storm onset and September 12, and it subsequently increased to near the original level for ~10 days. The similarity between the in-situ and entire disk measurements indicated that the variation of O abundance was a global feature in the Martian upper atmosphere. O depletion in the upper atmosphere during a global dust storm in 2018 has previously been reported14,15, and our observation shows that even a regional dust storm globally affects the upper atmosphere of Mars.
One possible explanation for the O depletion is that rapid atmospheric heating below 60âkm caused a rapid atmospheric expansion, forming H2O ice clouds even in the thermosphere that directly absorb O atoms. On Earth, direct uptake of O atoms has been observed near noctilucent clouds in the mesosphere31,32. Indeed, our observations show that H2O ice clouds rapidly formed at 60â80âkm (Fig. 2k), and they may exist even above 100 km10. If the H2O ice clouds formed near 130âkm altitude during this dust storm period, the direct uptake of O atoms could occur. After the O absorption, the H2O ice clouds may fall and eventually sublimate in the lower atmosphere, which indicates the transportation of O atoms to the lower atmosphere. Another explanation for the O depletion is that the loss of O was accelerated through odd hydrogen catalysis in the thermosphere31,32,33,34. The source of the odd hydrogen may be photolysis of water vapor or H2O ice clouds in the thermosphere4,11,31,32. At higher altitudes, the abundance of odd hydrogen decreases, and collisions for the catalysis would become ineffective. This is consistent with a greater decrease of O abundance observed at ~130âkm by OI 1356 airglow, than at 180â210âkm by NGIMS. Another possibility is that dust storms may play a role in transporting not only water vapor but also dust-driven turbulence. High-altitude turbulence causes eddy diffusion that transports O atoms downward from the upper atmosphere34,35. In fact, the air temperature at 60â80âkm altitude began to decrease gradually ~4 days after the storm onset (Fig. 2l), which may demonstrate that CO2 molecules were transported upward by atmospheric expansion while O atoms were transported downward, resulting in a larger O/CO2 at an altitude of 60â80âkm and radiative cooling of the CO2 atmosphere36,37.
Alternating H and O abundances in the upper atmosphere
In addition to the gradual hydrogen rise and rapid oxygen decline in the upper atmosphere, the hydrogen and oxygen airglow exhibited periodic variations. Figure 3a, b show the residual brightness and variation percentages of the Ly-β and 1356 airglow with respect to their four-day rolling averages (green lines in Fig. 2g, i, respectively) as a function of time. Clear periodic variations appear both in the hydrogen and oxygen residual brightness, exhibiting an alternating relationship. Using the Lomb-Scargle periodogram method38,39, we obtained a periodicity of 6.6 days in Ly-β (Fig. 3d). A peak was observed at ~8 days for the OI 1356 emission, although this peak was smaller than the confidence level. Yet, a 6.8-day periodicity was observed in NGIMS O density data at 180â210âkm altitude (Supplementary Fig. 3). The amplitude of the airglow variations was 20-50% (Fig. 3b), suggesting that the global abundances (i.e., column density) of H and O atoms alternately varied by 20â50% every 7â8 days. Figure 3c shows the Mars Curiosity rover observations of air pressure and air temperature averaged on the dayside. Measurements demonstrate that air pressure gradually increased over the period, which was consistent with that observed by MCS. Conversely, air temperature exhibited distinct periodic variations ranging between 250âK and 270âK, demonstrating a periodicity of ~6 days, as determined via the periodogram method (Fig. 3e). These features were consistent with those of baroclinic waves usually observed in this season40,41. This periodicity was reasonably comparable to the airglow periodicity, given that the rover and space telescope measurements were based on a single point on Mars and the global view of the entire disk, respectively. This implies that the baroclinic waves propagated horizontally and vertically, modulating the atmospheric circulation to affect dust and water transportations from the lower to the upper atmosphere, which periodically affected the formation of H2O ice, odd hydrogen productions, or downward transportations and alternately changed H and O abundances in the upper atmosphere.
Our findings offer direct evidence for a coupling between the Martian lower and upper atmospheres on the dayside via a major regional dust storm and coincident atmospheric waves. The dust storm event investigated in this study played a role in increasing the hydrogen abundance in the upper atmosphere by a factor of ~2 over 20 days, but temporarily reduced oxygen abundance by a factor of 3 over 6 days. In addition, vertically propagating baroclinic waves were found to modulate dust, water, and atmospheric circulation, resulting in alternate oscillations of hydrogen and oxygen abundances in the upper atmosphere. Water ice clouds, odd hydrogen catalysis, or dust-driven turbulence as well as high-altitude water may be fundamental in causing these variations in hydrogen and oxygen atmospheres. If these are common processes for every dust storm, more hydrogen atoms may equate to fewer oxygen atoms in the upper atmosphere. This potentially stipulates that with increasing hydrogen escape, oxygen will less efficiently escape from Mars. Thus, dust storms play a role in oxidizing the atmosphere, at least on present-day Mars where dust storms seasonally occur. If dust storms have repeatedly occurred and controlled the hydrogen and oxygen escape, past Mars may have had a less oxidizing atmosphere than today, although there are many uncertainties in the history of Mars atmospheric and dust storm conditions. Hence, our findings provide insights regarding the Martian history of water loss and its redox state, which is crucial for a habitable environment.
Methods
Hisaki Space Telescope dataset
In this study, we analyzed the extreme ultraviolet (EUV) spectra obtained from the Hisaki Space Telescope19,20. The Martian upper atmosphere was observed from the end of August to the end of September 2016, during which a regional dust storm occurred. We obtained Mars EUV spectra every two days (except for September 28). Examples of daily-averaged spectra of HI Ly-β (1025âà ), OI 1304âà , and OI 1356âà airglow are shown in Supplementary Figs. 4â6. We calculated the average disk brightness and its error using the following equations:
where Ctot and CBG are count rates of the entire Mars disk and background observations (integrations of data in panels d and e in Supplementary Figs. 4â6, respectively) and ÎCtot and ÎCBG are their Poisson errors (standard deviations). G is the geometric factor converting count rates to Rayleigh. Ω and Ωd are solid angles of a single pixel and the illuminating are of the Mars disk, respectively, and their ratio is calculated as
where ÎÏ and Îθ are the angular plate scales in the spectral and slit directions, corresponding to 10âarcsec and 4.2âarcsec, respectively. The value d is the apparent diameter of Mars and p is the fraction of the illuminating area. NRM is selected to be unity for OI 1304 and OI 1356 and two for HI Ly- β to scale the size of the illuminating area. Following the above data processing method, we calculated the 2-RM (RM denotes a Mars radius) disk brightness of the hydrogen Ly-β airglow and the 1-RM disk brightness of the oxygen 1304âà and 1356âà airglow every other day. We discarded data affected by cosmic rays or EUV stars, and only used those obtained when the telescope orbited inside the shadow of Earth (defined as the local time between 20âh and 4âh) to reduce the geocoronal effect. We then employed a data selection method to increase data quality by removing pointing errors13. Finally, using these data, we studied HI Ly-β (1025âà ), OI 1304âà , and OI 1356âà airglow variations every two days during the dust storm event. Note that HI Ly-α data was not used in this study because the detectorâs sensitivity near 1216âà significantly degraded in the middle of the mission due to its constant exposure to the extremely strong Ly-α emissions of geocorona and planetary upper atmospheres such as Venus42,43,44 and Jupiter45. So, Ly-α data is currently available only at the beginning of the mission46. In-flight calibration for this issue is currently ongoing for the future use of the whole Ly-α dataset.
Mars Reconnaissance Orbiter (MRO)/Mars Climate Sounder (MCS) dataset
MRO is orbiting Mars since March 2006 to study the geology and climate of Mars47, and the MCS instrument has measured thermal emission of the Martian atmosphere21. In this study, we used a dataset of dust opacity, H2O ice opacity, air temperature, and air pressure measurements retrieved by MCS, which were stored in the MRO/MCS derived data records Version 5 in the NASA Planetary Data System. To study the dynamics of Marsâ dayside lower and middle atmospheres, these measurements were averaged within the local time of 11âh and 13âh for every 20âkm altitude between 0âkm and 80âkm every Martian day. The dataset at an altitude of 0â20âkm during Mars Year 33 is shown in Supplementary Fig. 1, and MCS geometry information during Hisakiâs observation period is shown in Supplementary Fig. 7.
Mars Atmosphere and Volatile Evolution (MAVEN)/Solar Extreme Ultravioet Monitor (EUVM) and Neutral Gas and Ion Mass Spectrometer (NGIMS) dataset
MAVEN is orbiting Mars since September 2014 to investigate the Martian atmospheric escape processes controlled by solar and atmospheric forcings48. In this study, we used the solar ultraviolet irradiance of EUVM level 3 data49, obtained from the Laboratory for Atmospheric and Space Physics Interactive Solar Irradiance Data Center. We selected the 305âà line as a proxy for photoionization of the CO2 atmosphere to produce photoelectrons in the Martian ionosphere50 and the 1025âà and 1305âà lines as sources of resonantly scattered HI Ly-β and OI 1304âà airglow13, respectively. To display the three irradiance data in the same panel in Fig. 2a, the solar 1025âà and 1305âà irradiances were scaled to the 305âà irradiance value on August 31, 2016.
We also used O number density of NGIMS level 2 data. In this study, we only used inbound measurements near periapsis (150â300âkm altitude range). Supplementary Fig. 3a shows binned (5âkm by 5âh) O densities between 150â300âkm, which is identical to Fig. 2b. Supplementary Fig. 3bâd show O density variations of three different altitude ranges (180â190âkm, 190â200âkm, and 200â210âkm), which are also identical to Fig. 2câe. To analyze their average diurnal variations, we calculated their running mean values using a 1-day window (blue dots). Supplementary Fig. 3eâh show MAVEN geometry data. We can see that MAVEN positions developed from the nightside to the dayside across the terminator in this period.
Mars Express (MEX)/Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Mars (SPICAM) dataset
MEX is orbiting Mars since December 2003 and has collected comprehensive dataset, including atmospheric water contents with the SPICAM instrument22. In this study, we used the volume mixing ratio of water vapor in the Martian atmosphere retrieved from the SPICAM measurements by ref. 18. Depending on their tangential attitude, we divided the measurements into four altitude ranges (20â40âkm, 40â60âkm, 60â80âkm, and 80â90 km), as shown in Fig. 2n. Supplementary Fig. 2 also shows the observation coverage of the SPICAM. In 20â60âkm altitude ranges, SPICAM measurements cover both northern and southern hemispheres, whereas they are only derived from the southern hemisphere in the 60â90âkm range.
Curiosity/ Rover Environmental Monitoring Station (REMS) dataset
The Curiosity rover landed in the Gale crater at the longitude of 137.4°E and latitude of 4.6°S in early August 201251. Since then, the rover has performed various experiments and observations to study habitable environment on Mars, including its meteorology with the REMS instrument25. When the Hisaki telescope observed the Martian upper atmosphere, the rover was traveling over a region called Murrays Buttes, ~8âkm southwest of the landing site, and the REMS instrument monitored atmospheric conditions. In this study, we used a dataset of air temperature and pressure of the REMS-reduced data stored in the NASA Planetary Data System. We averaged these measurements within the local time of 6âh and 18âh to focus on their dayside variations as shown in Fig. 3c.
Data availability
The Hisaki level 2 data are publicly available in the JAXA Data Archives and Transmission System (DARTS, https://data.darts.isas.jaxa.jp/pub/hisaki/euv/l2/). The Mars Atmosphere and Volatile Evolution Extreme Ultraviolet Monitor l3 data are publicly available at the LASP Interactive Solar Irradiance Data Center (https://lasp.colorado.edu/lisird/). Neutral Gas and Ion Mass Spectrometer l2 data are publicly available at MAVEN Science Data Center (https://lasp.colorado.edu/maven/sdc/public/pages/datasets/ngims.html). Mars Climate Sounder and Curiosity data are publicly available in the NASA Planetary Data System (https://atmos.nmsu.edu/data_and_services/atmospheres_data/MARS/mcs.html and https://atmos.nmsu.edu/PDS/data/mslrem_1001/, respectively). The Mars Express/Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Mars data derived from ref. 16 are stored at https://data.mendeley.com/datasets/vx4gks6bx7/1, and we would like to thank the authors for making their dataset publicly available. The Hisaki disk average airglow brightness data generated in this study is provided in Supplementary Data 1.
References
Hunten, D. M. & McElroy, M. B. Production and escape of hydrogen on Mars. J. Geophys. Res. 75, 5989â6001 (1970).
Krasnopolsky, V. A., Mumma, M. J. & Gladstone, G. R. Detection of atomic deuterium in the upper atmosphere of Mars. Science 280, 1576â1580 (1998).
Aoki, S. et al. Water vapor vertical profiles on Mars in dust storms observed by TGO/NOMAD. J. Geophys. Res. Planets 124, 3482â3497 (2019).
Belyaev, D. A. et al. Revealing a high water abundance in the upper mesosphere of Mars with ACS onboard TGO. Geophys. Res. Lett. 48, e2021GL093411 (2021).
Clarke, J. T. et al. A rapid decrease of the hydrogen corona of Mars. Geophys. Res. Lett. 41, 8013â8020 (2014).
Bhattacharyya, D. et al. Seasonal changes in hydrogen escape from mars through analysis of HST observations of the Martian exosphere near perihelion. J. Geophys. Res. Space Phys. 122, 11756â11764 (2017).
Halekas, J. S. Seasonal variability of the hydrogen exosphere of Mars. J. Geophys. Res. Planets 122, 901â911 (2017).
Chaffin, M. S. et al. Unexpected variability of Martian hydrogen escape. Geophys. Res. Lett. 41, 314â320 (2014).
Chaffin, M. S. et al. Martian water loss to space enhanced by regional dust storms. Nat. Astron. 5, 1036â1042 (2021).
Stone, S. W. et al. Hydrogen escape from Mars is driven by seasonal and dust storm transport of water. Science 370, 824â831 (2020).
Fedorova, A. A. et al. Stormy water on Mars: The distribution and saturation of atmospheric water during the dusty season. Science 367, 297â300 (2020).
Heavens, N. G. et al. Hydrogen escape from Mars enhanced by deep convection in dust storms. Nat. Astron. 2, 126â132 (2018).
Masunaga, K. et al. Martian oxygen and hydrogen upper atmospheres responding to solar and dust storm drivers: Hisaki space telescope observations. J. Geophys. Res. Planets 125, e2020JE006500 (2020).
Elrod, M. K. et al. Structural and compositional changes in the upper atmosphere related to the PEDE-2018 dust event on Mars as observed by MAVEN NGIMS. Geophys. Res. Lett. 47, e2019GL084378 (2020).
Niu, D. D. et al. Species-dependent response of the Martian ionosphere to the 2018 global dust event. J. Geophys. Res. Planets 126, e2020JE006679 (2021).
Yoshida, N. et al. Seasonal and dust-related variations of the dayside thermospheric and ionospheric compositions of Mars observed by MAVEN/NGIMS. J. Geophys. Res. Planets 126, e2021JE006926 (2021).
Wu, Z., Li, T., Zhang, X., Li, J. & Cui, J. Dust tides and rapid meridional motions in the Martian atmosphere during major dust storms. Nat. Commun. 11, 614 (2020).
Fedorova, A. et al. Multi-annual monitoring of the water vapor vertical distribution on Mars by SPICAM on Mars Express. J. Geophys. Res. Planets 126, e2020JE006616 (2021).
Yoshioka, K. et al. The extreme ultraviolet spectroscope for planetary science, EXCEED. Planet. Space Sci. 85, 250â260 (2013).
Yoshikawa, I. et al. Extreme ultraviolet Radiation measurement for planetary atmospheres/magnetospheres from the earth-orbiting spacecraft (Extreme ultraviolet Spectroscope for Exospheric Dynamics: EXCEED). Space Sci. Rev. 184, 237â258 (2014).
McCleese, D. J. et al. Mars Climate Sounder: An investigation of thermal and water vapor structure, dust and condensate distributions in the atmosphere, and energy balance of the polar regions. J. Geophys. Res. 112, E05S06 (2007).
Bertaux, J. L. et al. SPICAM on Mars Express: Observing modes and overview of UV spectrometer data and scientific results. J. Geophys. Res. Planets 111, E10S90 (2006).
Eparvier, F. G., Chamberlin, P. C., Woods, T. N. & Thiemann, E. M. B. The solar extreme ultraviolet monitor for MAVEN. Space Sci. Rev. 195, 293â301 (2015).
Mahaffy, P. R. et al. The neutral gas and ion mass spectrometer on the mars atmosphere and volatile evolution mission. Space Sci. Rev. 195, 49â73 (2015).
Gómez-Elvira, J. et al. REMS: The environmental sensor suite for the Mars Science Laboratory rover. Space Sci. Rev. 170, 583â640 (2012).
Chaffin, M. et al. Elevated atmospheric escape of atomic hydrogen from Mars induced by high-altitude water. Nat. Geosci. 10, 174â178 (2017).
Chamberlain, J. W. Planetary coronae and atmospheric evaporation. Planet. Space Sci. 11, 901â960 (1963).
Chaffin, M. S. et al. Mars H escape rates derived from MAVEN/IUVS Lyman alpha brightness measurements and their dependence on model assumptions. J. Geophys. Res. Planets 123, 2192â2210 (2018).
Mayyasi, M. et al. The variability of atmospheric deuterium brightness at Mars: Evidence for seasonal dependence. J. Geophys. Res. Space Phys. 122, 10811â10823 (2017).
Ritter, B. et al. Characteristics of Mars UV dayglow emissions from atomic oxygen at 130.4 and 135.6 nm: MAVEN/IUVS limb observations and modeling. J. Geophys. Res. Space Phys. 124, 4809â4832 (2019).
Murray, B. J. & Plane, J. M. C. Atomic oxygen depletion in the vicinity of noctilucent clouds. Adv. Space Res. 31, 2075â2084 (2003).
Murray, B. J. & Plane, J. M. C. Modelling the impact of noctilucent cloud formation on atomic oxygen and other minor constituents of the summer mesosphere. Atmos. Chem. Phys. 5, 1027â1038 (2005).
Daerden, F. et al. Planet-wide ozone destruction in the middle atmosphere on Mars during global dust storm. Geophys. Res. Lett. 49, e2022GL098821 (2022).
Shimazaki, T. Photochemical stability of CO2 in the Martian atmosphere: Reevaluation of the eddy diffusion coefficient and the role of water vapor. J. Geomag. Geoelectr. 41, 273â301 (1989).
Nara, Y. et al. Vertical coupling between the cloud-level atmosphere and the thermosphere of Venus inferred from the simultaneous observations by Hisaki and Akatsuki. J. Geophys. Res. Planets 125, e2019JE006192 (2020).
Chaufray, J. Y. et al. Martian oxygen density at the exobase deduced from OI 130.4-nm observations by Spectroscopy for the Investigation of the Characteristics of the Atmosphere of Mars on Mars Express. J. Geophys. Res. 114, E02006 (2009).
Forget, F. et al. Density and temperatures of the upper Martian atmosphere measured by stellar occultations with Mars Express SPICAM. J. Geophys. Res. 114, E01004 (2009).
Lomb, N. R. Least-squares frequency analysis of unequally spaced data. Astrophys. Space Sci. 39, 447â462 (1976).
Scargle, J. D. Studies in astronomical time series analysis. II. Statistical aspects of spectral analysis of unevenly spaced data. Astrophys. J. 263, 835â853 (1982).
Barnes, J. R. Time spectral analysis of midlatitude disturbances in the Martian atmosphere. J. Atmos. Sci. 37, 2002â2015 (1980).
Barnes, J. R. Midlatitude disturbances in the Martian atmosphere: A second Mars year. J. Atmos. Sci. 38, 225â234 (1981).
Masunaga, K. et al. Periodic variations of oxygen EUV dayglow in the upper atmosphere of Venus: Hisaki/EXCEED observations. J. Geophys. Res. Planets 120, 2037â 2052 (2015).
Masunaga, K. et al. Dawn-dusk difference of periodic oxygen EUV dayglow variations at Venus observed by Hisaki. Icaurs 292, 102â110 (2017).
Nara, Y. et al. Extreme ultraviolet spectra of Venusian airglow observed by EXCEED. Icarus 307, 207â215 (2018).
Kimura, T. et al. Transient internally driven aurora at Jupiter discovered by Hisaki and the Hubble Space Telescope. Geophys. Res. Lett. 42, 1662â1668 (2015).
Kuwabara, M. et al. The geocoronal responses to the geomagnetic disturbances. J. Geophys. Res. Space Phys. 122, 1269â1276 (2017).
Zurek, R. W. & Smrekar, S. E. An overview of the Mars Reconnaissance Orbiter (MRO) science mission. J. Geophys. Res. Planets 112, E05S01 (2007).
Jakosky, B. M., Grebowsky, J. M., Luhmann, J. G. & Brain, D. A. Initial results from the MAVEN mission to Mars. Geophys. Res. Lett. 42, 8791â8802 (2015).
Thiemann, E. M. et al. The MAVEN EUVM model of solar spectral irradiance variability at Mars: Algorithms and results. J. Geophys. Res. Space Phys. 122, 2748â2767 (2017).
Frahm, R. A. et al. Carbon dioxide photoelectron energy peaks at Mars. Icarus 182, 371â382 (2006).
Grotzinger, J. P. et al. Mars Science Laboratory Mission and Science Investigation. Space Sci. Rev. 170, 5â56 (2012).
Acknowledgements
K.M. is supported by JSPS KAKENHI (Grant numbers JP21K20387 and JP22K03708) and the JAXA Hisaki and MMX projects. N.T. is supported by JSPS KAKENHI (Grant numbers JP18H05439, JP18KK0093, JP19H00707, JP20H00192, and JP22H00164). N.Y. is supported by JSPS KAKENHI Grant number JP21J13710 and the International Joint Graduate Program in Earth and Environmental Sciences, Tohoku University (GP-EES). Y.N. is supported by JSPS KAKENHI Grant Number JP22J14954 and the International Joint Graduate Program in Earth and Environmental Sciences, Tohoku University (GP-EES).
Author information
Authors and Affiliations
Contributions
K.M. conceived the study and analyzed the datasets of the Hisaki telescope and multiple spacecraft, and wrote the paper. N.T., N.Y., Y.N., T.K., K.Y., Y.S., H.N., and T.U. contributed to the scientific interpretation of the results. K.Y., G.M., F.T., and T.K. contributed to the calibration of the Hisaki data. A.Y. led the operation of the Hisaki space telescope for Mars observation. I.Y. directed the Hisaki mission as the principal investigator. All authors participated in the discussion of the results of this paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Juan Alday and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Masunaga, K., Terada, N., Yoshida, N. et al. Alternate oscillations of Martian hydrogen and oxygen upper atmospheres during a major dust storm. Nat Commun 13, 6609 (2022). https://doi.org/10.1038/s41467-022-34224-6
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-34224-6