ALBERT

Description: The TARFOX (Tropospheric Aerosol Radiative Forcing Observational Experiment) intensive field campaign was designed to reduce uncertainties in estimates of the effects of anthropogenic aerosols on climate by measuring direct radiative effects and the optical, physical, and chemical properties of aerosols [1]. TARFOX was conducted off the East Coast of the United States between July 10-31, 1996. Ground, aircraft, and satellite-based sensors measured the sensitivity of radiative fields at various atmospheric levels to aerosol optical properties (i.e., optical thickness, phase function, single-scattering albedo) and to the vertical profile of aerosols. The LASE (Lidar Atmospheric Sensing Experiment) instrument, which was flown on the NASA ER-2 aircraft, measured vertical profiles of total scattering ratio and water vapor during a series of 9 flights. These profiles were used in real-time to help direct the other aircraft to the appropriate altitudes for intensive sampling of aerosol layers. We have subsequently used the LASE aerosol data to derive aerosol backscattering and extinction profiles. Using these aerosol extinction profiles, we derived estimates of aerosol optical thickness (AOT) and compared these with measurements of AOT from both ground and airborne sun photometers and derived from the ATSR-2 (Along Track and Scanning Radiometer 2) sensor on ERS-2 (European Remote Sensing Satellite-2). We also used the water vapor mixing ratio profiles measured simultaneously by LASE to derive precipitable water vapor and compare these to ground based measurements.

Description: In recent years, the role of PSCs in the ozone depletion process has become better understood. PSCs provide the surfaces upon which heterogeneous reactions take place that affect the gas phase partitioning between active and reservoir chlorine and nitrogen species. Present methods of PSC detection include in situ measurements by lidar and various satellite-borne instruments such as the Stratospheric Aerosol Measurement II (SAM II) on the Nimbus 7 spacecraft, which produced PSC measurements from 1978 to 1994 and several instruments onboard the Upper Atmospheric Research Satellite (UARS) such as the Cryogenic Limb Array Etalon Spectrometer (CLAES) which provided measurements for 1991-1993. All of the PSC-detection methods devised so far have been hampered by incomplete sampling of the places and times in which PSCs are likely to form. There is a need to understand the climatology of PSCS, in particular the timing of their onset and duration, their vertical distribution, geographic extent, annual variability and responses to volcanic aerosol forcing. Poole and Pitts [1994] assembled a PSC climatology based on SAM II data, but this climatology is incomplete, as it is limited to the edge of the polar night due to the limitations of the solar occultation scan geometry. The Advanced Very High Resolution Radiometer (AVHRR) five- channel sensors onboard the NOAA polar-orbiting satellites have been collecting data over the polar regions continuously since 1979. These operational satellites provide unmatched coverage in space and time of both polar regions, but were not designed for the detection of optically-thin PSCS. However, the AVHRR data archive would be an invaluable source for the construction of a long-term climatology of PSCs if techniques can be developed and tested to detect PSCs in AVHRR data. In the last few years, the members of our group at San Francisco State University and NASA Ames Research Center have been engaged in the development of various PSC detection methods using AVHRR data. There is strong evidence that a subset of PSCS, those that are optically thick, can be readily identified in the AVHRR data set. Our group has also made significant progress in the identification of optically thinner PSCs using a variety of techniques.

Description: Lidar backscattering profiles available from the LITE data set have been used to estimate the optical depths of the Saharan dust layer over West Africa and E. Atlantic regions, in the context of validating the 3-D conceptual model of the Saharan dust plume proposed by Karyampudi and Carlson. The aerosol extinction profiles and optical depths were retrieved from LITE using the Fernald et al. (1972) method. An extinction-to-backscattering ratio, S(sub a), of 25 was selected for optical depth calculations. The spatial analysis of total column and Saharan dust layer optical depths show higher optical depths over W. Africa that decrease westward over E. Atlantic. The higher optical depths over W. Africa, in general, are associated with heavy dust being raised from the surface in dust source regions. Rapid depletion of these heavy dust particles, perhaps due to sedimentation, appear to decrease the dust loading within the dust layer as the plume leaves the west African continent. Higher optical depths are generally confined to the southern edge of the dust layer, where the middle level jet appears to transport the heavy dust concentrations that tend to mix downward from vertical mixing associated with the strong vertical shears underneath the middle jet. Thus, LITE measurements although, in general, validate the Saharan dust plume conceptual model, show maximum values of optical depths near the southern edge of the dust plume over the E. Atlantic region instead of near the center of the dust plume as described in the conceptual model.

Description: Water vapor plays an important role in the energetics of the boundary layer processes which in turn play a key role in regulating regional and global climate. It plays a primary role in Earth's hydrological cycle, in radiation balance as a direct absorber of infrared radiation, and in atmospheric circulation as a latent heat energy source, as well as in determining cloud development and atmospheric stability. Water vapor concentration, expressed as a mass mixing ratio (g kg(exp -l)), is conserved in all meteorological processes except condensation and evaporation. This property makes it an ideal choice for studying many of the atmosphere's dynamic features. Raman scattering measurements from lidar also allow retrieval of water vapor mixing ratio profiles at high temporal and vertical resolution. Raman lidars sense water vapor to altitudes not achievable with towers and surface systems, sample the atmosphere at much higher temporal resolution than radiosondes or satellites, and do not require strong vertical gradients or turbulent fluctuations in temperature that is required by acoustic sounders and radars. Analysis of highly-resolved water vapor profiles are used here to characterize two important mesoscale flows: thunderstorm outflows and a cold front passage. The data were obtained at the Atmospheric Radiation Measurement Site (CART) by the groundbased Department of Energy/Sandia National Laboratories lidar (CART Raman lidar or CARL) and Goddard Space Flight Center Scanning Raman Lidar (SRL). A detailed discussion of the SRL and CARL performance during the IOPs is given by others in this meeting.

Description: The NASA Langley Research Center's airborne UV Differential Absorption Lidar (DIAL) system participated in the Subsonic Assessment, Ozone and Nitrogen Oxide Experiment (SONEX) mission from October 13 to November 12, 1997. The purpose of the mission was to study the upper troposphere/lower stratosphere in and near the North Atlantic flight corridor to better understand this region of the atmosphere and how civilian air travel in the corridor might be affecting the atmospheric chemistry. Bases of operations included NASA Ames, California (37.4 deg N, 122.1 deg W); Bangor, Maine (44.8 deg N, 68.8 deg W); Shannon, Ireland (52.7 deg N, 8.9 deg W); and Lajes, Terceira Island, Azores (38.8 deg N, 27.1 deg W). Since the UV DIAL system observes in the nadir as well as the zenith, aerosol and ozone data were obtained from near the Earth's surface to the lower stratosphere. A number of interesting features were noted relating to both chemistry and dynamics of the troposphere, which are reported here.

Description: Worldwide, about ten Differential Absorption Lidars are used for long-term monitoring of stratospheric ozone. These systems are an important component of the Network for the Detection of Stratospheric Change. Although DIALs are self-calibrating in principle, regular intercomparisons with other ozone-lidars, microwave radiometers or ozone-sondes are highly desirable to ensure high data quality at a well known level. The Network for the Detection of Stratospheric Change (NDSC) validation policy suggests that such intercomparisons be "blind", meaning all participants submit their data to an impartial referee, without seeing results from the other participants. Here we report on the "blind" intercomparison taking place from January 20th to February 10th 1998 at Ny-Alesund, Spitsbergen (78.92 deg N, 11.95 deg E). Participating groups were from the Alfred Wegener Institute, Potsdam, operating the NDSC DIAL system at Ny-Alesund, from the University of Bremen operating the NDSC microwave radiometer for ozone profiling at Ny-Alesund, and the NASA Goddard Space Flight Center group with the "NDSC travelling standard" STROZ-LITE. The first author acted as the impartial referee. Also used for the intercomparison were data from ECC-6A/Vaisala RS80 ozone sondes routinely launched at Ny-Alesund by the AWI group. A 1% KI solution (3 ml) and the 1986 ECC pump correction (1.092 at 5 hPa) are used. The ECC-data were available to all participants during the campaign and thus were not "blind". Table 1 summarizes the expected performance of the instruments participating in the ozone intercomparison reported in this paper.

Description: Many modem and ancient carbonate deposits around the world have been recognized as microbial buildups or microbialites. Ancient microbialite structures have been divided into two basic categories based on their internal fabric or texture. They include stromatolites which have a predominantly laminated internal fabric and thrombolites which have an open-porous clotted fabric, that lacks laminae. The origin of these two basic microbial fabrics is still being debated in the literature. Understanding the origin and the various microorganisms involved in forming these modem fabrics is the key to the interpretation of similar fabrics in ancient and possibly Martian rocks. Therefore, detailed studies are needed on the microbiological makeup and origin of the fabrics in modem microbialites. Such studies may serve as analogs for ancient and Martian microbialites in the future. The purpose of this study is to examine the textures and carbon isotopic signatures of the following modem microbialites from the Bahamas: 1) a modem subtidal microbialite from Iguana Cay, Bahamas and 2) a modem microbial mat (stromatolite) from a hypersaline pond on Lee Stocking Island, Bahamas.

Description: Research has established the importance of global tropospheric wind measurements for large scale improvements in numerical weather prediction. In addition, global wind measurements provide data that are fundamental to the understanding and prediction of global climate change. These tasks are closely linked with the goals of the NASA Earth Science Enterprise and Global Climate Change programs. NASA Goddard has been actively involved in the development of direct detection Doppler lidar methods and technologies to meet the wind observing needs of the atmospheric science community. A variety of direct detection Doppler wind lidar measurements have recently been reported indicating the growing interest in this area. Our program at Goddard has concentrated on the development of the edge technique for lidar wind measurements. Implementations of the edge technique using either the aerosol or molecular backscatter for the Doppler wind measurement have been described. The basic principles have been verified in lab and atmospheric lidar wind experiments. The lidar measurements were obtained with an aerosol edge technique lidar operating at 1064 nm. These measurements demonstrated high spatial resolution (22 m) and high velocity sensitivity (rms variances of 0.1 m/s) in the planetary boundary layer (PBL). The aerosol backscatter is typically high in the PBL and the effects of the molecular backscatter can often be neglected. However, as was discussed in the original edge technique paper, the molecular contribution to the signal is significant above the boundary layer and a correction for the effects of molecular backscatter is required to make wind measurements. In addition, the molecular signal is a dominant source of noise in regions where the molecular to aerosol ratio is large since the energy monitor channel used in the single edge technique measures the sum of the aerosol and molecular signals. To extend the operation of the edge technique into the free troposphere we have developed a variation of the edge technique called the double edge technique. In this paper a ground based aerosol double edge lidar is described and the first measurements of wind profiles in the free troposphere obtained with this lidar will be presented.

Description: Earth's thin layer of soil is a fragile resource, made up of minerals, organic materials, air, water, and billions of living organisms. Soils plays a variety of critical roles that sustain life on Earth. If we think about soil, we tend to see it first as the source of most of the food we eat and the fibers we use, such as wood and cotton. Few students realize that soils also provide the key ingredients to many of the medicines (including antibiotics), cosmetics, and dyes that we use. Fewer still understand the importance of soils in integrating, controlling, and regulating the movement of air, water, materials, and energy between the hydrosphere, lithosphere, atmosphere, and biosphere. Because soil sustains life, it offers both a context and a natural laboratory for investigating these interactions. The enclosed poster, which integrates soil profiles with typical landscapes in which soils form, can also help students explore the interrelationships of Earth systems and gain an understanding of our soil resources. The poster, produced jointly by the American Geological Institute and the Soil Science Society of America, aims to increase awareness of the importance of soil, as does the GLOBE (Global Learning and Observations To Benefit the Environment) Program. Vice President Al Gore instituted the GLOBE Program on Earth Day of 1993 to increase environmental awareness of individuals throughout the world, contribute to a better scientific understanding of the Earth, and help all students reach higher levels of achievement in science and mathematics. GLOBE functions as a partnership between scientists, students, and teachers in which scientists design protocols for specific measurements they need for their research that can be performed by K-12 students. Teachers are trained in the GLOBE protocols and teach them to their students. Students make the measurements, enter data via the Internet to a central data archive, and the data becomes available to scientists and the general community. Students benefit by having a "hands-on"experience in science, math, and technology, using their local environment as a learning laboratory, as well as contact with scientists and other students around the world. Soil investigations have become an essential component of GLOBE. The protocols that have been developed so far within the GLOBE program include GPS Location, Atmosphere/Climate, Soil Characterization, Soil Moisture and Temperature, Land Cover/Biometry, Hydrology, and Satellite Image Classification. For the GLOBE Soil Characterization Protocol, students explore the physical. chemical, and morphological properties of the soil at their study site. They are asked to dig a pit or use an auger to about 1 meter at at least 2 sites.

Description: This study investigates the use of H2, p-synthesis, and mixed H2/mu methods to construct full-order controllers and optimized controllers of fixed dimensions. The benchmark problem definition is first extended to include uncertainty within the controller bandwidth in the form of parametric uncertainty representative of uncertainty in the natural frequencies of the design model. The sensitivity of H2 design to unmodelled dynamics and parametric uncertainty is evaluated for a range of controller levels of authority. Next, mu-synthesis methods are applied to design full-order compensators that are robust to both unmodelled dynamics and to parametric uncertainty. Finally, a set of mixed H2/mu compensators are designed which are optimized for a fixed compensator dimension. These mixed norm designs recover the H, design performance levels while providing the same levels of robust stability as the u designs. It is shown that designing with the mixed norm approach permits higher levels of controller authority for which the H, designs are destabilizing. The benchmark problem is that of an active tendon system. The controller designs are all based on the use of acceleration feedback.

Description: Vertical mixing ratio profiles of four relatively long-lives gases, HCN, C2H2, CO, and C2H6, have been retrieved from 0.01/cm resolution infrared solar occultation spectra recorded between latitudes of 5.3degN and 31.4degN. The observations were obtained by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer during the Atmospheric Laboratory for Applications and Science (ATLAS) 3 shuttle flight, 3-12 November 1994. Elevated mixing ratios below the tropopause were measured for these gases during several of the occultations. The positive correlations obtained between the simultaneously measured mixing ratios suggest that the enhancements are likely the result of surface emissions, most likely biomass burning and/or urban industrial activities, followed by common injection via deep convective transport of the gases to the upper troposphere. The elevated levels of HCN may account for at least part of the "missing NO," in the upper troposphere. Comparisons of the observations with values measured during a recent aircraft campaign are presented.

Description: An autonomous GPS array is being implemented in the Transantarctic Mountains, sponsored by NSF and NASA, for the purpose of measuring uplift resulting from post-glacial rebound (PGR). The rebound of the solid earth due to unloading of ice since the Last Glacial Maximum is expected to dominate the measured uplift for most of West Antarctica, dwarfing the signals due to present-day ice sheet mass balance changes and tectonic motion, as long as mantle viscosity is greater than about 10(exp 20) Pa-s. Predicted uplift patterns have been calculated for a range of model scenarios, which illustrate how the uplift pattern might distinguish between different-sized ice sheets and deglaciation histories as represented by the competing models. The scenarios considered by James and Ivins (1998) include ICE-3G, CLIMAP and a variation of the CLIMAP model by Denton et al. For these models, peak uplift rates occur in the Transantarctic Mountains, and differences between models is often large there. Thus, the Transantarctic Mountains are an ideal place to obtain uplift measurements to constrain deglaciation models.

Description: In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.

Description: As a result of photochemistry, some relationship between the stratospheric age-of-air and the amount of tracer contained within an air sample is expected. The existence of such a relationship allows inferences about transport history to be made from observations of chemical tracers. This paper lays down the conceptual foundations for the relationship between age and tracer amount, developed within a Lagrangian framework. In general, the photochemical loss depends not only on the age of the parcel but also on its path. We show that under the "average path approximation" that the path variations are less important than parcel age. The average path approximation then allows us to develop a formal relationship between the age spectrum and the tracer spectrum. Using the relation between the tracer and age spectra, tracer-tracer correlations can be interpreted as resulting from mixing which connects parts of the single path photochemistry curve, which is formed purely from the action of photochemistry on an irreducible parcel. This geometric interpretation of mixing gives rise to constraints on trace gas correlations, and explains why some observations are do not fall on rapid mixing curves. This effect is seen in the ATMOS observations.

Description: The Global Modeling Initiative (GMI) science team is developing a three dimensional chemistry and transport model (CTM) to be used in assessment of the atmospheric effects of aviation. Requirements are that this model be documented, be validated against observations, use a realistic atmospheric circulation, and contain numerical transport and photochemical modules representing atmospheric processes. The model must also retain computational efficiency to be tractable to use for multiple scenarios and sensitivity studies. To meet these requirements, a facility model concept was developed in which the different components of the CTM are evaluated separately. The first use of the GMI model will be to evaluate the impact of the exhaust of supersonic aircraft on the stratosphere. The assessment calculations will depend strongly on the wind and temperature fields used by the CTM. Three meteorological data sets for the stratosphere are available to GMI: the National Center for Atmospheric Research Community Climate Model (CCM2), the Goddard Earth Observing System Data Assimilation System (GEOS DAS), and the Goddard Institute for Space Studies general circulation model (GISS). Objective criteria were established by the GMI team to identify the data set which provides the best representation of the stratosphere. Simulations of gases with simple chemical control were chosen to test various aspects of model transport. The three meteorological data sets were evaluated and graded based on their ability to simulate these aspects of stratospheric measurements. This paper describes the criteria used in grading the meteorological fields. The meteorological data set which has the highest score and therefore was selected for GMI is CCM2. This type of objective model evaluation establishes a physical basis for interpretation of differences between models and observations. Further, the method provides a quantitative basis for defining model errors, for discriminating between different models, and for ready re-evaluation of improved models. These in turn will lead to a higher level of confidence in assessment calculations.

Description: Supported by numerical experiment results, the abrupt change of the location of the intertropical convergence zone (ITCZ), from the equatorial trough flow regime to the monsoon trough flow regime is interpreted as a subcritical instability. The existence of these multiple quasi-equilibria is due to the balance of two "forces" on the ITCZ: one toward the equator, due to the earth's rotation, has a nonlinear latitudinal dependence; and the other toward the latitude of the sea surface (or ground) temperature peak has a relatively linear latitudinal dependence. This work pivots on the finding that the ITCZ and Hadley circulation can still exist without the pole-to-equator gradient of radiative-convective equilibrium temperature.

Description: Solar radiation is the primary source of energy driving atmospheric and oceanic circulations. Concerned with the huge computing time required for computing radiative transfer in weather and climate models, solar heating in minor absorption bands has often been neglected. The individual contributions of these minor bands to the atmospheric heating is small, but collectively they are not negligible. The solar heating in minor bands includes the absorption due to water vapor in the photosynthetically active radiation (PAR) spectral region from 14284/cm to 25000/cm, the ozone absorption and Rayleigh scattering in the near infrared, as well as the O2 and CO2 absorption in a number of weak bands. Detailed high spectral- and angular-resolution calculations show that the total effect of these minor absorption is to enhance the atmospheric solar heating by approximately 10%. Depending upon the strength of the absorption and the overlapping among gaseous absorption, different approaches are applied to parameterize these minor absorption. The parameterizations are accurate and require little extra time for computing radiative fluxes. They have been efficiently implemented in the various atmospheric models at NASA/Goddard Space Flight Center, including cloud ensemble, mesoscale, and climate models.

Description: Volcanic eruptions loft gases and ash particles into the atmosphere and produce effects that are both short term (aircraft hazards, interference with satellite measurements) and long term (atmospheric chemistry, climate). Large (greater than 0.5mm) ash particles fall out in minutes [Rose et al, 1995], but fine ash particles can remain in the atmosphere for many days. This fine volcanic ash is a hazard to modem jet aircraft because the operating temperatures of jet engines are above the solidus temperature of volcanic ash, and because ash causes abrasion of windows and airframe, and disruption of avionics. At large distances(10(exp 2)-10(exp 4) km or more) from their source, drifting ash clouds are increasingly difficult to distinguish from meteorological clouds, both visually and on radar [Rose et al., 1995]. Satellites above the atmosphere are unique platforms for viewing volcanic clouds on a global basis and measuring their constituents and total mass. Until recently, only polar AVHRR and geostationary GOES instruments could be used to determine characteristics of drifting volcanic ash clouds using the 10-12 micron window [Prata 1989; Wen and Rose 1994; Rose and Schneider 1996]. The NASA Total Ozone Mapping Spectrometer (TOMS) instruments aboard the Nimbus-7, Meteor3, ADEOS, and Earth Probe satellites have produced a unique data set of global SO2 volcanic emissions since 1978 (Krueger et al., 1995). Besides SO2, a new technique has been developed which uses the measured spectral contrast of the backscattered radiances in the 330-380nm spectral region (where gaseous absorption is negligible) in conjunction with radiative transfer models to retrieve properties of volcanic ash (Krotkov et al., 1997) and other types of absorbing aerosols (Torres et al., 1998).

Description: This study investigates the distribution of clear-sky ultraviolet-B (UV-B, wavelengths 290-320 nm) trends in northern midlatitudes using 1979-1991 Nimbus 7 total ozone mapping spectrometer (TOMS) version 7 low-reflectivity (R〈0.2) total ozone footprint measurements. The incorporation of essentially cloud-free ozone data from TOMS provides a direct method for separating transient cloud effects from anthropogenic and other dynamical factors present in UV-B. This study has also included both National Oceanic and Atmospheric Administration (NOAA) microwave sounding unit channel 4 (MSU4) and National Centers for Environmental Prediction (NCEP) 500 hPa temperature (T500) fields in our trend models to improve UV-Index (UVI) trend statistics and to investigate the effects of interannual changes in UVI caused by synoptic-scale (horizontal wavelengths 4000-8000 km) and planetary-scale (horizontal wavelengths greater than 8000 km) dynamical events. Clear-sky UVI trends in the northern midlatitudes show large increases (exceeding 10 % per decade) and distinct regional variability especially during winter-spring months which can be attributed to topography and dynamical forcing effects. In the UV-important summer-autumn months, these trends are more uniformly distributed and still statistically significant, although smaller at around +2 to +3 % per decade. Specifically, during April largest increases in midlatitude UVI are seen to extend from near the dateline eastward across North America. In June months largest UVI increases occur over the east Asian continent with values around +5 to +6 % per decade. These increases in UVI over both the Pacific and Asian continent regions persist through summer into Autumn. In the the European sector, statistically significant increases in clear-sky UVI are found over central Europe with values around +2 to +3 % per decade and +8 to +9 % per decade during summer and winter-spring months, respectively. Over the nearby Mediterranean region these seasonal trends are around +2 to +3 and +5 to +6 % per decade.

Description: In August 1995 multiple instruments that measure the stratospheric ozone vertical distribution were intercompared at the Mauna Loa Observatory, Hawaii, under the auspices of the Network for the Detection of Stratospheric Change. The instruments included two UV lidar systems, one from JPL and the other from Goddard Space Flight Center, ECC balloon-sondes, a ground-based microwave instrument, Umkehr measurements, and a new ground-based FTIR instrument. The MLS instrument on the UARS satellite provided correlative profiles of ozone, and there was one close overpass of the SAGE II instrument. The results show that much better consistency among instruments is being achieved than even a few years ago, usually to within the instrument uncertainties. The different measurement techniques in this comparison agree to within +/-10% at almost all altitudes, and in the 20 km to 45 km region most agreed within +/-5%. The results show that the current generation of lidars are capable of accurate measurement of the ozone profile to a maximum altitude of 50 km. SAGE agreed well with both lidar and balloon-sonde down to at least 17 km. The ground-based microwave measurement agreed with other measurements from 22 km to above 50 km. One minor source of disagreement continues to be the pressure-altitude conversion needed to compare a measurement of ozone density versus altitude with a measurement of ozone mixing ratio versus pressure.

Description: The decadal, 12-14 year, cycle observed in the North Atlantic SST and tide gauge data was examined using the NCEP/NCAR reanalyses, COADS data and an ocean model simulation. Besides this decadal mode, a shorter, subdecadal period of about 8 years exists in tide gauge data north of 40N, in the subpolar SST and in the winter North Atlantic Oscillation (NAO) index and in subpolar winter heat flux values. The decadal cycle is a well separated mode in a singular spectrum analysis (SSA) for a time series of SST EOF mode 1 with a center over the Gulf Stream extension. Tide gauge and SST data are consistent in that both show a significant subdecadal periodicity exclusively in the subpolar gyre, but in subtropics the 12-14 year period is the prominent, but nonstationary, decadal signal. The main finding of this study is that this 12-14 year cycle can be constructed based on the leading mode of the surface heat flux. This connection to the surface heat flux implicates the participation of the thermohaline circulation in the decadal cycle. During the cycle starting from the positive index phase of NAO, SST and oceanic heat content anomalies are created in subtropics due to local heat flux and intensification of the thermohaline circulation. The anomalies advect to the subpolar gyre where they are amplified by local heat flux and are part of the negative feedback of thermohaline circulation on itself. Consequently the oceanic thermohaline circulation slows down and the opposite cycle starts. The oscillatory nature would not be possible without the active atmospheric participation in the cycle, because it provides the unstable interaction through heat flux, without it, the oceanic mode would be damped. This analysis suggests that the two principal modes of heat flux variability, corresponding to patterns similar to North Atlantic Oscillation (NAO) and Western Atlantic (WA), are part of the same decadal cycle and an indirect measure of the north-south movement of the storm tracks.

Description: HCl measurements from HALOE in the northern hemisphere during mid-May 1997 revealed vortex fragments in which the chlorine reservoir partitioning was strongly pushed toward HCl (approx. 90% HCl, approx. 10% ClONO2), similar to partitioning previously observed in the Antarctic vortex region. In contrast, observations of ClONO2 and HCl in the northern polar spring, 1992, and in other years, show these species established the balance typical for gas phase photochemical reactions in this region (approx. 60% HCl, approx. 40% ClONO2). Annually, chlorine reservoirs in the winter lower stratosphere polar vortex are converted to chlorine radicals via heterogeneous reactions on particle surfaces at very cold temperatures (less than about 200 K). As temperatures warm in spring, the heterogeneous processes become insignificant compared with gas phase reactions, and the chlorine reservoirs are reformed. Measurements through the northern winter/spring in 1992 show rapid formation of ClONO2, followed by steady loss of ClONO2 and increasing HCl. Although ClONO2 measurements are not available for 1997, the HCl increase in 1997 is observed to be much more rapid and the eventual HCl mixing ratio is about 50% greater than that of 1992. The observations are examined through comparison with the Goddard three-dimensional chemistry and transport model. This model utilizes winds and temperatures from the Goddard Earth Observing System Data Assimilation System and a complete integration scheme for stratospheric photochemistry. Analysis of the evolution of HCl and ClONO2 shows that the observed difference in the overall rate of HCl formation is explained by the sensitivity of the gas-phase chemistry to the ozone mixing ratio and the temperature. The results show that the model accurately simulates HCl and ClONO2 evolution during these two winters. Model validity is further supported by comparisons with O3 and reactive nitrogen species NO and NO2. This analysis provides a sensitive test of the lower stratospheric chlorine photochemistry, particularly because the analysis considers constituent evolution at a time when the HCl and ClONO2 are far from a photochemical stationary state.

Description: Vertical profiles of infrared cirrus extinction have been derived from tropical and subtropical upper tropospheric solar occultation spectra. The measurements were recorded by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer during the Atmospheric Laboratory for Applications and Sciences (ATLAS) 3 shuttle flight in November 1994. The presence of large numbers of small ice crystals is inferred from the appearance of broad extinction features in the 8-12 micron region. These features were observed near the tropopause and at lower altitudes. Vertical profiles of the ice extinction (/km) in microwindows at 831, 957, and 1204/cm have been retrieved from the spectra and analyzed with a model for randomly oriented spheroidal ice crystals. An area-equivalent spherical radius of 6 microns is estimated from the smallest ice crystals observed in the 8-12 gm region. Direct penetration of clouds into the lower stratosphere is inferred from observations of cloud extinction extending from the upper troposphere to 50 mbar (20 km altitude). Cloud extinction between 3 and 5 microns shows very little wavelength dependence, at least for the cases observed by the ATMOS instrument in the tropics and subtropics during ATLAS 3.

Description: A statistical description of Earth's broad scale, core-source magnetic field has been developed and tested. The description features an expected, or mean, spatial magnetic power spectrum that is neither "flat" nor "while" at any depth, but is akin to spectra advanced by Stevenson and McLeod. This multipole spectrum describes the magnetic energy range; it is not steep enough for Gubbins' magnetic dissipation range. Natural variations of core multipole powers about their mean values are to be expected over geologic time and are described via trial probability distribution functions that neither require nor prohibit magnetic isotropy. The description is thus applicable to core-source dipole and low degree non-dipole fields despite axial dipole anisotropy. The description is combined with main field models of modem satellite and surface geomagnetic measurements to make testable predictions of: (1) the radius of Earth's core, (2) mean paleomagnetic field intensity, and (3) the mean rates and durations of both dipole power excursions and durable axial dipole reversals. The predicted core radius is 0.7% above the 3480 km seismologic value. The predicted root mean square paleointensity (35.6 mu T) and mean Virtual Axial Dipole Moment (about 6.2 lx 1022 Am(exp 2)) are within the range of various mean paleointensity estimates. The predicted mean rate of dipole power excursions, as defined by an absolute dipole moment 〈20% of the 1980 value, is 9.04/Myr and 14% less than obtained by analysis of a 4 Myr paleointensity record. The predicted mean rate of durable axial dipole reversals (2.26/Myr) is 2.3% more than established by the polarity time-scale for the past 84 Myr. The predicted mean duration of axial dipole reversals (5533 yr) is indistinguishable from an observational value. The accuracy of these predictions demonstrates the power and utility of the description, which is thought to merit further development and testing. It is suggested that strong stable stratification of Earth's uppermost outer core leads to a geologically long interval of no dipole reversals and a very nearly axisymmetric field outside the core. Statistical descriptions of other planetary magnetic fields are outlined.

Description: This study shows that: 1) the NSCAT backscatter fields provide an estimate of the PIZ coverage of the Arctic Ocean; and, 2) the decrease in PIZ area over the winter gives an indication of the PIZ area exported through Fram Strait.

Description: The tidal signature in the middle atmospheric thermal structure is investigated using more than 140 hours of nighttime lidar measurements at Table Mountain (34.4 degrees north) during January 1997 and February 1998.

Description: In addition, the effective Henry's law constant and the associative enthalpy change of solvation of PNA in water are determined to be 39.95 mol kg^-1 atm^-1 and -69.84 kJ mol^-1 at 298.15K, respectively.

Description: The terrestrial location of the Earth's rotation pole has been under continuous observation since 1899 when the International Latitude Service (ILS) began conducting optical astrometric measurements of star positions to determine variations in station latitude and hence variations in the location of the rotation pole.

Description: The Ulysses mission explored for the first time our heliosphere at all latitudes up to +/-80 and therefore has been an ideal mission to study potential gradients in heliolatitude of discontinuity occurrence rates and types.

Description: Using more than 145 hours of nighttime lidar measurements obtained during October 1-16, 1996 and October 2-11, 1997, the tidal signature in the middle atmopheric thermal structure (15-95 km) at Mauna Loa, Hawaii, is investigated.

Description: One of the oldest mysteries in geomagnetism is the linkage between solar and geomagnetic activity. The 11-year cycles of both the numbers of sunspots and Earth geomagnetic storms were first noted by Sabine (1852).

Description: The circulation and upwelling processes (coastal and gyre-induced) that control the phytoplankton distribution in the Alboran sea are examined by analyzing monthly climatological patterns of Coastal Zone Color Scanner (CZCS) pigment concentrations, sea surface temperatures, winds, and seasonal geostrophic fields.

Description: Measurements of HCl in the lower stratosphere (15-21 km) from aircraft, balloon, Space Shuttle, and satellite reveal a growth in its mean abundance relative to that of total inorganic chlorine (Cly) from HCl/Cly = 57(+/-5)% in early 1993 to 75(+/-7)% by the end of 1997.

Description: The Ulysses mission explored for the first time our heliosphere at all latitudes up to +/-80 and therefore has been an ideal mission to study potential gradients in heliolatitude(and radial distance) of discontinuity occurrence rates and types.

Description: The radio signal from the Galileo prove to the orbiter experienced attenuation due to ammonia in Jupiter's atmosphere during the probe descent.

Description: Wind taken from the National Aeronautics and Space Administration (NASA) scatterometer (NSCAT) is compared with the operational analysis from European Center for Medium-Range Forecast (ECMWF) for the entire duration (about 9 months) of the NSCAT mission.

Description: We report the first finding of a pair of forward and reverse slow-mode shocks in the distant heliosphere using plasma and magnetic field data from the Ulysses spacecraft located at 5.3 AU and 9 degrees South heliolatitude.

Description: We estimate the winter sea ice export through the Fram Strait using ice motion from satellite passive microwave data. Sea ice motion (October through May)is obtained by tracking the displacement of common features in sequential 85 GHz and 37 GHz brightness temperature fields.

Description: The goal of the Global Rain Forest Mapping project (GRFM) is to acquire continuous, high resolution L-band Synthetic Aperture Radar (SAR) imagery of the major tropical rain forests of the world using the Japanese Earth Resources Satellite (JERS-1) SAR.

Description: The purpose of this paper is to study an interplanetary Bs feature ahead of a magnetic cloud and its related auroral and magnetospheric responses.

Description: A combination of Earth orientation measurements has been generated from space-geodetic observations spanning 1976-1997.

Description: A Kalman filter has been used to combine independent measurements of the Earth's orientation taken by optical astrometry and the space-geodetic techniques of LLR, SLR, VLBI, and GPS.

Description: Using a helium plus reversible fluid balloon system as the observing platform, multiple profiles of shortwave irradiance between 4 and 10 km were recently obtained over the Los Angeles basin.

Description: Results of an investigation of mesospheric temperature inversion layers using long-term lidar measurements at mid- and low-latitudes are reported. In this paper, new results from different lidar observations of the invasion layers will be presented.

Description: The ionospheric storm evolution process was monitored during the January 10, 1997 magnetic cloud event, through measurements of the inonospheric total electron content (TEC) from 150 GPS stations. The first significant response of the inonospheric TEC to the geomagnetic storm was at 0300 UT as an auroral/subauroral enhancement around the Alaskan evening sector.

Description: This paper presents a study of the equatorial diurnal tide at 27-55 km based on temperature data from the Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS).

Description: Our understanding of snow distributions in the polar regions is severly restricted due to the heterogeneity, both in space and time, of this solid precipitate. Processes such as vapor and mass fluxes across the interface are, to a large extent, controlled by the presence and geophysical state of the snow cover on sea ice.

Description: Radio occultation measurements, which probe electron density over a wide dynamic range with high sensitivity and high spatial and temporal resolution reveal a solar corona permeated by a hierarchy of filamentary structures.

Description: Upper tropospheric humidity (UTH) is a fundamental importance in understanding earth's atmosphere and climate. Water vapor is the most important greenhouse gas and it is in the upper troposphere that water vapor most strongly influences radiative forcing.

Description: The effect of Faraday rotation on linearly polarized backscatter measurements from space is addressed.

Description: Dual wavelength aerosol lidar observations of tropical high-altitude cirrus clouds were performed during the ALBATROSS 1996 campaign aboard the research vessel POLARSTERN on the Atlantic ocean in October-November 1996.

Description: This report describes the development of a three-dimensional scenario of aircraft fuel burn and emissions (fuel burned, NOx, CO, and hydrocarbons)for projected year 2015 scheduled air traffic. These emission inventories are available for use by atmospheric scientists conducting the Atmospheric Effects of Aviation Project (AEAP) modeling studies. Fuel burned and emissions of nitrogen oxides (NOx as NO2), carbon monoxides, and hydrocarbons have been calculated on a 1 degree latitude x 1 degree longitude x 1 kilometer altitude grid and delivered to NASA as electronic files.

Description: This report describes the development of a three-dimensional database of aircraft fuel burn and emissions (fuel burned, NOx, CO, and hydrocarbons) from projected fleets of high speed civil transports (HSCTs) on a universal airline network. Inventories for 500 and 1000 HSCT fleets, as well as the concurrent subsonic fleets, were calculated. The HSCT scenarios are calculated using the NASA technology concept airplane (TCA) and update an earlier report. These emissions inventories are available for use by atmospheric scientists conducting the Atmospheric Effects of Stratospheric Aircraft (AESA) modeling studies. Fuel burned and emissions of nitrogen oxides (NOx as NO2), carbon monoxide, and hydrocarbons have been calculated on a 1 degree latitude x 1 degree longitude x 1 kilometer pressure altitude grid and delivered to NASA as electronic files.

Description: The report describes development of databases estimating aircraft engine exhaust emissions for the years 1976 and 1984 from global operations of Military, Charter, historic Soviet and Chinese, Unreported Domestic traffic, and General Aviation (GA). These databases were developed under the National Aeronautics and Space Administration's (NASA) Advanced Subsonic Assessment (AST). McDonnell Douglas Corporation's (MDC), now part of the Boeing Company has previously estimated engine exhaust emissions' databases for the baseline year of 1992 and a 2015 forecast year scenario. Since their original creation, (Ward, 1994 and Metwally, 1995) revised technology algorithms have been developed. Additionally, GA databases have been created and all past NIDC emission inventories have been updated to reflect the new technology algorithms. Revised data (Baughcum, 1996 and Baughcum, 1997) for the scheduled inventories have been used in this report to provide a comparison of the total aviation emission forecasts from various components. Global results of two historic years (1976 and 1984), a baseline year (1992) and a forecast year (2015) are presented. Since engine emissions are directly related to fuel usage, an overview of individual aviation annual global fuel use for each inventory component is also given in this report.

Description: A 3-dimensional model was developed to support mechanistic studies. The model solves the global primitive equations in isentropic coordinates, which directly characterize diabatic processes forcing the Brewer-Dobson circulation of the middle atmosphere. It's numerical formulation is based on Hough harmonics, which partition horizontal motion into its rotational and divergent components. These computational features, along with others, enable 3D integrations to be performed practically on RISC computer architecture, on which they can be iterated to support mechanistic studies. The model conserves potential vorticity quite accurately under adiabatic conditions. Forced by observed tropospheric structure, in which integrations are anchored, the model generates a diabatic circulation that is consistent with satellite observations of tracer behavior and diabatic cooling rates. The model includes a basic but fairly complete treatment of gas-phase photochemistry that represents some 20 chemical species and 50 governing reactions with diurnally-varying shortwave absorption. The model thus provides a reliable framework to study transport and underlying diabatic processes, which can then be compared against chemical and dynamical structure observed and in GCM integrations. Integrations with the Langley GCM were performed to diagnose feedback between simulated convection and the tropical circulation. These were studied in relation to tropospheric properties controlling moisture convergence and environmental conditions supporting deep convection, for comparison against mechanistic integrations of wave CISK that successfully reproduce the Madden-Julian Oscillation (MJO) of the tropical circulation. These comparisons were aimed at identifying and ultimately improving aspects of the convective simulation, with the objective of recovering a successful simulation of the MJO in the Langley GCM, behavior that should be important to budgets of upper-tropospheric water vapor and chemical species.

Description: Sporadic layers, which appear in the region from 100 km to 120 km are thought to be formed by convergent Pedersen drifts induced by altitude gradients in the zonal neutral wind. In this altitude region NO+ and 02+ are the major ions produced by photoionization and charge exchange of atmospheric and ionospheric species. The relative composition of atmospheric ions and meteoric ions in sporadic layers is important in determining their persistence, the time scales for formation, and the electrical conductivity of the layers. This rocket investigation will include a diagnosis of the neutral wind field and the electric field distribution. Coupled with ion composition measurements we will be able to expose the relevant formation mechanisms and the electrodynamic consequences of their existence. A rocket trajectory has been chosen to provide substantial horizontal sampling of the layer properties and knowledge of the horizontal gradients in composition and density are essential to determine the polarization electric fields that may be associated with ionospheric layers. The University of Texas at Dallas (UTD) is responsible for designing, building, and operating the ion mass spectrometers included on these rockets. The following provides a summary of the UTD accomplishments in the second year of the project as well as a description of the plans for the third year's activities. The UTD mass spectrometer acronym has been coined as PRIMS for Puerto Rico Ion Mass Spectrometer.

Description: Since its' launch on August 10, 1992, the TOPEX/Poseidon satellite hs successfully observed the earth's ocean circulation using a combination of precision orbit determination (POD) and dual-frequency radar altimetry.

Description: We study cases of abrupt dayside auroral brightenings and very fast auroral propagation using the POLAR UV imaging data.

Description: The effects of decade-to-century (Dec-Cen) variations in total solar irradiance (TSI) on global mean surface temperature Ts during the pre-Pinatubo instrumental era (1854-1991) are studied by using two different proxies for TSI and a simplified version of the IPCC climate model.

Description: The Acoustic Doppler Current Profiler (ADCP) combined with accurate navigation provides absolute current velocities which include information from all the frequencies which have a dynamical presence in the ocean.

Description: Both surface and internal mass redistribution causes geocenter variations. The global surface mass redistribution can be measured or modeled and is considered as the primary contributor to the geocenter variations on seasonal time scale.

Description: Ulysses is a joint NASA/ESA mission that is exporing a three-dimensional structure of the heliosphere via an orbit over the poles of the Sun. Since its inception, Ulysses, both in NASA and ESA, has had a public and science outreach program; it is only within the last three years that the project has begun to expand its outreach activities to include education.

Description: In late 1997, the Ulysses spacecraft completed its first orbit around the Sun, observing the properties of the heliosphere at all latitudes between 80 degrees South and 80 degrees North. Because the mission occurred during a period of near-minimum solar activity, the configuration of the solar wind and interplanetary magnetic field were particularly simple, thus allowing confident comparisons between the properties of the polar corona observed by instruments of the Spartan and SOHO spacecraft and the resulting properties of the solar wind.

Description: Simultaneous in situ measurements of the long-lived trace gases N(sub 2)O and CH(sub 4) were made with a tunable diode laser spectrometer (ALIAS II) aboard the Observations from the Middle Stratosphere (OMS) balloon platform from New Mexico, Alaska, and Brazil during 1996 and 1997.

Description: Sources such as atmospheric or buried explosions and shallow earthquakes producing strong vertical ground displacements are known to produce infrasonic pressure waves in the atmosphere. Because of the coupling between neutral particles and electrons at ionospheric altitudes, these acoustic waves induce variations of the ionospheric electron density. The Global Positioning System provides a way of directly measuring the Total Electron Content in the ionosphere and, therefore. of detecting such perturbations in the upper atmosphere. In this work, we demonstrate the capabilities of the GPS technique to detect ionospheric perturbations caused by the January 17. 1994, M (sub w) =6.7, Northridge earthquake and the STS-58 Space Shuttle ascent. In both cases, we observe a perturbation of the ionospheric electron density lasting for about 30 m, with periods less than 10 m. The perturbation is complex and shows two sub-events separated by about 15 m. The phase velocities and waveform characteristics of the two sub-events lead us to interpret the first arrival as the direct propagation of 2 free wave, followed by oscillatory guided waves propagating along horizontal atmospheric interfaces at 120 km altitude and below.

Description: The contribution of VLBI to geophysics and geodesy arises from its ability to measure distances between stations in a network and to determine the orientation of stations in the network as well as the orientation of the network with respect to the external reference frame of extragalactic radio objects. Integrating nearly two decades of observations provides useful information about station positions and velocities and the orientation of the Earth, but the complications of the real world and the limitations of observing, modeling and analysis prevent recovery of all effects. Of the factors that limit the accuracy of seemingly straightforward geodetic parameters, the neutral propagation medium has been subject to the greatest scrutiny, but the treatment of the mapping function, the wet component and spatial/temporal inhomogeneities is still improving. These affect both the terrestrial scale and consistency over time. The modeling of non-secular site motions (tides and loading) has increased in sophistication, but there are some differences between the models and the observations. VLBI antennas are massive objects, so their behavior is quite unlike GPS monuments, but antenna deformations add some (generally) unmodeled signal. Radio sources used in geodetic VLBI observations are selected for strength and (relative) absence of structure, but apparent changes in position can leak into geodetic parameters. A linear rate of change of baseline or site parameters is the simplest model and its error improves with time span. However, in most cases the VLBI data distribution is insufficient to look for real non-linear behavior that might affect the average rate. A few sites have multiple VLBI antennas, and some show small differences in rate. VLBI intrinsically measures relative positions and velocities, but individual site positions and velocities are generally more useful. The creation of the VLBI terrestrial reference frame, which transforms relative information into individual results, is an empirical process that has intrinsic errors. While UT1 is uniquely measured by VLBI, the geographical distribution and availability of VLBI stations, especially in the southern hemisphere, and the consistency of the VLBI terrestrial reference frame may limit the accuracy of Earth orientation measurements. The effects of particular error sources on geodetic and geophysical parameters derived from VLBI data will be illustrated.

Description: The plasmasphere is a toroidal region around the Earth and is filled with cold dense plasma in which the magnetic field lines are mainly closed. The plasmasphere has been studied since the discovery of the plasmaspause by Carpenter in 19xx. Since that time an impressive array of satellites and rockets have been brought to bear on studies of the region. It is a region in which the science seems to be considered as mature. We will explore the known aspects of the plasmasphere, those things that appear to be known well enough to be modeled with some degree of confidence. The areas that are not as confidently modeled, mostly on the outer boundary, aspects of refilling, convection around to the dayside, and duskside phenomenon, will be examined with the object of showing the areas of research in the plasmasphere that need further investment of resources.

Description: The high-latitude ion velocities measured by the Thermal Ion Dynamics Experiment (TIDE) instrument on the Polar spacecraft will be examined in relation to magnetospheric circulation. TIDE derives ion velocities from moments of measured distribution functions. Hydrogen and oxygen ions are E X B drifting in the polar cap and cleft regions with a speed of about 5-20 km/s at apogee (approximately 9 Re) and a speed of 1-2 km/s at perigee (approximately 1.8 Re). At perigee 0+ is typically seen flowing down in the polar cap and outflowing from the cleft. At the transition from downflowing to upflowing there is also seen a reversal in the ion convection. The convection at perigee is consistent with standard ionospheric convection models for given Interplanetary Magnetic Field (IMF) conditions. Convection at high altitude (approximately 8.9 Re) polar regions has not been studied very much since there have not been many satellites in this region. Unlike previous missions to this region TIDE in conjunction the Plasma Source Instrument (PSI) can measure ions with as low an energy as several electron Volts. The outflowing ions observed by TIDE at apogee are believed to be important to the overall circulation of the magnetosphere. The convection of these outflowing ions at apogee will be related to the IMF. This study tries to answer the question of how the IMF response of the convection influences the overall circulation of the magnetosphere.

Description: The significant results obtained with support of this grant include the following: (1) Using VLBI data in combination with other geodetical, geophysical, and geological data to bound the present rotation of the Colorado Plateau, and to evaluate to its implications for the kinematics and seismogenic potential of the western half of the conterminous U.S. (2) Determining realistic estimates of uncertainties for VLBI data and then applying the data and uncertainties to obtain an upper bound on the integral of deformation within the "stable interior" of the North American and other plates and thus to place an upper bound on the seismogenic potential within these regions. (3) Combining VLBI data with other geodetic, geophysical, and geologic data to estimate the motion of coastal California in a frame of reference attached to the Sierra Nevada-Great Valley microplate. This analysis has provided new insights into the kinematic boundary conditions that may control or at least strongly influence the locations of asperities that rupture in great earthquakes along the San Andreas transform system. (4) Determining a global tectonic model from VLBI geodetic data that combines the estimation of plate angular velocities with individual site linear velocities where tectonically appropriate. and (5) Investigation of the some of the outstanding problems defined by the work leading to global plate motion model NUVEL-1. These problems, such as the motion between the Pacific and North American plates and between west Africa and east Africa, are focused on regions where the seismogenic potential may be greater than implied by published plate tectonic models.

Description: On May 29, 1996, under steady strong northward IMF and high solar wind dynamic pressure conditions both Polar and Interball cross field lines that pass through the northern cusp and apparently close to the post-cusp reconnection site. The magnetopause current observed by Interball consists of two quite distinct layers, an inner broad current that is quite turbulent and another current that is quite abrupt and quiet. Polar also crosses current layers, similar to the Interball inner one. These observations support a model in which cusp field lines experience essentially stochastic behavior but on average provide topological connection between the cusp and magnetosheath.

Description: In this study we investigate the transport of H+ ions that made up the complex ion distribution function observed by the Geotail spacecraft at 0740 UT on November 24, 1996. This ion distribution function, observed by Geotail at approximately 20 R(sub E) downtail, was used to initialize a time-dependent large-scale kinetic (LSK) calculation of the trajectories of 75,000 ions forward in time. Time-dependent magnetic and electric fields were obtained from a global magnetohydrodynamic (MHD) simulation of the magnetosphere and its interaction with the solar wind and the interplanetary magnetic field (IMF) as observed during the interval of the observation of the distribution function. Our calculations indicate that the particles observed by Geotail were scattered across the equatorial plane by the multiple interactions with the current sheet and then convected sunward. They were energized by the dawn-dusk electric field during their transport from Geotail location and ultimately were lost at the ionospheric boundary or into the magnetopause.

Description: The primary goal of this project is to investigate two aspects of upper stratosphere and mesosphere circulation dynamics using data from the Upper Atmosphere Research Satellite's (UARS) Microwave Limb Sounder (MLS) and Improved Stratospheric and Mesospheric Sounder (ISAMS) instruments. The first aspect involves observational analyses and comparison with model predictions for the polar night "4-day wave". The second project focus is the distribution of carbon monoxide (CO) in the middle atmosphere and circulation features that can be deduced from analyses of such a tracer and comparison with atmospheric circulation models.

Description: This grant covered the transition period between the ISTP/GGS main mission MAMI work (NAGW5-1097) to the extended mission of GGS/SolaxMax. During this period we concentrated on the following subjects: (1) Ionospheric modeling and in particular the coupling of the ionosphere and the magnetosphere require knowledge of the ionospheric conductances. We have an algorithm to derive these conductances from POLAR UVI or POLAR VIS images. The image analysis requires carefully calibrated measurements, separation of auroral brightness from scattered solar EUV and emissions caused by photoelectrons, and modeling of the auroral conditions, Under this grant we evaluated the various sources of uncertainty that contribute to the conductance calculation by comparing different model approaches and we conducted sensitivity studies. A joint modeling effort between our ionospheric models and a global MHD magnetospheric model with S. Slinker and J. Fedder lead to improved parameterizations for the MHD model. The predictions of that new model were validated using data from ground based instruments, POLAR UVI images, and DMSP satellites. These findings were presented at the GEM meeting in Snowmass, and at the Huntsville modeling workshop in Guntersville. A manuscript for publication is in preparation. A study to evaluate the thermospheric effects of large substorms and storms was con- ducted with the TIE-GCM model. Of particular interest are compositional changes of the thermosphere in the southern hemisphere. Southern hemisphere effects are stronger than those on the northern hemisphere because of the larger offset of the magnetic dipole from the rotational axis. These results were presented at the Huntsville modeling workshop in Guntersville.

Description: The results of the studies carried out under the grant are described in detail in the following publications: Atmospheric Response in Aurora experiment: Observations of E and F region neutral winds in a region of postmidnight diffuse aurora; Walterscheid, Modified geostrophy in the thermosphere, ARIA II neutral flywheel-driven field-aligned currents in the postmidnight sector of the auroral oval: A case study; Observations of unstable atmospheric shear layers in the lower E region in the post-midnight auroral oval; and Depletion of oxygen in aurora: Evidence for a local mechanism. The objective of the Atmospheric Response in Aurora (ARIA) experiment was to measure the response of the E-region neutral flow to substorm activity in the post-midnight sector of the auroral oval. To this end, neutral wind profiles were obtained from a series of four sets of rocket flights using the chemical release technique. The measurements covered conditions ranging from quiet to disturbed. A consistent feature of the observed winds in disturbed conditions is the presence of an E-region jet located between approximately 110 and 120 km altitude. The profiles presented here show that the neutral wind speeds near 110-120 km altitude increase and the hodographs become more elongated or linear in response to higher magnetic activity levels. The wind speeds decrease and the hodographs become more circular as the activity level decreases. The presence of a wind maximum implies shears both above and below. The observations consistently show larger shears below the peak and smaller shears above. In fact, the bottom-side shears during disturbed conditions have Richardson numbers close to or below the critical value of 0.25 suggesting that the flow is highly unstable in that height range. The instability is expected to produce waves, turbulence, and enhanced diffusivity.

Description: Magnetic merging at Earth's magnetopause produces distinct mixtures of ions and electrons as well as signatures in their distribution functions. High resolution measurements allow for the separation of the different distributions and quantification of their characteristics. This provides details on the temporal and spatial nature of the merging site and the resulting history of the merged fields. The event of May 29, 1996 resulted in remote observations of the effects of reconnection on both magnetosheath and magnetosphere populations for a period of approximately three hours. Three-dimensional ion distributions obtained by the Thermal Ion Dynamics Experiment on the Polar spacecraft show that field lines threading the spacecraft's location in the northern cusp region contained a mix of D-shaped ions from the magnetosheath and accelerated magnetospheric ions both moving parallel to the local magnetic field. This mix of ions resulted from transmission of magnetosheath ions across the magnetopause at speeds greater than the de-Hoffman-Teller speed and the reflection of cold, slow-moving plasmasphere-like ions at the magnetopause. These observations are used to conclude that these field lines were connected to the ionosphere in the northern hemisphere and, southward of the spacecraft, the interplanetary magnetic field and crossed the magnetopause in the equatorial region southward of the spacecraft.

Description: We have used GPS at 10 stations in southern Alaska with three epochs of measurements to estimate short-term uplift rates. A number of great earthquakes as well as recent large earthquakes characterize the seismicity of the region this century. To reliably estimate uplift rates from GPS data, numerical models that included both the slip distribution in recent large earthquakes and the general slab geometry were constructed.

Description: A new approach is presented to quantify upper-level moisture transport from geostationary satellite data. Daily time sequences of GOES-7 water vapor imagery were used to produce estimates of winds and water vapor mixing ratio in the upper-troposphere sensed by the 6.7-microns water vapor channel. The winds and mixing ratio values were gridded and then combined to produce a parameter called the Water Vapor Transport Index (WVTI) which represents the magnitude of the two dimensional transport of water vapor in the upper troposphere. Daily grids of WVTI, meridional moisture transport, mixing ratio, pressure and other associated parameters were averaged to produce monthly fields for June, July and August of 1987 and 1988 over the Americas and surrounding oceanic regions. The WVTI was used to compare upper-tropospheric moisture transport between the summers of 1987 and 1988, contrasting the latter part of the 1986/87 El Nino event and the La Nina period of 1988. A similar product derived from the NCEP reanalysis was used to help validate the index and understand interannual variability of moisture transport from the modeling perspective. Both datasets showed large upper-level water vapor transport associated with synoptic features over the Americas and with outflow from tropical convective systems. Minimal transport occurred over tropical and subtropical high pressure regions where winds were light. Index values from NCEP were 2-3 times larger than that determined from the GOES satellite This difference resulted from an over estimate of the zonal winds and upper-tropospheric moisture in the reanalysis model. A comparison of the satellite-derived monthly values between the summers of 1987 and 1988 provided some insight into the impact of the ENSO event on upper-level moisture and its transport during the period. During July 1987, a large portion of the tropics in the eastern Pacific Ocean and Caribbean Sea was dominated by strong vapor transport in excess of 5 g/kg/ms with relatively small amounts in the other months. In contrast, JJA 1988 showed an opposite trend with July 1988 being less dominated by tropical water vapor transport. The meridional transport of upper-level water vapor indicated significant poleward transport from the tropics to the higher latitudes. This transport favored the Southern Hemisphere with large transport occurring south of the ITCZ which extended across the eastern Pacific and northern South America. Zonally-averaged monthly transport values were shown to provide a simple way to quantify the monthly and interannual changes in water vapor transport. Zonally-averaged WVTI values peaked in the Southern Hemisphere subtropics during both Austral winters. In the tropics a single more pronounced peak located over the equator and south latitudes occurred in 1988 as opposed to a dual peak in 1987. The second peak around 20 deg.N latitude is consistent with findings of others where upper-tropospheric winds were noted to be stronger in this region during warm ENSO events. Zonally-averaged meridional transport was southward for all summer months and stronger in 1988. The asymmetric nature of the zonally-averaged meridional transport (more southerly water vapor transport) was enhanced during JJA 1998 thus indicating a stronger upper-level branch of the Hadley circulation during this notable strong La Nina period.

Description: We have been studying how substorms work by examining the global polar Ultraviolet Imager (UVI) images in correlation with observations from the ground, interplanetary space and the geomagnetic tail between 10-20 earth radii. One of the objectives of our study is to better understand the connection among many complex phenomena going on close to Earth and those in the distant plasma sheet. We have studied, for example, how the aurora[ and polar cap boundaries at different local times behave in relation to variations observed in the solar wind and plasma sheet during substorms. Preliminary results indicate that the polar cap and auroral oval boundaries expand and contract in a complicated but systematic way. These variations are correlated to solar wind parameters, and thinning and recovery phenomena in the plasma sheet. These results will be presented and interpreted in terms of directly driven and/or unloading substorm processes.