The Coma of Comet 9P/Tempel 1

As comet 9P/Tempel 1 approaches the Sun in 2004–2005, a temporary atmosphere, or “coma,” will form, composed of molecules and dust expelled from the nucleus as its component icy volatiles sublimate. Driven mainly by water ice sublimation at surface temperatures T > 200 K, this coma is a gravitationally unbound atmosphere in free adiabatic expansion. Near the nucleus (≤ 10² km), it is in collisional equilibrium, at larger distances (≥10⁴ km) it is in free molecular flow. Ultimately the coma components are swept into the comet’s plasma and dust tails or simply dissipate into interplanetary space. Clues to the nature of the cometary nucleus are contained in the chemistry and physics of the coma, as well as with its variability with time, orbital position, and heliocentric distance. The DI instrument payload includes CCD cameras with broadband filters covering the optical spectrum, allowing for sensitive measurement of dust in the comet’s coma, and a number of narrowband filters for studying the spatial distribution of several gas species. DI also carries the first near-infrared spectrometer to a comet flyby since the VEGA mission to Halley in 1986. This spectrograph will allow detection of gas emission lines from the coma in unprecedented detail. Here we discuss the current state of understanding of the 9P/Tempel 1 coma, our expectations for the measurements DI will obtain, and the predicted hazards that the coma presents for the spacecraft. An erratum to this article is available at .     

The History and Dynamics of Comet 9P/Tempel 1

Since its discovery in 1867, periodic comet 9P/Tempel 1 has been observed at 10 returns to perihelion, including all its returns since 1967. The observations for the seven apparitions beginning in 1967 have been fit with an orbit that includes only radial and transverse nongravitational accelerations that model the rocket-like thrusting introduced by the outgassing of the cometary nucleus. The successful nongravitational acceleration model did not assume any change in the comet’s ability to outgas from one apparition to the next and the outgassing was assumed to reach a maximum at perihelion. The success of this model over the 1967–2003 interval suggests that the comet’s spin axis is currently stable. Rough calculations suggest that the collision of the impactor released by the Deep Impact spacecraft will not provide a noticeable perturbation on the comet’s orbit nor will any new vent that is opened as a result of the impact provide a noticeable change in the comet’s nongravitational acceleration history. The observing geometries prior to, and during, the impact will allow extensive Earth based observations to complement the in situ observations from the impactor and flyby spacecraft.     

The Information Content of Synthetic Aperture Radar Images of Terrain

A statistical model is developed that portrays an imaging radar as a noisy communication channel with multiplicative noise, and the model is used to evaluate the average amount of information that can be extracted about a target from its radar image. The average information content is also used to define a measure of radiometric resolution for radar images. It is shown that the information content and the resolution capabilities of an imaging radar reach a limit beyond which an increase in scene dynamic range does not improve the information content or the resolution. This limitation results from the multiplicative nature of the noise introduced in the imaging process.     

A Comparison of Numerical Methods for Solving Optimal Control Problems

Linear systems with quadratic performance criteria typically lead to two-point boundary value problems. The numerical solution of such problems is often difficult because of roundoff and truncation errors. In this paper numerical methods for digital computer solutions are considered for a system described by a given second-order linear differential equation and the results are compared. The results show that where accurate analytic solutions cannot be obtained, the invariant imbedding technique yields more accurate results, for a given integration grid interval, than the other methods considered.     

MHD ANALYSIS OF PETSCHEK-TYPE RECONNECTION IN NON-UNIFORM FIELD AND FLOW GEOMETRIES

Analytical studies of reconnection have, for the most part, been confined to steady and uniform current sheet geometries. In contrast to these implifications, natural phenomena associated with the presence of current sheets indicate highly non-uniform structure and time-varying behaviour. Examples include the violent outbursts of energy on the Sun known as solar flares, and magnetospheric phenomena such as flux transfer events, plasmoids, and auroral activity. Unlike the theoretical models, reconnection therefore occurs in a highly dynamic and structured plasma environment. In this article we review the mathematical tools and techniques which are available to formulate models capable of describing the effects of reconnection in such situations. We confine attention to variants of the reconnection model first discussed by Petschek in the 1960s, in view of its successful application in predicting and interpreting phenomena in the terrestrial magnetosphere. The analysis of Petschek-type reconnection is based on the equations of ideal magnetohydrodynamics (MHD), which describe the large-scale behaviour of the magnetic field and plasma flow outside the diffusion region, which we determine as a localised part of the current sheet in which reconnection is initiated. The approach we adopt here is to transform the MHD equations into a Lagrangian or so-called 'frozen-in' coordinate system. In this coordinate system, the equation of motion transforms into a set of coupled nonlinear equations, in which the presence of inhomogeneous magnetic fields and/or plasma flows gives rise to a term similar to that which appears in the study of the ordinary string equation in a non-homogeneous medium. As demonstrated here, this approach not only clarifies and highlights the effects of such non-uniformities, it also simplifies the solution of the original set of MHD equations. In particular, this is true for those types of problem in which the total pressure can be considered as a known quantity from the outset. To illustrate the method, we solve several 2D problems involving magnetic field and flow non-uniformities: reconnection in a stagnation-point flow geometry with antiparallel magnetic fields; reconnection in a Y-type magnetic field geometry with and without velocity shear across the current sheet; and reconnection in a force-free magnetic field geometry with field lines of the form xy = const. These case examples, chosen for their tractability, each incorporate some aspects of the field and flow geomtries encountered in solar-terrestrial applications, and they provide a starting point for further analytical as well as numerical studies of reconnection.     

Resonant inverters for photovoltaic array to utility interface

A number of high-frequency link power converter schemes suitable for interconnecting a photovoltaic (PV) array to the utility line using resonant inverters are presented. Use of high-frequency resonant inverters (HFRI) results in reduced size, weight and cost. Scheme 1 uses the HFRI-cycloconverter, scheme 2 uses the HFRI-rectifier pulsewidth-modulated voltage source inverter, scheme 3 uses the HFRI-rectifier line commutated inverter working with a fixed maximum possible inverting fringe angle, and scheme 4 uses the HFRI-rectifier line-connected inverter (with line current modulation). These schemes are compared. Based on the comparison, scheme 3 was selected for an initial study. To overcome various problems associated with scheme 3, the line current is modulated sinusoidally in scheme 4. In this scheme, the HFRI is controlled appropriately to obtain a rectified sine wave of 120 Hz in the DC link and the function of the line-connected inverter is to unfold this waveform on alternate half cycles to obtain a sinusoidal line current with the line voltage. Methods of commutation of the line connected inverter are explained. Experimental results obtained with scheme 4 are presented     

Spatio-temporal pattern recognition using hidden Markov models

A spatio-temporal method for identifying objects contained in an image sequence is presented. The Hidden Markov Model (HMM) technique is used as the classification algorithm, making classification decisions based on a spatio-temporal sequence of observed object features. A five class problem is considered. Classification accuracies of 100% and 99.7%, are obtained for sequences of images generated over two separate regions of viewing positions. HMMs trained on image sequences of the objects moving in opposite directions showed a 98.1% successful classification rate by class and direction of movement. The HMM technique proved robust to image corruption with additive correlated noise and had a higher accuracy than a single-look nearest neighbor method. A real image sequence of one of the objects used was successfully recognized with the HMMs trained on synthetic data. This study shows the temporal changes that observed feature vectors undergo due to object motion hold information that can yield superior classification accuracy when compared with single-frame techniques     

Remote Sensing of H from Ulysses and Galileo

We model interplanetary H Lyman-α (Lα) observations from Galileo UVS (Ultraviolet Spectrometer) and EUVS (Extreme Ultraviolet Spectrometer) (Hord et al., 1992) and the Ulysses interstellar neutral gas (GAS) instrument (Witte et al., 1992). EUVS measurements near solar maximum (max) in 1990–1992 have a peaked brightness maximum upwind due to a rather isotropic solar wind charge-exchange ionization pattern (A=0–0.25). GAS measurements from solar minimum (min) in 1997 have a plateau in the upwind direction that we model using Ulysses SWOOPS (solar wind plasma experiment) solar min data on solar wind density and velocity at different heliographic latitudes. The isotropic ionization pattern deduced from EUVS at solar max may be consistent with recent SWOOPS results (McComas et al., 2000b, c) that high speed solar wind is absent at high latitudes during solar max. Galileo and Ulysses Lα data favor higher H temperatures (15 000–18 000 K) than previous models. This revised version was published online in August 2006 with corrections to the Cover Date.     

Proton phase space densities (0.5eV<Ep<5MeV) at midlatitudes from Ulysses SWICS/HI-SCALE measurements

Proton phase space densities in the solar wind frame from suprathermal velocities 10 km s^–1 to 30,000 km s^–1 (0.5 eV–5 MeV) were derived from combined SWICS and HISCALE measurements when Ulysses was at 5 AU and –24° heliolatitude. The period (19–23 January 1993) encompasses a forward/reverse shock pair (20 January, 0500 UT and 22 January, 0300 UT). Strong evidence is found for shock acceleration of pickup protons from interstellar hydrogen at all energies measured.     

THE CLUSTER ION SPECTROMETRY (CIS) EXPERIMENT

The Cluster Ion Spectrometry (CIS) experiment is a comprehensive ionic plasma spectrometry package on-board the four Cluster spacecraft capable of obtaining full three-dimensional ion distributions with good time resolution (one spacecraft spin) with mass per charge composition determination. The requirements to cover the scientific objectives cannot be met with a single instrument. The CIS package therefore consists of two different instruments, a Hot Ion Analyser (HIA) and a time-of-flight ion COmposition and DIstribution Function analyser (CODIF), plus a sophisticated dual-processor-based instrument-control and Data-Processing System (DPS), which permits extensive on-board data-processing. Both analysers use symmetric optics resulting in continuous, uniform, and well-characterised phase space coverage. CODIF measures the distributions of the major ions (H⁺, He⁺, He⁺⁺, and O⁺) with energies from ~0 to 40 keV/e with medium (22.5°) angular resolution and two different sensitivities. HIA does not offer mass resolution but, also having two different sensitivities, increases the dynamic range, and has an angular resolution capability (5.6° × 5.6°) adequate for ion-beam and solar-wind measurements.