Comparative study of the dust emission of 19P/Borrelly (Deep space 1) and 1P/Halley

Images obtained by the Miniature Integrated Camera and Imaging Spectrometer (MICAS) experiment onboard the Deep Space 1 spacecraft which encountered comet 19P/Borrelly on September 22nd 2001 show a dust coma dominated by jets. In particular a major collimated dust jet on the sunward side of the nucleus was observed. Our approach to analyse these features is to integrate the observed intensity in concentric envelopes around the nucleus. The same procedures has been used on the Halley Multicolour Camera images of comet 1P/Halley acquired on March 14th 1986. We are able to show that at Borrelly the dust brightness dependence as a function of radial distance is different to that of Halley. At large distances both comets show constant values as the size of the concentric envelopes increases (as one would expect for force free radial outflow). For Halley the integral decreases as one gets closer to the nucleus. Borrelly shows opposite behaviour. The main cause for Halley's intensity distribution is either high optical thickness or particle fragmentation. For Borrelly, we have constructed a simple model of the brightness distribution near the nucleus. This indicates that the influence of deviations from point source geometry is insufficient to explain the observed steepening of the intensity profile close to the nucleus. Dust acceleration or fragmentation into submicron particles appear to be required. We also estimate the dust production rate of Borrelly with respect to Halley and compare their dust to gas ratios.     

A quantitative model for comet nucleus topography

Topography of the surface of a comet nucleus is likely rough at all scales smaller than the mean effective radius. We present a flexible and easily scalable model for quantitative calculations simulating the effects of comet nucleus topography on gas release and dust mantle evolution. The topographic features we describe must be large enough (typically> 10 m) so that they will not erode in one orbit of the nucleus around the Sun. The maximum effective size of a hill is about 1/√2 times the effective radius of the nucleus. If it is larger, then an ellipsoidal shape of the nucleus is more appropriate. The procedure described here also permits for inhomogeneous composition of the topographic features, leading to locally different rates of gas production (e.g., jet-like features and filaments) or different thicknesses of the dust mantle. It also can give rise to different temperature patches, locally varying albedos and emissivities, and may explain the formation of permanent dust mantles.     

Modelling of cometary nuclei: Planetary missions preparation

The planning of planetary missions requires the knowledge as much as possible accurate of the target. This knowledge, which cannot be obtained only by ground based observations, can be supported by theoretical modeling. This is particularly true in the case of Rosetta, a cornerstone ESA mission that will be launched in January 2003, and of its target, the comet 46P/Wirtanen. In this paper we show how, using a nucleus thermal evolution model, it is possible to foresee the temperature of the nucleus surface and the activity level along the orbit and define the range of possible values for these properties. Activity level close to the aphelion depends on the presence of CO. Surface temperatures depend on the physical properties attributed to the dust: temperatures reached by a crust mainly composed by silicatic grains differ from the temperatures reached by a crust mainly composed by organic grains.     

The derivation of Halley parameters from observations

A set of nominal model parameters for P/Halley is derived from its light curve and spectra. In those cases where Halley observations are not sufficient, the average value derived from a large set of other comets has been used, or data from comet Bennett, Halley's best analogue has been taken. The derived parameters include nucleus mass, size, density, albedo, rotation period, axial inclination, and surface temperature, the composition of the parent molecules, the total gas and dust production rates, distributions for the dust size and bulk density as well as various other parameters.     

Coma morphology and dust emission pattern of comet Halley

This contribution starts with a short overview on cometary dust modelling and then focuses on the application of coma modelling with respect to in-situ measurements of cometary dust and ground based observations. The fountain model, valid for the dynamics of small cometary dust particles, is discussed. Models using Keplarian theory for the motion of the dust particles are outlined and the ESOC coma model is presented. Some direct applications of this model to analyse the results of the recent spacecraft flybys of comet Halley, as dust flux profiles, particle ground tracks and envelope positions, are shown. To compare the model with ground-based astronomical observations, the utilization of the ESOC coma model for the generation of synthetic images is demonstrated and some future prospects of this technique are outlined.     

Nucleus studies of comet Halley

Early results from, and research initiatives warranted by, the Earth-based observations of Halley's near-nucleus and related phenomena are reviewed. Where appropriate, this information is combined with spacecraft data obtained by the various flight projects. The basic objective is to gain a greater insight into the nature of the comet's nucleus and its environment. Among topics are the brightness variations at large heliocentric distances along the inbound leg of the orbit; the bulk and rotational properties of the nucleus, including possible precession; the surface morphology and the formation of dust jets; subfemtogram dust particles and their presence in a sunward spike and relation to CN jets; comparison of the hydrogen coma's “pulsation” pattern with the surface distribution of major dust vents; the events causing Giotto's wobbling near its closest approach to the comet; and the recent developments in theoretical modeling of the icy-conglomerate nucleus.     

Numerical simulations of the radiance from the comet 46P/Wirtanen in the various configurations of the measurements during “Rosetta” mission

The work we present deals with the spectrometric measurements of VIRTIS instrument of the Comet P/Wirtanen planned for the Rosetta mission. This spectrometer can monitor (VIRTIS M channel: 0.250μm – 0.980μm; Δκ=20cm⁻¹; 0.980 – 5.0 μm; Δκ = 5cm⁻¹; VIRTIS H channel: 2.0 μm – 5.0 μm; Δκ=2cm⁻¹) the nucleus and the coma in order to provide a general picture of coma's composition, the production of gas and dust, the relationship of coma production to surface composition and the structure and variation of mineralogy of the nucleus surface. During the mission the observation conditions of the spectroscopic investigation change due to different relative positions spacecraft/comet, and to the different illumination conditions of the surface at various distances of the comet to the Sun. The nucleus surface is continuously modified by the ice sublimation accompanied by gas and dust emission. Consequently the surface also its spectrophotometric properties changes and their monitoring can give a new insight. The important role of simulations is to predict the results of measurements in various experimental condition what, in the future, can help in interpretation of the measured data.

In this paper the first results of our simulation the radiance from the comet in the 0.25–5.0μm spectral range at two distances from the Sun (1AU and 3AU) are shown. The distance between the Rosetta orbiter and the nucleus surface as well as the sun zenith angles are taken into account according to the Rosetta mission phases. In fact the surface and coma properties vary along the comet orbit, and should be taken into account in our calculations. The optical parameters of the dust on the surface (e.g. reflectance) and in the coma (e.g. Q_ext) were calculated from optical constants of possible comet analogues. The thermodynamic parameters of the comet are taken from the models of comet evolution. Through this kind of modelling it is possible to identify the surface characteristics in spectra of the radiation from the surface of nucleus transmitted through the coma loaded with dust and gases.

Even if the “Rosetta mission” is postponed, with the consequence of a target change, we think that our idea and the method used for the simulations can be useful also for the new Rosetta target - the comet 67P/Churyumov Gerasimenko.     

A first step towards proton flux forecasting

We present a preliminary version of a potential tool for real time proton flux prediction which provides proton flux profiles and cumulative fluence profiles at 0.5 and 2 MeV of solar energetic particle events, from their onset up to the arrival of the interplanetary shock at the spacecraft position (located at 1 or 0.4 AU). Based on the proton transportation model by Lario et al. [Lario, D., Sanahuja, B., Heras, A.M. Energetic particle events: efficiency of interplanetary shocks as 50 keV E < 100 MeV proton accelerators. Astrophys. J. 509, 415–434, 1998] and the magnetohydrodynamic shock propagation model of Wu et al. [Wu, S.T., Dryer, M., Han, S.M. Non-planar MHD model for solar flare-generated disturbances in the Heliospheric equatorial plane. Sol. Phys. 84, 395–418, 1983], we have generated a database containing “synthetic” profiles of the proton fluxes and cumulative fluences of 384 solar energetic particle events. We are currently validating the applicability of this code for space weather forecasting by comparing the resulting “synthetic” flux profiles with those of several real events.     

A fine mist of very small comet dust particles

A comet nucleus considered as an aggregate of interstellar dust would produce a mist of very finely divided (radius ～ 0.01 μm) particles of carbon and metal oxides accompanying the larger dust grains. These small particles which are very abundant in the interstellar dust size spectrum would provide substantial physical effects because of their large surface area. They may show up strongly in particle detectors on the Halley probes. A strong basis for serious consideration of these particles comes from the other evidence that interstellar dust grains are the building blocks of comets; e.g. (1) the explanation of the “missing” carbon in comets; (2) The S₂ molecule detection which suggests that the comet solid ice materials have been previously subjected to ultraviolet radiation (as are interstellar grains) before aggregation into the comet; (3) the predicted dust to gas ratio.     

Applications of PVDF dust sensor systems in space

The advantages of ruggedness, no bias requirement, ease of large area sensor construction, high counting rate capability, and space reliability inherent in the Polyvinylidene Fluoride (PVDF) dust sensors which have been under development at the University of Chicago over the last decade have led to PVDF flux/mass/velocity/trajectory systems which have advantages over other systems and are well suited for a variety of dust studies in space. The thermal stability characteristics and flux/mass/velocity/trajectory determining characteristics of PVDF and Vinylidene Fluoride/Trifluoroethylene (PVDF copolymer) dust sensors are described. We summarize the objectives and designs of our earlier VEGA- comet Halley instruments, a PVDF velocity/trajectory dust instrument for launch on the Advanced Research and Global Observation Satellite (ARGOS) in January 1996, and a PVDF high flux dust instrument for launch on the CASSINI spacecraft to Saturn in October 1997.     

Thermal modeling of Halley's comet

The comet thermal model of Weissman and Kieffer is used to calculate gas production rates and other parameters for the 1986 perihelion passage of Halley's Comet. Gas production estimates are very close to revised pre-perihelion estimates by Newburn based on 1910 observations of Halley; the increase in observed gas production post-perihelion may be explained by a variety of factors. The energy contribution from multiply scattered sunlight and thermal emission by coma dust increases the total energy reaching the Halley nucleus at perihelion by a factor of 2.4. The high obliquity of the Halley nucleus found by Sekanina and Larson may help to explain the asymmetry in Halley's gas production rates around perihelion.     

Comet Halley: Nucleus and jets (results of the VEGA mission)

The VEGA-1 and VEGA-2 spacecraft made their closest approach to Comet Halley on 6 and 9 March, respectively. In this paper those results of the onboard imaging experiment which were obtained around closest approach are discussed. The nucleus of the comet was clearly identifiable as an irregularly shaped object, with overall dimensions of (16±1)×(8±1)×(8±1) km. The nucleus rotates in the prograde sense about an axis nearly perpendicular to the orbital plane with a period of 53±2 hours. Its albedo is only $\text{0.04±}\text{0.02}\text{0.01}$ Many of the jet features observed during the second fly-by have been spatially reconstructed. Their sources form a quasi-linear structure on the surface. The dust above the surface is shown to be generally optically thin with the exception of certain specific dust jets. Brightness features on the surface are clearly seen. Correlating our data with other measurements, we conclude that the dirty snow-ball model will probably need to be revised.     

A theoretical model for the magnetic helicity of solar active regions

Active regions on the solar surface are known to possess magnetic helicity, which is predominantly negative in the northern hemisphere and positive in the southern hemisphere. Choudhuri et al. [Choudhuri, A.R. On the connection between mean field dynamo theory and flux tubes. Solar Phys. 215, 31–55, 2003] proposed that the magnetic helicity arises due to the wrapping up of the poloidal field of the convection zone around rising flux tubes which form active regions. Choudhuri [Choudhuri, A.R., Chatterjee, P., Nandy, D. Helicity of solar active regions from a dynamo model. ApJ 615, L57–L60, 2004] used this idea to calculate magnetic helicity from their solar dynamo model. Apart from getting broad agreements with observational data, they also predict that the hemispheric helicity rule may be violated at the beginning of a solar cycle. Chatterjee et al. [Chatterjee, P., Choudhuri, A.R., Petrovay, K. Development of twist in an emerging magnetic flux tube by poloidal field accretion. A&A 449, 781–789, 2006] study the penetration of the wrapped poloidal field into the rising flux tube due to turbulent diffusion using a simple 1-d model. They find that the extent of penetration of the wrapped field will depend on how weak the magnetic field inside the rising flux tube becomes before its emergence. They conclude that more detailed observational data will throw light on the physical conditions of flux tubes just before their emergence to the photosphere.     

Impact induced plasma during a cometary fly-by

GIOTTO, the probe which is presently developed by the European Space Agency, will encounter comet Halley in March 1986 with a relative velocity of 69 km/s. The fore section of the surface will be submitted to the bombardment of dust grains and neutral molecules in the final phase of the mission, like that of an Earth orbiter during atmospheric re-entry. These particles have a kinetic energy of 24 eV per a.m.u.; they produce secondary ions and electrons which form a plasma cloud around the body and control the electric potential of its surface. This paper is a review of the work which has been performed on the subject by dedicated study groups; the purpose of their action was to gather information and produce new findings which might have an influence on the design of the spacecraft and help in the interpretation of the data collected by the scientific payload.

The effect of impact induced plasma may already be significant at 10⁵ km from the comet nucleus; at a distance of 1000 km the flux of ions and electrons produced by cometary dust and neutrals will possibly exceed that of the ambient plasma by more than three orders of magnitude. It is expected that the spacecraft surface potential will be positive and will reach at least a few tens of volts; coating the leading surface of the spacecraft with a thin layer of gold or silver will help reducing the emission of ions from neutral gas. Computer simulation models are used to predict the structure of the charged particle density distribution in the vicinity of the surface. Effects associated with the wake and differential charging are also discussed. The significance of these results is conditioned by the validity of the models and the largest source of uncertainty seems to be associated with the plasma generated by dust impact.     

The Halley dust model

The properties of dust ejecta from Comet Halley are studied on the basis of (a) evidence from the comet's past apparitions and (b) analogy with recent, physically similar comets. Specifically discussed are the light curve and spectrum, discrete phenomena in the head, the physical properties of the nucleus (size, albedo, rotation, surface temperature, and morphology), and an interaction between the nucleus and dust atmosphere. Also reviewed are constraints on the size and mass distributions of dust particles, information on submicron-size and submillimeter-size grains from the comet's dust tail and antitail, and the apparent existence of more than one particle type. Similarities between the jet patterns of Halley and the parent comet of the Perseid meteor stream are depicted, and effects of the surface heterogeneity (discrete active regions) on the dust flow are assessed. Current dust models for Halley are summarized and the existence of short-term variations in the dust content in the comet's atmosphere is suggested.     

Thermal emission from cometary dust

The thermal emission from the dust coma of a comet can be analyzed to yield the flux and size distribution of the dust grains and the relative abundance of silicate and absorbing grains.     

Cometary constraints on the planet forming environment.

Molecular elemental and isotopic abundances of comets provide sensitive diagnostics for models of the primitive solar nebula. New measurements of the N2, NH and NH2 abundances in comets together with the in situ Giotto mass spectrometer and dust analyzer data provide new constraints for models of the comet forming environment in the solar nebula. An inventory of nitrogen-containing species in comet Halley indicates that NH3 and CN are the dominant N carriers observed in the coma gas. The elemental nitrogen abundance in the gas component of the coma is found to be depleted by a factor approximately 75 relative to the solar photosphere. Combined with the Giotto dust analyzer results for the coma dust component, we find for comet Halley Ngas + dust approximately 1/6 the solar value. The measurement of the CN carbon isotope ratio from the bulk coma gas and dust in comet Halley indicates a significantly lower value, 12C/13C = 65 +/- 9 than the solar system value of 89 +/- 2. Because the dominant CN carrier species in comets remains unidentified, it is not yet possible to attribute the low isotope ratio predominantly to the bulk gas or dust components. The large chemical and isotopic inhomogeneities discovered in the Halley dust particles on 1 mu scales are indicative of preserved circumstellar grains which survived processing in the interstellar clouds, and may be related to the presolar silicon carbide, diamond and graphite grains recently discovered in carbonaceous chondrites. Less than 0.1% of the bulk mass in the primitive meteorites studied consists of these cosmically important grains. A larger mass fraction (approximately 5%) of chemically heterogeneous organic grains is found in the nucleus of comet Halley. The isotopic anomalies discovered in the PUMA 1 Giotto data in comet Halley are probably also attributable to preserved circumstellar grains. Thus the extent of grain processing in the interstellar environment is much less than predicted by interstellar grain models, and a significant fraction of comet nuclei (approximately 5%) may be in the form of preserved circumstellar matter. Comet nuclei probably formed in much more benign environments than primitive meteorites.     

A European probe to comet Halley

A European probe to comet Halley is proposed. The probe's model payload consists of 8 scientific instruments, viz. neutral, ion and dust impact mass spectrometers, magnetometer, medium energy ion and electron analyzer, camera, dust impact detectors and plasma wave experiment. Fly-by of the comet Halley nucleus will take place on November 28th, 1985, at about 500 km miss distance. The main spacecraft serves as relay link to transmit the observed data to Earth. As probe, a modified ISEE 2 design is proposed. Because of the cometary dust hazard expected in the coma a heavy dust shield (27 kg) is required, consisting of a thin front sheet and a 3 layer rear sheet. The probe is spin-stabilized (12 rpm), has no active attitude and orbit control capability and uses battery power only to provide about 1000 Wh for a measuring phase. A despun antenna transmits up to 20 kbit/s, in X-band. The total probe mass is estimated at 250 kg. The 3 model development programme should start in mid 1981 with Phase B.     

The Giotto mission to Halley''s comet

ESA's Giotto mission to Halley's comet is a fast flyby in March 1986, about four weeks after the comet's perihelion passage when it is most active. The scientific payload comprises 10 experiments with a total mass of about 60 kg: a camera for imaging the comet nucleus, three mass spectrometers for analysis of the elemental and isotopic composition of the cometary gas and dust environment, various dust impact detectors, a photopolarimeter for measurements of the coma brightness, and a set of plasma instruments for studies of the solar wind/comet interaction. In view of the high flyby velocity of 68 km/s the experiment active time is very short (only 4 hours) and all data are transmitted back to Earth in real time at a rate of 40 kbps. The Giotto spacecraft is spin-stabilised with a despun high gain parabolic dish antenna inclined at 44.3° to point at the Earth during the encounter while a specially designed dual-sheet bumper shield at the other end protects the spacecraft from being destroyed by hypervelocity dust impacts. The mission will probably end near the point of closest approach to the nucleus when the spacecraft attitude will be severely perturbed by impacting dust particles leading to a loss of the telecommunications link.