ALBERT

Description: The development of an integrated approach to the modeling of forest dynamics encompassing submodels of forest growth and succession, soil processes and radiation interactions, is reported. Remote sensing technology is a key element of this study in that it provides data for developing, initializing, updating, and validating the models. The objectives are reviewed, the data collected and models in use are discussed, and a framework for studying interactions between the forest growth, soil process and energy interaction components, is described. Remote sensing technology used in the study includes optical and microwave field, aircraft and satellite borne instruments. The types of data collected during intensive field and aircraft campaigns included bidirectional reflectance, thermal emittance and multifrequency, multipolarization synthetic aperture radar backscatter. Synthetic imagery of derived products such as forest biomass and NDVI (Normalized Difference Vegetative Index), and collections of ground data are being assembled in a georeferenced data base. These data are used to drive or test multidiscipline simulations of forested ecosystems. Enhancements to the modeling environment permit considerable flexibility in configuring simulations and selecting results for reporting and graphical display.

Description: The quantitative interpretation of satellite observations requires the use of mathematical tools to extract the desired information on terrestrial environments from the radiation data collected in space. A whole range of approaches can be pursued, from the development of models capable of explaining the nature of the physical signal being measured and of characterizing the state of the system under observation, to the empirical correlations between the variables of interest and the space measurements. The premises and implications of these approaches are outlined, paying special attention to the mathematical and numerical requirements. The role and specific applications of empirical bidirectional reflectance models is also discussed, even though these models do not contribute to the understanding of the theory of radiation transfer or to the assessment of the variables of interest. The advantages and drawbacks of these various approaches and the research priorities for the next few years are discussed in the context of the planned availability of new sensors.

Description: Management of crop residues, the portion of a crop left in the field after harvest, is an important conservation practice for minimizing soil erosion and for improving water quality. Quantification of crop residue cover is required to evaluate the effectiveness of conservation tillage practices. Methods are needed to quantify residue cover that are rapid, accurate, and objective. The fluorescence of crop residue was found to be a broadband phenomenon with emission maxima at 420 to 495 nm for excitations of 350 to 420 nm. Soils had low intensity broadband emissions over the 400 to 690 nm region for excitations of 300 to 600 nm. The range of relative fluorescence intensities for the crop residues was much greater than the fluorescence observed of the soils. As the crop residues decompose their blue fluorescence values approach the fluorescence of the soil. Fluorescence techniques are concluded to be less ambiguous and better suited for discriminating crop residues and soils than reflectance methods. If properly implemented, fluorescence techniques can be used to quantify, not only crop residue cover, but also photosynthetic efficiency in the field.

Description: The importance of the measurement of wind fields is discussed. Wind regime data can be used to infer the amount and type of wind induced (aerolian) transport of sand and dust, or to establish global circulation models, for example on other planets. Since local measurements are costly and often impossible, it is desired to infer such data from remotely sensed information. A potential mechanism for remotely inferring the wind regime by using synthetic aperture radar data to describe the roughness of the surface is described. A project to estimate the practicality of using such a mechanism is described. An experiment that extends the mechanism to vegetated sites, where the goal is to measure potential for erosion, is reported.

Description: Surface reflectance is required to quantitatively investigate molecular absorption and particle scattering properties of materials on the Earth's surface. Atmospheric aerosol optical depth, surface pressure and water vapor are required to constrain a radiative transfer code for the inversion of measured spectral radiance to apparent surface reflectance. A suite of algorithms using nonlinear least squares fitting techniques are described that directly estimate these atmospheric parameters from spectral radiance measured by the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS). The derived atmospheric parameters are used to constrain a radiative transfer code for the inversion of the imaging spectrometer radiance to apparent reflectance. The derived apparent reflectance is validated with respect to in situ measurement on the same target.

Description: Investigations designed to study land surface hydrologic-atmospheric interactions, showing the potential of L band passive microwave radiometry for measuring surface soil moisture over large areas, are discussed. Satisfying the data needs of these investigations requires the ability to map large areas rapidly. With aircraft systems this means a need for more beam positions over a wider swath on each flightline. For satellite systems the essential problem is resolution. Both of these needs are currently being addressed through the development and verification of Electronically Scanned Thinned Array Radiometer (ESTAR) technology. The ESTAR L band radiometer was evaluated for soil moisture mapping applications in two studies. The first was conducted over the semiarid rangeland Walnut Gulch watershed located in south eastern Arizona (U.S.). The second was performed in the subhumid Little Washita watershed in south west Oklahoma (U.S.). Both tests showed that the ESTAR is capable of providing soil moisture with the same level of accuracy as existing systems.

Description: A weather resistant automatic scanning Sun photometer system is assessed and demonstrated as practical for measurements of aerosol concentrations and properties at remote sites. Interfaced with a transmitter using the Geostationary Data Collection System (GDCS), the data are processed in near real time. The processing allows a time dependence of the aerosols and water vapor and an ongoing assessment of the health and calibration of the instruments. The system's automatic data acquisition, transmission, and processing offer immediate application to atmospheric monitoring and modeling on a regional to global scale and validation of satellite retrievals. It is estimated that under normal circumstances the retrieved aerosol optical thickness has a network wide accuracy of +/- 0.02 from 340 nm to 1020 nm, water vapor +/- 0.2 cm and size distribution from 0.1 to 3 micrometers.

Description: As part of a global program to validate the ocean surface sensors on board ERS-1, a joint experiment on the Grand Banks of Newfoundland was carried out in Nov. 1991. The principal objective was to provide a field validation of ERS-1 Synthetic Aperture Radar (SAR) measurement of ocean surface structure. The NASA-P3 aircraft measurements made during this experiment provide independent measurements of the ocean surface along the validation swath. The Radar Ocean Wave Spectrometer (ROWS) is a radar sensor designed to measure direction of the long wave components using spectral analysis of the tilt induced radar backscatter modulation. This technique greatly differs from SAR and thus, provides a unique set of measurements for use in evaluating SAR performance. Also, an altimeter channel in the ROWS gives simultaneous information on the surface wave height and radar mean square slope parameter. The sets of geophysical parameters (wind speed, significant wave height, directional spectrum) are used to study the SAR's ability to accurately measure ocean gravity waves. The known distortion imposed on the true directional spectrum by the SAR imaging mechanism is discussed in light of the direct comparisons between ERS-1 SAR, airborne Canadian Center for Remote Sensing (CCRS) SAR, and ROWS spectra and the use of the nonlinear ocean SAR transform.

Description: Terrain slopes, which can be measured with Synthetic Aperture Radar (SAR) interferometry either from a height map or from the interferometric phase gradient, were used to calculate the local incidence angle and the correct pixel area. Both are required for correct thematic interpretation of SAR data. The interferometric correlation depends on the pixel area projected on a plane perpendicular to the look vector and requires correction for slope effects. Methods for normalization of the backscatter and interferometric correlation for ERS-1 SAR are presented.

Description: The relationship between the gradient of the interferometric phase and the terrain slope, which, it is thought, would allow a derivation of the terrain slopes without phase unwrapping, is presented. A linear relationship between the interferometric phase gradient and the terrain slopes was found. A quantitative error analysis showed that only very small errors are introduced by these approximations for orbital Synthetic Aperture Radar (SAR) geometries. An example of a slope map for repeat pass interferometry from ERS-1 SAR data is given. A number of direct and indirect applications of the terrain slope are indicated: erosion and avalanche hazard studies, radiometric calibration of SAR data, and normalization of the interferometric correlation coefficient.

Description: Research on the use of active microwaves in remote sensing, presented during plenary and poster sessions, is summarized. The main highlights are: calibration techniques are well understood; innovative modeling approaches have been developed which increase active microwave applications (segmentation prior to model inversion, use of ERS-1 scatterometer, simulations); polarization angle and frequency diversity improves characterization of ice sheets, vegetation, and determination of soil moisture (X band sensor study); SAR (Synthetic Aperture Radar) interferometry potential is emerging; use of multiple sensors/extended spectral signatures is important (increase emphasis).

Description: Crop residues, the portion of the crop left in the field after harvest, can be an important management factor in controlling soil erosion. Methods to quantify residue cover are needed that are rapid, accurate, and objective. Scenes with known amounts of crop residue were illuminated with long wave ultraviolet (UV) radiation and fluorescence images were recorded with an intensified video camera fitted with a 453 to 488 nm band pass filter. A light colored soil and a dark colored soil were used as background for the weathered soybean stems. Residue cover was determined by counting the proportion of the pixels in the image with fluorescence values greater than a threshold. Soil pixels had the lowest gray levels in the images. The values of the soybean residue pixels spanned nearly the full range of the 8-bit video data. Classification accuracies typically were within 3(absolute units) of measured cover values. Video imaging can provide an intuitive understanding of the fraction of the soil covered by residue.

Description: Global study of land surface properties uses AVHRR channels 1 and 2, but channel 3 may be of interest, although its use requires preprocessing. It consists of both a reflective part and an emissive part, the former can be derived from T3, T4 and T5. Since the water vapor affects channel 3, its content is retrieved from the channel 4 and 5 using the split window technique. A formula of reflective part retrieval at 3.75 micrometers is tested in the case of sunglint observations where the emissivities of channels 4 and 5 can be set to the unity. The formula is adapted and validated to land surface using the FIFE-87 data set. Preliminary applications of the reflectance at 3.75 micrometers to the studies of surface properties retrieval, aerosol retrieval over land, and desertic aerosol retrieval, are addressed.

Description: The polarization of the sunlight scattered by atmospheric aerosols or cloud droplets and reflected from ground surfaces or plant canopies may convey much information when used for remote sensing purposes. The typical polarization features of aerosols, cloud droplets, and plant canopies, as observed by ground based and airborne sensors, are investigated, looking especially for those invariant properties amenable to description by simple models when possible. The question of polarization measurements from space is addressed. The interest of such measurements for remote sensing purposes is investigated, and their feasibility is tested by using results obtained during field campaigns of the airborne POLDER instrument, a radiometer designed to measure the directionality and polarization of the sunlight scattered by the ground atmosphere system.

Description: Papers focused on land surface, atmospheric, and ocean properties are reported. Specific comments pertaining to polarization, models and inversion, and measurements, are given. Recommendations are: continued research into the application potential of the BRDF (Bidirectional Reflectance Distribution Function) and polarization properties of ground surface and atmospheric targets; three dimensional models, which account for the statistical behavior of remotely sensed data, should be extended and inverted in order to support analysis of data potentially covering rolling terrain such that pixels represent heterogeneous mixtures of surface cover types and project ground footprints with sizes between 10 to 6 km, the ground pixel sizes of planned future sensors; available reflectance models should be further validated by means of multi dimensional (directional, spectral, temporal) field data and existing models should be intercompared in more depth to evaluate their performance and limitations; existing methods for model inversion should be validated in more depth in order to quantify the practical limitations and the expected accuracy of the parameters retrieved and new approaches should be developed based upon apriori knowledge of plant canopy development and spectral BRDF properties; there is a need to establish a protocol of validation and intercomparison of the indices and compositing techniques which have been proposed during these last years.

Description: Traditionally, the remote sensing community has relied totally on spectral knowledge to extract vegetation characteristics. However, there are other knowledge bases (KB's) that can be used to significantly improve the accuracy and robustness of inference techniques. Using AI (artificial intelligence) techniques a KB system (VEG) was developed that integrates input spectral measurements with diverse KB's. These KB's consist of data sets of directional reflectance measurements, knowledge from literature, and knowledge from experts which are combined into an intelligent and efficient system for making vegetation inferences. VEG accepts spectral data of an unknown target as input, determines the best techniques for inferring the desired vegetation characteristic(s), applies the techniques to the target data, and provides a rigorous estimate of the accuracy of the inference. VEG was developed to: infer spectral hemispherical reflectance from any combination of nadir and/or off-nadir view angles; infer percent ground cover from any combination of nadir and/or off-nadir view angles; infer unknown view angle(s) from known view angle(s) (known as view angle extension); and discriminate between user defined vegetation classes using spectral and directional reflectance relationships developed from an automated learning algorithm. The errors for these techniques were generally very good ranging between 2 to 15% (proportional root mean square). The system is designed to aid scientists in developing, testing, and applying new inference techniques using directional reflectance data.

Description: Most earth surfaces, particularly those supporting natural vegetation ecosystems, constitute structurally and spectrally complex surfaces that are distinctly non-Lambertian reflectors. Obtaining meaningful measurements of the directional radiances of landscapes and obtaining estimates of the complete bidirectional reflectance distribution functions of ground targets with complex and variable landscape and radiometric features are challenging tasks. Reasons for the increased interest in directional radiance measurements are presented, and the issues that must be addressed when trying to acquire directional radiances for vegetated land surfaces from different types of remote sensing platforms are discussed. Priority research emphases are suggested, concerning field measurements of directional surface radiances and reflectances for future research. Primarily, emphasis must be given to the acquisition of more complete and directly associated radiometric and biometric parameter data sets that will empower the exploitation of the 'angular dimension' in remote sensing of vegetation through enabling the further development and rigorous validation of state of the art plant canopy models.

Description: Synthetic aperture radar (SAR) images of the Greenland ice sheet collected by an airborne system clearly reveal the four melting facies of this sheet defined 30 years ago from snow stratigraphy studies by glaciologists. In particular, the radar echoes from the percolation facies have radiometric and polarimetric characteristics that are unique among terrestrial surfaces, but that resemble the exotic radar echoes recorded from the icy Galilean satellites. There, the radar signals interact with subsurface, massive ice features created in the cold, dry snow by seasonal melting and refreezing events. The subsurface features act as efficient reflectors of the incident radiation most likely via internal reflections. In the soaked-snow facies, the radar reflectivity is much lower because radar signals are attenuated by the wetter snow before they can interact with subsurface structures. Inversion algorithms to derive geophysical information from the SAR data are developed in both cases to estimate snow wetness in the soaked-snow facies and the mass of ice water retained in the percolation facies.

Description: An operational stratospheric correction scheme used after the Mount Pinatubo (Phillipines) eruption (Jun. 1991) is presented. The stratospheric aerosol distribution is assumed to be only variable with latitude. Each 9 days the latitudinal distribution of the optical thickness is computed by inverting radiances observed in the NOAA AVHRR channel 1 (0.63 micrometers) and channel 2 (0.83 micrometers) over the Pacific Ocean. This radiance data set is used to check the validity of model used for inversion by checking consistency of the optical thickness deduced from each channel as well as optical thickness deduced from different scattering angles. Using the optical thickness profile previously computed and radiative transfer code assuming Lambertian boundary condition, each pixel of channel 1 and 2 are corrected prior to computation of NDVI (Normalized Difference Vegetation Index). Comparison between corrected, non corrected, and years prior to Pinatubo eruption (1989 to 1990) NDVI composite, shows the necessity and the accuracy of the operational correction scheme.

Description: Aspects of aerosol studies and remote sensing are reviewed. Aerosol scatters solar radiation before it reaches the surface and scatters and absorbs it again after it is reflected from the surface and before it reaches the satellite sensor. The effect is spectrally and spatially dependent. Therefore atmospheric aerosol (dust, smoke and air pollution particles) has a significant effect on remote sensing. Correction for the aerosol effect was never achieved on an operational basis though several case studies were demonstrated. Correction can be done in a direct way by deriving the aerosol loading from the image itself and correcting for it using the appropriate radiative transfer model or by an indirect way, by defining remote sensing functions that are less dependent on the aerosol loading. To some degree this was already achieved in global remote sensing of vegetation where a composite of several days of NDVI (Normalized Difference Vegetation Index) measurements, choosing the maximal value, was used instead of a single cloud screened value. The Atmospheric Resistant Vegetation Index (ARVI) introduced recently for the NASA Earth Observing System EOS-MODIS is the most appropriate example of indirect correction, where the index is defined in such a way that the atmospheric effect in the blue spectral channel cancels to a large degree the atmospheric in the red channel in computations of a vegetation index. Atmospheric corrections can also use aerosol climatology and ground based instrumentation.

Description: A session dedicated to high spectral resolution in the solar spectrum, covering topics of calibration, atmospheric correction, geology/pedology, inland water, and vegetation, is reported. The session showed a high degree of diversity in the topics and the approaches used. It was highlighted that high spectral resolution data could provide atmospherically corrected ground level calibrated reflectance values. Important advances were shown in the use of radiative transfer models applied either on water bodies or vegetation. Several studies highlighted the high degree of redundancy contained in high spectral resolution data.

Description: A 54.95-MHz coherent backscatter radar, an ionosonde and the magnetometer located at Trivandrum in India (8.5 deg N, 77 deg E, 0.5 deg N dip angle) recorded large-amplitude ionospheric fluctuations and magnetic field fluctuations associated with a Pc5 micropulsation event, which occurred during an intense magnetic storm on 24 March 1991 (A(sub p) = 161). Simultaneous 100-n T-level fluctuations are also observed in the H-component at Brorfelde, Denmark (55.6 deg N gm) and at Narsarsuaq, Greenland (70.6 deg N gm). Our study of the above observations shows that the E-W electric field fluctuations in the E- and F-regions and the magnetic field fluctuations at Thumba are dominated by a near-sinusoidal oscillation of 10 min during 1730-1900 IST (1200-1330 UT), the amplitude of the electric field oscillation in the equatorial electrojet (EEJ) is 0.1-0.25 mV/m and it increases with height, while it is about 1.0 mV/m in the F-region, the ground-level H-component oscillation can be accounted for by the ionospheric current oscillation generated by the observed electric field oscillation in the EEJ and the H-component oscillations at Trivandrum and Brofelde are in phase with each other. The observations are interpreted in terms of a compressional cavity mode resonance in the inner magnetosphere and the assoicated ionospheric electric field penetrating from high latitudes to the magnetic equator.

Description: The 3.0-micrometers water of hydration absorption feature observed in the IR photometry of many low-albedo and some medium-albedo asteroids strongly correlates with the 0.7-micrometers Fe(+2) to Fe(+3) oxidized iron absorption feature observed in narrowband spectrophotometry of these asteroids. Using this relationship, an empirical algorithm for predicting the presence of water of hydration in the surface material of a Solar System body using photometry obtained through the Eight-Color Asteroid Survey nu (0.550 micrometers), w (0.701 micrometers), and x (0.853 micrometers) filters was developed and applied to the ECAS photometry of asteroids and outer planet satellites. The percentage of objects in low-albedo, outer main-belt asteroid classes that test positively for water of hydration increases from P to B to C to G class and correlates linearly with the increasing mean albedos of those objects testing positively. The medium-albedo M-class asteroids do not test positively in large number using this algorithm. Aqueously altered asteroids dominate the Solar System population between heliocentric distances of 2.6 to 3.5 AU, bracketing the Solar System region where the aqueous alteration mechanism operated most strongly. One jovian satellite, J VI Himalia, and one saturnian satellite. Phoebe, tested positively for water of hydration, supporting the hypothesis that these may be captured C-class asteroids from a postaccretional dispersion. The proposed testing technique could be applied to an Earth-based survey of asteroids or a space-probe study of an asteroid's surface characteristic in order to identify a potential water source.

Description: The Coast Ranges of the Cascadia margin are overriding the subducted Juan de Fuca/Gorda plate. We investigate the extent to which the latitudinal change in attributes related to the subduction process. These attributes include the varibale age of the subducted slab that underlies the Coast Ranges and average vertical crustal velocities of the western margin of the Coast Rnages for two markedly different time periods, the last 45 years and the last 100 kyr. These vertical crustal velocities are computed from the resurveying of highway bech marks and from the present elevation of shore platforms that have been uplifted in the late Quaternary, respectively. Topogarphy of the Coast Ranges is in part a function of the age and bouyancy of the underlying subducted plate. This is evident in the fact that the two highest topographic elements of the Coast Rnages, the Klamath Mountains and the Olympic Mountains, are underlain by youngest subducted oceanic crust. The subducted Blanco Fracture Zone in southernmost Oregon is responsible for an age discontinuity of subducted crust under the Klamath Mountains. The norhtern terminus of hte topographically higher Klamaths is offset to the north relative to the position of the underlying Blanco Fracture Zone, teh offset being in the direction of migration of the farcture zone, as dictated by relative plate motions. Vertical crustal velocities at the coast, derived from becnh mark surveys, are as much as an order of magnitude greater than vertical crustal velocities derived from uplifted shore platforms. This uplift rate discrepancy indicates that strain is accumulating on the plate margin, to be released during the next interplate earthquake. In a latitudinal sense, average Coast Rnage topography is relatively high where bench mark-derived, short-term vertical crustal velocities are highest. Becuase the shore platform vertical crustal velocities reflect longer-term, premanent uplift, we infer that a small percentage of the interseismic strain that accumulates as rapid short-term uplift is not recovered by subduction earthquakes but rather contributes to rock uplift of the Coast Ranges. The conjecture that permanent rock uplift is related to interseismic uplift is consistent with the observation that those segments of the subduction zone subject to greater interseismic uplift rates are at approximately the same latitudes as those segments of the Coast Ranges that have higher magnitudes of rock uplift over the long term.

Description: A comprehensive review is presented of the mathematical models used to represent magnetic fields in the Earth's magnetosphere, of the way existing data-based models use these methods and of the associated problems and concepts. The magnetic field has five main components: the internal field, the magnetopause, the ring current, the tail and Birkeland currents. Methods of representing separately each of these are discussed, as is the deformation of magnetic fields; Appendix B traces the connection between deformations and the Cauchy integral. A summary section lists the uses of data-based models and their likely future evolution, and Appendix A supplements the text with a set of problems.

Description: We present a map of the coseimic displacement field resulting from the Landers, California, June 28, 1992, earthquake derived using data acquired from an orbiting high-resolution radar system. We achieve results more accurate than previous space studies and similar in accuracy to those obtained by conventional field survey techniques. Data from the ERS 1 synthetic aperture radar instrument acquired in April, July, and August 1992 are used to generate a high-resolution, wide area map of the displacements. The data represent the motion in the direction of the radar line of sight to centimeter level precision of each 30-m resolution element in a 113 km by 90 km image. Our coseismic displacement contour map gives a lobed pattern consistent with theoretical models of the displacement field from the earthquake. Fine structure observed as displacement tiling in regions several kilometers from the fault appears to be the result of local surface fracturing. Comparison of these data with Global Positioning System and electronic distance measurement survey data yield a correlation of 0.96; thus the radar measurements are a means to extend the point measurements acquired by traditional techniques to an area map format. The technique we use is (1) more automatic, (2) more precise, and (3) better validated than previous similar applications of differential radar interferometry. Since we require only remotely sensed satellite data with no additioanl requirements for ancillary information. the technique is well suited for global seismic monitoring and analysis.

Description: The infrared transmission spectra and photochemical behavior of various organic compounds isolated in solid N2 ices, appropriate for applications to Triton ad Pluto, are presented. It is shown that excess absorption in the surface spectra of Triton and Pluto, i.e., absorption not explained by present models incorporating molecules already identified on these bodies (N2, CH4, CO, and CO2), that starts near 4450/cm (2.25 microns) and extends to lower frequencies, may be due to alkanes (C(n)H(2n+2)) and related molecules frozen in the nitrogen. Branched and linear alkanes may be responsible. Experiments in which the photochemstry of N2: CH4 and N2: CH4: CO ices was explored demonsrtrate that the surface ices of Triton and Pluto may contain a wide variety of additional species containing H, C, O, and N. Of these, the reactive molecule diazomethane, CH2N2, is particularly important since it may be largely responsible for the synthesis of larger alkanes from CH4 and other small alkanes. Diazomethane would also be expected to drive chemical reactions involving organics in the surface ices of Triton and Pluto toward saturation, i.e., to reduce multiple CC bonds. The positions and intrinsic strengths (A values) of many of the infrared absorption bands of N2 matrix-isolated molecules of relevance to Triton and Pluto have also been determined. These can be used to aid in their search and to place constraints on their abundances.

Description: High-resolution (0.01/cm) absorption spectra of lean mixtures of CH4 in dry air were recorded with the McMath-Pierce Fourier transform spectrometer (FTS) of the National Solar Observatory on Kitt Peak at various temperatures between 24 and -61 C. The spectra have been analyzed to determine the values at room temperature of pressure-broadened widths and pressure-induced shifts of more than 740 transitions. The temperature dependence of air-broadened widths and pressure-induced shifts was deduced for approx. 370 transitions in the nu(sub 1) + nu(sub 4), nu(sub 3) + nu(sub 4), and nu(sub 2) + nu(sub 3) bands of (12)CH4 located between 4118 and 4615/cm. These results were obtained by analyzing a total of 29 spectra simultaneously using a multi-spectral non-linear least-squares fitting technique. This new technique allowed the determination of correlated spectral line parameters (e.g. intensity and broadening coefficient) better than the procedure of averaging values obtained by fitting the spectra individually. This method also provided a direct determination of the uncertainties in the retrieved parameters due to random errors. For each band analysed in this study the dependence of the various spectral line parameters upon the tetrahedral symmetry species and the rotational quantum numbers of the transitions is also presented.

Description: This report presents the most recent spherical harmonic topography model of Venus developed at Jet Propulsion Laboratory. It was produced by a spherical harmonic analysis of the most complete set of Magellan altimetry data, augmented by Pioneer Venus and Venera data. The harmonic coefficients of the topography were computed to degree and order 360. Compared to previous topography models, this one has the highest correlation with the gravity field of Venus.

Description: The 500-Myr average crater retention age for Venus has raised questions about the present-day level of tectonic activity. In this study we examine the relationship between the gravity and topography of four large volcanic swells, Beta, Atla, Bell, and Western Eistla Regiones, for clues about their stage evolution. The Magellan line-of-sight gravity data are inverted using a point mass model of the anomalous mass to solve for the local vertical gravity field. Spectral admittance calculated from both the local gravity inversions and a spherical harmonic model is compared to three models of compensation: local compensation, a 'flexural' model with local and regional compensation of surface and subsurface loads, and a 'hotspot' model of compensation that includes top loading by volcanoes and subsurface loading due to a deep, low density mass anomaly. The coherence is also calculated in each region, but yields an elastic thickness estimate only at Bell Regio. In all models, the long wavelengths are compensated locally. Our results may indicate a relatively old, possibly inactive plume.

Description: New line parameters for two heavy odd nitrogen molecules HNO3 in the nu(sub 5)/2nu(sub 9) region, and ClONO2 in the nu(sub 4) region are incorporated in the analysis of high resolution i.r. atmospheric spectra. The line parameters are tested and renormalized vs laboratory spectra, and then applied to retrievals from balloon-borne and ground-based solar absorption spectra.

Description: Spectral absorption-coefficients (cross-sections) kappa(sub nu) (/cm/atm) have been measured in the 7.62, 8.97, and 12.3 micrometer bands of HCFC-22 (CHClF2) and the 10.6 micrometer bands of SF6 employing a high-resolution Fourier-transform spectrometer. Temperature and total pressure have been varied to simulate conditions corresponding to tropospheric and stratospheric layers in the atmosphere. The kappa(sub nu) are compared with values measured by us previously using a tunable diode laser spectrometer and with the appropriate entries in HITRAN and GEISA, two of the databases known to the atmospheric scientist. The measured absolute intensities of the bands are compared with previously published values.

Description: The spectroscopic identification for the HNO3 3 nu(sub 9) - nu(sub 9) band Q branch at 830.4/cm is reported based on 0.01/cm resolution solar occultation spectra of the lower stratosphere recorded by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer and a recent analysis of this band. Least-squares fits to 0.0025/cm resolution laboratory spectra in the Q branch region indicate an integrated intensity of 0.529 x 10(exp -18)/cm/mol/sq cm at 296 K for this weak band. Stratospheric HNO3 retrievals derived from the ATMOS data are consistent with this value within its estimated uncertainty of about +/- 30%. A set of spectroscopic line parameters suitable for atmospheric studies has been generated.

Description: About 200 i.r. solar spectra recorded at 0.01/cm resolution on 71 days between November 1991 and July 1993 at the Network for the Detection of Stratospheric Change (NDSC) station at Mauna Loa, Hawaii (latitude 19.53 deg N, longitude 155.58 deg W, elevation 3.459 km) have been analyzed with a nonlinear least-squares spectral fitting technique to study temporal variations in the total column of atmospheric ethane (C2H6) above the site. The results were derived from the analysis of the unresolved nu(sub 7) band (P)Q(sub 3) subbranch at 2976.8/cm. A distinct seasonal cycle is observed with a factor of 2 variation, a maximum total column of 1.16 x 10(exp 16) mol/sq cm at the end of winter, and a minimum total column of 0.53 x 10(exp 16) mol/sq cm at the end of summer. Our measurements are compared with previous observations and model predictions.

Description: Quantitative measurements of the wavelength dependence of aerosol extinction in the 750-3400/cm spectral region have been derived from 0.01/cm resolution stratospheric solar occultation spectra recorded by the ATMOS (Atmospheric Trace Molecule Spectroscopy) Fourier transform spectrometer about 9 1/2 months after the Mt Pinatubo volcanic eruption. Strong, broad aerosol features have been identified near 900, 1060, 1190, 1720, and 2900/cm below a tangent height of approximately 30 km. Aerosol extinction measurements derived from approximately 0.05/cm wide microwindows nearly free of telluric line absorption in the ATMOS spectra are compared with transmission calculations derived from aerosol size distribution profiles retrieved from correlative SAGE (Stratospheric Aerosol and Gas Experiment) II visible and near i.r. extinction measurements, seasonal and zonally averaged H2SO4 aerosol weight percentage profiles, and published sulfuric acid optical constants derived from room temperature laboratory measurements. The calculated shapes and positions of the aerosol features are generally consistent with the observations, thereby confirming that the aerosols are predominantly concentrated H2SO4-H2O droplets, but there are significant differences between the measured and calculated wavelength dependences of the aerosol extinction. We attribute these differences as primarily the result of errors in the calculated low temperature H2SO4-H2O optical constants. Errors in both the published room temperature optical constants and the limitations of the Lorentz-Lorenz relation are likely to be important.

Description: We present vertical column abundances of H2O, N2O, HNO3, NO2, O3, HF, HCl, and ClNO3, determined from solar absorption spectra measured by the JPL MkIV interferometer from the NASA DC-8 aircraft. These observations, taken in 1987 and 1992, covered latitudes ranging from 85 deg S to 85 deg N. Although most gases display latitude symmetry, large asymmetries in H2O, HNO3, and O3 are apparent, which can be ascribed to processes enhanced by the colder Antarctic winter temperatures.

Description: Simultaneous in situ measurements of NO2, NO, O3, ClO, pressure and temperature have been made for the first time, presenting a unique opportunity to test our current understanding of the photochemistry of the lower stratospere. Data were collected from several flights of the ER-2 aircraft at mid-latitudes in May 1993 during NASA's Stratospheric Photochemistry, Aerosols and Dynamics Expedition (SPADE). The daytime ratio of NO2/NO remains fairly constant at 19 km with a typical value of 0.68 and standard deviation of +/- 17. The ratio observations are compared with simple steady-state calculations based on laboratory-measured reaction rates and modeled NO2 photolysis rates. At each measurement point the daytime NO2/NO with its measurements uncertainty overlap the results of steady-state caculations and associated uncertainty. Possible sources of error are examined in both model and measurements. It is shown that more accurate laboratory determinations of the NO + 03 reaction rate and of the NO2 cross-sections in the 200-220 K temperature range characteristic of the lower stratosphere would allow for a more robust test of our knowledge of NO(X) phtochemistry by reducing significant sources if uncertainties in the interpretation of statospheric measurements.

Description: The ability to predict short-term variations in the Earth's rotation has gained importance in recent years owing to more precise spacecraft tracking requirements. Universal time (UT1), that component of the Earth's orientation corresponding to the rotation angle, can be measured by number of high-precision space geodetic techniques. A Kalman filter developed at the Jet Propulsion Laboratory (JPL) optimally combines these different data sets and generates a smoothed times series and a set of predictions for UT1, as well as for additional Earth orientation components. These UT1 predictions utilize an empirically derived random walk stochastic model for the length of the day (LOD) and require frequent and up-to-date measurements of either UT1 or LOD to keep errors from quickly accumulating. Recent studies have shown that LOD variations are correlated with changes in the Earth's axial atmospheric angular momentum (AAM) over timescales of several years down to as little as 8 days. AAM estimates and forecasts out to 10 days are routinely available from meteorological analysis centers; these data can supplement geodetic measurements to improve the short-term prediction of LOD and have therefore been incorporated as independent data types in the JPL Kalman filter. We find that AAM and, to a lesser extent, AAM forecast data are extremely helpful in generating accurate near-real-time estimates of UT1 and LOD and in improving short-term predictions of these quantities out to about 10 days.

Description: Recent mapping studies west of Elysium Mons, Mars, have pinpointed subice features that suggest the existence of a frozen paleolake in Utopia Planitia as recently as 1.8 billion years ago. The subice features are interpreted to be hyaloclastic ridges and hills, table moutains, associated joekulhalaup deposits, and fluvial channels. Photoclinometric studies of these features and of a basal scarp around the northwest flank of Elysium Mons interpreted to have been an ice-sheet boundary indicate that the maximum thickness of ice within the basin may have been about 180 m. This thickness of ice during a relatively late stage of Martian geologic history would have important implications concerning the atmospheric, the climatic, and possibly the exobiologic history of the planet.

Description: Disk-shaped current distributions are useful tools in modeling the magnetospheres of Earth and other planets and have also been adapted for modeling the geotail. Such models usually start with an axisymmetric vector potential but may be modified to account for observed asymmetries, variable thickness, and warping in response to an inclined orientation of the planetary dipole axis. Models of this type until now either have lacked the ability to simulate a sharp inner edge of the current and to control accurately its falloff with distance or did not allow a simple analytical representation. Here existing methods will be reviewed, after which a new class of models which overcomes the above deficiencies and also allows the modeling of current disks of finite thickness flanked by current-free regions will be presented.

Description: An improved model of Earth's gravitational field, Goddard Earth Model T-3 (GEM-T3), has been developed from a combination of satellite tracking, satellite altimeter, and surface gravimetric data. GEM-T3 provides a significant improvement in the modeling of the gravity field at half wavelengths of 400 km and longer. This model, complete to degree and order 50, yields more accurate satellite orbits and an improved geoid representation than previous Goddard Earth Models. GEM-T3 uses altimeter data from GEOS 3 (1975-1976), Seasat (1978) and Geosat (1986-1987). Tracking information used in the solution includes more than 1300 arcs of data encompassing 31 different satellites. The recovery of the long-wavelength components of the solution relies mostly on highly precise satellite laser ranging (SLR) data, but also includes Tracking Network (TRANET) Doppler, optical, and satellite-to-satellite tracking acquired between the ATS 6 and GEOS 3 satellites. The main advances over GEM-T2 (beyond the inclusion of altimeter and surface gravity information which is essential for the resolution of the shorter wavelength geoid) are some improved tracking data analysis approaches and additional SLR data. Although the use of altimeter data has greatly enhanced the modeling of the ocean geoid between 65 deg N and 60 deg S latitudes in GEM-T3, the lack of accurate detailed surface gravimetry leaves poor geoid resolution over many continental regions of great tectonic interest (e.g., Himalayas, Andes). Estimates of polar motion, tracking station coordinates, and long-wavelength ocean tidal terms were also made (accounting for 6330 parameters). GEM-T3 has undergone error calibration using a technique based on subset solutions to produce reliable error estimates. The calibration is based on the condition that the expected mean square deviation of a subset gravity solution from the full set values is predicted by the solutions' error covariances. Data weights are iteratively adjusted until this condition for the error calibration is satisfied. In addition, gravity field tests were performed on strong satellite data sets withheld from the solution (thereby ensuring their independence). In these tests, the performance of the subset models on the withheld observations is compared to error projections based on their calibrated error covariances. These results demonstrate that orbit accuracy projections are reliable for new satellites which were not included in GEM-T3.

Description: The ion and electron momentum equations, along with Ampere's law, are solved for the ion and electron drift velocities and the electric field in the subsolar Venus ionosphere, assuming a partially ionized gas and a single ion species having the ion mean mass. All collision terms among the ions, electrons and neutral particles are retained in the equations. A general expression for the evolution of the magnetic field is derived and compared with earlier expressions. Subsolar region data in the altitude range 150-300 km from the Pioneer Venus Orbiter are used to calculate altitude profiles of the components of the current due to the electric field, gradients of pressure, and gravity. Altitude profiles of the ion and electron velocities as well as the electric field, electrodynamic heating, and the energy density are determined. Only orbits having a complete set of measured plasma temperatures and densities, neutral densities, and magnetic field were considered for analysis; the results are shown only for orbit 202. The vertical velocity at altitudes above 220 km is upgoing for orbit 202. This result is consistent with observations of molecular ions at high altitudes and of plasma flow to the nightside, both of which require upward velocity of ions from the dayside ionosphere. Above about 230 km the momentum equations are extremely sensitive to the altitude profiles of density, temperature, and magnetic field.

Description: We present predictions of the signatures of magnetosheath particle precipitation (in the regions classified as open low-latitude boundary layer, cusp, mantle and polar cap) for periods when the interplanetary magnetic field has a southward component. These are made using the 'pulsating cusp' model of the effects of time-varying magnetic reconnection at the dayside magnetopause. Predictions are made for both low-altitude satellites in the topside ionosphere and for midaltitude spacecraft in the magnetosphere. Low-altitude cusp signatures, which show a continuous ion dispersion signature, reveal 'quasi-steady reconnection' (one limit of the pulsating cusp model), which persists for a period of at least 10 min. We estimate that 'quasi-steady' in this context corresponds to fluctuations in the reconnection rate of a factor of 2 or less. The other limit of the pulsating cusp model explains the instantaneous jumps in the precipitating ion spectrum that have been observed at low altitudes. Such jumps are produced by isolated pulses of reconnection: that is, they are separated by intervals when the reconnection rate is zero. These also generate convecting patches on the magnetopause in which the field lines thread the boundary via a rotational discontinuity separated by more extensive regions of tangential discontinuity. Predictions of the corresponding ion precipitation signatures seen by midaltitude spacecraft are presented. We resolve the apparent contradiction between estimates of the width of the injection region from midaltitude data and the concept of continuous entry of solar wind plasma along open field lines. In addition, we reevaluate the use of pitch angle-energy dispersion to estimate the injection distance.

Description: We report here analyses of olivines and pyroxenes, and petrofabrics of 27 chondritic interplanetary dust particles (IDPs), comparing those from anhydrous and hydrous types. Approximately 40% of the hydrous particles contain diopside, a probable indicator of parent body thermal metamorphism, while this mineral is rarely present in the anhydrous particles. Based on this evidence, we find that hydrous and anhydrous IDPs are, in general, not directly related, and we conclude that olivine and pyroxene major-element compositions can be used to help discriminate between IDPs that are (1) predominantly nebular condensates, and lately resided in anhydrous or icy (no liquids) primitive parent bodies, and (2) those originating from more geochemically active parent bodies (probably hydrous and anhydrous asteroids).

Description: Simultaneous profiles of aerosol backscatter ratio were measured over Lauder, New Zealand (45 deg S, 170 deg E) on the night of November 24, 1992. Instrumentation comprised two complementary lidar systems and a backscattersonde, to give measurements at wavelengths 351, 490, 532, and 940 nm. The data from the lidars and the backscattersonde were self-consistent, enabling the wavelength dependence of aerosol backscatter to be determined as a function of altitude. This wavelength-dependence is a useful parameter in radiative transfer calculations. In the stratosphere, the average wavelength exponent between 351 and 940 nm was -1.23 +/- 0.1, which was in good agreement with values derived from measured physical properties of aerosols.

Description: High-energy electrons have been measured systematically in a low-altitude (520 x 675 km), nearly polar (inclination = 82 deg) orbit by sensitive instruments onboard the Solar, Anomalous, and Magnetospheric Particle Explorer (SAMPEX). Count rate channels with electron energy thresholds ranging from 0.4 MeV to 3.5 MeV in three different instruments have been used to examine relativistic electron variations as a function of L-shell parameter and time. A long run of essentially continuous data (July 1992 - July 1993) shows substantial acceleration of energetic electrons throughout much of the magnetosphere on rapid time scales. This acceleration appears to be due to solar wind velocity enhancements and is surprisingly large in that the radiation belt 'slot' region often is filled temporarily and electron fluxes are strongly enhanced even at very low L-values (L aprroximately 2). A superposed epoch analysis shows that electron fluxes rise rapidly for 2.5 is approximately less than L is approximately less than 5. These increases occur on a time scale of order 1-2 days and are most abrupt for L-values near 3. The temporal decay rate of the fluxes is dependent on energy and L-value and may be described by J = Ke-t/to with t(sub o) approximately equals 5-10 days. Thus, these results suggest that the Earth's magnetosphere is a cosmic electron accelerator of substantial strength and efficiency.

Description: This paper reports on a magnetic field phenomenon, hereafter referred to as null fields, which were discovered during the inbound pass of the recent flyby of Jupiter by the Ulysses spacecraft. These null fields which were observed in the outer dayside magnetosphere are characterised by brief but sharp decreases of the field magnitude to values less than 1 nT. The nulls are distinguished from the current sheet signatures characteristic of the middle magnetosphere by the fact that the field does not reverse across the event. A field configuration is suggested that accounts for the observed features of the events.

Description: Rocket data have been used to evaluate the characteristics of precipitating relativistic electrons and their effects on the electrodynamic structure of the middle atmosphere. These data were obtained at Poker Flat, Alaska, on May 13 and 14, 1990, during a midday, highly relativistic electron (HRE) precipitation event. Solid state detectors were used to measure the electron fluxes and their energy spectra. An X ray scintillator was included on each flight to measure bremsstrahlung X rays produced by energetic electrons impacting on the upper atmosphere. However, these were found the be of negligible importance for this particular event. The energy deposition by the electrons has been determined from the flux measurements and compared with in situ measurements of the atmospheric electrical response. The electrodynamic measurements were obtained by the same rockets and additionally on May 13, with an accompanying rocket. The impact flux was highly irregular, containing short-lived bursts of relativistic electrons, mainly with energies below 0.5 MeV and with fluxes most enhanced between pitch angles of 0 deg - 20 deg. Although the geostationary counterpart of this measured event was considered to be of relatively low intensity and hardness, energy deposition peaked near 75 km with fluxes approaching an ion pair production rate in excess of 100/cu cm s. This exceeds peak fluxes in relativistic electron precipitation (REP) events as observed by us in numerous rocket soundings since 1976. Conductivity measurements from a blunt probe showed that negative electrical conductivities exceeded positive conductivities down to 50 km or lower, consistent with steady ionization by precipitating electrons above 1 MeV. These findings imply that the electrons from the outer radiation zone can modulate the electrical properties of the middle atmosphere to altitudes below 50 km. During the decline and activity minimum of the current solar cycle, we anticipate the occurence of similar events but with fluxes 1-2 orders of magnitude above that reported here, based on studies of earlier solar cycles (e.g., Baker et al., 1993).

Description: Nitrous oxide (N2O) measured on board the ER-2 aircraft during the Airborne Arctic Stratospheric Expedition 2 (AASE 2) has been used to monitor descent of air inside the Arctic vortex between October 1991 and March 1992. Monthly mean N2O fields are calculated from the flight data and then compared with mean fields calculated from the high-resolution Geophysical Fluid Dynamics Laboratory general circulation model SKYHI in order to evaluate the model's simulation of the polar vortex. From late fall through winter the model vortex evolves in much the same way as the 1991-1992 vortex, with N2O gradients at the edge becoming progressively steeper. The October to March trends in N2O profiles inside the vortex are used to verify daily net heating rates in the vortex that were computed from clear sky radiative heating rates and National Meteorological Center temperature observations. The computed heating rates successfully estimate the descent of vortex air from December through February but suggest that before December, air at high latitudes may not be isolated from the midlatitudes. SKYHI heating rates are in good agreement with the computed rates but tend to be slightly higher (i.e., less cooling) due to meteorological differences between SKYHI and the 1991-1992 winter. Three ER-2 flights measured N2O just north of the subtropical jet. These low-midlatitude profiles show only slight differences from the high-midlatitude profiles (45 deg - 60 deg N), indicating strong meridional mixing in the midlatitude 'surf zone.' Mean midwinter N2O profiles inside and outside the vortex calculated from AASE 2 data are shown to be nearly identical to 1989 AASE profiles, pointing to the N2O/potential temperature relationship as an excellent marker for vortex air.

Description: The combined Nimbus 7 solar backscattered ultraviolet (SBUV) and NOAA 11 SBUV/2 ozone data, covering a period of more than a solar cycle (about 15 years), are used to study the UV response of ozone in the stratosphere. The study shows that about 2% change in total column ozone and about 5-7% change in ozone mixing ratio in the upper stratosphere (0.7 to 2 hPa) may be attributed to the change in the solar UV flux over a solar cycle. In the upper stratosphere, where photochemical processes are expected to play a major role, the measured solar cycle variation of ozone is significantly larger than inferred either from the photochemical models or from the ozone response to the 27-day solar UV modulation. For example, the observed solar cycle related change in ozone mixing ratio at 2 hPa is about 1% for 1% change in the solar UV flux near 200 nm. The inferred change in ozone from either the photochemical models or from the 27-day ozone-UV response is about a factor of 2-3 lower than this value.

Description: Vertical profiles of N2O and NO(y) taken by the ER-2 outside the vortex are used to construct average vertical profiles of F(NO(y)) = NO(y)/(A-N2O), where A is the tropospheric content of N2O three years prior to the measurements. The southern hemisphere had less nitrous oxide in the range 400 less than Theta less than 470 K, by up to 25% relative to the northern hemisphere. F(NO(y)) is the ratio of NOy produced to N2O lost in a stratospheric air mass since entry from the troposphere. The profiles of F(NO(y)) have the following characteristics: (1) Relative to 1991-1992, a year without denitrification inside or outside the vortex, the northern hemisphere in 1988-1989 showed denitrification outside the vortex ranging up to 25% and averaging 17% above Theta = 425 K. (2) Relative to the northern hemisphere in 1991-1992, the southern hemisphere in 1987 showed denitrification outside the vortex ranging up to 32% and averaging 20% above Theta = 400 K. (3) Below Theta = 400 K the southern hemisphere showed enhancements of F(NO(y)) relative to the northern hemisphere in 1991-1992 ranging up to 200% at Theta = 375 K, outside the vortex. Corresponding profiles of residual water, R(H2O) = H2O - 2(1.6 - CH4), are considered and shown to be consistent with those of F(NO(y)) in the sense that they show deficits outside the Antarctic vortex, which was both dehydrated and denitrified, but not outside the 1988-1989 Arctic vortex, which was denitrified but not dehydrated. R(H2O) is the water content of stratospheric air with the contribution from methane oxidation subtracted. Comparison of F(NO(y)) and R(H2O) below 400 K outside the Antarctic vortex leads to the suggetion that dehydration in the Antarctic vortex occurs by the sedimentation of ice crystals large enough to fall out of the stratosphere, whereas denitrification occurs mainly on mixed nitric acid-water crystals which evaporate below the base of the vortex at Theta = 400 K but above the tropopause.

Description: During January and August 1985, the scanning radiometers of the Earth Radiation Budget Experiment(ERBE) aboard the Earth Radiation Budget Satellite (ERBS) and the NOAA-9 satellite were operated in along-track scanning modes. Along-track scanning permits the study of many measurement problems. It provides the data for developing a limb-darkening model for a single site over a short period of time and also permits the indentification of the scene from data taken at smaller nadir angles. The earth-emitted radiation measured by the scanners has been analyzed to produce limb-darkening models for a variety of scene types. Limb-darkening models relate the radiance in any given direction to the radiant flux. The scene types were computed using measurements within 10 deg of zenith. The models have values near zenith of 1.02-1.09. The typical zenith values of the model are 1.06 for both day and night for ERBS, and for NOAA-9, 1.06 for day and 1.05 for night. Mean models are formed for the ERBS and NOAA-9 results and are found to differ less than 1%, the ERBS results being the higher. The models vary about 1% with latitude near zenith and agree with earlier models that were used to analyze ERBE data typically to 2%.

Description: Geochemical profiles of surface units, impact, and volcanic features are studied in detail to determine the underlying structure in an area of extensive mare/highland interface, Sinus Amoris. This study region includes and surrounds the northeastern embayment of Mare Tranquillitatis. The concentrations of two major rock-forming elements (Mg and Al), which were derived from the Apollo 15 orbital geochemical measurements, were used in this study. Mapped units and deposits associated with craters in the northwestern part of the region tend to have correlated low Mg and Al concentrations, indicating the presence of Potassium (K)-Rare Earth Elements (REE)-Phosphorus (P) (KREEP)-enriched basalt. Found along the northeastern rim of Tranquillitatis were areas with correlated high Mg and Al concentration, indicating the presence of troctolite. Distinctive west/east and north/south trends were observed in the concentrations of Mg and Al, and, by implication, in the distribution of major rock components on the surface. Evidence for a systematic geochemical transition in highland or basin-forming units may be observed here in the form of distinctive differences in chemistry in otherwise similar units in the western and eastern portions of the study region.

Description: One of the oldest mysteries in geomagnetism is the linkage between solar and geomagnetic activity. In investigating the causes of geomagnetic storms occurring during solar maximum, the following topics are discussed: solar phenomena; types of solar wind; magnetic reconnection and magnetic storms; an interplanetary example; and future space physics missions.

Description: Atmospheric mass loading produces a primarily vertical displacement of the Earth's crust. This displacement is correlated with surface pressure and is large enough to be detected by very long baseline interferometry (VLBI) measurements. Using the measured surface pressure at VLBI stations, we have estimated the atmospheric loading term for each station location directly from VLBI data acquired from 1979 to 1992. Our estimates of the vertical sensitivity to change in pressure range from 0 to -0.6 mm/mbar depending on the station. These estimates agree with inverted barometer model calculations (Manabe et al., 1991; vanDam and Herring, 1994) of the vertical displacement sensitivity computed by convolving actual pressure distributions with loading Green's functions. The pressure sensitivity tends to be smaller for stations near the coast, which is consistent with the inverted barometer hypothesis. Applying this estimated pressure loading correction in standard VLBI geodetic analysis improves the repeatability of estimated lengths of 25 out of 37 baselines that were measured at least 50 times. In a root-sum-square (rss) sense, the improvement generally increases with baseline length at a rate of about 0.3 to 0.6 ppb depending on whether the baseline stations are close to the coast. For the 5998-km baseline from Westford, Massachusetts, to Wettzell, Germany, the rss improvement is about 3.6 mm out of 11.0 mm. The average rss reduction of the vertical scatter for inland stations ranges from 2.7 to 5.4 mm.

Description: A theoretical model is used to describe and investigate the effects of simultaneous crystallization, radiation loss, and entrainment of cooler material on the temperature of a well-mixed core of an active aa lava flow. Entrainment of crust, levee debris, and base material into the interior of active flows has been observed, but the degree of assimilation and the thermal consequences are difficult to quantify. The rate of entrainment can be constrained by supplementing the theoretical model with information on the crystallization along the path of the flow and estimation of the radiative loss from the flow interior. Application of the model is demonstrated with the 1984 Mauna Loa flow, which was erupted about 30 C undercooled. Without any entrainment of cooler material, the high crystallization rates would have driven temperatures in the core wall above temperatures measured by thermocouple and estimated from glass geothermometry. One plausible scenario for this flow, which agrees with available temperature and crystallinity measurements, has a high initial rate of entrainment during the first 8 hours of travel (a mass ratio of entrained material to fluid core of about 15% if the average temperature of the entrained material was 600 C), which counterbalances the latent heat from approximately 40% crystallization. In this scenario, the model suggests an additional 5% crystallization and a 5% entrainment mass ratio over the subsequent 16-hour period. Measurements of crystallization, radiative losses, and entrainment factors are necessary for understanding the detailed thermal histories of active lava flows.

Description: The data from the infrared telescope (IRT), which was flown on space shuttle Challenger Spacelab 2 mission (July 1985), were originally reported by Koch et al. (1987) as originating from near orbital emissions, primarily H2O. In this study, analysis of this data was extended to determine the collisional cross sections for the excitation of the low lying vibrational levels of H2O, present in the orbiter cloud, by atmospheric O(3P). The evaluation of the contribution to the measured signal from solar excitation and ram O excitation of outgassing H2O permits the determination of the H2O column density and the excitation cross section of the (101) level at an O(3P) velocity of approximately 7.75 km/s. Contributions to the radiation in the 1.7-3.0 micron band by transitions from the (100), (001), and multiquantum excited levels are discussed. The findings of the study are (1) the IRT data for the 4.5-9.5 micron and the nighttime data for the 1.7-3.0 micron sensors are consistent with being explained by collision excitation of H2O by O(3P), (2) diurnal variations of 4.5-9.5 micron intensities follow the model predicted O density for a full orbit, (3) daytime increases in the H2O cloud density were not evident, (4) the cross sections for the collisional excitation process are derived and compared to values computated by Johnson (1986) and Redmon et al. (1986), (5) theoretical investigation suggests greater than 60% of the radiation from H2O is a result of multiphoton emission resulting from collisional multiquanta excitation, and (6) the large daytime increase in the 1.7-3.0 micron intensity data suggests that O(+) may likely be instrumental in producing excited H2O(+) through charge exchange.

Description: The aerobraking orbital activities of Magelland during the gravity mapping of Venus are discussed. The goal of aerobraking was to circularize Magellan's orbit. By aerobraking the spacecraft into a nearly circula orbit, the Magellan team was able to provide scientists with a different data set to deepen their understanding of what is going on beneath Venus' surface. Before undertaking its gravity-mapping mission, Magellan completed three cycles of radar mapping. This repeated coverage allowed the spacecraft to see some of Venus' geologic features from different viewing angles. Various aspects of the mission are discussed, and maps of Venus are presented.

Description: The devolatilization of calcium sulfate, which is present in the target rock of the Chicxulub, Mexico impact structure, and dispersal in the stratosphere of the resultant sulfuric acid aerosol have been suggested as a possible mechanism for the Cretaceous-Tertiary extinctions. We measured the amount of SO2 produced from two shock-induced devolatilization reactions of calcium sulfate up to 42 GPa in the laboratory. We found both to proceed to a much lower extent than calculated by equilibrium thermodynamic calculations. Reaction products are found to be approx. 10(exp -2) times those calculated for equilibrium. Upon modeling the quantity of sulfur oxides degassed into the atmosphere from shock devolatilization of CaSO4 in the Chicxulub lithographic section, the resulting 9 x 10(exp 16) to 6 x 10(exp 17) g (in sulfur mass) is lower by a factor of 10-100 than previous upper limit estimates, the related environmental stress arising from the resultant global cooling and fallout of acid rain is insufficient to explain the widespread K-T extinctions.

Description: Energetic atomic (O(+1) and N(+1)) and molecular (O2(+1), NO(+1), and N2(+1)) ions of ionospheric origin were observed in Earth's magnetotail at X approximately -146 R(sub E) during two plasma sheet sunward/tailward flow-reversal events measured by instruments on the GEOTAIL spacecraft. These events were associated with concurrent ground-measured geomagnetic disturbance intensification at auroral-and mid-latitudes (Kp = 7(-)). Energetic ions in the sunward-component and tailward flows were from both the solar wind and ionosphere. Plasma and energetic ions participated in the flows. During tailward flow, ionospheric origin ion abundance ratios at approximately 200-900 km/s in the rest frame were N(+1)/O(+1) = approximately 25-30% and ((O2(+1), NO(+1), and N2(+1))/O(+1) = approximately 1-2%. We argue that tailward flow most likely initiated approximately 80-100 R(sub E) tailward of Earth and molecular ions were in the plasma sheet prior to geomagnetic intensification onset.

Description: Wave forms of BEN (Broadband Electrostatic Noise) in the geomagnetic tail were first detected by the Wave Form Capture reciever on the GEOTAIL spacecraft. The results show that most of the BEN in the plasma sheet boundary layer (PSBL) are not continuous broadband noise but are composed of a series of solitary pulses having a special form which we term 'Electrostatic Solitary Waves (ESW)'. A nonlinear BGK potential model is proposed as the generation mechanism for the ESW based upon a simple particle simulation which considers the highly nonlinear evolution of the electron beam instability. The wave forms produced by this simulation are very similar to those observed by GEOTAIL and suggest that the nonlinear dynamics of the electron beam play an essential role in the generation of ESW.

Description: We report on the observations of a number of quasi-dc electric field events associated with large-scale atmospheric weather formations. The observations were made by the electric field experiment onboard the San Marco D satellite, operational in an equatorial orbit from May to December 1988. Several theoretical studies suggest that electric fields generated by thunderstorms are present at high altitudes in the ionosphere. In spite of such favorable predictions, weather-related events are not often observed since they are relatively weak. We shall report here on a set of likely E field candidates for atmospheric-ionospheric causality, these being observed over the Indonesian Basin, northern South America, and the west coast of Africa; all known sites of atmospheric activity. As we shall demonstrate, individual events often be traced to specific active weather features. For example, a number of events were associated with spacecraft passages near Hurricane Joan in mid-October 1988. As a statistical set, the events appear to coincide with the most active regions of atmospheric weather.

Description: A radiation model, together with National Meteorological Center temperature observations, was used to compute daily net heating rates in the northern hemisphere (NH) for the Arctic late fall and winter periods of both 1988-1989 and 1991-1992 and in the southern hemisphere (SH) for the Antarctic fall and winters of 1987 and 1992. The heating rates were interpolated to potential temperature (theta) surfaces between 400 K and 2000 K and averaged within the polar vortex, the boundary of which was determined by the maximum gradient in potential vorticity. The averaged heating rates were used in a one-dimensional vortex interior descent model to compute the change in potential temperature with time of air parcels initialized at various theta values, as well as to compute the descent in log pressure coodinates. In the NH vortex, air parcels which were initialized at 18 km on November 1, descended about 6 km by March 21, while air initially at 25 km descended 9 km in the same time period. this represents an average descent rate in the lower stratosphere of 1.3 to 2 km per month. Air initialized at 50 km descended 27 km between November 1 and March 21. In the SH vortex, parcels initialized at 18 km on March 1, descended 3 km, while air at 25 km descended 5-7 km by the end of October. This is equivalent to an average descent in the lower stratosphere of 0.4 to 0.9 km per month during this 8-month period. Air initialized at 52 km descended 26-29 km between March 1 and October 31. In both the NH and the SH, computed descent rates increased markedly with height. The descent for the NH winter of 1992-1993 and the SH winter of 1992 computed with a three-dimensional trajectory model using the same radiation code was within 1 to 2 km of that calculated by the one-dimensional model, thus validating the vortex averaging procedure. The computed descent rates generally agree well with observations of long-lived tracers, thus validating the radiative transfer model.

Description: Spatial correlation among densely packed particles can substantially change their single-scattering properties, thus making questionable the applicability of the independent scattering approximation in calculations of light scattering by planetary regoliths. The same problem arises in geophysics in light scattering computations for snow, frosts, and bare soil. In this paper, we use a dense-medium light-scattering theory based on the introduction of the static structure factor to calculate asymmetry parameters of the phase function for densely packed particles with real refractive indices 1.31 and 1.66, approximating water ice and soil particles, respectively, and imaginary refractive indices 0, 0.01, and 0.3. For sparsely distributed, independently scattering grains, the calculated asymmetry parameters are always positive and always larger than those for densely packed particles. For densely packed grains, the asymmetry parameters may be negative but only for radius-to-wavelength ratios from about 0.1 to about 0.4. With decreasing particle size, the calculated asymmetry parameters tend to zero independently of the compaction state. In the geometrical optics regime, the asymmetry parameters for densely packed scatterers are positive and very close to those for independently scattering grains. These results may have important implications for remote sensing of the Earth and solid planetary surfaces. In particular, it is demonstrated that negative asymmetry parameters derived with some approximate multiple-scattering theories may be physically irrelevant and can be the result of using an inaccurate bidirectional reflection function combined with the ill-conditionally of the inverse scattering problem.

Description: The map of the coseismic displacement field generated by interferometric processing of synthetic aperture radar (SAR) images taken before and after the June 28, 1992, Landers earthquake sequence brings new insights into the nature of deformation caused by these earthquakes. We use the interferometric map generated by Massonnet et al. (1993) to analyze the surface displacement field in the vicinity of the fault trace. Complexities in the fringe pattern near the fault reflect short-wavelength variations of the surface rupture and slip distribution, and attest to large displacement gradients. Along two sections of the fault, characteristic fringe patterns can be recognized, contrasting in density and direction with patterns observed away from the rupture. In order to understand the observed fringe patterns, we compute synthetic interferograms in three simple cases: (1) rigid-body rotations about a vertical axis, (2) about a horizontal axis (tilt), and (3) distributed, simple shear. The orientation and spatial separation of interferometric fringes predicted by these models help constrain near-field deformation and rupture parameters. Where the Kickapoo fault connects with the Homestead Valley fault, the interferogram shows a clear pattern of parallel N20 deg W fringes separated by about 160 m. This pattern and vertical offsets measured along the Kickapoo fault suggest that the block between this fault and the Johnson Valley fault may have been tilted, down to the west. A 5-km block lifted by 1 m on one side would be tilted by an angle of 0.01 deg (190 microrad), producing fringes separated by about 160 m, parallel to the tilt axis. Such a tilt, parallel to a N20 deg W direction, would account for the gradual, northward increase of the vertical slip component observed along the Kickapoo fault. This tilt may also explain the 1 m of reverse slip observed along the 'slip gap' section of the Homestead Valley break. Between the southern end of the Johnson Valley fault and the Eureka Peak fault, where no surface rupture has been mapped, the dense pattern of fringes implies distributed shear, probably resulting from fault slip at depth. The density and direction of the fringes in the gap are consistent with a right-lateral slip of 1.2-3.8 m on a blind fault locked above the depth of 1.5-2 km. Such observations of small wavelength features in the SAR interferogram bring new insights into the near-field displacement gradient and thus on response of the uppermost crust to seismic rupture.

Description: Microwave radar and radiometer measurements of grasslands indicate a substantial reduction in sensor sensitivity to soil moisture in the presence of a thatch layer. When this layer is wet it masks changes in the underlying soil, making the canopy appear warm in the case of passive sensors (radiometer) and decreasing backscatter in the active case (scatterometer). A model for a grass canopy with thatch will be presented in this paper to explain this behavior and to compare with observations. The canopy model consists of three layers: grass, thatch, and the underlying soil. The grass blades are modeled by elongated elliptical discs and the thatch is modeled as a collection of disk shaped water droplets (i.e., the dry matter is neglected). The ground is homogeneous and flat. The distorted Born approximation is used to compute the radar cross section of this three layer canopy and the emissivity is computed from the radar cross section using the Peake formulation for the passive problem. Results are computed at L-band (1.4 GHz) and C-band (4.75 GHz) using canopy parameters (i.e., plant geometry, soil moisture, plant moisture, etc.) representative of Konza Prairie grasslands. The results are compared to C-band scatterometer measurements and L-band radiometer measurements at these grasslands.

Description: The acceleration of protons in a dynamically evolving magnetotail is investigated by tracing particles in the fields obtained from a three-dimensional resistive magnetohydrodynamic (MHD) simulation. The MHD simulation, representing plasmoid formation and ejection through a near-Earth reconnection process, leads to cross-tail electric fields of up to approximately 4 mV/m with integrated voltages across the tail of up to approximately 200 kV. Energization of particles takes place over a wide range along the tail, due to the large spatial extent of the increased electric field together with the finite cross-tail extent of the electric field region. Such accelerated particles appear earthward of the neutral line over a significant portion of the closed field line region inside of the separatrix, not just in the vicinity of the separatrix. Two different acceleration processes are identified: a 'quasi-potential' acceleration, due to particle motion in the direction of the cross-tail electric field, and a 'quasi-betatron' effect, which consists of multiple energy gains from repeated crossings of the acceleration region, mostly on Speiser-type orbits, in the spatially varying induced electric field. The major source region for accelerated particles in the hundreds of keV range is the central plasma sheet at the dawn flank outside the reconnection site. Since this source plasma is already hot and dense, its moderate energization by a factor of approximately 2 may be sufficient to explain the observed increases in the energetic particle fluxes. Particles from the tail are the source of beams at the plasma sheet/lobe boundary. The temporal increase in the energetic particle fluxes, estimated from the increase in energy gain, occurs on a fast timescale of a few minutes, coincident with a strong increase in B(sub z), despite the fact that the inner boundary ('injection boundary') of the distribution of energized particles is fairly smooth.

Description: The first comprehensive discussion of the south seasonal polar cap spectra obtained by the Mariner 7 infrared spectrometer in the short-wavelength region (2-4 microns) is presented. The infrared spectra is correlated with images acquired by the wide-angle camera. Significant spectral variation is noted in the cap interior and regions of varying water frost abundance, CO2 ice/frost cover, and CO2-ice path length can be distinguished. Many of these spectral variations correlate with heterogeneity noted in the camera images, but certain significant infrared spectral variations are not discernible in the visible. Simple reflectance models are used to classify the observed spectral variations into four regions. Region I is at the cap edge, where there is enhanced absorption beyond 3 microns inferred to be caused by an increased abundance of water frost. The increase in water abundance over that in the interior is on the level of a few parts per thousand or less. Region II is the typical cap interior characterized by spectral features of CO2 ice at grain sizes of several millimeters to centimeters. These spectra also indicate the presence of water frost at the parts per thousand level. A third, unusual region (III), is defined by three spectra in which weak CO2 absorption features are as much as twice as strong as in the average cap spectra and are assumed to be caused by an increased path length in the CO2. Such large paths are inconsistent with the high reflectance in the visible and at 2.2 microns and suggest layered structures or deposition conditions that are not accounted for in current reflectance models. The final region (IV) is an area of thinning frost coverage or transparent ice well in the interior of the seasonal cap. These spectra are a combination of CO2 and ground signatures.

Description: A high-precision radiometric satellite tracking system, Doppler Orbitography and Radio-positioning Integrated by Satellite system (DORIS), has recently been developed by the French space agency, Centre National d'Etudes Spatiales (CNES). DORIS was designed to provide tracking support for missions such as the joint United States/French TOPEX/Poseidon. As part of the flight testing process, a DORIS package was flown on the French SPOT 2 satellite. A substantial quantity of geodetic quality tracking data was obtained on SPOT 2 from an extensive international DORIS tracking network. These data were analyzed to assess their accuracy and to evaluate the gravitational modeling enhancements provided by these data in combination with the Goddard Earth Model-T3 (GEM-T3) gravitational model. These observations have noise levels of 0.4 to 0.5 mm/s, with few residual systematic effects. Although the SPOT 2 satellite experiences high atmospheric drag forces, the precision and global coverage of the DORIS tracking data have enabled more extensive orbit parameterization to mitigate these effects. As a result, the SPOT 2 orbital errors have been reduced to an estimated radial accuracy in the 10-20 cm RMS range. The addition of these data, which encompass many regions heretofore lacking in precision satellite tracking, has significantly improved GEM-T3 and allowed greatly improved orbit accuracies for Sun-synchronous satellites like SPOT 2 (such as ERS 1 and EOS). Comparison of the ensuing gravity model with other contemporary fields (GRIM-4C2, TEG2B, and OSU91A) provides a means to assess the current state of knowledge of the Earth's gravity field. Thus, the DORIS experiment on SPOT 2 has provided a strong basis for evaluating this new orbit tracking technology and has demonstrated the important contribution of the DORIS network to the success of the TOPEX/Poseidon mission.

Description: The presence of anisotropic plasma distributions, trapped at the Earth's magnetic equator, has consequences for the electric field structure which must exist in equilibrium along the magnetic field line. Data from SCATHA and Dynamics Explorer 1 indicate that the core ion distributions at the magnetic equator can be well described as bi-Maxwellian distributions, with a perpendicular temperature an order of magnitude larger than the parallel temperature. A collisionless model is developed for the variation in plasma parameters, following the forms developed by Whipple (1977). If the core electron anisotropy is low, the resulting electric field of approximately 0.1 microV/m is pointed away from the equator. Under these conditions the self-consistent electric field will not overcome the effects of magnetic trapping. The resulting potential distribution results in a local maximum in total plasma density at the equator. Only when the electron distribution is primarily field-aligned can there be a density minimum at the equator. Comparisons are made between this model and the observed variations in DE 1 plasma parameters with latitude.

Description: We report new measurements of the sodium emission intensity seen in a line of sight just above the surface of the Moon. These data show a strong dependence on lunar phase. The emission intensity decreases from a maximum around first quarter (phase angle 90 deg) to very small values near full Moon (phase angle 0 deg). This suggests that the rate of sodium vapor production from the lunar surface is largest at the subsolar point and becomes small near the terminator. However, the sodium emission near full Moon falls below that which would be expected for solar photon-driven processes. Since the solar wind flux decreases substantially when the Moon enters the Earth's magnetotail near full Moon, while the global solar photon flux is undiminished, we suggest that solar wind sputtering is the dominant process for sodium production.

Description: A combined 3-dimensional circulation model and aerosol microphysical/transport model is used to simulate the dispersion of the Mt. Pinatubo volcanic cloud in the stratosphere for the first few months following the eruption. Radiative heating of the cloud due to upwelling infrared radiation from the troposphere is shown to be an important factor affecting the transport. Without cloud heating, cloud material stays mostly north of the equator, whereas with cloud heating, the cloud is transported southward across the equator within the first two weeks following the eruption. Generally the simulations agree with Total Ozone Mapping Spectrometer (TOMS), Advanced Very High Resolution Radiometer (AVHRR), and Stratospheric Aerosol and Gas Experiment (SAGE) satellite observations showing the latitude distribution of cloud material to be between about 20 deg S and 30 deg N within the first few months. Temperature perturbations in the stratosphere induced by the aerosol heating are generally 1-4 K, in the range of those observed.

Description: We use geodetic data from Very Long Baseline Interferometry (VLBI) to determine the pre- and postseismic velocities of two sites. We then place limits on variations in interseismic strain buildup. The 1987 and 1988 Gulf of Alaska earthquakes (each Ms = 7.6) broke the Pacific plate interior. During the earthquakes the Cape Yakataga site moved 78 mm toward southwest. During the 1989 Loma Prieta earthquake (Ms = 7.1) the Fort Ord site moved 48 mm toward north. Baselines (a) from Fairbanks to Cape Yakataga and (b) from Mojave to Fort Ord change at nearly the same rate before and after the earthquakes. Postseismic transients, which we determine from differences between post- and preseismic rates, are minor: at Cape Yakataga the transient is 3 +/- 4 mm in a postseismic interval of 23 months, and at Fort Ord the transient is 6 +/- 5 mm in 21 months. The slip beneath the Loma Prieta rupture needed to generate the Fort Ord transient is 0.22 +/- 0.19 m, one-tenth the coseismic slip (2 m). We analyze elastic lithosphere-viscous asthenosphere models to determine that the characteristic time describing exponential decay in deep fault slip is longer than 6 years. The VLBI measurements are consistent with uniform interseismic strain buildup. They disagree with fast postseismic rates caused by an asthenosphere with very low viscosity.

Description: The AEPI (Atmospheric Emissions Photometric Imager) experiment on the ATLAS-1 shuttle mission (launched March 24, 1992) imaged the earth night airglow emission in O2 Atmospheric (0,0) bands, at 762.0 nm. Earthward views of O2 A bands show structure from gravity waves which exhibit extended horizontal structure with horizontal wavelenghts on the order of 50-100 km. These observations of the O2 A (0,0) bands are particularly interesting since in this wavelength the lower atmosphere absorbs all the earth-reflected emissions and most of the spectrally diffuse backgrounds. Herein we present observations of gravity waves using a topside airglow imaging technique.

Description: Fujii et al. (1994) obtained characteristics of the electrodynamic parameters, that is, field-aligned currents, electric fields, and electron precipitation, which are associted with auroral substorm events in the nighttime sector, through a unique analysis that places the ionospheric measurements of these parameters into the context of a generic substorm determined from global auroral images. In this paper we investigate in considerably more detail the characteristics of the field-aligned currents using data from the same set of passes as the previuos study. We show for the first time that the net upward field-aligned currents throughout the surge and surge horn are sufficient to account for most if not all of the converging currents of the auroral electrojets. Current densities are largest in the surge and surge horn. Current region continuity does not appear to exist across the substorm bulge region. Much of the auroral substorm field-aligned current is composed of filamentary currents and finite current segments at large angles to each other. The westward electrojet may contain large gradients in intensity both in local time and latitude due to sets of localized field-aligned currents. The net downward current for several hours to the west of the surge is insufficient to account for the eastward electrojet, consistent with the concept that this electrojet originates primarily on the dayside. Our pattern of field-aligned currents associated with the surge has common features and also differs significantly from the patterns previously derived from data from radars and ground-based magnetometer arrays. Our pattern is considerably more complex, probably due to the much higher resolution in latitude of the satellite data. It is also larger in area, since our average substorm is much larger than those pertaining to the previous patterns, giving a substorm wedge considerably wider than that obtained from the radar and array data.

Description: Infrared diffuse reflectance spectra (2.53-25 microns) of some carbonaceous (C) chondrites were measured. The integrated intensity of the absorption bands near 3 microns caused by hydrous minerals were compared with the modal content of hydrous minerals for the meteorites. The CM and CI chondrites show larger values of the intergated intensity than those of the unique C chondrites Y82162, Y86720 and B7904, suggesting that the amount of hydrous minerals in the CM and CI chondrites is larger, which supports the contention that hydrous minerals were dehydrated by thermal metamorphism in the unique chondrites. Orgueil (CI) has the largest value of the integrated intensity among the C chondrites we measured and shows a sharp absorption band at 3685/cm (2.71 microns) that is not seen in the spectra of the CM chondrites. There is an excellent correlation between the observed hydrogen content in C chondrites and the integrated intensity. The CM chondrites show a wide variation in the strength of absorption bands at 1470/cm (6.8 microns), despite the similarity in absorption features near 3 micron for all CM chondites. The 1470/cm band could be due to the presence of some hydrocarbons but may also be a result of terrestrial alteration processes.

Description: An analytical model of the bidirectional reflectance for optically semi-infinite plant canopies has been extended to describe the reflectance of finite depth canopies contributions from the underlying soil. The model depends on 10 independent parameters describing vegetation and soil optical and structural properties. The model is inverted with a nonlinear minimization routine using directional reflectance data for lawn (leaf area index (LAI) is equal to 9.9), soybeans (LAI, 2.9) and simulated reflectance data (LAI, 1.0) from a numerical bidirectional reflectance distribution function (BRDF) model (Myneni et al., 1988). While the ten-parameter model results in relatively low rms differences for the BRDF, most of the retrieved parameters exhibit poor stability. The most stable parameter was the single-scattering albedo of the vegetation. Canopy albedo could be derived with an accuracy of less than 5% relative error in the visible and less than 1% in the near-infrared. Sensitivity were performed to determine which of the 10 parameters were most important and to assess the effects of Gaussian noise on the parameter retrievals. Out of the 10 parameters, three were identified which described most of the BRDF variability. At low LAI values the most influential parameters were the single-scattering albedos (both soil and vegetation) and LAI, while at higher LAI values (greater than 2.5) these shifted to the two scattering phase function parameters for vegetation and the single-scattering albedo of the vegetation. The three-parameter model, formed by fixing the seven least significant parameters, gave higher rms values but was less sensitive to noise in the BRDF than the full ten-parameter model. A full hemispherical reflectance data set for lawn was then interpolated to yield BRDF values corresponding to advanced very high resolution radiometer (AVHRR) scan geometries collected over a period of nine days. The resulting parameters and BRDFs are similar to those for the full sampling geometry, suggesting that the limited geometry of AVHRR measurements might be used to reliably retrieve BRDF and canopy albedo with this model.

Description: There have been many significant improvements in the public access to the Space Shuttle Earth Observations Photography Database. New information is provided for the user community on the recently released videodisc of this database. Topics covered included the following: earlier attempts; our first laser videodisc in 1992; the new laser videodisc in 1994; and electronic database access.

Description: The Challenge Awards are designed to provide a unique perspective to students gifted in the arts and humanities from which to understand scientific endeavor by giving students an opportunity to participate in an ongoing research project. In the graduate program, seven students who had participated in previous Challenge Awards programs were selected to help develop the tools for Earth observations for the astronauts on the Space Radar Laboratory (SRL) missions. The goal of the Challenge Awards program was to prepare a training manual for the astronauts on the SRL missions. This paper describes the observations to be made by the astronauts on the SRL missions. The emphasis is on the dynamic seasonal features of the Earth's surface and atmosphere which justify the need for more than one flight of the SRL. Complete notebooks of the sites, global seasonal patterns, examples of radar and the Measurement of Air Pollution from Satellites data, and shuttle photographs have been given to each of the SRL crews.

Description: Oxygen production from a lunar rock has been experimentally demonstrated for the first time. A 10 g sample of high-Ti basalt 70035 was reduced with hydrogen in seven experiments at temperatures of 900-1050 C and pressures of 14.7-150 psia. In all experiments, water evolution began almost immediately and was essentially complete in tens of minutes. Oxygen yields ranged from 2.93 to 4.61% of the starting sample weight, and showed weak dependence on temperature and pressure. Analysis of the solid samples demonstrated total reduction of Fe(2+) in ilmenite and small degrees of reduction in olivine and pyroxene. Ti O2 was also partially reduced to one or more suboxides. Data from these experiments provide a basis for predicting the yield of oxygen from lunar basalt as well as new constraints on natural reduction in the lunar regolith.

Description: MacSigma0 is an interactive tool for the Macintosh which allows you to display and make computations from radar data collected by the following sensors: the JPL AIRSAR, ERS-1, JERS-1, and Magellan. The JPL AIRSAR system is a multi-polarimetric airborne synthetic aperture radar developed and operated by the Jet Propulsion Laboratory. It includes the single-frequency L-band sensor mounted on the NASA CV990 aircraft and its replacement, the multi-frequency P-, L-, and C-band sensors mounted on the NASA DC-8. MacSigma0 works with data in the standard JPL AIRSAR output product format, the compressed Stokes matrix format. ERS-1 and JERS-1 are single-frequency, single-polarization spaceborne synthetic aperture radars launched by the European Space Agency and NASDA respectively. To be usable by MacSigma0, The data must have been processed at the Alaska SAR Facility and must be in the "low-resolution" format. Magellan is a spacecraft mission to map the surface of Venus with imaging radar. The project is managed by the Jet Propulsion Laboratory. The spacecraft carries a single-frequency, single-polarization synthetic aperture radar. MacSigma0 works with framelets of the standard MIDR CD-ROM data products. MacSigma0 provides four basic functions: synthesis of images (if necessary), statistical analysis of selected areas, analysis of corner reflectors as a calibration measure (if appropriate and possible), and informative mouse tracking. For instance, the JPL AIRSAR data can be used to synthesize a variety of images such as a total power image. The total power image displays the sum of the polarized and unpolarized components of the backscatter for each pixel. Other images which can be synthesized are HH, HV, VV, RL, RR, HHVV*, HHHV*, HVVV*, HHVV* phase and correlation coefficient images. For the complex and phase images, phase is displayed using color and magnitude is displayed using intensity. MacSigma0 can also be used to compute statistics from within a selected area. The statistics computed depend on the image type. For JPL AIRSAR data, the HH, HV, VV, HHVV* phase, and correlation coefficient means and standard deviation measures are calculated. The mean, relative standard deviation, minimum, and maximum values are calculated for all other data types. A histogram of the selected area is also calculated and displayed. The selected area can be rectangular, linear, or polygonal in shape. The user is allowed to select multiple rectangular areas, but not multiple linear or polygonal areas. The statistics and histogram are displayed to the user and can either be printed or saved as a text file. MacSigma0 can also be used to analyze corner reflectors as a measure of the calibration for JPL AIRSAR, ERS-1, and JERS-1 data types. It computes a theoretical radar cross section and the actual radar cross section for a selected trihedral corner reflector. The theoretical cross section, measured cross section, their ratio in dBs, and other information are displayed to the user and can be saved into a text file. For ERS-1, JERS-1, and Magellan data, MacSigma0 simultaneously displays pixel location in data coordinates and in latitude, longitude coordinates. It also displays sigma0, the incidence angle (for Magellan data), the original pixel value (for Magellan data), and the noise power value (for ERS-1 and JERS-1 data). Grey scale computed images can be saved in a byte format (a headerless format which saves the image as a string of byte values) or a PICT format (a standard format readable by other image processing programs for the Macintosh). Images can also be printed. MacSigma0 is written in C-language for use on Macintosh series computers. The minimum configuration requirements for MacSigma0 are System 6.0, Finder 6.1, 1Mb of RAM, and at least a 4-bit color or grey-scale graphics display. MacSigma0 is also System 7 compatible. To compile the source code, Apple's Macintosh Programmers Workbench (MPW) 3.2 and the MPW C language compiler version 3.2 are required. The source code will not compile with a later version of the compiler; however, the compiled application which will run under the minimum hardware configuration is provided on the distribution medium. In addition, the distribution media includes an executable which runs significantly faster but requires a 68881 compatible math coprocessor and a 68020 compatible CPU. Since JPL AIRSAR data files can be very large, it is often desirable to reduce the size of a data file before transferring it to the Macintosh for use in MacSigma0. A small FORTRAN program which can be used for this purpose is included on the distribution media. MacSigma0 will print statistics on any output device which supports QuickDraw, and it will print images on any device which supports QuickDraw or PostScript. The standard distribution medium for MacSigma0 is a set of five 1.4Mb Macintosh format diskettes. This program was developed in 1992 and is a copyrighted work with all copyright vested in NASA. Version 4.2 of MacSigma0 was released in 1993.

Description: VICAR (Video Image Communication and Retrieval) is a general purpose image processing software system that has been under continuous development since the late 1960's. Originally intended for data from the NASA Jet Propulsion Laboratory's unmanned planetary spacecraft, VICAR is now used for a variety of other applications including biomedical image processing, cartography, earth resources, and geological exploration. The development of this newest version of VICAR emphasized a standardized, easily-understood user interface, a shield between the user and the host operating system, and a comprehensive array of image processing capabilities. Structurally, VICAR can be divided into roughly two parts; a suite of applications programs and an executive which serves as the interfaces between the applications, the operating system, and the user. There are several hundred applications programs ranging in function from interactive image editing, data compression/decompression, and map projection, to blemish, noise, and artifact removal, mosaic generation, and pattern recognition and location. An information management system designed specifically for handling image related data can merge image data with other types of data files. The user accesses these programs through the VICAR executive, which consists of a supervisor and a run-time library. From the viewpoint of the user and the applications programs, the executive is an environment that is independent of the operating system. VICAR does not replace the host computer's operating system; instead, it overlays the host resources. The core of the executive is the VICAR Supervisor, which is based on NASA Goddard Space Flight Center's Transportable Applications Executive (TAE). Various modifications and extensions have been made to optimize TAE for image processing applications, resulting in a user friendly environment. The rest of the executive consists of the VICAR Run-Time Library, which provides a set of subroutines (image I/O, label I/O, parameter I/O, etc.) to facilitate image processing and provide the fastest I/O possible while maintaining a wide variety of capabilities. The run-time library also includes the Virtual Raster Display Interface (VRDI) which allows display oriented applications programs to be written for a variety of display devices using a set of common routines. (A display device can be any frame-buffer type device which is attached to the host computer and has memory planes for the display and manipulation of images. A display device may have any number of separate 8-bit image memory planes (IMPs), a graphics overlay plane, pseudo-color capabilities, hardware zoom and pan, and other features). The VRDI supports the following display devices: VICOM (Gould/Deanza) IP8500, RAMTEK RM-9465, ADAGE (Ikonas) IK3000 and the International Imaging Systems IVAS. VRDI's purpose is to provide a uniform operating environment not only for an application programmer, but for the user as well. The programmer is able to write programs without being concerned with the specifics of the device for which the application is intended. The VICAR Interactive Display Subsystem (VIDS) is a collection of utilities for easy interactive display and manipulation of images on a display device. VIDS has characteristics of both the executive and an application program, and offers a wide menu of image manipulation options. VIDS uses the VRDI to communicate with display devices. The first step in using VIDS to analyze and enhance an image (one simple example of VICAR's numerous capabilities) is to examine the histogram of the image. The histogram is a plot of frequency of occurrence for each pixel value (0 - 255) loaded in the image plane. If, for example, the histogram shows that there are no pixel values below 64 or above 192, the histogram can be "stretched" so that the value of 64 is mapped to zero and 192 is mapped to 255. Now the user can use the full dynamic range of the display device to display the data and better see its contents. Another example of a VIDS procedure is the JMOVIE command, which allows the user to run animations interactively on the display device. JMOVIE uses the concept of "frames", which are the individual frames which comprise the animation to be viewed. The user loads images into the frames after the size and number of frames has been selected. VICAR's source languages are primarily FORTRAN and C, with some VAX Assembler and array processor code. The VICAR run-time library is designed to work equally easily from either FORTRAN or C. The program was implemented on a DEC VAX series computer operating under VMS 4.7. The virtual memory required is 1.5MB. Approximately 180,000 blocks of storage are needed for the saveset. VICAR (version 2.3A/3G/13H) is a copyrighted work with all copyright vested in NASA and is available by license for a period of ten (10) years to approved licensees. This program was developed in 1989.

Description: The Interactive Image Display Program (IMDISP) is an interactive image display utility for the IBM Personal Computer (PC, XT and AT) and compatibles. Until recently, efforts to utilize small computer systems for display and analysis of scientific data have been hampered by the lack of sufficient data storage capacity to accomodate large image arrays. Most planetary images, for example, require nearly a megabyte of storage. The recent development of the "CDROM" (Compact Disk Read-Only Memory) storage technology makes possible the storage of up to 680 megabytes of data on a single 4.72-inch disk. IMDISP was developed for use with the CDROM storage system which is currently being evaluated by the Planetary Data System. The latest disks to be produced by the Planetary Data System are a set of three disks containing all of the images of Uranus acquired by the Voyager spacecraft. The images are in both compressed and uncompressed format. IMDISP can read the uncompressed images directly, but special software is provided to decompress the compressed images, which can not be processed directly. IMDISP can also display images stored on floppy or hard disks. A digital image is a picture converted to numerical form so that it can be stored and used in a computer. The image is divided into a matrix of small regions called picture elements, or pixels. The rows and columns of pixels are called "lines" and "samples", respectively. Each pixel has a numerical value, or DN (data number) value, quantifying the darkness or brightness of the image at that spot. In total, each pixel has an address (line number, sample number) and a DN value, which is all that the computer needs for processing. DISPLAY commands allow the IMDISP user to display all or part of an image at various positions on the display screen. The user may also zoom in and out from a point on the image defined by the cursor, and may pan around the image. To enable more or all of the original image to be displayed on the screen at once, the image can be "subsampled." For example, if the image were subsampled by a factor of 2, every other pixel from every other line would be displayed, starting from the upper left corner of the image. Any positive integer may be used for subsampling. The user may produce a histogram of an image file, which is a graph showing the number of pixels per DN value, or per range of DN values, for the entire image. IMDISP can also plot the DN value versus pixels along a line between two points on the image. The user can "stretch" or increase the contrast of an image by specifying low and high DN values; all pixels with values lower than the specified "low" will then become black, and all pixels higher than the specified "high" value will become white. Pixels between the low and high values will be evenly shaded between black and white. IMDISP is written in a modular form to make it easy to change it to work with different display devices or on other computers. The code can also be adapted for use in other application programs. There are device dependent image display modules, general image display subroutines, image I/O routines, and image label and command line parsing routines. The IMDISP system is written in C-language (94%) and Assembler (6%). It was implemented on an IBM PC with the MS DOS 3.21 operating system. IMDISP has a memory requirement of about 142k bytes. IMDISP was developed in 1989 and is a copyrighted work with all copyright vested in NASA. Additional planetary images can be obtained from the National Space Science Data Center at (301) 286-6695.

Description: PC-SEAPAK is a user-interactive satellite data analysis software package specifically developed for oceanographic research. The program is used to process and interpret data obtained from the Nimbus-7/Coastal Zone Color Scanner (CZCS), and the NOAA Advanced Very High Resolution Radiometer (AVHRR). PC-SEAPAK is a set of independent microcomputer-based image analysis programs that provide the user with a flexible, user-friendly, standardized interface, and facilitates relatively low-cost analysis of oceanographic satellite data. Version 4.0 includes 114 programs. PC-SEAPAK programs are organized into categories which include CZCS and AVHRR level-1 ingest, level-2 analyses, statistical analyses, data extraction, remapping to standard projections, graphics manipulation, image board memory manipulation, hardcopy output support and general utilities. Most programs allow user interaction through menu and command modes and also by the use of a mouse. Most programs also provide for ASCII file generation for further analysis in spreadsheets, graphics packages, etc. The CZCS scanning radiometer aboard the NIMBUS-7 satellite was designed to measure the concentration of photosynthetic pigments and their degradation products in the ocean. AVHRR data is used to compute sea surface temperatures and is supported for the NOAA 6, 7, 8, 9, 10, 11, and 12 satellites. The CZCS operated from November 1978 to June 1986. CZCS data may be obtained free of charge from the CZCS archive at NASA/Goddard Space Flight Center. AVHRR data may be purchased through NOAA's Satellite Data Service Division. Ordering information is included in the PC-SEAPAK documentation. Although PC-SEAPAK was developed on a COMPAQ Deskpro 386/20, it can be run on most 386-compatible computers with an AT bus, EGA controller, Intel 80387 coprocessor, and MS-DOS 3.3 or higher. A Matrox MVP-AT image board with appropriate monitor and cables is also required. Note that the authors have received some reports of incompatibilities between the MVP-AT image board and ZENITH computers. Also, the MVP-AT image board is not necessarily compatible with 486-based systems; users of 486-based systems should consult with Matrox about compatibility concerns. Other PC-SEAPAK requirements include a Microsoft mouse (serial version), 2Mb RAM, and 100Mb hard disk space. For data ingest and backup, 9-track tape, 8mm tape and optical disks are supported and recommended. PC-SEAPAK has been under development since 1988. Version 4.0 was updated in 1992, and is distributed without source code. It is available only as a set of 36 1.2Mb 5.25 inch IBM MS-DOS format diskettes. PC-SEAPAK is a copyrighted product with all copyright vested in the National Aeronautics and Space Administration. Phar Lap's DOS_Extender run-time version is integrated into several of the programs; therefore, the PC-SEAPAK programs may not be duplicated. Three of the distribution diskettes contain DOS_Extender files. One of the distribution diskettes contains Media Cybernetics' HALO88 font files, also licensed by NASA for dissemination but not duplication. IBM is a registered trademark of International Business Machines. MS-DOS is a registered trademark of Microsoft Corporation. HALO88 is a registered trademark of Media Cybernetics, but the product was discontinued in 1991.

Description: The Land Analysis System (LAS) is an image analysis system designed to manipulate and analyze digital data in raster format and provide the user with a wide spectrum of functions and statistical tools for analysis. LAS offers these features under VMS with optional image display capabilities for IVAS and other display devices as well as the X-Windows environment. LAS provides a flexible framework for algorithm development as well as for the processing and analysis of image data. Users may choose between mouse-driven commands or the traditional command line input mode. LAS functions include supervised and unsupervised image classification, film product generation, geometric registration, image repair, radiometric correction and image statistical analysis. Data files accepted by LAS include formats such as Multi-Spectral Scanner (MSS), Thematic Mapper (TM) and Advanced Very High Resolution Radiometer (AVHRR). The enhanced geometric registration package now includes both image to image and map to map transformations. The over 200 LAS functions fall into image processing scenario categories which include: arithmetic and logical functions, data transformations, fourier transforms, geometric registration, hard copy output, image restoration, intensity transformation, multispectral and statistical analysis, file transfer, tape profiling and file management among others. Internal improvements to the LAS code have eliminated the VAX VMS dependencies and improved overall system performance. The maximum LAS image size has been increased to 20,000 lines by 20,000 samples with a maximum of 256 bands per image. The catalog management system used in earlier versions of LAS has been replaced by a more streamlined and maintenance-free method of file management. This system is not dependent on VAX/VMS and relies on file naming conventions alone to allow the use of identical LAS file names on different operating systems. While the LAS code has been improved, the original capabilities of the system have been preserved. These include maintaining associated image history, session logging, and batch, asynchronous and interactive mode of operation. The LAS application programs are integrated under version 4.1 of an interface called the Transportable Applications Executive (TAE). TAE 4.1 has four modes of user interaction: menu, direct command, tutor (or help), and dynamic tutor. In addition TAE 4.1 allows the operation of LAS functions using mouse-driven commands under the TAE-Facelift environment provided with TAE 4.1. These modes of operation allow users, from the beginner to the expert, to exercise specific application options. LAS is written in C-language and FORTRAN 77 for use with DEC VAX computers running VMS with approximately 16Mb of physical memory. This program runs under TAE 4.1. Since TAE 4.1 is not a current version of TAE, TAE 4.1 is included within the LAS distribution. Approximately 130,000 blocks (65Mb) of disk storage space are necessary to store the source code and files generated by the installation procedure for LAS and 44,000 blocks (22Mb) of disk storage space are necessary for TAE 4.1 installation. The only other dependencies for LAS are the subroutine libraries for the specific display device(s) that will be used with LAS/DMS (e.g. X-Windows and/or IVAS). The standard distribution medium for LAS is a set of two 9~track 6250 BPI magnetic tapes in DEC VAX BACKUP format. It is also available on a set of two TK50 tape cartridges in DEC VAX BACKUP format. This program was developed in 1986 and last updated in 1992.

Description: There is no simple and perfect way to measure residual stresses in metal parts that have been welded or deformed to make complex structures such as pressure vessels and aircraft, yet these locked-in stresses can contribute to structural failure by fatigue and fracture. However, one proven and tested technique for determining the internal stress of a metal part is to drill a test hole while measuring the relieved strains around the hole, such as the hole-drilling strain gage method described in ASTM E 837. The program HOLEGAGE processes strain gage data and provides additional calculations of internal stress variations that are not obtained with standard E 837 analysis methods. The typical application of the technique uses a three gage rosette with a special hole-drilling fixture for drilling a hole through the center of the rosette to produce a hole with very small gage pattern eccentricity error. Another device is used to control the drilling and halt the drill at controlled depth steps. At each step, strains from all three strain gages are recorded. The influence coefficients used by HOLEGAGE to compute stresses from relieved hole strains were developed by published finite element method studies of thick plates for specific hole sizes and depths. The program uses a parabolic fit and an interpolating scheme to project the coefficients to other hole sizes and depths. Additionally, published experimental data are used to extend the coefficients to relatively thin plates. These influence coefficients are used to compute the stresses in the original part from the strain data. HOLEGAGE will compute interior planar stresses using strain data from each drilled hole depth layer. Planar stresses may be computed in three ways including: a least squares fit for a linear variation with depth, an integral method to give incremental stress data for each layer, or by a linear fit to the integral data (with some surface data points omitted) to predict surface stresses before strain gage sanding preparations introduced additional residual stresses. Options are included for estimating the effect of hole eccentricity on calculations, smoothing noise from the strain data, and inputting the program data either interactively or from a data file. HOLEGAGE was written in FORTRAN 77 for DEC VAX computers under VMS, and is transportable except for system-unique TIME and DATE system calls. The program requires 54K of main memory and was developed in 1990. The program is available on a 9-track 1600 BPI VAX BACKUP format magnetic tape (standard media) or a TK50 tape cartridge. The documentation is included on the tape. DEC VAX and VMS are trademarks of Digital Equipment Corporation.

Description: Spaceborne Synthetic Aperture Radar (SAR) images are useful for planetary mapping and Earth sciences investigations. However, swath widths rarely exceed 100 Kilometers, and images must be patched together to create a mosaic in order to analyze larger areas. The primary function of this program is to generate large digital mosaics of SAR imagery without manually marked tiepoints. MOSK can produce multiframe mosaics by combining images in the along-track, adjacent cross-track swaths, or ascending and descending passes. Geocoded map registered images, such as the ones produced by MAPJTC (NPO-17718), are required as input. The output is a geocoded mosaic on a standard map grid which permits easy registration with other geocoded data sets. Mosaicking of geocoded SAR imagery involves three steps. First, a match point is selected at the center of the overlapping area, then an image patch around the match point is extracted from both images and cross-correlation is done on this area. Then, images with their refined match points are merged together to form a mosaic. To handle the large data volume of overlapping intermediate stages, large mosaics are divided into equal size quadrants with each quadrant cut from an intermediate mosaic. The full mosaic can then be assembled from the individual quadrants. Finally, radiometric disparities at the image seams are smoothed by a "feathering" technique. The automatic mosaic system generates output with minimal operator interaction. However, manual tiepointing is required in cases of a large registration error or two images with smooth surfaces such as ocean images. MOSK is implemented on a DEC VAX 11/785 running VMS 4.5. Most subroutines are in FORTRAN, but three are in MAXL and one is in APAL. The program requires 1 Mb of memory and a Floating Point Systems AP-5210 array processor. The system memory usage is approximately 1000 pages and the requirement of page file size is 2000 blocks. MOSK was developed in 1988.

Description: MAPJTC was designed to rectify and transform the standard image output of the digital Synthetic Aperture Radar (SAR) correlator into a geocoded map registered image without operator interaction or manual tiepointing. This is accomplished by modeling the distortions and predicting the pixel displacements based on platform and radar parameters. The map projection implemented in MAPJTC is the Universal Transverse Mercator. Since the re-sampling operation is independent of the transformation data generation, other cartographic projections can be implemented with few software modifications. MAPJTC makes a precise determination of the geodetic location of an arbitrary pixel within the image frame based on the simultaneous solution of a set of earth model equations, SAR Doppler equations, and SAR range equations that identify the slant range from the sensor to the target at a specific image pixel. Based on a table of geodetic coordinates of the image pixels, the image is then mapped onto the desired cartographic projection by applying the appropriate transformation equations. Typically, mapping involves a two-dimensional re-sampling and is very computationally intensive. MAPJTC reduces the procedure to two one-dimensional passes, which saves computer time. Geocoding transforms the rectified image into a grid defined by a specific map projection. (The image is rotated and rectified to match the map projection.) Again, the two dimensional re-sampling process can be separated into two one-dimensional re-sampling processes. Optionally, MAPJTC can correct terrain-induced distortions in SAR imagery when a digital elevation map is available. MAPJTC was developed on a DEC VAX 11/785 under VMS 4.5. The program is written in FORTRAN (84%), APAL (2%), and MAXL (14%). It requires 6Mb of memory and a Floating Point Systems AP-5210 Array Processor equipped with 1Mb of memory. MAPJTC can run interactively or as a batch job. MAPJTC was developed in 1987.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The Spectral Analysis Manager (SPAM) was developed to allow easy qualitative analysis of multi-dimensional imaging spectrometer data. Imaging spectrometers provide sufficient spectral sampling to define unique spectral signatures on a per pixel basis. Thus direct material identification becomes possible for geologic studies. SPAM provides a variety of capabilities for carrying out interactive analysis of the massive and complex datasets associated with multispectral remote sensing observations. In addition to normal image processing functions, SPAM provides multiple levels of on-line help, a flexible command interpretation, graceful error recovery, and a program structure which can be implemented in a variety of environments. SPAM was designed to be visually oriented and user friendly with the liberal employment of graphics for rapid and efficient exploratory analysis of imaging spectrometry data. SPAM provides functions to enable arithmetic manipulations of the data, such as normalization, linear mixing, band ratio discrimination, and low-pass filtering. SPAM can be used to examine the spectra of an individual pixel or the average spectra over a number of pixels. SPAM also supports image segmentation, fast spectral signature matching, spectral library usage, mixture analysis, and feature extraction. High speed spectral signature matching is performed by using a binary spectral encoding algorithm to separate and identify mineral components present in the scene. The same binary encoding allows automatic spectral clustering. Spectral data may be entered from a digitizing tablet, stored in a user library, compared to the master library containing mineral standards, and then displayed as a timesequence spectral movie. The output plots, histograms, and stretched histograms produced by SPAM can be sent to a lineprinter, stored as separate RGB disk files, or sent to a Quick Color Recorder. SPAM is written in C for interactive execution and is available for two different machine environments. There is a DEC VAX/VMS version with a central memory requirement of approximately 242K of 8 bit bytes and a machine independent UNIX 4.2 version. The display device currently supported is the Raster Technologies display processor. Other 512 x 512 resolution color display devices, such as De Anza, may be added with minor code modifications. This program was developed in 1986.

Description: The ORBSIM program was developed for the accurate extraction of geophysical model parameters from Doppler radio tracking data acquired from orbiting planetary spacecraft. The model of the proposed planetary structure is used in a numerical integration of the spacecraft along simulated trajectories around the primary body. Using line of sight (LOS) Doppler residuals, ORBSIM applies fast and efficient modelling and optimization procedures which avoid the traditional complex dynamic reduction of data. ORBSIM produces quantitative geophysical results such as size, depth, and mass. ORBSIM has been used extensively to investigate topographic features on the Moon, Mars, and Venus. The program has proven particulary suitable for modelling gravitational anomalies and mascons. The basic observable for spacecraft-based gravity data is the Doppler frequency shift of a transponded radio signal. The time derivative of this signal carries information regarding the gravity field acting on the spacecraft in the LOS direction (the LOS direction being the path between the spacecraft and the receiving station, either Earth or another satellite). There are many dynamic factors taken into account: earth rotation, solar radiation, acceleration from planetary bodies, tracking station time and location adjustments, etc. The actual trajectories of the spacecraft are simulated using least squares fitted to conic motion. The theoretical Doppler readings from the simulated orbits are compared to actual Doppler observations and another least squares adjustment is made. ORBSIM has three modes of operation: trajectory simulation, optimization, and gravity modelling. In all cases, an initial gravity model of curved and/or flat disks, harmonics, and/or a force table are required input. ORBSIM is written in FORTRAN 77 for batch execution and has been implemented on a DEC VAX 11/780 computer operating under VMS. This program was released in 1985.

Description: The PYROLASER package is an operating system for the Pyrometer Instrument Company's Pyrolaser. There are 6 individual programs in the PYROLASER package: two main programs, two lower level subprograms, and two programs which, although independent, function predominantly as macros. The package provides a quick and easy way to setup, control, and program a standard Pyrolaser. Temperature and emissivity measurements may be either collected as if the Pyrolaser were in the manual operations mode, or displayed on real time strip charts and stored in standard spreadsheet format for post-test analysis. A shell is supplied to allow macros, which are test-specific, to be easily added to the system. The Pyrolaser Simple Operation program provides full on-screen remote operation capabilities, thus allowing the user to operate the Pyrolaser from the computer just as it would be operated manually. The Pyrolaser Simple Operation program also allows the use of "quick starts". Quick starts provide an easy way to permit routines to be used as setup macros for specific applications or tests. The specific procedures required for a test may be ordered in a sequence structure and then the sequence structure can be started with a simple button in the cluster structure provided. One quick start macro is provided for continuous Pyrolaser operation. A subprogram, Display Continuous Pyr Data, is used to display and store the resulting data output. Using this macro, the system is set up for continuous operation and the subprogram is called to display the data in real time on strip charts. The data is simultaneously stored in a spreadsheet format. The resulting spreadsheet file can be opened in any one of a number of commercially available spreadsheet programs. The Read Continuous Pyrometer program is provided as a continuously run subprogram for incorporation of the Pyrolaser software into a process control or feedback control scheme in a multi-component system. The program requires the Pyrolaser to be set up using the Pyrometer String Transfer macro. It requires no inputs and provides temperature and emissivity as outputs. The Read Continuous Pyrometer program can be run continuously and the data can be sampled as often or as seldom as updates of temperature and emissivity are required. PYROLASER is written using the Labview software for use on Macintosh series computers running System 6.0.3 or later, Sun Sparc series computers running OpenWindows 3.0 or MIT's X Window System (X11R4 or X11R5), and IBM PC or compatibles running Microsoft Windows 3.1 or later. Labview requires a minimum of 5Mb of RAM on a Macintosh, 24Mb of RAM on a Sun, and 8Mb of RAM on an IBM PC or compatible. The Labview software is a product of National Instruments (Austin,TX; 800-433-3488), and is not included with this program. The standard distribution medium for PYROLASER is a 3.5 inch 800K Macintosh format diskette. It is also available on a 3.5 inch 720K MS-DOS format diskette, a 3.5 inch diskette in UNIX tar format, and a .25 inch streaming magnetic tape cartridge in UNIX tar format. An electronic copy of the documentation in Macintosh WordPerfect version 2.0.4 format is included on the distribution medium. Printed documentation is included in the price of the program. PYROLASER was developed in 1992.

Description: MSIS is an empirical model of the thermosphere based on Mass Spectrometer and Incoherent Scatter data. It provides a description of atmospheric temperature, density, and composition for altitudes higher than 85 kilometers. There are coefficients in MSIS to account for yearly and daily variations, geodetic latitude and longitude, altitude, and solar activity. Variations due to magnetic storms are represented by three-hour magnetic ap indices. The MSIS model enables a more timely prediction of aeronomic densities for specific events, such as rocket flights. The database for this model is a comprehensive summary of rocket flight, satellite, incoherent scatter radar, grenade, and falling sphere measurements. Subsets of data were formed by random selection after sorting on altitude, latitude, time of day, etc. Curve fitting was done with four to five thousand data points at a time. The resulting coefficients are presented in subroutines which calculate thermospheric composition and temperature for a user-supplied position and time. MSIS is written in FORTRAN 77 for use with batch or interactive programs and has been implemented on a DEC VAX series computer operating under VMS 4.3 with a central memory requirement of approximately 25K of 8 bit bytes. MSIS is based on a 1977 thermosphere model and was last updated in 1987 to reflect the CIRA 1986 Neutral Thermosphere Model.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: The Science and Technology Laboratory Applications Software (ELAS) was originally designed to analyze and process digital imagery data, specifically remotely-sensed scanner data. This capability includes the processing of Landsat multispectral data; aircraft-acquired scanner data; digitized topographic data; and numerous other ancillary data, such as soil types and rainfall information, that can be stored in digitized form. ELAS has the subsequent capability to geographically reference this data to dozens of standard, as well as user created projections. As an integrated image processing system, ELAS offers the user of remotely-sensed data a wide range of capabilities in the areas of land cover analysis and general purpose image analysis. ELAS is designed for flexible use and operation and includes its own FORTRAN operating subsystem and an expandable set of FORTRAN application modules. Because all of ELAS resides in one "logical" FORTRAN program, data inputs and outputs, directives, and module switching are convenient for the user. There are over 230 modules presently available to aid the user in performing a wide range of land cover analyses and manipulation. The file management modules enable the user to allocate, define, access, and specify usage for all types of files (ELAS files, subfiles, external files etc.). Various other modules convert specific types of satellite, aircraft, and vector-polygon data into files that can be used by other ELAS modules. The user also has many module options which aid in displaying image data, such as magnification/reduction of the display; true color display; and several memory functions. Additional modules allow for the building and manipulation of polygonal areas of the image data. Finally, there are modules which allow the user to select and classify the image data. An important feature of the ELAS subsystem is that its structure allows new applications modules to be easily integrated in the future. ELAS has as a standard the flexibility to process data elements exceeding 8 bits in length, including floating point (noninteger) elements and 16 or 32 bit integers. Thus it is able to analyze and process "non-standard" nonimage data. The VAX (ERL-10017) and Concurrent (ERL-10013) versions of ELAS 9.0 are written in FORTRAN and ASSEMBLER for DEC VAX series computers running VMS and Concurrent computers running MTM. The Sun (SSC-00019), Masscomp (SSC-00020), and Silicon Graphics (SSC-00021) versions of ELAS 9.0 are written in FORTRAN 77 and C-LANGUAGE for Sun4 series computers running SunOS, Masscomp computers running UNIX, and Silicon Graphics IRIS computers running IRIX. The Concurrent version requires at least 15 bit addressing and a direct memory access channel. The VAX and Concurrent versions of ELAS both require floating-point hardware, at least 1Mb of RAM, and approximately 70Mb of disk space. Both versions also require a COMTAL display device in order to display images. For the Sun, Masscomp, and Silicon Graphics versions of ELAS, the disk storage required is approximately 115Mb, and a minimum of 8Mb of RAM is required for execution. The Sun version of ELAS requires either the X-Window System Version 11 Revision 4 or Sun OpenWindows Version 2. The Masscomp version requires a GA1000 display device and the associated "gp" library. The Silicon Graphics version requires Silicon Graphics' GL library. ELAS display functions will not work with a monochrome monitor. The standard distribution medium for the VAX version (ERL~10017) is a set of two 9-track 1600 BPI magnetic tapes in DEC VAX BACKUP format. This version is also available on a TK50 tape cartridge in DEC VAX BACKUP format. The standard distribution medium for the Concurrent version (ERL-10013) is a set of two 9-track 1600 BPI magnetic tapes in Concurrent BACKUP format. The standard distribution medium for the Sun version (SSC-00019) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Masscomp version, (SSC-00020) is a .25 inch streaming magnetic tape cartridge in UNIX tar format. The standard distribution medium for the Silicon Graphics version (SSC-00021) is a .25 inch streaming magnetic IRIS tape cartridge in UNIX tar format. Version 9.0 was released in 1991. Sun4, SunOS, and Open Windows are trademarks of Sun Microsystems, Inc. MIT X Window System is licensed by Massachusetts Institute of Technology.

Description: IMAGEP is a FORTRAN computer algorithm containing various image processing, analysis, and enhancement functions. It is a keyboard-driven program organized into nine subroutines. Within the subroutines are other routines, also, selected via keyboard. Some of the functions performed by IMAGEP include digitization, storage and retrieval of images; image enhancement by contrast expansion, addition and subtraction, magnification, inversion, and bit shifting; display and movement of cursor; display of grey level histogram of image; and display of the variation of grey level intensity as a function of image position. This algorithm has possible scientific, industrial, and biomedical applications in material flaw studies, steel and ore analysis, and pathology, respectively. IMAGEP is written in VAX FORTRAN for DEC VAX series computers running VMS. The program requires the use of a Grinnell 274 image processor which can be obtained from Mark McCloud Associates, Campbell, CA. An object library of the required GMR series software is included on the distribution media. IMAGEP requires 1Mb of RAM for execution. The standard distribution medium for this program is a 1600 BPI 9~track magnetic tape in VAX FILES-11 format. It is also available on a TK50 tape cartridge in VAX FILES-11 format. This program was developed in 1991. DEC, VAX, VMS, and TK50 are trademarks of Digital Equipment Corporation.

Description: Scale-space filtering is used to screen information obtained from signals that produce a complex curve (such as that in geographic and thermal analysis) to gain a truer representation of the area under analysis. PSF extends this technique to extract non-periodic hills and valleys from a signal. Because the signal's information is sometimes too complex to determine with certainty if some features are real or artificial, PSF calculates probabilities, with the extracted features corresponding to real events, in order to aid in determining the signal's accuracy. Since the probabilities associated with the features are derived from domain-specific statistics, it is (most likely) necessary to modify the program code to correspond to the user's particular domain. PSF also provides a standard scale-space filtering algorithm for use when the desired features can be identified with certainty or when it is not practical to get the domain-specific statistics. The PSF algorithm is based on Witkin's scale-space filtering theory. The program detects signal variations by finding the points of inflection in the input signal. The number and position of these points are dependent upon the scale of the derivative operators used to detect them. Therefore, instead of assuming any single scale to be correct, PSF identifies points of inflection in a large number of different scales. It then describes the curve according to the groups of points of inflection, across all scales, caused by the same physical process. PSF provides an output table giving the following information: the abscissa of the first inflection of the peak, the type of peak, the distance between the first and second inflection points, the abscissa of the peak, and the probability of the feature corresponding to a real event in the curve. The program will also list points representing a graphical image of the signal and detected peaks. This data can be used with a standard plotting program (not included) to display the signal and its features graphically. PSF is written in C language (49%) and Common LISP (51%) for use on a Sun SPARC workstation running the UNIX operating system. PSF requires 4Mb of RAM. The standard distribution medium for this program is a .25 streaming magnetic tape cartridge in UNIX tar format. It is also available on a 3.5 inch diskette in UNIX tar format. PSF was developed in 1991. Sun and SPARC are trademarks of Sun Microsystems, Inc. UNIX is a registered trademark of AT&T Bell Laboratories.

Description: The complex environment of the typical research laboratory requires flexible process control. This program provides natural language process control from an IBM PC or compatible machine. Sometimes process control schedules require changes frequently, even several times per day. These changes may include adding, deleting, and rearranging steps in a process. This program sets up a process control system that can either run without an operator, or be run by workers with limited programming skills. The software system includes three programs. Two of the programs, written in FORTRAN77, record data and control research processes. The third program, written in Pascal, generates the FORTRAN subroutines used by the other two programs to identify the user commands with the user-written device drivers. The software system also includes an input data set which allows the user to define the user commands which are to be executed by the computer. To set the system up the operator writes device driver routines for all of the controlled devices. Once set up, this system requires only an input file containing natural language command lines which tell the system what to do and when to do it. The operator can make up custom commands for operating and taking data from external research equipment at any time of the day or night without the operator in attendance. This process control system requires a personal computer operating under MS-DOS with suitable hardware interfaces to all controlled devices. The program requires a FORTRAN77 compiler and user-written device drivers. This program was developed in 1989 and has a memory requirement of about 62 Kbytes.

Description: The Vibration Pattern Imager (VPI) system was designed to control and acquire data from laser vibrometer sensors. The PC computer based system uses a digital signal processing (DSP) board and an analog I/O board to control the sensor and to process the data. The VPI system was originally developed for use with the Ometron VPI Sensor (Ometron Limited, Kelvin House, Worsley Bridge Road, London, SE26 5BX, England), but can be readily adapted to any commercially available sensor which provides an analog output signal and requires analog inputs for control of mirror positioning. VPI's graphical user interface allows the operation of the program to be controlled interactively through keyboard and mouse-selected menu options. The main menu controls all functions for setup, data acquisition, display, file operations, and exiting the program. Two types of data may be acquired with the VPI system: single point or "full field". In the single point mode, time series data is sampled by the A/D converter on the I/O board at a user-defined rate for the selected number of samples. The position of the measuring point, adjusted by mirrors in the sensor, is controlled via a mouse input. In the "full field" mode, the measurement point is moved over a user-selected rectangular area with up to 256 positions in both x and y directions. The time series data is sampled by the A/D converter on the I/O board and converted to a root-mean-square (rms) value by the DSP board. The rms "full field" velocity distribution is then uploaded for display and storage. VPI is written in C language and Texas Instruments' TMS320C30 assembly language for IBM PC series and compatible computers running MS-DOS. The program requires 640K of RAM for execution, and a hard disk with 10Mb or more of disk space is recommended. The program also requires a mouse, a VGA graphics display, a Four Channel analog I/O board (Spectrum Signal Processing, Inc.; Westborough, MA), a break-out box and a Spirit-30 board (Sonitech International, Inc.; Wellesley, MA) which includes a TMS320C30 DSP processor, 256Kb zero wait state SRAM, and a daughter board with 8Mb one wait state DRAM. Please contact COSMIC for additional information on required hardware and software. In order to compile the provided VPI source code, a Microsoft C version 6.0 compiler, a Texas Instruments' TMS320C30 assembly language compiler, and the Spirit 30 run time libraries are required. A math co-processor is highly recommended. A sample MS-DOS executable is provided on the distribution medium. The standard distribution medium for this program is one 5.25 inch 360K MS-DOS format diskette. The contents of the diskettes are compressed using the PKWARE archiving tools. The utility to unarchive the files, PKUNZIP.EXE, is included. VPI was developed in 1991-1992.

Description: Continuous optical observations of cusp/cleft auroral activities within approximately equal to 09-15 MLT and 70-76 deg magnetic latitude are studied in relation to changes in solar wind dynamic pressure and interplanetary magnetic field (IMF) variability. The observed latitudinal movements of the cusp/cleft aurora in response to IMF B(sub z) changes may be explained as an effect of a variable magnetic field intensity in the outer dayside magnetosphere associated with the changing intensity of region 1 field-aligned currents and associated closure currents. Ground magnetic signatures related to such currents were observed in the present case (January 10, 1993). Strong, isolated enhancements in solar wind dynamic pressure (Delta p/p is greater than or equal to 0.5) gave rise to equatorward shifts of the cusp/cleft aurora, characteristic auroral transients, and distinct ground magnetic signatures of enhanced convection at cleft latitudes. A sequence of auroral events of approximately equal to 5-10 min recurrence time, moving eastward along the poleward boundary of the persistent cusp/cleft aurora in the approximately equal to 10-14 MLT sector, during negative IMF B(sub z) and B(sub y) conditions, were found to be correlated with brief pulses in solar wind dynamic pressure (0.1 is less than Delta p/p is less than 0.5). Simultaneous photometer observations from Ny Alesund, Svalbard, and Danmarkshavn, Greenland, show that the events often appeared on the prenoon side (approximately equal to 10-12 MLT), before moving into the postnoon sector in the case we study here, when IMF B(sub y) is less than 0. In other cases, similar auroral event sequences have been observed to move westward in the prenoon sector, during intervals of positive B(sub y). Thus a strong prenoon/postnoon asymmetry of event occurence and motion pattern related to the IMF B(sub y) polarity is observed. We find that this category of auroral event sequence is stimulated bursts of electron precipitation that originate from magnetosheath plasma that has accessed that dayside magnetosphere in the noon or near-noon sector, possibly at high latitudes, partly governed by the IMF orientation as well as by solar wind dynamic pressure pulses.