Simulation of source intensity variations from atmospheric dust for solar occultation Fourier transform infrared spectroscopy at Mars

https://doi.org/10.1016/j.jms.2015.11.008Get rights and content

Highlights

Abstract

A Fourier transform spectrometer observing in solar occultation mode from orbit is ideally suited to detecting and characterizing vertical profiles of trace gases in the Martian atmosphere. This technique benefits from a long optical path length and high signal strength, and can have high spectral resolution. The Martian atmosphere is often subject to large quantities of suspended dust, which attenuates solar radiation along the line-of-sight. An instrument making solar occultation measurements scans the limb of the atmosphere continuously, and the optical path moves through layers of increasing or decreasing dust levels during a single interferogram acquisition, resulting in time-varying signal intensity. If uncorrected, source intensity variations (SIVs) can affect the relative depth of absorption lines, negatively impacting trace gas retrievals. We have simulated SIVs using synthetic spectra for the Martian atmosphere, and investigated different techniques to mitigate the effects of SIVs. We examined high-pass filters in the wavenumber domain, and smoothing methods in the optical path difference (OPD) domain, and conclude that using a convolution operator in the OPD domain can isolate the SIVs and be used to correct for it. We observe spectral residuals of less than 0.25% in both high- and low-dust conditions, and retrieved volume mixing ratio vertical profile differences on the order of 0.5–3% for several trace gases known to be present in the Martian atmosphere. These differences are smaller than those caused by adding realistic noise to the spectra. This work thus demonstrates that it should be possible to retrieve vertical profiles of trace gases in a dusty Martian atmosphere using solar occultation if the interferograms are corrected for the effects of dust.

Introduction

Solar absorption spectroscopy is affected by airborne aerosols, which absorb and scatter incoming solar radiation. These often take the form of thin clouds, water vapour, pollution, and smog. In the case of ground-based observations, these conditions may change during the day, and lead to biases in retrieved volume mixing ratios (VMRs) that may vary between measurements. While making remote sensing observations from orbit, the optical path observed by the instrument changes during acquisition, and if the line-of-sight passes through atmospheric layers with varying aerosol loading, the aerosol optical depth will change during acquisition.

The Atmospheric Chemistry Experiment Fourier Transform Spectrometer (ACE-FTS) is a high-resolution Fourier transform spectrometer (FTS) in near-polar, low-Earth orbit on the Canadian Space Agency’s (CSA’s) SCISAT, launched in 2004 [1]. It operates in solar occultation geometry, measuring the absorption of solar radiation along the atmospheric limb and yielding transmission spectra using observations of the un-occulted Sun and deep space. An ACE-FTS-like instrument would be ideally suited for detecting trace gases on Mars, where the atmospheric chemistry and the existence and distribution of trace gases are not well known. ACE-FTS has a wide spectral range (750–4400 cm−1), allowing it to search for dozens of trace gases active in the infrared. It has a spectral resolution of 0.02 cm−1, orders of magnitude better than current Mars missions [8], [14], and capable of distinguishing isotopologues. Solar occultation geometry provides very high signal-to-noise ratios (SNRs) and long optical path lengths, and allows for self-calibration between each occultation.

A challenge of applying the ACE-FTS technique to the Martian atmosphere is the presence of suspended dust particles. Dust storms occur frequently on Mars, can be global in scale, and can elevate dust to altitudes above 50 km [17], [9]. With ACE-FTS, the treatment of interference from aerosols involves the use of retrievals from altitudes with clear skies, or specific studies of cloud properties (e.g., [7], [5]) or dust events (e.g., [21], [6]). However, on Mars, the extent of the dust layers can be too large to discount, while the duration of dust events can last the majority of a proposed mission length [4], [14], so retrieval algorithms for an ACE-FTS-like instrument at Mars must be able to derive trace gas VMR vertical profiles from a dusty atmosphere.

Keppel-Aleks et al. [12] proposed a now-widely-used technique to mitigate the effects of source intensity variation (SIV) for instruments in the Total Carbon Column Observing Network (TCCON) [25]. The Greenhouse gases Observing SATellite (GOSAT) Thermal And Near infrared Sensor for carbon Observation (TANSO) FTS uses a similar technique [15]. Both techniques Fourier transform a raw interferogram, apply a high-pass filter, perform an inverse Fourier transform, and divide the raw interferogram by the filtered interferogram. This requires knowledge of the DC signal level and cannot be applied to AC-coupled interferograms, which are commonly recorded to satisfy the requirements of specific analog-to-digital converters (ADCs) used on the ground (e.g., [24], [19]) and from orbit (e.g., [1], [13]). If operating an ACE-FTS-like instrument at Mars, DC coupling will be a necessary requirement to measure and mitigate changes in the incoming solar signal.

A solar occultation instrument tracks the centre of the solar disk as the spacecraft comes out of, or enters, the shadow of the planet. During an occultation, the location and altitude of the tangent point along the optical path changes continuously. ACE-FTS uses a double pendulum swing arm with a maximum optical path difference (OPD) of ±25 cm, and interferogram acquisition takes 2 s. How many interferograms are acquired during an occultation, and the altitude spacing between them, depends on the β angle (between the orbit plane and the vector from the Sun). With ACE-FTS, β angles between ±20° result in a mean tangent altitude spacing between measurements of 5.5–6 km above 20 km during an occultation.

On Mars, the amount of dust along the optical path can vary significantly over the altitude range tracked during a single interferogram acquisition (1–6 km, depending on β angle), especially at the boundary of a dust layer. We generated synthetic spectra to simulate Mars atmospheric conditions, transformed these spectra into interferograms, and added DC signals. To simulate continuous acquisition, each interferogram was perturbed using the interferograms and DC levels of the measurements from the previous and next tangent height. We then investigated three methods to recover transmission spectra and compared them to the original synthetic spectra.

In Section 2, we describe the creation of synthetic spectra for the Mars atmosphere, their transformation into interferograms, and the SIV perturbation applied. In Section 3, we present the SIV mitigation strategies we investigated, and in Sections 4 Results, 5 Discussion we discuss comparisons of spectra and gas retrievals between the original synthetic spectra and those influenced by SIVs.

Section snippets

Simulated spectra

Synthetic transmission spectra, with a range of 850–4320 cm−1 and resolution of 0.02 cm−1 were generated using the GGG software suite used for analysis of spectra from the MkIV balloon-borne FTS [22] and TCCON [25]. The full spectral range is divided into two channels representing an HgCdTe (MCT) detector between 850 and 2000 cm−1, and an InSb detector between 1900 and 4320 cm−1. A priori profiles were developed at NASA’s Jet Propulsion Laboratory (JPL), based on Viking mission results [11], [16],

Mitigation

The strategy to mitigate SIVs, originally suggested by Brault [3], is to obtain a smooth function with which to re-weight the interferogram as:Icorr(x)=ISIV(x)Ismooth(x).The resulting Icorr(x) will have a constant DC level of 1, preserving spectral information in the centreburst and high-OPD wings, but requiring re-normalization. We examined three methods to obtain Ismooth(x):

  • (i)

    using the known V(x) as Ismooth(x),

  • (ii)

    obtaining Ismooth(x) by high-pass filtering ISIV(x) in the wavenumber domain as in

Spectra

The effects of the SIV perturbation are shown in Fig. 4. The top panel shows the maximum OPD regions of an interferogram with the +25 cm−1 region’s x-axis reversed. This interferogram has been corrected using only the known V(x), such that the -25 cm−1 region and the +25 cm−1 have the same DC level. This illustrates the magnitude of the asymmetry caused by perturbing the interferograms with those from altitudes above and below, simulating interferogram acquisition beginning in a different optical

Discussion

Keppel-Aleks et al. [12] distinguish between grey (absorbed equally at all wavelengths) and non-grey SIVs, correctly identifying the limitation of Method (ii) when the SIV strength is wavenumber-dependent. They applied and evaluated their filter on non-grey SIVs measured by a ground-based interferometer, then simulated low-amplitude grey SIVs that may be encountered by a ground-based FTS on Earth. They found that the correction was less effective for the grey dataset than non-grey, but still

Conclusions

A solar occultation FTS similar to the Earth-observing ACE-FTS is ideally suited to detecting unknown trace gases in the Martian atmosphere and retrieving their VMR vertical profiles. Dust storms elevate the dust content of the Martian atmosphere, scattering and absorbing transmitted solar radiation. A problem faced by an ACE-FTS-like mission to Mars is that the dust level through the optical path can vary strongly during the acquisition of each interferogram. We simulated transmission spectra

Acknowledgements

Funding for this project was provided by the CSA and the Natural Sciences and Engineering Research Council of Canada (NSERC). We would like to thank the ACE Science Team for providing Level 1 data (spectra), Level 0 data (raw interferograms), and for their help and input throughput the project. We want to thank collaborators on MATMOS and members of TCCON for guidance with GGG.

References (25)

  • P.F. Bernath et al.

    Atmospheric chemistry experiment (ACE): mission overview

    Geophys. Res. Lett.

    (2005)
  • C.D. Boone et al.

    Retrievals for the atmospheric chemistry experiment Fourier-transform spectrometer

    Appl. Opt.

    (2005)
  • J.W. Brault

    Fourier transform spectroscopy

  • B.A. Cantor et al.

    MARCI and MOC observations of the atmosphere and surface cap in the north polar region of Mars

    Icarus

    (2010)
  • J. Dodion et al.

    Cloud detection in the upper troposphere-lower stratosphere region via ACE imagers: a qualitative study

    J. Geophys. Res.

    (2007)
  • D. Doeringer et al.

    Observation of sulfate aerosols and SO2 from the Sarychev volcanic eruption using data from the atmospheric chemistry experiment (ACE)

    J. Geophys. Res.

    (2012)
  • M.N. Eremenko et al.

    Properties of high-altitude tropical cirrus clouds determined from ACE FTS observations

    Geophys. Res. Lett.

    (2005)
  • V. Formisano et al.

    The planetary Fourier spectrometer (PFS) onboard the European Mars Express mission

    Planet Space Sci.

    (2005)
  • S.D. Guzewich et al.

    High-altitude dust layers on Mars: observations with the thermal emission spectrometer

    J. Geophys. Res.

    (2013)
  • N.G. Heavens et al.

    The vertical distribution of dust in the Martian atmosphere during northern spring and summer: observations by the Mars Climate Sounder and analysis of zonal average vertical dust profiles

    J. Geophys. Res.

    (2011)
  • M.N. Izakov

    The Martian upper atmosphere structure from the Viking spacecraft experiments

    Icarus

    (1978)
  • G. Keppel-Aleks et al.

    Reducing the impact of source brightness fluctuations on spectra obtained by Fourier-transform spectrometry

    Appl. Opt.

    (2007)
  • Cited by (0)

    View full text