ALBERT

Description: The time scales and phases of auroral substorm, activity are quantied in this study using the hemispheric power computed from Polar Ultraviolet Imager (UVI) images. We have applied this technique to several hundred substorm events and we are able to quantify how the characterist act, of substorms vary with season and IMF Bz orientation. We show that substorm time scales vary more strongly with season than with IMF Bz orientation. The recovery time for substorm. activity is well ordered by whether or not the nightside oral zone is sunlit. The recovery time scales for substorms occurring in the winter and equinox periods are similar and are both roughly a factor of two longer than in summer when the auroral oval is sunlit. Our results support the hypothesis that the ionosphere plays an active role in governing the dynamics of the aurora.

Description: Smyer is an H-chondrite impact-melt breccia containing approx.20 vol% 0.5- to 13-mm-thick silicate-rich melt veins surrounding unmelted subrounded chondritic clasts up to 7 cm in maximum dimension. At the interface between some of the melt veins and chondritic clasts, there are troilite-rich regions consisting of unmelted. crushed 0.2- to 140-micron-size angular silicate grains and chondrule fragments surrounded by troilite and transected by thin troilite veins. Troilite fills every available fracture in the silicates. including some as thin as 0.1 microns. Little metallic Fe-Ni is present in these regions: the FeS/Fe modal ratio ranges from -25: 1 to approx.500: 1, far higher than the eutectic weight ratio of 7.5: 1. The texture of these regions indicates that the sulfide formed from a fluid of very low viscosity. The moderately high viscosity (0.2 poise) and large surface tension of liquid FeS, its inability to wet silicate grain surfaces at low oxygen fugacities. and the supereutectic FeS/Fe ratios in the troilite-rich regions indicate that the fluid was a vapor. It seems likely that during the shock event that melted Smyer, many silicates adjacent to the melt veins were crushed. Upon release of shock pressure. some of the troilite evaporated and dissociated. Molecules of S2 were transported and condensed into fractures and around tiny silicate grains: there, they combined with Fe from small adjacent metallic Fe-Ni grains to form troilite. The Ni content at the edges of some of these metal grains increased significantly; Co from these Ni-rich grains diffused into nearby kamacite. Impact-induced S volatilization may have played a major role in depleting the surface of 433 Eros (and other chondritic asteroids) in S.

Description: We conducted gravity wave ray-tracing experiments within an atmospheric region centered near the ARCLITE lidar system at Sondrestrom, Greenland (67N, 310 deg E), in efforts to understand lidar observations of both upper stratospheric gravity wave activity and mesospheric clouds during August 1996 and the summer of 2001. The ray model was used to trace gravity waves through realistic three-dimensional daily-varying background atmospheres in the region, based on forecasts and analyses in the troposphere and stratosphere and climatologies higher up. Reverse ray tracing based on upper stratospheric lidar observations at Sondrestrom was also used to try to objectively identify wave source regions in the troposphere. A source spectrum specified by reverse ray tracing experiments in early August 1996 (when atmospheric flow patterns produced enhanced transmission of waves into the upper stratosphere) yielded model results throughout the remainder of August 1996 that agreed best with the lidar observations. The model also simulated increased vertical group propagation of waves between 40 km and 80 km due to intensifying mean easterlies, which allowed many of the gravity waves observed at 40 km over Sondrestrom to propagate quasi-vertically from 40-80 km and then interact with any mesospheric clouds at 80 km near Sondrestrom, supporting earlier experimentally-inferred correlations between upper stratospheric gravity wave activity and mesospheric cloud backscatter from Sondrestrom lidar observations. A pilot experiment of real-time runs with the model in 2001 using weather forecast data as a low-level background produced less agreement with lidar observations. We believe this is due to limitations in our specified tropospheric source spectrum, the use of climatological winds and temperatures in the upper stratosphere and mesosphere, and missing lidar data from important time periods.

Description: A number of recently published papers suggest that mountain-wave activity in the stratosphere, producing ice particles when temperatures drop below the ice frost point, may be the primary source of large NAT particles. In this paper we use measurements from the Advanced Very High Resolution Radiometer (AVHRR) instruments on board the National Oceanic and Atmospheric Administration (NOAA) polar-orbiting satellites to map out regions of ice clouds produced by stratospheric mountain-wave activity inside the Arctic vortex. Lidar observations from three DC-8 flights in early December 1999 show the presence of solid nitric acid (Type Ia or NAT) polar stratospheric clouds (PSCs). By using back trajectories and superimposing the position maps on the AVHRR cloud imagery products, we show that these observed NAT clouds could not have originated at locations of high-amplitude mountain-wave activity. We also show that mountain-wave PSC climatology data and Mountain Wave Forecast Model 2.0 (MWFM-2) raw hemispheric ray and grid box averaged hemispheric wave temperature amplitude hindcast data from the same time period are in agreement with the AVHRR data. Our results show that ice cloud formation in mountain waves cannot explain how at least three large scale NAT clouds were formed in the stratosphere in early December 1999.

Description: In two separate orbital campaigns (November, 1994 and August, 1997), the Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere (CRISTA) instrument acquired global stratospheric data of high accuracy and high spatial resolution. The standard limb-scanned CRISTA measurements resolved atmospheric spatial structures with vertical dimensions greater than or equal to 1.5 - 2 km and horizontal dimensions is greater than or equal to 100 - 200 km. A fluctuation analysis of horizontal temperature distributions derived from these data is presented. This method is somewhat complementary to conventional power-spectral analysis techniques.

Description: Earth rotation measurements were obtained using Global Positioning System (GPS) data for 11 days during the Epoch '92 campaign in the Summer of 1992. Earth orientation was measured simultaneously with several very long baseline interferornetry (VLBI) networks. These data were processed to yield both GPS and VLBI estimates of UT1 with 3-hour time resolution, which were then compared and analyzed. The high frequency behavior of both data sets is similar, although drifts between the two series of approx.0,1 ms over 2-5 days are evident, Models for tidally induced UT1 variations and estimates of atmospheric angular momentum (AAM) at 6-hour intervals were also compared with the geodetic data, These studies indicate that most of the geodetic signal in the diurnal and semidiurnal frequency bands can be attributed to tidal processes, and that UT1 variations over a few days are mostly atmospheric in origin.

Description: The ability to estimate oxygen fugacities for mare basalts and to extend these observations to the lunar mantle is limited using bulk analysis techniques based on buffering assemblages or the valence state of iron. These limitations are due to reequilibration of mineral assemblages at subsolidus conditions, deviations of mineral compositions from thermodynamic ideality, size requirements, and the limits of the iron valence at very low fO2. Still, these approaches have been helpful and indicate that mare basalts crystallized at fO2 between the iron-w stite buffer (IW) and the ilmenite breakdown reaction (ilmenite = rutile + iron). It has also been inferred from these estimates that the lunar mantle is also highly reduced lying at conditions below IW. Generally, these data cannot be used to determine if the mare basalts become increasingly reduced during transport from their mantle source and eruption at the lunar surface and if there are differences in fO2 among mare basalts or mantle sources. One promising approach to determining the fO2 of mare basalts is using the mean valence of vanadium (2+, 3+, 4+, 5+) determined on spots of a few micrometers in diameter using synchrotron x-ray absorption fine structure (XAFS) spectroscopy. The average valence state of V in basaltic glasses is a function of fO2, temperature, V coordination, and melt composition. Here, we report the initial results of this approach applied to lunar pyroclastic glasses.

Description: Our analysis concerns macro and meso-scale aspects of coupling between the IMF and the magnetosphere-ionosphere system, as opposed to the microphysics of determining how electron gyrotropy is broken and merging actually occurs. We correlate observed behaviors at Cluster and at Polar with temporal variations in other regions, such as in the ionosphere as measured by SuperDARN. Addressing problems with simultaneous observations from diverse locations properly constrains our interpretations.

Description: Carbon dioxide is a low-concentration, but important, component of the Earth's atmosphere. This gas absorbs electromagnetic radiation (EMR) in several regions of the spectrum. Absorption of energy by carbon dioxide adds heat to the atmosphere. In the world today, the burning of fossil fuels and other anthropogenic processes adds carbon dioxide to the atmosphere. Other natural processes in the Earth's system both add and remove carbon dioxide. Overall, measurements of atmospheric carbon dioxide at selected sites around the globe show an increased carbon dioxide concentration in the atmosphere. A figure shows the measured carbon dioxide from Mauna Loa, Hawaii, from 1958 to 2000. Overall, the concentration has increased from 315 to 365 ppm at this site over this period. (There is also a yearly cycle to the concentration that is timed with and hypothesized to be related to the vegetation growing season in the Northern Hemisphere.) The overall expected effect of this increase of atmospheric carbon dioxide is trapping of heat in the atmosphere and global warming. While this overall relationship between carbon dioxide and global warming seems straightforward, many of the specific details relating to regional and local sources and sinks and gradients of carbon dioxide are not well understood. A remote sensing capability to measure carbon dioxide could provide important inputs for scientific research to better understand the distribution and change in atmospheric carbon dioxide at detailed spatial and temporal levels. In pursuit of this remote sensing of carbon dioxide objective, this paper analyzes the expression of carbon dioxide in the spectral range measured by the Airborne Visible/Infrared Imagery Spectrometer (AVIRIS). Based on these analyses, a spectral-fitting algorithm that uses AVIRIS measured spectra and MODTRAN radiative-transfer code modeled spectra to derive total column carbon dioxide abundance has been developed. This algorithm has been applied to an AVIRIS data set acquired over Pasadena, California, in 1999 and a data set acquired over the Pacific Ocean near Hawaii in 2000 with promising results. This is ongoing research; the current initial analyses, measurements, and results are reported in this paper.

Description: Second quarter progress report for year 3 on contract NASW-99002 "What is the Relationship between heavy ion outflow and high latitude energetic particle precipitation". In this Project we are studying the relationship between the fluxes, mean energies, and field-aligned flow speeds of escaping suprathermal H+ and 0+ measured by the TEAMS instrument on FAST and the energy flux of precipitating electrons obtained from the LBHL images taken by the UVI camera on POLAR. In this portion of the project we are using UVI images to tell us when substorm onsets occur and how the auroral zone changes during the course of a substorm. We are correlating this information with TEAMS flux measurements made over the auroral zone at times close to these substorms. The goal is to understand how the flux of suprathermal ion outflow varies with substorm phase.

Description: We review solar geophysical data relating to the great magnetic storm of 14-15 May 1921, with emphasis on observations of the low-latitude visual aurora. From the reports we have gathered for this event the lowest geomagnetic latitude of definite overhead aurora (coronal form) was 40 deg and the lowest geomagnetic latitude from which auroras were observed on the poleward horizon in the northern hemisphere was 30 deg. For comparison, corresponding overhead/low-latitude values of 48 deg/32 deg and 41 deg/20 deg were reported for the great auroras on 28-29 August and 1-2 September 1859, respectively. However for the 1921 event, there is a report of aurora from Apia, Samoa, in the southern hemisphere, within 13 deg of the geomagnetic equator. This report by professional observers appears to be credible, based on the aurora description and timing, but is puzzling because of the discrepancy with the lowest latitude of observation in the northern hemisphere and the great implied aurora height (approximately 2000 km, assuming overhead aurora at Auckland, New Zealand). We discuss various possibilities that might account for this observation.

Description: During June to September 1992, a special campaign was held to measure rapid changes in Earth's rotation rate and to relate these measurements to variations in the atmosphere's angular momentum, due principally to changes in zonal winds. A strong rise in both length of day and atmospheric momentum during a particular 6-day subperiod is documented, and this example of a short-period perturbation is identified with a specific regional coupling mechanism. Mountain torques within the southern tropics appear to account for most of the rapid momentum transfer between the solid Earth and atmosphere, with those across South America especially important.

Description: In situ measurements of hydrogen, nitrogen, and chlorine radicals obtained through sunrise and sunset in the lower stratosphere during SPADE are compared to results from a photochemical model constrained by observed concentrations of radical precursors and environmental conditions. Models allowing for heterogeneous hydrolysis of N205 on sulfate aerosols agree with measured concentrations of NO, NO2, and ClO throughout the day, but fail to account for high concentrations of OH and H02 observed near sunrise and sunset. The morning burst of [OH] and [HO2] coincides with the rise of [NO] from photolysis of N02, suggesting a new source of HO, that photolyzes in the near UV (350 to 400 nm) spectral region. A model that allows for the heterogeneous production of HN02 results in an excellent simulation of the diurnal variations of [OH] and [HO2].

Description: IAEMIS (Integrated Automated Emergency Management Information System) is the principal tool of an earthquake preparedness program developed by Martin Marietta and the Mid-America Remote Sensing Center (MARC). It is a two-component set of software, data and procedures to provide information enabling management personnel to make informed decisions in disaster situations. The NASA-developed program ELAS, originally used to analyze Landsat data, provides MARC with a spatially-oriented information management system. Additional MARC projects include land resources management, and development of socioeconomic data.

Description: Imaging spectroscopy offers a framework based in physics and chemistry for scientific investigation of a wide range of phenomena of interest in the Earth environment. In the scientific discipline of volcanology knowledge of lava temperature and distribution at the surface provides insight into the volcano status and subsurface processes. A remote sensing strategy to measure surface lava temperatures and distribution would support volcanology research. Hot targets such as molten lava emit spectral radiance as a function of temperature. A figure shows a series of Planck functions calculated radiance spectra for hot targets at different temperatures. A maximum Lambertian solar reflected radiance spectrum is shown as well. While similar in form, each hot target spectrum has a unique spectral shape and is distinct from the solar reflected radiance spectrum. Based on this temperature-dependent signature, imaging spectroscopy provides an innovative approach for the remote-sensing-based measurement of lava temperature. A natural site for investigation of the measurement of lava temperature is the Big Island of Hawaii where molten lava from the Kilauea vent is present at the surface. In the past, Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data sets have been used for the analysis of hot volcanic targets and hot burning fires. The research presented here builds upon and extends this earlier work. The year 2000 Hawaii AVIRIS data set has been analyzed to derive lava temperatures taking into account factors of fractional fill, solar reflected radiance, and atmospheric attenuation of the surface emitted radiance. The measurements, analyses, and current results for this research are presented here.

Description: The formation mechanism of terrestrial planetary is still poorly understood, and has been the subject of numerous experimental studies. Several mechanisms have been proposed by which metal-mainly iron with some nickel-could have been extracted from a silicate mantle to form the core. Most recent models involve gravitational sinking of molten metal or metal sulphide through a partially or fully molten mantle that is often referred to as a'magma ocean. Alternative models invoke percolation of molten metal along an interconnected network (that is, porous flow) through a solid silicate matrix. But experimental studies performed at high pressures have shown that, under hydrostatic conditions, these melts do not form an interconnected network, leading to the widespread assumption that formation of metallic cores requires a magma ocean. In contrast, here we present experiments which demonstrate that shear deformation to large strains can interconnect a significant fraction of initially isolated pockets of metal and metal sulphide melts in a solid matrix of polycrystalline olivine. Therefore, in a dynamic (nonhydrostatic) environment, percolation remains a viable mechanism for the segregation and migration of core-forming melts in a solid silicate mantle.

Description: We use the first simultaneous in situ measurements of ClONO2, ClO, and HCl acquired using the NASA ER-2 aircraft during the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) mission to test whether these three compounds quantitatively account for total inorganic chlorine (Cly) in the lower stratosphere in 1997. We find (ClO + ClONO2 + HCl)/Cly = 0.92 +/- 0.10, where Cly is inferred from in situ measurements of organic chlorine source gases. These observations are consistent with our current understanding of the budget and partitioning of Cly in the lower stratosphere. We find no evidence in support of missing inorganic chlorine species that compose a significant fraction of Cly. We apply the analysis to earlier ER-2 observations dating from 1991 to investigate possible causes of previously observed discrepancies in the inorganic chlorine budget. Using space shuttle, satellite, balloon, and aircraft measurements in combination with ER-2 data, we find that the discrepancy is unlikely to have been caused by missing chlorine species or an error in the photolysis rate of chlorine nitrate. We also find that HCl/Cly is not significantly controlled by aerosol surface area density in the lower stratosphere.

Description: We examine inorganic chlorine (Cly) partitioning in the summer lower stratosphere using in situ ER-2 aircraft observations made during the Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) campaign. New steady state and numerical models estimate [ClONO2]/[HCl] using currently accepted photochemistry. These models are tightly constrained by observations with OH (parameterized as a function of solar zenith angle) substituting for modeled HO2 chemistry. We find that inorganic chlorine photochemistry alone overestimates observed [ClONO2]/[HCl] by approximately 55-60% at mid and high latitudes. On the basis of POLARIS studies of the inorganic chlorine budget, [ClO]/[ClONO2], and an intercomparison with balloon observations, the most direct explanation for the model-measurement discrepancy in Cly partitioning is an error in the reactions, rate constants, and measured species concentrations linking HCl and ClO (simulated [ClO]/[HCl] too high) in combination with a possible systematic error in the ER-2 ClONO2 measurement (too low). The high precision of our simulation (+/-15% 1-sigma for [ClONO2]/[HCl], which is compared with observations) increases confidence in the observations, photolysis calculations, and laboratory rate constants. These results, along with other findings, should lead to improvements in both the accuracy and precision of stratospheric photochemical models.

Description: The consequences of electric field acceleration and an inhomogencous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one- dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogencous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of 0(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function in investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogencous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Description: The Lidar In-space Technology Experiment (LITE) is a multi-wavelength backscatter lidar developed by NASA Langley Research Center to fly on the Space Shuttle. The LITE instrument is built around a three-wavelength ND:YAG laser and a 1-meter diameter telescope. The laser operates at 10 Hz and produces about 500 mJ per pulse at 1064 nm and 532 nm, and 150 mJ per pulse at 355 nm. The objective of the LITE program is to develop the engineering processes required for space lidar and to demonstrate applications of space-based lidar to remote sensing of the atmosphere. The LITE instrument was designed to study a wide range of cloud and aerosol phenomena. To this end, a comprehensive program of scientific investigations has been planned for the upcoming mission. Simulations of on-orbit performance show the instrument has sufficient sensitivity to detect even thin cirrus on a single-shot basis. Signal averaging provides the capability of measuring the height and structure of the planetary boundary layer, aerosols in the free troposphere, the stratospheric aerosol layer, and density profiles to an altitude of 40 km. The instrument has successfully completed a ground-test phase and is scheduled to fly on the Space Shuttle Discovery for a 9-day mission in September 1994.

Description: Improved techniques have been developed for empirical modeling of the high-latitude electric potentials and magnetic field aligned currents (FAC) as a function of the solar wind parameters. The FAC model is constructed using scalar magnetic Euler potentials, and functions as a twin to the electric potential model. The improved models have more accurate field values as well as more accurate boundary locations. Non-linear saturation effects in the solar wind-magnetosphere coupling are also better reproduced. The models are constructed using a hybrid technique, which has spherical harmonic functions only within a small area at the pole. At lower latitudes the potentials are constructed from multiple Fourier series functions of longitude, at discrete latitudinal steps. It is shown that the two models can be used together in order to calculate the total Poynting flux and Joule heating in the ionosphere. An additional model of the ionospheric conductivity is not required in order to obtain the ionospheric currents and Joule heating, as the conductivity variations as a function of the solar inclination are implicitly contained within the FAC model's data. The models outputs are shown for various input conditions, as well as compared with satellite measurements. The calculations of the total Joule heating are compared with results obtained by the inversion of ground-based magnetometer measurements. Like their predecessors, these empirical models should continue to be a useful research and forecast tools.

Description: During the 1997 Photochemistry of Ozone Loss in the Arctic Region in Summer (POLARIS) mission, simultaneous in situ observations of NOx and HOx radicals, their precursors, and the radiation field were obtained in the lower stratosphere. We use these observations to evaluate the primary mechanisms that control NOx-HNO3 exchange and to understand their control over the partitioning between NO2 and HNO3 in regions of continuous sunlight. We calculate NOx production (PNOx) and loss (LNOx) in a manner directly constrained by the in situ measurements and current rate constant recommendations, using approaches for representing albedo, overhead O3 and [OH] that reduce model uncertainty. We find a consistent discrepancy of 18% between modeled rates of NOx production and loss (LNOx = 1.18P(sub NOx)), which is within the measurement uncertainty of +/- 27%. The partitioning between NOx production processes is [HNO3 + OH (41 +/- 2)%; HNO3 + hv (59 +/- 2)%] and between NOx loss processes is [NO2 + OH, 90% to 〉97%; BrONO2 + H2O, 10% to 〈3%]. The steady-state description of NOx-HNO3 exchange reveals the significant influence of the tight correlation between the photolysis rate of HNO3 and [OH] established by in situ measurements throughout the lower stratosphere. Parametrizing this relationship, we find: (1) the steady-state value of [NO2](sub 24h-avg)/[HNO3] in the continuously sunlit, lower stratosphere is a function only of temperature and number density; and (2) the partitioning of NOx production between HNO3 + OH and HNO3 + hv is nearly constant throughout most of the lower stratosphere. We describe a methodology (functions of latitude, day, temperature, and pressure) for accurately predicting the steady-state value of [NO2](sub 24h-avg)/[HNO3] and the partitioning of NOx production within these regions. The results establish a metric to compare observations of [NO2](sub 24h-avg)/[HNO3] within the continuously sunlit region and provide a simple diagnostic for evaluating the accuracy of models that attempt to describe the coupled NOx-HOx photochemistry in the lower stratosphere.

Description: This grant is an extension to our previous NASA Grant NAG5-3461, providing incremental funding to continue GOME (Global Ozone Monitoring Experiment) and SCIAMACHY (SCanning Imaging Absorption SpectroMeter for Atmospheric CHartographY) studies. This report summarizes research done under these grants through December 31, 2002. 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, consultation and development for SCIAMACHY near-real-time (NRT) and off-line (OL) data products, and participation in initial SCIAMACHY validation studies. The Global Ozone Monitoring Experiment was successfully launched on the ERS-2 satellite on April 20, 1995, and remains working in normal fashion. SCIAMACHY was launched March 1, 2002 on the ESA Envisat satellite. Three GOME-2 instruments are now scheduled to fly on the Metop series of operational meteorological satellites (Eumetsat). K. Chance is a member of the reconstituted GOME Scientific Advisory Group, which will guide the GOME-2 program as well as the continuing ERS-2 GOME program.

Description: Eighteen long-range flights over the Pacific Ocean between 38 S to 20 N and 166 E to 90 W were made by the NASA DC-8 aircraft during the NASA Pacific Exploratory Mission (PEM) Tropics B conducted from March 6 to April 18, 1999. Two lidar systems were flown on the DC-8 to remotely measure vertical profiles of ozone (O3), water vapor (H2O), aerosols, and clouds from near the surface to the upper troposphere along their flight track. In situ measurements of a wide range of gases and aerosols were made on the DC-8 for comprehensive characterization of the air and for correlation with the lidar remote measurements. The transition from northeasterly flow of Northern Hemispheric (NH) air on the northern side of the Intertropical Convergence Zone (ITCZ) to generally easterly flow of Southern Hemispheric (SH) air south of the ITCZ was accompanied by a significant decrease in O3, carbon monoxide, hydrocarbons, and aerosols and an increase in H2O. Trajectory analyses indicate that air north of the ITCZ came from Asia and/or the United States, while the air south of the ITCZ had a long residence time over the Pacific, perhaps originating over South America several weeks earlier. Air south of the South Pacific Convergence Zone (SPCZ) came rapidly from the west originating over Australia or Africa. This air had enhanced O3 and aerosols and an associated decrease in H2O. Average latitudinal and longitudinal distributions of O3 and H2O were constructed from the remote and in situ O3 and H2O data, and these distributions are compared with results from PEM-Tropics A conducted in August-October 1996. During PEM-Tropics B, low O3 air was found in the SH across the entire Pacific Basin at low latitudes. This was in strong contrast to the photochemically enhanced O3 levels found across the central and eastern Pacific low latitudes during PEM-Tropics A. Nine air mass types were identified for PEM-Tropics B based on their O3, aerosols, clouds, and potential vorticity characteristics. The data from each flight were binned by altitude according to air mass type, and these results showed the relative observational frequency of the different air masses as a function of altitude in seven regions over the Pacific. The average chemical composition of the major air mass types was determined from in situ measurements in the NH and SH, and these results provided insight into the origin, lifetime, and chemistry of the air in these regions.

Description: Accurate determination of spatial soil moisture distribution and monitoring its temporal variation have a significant impact on the outcomes of hydrologic, ecologic, and climatic models. Development of a successful remote sensing instrument for soil moisture relies on the accurate knowledge of the soil dielectric constant (epsilon(sub soil)) to its moisture content. Two existing methods for measurement of dielectric constant of soil at low and high frequencies are, respectively, the time domain reflectometry and the reflection coefficient measurement using an open-ended coaxial probe. The major shortcoming of these methods is the lack of accurate determination of the imaginary part of epsilon(sub soil). In this paper a microstrip ring resonator is proposed for the accurate measurement of soil dielectric constant. In this technique the microstrip ring resonator is placed in contact with soil medium and the real and imaginary parts of epsilon(sub soil) are determined from the changes in the resonant frequency and the quality factor of the resonator respectively. The solution of the electromagnetic problem is obtained using a hybrid approach based on the method of moments solution of the quasi-static formulation in conjunction with experimental data obtained from reference dielectric samples. Also a simple inversion algorithm for epsilon(sub soil) = epsilon'(sub r) + j(epsilon"(sub r)) based on regression analysis is obtained. It is shown that the wide dynamic range of the measured quantities provides excellent accuracy in the dielectric constant measurement. A prototype microstrip ring resonator at L-band is designed and measurements of soil with different moisture contents are presented and compared with other approaches.

Description: In the fall of 1997, during an Intensive Observation Period (IOP), the Atmospheric Radiation Measurement (ARM) program conducted a study of water vapor abundance measurement at its Southern Great Plains (SGP) site. Among a large number of instruments, four sun-tracking radiometers were present to measure the columnar water vapor (CWV). All four solar radiometers retrieve CWV by measuring total solar transmittance in the 0.94-gm water vapor absorption band and subtracting contributions due to Rayleigh, ozone and aerosol transmittances. The aerosol optical depth comparisons among the same four radiometers has been presented elsewhere (Geophys. Res. Lett., 26, 17, 2725-2728, 1999). We have used three different methods to retrieve CWV. In a first round of comparison no attempt was made to standardize on the same radiative transfer model and its underlying water vapor spectroscopy. In the second round of comparison we used the same line-by-line code (which includes recently corrected H2O spectroscopy) to retrieve CAN from all four suntracking radiometers. This decreased the mean CWV by 8% or 13%. The spread of 8% in the solar radiometer results found when using the same model is an indication of the other-than-model uncertainties involved in determining CWV from solar transmittance measurements with current instrumentation.

Description: Trajectory hunting (i.e., a technique to find air parcels sampled at least twice over the course of a few days) is applied to analyze Upper Atmosphere Research Satellite (UARS) measurements in conjunction with the Atmospheric and Environmental Research, Inc. (AER) photochemical box model. As a case study, we investigate rapid chlorine activation in the Arctic lower stratosphere on December 29, 1992 associated with a polar stratospheric cloud (PSC) event. Eleven air parcels that have been sampled several times along five-day trajectories at the 465 K (approx. 46 hPa), 520 K (approx. 31 hPa), and 585 K (approx. 22 hPa) levels were investigated. For the first time, the latest versions of the Cryogenic Limb Array Etalon Spectrometer (CLAES, version 9) and Microwave Limb Sounder (MLS, version 5) data sets are analyzed, and their consistency is assessed. A detailed sensitivity study with the AER photochemical box model along these trajectories leads to the conclusion that for the December 24-29, 1992 episode: (1) the individual CLAES ClONO2 and MLS ClO measurements are self-consistent within their uncertainties; and (2) most of the time, UARS measurements of ClO, ClONO2, HNO3, and aerosol extinction at 780 cm(exp -1) agree within the range of their uncertainties with the model calculations. It appears that the HNO3 and aerosol extinction measurements for four parcels at 520 K look more supportive for the nitric acid trihydrate (NAT) scheme, However, the uncertainties in the individual UARS measurements and UK Meteorological Office temperature do not allow a definite discrimination between the NAT and supercooled ternary solution (STS) PSC schemes for this chlorine activation episode in December 1992.

Description: Sonoluminescence is the term used to describe the emission of light from a violently collapsing bubble. Sonoluminescence ("light from sound") is the result of extremely nonlinear pulsations of gas/vapor bubbles in liquids when subject to sufficiently high amplitude acoustic pressures. In a single collapse, a bubble's volume can be compressed more than a thousand-fold in the span of less than a microsecond. Even the simplest consideration of the thermodynamics yields pressures on the order of 10,000 ATM, and temperatures of at least 10,000K. On the face of things, it is not surprising that light should be emitted from such an extreme process. Since 1990 (the year that Gaitan discovered light from a single bubble) there has been a tremendous amount of experimental and theoretical research in stable, single-bubble sonoluminescence (SBSL), yet there remain at least four unexplained phenomena associated with SBSL in 1g: the light emission mechanism itself, the existence of anisotropies in the emitted light, the disappearance of the bubble at some critical acoustic pressure, and the appearance of quasiperiodic and chaotic oscillations in the flash timing. Gravity, in the context of the buoyant force, is implicated in all four of these. We are developing KC-135 experiments probing the effect of gravity on single bubble sonoluminescence. By determining the stability boundaries experimentally in microgravity, and measuring not only light emission but mechanical bubble response, we will be able to directly test the predictions of existing theories.

Description: To demonstrate the potential scientific and technical merits of in situ microscopy on Mars, we analyzed a possible Martian regolith analog - an acolian red dune sand from the central Australian desert (near Mt. Olga). This sand was chosen for its ubiquitous red coating and the desert environment in which is it found. Grains of this sand were analyzed using a variety of microanalytical techniques. A database of detailed studies of such terrestrial analogs would assist the study of geological and astrobiological specimens in future missions to Mars. Potential instrument concepts for in situ deployment on Mars include local electrode atom probe nanoanalysis (LEAP), vertical scanning white light interferometry (VSWLI), scanning electron microscopies, energy dispersive x-ray microanalysis (EDX), atomic force microscopy (AFM) and X-ray diffraction (XRD). While in situ deployment of these techniques is many years away, ground-based studies using these analytical techniques extend our understanding of the data obtained from instruments to be flown in the near future.

Description: The overall objective of this grant was to: 1) Quantify thermospheric neutral wind behavior in the ionosphere. This was to be achieved by developing an improved empirical wind model. 2) Validating the procedure for obtaining winds from the height of the peak density. 3) Improving the model capabilities and making updated versions of the model available to other scientists. The approach is to use neutral winds derived from ionosonde measurements of the height of the peak electron density (h(sub m)F(sub 2)). One of the proposed first year tasks was to perform some validation studies on the method. Substantial progress has been made with regard to both the empirical model and the validation study. Funding from this grant has also enabled a number of fruitful collaborations with other researchers; one of the stated aims in the proposal. Graduate student Mayra Martinez has developed the mathematical formulation for the empirical wind model as part of her dissertation. As proposed, authors continued validation studies of the technique for determining winds from h(sub m)F(sub 2). They are submitted a paper to the Journal of Geophysical Research in December 1996 entitled "Therinospheric neutral winds at southern mid-latitudes: comparison of optical and ionosonde h(sub m)F(sub 2) methods. A second paper entitled "Ionospheric behavior at a southern mid-latitude in March 1995" has come out of the March 1995 data set and was published in The Journal of Geophysical Research. A new algorithm was developed. The ionosphere also have been modeled.

Description: The purpose of this grant was to study the electrical current flow along and across the magnetic field lines in the middle and outer magnetosphere. These currents are both time stationary and fluctuating. Our study of the time stationary currents formed part of the Masters dissertation of graduate student Tom Meseroll. In the bibliography below we list the papers in journals and at meetings arising from these studies. In addition to these studies, we were asked to assist in several studies of the currents in the magnetotail, in the dayside magnetosphere, and near the polar cusp (with Hawkeye data). Also we undertook a study of the currents associated with sudden impulses caused by sudden solar wind pressure changes. Below we say a few words about these studies. This is followed by a bibliography listing the papers. There were no patents issued or inventions associated with this study.

Description: This Handbook summarizes the radiation data from the Soviet 'North Pole' drifting research stations operated in the Arctic from 1950 to 1991. The Handbook contains reduced mean monthly values of the fluxes of solar radiation for solar elevations up to 35 deg at 5 deg intervals, monthly totals of the net radiation and its constituents, and mean monthly values of the albedo of the sea-ice surfaces. The Handbook presents information of interest to specialists studying the physics of the atmosphere, climatology, geography, and related disciplines.

Description: We describe a method for retrieving neutral thermospheric composition and solar EUV flux from optical measurements of the O(+)(P-2) 732 nm and O(D-1) 630 nm airglow emissions. The parameters retrieved are the neutral temperature, the O, L2, and N2 density profiles, and a scaling factor for the solar EUV flux spectrum. The temperature, solar EUV flux scaling factor, and atomic oxygen density are first retrieved from the 732 nm emission, which are then used with the 630 nm emission to retrieve the O2 and N2 densities. Between the altitudes of 200 and 400 km the retrieval technique is able to statistically retrieve values to within 3.1% for thermospheric temperature, 3.3% for atomic oxygen, 2.3% for molecular oxygen, and 2.4% for molecular nitrogen. The solar EUV flux scaling factor has a retrieval error of 5.1%. We also present the results of retrievals using existing data taken from both groundbased and spacebased instruments. These include airglow data taken by the Visible Airglow Experiment on the Atmospheric Explorer spacecraft and the Imaging Spectrometric Observatory flown on the ATLAS 1 shuttle mission in 1992.

Description: The invention is a method for generating an image of a parametric surface, such as the compass direction toward which each surface element of terrain faces, commonly called the slope-aspect azimuth of the surface element. The method maximizes color contrast to permit easy discrimination of the magnitude, ranges, intervals or classes of a surface parameter while making it easy for the user to visualize the form of the surface, such as a landscape. The four pole colors of the opponent process color theory are utilized to represent intervals or classes at 90 degree angles. The color perceived as having maximum measured luminance is selected to portray the color having an azimuth of an assumed light source and the color showing minimum measured luminance portrays the diametrically opposite azimuth. The 90 degree intermediate azimuths are portrayed by unique colors of intermediate measured luminance, such as red and green. Colors between these four pole colors are used which are perceived as mixtures or combinations of their bounding colors and are arranged progressively between their bounding colors to have perceived proportional mixtures of the bounding colors which are proportional to the interval's angular distance from its bounding colors.

Description: This research program investigated changes of dynamical and chemical structure of the lower stratosphere and how they are related to elements of the tropospheric general circulation. These considerations were explored in total ozone data from TOMS on board the Nimbus-7 satellite. It was shown that most of the daily variance of total ozone was accounted for by quasi-horizontal transport of ozone along isentropic surfaces in the lower stratosphere. Air descending along theta surfaces experiences compression that increases the local ozone number density and column abundance. Just the reverse is experienced by air ascending along isentropics surfaces. Together, these mechanisms provide an explanation for ozone "mini-hole" phenomenon, which punctuates the circulation of the Southern Hemisphere.

Description: Previous studies have demonstrated the usefulness of major and minor elements in silicate phases to understand differences among basaltic systems and the influence of different planetary environments on basalt chemistry (e.g., Papike [1]). Intriguing data displays presented by Papike [1] include a plot of Mn vs. Fe (atoms per formula unit, afu) for pyroxene and olivine and a plot showing the anorthite content of plagioclase from different planetary basalts. Here we combine portions of these plots (Fig. 4) and provide all new data for olivine and plagioclase.

Description: THz spectroscopy of the atmosphere has been driven by the need to make remote sensing measurements of OH. While the THz region can be used for sensitive detection on many atmospheric molecules, the THz region is the best region for measuring the diurnal behavior of stratospheric OH by remote sensing. The infrared region near 3 microns suffers from chemiluminescence and from spectral contamination due to water. The ultraviolet region near 300 nm requires solar illumination. The three techniques for OH emission measurements in the THz region include Fourier Transform interferometry, Fabry-Perot interferometry, and heterodyne radiometry. The first two use cryogenic direct detectors while the last technique uses a local oscillator and a mixer to down convert the THz signal to GHz frequencies. All techniques have been used to measure stratospheric OH from balloon platforms. OH results from the Fabry-Perot based FILOS instrument will be given. Heterodyne measurement of OH at 2.5 THz has been selected to be a component of the Microwave Limb Sounder on the Earth Observing System CHEM-1 polar satellite. The design of this instrument will be described. A balloon-based prototype heterodyne 2.5 THz radiometer had its first flight on, 24 May 1998. Results form this flight will be presented.

Description: This workshop, held in April 2002, brought together various Earth Sciences experts to focus on the subseasonal prediction problem. While substantial advances have occurred over the last few decades in both weather and seasonal prediction, progress in improving predictions on these intermediate time scales (time scales ranging from about two weeks to two months) has been slow. The goals of the workshop were to get an assessment of the "state of the art" in predictive skill on these time scales, to determine the potential sources of "untapped" predictive skill, and to make recommendations for a course of action that will accelerate progress in this area. One of the key conclusions of the workshop was that there is compelling evidence for predictability at forecast lead times substantially longer than two weeks. Tropical diabatic heating and soil wetness were singled out as particularly important processes affecting predictability on these time scales. Predictability was also linked to various low-frequency atmospheric "phenomena" such as the annular modes in high latitudes (including their connections to the stratosphere), the Pacific/North American (PNA) pattern, and the Madden Julian Oscillation (MJO). The latter, in particular, was highlighted as a key source of untapped predictability in the tropics and subtropics, including the Asian and Australian monsoon regions.

Description: Using measurements from the High Altitude Plasma Instrument (HAPI) on the Dynamics-Explorer 1 (DE-1) spacecraft and the Low Altitude Plasma Instrument (LAPI) on Dynamics Explorer 2 (DE 2), we investigate both die high altitude and low altitude extents of the auroral acceleration region. To infer the high altitude limit, we searched the HAPI data base for evidence of upward-directed auroral electric fields located above the spacecraft when the HAPI spacecraft is above 9000 km altitude. We find that such acceleration is common when DE-1 flies through die auroral oval at an altitude of 9,000-11,000 km. At altitudes above 11,000 km, the fraction of the orbits with evidence of at least a 1000 V potential drop above the spacecraft falls, becoming essentially zero above an altitude of 15,000 km. Above that altitude, small (100 V) potential drops are frequently observed, but only rarely are approx. 1 kV potentials observed, typically associated with polar cap or 'theta' arcs or westward traveling surges. To investigate the low-altitude limit of the auroral acceleration region, we use conjunctions of DE 1 and DE 2 along auroral field lines and match the upgoing fluxes of ionospheric ions observed by DE 2 with the flux of accelerated upgoing ions observed at DE 1. Calculating the ionospheric scale height from the ion and electron temperatures and assuming that the parallel flow velocity is independent of height above 800 km, we calculate the altitude at which the upwelling ionospheric ions are effectively completely lost to upward acceleration. The initial lowest-altitude acceleration process could be either a perpendicular acceleration or a parallel electric field, but it must be sufficient to give the entire distribution escape energy. We find that in the two cases studied, near the region of peak auroral potential drop the altitude of this acceleration was around 1700 km (near the O/H neutral crossover altitude), but was significantly higher (approx. 2000 km) near the edges of the arc, where the potential was lower. The composition of the upgoing ion beam was consistent with these heights, being predominately H(+) near the edges and O(+) near the peak.

Description: During the last decade there has been an unprecedented improvement in both the accuracy, the temporal resolution of Earth's rotation measurements. Determination of the position of the Earth's rotation axis both in inertial space and with respect to the crust with accuracies of about 0.3 milliarcseconds (mas) are now routine. In recent years, there has been and emphasis on the determination of short-period (daily and less) variations in Earth rotation. Two space based geodetic systems, very long baseline interferometry (VLBI) and the global positioning system (GPS) have proved to be very successful in this endeavor. Results for the tidally coherent part of the subdaily Earth rotation variations determined from the analysis of VLBI data are discussed. The magnitude of other subdaily variations are also considered.

Description: Measurements from the Halogen Occultation Experiment (HALOE) revealed the infrared signature of polar mesospheric clouds (PMCs), for the first time, HALOE PMC observations at eight wavelengths (2.45 - 10 microns) show remarkable agreement with model PMC spectra based on ice particle extinction, and thus provide the first confirmation that water ice is the primary component of PMCs. Because PMCs respond to changes in temperature and water vapor, they are considered an indicator of global climate change. We propose to further the understanding of PMCs using a decade of infrared measurements form HALOE. This effort will characterize PMC spectral properties, extinction profiles, and size distributions. Using this information, HALOE measurements will be used to make simultaneous retrievals of H2O, O3, and temperature, in the presence of PMCs. The simultaneous retrievals of particle properties, H2O, O3, and temperature will be used with HALOE NO data to provide a significant step forward in the knowledge of PMC characteristics and formation conditions. We will challenge fundamental theories of PMC formation, and investigate changes in PMC properties and related conditions over the length of the HALOE measurement record. HALOE has been operating without flaw since it was launched on October 11, 1991. Consequently, ten southern and ten northern PMC seasons have been observed thus far, providing a wealth of data for the study of PMCs and related parameters.

Description: This project involved analyses of atmospheric constituent data fields, carbon monoxide in the upper stratospheric/lower mesosphere, and water vapor in the upper troposphere. The observational data analyses were compared with atmospheric models.

Description: The specific objectives of this research effort included the following: 1) Quantification of gravity wave propagation throughout the lower and middle atmosphere in order to define the roles of topographic and convective sources and filtering by mean and low-frequency winds in defining the wave field and wave fluxes at greater altitudes; 2) The influences of wave instability processes in constraining wave amplitudes and fluxes and generating turbulence and transport; 3) Gravity wave forcing of the mean circulation and thermal structure in the presence of variable motion fields and wave-wave interactions, since the mean forcing may be a small residual when wave interactions, anisotropy, and momentum and heat fluxes are large; 4) The statistical forcing and variability imposed on the thermosphere at greater altitudes by the strong wave forcing and interactions occurring in the MLTI.

Description: This report summarizes our investigations performed under NASA Grant NAG5-8058. The three-year research supported by the Geospace Sciences SR&T program (Ionospheric, Thermospheric, and Mesospheric Physics) has been designed to investigate fluxes of energetic oxygen and nitrogen atoms in the terrestrial thermosphere. Fast atoms are produced due to absorption of the solar radiation and due to coupling between the ionosphere and the neutral thermospheric gas. We have investigated the impact of hot oxygen and nitrogen atoms on the thermal balance, chemistry and radiation properties of the terrestrial thermosphere. Our calculations have been focused on the accurate quantitative description of the thermalization of O and N energetic atoms in collisions with atom and molecules of the ambient neutral gas. Upward fluxes of oxygen and nitrogen atoms, the rate of atmospheric heating by hot oxygen atoms, and the energy input into translational and rotational-vibrational degrees of atmospheric molecules have been evaluated. Altitude profiles of hot oxygen and nitrogen atoms have been analyzed and compared with available observational data. Energetic oxygen atoms in the terrestrial atmosphere have been investigated for decades, but insufficient information on the kinetics of fast atmospheric atoms has been a main obstacle for the interpretation of observational data and modeling of the hot geocorona. The recent development of accurate computational methods of the collisional kinetics is seen as an important step in the quantitative description of hot atoms in the thermosphere. Modeling of relaxation processes in the terrestrial atmosphere has incorporated data of recent observations, and theoretical predictions have been tested by new laboratory measurements.

Description: Although the exact site for the origin of the r-process isotopes remains mysterious, most thinking has centered on matter ejected from the cores of massive stars in core-collapse supernovae [13]. In the 1970's and 1980's, however, difficulties in understanding the yields from such models led workers to consider the possibility of r-process nucleosynthesis farther out in the exploding star, in particular, in the helium burning shell [4,5]. The essential idea was that shock passage through this shell would heat and compress this material to the point that the reactions 13C(alpha; n)16O and, especially, 22Ne(alpha; n)25Mg would generate enough neutrons to capture on preexisting seed nuclei and drive an "n process" [6], which could reproduce the r-process abundances. Subsequent work showed that the required 13C and 22Ne abundances were too large compared to the amounts available in realistic models [7] and recent thinking has returned to supernova core material or matter ejected from neutron star-neutron star collisions as the more likely r-process sites.

Description: We compared the version 5 Microwave Limb Sounder (MLS) aboard the Upper Atmosphere Research Satellite (UARS), version 3 Polar Ozone and Aerosol Measurement-III (POAM-111) aboard the French satellite SPOT-IV, version 6.0 Stratospheric Aerosol and Gas Experiment 11 (SAGE-II) aboard the Earth Radiation Budget Satellite, and NASA ER-2 aircraft measurements made in the northern hemisphere in January-February 2000 during the SAGE III Ozone Loss and Validation Experiment (SOLVE). This study addresses one of the key scientific objectives of the SOLVE campaign, namely, to validate multi-platform satellite measurements made in the polar stratosphere during winter. This intercomparison was performed using a traditional correlative analysis (TCA) and a trajectory hunting technique (THT). Launching backward and forward trajectories from the points of measurement, the THT identifies air parcels sampled at least twice within a prescribed match criterion during the course of 5 days. We found that the ozone measurements made by these four instruments agree most of the time within 110% in the stratosphere up to 1400 K (approximately 35 km). The water vapor measurements from POAM-III and the ER-2 Harvard Lyman-alpha hygrometer and JPL laser hygrometer agree to within 10.5 ppmv (or about +/-10%) in the lower stratosphere above 380 K. The MLS and ER-2 ClO measurements agree within their error bars for the TCA. The MLS and ER-2 nitric acid measurements near 17-20 km altitude agree within their uncertainties most of the time with a hint of a positive offset by MLS according to the TCA. We also applied the AER box model constrained by the ER-2 measurements for analysis of the ClO and HN03 measurements using the THT. We found that: (1) the model values of ClO are smaller by about 0.3-0.4 (0.2) ppbv below (above) 400 K than those by MLS and (2) the HN03 comparison shows a positive offset of MLS values by approximately 1 and 1-2 ppbv below 400 K and near 450 K, respectively. It is hard to quantify the HN03 offset in the 400-440 K range because of the high sensitivity of nitric acid to the PSC schemes. Our study shows that, with some limitations (like HN03 comparison under PSC conditions), the THT is a more powerful tool for validation studies than the TCA, making conclusions of the comparison statistically more robust.

Description: Extensive measurement campaigns by the NASA ER-2 research aircraft have obtained a nearly pole-to-pole database of the species that control HOx (OH + HO2) chemistry. The wide dynamic range of these in situ measurements provides an opportunity to demonstrate empirically the mechanisms that control the HOx system. Measurements in the lower stratosphere show a remarkably tight correlation of OH concentration with the solar zenith angle (SZA). This correlation is nearly invariant over latitudes ranging from 70 S to 90 N and all seasons. An analysis of the production and loss of HOx in terms of the rate determining steps of reaction sequences developed by Johnston and Podolske and Johnston and Kinnison is used to clarify the behavior of the system and to directly test our understanding of the system with observations. Calculations using in situ measurements show that the production rate of HOx is proportional to O3 and ultraviolet radiation flux. The loss rate is proportional to the concentration and the partitioning of NOy (reactive nitrogen) and the concentration of HO2. In the absence of heterogeneous reactions, the partitioning of NOy is controlled by O3 and NOx and the concentration of HO2 is controlled by NOy and O3, so that the removal rate of OH is buffered against changes in the correlation of O3 and NOy. The heterogeneous conversion of NO2 to HNO3 is not an important net source of HO, because production and removal sequences are nearly balanced. Changes in NOy partitioning resulting from heterogeneous chemistry have a large effect on the loss rates of HOx, but little or no impact on the measured abundance of OH. The enhanced loss rates at high NO2/HNO3 are offset in the data set examined here by enhanced production rates resulting from increased photolysis rates resulting from the decreased O3 column above the ER-2.

Description: This report provides a broad outline of the total body of research conducted during the past three years. We report on detailed model studies of the precise way in which infrared limb scanning satellites explicitly detect gravity waves, and use these insights to resolve earlier discrepant zonal mean estimates of gravity wave variances from satellite limb scanners. Detailed analysis of CRISTA temperature fluctuations are outlined, which provide new global information on long-wavelength stratospheric gravity waves generated by mountains, tropical convection and the mid-latitude jet stream-vortex system, as well as interactions higher in the mesosphere with the diurnal tide. The detailed insights gained from analyzing CRISTA data are applied to provide a multiyear analysis of stratospheric mountain waves over the Andes evident in MLS limb-track data. We also demonstrate for the first time that stratospheric temperature data from the CLAES instrument on UARS resolved gravity waves, and we highlight mountain wave information in a subset of these data. Final conclusions and recommendations are set forth, and the present findings are directly related to the original goals of the research contract. A full list of publications that resulted from this research is provided.

Description: We model spatial variations in horizontal displacements of 117 geodetic sites measured during annual surveys in 1989-1996 with the Global Positioning System (GPS) as elastic strain across a locked strike-slip fault to infer the contemporary slip rate, locking depth, and location of the Sumatran fault (SF) in northern Sumatra (1 S-3 N). GPS-derived slip rate estimates increase slightly northward from 23 plus or minus 3 mm/yr at 0.8 deg S to 26 plus or minus 2mm/yr at 2.7 N. They agree with geologic estimates north of the Equator, but at 0.5 S they are about 10 mm/yr higher. Strain appears to be distributed asymmetrically about the fault. South of 2 N, about 5 mm/yr of shear is required within the offshore forearc, west of the fault, to achieve a closer agreement of fault locations inferred from GPS velocities with geologically identified traces of the SF. Locking depth estimates are on the order of 10-20 km. The western branch of the major fault bifurcation near 1 N slips at a rate five times higher than the eastern branch. The two main strands of the fault at the northwestern tip of Sumatra (5.5 N) appear to be nearly free of horizontal strain; significant slip must occur away from the two strands, probably further east at two other geologically active branches. The Banda Aceh embayment is extruded to the northwest at a rate of 5 plus or minus 2 mm/yr. Within the estimated velocity uncertainties of several mm/yr, fault-normal deformation along the SF is insignificant. Almost strain free, the northern part of the back-arc basin is part of a rigid Sunda shelf, while the northern forearc is subjected to 8 plus or minus 5 x 10 (exp -8)/yr of extension nearly parallel to the arc.

Description: An analysis combining historical triangulation and recent Global Positioning System (GPS) survey measurements in West and North Sumatra, Indonesia, reveals a detailed slip history along the central part of the Sumatran Fault. The arc-parallel components of the combined velocity field are consistent with slip rates inferred from GPS data, ranging from 23 to 24 mm/yr. Between 1.0 S and 1.3 N the Sumatran Fault appears to be characterized by deep locking depths, on the order of 20 km, and the occurrence of large (M (sub w) approximately 7) earthquakes. The long-term (1883-1993) strains show simple right-lateral shear, with rates similar to GPS-measured, 1989-1993 strain rates. Coseismic deformation due to the 1892 Tapanuli and 1926 Padang Panjang earthquakes, estimated from triangulation measurements taken before and after the events, indicate that the main shocks were significantly larger than previously reported. The 1892 earthquake had a likely magnitude of M (sub w) approximately equal to 7.6, while the 1926 events appear to be comparable in size to the subsequent (M approximately 7) 1943 events, and an order of magnitude higher than previously reported.

Description: The solar magnetic field plays a key role in determining coronal. The principal input to MHD models is the observed solar magnetic field. 3D MHD models can be used to compare with eclipse and coronograph images, SOHO images (LOSCO, EIT), Ulysses and WIND spacecraft data, and interplanetary scintillation (IPS) measurements. MHD computations can tell us about the structure of the corona. Eclipses can help us to verify the accuracy of the models. 4 December, 2002 total eclipce: visible in the southern hemisphere (South Atlantic, southern Africa, Indian Ocean, and Australia). Total in center Angola is at 06:00 UT.

Description: Stress-strain tests were performed both on shock-damaged gabbro and limestone. The effective Young's modulus decreases with increasing initial damage parameter value, and an apparent work-softening process occurs prior to failure. To further characterize shock-induced microcracks, the longitudinal elastic wave velocity behavior of shock-damaged gabbro in the direction of compression up to failure was measured using an acoustic transmission technique under uniaxial loading. A dramatic increase in velocity was observed for the static compressive stress range of 0-50 MPa. Above that stress range, the velocity behavior of lightly damaged (D(sub 0) less than 0.1) gabbro is almost equal to unshocked gabbro. The failure strength of heavily-damaged (D(sub 0) greater than 0.1) gabbro is approx. 100-150 MPa, much lower than that of lightly damaged and unshocked gabbros (approx. 230-260 MPa). Following Nur's theory, the crack shape distribution was analyzed. The shock-induced cracks in gabbro appear to be largely thin penny-shaped cracks with c/a values below 5 x 10(exp -4). Moreover, the applicability of Ashby and Sammis's theory relating failure strength and damage parameter of shock-damaged rocks was examined and was found to yield a good estimate of the relation of shock-induced deficit in elastic modulus with the deficit in compressive strength.

Description: We report paleomagnetic results for core samples of the breccia and andesitic rocks recovered from the Yucatan-6 Petrolcos Mexicanos exploratory well within the Chicxulub structure (about 60 km SSW from its center), northern Yucatan, Mexico. A previous study has shown that the rocks studied contain high iridium levels and shocked breccia clasts and an Ar/Ar date of 65.2 +/- 0.4 Ma. Andesitic rocks are characterized by stable single-component magnetizations with a mean inclination of -42.6 deg +/- 2.4 deg. Breccias present a complex paleomagnetic record characterized by multivectorial magnetizations with widely different initial NRM inclinations. However, after alternating field demagnetization, well defined characteristic components with upward inclinations are defined. IRM acquisition experiments, comparison of IRM and NRM coercivity spectra and the single component magnetization of the andesitic rocks indicate the occurrence of iron-rich titanomagnetites of single or pseudo-single domain states as the dominant magnetic carriers. Mean inclinations from the andesitic rocks and most of the breccia samples give a mean inclination of about -40 deg to -45 deg, indicating a reverse polarity for the characteristic magnetization that is consistent with geomagnetic chron 29R, which spans the Cretaceous/Tertiary (K/T) boundary. The inclination is also consistent with the expected value (and corresponding paleolatitude) for the site estimated from the reference polar wander curve for North America. We suggest that the characteristic magnetizations for the andesitic and breccia rocks are the result of shock heating at the time of formation of the impact structure and that the age, polarity and pateolatitude are consistent with a time at the K/T boundary.

Description: Field aligned currents play an important role in the global coupling between the magnetosphere and the ionosphere and in their relationship to the auroral phenomena. Moreover, there exists evidence that resonant oscillations are related to large-scale Birkeland currents. The spatial confinement of the field-aligned currents forms an inhomogeneous system susceptible to low-frequency oscillations, which can be excited due to periodic variations in the solar wind pressure or to the Kelvin-Helmholtz (KH) instability. In this paper we present a study of ultralow-frequency (ULF) oscillations in an inhomogeneous magnetic field formed by a large-scale current. We investigate the effects of the field-aligned currents on the generation of localized Alfven waves. The field oscillations are described by an eigenvalue wave equation which includes the effects of the field aligned currents and which produces a discrete spectrum of Alfven waves. These waves are observed mainly in three regions of the magnetosphere: in the magnetosheath, in the polar cusp, and in the plasmasphere. In the present study we limit our investigation to the auroral region.

Description: Normally the gravitationally bound heavy O(+) ions in the Earth's ionosphere are in a diffusive equilibrium. However, when energized to superthermal energies of a few eV transverse to the geomagnetic field, the combined effects of the downward gravitational and the upward electric and mirror forces produce interesting flow patterns in the vertical direction like in a pulsating fountain. This flow pattern is studied by means of a particle-in-cell code.

Description: We list a few things that we do not understand about stars and that most people ignore. These are all hard problems. We can learn more cosmology by working on them to reduce the systematic errors they introduce than by trying to derive cosmological results that are highly uncertain.

Description: Studies of Galileo Near Infrared Mapping Spectrometer (NIMS) data and ground based data of volcanism at Prometheus and Loki Patera on Io reveal very different mechanisms of lava emplacement at these two volcanoes. Data analyses show that the periodic nature of Loki Patera s volcanism from 1990 to 2001 is strong evidence that Loki s resurfacing over this period resulted from the foundering of a crust on a lava lake. This process is designated passive , as there is no reliance on sub-surface processes: the foundering of the crust is inevitable. Prometheus, on the other hand, displays an episodicity in its activity which we designate active . Like Kilauea, a close analog, Prometheus s effusive volcanism is dominated by pulses of magma through the nearsurface plumbing system. Each system affords views of lava resurfacing processes through modelling.

Description: Stress accumulation rates along the primary segments of the San Andreas Fault system are computed using a three-dimensional (3-D) elastic half-space model with realistic fault geometry. The model is developed in the Fourier domain by solving for the response of an elastic half-space due to a point vector body force and analytically integrating the force from a locking depth to infinite depth. This approach is then applied to the San Andreas Fault system using published slip rates along 18 major fault strands of the fault zone. GPS-derived horizontal velocity measurements spanning the entire 1700 x 200 km region are then used to solve for apparent locking depth along each primary fault segment. This simple model fits remarkably well (2.43 mm/yr RMS misfit), although some discrepancies occur in the Eastern California Shear Zone. The model also predicts vertical uplift and subsidence rates that are in agreement with independent geologic and geodetic estimates. In addition, shear and normal stresses along the major fault strands are used to compute Coulomb stress accumulation rate. As a result, we find earthquake recurrence intervals along the San Andreas Fault system to be inversely proportional to Coulomb stress accumulation rate, in agreement with typical coseismic stress drops of 1 - 10 MPa. This 3-D deformation model can ultimately be extended to include both time-dependent forcing and viscoelastic response.

Description: The data was analyzed and two papers were written and published in the refereed journal: Meteoritics and Planetary Science. These paper describes the results of the study of the surface chemical composition of the asteroid Eros by the NEAR X-ray Fluorescence Spectrometer.

Description: The abundances of V, Cr, and Mn in the mantles of the Earth and Moon are very similar, and are distinct from all other sampled solar system bodies. This similarity has led to the suggestion that there is a genetic relationship between the Earth and Moon, more specifically, that the Moon was derived from the mantle of the Earth after formation of the core. This proposal is at odds with the Giant Impact paradigm of lunar origin and our current understanding of the formation of the terrestrial planets. During the award period, some of the necessary experimental procedures were developed. Specifically, an assembly design has been developed to permit heating of metal by insulating it from highly conducting diamond. A Nd-YAG laser was assembled, aligned and tested. The technique for preparation of the run products for electron microprobe analysis was developed, and an alternative approach has been identified. A scheme to carry out polybaric experiments on a single charge was developed. The experiments carried out did not lead to conclusive results, and the arguments advanced by Drake et aL (1989) and Drake (1990) remain the best summary of the status of the problem.

Description: The origin of the rough-textured aureoles that surround the immense Olympus Mons volcano on Mars is controversial. We present data from the Mars Global Surveyor and Mars Odyssey missions to demonstrate that at least two of the aureole lobes are derived from the volcano's flanks in large and probably catastrophic mass movement events, leaving behind headwalls that constitute the basal scarp. This evidence stems from the morphology and internal structure of aureole blocks, which exhibit remnants of volcanic flow units on their surfaces. Our claim is supported by plausible reconstructions of the prefailure flanks. Structural analogs to known flank failure events at Hawaiian volcanoes suggest that repeated cycles of flank growth and collapse at Olympus Mons allow generation of the observed aureoles from a protoedifice similar in size and shape to the present one.

Description: We demonstrate that high-resolution 557.7-nm all-sky images are useful tools for investigating the spatial and temporal evolution of merging on the dayside magnetopause. Analysis of ground and satellite measurements leads us to conclude that high-latitude merging events can occur at multiple sites simultaneously and vary asynchronously on time scales of 30 s to 3 min. Variations of 557.7 nm emissions were observed at a 10 s cadence at Ny-Alesund on 19 December 2001, while significant changes in the IMF clock angle were reaching the magnetopause. The optical patterns are consistent with a scenario in which merging occurs around the rim of the high-latitude cusp at positions dictated by the IMF clock angle. Electrons energized at merging sites represent plausible sources for 557.7 nm emissions in the cusp. Polar observations at the magnetopause have directly linked enhanced fluxes of 〉 or = 0.5 keV electrons with merging. Spectra of electrons responsible for some of the emissions, measured during a DMSP F15 overflight, exhibit "inverted-V" features, indicating further acceleration above the ionosphere. SuperDARN spectral width boundaries, characteristic of open-closed field line transitions, are located at the equatorward edge of the 557.7nm emissions. Optical data suggest that with IMF B(sub Y) 〉 0, the Northern Hemisphere cusp divides into three source regions. When the IMF clock angle was approx. 150 deg structured 557.7-nm emissions came from east of the 13:00 MLT meridian. At larger clock angles the emissions appeared between 12:00 and 13:00 MLT. No significant 557.7-nm emissions were detected in the prenoon MLT sector. MHD simulations corroborate our scenario, showing that with the observed large dipole-tilt and IMF clock angles, merging sites develop near the front and eastern portions of the high-altitude cusp rim in the Northern Hemisphere and near the western part of the cusp rim in the Southern Hemisphere.

Description: Magnetic merging on the dayside magnetopause often occurs at high latitudes. Polar measured fluxes of accelerated ions and wave Poynting vectors while skimming the subsolar magnetopause. The measurements indicate that their source was located to the north of the spacecraft, well removed from expected component merging sites. This represents the first use of wave Poynting flux as a merging discriminator at the magnetopause. We argue that wave Poynting vectors, like accelerated particle fluxes and the Walen tests, are necessary, but not sufficient, conditions, for identifying merging events. The Polar data are complemented with nearly simultaneous measurements from Cluster in the northern cusp, with correlated observations from the SuperDARN radar, to show that the locations and rates of merging vary. Magnetohydrodynamic (MHD) simulations are used to place the measurements into a global context. The MHD simulations confirm the existence of a high-latitude merging site and suggest that Polar and SuperDARN observed effects are attributable to both exhaust regions of a temporally varying X-line. A survey of 13 merging events places the location at high latitudes whenever the interplanetary magnetic field (IMF) clock angle is less than approximately 150 degrees. While inferred high-latitude merging sites favor the antiparallel merging hypothesis, our data alone cannot exclude the possible existence of a guide field. Merging can even move away from equatorial latitudes when the IMF has a strong southward component. MHD simulations suggest that this happens when the dipole tilt angle increases or when IMF B(sub X) increases the effective dipole tilt.

Description: A series of VLBA experiments were carried out at K and Q bands for astrometry and imaging within the KQ VLBI Survey Collaboration. The paired K and Q observations of each source are separated by approximately 3 minutes of time. We investigate the delay effect of the ionosphere between K and Q bands involving the interscan separation. This differential delay effect is intermixed with the differential fluctuation effect of the troposphere.

Description: Radio sounding of the Earth's topside ionosphere and magnetosphere is a proven technique from geospace missions such as the International Satellites for Ionospheric Studies (ISIS) and the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE). Application of this technique to Jupiter's icy moons and the surrounding Jovian magnetosphere will provide unique remote sensing observations of the plasma and magnetic field environments and the subsurface conductivities, of Europa, Ganymede, and Callisto. Spatial structures of ionospheric plasma above the surfaces of the moons vary in response to magnetic-field perturbations from (1) magnetospheric plasma flows, (2) ionospheric currents from ionization of sputtered surface material, and (3) induced electric currents in salty subsurface oceans and from the plasma flows and ionospheric currents themselves. Radio sounding from 3 kHz to 10 MHz can provide the global electron densities necessary for the extraction of the oceanic current signals and supplements in-situ plasma and magnetic field measurements. While radio sounding requires high transmitter power for subsurface sounding, little power is needed to probe the electron density and magnetic field intensity near the spacecraft. For subsurface sounding, reflections occur at changes in the dielectric index, e.g., at the interfaces between two different phases of water or between water and soil. Variations in sub-surface conductivity of the icy moons can be investigated by radio sounding in the frequency range from 10 MHz to 50 MHz, allowing the determination of the presence of density and solid-liquid phase boundaries associated with oceans and related structures in overlying ice crusts. The detection of subsurface oceans underneath the icy crusts of the Jovian moons is one of the primary objectives of the Jupiter Icy Moons Orbiter (JIMO) mission. Preliminary modeling results show that return signals are clearly distinguishable be&een an ice crust with a thickness of 7 km on 1) an ocean and 2) a layer of bedrock. Knowledge of the ionospheric contributions to the time delay of the low-frequency subsurface radar is shown to be important in obtaining accurate depth information.

Description: The origin of whistler mode radiation in the plasmasphere is examined from three years of plasma wave observations from the Dynamics Explorer and three years from the Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) spacecraft. These data are used to construct plasma wave intensity maps of whistler mode radiation in the plasmasphere. The highest average intensities of the radiation in the wave maps show source locations and/or sites of wave amplification. Each type of emission is classified based on its magnetic latitude and longitude rather than any spectral feature. Equatorial electromagnetic (EM) emissions (approx. 30-330 Hz), plasmaspheric hiss (approx. 330 Hz - 3.3 kHz), chorus (approx. 2 kHz - 6 kHz), and VLF transmitters (approx. 10-50 kHz) are the main types of waves that are clearly delineated in the plasma wave maps. Observations of the equatorial EM emissions show that the most intense region is on or near the magnetic equator in the afternoon sector and that during times of negative B(sub z) (interplanetary magnetic field),the maximum intensity moves from L values of 3 to less than 2. These observations are consistent with the origin of this emission being particle-wave interactions in or near the magnetic equator. Plasmaspheric hiss shows high intensity at high latitudes and low altitudes (L shells from 2 to 4) and in the magnetic equator over L values from 2 to 3 in the early afternoon sector. The longitudinal distribution of the hiss intensity (excluding the enhancement at the equator) is similar to the distribution of lightning: stronger over continents than over the ocean, stronger in the summer than winter, and stronger on the dayside than nightside. These observations strongly support lightning as the dominant source for plasmaspheric hiss, which through particle-wave interactions, maintains the slot region in the radiation belts. The enhancement of hiss at the magnetic equator is consistent with particle-wave interactions. The chorus emissions are most intense on the morning side as previously reported. At frequencies from approx. 10-50 kHz VLF transmitters dominate the spectrum. The maximum intensity of the VLF transmitters is in the late evening or early morning with enhancements all along L shells from 1.8 to 3.

Description: Results of an experimental and theoretical investigation into the effects of the time dependent neutral wind flywheel on high-latitude ionospheric electrodynamics are presented. The results extend our previous work which used the National Center for Atmospheric Research Thermosphere/Ionosphere General Circulation Model (NCAR TIGCM) to theoretically simulate flywheel effects in the aftermath of a geomagnetic storm. The previous results indicated that the neutral circulation, set up by ion-neutral momentum coupling in the main phase of a geomagnetic storm, is maintained for several hours after the main phase has ended and may dominate height-integrated Hall currents and field-aligned currents for up to 4-5 hours. We extend the work of Deng et al. to include comparisons between the calculated time-dependent ionospheric Hall current system in the storm-time recovery period and that measured by instruments on board the Dynamics Explorer 2 (DE 2) satellite. Also, comparisons are made between calculated field-aligned currents and those derived from DE 2 magnetometer measurements. These calculations also allow us to calculate the power transfer rate (sometimes called the Poynting flux) between the magnetosphere and ionosphere. The following conclusions have been drawn: (1) Neutral winds can contribute significantly to the horizontal ionospheric current system in the period immediately following the main phase of a geomagnetic storm, especially over the magnetic polar cap and in regions of ion drift shear. (2) Neutral winds drive Hall currents that flow in the opposite direction to those driven by ion drifts. (3) The overall morphology of the calculated field-aligned current system agrees with previously published observations for the interplanetary magnetic field (IMF) B(sub Z) southward conditions, although the region I and region 2 currents are smeared by the TI(ICM model grid resolution. (4) Neutral winds can make significant contributions to the field-aligned current system when B(sub Z) northward conditions prevail following the main phase of a storm, but can account for only a fraction of the observed currents. (5) DE 2 measurements provide a demonstration of "local" (satellite-altitude) flywheel effects. (6) On the assumption that the magnetosphere acts as an insulator, we calculate neutral-wind-induced polarization electric fields of approx. 20-30 kV in the period immediately following the geomagnetic storm.

Description: In situ measurements of hydrogen, nitrogen, and chlorine radicals obtained in the lower stratosphere during the Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE) are compared to results from a photochemical model that assimilates measurements of radical precursors and environmental conditions. Models allowing for heterogeneous hydrolysis of N2O5 agree well with measured concentrations of NO and ClO, but concentrations of HO2 and OH are underestimated by 10 to 25%, concentrations of NO2 are overestimated by 10 to 30%, and concentrations of HCl are overestimated by a factor of 2. Discrepancies for [OH] and [HO2] are reduced if we allow for higher yields of O(sup 1)D) from 03 photolysis and for heterogeneous production of HNO2. The data suggest more efficient catalytic removal of O3 by hydrogen and halogen radicals relative to nitrogen oxide radicals than predicted by models using recommended rates and cross sections. Increases in [O3] in the lower stratosphere may be larger in response to inputs of NO(sub y) from supersonic aircraft than estimated by current assessment models.

Description: Compositions and textures of melt rocks from the upper part of the Chicxulub structure are typical of melt rocks at other large terrestrial impact structures. Apart from variably elevated iridium concentrations (less than 1.5 to 13.5 +/- 0.9 ppb) indicating nonuniform dissemination of a meteoritic component, bulk rock and phenocryst compositions imply that these melt rocks were derived exclusively from continental crust and platform-sediment target lithologies. Modest differences in bulk chemistry among samples from wells located approximately 40 km apart suggest minor variations in relative contributions of these target lithologies to the melts. Subtle variations in the compositions of early-formed pyroxene and plagioclase also support minor primary differences in chemistry between the melts. Evidence for pervasive hydrothermal alteration of the porous mesostasis includes albite, K-feldspar, quartz, epidote, chlorite, and other phyllosilicates, as well as siderophile element-enriched sulfides, suggesting the possibility that Chicxulub, like Sudbury, may host important ore deposits.

Description: Recent satellite observations at high altitudes (greater than 5000 km) in the auroral zone have shown the existence of hybrid or bimodal ion beam distributions that are evidence of both parallel and perpendicular ion acceleration. Acceleration parallel to the magnetic field is most likely due to quasi-static electric fields (double layers) which can create outflowing ion beams; since ions of different mass will have different drift speeds due to this acceleration, a plasma configuration unstable to the ion-ion two-stream acoustic mode develops. When the net drift velocity (U) between the two ion species is greater than the sound speed (C(sub 0)), the ion-ion instability has maximum growth at oblique wave propagation. To study the nonlinear effects of the ion-ion instability in terms of plasma heating, a numerical simulation parametric study has been performed. It was found that the parallel acceleration that forms the ion beams occurs on a time scale faster than ion- ion wave growth at low drifts; thus ion-ion wave growth is expected to occur primarily for higher drift speeds (U greater than C(sub 0)) which results in strong oblique heating of the ions (both hydrogen and oxygen) forming elevated ion conics (sometimes called 'bowl' distributions). Also, strong parallel electron heating in the direction of the ion beams can occur, and electrons near the top of the acceleration region may attain a net upward drift along with the elevated ion conics. Variation of the oxygen density greatly affects the ion heating due to the ion-ion instability; as the oxygen density decreases, oxygen heating increases, in agreement with observations (Collin et al., 1987). Ion-ion electrostatic wave properties and the plasma heating that results over a wide range of auroral zone parameters are included.

Description: The magnetospheric convection electric field contributes to Birkeland currents. The effects of the field are to polarize the plasma by displacing the bounce paths of the ions from those of electrons, to redistribute the pressure so that it is not constant along magnetic field lines, and to enhance the pressure gradient by the gradient of the bulk speed. Changes in the polarization charge during the convection of the plasma are neutralized by electrons in the form of field-aligned currents that close through the ionosphere. The pressure drives field-aligned currents through its gradient in the same manner as in quasi-static plasma, but with modifications that are important if the bulk speed is of the order of the ion thermal speed; the variations in the pressure along field lines are maintained by a weak parallel potential drop. These effects are described in terms of the field-aligned currents in steady state, isotropic, MED plasma. Solutions are developed by taking the MHD limit of two-fluid solutions and illustrated in the special case of Maxwellian plasma for which the temperature is constant along magnetic field lines. The expression for the Birkeland current density is a generalization of Vasyliunas' expression for the field-aligned current density in quasi-static plasma and provides a unifying expression when both pressure gradients and ion inertia operate simultaneously as sources of field-aligned currents. It contains a full account of different aspects of the ion flow (parallel and perpendicular velocity and vorticity) that contribute to the currents. Contributions of ion inertia to field-aligned currents will occur in regions of strong velocity shear, electric field reversal, or large gradients in the parallel velocity or number density, and may be important in the low-latitude boundary layer, plasma sheet boundary layer, and the inner edge region of the plasma sheet.

Description: Direct imaging of the magnetosphere by instruments on the IMAGE spacecraft is supplemented by simultaneous observations of the global aurora in three far ultraviolet (FUV) wavelength bands. The purpose of the multi-wavelength imaging is to study the global auroral particle and energy input from thc magnetosphere into the atmosphere. This paper describes provides the method for quantitative interpretation of FUV measurements. The Wide-Band Imaging Camera (WIC) provides broad band ultraviolet images of the aurora with maximum spatial and temporal resolution by imaging the nitrogen lines and bands between 140 and 180 nm wavelength. The Spectrographic Imager (SI), a dual wavelength monochromatic instrument, images both Doppler-shifted Lyman alpha emissions produced by precipitating protons, in the SI-12 channel and OI 135.6 nm emissions in the SI-13 channel. From the SI-12 Doppler shifted Lyman alpha images it is possible to obtain the precipitating proton flux provided assumptions are made regarding the mean energy of the protons. Knowledge of the proton (flux and energy) component allows the calculation of the contribution produced by protons in the WIC and SI-13 instruments. Comparison of the corrected WIC and SI-13 signals provides a measure of the electron mean energy, which can then be used to determine the electron energy fluxun-. To accomplish this reliable modeling emission modeling and instrument calibrations are required. In-flight calibration using early-type stars was used to validate the pre-flight laboratory calibrations and determine long-term trends in sensitivity. In general, very reasonable agreement is found between in-situ measurements and remote quantitative determinations.

Description: The EOS Aura Mission is designed to make comprehensive chemical measurements of the troposphere and stratosphere. In addition the mission will make measurements of important climate variables such as aerosols, and upper tropospheric water vapor and ozone. Aura will launch in late 2003 and will fly 15 minutes behind EOS Aqua in a polar sun synchronous ascending node orbit with a 1:30 pm equator crossing time.

Description: An embodiment of the present invention uses a single detection system to approximate a location of lightning strikes. This system is triggered by a broadband RF detector and measures a time until the arrival of a leading edge of the thunder acoustic pulse. This time difference is used to determine a slant range R from the detector to the closest approach of the lightning. The azimuth and elevation are determined by an array of acoustic sensors. The leading edge of the thunder waveform is cross-correlated between the various acoustic sensors in the array to determine the difference in time of arrival, AT. A set of AT S is used to determine the direction of arrival, AZ and EL. The three estimated variables (R, AZ, EL) are used to locate a probable point of the lightning strike.

Description: In string theories, extra dimensions must be compactified. The possibility that gravity can have large radii of compactification leads to a violation of the inverse square law at submillimeter distances. The objective of ISLES is to perform a null test of Newton s law in space with a resolution of one part in 10(exp 5) or better at 100 microns. The experiment will be cooled to less than or equal to 2 K, which permits superconducting magnetic levitation of the test masses. To minimize Newtonian errors, ISLES employs a near null source, a circular disk of large diameter-to-thickness ratio. Two test masses, also disk-shaped, are suspended on the two sides of the source mass at a nominal distance of 100 microns. The signal is detected by a superconducting differential accelerometer. A ground test apparatus is under construction.

Description: Accurate analysis of precision ranges to the Moon have provided several tests of gravitational theory: the equivalence principle, geodetic precession, PPN parameters beta and gamma, and the constancy of the gravitational constant G. Other possible tests include the inverse square law at 20,000 km length scales and the PPN parameter 1. The uncertainties of these tests have decreased as data accuracies have improved and data time span has lengthened. We are exploring the modeling improvements necessary to proceed from cm to mm range accuracies. Looking to future exploration, what characteristics are desired for the next generation of ranging devices, what fundamental questions can be investigated, and what are the challenges for modeling and data analysis?

Description: A 2710 g meteorite, Sahara 99555 (Sah99), was recently recovered from the Sahara and reported to be the 5th angrite. It is the largest angrite ever found and may offer useful information to better understand the unusual petrogeneses of this rare achondrite group. It may also allow us to examine the chronological record of igneous activity in the very early solar system. We obtained a 2.6 g chip of Sah99 and here present a preliminary report of its petrology and mineralogy in conjunction with a crystallization experiment on an analogue composition.

Description: In situ observations indicate that the dynamical processes in the geoplasma environment generally entail localized intermittent processes and anomalous global transports. It was suggested by T. Chang that instead of considering the turbulence as a mixture of interacting waves, such type of patchy intermittency could be more easily understood in terms of the development, interaction, merging, preferential acceleration and evolution of coherent magnetic structures. In this three-year project, we have used direct numerical MHD simulations to study some aspects of the MHD dynamics in Chang's model. Our large-scale numerical calculations and simulations have been supplemented by and coordinated with theoretical studies conducted by Chang and his colleagues.

Description: This volume is the proceedings of the third General Meeting of the International VLBI Service for Geodesy and Astromctry IVS), held in Otlawa, Canada, February 9-11,2004. The keynote of the third GM was visions for the next decade following the main theme of "Today's Results and Tomorrow's Vision". with a recognition that the outstanding VLBI results available today are the foundation and motivation for the next generation VLBI system requirements. The goal of the meeting was to provide an interesting and informative program for a wide cross section of IVS members, including station operators, program managers, and analysts.

Description: On January 12, 2003 the Ice, Cloud and land Elevation Satellite (ICESat) was successfUlly placed into orbit. The ICESat mission carries the Geoscience Laser Altimeter System (GLAS), which consists of three near-infrared lasers that operate at 40 short pulses per second. The instrument has collected precise elevation measurements of the ice sheets, sea ice roughness and thickness, ocean and land surface elevations and surface reflectivity. The accurate geolocation of GLAS's surface returns, the spots from which the laser energy reflects on the Earth's surface, is a critical issue in the scientific application of these data Pointing, ranging, timing and orbit errors must be compensated to accurately geolocate the laser altimeter surface returns. Towards this end, the laser range observations can be fully exploited in an integrated residual analysis to accurately calibrate these geolocation/instrument parameters. Early mission ICESat data have been simultaneously processed as direct altimetry from ocean sweeps along with dynamic crossovers resulting in a preliminary calibration of laser pointing, ranging and timing. The calibration methodology and early mission analysis results are summarized in this paper along with future calibration activities

Description: This CD-ROM contains the abstracts and associated files for the Eleventh Annual V.M Goldschmidt Conference. Topics include: Organic Geochemistry; Metamorphic Processes; Igneous Processes; Stable and Radiogenic Isotopes; Planetary Geochemistry and Mineralogy; Mineralogy and Crystallography; Ore Deposits; and Aqueous Geochemistry.

Description: The IGS network today consists of 364 stations managed by about, 100 different agencies worldwide. The RINEX observations are contributed to the IGS Data Centers, which permanently archive the data and make it freely available to all users. The primary customer of the data set is the IGS Analysis Centers, which acquire the data for generation of precise GPS products such as ephemerides, clocks, earth orientation parameters, and station position and velocities.

Description: The Hydrosphere State (HYDROS) Mission has been selected for the National Aeronautics and Space Administration (NASA) Earth System Science Pathfinder (ESSP) program. The objectives of HYDROS are to provide frequent, global measurements of surface soil moisture and surface freeze/thaw state. In order to adequately measure these geophysical parameters, a system capable of simultaneously measuring L-Band radiometer brightness temperatures at 40 km resolution and L-Band radar backscatter at 3 km resolution over a very wide swath is required. In addition, these science requirements must be satisfied under the stringent cost-cap imposed on all ESSP missions. As a solution to this challenging set of requirements, a relatively large, six meter, conically-scanning reflector antenna architecture was selected for the mission design. The HYDROS instrument will fly on a General Dynamics SA-200HP spacecraft bus. Although large deployable mesh antennas have been used in communication applications, this will mark the first time such technology is applied in a rotating configuration for high-resolution remote sensing.

Description: The science objectives of the HYDROS mission are to provide frequent, global measurements of surface soil moisture and surface freeze/thaw state. In order to adequately measure these geophysical quantities, the key instrument requirements were determined by the HYDROS science team to be: 1) Dual-polarization L-Band passive radiometer measurements at 40 km resolution, 2) Dual-polarization L-Band active radar measurements at 3 km resolution, and 3) A wide swath to insure global three day refresh time for these measurements (1000 km swath at the selected orbit altitude of 670 km). As a solution to this challenging set of instrument requirements, a relatively large, 6 meter, conically-scanning reflector antenna architecture was selected for the instrument design. The deployable mesh antenna is shared by both the radiometer and radar electronics by employing a single L-Band feed.

Description: Salt marshes are constant depositional environments and as a result contain accurate indicators of past relative sea level rise and salinity. The Hudson River marshes are at least twice as deep when compared to coastal marshes on either side of the mouth of the Hudson. The reason for this difference in sedimentation is unclear. This study uses macrofossil data as well as sediment stratigraphy in order to understand the formation and evolution of these marshes. The composition of seeds, roots, shoots and foraminifera, are used to indicate past sea levels. The four sites involved in this study are, from south to north, the Arthur Kill Marsh in Staten Island (40 36 N, 74 77W), Piermont marsh (N 4100; 73 55W) Croton Point (41 14 N; 73 50W) and Iona Island (41 18N, 73 58W). These are all tidally influenced but with increasing distances from the New York Bight, which gives a good spectrum of tidal influence. AMS-C14 dates on basal macrofossils will document the time of each marsh formation. Basal material from Arthur Kill (8 m) includes freshwater seeds such as Viola, Potomageton and Alnus along with Salix buds. Basal material from Croton Point (10 m) includes fibrous woody material, foraminifera and Zanichellia seeds and other brackish vegetational components. The basal material from Piermont (13.77 m) is lacking any identifiable macrofossils between 150 and 500 microns. The basal material from Iona Island (10 m) has vegetation such as Scirpus and Cyperus seeds, probably implying a brackish environment. The freshwater origin of the Arthur Kill marsh in Staten Island is significant because it predates either sea level rise or the western channel incision. Additional implications for this study include evidence for changes in river channel geomorphology. Reasons for the relatively deeper river marshes include possible basal clay compaction, high production due to river and marine nutrients as well as tectonic activity. This study provides the groundwork for more high-resolution studies of these marshes to understand the fluctuations in salinity caused by relative sea level rise, tectonic faulting and/or changes in precipitation/evaporation.

Description: The principal objective of this grant and this research were to investigate the topographic development of an active glaciated orogenic belt in southern Alaska as that development relates to patterns of erosion and crustal deformation. A specific objective of the research was to investigate feedbacks between mountain building, orographic affects on climate, and patterns of exhumation and rock uplift. To that end, an orogen-scale analysis of topography was conducted with the aid of digital elevation models, magnitudes and patterns of crustal deformation were compiled from existing literature, present and past climate patterns were constrained using the modern and past distribution of glaciers, and styles, magnitudes, and extent of erosion were constrained with observations from the 1998 field season.

Description: Phosphorus incorporation in major rock-forming silicate minerals has the following implications: (1) Reactions between phosphorus-hosting major silicates and accessory phosphates, which are also major trace element carriers, may control the stability of the latter and thus may affect the amount of phosphorus and other trace elements released to the coexisting melt or fluid phase. (2) Less of a phosphate mineral is needed to account for the bulk phosphorus of planetaty mantles. (3) During partial melting of mantle mineral assemblages or equilibrium fractional crystallization of basaltic magmas, and in the absence or prior to saturation with a phosphate mineral, silicate melts may become enriched in phosphorus, especially in the geochemically important low melt fraction regime, Although the small differences in the ionic radii of IVp5+, IVSi4+, and IV Al3+ makes phosphoms incorporation into crystalline silicates perhaps unsurprising, isostructural silicate and phosphate crystalline solids do not readily form solutions, e.g., (Fe, Mg)2SiO4 vs. LiMgPO4, SiO)2 VS. AlPO4. Nonetheless, there are reports of, poorly characterized silico-phosphate phases in angrites , 2-4 wt% P2O5 in olivine and pyroxene grains in pallasites and reduced terestrial basalts which are little understood but potentially useful, and up to 17 wt% P2O5 in olivine from ancient slags. However, such enrichments are rare and only underscore the likelihood of phosphoms incorporation in silicate minerals. The mechanisms that allow phosphorus to enter major rock-forming silicate minerals (e.g., Oliv, Px, Gt) remain little understood and the relevant data base is limited. Nonetheless, old and new high-pressure (5-10 GPa) experimental data suggest that P2O5 wt% decreases from silica-poor to silica-rich compositions or from orthosilicate to chain silicate structures (garnet 〉 olivine 〉 orthopyroxene) which implies that phosphorus incorporation in silicates is perhaps more structure-than site-specific. The data also indicate that DXVliQP2O5 decrease in the same order, but DOVLiQP2O5 and DOpx/LiQP205 are likely constant, respectively equal to 0.08(3) and 0.007(4), in contrast, DG1ILiQP205 increases from 0.15(3) to 0.36(10) as garnet becomes majoritic, thus silica-enriched, and may also depend on liquid composition (SiO2, P2O5 and Na2O wt%).

Description: Basaltic and anorthositic glasses were subjected to aqueous weathering conditions in the laboratory where the variables were pH, temperature, glass composition, solution composition, and time. Leached layers formed at the surfaces of glasses followed by the precipitation of X-ray amorphous iron and titanium oxides in acidic and neutral solutions at 25 C over time. Glass under oxidative hydrothermal treatments at 150 C yielded a three-layered surface; which included an outer smectite layer, a Fe-Ti oxide layer and an innermost thin leached layer. The introduction of Mg into solutions facilitated the formation of phyllosilicates. Aqueous hydrothermal treatment of anorthositic glasses (high Ca, low Ti) at 200 C readily formed smectite, whereas, the basaltic glasses (high Ti) were more resistant to alteration and smectite was not observed. Alkaline hydrothermal treatment at 2000e produced zeolites and smectites; only smectites formed at 200 C in neutral solutions. These mineralogical changes, although observed under controlled conditions, have direct applications in interpreting planetary (e.g., meteorite parent bodies) and terrestrial aqueous alteration processes.

Description: The Athena Science Instrument Payload is providing geochemical and mineralogical information for determining the properties of rocks, soils, and outcrops at the Mars Exploration Rovers landing sites. These measurements indicate that a variety of aqueous processes as well as various degrees of alteration occurred at the two landing sites. Light-toned rocks around the Spirit landing site appear to have coatings or alteration rinds that may have resulted from limited aqueous alteration on the surfaces of basaltic rocks. Hematite and high Fe(III)lFe(total) occur at the surfaces of these rocks. High concentrations of elements highly mobile in water (i.e., S, Cl, and Br) occur in rock veins, vugs, and coatings and at the bottom of soil trenches in the "intercrater plains." One scenario for the formation of rock coatings or rinds and translocation of mobile elements is that water might have occurred briefly at the Martian surface during periods of high obliquity and thin films of water may have mobilized elements and altered the surfaces of rocks. Outcrops on the slopes of the Columbia Hills appear to be extensively altered as suggested by their relative "softness" (measured as resistance to abrasion) as compared to basalts on the adjacent plains, high Fe(III)lFe(total), iron mineralogy dominated by nanophase Fe(III) oxides and hematite, and high Br and CI concentrations beneath outcrop surfaces. These outcrops may have formed by the alteration of basaltic rocks and/or volcaniclastic materials by solutions that were rich in volatile elements (e.g., Br, CI, S). However, it is not clear whether aqueous alteration occurred at depth (e.g., metasomatism), by hydrothermal solutions (e.g., associated with volcanic or impact processes), by vapors rich in volcanic gases, or by low-temperature solutions. The occurrence of jarosite, hematite, and other sulfates (e.g., Mg sulfates) in Eagle and Endurance crater outcrops are strong indicators of aqueous processes at Meridiani Planum. These phases occur with siliciclastic materials in outcrops. Jarosite can only form by aqueous processes under very acidic conditions; e.g., acid sulfate weathering conditions resulting from the oxidation of Fe sulfides or by sulfuric acid alteration of basalts by solutions associated with S02-rich volcanic gases. It is plausible that acidic solutions rich in sulfur (and Fe(II)) reacted with basaltic sediments (which provided a host of soluble cations) under oxidizing conditions and then, through evaporation, formed sediments rich in jarosite and other sulfates along with siliciclastic materials. Hematite-rich spherules in outcrops may have formed by aqueous processes within the sedimentary layers, which promoted transport of Fe (II) solutions to nucleation sites where oxidation and precipitation occurred to form hematite-rich spherules.

Description: This paper presents the prototype of a predictive model capable of describing both magnitudes of deforestation and its spatial articulation into patterns of forest fragmentation. In a departure from other landscape models, it establishes an explicit behavioral foundation for algorithm development, predicated on notions of the peasant economy and on household production theory. It takes a 'bottom-up' approach, generating the process of land-cover change occurring at lot level together with the geography of a transportation system to describe regional landscape change. In other words, it translates the decentralized decisions of individual households into a collective, spatial impact. In so doing, the model unites the richness of survey research on farm households with the analytical rigor of spatial analysis enabled by geographic information systems (GIs). The paper describes earlier efforts at spatial modeling, provides a critique of the so-called spatially explicit model, and elaborates a behavioral foundation by considering farm practices of colonists in the Amazon basin. It then uses, insight from the behavioral statement to motivate a GIs-based model architecture. The model is implemented for a long-standing colonization frontier in the eastern sector of the basin, along the Trans-Amazon Highway in the State of Para, Brazil. Results are subjected to both sensitivity analysis and error assessment, and suggestions are made about how the model could be improved.

Description: Any large mass transport in the Earth system produces changes in the gravity field. Via the space geodetic technique of satellite-laser ranging in the last quarter century, the Earth's dynamic oblateness J2 (the lowest-degree harmonic component of the gravity field) has been observed to undergo a slight decrease -- until around 1998, when it switched quite suddenly to an increase trend which has continued to 2001 before sharply turning back to the value which it is "supposed to be"!. The secular decrease in J2 has long been attributed primarily to the post-glacial rebound in the mantle; the present increase signifies an even larger change in global mass distribution whose J2 effect overshadows that of the post-glacial rebound, at least over interannual timescales. Intriguing evidences have been found in the ocean water distribution, especially in the extratropical Pacific basins, that may be responsible for this J2 change. New techniques based on satellite-to-satellite tracking will yield greatly improved observations for time-variable gravity, with much higher precision and spatial resolution (i.e., much higher harmonic degrees). The most important example is the GRACE mission launched in March 2002, following the success of the CHAMP mission. Such observations are becoming a new and powerful tool for remote sensing of geophysical fluid processes that involve larger-scale mass transports.

Description: Orthorhombic, (Bbmm), (Al, Fe, Cr, Ti)(sub 2) TiO5-(Mg, Fe)Ti2O5 solid solutions (pseudobrookites, s.l.) are found either as an oxidation product of ilmenite and/or spinel or a primary crystallizing phase in igneous and metamorphic rocks on Earth (e.g., basalt flows, crustal and mantle xenoliths, hornfels), and basaltic rocks on the Moon. Moreover, orthorhombic oxides are often part of the crystalline matrix in glass/ceramics with useful applications, and play a major role in the industrial production of TiO2. To fully exploit the potential of these compounds as petrogenetic indicators and/or useful materials we need to quantitatively understand the factors controlling their properties and stability, and thus, to extrapolate beyond the calibrating experiments. For that purpose, we need to combine thermochemistry, phase equilibrium, and in situ P-V-T-cation disorder experimental data that presently either are incomplete or lacking. Perhaps, the most complete data set is that for MgTi2O5 (karrooite) which allows the calibration of models for the Gibbs free energy of the MgTi2O5 as a function of pressure, temperature, and the Mg2+-Ti4+ distribution between the two nonequivalent octahedral sites. Consequently, the effect of cation disorder on MgTi2O5 stability, and the phase relations among MgTi2O5, other titanium oxides, and silicate minerals can be examined. Calculated phase relations in the Mg-Ti-Si-O system and phase equilibrium experiments in Fe-bearing compositions suggest that pseudobrookite-type oxides may be a more common in rocks than previously realized. However, homogeneous and heterogeneous equilibria, and crystallization paths likely affect their stability. For example, isobaric increases in temperature favor disordering and thus entropy-stabilization, in contrast, isothermal increases in pressure have the opposite effect. Although, currently, the potential effect of composition to cation disorder cannot be fully explored, it appears that enrichment in trivalent cations probably enhances entropy-stabilization and thus may increase the stability of (Al, Fe, Cr, Ti)-rich pseudobrookites relative to that of (Mg, Fe)-rich ones. In addition, high-temperature, nearly isothermal, decompression paths of olivine+orthopyroxene+oxide assemblages may favor pseudobrookites (s.l.) over rutile and/or ilmenite, in contrast, cooling at low pressures seems to favor ilmenite and/or rutile. In the case of crustal and mantle xenoliths, the presence or absence of orthorhombic oxides is probably controlled by reactions with olivine, orthopyroxene, ilmenite, and rutile. In oceanic mantle xenoliths such reactions may also involve a TiO2-enriched but not SiO2-enriched melt/fluid, because pseudobrookites (s.l.) would react with the SiO2-enriched melt/fluid to form orthopyroxene and rutile. Parenthetically, experiments and model calculations in the Mg-Ti-Si-O system suggest that low degree partial melting of low-TiO2 bulk compositions may produce Ti-enriched liquids in equilibrium with olivine, orthopyroxen ad=nd MgTi2O5, rutile or ilmenite.

Description: Any large mass transport in the Earth system produces changes in the gravity field. Low harmonic degree components of such variations have been observed by the satellite-laser-ranging (SLR) technique, particularly in 52, the Earth's dynamic oblateness. 52 has long been observed to undergo a slight decrease due to the post-glacial rebound of the mantle -- until around 1998, when it switched quite suddenly to an increase trend which continued to 2001 before sharply turning back to normal , signifying a large change in global mass distribution whose 52 effect overshadows that of the post-glacial rebound over interannual timescales. Intriguing evidences have been found in the Ocean water distribution, especially related to the Pacific Decadal Oscillation in the extratropical Pacific basins, that may be responsible for this 52 anomaly. Besides the lowest-degree 52, recent updates in the SLR-derived time series of the Earth's low-degree gravity components also show shorter wavelength zonal and other longitudinal signals. While the formal uncertainty of these terms is significantly higher than that for 52, some of these series have significant signal that appears to be climatic in origin. For example, there is a significant correlation of the sectoral S2,2 with the Tahiti-Darwin Southern Oscillation Index (SOI), but preceding SO1 by about 1 year. Similar cases demonstrate the utility of assessing the mass component of climate variations. We will present these intriguing results.

Description: Conjugate studies of high-latitude geomagnetic activity are central to understanding the (truly) global magnetospheric response to external perturbations, i.e. the solar wind, and the role of the ionosphere in modulating and coupling with the magnetosphere. Interhemispheric asymmetries as manifested in auroral emissions have been observed for over 40 years. Unfortunately, the nature of the problem has limited the type and extent of studies that can be performed to either conjugate ground based observations or comparison of space based images with all-sky cameras. Recently Frank and Sigwarth (2002) published results using unique simultaneous images from Polar of the northern and southern aurora, but these are necessarily limited to the nighttime and at oblique angles. However, conjunctive images made with the cameras from the Polar and IMAGE missions contain simultaneous conjugate images of the large-scale aurora under many and varying conditions. This rich data set provides an opportunity to study interhemispheric auroral asymmetries and investigate their occurrence as a function of solar wind conditions and ionospheric parameters such as conductivity. In order to use images of different scenes from each of these cameras, knowledge of their relative response is required. To that end, this paper will present preliminary results of comparing the images with emphasis on data from the Polar UVI LBH filters and the IMAGE FUV WIC.

Description: The dominant role of the latitudinal peak of the sea surface temperature (SST) in determining the latitudinal location of the intertropical convergence zone (ITCZ) is well-known. However, the roles of the other factors are less well-known and are the topic of this study. These other factors include the inertial stability, the interaction between convection and surface fluxes and the interaction between convection and radiation. Since these interactions involve convection, in a model they involve the cumulus parameterization scheme. These factors are studied with a general circulation model with uniform SST and solar angle. Under the aforementioned model settings, the latitudinal location of the ITCZ is the latitude where the balance of two types of attraction on the ITCZ, both due to earth's rotation, exists. Directly related to the Coriolis parameter, the first type pulls the ITCZ toward the equator and is not sensitive to model design changes. Related to the convective circulation, the second type pulls the ITCZ poleward and is sensitive to model design changes. Due to the shape and the magnitude of the attractors, the balance of the two types of attractions is reached either at the equator or more than 10 degrees away from the equator. The former case results in a single ITCZ over the equator and the latter case a double ITCZ straddling the equator.

Description: One of the most promising methods to test the representation of cloud processes used in climate models is to use observations together with Cloud Resolving Models (CRMs). The CRMs use more sophisticated and realistic representations of cloud microphysical processes, and they can reasonably well resolve the time evolution, structure, and life cycles of clouds and cloud systems (size about 2-200 km). The CRMs also allow explicit interaction between out-going longwave (cooling) and incoming solar (heating) radiation with clouds. Observations can provide the initial conditions and validation for CRM results. The Goddard Cumulus Ensemble (GCE) Model, a cloud-resolving model, has been developed and improved at NASA/Goddard Space Flight Center over the past two decades. Dr. Joanne Simpson played a central role in GCE modeling developments and applications. She was the lead author or co-author on more than forty GCE modeling papers. In this paper, a brief discussion and review of the application of the GCE model to (1) cloud interactions and mergers, (2) convective and stratiform interaction, (3) mechanisms of cloud-radiation interaction, (4) latent heating profiles and TRMM, and (5) responses of cloud systems to large-scale processes are provided. Comparisons between the GCE model's results, other cloud-resolving model results and observations are also examined.

Description: The Geoscience Laser Altimeter System launched in early 2003 is the first satellite instrument in space to globally observe the distribution of clouds and aerosol through laser remote sensing. The instrument is a basic backscatter lidar that operates at two wavelengths, 532 and 1064 nm. The mission data products for atmospheric observations include the calibrated, observed, attenuated backscatter cross section for cloud and aerosol; height detection for multiple cloud layers; planetary boundary layer height; cirrus and aerosol optical depth and the height distribution of aerosol and cloud scattering cross section profiles. The data is expected to significantly enhance knowledge in several areas of atmospheric science, in particular the distribution, transport and influence of atmospheric aerosol. Measurements of the coverage and height of polar and cirrus cloud should be significantly more accurate than previous global measurement. Initial result from the first several months of operation will be presented.

Description: The April 2002 solar storm event provides a unique opportunity to study the resulting effects on upper atmospheric constituents such as NO and O3. Two sources may perturb these constituents. For a magnetic-storm-source, energetic particles collide with and dissociate N2 in the lower thermosphere to produce excited nitrogen atoms which then combine with O to produce NO. The NO subsequently is transported down to lower altitudes where it reacts with and thus destroys mesospheric and possibly stratospheric ozone. For a solar-particle-event-source, high energy particles penetrate directly into the mesosphere, break apart N2 and water vapor, creating NOx and HOx to destroy ozone in the middle atmosphere. We present perturbations in high northern latitude NO and O3 as measured by the Halogen Occultation Experiment (HALOE) aboard the Upper Atmosphere Research Satellite between April 20-27. HALOE observations show an order of magnitude increase in mesospheric NO and a factor of approximately 2 decrease in mesospheric O3. We will also compare these observations with NASA GSFC 2D model computations.

Description: SAFARI 2000 provided the opportunity to validate MODIS aerosol retrievals and to correct any assumptions in the retrieval process. By comparing MODIS retrievals with ground-based sunphotometer data, we quantified the degree to which the MODIS algorithm underestimated the aerosol optical thickness. This discrepancy was attributed to underestimating the degree of light absorption by the southern African smoke aerosol. Correcting for this underestimation of absorption, produces more realistic aerosol retrievals that allow various applications of the MODIS aerosol products. One such application is the calculation of the aerosol radiative forcing at the top and bottom of the atmosphere. The combination of MODIS accuracy, coverage, resolution and the ability to separate fine and coarse mode make this calculation substantially advanced over previous attempts with other satellites. We focus on the oceans adjacent to southern Africa and use a solar radiative transfer model to perform the flux calculations. The forcing at the top of atmosphere is calculated to be 10 W/sq m, while the forcing at the surface is -26 W/sq m. These results resemble those calculated from INDOEX data, and are most sensitive to assumptions of aerosol absorption, the same parameter that initially interfered with our retrievals.

Description: Stratospheric temperature trends derived from United States meteorological rocketsonde measurements obtained between the late 1960's and mid 1990's are examined at the 50-, 40-, and 25-km altitude levels. Although the trends are different at each of the launch locations there is an unequivocal downward slope of about -0.1 K to -0.3K per year at many of the launch sites. Distances between launch sites and, in some cases unequal data-record lengths, inhibit determination of trend inter-relationship among the various sites. Although the data only provide a 'snapshot' of atmospheric behavior for the specific location, a particular advantage resulting from using these particular rocketsonde observations is their consistency over time. Thus, using the same rocketsonde type over the data period, i.e., Datasonde insures a significant reduction of instrument induced anomalies in the temperature profiles. Trends at the 25-km altitude level, approximately 30 hPa, are compared with trends from radiosonde observations. Both rocketsonde and radiosonde measurements were obtained at approximately the same local times and within less than 100 km of each other. Rocketsonde temperatures from the Former Soviet Union are also used and often complement trends from US data.