ALBERT

Description: The Grant was a three year grant funded under the Space Physics Supporting Research and Technology and Suborbital Program. Our objective was to develop automated techniques needed to unfold or "invert" global images of the magnetospheric ion populations obtained by the new magnetospheric imaging techniques (ENA, EUV) in anticipation of future missions such as the Magnetospheric Imager and, now, IMAGE. Our focus on the present three year grant is to determine the degree to which such images can quantitatively constrain the global electromagnetic properties of the magnetosphere. In a previous three year grant period we successfully automated a forward modeling inversion algorithm, demonstrated that these inversions are robust in the face of realistic instrumental considerations such as counting statistics and backgrounds, applied error analysis techniques to the extracted parameters using variational procedures, implemented very realistic magnetospheric test images to test the inversion algorithms using the Rice University Magnetospheric Specification Model, and began the process of generating parametric models with the flexibility to handle the realistic magnetospheric images (e.g. Roelof et al, 1992; 1993). Our plan for the present 3 year grant period was to complete the development of the inversion tools needed to handle realistic magnetospheric images, assess the degree to which global electrodynamics is quantitatively constrained by ENA images of the magnetosphere, and bring the inversion of EUV images up to the maturity that we will have achieved for the ENA imaging. Below the accomplishments of our three year effort are present followed by a list of our presentations and publications. The accomplishments of all three years are presented here, and thus some of these items appeared on interim progress reports.

Description: The Swedish microsatellite ASTRID was launched by a Russian Cosmos rocket on January 24, 1995 into a 1000 km circular orbit with 83 deg inclination. Besides the main objective of technological demonstration, imaging of energetic neutral atoms (ENAS) was attempted. The imager detected ENA in the energy range 0.1 - 140 keV utilizing two different techniques. Neutrals of the energy 13 - 140 keV were recorded by 14 solid state detectors with the total field of view 5 deg x 322 deg. For half a spin (approx. 1.5 s) of the ASTRID spacecraft, almost all of space was covered with an angular resolution 2.5 deg x 25 deg. Less energetic neutrals of approx. 0.1 - 70 keV were converted on a graphite target into secondary particles which then were detected by a microchannel plate with 32 anodes. A fraction of primary neutrals was directly reflected towards the sensor. This technique provided the total ENA flux with an angular resolution 4.6 deg x 11.5 deg. The instrument weight is 3.13 kg. Successful operation of the instrument during the first 5 weeks of the mission provided the first ENA images of the ring current at low altitudes.

Description: Two aspects of the cloud ice parameterization in the Goddard Cumulus Ensemble Model cloud physics parameterization are examined: the conversion of cloud ice to snow by depositional growth, designated PSFI, and the saturation adjustment scheme. The original formulation of PSFI is shown to produce excessive conversion of cloud ice to snow because of an implicit assumption that the relative humidity is 100% with respect to water even though the air may actually be quite less humid. Two possible corrections to this problem are proposed, the first involving application of a relative humidity dependent correction factor to the original formulation of PSFI, and the second involving a new formulation of PSFI based on the equation for depositional growth of cloud ice. The sensitivity of these formulations of PSFI to the assumed masses of the ice particles is examined. Possible problems associated with using a saturation adjustment scheme for cloud ice are discussed and simulations of a squall line with and without application of the adjustment scheme for ice are compared.

Description: On May 23, 1995, the Comprehensive Plasma Instrumentation (CPI) onboard the Geotail spacecraft observed a complex and structured ion distribution function near the magnetotail midplane at x approximately -10 R(sub E). On the same day, the Wind spacecraft observed a very high density (approximately 40/cubic cm) solar wind and an interplanetary magnetic field (IMF) that was predominantly northward but had several southward turnings. We have inferred the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time using time-dependent electric and magnetic fields obtained from a global MHD (magnetohydrodynamic) simulation. Wind data were used as input for the MHD model. We found that three sources contributed to this distribution: the ionosphere, the plasma mantle which had near-Earth and distant tail components, and the low latitude boundary layer (LLBL). Moreover, distinct structures in the low energy part of the distribution function were found to be associated with individual sources. Structures near 0 deg pitch angle were made up of either ionospheric or plasma mantle ions, while structures near 90 deg pitch angle were dominated by ions from the LLBL source. Particles that underwent nonadiabatic acceleration were numerous in the higher energy part of the ion distribution function, whereas ionospheric and LLBL ions were mostly adiabatic. A large proportion of the near-Earth mantle ions underwent adiabatic acceleration, while most of the distant mantle ions experienced nonadiabatic acceleration.

Description: Ion density and velocity measurements from the Dynamics Explorer 2 (DE 2) spacecraft are used to obtain the average magnetic local time versus invariant latitude distribution of irregularities in the high-latitude F region ionosphere. To study the small-scale structure and its relationship to background conditions in the ionosphere, we have formed a reduced database using 2-s (approx. = 16 km) segments of the ion density and velocity data. The background gradients associated with each 2-s segment and the spectral characteristics, such as power at 6 Hz (approx. = 1.3 km) and spectral index, are among the reduced parameters used in this study. The relationship between the observed plasma structure and its motion is complex and dependent on the externally applied fields as well as locally generated plasma structure. The evolution of plasma structures also depends critically on the conductivity of the underlying ionosphere. Observations indicate an enhancement of irregularity amplitudes in two spatially isolated regions in both the ion density and the velocity. Convective properties seem to play a more important role in winter hemisphere where smaller-scale structures are maintained outside the source regions. (Delta)V irregularity amplitudes are enhanced in the cusp and the polar cap during northward interplanetary magnetic field regardless of season. The power in (Delta)V is usually higher than that associated with local polarization electric fields, suggesting that the observed structure in (Delta)N/N is strongly influenced by (Delta)V structure applied to large density gradients.

Description: Analysis of meteorological, chemical and microphysical data from the airborne SUCCESS (SUbsonic aircraft Contrail and Cloud Effects Special Study) mission is reported. Careful analysis of the complex DC-8 flight pattern of May 2, 1996 reveals 19 linear flight segments within six main geographical areas, which we have analyzed. Significant mountain wave activity is revealed in the data from the MMS (Meteorology Measurement System) and MTP (Microwave Temperature Profiler) instruments on the DC-8, which resembles previous observations of mountain wave structures near Boulder, Colorado. Strong mountain-wave-induced upwelling downwind of the Rockies is noted. Turbulence is also noted in regions of the mountain wave consistent with overturning near the tropopause. Zonal winds recorded on the ER-2 are shown to be consistent with mountain wave breaking at or near critical levels in the stratosphere, consistent with the strong turbulence reported by the pilot during the ER-2 flight. These observations have been supported with spectral analyses and modeling studies. 'Postcasts' of mountain wave activity on May 2, 1996 using the Naval Research Laboratory Mountain Wave Forecast Model predicts both strong mountain wave activity near the tropopause and strong mountain-wave-induced turbulence in the stratosphere.

Description: Observations of energetic neutral atoms (ENA) in the energy range 26- 52 keV are reported from four occasions during geomagnetically disturbed periods. The data were acquired by the ENA imager flown on the Swedish microsatellite Astrid in a 1000 km circular orbit with 83 deg inclination. The ENA imager separates charged particles from neutrals through an electrostatic deflection system in the energy range between 0.1 and 114 keV. ENA images obtained from vantage points in the polar cap and in the afternoon magnetic local time (MLT) hours looking into the antisunward hemisphere show intense ENA fluxes (approx. 10(exp 4)/sq cm sr s over 26-37 keV) coming from the dusk region and low altitudes (approx. 300 km). The morphology shows no relation to local magnetic field excluding the possibility of charged particle detection. It is concluded that the source of these ENAs are precipitating/mirroring ions from the ring current/trapped radiation interacting with the exobase on auroral L-shells and in the dusk region. The observed ENA fluxes show a relation with Kp and Dst geomagnetic indices. The observed ENA spectrum from a geomagnetic storm on February 8, 1995, is investigated in more detail and compared to the parent ion spectrum obtained by the Defense Meteorological Satellite Project (DMSP) satellite, Fl2, during the same period on L = 6 +/- 2 around dusk. The observed ENA spectral slope is used to derive the parent ion spectral temperature. The derived ion temperatures range is 3.0 - 6.0 keV for H and 4.5 - 8.5 keV for O. The higher of these ion temperatures comes closest in agreement to the extrapolated DMSP spectrum leading us to favor O over H as the species of the detected ENAS. It is shown that the detected ENAs must have been produced at L greater than or equal to 6 to reach the detector without atmospheric attenuation and that the main energy dependence of the ENA spectrum, apart from the parent ion spectrum, is governed by the energy dependence of the charge exchange cross section between ions and exospheric oxygen.

Description: The purpose of this project was to determine the suitability of measuring active deformation of volcanoes in Alaska using Interferometric Synthetic Aperture Radar (INSAR) techniques. Work sponsored by this grant supported one graduate student (for almost 2 years) and one postdoc (for several months), and has resulted in two published peer-reviewed papers and a front-page article in EOS. An additional paper is in review and a fourth is in preparation. An additional paper in preparation was based in part on research supported by this grant and in part by a successor grant from NASA's Solid Earth Natural Hazards program. Over the course of this research, we documented measurable uplift of Trident volcano in the Katmai group, conducted a systematic study of the change in phase coherence over time on volcanic surfaces, and measured and modeled the spectacular 1.5 m deflation of Okmok caldera associated with its 1997 eruption. We also generated initial interferograms spanning the 1996 seismic swarm of Akutan volcano; however, during the period covered by this project we were not able to remove topography. That has been done under the subsequent funding and a paper is now in preparation. This report summarizes work done under two separate contracts because both were based on the same proposal to NASA's ADRO (Application Development and Research Opportunity) program. The first year was funded out of a grant from NASA Headquarters and the second and third years out of a grant through Goddard. The work, however, was a continuous three year effort.

Description: This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model in which the ion current periodically updates the magnetic field. The plasma mantle supplies particles continuously to the magnetotail, and the perturbation magnetic field is calculated from the ion current using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the near-Earth X-line in the model. This variability is caused by the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. Particularly noteworthy is the value found for the characteristic time scale of variability in the magnetotail. on the order of 4 - 5 minutes.

Description: This study investigates the evolution of the magnetotail's magnetic field with the aid of a self-consistent two-dimensional model. In this model the plasma mantle continuously supplies particles to the magnetotail, the ion current periodically updates the magnetic field using the Biot-Savart law. The simulated magnetotail evolves into a quasi-steady state, characterized by the periodic motion of the model's near-Earth X-line. This variability results from the nonadiabatic acceleration of ions in the current sheet and their rapid loss from the tail. The characteristic time scale of variability in the magnetotail is on the order of 4 - 5 minutes. We also investigate how the magnetotail's topology responds to increased convection electric fields, and show examples of observations of variability in the magnetotail.