A&A 380, 727-733 (2001)
DOI: 10.1051/0004-6361:20011412
A. Vienne1,2 - W. Thuillot1 - C. H. Veiga3 - J.-E. Arlot1 - R. Vieira Martins3
1 - Institut de mécanique céleste et de calcul des éphémérides - Observatoire
de Paris, UMR 8028 du CNRS, 77 avenue Denfert-Rochereau, 75014 Paris, France
2 -
Université de Lille, 9000 Lille, France
3 -
Observatório Nacional, Rua Gal. José Cristino 77, 20921-400 Rio de Janeiro,
Brazil
Received 7 May 2001 / Accepted 4 October 2001
Abstract
This paper provides an analysis of astrometric measurements
of the main Saturnian satellites made thanks to CCD observations performed in
1995 at the Laboratório Nacional de Astrofísica at Itajubá in Brazil.
The astrometric reduction is discussed especially the small corrections
done here, but most of time neglected elsewhere. A catalog of 6006 differential
positions has been obtained. They have been compared to different ephemerides,
the Vienne & Duriez ephemerides (TASS 1.7), the Harper & Taylor ephemerides
and the Dourneau ephemerides. These observations provide a large set of modern
observations, and appear to be of good precision. This accuracy is needed for
future use of these data to improve the dynamical models. These positions are
included in the data base NSDC dedicated to the natural satellites (ftp://ftp.bdl.fr/pub/NSDC/saturn/raw_data/position/).
They are also available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr
(130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/380/727
Key words: Saturn - satellites of Saturn - astrometry
In 1995, an international campaign was organized in order to coordinate the observations of phenomena among the Saturnian satellite system (Arlot & Thuillot 1993). During this period the Earth and Sun crossed the orbital planes of these satellites and mutual occultations and eclipses occurred. Many observers then planned to observe these events in order to catch the short light-flux drops most of times thanks to CCD cameras. Sixty-six such observations of events were successfully performed in several countries (Thuillot et al. 2001). But actually many CCD frames were then obtained without any detection of events because of the difficulty to measure a faint magnitude drop closely to the planet. Here we intend to extract the astrometric data from these observations in order to get a catalog of positions jointly with the observations of the events themselves.
These positions could subsequently be applied to the improvement of the dynamical models and to explore some specific problems. For example accurate positions of both Mimas and Tethys are necessary to well determine the eccentricity of Tethys which has a major influence on the evolution of the Mimas-Tethys resonance and lead to a better knowledge of its dynamical history (Champenois & Vienne 1999a, 1999b).
First we present a description of these measurements issued from the observations made at the Laboratório Nacional de Astrofísica at Itajubá in Brazil.We present our reduction, and after we show the comparison of the obtained positions with those deduced from the ephemerides. And at last, we give the catalog of the corresponding astrometric data. These data could also be accessed under a standard form on the data base NSDC dedicated to the natural satellites (ftp://ftp.bdl.fr/pub/NSDC/saturn/raw_data/position/).
In Brazil, campaigns of astrometric observations of the natural satellites are performed (see for example Veiga & Vieira Martins 1999). During the 1995 international campaign of observation of the Saturnian mutual events, several nights were scheduled at the Laboratório Nacional de Astrofísica at Itajubá to observe the Saturnian satellites. The number of observations and its distribution over the period 1995-1996 are given in Fig. 1.
![]() |
Figure 1: Histogram of the observations with respect to times. Each bar corresponds to one night. |
Open with DEXTER |
The two telescopes used at Itajubá (32
19
East in longitude,
South in latitude, altitude 1872 m) were the 1.60 m and the 0.60 m reflectors.
They were equipped with CCD cameras the characteristics of which are given in
Table 1.
Camera | Description | series |
CCD0048 | EEV P88231 (770X1152) | 1 |
CCD009 | EEV-05.20.0.202 (770X1200) | 2-3-4 |
CAM2 | SITe SI003AB (1024X1024) | 6 |
All the frames were analyzed thanks to the ASTROL software (Colas 1996) which allowed us to use a centering algorithm based on the adjustment of a Point Spread Function (PSF).
No stars could generally be recorded, but in several frames some stars were although visible.
Several frames where stars are visible close to the satellites have been reduced thanks to the ESO software MIDAS (Warmels 1991) and a new algorithm for astrometric reduction named PRIAM recently developed by Fienga et al. (in preparation). This algorithm of reduction leads in these case to topocentric equatorial J2000 coordinates of the satellites and to the determination of the scale and orientation of the frames. Unfortunately, these frames are too few, so we have used these frames and the corresponding reductions only to validate our inter-satellite reduction.
We give in this section the main points of the reduction, we can find in Vienne et al. (2001) more details and a discussion of the method for the case of inter-satellite measurements.
The frames registered and measured in pixels are not directly comparable with
ephemeris, because they have been affected by some local effects. Before
any comparison, it is necessary to take into account these effects due to refraction,
stellar aberration, central projection (projection of the celestial sphere on
the tangential plane of the focal point), light-travel time between satellites
and topocentric parallax. These differential astrometric corrections are classical
and we have made them in the present work. In the past, the precision level
was about
for inter-satellite astrometry. It is surely the reason
that until now, concerning the inter-satellite observations of satellites of
Saturn, all these astrometric corrections were not automatically taken into
account by previous authors:
effect | correction |
refraction | 1'' (
![]() |
stellar aberration |
![]() |
central projection (typical value) |
![]() |
light-travel time between satellites |
![]() ![]() |
topocentric parallax |
![]() |
relativistic deflexion from the Sun |
![]() ![]() ![]() |
relativistic deflexion from the planet |
![]() ![]() |
The list of the differential corrections, to be applied to the five outer planetary
sytems, is given in Table 2. The values given in the right
column are indicative; they correspond to a maximum value, for a differential
effect applied to a standard frame 400'' large. In the present set of
observations, the geocentric elongation of Saturn from the Sun is greater than
85,
so the differential correction due to the relativistic deflexion
of the Sun is less than 0.009 mas. Due to Saturn, this effect is about 0.02 mas in the position of Iapetus. So, in the present reduction, we have neglected
these two effects.
The positions of Saturn are given by the ephemerides SLP96 from the "Institut
de mécanique céleste (IMCCE)'' (available at ftp://ftp.bdl.fr/pub/ephem/sun/slp96/)
found on the VSOP87 planetary theory (Bretagnon & Francou 1988). The
precision on the positions of Saturn is about
which is enough here.
We use three different ephemerides for the saturnicentric positions of the satellites:
TASS1.7 (Vienne & Duriez 1995; Duriez & Vienne 1997), H-T (Harper
& Taylor 1993) and Dourneau (1993). The time argument is the
Terrestrial Time (TT). At a given date corrected for the light-travel
time between the satellites, we compute the apparent coordinates, and
,
of the objects. These coordinates take into account all the
corrections listed above. We use the observed frame as a table of positions
,
given in pixels (see Sect. 2). So, with
the relation
, we compute an estimation
of a and
b by a least square procedure. We deduce the scale factor
and the orientation of the receptor (the angle from the
-axis of
the camera to the true equator)
.
Note that, because of the linearity of the relation between
and
,
no iterative procedure is necessary. It is not true
for the direct determination of
and
as it is done
by Shen et al. (2001). Using
and
,
and
applying once more time, but inversely, the local corrections, we get the coordinates
of each measured object. The coordinates
are astrometric, but there is an absolute part
and
which does not come from the observations. C is defined such as the frame
is the tangential plane of the celestial sphere at C. We have to know
the position of this point in order to make the corrections of the central projection.
As C does not correspond to a physical object, we do not know
.
We estimate
by an iterative procedure beginning
by applying the reduction without the present correction. Two iterations are
enough to obtain a good estimation. In order to avoid to interpret
the observations as absolute ones, we publish them in inter-satellite form
,
or more precisely,
and
.
The index
is for the satellite object, and the index
for the reference satellite. The content of the catalog is described in Sect. 5.
The procedure is applied once only for all frames of a given series. Generally, a series of frames corresponds to one campaign of observations and covers several nights. We then suppose that the receptor has been mounted in the same way for all the frames of the series.
This procedure is efficient for the eight major satellites of Saturn, or more generally for all objects for which we have ephemerides. But, in the least square procedure, only the positions of Tethys, Dione, Rhea and Titan are used to calibrate because these satellites have the best ephemerides, probably affected by the smallest systematic effects.
The positions
we have obtained
are really astrometric ones because they are given in the J2000 system and all
significant astrometric corrections have been done. But, for a given series,
they are given apart from a scale factor and from a rotation. As explained in
Vienne et al. (2001), if we want to compute the astrometric coordinates
in any other way, for example with other ephemerides, we have only to touch
up
and
.
![]() |
Figure 2:
Histogram for
![]() ![]() |
Open with DEXTER |
We use three different ephemerides for the saturnicentric positions of the satellites: TASS1.7, H-T and Dourneau. The position of Saturn itself is given by the ephemerides SLP96 (see above). The time argument is the Terrestrial Time (TT).
The observations rejected are the same for the three theories: an observation
is rejected when the residual is greater than
for each of the three
theories. In this way, we have rejected 6 observations of Mimas, 2 of Tethys,
1 of Titan and 12 of Hyperion.
![]() |
Figure 3:
Histogram for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Histogram for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 5:
Histogram for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 6:
Histogram for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Histogram for
![]() ![]() |
Open with DEXTER |
![]() |
Figure 8:
Histogram for
![]() ![]() |
Open with DEXTER |
TASS |
![]() |
![]() |
|||
satellites | N | means | rms | means | rms |
S1-Mimas | 216 | -23 | 85 | -1 | 78 |
S2-Enceladus | 865 | +11 | 92 | -6 | 68 |
S3-Tethys | 2151 | +3 | 79 | +3 | 65 |
S4-Dione | 1466 | -8 | 64 | -1 | 57 |
S6-Titan | 460 | +19 | 81 | +11 | 67 |
S7-Hyperion | 324 | -93 | 152 | -39 | 128 |
S8-Iapetus | 524 | -108 | 141 | +10 | 69 |
S3 S4 S5 S6 | 4077 | +1 | 74 | +2 | 62 |
H-T |
![]() |
![]() |
|||
satellites | N | means | rms | means | rms |
S1-Mimas | 216 | -54 | 108 | +12 | 80 |
S2-Enceladus | 865 | -3 | 90 | -10 | 69 |
S3-Tethys | 2151 | 0 | 79 | 0 | 65 |
S4-Dione | 1466 | -8 | 62 | +12 | 58 |
S6-Titan | 460 | 14 | 84 | +9 | 71 |
S7-Hyperion | 324 | -295 | 356 | +115 | 170 |
S8-Iapetus | 524 | +422 | 432 | +149 | 163 |
S3 S4 S5 S6 | 4077 | -2 | 74 | +5 | 63 |
Dourneau |
![]() |
![]() |
|||
satellites | N | means | rms | means | rms |
S1-Mimas | 216 | -63 | 106 | +6 | 78 |
S2-Enceladus | 865 | -2 | 90 | -6 | 71 |
S3-Tethys | 2151 | +8 | 81 | +2 | 67 |
S4-Dione | 1466 | -1 | 64 | +8 | 59 |
S6-Titan | 460 | +3 | 89 | +8 | 74 |
(S7-Hyperion) | 324 | +1892 | 1924 | +98 | 167 |
S8-Iapetus | 524 | +206 | 225 | -34 | 78 |
S3 S4 S5 S6 | 4077 | +4 | 76 | +5 | 65 |
The statistics of the comparison is shown in Table 3.
Figures 2 to 8 show the corresponding histograms
(for TASS only). From Table 3 and whatever the theory used,
we can deduce that the dispersion of the observations is about
,
and the bias of some milli-arcseconds. As Tethys, Dione, Rhea and Titan have
been used in the fit, this fact is true for these satellites taken separately
or globally. For Enceladus, the r.m.s of the residuals reaches
,
also whatever the used theory. It is not the case for Mimas, Hyperion and Iapetus.
For these satellites, TASS gives the lowest residuals:
(Mimas),
(Hyperion),
(Iapetus) for
,
and
(Mimas),
(Hyperion),
(Iapetus)
for
.
Each one corresponds to several hundred positions:
216 for Mimas, 324 for Hyperion and 524 for Iapetus.
For comparison, Shen et al. (2001) have analyzed CCD observations over
the 1990-1997 period from Qiao et al. (1999), Harper et al. (1997,
1999). From their analysis, they conclude that TASS is significantly
more accurate (
)
than other theories (
). They give
the following residuals:
(Mimas),
(Hyperion),
(Iapetus), corresponding respectively to 57, 218 and 230 positions. So, the
observations presented in our paper cover a large set of modern observations,
and are of good precision.
Many observations of Hyperion have bad residuals when the theory of Dourneau has been used. The reason is that Dourneau did not include any short period term in the expressions of Hyperion's theory. Note also, that Shen et al. (2001) use for the theory of Harper & Taylor (H-T) a slight different version not yet published.
opp | year | m | day (utc) | tt-utc | obs | ref | t | obj | fl | obs1 | obs2 | s | f | O-C1 | O-C2 | r | - | series | xpix | ypix |
s | arcsec (degree for C*) | arcsec | pixel | |||||||||||||||||
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 1 | 65 | 11 | -30.2330 | 0.0060 | 2 | 1 | 0.051 | 0.078 | 5 | 0 | 1 | -102.657 | -0.464 |
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 1 | 25 | 11 | 43.7666 | -4.3149 | 2 | 1 | -0.048 | 0.001 | 5 | 0 | 1 | 148.679 | -13.953 |
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 1 | 45 | 11 | 94.1876 | -8.5286 | 2 | 1 | -0.107 | 0.021 | 5 | 0 | 1 | 319.951 | -27.459 |
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 1 | 35 | 11 | 110.1468 | -8.6285 | 2 | 1 | -0.123 | -0.040 | 5 | 0 | 1 | 374.142 | -27.545 |
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 1 | C5 | 11 | 40.4228 | -4.6831 | 2 | 1 | 999.999 | 999.999 | 5 | 0 | 1 | 137.386 | -15.259 |
118 | 1995 | 6 | 10.3531230 | 61.184 | 874 | 301 | 0 | C* | 11 | 355.5182 | -4.1346 | 2 | 1 | 999.999 | 999.999 | 5 | 0 | 1 | 99999.999 | 99999.999 |
118 | 1995 | 6 | 10.3540727 | 61.184 | 874 | 301 | 1 | 65 | 11 | -30.1660 | -0.0413 | 2 | 1 | 0.026 | 0.038 | 5 | 0 | 1 | -102.429 | -0.624 |
... | ||||||||||||||||||||
118 | 1995 | 8 | 23.1517626 | 61.184 | 874 | 301 | 0 | C* | 11 | 354.5334 | -4.8798 | 2 | 1 | 999.999 | 999.999 | 6 | 0 | 6 | 99999.999 | 99999.999 |
118 | 1995 | 8 | 23.1524561 | 61.184 | 874 | 301 | 1 | 26 | 11 | -2.7989 | 1.3351 | 2 | 1 | -0.020 | 0.001 | 6 | 0 | 6 | -5.638 | -0.481 |
118 | 1995 | 8 | 23.1524561 | 61.184 | 874 | 301 | 1 | 16 | 11 | 63.7366 | -4.3066 | 2 | 1 | -0.003 | -0.076 | 6 | 0 | 6 | 104.309 | 52.014 |
118 | 1995 | 8 | 23.1524561 | 61.184 | 874 | 301 | 1 | 36 | 11 | 67.5644 | -4.1071 | 2 | 1 | 0.058 | -0.052 | 6 | 0 | 6 | 110.151 | 55.859 |
118 | 1995 | 8 | 23.1524561 | 61.184 | 874 | 301 | 1 | C6 | 11 | 33.9379 | 0.8566 | 2 | 1 | 999.999 | 999.999 | 6 | 0 | 6 | 52.653 | 32.667 |
118 | 1995 | 8 | 23.1524561 | 61.184 | 874 | 301 | 0 | C* | 11 | 354.5334 | -4.8798 | 2 | 1 | 999.999 | 999.999 | 6 | 0 | 6 | 99999.999 | 99999.999 |
The catalog is available in electronic form at ftp://ftp.bdl.fr/pub/NSDC/saturn/raw_data/position/. It is also available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/380/727. It contains astrometric coordinates, in J2000 system, of the observations analyzed in Sect. 4. Table 4 shows an extract of it. They are satellite/satellite positions expressed in arcsec, gathered by series with the corresponding scale (in pixels/arcsec) and orientation (in radians). The dates correspond to the mid-time of the pose. This date is not light-time corrected.
The format and the conventions of our catalog are very near to the one of Strugnell
& Taylor (1990) which contains most of the Earth-based observations.
In a FORTRAN code, one line is read with the format:
(i3,i5,i3,f11.7,f7.3,2i4,i2,1x,a2,a1,i2,i1,
2(1x,f13.7),2i2,2(1x,f7.3),3i3,2f10.3).
Mainly, we have changed the format in order to have more digits in the date, in the coordinates and in the O-C. In this way, the precision reaches the mas level. Of course, this precision is only the precision of the record.
series | pixel | angle of orientation to equator |
(arcsec/pix) | (rad) | |
1 | 0.294384053 | 0.004800195 |
2 | 0.547902822 | 0.531254729 |
3 | 0.547682483 | 2.625230525 |
4 | 0.547603319 | 0.529115226 |
6 | 0.547880584 | 0.530138505 |
The meanings of each parameter can be found in (Strugnell & Taylor 1990), but we recall them below with eventually, in parenthesis, the differences with the original meaning:
For any catalog, it is important to know the method of reduction to allow future uses of the observations. It is why, we add at the end of our catalog, some lines which indicate how the reduction have been done and which corrections have been applied to get the coordinates. In particular, we give the scale and the orientation of each series of observations. These values are given in Table 5. Note that the orientation refers to the true equator of the date because the calibration parameters are issued from a direct comparison between the observed frame and a "computed frame''. But as explained in Vienne et al. (2001), it is possible to touch up the scale and the orientation, directly by comparing the observed astrometric coordinates and the computed ones (then given in the J2000 system). For example, a correction of these parameters could be useful to adjust dynamical theory. Another example, a new value for the mass of Saturn can be considered throughout the scale factor.
We give 6006 differential positions of the main Saturnian satellites issued from astrometric measurements of CCD observations performed in 1995 at the Laboratório Nacional de Astrofísica at Itajubá in Brazil. Most of these frames have no reference stars, then we applied an inter-satellite reduction.
We have used the positions of Tethys, Dione, Rhea and Titan given by TASS1.7 to determine the scale factor and the orientation of the receptor.
The comparison of the observed positions and the computed positions shows that
the dispersion of the observations is
.
Furthermore, we point out
the quality and the number of the positions of the satellites not used in the
calibration:
for the 216 positions of Mimas,
for
Hyperion (324),
for Iapetus (524).
The catalog of the corresponding astrometric data could also be accessed under a standard form on the data base NSDC dedicated to the natural satellites ( ftp://ftp.bdl.fr/pub/NSDC/saturn/raw_data/position/). The format and the conventions of our catalog are very near to the one of (Strugnell & Taylor 1990). The coordinates are given apart from a scale factor and from a rotation, but all astrometric corrections are done. So, these differential positions are really astrometric ones in that meaning that, no astrometric consideration is necessary to use them, even if one wants to touch up the calibration. Nevertheless, the raw pixels are also given in order to allow anyone to reduce again the frames.
Acknowledgements
We thank all the colleagues and students who joined the campaign of observation of the Saturnian satellites at Itajubá, especially Marcelo Assafin, Serguei Puliaev, Marildo G. Pereira and Daniela Lazzaro.