ALBERT

Description: Ions of galactic origin are modified but not attenuated by the presence of shielding materials. Indeed, the number of particles and the absorbed energy behind most shield materials increases as a function of shield thickness. The modification of the galactic cosmic ray composition upon interaction with shielding is the only effective means of providing astronaut protection. This modification is intimately conntected with the shield transport porperties and is a strong function of shield composition. The systematic behavior of the shield properites in terms of microscopic energy absorption events will be discussed. The shield effectiveness is examined with respect to convectional protection practice and in terms of a biological endpoint: the efficiency for reduction of the probability of transformation of shielded C3H1OT1/2 mouse cells. The relative advantage of developing new shielding technologies is discussed in terms of a shield performance as related to biological effect and the resulting uncertainty in estimating astronaut risk.

Description: Variations in the Earth's trapped (Van Allen) belts produced by solar flare particle events are not well understood. Few observations of increases in particle populations have been reported. This is particularly true for effects in low Earth orbit, where manned spaceflights are conducted. This paper reports the existence of a second proton belt and it's subsequent decay as measured by a tissue-equivalent proportional counter and a particle spectrometer on five Space Shuttle flights covering an eighteen-month period. The creation of this second belt is attributed to the injection of particles from a solar particle event which occurred at 2246 UT, March 22, 1991. Comparisons with observations onboard the Russian Mir space station and other unmanned satellites are made. Shuttle measurements and data from other spacecraft are used to determine that the e-folding time of the peak of the second proton belt. It was ten months. Proton populations in the second belt returned to values of quiescent times within eighteen months. The increase in absorbed dose attributed to protons in the second belt was approximately 20%. Passive dosimeter measurements were in good agreement with this value.

Description: It has been known for some time that adequate assessment of spacecraft requirements and concomitant estimates of astronaut radiation exposures from galactic cosmic radiation requires accurate, quantitative methods for characterizing these radiation fields as they pass through thick absorbers. The main nuclear interaction processes involved are nuclear elastic an inelastic collisions, and nuclear breakup (fragmentation) and electromagnetic dissociation (EMD). Nuclear fragmenation and EMD are important because they alter the elemental and isotopic composition of the transported radiation fields. At present, there is no suitable accurate theory for predicting nuclear fragmentation cross sections for all collision pairs and energies of interest in space radiation protection. Typical cross-section differences between theory and experiment range from about 25 percent to a factor of two. The resulting errors in transported flux, for high linear energy transfer (LET) particles, are compared to these cross-sections errors. In this overview, theoretical models of heavy ion fragmentation currently used to generate input data bases for cosmic-ray transport and shielding codes are reviewed. Their shortcomings are discussed. Further actions needed to improve their accuracy and generality are presented.

Description: Both man and technological equipment must survive the near-earth space radiation environment, which can, under specific conditions, be extremely severe. This conference produced 17 papers on the dynamic space radiation environment covering: galactic, solar and trapped particles; nuclear fragmentation; nuclear interactions and transport theory; solar proton events; radiation shielding; and heavy ion fluences. Several papers present results from the recent SAMPEX mission.

Description: The fragmenting of high energy, heavy ions (HZE particles) by hydrogen targets is an important, physical process in several areas of space radiation research. In this work quantum mechanical optical model methods for estimating cross sections for HZE particle fragmentation by hydrogen targets are presented. The cross sections are calculated using a modified abrasion-ablation collision formalism adapted from a nucleus-nucleus collision model. Elemental and isotopic production cross sections are estimated and compared with reported measurements for the breakup of neon, sulphur, and iron, nuclei at incident energies between 400 and 910 Mev/nucleon. Good agreement between theory and experiment is obtained.

Description: Light nuclei are present in the primary galactic cosmic rays (GCR) and are produced in thick targets due to projectile or target fragmentation from both nucleon and heavy induced reactions. In the primary GCR, He-4 is the most abundant nucleus after H-1. However, there are also a substantial fluxes of H-2 and He-3. In this paper we describe theoretical models based on quantum multiple scattering theory for the description of light ion nuclear interactions. The energy dependence of the light ion fragmentation cross section is considered with comparisons of inclusive yields and secondary momentum distributions to experiments described. We also analyze the importance of a fast component of lights ions from proton and neutron induced target fragementation. These theoretical models have been incorporated into the cosmic ray transport code HZETRN and will be used to analyze the role of shielding materials in modulating the production and the energy spectrum of light ions.

Description: A model for the differential energy spectra of galactic cosmic radiation as a function of solar activity is described. It is based on the standard diffusion-convection theory of solar modulation. Estimates of the modulation potential based on fitting this theory to observed spectral measurements from 1954 to 1989 are correlated to the Climax neutron counting rates and to the sunspot numbers at earlier times taking into account the polarity of the interplanetary magnetic field at the time of observations. These regression lines then provide a method for predicting the modulation at later times. The results of this model are quantitatively compared to a similar Moscow State University (MSU) model. These model cosmic ray spectra are used to predict the linear energy transfer spectra, differential energy spectra of light (charge less than or = 2) ions, and single event upsets rates in memeory devices. These calculations are compared to observations made aboard the Space Shuttle.

Description: The large solar energetic particle (SEP) events and simultaneous large geomagnetic disturbances observed during October 1989 posed a significant, rapidly evolving space radiation hazard. Using data from the GOES-7, NOAA-10, IMP-8 and LDEF satellites, we determined the geomagnetic transmission, heavy ion fluences, mean Fe ionic charge state, and effective radiation hazard observed in low Earth orbit (LEO) for these SEPs. We modeled the geomagneitc transmission by tracing particles through the combination to the internal International Geomagnetic Reference Field (IGRF) and the Tsyganenko (1989) magnetospheric field models, extending the modeling to large geomagnetic disturbances. We used our results to assess the radiation hazard such very large SEP events would pose in the anticipated 52 deg inclination space station orbit.

Description: Observations aboard Cosmos satelites discovered trapped anomalous cosmic rays (ACRs), tracked the variation in their intensity in 1986-1988, and measured their fluence, spectrum, and composition at solar minimum in the previous solar cycle. The MAST instrument aboard the SAMPEX satellite has observed trapped anomalous cosmic rays in the present solar cycle, confirmed the general features of the Cosmos data, and provided the first detailed observations of trapped ACRs. In this paper we apply theoretical modeling of trapped ACRs, which is shown to provide a reasonably good description of both the Cosmos and SAMPEX data, to calculate the integral linear-energy-transfer (LET) spectra due to trapped ACRs in typical low-Earth orbits. We compare these calculations with the LET spectra produced by galactic cosmic rays (GCRs) and non-trapped ACRs in order to assess the relative radiation hazard posed by trapped ACRs.

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: 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: A number of planetary objects exhibit large radar reflectivity and polarization ratios, and more recently, a similar behavior has been observed over a vast portion of the Earth's surface: the percolation facies of the Greenland Ice Sheet. Surface-based ranging radar data and snow stratigraphy studies demonstrated that the radar properties of that portion of Greenland are caused by enhanced scattering from massive, large, solid-ice bodies buried in the top few meters of the dry, cold, clean snowy surface of the ice sheet and created by seasonal melting and refreezing events. Here, we model the icy inclusions as randomly oriented, discrete, noninteracting, dielectric cylinders embedded in a transparent snow medium. An exact analytical solution is used to compute the scattered field from the cylinders. Using this model, we correctly predict the polarimetric radar observations gathered by an airborne imaging system at three wavelengths (5.6, 24, and 68 cm), between 19 deg and 65 deg incidence angle. The diameter and number density of the cylinders that are inferred from the radar data using the model are consistent with in situ observations of the icy inclusions. The large radar reflectivity and polarization ratios are interpreted as arising internal reflections of the radar signals in the icy inclusions that first-order external reflection models fail to predict. The results compare favorably with predictions from the coherent backscatter or weak localization theory and may provide a complementary framework for interpreting exotic radar echoes from other planetary objects.

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 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 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: Synthetic-aperture radar (SAR) interferograms produced from ESA's ERS-1 satellite, provide the first synoptic view of ice flow dynamics of the western sector of the Greenland Ice Sheet. Glacial motion is detected in the radar ranging direction at millimetric scales, across a complete sequence of snow accumulation and melting regimes, despite significant varations in their radar scattering properties. Ice flow evolves from a slow, regular motion at the higher elevations. At lower elevations, motion is strongly convoluted by meter-scale undulations in surface topography, which have a unique interferometric signature that enables a novel approach for retrieving flow direction. Inferred flow directions, combined with surface displacements in the radar ranging direction, yield ice velocity estimates that are within 6% of in-situ measurements gathered along a 40 km survey line. Application of repeat-pass SAR interferometry to the entire Greenland Ice Sheet should enable precise mapping of its ice flow dynamics at an unprecedented level of spatial detail.

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 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: 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 present interferometric observations of ice-sheet motion in western Greenland based on pairs of ERS-1 synthetic aperture radar (SAR) images. These observations provide the first detailed regional view of ice motion associated with dynamically supported topography near the margin of an ice sheet. The interferograms of this area are much more complicated than other interferograms of ice sheets presented to date. We devote the largest part of this paper to explaining the source of the complexity in these interferograms. A synthetic interferogram based on a simple model helps to illustrate the effects of different components of the ice velocity field in interferometric data and suggests a method for estimating the large-scale ice velocity field from such data.

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: 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: 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 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: 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: 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: The new information on galactic cosmic rays (GCR) derived from the Spacelab-3 cosmic ray experiment 'Anuradha' shows that at 25-125 MeV/N GCR sub-iron and iron (Z = 21-28) particles consists of a mixture of partially ionized and fully ionized ions. Computation of electron capture and loss cross sections in hydrogen in 1-50 MeV/N energy range are made for Fe, Cr, Ti and Ni. From these it is concluded that: (1) these GCR particles must have captured orbital electrons at energies of about 1-5 MeV/N and (2) these particles are then reaccelerated to 300-500 MeV/N most probably in interstellar medium by collision with SNR shock fronts. Some reacceleration may take place also in heliospheric boundary region. It is suggested that these observations of partially ionized GCR ions of about 100 MeV/N in Spacelab-3 provide a direct evidence of reacceleration of GCR.

Description: Get-Away Special (GAS) G-301, named the Flying Falcon and scheduled for launch on the STS-77 Space Shuttle in April, 1996, is being prepared to perform an experiment designed by the Department of Geology, Bowling Green State University (BGSU). The experiment will employ a new type of infrared imager designed and built by a consortium of Teltron Technologies Inc., Hudson Research Inc., and BGSU that is an uncooled, quantrum ferro-electric, infrared return beam vidicon (IRBV) camera capable of detecting thermal infrared radiation throughout the 2.0-50.0 micron wavelength region, and to which an integral, unable Fabry-Perot filter and a telescopic lens have been added. The primary objectives in the experiment include the mapping of methane plumes from solid waste landfills and wetlands in the midwestern U.S., the mapping of methane plumes offshore in the Gulf of Mexico and in the Middle East, brief monitoring for precursors of volcanoes or earthquakes in the South China sea and the East Pacific Rise (about 300 km west of Easter Island), and the mapping of silica content in exposed outcrops and residual soils of the southwestern U.S. and Middle East.

Description: The GAMCIT (Gamma-ray Astrophysics Mission, California Institute of Technology) payload is a Get-Away-Special payload designed to search for high-energy gamma-ray bursts and any associated optical transients. This paper presents details on the development and construction of the GAMCIT payload. In addition, this paper will reflect upon the unique challenges involved in bringing the payload close to completion, as the project has been designed, constructed, and managed entirely by undergraduate members of the Caltech SEDS (Students for the Exploration and Development of Space). Our experience will definitely be valuable to other student groups interested in undertaking a challenge such as a Get-Away-Special payload.

Description: Monitoring spatial and temporal changes of soil moisture are of importance to hydrology, meteorology, and agriculture. This paper reports a result on study of using L-band SAR imagery to estimate soil moisture and surface roughness for bare fields. Due to limitations of the Small Perturbation Model, it is difficult to apply this model on estimation of soil moisture and surface roughness directly. In this study, we show a simplified model derived from the Integral Equation Model for estimation of soil moisture and surface roughness. We show a test of this model using JPL L-band AIRSAR data.

Description: In the tropical rain forests of Manu, in Peru, where forest biomass ranges from 4 kg/sq m in young forest succession up to 100 kg/sq m in old, undisturbed floodplain stands, the P-band polarimetric radar data gathered in June of 1993 by the AIRSAR (Airborne Synthetic Aperture Radar) instrument separate most major vegetation formations and also perform better than expected in estimating woody biomass. The worldwide need for large scale, updated biomass estimates, achieved with a uniformly applied method, as well as reliable maps of land cover, justifies a more in-depth exploration of long wavelength imaging radar applications for tropical forests inventories.

Description: Radar backscatter intensity as measured by calibrated synthetic aperture radar (SAR) systems is primarily controlled by three factors: local incidence angle, wavelength-scale roughness, and dielectric permittivity of surface materials. In order to make adequate use of radar observations for geological investigations of surface type, the relationships between lithology and the above characteristics must be adequately understood. In arid terrains weathering signatures (e.g. fracturing, debris grain size and shape, slope characteristics) are controlled to some extent by lithologic characteristics of the parent bedrock. These textural features of outcrops and their associated debris control radar backscatter to varying degrees. The quad-polarization JPL AIRSAR system allows sampling of textures at three distinct wavelength scales: C-band (5.66 cm), L-band (23.98 cm), and P-band (68.13 cm). This paper presents a discussion of AIRSAR data using recent field observations of weathered felsic and basaltic volcanic rock units exposed in the southern part of the Lunar Crater Volcanic Field, in the Pancake Range of central Nevada. The focus is on the relationship of radar backscatter at multiple wavelengths to weathering style and parent bedrock lithology.

Description: Silicic (silica-rich) lava flows, such as rhyolite, rhyodacite, and dacite, possess unique physical properties primarily because of the relatively high viscosity of the molten lava. Silicic flows tend to be thicker than basaltic flows, and the resulting large-scale morphology is typically a steep-sided dome or flow lobe, with aspect ratios (height/length) sometimes approaching unity. The upper surfaces of silicic domes and flows are normally emplaced as relatively cool, brittle slabs that fracture as they are extruded from the central vent areas, and are then rafted away toward the flow margin as a brittle carapace above a more ductile interior layer. This mode of emplacement results in a surface with unique roughness characteristics, which can be well-characterized by multiparameter synthetic aperture radar (SAR) observations. In this paper, we examine the scattering properties of several silicic domes in the Inyo volcanic chain in the Eastern Sierra of California, using AIRSAR and TOPSAR data. Field measurements of intermediate-scale (cm to tens of m) surface topography and block size are used to assess the mechanisms of the scattering process, and to quantify the unique roughness characteristics of the flow surfaces.

Description: The replacement of semidesert grassland by woody shrubland is a widespread form of desertification. This change in physiognomy and species composition tends to sharply reduce the productivity of the land for grazing by domestic livestock, increase soil erosion and reduce soil fertility, and greatly alter many other aspects of ecosystem structure and functioning. Remote sensing methods are needed to assess and monitor shrubland encroachment. Detection of woody shrubs at low density would provide a particularly useful baseline on which to access changes, because an initially low shrub density often tends to increase even after cessation of the disturbance (e.g., overgrazing, drought, or fire suppression) responsible for triggering the initial stages of the invasion (Grover and Musick, 1990). Limited success has been achieved using optical remote sensing. In contrast to other forms of desertification, biomass does not consistently decrease with a shift from grassland to shrubland. Estimation of green vegetation amount (e.g., by NDVI) is thus of limited utility, unless the shrubs and herbaceous plants differ consistently in phenology and the area can be viewed during a season when only one of these is green. The objective of this study was to determine if the potential sensitivity of active microwave remote sensing to vegetation structure could be used to assess the degree of shrub invasion of grassland. Polarimetric Airborne Synthetic Aperture Radar (AIRSAR) data were acquired for a semiarid site containing varied mixtures of shrubs and herbaceous vegetation and compared with ground observations of vegetation type and other landsurface characteristics. In this preliminary report we examine the response of radar backscatter intensity to shrub density. The response of other multipolarization parameters will be examined in future work.

Description: Conventional representations of polarization response are referred to a horizontally and vertically polarized basis. Recent studies by Freeman and Durden, van Zyl, and others suggest that alternative polarimetric features which more easily resolve the contributions of simple scattering mechanisms such as odd-bounce, even-bounce, and diffuse scattering could offer several advantages in terrain classification. The circular polarization covariance matrix is a potential source of such features. In this paper, we derive its relationship to the Stokes matrix, describe some of its properties, and compare the utility of linear and circular polarimetric features in classifying an AIRSAR scene containing urban, park, and ocean terrain.

Description: The relationship between habitat area, spatial dynamics of the landscape, and species diversity is an important theme in population and conservation biology. Of particular interest is how populations of various species are affected by increasing habitat edges due to fragmentation. Over the last decade, assumptions regarding the effects of habitat edges on biodiversity have fluctuated wildly, from the belief that they have a positive effect to the belief that they have a clearly negative effect. This change in viewpoint has been brought about by an increasing recognition of the importance of geographic scale and a reinterpretation of natural history observations. In this preliminary report from an ongoing project, we explore the use of remote sensing technology and geographic information systems to further our understanding of how species diversity and population density are affected by habitat heterogeneity and landscape composition. A primary feature of this study is the investigation of SAR for making more rigorous investigations of habitat structure by exploiting the interaction between radar backscatter and vegetation structure and biomass. A major emphasis will be on the use of SAR data to define relative structural types based on measures of structural consolidation using the vegetation surface area to volume ratio (SA/V). Past research has shown that SAR may be sensitive to this form of structural expression which may affect biodiversity.

Description: Landforms in arid regions record the interplay between tectonic forces and climate. Alluvial fans are a common landform in desert regions where the rate of uplift is greater than weathering or sedimentation. Changes in uplift rate or climatic conditions can lead to isolation of the currently forming fan surface through entrenchment and construction of another fan either further from the mountain front (decreased uplift or increased runoff) or closer to the mountain front (increased uplift or decreased runoff). Thus, many alluvial fans are made up of a mosaic of fan units of different age, some older than 1 million years. For this reason, determination of the stages of fan evolution can lead to a history of uplift and runoff. In an attempt to separate the effects of tectonic (uplift) and climatic (weathering, runoff, sedimentation) processes on the shapes of alluvial fan units, a modified conic equation developed by Troeh (1965) was fitted to TOPSAR digital topographic data for the Trail Canyon alluvial fan in Death Valley, California. This allows parameters for the apex position, slope, and radial curvature to be compared with unit age.

Description: During recent years signature analysis, classification, and modeling of Synthetic Aperture Radar (SAR) data as well as estimation of geophysical parameters from SAR data have received a great deal of interest. An important requirement for the quantitative use of SAR data is the accurate estimation of the backscattering coefficient sigma(exp 0). In terrain with relief variations radar signals are distorted due to the projection of the scene topography into the slant range-Doppler plane. The effect of these variations is to change the physical size of the scattering area, leading to errors in the radar backscatter values and incidence angle. For this reason the local incidence angle, derived from sensor position and Digital Elevation Model (DEM) data must always be considered. Especially in the airborne case, the antenna gain pattern can be an additional source of radiometric error, because the radar look angle is not known precisely as a result of the the aircraft motions and the local surface topography. Consequently, radiometric distortions due to the antenna gain pattern must also be corrected for each resolution cell, by taking into account aircraft displacements (position and attitude) and position of the backscatter element, defined by the DEM data. In this paper, a method to derive an accurate estimation of the backscattering coefficient using NASA/JPL AIRSAR data is presented. The results are evaluated in terms of geometric accuracy, radiometric variations of sigma(exp 0), and precision of the estimated forest biomass.

Description: An empirical model was developed to infer soil moisture and surface roughness from radar data. The accuracy of the inversion technique is assessed by comparing soil moisture obtained with the inversion technique to in situ measurements. The effect of vegetation on the inversion is studied and a method to eliminate the areas where vegetation impairs the algorithm is described.

Description: In 1993 NASA's Jet Propulsion Laboratory Airborne Synthetic Aperture Radar system (AIRSAR) was deployed to South America to collect multi-parameter radar data over pre-selected targets. Among the sites targeted was a series of wind streaks located in the Altiplano of Bolivia. The objective of this investigation is to study the effect of wavelength, polarization, and incidence angle on the visibility of wind streaks in radar data. Because this is a preliminary evaluation of the recently acquired data we will focus on one scene and, thus, only on the effects of wavelength and polarization. Wind streaks provide information on the near-surface prevailing winds and on the abundance of winderodible material, such as sand. The potential for a free-flyer radar system that could provide global radar images in multiple wavelengths, polarizations, and incidence angles requires definition of system parameters for mission planning. Furthermore, thousands of wind streaks were mapped from Magellan radar images of Venus; their interpretation requires an understanding of the interaction of radar with wind streaks and the surrounding terrain. Our experiment was conducted on wind streaks in the Altiplano of Bolivia to address these issues.

Description: The Thermal Infrared Multispectral Scanner (TIMS) was flown on the NASA C-130 aircraft for a series of 12 flights during HAPEX-Sahel at altitudes ranging from 0.25 to 6 km (0.6 to 15 m resolution). TIMS provides coverage of the 8 to 12 micrometer thermal infrared band in 6 contiguous channels. Thus it is possible to observe the spectral behavior of the surface emissivity over this wavelength interval.

Description: Thermal Infrared Multispectral Scanner (TIMS) data were acquired over the McDowell Mountains northeast of Scottsdale, Arizona during August 1994. The raw data were processed to emphasize lithologic differences using a decorrelation stretch and assigning bands 5, 3, and 1 to red, green, and blue, respectively. Processed data of alluvium flanking the mountains exhibit moderate color variation. The objective of this study was to determine, using a quantitative approach, what environmental variable(s), in the absence of bedrock, is/are responsible for influencing the spectral properties of the desert alluvial surface.

Description: In mineral exploration, the ability to distinguish and map petrochemical variations of magmatic rocks can be a useful reconnaissance tool. Alkalinity is one such petrochemical parameter and is used in the characterization of granitoid rocks. In quartz normative plutonic rocks, alkalinity is related to the composition and abundance of feldspars. Together with quartz abundance, knowledge of feldspar modes allows the classification of these igneous rocks according to the Streckeisen diagram. Alternative classification schemes rely on whole rock geochemistry instead of mineral identifications. The relative ease of obtaining whole rock analyses means that geochemical classifications tend to be favored in exploration geology. But the technique of thermal infrared spectroscopy of rocks yields information on mineralogy and is one that can be applied remotely. The goal of the current work then is to establish whether data from TIMS can be used to distinguish the mineralogical variations that relate to alkalinity. An ideal opportunity to test this thesis arises from the work presented in a paper by Dewitt (1989). This paper contains the results of mapping and analysis of Proterozoic plutonic rocks in north-central Arizona. The map resulting from this work delineates plutons according to alkalinity in an effort to establish a trend or polarity in the regional magmatism. Also contained within this paper are brief descriptions of the mineralogy of half of the region's plutons. This combination of mineralogical and geochemical information was the rationale behind choosing this area as a site for TIMS over flights. A portion of the region centered on the northern Bradshaw Mountains was selected because it contains plutons of all three alkalinity classifications (alkali-calcic, calc-alkalic, and calic) present on DeWitt's map within a relatively small area. The site was flown in August of 1994 and the data received a few days before the writing of this manuscript. Most of this paper is devoted to the description of laboratory based spectroscopy and spectral simulations. These are required to gain insight into the correct procedures for enhancing the relatively small differences in the low spectral resolution TIMS data.

Description: This paper describes an empirical method to correct TIMS (Thermal Infrared Multispectral Scanner) data for atmospheric effects by transferring calibration from a laboratory thermal emission spectrometer to the TIMS multispectral image. The method does so by comparing the laboratory spectra of samples gathered in the field with TIMS 6-point spectra for pixels at the location of field sampling sites. The transference of calibration also makes it possible to use spectra from the laboratory as endmembers in unmixing studies of TIMS data.

Description: Remote sensing is a tool that, in the context of aeolian studies, offers a synoptic view of a dune field, sand sea, or entire desert region. Blount et al. (1990) presented one of the first studies demonstrating the power of multispectral images for interpreting the dynamic history of an aeolian sand sea. Blount's work on the Gran Desierto of Mexico used a Landsat TM scene and a linear spectral mixing model to show where different sand populations occur and along what paths these sands may have traveled before becoming incorporated into dunes. Interpretation of sand transport paths and sources in the Gran Desierto led to an improved understanding of the origin and Holocene history of the dunes. With the anticipated advent of the EOS-A platform and ASTER thermal infrared capability in 1998, it will become possible to look at continental sand seas and map sand transport paths using 8-12 mu m bands that are well-suited to tracking silicate sediments. A logical extension of Blount's work is to attempt a similar study using thermal infrared images. One such study has already begun by looking at feldspar, quartz, magnetite, and clay distributions in the Kelso Dunes of southern California. This paper describes the geology and application of TIMS image analysis of a less-well known Holocene dune field in south central Oregon using TIMS data obtained in 1991.

Description: A research program has been initiated between Arizona State University and the City of Scottsdale, Arizona to study the potential applications of TIMS (Thermal Infrared Multispectral Scanner) data for urban scene classification, desert environmental assessment, and change detection. This program is part of a long-term effort to integrate remote sensing observations into state and local planning activities to improve decision making and future planning. Specific test sites include a section of the downtown Scottsdale region that has been mapped in very high detail as part of a pilot program to develop an extensive GIS database. This area thus provides excellent time history of the evolution of the city infrastructure, such as the timing and composition of street repavement. A second area of study includes the McDowell intensive study by state and local agencies to assess potential sites for urban development as well as preservation. These activities are of particular relevance as the Phoenix metropolitan area undergoes major expansion into the surrounding desert areas. The objectives of this study in urban areas are aimed at determining potential applications of TIMS data for classifying and assessing land use and surface temperatures. Land use centers on surface impermeability studies for storm runoff assessment and pollution control. These studies focus on determining the areal abundance of urban vegetation and undeveloped soil. Highly experimental applications include assessment and monitoring of pavement condition. Temperature studies focus on determining swimming pool area and temperature for use in monitoring evaporating and urban water consumption. These activities are of particular relevance as the Phoenix metropolitan area undergoes major expansion into the surrounding desert area.

Description: Remote sensing of sand dunes helps in the understanding of aeolian process and provides important information about the regional geologic history, environmental change, and desertification. Remotely sensed data combined with field studies are valuable in studying dune morphology, regional aeolian dynamics, and aeolian depositional history. In particular, active and inactive sands of the Kelso Dunes have been studied using landsat TM and AIRSAR. In this report, we describe the use of AVIRIS data to study the Kelso dunes and to compare the AVIRIS information with that from TM and AIRSAR.

Description: Vulcano Island is part of the Eolian archipelago, located about 25 km from the northeast coast of Sicily. The archipelago comprises seven major volcanic islands, two of which are active volcanoes (Vulcano and Stromboli). Vulcano covers an area of about 50 square km, and is about 10 km long. Explosive volcanic activity has predominated in the geological evolution of Vulcano Island, and there is no evidence that this pattern has ceased. Rather, the current situation is one of unrest, so a strict regimen of continuous geophysical and geochemical monitoring has been undertaken over the last decade. Though the year-round population of Vulcano is small (under 1000), during the summer the island becomes a very popular resort, and has thousands of additional tourists at any time throughout the high season, thus substantially increasing the number of people potentially at risk from an explosive eruption or other hazards such as noxious gas emissions (e.g., CO2, H2S, SO2). During the past ten years, remote sensing data have been repetitively acquired with optical and microwave airborne sensors. The present work shows the preliminary results of a study based on the integration of various remote sensing data sets with field spectroscopy, and other laboratory analyses, for the geological and geomorphological mapping of the island. It is hoped that such work will also usefully contribute to the evaluation of the volcanic hazard potential of the islands as well as to the evaluation of the status of its current activity.

Description: The rapid development of sophisticated imaging spectrometers and resulting flood of imaging spectrometry data has prompted a rapid parallel development of spectral-information extraction technology. Even though these extraction techniques have evolved along different lines (band-shape fitting, endmember unmixing, near-infrared analysis, neural-network fitting, and expert systems to name a few), all are limited by the spectrometer's signal to noise (S/N) and spectral resolution in producing useful information. This study grew from a need to quantitatively determine what effects these parameters have on our ability to differentiate between mineral absorption features using a band-shape fitting algorithm. We chose to evaluate the AVIRIS, HYDICE, MIVIS, GERIS, VIMS, NIMS, and ASTER instruments because they collect data over wide S/N and spectral-resolution ranges. The study evaluates the performance of the Tricorder algorithm, in differentiating between mineral spectra in the 0.4-2.5 micrometer spectral region. The strength of the Tricorder algorithm is in its ability to produce an easily understood comparison of band shape that can concentrate on small relevant portions of the spectra, giving it an advantage over most unmixing schemes, and in that it need not spend large amounts of time reoptimizing each time a new mineral component is added to its reference library, as is the case with neural-network schemes. We believe the flexibility of the Tricorder algorithm is unparalleled among spectral-extraction techniques and that the results from this study, although dealing with minerals, will have direct applications to spectral identification in other disciplines.

Description: Surface reflectance retrieval from imaging spectrometer data as acquired with the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) has become important for quantitative analysis. In order to calculate surface reflectance from remotely measured radiance, radiative transfer codes such as 5S and MODTRAN2 play an increasing role for removal of scattering and absorption effects of the atmosphere. Accurate knowledge of the exo-atmospheric solar irradiance (E(sub 0)) spectrum at the spectral resolution of the sensor is important for this purpose. The present study investigates the impact of differences in the solar irradiance function, as implemented in a modified version of 5S (M5S), 6S, and MODTRAN2, and as proposed by Green and Gao, on the surface reflectance retrieved from AVIRIS data. Reflectance measured in situ is used as a basis of comparison.

Description: The potential of airborne imaging spectrometry for assessing and monitoring natural resources is studied. Therefore, an AVIRIS scene of the NASA's MacEurope 1991 campaign - acquired in Central Switzerland - is used. The test site consists of an urban area, the Lake Zug with its surrounding fields, the Rigi mountain in the center of the test site, and the Lake of Four Cantons. The region is covered by the AVIRIS flight #910705, run 6 and 7 of the NASA ER-2 aircraft resulting in an average nominal pixel size of about 18 m. Simultaneous to the ER-2 overflight spectroradiometric measurements have been taken in various locations. Preselected reference targets were measured in the field with a GER Mark V spectroradiometer, and radiance measurements were taken to the lake using a Li-Cor LI 1800UW specroradiometer below and above the water surface. A comprehensive meteorological data set was obtained by joining the POLLUMET experiment which carried out measurements to investigate the summer smog in Switzerland on the same day. The quality assessment for the actual data set can be found in detail in Meyer et al. A parametric approach calculating the location of the airplane was used to simulate the observation geometry. This parametric preprocessing procedure, which takes care of effects of flight line and attitude variations as well as the pixel-by-pixel topographic corrections is described in Meyer.

Description: Roads and highways show up clearly in many bands of AVIRIS images. A typical lane in the U.S. is 12 feet wide, and the total width of a four lane highway, including 18 feet of paved shoulders, is 19.8 m. Such a highway will cover only a portion of any 20x20 m AVIRIS pixel that it traverses. The other portion of these pixels wil be usually covered by vegetation. An interesting problem is to precisely determine the location of a highway within the AVIRIS pixels that it traverses. This information may be used for alignment and spatial calibration of AVIRIS images. Also, since the reflection properties of highway surfaces do not change with time, and they can be determined once and for all, such information can be of help in calculating and filtering out the atmospheric noise that contaminates AVIRIS measurements. The purpose of this report is to describe a method for sub-pixel localization of highways.

Description: Heterogeneity in ecological phenomena are scale dependent and affect the hierarchical structure of image data. AVIRIS pixels average reflectance produced by complex absorption and scattering interactions between biogeochemical composition, canopy architecture, view and illumination angles, species distributions, and plant cover as well as other factors. These scales affect validation of pixel reflectance, typically performed by relating pixel spectra to ground measurements acquired at scales of 1m(exp 2) or less (e.g., field spectra, foilage and soil samples, etc.). As image analysis becomes more sophisticated, such as those for detection of canopy chemistry, better validation becomes a critical problem. This paper presents a methodology for bridging between point measurements and pixels using geostatistics. Geostatistics have been extensively used in geological or hydrogeolocial studies but have received little application in ecological studies. The key criteria for kriging estimation is that the phenomena varies in space and that an underlying controlling process produces spatial correlation between the measured data points. Ecological variation meets this requirement because communities vary along environmental gradients like soil moisture, nutrient availability, or topography.

Description: A simple mixing model employing reference endmembers (green vegetation, non-photosynthetic vegetation, soil and shade), and using 180 AVIRIS bands, was used to establish an interpretive framework for a forested area in the Pacific Northwest. A regrowth trend, based on changes in the endmember proportions, was defined for conifers that extends from clearcuts to mature forest, and by implication to old growth. Deciduous species within replanted forest plots caused the fractions to be displaced from the main coniferous regrowth trend and to move toward the green vegetation fraction. The results indicate that the spectral information in AVIRIS can be inverted to estimate approximate stand age and relative proportion of deciduous species in the context of the area studied. Using AVIRIS we measured a 3 to 5 percent increase in woody material in old-growth forest, as distinct from other mature forest. This result is consistent with a predicted increase in NPV in old-growth forest, based on field observations. Previous application of the mixing analysis to a TM image of the same area separated old growth based solely on the shade fraction; however the approach required successful removal of shade introduced by topography. Our new results suggest that with the high spectral resolution and high signal-to-noise of AVIRIS images it may be possible to characterize and map old-growth forests in the Northwest using both the NPV fraction and shade.

Description: Spectral Mixture Analysis (SMA) has become a well established procedure for analyzing imaging spectrometry data, however, the technique is relatively insensitive to minor sources of spectral variation (e.g., discriminating stressed from unstressed vegetation and variations in canopy chemistry). Other statistical approaches have been tried e.g., stepwise multiple linear regression analysis to predict canopy chemistry. Grossman et al. reported that SMLR is sensitive to measurement error and that the prediction of minor chemical components are not independent of patterns observed in more dominant spectral components like water. Further, they observed that the relationships were strongly dependent on the mode of expressing reflectance (R, -log R) and whether chemistry was expressed on a weight (g/g) or are basis (g/sq m). Thus, alternative multivariate techniques need to be examined. Smith et al. reported a revised SMA that they termed Foreground/Background Analysis (FBA) that permits directing the analysis along any axis of variance by identifying vectors through the n-dimensional spectral volume orthonormal to each other. Here, we report an application of the FBA technique for the detection of canopy chemistry using a modified form of the analysis.

Description: Surface albedo and snow-covered-area (SCA) are crucial inputs to the hydrologic and climatologic modeling of alpine and seasonally snow-covered areas. Because the spectral albedo and thermal regime of pure snow depend on grain size, areal distribution of snow grain size is required. Remote sensing has been shown to be an effective (and necessary) means of deriving maps of grain size distribution and snow-covered-area. Developed here is a technique whereby maps of grain size distribution improve estimates of SCA from spectral mixture analysis with AVIRIS data.

Description: Experimental results have shown the existence of a strong relationship between chlorophyll alpha concentration and remote sensing reflectance measured at lake level with a high resolution spectroradiometer. The objective of our study was to investigate the relationship between surface chlorophyll alpha concentration at Mono Lake and water reflectance retrieved from Airborne Visible - Infrared Imaging Spectrometer (AVIRIS) data obtained in october 7, 1992. AVIRIS data were atmospherically corrected as described by Green et al. A description of the lake-level sampling is found in Melack and Gastil. The relationship between chlorophyll concentration and both the single band reflectance and the first difference transformation of the reflectance spectra for the first 40 AVIRIS spectral bands (400 nm to 740 nm) was examined. The relationship was then used to produce a map of the surface chlorophyll distribution.

Description: We have demonstrated the unique utility of imaging spectroscopy in mapping mineral distribution. In the Summitville mining region we have shown that the mine site does not contribute clay minerals to the Alamosa River, but does contribute Fe-bearing minerals. Such minerals have the potential to carry heavy metals. This application illustrates only one specific environmental application of imaging spectroscopy data. For instance, the types of minerals we can map with confidence are those frequently associated with environmental problems related to active and abandoned mine lands. Thus, the potential utility of this technology to the field of environmental science has yet to be fully explored.

Description: Melting of the snowpack is a critical parameter that drives aspects of the hydrology in regions of the earth where snow accumulates seasonally. New techniques for measurement of snow melt over regional scales offer the potential to improve monitoring and modeling of snow-driven hydrological processes. We present the results of measuring the spectral absorption of liquid water in a melting snowpack with the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). AVIRIS data were acquired over Mammoth Mountain, in east central California on 21 May 1994 at 18:35 UTC. The air temperature at 2926 m on Mammoth Mountain at site A was measured at 15-minute intervals during the day preceding the AVIRIS data acquisition. At this elevation, the air temperature did not drop below freezing the night of May 20 and had risen to 6 degrees Celsius by the time of the overflight on May 21. These temperature conditions support the presence of melting snow at the surface as the AVIRIS data were acquired.

Description: The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) acquired data as part of the Boreal Ecosystem-Atmosphere Study (BOREAS) in 1994. Flights occurred over the northern study area (NSA) in the region of 56 degrees north latitude and 98.5 degrees west longitude and over the southern study area (SSA) at 54 degrees north latitude and 105 degrees west longitude. These data will be used to directly derive spectral properties of the surface and atmosphere and to provide supporting data for other instruments, models, and experiments in support of the BOREAS objectives. We present a preliminary evaluation of the AVIRIS data collected in BOREAS in terms of the AVIRIS-derived parameters: water vapor, leaf water, and apparent spectral reflectance.

Description: The remote estimation of leaf biochemical content from spaceborne platforms has been the subject of many studies aimed at better understanding of terrestrial ecosystem functioning. The major ecological processes involved in exchange of matter and energy, like photosynthesis, primary production, evaportranspiration, respiration, and decomposition can be related to plant properties e.g., chlorophyll, water, protein, cellulose and lignin contents. As leaves represent the most important plant surfaces interacting with solar energy, a top priority has been to relate optical properties to biochemical constituents. Two different approaches have been considered: first, statistical correlations between the leaf reflectance (or transmittance) and biochemical content, and second, physically based models of leaf scattering and absorption developed using the laws of optics. Recently reviewed by Verdebout et al., the development of models of leaf optical properties has resulted in better understanding of the interaction of light with plant leaves. Present radiative transfer models mainly use chlorophyll and/or water contents as input parameters to calculate leaf reflectance. Inversion of these models allows to retrieve these constituents from spectrophotometric measurements. Conel et al. recently proposed a two-stream Kubelka-Munk model to analyze the influence of protein, cellulose, lignin, and starch on leaf reflectance, but in fact, the estimation of leaf biochemistry from remote sensing is still an open question. In order to clarify it, a laboratory experiment associating visible/infrared spectra of plan leaves both with physical measurements and biochemical analyses was conducted at the Joint Research Center during the summer of 1993. This unique data set has been used to upgrade the PROSPECT model, by including leaf biochemistry.

Description: Plant species composition and plant architectural attributes are critical parameters required for the measuring, monitoring, and modeling of terrestrial ecosystems. Remote sensing is commonly cited as an important tool for deriving vegetation properties at an appropriate scale for ecosystem studies, ranging from local to regional and even synoptic scales. Classical approaches rely on vegetation indices such as the normalized difference vegetation index (NDVI) to estimate biophysical parameters such as leaf area index or intercepted photosynthetically active radiation (IPAR). Another approach is to apply a variety of classification schemes to map vegetation and thus extrapolate fine-scale information about specific sites to larger areas of similar composition. Imaging spectrometry provides additional information that is not obtainable through broad-band sensors and that may provide improved inputs both to direct biophysical estimates as well as classification schemes. Some of this capability has been demonstrated through improved discrimination of vegetation, estimates of canopy biochemistry, and liquid water estimates from vegetation. We investigate further the potential of leaf water absorption estimated from Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data as a means for discriminating vegetation types and deriving canopy architectural information. We expand our analysis to incorporate liquid water estimates from two spectral regions, the 1000-nm region and the 2200-nm region. The study was conducted in the vicinity of Jasper Ridge, California, which is located on the San Francisco peninsula to the west of the Stanford University campus. AVIRIS data were acquired over Jasper Ridge, CA, on June 2, 1992, at 19:31 UTC. Spectra from three sites in this image were analyzed. These data are from an area of healthy grass, oak woodland, and redwood forest, respectively. For these analyses, the AVIRIS-measured upwelling radiance spectra for the entire Jasper Ridge scene were transformed to apparent surface reflectance using a radiative transfer code-based inversion algorithm.

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: Using spectral imaging data acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS) from an ER-2 aircraft at 20 km altitude during various field programs, it was found that narrow channels near the center of the strong 1.38-micrometer water vapor band are very effective in detecting think cirrus clouds. Based on this observation from AVIRIS data, Gao and Kaufman proposed to put a channel centered at 1.375 micrometers with a width of 30 nm on the Moderate Resolution Imaging Spectrometer (MODIS) for remote sensing of cirrus clouds from space. The sensitivity of the 1.375-micrometer MODIS channel to detect thin cirrus clouds during the day time is expected to be one to two orders of magnitude better than the current infrared emission techniques. As a result, much larger fraction of the satellite data is expected to be identified as being covered by cirrus clouds, some of them so thin that their obscuration of the surface is very small. In order to make better studies of surface reflectance properties, thin cirrus effects must be removed from satellite images. Therefore, there is a need to study radiative properties of thin cirrus clouds, so that a strategy for correction or removal of the thin cirrus effects, similar to the correction of atmospheric aerosol effect, can be formed. In this extended abstract, we describe an empirical approach for removing/correcting thin cirrus effects in AVIRIS images using channels near 1.375 microns - one step beyond the detection of cirrus clouds using these channels.

Description: During the past two years, researchers from several institutes joined together to take part in two SCAR experiments. The SCAR-A (Sulfates, Clouds And Radiation - Atlantic) took place in the mid-Atlantic region of the United States in July, 1993. remote sensing data were acquired with the Airborne Visible Infrared Imaging Spectrometer (AVIRIS), the MODIS Airborne Simulator (MAS), and a RC-10 mapping camera from an ER-2 aircraft at 20 km. In situ measurements of aerosol and cloud microphysical properties were made with a variety of instruments equipped on the University of Washington's C-131A research aircraft. Ground based measurements of aerosol optical depths and particle size distributions were made using a network of sunphotometers. The main purpose of SCAR-A experiment was to study the optical, physical and chemical properties of sulfate aerosols and their interaction with clouds and radiation. Sulfate particles are believed to affect the energy balance of the earth by directly reflecting solar radiation back to space and by increasing the cloud albedo. The SCAR-C (Smoke, Clouds And Radiation - California) took place on the west coast areas during September - October of 1994. Sets of aircraft and ground-based instruments, similar to those used during SCAR-A, were used during SCAR-C. Remote sensing of fires and smoke from AVIRIS and MAS imagers on the ER-2 aircraft was combined with a complete in situ characterization of the aerosol and trace gases from the C-131A aircraft of the University of Washington and the Cesna aircraft from the U.S. Forest Service. The comprehensive data base acquired during SCAR-A and SCAR-C will contribute to a better understanding of the role of clouds and aerosols in global change studies. The data will also be used to develop satellite remote sensing algorithms from MODIS on the Earth Observing System.

Description: In this paper, we evaluate the potential for extracting the 'photochemical reflectance index' (PRI; previously called the 'physiological reflectance index') from AVIRIS data. This index, which is derived from narrow-band reflectance at 531 and 570 nm, has proven to be a useful indicator of photosynthetic function at the leaf and canopy scales. At the leaf level, PRI varies with photosynthetic capacity, radiation-use efficiency, and vegetation type (unpublished data). This finding is consistent with the hypothesis that vegetation types exhibiting chronically reduced photosynthesis during periods of stress (e.g. drought-tolerant evergreens) invest proportionally more in photoprotective processes than vegetation with high photosynthetic capacity (e.g. crops or deciduous perennials). Vertical transects in tropical and boreal forest canopies have indicated declines in PRI associated with downregulation of photosynthesis at the canopy tops under sunny, dry midday conditions (unpublished data). This reduced PRI in upper canopy levels provides a further basis for examining this signal with the 'view from above' afforded by aircraft overflights. Although many factors could confound interpretation of a subtle physiological signal at the landscape scale, we conducted a preliminary examination of PRI extracted from existing, AVIRIS imagery of Stanford University's Jasper Ridge Biological Preserve obtained on the June 2nd, 1992, overflight. The goal was to use the hyperspectral capabilities of AVIRIS to evaluate the potential of this index for obtaining useful physiological data at the landscape scale. The expectation based on leaf- and canopy-level studies was that regions containing vegetation of reduced photosynthetic capacity (e.g. chaparral or evergreen woodland) would exhibit lower PRI values than regions of high capacity (e.g. deciduous woodland).

Description: One of the challenges of Imaging Spectroscopy is the identification, mapping and abundance determination of materials, whether mineral, vegetable, or liquid, given enough spectral range, spectral resolution, signal to noise, and spatial resolution. Many materials show diagnostic absorption features in the visual and near infrared region (0.4 to 2.5 micrometers) of the spectrum. This region is covered by the modern imaging spectrometers such as AVIRIS. The challenge is to identify the materials from absorption bands in their spectra, and determine what specific analyses must be done to derive particular parameters of interest, ranging from simply identifying its presence to deriving its abundance, or determining specific chemistry of the material. Recently, a new analysis algorithm was developed that uses a digital spectral library of known materials and a fast, modified-least-squares method of determining if a single spectral feature for a given material is present. Clark et al. made another advance in the mapping algorithm: simultaneously mapping multiple minerals using multiple spectral features. This was done by a modified-least-squares fit of spectral features, from data in a digital spectral library, to corresponding spectral features in the image data. This version has now been superseded by a more comprehensive spectral analysis system called Tricorder.

Description: The success of imaging spectrometry in mineralogic mapping of natural terrains indicates that the technology can also be used to assess the environmental impact of human activities in certain instances. Specifically, this paper describes an investigation into the use of data from the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) for mapping the spread of, and assessing changes in, the mineralogic character of tailings from a major silver and base metal mining district. The area under investigation is the Coeur d'Alene River Valley in northern Idaho. Mining has been going on in and around the towns of Kellogg and Wallace, Idaho since the 1880's. In the Kellogg-Smelterville Flats area, west of Kellogg, mine tailings were piled alongside the South Fork of the Coeur d'Alene River. Until the construction of tailings ponds in 1968 much of these waste materials were washed directly into the South Fork. The Kellogg-Smelterville area was declared an Environmental Protection Agency (EPA) Superfund site in 1983 and remediation efforts are currently underway. Recent studies have demonstrated that sediments in the Coeur d'Alene River and in the northern part of Lake Coeur d'Alene, into which the river flows, are highly enriched in Ag, Cu, Pb, Zn, Cd, Hg, As, and Sb. These trace metals have become aggregated in iron oxide and oxyhydroxide minerals and/or mineraloids. Reflectance spectra of iron-rich tailing materials are shown. Also shown are spectra of hematite and goethite. The broad bandwidth and long band center (near 1 micron) of the Fe(3+) crystal-field band of the iron-rich sediment samples combined with the lack of features on the Fe(3+) -O(2-) charge transfer absorption edge indicates that the ferric oxide and/or oxyhydroxide in these sediments is poorly crystalline to amorphous in character. Similar features are seen in poorly crystalline basaltic weathering products (e.g., palagonites). The problem of mapping and analyzing the downriver occurrences of iron rich tailings in the Coeur d'Alene (CDA) River Valley using remotely sensed data is complicated by the full vegetation cover present in the area. Because exposures of rock and soil were sparse, the data processing techniques used in this study were sensitive to detecting materials at subpixel scales. The methods used included spectral mixture analysis and a constrained energy minimization technique.

Description: A major goal of airborne imaging spectrometry is to estimate the biochemical composition of vegetation canopies from reflectance spectra. Remotely-sensed estimates of foliar biochemical concentrations of forests would provide valuable indicators of ecosystem function at regional and eventually global scales. Empirical research has shown a relationship exists between the amount of radiation reflected from absorption features and the concentration of given biochemicals in leaves and canopies (Matson et al., 1994, Johnson et al., 1994). A technique commonly used to determine which wavelengths have the strongest correlation with the biochemical of interest is unguided (stepwise) multiple regression. Wavelengths are entered into a multivariate regression equation, in their order of importance, each contributing to the reduction of the variance in the measured biochemical concentration. A significant problem with the use of stepwise regression for determining the correlation between biochemical concentration and spectra is that of 'overfitting' as there are significantly more wavebands than biochemical measurements. This could result in the selection of wavebands which may be more accurately attributable to noise or canopy effects. In addition, there is a real problem of collinearity in that the individual biochemical concentrations may covary. A strong correlation between the reflectance at a given wavelength and the concentration of a biochemical of interest, therefore, may be due to the effect of another biochemical which is closely related. Furthermore, it is not always possible to account for potentially suitable waveband omissions in the stepwise selection procedure. This concern about the suitability of stepwise regression has been identified and acknowledged in a number of recent studies (Wessman et al., 1988, Curran, 1989, Curran et al., 1992, Peterson and Hubbard, 1992, Martine and Aber, 1994, Kupiec, 1994). These studies have pointed to the lack of a physical link between wavelengths chosen by stepwise regression and the biochemical of interest, and this in turn has cast doubts on the use of imaging spectrometry for the estimation of foliar biochemical concentrations at sites distant from the training sites. To investigate this problem, an analysis was conducted on the variation in canopy biochemical concentrations and reflectance spectra using forced entry linear regression.

Description: We analyzed AVIRIS data obtained over agricultural areas in the San Luis Valley of Colorado. The data were acquired on September 3, 1993. A combined method of radiative transfer modeling and ground calibration site reflectance was used to correct the flight data to surface reflectance. This method, called Radiative Transfer Ground Calibration, or RTGC, corrects for variable water vapor in the atmosphere and produces spectra free of artifacts with spectral channel to channel noise approaching the signal to noise of the raw data. The calibration site soil samples were obtained on the day of the overflight and measured on our laboratory spectrometer. The site was near the center of the AVIRIS scene and the spectra of the soil is spectrally bland, especially in the region of the chlorophyll absorption in the visible portion of the spectrum. The center of the scene is located at approximately 106 deg 03' longitude, 37 deg 23' latitude, and the scene covers about 92 square kilometers. This scene is one of 28 in the area for a general project to study the Summitville abandoned mine site, located in the mountains west of the San Luis Valley, and its effects on the surrounding environment.

Description: There is a user need for increasing spatial and spectral resolution in Earth Observation (EO) optical instrumentation. Higher spectral resolution will be achieved by the introduction of spaceborne imaging spectrometers. Higher spatial resolutions of 1 - 3m will be achieved also, but at the expense of sensor redesign, higher communications bandwidth, high data processing volumes, and therefore, at the risk of time delays due to large volume data-handling bottlenecks. This paper discusses a design concept whereby the hyperspectral properties of a spaceborne imaging spectrometer can be used to increase the image spatial resolution, without such adverse cost impact.

Description: A complete spectral unmixing of a complicated AVIRIS scene may not always be possible or even desired. High quality data of spectrally complex areas are very high dimensional and are consequently difficult to fully unravel. Partial unmixing provides a method of solving only that fraction of the data inversion problem that directly relates to the specific goals of the investigation. Many applications of imaging spectrometry can be cast in the form of the following question: 'Are my target signatures present in the scene, and if so, how much of each target material is present in each pixel?' This is a partial unmixing problem. The number of unmixing endmembers is one greater than the number of spectrally defined target materials. The one additional endmember can be thought of as the composite of all the other scene materials, or 'everything else'. Several workers have proposed partial unmixing schemes for imaging spectrometry data, but each has significant limitations for operational application. The low probability detection methods described by Farrand and Harsanyi and the foreground-background method of Smith et al are both examples of such partial unmixing strategies. The new method presented here builds on these innovative analysis concepts, combining their different positive attributes while attempting to circumvent their limitations. This new method partially unmixes AVIRIS data, mapping apparent target abundances, in the presence of an arbitrary and unknown spectrally mixed background. It permits the target materials to be present in abundances that drive significant portions of the scene covariance. Furthermore it does not require a priori knowledge of the background material spectral signatures. The challenge is to find the proper projection of the data that hides the background variance while simultaneously maximizing the variance amongst the targets.

Description: Starting in 1994, all AVIRIS data distributions include a new product useful for quantification and modeling of the noise in the reported radiance data. The 'postcal' file contains approximately 100 lines of dark current data collected at the end of each data acquisition run. In essence this is a regular spectral-image cube, with 614 samples, 100 lines and 224 channels, collected with a closed shutter. Since there is no incident radiance signal, the recorded DN measure only the DC signal level and the noise in the system. Similar dark current measurements, made at the end of each line are used, with a 100 line moving average, to remove the DC signal offset. Therefore, the pixel-by-pixel fluctuations about the mean of this dark current image provide an excellent model for the additive noise that is present in AVIRIS reported radiance data. The 61,400 dark current spectra can be used to calculate the noise levels in each channel and the noise covariance matrix. Both of these noise parameters should be used to improve spectral processing techniques. Some processing techniques, such as spectral curve fitting, will benefit from a robust estimate of the channel-dependent noise levels. Other techniques, such as automated unmixing and classification, will be improved by the stable and scene-independence noise covariance estimate. Future imaging spectrometry systems should have a similar ability to record dark current data, permitting this noise characterization and modeling.

Description: Mineral-bound ammonium (NH4+) was discovered by the U.S. Geological Survey in the southern Cedar Mountains of Esmeralda County, Nevada in 1989. At 10 km in length, this site is 100 times larger than any previously known occurrence in volcanic rocks. The ammonium occurs in two hydrothermally altered, crystal-rich rhyolitic tuff units of Oligocene age, and is both structurally and stratigraphically controlled. This research uses Advanced Visible/Infrared Imaging Spectrometer (AVIRIS) data to quantitatively map the mineral-bound ammonium (buddingtonite) concentration in the altered volcanic rocks. Naturally occurring mineral-bound ammonium is fairly rare; however, it has been found to occur in gold-bearing hydrothermal deposits. Because of this association, it is thought that ammonium may be a useful too in exploration for gold and other metal deposits. Mineral-bound ammonium is produced when an ammonium ion (NH4+) replaces the alkali cation site (usually K+) in the crystal structure of silicate minerals such as feldspars, micas and clays. Buddingtonite is an ammonium feldspar. The ammonium originates in buried organic plant matter and is transported to the host rock by hydrothermal fluids. Ammonium alteration does not produce visible changes in the rock, and it is barely detectable with standard x-ray diffraction methods. It is clearly identified, however, by absorption features in short wave-infrared (SWIR) wavelengths (2.0 - 2.5 micrometers). The ammonium absorption features are believed to be caused by N-H vibrational modes and are analogous to hydroxyl (O-H) vibrational modes, only shifted slightly in wavelength. Buddingtonite absorption features in the near- and SWIR lie at 1.56, 2.02 and 2.12 micrometers. The feature at 2.12 micrometer is the strongest of the three and is the only one used in this study. The southern Cedar Mountains are sparsely vegetated and are an ideal site for a remote sensing study.

Description: Near infrared laboratory spectra have been used for many years to determine nitrogen and lignin concentrations in plant materials. In recent years, similar high spectral resolution visible and infrared data have been available via airborne remote sensing instruments. Using data from NASA's Airborne visible/Infrared Imaging Spectrometer (AVIRIS) we attempt to identify spectral regions correlated with foliar chemistry at the canopy level in temperate forests.

Description: MODTRAN, the Moderate Resolution Atmospheric Radiance and Transmittance Model, encompasses all the capabilities of LOWTRAN 7, the widely used 20 cm(exp -1) resolution radiance code, but incorporates a much more sensitive molecular band model with 2 cm(exp -1) resolution. MODTRAN contains many important elements that other band model based radiative transfer codes do not incorporate. It shares with FASCODE: spherical geometry, single and multiple scattering default atmospheric profile descriptors (gases, aerosols, clouds, fogs, and rain), and molecular continua (H2O, CO2, O3, O2, N2). In addition, it can calculate the solar/lunar direct and scattered radiation. MODTRAN3 was released to the general public in November 1994. It has several important features that the previous version, MODTRAN2, does not have. Chloro-fluorocarbon (CFC) and related heavy molecules (whose spectroscopic properties first appear on the HITRAN92 data base as temperature-dependent cross sections) have been incorporated into pseudo-band models, with provision for using both default and user supplied profiles. The addition of SO2 and O2 in the UV, along with upgraded ozone Chappuis bands in the visible is also part of MODTRAN3. An improved multiple scattering algorithm, the DIScrete Ordinate Radiative Transfer (DISORT) has also been incorporated into MODTRAN3. MODTRAN is very fast: simple timing runs of MODTRAN3 vs. FASCOD3 show an improvement of more than a factor of 100 for a typical 500 cm(exp -1) spectral interval and comparable vertical layering. Speed is an important consideration in heating/cooling rates calculations, where a large number of radiative transfer calculations are needed. The MODTRAN3 used in this study is based on HITRAN92, but as mentioned, above, it will be upgraded to HITRAN94 upon its release at the end of 1994. MODTRAN has been adopted by: some researchers in the AVIRIS program as one radiative transfer code to derive surface reflectance from AVIRIS measurements. The accuracy of the code is very important because any errors in the radiative transfer calculation will directly translate into errors in the derived surface reflectance. In this paper, the new solar irradiance calculated by Kurucz, which is adopted in MODTRAN3, will be presented. Recent validations of MODTRAN3 with airborne high resolution interferometer measurements over ocean will be discussed. Good agreeement between model calculations and measurements was achieved.