ALBERT

Description: We have investigated the interaction of Io, Jupiter's innermost Galilean satellite, with the Io plasma torus. The interaction of Io with the plasma surrounding it has been a subject of interest for almost 30 years, dating from the discovery by Bigg (1964) that radio emissions from the Jovian magnetosphere are controlled by Io's position. Since that time, both ground-based and spacecraft observations have shown that Io is a unique satellite that influences the Jovian magnetosphere in important ways. In particular, material from Io is a major source of plasma for the magnetosphere, and the energy that this plasma harnesses from Jupiter's co-rotating magnetic field is an important power source for the magnetosphere. It is apparent that the local interaction of the torus plasma with Io plays a key role in the formation, composition, and energetics of the Io torus; the interaction is also highly nonlinear. We have modeled this interaction using time-dependent three-dimensional magnetohydrodynamic (MHD) simulations. During this past year, we have used NASA support to develop a new MHD code to study the interaction. As part of the Galileo spacecraft's recent successful insertion into orbit around Jupiter, the spacecraft passed within 900 km of Io's surface. Our calculations have focused on using Galileo particles and fields data to examine a question that was not resolved by the Voyager observations: Does Io have an intrinsic magnetic field? In this progress summary, we describe our efforts on this problem to date.

Description: This paper reports a new balance for the measurement of three components of force - lift, drag and pitching moment - in impulsively starting flows which have a duration of about one millisecond. The basics of the design of the balance are presented and results of tests on a 15 deg semi-angle cone set at incidence in the T4 shock tunnel are compared with predictions. These results indicate that the prototype balance performs well for a 1.9 kg, 220 mm long model. Also presented are results from initial bench tests of another application of the deconvolution force balance to the measurement of thrust produced by a 2D scramjet nozzle.

Description: The eruption of Mt. Pinatubo in June 1991 produced the largest enhancement of stratospheric aerosol loading ever observed by lidar over Hampton, Virginia. Low altitude layers (less than 20 km) were the first to arrive over Hampton in early August, the result of transport associated with a tropospheric anticyclonic cell over North America. The maximum peak scattering ratio, 34 at 22.4 km, and the maximum stratospheric integrated backscatter of 0.0053 sr(exp -1), both at 694 nm, observed since the eruption were measured on February 20, 1992. After decreasing during the spring and summer of 1992, the aerosol burden increased significantly during the winter of 1992-93, evidence of a poleward winter transport from the equatorial reservoir. Over the period from February 1992 to February 1994, the stratospheric aerosol loading decreased with an average 1/e decay time of 10.1 months. The vertical distribution, intensity, and transport of Pinatubo aerosols over this site are described and compared with similar measurements after El Chichon.

Description: A preliminary comparison of the GEOS-1 (Goddard Earth Observing System) data assimilation system convective cloud mass fluxes with fluxes from a cloud-resolving model (the Goddard Cumulus Ensemble Model, GCE) is reported. A squall line case study (10-11 June 1985 Oklahoma PRESTORM episode) is the basis of the comparison. Regional (central U. S.) monthly total convective mass flux for June 1985 from GEOS-1 compares favorably with estimates from a statistical/dynamical approach using GCE simulations and satellite-derived cloud observations. The GEOS-1 convective mass fluxes produce reasonable estimates of monthly-averaged regional convective venting of CO from the boundary layer at least in an urban-influenced continental region, suggesting that they can be used in tracer transport simulations.

Description: T-matrix computations of light scattering by polydispersions of randomly oriented nonspherical aerosols and Mie computations for equivalent spheres are compared. Findings show that even moderate nonsphericity results in suubstantial errors in the retrieved aerosol optical thickness if satellite reflectance measurements are analyzed using Mie theory. On the other hand, the use of Mie theory for nonspherical aerosols produces negligible errors in the computation of albedo and flux related quantities, provided that the aerosol size distribution and optical thickness are known beforehand. The first result can be explained by large nonspherical-spherical differences in scattering phase function, while the second result follows from small nonspherical-spherical differences in single-scattering albedo and asymmetry parameter. No cancellation of errors occurs if one consistently uses Mie theory in the retrieval algorithm and then in computing the albedo for the retrieved aerosol optical thickness.

Description: Azimuthal asymmetries in the atmospheric refractive index can lead to errors in estimated vertical and horizontal station coordinates. Daily average gradient effects can be as large as 50 mm of delay at a 7 deg elevation. To model gradients, the constrained estimation of gradient paramters was added to the standard VLBI solution procedure. Here the analysis of two sets of data is summarized: the set of all geodetic VLBI experiments from 1990-1993 and a series of 12 state-of-the-art R&D experiments run on consecutive days in January 1994. In both cases, when the gradient parameters are estimated, the overall fit of the geodetic solution is improved at greater than the 99% confidence level. Repeatabilities of baseline lengths ranging up to 11,000 km are improved by 1 to 8 mm in a root-sum-square sense. This varies from about 20% to 40% of the total baseline length scatter without gradient modeling for the 1990-1993 series and 40% to 50% for the January series. Gradients estimated independently for each day as a piecewise linear function are mostly continuous from day to day within their formal uncertainties.

Description: Particles leving the neutral sheet in the distant magnetotail at times display adiabatic trajectory sequences characterized by an inflection toward the equator and subsequent mirroring in its vicinity. We demonstrate that this low-latitude mirroring results primarily from a centrifugal deceleration due to the fast direction-changing E x B drift. This effect which we refer to as 'centrifugal trapping' appears both in guiding centre and full particle treatments. It thus does not directly relate to nonadiabatic motion. However, pitch angle scattering due to nonadiabatic neutral sheet interaction does play a role in reducing the parallel speed of the particles. We show that centrifugal trapping is an important mechanism for the confinement of the slowest (typically below the equatorial E x B drift speed) plasma sheet populations to the midplane vicinity.

Description: Estimates of the mass of dust suspended in the Martian atmosphere are derived from global and regional 9-micrometer opacity maps produced from Viking Infrared Thermal Mapper data. During the peak of the 1977b storm, a total dust mass of approximately 4.3 x 10(exp 14) g was suspended, equivalent to 4.3 x 10(exp -4) g/sq cm, or a layer 1.4 micrometers thick. During a local dust storm near Solis Planum at L(sub s) 227 deg, approximately 1.3 x 10(exp 13) g of dust were lofted, equal to about a 6-micrometer layer in that vicinity.

Description: The U.S. upper Midwest was subjected to severe flooding during the summer of 1993. Heavy rainfall in the Mississippi River basin from April through July caused flooding on many Midwest rivers, including the Mississippi, Illinois, Missouri, and Kansas Rivers. The flood crest of 15.1 m at St. Louis, Missouri, on 1 August 1993 was the highest ever measured, surpassing the previous record of 13.2 m set on 28 April 1973. Damage estimates include at least 47 flood-related deaths and a total damage cost of $12 billion. Remotely sensed imagery of severe flooding in the U.S. Midwest was obtained under cloud-free skies on 29 July 1993 by the MODIS (Moderate Resolution Imaging Spectroradiometer) Airborne Simulator (MAS). The MAS is a newly developed scanning spectrometer with 50 spectral bands in the wavelength range 0.55-14.3 micrometers. Estimation of the total flooded area in the MAS scene acquired near St. Louis was accomplished by comparing the MAS scene to a Landsat-5 thematic mapper (TM) scene of the same area acquired on 14 April 1984 in nonflood conditions. For comparison, the MAS band centered at 0.94 micrometers and the TM band centered at 1.65 micrometers were selected because of the high contrast seen in these bands between land and water-covered surfaces. An estimate of the area covered by water in the MAS and TM scenes was obtained by developing land/water brightness thresholds from histograms of the MAS and TM digital image data. Afetr applying the thresholds, the difference between the area covered by water in the MAS and TM scenes, and hence the flooded area in the MAS scene, was found to be about 396 sq km, or about 153 square miles.

Description: Shortwave radiative fluxes that reach the earth's surface are key factors that influence atmospheric and oceanic circulations as well as surface climate. Yet, information on these fluxes is meager. Surface site data are generally available from only a limited number of observing stations over land. Much less is known about the large-scale variability of the shortwave radiative fluxes over the oceans, which cover most of the globe. Recognizing the need to produce global-scale fields of such fluxes for use in climate research, the World Climate Research Program has initiated activities that led to the establishment of the Surface Radiation Budget Climatology Project with the ultimate goal to determine various components of the surface radiation budget from satellite data. In this paper, the first global products that resulted from this activity are described. Monthly and daily data on a 280-km grid scale are available. Samples of climate parameters obtainable from the dataset are presented. Emphasis is given to validation and limitations of the results. For most of the globe, satellite estimates have bias values between +/- 20 W/sq m and root mean square (rms) values are around 25 W/sq m. There are specific regions with much larger uncertainties however.

Description: The deployment of a space-based Doppler lidar would provide information that is fundamental to advancing the understanding and prediction of weather and climate. This paper reviews the concepts of wind measurement by Doppler lidar, highlights the results of some observing system simulation experiments with lidar winds, and discusses the important advances in earth system science anticipated with lidar winds. Observing system simulation experiments, conducted using two different general circulation models, have shown (1) that there is a significant improvement in the forecast accuracy over the Southern Hemisphere and tropical oceans resulting from the assimilation of simulated satellite wind data, and (2) that wind data are significantly more effective than temperature or moisture data in controlling analysis error. Because accurate wind observations are currently almost entirely unavailable for the vast majority of tropical cyclones worldwide, lidar winds have the potential to substan- tially improve tropical cyclone forecasts. Similarly, to improve water vapor flux divergence calculations, a direct measure of the ageostrophic wind is needed since the present level of uncer- tainty cannot be reduced with better temperature and moisture soundings alone.

Description: We model an infrared outburst on Io as being due to a large, erupting lava flow which increased its area at a rate of 1.5 x 10(exp 5)/sq m and cooled from 1225 to 555 K over the 2.583-hr period of observation. The inferred effusion rate of 3 x 10(exp 5) cu m/sec for this eruption is very high, but is not unprece- dented on the Earth and is similar to the high eruption rates suggested for early lunar volcanism. Eruptions occur approxi- mately 6% of the time on Io. These eruptions provide ample resurfacing to explain Io's lack of impact craters. We suggest that the large total radiometric heat flow, 10(exp 14) W, and the size and temperature distribution of the thermal anomalies (McEwen et al. 1992; Veeder et al. 1994) can be accounted for by a series of silicate lava flows in various stages of cooling. We propose that the whole suite of Io's currently observed thermal anomalies was produced by multiple, high-eruptive-rate silicate flows within the past century.

Description: The 'satellite-gauge model' (SGM) technique is described for combining precipitation estimates from microwave satellite data, infrared satellite data, rain gauge analyses, and numerical weather prediction models into improved estimates of global precipitation. Throughout, monthly estimates on a 2.5 degrees x 2.5 degrees lat-long grid are employed. First, a multisatellite product is developed using a combination of low-orbit microwave and geosynchronous-orbit infrared data in the latitude range 40 degrees N - 40 degrees S (the adjusted geosynchronous precipitation index) and low-orbit microwave data alone at higher latitudes. Then the rain gauge analysis is brougth in, weighting each field by its inverse relative error variance to produce a nearly global, observationally based precipitation estimate. To produce a complete global estimate, the numerical model results are used to fill data voids in the combined satellite-gauge estimate. Our sequential approach to combining estimates allows a user to select the multisatellite estimate, the satellite-gauge estimate, or the full SGM estimate (observationally based estimates plus the model information). The primary limitation in the method is imperfections in the estimation of relative error for the individual fields. The SGM results for one year of data (July 1987 to June 1988) show important differences from the individual estimates, including model estimates as well as climatological estimates. In general, the SGM results are drier in the subtropics than the model and climatological results, reflecting the relatively dry microwave estimates that dominate the SGM in oceanic regions.

Description: Measurements of superthermal electron fluxes in the solar wind indicate that field lines within coronal mass ejections, CMEs, near and beyond 1 AU are normally connected to the Sun at both ends. However, on occasion some field lines embedded deep within CMEs appear to be connected to the Sun at only one end. Here we propose an explanation for how such field lines arise in terms of 3-dimensional reconnection close to the Sun. Such reconnection also provides a natural explanation for the flux rope topology characteristic of many CMEs as well as the coronal loops formed during long-duration, solar X-ray events. Our consideration of the field topologies resulting from 3-dimensional reconnection indicates that field lines within and near CMEs may on occasion be connected to the outer heliosphere at both ends.

Description: We investigate the orbital dynamics of small dust particles generated via the continuous micrometeoroid bombardment of the Martian moons. In addition to Mar's oblateness, we also consider the radiation pressure perturbation that is complicated by the planet's eccentric orbit and tilted rotational axis. Considering the production rates and the lifetimes of dust grains, we show that particles from Deimos with radii of about 15 micrometers are expected to dominate the population of a permanently present and tilted dust torus. This torus has an estimated peak number density of approximately equals 5 x 10(exp -12)/cu cm and an optical depth of approximately equals 4 x 10(exp -8).

Description: The lower stratospheric variability of equatorial water vapor, measured by the Microwave Limb Sounder (MLS), follows an annual cycle modulated by the quasi-biennial oscillation. At levels higher in the stratosphere, water vapor measurements exhibit a semi-annual oscillatory signal with the largest amplitudes at 2.2 and 1hPa. Zonal-mean cross sections of MLS water vapor are consistent with previous satellite measurements from the limb infrared monitor of the stratosphere (LIMS) and the stratospheric Aerosol and Gas Experiment 2 (SAGE 2) instruments in that they show water vapor increasing upwards and the polewards from a well defined minimum in the tropics. The minimum values vary in height between the retrieved 46 and 22hPa pressure levels.

Description: We examined 11 cases when the interplanetary magnetic field (IMF) was intensely northward (greater than 10 nT) for long durations of time (greater than 3 hours), to quantitatively determine an uppler limit on the efficiency of solar wind energy injection into the magnetosphere. We have specifically selected these large B(sub N) events to minimize the effects of magnetic reconnection. Many of these cases occurred during intervals of high-speed streams associated with coronal mass ejections when viscous interaction effects might be at a maximum. It is found that the typical efficiency of solar wind energy injection into the magnetosphere is 1.0 x 10(exp -3) to 4.0 x 10(exp -3), 100 to 30 times less efficient than during periods of intense southward IMFs. Other energy sinks not included in these numbers are discussed. Estimates of their magnitudes are provided.

Description: The valley network channels on the heavily cratered ancient surface of Mars suggest the presence of liquid water approximately 3.8 Gyr ago. However, the implied warm climate is difficult to explain in the context of the standard solar model, even allowing for the maximum CO2 greenhouse heating. In this paper we investigate the astronomical and planetary implications of a nonstandard solar model in which the zero-age, main-sequence Sun had a mass of 1.05 +/- 0.02 Solar Mass. The excess mass was subsequently lost in a solar wind during the first 1.2(-0.2, +0.4) Gyr of the Sun's main sequence phase. The implied mass-loss rate of 4(+3, -2) x 10(exp -11) M/yr, or about 10(exp 3) x that of the current Sun, may be detectable in several nearby young solar type stars.

Description: Hellas basin on Mars has been the site of volcanism, tectonism, and modification by fluvial, mass-wasting, and eolian processes over its more than 4-b.y. existence. Our detailed geologic mapping and related studies have resulted in the following new interpretations. The asymmetric distribution of highland massifs and other structures that define the uplifted basin rim suggest a formation of the basin by the impact of a low-angle bolide having a trajectory heading S60E. During the Late Noachian, the basin was infilled, perhaps by lava flows, that were sufficiently thick (greater than 1 km) to produce wrinkle ridges on the fill material and extensional faulting along the west rim of the basin. At about the same time, deposits buried northern Malea Planum, which are interpreted to be pyroclastic flows from Amphitrites and Peneus Paterae on the basis of their degraded morphology, topology, and the application of a previous model for pyroclastic volcanism on Mars. Peneus forms a distinctive caldera structure that indicates eruption of massive volumes of magma, whereas Amphitrites is a less distinct circular feature surrounded by a broad, low, dissected shield that suggests generally smaller volume eruptions. During the Early Hesperian, an approximately 1-to 2km-thick sequence of primarily fined-grained, eolian material was deposited on the floor of Hellas basin. Subsequently, the deposit was deeply eroded, except where armored by crater ejecta, and it retreated as much as 200-300 km along its western margin, leaving behind pedestal craters and knobby outliers of the deposit. Local debris flows within the deposit attest to concentrations of groundwater, perhaps in part brought in by outflow floods along the east rim of the basin. These floods may have deposited approximately 100-200m of sediment, subduing wrinkle ridges in the eastern part of the basin floor. During the Late Hesperian and Amazonian, eolian mantles were emplaced on the basin rim and floor and surrounding highlands. Their subsequent erosion resulted in pitted and etched plains and crater fill, irregular mesas, and pedestal craters. Local evidence occurs for the possible former presence of ground ice or ice sheets approximately 100 km across; however, we disagree with a hypothesis that suggest that the entire south rim and much of the floor of Hellas have been glaciated. Orientations of dune fields and yardangs in lower parts of Hellas basin follow directions of the strongest winds predicted by a recently published general circulation model (GCM). Transient frost and dust splotches in the region are, by contrast, related to the GCM prediction for the season in which the images they appear in were taken.

Description: Visible and near-IR refectivity, Moessbauer, and X ray diffraction data were obtained on powders of impact melt rock from the Manicouagan Impact Crater located in Quebec, Canada. The iron mineralogy is dominated by pyroxene for the least oxidized samples and by hematite for the most oxidized samples. Phyllosilicate (smectite) contents up to approximately 15 wt % were found in some heavily oxidized samples. Nanophase hematite and/or paramagnetic ferric iron is observed in all samples. No hydrous ferric oxides (e.g., goethite, lepidocrocite, and ferrihydrite) were detected, which implies the alteration occurred above 250 C. Oxidative alteration is thought to have occurred predominantly during late-stage crystallization and subsolidus cooling of the impact melt by invasion of oxidizing vapors and/or solutions while the impact melt rocks were still hot. The near-IR band minimum correlated with the extent of aleration Fe(3+)/Fe(sub tot) and ranged from approximately 1000 nm (high-Ca pyroxene) to approximately 850 nm (bulk, well-crystalline hematite) for least and most oxidized samples, respectively. Intermediate band positions (900-920 nm) are attributed to low-Ca pyroxene and/or a composite band from hematite-pyroxene assemblages. Manicouagan data are consistent with previous assignments of hematite and pyroxene to the approximately 850 and approximately 1000nm bands observed in Martian reflectivity spectra. Manicouagan data also show that possible assignments for intermediate band positions (900-920 nm) in Martian spectra are pyroxene and/or hematite-pyroxene assemblages. By analogy with impact melt sheets and in agreement with observables for Mars, oxidative alteration of Martian impact melt sheets above 250 C and subsequent erosion could produce rocks and soils with variable proportions of hematite (both bulk and nanophase), pyroxene, and phyllosilicates as iron-bearing mineralogies. If this process is dominant, these phases on Mars were formed rapidly at relativly high temperatures on a sporadic basis throughout the history of the planet. The Manicouagan samples also show that this mineralogical diversity can be accomplished at constant chemical composition, which is also indicated for Mars from the analyses of soil at the two Viking landing sites.

Description: The subject of this paper is a self-consistent, magnetohydrodynamic numerical realization for the Earth's magnetosphere which is in a quasi-steady dynamic equilibrium for a due northward interplanetary magnetic field (IMF). Although a few hours of steady northward IMF are required for this asymptotic state to be set up, it should still be of considerable theoretical interest because it constitutes a 'ground state' for the solar wind-magnetosphere interaction. Moreover, particular features of this ground state magnetosphere should be observable even under less extreme solar wind conditions. Certain characteristics of this magnetosphere, namely, NBZ Birkeland currents, four-cell ionospheric convection, a relatively weak cross-polar potential, and a prominent flow boundary layer, are widely expected. Other characteristics, such as no open tail lobes, no Earth-connected magnetic flux beyond 155 R(sub E) downstream, magnetic merging in a closed topology at the cusps, and a 'tadpole' shaped magnetospheric boundary, might not be expected. In this paper, we will present the evidence for this unusual but interesting magnetospheric equilibrium. We will also discuss our present understanding of this singular state.

Description: The stability of collisionless tearing modes is examined in the presence of curvature drift resonances and the trapped particle effects. A kinetic description for both electrons and ions is employed to investigate the stability of a two-dimensional equilibrium model. The main features of the study are to treat the ion dynamics properly by incorporating effects associated with particle trajectories in the tail fields and to include the linear coupling of trapped particle modes. Generalized dispersion relations are derived in several parameter regimes by considering two important sublayers of the reconnecting region. For a typical choice of parameters appropriate to the current sheet region, we demonstrate that localized tearing modes driven by ion curvature drift resonance effects are excited in the current sheet region with growth time of the order of a few seconds. Also, we examine nonlocal characteristics of tearing modes driven by curvature effects and show that modes growing in a fraction of a second arise when mode widths are larger than the current sheet width. Further, we show that trapped particle effects, in an interesting frequency regime, significantly enhance the growth rate of the tearing mode. The relevance of this theory for substorm onset phase and other features of the substorms is briefly discussed.

Description: The second pulsating aurora rocket (PULSAUR 2) rocket was launched into a pulsating aurora on 9 February 1994 and carried a broad range of instruments in order to perform a study of this type of aurora. The rocket measurements were complemented with a set of ground-based measurements. The particle measurements performed in the rocket are related to the ground-based optical measurements performed along the rocket trajectory. It is found that the high energy electrons are largely in phase with the measured luminosity. The EISCAT measurements carried out during the flight are reviewed. The PULSAUR 2 campaign is described. Results concerning the auroral conditions and the particle measurements are presented.

Description: The high-speed correction factor to the O(+)-O collision frequency, resulting from drift velocities between ions and neutrals, is calculated by solving the integral expression in this factor both numerically and analytically. Although the analytic solution is valid for either small or large drift velocities between ions and neutrals, for temperatures of interest and all drift velocities considered, agreement is found between analytic and detailed numerical integration results within less than 1% error. Let T(sub r) designate the average of the ion and neutral temperatures in K, and u = nu(sub d)/alpha, where nu(sub d) is the relative drift velocity in cm/s, and alpha = 4.56 x 10(exp 3) square root of T(sub r) cm/s is the thermal velocity of the O(+)-O system. Then, as u ranges from 0 to 2, the correction factor multiplying the collision frequency increases monotonically from 1 to about 1.5. An interesting result emerging from this calculation is that the correction factor for temperatures of aeronomical interest is to a good approximation independent of the temperature, depending only on the scaled velocity u.

Description: This paper describes an efficient Monte Carlo algorithm for choosing a new direction of a photon after a scattering interaction. The algorithm chooses a scattering angle by linear interpolation in a table of the inverse cumulative scattering probability. A Legendre expansion of the phase function makes it easy to apply Clenshaw's algorithm to build the interpolation table. The points in the table are close enough together that linear interpolation is accurate. With a table of 100,000 entries, we can keep the absolute and relative errors in matching the probability distribution below 10(exp -5).

Description: The absolute rate constant for the reaction O((3)P) + HOBr has been measured between T = 233K and 423K using the discharge-flow kinetic technique coupled to mass spectrometric detection. The value of the rate coefficient at room temperature is (2.5 +/- 0.6) x 10(exp -11)cu cm/molecule/s and the derived Arrhenius expression is (1.4 +/- 0.5) x 10(exp -10) exp((-430 +/- 260)/T)cu cm/molecule/s. From these rate data the atmospheric lifetime of HOBr with respect to reaction with O((3)P) is about 0.6h at z = 25 km which is comparable to the photolysis lifetime based on recent measurements of the UV cross section for HOBr. Implications for HOBr loss in the stratosphere have been tested using a 1D photochemical box model. With the inclusion of the rate parameters and products for the O + HOBr reaction, calculated concentration profiles of BrO increase by up to 33% around z = 35 km. This result indicates that the inclusion of the O + HOBr reaction in global atmospheric chemistry models may have an impact on bromine partitioning in the middle atmosphere.

Description: The relation between the solar wind input to the magetosphere, VB(sub South), and the auroral geomagnetic index AL is modeled with two linear moving-average filtering methods: linear prediction filters and a driven harmonic oscillator in the form of an electric circuit. Although the response of the three-parameter oscillator is simpler than the filter's, the methods yield similar linear timescales and values of the prediction-observation correlation and the prediction Chi(exp 2). Further the filter responses obtained by the two methods are similar in their long-term features. In these aspects the circuit model is equivalent to linear prediction filtering. This poses the question of uniqueness and proper interpretation of detailed features of the filters such as response peaks. Finally, the variation of timescales and filter responses with the AL activity level is discussed.

Description: Empirical data-based models of the magnetosphereic magnetic field have been widely used during recent years. However, the existing models (Tsyganenko, 1987, 1989a) have three serious deficiencies: (1) an unstable de facto magnetopause, (2) a crude parametrization by the K(sub p) index, and (3) inaccuracies in the equatorial magnetotail B(sub z) values. This paper describes a new approach to the problem; the essential new features are (1) a realistic shape and size of the magnetopause, based on fits to a large number of observed crossing (allowing a parametrization by the solar wind pressure), (2) fully controlled shielding of the magnetic field produced by all magnetospheric current systems, (3) new flexible representations for the tail and ring currents, and (4) a new directional criterion for fitting the model field to spacecraft data, providing improved accuracy for field line mapping. Results are presented from initial efforts to create models assembled from these modules and calibrated against spacecraft data sets.

Description: We have investigated the middle atmospheric response to the 27-day and 11-yr solar UV flux variations at low to middle latitudes using a two-dimensional photochemical model. The model reproduced most features of the observed 27-day sensitivity and phase lag of the profile ozone response in the upper stratosphere and lower mesosphere, with a maximum sensitivity of +0.51% per 1% change in 205 nm flux. The model also reproduced the observed transition to a negative phase lag above 2 mb, reflecting the increasing importance with height of the solar modulated HO(x) chemistry on the ozone response above 45 km. The model revealed the general anti-correlation of ozone and solar UV at 65-75 km, and simulated strong UV responses of water vapor and HO(x) species in the mesosphere. Consistent with previous 1D model studies, the observed upper mesospheric positive ozone response averaged over +/- 40 was simulated only when the model water vapor concentrations above 75 km were significantly reduced relative to current observations. In agreement with observations, the model computed a low to middle latitude total ozone phase lag of +3 days and a sensitivity of +0.077% per 1% change in 205 nm flux for the 27-day solar variation, and a total ozone sensitivity of +0.27% for the 11-yr solar cycle. This factor of 3 sensitivity difference is indicative of the photochemical time constant for ozone in the lower stratosphere which is comparable to the 27-day solar rotation period but is much shorter than the 11-yr solar cycle.

Description: The accuracy of Geosat satellite altimetry over the Greenland ice sheet is evaluated by comparing the measured heights to radar elevations from the airborne Greenland Aerogeophysics Project. At the center of the ice sheet where the ice surface is nearly level, surface comparisons show a fit at the 1 to 3 m level as expected, but even at moderately sloping ice regions (0.3 deg-0.6 deg), satellite altimetry mean errors in the range of 10 to 35 m are observed. These errors are found for slope-corrected and waveform-retracked data, so most previous accuracy estimates of current satellite altimetry ice sheet elevations in regions of slopping or undulating ice appear to be too optimistic.

Description: We report here on a number of examples of anomalous enhancements of eastward electric fields near sunrise in the equatorial ionospheric F-region. These examples were selected from the data base of the equatorial satellite, San Marco D (1988), which measured ionospheric electric fields during a period of solar minimum. The eastward electric fields reported correspond to vertical plasma drifts. The examples studied here are similar in signature and polarity to the pre-reversal electric field enhancements seen near sunset from ground-based radar systems. The morphology of these sunrise events, which are observed on about 14% of the morning-side satellite passes, are studied as a function of local zonal velocity, magnetic activity, geographic longitude and altitude. The nine events studied occur at locations where the zonal plasma flow is generally measured to be eastward, but reducing as a function of local time and at satellite longitudes where the magnetic declination has the opposite polarity as the declination of the sunrise terminator.

Description: A diagnostic analysis of the VVP (volume velocity processing) retrieval method is presented, with emphasis on understanding the technique as a linear, multivariate regression. Similarities and differences to the velocity-azimuth display and extended velocity-azimuth display retrieval techniques are discussed, using this framework. Conventional regression diagnostics are then employed to quantitatively determine situations in which the VVP technique is likely to fail. An algorithm for preparation and analysis of a robust VVP retrieval is developed and applied to synthetic and actual datasets with high temporal and spatial resolution. A fundamental (but quantifiable) limitation to some forms of VVP analysis is inadequate sampling dispersion in the n space of the multivariate regression, manifest as a collinearity between the basis functions of some fitted parameters. Such collinearity may be present either in the definition of these basis functions or in their realization in a given sampling configuration. This nonorthogonality may cause numerical instability, variance inflation (decrease in robustness), and increased sensitivity to bias from neglected wind components. It is shown that these effects prevent the application of VVP to small azimuthal sectors of data. The behavior of the VVP regression is further diagnosed over a wide range of sampling constraints, and reasonable sector limits are established.

Description: Analysis of the version 16 Halogen Occultation Experiment (HALOE) CH4 data shows that this long-lived trace gas is well correlated with potential vorticity (PV) computed from National Meteorological Center balanced winds. Analyzing late September and October 1992 data, we show that very low CH4 values are confined to the interior of a vortex edge defined by the maximum gradient in PV. The CH4 and HF time tendency is used to estimate the descent rate in the Antarctic vortex. After removing a component of the trend correlated with the HALOE sampling pattern, we compute the lower stratosphere vertical descent rates and net heaing rates in the spring Antarctic vortex. Our computations of the spring Antarctic vortex heating rates give -0.5 to -0.1 K/day. Over the winter season, the overall lower stratospheric descent rate averages about 1.8-1.5 km/month. These computations are in line with radiative transfer estimates of the heating and descent rate. The HALOE data thus appear to be consistent with the picture of an isolated lower stratospheric Antarctic vortex.

Description: Using delta C-13 measurements in atmospheric CO2 from a cooperative global air sampling network, we determined the partitioning of the net uptake of CO2 between ocean and land as a function of latitude and time. The majority of delta C-13 measurements were made at the Institute of Arctic and Alpine Research (INSTAAR) of the University of Colorado. We perform an inverse deconvolution of both CO2 and delta C-13 observations, using a two-dimensional model of atmospheric transport. Also, the discrimination against C-13 by plant photosynthesis, as a function of latitude and time, is calculated from global runs of the simple biosphere (SiB) model. Uncertainty due to the longitudinal structure of the data, which is not represented by the model, is studied through a bootstrap analysis by adding and omitting measurement sites. The resulting error estimates for our inferred sources and sinks are of the order of 1 GTC (1 GTC = 10(exp 15) gC). Such error bars do not reflect potential systematic errors arising from our estimates of the isotopic disequilibria between the atmosphere and the oceans and biosphere, which are estimated in a separate sensitivity analysis. With respect to global totals for 1992 we found that 3.2 GTC of carbon dissolved into the ocean and that 1.5 GTC were sequestered by land ecosystems. Northern hemisphere ocean gyres north of 15 deg N absorbed 2.7 GTC. The equatorial oceans between 10 deg S and 10 deg N were a net source to the atmosphere of 0.9 GTC. We obtained a sink of 1.6 GTC in southern ocean gyres south of 20 deg S, although the deconvolution is poorly constrained by sparse data coverage at high southern latitudes. The seasonal uptake of CO2 in the northern gyres appears to be correlated with a bloom of phytoplankton in surface waters. On land, northern temperate and boreal ecosystems between 35 deg N and 65 deg N were found to be a major sink of CO2 in 1992, as large as 3.5 GTC. Northern tropical ecosystems (equator-30 deg N) appear to be a net source to the source to the atmosphere of 2 GTC which could reflect biomass burning. A small sink, 0.3 GTC, was inferred for southern tropical ecosystems (30 deg S-equator).

Description: Global sets of surface radiation budget (SRB) have been obtained from satellite programs. These satellite-based estimates need validation with ground-truth observations. This study validates the estimates of monthly mean surface insolation contained in two satellite-based SRB datasets with the surface measurements made at worldwide radiation stations from the Global Energy Balance Archive (GEBA). One dataset was developed from the Earth Radiation Budget Experiment (ERBE) using the algorithm of Li et al. (ERBE/SRB), and the other from the International Satellite Cloud Climatology Project (ISCCP) using the algorithm of Pinker and Laszlo and that of Staylor (GEWEX/SRB). Since the ERBE/SRB data contain the surface net solar radiation only, the values of surface insolation were derived by making use of the surface albedo data contained GEWEX/SRB product. The resulting surface insolation has a bias error near zero and a root-mean-square error (RMSE) between 8 and 28 W/sq m. The RMSE is mainly associated with poor representation of surface observations within a grid cell. When the number of surface observations are sufficient, the random error is estimated to be about 5 W/sq m with present satellite-based estimates. In addition to demonstrating the strength of the retrieving method, the small random error demonstrates how well the ERBE derives from the monthly mean fluxes at the top of the atmosphere (TOA). A larger scatter is found for the comparison of transmissivity than for that of insolation. Month to month comparison of insolation reveals a weak seasonal trend in bias error with an amplitude of about 3 W/sq m. As for the insolation data from the GEWEX/SRB, larger bias errors of 5-10 W/sq m are evident with stronger seasonal trends and almost identical RMSEs.

Description: A nonlinear filtering method is introduced for the study of the solar wind -- magnetosphere coupling and related to earlier linear techniques. The filters are derived from the magnetospheric state, a representation of the magnetospheric conditions in terms of a few global variables, here the auroral electrojet indices. The filters also couple to the input, a representation of the solar wind variables, here the rectified electric field. Filter-based iterative prediction of the indices has been obtained for up to 20 hours. The prediction is stable with respect to perturbations in the initial magnetospheric state; these decrease exponentially at the rate of 30/min. The performance of the method is examined for a wide range of parameters and is superior to that of other linear and nonlinear techniques. In the magnetospheric state representation the coupling is modeled as a small number of nonlinear equations under a time-dependent input.

Description: Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Description: Participants in this workshop, which convened in Venice, Italy, 6-8 May 1993, met to consider the current state of climate monitoring programs and instrumentation for the purpose of climatological prediction on short-term (seasonal to interannual) timescales. Data quality and coverage requirements for definition of oceanographic heat and momentum fluxes, scales of inter- and intra-annual variability, and land-ocean-atmosphere exchange processes were examined. Advantages and disadvantages of earth-based and spaceborne monitoring systems were considered, as were the structures for future monitoring networks, research programs, and modeling studies.

Description: The work described here makes it possible to identify anomalous wind behavior such as the nighttime meridional wind abatements that occur at F-region heights. A new analysis technique uses a simple empirical wind model to simulate measurements of 'normal' winds (as measured by the Neutral Atmosphere and Temperature Experiment (NATE) that flew on the Atmosphere Explorer-E (AE-E)) to highlight anomalous wind measurements made by the satellite while in circular orbits at 270-290 km altitude. Our approach is based on the recognition that the 'in orbit' wind variation must show the combined effects of the diurnal wind variation as seen from the ground with the latitude variation of the satellite orbit. For the data period 77250-78035 examined thus far, the wind abatement always occurred with a corresponding pressure or temperature maximum, and was detected on 12 out of the 36 nights with data. This study has revealed that the wind abatement occur only during or shortly after increases in solar extreme ultraviolet (EUV) flux, as indicated by daily radio flux measurements. In the past, nighttime wind reversals at mid-latitudes have been associated with increased geomagnetic activity. This study indicates that intensified solar EUV heating may be responsible for anomalous thermospheric nighttime winds at mid-latitudes.

Description: The concept of the well-known Langley plot technique, used for the calibration of ground-based instruments, has been generalized for application to satellite instruments. In polar regions, near summer solstice, the solar backscattered ultraviolet (SBUV) instrument on the Nimbus 7 satellite samples the same ozone field at widely different solar zenith angles. These measurements are compared to assess the long-term drift in the instrument calibration. Although the technique provides only a relative wavelength-to-wavelength calibration, it can be combined with existing techniques to determine the drift of the instrument at any wavelength. Using this technique, we have generated a 12-year data set of ozone vertical profiles from SBUV with an estimated accuracy of +/- 5% at 1 mbar and +/- 2% at 10 mbar (95% confidence) over 12 years. Since the method is insensitive to true changes in the atmospheric ozone profile, it can also be used to compare the calibrations of similar SBUV instruments launched without temporal overlap.

Description: An automated scheme to characterize precipitation echoes within small windows in the radar field is presented and applied to previously subjectively classified tropical rain cloud systems near Darwin, Australia. The classification parameters are (a) E(sub e), effective efficiency, as determined by cloud-top and cloud-base water vapor saturation mixing ratios; (b) BBF, brightband fraction, as determined by the fraction of the radar echo area in which the maximal reflectivity occurs within +/- 1.5 km of the 0 C isotherm level; and (c) del(sub r) Z, radial reflectivity gradients (dB/km). These classification criteria were applied to tropical rain cloud systems near Darwin, Australia, and to winter convective rain cloud systems in Israel. Both sets of measurements were made with nearly identical networks of C-band radars and rain gauge networks. The results of the application of these objective classification criteria to several independently predetermined rain regimes in Darwin have shown that better organized rain systems have smaller del(sub r) Z and larger BBF. Similarly, smaller del(sub r)Z and larger BBF were also observed from maritime rain cloud systems, as compared to continental rain cloud systems with the same degree of organization. Continental rain cloud system, regardless of their degree of organization, have larger depths, as expressed by E(sub e). The rainfall analyses presented in this study are based exclusively on rain gauge measurements, while radar information was used only to classify the individual gauge measurements.

Description: Many problems in geophysical and astrophysical convection systems are characterized by fast rotation and spherical shell geometry. The combined effects of Coriolis forces and spherical shell geometry produce a unique spatial symmetry for the convection pattern in a rapidly rotating spherical shell. In this paper, we first discuss the general spatial symmetries for rotating spherical shell convection. A special model, a spherical shell heated from below, is then used to illustrate how and when the spatial symmetries are broken. Symmetry breaking occurs via a sequence of spatial transitions from the primary conducting state to the complex multiple-layered columnar structure. It is argued that, because of the dominant effects of rotation, the sequence of spatial transitions identified from this particular model is likely to be generally valid. Applications of the spatial symmetry breaking to planetary convection problems are also discussed.

Description: Initial observations of a newly documented type of optical emission above thunderstorms are reported. 'Blue jets', or narrowly collimated beams of blue light that appear to propagate upwards from the tops of thunderstorms, were recorded on B/W and color video cameras for the first time during the Sprites94 aircraft campaign, June-July, 1994. The jets appear to propagate upward at speeds of about 100 km/s and reach terminal altitudes of 40-50 m. Fifty six examples were recorded during a 22 minute interval during a storm over Arkansas. We examine some possible mechanisms, but have no satisfactory theory of this phenomenon.

Description: It is shown from flux transfer event (FTE) occurrence statistics, observed as a function of MLT by the ISEE satellites, that recent 2-dimensional analytic theories of the effects of pulsed Petschek reconnection predict FTEs to contribute between 50 and 200 kV to the total reconnection voltage when the magnetosheath field points southward. The upper limit (200 kV) allows the possibility that FTEs provide all the antisunward transport of open field lines into the tail lobe. This range is compared with the voltages associated with series of FTEs signatures, as inferred from ground-based observations, which are in the range 10-60 kV. We conclude that the contribution could sometimes be made by a series of single, large events; however, the voltage is often likely to be contributed by several FTEs at different MLT.

Description: New investigations of the core ion motion within high-latitude topside ionosphere and near-Earth magnetosphere, using data from the Dynamics Explorer (DE) retarding ion mass spectrometer (RIMS), reveal the existence of significant regions of downward moving O(+). The occurences of downgoing versus upgoing O(+) are not clearly separable in terms of either polar zenith angle or Kp but are well distinquished by the direction of the z component of the interplanetary magnetic field (IMF). On the average, down flow dominates when IMFB2 less than O, while upflow dominates for IMFB2 greater than 0. Combining cross-field convection velocities derived from an empirical convection electric field model with the observed parallel velocities yields a two dimensional ion velocity field. This velocity field is consistent with a senario which has O(+) of cusp/cleft and auroral zone origin concvecting into the polar cap and, because of the dominance of gravitational energy over the upward kinetic energy, falling back into the inonsphere. This provides additional confirmation of the results of studies of the cleft ion fountain. Estimates of the flux of O(+) in the upflowing and downflowing regions for Lambda greater than or equal to 60 deg give a total upflow of approximately 6 X 10(exp 25) ions/sec for IMFB2 greater than 0 and total upflow and downflow of approximately 4 X 10 (exp 25) ions/sec and 1 x 10(exp 25) ions/sec, respectively, for IMFB less than 0. In all cases the magnitude of the dayside outflow is consistent with previous work on upwelling ions. While the magnitudes vary for high and low Kp the ratios of upward to downward flow are roughly the same at approximately 1.7. The downflowing O(+) shows a correlation with the magnitude of the outflow of light ions in the same region but the cause and effect of this relationship is not distinquishable. hable.

Description: We have developed the numerical algorithm for the computation of transient viscoelastic responses in the time domain for a radially stratified Earth model. Stratifications in both the elastic parameters and the viscosity profile have been considered. The particular viscosity profile employed has a viscosity maximum with a constrast of O(100) in the mid lower mantle. The distribution of relaxation times reveals the presence of a continuous spectrum situated between O(100) and O(exp 4) years. The principal mode is embedded within this continuous spectrum. From this initial-value approach we have found that for the low degree harmonics the non-modal contributions are comparable to the modal contributions. For this viscosity model the differences between the time-domain and normal-mode results are found to decrease strongly with increasing angular order. These calculations also show that a time-dependent effective relaxation time can be defined, which can be bounded by the relaxation times of the principal modes.

Description: The effects on double-probe electric field measurements induced by electron density and temperature gradients are investigated. We show that on some occasions such gradients may lead to marked spurious electric fields if the probes are assumed to lie at the same probe potential with repect to the plasma. The use of a proper bias current will decrease the magnitude of such an error. When the probes are near the plasma potential, the magnitude of these error signals, delta Epsilon, can vary as delta Epsilon approx. T(sub e)(Delta n(sub e)/n(sub e)) + 0.5 Delta T(sub e), where T(sub e) is the electron temperature, Delta n(sub e)/ n(sub e) the relative electron density variation between the two sensors, and Delta T(sub e) the electron temperature difference between the two sensors. This not only implies that the error signals will increase linearly with the density variations but also that such signatures grow with Delta T(sub e) i.e., such effects are 10 times larger in a 10-eV plasma than in a 1-eV plasma. This type of error is independent of the probe separation distance provided the gradient scale length is much larger than the distance. The largest errors occur when the probes are near to the plasma potential. During the crossing of a small structure (e.g, a double layer) the error signal appears as a bipolar signature. Our analysis shows that errors in double-probe measurements caused by plasma gradients are not significant at large scale (much greater than 1 km) plasma boundaries, and may only be important in cases where small-scale (less than 1 km), internal gradient structures exist. Bias currents tailored for each plasma parameter regime (i.e., variable bias current) would improve the double-probe response to gradient effects considerably.

Description: The evolution of the volcanic debris plume originating from the June 1991 eruption of Mt. Pinatubo has been monitored since its genesis using a ground-based backscatter lidar facility sited at the Jet Propulsion Laboratory (JPL). Both absolute and relative pre- and post-Pinatubo backscatter observations are in accord with Mie scattering projections based on measured aerosol particle size distributions reported in the literature. The post-Pinatubo column-integrated backscatter coefficient peaked approximately 400 days after the eruption, and the observed upper boundary of the aerosol column subsided at a rate of approximately 200 m/mon.

Description: Using Global Positioning System (GPS) receivers, we reoccupied several leveling benchmarks on the Kenai Peninsula of Alaska which had been surveyed by conventional leveling immediately following the March 27, 1964, Prince William Sound earthquake (M(sub w) = 9.3). By combining the two sets of measurements with a new, high-resolution model of the geoid in the region, we were able to determine the cumulative 1993-1964 postseismic vertical displacement. We find uplift at all of our benchmarks, relative to Seward, Alaska, a point that is stable according to tide gauge data. The maximum uplift of about 1 m occurs near the middle of the peninsula. The region of maximum uplift appears to be shifted northwest relative to the point of maximum coseismic subsidence. If we use tide gauge data at Nikishka and Seward to constrain the vertical motion, then the observed uplift has a trenchward tilt (down to the southeast) as well as an arching component. To explain the observations, we use creep-at-depth models. Most acceptable models require a fault slip of about 2.75 m, although this result is not unique. If the slip has been continuous since the 1964 earthquake, then the average slip rate is nearly 100 mm/yr, twice the plate convergence rate. Comparing the net uplift achieved in 29 years with that observed over 11 years in an adjacent region southeast of Anchorage, Alaska, we conclude that the rate of uplift is decreasing. A further decrease in the uplift rate is expected as the 29-year averaged displacement rate is about twice the plate convergence rate and therefore cannot be sustained over the entire earthquake cycle.

Description: Narrow enhancements of electron precipitation, with energy and flux well above typical values, have been observed with Dynamics Explorer 2 (DE 2) in the cusp/cleft region. The electron flux in the energy range 0.2-1 keV was 2 orders of magnitude higher in these structures than in the magnetosheath and were seen in approximaetly 80% of DE 2 cusp crossings at ionospheric altitudes. Typically, there was more than one electron structure in each cusp crossing. The position of these structures showed a systematic variation: for poleward ion dispersion (energy decreases with increasing latitude), electron structures were seen more often on the equatorial boundary of the cusp, while for equatorward ion dispersion (energy decreases with decreasing latitude), electron structures were more often seen on the poleward boundary. This suggests that the electron structures are associated with newly reconnected field lines. The electron spectra suggest that field-aligned acceleration processes could produce the electron structures, first near the boundary of the cusp/cleft during the reconnection of field lines and then in the cusp/cleft during the motion of reconnected flux tubes through the polar ionosphere.

Description: Results of a simulation study on the effects of optically thick stratospheric sulfate aerosol layers on the backscattered ultraviolet radiation (buv) in the range 256-340 nm are presented. In general, the increased Mie scattering produced by the aerosols results in an enhancement of the buv radiation. The increase is approximately linear with optical depth and strongly depends on solar zenith angle and aerosol layer altitude in relation to the ozone maximum. The effect is greatest at those wavelenghts whose contribution functions peak in the vicinity of the densest part of the aerosol layer. The aerosol induced perturbation of the buv field affects the ozone profile retrieval from space measurements by the Solar Backscatter Ultraviolet Experiment (SBUV) experiment. In the tropical stratosphere, the retrieved ozone between 25 and 45 km is underestimated as a result of increased Mie scattering. On the other hand, an algorithm related effect causes the retrieved ozone below 25 km to be overestimated by an amount similar to the stratospheric deficit.

Description: A permanent Global Positioning System (GPS) receiver at Casa Diablo Hot Springs, Long Valley Caldera, California was installed in January, 1993, and has operated almost continuously since then. The data have been transmitted daily to the Jet Propulsion Laboratory (JPL) for routine analysis with data from the Fiducial Laboratories for an International Natural sciences Network (FLINN) by the JPL FLINN analysis center. Results from these analyses have been used to interpret the on going deformation at Long Valley, with data excluded from periods when the antenna was covered under 2.5 meters of snow and from some periods when Anti Spoofing was enforced on the GPS signal. The remaining time series suggests that uplift of the resurgent dome of Long Valley Caldera during 1993 has been 2.5 +/- 1.1 cm/yr and horizontal motion has been 3.0 +/- 0.7 cm/yr at S53W in a no-net-rotation global reference frame, or 1.5 +/- 0.7 cm/yr at S14W relative to the Sierra Nevada block. These rates are consistent with uplift predicted from frequent horizontal strain measurements. Spectral analysis of the observations suggests that tidal forcing of the magma chamber is not a source of the variability in the 3 dimensional station location. These results suggest that remotely operated, continuously recording GPS receivers could prove to be a reliable tool for volcanic monitoring throughout the world.

Description: A comparison is made of the vertical distribution of high-level cloud tops derived from the Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements and from the International Satellite Cloud Climatology Project (ISCCP) for all Julys and Januarys in 1985 to 1990. The results suggest that ISCCP overestimates the pressure of high-level clouds by up to 50-150 mbar, particularly at low latitudes. This is caused by the frequent presence of clouds with diffuse tops (greater than 50% time when cloudy events are observed). The averaged vertical extent of the diffuse top is about 1.5 km. At midlatitudes where the SAGE II and ISCCP cloud top pressure agree best, clouds with distinct tops reach a maximum relative proportion of the total level cloud amount (about 30-40%), and diffuse-topped clouds are reduced to their minimum (30-40%). The ISCCP-defined cloud top pressure should be regarded not as the material physical height of the clouds but as the level which emits the same infrared radiance as observed. SAGE II and ISCCP cloud top pressures agree for clouds with distinct tops. There is also an indication that the cloud top pressures of optically thin clouds not overlying thicker clouds are poorly estimated by ISCCP at middle latitudes. The average vertical extent of these thin clouds is about 2.5 km.

Description: Global high-level clouds identified in Stratospheric Aerosol and Gas Experiment II (SAGE II) occultation measurements for January and July in the period 1985 to 1990 are compared with near-nadir-looking observations from the International Satellite Cloud Climatology Project (ISCCP). Global and zonal mean high-level cloud amounts from the two data sets agree very well, if clouds with layer extinction coefficients of less than 0.008/km at 1.02 micrometers wavelength are removed from the SAGE II results and all detected clouds are interpreted to have an average horizontal size of about 75 km along the 200 km transimission path length of the SAGE II observations. The SAGE II results are much more sensitive to variations of assumed cloud size than to variations of detection threshold. The geographical distribution of cloud fractions shows good agreement, but systematic regional differences also indicate that the average cloud size varies somewhat among different climate regimes. The more sensitive SAGE II results show that about one third of all high-level clouds are missed by ISCCP but that these clouds have very low optical thicknesses (less than 0.1 at 0.6 micrometers wavelength). SAGE II sampling error in monthly zonal cloud fraction is shown to produce no bias, to be less than the intraseasonal natural variability, but to be comparable with the natural variability at longer time scales.

Description: Enhancements in the fluxes of relativistic electrons trapped within the Earth's magnetosphere have been measured by the high-energy particle spectrometer, part of the particle environment monitor on the upper atmosphere research satellite (UARS). The largest increase in the electron fluxes with energies greater than 1 MeV observed on UARS from October 1991 through July 1994 was in early May 1992. The fluxes of trapped electrons in the drift loss cone and locally precipitating electrons showed differing buildup and decay rates as a function of invariant latitude. Increases of more than 2 orders of magnitude were observed in drift loss cone fluxes at magnetic latitudes of 40 deg-66 deg and in precipitating fluxes from 48 deg to 66 deg. The energy flux contained in the most intense local precipitation observed was approximately 0.1 erg/sq cm/s, entering the atmosphere and creating up to 1000 ion pairs/cu cm/s at 55-km altitude. The daily averaged energy flux from directly precipitating electrons with energies greater than 1 MeV deposited greater than 10(exp 20) erg/d worldwide into the atmosphere for the period May 12-21, 1992, producing greater than 10(exp 31) odd nitrogen molecules below 60-km altitude.

Description: Application of the window probability matching method to radar and rain gauge data that have been objectively classified into different rain types resulted in distinctly different Z(sub e)-R relationships for the various classifications. These classification parameters, in addition to the range from the radar, are (a) the horizontal radial reflectivity gradients (dB/km); (b) the cloud depth, as scaled by the effective efficiency; (c) the brightband fraction within the radar field window; and (d) the height of the freezing level. Combining physical parameters to identify the type of precipitation and statistical relations most appropriate to the precipitation types results in considerable improvement of both point and areal rainfall measurements. A limiting factor in the assessment of the improved accuracy is the inherent variance between the true rain intensity at the radar measured volume and the rain intensity at the mouth of the rain guage. Therefore, a very dense rain gauge network is required to validate most of the suggested realized improvement. A rather small sample size is required to achieve a stable Z(sub e)-R relationship (standard deviation of 15% of R for a given Z(sub e)) -- about 200 mm of rainfall accumulated in all guages combined for each classification.

Description: A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.

Description: Radiation resulting from interaction between the effluent cloud of a space shuttle thruster and the ambient atmosphere was observed with a spectograph aboard the shutttle. The spectral measurements were made between 400 and 800 nm with a resolutoion of 3 nm. The primary emissions are identified as NO2, HNO, O(1)D, and O(1)S. These are the first observations od O(1)S emission in the shuttle plume. These data are compared with the previous measurements, and possible excitation mechanisms are discussed. The results are also compared with a Monte Carlo simulation of thruster plume-atmosphere interaction radiation.

Description: We have studied the implications of recent trends in the annual mean and the amplitude of the annual harmonic of ozone in the upper stratosphere from the 15 years of the combined data from the Nimbus-7 SBUV and the NOAA-11 SBUV/2 instruments. This was done in the context of the GSFC 2D model predictions of these trends which are based on plausible scenarios of anthropogenic Cly increase in the atmosphere. The comparison of the observed and model-estimated annual mean ozone trends show some similarity in their latitude and altitude characteristics. Both the model and data show a maximum ozone decrease of -6 to -10 % per decade at high latitudes in the upper stratosphere. However, there are also significant differences between the observed and computed trends which may be related to both the model uncertainty and the uncertainty in correcting for the long term instrument drift. The observations also suggest a decrease of 10-25 % per decade in the annual amplitude of ozone at 2 mb between 40 deg - 60 deg in both hemispheres, with a relatively larger interannual variability in the northern hemisphere. These values are in general agreement with the model predictions and thus provide additional support in favor of the chlorine induced changes in ozone in the upper stratosphere.

Description: BrO measured from the NASA ER-2 during Airborne Arctic Stratospheric Expedition (AASE) II exhibited a mean value (for 20-minute averages) of 5.4 +/- 0.3 pptv, with a standard deviation of 3.1 pptv. Ratios of BrO to available inorganic bromine (Br(sub y)) show only slight increases in polar regions relative to midlatitudes. A comparison between observed latitudinal and diurnal variations of this same ratio and that calculated by photochemical models shows reasonable agreement in behavior, but significant discrepancies in magnitude. It is unclear whether this difference is due to errors in measurements, models or both.

Description: Measurements of stratospheric hydrofluoric acid (HF) have been made by the JPL MkIV interferometer during high-altitude balloon flights. Infrared solar absorption spectra were acquired near 35 deg N at altitudes between local tropopause and 38 km. Volume mixing ratio profiles of HF derived from 4 flights (1990-93), in conjunction with simultaneously observed N2O profiles, indicate an average rate of HF increase of (5.5 +/- 0.3)% per year, in agreement with time-dependent, two-dimensional model simulations (6% per year) and ATMOS measurements.

Description: Upper Atmosphere Research Satellite (UARS) Microwave Limb Sounder (MLS) measurements of lower stratospheric ClO during 1992-93 and 1993-94 Arctic winters are presented. Enhanced ClO in the 1992-93 winter was first observed in early December, and extensively during February when temperatures were continually low enough for polar stratospheric cloud (PSCs). Sporadic episodes of enhanced ClO were observed for most of the 1993-94 winter as minimum temperatures hovered near the PSC threshold, with largest ClO amounts occurring in early March after a sudden deep cooling in late February.

Description: We report multi-instrument observations during an isolated substorm on 17 October 1989. The European Incoherent Scatter (EISCAT) radar operated in the SP-UK-POLI mode measuring ionospheric convection at latitudes 71 deg Lambda - 78 deg Lambda. Sub-Auroral Magnetometer Network (SAMNET) and the EISCAT Magnetometer Cross provide information on the timing of substorm expansion phase onset and subsequent intensifications, as well as the location of the field aligned and ionospheric currents associated with the substorm current wedge. Interplanetary Monitoring Platform-8 (IMP-8) magnetic field data are also included. Evidence of a substorm growth phase is provided by the equatorward motion of a flow reversal boundary across the EISCAT radar field of view at 2130 MLT, following a southward turning of the interplanetary magnetic field (IMF). We infer that the polar cap expanded as a result of the addition of open magnetic flux in the tail lobes during this interval. The flow reversal boundary, which is a lower limit to the polar cap boundary, reached an invariant latitude equatorward of 71 deg Lambda by the time of the expansion phase onset. We conclude that the substorm onset region in the ionosphere, defined by the westward electrojet, mapped to a part of the tail radially earthward of the boundary between open and closed magnetic flux, the distant neutral line. Thus the substorm was not initiated at the distant neutral line, although there is evidence that it remained active during the expansion phase.

Description: The impact of the extraterrestrial object that formed the Chicxulub crater in the northwestern Yucatan peninsula of Mexico is the leading suspect for the extinction of the dinosaurs. This article reports on a Planetary Society expedition to Albion Island in the Rio Hondo region of Belize to investigate evidence supporting the impact theory.

Description: The thermal-infrared (longwave) emission from a vegetated terrain is generally anisotropic, i.e., the emission temperature varies with the view direction. If a directional measurement of temperature is considered to be equal to the effective temperature of the hemispheric emission, then the estimate of the latter can be significantly in error. The view-direction (zenith angle theta(sub eq) at which the emission equivalence does hold is determined in our modeling study. In a two-temperature field-of-view (soil and plants), theta(sub eq) falls in a narrow range depending on plant density and canopy architecture. Theta(sub eq) does not depend on soil and (uniform) plant temperatures nor on their ratio, even though the pattern of emission vs. the view direction depends crucially on this ratio. For a sparse canopy represented as thin, vertical cylindrical stalks (or vertical blades uniformly distributed in azimuth) with horizontal facets, theta(sub eq) ranges from 48 to 53 deg depending on the optical density of the vertical elements alone. When plant elements are modeled as small spheres, theta(sub eq) lies between 53 to 57 deg (for the same values of the canopy optical density). Only for horizontal leaves (a truly planophile canopy) is the temperature measured from any direction equal to the temperature of the hemispheric emission. When the emission temperature changes with optical depth within the canopy at a specified rate, theta(sub eq) depends to some extent on that rate. For practically any sparsely vegetated surface, a directional measurement at the zenith angle of 50 deg offers an appropriate evaluation of the hemispheric emission, since the error in the estimate will, at most, only slightly exceed 1% (around 4 W/sq m). Estimates of the hemispheric emission through a nadir measurement, on the other hand, can be in error in some cases by about 10%, i.e., on the order of 40 W/sq m.

Description: Direct current (DC) electric field and ion density measurements near density depletion regions (that is, equatorial plasma bubbles) are used to estimate the vertical neutral wind speed. The measured zonal electric field in a series of density depletions crossed by the San Marco D satellite at 01.47-01.52 Universal Time (UT) on 25 October 1988, can be explained if a downward neutral wind of 15-30 m/s exists. Simultaneously, the F-region plasma was moving downward at a speed of 30-50 m/s. These events appear in the local time sector of 23.00-23.15 in which strong downward neutral winds may occur. Indeed, airglow measurements suggest that downward neutral velocities of 25-50 m/s are possible at time near midnight in the equatorial F-region.

Description: An overview of the observations of backstreaming electrons in the foreshock and the mechanisms that have been proposed to explain their properties will be presented. A primary characteristic of observed foreshock electrons is that their velocity distributions are spatially structured in a systematic way depending on distance from the magnetic field line which is tangent to the shock. There are two interrelated aspects to explaining the structure of velocity distributions in the foreshock, one involving the acceleration mechanism and the other, propagation from the source to the observing point. First, the source distribution of electrons energized by the shock must be determined along the shock surface. Proposed acceleration mechanisms include magnetic mirroring of incoming solar wind particles and mechanisms involving transmission of particles through the shock. Secondly, the kinematics of observable electrons streaming away from a curved shock with an initial parallel velocity and a downstream perpendicular velocity component due to the motional electric field must be determined. This is the context in which the observations and their explanations will be reviewed.

Description: Based on the conservation of chemical elements in chemical reactions, a rule is proved that the number of boundary conditions given by densities and/or nonzero velocities should not be less than the number of chemical elements in the system, and the boundary conditions for species given by densities and velocities should include all elements in the system. Applications of this rule to Mars are considered. It is shown that the problem of the CO2-H2O chemistry in the lower and middle atmosphere of Mars, say, in the range of 0-80 km does not have a unique solution, if only CO2 and H2O densities are given at the lower boundary, and the remaining boundary conditions are fluxes. Two examples of models of this type are discussed. Two models of the photochemistry of the Martian atmosphere, with and without nitrogen chemistry, are considered. The oxygen nonthermal escape ratio of 1.2 x 10(exp 8)/cu cm/s is given at 240 km and is balanced with the total hydrogen escape rate within an uncertainty of 1% for both models. Both models fit the measured O2 and CO mixing ratios, the O3 abundance, and the O2 1.27-micrometer dayglow almost within the uncertainties of the measured values, though the model without nitrogen chemistry fits better. The importance of nitrogen chemistry in the lower and middle atmosphere of Mars depends on a fine balance between production of NO and N in the upper atmosphere which is not known within the required accuracy.

Description: Goethite-bearing samples with values of Mn(s) (Mn/(Mn+Fe) mole fraction) up to 0.206 were synthesized by precipitation from alkaline solution. Samples with Mn(s) less than or equal 0.061 were single-phase Mn-goethites: samples with higher Mn(s) values contained another Mn-bearing phase (probably jacobsite). Mn-hematites were prepared by dehydroxylation of corresponding Mn-goethites at 500 C. Orthorhombic a and b unit cell dimensions of Mn-goethites changed in a linear manner with Mn(s), but not at rates predicted by the Vegrad law. Hexagonal unit cell dimensions of Mn-hematites did not vary with Mn(s). Moessbauer parameters isomer shift (IS), quadrupole splitting (QS), and hyperfine field (B(sub hf)) were measured at 293 and 15 K. For all single-phase Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061), magnetic splitting was observed at both temperatures. At 293 K, small but systematic decreases in B(sub hf) were observed with increasing Mn substitution; IS and QS were not dependent on Mn(s). Mn substitution strongly lowered the Morin transition temperature of hematite. At 15 K, the Morin transition was not present for Mn(s) greater than 0.020(4). The saturation magnetization of Mn-goethites and Mn-hematites (Mn(s) less than or equal 0.061) was the expected zero (within error) for antiferromagnetic goethite and for hematites obtained from dehydroxylation of goethites. Mn-geothites with Mn(s) greater than 0.061 were magnetic because of the presence of strongly magnetic jacobsite. For reflectivity spectra, bands resulting from MN(3+) were centered near 454 and 596 nm for Mn-goethites and near 545 and 700 nm for Mn-hematites. There is evidence for a approximately 700 nm band in spectral data for Martian bright regions, but association of it with Mn(3+) is not a unique interpretation. Comparison of laboratory and Martian spectral data implies that Mn(s) less than 0.032 for the Mn(3+) content of Martian hematites.

Description: Ferric-iron-bearing materials play an important role in the interpretation of visible to near-IR Mars spectra, and they may play a similarly important role in the analysis of new mid-IR spacecraft spectral observations to be obtained over the next decade. We review exisiting data on mid-IR transmission spectra of ferric oxides/oxyhydroxides and present new transmission spectra for ferric-bearing materials spanning a wide range of mineralogy and crystallinity. These materials include 11 samples of well-crystallized ferric oxides (hematite, maghemite, and magnetite) and ferric oxyhydroxides (goethite, lepidocrocite). We also report the first transmission spectra for purely nanophase ferric oxide samples that have been shown to exhibit spectral similarities to Mars in the visible to near-IR and we compare these data to previous and new transmission spectra of terrestial palagonites. Most of these samples show numerous, diagnostic absorption features in the mid-IR due to Fe(3+) - 0(2-) vibrational transitions, structural and/or bound OH, and/or silicates. These data indicate that high spatial resolution, moderate spectral resolution mid-IR ground-based and spacecraft observations of Mars may be able to detect and uniquely discriminate among different ferric-iron-bearing phases on the Martian surface or in the airborne dust.

Description: Numerical simulation experiments were conducted to delineate the influence of in situ deforestation data on episodic rainfall by comparing two ensembles of five 5-day integrations performed with a recent version of the Goddard Laboratory for Atmospheres General Circulation Model (GCM) that has a simple biosphere model (SiB). The first set, called control cases, used the standard SiB vegetation cover (comprising 12 biomes) and assumed a fully forested Amazonia, while the second set, called deforestation cases, distinguished the partially deforested regions of Amazonia as savanna. Except for this difference, all other initial and prescribed boundary conditions were kept identical in both sets of integrations. The differential analyses of these five cases show the following local effects of deforestation. (1) A discernible decrease in evapotranspiration of about 0.80 mm/d (roughly 18%) that is quite robust in the averages for 1-, 2-, and 5-day forecasts. (2) A decrease in precipitation of about 1.18 mm/d (roughly 8%) that begins to emerge even in 1-2 day averages and exhibits complex evolution that extends downstream with the winds. (3) A significant decrease in the surface drag force (as a consequence of reduced surface roughness of deforested regions) that, in turn, affects the dynamical structure of moisture convergence and circulation. The surface winds increase significantly during the first day, and thereafter the increase is well maintained even in the 2- and 5-day averages.

Description: A steady-state scheme for data assimilation in the context of a single, short period (relative to a day), sun-synchronous, polar-orbiting satellite is examined. If the satellite takes observations continuously, the gains, which are the weights for blending observations and predictions together, are steady in time. For a linear system forced by random noise, the optimal steady-state gains (Wiener gains) are equivalent to those of a Kalman filter. Computing the Kalman gains increases the computational cost of the model by a large factor, but computing the Wiener gains does not. The latter are computed by iteration using prior estimates of the gains to assimilate simulated observations of one run of the model, termed 'truth' into another run termed 'prediction'. At each stage, the prediction errors form the basis for the next estimate of the gains. Steady state is achieved after three or four iterations. Further simplification is achieved by making the gains depend on longitudinal distance from the observation point, not on absolute longitude. For a single-layer primitive equation model, the scheme works well even if only the mass field is observed but not the velocity field. Although the scheme was developed for Mars Observer, it should be applicable to data retrieved from Earth atmosphere satellites, for example, UARS.

Description: Brightness temperature difference (BTD) values are calculated for selected Geostationary Operational Environmental Satellite (GOES-6) channels (3.9, 12.7 micrometer) and Advanced Very High Resolution Radiometer channels (3.7, 12.0 micrometer). Daytime and nighttime discrimination of particle size information is possible given the infrared cloud extinction optical depth and the BTD value. BTD values are presented and compared for cirrus clouds composed of equivalent ice spheres (volume, surface area) versus randomly oriented hexagonal ice crystals. The effect of the hexagonal ice crystals is to increase the magnitude of the BTD values calculated relative to equivalent ice sphere (volume, surface area) BTDs. Equivalent spheres (volume or surface area) do not do a very good job of modeling hexagonal ice crystal effects on BTDs; however, the use of composite spheres improves the simulation and offers interesting prospects. Careful consideration of the number of Legendre polynomial coefficients used to fit the scattering phase functions is crucial to realistic modeling of cirrus BTDs. Surface and view-angle effects are incorporated to provide more realistic simulation.

Description: An examination and analysis of video images of lightning, captured by the payload bay TV cameras of the space shuttle, provided a variety of examples of lightning in the stratosphere above thunderstorms. These images were obtained on several recent shuttle flights while conducting the Mesoscale Lightning Experiment (MLE). The images of stratospheric lightning illustrate the variety of filamentary and broad vertical discharges in the stratosphere that may accompany a lightning flash. A typical event is imaged as a single or multiple filament extending 30 to 40 km above a thunderstorm that is illuminated by a series of lightning strokes. Examples are found in temperate and tropical areas, over the oceans, and over the land.

Description: This paper focuses on the comparison of cloud amounts derived from an atmospheric general circulation model (AGCM), Satellite-observed clouds, and Ground-based cloud observations. Unlike Earth Radiation Budget Experiment (ERBE)-type comparisons it does not mix potential errors in the cloud amount with those in the radiation code embedded in the model. Long term cloud climatologies were used to compare global cloud amounts and regional seasonal cycles. The AGCM successfully reproduced the signatures of the warm pool and North Pacific seasonal cycle cloudiness but failed in the low stratus region off the coast of South America, a known problem for AGCMs. The data sets also reproduced the anomaly signature associated with El Nino in the warm pool region, but the model amounts were lower. Global results had a similar success rate, with the model generally producing lower total cloud compared to the satellite and in situ measurements. To compare cloud vertical distributions the cloud height may need to be validated using the corresponding radiation fields. Unfortunately there were also some large discrepancies between the two observed cloud data sets. While tremendously improved over the last decade the character of the observed cloud data sets, must be substantially enhanced before they will be useful in validating AGCMs by any but the crudest levels of comparison.

Description: A prototype of the Josephson-effect spectrum analyzer developed for the millimeter wave band is described. The measurement results for spectra obtained in the frequency band from 50 to 250 GHz are presented.

Description: This summary describes the spatial, spectral, and radiometric calibration of Thermal Infrared Multispectral Scanner (TIMS) performed at the Jet Propulsion Laboratory (JPL) Thermal Infrared Calibration Facility (TIRCAL) between May and August, 1994. The 1994 calibration of TIMS was the first to make use of the new EXABYTE (8mm helical-scan tape) recording system. With the new recorder, the TIMS data tapes may be read directly on any computer system that has an EXABYTE tape drive. We analyzed the calibration data sets using image processing procedures written in Interactive Data Language.

Description: This summary describes the 9 March 1994 Thermal Infrared Multispectral Scanner (TIMS) airborne calibration experiment conducted at Castaic Lake, California. This experiment was a collaborative effort between the TIMS and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) science teams at the Jet Propulsion Laboratory (JPL). TIMS was flown on the NASA/Ames Research Center C130 with the new retractable air fence installed in the TIMS instrument bay. The purpose of this experiment was to determine if the fence would reduce the air turbulence in the TIMS instrument bay, thereby reducing the errors in calibration caused by wind-blast cooling of the blackbody reference sources internal to TIMS. Previous experiments have indicated that the wind blast effect could cause TIMS to over-estimate surface temperatures by more than 10 C. We have examined the TIMS data from twelve lines flown over Castaic Lake. Four of the lines were flown at an altitude of approximately 2.5 km (MSL), four at an altitude of approximately 6.7 km, and four at approximately 8.3 km. At each altitude there were flights with northern and southern headings, with the aircraft level and at a positive pitch (nose-up attitude). The suite of twelve flights was designed to subject the TIMS/air fence system to different wind conditions and air temperatures. The TIMS flights were supported by a ground-truth team, who measured lake surface temperatures from a boat, and an atmosphere characterization team, who launched an airsonde and measured solar irradiance with a Reagan Sun Photometer. The Reagan measurements were used to construct a time-series of estimates of the total abundance of water vapor in the atmospheric column. These estimates were used to constrain modifications of the airsonde water vapor profile measurements made when processing the TIMS data with a customized version of the MODTRAN radiative transfer code.

Description: The sedimentary geology of the western Grand Canyon consists of gently northeast dipping sandstones, shales, and carbonates. However, due to facies changes within the units, the geomorphology varies from that seen by visitors at the National Park Headquarters. There, the cliff and slope expression of the rocks is replaced in the west by a series of mesas, ridges, and horizontal platforms. The largest of these occurs on the Esplanade Sandstone within the Supai Formation. The Esplanade is formed by slope retreat of the overlying units and resistance to erosion by the underlying limestones. It is onto this platform that the lavas of the Uinkaret Plateau were emplaced. The Uinkaret lava field lies 120 km south of St. George, Utah and is tectonically defined by two major normal faults -- the Hurricane to the west and the Toroweap to the east. The purpose of this investigation was to collect visible, near and thermal infrared data at different periods of the day and year. It is expected that these data will provide the ability to retrieve water temperatures; monitor sediment loads; map and examine any changes in the near shore vegetation communities and understand some of the intricacies of the geology. This paper will serve, to some degree, as a progress report on the Grand Canyon study, since only a fraction of the data has been received and processed thus far. Data from the Thermal Infrared Multispectral Scanner (TIMS) and the Landsat Thematic Mapper simulator (NS001) were acquired simultaneously on April 4, 1994. A second data acquisition occurred on August 27, 1994. Initial analysis of the TIMS data indicates a remarkably noise-free data set with minimal atmospheric attenuation. environments is evident.

Description: The Jet Propulsion Laboratory (JPL) Field Emission Spectrometer (FES) was built by Designs and Prototypes based on a set of functional requirements supplied by JPL. The instrument has a spectral resolution of 6 wavenumbers (wn) and can acquire spectra from either the Mid Infrared (3-5 mu m) or the Thermal Infrared (8-12 pm) depending on whether the InSb or HgCdTe detector is installed respectively. The instrument consists of an optical head system unit and battery. The optical head which is tripod mounted includes the interferometer and detector dewar assembly. Wavelength calibration of the interferometer is achieved using a Helium-Neon laser diode. The dewar needs replenishing with liquid Nitrogen approximately every four hours. The system unit includes the controls for operation and the computer used for acquiring viewing and processing spectra. Radiometric calibration is achieved with an external temperature-controlled blackbody that mounts on the fore-optics of the instrument. The blackbody can be set at 5 C increments between 10 and 55 C. The instrument is compact and weighs about 33 kg. Both the wavelength calibration and radiometric calibration of the instrument have been evaluated. The wavelength calibration was checked by comparison of the position of water features in a spectrum of the sky with their position in the output from a high resolution atmospheric model. The results indicatethat the features in the sky spectrum are within 6-8 wn of their position ill the model spectrum. The radiometric calibration was checked by first calibrating the instrument using the external blackbody supplied with the instrument and then measuring the radiance from another external blackbody at a series of temperatures. The temperatures of these radiance spectra were then recovered by inventing Planck's law and the recovered temperatures compared lo the measured blackbody temperature. These results indicate that radiometric calibration is good to 0.5 C over the range of temperatures 10 to 55 C. The results also indicate that the instrument drifts slowly over time and should be recalibrated every 20 to 30 minutes in the field to ensure good radiometric fidelity. The instrument has now been extensively tested in the field in the United States and Australia. These in situ field measurements are being used to validate emissivity spectra recovered from the Thermal Infrared Multispectral Scanner (TIMS) and also the Australian CO2 Laser. The availability of in situ measurements is proving crucial to validation of the spectra derived from the airborne instruments since many natural surfaces cannot be easily transported back to the laboratory.

Description: AVIRIS is a NASA-sponsored Earth-remote-sensing imaging spectrometer designed, built and operated by the Jet Propulsion Laboratory (JPL). While AVIRIS has been operational since 1989, major improvements have been completed in most of the sensor subsystems during the winter maintenance cycles. As a consequence of these efforts, the capabilities of AVIRIS to reliably acquire and deliver consistently high quality, calibrated imaging spectrometer data continue to improve annually, significantly over those in 1989. Improvements to AVIRIS prior to 1994 have been described previously. This paper details recent and planned improvements to AVIRIS in the sensor task.

Description: This paper serves as a brief overview of the AVIRIS instrument (Airborne Visible/Infrared Imaging Spectrometer). The AVIRIS sensor collects data that will be used for quantitative characterization of the Earth's surface and atmosphere from geometrically coherent spectroradiometric measurements. This data can be applied to studies in the fields of oceanography, environmental science, snow hydrology, geology, volcanology, soil and land management, atmospheric and aerosol studies, agriculture, and limnology. Applications under development include the assessment and monitoring of environmental hazards such as toxic waste, oil spills, and land/air/water pollution. Mission planning and flight operations are discussed, and recommendations are given regarding the deployment of ground truth experiments.

Description: The AVIRIS sensor must be calibrated at the time it measures spectra from the ER-2 airborne platform in order to achieve research and application objectives that are both quantitative and physically based. However, the operational environment inside the Q-bay of the ER-2 at 20 km altitude differs from that in the AVIRIS laboratory with respect to temperature, pressure, vibration, and high-frequency electromagnetic fields. Experiments at surface calibration targets are used in each flight season to confirm the accuracy of AVIRIS in-flight radiometric calibrations. For these experiments, the MODTRAN radiative transfer code is constrained by using in situ measurements to independently predict the upwelling spectral radiance arriving at AVIRIS for a specific calibration target. AVIRIS calibration is validated in flight by comparing the MODTRAN-predicted radiance to the laboratory-calibrated radiance measured by the AVIRIS sensor for the same time over the calibration target. We present radiometric calibration results for the AVIRIS in-flight calibration experiment held at the beginning of the 1994 flight season.

Description: The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) measures spectra from 400 to 2500 nm through 224 contiguous spectral channels. These spectra are used to identify and determine the abundances of constituents of the earth's surface and atmosphere through the analysis of spectral molecular absorptions and scattering properties of materials. Spectral calibration of the AVIRIS is required to pursue these determinations, which are based on these fundamental spectral characteristics. The spectral calibration of AVIRIS consists of a description of the spectral response function for each of the 224 spectral channels. The spectral response function of AVIRIS channels was measured in the laboratory and shown to be Gaussian in shape. A subset of nine AVIRIS spectral channels is shown with channel centers at 10 nm intervals and channel shapes of Gaussian form with 10 nm full width and half maximum (FWHM). For AVIRIS the spectral calibration is described by the spectral center and FWHM of a Gaussian function for each channel. Uncertainty in the determination of the spectral center and FWHM for each channel are also reported. This paper presents results from an investigation of the sensitivity of AVIRIS-measured upwelling radiance to errors in spectral calibration. Based on this sensitivity analysis in conjunction with earlier work, an improved spectral calibration requirement for AVIRIS is proposed.

Description: Movement of water as vapor in the atmosphere is a fundamental process in the Earth's hydrological cycle. Investigations of spatial and time scales of water vapor transport in the atmosphere are important areas of research. Water vapor transmits energy as a function of its abundance across the spectrum. This is shown in the 400- to 2500-nm spectral region where the transmission of the terrestrial atmosphere has been modeled using the MODTRAN radiative transfer code for a range of water vapor abundances. Based on these model results, spectra measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) have been used to investigate the movement of water vapor at 20-m spatial resolution over an 11-by-30-km area at approximately 15-minute time intervals (1.25 hours total). AVIRIS measures the upwelling spectral radiance from 400 to 2500 nm at 10-nm spectral intervals and collects images of 11-by-up-to-1000 km at 20-m spatial resolution. Data are collected at a rate of 1 km of flight line per 4.5 seconds. A set of five AVIRIS flight lines was acquired in rapid succession over Rogers Dry Lake, CA on May 18, 1993 at 18:59, 19:13, 19:29, 19:47, and 19:59 UTC. Rogers Dry lake is located 2 hours north of Los Angeles, CA at 34.84 degrees north latitude and 117.83 degrees west longitude in the Mojave Desert.

Description: Spectral calibration of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) as data are acquired in flight is essential to quantitative analysis of the measured upwelling spectral radiance. In each spectrum measured by AVIRIS in flight, there are numerous atmospheric gas absorption bands that drive this requirement for accurate spectral calibration. If the surface and atmospheric properties are measured independently, these atmospheric absorption bands may be used to deduce the in-flight spectral calibration of an imaging spectrometer. Both the surface and atmospheric characteristics were measured for a calibration target during an in-flight calibration experiment held at Lunar Lake, Nevada on April 5, 1994. This paper uses upwelling spectral radiance predicted for the calibration target with the MODTRAN radiative transfer code to validate the spectral calibration of AVIRIS in flight.

Description: Since initial contact between Earth Search Sciences, Inc. (ESSI) and the Idaho National Engineering Laboratory (INEL) in February, 1994, at least seven proposals have been submitted in response to a variety of solicitations to commercialize and improve the AVIRIS instrument. These proposals, matching ESSI's unique position with respect to agreements with the National Aeronautics and Space Administration (NASA) and the Jet Propulsion Laboratory (JPL) to utilize, miniaturize, and commercialize the AVIRIS instrument and platform, are combined with the applied engineering of the INEL. Teaming ESSI, NASA/JPL, and INEL with diverse industrial partners has strengthened the respective proposals. These efforts carefully structure the overall project plans to ensure the development, demonstration, and deployment of this concept to the national and international arenas. The objectives of these efforts include: (1) developing a miniaturized commercial, real-time, cost effective version of the AVIRIS instrument; (2) identifying multiple users for AVIRIS; (3) integrating the AVIRIS technology with other technologies; (4) gaining the confidence/acceptance of other government agencies and private industry in AVIRIS; and (5) increasing the technology base of U.S. industry.

Description: Recent laboratory calibrations of the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) include new methods for the characterization of the geometric, spectral, temporal and radiometric properties of the sensor. New techniques are desired in order to: (1) increase measurement accuracy and precision, (2) minimize measurement time and expense, (3) prototype new field and inflight calibration systems, (4) resolve measurement ambiguities, and (5) add new measurement dimensions. One of the common features of these new methods is the use of the full data collection and processing power of the AVIRIS instrument and data facility. This allows the collection of large amounts of calibration data in a short period of time and is well suited to modular data analysis routines.

Description: The AVIRIS instrument uses an onboard calibration system to provide auxiliary calibration data. The system consist of a tungsten halogen cycle lamp imaged onto a fiber bundle through an eight position filter wheel. The fiber bundle illuminates the back side of the foreoptics shutter during a pre-run and post-run calibration sequence. The filter wheel contains two neutral density filters, five spectral filters and one blocked position. This paper reviews the general workings of the onboard calibrator system and discusses recent modifications.

Description: The JPL DSN Microwave Antenna Holography System (MAHST) was applied to the newly constructed DSS-24 34-m beam-waveguide antenna at Goldstone, California. The application of MAHST measurements and corrections at DSS 24 provided the critical RF performance necessary to not only meet the project requirements and goals, but to surpass them. A performance increase of 0.35 dB at X-band (8.45 GHz) and 4.9 dB at Ka-band (32 GHz) was provided by MAHST, resulting in peak efficiencies of 75.25 percent at X-band and 60.6 percent at Ka-band (measured from the Cassegrain focus at f1). The MAHST enabled setting the main reflector panels of DSS 24 to 0.25-mm rms, making DSS 24 the highest precision antenna in the NASA/JPL DSN. The precision of the DSS-24 antenna (diameter/rms) is 1.36 x 10(exp 5), and its gain limit is at 95 GHz.

Description: The Voyager 1 spacecraft flew by Jupiter on March 5, 1979. Spacecraft navigation was performed with radio tracking data from NASA's Deep Space Network. In the years since then, there has been a great deal of progress in the definition of celestial reference frames and in determining the orbit and orientation of the Earth. Using these improvements, the radio metric range and Doppler data acquired from the Voyager 1 spacecraft near its encounter with Jupiter have been reanalyzed to determine the plane-of-sky position of Jupiter with much greater accuracy than was possible at the time of the encounter. The position of Jupiter at the time of encounter has been determined with an accuracy of 40 nrad in right ascension and 140 nrad in declination with respect to the celestial reference frame defined by the International Earth Rotation Service. This position estimate has been done to improve the ephemeris of Jupiter prior to the upcoming encounter of the Galileo spacecraft with Jupiter.

Description: Frequency measurements of electromagnetic oscillations of millimeter and submillimeter wavebands with frequency growth due to a number of reasons become more and more difficult. First, these frequencies are considered to be cutoffs for semiconductor converting devices and one has to use optical measurement methods instead of traditional ones with frequency transfer. Second, resonance measurement methods are characterized by using relatively narrow bands and optical ones are limited in frequency and time resolution due to the limited range and velocity of movement of their mechanical elements as well as the efficiency of these optical techniques decrease with the increase of wavelength due to diffraction losses. That requires a priori information on the radiation frequency band of the source involved. Method of measuring frequency of harmonic microwave signals in millimeter and submillimeter wavebands based on the ac Josephson effect in superconducting contacts is devoid of all the above drawbacks. This approach offers a number of major advantages over the more traditional measurement methods, that is one based on frequency conversion, resonance and interferometric techniques. It can be characterized by high potential accuracy, wide range of frequencies measured, prompt measurement and the opportunity to obtain a panoramic display of the results as well as full automation of the measuring process.

Description: As data communications rates climb toward 10 Gbits/s, clock recovery and synchronization become more difficult, if not impossible, using conventional electronic circuits. The high-speed photonic clock regenerator described in this article may be more suitable for such use. This photonic regenerator is based on a previously reported photonic oscillator capable of fast acquisition and synchronization. With both electrical and optical clock inputs and outputs, the device is easily interfaced with fiber-optic systems. The recovered electrical clock can be used locally and the optical clock can be used anywhere within a several kilometer radius of the clock/carrier regenerator.

Description: Industry and university participants have joined together to form the IMPA:Ct consortium (In-situ Monitors of the Particulate Ambient: Circumterrestrial) which offers a broad range of flight qualified instruments for monitoring the small particle (0.1 micron to 10 cm) environment in space. Instruments are available in 12 months or less at costs ranging from 0.5 to 1.5 million dollars (US) for the total program. Detector technologies represented by these groups are: impact-induced capacitor-discharge (MOS, metal-oxide-silicon), cratering or penetration of electroactive thin film (polyvinylidene fluoride (PVDF)), impact-plasma detection, acoustic detection, CCD tracking of optical scatter of sunlight, and photodiode detection of optical scatter of laser light. The operational characteristics, general spacecraft interface and resource requirements (mass/power/telemetry), cost and delivery schedules, and points of contact for seven different instruments are presented.

Description: The Tropical Rainfall Measuring Mission (TRMM) is an earth observing satellite that will be in a low earth orbit (350 kilometers) during the next period of maximum solar activity. The TRMM observatory is expected to experience an atomic oxygen fluence of 8.9 x 10(exp 22) atoms per square centimeter. This fluence is ten times higher than the atomic oxygen impingement incident to the Long Duration Exposure Facility (LDEF). Other environmental concerns on TRMM include: spacecraft glow, silicon oxide contaminant build-up, severe spacecraft material degradation, and contamination deposition resulting from molecular interactions with the dense ambient atmosphere. Because of TRMM's predicted harsh environment, TRMM faces many unique material concerns and subsystem design issues. The LDEF data has influenced the design of TRMM and the TRMM material selection process.

Description: This paper describes the requirements, design, assembly and testing of the linear Scan Mechanism (SM) of the Composite Infrared Spectrometer (CIRS) Instrument. The mechanism consists of an over constrained flexible structure, an innovative moving magnet actuator, passive eddy current dampers, a Differential Eddy Current (DEC) sensor, Optical Limit Sensors (OLS), and a launch lock. Although all the components of the mechanism are discussed, the flexible structure and the magnetic components are the primary focus. Several problems encountered and solutions implemented during the development of the scan mechanism are also described.

Description: On July 4,1997, the Mars Pathfinder spacecraft lands on Mars and starts conducting technological and scientific experiments. One experiment, the Alpha-Proton-X-ray Spectrometer, uses a sensor head placed against rocks and soil to determine their composition. To guarantee proper placement, a deployment mechanism mounted on the Mars Rover aligns the sensor head to within 20 deg of the rock and soil surfaces. In carrying out its task, the mechanism mimics the action of a human hand and arm. Consisting of a flexible wrist, a parallel link arm, a brush dc motor actuator, and a revolutionary non-pyrotechnic fail-safe release device, the mechanism correctly positions the sensor head on rocks as high as 0.29 m and on targets whose surfaces are tilted as much as 45 deg from the nominal orientation of the sensor head face. The mechanism weighs less than 0.5 kg, can withstand 100 g's, and requires less than 2.8 N x m of actuation torque. The fail-safe coupler utilizes Cerrobend, a metal alloy that melts at 60 C, to fuse the actuator and the rest of the mechanism together. A film heater wrapped around the coupler melts the metal, and Negator springs drive the mechanism into its stowed position. The fail-safe actuates using 6.75 Watts for 5 minutes in the event of an actuator failure.

Description: The University of California at Santa Cruz and the NASA Ames Research Center have developed a mid-infrared Si:Ga 128 x 128 array camera. The camera has been used on several flights on the KAO to obtain data in the mid-infrared wavelength (4-8 micron) range. At the f/17 bent Cassegrain focus, the pixel spacing is 1.6 arcsec, resulting in a FOV of 3.4 arcmin. The bandpass of the camera is determined by a set of 4 fixed filters and a 4-8 micron CVF. Camera performance of 92mJy/square arcsec/minute has been achieved, but is far from optimal. Realistic improvements to the camera and pointing stability on the KAO in the near future should increase the camera's sensitivity on future flights by 5-10x.

Description: We describe a focal plane guider for the Kuiper Airborne Observatory which consists of a CCD camera interfaced to an AMIGA personal computer. The camera is made by Photometrics Ltd. and utilizes a Thomson 576 x 384 pixel CCD chip operated in Frame Transfer mode. Custom optics produce a scale of 2.4 arc-sec/pixel, yielding an approx. 12 ft. diameter field of view. Chopped images of stars with HST Guide Star Catalog magnitudes fainter than 14 have been used for guiding at readout rates greater than or equal to 0.5 Hz. The software includes automatic map generation, subframing and zooming, and correction for field rotation when two stars are in the field of view.

Description: The High-efficiency, Infrared Faint Object Grating Spectrometer (HIFOGS) provides spectral coverage of selectable portions of the 3 to 18 micron range at resolving powers from 00 to 1000 using 120 Si/Bi detectors. Additional coverage to 30 microns is provided by a bank of 32 Si:P detectors. Selectable apertures, gratings and band-pass filters provide flexibility to this system. Software for operation of HIFOGS and reduction of the data runs on a MacIntosh computer. HIFOGS has been used to establish celestial flux standards using 3 independent approaches: comparison to star models, comparisons to asteroid models and comparisons to laboratory blackbodies. These standards are expected to have wide application in astronomical thermal-infrared spectroscopy.