A&A 368, 292-297 (2001)
J. T. T. Mäkinen1 - J.-L. Bertaux2 - T. I. Pulkkinen1 - W. Schmidt1 - E. Kyrölä1 - T. Summanen1 - E. Quémerais2 - R. Lallement2
1 - Finnish Meteorological Institute, Geophysics Research, PO Box 503, 00101 Helsinki, Finland
2 - Service d'Aéronomie, BP 3, 91371 Verrières-le-Buisson, France
Received 29 September 2000 / Accepted 8 December 2000
The SWAN instrument onboard the SOHO spacecraft is a Lyman scanning photometer cabable of mapping the whole sky with resolution. Since January 1996 the instrument has produced on average three full sky maps a week with the principal scientific objective of observing the distribution of heliospheric neutral hydrogen. In addition, these systematic observations are a valuable source for studying comets brighter than a visual magnitude of 7-11, the observing limit depending on the abundance ratios of produced radicals and the location of the comet relative to the galactic plane. When the data before the temporary loss of control of SOHO at the end of June 1998 were processed, altogether 18 comets were positively identified, of which one is a new discovery and another 5 can be detected on SWAN images before their actual discovery date. This demonstrates the feasibility of SWAN as an instrument for cometary surveys. The observations are used to estimate the water production rates of the detected comets near their perihelion passages.
Key words: data analysis - surveys - comets - ultraviolet: solar system
The SOHO (Solar and Heliospheric Observatory) spacecraft has already been recognized as the most successful comet finder ever, the total number of discoveries exceeding 200 during less than five years of operation. Almost all new comets belong to the Kreutz sungrazer family and can be detected from LASCO (Large Angle and Spectrometric Coronagraph) images. While most of the SOHO instruments are studying the immediate surroundings of the Sun, one of them covers the rest of the sky. The SWAN (Solar Wind Anisotropies) instrument (Bertaux et al. 1995) onboard SOHO is a Lyman multianode scanning photometer with an instantaneous field of view (FOV) of with 25 pixels each. The instrument consists of two sensor heads each of which has an overall FOV of over steradians covering northern and southern ecliptic hemisphere, respectively. In the normal operation mode the instrument is capable of mapping the whole sky in one day but since the time has to be shared with other kinds of observations, the instrument has produced a complete sky map (Fig. 1) every three days on average. The observing activity has varied over the operational period (Fig. 2). In this respect the period from December 1996 to June 1998 represents the full capability of the instrument. At the end of the period the control of the spacecraft was lost for several months, and data gathered after the recovery show that the instrument was degraded through direct exposure to sunlight during the period of inactivity, decreasing its sensitivity and spectral resolution. Despite additional setbacks like the subsequent loss of attitude control gyroscopes of the spacecraft the observing campaign proceeds at full scale.
|Figure 1: SWAN full sky map in ecliptic coordinates. This resolution UV image recorded April 9, 1996, shows the spatial distribution of interplanetary hydrogen, the area of bright stars near the galactic plane and the comet C/1996 B2 Hyakutake, denoted by the white arrow. The view of the instrument is always obstructed in two directions: pointing the instrument too close to the Sun (s) might damage the sensors, and the antisolar direction (a) is partly obstructed by the spacecraft itself. The locations of the unobservable zones change as the spacecraft moves on its orbit around the Sun|
|Open with DEXTER|
The primary use of the SWAN full sky UV maps is to study the latitudinal distribution of solar wind deductible from asymmetries in the cavity that it carves in the passing cloud of interstellar neutral hydrogen which resonantly scatters solar light (Bertaux et al. 1999b). Another contribution to neutral hydrogen in the solar system comes from the photodissociation of , the major volatile component of cometary nuclei. Several known comets have been observed separately, obtaining valuable results (Bertaux et al. 1998, 1999a; Combi et al. 2000). Since these observations cannot create a complete track record of comets possibly detectable by the instrument, especially new comets only discovered near their perihelion, a cometary survey was undertaken using the full sky images which constitute the most systematic set of measurements with best available coverage.
Since the SWAN instrument was not designed primarily to detect comets, its performance in this respect is far from ideal. The spatial resolution is restricted to the
of sensor pixels although this shortcoming is somewhat compensated by the fact that the hydrogen cloud of a comet is almost spherically symmetric and orders of magnitude larger than the visible dust tail. The point spread function (PSF) of the instrument has a standard deviation comparable to the pixel size near the
line but grows significantly towards the limits of the observing window of 115-170 nm and has a slight dependence on observing geometry. This spreads out the images of hot UV stars on the sky maps. The situation is further complicated by line of sight (LOS) retrieval inaccuracies which contribute another
of random diffusion to the observed signal. These effects together pose great difficulties for observing comets especially near the galactic plane where UV stars are abundant. Also the short exposure time of the normal mode - 13 s in any particular direction - gives only a modest signal to noise ratio despite the relatively high photometric sensitivity of about 0.84 counts per second per Rayleigh per pixel. All these complications are further accentuated on the post-recovery data but the systematic nature of observations combined with very high coverage and the location of the spacecraft at the L1 point between the Earth and the Sun,
unpolluted by the Earth's exospheric
emission still makes the SWAN instrument an important tool for cometary studies.
|Figure 2: Number of produced SWAN full sky maps per month from January 1996 to June 1998|
|Open with DEXTER|
|Figure 3: Detecting comets with neural network. The trail of C/1995 O1 (Hale-Bopp) as denoted by the white arrow, is easily discernible on the combined sky image for the period of May to July 1997 (above). The neural network tracks minute differences in subsequent images and highlights the trails of other comets as well (below). Denoted by arrows are the trails of (a) C/1997 O1 (Tilbrook), (b) 2P/Encke, (c) C/1995 O1, (d) C/1997 N1 (Tabur) and (e) C/1997 K2. The random noise around the data gap is caused by reflections from the spacecraft|
|Open with DEXTER|
The aforementioned shortcomings in mind, a heuristic combined filtering and neural network method (Mäkinen et al. 2000b) was developed to detect all possible candidates for cometary objects. The single exposures of one observing session are binned into a
grid in ecliptic coordinates. A variable-scale median difference filtering removes the smooth background caused by the interstellar neutral hydrogen cloud. The method is similar to the conventional overlaying of an image with its out-of-focus negative to reveal fine structure. The median filter was chosen because it preserves sharp changes in base intensity, e.g., around the borders of a data gap. Since the remaining comets and stars appear identical in shape and because of the mentioned inaccuracies the random deformations of stars prevent direct image subtraction, a 3-layer bidirectional neural network was developed to process several months worth of images in a batch, highlighting consistent motion of any image feature. If Skij is a set of T sky maps where the indices i and j denote longitude and latitude, respectively, and k the observing time, then the nodes in the successive layers A, B and C of the neural network are defined as
The implementation of the network contains several parameters for restricting the evaluation of layer C nodes shortening the processing time considerably. After the layer C is completed the data need to be visualized. Because the node space is multidimensional, the flow is reversed by feeding the obtained weights through established connections back to the layer A from which the results can be read. The output from the neural network (Fig. 3) depicts cometary trails which can be immediately compared to the orbits of known comets. The sensitivity of the neural network has been tested with simulated data and found to be comparable or in some cases even better than visual inspection of time lapse series of filtered images. The largest problem with the algorithm so far has been that the observation times are only known to an accuracy of one day but this problem could be eliminated by improving the data preprocessing tools.
The SWAN images were processed by the neural network in quarterly sets. Altogether 18 comets as listed in Table 1 were identified from the SWAN full sky images from January 1996 to June 1998 of which C/1997 K2 proved to be a new discovery (Mäkinen et al. 2000a). The visibility period of every comet on SWAN images was determined by visually estimating the first and last images where the existence of a comet could be confirmed without a priori knowledge of its location. When the known orbital elements are used in combination with advanced processing methods these limits can be extended to yield valuable data. The list contains 6 short-period comets and 12 long-period comets. When the actual discovery dates of the long-period comets are compared to their visibility on SWAN it can be noticed that half of the long-period comets are visible on SWAN before their discovery. The situation is especially clear with comets C/1997 O1 (Tilbrook) and C/1998 J1 (SOHO) both of which were found near the perihelion but the SWAN instrument has recorded months of prediscovery data.
The visual magnitude of a comet at the last time when it is detectable in SWAN images m1L is given in Table 1 where it is relevant. From the listed values it can be estimated that the limiting magnitude for the SWAN instrument is 7-8 in areas of high star density and 10-11 in voids. Furthermore, the southern hemisphere sensor is less sensitive than the northern one by a factor of about 2.6 and this affects the detection limit. Considering these values, the early disappearance of C/1995 O1 (Hale-Bopp) is quite peculiar but it can be explained by the fact that at the moment the comet was located in a densely populated part of the galactic plane in the direction of Puppis and moving relatively little and it was, at least by visual observation, lost among background star contamination.
|Figure 4: Spatial coverage of all comets detected by SWAN. The trails of detected comets are plotted on ecliptic coordinates with one day stepping.|
|Open with DEXTER|
The water production rate
of a comet is important because it is the most abundant product of cometary activity to which other production ratios are scaled. It can be derived from the observed neutral H distribution by assuming that the photodissociation of
outside the collision sphere is the only noteworthy process. A simple stationary and monokinetic Haser model (Haser 1957; Festou 1981) is used to calculate the H column densities for a
reference value of
The equations follow directly from recursively applying the relation
for the density n where the production term P(r) for daughter products is determined by the destruction rate of parent molecules. Thus for the
In the recent years the search for near-Earth objects (NEOs) has received considerable attention. The existing surveys, however, concentrate primarily on cataloguing all the potentially hazardous asteroids which is arguably a feasible objective since these objects are on orbits which bring them near enough to be observed every few years. The survey coverage is fairly limited which is demonstrated by the fact that amateurs still have a fair chance of discovering a new comet. Comets have much more variation in their orbital parameters and a large part of them may visit the inner solar system just once. Because of their higher kinetic energy and virtual invisibility before the nucleus is activated they pose a direct long term global threat. Collision probability estimates depend on the size distribution of comets which is still not adequately known.
The late discovery of C/1997 K2 and other prediscovery data underline the advantage of an instrument with full sky coverage in detecting new comets. An instrument looking for OH emission as suggested by Brandt et al. (1996a, 1996b) would not be affected by the interstellar neutral hydrogen. On the other hand this is not the largest problem in the SWAN data. The binning method applied in producing full sky maps contributes to the degradation of spatial resolution since it is optimized for large bin sizes. With advanced processing methods higher resolution will be achieved.
The use of the simple Haser model is justifiable in most situations. The largest comets provide ample data so that a more complex model can be used like with C/1995 O1 (Combi et al. 2000). Such models can also benefit from the rudimentary spectral measurement capacity of the instrument given by a H cell filter which can be used to derive the velocity of neutral H atoms. Full sky maps with H cell active are made as well although not as often as ordinary observations. Another shortcoming of the Haser model is apparent with comets whose undergoes rapid fluctuations. A dynamical model could not use the currently available full sky maps since their time resolution is too coarse. Thus, in combination with developing a time-dependent model one should use single exposures directly.
The SWAN full sky maps are very useful for calculating systematic water production rates with some caveats. Besides that the random spreading of stars on SWAN full sky maps introduces uncertainties to the determination of cometary water production rates, other sources of error exist which are not compensated for in this study. The spatial and temporal solar intensity variations are considerable but in principle they could be tracked by observing the apparent background intensity changes at suitable areas over time. Furthermore, the instrument still has some calibration issues which must be addressed before a systematic record can be constructed. Especially the post-recovery data will need considerable calibrating effort. Once these issues have been adequately addressed, more comprehensive reviews can be produced concerning each major data set: the initial full sky observations, the post-recovery full sky observations and the comet specific observations.
SOHO is an international co-operative mission of ESA and NASA. SWAN was financed in France by CNES with support from CNRS and in Finland by TEKES and the Finnish Meteorological Institute. The work of J.T.T.M. and T.S. was supported by the Academy of Finland.