ALBERT

Description: Helmet streamers on the sun were observed to be the site of coronal mass ejections, dynamic events that eject coronal plasma and magnetic fields into the solar wind. A two dimensional (azimuthally symmetric) helmet streamer configuration was developed by computing solutions of the time dependent magnetohydrodynamic (MHD) equations, for a specified magnetic flux distribution on the sun. The helmet streamer is not symmetric about the equator. The evolution of the configuration, when differential rotation is applied, was investigated. It was found that after many rotations the configuration does not reach a steady state, but disrupts recurrently with the ejection of a plasmoid. These results suggest that differential rotation may be one of the mechanisms by which mass ejections are initiated.

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: The ground-based observing facilities of the National Solar Observatory (NSO) are reviewed from the perspective of joint observations with SOHO. A specific proposal is presented for observations of the HE-I 1083.0 nm line with the NASA/NSO spectromagnetograph and He 10830 video filtergraph/magnetograph in coordination with ultraviolet sensitive instruments on SOHO. The first task will be to look for associations of low-temperature transition-region lines with He 1083 nm absorption to investigate Andretta's conjecture, i.e. that the He 1083 nm line is formed in two layers where extreme ultraviolet radiation produced both in the low-temperature transition region (the upper layer) and in the surrounding corona products - a lower layer of absorption in the upper chromosphere.

Description: The white light coronograph (WLC) on Skylab provided an opportunity to study the corona at high spatial and temporal resolution. The spatial resolution of the instrument was approximately 25 cm with images taken approximately one per min. One set of images taken over a 10 min period was digitized, providing ten high spatial resolution images for analysis. The progress in data processing techniques available at the time was not sufficient to permit a reliable study of the fine structure in these images. Using current techniques an investigation of the sizes and lifetimes of the smallest scale features in the data was carried out. A preliminary analysis of an area between 2 and 3 Ro was completed. The results show that very narrow rays extend from at least 2 to 3 Ro. The narrowest of these rays has a thickness of approximately 75 cm. The contrast is so low that they are very close to the noise limit of the data. Most of the rays observed become unrecognizable after 10 min, although some remain visible over the entire time. Some notion seems to be detectable in the fine structure rays, but analysis of more frames will be needed to quantify these results.

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: There is a frequent need to measure the frequency stability and phase noise levels of very high performance signal sources that are required for certain spacecraft missions. These measurements need to be done at different locations as the spacecraft subsystems progress through the various stages of development, assembly, test, and integration. Allan Deviation and Phase Noise of high performance sources are generally measured by comparing the unit under test to a reference standard. Five basic requirements are associated with making these kind of measurements: (1) the reference standard performance needs to be equal or better than the unit under test; (2) the measurement system needs to accommodate odd, nonstandard measurement frequencies that can range from 4 MHz to 35 GHz; (3) warm-up frequency drift and aging can corrupt a measurement and must be dealt with; (4) test equipment generated noise must be understood and prevented from limiting the measurements; (5) test equipment noise performance must be verifiable in the field as needed. A portable measurement system that was built by JPL and used in the field is described. The methods of addressing the above requirements are outlined and some measurement noise floor values are given. This test set was recently used to measure state of the art crystal oscillator frequency standards on the TOPEX and MARS OBSERVER spacecraft during several stages of acceptance tests.

Description: A study was conducted to assess the effects of transmitting a precision clock synchronization signal over a commercial multiplexed fiber optic communication system. This study is an evaluation of the distortion and jitter introduced into the signal by this type of transmission system. An analysis comparing signal quality at the multiplexing and demultiplexing ends of the fiber optic communication system shows that the amplitude and phase distortion added to the clock synchronization signal by the transmission system is minimal.

Description: Historically, Maximum Time Interval Error (MTIE) and Maximum Relative Time Interval Error (MRTIE) have been the main measurement techniques used to characterize timing performance in telecommunications networks. Recently, a new measurement technique, Time Variance (TVAR) has gained acceptance in the North American (ANSI) standards body. TVAR was developed in concurrence with NIST to address certain inadequacies in the MTIE approach. The advantages and disadvantages of each of these approaches are described. Real measurement examples are presented to illustrate the critical issues in actual telecommunication applications. Finally, a new MTIE measurement is proposed (ZTIE) that complements TVAR. Together, TVAR and ZTIE provide a very good characterization of network timing.

Description: For proper operation of today's large digital private networks, high quality synchronization must be achieved. A general telecommunication performance objective is to maintain long-term frequency accuracy of ten parts per trillion at all synchronous digital equipment in the network. Many times, however, this is not achieved in private networks. Low quality clocks, errored transmission facilities, and incorrectly designed synchronization plans are often cause for poor performance. It is shown that properly-designed private networks can operate with long term frequency averages between ten parts per trillion to ten parts per million. These performance levels can adversely impact customer applications. The most demanding applications are digital and voice band data, encrypted voice, facsimile, and video. In a typical private network operating at 0.01 parts per million, the user would experience reduced data throughput, dropped encrypted calls, unreadable facsimile pages, or interrupted video transmission dozens of times per day. The major contribution to poor private network synchronization performance is the interaction of Customer Premises Equipment (CPE) clocks and the network facilities used to distribute timing. The performance of typical CPE clocks and facilities, and their impact on customer applications, are discussed. CPE clock performance issues, along with private network architectural constraints, make synchronization planning extremely difficult. Planning is usually costly and requires specialized expertise.

Description: A new method that uses round-trip paths to accurately measure transmission delay for time synchronization is proposed. The performance of the method in Synchronous Digital Hierarchy networks is discussed. The feature of this method is that it separately measures the initial round trip path delay and the variations in round-trip path delay. The delay generated in SDH equipment is determined by measuring the initial round-trip path delay. In an experiment with actual SDH equipment, the error of initial delay measurement was suppressed to 30ns.

Description: A coded time/date information dissemination system (CTD), based on telephone lines and commercial modems, is now in its experimental phase in Italy at IEN. This service, born from a cooperation with other metrological laboratories (TUG, Austria, SNT, Sweden, VSL, The Netherlands), represents an attempt towards an European standardization. Some results of an experimental analysis in which a few modems were tested, both in laboratory conditions and connected to the telephone network, in order to evaluate the timing capability of the system are given. When the system is used in a one-way mode, in many practical cases the modems delay turns out to be the main factor which limits the accuracy, even more than the telephone line delays. If the two-way mode is used, the modems asymmetry, i.e., the delay difference between transmission and reception, is almost always the most important source of uncertainty, provided the link is not including a space segment. Comparing the widely used V.22 modems to the old V.21 ones, the latters turn out to be better both in delay time (30-100 ms V.22, and 7-15 ms V.21) and asymmetry (10-50 micro-s V.22, and 10 ms V.22). Time transfer accuracies of 10 micron-s (same turn) to 100 micro-s (long distance calls) were obtained in two-way mode with commercial V.21 modems.

Description: Filaments, flare sprays, prominences and 'post-flare' loops are familiar to H alpha observers in their frequent appearances 'in absorption', dark against the chromospheric background or plages. Observations of the X-ray corona are generally interpreted as due to emission via optically thin thermal bremsstrahlung. Several cases of X-ray coronal structures in Yohkoh images, due to high opacity, absorbing matter in coronograph loops, are presented. The presence of the absorbing matter, mixed with emitting matter, complicates inference of physical parameters such as emission measures in X-ray sources. In the case of well defined features, absorption provides an opportunity to infer density. Quantitative estimates of the attenuation due to the absorption in example features are presented.

Description: The causal association of major solar particle events seen at earth with coronal mass ejections (CME's), and not with solar flares, is discussed. Evidence that led to the demise of the flare dominated paradigm for major solar energetic particle events are described. The possibility of distinguishing particles from impulsive and gradual events using only observations is described. Particle acceleration at the CME level is discussed. Multi-spacecraft observations of CME events are described. Concerning the interplanetary CME, bidirectional proton events are discussed. Conclusions from progress in understanding the characteristics of solar energetic particles and their relation to the physical mechanisms of acceleration are given.

Description: Two dimensional magnetohydrodynamic simulations of the distortion of a magnetic flux tube, accelerated through ambient solar wind plasma, are presented. Vortices form on the trailing edge of the flux tube, and couple strongly to its interior. If the flux tube azimuthal field is weak, it deforms into an elongated banana-like shape after a few Alfven transit times. A significant azimuthal field component inhibits this distortion. In the case of magnetic clouds in the solar wind, it is suggested that the shape observed at 1 AU was determined by distortion of the cloud in the inner heliosphere. Distortion of the cloud beyond 1 AU takes many days. It is estimated that effective drag coefficients slightly greater than unity are appropriate for modeling flux tube propagation. Synthetic magnetic field profiles as would be seen by a spacecraft traversing the cloud are presented.

Description: Since radio propagation measurements using either natural or spacecraft radio signals are used for probing the solar wind in the vicinity of the sun, they represent a key tool for studying the interplanetary consequences of solar structure and dynamic phenomena. New information on the near sun consequences was obtained from radio scintillation observations of coherent spacecraft signals. The results covering density fluctuations, fractional density fluctuations, coronal streamers, heliospheric current sheets, coronal mass ejections and interplanetary shocks are reviewed. A joint ICE S-band (13 cm wavelength) Doppler scintillation measurement with the SOHO white-light coronograph (LASCO) is described.

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: The knowledge of vegetation dielectric behavior is important in studying the scattering properties of the vegetation canopy and radar backscatter modelling. Until now, a limited number of studies have been published on the dielectric properties in the boreal forest context. This paper presents the results of the dielectric constant as a function of depth in the trunks of two common boreal forest species: black spruce and trembling aspen, obtained from field measurements. The microwave penetration depth for the two species is estimated at C, L, and P bands and used to derive the equivalent dielectric constant for the trunk as a whole. The backscatter modelling is carried out in the case of black spruce and the results are compared with the JPL AIRSAR data. The sensitivity of the backscatter coefficient to the dielectric constant is also examined.

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: The Galileo Ultravilet Spectrometer Experiment (UVS) obtained a partial celestial sphere map of interplanetary Lyman-alpha (IP L alpha) on 13-14 December 1990 during the first Earth encounter. The Galileo spacecraft was near the downwind axis of the local interstellar medium flow. These UVS measurements sampled the downwind, anti-sunward hemisphere. The data were modeled using a hot model of the interplanetary hydrogen density distribution with the goal of studying multiple scattering effects in the inner solar system. The derived ratio in the downwind direction of the observed brightness and a single scattering model brightness, both normalized to unity in the upwind direction, is 1.82 +/- 0.2. This brightness ratio requires a multiple scattering correction which is 36% larger than can be accounted for by theoretical calculations. The hot model may require: (1) a temperature perturbation of the interstellar wind velocity distribution or (2) an additional downstream source of interplanetary hydrogen. However, a more likely exlanation which affects the hot model is the latitude dependence of the radiation pressure. This dependence, based on the known solar L alpha flux latitude variation at solar maximum, causes a downwind brightness enhancement by preferential focusing of H-atoms with trajectory planes containing the solar poles. This result implies that radiation pressure near the solar poles is nearly independent of solar cycle and is insufficient to lead to a net repulsion of hydrogen atoms by the sun, as can occur near the ecliptic plane during the solar maximum. In addition, the UVS performed 13 observations of IP L alpha while in cruise between Venus and the Earth in 3 directions fixed in ecliptic coordinates.

Description: In this paper we generalize earlier gasdynamic analyses of the motion of the heliospheric termination shock in response to upstream disturbances (Barnes, 1993, 1994; Naidu and Barnes, 1994), to include magnetohydrodynamic (MHD) phenomena. We assume that the termination shock is a strong, perpendicular shock and that the initial upstream disturbance is a tangential discontinuity. The resulting configuration after the interaction is very similar to that in the gasdynamic models after an interaction with a contact discontinuity or interplanetary shock, and for an increase (decrease) in dynamic pressure consists of an outward (inward) propagating termination shock and an outward propagating shock (MHD rarefraction wave) that carries the signal of the disturbance into the far downstream plasma. The plasma immediately behind the new termination shock is separated from the downstream signal by a tangential discontinuity. The results of the model show that the speed of the new termination shock depends mainly on the magnitude of the change in dynamic pressure and are typically of order approximately 100 km/s, comparable to the results of the gasdynamic models.

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 term 'magnetic hole' has been used to denote isolated intervals when the magnitude of the interplanetary magnetic field drops to a few tenths, or less, of its ambient value for a time that corresponds to a linear dimension of tens to a few hundreds of proton gyro-radii. Data obtained by the Ulysses magnetometer and solar wind anlayzer have been combined to study the properties of such magnetic holes in the solar wind between 1 AU and 5.4 AU and to 23 deg south latitude. In order to avoid confusion with decreases in field strength at interplanetary discontinuities, the study has focused on linear holes across which the field direction changed by less than 5 deg. The holes occurred preferentially, but not without exception, in the interaction regions on the leading edges of high-speed solar wind streams. Although the plasma surrounding the holes was generally stable against the mirror instability, there are indications that the holes may have been remnants of mirror-mode structures created upstream of the points of observation. Those indications include the following: (1) For the few holes for which proton of alpha-particle pressure could be measured inside the hole, the ion thermal pressure was always greater than in the plasma adjacent to the holes. (2) The plasma surrounding many of the holes was marginally stable for the mirror mode, while the plasma environment of all holes was significantly closer to mirror instability than was the average solar wind. (3) The plasma containing trains of closely spaced holes was closer to mirror instability than was the plasma containing isolated holes. (4) The near-hole plasma had much higher ion beta (ratio of thermal to magnetic pressure) than did the average solar wind. (5) Near the holes, T(sub perp)/T(sub parallel) tended to be either greater than 1 or larger than in the average wind. (6) The proton and alpha-particle distribution functions measured inside the holes occasionally exhibited the flattened phase-space-density contoures in nu(sub perp)/nu(sub parallel) space found in some numerical simulations of the mirror instability.

Description: During the international campaign of June 1991, the active region AR 6659 produced six very large, long-duration flares (X10/12) during its passage across the solar disk. We present the characteristics of four of them (June 4, 6, 9, 15). Precise measurements of the spot motions from Debrecen and Tokyo white-light pictures are used to understand the fragmentation of the main sunspot group with time. This fragmentation leads to a continuous restructuring of the magnetic field pattern while rapid changes are evidenced due to fast new flux emergence (magnetograms of Marshall Space Flight Center (MSFC), Huairou). The first process leads to a shearing of the field lines along which there is energy storage; the second one is the trigger which causes the release of energy by creating a complex topology. We conjecture that these two processes with different time scales are relevant to the production of flares.

Description: Because changes in the Earth's environment have become major global issues, continuous, longterm scientific information is required to assess global problems such as deforestation, desertification, greenhouse effects and climate variations. Global change studies require understanding of interactions of complex processes regulating the Earth system. Space-based Earth observation is an essential element in global change research for documenting changes in Earth environment. It provides synoptic data for conceptual predictive modeling of future environmental change. This paper provides a brief overview of remote sensing technology from the perspective of global change research.

Description: Interplanetary disturbances characterized by plasma that is more turbulence and/or moves faster than the background solar wind are readily defected as transients in Doppler scintillation measurements of the near-Sun solar wind. Systematic analysis of over 23,000 hours of Pioneer Venus Orbiter Doppler measurements obtained inside 0.5 AU during 1979-1987 have made it possible for the first time to investigate the frequency of occurrence of Doppler scintillation transients under solar minimum conditions and to determine its dependence on solar cycle. On the basis of a total of 142 transients, Doppler scintillation transient rates vary from a high of 0.22 in 1979 (one every 4.6 days) to a low of 0.077 transients/d in 1986 (one every 13 days), a decrease by almost a factor of 3 from solar maximum to solar minimum. This solar cycle variation, the strongest yet of any solar wind Doppler scintillation property, is highly correlated with both solar activity characterized by sunspot number and the coronal mass ejection rates deduced from Solswind and Solar Maximum Mission (SMM) coronagraph observations. These results indicate that coronal mass ejections and Doppler scintillation transients are closely related not just during solar maximum, as occasional individual comparisons have shown in the past, but throughout the entire solar cycle, and strengthen the notation that the Doppler scintillation and optical transients are different manifestations of the same physical phenomenon. The magnitudes of the transients, as described by the ratio of peak to pretransient scintillation levels (EF for enhancement factor), and their distribution iwth heliocentric distance also vary with solar cycle. While EF tends to diminish with increasing heliocentric distance during high solar activity, it is more evenly distributed during low solar activity. EF is also lower during solar minimum, as 13% of the transients during solar maximum have values exceeding 23, the highest EF observed during solar minimum. These results are consistent with the fact that occasional major fast-moving interplanetary shocks that are observed during solar maximum are very rate during solar minimum.

Description: It has recently been reported that the total radiative emission variations from solar type stars exceeds the currently solar constant variations (from spacecraft over the last decade) by a factor near 4. Aside from other remote alternatives, this suggests three clear possibilities: (1) the Sun may undergo irradiance variations several times larger than any we have seen; (2) our Sun is highly unusual with regard to its radiative output; or (3) our terrestrial position in the heliosphere provides a special vantage point which reduces the observed solar irradiance variations. We investigate the last possibility by considering the influence of observer latitude upon calculated irradiance variations using a simple model for emission from solar contrast features. We consider modeled sunspots, faculae, and network structures. As the latitude angle of the observer rises relative to the heliographic equator, sunspot deficit contributions diminish and facular plus network contributions escalate. We find that the observing latitude can influence the irradiance variations by a factor near 6. When we integrate the irradiance variations, over the celestial sphere, they average to 3 times the terrestrial effect, suggesting that the solar cycle luminosity variations are proportionally, 3 times larger than the solar constant variations. Thus we suggest the Sun's luminosity output varies even more strongly with the solar cycle than is apparent in the solar constant variations. The influence of the observer viewing angle relative to stellar spin axis, studied here, may be possible to investigate with a thorough statistical examination of other type stars. Additionally, the rotational modulation due to active regions (as a function of observer viewing angle) may also be a valuable are for future investigation.

Description: For more than 90 years, solar extreme ultraviolet (EUV) irradiance modeling has progressed from empirical blackbody radiation formulations, through fudge factors, to typically measured irradiances and reference spectra was well as time-dependent empirical models representing continua and line emissions. A summary of recent EUV measurements by five rockets and three satellites during the 1980s is presented along with the major modeling efforts. The most significant reference spectra are reviewed and threee independently derived empirical models are described. These include Hinteregger's 1981 SERF1, Nusinov's 1984 two-component, and Tobiska's 1990/1991/SERF2/EUV91 flux models. They each provide daily full-disk broad spectrum flux values from 2 to 105 nm at 1 AU. All the models depend to one degree or another on the long time series of the Atmosphere Explorer E (AE-E) EUV database. Each model uses ground- and/or space-based proxies to create emissions from solar atmospheric regions. Future challenges in EUV modeling are summarized including the basic requirements of models, the task of incorporating new observations and theory into the models, the task of comparing models with solar-terrestrial data sets, and long-term goals and modeling objectives. By the late 1990s, empirical models will potentially be improved through the use of proposed solar EUV irradiance measurements and images at selected wavelengths that will greatly enhance modeling and predictive capabilities.

Description: Enhanced fluxes of suprathermal electrons are commonly observed upstream of corotating forward and reverse shocks in the solar wind at heliocentric distances beyond approximately 2 AU by the Los Alamos plasma experiment on Ulysses. The average duration of these events, which are most intense immediately upstream from the shocks and which fade with increasing distance from them, is approximately 2.4 days near 5 AU. These events are caused by the leakage of shock-heated electrons into the upstream region. The upstream regions of these shocks face back toward the Sun along the interplanetary magnetic field, so these leaked electrons commonly counterstream relative to the normal solar wind electron heat flux. The observations suggest that conservation of magnetic moment and scattering typically limit the sunward propagation of these electrons as beams to field-aligned distances of approximately 15 AU. Although it seems unlikely that these shock-associated events are an important source of counterstreaming events near 1 AU, remnants of the backstreaming beams may contribute importantly to the diffuse solar wind halo electron population there.

Description: In May, 1993, the heliospheric current sheet (HCS) ceased to be seen by the Ulysses spacecraft at a heliocentric latitude of approximately 30 deg S and distance of 4.7 AU. The disappearance of the HCS coincided with the solar wind speed remaining greater than 560 km/s and with the disappearance of one of four interaction regions previously seen on each solar rotation. The heliographic latitude of the disappearance of the HCS at Ulysses was 11 deg equatorward of the latitude of the magnetic neutral sheet computed at the source surface at 2.5 solar radii, and it occurred a half year earlier than predicted on the basis of the persistance of the time profile of the neutral sheet tilt from one solar cycle to the next.

Description: The unique vantage point of the Ulysses spacecraft throughout 1992 and the beginning of 1993, at a close to constant heliocentric distance of about 5 AU and a slowly varying heliographic latitude from 5 deg to 30 deg south is used to describe and discuss the evolution of the sector structure of the interplanetary magnetic field during the declining phase of the solar cycle. From the end of 1990 to the beginning of 1992 the sector structure changed from a four sector to a two sector structure, but remained constant in solar longitude. From about June-July 1992, the structure, matching the evolution in the computed coronal magnetic fields, drifted eastwards, with a recurrence period of about 28 days. This result may indicate a slower rotation rate for the dipolar component of the solar magnetic field which becomes dominant about this time in the solar cycle.

Description: Near-ecliptic solar wind observations by Ulysses on its way to the polar regions of the Sun, compared with those from IMP 8 at 1 AU, showed that high-speed streams decay and broaden with heliocentric distance from IMP 8 to Ulysses, as expected. In July 1992 while traveling south at approximately 13 deg S and 5.3 AU, Ulysses encountered a recurrent high-speed stream, that may also have been observed at IMP 8. The stream has been observed a total of 14 times, once in each solar rotation through June 1993 at approximately 34 deg S. The source of the high-speed stream is an equatorward extension of the south polar coronal hole. From July 1992 through June 1993, averages of solar wind peak speed increased while density decreased with heliographic latitude. Both the stream and a low-speed, high-density flow, presumably associated with the heliomagnetic (coronal) streamer belt encircling the heliomagnetic equator, crossed Ulysses with the solar rotation period until April 1993 when the spacecraft was at approximately 29 deg S heliographic latitude. After this time, as the spacecraft climbed to higher latitudes, the central portion of the streamer belt with lowest speed and highest density disappeared. Therefore, at its maximum inclination, the belt was tilted at approximately 29 deg to the heliographic equator at this point in the solar cycle.

Description: During the late declining phase of the solar cycle, the tilt of the solar magnetic dipole with respect to the Sun's rotation axis leads to large-scale organization of the solar wind, such that alternating regions of high- and low-speed solar wind are observed in the ecliptic plane. In this paper, we use Doppler scintillation measurements to investigate mass flux of these two types of solar wind in the ecliptic plane and inside 0.3 AU, where in situ measurements have not been possible. To the extent that Doppler scintillation reflects mass flux, we find that mass flux in high-speed streams: (1) is lower (by a factor of approximately 2.2) than the mass flux of the average solar wind in the heliocentric distance range of 0.3-0.5 AU; (2) is lower still (by as much as a factor of about 4) than the mass flux of the slow solar wind associated with the streamer belt; and (3) appears to grow with heliocentric distance. These Doppler scintillation results are consistent with the equator to pole decrease in mass flux observed in earlier spectral broadening measurements, and with trends and differences between high- and low-speed solar wind observed by in situ measurements in the range of 0.3-0.1 AU. The mass flux results suggest that the solar wind flow in high-speed streams is convergent towards the ecliptic near the Sun, becoming less convergent and approaching radial with increasing heliocentric distance beyond 0.3 AU. The variability of mass flux observed within equatorial and polar high-speed streams close to the Sun is strikingly low. This low variability implies that, as Ulysses currently ascends to higher latitudes and spends more time in the south polar high-speed stream after crossing the heliocentric current sheet, it can expect to observe a marked decrease in variations of both mass flux and solar wind speed, a trend that appears to have started already.

Description: Highly structured Langmuir waves, also known as electron plasma oscillations, have been observed in the foreshock of Venus using the plasma wave experiment on the Galileo spacecraft during the gravity assist flyby on February 10, 1990. The Galileo wideband sampling system provides digital electric field waveform measurements at sampling rates up to 201,600 samples per second, much higher than any previous instrument of this type. The main Langmuir wave emission band occurs near the local electron plasma frequency, which was approximately 43 kHz. The Langmuir waves are observed to shift above and below the plasma frequency, sometimes by as much as 20 kHz. The shifts in frequency are closely correlated with the downstream distance from the tangent field line, implying that the shifts are controlled by the electron beam velocity. Considerable fine structure is also evident, with timescales as short as 0.15 ms, corresponding to spatial scales of a few tens of Debye lengths. The frequency spectrum often consists of beat-type waveforms, with beat frequencies ranging from 0.2 to 7 kHz, and in a few cases, isolated wave packets. The peak electric field strengths are approximately 1 mV/m. These field strengths are too small for strongly nonlinear processes to be important. The beat-type waveforms are suggestive of a parametric decay process.

Description: Four different arraying schemes applicable to deep space communications are discussed and analyzed. These include symbol stream combining (SSC), baseband combining (BC), carrier arraying (CA) and full spectrum combining (FSC). Complexity versus performance is traded off throughout the paper and benefits to the reception of existing spacecraft signals are discussed.

Description: The scattering properties of perfectly conducting and resistive strips are predicted for strips which are located on a dielectric slab backed by a perfectly conducting ground plane. The spectral domain Green's function is used to relate the currents and fields on the strip, and the resulting integral equation is solved using the method of moments. Both TE and TM strips are examined using piecewise linear and pulse subdomain basis functions, respectively, to model the current on the strip. Calculated results are compared with results measured at the NASA Langley Research Center.

Description: The irregular polygonal pattern of solar granulation is analyzed for size-shape relations using statistical crystallography. In contrast to previous work which has assumed perfectly hexagonal patterns for granulation, more realistic accounting of cell (granule) shapes reveals a broader basis for quantitative analysis. Several features emerge as noteworthy: (1) a linear correlation between number of cell-sides and neighboring shapes (called Aboav-Weaire's law); (2) a linear correlation between both average cell area and perimeter and the number of cell-sides (called Lewis's law and a perimeter law, respectively) and (3) a linear correlation between cell area and squared perimeter (called convolution index). This statistical picture of granulation is consistent with a finding of no correlation in cell shapes beyond nearest neighbors. A comparative calculation between existing model predictions taken from luminosity data and the present analysis shows substantial agreements for cell-size distributions. A model for understanding grain lifetimes is proposed which links convective times to cell shape using crystallographic results.

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: Dual-frequency ranging and Doppler measurements were conducted in support of the Ulysses Solar Corona Experiment (SCE) at and around the spacecraft's first solar conjunction in 1991 August. The differential group delay time between range codes on the two downlink carrier signals at the wavelengths 13.1 and 3.6 cm, a direct measure of the total electron content between spacecraft and ground station, was used to derive the electron density distribution in the solar corona. Linear power-law representations of the coronal electron density were derived for the range of solar distances from 4 solar radii to 40 solar radii on both sides of the Sun. The corona was found to be very nearly symmetric; the radial falloff exponent being 2.54 +/- 0.05 for occultation ingress (east solar limb) and 2.42 +/- 0.05 for egress (west limb), respectively. The departure of these exponents from the inverse equare relation implies that significant solar wind acceleration is occurring within the radial range of the observations. The electron density level was found to be considerably lower than that observed during the 1988 December solar occultation of Voyager 2. Although the smoothed sunspot number R(sub z) (a standard indicator of solar activity) was almost the same in 1988 December and 1991 August, the mean electron density at 20 solar radii was found to be 1.7 +/- 0.1 x 10(exp 3)/cu cm during the Ulysses conjunction, a decline by almost a factor of 4 from the value obtained during the Voyager conjunction.

Description: One of the long-standing uncertainties in the wave-resonance theory of coronal heating is the stability of the resonance layer. The wave motions in the resonance layer produce highly localized shear flows which vary sinusoidally in time with the resonance period. This configuration is potentially susceptible to the Kelvin-Helmholtz instability (KHI), which can enhance small-scale structure and turbulent broadening of shear layers on relatively rapid ideal timescales. We have investigated numerically the response of a characteristic velocity profile, derived from resonance absorption models, to finite fluid perturbations comparable to photospheric fluctuations. We find that the KHI primarily should affect long (approximately greater than 6 x 10(exp 4) km) loops where higher velocity flows (M approximately greater than 0.2) exist in resonance layers of order 100 km wide. There, the Kelvin-Helmholtz growth time is comparable to or less than the resonance quarter-period, and the potentially stabilizing magnetic effects are not felt until the instability is well past the linear growth stage. Not only is the resonance layer broadened by the KHI, but also the convective energy transport out of the resonance layer is increased, thus adding to the efficiency of the wave-resonance heating process. In shorter loops, e.g., those in bright points and compact flares, the stabilization due to the magnetic field and the high resonance frequency inhibit the growth of the Kelvin-Helmholtz instability beyond a minimal level.

Description: The solar flare and coronal mass ejection (CME) events associated with the large and complex March 1989 active region are discussed. This active region gave us a chance to study the relation of CME with truly major solar flares. The work concentrates on questions of the relation of CMEs and flares to one another and to other types of activity on the Sun. As expected, some major (X-3B class) flares had associated CMEs. However, an unexpected finding is that others did not. In fact, there is strong evidence that the X4-4B flare of March 9th had no CME. This lack of a CME for such an outstanding flare event has important implications to theories of CME causation.Apparently, not all major flares cause CMEs or are caused by CMEs. The relations between CMEs and other types of solar activity are also discussed. No filament disappearances are reported for major CMEs studied here. Comparing these results with other studies, CMEs occur in association with flares and with erupting prominences, but neither are required for a CME. The relation between solar structures showing flaring without filament eruptions and structures showing filament eruptions without flares becomes important. The evolutionary relation between an active flaring sunspot region and extensive filaments without sunspots is reviewed, and the concept of an 'evolving magnetic structure' (EMS) is introduced. It is suggested that all CMEs arise in EMSs and that CMEs provide a major path through which azimuthal magnetic fields escape form the Sun during the solar cycle.

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: This paper presents the results of the correlation analysis of the Skylab S-193 13.9 GHz Radiometer/Scatterometer data. Computer analysis of the S-193 data shows more than 50 percent of the radiometer and scatterometer data are uncorrelated. The correlation coefficients computed for the data gathered over various ground scenes indicates the desirability of using both active and passive sensors for the determination of various Earth phenomena.

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: The large gradual solar-energetic-particle (SEP) events, where abundances are commonly measured, are produced when coronal mass ejections (CMEs) drive shock waves through the corona and the interplanetary medium. The shock accelerates particles from the highly-ionized, approximately 1.5 MK, plasma in a manner that depends only weakly upon the Q/A of the ion, except at very high energies. Averaging the approximately 1 MeV/amu abundances over many events compensates for the acceleration effects to produce abundances that appear to correspond directly to those in the coronal source for all observed elements, including H. The resulting abundances reflect the 4 x enhancement of ions with low values of first ionization potential (FIP) arising from ion-neutral fractionation that occurs as the atoms are transported up from the photosphere. A different pattern of fractionation is found for ions that are shock-accelerated from the high speed solar wind emerging from coronal holes.

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: POMESH is a computer program capable of predicting the performance of reflector antennas. Both far field pattern and gain calculations are performed using the Physical Optics (PO) approximation of the equivalent surface currents. POMESH is primarily intended for relatively small reflectors. It is useful in situations where the surface is described by irregular data that must be interpolated and for cases where the surface derivatives are not known. This method is flexible and robust and also supports near field calculations. Because of the near field computation ability, this computational engine is quite useful for subreflector computations. The program is constructed in a highly modular form so that it may be readily adapted to perform tasks other than the one that is explicitly described here. Since the computationally intensive portions of the algorithm are simple loops, the program can be easily adapted to take advantage of vector processor and parallel architectures. In POMESH the reflector is represented as a piecewise planar surface comprised of triangular regions known as facets. A uniform physical optics (PO) current is assumed to exist on each triangular facet. Then, the PO integral on a facet is approximated by the product of the PO current value at the center and the area of the triangle. In this way, the PO integral over the reflector surface is reduced to a summation of the contribution from each triangular facet. The source horn, or feed, that illuminates the subreflector is approximated by a linear combination of plane patterns. POMESH contains three polarization pattern definitions for the feed; a linear x-polarized element, linear y-polarized element, and a circular polarized element. If a more general feed pattern is required, it is a simple matter to replace the subroutine that implements the pattern definitions. POMESH obtains information necessary to specify the coordinate systems, location of other data files, and parameters of the desired calculation from a user provided data file. A numerical description of the principle plane patterns of the source horn must also be provided. The program is supplied with an analytically defined parabolic reflector surface. However, it is a simple matter to replace it with a user defined reflector surface. Output is given in the form of a data stream to the terminal; a summary of the parameters used in the computation and some sample results in a file; and a data file of the results of the pattern calculations suitable for plotting. POMESH is written in FORTRAN 77 for execution on CRAY series computers running UNICOS. With minor modifications, it has also been successfully implemented on a Sun4 series computer running SunOS, a DEC VAX series computer running VMS, and an IBM PC series computer running OS/2. It requires 2.5Mb of RAM under SunOS 4.1.1, 2.5Mb of RAM under VMS 5-4.3, and 2.5Mb of RAM under OS/2. The OS/2 version requires the Lahey F77L compiler. The standard distribution medium for this program is one 5.25 inch 360K MS-DOS format diskette. It is also available on a .25 inch streaming magnetic tape cartridge in UNIX tar format and a 9-track 1600 BPI magnetic tape in DEC VAX FILES-11 format. POMESH was developed in 1989 and is a copyrighted work with all copyright vested in NASA. CRAY and UNICOS are registered trademarks of Cray Research, Inc. SunOS and Sun4 are trademarks of Sun Microsystems, Inc. DEC, DEC FILES-11, VAX and VMS are trademarks of Digital Equipment Corporation. IBM PC and OS/2 are registered trademarks of International Business Machines, Inc. UNIX is a registered trademark of Bell Laboratories.

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: The Interactive Digital Signal Processor, IDSP, consists of a set of time series analysis "operators" based on the various algorithms commonly used for digital signal analysis work. The processing of a digital time series to extract information is usually achieved by the application of a number of fairly standard operations. However, it is often desirable to "experiment" with various operations and combinations of operations to explore their effect on the results. IDSP is designed to provide an interactive and easy-to-use system for this type of digital time series analysis. The IDSP operators can be applied in any sensible order (even recursively), and can be applied to single time series or to simultaneous time series. IDSP is being used extensively to process data obtained from scientific instruments onboard spacecraft. It is also an excellent teaching tool for demonstrating the application of time series operators to artificially-generated signals. IDSP currently includes over 43 standard operators. Processing operators provide for Fourier transformation operations, design and application of digital filters, and Eigenvalue analysis. Additional support operators provide for data editing, display of information, graphical output, and batch operation. User-developed operators can be easily interfaced with the system to provide for expansion and experimentation. Each operator application generates one or more output files from an input file. The processing of a file can involve many operators in a complex application. IDSP maintains historical information as an integral part of each file so that the user can display the operator history of the file at any time during an interactive analysis. IDSP is written in VAX FORTRAN 77 for interactive or batch execution and has been implemented on a DEC VAX-11/780 operating under VMS. The IDSP system generates graphics output for a variety of graphics systems. The program requires the use of Versaplot and Template plotting routines and IMSL Math/Library routines. These software packages are not included in IDSP. The virtual memory requirement for the program is approximately 2.36 MB. The IDSP system was developed in 1982 and was last updated in 1986. Versaplot is a registered trademark of Versatec Inc. Template is a registered trademark of Template Graphics Software Inc. IMSL Math/Library is a registered trademark of IMSL Inc.

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: This collection of programs comprises The UNIX Based Data Management System for the Pilot Field Experiment (PiFEx) which is an attempt to mimic the Mobile Satellite (MSAT) scenario. The major purposes of PiFEx are to define the mobile communications channels and test the workability of new concepts used to design various components of the receiver system. The results of the PiFex experiment are large amounts of raw data which must be accessed according to a researcher's needs. This package provides a system to manage the PiFEx data in an interactive way. The system not only provides the file handling necessary to retrieve the desired data, but also several FORTRAN programs to generate some standard results pertaining to propagation data. This package assumes that the data file initially generated from the experiment has been already converted from binary to ASCII format. The Data Management system described here consists of programs divided into two categories: those programs that handle the PiFEx generated files and those that are used for number-crunching of these files. Five FORTRAN programs and one UNIX shell script file are used for file manipulation purposes. These activities include: calibration of the acquired data; and parsing of the large data file into datasets concerned with different aspects of the experiment such as the specific calibrated propagation data, dynamic and static loop error data, statistical data, and temperature and spatial data on the hardware used in the experiment. The five remaining FORTRAN programs are used to generate usable information about the data. Signal level probability, probability density of the signal fitting the Rician density function, frequency of the data's fade duration, and the Fourier transform of the data can all be generated from these data manipulation programs. In addition, a program is provided which generates a downloadable file from the signal levels and signal phases files for use with the plotting routine AKPLOT (NPO-16931). All programs in this package are written in either FORTRAN-77 or UNIX shell-scripts. The package does not include test data. The programs were developed in 1987 for use with a UNIX operating system on a DEC MicroVAX computer.

Description: The Optical Communication Link Analysis Program, OPTI, analyzes optical and near-infrared communication links that use pulse position modulation (PPM) and direct detection. The program prompts for inputs of system component parameters, modulation format and other operational parameters, background noise sources, and the desired link bit error rate. From these inputs, link margin is determined and a link design control table (DCT) is generated. The program also allows the user to save sets of input parameters defining a given link and read them back into the program at a later time. Further, the program has the capability of altering automatically any of the input parameters to achieve a desired link margin. The program provides a table of extended background sources, e.g. planets, the moon, and the sun. To compute background noise, only the distance from the receiver to the noise source(s) must be entered. This determines whether or not the whole object is in the field-of-view. The program assumes that each object is a blackbody (at 5900K) with an overall visible magnitude scaled to match the tabulated data. Also provided is a table of 19 bright stars. If the noise source is one or several of these, then only the name(s) of the star(s) are required. Noise sources that are not among those contained in the program can be entered as "additional noise sources". In this case required information includes whether or not the source is a point or extended source, its radiance (for extended sources) or irradiance, and receiver to source distance in A.U. (for point sources). The OPTI program is written in FORTRAN-77 and was designed to be used on the IBM PC and PC/AT personal computers. (Note: The 8087/80287 math coprocessor option is highly recommended for use with this program.) The program will also compile under UNIX 4.3 BSD FORTRAN-77 with minor changes. OPTI was developed in 1987.

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 "Data Management System for Mobile Satellite Propogation" package is a collection of FORTRAN programs and UNIX shell scripts designed to handle the huge amounts of data resulting from Mobile Satellite propogation experiments. These experiments are designed to assist in defining channels for mobile satellite systems. By understanding multipath fading characteristics of the channel, doppler effects, and blockage due to manmade objects as well as natural surroundings, characterization of the channel can be realized. Propogation experiments, then, are performed using a prototype of the system simulating the ultimate product environment. After the data from these experiments is generated, the researcher must access this data with a minimum of effort and to derive some standard results. The programs included in this package manipulate the data files generated by the NASA/JPL Mobile Satellite propogation experiment on an interactive basis. In the experiment, a transmitter operating at 869 MHz was carried to an altitude of 32Km by a stratospheric balloon. A vehicle within the line-of-sight of the transmitter was then driven around, splitting the incoming signal into I and Q channels, and sampling the resulting signal strength at 1000 samples per second. The data was collected at various antenna elavation angles and different times of day generating the ancillary data for the experiment. This package contains a program to convert the binary format of the data generated into standard ASCII format suitable for use with a wide variety of machine architectures. Also included is a UNIX shell-script designed to parse this ASCII file into those records of data that match the researcher's desired values for the ancillary data parameters. In addition, four FORTRAN programs are included to obtain standard quantities from the data. Quantities such as probability of signal level greater than or equal to a specified signal level, probability density of the signal levels, frequency of fade duration, and Fourier Transforms of the sampled data can be generated from the propogation experiment data. All programs in this package are written in either FORTRAN 77 or UNIX shell-scripts. The package does not include test data. The programs were developed in 1987 for use with a UNIX operating system on a DEC MicroVAX computer.

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 Distributed Architecture for Phased Array Antennas program, DISTAR, is a simulation tool used to study the implementation of distributed phased array hardware. DISTAR allows the placement of possibly faulty transmit/receive modules (T/R) at locations throughout the array. Variations in amplifiers and phase shifters may degrade antenna performance, depending on environmental conditions and array architecture. DISTAR enables the antenna designer to examine array characteristics and how they effect both type and extent of antenna failures. General specifications for amplifier and phase shifter tolerances may also be determined for various architectures. The array geometry is assumed to be rectangular with optional staggering in the X or Y direction. The array is divided into various subgroups both geometrically and electronically. Each T/R drives a subgroup of elemental radiators. Both hard and soft failures of the T/R amplifiers are modeled. Hard failures are catastrophic as no power is transmitted to those elements controlled by that particular amplifier. Soft, or non-catastrophic, failures are simulated by a modified Gaussian distribution. Flat, warped, and parabolic arrays can be modeled. The user must supply array geometry, amplitude and phase for each row and column, and location of all T/Rs. Specific or randomly chosen amplifiers may be deemed faulty. DISTAR outputs the normalized antenna gain patterns in the form of tables, and 2-D or 3-D graphs. DISTAR is written in FORTRAN 77 for interactive execution and has been implemented on a DEC VAX 11/780 computer operating under VMS 4.4 with a central memory requirement of approximately 35K of 8 bit bytes (without graphics). DISTAR requires the commercial DISSPLA plot library for full graphics implementation. This program was developed in 1986.

Description: RHCP, the Right-Handed, Circularly Polarized Microstrip Antenna program, aids in the design of a rectangular antenna element, given the desired frequency of operation and substrate characteristics. RHCP begins the design calculations based on a square element with linear polarization. The effective dielectric constant and changes in electrical length due to fringing at the radiating element edges are taken into account. A coaxial feed is inset with 50 ohms input impedance. By placing the feed such that two orthonormal modes are produced in the antenna cavity, right- or left-handed circular polarization is obtained. Input to RHCP consists of desired frequency, dielectric constant, and substrate thickness. Output consists of the final rectangular geometry, the proposed feed inset placement, and actual input impedance. RHCP has been used successfully for frequencies between 2 and 15 GHz for thin substrates. This program was used to fabricate antenna elements for the S-band quad antennas on board the Space Shuttle, and is a part of the design project for the S-band phased array antenna radiating aperture. RHCP is written in FORTRAN 77 for interactive execution and has been implemented on a DEC VAX series computer operating under VMS. This program was developed in 1985.

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: MOM3D (LAR-15074) is a FORTRAN method-of-moments electromagnetic analysis algorithm for open or closed 3-D perfectly conducting or resistive surfaces. Radar cross section with plane wave illumination is the prime analysis emphasis; however, provision is also included for local port excitation for computing antenna gain patterns and input impedances. The Electric Field Integral Equation form of Maxwell's equations is solved using local triangle couple basis and testing functions with a resultant system impedance matrix. The analysis emphasis is not only for routine RCS pattern predictions, but also for phenomenological diagnostics: bistatic imaging, currents, and near scattered/total electric fields. The images, currents, and near fields are output in form suitable for animation. MOM3D computes the full backscatter and bistatic radar cross section polarization scattering matrix (amplitude and phase), body currents and near scattered and total fields for plane wave illumination. MOM3D also incorporates a new bistatic k space imaging algorithm for computing down range and down/cross range diagnostic images using only one matrix inversion. MOM3D has been made memory and cpu time efficient by using symmetric matrices, symmetric geometry, and partitioned fixed and variable geometries suitable for design iteration studies. MOM3D may be run interactively or in batch mode on 486 IBM PCs and compatibles, UNIX workstations or larger computers. A 486 PC with 16 megabytes of memory has the potential to solve a 30 square wavelength (containing 3000 unknowns) symmetric configuration. Geometries are described using a triangular mesh input in the form of a list of spatial vertex points and a triangle join connection list. The EM-ANIMATE (LAR-15075) program is a specialized visualization program that displays and animates the near-field and surface-current solutions obtained from an electromagnetics program, in particular, that from MOM3D. The EM-ANIMATE program is windows based and contains a user-friendly, graphical interface for setting viewing options, case selection, file manipulation, etc. EM-ANIMATE displays the field and surface-current magnitude as smooth shaded color fields (color contours) ranging from a minimum contour value to a maximum contour value for the fields and surface currents. The program can display either the total electric field or the scattered electric field in either time-harmonic animation mode or in the root mean square (RMS) average mode. The default setting is initially set to the minimum and maximum values within the field and surface current data and can be optionally set by the user. The field and surface-current value are animated by calculating and viewing the solution at user selectable radian time increments between 0 and 2pi. The surface currents can also be displayed in either time-harmonic animation mode or in RMS average mode. In RMS mode, the color contours do not vary with time, but show the constant time averaged field and surface-current magnitude solution. The electric field and surface-current directions can be displayed as scaled vector arrows which have a length proportional to the magnitude at each field grid point or surface node point. These vector properties can be viewed separately or concurrently with the field or surface-current magnitudes. Animation speed is improved by turning off the display of the vector arrows. In RMS modes, the direction vectors are still displayed as varying with time since the time averaged direction vectors would be zero length vectors. Other surface properties can optionally be viewed. These include the surface grid, the resistance value assigned to each element of the grid, and the power dissipation of each element which has an assigned resistance value. The EM-ANIMATE program will accept up to 10 different surface current cases each consisting of up to 20,000 node points and 10,000 triangle definitions and will animate one of these cases. The capability is used to compare surface-current distribution due to various initial excitation directions or electric field orientations. The program can accept up to 50 planes of field data consisting of a grid of 100 by 100 field points. These planes of data are user selectable and can be viewed individually or concurrently. With these preset limits, the program requires 55 megabytes of core memory to run. These limits can be changed in the header files to accommodate the available core memory of an individual workstation. An estimate of memory required can be made as follows: approximate memory in bytes equals (number of nodes times number of surfaces times 14 variables times bytes per word, typically 4 bytes per floating point) plus (number of field planes times number of nodes per plane times 21 variables times bytes per word). This gives the approximate memory size required to store the field and surface-current data. The total memory size is approximately 400,000 bytes plus the data memory size. The animation calculations are performed in real time at any user set time step. For Silicon Graphics Workstations that have multiple processors, this program has been optimized to perform these calculations on multiple processors to increase animation rates. The optimized program uses the SGI PFA (Power FORTRAN Accelerator) library. On single processor machines, the parallelization directives are seen as comments to the program and will have no effect on compilation or execution. MOM3D and EM-ANIMATE are written in FORTRAN 77 for interactive or batch execution on SGI series computers running IRIX 3.0 or later. The RAM requirements for these programs vary with the size of the problem being solved. A minimum of 30Mb of RAM is required for execution of EM-ANIMATE; however, the code may be modified to accommodate the available memory of an individual workstation. For EM-ANIMATE, twenty-four bit, double-buffered color capability is suggested, but not required. Sample executables and sample input and output files are provided. Electronic documentation is provided for both EM-ANIMATE and MOM3D in PostScript format. Documentation for EM-ANIMATE is also provided in the form of IRIX man pages. The standard distribution medium for COS-10048 is a .25 inch streaming magnetic IRIX tape cartridge in UNIX tar format. MOM3D and EM-ANIMATE are also available separately as LAR-15074 and LAR-15075, respectively. MOM3D was developed in 1992. EM-ANIMATE was developed in 1993.

Description: The signal from a stable periodicity can seem to be intermittent when it is partially masked by an unmodelled window function or when the data set is too short to resolve closely spaced periodicities. By taking this into account, short-lived periodicities in solar data can be reinterpreted as evidence for continuously periodic behavior. The periodic sources are located in the solar interior and caused by global oscillation modes. The convective envelope acts as the window for these sources. Recent reports of seven periodicities from 100 to 1000 days are compared with this model. Precise long-term values for the periodicities are predicted and they agree closely with observations. Some elements are suggested that might explain the well-documented 155-day periodicity. Conventional filtering methods to suppress effects of the 11-year cycle are criticized as inadequate.

Description: Data on solar wind electrons at the ISEE 3 spacecraft located 0.01 AU upstream from the earth (McComas et al., 1989) showed periods of time when the flux of antisunward suprathermal electrons would decrease suddenly, leading to heat flux dropouts (HFDs). This paper examines data from ISEE 1 at the 1.5 x 10 exp 6 km downstream location to determine whether HFDs identified at ISEE 3 by McComas et al. can be detected at this location and whether the ISEE 1 observations can provide information to one or the other possible interpretations of HFDs: that HFDs are due to enhanced Coulomb scattering, or to disconnection from the sun of the magnetic flux tube. The results of the examination identified the presence of HFD events in the ISEE 1 data, and the findings indicate that Coulomb scattering plays a substantial role in at least some HFD events.

Description: This paper discusses the biophysical stratification of the FIFE site, implementation of the stratification utilizing geographic information system methods, and validation of the stratification with respect to field measurements of biomass, Bowen ratio, soil moisture, and the greenness vegetation index (GVI) derived from TM satellite data. Maps of burning and topographic position were significantly associated with variation in GVI, biomass, and Bowen ratio. The stratified design did not significantly alter the estimated site-wide means for surface climate parameters but accounted for between 25 and 45 percent of the sample variance depending on the variable.

Description: During the intensive field campaigns of the First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment (FIFE) in May-October of 1987, several nearly simultaneous measurements were made with low-altitude flights of the L-band radiometer and C- and X-band scatterometers over two transects in the Konza Prairie Natural Research Area, some 8 km south of Manhattan, Kansas. These measurements showed that although the scatterometers were sensitive to soil moisture variations in most regions under the flight path, the L-band radiometer lost most of its sensitivity in regions unburned for many years. The correlation coefficient derived from the regression between the radar backscattering coefficient and the soil moisture was found to improve with the increase in antenna incidence angle. This is attributed to a steeper falloff of the backscattering coefficient as a function of local incidence at angles near nadir than at angles greater than 30 deg.

Description: The phased array has become a key component in the design of advanced antenna systems. This computer program was developed to examine the radiation characteristics of a generalized phased array antenna. Using a very efficient numerical technique, this program calculates the array's radiated power and its directivity. The results can be used to determine the radiation pattern of a generalized phased array at near- or far-field observation points. This program is a key research tool at the NASA Lewis Research Center for analyzing advanced space communication antenna systems. Results from this program compare favorably with experimental Lewis results for arrays of 2x2 and 3x3 elements. Given the array geometry and element characteristics, generalized phased array attributes can be broken into two areas: 1) the power radiated and its directivity at any given point, and 2) the co- and cross-polarization field components. This program allows arbitrarily located source elements with an analytically described cosine pattern. The formulation is based on a Romberg integration scheme and takes into account arbitrary element polarization, E and H plane element patterns, and mutual coupling. The input consists of the array geometry; phase, amplitude, linear and circular polarization of each source element; and the cosine functions of the element pattern. The output is a series of observation angles with their associated field magnitude and polarizations. Total radiated power and peak directivity are also calculated. This program is written in FORTRAN IV for batch execution and has been implemented on an IBM 370 computer operating under TSS with a central memory requirement of approximately 22K of 8 bit bytes. The IBM Scientific Subroutine Package (SSP) is required to run the program. This program was developed in 1986.

Description: The EM-ANIMATE program is a specialized visualization program that displays and animates the near-field and surface-current solutions obtained from an electromagnetics program, in particular, that from MOM3D (LAR-15074). The EM-ANIMATE program is windows based and contains a user-friendly, graphical interface for setting viewing options, case selection, file manipulation, etc. EM-ANIMATE displays the field and surface-current magnitude as smooth shaded color fields (color contours) ranging from a minimum contour value to a maximum contour value for the fields and surface currents. The program can display either the total electric field or the scattered electric field in either time-harmonic animation mode or in the root mean square (RMS) average mode. The default setting is initially set to the minimum and maximum values within the field and surface current data and can be optionally set by the user. The field and surface-current value are animated by calculating and viewing the solution at user selectable radian time increments between 0 and 2pi. The surface currents can also be displayed in either time-harmonic animation mode or in RMS average mode. In RMS mode, the color contours do not vary with time, but show the constant time averaged field and surface-current magnitude solution. The electric field and surface-current directions can be displayed as scaled vector arrows which have a length proportional to the magnitude at each field grid point or surface node point. These vector properties can be viewed separately or concurrently with the field or surface-current magnitudes. Animation speed is improved by turning off the display of the vector arrows. In RMS modes, the direction vectors are still displayed as varying with time since the time averaged direction vectors would be zero length vectors. Other surface properties can optionally be viewed. These include the surface grid, the resistance value assigned to each element of the grid, and the power dissipation of each element which has an assigned resistance value. The EM-ANIMATE program will accept up to 10 different surface current cases each consisting of up to 20,000 node points and 10,000 triangle definitions and will animate one of these cases. The capability is used to compare surface-current distribution due to various initial excitation directions or electric field orientations. The program can accept up to 50 planes of field data consisting of a grid of 100 by 100 field points. These planes of data are user selectable and can be viewed individually or concurrently. With these preset limits, the program requires 55 megabytes of core memory to run. These limits can be changed in the header files to accommodate the available core memory of an individual workstation. An estimate of memory required can be made as follows: approximate memory in bytes equals (number of nodes times number of surfaces times 14 variables times bytes per word, typically 4 bytes per floating point) plus (number of field planes times number of nodes per plane times 21 variables times bytes per word). This gives the approximate memory size required to store the field and surface-current data. The total memory size is approximately 400,000 bytes plus the data memory size. The animation calculations are performed in real time at any user set time step. For Silicon Graphics Workstations that have multiple processors, this program has been optimized to perform these calculations on multiple processors to increase animation rates. The optimized program uses the SGI PFA (Power FORTRAN Accelerator) library. On single processor machines, the parallelization directives are seen as comments to the program and will have no effect on compilation or execution. EM-ANIMATE is written in FORTRAN 77 for implementation on SGI IRIS workstations running IRIX 3.0 or later. A minimum of 55Mb of RAM is required for execution of this program; however, the code may be modified to accommodate the available memory of an individual workstation. For program execution, twenty-four bit, double-buffered color capability is suggested, but not required. Sample input and output files and a sample executable are provided on the distribution medium. Electronic documentation is provided in PostScript format and in the form of IRIX man pages. The standard distribution medium for EM-ANIMATE is a .25 inch streaming magnetic IRIX tape cartridge in UNIX tar format. EM-ANIMATE is also available as part of a bundled package, COS-10048 that includes MOM3D, an IRIS program that produces electromagnetic near field and surface current solutions. This program was developed in 1993.

Description: Four radiometric correction methods for the reduction of slope-aspect effects in a Landsat TM data set are tested in a mountainous test site with regard to their physical soundness and their influence on forest classification, as well as on the visual appearance of the scene. Excellent ground reference information and a fine-resolution DEM allowed precise assessment of the applicability of the methods under investigation. The results of the study presented here demonstrate the weakness of the classical cosine correction method for radiometric correction in rugged terrain. The statistical, Minnaert and C-correction approaches, however, yielded an improvement of the forest classification and an impressive reduction of the visual topography effect.

Description: Development of a first-order radiative transfer model for predicting backscatter from tree canopies has been underway at the University of Michigan Radiation Laboratory for some time. This model is known as the Michigan Microwave Canopy Scattering (MIMICS) model. This article presents the second-generation MIMICS model (MIMICS II) which accounts for canopies with discontinuous (open) crown layer geometries. MIMICS II models open crown layers by treating the location, size, and shape of the individual tree crowns as random variables. The backscattering coefficients for the canopy are then determined by introducing statistics derived from these parameters into the radiative transfer solution. Application of the radiative transfer equations to the discontinuous canopy geometry is presented. The resulting model is a robust fully polarimetric solution that is applicable over a wide variety of canopy architectures. Model simulations are compared to results generated with the continuous canopy model. The effect of the open crown geometry is found to be most significant at shallow incidence angles and at high frequencies for trees with well-developed crowns. Under these conditions, the gaps in the crown layer give rise to a notable increase in crown layer transmissivity which allows the radar to see through to the lower layers of the canopy more easily, thereby directly affecting the backscatter contribution of the trunks and ground.

Description: Experimental studies are performed on some coniferous trees (Austrian pine, Nordmann spruce, and Norway spruce) to investigate the relation between the tree architecture and radar signal at X-band. For a single tree, the RCS is measured as a function of the scatterer location at 90 deg incidence. It is found that the main scatterers are the leafy branches and the difference between sigma(vv) and sigma(hh) is significant at the upper portion of the tree. At the lower portion of the tree, sigma(vv) and sigma(hh) have almost the same level. For a group of trees the angular trends of sigma(vv) and sigma(hh) are measured. It is found that the levels of sigma(vv) and sigma(hh) are of the same order, but their angular trends vary from one tree species to the other depending on the tree species structure. The interpretation of these experimental results is carried out with the help of a theoretical model which accounts for the structure of the tree. According to this theoretical study, the major scattering trend is due to the leaves, while the perturbation to the angular trend and the level difference between sigma(vv) and sigma(hh) are due to the branch orientation distributions (i.e., the tree architecture).

Description: The dynamics of radially expanding magnetic clouds is rigorously analyzed within the framework of ideal MHD. The cloud is modelled as a cylindrically symmetric magnetic flux rope. In the force balance we include the gas pressure gradient and the Lorentz force. Interaction with the ambient solar wind due to expansion of the magnetic cloud is represented by a drag force proportional to the bulk velocity. We consider the self-similar expansion of a polytrope, and reduce the problem to an ordinary nonlinear differential equation for the evolution function. Analyzing the asymptotic behavior of the evolution function, we formulate theoretical expectations for the long-term behavior of cloud parameters. We focus on the temporal evolution of (1) the magnetic field strength; (2) the twist of the field lines; (3) the asymmetry of the total field profile; and (4) the bulk flow speed. We present data from two magnetic clouds observed at 1 AU and 2 AU, respectively, and find good agreement with theoretical expectations. For a peak magnetic field strength at 1 AU of 25 nT and a polytropic index of 0.5, we find that a magnetic cloud can be distinguished from the background interplanetary field up to a distance of about 5 AU. Taking larger magnetic fields and bigger polytropic indices this distance can double.

Description: Using unique 3-d velocity space measurements by the Ulysses solar wind plasma experiment from 1.15 to 5.34 AU, we assess the radial gradient in thermal electron temperature. Until 3.8 AU, the gradient was steeper than previously reported but flatter than adiabatic; after 3.8 AU the gradient flattened. Trends in the observed electron distribution shapes qualitatively support predictions for regulation by Coulomb collisions and by expansion in a spiral IMF.

Description: A radar system based on a network analyzer has been developed to study the backscatter from vegetation. The radar is operated at L-band. Radar measurements of a grass field were made in 1991. The radar returns from the grass were measured at three incidence angles. Ground truth and canopy parameters such as blade and stem dimensions, moisture content of the grass and the soil, and blade and stem density, were measured. These parameters are used in a distorted Born approximation model to compute the backscatter coefficients from the grass layer. The model results are compared with the radar data.

Description: The direction angle sensitivity of agricultural field backscatter is studied. The direction angle is defined as the angle between the incident plane and the perpendicular to the row direction. Maximum backscatter power from an angricultural field is expected to occur when the furrow induced slopes are oriented towards the radar, i.e., for a 0 deg direction angle. This effect is known as the cardinal effect. Because of the way the looks are formed in the NASA/JPL airborne synthetic aperture radar (AIRSAR) processor, each look corresponds to a slightly different squint angle. This provides a unique data set to analyze the cardinal effect, as it allows simultaneous observations of the backscatter of a field for sixteen different direction angles. The backscatter variations of the agricultural fields with direction angle at P-, L-, and C-bands is described. The observed variations in backscatter are compared with model predictions. The model predicts that the maximum backscatter occurs for a 0 deg direction angle, but underestimates the backscatter variations with direction angle by more than 10 dB.

Description: C-, L-, and P-band polarimetric signatures of wet snow surfaces have been analyzed, based on airborne synthetic aperture radar (AIRSAR) surveys of an Alpine test site. The importance of surface roughness is evident in the C- and L-band signatures, whereas the diffuse scattering contribution by internal inhomogeneities in the snowpack increases from the C- to the P-band at incidence angles below 50 deg due to increasing penetration.

Description: The Jet Propulsion Laboratory (JPL) airborne synthetic aperture radar (AIRSAR) L-band along-track interferometer images currents and waves on the ocean surface. By modifying the operating procedure of this two antenna interferometer, a technique has been developed to enable interferometric measurements to be made simultaneously at two different baselines. The availability of such data allows measurement of the decorrelation process of the ocean in greater detail. The coherence time of the ocean surface can be measured at high resolution over large areas. In addition to the L-band interferometer, a C-band along-track interferometer has been developed. It allows C-band dual-baseline measurements to be made simultaneous with the L-band measurements. The dual-baseline technique and AIRSAR implementation are described, and some example data are presented.

Description: The authors discuss the feasibility of determining the surface flux of sensible heat from forests with surface temperatures measured by satellites together with temperature soundings in the unstable atmospheric boundary layer (ABL). The latent heat flux is derived from the sensible heat flux by means of the energy budget. The study makes use of data collected during HAPEX-MOBILHY (Hydrologic Atmospheric Pilot Experiment-Modelization du Bilan Hydrique). The methodology is based on turbulence similarity for the unstable ABL. The surface temperature data were derived from measurements by the advanced very high resolution radiometer (AVHRR) aboard the NOAA-9 satellite; the atmospheric profiles were obtained by radiosondes.

Description: The authors examine the hypothesis that some part of the ecosystem-dependent variability of the vegetation indices is attributable to the effects of specular reflection of sunlight by leaves. A new class of vegetation indices, or 'minus specular' vegetation indices, is defined to account for the effects of specularly reflected light. Results show that the 'minus specular' indices, when compared to the traditional vegetation indices, potentially provide better estimates of the photosynthetic activity of a canopy than the traditional vegetation indices, particularly as a function of sun and view angles.

Description: The field experiment described took place in the Sayani Mountains of Siberia. The purpose of the joint field campaign was to observe and exchange methodologies with Russian scientists with regard to the development of remote sensing techniques for the early detection and assessment of forest decline damage believed to be associated with atmospheric deposition and/or insect and disease infestations. Several types of passive and active remote sensing measurements were made in conjunction with biophysical measurements on vegetative samples collected from four study sites representing a strong elevational gradient. Relatively cloud-free SPOT data were also acquired over the study area. Moderate canopy damage was recorded at the mid-elevation site (3400 ft/1037 m). The lowest levels of damage were recorded at the lowest elevation site (2300 ft/701 m.) At all sites, east versus west flagging of the canopy was noted (i.e., full canopy on the west-facing side of the canopy, significantly less foliage on the east-facing side).

Description: Earth photography from the Space Shuttle is used to examine the ice cover on Lake Baikal and correlate the patterns of weakened and melting ice with known hydrothermal areas in the Siberian lake. Particular zones of melted and broken ice may be surface expressions of elevated heat flow in Lake Baikal. The possibility is explored that hydrothermal vents can introduce local convective upwelling and disrupt a stable water column to the extent that the melt zones which are observed in the lake's ice cover are produced. A heat flow map and photographs of the lake are overlaid to compare specific areas of thinned or broken ice with the hot spots. The regions of known hydrothermal activity and high heat flow correlate extremely well with circular regions of thinned ice, and zones of broken and recrystallized ice. Local and regional climate data and other sources of warm water, such as river inlets, are considered.