ALBERT

Description: The He I 1083 nm multiplet is a powerful tool for observing the outer solar atmosphere but is difficult to analyze because the lines are weak, highly variable, and spectrally blended, both internally and with other neighboring solar and telluric lines. After separation from nearby spectral features, two components of the He I multiplet are resolved. Fitting these lines with two unconstrained Gaussian profiles always gives a ratio of major to minor component of less than half the value which would be expected for optically thin lines. One possibility for explaining the discrepancy between the weakness of the lines and the ratio of the spectral components is to assume that the line formation regions are concentrated in laterally unresolved, optically thick structures with small filling factor. However, we present here a least-squares fitting technique using cubic splines with fixed breakpoints with the constraint that the blend is the sum of three identically shaped profiles shifted in wavelength according to the atomic structure of the blend and weighted by the corresponding statistics weights, in agreement with optically thin line formation. The basis functions for the fitting procedure have no built-in spectral symmetry or shape. The resulting underlying profiles tend to be asymmetric with excess absorption to the blue, consistent with formation by "hot" and "cool" spatial elements within the observed volume, with the hotter regions having differential motion toward the observer. The results agree well with NASA/XSO Spectromagnetograph observations in quiet sun and coronal holes. Partial funding of this research was provided through the NASA Sun-Earth Connection SR&T program.

Description: GeoSpec will support several possible future mission concepts in the Atmospheric Sciences and in Land and Ocean Sciences by measurement of both chemically linked atmospheric trace gas concentrations and profiles of important molecules such as OS, N02, CH20 and SO2 and at the same time coastal and ocean pollution events, tidal effects, and the origin and evolution of aerosol plumes. The instrument design concept we will describe is a dual spectrograph covering the WMS wavelength region of 310- 481 nm and the VIS/NIR wavelength region of 500-900 nm. A third channel in the short- wave infrared (SWIR) region between 2.2 p and 2.4 pn for total column measurements of CO, CH4, and N20 will also be described. The goal is to design a system capable of making moderate spatial resolution (750 meters at nadir) hyperspectral measurements (0.2 to 1.2 nm resolution) from a geostationary orbit. This would enable studies of time- varying pollution and coastal change processes with a temporal resolution of 5 minutes on a regional scale to 1 hour on a continental scale. Technological advances in the design and fabrication of convex holographic gratings and large format, high dynamic range PIN/CMOS detectors at the focal plane will be exploited. By simply changing the focal length of the front-end telescope GeoSpec can accommodate different orbital altitudes, including low Earth orbit, the Sun-side Lagrangian point L1, and/or different spatial resolutions.

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 J2, the Earth's dynamic oblateness. J2 has long been observed to undergo a slight decrease due primarily 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 turning back, signifymg a large, anomalous mass redistribution whose J2 effect overshadows that of the post-glacial rebound over interannual timescales. Intriguing evidences have been found in the extratropical Pacific basins, especially related to the Pacific Decadal Oscillation, and in land hydrology. We will examine the latest results based on ocean altimetry, sea-surface temperature, and ocean and hydrology model outputs. Besides 52, 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, some of these series have significant signal that show correlation to various climatic signals. For example, the SLR-observed S2,l tesseral (and 53 zonal) coefficients compare favorably with those computed from the NCEP atmospheric pressure field. There is a significant correlation of the sectoral S2,2 with the Southern Oscillation Index signifying the influence of El NinoLa Nina, but preceding SO1 by about 1 year. Cases such as these demonstrate the utility of assessing the mass component of climate variations.

Description: In late February and March, 2003, the Ice, Cloud, and land Elevation Satellite (ICESat) measured ice and land elevations along profiles across southern Alaska. During this initial data acquisition stage ICESat observations were made on 8-day repeat tracks to enable calibration and validation of the ICESat data products. Each profile consists of a series of single point values derived from centroid elevations of an $\approx$70 m diameter laser footprint. The points are s4pakated by $\approx$172 m along track. Data siets of 8-day observations (an ascending and descending ground track) crossed the Bering and Malaspina Glacier. Following its 1993--1995 surge; the Bering Glacier has undergone major terminus retreat as well as ike thinning in the abtation zone. During the later part of the 20th century, parts of the Malaspina thinned by about 1 m/yr. The multiple observation profiles across the Bering and Malaspina piedmont lobes obtained in February/March are being geolocated on Landsat images and the elevation profiles will be used for a number o scientific objectives. Based on our simulations of ICESat performance over the varied ice surface of the Jakobshavn Glacier of GReenland, 2003, we expect to measure annual, and possibly seasonal, ice elevation changes on the large Alaskan glaciers. Using elevation data obtained from a second laser, we plan to estimate ice elevation changes on the Bering Glacier between March and October 2003.

Description: Over two decades of geodetic satellite-laser-ranging (SLR) data show that the variation of the Earth's oblateness parameter J2 has a clear seasonal signal of amplitude of about 3e-10 and a secular decrease of about -2.8e-11/year, superimposed on some interesting interannual fluctuations. Physically, any change in mass distribution or/inside the Earth will be reflected in the time-variable gravity signal obtained outside the Earth, according to Newton s gravitational law. Therefore, such signal contains contributions from all geophysical sources that redistribute mass, on all temporal and spatial scales, including those from the core. Besides Earth rotation and geomagnetic field variations, the time-variable gravity also contains information linking Earth surface observations with internal core dynamical processes. The time scales of the gravity signal are critical in helping differentiate different contributions. The atmosphere and hydrosphere are responsible for the seasonal and much of the interannual and intraseasoanl fluctuations, while the secular trend is due mainly to the post-glacial rebound but possibly core mass flow. To estimate the latter effect, we use our MoSST (Modular, Scalable, Self-consistent, Three-dimensional) core dynamics model to forward simulate the core flow, and density variation due to the core convection. Our results suggest that, when upward continued to the surface, the J2 component of the core mass redistribution can reach an overall amplitude of e-11/year, approaching the SLR detectability and significant in geophysical terms. We also find a general westward drift of the mass flow, with a speed comparable to that of the geomagnetic westward drift.

Description: This presentation will discuss the sensitivity of assimilated ozone fields in the upper troposphere and lower stratosphere (UTLS) to a number of factors, focusing mainly on aspects of data selection and the prediction model. This is important, because assimilation represents an attempt to construct our best estimates of the true ozone field; however, inaccuracies in the UTLS ozone distribution translate into an uncertainty in factors such as the calculated radiative forcing of climate or the inferred stratosphere-troposphere exchange (STE) of ozone. The 3D ozone data assimilation system, from NASA's Global Modeling and Assimilation Office (GMAO), combines observations of total ozone column and stratospheric profiles with predictions from an off-line, parameterized chemistry and transport model (pCTM) to produce six-hourly, global analyses. The first experiments discussed assimilate ozone retrievals from the Earth-Probe Total Ozone Mapping Spectrometer (EPTOMS) and stratospheric profiles from the Solar Backscatter UltraViolet/2 (SBUV/2) instrument. The SBUV/2 ozone data have a coarse vertical resolution, with increased uncertainty below the ozone maximum, and TOMS provides only total ozone columns. Thus, the assimilated ozone profiles in the UTLS region are only weakly constrained by the incoming SBUV and TOMS data. Consequently, the assimilated ozone distribution should be sensitive to changes in inputs to the statistical analysis scheme. Sensitivity studies have been conducted to examine the responses to TOMS and SBUV/2 data selection, modifications of the forecast and observation error covariance models, and the model formulation (turning off chemistry or using different wind analyses in the pCTM). The second set of experiments includes an additional data type: ozone retrieved from infrared limb-emission by MIPAS on Envisat. These data offer not only improved vertical resolution in the stratosphere, but also give measurements in the polar night. Comparisons of the assimilated ozone fields from both sets of experiments with independent observations, primarily ozone sondes, are used to determine the impact of each of these changes. It is shown that many of the changes have a significant impact on the UTLS ozone estimates. Implications for interpretation of STE and radiative forcing of climate are discussed.

Description: A network of 12 southern hemisphere tropical and subtropical stations in the Southern Hemisphere ADditional OZonesondes (SHADOZ) project has provided over 2000 profiles of stratospheric and tropospheric ozone since 1998. Balloon-borne electrochemical concentration cell (ECC) ozonesondes are used with standard radiosondes for pressure, temperature and relative humidity measurements. The archived data are available at:http: //croc.gsfc.nasa.gov/shadoz. In Thompson et al., accuracies and imprecisions in the SHADOZ 1998- 2000 dataset were examined using ground-based instruments and the TOMS total ozone measurement (version 7) as references. Small variations in ozonesonde technique introduced possible biases from station-to-station. SHADOZ total ozone column amounts are now compared to version 8 TOMS; discrepancies between the two datasets are reduced 2\% on average. An evaluation of ozone variations among the stations is made using the results of a series of chamber simulations of ozone launches (JOSIE-2000, Juelich Ozonesonde Intercomparison Experiment) in which a standard reference ozone instrument was employed with the various sonde techniques used in SHADOZ. A number of variations in SHADOZ ozone data are explained when differences in solution strength, data processing and instrument type (manufacturer) are taken into account.

Description: Besides operating its own energy budget, an earthquake acts as an agent transferring a much greater amount of energy among the Earth's rotation, elastic field, gravitational field and internal heat. We compute the co-seismic, globally integrated gravitational and rotation changes induced by some 20,000 large earthquakes that occurred in the last quarter century, according to Chao et al. (1995, GJI, 122,776- 783,784-789) and using the Harvard CMT catalog. The result confirms an extremely strong tendency for the earthquakes to decrease the global gravitational energy and to increase the spin energy. It is found that energy is being extracted from the Earth's gravitational field by the action of earthquakes at an average rate of about approx. 2 TeraW during the studied period, larger by far than the approx. 7 GigaW for the average rate of the earthquake-induced rotational energy increase and the approx. 5 GigaW for the seismic energy release. Based on energetics considerations and assuming the inability of the Earth to build up elastic energy continuously over time, it is argued that earthquakes, by converting gravitational energy, may make a significant contribution to the global hedflow.

Description: The report present the atmosphere, ocean color, land and solar irradiation data sets. The data presented: total ozone, aerosol, cloud optical and physical parameters, temperature and humidity profiles, radiances, rain fall, drop size distribution.

Description: Water vapor measurements from the GSFC Stratospheric Ozone Lidar were made for the first time during a campaign at NOAA's Mauna Loa Observatory. Comparisons were made among the GSFC lidar, the NOAA Lidar and water vapor sondes which were flown from the observatory at times coincident with the lidar measurements.

Description: NASA's Ice, Cloud and Land Elevation Satellite (ICESat) mission is designed to measure changes in elevation of the Greenland and Antarctic ice sheets beginning in January 2003. Time-series of elevation changes will enable determination of the present- day mass balance of the ice sheets, study of associations between observed ice changes and polar climate, and estimation of the present and future contributions of the ice sheets to global sea level rise. Other scientific objectives of ICESat include: global measurements of cloud heights and the vertical structure of clouds and aerosols; precise measurements of land topography and vegetation canopy heights; and measurements of sea ice roughness, sea ice thickness, ocean surface elevations, and surface reflectivity. The Geoscience Laser Altimeter System (GLAS) on ICESat has a 1064 nm laser channel for surface altimetry and dense cloud heights and a 532 nm lidar channel for the vertical distribution of clouds and aerosols. Differences between the characteristics of laser and radar altimetry, such as effective depth of the backscattered signal, elevation accuracy, and footprint location, and their relevance to inter-relating measurements from ERS, Envisat, ICESat, and Cryosat are discussed. Preliminary ICESat results obtained during the calibration and validation period of ICESat are described. ICESat is designed to operate for 3 to 5 years and should be followed by successive missions to measure ice changes for at least 15 years.

Description: Over the past several years the Pennsylvania Departments of Environmental Protection (DEP), Health (DOH), and Agriculture (PDA) built the GIs-based Pennsylvania West Nile Surveillance System. That system has become a model for collecting data that has a field component, laboratory component, reporting and mapping component, and a public information component. Given the success of the West Nile Virus System and the events of September 11, 2001, DEP then embarked on the development of the Pennsylvania Incident Response System, or PAIRS. PAIRS is an effective GIs-based approach to providing a system for response to incidents of any kind, including terrorism because it is building upon the existing experience, infrastructure and databases that were successfully developed to respond to the West Nile Virus by DEP, DOH, and PDA. The proposed system can be described as one that supports data acquisition, laboratory forensics, decision making/response, and communications. Decision makers will have tools to view and analyze data from various sources and, at the same time, to communicate with the large numbers of people responding to the same incident. Recent collaborations with NASA partners are creating mechanisms for the PAIRS system to incorporate space-based and other remote sensing geophysical parameters relevant to public health assessment and management, such as surface temperatures, precipitation, land cover/land use change, and humidity. This presentation will describe the PAIRS system and outline the Pennsylvania-NASA collaboration for integration of space-based data into the PAIRS system.

Description: Simulations were conducted to investigate the influence of rapid electric field fluctuations on electron energization in the inner magnetosphere based on the assimilative mapping of ionospheric electrodynamics (AMIE) technique. Simulations for four different magnetic storms were run, namely those that occurred on May 15,1997, May 4, 1998, September 25, 1998, and October 19, 1998. Here, we have examined the formation of high energy electrons in the inner magnetosphere during these storm events with our recently-developed relativistic radiation belt transport code. The point of this numerical experiment is to show that a simulation of a real event must have the high time resolution electric field input files in order to produce the seed population for the radiation belts, which are often observed to increase in the days following a magnetic storm. Specifically, a cadence of the global electric field pattern of 5 minutes or less produces inner magnetospheric fluxes that are larger (by up to 5 orders of magnitude) than fluxes produced with a longer cadence. Differences were particularly large relative to simulation results with a 3-hour time cadence, analogous to a Kp-driven electric field model.

Description: The Magnetospheric Multiscale Mission (MMS) is a NASA mission intended to make fundamental advancements in our understanding of the Earth's Magnetosphere. There are three processes that MMS is intended to study including magnetic recon- nection, charged particle acceleration, and turbulence. There are four phases of the MMS mission and each phase is designed to study a particular region of the Earth's magnetosphere. The mission is composed of a formation of four spacecraft that are nominally in a regular tetrahedron formation. In this work, we present optimal orbit designs for Phase I and II. This entails designing optimal reference orbits so that the spacecraft dwell-time in the region of interest is a maximum. This is non-trivial because the Earth's magnetosphere is dynamic and its shape and position are not constant in inertial space. Optimal orbit design for MMS also entails designing the formation so that the relative motion of the four spacecraft yields the greatest science return. We develop performance metrics that are related to the science return, and use Sequential Quadratic Programming (SQP) to determine optimal relative motion solutions. We also ensure that practical constraints such as maximum eclipse time and minimum inter-spacecraft separation distances are not violated.

Description: At altitudes near 5000 km over the Southern polar cap region of the terrestrial magnetospherehonosphere, the Thermal Ion Dynamics Experiment (TIDE) onboard the Polar satellite has observed O+ ion density trough regions, in which the O+ densities were at least one order of magnitude lower than the surrounding O+ densities. In the O+ demify trough regions, the estimated O+ densities were generally lower than 0.01 per cc. The boundaries between normal density level regions and the trough density regions were usually abrupt transitions. From December 1, 1997 to November 30, 1998, polar cap O+ troughs in Polar/TIDE observations occurred at a frequency of about 48%. Statistical examination of the Polar perigee observations from December 1 , 1997 to November 30, 1998 shows that the Polar perigee passes evenly covered the southern polar cap region, while the O+ density trough was always located on the nightside portion of the polar cap magnetospherehonosphere, and that invariant latitude spans of such troughs could be as large as 230 in extent. The trough occurrence displayed strong seasonal dependence; in the winter season (e.g. for July in the southern hemisphere) the O+ ion density trough occurrence frequency ranged up to 92%, while in the summer season (e.g. for January in the southern hemisphere) it decreased to as infrequent as 15%. The O+ ion density trough occurrence appeared relatively independent of the geomagnetic Kp index, and IMF Bz, By conditions. However, as suggested by the seasonal dependence, the O+ ion density trough occurrence was strongly related to the solar zenith angle (SZA). In the SZA range 500 to 1250, the trough occurrence increased monotonically with SZA. Also, case-by-case examinations of near-simultaneous O+ densities and vertical velocities observed by the DMSP satellite group orbiting at 840 km altitude indicate that the O+ density troughs observed at 5000 km altitude exhibit moderate correlation or anti-correlation with topside ionosphere density or velocity features at 840 km altitude for approximately half of the near-conjunction cases examined.

Description: Magnetic flux tubes containing plasmaspheric ion density and composition have been observed in the region between the classical plasmapause and the magnetopause. New observations show that these ion distributions exist at the equatorial, post-noon magnetopause. Comparison to observations of similar distributions at geosynchronous orbit and to simulations leads to the conclusion that these ions are convected from these regions to the magnetopause. This represents an extension of the geopause to the outer edge of the magnetosphere on the dayside. The presence of ion densities 〉 10 cu cm in this region must have profound impact on the nature of plasma processes that occur there.

Description: We present EUV observations of the plasmasphere-plasmapause from 19:38-22:11 UT on 28 June 2000 characterized by the presence of bifurcated radial enhancements of the He(+) plasma distribution in the nighside sector. These features remain stable throughout the period of observation and are found to co-rotate at 67% of the expected rate. Two-dimensional simulation of the plasmasphere assuming the presence of field lines resonances at L = 1.8 and 2.5 suggest that the organization of the outer plasmasphere and plasmapause is the result convective motion driven by a standing ULF-wave. Preliminary analysis of ground-based magnetometer data provided by the IMAGE magnetometer network during the period of EUV observation indicates the presence of a discrete spectrum of field line resonances extending down to 0.68-mHz.

Description: We have produced organic material simulating a methane photochemical haze in a CO2- rich atmosphere of the early Earth by irradiating gas mixtures in an inductively coupled cold plasma chamber with pressure approx. 0.25 mbar at 100 W total power. The flow rate was 24 cm3 min. We added progressively higher levels of CH, by combining gas mixtures of N2/CH4 (9/1) and N2/CO2 (9/1) to change the ratio of CH4/CO2. Tholin was accumulated for 5 hours in each experiment; the onset of tholin formation is in the range CH4/CO2 = 0.5 to 1. As the mixing ratio of CH, is increased, the production rate of the brownish tholin film increases. IR spectra showed the C-H and N-H bands similar to that of Titan tholin and closely resemble Titan tholin made at 0.13 mbar pressure. A decrease in the CH bonds on decreasing CH4/CO2 is noted. Ether bands (-(2-O-C) were tentatively detected, but no detectable carbonyl (C=O) band was found. The absorption in the UV region for the early Earth tholin is found to be substantially greater than the Titan tholin. Quantitative values of the optical constants of early Earth tholin are currently being measured.

Description: Between 2400A to 3160A, we have identified approximately 500 absorption line complexes, each with up to 20 velocity components. Lines of Fe I, Fe II, Ti II, V II, Ni II, Co II, Mn II, Mg I, Mg II and Na I have been identified. Surprisingly, most of the lines originate from energy levels significantly above the ground level. This is indcative of optical pumping from the Central Source. Line widths and population of various levels are non-thermal. The relative column densities change with velocity. For example, Fe 11 column densities for one transition arising from approx. 0.l ev increase with ejecta velocity while Fe 11 column densities for another transition decreases with velocity. This may be due to softening of the ultraviolet radiation that pumps the various ions (neutrals) with velocity. If we assume that the distance of each system scales with distance from the Central Source, only a thirty percent change in distance is noted; yet the ratio of column densities for the Fe I1 examples given above changes by nearly thirty-fold. If the ejecta distance scales with distance from Eta Carinae, then it is likely that this ejecta originated at nearly the same time. The geometry of the Homunculus has been determined to be a double-lobed structure tilted out of the plane of the sky. We interpret the ejecta as being in the wall of the Southwest lobe, and that this wall just happens to be in line of sight from Eta Carinae to the observer. As Eta Carinae enters into the upcoming minimum, we are already seeing some evidence for changes in column densities due to changes in ultraviolet fluxes. This is reinforced by IUE observations that we have recently re-analyzed with respect to the 5.52 year (2020 +/- 10 days) spectroscopic period. Observations were done through STScI and funding was through the STIS GTO resources.

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 date. 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. In addition, although less precise than GRACE, the GPS/Meteorology constellation mission COSMIC, with 6 mini-satellites to be launched in late 2005, is expected to provide continued and complementary time-variable gravity observations. Such observations are becoming a new and powerful tool for remote sensing of geophysical fluid processes that involve larger-scale mass transports.

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: The unique instrumentation on the Polar satellite combined with the simultaneous measurement of different parts of the magnetosphere with multiple satellites make possible the study of magnetospheric processes in a special way. In particular, the study of the ionospheric supply of plasma to the magnetosphere can by accomplished to give important results on the plasmas which drive the dynamics of the magnetosphere. This study concentrated on the period of September to December, 2001 in which the Polar orbit had precessed to the point that the line of apsides was near the equatorial plane. This unique orbital configuration enabled the TIDE instrument to measure outflowing ions across the polar cap and then transit the magnetotail lobes and observe the dramatic change in plasma characteristics as the satellite entered the plasma sheet. Contact was made with investigators on the Cluster and Geotail satellite missions and corresponding time frames were studied in the data. There were two approximate conjunctions between Polar and Geotail and data were compared to look for features which might be related. The higher energy concentration of the Geotail instrument made direct comparisons with TIDE difficult, and the Cluster measurements did not surface any cases that corresponded closely in space and time. There were, however, many interesting aspects of the Polar orbits which permitted the observation of the changing ionospheric outflowing plasma characteristics. As in earlier measurements, the ionospheric plasma could be seen flowing up the magnetic field lines out of the northern and southern polar caps. Its energy suggested a polar wind origin energized by the centrifugal acceleration of flow through the polar cusp. The roughly 10eV ions then moved out into the lobes of the magnetotail where they could be seen flowing toward the plasma sheet in both the northern and southern magnetotail lobes. There was also a double field-aligned region of warm ions observed just outside the plasmasphere, stretching toward the auroral zone and inner plasma sheet boundary. Upon entering the plasma sheet, the plasma energy jumped from 10 s of eV to greater than 1 keV. The single field aligned flows transitioned to highly variable spatially choppy energized ion distributions with a variety of pitch angle configurations. This pattern was quite repeatable in all of the Polar orbits that were examined and are compatible with the source of ions being the polar wind which is then energized to 10 s of eV by the centrifugal acceleration. These modestly energized polar wind ions are then carried to the magnetotail where they are substantially energized by the curvature drift-induced movement across the cross-tail potential of the magnetotail. This latter drift energizes the ions to the energies typically found in the plasma sheet. Subsequent drift and energization can cause the ions to become part of the ring current. The results of this study were presented at the Spring AGU meeting in 2002 and the GEM meeting in June 2003. They are the foundation for a paper that has been submitted by Matthew Huddleston to the Journal of Geophysical Research in December 2003. This work was part of the thesis that Matthew completed in finishing his Ph.D. in Physics at Vanderbilt University.

Description: An improved specification of the plasma environment has been developed for use in modeling spacecraft charging. It was developed by statistically analyzing a large part of the LANL Magnetospheric Plasma Analyzer (MPA) data set for ion and electron spectral signature correlation with spacecraft charging, including anisotropies. The objective is to identify a relatively simple characterization of the full particle distributions that yield an accurate predication of the observed charging under a wide variety of conditions.

Description: There is abundant evidence for the existence of 146Sm in the early solar system and for preservation of effects in the 146Sm-142Nd system in differentiated meteorites ([1]; see recent discussion in Stewart et al. 1994). Information from the 182Hf-182W system, as revised by new careful work [2-3] also indicates that the Earth s core formed relatively early. It is in principle possible for early-formed crust and mantle reservoirs on Earth to have preserved evidence for 146Sm if such reservoirs were produced with high Sm/Nd fractionation and if they have remained isolated and closed since ~4.3 Ga. The mean life of 146Sm of 149 Ma is sufficiently long to make this an intriguing possibility.

Description: Fluid inclusions are microscopic volumes of fluid trapped within minerals as they precipitate. Fluid inclusions are common in terrestrial minerals formed under a wide array of geological settings from surface evaporite deposits to kimberlite pipes. While fluid inclusions in terrestrial rocks are the rule rather than the exception, only few fluid inclusion-bearing meteorites have been documented. The rarity of fluid inclusions in meteoritic material may be explained in two ways. First, it may reflect the absence of fluids (water?) on meteorite parent bodies. Alternatively, fluids may have been present when the rock formed, but any fluid inclusions originally trapped on the parent body were destroyed by the extreme P-T conditions meteorites often experience during impact events. Distinguishing between these two possibilities can provide significant constraints on the likelihood of life on the parent body. Just as textures, structures, and compositions of mineral phases can be significantly altered by shock metamorphism upon hypervelocity impact, fluid inclusions contained within component minerals may be altered or destroyed due to the high pressures, temperatures, and strain rates associated with impact events. Reequilibration may occur when external pressure-temperature conditions differ significantly from internal fluid isochoric conditions, and result in changes in fluid inclusion properties and/or textures. Shock metamorphism and fluid inclusion reequilibration can affect both the impacted target material and the meteoritic projectile. By examining the effects of shock deformation on fluid inclusion properties and textures we may be able to better constrain the pressure-temperature path experienced by shocked materials and also gain a clearer understanding of why fluid inclusions are rarely found in meteoritic samples.

Description: The idea of impact induced volcanism continues to blossom ([1-3] and other references). However, this appealing idea is seldom supported with an appropriate physical mechanism. The aim of this publication is to critically examine some frequently cited mechanisms of impact energy transformation into a trigger for terrestrial volcanism and magmatism.

Description: Based on the lunar cratering record, impacts were larger and more frequent on the early Earth than they are today. There is no persevered record of these early terrestrial impacts because rocks of this age have been obliterated by tectonism and erosion. The oldest known evidence of impacts on Earth lies in four beds (S1, S2, S3 and S4) in the Barberton Greenstone Belt (BGB), South Africa, ranging in age from about 3.24 to 3.47 Ga. These beds are composed in large part of sand-sized spherical particles, termed spherules, that are thought to have formed by the condensation of rock vapor clouds ejected above the atmosphere as a result of large impacts. Spherule beds S2 and S3 are both about 20 cm thick where composed entirely of fall-deposited spherules and up to a meter thick where spherules are mixed with locally derived debris. The diameters the bolides have been estimated to be between 20 and 50 km, based on bed thickness, size of the largest spherules, Ir fluence and extraterrestrial Cr.

Description: Regional to continental scale magnetic anomaly maps are becoming increasingly available from airborne, shipborne, and terrestrial surveys. Satellite data are commonly considered to fill the coverage gaps in regional compilations of these near-surface surveys. For the near-surface Antarctic magnetic anomaly map being produced by the Antarctic Digital Magnetic Anomaly Project (ADMAP), we show that near-surface magnetic anomaly estimation is greatly enhanced by the joint inversion of the near-surface data with the satellite observations relative to the conventional technique such as minimum curvature. Orsted observations are especially advantageous relative to the Magsat data that have order-of-magnitude greater measurement errors, albeit at much lower orbital altitudes. CHAMP is observing the geomagnetic field with the same measurement accuracy as the Orsted mission, but at the lower orbital altitudes covered by Magsat. Hence, additional significant improvement in predicting near-surface magnetic anomalies can result as these CHAMP data are available. Our analysis also suggests that considerable new insights on the magnetic properties of the lithosphere may be revealed by a further order-of-magnitude improvement in the accuracy of the magnetometer measurements at minimum orbital altitude.

Description: Seismic refraction studies have provided critical, but spatially restricted constraints on the structure of the Icelandic crust. To obtain a more comprehensive regional view of this tectonically complicated area, we spectrally correlated free-air gravity anomalies against computed gravity effects of the terrain for a crustal thickness model that also conforms to regional seismic and thermal constraints. Our regional crustal thickness estimates suggest thickened crust extends up to 500 km on either side of the Greenland-Scotland Ridge with the Iceland-Faeroe Ridge crust being less extended and on average 3-5 km thinner than the crust of the Greenland-Iceland Ridge. Crustal thickness estimates for Iceland range from 25-35 km in conformity with seismic predictions of a cooler, thicker crust. However, the deepening of our gravity-inferred Moho relative to seismic estimates at the thermal plume and rift zones of Iceland suggests partial melting. The amount of partial melting may range from about 8% beneath the rift zones to perhaps 20% above the plume core where mantle temperatures may be 200-400 C above normal. Beneath Iceland, areally limited regions of partial melting may also be compositionally and mechanically layered

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: 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: 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: 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: The discrete aurora results when energized electrons bombard the Earth's atmosphere at high latitudes. This paper examines the physical processes that can cause field-aligned acceleration of plasma particles in the auroral region. A data and theoretical study has been carried out to examine the acceleration mechanisms that operate in the auroral zone and to identi@ the magnetospheric drivers of these acceleration mechanisms. The observations used in the study were collected by the Fast Auroral Snapshot (FAST) and Polar satellites when the two satellites were in approximate magnetic conjunction in the auroral region. During these events FAST was in the middle of the auroral zone and Polar was above the auroral zone in the near-Earth plasma sheet. Polar data were used to determine the conditions in the magnetotail at the time field-aligned acceleration was measured by FAST in the auroral zone. For each of the magnetotail drivers identified in the data study, the physics of field-aligned acceleration in the auroral region was examined using existing theoretical efforts and/or a long-system particle in cell simulation to model the magnetically connected region between the two satellites. Results from the study indicate that there are three main drivers of auroral acceleration: (1) field-aligned currents that lead to quasistatic parallel potential drops (parallel electric fields), (2) earthward flow of high-energy plasma beams from the magnetotail into the auroral zone that lead to quasistatic parallel potential drops, and (3) large-amplitude Alfven waves that propagate into the auroral region from the magnetotail. The events examined thus far confm the previously established invariant latitudinal dependence of the drivers and show a strong dependence on magnetic activity. Alfven waves tend to occur primarily at the poleward edge of the auroral region during more magnetically active times and are correlated with intense electron precipitation. At lower latitudes away from the poleward edge of the auroral zone is the primary field-aligned current region which results in the classical field- aligned acceleration associated with the auroral zone (electrons earthward and ion beams tailward). During times of high magnetic activity, high-energy ion beams originating from the magnetotail are observed within, and overlapping, the regions of primary and return field-aligned current. Along the field lines where the high-energy magnetotail ion beams are located, field-aligned acceleration can occur in the auroral zone leading to precipitating electrons and upwelling ionospheric ion beams. Field-aligned currents are present during both quiet and active times, while the Alfven waves and magnetotail ion beams were observed only during more magnetically active events.

Description: We report measurement of whistler-mode chorus by the four Cluster spacecraft at close separations. We focus our analysis on the generation region close to the magnetic equatorial plane at a radial distance of 4.4 Earth's radii. We use both linear and rank correlation analysis to define perpendicular dimensions of the sources of chorus elements below one half of the electron cyclotron frequency. Correlation is significant throughout the range of separation distances of 60-260 km parallel to the field line and 7-100 km in the perpendicular plane. At these scales, the correlation coefficient is independent for parallel separations, and decreases with perpendicular separation. The observations are consistent with a statistical model of the source region assuming individual sources as gaussian peaks of radiated power with a common half-width of 35 km perpendicular to the magnetic field. This characteristic scale is comparable to the wavelength of observed waves.

Description: The radio plasma imager (RPI) on the IMAGE satellite performs radio sounding in the magnetosphere, transmitting coded signals stepping through the frequency range of interest and receiving the returned echoes. The measurements provide the echo amplitude as a function of frequency and echo delay time on a so-called plasmagram. A newly developed algorithm inverts THE echo traces on a plasmagram to electron density spatial distributions. Rased on these observed density distributions, an empirical model is constructed to describe the two-dimensional density distribution in the inner magnetosphere.

Description: Dartmouth College provided a single-channel high-frequency wave receiver to the Sounding of the Ion Energization Region: Resolving Ambiguities (SIERRA) rocket experiment launched from Poker Flat, Alaska, in January 2002. The receiver used signals from booms, probes, preamplifiers, and differential amplifiers provided by Cornell University coinvestigators. Output was to a dedicated 5 MHz telemetry link provided by WFF, with a small amount of additional Pulse Code Modulation (PCM) telemetry required for the receiver gain information. We also performed preliminary analysis of the data. The work completed is outlined below, in chronological order.

Description: The primary goal of the work performed under this grant was to locate, if possible, the onset of subphotospheric convection zones in normal main sequence stars by using the presence of emission in high temperature lines in far ultraviolet spectra from the FUSE spacecraft as a proxy for convection. The change in stellar structure represented by this boundary between radiative and convective stars has always been difficult to find by other empirical means. A search was conducted through observations of a sample of A-type stars, which were somewhat hotter and more massive than the Sun, and which were carefully chosen to bridge the theoretically expected radiative/convective boundary line along the main sequence.

Description: We present X-ray Luminosity Distributions (XLDs) of late-type members (dF, dG, dK, dM) of the Blanco 1 cluster, based on ROSAT-HRI data and new astrometric-photometric membership obtained from the GSC-II project. For the first time we present the XLD of dM stars of this cluster. The high metallicity of Blanco 1 allows us to investigate the role of chemical composition on the coronal emission of late-type stars. Comparison between X-ray Luminosity Distributions of Blanco 1 and Pleiades, NGC2516 and alpha Per suggests a possible metallicity effect in dM stars.

Description: Global pressure data are often needed for processing or interpreting modern geodetic and oceanographic measurements. The most common source of these data is the analysis or reanalysis products of various meteorological centers. Tidal signals in these products can be problematic for several reasons, including potentially aliased sampling of the semidiurnal solar tide as well as the presence of various modeling or timing errors. Building on the work of Van den Dool and colleagues, we lay out a strategy for handling atmospheric tides in (re)analysis data. The procedure also offers a method to account for ocean loading corrections in satellite altimeter data that are consistent with standard ocean-tide corrections. The proposed strategy has immediate application to the on-going Jason-1 and GRACE satellite missions.

Description: We give measurement results of electron temperature and electron density from the Floating Potential Probe (FPP) on the International Space Station (ISS), and relate them to the electron current collection of the ISS solar arrays and the degree of charging of ISS when its Plasma Contacting Units (PCUs) are not operating. We show that on days of high solar activity index Kp, high levels of ISS charging are significantly more probable than on days of low solar activity, due to some abnormally low morning electron temperatures. Although the FPP electron temperatures measured are almost always higher than predicted by the International Reference Ionosphere 90 model (IRI-90), it is shown that the CHAMP satellite Langmuir Probe (PLP) also shows low dawn electron temperatures on the same day as those found by FPP. It is further shown that similar high levels of predicted charging, accompanied by vxB charging on the ISS structure, could exceed the -40 V specification on ISS charging, and could be dangerous to ISS astronauts if the PCUs fail to operate.

Description: The dependency of non-sphericity on gravitational settling of mineral dust particles is parameterized for prolate ellipsoids and Reynolds number lower than 2. The settling speed is numerically solved from the momentum equation as a function of particle diameter and aspect ratio. The reduction of settling speed due to non-sphericity is included in the GO-CART model to simulate dust size distribution for April 2001. Two numerical schemes for solving sedimentation are compared. For particles of diameter greater than 5 micron the simulated size distribution is sensitive to the numerical sedimentation scheme. Changing the particle shape from spherical to non-spherical with lambda=2, makes little difference to the simulated surface concentration and size distribution except at the periphery of the dust sources. However, when very elongated particles (lambda=5) are simulated the differences between non-spherical and spherical particles are significant. With limited in-situ measurements reporting most frequent lambda around 1.5, the overall effects on global modeling is rather negligible and the essential benefit is to relax the CFL condition of Eulerian settling schemes.

Description: NASA s vision for Earth science is to build a "sensor web": an adaptive array of heterogeneous satellites and other sensors that will track important events, such as storms, and provide real-time information about the state of the Earth to a wide variety of customers. Achieving this vision will require automation not only in the scheduling of the observations but also in the processing of the resulting data. To address this need, we are developing a planner-based agent to automatically generate and execute data-flow programs to produce the requested data products.

Description: A central goal of extragalactic observational astronomy is to understand how normal galaxies evolve with redshift, and particularly when galaxies formed their stars. While optical and rest-frame UV observations have begun to address these issues, the interpretation of such data is particularly challenging because of the sensitivity to dust obscuration (at optical and UV wavelengths). The absorbed light is re-radiated at IR wavelengths, hence the optimal indicators of the star formation rate (SFR) is at a rest-frame wavelength of approx. 60 microns. The SIRTF mission will revolutionize the study of the global evolution of the SFR by providing mass-selected, complete samples of galaxies and far-IR estimators of the SFR. SIRTF will be two orders of magnitude more sensitive than previous IR space missions, and therefore will open up a completely new region of parameter space in sensitivity and angular resolution. This research program will study the SFR using statistical samples of galaxies in the local universe, at intermediate redshifts, and set the stage for continuing studies up to z = 5. The overall research program is divided into three main investigations: a Mid-IR Hubble Atlas and SFR estimators in the local universe, evolution of the SFR at 0 less than z less that 1 using pencil beam redshift surveys, and Galaxy formation and evolution at 1 less than z less than 5.

Description: Overview of project: we combined and jointly analysed geophysical, remote sensing and glaciological data for investigating the temporal changes in ice flow and the role of geologic control on glacial drainage. The project included two different studies, the investigation of recent changes of the Kangerlussuaq glacier and the study of geologic control of ice flow in NW Greenland, around the Humboldt, Petermann and Ryder glaciers.

Description: Ion charge states measured in situ in interplanetary space are formed in the inner coronal regions below 5 solar radii, hence they carry information on the properties of the solar wind plasma in that region. The plasma parameters that are important in the ion forming processes are the electron density, the electron temperature and the flow speeds of the individual ion species. In addition, if the electron distribution function deviates from a Maxwellian already in the inner corona, then the enhanced tail of that distribution function, also called halo, greatly effects the ion composition. The goal of the proposal is to make use of ion fractions observed in situ in the solar wind to learn about both, the plasma conditions in the inner corona and the expansion and ion formation itself. This study is carried out using solar wind models, coronal observations, and ion fraction calculations in conjunction with the in situ observations.

Description: The ultimate goal of this project is to better understand the current transfer of mass between the Greenland Ice Sheet, the world's oceans and the atmosphere, and to identify processes controlling the rate of this transfer, to be able to predict with greater confidence future contributions to global sea level rise. During the first year of this project, we focused on establishing longer-term records of change of selected outlet glaciers, reevaluation of mass input to the ice sheet and analysis of climate records derived from ice cores, and modeling meltwater production and runoff from the margins of the ice sheet.

Description: The grant has supported the completion of 16 papers and 4 conference proceedings to date. During the first year of the project we completed five papers, each of which represents a new direction in the theory and interpretation of gravitational lenses. In the first paper, "The Importance of Einstein Rings", we developed the first theory for the formation and structure of the Einstein rings formed by lensing extended sources like the host galaxies of quasar and radio sources. We applied the theory to three lenses with lensed host galaxies. For the time delay lens PG 1115+080 we found that the structure of the Einstein ring ruled out models of the gravitational potential which permitted a large Hubble constant (70 km/s Mpc). In the second paper, :Cusped Mass Models Of Gravitational Lenses", we introduced a new class of lens models where the central density is characterized by a cusp ( rho proportional to tau(sup -gamma), 1 less than gamma less than 2) as in most modern models and theories of galaxies rather than a finite core radius. In the third paper, "Global Probes of the Impact of Baryons on Dark Matter Halos", we made the first globally consistent models for the separation distribution of gravitational lenses including both galaxy and cluster lenses. We show that the key physics for the origin of the sharp separation cutoff in the separation distribution near 3 arc sec is the effect of the cooling baryons in galaxies on the density structure of the system.

Description: The Montreal Protocol restricts production of ozone-depleting halocarbons worldwide. Enforcement of the protocol has relied mainly on annual government statistics of production and consumption of these compounds (bottom-up approach). We show here that aircraft observations of ha1ocarbon:CO enhancement ratios on regional to continental scales can be used to infer halocarbon emissions, providing independent verification of the bottom-up approach. We apply this topdown approach to aircraft observations of Asian outflow &om the TRACE-P mission over the western Pacific (March-April 2001) and derive emissions from eastern Asia (China, Japan, and Korea). We derive an eastern Asian carbon tetrachloride (CCl4) source of 21.5 Gg yr(sup -1), several-fold larger than previous estimates and amounting to -30% of the global budget for this gas. Our emission estimate for CFC-11 from eastern Asia is 50% higher than inventories derived from manufacturing records. Our emission estimates for methyl chloroform (CH3CC13) and CFC-12 are in agreement with existing inventories. For halon 1211 we find only a strong local source originating from the Shanghai area. Our emission estimates for the above gases result in a approximately equal to 40% increase in the ozone depletion potential (ODP) of Asian emissions relative to previous estimates, corresponding to a approximately equal to 10% global increase in ODP.

Description: We analyze the Asian outflow of CO, ozone, and nitrogen oxides (NOx) to the Pacific in spring by using the GEOS-CHEM global three-dimensional model of tropospheric chemistry and simulating the Pacific Exploratory Mission-West (PEM-West B) aircraft mission in February-March 1994. The GEOS-CHEM model uses assimilated meteorological fields from the NASA Goddard Earth Observing System (GEOS). It reproduces relatively well the main features of tropospheric ozone, CO, and reactive nitrogen species observed in PEM-West B, including latitudinal and vertical gradients of the Asian pollution outflow over the western Pacific although simulated concentrations of CO tend to be too low (possibly because biogenic sources are underestimated). We use CO as a long-lived tracer to diagnose the processes contributing to the outflow. The highest concentrations in the outflow are in the boundary layer (0-2 km), but the strongest outflow fluxes are in the lower free troposphere (2-5 km) and reflect episodic lifting of pollution over central and eastern China ahead of eastward moving cold fronts. This frontal lifting, followed by westerly transport in the lower free troposphere, is the principal process responsible for export of both anthropogenic and biomass burning pollution from Asia. Anthropogenic emissions from Europe and biomass burning emissions from Africa make also major contributions to the Asian outflow over the western Pacific; European sources dominate in the lower troposphere north of 40 degrees N, while African sources are important in the upper troposphere at low latitudes. For the period of PEM-West B (February-March) we estimate that fossil fuel combustion and biomass burning make comparable contributions to the budgets of CO, ozone, and NO, in the Asian outflow. We find that 13% of NO, emitted in Asia is exported as NO, or PAN, a smaller fraction than for the United States because of higher aerosol concentrations that promote heterogeneous conversion of NOx to HNO3. Production and export of ozone from Asia in spring is much greater than from the United States because of the higher photochemical activity.

Description: This short program of LIDORT-based research in atmospheric trace gas retrieval was conducted over the 1 year period 01 July 2002 to 30 June 2003. After consultation with the NASA reporting officer, the first of the two original proposal activities (development of a direct-fitting total O3 column retrieval algorithm with operational capability for GOME data) was replaced by other tasks. The three activities addressed were: (1) Sensitivity studies for column and profile retrieval of NO2 distributions from a new generation of multi-axis ground-based spectrometers; (2) use of the LIDORT-RRS model to determine the effect of inelastic rotational Raman scattering at SBUV wavelengths; (3) an examination of ozone profile weighting functions in the presence of optically thick tropospheric clouds.

Description: We use measurements and models to develop aerosol models for use in the inversion algorithms for the Cloud Aerosol Lidar and Imager Pathfinder Spaceborne Observations (CALIPSO). Radiance measurements and inversions of the AErosol RObotic NETwork (AERONET1, 2) are used to group global atmospheric aerosols using optical and microphysical parameters. This study uses more than 105 records of radiance measurements, aerosol size distributions, and complex refractive indices to generate the optical properties of the aerosol at more 200 sites worldwide. These properties together with the radiance measurements are then classified using classical clustering methods to group the sites according to the type of aerosol with the greatest frequency of occurrence at each site. Six significant clusters are identified: desert dust, biomass burning, urban industrial pollution, rural background, marine, and dirty pollution. Three of these are used in the CALIPSO aerosol models to characterize desert dust, biomass burning, and polluted continental aerosols. The CALIPSO aerosol model also uses the coarse mode of desert dust and the fine mode of biomass burning to build a polluted dust model. For marine aerosol, the CALIPSO aerosol model uses measurements from the SEAS experiment 3. In addition to categorizing the aerosol types, the cluster analysis provides all the column optical and microphysical properties for each cluster.

Description: This report covers the time period 1 January 2002 to 31 October 2003. During this period we had two balloon flights, continued analyzing data from past and recent flights, exploring issues such as radical partitioning, stratospheric transport, and molecular spectroscopy and further developed our beamsplitter technology.

Description: The project has resulted in four separate investigations, which are each in various stages of publication in the refereed scientific journals. The first investigation was of the generation of electrostatic electron cyclotron waves observed by the Polar spacecraft throughout the auroral regions, dayside cusp, and polar magnetosphere. We have since discovered that these waves are also present within the magnetopause and magnetosheath, which is one of the topics of a second study, entitled: 'Polar observations of plasma waves in and near the dayside magnetopause/magnetosheath.' A third study of plasma waves focussed on kilometric continuum (KC) emission. This work is reported in a paper entitled 'Near-source and Remote Observations of Kilometric Continuum Radiation From Multi-spacecraft Observations'.The final investigation of this program concerns the possible transverse heating of auroral ions by impulsive wave structures. We summarize that substantial transverse ion heating has already occurred at lower altitudes. Abstracts of the above four studies are included in the Appendix to this final report.

Description: This paper focuses on biosphere-climate interactions including the influences of human activities. Recognizing this is only one aspect of biospheric processes, this places an emphasis of those biogeochemical processes that have a profound effect on numerous other aspects of the biosphere and the services it provides, services which are critical to sustaining life on Earth. And, the paper will focus on the various scientific aspects of assessing the availability of fresh water, including its sensitivity to climate variance and land use changes. Finally, this paper hopes to emphasize the potential role that greatly expanded space observations and interactive modeling can play in developing our understanding of Earth and its the living systems.

Description: The Third International Conference on Large Meteorite Impacts (formerly Large Meteorite Impacts and Planetary Evolution) was held August 5-7, 2003, at the Kloesterle, Noerdlingen. The conference addressed a wide range of topics relating to planetary impacts and their effects on planetary crusts, climate, life on Earth, rocks, and rock-forming minerals. This CD-ROM contains the preface, table of contents, program, abstracts and indexes.

Description: The progress on developing models of the radiation environment since the 1960s is reviewed with emphasis on models that can be applied to predicting the performance of microelectronics used in spacecraft and instruments. Space, atmospheric, and ground environments are included. It is shown that models must be adapted continually to account for increased understanding of the dynamics of the radiation environment and the changes in microelectronics technology. The IEEE Nuclear and Space Radiation Effects Conference is a vital forum to report model progress to the radiation effects research community.

Description: The NASA Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft mission was conducted in February-April 2001 over the NW Pacific (1) to characterize the Asian chemical outflow and relate it quantitatively to its sources and (2) to determine its chemical evolution. It used two aircraft, a DC-8 and a P-3B, operating out of Hong Kong and Yokota Air Force Base (near Tokyo), with secondary sites in Hawaii, Wake Island, Guam, Okinawa, and Midway. The aircraft carried instrumentation for measurements of long-lived greenhouse gases, ozone and its precursors, aerosols and their precursors, related species, and chemical tracers. Five chemical transport models (CTMs) were used for chemical forecasting. Customized bottom-up emission inventories for East Asia were generated prior to the mission to support chemical forecasting and to serve as a priori for evaluation with the aircraft data. Validation flights were conducted for the Measurements Of Pollution In The Troposphere (MOPITT) satellite instrument and revealed little bias (6 plus or minus 2%) in the MOPITT measurements of CO columns. A major event of transpacific Asian pollution was characterized through combined analysis of TRACE-P and MOPITT data. The TRACE-P observations showed that cold fronts sweeping across East Asia and the associated warm conveyor belts (WCBs) are the dominant pathway for Asian outflow to the Pacific in spring. The WCBs lift both anthropogenic and biomass burning (SE Asia) effluents to the free troposphere, resulting in complex chemical signatures. The TRACE-P data are in general consistent with a priori emission inventories, lending confidence in our ability to quantify Asian emissions from socioeconomic data and emission factors. However, the residential combustion source in rural China was found to be much larger than the a priori, and there were also unexplained chemical enhancements (HCN, CH3Cl, OCS, alkylnitrates) in Chinese urban plumes. The Asian source of CCl4 was found to be much higher than government estimates. Measurements of HCN and CH3CN indicated a dominant biomass burning source and ocean sink for both gases. Large fractions of sulfate and nitrate were found to be present in dust aerosols. Photochemical activity in the Asian outflow was strongly reduced by aerosol attenuation of UV radiation, with major implications for the concentrations of HOx, radicals. New particle formation, apparently from ternary nucleation involving NH3, was observed in Chinese urban plumes.

Description: A complete description of a self-consistent model of magnetospheric ring current interacting with electromagnetic ion cyclotron waves is presented. The model is based on the system of two kinetic equations; one equation describes the ring current ion dynamics, and another equation describes the wave evolution. The effects on ring current ions interacting with electromagnetic ion cyclotron waves and back on waves are considered self-consistently by solving both equations on a global magnetospheric scale under nonsteady state conditions. The developed model is employed to simulate the entire 2-7 May 1998 storm period. First, the trapped number fluxes of the ring current protons are calculated and presented along with comparison with the data measured by the three- dimensional hot plasma instrument Polar/HYDRA. Incorporating in the model the wave-particle interaction leads to much better agreement between the experimental data and the model results. Second, examining of the wave (MLT, L shell) distributions produced by the model during the storm progress reveals an essential intensification of the wave emission about 2 days after the main phase of the storm. This result is well consistent with the earlier ground-based observations. Finally, the theoretical shapes and the occurrence rates of the wave power spectral densities are studied. It is found that about 2 days after the storm s main phase on 4 May, mainly non-Gaussian shapes of power spectral densities are produced.

Description: The main objective of this project was to better understand the geodynamics controlling the mountain building and topographic evolution in the central Andes using an integrated approach that combines GPS, seismological, and numerical studies.

Description: Tn determine the effect of crustal thickness variation on satellite-altitude geopotential anomalies we compared two regions of Europe with vastly different values, South and Central Finland and the Pannonian Basin. In our study regions, crustal thickness exceeds 44 km in Finland and is less than 26 km in the Pannonian Basin. Heat-flow data indicate that the thinner and more active crust of the Pannonian Basin has a value nearly three times that of the Finnish Svecofennian Province. An ovoid positive CHAMP gravity anomaly (-4 mGal) is quasi-coincidental with the CHAMP magnetic anomaly traverses the Pannonian Basin while ground based gravity mapping in Hungary shows that the free-air gravity anomalies across the Pannonian Basin are near 0 to +20 mGal with shorter wavelength anomalies from +40 to less than +60 mGal and some 0 to greater than -20 mGal. Larger anomalies are detected in the mountainous areas. The minor value anomalies can indicate the isostatic equilibrium for Hungary (the central part of the Pannonian Basin). Gravity data over Finland bear overprint of de-glaciation. CHAMP gravity data indicates a west-east positive gradient of less than 4 mGal across South and Central Finland. CHAMP magnetic data (400 km) reveal elongated semi-circular negative anomalies for both regions with South-Central Finland having larger amplitude (less than -6 nT) than that over the Pannonian Basin, Hungary (less than -5 nT). In the latter subducted oceanic lithosphere has been proposed as the anomalous body.

Description: CHAMP is recording state-of-the-art magnetic and gravity field observations at altitudes ranging over roughly 300 - 550 km. However, anomaly continuation is severely limited by the non-uniqueness of the process and satellite anomaly errors. Indeed, our numerical anomaly simulations from satellite to airborne altitudes show that effective downward continuations of the CHAMP data are restricted to within approximately 50 km of the observation altitudes while upward continuations can be effective over a somewhat larger altitude range. The great unreliability of downward continuation requires that the satellite geopotential observations must be analyzed at satellite altitudes if the anomaly details are to be exploited most fully. Given current anomaly error levels, joint inversion of satellite and near- surface anomalies is the best approach for implementing satellite geopotential observations for subsurface studies. We demonstrate the power of this approach using a crustal model constrained by joint inversions of near-surface and satellite magnetic and gravity observations for Maude Rise, Antarctica, in the southwestern Indian Ocean. Our modeling suggests that the dominant satellite altitude magnetic anomalies are produced by crustal thickness variations and remanent magnetization of the normal polarity Cretaceous Quiet Zone.

Description: By modifying the IGRF it is possible to learn what may happen to the geomagnetic field during a geomagnetic reversal. If the entire IGRF reverses then the declination and inclination only reverse when the field strength is zero. If only the dipole component of the IGRF reverses a large geomagnetic field remains when the dipole component is zero and he direction of the field at the end of the reversal is not exactly reversed from the directions at the beginning of the reversal.

Description: Observations and analyses of mesoscale gravity waves in the stratosphere from the Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) are summarized, with focus on global distribution of topography related wave activities. We found most of the orographical wave activities occur during the winter seasons over high latitude mountain ridges. In the northern hemisphere, the strongest waves are those over Scandinavia, Central Eurasia, and southern Greenland, whereas in the southern hemisphere, wave activities are outstanding over the Andes, New Zealand, and Antarctic rim;, MLS observations suggest that these orographic waves are located mostly on the down stream side of the mountain ridge with downward phase progression and have horizontal phase velocities opposite to the stratospheric jet-stream. Future studies using MLS data and numerical modeling will lead to better understanding of gravity wave effects on dynamics and chemistry in the middle atmosphere.

Description: The Galileo spacecraft Energetic Particle Detector (EPD) has been used to map Jupiter's trapped electron radiation in the jovian equatorial plane for the range 8 to 16 Jupiter radii.

Description: The goal of this study is in improvement of the short-term prediction of the polar motion in the context of state-space modeling suitable for the Kalman filter operations, such as the Kalman Earth Orientation Filter (KEOF) at Jet Propulsion Laboratory.

Description: The vast illuminating power of the Global Positioning System (GPS), which transformed space geodesy in the 199Os, is now serving to probe the earth's fluid envelope in unique ways. Three distinct techniques have emerged: ground-based sensing of the integrated atmospheric moisture; space-based profiling of atmospheric refractivity, pressure, temperature, moisture, and other properties by active limb sounding; and surface (ocean and ice) altimetry and scatterometry with reflected signals detected from space. Ground-based GPS moisture sensing is already in provisional use for numerical weather prediction. Limb sounding, while less mature, offers a bevy of attractions, including high accuracy, stability, and vertical resolution; all-weather operation; and exceptionally low cost. GPS bistatic radar, r 'reflectometry,' is the least advanced but shows promise for a number of niche applications.

Description: The Jet Propulsion Laboratory (JPL) Microwave Temperature Profiler (MTP) was the first United States instrument to fly on the Russian Geophysica high-altitude research aircraft. Careful comparison of MTP measurements with radiosondes launched near the Geophysica flight track has allowed us to establish the flight level temperature to an accuracy of 0.2K.

Description: In the interest of improving the performance and efficiency of space geodesy a diverse group in the U.S., in collaboration with IGGOS, has begun to establish a unified National Geodetic Observatory (NGO).

Description: Hydrogen emissions are the signature of proton aurora. The Doppler-shifted hydrogen emission lines can be interpreted in terms of the mean energy of the precipitating protons. A red shifted component of the line profiles observed from the ground indicates upward going hydrogen atoms due to angular redistribution of the precipitation. Secondary electrons from ionization and stripping collisions also contribute to the auroral emissions. Since the energy distribution of these secondaries has a lower mean energy than secondary electrons in electron aurora, the relative brightness of eniission features differs from that in electron aurora. The secondaries contribute little to additional ionization. These differences between proton and electron aurora can lead to misinterpretation when brightness ratios are used to derive ionospheric conductances with parameterizations that are based on electron aurora.

Description: Recent drilling has produced the samples of impact melt breccias examined here. All samples examined here are from the YAX-1 borehole [7]. Two main types of impact breccias have been studied. The first type is a green altered impact meltrock found in the lower portion of the impact sequence. The texture of the rock is microcrystalline and is composed of pyroxene, plagioclase, and alkali feldspars. Its composition is consistent with continental margin rocks. It is generally massive with some flow structure. The rock was brecciated and altered after solidfication and contains small amounts of both shocked and unshocked clasts of the impacted lithologies. These lithologies include lithic quartzite, and isolated feldspar crystals. The compositions of these rocks are similar to those seen in meltrocks sampled by the Yucatan-6 borehole [8-9]. Our study includes samples YAX-1_861.4, YAX-1_863.51, and YAX- 1_876.46, which represent both the top and lower portion of the green impact meltrock. The middle sample in the sequence has the least amount of (mineralogical) alteration [10]. The second type of melt breccia under study is a brown altered impact meltrock. It also has a microcrystalline texture and both shocked and unshocked clasts of the target material. Even though this rock type has been altered, remnant schleiren, metaquartzite, and micritic calcite have been identified. Sample YAX-1_841.32 is representative of this type of rock. It was recovered from a polymict breccia in the middle of the impact sequence.

Description: The Chicxulub Scientific Drilling Project (www.icdp-online.de) recovered a continuous core from a depth of 404 m (in Tertiary cover) to 1511 m (in a megablock of Cretaceous target sediments), penetrating approx. 100 m of melt-bearing impactites between 794 and 895 m. The Yaxcopoil-1 (YAX-1) borehole is approx. 60-65 km from the center of the Chicxulub structure, which is approx. 15 km beyond the limit of the estimated approx. 50 km radius transient crater (excavation cavity), but within the rim of the estimated approx. 90 km radius final crater. In general, the impactite sequence is incredibly rich in impact melts of unusual textural variety and complexity, quite unlike melt-bearing impact formations from other terrestrial craters.

Description: The Silverpit crater is a recently discovered, 60-65 Myr old complex crater, which lies buried beneath the North Sea, about 150 km east of Britain. High-resolution images of Silverpit's subsurface structure, provided by three-dimensional seismic reflection data, reveal an inner-crater morphology similar to that expected for a 5-8 km diameter terrestrial crater. The crater walls show evidence of terrace-style slumping and there is a distinct central uplift, which may have produced a central peak in the pristine crater morphology. However, Silverpit is not a typical 5-km diameter terrestrial crater, because it exhibits multiple, concentric rings outside the main cavity. External concentric rings are normally associated with much larger impact structures, for example Chicxulub on Earth, or Orientale on the Moon. Furthermore, external rings associated with large impacts on the terrestrial planets and moons are widely-spaced, predominantly inwardly-facing, asymmetric scarps. However, the seismic data show that the external rings at Silverpit represent closely-spaced, concentric faultbound graben, with both inwardly and outwardly facing fault-scarps. This type of multi-ring structure directly analogous to the Valhalla-type multi-ring basins found on the icy satellites. Thus, the presence and style of the multiple rings at Silverpit is surprising given both the size of the crater and its planetary setting. A further curiosity of the Silverpit structure is that the external concentric rings appear to be extensional features on the West side of the crater and compressional features on the East side. The crater also lies in a local depression, thought to be created by postimpact movement of a salt layer buried beneath the crater.

Description: By combining complementary monitoring observations spanning long, medium and short time scales, we have constructed power spectral densities (PSDs) of six Seyfert 1 galaxies. These PSDs span approx. greater than 4 orders of magnitude in temporal frequency, sampling variations on time scales ranging from tens of minutes to over a year. In at least four cases, the PSD shows a "break," a significant departure from a power law, typically on time scales of order a few days. This is similar to the behavior of Galactic X-ray binaries (XRBs), lower mass compact systems with breaks on time scales of seconds. NGC 3783 shows tentative evidence for a doubly-broken power law, a feature that until now has only been seen in the (much better-defined) PSDs of low-state XRBs. It is also interesting that (when one previously-observed object is added to make a small sample of seven), an apparently significant correlation is seen between the break time scale T and the putative black hole mass M(sub BH), while none is seen between break time scale and luminosity. The data are consistent with the linear relation T = M(sub BH) /10(exp 6.5) solar mass; extrapolation over 6-7 orders of magnitude is in reasonable agreement with XRBs. All of this strengthens the case for a physical similarity between Seyfert 1s and XRBs.

Description: The Polar Cap is the Upper-Atmosphere cum Mag-netosphere region which is enclosed by the poleward boundary of the Auroral Oval and is threaded by open geomagnetic tield lines. In this region, there is normally a steady precipition (Polar "drizzle") of low energy (w 300eV) electrons that excite optical emissions from the ionosphere. At times, enhanced ionization patches are formed near the Dayside Cusp regions that drift across the Polar Cap towards the Night Sector of the Auroral Oval. Discrete auroral arcs and auroras formed during Solar Magnetic Cloud (SMC)/Coronal Mass Ejection (CME) events are also observed in the Polar Cap. Spectrophotometric observations of all these Polar Cap phenomena provide a measure of the average energy as well a energy flux of the electrons precipitating in the Polar Cap region during these disturbances. Such measurements also point to modulations of the Polar Cap Mesosphere-Lower Thermosphere (MLT) air density and temperature by zonally symmetric tides whose Hough functions peak in the Polar region. MLT cooling during Stratospheric Warming events and their relation to Polar Vortex and associated Gravity wave activities are also observed at the Polar Cap sites.

Description: For over a decade, the oldest evidence for life on this planet has been microfossils in the 3.5 Ga Apex Chert in Western Australia. Recently, the biogenicity of these carbon-rich structures has been called into question through reanalysis of the local geology and reinterpretation of the original thin sections. Although initially described as a stratiform, bedded chert of siliceous clasts, the unit is now thought to be a brecciated hydrothermal vein chert. The high temperatures of a hydrothermal environment would probably have detrimental effects to early non-hyperthermophilic life, compared to that of a shallow sea. Conversely, a hydrothermal origin would suggest that if the microfossils were valid, they might have been hyperthermophilic. Apex Chert controversy. The Apex Chert microfossils were originally described as septate filaments composed of kerogen similar in morphology to Proterozoic and modern cyanobacteria. However new thin section analysis shows that these carbonaceous structures are not simple filaments. Many of the original microfossils are branched and have variable thickness when the plane of focus is changed. Hydrothermal alteration of organic remains has also been suggested for the creation of these strange morphologies. Another point of contention lies with the nature of the carbon material in these proposed microfossils. Kerogen is structurally amorphous, but transforms into well-ordered graphite under high pressures and temperatures. Raman spectrometry of the carbonaceous material in the proposed microfossils has been interpreted both as partially graphitized kerogen and amorphous graphite. However, these results are inconclusive, since Raman spectrometry cannot adequately discriminate between kerogen and disordered graphite. There are also opposing views for the origin of the carbon in the Apex Chert. The carbon would be biogenic if the proposed microfossils are indeed the remains of former living organisms. However, an inorganic Fischer- Tropsch-type synthesis is also a possible explanation for the formation of large-aggregate carbonaceous particles and could also account for the depletion of (13)C observed.

Description: Experiments investigating the partitioning of V, Cr, and Si between metal and silicate at various pressures, temperatures, redox state, and composition demonstrate that V and Cr are always more siderophile than Si. The relatively high abundances of V and Cr in the Earth's upper mantle indicate that the high Mg/Si ratio of the Earth's upper mantle cannot be attributed to extraction of Si into the core and must be an intrinsic bulk property of the silicate Earth.

Description: Interferometric synthetic aperture radar (InSAR) provides a practical means of mapping creep along major strike-slip faults. The small amplitude of the creep signal (less than 10 mm/yr), combined with its short wavelength, makes it difficult to extract from long time span interferograms, especially in agricultural or heavily vegetated areas. We utilize two approaches to extract the fault creep signal from 37 ERS SAR images along the southem San Andreas Fault. First, amplitude stacking is utilized to identify permanent scatterers, which are then used to weight the interferogram prior to spatial filtering. This weighting improves correlation and also provides a mask for poorly correlated areas. Second, the unwrapped phase is stacked to reduce tropospheric and other short-wavelength noise. This combined processing enables us to recover the near-field (approximately 200 m) slip signal across the fault due to shallow creep. Displacement maps fiom 60 interferograms reveal a diffuse secular strain buildup, punctuated by localized interseismic creep of 4-6 mm/yr line of sight (LOS, 12-18 mm/yr horizontal). With the exception of Durmid Hill, this entire segment of the southern San Andreas experienced right-lateral triggered slip of up to 10 cm during the 3.5-year period spanning the 1992 Landers earthquake. The deformation change following the 1999 Hector Mine earthquake was much smaller (4 cm) and broader than for the Landers event. Profiles across the fault during the interseismic phase show peak-to-trough amplitude ranging from 15 to 25 mm/yr (horizontal component) and the minimum misfit models show a range of creeping/locking depth values that fit the data.

Description: A 30-year satellite record of sea ice extents derived mostly from satellite microwave radiometer observations reveals that the Arctic sea ice extent decreased by 0.30+0.03 x 10(exp 6) square kilometers per 10 yr from 1972 through 2002, but by 0.36 plus or minus 0.05 x 10(exp 6) square kilometers per 10yr from 1979 through 2002, indicating an acceleration of 20% in the rate of decrease. In contrast, the Antarctic sea ice extent decreased dramatically over the period 1973-1977, then gradually increased. Over the full 30-year period, the Antarctic ice extent decreased by 0.15 plus or minus 0.08 x 10(exp 6) square kilometers per 10 yr. The trend reversal is attributed to a large positive anomaly in Antarctic sea ice extent in the early 1970's, an anomaly that apparently began in the late 1960's, as observed in early visible and infrared satellite images.

Description: The 50,000 year old, 1.8km diameter Lonar crater is located in Maharashtra, India. This relatively small crater is of particular interest because of its unique morphological and mineralogical properties, which make it a valid analogue for similar craters on the surface of Mars. We show that even in this relatively small crater, substantial hydrothermal alteration of shocked breccias in the floor of the crater has occurred, probably due to the thermal effects of the impact event. The 38 my old, 28 km diameter, Mistastin crater contains an 80 m thick impact melt sheet. We have also documented the presence of alteration phases in the material from this larger crater.

Description: This 12-part report documents the data obtained from various sensor measurements taken aboard the Russian Express-A2 and Express-A3 spacecraft in Geosynchronous Earth Orbit (GEO). These GEO communications satellites, which were designed and built by NPO Prikladnoy Mekhaniki (NPO PM) of Zheleznogorsk, Russia, utilize Hall thruster propulsion systems for north-south and east-west stationkeeping and as of June 2002, were still operating at 80 deg. E. and 11 deg. W respectively. Express-A2 was launched on March 12, 2000, while Express-A3 was launched on June 24, 2000. The diagnostic equipment from which these data were taken includes electric field strength sensors, ion current and energy sensors, and pressure sensors. The diagnostics and the Hall thruster propulsion systems are described in detail along with lists of tabular data from those diagnostics and propulsion system and other satellite systems. Space Power, Inc., now part of Pratt & Whitney's Chemical Systems Division, under contract NAS3 99151 to the NASA Glenn Research Center, obtained these data over several periods from March 12, 2000, through September 30, 2001. Each of the 12 individual reports describe, in detail, the propulsion systems as well as the diagnostic sensors utilized. Finally, parts 11 and 12 include the requirements to which NPO PM prepared and delivered these data.

Description: The study of extraterrestrial matter in sediments began with the discovery of cosmic spherules during the HMS Challenger Expedition (1873-1876), but has evolved into a multidisciplinary study of the chemical, physical, and isotopic study of sediments. Extraterrestrial matter in sediments comes mainly from dust and large impactors from the asteroid belt and comets. What we know of the nature of these source materials comes from the study of stratospheric dust particles, cosmic spherules, micrometeorites, meteorites, and astronomical observations. The most common chemical tracers of extraterrestrial matter in sediments are the siderophile elements, most commonly iridium and other platinum group elements. Physical tracers include cosmic and impact spherules, Ni-rich spinels, meteorites, fossil meteorites, and ocean-impact melt debris. Three types of isotopic systems have been used to trace extraterrestrial matter. Osmium isotopes cannot distinguish chondritic from mantle sources, but provide a useful tool in modeling long-term accretion rates. Helium isotopes can be used to trace the long-term flux of the fine fraction of the interplanetary dust complex. Chromium isotopes can provide unequivocal evidence of an extraterrestrial source for sediments with high concentrations of meteoritic Cr. The terrestrial history of impacts, as recorded in sediments, is still poorly understood. Helium isotopes, multiple Ir anomalies, spherule beds, and craters all indicate a comet shower in the late Eocene. The Cretaceous-Tertiary boundary impact event appears to have been caused by a single carbonaceous chondrite projectile, most likely of asteroid origin. Little is known of the impact record in sediments from the rest of the Phanerozoic. Several impact deposits are known in the Precambrian, including several possible mega-impacts in the Early Archean.

Description: The Grant supported research addressing the question of how the NASA Solar Terrestrial Probes (STP) Mission called Geospace electrodynamics Connections (GEC) will resolve space-time structures as well as collect sufficient information to solve the coupled thermosphere-ionosphere- magnetosphere dynamics and electrodynamics. The approach adopted was to develop a high resolution in both space and time model of the ionosphere-thermosphere (I-T) over altitudes relevant to GEC, especially the deep-dipping phase. This I-T model was driven by a high- resolution model of magnetospheric-ionospheric (M-I) coupling electrodynamics. Such a model contains all the key parameters to be measured by GEC instrumentation, which in turn are the required parameters to resolve present-day problems in describing the energy and momentum coupling between the ionosphere-magnetosphere and ionosphere-thermosphere. This model database has been successfully created for one geophysical condition; winter, solar maximum with disturbed geophysical conditions, specifically a substorm. Using this data set, visualizations (movies) were created to contrast dynamics of the different measurable parameters. Specifically, the rapidly varying magnetospheric E and auroral electron precipitation versus the slower varying ionospheric F-region electron density, but rapidly responding E-region density.

Description: The grant funds were provided to carry out chemical tracer wind and turbulence measurements in the mesosphere and lower thermosphere as part of the Turbulent Oxygen Mixing Experiment (TOMEX) for which Dr. James Hecht of the Aerospace Corp. was the Principal Investigator. Clemson University designed, built, and tested two chemical tracer release payload sections for the 21.126 and 21.127 payloads which, in addition, had photometer and ionization gauge instrumentation. The tracer chemical was trimethyl aluminum (TMA). The experiment was carried out at White Sands Missile Range in New Mexico on October 26,2000. The location was chosen because of the proximity to the Starfire Optical Range (SOR) which hosted the powerful University of Illinois sodium lidar for an extended period prior to and also during the launch window. Since the SOR telescope is fully steerable, lidar measurements in the same volume sampled by the rocket were possible. The primary objective of the experiment was to measure the turbulent diffusion and mixing in the mesosphere and lower thermosphere, especially in layers characterized by convective and/or dynamical instabilities. The lidar, which is capable of measuring the sodium density, temperatures, and winds with good range and time resolution, provided the launch criteria, as well as context measurements for the in situ rocket observations.

Description: We propose a scheme to model the absorbing aerosol index and improve the biomass carbon inventories by optimizing the difference between TOMS aerosol index (AI) and modeled AI with an inverse model. Two absorbing aerosol types are considered, including biomass carbon and mineral dust. A priori biomass carbon source was generated by Liousse et al [1996]. Mineral dust emission is parameterized according to surface wind and soil moisture using the method developed by Ginoux [2000]. In this initial study, the coupled CCM1 and GRANTOUR model was used to determine the aerosol spatial and temporal distribution. With modeled aerosol concentrations and optical properties, we calculate the radiance at the top of the atmosphere at 340 nm and 380 nm with a radiative transfer model. The contrast of radiance at these two wavelengths will be used to calculate AI. Then we compare the modeled AI with TOMS AI. This paper reports our initial modeling for AI and its comparison with TOMS Nimbus 7 AI. For our follow-on project we will model the global AI with aerosol spatial and temporal distribution recomputed from the IMPACT model and DAO GEOS-1 meteorology fields. Then we will build an inverse model, which applies a Bayesian inverse technique to optimize the agreement of between model and observational data. The inverse model will tune the biomass burning source strength to reduce the difference between modelled AI and TOMS AI. Further simulations with a posteriori biomass carbon sources from the inverse model will be carried out. Results will be compared to available observations such as surface concentration and aerosol optical depth.

Description: Seismic studies suggest that the bulk of the mantle is heterogeneous, with density variations in depth as well as in horizontal directions (latitude and longitude). This density variation produces a three- dimensional gravity field throughout the Earth. On the other hand, the core density also varies in both time and space, due to convective core flow. Consequently, the fluid outer core and the solid mantle interact gravitationally due to the mass anomalies in both regions. This gravitational core-mantle interaction could play a significant role in exchange of angular momentum between the core and the mantle, and thus the change in Earth's rotation on time scales of decades and longer. Aiming at estimating the significance of the gravitational core-mantle interaction on Earth's rotation variation, we introduce in our MoSST core dynamics model a heterogeneous mantle, with a density distribution derived from seismic results. In this model, the core convection is driven by the buoyancy forces. And the density variation is determined dynamically with the convection. Numerical simulation is carried out with different parameter values, intending to extrapolate numerical results for geophysical implications.

Description: Aerosol absorption in the atmosphere poses a major uncertainty in assessing the aerosol climate effects. This uncertainty arises from the poorly quantified aerosol sources, especially black carbon emissions, and our limited knowledge of aerosol mixing state and optical properties. Here we use a global model GOCART to simulate atmospheric aerosols, including sulfate, black carbon, organic carbon, dust, and sea salt. We compare the model calculated total aerosol optical thickness, extinction, and absorption with those quantities from the ground-based sun photometer measurements from AERONET at several different wavelengths and the field observations from ACE-Asia. We will examine what are the most sensitive factors in determining the aerosol absorption, and the consequences of assessing the aerosol radiative forcing and atmospheric heating associated with those factors.

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

Description: Stress 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: c1inoptilolite-rich tuff-hydroxyapatite mixture (zeoponic substrate) has the potential to serve as a synthetic soil-additive for plant growth. Essential plant macro-nutrients such as calcium, phosphorous, magnesium, ammonium and potassium are released into solution via dissolution of the hydroxyapatite and cation exchange on zeolite charged sites. Plant growth experiments resulting in low yield for wheat have been attributed to a Ca deficiency caused by a high degree of cation exchange by the zeolite. Batch-equilibration experiments were performed in order to determine if the Ca deficiency can be remedied by the addition of a second Ca-bearing, soluble, mineral such as calcite, dolomite or wollastonite. Variations in the amount of calcite, dolomite or wollastonite resulted in systematic changes in the concentrations of Ca and P. The addition of calcite, dolomite or wollastonite to the zeoponic substrate resulted in an exponential decrease in the phosphorous concentration in solution. The exponential rate of decay was greatest for calcite (5.60 wt. % -I), intermediate for wollastonite (2.85 wt.% -I) and least for dolomite (1.58 wt.% -I). Additions of the three minerals resulted in linear increases in the calcium concentration in solution. The rate of increase was greatest for calcite (3.64), intermediate for wollastonite (2.41) and least for dolomite (0.61). The observed changes in P and Ca concentration are consistent with the solubilities of calcite, dolomite and wollastonite and with changes expected from a common ion effect with Ca. Keywords: zeolite, zeoponics, common-ion effect, clinoptilolite, hydroxyapatite

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: In this paper we present the first early dry season (early June-early August) emission factor measurements for carbon dioxide (CO2), carbon monoxide (CO), methane (Ca), nonmethane hydrocarbons (NMHC), and particulates with a diameter less than 2.5 microns (pM2.5) for southern African grassland and woodland fires. Seasonal emission factors for grassland fires correlate linearly with the proportion of green grass, used as a surrogate for the fuel moisture content, and are higher for products of incomplete combustion in the early part of the dry season compared with later in the dry season. Models of emission factors for NMHC and PM(sub 2.5) versus modified combustion efficiency (MCE) are statistically different in grassland compared with woodland ecosystems. We compare predictions based on the integration of emissions factors from this study, from the southern African Fire-Atmosphere Research Initiative 1992 (SAFARI-92), and from SAFARI-2000 with those based on the smaller set of ecosystem-specific emission factors to estimate the effects of using regional-average rather than ecosystem-specific emission factors. We also test the validity of using the SAFARI-92 models for emission factors versus MCE to predict the early dry season emission factors measured in this study. The comparison indicates that the largest discrepancies occur at the low end (0.907) and high end (0.972) of MCE values measured in this study. Finally, we combine our models of MCE versus proportion of green grass for grassland fires with emission factors versus MCE for selected oxygenated volatile organic compounds measured in the SAFARI-2000 campaign to derive the first seasonal emission factors for these compounds. The results of this study demonstrate that seasonal variations in savanna fire emissions are important and should be considered in modeling emissions at regional to continental scales.

Description: Airborne measurements made on initial smoke from 10 savanna fires in southern Africa provide quantitative data on emissions of 50 gaseous and particulate species, including carbon dioxide, carbon monoxide, sulfur dioxide, nitrogen oxides, methane, ammonia, dimethyl sulfide, nonmethane organic compounds, halocarbons, gaseous organic acids, aerosol ionic components, carbonaceous aerosols, and condensation nuclei (CN). Measurements of several of the gaseous species by gas chromatography and Fourier transform infrared spectroscopy are compared. Emission ratios and emission factors are given for eight species that have not been reported previously for biomass burning of savanna in southern Africa (namely, dimethyl sulfide, methyl nitrate, five hydrocarbons, and particles with diameters from 0.1 to 3 microns). The emission factor that we measured for ammonia is lower by a factor of 4, and the emission factors for formaldehyde, hydrogen cyanide, and CN are greater by factors of about 3, 20, and 3 - 15, respectively, than previously reported values. The new emission factors are used to estimate annual emissions of these species from savanna fires in Africa and worldwide.

Description: Retrieved tropospheric carbon monoxide (CO) column densities are presented for more than 9000 spectra obtained by the University of Wisconsin-Madison (UWis) Scanning High-Resolution Interferometer Sounder (SHIS) during a flight on the NASA ER-2 on 7 September 2000 as part of the Southern African Regional Science Initiative (SAFARI 2000) dry season field campaign. Enhancements in tropospheric column CO were detected in the vicinity of a controlled biomass burn in the Timbavati Game Reserve in northeastern South Africa and over the edge of the river of smoke in south central Mozambique. Relatively clean air was observed over the far southern coast of Mozambique. Quantitative comparisons are presented with in situ measurements from five different instruments flying on two other aircraft: the University of Washington Convair-580 (CV) and the South African Aerocommander JRB in the vicinity of the Timbavati fire. Measured tropospheric CO columns (extrapolated from 337 to 100 mb) of 2.1 x 10(exp 18) per square centimeter in background air and up to 1.5 x 10(exp 19) per square centimeter in the smoke plume agree well with SHIS retrieved tropospheric CO columns of (2.3 plus or minus 0.25) x 10(exp 18) per square centimeter over background air near the fire and (1.5 plus or minus 0.35) x 10(exp 19) per square centimeter over the smoke plume. Qualitative comparisons are presented with three other in situ CO profiles obtained by the South African JRA aircraft over Mozambique and northern South Africa showing the influence of the river of smoke.

Description: During the flights of the University of Washington's Convair-580 in the Southern African Regional Science Initiative (SAFARI 2000) in southern Africa, a phenomenon was observed that has not been reported previously. This was the occurrence of thin layers of remarkably clean air, sandwiched between heavily polluted air, which persisted for many hours during the day. Photographs are shown of these clean air slots (CAS), and particle concentrations and light scattering coefficients in and around such slot are presented. An explanation is proposed for the propensity of CAS to form in southern Africa during the dry season.

Description: During the dry season airborne campaign of the Southern African Regional Science Initiative (SAFARI 2000), coordinated observations were made of massive thick aerosol layers. These layers were often dominated by aerosols from biomass burning. We report on airborne Sun photometer measurements of aerosol optical depth (lambda = 0.354- 1.557 microns), columnar water vapor, and vertical profiles of aerosol extinction and water vapor density that were obtained aboard the University of Washington's Convair-580 research aircraft. We compare these with ground-based AERONET Sun/sky radiometer results, with ground based lidar data (MPL-Net), and with measurements from a downward pointing lidar aboard the high-flying NASA ER-2 aircraft. Finally, we show comparisons between aerosol optical depths fiom the Sun photometer and those retrieved over land and over water using four spaceborne sensors (TOMS, MODIS, MISR, and ATSR-2).