ALBERT

Description: Our studies elucidated the relationship between the auroral arcs and magnetotail phenomena. One paper examined particle energization in the source region of the field-aligned currents that intensify at substorm onset when the arc brightens to form the westward electrojet. A second paper examined the relationship between the precipitating particles in the arcs, the location of the westward electrojet, and magnetospheric source regions. Two earlier papers also investigated the roles that field aligned currents and particle acceleration have during substorms.

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: This final report describes the work done by Dr. Marc Hairston and Dr. Rod Heelis on NASA SR&T grant NAGW-4411 studying the theta aurora using DE-1 ultraviolet imager data and DMSP particle data. This report covers the period from summer 1995 through summer 1996 along with a review of the previous work. Previous work on this grant looked at the time period from the launch of DMSP-F8 in June 1987 through the end of mission of DE-1 in summer 1991. Despite the sporadic and decreasing frequency of observations from DE-1 over this time period, we were able to identify six events for study where the DE-1 imager observed a theta aurora occurring during a period in which the DMSP-F8 satellite flew through the region of the aurora at an altitude of 800 km. We focused on the best two events where the theta aurora persisted for an extended period so that we could observe the DMSP particle signatures in both the hemisphere observed by DE-1 and in the other polar hemisphere immediately before or afterwards. These results were presented at the Fall 1994 meeting of the AGU. Initially we had hoped to expand on this work. However, further work showed that none of these events gave a clear enough signature in the DMSP data for us to identify a theta aurora in the hemisphere opposite to the hemisphere imaged by DE-1. Without that, there were no results from this work that were new enough to warrant publication. So instead we used the final year of the grant to work with our colleague, J. A. Cumnock on a similar project using DE data to study the evolution of theta auroras as a function of the IMF A paper from that work was published which acknowledged this grant and a copy of that paper is included with this final report.

Description: The accuracy of a new local gravity field model, GEOID94A, is examined at a site on the western Greenland ice sheet. The model, developed by the Danish National Survey and Cadastre, incorporates several new gravity data sets including an extensive amount of airborne gravity data. Model-derived geoid undulations were compared to independently determined undulations found by differencing the elevations from Global Positioning System controlled airborne laser altimetry and optical leveling surveys. Differences between the two sets of undulations were less than +/- 6 cm RMS. The comparison improved (+/- 5 cm RMS) when GEOID94A undulations were adjusted by local gravity observations also acquired at the site. Our comparisons demonstrate that GEOID94A adequately models the long to intermediate wavelengths of the gravity field. We conclude that GEOID94A constitutes a reliable reference model for studies of Greenland's gravity field.

Description: Under the mandate contained in the FY 1976 NASA Authorization Act, the National Aeronautics and Space Administration (NASA) has developed and is implementing a comprehensive program of research, technology, and monitoring of the Earth's upper atmosphere, with emphasis on the stratosphere. This program aims at expanding our understanding to permit both the quantitative analysis of current perturbations as well as the assessment of possible future changes in this important region of our environment. It is carried out jointly by the Upper Atmosphere Research Program (UARP) and the Atmospheric Chemistry Modeling and Analysis Program (ACMAP), both managed within the Science Division in the Office of Mission to Planet Earth at NASA. Significant contributions to this effort are also provided by the Atmospheric Effects of Aviation Project (AEAP) of NASA's Office of Aeronautics. The long-term objectives of the present program are to perform research to: understand the physics, chemistry, and transport processes of the upper atmosphere and their effect on the distribution of chemical species in the stratosphere, such as ozone; understand the relationship of the trace constituent composition of the lower stratosphere and the lower troposphere to the radiative balance and temperature distribution of the Earth's atmosphere; and accurately assess possible perturbations of the upper atmosphere caused by human activities as well as by natural phenomena. In compliance with the Clean Air Act Amendments of 1990, Public Law 101-549, NASA has prepared a report on the state of our knowledge of the Earth's upper atmosphere, particularly the stratosphere, and on the progress of UARP and ACMAP. The report for the year 1996 is composed of two parts. Part 1 summarizes the objectives, status, and accomplishments of the research tasks supported under NASA UARP and ACMAP in a document entitled, Research Summary 1994-1996. Part 2 is entitled Present State of Knowledge of the Upper Atmosphere 1996.- An Assessment Report. It consists primarily of the Executive Summary and Chapter Summaries of the World Meteorological Organization Global Ozone Research and Monitoring Project Report No. 37, Scientific Assessment of Ozone Depletion: 1994, sponsored by NASA, the National Oceanic and Atmospheric Administration (NOAA), the UK Department of the Environment, the United Nations Environment Program, and the World Meteorological Organization. Other sections of Part 11 include summaries of the following: an Atmospheric Ozone Research Plan from NASA's Office of Mission to Planet Earth; summaries from a series of Space Shuttle-based missions and two recent airborne measurement campaigns; the Executive Summary of the 1995 Scientific Assessment of the Atmospheric Effects of Stratospheric Aircraft, and the most recent evaluation of photochemical and chemical kinetics data (Evaluation No. 12 of the NASA Panel for Data Evaluation) used as input parameters for atmospheric models.

Description: Horizontal wavenumber power spectra of vertical and horizontal wind velocities, potential temperatures, and ozone and N(2)O mixing ratios, as measured in the mid-stratosphere during 73 ER-2 flights (altitude approx. 20km) are presented. The velocity and potential temperature spectra in the 100 to 1-km wavelength range deviate significantly from the uniform -5/3 power law expected for the inverse energy-cascade regime of two-dimensional turbulence and also for inertial-range, three-dimensional turbulence. Instead, steeper spectra approximately consistent with a -3 power law are observed at horizontal scales smaller than 3 km for all velocity components as well as potential temperature. Shallower spectra are observed at scales longer than 6 km. For horizontal velocity and potential temperature the spectral indices at longer scales are between -1.5 and -2.0. For vertical velocity the spectrum at longer scales become flat. It is argued that the observed velocity and potential temperature spectra are consistent with gravity waves. At smaller scales, the shapes are also superficially consistent with a Lumley-Shur-Weinstock buoyant subrange of turbulence and/or nonlinear gravity waves. Contemporaneous spectra of ozone and N(sub 2)O mixing ratio in the 100 to 1-km wavelength range do conform to an approximately uniform -5/3 power law. It is argued that this may reflect interactions between gravity wave air-parcel displacements and laminar or filamentary structures in the trace gas mixing ratio field produced by enstropy-cascading two-dimensional turbulence.

Description: Cryptoblemes are subtle impact shock signatures imprinted by cosmic debris on the crustal surfaces of lunar planetary bodes. These signatures constitute a complex cumulative overprinting of topographic, structural geophysical, and tectonic patterns that have a conspicuous radial centric multiringed symmetry. The geometry and distribution of cryptoblemes on Earth is comparable to the size and density of impact features on lunar planetary surfaces. Analysis of satellite imagery, sea-floor sonar, side-looking radar and aerial photographs of specific sites reveals new criteria for the identification and confirmation of impact-shock signatures. These criteria include joint and foliation patterns with asbestiform minerals, ribbon-quartz, spheroidal weathering, domal exfoliation, pencil shale, and shock spheres, which may originate from hydrocavitation of water-saturated sedimentary rocks. Cryptoblemes may also be associated with breccia pipes, sinkholes, buttes, mesas, and bogs, high-Rn anomalies, nodular concentrations, and earthquake epicenters. Major implications of cryptobleme identification include exploratory targeting of hydrocarbon and mineral deposits and the explanation of their origins. Analysis of known mineral deposits, structural traps and sedimentary basins show a direct correlation with cryptobleme patterns. Significant geologic paradigm shifts related to cryptoblemes include mountain building processes, structural orogenies, induced volcanism, earthquake origins, hydrocarbon diagenesis, formation mineral deposits, continental rifting, and plate movements, magnetic overprinting and local regional, and global geologic extinction and speciation patterns. Two figures provide a comparison between a multiring impact overprint in water and multiring cryptobleme in the U.S. basin range. (Additional information is contained in the original document).

Description: In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Description: Twenty minerals that were not included in the most recent list of meteoritic minerals have been reported as occurring in meteorites. Extraterrestrial anhydrous Ca phosphate should be called menillite, not whitlockite.

Description: Galim is a polymict breccia consisting of a heavily shocked (shock stage S6) LL6 chondrite, Galim (a), and an impact-melted EH chondrite, Galim (b). Relict chondrules in Galim (b) served as nucleation sites for euhedral enstatite grains crystallizing from the impact melt. Many of the reduced phases typical of EH chondrites (e.g., Si-bearing metallic Fe-Ni; Ti-bearing troilite) are absent. Galim (b) was probably shock-melted while in contact with a more oxidized source, namely, Galim (a); during this event, Si was oxidized from the metal and Ti was oxidized from troilite. Galim (a) contains shock veins and recrystallized, unzoned olivine. The absence of evidence for reduction in Galim (a) may indicate that the amount of LL material greatly exceeded that of EH material; shock metamorphism may have taken place on the LL parent body. Shock-induced redox reactions such as those inferred for the Galim breccia appear to be restricted mainly to asteroids because the low-end tail of their relative-velocity distribution permits mixing of intact disparate materials (including accretion of projectiles of different oxidation states), whereas the peak of the distribution leads to high equilibration shock pressures (allowing impact-induced exchange between previously accreted, disequilibrated materials). Galim probably formed by a two-stage process: (I) accretion to the LL parent body of an intact EH projectile at low relative velocities, and (2) shock metamorphism of the assemblage by the subsequent impact of another projectile at significantly higher relative velocities.

Description: Approximately 275 mineral species have been identified in meteorites, reflecting diverse redox environments, and, in some cases, unusual nebular formation conditions. Anhydrous ordinary, carbonaceous and R chondrites contain major olivine, pyroxene and plagioclase; major opaque phases include metallic Fe-Ni, troilite and chromite. Primitive achondrites are mineralogically similar. The highly reduced enstatite chondrites and achondrites contain major enstatite, plagioclase, free silica and kamacite as well as nitrides, a silicide and Ca-, Mg-, Mn-, Na-, Cr-, K- and Ti-rich sulfides. Aqueously altered carbonaceous chondrites contain major amounts of hydrous phyllosilicates, complex organic compounds, magnetite, various sulfates and sulfides, and carbonates. In addition to kamacite and taenite, iron meteorites contain carbides, elemental C, nitrides, phosphates, phosphides, chromite and sulfides. Silicate inclusions in IAB/IIICD and lIE iron meteorites consist of mafic silicates, plagioclase and various sulfides, oxides and phosphates. Eucrites, howardites and diogenites have basaltic to orthopyroxenitic compositions and consist of major pyroxene and calcic plagioclase and several accessory oxides. Ureilttes .are made up mainly of calcic, chromian olivine and low-Ca clinopyroxene embedded in a carbonaceous matrix; accessory phases include the C polymorphs graphite, diamond, lonsdaleite and chaoite as well as metallic Fe-Ni, troilite and halides. Angrites are achondrites rich in fassaitic pyroxene (i.e. , AI-Ti diopside); minor olivine with included magnesian kirschsteinite is also present. Martian meteorites comprise basalts, Iherzolites, a dunite and an orthopyroxenite. Major phases include various pyroxenes and olivine; minor to accessory phases include various sulfides, magnetite, chromite and Ca-phosphates. Lunar meteorites comprise mare basalts with major augite and calcic plagioclase and anorthositic breccias with major calcic plagioclase. Several meteoritic phases were formed by shock metamorphism. Martensite (a2-fe,Ni) has a distorted body-centered-cubic structure and formed by a shear transformation from taenite during shock reheating and rapid cooling. The C polymorphs diamond, lonsdaleite and chaoite formed by shock from graphite. Suessite formed in the North Haig ureilite by reduction of Fe and Si (possibly from olivine) via reaction with carbonaceous matrix material. Ringwoodite, the spinel form of (Mg,Fe)2Si04, and majorite, a polymorph of (Mg,Fe)Si03 with the garnet structure, formed inside shock veins in highly shocked ordinary chondrites. Secondary minerals in meteorite finds that formed during terrestrial weathering include oxides and hy-. droxides formed directly from metallic Fe-Ni by oxidation, phosphates formed by the alteration of schreibersite, and sulfates formed by alteration of troilite.

Description: Igneous graphite. a rare constituent in terrestrial mafic and ultramafic rocks. occurs in three EH and one EL enstatite chondrite impact-melt breccias as 2-150 Ilm long euhedrallaths. some with pyramidal terminations. In contrast. graphite in most enstatite chondrites exsolved from metallic Fe-Ni as polygonal. rounded or irregular aggregates. Literature data for five EH chondrites on C combusting at high temperatures show that Abee contains the most homogeneous C isotopes (i.e. delta(sup 13)C = -8.1+/-2.1%); in addition. Abee's mean delta(sup l3)C value is the same as the average high-temperature C value for the set of five EH chondrites. This suggests that Abee scavenged C from a plurality of sources on its parent body and homogenized the C during a large-scale melting event. Whereas igneous graphite in terrestrial rocks typically forms at relatively high pressure and only moderately low oxygen fugacity (e.g., approx. 5 kbar. logfO2, approx. -10 at 1200 C ). igneous graphite in asteroidal meteorites formed at much lower pressures and oxygen fugacities.

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: The power of mineralogical analysis as a descriptive or predictive technique stems from the fact that only a few thousand minerals are known to occur in nature as compared to several hundred thousand inorganic compounds. Further, all of the known minerals have specific stability ranges in pressure, temperature, an composition. A specific knowledge of the mineralogy of a planets surface or interior therefore allows one to characterize the present or past conditions under which the minerals were formed or have existed. For the purposes of this paper, a slightly broader definition of mineralogy was adopted by including not only crystalline materials found on planetary surfaces, but also ices and classes that can benefit from in situ types of analyses. Both visual examination and the various spectroscopies available for robotic probes to planetary surfaces are discussed.

Description: A conceptual design study of the ARIM-1 sounding rocket mission, whose goal is to study atmospheric turbulence in the tropopause region of the atmosphere, is presented. The study was conducted by an interdisciplinary team of students at the University of Alaska Fairbanks who were enrolled in a Space Systems Engineering course. The implementation of the ARIM-1 mission will be carried out by students participating in the Alaska Student Rocket Program (ASRP), with a projected launch date of August 1997. The ARIM-1 vehicle is a single stage sounding rocket with a 3:1 ogive nose cone, a payload diameter of 8 in., a motor diameter of 7.6 in., and an overall height of 17.0 ft including the four fins. Emphasis is placed on standardization of payload support systems. The thermosonde payload will measure the atmospheric turbulence by direct measurement of the temperature difference over a distance of one meter using two 3.45-micron 'hot-wire' probes. The recovery system consists of a 6 ft. diameter ribless guide surface drogue chute and a 33 ft. diameter main cross parachute designed to recover a payload of 31 pounds and slow its descent rate to 5 m/s through an altitude of 15 km. This document discusses the science objectives, mission analysis, payload mechanical configuration and structural design, recovery system, payload electronics, ground station, testing plans, and mission implementation.

Description: We present a ray-tracing study of the propagation of Pc 3 Alfven mode waves originating at the dayside magnetopause. This study reveals interesting features of magnetospheric filter effect for these waves. Pc 3 Alfven mode waves cannot penetrate to low Earth altitudes unless the wave frequency is below approximately 30 mHz. Configurations of the dispersion curves and the refractive index show that the gyroresonance and pseudo-cutoff introduced by the heavy ion O(+) block the waves. When the O(+) concentration is removed from the plasma composition, the barriers caused by the O(+) no longer exist, and waves with much higher frequencies than 30 mHz can penetrate to low altitudes. The result that the 30 mHz or lower frequency Alfven waves can be guided to low altitudes agrees with ground-based power spectrum observation at high altitudes.

Description: The Hatizyo Hydrographic Observatory, which is one of the essential magnetic observatories in Japan, was established in 1979 and is currently operated by the Hydrographic Department, Maritime Safety Agency. This is the annual report compiled from the results of magnetic observations carried out at the observatory in 1993. As to the instruments used for magnetic observations, the digital recording variometer was replaced by a fluxgate magnetometer in 1986, and one set each of the proton and fluxgate magnetometers was additionally installed in January and October 1992, respectively.

Description: Ion outflow from the ionosphere plays a fundamental but poorly defined role in magnetospheric processes. The purpose of the research is to better understand the mass coupling between the Earth's ionosphere and Magnetosphere. The work performed under this grant falls in three areas: (1) event studies using archived data from the DE-1/2 satellites; (2) investigations using Data from the ISTP satellites; and (3) work supporting a Space Physics Educational Outreach (SPEO) grant supplement.

Description: The purpose of this project was to quantify the rates of two processes which are crucial to our understanding of radiative energy balance in the upper atmosphere. The first process is radiative emission from vibrationally hot OH radicals following the H + O3 reaction in the upper mesosphere. The importance of this process depends strongly on the OH radiative emission coefficients. Our goal was to measure the OH permanent dipole moment in excited vibrational states and to use these measurements to construct an improved OH dipole moment function and improved radiative emission coefficients. Significant progress was made on these experiments including the construction of a supersonic jet source for vibrationally excited OH radicals. Unfortunately, our efforts to transport the OH radicals into a second lower pressure vacuum chamber were not successful, and we were unable to make improved dipole moment measurements for OH. The second key kinetic process which we attempted to quantify during this project is the rate of relaxation of bend-excited CO2 by oxygen atoms. Since excitation of the bending vibrational mode of CO2 is the major cooling mechanism in the upper mesosphere/lower thermosphere, the cooling rate of this region depends crucially on the rate of energy transfer out of this state. It is believed that the most efficient transfer mechanism is via atomic oxygen but the rate for this process has not been directly measured in the laboratory at appropriate temperatures and even the room temperature rate remains controversial. We attempted to directly measure the relaxation rate Of CO2 (010) by oxygen atoms using the discharge flow technique. This experiment was set up at Aerodyne Research. Again, significant progress was achieved in this experiment. A hot CO2 source was set up, bend excited CO2 was detected and the rate of relaxation of bend excited CO2 by He atoms was measured. Unfortunately, the project ran out of time before the oxygen atom kinetic studies could be implemented.

Description: There is an increasing concern of land owners to protect and maintain healthy and sustainable agroecosystems through the implementation of best management practices (BMP). The objectives of this study were: (1) To develop and evaluate the use of a Geographic Information System (GIS) technology for enhancing field-scale management practices; (2) evaluate the use of 2-dimensional displays of the landscape and (3) define spatial classes of variables from interpretation of geostatistical parameters. Soil samples were collected to a depth of 2 m at 15 cm increments. Existing data from topographic, land use, and soil survey maps of the Winfred Thomas Agricultural Research Station were converted to digital format. Additional soils data which included texture, pH, and organic matter were also generated. The digitized parameters were used to create a multilayered field-scale GIS. Two dimensional (2-D) displays of the parameters were generated using the ARC/INFO software. The spatial distribution of the parameters evaluated in both fields were similar which could be attributed to the similarity in vegetation and surface elevation. The ratio of the nugget to total semivariance, expressed as a percentage, was used to assess the degree of spatial variability. The results indicated that most of the parameters were moderate spatially dependent Biophysical constraint maps were generated from the database layers, and used in multiple combination to visualize results of the BMP. Understanding the spatial relationships of physical and chemical parameters that exists within a field should enable land managers to more effectively implement BMP to ensure a safe and sustainable environment.

Description: We are investigating the use of nitron as a potential chemical sensor for nitric acid and other electron deficient nitrogen oxides. Solutions of nitron in 1-propanol, toluene, and chloroform have been tested for use on a piezoelectric quartz crystal microbalance. We are testing various solvents and metal cations which can maximize the lifetime and reaction specificity of nitron so that they may be used as chemical coatings for stratospheric measurement of trace gases. Results of the work to date will be shown, and future direction discussed.

Description: There has been considerable recent interest in the relationship between the cross-tail magnetic field component B(sub y) and tail dynamics. The purpose of this paper is to give an overall description of the penetration of the interplanetary magnetic field (IMF) B(sub y) into the near-Earth plasma sheet. We show that plasma sheet B(sub y) may be generated by the differential shear motion of field lines and enhanced by flux tube compression. The latter mechanism leads to a B(sub y) analogue of the pressure-balance inconsistency as flux tubes move from the far tail toward the Earth. The growth of B(sub y), however, may be limited by the dawn-dusk asymmetry in the shear velocity as a result of plasma sheet tilting. B(sub y) penetration into the plasma sheet implies field-aligned currents flowing between hemispheres. These currents together with the IMF B(sub y) related mantle field-aligned currents effectively shield the lobe from the IMF B(sub y).

Description: Distinctive wave forms in the distributions of vertical velocity and temperature of both neutral particles and ions are frequently observed from Dynamics Explorer 2 at altitudes above 250 km over the polar caps. These are interpreted as being due to internal gravity waves propagating in the neutral atmosphere. The disturbances characterized by vertical velocity perturbations of the order of 100 m/s and horizontal wave lengths along the satellite path of about 500 km. They often extend across the entire polar cap. The associated temperature perturbations indicate that the horizontal phase progression is from the nightside to the dayside. Vertical displacements are inferred to be of the order of 10 km and the periods to be of the order of 10(exp 3) s. The waves must propagate in the neutral atmosphere, but they usually are most clearly recognizable in the observations of ion vertical velocity and ion temperature. By combining the neutral pressure calculated from the observed neutral concentration and temperature with the vertical component of the neutral velocity, an upward energy flux of the order of 0.04 erg/sq cm-s at 250 km has been calculated, which is about equal to the maximum total solar ultraviolet heat input above that altitude. Upward energy fluxes calculated from observations on orbital passes at altitudes from 250 to 560 km indicate relatively little attenuation with altitude.

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 geomagnetic field at a given epoch can be used to partition the surface of the liquid outer core into a finite number of contiguous regions in which the radial component of the magnetic flux density, B (sub r), is of one sign. These flux patches are instrumental in providing detail to surface fluid flows inferred from the changing geomagnetic field and in evaluating the validity of the frozen-flux approximation on which such inferences rely. Most of the flux patches in models of the modem field are small and enclose little flux compared to the total unsigned flux emanating from the core. To demonstrate that such patches are not required to explain the most spatially complete and accurate data presently available, those from the Magsat mission, I have constructed a smooth core field model that fits the Magsat data but does not possess small flux patches. I conclude that our present knowledge of the geomagnetic field does not allow us to resolve these features reliably at the core-mantle boundary; thus we possess less information about core flow than previously believed.

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: In this paper we derive the average configuration of the ring current as a function of the state of the magnetosphere as indicated by the Dst index. We sort magnetic field data from the Combined Release and Radiation Effects Satellite (CRRES) by spatial location and by the Dst index in order to produce magnetic field maps. From these maps we calculate local current systems by taking the curl of the magnetic field. NN7e find both the westward (outer) and the eastward (inner) components of the ring current. We find that the ring current intensity varies linearly with D.St as expected, and that the ring current is asymmetric for all Dst values. The azimuthal peak of the ring current is located in the afternoon sector for quiet conditions, and near midnight for disturbed conditions. The ring current also moves closer to the Earth during disturbed conditions. We are able to recreate the Dst index by integrating the magnetic perturbations caused by the ring current. We find that we needed to apply a 20 nT offset to Dst, and assume a perfectly conducting Earth to obtain an optimal agreement between the computed and the observed Dst. We interpret the 20 nT offset as the magnetic field generated by the quiet time ring current used as baseline in computing Dst.

Description: This project resulted in the first 2-D maps of magnetotail pressure, density and temperature. The results were published in JGR. A copy of this paper is attached. Also a magnetotail viewer was developed to allow the user to examine magnetotail plasma from different vantages. We hope to have this viewer online soon (at our web site http://sd-www.jhuapl.edu/Aurora).

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: Dimethylsulfide (DMS) is important in influencing the formation of aerosols in the troposphere over large areas of the world's oceans. Understanding the dynamics of aerosols is important to understanding the earth's radiation balance. In evaluating the factors controlling DMS in the troposphere it is vital to understand the dynamics of DMS in the surface ocean. The biogeochemical processes controlling DMS concentration in seawater are myriad; modeling and theoretical estimation are problematic. At the beginning of this project we believed that we were on the verge of simplifying the ship-track measurement of DMS, and we proposed to deploy such a system to develop a database relating high frequency DMS measurements to biological and physicochemical and optical properties of surface water that can be quantified by remote sensing techniques. We designed a system to measure DMS concomitantly with other basic chemical and biological data in a flow-through system. The project was collaborative between Woods Hole Oceanographic Institution (WHOI) and Bermuda Biological Station for Research (BBSR). The project on which we are reporting was budgeted for only one year with a one year no-cost extension. At WHOI our effort was directed towards designing traps which would be used to concentrate DMS from seawater and allow storage for subsequent analysis. At that time, GC systems were too large for easy long-term deployment on a research vessel like R/V Weatherbird, so we focused on simplifying the shipboard sampling procedure. Initial studies of sample recovery with high levels of DMS suggested that Carboxen 1000, a relatively new carbon molecular sieve, could be used as a stable storage medium. The affinity of Carboxen for DMS is several orders of magnitude higher than gold wool (another adsorbent used for DMS collection) on a weight or volume basis. Furthermore, Carboxen's affinity for DMS is also far less susceptible to humidity than gold wool. Unfortunately, further experiments with low level DMS indicated that recovery of DMS after storage was not quantitative. The material has proven to be completely acceptable for short term storage and has been incorporated into a micro-GC system. Since working on this project, we have collaborated with RVM Scientific in Santa Barbara in the design and construction of small portable micro-GC's that will make feasible at-sea measurement in moving ships, making rapid gas analysis and quantification feasible in a ship-track mode. Throughout this period at both WHOI and BBSR, we continued to analyze field data to understand that patterns of time and space variability in DMS and the processes that govern it. These insights will be crucial to determining the specifications for our automated sampling program. The data from this, the longest continuous sampling program for ocean DMS, provided insights into year to year and short-term variability.

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.

Description: This study employs Geotail plasma observations and numerical modeling to determine sources of the ions observed in the near-Earth magnetotail near midnight during a substorm. The growth phase has the low-latitude boundary layer as its most important source of ions at Geotail, but during the expansion phase the plasma mantle is dominant. The mantle distribution shows evidence of two distinct entry mechanisms: entry through a high latitude reconnection region resulting in an accelerated component, and entry through open field lines traditionally identified with the mantle source. The two entry mechanisms are separated in time, with the high-latitude reconnection region disappearing prior to substorm onset.

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 transmission electron microscope study of experimental shock metamorphism in natural pre-graphitic carbon simulates the response of the most common natural carbons to increased shock pressure. The d-spacings of this carbon are insensitive to the shock pressure and have no apparent diagnostic value, but progressive comminution occurs in response to increased shock pressure up to 59.6 GPa. The function, P = 869.1 x (size(sub minimum )(exp -0.83), describes the relationship between the minimum root-mean-square subgrain size (nm) and shock pressure (GPa). While a subgrain texture of natural pregraphitic carbons carries little information when pre-shock textures are unknown, this texture may go unnoticed as a shock metamorphic feature.

Description: Platey grains of cubic Bi2O3, alpha-Bi2O3, and Bi2O(2.75), nanograins were associated with chondritic porous interplanetary dust particles W7029C1, W7029E5, and 2011C2 that were collected in the stratosphere at 17-19 km altitude. Similar Bi oxide nanograins were present in the upper stratosphere during May 1985. These grains are linked to the plumes of several major volcanic eruptions during the early 1980s that injected material into the stratosphere. The mass of sulfur from these eruptions is a proxy for the mass of stratospheric Bi from which we derive the particle number densities (p/cu m) for "average Bi2O3 nanograins" due to this volcanic activity and those necessary to contaminate the extraterrestrial chondritic porous interplanetary dust particles via collisional sticking. The match between both values supports the idea that Bi2O3 nanograins of volcanic origin could contaminate interplanetary dust particles in the Earth's stratosphere.

Description: The top-of-atmosphere (TOA) path radiance generated by an aerosol mixture can be synthesized by linearly adding the contributions of the individual aerosol components, weighted by their fractional optical depths. The method, known as linear mixing, is exact in the single-scattering limit. When multiple scattering is significant, the method reproduces the atmospheric path radiance of the mixture with less than 3% errors for weakly absorbing aerosols up to optical thickness of 0.5. However, when strongly absorbing aerosols are included in the mixture, the errors are much larger. This is due to neglecting the effect of multiple interactions between the aerosol components, especially when the values of the single-scattering albedos of these components are so different that the parameter epsilon = (Sigma)f(sub i) absolute value of bar omega(sub i) - bar omega(sub mix)/bar omega(sub i), is larger than approx. 0.1, where bar omega(sub i) and f(sub i) are the single-scattering albedo and the fractional abundance of the i th component, and bar omega(sub i) is the effective single-scattering albedo of the mixture. We describe an empirical, modified linear-mixing method which effectively accounts for the multiple interactions between aerosol components. The modified and standard methods are identical when epsilon = 0.0 and give similar results when epsilon less than or equal to 0.05. For optical depths larger than approx. 0.5, or when epsilon greater than 0.05, only the modified method can reproduce the radiances within 5% error for common aerosol types up to optical thickness of 2.0. Because this method facilitates efficient and accurate atmospheric path radiance calculations for mixtures of a wide variety of aerosol types, it will be used as part of the aerosol retrieval methodology for the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), scheduled for launch into polar orbit in 1998.

Description: Quasi-electrostatic (QE) fields which exist above thunderclouds after lightning discharges can lead to the formation of columnar channels of breakdown ionization and carrot-like vertical luminous structures with typical transverse dimension approximately 5-10 km spaning an altitude range from approximately 80 km to well below approximately 50 km. The carrot-like forms closely resemble those observed in sprites. Results indicate that the appearance of optical emissions can be significantly delayed in time (approx. 1-20 ms) with respect to the causative lightning discharge.

Description: Trajectory calculations are used to examine ozone transport in the polar winter stratosphere during periods of the Upper Atmosphere Research Satellite (UARS) observations. The value of these calculations for determining mass transport was demonstrated previously using UARS observations of long-lived tracers, In the middle stratosphere, the overall ozone behavior observed by the Microwave Limb Sounder in the polar vortex is reproduced by this purely dynamical model. Calculations show the evolution of ozone in the lower stratosphere during early winter to be dominated by dynamics in December 1992 in the Arctic. Calculations for June 1992 in the Antarctic show evidence of chemical ozone destruction and indicate that approx. 50% of the chemical destruction may be masked by dynamical effects, mainly diabatic descent, which bring higher ozone into the lower-stratospheric vortex. Estimating differences between calculated and observed fields suggests that dynamical changes masked approx. 20% - 35% of chemical ozone loss during late February and early March 1993 in the Arctic. In the Antarctic late winter, in late August and early September 1992, below approx. 520 K, the evolution of vortex-averaged ozone is entirely dominated by chemical effects; above this level, however, chemical ozone depletion can be partially or completely masked by dynamical effects. Our calculations for 1992 showed that chemical loss was nearly completely compensated by increases due to diabatic descent at 655 K.

Description: We use geodetic data taken over four years with the Global Positioning System (GPS) to estimate: (1) motion between six major plates and (2) motion relative to these plates of ten sites in plate boundary zones. The degree of consistency between geodetic velocities and rigid plates requires the (one-dimensional) standard errors in horizontal velocities to be approx. 2 mm/yr. Each of the 15 angular velocities describing motion between plate pairs that we estimate with GPS differs insignificantly from the corresponding angular velocity in global plate motion model NUVEL-1A, which averages motion over the past 3 m.y. The motion of the Pacific plate relative to both the Eurasian and North American plates is observed to be faster than predicted by NUVEL-1A, supporting the inference from Very Long B ase- line Interferometry (VLBI) that motion of the Pacific plate has speed up over the past few m.y. The Eurasia-North America pole of rotation is estimated to be north of NUVEL-1A, consistent with the independent hypothesis that the pole has recently migrated northward across northeast Asia to near the Lena River delta. Victoria, which lies above the main thrust at the Cascadia subduction zone, moves relative to the interior of the overriding plate at 30% of the velocity of the subducting plate, reinforcing the conclusion that the thrust there is locked beneath the continental shelf and slope.

Description: A composite model is proposed to describe the time-dependent response of the Earth's lower crust. The motivation for such it model is twofold: First, new observations of widespread postseismic deformation indicate that the deep continental crust responds viscoelastically, having both long-and short-term decay times. Second, by any number of observationally based rationales, the lower crust is compositionally and structurally heterogeneous over many length scales. For heterogeneities that have much smaller characteristic lengths than the minimum deformation wavelength of interest, the aggregate rheology can be described by composite media theory. For wavelengths of the order of the thickness of the lower crust (approx. = 25-40 km) and larger, composite theory may be applied to heterogeneities that are smaller than about several hundred meters, or equivalent to the vertical extent of a thick lower crustal mylonitic shear zone. The composite media theory developed here is constructed using both Eshelhy-Mori-Tanaka theory for aligned generalized spheroidal inclusions and a generalized self-consistent method. The inclusions and matrix are considered to be Maxwellian viscoelastic: a rheology that is consistent with past homogeneous models of postseismic stress relaxation. The composite theory presented here introduces a transient response to a suddenly imposed stress field which does not appear in homogeneous Maxwell models. Analytic expressions for the amplitude and duration of the transient and for the effective long-and short-term viscosities of the composite are given which describe the sensitivity to inclusion concentration (phi), to shape, and to ratio of inclusion-to-matrix viscosity (R).

Description: The ability of water vapor radiometers (WVRs) to calibrate changes in tropospheric delay was demonstrated during very long baseline radio interferometer (VLBI) observations at Goldstone, California. WVR measurements reduced the observed VLBI delay variations over a 13 hr period by a factor of approx. = 2.5. When applied to shorter time scales, a approx. = 50% reduction in 100-700 s delay variations was achieved during conditions of high tropospheric activity. Thermal WVR noise precluded calibration of short time scale delay fluctuations during quiet tropospheric conditions.

Description: Measurements of hydrogen, nitrogen and chlorine radicals from a balloon flight on 25 September 1993 from Ft. Sumner, NM provide an opportunity to quantify photochemical production and loss of stratospheric ozone. Ozone loss rates determined using measured radical concentrations agree fairly well with loss rates calculated using a photochemical model. Catalytic cycles involving OH and HO2 are shown to dominate photochemical loss of ozone for altitudes between 44 and 50 km. Reactions involving NO and NO2 are the dominant sink for ozone between 25 and 38 km. The total ozone loss rate determined from the measurements balances calculated production rates for altitudes between 30 and 40 km. However, loss of ozone exceeds production by -35% between 42 and 50 km. The imbalance between production and loss of ozone above 42 km is larger than the uncertainty of any one of the critical kinetic parameters or species concentrations. No single adjustment to any of these parameters can simultaneously resolve the imbalance and satisfy constraints imposed by measured OH, HO2, NO2 and ClO. Our results are consistent with an additional mechanism for ozone production above 40 km other than photolysis of ground state O2.

Description: An anelastic approximation is used with a time-variable coordinate transformation to formulate a two-dimensional numerical model that describes the evolution of gravity waves. The model is solved using a semi-Lagrangian method with monotone (nonoscillatory) interpolation of all advected fields. The time-variable transformation is used to generate disturbances at the lower boundary that approximate the effect of a traveling line of thunderstorms (a squall line) or of flow over a broad topographic obstacle. The vertical propagation and breaking of the gravity wave field (under conditions typical of summer solstice) is illustrated for each of these cases. It is shown that the wave field at high altitudes is dominated by a single horizontal wavelength; which is not always related simply to the horizontal dimension of the source. The morphology of wave breaking depends on the horizontal wavelength; for sufficiently short waves, breaking involves roughly one half of the wavelength. In common with other studies, it is found that the breaking waves undergo "self-acceleration," such that the zonal-mean intrinsic frequency remains approximately constant in spite of large changes in the background wind. It is also shown that many of the features obtained in the calculations can be understood in terms of linear wave theory. In particular, linear theory provides insights into the wavelength of the waves that break at high altitudes, the onset and evolution of breaking. the horizontal extent of the breaking region and its position relative to the forcing, and the minimum and maximum altitudes where breaking occurs. Wave breaking ceases at the altitude where the background dissipation rate (which in our model is a proxy for molecular diffusion) becomes greater than the rate of dissipation due to wave breaking, This altitude, in effect, the model turbopause, is shown to depend on a relatively small number of parameters that characterize the waves and the background state.

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.

Description: This report describes the work done with funding from NASA Grant NAGW 3525 during 1995. The grant was initiated in March, 1992 and has a planned duration of three years. This report covers the time period March 1995 to Present. In this report I present a short description of the projects carried out and documentation of the work done in terms of papers presented, etc.

Description: Two rockets were flown to peak altitudes of 220 km in September 1959 to test various methods planned for future measurements of ionization parameters in the ionosphere, exosphere, and interplanetary plasma. The experiments used techniques which sample the ambient environment in the immediate vicinity of the research vehicle. Direct methods were chosen since indirect propagation techniques do not provide the temperatures of charged particles, are insensitive to ion densities, and cannot measure local electron densities under all conditions. Very encouraging results have been obtained from a preliminary analysis of data provided by one of the two flights. A new rf probe technique was successfully used to determine the electron density profile. This was indicated by its agreement with the results of a companion cw propagation experiment, particularly when the probe data were corrected for the effects of the ion sheath which surrounds the vehicle. The characteristics of this sheath were determined directly in flight by an electric field meter which provided the sheath field, and by a Langmuir probe which measured the total potential across the sheath. The electron temperatures deduced from the Langmuir probe data are greater than the neutral gas temperatures previously measured for the same location and season, but these measurements possibly were taken under different atmospheric conditions. Ion densities were calculated from the ion trap data for several altitudes ranging from 130 to 210 km and were found to be within 20 percent of the measured electron densities.

Description: Measurements of clear-air turbulence by use of airplane-borne instrumentation have been obtained from NASA VGH recorders during research flights of Lockheed U-2 airplanes to altitudes of 75,000 feet over several areas of the Northern Hemisphere. An analysis of these data has indicated that for the higher altitudes (50,000 to 75,000 feet), turbulence is both less frequent and less severe than for the lower altitudes (20,000 to 50,000 feet). Turbulence appears to be present at the high altitudes (60,000 to 75,000 feet) less than 1 percent of the time. Moderately heavy turbulence appears to exist on occasion at altitudes of about 50,000 feet over Japan. As a consequence, the gust experience appears to be more severe for operations over Japan than for the other areas. Less than 50 percent of the turbulent areas exceeded 10 miles in length.

Description: Abundant evidence now shows that the buried Chicxulub structure in northern Yucatan, Mexico, is indeed the intensely sought-after source of the ejecta found world-wide at the Cretaceous-Tertiary (K/T) boundary. In addition to large-scale concentric patterns in gravity and magnetic data over the structure, recent analyses of drill-core samples reveal a lithological assemblage similar to that observed at other terrestrial craters. This assemblage comprises suevite breccias, ejecta deposit breccias (Bunte Breccia equivalents), fine-grained impact melt rocks, and melt-matrix breccias. All these impact-produced lithologies contain diagnostic evidence of shock metamorphism, including planar deformation features in quartz, feldspar, and zircons; diaplectic glasses of quartz and feldspar; and fused mineral melts and whole-rock melts. In addition, elevated concentrations of Ir, Re, and Os, in meteoritic relative proportions, have been detected in some melt-rock samples from the center of the structure. Isotopic analyses, magnetization of melt-rock samples, and local stratigraphic constraints identify this crater as the source of K/T boundary deposits.

Description: I discuss errors in theory and in interpreting observations that are produced by the failure to consider resolution in space, time, and energy. I discuss convection in stellar model atmospheres and in stars. Large errors in abundances are possible such as the factor of ten error in the Li abundance for extreme Population II stars. Finally I discuss the variation of microturbulent velocity with depth, effective temperature, gravity and abundance. These variations must be dealt with in computing models and grids and in any type of photometric calibration.

Description: On Venus, global topography shows the presence of highs and lows including regional highly deformed plateaus (tesserae), broad rifted volcanic rises, linear lows flanking uplands, and more equidimensional lowlands (e.g. Lavinia and Atalanta planitiae) Each of these terrain types on Venus has relatively distinctive characteristics, but origins are uncertain in terms of mode of formation, time of formation, and potential evolutionary links. There is a high level of uncertainty about the formation and evolution of lowlands on Venus. We have undertaken the mapping of a specific lowlands region of Venus to address several of these major questions. Using geologic mapping we have tried to establish: What is the sequence of events in the formation and evolution of large-scale equidimensional basins on Venus? When do the compressional features typical of basin interiors occur? What is the total volume of lava that occurs in the basins and is this similar to other non-basin areas? How much subsidence and downwarping has occurred after the last major plains units? WE have undertaken an analysis of the geology of the V55 Lavinia Planitia quadrangle in order to address many of these issues and we report on the results here.

Description: The dynamics and structure of the magnetosheath and the magnetopause, a new theory of MHD shock waves, evolution of intermediate shocks, and shock waves in an anisotropic plasma, are reported.

Description: Most of the work carried out to date on this project is summarized in the enclosed reprints of two papers that were just published. The earlier paper that is also enclosed, Structure of the Magnetotail, by D. L. Larson and R. L. Kaufmann, was primarily intended to describe our Consistent Orbit Tracing (COT) technique and to show that the resulting magnetotail models were in good agreement with published experimental observations. The following are the most important results from the two new papers.

Description: A final report for the NASA-supported project on laboratory investigations of stratospheric halogen chemistry is presented. In recent years, this project has focused on three areas of research: (1) kinetic, mechanistic, and thermochemical studies of reactions which produce weakly bound chemical species of atmospheric interest; (2) development of flash photolysis schemes for studying radical-radical reactions of stratospheric interest; and (3) photochemistry studies of interest for understanding stratospheric chemistry. The first section of this paper contains a discussion of work which has not yet been published. All subsequent chapters contain reprints of published papers that acknowledge support from this grant.

Description: The Slate Islands in northern Lake Superior represent the eroded remains of a complex impact crater, originally approximately 32 km in diameter. New field studies there reveal allogenic crater fill deposits along the eastern and northern portions of the islands indicating that this 500-800 Ma impact structure is not as heavily eroded as previously thought. Near the crater center, on the western side or Patterson Island, massive blocks of target rocks, enclosed within a matrix of fine-grained polymict breccia, record the extensive deformation associated with the central uplift. Shatter cones are a common structural feature on the islands and range from less than 3 cm to over 10 m in length. Although shatter cones are powerful tools for recognizing and analyzing eroded impact craters, their origin remains poorly constrained.

Description: The near earth reconnection model of substorms represents an attempt to place a broad range of observations into a consistent framework. The roles and requirements of reconnection are discussed. High speed plasma sheet flows, thin current sheet instability, substorm triggering, plasmoids and flux ropes in the distant tail, and magnetohydrodynamic simulations are discussed. Substorms are global, coherent sequences of processes involving solar wind/magnetosphere/ionosphere interaction. Magnetic reconnection is required to explain different dayside and polar cap phenomena, which required nightside reconnection. The modification and expansion of the standard near earth neutral line (NENL) model can integrate breakup arcs, current disruption, current wedge features, and localized plasma flows into the magnetic reconnection framework.

Description: Recent results from investigations carried out using the active magnetospheric particle tracer explorer/ion release module (AMPTE/IRM) satellite and International Sun-Earth Explorer (ISEE) show that bursty, near-neutral sheet fast flows are an important part of the transport in the near-earth tail. The results related to such flows are reviewed, and information is presented on their relation to the average plasma sheet characteristics and substorms. The average ion density, temperature and flow variability in the quiet inner plasma sheet exhibit a spatial dependence which suggests a dependence on Bursty Bulk Flows (BBFs). The energy and flux transport in the tail is discussed and it is stated that BBFs are responsible for such transport. The tailward progression of activity is considered.

Description: Many substorm expansions are triggered by interplanetary magnetic field changes that reduce magnetospheric convection. This suggests that expansion onsets are a result of a reduction in the large-scale electric field imparted to the magnetosphere from the solar wind. Such a reduction disrupts the inward motion and energization of plasma sheet particles that occur during the growth phase. It is proposed that the resulting magnetic drift of particles and a large dawn to dusk gradient in the ion energies leads to a longitudinally localized reduction in the plasma pressure, and thus, to the current wedge formation. This theory accounts for the rapid development of the expansion phase relative to growth phase, the magnitude of the wedge currents, the speeds of tailward and westward expansion of the current reduction region in the equatorial plane, and the speeds of the poleward and westward motion of active aurora in the ionosphere.

Description: Over the past four years of funding, SRI, in collaboration with the University of Texas at Dallas, has been involved in assessing the influence of thermospheric neutral winds on the electric field and current systems at high latitudes. The initial direction of the project was to perform a set of numerical experiments concerning the contribution of the magnetospheric and neutral wind dynamo processes, under specific boundary conditions, to the polarization electric field and/or the field-aligned current distribution at high latitudes. To facilitate these numerical experiments we developed a numerical scheme that relied on using output from the NCAR Thermosphere-Ionosphere General Circulation Model (NCAR-TIGCM), expanding them in the form of spherical harmonics and solving the dynamo equations spectrally. Once initial calculations were completed, it was recognized that the neutral wind contribution could be significant but its actual contribution to the electric field or currents depended strongly on the generator properties of the magnetosphere. Solutions to this problem are not unique because of the unknown characteristics of the magnetospheric generator, therefore the focus was on two limiting cases. One limiting case was to consider the magnetosphere as a voltage generator delivering a fixed voltage to the high-latitude ionosphere and allowing for the neutral wind dynamo to contribute only to the current system. The second limiting case was to consider the magnetosphere as a current generator and allowing for the neutral wind dynamo to contribute only to the generation of polarization electric fields. This work was completed and presented at the l994 Fall AGU meeting. The direction of the project then shifted to applying the Poynting flux concept to the high-latitude ionosphere. This concept was more attractive as it evaluated the influence of neutral winds on the high-latitude electrodynamics without actually having to determine the generator characteristics of the magnetosphere. The influence of the neutral wind was then determined not by estimating how much electric potential or current density it provides, but by determining the contribution of the neutral wind to the net electromagnetic energy transferred between the ionosphere and magnetosphere. The estimate of the net electromagnetic energy transfer and the role of the neutral winds proves to be a more fundamental quantity in studies of magnetosphere- ionosphere coupling also showed that by using electric and magnetic field measurements from the HILAT satellite, the Poynting flux could be a measurable quantity from polar-orbiting, low- altitude spacecraft. Through collaboration with Dr. Heelis and others at UTD and their expertise of the electric field measurements on the DE-B satellite, an extensive analysis was planned to determine the Poynting flux from the DE-B measurements in combination with a modeling effort to help interpret the observations taking into account the coupled magnetosphere-ionosphere.

Description: The basaltic Kaupulehu 1800-1801 lava flow of Hualalai Volcano, Hawaii contains abundant ultramafic xenoliths. Many of these xenoliths occur as bedded layers of semi-rounded nodules, each thinly coated with a veneer (typically 1 mm thick) of lava. The nodule beds are analogous to cobble deposits of fluvial sedimentary systems. Although several mechanisms have been proposed for the formation of the nodule beds, it was found that, at more than one locality, the nodule beds are overbank levee deposits. The geological occurrence of the nodules, certain diagnostic aspects of the flow morphology and consideration of the inferred emplacement process indicate that the Kaupulehu flow had an exceptionally low viscosity on eruption and that the flow of the lava stream was extremely rapid, with flow velocities of at least 10 m/s (more than 40 km/h. This flow is the youngest on Hualalai Volcano and future eruptions of a similar type would pose considerable hazard to life as well as property.

Description: The proposal had outlined a year for program conversion, a year for testing and debugging, and two years for numerical experiments. We kept to that schedule. In first (partial) year, author designed a finite element for isostatic thin-shell deformation on a sphere, derived all of its algebraic and stiffness properties, and embedded it in a new finite element code which derives its basic solution strategy (and some critical subroutines) from earlier flat-Earth codes. Also designed and programmed a new fault element to represent faults along plate boundaries. Wrote a preliminary version of a spherical graphics program for the display of output. Tested this new code for accuracy on individual model plates. Made estimates of the computer-time/cost efficiency of the code for whole-earth grids, which were reasonable. Finally, converted an interactive graphical grid-designer program from Cartesian to spherical geometry to permit the beginning of serious modeling. For reasons of cost efficiency, models are isostatic, and do not consider the local effects of unsupported loads or bending stresses. The requirements are: (1) ability to represent rigid rotation on a sphere; (2) ability to represent a spatially uniform strain-rate tensor in the limit of small elements; and (3) continuity of velocity across all element boundaries. Author designed a 3-node triangle shell element which has two different sets of basis functions to represent (vector) velocity and all other (scalar) variables. Such elements can be shown to converge to the formulas for plane triangles in the limit of small size, but can also applied to cover any area smaller than a hemisphere. The difficult volume integrals involved in computing the stiffness of such elements are performed numerically using 7 Gauss integration points on the surface of the sphere, beneath each of which a vertical integral is performed using about 100 points.

Description: The research performed during this reporting period includes development and maintenance of scientific software for the GOME retrieval algorithms, consultation on operational software development for GOME, further sensitivity and instrument studies to help finalize the definition of the SCIAMACHY instrument, and consultation on optical and detector issues for both GOME and SCIAMACHY. The Global Ozone Monitoring Experiment was successfully launched on the ERS-2 satellite on April 20, 1995, during this reporting period, and is working in the expected fashion. The European Space Agency has made their selections from responses to the Announcement of Opportunity for GOME validation and science studies, part of the overall ERS AO. The Smithsonian Astrophysical Observatory (SAO) proposal has been selected. These proposals are primarily for access to the data; ESA does not provide research funding for the selected investigations. The SAO activities that are carried out as a result of selection by ESA are funded by the present grant, to the limit that can be accomplished at the present level of funding. SCIAMACHY is currently in Phase C/D. Instrument design is almost finalized and selection of infrared detectors from the initial production run has been made.

Description: The SUNDIAL-ATLAS effort was a global-scale investigation which responded to the science priorities of the ITM Panel, the Integrated SPD Strategy Implementation Plan as a whole, and the need for potential cost-saving design criteria for the TIMED mission. The investigation focused on coupling processes in the ionospheric-thermospheric system, taking advantage of the timelines of the ATLAS-1 mission (March 1992), and the global-scale ground-based measurement and modeling activities of the SUNDIAL program. The collaborative SUNDIAL-ATLAS activity was the first opportunity for global measurements of the chemistry, kinetics, and electrodynamics which couple the E-, Fl-, and F2-regions into a single interactive system. As such, the program represented an important first step in studying global issues; and accordingly, was an important proof of concept experiment relevant to the strategic mission plans for the ITM community and the upcoming intermediate class satellite program called TIMED. To meet its projected goals, TIMED must perform a number of critical measurements and execute a number of correlations that were to be tried and tested for the first time in the SUNDIAL-ATLAS investigation. This was designed to include global correlations of thermospheric and ionospheric composition during quiet and disturbed conditions and the co-registration of global-scale ground-based measurements with along-track satellite diagnostics. The SUNDIAL component of the current investigation addressed this need by acquiring, reducing, and analyzing a multi-sensor database that complemented and extended that which was generated in the ATLAS mission (Atmospheric Laboratory for Applications and Science). The SUNDIAL data defined the state and condition of the global-scale ionosphere in the altitude range from 100 km to the F2-peak. These data specified the peak heights and densities of the E-, Fl-, and F2-regions, along with the global distributions of intermediate, descending, and sequential layers which play a critical role in the dynamo region of the lower ionospheric-thermospheric domain. The data were collected by the SUNDIAL network of more than 50 ground-based stations utilizing ionosondes, radars, photometers, Fabry-Perot interferometers, and total electron content measurements. The data were acquired during a three-week period centered on the eight-day ATLAS-1 mission, which provided image and photometric sensing of the altitude distributions of the major and minor ions and neutrals in the ITM system. This report focuses on the scientific contributions of the SUNDIAL component of the overall investigation. Specific findings are described in seven papers (attached) published in the Journal of Geophysical Research.

Description: The work proposed on this project included both field and laboratory studies. The laboratory studies were to consist of measurements of H2SO4 uptake and evaporation from aerosols of varying chemical composition, while the field component would include measurements of H2SO4 and other compounds which would be conducted as part of a large field campaign. By chance, the opportunity to conduct such an H2SO4/aerosol/ultrafine particle study in conjunction with an OH intercomparison/photochemistry study became available very early in this project (September 1993). This study was conducted at Caribou, Colorado in conjunction with several other groups from NCAR, NOAA and a number of universities. Our group measured OH, H2SO4, SO2, and H20, while Dr. McMurfy's group measured ultrafine particles, and total particle number and size distribution. In addition measurements of HO2/RO2, O3, NO, NO2, NO(y) CO, hydrocarbons, CH2O, and other chemical compounds and meteorological parameters were performed by the other participants and a new laser oblation/mass spectrometry technique was also employed by the NOAA Aeronomy Laboratory to study aerosol composition. The study of aerosol production and growth in conjunction with photochemical measurements is highly advantageous because particle growth precursors such as H2SO4.or MSA are formed by OH initiated sulfur oxidation. The large number of hydrocarbon measurements included in this study were also important in understanding particle growth.

Description: The thermal emission of two palagonitic soils, common visible and near infrared spectral analogs for bright soils on Mars, was measured over the wavelength range of 5 to 25 microns (2000 to 400 cm(exp -1) for several particle size separates. All spectra exhibit emissivity features due to vibrations associated with H2O and SiO. The maximum variability of emissivity is approx. = 20% in the short wavelength region (5 to 6.5 microns, 2000 to 1500 cm(exp -1)), and is more subdued, 〈 4%, at longer wavelengths. The strengths of features present in infrared spectra of Mars cannot be solely provided by emissivity variations of palagonite; some other material or mechanism must provide additional absorption(s).

Description: Fifteen polar ice samples-fourteen from Greenland and one from Antarctica-have been analyzed for the extraterrestrial amino acid alpha-aminoisobutyric acid (AIB) in an effort to estimate the flux of interplanetary organic material to the Earth's surface. Only one sample (Greenland GISP II, 4270-4440 years old) contains detectable amounts of AIB, apparently the signature of a transient delivery event. The maximum oceanic concentration of AIB from such an event would be less than 10(exp-9) M.

Description: The geoid and the boundary topography caused by mass loads inside the earth were estimated. It is shown that the estimates are affected by compressibility, by a radially varying density distribution, and by the presence of phase boundaries with density discontinuities. The geoid predicted in the chemical boundary case is 30 to 40 percent smaller than that predicted in the phase case. The effects of compressibility and radially varying density are likely to be small. The inner core-outer core topography for loading inside the mantle and for loading inside the inner core were computed.

Description: This document provides a summary of the launch and post-launch activities of Virginia Space Grant Consortium Upper Atmospheric Payload Balloon System, V(ps). It is a comprehensive overview covering launch activities, post-launch activities, experimental results, and future flight recommendations.

Description: The initial substorm interval to be studied has been selected to be the magnetic storm of June 4, 1991, and following. The CRRES electric and magnetic field data has been processed. We have added DMSP ion drift and energetic particle data and further refined the CRRES data for intercomparison. The DMSP data increase the frequency of monitoring of the temporal response of the penetration electric fields to every 100 min. Energy is seen to flow between the ionosphere and magnetosphere at low L values during the main phase of the magnetic storm in the form of field-aligned Poynting flux. This indicates electrodynamic coupling of the regions with Alfven waves. The first comparisons of the data with outputs of the Rice Magnetospheric Specification Model (MSM) were made. Both positive and negative correlation were seen as might be expected. Differences were especially evident in the time constants of the processes. Comparisons with the more physically self-consistent Rice Convection Model (RCM) with both electric fields and particle data are in progress to suggest physical constraints for our understanding of the phenomena.

Description: A full-wave propagation model was developed that describes the propagation of gravity waves from the Earth's surface to the upper boundary, which can be placed anywhere between 150 and 500 km altitude. The model includes a realistic background atmosphere, and includes the effects of mean horizontal winds and their vertical shears, mean vertical temperature gradients, the eddy and molecular diffusion of heat and momentum, and the effects of ion-drag. This model solves five coupled second-order differential equations (continuity, momentum, and energy) in the vertical coordinate to derive the perturbation variables u', v', w' (horizontal and vertical velocity components), T' (temperature) and p' (pressure). The upper boundary can be automatically selected based on tests using the radiation condition at the upper boundary, wherein the height is increased until the wave is experiencing severe dissipation at the upper boundary, ensuring that substantial absorption occurs for any waves reflected from the upper boundary. The determination of wave amplitude is a key requirement of wave energetics. Therefore, the fullwave model has been applied to airglow observations in order to determine wave amplitudes as a function of altitude. This was accomplished by using the full-wave model output to drive a chemistry perturbation module that describes minor species perturbations and the resulting airglow perturbations. The full-wave output was multiplied by an altitude-independent factor such that the modeled and observed relative airglow intensity perturbations were equal. The effects of mean winds were included in these studies, and found to be the most important model input affecting the calculations (being more important than the choice of eddy diffusion profiles and chemical kinetic coefficients). In one study (Hickey et al., 1997a) these winds could not be well estimated from the measurements, whereas in the second study (Hickey et al.,1997b) the mean were well defined with a sodium wind-temperature lidar.

Description: The frozen-flux hypothesis for the Earth's liquid core assumes that convective terms dominate diffusive terms in the induction equation governing the behaviour of the magnetic field at the surface of the core. While highly plausible on the basis of estimates of physical parameters, the hypothesis has been questioned. To study this hypothesis, we improve the method which tests the consistency of magnetic observations with the hypothesis by constructing simple, flux-conserving core-field models fitting the data at pairs of epochs. We introduce a new approach that fixes the patch configurations at each of the two epochs before inversion, so that each configuration is consistent with its respective data set but possesses the same patch topology. We expand upon the inversion algorithm, using quadratic programming to maintain the proper flux sign within patches; the modelling calculations are also extended to include data types that depend non-linearly on the model. Every test of a hypothesis depends on the characterization of the observational uncertainties; we undertake a thorough review of this question. For main-field models, the primary source of uncertainty comes from the crustal field. We base our analysis on statistical models of the crustal magnetization, adjusted to bring it into better conformity with our data set. The noise model permits us to take into account the correlations between the measurements and requires that a different weighting be given to horizontal and vertical components. It also indicates that the observations should be fit more closely than has been the practice heretofore. We apply the revised method to Magsat data from 1980 and survey and observatory data from 1915.5, two data sets believed to be particularly difficult to reconcile with the frozen-flux hypothesis. We compute a pair of simple, flux-conserving models that fit the averaged data from each epoch. We therefore conclude that present knowledge of the geomagnetic fields of 1980 and 1915.5 is consistent with the frozen-flux hypothesis.

Description: One of the oldest mysteries in geomagnetism is the linkage between solar and geomagnetic activity. The 11-year cycles of both the numbers of sunspots and Earth geomagnetic storms were first noted by Sabine. A few years later, speculation on a causal relationship between flares and storms arose when Carrington reported that a large magnetic storm followed the great September 1859 solar flare. However, it was not until this century that a well-accepted statistical survey on large solar flares and geomagnetic storms was performed, and a significant correlation between flares and geomagnetic storms was noted. Although the two phenomena, one on the Sun and the other on the Earth, were statistically correlated, the exact physical linkage was still an unknown at this time. Various hypotheses were proposed, but it was not until interplanetary spacecraft measurements were available that a high-speed plasma stream rich in helium was associated with an intense solar flare. The velocity of the solar wind increased just prior to and during the helium passage, identifying the solar ejecta for the first time. Space plasma measurements and Skylab's coronagraph images of coronal mass elections (CMES) from the Sun firmly established the plasma link between the Sun and the Earth. One phenomenon associated with magnetic storms is brilliant "blood" red auroras, as shown.

Description: On 9 April 1994 the Clementine spacecraft high-gain antenna was aimed toward the Moon's surface and the resulting 13-cm wavelength radio echoes were received on Earth. Using these data, we have found that the lunar surface generally follows a Lambertian bistatic scattering function sigma(sub 0) = K(sub D)cos(theta(sub i) with K(sub D) approx. 0.003 for the opposite (expected) sense of circular polarization and K(sub D) approx. 0.001 for the same (unexpected) sense. But there are important deviations-of up to 50% in some parts of the echo spectrum-from this simple form. Based on an earlier analysis of these same data, Nozette et al. claimed detection of an enhancement in echoes with right circular polarization from regions near the South Pole in a near-backscatter geometry. Such behavior would be consistent with presence of perhaps large quantities of water ice near the Pole. We have been unable to reproduce that result. Although we find weak suggestions of enhanced echoes at the time of South Pole backscatter, similar features are present at earlier and later times, adjacent frequencies, and in left circular polarization. If enhanced backscatter is present, it is not unique to the South Pole; if not unique to the Pole, then ice is less likely as an explanation for the enhancement.

Description: The research performed during this reporting period has been a collaboration between institutions including the Smithsonian Astrophysical Observatory, the National Institute of Standards and Technology, the University of Oregon, and the NASA Langley Research Center. The program has included fully line-resolved measurements of submillimeter and far infrared spectroscopic line parameters (pressure broadening coefficients and their temperature dependences, and fine positions) for the analysis of field measurements of stratospheric constituents, far infrared database improvements, and studies for improved satellite measurements of the Earth's atmosphere. This research program is designed to enable the full utilization of spectra obtained in far infrared/submillimeter field measurements, such as FIRS-2, FILOS, IBEX, SLS, EosMLS, and proposed NASA and European Space Agency measurements of ClO and OH (e.g., PIRAMHYD) for the retrieval of accurate stratospheric altitude profiles of key trace gases involved in ozone layer photochemistry.

Description: Three-dimensional, global MHD simulations of solar wind flow onto a prescribed magnetopause obstacle are used to show that a bow shock's nose location a(sub s), and the relative subsolar magnetosheath thickness Delta(sub ms)/a(sub mp) are strong functions of the IMF cone angle theta (between v(sub sw) and B(sub sw)) and the Alfven Mach number M(sub A). For a given M(sub A) the shock is more distant for higher theta (restricted to the interval 0-90deg by symmetries), while a(sub s)/a(sub mp) and Delta(sub ms/a(sub mp) increase with decreasing M(sub A) for theta greater than or approximately 20deg but decrease with decreasing M(sub A) for theta approximately Odeg. Large differences in Delta(sub ms/a(sub mp) are predicted between theta = Odeg and 90deg at low M(sub A), with smaller differences remaining even at M(sub A) approximately 10. The theta = Odeg results confirm and extend the previous work of Spreiter and Rizzi [1974]. The simulations show that successful models for the subsolar shock location cannot subsume the dependences on M(sub A) and theta into a sole dependence on M(ms). Instead, they confirm a recent prediction [Cairns and Grabbe, 1994] that a(sub s)/a(sub mp) and Delta(sub ms)/a(sub mp) should depend strongly on theta and M(sub A) for M(sub A) less than or approximately 10 (as well as on other MHD variables). Detailed comparisons between theory and data remain to be done. However, preliminary comparisons show good agreement, with distant shock locations found for low M(sub A) and large theta greater than or approximately 45deg and closer locations found for theta less than or approximately 20deg even at M/A approximately 8.

Description: Observational studies were performed of westward propagating synoptic scale waves in the tropical troposphere, the structure of monsoon circulations in the upper troposphere and lower stratosphere, and zonally propagating features in deep tropical convection. The effect of the quasi-biennial oscillation (QBO) were investigated, and a numerical study of the QBO was performed using a two-dimensional model, highlighting the role of gravity waves in the momentum balance of the QBO. Vertical coupling of the troposphere and stratosphere was examined in polar regions on intraseasonal and interannual timescales. A deep circumpolar mode was discovered, now known as the Arctic Oscillation.

Description: Conservation of angular momentum dictates that as the wind-driven axial atmospheric angular momentum changes, so will the length-of-day (LOD). In particular, as the strength of the seasonal zonal winds change, so should the strength of the seasonal LOD signals. Here, observed changes in the strengths of the annual and semiannual LOD signals during 1963-1991 are analyzed and shown to be both significantly correlated (at the 99% significance level) with the Southern Oscillation Index (SOI), and to exhibit trends of comparable magnitude but opposite signs. This reported correlation between the SOI and changes in the amplitude of the seasonal LOD signals demonstrates a linkage between seasonal LOD (and hence seasonal zonal wind) variability and the El Nino / Southern Oscillation (ENSO) phenomenon. Furthermore, this study suggests that observed variations in the amplitudes of the seasonal LOD signals can be used to study changes in the strengths of the seasonal atmospheric zonal winds on interannual to decadal and longer time scales.

Description: The Upper Atmosphere Research Satellite (UARS) microwave limb sounder (MLS) makes measurements of thermal emission at 183.3 GHz which are used to infer the concentration of water vapor over a pressure range of 46-0.2hPa (approximately 20-60 km). We provide a validation of MLS H2O by analyzing the integrity of the measurements, by providing an error characterization, and by comparison with data from other instruments. It is estimated that version 3 MLS H2O retrievals are accurate to within 20-25% in the lower stratosphere and to within 8-13% in the upper stratosphere and lower mesosphere. The precision of a single profile is estimated to be approximately 0.15 parts per million by volume (ppmv) in the midstratosphere and 0.2 ppmv in the lower and upper stratosphere. In the lower mesosphere the estimate of a single profile precision is 0.25-0.45 ppmv. During polar winter conditions, H2O retrievals at 46 hPa can have a substantial contribution from climatology. The vertical resolution of MLS H2O retrievals is approximately 5 km.

Description: On February 9, 1995, the Comprehensive Plasma Instrumentation (CPI) on the Geotail spacecraft observed a complex, structured ion distribution function near the magnetotail midplane at x approximately -30 R(sub E). On this same day the Wind spacecraft observed a quiet solar wind and an interplanetary magnetic field (IMF) that was northward for more than five hours, and an IMF B(sub y) component with a magnitude comparable to that of the RAF B(sub z) component. In this study, we determined the sources of the ions in this distribution function by following approximately 90,000 ion trajectories backward in time, using the time-dependent electric and magnetic fields obtained from a global MHD simulation. The Wind observations were used as input for the MHD model. The ion distribution function observed by Geotail at 1347 UT was found to consist primarily of particles from the dawn side low latitude boundary layer (LLBL) and from the dusk side LLBL; fewer than 2% of the particles originated in the ionosphere.

Description: Understanding the large-scale dynamics of the magnetospheric boundary is an important step towards achieving the ISTP mission's broad objective of assessing the global transport of plasma and energy through the geospace environment. Our approach is based on three-dimensional global magnetohydrodynamic (MHD) simulations of the solar wind-magnetosphere- ionosphere system, and consists of using interplanetary magnetic field (IMF) and plasma parameters measured by solar wind monitors upstream of the bow shock as input to the simulations for predicting the large-scale dynamics of the magnetospheric boundary. The validity of these predictions is tested by comparing local data streams with time series measured by downstream spacecraft crossing the magnetospheric boundary. In this paper, we review results from several case studies which confirm that our MHD model reproduces very well the large-scale motion of the magnetospheric boundary. The first case illustrates the complexity of the magnetic field topology that can occur at the dayside magnetospheric boundary for periods of northward IMF with strong Bx and By components. The second comparison reviewed combines dynamic and topological aspects in an investigation of the evolution of the distant tail at 200 R(sub E) from the Earth.

Description: We present Geotail plasma and field observations from the middle magnetotail near X(sub GSE) = -46 R(sub E) for the time period 1400 to 1800 UT on December 14, 1994. During that period, the Wind satellite monitored the solar wind plasma and interplanetary magnetic field (IMF) upstream of the bow shock. The IMF was northward and the plasma parameters near average. Geotail observed slow tailward flows and a northward field. The plasma and field parameters indicate that Geotail is either in the plasma sheet or in a boundary layer. We used the Wind solar wind plasma and IMF data as input for a global simulation of that time interval. Comparison of the simulation results with the observational data show very good overall agreement of the magnitudes of the plasma and field parameters. In particular, the simulation reproduces the slow tailward flows and northward field found at Geotail. Small scale temporal, variations are less well reproduced. The simulation shows the formation of a broad boundary layer (which we call tail flank boundary layer, TFBL) that consists of closed flux which is formed by magnetic magnetic reconnection of IMF and lobe field lines. The simulation results indicate that Geotail is located very close to the TFBL and may have entered the TFBL proper. We show that the TFBL plays an important role in energy transport from the solar wind into the magnetosphere during northward IMF conditions.

Description: Multiangle, multispectral remote sensing observations, such as those anticipated from the Earth Observing System (EOS) multiangle imaging spectroradiometer (MISR), can distinguish spherical from nonspherical particles over calm ocean for mineral-dust-like particles with the range of sizes and column amounts expected under natural conditions. The ability to make such distinctions is critical if remote sensing of atmospheric aerosol properties is to provide significant new contributions to our understanding of the global-scale, clear-sky solar radiation balance. According to theoretical simulations the measurements can retrieve column optical depth for nonspherical particles to an accuracy of at least 0.05 or 10%, whichever is larger. In addition, three to four distinct size groups between 0.1 and 2.0 microns effective radius can be identified at most latitudes.

Description: Remote airborne measurements of the vertical and horizontal structure of the atmospheric pressure field in the lower troposphere are made with an oxygen differential absorption lidar (DIAL). A detailed analysis of this measurement technique is provided which includes corrections for imprecise knowledge of the detector background level, the oxygen absorption fine parameters, and variations in the laser output energy. In addition, we analyze other possible sources of systematic errors including spectral effects related to aerosol and molecular scattering interference by rotational Raman scattering and interference by isotopic oxygen fines.

Description: A new method for estimating downwelling shortwave and longwave radiation fluxes in the Arctic from TOVS brightness temperatures has been developed. The method employs a neural network to bypass computationally intensive inverse and for-ward radiative transfer calculations, Results from two drifting ice camps (CEAREX, LeadEx) and from one coastal station show that downwelling fluxes can be estimated with r.m.s. errors of 20 W/sq m for longwave radiation and 35 W/sq m for shortwave radiation. Mean errors are less than 4 W/sq m and are well within the bounds required for many climate process studies.

Description: The rate of slow deformation processes along faults activated during the Landers 1992 earthquake was investigated. The analysis was performed by combining pairs of ERS-1/2 synthetic aperture radar (SAR) of Southern California (U.S.). The interferograms revealed several centimeters of post-seismic rebound in step-overs of the 1992 break. The southern branches of the 1992 break experienced surface creep, producing sharp phase cuts in the interferometric maps. The same approach was used in the Los Angeles (CA) basin. The tectonic signal in the interferograms of the Los Angeles basin was intermingled with signals due to ground subsidence caused by oil and water withdrawal.

Description: We have continued our modeling of the solar wind/magnetosphere interaction for steady solar Wind conditions by carrying out a series of 3 D global magneto-hydrodynamics (MHD) simulations using a set of predetermined solar wind parameters as input to the code. With the results from the simulations that we ran last year for two other dipole tilts (0 and 35 deg), we have 36 data sets to study the dynamics of mesoscale structures at the magnetospheric boundary. We have started the analysis of these runs to investigate geometrical properties and the topology of the magnetic and electric fields for the different solar wind regimes and orientations considered. Preliminary results indicate that merging sites are consistent with patterns proposed for antiparallel merging at the dayside magnetopause. Another goal of this investigation is to establish the displacement of the cusp region as a function of the solar wind dynamic pressure, IMF direction and magnetic field dipole tilt. One of the difficulties of that study is to locate precisely the cusp in the simulation results. We have used several case studies to establish a series of criteria in to permit a routine identification of that location. We need now to process the rest of the simulation results using these criteria to obtain statistical results.

Description: Over the past three years, the primary foci of our research have been on studies of stratospheric processes and the potential effects of natural perturbations and human related activities on global stratospheric ozone, on interactions between atmospheric chemistry, atmospheric composition and concerns about climate change, and, to a lesser degree, on studies of tropospheric processes. Published research articles are listed.

Description: Our research activities under this NASA grant have focused on two broad topics associated with the Earth Radiation Budget Experiment (ERBE): (1) the role of clouds and the surface in modifying the radiative balance; and (2) the spatial and temporal variability of the earth's radiation budget. Each of these broad topics is discussed separately in the text that follows. The major points of the thesis are summarized in section 3 of this report. Other dissertation focuses on deriving the radiation budget over the TOGA COARE region.

Description: This final report describes the activities under NASA contract to Lockheed Missiles and Space Company. It covers the period from 10-1-94 to 12-31-97. The objective of this investigation is to identify and characterize the spatial relationships of auroral particle acceleration features relative to the characteristic transition features in the surrounding polar ionospheric plasmas. Due to the reduced funding level approved for this contract, the original scope of the proposed work was readjusted with the focus placed on examining spatial relationships with respect to particle structures.

Description: This section outlines those tasks undertaken in the final year that contribute integrally to the overarching project goals. Fast, during the final year, it is important to note that the project benefited greatly with the addition of a Boston University graduate student, Ms. Karen Hirsch. Jointly, we made substantial progress on the development of and improvements to magnetotail magnetic field and plasma models. The ultimate aim of this specific task was to assess critically the utility of such models for mapping low-altitude phenomena into the magnetotail (and vice-versa). The bulk of this effort centered around the finite-width- magnetotail convection model developed by and described by Spence and Kivelson (J. Geophys. Res., 98, 15,487, 1993). This analytic, theoretical model specifies the bulk plasma characteristics of the magnetotail plasma sheet (number density, temperature, pressure) across the full width of the tail from the inner edge of the plasma sheet to lunar distances. Model outputs are specified by boundary conditions of the source particle populations as well as the magnetic and electric field configuration. During the reporting period, we modified this code such that it can be interfaced with the auroral particle precipitation model developed by Dr. Terry Onsager. Together, our models provide a simple analytic specification of the equatorial distribution of fields and plasma along with their low-altitude consequences. Specifically, we have built a simple, yet powerful tool which allows us to indirectly 'map' auroral precipitation signatures (VDIS, inverted-V's, etc.) measured by polar orbiting spacecraft in the ionosphere, to the magnetospheric equatorial plane. The combined models allow us to associate latitudinal gradients measured in the ion energy fluxes at low-altitudes with the large-scale pressure gradients in the equatorial plane. Given this global, quasi-static association, we can then make fairly strong statements regarding the location of discrete features in the context of the global picture. We reported on our initial study at national and international meetings and published the results of our predictions of the low-altitude signatures of the plasma sheet. In addition, the PI was invited to contribute a publication to the so-called 'Great Debate in Space Physics' series that is a feature of EOS. The topic was on the nature of magnetospheric substorms. Specific questions of the when and where a substorm occurs and the connection between the auroral and magnetospheric components were discussed in that paper. This paper therefore was derived exclusively from the research supported by this grant. Attachment: Empirical modeling of the quite time nightside magnetosphere.' 'CRRES observations of particle flux dropout event.' The what, where, when, and why of magnetospheric substorm triggers'. and 'Low altitude signature of the plasma sheet: model prediction of local time dependence'.

Description: Two tools for the solution of radiative transfer problems are presented. Streamer is a highly flexible medium spectral resolution radiative transfer model based on the plane-parallel theory of radiative transfer. Capable of computing either fluxes or radiances, it is suitable for studying radiative processes at the surface or within the atmosphere and for the development of remote-sensing algorithms. FluxNet is a fast neural network-based implementation of Streamer for computing surface fluxes. It allows for a sophisticated treatment of radiative processes in the analysis of large data sets and potential integration into geophysical models where computational efficiency is an issue. Documentation and tools for the development of alternative versions of Fluxnet are available. Collectively, Streamer and FluxNet solve a wide variety of problems related to radiative transfer: Streamer provides the detail and sophistication needed to perform basic research on most aspects of complex radiative processes while the efficiency and simplicity of FluxNet make it ideal for operational use.

Description: An orbit tracing technique was used to generate current sheets for three magnetotail models. Groups of ions were followed to calculate the resulting cross-tail current. Several groups then were combined to produce a current sheet. The goal is a model in which the ions and associated electrons carry the electric current distribution needed to generate the magnetic field B in which ion orbits were traced. The region -20 R(E) less than x less than -14 R(E) in geocentric solar magnetospheric coordinates was studied. Emphasis was placed on identifying the categories of ion orbits which contribute most to the cross-tail current and on gaining physical insight into the manner by which the ions carry the observed current distribution. Ions that were trapped near z = 0, ions that magnetically mirrored throughout the current sheet, and ions that mirrored near the Earth all were needed. The current sheet structure was determined primarily by ion magnetization currents. Electrons of the observed energies carried relatively little cross-tail current in these quiet time current sheets. Distribution functions were generated and integrated to evaluate fluid parameters. An earlier model in which B depended only on z produced a consistent current sheet, but it did not provide a realistic representation of the Earth's middle magnetotail. In the present study, B changed substantially in the x and z directions but only weakly in the y direction within our region of interest. Plasmas with three characteristic particle energies were used with each of the magnetic field models. A plasma was found for each model in which the density, average energy, cross-tail current, and bulk flow velocity agreed well with satellite observations.

Description: This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.