ALBERT

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

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

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

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

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

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

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

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

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

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

Description: 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 52 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 @e., much higher harmonic degrees). The most important example is the GRACE mission launched in March 2002, following the success of the CHAMP mission. Such observations are becoming a new and powerful tool for remote sensing of geophysical fluid processes that involve larger-scale mass transports.

Description: During the SAFARI 2000 field campaign, both smoke aerosols from savanna fires and haze aerosols in the boundary layer and in the free troposphere were collected from an aircraft in southern Africa. These aerosol samples were analyzed for their water-soluble chemical components, particularly the organic species. A novel technique, electrospray ionization-ion trap mass spectrometry, was used concurrently with an ion chromatography system to analyze for carbohydrate species. Seven carbohydrates, seven organic acids, five metallic elements, and three inorganic anions were identified and quantified. On the average, these 22 species comprised 36% and 27% of the total aerosol mass in haze and smoke aerosols, respectively. For the smoke aerosols, levoglucosan was the most abundant carbohydrate species, while gluconic acid was tentatively identified as the most abundant organic acid. The mass abundance and possible source of each class of identified species are discussed, along with their possible formation pathways. The combustion phase of a fire had an impact on the chemical composition of the emitted aerosols. Secondary formation of sulfate, nitrate, levoglucosan, and several organic acids occurred during the initial aging of smoke aerosols. It is likely that under certain conditions, some carbohydrate species in smoke aerosols, such as levoglucosan, were converted to organic acids during upward transport.

Description: During the dry season component of the Southern African Regional Science Initiative (SAFARI) in late winter 2000, the net solar spectral irradiance was measured at flight levels throughout biomass burning haze layers. From these measurements, the flux divergence, fractional absorption, instantaneous heating rate, and absorption efficiency were derived. Two cases are examined: on 24 August 2000 off the coast of Mozambique in the vicinity of Inhaca Island and on 6 September 2000 in a very thick continental haze layer over Mongu, Zambia. The measured absolute absorption was substantially higher for the case over Mongu where the measured midvisible optical depth exceeded unity. Instantaneous heating from aerosol absorption was 4 K d(sup -1) over Mongu, Zambia and 1.5 K d(sup -1) near Inhaca Island, Mozambique. However, the spectral absorption efficiency was nearly identical for both cases. Although the observations over Inhaca Island preceded the river of smoke from the southern African continent by nearly 2 weeks, the evidence here suggests a continental influence in the lower tropospheric aerosol far from source regions of burning.

Description: Individual aerosol particles collected over southern Africa during the SAFARI 2000 field study were studied using transmission electron microscopy and field-emission scanning electron microscopy. The sizes, shapes, compositions, mixing states, surface coatings, and relative abundances of aerosol particles from biomass burning, in boundary layer hazes, and in the free troposphere were compared, with emphasis on aging and reactions of inorganic smoke particles. Potassium salts and organic particles were the predominant species in the smoke, and most were internally mixed. More KCl particles occur in young smoke, whereas more K2SO4 and KNO3 particles were present in aged smoke. This change indicates that with the aging of the smoke, KCl particles from the fires were converted to K2SO4 and KNO3 through reactions with sulfur- and nitrogen- bearing species from biomass burning as well as other sources. More soot was present in smoke from flaming grass fires than bush and wood fires, probably due to the predominance of flaming combustion in grass fires. The high abundance of organic particles and soluble salts can affect the hygroscopic properties of biomass-burning aerosols and therefore influence their role as cloud condensation nuclei. Particles from biomass burning were important constituents of the regional hazes.

Description: Individual aerosol particles in smoke plumes from biomass fires and in regional hazes in southern Africa were studied using analytical transmission electron microscopy (TEM), which allowed detailed characterization of carbonaceous particle types in smoke and determination of changes in particle properties and concentrations during smoke aging. Based on composition, morphology, and microstructure, three distinct types of carbonaceous particles were present in the smoke: organic particles with inorganic (K-salt) inclusions, tar ball particles, and soot. The relative number concentrations of organic particles were largest in young smoke, whereas tar balls were dominant in a slightly aged (1 hour) smoke from a smoldering fire. Flaming fires emitted relatively more soot particles than smoldering fires, but soot was a minor constituent of all studied plumes. Further aging caused the accumulation of sulfate on organic and soot particles, as indicated by the large number of internally mixed organic/sulfate and soot/sulfate particles in the regional haze. Externally mixed ammonium sulfate particles dominated in the boundary layer hazes, whereas organic/sulfate particles were the most abundant type in the upper hazes. Apparently, elevated haze layers were more strongly affected by biomass smoke than those within the boundary layer. Based on size distributions and the observed patterns of internal mixing, we hypothesize that organic and soot particles are the cloud-nucleating constituents of biomass smoke aerosols. Sea-salt particles dominated in the samples taken in stratus clouds over the Atlantic Ocean, off the coast of Namibia, whereas a distinct haze layer above the clouds consisted of aged biomass smoke particles.

Description: During a polarity transition of the Earth's magnetic field, the structure and strength of the field change significantly from their present values. This will alter the global pattern of charged particle precipitation into the atmosphere. Thus, particle precipitation is possible into regions that are at the moment effectively shielded by the Earth's magnetic field. A two-dimensional global fbHy coupled chemistry, radiation, and transport model of the atmosphere has been used to investigate how the increased particle precipitation affects the chemical composition of the middle and lower atmosphere. Ozone losses resulting from large energetic particle events are found to increase significantly, with resultant losses similar to those observed m the Antarctic ozone hole of the 1990s. This results in significant increases in surface UV-B radiation as well as changes in stratospheric temperature and circulation over a period of several months after large particle events.

Description: We analyze two LLBL crossings made by Interball Tail satellite under southward or variable magnetosheath magnetic field: one crossing on the flank of the magnetosphere, and another one closer to the subsolar point. Three different types of ion velocity distributions within LLBL are observed: (a) D-shaped distributions, (b) ion velocity distributions consisting of two counter-streaming components of magnetosheath-ty and (c) distributions with three components one of which has nearly zero parallel velocity and two counter-streaming components. Only the (a) type fits to the single magnetic flux tube formed by reconnection between magnetospheric and magnetosheath magnetic fields. We argue that two counter-streaming magnetosheath-like ion components observed by Interball within LLBL cannot be explained by the reflection of the ions from the magnetic mirror deeper within magnetosphere. Types (b) and (c) ion velocity distributions would form within spiral magnetic flux tube consisting of a mixture of alternating segments originating from the magnetosheath and from magnetospheric plasma. The shapes of ion velocity distributions and their evolution with decreasing number density in LLBL indicate that a significant part of LLBL is located on magnetic field lines of long spiral flux tube islands at the magnetopause, as has been proposed and found to occur in magnetopause simulations. We consider these observations as evidence for multiple reconnection X-ray lines between magnetosheath and magnetospheric flux tubes.

Description: The marine record of platinum group elements (PGEs) and Os isotopic compositions provides information on different inputs of PGEs into the oceans. Some studies based on a smaller subset of the PGEs suggest that the PGEs may suffer post-depositional mobility during diagenesis. In some K/T boundary clays, Kyte and others showed that the relative abundances of Pt, Pd, Ir, and Os can differ significantly from chondritic, which is the signature expected from fallout of the meteorite impact. In some K/T boundary sections, elevated Ir concentrations are observed as far as 1 meter from the cm-thick boundary clay containing the meteoritic ejecta. The purpose of this study was to characterize Pt, Pd, Ir, Os, and Ru abundances in zones including the K/T boundary. We determined PGE abundances of boundary clays at two hemipelagic sites (Stevns Klint, Denmark and Caravaca, Spain) in which previous studies by Kyte and others showed that the Ir anomaly is confined to within a few cm. We also analyzed two pelagic Pacific sites: a boundary clay from the north Pacific (Hole 465A) characterized by a 0.5 m thick Ir anomaly and a transect across the K/T boundary from the south Pacific (Hole 596) where the Ir anomaly spans 2 m. The Stevns Klint, Caravaca, and north Pacific sites are characterized by abundant marls and limestones in the section, whereas the south Pacific site is dominated by clays. Samples were spiked with isotopic tracers, mixed with a flux, S and Ni, and equilibrated by fusion. PGEs were extracted from the Ni and analyzed on a Finnigan Element ICP-MS. We find that the narrow Caravaca and Stevns Klint boundary clays have relative PGE abundance patterns indistinguishable from chondritic values. The two Pacific sites were found to have nearly identical PGE patterns but have ratios at the peak, which differ from chondritic values as found earlier by Evans et al. The Pacific sites were found to have nearly identical PGE patterns but are extremely depleted in OS (Os/Ir = 0.07-0.15) and slightly enriched in Pd and Pt relative to Ir.

Description: Observations have demonstrated that length of day (LOD) variation on decadal time scales results from exchange of axial angular momentum between the solid mantle and the core. There are in general four core-mantle interaction mechanisms that couple the core and the mantle. Of which, three have been suggested likely the dominant coupling mechanism for the decadal core-mantle angular momentum exchange, namely, gravitational core-mantle coupling arising from density anomalies in the mantle and in the core (including the inner core), the electromagnetic coupling arising from Lorentz force in the electrically conducting lower mantle (e.g. D-layer), and the topographic coupling arising from non-hydrostatic pressure acting on the core-mantle boundary (CMB) topography. In the past decades, most effort has been on estimating the coupling torques from surface geomagnetic observations (kinematic approach), which has provided insights on the core dynamical processes. In the meantime, it also creates questions and concerns on approximations in the studies that may invalidate the corresponding conclusions. The most serious problem is perhaps the approximations that are inconsistent with dynamical processes in the core, such as inconsistencies between the core surface flow beneath the CMB and the CMB topography, and that between the D-layer electric conductivity and the approximations on toroidal field at the CMB. These inconsistencies can only be addressed with numerical core dynamics modeling. In the past few years, we applied our MoSST (Modular, Scalable, Self-consistent and Three-dimensional) core dynamics model to study core-mantle interactions together with geodynamo simulation, aiming at assessing the effect of the dynamical inconsistencies in the kinematic studies on core-mantle coupling torques. We focus on topographic and electromagnetic core-mantle couplings and find that, for the topographic coupling, the consistency between the core flow and the CMB topography is critical for correct evaluation of the coupling torque.

Description: Modeling of the Earth's core has reached a level of maturity to where the incorporation of observations into the simulations through data assimilation has become feasible. Data assimilation is a method by which observations of a system are combined with a model output (or forecast) to obtain a best guess of the state of the system, called the analysis. The analysis is then used as an initial condition for the next forecast. By doing assimilation, not only we shall be able to predict partially secular variation of the core field, we could also use observations to further our understanding of dynamical states in the Earth's core. One of the first steps in the development of an assimilation system is a comparison between the observations and the model solution. The highly turbulent nature of core dynamics, along with the absence of any regular external forcing and constraint (which occurs in atmospheric dynamics, for example) means that short time comparisons (approx. 1000 years) cannot be made between model and observations. In order to make sensible comparisons, a direct insertion assimilation method has been implemented. In this approach, magnetic field observations at the Earth's surface have been substituted into the numerical model, such that the ratio of the multiple components and the dipole component from observation is adjusted at the core-mantle boundary and extended to the interior of the core, while the total magnetic energy remains unchanged. This adjusted magnetic field is then used as the initial field for a new simulation. In this way, a time tugged simulation is created which can then be compared directly with observations. We present numerical solutions with and without data insertion and discuss their implications for the development of a more rigorous assimilation system.

Description: Using the satellite-laser-ranging (SLR) data, Cox and Chao [2002] reported the detection of a large post-1998 anomaly (in the form of a positive jump) in the time series of Earth s lowest-degree gravity harmonic 52, or the dynamic oblateness. Among several groups now examining the mass redistribution in the global geophysical fluids in search of the cause(s), we report here a temporally coinciding anomalies found in the extratropic north + south Pacific basins. Clearly seen in the leading EOFPC mode for extratropic Pacific, these anomalies occurred in sea-surface height, sea-surface temperature, and temperature- and salinity-depth profiles. We based our analysis on two different data sources: TOPEX/Poseidon altimetry, and the ECCO ocean general circulation model output assimilating T/P data. The magnitude of these changes, when converted to equivalent J2 change, appears to be a few times too small to explain the observed J2 directly. These findings, and the fact that the anomalies occurred following the strong 1997-98 El Nino, suggest strong geophysical connection of the interannual-to-decadal variation of 52 with the Pacific Decadal Oscillation (PDO) and the ultimate global-change processes that cause PDO. More work is underway, and additional independent data sources are examined, paying close attention to the fact that the J2 anomaly has been reversing back to normal since 2001. These include: (1) cryospheric contributions (melting of glaciers and ice sheets); (2) land hydrological contributions; (3) polar sea influences ( e g , via deep flow); (4) fluid flow in Earth's core; (5) time-variable gravity signals from SLR in higher harmonic degree/order, including J3,J4, (2,1), and (2,2) coefficients, considering their lower signal-to-noise ratios; (6) Earth rotation data in terms of length-of-day and polar motion.

Description: Tidally-induced stick-slip motion in the mouth of Whillans Ice Stream provides a unique natural experiment in ice-stream response behavior and fiom which we might learn a great deal about subglacial till properties and sub-ice-stream conditions. At the IGS Symposium on Fast Glacier Flow (Yakutat, 2002), we reported our observations of stick- slip motion and demonstrated its synchronicity with tidal forcing. Recently, we have completed additional processing of our GPS data in differential mode. It reveals more details of the stick-slip events and illustrates that within 30 seconds, the temporal interval of our data, the ice stream accelerates to a speed corresponding to a completely lubricated bed. While details of individual events vary, there seems to be strong evidence of an elastic rebound on the time scale of one hour following most events. This suggests the event involves the release of stored elastic strain energy in the ice. The similar displacements of events suggest further that till or subglacial hydrologic properties limit the amount of elastic strain released in any single event. We follow a line of reasoning that dilatant strengthening limits the slip displacement and present model of the stick-slip process. To match the observed delay between the peak ocean tide and stick-slip events, our model includes a propagating pressure wave in the subglacial hydrologic system between the grounding line, where the rising tide first increases the subglacial water pressure and regions upstream where stored elastic strain increases the basal shear stress. This high-tide event is released when the increased water pressure reaches the region of increased shear stress. Dilatant strengthening stops the event by increasing pore volume and lowering the water pressure. Following this event, falling tide increases the normal forces, compresses the till and increases pore pressure again, leading to the second falling-tide event we observe every tidal cycle.

Description: The excitation of lower hybrid waves (LHWs) is a widely discussed mechanism of interaction between plasma species in space, and is one of the unresolved questions of magnetospheric multi-ion plasmas. In this paper we present the morphology, dynamics, and level of LHW activity generated by electromagnetic ion cyclotron (EMIC) waves during the May 2-7, 1998 storm period on the global scale. The LHWs were calculated based on our newly developed self-consistent model that couples the system of two kinetic equations: one equation describes the ring current (RC) ion dynamic, and another equation describes the evolution of EMIC waves. It is found that the LHWs are excited by helium ions due to their mass dependent drift in the electric field of EMIC waves. The level of LHW activity is calculated assuming that the induced scattering process is the main saturation mechanism for these waves. The calculated LHWs electric fields are consistent with the observational data.

Description: During the Second SAGE 111 Ozone Loss and Validation Experiment (SOLVE II), the 14- channel NASA Ames Airborne Trackmg Sunphotometer (AATS-14) was mounted on the NASA DC-8 and successfully measured spectra of total and aerosol optical depth (TOD and AOD) during the sunlit portions of eight science flights. Values of ozone column content above the aircraft have been derived from the AATS-14 data by using a linear least squares method. For each AATS-14 measured TOD spectrum, this method iteratively finds the ozone column content that yields the best match between measured and calculated TOD. The calculations assume the known Chappuis ozone band shape and a three-parameter AOD shape (quadratic in log-log space). Seven of the AATS-14 channels (each employing an interference filter with a nominal full-width at half maximum bandpass of -5 nm) are within the Chappuis band, with center wavelengths between 452.9 nm and 864.5 nm. One channel (604.4 nm) is near the peak, and three channels (499.4, 519.4 and 675.1 nm) have ozone absorption within 30-40% of that at the peak. For the typical DC-8 SOLVE II cruising altitudes of approx. 8-12 km and the background stratospheric aerosol conditions that prevailed during SOLVE 11, absorption of incoming solar radiation by ozone comprised a significant fraction of the aerosol-plus-ozone optical depth measured in the four AATS-14 channels centered between 499.4 and 675.1 nm. Typical AODs above the DC-8 ranged from 0.003-0.008 in these channels. For comparison, an ozone overburden of 0.3 atm-cm (300 DU) translates to ozone optical depths of 0.009,0.014, 0.041, and 0.012, respectively, at these same wavelengths. In this paper, we compare AATS-14 values of ozone column content with temporally and spatially near-coincident values derived from measurements acquired by the Stratospheric Aerosol and Gas Experiment III (SAGE III) and the Polar Ozone and Aerosol Measurement 111 (POAM III) satellite sensors. We also compare AATS-14 ozone retrievals during selected DC-8 latitudinal and longitudinal transects with total column ozone data acquired by the Total Ozone Mapping Spectrometer (TOMS) and the Global Ozone Monitoring Experiment (GOME) satellite sensors. To enable this comparison, the amount of ozone in the column below the aircraft is estimated by combining SAGE and/or POAM data with high resolution, fast response in-situ ozone measurements acquired during the DC-8 ascent at the start of each science flight.

Description: Flux transfer events (FTE) have been postulated to result from transient magnetic merging. If so, the ion distributions within an event should exhibit features known to result from merging. Observations of a FTE by instruments on the Polar spacecraft revealed classical merging signatures that included: 1) D-shaped, accelerated, magnetosheath ion distributions, 2) a well defined de Hoffman-Teller frame, 3) local stress balance, and 4) a P-N magnetic field signature. This FTE was observed near the magnetic equator at approx. 13 MLT under conditions of a moderately northward interplanetary magnetic field (IMF) (clock angle of less than 10 deg). The nature of the ion distributions and the consistency of the measured cutoff speed with that calculated from the measured local magnetic field and the derived de Hoffman-Teller speed show the ion injection to be local. Coupled with the northward IMF these results lead to the conclusion that component merging in the low latitude region was responsible for the FTE.

Description: To determine the effect of crustal thickness variation on satellite-altitude geopotential anomalies we compared two regions of Europe with vastly different values, Central/Southern Finland and the Pannonian Basin. Crustal thickness exceeds 62 km in Finland and is less than 26 km in the Pannonian Basin. Heat-flow maps indicate that the thinner and more active crust of the Pannonian Basin has a value nearly three times that of the Finnish Svecofennian Province. 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 are complicated by de-glaciation. CHAMP gravity data (400 km) indicates a west-east positive gradient of greater than 4 mGal across Central/Southern Finland and an ovoid positive anomaly (approximately 4 mGal) quasi-coincidental with the magnetic anomaly traversing the Pannonian Basin. CHAMP magnetic data (425 km) reveal elongated semicircular 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 both regions subducted oceanic lithosphere has been proposed as the anomalous body.

Description: Complete description of a self-consistent model for 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 non steady-state conditions. In the paper by Khazanov et al. [2002] this self-consistent model has only been shortly outlined, and discussions of many the model related details have been omitted. For example, in present study for the first time a new algorithm for numerical finding of the resonant numbers for quasilinear wave-particle interaction is described, or it is demonstrated that in order to describe quasilinear interaction in a multi-ion thermal plasma correctly, both e and He(+) modes of electromagnetic ion cyclotron waves should be employed. The developed model is used to simulate the entire May 2-7, 1998 storm period. Trapped number fluxes of the ring current protons are calculated and presented along with their comparison with the data measured by the 3D hot plasma instrument Polar/HYDRA. Examining of the wave (MLT, L shell) distributions produced during the storm progress reveals an essential intensification of the wave emissions in about two days after main phase of storm. This result is well consistent with the earlier ground-based observations. Also the theoretical shapes and the occurrence rates for power spectral densities of electromagnetic ion cyclotron waves are studied. It is found that in about 2 days after the storm main phase on May 4, mainly non Gaussian shapes of power spectral densities are produced.

Description: The analysis of global VLBI observations provides vertical rates for 50 sites with formal errors less than 2 mm/yr and median formal error of 0.4 mm/yr. These sites are largely in Europe and North America with a few others in east Asia, Australia, South America and South Africa. The time interval of observations is up to 20 years. The error of the velocity reference frame is less than 0.5 mm/yr, but results from several sites with observations from more than one antenna suggest that the estimated vertical rates may have temporal variations or non-geophysical components. Comparisons with GPS rates and corresponding site position time series will be discussed.

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

Description: Our research was targeted at three issues: (1) the factors controlling ozone in the tropical troposphere, (2) the Asian outflow of ozone and its precursors, and (3) the causes of decadal trends observed in ozone and CO. We have also used support from this ACMAP grant to (1) work with Kelly Chance on the retrieval and interpretation of HCHO and NO2 observations from GOME, and (2) develop GEOS-CHEM into a versatile model supporting the work of a large number of users including outside Harvard. ACMAP has provided the core support for GEOS-CHEM development. Applications of the GEOS-CHEM model with primary support from ACMAP are discussed below. A list of publications resulting from this grant is given at the end of the report.

Description: The meteorological pathways contributing to Asian pollution outflow over the Pacific are examined with a global three-dimensional model analysis of CO observations from the Transport and Chemical Evolution over the Pacific (TRACE-P) aircraft mission (February-April 2001). The model is used also to place the TRACE-P observations in an interannual (1994-2001) and seasonal context. The major process driving Asian pollution outflow in spring is frontal lifting ahead of southeastward-moving cold fronts (the leading edge of cold surges) and transport in the boundary layer behind the cold fronts. Orographic lifting over central and eastern China combines with the cold fronts to promote the transport of Chinese pollution to the free troposphere. Outflow of seasonal biomass burning in Southeast Asia during spring takes place mostly by deep convection but also by northeastward transport and frontal lifting, mixing with the anthropogenic outflow. Boundary layer outflow over the western Pacific is largely devoid of biomass burning influence. European and African (biomass burning) plumes in Asian outflow during TRACE-P were weak (less than 60 ppbv and 20 ppbv CO, respectively) and were not detectable in the observations because of superposition of the much larger Asian pollution signal. Spring 2001 (La Nina) was characterized by unusually frequent cold surge events in the Asian Pacific rim and strong convection in Southeast Asia, leading to unusually strong boundary layer outflow of anthropogenic emissions and convective outflow of biomass burning emissions in the upper troposphere. The Asian outflow flux of CO to the Pacific is found to vary seasonally by a factor of 3-4 (maximum in March and minimum in summer). The March maximum results from frequent cold surge events and seasonal biomass burning emissions.

Description: This paper describes the data analysis technique used for magnetic testing at the NASA Goddard Space Flight Center (GSFC). Excellent results have been obtained using this technique to convert a spacecraft s measured magnetic field data into its respective magnetic dipole moment model. The model is most accurate with the earth s geomagnetic field cancelled in a spherical region bounded by the measurement magnetometers with a minimum radius large enough to enclose the magnetic source. Considerably enhanced spacecraft magnetic testing is offered by using this technique in conjunction with a computer-controlled magnetic field measurement system. Such a system, with real-time magnetic field display capabilities, has been incorporated into other existing magnetic measurement facilities and is also used at remote locations where transport to a magnetics test facility is impractical.

Description: On January 12, 2003 NASA launched the first satellite-based lidar, the Geoscience Laser -Altimeter System (GLAS), onboard the ICESat spacecraft. The GLAS atmospheric measurements introduce a fundamentally new and important tool for understanding the atmosphere and climate. In the past, aerosols have only been studied from space using images gathered by passive sensors. Analysis of this passive data has lead to an improved understanding of aerosol properties, spatial distribution, and their effect on the earth's climate. However, these images do not show the aerosol's vertical distribution. As a result, a key piece of information has been missing. The measurements now obtained by GLAS will provide information on the vertical distribution of aerosols and clouds, and improve our ability to study their transport processes and aerosol-cloud interactions. Here we show an overview of GLAS, provide an update of its current status, and present initial observations of dust profiles. In particular, a strategy of characterizing the height profile of dust plumes over source regions will be presented.

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 Coordinated Observations of Dynamics in Aurora (CODA) experiment for which Dr. Andy Christensen of the Aerospace Corp. was the Principal Investigator. Clemson University was a co-investigator institution on the project and was tasked with designing, building, and testing a total of four chemical tracer release payloads which were used to release trimethyl aluminum (TMA) to measure the winds and turbulent diffusion. All four payloads were flown from the rocket range at Poker Flat, Alaska, in support of instrumented payloads flown from the same location. The CODA experiment was designed to investigate the changes in the chemical composition in the lower thermosphere during an auroral substorm event. In particular, the objective was to determine the relative importance of two competing mechanisms. The first mechanism is associated with the advection of air from the polar cap into the auroral oval. A flow in that direction is typical in the postmidnight sector of the oval. As the air moves from the region of negligible forcing in the polar cap to the region of strong forcing and heating in the oval, changes in the composition are expected to occur. Alternatively, enhanced local mixing will also cause such changes. Since the auroral forcing accelerates the neutral winds in the lower E region, but not at lower altitudes, the forcing results in large shears that can be large enough to be unstable in the Richardson number sense. The unstable shears are expected to generate turbulence and therefore enhanced mixing. Such mixing will also create composition changes similar to what is observed. In the CODA experiment an instrumented payload designed and built by the Aerospace Corporation was used to measure the composition and plasma parameters needed to define the auroral forcing. In addition, three chemical tracer payloads were to be flown along three trajectories widely separated in azimuth to map out the winds and wind gradients. The measurements would thus provide information about the local composition changes above the launch site as well as detailed information about the advection and location of stable or unstable shears both locally and over a broader area in the vicinity of the launch site.

Description: The Advanced Thin Ionization Calorimeter (ATIC) balloon experiment completed a 19.7 day science flight, 29 Dec. 02 to 18 Jan. 03, from McMurdo, Antarctica, yielding 65 Gbytes of science data. The instrument was fully recovered just over the trans-Antarctic mountain range. ATIC is designed to make new measurements of the composition and energy spectra of cosmic rays in the region up to 100 TeV total energy. The flight and the instrument performance are described and compared to the ATIC-1 flight from 2000-01. Although the data from this new flight will not be fully processed and analyzed, some preliminary results are discussed.

Description: The goal of this project is to determine precursors and signatures of local substorm onset and how they evolve in the plasma sheet using the Geotail near-Earth database. This project is part of an ongoing investigation involving this PI, Nelson Maynard (Mission Research Corporation), and William Burke (AFRL) toward an empirical understanding of the onset and evolution of substorms. The first year began with dissemination of our CRRES findings, which included an invited presentation and major publication. The Geotail investigation began with a partial survey of onset signature types at distances X less than 15 R(sub E) for the first five months (March-July 1995) of the Geotail near-Earth mission. During the second year, Geotail data from March 1995 to present were plotted. Various signatures at local onset were catalogued for the period through 1997. During this past year we performed a survey of current-disruption-like (CD-like) signatures at distances X less than or equal to 14 R(sub E) for the three years 1995-1997.

Description: We have investigated the energy distributions of the metastable oxygen atoms in the terrestrial thermosphere. Nascent O(lD) atoms play a fundamental role in the energy balance and chemistry of the terrestrial atmosphere, because they are produced by photo-chemical reactions in the excited electronic states and carry significant translational energies.

Description: Silicate inclusions in nine known IIE irons show diversity in mineralogy, and Colomera, Kodaikanal, Elga and Miles contain alkali-rich silicate inclusions. Bogard et al. showed evidence of a complex parent body evolution for IIE irons based on Ar-39-Ar-40 ages. Colomera contained a sanidine-rich surface inclusion and the K-enrichment trends in the Na-rich inclusions are different from those of other IIEs. To elucidate the origin of K-rich materials, we studied the mineralogy and Ar-Ar age of silicate inclusions from the Tarahumara IIE iron meteorite.

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.