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: 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: 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: 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: 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: 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: 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: 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 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: 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.

Description: Recent more accurate antiproton data obtained by the BESS team during the last solar minimum pose a challenge to conventional propagation models of cosmic rays. In particular, the diffusive reacceleration model, which matches well key secondary/primary isotope ratios in cosmic rays, fails to reproduce the secondary antiproton spectrum. Tuning both secondary/primary isotope ratios and antiprotons is possible, but requires artificial breaks in the diffusion coefficient and the injection spectrum of primaries. We will discuss some possibilities to overcome these difficulties in the propagation models. We will present new results of our calculation of CR propagation in the Galaxy using the GALPROP code.

Description: In this talk, five specific major GCE improvements: (1) ice microphysics, (2) longwave and shortwave radiative transfer processes, (3) land surface processes, (4) ocean surface fluxes and (5) ocean mixed layer processes are presented. The performance of these new GCE improvements will be examined. Observations are used for model validation.

Description: Evidence has been found in IMAGE observations that overlap of the plasmasphere and the ring current may lead to enhanced loss of plasma into the ionosphere. It has long been anticipated that this mixing of plasma leads to coupling and resulting consequences on both populations. Wave generation, pitch angle scattering, and heating are some of the consequences that are anticipated. IMAGE plasmasphere ring current, and auroral observations will be presented and used to explore these interactions and their effects.

Description: The Extreme Universe Space Observatory (EUSO) mission will measure the uorescence signal produced by Extensive Air Showers (EAS). To determine the intrinsic shower luminance and hence the energy of the incident cosmic ray, a correction must be made for attenuation in the atmosphere. The correction will be determined on an event-by-event basis by atmospheric sounding with a LIDAR. We report here on a "ground truth" concept to verify EUSO's measurements of intrisic shower luminance.

Description: The Advanced Thin Ionization Calorimeter (ATIC) is a balloon borne experiment to measure the composition and energy spectra of Z = 1 to 28 cosmic rays over the energy range approx. 30 GeV - 100 TeV. The instrument consists of a fully active 320-crystal Bismuth Germanate (BGO) calorimeter, 202 scintillator strips in 3 hodoscopes interleaved with a graphite target, and a 4480-pixel silicon matrix charge detector. ATIC has had two successful Long Duration Balloon flights from McMurdo, Antarctica: from 12/28/00 to 01/13/01 and from 12/29/02 to 01/18/03. We have developed the ATIC Data Processing System (ADPs), which is an Object Oriented data processing program based on ROOT. In this paper, we describe the processing scheme used in handling the flight data, especially the calibration method and the event reconstruction algorithm.

Description: Formed in the gravity regime, complex craters are larger than their simple crater equivalents, due to a combination of slumping and uplift. Just how much larger is a matter of great interest for, for example, age dating studies. We examine three empirical scaling laws for complex crater size, examining their strengths and weaknesses, as well as asking how well they accord with previously published and new data from lunar, terrestrial, and venusian craters.

Description: The goal of this project is to identify the process responsible for the formation of geomorphic features in the Black Point Drive area of Merritt Island National Wildlife Refuge/Kennedy Space Center (MINWR/KSC), northwest Cape Canaveral. This study confirms the principal landscape components (geomorphology) of Black Point Drive reflect interaction between surficial sediments deposited in association with late-Quaternary sea-level highstands and the chemical evolution of late-Cenozoic sub-surface limestone formations. The Black Point Drive landscape consists of an undulatory mesic terrain which dips westward into myriad circular and channel-like depression marshes and lakes. This geomorphic gradient may reflect: (1) spatial distinctions in the elevation, character or age of buried (pre-Miocene) limestone formations, (2) dissolution history of late-Quaternary coquina and/or (3) thickness of unconsolidated surface sediment. More detailed evaluation of subsurface data will be necessary before this uncertain0 can be resolved.

Description: Global Reference Atmospheric Model (GRAM-99) is an engineering-level model of the Earth's atmosphere. It provides both mean values and perturbations for density, temperature, pressure, and winds, as well as monthly- and geographically-varying trace constituent concentrations. From 0-27 km, thermodynamics and winds are based on National Oceanic and Atmospheric Administration Global Upper Air Climatic Atlas (GUACA) climatology. Above 120 km, GRAM is based on the NASA Marshall Engineering Thermosphere (MET) model. In the intervening altitude region, GRAM is based on Middle Atmosphere Program (MAP) climatology that also forms the basis of the 1986 COSPAR Intemationa1 Reference Atmosphere (CIRA). MAP data in GRAM are augmented by a specially-derived longitude variation climatology. Atmospheric composition is represented in GRAM by concentrations of both major and minor species. Above 120 km, MET provides concentration values for N2, O2, Ar, O, He, and H. Below 120 km, species represented also include H2O, O3, N2O, CO, CH, and CO2. Water vapor in GRAM is based on a combination of GUACA, Air Force Geophysics Laboratory (AFGL), and NASA Langley Research Center climatologies. Other constituents below 120 km are based on a combination of AFGL and h4AP/CIRA climatologies. This report presents results of comparisons between GRAM Constituent concentrations and those provided by the Naval Research Laboratory (NRL) climatology of Summers (NRL,/MR/7641-93-7416, 1993). GRAM and NRL concentrations were compared for seven species (CH4, CO, CO2, H2O, N2O, O2, and O3) for months January, April, July, and October, over height range 0-115 km, and latitudes -90deg to + 90deg at 10deg increments. Average GRAM-NRL correlations range from 0.878 (for CO) to 0.975 (for O3), with an average over all seven species of 0.936 (standard deviation 0.049).

Description: The hierarchical organization of microbial ecosystems determines the rates of processes that shape Earth's environment, define the stage upon which major evolutionary events occurred, and create biosignatures in sediments and atmospheres. In cyanobacterial mats, oxygenic photosynthesis provides energy, organic substrates and oxygen to the ecosystem. Incident light changes with depth in the mat, both in intensity and spectral composition, and counteracting gradients of oxygen and sulfide shape the chemical microenvironment. A combination of benefits and hazards of light, oxygen and sulfide promotes the allocation of the various essential mat processes between light and dark periods and to various depths in the mat. Microliters produce hydrogen, small organic acids, nitrogen and sulfur species. Such compounds fuel a flow of energy and electrons in these ecosystems and thus shape interactions between groups of microorganisms. Coordinated observations of population distribution, abundance, and activity for an entire community are making fundamental questions in ecology accessible. These questions address those factors that sustain the remarkable diversity of microorganisms that are now being revealed by molecular techniques. These questions also target the processes that shape the various kinds of biosignatures that we will seek, both in ancient rocks from Earth and Mars, and in atmospheres of distant planets beyond our Solar System.

Description: INTEX-NA is an integrated atmospheric chemistry field experiment to be performed over North America using the NASA DC-8 and P-3B aircraft as its primary platforms. It seeks to understand the exchange of chemicals and aerosols between continents and the global troposphere. The constituents of interest are ozone and its precursors (hydrocarbons, NOX and HOX), aerosols, and the major greenhouse gases (CO2, CH4, N2O). INTEX-NA will provide the observational database needed to quantify inflow, outflow, and transformations of chemicals over North America. INTEX-NA is to be performed in two phases. Phase A will take place during the period of May-August 2004 and Phase B during March-June 2006. Phase A is in summer when photochemistry is most intense and climatic issues involving aerosols and carbon cycle are most pressing, and Phase B is in spring when Asian transport to North America is at its peak. INTEX-NA will coordinate its activities with concurrent measurement programs including satellites (e. g. Terra, Aura, Envisat), field activities undertaken by the North American Carbon Program (NACP), and other U.S. and international partners. However, it is being designed as a 'stand alone' mission such that its successful execution is not contingent on other programs. Synthesis of the ensemble of observation from surface, airborne, and space platforms, with the help of global/regional models is an important It is anticipated that approximately 175 flight hours for each of the aircraft (DC-8 and P-3B) will be required for each Phase. Principal operational sites are tentatively selected to be Bangor, ME; Wallops Island, VA; Seattle, WA; Rhinelander, WI; Lancaster, CA; and New Orleans, LA. These coastal and continental sites can support large missions and are suitable for INTEX-NA objectives. The experiment will be supported by forecasts from meteorological and chemical models, satellite observations, surface networks, and enhanced O3,-sonde releases. In addition to characterizing Atlantic-outflow and Pacific-inflow, INTEX-NA will characterize air masses transported between the U.S., Canada, and Mexico. INTEX-NA will be the first continental scale inflow, outflow, and transformation experiment to be performed over North America. It will provide the most comprehensive observational data set to date to understand the O3/NOX/HOX/aerosol photochemical system and the carbon cycle. One of the critical needs of the carbon cycle research is to obtain large-scale vertical and horizontal concentration gradients of CO2, throughout the troposphere over continental source/sink regions. INTEX-NA is ideally suited to perform this role. Coastal and continental operational sites will allow us to develop a curtain profile of greenhouse gases (e. g. CO2,) and other key pollutants across North America. Such information is central to our quantitative understanding of chemical budgets on the continental scale. We expect to provide a number of satellite under-flights over land and water to test and validate observations from the appropriate satellite platform (e. g. Aura). We plan to develop strong collaborations with other national and international observational programs. Results from INTEX-NA should directly benefit the development of environmental policy for air quality and climate change.

Description: The Living With a Star Geospace Investigations is established to effectively address those phenomena the Geospace environment that directly affect life and society. The priority science questions focus on two broad areas: (1) ionospheric variability, especially at mid-latitudes, that affects navigation and communications and (2) the source, acceleration mechanisms, and sinks of the radiation belts that degrade satellite lifetimes, produce surface charging, and threaten manned space flight. Candidate missions to address these science foci will be presented as well as possible additional investigations and experiments that would enable an understanding of the Geospace at the system level.

Description: An evolution of ionospheric convection was studied for a double transpolar arc phenomenon on February 11, 1999. While one transpolar arc split from the auroral oval in the morning sector and drifted duskward, another arc appeared in the evening sector. The convection was investigated with three velocity data sets: E B drift velocities from the ASTRID-2 satellite; Ion Driftmeter data from the DMSP satellites; and Doppler-shift data from the Super-DARN radars. We inferred convection cells from these data sets and found that the number of convection cells changed from three to four as the dominance of IMF changed from a negative By to a positive Bz. Our result suggests that the ionospheric convection that has been so far discussed for various conditions of IMF may be applied even to the cases accompanied by transpolar arcs.

Description: Magnetosphere-Ionosphere (MI) coupling and associated with this process electron and ion energization processes have interested scientists for decades and, in spite of experimental and theoretical research efforts, are still ones of the least well known dynamic processes in space plasma physics. The reason for this is that the numerous physical processes associated with MI coupling occur over multiple spatial lengths and temporal scales. One typical example of MI coupling is large scale ring current (RC) electrodynamic coupling that includes calculation of the magnetospheric electric field that is consistent with the ring current (RC) distribution. A general scheme for numerical simulation of such large-scale magnetosphere-ionosphere coupling processes has been presented earlier in many works. The mathematical formulation of these models are based on "modified frozen-in flux theorem" for an ensemble of adiabatically drifting particles in the magnetosphere. By tracking the flow of particles through the inner magnetosphere, the bounce-averaged phase space density of the hot ions and electrons can be reconstructed and the magnetospheric electric field can be calculated such that it is consistent with the particle distribution in the magnetosphere. The new a self-consistent ring current model has been developed that couples electron and ion magnetospheric dynamics with calculation of electric field. Two new features were taken into account in addition to the RC ions, we solve an electron kinetic equation in our model, self-consistently including these results in the solution. Second, using different analytical relationships, we calculate the height integrated ionospheric conductances as the function of precipitated high energy magnetospheric electrons and ions as produced by our model. This results in fundamental changes to the electric potential pattern in the inner magnetosphere, with a smaller Alfven boundary than previous potential formulations would predict but one consistent with recent satellite observations. This leads to deeper penetration of the plasma sheet ions and electrons into the inner magnetosphere and more effective ring current ions and electron energization.

Description: The ALTUS Cumulus Electrification Study (ACES) is an uninhabited aerial vehicle (UAV)-based project that will investigate thunderstorms in the vicinity of the Florida Everglades in August 2002. ACES is being conducted to both investigate storm electrical activity and its relationship to storm morphology, and validate Tropical Rainfall Measurement Mission (TRMM) satellite measurements. In addition, as part of NASA's UAV-based science demonstration program, this project will provide a scientifically useful demonstration of the utility and promise of UAV platforms for Earth science and applications observations. Part of the demonstration involves getting approvals from the Federal Aviation Administration and the NASA airworthiness flight safety review board. ACES will employ the ALTUS II aircraft, built by General Atomics - Aeronautical Systems, Inc. Key science objectives simultaneously addressed by ACES are to: (1) investigate lightning-storm relationships, (2) study storm electrical budgets, and (3) provide Lightning Imaging Sensor validation. The ACES payload, already developed and flown on ALTUS, includes electrical, magnetic, and optical sensors to remotely characterize the lightning activity and the electrical environment within and around thunderstorms. ACES will contribute important electrical and optical measurements not available from other sources. Also, the high altitude vantage point of the UAV observing platform (up to 55,000 feet) offers a useful 'cloud-top' perspective. By taking advantage of its slow flight speed (70 to 100 knots), long endurance, and high altitude flight, the ALTUS will be flown near, and when possible, above (but never into) thunderstorms for long periods of time, allowing investigations to be conducted over entire storm life cycles. In addition, concurrent ground-based observations will enable the UAV measurements to be more completely interpreted and evaluated in the context of the thunderstorm structure, evolution, and environment.

Description: Stratospheric water vapor is important not only for its greenhouse forcing, but also because it plays a significant role in stratospheric chemistry. Several recent studies have focused on the potential for dehydration due to ice cloud formation in air rising slowly through the tropical tropopause layer (TTL). Holton and Gettelman showed that temperature variations associated with horizontal transport of air in the TTL can drive ice cloud formation and dehydration, and Gettelman et al. recently examined the cloud formation and dehydration along kinematic trajectories using simple assumptions about the cloud properties. In this study, a Lagrangian, one-dimensional cloud model has been used to further investigate cloud formation and dehydration as air is transported horizontally and vertically through the TTL. Time-height curtains of temperature are extracted from meteorological analyses. The model tracks the growth, advection, and sedimentation of individual cloud particles. The regional distribution of clouds simulated in the model is comparable to the subvisible cirrus distribution indicated by SAGE II. The simulated cloud properties and cloud frequencies depend strongly on the assumed supersaturation threshold for ice nucleation. The clouds typically do not dehydrate the air along trajectories down to the temperature minimum saturation mixing ratio. Rather the water vapor mixing ratio crossing the tropopause along trajectories is 10-50% larger than the saturation mixing ratio. I will also discuss the impacts of Kelvin waves and gravity waves on cloud properties and dehydration efficiency. These simulations can be used to determine whether observed lower stratospheric water vapor mixing ratios can be explained by dehydration associated with in situ TTL cloud formation alone.

Description: Empirical models of plasmaspheric properties date from the pioneering work of Storey where he developed the analysis of ground whistler observations that lead to his estimate for the equatorial plasma density at L=3. The most recent in situ satellite study takes us to 1000 CRRES satellite passes and a statistical analysis of the plasmapause location at all local times and for varying geomagnetic conditions by Moldwin et al. These and many other studies over the intervening 49 years have given us a strong familiarity with the distribution of cold plasmaspheric ions throughout the magnetosphere. The major components of inner plasmasphere, nightside bulge, sunward convection tail, and plasmapause are all well established. Storm-time erosion and the resulting ionospheric refilling has been encompassed, even if not completely understood. Small-scale density variations near the plasmapause and extending at least to geosynchronous orbit have been characterized in a variety of ways, even though we do not yet understand their origin. This paper will present early empirical modeling results from the inversion of IMAGE/EW global intensity images to density distributions. Densities are obtained in this initial study through use of forward image modeling with a simple 3-parameter plasmaspheric and plasmapause mathematical model. Individual interior plasmaspheric density profiles and plasmapause locations are obtained every 10 degrees in magnetic local time for each E W image analyzed. Derived profile parameters are statistically characterized in the context of storm magnitude and evolution. Identified patterns in the appearance of plasmaspheric structures, plasmapause erosion, and refilling will be presented. Comparisons to existing empirical plasmaspheric models and the implications for new modeling will be presented. Additional information is included in the original extended abstract.

Description: Water vapor in the winter arctic tropopause region is important because, after the tropical tropopause region, the winter arctic tropopause has the coldest temperatures in the tropospheric northern hemisphere. This suggests the potential for cloud formation that can remove water vapor from a part of the atmosphere where radiatively active gases (such as water) exert a disproportionate influence on the earth's radiation budget. Previous work by the same authors has shown that this cloud formation extends into the stratosphere, with 20% of the parcels having ozone values of 300-350 ppbv experiencing ice saturation in any given 10 day period period during the late winter. In fact, temperatures are cold enough that 5-10% of the parcels experience saturation even if the water content is below the prevailing stratospheric value of 5 ppmv. This work describes a case study of clouds observed by aircraft near the winter arctic tropopause during the SAGE Ozone Loss and Validation Experiment (SOLVE). This provided a unique opportunity to examine dehydration processes in this region since in situ water, tracer, cloud particle, and meteorological data were all available simultaneously. During this period, temperatures were cold enough at the tropopause to produce saturation mixing ratios of 3-4 ppmv. Thus, clouds were actually observed within the stratosphere. Back trajectories indicate that the air in these clouds came from lower latitudes and altitudes. The study describes the nature of the clouds, the history of the air, and the possible implications for the upper tropospheric water budget.

Description: Spacecraft designers and operations support personnel involved in space environment analysis for low Earth orbit missions require ionospheric specification and forecast models that provide not only average ionospheric plasma parameters for a given set of geophysical conditions but the statistical variations about the mean as well. This presentation describes the development of a prototype empirical model intended for use with the International Reference Ionosphere (IRI) to provide ionospheric Ne and Te variability. We first describe the database of on-orbit observations from a variety of spacecraft and ground based radars over a wide range of latitudes and altitudes used to obtain estimates of the environment variability. Next, comparison of the observations with the IRI model provide estimates of the deviations from the average model as well as the range of possible values that may correspond to a given IRI output. Options for implementation of the statistical variations in software that can be run with the IRI model are described. Finally, we provide example applications including thrust estimates for tethered satellites and specification of sunrise Ne, Te conditions required to support spacecraft charging issues for satellites with high voltage solar arrays.

Description: Tropospheric measurements of ozone from SAGE II (version 6.1) in the tropics have been analyzed using 12 years of data (1985-1990, 1994-1999). The seasonally averaged vertical profiles of the ozone mixing ratio in the upper troposphere have been presented for the first time from satellite measurements. These profiles show qualitative similarities with corresponding seasonal mean ozonesonde profiles at northern and southern tropical stations and are about 40-50% less than the sonde values. Despite this systematic offset, the measurements appear to be consistent with a zonal wave one pattern in the upper tropospheric column ozone and with the recently predicted summertime ozone enhancement over the Middle East. These results thus affirm the usefulness of the occultation method in studying tropospheric ozone.

Description: Long term X-ray monitoring data from the RXTE All Sky Monitor (ASM) reveal that the third (superorbital) period in SMC X-1 is not constant but varies between 40-60 days. A dynamic power spectrum analysis indicates that the third period has been present continuously throughout the five years of ASM observations. This period changed smoothly from 60 days to 45 days and then returned to its former value, on a timescale of approximately 1600 days. During the nearly 4 years of overlap between the CGRO & RXTE missions, the simultaneous BATSE hard X-ray data confirm this variation in SMC X-1. Sources of systematic error and possible artefacts are investigated and found to be incapable of reproducing the results reported here. Our disco cry of such an instability in the superorbital period of SMC X-1 is interpreted in the context of recent theoretical studies of warped, precessing accretion discs. We find that the behaviour of SMC X-1 is consistent with a radiation - driven warping model.

Description: The Upper Atmosphere Research Satellite (UARS) was launched in September 1991 by the Space Shuttle Discovery. Seven of the original ten instruments aboard the UARS are functional and six instruments regularly take measurements. The UARS is in a stable observing configuration, in spite of experiencing several anomalies which have impacted the data tape storage and power available. Power sharing and "real-time" operations using the Tracking and Data Relay Satellites have reduced the data collection over the past five years. Although the UARS measurements have not been continuous, the UARS instruments HALOE, SUSIM, SOLSTICE, HRDI, WINDII, and PEM have provided important observations over an entire solar cycle. The UARS HALOE instrument measures ozone and other constituents important in understanding ozone variations. The UARS SUSIM and SOLSTICE instruments observe ultraviolet light between 120 and 420 nm, which influence middle atmospheric constituent fluctuations. The UARS HRDI and WINDII instruments measure winds in the stratosphere, mesosphere, and thermosphere, which move constituents among the various atmospheric regions. The UARS PEM instrument provides observations of input precipitating charged particles, including both electrons and protons. This paper will provide an overview about the UARS mission and its relevant atmospheric measurements.

Description: Airborne measurements of a large number of oxygenated organics were carried out in the Pacific troposphere (to 12 km) in the Spring of 2001 (Feb. 24-April 10). Specifically these measurements included acetaldehyde, propanaldehyde, acetone, methylethyl ketone, methanol, ethanol, PAM and organic nitrates. Independent measurements of formaldehyde, peroxides, and tracers were also available. Highly polluted as well as pristine air masses were sampled. Oxygenated organics were abundant in the clean In troposphere and were greatly enhanced in the outflow regions from Asia. Extremely high concentrations of aldehydes could be measured in the troposphere. It is not possible to explain the large abundances of aldehydes in the background troposphere without invoking significant oceanic sources. A strong correlation between the observed mixing ratios of formaldehyde and acetaldehyde is present. We infer that higher aldehydes (such as acetaldehyde and propanaldehyde) may provide a large source of formaldehyde and sequester Cox throughout the troposphere. The atmospheric behavior of acetone, methylethyl ketone, and methanol is generally indicative of their common terrestrial sources with a Image contribution from biomass/biofuel burning. A vast body of data has been collected and it is being analyzed both statistically and with the help of models to better understand the role that oxygenated organics play in the atmosphere and to unravel their sources and sinks. These results will be presented.

Description: We report the first in-situ measurements of hydrogen cyanide (HCN) and acetonitrile (CH3CN) from the Pacific troposphere (0-12 km) obtained during the NASA/Trace-P mission (Feb.-April, 2001). Mean HCN and CH3CN mixing ratios of 243 (+/-118) ppt and 149 (+/-56) ppt respectively, were measured. The in-situ observations correspond to a total HCN column of 4.4-4.9 x 10(exp 15) molec. cm(exp -2) and a CH3CN column of 2.8-3.0 x 10(exp 15) molec. cm(exp -2). This HCN column is in good agreement with available spectroscopic observations. The atmospheric concentrations of HCN and CH3CN were greatly influenced by outflow of pollution from Asia. There is a linear relationship between the mixing ratios of HCN and CH3CN, and in turn these are well correlated with tracers of biomass combustion (e.g. CH3Cl, CO). Relative enhancements with respect to known tracers of biomass combustion within selected plumes in the free troposphere, and pollution episodes in the boundary layer allow an estimation of a global biomass burning source of 0.8+/-0.4 Tg (N)/y for HCN and 0.4+/-0.1 Tg (N)/y for CH3CN. In comparison, emissions from automobiles and industry are quite small (〈0.05 Tg (N)/y). The vertical structure of HCN and CH3CN indicated reduced mixing ratios in the MBL (Marine Boundary Layer). Using, a simple box model, the observed gradients across the top of the MBL are used to derive an oceanic flux of 6.7 x 10(exp -15) g (N) cm(exp -2)/s for HCN and 4.8 x 10(exp -15) g (N) cm(exp -2)/s for CH3CN. An air-sea exchange model is used to conclude that this flux can be maintained if the oceans are under-saturated in HCN and CH3CN by 23% and 17%, respectively. It is inferred that oceanic loss is a dominant sink for these nitrites, and they deposit some 1.3 Tg (N) of nitrogen annually to the oceans. Assuming reaction with OH radicals and loss to the oceans as the major removal processes, a mean atmospheric residence time of 4.7 months for HCN and 5.1 months for CH3CN is calculated. A global budget analysis shows that the sources and sinks of HCN and CH3CN are roughly in balance. There are indications that biogenic sources may also be present. Mechanisms involved in nitrate formation during combustion and removal in the oceans are poorly understood.