ALBERT

Description: Sources such as atmospheric or buried explosions and shallow earthquakes producing strong vertical ground displacements produce pressure waves that propagate at infrasonic speeds in the atmosphere. At ionospheric altitudes low frequency acoustic waves are coupled to ionispheric gravity waves and induce variations in the ionoispheric electron density. Global Positioning System (GPS) data recorded in Southern California were used to compute ionospheric electron content time series for several days preceding and following the January 17, 1994, M(sub w) = 6.7 Northridge earthquake. An anomalous signal beginning several minutes after the earthquake with time delays that increase with distance from the epicenter was observed. The signal frequency and phase velocity are consistent with results from numerical models of atmospheric-ionospheric acoustic-gravity waves excited by seismic sources as well as previous electromagnetic sounding results. It is believed that these perturbations are caused by the ionospheric response to the strong ground displacement associated with the Northridge earthquake.

Description: Observations of waves at frequencies below approximately 200 Hz obtained near the magnetopause are presented. For one case identified in the ISEE 1 data as a period when steady state reconnection was occurring, there were waves below the lower hybrid frequency with amplitudes up to approximately 7 mV/m. Intense low-frequency waves with amplitudes up to approximately 20 mV/m at the subsolar magnetopause have also been observed by the Geotail electric field instrument. In some cases, large spiky fields were embedded in the waves. The waves observed by ISEE 1 and Geotail were large enough to provide the dissipation required for reconnection to occur.

Description: Pyranometers have been used for many years to measure broadband surface incoming solar irradiance, data that is necessary for surface energy budget, cloud forcing, and satellite validation research. Because such measurements are made at a specific location, it is unclear how representative they may be of a larger area. This study attempts to determine a reasonable spacing between measurement sites for such research by computing the correlation, and standard deviation from perfect correlation, between simultaneous measurements of incoming solar irradiance for a network of surface measurement sites covering a 75 km x 75 km area. Using 1-min data collected from this network of 11 sites during the NASA First ISSCP Radiation Experiment/Surface Radiation Budget (FIRE/SRB) Project temporal averages were calculated. The correlation between any two of these sites was determined by comparing simultaneous measurement averages for the 55 possible combinations of site pairs, along with the distances between them. In an attempt to remove the effect of the diurnal cycle, thus leaving clouds as the primary influence on correlation of the radiation field, model results for a clear day were used to normalize measured irradiances and correlations were again calculated.

Description: We examine the consequences of magnetic reconnection at the high-latitude magnetopause using a three-dimensional global magnetohydrodynamic simulation of the solar wind interaction with the Earth's magnetosphere. Magnetic field lines from the simulation reveal the formation of magnetic flux ropes during periods with northward interplanetary magnetic field. These flux ropes result from multiple reconnection processes between the lobes field lines and draped magnetosheath field lines that are convected around the flank of the magnetosphere. The flux ropes identified in the simulation are consistent with features observed in the magnetic field measured by Hawkeye-1 during some high-latitude magnetopause crossings.

Description: Analysis of partical data from the CRIT 2 experiment, studying Alfven's critical ionization velocity (CIV) effect, shows that the density of newly created ions (presumably Ba(+) from the shaped-charge beam) is consistent with the increase in total plasma density measured by the independent RF plasma probe on board (Swenson et al., 1990) at the most active time period. We model this ion production using the measured electron flux data and the neutral barium model of Stenbaek-Nielsen et al. (1990a). To identify the main source mechanisms which may contribute most to the barium ionization, a simple model for the barium ion density at the payload location is developed based on Liouvilles theorem. We estimate that the electron impact ionization is responsible for 90% of the barium ion production observed by CRIT 2 in the first release and up to 45% in the second release. By employing a two-state approximation calculation (Rapp and Francis, 1962), the Ba-O(+) charge exchange cross section is found to range from about 2.0 X 10(exp -17) sq cm at a velocity of 4 km/s to 2.0 X 10(exp -15) sq cm at a velocity of 20 km/s. This result suggests that the Ba-O(+) charge exchange is probably dominant among all the non-CIV ionization processes. By considering the charge exchange process in our density model, the barrium ion densities are calculated for the two releases on CRIT II. The comparison between the model results and the observed data is found to be resonably consistent if the cross sections, as calculated above, are multiplied by 0.3 for the first release and 1.0 for the second release. Our result suggests that the charge exchange process could be the most important non-CIV ionization mechanism in the CRIT II experiment and it should be considered carefully case by case in CIV experiments.

Description: A selected set of far ultra violet images of Earth have been analyzed quanitatively to establish their validity for studying thermospheric weather. The set of images chosen for the study was restricted to mostly geomagnetically quiet conditions in order to obtain a baseline understanding of the relationship between the observations and thermospheric phenomenology. The images included low to modrerate solar activity levels. A new model was developed to generate global dayglow images using first principles methods. The mass Spectrometer/incoherent scatter (MSIS-86) model was used to predict the thermospheric concentrations. The analyses of thermospheric images observed in the 123 to 160-nm nominal passband show that the spectral composition for observations on the projected earth disk is dominated by O I 130.4-nm radiation (85-90%), with concentrations from O I 135.6-nm and N2 Lyman-Birge-Hopefield (LBH) bands of about 5-8% each. The synthetic images reproduce the global features of the observed images rather well. Differences between the model and the data are attributed to real atmospheric effects, such as atomic oxygen depletions which are not well reproduced by the MSIS model when geomagnetic activity is elevated. The absoute values recorded were 38-54% higher than predicted. We attribute this discrepancy to low values of the solar extreme ultraviolet irradiances used in the model. Images obtained in the 136 to 165-nm nominal passband are a factor of 2.7 greater than the model. The excess signal observed is most likely due to a long wavelength tail in the instrument sensitivity which allowed Rayleigh scattered sunlight between 180and 250nm to be detected. The understanding of the Dynamics Explorer (DE 1) images gained by this study provides the basis for future work on the global response of the theremosphere to geomagnetic forcing.

Description: There are four low-frequency modes which may propagate in a high-beta nearly bi-Maxwellian plasma. These are the magnetosonic, Alfven, ion acoustic, and mirror modes. This manuscript defines a procedure based on linear Vlasov theory for the unique identification of these modes by use of transport ratios, dimensionless ratios of the fluctuating field and plasma quantities. A single parameter, the mode deviation is calculated using the plasma and magnetic field data gathered by the Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM) spacecraft to identify the modes observed in the terrestial magnetosheath near the magnetopause. As well as determining the mode which best describes the observed fluctuations, it gives us a measure of whether or not the resulting identification is unique. Using 17 time periods temporally close to a magnetopause crossing, and confining our study to the frequency range from 0.01 to 0.04 Hz, we find that the only clearly identified mode in this frequency range is the mirror mode. Most commonly, the quasi-perpendicular mirror mode (with wave vector k roughly perpendicular to the background magnetic field B(sub zero) is observed. In two events the quasi-parallel mirror mode k parallel B(sub zero) was identified.

Description: Low-altitude spacecraft on magnetotail field lines often detect a distinctive signature in the precipitating ion flux. A velocity-dispersed ion structure is often observed near the poleward boundary of the auroral oval. At the low-latitude edge of this structure an absence of precipitating ions is seen, previously referred to as 'the gap,' separating the velocity -- dispersed ions at the higher latitudes from the more diffuse, plasma sheet-like ions at lower latitudes. We present a model of low-altitude particle precipitation that reproduces these observed features in the ion spectra and provides a quantitative estimate of the downtail plasma sheet properties. The model calculations are compared with observations from the Akebono spacecraft. In this model, the dispersed ion velocity signature maps to a region in the distant plasma sheet where the plasma has a field-aligned bulk flow. The gap maps to a region in the distant magnetotail where the ion fluxes are below the detection threshold of the instrument, due to the low plasma sheet density and temperature in that region.

Description: Earthquake recurrence data from the Pallett Creek and Wrightwood paleoseismic sites on the San Andreas fault appear to show temporal variations in repeat interval. We investigate the interaction between strike-slip faults and auxiliary reverse and normal faults as a physical mechanism capable of producing such variations. Under the assumption that fault strength is a function of fault-normal stress (e.g. Byerlee's Law), failure of an auxiliary fault modifies the strength of the strike-slip fault, thereby modulating the recurrence interval for earthquakes. In our finite element model, auxiliary faults are driven by stress accumulation near restraining and releasing bends of a strike-slip fault. Earthquakes occur when fault strength is exceeded and are incorporated as a stress drop which is dependent on fault-normal stress. The model is driven by a velocity boundary condition over many earthquake cycles. Resulting synthetic strike-slip earthquake recurrence data display temporal variations similar to observed paleoseismic data within time windows surrounding auxiliary fault failures. Our simple model supports the idea that interaction between a strike-slip fault and auxiliary reverse or normal faults can modulate the recurrence interval of events on the strike-slip fault, possibly producing short term variations in earthquake recurrence interval.

Description: Global Positioning System (GPS) measurements spanning approximately 3 years have been used to determine velocities for 7 sites on the Australian, Pacific and Antarctic plates. The site velocities agree with both plate model predictions and other space geodetic techniques. We find no evidence for internal deformation of the interior of the Australian plate. Wellington, New Zealand, located in the Australian-Pacific plate boundary zone, moves 20 +/- 5 mm/yr west-southwest relative to the Australian plate. Its velocity lies midway between the predicted velocities of the two plates. Relative Euler vectors for the Australia-Antarctica and Pacific-Antarctica plates agree within one standard deviation with the NUVEL-1A predictions.

Description: Coseismic surface deformation associated with the M(sub w) 6.1, April 23, 1992, Joshua Tree earthquake is well represented by estimates of geodetic monument displacements at 20 locations independently derived from Global Positioning System and trilateration measurements. The rms signal to noise ratio for these inferred displacements is 1.8 with near-fault displacement estimates exceeding 40 mm. In order to determine the long-wavelength distribution of slip over the plane of rupture, a Tikhonov regularization operator is applied to these estimates which minimizes stress variability subject to purely right-lateral slip and zero surface slip constraints. The resulting slip distribution yields a geodetic moment estimate of 1.7 x 10(exp 18) N m with corresponding maximum slip around 0.8 m and compares well with independent and complementary information including seismic moment and source time function estimates and main shock and aftershock locations. From empirical Green's functions analyses, a rupture duration of 5 s is obtained which implies a rupture radius of 6-8 km. Most of the inferred slip lies to the north of the hypocenter, consistent with northward rupture propagation. Stress drop estimates are in the range of 2-4 MPa. In addition, predicted Coulomb stress increases correlate remarkably well with the distribution of aftershock hypocenters; most of the aftershocks occur in areas for which the mainshock rupture produced stress increases larger than about 0.1 MPa. In contrast, predicted stress changes are near zero at the hypocenter of the M(sub w) 7.3, June 28, 1992, Landers earthquake which nucleated about 20 km beyond the northernmost edge of the Joshua Tree rupture. Based on aftershock migrations and the predicted static stress field, we speculate that redistribution of Joshua Tree-induced stress perturbations played a role in the spatio-temporal development of the earth sequence culminating in the Landers event.

Description: Ions that are energized at quasi-parallel collisionless shocks and move back upstream generate low-frequency waves, largely on the fast/magnetosonic branch. At sufficient Mach number, the waves are convected back into the shock, lead to shock re-formation, and are mode-converted into downstream (magnetosheath) Alfvenic turbulence. Other waves are generated more locally at the interface of the incoming solar wind and the partially thermalized plasma. This paper reviews how recent simulation studies of collisionless shocks in conjunction with linear kinetic theory and proper wave diagnostics have aided in our understanding of the upstream and magnetosheath waves.

Description: A variety of suprathermal and energetic ion distributions are found upstream from shocks. Some distributions, such as field-aligned beams, are generated directly at the shock either through reflection processes or through leakage from the hotter downstream region. Other distributions, such as intermediate distributions, evolve from these parent distributions through wave-particle interactions. This paper reviews our current understanding of the creation and evolution of suprathermal distributions at shocks. Examples of suprathermal ion distributions are taken from observations at the Earth's bow shock. Particular emphasis is placed on the creation of field-aligned beams and specularly reflected ion distributions and on the evolution of these distributions in the Earth's ion foreshock. However, the results from this heavily studied region are applicable to interplanetary shocks, bow shocks at other planets, and comets.

Description: In this study, Atmosphere Explorer data and model results for the ion and electron temperature and the density of N(+), O(+), H(+), and He(+) between 120 and 1400 km altitude are compared for two midlatitude ranges (L=2 and L=4), noon and midnight local time, winter and summer, at solar minimum. The data for the heavy atomic ions (O(+) and N(+)) show that their densities are greater at noon than at midnight for a given season and greater in summer than winter for a given local time. There is only a weak latitudinal variation in the density of these ions. The data show that the light ion (H(+) and He(+)) densities are greater at midnight than at noon and are generally greater in winter than summer. There is a strong latitudinal variation of the light ion densities, with the densities decreasing with increasing latitude. The model densities are in good agreement with the AE densities for N(+), O(+), and H(+). Model He(+) densities are lower, by a factor of 2 or more, than the measured densities. Model ion and electron temperatures agree well with the measured temperatures with only a modest increase in plasmapheric heating.

Description: From a study of the 21 largest geomagnetic storms during solar cycle 21, a strong correlation is established between the ring current index Dst and the time-weighted accumulation of the 1-hour auroral electrojets indices, AE and AL. The time-weighted accumulation corresponds to convolution of the auroral electrojet indices with an exponential weighting function with an e-folding time of 9.4 hours. The weighted indices AE(sub w) and AL(sub w) have correltation coefficients against Dst ranging between 0.8 and 0.95 for 20 of the 21 storms. Correlation over the entire solar cycle 21 database is also strong but not as strong as for an individual storm. A set of simple Dst prediction functions provide a first approximation of the inferred dependence, but the specific functional relationship of Dst (AL(sub w)) or Dst (AL(sub w)) varies from one storm to the next in a systematic way. This variation reveals a missing parametric dependence in the transfer function. However, our results indicate that auroral electroject indices are potentially useful for predicting storm time enhancements of ring current intensity with a few hours lead time.

Description: The poloidal mode field line resonance in the Earth's dipole magnetic field is investigated using cold plasma ideal MHD simulations in dipole geometry. In order to excite the poloidal mode resonance, we use either an initial or a continuous velocity perturbation to drive the system. The perturbation is localized at magnetic shell L = 7 with plasma flow in the radial direction (electric field component in the azimuthal direction). It is found that with the initial perturbation alone, no polodial mode resonance can be obtained and the initially localized perturbation spreads out across all magnetic L shells. With the continuous perturbation, oscillating near the poloidal resonance frequency, a global-scale poloidal cavity mode can be obtained. For the first time, a localized guided poloidal mode resonance is obtained when a radial component of electric field is added to the initial perturbation such that the curl of the electric field is everywhere perpendicular to the background dipole magnetic field. During the localized poloidal resonance, plasma vortices parallel/antiparallel to the background dipole magnetic field B(sub 0). This circular flow, elongated radially, results in twisting of magnetic field flux tubes, which, in turn, leads to the slowdown of the circular plasma flow and reversal of the plasma vortices. The energy associated with the localized poloidal resonance is conserved as it shifts back and forth between the oscillating plasma vortices and the alternately twisted magnetic flux tubes. In the simulations the eigenfunctions associated with the localized poloidal resonance are grid-scale singular functions. This result indicates that ideal MHD is inadequate to describe the underlying problem and nonideal MHD effects are needed for mode broadening.

Description: The long-term time series of global ozone from the Nimbus-7 Solar Backscatter Ultraviolet instrument (SBUV) (Nov. 1978-June 1990) are extended through June 1994 by using measurements from the NOAA-11 SBUV/2. The data sets are merged based upon comparisons during the 18-month overlap period in which both instruments were operational. During this period, the average offset between the two time series is less than 2% in total ozone, and less than 6% in Umkehr layers 1-10. A linear-regression trend model is applied to the extended time series to calculate updated trends as a function of latitude and altitude. Trends through June 1994 are 1.5-2% per decade less negative than through June 1990 in the tropical middle stratosphere (35-40 km) and in the upper stratosphere (45-50 km) at mid-latitudes. In the lower stratosphere, the trends are nearly 1.5% per decade more negative in the southern hemisphere tropical regions to 25 deg S, but remain relatively unchanged elsewhere. The seasonal structure of the total ozone trends is similar to past trend study results, but the magnitude of the seasonal trend can vary by 2% per decade depending on the length of the time series. Both Total Ozone Mapping Spectrometer (TOMS) (through April 1993) and SBUV total ozone time series show small negative trends in the equatorial region, though they are not statistically at the 2-sigma level.

Description: A new technique for estimating electron density perturbation amplitudes of traveling ionospheric disturbances (TIDs), using HF radar data, is presented. TIDs are observed in HF radar data as enhancements of the ground-scattered power which propagate through the radar's field of view. These TIDs are the ionospheric manifestation of atmospheric acoustic-gravity waves. TID electron density perturbation amplitudes were estimated by simulating the radar returns, using HF ray tracing through a model ionosphere perturbed by a model gravity wave. The simulation determined the return power in the ground-scattered portion of the signal as a function of range, and this was compared to HF radar data from the Goose Bay HF radar at a time when evidence of gravity waves was present in the data. By varying the amplitude of the electron density perturbation in the model it was possible to estimate the perturbation of the actual wave. It was found that the perturbations that are observed by the Goose Bay HF radar are of the order of 20% to 35%. It was also found that the number of observable power enhancements, and the relative amplitudes of these enhancements, depended on the vertical thickness of the gravity wave's source region. From the simulations and observations it was estimated that the source region for the case presented here was approximately 20 km thick. In addition, the energy in the wave packet was calculated and compared to an estimate of the available energy in the source region. It was found that the wave energy was about 0.2% of the estimated available source region energy.

Description: Distribution functions of ions heated in quasi-perpendicular bow shocks have a large perpendicular temperature anisotropy that provides free energy for the growth of Alfven ion cyclotron (AIC) waves and mirror waves. Both types of waves have been observed in the Earth's magnetosheath downstream of quasi-perpendicular shocks. We use a two-dimensional hybrid simulations to give a self-consistent description of the evolution of the wave spectra downstream of quasi-perpendicular shocks. Both mirror and AIC waves are identified in the simulated magnetosheath. They are generated at or near the shock front and convected away from it by the sheath plasma. Near the shock, the waves have a broad spectrum, but downstream of the shock, shorter-wavelength modes are heavily damped and only longer-wavelength modes persist. The characteristics of these surviving modes can be predicted with reasonable accuracy by linear kinetic theory appropriate for downstream conditions. We also follow the evolution of the ion distribution function. The shocked ions that provide the free energy for wave growth have a two-component distribution function. The halo is initially gyrophase-bunched and extremely anisotropic. Within a relatively short distance downstream of the shock (of the order of 10 ion inertial lengths), wave-particle interactions remove these features from the halo and reduce the anisotropy of the distribution to near-threshold levels for the mirror and AIC instabilities. A similar evolution has been observed for ions at the Earth's bow shock.

Description: The following areas of interest in this progress report are: (1) the continuation of software development in the examination of F-region gravity-wave power using in-situ data from the Atmosphere Explorer (AE-E); (2) the inquiry into the use of the San Marco data for the study of the initiation and growth of bubbles, particularly when the satellite passes through the early evening hours at relatively high altitudes, and the development of bubbles using not only the San Marco data but includes the use of airglow observations made in Hawaii; and (3) the promising development in the observation of distinct well formed waves at about 400 km altitude in the equatorial region. These waves look very much like waves seen over the polar cap that are attributed to internal gravity waves in the neutral atmosphere driving ionization up and down the magnetic field lines. These equatorial waves show no modulation of the total ion concentration.

Description: This final report was concerned with the ideas that: (1) magnetospheric boundary layers link disparate regions of the magnetosphere-solar wind system together; and (2) global behavior of the magnetosphere can be understood only by understanding its internal linking mechanisms and those with the solar wind. The research project involved simultaneous research on the global-, meso-, and micro-scale physics of the magnetosphere and its boundary layers, which included the bow shock, the magnetosheath, the plasma sheet boundary layer, and the ionosphere. Analytic, numerical, and simulation projects were performed on these subjects, as well as comparisons of theoretical results with observational data. Other related activity included in the research included: (1) prediction of geomagnetic activity; (2) global MHD (magnetohydrodynamic) simulations; (3) Alfven resonance heating; and (4) Critical Ionization Velocity (CIV) effect. In the appendixes are list of personnel involved, list of papers published; and reprints or photocopies of papers produced for this report.

Description: A brief description of the major activities pursued during the last year (March 1994 - February 1995) of this grant are: (1) the development of a 200 km to 1 Re, O(+) H(+) Model; (2) the extension of the E x B convection heating study to include centrifugal effects; (3) the study of electron precipitation effects; (4) the study of wave heating of O(+); and (5) the polar wind acceleration study. A list of both papers published and papers submitted, along with a proposal for next year's study and a copy of the published paper is included.

Description: The Langmuir probe flown as part of the Solar Array Module Plasma Interactions Experiment (SAMPIE) package aboard the space shuttle flight STS-62 was used to determine plasma potential fluctuations in the vicinity of the shuttle. The broadband noise was observed at frequencies 250 - 20,000 Hz. Measurements were performed in ram conditions; thus, it seems reasonable to believe that the influence of spacecraft operations on plasma parameters was absolutely negligible. The average spectrum of fluctuations is in agreement with theoretical predictions. The influence on the observed spectra of arcing generated by high negative bias voltages applied to solar cell samples is briefly discussed.

Description: The SBUV/2 instrument onboard the NOAA-11 satellite made daily solar spectral irradiance measurements in the wavelength region 160405 nm at 1.1 nm resolution between January 1989 and October 1994. These observations continued the uninterrupted series of solar measurements begun by the Nimbus-7 SBUV in 1978 and continued by NOAA-9 SBUV/2. While the measurements made by the SBUV-series instruments furnish an excellent data base for studies of solar UV variability, these instruments do not have an internal mew to evaluate and correct for long-term instrument sensitivity degradation, needed to evaluate solar cycle timescale irradiance change. During yearly Shuttle flights the Shuttle SBUV (SSBUV) also performs solar spectral irradiance measurements in the wavelength region 200 to 400 nm with an instrument that is calibrated preflight, inflight, and postflight. Comparisons between the simultaneous NOAA-11 SBUV/2 and SSBUV solar measurements are used to identify and correct long term sensitivity changes in the satellite instrument. The NOAA-11 data will then be used to evaluate long-term solar change. We present a progress report on the above process. At this preliminary stage uncertainties in the calibration transfer between SSBUV and NOAA-11 SBUV/2 are too large to accurately evaluate long-term solar change near the A1 edge, but solar rotational activity variations can be evaluated. We find that rotational activity declined from roughly 6% peak-to-peak (p-p) near the maximum of solar cycle 22 in 1989-1991 to approximately 3% p-p in mid 1992 and 2% p-p by mid 1994. Emphasizing rotational variations, comparisons between the 200 nm data and the NOAA-11 Mg II proxy index are presented.

Description: A far-wing line shape theory based on the binary collision and quasistatic approximations that is applicable for both the low- and high-frequency wings of the vibration-rotational bands has been developed. This theory has been applied in order to calculate the frequency and temperature dependence of the continuous absorption coefficient for frequencies up to 10,000 cm(exp -1) for pure H2O and for H2O-N2 mixtures. The calculations were made assuming an interaction potential consisting of an isotropic Lennard-Jones part with two parameters that are consistent with values obtained from other data, and the leading long-range anisotropic part, together with the measured line strengths and transition frequencies. The results, obtained without the introduction of adjustable parameters, compare well with the existing laboratory data, both in magnitude and in temperature dependence. This leads us to the conclusion that the water continuum can be explained in terms of far-wing absorption. Current work in progress to extend the theory and to validate the theoretically calculated continuum will be discussed briefly.

Description: The Argon Release for Controlled Studies (ARCS) 4 sounding rocket was launched northward into high altitude from Poker Flat Research Range on February 23, 1990. The vehicle crossed geomagnetic field lines containing discrete auroral activity. An instrumented subpayload released 100-eV and 200-eV Ar(+) ion beams sequentially, in a direction largely perpendicular to both the local geomagnetic field and the subpayload spin axis. The instrumented main payload was separated along field lines from the beam emitting subpayload by a distance which increased at a steady rate of approximately 2.4 m/s. Three dimensional mass spectrometric ion observations of ambient H(+) and O(+) ions, obtained on board the main payload, are presented. Main payload electric field observations in the frequency range 0-16 kHz, are also presented. These observations are presented to demonstrate the operation of transverse ion acceleration, which was differential with respect to ion mass, primarily during 100-eV beam operations. The preferential transverse acceleration of ambient H(+) ions, as compared with ambient O(+) ions, during the second, third, fourth, and fifth 100-eV beam operations, is attributed to a resonance among the injected Ar(+) ions, beam-generated lower hybrid waves, and H(+) ions in the tail of the ambient thermal distribution. This work provides experimental support of processes predicted by previously published theory and simulations.

Description: Spectral and spatial images obtained with the Imaging Spectrometric Observatory on the ATLAS 1 and Spacelab 1 missions are used to study the ultraviolet emissions of nitric oxide in the thermosphere. By synthetically fitting the measured NO gamma bands, intensities are derived as a function of altitude and latitude. We find that the NO concentrations inferred from the ATLAS 1 measurements are higher than predicted by our thermospheric airglow model and tend to lie to the high side of a number of earlier measurements. By comparison with synthetic spectral fits, the shape of the NO gamma bands is used to derive temperature as a function of altitude. Using the simultaneous spectral and spatial imaging capability of the instrument, we present the first simultaneously acquired altitude images of NO gamma band temperature and intensity in the thermosphere. The lower thermospheric temperature images show structure as a function of altitude. The spatial imaging technique appears to be a viable means of obtaining temperatures in the middle and lower thermosphere, provided that good information is also obtained at the higher altitudes, as the contribution of the overlying, hotter NO is nonnegligible. By fitting both self-absorbed and nonabsorbed bands of the NO gamma system, we show that the self absorption effects are observable up to 200 km, although small above 150 km. The spectral resolution of the instrument (1.6 A) allows separation of the N(+)(S-5) doublet, and we show the contribution of this feature to the combination of the NO gamma (1, 0) band and the N(+)(S-5) doublet as a function of altitude (less than 10% below 200 km). Spectral images including the NO delta bands support previous findings that the fluorescence efficiency is much higher than that determined from laboratory measurements. The Spacelab 1 data indicate the presence of a significant population of hot NO in the vehicle environment of that early shuttle mission.

Description: This work describes a method for estimating vertical fluxes of horizontal momentum carried by short horizontal scale gravity waves (lambda(sub x) = 10-100 km) using aircraft measured winds in the lower stratosphere. We utilize in situ wind vector and pressure altitude measurements provided by the Meteorological Measurement System (MMS) on board the ER-2 aircraft to compute the momentum flux vectors at the flight level above deep convection during the tropical experiment of the Stratosphere Troposphere Exchange Project (STEP-Tropical). Data from Flight 9 are presented here for illustration. The vertical flux of horizontal momentum these observations points in opposite directions on either side of the location of a strong convective updraft in the cloud shield. This property of internal gravity waves propagating from a central source compares favorably with previously described model results.

Description: The major activities of the Global Tropospheric Experiment at the Hong Kong Atmospheric Chemistry Measurement Station are presented for the period 1 January - 31 December 1995. Activities included data analysis, reduction, and archiving of atmospheric measurements and sampling. Sampling included O3, CO, SO2, NO, TSP, RSP, and ozone column density. A data archive was created for the surface meteorological data. Exploratory data analysis was performed, including examination of time series, frequency distributions, diurnal variations and correlation. The major results have been or will be published in scientific journals as well as presented at conferences/workshops. Abstracts are attached.

Description: The work performed under this contract provides for the improvement of modeling of troposphere propagation delay in very long base interferometry (VLBI) data, defining of required accuracy of ancillary data type(s) to measure tropospheric behavior and to develop a clear correspondence between tropospheric behavior and the quality of geodetic VLBI measurements.

Description: The main focus of this work is to understand the dynamics of mass exchange between the tropics and the midlatitudes and to determine any links between tropospheric exchange and that in the stratosphere. We have approached this problem from two different perspectives. The first is aimed towards understanding the troposphere's role in inducing lower stratospheric tropical/midlatitude exchange. For this project we focus on observational analysis of the lower stratosphere to assess the key regions of transport in/out of the tropics and to what extent this transport is driven by tropospheric processes. The second approach is to determine the extent to which stratospheric processes influence the troposphere. In this project we are performing potential vorticity (PV) inversions to assess the winds induced near the tropopause when the stratospheric polar vortex is displaced equatorward. These are each discussed in more detail in the subsections below. Also, we have organized a session on Tropical/Midlatitude Interaction and Transport at the Fall AGU where we will be showing our latest results.

Description: The major focus of the subject contract was on helping to resolve one of the more notable discrepancies still existing in the axial momentum budget of the solid Earth-atmosphere system, namely the disappearance of coherence between length-of-day (l.o.d.) and atmospheric angular momentum (AAM) at periods shorter than about a fortnight. Recognizing the importance of identifying the source of the high-frequency momentum budget anomaly, the scientific community organized two special measurement campaigns (SEARCH '92 and CONT '94) to obtain the best possible determinations of l.o.d. and AAM. An additional goal was to analyze newly developed estimates of the torques that transfer momentum between the atmosphere and its underlying surface to determine whether the ocean might be a reservoir of momentum on short time scales. Discrepancies between AAM and l.o.d. at sub-fortnightly periods have been attributed to either measurement errors in these quantities or the need to incorporate oceanic angular momentum into the planetary budget. Results from the SEARCH '92 and CONT '94 campaigns suggest that when special attention is paid to the quality of the measurements, better agreement between l.o.d. and AAM at high frequencies can be obtained. The mechanism most responsible for the high-frequency changes observed in AAM during these campaigns involves a direct coupling to the solid Earth, i.e, the mountain torque, thereby obviating a significant oceanic role.

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: The NASA Combined Release and Radiation Effects Satellite (CRRES) El Coqui rocket campaign was successfully carried out in Puerto Rico during the period 18 May through 12 July 1992. This report describes five chemical release experiments in the upper ionosphere supported by Geospace Research, Inc. during the El Coqui campaign. Additional spin-off science is also discussed. The El Coqui releases are designated AA-1 (rocket 36-082), AA-2 (rocket 36-081), AA-3b (rocket 36-064), AA-4 (rocket 36-065), and AA-7 (rocket 36-083). Particular attention is paid to releases AA-2 and AA-4. These two experiments involved the illumination of ionospheric release regions with powerful high-frequency (HF) radio waves transmitted from the Arecibo HF facility. In the AA-2 experiment, microinstabilities excited by the HF wave in a Ba(+) plasma were examined. This release yielded a smooth plasma cloud that helped clarify several fundamental issues regarding the physics of wave plasma instabilities. During AA-2 extremely strong HF-induced Langmuir turbulence was detected with the Arecibo 430 MHz radar. CF3Br was released in the AA-4 study to create an ionospheric hole that focused the HF beam. This experiment successfully explored wave-plasma coupling in an O(+) ionosphere under conditions of very high HF electric field strengths.

Description: The primary objective of the three-day workshop on results from the Data Assimilation Office (DAO) five-year assimilation was to provide timely feedback from the data users concerning the strengths and weaknesses of version 1 of the Goddard Earth Observing System (GEOS-1) assimilated products. A second objective was to assess user satisfaction with the current methods of data access and retrieval. There were a total of 49 presentations, with about half (23) of the presentations from scientists from outside of Goddard. The first two days were devoted to applications of data: studies of the energy diagnostics, precipitation and diabatic heating, hydrological modeling and moisture transport, cloud forcing and validation, various aspects of intraseasonal, seasonal, and interannual variability, ocean wind stress applications, and validation of surface fluxes. The last day included talks from the National Meteorological Center (NMC), the National Center for Atmospheric Research (NCAR), the Center for Ocean-Land-Atmosphere Studies (COLA), the United States Navy, and the European Center for Medium Range Weather Forecasts (ECMWF).

Description: The fast neutron flux in near-Earth orbit has been measured with the COMPTEL instrument on the Compton Gamma Ray Observatory (CGRO). For this measurement one of COMPTEL's seven liquid scintillator modules was used as an uncollimated neutron detector with threshold of 12.8 MeV. The measurements cover a range of 4.8 to 15.5 GV in vertical cutoff rigidity and 3 deg to 177 deg in spacecraft geocenter zenith angle. One of the measurements occurred near the minimum of the deepest Forbush decrease ever observed by ground-level neutron monitors. After correction for solar modulation, the total flux is well fitted by separable functions in rigidity and zenith angle. With the spacecraft pointed near the nadir the flux is consistent with balloon measurements of the atmospheric neutron albedo. The flux varies by about a factor of 4 between the extremes of rigidity and a factor of 2 between the extremes of zenith angle. The effect of the spacecraft mass in shielding the detector from the atmospheric neutron albedo is much more important than its role as a source of additional secondary neutrons. The neutron spectral hardness varies little with rigidity or zenith angle and lies in the range spanned by earlier atmospheric neutron albedo measurements.

Description: Forty-three examples of ISEE 1 tailward flank side magnetopause crossings are examined and directly compared with upstream solar wind parameters. The crossings are classified into two groups. In the first group, a few sudden magnetopause crossings are observed, whereas repeated magnetopause crossings and oscillatory motions, often with boundary layer signatures, are observed in the second group. These distinctive characteristics of the two groups are interpreted in terms of the surface waves due to the Kelvin-Helmholtz instability. It is found that low solar wind speed tends to favor characteristics of the first group, whereas high solar wind speed yields those of the second group. However, no evident correlations between the groups and the interplanetary magnetic field directions are found.

Description: During periods of southward interplanetary magnetic field (IMF) orientation the magnetic field geometry at the dayside magnetopause is susceptible to magnetic reconnection. It has been suggested that reconnection may occur in a localized manner at several patches on the magnetopause. A major problem with this picture is the interaction of magnetic flux ropes which are generated by different reconnection processes. An individual flux rope is bent elbowlike where it intersects the magnetopause and the magnetic field changes from magnetospheric to interplanetary magnetic field orientation. Multiple patches of reconnection can lead to the formation of interlinked magnetic flux tubes. Although the corresponding flux is connected to the IMF the northward and southward connected branches are hooked into each other and cannot develop independently. We have studied this problem in the framework of three-dimensional magnetohydrodynamic simulations. The results indicate that a singular current sheet forms at the interface of two interlinked flux tubes if no resistivity is present in the simulation. This current sheet is strongly tilted compared to the original current sheet. In the presence of resistivity the interaction of the two flux tubes forces a fast reconnection process which generates helically twisted closed magnetospheric flux. This linkage induced reconnection generates a boundary layer with layers of open and closed magnetospheric flux and may account for the brightening of auroral arcs poleward of the boundary between open and closed magnetic flux.

Description: Classical models of magnetic reconnection consist of a small diffusion region bounded by two symmetric slow shocks, across which the plasma is accelerated. Asymmetries often present in space plasmas are sheared plasma flow and dissimilar plasma densities on the two sides of current sheets. In this paper, we investigate magnetic reconnection in the presence of a shear flow and an asymmetric density across the current sheet using two-dimensional magnetohydrodynamic (MHD) simulations. The results demonstrate that magnetic reconnection can occur only for a plasma flow velocity (in the frame of the X line) which is below the Alfven speed in each inflow region. This limits the velocity of the X line to a certain range for a given flow shear and provides an upper limit to the total velocity shear at which reconnection ceases to operate. Depending on the direction of the flow in the adjacent inflow region, the effects from the sheared flow and from the density asymmetry will compete with or enhance each other in respect to the magnitude and location of the currents which bound the outflow regions. The results are applied to the dayside and flank regions of the magnetosphere. For the dayside region where the magnetosheath flow is slow, the magnetic field transition region is thin and the accelerated flow is earthward of the sharp current layer (magnetopause). At the flanks tailward of the X line, shear flow and density asymmetry effects compete making the magnetic field transition layer broad with the high-speed flow contained within the transition region which explains corresponding observations. At the flanks sunward of the X line, shear flow and density asymmetry effects enhance each other and lead to a strong current sheet on the magnetosheath side of the accelerated flow. The total volume affected by magnetic reconnection is much larger than the steady state region. A large bulge region precedes the steady state region. Qualitatively, the bulge and the steady state region have similar signatures and both can explain observations. We provide criteria in order to distinguish between the bulge and the steady state region in observations.

Description: A magnetofricional method is used to construct two-dimensional MHD equilibria of the Earth's magnetosphere for a given distribution of entropy functions(S = pV(exp gamma), where p is the plasma pressure and V is the tube volume per unit magnetic flux. It is found that a very thin current sheet with B (sub zeta) is less than 0.5 nu T and thickness less than 1000 km can be formed in the near-earth magnetotail (x is approximately -8 to -20R(sub e) during the growth phase of substorm. The tail current sheets are found to become thinner as the entropy or the entropy gradient increases. It is suggested that the new entropy anti-diffusion instability associated with plasma transport across field lines leads to magnetic field dipolarization and accelerates the formation of thin current sheet, which may explain the observed explosive growth phase of substorms.

Description: Strong interplanetary magnetic field (IMF) and seasonal effects in the convection of nightside ionospheric plasma are described. The findings are based on a statistical analysis of observations made with the Johns Hopkins University/ Applied Physics Lab (JHU/APL) HF radar located at Goose Bay, Labrador. For positive sign of the IMF dusk-dawn component, By greater than 0 the dawn cell is more crescent shaped and the dusk cell more round while for BY less than 0 these pairings of size and shape are reversed. The more extreme crescent /round cell dichotomy is obtained for BY greater than 0. The return flows associated with the crescent-shaped cell dominate at midnight MLT (magnetic local time); the reversal in the zonal velocity in the 67 deg-69 deg lambda (magnetic latitude) interval occurs 2.5 hr earlier in summer than in winter. The maximum effects are obtained on the nightside for the pairings By greater than 0, summer and BY less than 0, winter; the first produces the more structured cell in the morning, the second in the evening, and this cell dominates the return flow at midnight. The difference in the zonal flow reversals for these pairings exceeds 4 hr in MLT.

Description: Impulsive ELF/VLF electric field bursts observed by the vector electric field instrument (VEFI) on the Dynamics Explorer 2 (DE 2) satellite on almost every crossing of the geomagnetic equator in the evening hours are interpreted as originating in lightning discharges. These signals that peak in intensity near the magnetic equator are observed within 5-20 deg latitude of the geomagnetic equator at altitudes of 300-500 km with amplitudes of the order of approximately mV/m in the 512- or 1024-Hz frequency band of the VEFI instrument. Whistler-mode ELF/VLF wave propagation through a horizontally stratified ionosphere predicts strong attenuation of subionospheric signals reaching the equator at low altitudes. However, ray tracing analysis shows that the presence of the equatorial density anomaly, commonly observed in the upper ionosphere during evening hours, leads to the focusing of the wave energy from lightning near the geomagnetic equator at low altitudes, thus accounting for all observed aspects of the phenomenon. The observations presented here indicate that during certain hours in the evening, almost all the energy input from lightning discharges entering the ionosphere at less than 30 deg latitude remains confined to a small region (in altitude and latitude) near the geomagnetic equator. The net wideband electric field, extrapolated from the observed electric field values in the 512- to 1024-Hz band, can be approximately 10 mV/m or higher. These strong electric fields generated in the ionosphere by lightning at local evening times may be important for the equatorial electrodynamics of the ionosphere.

Description: The tail currents predicted by empirical magnetic field models and global MHD simulations are compared. It is shown that the near-Earth currents obtained from the MHD simulations are much weaker than the currents predicted by the Tsyganenko models, primarily because the ring current is not properly represented in the simulations. On the other hand, in the mid-tail and distant tail the lobe field strength predicted by the simulations is comparable to what is observed at about 50 R(sub E) distance, significantly larger than the very low lobe field values predicted by the Tsyganenko models at that distance. Ways to improve these complementary approaches to model the actual magnetospheric configuration are discussed.

Description: A new concept for an active experiments mission is presented: The Active Magnetospheric Particle Acceleration Satellite (AMPAS). It is proposed to fly a dual-payload tethered satellite to study three important areas of magnetopsheric physics: (1) electric and magnetic field structures of the Earth's magnetosphere, (2) electron beam-plasma interactions, and (3) far ultraviolet (FUV) signatures of energetic electron precipitation. The tethered system is to be flown in a gravity-gradient stabilized configuration with a 70 deg circular orbit at altitudes in the range 300-800 km. The upper payload will carry two electron beam accelerators (1-10 keV, 1 A) and a plasma contractor to electrically neutralize the payload, while the lower will carry a complement of optical imaging and plasma diagonstics instrumentation. Upward directed beams are injected to 'sound' electric and magnetic field structures, including parallel electric fields in the auroral region, while simultaneous optical instruments aimed downwards measure characteristics of reflected beam pulses precipitating into the upper atmosphere below the satellite. Downward directed beams, and simultaneous plasma wave, particle and optical measurements, are used for the investigation of beam-plasma interactions and electron precipitation studies. The payloads are to be connected by a 1-6 km long tether to be deployed in stages. The configuration of the beam source and the plasma and optical diagnostics on separate, tethered payloads allows the diagnostics to be performed outside of the disturbed region around the source payload, and plasma wave and particle observations of beam-plasma interactions to be performed at well-defined locations relative to the source.

Description: As part of the SAMPIE (The Solar Array Module Plasma Interaction Experiment) program, the Langmuir probe (LP) was employed to measure plasma characteristics during the flight STS-62. The whole set of data could be divided into two parts: (1) low frequency sweeps to determine voltage-current characteristics and to find electron temperature and number density; (2) high frequency turbulence (HFT dwells) data caused by electromagnetic noise around the shuttle. The broadband noise was observed at frequencies 250-20,000 Hz. Measurements were performed in ram conditions; thus, it seems reasonable to believe that the influence of spacecraft operations on plasma parameters was minimized. The average spectrum of fluctuations is in agreement with theoretical predictions. According to purposes of SAMPIE, the samples of solar cells were placed in the cargo bay of the shuttle, and high negative bias voltages were applied to them to initiate arcing between these cells and surrounding plasma. The arcing onset was registered by special counters, and data were obtained that included the amplitudes of current, duration of each arc, and the number of arcs per one experiment. The LP data were analyzed for two different situations: with arcing and without arcing. Electrostatic noise spectra for both situations and theoretical explanation of the observed features are presented in this report.

Description: The cause of the apparent small friction exhibited by long runout landslides has long been speculated upon. In an attempt to provide some insight into the matter, this paper describes results obtained from a discrete particle computer simulation of landslides composed of up to 1,000,000 two-dimensional discs. While simplified, the results show many of the characteristics of field data (the volumetric effect on runout, preserved strata, etc.) and with allowances made for the two-dimensional nature of the simulation, the runouts compare well with those of actual landslides. The results challenge the current view that landslides travel as a nearly solid block riding atop a low friction basal layer. Instead, they show that the mass is completely shearing and indicate that the apparent friction coefficient is an increasing function of shear rate. The volumetric effect can then be understood. With all other conditions being equal, different size slides appear to travel with nearly the same average velocity; however, as the larger landslides are thicker, they experience smaller shear rates and correspondingly smaller frictional resistance.

Description: The extended flight of the Airborne Ionospheric Observatory during the Geospace Environment Modeling (GEM) Pilot program on January 16, 1990, allowed continuous all-sky monitoring of the two-dimensional ionospheric footprint of the northward interplanetary magnetic field (IMF) cusp in several wavelengths. Especially important in determining the locus of magnetosheath electron precipitation was the 630.0-nm red line emission. The most striking morphological change in the images was the transient appearance of zonally elongated regions of enhanced 630.0-nm emission which resembled 'rays' emanating from the centroid of the precipitation. The appearance of these rays was strongly correlated with the Y component of the IMF: when the magnitude of B(sub y) was large compared to B(sub z), the rays appeared; otherwise, the distribution was relatively unstructured. Late in the flight the field of view of the imager included the field of view of flow measurements from the European incoherent scatter radar (EISCAT). The rays visible in 630.0-nm emission exactly aligned with the position of strong flow jets observed by EISCAT. We attribute this correspondence to the requirement of quasi-neutrality; namely, the soft electrons have their largest precipitating fluxes where the bulk of the ions precipitate. The ions, in regions of strong convective flow, are spread out farther along the flow path than in regions of weaker flow. The occurrence and direction of these flow bursts are controlled by the IMF in a manner consistent with newly opened flux tubes; i.e., when absolute value of B(sub y) greater than absolute value of B(sub z), tension in the reconnected field lines produce east-west flow regions downstream of the ionospheric projection of the x line. We interpret the optical rays (flow bursts), which typically last between 5 and 15 min, as evidence of periods of enhanced dayside (or lobe) reconnection when absolute value of B(sub y) greater than absolute value of B(sub z). The length of the reconnection pulse is difficult to determine, however, since strong zonal flows would be expected to persist until the tension force in the field line has decayed, even if the duration of the enhanced reconnection was relatively short.

Description: A general study of the structure and stability of intermediate shocks (IS) in an isotropic plasma is presented using a hybrid as well as a resistive Hall MHD code. Special emphasis is put on the question of whether the rotational layers observed at the magnetopause can be intermediate shocks. The shocks are formed dynamically by the interaction between a flowing plasma and a stationary piston. Coplanar ISs (both strong and weak) are found to be stable in a collisionless plasma. The existence of slow shocks in a high beta plasma is also established for the first time. Noncoplanar ISs are found to be time-dependent, evolving toward a rotational discontinuity (RD) after some characteristic time tau which can be quite long (1000 Omega(exp -1), where Omega is the ion gyrofrequency). The value tau is larger the closer the rotation angle is to 180 deg. Rotations larger than 180 deg are found to be unstable, decaying into a state of minimum shear (i.e., rotation angle less than 180 deg). There are various length scales associated with an IS in the kinetic regime. The shortest scale is found to be the length scale over which rotation of the transverse component of the magnetic field takes place. This scale can have a half width as small as one ion inertial length (c/omega(sub p)) for electron sense rotations and 3c/omega(sub p) for ion sense rotations, for an upstream ion beta of unity. Both of these scales are consistent with the observed thickness at the magnetopause and identical to the corresponding RD scales. A detailed study of the mode conversion of the Alfven ion cyclotron waves (A/IC) waves across both slow and intermediate shocks and the resulting downstream wave spectrum are presented. The possibility that the large number of relfected ions observed at the magnetopause may be due to the presence of strong ISs is considered. The identification of strong ISs and their distinction from RDs should be possible in observations due to significant differences that exist between jump conditions and overall structure of the two discontinuities. The jumps in the plasma parameters across a weak IS are typically small. This together with the fact that the weak ISs and RDs have very similar thickness and other overall properties makes the distinction between weak ISs and RDs in the observations largely inconsequential. However, at large noncoplanarity angles the weak IS approaches the RD limit in a relatively short time (approximately less than 100 Omega (exp -1)). Thus, magnetopause rotations with large noncoplanarity angles are most likely either RDs or strong ISs. Finally, direct comparisons between fluid (resistive Hall MHD) and kinetic simulations show that fluid theory is not applicable to study of ISs in a collisionless plasma.

Description: Orbit-averaged geomagnetic transmission measurements during the large solar energetic particle events of October 1989 are presented using proton data from the NOAA-10 and GOES-7 satellies. The measurements are compared to geomagnetic transmission calculations determined by tracing particle trajectories through the combination of the International Geomagnetic Reference Field (IGRF) model and the 1989 Tsyganenko magnetospheric magnetic field model. The effective 'ring current' parameter in the 1989 Tsyganenko model based on the Dst data. Results are compared to calculations employing only the IGRF and to a parameterization of geomagnetically quiet-time cutoff rigidities derived from Cosmos/intercosmos observations. The 3-hour orbit-averaged results have approximately 15% accuracy during the October 1989 events.

Description: Low altitude (less than 1000 km) measurements of ions precipitating into the morning auroral region are presented and analyzed. The ion fluxes exhibited time-energy signatures consistent with impulsive injection onto high-altitude field lines, followed by time-of-flight dispersion. The origin of these ions is investigated through the detailed examination of these signatures in conjunction with simultaneous measurements of precipitating electrons and a magnetic field model. A model is developed which indicates that the source for these particles was located in or near the magnetopause boundary layer, with the position deduced to be in the midlatitude flank region about 20-30 R(sub E) tailward of the Earth. The model explains the existence of multiple injections on a given field line as due to a quasi-periodic source, with the periodicity being about 100-200 s at the source. Several mechanisms are examined in an attempt to explain the injections, with a mechanism related to the propagation of waves on the surface of the boundary layer found to be the most plausible. The observations and results are compared to those of similar experiments and some unifying ideas are discussed.

Description: The assimilative mapping of ionospheric electrodynamics (AMIE) technique has been used to estimate global distributions of high-latitude ionospheric convection and field-aligned current by combining data obtained nearly simultaneously both from ground and from space. Therefore, unlike the statistical patterns, the 'snapshot' distributions derived by AMIE allow us to examine in more detail the distinctions between field-aligned current systems associated with separate magnetospheric processes, especially in the dayside cusp region. By comparing the field-aligned current and ionospheric convection patterns with the corresponding spectrograms of precipitating particles, the following signatures have been identified: (1) For the three cases studied, which all had an IMF with negative y and z components, the cusp precipitation was encountered by the DMSP satellites in the postnoon sector in the northern hemisphere and in the prenoon sector in the southern hemisphere. The equatorward part of the cusp in both hemispheres is in the sunward flow region and marks the beginning of the flow rotation from sunward to antisunward. (2) The pair of field-aligned currents near local noon, i.e., the cusp/mantle currents, are coincident with the cusp or mantle particle precipitation. In distinction, the field-aligned currents on the dawnside and duskside, i.e., the normal region 1 currents, are usually associated with the plasma sheet particle precipitation. Thus the cusp/mantle currents are generated on open field lines and the region 1 currents mainly on closed field lines. (3) Topologically, the cusp/mantle currents appear as an expansion of the region 1 currents from the dawnside and duskside and they overlap near local noon. When B(sub y) is negative, in the northern hemisphere the downward field-aligned current is located poleward of the upward current; whereas in the southern hemisphere the upward current is located poleward of the downward current. (4) Under the assumption of quasi-steady state reconnection, the location of the separatrix in the ionosphere is estimated and the reconnection velocity is calculated to be between 400 and 550 m/s. The dayside separatrix lies equatorward of the dayside convection throat in the two cases examined.

Description: The existence of localized regions of intense lower hybrid waves in the auroral ionosphere recently observed by rocket and satellite experiments can be understood by the study of a non-linear two-timescale coupling process. In this Letter, we demonstrate that the leading non-linear term in the standard Musher-Sturman equation vanishes identically in strict two-dimensions (normal to the magnetic field). Instead, the new two-dimensional equation is characterized by a much weaker non-linear term which arises from the ponderomotive force perpendicular to the magnetic field, particularly that due to the ions. The old and new equations are compared by means of time-evolution calculations of wave fields. The results exhibit a remarkable difference in the evolution of the waves as governed by the two equations. Such dissimilar outcomes motivate our investigation of the limitation of Musher-Sturman equation in quasi-two-dimensions. Only within all these limits can Musher-Sturman equation adequately describe the collapse of lower hybrid waves.

Description: The field line interhemispheric plasma (FLIP) model is used to study the 6300 A line intensity measured during three morning twilights from the McDonald Observatory in Texas. The Imaging Spectrometric Observatory (ISO) measured the 6300 A intensity during the winter of 1987 and the spring and summer of 1988. The FLIP model reproduces the measured intensity and its variation through the twilight well on each day using neutral densities from the MSIS-86 empirical model. This is in spite of the fact that different component sources dominate the integrated volume emission rate on each of the days analyzed. The sensitivity of the intensity to neutral composition is computed by varying the N2, O2, and O densities in the FLIP model and comparing to the intensity computed with the unmodified MSIS-86 densities. The ion densities change self-consistently. Thus the change in neutral composition also changes the electron density. The F2 peak height is unchanged in the model runs for a given day. The intensity changes near 100 deg SZA are comparable to within 10% when either (O2), (N2), or (O) is changed, regardless of which component source is dominant. There is strong sensitivity to changes in (N2) when dissociative recombination is dominant, virtually no change in the nighttime (SZA greater than or equal to 108 deg) intensity with (O2) doubled, and sensitivity of over 50% to doubling or halving (O) at night. When excitation by conjugate photoelectrons is the dominant nighttime component source, the relative intensity change with (O) doubled or halved is very small. This study shows the strong need for simultaneous measurements of electron density and of emissions proportional to photoelectron fluxes if the 6300 A twilight airglow is to be used to retrieve neutral densities.

Description: Ion drift meter observations from the Atmosphere Explorer E (AE-E) satellite during the period of January 1977 to December 1979 are used to study the dependence of equatorial (dip latitudes less than or equal to 7.5 deg) F region vertical plasma drifts (east-west electric fields) on solar activity, season, and longitude. The satellite-observed ion drifts show large day-to-day and seasonal variations. Solar cycle effects are most pronounced near the dusk sector with a large increase of the prereversal velocity enhancement from solar minimum to maximum. The diuurnal, seasonal, and solar cycle dependence of the logitudinally averaged drifts are consistent with results from the Jicamarca radar except near the June solstice when the AE-E nighttime downward velocities are significantly smaller than those observed by the radar. Pronounced presunrise downward drift enhancements are often observed over a large longituudinal range but not in the Peruvian equatorial region. The satellite data indicate that longitudinal variations are largest near the June solstice, particularly near dawn and dusk but are virtually absent during equinox. The longitudinal dependence of the AE-E vertical drifts is consistent with results from ionosonde data. These measurements were also used to develop a description of equatorial F region vertical drifts in four longitudinal sectors.

Description: By comparing data from the Special Sensor Microwave Imager (SSM/I) to field data, a melt threshold of the cross-polarized gradient ratio (XPGR), which is a normalized difference between the 19 GHz horizontally-polarized and 37 GHz vertically polarized brightness temperatures, is determined. This threshold, XPGR = -0.025, is used to classify dry and wet snow. The annual areal extent of melt is mapped for the years 1988 through 1991, and inter-annual variations of melt extent are examined. The results show that the melt extent varied from a low of 38.3% of the ice sheet (1990) to a high of 41.7% (1991) during the years 1988-1991.

Description: Dynamic isolation of the winter Arctic circumpolar vortex during 1992-1993 and 1993-1994 (the second and third northern hemisphere winters of the UARS mission) is studied using quasi-horizontal isentropic trajectories. Ejection of vortex air and entrainment of mid-latitude air into the vortex are quantified and compared with climatological values obtained from the analysis of 16 Arctic winters. A number of unusual features of both winters are discussed. The most notable features are the anomalous isolation experienced by the vortex during December 1992 and the unusual degree of isolation and persistence of the vortex during February and March of both years. the 1992-1993 winter season is the most consistently isolated vortex on record. Only during January 1993, when entrainment is large, is this pattern of extreme isolation broken.

Description: We present a statisical survey of Prognoz 10 solar wind observations at the times of transient (step function and impulsive) variations in the dayside magnetospheric magnetic field strength measured by the GOES 5 and 6 geosynchronous satellites. The results indicate that 51% of the magnetospheric events can be associated with corresponding variations in the solar wind dynamic pressure. A further 17% of the events can be associated with fluctuations in the interplanetary magnetic field orientation in the sense previously associated with foreshock pressure pulses. We find no tendency for impulsive events at dayside geosynchronous orbit to be associated with north/south fluctuations in the interplanetary magnetic field (IMF) orientation, nor for the events to occur primarily during intervals of southward IMF. The success rate for associating transient events at dayside geosynchronous orbit with solar wind features decreases as Prognoz 10 moves farther from the Earth-Sun line. The observations indicate that variations in the solar wind dynamic pressure and foreshock pressure pulses associated with variations in the IMF cone angle are the predominant causes of large-amplitude transient events observed at dayside geosynchronous orbit.

Description: Plasmoids are thought to occur as a consequence of the formation of a near-Earth neutral line during the evolution of a geomagnetic substorm. Using a 3D, global MHD simulation of the interaction of the Earth's magnetosphere with the solar wind, we initiate a substorm by a southward turning of the Interplanetary Magnetic Field (IMF) after a long period of steady northward field. A large plasmoid is formed and ejected. We show field line maps of its shape and relate its formation time to the progress of the substorm as indicated by the cross polar potential. Because of the large region of closed field in the magnetotail at the time of the substorm, this plasmoid is longer in axial dimension than is typically observed. We compare the simulation results with the type of satellite observations which have been used to argue for the existence of plasmoids or of traveling compression regions (TCRs) in the lobes or magnetosheath. The simulation predicts that plasmoid passage would result in a strong signal in the cross tail electric field.

Description: Plasma and magnetic field data from the International Sun Earth Explorer (ISEE) 2 spacecraft recorded on 29 Oct 1979 provide evidence for a slow shock (SS) in the reconnection layer of the dayside magnetopause. This layer is bounded on the magnetosheath side by the SS and on the magnetospheric side by a rotational discontinuity (RD). The direction of the accelerated plasma flow, the earthward sense of the normal magnetic field across both discontinuities, and the relative orientation of the SS and the RD all indicate that the reconnection site was located south of the spacecraft. Examination of the substantial pressure anisotropy downstream of the SS explains two unusual properties of the shock: (1) the slow-mode and intermediate-mode phase speeds are inverted downstream of the SS such that the RD propagates behind the SS rather than ahead of it; (2) the magnetic wave polarization reserves such that the SS initially displays a left-handed polarization and then switches to a right-handed polarization inside the shock structure.

Description: Measurements of crustal uplift from bedrock around the edges of Antarctica or Greenland could help constrain the mass balance of those ice caps. Present-day changes in ice could cause vertical displacement rates of several mm/yr around Antarctica and up to 10-15 mm/yr around Greenland. Horizontal displacement rates are likely to be about 1/3 the vertical rates. The viscoelastic response of the earth to past changes in ice could cause uplift rates that are several times larger. By measuring both gravity and vertical displacements, it is possible to remove the viscoelastic effects, so that the observations can be used to constrain present-day thickness changes.

Description: Recent updates and extensions to a steady-state two-dimensional linearized model of global-scale atmospheric waves have facilitated improved calculations of those which are subharmonics of a solar day and propagate with the apparent motion of the sun. The model improvements are briefly described and some updated predictions of the migrating solar diurnal component are highlighted. The latter represent the first numerical modeling effort to examining the seasonal variability of the migrating diurnal harmonic as it propagates into the mesosphere and lower thermosphere.

Description: Numerical simulations of mantle convection with a composite temperature-dependent, Newtonian and non-Newtonian creep law have revealed a transition in the dominant creep mechanism with the increasing vigour of convection. Newtonian creep is found to dominate in the low Rayleigh number regime. With sufficiently high effective Rayleigh number, the overall creep mechanism in the convective flow becomes non-Newtonian. The transitional Rayleigh number increases strongly with the activation energy. These results would suggest a scenario that in the early epochs of Earth the flow in the mantle would have been governed by non-Newtonian rheology and would have exhibited both strong spatial and temporal fluctuations. With time the flow mechanism would behave like a Newtonian fluid and would have a different time-dependent character. In time-dependent Newtonian-dominated flows there are still localized features with distinctly non-Newtonian character. Our analysis of the relative contributions to the lateral viscosity field supports the idea that the inference of the nature of lateral viscosity heterogeneities by seismic tomography may be strongly contaminated by the dominant non-Newtonian contributions to the total lateral viscosity field.

Description: Test particle calculations are used to compute the effects of gravity and ponderomotive acceleration by shear Alfven wave oscillations on the distribution function of O(+) ions along auroral field lines, assuming an ionospheric Maxwellian source of the ions at 2000 km altitude with approximately 0.5 eV of thermal energy in the parallel component of velocity. The electric field model corresponds to a standing wave oscillation with a frequency approximately 1 Hz in the azimuthal direction superimposed on the background dipole field, in which the wave amplitude is either increasing or decreasing in time. The electric field is taken to be primarily in the perpendicular direction. The time varying wave produces broad distributions with widths of 2 to 10 times the initial 0.5-eV thermal energy of the Maxwellian source, and the density and flux of upward going O(+) ions at one Earth radius are both enhanced in this model. The oxygen ion distribution functions at 1 R(sub E) altitude resulting from interaction with waves whose amplitudes are increasing in time have a more gradual lower energy cutoff than do the distribution functions resulting from decaying waves. The high-energy part of the distribution functions in growing waves reflects the temperature of the Maxwellian source, while the high-energy part of the distributions resulting from decaying waves steepens with time, independent of the source temperature.

Description: The structure of the magnetopause is investigated by studying the interaction between two plasmas with solar wind and magnetospheric properties. Both Hall-MHD and hybrid (particle ions and fluid electrons) simulations are performed to compare and contrast the nature of the solutions in the fluid and the kinetic limits. It is shown that, in accordance with previous studies, the fluid solutions consist of multiple discontinuities and waves such as slow shocks and waves, as well as contact and rotational discontinuities. In contrast, the kinetic solutions consist of fewer discontinuities and include non-MHD boundaries. The difference between the two types of solutions are attributed to the absence of contact surfaces in collisionless plasmas and the possible Landau damping of slow waves. The kinetic solutions are found to be in a much better agreement with magnetopause observations, which have shown little evidence for the presence of slow shocks or contact discontinuities. The results of kinetic calculations suggest that the presence of a small but finite normal component of the magnetic field allows for the mixing of magnetosheath and magnetospheric plasmas and may, in part, be responsible for the formation of the boundary layer even during the periods of northward interplanetary magnetic field. These results also show that the necessary changes in the fields and the plasma parameters across the magnetopause do not occur simultaneously (i.e., in the same spatial location). As such, the magnetopause cannot be viewed as a single boundary but instead is a region which includes both abrupt and gradual changes in the fields and plasma parameters.

Description: Electron distributions at energies above 50 eV have been found to be a sensitive indicator of magnetic topology for magnetopause crossings of the AMPTE/CCE spacecraft. Progressing from the magnetosheath to the magnetosphere two abrupt transitions occur. First, the magnetosheath electron population directed either parallel or antiparallel to the magnetic field is replaced by a streaming, heated magnetosheath electron population. The other half of the distribution is unchanged. The region with unidirectional, heated magnetosheath electrons is identified as the magnetosheath boundary layer (MSBL). Second, the unheated magnetosheath electron population is replaced by a heated population nearly identical to the population encountered in the MSBL, resulting in a symmetric counterstreaming distribution. The region populated by the bidirectional heated magnetosheath electrons is identified as the low-latitude boundary layer (LLBL). The MSBL and LLBL identified by the electron transitions are the same as the regions identified using ion composition measurements. The magnetosheath-MSBL transition reflects a change in magnetic topology from a solar wind field line to one that threads the magnetopause, and the existence of a magnetosheath-MSBL transition implies that the magnetopause is open. When the current layer is easily identified, the MSBL-LLBL transition coincides with the magnetopause current layer, indicating that the magnetosheath electrons are heated in the current layer. Both magnetosheath-MSBL and MSBL-LLBL transitions are observed for low as well as high magnetic shears. Moreover, the transitions are particularly clear for low shear implying that magnetic topology boundaries are sharp even when abrupt changes in the field and other plasma parameters are absent. Furthermore, for low magnetic shear, solar wind ions with low parallel drift speeds make up the majority of the LLBL population indicating that the magnetosheath plasma has convected directly across the magnetosheath plasma has converted directly across the magnetopause. These observations are consistent with quasi-steady, high-latitude reconnection and indicate that the signatures of this reconnection geometry are commonly present in the subpolar region.

Description: We have developed the linear theory of collisionless ion tearing in a two-dimensional magnetotail equilibrium for a single resonant species. We have solved the normal mode problem for tearing instability by an algorithm that employs particle-in-cell simulation to calculate the orbit integrals in the Maxwell-Vlasov eigenmode equation. The results of our single-species tearing analysis can be applied to ion tearing where electron effects are not included. We have calculated the tearing growth rate as a function of the magnetic field component B(sub n) normal to the current sheet for thick and thin current sheets, and we show that marginal stability occurs when the normal gyrofrequency Omega(sub n) is comparable to the Harris neutral sheet growth rate. A cross-tail B(sub y) component has little effect on the growth rate for B(sub y) approximately = B(sub n). Even in the limit B(sub y) much greater than B(sub n), the mode is strongly stabilized by B(sub n). We report than random pitch angle scattering can overcome the stabilizing effect of B(sub n) and drive the growth rate up toward the Harris neutral sheet (B(sub n) = 0) value when the pitch angle diffusion rate is comparable to Omega(sub n).

Description: A model of mantle convection which generates plate tectonics requires strain rate- or stress-dependent rheology in order to produce strong platelike flows with weak margins as well as strike-slip deformation and plate spin (i.e., toroidal motion). Here, we employ a simple model of source-sink driven surface flow to determine the form of such a rheology that is appropriate for Earth's present-day plate motions. In this model, lithospheric motion is treated as shallow layer flow driven by sources and sinks which correspond to spreading centers and subduction zones, respectively. Two plate motion models are used to derive the source sink field. As originally implied in the simpler Cartesian version of this model, the classical power law rheologies do not generate platelike flows as well as the hypothetical Whitehead-Gans stick-slip rheology (which incorporates a simple self-lubrication mechanism). None of the fluid rheologies examined, however, produce more than approximately 60% of the original maximum shear. For either plate model, the viscosity fields produced by the power law rheologies are diffuse, and the viscosity lows over strike-slip shear zones or pseudo-margins are not as small as over the prescribed convergent-divergent margins. In contrast, the stick-slip rheology generates very platelike viscosity fields, with sharp gradients at the plate boundaries, and margins with almost uniformly low viscosity. Power law rheologies with high viscosity contrasts, however, lead to almost equally favorable comparisons, though these also yield the least platelike viscosity fields. This implies that the magnitude of toroidal flow and platelike strength distributions are not necessarily related and thus may present independent constraints on the determination of a self-consistent plate-mantle rheology.

Description: We present measurements of the middle atmospheric water vapor mixing ratio profile obtained using the ground-based Naval Research Laboratory (NRL) water vapor millimeter-wave spectrometer (WVMS) instrument at the Jet Propulsion Laboratory (JPL) Table Mountain Observatory. The measurements cover a period of 262 days from January 23, 1992, to October 13, 1992. During this campaign it was possible to retrieve useful daily mixing ratio profiles for 186 days. We thus have a nearly continuous record of water vapor mixing ratios for altitudes from approximately equals 35 to 75 km. The retrievals are obtained using the optimal estimation method. Details of the error analysis are presented, and a technique is introduced that reduces baseline effects and helps to estimate the baseline error. The high-altitude (greater than or approx. equal to 65 km) data show a sharp rise prior to the expected maximum near the summer solstice and a gradual decline in the following months. The mixing ratios generally peak between 55 and 65 km, at which point the mixing ratios are 6 - 7 parts per million by volume. The highest peaks occur in January, May, and October.

Description: During December 1988, 24 hours of darkness and clear sky conditions permitted continuous observations of the O I(6300 A) airglow by a Fabry-Perot interferometer located at Thule Air Base, Greenland. Thus a continuous record of the F region neutral winds was obtained for that month. During this same time period, a digital ionosonde located at Qanaq, Greenland (110 km north of Thule Air Base), was in operation measuring electron density profiles and F region ion drifts. This combination of ground-based observations allowed the investigation of ion/neutral coupling at a temporal resolution of about 15 min. Interplanetary magnetic field (IMF) data from the IMP 8 satellite were also available from December 16 to 24 and indicated intervals of B(sub z) northward IMF conditions during this period. Here we investigate the observed response of the neutral wind to convection changes in the ion drift inside the polar cap for southward and northward IMF B(sub z) conditions. In particular, we establish a control day illustrating the typical antisunward neutral wind and ion drift patterns observed for southward B(sub z) over Thule and Qanaq, and we compare it with observations made when the IMF B(sub z) is directed northward. The observations during periods of northward B(sub z) display sunward directed ion drifts over the polar cap accompanied by decreasing antisunward directed neutral winds. We investigate these times of northward B(sub z) further and demonstrate that the ion drag term alone cannot describe the observed response in the neutral wind during northward IMF.

Description: Plasmaspheric densities have been observed previously to be higher in December than in June, with the ratio varying between 1.5 and 3.0 and with larger variations at lower L shells. In order to search for the cause of the observed annual variations, we have modeled plasmaspheric density, using a time-dependent hydrodynamic model. On an L = 2 field line with geomagnetic longitude equal to 300 deg, the modeled plasmaspheric densities were a factor of 1.5 times higher in December than in June. The modeled December to June density ratio was found to increase slightly with L shell, in contrast to observations; this discrepancy may be due to the fact that outer plasmaspheric flux tubes are never completely full. In addition, for an L = 2 field line with geomagnetic longitude equal to 120 deg, the modeled plasmaspheric density was higher in June than in December by a factor of about 1.2. Various numerical tests were also performed in order to examine the sensitivity of plasmaspheric density to various parameters. In particular, a large vertical neutral wind was applied in order to raise the O(+) profile which had the effect of raising plasmaspheric density by a factor of 6. This in conjunction with a theoretical analysis suggests that plasmaspheric density levels are very sensitive to O(+) levels in the upper ionosphere. We conclude that annual variations in plasmaspheric density are due to similar variations in ionospheric O(+).

Description: Continuous ground-based observations of ionospheric and magnetospheric regions are critical to the Geospace Environmental Modeling (GEM) program. It is therefore important to establish clear intercalibrations between different ground-based instruments and satellites in order to clearly place the ground-based observations in context with the corresponding in situ satellite measurements. HF-radars operating at high latitudes are capable of observing very large spatial regions of the ionosphere on a nearly continuous basis. In this paper we report on an intercalibration study made using the Polar Anglo-American Conjugate Radar Experiment radars located at Goose Bay, Labrador, and Halley Station, Antarctica, and the Defense Meteorological Satellite Program (DMSP) satellites. The DMSP satellite data are used to provide clear identifications of the ionospheric cusp and the low-latitude boundary layer (LLBL). The radar data for eight cusp events and eight LLBL events have been examined in order to determine a radar signature of these ionospheric regions. This intercalibraion indicates that the cusp is always characterized by wide, complex Doppler power spectra, whereas the LLBL is usually found to have spectra dominated by a single component. The distribution of spectral widths in the cusp is of a generally Gaussian form with a peak at about 220 m/s. The distribution of spectral widths in the LLBL is more like an exponential distribution, with the peak of the distribution occurring at about 50 m/s. There are a few cases in the LLBL where the Doppler power spectra are strikingly similar to those observed in the cusp.

Description: Ten transient magnetic structures in Earth's magnetotail, as observed in GEOTAIL measurements, selected for early 1993 (at (-) X(sub GSM) = 90 - 130 Earth radii), are shown to have helical magnetic field configurations similar to those of interplanetary magnetic clouds at 1 AU but smaller in size by a factor of approximately = 700. Such structures are shown to be well approximated by a comprehensive magnetic force-free flux-rope model. For this limited set of 10 events the rope axes are seen to be typically aligned with the Y(sub GSM) axis and the average diameter of these structures is approximately = 15 Earth radii.

Description: The problem of the exitation of obliquely propagating magnetosonic waves which can steepen up (also known as shocklets) is considered. Shocklets have been observed upstream of the Earth's bow shock and at comets Giacobini-Zinner and Grigg-Skjellerup. Linear theory as well as two-dimensional (2-D) hybrid (fluid electrons, particle ions) simulations are used to determine the properties of waves generated by ring-beam velocity distributions in great detail. The effects of both proton and oxygen ring-beams are considered. The study of instabilities excited by a proton ring-beam is relevant to the region upstream of the Earth's bow shock, whereas the oxygen ring-beam corresponds to cometary ions picked up by the solar wind. Linear theory has shown that for a ring-beam, four instabilities are found, one on the nonresonant mode, one on the Alfven mode, and two along the magnetosonic/whistler branch. The relative growth rate of these instabilities is a sensitive function of parameters. Although one of the magnetosonic instabilities has maximum growth along the magnetic field, the other has maximum growth in oblique directions. We have studied the competition of these instabilities in the nonlinear regime using 2-D simulations. As in the linear limit, the nonlinear results are a function of beam density and distribution function. By performing the simulations as both initial value and driven systems, we have found that the outcome of the simulations can vary, suggesting that the latter type simulations is needed to address the observations. A general conclusion of the simulation results is that field-aligned beams do not result in the formation of shocklets, whereas ring-beam distributions can.

Description: In order to facilitate bounce-averaged guiding center simulations of geomagnetically trapped particles, we express the kinetic energy of a particle with magnetic coordinates (L,phi) as an analytic function of the first two adiabatic invariants (M,J) and the L value of the field line. The magnetic field model is axisymmetric, consisting of a dipolar vector-B field plus a uniform southward magnetic field parallel to the dipole moment mu(sub E). This model magnetosphere is surrounded by a circular equatorial neutral line whose radius b is an adjustable parameter. Our formulation provides a computationally efficient method for tracing the bounce-averaged adiabatic motion (conserving all three invariants) and nonadiabatic transport (violating the third invariant while conserving the first two invariants) of geomagnetically trapped particles in the model magnetosphere.

Description: The Earth's bow shock is a bountiful generator of waves. Some of these waves have group velocities that exceed the solar wind velocity directed into the shock and can propagate upstream against the flow. Upstream whistlers observed close to one Hertz in the spacecraft frame have been seen many Earth radii upstream. A second whistler mode wave, called the precursor, propagates upstream along the shock normal but is phase standing in the solar wind flow. The damping of both whistler mode waves is consistent with Landau damping. At low Mach numbers the precursor is connected to the non-coplanarity component in the shock ramp. At higher Mach numbers the upstream waves cannot propagate upstream and ion reflection becomes more important in providing free energy for wave particle interactions. The non-coplanarity component is still present but it now initiates a downstream wave train. Generally the waves just downstream from the bow shock are left hand circularly polarized ion cyclotron waves propagating along the magnetic field at the Alfven velocity. When the upstream Mach number is high and the helium content of the plasma is high, mirror mode waves are observed.

Description: We review our current knowledge of ULF waves in planetary foreshocks. Most of this knowledge comes from observations taken within a few Earth radii of the terrestrial bow shock. Terrestrial foreshock ULF waves can be divided into three types, large amplitude low frequency waves (approximately 30-s period), upstream propagating whistlers (1-Hz waves), and 3-s waves. The 30-s waves are apparently generated by back-streaming ion beams, while the 1-Hz waves are generated at the bow shock. The source of the 3-s waves has yet to be determined. In addition to issues concerning the source of ULF waves in the foreshock, the waves present a number of challenges, both in terms of data acquisition, and comparison with theory. The various waves have different coherence scales, from approximately 100 km to approximately 1 Earth radius. Thus multi-spacecraft separation strategies must be tailored to the phenomenon of interest. From a theoretical point of view, the ULF waves are observed in a plasma in which the thermal pressure is comparable to the magnetic pressure, and the rest-frame wave frequency can be moderate fraction of the proton gyro-frequency. This requires the use of kinetic plasma wave dispersion relations, rather than multi-fluid MHD. Lastly, and perhaps most significantly, ULF waves are used to probe the ambient plasma, with inferences being drawn concerning the types of energetic ion distributions within the foreshock. However, since most of the data were acquired close to the bow shock, the properties of the more distant foreshock have to be deduced mainly through extrapolation of the near-shock results. A general understanding of the wave and plasma populations within the foreshock, their interrelation, and evolution, requires additional data from the more distant foreshock.

Description: Whether a dike can propagate far from a magma reservoir depends upon the competition between the rate at which propagation widens the dike and the rate at which freezing constricts the aperture available for magma flow. Various formulations are developed for a viscous fluid at temperature T(sub m) intruding a growing crack in an elastic solid. The initial solid temperature equals T(sub m) at the source and decreases linearly with distance from the source. If T(sub m) is the unique freezing temperature of the fluid, dike growth is initially self-similar and an essentially exact solution is obtained; if T(sub m) is above the solidus temperature, the solution is approximate but is designed to overestimate the distance the dike may propagate. The ability of a dike to survive thermally depends primarily upon a single parameter that is a measure of the ratio of the dike frozen margin thickness to elastic thickness. Perhaps more intuitively, one may define a minimum distance from the essentially solid reservoir wall to the point at which the host rock temperature drops below the solidus, necessary for dikes to propagate far into subsolidus rock. It is concluded that for reasonable material properties and source conditions, most basalt dikes will have little difficulty leaving the source region, but most rhyolite dikes will be halted by freezing soon after the magma encounters rock at temperatures below the magma solidus. While these results can explain why granitic dikes are common near granitic plutons but rare elsewhere, the potentially large variation in magmatic systems makes it premature to rule out the possibility that most granites are transported through the crust in dikes. Nonetheless, these results highlight difficulties with such proposals and suggest that it may also be premature to rule out the possibility that most granite plutons ascend as more equidimensional bodies.

Description: The early time evolution of an ionospheric electron depletion produced by a radially expanding electron attachment chemical release is studied with a two-dimensional simulation model. The model includes electron attachment chemistry, incorporates fluid electrons, particle ions and neutrals, and considers the evolution in a plane perpendicular to the geomagnetic field for a low beta plasma. Timescales considered are of the order of or less than the cyclotron period of the negative ions that result as a by-product of the electron attacment reaction. This corresponds to time periods of tenths of seconds during recent experiemts. Simulation results show that a highly sheared azimuthal electron flow velocity develops in the radially expanding depletion boundary. This sheared electron flow velocity and the steep density gradients in the boundary give rise to small-scale irregulatities in the form of electron density cavities and spikes. The nonlinear evolution of these irregularities results in trapping and ultimately turbulent heating of the negative ions.

Description: We study gradients of the energetic ion intesity observed at the edge of the plasma sheet boundary layer (PSBL) by the energetic ion anisotropy spectrometer (EPAS) on International Sun Earth Explorer 3 (ISEE 3). In particular, we have determined the velocity of the boundary relative to the spacecraft in the direction perpendicular to the tail axis and the angle which the boundary normal makes to the spacecraft spin axis for 1160 PSBL encounters at X(sub GSM) greater than -240 R(sub E). By asuming that, on average, the edge of the PSBL is parallel to the cross-tail current sheet, we are then able to determine a number of properties of the structure, orientation and motion of the deep geomagnetic tail. We conclude the following: (1) Most crossings of the edge of the PSBL are caused by transverse motuion of the entire tail induced by solar wind direction variations, although some are caused by reconfiguration of the tail due to geomagnetic activity. (2) The typical velocity of the PSBL (and hence of the tail) in the direction perpendicular to the tail axis is 50-85 km/s. (3) The average twist of the tail is near zero, with the edge of the PSBL (and by inference the cross-tail current sheet) lying parallel to the ecliptic plane (however, large twists are found in individual events and the distribution of twists is broad, with one standard deviation of approximately 50 deg. (4) The width of the distribution decreases with downtail distance. (5) The variation of the distributions with cross-tail position reveals that this decreas in width is most likely due to the edge of the PSBL being concave, or significantly flared at the tail flanks, in the near-Earth region. (6) During days on which the Interplanetary Magnetic Field (IMF) has 'away' sector structure, the north lobe of the trail is twisted on average towards dawn by 7.0 +/-2.4 deg. (7) During days on which the IMF has 'toward' sector structure, the north lobe is tilted towards duskby 3.8 +/- 2.3. (8) A subset of events for which IMP 8 solar wind data are available show that, for southward IMF BH(sub z) the tail has a mean twist of -12.3 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 5.5 +/- 3.8 deg for IMF B(sub Y) less than 0 (positive twist angles correspond to a tilt of the northern lobe towards dusk). (9) For northward IMF B(sub z) the tail has a twist of -23.9 +/- 5.0 deg for IMF B(sub Y) greater than 0 and 13.4 +/- 6.0 deg for IMF B(sub y) less than 0. Hence the tail appears more twisted on average for the IMF B(sub Z) northward case. (10) The distribution of tail twist is wider for lower levels of geomagnetic activity, indicating that the tail is able to twist more at lower levels of activity. (11) The data set reveals no evident effect of the earth's dipole wobble; tail orientation appears to be controlled by the solar wind and IMF, such that the GSE coordinate system may be appropriate for the study of field and plasma structures in the distant tail region.

Description: Intermediate shocks (ISs) lead to a transition from super-Alfvenic to sub-Alfvenic flow and are different from slow and fast shocks in that an IS rotates the component of the magnetic field tangent to the shock plane by 180 deg. Another peculiarity of ISs is that for the same upstream conditions an IS can have two different downstream states. There also exist a second class of ISs which rotate the magnetic field by an angle other than 180 deg. Due to their noncoplanar nature they cannot be time-stationary and are referred to as time-dependent intermediate shocks (TDIS). The existence of ISs has been the subject of much controversy over the years. Early studies questioned the physical reality of ISs. However, the studies of ISs found a new impetus when C.C. Wu showed that ISs do exist and are stable within the resistive MHD framework. In this paper, after a brief historical overview of the subject, we will review the latest developments in the study of ISs. In particular, we will address the questions of stability and structure of ISs and the relationship between ISs and other discontinuities. One of the recent developments has been the finding that ISs can be unsteady, reforming in time. Details of this process will be discussed. Finally, we examine the effect of anisotropy on the resolutions and discuss the relevance of ISs to the observed field rotations at the Earth's magnetopause.

Description: Previous observations have shown that the solar wind is decelerated and deflected in the earth's upstream region populated by long-period waves. This deceleration is corelated with the 'diffuse' but not with the 'reflected' ion population. The speed of the solar wind may decrease tens of km/s in the foreshock region. The solar wind dynamic pressure exerted on the magnetopause may vary due to the fluctuation of the solar wind speed and density in the foreshock region. In this study, we examine this solar wind deceleration and determine how the solar wind deceleration varies in the foreshock region.

Description: It is generally believed that flux transfer events (FTEs) in the outer dayside magneosphere, usually identified by transient (approximately 1 min) bipolar magneitc field perturbations in the direction normal to the nominal magnetopause, occur when the magnetosheath magetic field has a southward component. We compare the results of three methods for determining the magnetosheath magnetic field orientationat the times of previously identified UKS/IRM events: (1) the average magnetosheath magnetic field orientation in the 30-min period adjacent to the nearest magnetopause crossing, (2) the magnetosheath magnetic field orientation observed just outside the magnetopause, and (3) the lagged interplanetary magnetic field (IMF) orientation at the time of the transient events. Whereas the results of method 2 indicate that the events tend to occur for a southward magnetosheath magnetic field, the results of methods 1 and 3 show no such tnedency. The fact that the three methods yield significantly diffeent results emphasizes the need for caution in future studies.

Description: The spatial structure of dayside large-scale field-aligned current (FAC) systems is examined by using Viking and Defense Meteorological Satellite Program-F7 (DMSP-F7) data. We focus on four events in which the satellites simultaneously observed postnoon and prenoon three FAC systems: the region 2, the region 1, and the mantle (referred to as midday region O) systems, from equatorward to poleward. These events provide the most solid evidence to date that the midday region O system is a separate and unique FAC system, and is not an extension of the region 1 system from other local times. The events are examined comprehensively by making use of a mulit-instrumental data set, which includes magnetic field, particle flux, electric field, auroral UV image data from the satellites, and the Sondrestrom convection data. The results are summarized as follows: (1) Region 2 currents flow mostly in the central plasma sheet (CPS) precipitation region, often overlapping with the boundary plasma sheet (BPD) at their poleward edge. (2) The region 1 system is located in the core part of the auroral oval and is confined in a relatively narrow range in latitude which includes the convection reversal. The low-latitude boundary layer, possibly including the outer part of the plasma sheet, and the external cusp are the major source regions of dayside region 1 currents. (2) Midday region O currents flow on open field lines and are collocated with the shear of antisunward convection flows with velocites decreasing poleward. On the basis of these results we support the view that both prenoon and postnoon current systems consist of the three-sheet structure when the disctortion ofthe convection pattern associated with interplanetary magnetic field (IMF) B(sub Y) is small and both morningside and eveningside convection cells are crescent-shaped. We also propose that the midday region O and a part of the region 1 systems are closely coupled to the same source.

Description: A one-dimensional hybrid satellite track model has been developed to calculate the high-latitude thermospheric/ionospheric structure below the satellite altitude using Dynamics Explorer 2 (DE 2) satellite measurements and theory. This model is based on Emery et al. (1985) satellite track code but also includes elements of Roble et al. (1987b) global mean thermosphere/ionosphere model. A number of parameterizations and data handling techniques are used to input satellite data from several DE 2 instruments into this model. Profiles of neutral atmospheric densities are determined from the Mass Spectrometer Incoherent Scatter 1990 (MSIS-90) model and measured neutral temperatures. Measured electron precipitation spectra are used in an auroral model to calculate particle impact ionization rates below the satellite. These rates are combined with a solar ionization rate profile and used to solve the O(+) diffusion equation, with the measured electron density as an upper boundary condition. The calculated O(+) density distribution, as well as the ionization profiles, are then used in a photochemical equilibrium model to calculate the electron and molecular ion densities. The electron temperature is also calculated by solving the electron energy equation with an upper boundary condition determined by the DE 2 measurement. The model enables calculations of altitude profiles of conductivity and Joule heating rate along and below the satellite track. In a first application of the new model, a study is made of thermospheric and ionospheric structure below the DE 2 satellite for a single orbit which occurred on October 25, 1981. The field-aligned Poynting flux, which is independently obtained for this orbit, is compared with the model predictions of the height-integrated energy conversion rate. Good quantitative agreement between these two estimates has been reached. In addition, measurements taken at the incoherent scatter radar site at Chatanika (65.1 deg N, 147.4 deg W) during a DE 2 overflight are compared with the model calculations. A good agreement was found in lower thermospheric conductivities and Joule heating rate.

Description: Melting of the upper mantle and extraction of melt result in the formation of a less dense depleted mantle. This paper describes series of two-dimensional models that investigate the effects of chemical buoyancy induced by these density variations. A tracer particles method has been set up to follow as closely as possible the chemical state of the mantle and to model the chemical buoyant force at each grid point. Each series of models provides the evolution with time of magma production, crustal thickness, surface heat flux, and thermal and chemical state of the mantle. First, models that do not take into account the displacement of plates at the surface of Earth demonstrate that chemical buoyancy has an important effect on the geometry of convection. Then models include horizontal motion of plates 5000 km wide. Recycling of crust is taken into account. For a sufficiently high plate velocity which depends on the thermal Rayleigh number, the cell's size is strongly coupled with the plate's size. Plate motion forces chemically buoyant material to sink into the mantle. Then the positive chemical buoyancy yields upwelling as depleted mantle reaches the interface between the upper and the lower mantle. This process is very efficient in mixing the depleted and undepleted mantle at the scale of the grid spacing since these zones of upwelling disrupt the large convective flow. At low spreading rates, zones of upwelling develop quickly, melting occurs, and the model predicts intraplate volcanism by melting of subducted crust. At fast spreading rates, depleted mantle also favors the formation of these zones of upwelling, but they are not strong enough to yield partial melting. Their rapid displacement toward the ridge contributes to faster large-scale homogenization.

Description: Electromagnetic wave generation and resulting cross-field diffusion of plasma are considered at a tangential discontinuity, which characterizes the magnetopause for northward interplanetary magnetic field. Two-dimensional hybrid (particle ions, massless fluid electrons) simulations, in which the tangential discontinuity is generated self-consistently via a stream-stream interaction, are used to show that wave growth occurs when the ambient magnetic field is predominantly perpendicular to the direction of the density gradient. Low-frequency (much less than ion gyrofrequency) waves, with amplitudes delta B/B less than or equal to 0.2 and anticorrelated density fluctuations delta n/n less than or equal to 0.6, are generated at the discontinuity, resulting in cross-field diffusion that is comparable to the Bohm rate. Both the fluctuation level and the lack of ion heating in the calculations are consistent with observations at the magnetopause. The magnitude of the diffusion is considered in the presence of numerical effects and in the context of the inferred diffusion rate at the magnetopause. The relation of the low-frequency waves and their consequences to faster growing, short-wavelength waves due to the lower hybrid drift instability is also addressed. The overall conclusion of this initial study is that diffusion due to low frequency waves is not likely to be a major effect at the magnetopause.

Description: High spatial resolution numerical simulations of mantle plumes impinging from below on the endothermic phase change at 660-km depth are used to investigate the effects of latent heat release on the plume-phase change interaction. Both axisymmetric and planar upflows are considered, and the strong temperature dependence of mantle viscosity is taken into account. For plume strengths considered, a Clapeyron slope of -4 MPa/K prevents plume penetration of the phase change. Plumes readily penetrate the phase change for a Clapeyron slope of -2 MPa/K and arrive in the upper mantle considerably hotter than if they had not traversed the phase change. For the same amount of thermal drive, i.e., the same excess basal temperature, axisymmetric plumes are hotter upon reaching the upper mantle than are planar upwellings. Heating of plumes by their passage through the spinel-perovskite endothermic phase change can have important consequences for the ability of the plume to thermally thin the lithosphere and cause melting and volcanism.

Description: We have studied 2-D time-dependent convection for a rheology which is both non-Newtonian and temperature-dependent. Strong effects associated with viscous heating are found in the downwelling sheets, which are heated on both sides with an intensity around O(100) times the chondritic value. The magnitude of viscous heating increases strongly with the subduction speed. The slab interior is weakened by viscous heating and slab breakoff then takes place. This process provides a self-regulating mechanism for governing the speed of intact slabs able to reach the deep mantle. Timescales associated with viscous heating are quite short, a few million years. Internal heating by radioactivity decreases the amount of shear heating.

Description: Electric field and conductivity measurements in the stratosphere between November 1992 and March 1993 have been made using superpressure balloons in the southern hemisphere. Over 400 payload-days of data have been made during a record setting experiment called ELBBO (Extended Life Balloon Borne Observatories). This experiment resulted in 4 flights aloft simultaneously for over 2 months including one flight which lasted over 4 months. Electrodynamical coupling between the atmosphere and ionosphere is studied using the measured electric fields, and a simple empirical model of the stratospheric conductivity. Altitude profiles of conductivity have been obtained from several superpressure balloon flights using the large end-of-flight altitude swings on the last few days of each flight (as the balloon begins to loose superpressure). Coupling between the fields and atmospheric inertial waves has been observed. Effects and dynamics of the global circuit suggest that standard models are missing significant phenomena. Large scale ionospheric convection activity has been studied from the polar cap to the middle latitudes. Cusp latitude fields have been continuously measured for many days in a row.

Description: Large-scale ocean observing programs such as the Joint Global Ocean Flux Study (JGOFS) and the World Ocean Circulation Experiment (WOCE) today, must face the problem of designing an adequate sampling strategy. For ocean chemical variables, the goals and observing technologies are quite different from ocean physical variables (temperature, salinity, pressure). We have recently acquired data on the ocean CO2 properties on WOCE cruises P16c and P17c that are sufficiently dense to test for sampling redundancy. We use linear and quadratic interpolation methods on the sampled field to investigate what is the minimum number of samples required to define the deep ocean total inorganic carbon (TCO2) field within the limits of experimental accuracy (+/- 4 micromol/kg). Within the limits of current measurements, these lines were oversampled in the deep ocean. Should the precision of the measurement be improved, then a denser sampling pattern may be desirable in the future. This approach rationalizes the efficient use of resources for field work and for estimating gridded (TCO2) fields needed to constrain geochemical models.

Description: This paper examines an isolated magnetospheric VLF/radio noise event that is highly suggestive of the triggering of terrestrial auroral kilometric radiation (AKR) bu solar type III radio emission and of a close relation between AKR and broadband hiss. The solar type III burst was measured on polar HF riometers and was coincident with local dayside VLF/LF noise emission bursts at South Pole station. It was also coincident with AKR bursts detected onthe AMPTE/IRM satellite, at the same magnetic local time as South Pole. On the basis of the close association of AKR and VLF bursts, and from geometric considerations relating to wave propagation, it is likely that the AKR source was on the dayside and on field lines near South Pole station. The general level of geomagnetic activity was very low. However, an isolated magnetic impulse event (MIE) accompanied by a riometer absorption pulse was in progress when all of the VLF/radio noise bursts occurred. The very close association of the typew III burst at HF with the AKR is consistent with external stimulation of the AKR, is different, more immediate,triggering process than that implied by Calvert (1981) is invoked. It is suggested here that some of the HF solar radiant energy may decay into waves with frequences comparable to those of the AKR by paraetric excitation or some other process, thus providing the few background photons required for the generation of AKR by the WU and Lee (1979) cyclotron maser instability. The AKR, perhaps by modifying the magnetospheric electron velocity distribution, might have produced the observed VLF emissions. Alternatively, the VLF emissions may have arisen from the same anisotropic and unstable electron distribution function responsible for the AKR.

Description: A system of four current sheets of large-scale field-aligned currents (FACs) was discovered in the data set of simultaneous Viking and Defense Meteorological Satellire Program-F7 (DMSP-F7) crossing of the dayside high-latitude region. This paper reports four examples of this system that were observed in the prenoon sector. The flow polarities of FACs are upward, downward, upward, and downward, from equatorward to poleward. The lowest-latitude upward current is flowing mostly in the central plasma sheet (CPS) precipitation region, often overlapping with the boundary plasma sheet (BPS) at its poleward edge, andis interpreted as a region 2 current. The pair of downward and upward FACs in the middle of te structure are collocated with structured electron precipitation. The precipitation of high-energy (greater than 1 keV) electrons is more intense in the lower-latitude downward current sheet. The highest-latitude downward flowing current sheet is located in a weak, low-energy particle precipitation region, suggesting that this current is flowing on open field lines. Simulaneous observations in the postnoon local time sector reveal the standard three-sheet structure of FACs, sometimes described as region 2, region 1, and mantle (referred to the midday region O) currents. A high correlation was found between the occurrence of the four FAC sheet structure and negative interplanetary magnetic field (IMF) B(sub Y). We discuss the FAC structurein terms of three types of convection cells: the merging, viscous, andlobe cells. During strongly negative IMF B(sub Y), two convection reversals exist in the prenoon sector; one is inside the viscous cell, and the other is between the viscous cell and the lobe cell. This structure of convection flow is supported by the Viking electric field and auroral UV image data. Based on the convection pattern, the four FAC sheet structure is interpreted as the latitude overlap of midday and morning FAC systems. We suggest that the for-current sheet structure is common in a certain prenoon localtime sector during strongly negative IMF B(sub Y).

Description: We investigate the ionization of interstellar hydrogen and helium due to electron impact by shock-heated electrons. Taking the electron distributions measured at four interplanetary shocks at 1 AU, we show that the electrons in the downstream region of strong shocks can ionize interstellar atoms at rates matching or exceeding the nominal photoionization or charge-exchange rates. We suggest that this process may explain some puzzling observations of interstellar pickup ions by the Ulysses spacecraft.

Description: Onsager et al. (1991) have put forward a model of the formation of the plasma sheet boundary layer (PSBL) which relies on a steady source of plasma from a spatially extended plasma sheet, together with steady equatorward and earthward ExB convection of field lines due to reconnection at a downtail neutral line. This model is a synthesis of earlier proposals and it explains such features as an electron layer exterior to the ion boundary layer, ion velocity dispersion, counter streaming beams, low-speed cutoffs in the beams. It also explains the apparent evolution of the ion beams through 'kidney bean' shaped velocity-space distributions toward quasi-isotropic shells without invoking pitch angle scattering or energy diffusion. In this paper we explore two ramifications of the model. In principle we can map, as a function of time, the downtail neutral line distance and establish whether or not it is retreating during substorm recovery. We can also reconstruct the plasma distribution function near the neutral line to see if it is most consistent with mantle or plasma sheet plasma. We perform this analysis using International Sun Earth Explorer (ISEE) Fast Plasma Experiment (FPE) data for two plasma sheet recovery events, one on March 1, 1978, and the other on April 18, 1978. On March 1, 1978, we find evidence for an initial retreat from around 110 to 160 R(sub E) in the first 15 min; little further retreat occurs thereafter. On April 18, 1978, the neutral line location ranges from as little as 40 R(sub E) tailward of the satellite to as much as 200 R(sub E), but there is no evidence for a systematic retreat. The reconstructed ion distributions for these events are most consistent with a plasma sheet origin for the March 1 case and possibly plasma mantle or low-latitude boundary layer for the April 18 case.

Description: We have made a statistical study of the spatial distribution of low frequency waves (approx. 0.01-0.1 Hz) in the region upstream of the pre-dawn to dawn side bow shock (-50 Re less than X less than 15 Re) using both GEOTAIL and international sun earth explorer 3 (ISEE-3) magnetometer data. We have found that the wave amplitude dependence on D and X(sub s), where D is the distance from the bow shock and X(sub s) the x-coordinate position of shock foot point of the IMF, can be described by a functional form of A exp (X(sub s)/L(sub X)-D/L(sub D), with the characteristic attenuation distances, L(sub X) = 62 +/- 12 Re and L(sub D) = 59 +/- 38 Re.

Description: WINDII, the Wind Imaging Interferometer on the Upper Atmosphere Research Satellite, observes winds, temperatures and emission rates in the upper mesosphere and thermosphere. In this paper we report on nighttime observations of the vertical distribution of the O(1S) 557.7 nm emission near the geographic equator for March/April, 1993. The airglow volume emission rate distribution is found to be strongly dependent on local time. Beginning at dusk, an intense airglow emission layer descends from a mean altitude of 95 km, reaching 89 km by midnight after which the emission rapidly decays. Shortly after midnight it reappears weakly at a higher altitude and remains at this level as the emission rate gradually increases towards dawn. This strong local time dependence leads us to conclude that the effect is tidally driven. Comparison with the Forbes (1982a,b) model suggest that total density perturbations and changes in the atomic oxygen mixing ratio may the cause of the changes in emission rate distribution between dusk and midnight. The reappearance of the emission after midnight may be caused by downward winds bringing oxygen-rich air from above.

Description: Numerical simulations of the Farley-Buneman instability in 2-1/2 dimensions using particle ions and fluid electrons show the growth, saturation and nonlinear behavior of two-stream waves. This hybrid technique models the saturated state of the instability for a much longer period of time than the pure particle codes that preceded it. While focusing principally on modeling the topside E region equatorial electrojet, many of these results apply to the auroral two-stream instability as well. The following features are seen in all our hybrid simulations: (1) wave growth at an angle offset from the electron drift direction where the angle depends on the strength of the driving electric field, (2) nonlinear coupling to waves traveling perpendicular to the propagation direction of the principal two-stream waves, (3) a saturated wave phase velocity at or above the sound speed but well below the velocity predicted by linear theory and (4) phase velocities which remain almost constant as a simulated radar sweeps from a horizontal direction to nearly vertical. The nonlinear electron motion dominates the behavior of these waves. Further, these simulations indicate that ion kinetic effects are not essential for the saturation of the instability and that electron temperature effects have a minor impact on the final saturated state.

Description: We have used a newly developed, parallelized, global MHD magnetosphere - ionosphere simulation model with a 400 R(sub E) long tail to study the evolution, structure, and dynamics of the distant magnetotail during extended periods of northward interplanetary magnetic field (IMF). We find that the tail evolves to a nearly time stationary structure about one solar wind transit time after the IMF turns northward. Four regions of different magnetic topology can be distinguished which extend at least to the end of the simulation box at 400 R(sub E). Besides lobe field lines and open solar wind field lines tailward of an X-line, there is a broad boundary layer of closed field lines which we call the tail flank boundary layer (TFBL). Just inside the TFBL there is a region of closed field loops. Besides the X-line we find two O-lines which are enclosed by the closed field loops and are roughly aligned with the tail axis. Together they form a U shaped separator between the northward and the southward plasma sheet fields.

Description: We demonstrate the use of Geostationary Operational Environmental Satellite (GOES) images of the May 18, 1980, Mount St. Helens volcanic plume in providing details of the dynamics and changing character of this major explosive eruption. Visible and thermal infrared (IR) data from a sequence of images at 30-min intervals from 0850 to 1720 Local Time (LT) give information on dispersal and plume top temperature. Initial visible and IR images at 0850 show the top of a spreading co-ignimbrite-like umbrella plume and an overshooting column emerging from it, both rising off the ground-hugging pyroclastic gravity flow generated by the opening directed blast. The overshooting column had a minimum temperature significantly colder than local ambient atmosphere, indicating substantial undercooling, and a maximum altitude of 31 +/- 2 km at 0920. This large plume system then formed a high-velocity, radially spreading, gravitationally driven current before becoming advected in the wind field at an average downwind velocity of 29 m/s. Reflectance values from visible GOES data change from lower to higher during periods of transition from darker toned Plinian to lighter toned co-ignimbrite plumes indicating that in this case satellite data resolved changes in eruptive style from plumes with a coarser to a finer dominant particle size.

Description: Multi-instrument data sets from the ground and satellites at both low and high altitude have provided new results concerning substorm onset and its source region in the magnetosphere. Twenty-six out of 37 substorm onset events showed evidence of azimuthally spaced auroral forms (AAFs) prior to the explosive poleward motion associated with optical substorm onset. AAFs can span 8 hours of local time prior to onset and generally propagate eastward in the morning sector. Onset itself is, however, more localized spanning only about 1 hour local time. AAF onset occur during time periods when the solar wind pressure is relatively high. AAFs brighten in conjunction with substorm onset leading to the conclusion that they are a growth phase activity casually related to substorm onset. Precursor activity associated with these AAFs is also seen near geosynchronous orbit altitude and examples show the relationship between the various instrumental definitions of substorm onset. The implied mode number (30 to 135) derived from this work is inconsistent with cavity mode resonances but is consistent with a modified flute/ballooning instability which requires azimuthal pressure gradients. The extended source region and the distance to the open-closed field line region constrain reconnection theory and local mechanisms for substorm onset. It is demonstrated that multiple onset substorms can exist for which localized dipolarizations and the Pi 2 occur simultaneously with tail stretching existing elsewhere. These pseudobreakups can be initiated by auroral streamers which originate at the most poleward set of arc systems and drift to the more equatorward main UV oval. Observations are presented of these AAFs in conjunction with low- and high-altitutde particle and magnetic field data. These place the activations at the interface between dipolar and taillike field lines probably near the peak in the cross-tail current. These onsets are put in the context of a new scenario for substorm morphology which employs individual modules which operate independently or couple together. This allows particular substorm events to be more accurately described and investigated.