ALBERT

Description: The delivery of amino acids by micrometeorites to the early Earth during the period of heavy bombardment could have been a significant source of the Earth's prebiotic amino acid inventory provided that these organic compounds survived atmospheric entry heating. To investigate the sublimation of amino acids from a micrometeorite analog at elevated temperature, grains from the CM-type carbonaceous chondrite Murchison were heated to 550 C inside a glass sublimation apparatus (SA) under reduced pressure. The sublimed residue that had collected on the cold finger of the SA after heating was analyzed for amino acids by HPLC. We found that when the temperature of the meteorite reached approx. 150 C, a large fraction of the amino acid glycine had vaporized from the meteorite, recondensed onto the end of the SA cold finger, and survived as the rest of the grains heated to 550 C. alpha-Aminoisobutryic acid and isovaline, which are two of the most abundant non-protein amino acids in Murchison, did not sublime from the meteorite and were completely destroyed during the heating experiment. Our experimental results suggest that sublimation of glycine present in micrometeorite grains may provide a way for this amino acid to survive atmospheric entry heating at temperatures less than 550 C; all other amino acids apparently are destroyed. Key Words: Amino acids-Exogenous delivery-Micrometeorites-Sublimation.

Description: Many field studies have been conducted that document the morphology of ventifacts and the directionality of their features relative to current and past wind regimes. Field plots and wind tunnel studies have identified heights and particle concentrations above the surface where maximum abrasion occurs. However, as of yet, the rates and detailed methods by which rocks abrade and evolve into ventifacts are poorly documented and understood. This abstract addresses this gap in knowledge by interpreting controlled laboratory and field analog studies. We begin with an overview of the methods by which the wind tunnel experiments and field studies were done, followed by how the resulting data were analyzed and interpreted. A presentation of the results comes next, after which the implications for rock abrasion and ventifact formation on Earth and Mars are discussed. We show that initial rock shape and texture play important roles in determining both rate and style of abrasion, with steep-sided, rough rocks eroding the fastest but with intermediate-angled faces exhibiting the greatest shape change. Most rocks tend to evolve toward an equilibrium shape whose form is poorly conducive to further abrasion. Most rocks on Mars and in terrestrial ventifact localities never reach this mature state, with erosion ceasing or slowing down due to exhaustion of the sand supply and other factors.

Description: New analyses of rocks and soils at the Mars Pathfinder landing site have been completed using the full Imager for Mars Pathfinder (IMP) 12- color SuperPan panorama. These revise early conclusions that rocks at the landing site are a single lithology coated only by windblown dust. We conclude instead that there is also a second lithology in addition to the dominant gray rock, and that it is consistent with highlands material excavated from beneath a thin veneer of northern plains; that many rocks have cemented coatings that formed during an early, probably wetter climate; and that young rocks excavated after coating formation ceased are mainly breccias or conglomerates.

Description: We investigate the source of poleward propagating atmospheric zonal wind anomalies, originating at the equator and penetrating to high latitudes in both hemispheres in conjunction with ENSO [Dickey et al.,1992], and report the discovery of similar variability on decadal and longer timescales. Since atmospheric dissipation times are generally on the order of a month or less, we examine the ocean as a 'memory' source for these globally coherent anomalies. This hypothesis is substantiated by the observation of complementaryoscillation in the sea surface temperature (SST) field; further, we detect a robust decadal variability (~1012 yrs)in both the SST and contemporaneous atmospheric angular momentum (AAM) series. Analyzing GISST SST data beginning in 1902, we confirm this decadal mode and find signatures of longer (multidecadal) SST variability centered in the equatorial and North Pacific.

Description: As a result of the research leading to the 1998 AIR workshop and the subsequent analysis, the neutron issues posed by Foelsche et al. and further analyzed by Hajnal have been adequately resolved. We are now engaged in developing a new atmospheric ionizing radiation (AIR) model for use in epidemiological studies and air transportation safety assessment. A team was formed to examine a promising code using the basic FLUKA software but with modifications to allow multiple charged ion breakup effects. A limited dataset of the ER-2 measurements and other cosmic ray data will be used to evaluate the use of this code.

Description: We have investigated the importance of intercontinental transport using source-receptor relationship. A global radon-like and seven regional tracers were used in three-dimensional model simulations to quantify their contributions to column burdens and vertical profiles at world-wide receptors. Sensitivity of these contributions to meteorological input was examined using different years of meteorology in two atmospheric simulations. Results show that Asian emission influences tracer distributions in its eastern downwind regions extending as far as Europe with major contributions in mid- and upper troposphere. On the western and eastern sides of the US, Asian contribution to annual average column burdens are 37% and 5% respectively with strong monthly variations. At an altitude of 10 km, these contributions are 75% and 25% respectively. North American emissions contribute more than 15% to annual average column burden and about 50% at 8 km altitude over the European region. Contributions from tropical African emissions are wide-spread in both the hemispheres. Differences in meteorological input cause non-uniform redistribution of tracer mass throughout the troposphere at all receptors. We also show that in model-model and model-data comparison, correlation analysis of tracer's spatial gradients provides an added measure of model's performance.

Description: We demonstrated that, in our model, non-linear interactions between planetary waves (PW) and migrating tides could generate in the upper mesosphere non-migrating tides with amplitudes comparable to those observed. The Numerical Spectral Model (NSM) we employ incorporates Hines Doppler Spread Parameterization for small-scale gravity waves (GW), which affect in numerous ways the dynamics of the mesosphere. The latitudinal (seasonal) reversals in the temperature and zonal circulation, which are largely caused by GWs (Lindzen, 198l), filter the PWs and contribute to the instabilities that generate the PWs. The PWs in turn are amplified by the momentum deposition of upward propagating GWs, as are the migrating tides. The GWs thus affect significantly the migrating tides and PWs, the building blocks of non-migrating tides. In the present paper, we demonstrate that GW filtering also contributes to the non-linear coupling between PWs and tides. Two computer experiments are presented to make this point. In one, we simply turn off the GW source to show the effect. In the second case, we demonstrate the effect by selectively suppressing the momentum source for the m = 0 non-migrating tides.

Description: We discuss results from a modeling study with our Numerical Spectral Model (NSM) that specifically deals with the non-migrating tides generated in the mesosphere. The NSM extends from the ground to the thermosphere, incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GWs), and it describes the major dynamical features of the atmosphere including the wave driven equatorial oscillations (QBO and SAO), and the seasonal variations of tides and planetary waves. Accounting solely for the excitation sources of the solar migrating tides, the NSM generates through dynamical interactions also non-migrating tides in the mesosphere that are comparable in magnitude to those observed. Large non-migrating tides are produced in the diurnal and semi-diurnal oscillations for the zonal mean (m = 0) and in the semidiurnal oscillation for m = 1. In general, significant eastward and westward propagating tides are generated for all the zonal wave numbers m = 1 to 4. To identify the cause, the NSM is run without the solar heating for the zonal mean (m = 0), and the amplitudes of the resulting non-migrating tides are then negligibly small. In this case, the planetary waves are artificially suppressed, which are generated in the NSM through instabilities. This leads to the conclusion that the non-migrating tides are generated through non-linear interactions between planetary waves and migrating tides, as Forbes et al. and Talaat and Liberman had proposed. In an accompanying paper, we present results from numerical experiments, which indicate that gravity wave filtering contributes significantly to produce the non-linear coupling that is involved.

Description: This paper presents early results from aircraft measurements made by a Low-LET Radiation Spectrometer (LoLRS), as part of a long-range effort to study the complex dynamics of the atmospheric radiation field. For this purpose, a comprehensive data base is being generated to enable a multivariable global mapping (and eventually modeling) of doses and Linear-Energy-Transfer (LET) spectra at aviation altitudes. To accomplish this, a methodical collection of data from the LoLRS (and other instruments), is planned over extended periods of time, in a manner that complements some previous isolated and sporadic measurements by other workers, with the objective to generate a detailed long-range description of the cosmic-ray induced particle environment and to study its variability and dependence on atmospheric thickness, magnetic latitude, L-shell or rigidity, space weather, solar particle events, solar cycle effects, magnetic field variation, diurnal and seasonal effects, and atmospheric weather. Analysis of initial data indicates that the dose is rising with increasing altitude and increasing magnetic latitude. Comparison of total doses with predictions is in good agreement.

Description: Enormous advances have been made in the last quarter century in all of these needed areas, covering the two essential halves of the Sun-Climate question: in what we know of solar variations and, equally important, in what we know of the climate system and of climatic changes. These research achievements allow us to examine all aspects of the question more directly and quantitatively than was ever possible before, and in the brighter light and more objective context of other known or suspected climate change mechanisms, including human-induced global greenhouse warming. Brief summaries of present status and current understanding are given below for nine facets of Sun-Climate science in which major progress has been made in recent years. At the same time it will be seen that in every instance, significant elements of uncertainty still remain, Some of the most important of these unanswered questions are considered later, in Section IV.

Description: We review research developments on satellite detection of the convection generated stresses in the Earth for seismic hazard assessment and Earth resource survey. Particular emphasis is laid upon recent progress and results of stress calculations from which the origin and evolution of the tectonic features on Earth's surface can be scientifically addressed. An important aspect of the recent research development in tectonic stresses relative to earthquakes is the implications for earthquake forecasting and prediction. We have demonstrated that earthquakes occur on the ring of fire around the Pacific in response to the tectonic stresses induced by mantle convection. We propose a systematic global assessment of the seismic hazard based on variations of tectonic stresses in the Earth as observed by satellites. This space geodynamic approach for assessing the seismic hazard is unique in that it can pinpoint the triggering stresses for large earthquakes without ambiguities of geological structures, fault geometries, and other tectonic properties. Also, it is distinct from the probabilistic seismic hazard assessment models in the literature, which are based only on extrapolations of available earthquake data.

Description: We assess the consistency (closure) between solar beam attenuation by aerosols and water vapor measured by airborne sunphotometry and derived from airborne in-situ, and ship-based lidar measurements during the April 2001 Asian Pacific Regional Aerosol Characterization Experiment (ACE-Asia). The airborne data presented here were obtained aboard the Twin Otter aircraft. Comparing aerosol extinction o(550 nm) from four different techniques shows good agreement for the vertical distribution of aerosol layers. However, the level of agreement in absolute magnitude of the derived aerosol extinction varied among the aerosol layers sampled. The sigma(550 nm) computed from airborne in-situ size distribution and composition measurements shows good agreement with airborne sunphotometry in the marine boundary layer but is considerably lower in layers dominated by dust if the particles are assumed to be spherical. The sigma(550 nm) from airborne in-situ scattering and absorption measurements are about approx. 13% lower than those obtained from airborne sunphotometry during 14 vertical profiles. Combining lidar and the airborne sunphotometer measurements reveals the prevalence of dust layers at altitudes up to 10 km with layer aerosol optical depth (from 3.5 to 10 km altitude) of approx. 0.1 to 0.2 (500 nm) and extinction-to-backscatter ratios of 59-71 sr (523 nm). The airborne sunphotometer aboard the Twin Otter reveals a relatively dry atmosphere during ACE- Asia with all water vapor columns less than 1.5 cm and water vapor densities w less than 12 g/cu m. Comparing layer water vapor amounts and w from the airborne sunphotometer to the same quantities measured with aircraft in-situ sensors leads to a high correlation (r(sup 3)=0.96) but the sunphotometer tends to underestimate w by 7%.

Description: The Magnetospheric Constellation Mission Dynamic Response and Coupling Observatory (MC-DRACO), the logical outgrowth of a sequence of STP missions, will explore plasma transport and energy conversion processes over a broad range of spatial sizes. Designed to be a meso-/macroscope for the magnetotail, it will resolve persistent controversies and yield a new understanding on which to build a predictive science of next generation magnetospheric meteorology.

Description: We present here some of the scientific objectives which can be achieved by a single seismic station on Mars, equipped with a 3 axis VBB seismometer and a 3 axis Short Period Seismometer. We assume that this station is also equipped with meteorological sensors, including infra-sound and pressure, in order to perform a complete meteorological noise correction. The science objectives are listed in order of increasing difficulty.

Description: The severe mass extinction of marine and terrestrial organisms at the end of the Permian Period (approx. 251 Ma) was accompanied by a rapid negative excursion of approx. 3 to 4 per mil in the carbon-isotope ratio of the global surface oceans and atmosphere that persisted for some 500,000 into the Early Triassic. Simulations with an ocean-atmosphere/carbon-cycle model suggest that the isotope excursion can be explained by collapse of ocean primary productivity (a Strangelove Ocean) and changes in the delivery and cycling of carbon in the ocean and on land. Model results also suggest that perturbations of the global carbon cycle resulting from the extinctions led to short-term fluctuations in atmospheric pCO2 and ocean carbonate deposition, and to a long-term (〉1 Ma) decrease in sedimentary burial of organic carbon in the Triassic. Deposition of calcium carbonate is a major sink of river-derived ocean alkalinity and for CO2 from the ocean/atmosphere system. The end of the Permian was marked by extinction of most calcium carbonate secreting organisms. Therefore, the reduction of carbonate accumulation made the oceans vulnerable to a build-up of alkalinity and related fluctuations in atmospheric CO2. Our model results suggest that an increase in ocean carbonate-ion concentration should cause increased carbonate accumulation rates in shallow-water settings. After the end-Permian extinctions, early Triassic shallow-water sediments show an abundance of abiogenic and microbial carbonates that removed CaCO3 from the ocean and may have prevented a full 'ocean-alkalinity crisis' from developing.

Description: In this presentation, we will discuss what we can expect for the future of Earth rotation modeling, and the great challenges ahead.

Description: In this session we will describe and demonstrate the interface to the PDS access tools and functions developed for the scientific community, and discuss the potential for its utilization in K-14 formal and informal settings.

Description: It has been suggested that Magnetic Holes (MHs) are caused by perpendicular proton energization due to the Ponderomotive Force (PF) associated with phase steepened Alfven waves.

Description: Magnetic component anomaly maps were made from five mapping cycles of the Mars Global Surveyor s magnetometer data. Our goal was to find and isolate positive and negative anomaly pairs which would indicate magnetization of a single source body. From these anomalies we could compute the direction of the magnetizing vector and subsequently the location of the magnetic pole existing at the time of magnetization. We found nine suitable anomaly pairs and from these we computed four North and 3 South poles with two at approximately 60 degrees north latitude. These results suggest that during the existence of the Martian main magnetic field it experienced several reversals.

Description: The Upper Atmosphere Research Satellite (UARS) was launched in September 1991 by the Space Shuttle Discovery and continues to make relevant atmospheric measurements (as of October 2002). This successful satellite has fostered a better understanding of the middle atmospheric processes, especially those important in the control of ozone. Seven of the original ten instruments aboard the UARS are still functional and six instruments regularly make measurements. The UARS is in a stable observing configuration, in spite of experiencing several anomalies over its lifetime. It is expected that the UARS will overlap the Earth Observing System (EOS) Aura satellite (scheduled launch in January 2004) for several months before the end of the UARS mission.

Description: This paper provides information on the products available at the NASA Goddard Earth Sciences (GES) Distributed Active Archive Center (DAAC) from the Upper Atmosphere Research Satellite (UARS) mission. The GES DAAC provides measurements from the primary UARS mission, which extended from launch in September 1991 through September 2001. The ten instruments aboard UARS provide measurements of atmospheric trace gas species, dynamical variables, solar irradiance input, and particle energy flux. All standard Level 3 UARS products from all ten instruments are offered free to the public and science user community. The Level 3 data are geophysical parameters, which have been transformed into a common format and equally spaced along the measurement trajectory. The UARS data have been reprocessed several times over the years following improvements to the processing algorithms. The UARS data offered from the GES DAAC are the latest versions of each instrument. The UARS data may be accessed through the GES DAAC website at

Description: A novel method of calculating the downward ozone flux across the midlatitude (30 deg-60 deg) tropopause shows the Northern Hemisphere (NH) ozone flux to be significantly larger (approximately 24%) than that calculated in the Southern Hemisphere (SH) during the year 2000. This diagnostic method makes it possible to separate dynamical aspects of transport from the seasonal cycle of ozone in the lowermost stratosphere and explain the hemispheric difference in the ozone flux. The SH total horizontal area of exchange is equal to or slightly greater than the area of exchange in the NH throughout an annual cycle. The mean changes in potential vorticity of parcels near the tropopause are also similar or slightly greater in the SH, suggesting that NH and SH downward total mass transport to the troposphere are comparable. These results imply that the greater NH ozone flux is mostly due to the amount of ozone available for exchange rather than net hemispheric dynamical differences near the tropopause level.

Description: This study examines the winter southern hemisphere vortex of 1998 using four times daily output from a data assimilation system to focus on the polar 2-day, wave number 2 component of the 4-day wave. The data assimilation system products are from a test version of the finite volume data assimilation system (fvDAS) being developed at Goddard Space Flight Center (GSFC) and include an ozone assimilation system. Results show that the polar 2-day wave dominates during July 1998 at 70 degrees. The period of the quasi 2-day wave is somewhat shorter than 2 days (about 1.7 days) during July 1998 with an average perturbation temperature amplitude for the month of over 2.5 K. The 2-day wave propagates more slowly than the zonal mean zonal wind, consistent with Rossby wave theory, and has EP flux divergence regions associated with regions of negative horizontal potential vorticity gradients, as expected from linear instability theory. Results for the assimilation-produced ozone mixing ratio show that the 2-day wave represents a major source of ozone variation in this region. The ozone wave in the assimilation system is in good agreement with the wave seen in the POAM (Polar Ozone and Aerosol Measurement) ozone observations for the same time period. Some differences with linear instability theory are noted as well as spectral peaks in the ozone field, not seen in the temperature field, that may be a consequence of advection.

Description: The influence of the uniform large-scale flow, beta effect, and vertical shear of the environmental flow on hurricane intensity is investigated in the context of the induced convective or potential vorticity asymmetries with a hydrostatic primitive equation hurricane model. In agreement with the previous studies, imposing of one of these environmental effects can substantially weaken the simulated tropical cyclones. In response t o the environmental influence, significant asymmetries develop with a structure similar to the spiral bands in real hurricanes, which are dominated by wavenumber-one components. The tendencies of the mean radial, azimuthal winds and temperature associated with the environment-induced convective asymmetries are evaluated respectively. The resulting asymmetries can effectively reduce hurricane intensity by directly producing the negative tendency of the mean tangential wind in the vicinity of the radius of maximum wind, and by weakening the mean radial circulation. The reduction effects are closely associated with the spiral structure of the induced asymmetries. The time lag observed between the imposition of the environmental influence and the resulting rise in the minimum central pressure is the time required for developing the spiral structure. This study also confirms the axisymmetrization process associated with the induced wavenumber-one components of potential vorticity asymmetries, but it exists only within the radius of maximum wind.

Description: This article, for Advances in Geophysics, is a summary of crustal deformation studies in southcentral Alaska. In 1964, southcentral Alaska was struck by the largest earthquake (moment magnitude 9.2) occurring in historical times in North America and the second largest earthquake occurring in the world during the past century. Conventional and space-based geodetic measurements have revealed a complex temporal-spatial pattern of crustal movement. Numerical models suggest that ongoing convergence between the North America and Pacific Plates, viscoelastic rebound, aseismic creep along the tectonic plate interface, and variable plate coupling all play important roles in controlling both the surface and subsurface movements. The geodetic data sets include tide-gauge observations that in some cases provide records back to the decades preceding the earthquake, leveling data that span a few decades around the earthquake, VLBI data from the late 1980s, and GPS data since the mid-1990s. Geologic data provide additional estimates of vertical movements and a chronology of large seismic events. Some of the important features that are revealed by the ensemble of studies that are reviewed in this paper include: (1) Crustal uplift in the region that subsided by up 2 m at the time of the earthquake is as much as 1 m since the earthquake. In the Turnagain Arm and Kenai Peninsula regions of southcentral Alaska, uplift rates in the immediate aftermath of the earthquake reached 150 mm/yr , but this rapid uplift decayed rapidly after the first few years following the earthquake. (2) At some other locales, notably those away the middle of the coseismic rupture zone, postseismic uplift rates were initially slower but the rates decay over a longer time interval. At Kodiak Island, for example, the uplift rates have been decreasing at a rate of about 7mm/yr per decade. At yet other locations, the uplift rates have shown little time dependence so far, but are thought not to be sustainable throughout the several hundred year recurrence time for great earthquake. The nearly 10 mm/yr uplift rate at Seldovia on the Kenai Peninsula is an example.

Description: We have processed all available DORIS data from all available satellites, except JASON over the past 10 years. Weekly solutions have been produced for stations positions coordinates, geocenter motion and scale factor stability. We present here the current accuracy presently achievable for all types of potential geodetic products.

Description: The goal of this paper is to discuss the status of the DORIS Pilot experiment, to analyze the progress made during the past two years and to propose the official creation within the IAG of an International DORIS Service (IDS).

Description: Geophysics research has been faced with a growing need for automated techniques with which to process large quantities of data. A successful tool must meet a number of requirements: it should be consistent, require minimal parameter tuning, and produce scientifically meaningful results in reasonable time. We introduce a hidden Markov model (HMM)-based method for analysis of geophysical data sets that attempts to address these issues.

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