***** File GIOGRE.TXT NOTE: This file was created by scanning the original hardcopy article and only the Figure captions are included. Table 1 was missing from the printed version and is not included. The Giotto Radio- Science Experiment P. Edenhofer Institut fur Hoch- und Hochstfrequenztechnik, University of Bochum, W. Germany M.K. Bird Radioastronomisches Institut, University of Bonn, W. Germany H. Buschert Institut fur Hoch- und Hochstfrequenztechnik, University of Bochum, W. Germany P.B. Esposito Jet Propulsion Laboratory, Pasadena, USA H. Porsche Institut fur Physik der Atmosphare, Deutsche Forschungs- und Versuchsanstalt fur Luft- und Raumfahrt (DFVLR), Oberpfaffenhofen, W. Germany H. Volland Radioastronomisches Institut, University of Bonn, W. Germany Abstract The scientific objectives of the Giotto Radio-Science Experiment (GRE) are to determine both the columnar electron content of Comet Halley's ionosphere and the cometary mass fluence from atmospheric drag by using the radio signals from Giotto during the Halley encounter on 13 - 14 March 1986. The radio-science data (S- and X-band Doppler and range measurements) will be collected at NASA's deep-space 64 m tracking antenna at Tidbinbilla near Canberra, in Australia. In order to separate the effects of the terrestrial ionosphere and the interplanetary plasma, S-band Doppler measurements will also be taken at Tidbinbilla along the line-of-sight of Japan's com- etary probe Sakigake during the Giotto - Halley encounter. The measurements of com- etary electron content and mass fluence will be inverted to derive the spatial distribution of the electron and mass (dust and gas) density within Halley's coma. The GRE is the only experiment on Giotto capable of measuring the low-energy (~ 10 eV) electron bulk population of Halley's ionosphere and the total cometary mass flow im- pacting upon the spacecraft. 1. Scientific Objectives The Giotto spacecraft is equipped with redundant radio transponders capable of downlink transmission at S-band (2.3 GHz) and X-band (8.4 GHz). It is expected to travel through the dayside atmosphere of Comet Halley (Reinhard, 1981), with a relative velocity of 69 km/s. The target distance for the cometary flyby will probably be 500 +/- 100 km. The goal of the GRE (Edenhofer, 1980) is to take measurements of the phase (Doppler) shifts of the downlink carrier signals both at S- and X-band as a function of time during the Halley encounter. The specific scientific objectives are determination of: (a) the columnar electron content of the ionosphere of Comet Halley (b) the cometary mass fluence and total mass flow due to atmospheric drag. During encounter, the carriers of the telemetry signals from the Giotto spacecraft will be subject to: (a) dispersion by cometary plasma influencing the phase and group velocities of signal propagation (strongest at S-band); and (b) Doppler frequency shift induced by the deceleration of the spacecraft due to atmospheric drag (strongest at X- band). Since both effects will be intermingled in the data collected, joint S/X-band data analysis and evaluation will be required to determine both the columnar electron con- tent and the mass fluence. By applying appropriate models of the cometary atmosphere and using integral inversion techniques, the spatial distributions of electrons as well as mass (dust and gas) density will be derived for the comet's coma. The GRE investigations will complement the in-situ particle measurements to be made onboard Giotto. The GRE determination of the total cometary electron content (along the Giotto -Earth line-of-sight) and mass flow (along the spacecraft trajectory) will provide useful reference values for the on-board plasma analyzer (RPA experi- ment) and dust detectors (DID, PIA). In view of the sensitivity threshold of the RPA, only the GRE is capable of measuring the low-energy electrons (~ 10 eV) which con- stitute the bulk of the charged-particle population within Halley's coma. 2. Scientific Background 2.1 Electron content The GRE is an occultation-type, remote-sensing experiment which relies on the pro- pagation effects of electromagnetic waves. Radio-science experiments of this type have been performed in the past to probe the inner and outer coma of the Sun (Edenhofer et al., 1977; Esposito et al., 1980; Volland et al., 1977; Bird, 1982), or the planetary atmospheres (Fjelbo et al., 1975a: Fjelbo et al., 1975b; Levy et al., 1981; Vasilyev et al., 1980), using occulting space probes such as Helios, Pioneer, Mariner, Voyager, Venera, etc. A cometary dual-frequency radio-science experiment is also being undertaken in conjunction with the Vega- 1 and -2 missions (Armand et al., 1982). As a function of time, the ground-received phase phi(t) (in cycles) of a stable, monochromatic radio signal of frequency f is proportional to the (downlink) columnar electron content I(t), i.e. the integrated electron density N_e along the ray path from the spacecraft to the tracking antenna on Earth. The signal's phase will shift in time with respect to that of a local oscillator tuned to the undisturbed (i.e. vacuum wave propagation) Doppler frequency f' = f(1 -(dR(0)/dt)/c) by an amount delta phi(t) = -f(rho(t))/c+40.3/(cf) I(t), where the quantity dR(0)/dt (known from celestial mechanics and orbit deter- mination) is the radial component of the relative velocity between Earth and probe at times close to encounter (t = 0), and the distance rho(t) is associated with the drag- induced deceleration delta v(t) of the probe, via rho(t) = integral(delta v(tau) dtau (c = velocity of light; MKS units; phase constants omitted). This corresponds to a Doppler frequency shift delta f(t) = [-f delta v(t)/c] + [40.3/(cf)*I(t)] (Bird et al., 1982; Porsche & Edenhofer, 1982: Bird et al., 1983; Bird et al., 1984). During the encounter, the first term of this equation is expected to result in a 'blue shift' (delta v<0), caused by the deceleration of the spacecraft in the cometary atmosphere (frequency increase). The second term yields an increase in frequency for cometary plasma entering the line-of-sight (time derivative I>0; immersion), but reverses its sign for maximum electron content and then produces a 'red shift' or frequency decrease, respectively (emersion). Since the Giotto radio subsystem features two downlink carrier signals at frequen- cies f(s) and f(x), two such phase shifts delta phi(s) and delta phi(x) can be measured simultaneously. The plasma and drag effects due to the cometary atmosphere (ionized and neutral part) can be separated unambiguously by a linear combination of these measured dual- frequency phase shifts: delta psi(s) = delta phi(s) - (f(s)/f(x)) delta phi(x) = 40.3/(cf(s)) [1 -(f(s)/f(x))**2] I(t) delta psi(x) = (f(s)/f(x)) delta phi(s) - delta phi(x) = f(x)/c [1 -(f(s)/f(x))**2] rho(t) With just one downlink frequency, plasma and drag effects cannot be separated, even if an uplink is maintained. The only possibility is then a least-squares fit of the record- ed data to model time profiles of the unknown quantities I(t) and rho(t). To determine the deceleration delta v reliably, the ratio between the Doppler frequency shift delta f and the nominal downlink frequency f(0) must be less than about 10**- 10, which is the expected order of magnitude for delta v/c at the time of encounter, e.g. delta psi(x)= f(xo)(delta v/c - delta f/f(so)) [1 - (f(so)/f(xo))**2]. The quantity delta f(t) represents small frequency perturba- tions due to systematic drifts, randomly distributed instabilities, etc. The most important parameter controlling the electron density and therefore the col- umnar electron content in the cometary ionosphere is the total gas production rate Q. It depends, for example, on the size of the cometary nucleus and the distance from the Sun (Reinhard, 1981; Mendis & Ip, 1979). A postulated balance between photo- ionization and losses by recombination governs the physics of cometary plasma, yielding a model representation of the cometary electron density (Mendis & Houpis, 1982) of the form N_c(r) = (A/r) m**-3 (r in meters; A ~ Q**1/2) For Comet Halley, a reasonable value for the quantity Q is taken to be Q(0)= 1 X 10**29 s**-1 at 1 AU (Ip, 1980). Out to about 10**4 km, the electron density is found to be well represented by A = 3 X 10**15 m**-2. Further out, the radial dependence may be somewhat different because of the compressing solar-wind plasma. The outer boundary of the cometary ionosphere is estimated to be about r(c) = 5 X 10**5 km. In a first approximation, the electron density distribution is considered to be spherically symmetric. Thus, in a nucleus-centred coordinate system, the cometary content is denoted by I(c) = Integral(from S(g) to (S(c)) N_c(r) ds = A ln [s(c) + r(c)/(s(g) + r(g))] where s(g)**2 = r(g)**2 - p**2; s(c)**2 = r(c)**2 - p**2, and the quantities r(g), r(c), and p represent the distances to Giotto. to the outer boundary of Halley's ionosphere, and to the Earth-Giotto (impact parameter) ray path, respectively (Edenhofer et al., 1980; Bird et al., 1982; Bird et al., 1983). For a flyby distance of 500 km, Figure 1 shows the temporal variation of I(c) (time profile) to be expected during encounter (S-band). The maximum columnar electron content I(comax) = 3 X 10**16 m**-2 will be obtained about 3.5 s after encounter, where this value may increase to some 7 X 10**16 m**-2, depending on whether or not there are locally enhanced regions from internal sources of additional ionization (e.g. jet-stream configurations) within the cometary coma intersected by the ray path. Figure 1. (a) Columnar electron content of Comet Halley during Giotto encounter for a flyby distance of 500 km at S-band (1 hexem=10**16 m**-2). Models of cometary electron density distribution vary with r**-1 (photo-ionization only) and r**-2, respectively. (b) Electron content I(c) (---) and rate of change dI(c)/dt (---) versus time, corresponding to variations in phase and Doppler frequency shift. A noise level of +/-3x10**14 m**-2 s**-1 (+/-15 mHz) is indicated, as derived from Giotto Doppler residuals collected at the NASA DSS 42 ground station on 31 July 1985 The taking into account of such kinds of phenomena can result in significantly dif- ferent profiles for I(c). For a gas production rate Q = 4 Q0, the maximum electron con- tent will be raised by a factor of two. Computer simulations indicate, however, that I(comax) may vary only about +20% and -9%, if the distance changes to 100 and 1000 km, respectively. For the mission baseline (X-band downlink only during encounter), the electron content would be reduced by a factor of (2.3/8.4)**2 = 1:13, corresponding to a loss in signal strength for I(c) of about 11 dB. Obviously, neither optimistic variations of parameters for encounter geometry nor for cometary plasma physics can balance such large a loss in signal strength (see Sections 3 and 4 for intrinsic noise levels of GRE measurements). Figure 2 shows some details of the geometry of encounter and the resulting com- etary electron content. The temporal behaviour of the electron content during en- counter is expected to be specifically different from the steady-state time signature of those portions of the total electron content measured that are due to the terrestrial ionosphere and interplanetary plasma. Figure 2. Schematic of the geometry of Giotto's Halley encounter on 13/14 March 1986. The cometary electron content represents the integrated electron density of Halley's ionosphere along the ray path C - A, the cometary mass fluence being accumulated from atmospheric drag along the trajectory C - B 2.2 Cometary mass fluence: superposition of plasma and mass effects As Giotto travels through the atmosphere of Comet Halley, dust and gas particles will hit its bumper shield. Cometary dust and gas particles typically consist of grain materials such as magnetite, graphite, quartz, water ice (grain radii roughly 10**-7- l0**-2 m; mass spectrum 10**-17- 10**-3 kg) and sublimation products of nuclear constituents in terms of molecules such as H2O, CO2, NH3 (e.g. mean thermal speed 400 m/s, outflow speed at surface of nucleus 180 m/s) (Newburn, 1981; Divine, 19,81; Hubner & Keady, 1982; Hellmich & Keller, 1981; Schwehm & Kneissel, 1981). The impacts of these particles build up an atmospheric drag effect, causing a deceleration delta v of the spacecraft (Edenhofer et al., 1980; Porsche & Edenhofer, 1982; Bird et al., 1983). This deceleration can be estimated, from momentum and energy- transfer relationships, to be delta v = m(0)*v(0)*(1 + q**1/2 cos(phi)) / M where m(0) and v(0) are the mass and the velocity of the impacting (recoiling) particles, M is the mass of the spacecraft (512 kg during encounter) and phi is the angle of intersection between the velocity vectors v and - v(0). Here inelastic energy losses delta Q (e.g. heat generation and ionization processes) are taken to be negligible, according to delta Q << m(0)v(0)/2. Since the quantity delta v can be measured from the observables of phase (and group) delay shift (e.g. delta psi(x)) via rho(t), the cometary mass fluence is finally derivable as the mass flow integrated along the trajectory of the spacecraft. Multipli- cation by the corresponding spacecraft area (2.64 m**2) then yields the total impacting cometary mass. The mass-fluence determination is complicated by the fact that the ratio q of the mass of the recoiling particles to that of the impacting particles is also a function of the energy of the impacting particles (Hirao, 1982). High-energy gas molecules in the cometary coma, for instance, may generate multiple secondary ions. At the encounter velocity of v(0)=69 km/s, energy levels higher than 1 keV apply for molecules of atomic weight as low as 42. As far as dust particles are concerned, the production rate for emission of secondary ions is expected to be even higher than in the case of gas particles. In laboratory experiments, q-values as high as 40 have been observed for velocities much lower than v(0) (Porsche & Edenhofer, 1982). Following refined calculations of the expected dust-particle fluence based on data from the last apparition of Comet Halley in 1910 (Newburn, 1981), and on the com- putations by Hubner & Keady (1982) for the cometary plasma mass which will also impact on the bumper shield of the spacecraft, the total integrated mass flow due to Halley's atmosphere hitting Giotto during the encounter is estimated to be of the order of m(c)=9 X 10**-5 kg. Allowing for an uncertainty of 1