Composition and size of Martian aerosols as seen in the IR from solar occultation measurements by NOMAD onboard TGO
Authors: stolzenbach aurélien , López Valverde Miguel-Angel , Brines Adrian , Modak Ashimananda , Funke Bernd , González-Galindo Francisco , Thomas Ian , Liuzzi Giuliano , Villanueva Geronimo , Luginin Mikhail , Aoki Shohei 
Published: September 2022

&lt;p&gt;&lt;strong&gt;Introduction&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;The nature, size and content of aerosols in the atmosphere affect the energy budget on all planets, hence the atmospheric dynamic of the planet. Mars exhibits three types of atmospheric aerosol. Mineral dust, water ice and carbon dioxide ice. Martian aerosols nature and size distribution were observed using many different methods and experiments, from rovers to satellites. Exhaustive review scan be found in [1] and in [2]. Usually, dust effective radius, r&lt;sub&gt;eff&lt;/sub&gt;, ranges from 1 to 2 &amp;#956;m and its effective variance, &amp;#957;&lt;sub&gt;eff&lt;/sub&gt;, from 0.2 to 0.4. H&lt;sub&gt;2&lt;/sub&gt;O ice r&lt;sub&gt;eff&lt;/sub&gt; ranges from 1 to 5 &amp;#956;m and its &amp;#957;&lt;sub&gt;eff&lt;/sub&gt; from 0.1 to 0.4. However, these two parameters and their variability are poorly constraint in the vertical to date. ExoMars TGO mission (ESA/Roscosmos) was primarily designed to study trace gases, thermal structure and aerosol content in Mars atmosphere with unprecedented vertical resolution [3]. &lt;/p&gt;&lt;p&gt;&lt;strong&gt;NOMAD-SO Data processing&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;NOMAD (Nadir and Occultation for MArs Discovery) is suite of two infrared spectrometers onboard the ExoMars 2016 Trace Gas Orbiter (TGO) orbiter, covering the spectral range of 0.2 to 4.3 &amp;#956;m [4]. An Acousto-Optical Tunable Filter (AOTF) is used to select different spectral windows. The sampling of this channel is approximately of 1 second, allowing a vertical sampling about 1km. the SO channel is able to observe the atmosphere at a given altitude with 6 different diffraction orders. For this study, we selected a configuration of 5 diffraction orders (121,134,149,168,190) effectively spanning the overall spectral range of NOMAD.&lt;/p&gt;&lt;p&gt;In order to evaluate the local extinction due to aerosols, we use an inversion program called Retrieval Control Program (RCP). It is a multi-parameter non-linear least squares fitting of measured and modelled spectra [5]. Its forward model, KOPRA, was recently adapted to limb emissions on Mars [6] and for solar occultation data on Mars for the first time. RCP solves iteratively the inverse problem [7] and is described in details in [8]. The regularization matrix is build from Tikhonov-type terms of different orders which can be combined to obtain a custom-tailored regularization for any particular retrieval problem.&lt;/p&gt;&lt;p&gt;An example of the retrieved extinction profile is shown in Fig 1. The retrieved extinctions differs from previous work on aerosols using ACS data [9,10]&amp;#160; using the Onion-peeling or Abel's transform method since this global fit is less affected by the large error propagation to low altitudes typical of those methods, and the lower Martian atmosphere is precisely where aerosols are particular relevant.&lt;/p&gt;&lt;p&gt;&lt;img src=&quot;https://contentmanager.copernicus.org/fileStorageProxy.php?f=gnp.9dee3cb4b48268814682561/sdaolpUECMynit/2202CSPE&amp;app=m&amp;a=0&amp;c=c3ce34094b99922ae3b0d1b323aad8d3&amp;ct=x&amp;pn=gnp.elif&amp;d=1&quot; alt=&quot;&quot; width=&quot;486&quot; height=&quot;468&quot;&gt;&lt;/p&gt;&lt;p&gt;Fig 1.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Mean extinction cross-section ratio modelling&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;In order to model the optical behavior of the Martian aerosol we chose the log-normal distribution which is widely used in atmospheric sciences. It is a function of two parameters (r&lt;sub&gt;g&lt;/sub&gt;, &amp;#963;&lt;sub&gt;g&lt;/sub&gt;). In optics, we change those parameters to more suitable ones, the effective radius, r&lt;sub&gt;eff &lt;/sub&gt;and its corresponding effective variance &amp;#957;&lt;sub&gt;eff&lt;/sub&gt;. For any aerosol size distribution, the extinction k is km&lt;sup&gt;-1&lt;/sup&gt;&amp;#160;is k(&amp;#955;) = N . &amp;#963;&lt;sub&gt;ext&lt;/sub&gt;&amp;#160;(&amp;#955;r&lt;sub&gt;eff&lt;/sub&gt;,&amp;#957;&lt;sub&gt;eff&lt;/sub&gt;). N is the aerosol number density and &amp;#963;&lt;sub&gt;ext&lt;/sub&gt; (&amp;#955;,r&lt;sub&gt;eff&lt;/sub&gt;,&amp;#957;&lt;sub&gt;eff&lt;/sub&gt;) is the mean average extinction cross-section at a wavelength &amp;#955;, a specific aerosol distribution defined by (r&lt;sub&gt;eff&lt;/sub&gt;,&amp;#957;&lt;sub&gt;eff&lt;/sub&gt;). We build a look-up table of dust and water ice &amp;#963;&lt;sub&gt;ext&lt;/sub&gt; at the selected NOMAD order's wavelengths for different sets of (r&lt;sub&gt;eff&lt;/sub&gt;,&amp;#957;&lt;sub&gt;eff&lt;/sub&gt;). The extinction are evaluated with a Lorenz-Mie code for polydisperse spherical particle from [11].&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Aerosol composition and size distribution evaluation&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;We will detail the process of evaluating the aerosol composition and size distribution that consists of a mix of non-linear least square and brute force in order to evaluate the best set of parameters (r&lt;sub&gt;eff&lt;/sub&gt;,&amp;#957;&lt;sub&gt;eff&lt;/sub&gt; ,&amp;#947;) where &amp;#947; represent a mixture of dust and H&lt;sub&gt;2&lt;/sub&gt;O ice&lt;sub&gt;.&lt;/sub&gt; The NLSQ algorithm is provided by the SciPy Python package [12]. To assess the robustness and limitations of our evaluation procedure, we will present results against synthetic extinction signal. We will discuss our main results, especially for the period covering the Global Dust Storm of MY34 (Fig 2.).&lt;/p&gt;&lt;p&gt;&lt;img src=&quot;https://contentmanager.copernicus.org/fileStorageProxy.php?f=gnp.a8e3e9d4b48267124682561/sdaolpUECMynit/2202CSPE&amp;app=m&amp;a=0&amp;c=2a5bd4f384692b29ec745b0c68ef1ecb&amp;ct=x&amp;pn=gnp.elif&amp;d=1&quot; alt=&quot;&quot; width=&quot;801&quot; height=&quot;570&quot;&gt;&lt;/p&gt;&lt;p&gt;Fig 2.&lt;/p&gt;&lt;p&gt;&lt;strong&gt;Acknowledgments&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;The IAA/CSIC team acknowledges financial support from the State Agency for Research of the Spanish MCIU through the \emph{&quot;Center of Excellence Severo Ochoa&quot;} award for the Instituto de Astrof&amp;#237;sica de Andaluc&amp;#237;a (SEV-2017-0709) and funding by grant PGC2018-101836-B-100 (MCIU/AEI/FEDER, EU). ExoMars is a space mission of the European Space Agency (ESA) and Roscosmos. The NOMAD experiment is led by the Royal Belgian Institute for Space Aeronomy (IASB-BIRA), assisted by Co-PI teams from Spain (IAA-CSIC), Italy (INAF-IAPS), and the United Kingdom (Open University).&lt;/p&gt;&lt;p&gt;&lt;strong&gt;References&lt;/strong&gt;&lt;/p&gt;&lt;p&gt;[1] Robert M. Haberle et al., eds. The Atmosphere and Climate of Mars. Cambridge University Press, 2017.&lt;/p&gt;&lt;p&gt;[2] R. Todd Clancy et al. &amp;#8220;The distribution, composition, and particle properties of Mars meso-spheric aerosols: An analysis of CRISM visible/near-IR limb spectra with context from near-coincident MCS and MARCI observations&amp;#8221;. Icarus 328 (2019).&lt;/p&gt;&lt;p&gt;[3] J. Vago et al. &amp;#8220;ESA ExoMars program: The next step in exploring Mars&amp;#8221;. SSR 49.7 (2015).&lt;/p&gt;&lt;p&gt;[4] A. C. Vandaele et al. &amp;#8220;NOMAD, an Integrated Suite of Three Spectrometers for the ExoMarsTrace Gas Mission: Technical Description, Science Objectives and Expected Performance&amp;#8221;. SSR 214.5 (2018).&lt;/p&gt;&lt;p&gt;[5] T. von Clarmann et al. &amp;#8220;Retrieval of temperature and tangent altitude pointing from limb emission spectra recorded from space by the Michelson Interferometer for Passive Atmospheric Sounding (MIPAS)&amp;#8221;. JGR: Atmospheres 108.D23 (2003).&lt;/p&gt;&lt;p&gt;[6] Sergio Jim&amp;#233;nez-Monferrer et al. &amp;#8220;CO2 retrievals in the Mars daylight thermosphere from its 4.3&amp;#956;m limb emission measured by OMEGA/MEx&amp;#8221;. Icarus 353 (2021).&lt;/p&gt;&lt;p&gt;[7] Clive D Rodgers. Inverse Methods for Atmospheric Sounding. WORLD SCIENTIFIC, 2000.&lt;/p&gt;&lt;p&gt;[8] Jurado Navarro et al. Retrieval of CO2 and collisional parameters from the MIPAS spectra in the Earth atmosphere. Universidad de Granada, 2016.&lt;/p&gt;&lt;p&gt;[9] M. Luginin et al. &amp;#8220;Properties of Water Ice and Dust Particles in the Atmosphere of Mars During the 2018 Global Dust Storm as Inferred From the Atmospheric Chemistry Suite&amp;#8221;. JGR: Planets 125.11 (2020).&lt;/p&gt;&lt;p&gt;[10] A. Stcherbinine et al. &amp;#8220;Martian Water Ice Clouds During the 2018 Global Dust Storm as Observed by the ACS-MIR Channel Onboard the Trace Gas Orbiter&amp;#8221;. JGR: Planets 125.3 (2020).&lt;/p&gt;&lt;p&gt;[11] Michael I Mishchenko et al. Scattering, absorption, and emission of light by small particles. Cambridge university press, 2002.&lt;/p&gt;&lt;p&gt;[12] Pauli Virtanen et al. &amp;#8220;SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python&amp;#8221;. Nature Methods 17 (2020).&lt;/p&gt;