ALBERT

Description: We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus. The interaction of Io with the plasma surrounding it has been a subject of interest for almost 30 years, dating from the discovery by Bigg (1964) that radio emissions from the Jovian magnetosphere are controlled by Io's position. Since that time, both ground-based and spacecraft observations have shown that Io is a unique satellite that influences the Jovian magnetosphere in important ways. In particular, material from Io is a major source of plasma for the magnetosphere, and the energy that this plasma harnesses from Jupiter's co-rotating magnetic field is an important power source for the magnetosphere. It is apparent that the local interaction of the torus plasma with Io plays a key role in the formation, composition, and energetics of the Io torus; the interaction is also highly nonlinear. We have modeled this interaction using time-dependent three-dimensional magnetohydrodynamic (MHD) simulations. During this past year, we have used NASA support to develop a new MHD code to study the interaction. As part of the Galileo spacecraft's recent successful insertion into orbit around Jupiter, the spacecraft passed within 900 km of Io's surface. Our calculations have focused on using Galileo particles and fields data to examine a question that was not resolved by the Voyager observations: Does Io have an intrinsic magnetic field? In this progress summary, we describe our efforts on this problem to date.

Description: The eruption of Mt. Pinatubo in June 1991 produced the largest enhancement of stratospheric aerosol loading ever observed by lidar over Hampton, Virginia. Low altitude layers (less than 20 km) were the first to arrive over Hampton in early August, the result of transport associated with a tropospheric anticyclonic cell over North America. The maximum peak scattering ratio, 34 at 22.4 km, and the maximum stratospheric integrated backscatter of 0.0053 sr(exp -1), both at 694 nm, observed since the eruption were measured on February 20, 1992. After decreasing during the spring and summer of 1992, the aerosol burden increased significantly during the winter of 1992-93, evidence of a poleward winter transport from the equatorial reservoir. Over the period from February 1992 to February 1994, the stratospheric aerosol loading decreased with an average 1/e decay time of 10.1 months. The vertical distribution, intensity, and transport of Pinatubo aerosols over this site are described and compared with similar measurements after El Chichon.

Description: T-matrix computations of light scattering by polydispersions of randomly oriented nonspherical aerosols and Mie computations for equivalent spheres are compared. Findings show that even moderate nonsphericity results in suubstantial errors in the retrieved aerosol optical thickness if satellite reflectance measurements are analyzed using Mie theory. On the other hand, the use of Mie theory for nonspherical aerosols produces negligible errors in the computation of albedo and flux related quantities, provided that the aerosol size distribution and optical thickness are known beforehand. The first result can be explained by large nonspherical-spherical differences in scattering phase function, while the second result follows from small nonspherical-spherical differences in single-scattering albedo and asymmetry parameter. No cancellation of errors occurs if one consistently uses Mie theory in the retrieval algorithm and then in computing the albedo for the retrieved aerosol optical thickness.

Description: Azimuthal asymmetries in the atmospheric refractive index can lead to errors in estimated vertical and horizontal station coordinates. Daily average gradient effects can be as large as 50 mm of delay at a 7 deg elevation. To model gradients, the constrained estimation of gradient paramters was added to the standard VLBI solution procedure. Here the analysis of two sets of data is summarized: the set of all geodetic VLBI experiments from 1990-1993 and a series of 12 state-of-the-art R&D experiments run on consecutive days in January 1994. In both cases, when the gradient parameters are estimated, the overall fit of the geodetic solution is improved at greater than the 99% confidence level. Repeatabilities of baseline lengths ranging up to 11,000 km are improved by 1 to 8 mm in a root-sum-square sense. This varies from about 20% to 40% of the total baseline length scatter without gradient modeling for the 1990-1993 series and 40% to 50% for the January series. Gradients estimated independently for each day as a piecewise linear function are mostly continuous from day to day within their formal uncertainties.

Description: Particles leving the neutral sheet in the distant magnetotail at times display adiabatic trajectory sequences characterized by an inflection toward the equator and subsequent mirroring in its vicinity. We demonstrate that this low-latitude mirroring results primarily from a centrifugal deceleration due to the fast direction-changing E x B drift. This effect which we refer to as 'centrifugal trapping' appears both in guiding centre and full particle treatments. It thus does not directly relate to nonadiabatic motion. However, pitch angle scattering due to nonadiabatic neutral sheet interaction does play a role in reducing the parallel speed of the particles. We show that centrifugal trapping is an important mechanism for the confinement of the slowest (typically below the equatorial E x B drift speed) plasma sheet populations to the midplane vicinity.

Description: Estimates of the mass of dust suspended in the Martian atmosphere are derived from global and regional 9-micrometer opacity maps produced from Viking Infrared Thermal Mapper data. During the peak of the 1977b storm, a total dust mass of approximately 4.3 x 10(exp 14) g was suspended, equivalent to 4.3 x 10(exp -4) g/sq cm, or a layer 1.4 micrometers thick. During a local dust storm near Solis Planum at L(sub s) 227 deg, approximately 1.3 x 10(exp 13) g of dust were lofted, equal to about a 6-micrometer layer in that vicinity.

Description: We model an infrared outburst on Io as being due to a large, erupting lava flow which increased its area at a rate of 1.5 x 10(exp 5)/sq m and cooled from 1225 to 555 K over the 2.583-hr period of observation. The inferred effusion rate of 3 x 10(exp 5) cu m/sec for this eruption is very high, but is not unprece- dented on the Earth and is similar to the high eruption rates suggested for early lunar volcanism. Eruptions occur approxi- mately 6% of the time on Io. These eruptions provide ample resurfacing to explain Io's lack of impact craters. We suggest that the large total radiometric heat flow, 10(exp 14) W, and the size and temperature distribution of the thermal anomalies (McEwen et al. 1992; Veeder et al. 1994) can be accounted for by a series of silicate lava flows in various stages of cooling. We propose that the whole suite of Io's currently observed thermal anomalies was produced by multiple, high-eruptive-rate silicate flows within the past century.

Description: Measurements of superthermal electron fluxes in the solar wind indicate that field lines within coronal mass ejections, CMEs, near and beyond 1 AU are normally connected to the Sun at both ends. However, on occasion some field lines embedded deep within CMEs appear to be connected to the Sun at only one end. Here we propose an explanation for how such field lines arise in terms of 3-dimensional reconnection close to the Sun. Such reconnection also provides a natural explanation for the flux rope topology characteristic of many CMEs as well as the coronal loops formed during long-duration, solar X-ray events. Our consideration of the field topologies resulting from 3-dimensional reconnection indicates that field lines within and near CMEs may on occasion be connected to the outer heliosphere at both ends.

Description: We investigate the orbital dynamics of small dust particles generated via the continuous micrometeoroid bombardment of the Martian moons. In addition to Mar's oblateness, we also consider the radiation pressure perturbation that is complicated by the planet's eccentric orbit and tilted rotational axis. Considering the production rates and the lifetimes of dust grains, we show that particles from Deimos with radii of about 15 micrometers are expected to dominate the population of a permanently present and tilted dust torus. This torus has an estimated peak number density of approximately equals 5 x 10(exp -12)/cu cm and an optical depth of approximately equals 4 x 10(exp -8).

Description: We examined 11 cases when the interplanetary magnetic field (IMF) was intensely northward (greater than 10 nT) for long durations of time (greater than 3 hours), to quantitatively determine an uppler limit on the efficiency of solar wind energy injection into the magnetosphere. We have specifically selected these large B(sub N) events to minimize the effects of magnetic reconnection. Many of these cases occurred during intervals of high-speed streams associated with coronal mass ejections when viscous interaction effects might be at a maximum. It is found that the typical efficiency of solar wind energy injection into the magnetosphere is 1.0 x 10(exp -3) to 4.0 x 10(exp -3), 100 to 30 times less efficient than during periods of intense southward IMFs. Other energy sinks not included in these numbers are discussed. Estimates of their magnitudes are provided.

Description: The valley network channels on the heavily cratered ancient surface of Mars suggest the presence of liquid water approximately 3.8 Gyr ago. However, the implied warm climate is difficult to explain in the context of the standard solar model, even allowing for the maximum CO2 greenhouse heating. In this paper we investigate the astronomical and planetary implications of a nonstandard solar model in which the zero-age, main-sequence Sun had a mass of 1.05 +/- 0.02 Solar Mass. The excess mass was subsequently lost in a solar wind during the first 1.2(-0.2, +0.4) Gyr of the Sun's main sequence phase. The implied mass-loss rate of 4(+3, -2) x 10(exp -11) M/yr, or about 10(exp 3) x that of the current Sun, may be detectable in several nearby young solar type stars.

Description: Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Description: Visible and near-IR refectivity, Moessbauer, and X ray diffraction data were obtained on powders of impact melt rock from the Manicouagan Impact Crater located in Quebec, Canada. The iron mineralogy is dominated by pyroxene for the least oxidized samples and by hematite for the most oxidized samples. Phyllosilicate (smectite) contents up to approximately 15 wt % were found in some heavily oxidized samples. Nanophase hematite and/or paramagnetic ferric iron is observed in all samples. No hydrous ferric oxides (e.g., goethite, lepidocrocite, and ferrihydrite) were detected, which implies the alteration occurred above 250 C. Oxidative alteration is thought to have occurred predominantly during late-stage crystallization and subsolidus cooling of the impact melt by invasion of oxidizing vapors and/or solutions while the impact melt rocks were still hot. The near-IR band minimum correlated with the extent of aleration Fe(3+)/Fe(sub tot) and ranged from approximately 1000 nm (high-Ca pyroxene) to approximately 850 nm (bulk, well-crystalline hematite) for least and most oxidized samples, respectively. Intermediate band positions (900-920 nm) are attributed to low-Ca pyroxene and/or a composite band from hematite-pyroxene assemblages. Manicouagan data are consistent with previous assignments of hematite and pyroxene to the approximately 850 and approximately 1000nm bands observed in Martian reflectivity spectra. Manicouagan data also show that possible assignments for intermediate band positions (900-920 nm) in Martian spectra are pyroxene and/or hematite-pyroxene assemblages. By analogy with impact melt sheets and in agreement with observables for Mars, oxidative alteration of Martian impact melt sheets above 250 C and subsequent erosion could produce rocks and soils with variable proportions of hematite (both bulk and nanophase), pyroxene, and phyllosilicates as iron-bearing mineralogies. If this process is dominant, these phases on Mars were formed rapidly at relativly high temperatures on a sporadic basis throughout the history of the planet. The Manicouagan samples also show that this mineralogical diversity can be accomplished at constant chemical composition, which is also indicated for Mars from the analyses of soil at the two Viking landing sites.

Description: The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Description: The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Description: The second pulsating aurora rocket (PULSAUR 2) rocket was launched into a pulsating aurora on 9 February 1994 and carried a broad range of instruments in order to perform a study of this type of aurora. The rocket measurements were complemented with a set of ground-based measurements. The particle measurements performed in the rocket are related to the ground-based optical measurements performed along the rocket trajectory. It is found that the high energy electrons are largely in phase with the measured luminosity. The EISCAT measurements carried out during the flight are reviewed. The PULSAUR 2 campaign is described. Results concerning the auroral conditions and the particle measurements are presented.

Description: The high-speed correction factor to the O(+)-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T(sub r) designate the average of the ion and neutral temperatures in K, and u = nu(sub d)/alpha, where nu(sub d) is the relative drift velocity in cm/s, and alpha = 4.56 x 10(exp 3) square root of T(sub r) cm/s is the thermal velocity of the O(+)-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.

Description: This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Description: The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Description: Empirical data-based models of the magnetosphereic magnetic field have been widely used during recent years. However, the existing models (Tsyganenko, 1987, 1989a) have three serious deficiencies: (1) an unstable de facto magnetopause, (2) a crude parametrization by the K(sub p) index, and (3) inaccuracies in the equatorial magnetotail B(sub z) values. This paper describes a new approach to the problem; the essential new features are (1) a realistic shape and size of the magnetopause, based on fits to a large number of observed crossing (allowing a parametrization by the solar wind pressure), (2) fully controlled shielding of the magnetic field produced by all magnetospheric current systems, (3) new flexible representations for the tail and ring currents, and (4) a new directional criterion for fitting the model field to spacecraft data, providing improved accuracy for field line mapping. Results are presented from initial efforts to create models assembled from these modules and calibrated against spacecraft data sets.

Description: We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO(x) chemistry on the ozone response above 45 km. The model revealed the general anti-correlation of ozone and solar UV at 65-75 km, and simulated strong UV responses of water vapor and HO(x) species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over +/- 40 was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.

Description: The accuracy of Geosat satellite altimetry over the Greenland ice sheet is evaluated by comparing the measured heights to radar elevations from the airborne Greenland Aerogeophysics Project. At the center of the ice sheet where the ice surface is nearly level, surface comparisons show a fit at the 1 to 3 m level as expected, but even at moderately sloping ice regions (0.3 deg-0.6 deg), satellite altimetry mean errors in the range of 10 to 35 m are observed. These errors are found for slope-corrected and waveform-retracked data, so most previous accuracy estimates of current satellite altimetry ice sheet elevations in regions of slopping or undulating ice appear to be too optimistic.

Description: We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.

Description: A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Description: Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Description: Many problems in geophysical and astrophysical convection systems are characterized by fast rotation and spherical shell geometry. The combined effects of Coriolis forces and spherical shell geometry produce a unique spatial symmetry for the convection pattern in a rapidly rotating spherical shell. In this paper, we first discuss the general spatial symmetries for rotating spherical shell convection. A special model, a spherical shell heated from below, is then used to illustrate how and when the spatial symmetries are broken. Symmetry breaking occurs via a sequence of spatial transitions from the primary conducting state to the complex multiple-layered columnar structure. It is argued that, because of the dominant effects of rotation, the sequence of spatial transitions identified from this particular model is likely to be generally valid. Applications of the spatial symmetry breaking to planetary convection problems are also discussed.

Description: Initial observations of a newly documented type of optical emission above thunderstorms are reported. 'Blue jets', or narrowly collimated beams of blue light that appear to propagate upwards from the tops of thunderstorms, were recorded on B/W and color video cameras for the first time during the Sprites94 aircraft campaign, June-July, 1994. The jets appear to propagate upward at speeds of about 100 km/s and reach terminal altitudes of 40-50 m. Fifty six examples were recorded during a 22 minute interval during a storm over Arkansas. We examine some possible mechanisms, but have no satisfactory theory of this phenomenon.

Description: It is shown from flux transfer event (FTE) occurrence statistics, observed as a function of MLT by the ISEE satellites, that recent 2-dimensional analytic theories of the effects of pulsed Petschek reconnection predict FTEs to contribute between 50 and 200 kV to the total reconnection voltage when the magnetosheath field points southward. The upper limit (200 kV) allows the possibility that FTEs provide all the antisunward transport of open field lines into the tail lobe. This range is compared with the voltages associated with series of FTEs signatures, as inferred from ground-based observations, which are in the range 10-60 kV. We conclude that the contribution could sometimes be made by a series of single, large events; however, the voltage is often likely to be contributed by several FTEs at different MLT.

Description: New investigations of the core ion motion within high-latitude topside ionosphere and near-Earth magnetosphere, using data from the Dynamics Explorer (DE) retarding ion mass spectrometer (RIMS), reveal the existence of significant regions of downward moving O(+). The occurences of downgoing versus upgoing O(+) are not clearly separable in terms of either polar zenith angle or Kp but are well distinquished by the direction of the z component of the interplanetary magnetic field (IMF). On the average, down flow dominates when IMFB2 less than O, while upflow dominates for IMFB2 greater than 0. Combining cross-field convection velocities derived from an empirical convection electric field model with the observed parallel velocities yields a two dimensional ion velocity field. This velocity field is consistent with a senario which has O(+) of cusp/cleft and auroral zone origin concvecting into the polar cap and, because of the dominance of gravitational energy over the upward kinetic energy, falling back into the inonsphere. This provides additional confirmation of the results of studies of the cleft ion fountain. Estimates of the flux of O(+) in the upflowing and downflowing regions for Lambda greater than or equal to 60 deg give a total upflow of approximately 6 X 10(exp 25) ions/sec for IMFB2 greater than 0 and total upflow and downflow of approximately 4 X 10 (exp 25) ions/sec and 1 x 10(exp 25) ions/sec, respectively, for IMFB less than 0. In all cases the magnitude of the dayside outflow is consistent with previous work on upwelling ions. While the magnitudes vary for high and low Kp the ratios of upward to downward flow are roughly the same at approximately 1.7. The downflowing O(+) shows a correlation with the magnitude of the outflow of light ions in the same region but the cause and effect of this relationship is not distinquishable. hable.

Description: We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Description: The effects on double-probe electric field measurements induced by electron density and temperature gradients are investigated. We show that on some occasions such gradients may lead to marked spurious electric fields if the probes are assumed to lie at the same probe potential with repect to the plasma. The use of a proper bias current will decrease the magnitude of such an error. When the probes are near the plasma potential, the magnitude of these error signals, delta Epsilon, can vary as delta Epsilon approx. T(sub e)(Delta n(sub e)/n(sub e)) + 0.5 Delta T(sub e), where T(sub e) is the electron temperature, Delta n(sub e)/ n(sub e) the relative electron density variation between the two sensors, and Delta T(sub e) the electron temperature difference between the two sensors. This not only implies that the error signals will increase linearly with the density variations but also that such signatures grow with Delta T(sub e) i.e., such effects are 10 times larger in a 10-eV plasma than in a 1-eV plasma. This type of error is independent of the probe separation distance provided the gradient scale length is much larger than the distance. The largest errors occur when the probes are near to the plasma potential. During the crossing of a small structure (e.g, a double layer) the error signal appears as a bipolar signature. Our analysis shows that errors in double-probe measurements caused by plasma gradients are not significant at large scale (much greater than 1 km) plasma boundaries, and may only be important in cases where small-scale (less than 1 km), internal gradient structures exist. Bias currents tailored for each plasma parameter regime (i.e., variable bias current) would improve the double-probe response to gradient effects considerably.

Description: The evolution of the volcanic debris plume originating from the June 1991 eruption of Mt. Pinatubo has been monitored since its genesis using a ground-based backscatter lidar facility sited at the Jet Propulsion Laboratory (JPL). Both absolute and relative pre- and post-Pinatubo backscatter observations are in accord with Mie scattering projections based on measured aerosol particle size distributions reported in the literature. The post-Pinatubo column-integrated backscatter coefficient peaked approximately 400 days after the eruption, and the observed upper boundary of the aerosol column subsided at a rate of approximately 200 m/mon.

Description: Using Global Positioning System (GPS) receivers, we reoccupied several leveling benchmarks on the Kenai Peninsula of Alaska which had been surveyed by conventional leveling immediately following the March 27, 1964, Prince William Sound earthquake (M(sub w) = 9.3). By combining the two sets of measurements with a new, high-resolution model of the geoid in the region, we were able to determine the cumulative 1993-1964 postseismic vertical displacement. We find uplift at all of our benchmarks, relative to Seward, Alaska, a point that is stable according to tide gauge data. The maximum uplift of about 1 m occurs near the middle of the peninsula. The region of maximum uplift appears to be shifted northwest relative to the point of maximum coseismic subsidence. If we use tide gauge data at Nikishka and Seward to constrain the vertical motion, then the observed uplift has a trenchward tilt (down to the southeast) as well as an arching component. To explain the observations, we use creep-at-depth models. Most acceptable models require a fault slip of about 2.75 m, although this result is not unique. If the slip has been continuous since the 1964 earthquake, then the average slip rate is nearly 100 mm/yr, twice the plate convergence rate. Comparing the net uplift achieved in 29 years with that observed over 11 years in an adjacent region southeast of Anchorage, Alaska, we conclude that the rate of uplift is decreasing. A further decrease in the uplift rate is expected as the 29-year averaged displacement rate is about twice the plate convergence rate and therefore cannot be sustained over the entire earthquake cycle.

Description: Narrow enhancements of electron precipitation, with energy and flux well above typical values, have been observed with Dynamics Explorer 2 (DE 2) in the cusp/cleft region. The electron flux in the energy range 0.2-1 keV was 2 orders of magnitude higher in these structures than in the magnetosheath and were seen in approximaetly 80% of DE 2 cusp crossings at ionospheric altitudes. Typically, there was more than one electron structure in each cusp crossing. The position of these structures showed a systematic variation: for poleward ion dispersion (energy decreases with increasing latitude), electron structures were seen more often on the equatorial boundary of the cusp, while for equatorward ion dispersion (energy decreases with decreasing latitude), electron structures were more often seen on the poleward boundary. This suggests that the electron structures are associated with newly reconnected field lines. The electron spectra suggest that field-aligned acceleration processes could produce the electron structures, first near the boundary of the cusp/cleft during the reconnection of field lines and then in the cusp/cleft during the motion of reconnected flux tubes through the polar ionosphere.

Description: A permanent Global Positioning System (GPS) receiver at Casa Diablo Hot Springs, Long Valley Caldera, California was installed in January, 1993, and has operated almost continuously since then. The data have been transmitted daily to the Jet Propulsion Laboratory (JPL) for routine analysis with data from the Fiducial Laboratories for an International Natural sciences Network (FLINN) by the JPL FLINN analysis center. Results from these analyses have been used to interpret the on going deformation at Long Valley, with data excluded from periods when the antenna was covered under 2.5 meters of snow and from some periods when Anti Spoofing was enforced on the GPS signal. The remaining time series suggests that uplift of the resurgent dome of Long Valley Caldera during 1993 has been 2.5 +/- 1.1 cm/yr and horizontal motion has been 3.0 +/- 0.7 cm/yr at S53W in a no-net-rotation global reference frame, or 1.5 +/- 0.7 cm/yr at S14W relative to the Sierra Nevada block. These rates are consistent with uplift predicted from frequent horizontal strain measurements. Spectral analysis of the observations suggests that tidal forcing of the magma chamber is not a source of the variability in the 3 dimensional station location. These results suggest that remotely operated, continuously recording GPS receivers could prove to be a reliable tool for volcanic monitoring throughout the world.

Description: A comparison is made of the vertical distribution of high-level cloud tops derived from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements and from the International Satellite Cloud Climatology Project (ISCCP) for all Julys and Januarys in 1985 to 1990. The results suggest that ISCCP overestimates the pressure of high-level clouds by up to 50-150 mbar, particularly at low latitudes. This is caused by the frequent presence of clouds with diffuse tops (greater than 50% time when cloudy events are observed). The averaged vertical extent of the diffuse top is about 1.5 km. At midlatitudes where the SAGE II and ISCCP cloud top pressure agree best, clouds with distinct tops reach a maximum relative proportion of the total level cloud amount (about 30-40%), and diffuse-topped clouds are reduced to their minimum (30-40%). The ISCCP-defined cloud top pressure should be regarded not as the material physical height of the clouds but as the level which emits the same infrared radiance as observed. SAGE II and ISCCP cloud top pressures agree for clouds with distinct tops. There is also an indication that the cloud top pressures of optically thin clouds not overlying thicker clouds are poorly estimated by ISCCP at middle latitudes. The average vertical extent of these thin clouds is about 2.5 km.

Description: Global high-level clouds identified in Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements for January and July in the period 1985 to 1990 are compared with near-nadir-looking observations from the International Satellite Cloud Climatology Project (ISCCP). Global and zonal mean high-level cloud amounts from the two data sets agree very well, if clouds with layer extinction coefficients of less than 0.008/km at 1.02 micrometers wavelength are removed from the SAGE II results and all detected clouds are interpreted to have an average horizontal size of about 75 km along the 200 km transimission path length of the SAGE II observations. The SAGE II results are much more sensitive to variations of assumed cloud size than to variations of detection threshold. The geographical distribution of cloud fractions shows good agreement, but systematic regional differences also indicate that the average cloud size varies somewhat among different climate regimes. The more sensitive SAGE II results show that about one third of all high-level clouds are missed by ISCCP but that these clouds have very low optical thicknesses (less than 0.1 at 0.6 micrometers wavelength). SAGE II sampling error in monthly zonal cloud fraction is shown to produce no bias, to be less than the intraseasonal natural variability, but to be comparable with the natural variability at longer time scales.

Description: Enhancements in the fluxes of relativistic electrons trapped within the Earth's magnetosphere have been measured by the high-energy particle spectrometer, part of the particle environment monitor on the upper atmosphere research satellite (UARS). The largest increase in the electron fluxes with energies greater than 1 MeV observed on UARS from October 1991 through July 1994 was in early May 1992. The fluxes of trapped electrons in the drift loss cone and locally precipitating electrons showed differing buildup and decay rates as a function of invariant latitude. Increases of more than 2 orders of magnitude were observed in drift loss cone fluxes at magnetic latitudes of 40 deg-66 deg and in precipitating fluxes from 48 deg to 66 deg. The energy flux contained in the most intense local precipitation observed was approximately 0.1 erg/sq cm/s, entering the atmosphere and creating up to 1000 ion pairs/cu cm/s at 55-km altitude. The daily averaged energy flux from directly precipitating electrons with energies greater than 1 MeV deposited greater than 10(exp 20) erg/d worldwide into the atmosphere for the period May 12-21, 1992, producing greater than 10(exp 31) odd nitrogen molecules below 60-km altitude.

Description: A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Description: Radiation resulting from interaction between the effluent cloud of a space shuttle thruster and the ambient atmosphere was observed with a spectograph aboard the shutttle. The spectral measurements were made between 400 and 800 nm with a resolutoion of 3 nm. The primary emissions are identified as NO2, HNO, O(1)D, and O(1)S. These are the first observations od O(1)S emission in the shuttle plume. These data are compared with the previous measurements, and possible excitation mechanisms are discussed. The results are also compared with a Monte Carlo simulation of thruster plume-atmosphere interaction radiation.

Description: We report multi-instrument observations during an isolated substorm on 17 October 1989. The European Incoherent Scatter (EISCAT) radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71 deg Lambda - 78 deg Lambda. Sub-Auroral Magnetometer Network (SAMNET) and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. Interplanetary Monitoring Platform-8 (IMP-8) magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux in the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71 deg Lambda by the time of the expansion phase onset. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the distant neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase.

Description: The impact of the extraterrestrial object that formed the Chicxulub crater in the northwestern Yucatan peninsula of Mexico is the leading suspect for the extinction of the dinosaurs. This article reports on a Planetary Society expedition to Albion Island in the Rio Hondo region of Belize to investigate evidence supporting the impact theory.

Description: Direct current (DC) electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 Universal Time (UT) on 25 October 1988, can be explained if a downward neutral wind of 15-30 m/s exists. Simultaneously, the F-region plasma was moving downward at a speed of 30-50 m/s. These events appear in the local time sector of 23.00-23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25-50 m/s are possible at time near midnight in the equatorial F-region.

Description: An overview of the observations of backstreaming electrons in the foreshock and the mechanisms that have been proposed to explain their properties will be presented. A primary characteristic of observed foreshock electrons is that their velocity distributions are spatially structured in a systematic way depending on distance from the magnetic field line which is tangent to the shock. There are two interrelated aspects to explaining the structure of velocity distributions in the foreshock, one involving the acceleration mechanism and the other, propagation from the source to the observing point. First, the source distribution of electrons energized by the shock must be determined along the shock surface. Proposed acceleration mechanisms include magnetic mirroring of incoming solar wind particles and mechanisms involving transmission of particles through the shock. Secondly, the kinematics of observable electrons streaming away from a curved shock with an initial parallel velocity and a downstream perpendicular velocity component due to the motional electric field must be determined. This is the context in which the observations and their explanations will be reviewed.

Description: Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Description: Goethite-bearing samples with values of Mn(s) (Mn/(Mn+Fe) mole fraction) up to 0.206 were synthesized by precipitation from alkaline solution. Samples with Mn(s) less than or equal 0.061 were single-phase Mn-goethites: samples with higher Mn(s) values contained another Mn-bearing phase (probably jacobsite). Mn-hematites were prepared by dehydroxylation of corresponding Mn-goethites at 500 C. Orthorhombic a and b unit cell dimensions of Mn-goethites changed in a linear manner with Mn(s), but not at rates predicted by the Vegrad law. Hexagonal unit cell dimensions of Mn-hematites did not vary with Mn(s). Moessbauer parameters isomer shift (IS), quadrupole splitting (QS), and hyperfine field (B(sub hf)) were measured at 293 and 15 K. For all single-phase Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061), magnetic splitting was observed at both temperatures. At 293 K, small but systematic decreases in B(sub hf) were observed with increasing Mn substitution; IS and QS were not dependent on Mn(s). Mn substitution strongly lowered the Morin transition temperature of hematite. At 15 K, the Morin transition was not present for Mn(s) greater than 0.020(4). The saturation magnetization of Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061) was the expected zero (within error) for antiferromagnetic goethite and for hematites obtained from dehydroxylation of goethites. Mn-geothites with Mn(s) greater than 0.061 were magnetic because of the presence of strongly magnetic jacobsite. For reflectivity spectra, bands resulting from MN(3+) were centered near 454 and 596 nm for Mn-goethites and near 545 and 700 nm for Mn-hematites. There is evidence for a approximately 700 nm band in spectral data for Martian bright regions, but association of it with Mn(3+) is not a unique interpretation. Comparison of laboratory and Martian spectral data implies that Mn(s) less than 0.032 for the Mn(3+) content of Martian hematites.

Description: Ferric-iron-bearing materials play an important role in the interpretation of visible to near-IR Mars spectra, and they may play a similarly important role in the analysis of new mid-IR spacecraft spectral observations to be obtained over the next decade. We review exisiting data on mid-IR transmission spectra of ferric oxides/oxyhydroxides and present new transmission spectra for ferric-bearing materials spanning a wide range of mineralogy and crystallinity. These materials include 11 samples of well-crystallized ferric oxides (hematite, maghemite, and magnetite) and ferric oxyhydroxides (goethite, lepidocrocite). We also report the first transmission spectra for purely nanophase ferric oxide samples that have been shown to exhibit spectral similarities to Mars in the visible to near-IR and we compare these data to previous and new transmission spectra of terrestial palagonites. Most of these samples show numerous, diagnostic absorption features in the mid-IR due to Fe(3+) - 0(2-) vibrational transitions, structural and/or bound OH, and/or silicates. These data indicate that high spatial resolution, moderate spectral resolution mid-IR ground-based and spacecraft observations of Mars may be able to detect and uniquely discriminate among different ferric-iron-bearing phases on the Martian surface or in the airborne dust.

Description: This paper focuses on the comparison of cloud amounts derived from an atmospheric general circulation model (AGCM), Satellite-observed clouds, and Ground-based cloud observations. Unlike Earth Radiation Budget Experiment (ERBE)-type comparisons it does not mix potential errors in the cloud amount with those in the radiation code embedded in the model. Long term cloud climatologies were used to compare global cloud amounts and regional seasonal cycles. The AGCM successfully reproduced the signatures of the warm pool and North Pacific seasonal cycle cloudiness but failed in the low stratus region off the coast of South America, a known problem for AGCMs. The data sets also reproduced the anomaly signature associated with El Nino in the warm pool region, but the model amounts were lower. Global results had a similar success rate, with the model generally producing lower total cloud compared to the satellite and in situ measurements. To compare cloud vertical distributions the cloud height may need to be validated using the corresponding radiation fields. Unfortunately there were also some large discrepancies between the two observed cloud data sets. While tremendously improved over the last decade the character of the observed cloud data sets, must be substantially enhanced before they will be useful in validating AGCMs by any but the crudest levels of comparison.

Description: Movement of water as vapor in the atmosphere is a fundamental process in the Earth's hydrological cycle. Investigations of spatial and time scales of water vapor transport in the atmosphere are important areas of research. Water vapor transmits energy as a function of its abundance across the spectrum. This is shown in the 400- to 2500-nm spectral region where the transmission of the terrestrial atmosphere has been modeled using the MODTRAN radiative transfer code for a range of water vapor abundances. Based on these model results, spectra measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) have been used to investigate the movement of water vapor at 20-m spatial resolution over an 11-by-30-km area at approximately 15-minute time intervals (1.25 hours total). AVIRIS measures the upwelling spectral radiance from 400 to 2500 nm at 10-nm spectral intervals and collects images of 11-by-up-to-1000 km at 20-m spatial resolution. Data are collected at a rate of 1 km of flight line per 4.5 seconds. A set of five AVIRIS flight lines was acquired in rapid succession over Rogers Dry Lake, CA on May 18, 1993 at 18:59, 19:13, 19:29, 19:47, and 19:59 UTC. Rogers Dry lake is located 2 hours north of Los Angeles, CA at 34.84 degrees north latitude and 117.83 degrees west longitude in the Mojave Desert.

Description: The Voyager 1 spacecraft flew by Jupiter on March 5, 1979. Spacecraft navigation was performed with radio tracking data from NASA's Deep Space Network. In the years since then, there has been a great deal of progress in the definition of celestial reference frames and in determining the orbit and orientation of the Earth. Using these improvements, the radio metric range and Doppler data acquired from the Voyager 1 spacecraft near its encounter with Jupiter have been reanalyzed to determine the plane-of-sky position of Jupiter with much greater accuracy than was possible at the time of the encounter. The position of Jupiter at the time of encounter has been determined with an accuracy of 40 nrad in right ascension and 140 nrad in declination with respect to the celestial reference frame defined by the International Earth Rotation Service. This position estimate has been done to improve the ephemeris of Jupiter prior to the upcoming encounter of the Galileo spacecraft with Jupiter.

Description: On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Description: We have carried out three (piggyback) radon-related projects aboard the KAO. The first, which was limited to upper tropospheric measurements while in level flight, revealed the systematic occurrence of unexpectedly high radon concentrations in this region of the atmosphere. The second project was an instrument development project, which led to the installation of an automatic radon measurement system aboard the NASA ER-2 High Altitude Research Aircraft. In the third, we installed a new system capable of collecting samples during the normal climb and descent of the KAO. The results obtained in these projects have resulted in significant contributions to our knowledge of atmospheric transport processes, and are currently playing a key role in the validation of global circulation and transport models.

Description: A data base comprised of all available ionosphere satellite ion composition measurements - the Goddard Comprehensive Ionosphere Data Base (GCID) - has been set upon optical disks for convenient merging and accessing of data from different satellites. This data has recentely been expanded to include all accessible satellite electron density and plasma temperature measurements. This paper demonstrates with a couple examples, the potential of GCID as an alternative to empirical models for undertanding ionsophere physics and chemistry. Through binning techniques GCID provdes a tool that complements and improves on the International Reference Ionosphere model (IRI) in delineating the topology of high latitude ion composition. GCID provides at a glance a measure of the local statistical variabilty of the ion compostions that is not in empirical models such as IRI. Furthermusing the data to determine the spatial and geophysical parameter range over which the minor ion species are approximately in chemical equilibrium, the number of data points available for empirical ion composition models can be increased by using the statiscally more reliable neutral and electron empirical models to derive ions composition. Currently available empirical models are not capable of adequately defining the comples high latitude distribution - a measure of the variablity is needed.

Description: An international conference on high-latitude ionospheric modeling produced 27 papers in the areas of ionospheric mapping, electron density and distribution, ion density and distribution, ionospheric storems, ionospheric composition, and ionospheric sounding techniques. Upgrades to the International Reference Ionosphere (IRI) model were proposed in several papers.

Description: Two symposia related to the exploration of the Venus and Mars atmospheres produced papers which described the results of spacecraft-based and ground-based observations of these two planets. The topics addressed included solar wind interactions, molecular ions, planetary magnetospheres, planetary ionospheres, planetary thermospheres, gravity waves, the homopause, radiative transfer, radiative balance, atmospheric water, planetary lightning, very low frequency bursts, cosmic-ray ionization, planetary loss, and future missions to study Mars aeronomy, energetic particles near Mars, radar exploration of Mars, and the Mars Environmental Survey (MESUR) mission.

Description: At the last International Reference Ionosphere (IRI) Workshop it was decided that future editions of the model should include a representationof the aurroral oval boundaries. In this paper we review the different existing parameterizations of the auroral oval discussing their data base, boundary criteria, matematical formation, and overall usefulness for IRI. As a first candidate for incorporation into IRI we recomment the parameterization of the Feldstein ovals by Holzworth and Meng. Ways of implementing this model into IRI are discussed. We will also address adjustability with user-provided boundaries or boundary-related parameters, to better support strom-related studies.

Description: A newly discovered source of energetic particle streams directed at Venus' nightside may cause several exosphere phenomena. The streams evolve from dispersed cosmoids, a surreptitious population of meteoroids in nearly hyperbolic orbits, measured with three dust experiments on Pioneer 10/11. Loose fragile comet-like ensembles of frozen volatile material, they are orders of magnitude more populous than short-period meteoroids. Dispersion is forced near a planet masking the meteor signature at Earth, but recently terrestrial exosphere interactions have been detected principally from VLF radar returns and as suddenly formed layers of neutral sodium and iron probed by lidar. Venus influx is estimated greater than 10(exp -14)(g)/((sq cm)(s)) with nightside directed kinetic power of greater than 0.3 (erg)/((sq cm)s)). Compared to Earth, Venus, lacking a dipole field with almost no rotation, has a readily recognized set of nightside interaction signatures of these downward energetic particle streams: i.e., density depressions of the neutral thermosphere, hydrogen and helium bulges, a non-disappearing ionosphere, electron holes, nightglows, VLF bursty signals, and local magnetic fields. From Pioneer Venus Orbiter (PVO) measurements, the cosmoid population, its temporal and solar azimuthal variation may be determined.

Description: High-resolution, near-infrared (1.09 to 2.5 micrometers) spectra of the night side of Venus have been obtained in 1990 and 1991 using the Fourier Transform Spectrometer at the 3.6-m Canada-France-Hawaii telescope. Absorptions due to H2O were detected in spectral windows near 2.3, 1.74, and 1.18 micrometers. Our analysis of these absorptions constrains the abundance of water vapor in three different altitude ranges located between the clouds and the surface: 30-40 km, 15-25 km and 0-15 km. A constant water vapor mixing ratio of 30 +/- 15 ppm below the clouds can fit the observations. These values are consistent with recent near-infrared studies of the night side of Venus at lower spectral resolution. The atmosphere of Venus appears to be dryer than originally suggested by the in-situ measurements made by the Pioneer Venus and Venera mass-spectrometers and gas-chromatographs.

Description: Using Pioneer Venus data, hydrogen and deuterium ions are shown to escape from the hydrogen bulge region in the nightside ionosphere. The polarization electric field propels these light ions upward through the ionosphere and into the ion-exosphere, where H(+) and D(+) continue to be accelerated away from Venus and move into the ionotail and beyond. The vertical flow speeds of H(+) and D(+) are found to be about the same; therefore, selective escape between H(+) and D(+) is negligible for this mechanism. Present day planetary loss rates of about 8.6 x 10(exp 25)/s and 3.2 X 10(exp 23)/s were obtained for H(+) and D(+), respectively. Such rates, persisting over a few billion years, should have significantly affected the planetary water budget.

Description: This paper was presented during the special session that was held at the 1993 International Reference Ionosphere (IRI) Workshop in honor of Karl Rawer's 80th birthday. It retraces the steps that led from the start of the IRI project to the present edition of the model highlighting the important role that the honoree played in guiding IRI from infancy to maturity. All summary view graphs are reproduced at the end of the article.

Description: A description of the electric structure as observed during a rocket sounding noctilucent cloud (NLC) observation program in the summer of 1991, is presented. Both NLC and polar mesosphere summer echo (PMSE) conditions were determined. The observable effects on ions and electrons were measured in the NLC and PMSE regions and associated electric field (E-fields) measurements indicated small alternating current vertical fields. The following findings are reported: small and large scale electron density structure is noticeable in the NLC regions; above 82 km, the occurrence of very low-mobility positive ions directly demonstrates the existence of charged aerosols, and small alternating current vertical electric fields were observed in the NLC and PMSE region.

Description: The mesospheric and lower thermospheric equatorial dynamics program (MALTED), conducted from the Alcantara, Brazil rocket site as part of the international Guara rocket campaign to study equatorial dynamics, irregularities, and instabilities, is reported. The MALTED program was concerned with the 16-day modulation of the diurnal tidal amplitude, which shows high extremes during August in the equatorial belt. The interaction of this global phenomenon with locally produced gravity waves and turbulence in the mesopause region is studied in order to gain a better understanding of the dynamic influences on the equatorial middle atmosphere. Four identical payloads designed to investigate small-scale turbulence and irregularities, were coordinated with 20 falling sphere rockets designed to measure meteorological parameters. The prediction and monitoring of global mesospheric effects were obtained through coordination with various ground-based radar observatories. The campaign logistics, the instrumentation, and the preliminary results are described.

Description: A spaceborne solar UV irradiance observation from NOAA-14 polar orbiting satellite by using an 1/4 meter double-ebert spectrometer, is reported.

Description: Terrestrial far ultraviolet (FUV) airglow emissions have been suggested as a means for remote sensing the structure of the upper atmosphere. The energy which leads to the excitation of FUV airglow emissions is solar irradiance at extreme ultraviolet (EUV) and soft x-ray wavelengths. Solar irradiance at these wavelengths is known to be highly variable; studies of nitric oxide (NO) in the lower thermosphere have suggested a variability of more than an order of magnitude in the solar soft x-ray irradiance. To properly interpret the FUV airflow, the magnitude of the solar energy deposition must be known. Previous analyses have used the electron impact excited Lyman-Birge-Hopfield (LBH) bands of N2 to infer the flux of photoelectrons in the atmosphere and thus to infer the magnitude of the solar irradiance. This dissertation presents the first simultaneous measurements of the FUV airglow, the major atmospheric constituent densities, and the solar EUV and soft x-ray irradiances. The measurements were made on three flights of an identical sounding rocket payload at different levels of solar activity. The linear response in brightness of the LBH bands to variations in solar irradiance is demonstrated. In addition to the N2 LBH bands, atomic oxygen lines at 135.6 and 130.4 nm are also studied. Unlike the LBH bands, these emissions undergo radiative transfer effects in the atmosphere. The OI emission at 135.6 nm is found to be well modeled using a radiative transfer calculation and the known excitation processes. Unfortunately, the assumed processes leading to OI 130.4 nm excitation are found to be insufficient to reproduce the observed variability of this emission. Production of NO in the atmosphere is examined; it is shown that a lower than previously reported variability in the solar soft x-ray irradiance is required to explain the variability of NO.

Description: Simultaneous very long baseline interferometry (VLBI) and water vapor radiometer (WVR) measurements on a 21 km baseline showed that calibration by WVRs removed a significant fraction of the effect of tropospheric delay fluctuations for these experiments. From comparison of the residual delay variations within scans and between scans, the total tropospheric contribution t the delay residuals for each of the three 5 to 20 hour sessions was estimated as 1, 17, and 10%, with the first value being uncertain. The observed improvement in rms residual delay from WVR calibration during these three sessions was 4, 16, and 2%, respectively. The improvement is consistent with the estimated 2 to 3 mm path delay precision of current WVRs. The VLBI measurements, of natural radio sources, were conducted in April and May 1993 at Goldstone, California. Dual-frequency (2.3 and 8.4 GHz) observations were employed to remove the effects of charged particles from the data. Measurements with co-pointed WVRs, located within 50 m of the axis of each antenna, were performed to test the ability of the WVRs to calibrate line-of-sight path delays. Factors that made WVR performance assessment difficult included (1) the fact that the level of tropospheric fluctuations was smaller than is typical for Goldstone during these experiments and (2) VLBI delay variations on longer time scales (i.e., over multiple scans) contained uncalibrated instrumental effects (probably a result of slow temperature variations in the VLBI hardware) that were larger than the tropospheric effects.

Description: Substantial evidence suggests that a UV spectrally Absorbing Material (UV-SAM) exists on Triton's surface. This evidence is found in the positive slope in Triton's spectrum from the UV to the near-IR, and the increasing contrast in Triton's light curve in the blue and UV. Although it is now widely-thought that UV-SAMs exist on Triton, little is known about their distribution and spectral properties. The goal of this NDAP Project is to determine the spatial distribution and geological context of the UV-SaM material. We hope to determine if UV-SAMs on Triton are correlated with geologic wind streaks, craters, calderas, geomorphic/topographic units, regions containing (or lacking) volatile frosts, or some other process (e.g., magnetospheric interactions). Once the location and distribution of UV-SAMs has been determined, further constraints on their composition cable made by analyzing the spectrographic data set. To accomplish these goals, various data sets will be used, including Voyager 2 UV and visible images of Triton's surface, IUE and HST spectra of Triton, and a geologic map of the surface based on voyager 2 and spectrophotometric data. The results of this research will be published in the planetary science literature.

Description: For three decades, magnetospheric field and plasma measurements have been made by diverse instruments flown on spacecraft in may different orbits, widely separated in space and time, and under various solar and magnetospheric conditions. Scientists have used this information to piece together an intricate, yet incomplete view of the magnetosphere. A simultaneous global view, using various light wavelengths and energetic neutral atoms, could reveal exciting new data nd help explain complex magnetospheric processes, thus providing a clear picture of this region of space. This report documents the scientific rational for such a magnetospheric imaging mission and provides a mission concept for its implementation.

Description: The goal of the Atmospheric Response in Aurora (ARIA) experiment carried out at Poker Flat, Alaska, on March 3, 1992, was to determine the response of the neutral atmosphere to the long-lived, large-scale forcing that is characteristic of the diffuse aurora in the post midnight sector. A combination of chemical release rocket wind measurements, instrumented rocket composition measurements, and ground-based optical measurements were used to characterize the response of the neutral atmosphere. The rocket measurements were made at the end of a 90-min period of strong Joule heating. We focus on the neutral wind measurements made with the rocket. The forcing was determined by running the assimilated mapping of ionospheric electrodynamics (AMIE) analysis procedure developed at the National Center for Atmospheric Research. The winds expected at the latitude and longitude of the experiment were calculated using the spectral thermospheric general circulation model developed at the Danish Meteorological Institute. Comparisons of the observations and the model suggest that the neutral winds responded strongly in two height ranges. An eastward wind perturbation of approximately 100 m/s developed between 140 and 200 km altitude with a peak near 160 km. A southwestward wind with peak magnitude of approximately 150 m/s developed near 115 km altitude. The large amplitude winds at the lower altitude are particularly surprising. They appear to be associated with the upward propagating semidiurnal tide. However, the amplitude is much larger than predicted by any of the tidal models, and the shear found just below the peak in the winds was nominally unstable with a Richardson number of approximately 0.08.

Description: Work on completing our analysis of the Voyager UVS solar occultation data acquired during Neptune encounter is essentially complete, as testified by the attached poster materials. The photochemical modeling addresses the recent revision in branching ratios for radical production in the photolysis of methane at H Lyman alpha implied by the lab measurements of Mordaunt et al. (1993). The software generated in this effort has been useful for checking the degree to which photochemical models addressing other datasets (mainly infrared) are consistent with the UVS data. This work complements the UVS modeling results in that the IR data refer to deeper pressure levels; as regards the modeling of UVS data, the most significant result is the convincing support for the presence of a stagnant lower stratosphere. Evidence for strong dynamical (mixing) transport of minor constituents at shallower pressures is provided by the UVS data analysis.

Description: Moments of plasma distributions observed in the magnetotail vary with different time scales. In this paper we attempt to explain the observed variability on intermediate timescales of approximately 10-20 min that result from the simultaneous energization and spatial structuring of solar wind plasma in the distant magnetotail. These processes stimulate the formation of a system of spatially disjointed. highly accelerated filaments (beamlets) in the tail. We use the results from large-scale kinetic modeling of magnetotail formation from a plasma mantle source to calculate moments of ion distribution functions throughout the tail. Statistical restrictions related to the limited number of particles in our system naturally reduce the spatial resolution of our results, but we show that our model is valid on intermediate spatial scales Delta(x) x Delta(z) equal to approximately 1 R(sub E) x 1000 km. For these spatial scales the resulting pattern, which resembles a mosaic, appears to be quite variable. The complexity of the pattern is related to the spatial interference between beamlets accelerated at various locations within the distant tail which mirror in the strong near-Earth magnetic field. Global motion of the magnetotail results in the displacement of spacecraft with respect to this mosaic pattern and can produce variations in all of the moments (especially the x-component of the bulk velocity) on intermediate timescales. The results obtained enable us to view the magnetotail plasma as consisting of two different populations: a tailward-Earthward system of highly accelerated beamlets interfering with each other, and an energized quasithermal population which gradually builds as the Earth is approached. In the near-Earth tail, these populations merge into a hot quasi-isotropic ion population typical of the near-Earth plasma sheet. The transformation of plasma sheet boundary layer (PSBL) beam energy into central plasma sheet (CPS) quasi-thermal energy occurs in the absence of collisions or noise. This paper also clarifies the relationship between the global scale where an MHD description might be appropriate and the lower intermediate scales where MHD fails and large-scale kinetic theory should be used.

Description: Spectral observations providing evidence for the presence of volatile-bearing minerals on the surface of Mars were obtained in 1988 and 1990 from the KAO. The 1988 data suggest the presence of 1-3 weight percent (wt%) of carbonate/bicarbonate and 10-15 wt% sulfate/bisulfate associated with martian atmospheric dust. Estimates of the optical depths are approximately 0.60 and approximately 0.35 in 1988 and 1990, respectively.

Description: We use a dynamical guiding-center model to investigate the stormtime transport of ring current and radiation-belt ions. We trace the motion of representative ions' guiding centers in response to model substorm-associated impulses in the convection electric field for a range of ion energies. Our simple magnetospheric model allows us to compare our numerical results quantitatively with analytical descriptions of particle transport, (e.g., with the quasilinear theory of radial diffusion). We find that 10-145-keV ions gain access to L approximately 3, where they can form the stormtime ring current, mainly from outside the (trapping) region in which particles execute closed drift paths. Conversely, the transport of higher-energy ions (approximately greater than 145 keV at L approximately 3) turns out to resemble radial diffusion. The quasilinear diffusion coefficient calculated for our model storm does not vary smoothly with particle energy, since our impulses occur at specific (although randomly determined) times. Despite the spectral irregularity, quasilinear theory provides a surprisingly accurate description of the transport process for approximately greater than 145-keV ions, even for the case of an individual storm. For 4 different realizations of our model storm, the geometric mean discrepancies between diffusion coefficients D(sup sim, sub LL) obtained from the simulations and the quasilinear diffusion coefficient D(sup ql, sub LL) amount to factors of 2.3, 2.3, 1.5, and 3.0, respectively. We have found that these discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) can be reduced slightly by invoking drift-resonance broadening to smooth out the sharp minima and maxima in D(sup ql, sub LL). The mean of the remaining discrepancies between D(sup sim, sub LL) and D(sup ql, sub LL) for the 4 different storms then amount to factors of 1.9, 2.1, 1.5, and 2.7, respectively. We find even better agreement when we reduce the impulse amplitudes systematically in a given model storm (e.g., reduction of all the impulse amplitudes by half reduces the discrepancy factor by at least its square root) and also when we average our results over an ensemble of 20 model storms (agreement is within a factor of 1.2 without impulse-amplitude reduction). We use our simulation results also to map phase-space densities f in accordance with Liouville's theorem. We find that the stormtime transport of approximately greater than 145-keV ions produces little change in f-bar the drift-averaged phase-space density on any drift shell of interest. However, the stormtime transport produces a major enhancement from the pre-storm phase-space density at energies approximately 30-145 keV, which are representative of the stormtime ring current.

Description: This final report describes the activities of the Lockheed Martin Palo Alto Research Laboratory in studying the measurement of ion and electron precipitation induced Far Ultra-Violet (FUV) emissions and Geocoronal Lyman Alpha for the NASA Inner Magnetospheric Imager (IMI) mission. this study examined promising techniques that may allow combining several FUV instruments that would separately measure proton aurora, electron aurora, and geocoronal Lyman alpha into a single instrument operated on a spinning spacecraft. The study consisted of two parts. First, the geocoronal Lyman alpha, proton aurora, and electron aurora emissions were modeled to determine instrument requirements. Second, several promising techniques were investigated to determine if they were suitable for use in an IMI-type mission. Among the techniques investigated were the Hydrogen gas cell for eliminating cold geocoronal Lyman alpha emissions, and a coded aperture spectrometer with sufficient resolution to separate Doppler shifted Lyman alpha components.

Description: Near simultaneous observations of the Saturnian H2 north ultraviolet aurora and the polar haze were made at 153 nm and 210 nm respectively with the Faint Object Camera on board the Hubble Space Telescope. The auroral observations cover a complete rotation of the planet and, when co-added, reveal the presence of an auroral emission near 80 deg N with a peak brightness of about 150 kR of total H2 emission. The maximum optical depth of the polar haze layer is found to be located approximately 5 deg equatorward of the auroral emission zone. The haze particles are presumably formed by hydrocarbon aerosols initiated by H2+ auroral production. In this case, the observed haze optical depth requires an efficiency of aerosol formation of about 6 percent, indicating that auroral production of hydrocarbon aerosols is a viable source of high-latitude haze.

Description: The second pulsating aurora (PULSAUR 2) rocket was a comprehensively instrumented rocket for the study of the generating mechanisms of pulsating auroras and their related atmospheric and ionospheric effects. It carried instruments to measure electrons and ions over a wide range of energies and pitch angles, optical emissions and X-rays, direct currents and alternating currents, magnetic fields, electron densities and temperatures. The rocket was flown from Andoeya, Norway on 9 February 1994, during good auroral and optical conditions, and reached an altitude of 291.5 km. Coordinated measurements were made from the rocket range by television cameras, photometers, magnetometers, riometers and very low frequency receivers. Ionospheric parameters were measured by the EISCAT radar. The main objectives and performance of the project are discussed, and some of the data obtained are mentioned.

Description: The goals of this project are to observe how the open water/thin ice fraction in a high-concentration ice pack responds to different short-period atmospheric forcings, and how this response is represented in different scales of observation. The objectives can be summarized as follows: determine the feasibility and accuracy of ice concentration and ice typing by ERS-1 SAR backscatter data, and whether SAR data might be used to calibrate concentration estimates from optical and massive-microwave sensors; investigate methods to integrate SAR data with other satellite data for turbulent heat flux parameterization at the ocean/atmosphere interface; determine how the development and evolution of open water/thin ice areas within the interior ice pack vary under different atmospheric synoptic regimes; compare how open-water/thin ice fractions estimated from large-area divergence measurements differ from fractions determined by summing localized openings in the pack; relate these questions of scale and process to methods of observation, modeling, and averaging over time and space.

Description: The use of radio sounding techniques for the study of the ionospheric plasma dates back to G. Briet and M. A. Tuve in 1926. Ground based swept frequency sounders can monitor the electron number density (N(sub e)) as a function of height (the N(sub e) profile). These early instruments evolved into a global network that produced high-resolution displays of echo time delay vs frequency on 35-mm film. These instruments provided the foundation for the success of the International Geophysical Year (1958). The Alouette and International Satellites for Ionospheric Studies (ISIS) programs pioneered the used of spaceborne, swept frequency sounders to obtain N(sub e) profiles of the topside of the ionosphere, from a position above the electron density maximum. Repeated measurements during the orbit produced an orbital plane contour which routinely provided density measurements to within 10%. The Alouette/ISIS experience also showed that even with a high powered transmitter (compared to the low power sounder possible today) a radio sounder can be compatible with other imaging instruments on the same satellite. Digital technology was used on later spacecraft developed by the Japanese (the EXOS C and D) and the Soviets (Intercosmos 19 and Cosmos 1809). However, a full coherent pulse compression and spectral integrating capability, such as exist today for ground-based sounders (Reinisch et al., 1992), has never been put into space. NASA's 1990 Space Physics Strategy Implementation Study "The NASA Space Physics Program from 1995 to 2010" suggested using radio sounders to study the plasmasphere and the magnetopause and its boundary layers (Green and Fung, 1993). Both the magnetopause and plasmasphere, as well as the cusp and boundary layers, can be observed by a radio sounder in a high-inclination polar orbit with an apogee greater than 6 R(sub e) (Reiff et al., 1994; Calvert et al., 1995). Magnetospheric radio sounding from space will provide remote density measurements of unprecedented precision and coverage in the plasmasphere, inner magnetosphere and magnetopause, from which the structure, inter-relationship, and variations of different plasma regions can be determined (Armstrong Johnson, 1995). A space-borne Radio Plasma Imager (RPI) could provide a unique global view of the magnetosphere revealing the underlying structure of remote plasma regions, thereby providing a framework for the interpretation of images obtained by other techniques as identified in the technical areas TA1 to TA4 in the MSFC NRA8-8.

Description: One year of near-continuous water vapor radiometer (WVR) measurements at DSS 13 has provided a database for characterizing the Goldstone tropospheric delay properties in a statistical sense. The results have been expressed in terms of the Allan standard deviation of delay and compared to a previous model for Goldstone fluctuations and the specifications of the Cassini Gravitational Wave Experiment (GWE). The new WVR data indicate that average fluctuation levels at hour time scales or less are approximately 30 percent lower than the earlier Goldstone model predictions. At greater than 1 h time scales, the WVR indicated fluctuation levels are in closer agreement with the model, although noise floor limitations may be artificially raising the average WVR-derived atmospheric fluctuation levels at the longer time scales. When scaled to two-way Doppler tracking at 20 deg elevation, as will occur for the GWE, these results indicate that Goldstone winter tropospheric delay fluctuations will typically be a factor of 10 larger than the GWE requirements at 1000 s and a factor of 4 larger at 10,000 s.

Description: The main results of the grant were (1) finishing the manuscript of a proof of completeness of the Poincare modes in an incompressible nonviscous fluid corotating with a rigid ellipsoidal boundary, (2) partial completion of a manuscript describing a definition of helicity that resolved questions in the literature about calculating the helicities of vector fields with complicated topologies, and (3) the beginning of a reexamination of the inverse problem of inferring properties of the geomagnetic field B just outside the core-mantle boundary (CMB) from measurements of elements of B at and above the earth's surface. This last work has led to a simple general formalism for linear and nonlinear inverse problems that appears to include all the inversion schemes so far considered for the uniqueness problem in geomagnetic inversion. The technique suggests some new methods for error estimation that form part of this report.

Description: The various regions of the magnetosphere-ionosphere system are coupled by flows of charged particle beams and electromagnetic waves. This coupling gives rise to processes that affect both technical and non-technical aspects of life on Earth. The CRRES Program sponsored experiments which were designed to produce controlled and known input to the space environment and the effects were measured with arrays of diagnostic instruments. Large amounts of material were used to modify and perturb the environment in a controlled manner, and response to this was studied. The CRRES and PEGSAT satellites were dual-mission spacecraft with a NASA mission to perform active chemical-release experiments, grouped into categories of tracer, modification, and simulation experiments. Two sounding rocket chemical release campaigns completed the study.

Description: Among the hazards that must be negotiated by lunar-landing spacecraft are blocks on the surface of the Moon. Unfortunately, few data exist that can be used to evaluate the threat posed by such blocks to landing spacecraft. Perhaps the best information is that obtained from Surveyor photographs, but those data do not extend to the dimensions of the large blocks that would pose the greatest hazards. Block distributions in the vicinities of the Surveyor 1, 3, 6, and 7 sites have been determined from Lunar Orbiter photography and are presented here. Only large (i.e., greater than or equal to 2.5 m) blocks are measurable in these pictures, resulting in a size gap between the Surveyor and Lunar Orbiter distributions. Nevertheless, the orbital data are self-consistent, a claim supported by the similarity in behavior between the subsets of data from the Surveyor 1, 3, and 6 sites and by the good agreement in position (if not slopes) between the data obtained from the Surveyor 3 photography and those derived from the Lunar Orbiter photographs. Confidence in the results is also justified by the well-behaved distribution of large blocks at the surveyor site. Comparisons between the Surveyor distributions and those derived from the orbital photography permit these observations: (1) in all cases but that for Surveyor 3, the density of large blocks is overestimated by extrapolation of the Surveyor-derived trends; (2) the slopes of the Surveyor-derived distributions are consistently lower than those determined for the large blocks; and (3) these apparent disagreements could be mitigated if the overall shapes of the cumulative lunar block populations were nonlinear, allowing for different slopes over different size intervals. The relatively large gaps between the Surveyor-derived and Orbiter-derived data sets, however, do not permit a determination of those shapes.

Description: The Space Studies Board of the National Research Council charged its Committee on Planetary and Lunar Exploration (COMPLEX) to (1) examine the degree to which small missions, such as those fitting within the constraints of the Discovery program, can achieve priority objectives in the lunar and planetary sciences; (2) determine those characteristics, such as level of risk, flight rate, target mix, university involvement, technology development, management structure and procedures, and so on, that could allow a successful program; (3) assess issues, such as instrument selection, mission operations, data analysis, and data archiving, to ensure the greatest scientific return from a particular mission, given a rapid deployment schedule and a tightly constrained budget; and (4) review past programmatic attempts to establish small planetary science mission lines, including the Planetary Observers and Planetary Explorers, and consider the impact management practices have had on such programs. A series of small missions presents the planetary science community with the opportunity to expand the scope of its activities and to develop the potential and inventiveness of its members in ways not possible within the confines of large, traditional programs. COMPLEX also realized that a program of small planetary missions was, in and of itself, incapable of meeting all of the prime objectives contained in its report 'An Integrated Strategy for the Planetary Sciences: 1995-2010.' Recommendations are provided for the small planetary missions to fulfill their promise.

Description: Recently, extremely interesting high-altitude cloud-ionosphere electrical discharges, like lightning above thunderstorms, have been observed from NASA's space shuttle missions and during airborne and ground-based experiments. To understand these discharges, a new experiment was conceived to simulate a thundercloud in a vacuum chamber using a dielectric in particulate form into which electrodes were inserted to create charge centers analogous to those in an electrified cloud. To represent the ionosphere, a conducting medium (metallic plate) was introduced at the top of the chamber. It was found that for different pressures between approximately 1 and 300 mb, corresponding to various upper atmospheric altitudes, different discharges occurred above the simulated thundercloud, and these bore a remarkable similarity to the observed atmospheric phenomena. At pressures greater than 300 mb, these discharges were rare and only discharges within the simulated thundercloud were observed. Use of a particulate dielectric was critical for the successful simulation of the high-altitude lightning.

Description: The objectives are to look at time scales of overpressure compared to cratering and to determine: what are the transient pressure and density due to atmospheric entry; do shock waves evacuate ambient gas; do transient atmospheric disturbances 'settle down' during cratering; can the pressure/density field be approximated as quasi-static; how does disturbance scale with impactor size; and what is the role of atmospheric thickness. The general approach is to perform inexpensive exploratory calculations, perform experiments to validate code and observe crater growth, and to follow up with more realistic coupling calculations. This viewgraph presentation presents progress made with the objective to obtain useful scaling relationships for crater formation when atmospheric effects are important.

Description: The past few decades have witnessed significant advances in our understanding of how stars form, and there has been an associated increase in our knowledge of conditions and phenomena in the early solar system. These have led to the formulation of a paradigm for the origin of the solar system that is sufficiently complete that its basic elements can be tested directly through observations. A simple, but profound, consequence of the paradigm is that most if not all stars should be accompanied by planetary systems. The accuracy of instruments that can be used in such searches has improved to the point that Jupiter-like companions to a number of nearby stars could be detected. However, the results to date are that no other planetary systems have been detected, and the absence of detection is becoming statistically significant, particularly as it relates to the existence of brown dwarf companions to main-sequence stars.

Description: The latest version of the Global Reference Atmospheric Model (GRAM-95) is presented and discussed. GRAM-95 uses the new Global Upper Air Climatic Atlas (GUACA) CD-ROM data set, for 0- to 27-km altitudes. As with earlier versions, GRAM-95 provides complete geographical and altitude coverage for each month of the year. Individual years 1985 to 1991 and a period-of-record (1980 to 1991) can be simulated for the GUACA height range. GRAM-95 uses a specially developed data set, based on Middle Atmosphere Program (MAP) data, for the 20- to 120-km height range, and the NASA Marshall Engineering Thermosphere (MET) model for heights above 90 km. Fairing techniques assure a smooth transition in the overlap height ranges (20 to 27 km and 90 to 120 km). In addition to the traditional GRAM variables of pressure, density, temperature and wind components, GRAM-95 now includes water vapor and 11 other atmospheric constituents (O3, N2O, CO, CH4, CO2, N2, O2, O, A, He, and H). A new, variable-scale perturbation model provides both large-scale and small-scale deviations from mean values for the thermodynamic variables and horizontal and vertical wind components. The perturbation model includes new features that simulate intermittency (patchiness) in turbulence and small-scale perturbation fields. The density perturbations and density gradients (density shears) computed by the new model compare favorably in their statistical characteristics with observed density perturbations and density shears from 32 space shuttle reentry profiles. GRAM-95 provides considerable improvement in wind estimates from the new GUACA data set, compared to winds calculated from the geostrophic wind relations previously used in the 0- to 25-km height range. The GRAM-95 code has been put into a more modular form, easier to incorporate as subroutines in other programs (e.g., trajectory codes). A complete user's guide for running the program, plus sample input and output, is provided.

Description: Measurements by single-frequency satellite altimeter (Geosat, ERS-1) require a ionospheric correction to account for the signal time delay in the ionosphere. We propose using the International Reference Ionosphere (IRI) for the determination of this time delay. To investigate the effectiveness of an IRI correction, we have compared the IRI values with ionospheric corrections deduce from measurements by the dual-frequency Topex altimeter. By measuring at two frequencies, the Topex instrument can record (and thus eliminate) the ionospheric influence. We find that IRI agrees with the Topex data much better than the model that is currently used in Geosat data analysis. In particular the earlier model does not represent the equatorial double-peak (equator anomaly) clearly seen in the Topex data. Overall, the use of IRI results in a 30% improvement (over the older model) in the accuracy of ionospheric corrections computed for the first year of the Topex mission.

Description: Papers from the conference are presented and cover the following topics: ion and electron beams; ionospheric modification; spacecraft interactions; chemical releases; and plasma waves. Auroras and plasma emissions are reported from electron beam injection experiments on the EXCEDE 3 rocket and APEX satellite respectively. The important parameters affecting the charging of spacecraft during the operation of electron guns is covered. The Active Magnetospheric Particle Acceleration Satellite (AMPAS) mission utilizing dual-payload tethered satellites and both up and downward directed electron beams is proposed to study the magnetosphere. Recent results and associated theories from the Sura, Arecibo and Troms ionospheric heating facitlites are presented. The effects of neutral gases on spacecraft charging are examined in a series of rocket flights. Many results from the Combined Release and Radiation Effects Satellite chemical release experiments are presented.

Description: The main subject of the symposium was 'off-median phenomena'. The title denotes a range of problems that are rarely considered in 'pure science' studies of the ionosphere. The appearance of regular ionospheric variations is well known. Most of these depend on evident solar-geophysical influences like day and night, solar zenith angle, the seasons, geomagnetic control, solar activity, etc. Applicants and theoreticians as well used to work with monthly medians so that the in fact existing and quite important day-to-day variability is systematically overlooked. It is evident that a descriptive model like the International Reference Ionosphere (IRI) would be inadequate if this variability were denied. Interesting contributions from the symposium on 'off-median phenomena' and the IRI are presented.

Description: This serves as a final report entitled Energization and Transport of Ions of Ionospheric Origin in the Terrestrial Magnetosphere. The work has been predominantly focused on ion outflows identified in two data sets: (1) Prognoz 7; and (2) Dynamics Explorer. The study analyzed ion densities, temperatures, and flow velocities in the magnetotail. The work performed under this contract consisted of developing a program to load the raw data, compute the background subtraction of a strong sun pulse, and use the net counts to calculate the low order moments of the distribution function. The study confirms the results of ISEE that the the cusp is a major source of plasmasheet plasma and goes beyond this to discuss the use of ion velocities as a way to examine the motions of the magnetotail.

Description: The Compton Observatory commonly encounters fluxes of energetic electrons which have been scattered from the inner radiation belt to the path of the satellite by resonant interactions with VLF waves from powerful man-made transmitters. The present investigation was motivated by the fact that in the Fall of 1993, the Gamma Ray Observatory was boosted from a 650 km altitude circular orbit to a 750 km altitude circular orbit. This was an opportunity, for the first time, to make observations at two different altitudes using the same instrument. We have examined DISCLA data from the Burst & Transient Source Experiment (BATSE) experiment from 1 Sep. 1993 to 29 Jan. 1994. During the period of study we identified 48 instances of the satellite encountering a cloud of energetic electrons which had been scattered by VLF transmitters. We find that boosting the altitude of the circular orbit from 650 km to 750 km increased the intensity of cyclotron resonance scattered electrons by a factor of two. To search for long term changes in the cyclotron resonance precipitation, we have compared the approx. 750 km altitude data from 106 days at the end of 1993 with data at the same altitudes and time of year in 1991. The cyclotron resonance events in 1991 were three times more frequent and 25% of those cases were more intense than any seen in the 1993 data. We attribute this difference to increased level of geomagnetic activity in 1991 near the Solar Maximum.

Description: Mesoscale model simulations provide insight into the complex jet streak adjustments on 11-12 July 1981 that preceded the first of two significant gravity wave events to have been generated over the Rocky Mountains in Montana. Simulations employing a variety of terrain treatments indicate that prior to wave formation, geostrophic adjustment processes modified the structure of the mid-upper tropospheric jet streak by creating secondary jetlets to the southeast of the polar jet streak in proximity to the gravity wave generation region. This simulated restructuring of the mid-upper tropospheric jet streak is the result of a four stage process. During stage 1, the wind adjusts to the mass field as the jet streak exit region propagates into the inflection point between the upstream trough and downstream ridge in the height field. Stage 2 is initiated as the mass field is forced to adjust to the new ageostrophic wind field created during stage 1. Stage 3 is defined by a second geostrophic adjustment process occurring in a similar manner but to the south and east of the adjustment which occurs during stage 1. A low-level mesoscale jetlet is formed during stage 4 in response to the low-level pressure falls that are established during stage 3. The perturbation of this jetlet, caused by orographically-induced adiabatic and diabatic physical processes, is the likely mechanism responsible for the generation of the first and second episode of observed gravity waves. The dynamics responsible for this wave episode are discussed as differential surface sensible heating inducing an orographically-forced mountain-plains solenoid, resulting in the formation of additional mesoscale jetlets and internal gravity waves. Also discussed is how convective latent heating modifies the numerically simulated terrain-induced internal gravity waves, especially their amplitude and phase velocities, which provide better agreement with those wave characteristics observed in nature. Finally, the three-dimensional linear response of a zonally uniform barotropic flow in a vertically unbounded, continuously stratified, Boussinesq atmosphere which is perturbed from geostrophic equilibrium is investigated.

Description: The retarding ion mass spectrometer (RIMS) experiment onboard the Dynamics Explorer 1 (DE 1) satellite was designed to perform energy and mass-per-charge analysis on low-energy ions (less than 50 eV) with mass/charge ratios ranging from 1 to 40 amu/Z. the DE 1 satellite, carrying the RIMS experiment, was launched into an elliptical polar orbit on August 3, 1981. The approximately 7.5 hour orbit has perigee of 675 km altitude and apogee of 24,875 km altitude. This document and those that follow in this series, contain summary RIMS data spectrograms for each orbit for which RIMS data are available. The RIMS instrument began returning science data on day 280 of 1981 and continued to return usable data until the end of the DE mission in March 1991. It should be noted that studies of the RIMS data set should be conducted only with a thorough awareness of the material described in the introduction section presented here, or in collaboration with a scientist familiar with RIMS data analysis.

Description: If current plans are realized, within the next few years, an extraordinary set of coordinated research efforts focusing on energy flows in the Arctic will be implemented. All are motivated by the prospect of global climate change. SHEBA (Surface Energy Budget of the Arctic Ocean), led by the National Science Foundation (NSF) and the Office of Naval Research (ONR), involves instrumenting an ice camp in the perennial Arctic ice pack, and taking data for 12-18 months. The ARM (Atmospheric Radiation Measurement) North Slope of Alaska and Adjacent Arctic Ocean (NSA/AAO) Cloud and Radiation Testbed (CART) focuses on atmospheric radiative transport, especially in the presence of clouds. The NSA/AAO CART involves instrumenting a sizeable area on the North Slope of Alaska and adjacent waters in the vicinity of Barrow, and acquiring data over a period of about 10 years. FIRE (First ISCCP (International Satellite Cloud Climatology Program) Regional Experiment) Phase 3 is a program led by the National Aeronautics and Space Administration (NASA) which focuses on Arctic clouds, and which is coordinated with SHEBA and ARM. FIRE has historically emphasized data from airborne and satellite platforms. All three program anticipate initiating Arctic data acquisition during spring, 1997. In light of his historic opportunity, the authors discuss a strawman atmospheric radiative transfer experimental plan that identifies which features of the radiative transport models they think should be tested, what experimental data are required for each type of test, the platforms and instrumentation necessary to acquire those data, and in general terms, how the experiments could be conducted. Aspects of the plan are applicable to all three programs.

Description: A survey of bulk parameters of analyzable O(+) outward streams in the mid-altitude (3-4.7 R(sub E) geocentric distance) polar cap magnetosphere is obtained from measurements by the Retarding Ion Mass Spectrometer (RIMS) aboard the Dynamics Explorer-1 (DE-1) spacecraft. There is wide scatter in the obtained densities, but they do display discernible trends: the average O(+) density in these streams decreases from about 60 ions/cc at 3.5 R(sub E) to about 1 ion/cc at 4.6 R(sub E). The streaming velocities are somewhat more defined, and their average increases from about 8 km/s at 3.5 R(sub E) to about 12 km/s at 4.6 R(sub E). The densities and bulk velocities are inversely correlated. We have further compared these observational trends with model profiles for the centrifugally-accelerated polar wind as recently described by Horwitz et al. (1994). The large outflow velocities observed can be understood in part as centrifugally-driven by convection with ionospheric electric field magnitudes of the order 50-70 mV/m, perhaps including plasma expansion effects.

Description: This report summarizes the current JPL efforts of generating a Mars gravity field from Viking 1 and 2 and Mariner 9 Doppler tracking data. The Mars 50c solution is a complete gravity field to degree and order 50 with solutions as well for the gravitational mass of Mars, Phobos, and Deimos. The constants and models used to obtain the solution are given and the method for determining the gravity field is presented. The gravity field is compared to the best current gravity GMM1 of Goddard Space Flight Center.

Description: This report describes the analysis component of the Goddard Earth Observing System, Data Assimilation System, Version 1 (GEOS-1 DAS). The general features of the data assimilation system are outlined, followed by a thorough description of the statistical interpolation algorithm, including specification of error covariances and quality control of observations. We conclude with a discussion of the current status of development of the GEOS data assimilation system. The main components of GEOS-1 DAS are an atmospheric general circulation model and an Optimal Interpolation algorithm. The system is cycled using the Incremental Analysis Update (IAU) technique in which analysis increments are introduced as time independent forcing terms in a forecast model integration. The system is capable of producing dynamically balanced states without the explicit use of initialization, as well as a time-continuous representation of non- observables such as precipitation and radiational fluxes. This version of the data assimilation system was used in the five-year reanalysis project completed in April 1994 by Goddard's Data Assimilation Office (DAO) Data from this reanalysis are available from the Goddard Distributed Active Center (DAAC), which is part of NASA's Earth Observing System Data and Information System (EOSDIS). For information on how to obtain these data sets, contact the Goddard DAAC at (301) 286-3209, EMAIL daac@gsfc.nasa.gov.

Description: Substantial evidence suggests that a UV (ultraviolet) Spectrally Absorbing Material (UV-SAM) exists on Triton's surface. This evidence is found in the positive slope in Triton's spectrum from the UV to the near-IR, and the increasing contrast in Triton's light curve in the blue and UV. Although it is now widely-thought that UV-SAM's exist on Triton, little is known about their distribution and spectral properties. The goal of this NDAP Project is to determine the spatial distribution and geological context of the UV-SAM material. We hope to determine if UV-SAM's on Triton are correlated with geologic wind streaks, craters, calderas, geomorphic/topographic units, regions containing (or lacking) volatile frosts, or some other process (e.g., magnetospheric interactions). Once the location and distribution of UV-SAM's has been determined, further constraints on their composition can be made by analyzing the spectrographic data set. To accomplish these goals, various data sets will be used, including Voyager 2 UV and visible images of Triton's surface, IUE and HST spectra of Triton, and a geologic map of the surface based on Voyager 2 and spectrophotometric data. The results of this research will be published in the planetary science literature.

Description: Consensus on the timing and mapping of substorm features has permitted a synthesis of substorm models. Within the synthesis model the mechanism for onset of substorm expansion is still unknown. Possible mechanisms are: growth of an ion tearing mode, current disruption by a cross-field current instability, and magnetosphere-ionosphere coupling. While the synthesis model is consistent with overall substorm morphology, including near-Earth onset, none of the onset theories, taken individually, appear to account for substorm expansion onset. A grand synthesis with unification of the underlying onset theories appears necessary.

Description: High resolution numerical simulations of thermal convection in a rapidly rotating channel with gravity perpendicular to the rotation vector are described. The convecting columns are subject to a beta-effect resulting from cross-channel topographic vortex stretching. The symmetries of the problem allow many invariant wavenumber sets, and this property is associated with the existence of stable multiple-equilibria at modest supercriticality. The transition to chaotic behavior involves the production of intermittent unstable orbits off a two-torus in energy space. At very high Rayleigh number (of order 10(exp 6) to 10(exp 7)) the motion can be turbulent, depending on the size of beta. However, the turbulence is usually characterized by an almost-periodic formation of patches of small scale convection that cause regular pulsations in the accompanying strong zonal jets. The processes maintaining these flows may be related to those responsible for the zonal currents on Jupiter and for cyclic variability on the Sun.