Articles

VIEWS FROM EPOXI: COLORS IN OUR SOLAR SYSTEM AS AN ANALOG FOR EXTRASOLAR PLANETS

, , , , , , , , , and

Published 2011 February 18 © 2011. The American Astronomical Society. All rights reserved.
, , Citation Carolyn A. Crow et al 2011 ApJ 729 130 DOI 10.1088/0004-637X/729/2/130

0004-637X/729/2/130

ABSTRACT

The first visible-light studies of Earth-sized extrasolar planets will employ photometry or low-resolution spectroscopy. This work uses EPOCh medium-band filter photometry between 350 and 950 nm obtained with the Deep Impact (DI) High Resolution Instrument (HRI) of Earth, the Moon, and Mars in addition to previous full-disk observations of the other six solar system planets and Titan to analyze the limitations of using photometric colors to characterize extrasolar planets. We determined that the HRI 350, 550, and 850 nm filters are optimal for distinguishing Earth from the other planets and separating planets to first order based on their atmospheric and surface properties. Detailed conclusions that can be drawn about exoplanet atmospheres simply from a color–color plot are limited due to potentially competing physical processes in the atmosphere. The presence of a Rayleigh scattering atmosphere can be detected by an increase in the 350–550 nm brightness ratio, but the absence of Rayleigh scattering cannot be confirmed due to the existence of atmospheric and surface absorbing species in the UV. Methane and ammonia are the only species responsible for strong absorption in the 850 nm filter in our solar system. The combination of physical processes present on extrasolar planets may differ from those we see locally. Nevertheless, a generation of telescopes capable of collecting such photometric observations can serve a critical role in first-order characterization and constraining the population of Earth-like extrasolar planets.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

With extrasolar planet detection becoming more commonplace, the frontiers of extrasolar planet science have moved beyond detection to the observations required to determine planetary properties. Once the existing observational challenges have been overcome, the first visible-light studies of extrasolar Earth-sized planets will likely employ filter photometry or low-resolution spectroscopy to observe disk-integrated radiation from the unresolved planet. While spectroscopy of these targets is highly desirable (e.g., Des Marais et al. 2002) and provides the most robust form of characterization, signal-to-noise ratio (S/N) considerations presently limit spectroscopic measurements of extrasolar worlds. Broadband filter photometry will thus serve as a first line of characterization.

The most ambitious mission concepts for visible-light spectroscopic and photometric characterization of extrasolar planets include large space-based telescopes that will employ coronographic imagers or occulters (e.g., Levine et al. 2009; Cash 2006) to greatly reduce the brightness of the parent star and reveal the fainter planet (Trauger & Traub 2007). The James Webb Space Telescope, currently under construction and planned to launch in 2014, may provide a more immediate coronographic capability for extrasolar planet studies if it can be paired with a free-flying occulter spacecraft (Soummer et al. 2009).

Visual inspection of solar system planets reveals an array of planetary colors, resulting from the diversity of atmospheric and surface composition and structure. While this handful of planets will not fully represent the diversity to be found even in the 430+ planets now known, these are the only ones for which we can currently obtain multi-wavelength high S/N visible photometry.

Previous attempts to obtain disk-integrated spectra or photometry of Earth and Mars have been made by interplanetary spacecraft en route to other targets in the solar system. The first spacecraft mission dedicated to whole-disk Earth observations as a principal scientific goal is the Extrasolar Observation and Characterization (EPOCh) investigation, which is part of NASA's EPOXI mission (Ballard et al. 2010; Cowan et al. 2009). EPOCh utilized the Deep Impact (DI) flyby spacecraft to provide time-resolved multi-wavelength medium-band filter photometry and moderate-resolution spectroscopy of Earth as an analog experiment for photometric observations of an extrasolar Earth-like planet. These data were taken from distances of 0.18–0.3 AU, providing multi-wavelength coverage of five full rotations of Earth from equatorial and polar views, and over several seasons (Livengood et al. 2009). EPOCh observed a transit of Earth by the Moon in 2008 May, facilitating a comparison between two terrestrial bodies with significantly differing surface properties (Deming et al. 2007). Near the end of the mission, medium-band filter photometry of Mars was also obtained.

In this paper, we use EPOCh filter photometry of Earth, the Moon and Mars, model spectra, and previous photometric and spectroscopic observations of a range of the solar system planets, Titan, and the Moon to explore the limitations of using color as a baseline for characterizing extrasolar planets. We expand upon the work of Traub (2003) to determine the optimal filters for distinguishing Earth-like, terrestrial-like, and Jovian-like bodies.

2. EARTH, MOON, AND MARS OBSERVATIONS

2.1. Previous Observations

One of the earliest studies of Earth for the purpose of understanding whether its habitability could be remotely determined was undertaken by the Galileo spacecraft as it flew by Earth during a gravity assist maneuver on the way to its primary mission at Jupiter (Sagan et al. 1993). The data collected on this flyby included spatially resolved UV and NIR spectra and visible photometry, as well as observations at radio wavelengths. Whole-disk high-resolution spectroscopy at visible wavelengths has since been acquired by the Visible and Infrared Thermal Imaging Spectrometer (VIRTIS) aboard Venus Express, but has yet to be analyzed (Coradini et al. 1998; Grinspoon et al. 2008).

Photometric studies of the Moon have until recently been restricted to ground-based observations and thus limited to only the lunar near side. McCord et al. (1972), using low-resolution filter photometry, studied the spectral characteristics of various near side lunar regions including mare, upland, and bright craters using ground-based telescopes calibrated relative to laboratory-reflected light measurements of Apollo 11 returned samples. Observations of the lunar surface were collected with the 21 inch and 60 inch telescopes at the Mount Wilson Observatory using 24 medium-band filters, between 0.3 and 1.1 μm. The McCord data are characterized by variations in the lunar surface reflectance between 800 and 1000 nm. See Section 4 for more details.

Three spacecraft missions observed the lunar far side through narrow- to medium-band filters similar to those of the High Resolution Instrument (HRI) used in the data discussed herein. The data from Galileo's solid state imaging cameras (Belton et al. 1992), Clementine's UV/Visible instrument (Hillier et al. 1999), and SMART-1 orbiter's visible and NIR camera AMIE (Advanced Moon micro-Imager Experiment; Josset 2008) were all at higher spatial resolution than the EPOCh images though no single data set spanned the UV to NIR spectral regions as completely as the HRI.

Multiple full-disk brightness measurements of Mars have been collected with ground-based telescopes (e.g., Younkin 1966; Irvine et al. 1968a, 1968b; McCord & Westphal 1971), but they are limited by atmospheric interference and do not cover a complete rotation of Mars. Visible spectroscopic observations of Mars have been collected by Mars Global Surveyor and are not contaminated by terrestrial atmospheric effects (e.g., Rogers et al. 2007), but they are at higher spatial resolution than the EPOCh Mars data and are not full-disk observations. The EPOCh data set therefore afforded the unique opportunity to observe the full disk of Mars over a full revolution without the effects of neither Earth's nor Mars' atmosphere.

2.2. EPOCh Observations

The EPOCh investigation planned five observations of Earth for 2008. EarthObs4, presented here, was conducted over 24 hr beginning approximately on 20:00 UT 2008 May 28 and included a transit of Earth by the Moon; see Figure 1(a). This sequence provided contemporaneous photometry and spectroscopy of the Moon and Earth. It also allowed for observations of the far side of the Moon, which is not visible from Earth due to the Moon's synchronous rotation. The illuminated region of the Moon in the 2008 May images spanned from 160° W to 260° W longitudes (approximately half of the lunar far side) and the phase angle for both Earth and Moon averaged 75fdg1. An additional set of calibration data from 2005 January 16 was used to analyze the near side lunar reflectance for comparison with the far side data. Approximately 2/3 of the illuminated surface of the front side of the Moon was visible in the calibration images and encompassed the southeast near side Balmer crater regions; see Figure 1(b) (Maxwell & Andre 1981). The phase angle for these images was 97fdg6. In addition to lunar and Earth data, we analyzed 64 sets of full-disk observations of Mars collected over a 24 hr period starting on 2009 November 20. The observations covered a nearly complete rotation period of Mars at a phase angle of 37° (Figure 1(c)). Figure 2 shows the individual data sets for each body and Table 1 contains details on observation parameters for all data sets. The variations in the wavelength-dependent brightness of the Moon from the different observations shown in Figure 2 are relatively small. For Mars, the variations are larger and are due to the planet's rotation causing different surfaces with distinct compositions rotating into and out of view. Similarly, for Earth the variations are due to the changing scene including cloud coverage and structures, continents, and oceans on the illuminated disk as Earth rotates. The variability demonstrated in Figure 2 does not affect the relative reflectance values used in this work to characterize the worlds.

Figure 1.

Figure 1. (a) Color composite image of lunar transit of Earth using the HRI 450, 550, and 650 nm filters. Half of the lunar far side is illuminated, and the Sahara Desert of northern Africa is identified by the orange region on the illuminated Earth. Credit: NASA/JPL-Caltech/UMD/GSFC. (b) Calibration image from 2005 January 16 used to calculate the near side reflectance of the Moon. The region visible encompasses the Balmer crater region located in the southeast of the lunar near side. (c) HRI image of Mars acquired on 2009 November 20.

Standard image High-resolution image
Figure 2.

Figure 2. Unnormalized, summed, I/F (intensity divided by solar flux) for each image of the Earth, Moon, and Mars. These are displayed to show the variance in I/F during the observed time interval, the Earth and Mars being observed over one full rotation. The value, I/F, is a function of the reflectance characteristics of the body at the observed heliocentric distance, spacecraft to target distance, and phase angle. Since these are different for all three bodies, their unnormalized data cannot be compared.

Standard image High-resolution image

Table 1. EPOCh Data Set Details

Data Set Date Time Phase Angle Geocentric Distance Heliocentric Distance
    (UT) (deg) (AU) (AU)
Lunar calibration 2005 Jan 16 20:10:20 97.6 0.0110 0.9826
EarthObs4: Initial image 2008 May 29 2:03:45 75.0 0.3308 0.9811
EarthObs4: Final image 2008 May 29 13:03:45 75.1 0.3313 0.9818
Mars Obs: Initial image 2009 Nov 20 11:53:18 37.0 0.0717 0.9814
Mars Obs: Final image 2009 Nov 21 11:54:38 37.0 0.0701 0.9811

Note. Columns 5 and 6 are spacecraft geocentric and heliocentric distances, respectively.

Download table as:  ASCIITypeset image

The EPOCh observations were taken using the HRI on the DI flyby spacecraft (Hampton et al. 2005), a 0.3 m f/35 optical Cassegrain telescope equipped with a 9-position filter wheel and a CCD camera. Two of the positions have broadband white-light filters centered at 650 nm and the other seven are medium-band filters with 100 nm bandwidths evenly spaced between 300 and 1000 nm (Hampton et al. 2005). Images of Earth and the Moon were taken through the 350, 750, and 950 nm filters hourly and through the remaining medium-band filters every 15 minutes, with exposure times ranging from 8 to 73 ms depending on the filter throughputs. See Table 2 for a list of filter center wavelengths and exposure times. We used 8 of the 12 sets of lunar transit observations for our analysis. The remaining four sets were discarded because either the Moon was transiting Earth or the gradient of scattered Earthlight across the Moon was too large for a single background value subtraction.

Table 2. HRI Filter Details for Earth Observations

Filter Number Wavelength Exposure Time
  (nm) (ms)
2 450 13
3 550 8
4 350 73
5 950 61
7 750 13
8 850 26
9 650 9

Download table as:  ASCIITypeset image

3. EPOXI DATA PROCESSING

The raw data from the DI spacecraft are processed through a calibration pipeline that converts the raw data numbers into measured flux in units of W m−2 sr−1 for a given bandpass (Klaasen et al. 2008). In the pipeline, the images are first decompressed, saturated pixels are flagged, and a dark frame is subtracted. Crosstalk functions are then applied to eliminate ghosts of the other three quadrants of the CCD, and the images are divided by the flat field to mitigate pixel-to-pixel variations. The HRI uses a frame-transfer CCD, so a desmearing process is necessary to correct for residual persistent smear left by bright objects during the frame transfer from the illuminated portion of the chip to the shielded part of the chip. The images are next interpolated over bad pixels and gaps, and then median filtered to produce a sharper image (Klaasen et al. 2008). We used I/F images that are produced by dividing the pipeline-calibrated images by the solar spectrum. The resulting images are therefore a measure of the reflectance of the planetary bodies.

Additional processing was needed in order to analyze the EPOXI lunar data. We subtracted the background of scattered Earthlight from the lunar data to get a more accurate measurement of lunar reflectance. The background level was determined by taking multiple profiles of the Moon in each pipeline-calibrated image. For most of the data, the background level was relatively constant across the disk of the Moon, so a single background value was used for subtraction. The background level decreased with increasing distance from Earth and averaged around 2%. The highest background levels were in the 350 nm filter when the Moon was close to Earth, reaching a maximum of 13% of the total signal.

Following background determination, we calculated the integrated brightness of the lunar disk by defining a region around the Moon and summing the background-subtracted pixel values within the region. This process was also used to calculate the disk-integrated reflectance of Earth and Mars, but no background was subtracted for these data. To calculate the lunar near side colors, we used calibration data from 2005 January. Only the Balmer crater region on the southeast of the lunar near side, not the entire illuminated surface, was visible in the image. For analysis we defined a region to encompass the portion that was visible, and no background was subtracted from these data either.

4. RESULTS AND ANALYSIS

4.1. The Moon and Mars

The reflectance of the lunar far and southeast region of the near side, shown in Figure 3, increases with a relatively linear slope from the UV to the NIR except in the 950 nm filter of the far side. In this bandpass, the disk-integrated far side exhibits an absorption feature not present in the summed images of the near side region observed by the DI spacecraft. The variation in the 950 nm filter is consistent with narrowband photometry of lunar regions from McCord et al. (1972) presented in Figure 4. Absorptions are seen in the McCord data for the lunar standard region Mare Serenitatis 2 (d) and are also prominent in other bright crater regions (e, f, and g). The disk-integrated reflectance of Mars shown in Figure 2 is also characterized by a nonlinear increase from the UV to the NIR with an absorption feature at 950 nm.

Figure 3.

Figure 3. Normalized UV to NIR reflectance (W m−2 sr−1) of far side and southeast region of the lunar near side obtained with the HRI and Sea of Serenity 2 (also known as Mare Serenitatis 2) from McCord et al. (1972) scaled to 550 nm. The near and far side HRI reflectance vary in the 950 nm filter. The feature at 950 nm in both the HRI and Mare Serenitatis 2 data is due to Fe2+ absorption. The linear slope between UV and NIR in all three spectra is due to charge transfer absorptions within mafic silicates on the lunar surface.

Standard image High-resolution image
Figure 4.

Figure 4. Figure 1 from McCord et al. (1972) (Copyright 1972 American Geophysical Union. Reproduced by permission of American Geophysical Union.) showing ground-based reflectance measurements of various lunar regions. Absorption at 950 nm is also in McCord spectra of bright crater regions (e), (f), and (g) and minimal in the spectra of upland and mare regions (a) and (b). These variations are consistent with the EPOCh lunar far and southeast near side region data.

Standard image High-resolution image

The Moon and Mercury are airless bodies and Mars' atmosphere is transparent in the visible and NIR wavelengths. Consequently, their reflected light is dominated by the surface composition and the texture of their regolith (Soderblom 1992; Roush et al. 1993; McKay 1996). Each body experienced differentiation soon after it was formed, followed by an extended period of regolith formation primarily from repeated impacts during the late heavy bombardment period of solar system evolution (3.8–4.1 Gyr ago; Tera et al. 1974). Photons hitting the regolith surface are both scattered from grain boundaries and absorbed and reradiated after interacting with the minerals in the soil (Burns 1993; Hapke 2005).

In the case of the Moon, the regolith contains minerals, glass and nano-phase iron (Pieters et al. 2000). The steep slope toward the NIR in the lunar data is primarily due to the continuously varying absorption coefficient of nano-phase iron that is produced by reduction of Fe and vapor deposition resulting from micrometeorite bombardment of the lunar surface. A large portion of the observed region of the lunar far side encompasses the South Pole Aitken Basin, which is enriched in mafic minerals relative to the near side (Nakamura et al. 2009). Charge transfer interactions with Fe2+ on the surface are responsible for the 950 nm absorption in the lunar spectrum (Lucey et al. 1998). It is important to note here that the 950 nm absorption cannot be due to the hydrated materials on the surface of the Moon detected by Sunshine et al. (2009) using the DI Medium Resolution Instrument (MRI) and the high-resolution near-IR spectrometer. The amount of lunar water is too small for its absorption effects to be measurable in the low spectral resolution HRI filter data.

In the case of Mars, the surface likely consists of fine-grained loess-like material as discussed in Singer et al. (1982). Surface mineral components containing optically active absorption bands in the visible are a combination of charge transfer absorption and spin-forbidden crystal field transitions due to Fe3+ and Fe2+ cations in ilmenite and pyroxene minerals. If the linear nature of the UV absorption in the lunar spectrum is indicative of the presence of nano-phase iron, it would appear that the disk-integrated Mars spectrum is not influenced by nano-phase iron. It is important to note here that the large vertical bars placed on the Mars data at longer wavelengths in Figure 5 do not reflect errors within the measurement, but instead reflect variations in the reflectance of the Martian surface observed during one full rotation. Irvine et al. (1968a, 1968b) also collected photometric data of Mars at various phases and rotational angles over a two-year time period (see Section 4.3). His spectrum of Mars at a 36fdg7 phase angle falls within the range of values observed with the HRI and is therefore consistent with our data.

Figure 5.

Figure 5. Disk-integrated photometry of Earth, Moon, and Mars from EPOCh HRI observations. The Moon's reflectance increases from UV to NIR due to photon interactions with mafic silicates and nano-phase iron on the lunar surface; the absorption at 950 nm is a signature of Fe2+. Mars exhibits a steeper UV to NIR slope than the Moon due to the presence of both iron oxides and mafic silicates on its surface. The iron oxides are stronger absorbers of UV wavelengths and cause Mars' red reflectance to be relatively higher. Earth shows a turn-up toward the UV due to Rayleigh scattering and absorption at 950 nm due to water. Vertical bars for EPOCh data reflect the range of observed full-disk reflectance values and not the instrument error.

Standard image High-resolution image

4.2. Earth

The disk-integrated reflectance spectrum of Earth shown in Figure 5 exhibits a steep incline toward UV wavelengths and is relatively flat at green to NIR wavelengths. The flat shape of the spectrum at longer wavelengths is attributed to clouds. Although, a slight rise in Earth's reflectivity from the 650 nm filter to the 850 nm filter is due to the fact that many surfaces on Earth (e.g., soils, sands, and vegetation) are more reflective at NIR wavelengths than at visible wavelengths. An absorption feature can be seen in the 950 nm filter due to a water overtone band located within this wavelength range (Traub 2003).

The increase in Earth's reflectance toward the UV is the result of Rayleigh scattering in the atmosphere. To investigate the significance of Rayleigh scattering in our Earth data, we produced two simulated high-resolution Earth spectra using the NASA Astrobiology Institute's Virtual Planetary Laboratory (VPL) three-dimensional spectral Earth model (Robinson et al. 2010) and convolved them with the HRI filter transmission curves, see Figure 6. In the first spectrum (solid line), realistic Rayleigh scattering is assumed and a rise in reflectivity at blue wavelengths is reproduced. The model used to produce the second spectrum (dashed line) is identical to the first except that Rayleigh scattering has been removed from the radiative transfer calculations. A strong rise in reflectivity toward blue wavelengths is no longer seen, although a slight increase is attributed to the blue color of Earth's oceans. These models agree with our assumption that the high near-UV reflectance of Earth is primarily due to Rayleigh scattering.

Figure 6.

Figure 6. Model spectra produced using NASA Astrobiology Institute's Virtual Planetary Laboratory three-dimensional spectral Earth model convolved with HRI filter transmission curves. The model with Rayleigh scattering (solid line) has significantly higher UV values than the model with scattering removed (dashed line); a slight increase in UV values of the latter is due to the blue color of Earth's oceans. These models verify that Rayleigh scattering is the dominant process responsible for the high UV values in the EPOCh HRI observations.

Standard image High-resolution image

Traub (2003) discusses the prospect of detecting vegetation on a planetary body by looking for the vegetation red edge (VRE), which he was unable to detect in his analysis of Earthshine data. The VRE is defined as

Equation (1)

where rIR and rR are the reflectivities at NIR and red wavelengths, respectively. Narrowband filters designed to bound the steep rise in the reflectivity of planets around 700 nm are usually employed to measure the VRE. For example, a pair of 10 nm wide filters centered at 675 nm and 775 nm give a VRE value of ∼100% for a typical reflectance spectrum of a planet, ∼5% for a high-resolution reflectance spectrum of Mars, and between 0% and 10% for Earth (Arnold 2008), where the vegetation signature is washed out by other features on the planet (e.g., clouds). Cowan et al. (2009) used the EPOXI data set to correlate reddening in Earth's disk-integrated reflectance spectrum with the presence of landmass on the visible, illuminated disk. Vegetation-covered land was not distinguished from other land types because the medium-band HRI filters are nonideal for characterizing a steep rise in reflectivity and cannot properly distinguish between a gradual rise in reflectivity toward red wavelengths (as is the case for Mars, Mercury, and the Moon) and a steep rise in reflectivity near 700 nm, which is seen in some high-resolution, disk-integrated Earth spectra (Arnold et al. 2002).

4.3. Other Solar System Worlds

To determine how the EPOCh data compare with the reflectance spectra of other solar system bodies, we expanded our analysis to include other terrestrial and Jovian worlds. Observations for these bodies were not collected with the DI HRI, so full-disk reflectance data from previously published studies were used for our analysis. We reproduced reflectance spectra of Mercury and Venus from Irvine et al. (1968a, 1968b) collected at Le Houga Observatory, France and Boyden Observatory, South Africa. The observations were taken through 10 medium-band filters centered between 315 nm and 1060 nm, and the calculated magnitudes of the planets were normalized to a distance of 1 AU from the Sun. We chose the observations with phase angles closest to that of the May EPOCh observations of 75°; the angles for the Mercury and Venus were 77° and 76°, respectively. The band centers of the filters used by Irvine et al. are similar to the HRI filters, so a linear interpolation and conversion to reflectance was necessary for comparison with the EPOCh data. Figure 7 shows a comparison between the EPOCh and Irvine data sets.

Figure 7.

Figure 7. EPOCh data for Earth, Moon, and Mars overlaid with Irvine et al. (1968a, 1968b) data for Venus and Mercury. Irvine data, which were originally reported in magnitudes with an error of ±0.02 mag, were converted to normalized reflectance for comparison with EPOCh data. The UV downturn in the Venusian atmosphere is due to atmospheric aerosol absorbers paired with a high cloud deck prohibiting appreciable Rayleigh scattering. The Venusian cloud deck is relatively non-absorptive at longer wavelengths causing the reflectance to be flat in the NIR. Mercury's reflectance increases from UV to NIR due to charge transfer of low-iron silicates. This figure demonstrates that Earth is the only body with an increase in reflectivity toward UV wavelengths due to Rayleigh scattering.

Standard image High-resolution image

The terrestrial worlds, other than Earth, are all less reflective in the UV than in the NIR. Mercury, an airless body, has a reflectance spectrum increasing steadily from the UV to the IR due to light interactions with surface minerals (Vilas 1985). The Mercurian soil is low in Fe with <2–3 wt % at its surface, so Fe is not a major contributor to charge transfer slope in the reflectance of Mercury (McClintock et al. 2008). The visible colors of Mercury are not diagnostic of the surface mineral composition.

The reflectivity of Venus shows a depression in the UV and a relatively flat reflectance from the green to the NIR filters. Venus has a thick atmosphere unlike the airless bodies discussed previously, so its atmospheric structure and composition are solely responsible for the planet's colors. Venus' opaque atmosphere contains several distinct layers of sulfuric acid clouds, which are distinguished by their aerosol particle size. The cloud deck which is primarily responsible for Venus' high reflectivity is located between about 50 and 80 km in altitude (Marov et al. 1973; Marov 1978). Extinction optical depths through the Venusian haze (which extends up to about 90 km in altitude) and the upper cloud deck reach unity by about 30 mbar (Knollenberg et al. 1980), so no appreciable amount of Rayleigh scattering can occur. The effects of a "UV absorber" can be seen in the reflectance spectrum of Venus blueward of about 650 nm. The absorption efficiency of this substance increases from 650 to 300 nm, explaining the drop in Venus' reflectivity at blue wavelengths (Pollack et al. 1980). The flat reflectance of Venus' reflectivity curve beyond about 650 nm is explained by the sulfuric acid droplets' essentially non-absorptive behavior at visible wavelengths (Crisp 1986).

In addition to Irvine's photometric data, we used high spectral resolution, full-disk reflectance data for Jupiter, Saturn, Uranus, Neptune, and Titan from Karkoschka (1994). The observations were collected at the European Southern Observatory and span the wavelength range of 300–1000 nm at phase angles of 9fdg8, 2fdg7, 0fdg6, 0fdg4, and 2fdg7, respectively. We convolved these spectra with the HRI filter transmission curves to determine the medium-band reflectance spectra of the five bodies. An overlay of the high-resolution reflectance spectrum of Jupiter from Karkoschka (1994) and the HRI filter transmission curves are seen in Figure 8, and the computed broadband photometric reflectance values for the Jovian worlds are shown in Figure 9.

Figure 8.

Figure 8. Karkoschka (1994) high-resolution spectrum of Jupiter overlaid with HRI–VIS filter transmission curves (Hampton et al. 2005).

Standard image High-resolution image
Figure 9.

Figure 9. Karkoschka (1994) high-resolution spectra and equivalent HRI reflectance values computed by convolution with HRI–VIS filter transmission functions. The reflectance of Jupiter, Saturn, and Titan has strong absorption due to aerosols in the UV. Uranus and Neptune have competing aerosol absorption and Rayleigh scattering causing their reflectivity to be flat UV spectra. All five bodies exhibit varying amounts of methane absorption in the NIR.

Standard image High-resolution image

Uranus and Neptune are unique in that their reflectance spectra are relatively flat at blue and UV wavelengths while containing strong absorption features in the red and NIR. A perfectly scattering, semi-infinite atmosphere consisting of pure H2 would have a wavelength-independent, zero-phase reflectance of about 0.8 (Prather 1974), which is about 30%–40% larger than the reflectance values reported by Karkoschka (1994, 1998) for Uranus and Neptune at short wavelengths. The discrepancy between the theoretical value and the measured value is due to Raman scattering, which decreases the zero-phase reflectance values for these worlds by about 5%, and the presence of aerosols (Karkoschka 1994). In the case of Uranus, a methane aerosol haze near 1.3 bar and a thick cloud of H2S aerosols significantly reduce the reflectance values at UV and near-UV wavelengths (Baines & Bergstralh 1986). The methane haze preferentially absorbs blue wavelengths, which is possibly due to contamination from a thin stratospheric haze (Rages et al. 1991), leading to a drop in Uranus' reflectance between the 550 nm and 350 nm filters. Similarly for Neptune, a methane haze near 1 bar in the atmosphere and photochemically produced stratospheric hazes of hydrocarbons reduce the reflectance values at blue and UV wavelengths (Baines & Smith 1990). Both Uranus and Neptune experience strong CH4 absorption features at wavelengths longer than about 500 nm (Karkoschka 1998). The wings of these absorption features deplete most of the continuum at red and NIR wavelengths, and are especially strong in the 850 nm and 950 nm HRI filters. The single-scattering albedos of aerosols in the atmospheres of Uranus and Neptune decrease between green and red wavelengths, causing a slight drop in the continuum reflectivity of these worlds at red wavelengths (Baines & Bergstralh 1986; Baines & Smith 1990).

The reflectance spectra of Jupiter, Saturn, and Titan are different from the ice giants as they exhibit weaker absorption in the NIR and stronger aerosol absorption at blue and UV wavelengths. High-resolution reflectance spectra of Jupiter and Saturn are similar in appearance. Both are dominated by methane absorption at wavelengths longer than about 700 nm, especially in the 850 nm and 950 nm HRI filters (Karkoschka 1994, 1998). At shorter wavelengths, the reflectance spectra of these worlds begin to decrease due to aerosol absorption. The composition of these aerosols is currently unknown, but is commonly referred to as "chromophores." Compounds that may explain the blue absorption in the atmospheres of Jupiter and Saturn include photochemically produced hydrocarbons and/or sulfur-bearing compounds (see West 2007; Simon-Miller et al. 2001). Models predict multiple cloud decks in the atmospheres of Jupiter and Saturn (Weidenschilling & Lewis 1973), which was verified in the case of Jupiter by the Galileo probe (Ragent et al. 1998). The uppermost cloud deck of these worlds is composed of ammonia ice, which tends to form at higher pressures (lower altitudes) in the Saturnian atmosphere due to Saturn's colder temperature profile (Atreya et al. 1999; Gierasch & Conrath 1993). A slight upturn in reflectance at UV wavelengths for Saturn is attributed to Rayleigh scattering (McCord et al. 1971). The spectral behavior of Saturn in this wavelength region is controlled by competition between aerosol absorption and Rayleigh scattering.

Titan's reflectance spectrum is similar to that of Jupiter and Saturn, and exhibits absorption at both UV and NIR regions. At visible wavelengths, Titan's reflectance spectrum is dominated by absorption from a stratospheric haze at short wavelengths, scattering and absorption by nitrogen-rich hydrocarbon aerosols (tholins) at intermediate wavelengths, and gas absorption at longer wavelengths. Spectroscopic and laboratory studies have revealed that the tholins which make up Titan's haze generally have an increasing absorption efficiency with decreasing wavelength through the visible (Rages & Pollack 1980; McKay 1996; Tomasko et al. 2008). As a result, Titan's reflectance spectrum decreases steadily from 650 nm to UV wavelengths. Methane absorption dominates Titan's reflectance spectrum longward of 650 nm. Although these features are not as strong as those in the spectra of Uranus and Neptune, the 850 nm and 950 nm HRI filters show marked decreases in Titan's reflectance due to methane absorption.

5. DISCUSSION

The objective of our study was to analyze the colors of the planets within our solar system and use them to create a baseline for characterizing extrasolar planets. Traub (2003) discusses the benefits of using color to broadly characterize the types of planets detected by Terrestrial Planet Finder and other similar missions. Similar to his analysis, we used data from Irvine et al. (1968a, 1968b) and Karkoschka (1994) to reproduce the colors of terrestrial and Jovian worlds. We additionally presented photometric observations of Earth, Moon, and Mars taken with the DI spacecraft. Our data improve upon Traub's study, in which Earth's colors were from Earthshine and modeled Earth spectra and Mars' colors were from ground-based data that are contaminated by the terrestrial atmosphere.

Traub proposed using a color–color diagram and defined three broadband filters: 400–600, 600–800, and 800–1000 nm. Although his choice of bandpasses separates planets into groups, he acknowledges that the filter selection could be improved upon. The DI HRI afforded us the opportunities to observe Earth, Moon, and Mars from space, explore a range of filter combinations, and to determine the optimal filters for distinguishing between different types of planets. We chose the three filter combinations that reveal Earth's unique characteristics and partitioned the solar system bodies into color groups. Figure 10 shows the resulting color–color diagram of the reflectance of eight planets plus Titan and the Moon through the HRI 350, 550, and 850 nm filters. The ratio between the reflectivity in the 350 and 550 nm filters is plotted on the vertical axis, and the ratio between the reflectivity in the 850 and 550 nm filters is plotted on the horizontal axis. The lines where these ratios are equal to unity are also shown, and a perfectly reflective body would fall where they intersect.

Figure 10.

Figure 10. Color–color plot with ratio of reflectance in HRI 350 and 550 nm filters on the vertical axis and 850 and 550 nm filters on the horizontal axis. The lines of unity are plotted for reference and the Sun's radiance resides at the intersection of the two lines. Bodies with 350:550 > 1 have atmospheres dominated by Rayleigh scattering and those with 350:550 < 1 are dominated by atmospheric or charge-transfer absorption. Although 350:550 is diagnostic of the presence of Rayleigh scattering it is not sufficient to determine its absence. The 850:550 groups the bodies into three regions: (1) airless bodies on the right, (2) intermediate cloudy atmospheres near unity, and (3) strong NIR absorbing atmospheres on the left. Earth is the only body that resides in the upper right quadrant of the diagram due to Rayleigh scattering in its intermediate cloudy atmosphere. Bars for the EPOCh targets demonstrate the range of observed full-disk reflectance values where as the bars for the other bodies represent instrument error.

Standard image High-resolution image

The ratio between the reflectivity in the 350 and 550 nm filter in Figure 10 characterizes blue and UV reflectance. Earth, Uranus, and Neptune fall above or near the line of unity, meaning that they are blue or white. The increase toward UV wavelengths in the reflectance of these worlds is due to Rayleigh scattering in their atmospheres. Uranus and Neptune also have Raman scattering and H2S or hydrocarbon absorption in their atmospheres consequently reducing the effect of Rayleigh scattering. If H2S and other hydrocarbons were not present in the atmospheres of Uranus and Neptune, the planets would have 30%–40% higher reflectance at shorter wavelengths. Earth is therefore the bluest of the planets because it has no absorbing species that counter Rayleigh scattering at these wavelengths. This suggests that the ratio of reflectance at 350 and 550 nm could be used to detect the presence of a Rayleigh scattering atmosphere. However, the ratio is limited in determining the absence of Rayleigh scattering. As seen in the spectrum of Saturn, atmospheric absorption can overpower the effects of Rayleigh scattering and result in a decrease in UV reflectance. When analyzing the colors of extrasolar planet reflectance, higher resolution data will be needed to determine whether a UV color value near unity is due to absorption or an absence of Rayleigh scattering.

On the horizontal axis, the 850–550 nm reflectivity ratio is dependent on the NIR reflectance. Mercury, Venus, Earth, Moon, and Mars all have 850:550 > 1, with Mars being the reddest planet, and Earth and Venus being relatively white. The airless bodies have high relative reflectance in the NIR due to charge transfer interactions and crystal field interactions with minerals, mainly iron-bearing silicates and glass, on their surfaces. The two terrestrial worlds with atmospheres appear white due to clouds and an absence of atmospheric absorbing species at longer wavelengths. Conversely, Jupiter, Saturn, Uranus, Neptune, and Titan all have 850:550 < 1 as a result of strong methane absorption in their atmospheres. Ammonia also absorbs radiation at longer wavelengths, but is only attributed to a small fraction of the NIR absorption of these atmospheres. The most common physical process in our solar system responsible for relatively low NIR reflectivity is strong methane absorption, so similarly to Rayleigh scattering, our analysis suggests that the presence of a CH4 absorbing atmosphere is evident from planetary colors.

Traub (2003) produced a similar color–color plot of the planets in our solar system using their broadband colors. His results vary from Figure 10 because his filters were broader than those of the HRI and consequently could not detect slight variations in color due to differences in surface processes. For instance, Traub's filters average out the increase at short wavelengths due to Rayleigh scattering causing Uranus and Neptune to fall further to the blue end of the color–color plot than their positions in Figure 10. This effect may also be due to the fact that the HRI filters extend further into the UV where the effect of Rayleigh scattering is more prominent as a result of the λ−4 wavelength dependence. Another difference in Traub's plot is that Mercury appears redder than Mars. This could be due to either the higher quality of HRI data and/or to selection effects of the filters. In both cases, the effectiveness of grouping the bodies based on surface phenomena is increased by the narrower bandpasses of the HRI filters.

The general trend along the 850:550 nm axis places airless bodies toward the right with high NIR reflectivity, cloudy bodies near unity, and worlds with CH4 absorbing atmospheres toward the left. The pairing of this separation with the 350:550 ratio, which is dependent on Rayleigh scattering, shows that Earth is unique in our solar system. It is the only major solar system body with a cloudy, Rayleigh scattering atmosphere and consequently the only world in the upper right quadrant of the color–color diagram.

Our study finds that the HRI filters are successful in distinguishing Earth from other major solar system bodies and give evidence of possible processes responsible for planetary colors. As such, this filter combination has the potential for defining a more cost-efficient test and initial characterization of Earth-like exoplanets. The caveat is that currently known exoplanets are different from the worlds in our solar system. Thus, if future photometric observations of extrasolar planets were to place a world close to Earth on the color–color plot, it would indicate a necessary condition but is not sufficient to indicate a positive detection of an Earth-like planet. It could, however, suggest that some of the same physical processes responsible for Earth's colors are present on the extrasolar world. Future work will explore the effects of different combinations of physical processes on the colors of extrasolar planets. Model spectra of various potential planets may help determine if our filter choices are optimal for extrasolar planet characterization.

6. CONCLUSION

This work uses EPOCh HRI observations of Earth, the Moon, and Mars in addition to previous full-disk observations of the other six solar system planets and Titan to analyze the limitations of using photometric colors to characterize extrasolar planets expanding on work by Traub (2003). We determined that the 350, 550, and 850 nm HRI filters were optimal for characterizing planetary colors in our solar system by producing a color–color plot using relative full-disk reflectance through these filters. The ratio between the reflectivity in the 350 and 550 nm filters can be used to determine the presence of Rayleigh scattering with all the bodies containing appreciable Rayleigh scattering having 350:550 > 1. The ratio between the reflectivity in the 850 and 550 nm filters divides the worlds into three groups: those with strong methane absorbing atmospheres, with non-absorbing cloudy atmospheres, and with no atmospheres. The effect of using all three filters is to distinguish Earth as unique because it is the only body with a Rayleigh scattering non-absorbing cloudy atmosphere in our solar system as gleaned from disk-integrated observations.

Using these photometric colors, we find that Earth is easily distinguished from other solar system planets. As in this study, detailed conclusions that can be drawn about exoplanet atmospheres simply from a color–color plot are limited due to potentially competing physical processes in the atmosphere. For example, there are multiple mechanisms that can produce low UV reflectance and high IR reflectance. In the UV, both aerosol absorption and charge transfer with surface materials cause a decrease in relative reflectance. In the IR, both cloudy, non-absorbing atmospheres and crystal field absorptions on planetary surfaces produce 850:550 ratios greater than 1. We are able to distinguish between the different physical processes responsible for planetary colors only through analysis of higher resolution spectroscopy. Nevertheless, a generation of telescopes capable of collecting such photometric observations can serve a critical role in first-order characterization and constraining the population of Earth-like extrasolar planets. Such missions will function as a cost-efficient precursor for guiding the next generation of missions with advanced spectroscopic capabilities.

APPENDIX

EPOCh observations, convolved Karkoschka data, and Irvine relative reflectance values are listed in Tables A1A3, respectively.

Table A1. EPOCh Observations

Filter Center Earth Lunar Far Lunar Near Mars
(nm) (I/F) Side (I/F) Side (I/F) (I/F)
350 1610.47 13.83 1894.24 4.52
450 1342.71 18.06 2387.17 6.87
550 1059.02 22.54 2832.17 13.09
650 1014.27 26.68 3269.35 21.60
750 1062.25 30.74 4091.49 25.32
850 1168.00 34.86 4787.91 26.55
950 853.79 34.14 5506.36 24.31
Phase angle (deg) 75.1 75.1 97.6 37.0

Download table as:  ASCIITypeset image

Table A2. Convolved Karkoschka Data

Filter Center Jupiter Saturn Uranus Neptune Titan
(nm) (I/F) (I/F) (I/F) (I/F) (I/F)
350 0.604 0.445 0.984 1.252 0.344
450 0.824 0.684 1.067 1.235 0.584
550 1.000 1.000 1.000 1.000 1.000
650 1.006 1.144 0.647 0.555 1.248
750 0.832 1.001 0.292 0.237 1.101
850 0.636 0.776 0.148 0.132 0.879
Phase angle (deg) 9.8 2.7 0.6 0.4 2.7

Note. Values from convolution of Karkoschka (1994) data with HRI filter transmission functions.

Download table as:  ASCIITypeset image

Table A3. Irvine Relative Reflectance Values

Filter Center Mercury Venus
(nm)    
314.7 0.61 0.59
359.0 0.56 0.64
392.6 0.66 0.72
415.5 0.71 0.58
457.5 0.86 0.95
501.2 1.00 1.00
626.4 1.37 1.16
729.7 1.50 1.08
959.5 1.77 1.16
1063.5 2.01 1.12
Phase angle (deg) 77 76

Note. Calculated from magnitude data (Irvine et al. 1968a, 1968b).

Download table as:  ASCIITypeset image

Please wait… references are loading.
10.1088/0004-637X/729/2/130