next previous
Up: A catalog of observed comets


Subsections

5 Sources of errors and uncertainties

Several sources of uncertainties can be identified in the determination of the adopted nuclear magnitudes and in the computation of the effective nuclear radius. Let us now analyze these in some detail.

We may broadly distinguish between methodological errors due to shortcomings of the observing techniques and intrinsic errors due to unknown properties of the nuclei themselves. The former include all random and systematic errors of photometry, which are often particularly bad for attempted measurements on cometary nuclei. Some of the reasons are that many of the data were taken for non-photometric purposes where magnitudes resulted as a spinoff without much attention to the photometric conditions of the sky, and that the data often refer to objects near the plate limit without appropriate standard sequences, or uncertainty of the CCD extrapolation.

It must be kept in mind that most of the magnitudes in the literature, even if referred to as nuclear, are not truly nuclear magnitudes according to our physical definition. This is the main concern of our investigation, and our most important aim is to minimize these errors, which are due to coma contamination. Nonetheless, the discrepancies remaining between our adopted nuclear magnitudes and the true values have an important component that can be referred to as residual nuclear activity. This component borders between the methodological and intrinsic errors, because:

(1)
the angular resolution of the observation has a major influence on the amount of coma contamination for an active comet (methodological error);

(2)
even if the best possible angular resolution is achieved, or if the comet appears perfectly inactive, the nucleus is never resolved and the apparent nucleus may always be contaminated by circumnuclear dust (intrinsic error).

To the above errors, we should finally add some further intrinsic errors due to the nucleus itself: variations of the photometric nucleus cross-section S due to rotation; and a phase angle effect. Furthermore, the computation of the nuclear effective radius requires precise knowledge of the geometric albedo, which is generally not available.

5.1 The adopted nuclear magnitudes

5.1.1 The data values

The error of each individual observed magnitude may be of the order of 0.1 magnitudes or more, except for the case of CCD observations obtained from a few papers where photometry was done carefully; errors as low as a few hundredths of a magnitude can be attained in this case. A few tenths of a magnitude is the usual error of photographic observations. While most CCD magnitudes reported in the MPCs are more indicative values than precise determinations, the great improvement of the recent CCD data with respect to the old photographic data is the huge number of magnitude determinations at larger heliocentric distances. Unfortunately, in many cases the color of the magnitude is not clearly established, or it can only be crudely approximated by a standard color, as occurs, for instance, with Scotti's unfiltered magnitudes.

5.1.2 Nuclear rotation

Cometary nuclei seem to be very irregular and elongated objects. Therefore, the photometric nucleus cross-section may vary with the position of the nucleus with respect to the line of sight. If the spin axis of the nucleus is not aligned with the line of sight, rotation produces fast variations of the photometric cross-section S. Rotational lightcurves have been used to determine the rotation period of cometary nuclei (see e.g. Jewitt & Meech 1987). Assuming that the nucleus is a triaxial ellipsoid with axes a > b > c, that it has a homogeneous surface, and that its spin is relaxed to pure rotation around the fixed axis c, it is easy to see that the maximum amplitude of its rotational lightcurve is


\begin{displaymath}B_{\rm N2}/B_{\rm N1}=b/a=10^{-0.4 \Delta m}
\end{displaymath} (3)

where $B_{\rm N1}$ and $B_{\rm N2}$ are the maximum and minimum brightnesses. For instance, the lightcurve of a nucleus with a/b=2 fulfilling the previous conditions will show a maximum amplitude of 0.75 mag. So, for typical complex shapes, variations up to the order of one magnitude may be observed due only to the rotation of the comet (see e.g. Lamy et al. 1998b).

5.1.3 Phase coefficient $\beta $

As regards the phase coefficient $\beta $, Scotti (private communication) has pointed out that it might depart significantly from our adopted value. Since we have taken a limit for the phase angle of $30^\circ$, an uncertainty on $\beta $ of $\pm\ 0.02$ would introduce a maximum uncertainty of $\pm\ 0.6$ in the adopted nuclear magnitude.

The three effects discussed so far imply that, even for a non-active nucleus, a dispersion of the order of at least one magnitude in the plots should be considered as normal.

5.1.4 Non-detected activity

Some comets present much larger dispersions or a trend in the plots (see Sect. 8.3) that suggests some non-detected activity. The detection of a small coma, especially when the comet is far from the Earth, is not a simple task. Licandro et al. (1999a) compare the comet brightness profile to the profile of field stars, which proves to be very useful, but even with such methods some faint coma may remain undetected. But for most of the data included in the catalog the presence of a coma is reported only when it is easily detected by eye or, in the best case, by comparing the FWHM of the comet with respect to the one of the stars in CCD images.

The only way to determine if the comet is active with the data of the catalog is to compare the nuclear magnitudes obtained at a wide range of heliocentric distances. If the trend of the data becomes horizontal beyond a certain heliocentric distance, and the dispersion of the data is not large (say, less than one magnitude), we should expect that the corresponding average magnitude will correspond to the one of the bare nucleus. Yet, some residual dust coma might remain at some constant level over the whole observed part of the cometary orbit, even if the comet does not present any signs of activity. The main problem with this uncertainty, with respect to the other ones, is that it introduces a systematic effect that produces brighter "nuclear'' magnitudes.

5.1.5 Coma subtraction

A method was already explained in Sect. 3.2 as developed and applied by James Scotti. We showed that it involves in itself several sources of errors and uncertainties.

The HST team also applied a method of coma subtraction to their observations. A narrow PSF and a very good sampling enables them to obtain a more reliable estimate of the nucleus contribution, even in cases similar to the one sketched in Fig. 1a.

5.2 The computation of the nuclear radius

Let us also consider the uncertainties in the determination of the effective nuclear radius. Let us consider an uncertainty in the geometric albedo of $\Delta p_{\rm v}$ and in the absolute magnitude of $\Delta H_{\rm N}$, leading to a combined uncertainty in the radius of $\Delta R_{\rm N}$. Applying Eq. (1) to $(R_{\rm N}, H_{\rm N}, p_{\rm v})$ and $(R_{\rm N}+\Delta R_{\rm N}, H_{\rm N}+\Delta H_{\rm N}, p_{\rm v}+\Delta p_{\rm v})$ and taking the difference, we obtain


\begin{displaymath}\log \left( 1\! +\! \frac{\Delta R_{\rm N}}{R_{\rm N}} \right...
...m v}}{p_{\rm v}} \right)\! +\! 0.4 \ \Delta H_{\rm N} \right].
\end{displaymath} (4)

For example, an uncertainty of $\Delta H_{\rm N} = \pm 1$ mag would lead to a -37%/+58% relative error in the radius estimate. Adopting a value $p_{\rm v} = 0.04$ and considering an uncertainty of $\pm\ 0.02$, the uncertainty in the radius estimate would be -44%/+78%. Note that, for radius determination, a one magnitude error in the nuclear magnitude is even less important than the uncertainty due to the bad knowledge of cometary albedo. Individual nuclear magnitudes to a precision much better than one magnitude are not needed for a first characterization of the size or mass distribution either. We are thus justified in considering all $H_{\rm N}$ values with errors within roughly one magnitude as meaningful data.


next previous
Up: A catalog of observed comets

Copyright The European Southern Observatory (ESO)