Issue 
A&A
Volume 554, June 2013



Article Number  A71  
Number of page(s)  16  
Section  Planets and planetary systems  
DOI  https://doi.org/10.1051/00046361/201220680  
Published online  05 June 2013 
Physical properties of Btype asteroids from WISE data^{⋆}
^{1} Instituto de Astrofísica de Canarias (IAC), c/ Vía Láctea s/n, 38205 La Laguna, Tenerife, Spain
email: vali@iac.es
^{2} Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
^{3} Departamento de Edafología y Geología, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
^{4} UNS–CNRS–Observatoire de la Côte d’Azur, BP 4229, 06304 Nice Cedex 4, France
^{5} Physics Department, University of Central Florida, PO Box 162385, Orlando, FL 32816.2385, USA
^{6} Instituto de Astrofísica de Andalucía (IAA), 18008 Granada, Spain
^{7} Department of Earth and Planetary Sciences, University of Tennessee, 1412 Circle Dr, Knoxville TN 379961410, USA
^{8} Department of Astronomy, University of Maryland, College Park, MD 204722421, USA
Received: 31 October 2012
Accepted: 19 March 2013
Aims. Our aim is to obtain more information about the physical nature of Btype asteroids and extend previous work by studying their physical properties as derived from fitting an asteroid thermal model to their NASA’s Widefield Infrared Survey Explorer (WISE) data. We also examine the Pallas collisional family, a Btype family with a moderately high albedo in contrast to the large majority of Btypes.
Methods. We applied a combination of the nearEarth asteroid thermal model and a model of the reflected sunlight to WISE asteroid data in order to derive up to four parameters: effective diameter (D), the socalled infrared beaming parameter (η), ratio of infrared to visible albedo (R_{p} = p_{IR}/p_{V}), and visible geometric albedo (p_{V}).
Results. We obtained the effective diameter, geometric visible albedo, infraredtovisible albedo ratio, and beaming parameter for ≳ 100 Btypes asteroids and plotted the value distributions of p_{V}, R_{p}, and η (p̅_{V} = 0.07 ± 0.03, R̅_{p} = 1.0 ± 0.2, and η̅ = 1.0 ± 0.1). By combining the IR and visible albedos with 2.5 μm reflectances from the literature we obtained the ratio of reflectances at 3.4 and 2.5 μm, from which we found statistically significant indications that the presence of a 3μm absorption band related to water may be commonplace among the Btypes. Finally, the Pallas collisional family members studied (~50 objects) present moderately high values of p_{V}, p̅_{V} = 0.14 ± 0.05, which is significantly higher than the average albedo of Btypes. In addition, this family presents the lowest and most homogeneously distributed R_{p}values of our whole sample, which shows that this group is clearly different from the other Btypes, probably because its members are fragments likely originating from the same region of (2) Pallas, a particularly highalbedo Btype asteroid.
Key words: minor planets, asteroids: general / surveys / infrared: planetary systems
Tables 1, 4, and Appendices A and B are available in electronic form at http://www.aanda.org
© ESO, 2013
1. Introduction
The study of asteroids is important for gaining knowledge about the origin and evolution of our planetary system. Asteroids are relics of the solar system’s formation and the building blocks of the terrestrial planets. Primitive asteroids, i.e. those belonging to the socalled spectroscopic Ccomplex and having in general visible geometric albedo p_{V} ≲ 0.1 and featureless, flat visible spectra, are particularly relevant in this context. Thought to have formed farther away from the Sun than the other asteroid classes, primitive asteroids have experienced less heating and alteration processes and have a more pristine composition, potentially preserving crucial information about the early solar system. In addition, primitive asteroids play an important role in current exobiological scenarios since they delivered complex organic molecules to the early Earth. This organic matter is a prerequisite for the synthesis of prebiotic biochemical compounds that would subsequently lead to the emergence of life (Maurette 2006, and references therein). For these and other reasons, upcoming sample return space missions have selected primitive asteroids as primary targets: NASA’s OSIRISRex (Campins et al. 2010b; Lauretta et al. 2010), ESA’s Marco PoloR (Barucci et al. 2012; de León et al. 2011), and JAXA’s Hayabusa2^{1}.
The taxonomic classification of primitive asteroids has been traditionally based on their low visible albedo ( ≲ 0.08–0.1), relatively flat or slightly blue visible spectra with weak or no absorption features (for a detailed review, see Clark et al. 2010). Several primitive classes have been defined in Tholen’s taxonomy, e.g. B, C, F, G, D, and P (Tholen 1984, 1989). Bus’s featurebased classification, independent of the albedo, merged some of these and defined new primitive taxons that were extended with minor changes into the BusDeMeo taxonomy (Bus & Binzel 2002; DeMeo et al. 2009).
Btypes are of particular interest among the primitive asteroids for a number of reasons: (i) there is as yet no compelling explanation for their defining feature, i.e. their slightly blue spectral slope in the visible range; (ii) Btype asteroids constitute the only primitive class that presents a wide range of spectral slopes in the 0.8–2.5 μm nearinfrared (NIR) region (de León et al. 2012); (iii) the few Btypes studied present the 3μm absorption feature related to hydrated minerals; (iv) water ice has been detected on the surface of (24) Themis (Campins et al. 2010a; Rivkin & Emery 2010); (v) the majority of asteroids that have been observed to display cometarylike activity are Btypes (Licandro et al. 2012, and references therein); (vi) the target of NASA OSIRISRex mission, 2006 RQ_{36}, is a Btype asteroid.
The Btype asteroids have been widely related to carbonaceous chondrites, composed of carbonaceous minerals and phyllosilicates, in terms of their generally low albedo and broad spectral properties (Gaffey et al. 1989; Vilas & Gaffey 1989; Vilas et al. 1994). More recently, de León et al. (2012) have examined visible to nearinfrared (VNIR) spectra of a sample of 45 Btypes and found that the characteristic negative spectral slope in visible wavelengths diverges into a continuum of gradually varying NIR spectral slopes, from a monotonic negative (blue) slope to a positive (red) slope. De León et al. classified their spectra into six “average spectra” or “centroids” representative of the whole sample by means of statistical clustering analysis (Marzo et al. 2009). These centroids were compared against meteorite spectra from the RELAB database (Pieters & Hiroi 2004). The best meteorite analogues found for the six “clusters” were all carbonaceous chondrites with a gradual change in their degree of hydration, from aqueously altered CM2 chondrites for the reddest cluster to the heated/thermally metamorphosed CK4 chondrites for the bluest one.
This work is an extension of the results obtained by de León et al. (2012), which are part of a ongoing programme devoted to improving our knowledge of Btypes. Our aim here is to study the physical properties of Btype asteroids that can be derived by fitting a thermal model to their NASA’s Widefield Infrared Explorer (WISE) observations, i.e. effective diameter, beaming parameter, and p_{IR}/p_{V}, where p_{IR} is the albedo at 3.4–4.6 μm as defined in Mainzer et al. (2011b). We closely follow the methodology of Mainzer et al. (2011b), though with a number of differences, as described in Sects. 2 and 3.
Mainzer et al. (2011c) studied all groups of spectrophotometrically classified asteroids in the Tholen, Bus, and BusDeMeo taxonomies observed by WISE, including Btypes. Those authors present visible and NIR geometric albedo distributions and median values of Btypes and conclude that, in spite of having analogously low albedos, B, C, D, and Ttype asteroids can be distinguished from their values of NIR reflectance. In particular, Mainzer et al. (2011c) point out that Btypes have a lower p_{IR}/p_{V} ratio than Ctypes and attribute this to their blue VNIR slopes likely extending out to 3–4 μm. Our definition of Btype asteroid in this work is different: following Clark et al. (2010) and de León et al. (2012), we consider all objects that have a flat to slightly blue spectral slope in the visible range, i.e. any object that has ever been classified as Btype, including Tholen’s Ftypes and ambiguous designations. This criterion produces a total of 162 asteroids classified as Btypes.
We also study the collisional family of (2) Pallas (hereafter PCF and Pallas, respectively). The PCF is interesting for the following reasons: (i) it is a Btype family, given that Pallas and the very few family members that have been taxonomically classified are Btypes (GilHutton 2006; de León et al. 2010), and the five members studied in de León et al. (2012) were spectrally related to carbonaceous chondrites, which establishes their primitive nature; (ii) the average values of geometric albedo of members of the family calculated by Masiero et al. (2011) are roughly 0.15 (see their Fig. 19), significantly greater than expected for primitive bodies (<0.1), though no explicit comment is made by these authors on this intriguing result; (iii) the nearEarth asteroid (3200) Phaethon, an activated asteroid parent of the Geminid meteor shower, likely originated in the Pallas family (de León et al. 2010); (iv) this family is isolated well in (proper) element space, thus the potential identification of interlopers as members is greatly reduced.
The paper is organised as follows. In Sect. 2, we briefly describe the WISE data set and our selection criteria. The thermal modelling of the data is explained in detail in Appendix A, whereas Sect. 3 includes relevant comments on the very few differences introduced in this work. In Appendix B our parameter determinations are compared to those by Masiero et al. (2011). We present our results in Sect. 4, a discussion of the implications of this work is put forward in Sect. 5, and our conclusions are enumerated in Sect. 6.
2. Data
A general introduction to WISE can be found in Wright et al. (2010) and references therein. Of particular interest to solar system science is the NEOWISE project. This acronym collectively refers to two enhancements to the WISE data processing system that were designed to allow detection and archiving of solar system objects (for details, see Mainzer et al. 2011a).
WISE used four broadband filters with approximate isophotal wavelengths at 3.4, 4.6, 12, and 22 μm, referred to as W1, W2, W3, and W4, respectively (Wright et al. 2010). The WISE AllSky Single Exposure L1b Working Database, published in April 2012 and available via the IRSA/IPAC archive^{2}, includes the corresponding magnitudes and uncertainties in the Vega system, as well as quality and contamination and confusion flags that enable us to reject defective data (Cutri et al. 2012).
We follow a combination of criteria found in Mainzer et al. (2011b,c), Masiero et al. (2011), and Grav et al. (2012) to ensure the reliability of the data. We implement the correction to the red and blue calibrator discrepancy in W3 and W4, and we use a cone search radius of 0.3′′ centred on the MPC ephemeris of the object in our queries. All artefact flags other than p, P, and 0 and quality flags other than A, B, and C are rejected, and we require the modified Julian date to be within four seconds of the time specified by the MPC and split groups of epochs separated more than three days (see the end of this section). We ensure that the data is not contaminated by inertial sources by removing those points that return a positive match from the WISE Source Catalog within 6′′. Finally, all remaining observations in a given band are rejected if they are fewer than 40% of the data in the band with the maximum number of detections.
On the other hand, we do not use data saturated to any extent. The onset of saturation is reported to correspond to magnitudes M_{W1} < 6, M_{W2} < 6, M_{W3} < 4, M_{W4} < 3 (Cutri et al. 2012). We found that enlarging the error bar of partially saturated data to 0.2 magnitudes (which translates into a relative error of 20% in fluxes) renders the corresponding band unable to play any effective role in the thermal model fit by not contributing significantly to the χ^{2}.
The application of the above criteria results in a sample of 111 Btype mainbelt asteroids with WISE observations usable for our purposes. Some of these have been observed by WISE in more than one uninterrupted group of epochs with different observation geometries. We also model such groups of observations separately if they are more than three days apart (see Appendix A and Mainzer et al. 2011b). Consequently, we have a larger set of parameter determinations than asteroids in the sample.
3. Thermal modelling
The modelling of WISE asteroid data implemented in this work closely follows Mainzer et al. (2011b,c), and Masiero et al. (2011) and is based on the nearEarth asteroid thermal model (NEATM, Harris 1998) and the IAU phase curve correction to the visible magnitude (Bowell et al. 1989). For the sake of reproducibility, we include a detailed account of our procedure in Appendix A and enumerate the few differences with respect to Masiero et al. and Mainzer et al. below.
The number of parameters we can fit for each object depends upon how many and which WISE bands are present in its data set. Parameter default values are chosen based on the peak of their respective fitted value distributions of main belt asteroids presented in Masiero et al. (2011). Whenever there is one single or no thermal band available (W2, W3, or W4) we assume η = 1.0; R_{p} is fixed to 1.5 unless we have at least a 50% contribution of reflected sunlight in W1 data. The last criterion is based on the consistency of our parameter determinations for objects with double detections^{3}; namely, both groups of observations of asteroids (1076) and (2446) have >50% sunlight and their R_{p} values are consistent within the error bar. On the other hand, (3579) has incompatible R_{p} determinations from W1 data with >70% and ~25% reflected sunlight, respectively. We thus reject three R_{p} values belonging to asteroids (288), (1493), and (3579).
These considerations allow asteroid size to be fitted in all cases and, by means of the relation (1)the geometric visible albedo can be computed. In contrast, W1 and W2 data are more often rejected based on the data requirements (see Sect. 2) than the purely thermal bands, and one will usually be able to obtain fewer R_{p} determinations than η or indeed D.
It is important to point out that we do not use physical data previously determined by direct measurements, such as radar diameters or albedos, to constrain our fits; i.e., we limit ourselves to using radiometrically derived sizes. This will introduce variations in the parameter determinations of some individual asteroids – especially the largest ones, since it is more likely that more direct measurements have been performed – as compared to Masiero et al. (2011), but should not affect the result of statistical analyses if the populations studied consist of a significant number of objects. In Appendix B we show that, given the same input values of H, our bestfit parameter values are consistent within the error bars, though we find that our R_{p} determinations are systematically lower by ~10%. On the other hand, the update of ~45% of the MPC Hvalues in our sample does change the values of p_{V} (see Eq. (1)) and R_{p}. As shown in Fig. 1, the updates tend to be toward higher values of H, which will result in lower p_{V} and higher values of R_{p} than those of Masiero et al. (2011). For more details, see Appendix B.
Fig. 1 Differences in absolute magnitude values between those used by Masiero et al. (2011), H_{M}, and the most updated ones (as of May 2012) used in this paper, H_{U}. Note that the cases verifying ΔH = 0 do not contribute to this histogram. 

Open with DEXTER 
4. Results
4.1. Value distributions of η, p_{V}, and R_{p}
The distributions of η, p_{V}, R_{p}, and p_{IR} obtained for the Btypes are shown in Fig. 2 (complete set of parameter determinations in Table 1). For comparison, we overplotted the corresponding histograms with bestfit parameter values from Masiero et al. (2011) (see Appendix B for a detailed comparison).
Fig. 2 Histograms of value distributions obtained for the Btype asteroids observed by WISE: η (top left), p_{V} (bottom left), R_{p} (top right), and p_{IR} (bottom right). The corresponding histograms using the results from Masiero et al. (2011) are overplotted. 

Open with DEXTER 
The parameter median and mean values and standard deviations of this work, as well as the number of parameter determinations obtained in each case (N), are presented in Table 2. Note that p_{V} and p_{IR} are not fitted, but computed. The former is obtained from Eq. (1) with the bestfit value of D as input, whereas p_{IR} = R_{p}p_{V}. These results are consistent with previous work by Mainzer et al. (2011c): if we take the weighted mean of median p_{V} and R_{p}values corresponding to their Tholen, Bus, and BusDeMeo Btypes and Tholen Ftypes, we obtain the same median values.
We find a ηvalue distribution centered at unity, consistent with the average value obtained for the whole main belt (Masiero et al. 2011). The broad and asymmetrical p_{V} distribution extends to p_{V} > 0.1. The R_{p} distribution is also broad, whereas the values of p_{IR} are more compactly distributed around the mean.
Median and mean values and standard deviations of η, R_{p}, p_{V}, and p_{IR} derived for the Btype asteroids observed by WISE.
4.2. Albedo ratio
Figure 3 shows a plot of R_{p} versus p_{V}. Similar plots including all taxonomic classes in different classification schemes are presented in Figs. 14 and 15 by Mainzer et al. (2011c) to show how clearly different taxons may be distinguished. Here we concentrate on the p_{V} < 0.18 range, with all Btypes with WISE data for which R_{p}values could be derived are plotted (see Table 2); we also include all main belt asteroids taken from Table 1 of Masiero et al. (2011). The cloud of points exhibits a characteristic “waningmoon” shape, with no points in either the upper right or lower lefthand regions of the plot. Mainzer et al. caution that, while WISE is essentially unbiased against p_{V}, spectroscopic surveys conducted to create the classification schemes are inherently biased against small, lowp_{V} objects; in addition, the computation of R_{p} from WISE data requires sufficient reflected sunlight contribution in bands W1 and W2, which will tend to exclude objects with low enough values of p_{V} and R_{p}. This could explain the lack of points in the lower lefthand part of the plot. However, these biases cannot be solely responsible for this characteristic shape since other taxonomic classes with higher values of p_{V} and R_{p} also cluster similarly. Furthermore, if we plot all main belt objects irrespective of whether they have a taxonomic classification or not, distinct clouds of points with the same shape become apparent. Thus, we emphasise that because of their characteristic p_{V} and p_{IR}value distributions, there are no highR_{p} objects among the highp_{V} Btype asteroids.
Fig. 3 Albedo ratio versus visible geometric albedo. The 20% error bars in p_{V} and R_{p} are not shown to facilitate visualisation. Btypes observed by WISE are plotted in black circles; all main belt objects featured in Table 1 of Masiero et al. (2011) are in empty grey circles. The horizontal line is an artificial feature corresponding to objects with the default fixed value of R_{p} = 1.5. 

Open with DEXTER 
We have also analysed the R_{p}values of the sample of 45 Btypes studied by de León et al. (2012) separately. The spectra of these asteroids were classified into six “average spectra” or “centroids” referred to as G1, G2, ... G6 (see Sect. 1). These show a progressive decrease in spectral gradient in the NIR interval (0.8–2.5 μm), ranging from a positive (red) slope for G1 to a negative (blue) slope for G6. In Fig. 4, we plot R_{p} versus p_{V} labelling the objects in the different centroids G1, ... G5. Cluster G6 is not included since it is only composed of one member, (3200) Phaethon, which did not have enough WISE observations to perform a reliable fit.
Fig. 4 R_{p} vs. p_{V} for the asteroids in de León et al. (2012) distinguishing the clusters to which they belong. Within parentheses, the average p_{V} of each cluster is shown. On average, p_{V} increases from G1 to G3, keeping the same value from G3 to G5. 

Open with DEXTER 
The average values of p_{V} for each cluster increases from G1 to G3, retaining the same value from G3 to G5. This might suggest an inverse correlation between the cluster NIR slope and p_{V}, though the small number of objects per cluster with WISE observations (2 objects in G1, 12 in G2, 6 in G3, 5 in G4, and 3 in G5) prevents establishing a firm conclusion.
In Fig. 5, the average R_{p}value is plotted for the different clusters.
Fig. 5 Average values of R_{p} for each cluster defined by de León et al. (2012). In terms of spectral slope up to 2.5 μm, cluster G5 is the bluest, whereas G1 is the reddest. The error bar is the standard error of the mean. 

Open with DEXTER 
This figure suggests a correlation between the average R_{p} and the NIR slope of the relative reflectance of the clusters: objects with higher R_{p} belong (on average) to clusters with higher NIR spectral slope in the ~ 1.0–2.5 μm range (see also Fig. 5 of de León et al. 2012). This indicates that the reflectivity at 3.4 μm tends to continue the trend observed at shorter IR wavelengths, as hypothesised by Mainzer et al. (2011c).
The W1 band pass spans the 2.8 to 3.8 μm range (Wright et al. 2010); therefore, the values of R_{p} may also indicate the presence of the 3μm absorption feature attributed to hydrated minerals or water ice detected on many asteroids (Rivkin et al. 2000; Gaffey et al. 2002; Campins et al. 2010b; Rivkin & Emery 2010; Licandro et al. 2011). Mainzer et al. (2011c) ruled out the possibility of detecting the hydration band from WISE data based on the fact that the average R_{p} for a sample of 7 Mtypes with positive detections of the band (Rivkin et al. 2000) cannot be distinguished from those corresponding to other 33 Mtypes. However, this test might not be meaningful given that the R_{p}values have error bars that are at least larger than the characteristic depth of the absorption feature. Below we provide evidence that the W1 may be sensitive to the 3μm feature.
Our first step was to assemble VNIR spectra up to ~ 3.6 μm of a list of nine Ccomplex control asteroids. All of these asteroids are primitive, including some Btypes, such as (2) Pallas or (45) Eugenia, and some of them show a distinct absorption feature. The assembled spectra are plotted along with the bestfit value of R_{p} in Fig. 6. We collected or digitised data from Hiroi et al. (1996), SMASSII (Bus & Binzel 2002), the 52Color Survey (Bell et al. 2005), Rivkin et al. (2003) and de León et al. (2012).
Fig. 6 Assembled spectra of the control objects chosen to study the reflectivity at 3.4 μm derived from WISE data as a diagnostic of the presence of the 3μm feature. The dashed lines approximately enclose the wavelength integration range of the W1 filter. 

Open with DEXTER 
Taking into account the errors associated with superimposing spectra obtained at different epochs of observation with different equipment and the uncertainties in the spectra and in the value of the albedo ratio, the latter is not expected to exactly match the value of the spectra at 3.4 μm. Owing to these deviations, seen in Fig. 6, one can neither confirm nor rule out the presence of the absorption feature based on the value of R_{p} alone.
Next, we combined the values of R_{p} with the observed relative reflectances at 2.5 μm available from other datasets. Rivkin et al. (2003) use the parameter 1 − R_{λ}/R_{2.5} as a rough measure of band depth. Because the band minima are usually near 3.0 μm and given that W1 results from an average over 2.8 to 3.8 μm, the parameter b ≡ 1 − R_{p}/R_{2.5} is not to be taken as a measure of band depth but as a helpful parameter to quantitatively compare the values of relative reflectances at 2.5 and 3.4 μm. If a deep absorption band is present, R_{p} is in general expected to be ≲ R_{2.5}, hence b ≳ 0, as is the case for (19) Fortuna. There is one general case for which this interpretation would be wrong: if the NIR slope is negative up to 3.4 μm, we would have b ≳ 0 even if no band was present, though blue NIR slopes such as Pallas’ have only been measured in a low percentage of cases (e.g. 2 out of 45 in the sample of de León et al. 2012).
A plot of b vs. p_{V} is shown in Fig. 7 with an estimated b error bar of 0.2 to account for the large uncertainties in R_{p} and the assembly of the spectra. We see that those asteroids with a weak or nonexistent absorption band tend to have b ≲ 0, whereas b ≳ 0 for those with a higher contrast feature. Only in the case of (13) Egeria is the value of b clearly inconsistent with the absortion band observed. Indeed, the large error bar in the value of b still prevents irrefutable detection of the band on a casebycase basis, but from Fig. 7, the correlation between the sign of b and the spectra appears robust, so that systematically obtaining b ≳ 0 for a given population may be statistically significant.
In Fig. 8 (top panel) we show a bvalue histogram of those Btype asteroids for which the value of b could be determined, and R_{2.5} was taken from the Btype spectra presented in de León et al. (2012), after normalising to unity at 0.55 μm. To test for statistical significance, we resort to the KolmogorovSmirnov (KS) test, which enables one to reject the null hypothesis that a given set of unbinned values is compatible with having been drawn from a given distribution function (see e.g. Press et al. 1986). As the null hypothesis we take a Gaussian distribution with zero mean. This choice is based on the mean value of b that we would expect considering that interpolating between R_{0.55} = 1 and gives (see Fig. 9). The KS test amply rules out the null hypothesis that the Btypes bvalues are drawn from a Gaussian distribution of zero mean regardless of its width (σ).
To demonstrate the robustnest of this result even further, we carried out the same procedure for a list of asteroids belonging to the Scomplex, including all BusDeMeo pure Stypes for which we found that the value of R_{2.5} was available from the literature and for which R_{p} could be computed, plus enough randomly selected S subtypes to get the same number of bvalues as we derived for the Btypes (see Table 5). Stype asteroids are “anhydrous” and have positive spectral slopes up to 2.5 μm so that in the absence of a 3μm absorption feature, one would expect a negative value of b if the spectral slope maintains its trend up to 3.4 μm. Taking into account that for this sample of Stypes, from the interpolated value (see Table 3), the expected mean value of b would be (see Fig. 9). The distribution of bvalues obtained for the Stypes, shown in the bottom panel of Fig. 8, presents a negative mean value. In this case, the KS test does not rule out the null hypothesis that the bvalues are drawn from a Gaussian distribution centered at − 0.14 with σ = 0.2 with pvalue >0.9. This rules out the possibility that a systematic error is causing the bvalues of the Btype asteroids to be >0.
In the preceeding section we have shown that our values of R_{p} are systematically ~ 10% lower than those of Masiero et al. (2011). We carried out the same procedure enlarging our R_{p} values 10% and the conclusions still hold.
To conclude, we find that the majority of Btypes with computed bvalues verify b > 0 and that very few present a clearly negative value of b, which means that even for those Btype spectra with a positive slope in the 2.5 μm region (approximately half of the objects in de León et al. 2012), there is a reduction in the reflectivity around 3.4 μm. As discussed above, for asteroids of the Ccomplex, the 3μm absorption feature has been attributed to hydrated minerals or water ice. The presence of goethite has been also proposed as an alternative explanation for this band (Beck et al. 2011). Nonetheless, while other closely related minerals have been found in both meteorite and asteroid spectra, extrarrestrial goethite has never been identified within the meteorite inventory, so the possibility that putative goethitecontaining asteroids never found a dynamical collisional pathway to Earth is less likely than the simpler interpretation: goethite is not present in asteroidal surfaces (Jewitt & GuilbertLepoutre 2012). Therefore, from the distribution of bvalues, we conclude that most asteroids in this sample (which constitutes ≲ 40% of the Btype population with computed R_{p}) present this absorption and that “water” (be it bound or free) may be common among the Btype asteroids.
4.3. The Pallas collisional family (PCF)
We use the most uptodate Pallas family list by (Nesvorny 2012). WISE has observed 46 of the objects in this list. Histograms of beaming parameter and albedo determinations are shown in Fig. 10 (complete parameter set in Table 4). Given that the R_{p} value could only be fitted for seven objects, we do not include histograms for R_{p} and p_{IR}. In Table 6 we present the mean values of η, p_{V}, and R_{p}, along with their corresponding standard deviations and contributing number of determinations.
Fig. 7 b = 1 − R_{p}/R_{2.5} versus p_{V} for a set of primitive control objects. The values of b obtained from the combination of our R_{p} and the R_{2.5} taken from the assembly of spectra published by other authors (cf. Fig. 6) and the shape of the 3 μm absorption band are consistent except for (13) Egeria. 

Open with DEXTER 
Fig. 8 Normalised histograms showing the distribution of bvalues. Top panel: Btypes; the fact that b > 0 in the majority of cases suggests that a high percentage of these asteroids have absorption features in the ~ 3μm region. Bottom panel: Stypes; the opposite conclusion is reached, consistent with the anhydrous nature of Stypes. 

Open with DEXTER 
The average η value of the PCF is consistent with that of the Btype population and that of the main belt (Masiero et al. 2011). However, the PCF has a moderately high albedo of ~ 14%, significantly higher than the rest of the Btypes and indeed higher than the value expected for primitive asteroids, whereas their average R_{p} is lower than the average value for the Btypes (cf. Table 2). The seven R_{p}values computed are quite homogeneous compared to the distribution observed for the rest of the Btype population. As we saw in Sect. 4.1, our parameter determinations are consistent within the error bars with those of Masiero et al. (2011), so we do not perform the same comparison for the PCF members. Masiero et al. show a histogram of log p_{V} with peak between log p_{V} = ( − 1) – ( − 0.8), i.e. p_{V} = 0.10–0.16, with which our results are consistent. (They provide no mean value to compare with.) On the other hand, though they provide the bestfit parameters in their Table 1, Masiero et al. do not discuss the R_{p} values of the PCF.
The average p_{V} and R_{p} values continue the trend observed for the clusters of de León et al. (2012) in Figs. 4 and 5: they have higher visible albedos on average than the G5 cluster and are bluer at 3.4 μm than the bluest cluster, G4. (The average R_{p} value of cluster G5 is obtained from only three objects, which could explain its deviation from the trend.) Members of the Pallaslike group of objects of Clark et al. (2010) were concentrated in clusters G4 and G5 in de León et al. (2012) (see Sect. 1). In this work, the list of Pallas family members observed by WISE includes still more asteroids. The four additional objects with R_{p} determinations that were not present in de León et al. (2012) are also located in the lower righthand region of Fig. 4. This property would be consistent with the characteristics of an extrapolated G6 centroid (see the right panel of Fig. 11). Interestingly, the G6 centroid of de León et al. (2012) contained asteroid (3200) Phaethon alone. Therefore, the connection of the NEA (3200) Phaethon to the PCF established by de León et al. (2010) based on spectroscopical and dynamical arguments is also supported by the values of relative reflectances at 3.4 μm.
5. Discussion
In Sect. 4.2 we present statistically significant indications that water may be common within our sample, and in Sect. 4.3 we point out how the PCF have higher geometric visible albedos and more homogeneously distributed R_{p}values than the rest of the Btypes. The purpose of this section is to put these results in the context of other studies.
In the context of the geophysical models of the Themis and Pallas parent bodies by CastilloRogez & Schmidt (2010) and Schmidt & CastilloRogez (2012), the detections of water ice and organics on the surface of (24) Themis (Rivkin & Emery 2010; Campins et al. 2010b) and (65) Cybele (Licandro et al. 2011) have been invoked as provocative indications that water has played an important role in the accretion and evolution of these asteroids’ parent bodies in the midouter asteroid belt. Finding indications that a 3μm is frequently present in Btypes spectra adds further support to this conceptual framework, in which the gradual differences in the observed NIR spectral slopes of Btypes might also be explained as resulting from different subsequent evolution and processing of their parent bodies.
Schmidt & CastilloRogez (2012) discuss that (1) Ceres, being the largest body of the asteroid belt and plausibly a waterrich asteroid, mostly preserved its integrity, whereas the Themis parent body was catastrophically disrupted. Pallas would be an intermediate case, still intact but showing evidence of heating, internal evolution, impact, and loss of water. In accordance with this scenario, the Themis family members included in the de León et al. (2012) sample were distributed more or less homogeneously among the representative centroids of the complete sample, spanning from neutral to red slopes and matching several meteorite analogues, whereas the PCF members were concentrated in the bluest centroids. As our results show, these trends are also verified longwards of 2.5 μm: the Themis family members show more heterogeneous R_{p} values, whereas the Pallas family members have distinctly higher albedos and lower, more homogeneous values of R_{p} than the rest of Btypes. This homogeneity is consistent with the collisional family being the result of a nondisruptive cratering event on Pallas, which has a significantly higher albedo than the rest of large Btypes.
Fig. 9 Schematic diagram illustrating the interpolation of R_{2.5} (empty points) from the average values of R_{p} (filled points) for the B and S asteroids studied. 

Open with DEXTER 
Bestfitting values of physical parameters determined for the Stypes with sufficient W1 WISE observations and published 2.5 μm reflectances.
Fig. 10 Beaming parameter (left) and geometric albedo (right) distributions of the values derived for the PCF members observed by WISE. The albedo distributions are normalised, i.e. divided by their respective total number of counts. 

Open with DEXTER 
Fig. 11 Left panel: PCF members are on average brighter in the visible and bluer in the 3.4 μm region than the Btypes sample (cf. Fig. 3). The 20% error bars are not plotted for clarity. Right panel: the average R_{p} of the PCF continues the decreasing trend of the de León et al. (2012) clusters (cf. Fig. 5). 

Open with DEXTER 
The geophysical models by Schmidt & CastilloRogez (2012) favour a waterrich past for Pallas. Its near spherical shape may be explained as a consequence of early melting of a substantial initial icy component that is subsequently removed by a combination of thermal and impact processes. The waterloss processes that took place on the surface of Pallas may have taken place on other active bodies, such as main belt comets or indeed (3200) Phaethon, linked to the PCF by de León et al. (2010), and this activity has already been proposed to explain the characteristic NIR blue spectral slope (Schmidt & CastilloRogez 2012, and references therein). According to this model, waterrich materials are also expected to be important components within some members of Pallas’ family. The fact that the PCF members in this study are found to match the properties of the least hydrated clusters of de León et al. (2012) would be inconsistent with this view, but caution must be exercised when associating NIR slopes of meteorite analogues to hydration. For instance, while having suffered intense thermal metamorphism, CK4 chondrites (the best matches for Pallas) still show an absorption feature at 2.9 μm, albeit less prominent than e.g. CM chondrites.
From the dynamical standpoint, the simulations of Walsh et al. (2011) also provide a congruous context for our results. Their model suggests that Ccomplex asteroids were formed in the giantplanetforming region and that, for every Ctype planetesimal from beyond 8 AU that would later be located in the outer main belt, more than ten would have ended up in the region where terrestrial planets formed. Assuming that the composition of these objects is 10% water by mass, this may account for the minimal mass required to bring the current amount of water to our planet by a factor of 6–22.
On the other hand, the question of the moderately high values of p_{V} obtained for the PCF remains unanswered. Finding members of the Pallas and Themis collisional families distributed differently among the de León et al. (2012) clusters and in Fig. 11 (as discussed above) also leads us to ponder the possibility that a sequential or progressive physical process could explain the differences in the value distributions of p_{V} and R_{p}. Unfortunately, there are many possible mechanisms underlying the observed NIR spectral variability of primitive asteroids (e.g. regolith particle sizes, space weathering as a function of asteroidfamily age and composition, thermal processing) and these are difficult to disentangle (Ziffer et al. 2011).
Median and standard deviation of bvalues for Btypes and Stypes.
Mean values and standard deviations of η, R_{p}, and p_{V} derived for the members of the Pallas collisional family observed by WISE.
6. Conclusions
We have derived values of D, η, and R_{p} of 111 Btype asteroids by thermal modelling WISE data and updated H values (most of which have increased since the work of Masiero et al. 2011,as shown in Fig. 1). Additionally, from H and D we computed the corresponding values of p_{V} (Table 2). Our results agree within the error bars of the model with those previously published by Masiero et al. (2011). However, we obtained a systematic trend of approximately − 10% discrepancies in the determinations of R_{p} that we cannot explain, but it is most probably attributable to small differences in the tabulated solar flux data needed to estimate the reflected light component at 3.4 μm. This work led to the following conclusions.
 1.
We derived the distribution of η, p_{V}, and R_{P} fot the Btype asteroids (see Fig. 2) and obtained the following mean values: , , and .
There are no highp_{V}, highR_{p} Btype asteroids (p_{V} > 0.10, R_{p} > 1.0; see Fig. 3). The average R_{p}values of the centroids of de León et al. (2012) clearly decrease from G1 to G5, which implies a relationship between the IR slope of the asteroid spectra up to 2.5 μm (Figs. 4 and 5).
We computed bvalues for a set of Btype asteroids and Scomplex asteroids, which present and , respectively (Fig. 8). While the latter result is consistent with objects of the Scomplex being anhydrous, the former indicates that the majority of Btype asteroids bvalues in this study are consistent with the presence of a 3μm absorption feature usually attributed to hydrated minerals or water ice; therefore, water must have played a key role in the evolution of a large fraction of the Btypes, supporting recent works discussed in Sect. 5.
We also studied the Pallas collisional family. On the one hand, the average albedo ( 0.14 ± 0.05) of this familly is significantly higher than the average albedo of Btypes ( 0.07 ± 0.03) and moderately high compared to what is traditionally considered to be the albedo of primitive asteroids (<0.1). On the other hand, the albedo ratio values of the PCF members are very low and homogeneous ( 0.5 ± 0.1). These results clearly show the fundamental differences between this family and the rest of the Btypes.
Our results support the scenarios by the geophysical models by CastilloRogez & Schmidt (2010) and Schmidt & CastilloRogez (2012) and the simulations of Walsh et al. (2011), which suggest that water played an important role in the origin of primitive asteroid parent bodies of the midouter belt. In addition, the connection of the NEA (3200) Phaethon to the PCF established by de León et al. (2010) based on spectroscopical and dynamical arguments is also supported by the values of relative reflectances at 3.4 μm.
Acknowledgments
We thank the referee for a careful and constructive revision. V.A.L. acknowledges support from the project AYA201129489C0302 (MEC). J.L. acknowledges support from the projects AYA201129489C0302 and AYA201239115C0303 (MINECO). J.d.L. is grateful for financial support via a “Juan de la Cierva” contract from the Spanish “Secretaría de Estado de Investigación, Desarrollo e Innovación”. M.D.B. thanks the Space Situation Awareness programme of the European Space Agency (ESASSA) for financial support. N.P.A. was financed by the Spanish Ministry of Economy and Competitiveness through the “Juan de la Cierva” programme. This publication 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. This research also 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.
References
 Barucci, M. A., Cheng, A. F., Michel, P., et al. 2012, Exp. Astron., 33, 645 [NASA ADS] [CrossRef] [Google Scholar]
 Beck, P., Quirico, E., Sevestre, D., et al. 2011, A&A, 526, A85 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Bell, J. F., Owensby, P. D., Hawke, B. R., et al. 2005, NASA Planetary Data System, 28 [Google Scholar]
 Bowell, E., Hapke, B., Domingue, D., et al. 1989, in Asteroids II, eds. R. P. Binzel, T. Gehrels, & M. S. Matthews, 524 [Google Scholar]
 Bus, S. J., & Binzel, R. P. 2002, Icarus, 158, 146 [NASA ADS] [CrossRef] [Google Scholar]
 Campins, H., Hargrove, K., PinillaAlonso, N., et al. 2010a, Nature, 464, 1320 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Campins, H., Morbidelli, A., Tsiganis, K., et al. 2010b, ApJ, 721, L53 [NASA ADS] [CrossRef] [Google Scholar]
 CastilloRogez, J. C., & Schmidt, B. E. 2010, Geophys. Res. Lett., 37, 10202 [NASA ADS] [CrossRef] [Google Scholar]
 Clark, B. E., Ziffer, J., Nesvorny, D., et al. 2010, J. Geophys. Res. (Planets), 115, 6005 [NASA ADS] [CrossRef] [Google Scholar]
 Cutri, R. M., et al. 2012, VizieR Online Data Catalog, II/311 [Google Scholar]
 de León, J., Campins, H., Tsiganis, K., Morbidelli, A., & Licandro, J. 2010, A&A, 513, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de León, J., MothéDiniz, T., Licandro, J., PinillaAlonso, N., & Campins, H. 2011, A&A, 530, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 de León, J., PinillaAlonso, N., Campins, H., Licandro, J., & Marzo, G. A. 2012, Icarus, 218, 196 [NASA ADS] [CrossRef] [Google Scholar]
 Delbó, M., & Harris, A. W. 2002, Meteor. Planet. Sci., 37, 1929 [NASA ADS] [CrossRef] [Google Scholar]
 Delbó, M., Dell’Oro, A., Harris, A. W., Mottola, S., & Mueller, M. 2007, Icarus, 190, 236 [NASA ADS] [CrossRef] [Google Scholar]
 DeMeo, F. E., Binzel, R. P., Slivan, S. M., & Bus, S. J. 2009, Icarus, 202, 160 [NASA ADS] [CrossRef] [Google Scholar]
 Gaffey, M. J., Bell, J. F., & Cruikshank, D. P. 1989, in Asteroids II, eds. R. P. Binzel, T. Gehrels, & M. S. Matthews, 98 [Google Scholar]
 Gaffey, M. J., Cloutis, E. A., Kelley, M. S., & Reed, K. L. 2002, in Asteroids III, eds. W. F. Bottke Jr., A. Cellino, P. Paolicchi, & R. P. Binzel, 183 [Google Scholar]
 GilHutton, R. 2006, Icarus, 183, 93 [NASA ADS] [CrossRef] [Google Scholar]
 Grav, T., Mainzer, A. K., Bauer, J., et al. 2012, ApJ, 744, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Harris, A. W. 1998, Icarus, 131, 291 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Harris, A. W. 2006, in Asteroids, Comets, Meteors, eds. L. Daniela, M. Sylvio Ferraz, & F. J. Angel, IAU Symp., 229, 449 [Google Scholar]
 Hiroi, T., Zolensky, M. E., Pieters, C. M., & Lipschutz, M. E. 1996, Meteor. Planet. Sci., 31, 321 [NASA ADS] [CrossRef] [Google Scholar]
 Jewitt, D., & GuilbertLepoutre, A. 2012, AJ, 143, 21 [NASA ADS] [CrossRef] [Google Scholar]
 Lauretta, D. S., Drake, M. J., Benzel, R. P., et al. 2010, Meteor. Planet. Sci. Suppl., 73, 5153 [NASA ADS] [Google Scholar]
 Licandro, J., Campins, H., Kelley, M., et al. 2011, A&A, 525, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Licandro, J., Hargrove, K., Kelley, M., et al. 2012, A&A, 537, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Mainzer, A., Bauer, J., Grav, T., et al. 2011a, ApJ, 731, 53 [NASA ADS] [CrossRef] [Google Scholar]
 Mainzer, A., Grav, T., Masiero, J., et al. 2011b, ApJ, 736, 100 [NASA ADS] [CrossRef] [Google Scholar]
 Mainzer, A., Grav, T., Masiero, J., et al. 2011c, ApJ, 741, 90 [NASA ADS] [CrossRef] [Google Scholar]
 Marzo, G. A., Roush, T. L., & Hogan, R. C. 2009, J. Geophys. Res. (Planets), 114, 8001 [NASA ADS] [CrossRef] [Google Scholar]
 Masiero, J. R., Mainzer, A. K., Grav, T., et al. 2011, ApJ, 741, 68 [NASA ADS] [CrossRef] [Google Scholar]
 Maurette, M. 2006, Micrometeorites and the Mysteries of Our Origins (Berlin: SpringerVerlag) [Google Scholar]
 Neckel, H., & Labs, D. 1984, Sol. Phys., 90, 205 [NASA ADS] [CrossRef] [Google Scholar]
 Nesvorny, D. 2012, NASA Planet. Data Syst., 189 [Google Scholar]
 Pieters, C. M., & Hiroi, T. 2004, in Lunar and Planetary Institute Science Conference Abstracts, eds. S. Mackwell, & E. Stansbery, 35, 1720 [Google Scholar]
 Press, W. H., Flannery, B. P., & Teukolsky, S. A. 1986, Numerical recipes. The art of scientific computing [Google Scholar]
 Rivkin, A. S., & Emery, J. P. 2010, Nature, 464, 1322 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Rivkin, A. S., Howell, E. S., Lebofsky, L. A., Clark, B. E., & Britt, D. T. 2000, Icarus, 145, 351 [NASA ADS] [CrossRef] [Google Scholar]
 Rivkin, A. S., Davies, J. K., Johnson, J. R., et al. 2003, Meteor. Planet. Sci., 38, 1383 [NASA ADS] [CrossRef] [Google Scholar]
 Schmidt, B. E., & CastilloRogez, J. C. 2012, Icarus, 218, 478 [NASA ADS] [CrossRef] [Google Scholar]
 Tholen, D. J. 1984, Ph.D. Thesis, Arizona Univ., Tucson [Google Scholar]
 Tholen, D. J. 1989, in Asteroids II, eds. R. P. Binzel, T. Gehrels, & M. S. Matthews, 1139 [Google Scholar]
 Vilas, F., & Gaffey, M. J. 1989, Science, 246, 790 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Vilas, F., Jarvis, K. S., & Gaffey, M. J. 1994, Icarus, 109, 274 [NASA ADS] [CrossRef] [Google Scholar]
 Walsh, K. J., Morbidelli, A., Raymond, S. N., O’Brien, D. P., & Mandell, A. M. 2011, Nature, 475, 206 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 [NASA ADS] [CrossRef] [Google Scholar]
 Ziffer, J., Campins, H., Licandro, J., et al. 2011, Icarus, 213, 538 [NASA ADS] [CrossRef] [Google Scholar]
Online material
Bestfitting values of physical parameters determined for the Btypes with WISE observations.
Bestfitting values of physical parameters determined for the Pallas collisional family asteroids excluding (2) Pallas with WISE observations.
Appendix A: Thermal modelling of WISE asteroid data
Our aim is to model the observed asteroid flux as a function of several physical parameters and derive the set of parameter values that most closely reproduce the actually measured fluxes. In this work we follow the method described by Mainzer et al. (2011b). The set of wavelengths covered by WISE (specified in Sect. 2) allow us to derive up to three parameters by fitting a thermal model to asteroid WISE data: asteroid effective diameter, beaming parameter, and reflectance at 3.4 μm (defined below). Within the wavelength range covered, the observed asteroid flux consists of two components: (A.1)The thermal flux component (f_{th,λ}) is the main contribution to W3 and W4, whereas the reflected sunlight component (r_{s,λ}) dominates in band W1. In general, W2 will have nonnegligible contributions from both components (Mainzer et al. 2011b).
The computation of f_{th,λ} is based on the Near Earth Asteroid Thermal Model (NEATM; see Harris 1998; Delbó & Harris 2002). The asteroid is assumed to be spherical, and its surface is divided into triangular facets that contribute to the total thermal flux observed by WISE in accordance with the facet temperature (T_{i}), the geocentric distance (Δ), and the phase angle (α_{⊙}). In turn, the temperature of each facet depends on the asteroid heliocentric distance (r_{⊙}) and its orientation with respect to the direction towards the sun. It is given by (A.2)which results from assuming that each surface element δa_{i} is in instantaneous equilibrium with solar radiation. S_{⊙} is the solar power at a distance of 1 AU, A the bolometric Bond albedo, ϵ the emissivity (usually taken to be 0.9; see Delbó et al. 2007, and references therein), σ the StefanBoltzmann constant, and μ_{i} = cosθ_{i}, where θ_{i} is the angle between the normal to the surface element i and the direction towards the Sun. Nonilluminated facets will be instantaneously in equilibrium with the very low temperatures of the surroundings (~0 K), and thus their contribution to f_{th,λ} is neglected in the NEATM. Finally, the beaming parameter (η) can be thought of as a normalisation or calibration factor that accounts for the different effects that would change the apparent dayside temperature distribution of the asteroid compared to that of a perfectly smooth, nonrotating sphere (Harris 1998). These include, for example, the enhanced sunward thermal emission due to surface roughness (η < 1), or the nonnegligible nightside emission of surfaces with high thermal inertia that, in order to conserve energy, causes the dayside temperature to be lower than that compared to the ideal case with zero thermal inertia (η > 1).
The asteroid thermal flux component is then given by (A.3)where f_{i,λ} is the contribution from each illuminated facet of a 1km sphere; Ω ≡ (D/1 km)^{2} scales the crosssection of the latter to the corresponding value of an asteroid of diameter D. The colour correction associated with each value of T_{i} and each WISE band is applied to the facet flux. By definition, it is the quotient of the inband flux of the black body at the given temperature to that of Vega (Wright et al. 2010). A colour correction table was generated for all integer temperatures from 70 K up to 1000 K using the filter profiles available from Cutri et al. (2012).
The reflected light component, the second term on the righthand side of Eq. (A.1), is calculated as follows. First, the asteroid visible magnitude (V) that would be observed at a given geometry (r_{⊙}, Δ and α_{⊙}) can be estimated using the IAU phase curve correction (Bowell et al. 1989), along with the tabulated values of asteroid absolute magnitude (H) and slope parameter (G) from the Minor Planet Center. Secondly, knowledge of the solar visible magnitude and flux at 0.55 μm (V_{⊙} and f_{V⊙}, respectively) allows us to calculate the sunlight reflected from the asteroid at that particular wavelength: (A.4)If we assume that the Sun is approximated well by a blackbody emitter at the solar effective temperature (T_{⊙} = 5778 K), the estimated reflected flux at any other desired wavelength (r_{λ}) can be computed by normalising the black body emission B_{λ}(T_{⊙}) to verify r_{V}, i.e. (A.5)In this approximation, we can also consider (A.6)from which we arrive at the following expression: (A.7)where the subscript IR denotes 3.4 μm. We do not colourcorrect this component given the small correction to the flux of a G2V star (see Table 1 of Wright et al. 2010). Finally, to account for possible differences in the reflectivity at wavelengths longward of 0.55 μm, a prefactor to r_{λ} is included in the model, such that (A.8)This prefactor, R_{p}, is by definition equivalent to the ratio of p_{IR} and the the visible geometric albedo, so we will refer to it as the “albedo ratio”. The paremeter p_{IR} is the reflectivity at 3.4 and 4.6 μm defined by Mainzer et al. (2011b).
To sum up, the observed model flux can then be written as (A.9)We use the LevenbergMarquardt algorithm (Press et al. 1986) in order to find the values of asteroid size (, in km), beaming parameter (η) and albedo ratio (R_{p}) that minimise the χ^{2} of the asteroid’s WISE data set, namely (A.10)where F_{j,λ} and σ_{j,λ} are the measured fluxes and corresponding uncertainties, j runs over the observation epochs, and λ labels the WISE bands. The implementation of this technique involves calculating the partial derivatives of with respect to the fitting parameters, which is straightforward in the case of Ω and R_{p}. The partial derivative with respect to η can be derived from (A.11)
Appendix B: Comparison with Masiero et al. (2011)
Fig. B.1 Fractional difference histograms of D, η, p_{V}, and R_{p}. We define ε = 100(x − x_{M})/x, where x is the parameter value in this work and x_{M} the correspoding value taken from Table 1 by Masiero et al. (2011). The vertical lines mark the corresponding average values. Only parameters resulting from the same input values of H contribute to these histograms. 

Open with DEXTER 
Figure 2 shows that our parameter determinations and those of Masiero et al. (2011) are compatible in spite of the slight differences in the data set and the thermal modelling used in this work (refer to Sect. 2 and Appendix A), from which we do not expect to obtain exactly the same bestfit parameters for each object. In order to carry out a detailed comparison between our results and those of Masiero et al. (2011), we computed the mean fractional difference (ε) and corresponding standard deviations of D, η, p_{V}, and R_{p}. Let ε = 100(x − x_{M})/x, where x is the parameter value for a given object in this work, and x_{M} the correspoding value taken from Table 1 by Masiero et al. (2011). The distributions of ε values are plotted in Fig. B.1. These histograms only include parameter determinations that have the same H as input in order to identify possible discrepancies in results not caused by different values of H. We find that our values of D and η tend to be slightly higher by 1% and 3%, respectively, whereas our p_{V} values are lower by 2%, though these deviations are small compared to the error bars. On the other hand, there is a large bias towards lower values of R_{p} that, while still being within the error bar, must be addressed.
Most probably, the R_{p} discrepancy is associated with how the reflected flux r_{λ} is calculated. In particular, we take the solar flux at 3.4 μm (f_{IR⊙} in Eq. (A.7)) from the solar power spectrum at zero air mass of Wehrli^{4}, based on the one by Neckel & Labs (1984). Any differences in input, including solar visible magnitude, taken from tabulated data sources that may cause our r_{λ} to be systematically 10% greater than that of Masiero et al. (2011) would explain our higher values of R_{p}. For instance, considering that there is only one optimum value of r_{s,λ} to fit a given W1 data set, from Eq. (A.8) it is clear that larger r_{λ} will have associated a lower bestfit value of R_{p}.
Fig. B.2 Differences in albedo ratio determinations versus difference in absolute magnitude corresponding to the Btypes in this paper and those by Masiero et al. (2011). 

Open with DEXTER 
The Monte Carlo estimations by the NEOWISE team show that the error bars associated to the fitting of the data are always small compared to the errors inherent to the thermal model itself. The relative errors in diameters derived from the NEATM have been characterised as ~10%–15% (Harris 2006). From these facts and the widths of the εvalue distributions of Fig. B.1, we consider it safe to assume a minimum relative error of 10% in diameter and 20% in beaming parameter, p_{V} and R_{p}. On the other hand, large uncertainties in the absolute magnitude (sometimes as large as ~0.3 mag) will also affect the values of p_{V}, so 20% is probably an optimistic assumption in some cases.
Finally, we also evaluate how differences in the values of H result in different values of p_{V} and R_{p}. We downloaded the MPC orbital element file as of May 2012 and compared the values of absolute magnitude (H_{U}) to those used by Masiero et al. (2011), H_{M}. About 50000 Hvalues have been updated between these two works, and ~38 000 have been enlarged. Figure 1 shows a histogram of ΔH ≡ H_{U} − H_{M} for the Btypes in this work. Out of the 52 objects with ΔH ≠ 0, as many as 43 of them have ΔH > 0. Our size determinations agree to within 10%, therefore higher updated values of H will result in lower values of geometric albedos.
In Fig. B.2 we show a plot of ΔR_{p} ≡ R_{p} − (R_{p})_{M} versus ΔH for all the Btypes with determined values of R_{p}. The notation (R_{p})_{M} refers to the corresponding albedo ratios by Masiero et al. (2011). There are three features to note in this plot: (1) our values of R_{p} tend to be ~10% systematically lower, as we already noted (see Fig. B.1); (2) most points off the ΔH = 0
axis show a direct correlation between ΔR_{p} and ΔH, as expected from the discussion above; (3) some points show ΔR_{p} < −0.5, even though ΔH = 0. The points of feature (3) are explained by an inconsistency in the p_{V} values of Masiero et al. (2011) with their corresponding values of D and H: they do not verify Eq. (1) and are always lower than the predicted p_{V}.
To sum up, we have shown that if the input values of H are equal, our model fits are consistent within the model error bars with those presented in Table 1 of Masiero et al. (2011). The tendency to 10% lower values of R_{p} is likely caused by differences in solar power spectra data taken to estimate the reflected light component at NIR wavelengths (see Eq. (A.7)). We have also examined how updated input values of H affect the bestfit parameter values and showed how increasing the value of H results in greater values of R_{p} and vice versa.
All Tables
Median and mean values and standard deviations of η, R_{p}, p_{V}, and p_{IR} derived for the Btype asteroids observed by WISE.
Bestfitting values of physical parameters determined for the Stypes with sufficient W1 WISE observations and published 2.5 μm reflectances.
Mean values and standard deviations of η, R_{p}, and p_{V} derived for the members of the Pallas collisional family observed by WISE.
Bestfitting values of physical parameters determined for the Btypes with WISE observations.
Bestfitting values of physical parameters determined for the Pallas collisional family asteroids excluding (2) Pallas with WISE observations.
All Figures
Fig. 1 Differences in absolute magnitude values between those used by Masiero et al. (2011), H_{M}, and the most updated ones (as of May 2012) used in this paper, H_{U}. Note that the cases verifying ΔH = 0 do not contribute to this histogram. 

Open with DEXTER  
In the text 
Fig. 2 Histograms of value distributions obtained for the Btype asteroids observed by WISE: η (top left), p_{V} (bottom left), R_{p} (top right), and p_{IR} (bottom right). The corresponding histograms using the results from Masiero et al. (2011) are overplotted. 

Open with DEXTER  
In the text 
Fig. 3 Albedo ratio versus visible geometric albedo. The 20% error bars in p_{V} and R_{p} are not shown to facilitate visualisation. Btypes observed by WISE are plotted in black circles; all main belt objects featured in Table 1 of Masiero et al. (2011) are in empty grey circles. The horizontal line is an artificial feature corresponding to objects with the default fixed value of R_{p} = 1.5. 

Open with DEXTER  
In the text 
Fig. 4 R_{p} vs. p_{V} for the asteroids in de León et al. (2012) distinguishing the clusters to which they belong. Within parentheses, the average p_{V} of each cluster is shown. On average, p_{V} increases from G1 to G3, keeping the same value from G3 to G5. 

Open with DEXTER  
In the text 
Fig. 5 Average values of R_{p} for each cluster defined by de León et al. (2012). In terms of spectral slope up to 2.5 μm, cluster G5 is the bluest, whereas G1 is the reddest. The error bar is the standard error of the mean. 

Open with DEXTER  
In the text 
Fig. 6 Assembled spectra of the control objects chosen to study the reflectivity at 3.4 μm derived from WISE data as a diagnostic of the presence of the 3μm feature. The dashed lines approximately enclose the wavelength integration range of the W1 filter. 

Open with DEXTER  
In the text 
Fig. 7 b = 1 − R_{p}/R_{2.5} versus p_{V} for a set of primitive control objects. The values of b obtained from the combination of our R_{p} and the R_{2.5} taken from the assembly of spectra published by other authors (cf. Fig. 6) and the shape of the 3 μm absorption band are consistent except for (13) Egeria. 

Open with DEXTER  
In the text 
Fig. 8 Normalised histograms showing the distribution of bvalues. Top panel: Btypes; the fact that b > 0 in the majority of cases suggests that a high percentage of these asteroids have absorption features in the ~ 3μm region. Bottom panel: Stypes; the opposite conclusion is reached, consistent with the anhydrous nature of Stypes. 

Open with DEXTER  
In the text 
Fig. 9 Schematic diagram illustrating the interpolation of R_{2.5} (empty points) from the average values of R_{p} (filled points) for the B and S asteroids studied. 

Open with DEXTER  
In the text 
Fig. 10 Beaming parameter (left) and geometric albedo (right) distributions of the values derived for the PCF members observed by WISE. The albedo distributions are normalised, i.e. divided by their respective total number of counts. 

Open with DEXTER  
In the text 
Fig. 11 Left panel: PCF members are on average brighter in the visible and bluer in the 3.4 μm region than the Btypes sample (cf. Fig. 3). The 20% error bars are not plotted for clarity. Right panel: the average R_{p} of the PCF continues the decreasing trend of the de León et al. (2012) clusters (cf. Fig. 5). 

Open with DEXTER  
In the text 
Fig. B.1 Fractional difference histograms of D, η, p_{V}, and R_{p}. We define ε = 100(x − x_{M})/x, where x is the parameter value in this work and x_{M} the correspoding value taken from Table 1 by Masiero et al. (2011). The vertical lines mark the corresponding average values. Only parameters resulting from the same input values of H contribute to these histograms. 

Open with DEXTER  
In the text 
Fig. B.2 Differences in albedo ratio determinations versus difference in absolute magnitude corresponding to the Btypes in this paper and those by Masiero et al. (2011). 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.