A&A 381, 324-339 (2002)
DOI: 10.1051/0004-6361:20011491

Asteroids as calibration standards in the thermal infrared for space observatories[*],[*]

T. G. Müller1 - J. S. V. Lagerros2


1 - ISO Data Centre, Astrophysics Division, Space Science Department of ESA, Villafranca, PO Box 50727, 28080 Madrid, Spain
2 - Astronomiska observatoriet, Box 515, 75237 Uppsala, Sweden

Received 5 June 2001 / Accepted 19 October 2001

Abstract
Asteroids have been used extensively as calibration sources for the Infrared Space Observatory (ISO) and are planned to be used by future groundbased, airborne and space-based projects in the thermal infrared (IR) and in the sub-millimetre. We summarize the general IR observational parameters with a focus on space observatories and discuss brightness variations, apparent velocities and background influences. During the ISO mission ten well-studied asteroids were used for the photometric calibration of ISOPHOT, but additionally the bright asteroids turned out to be of great interest for many technical tests and calibration aspects. We evaluated the different applications, like testing the photometry of the spectrometers, validation of relative spectral response functions, determination of beam profiles or colour correction tests. The description of the asteroids' thermal emission has been obtained by a recent thermophysical model (TPM). The important model aspects are size, albedo, shape together with the spin vector, a beaming model, thermal inertia and a wavelength-dependent emissivity. With a large sample of observational data provided by three different ISO instruments we had for the first time the possibility to study the thermal emission of several asteroids in detail. The intercomparison between results from different instruments allowed us to distinguish between observational errors and model shortcomings. It turned out that the accuracy of TPM predictions is in many cases strongly related to the limited knowledge of the asteroid shapes. The concepts of beaming, thermal inertia and wavelength dependent emissivities were nicely confirmed for a wide range of observing and illumination geometries under many aspect angles for different asteroids. The TPM predictions for Ceres, Pallas and Vesta are accurate within 5% over the full wavelength range from 5 to 200 ${\rm\mu m}$, for Hygiea and a few other asteroids the predictions and observations agree within 10 to 15%. We found similar emissivity behaviour for the four large asteroids over the full ISO wavelength range. Up to now, no clear spectral features have been seen in the asteroid far-IR spectra.

Key words: minor planets, asteroids - radiation mechanisms: thermal - infrared: solar system


1 Introduction

With the availability of the thermal IR wavelength range (from a few micron to the sub-millimetre range) through balloon, airborne and spaceborne instruments, it became necessary to establish new calibration standards and to develop new calibration strategies. In the early phases, instruments working in these wavelengths were usually calibrated against planets (mainly Mars, Uranus and Neptune), but with the development of more sensitive instruments these objects were too bright. Sandell (1994) made a first attempt to establish a set of secondary calibrators at sub-millimetre (submm) wavelengths. The emission from stellar photospheres is too faint in the far-IR and submm and therefore all calibrators were chosen from dust-rich sources, i.e. Ultra-Compact $H {\sc ii}$ regions, protostars, protoplanetary nebulae and AGB-stars surrounded by massive dust envelopes. But often these sources are embedded in dust clouds which provide a strong and variable background (Sandell & Aspin 1998). Additionally, modeling proves to be difficult and accurate far-IR extrapolations are almost impossible. In the mid-IR, the connection from ground-based N and Q band photometry to the full mid-IR wavelength available from space was done via stellar models. In this way also the IRAS 12, 25 and 60 ${\rm\mu m}$ bands were calibrated. But neither the stellar models nor planet model extrapolations were considered reliable for the 100 ${\rm\mu m}$ band. Here, for the first time, asteroid models were used to "transport'' the calibration from 60 to 100 ${\rm\mu m}$ (Beichman et al. 1988).

The next big step in terms of sensitivity and thermal IR wavelength coverage was then achieved with the Infrared Space Observatory (ISO) (Kessler et al. 1996). ISO was equipped with four scientific instruments, covering the IR wavelengths from 2.5-240 ${\rm\mu m}$. At a wavelength of 12 ${\rm\mu m}$, ISO was one thousand times more sensitive and had one hundred times better angular resolution than its predecessor. During its routine operational phase (4 February 1996 to 8 April 1998) ISO successfully made over 26450 individual scientific observations ranging from objects in our own solar system right out to the most distant extragalactic sources, and approximately 4000 calibration observations. Especially in the far-IR beyond 50 ${\rm\mu m}$ ISO was lacking reliable photometric standards. Uranus and Neptune at the bright end and the brightest stellar calibrators at the faint end left a gap of two orders of magnitude in flux. And, a special in-flight programme to establish secondary standards was not foreseen. Müller & Lagerros (1998, in the following Paper I) provided a set of 10 asteroids, based on a recent thermophysical model code by Lagerros (1996a,1997,1998). The model input parameters were derived from ground-based mid-IR and submm observations and IRAS data. These sources have been extensively observed by ISO for photometric far-IR calibration, for testing relative spectral response functions and for many technical and calibration purposes.

In the following we summarize the important aspects of IR observational parameters and constraints for asteroids (Sect. 2). The crucial modeling issues are explained briefly (Sect. 3). In these two sections our goal is to provide useful information on the possibilities and limitations of using asteroids as calibration standards for space observatories. Section 4 gives an overview over the thermophysical model applications for ISO, including spectroscopic and photometric calibration aspects for different ISO instruments. It is meant as an illustration in what context and how asteroids can be utilized for other projects. With the ISO data it was for the first time possible to test model predictions at far-IR wavelengths in a direct way. Photometrically useful observations of asteroids from ISOPHOT (Lemke et al. 1996), SWS (de Graauw et al. 1996) and LWS (Clegg et al. 1996) are presented and discussed. We summarize these new observations and the special data reduction procedures which were necessary to calibrate the observations exclusively against stars and planets (Sect. 5). The results clearly demonstrate the excellent quality of the models, but also the possible areas of improvement of the models for future missions with even higher requirements (Sect. 6). The new observational results from ISO photometry with ISOPHOT and LWS are listed in the Appendix.

   
2 Asteroid IR observational aspects

Asteroids are point-like (<1 $\hbox{$^{\prime\prime}$ }$ angular diameter) and emit the largest fraction of their energy at thermal IR wavelengths. The spectrum of an asteroid is dominated by reflected solar radiation in the visual and near IR, and by the thermal emission at longer wavelength. The transition region is narrow due to the steep slope in the Wien part of the thermal spectrum. Outside the transition region the contribution from either the reflected light or the thermal emission is completely negligible with respect to the other.

So far, no prominent features are known in the thermal spectra of asteroids at far-IR wavelength. At shorter wavelength, between 8 and 11 ${\rm\mu m}$, broadband silicate features have been seen (Cohen et al. 1998; Dotto et al. 2000). Heras et al. (2000) report mid-IR emission signatures on Vesta, tentatively identified with olivine and pyroxene silicate groups. These features do influence ground-based N-band photometry but they are not relevant in the far-IR.

2.1 Brightness variations

The modeling of the thermal spectrum is reviewed in Paper I, but to a zeroth order the IR brightness is determined by the apparent size and the temperature across the surface. Over the orbital period (3-6yr in the main-belt) the changing distance to the Sun will induce temperature variations covering the dynamic range approximately given by $\sqrt{(1+e)/(1-e)}$, where e is the eccentricity of the orbit. From this the temperatures are typically $(17\pm13)\%$ higher in perihelion than in aphelion.

In principle the varying distance between the Earth and a main-belt asteroid (MBA), typically yields a factor of 5 between minimum and maximum flux density. In practice, however, there are constraints on the observing geometry, such as the solar elongation, which significantly reduces the dynamic range resulting from distance scaling. For ISO a factor of 2 was possible for MBAs in some extreme cases, but typical values were significantly smaller.

Asteroids rotate which produces a light curve, mainly due to the elongated shape and in some cases albedo variations across the surface. For the asteroids larger than 200km in diameter, the range of periods is 4.2- $29\,{\rm h}$ with a mean of $8.7\,{\rm h}$, and the mean amplitude of the visual light curve $0.22\pm0.15$ mag (Lagerkvist, private communication). The amplitude in the thermal IR is comparable to that in the visual, since both are approximately proportional to the projected surface area. Depending on the target and the application it may be important to consider the flux changes due to rotation, as given by the TPM, for integration times longer than about $0.5\,{\rm h}$.

2.2 Apparent velocities

Typical MBAs had between 0 and about 80arcsec/hour apparent speed, as seen from the Earth (vector addition of Earth and MBA movement). The ISO motion added in the worst case about 20arcsec/hour at certain points of the 24 h orbit (see Fig. 1). Near Earth objects can have up to a few degrees/hour for close encounters with Earth.

  \begin{figure}
\par\includegraphics[angle=90,width=8.8cm,clip]{MS1543f01.eps} \end{figure} Figure 1: The apparent sky movement [ $^{\prime \prime }$/hour] of the asteroid Vesta with respect to the geocentre and ISO (24h periodicity). ISO visibility periods are for $\lambda {-}\lambda _\odot $ between 60 and 120$^{\circ }$, corresponding to days 151 to 237 and 343 to 418 in the time scale of the graphics.
Open with DEXTER

Observations of solar system objects require therefore accurate, and in some cases, fast tracking. In case of ISO, tracking was done as a spacecraft raster with 2arcsec step size, limited to apparent velocities of less then than 120arcsec/h with respect to equatorial coordinates. This excluded fast moving objects like a few near Earth asteroids and comet Hyakutake from the ISO programme. The relative velocity between typical MBAs and Earth varies between 10 and 50 kms-1. The ISO orbit added an amplitude of about 10 kms-1 with a 24h periodicity. The apparent sky motion of solar system objects has advantages and disadvantages, depending on the purpose of a measurement (see Sect. 4.4). In general, tracked satellite observations are usually more difficult than similar ones from the ground due to the involvement of gyros and reaction wheels.

   
2.3 Infrared background

The IR background in the ecliptic plane as seen from space is dominated by zodiacal light and therefore relatively smooth. It can be approximated as diluted black body radiation with temperatures ranging from 259K (at $\lambda-\lambda_\odot=120^\circ$) to 280K (at $\lambda-\lambda_\odot=60^\circ$) (Abraham et al. 1999). Only at around $\lambda _{\rm ecl.}=90^{\circ }$ and $\lambda _{\rm ecl.}=270^{\circ }$, where the ecliptic crosses the galactic plane, the background strongly increases (see Fig. 2).

  \begin{figure}
\par\includegraphics[angle=90,width=13cm,clip]{MS1543f02.eps} \end{figure} Figure 2: The sky background in the ecliptic plane for a solar elongation of $\lambda _{\rm ecl.}-\lambda _\odot =90^\circ $ at IR wavelengths. At $\lambda _{\rm ecl.}=90^{\circ }$ and $\lambda _{\rm ecl.}=270^{\circ }$ the ecliptic crosses the galactic plane. Note: the values are given in logarithmic scale.
Open with DEXTER

At 25 ${\rm\mu m}$ the background along the ecliptic plane varies between 60 and 100MJy/sr (as given by the COBE/DIRBE annual maps, corresponding to $\lambda-\lambda_\odot=90^\circ$, Hauser et al. 1998) at 140 ${\rm\mu m}$ the background values are around 10MJy/sr, but reaching up to 60MJy/sr at $\lambda=90^{\circ}$ and more than 1500MJy/sr at $\lambda=270^{\circ}$ (close to the galactic centre). In these places of high background the interstellar cirrus dominates at wavelengths longer than 100 ${\rm\mu m}$.

The background of a given position in the ecliptic plane changes with time, depending on the angular distance from the Sun. Under the ISO visibility constraints, between $\lambda-\lambda_\odot=60^\circ$ and $\lambda-\lambda_\odot=120^\circ$, the background varied up to a factor of 3 at the mid-IR and about 10 to 30% at the far-IR. Preferable, low background configurations are therefore between $\lambda_{\rm ecl.}=0...70^{\circ}$, 100...250$^{\circ }$ and at $\lambda_{\rm ecl.}=280...360^{\circ}$ for $\lambda-\lambda_\odot > 90^\circ$. Background values of 10, 100 and 1000MJy/sr correspond to approximately 2, 20 and 200Jy in a 100 $^{\prime \prime }$ circular aperture and to 0.02, 0.2 and 2Jy in a 10 $^{\prime \prime }$ aperture. Although sky backgrounds can be high in some cases they only play a role for "cold" space observatories. Groundbased, airborne and passively cooled space instruments are limited by atmospheric and instrument self-emission. Sources are therefore measured in chopped mode and the sky background is negligible.

   
3 Thermophysical modeling

Using asteroids as calibrators requires a good understanding of the temporal flux variations on orbital and rotational time scales, as discussed above. Paper I and Müller et al. (1999) applied the thermophysical model developed by Lagerros (1996a,1997,1998), to a large data base of IR observations of a selection of previously well studied asteroids, in order to establish the asteroid standards used for ISOPHOT (Lemke et al. 1996). The model produces accurate thermal IR spectra and thermal lightcurves, taking into account a number of physical processes. Subjectively ordered by their relative importance, the main components of the model are:

The fundamental parameter for far-IR model predictions is the object cross section at the time of the observation. Therefore, the size and shape and spin state uncertainties clearly dominate the model uncertainties. Tied to the size is the albedo because it determines the energy balance between reflected and absorbed solar radiation. In the model the emissivity was fitted to mid-IR and submm data and interpolated in the unexplored far-IR. Theoretically the beaming effect is much stronger in the mid-IR. For example, if the effect is constrained to a 10%-level by data at $10\,{\rm\mu m}$ it is typically within a 3-5%-level at $50\,{\rm\mu m}$, and 1-3%-level at $200\,{\rm\mu m}$. The thermal inertia causes the morning side to appear colder than the afternoon side, and some of the emission to occur on the night side (Morrison 1977; Hansen 1977; Lebofsky et al. 1986). Increasing the inertia decreases the amplitude of the thermal light curve. Varying the inertia introduces a significantly larger effect at shorter wavelengths, which are more sensitive to changes in surface temperature. Individual values of TPM input parameters are given in Paper I for 10 asteroids. TPM spectral energy distributions and lightcurve predictions at any thermal wavelength can be requested from the authors via electronic mail. But reliable predictions are only possible for cases where the sizes and shapes are well known.

   
4 TPM applications for ISO

The asteroids were utilized in a wide range of calibration aspects ranging from absolute photometry, relative spectral response function checks to cross calibration and calibration of special modes. See also Müller & Lagerros (2001) for a more technically oriented overview of all asteroid related calibration of ISO. The brief description of the main calibration aspects also illustrate the potential use of asteroids for future far-IR and submillimetre projects.

   
4.1 Photometric calibrations of ISOPHOT in the far-IR

A set of 10 asteroids was used during the ISO mission, together with bright standard stars and the planets Uranus and Neptune, to establish the photometric calibration of ISOPHOT at far-IR wavelength beyond 45 ${\rm\mu m}$ (Schulz et al. 1999). Pallas was used also at shorter wavelength, down to 12 ${\rm\mu m}$.


  \begin{figure}
\par\includegraphics[angle=90,width=12.8cm,clip]{MS1543f03.eps} \end{figure} Figure 3: The ISOPHOT calibration scheme. The stars cover the lowest flux range, the asteroids the intermediate one and the planets the high flux range. Individual measurements of celestial standards were used to calibrate the corresponding flux output of the internal calibration source (FCS), expressed in heating power.
Open with DEXTER

Figure 3 shows an example of the photometric calibration of ISOPHOT's internal reference source (FCS) against celestial standards for one pixel of the C200 camera in the 120 ${\rm\mu m}$ filter. The signals of the FCS (expressed in terms of heating power) were ratioed to the signals produced by celestial standards, where the fluxes (or inband powers) were predicted from models. The model predictions for the asteroids were provided through Paper I. They established the calibration of the ISOPHOT FCS at intermediate flux level, whereas the stars served in the low range and the planets at the bright end of the detectors. In this way an absolute photometric calibration accuracy of 10% (P3), 15% (C100) and 10% (C200) was achieved for point sources close to the brightness level of the 10 asteroids (Klaas et al. 2000). Taking the general difficulties of opening a so far unexplored wavelength region beyond 100 ${\rm\mu m}$ into account and considering the asteroid shape and thermal modeling aspects, the concept of using asteroids for the absolute calibration worked extremely well. Although the always changing background conditions for the asteroids required special attention and many dedicated measurements.

   
4.2 LWS photometry and spectral response function

Both the LWS absolute flux calibration and the relative spectral response function (RSRF) have been established using observations of Uranus in combination with a model (Griffin & Orton 1993).

The photometry of LWS has been tested against Ceres, Pallas, Vesta and Hygiea repeatedly. The asteroids were considered as a viable possibility to ensure the calibration, since the brightest stellar standards are already too faint to obtain sufficient S/N values. The bright asteroids provide approximately one tenth of the Uranus flux level (beyond 100 ${\rm\mu m}$), whereas $\alpha$Boo and $\alpha$Tau are only at one hundredth of the Uranus brightness. The better flux coverage revealed non-linearities in the detector responses and allowed their empirical correction (Gry et al. 2001). Background contributions are not important for the bright planets, but had to be considered for the asteroids.

No prominent lines are known for the asteroids in the LWS wavelengths range. They therefore provide means of an independent RSRF determination. Photometric and spectroscopic asteroid measurements are currently analysed and the findings will be implemented in the LWS calibration scheme.

   
4.3 SWS photometric and spectroscopic calibration

Most of the SWS calibration programme were done on stars which perfectly cover the dynamic range of the detectors, except the highest brightness level beyond 12 ${\rm\mu m}$. Here, bright asteroids and planets were used to check the final photometric and spectroscopic calibration (Morris 1999).

The grating flux calibration programme included Ceres, Pallas, Juno, Vesta, Hygiea together with the planets Saturn, Uranus and Neptune. For the Fabry-Perot interferometer calibration only measurements on Ceres and Pallas were used in addition to the stars.

Further qualitative comparisons related to spectral features or broadband shapes on the grating relative spectral response function included measurements on Ceres, Pallas, Juno and Hygiea. Also for the validation of different instrument configurations (SWS01/SWS06) a few asteroids were used (Leech et al. 2001).

   
4.4 Miscellaneous applications and results

Tracking and beam profiles:

the apparent sky movements of three asteroids were used to validate the pointing-tracking system of ISO. Additionally, the apparent sky movement of Vesta served in a special calibration programme to characterize the beam profile in dispersion direction on subarcsecond level. This complemented the extremely time consuming beam profile measurements (typically done with a series of satellite repointings on bright stars). In this way also the pointing uncertainties of the telescope could be minimized.

Colour correction and filter leaks:

the accuracy of the photometric calibration depends strongly on the quality of the blocking elements and the characterization of the band passes. Ceres (in December 1997) and HR6705 ($\gamma$ Dra) had both approximately 150Jy flux density at 10 ${\rm\mu m}$. In case of the asteroid, 10 ${\rm\mu m}$ is still on the Wien-part of the spectrum with its energy peak at around 15 ${\rm\mu m}$. For HR6705 10 ${\rm\mu m}$ lies on the Rayleigh-Jeans part of the spectrum with its peak wavelength in the visible. The measured in-band fluxes are therefore different and depend strongly on the width and shape of the band pass. These differences are clearly visible in the colour correction factors, which are calculated for a wide range of sources in comparison to a reference spectrum (IRAS and ISO used a constant energy spectrum $\nu F_{\nu} = \rm const.$ as reference).


 

 
Table 1: Colour correction factors for blue ( $T_{\rm BB}=5000$K) and red sources ( $T_{\rm BB}=200$K) for 5 selected ISOPHOT filters.

Detector-Filter
$T_{\rm BB}$ P1_3.6 P1_12 P2_25 P3_60 C100_120

5000K
1.05 1.30 1.28 1.13 1.21
200K 1.69 0.89 1.05 1.07 1.17


Table 1 gives an example of colour correction values between typical stellar and asteroid temperatures for different ISOPHOT filters. Indications for filter leaks are given when colour corrected monochromatic fluxes of one source type agree nicely with model predictions, but for the "opposite'' source type they do not agree. ISOPHOT and also ISOCAM used asteroids for the blue versus red sources tests, to validate the colour correction factors and to investigate filter leaks.
ISOPHOT Serendipity Mode:

Whenever possible, the long-wavelength channel of the photometer was used during satellite slews. This ISOPHOT Serendipity Mode (e.g. Stickel et al. 1999) led to a partial sky survey, covering approximately 15% of the sky, at wavelengths of 170 ${\rm\mu m}$, a spectral region not covered by the IRAS survey.

The bright asteroids, seen serendipitously during the ISO slews, allowed a flux calibration in a direct way, using TPM predictions (Müller et al., in preparation). This procedure improved the calibration significantly and extended the calibrated flux range to higher brightness levels, where many new galactic and extragalactic sources have been found (Stickel et al. 2000).

   
5 New ISO observations of asteroids

ISO used asteroids for calibration purposes in many contexts, but it also allowed an independent testing of TPM predictions in the unexplored far-IR. To avoid "self-calibration" of the asteroid observations with a calibration based on asteroids, we developed special procedures to establish "asteroid free" calibration for the following data samples.

5.1 Independent ISOPHOT observations

With the method described in Paper I it was possible to exclude all asteroid observations from the ISOPHOT calibration scheme and to extract independently calibrated asteroid fluxes via the internal calibration source. Asteroid observations were in this way treated as normal scientific measurements, but now calibrated only against stars and planets. In total we derived 94 individual photometric data points between 50 and $200\,{\rm\mu m}$ (Tables B.1-B.3). All values are background subtracted and colour corrected. The uncertainties are rms values of observational errors, calibration source model uncertainties and estimated method errors. Additionally, the TPM predictions are given in the tables. There are many more ISOPHOT observations of asteroids where our method could not be applied because of a lack of stellar and planetary calibrators in the specific instrument configuration. A detailed discussion per object follows in Sect. 6.

5.2 LWS AOT observations of asteroids

5.2.1 Full grating scans (LWS01)

The ISO Data Archive contains two full grating scans (observing mode LWS01) on Ceres and one on Hygiea (see Table 2). Additional grating scans exist on 4 asteroids, but all taken in the photometrically unreliable LWS99 mode.


 

 
Table 2: Summary of LWS01 observations on asteroids.

TDT
JD (mid time) Object Scans

74803403
2450785.67361 Ceres 22 up
76903203 2450806.59514 Ceres 12 up & 12 down
83201702 2450869.41667 Hygiea 10 up & 10 down


The observations (OLP10 products) have been analysed with ISAP2.0[*]. An equal number of up and down scans for each detector were deglitched, sigma clipped (2.5$\times$rms) and averaged. Detectors SW1 to SW5 were then smoothed with the nominal resolution element of 0.29 ${\rm\mu m}$, LW1 to LW5 with a resolution element of 0.6 ${\rm\mu m}$ and both scan directions averaged. Observations on a large number of standard sources (planets, asteroids, stars) showed that the responsivities of the detectors LW1, LW2 and LW3 are flux dependent. Empirical correction values have been derived and were applied to the asteroid scans. The background values have been taken from COBE-DIRBE (Hauser et al. 1998) weekly maps (25-100 ${\rm\mu m}$) and yearly maps (140-240 ${\rm\mu m}$), together with the LWS solid angles and the correction factors for extended sources. Maximal background contributions were 4% for Ceres and up to 50% for Hygiea with the largest influence at long wavelengths. One Ceres observation (74803403) was affected by straylight, which is noticeable in the noise of the long wavelength detector signals (Fig. 4).

The jumps are produced by the individual detector calibration. Additional artifacts of the RSRF can be seen at the edges of the detector bands. At longer wavelength, where the source signals are already comparable with the dark current, the noise level is much higher, but the overall flux levels agree still quite well with the model predictions (the LW4 and LW5 are considered not to be photometrically reliable at this flux range). See also Sect. 6 for quantitative analysis of the observation to model ratios.


  \begin{figure}
\par\includegraphics[angle=90,width=8.8cm,clip]{MS1543f04.eps} \end{figure} Figure 4: Calibrated and background subtracted LWS spectrum of Ceres (3 Dec. 1997, revolution 748). Solid line: TPM prediction.
Open with DEXTER

5.2.2 Fixed grating positions (LWS02)

Ceres, Pallas, Vesta and Hygiea have been observed in the L02 narrow band mode, where 10 wavelengths were taken at the same time - one in each channel - with the grating remaining at a fixed position.

The data reduction was done manually by taking only the last 100s of the measurements where the signals were usually stabilized. Dark current subtraction and flux calibration was done in the ISAP command mode. The detector non-linearities and the backgrounds have been corrected in the same way as for the LWS01 mode. Maximal background contributions were up to 16% for Ceres, up to about 10% for Pallas and Vesta and up to 40% for Hygiea. The dark current signals in LW4 (160.6 ${\rm\mu m}$) and LW5 (178.0 ${\rm\mu m}$) are in many cases unstable during the integration times and the corresponding fluxes could therefore not be used. Four Ceres and one Pallas L02 observations took place before ISO revolution 236 when the illuminator sequence for the instrument calibration was different. In these cases the calibration of the SW1 detector seems to be systematically off by 20 to 25% and fluxes are too low. No empirical correction has been done so far, but based on the asteroid results, this might be possible in the future. Tables B.1-B.4 summarize the observational results together with the model predictions. See also Sect. 6 for further discussions and the implications of these findings.

  \begin{figure}
\par\includegraphics[angle=90,width=8.8cm,clip]{MS1543f05.eps} \end{figure} Figure 5: Calibrated and background subtracted LWS spectrum of Ceres (24 Dec. 1997, revolution 769). Solid line: TPM prediction.
Open with DEXTER

5.3 SWS observations of 1 Ceres and 10 Hygiea

Two sets of measurements (see Table 3) of Ceres and Hygiea have been analysed in a standard way, including the latest tracking/beam profile corrections (OLP10 products). No emphasis was put on the verification of spectral features, only the overall flux levels per band were of interest. Background values are usually below the 1-2% level for the bright asteroids (see also Sect. 2.3) and were therefore neglected.
 

 
Table 3: Summary of 2 sets of SWS observations (SWS06) on Ceres and Hygiea. The nominal wavelength ranges of the SWS bands are [in ${\rm\mu m}$]: 1A: 2.38-2.60; 1B: 2.60-3.02; 1D: 3.02-3.52; 1E: 3.52-4.08; 2A: 4.08-5.30; 2B: 5.30-7.00; 2C: 7.0-12.0; 3A: 12.0-16.5; 3C: 16.5-19.5; 3D: 19.5-27.5; 3E: 27.5-29.0; 4: 29.0-45.2.

TDT
JD (mid time) Object Bands

71401706
2450751.78889 Ceres 2C, 3E, 4
71401805 2450751.82361 Ceres 1D, 1E, 2A, 2B, 3C, 3D
71401904 2450751.85486 Ceres 1A, 1B, 3A

67000903
2450707.67153 Hygiea 1B, 3A
67000904 2450707.72500 Hygiea 1E, 2A, 2B, 3C, 3D
67000905 2450707.77639 Hygiea 2C, 3E, 4


Figure 7 shows the observational results together with the model predictions on absolute scales. At short wavelengths the transition region between reflected light and thermal emission can be seen. A discussion on the quality of the TPM predictions per object is given in Sect. 6.

   
6 Discussion

6.1 Asteroids as technical calibrators for spacebased IR projects

The brightness level of the point-like asteroids, together with the high accuracy of the model predictions allow their usage as technical calibrators for space observatories in many aspects where high S/N is needed. Even their apparent sky movement can be employed for the characterization of point-spread function, beam and aperture profile scans without moving the satellite. Satellite tracking in combination with instrument apertures significantly smaller than the point-spread function can cause severe flux changes depending on the source position within the aperture. In some ISO configurations tracking was done with apertures smaller than 20 $^{\prime \prime }$ diameter which produced repeated flux changes for each tracking step.

  \begin{figure}
\par\includegraphics[angle=90,width=8.8cm,clip]{MS1543f06.eps} \end{figure} Figure 6: Calibrated and background subtracted LWS spectrum of Hygiea (24 Feb. 1998, revolution 832). Solid line: TPM prediction.
Open with DEXTER

The spectral energy distribution of main-belt asteroids with energy peaks at around 15 ${\rm\mu m}$ makes them perfect counterparts to stellar calibrators, ideal to perform filter leak tests and to validate colour correction factors in the mid-IR. The high ephemeris precision needed for modern state-of-the-art instruments is not an obstacle anymore. Recent N-body ephemeris software even allow satellite centric position calculations for solar system objects, including parallax corrections due to the satellite orbit. For the well known bright main-belt asteroids ephemeris accuracies of 1arcsec are reached.
  \begin{figure}
\par\includegraphics[angle=90,width=13cm,clip]{MS1543f07.eps}
\par\includegraphics[angle=90,width=13cm,clip]{MS1543f08.eps} \end{figure} Figure 7: SWS observations of Ceres and Hygiea. The TPM predictions are plotted as solid lines.
Open with DEXTER

The always changing sky background conditions are clearly a disadvantage of asteroids. The high background areas, especially where the ecliptic and the galactic plane intersect, should be avoided. Depending on the size of the instrument aperture, far-IR space observations require investigations of the actual background conditions also outside the galactic plane. All asteroids presented here are known to have only small amplitudes in their visual light curve (less than 0.23mag in the most extreme case, Lagerkvist et al. 1989). Depending on albedo, thermal parameters and observing geometry the thermal light curves have even smaller amplitudes. But long calibration observations of several hours might be affected by brightness changes of a few percent. But with the TPM it is possible to predict the light curve variations at any wavelength, based on shape and rotational aspects, and to take it into account for calibration observations. The lack of prominent lines in the asteroid spectra make them also useful calibrators for the RSRF determination. Especially at far-IR and submillimetre wavelengths where the behaviour of stellar photospheres and planetary atmospheres is not perfectly well known, asteroids provide independent means to verify the spectroscopic calibration.

6.2 The quality of the asteroid models

6.2.1 Observation to model ratios

The independently derived ISO data allowed us to test the TPM quality. The observations have been brought together by calculating the corresponding TPM predictions and plotting the observation to model ratios.

  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f09.eps} \end{figure} Figure 8: Ceres SWS observation divided by the TPM prediction. Jumps between SWS bands and the slope of band 4 are calibration artifacts. At around 10 ${\rm\mu m}$ broad band Silicate features are visible.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f10.eps} \end{figure} Figure 9: Ceres LWS (fixed grating: squares; grating scans: dots) and ISOPHOT (triangles) observations divided by the corresponding TPM predictions. The LWS LW4 detector scan is off by about 25% for unknown reasons. No systematic differences between LWS and ISOPHOT can be seen.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f11.eps} \end{figure} Figure 10: Pallas LWS (fixed grating: squares) and ISOPHOT (triangles) observations divided by the corresponding TPM predictions. No systematic differences between LWS and ISOPHOT can be seen.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f12.eps} \end{figure} Figure 11: Vesta LWS (fixed grating: squares) and ISOPHOT (triangles) observations divided by the corresponding TPM predictions. No systematic differences between LWS and ISOPHOT can be seen.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f13.eps} \end{figure} Figure 12: Hygiea SWS observation divided by the TPM prediction. Jumps between SWS bands and the slope of band 4 are calibration artifacts.
Open with DEXTER


  \begin{figure}
\par\includegraphics[angle=90,width=12.2cm,clip]{MS1543f14.eps} \end{figure} Figure 13: Hygiea LWS (fixed grating: squares; grating scans: dots) and ISOPHOT (triangles) observations divided by the corresponding TPM predictions. Jumps between LWS bands are calibration artifacts. Especially the high quality ISOPHOT C200 measurements are consistent with the TPM predictions.
Open with DEXTER

The observation to model ratio representation (Figs. 8-13) include a large variety of observing geometries, aspect and phase angles, typically 2-3 orders of magnitude in flux and different background conditions for each asteroid. The observation from the 3 instruments (SWS, LWS, ISOPHOT) are taken in different observing modes (SWS06, LWS01, LWS02, PHT22, PHT25, PHT99) and different integration times, ranging from a few seconds (PHT) to hours (SWS). They underwent completely different data processing and calibration schemes. Nevertheless, the observation to model ratios give a consistent picture, expressed in ratios close to one. Deviations from the model are caused by different effects. In case of the ISOPHOT, the uncertainties in the background determination and in the calibration dominate the final errors. For SWS and LWS the background contributions are of minor importance. Instead problems at the band edges appear. The SWS band 4 in addition has a wrong slope and LWS detectors LW4 and LW5 are affected by dark current problems. Some deviations are due to model uncertainties, mainly by shape and size approximations. Systematical offsets in Figs. 8-13 indicate incorrect sizes (for example Juno), standard deviations (Tables 4, 5) larger than typical observational errors indicate problems with the shape approximation by ellipsoids (for example Hygiea). Within the error bars the observations agree extremely well with the model predictions over all wavelengths and at all flux levels. Clearly visible features (up to 5% over the continuum) exist only at around 10 ${\rm\mu m}$ (see also Dotto et al. 2000). More detailed spectroscopic analysis in the future might reveal additional structures at levels below 5%.

   
6.2.2 Results per instrument
Before the quality of individual asteroid models are discussed, we analysed the results per instrument to see if different calibration concepts lead to the same results.
ISOPHOT:
the statistical analysis of the ratios between observations and model predictions, are given in Table 4.
 

 
Table 4: Averaged ratios between ISOPHOT observations and model predictions. The total number of observations varies between 3 and 16 per object. All measurements are taken between 60 ${\rm\mu m}$ and 200 ${\rm\mu m}$.

Asteroid
No of Meas. Mean Ratio Std.Dev.

1 Ceres
16 1.03 0.08
2 Pallas 13 0.97 0.08
3 Juno 11 1.13 0.10
4 Vesta 12 0.97 0.11
10 Hygiea 16 1.00 0.16
54 Alexandra 3 0.98 0.11
65 Cybele 13 1.04 0.10
532 Herculina 10 0.93 0.07


In general observations and model predictions agree on absolute level very well over the whole ISOPHOT wavelength range. In case of Juno the model predictions seem to underestimate and in the case of Herculina to overestimate the observed fluxes. The difficulty to calibrate these measurements with stars and planets only (see Sect. 5) is reflected in the large scatter of the ratios. Additionally, the uncertainties in the shape model might contribute to the scatter for some objects.

LWS:
the ratios between LWS01 observations of Ceres and the corresponding TPM predictions per detector range between 0.94 and 1.04, with a weighted mean of $0.99\pm0.03$. For Hygiea the individual ratios vary between 0.83 and 1.22, with a weighted mean of $1.04\pm0.15$ (Figs. 4-69 and 13). Here, as in case of the SWS measurements, the band average are the relevant numbers for the comparison. The band slope and band edge effects are of minor interest for the model tests. The photometry of the LW4 and LW5 detectors was not reliable due to dark current problems.
 

 
Table 5: Ratios between LWS02 observations and TPM predictions, weighted mean values and standard deviations per detector (SW1, 2, 3, 4, 5 and LW1, 2, 3) and per asteroid. The photometry of the LW4 and LW5 detectors was not reliable due to dark current problems.

Ceres Pallas Vesta Hygiea
$\lambda$ 15 Obs. 4 Obs. 14 Obs. 2 Obs.
[$\mu$m] Obs/Mod Obs/Mod Obs/Mod Obs/Mod

46.2
$0.94\pm0.11$ $0.91\pm0.12$ $0.99\pm0.05$ $0.86\pm0.03$
56.2 $0.97\pm0.03$ $0.98\pm0.05$ $1.01\pm0.07$ $1.09\pm0.16$
66.1 $1.01\pm0.02$ $1.04\pm0.04$ $1.02\pm0.04$ $1.23\pm0.11$
75.7 $0.99\pm0.02$ $0.94\pm0.03$ $1.01\pm0.06$ $0.98\pm0.08$
84.8 $1.00\pm0.04$ $0.99\pm0.11$ $0.94\pm0.12$ $1.22\pm0.07$
102.4 $1.00\pm0.04$ $0.98\pm0.03$ $1.02\pm0.06$ $1.08\pm0.05$
122.2 $1.01\pm0.05$ $1.01\pm0.05$ $1.05\pm0.11$ $1.10\pm0.07$
141.8 $0.99\pm0.06$ $1.04\pm0.09$ $1.04\pm0.22$ $1.05\pm0.23$

Mean
$1.00\pm0.02$ $0.98\pm0.04$ $1.01\pm0.03$ $0.98\pm0.11$


The general agreement between the photometric observations (L02) and models (Table 5) is better than 5% for Ceres, Pallas and Vesta and better than 15% for Hygiea. The two fixed grating measurements (squares in Fig. 13) follow the same jumps as the individually scanned bands (dots). We believe that this is due to dark current effects at the low flux level for Hygiea.

SWS:

the ratios per detector between SWS observations of Ceres and TPM predictions range between 1.01 and 1.05, with a weighted mean of $1.02\pm0.01$. For Hygiea the ratios vary between 0.85 and 1.18, with an weighted mean of $1.03\pm0.10$ (Figs. 7, 8 and 12). SWS band 4 is photometrically difficult due to the deviating band slope, which seems to occur for all solar system objects. The Hygiea band 2b is affected by the low level flux.

No systematic offsets on absolute terms between the instruments or observing modes can be seen. The jumps between SWS and LWS detector bands are instrumental artifacts which are not yet perfectly understood, but photometrically reliable are in general only the values at the key wavelengths of the bands.

6.2.3 Model quality per asteroid

Not all individual asteroid models have the same quality. Already in Paper I the 10 asteroids were divided in primary and secondary calibrators, depending on the quality and coverage of the available observational data and on the reliability of the size/shape values. A similar separation can also be seen in Figs. 8-13 with observations from ISOPHOT, SWS and LWS.

1 Ceres:

based on a large set of ISO observations from 3 different instruments the model predictions agree within 3% at wavelength between 5 and 200 ${\rm\mu m}$ (Tables 4, 5 and Figs. 8, 9). The larger scatter in the ISOPHOT data is related to the calibration procedure. Submm and millimetre data proved that the Ceres TPM works also out to about 2mm (Paper I). Ceres is therefore the best modeled asteroid and an excellent far-IR calibration standard. The TPM concept (Lagerros 1996a,1997,1998) and the model input parameters (Müller et al. 1999) are confirmed. Caution is only necessary for N-band photometry (see SWS band 2C in Fig. 8) where the spectral behaviour from 8 to 11 ${\rm\mu m}$ suggests the presence of silicates on the surface (Dotto et al. 2000). Detailed spectroscopic analysis might reveal other low level spectral features at longer wavelengths.

2 Pallas:

the 13 ISOPHOT and the 4 LWS fixed grating observations are within 3% of the predictions. Individual measurements with large deviations can be tracked down in most cases to difficulties in the calibration (see for example LWS SW1 problem mentioned in Sect. 5). The shape values of Pallas are uncertain and our parameters deviate from the occultation and speckle results (Paper I). The observation to model ratios close to one seem to confirm our parameters, but some scatter might be due to deviations from the real physical shape of Pallas. For this asteroid a 3-dimensional shape from time delay and Doppler frequency radar measurements (Ostro 1993) would be of great importance.

4 Vesta:

the 12 ISOPHOT and the 14 LWS grating observations are within 3% of the predictions. The scatter in Fig. 11 is large and a few individual observations are off by 25%. This might be an indication that the physical shape model could be improved, although the highest quality 3 dimensional HST shape (Stooke 1997) has been used. The albedo variations on the surface (Binzel et al. 1997) have not been taken into account (Müller & Lagerros 1998) and might cause deviations in one or the other direction, depending on the rotational phase. The ISOPHOT data points beyond 150 ${\rm\mu m}$ are systematically below 1.0. This indicates that the emissivity drop to 0.6 in the submm (Redman et al. 1998; Redman et al. 1992) might start already in the far-IR.

10 Hygiea:

all data sets agree on average within a few percent with the model predictions. Nevertheless, the model predictions are only accurate to about 15% for a given time. The reason lies mainly in the uncertainties of the shape model. Based on lightcurve observations the axis ratio of the ellipsoidal shape in Paper I was assumed to be a/b=1.29 and b/c=1.18. Erikson (2000) stated now that a meaningful determination of b/c was not possible. Based on IRAS data the effective diameter would be 407.1 km (Tedesco et al. 1992), we found in Paper I 429.9 km. Our shape model produced at the time of the SWS observation an apparent effective diameter of 448km and in consequence a too high model flux. Forcing the apparent effective diameter to our previous value of 429.9 km gave much better agreement with the SWS observations (Fig. 7), except for band 4 where the calibration is problematic (Barucci et al. 2001). Here again, a radar based shape model could improve the situation significantly.

Other asteroids:

for the asteroids Juno, Alexandra, Cybele and Herculina we have only ISOPHOT far-IR observations to test the model quality. Additionally, their diameters are based on radiometric methods only (Paper I). Radar observations and direct imaging (HST, Speckle) would improve the size and shape models and put the models on more solid ground. For Juno the eleven ISOPHOT measurements show that the model systematically underestimated the derived fluxes by more than 10%. Since the direct diameter measurements and the radiometric diameters agree almost perfectly (Paper I), the source of the deviations is not known. Independent fluxes from the submm would allow to test our model values, especially the emissivity behaviour and the radiometric diameter. Erikson (2000) published an ellipsoidal shape model which differs from the one we have used. Additional lightcurve, thermal and/or radar observations could also reveal the cause of the deviations. The 3 ISOPHOT observations of Alexandra agree nicely, but a direct size measurement and an improved shape model would lower the model uncertainties. The Cybele model agrees on average within 4% with all 13 ISOPHOT measurements, but has a 10% uncertainty for individual measurements. The radiometric size (Paper I) and the shape and spin vector solution (Erikson 2000) are not very reliable and could be improved by additional measurements. Using the 9% lower IRAS diameter in the TPM would produce an average ratio of 0.85 for the 13 ISO observations. Here again, an improved shape from other observing techniques would lower the uncertainties. Herculina's shape was derived from lightcurve observations (Michalowski 1996). The combination of occultation measurements (Bowell et al. 1978), speckle observations (Drummond et al. 1985) and radiometric analysis (Paper I) agreed resonably well. Nevertheless, the model predictions are systematically too high. Our previous classification as primary calibrator cannot be kept anymore and individual model predictions can be off by 15% in some extreme cases.

In summary, we classify now only Ceres, Pallas and Vesta as primary calibrators with an estimated model uncertainty of 5%, although individual predictions can be off by as much as 10%. All other object models have larger uncertainties. We estimate the TPM prediction accuracy to about 10 to 15% with the possibility that individual calculations might be off by as much as 20%. The reliable TPM wavelengths range can be given with 5 to 200 ${\rm\mu m}$ for well known main-belt asteroids. At the submm and millimetre range the TPM can be used, but with lower accuracy (Paper I).

6.3 Model input parameters

Size/shape:
the sizes and shapes of the asteroids determine the absolute flux level. Since the observation to model ratios are all close to 1.0 we conclude that the specified asteroid sizes are close to the true values. Nevertheless, some of the diameters need confirmation from direct size methods, like imaging, speckle, occulation or radar techniques. The object shapes, and therefore also the cross sections at a given time, are more complex. Approximating the shape as a 3-dimensional ellipsoid reproduces the visible and thermal light curves to a large extent, but some deviating model predictions might be due to the simplification of the shapes in combination with uncertain ellipsoid ratios of the semi-major axis.
Albedo variations:
albedo changes on the asteroid surface of $\pm 10$% between different hemispheres would only produce far-IR flux changes of about $\pm$1%. Thus, even from high quality taken over a full rotation period it would be difficult to determine these variations. Rotationally resolved mid-IR observations would be better suited. Although the influence of albedo variations in the far-IR is of minor importance, in case of Vesta it would improve the model quality.
Emissivity:
in Paper I the wavelength dependent emissivity was derived from mid-IR and submm observations. All asteroids show a similar behaviour in the mid-IR, with emissivities between 0.9 and 1.0. At longer wavelength a clear trend of decreasing emissivity can be seen, reaching values of 0.83 at 200 ${\rm\mu m}$ and 0.8 in the submm. In case of Vesta, submm emissivities as low as 0.6 are found (Redman et al. 1992,1998). But due to a lack of observations it was not clear where and in what way the emissivity changes in the far-IR. To have better possibilities to investigate the emissivity behaviour in the unknown far-IR range we used exactly the same emissivity model for all asteroids. And in fact it is difficult to distinguish Vesta from the others by looking at the ratio figures in the ISO far-IR wavelength range. There is a trend in Figs. 10 and 11 with lower emissivities for Pallas and Vesta beyond 130 ${\rm\mu m}$, but more observations are needed to clarify that aspect. We conclude therefore that Vesta's emissivity has to drop within a short wavelength range closer to the submm from 0.8 to 0.6, possibly related to grain size effects and/or Vesta's composition. Laboratory measurements of different possible asteroid surface materials are currently underway (Brucato, priv. communications). This might be an additional possibility to investigate the emissivity behaviour in the thermal-IR.
Light curve variations:
the thermal light curve variations are predicted by the TPM based on shape models and the thermal behaviour. But only limited observational material is available to test the model predictions. The here presented ISO observations contain also 13 consecutive LWS L02 measurements over a full rotation period (5.34hours) of Vesta. The TPM predicts a light curve amplitude of 3.8% for this epoch at the LWS wavelengths range. The normed mean values of all 10 LWS detectors varied at the same time between $1.01\pm0.03$. But distinguishing LWS observational uncertainties from effects due to the asteroid rotation is difficult at this low level. The same is true for the analysis of the light curve phase: the large scatter between the results from individual detectors does not give a consistent result. A combination with ISOPHOT observations taken simultaneously might give a clearer picture (Schulz et al. in preparation).

6.4 Outlook

The ISO project and recent requests from far-IR and submm projects clearly show the need for new calibration sources. The large asteroids can close this gap under the above specified limitations of variability, background influences and the individual TPM uncertainties. There are still many open questions and aspects which have to be addressed in the future. The asteroid shape models together with the absolute size and albedo values are crucial for the TPM and can be improved. First indications of small features in the asteroid spectra have already been seen. But what materials do they correspond to? Would it be possible to include them in a radiative transfer code to be able to explain the empirical emissivity values? Work on identifying the spectral features in the mid-IR is on going and first empirical results have been published (Dotto et al. 2000; Barucci et al. 2001). In the far-IR beyond 30 ${\rm\mu m}$ the extraction of spectral features below a 5% level is more difficult and hardly any laboratory data exist for comparison. For the fainter asteroids in the ISO programme there are more unknowns and larger uncertainties in the TPM input parameters. Instead of connecting them only to stars and planets, one could use in addition the bright asteroids to recalibrate observation and to improve our knowledge about them. Establishing new sources as a reference with better known model parameters would also be an alternative and a possibility to enlarge our small sample of new far-IR standards. The radar techniques are especially promising to improve the quality of existing standards and to establish new ones.

Acknowledgements
We thank Cecile Gry, Martin Burgdorf, Bernhard Schulz (ISO Data Centre in Villafranca), Fred Lahuis (Dutch ISO Data Analysis Centre) and Tanja Lim (Rutherford Appleton Laboratory) for the helpful discussions and support in LWS, SWS and ISOPHOT data reduction. We also would like to express our thanks to the referee, Dr. G. Sandell. His comments were very helpful and improved the paper significantly.

Appendix A: Independent ISOPHOT Observations

A.1 ISOPHOT P3 Observations

A.2 ISOPHOT C100 Observations

A.3 ISOPHOT C200 Observations

Appendix B: Independent LWS Observations

B.1 LWS02 Observations on 1 Ceres

B.2 LWS02 Observations on 2 Pallas

B.3 LWS02 Observations on 4 Vesta

B.4 LWS02 Observations on 10 Hygiea

References

 

Online Material


   
Table A.1: Summary of ISOPHOT P3 observations of asteroids and the corresponding model predictions. The filter photometry was converted to monochromatic fluxes at the specified wavelengths. The estimated errors include observational uncertainties and the estimated method uncertainty (between 5 and 10%), depending on the number of calibrators used and the flux differences between calibrators and asteroids.

Asteroid
Julian Date ${\rm\lambda}$ FD rms TPM
Name (mid time) [ ${\rm\mu m}$] [Jy] [Jy] [Jy]

Pallas
2450289.60547 60 112.72 12.60 100.5
Vesta 2450282.40928 60 327.55 36.62 321.9
Alexandra 2450331.60832 60 26.44 2.96 23.4
Herculina 2450520.05851 60 12.80 1.43 16.0
Pallas 2450289.59131 100 43.47 4.86 48.8
Vesta 2450282.39512 100 127.33 14.24 141.9
Hygiea 2450457.26032 100 15.74 1.76 13.3
Hygiea 2450869.36059 100 17.15 1.92 21.3


 

 
Table A.2: Summary of ISOPHOT C100 observations of asteroids and the corresponding model predictions.

Asteroid
Julian Date ${\rm\lambda}$ FD rms TPM
Name   [ ${\rm\mu m}$] [Jy] [Jy] [Jy]

Ceres
2450575.88955 65 247.84 33.34 206.0
Pallas 2450303.66288 65 86.18 9.64 78.2
Juno 2450275.67887 65 46.20 3.27 39.5
Vesta 2450652.63162 65 187.92 29.35 148.6
Hygiea 2450450.32579 65 37.31 2.64 27.1
Cybele 2450415.29284 65 12.94 1.33 11.2
Pallas 2450303.70392 60 91.88 10.27 89.0
Juno 2450275.69965 60 52.76 3.73 43.4
Vesta 2450652.59831 60 172.89 16.31 170.5
Hygiea 2450450.27524 60 39.62 3.41 33.5
Cybele 2450415.27231 60 15.09 1.30 12.9
Herculina 2450220.66045 60 28.29 2.91 32.7
Ceres 2450575.87253 80 129.70 27.85 151.2
Pallas 2450303.68341 80 47.91 6.45 57.8
Juno 2450275.72060 80 29.07 3.55 25.1
Vesta 2450652.66490 80 100.83 13.57 111.1
Hygiea 2450450.30159 80 25.07 2.80 20.9
Hygiea 2450869.29659 80 25.98 3.03 27.8
Cybele 2450415.25179 80 9.56 1.29 8.6
Herculina 2450255.79757 80 11.54 1.24 11.8
Pallas 2450303.72448 90 48.98 6.37 49.4
Juno 2450275.74111 90 24.22 3.15 19.5
Hygiea 2450444.42684 90 20.53 1.77 22.6
Hygiea 2450869.27307 90 23.30 2.20 23.3
Alexandra 2450331.58479 90 10.26 0.97 11.0
Cybele 2450387.32012 90 10.10 1.04 9.7
Ceres 2450575.85558 100 125.38 21.02 109.2
Pallas 2450303.62181 100 45.54 4.69 43.1
Juno 2450275.76163 100 22.21 2.29 16.6
Vesta 2450652.69807 100 82.70 9.25 77.8
Hygiea 2450444.40270 100 19.32 2.33 19.9
Cybele 2450387.36479 100 9.18 1.03 8.0
Herculina 2450255.75420 100 9.70 0.83 9.9
Herculina 2450527.00830 100 8.25 1.15 8.0
Pallas 2450303.64233 105 38.91 3.67 40.4
Juno 2450275.78214 105 17.47 1.80 16.6
Vesta 2450652.73117 105 63.27 6.51 69.6
Hygiea 2450444.37853 105 19.09 1.80 18.7
Alexandra 2450331.56427 105 7.71 1.00 8.9
Cybele 2450387.34427 105 8.43 1.17 7.6
Herculina 2450527.03248 105 7.01 0.66 7.1



 

 
Table A.3: Summary of ISOPHOT C200 observations of asteroids and the corresponding model predictions.

Asteroid
Julian Date ${\rm\lambda}$ FD rms TPM
Name   [ ${\rm\mu m}$] [Jy] [Jy] [Jy]

Ceres
2450313.65701 120 128.21 12.82 126.3
Ceres 2450616.79198 120 125.06 10.69 118.5
Pallas 2450289.52763 120 34.95 3.50 37.5
Juno 2450457.32374 120 18.85 2.27 17.9
Vesta 2450282.33141 120 101.17 9.05 104.8
Hygiea 2450444.32073 120 13.35 1.34 14.9
Cybele 2450244.99559 120 7.19 0.61 7.4
Herculina 2450220.71909 120 8.45 0.76 9.6
Ceres 2450149.86646 150 109.58 7.02 94.4
Ceres 2450313.67502 150 86.39 6.11 84.7
Ceres 2450631.77793 150 94.93 7.41 92.8
Pallas 2450289.54560 150 23.12 1.35 25.3
Juno 2450457.34175 150 11.93 0.76 10.9
Vesta 2450282.34939 150 59.85 3.49 67.7
Hygiea 2450444.33875 150 9.58 0.49 10.2
Cybele 2450245.01358 150 5.48 0.35 5.1
Herculina 2450220.73706 150 5.41 0.32 6.4
Ceres 2450149.89059 170 77.45 8.66 73.7
Ceres 2450313.69306 170 65.94 3.84 66.2
Ceres 2450616.81978 170 62.05 3.97 62.2
Pallas 2450289.56360 170 18.15 2.03 19.6
Juno 2450457.35978 170 8.58 0.81 8.1
Vesta 2450282.36740 170 45.02 4.64 52.3
Hygiea 2450444.35679 170 7.36 0.63 8.1
Hygiea 2450869.34271 170 5.85 0.60 8.4
Cybele 2450245.03159 170 3.37 0.24 4.0
Ceres 2450313.63899 180 57.23 3.34 58.9
Ceres 2450631.74854 180 60.17 5.68 64.5
Pallas 2450289.50964 180 16.71 1.18 17.6
Juno 2450457.30569 180 8.69 0.56 9.0
Vesta 2450282.31343 180 43.85 4.14 48.7
Hygiea 2450444.30271 180 6.82 0.40 7.0
Cybele 2450244.97759 180 3.33 0.21 3.4
Cybele 2450415.23097 180 2.18 0.15 2.3
Herculina 2450220.70110 180 4.45 0.42 4.7
Ceres 2450149.91472 200 56.09 5.29 52.4
Ceres 2450313.62098 200 50.23 4.74 47.1
Ceres 2450616.76421 200 42.45 3.32 44.2
Pallas 2450289.49167 200 12.97 0.92 13.9
Juno 2450457.28769 200 7.61 0.59 7.4
Vesta 2450282.29545 200 35.64 2.28 39.3
Hygiea 2450444.28470 200 5.47 0.32 5.5
Cybele 2450244.95961 200 2.69 0.23 2.7
Cybele 2450415.21300 200 1.72 0.13 1.8
Herculina 2450220.68312 200 3.76 0.32 4.0



 

 
Table B.1: Summary of LWS02 observations on Ceres. The calibrated monochromatic flux densities, the errors and the model predictions are given in [Jy]. Together with the object name and the ISO specific number (TDT), the Julian Dates of observation mid time are tabulated.

Wavelengths [ ${\rm\mu m}$] / Flux [Jy]
  46.2 56.2 66.1 75.7 84.8 102.4 122.2 141.8 160.6 178.0

1 Ceres

TDT09300401 JD2450131.68839          
Obs. 395.0 399.2 329.0 260.3 227.5 155.4 123.3 85.7 66.0 64.6
rms. 28.8 25.6 20.3 16.7 14.6 9.6 8.2 5.8 4.6 6.4
TPM 512.0 401.6 320.6 261.9 219.1 160.6 117.8 89.7 70.8 58.0

1 Ceres

TDT10500402 JD2450143.72387          
Obs. 447.1 427.9 377.5 297.5 235.7 177.8 134.3 104.9 67.8 64.5
rms. 32.5 26.4 23.0 18.3 15.1 11.0 8.6 7.1 4.6 6.0
TPM 577.2 452.8 361.6 295.3 247.1 181.1 132.9 101.2 79.9 65.4

1 Ceres

TDT11900214 JD2450157.61316          
Obs. 498.1 473.9 417.0 333.9 273.8 201.0 151.7 111.8 87.4 81.0
rms. 36.3 29.2 25.7 20.6 17.5 12.4 10.1 7.8 6.1 7.5
TPM 657.0 515.7 411.9 336.5 281.6 206.5 151.5 115.3 91.1 74.6

1 Ceres

TDT12600114 JD2450164.58605          
Obs. 551.5 509.1 437.5 349.8 315.2 203.3 156.0 116.1 92.1 74.2
rms. 39.8 31.4 27.0 21.6 20.2 12.5 10.0 7.9 6.2 8.8
TPM 713.6 560.1 447.3 365.4 305.8 224.2 164.5 125.2 98.9 81.0

1 Ceres

TDT25800302 JD2450296.27109          
Obs. 676.9 511.5 425.9 336.7 298.9 207.3 172.2 121.5 82.7 61.4
rms. 48.8 31.5 26.3 20.8 19.1 12.6 11.0 8.2 5.8 5.7
TPM 665.9 525.7 421.5 345.3 289.6 213.0 156.6 119.5 94.5 77.4

1 Ceres

TDT26500301 JD2450303.25720          
Obs. 596.4 467.7 386.1 309.1 265.1 187.5 156.5 112.4 76.9 58.9
rms. 43.0 28.8 23.8 19.8 18.0 11.6 10.0 7.9 5.4 5.2
TPM 607.5 479.8 384.9 315.5 264.6 194.7 143.2 109.3 86.4 70.8

1 Ceres

TDT32100204 JD2450359.10520          
Obs. 331.3 254.0 220.8 174.4 149.4 109.5 79.8 56.7 44.2 32.7
rms. 24.1 17.2 13.6 11.8 9.6 6.7 5.1 4.0 3.5 4.1
TPM 346.5 273.7 219.6 179.9 150.9 111.0 81.6 62.3 49.2 40.4

1 Ceres

TDT53802209 JD2450575.96948          
Obs. 398.3 278.7 230.5 174.6 146.4 120.7 81.8 59.6 47.3 31.1
rms. 29.8 17.8 14.8 11.2 10.7 7.7 5.2 4.2 3.7 5.0
TPM 354.2 279.9 224.5 183.9 154.2 113.4 83.4 63.6 50.3 41.2

1 Ceres

TDT57902409 JD2450616.87244          
Obs. 502.7 407.8 324.9 269.7 225.1 167.5 124.3 93.5 69.5 48.8
rms. 37.6 25.1 20.0 16.6 15.3 10.3 8.2 6.5 4.9 4.3
TPM 512.8 405.4 325.3 266.6 223.6 164.5 121.0 92.3 72.9 59.8

1 Ceres

TDT59401908 JD2450631.73749          
Obs. 595.5 475.4 389.5 319.0 270.4 199.4 140.7 117.0 85.1 78.7
rms. 43.4 29.3 24.0 20.4 17.3 12.8 9.0 8.2 6.0 7.3
TPM 601.2 474.9 380.9 312.0 261.7 192.4 141.4 107.8 85.3 69.8

1 Ceres

TDT72001901 JD2450757.54860          
Obs. 521.3 398.4 319.2 263.3 216.3 164.9 117.9 89.8 70.8 59.2
rms. 38.0 24.6 19.7 16.9 13.9 10.0 7.6 6.3 5.2 6.2
TPM 511.5 405.3 325.8 267.4 224.5 165.4 121.8 93.0 73.6 60.3

1 Ceres

TDT74803304 JD2450785.62839          
Obs. 389.7 301.6 243.0 209.3 166.7 126.6 99.0 75.7 -- --
rms. 28.4 20.5 15.0 13.4 11.3 7.7 6.3 5.3 -- --
TPM 389.5 308.7 248.2 203.7 171.1 126.0 92.8 70.9 56.1 46.0

1 Ceres

TDT75503003 JD2450792.56800          
Obs. 357.5 286.2 237.7 192.6 162.9 120.2 83.3 62.7 -- --
rms. 26.0 18.3 14.7 13.1 13.9 7.4 6.2 4.4 -- --
TPM 361.2 286.2 230.1 188.9 158.6 116.9 86.1 65.7 52.0 42.6

1 Ceres

TDT76200502 JD2450799.06191          
Obs. 336.5 277.2 225.5 181.0 150.3 111.5 80.5 56.9 -- --
rms. 24.5 17.1 13.9 12.3 9.6 7.1 5.3 4.0 -- --
TPM 348.2 275.9 221.7 181.9 152.8 112.5 82.8 63.2 50.0 41.0

1 Ceres

TDT76903102 JD2450806.54772          
Obs. 318.8 257.8 207.2 171.9 151.8 109.5 80.1 61.4 40.3 35.5
rms. 23.2 15.9 12.8 12.5 11.1 7.0 5.3 4.8 3.1 4.5
TPM 323.8 256.4 206.1 169.1 141.9 104.5 77.0 58.7 46.5 38.1



 

 
Table B.2: Summary of LWS02 observations on Pallas.

Wavelengths [ ${\rm\mu m}$] / Flux [Jy]
  46.2 56.2 66.1 75.7 84.8 102.4 122.2 141.8 160.6 178.0

2 Pallas

TDT23000306 JD2450268.35000          
Obs. 141.1 156.3 131.4 95.9 91.6 57.0 47.4 37.8 25.1 23.3
rms. 11.0 10.6 9.6 8.2 5.9 3.5 3.3 2.8 2.3 4.1
TPM 193.2 151.4 120.7 98.5 82.3 60.3 44.2 33.5 26.4 21.6

2 Pallas

TDT25100202 JD2450289.28410          
Obs. 149.6 118.4 99.4 74.4 68.2 48.7 35.6 28.9 19.4 --
rms. 11.7 8.6 6.4 5.9 6.3 3.1 2.9 2.3 2.0 --
TPM 153.8 120.8 96.5 78.8 65.9 48.3 35.4 26.9 21.2 17.4

2 Pallas

TDT26500503 JD2450303.28000          
Obs. 132.9 107.2 83.6 62.3 50.1 42.5 30.4 21.7 17.9 --
rms. 11.0 7.8 5.7 4.5 6.3 3.1 2.4 1.7 2.2 --
TPM 135.1 106.2 84.9 69.3 58.0 42.6 31.2 23.7 18.7 15.3

2 Pallas

TDT27200203 JD2450310.22194          
Obs. 123.3 92.2 85.0 63.3 52.8 38.8 28.8 23.2 14.8 --
rms. 10.2 9.3 6.7 5.4 4.5 2.6 2.2 1.9 1.4 --
TPM 128.2 100.8 80.5 65.8 55.0 40.4 29.6 22.5 17.7 14.5



 

 
Table B.3: Summary of LWS02 observations on Vesta.

Wavelengths [ ${\rm\mu m}$] / Flux [Jy]
  46.2 56.2 66.1 75.7 84.8 102.4 122.2 141.8 160.6 178.0

4 Vesta

TDT24402202 JD2450282.43072          
Obs. 424.2 339.4 280.4 212.5 179.2 124.7 105.0 75.5 46.8 33.5
rms. 30.6 21.7 18.0 13.1 12.2 7.7 7.0 5.3 3.7 3.3
TPM 457.6 355.8 282.2 229.4 191.2 139.6 102.1 77.3 60.8 49.7

4 Vesta

TDT80500101 JD2450841.94186          
Obs. 149.7 116.0 94.8 78.2 54.8 47.2 36.7 27.5 18.5 15.6
rms. 11.7 7.9 6.1 5.7 5.5 3.2 3.0 2.0 2.2 3.3
TPM 147.0 116.0 92.9 76.0 63.7 46.8 34.4 26.2 20.6 16.9

4 Vesta

TDT80500104 JD2450841.96022          
Obs. 156.6 112.5 93.2 77.5 64.9 49.5 38.2 31.7 -- --
rms. 11.7 8.2 6.3 6.1 5.5 3.1 3.1 2.3 -- --
TPM 146.3 115.5 92.6 75.8 63.5 46.7 34.3 26.1 20.6 16.8

4 Vesta

TDT80500107 JD2450841.97855          
Obs. 142.4 111.3 99.7 76.6 72.0 49.1 39.2 32.5 -- --
rms. 10.7 7.1 6.4 6.5 9.1 3.3 2.6 2.7 -- --
TPM 148.6 117.2 94.0 76.9 64.4 47.4 34.8 26.5 20.9 17.1

4 Vesta

TDT80500110 JD2450841.99691          
Obs. 142.2 114.9 93.9 76.0 70.7 48.8 33.8 33.5 -- --
rms. 10.6 7.8 5.8 6.0 10.2 3.3 2.9 2.5 -- --
TPM 151.4 119.4 95.7 78.3 65.6 48.2 35.4 26.9 21.2 17.4

4 Vesta

TDT80500113 JD2450842.01524          
Obs. 151.5 108.6 98.2 68.0 61.2 48.5 34.8 29.5 -- --
rms. 11.8 8.6 6.3 5.4 6.2 3.1 2.8 2.7 -- --
TPM 152.2 120.0 96.1 78.7 65.9 48.4 35.6 27.0 21.3 17.4

4 Vesta

TDT80500116 JD2450842.03362          
Obs. 148.3 120.7 92.6 78.6 64.9 48.9 40.0 31.0 -- --
rms. 12.2 10.3 5.9 6.2 7.0 3.6 2.8 2.4 -- --
TPM 150.1 118.3 94.8 77.5 64.9 47.7 35.1 26.7 21.0 17.2

4 Vesta

TDT80500119 JD2450842.05196          
Obs. 158.0 128.0 95.7 75.9 64.6 52.1 39.2 28.4 -- --
rms. 13.0 8.2 6.1 4.9 5.1 3.3 2.9 2.4 -- --
TPM 147.1 116.0 93.0 76.1 63.7 46.8 34.4 26.2 20.6 16.9

4 Vesta

TDT80500122 JD2450842.07031          
Obs. 148.4 130.2 96.4 80.0 59.9 47.5 40.5 29.8 -- --
rms. 11.1 8.8 6.2 5.8 6.0 3.0 2.8 2.3 -- --
TPM 146.1 115.3 92.4 75.7 63.4 46.6 34.3 26.1 20.5 16.8

4 Vesta

TDT80500125 JD2450842.08865          
Obs. 157.2 127.8 97.2 81.9 48.8 48.5 39.5 24.7 -- --
rms. 11.8 10.1 7.7 6.0 5.3 3.0 3.0 2.1 -- --
TPM 148.1 116.9 93.7 76.7 64.2 47.2 34.7 26.4 20.8 17.0

4 Vesta

TDT80500128 JD2450842.10700          
Obs. 149.9 117.3 96.8 85.5 48.2 46.0 26.2 -- -- --
rms. 12.4 9.2 6.2 7.3 8.2 2.9 2.0 -- -- --
TPM 151.1 119.2 95.5 78.2 65.5 48.1 35.4 26.9 21.2 17.4

4 Vesta

TDT80500131 JD2450842.12534          
Obs. 140.2 116.9 89.3 82.6 57.1 49.1 36.3 27.8 -- --
rms. 11.0 7.9 5.7 7.1 5.3 3.3 2.9 2.2 -- --
TPM 152.1 119.9 96.1 78.6 65.8 48.4 35.5 27.0 21.3 17.4

4 Vesta

TDT80500134 JD2450842.14369          
Obs. 147.6 119.8 100.2 82.2 70.7 48.9 35.7 26.9 -- --
rms. 11.5 8.1 6.8 6.0 8.3 3.0 2.7 2.8 -- --
TPM 150.0 118.3 94.7 77.5 64.9 47.7 35.1 26.6 21.0 17.2

4 Vesta

TDT80500137 JD2450842.16203          
Obs. 142.7 123.5 99.7 74.0 50.1 52.9 38.3 26.6 -- --
rms. 10.7 8.4 7.3 5.0 4.6 4.2 3.1 2.0 -- --
TPM 147.0 116.0 92.9 76.0 63.7 46.8 34.4 26.2 20.6 16.9



 

 
Table B.4: Summary of LWS02 observations on Hygiea.

Wavelengths [ ${\rm\mu m}$] / Flux [Jy]
  46.2 56.2 66.1 75.7 84.8 102.4 122.2 141.8 160.6 178.0

10 Hygiea

83201803 2450869.44722          
Obs. 57.1 52.9 50.7 37.4 35.4 23.3 17.2 11.1 -- --
rms. 5.7 5.3 5.9 3.2 4.1 2.0 1.3 1.3 -- --
TPM 68.4 54.5 43.9 36.1 30.3 22.4 16.5 12.6 9.9 8.1

10 Hygiea

85303402 2450890.42337          
Obs. 46.7 50.8 44.4 25.9 29.7 19.4 14.7 11.8 -- --
rms. 3.9 4.3 4.5 3.0 7.3 1.4 1.2 1.3 -- --
TPM 53.0 42.3 34.1 28.0 23.5 17.4 12.8 9.8 7.7 6.3



Copyright ESO 2001