ALBERT

Description: In March of 1994, the GSFC Stratospheric Ozone Lidar was deployed to the Network for the Detection of Stratospheric Change (NDSC) site at Lauder, NZ. This was in conjunction with a series of NASA ER-2 flights from Christchurch, NZ south to the Antarctic Circle. These flights were organized to study the chemistry of the stratosphere before, during and after the formation of the well-known 'ozone hole'. Lidar measurements were made at four different time periods corresponding to the times of the ER-2 flights. Lauder is situated nearly along the flight path as the aircraft flew south and so the lidar measurements provide a checkpoint for the ozone, aerosol and temperature instruments onboard the aircraft. Whenever the weather permitted, lidar measurements were made as near to dawn, prior to the flight, and as near to sunset, after the flight. This provided data as close to the aircraft transit time as possible. More than 70 individual lidar measurements were made, each consisting of a vertical profile of ozone, temperature, and aerosol. These were made over three different seasons and show seasonal variation. Of particular interest in the lidar data base is the wintertime stratospheric - mesospheric temperature profiles, which show large variations at the stratopause and also some significant wave activity.

Description: The HALogen Occultation Experiment (HALOE) was launched on the Upper Atmosphere Research Satellite (UARS) by the Space Shuttle Discovery at 7:11:04 EDT on September 12, 1991. After allowing for a period of outgassing, HALOE began taking routine science observations on October 11. HALOE uses the experiment approach of solar occultation and the gas filter and broad band radiometer instruments techniques to measure vertical profiles of HCl, HF, CH4, NO, NO2, H2O, O3, aerosol, and temperature versus pressure. The measurements cover a broad altitude range from the upper troposphere in some cases to the lower thermosphere in the case of nitric oxide. Latitude coverage provided by the occultation geometry ranges from 80 deg S to 80 deg N over the course of one year. The experiment has operated essentially without flaw for more than three years. Instrument stability over this time, as judged by the maximum signal change when viewing the sun exoatmospherically is less than or equal to 2 to 3%.

Description: To model radiative transfer through the atmosphere with high accuracy, one must resort to the calculation of spectral absorption coefficients on a line-by-line basis. The calculation of these coefficients is computationally expensive for three reasons: (1) thousands of spectral lines can contribute to absorption at a single frequency; (2) the tails of spectral line profiles are long (i.e., a given line can contribute to absorption over a wide range of frequencies); and (3) the sampling frequencies at which monochromatic radiances are to be calculated must be spaced sufficiently close together to resolve the thinnest lines of interest (e.g., those that arise in the stratosphere). We have developed a new algorithm to accelerate the calculation of spectral absorption coefficients while retaining high numerical accuracy.

Description: The Lidar In-Space Technology Experiment (LITE) is a backscatter lidar built by NASA Langley Research Center to fly on the Space Shuttle. The purpose of the program was to develop the engineering processes required for space lidar and to demonstrate applications of space lidar to remote sensing of the atmosphere. The instrument was flown on Discovery in September 1994. Global observations of clouds and aerosols were made between the latitudes of 57 deg N and 57 deg S during 10 days of the mission.

Description: The Sounding of the Atmosphere Using Broadband Emission Radiometry (SABER) experiment has been selected for flight on the Thermosphere-Ionosphere-Mesosphere Energetics and Dynamics (TIMED) mission expected to fly in the latter part of this decade. The primary science goal of SABER is to achieve fundamental and important advances in understanding of the energetics, chemistry, and dynamics, in the atmospheric region extending from 60 km to 180 km altitude, which has not been comprehensively observed on a global basis. This will be accomplished using the space flight proven experiment approach of broad spectral band limb emission radiometry. SABER will scan the horizon in 12 selected bands ranging from 1.27 microns to 17 microns wavelength. The observed vertical horizon emission profiles will be mathematically inverted in ground data processing to provide vertical profiles with 2 km vertical resolution, of temperature, O3, H2O, NO, NO2, CO, and CO2. SABER will also observe key emissions needed for energetics studies at 1.27 microns (O2((sup 1)delta)), 2 microns (OH(v = 7,8,9)) 1.6 microns (OH(v = 3,4,5)), 4.3 microns (CO2(v(sub 3))) 5.3 microns (NO) 9.6 microns (O3), and 15 microns (CO2(v(sub 2))). These measurements will be used to infer atomic hydrogen and atomic oxygen, the latter inferred three different ways using only SABER observations. Measurements will be made both night and day over the latitude range from the southern to northern polar regions.

Description: The purpose of the chemistry component of the model comparison is to assess to what extent differences in the formulation of chemical processes explain the variance between model results. Observed concentrations of chemical compounds are used to estimate to what degree the various models represent realistic situations. For readability, the materials for the chemistry experiment are reported in three separate sections. This section discussed the data used to evaluate the models in their simulation of the source gases and the Nitrogen compounds (NO(y)) and Chlorine compounds (Cl(y)) species.

Description: Space-based and airborne coherent Doppler lidars designed for measuring global tropospheric wind profiles in cloud-free air rely on backscatter, beta from aerosols acting as passive wind tracers. Aerosol beta distribution in the vertical can vary over as much as 5-6 orders of magnitude. Thus, the design of a wave length-specific, space-borne or airborne lidar must account for the magnitude of 8 in the region or features of interest. The SPAce Readiness Coherent Lidar Experiment under development by the National Aeronautics and Space Administration (NASA) and scheduled for launch on the Space Shuttle in 2001, will demonstrate wind measurements from space using a solid-state 2 micrometer coherent Doppler lidar. Consequently, there is a critical need to understand variability of aerosol beta at 2.1 micrometers, to evaluate signal detection under varying aerosol loading conditions. Although few direct measurements of beta at 2.1 micrometers exist, extensive datasets, including climatologies in widely-separated locations, do exist for other wavelengths based on CO2 and Nd:YAG lidars. Datasets also exist for the associated microphysical and chemical properties. An example of a multi-parametric dataset is that of the NASA GLObal Backscatter Experiment (GLOBE) in 1990 in which aerosol chemistry and size distributions were measured concurrently with multi-wavelength lidar backscatter observations. More recently, continuous-wave (CW) lidar backscatter measurements at mid-infrared wavelengths have been made during the Multicenter Airborne Coherent Atmospheric Wind Sensor (MACAWS) experiment in 1995. Using Lorenz-Mie theory, these datasets have been used to develop a method to convert lidar backscatter to the 2.1 micrometer wavelength. This paper presents comparison of modeled backscatter at wavelengths for which backscatter measurements exist including converted beta (sub 2.1).

Description: The typical fair weather electric field at the ground is between -100 and -300 V/m. At the NASA Kennedy Space Center and US Air Force Cape Canaveral Air Station (KSC) the electric field at the ground sometimes reaches -400 to -1200 V/m within an hour or two after sunrise on days that otherwise seem to be fair weather. We refer to the enhanced negative electric fields as the "sunrise enhancement." To investigate the sunrise enhancement at KSC we measured the electric field (E) in the first few hundred meters above the ground before and during several sunrise enhancements. From these E soundings we can infer the presence of charge layers and determine their thickness and charge density.

Description: Progress, significant results, and future plans regarding the following project objectives are presented: (a) Develop techniques for optimizing structural analysis of basement trends in arid regions with extremely subdued topography and/or thin aeolian cover. b) Apply results of (a) to map the southern extension of the Hamisana Shear Zone and the western extension of Nakasib Suture. c) Apply results of (b) to constrain the roles of terrane accretion and strike-slip re-organization for late Precambrian crustal evolution in NE Africa.

Description: Progress and future plans for the following objectives are presented: (1) To develop a technique to obtain values of aeolian roughness for geologic surfaces from values of surface roughness determined from calibrated L- and C-band, like- and cross-polarized, multiple incidence angle radar data from SIR-C; (2) To define the optimal combination of radar parameters from which aeolian roughness can be derived; and (3) To gain an understanding of the physical processes behind the empirical relationship.

Description: The activities and aims of the Inter-Agency Consultative Group (IACG) for space science are reported on. The principle aim is to enhance the scientific return among the members through the coordination of their fleets of current and future spacecraft. The four current campaigns are: magnetotail energy flow and nonlinear dynamics; boundaries in the collisionless plasma; solar events and their manifestations in interplanetary space and geoscience, and solar sources of heliospheric structure observed out of the ecliptic. The first of these campaigns and its implementation are reviewed.

Description: The characteristics of traveling compression regions (TCRs) in the midtail lobes are examined. Through the use of the AL index, isolated substorm events with well developed expansion phases are selected. The TCR events which feature a field compression coincident with modified Bz variations are categorized into different types, and the magnetic variations are interpreted in terms of the relative location of the point of observation to the plasmoid at the time of release and the effects of tail flaring. In order to understand the relationship between the plasmoid release time and the substorm onset time, the time difference between the different types of TCR and the substorm onset determined by Pi 2 pulsations at mid-latitude ground stations, is examined. The results suggest that the downtail release of most of the plasmoids created earthwards of -38 earth radii occurs at almost the same distance as the substorm onset.

Description: The Global Geospace Science (GGS) program's Polar satellite is reported on. The satellite aims to measure: the plasma flux in the polar magnetosphere and the geomagnetic tail; the plasma flux to and from the ionosphere, and the deposition of particle energy in the upper atmosphere. To accomplish these objectives, the satellite was placed on a 86 deg inclination, elliptical orbit whose apogee is located over the northern polar region. The spacecraft carries instruments for observing and measuring the magnetic field and charged particles as well as the imaging instruments.

Description: This section contains a number of special diagnostics that are designed to examine certain mechanisms. Section 1 reports on the method used to test the photochemical partitioning in the models. Sections 2 and 3 represent efforts to examine the model calculated production and removal rates for ozone and how the values are combined with transport rates in the models to produce the simulated ozone distributions. Sections 4 and 5 concentrate on polar processes including the dynamics aspect of vortex confinement and the chemical aspects of chlorine activation.

Description: MM II defined a series of experiments to better understand and characterize model transport and to assess the realism of this transport by comparison to observations. Measurements from aircraft, balloon, and satellite, not yet available at the time of MM I [Prather and Remsberg, 1993], provide new and stringent constraints on model transport, and address the limits of our transport modeling abilities. Simulations of the idealized tracers the age spectrum, and propagating boundary conditions, and conserved HSCT-like emissions probe the relative roles of different model transport mechanisms, while simulations of SF6 and C02 make the connection to observations. Some of the tracers are related, and transport diagnostics such as the mean age can be derived from more than one of the experiments for comparison to observations. The goals of the transport experiments are: (1) To isolate the effects of transport in models from other processes; (2) To assess model transport for realistic tracers (such as SF6 and C02) for comparison to observations; (3) To use certain idealized tracers to isolate model mechanisms and relationships to atmospheric chemical perturbations; (4) To identify strengths and weaknesses of the treatment of transport processes in the models; (5) To relate evaluated shortcomings to aspects of model formulation. The following section are included:Executive Summary, Introduction, Age Spectrum, Observation, Tropical Transport in Models, Global Mean Age in Models, Source-Transport Covariance, HSCT "ANOY" Tracer Distributions, and Summary and Conclusions.

Description: A global lightning model that includes diurnal and annual lightning variation, and total flash density versus latitude for each major land and ocean, has been used as the basis for simulating the global electric circuit charging rate. A particular objective has been to reconcile the difference in amplitude ratios [AR=(max-min)/mean] between global lightning diurnal variation (AR approx. = 0.8) and the diurnal variation of typical atmospheric potential gradient curves (AR approx. = 0.35). A constraint on the simulation is that the annual mean charging current should be about 1000 A. The global lightning model shows that negative ground flashes can contribute, at most, about 10-15% of the required current. For the purpose of the charging rate simulation, it was assumed that each ground flash contributes 5 C to the charging process. It was necessary to assume that all electrified clouds contribute to charging by means other than lightning, that the total flash rate can serve as an indirect indicator of the rate of charge transfer, and that oceanic electrified clouds contribute to charging even though they are relatively inefficient in producing lightning. It was also found necessary to add a diurnally invariant charging current component. By trial and error it was found that charging rate diurnal variation curves in Universal time (UT) could be produced with amplitude ratios and general shapes similar to those of the potential gradient diurnal variation curves measured over ocean and arctic regions during voyages of the Carnegie Institute research vessels.

Description: Observations of ionospheric convection flows at a range of local times during the various phases of the substorm cycle are reported on with the aim of investigating the convection behavior during a range of times and phases. The ionospheric flow observations are from the EISCAT and DMSP satellites. The substorm phases are identified from energetic particle measurements from geosynchronous satellites. The growth phase convection indicates an initial expansion of the polar cap. There is an unexplained poleward motion of the flow reversal boundary (FRB). It is concluded that this motion does not necessarily provide a true representation of the balance between reconnection at the dayside and in the tail. The expansion phase flows do not show any evidence for tail reconnection until late in the phase. The convection during the recovery phase is indicative of tail reconnection as there is evidence that there is only a lobe cell driving convection on the dayside.

Description: GOME radiance, irradiance, and ozone products were validated by NASA Goddard Space Flight Center through three tasks which included, pre-launch calibration comparisons with SBUV and TOMS radiometric standards, validation of GOME Level-1 irradiance and radiance and Level 2 total ozone data products using SBUV/2 and TOMS algorithms and data, and studies of GOME data using the Goddard radiative transfer code. The prelaunch calibration using the NASA large aperture integrating sphere was checked against that provided by TPD. Agreement in the calibration constants, derived in air, between the Goddard and TPD system were better than 3%. Validation of Level-1 irradiance data included comparison of GOME and SSBUV and the UARS solar irradiances measurements. Large wavelength dependent differences, as high as 10%, were noted between GOME and the US instruments. This discrepancy has now been attributed to radiometric sensitivity changes experienced by GOME when operating in a vacuum. GOME Earth radiance data were then compared to the NOAA-14 SBUV/2 radiances. These results show that between 340 and 400 nm the differences in GOME and SBUV/2 data are less than 5% with some wavelength dependence. At wavelengths shorter than 300 nm, differences are of the order of 10% or more where the GOME radiances are larger. To test GOME DOAS retrieved total ozone values, these values were compared with ozone amounts retrieved using GOME radiances in the TOMS version-7 algorithm. The differences showed a solar zenith angle dependence ranging from 0 to 10% where the TOMS algorithm values were higher. GOME radiances below 300 nm were further validated by selecting radiances at wavelengths normally used by SBUV and processing them through the SBUV ozone profile algorithm and then compared to climatological values. The GOME ozone profiles ranged from 10-30% lower over altitude compared to climatological values. This is consistent with the offsets detected in the SBUV/2 radiance comparisons at wavelengths shorter than 300 nm.

Description: The importance of the OH radical as an intermediate in many combustion reactions and in atmospheric photochemistry has led many researchers to use it as a diagnostic tool in these processes. The amount of data that has been acquired over the years for this radical is quite considerable. However, the quenching rate of OH with water molecules as a function of temperature and the rotational level of the excited state is not very well understood. The motivation of the studies undertaken is to bridge the gap between the low temperature measurements and the high temperature ones reported in the literature. The technique generally employed in these diagnostics is laser-induced fluorescence (LIF), through which rotational state selective excitation of the radical is possible. Furthermore, in a combustion medium, water is produced in abundance so that knowledge of the quenching rate of OH due to water molecules plays a crucial role in interpreting the data. In general, the precursor to an understanding of the collisional quenching rates of OH involves a characterization of the mode in which the radical is produced; the resulting rotational and translational distribution, followed by a measurement of the OH temperature; and ultimately obtaining the rate constants from the pressure dependence of the fluorescence signal. The experimental implementation of these measurements therefore involved, as a first step, the production of the OH radicals in a microwave discharge cell using water vapor as the source, wherein a hydrogen atom is abstracted from H2O. The second step involved the absorption of photons from the frequency-doubled output of a pulsed amplified, single-frequency cw ring dye laser. By tuning the laser to the peak of the transition and observing the fluorescence decay after the laser pulse, the lifetime of the OH in a particular rotational electronic state was determined (tau = 1.4 microseconds for Q(sub 1)(3)). Knowledge of this parameter led to a determination of the quenching rate. By varying the water vapor pressure in the cell and measuring the lifetime as a function of pressure a linear plot of the quenching rate as a function of pressure was obtained. Using this plot, the quenching cross section was deduced. It has therefore been possible to measure the local translational temperature and the quenching cross section with one laser system.

Description: The application of remote sensing techniques to the analysis of the dynamics of storms and substorm processes is discussed. The Imager for Magnetopause-to-Aurora Global Exploration (IMAGE) mission, its mission and instruments are presented. The following are discussed: neutral atom imaging; radio plasma techniques; photon imaging, and substorm observations.

Description: The energy and magnetic flux budgets of the magnetotail plasma sheet during substorm expansion are investigated. The possible mechanisms that change the energy content of the closed field line region which contains all the major dissipation mechanisms of relevance during substorms, are considered. The compression of the plasma sheet mechanism and the diffusion mechanism are considered and excluded. It is concluded that the magnetic reconnection mechanism can accomplish the required transport. Data-based empirical magnetic field models are used to investigate the magnetic flux transport required to account for the observed magnetic field dipolarizations in the inner magnetosphere. It is found that the magnetic flux permeating the current sheet is typically insufficient to supply the required magnetic flux. It is concluded that no major substorm-type magnetospheric reconfiguration is possible in the absence of magnetic reconnection.

Description: Hybrid simulations are used to investigate the formation of a thin current sheet inside the plasma sheet of a magnetotail-like configuration. The initial equilibrium is subjected to a driving electric field which is qualitatively similar to what would be expected from solar wind driving. As a result, a new current sheet with the thickness of approximately the ion inertial length is formed. The current density inside the current sheet region is supplied largely by the electrons. Ion acceleration in the cross-tail direction is absent as the driving electric field fails to penetrate into the equatorial region.

Description: Ion acceleration and flux increase associated with substorm energetic particle injections are investigated on the basis of geosynchronous observations and test proton orbits in the dynamic fields of a three dimensional magnetohydrodynamic simulation of neutral line formation and dipolarization in the magnetotail. The energetic particle flux changes obtained from the test particle orbits agree well with observations that demonstrate rapid ion flux increases at energies of above 20 keV. The injection region inferred from the test particles has a sharp earthward boundary and a sharp ragged tailward boundary. The earthward portion of the enhanced ion flux can be traced to the enhanced cross-tail electric field associated with the near-earth x-type neutral line. Due to the rapid earthward motion of accelerated ions away from the neutral line, this boundary is displaced earthward to where the energetic ions become more adiabatic in the stronger dipolar field.

Description: The theoretical modeling of atmospheric spectra is important for a number of different applications: for instance, in the determination of minor atmospheric constituents such as ozone, carbon dioxide, CFC's etc.; in monitoring the temperature profile for climate studies; and in measuring the incoming and outgoing radiation to input into global climate models. In order to accomplish the above mentioned goal, one needs to know the spectral parameters characterizing the individual spectral lines (frequency, width, strength, and shape) as well as the physical parameters of the atmosphere (temperature, abundances, and pressure). When all these parameters are known, it is usually assumed that the resultant spectra and concomitant absorption coefficient can then be calculated by a superposition of individual profiles of appropriate frequency, strength and shape. However, this is not true if the lines are 'coupled'. Line coupling is a subtle effect that takes place when lines of a particular molecule overlap in frequency. In this case when the initial states and the final states of two transitions are connected by collisions, there is a quantum interference resulting in perturbed shapes. In general, this results in the narrowing of Q-branches (those in which the rotational quantum number does not change), and vibration-rotational R- and P branches (those in which the rotational quantum number changes by +/- 1), and in the spectral region beyond band heads (regions where the spectral lines pile up due to centrifugal distortion). Because these features and spectral regions are often those of interest in the determination of the abundances and pressure-temperature profiles, one must take this effect into account in atmospheric models.

Description: Gas filter correlation radiometer (GFCR) is a passive remote sensing technique used in a variety of atmospheric measurements. In recent years, a nonmechanical optical switching GFCR has been invented and developed at NASA Langley Research Center. The use of a polarization modulator, in conjunction with a polarization beamsplitter, enables rapid optical switching without mechanically moving parts. In comparison with the conventional GFCR, which involves mechanical chopping or switching between two optical paths, the nonmechanical GFCR possesses some very attractive advantages such as fast sampling rate, high reliability, low weight, and long operational life time. In a recent study, we have developed a new GFCR configuration and have fabricated a compact, nonmechanical breadboard instrument. Using this instrument, we have carried out atmospheric methane measurements in the 2.3 micron region. Measurement results are compared with theoretical predictions using the HITRAN database.

Description: Knowledge of the global scale distribution of atmospheric ozone and its temporal variability can be achieved using a satellite-based nadir-viewing device making high spectral resolution measurements with high signal-to-noise ratios. This would enable observation in the pressure-broadened wings of strong O3 lines while minimizing the impact of undesirable signal contributions associated with, for example, the terrestrial surface and interfering species. The Fabry-Perot interferometer (FPI) provides high spectral resolution and high throughput capabilities that are essential for this measurement task. The periodic nature of the Fabry-Perot instrument function can be advantageous when observation of periodic spectra is desired. However, for most applications, additional optical elements are necessary to reduce the effect of unwanted passbands. This is frequently accomplished using additional Fabry-Perot etalons in a series configuration in conjunction with a bandpass filter. This paper discusses a Fabry-Perot interferometer conceptual instrument design to achieve tropospheric and total ozone monitoring capability from a satellite-based nadir-viewing geometry. The design involves a double-etalon fixed-gap series configuration FPI along with an ultra-narrow bandpass filter to achieve single-order operation with an overall spectral resolution of approximately .068 cm(exp -1). The impact of inter-etalon reflections has been reduced to acceptable levels by placement of a slightly attenuating medium in between the etalons. A passive device is selected for low power consumption, and continuous day/night coverage, independent of solar zenith angle, is enabled by observing within the strong 9.6 micron ozone infrared band. The IR-FPI detection will be performed through implementation of the new Circle to Line Interferometer Optical (CLIO) system, developed by researchers at the Space Physics Research Laboratory (SPRL) of the University of Michigan, to accomplish focal plane fringe detection; the CLIO system converts the circular interferometric fringes into a linear pattern which then can be detected by conventional linear array detectors. A multiplex signal advantage is achievable as all necessary frequencies can be measured simultaneously using a multichannel configuration. Through proper selection of channel spectral regions, the FPI optimized for tropospheric O3 measurements can simultaneously observe a stratospheric component and thus the total O3 column abundance.

Description: A practical procedure for the retrieval of ozone vertical profiles from ground-based high resolution Fourier transform infrared solar spectra has been developed. The analysis is based on a multilayer line-by-line forward model and a semi-empirical version of the optimal estimation inversion method of Rodgers. The 1002.6-1003.2 cm(exp -1) spectral interval has been selected for the analysis on the basis of synthetic spectrum calculations. This interval contains numerous ozone lines covering a range of intensities and providing retrieval sensitivity from ground level to about 35 km. Characterization of the method and an error analysis have been performed. For a spectral resolution of 0.05-0.01 cm(exp -1) and a signal-to-noise ratio greater than or equal to 100 the retrieval is stable with a vertical resolution of approximately 5 km attainable near the surface degrading to approximately 10 km in the stratosphere. Synthetic spectra studies show that the a priori profile and weak constraints selected for the retrievals introduce no significant biases for a wide range of ozone profiles.

Description: Estimates of the effect of pulse stretching on satellite laser altimetry in particular the Geoscience Laser Altimeter System (GLAS), by cloud multiple scattering were made from an analytical method and from Monte Carlo calculations. The path delay of the return pulse was found to be largest for low-level clouds with particle radii (3-20 microns). The magnitude of the path delay was affected by several factors including cloud height, cloud optical depth, cloud particle size, particle shape, and receiver field of view. Polar aerosol and Rayleigh scattering usually added less than 1 cm to the overall path delay. Path delay estimates for all cloud conditions would be less if a simple Gaussian fit of the return pulse peak were used to measure the pulse's centroid. However, care must be taken in determining the centroid as factors such as pulse width, surface slope and the fitting method used will affect the estimate. A planned application for laser altimetry is high precision monitoring of the height change of polar ice sheets. In the absence of a correction algorithm, the required GLAS altimetry accuracies will not be achieved unless atmospheric channel information is used to remove profiles that are likely to be heavily contaminated by multiple scattering.

Description: The zonal mean eddy heat flux is directly proportional to the wave activity that propagates from the troposphere into the stratosphere. This quantity is a simple eddy diagnostic which is easily calculated from conventional meteorological analyses. Because this "wave driving" of the stratosphere has a strong impact on the stratospheric temperature, it is necessary to compare the impact of the flux with respect to stratospheric radiative changes caused by greenhouse gas changes. Hence, we must understand the precision and accuracy of the heat flux derived from our global meteorological analyses. Herein, we quantify the stratospheric heat flux using five different meteorological analyses, and show that there are 30% differences between these analyses during the disturbed conditions of the northern hemisphere winter. Such large differences result from the planetary differences in the stationary temperature and meridional wind fields. In contrast, planetary transient waves show excellent agreement amongst these five analyses, and this transient heat flux appears to have a long term downward trend.

Description: We examine concurrent upper tropospheric measurements of CN (diameter greater than 4 nm). NO, and NO(Y) during the SONEX Experiment over the North Atlantic (Oct.-Nov., 1997). Elevated CN and NO(Y) concentrations observed in the upper troposphere are attributed largely to enhancements in convective outflows. We estimate that less than 7% of observed high-CN plumes (greater than 10000 /cc) may be attributed to aircraft emissions. Dilution of high-CN convective and aircraft plumes appears to be much more rapid than losses of NO(X) and CN by oxidation and coagulation, respectively, and accounts for much of observed CN concentrations. When taking into account of different time scales against dilution for observable aircraft and convective high-CN plumes (estimated to be 1:4), the contribution by aircraft emissions to CN concentrations is significant, about 20% of the convective source. We find no evidence that particle formation in convective plumes is limited by OH oxidation of SO2.

Description: This paper points out that the nutation terms of Zhu & Groten due to the tidal potential of degree 3 are erroneous. Correct values are deduced here and they coincide very well with those given in Kinoshita & Souchay. These errors explain the discrepancies between the evaluation of the theories of Zhu & Groten and Kinoshita & Souchay, which can reach values up to 165 micro arc second. Also the two leading nutation terms due to the tidal potential of degree 4 are given. Finally, some of the computer programs (NUTC.F, KSV_1994.F) for evaluating the,13 nutation terms are also erroneous.

Description: This study examines a unique data set returned by IMP8 and Geotail on January 29, 1995 during a substorm which resulted in the ejection of a plasmoid. The two spacecraft (s/c) were situated in the north lobe of the tail and both observed a traveling compression region (TCR). From single s/c observations only the length of the plasmoid in X and an estimate of its height in Z can be determined. However, we show that dual s/c measurements of TCRs can be used to model all three dimensions of the underlying plasmoid and to estimate of its rate of expansion or contraction. For this event plasmoid dimensions of Delta(X) approximates 18, Delta(Y) approximates 30, and Delta(Z) approximates 10 R(sub e) are inferred from the IMP8 and Geotail lobe magnetic field measurements. The earthward end of the plasmoid was inferred to be near the mean location of the near-earth neutral line, X approximates -26 R(sub e). Its center was underneath IMP 8 at X approximates -34 R(sub e) and its tailward end appeared to be near X approximates -44 R(sub e). Furthermore, a factor of approximately 2 increase in the amplitude of the TCR occurred in the 1.5 min it took to move from IMP 8 to Geotail. Modeled using conservation of the magnetic flux, this increase in lobe compression implies that the underlying plasmoid was expanding at a rate of approximately 140 km/s. Such an expansion is comparable to recently reported V(sub y) speeds in "young" plasmoids in this region of the tail. Finally, the Geotail measurements indicate that a reconfiguration of the lobe magnetic field closely followed the ejection of the plasmoid which moved magnetic flux tubes into the wake behind the plasmoid where they would convect into the near-earth neutral line and reconnect.

Description: Prior to 1991, major warmings (defined by increasing zonal mean temperatures and zonal mean easterly winds from 60degN to the pole at 10 hPa) typically occurred approximately once every two Arctic winters; a major warming in mid-Dec. 1998 was the first since Feb. 1991. The Dec. 1998 warming was also the second earliest on record. The earliest, and the only other major warming on record before the end of Dec. was in early Dec 1987; prior to that, the earliest was in late Dec./early Jan. 1984-85. The 1984-85 and 1987 warmings resulted in the warmest and weakest lower stratospheric polar vortices in the 20 years before 1998-99. Fig. 1 compares temperatures and vortex strength in 1998-99 with those in the previous 20 years, using the US National Center for Environmental Prediction (NCEP) record; 1987-88 and 1984-85 are also highlighted. The Dec. 1998 warming had a more pronounced effect on mid-stratospheric temperatures than the Dec. 1987 warming (Fig. 1a), although smaller than that of warmings later in winter (e.g., 1984-85). 10-hPa temperatures fell well below average again in late Jan. 1999 and remained unusually low until an early final warming began in late Feb. 840 K PV gradients (Fig. 1c) set a record minimum in Jan. 1999, but were near average in Feb before the final warming. The effect of the Dec. 1998 warming on lower stratospheric temperatures was comparable to that of other major warmings; there was a brief period of record-high minimum 46-hPa temperatures in early Jan 1999 (Fig. 1b), and temperatures then fell to near average for a short period in mid-Feb. Lower stratospheric PV gradients were the weakest on record during the 1998-99 winter (Fig. 1d). The evolution of the vortex and minimum temperatures during 1998-99 was remarkably similar to that during 1987-88, the only previous year when a major warming was observed before the end of Dec.

Description: The study of the magnetosphere is a 400 year old science that began with the publication by Gilbert, in 1600, of his hypotheses that the Earth was a giant magnet. Since then we have learned many things about the magnetosphere, particularly in the last 40 years of the space age, but we still have many unanswered questions. In spite of the many thousands of observations of this system we still lack a global understanding of how it works. This is due to its large size and tenuous nature that mean that any measurement made of the fields or particles involved only give one a knowledge of the local conditions at a given time. To gain a global perspective through such observations would require the simultaneous operation of thousands of satellites spread throughout the magnetospheric system in addition to observations made on the ground. Such a program would be impractical at least from financial considerations. What is needed for the advancement of magnetospheric physics is to develop the same capabilities that astrophysicists, solar physicists and meteorologists have been using for years --- the ability to stand back from the object under study and see it in its entirety. The challenge for doing this for the magnetosphere is that the particle densities are very low and the material is, for the most part, not luminous. In the last 25 years several ideas have been proposed that would allow at least the imaging of certain portions of the magnetosphere. These include imaging of the plasmasphere through the resonant scattering of solar 304 A from He+ ions, imaging of various hot plasma populations (i.e. the ring current, plasmasheet, upflowing ionospheric ions, etc.) from the neutral atoms that result when ions of these populations charge exchange with the hydrogen geocorona, and imaging the aurora at various wavelengths in the far ultraviolet. In addition, a novel technique for probing various boundaries in the magnetosphere by bouncing low frequency radio waves off of them has been extensively studied. Such a technique is analogous to the way the under water world can be probed with sonar. About five years ago NASA convened a science working group to study the possibility of flying a magnetospheric imaging mission. This resulted in a number of proposals for such a mission, one of which was selected to be the first MIDEX mission, to be launched in early 2000. The mission is called IMAGE (Imager for Magnetopause to Aurora Global Exploration) and its P.I. is J. Burch at SwRI. The IMAGE spacecraft will carry imagers to view the plasmasphere, aurora, ring current, inner plasmasheet, and upflowing ionospheric ions as well as a radio sounder to probe the location, shape and dynamics of the magnetopause, plasmapause, etc. Between its selection last April and the non advocacy mission review, which takes place next spring, the IMAGE teams needs to further refine the design of the mission and its instruments. The theory and modeling (T&M) subgroup of this team has the task of demonstrating what kind of images the instruments on IMAGE will see as well as showing that useful scientific information can be extracted from such images. As a central element to the efforts of the T&M subgroup we have decided to simulate and create synthetic images for the magnetic cloud event of October, 1995. In this event a large cloud, with high plasma densities and strong magnetic fields, ejected from the sun collided with the earth's magnetosphere triggering a three day period of intense magnetic storms and substorms. This event was observed from a number of different spacecraft and on the ground so we have a good data set to work with. In our work we will place the IMAGE spacecraft in the magnetosphere on its proposed orbit, with its proposed instruments, to see what it would see had it been there. Existing models of the plasmasphere, ring current and magnetopause will be run for this event to give the structures for the imaging instruments. There are several models which are lacking and which need to be developed. These include a model for the cusp, the inner plasmasheet and the upflowing ions. My task this summer was to develop the upflowing ion model and use it to create synthetic images.

Description: The objective of my NASA summer research project was to create a homepage to describe and present results from the NASA Global Tropospheric Experiment (GTE). The GTE is a major component of NASA's Tropospheric Chemistry Program and is managed in the Atmospheric Studies Branch, Atmospheric Sciences Division at the NASA Langley Research Center.

Description: The Clouds and The Earth's Radiant Energy System (CERES) simulator will allow flight operational familiarity with the CERES instrument prior to launch. It will provide a CERES instrument simulation facility for NASA Langley Research Center. NASA Goddard Space Flight Center and TRW. One of the objectives of building this simulator would be for use as a testbed for functionality checking of atypical memory uploads and for anomaly investigation. For instance, instrument malfunction due to memory damage requires troubleshooting on a simulator to determine the nature of the problem and to find a solution.

Description: Two different methods for retrieving Upper Tropospheric Humidities (UTH) from the TOVS (TIROS Operational Vertical Sounder) instruments aboard NOAA polar orbiting satellites are presented and compared. The first one, from the Environmental Technology Laboratory, computed by J. Bates and D. Jackson (hereafter BJ method), estimates UTH from a simplified radiative transfer analysis of the upper tropospheric infrared water vapor channel at wavelength measured by HIRS (6.3 micrometer). The second one results from a neural network analysis of the TOVS (HIRS and MSU) data developed at, the Laboratoire de Meteorologie Dynamique (hereafter the 3I (Improved Initialization Inversion) method). Although the two methods give very similar retrievals in temperate regions (30-60 N and S), an absolute bias up to 16% appears in the convective zone of the tropics. The two datasets have also been compared with UTH retrievals from infrared radiance measurements in the 6.3 micrometer channel from the geostationary satellite METEOSAT (hereafter MET method). The METEOSAT retrievals are systematically drier than the TOVS-based results by an absolute bias between 5 and 25%. Despite the biases, the spatial and temporal correlations are very good. The purpose of this study is to explain the deviations observed between the three datasets. The sensitivity of UTH to air temperature and humidity profiles is analysed as are the clouds effects. Overall, the comparison of the three retrievals gives an assessment of the current uncertainties in water vapor amounts in the upper troposphere as determined from NOAA and METEOSAT satellites.

Description: Carbonados are porous polycrystalline (with crystal sizes up to 100 micrometer) diamonds. Carbonado is found only in alluvial deposits in Bahia, Brazil and in the Central African Republic (CAR). Alluvial deposit host is 1.5 Ga while the carbonados are between 2.6 - 3.8 Ga. The process of fusing the carbonado microcrystals together is not fully understood, partly due to fact that the origin of these carbonado, is not known. Several modes of origins are proposed for carbonado. First, a crustal origin. Carbonados have a light carbon and helium isotopic signature. They contain an enrichment of the rare-earth elements (REE). Carbonados have tightly trapped atmospheric noble gases and contain an evidence of high He content despite the carbonado expected depletion of He at mantle temperatures. Carbonados have high porosity incompatible with high pressure mantle conditions. Second, a mantle origin is proposed based on similar REE pattern to kimberlites. The presence of nitrogen platelet (by IR spectra) indicates high temperature origin and syngenetic inclusions of rutile, ilmenite, and magnetite indicates high pressure and high temperature conditions consistent with mantle origin as well. Third, it is proposed that carbonado diamonds are a result of early impacts into crustal rocks. This is indicated by the rare and controversial occurrence of high pressure diamond polymorph, londsdaleite, frequently found in diamonds from meteorite impact sites, and by observation of planar deformation features, possibly indicating shock events. Finally, it is suggested that carbonados have an extraterrestrial origin, as indicated by a long term annealing based on observation of a zero-phonon line, attributed to paired nitrogen atoms in association with a vacancy.

Description: Vertical electron-density (N (sub e)) profiles, deduced from newly-available ISIS-II digital ionospheric topside-sounder data, are used to investigate the "polar-hole" region within the winter, nighttime polar cap ionosphere during solar minimum. The hole region is located around 0200 MLT near the poleward side of the auroral oval. Earlier investigations had revealed very low N (sub e) values in this region (down to 200/cu cm near 300 km). In the present study, such low N, values (approx. 100/cu cm) were only found near the ISIS (International Satellite for Ionospheric Study)-II altitude of 1400 km. The peak ionospheric concentration below the spacecraft remained fairly constant (approx. 10 (exp 5)/cu cm across the hole region but the altitude of the peak dropped dramatically. This peak dropped, surprisingly, to the vicinity of 100 km. These observations suggest that the earlier satellite in situ measurements, interpreted as deep holes in the ionospheric F-region concentration, could have been made during conditions of an extreme decrease in the altitude of the ionospheric N (sub e) peak. The observations, in combination with other data, indicate that the absence of an F-layer peak may be a frequent occurrence at high latitudes.

Description: We have examined the sea level height tide records at seven tide gauge sites in the region of southcentral Alaska that were affected by the 1964 Prince William Sound earthquake to determine the history of crustal uplift subsequent to the earthquake. There is considerable variation in the behavior depending on the location of the site relative to the 1964 rupture. At Seward, on the eastern side of the Kenai Peninsula we find a slow uplift that is consistent with elastic strain accumulation while at Seldovia and Nikiski on the western side of the Kenai we find a persistent rapid uplift of about 1 cm/yr that most likely represents a long term transient response to the earthquake, but which cannot be sustained over the expected recurrence interval for a great earthquake of several hundred years. Further to the southwest, at Kodiak, we find evidence that the rate of uplift, which is still several mm/yr, has slowed significantly over the past three and a half decades. To the east of the Kenai Peninsula we find subsidence at Cordova and an uncertain behavior at Valdez. At both of these sites there is a mathematically significant time-dependence to the uplift behavior, but the data confirming this time dependence are not as convincing as at Kodiak. At Anchorage, to the north there is little evidence of vertical motion since the earthquake. We compare these long term tide gauge records to recent GPS observations. In general there is reasonable consistency except at Anchorage and Cordova where the GPS measurement indicate somewhat more rapid uplift and subsidence, respectively.

Description: The geographic limits of the South Atlantic Anomaly, as defined by radiation damage, are compared to contours of geomagnetic total field intensity, as defined by the 1995 IGRF, for the present and recent past. The most likely secular variation of the geomagnetic field is estimated and used to extrapolate the geomagnetic field to the year 2100. This indicates that radiation damage to spacecraft and humans in space will likely increase and to cover a much larger geographic area.

Description: We have adopted the transport scenarios used in Part 1 to examine the sensitivity of stratospheric aircraft perturbations to transport changes in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric residence time and reduced the magnitude of the negative perturbation response in total ozone. Increasing the stratospheric K(sub yy) increased the residence time and enhanced the global scale negative total ozone response. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and results in a significantly weaker perturbation response, relative to the base case, throughout the stratosphere. We found a relatively minor model perturbation response sensitivity to the magnitude of K(sub yy) in the tropical stratosphere, and only a very small sensitivity to the magnitude of the horizontal mixing across the tropopause and to the strength of the mesospheric gravity wave drag and diffusion. These transport simulations also revealed a generally strong correlation between passive NO(sub y) accumulation and age of air throughout the stratosphere, such that faster transport rates resulted in a younger mean age and a smaller NO(y) mass accumulation. However, specific variations in K(sub yy) and mesospheric gravity wave strength exhibited very little NO(sub y)-age correlation in the lower stratosphere, similar to 3-D model simulations performed in the recent NASA "Models and Measurements" II analysis. The base model transport, which gives the most favorable overall comparison with inert tracer observations, simulated a global/annual mean total ozone response of -0.59%, with only a slightly larger response in the northern compared to the southern hemisphere. For transport scenarios which gave tracer simulations within some agreement with measurements, the annual/globally averaged total ozone response ranged from -0.45% to -0.70%. Our previous 1995 model exhibited overly fast transport rates, resulting in a global/annually averaged perturbation total ozone response of -0.25%, which is significantly weaker compared to the 1999 model. This illustrates how transport deficiencies can bias model simulations of stratospheric aircraft.

Description: According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Description: The subject of neotectonics, covering the structures and structural activity of the last 5 million years (i.e., post-Miocene) is a well-recognized field, including "active tectonics," focussed on the last 500,000 years in a 1986 National Research Council report of that title. However, there is a cartographic gap between tectonic maps, generally showing all features regardless of age, and maps of current seismic or volcanic activity. We have compiled a map intended to bridge this gap, using modern data bases and computer-aided cartographic techniques. The maps presented here are conceptually descended from an earlier map showing tectonic and volcanic activity of the last one million years. Drawn by hand with the National Geographic Society's 1975 "The Physical World" map as a base, the 1981 map in various revisions has been widely reproduced in textbooks and various technical publications. However, two decades of progress call for a completely new map that can take advantage of new knowledge and cartographic techniques. The digital tectonic activity map (DTM), presented in shaded relief (Fig. 1) and schematic (Fig. 2) versions, is the result. The DTM is intended to show tectonism and volcanism of the last one million years, a period long enough to be representative of global activity, but short enough that features such as fault scarps and volcanos are still geomorphically recognizable. Data Sources and Cartographic Methods The DTM is based on a wide range of sources, summarized in Table 1. The most important is the digital elevation model, used to construct a shaded relief map. The bathymetry is largely from satellite altimetry, specifically the marine gravity compilations by Smith and Sandwell (1996). The shaded relief map was designed to match the new National Geographic Society world physical map (1992), although drawn independently, from the digital elevation model. The Robinson Projection is used instead of the earlier Van der Grinten one. Although neither conformal nor equal-area, the Robinson Projection provides a reasonable compromise and retains useful detail at high latitudes.

Description: There are procedures and methods for verification of coding algebra and for validations of models and calculations that are in use in the aerospace computational fluid dynamics (CFD) community. These methods would be efficacious if used by the glacier dynamics modelling community. This paper is a presentation of some of those methods, and how they might be applied to uncertainty management supporting code verification and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modelling are discussed. After establishing sources of uncertainty and methods for code verification, the paper looks at a representative sampling of verification and validation efforts that are underway in the glacier modelling community, and establishes a context for these within overall solution quality assessment. Finally, an information architecture and interactive interface is introduced and advocated. This Integrated Cryospheric Exploration (ICE) Environment is proposed for exploring and managing sources of uncertainty in glacier modelling codes and methods, and for supporting scientific numerical exploration and verification. The details and functionality of this Environment are described based on modifications of a system already developed for CFD modelling and analysis.

Description: The behavior of mesospheric odd nitrogen species during and following relativistic and diffuse auroral precipitation events is simulated, Below 75 km nitric oxide is enhanced in proportion to the ion pair production function associated with the electron precipitation and the length of the event. Nitrogen dioxide and nitric acid are also enhanced. At 65 km the percentage of odd nitrogen for N is 0.1%, HNO3 is 1.6%, NO2 is 15%, and NO is 83.3%. Between 75 and 85 km NO is depleted during particle events due to the faster destruction of NO by N relative to the production of NO by N reacting with O2. Recovery of NO depends on transport from the lower thermosphere, where NO is produced in abundant amounts during particle events.

Description: Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Description: Magnetic remanence of crustal rocks can reside in three common rock-forming magnetic minerals: magnetite, pyrrhotite, and hematite. Thermoremanent magnetization (TRM) of magnetite and pyrrhotite is carried mostly by single domain (SD) grains. The TRM of hematite grains, however, is carried mostly by multidomain (NM) grains. This characteristic is illustrated by TRM acquisition curves for hematite of variable grainsizes. The transition between truly NM behavior and tendency towards SD behavior his been established between hematite grainsizes of 0. 1 and 0.05 mm. Coarse grainsize of lower crustal rocks and the large sensitivity of MD hematite grains to acquire TRM indicates that hematite could be a significant contributor to long-wavelength magnetic anomalies.

Description: We present an extension for the 2D (zonal mean) version of our Numerical Spectral Mode (NSM) that incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations and the semi-annual and quasi-biennial oscillations (SAO and QBO). Our earlier model reproduced the salient features of the mean zonal circulation in the middle atmosphere, including the QBO extension into the upper mesosphere inferred from UARS measurements. In the present model we incorporate also tropospheric heating to reproduce the upwelling at equatorial latitudes associated with the Brewer-Dobson circulation that affects significantly the dynamics of the stratosphere as Dunkerton had pointed out. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. The QBO period in the model is 30 months (mo), which is conducive to synchronize this oscillation with the seasonal cycle of solar forcing. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. Quadratic non-linearities generate interseasonal variations to produce a complicated pattern of variability associated with the QBO. The computed temperature amplitudes for the SAO and QBO are in substantial agreement with observations at equatorial and extratropical latitudes. At high latitudes, however, the observed QBO amplitudes are significantly larger, which may be a signature of propagating planetary waves not included in the present model. The assumption of hydrostatic equilibrium not being imposed, we find that the effects from the vertical Coriolis force associated with the equatorial oscillations are large for the vertical winds and significant for the temperature variations even outside the tropics but are relatively small for the zonal winds.

Description: This paper reviews observations of stratospheric temperatures that have been made over a period of several decades. Those observed temperatures have been used to assess variations and trends in stratospheric temperatures. A wide range of observation datasets have been used, comprising measurements by radiosonde (1940s to the present), satellite (1979 - present), lidar (1979 - present) and rocketsonde (periods varying with location, but most terminating by about the mid-1990s). In addition, trends have also been assessed from meteorological analyses, based on radiosonde and/or satellite data, and products based on assimilating observations into a general circulation model. Radiosonde and satellite data indicate a cooling trend of the annual-mean lower stratosphere since about 1980. Over the period 1979-1994, the trend is 0.6K/decade. For the period prior to 1980, the radiosonde data exhibit a substantially weaker long-term cooling trend. In the northern hemisphere, the cooling trend is about 0.75K/decade in the lower stratosphere, with a reduction in the cooling in mid-stratosphere (near 35 km), and increased cooling in the upper stratosphere (approximately 2 K per decade at 50 km). Model simulations indicate that the depletion of lower stratospheric ozone is the dominant factor in the observed lower stratospheric cooling. In the middle and upper stratosphere both the well-mixed greenhouse gases (such as CO) and ozone changes contribute in an important manner to the cooling.

Description: This paper presents a detailed characterization of seasonal and interannual variability in tropical tropospheric column ozone (TCO). TCO time series are derived from 20 years (1979-1998) of total ozone mapping spectrometer (TOMS) data using the convective cloud differential (CCD) method. Our study identifies three regions in the tropics with distinctly different zonal characteristics related to seasonal and interannual variability. These three regions are the eastern Pacific, Atlantic, and western Pacific. Results show that in both the eastern and western Pacific seasonal-cycle variability of northern hemisphere (NH) TCO exhibits maximum amount during NH spring whereas largest amount in southern hemisphere (SH) TCO occurs during SH spring. In the Atlantic, maximum TCO in both hemispheres occurs in SH spring. These seasonal cycles are shown to be comparable to seasonal cycles present in ground-based ozonesonde measurements. Interannual variability in the Atlantic region indicates a quasi-biennial oscillation (QBO) signal that is out of phase with the QBO present in stratospheric column ozone (SCO). This is consistent with high pollution and high concentrations of mid-to-upper tropospheric O3-producing precursors in this region. The out of phase relation suggests a UV modulation of tropospheric photochemistry caused by the QBO in stratospheric O3. During El Nino events there is anomalously low TCO in the eastern Pacific and high values in the western Pacific, indicating the effects of convectively-driven transport of low-value boundary layer O3 (reducing TCO) and O3 precursors including H2O and OH. A simplified technique is proposed to derive high-resolution maps of TCO in the tropics even in the absence of tropopause-level clouds. This promising approach requires only total ozone gridded measurements and utilizes the small variability observed in TCO near the dateline. This technique has an advantage compared to the CCD method because the latter requires high-resolution footprint measurements of both reflectivity and total ozone in the presence of tropopause-level cloud tops.

Description: Tropospheric column ozone (TCO) and stratospheric column ozone (SCO) gridded data in the tropics for 1979-present are now available from NASA Goddard Space Flight Center via either direct ftp, world-NN,ide-NN,eb, or electronic mail. This note provides a brief overview of the method used to derive the data set including validation and adjustments.

Description: In this study, we examine the sensitivity of long lived tracers to changes in the base transport components in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric mean age, and increased the rate of removal of material from the stratosphere. Increasing the stratospheric K(sub yy) increased the mean age due to the greater recycling of air parcels through the middle atmosphere, via the residual circulation, before returning to the troposphere. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and produces significantly younger ages throughout the stratosphere. Simulations with very small tropical stratospheric K(sub yy) decreased the globally averaged age of air by as much as 25% in the middle and upper stratosphere, and resulted in substantially weaker vertical age gradients above 20 km in the extratropics. We found only very small stratospheric tracer sensitivity to the magnitude of the horizontal mixing across the tropopause, and to the strength of the mesospheric gravity wave drag and diffusion used in the model. We also investigated the transport influence on chemically active tracers and found a strong age-tracer correlation, both in concentration and calculated lifetimes. The base model transport gives the most favorable overall comparison with a variety of inert tracer observations, and provides a significant improvement over our previous 1995 model transport. Moderate changes to the base transport were found to provide modest agreement with some of the measurements. Transport scenarios with residence times ranging from moderately shorter to slightly longer relative to the base case simulated N2O lifetimes that were within the observational estimates of Volk et al. [1997]. However, only scenarios with rather fast transport rates were comparable with the Volk et al. estimates of CFCl3 lifetimes. This is inconsistent with model-measurement comparisons of mean age in which the base model or slightly slower transport rates compared the most favorably with balloon SF6 data. For all comparisons shown, large transport changes away from the base case resulted in simulations that were outside the range of measurements, and in many cases, far outside this range.

Description: An intensive field campaign involving measurement of various aerosol physical, chemical, and radiative properties was conducted at Sde Boker in the Negev Desert of Israel, from 18 February to 15 March 1997. Nephelometer measurements gave average background scattering coefficient values of about 25 M/m at 550 nm wavelength, but strong dust events caused the value of this parameter to rise up to about 800 M/m Backscattering fractions did not depend on aerosol loading, and generally fell in the range of 0.1 to 0.25, comparable to values reported for marine and Arctic environments. Chemical analysis of the aerosol revealed that, in the coarse size range (2 - 10 micrometer equivalent aerodynamic diameter (EAD)), calcium (Ca) was by far the most abundant element followed by silicon (Si), both of which are indicators for mineral dust. In the fine size fraction (〈 2 micrometers EAD), sulfur (S) generally was the dominant element, except during high dust episodes when Ca and Si were again the most abundant. Furthermore, fine black carbon (BC) correlates with S, suggesting that they may have originated from the same sources or source regions. An indication of the short-term effect of aerosol loading on radiative forcing was provided by measurements of global and diffuse solar radiation, which showed that during high turbidity periods (strong dust events) almost all of the solar radiation reaching the area is scattered or absorbed.

Description: In 1891, the Virunga Mountains of Eastern Zaire were first acknowledged as volcanoes, and since then, the Virunga Mountain chain has demonstrated its potentially violent volcanic nature. The Virunga Mountains lie across the Eastern African Rift in an E-W direction located north of Lake Kivu. Mt. Nyamuragira and Mt. Nyiragongo present the most hazard of the eight mountains making up Virunga volcanic field, with the most recent activity during the 1970-90's. In 1977, after almost eighty years of moderate activity and periods of quiescence, Mt. Nyamuragira became highly active with lava flows that extruded from fissures on flanks circumscribing the volcano. The flows destroyed vast areas of vegetation and Zairian National Park areas, but no casualties were reported. Mt. Nyiragongo exhibited the same type volcanic activity, in association with regional tectonics that effected Mt. Nyamuragira, with variations of lava lake levels, lava fountains, and lava flows that resided in Lake Kivu. Mt. Nyiragongo, recently named a Decade volcano, presents both a direct and an indirect hazard to the inhabitants and properties located near the volcano. The Virunga volcanoes pose four major threats: volcanic eruptions, lava flows, toxic gas emission (CH4 and CO2), and earthquakes. Thus, the volcanoes of the Eastern African volcanic field emanate harm to the surrounding area by the forecast of volcanic eruptions. During the JSC Summer Fellowship program, we will acquire and collate remote sensing, photographic (Space Shuttle images), topographic and field data. In addition, maps of the extent and morphology(ies) of the features will be constructed using digital image information. The database generated will serve to create a Geographic Information System for easy access of information of the Eastem African volcanic field. The analysis of volcanism in Eastern Africa will permit a comparison for those areas from which we have field data. Results from this summer's work will permit further study and monitoring of the volcanic activity in the area. This is of concern due to the large numbers of refugees fleeing into Zaire where they are being positioned at the base of Mt. Nyiragongo. The refugees located at the base of the volcano are in direct hazard of suffocation by gas emission and destruction by lava flow. The results from this summer study will be used to secure future funding to enable continuation of this project.

Description: The orbital decay curve of a satellite having constant cross-sectional area and in circular orbit can furnish valuable information regarding the dynamical state of the atmosphere. It is shown that a rectilinear decay curve having constant downward slope (zero curvature) should indicate that the atmosphere was undergoing compression during that period. A decay curve having concavity upwards (positive curvature) will strongly indicate that the atmosphere was in a contracting phase. A decay curve with downward concavity (negative curvature) may indicate an expanding, a stationary or a contracting atmosphere. This theory, when applied to the orbital decay of the Orbital Debris Radar Calibration Spheres (ODERACS) satellites, indicates that during the period from Day 90 through Day 240 in the year 1994, the atmosphere was very definitely in a compression mode. During this period, ODERACS Sphere 1 faced nearly constant densities while Sphere 6 actually encountered progressively smaller air densities as they descended. The atmospheric scale height as calculated from the orbital data of Spheres 1 and 6 diminished steadily during the same period. It is shown that Spheres 1 and 6 descended faster and slower respectively, than the level of constant air density equal to 5 x 10 kg/m . During a brief period from Day 240 through Day 290, the atmosphere reversed to a strongly expanding mode. Thereafter, the atmosphere reverted back to a compression mode from Day 290 through Day 390, 1994.

Description: This presentation provides a brief overview of the scientific activities of the International GPS Service (IGS). This was an approved activity of the International Association of Geodesy (IAG) with official start of service on 1 Jan 1994. The mission of the IGS is "To provide a service to support geodetic and geophysical research activities, through GPS data and data products." The presentation explains the concept of the IGS working group, and pilot projects, and reviews the current working groups and pilot projects.

Description: Mantle-wide heterogeneity is largely controlled by deeply penetrating thermal convective currents. These thermal currents are likely to produce significant lateral variation in rheology, and this can profoundly influence overall material behaviour. How thermally related lateral viscosity variations impact models of glacio-isostatic and tidal deformation is largely unknown. An important step towards model improvement is to quantify, or bound, the actual viscosity variations that characterize the mantle. Simple scaling of viscosity to shear-wave velocity fluctuations yields map-views of long- wavelength viscosity variation. These give a general quantitative description and aid in estimating the depth dependence of rheological heterogeneity throughout the mantle. The upper mantle is probably characterized by two to four orders of magnitude variation (peak-to-peak). Discrepant time-scales for rebounding Holocene shorelines of Hudson Bay and southern Iceland are consistent with this characterization. Results are given in terms of a local average viscosity ratio, (Delta)eta(bar)(sub i), of volumetric concentration, phi(sub i). For the upper mantle deeper than 340 km the following reasonable limits are estimated for (delta)eta(bar) approx. equal 10(exp -2): 0.01 less than or equal to phi less than or equal to 0.15. A spectrum of ratios (Delta)eta(bar)(sub i) less than 0.1 at concentration level eta(sub i) approx. equal 10(exp -6) - 10(exp -1) in the lower mantle implies a spectrum of shorter time-scale deformational response modes for second-degree spherical harmonic deformations of the Earth. Although highly uncertain, this spectrum of spatial variation allows a purely Maxwellian viscoelastic rheology simultaneously to explain all solid tidal dispersion phenomena and long-term rebound-related mantle viscosity. Composite theory of multiphase viscoelastic media is used to demonstrate this effect.

Description: Within the last 7-8 years, there has been a substantial growth in out knowledge of the solar and interplanetary causes of geomagnetic storms at Earth. This review article will not attempt to cover all of the work done during this period. This can be found elsewhere. Our emphasis here will be on recent efforts that expose important, presently unanswered questions that must be addressed and solved before true predictability of storms can be possible. Hopefully, this article will encourage some readers to join this effort and perhaps make major contributions to the field.

Description: Stratospheric mixing of CH3D from 100 mb to 17 mb (?to 28 km) and HDO from 100 mb to 10.

Description: A review of the geomagnetic response to large-amplitude interplanetary Alfven wave trains is presented, highlighting its dependence on solar activity level and its role in the storm/substorm relationship problem. Also discussed are some recent observations obtained by the Ulysses spacecraft at high heliospheric latitudes dealing with the origin and dynamics of these wave trains.

Description: The secular drift of the precession of the ascending node of the LAGOES -1 satellite is apparently linked to the Earth s paleoclimate through the slow viscous response of the mantle to ice sheet/ocean mass transfer during the last great continental deglaciation . The secular node acceleration is particularly sensitive to the longest wavelengths of the paleo -surface loading that have been memorized by the mantle glacio -isostatic flow. Tide gauge records for the last 130 years show a post-glacial rebound-corrected sea-level rise of 2.4 n 0.9 mm yr-1.

Description: Efflorescent salt crusts and associated sediments in Death Valley, California, were studied with.

Description: Slow progress is being made in quantifying the mechanical and dynamical states of the lower mantle. The large plastic strength of the mantle above the D layer exhibits a non-tidal gravity field coefficient that is driven by isostatic flow of the deep mantle, as suggested by Yoder et al.

Description: A description of the DSN VLBI data set and of last year's analysis can be found in last year's report. Other than including another year's data, the main changes in this year's analysis from last year's are in the use of meteorological data for determining tropospheric parameters and in the weighting of the data to account for the uncertainty in the observables caused by tropospheric effects and source structure. A priori dry zenith tropospheric delays were determined from barometric pressure measurements at the DSN sites, corrected for height differences between the pressure sensor and the antennas. A priori wet zenith tropospheric delays were derived from tables of monthly average wet zenith delays for each station, which are based on historical radiosonde data. The Lanyi function was used for mapping zenith tropospheric delays to observed elevations. the temperature at the top of the boundary layer, a parameter in the Lanyi function, was taken to be the 24-hour average of the surface temperature at the station. Adjustments to the wet troposphere zenith delays were estimated every two to three hours.

Description: The Galileo spacecraft's Near-Infrared Mapping Spectrometer.

Description: High resolution trajectory calculations are used to show the patterns of the.

Description: In an effort to formulate the high frequency coupling between antennas.

Description: Ground-based solar absorption spectra were measured from Fairbanks, Alaska (65 degrees North, 148 degrees West) and from Ny Alesund, Spitsbergen, (79 degrees North, 12 degrees East) from March to September 1997 by Fourier Transform Infra-Red (FTIR) spectrometers.

Description: Remote sensing of chlorine monoxide (CIO) by the Microwave Limb Sounder (MLS) experiment aboard the Upper Atmosphere Research Satellite (UARS) has provided global measurements of variations in stratospheric free chlorine (for 1991 to 1997).

Description: There has been a question on whether it is possible to measure the surface wind speeds of tropical cyclones with a spaceborne Ku-band scatterometer.

Description: We examine energetic particle transport at or near Corotating Interaction Regions (CIRs).

Description: In global geodetic solutions vertical rates of site motion are usually estimated relative to the geocenter (center of figure) of the solid earth. The velocity of the geocenter is estimated assuming that the plates are rigid, that the velocities of the plates equal those in NUVEL-1A (DeMets et al. 1990, 1994) and that the uplift, subsidence, and intraplate deformation due to glacial isostatic adjustment is negligible.

Description: With aircraft-mounted in-situ and remote sensing instruments for dynamical, thermal, and chemical measurements, we studied two cases of tropopause folding.

Description: Current passive-microwave rain-retrieval methods are largely based on databases built off-line using cloud models.

Description: Geodetic measurements using the Global Positioning System and other techniques show north-south shortening near Los Angeles to be fastest across the northern part of the metropolitan area, where an ESE-striking, 5- to 40-km-wide belt lying to the south of San Gabriel Mountains and to the north of downtown and West Los Angeles is shortening at 5 mm/yr.

Description: Data from terrestrial GPS receivers are being used in growing numbers of applications requiring precise tropospheric sensing. One emerging application is the calibration of water vapor measurements from spaceborne microwave radiometers.

Description: Trace gases measured by the Atmospheric Trace Molecule Spectroscopy (ATMOS) instrument during the Mar/Apr 1992(AT-1), Apr 1993(AT-2), and Nov 1994(AT-3) space-shuttle missions have been mapped into equivalent latitude/potential temperature (EqL/0) coordinates.

Description: The late Holocene glacial moraine chronology in the southernmost South American Andes includes four (sup 14)C dated Neoglacial advances and retreats.

Description: Airborne measurements of NO, total reactive nitrogen(NO(sub y)), O(sub 3), and condensation nuclei (CN) were made, together with other tracers such as CO and H(sub 2)O, within air traffic corridors over the U.S. and North Atlantic regions in October and November 1997.

Description: In this paper, we study numerically the excitation of plasma waves by electron beams, in the auroral region above 2000 km of altitude. We have solved the fully kinetic dispersion relation, using numerical method and found the real frequency and the growth rate of the plasma wave modes. We have examined the instability properties of low-frequency waves such as the Electromagnetic Ion Cyclotron (EMIC) wave as well as Lower-Hybrid (LH) wave in the range of high-frequency. In all cases, the source of free energy are electron beams propagating parallel to the geomagnetic field. We present some features of the growth rate modes, when the cold plasma parameters are changed, such as background electrons and ions species (H(+) and O(+)) temperature, density or the electron beam density and/or drift velocity. These results can be used in a test-particle simulation code, to investigate the ion acceleration and their implication in the auroral acceleration processes, by wave-particle interaction.

Description: (See scanned abstract).