A strategy for investigation of the outer solar system

The requirements of systematic exploration of the outer solar system have been intensively studied by a Science Advisory Group (SAG) of consulting scientists for the National Aeronautics and Space Administration (NASA). Comets and Asteroids were excluded from this study, as a separate group is planning missions to these bodies. This paper and accompanying articles on specific related scientific subjects written by members of the SAG, summarize the findings and recommendations of this group. These recommendations should not be interpreted as official NASA policy. Following some general introductory remarks, a brief sketch is given of the development and current status of scientific missions to the inner planets by the U.S. and the U.S.S.R. With this perspective, the development of the U.S. program for investigation of the outer solar system is described. The scientific focus of outer solar system exploration has been studied in detail. The relationship of the outer planetary bodies to one another and to the inner planets, as parts in a unified solar system evolved from a primitive solar nebula, is emphasized. Deductions from outer solar system investigations regarding the conditions of the solar nebula at the time of planetary formation have been considered. Investigations have been proposed that are relevant to studies of the atmospheric structure and dynamics, internal structure of the planets, satellite composition and morphology, and planetary and interplanetary fields and energetic particles. The mission type and sequence required to conduct a systematic exploration of the outer solar system has been developed. Technological rationales for the suggested missions are discussed in general terms. The existing NASA program for outer solar system exploration is comprised of four missions:

1. Pioneer 10 fly-by mission to Jupiter and beyond, currently underway, with launch on 3 March 1972;
2. Pioneer G, intended for a similar mission with planned launch 2–22 April 1973; and
3. Two Mariner Jupiter/Saturn fly-bys in 1977, with experiment selection scheduled for late 1972 and detailed engineering design during 1972–74.

The Science Advisory Group advocates that detailed mission planning be undertaken on the following additional missions for launches during the late 1970's and early 1980's. Together with existing plans, these would provide a balanced, effective exploration program.

1. 1976 Pioneer Jupiter/Out-of-Ecliptic (One Mission)
2. 1979 Mariner Jupiter/Uranus Fly-bys (Two Missions)
3. 1979 Pioneer Entry Probe to Saturn 1980 Pioneer Entry Probe to Uranus via Saturn Fly-by (Three Missions)
4. 1981/1982 Mariner Jupiter Orbiter (Two Missions).

     

Solar wind interaction with comets

The planned missions to Comet Halley, which will arrive at the nearest space of the Sun in 1986, have recently revived interest in studying solar wind interaction with comets. Several unsolved problems exist and the most urgent of them are as follows:

1. The character of the solar wind interaction with comets: bow shocks and contact surface formation near comets; similarities and differences of solar- wind interaction with comets and with Venus. The differences are probably associated with a great extension of neutral atmospheres of comets (due to a practical lack of cometary gravitation) and the ‘loading’ of the solar wind flux by cometary ions during the interaction.
2. The anomalous ionization in cometary heads.
3. The problem of the anamalously high accelerations of ions in the plasma tails of comets.
4. The variability of plasma structures observed in cometary tails.

     

Satellite observations of power line harmonic radiation

In-situ spectral observations of power-line harmonic radiation (PLHR) are still quite rare and almost all the detailed characteristics have been derived from studies at Antarctic stations such as Siple and Halley, and their conjugates in North America. Because of the lack of more direct satellite evidence of PLHR and the difficulties in interpretation of morphological studies, such as those of Ariel 3 and 4, there is considerable controversy concerning the relative importance of PLHR and its contribution to wave-particle interactions (WPI) in the magnetosphere. The early Ariel 3 and 4 global surveys indicated that, in terms of true mean wave energy, there is no longitudinal localisation, the contribution of world-wide intense VLF emissions, associated with magnetic storms, being dominant. Also, the most intense wave emission, that of plasmaspheric hiss at ELF (< 1 kHz) exhibits little evidence of localisation. The PLHR phenomenon is most conspicuous by its persistence in quiet times (Kp ≤ 2+) at 45° < Λ < 55° over North America and its conjugate region, even though the less frequent strongest emissions, to which it gives rise in the summer, are located polewards at 3 < L < 5. In the northern winter, when spheric activity over both North America and its conjugate are low, there is a high occurrence of strong discrete emissions, which are more sharply localised than in the summer, on the NE industrial U.S.A. field line. The most recent Ariel 4 studies, particularly on the spheric wavefield over North America (using data from the Appleton Laboratory impulse counters) and on the character of the wavefield over the mainland and over the Atlantic immediately to the east (where the spheric contribution is similar) throw new light on the problem. It appears that the principal role of the PLHR may be to sustain duct structure and multihop propagation which is relatively much rarer over the Atlantic. Typical industrial PLHR consists of a series of narrow pulses at twice the mains frequency. It is suggested that these ‘artificial spherics’ may help to sustain the WPI and multihop duct structure. At L = 4, Yoshida et al. (1980) have shown that there is a strong, sharp maximum for WPIs originating in spherics, at f ? 3 kHz, in good agreement with Siple observations.     

Lower and Upper Ionosphere of Mars

The ionosphere of Mars has been explored mostly with the radio occultation experiment onboard Mariners 6, 7, 9; Mars 2, 3, 4, 6; Viking 1, 2, and more recently on Mars Global Surveyor (MGS) and Mars Express (MEX). In addition to the radio occultation experiment, MEX also carried Mars Advanced Radar for the Subsurface and Ionosphere Sounding (MARSIS) experiment which provided electron density profiles well above the main ionospheric peak. The atmosphere of Mars was measured directly by the neutral mass spectrometer onboard Viking 1 and 2 Landers. Later, an accelerometer and radio occultation experiment on MGS provided large data sets of atmospheric density at various locations in the upper and lower atmospheres of Mars, respectively. In this paper we review results of these upper and lower atmospheric/ionospheric measurements. Results of these measurements have been compared with theoretical models by several workers; therefore, we also review various atmospheric and ionospheric models of Mars.     

Voyager mission description

The Voyager Project, managed by the Jet Propulsion Laboratory, involves the lauching of two advanced spacecraft to explore the Jovian and Saturnian systems, as well as interplanetary space. The one-month lauch period opens on August 20, 1977, with arrivals at Jupiter in March and July of 1979, and at Saturn in November of 1980 and August of 1981. Gravity-assist swingbys of Jupiter are utilized in order to reduce the lauch energy demands needed to reach Saturn. In addition, a gravity-assist targeting option at Saturn will be maintained on the second-arriving Voyager for a possible continuation on to Uranus, with arrival in January of 1986. Flight through the Jovian and Saturnian systems will achieve close to moderate flyby encounters with several of the natural satellites, including special flyby geometry conditions for Io and Titan, as well as an Earth occultation of the spacecraft's radio signal by the rings of Saturn. The purpose of this paper is to describe the Voyager mission characteristics in order to establish a framework upon which to better understand the objectives and goals of the eleven scientific investigations which are described in subsequent papers.     

Searching for Black Holes in Space

Although General Relativity had provided the physical basis of black holes, evidence for their existence had to await the Space Era when X-ray observations first directed the attention of astronomers to the unusual binary stars Cygnus X-1 and A0620-00. Subsequently, a number of faint Ariel 5 and Uhuru X-ray sources, mainly at high Galactic latitude, were found to lie close to bright Seyfert galaxies, suggesting the nuclear activity in AGN might also be driven by accretion in the strong gravity of a black hole. Detection of rapid X-ray variability with EXOSAT later confirmed that the accreting object in an AGN is almost certainly a supermassive black hole.     

Global Muon Detector Network Used for Space Weather Applications

In this work, we summarize the development and current status of the Global Muon Detector Network (GMDN). The GMDN started in 1992 with only two muon detectors. It has consisted of four detectors since the Kuwait-city muon hodoscope detector was installed in March 2006. The present network has a total of 60 directional channels with an improved coverage of the sunward Interplanetary Magnetic Field (IMF) orientation, making it possible to continuously monitor cosmic ray precursors of geomagnetic storms. The data analysis methods developed also permit precise calculation of the three dimensional cosmic ray anisotropy on an hourly basis free from the atmospheric temperature effect and analysis of the cosmic ray precursors free from the diurnal anisotropy of the cosmic ray intensity.     

Observational Appearance of Black Holes in X-Ray Binaries and AGN

Accretion onto black holes powers most luminous compact sources in the Universe. Black holes are found with masses extending over an extraordinary broad dynamic range, from several to a few billion times the mass of the Sun. Depending on their position on the mass scale, they may manifest themselves as X-ray binaries or active galactic nuclei. X-ray binaries harbor stellar mass black holes—endpoints of the evolution of massive stars. They have been studied by X-ray astronomy since its inception in the early 60-ies, however, the enigma of the most luminous of them—ultra-luminous X-ray sources, still remains unsolved. Supermassive black holes, lurking at the centers of galaxies, are up to hundreds of millions times more massive and give rise to the wide variety of different phenomena collectively termed “Active Galactic Nuclei”. The most luminous of them reach the Eddington luminosity limit for a few billions solar masses object and are found at redshifts as high as z≥5–7. Accretion onto supermassive black holes in AGN and stellar- and (possibly) intermediate mass black holes in X-ray binaries and ultra-luminous X-ray sources in star-forming galaxies can explain most, if not all, of the observed brightness of the cosmic X-ray background radiation. Despite the vast difference in the mass scale, accretion in X-ray binaries and AGN is governed by the same physical laws, so a degree of quantitative analogy among them is expected. Indeed, all luminous black holes are successfully described by the standard Shakura-Sunyaev theory of accretion disks, while the output of low-luminosity accreting black holes in the form of mechanical and radiative power of the associated jets obeys to a unified scaling relation, termed as the “fundamental plane of black holes”. From that standpoint, in this review we discuss formation of radiation in X-ray binaries and AGN, emphasizing their main similarities and differences, and examine our current knowledge of the demographics of stellar mass and supermassive black holes.     

Evolution of the Moon: The 1974 model

The geology of the decade of Apollo and Luna probably will become one of the fundamental turning points in the history of all science. For the first time, the scientists of the Earth have been presented with the opportunity to interpret their home planet through the direct investigations of another. Mankind can be proud and take heart in this fact. The interpretive evolution of the Moon can be divided now into seven major stages beginning sometime near the end of the formation of the solar system. These stages and their approximate durations in time are as follows:

1. The Beginning — 4.6 billion years ago.
2. The Melted Shell — 4.6–4.4 billion years ago.
3. The Cratered Highlands — 4.4–4.1 billion years ago.
4. The Large Basins — 4.1–3.9 billion years ago.
5. The Light-colored Plains — 3.9–3.8 billion years ago.
6. The Basaltic Maria — 3.8–3.0 (?) billion years ago.
7. The Quiet Crust — 3.0 (?) billion years ago to the present.

The Apollo and Luna explorations that permit us to study these stages of evolution each have contributed in progressive and significant ways. Through them we now can look with new insight into the early differentiation of the Earth, the nature of the Earth's protocrust, the influence of the formation of large impact basins in that crust, the effects of early partial melting of the protomantle and possibly the earliest stages of the breakup of the protocrust into continents and ocean basins.     

The atmosphere of Venus

The investigations of Venus take a special position in planetary researches. It was just the atmosphere of Venus where first measurements in situ were carried out by means of the equipment delivered by a space probe (Venera 4, 1967). Venus appeared to be the first neighbor planet whose surface had been seen by us in the direct nearness made possible by means of the phototelevision device (Venera 9 and Venera 10, 1975). The reasons for the high interest in this planet are very simple. This planet is like the Earth by its mass, size and amount of energy obtained from the Sun and at the same time it differs sharply by the character of its atmosphere and climate. We hope that the investigations of Venus will lead us to define more precisely the idea of complex physical and physical-chemical processes which rule the evolution of planetary atmospheres. We hope to learn to forecast this evolution and maybe, in the far future, to control it. The last expeditions to Venus carried out in 1978 — American (Pioneer-Venus) and Soviet (Venera 11 and 12) — brought much news and it is interesting to sum up the results just now. The contents of this review are:

1. The planet Venus — basic astronomical data.
2. Chemical composition.
3. Temperature, pressure, density (from 0 to 100 km).
4. Clouds.
5. Thermal regime and greenhouse effect.
6. Dynamics.
7. Chemical processes.
8. Upper atmosphere.
9. Origin and evolution.
10. Problems for future studies

Here we have attempted to review the data published up to 1979 and partly in 1980. The list of references is not exhaustive. Publications of special issues of magazines and collected articles concerning separate space expeditions became traditional last time. The results obtained on the Soviet space probes Venera 9, 10 (the first publications) are collected in the special issues of Kosmicheskie issledovanija (14, Nos. 5, 6, 1975), analogous material about Venera 11, 12 is given at Pis'ma Astron. Zh. (5, Nos. 1 and 5, 1978), and in Kosmicheskie issledovanija (16, No. 5, 1979). The results of Pioneer-Venus mission are represented in two Science issues (203, No. 4382; 205, No. 4401) and special issue of J. Geophys. Res. (1980). We shall mention some articles to the same topic among previous surveys: (Moroz, 1971; Sagan, 1971; Marov, 1972; Hunten et al., 1977; Hoffman et al., 1977) and also the books by Kuzmin and Marov (1974) and Kondrat'ev (1977). Some useful information in the part of ground-based observations may be found in the older sources (for example, Sharonov, 1965; Moroz, 1967). For briefness we shall use as a rule the abbreviations of space missions names: V4 instead of Venera 4, M10 instead of Mariner 10 and so on. The first artificial satellites of Venus in the world (orbiters Venera 9 and 10) we shall mark as V9-O, V10-O unlike the descent probes V9, V10. Fly-by modules of Venera 11 and Venera 12 we shall mark as V11-F and V12-F. Pioneers descent probes — Large (Sounder), Day, Night and North — will be marked as P-L, P-D, P-Ni, P-No, orbiter as P-O, and bus as P-B.     

Measurements of quasi-static electric fields between 3 and 7 Earth radii on GEOS-1

Quasi-static electric fields have been measured with two spherical probes supported by cable booms providing a baseline of 42 m for the measurement. The performance of the experiment is outlined to demonstrate that electric fields can be measured with accuracies of ±0.7 mV m^-1 and ±1.0 mV m^-1 in the dawn-dusk and satellite-sun directions respectively. These uncertainties can be considerably reduced under favourable plasma conditions. Examples of typical observations are described.

1. The average electric field is always characterized by an irregular structure with time scales 0.5–5 min and with amplitudes of a few mV m^-1.
2. During substorms dawn-dusk electric fields up to 20–30 mV m^-1 have been observed over intervals of 30–60 s.
3. Oscillating electric fields with peak-to-peak amplitudes up to 10 mV m^-1 and periods of 3–10 min have been observed following magnetospheric disturbances.

The observations are discussed in terms of plasma motions and possible spatial scale sizes of the phenomena, standing magnetospheric wave modes and electrostatic potentials.     

Conjugate features of magnetospheric electron dynamics observed at balloon altitudes

Analysis of recent observations (from balloons, spacecraft, and surface observatories) demonstrate regional, shell, and nearpoint conjugacy at L ～ 7 during precipitative events which were characterized by local acceleration as well as release of gradient-drifted electrons injected during substorms. A number of new features of magnetospheric dynamics relating to substorm development and sudden-commencement effects, have been brought to light which, though poorly understood at present, may prove of considerable importance and are worthy of further investigation.

1. During the initial period of instability in substorm evolution, preceding the slower magnetotail convective injection, precipitation of waves of electrons in rapid polewards motion exhibit L-shell conjugacy near midnight.
2. Transient, large scale expansions of the magnetospheric electron population accompanied by temporally imbedded substorms display large scale regional conjugacy and are simultaneously observed as similarly transient intensity dropouts at balloon altitudes.
3. Precipitation from gradient-drifting electrons in the dayside magnetosphere exhibits near point-conjugacy, at least down to the order of 50 km and quite probably less.

Similarly tight conjugacy applies to the release of electrons showing a specific local response to sudden commencements.

1. Analysis of the approach to and attainment of spectral equilibrium in the precipitation observed from drifting electrons may provide information about either, or both, the source spectrum at injection and the process of local release.
2. The specific precipitation effect sometimes observed at the time of an SC remains a rather puzzling feature, although it seems clear now that the acceleration and/or release process responsible is of a highly local nature and works selectively at small pitch angles well within the magnetospheric boundary. Coupling of the interplanetary shock with the magnetosphere must be an important aspect, but the details are not clear as yet.
3. On at least one occasion, a large part (perhaps all) of the magnetospheric electron population varied in a nearly synchronous manner in response to solar wind induced distortions during the variable compressive phase of a sudden commencement geomagnetic storm.

In the ongoing effort to identify and understand acceleration and release mechanisms involved in magnetospheric dynamics, balloon-borne experiments will continue to be useful, providing essential information presently unattainable by other means.     

A New Morphology of Solar Activity and Recurrent Geomagnetic Disturbances: The Late-Declining Phase of the Sunspot Cycle

Certain aspects of the Sun and resulting geomagnetic disturbances can be studied better on the source surface, an imaginary spherical surface of 3.5 solar radii, than on the photospheric surface. This paper presents evidence that the Sun exhibits one of the most fundamental aspects of activities most clearly during the late-declining phase of the sunspot cycle. It is the period when 27-day average values of the solar wind speed and of geomagnetic disturbances tend to be highest during the sunspot cycle. Important findings of this study on the late-declining phase of the sunspot cycle are the following:

1. By introducing a new coordinate system, modifying the Carrington coordinates, it is shown that various solar activity phenomena, solar flares, the brightest coronal regions, and also the lowest solar wind speed region, tend to concentrate in two quadrants, one around 90^° in longitude in the northern hemisphere (NE) and the other around 270^° in longitude in the southern hemisphere (SW). For this reason, the new coordinate system is referred to as the NESW coordinate system.
2. It is shown that the above results are closely related to the fact that the neutral line exhibits a single wave (sinusoidal or rectangular) in both the Carrington coordinates and the NESW coordinate system during the late-declining phase. The shift of the neutral line configuration during successive solar rotations during the late-declining phase causes longitudinal scatter of the location of solar flares with respect to the neutral line in a statistical study. The NESW coordinate system is designed to suppress the shift, so that the single wave location is fixed and thus a ‘nest’ of solar flares emerges in the NE and SW quadrants.
3. It is also shown that the single wave is the source of the double peak of the solar wind speed and two series of recurrent geomagnetic disturbances in each solar rotation, making the 27-day average solar wind and geomagnetic disturbances highest during the sunspot cycle. The double peak is a basic feature during the late-declining phase, but is obscured by several complexities which we identified in this paper; see item 8.
4. The single wave of the neutral line configuration can be approximated by three dipole fields, one which can be represented by a central dipole (parallel or anti-parallel to the rotation axis) and two hypothetical dipoles on the photosphere. This configuration is referred to as the triple dipole model.
5. The location of the two hypothetical photospheric dipoles coincide with the two active regions (solar flares, the brightest coronal region) and also the lowest solar wind speed region in the NESW coordinate system; the lowest solar wind regions are the cause of the valleys of the double peak of the solar wind speed.
6. The two hypothetical dipole fields actually do exist at the location of the two active regions in a coarse magnetic map (5 × 5^°). The two dipoles follow the Hale–Nicholson polarity law. Thus, they are real physical entities.
7. The apparent meridional rotation of the dipolar field on the source surface during the sunspot cycle results from combined changes of both the central dipole field and of the two photospheric dipoles, although the central dipole remains axially parallel or anti-parallel. Thus, the Sun has a general field that can be represented by an axially aligned dipole located at the center of the Sun throughout the sunspot cycle, except for the sunspot maximum period when the polarization reversal occurs.
8. The complexity of recurrent geomagnetic disturbances can also be understood by having the NESW coordinate system for various solar phenomena and the relative location of the earth with respect to the solar equatorial plane.
9. As the intensity of the two dipoles decreases toward the end of the sunspot cycle, the amplitude of the single wave decreases, and the neutral line tends to align with the heliographic equator.
10. The neutral line shows a double wave structure during certain epochs of the sunspot cycle. In such a situation, it can be considered that two NESW coordinate systems are present in one Carrington coordinate, resulting in four active regions.
11. The so-called classical “sector boundary” arises when the peaks (top and bottom) of the single wave reached 90^° in latitude in both hemispheres.
12. In summary: A study of the late-declining period of the sunspot cycle is very important compared with the sunspot maximum period. In the late-declining period, the Sun shows its activities in the simplest form. It is suggested that some of the basic features of solar activities and recurrent geomagnetic disturbances that have been studied by many researchers in the past can be synthesized in a simplest way by introducing the NESW coordinate system and the triple dipole model. There is a possibility that the basic results we learned during the late phase of the sunspot cycle can be applicable to the rest of the sunspot cycle.

     

Particle code simulations with injected particles

As problems we are interested in become more complex, we often find our simulations stretching the limits of available computer resources. For example, an interesting problem is simulation of dissipation processes in sub-critical collisionless shocks. To simulate this system our simulation box must contain the shock and its upstream and downstream regions over the entire length of a run. If the shock moves with any appreciable speed the box must then be considerably larger than the shock thickness making it hard to resolve the shock front itself with a reasonable number of grid points. A solution to this problem is to run the simulation in the frame of reference of the shock. Particles are injected upstream of the shock and leave the simulation box downstream. With the shock stationary in the simulation box, we only need to contain enough of the up and downstream regions for the fields, etc., to settle down and separate the shock from the box boundaries. In this tutorial we consider some basic algorithms used in a practical particle injection code, such as the two dimensional WAVE code used at Los Alamos. We will try to present these ideas in a simple format general enough to be easily included in any particle code. Topics covered are:

1. Smoothly Injecting Particles.
2. Generating the Distribution Functions.
3. Time Dependent Injection Density.
4. Boundary Conditions on Fields and Particles.

(Flux and Charge Conservation)     

Hard rock cosmic ray archaeology

Recent examinations of extraterrestrial materials exposed to cosmic rays for different intervals of time during the geological history of the solar system have generated a wealth of new information on the history of cosmic radiation. This information relates to the temporal variations in

1. the flux and energy spectrum of low energy (solar) protons of ? 10 MeV kinetic energy;
2. the flux and energy spectrum of (solar) heavy nuclei of Z > 20 of kinetic energy, 0.5–10 MeV/n;
3. the integrated flux of protons and heavier nuclei of ? 0.5 GeV kinetic energy, and
4. the flux and energy spectrum of nuclei of Z > 20 of medium energy — 100–2000 MeV/n kinetic energy.

The above studies are entirely based on the natural detector method which utilises two principal cosmogenic effects observed in rocks, (i) isotopic changes and (ii) changes in the crystalline structure of rock constituents, due to cosmogenic interactions. The information available to date in the field of hard rock cosmic ray archaeology refers to meteorites and lunar rocks/soil. Additional information based on study of cosmogenic effects in man-made materials exposed to cosmic radiation in space is also discussed. It is shown that the natural detectors inspite of their extreme simplicity have begun to provide cosmic ray information in a very quantitative and precise manner comparable to the most sophisticated electronic particle detectors. The single handicap in using the hard rock detectors is however the uncertainty regarding their manner of exposure, geometry etc. At present, a variety of techniques are being used to study the evolutionary history of extraterrestrial materials and as this field grows, uncertainties in cosmic ray archaeology will correspondingly decrease.     

Evolution into contact of the low-mass close binary systems

We investigated the effect of mass accretion on the secondary components in close binomy systems (M _total ≤ 2.5 M _⊙ M _2,0 ≤ 0.75 M _⊙) exchanging mass in the case A. The evolution of the low-mass close binary systems (M _total ≤ 2.5 M _⊙) exchanging the mass in the case A depends on the three main factors:

-the initial mass ratio (q ₀ = M _2,0/M _1,0), which determines the rate of mass transfer between components;

-the inital mass of the secondary component (M _2,0) and

-the effectiveness of the heating of the photosphere of the secondary component, by infalling matter.

The second factor allows to divide all systems into two essentially different groups:

1. systems in which the secondary component is a star with a radiative envelope, or with a thin convection zone in the uppermost layers;
2. and systems in which secondary component has a thick convective envelope or is fully convective.

The systems from the first group evolve into contact in a characteristic time scale 10⁵ – 10⁷ years, and reach contact after transfering of 0.03 – 0.3 M _⊙. The mass exchange proceeds only in a thermal time scale. For the systems from the group b the effectiveness of the heating of the stellar surface is the most important. In the case when the entropy of the newly accreted matter is the same as the surface entropy of the secondary, a convective star should shrink upon accretion. Then contact binaries are not formed. In the case when the entropy of the infalling matter is greater then that on the surface, the reaction of the secondary is different. The radius of the secondary component grows rapidly in response to accretion, and the systems reaches contact after the 10³ – 3 10⁶ years, and after transfer of 0.002 – 0.2. M _⊙. The reaction of the secondary is determined by the formation of the temperature inversion layer below the stellar surface. Full references in: Sarna, M.J. and Fedorova, A.V. (1988) “Evolutionary status of W UMa-type Binaries — Evolution into contact”, Astron. Astrophys., in press.