arxiv: v1 [astro-ph.ep] 10 Nov 2016

Transcription

1 Near-infrared thermal emission from near-earth asteroids: Aspect-dependent variability arxiv: v1 [astro-ph.ep] 10 Nov 2016 Nicholas A. Moskovitz a,c,1 David Polishook b Francesca E. DeMeo c Richard P. Binzel c Thomas Endicott d Bin Yang e,f Ellen S. Howell g Ronald J. Vervack, Jr. h Yanga R. Fernández i a Lowell Observatory, 1400 West Mars Hill Road, Flagstaff, AZ (U.S.A.) b Department of Earth and Planetary Sciences, Weizmann Institute of Science, Rehovot (Israel) c Massachusetts Institute of Technology, Department of Earth, Atmospheric and Planetary Sciences, 77 Massachusetts Avenue, Cambridge, MA (U.S.A.) d University of Massachusetts, Boston, 100 Morrissey Blvd., Boston, MA (U.S.A.) e Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI (U.S.A.) f European Southern Observatory, Alonso de Cordova, 3107, Vitacura, Casilla 19001, Santiago de Chile (Chile) g Lunar and Planetary Lab, University of Arizona, Tucson AZ h JHU/Applied Physics Lab, Johns Hopkins Road, Laurel MD (U.S.A.) i University of Central Florida, Dept. of Physics, 4000 Central Florida Blvd. Preprint submitted to Elsevier 3 April 2018

2 Orlando FL (U.S.A.) Copyright c 2016 Nicholas A. Moskovitz Number of pages: 28 Number of tables: 2 Number of figures: 5 1 Observations conducted while at MIT, current affiliation is Lowell Observatory. 2

3 Proposed Running Head: NEA Thermal Emission Please send Editorial Correspondence to: Nicholas A. Moskovitz Lowell Observatory 1400 West Mars Hill Road Flagstaff, AZ 86001, USA. Phone: (928)

4 1 Introduction Obliquity (a scalar) or spin vector is a fundamental property of all celestial bodies. In the Solar System, spin vectors directly influence the efficiency in which the Yarkovsky and YORP effects (i.e. thermal radiation forces and torques) change the orbital and rotational properties of minor planets (Rubincam, 2000; Bottke et al., 2006; Vokrouhlicky et al, 2015). The Yarkovsky effect plays a key role in the dynamical transport of small asteroids from the Main Belt into the near-earth population (Farinella & Vokrouhlický, 1999). The YORP effect directly influences rotation state and is frequently invoked to explain the formation of binary and multiple asteroid systems through rotational fission or mass shedding (Margot et al., 2002; Scheeres, 2007; Walsh et al., 2008; Pravec et al., 2010; Jacobson & Scheeres, 2011). Furthermore, spin states are a direct consequence of asteroid-asteroid impacts and thus trace collisional environments across dynamical populations (Paolicchi et al, 2002). Despite the importance of spin vectors, only a few hundred out of three quarters of a million known asteroids have constrained pole orientations. This is a direct consequence of the many years or even decades required to constrain spin vectors through astrometric or photometric observations (e.g. Kaasalainen et al., 2002; Slivan, 2002; Lowry et al., 2007; Nugent et al., 2012). Analyses across shorter temporal baselines often suffer from retrograde-prograde degeneracies. Doppler delay radar imaging can determine shape and pole solutions, but requires sufficient sky motion during the observing window and is limited to the closest near-earth and largest Main Belt asteroids (Ostro et al, 2002). Here we focus on an alternate technique for constraining spin orientation, namely by monitoring thermal emission as a function of solar phase angle. 4

5 This technique does not solve for the obliquity or specific spin vector, but can constrain retrograde versus prograde spin orientation. Measuring thermal emission from a body at multiple viewing geometries on either side of opposition can reveal differences in local morning versus afternoon temperature distributions across the surface of an asteroid (Figure 1). The rotational sense of the body directly influences the local time and temperature distribution on the asteroid surface as perceived from Earth, particularly when observing at geometries of high phase angle. This afternoon-morning temperature dichotomy can observationally manifest as a difference in thermal emission. For the case depicted in Figure 1, preferential viewing of the cold hemisphere before opposition and the warm hemisphere afterwards would be reversed if this object were in retrograde rotation. Therefore, variations in thermal flux on either side of opposition can serve to constrain pole orientation. Thermal observations across a range of phase angles is a standard method for determining asteroid thermal inertia (Delbo et al, 2015). This technique, originally described (Matson, 1971) and pioneered in the 1970 s, has not been widely used to constrain pole orientation, but has seen some success. Morrison (1976) confirmed that the near-earth asteroid 433 Eros is in a prograde rotation state based on enhanced 10- and 20-µm thermal emission prior to opposition, consistent with preferential viewing of local afternoon. Even though Eros pole only points 11 above the ecliptic (Zuber et al., 2000) its moderate orbital inclination (10.8 ) and observations by Morrison (1976) at moderate to high solar phase angles (-41, -24, and +42 ) ensured that the afternoon and morning sides of Eros were preferentially viewed prior to and after opposition respectively. Morrison (1977) and Hansen (1977) extended this technique to several large Main Belt asteroids. They were able to correctly predict prograde rotation for 1 Ceres (Johnson et al., 1983), 4 5

6 June 15! = -39 July 5! = -36 June 8! = -30 June 2! = +18 May 30! = +34 Afternoon May 11! = +71 Morning Fig. 1. Schematic of the opposition ( ) centered orbital longitude of 1998 QE2 during its 2013 apparition. The dates and corresponding solar phase angles (α, signed relative to opposition) are indicated for each of our observations. The asteroid cartoons are shaded to indicate illumination (black = nighttime hemisphere) and surface temperate distribution for a body with finite thermal inertia (blue = cold, red = hot). The predicted prograde rotation of the asteroid (Springmann et al., 2014) is indicated by the black arrows. The cartoons of the asteroid at maximal phase angles (+71 and 39 ) are larger to emphasize the preferential viewing of cooler morning temperatures prior to opposition and warmer afternoon temperatures post-opposition. Note the phase angle actually decreased after the June 15 observation. Earth image credit: NASA/GSFC. Vesta (Gehrels, 1967) and 19 Fortuna (Torppa et al., 2003), and a retrograde state for 10 Hygeia (Kaasalainen et al., 2002). However they did not correctly predict the retrograde rotation of 2 Pallas (Torppa et al., 2003; Carry et al., 2010). There are likely several reasons for this. First, the prediction for Pallas was tenuous anyways due to it having the fewest observations of the objects under investigation. Second, the standard thermal model used to derive physical properties (including pole orientation) assumed non-rotating spherical bodies. Third, the range of accessible phase angles for objects in the Main Belt is limited relative to near-earth objects. Finally, the possibility of low surface thermal inertia for large (> 100 km) asteroids (e.g. Delbo et al, 6

7 2015) acts to dampen morning-afternoon temperature gradients. In light of these issues, it is impressive that these early investigations resulted in meaningful constraints on the rotation states of any of the targeted asteroids. Since these initial studies, the use of a morning-afternoon temperature dichotomy to inform pole orientation has only been discussed and applied in a few cases (e.g. Lebofsky et al., 1986; Harris et al., 2007; Muller et al., 2012). Note that this method is not sensitive to objects with poles aligned near to the line of sight, nor does it provide a measure of the object obliquity, instead it simply provides a means to constrain a prograde versus retrograde spin state. Traditionally, the measurement of asteroid thermal emission and subsequent modeling of physical properties has relied on mid-infrared (3-25 µm) photometry. An asteroid thermo-physical model called SHERMAN has been developed for application to airless bodies of arbitrary shape and has primarily been applied to (though is not restricted to) observations at wavelengths between µm (Howell et al., 2015). However, low albedo (< 10%) asteroids at small heliocentric distances (< 1.2 AU) have high enough surface temperatures that the Wien tail of their thermal emission profiles can be detected at near-infrared wavelengths ( µm). At such short wavelengths the measured signal is a combination of solar reflectance and thermal emission. This thermal excess or thermal tail can be modeled to retrieve surface temperature and albedo (Rivkin et al., 2005; Reddy et al., 2009, 2012). However, models of this excess from normalized reflectance spectra are not directly sensitive to diameter because such data are not flux calibrated. Here we present case studies of the near-earth asteroids (285263) 1998 QE2 ( 3) and (175706) 1996 FG3 ( 4) based on observations during apparitions favorable to the detection of phase-dependent thermal variability. This analysis 7

8 represents the first systematic study of thermal variability in the near-ir K- band as a function of phase angle. 2 Observations and Thermal Modeling Both 1998 QE2 and 1996 FG3 have low geometric albedos of about 3% and 4% respectively (Trilling et al., 2010; Wolters et al., 2011), and thus display significant thermal emission in the near-infrared when they are at heliocentric distances of 1 AU. We present near-infrared spectra of 1998 QE2 taken between May and July of 2013 and of 1996 FG3 taken during apparitions in 2009 and 2011 (Tables 1 and 2). All data analyzed here were obtained with SpeX at NASA s Infrared Telescope Facility, which was configured in its low resolution (R 100) prism mode with a 0.8 slit (Rayner et al., 2003). All observations were conducted with the telescope operating in a standard ABBA nod pattern. Commonly-used solar analogs (DeMeo et al., 2009) were observed to correct for telluric absorption and to remove the solar contribution from the measured reflectance. Reduction protocols followed Moskovitz et al. (2013) and DeMeo et al. (2009). 8

9 Observation Date V r α η 3% η 6% η 1% (UT) [mag] [AU] [AU] [deg] 2013 May ± May ± June ± June ± June ± July ± Table 1 Summary of (285263) 1998 QE2 observations and best-fit NEATM thermal beaming parameter (η) values for the full range of possible albedos from 1-6%. 9

10 Observation Date V r α η3.9% η5.1% η2.7% (UT) [mag] [AU] [AU] [deg] Mutual Event? 2009 March ± No 2009 April ± Secondary occulted by primary 2011 December ± Secondary occulted by primary 2011 December ± Primary eclipsed and occulted by secondary Table 2 Summary of (175706) 1996 FG3 observations and best-fit NEATM η values for the full range of possible albedos from %. 10

11 We model the thermal tails of our targets with the near-earth asteroid thermal model (NEATM; Harris, 1998). A number of parameters go into this model including slope parameter (G), visible to near-ir reflectance ratio (0.55 µm / 2.5 µm), emissivity, and the coupled parameters geometric albedo (p V ), diameter, and absolute magnitude (H). The NEATM assumes that the asteroid is spherical and that the hottest point on the surface is the sub-solar point (i.e. the surface is in instantaneous equilibrium with the incident solar flux). The NEATM thermal beaming parameter η compensates for a lack of explicit treatment of thermophysical properties like thermal inertia and surface roughness. In the NEATM, η is inversely proportional to temperature to the fourth power 2 and thus can indicate whether the measured thermal emission is higher or lower than expected from an idealized, smooth surface with zero thermal inertia at zero phase angle. Previous application of the NEATM to multi-epoch observations has successfully produced results that are broadly consistent with those from more sophisticated thermophysical models (e.g. Mueller et al., 2006; Harris et al., 2007; Delbo et al, 2015). Delbo (2004) presented a comprehensive suite of NEATM models that quantitatively describe thermal variability as a function of solar phase angle as well as other physical parameters. These models suggest several possible caveats to the detectability of aspect-dependent thermal variability. (1) The dimensionless thermal parameter Θ (Spencer et al., 1989) quantifies the ratio of radiative to rotational timescales, and is a function of spin period, thermal inertia (Γ), and subsolar temperature as dictated by heliocentric distance. Though the smallest near-earth asteroids can spin rapidly, the vast majority of NEAs with sizes bigger than about Technically, in NEATM formalism ηt 4 SS is a constant, where T SS is the sub-solar equilibrium temperature on the surface of the asteroid. 11

12 meters have rotation periods on the order of several hours (e.g. Polishook et al., 2012). Thus, for typical NEA rotation periods and heliocentric distances ( 1 AU), a thermal inertia greater than about 200 in MKS units (J m 2 s 0.5 K 1 ) ensures that heat deposited by solar insolation does not radiate away so quickly as to prevent morning-afternoon temperature differences. The typical thermal inertia for asteroids with diameters < 5km would be in that range (Delbo et al., 2007). While high thermal inertia certainly facilitates a morning-afternoon temperature dichotomy, thermally-constrained poles for several large Main Belt asteroids (Morrison, 1977; Hansen, 1977) suggests that either these early determinations were simply lucky or that high thermal inertia is not an exclusive requirement. (2) The asteroid must be observable across a wide range of solar phase angles, particularly at phases in excess of 50. In general, such phase angle ranges are only accessible for NEAs. (3) If thermal emission is to be detected in the near-ir K-band, then the object needs to be at a heliocentric distance inside of about 1 AU and have an albedo < 10%. (4) The asteroid should be observed during a single apparition to avoid significant changes in viewing aspect that could occur across multi-epoch observations. Even though all of these properties were not necessarily known in advance for QE2 and FG3, their low albedos and wide range of observable phase angles made them good candidates for monitoring thermal variability in the near-ir as a function of viewing aspect. 12

13 3 Case Study: (285263) 1998 QE2 Radar observations of 1998 QE2 in 2013 revealed a diameter of approximately 3 km and a prograde rotation state with the pole pointing approximately 40 above the ecliptic 3 (Springmann et al., 2014). From the JPL Horizons system the absolute magnitude of QE2 is H = 17.3, which is revised from a value of H = 16.4 used in previous works (e.g. Trilling et al., 2010). An updated analysis of warm Spitzer observations using this new H magnitude yields a low albedo for QE2 of % (Trilling et al., 2016). We adopt H = 17.3, even though it may be less accurate than the value of H = measured by Hicks et al. (2013), because it is fully consistent with the formal log-normal albedo error bars as computed by Trilling et al. (2016). As we discuss below, these are important for investigating systematic errors in our NEATM models due to uncertainties on the input parrameters. Figure 1 depicts a schematic of the opposition-centered orbital geometry of QE2 during its 2013 apparition. To facilitate our analysis we employ a signed phase angle (α) based on the geocentric ecliptic longitude of asteroid opposition, specifically α < 0 when the asteroid s ecliptic longitude minus the Sun s ecliptic longitude is greater than 180. The putative prograde rotation of QE2 suggests that prior to opposition (α > 0 ) Earth-based observers would preferentially see the cooler morning side of the asteroid, whereas post-opposition (α < 0 ) observations would preferentially access the warmer afternoon side. Near-IR spectra of QE2 spanning opposition are shown in Figure 2. A composite visible/near-infrared spectrum of QE2 is shown at the bottom 3 Incidentally these observations also discovered a small satellite roughly 750 meters in size. Unless the satellite has a dramatically different albedo or thermal properties relative to the primary, its contribution to the system s measured flux should not significantly influence any of the results presented here. 13

14 4 May 11 May 30 Normalized Reflectance 3 2 June 2 June 8 June 15 1 July 5 Nogoya (CM2 Chondrite) Wavelength (microns) Fig. 2. IRTF/SpeX spectra of 1998 QE2 taken in the months surrounding its opposition passage in mid The asteroid spectrum from July 5 also includes visible data ( µm) from Hicks et al. (2013). QE2 is taxonomically classified as a Ch type (DeMeo et al., 2009) and is a close spectral analog to the CM2 carbonaceous chondrite Nagoya. The measured spectra are shown in grey, the NEATM-corrected spectra with thermal tails removed are overplotted in black. Spectra have been vertically offset for clarity. of Figure 2. Due to the presence of a weak 0.7 µm absorption feature the asteroid is best fit by a Ch-type taxonomy in the DeMeo et al. (2009) system. Using the technique outlined in Moskovitz et al. (2013) we find that QE2 has a 14

15 spectrum and albedo most closely matched to RELAB (Pieters, 1983) spectra of carbonaceous chondrites. The CM2 Nagoya is a particularly good spectral analog (Fig. 2). A potential link between Ch-asteroids and CM chondrites has been noted previously (e.g. Burbine, 1998). This taxonomic classification is consistent with the low Spitzer albedo. The spectra in Figure 2 clearly show thermal emission long-wards of 2µm. We model these thermal tails with the following NEATM parameters: geometric albedo p V = % (Trilling et al., 2016); slope parameter G = 0.086, equal to the mean for C-complex asteroids (Warner et al., 2009); visible to near- IR reflectance ratio = 1.11, measured from our combined visible and near-ir spectral data (Fig. 2); absolute magnitude H = 17.3 from the JPL Horizons system; and emissivity ɛ = 0.9. These values are either consistent with observations of QE2 or are typical for objects like QE2, and are held constant for all of the modeling conducted here. Again we adopt H from JPL as opposed to Hicks et al. (2013) so that we can propogate the albedo error bars from Trilling et al. (2016) who also used the JPL value. With each of these parameters held fixed, the beaming parameter η is left as a free parameter. We vary η to best remove the thermal tails from the spectra in Figure 2. The best fit η for each spectrum is determined by minimizing the error-weighted RMS residual between the thermally corrected spectrum and a linear fit to the data between 1.9 and 2.49 µm. This approach thus requires three free parameters: η and the two parameters associated with the linear fit. Uncertainties on the fitted η values reflecting the signal-to-noise of the data (i.e. statistical errors) are conservatively determined by bracketing the range of η values that clearly over- or under-compensate the thermal emission. These uncertainties are somewhat subjective but do correspond to deviations from 15

16 the best fit by approximately one-sigma as dictated by the noise level in the spectral data. Best-fit η values and signal-to-noise related uncertainties are given in Table 1 for the nominal assumed albedo of 3%. A similar approach to modeling of near-ir spectra and estimating uncertainties on the model parameters has been employed elsewhere (e.g. Rivkin et al., 2005; Reddy et al., 2012). We further investigate uncertainty on η by considering systematic errors on the NEATM input parameters which would produce an overall offset to the best-fit η values. The assumed albedo is the dominant source of systematic uncertainty in η, which for QE2 ranges from 1-6% (Trilling et al., 2016). We fit η assuming these extreme values and include those additional fits in Table 1. The uncertain albedo for QE2 results in large systematic errors ( 50%) in the absolute η values. However, any change in the albedo from the nominal 3% affects all η values in the same way (increase or decrease) and by approximately the same amount. The assumed value of G = also represents a potential source of systematic error. However, considering a wide range of G values from -0.1 to results in η fits that span less than half the range produced by considering the full extent of possible albedos. In other words the systematic errors induced by uncertain G are approximately 20% on the best fit η values. The NEATM results for QE2 are shown in Figure 3. An asymmetric trajectory in this η α parameter space can indicate physical properties such as sense of rotation or thermal inertia, however the use of NEATM to diagnose these physical traits must be treated with caution due to known systematic errors inherent to NEATM (Wolters & Green, 2009) and known correlations between η and α (Delbo et al., 2003; Wolters et al., 2008; Masiero et al., 2011). Nevertheless, the results in Figure 3 are suggestive that a small number of carefully 16

17 Γ= May 11 NEATM η June 15 July 5 June 8 June 2 May 30 Γ= PM Phase Angle (deg) AM Fig. 3. NEATM η as a function of solar phase angle, α, for 1998 QE2 during its 2013 apparition. As QE2 moved through opposition it tracks from the upper right corner to the lower left of this η α space (black circles). NEATM models for asteroids with prograde rotation, surface roughness θ = 36, and thermal inertia (Γ) of 200 and 400 in MKS units are over-plotted as the solid and dashed grey lines respectively (Delbo, 2004). For a prograde rotator the predicted η α relationship is asymmetric with higher η (lower morning temperatures) at positive phase angles and lower η (higher afternoon temperatures) at negative phase angles. The curves for retrograde rotators would be flipped about α = 0. Statistical errors based on the signal-to-noise of our data are in black, while systematic error bars accounting for uncertainty in albedo are shown in grey. timed observations in the near-infrared, particularly at phase angles > 50, could be used to constrain physical properties of small NEAs. Unfortunately the lack of repeat observations of QE2 at high phase angles means that the results for this specific object are tenuously hinged to the single May 11 data point. In consideration of this technique, we compare our QE2 results to model curves extracted from Delbo (2004). Specifically, we compare our η values to a suite 17

18 of models that correspond to bodies with a macroscopic surface roughness parameter θ = 36. This particular value was selected as an intermediate case in the range of models presented by Delbo (2004) and is consistent with expected roughness parameters from more sophisticated photometric models of other well studied asteroids (e.g. Delbo & Tanga, 2009). The assumption of this particular roughness does not significantly influence the results presented here. If we were to adopt other roughness values (10, 20, 58 ) then the curves in Figure 3 would shift by no more than 10 20%. More specifically, the afternoon-side curves (α < 0 ) remain largely unchanged across all roughness values, whereas the morning-side (α > 0 ) curves either get steeper by 20% for the largest roughness or shallower by 10% for the smallest roughness value. As seen in Figure 3, the large η prior to opposition (May 11) suggests a cooler relative surface temperature, consistent with (though not uniquely diagnostic of) viewing of local morning on the asteroid surface and hence prograde rotation. Furthermore, the η values and magnitude of the η-enhancement at high positive phase angles are consistent with models for a prograde rotator with thermal inertial Γ (Delbo, 2004). Interestingly, a thermal inertia around 200 is the expectation for a 3-km asteroid like QE2 (Delbo et al, 2015). Again, the assumption of surface roughness = 36 does not significantly affect the rough fit to our data by theoretical curves corresponding to thermal inertia values around For the full range of roughness values considered by Delbo (2004), a thermal inertia of 50 is too low to produce large variations in η, and thermal inertia of 900 results in values of η that are much larger than those retrieved from our data. The systematic uncertainty in our η values due to the range of possible albedos could shift the points in Figure 3 up or down by several tens of percent, however the general result that the 18

19 data are indicative of a thermal inertia of a few hundred would still hold. Furthermore, the relative orientation and increase at high positive phase angle of the η α points in Figure 3 would be preserved. 4 Case Study: (175706) 1996 FG3 The NEA 1996 FG3 is also a binary system (Pravec & Harris, 2000) with a primary about 1.7 km in size (Wolters et al., 2011; Mueller et al., 2011; Walsh et al., 2012; Yu et al., 2014), a spacecraft accessible orbit (i.e. v = 6.6 km/s), and a putative primitive composition (Binzel et al., 2001; de Leon et al., 2011; Walsh et al., 2012; de Leon et al., 2013). FG3 s satellite is roughly 450 meters in size (Yu et al., 2014) and is in a tidaly-locked retrograde orbit (Benner et al., 2012; Scheirich et al., 2015), which suggests that the primary is also in a retrograde spin state. Again, the satellite is unlikely to have a significant influence on the system s thermal flux unless it has a dramatically different albedo or thermal properties relative to the primary. A low albedo around 4% (Wolters et al., 2011) makes FG3 another ideal object for monitoring thermal variability in the near-ir. As with QE2, the flux contribution from the small satellite should not have a significant effect on any of the measurements presented here. Numerous near-ir spectra of FG3 have been obtained during apparitions in 2009 and 2011 (de Leon et al., 2011; Walsh et al., 2012; de Leon et al., 2013). We analyze four individual spectra from these apparitions (Table 2; Figure 4). These spectra can be qualitatively divided into two groups: in 2009 the spectra have a neutral spectral slope and minimal thermal excess, whereas in 2011 the spectra show a pronounced red spectral slope and clear thermal emission. 19

20 Additional near-ir spectra not analyzed here (de Leon et al., 2011; Rivkin et al., 2012) are consistent with this qualitative description. Taxonomically (DeMeo et al., 2009), the 2009 data suggest that FG3 is a member of the C-complex while the 2011 data suggest membership in the X-complex. The repeatability of these spectra on different nights in each of the 2009 and 2011 epochs suggests that this variability is real and not an observational artifact. We have compared the timing of these observations to predictions for mutual events (P. Scheirich, private communication). Those predictions are included in Table 2. In short, the occurrence of mutual events does not seem to have any major influence on the spectral data, which is not surprising considering the significant size difference between the primary and secondary. The two spectra that differ the most, i.e. the 2009 April 27 and 2011 December 1 data, were both obtained when the secondary was occulted by the primary. Furthermore, within each epoch the pairs of spectra were obtained during different system configurations, again suggesting that the secondary had little optical and thermal influence on the measurements. It is worth noting however that occultation of the primary by the secondary, particular if the secondary shadows the hottest region on the primary s surface, could have significant thermal implications. This occultation would lead to reduced thermal emission and a correspondingly elevated η. The only data analyzed here that corresponded to occultation of the primary were from 2011 December 30. This data point is further discussed below. As with 1998 QE2 ( 3) we model the thermal tails with the NEATM. The following parameters are employed: geometric albedo p V = 3.9 ± 1.2% (Walsh et al., 2012); absolute magnitude H = ± 0.03 and slope parameter G = 0.07±0.02 as measured by Pravec et al. (2006); visible to near-ir reflectance 20

21 Normalized Reflectance Wavelength (microns) Fig. 4. Reflectance spectra of 1996 FG3 (grey lines), offset for clarity. In general the data from 2009 display a neutral spectral slope and limited thermal emission, whereas the data from 2011 show a red spectral slope and pronounced thermal emission. Visible data from Binzel et al. (2001) have been appended to the bottom spectrum. NEATM corrected spectra are over-plotted in black. ratio ranging from as measured from our combined visible and near-ir spectra; and emissivity ɛ = 0.9. These values are held constant for all of the FG3 models and again η is left as a free parameter in the NEATM fits. Our NEATM results for FG3 are shown in Figure 5. As with QE2 we compute error bars on the derived η values that account for uncertainty due to the signal-tonoise of the data and we compute systematic offsets in η attributable to the range of possible albedos for FG3 (Walsh et al., 2012). Again the uncertainty on albedo is the dominant source of error for these models. Previous results based on observations at longer wavelengths are also included in Figure 5; in particular, NEATM η values from 2 4 µm spectral observations by Rivkin et al. (2013) and thermal infrared ( 8 20 µm) photometry by Walsh et al. 21

22 (2012). NEATM η Ecliptic lat. (deg) Dec 01 AM Γ=400 Γ=200 Dec 06 Mar 30 Dec 24 Dec 30 Apr 27 May Phase Angle (deg) PM Fig. 5. Top panel: The η α relationship for 1996 FG3 based on 2009 (black) and 2011 (grey) observations. Filled circles are the data analyzed here (Figure 4). Statistical errors based on the signal-to-noise of the data are shown in black, systematic errors due to the range of possible albedos are in grey. Open circles are from Rivkin et al. (2013) and the open triangle is from Walsh et al. (2012). NEATM models for a retrograde rotator with thermal inertia (Γ) of 200 and 400 in MKS units, and surface roughness θ = 36 are over-plotted (Delbo, 2004). The curves for prograde rotators would be flipped about α = 0. The 2011 data are broadly consistent with the models for retrograde rotation, however the 2009 data are not. Bottom panel: the geocentric ecliptic latitude of the asteroid at the time of each observation. The discrepancy between the data and the models could be explained by a texturally and spectrally heterogeneous surface revealed by the different viewing geometries, i.e. preferential viewing of the southern hemisphere in 2011 and of the northern hemisphere in Unlike the η α results for 1998 QE2, the fits in Figure 5 do not clearly follow the NEATM predictions for a retrograde rotator. In fact the highest values of η for FG3 correspond to viewing geometries that would have preferentially viewed the afternoon hemisphere for a retrograde rotator. Our data thus suggest that FG3 is in a prograde rotation state. This apparent contradiction is 22

23 discussed in the following section. Values of Γ most closely overlap our data points and are broadly consistent, though higher by a factor of 2, with other estimates for FG3 s surface (Wolters et al., 2011). Interestingly the η value for December 30 is slightly higher than expected for a low phase angle observation. As mentioned above this could be attributed to observation of the system while the primary was occulted by the secondary. Modeling of this scenario is beyond the scope of this work, and given the unexplained variability in η for this object we do not pursue this idea any further. 5 Discussion We have demonstrated the detection of phase dependent variations in thermal beaming parameter for the NEAs (285263) 1998 QE2 and (175706) 1996 FG3. While such variability is not novel correlations between phase angle and η have been shown previously (e.g Delbo et al., 2003; Wolters et al., 2008; Masiero et al., 2011) this work represents the first systematic analysis of such variability at near-ir wavelengths. Ultimately the goal of any such investigation would be to correlate actual physical properties with derived model parameters. It is in this vein that we discuss the thermal variability seen for QE2 and FG3. As discussed above, the NEATM modeling of the QE2 data indicate a significant increase in η at high positive phase angles (Figure 3). This is reasonably well fit with models for a prograde pole orientation and a surface thermal inertia of roughly J m 2 s 0.5 K 1, properties which are consistent with other observations (Springmann et al., 2014) and expectations for objects of this size (Delbo et al., 2007). It is unclear whether the pole, pointing 23

24 at ecliptic latitude +40 is near enough to purely prograde to fully explain the observed thermal effects as a function of solar phase angle. While our QE2 results are suggestive that the thermal variability in the near- IR is indicative of physical properties, a data point at large negative phase angles (α < 50 ) would have helped to definitively determine whether the apparent asymmetry in this η α space is a signature of pole orientation or simply a side effect of the NEATM method. The diagnostic importance of such high phase angle measurements on both sides of opposition is an important lesson for guiding future observational campaigns. In general, the models of Delbo (2004) suggest that objects with Γ << 200 would not display significant differences in η as a function of phase angle, however future targeting of low albedo NEAs at near-ir wavelengths across a range of sizes are needed to test this model prediction. The interpretation of the FG3 data is less clear, which is perhaps unsurprising given the numerous previous attempts to understand its variability in reflectance and thermal emission (de Leon et al., 2011; Binzel et al., 2012; de Leon et al., 2013). We have built upon these previous works and show that there are clear differences in η as a function of phase angle (Figure 5), and that those differences are asymmetric about opposition. We argue that these differences are not an artifact of any observational issues such as signal-tonoise or data calibration. The 2009 data are of lower signal-to-noise relative to the 2011 data and they do show the greatest offsets from model predictions (Fig. 5), but there are no clear indicators of poor calibration. For example, the near complete removal of any residual telluric absorption features around 2.0 µm suggests that the dispersion solution is good and that the atmosphere was stable enough during those nights to enable removal of the telluric bands. Fur- 24

25 thermore, there are no odd continuum or slope changes in the 2009 data that might indicate use of a bad solar analog. Certainly, the repeatability of spectral properties in each of the two epochs suggests that the spectral changes are real. Furthermore, the results of Walsh et al. (2012) and Rivkin et al. (2013) are in broad agreement with our NEATM η values (Fig. 5), thus providing further support for a scenario in which spectral and thermal changes did occur. However, these detected changes, particularly the η α relationship, are not consistent with an expected retrograde spin state (Scheirich et al., 2015; Yu et al., 2014). As such, we offer three possible explanations for why our results may not be consistent with retrograde models, however in all three cases additional data are needed to better understand the unusual spectral and thermal behavior displayed by FG3. First, the variability in the data could be a consequence of the known correlation between η and phase angle in NEATM models (Spencer et al., 1989; Wolters & Green, 2009; Reddy et al., 2012). Values of η derived from NEATM reflect the observed surface temperature. At high phase angles surface roughness causes increased shadowing and an apparent drop in surface temperature (i.e. higher η). Such correlations are typically represented as linear fits to an ensemble of data. Table 3 presents a compilation of coefficients based on linear fits to various NEA data sets. Also included are fit parameters to the FG3 data from the 2009 and 2011 apparitions with errors indicating 1σ uncertainties. These linear fits are error weighted by the systematic uncertainties on each of the η values. In comparison to the ensemble fits our FG3 values are all shifted to higher η, i.e. larger intercepts, and shallower slopes. The η values from the QE2 data are also slightly higher than the mean values from the studies represented in Table 3. Though it is not entirely correct to compare our two objects to the diverse ensemble of properties (surface roughness, thermal inertia, mor- 25

26 phology, size, spin vector) represented by the population studies in Table 3, it is interesting to consider that these offsets may be due to currently unquantified differences in the η values returned from fitting near-ir data (< 2.5µm) versus those obtained by fitting longer wavelength data (> 2.5µm). If this is the case then our results provide insight into the magnitude of these offsets of order several 10 s of percent. However, our limited sample makes it difficult to determine whether this offset is simply a consequence of small number statistics or due to phenomena not treated in the model such as wavelength dependent emissivity. Table 3 NEATM η α linear-fit coefficients Reference Intercept Slope Delbo (2004) 0.92 ± ± Wolters et al. (2008) 0.91 ± ± Mainzer et al. (2011) ± ± This work; FG3 from ± ± This work; FG3 from ± ± A second explanation for the mismatch between our data and the retrograde models is the unlikely possibility that FG3 is actually in a prograde state. The analysis of light curve photometry strongly suggests that FG3 s satellite is in a retrograde orbit (Scheirich et al., 2015) and the satellite seems to be tidally locked (Benner et al., 2012). Attempts to model the spin pole of the primary have most recently been inconclusive (Scheirich et al., 2015), but previous analyses suggest a retrograde state (Yu et al., 2014). While a retrograde satellite orbit and a prograde primary spin state seem unlikely from the perspective of binary formation models involving YORP spin-up and mass loss (e.g. Jacobson & Scheeres, 2011), the current ambiguity in the interpretation of the lightcurve data leaves this as an open possibility. 26

27 Future observations to unambiguously determine the primary spin state would directly test this scenario. Interestingly, the linear fits to our FG3 data are significantly different for each of the two epochs of observation (Table 3). Furthermore, the reflectance spectra of FG3 are distinctly different in each of the two epochs (Figure 4). Perhaps the most interesting comparison to highlight these differences is between the 2009 April 27 and 2011 December 1 data. These data clearly fall into the categories of neutral slope, low thermal emission in 2009, and red slope, high thermal emission in However the heliocentric distance of FG3 was almost identical on these two dates. The primary difference between these observations is viewing geometry: the solar phase angles are essentially ±50 on either side of opposition (Table 2) and the sub-earth latitudes on the surface of the body were at their maximum separation amongst the data considered here (Fig. 5). This raises a third and ultimately testable interpretation of the FG3 results. As shown in the bottom panel of Figure 5 the two epochs corresponded to two different viewing aspects: in 2009 the asteroid was at negative ecliptic latitude and thus observations preferentially viewed its northern hemisphere, whereas in 2011 the asteroid was at positive ecliptic latitude and thus observations preferentially accessed its southern hemisphere. The variability in optical and thermal properties could thus be suggestive of a heterogeneous surface. The top-shaped morphology of FG3 with its pronounced equatorial ridge (Benner et al., 2012; Yu et al., 2014) could serve to isolate the optical and thermal influence of the northern and southern hemispheres and thus accentuate surface heterogeneity even though the difference in viewing aspect was only across the two epochs. If the surface of FG3 is hetero- 27

28 geneous, then we predict that the northern hemisphere would have optically neutral reflectance properties and a higher relative albedo, which would lower surface temperature and yield higher η. Conversely, the southern hemisphere would have optically red reflectance properties and a lower relative albedo, which would result in a higher relative surface temperature and correspondingly lower η. FG3 will be a strong optical target in 2022 and will be at negative ecliptic latitudes, thus providing a vantage of its northern hemisphere. Additional long term monitoring of this object may ultimately determine whether its measured variability is due to a heterogeneous surface (de Leon et al., 2011, 2013) or some other physical effect (Binzel et al., 2012). Spectral heterogeneity is not commonly seen on asteroids, particularly for km-scale objects like FG3, however near-ir observations of the 500-meter near-earth object Bennu indicate variations in spectral slope that may be related to differences in regolith grain size across the surface (Binzel et al., 2015). In light of the uncertainty surrounding our FG3 results it remains unclear whether physical properties like thermal inertia and pole orientation can be robustly derived from NEATM modeling of thermal emission in K-band. The highly non-uniform distribution of NEA obliquity, specifically a preference toward purely retrograde states (La Spina et al., 2004), would suggest that this technique should be well suited to detection of aspect-dependent thermal variability in the NEA population. Following on the work of Delbo (2004), new models to quantify the detectability of thermal asymmetry as a function of obliquity would provide new insight into the feasibility of such observational studies. Additionally, a larger sample size is needed to definitively determine whether NEATM modeling at near-ir wavelengths can robustly retrieve physical properties like pole orientation and thermal inertia. 28

29 It is clear that any future near-ir-based study of thermal inertia and pole orientation would require targets that span the widest possible range of phase angles on either side of opposition. On any given night 4, approximately 20% of observable NEAs (V < 20, solar elongation > 60, galactic latitude > 20 ) are at phase angles greater than 50. Some of this uncertainty could be resolved by extending such phase dependent measurements to longer wavelengths, for instance photometry at L- and M-bands (3.5 and 4.5 µm respectively) or in the N- and Q-bands ( 10 and 20 µm respectively) which are closer to the thermal emission peaks typical of NEOs. However, a systematic comparison of NEATM results based on data in different wavelength regimes is needed to more completely understand the limitations of this model. Extending measurements to longer wavelengths would have a primary benefit of enabling detection and modeling of thermal emission for any NEA, not just low albedo objects like QE2 and FG3. A targeted spectro-photometric survey to monitor phase-dependent thermal emission from select NEAs could produce an unprecedented sample of spin vector orientations. Such a sample would provide directs tests for models of NEA source regions in the Main Belt that depend on spin orientation and the subsequent direction of semi-major axis drift due to the Yarkovsky effect (e.g. Bottke et al., 2002). Acknowledgements We are grateful to Eric Volquardsen for use of his NEATM thermal model code. This manuscript was significantly improved thanks to insightful reviews from Michael Mueller, Al Harris, and an anonymous referee. This work in- 4 Based on a filtering of NEAs with the observability tool at 29

30 cludes data obtained at NASA s IRTF located on Mauna Kea in Hawaii, which is operated by the University of Hawaii under contract NNH14CK55B with the National Aeronautics and Space Administration. Support for this project was provided through the National Science Foundation Astronomy and Astrophysics Postdoctoral Fellowship awarded to N.A.M and by NASA grant number NNX14AN82G issued through the Near-Earth Object Observations program. F.E.D. acknowledges support by the National Science Foundation under Grant and by NASA under Grant NNX12AL26G and through Hubble Fellowship Grant HST-HF A awarded by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., for NASA, under Contract NAS D.P. is thankful to the AXA research fund for support through their postdoctoral fellowship. E.S.H., R.J.V., abd Y.R.F. acknowledge partial support from National Science Foundation grant AST References Benner, L. A. M. and 9 co-authors, Arecibo and Goldstone Radar Observations of Binary Near-Earth Asteroid and Marco Polo-R Mission Target (175706) 1996 FG3. American Astronomical Society, DPS meeting #44, abstract Binzel, R. P., Harris, A. W., Bus, S. J., and Burbine, T. H., Spectral Properties of Near-Earth Objects: Palomar and IRTF Results for 48 Objects Including Spacecraft Targets (9969) Braille and (10302) 1989 ML. Icarus 151, Binzel, R. P., Polishook, D., DeMeo, F. E., Emery, J. P., Rivkin, A. S., Marco Polo-R target asteroid (175706) 1996 FG3: Possible evidence for an 30

31 annual thermal wave. Lunar Planet. Sci. 43. Abstract Binzel, R. P. and 16 co-authors, Spectral slope variations for OSIRIS- REx target Asteroid (101955) Bennu: Possible evidence for a fine-grained regolith equatorial ridge. Icarus 256, Bottke, W. F., Morbidelli, A., Jedicke, R., Petit, J. M., Levison, H. F., Michel, P., Metcalfe, T., Debiased orbital and absolute magnitude distribution of near-earth objects. Icarus 156, Bottke, W. F., Vokroulický, D., Rubincam, D. P., Nesvorný, D., The Yarkovsky andyorp Effects: Implications for asteroid dynamics. AREPS 34, Burbine, T. H., Could G-class asteroids be the parent bodies of CM chondrites? MAPS 33, Carry, B. and 10 co-authors, Physical properties of (2) Pallas. Icarus 205, Delbo, M., Harris, A. W., Binzel, R. P., Pravec, P., and Davies, J. K., Keck observations of near-earth asteroids in the thermal infrared. Icarus 166, Delbo, M., The nature of near-earth asteroids from the study of their thermal infrared emission. Doctoral thesis, Freie Universitat Berlin. Delbo, M., dell Oro, A., Harris, A. W., Mottola, S., Mueller, M., Thermal inertia of near-earth asteroids and implications for the magnitude of the Yarkovsky effect. Icarus 190, Delbo, M. and Tanga, P., Thermal inertia of main belt asteroids smaller than 100 km from IRAS data. Planetary and Space Science 57, Delbo, M., Mueller, M., Emery, J. P., Rozitis, B. E. and Capria, M. T., Asteroid thermophysical modeling. In: Michel, P., DeMeo, F. E., Bottke, W. F. (Eds.), Asteroids IV. Univ. Arizona Press, Tucson, pp

32 de Leon, J., Mothe-Diniz, T., Licandro, J., Pinilla-Alonso, N., and Campins, H., New observations of asteroid (175706) 1996 FG3, primary target of the ESA Marco Polo-R mission. A&A 530, L12. de Leon, J., Lorenzi, V., Ali-Lagoa, V., Licandro, J., Pinilla-Alonso, N., and Campins, H., Additional spectra of asteroid 1996 FG3, backup target of the ESA MarcoPolo-R mission. A&A 556, A33. DeMeo, F. E., Binzel, R. P., Slivan, S. M., and Bus, S. J., An extension of the Bus asteroid taxonomy into the near-infrared. Icarus 202, Farinella, P. and Vokrouhlický, D., Semimajor axis mobility of asteroidal fragments. Science 283, Gehrels, T., Minor Planets. I. The rotation of Vesta. AJ 72, Hansen, O., On the prograde rotation of asteroids. Icarus 32, Harris, A. W., A thermal model for near-earth asteroids. Icarus 131, Harris, A. W., Mueller, M., Delbo, M., and Bus, S. J., Physical characterization of the potentially hazardous high-albedo Asteroid (33342) 1998 WT24 from thermal-infrared observations. Icarus 188, Hicks, M., Lawrence, K., Chesley, S., Chesley, J., Rhoades, H., Elberhar, S., Carcione, A., Borlase, R., Palomar Spectroscopy of Near-Earth Asteroids (1999 KX4), (1999 JV3), (2002 GT), (2002 OD20), and (1998 QE2). The Astronomer s Telegram 5132, 1. Howell, E. S. and 7 co-authors, SHERMAN: A Shape-based Thermophysical Model for Near-Earth Asteroids. American Astronomical Society, DPS meeting #47, abstract Jacobson, S. A. and Scheeres, D. J., Dynamics of rotationally fissioned asteroids: Source of observed small asteroid systems. Icarus 21, Johnson, P. E., Kemp, J. C., Lebofsky, M. J., and Rike, G. H., µm 32