Issue |
A&A
Volume 503, Number 2, August IV 2009
|
|
---|---|---|
Page(s) | L17 - L20 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200912374 | |
Published online | 27 July 2009 |
LETTER TO THE EDITOR
Spitzer observations of spacecraft target 162173 (1999 JU3)
H. Campins1,2 - J. P. Emery3 - M. Kelley4 - Y. Fernández1 - J. Licandro2 - M. Delbó5 - A. Barucci6 - E. Dotto7
1 - Instituto de Astrofísica de Canarias, c/Vía Láctea s/n, 38200 La Laguna, Tenerife, Spain
2 - Physics Department, University of Central Florida, Orlando, FL 32816, USA
3 - University of Tennessee, Knoxville, Tennessee
4 - Department of Astronomy, University of Maryland, College Park, MD 20742-2421, USA
5 - UNS, CNRS, Observatoire de la Côte d'Azur, Nice, France
6 - LESIA, Observatoire de Paris, France
7 - INAF, Osservatorio Astronomico di Roma, Roma, Italy
Received 22 April 2009 / Accepted 20 July 2009
Abstract
Context. Near-Earth asteroid 162173 (1999 JU3) is the primary target of the Japanese Aerospace Exploration Agency (JAXA) Hayabusa-2 sample return mission, and is also on the list of potential targets for the European Space Agency (ESA) Marco Polo sample return mission. Earth-based studies of this object are fundamental to these missions.
Aims. Our aim is to provide new constraints on the surface properties of this asteroid.
Methods. We present a mid-infrared spectrum (5-38 m) obtained with NASA's Spitzer Space Telescope in May 2008 and results from the application of thermal models.
Results. These observations place new constraints on the surface properties of this asteroid. To fit our spectrum we used the near-Earth asteroid thermal model (NEATM) and the more complex thermophysical model (TPM). However, the position of the spin-pole, which is uncertain, is a crucial input parameter for constraining the thermal inertia with the TPM; hence, we consider two pole orientations. First is the extreme case of an equatorial retrograde geometry from which we derive a rigorous lower limit to the thermal inertia of 150 Jm-2 s
.
Second, when we adopt the pole orientation of Abe et al. (2008a,
COSPAR Scientific Assembly) our best-fit thermal model yields a value for the thermal inertia of
Jm-2 s-0.5 K-1 and even higher values are allowed by the uncertainty in the spectral shape due to the absolute flux calibration. Our best estimates of the diameter (
km) and geometric albedo (
)
of asteroid 162173 are consistent with values based on previous mid-infrared observations.
Conclusions. We establish a rigorous lower limit to the thermal inertia, which is unlikely but possible, and would be consistent with a fine regolith similar to wthat is found for asteroid 433 Eros. However, the thermal inertia is expected to be higher, possibly similar to or greater than that on asteroid 25143 Itokawa. An Accurately determining the spin-pole of asteroid 162173 will narrow the range of possible values for its thermal inertia.
Key words: minor planets, asteroids - infrared: solar system - space vehicles
1 Introduction
Near-Earth asteroid 162173 (1999 JU3) is the primary target of the Japanese Aerospace Exploration Agency (JAXA) Hayabusa-2 sample return mission. This object is also on the list of potential targets for the European Space Agency (ESA) Marco Polo mission, which was selected for an assessment phase study in late 2007. The goal of both missions is to return to Earth samples of a primitive near-Earth asteroid for laboratory analysis.
Since these space missions seek to understand the origin and nature of volatile and organic material in the early Solar System, samples of a primitive asteroid (C, P, or D asteroids in the Tholen classification system; Tholen & Barucci 1989) are highly desirable. Asteroid 162173 is a C-class object, more specifically a Cg-type in the SMASS classification system (Binzel et al. 2001). In the asteroid population, C-class objects are abundant in the middle to outer Main Belt, and are primitive, volatile-rich remnants of the early solar system. Spectroscopy of C-class objects and comparisons with potential meteorite analogs suggest surface constituents of anhydrous silicates, hydrated minerals, organic polymers, magnetite, and sulfides (e.g., Rivkin et al. 2002; Gaffey et al. 2002). Dynamically, asteroid 162173 is most likely to have come from the main asteroid belt through the 3:1 jovian mean motion resonance (e.g., Bottke et al. 2002). At visible and near infrared wavelengths (0.5 to 2.0 m), this asteroid shows a spectrum with a neutral to bluish spectral slope (Binzel et al. 2001; Abe et al. 2008b). A broad and shallow absorption, centered near 0.7
m and associated with hydrated minerals (iron-bearing phyllosilicates), has been reported at one rotational phase by Vilas (2008). This asteroid appears to be nearly spherical with a rotational amplitude of about 0.1 mag (Abe et al. 2008b) and a rotation period of 0.3178 days (Hasegawa et al. 2008). Ground-based and space-based narrowband mid-IR photometry were used by Hasegawa et al. (2008) to derive estimates of the geometric albedo (0.063 +0.020 -0.015), diameter (
km), and thermal inertia (
500 Jm-2 s-0.5 K-1). The small rotational amplitude makes determining a spin-pole orientation for this asteroid a challenging task; nevertheless, a spin-pole orientation has been estimated at ecliptic longitude and latitude of 331 and 20
,
respectively (Abe et al. 2008a). In this work, we present new mid-infrared spectroscopy (5-38
m) obtained with NASA's Spitzer Space Telescope, which has allowed us to place new constraints on the surface properties of this asteroid.
2 Observations and data reduction
The spectrum was obtained with the Infrared Spectrograph (IRS) instrument (Houck et al. 2004) on NASA's Spitzer Space Telescope (Werner et al. 2004), on UT May 2.084, 2008. The asteroid was observed at a heliocentric distance of 1.202 AU, at a Spitzer distance of 0.416 AU, a phase angle of 52.6,
and a galactic latitude of 3.2
.
The IRS measures spectra over the range 5.2-38
m. The observations presented here used the low spectral resolution mode (
-128), which covers the full spectral range in four long-slit segments, and the final spectrum corresponds to a total on-object exposure time of 738.2 s. These four segments are short wavelength, 2nd order (SL2; 5.2-8.5
m), short wavelength, 1st order (SL1; 7.4-14.2
m), long wavelength, 2nd order (LL2; 14.0-21.5
m), and long wavelength, 1st order (LL1; 19.5-38.0
m). When each of the low-resolution orders is centered on the target, the other is offset (
50'' for SL and
108'' for LL) from the target, providing a measurement of the sky background. We anticipated a crowded field, so we also obtained the sky background from a ``shadow'' observation, which duplicated the first observation but with the asteroid out of the field of view. We did not use the peak-up routine because (a) the ephemeris uncertainty was only about 0.3 arcsec, much less than the slit widths; (b) the Spitzer pointing accuracy is about 1 arcsec, also less than the slit widths; and (c) the object's Galactic latitude was low so we did not want IRS to inadvertently peak-up on a background source. Our spectra were cleaned, extracted, and co-added using the Spitzer IRS Custom Extraction tool (SPICE), and our own routines. More details on our data reduction can be found in Kelley et al. (2009).
When the reduction was completed, we found a discrepancy of 10% between the fluxes of overlapping wavelengths in the SL and LL orders, which is equal to the expected uncertainty in the absolute flux calibration (Houck et al. 2004), and small offsets like this one are not unusual. We scaled the LL and SL fluxes to match each other in the overlapping wavelengths, the SL flux by 1.05 and the LL by 0.95, and we include this 10% systematic absolute calibration uncertainty in the modeled parameters discussed in the next section. This spectrum is plotted in Fig. 1 along with the best-fit thermophysical model (TPM, Sect. 3).
![]() |
Figure 1: The Spitzer spectrum of near-Earth asteroid 162173 (1999 JU3) plotted along with our best-fit thermophysical model (TPM). The model parameters are given in Table 1 and discussed in Sect. 3. The LL and SL fluxes have been scaled to match each other in the overlapping wavelengths (see Sect. 2). |
3 Thermal modeling
The flux from near-Earth asteroids in the wavelength range covered by the IRS is dominated by thermal emission. The measured spectral energy distribution (SED) depends on the object's size, composition, and temperature distribution. This last term is dependent on several factors, including distance from the Sun, albedo, thermal inertia, surface roughness, rotation rate, shape, and spin-pole orientation. We used two models to fit our 5 to 38 m spectrum, the relatively simple near-Earth asteroid thermal model (NEATM, Harris 1998) and the more complex thermophysical model (TPM, e.g., Harris & Lagerros 2002 and references therein).
The NEATM is a refinement of the standard thermal model (STM, Lebofsky et al. 1986; Lebofsky & Spencer 1989), which was developed and calibrated for main-belt asteroids. The NEATM accounts for observations at larger phase angles and contemplates possible deviations from the STM in the thermal behavior of smaller asteroids. Unlike the STM, the NEATM requires observations at multiple wavelengths and uses this information to force the model temperature distribution to be consistent with the apparent color temperature of the asteroid. The NEATM solves simultaneously for the beaming parameter ()
and the diameter (D). The beaming parameter was originally introduced in the STM to allow the model temperature distribution to fit the observed enhancement of thermal emission at small solar phase angles due to surface roughness. In practice,
can be thought of as a modeling parameter that allows a first-order correction for any effect that influences the observed surface temperature distribution (such as beaming, thermal inertia, and rotation). The values of
derived for asteroid 162173 using the NEATM fits to our spectrum are given in Table 1; for all our thermal models we assumed a bolometric emissivity (
)
of 0.9. These
values can be compared with those of other asteroids studied in a similar manner, and in fact fit well the trend of increasing
(lower thermal fluxes and cooler color temperatures) with increasing solar phase angle, observed for about 40 near-Earth asteroids so far (e.g., Delbó 2004; Campins et al. 2009). We did not make any observations of this asteroid's reflected continuum, so to estimate the geometric albedo (
), we adopted the values given by Hasegawa et al. (2008) of
and
,
where
is the visible absolute magnitude and G is the slope parameter.
Table 1: NEATM and TPM best-fit models.
In Table 1, scaling means that a correction was applied to the fluxes in the LL and SL orders (Sect. 2). Diameter, ,
,
and
are the derived diameter, beaming parameter, thermal inertia, and geometric albedo, respectively. TPM1 yields a rigorous lower limit to the thermal inertia (Sect. 3) and TPM2 assumes the pole orientation of Abe et al. (2008a) given in Sect. 1. Uncertainties in the diameter estimated with thermal modeling usually exceed the formal errors and are typically 10-15% (e.g., Delbó et al. 2003). The listed error in the geometric albedo includes the uncertainty in
,
with a 0.05 mag component due to the rotational amplitude of 0.1 mag, and the diameter uncertainty propagated into the albedo.
The TPM is applied to asteroids with known shape and spin state, and is generally used to constrain the thermal inertia (), which is particularly relevant to the proposed sample-return missions since
can be diagnostic of the regolith. For example, typical values of
(in S.I. units Jm
)
are 50 for lunar regolith, 400 for coarse sands, and 2500 for bare solid rocks (e.g., Delbó & Tanga 2009, and references therein). As mentioned in Sect. 1, there are estimates of the rotation period, approximate spin-pole orientation and the shape of asteroid 162173, so in principle we can apply the TPM. However, because of the uncertainties in the pole position, the TPM does not uniquely constrain the thermal inertia of this asteroid; we elaborate on this point in the next paragraph. Another input parameter needed in this case is the surface roughness of the asteroid. This roughness can also be constrained by the TPM if observations are made at several phase angles, but in our case it has to be assumed because we obtained only one spectrum. Little is known about the typical surface roughness of NEAs, and there is no consensus about what roughness parameters should be used (e.g., Harris et al. 2007 and references therein). We explored the effect of surface roughness by incorporating two values of surface roughness into our models (see below). The difference in the derived thermal inertia between the smooth surface case and the rough case was less than 30%. Because this uncertainty is considerably smaller than that due to the spin-pole orientation, it is not meaningful to try to refine the roughness parameters at this time.
As mentioned, the spin-pole orientation is a crucial input parameter for constraining the thermal inertia. However, the estimate of the spin-pole position for asteroid 162173 is uncertain because the amplitude of the rotational light curve is small and the range of geometries for the light curve observations has been limited. We, therefore, consider two cases for the orientation of the spin-pole. The first is the extreme case of an equatorial retrograde geometry, i.e., the spin-pole pointing at the South ecliptic pole. As input parameters for this TPM, we used a spherical shape (a reasonable assumption given the small amplitude of the light curve), the equatorial retrograde geometry, the rotation period of 0.3178 days, and a perfectly smooth surface (this last parameter yields a lower thermal inertia than in the rough surface case). Then we allowed the unconstrained parameters in the model (D, ,
and
)
to vary, an approach that worked particularly well because our spectrum spans a broad wavelength range that includes the peak in thermal emission. We generated a grid of models with
ranging from 0 to 2500 Jm-2 s-0.5 K-1, and the best-fit model was determined by minimizing the reduced
between the model fluxes and our data. This case (equatorial retrograde geometry and smooth surface) yields a rigorous lower limit to the thermal inertia of 150 Jm-2 s-0.5 K-1.
![]() |
Figure 2:
The the reduced |
For the second case, we adopted the pole orientation of Abe et al. (2008a) given in Sect. 1, which implies a subsolar latitude of approximately 30
at the time of the Spitzer observations. As in the first case, we generated a grid of models, and the best-fit model was determined by minimizing the reduced
.
In Fig. 1, we plot our best-fit TPM model along with the observed spectrum, and the corresponding parameters are listed in Table 1, again assuming
and the
and G from Hasegawa et al. (2008). Figure 2 shows the reduced
vs. thermal inertia and the estimated diameter and geometric albedo vs. thermal inertia for the rough-surface case (the parameters for the rough-surface model are
[half-angle opening of spherical section craters] and f = 0.77 [fractional surface coverage of craters]). The smooth surface model is not plotted because of space limitations, but fits equally well and yields about 30% lower values for the thermal inertia.
Uncertainties in thermal modeling usually exceed the formal errors resulting from the scatter of the flux measurements. The main sources of uncertainties in our derived physical parameters are the absolute photometric uncertainty in the thermal flux (diameter), the relative spectral uncertainty (temperature, beaming parameter and thermal conductivity), and the uncertainty in
(albedo), including the rotational light curve. Comparison of diameters derived from thermal models with those derived from other sources, such as radar, indicates that overall errors are normally less than 15% in diameter (Delbó et al. 2003). The contribution of relative spectral uncertainty to errors in the thermal continuum and derived parameters is dominated by the absolute-flux-calibration adjustment of the LL and SL segments of the spectrum (Sect. 2). The estimated uncertainties are given in Table 1, and in the case of the (second case) TPM, they are also illustrated in Fig. 2.
Spectral features superposed on the mid-infrared thermal continuum have been detected in other asteroids and can be compositionally diagnostic. However, except for a possible weak feature near 15 m (Fig. 1), no obvious structure is apparent in our spectrum, and a detailed analysis is deferred to a later publication.
4 Results and conclusion
Our results are summarized in Table 1, and the most relevant to the proposed sample return missions are the constraints our observations impose on the thermal inertia, .
When we adopt the extreme case of an equatorial retrograde geometry and a smooth surface for this asteroid, we derive a rigorous lower limit to the thermal inertia of 150 Jm-2 s-0.5 K-1. This value would be consistent with a fine regolith, similar to the one found in asteroid 433 Eros. Adopting the pole orientation of Abe et al. (2008a) yields a best-fit model with thermal inertia at 700 +/-200 Jm-2 s-0.5 K-1. Even higher values are allowed by the uncertainty in the spectral shape due to absolute flux calibration. This best-fit thermal inertia and its size place asteroid 162173 about two-sigma above the trend toward increasing thermal inertia with decreasing asteroid diameter as discussed by Delbó et al. (2007). The range in possible thermal inertias for asteroid 25143 Itokawa adopted by Delbó et al. (2007) (based on Müller et al. 2005 and Mueller 2007) is 350 to 800 Jm-2 s-0.5 K-1, which is within the range of possible values for asteroid 162173. If the spin-pole of Abe et al. (2008a) turns out to be correct, these last two results would suggest that the surface of asteroid 16217 has a texture rockier than or similar to that of asteroid 25143 Itokawa. However, given the uncertainties in the spin-pole position, we hesitate to even make this analogy.
Our estimates of the diameter and albedo of asteroid 162173 are consistent with those of Hasegawa et al. (2008). However, there are marginally significant discrepancies that may deserve further attention. There are temperature differences between their space-based (Akari) and ground-based (Subaru) data sets and between their data and ours. As a way to illustrate these differences we give our best-fit NEATM parameters for their Akari and Subaru data, which are
km,
,
and
km,
,
,
respectively. Only for
are the differences slightly larger than one sigma. It may be possible to constrain the spin-pole orientation better by modeling all three data sets (ours and the two sets in Hasegawa et al. 2008) simultaneously, allowing the pole solution to be a free parameter; however, this is work beyond the scope of this letter.
In summary, thermal models of our Spitzer spectrum of asteroid 162173 (1999 JU3) provide a rigorous lower limit to the thermal inertia of 150 Jm-2 s-0.5 K-1, which is unlikely but possible. This value would be consistent with a fine regolith, similar to the one found in asteroid 433 Eros. However, if we adopt the pole orientation of Abe et al. (2008), the best-fit model yields a thermal inertia of
Jm-2 s-0.5 K-1. Even higher values are allowed by the uncertainty in the spectral shape due to absolute flux calibration. Our best estimates of the diameter (
km) and geometric albedo (
)
of asteroid 162173 are consistent with those of Hasegawa et al. (2008).
Acknowledgements
H.C. acknowledges support from NASA's Spitzer Science Center, Jet Propulsion Laboratory, and Planetary Astronomy program. H.C. was a visiting Fulbright Scholar at the ``Instituto de Astrofísica de Canarias'' in Tenerife, Spain. J.L. acknowledges support from the Spanish ``Ministerio de Ciencia e Innovación'' projects AYA2005-07808-C03-02 and AYA2008-06202-C03-02. We also thank an anonymous referee for helpful comments.
References
- Abe, M., Kawakami, K., Hasegawa, S., et al. 2008a, 37th COSPAR Scientific Assembly, Montréal, Canada, 12 (In the text)
- Abe, M., Kawakami, K., Hasegawa, S., et al. 2008b, LPSC XXXIX, abstract #1594
- Binzel, R. P., Harris, A. W., Bus, S. J., et al. 2001, Icarus, 151, 139 [NASA ADS] [CrossRef] (In the text)
- Bottke, W. F., Morbidelli, A., Jedicke, R., et al. 2002, Icarus, 156, 399 [NASA ADS] [CrossRef] (In the text)
- Campins, H., Kelley, M. S., Fernández, Y., Licandro, J., & Hargrove, K. 2009, Earth, Moon & Planets, Aug. 2009, in press (In the text)
- Delbó, M. 2004. Ph. D. Thesis, Freie Universität, Berlin (In the text)
- Delbó, M., & Tanga, P. 2009, Planet. Space Sci., 57, 259 [NASA ADS] [CrossRef] (In the text)
- Delbó, M., Harris, A. W., Binzel, R. P., Pravec, P., & Davies, J. K. 2003. Icarus, 166, 116 (In the text)
- Delbó, M., dellOro, A., Harris, A. W., Mottola, S., & Mueller, M. 2007, Icarus, 190, 236 [NASA ADS] [CrossRef] (In the text)
- Emery, J. P. et al. 2009, in preparation
- Gaffey, M. J., Cloutis, E. A., Kelley, M. S., et al. 2002, in Asteroids III, ed. W. F. Bottke Jr., A. Cellino, P. Paolicchi, & R. P. Binzel (Tucson: Univ. of Arizona Press), 183 (In the text)
- Harris. A. W. 1998, Icarus, 131, 291 [NASA ADS] [CrossRef] (In the text)
- Harris, A. W., & Lagerros, J. S. V. 2002, in Asteroids III, ed. W. F. Bottke Jr., A. Cellino, P. Paolicchi, & R. P. Binzel (Tucson: Univ. of Arizona Press), 205 (In the text)
- Harris, A. W., Mueller, M., Delbó, M., et al. 2007, Icarus, 188, 414 [NASA ADS] [CrossRef] (In the text)
- Hasegawa, S., Müller, T. G., Kawakami, K., et al. 2008, Pub. Astron. Soc. Japan, 60, S399 [NASA ADS] (In the text)
- Houck, J. R., Roellig, T. L., van Cleve, J., et al. 2004, Astrophys. J. Suppl. Ser. 154, 18 (In the text)
- Kelley, M. S., Wooden, D. H., Tubiana, C., et al. 2009, AJ, 137, 4633 [NASA ADS] [CrossRef] (In the text)
- Lebofsky, L. A., & Spencer, J. R. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews (Tucson: Univ. of Arizona Press), 128 (In the text)
- Lebofsky, L. A., Sykes, M. V., Tedesco, E. F., et al. 1986, Icarus, 68, 239 [NASA ADS] [CrossRef] (In the text)
- Mueller, M. 2007. Ph.D. Thesis, Freie Universität, Berlin (In the text)
- Müller, T. G., Sekiguchi, T., Kaasalainen, M., Abe, M., & Hasegawa, S. 2005, A&A, 443, 347 [NASA ADS] [CrossRef] [EDP Sciences] (In the text)
- Rivkin, A. S., Howell, E. S., Vilas, F., et al. 2002, in Asteroids III, ed. W. F. Bottke Jr., A. Cellino, P. Paolicchi, & R. P. Binzel (Tucson: Univ. of Arizona Press), 235 (In the text)
- Tholen, D. J., & Barucci, M. A. 1989, in Asteroids II, ed. R. P. Binzel, T. Gehrels, & M. S. Matthews (Tucson: Univ. of Arizona Press), 298 (In the text)
- Vilas, F. 2008, AJ, 135, 1101 [NASA ADS] [CrossRef] (In the text)
- Werner, M. W., Roellig, T. L., Low, F. J., et al. 2004, ApJS, 154, 1 [NASA ADS] [CrossRef] (In the text)
Footnotes
- ... unusual
- We considered and ruled out other potential sources for the SL and LL flux difference, such as placement of the object within the slit and rotational variability of the asteroid.
All Tables
Table 1: NEATM and TPM best-fit models.
All Figures
![]() |
Figure 1: The Spitzer spectrum of near-Earth asteroid 162173 (1999 JU3) plotted along with our best-fit thermophysical model (TPM). The model parameters are given in Table 1 and discussed in Sect. 3. The LL and SL fluxes have been scaled to match each other in the overlapping wavelengths (see Sect. 2). |
In the text |
![]() |
Figure 2:
The the reduced |
In the text |
Copyright ESO 2009
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.