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: 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: We present a new gap-free version of the seasonal and longitudinal 0 (s/l) variations of P(sub EFI), the equatorial F region irregularity (EFI) occurrence probability, based on data from the AE-E spacecraft. The agreement of this and three earlier partial P(sub EFI) patterns verifies all four. We reinterpret another earlier gap-ridden pattern, that of D(bar)(sub RSF), a topside ionogram index of average darkening by range spread F. We compare it with P(sub EFI) and, using ionosonde radio science considerations, we conclude that D(bar)(sub RSF) = P(sub EFI) times a factor depending on the average number of topside plasma bubbles visible to the ionosonde. The s/l variations of D(baar)(sub RSF) thus imply s/l variations in the average spacing of bubbles, whose seeds have an occurrence probability pattern P(sub seed). For discussion we assume P(sub EFI) = P(sub inst)P(sub seed) is the pattern of F region instability. The P(sub EFI) pattern, which is by definition independent of seed and/or bubble spacing, is far too complex to be explained by the dominant paradigm, that of changes in P(sub inst) by simple changes in the F region altitude and/or north-south asymmetry. We examine evidence behind this dominance, and find it unconvincing. Both the asymmetry and sunset-node/altitude hypotheses of 1984 and 1985, respectively, seem to be partly based on misunderstood data, and their features appear displaced in time and space from those of our repeatable P(sub EFI) pattern. In contrast, if P(sub seed) variations influence the P(sub EFI) pattern and depend on thermospheric gravity waves from tropospheric convection near the dip equator, then the seasonal maxima (minima) Of P(sub EFI) could be explained, since they all occur above relatively warm (cold) surface features, where convection is maximal (minimal). Also, the hypothesis of the dominance of the P(sub seed) term could explain an unusual December/January P(sub EFI) maximum in the deep, wide, normal Pacific minimum in the one data set obtained in El Nino years. Based on the experiments we consider, we predict that the s/l variations Of P(sub seed) will be found to be similar to those of P(sub EFI) and largely to explain them. Finally, we find reasons, based on the similarity of the D(sub RSF) variations to s/l patterns of the average scintillation index, for not using, as is commonly done, such scintillation patterns as substitutes for P(sub EFI) or P(sub inst) patterns.

Description: The concentrations of the hydrogen radicals OH and HO2 in the middle and upper troposphere were measured simultaneously with those of NO, O3, CO, H2O, CH4, non-methane hydrocarbons, and with the ultraviolet and visible radiation field. The data allow a direct examination of the processes that produce O3 in this region of the atmosphere. Comparison of the measured concentrations of OH and HO2 with calculations based on their production from water vapor, ozone, and methane demonstrate that these sources are insufficient to explain the observed radical concentrations in the upper troposphere. The photolysis of carbonyl and peroxide compounds transported to this region from the lower troposphere may provide the source of HO, required to sustain the measured abundances of these radical species. The mechanism by which NO affects the production Of O3 is also illustrated by the measurements. In the upper tropospheric air masses sampled, the production rate for ozone (determined from the measured concentrations of HO2 and NO) is calculated to be about I part per billion by volume each day. This production rate is faster than previously thought and implies that anthropogenic activities that add NO to the upper troposphere, such as biomass burning and aviation, will lead to production of more 03 than expected.

Description: Vertical profiles of infrared cirrus extinction have been derived from tropical and subtropical upper tropospheric solar occultation spectra. The measurements were recorded by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer during the Atmospheric Laboratory for Applications and Sciences (ATLAS) 3 shuttle flight in November 1994. The presence of large numbers of small ice crystals is inferred from the appearance of broad extinction features in the 8-12 micron region. These features were observed near the tropopause and at lower altitudes. Vertical profiles of the ice extinction (/km) in microwindows at 831, 957, and 1204/cm have been retrieved from the spectra and analyzed with a model for randomly oriented spheroidal ice crystals. An area-equivalent spherical radius of 6 gm is estimated from the smallest ice crystals observed in the 8-12 micron region. Direct penetration of clouds into the lower stratosphere is inferred from observations of cloud extinction extending from the upper troposphere to 50 mbar (20 km altitude). Cloud extinction between 3 and 5 micron shows very little wavelength dependence, at least for the cases observed by the ATMOS instrument in the tropics and subtropics during ATLAS 3.

Description: This study uses Global magnetohydrodynamic (MHD) simulations driven by solar wind data along with Geotail observations of the magnetotail to investigate the magnetotail's response to changes in the interplanetary magnetic field (IMF); observed events used in the study occurred on March 29, 1993 and February 9, 1995. For events from February 9, 1995, we also use the time-dependent MHD magnetic and electric fields and the large-scale kinetic (LSK) technique to examine changes in the Geotail ion velocity distributions. Our MHD simulation shows that on March 29, 1993, during a long period of steady northward IMF, the tail was strongly squeezed and twisted around the Sun-Earth axis in response to variations in the IMF B(sub y) component. The mixed (magnetotail and magnetosheath) plasma observed by Geotail results from the spacecraft's close proximity to the magnetopause and its frequent crossings of this boundary. In our second example (February 9, 1995) the IMF was also steady and northward, and in addition had a significant B(sub y) component. Again the magnetotail was twisted, but not as strongly as on March 29, 1993. The Geotail spacecraft, located approximately 30 R(sub E) downtail, observed highly structured ion distribution functions. Using the time-dependent LSK technique, we investigate the ion sources and acceleration mechanisms affecting the Geotail distribution functions during this interval. At 1325 UT most ions are found to enter the magnetosphere on the dusk side earthward of Geotail with a secondary source on the dawn side in the low latitude boundary layer (LLBL). A small percentage come from the ionosphere. By 1347 UT the majority of the ions come from the dawn side LLBL. The distribution functions measured during the later time interval are much warmer, mainly because particles reaching the spacecraft from the dawnside are affected by nonadiabatic scattering and acceleration in the neutral sheet.

Description: A method of inferring central plasma sheet (CPS) temperature, density, and pressure from ionospheric observations is developed. The advantage of this method over in situ measurements is that the CPS can be studied in its entirely, rather than only in fragments. As a result, for the first time, comprehensive two-dimensional equatorial maps of CPS pressure, density, and temperature within the isotropic plasma sheet are produced. These particle properties are calculated from data taken by the Special Sensor for Precipitating Particles, version 4 (SSJ4) particle instruments onboard DMSP F8, F9, F10, and F11 satellites during the entire year of 1992. Ion spectra occurring in conjunction with electron acceleration events are specifically excluded. Because of the variability of magnetotail stretching, the mapping to the plasma sheet is done using a modified Tsyganenko [1989] magnetic field model (T89) adjusted to agree with the actual magnetotail stretch at observation time. The latter is inferred with a high degree of accuracy (correlation coefficient -0.9) from the latitude of the DMSP b2i boundary (equivalent to the ion isotropy boundary). The results show that temperature, pressure, and density all exhibit dawn-dusk asymmetries unresolved with previous measurements. The ion temperature peaks near the midnight meridian. This peak, which has been associated with bursty bulk flow events, widens in the Y direction with increased activity. The temperature is higher at dusk than at dawn, and this asymmetry increases with decreasing distance from the Earth. In contrast, the density is higher at dawn than at dusk, and there appears to be a density enhancement in the low-latitude boundary layer regions which increases with decreasing magnetic activity. In the near-Earth regions, the pressure is higher at dusk than at dawn, but this asymmetry weakens with increasing distance from the Earth and may even reverse so that at distances X less than approx. 10 to -12 R(sub E), depending on magnetic activity, the dawn sector has slightly higher pressure. The temperature and density asymmetries in the near-Earth region are consistent with the ion westward gradient/curvature drift as the ions ExB convect earthward. When the solar wind dynamic pressure increases, CPS density and pressure appear to increase, but the temperature remains relatively constant. Comparison with previously published work indicates good agreement between the inferred pressure, temperature, and density and those obtained from in situ data. This new method should provide a continuous mechanism to monitor the pressure, temperature, and density in the magnetotail with unprecedented comprehensiveness.

Description: For over 20 years researchers have been investigating the feasibility of profiling tropospheric vector wind velocity from space with a pulsed Doppler lidar. Efforts have included theoretical development, system and mission studies, technology development, and ground-based and airborne measurements. Now NASA plans to take the next logical step towards enabling operational global tropospheric wind profiles by demonstrating horizontal wind measurements from the Space Shuttle in early 2001 using a coherent Doppler wind lidar system.