Mars Atmosphere: Modeling and Observations (2008) 9078.pdf

Polarimetry of Mars with SPEX, an Innovative Spectropolarimeter

Daphne M. Stam1,2 ([email protected]), Erik Laan3, Frans Snik4, Theodora Karalidi4, Christoph Keller4, Rik ter Horst5, Ramon Navarro5, Christina Aas1, Johan de Vries3, Gijs Oomen6, and Ruud Hoogeveen2, 1 DEOS, Aerospace Engineering, Technical University, Delft, 2 SRON Netherlands Institute for Space Research, Utrecht, 3 TNO Science & Industry, Delft, 4 Sterrekundig Instituut Utrecht, Utrecht, 5 ASTRON Netherlands Foundation for Research in , Dwingeloo, 6 Dutch Space, Leiden, all in the Netherlands

We present SPEX, an innovative, compact, and ro- bust Spectropolarimeter for Planetary EXploration, that measures fluxes F and P of sunlight re- flected by Mars from 400 to 800 nm with 2 nm (F ) to 20 nm (P ) spectral resolution, in 9 fixed viewing direc- tions. With SPEX, dust and ice cloud particles in the atmosphere and the surface can be studied. To illustrate the power of spectropolarimetry, we present sample sim- ulations of F and P of Mars for different dust scenarios. Introduction Polarimetry of sunlight that is reflected by a has proven to be a powerful tool to characterize the planetary atmosphere and surface. The main reason is that the state of polarization of scattered sunlight is more sensitive to the microphysical properties of particles (size, shape, and Figure 1: With 7 fixed viewing directions, SPEX will composition) than the flux [1, 2, 3]. An early example of measure flux and polarization phase functions of scenes the power of polarimetry was the derivation of the size along its ground track (the 2 limb viewers are not shown). and composition of Venus cloud particles from Earth- based polarimetry at two wavelengths [4]. Since then, have flown on the Pioneer Venus, Galileo, SPEX instrument concept and Cassini missions. The POLDER has Sunlight that has been scattered within the Mars atmo- flown on various Earth observing missions, and NASA’s sphere and/or that has been reflected by the surface can Glory mission with the Aerosol Polarimetry Sensor on- be described by a (column) vector F [see e.g. 1] board will be launched at the end of 2008. Mars, with its dust storms and ice clouds and dusty F = [F, Q, U, V ] , (1) surface, is an ideal polarimetry target. Indeed, HST- with F the total flux, Q and U the linearly polarized flux, observations show polarization variations correlating and V the circularly polarized flux. We’ll ignore V be- with clouds and dust [5]. However, these observations cause it is very small [see e.g. 8]. The degree of (linear) were done with Mars at opposition; and the degrees of polarization P of the reflected sunlight is defined as polarization are therefore very small. To truly enjoy the p advantages of polarimetry for Mars remote-sensing, a P = Q2 + U 2/F. (2) polarimeter should orbit the planet, because only then the angles that hold most of the information The polarization direction χ is given by tan 2χ = U/Q are within reach. Two polarimeters have flown onboard (defined with respect to the local meridian plane) [see 1]. the Soviet spacecraft MARS-5 [6, 7] providing informa- In planetary (including Earth) remote-sensing, fluxes tion about particle shapes, sizes, and composition, even F , Q, and U of scattered sunlight are usually determined though they encountered a very clear Mars atmosphere. from flux measurements through 3 polarization filters We present SPEX, an innovative (patented), compact with their optical axes under different angles, or through (1 liter), and robust (no moving parts) spectropolarime- a single rotating filter [see e.g. 1]. The advantage of this ter in a design tailored for observing Mars from orbit, method is its simplicity. There are, however, also dis- together with sample simulations of the observable flux advantages: 1. the time lapse between the flux measure- and degree of polarization of reflected sunlight. ments can yield polarization errors of several percents, Mars Atmosphere: Modeling and Observations (2008) 9078.pdf

2

depending on the observed scene and on the change of Sample simulations the geometries during the integration time; 2. different Both F and P (and χ) of the reflected sunlight depend on filters have different transmission properties and may de- the composition of Mars’ atmosphere, the optical prop- grade differently during flight, which will lead to errors; erties of the dust and ice cloud particles, the reflective 3. when combining color filters with polarization filters, properties of the surface, and the illumination and view- P as is usually done, you have limited information on ’s ing angles. Flux F of course also depends on the incom- spectral behaviour; and 4. filter wheels can get stuck. ing solar flux F0, hence on the distance between Mars With SPEX, fluxes F , Q and U in a given view- and the sun, but in the following, we assume F0 = π, ing direction are measured simultaneously, from 400 to independent of the wavelength. Both P and χ of the re- 800 nm. As explained by Snik et al. [9], P and χ of the flected sunlight are independent of F0, since Q, U, and reflected sunlight are encoded in the flux spectrum of the 2 F are all proportional to F0. . The encoding is a spectral modulation: the flux For our sample simulations of F and P of reflected spectrum is multiplied with a sinusoidal function. The sunlight, we describe the atmosphere as a stack of ho- modulation amplitude yields P , and the phase, χ. The mogeneous layers containing gas and dust particles. The amplitude and phase will depend on λ, like P (see the surface pressure is 6 mbars, and an average temperature next section) and χ. The advantage of the modulation is profile is chosen. We ignore the few, shallow absorp- that the full polarization information is obtained with just tion bands of CO2 in the SPEX wavelength region. The one flux spectrum measurement. surface reflects Lambertian (isotropically and depolariz- The spectral modulation is achieved with [see 9] an ing), and has a realistic, wavelength dependent albedo. achromatic quarter-wave retarder, followed by an ather- The dust particles are homogeneous palagonite spheres, mal multiple-order retarder and a polarizing beamsplit- sized according to the standard distribution of Hansen ter. The modulation frequency (i.e. the width of the mod- and Travis [1], with an effective variance of 1.5, and ef- ulations) depends on the of the com- fective radii of 0.5 or 1.0 µm. For each wavelength, the ponents, such as their refractive index and geometrical optical properties of the dust particles are calculated us- thickness. SPEX’ components have been chosen such ing Mie-theory. Although real dust particles probably [see 9] that the modulated flux is measured with a reso- have irregular shapes, spherical particles are fine to illus- lution of about 2 nm (by combining two simultaneously trate the effects of dust on the flux and polarization. We measured spectra, the flux itself can be retrieved with this choose two values for the dust optical thickness τD: 0.2 same resolution), The polarization can be retrieved with and 1.0 (at 630 nm). The dust is located in the lower part a resolution of about 20 nm. of the atmosphere (up till 1.4 mbars). SPEX has 9 apertures and measures modulated flux The simulations are performed using an adding- spectra simultaneously in 9 fixed viewing directions doubling radiative transfer code, which fully includes 1 along its flight track . Two apertures look at the limb multiple scattering and polarization [12, 13, 14]. (forward and backward) for studying ice clouds. The Figure 2 shows F and P for the 2 particle sizes and ◦ other viewing angles are approximately: 0 (nadir), the 2 dust optical thicknesses for the 7 viewing directions ◦ ◦ ◦ ±20 , ±40 , and ±55 (the + indicates the forward, of SPEX (the limb directions have been excluded). The and the − the backward direction). With these apertures, solar zenith angle is 60◦ and the sun is in the orbital plane SPEX obtains the flux and polarization phase functions of the satellite. The surface albedo strongly determines of the scenes it flies over (see Fig. 1). Each apertures ◦ F , especially when τD = 0.2 and for nadir and ±20 ◦ ◦ field-of-view is 1 × 7 (along × across track). Ground viewing angles. With increasing viewing angle and/or for pixel sizes will depend on the satellite’s orbit. τ = 1.0, F increases because of the increased scattering A provisional patent application for SPEX’ passive of light by the dust. Some sensitivity of F to the dust polarimetric method (the concept of which has been particle size is seen, especially for the largest viewing proven with a single aperture optical train in an optical angles, but the spectral shape of F changes little. laboratory) has been submitted. The project to build a P appears to be more sensitive to the particle size breadboard of SPEX with its 9 apertures has started in than F is. The variation in P with the viewing angle September 2008. The flight model will have a volume reflects the dust particles’ single scattering polarization of about 1 liter, a mass smaller than 5 kg, and a power phase function, which depends stronger on the particle budget smaller than 5 W. size than the flux phase function. The spectral features

1 2 with a slight off-set to compensate for Mars’ rotation except possibly across Fraunhofer lines in F0 [see e.g. 10, 11] Mars Atmosphere: Modeling and Observations (2008) 9078.pdf

3

of the surface albedo also show up in P , especially at the References longer wavelengths, where the albedo is highest: because [1] J. E. Hansen and L. D. Travis. Light scattering in plane- the model surface completely depolarizes the light it re- tary atmospheres. Space Sci. Rev., 16:527–610, 1974. flects, P of the observed light decreases with increasing [2] M. I. Mishchenko and L. D. Travis. Satellite retrieval albedo. This effect is less strong when τD is large. How- of aerosol properties over the ocean using polarization as ever, a large τD will usually also lead to an increase of well as intensity of reflected sunlight. J. Geophys. Res., the multiple scattering (depending on the albedo of the 102:16989–17014, 1997. doi: 10.1029/96JD02425. dust particles) and hence to a decrease of P . [3] Y. Shkuratov, N. Opanasenko, E. Zubko, Y. Grynko, Our sample simulations illustrate the power of po- V. Korokhin, C. Pieters, G. Videen, U. Mall, and A. Opanasenko. Multispectral polarimetry as a tool to larimetry, and show that measuring both F and P across investigate texture and chemistry of lunar regolith parti- a broad spectral region gives better information about a cles. Icarus, 187:406–416, 2007. doi: 10.1016/j.icarus. planet than flux measurements alone, as was also con- 2006.10.012. vincingly shown by e.g. Mishchenko and Travis [2]. The [4] J. E. Hansen and J. W. Hovenier. Interpretation of the Po- simulations also show that the degree of polarization of larization of Venus. J. Atmos. Sci., 31:1137–1160, 1974. Mars can be several percents, which could induce errors [5] Y. Shkuratov, M. Kreslavsky, V. Kaydash, G. Videen, in polarization sensitive spectrometers [14]. J. Bell, M. Wolff, M. Hubbard, K. Noll, and A. Lubenow. Our simulations do not yet include realistically re- Hubble Space Telescope imaging polarimetry of Mars flecting Martian surfaces. Such surfaces are expected during the 2003 opposition. Icarus, 176:1–11, 2005. doi: 10.1016/j.icarus.2005.01.009. to have non-isotropic reflection, e.g. with an opposition [6] R. Santer, M. Deschamps, L. V. Ksanfomaliti, and surge, and to polarize the light. The reflection charac- A. Dollfus. Photopolarimetric analysis of the Martian at- teristics of a surface give information about the micro- mosphere by the Soviet MARS-5 orbiter. I - White clouds physical properties of the dust particles covering it [e.g. and dust veils. Astron. Astrophys., 150:217–228, 1985. 3]. In particular when the atmospheric optical thickness [7] R. Santer, M. Deschamps, L. V. Ksanfomaliti, and is small, SPEX can be used to study the surface. A. Dollfus. Photopolarimetry of Martian aerosols. II - Limb and terminator measurements. Astron. Astrophys., Summary and discussion 158:247–258, 1986. SPEX, the Spectropolarimeter for Planetary EXplo- [8] W. B. Sparks, J. H. Hough, and L. E. Bergeron. A Search for Chiral Signatures on Mars. , 5:737–748, ration, is an innovative, robust, and compact instrument 2005. doi: 10.1089/ast.2005.5.737. to measure the flux and state of polarization of sunlight [9] F. Snik, Th. Karalidi, Ch. Keller, E. Laan, R. ter Horst, reflected by a planet from about 400 to 800 nm. SPEX’ R. Navarro, D.M. Stam, Ch. Aas, J. de Vries, G. Oomen, current design (with 9 fixed viewing directions) is tai- and R. Hoogeveen. SPEX: An in-orbit spectropolarimeter lored for remote-sensing of Mars from an orbiter. With for planetary exploration. In SPIE Astronomical Instru- SPEX, the microphysical properties (size, shape, and mentation, 23-28 June, Marseille, France, 2008. composition) of dust and ice cloud particles in the at- [10] I. Aben, D. M. Stam, and F. Helderman. The ring effect mosphere and on the surface of Mars can be derived. in skylight polarization. Geophys. Res. Lett., 28:519–522, The project to build a breadboard model for SPEX 2001. doi: 10.1029/2000GL011901. has started in September 2008. Because of its small size [11] D. M. Stam, I. Aben, and F. Helderman. Skylight polar- ization spectra: Numerical simulation of the Ring effect. and its small power budget, SPEX could be added to al- J. Geophys. Res., 107:AAC 1–1 – 15, 2002. most any orbiter. An adapted version of SPEX (with less [12] J. F. de Haan, P. B. Bosma, and J. W. Hovenier. The and/or different viewing directions) could also be placed adding method for multiple scattering calculations of po- on the surface of Mars to study the sunlight that is dif- larized light. Astron. Astrophys., 183:371–391, 1987. fusely transmitted through the atmosphere. Indeed, a [13] D. M. Stam, J. F. De Haan, J. W. Hovenier, and SPEX has been proposed as payload for the balloon fore- P. Stammes. Degree of linear polarization of light emerg- seen in the Tandem mission to Titan to characterize the ing from the cloudless atmosphere in the oxygen A band. haze particles on this moon, and two SPEXes have been J. Geophys. Res., 104(13):16843–16858, 1999. doi: chosen as payload for a formation flying microsats mis- 10.1029/1999JD900159. sion to study aerosol and clouds on Earth. [14] D. M. Stam, J. F. De Haan, J. W. Hovenier, and I. Aben. Detecting radiances in the O2 A band using polarization- Acknowledgements The design and breadboarding sensitive satellite instruments with application to the study for SPEX is supported by a PEP-grant from NIVR Global Ozone Monitoring Experiment. J. Geophys. Res., (the Netherlands Agency for Aerospace Programmes). 105(14):22379–22392, 2000. Mars Atmosphere: Modeling and Observations (2008) 9078.pdf

4

Figure 2: Model calculations of F (left) and P (right) of sunlight reflected by Mars with different dust optical thick- nesses τD (at 630 nm) and dust particles with different effective radii reff . The surface reflects Lambertian with a wavelength dependent albedo. From top to bottom: τD= 0.2, reff = 1.0 µm; τD= 0.2, reff = 0.5 µm; τD= 1.0, ◦ reff = 1.0 µm; τD= 1.0, reff = 0.5 µm. Solar zenith angle θ0 is 60 , different lines represent 7 different SPEX view- ing angles θ (measured with respect to the downward vertical), plus-signs indicate azimuthal angles of 0◦, minus-signs azimuthal angles of 180◦ (i.e. the satellite has the sun in its back).