Articles

PRELIMINARY RESULTS FROM NEOWISE: AN ENHANCEMENT TO THE WIDE-FIELD INFRARED SURVEY EXPLORER FOR SOLAR SYSTEM SCIENCE

, , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , , and

Published 2011 March 23 © 2011. The American Astronomical Society. All rights reserved.
, , Citation A. Mainzer et al 2011 ApJ 731 53 DOI 10.1088/0004-637X/731/1/53

0004-637X/731/1/53

ABSTRACT

The Wide-field Infrared Survey Explorer (WISE) has surveyed the entire sky at four infrared wavelengths with greatly improved sensitivity and spatial resolution compared to its predecessors, the Infrared Astronomical Satellite and the Cosmic Background Explorer. NASA's Planetary Science Division has funded an enhancement to the WISE data processing system called "NEOWISE" that allows detection and archiving of moving objects found in the WISE data. NEOWISE has mined the WISE images for a wide array of small bodies in our solar system, including near-Earth objects (NEOs), Main Belt asteroids, comets, Trojans, and Centaurs. By the end of survey operations in 2011 February, NEOWISE identified over 157,000 asteroids, including more than 500 NEOs and ∼120 comets. The NEOWISE data set will enable a panoply of new scientific investigations.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

A full understanding of the history of our solar system requires study of its most primitive and dynamic members, the asteroids and comets. These different classes of small bodies are related; recent models of planetary migration, small body collisional history, and the role that primitive bodies have played in shaping our planet's evolution predict that these signatures are imprinted on the numbers, distribution, transport, and physical properties of these objects (Gomes et al. 2005; Tsiganis et al. 2005; Morbidelli et al. 2005). Only through detailed physical studies of all the major populations of our solar system's small bodies can their origins, evolution, and fate be learned.

To date, the most prolific surveys to discover and characterize asteroids and comets have been conducted at visible wavelengths (e.g., the Catalina Sky Survey, the Lincoln Near-Earth Asteroid Research Program, the Near-Earth Asteroid Tracking Program, and Spacewatch; Larsøn 2007; Stokes et al. 2000; Helin et al. 1997; McMillan 2007). The number of in situ missions to primitive bodies is necessarily small (the Near-Earth Asteroid Rendevous mission, Giotto, Galileo, Deep Space 1, Deep Impact, Hyabusa, Dawn, Stardust, and Rosetta; Cheng 2002; Oberst et al. 2001; Reinhard 1986; Thomas et al. 1994; A'Hearn et al. 2005; Fujiwara et al. 2006; Russell et al. 2007; Brownlee et al. 2003; Barucci et al. 2005). Radar observations (Ostro et al. 2007) are a powerful means of obtaining detailed physical characterization remotely, but the objects must pass quite close to Earth. Hence, although, approximately 7600 near-Earth objects (NEOs), hundreds of thousands of Main Belt asteroids (MBAs), and ∼3000 comets have been discovered to date, detailed physical characterizations have been made for only a small fraction of these objects (e.g., Binzel et al. 2004; Wolters et al. 2008; Harris et al. 2009; Bhattacharya et al. 2010).

1.1. Infrared Properties of Primitive Bodies

Thermal infrared observations offer a powerful means of characterizing physical properties of primitive bodies, and all-sky infrared surveys offer an opportunity to obtain these data for large numbers of objects in a relatively short time. Once an orbit has been calculated for an asteroid, its distance from the Earth at any given time is known, allowing the measured flux to be related directly to the flux emitted at the body's surface. Likewise, the known distance to the Sun allows derivation of the incident solar flux. The physical surface properties of the asteroid can then be modeled to reproduce the measured emitted flux as a function of the incident radiation (Harris & Lagerros 2002; Lebofsky & Spencer 1989). Diameters derived from infrared measurements have smaller uncertainties than diameters derived from visible-light observations. Thermal flux emitted by an asteroid has a weaker dependency on a body's unknown albedo value (pv) than on reflected sunlight. With both visible and IR data, the degeneracy between diameter and albedo can be broken, and both quantities can be independently fit. Current estimates for the range of albedos for NEOs result in a factor of five uncertainty in diameter for a given reflected optical flux (Stuart & Binzel 2004), since optical flux is linearly proportional to pv. Yet, the variation of an asteroid's thermal flux is proportional to (1 − qpv)1/4, where q is the phase integral (Lebofsky & Spencer 1989). Strong constraints on diameter can still be obtained even for objects which have thermal fluxes but no visible measurements. For example, for a 0.13 km NEO with a well-known orbit observed to 5σ at 12 μm at a distance of 0.5 AU from Earth and 1.1 AU from the Sun, the diameter can be determined within 10% (Wright et al. 2010; Mainzer et al. 2011).

Because thermal flux is more weakly dependent on albedo than on visible flux, infrared surveys are less subject to the bias against low-albedo objects than against visible light surveys. Although considerable work has been done to debias visible surveys (e.g., Stuart & Binzel 2004; Jedicke & Metcalfe 1998; Harris 2008), IR surveys offer a complementary approach in determining the size and albedo distributions of small bodies in the solar system.

The Infrared Astronomical Satellite (IRAS) observed ∼1800 asteroids, most in the Main Belt (Tedesco et al. 1988; Matson 1986). Ground-based thermal infrared measurements of asteroids have been made of ∼100 objects (Delbó et al. 2003). The Spitzer Space Telescope has undertaken a number of programs to observe solar system objects, including 700 NEOs (Trilling et al. 2010). However, almost all objects observed by IRAS, Spitzer, and ground-based infrared telescopes have been optically selected, so the bias against low-albedo objects is still present, complicating interpretation of the size and albedo distributions derived from these data. A large area infrared survey capable of independent detection and discovery of a statistically significant number of objects is needed. With the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010), such a resource is now available; the WISE bandpasses cover the thermal peak of most solar system objects out to ∼9 AU. The WISE data processing pipeline has been augmented to allow identification and discovery of new moving objects independently of previous visible asteroid surveys; these modifications are collectively known as "NEOWISE." Since it is capable of independent discovery of new objects, NEOWISE is unique in its ability to carry out a survey that is much less dependent upon albedo than other surveys. The purpose of this paper is to introduce the community to the rich NEOWISE data set, including the pipeline that was used to detect moving objects, and to give an overview of the science enabled by NEOWISE. More detailed results will follow in separate papers.

1.2. The Wide-field Infrared Survey Explorer

WISE is a NASA medium-class explorer mission designed to survey the entire sky in four infrared wavelengths, 3.4, 4.6, 12, and 22 μm (denoted by W1, W2, W3, and W4, respectively; Wright et al. 2010; Liu et al. 2008; Mainzer et al. 2005). WISE is producing a multi-epoch image atlas and source catalog that will serve as an important legacy for future research. The mission's two primary science goals are to find the most luminous galaxies in the entire universe and to find the closest and coolest stars; a secondary goal is to study the nature of Earth's nearest neighbors, the asteroids and comets.

WISE consists of a 40 cm telescope in a Sun-synchronous low-Earth orbit. The telescope scans the sky continuously, pointing near the zenith while a scan mirror freezes the sky on the focal planes for approximately 11 s. While the sky is frozen on the focal plane, an 8.8 s exposure is taken of the same 47' × 47' field of view in the four WISE bands. The scan plane precesses by ∼1° day−1. Each image overlaps by 10% in the in-scan direction, with 90% overlap in the cross-scan direction. This orbit allows WISE to cover nearly every part of the sky, a minimum of eight times, with more coverage at the ecliptic poles. The survey began on 2010 January 14 and full coverage of the sidereal sky was achieved in 2010 mid-July. The telescope and detectors were maintained at their operating temperatures through the use of a dual-stage solid hydrogen cryostat. The secondary cryogen tank depleted on 2010 August 5 and the primary tank was exhausted on 2010 September 30. After that time, the first two channels (W1 and W2) continued to operate, and a four month survey (known as the NEOWISE Post-Cryogenic Mission) using these two channels with the same survey strategy was undertaken. The NEOWISE Post-Cryogenic Mission was completed on 2011 February 1.

Astrometric errors are less than 0farcs5 with respect to Two Micron All Sky Survey for objects with high signal-to-noise measurements. The preliminary estimated S/N = 5 point-source sensitivity on the ecliptic is 0.08, 0.1, 0.85, and 5.5 mJy in the four bands (assuming eight individual exposures per band; Wright et al. 2010). Sensitivity improves away from the ecliptic due to denser coverage and lower zodiacal background. For an asteroid with a bolometric albedo of 0.05, the W3 sensitivity of 0.85 mJy corresponds to an visible magnitude of 23 (Wright et al. 2010). However, as described below, the search for new minor planets was carried out using individual exposures.

The WISE survey cadence yielded excellent repeated coverage of moving objects in the solar system, with an average of 10 detections of a typical asteroid or comet spaced over ∼36 hr. This is a slightly higher number of detections than the corresponding value of eight for an inertially fixed object on the ecliptic plane because asteroids tend to move in the same direction as the WISE scan plane.

As of 2011 February, WISE detected more than 157,000 moving solar system objects. The data set contains at least 584 NEOs, ∼2000 Jupiter Trojans, 18 Centaurs and scattered-disk objects, and ∼120 comets. WISE has also detected a handful of irregular satellites of Jupiter and Saturn, including Iapetus and Phoebe, as well as cometary dust trails and dust bands in the zodiacal cloud.

2. NEOWISE: AN ENHANCEMENT TO THE BASELINE WISE MISSION

Data processing in the baseline WISE mission does not include any provision for identifying previously unknown solar system objects. Furthermore, only the co-added Atlas images that combine many exposures covering each point on the sky were planned to be publicly released. Outlier rejection algorithms designed to reduce susceptibility to transient phenomena such as cosmic rays suppress most moving objects from the co-added images. To facilitate solar system science with WISE, two enhancements were made to the baseline pipeline; these tasks are collectively called NEOWISE ("near-Earth object + WISE").

2.1. NEOWISE Task 1: Creation of the Archive of Single-exposure Images Plus Searchable Interface

NEOWISE Task 1 allows researchers to query the WISE data archive for pre-discovery astrometry and physical data for solar system objects that are discovered after the final data processing is complete. It also enables users to download the individual single-epoch images of these objects.

Recovery of small bodies in the WISE data archive requires long-term access to the time-tagged calibrated images and the database of extracted sources from the individual (or single-epoch) WISE exposures. These are referred to as the Level 1b image and source products. The Level 1b images are FITS format images in each of the four bands that have instrumental signatures removed (e.g., linearized, droop-corrected, dark-subtracted, flat-fielded, illumination-corrected, etc.) and that have the reconstructed world coordinate system and photometric zero point information in their headers. The Level 1b extracted source database contains the reconstructed equatorial positions, calibrated brightness in four bands and associated uncertainties, and assorted quality and characterization flagging for sources detected on the Level 1b images. Latent images and other instrumental artifacts are flagged in the source database. These products are produced by the scan/frame pipeline component of the WISE Science Data System.

The atlas, catalog, and assorted metadata products are archived for the long-term and served to the research community and public via the online and machine-friendly interfaces of the NASA Infrared Science Archive (IRSA). An enhanced interface to the IRSA catalog and image query services that allows searches based on solar system object name (via Horizons name resolution) or orbital elements has been developed as part of NEOWISE Task 1.

2.2. NEOWISE Task 2: Rapid Identification and Posting of Previously Unknown Solar System Object Candidates

NEOWISE Task 2 consisted of an augmentation to the ground data processing system that exploits the WISE observing cadence to identify moving objects by comparing source extractions made from imaging of the same region of sky on successive orbits. New moving object candidates were reported to the Minor Planet Center (MPC) within 10 days of the midpoint of their observations on board the flight system. The MPC serves as the central clearinghouse for NEOWISE astrometry and associates it with all other observations.

2.3. WMOPS: The WISE Moving Object Processing Software

Candidate moving object tracklets (sets of time-tagged positions) were identified using the WISE Moving Object Processing System (WMOPS). WMOPS was developed using elements of the PanSTARRS Moving Object Processing System (MOPS; Denneau et al. 2007; Grav et al. 2009; Jedicke et al. 2009) that linked detections of candidate moving objects from multiple exposures covering the same regions of the sky. WISE-specific elements of WMOPS included optimizing detection for the WISE survey cadence, adapting to WISE detector and source extractor characteristics, and specialized quality assurance filters and tools to allow vetting of moving objects prior to delivery to the MPC.

WMOPS was integrated into the general WISE data processing pipeline system and followed scan/frame processing. The scan/frame pipeline performed instrumental image calibration (linearization, dark-subtraction, flat-fielding, sky subtraction, etc.) on individual WISE images from one scan (where a scan contains ∼260 four-band exposures taken during 1/2 of a WISE orbit), detected and measured position and fluxes for sources in the calibrated images, and flagged sources that may be spurious detections of image artifacts produced by bright stars. Scan/frame processing produced photometrically and astrometrically calibrated single-exposure (Level 1b) images and extracted source lists.

WMOPS processing steps are illustrated in Figure 1. Input was the source lists from ∼90 to 120 scans (spanning approximately 3–4 days). The input source lists were first filtered to exclude low-reliability extractions such as flagged artifacts, very low signal-to-noise ratio sources, and extractions with pathologically large values of the point-spread function fitting chi-squared values. Stationary object rejection (SOR) was then performed by positionally associating the filtered source lists from overlapping images and identifying pairs of detections that do not move by more than 2 arcsec between adjacent exposures, and eliminating these detections from further processing. After SOR, the remaining moving object candidate detections are processed by the findtracklets routine that creates pairs of detections within a region of sky and correlates each pair's velocity and direction of motion, along with their sky location. The pair lists are stored in a database and then fed through the collapsetracklets routine that creates chains of pair-correlations by matching pairs according to velocity and position within a specified uncertainty threshold at each point. Candidate tracklets were required to have a minimum of five independent detections to be considered reliable.

Figure 1.

Figure 1. Flowchart summarizing WMOPS.

Standard image High-resolution image

WMOPS efficiently detected all moving objects whether they were previously known or unknown. Approximately, 80% of the tracklets were previously known as solar system objects and were automatically accepted. The remaining 20% of tracklets were validated by visual inspection by NEOWISE scientists. Figure 2 illustrates the web-based quality assurance tool that enabled rapid assessment of a tracklet's reliability. Tracklets for known objects and reliable new detections were reported to the MPC.

Figure 2.

Figure 2. Sample of the quality assurance pages used to validate a tracklet for a previously unknown moving object. The top two rows show cutouts of the images in bands W3 and W4, respectively; the following rows show W3 images of the same inertial position on the sky before and after the exposures containing the target. These rows allow verification that the source observed is not inertially fixed. The final two rows show the exposures immediately preceding the top row; this allows verification that the transient source observed in the top frame is indeed a moving astronomical object and not an artifact fixed in pixel space but not in time.

Standard image High-resolution image

WMOPS was run approximately every four days to ensure that the accumulated position uncertainties for the fastest moving objects such as NEOs remained small enough to enable efficient recovery by ground-based observers. Tracklets from each WMOPS run were sent to the MPC for publication and linkage to other observations. Objects that the MPC determined to have a high probability of being NEOs or comets were posted to the NEO confirmation page (http://www.minorplanetcenter.net/iau/NEO/ToConfirm.html) to be followed up by the world-wide network of professional and amateur astronomers. Ground-based follow-up was essential for NEOs in particular because the WISE observational arcs, while generally long enough to assign a designation to new objects, do not span a sufficient timebase to allow recovery at subsequent apparitions. Ground-based follow-up within two to four weeks allows well-determined orbits to be computed. Objects that are confirmed to be NEOs or comets are published by the MPC via the Minor Planet Electronic Circulars (MPECs) or IAU Circulars. The astrometry and computed orbits of slower moving objects were published on a monthly basis in the MPC's Minor Planet Orbit Supplement.

2.3.1. Range of Motion Limits

The minimum tracklet length validation requirement and the minimum sky-plane projected velocity imposed by the SOR define the range of motion to which WMOPS is sensitive. Parallax motions below 2 arcsec over 90 minutes (approximate time between WISE coverages near the ecliptic) correspond to a distance of roughly 70 AU. The limits of five detections impose an upper limit in that objects which may outpace the survey cadence before the minimum number of detections can be obtained. Motions greater than the upper limits perpendicular to the scan direction are roughly 3.fdg25 day−1, considering both an irregular cadence pattern as well as an evenly progressing survey. The maximum speed of an object detected by WMOPS was 3.fdg3 day−1, but the component perpendicular to the scan direction was 3.fdg22 day−1.

3. RESULTS

Figure 3 shows a top–down view of all of the objects detected by NEOWISE as of 2011 February. Figures 4(a) and (b) show the distribution of orbital elements of all the solar system objects identified by NEOWISE to date, including NEOs, MBAs, comets, Trojans, and Centaurs. By 2011 February, more than 252,000 tracklets were reported to the MPC during the cryogenic and post-cryogenic missions, representing observations of ∼157,000 unique minor planets.

Figure 3.

Figure 3. WISE always surveys near 90° solar elongation; this figure shows a top–down view of the objects detected by NEOWISE as of 2011 February (distances are given in AU). The outermost circle represents Jupiter's orbit; the interior circles represent the terrestrial planets. Previously known NEOs are shown as green circles; new NEOs discovered by NEOWISE are shown as red circles; previously known comets observed by WISE are shown as cyan squares, and comets discovered by NEOWISE are shown as yellow squares. All other objects are shown as black points. The drop in density of objects observed near (+2, +2) AU in the figure is due to the exhaustion of the secondary tank's cryogen on 2010 August 5, resulting in the loss of band W4. The dashed red line indicates the exhaustion of the primary tank and the start of the NEOWISE Post-Cryogenic Mission on 2010 October 1; the survey was completed on 2011 February 1. The drop in detections near (+2, −2) AU in the figure is due to the intersection of the galactic plane with the ecliptic plane; the higher backgrounds and confusion caused by galactic cirrus resulted in the identification of fewer sources.

Standard image High-resolution image
Figure 4.

Figure 4. (a) NEOWISE has detected objects throughout the solar system. The orbital elements of all objects that NEOWISE has detected in the inner solar system with arcs longer than 30 days are shown as black dots; objects with arcs less than 30 days are shown in cyan. The quantization seen in the cyan points is due to the method by which the MPC computes orbits for objects with short observational arcs. (b) a–i plot of objects observed by NEOWISE in the inner solar system; color scheme is the same as in Figure 4(a).

Standard image High-resolution image

3.1. Near-Earth Objects

NEOs are thought to be a transient population, persisting in their orbits for only a few million years before being ejected into the outer solar system, crashing into the Sun, or colliding with other bodies (Gladman et al. 2000). Therefore, they must be continually resupplied from more permanent source reservoirs such as the Main Belt and comets. Recent models have improved our understanding of the transport mechanisms and sources of NEOs (e.g., Bottke et al. 2002) based on observations of their size frequency distribution (SFD). However, current NEO SFD models rely on the absolute magnitude (H) through an assumed albedo to derive size. By observing and discovering a sample of NEOs at thermal infrared wavelengths, we can obtain an improved computation of the SFD because (1) IR-derived diameters are intrinsically more accurate; (2) we can obtain albedos and further constrain diameters for all individual objects that have measured optical magnitudes; and (3) an IR survey is less biased against low-albedo objects. The improved SFD offered by NEOWISE will allow us to test NEO population models, to study the collisional histories of NEO and MBA populations, as well as to examine how SFD and albedo vary with asteroid collisional families.

Our understanding of the hazard posed by NEOs also depends on knowledge of their SFD, and efforts to plan mitigation strategies depend on the physical properties of NEOs. To the extent that albedo can be used as a proxy for asteroid composition, it will be possible to derive an approximate density, mass, and estimated impact energy for previously known and new objects observed by WISE.

As of 2011 February, NEOWISE has identified more than 584 NEOs, including discovery of 135 previously unknown NEOs. Figure 5 shows the preliminary distribution of H magnitudes for all NEOs known to have been observed by NEOWISE during the cryogenic mission. The distribution peaks at H ∼ 18, corresponding to diameters between ∼470 and 1500 m for albedos between 0.05 and 0.5 (H values for previously known objects were taken from the MPC observation file; H values for newly discovered objects were obtained by visible follow-up observations from a world-wide network of amateur and professional observers). The survey cadence includes many repeated coverages at the ecliptic poles, regions which most ground-based surveys preferentially avoid. NEOWISE is thus sensitive to objects at extreme declinations, which are more likely to have high orbital inclinations or have orbits that bring them very close to the Earth (Figure 6). Furthermore, unlike ground-based telescopes, the survey is unimpeded by daytime, weather, or variable seeing, so it is easier to debias (cf. Jedicke & Metcalfe 1998). However, of the 135 new NEOs discovered, ∼15 NEOs were designated by the MPC, but received no optical follow-up either due to their optical faintness or unavailability of suitable ground-based facilities. Due to the ∼36 hr arc length typical of most NEOWISE detections, ground-based follow-up was necessary in most cases to obtain a secure orbit for the newly discovered NEOs. An additional ∼15 NEO candidates were placed on the MPC's NEO confirmation page, but received neither designations nor follow-up. A radiometric diameter can be determined for these objects (subject to the uncertainties in Sun–object–observer distances due to orbital uncertainties), but albedo cannot be computed without optical fluxes. It is hoped that these objects will be recovered either from archival data or future surveys. Without this, these objects will result in an additional uncertainty in the derived NEO albedo distribution.

Figure 5.

Figure 5. H magnitude distribution for the >500 NEOs detected by WISE as of 2011 February.

Standard image High-resolution image
Figure 6.

Figure 6. WISE's survey cadence results in increased coverage at the ecliptic poles, making it easier to detect NEOs at high inclinations. The preliminary, non-debiased NEO inclination distribution is shown here for all NEOs observed by WISE. The data set must be debiased before conclusions about the NEO inclination distribution can be drawn.

Standard image High-resolution image

3.2. Main Belt Asteroids

NEOWISE has observed over 154,000 MBAs as of 2011 February, including ∼33,000 new discoveries. This represents an advance of nearly two orders of magnitude over the IRAS Minor Planet Survey (Tedesco et al. 2002), allowing a correspondingly more detailed view of the size and albedo distributions throughout the Main Belt. Figure 7 shows that the preliminary distribution of H magnitudes for all MBAs observed by WISE by the end of the cryogenic mission peaks at ∼15, corresponding to diameters between 1.9 and 5.9 km for albedos between 0.05 and 0.5. H values were taken from the MPC observation files for previously known objects; however, the H values for the newly discovered objects are approximated from the objects' thermal fluxes since no dedicated visible follow-up for these objects was carried out. It is expected that linkages will eventually be made between visible survey data and WISE observations, resulting in increasing numbers of NEOWISE-discovered objects with well-measured H values and refined orbits. Figure 8 shows the preliminary H magnitude distribution as a function of semimajor axis within the Main Belt. With thermal infrared fluxes available for all of these Main Belt asteroids, it is now possible to derive diameters and (where optical fluxes are available) albedos. This will facilitate study of the distribution of sizes and albedos within dynamical families (Hirayama 1918) as has been done with multi-band optical photometry (Parker et al. 2008; Ivezić et al. 2002).

Figure 7.

Figure 7. Preliminary distribution of absolute magnitudes of ∼150,000 Main Belt asteroids detected by NEOWISE peaks at H ∼ 15; this corresponds to diameters between 1.9 and 5.9 km for albedos of 0.5–0.05.

Standard image High-resolution image
Figure 8.

Figure 8. Peak of the Main Belt asteroid H magnitude distributions shifts as a function of semimajor axis. The plot shows the distributions of four categories of Main Belt asteroids, as defined by Zellner (1979): MBA I (2.06 < a ⩽ 2.5), MBA II (2.5 < a ⩽ 2.82), MBA III (2.82 < a ⩽ 3.27), and MBA IV (3.27 < a ⩽ 3.65). The solid line represents the histogram of MBA I, the dotted line represents MBA II, the dashed line is MBA III, and the dash-dotted line represents MBA IV. The histogram numbers of MBA IV have been multiplied by 10 for clarity.

Standard image High-resolution image

The NEOWISE data set contains observations of asteroids that have been previously associated with one of the 54 dynamical families defined by Nesvorny et al. (2006a). Figure 9 shows the number of objects observed by NEOWISE for each family; for some families, NEOWISE has observed thousands of known members. The combination of WISE infrared fluxes with multi-band optical photometry should allow new family members to be identified on the basis of their albedos and optical colors as well as their dynamical properties. Parker et al. (2008) have shown that it is possible to use Sloan Digital Sky Survey (SDSS) colors to find two to three times more family members that make up "halos" surrounding existing families. These are family members that have moved far enough into the background population that they can no longer be identified by hierarchical clustering methods, yet they have colors identical to the cores of the families (e.g., Vesta and Eos). In particular, this technique may facilitate the association of the smallest, least massive objects with various families. These objects are of interest because their small sizes make them more susceptible to the non-gravitational forces that cause the most spread in the orbital elements which are used to compute family ages (cf. Nesvorny et al. 2006a). By incorporating size and albedo distributions derived from infrared fluxes, as well as by finding additional family members using their albedos and colors, the NEOWISE data set should facilitate an improved understanding of the ages and timing of the breakups that created the present-day Main Belt families.

Figure 9.

Figure 9. Approximately, 26,000 of the Main belt asteroids observed by NEOWISE during the cryogenic mission are objects that are known to belong to dynamical families Nesvorny et al. (2006a). In addition to these objects, it is highly likely that other previously unrecognized family members are distributed among the more than 150,000 other MBAs detected by NEOWISE. These objects may be identified by their IR-derived albedos and optical colors. This will facilitate a wide range of studies aiming to improve the precision of calculations of parent body breakups.

Standard image High-resolution image

3.3. Comets

As reservoirs of some of the most pristine primordial material in the solar system, comets serve as Rosetta stones for understanding the transport and history of volatiles. Comets sample the most distant regions of the solar system: the Kuiper Belt just outside of Neptune's orbit, and the Oort cloud, which may extend halfway to the nearest star. While we cannot directly observe the earliest stages of comet formation in our own solar system, the SFD and orbital distribution of present-day comets preserve a record of the outer solar nebula's mass distribution that can be used to model the collisional and dynamical evolution of comet source regions over the last 4.5 Gyr (cf. Tsiganis et al. 2005; Levison et al. 2009). The cometary size and orbital element distributions allow a test of models of the origins of Halley family comets (HFCs) and long-period comets (LPCs) as well as Jupiter family comets (JFCs). Low-inclination JFCs are thought to originate in the Kuiper Belt and are expected to have smaller sizes than the HFCs and LPCs due to collisional processes and mass loss caused by their numerous passages by the Sun (Weissman & Levison 1997). The HFCs and LPCs, believed to originate in the Oort cloud, should have retained their original size distribution (Stern & Weissman 2001) and be larger (Meech et al. 2004). Furthermore, the distribution of perihelia can set constraints on the number and mass of objects in the Oort cloud (Francis 2005).

WISE has observed ∼120 comets as of 2011 February, of which 20 are new discoveries (including three previously known asteroids on which WISE observed activity; these objects were subsequently redesignated as comets). The NEOWISE survey of comets offers a data set that is less biased against low-albedo objects than previous optical surveys, and as an all-sky survey it has covered high-declination regions as well as the ecliptic plane with a well-determined set of survey pointings. Detections of previously known comets were identified by WMOPS, and new discoveries were made by inspecting the ∼20% of tracklets that were not associated with known objects to look for evidence of extended emission. New discoveries suspected of being comets were posted to the MPC's NEO confirmation page. All of the 20 new discoveries received visible follow-up, allowing orbits to be refined beyond the initial orbits computed solely from NEOWISE observations. As with the other major populations detected by WISE, infrared fluxes can be used to compute nucleus size when cometary activity does not significantly obscure it. While the Spitzer Space Telescope has undertaken a survey of 100 comets (Fernandez et al. 2008), they are all JFCs that have been selected by optical surveys. The NEOWISE survey includes both long- and short-period comets, as well as independent discoveries, allowing constraints to be set on the ratio of long- to short-period comets. Of the ∼120 comets observed by the end of the cryogenic portion of the mission, 35% are LPCs. Figures 10(a) and (b) show the distribution of orbital elements of comets detected by NEOWISE during the cryogenic mission. Many of the comets have high inclination orbits; however, the survey cadence resulted in increased coverage at the ecliptic poles, leading to a bias in favor of detecting comets at high inclinations (or closer to Earth). Of the 20 comets discovered by NEOWISE as of 2011 February, approximately one-third have perihelia greater than 4 AU, and the remaining two-thirds are LPCs, four of which have nearly parabolic orbits.

Figure 10.

Figure 10. (a) and (b) Preliminary, non-debiased distribution of orbital elements of NEOWISE-detected comets reveals many with high inclinations.

Standard image High-resolution image

For objects where significant cometary activity is observed, it is possible to compute particle size distributions, allowing the strength of non-gravitational forces to be determined as a function of particle size. For comets with sufficiently high fluxes at bands W1 (3.4 μm) and W2 (4.6 μm), gas ratios of CO and CO2 can be computed.

NEOWISE has also observed four objects that have been targeted by spacecraft for in situ visits, 19P/Borrelly (Deep Space 1), 9P/Tempel 1 (Deep Impact), 103P/Hartley (EPOXI; Bauer et al. 2010), and 67P/Churyumov–Gerasimenko, which the Rosetta mission will encounter in 2014. Figures 11(a)–(c) show the WISE views of 67P/Churyumov–Gerasimenko at two apparitions. A preliminary fit of a blackbody function to the first epoch of observations on 2010 January 19 yields an effective temperature of 179 K using an aperture radius of 22 arcsec. A similar fit to the second epoch observations on 2010 July 4 yields 155 K (Figure 11(c)). These fitted temperatures exceed ambient blackbody temperatures by about 20 K. It is not clear that this is due to the dust grain temperature, as there is a silicate emission band near the center of WISE band W3 at ∼10–12 μm.

Figure 11.

Figure 11. (a)–(c) Rosetta target 67/P Churyumov-Gerasimenko. (a) and (b) Three-color WISE images of the comet observed at two epochs: 2010 January 19 (left panel a) at a heliocentric distance R = 3.32 AU and an observer distance Δ = 3.14 AU, and 2010 June 29 (right panel b) with R = 4.18 AU, Δ = 3.97 AU. In these images, the blue channel is band W2 (4.6 μm), green is band W3 (12 μm), and red is band W4 (22 μm). (c) Blackbody fits to 67/P for the two epochs of WISE observations. Preliminary corrections to fluxes were made by interpolating between the blackbody color corrections provided in Wright et al. (2010). The blackbody fit for the first epoch of observations is shown as the solid black line, and the WISE data points are red diamonds; the blackbody fit for the second observational epoch is the dashed black line, and the data points are blue triangles. 67/P was undetected in bands W1 and W2 in both epochs.

Standard image High-resolution image

3.4. Jupiter Trojans

Trojans are populations of minor planets that are co-orbital with a planet. The known Jupiter Trojans (there are also a handful of objects that share orbits with Mars and Neptune) are separated into two clouds of objects that librate around the L4 and L5 Lagrangian points in Jupiter's orbit. These two points are found 60° of heliocentric ecliptic longitude ahead (also called the leading cloud) and behind (the trailing cloud) the planet, and both clouds possess enough dynamical stability to survive over the age of the solar system (Marzari & Scholl 1998). Reviews of the Trojan populations can be found in Jewitt et al. (2004) and separately in Dotto et al. (2008). The Trojans are an important population for understanding the origin and evolution of Jupiter, and they also provide a potential window into the solar nebula at the distances of the giant planets, since they are the only major population that might have survived the dynamical clearing of minor planets in this region during the early formation of the solar system. But how and when the Trojans were trapped in Jupiter's 1:1 mean-motion resonance remains unknown. Whether the capture happened during a very early epoch while the planet was forming or much later during the dynamical clearing phase is still under discussion (Morbidelli et al. 2005; Marzari & Scholl 2007).

There are currently about 4700 known Jupiter Trojans, and NEOWISE has reported observations of ∼1600 of these, increasing the number of thermal emission measurements of this population by more than an order of magnitude compared to earlier surveys (Cruikshank et al. 1977; Tedesco et al. 2002; Fernandez et al. 2003, 2009). In addition, NEOWISE detected an additional 400 objects that the MPC has designated with Trojan-like orbits, but these have arcs that are too short for their identification to be secure. WISE observes the sky in scans of constant ecliptic longitude, and as the survey began, the first scans were already in the middle of the leading Trojan cloud. The scans moved toward the tail of the cloud, resulting in about 20%–30% of the leading cloud closest to Jupiter remaining unobserved during the cryogenic survey (although it was observed about seven months later with bands W1 and W2 during the post-cryogenic mission). The survey exited the leading cloud at the end of 2010 February and almost immediately entered the trailing cloud on the opposite side of the solar system (Figure 3). Observations of the trailing cloud were completed in late 2010 May. The tail end of the trailing cloud was close to the galactic center when it was observed by NEOWISE, leading to noticeable loss of detection efficiency for less than 10% of the trailing cloud. While the two clouds were not observed with perfect uniformity, NEOWISE has provided the most uniform survey of the Jupiter Trojan population to date. Of the ∼1600 objects that have been conclusively identified as Trojans, ∼900 were observed in the leading cloud and ∼700 in the trailing cloud. This yields a non-debiased ratio of the number of objects in the L4 and L5 cloud of ∼1.3. The 400 probable Trojans are distributed roughly equally among the two clouds. The H-distribution (Figure 12) shows that the survey's limit is H ∼ 13 (or ∼12 km for an assumed albedo of 0.07).

Figure 12.

Figure 12. NEOWISE distribution of observed Trojans peaks at H ∼ 13. Of the 4237 Trojans known, WISE has observed ∼1600 during the one year survey (solid black line), along with ∼400 previously unknown Trojans for which orbits are presently poorly constrained. More Trojans were observed in the L4 cloud that leads Jupiter (red dashed line) than the L5 trailing cloud (blue dash-dotted line), despite the fact that the cryogenic portion of the survey ended before the survey of the L4 cloud was complete.

Standard image High-resolution image

The survey provides a unique opportunity to not only characterize the physical properties of a large fraction of the Trojans, but it also provides a uniform survey of both clouds. The survey must be debiased in order to test the initial results that show an overabundance of objects in the leading cloud compared to the trailing cloud. In addition, a proper vetting of the short arc sample using the on-sky velocity (similar to the method used on the SDSS data by Szabo et al. 2007) and thermal colors is needed to define the Trojan sample. As with the Main Belt asteroids, NEOWISE detections of new Trojan candidates may be linked to other observations in the future, allowing refinement of their orbits.

3.5. Unusual Objects and Outer Solar System Objects

WISE's wavelengths are not long enough to detect large number of objects in the outer solar system via their thermal emission. Nonetheless, several dozen objects were observed and discovered in this region with the NEOWISE pipeline. A number of unusual objects with high inclinations have been discovered by NEOWISE such as 2010 OA101, an outer Main Belt asteroid with an inclination of 84°, and 2010 CR140, an object with a semimajor axis of 5.6 AU and an inclination of 74.fdg7. 2010 CR140 may represent an object in transition between the Jupiter Trojan and Centaur populations. In total, 11 objects were observed in the outer solar system with inclinations larger than 80°. Figure 13 shows an a–e plot of WISE-observed objects, highlighting the outer solar system. WISE has observed ∼20 Centaurs to date, including five new discoveries.

Figure 13.

Figure 13. a–e distribution of objects observed by NEOWISE in the outer solar system.

Standard image High-resolution image

3.6. Light Curve Studies and Thermophysical Modeling

By obtaining an average of 10 observations spaced uniformly over ∼36 hr for each solar system object, the WISE cadence lends itself to the study of small body light curves. Ground-based observers are interrupted by daytime, clouds, and seeing variations; however, WISE is able to provide an extremely uniform, stable platform with regular observations spanning days or even, in a few cases, weeks. An example of this type of long-period light curve for a near-Earth asteroid is shown in Figure 14. This makes the WISE data set unique in comparison with ground-based light curve surveys that are generally limited to studying rotational periods faster than ∼1 day (Masiero et al. 2009; Pravec et al. 2008, 2006). The WISE observations offer a unique window into long-period light curve studies, and with suitable debiasing, the fraction of objects with rotational periods longer than a day may be determined. Figure 15 shows the distribution of the number of observations for objects detected by WMOPS; on average, a typical solar system object was detected every three hours by NEOWISE.

Figure 14.

Figure 14. 12 μm time series of the near-Earth asteroid 230118 (2001 DB3) extends over two days. The rotational modulation is readily visible due to WISE's unique observing cadence.

Standard image High-resolution image
Figure 15.

Figure 15. On average, WMOPS detected most moving solar system objects ∼10–12 times over 1.5 days; however, some NEOs have been observed more than a hundred times.

Standard image High-resolution image

Binary asteroids have been found via photometric observations, direct imaging (ground-based and space-based), adaptive optics, and radar. Binary systems provide a useful laboratory because they allow masses to be determined; when combined with infrared diameters, density and porosity may be computed. It is possible to search for binary systems in the WISE data set by searching for objects with light curve amplitude variations ≳1.5 mag; such objects are the most likely to be binary systems (Chandrasekhar 1969; Holsapple 2001). An example of a Main Belt asteroid displaying the large-amplitude light curve indicative of a candidate binary system is shown in Figure 16. With a period of more than 2.5 days and an amplitude of ∼1.5 mag, this object, (49137) 1998 SC35, will require a concerted follow-up observing campaign to confirm its binarity. Preliminary fits to the individual WISE data points using a thermophysical model assuming a single spherical object and the Near-Earth Asteroid Thermal Model (NEATM; Harris 1998) yield diameters ranging from 3.8 to 7.5 km; the wide range of diameters is due to rotational effects and is much larger than the measurement uncertainties.

Figure 16.

Figure 16. 12 μm light curve of Main Belt Asteroid 49137 extends over 2.5 days and shows an amplitude variation large enough (∼1.5 mag) to indicate that the object may be a binary. The bottom of the light curve shows a flattened region that may indicate the signature of the eclipsing secondary. The object was detected by WMOPS beginning on MJD 55306.640.

Standard image High-resolution image

As with visible data, shape and rotation state can be derived from WISE observations (cf. Delbó et al. 2003; Wright et al. 2007a, 2007b). If sufficient data are available at a range of different phase angles, it is possible to compute the pole orientation and rotation state of an object in addition to its shape model through light curve inversion (Kaasalainen et al. 2002; Müller et al. 2005; Ďurech et al. 2010). However, unlike models computed with visible data alone, the infrared fluxes allow the absolute size of the object to be determined as well. For some WISE-observed objects, it may be possible to derive an albedo map across an object's surface.

3.7. Comet Trails and Dust Bands

WISE thermal observations are sensitive to extended structures related to the principal mechanisms of dust production in the inner solar system: comet dust trails and zodiacal dust bands. Dust trails are narrow (arcminute) structures extending from degrees to tens of degrees in the sky, and they consist of low-velocity millimeter to centimeter diameter particle emissions from short-period comets (Sykes et al. 1986b; Sykes & Walker 1992). Spitzer Space Telescope observations have shown that 80% of the short-period comets surveyed have trails (Reach et al. 2007). A cratering or disruptive collisional event between two asteroids would also result in the formation of a dust trail, though the ejection velocities (compared to trail-forming comet emissions) are expected to be higher, resulting in a broader trail (e.g., Sykes 1988; Nesvorny et al. 2006b). Over time, gravitational perturbations acting on the trail particles result in the differential precession of their nodes, resulting in the formation of a torus of collisional material extending around the Sun. Some of these bands have been associated with the breakups of asteroids in the Main Belt that lead to the formation of young asteroid families (Nesvorny et al. 2003, 2006b). Identifying the dust bands associated with young families will serve to confirm their age and will allow quantitative studies of the dispersal time for the bands (Espy et al. 2009). WISE has already detected a number of previously known trails such as that created by Comet 65/P Gunn (Figure 17), which was discovered with IRAS and observed by Infrared Space Observatory (Sykes et al. 1986a; Colangeli et al. 1998). The WISE image is approximately a hundred times more sensitive than IRAS at 22 μm and has significantly better spatial resolution. Future work should allow discovery of more trails that can potentially be associated with either asteroid families or terrestrial meteor streams.

Figure 17.

Figure 17. Jupiter-family Comet 65P/Gunn was observed by WISE on 2010 April 24. This single-frame observation shows 3.4 μm as blue, 4.6 μm as green, 12 μm as orange, and 22 μm as red. The nucleus and tail of the comet are clearly visible, along with the trail of debris that has been spread around the comet's orbit.

Standard image High-resolution image

4. CONCLUSIONS

With NEOWISE, we have detected over 157,000 minor planets during a year long survey, including discovery of ∼33,000 new minor planets. We have shown that the WMOPS system is capable of rapid identification of moving objects over a wide range of apparent velocities when provided with appropriately calibrated input images and source detections. The unique nature of the NEOWISE data set will facilitate a wide range of small body studies from fields as diverse as population studies, infrared light curve shape modeling, thermophysical studies, comets, and dust bands and comet trails. With close to 100 times more objects observed than its predecessor, IRAS, NEOWISE will leave an even greater legacy for future researchers.

This publication makes use of data products from WISE, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. This publication also makes use of data products from NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology, funded by the Planetary Science Division of the National Aeronautics and Space Administration. We gratefully acknowledge the extraordinary services specific to NEOWISE contributed by the International Astronomical Union's Minor Planet Center, operated by the Harvard-Smithsonian Center for Astrophysics, and the Central Bureau for Astronomical Telegrams, operated by Harvard University. We appreciate the contributions of the referee, J. Emery, which greatly improved the manuscript. We also thank the worldwide community of dedicated amateur and professional astronomers devoted to minor planet follow-up observations. This research has made use of the NASA/IPAC Infrared Science Archive, which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. We thank the Auton Lab of Carnegie Mellon University for use of their libraries, and the PanSTARRS project for MOPS.

Please wait… references are loading.
10.1088/0004-637X/731/1/53