ALBERT

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

Description: Landsat Thematic Mapper (TM) data were used to evaluate young conifer stands in the western Cascade Mountains of Oregon. Regression and correlation analyses were used to describe the relationships between TM band values and age of young Douglas-fir stands (2 to 35 years old). Spectral data from well regenerated Douglas-fir stands were compared to those of poorly regenerated conifer stands. TM bands 1, 2, 3, 5, 6, and 7 were inversely correlated with the age (r greater than or equal to -0.80) of well regenerated Douglas-fir stands. Overall, the 'structural index' (TM 4/5 ratio) had the highest correlation to age of Douglas-fir stands (r = 0.96). Poorly regenerated stands were spectrally distinct from well regenerated Douglas-fir stands after the stands reached an age of approximately 15 years.

Description: NOAA-9 satellite data from the Advanced Very High Resolution Radiometer (AVHRR) were used in conjunction with Landsat Multispectral Scanner (MSS) data to determine the proportion of closed canopy conifer forest cover in the Cascade Range of Oregon. A closed canopy conifer map, as determined from the MSS, was registered with AVHRR pixels. Regression was used to relate closed canopy conifer forest cover to AVHRR spectral data. A two-variable (band) regression model accounted for more variance in conifer cover than the Normalized Difference Vegetation Index (NDVI). The spectral signatures of various conifer successional stages were also examined. A map of Oregon was produced showing the proportion of closed canopy conifer cover for each AVHRR pixel. The AVHRR was responsive to both the percentage of closed canopy conifer cover and the successional stage in these temperate coniferous forests in this experiment.

Description: Thematic Mapper (TM) digital imagery was used to map forest successional stages and to evaluate spectral differences between old-growth and mature forests in the central Cascade Range of Oregon. Relative sun incidence values were incorporated into the successional stage classification to compensate for topographic induced variation. Relative sun incidence improved the classification accuracy of young successional stages, but did not improve the classification accuracy of older, closed canopy forest classes or overall accuracy. TM bands 1, 2, and 4; the normalized difference vegetation index (NDVI); and TM 4/3, 4/5, and 4/7 band ratio values for old-growth forests were found to be significantly lower than the values of mature forests (P less than or equal to 0.010). Wetness and the TM 4/5 and 4/7 band ratios all had low correlations to relative sun incidence (r(exp 2) less than or equal to 0.16). The TM 4/5 band ratio was named the 'structural index' (SI) because of its ability to distinguish between mature and old-growth forests and its simplicity.

Description: Digital Landsat Thematic Mapper (TM) and Satellite Probatoire d'Observation de la Terre (SPOT) High Resolution Visible (HRV) images of coniferous forest canopies were compared in their relationship to forest wood volume using correlation and regression analyses. Significant inverse relationships were found between softwood volume and the spectral bands from both sensors (P less than 0.01). The highest correlations were between the log of softwood volume and the near-infrared bands (HRV band 3, r = -0.89; TM band 4, r = -0.83).

Description: Forest and wildlife habitat analyses were conducted at the H.J. Andrews Experimental Forest in the Central Cascade Mountains of Oregon using remotely sensed data and a geographic information system (GIS). Landsat Thematic Mapper (TM) data were used to determine forest successional stages, and to analyze the structure of both old and young conifer forests. Two successional stage maps were developed. One was developed from six TM spectral bands alone, and the second was developed from six TM spectral bands and a relative sun incidence band. Including the sun incidence band in the classification improved the mapping accuracy in the two youngest successional stages, but did not improve overall accuracy or accuracy of the two oldest successional stages. Mean spectral values for old-growth and mature stands were compared in seven TM bands and seven band transformations. Differences between mature and old-growth successional stages were greatest for the band ratio of TM 4/5 (P = 0.00005) and the multiband transformation of wetness (P = 0.00003). The age of young conifer stands had the highest correlation to TM 4/5 values (r = 0.9559) of any of the TM band or band transformations used. TM 4/5 ratio values of poorly regenerated conifer stands were significantly different from well regenerated conifer stands after age 15 (P = 0.0000). TM 4/5 was named a 'Successional Stage Index' (SSI) because of its ability to distinguish forest successional stages. The forest successional stage map was used as input into a vertebrate richness model using GIS. The three variables of (1) successional stage, (2) elevation, and (3) site moisture were used in the GIS to predict the spatial occurrence of small mammal, amphibian, and reptile species based on primary and secondary habitat requirements. These occurrence or habitat maps were overlayed to tally the predicted number of vertebrate at any given point in the study area. Overall, sixty-three and sixty-seven percent of the model predictions for vertebrate occurrence matched the vertebrates that were trapped in the field in eight forested stands. Of the three model variables, site moisture appeared to have the greatest influence on the pattern of high vertebrate richness in all vertebrate classes.

Description: We compared two digital elevation models (DEM's) for the Echo Mountain SE quadrangle in the Cascade Mountains of Oregon. Comparisons were made between 7.5-minute (1:24,000-scale) and 1-degree (1:250,000-scale) images using the variables of elevation, slope aspect, and slope gradient. Both visual and statistical differences are presented.

Description: This paper describes the use of a set of spatial statistics to quantify the landscape pattern caused by the patchwork of clearcuts made over a 15-year period in the western Cascades of Oregon. Fifteen areas were selected at random to represent a diversity of landscape fragmentation patterns. Managed forest stands (patches) were digitized and analyzed to produce both tabular and mapped information describing patch size, shape, abundance and spacing, and matrix characteristics of a given area. In addition, a GIS fragmentation index was developed which was found to be sensitive to patch abundance and to the spatial distribution of patches. Use of the GIS-derived index provides an automated method of determining the level of forest fragmentation and can be used to facilitate spatial analysis of the landscape for later coordination with field and remotely sensed data. A comparison of the spatial statistics calculated for the two years indicates an increase in forest fragmentation as characterized by an increase in mean patch abundance and a decrease in interpatch distance, amount of interior natural forest habitat, and the GIS fragmentation index. Such statistics capable of quantifying patch shape and spatial distribution may prove important in the evaluation of the changing character of interior and edge habitats for wildlife.

Description: The flight of the Charleston County School District Can Do Project GeoCam payload on STS-57 was the climax of a decade long endeavor to bring the promise and excitement of the space program directly into the classroom. The payload carried four cameras designed to take high resolution photographs of the Earth under the direction of children operating the first ever student control room. During the course of the flight, the students followed the Shuttle's orbital tract, satellite weather images and selected a target list that was sent up to the crew each night as part of the execute package. Targets from this list, as well as ones chosen by the crew visually, resulted in the successful collection of photographic runs at many interesting sites on three on three continents.

Description: We examined the process of fragmentation in a managed forest landscape by comparing rates and patterns of disturbance (primarily clear-cutting) and regrowth between 1972 and 1988 using Landsat imagery. A 2589-km(exp 2) managed forest landscape in western Oregon was classified into two forest types, closed-canopy conifer forest (CF) (typically, greater than 60% conifer cover) and other forest and nonforest types (OT) (typically, less than 40 yr old or deciduous forest). The percentage of CF declined from 71 to 58% between 1972 and 1988. Declines were greatest on private land, least in wilderness, and intermediate in public nonwilderness. High elevations (greater than 914 m) maintained a greater percentage of CF than lower elevations (less than 914 m). The percentage of the area at the edge of the two cover types increased on all ownerships and in both elevational zones, whereas the amount of interior habitat (defined as CF at least 100 m from OT) decreased on all ownerships and elevational zones. By 1988 public lands contained approximately 45% interior habitat while private lands had 12% interior habitat. Mean interior patch area declined from 160 to 62 ha. The annual rate of disturbance (primarily clear-cutting) for the entire area including the wilderness was 1.19%, which corresponds to a cutting rotation of 84 yr. The forest landscape was not in a steady state or regulated condition which is not projected to occur for at least 40 yr under current forest plans. Variability in cutting rates within ownerships was higher on private land than on nonreserve public land. However, despite the use of dispersed cutting patterns on public land, spatial patterns of cutting and remnant forest patches were nonuniform across the entire public ownership. Large remaining patches (less than 5000 ha) of contiguous interior forest were restricted to public lands designated for uses other than timber production such as wilderness areas and research natural areas.

Description: We investigated how the amount of old-growth and mature forest influences the selection of nest sites by northern spotted owls (Strix occidentalis caurina) in the Central Cascade Mountains of Oregon. We used 7 different plot sizes to compare the proportion of mature and old-growth forest between 30 nest sites and 30 random sites. The proportion of old-growth and mature forest was significantly greater at nests sites than at random sites for all plot sizes (P less than or equal to 0.01). Thus, management of the spotted owl might require setting the percentage of old-growth and mature forest retained from harvesting at least 1 standard deviation above the mean for the 30 nest sites we examined.

Description: There has been an increased interest in the quantification of pattern in ecological systems over the past years. This interest is motivated by the desire to construct valid models which extend across many scales. Spatial methods must quantify pattern, discriminate types of pattern, and relate hierarchical phenomena across scales. Wavelet analysis is introduced as a method to identify spatial structure in ecological transect data. The main advantage of the wavelet transform over other methods is its ability to preserve and display hierarchical information while allowing for pattern decomposition. Two applications of wavelet analysis are illustrated, as a means to: (1) quantify known spatial patterns in Douglas-fir forests at several scales, and (2) construct spatially-explicit hypotheses regarding pattern generating mechanisms. Application of the wavelet variance, derived from the wavelet transform, is developed for forest ecosystem analysis to obtain additional insight into spatially-explicit data. Specifically, the resolution capabilities of the wavelet variance are compared to the semi-variogram and Fourier power spectra for the description of spatial data using a set of one-dimensional stationary and non-stationary processes. The wavelet cross-covariance function is derived from the wavelet transform and introduced as a alternative method for the analysis of multivariate spatial data of understory vegetation and canopy in Douglas-fir forests of the western Cascades of Oregon.

Description: In contrast to traditional Synthetic Aperture Radar (SAR), an Interferometric SAR (INSAR) yields two additional measurements: the phase difference and the correlation between the two interferometric channels. The phase difference has been used to estimate topographic height. For homogeneous surfaces, the correlation depends on the system signal-to-noise (SNR) ratio, the interferometer parameters, and the local slope. In the presence of volume scattering, such as that encountered in vegetation canopies, the correlation between the two channels is also dependent on the degree of penetration of the radiation into the scattering medium. In this paper, we propose a method for removing system and slope effects in order to obtain the decorrelation due to penetration alone. The sensitivities and accuracy of the proposed method are determined by Monte Carlo experiments, and we show that the proposed technique has sufficient sensitivity to provide penetration measurements for airborne SAR systems. Next, we provide a theoretical model to estimate the degree of penetration in a way which is independent of the details of the scattering medium. We also present a model for the correlation from non-homogeneous layers. We assess the sensitivity of the proposed inversion technique to these inhomogeneous situations. Finally, we present a comparison of the interferometric results against in situ data obtained by an airborne laser profilometer which provides a direct measurement of tree height and an estimate of the vegetation density profile in the forested areas around Mt. Adams, WA.

Description: In recent years, monitoring vegetation biomass over various climate zones has become the primary focus of several studies interested in assessing the role of the ecosystem responses to climate change and human activities. Airborne and spaceborne Synthetic Aperture Radar (SAR) systems provide a useful tool to directly estimate biomass due to its sensitivity to structural and moisture characteristics of vegetation canopies. Even though the sensitivity of SAR data to total aboveground biomass has been successfully demonstrated in many controlled experiments over boreal forests and forest plantations, so far, no biomass estimation algorithm has been developed. This is mainly due to the fact that the SAR data, even at lowest frequency (P-band) saturates at biomass levels of about 200 tons/ha, and the structure and moisture information in the SAR signal forces the estimation algorithm to be forest type dependent. In this paper, we discuss the development of a hybrid forest biomass algorithm which uses a SAR derived land cover map in conjunction with a forest backscatter model and an inversion algorithm to estimate forest canopy water content. It is shown that unlike the direct biomass estimation from SAR data, the estimation of water content does not depend on the seasonal and/or environmental conditions. The total aboveground biomass can then be derived from canopy water content for each type of forest by incorporating other ecological information. Preliminary results from this technique over several boreal forest stands indicate that (1) the forest biomass can be estimated with reasonable accuracy, and (2) the saturation level of the SAR signal can be enhanced by separating the crown and trunk biomass in the inversion algorithm. We have used the JPL AIRSAR data over BOREAS southern study area to test the algorithm and to generate regional scale water content and biomass maps. The results are compared with ground data and the sources of errors are discussed. Several SAR images in synoptic modes are used to generate the parameter maps. The maps are then combined to generate mosaic maps over the BOREAS modeling grid.

Description: We report the Bayes classification of terrain types at different sites using airborne interferometric synthetic aperture radar (INSAR) data. A Gaussian maximum likelihood classifier was applied on multidimensional observations derived from the SAR intensity, the terrain elevation model, and the magnitude of the interferometric correlation. Training sets for forested, urban, agricultural, or bare areas were obtained either by selecting samples with known ground truth, or by k-means clustering of random sets of samples uniformly distributed across all sites, and subsequent assignments of these clusters using ground truth. The accuracy of the classifier was used to optimize the discriminating efficiency of the set of features that was chosen. The most important features include the SAR intensity, a canopy penetration depth model, and the terrain slope. We demonstrate the classifier's performance across sites using a unique set of training classes for the four main terrain categories. The scenes examined include San Francisco (CA) (predominantly urban and water), Mount Adams (WA) (forested with clear cuts), Pasadena (CA) (urban with mountains), and Antioch Hills (CA) (water, swamps, fields). Issues related to the effects of image calibration and the robustness of the classification to calibration errors are explored. The relative performance of single polarization Interferometric data classification is contrasted against classification schemes based on polarimetric SAR data.

Description: Many algorithms for spectral analysis of imaging spectroscopy data of the Earth's surface require that the data be calibrated to surface reflectance. Calibration requires removing instrumental response, solar irradiance, atmospheric transmittance, and atmospheric scattering from the radiance detected at the sensor. Depending on the amount of support data, this can be a formidable task. This paper examines four methods of calibration: (1) a radiative transfer model from the University of Colorado (ATREM: Gao and Goetz, 1990; Gao et al., 1992), (2) a MODTRAN-based method developed at the Jet Propulsion Lab by Green et al., (1191), (3) a ground calibration using known sites as standards, and (4) a combined approach using radiative transfer methods and ground calibration. Data from the Airborne Visual and Infra-Red Imaging Spectrometer (AVIRIS) instrument were evaluated from data sets obtained over multiple years and multiple sites.

Description: This report documents the final reliability prediction performed on the Earth Observing System/Advanced Microwave Sounding Unit-A (EOS/AMSU-A). The A1 Module contains Channels 3 through 15, and is referred to herein as 'EOS/AMSU-A1'. The A2 Module contains Channels 1 and 2, and is referred herein as 'EOS/AMSU-A2'. The 'specified' figures were obtained from Aerojet Reports 8897-1 and 9116-1. The predicted reliability figure for the EOS/AMSU-A1 meets the specified value and provides a Mean Time Between Failures (MTBF) of 74,390 hours. The predicted reliability figure for the EOS/AMSU-A2 meets the specified value and provides a MTBF of 193,110 hours.

Description: A suborbital experiment was designed to study the photochemistry of the mesosphere by observing simultaneously the airglow emissions with in-situ minor species number density profiles. The experiment was very successful and some preliminary results have already been reported in various scientific meetings. Two scientific papers are currently in the process of final preparation for submission for publication. In this final project report, we will first give a background description of the experiment and follow by the summaries of the scientific papers currently being prepared.

Description: This document constitutes the final technical report for the National Aeronautics and Space Administration (NASA) Grant NAGW-3172. This grant was instituted to provide for the conduct of research under the Universities Space Research Association's (USRA's) Universities Earth System Scientist Program (UESSP) for the Office of Mission to Planet Earth (OMTPE) at NASA Headquarters. USRA was tasked with the following requirements in support of the Universities Earth System Scientists Programs: (1) Bring to OMTPE fundamental scientific and technical expertise not currently resident at NASA Headquarters covering the broad spectrum of Earth science disciplines; (2) Conduct basic research in order to help establish the state of the science and technological readiness, related to NASA issues and requirements, for the following, near-term, scientific uncertainties, and data/information needs in the areas of global climate change, clouds and radiative balance, sources and sinks of greenhouse gases and the processes that control them, solid earth, oceans, polar ice sheets, land-surface hydrology, ecological dynamics, biological diversity, and sustainable development; (3) Evaluate the scientific state-of-the-field in key selected areas and to assist in the definition of new research thrusts for missions, including those that would incorporate the long-term strategy of the U.S. Global Change Research Program (USGCRP). This will, in part, be accomplished by study and evaluation of the basic science needs of the community as they are used to drive the development and maintenance of a global-scale observing system, the focused research studies, and the implementation of an integrated program of modeling, prediction, and assessment; and (4) Produce specific recommendations and alternative strategies for OMTPE that can serve as a basis for interagency and national and international policy on issues related to Earth sciences.

Description: A key drawback to estimating geodetic and geodynamic parameters over time based on satellite laser ranging (SLR) observations is the inability to accurately model all the forces acting on the satellite. Errors associated with the observations and the measurement model can detract from the estimates as well. These 'model errors' corrupt the solutions obtained from the satellite orbit determination process. Dynamical models for satellite motion utilize known geophysical parameters to mathematically detail the forces acting on the satellite. However, these parameters, while estimated as constants, vary over time. These temporal variations must be accounted for in some fashion to maintain meaningful solutions. The primary goal of this study is to analyze the feasibility of using a sequential process noise filter for estimating geodynamic parameters over time from the Laser Geodynamics Satellite (LAGEOS) SLR data. This evaluation is achieved by first simulating a sequence of realistic LAGEOS laser ranging observations. These observations are generated using models with known temporal variations in several geodynamic parameters (along track drag and the J(sub 2), J(sub 3), J(sub 4), and J(sub 5) geopotential coefficients). A standard (non-stochastic) filter and a stochastic process noise filter are then utilized to estimate the model parameters from the simulated observations. The standard non-stochastic filter estimates these parameters as constants over consecutive fixed time intervals. Thus, the resulting solutions contain constant estimates of parameters that vary in time which limits the temporal resolution and accuracy of the solution. The stochastic process noise filter estimates these parameters as correlated process noise variables. As a result, the stochastic process noise filter has the potential to estimate the temporal variations more accurately since the constraint of estimating the parameters as constants is eliminated. A comparison of the temporal resolution of solutions obtained from standard sequential filtering methods and process noise sequential filtering methods shows that the accuracy is significantly improved using process noise. The results show that the positional accuracy of the orbit is improved as well. The temporal resolution of the resulting solutions are detailed, and conclusions drawn about the results. Benefits and drawbacks of using process noise filtering in this type of scenario are also identified.

Description: Given the possibility that the passive Microwave Humidity Sounder (MHS) for NASA's Earth Observing System (EOS) may have to be replaced, two alternate configurations are analyzed in this document. One option mirrors AMSU-B, the Advanced Microwave Sounding Unit (AMSU) moisture sounder; which consists of five channels at 89, 150, 183 +/- 1, 183 +/- 3, and 183 +/- 7 GHz. The second option contains an additional channel at 183 +/- 14 GHz and replaces the channels at 89 and 150 GHz with three channels at 118 +/- 0.5, 118 +/- 1.4, and 118 +/- 3 GHz. The latter configuration is considered to be superior due to its greater scientific benefits and reduced cost and complexity -- it would require only two local oscillators and a smaller antenna size for a given resolution on the ground. As shown by means of simulations, humidity profile retrieval accuracy for the second option is superior in most cases and only slightly degraded relative to the first option in the worst-case scenario. The results are summarized. Including the 118-GHz channels also offers the possibility of cell-top altitude retrievals and improved temperature profile retrievals when used in conjunction with a temperature sounder such as the 15-channel AMSU-A.

Description: In this report the Committee on Earth Studies (CES), a standing committee of the Space Studies Board (SSB) within the National Research Council (NRC), reviews the recent history (nominally from 1981 to 1995) of the U.S. earth observations programs that serve civilian needs. The principal observations programs examined are those of NASA and the National Oceanic and Atmospheric Administration (NOAA). The Air Force' s Defense Meteorological Satellite Program (DMSP) is discussed, but only from the perspective of its relationship to civil needs and the planned merger with the NOAA polar-orbiting system. The report also reviews the interfaces between the earth observations satellite programs and the major national and international environmental monitoring and research programs. The monitoring and research programs discussed are the U.S. Global Change Research Program (USGCRP), the International Geosphere-Biosphere Program (IGBP), the World Climate Research Program (WCRP), related international scientific campaigns, and operational programs for the sharing and application of environmental data. The purpose of this report is to provide a broad historical review and commentary based on the views of the CES members, with particular emphasis on tracing the lengthy record of advisory committee recommendations. Any individual topic could be the subject of an extended report in its own right. Indeed, extensive further reviews are already under way to that end. If the CES has succeeded in the task it has undertaken. This report will serve as a useful starting point for any such more intensive study. The report is divided into eight chapters: ( I ) an introduction, (2) the evolution of the MTPE, (3) its relationship to the USGCRP, (4) applications of earth observations data, (5) the role that smaller satellites can play in research and operational remote sensing, (6) earth system modeling and information systems, (7) a number of associated activities that contribute to the MTPE and the USGCRP, and (8) organizational issues in the conduct of civil earth observations programs. Following the body of the report is a series of appendixes: after a list of acronyms and abbreviations and collected short biographies of CES members, six brief tutorials discuss several scientific topics important to the science and applications of earth observations. Highlights from the eight chapters follow.

Description: The continuation of the EOSDIS testbed ('Testbed') has materialized from a multi-task system to a fully functional stand-alone data archive distribution center that once was only X-Windows driven to a system that is accessible by all types of users and computers via the World Wide Web. Throughout the past months, the Testbed has evolved into a completely new system. The current system is now accessible through Netscape, Mosaic, and all other servers that can contact the World Wide Web. On October 1, 1995 we will open to the public and we expect that the statistics of the type of user, where they are located, and what they are looking for will drastically change. What is the most important change in the Testbed has been the Web interface. This interface will allow more users access to the system and walk them through the data types with more ease than before. All of the callbacks are written in such a way that icons can be used to easily move around in the programs interface. The homepage offers the user the opportunity to go and get more information about each satellite data type and also information on free programs. These programs are grouped into categories for types of computers that the programs are compiled for, along with information on how to FTP the programs back to the end users computer. The heart of the Testbed is still the acquisition of satellite data. From the Testbed homepage, the user selects the 'access to data system' icon, which will take them to the world map and allow them to select an area that they would like coverage on by simply clicking that area of the map. This creates a new map where other similar choices can be made to get the latitude and longitude of the region the satellite data will cover. Once a selection has been made the search parameters page will appear to be filled out. Afterwards, the browse image will be called for once the search is completed and the images for viewing can be selected. There are several other option pages, but once an order has been selected the Testbed will bring up the order list page and the user will then be able to place their order. After the order has been completed, the Testbed will mail the user to notify them of the completed order and how the images can be picked up.

Description: Combining quantitative geophysical information extracted from the optical and microwave wavelengths provides complementary information about both the surface mineralogy and morphology. This study combines inversion results from two remote sensing instruments, a polarimetric synthetic aperture radar, AIRSAR, and an imaging spectrometer, AVIRIS, for Trail Canyon alluvial fan in Death Valley, California. The NASA/JPL Airborne Synthetic Aperture Radar (AIRSAR) is a quad-polarization, three frequency instrument. AIRSAR collects data at C-band = 5.66 cm, L-band = 23.98 cm, and P-band = 68.13 cm. The data are processed to four-looks and have a spatial resolution of 10 m and a swath width of 12 km. The AIRSAR data used in this study were collected as part of the Geologic Remote Sensing Field Experiment (GRSFE) over Death Valley on 9/14/89. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is a NASA/JPL instrument that flies in an ER-2 aircraft at an altitude of 20 km. AVIRIS uses four spectrometers to collect data in 224 spectral channels from 0.4 micrometer to 2.45 micrometer. The width of each spectral band is approximately 10 nm. AVIRIS collects data with a swath width of 11 km and a pixel size of 20 m. The AVIRIS data used in this study were collected over Death Valley on 5/31/92.

Description: Accurate wavelength calibration of imaging spectrometer data is essential if proper atmospheric transmission corrections are to be made to obtain apparent surface reflectance. Accuracies of 0.1 nm are necessary for a 10 nm-sampling instrument in order to match the slopes of the deep atmospheric water vapor features that dominate the 0.7-2.3 micrometer regions. The Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) is calibrated in the laboratory to determine the wavelength position and full-width-half-maximum (FWHM) response for each of the 224 channels. The accuracies are limited by the quality of the monochromator used as a source. The accuracies vary from plus or minus to plus or minus 1.5 nm depending on the wavelength region, in general decreasing with increasing wavelength. Green et al. make corrections to the wavelength calibrations by using the known positions of 14 atmospheric absorption features throughout the 0.4-2.5 micrometer wavelength region. These features, having varying width and intensity, were matched to the MODTRAN model with a non-linear least squares fitting algorithm. A complete calibration was developed for all 224 channels by interpolation. Instrument calibration cannot be assumed to be stable to 0.1 nm over a flight season given the rigors of thermal cycling and launch and landing loads. The upcoming sensor HYDICE will require a means for in-flight spectral calibration of each scene because the calibration is both temperature and pressure sensitive. In addition, any sensor using a two-dimensional array has the potential for systematic wavelength shifts as a function of cross-track position, commonly called 'smile'. Therefore, a rapid means for calibrating complete images will be required. The following describes a method for determining instrument wavelength calibration using atmospheric absorption features that is efficient enough to be used for entire images on workstations. This study shows the effect of the surface reflectance on the calibration accuracy and the calibration history for the AVIRIS B spectrometer over the 1992 flight season.

Description: This study examines various aspects of the Microwave Vegetation Index (MVI). MVI is a derived signal created by differencing the spectral response of the 37 GHz horizontally and vertically polarized passive microwave signals. The microwave signal employed to derive this index is thought to be primarily influenced by vegetation structure, vegetation growth, standing water, and precipitation. The state of California is the study site for this research. Imagery from the Special Sensor Microwave/Imager (SSM/I) is used for the creation of MVI datasets analyzed in this research. The object of this research is to determine whether MVI corresponds with some quantifiable vegetation parameter (such as vegetation density) or whether the index is more affected by known biogeophysical parameters such antecedent precipitation. A secondary question associated with the above is whether the vegetation attributes that MVI is employed to determine can be more easily and accurately evaluated by other remote sensing means. An important associated question to be addressed in the study is the effect of different multi-temporal composting techniques on the derived MVI dataset. This work advances our understanding of the fundamental nature of MVI by studying vegetation as a mixture of structural types, such as forest and grassland. The study further advances our understanding by creating multitemporal precipitation datasets to compare the affects of precipitation upon MVI. This work will help to lay the groundwork for the use of passive microwave spectral information either as an adjunct to visible and near infrared imagery in areas where that is feasible or for the use of passive microwave alone in areas of moderate cloud coverage. In this research, an MVI dataset, spanning the period February 15, 1989 through April 25, 1990, has been created using National Aeronautic and Space Administration (NASA) supplied brightness temperature data. Information from the DMSP satellite 37 GHz wavelength SSM/I sensor in both horizontal and vertical polarization has been processed using the MVI algorithm. In conjunction with the MVI algorithm a multitemporal compositing technique was used to create datasets that correspond to 14 day periods. In this technical report, Section Two contains background information on the State of California and the three MVI study sites. Section Three describes the methods used to create the MVI and independent variables datasets. Section Four presents the results of the experiment. Section Five summarizes and concludes the work.

Description: The object of this research is to use indirect measurements, notably thermal infrared, to describe urbanization and deforestation with parameters that can be used to assess, as well as predict, the effects of land use changes on local microclimate. More specifically, we use a new approach for the treatment of remotely sensed data; this is referred to as the 'triangle' method. The name triangle is given because the envelope of data points, when plotted as a function of surface radiant temperature versus vegetation index or fractional vegetation cover, exhibits the shape of a triangle. From the information contained on these 'scatter plots', land use changes can be related to two intrinsic surface variables, the surface moisture availability (M(sub 0))(sup 1) and fractional vegetation cover. Recent work by Carlson et al. indicate that the triangle shape on the scatter plots may be scale similar, suggesting that these two parameters are subject to the same interpretation on differing scales. A second objective in this research is to determine if historical data for Advanced Very High Resolution Radiometer (AVHRR) (NOAA satellite; 1.1 km resolution at nadir) can be used to assess changes in regional land use over time. To this end, two target areas were chosen for the investigation of urbanization and two for deforestation. The former comprise tow areas in Pennsylvania, one a small but rapidly growing population center (State College) and the other a medium-sized urban area which continues to undergo development (Chester County). The two deforestation sites consist of rain forest areas in western and central Costa Rica and a region in the Brazilian Amazon.

Description: The goal of this project is to develop satellite-based observational techniques for measuring both oceanic and atmospheric variables using passive polarimetric radiometry. Polarimetric radiometry offers a potential alternative to radar scatterometry in observing global ocean surface wind direction from satellites. Polarimetric radiometry might also provide a means of detecting cell-top ice in convective storms by virtue of the polarizing properties of oriented ice particles, and thus facilitate estimation of the phase of the storm. The project focuses on the development of polarimetric microwave radiometers using digital cross-correlators for obtaining precise measurements of all four Stokes' parameters. As part of the project a unique four-band polarimetric imaging radiometer, the Polar Scanning Radiometer (PSR), is being designed for use on the NASA DC-8 aircraft. In addition to providing an aircraft-based demonstration of digital correlation technology the PSR will significantly enhance the microwave imaging capability of the existing suite of DC-8 instruments. During the first grant year excellent progress has been made in the following areas: (1) demonstrating digital correlation radiometry, (2) fabricating aircraft-qualified correlators for use in the PSR, and (3) modeling observed SSM/I brightness signatures of ocean wind direction.

Description: The research on model functions for ADEOS and ERS-1 are summarized and an analysis of the differences between the three kinds of models is provided in this final report. The success of the AMI on ERS-1 obtained at GSFC and NMC is highlighted. The problem of wind stress description is reviewed within and the scatterometer model being developed for high winds monitoring for the AMI on ERS-1 and ERS-2 is described.

Description: The data processing analysis and conversion of polarization measurements to polarization parameters from the Sunflower multiband radiometer is presented in this final report. Included is: (1) the actual data analysis; (2) the comparison of the averaging techniques and the percent polarization derived from the original and averaged I, Q, U parameters; (3) the polarizer angles used in conversion; (4) the Matlab files; (5) the relative ground size, field of view location, and view zenith angles, and (6) the summary of all the sky data for all dates.

Description: During the second phase project year we have made progress in the development and refinement of surface temperature retrieval algorithms and in product generation. More specifically, we have accomplished the following: (1) acquired a new advanced very high resolution radiometer (AVHRR) data set for the Beaufort Sea area spanning an entire year; (2) acquired additional along-track scanning radiometer(ATSR) data for the Arctic and Antarctic now totalling over eight months; (3) refined our AVHRR Arctic and Antarctic ice surface temperature (IST) retrieval algorithm, including work specific to Greenland; (4) developed ATSR retrieval algorithms for the Arctic and Antarctic, including work specific to Greenland; (5) developed cloud masking procedures for both AVHRR and ATSR; (6) generated a two-week bi-polar global area coverage (GAC) set of composite images from which IST is being estimated; (7) investigated the effects of clouds and the atmosphere on passive microwave 'surface' temperature retrieval algorithms; and (8) generated surface temperatures for the Beaufort Sea data set, both from AVHRR and special sensor microwave imager (SSM/I).

Description: Synoptic observations of ice motion in the Arctic Basin are currently limited to those acquired by drifting buoys and, more recently, radar data from ERS-1. Buoys are not uniformly distributed throughout the Arctic, and SAR coverage is currently limited regionally and temporally due to the data volume, swath width, processing requirements, and power needs of the SAR. Additional ice-motion observations that can map ice responses simultaneously over large portions of the Arctic on daily to weekly time intervals are thus needed to augment the SAR and buoys data and to provide an intermediate-scale measure of ice drift suitable for climatological analyses and ice modeling. Principal objectives of this project were to: (1) demonstrate whether sufficient ice features and ice motion existed within the consolidated ice pack to permit motion tracking using AVHRR imagery; (2) determine the limits imposed on AVHRR mapping by cloud cover; and (3) test the applicability of AVHRR-derived motions in studies of ice-atmosphere interactions. Each of these main objectives was addressed. We conclude that AVHRR data, particularly when blended with other available observations, provide a valuable data set for studying sea ice processes. In a follow-on project, we are now extending this work to cover larger areas and to address science questions in more detail.

Description: An improved method for detecting, measuring, and distinguishing crop residue, live vegetation, and mineral soil is presented. By measuring fluorescence in multiple bands, live and dead vegetation are distinguished. The surface of the ground is illuminated with ultraviolet radiation, inducing fluorescence in certain molecules. The emitted fluorescent emission induced by the ultraviolet radiation is measured by means of a fluorescence detector, consisting of a photodetector or video camera and filters. The spectral content of the emitted fluorescent emission is characterized at each point sampled, and the proportion of the sampled area covered by residue or vegetation is calculated.

Description: During Phase 2 airborne Gamma Radiation measurements were conducted over 48 BOREAS flight lines. Ground measurements of the soil moisture (SM) of the mineral soil and of the depth and water content (WC) of the moss/humus layer were collected for calibration of the flight lines. Special attention has been given to the flight lines over and near the primary tower sites in the southern study area (SSA): Old Black Spruce (OBS), Old Jack Pine, Old Aspen, and Young Aspen. Multiple ground surveys (September 1993; July, August, and September 1994) show the variation of the water content of the moss/humus layers and changes in the amount of standing water near the SSA, OBS tower during the period September 1993 to September 1994. All ground data, airborne estimates and locations of flight lines and ground sampling points have been submitted to BORIS (Boreas Information System). On 8-10 September 1994 exceptionally high values of cosmic radiation were observed by the airborne gamma radiation system over the SSA. Follow-up investigation has not determined the source of, or what caused, the high cosmic count rates.

Description: This progress report covers the following efforts initiated for the year: year-round monthly soil CO2 flux measurements were started in both primary and secondary forests and in managed and degraded pastures; root sorting and weighing has begun and all four ecosystems at Paragominas have been analyzed through samples; regional modeling of soil water dynamics and minimum rooting depth has been done and the RADAMBRASIL soils database has been digitized and a 20 year record of the precipitation for the region has been produced, along with a hydrological ('bucket-tipping') model that will run within a GIS framework; prototype tension lysimeters have been designed and installed in soil pits to begin assessing the importance of DOC as a source of organic matter in deep soils; and many publications, listed in this document, have resulted from this year's research. Two of the papers published are included with this annual report document.

Description: A significant progress has been made in TIR instrumentation which is required to establish the spectral BRDF/emissivity knowledge base of land-surface materials and to validate the land-surface temperature (LST) algorithms. The SIBRE (spectral Infrared Bidirectional Reflectance and Emissivity) system and a TIR system for measuring spectral directional-hemispherical emissivity have been completed and tested successfully. Optical properties and performance features of key components (including spectrometer, and TIR source) of these systems have been characterized by integrated use of local standards (blackbody and reference plates). The stabilization of the spectrometer performance was improved by a custom designed and built liquid cooling system. Methods and procedures for measuring spectral TIR BRDF and directional-hemispheric emissivity with these two systems have been verified in sample measurements. These TIR instruments have been used in the laboratory and the field, giving very promising results. The measured spectral emissivities of water surface are very close to the calculated values based on well established water refractive index values in published papers. Preliminary results show that the TIR instruments can be used for validation of the MODIS LST algorithm in homogeneous test sites. The beta-3 version of the MODIS LST software is being prepared for its delivery scheduled in the early second half of this year.

Description: An investigation of the influence of stratospheric aerosol on the performance of the atmospheric correction algorithm was carried out. The results indicate how the performance of the algorithm is degraded if the stratospheric aerosol is ignored. Use of the MODIS 1380 nm band to effect a correction for stratospheric aerosols was also studied. The development of a multi-layer Monte Carlo radiative transfer code that includes polarization by molecular and aerosol scattering and wind-induced sea surface roughness has been completed. Comparison tests with an existing two-layer successive order of scattering code suggests that both codes are capable of producing top-of-atmosphere radiances with errors usually less than 0.1 percent. An initial set of simulations to study the effects of ignoring the polarization of the the ocean-atmosphere light field, in both the development of the atmospheric correction algorithm and the generation of the lookup tables used for operation of the algorithm, have been completed. An algorithm was developed that can be used to invert the radiance exiting the top and bottom of the atmosphere to yield the columnar optical properties of the atmospheric aerosol under clear sky conditions over the ocean, for aerosol optical thicknesses as large as 2. The algorithm is capable of retrievals with such large optical thicknesses because all significant orders of multiple scattering are included.

Description: During Phase 2 airborne Gamma Radiation measurements were conducted over 48 BOREAS (boreal ecosystem atmosphere study) flight lines. Ground measurements of the soil moisture (SM) of the mineral soil and the depth and water content (WC) of the moss/humus layer were collected for calibration of the flight lines. Special attention has been given to the flight lines over, and near, the primary tower sites in the southern study area (SSA); Old Black Spruce (OBS), Old Jack Pine, Old Aspen, and Young Aspen. Multiple ground surveys (September 1993, July, August, and September 1994) show the variation of the water content of the moss/humus layers and changes in the amount of standing water near the SSA OBS tower during the period September 1993 to September 1994. All ground data, airborne estimates and locations of flight lines and ground sampling points have been submitted to BORIS (BOREAS Information System). On 8-10 September 1994 exceptionally high values of cosmic radiation were observed by the airborne gamma radiation system over the SSA. Follow up investigation has not determined the source of, or what caused, the high cosmic count rates.

Description: This document is a summary of the project funded by NAGw-1460 as part of the Earth Observation Commericalization/Applications Program (EOCAP) directed by NASA's Earth Science and Applications Division. The goal was to work with several agencies to focus on forest structure and landscape characterizations for wildlife habitat applications. New analysis techniques were used in remote sensing and landscape ecology with geographic information systems (GIS). The development of GIS and the emergence of the discipline of landscape ecology provided us with an opportunity to study forest and wildlife habitat resources from a new perspective. New techniques were developed to measure forest structure across scales from the canopy to the regional level. This paper describes the project team, technical advances, and technology adoption process that was used. Reprints of related refereed journal articles are in the Appendix.

Description: Several significant accomplishments were made during the present reporting period. (1) Initial simulations to understand the applicability of the MODerate Resolution Imaging Spectrometer (MODIS) 1380 nm band for removing the effects of stratospheric aerosols and thin cirrus clouds were completed using a model for an aged volcanic aerosol. The results suggest that very simple procedures requiring no a priori knowledge of the optical properties of the stratospheric aerosol may be as effective as complex procedures requiring full knowledge of the aerosol properties, except the concentration which is estimated from the reflectance at 1380 nm. The limitations of this conclusion will be examined in the next reporting period; (2) The lookup tables employed in the implementation of the atmospheric correction algorithm have been modified in several ways intended to improve the accuracy and/or speed of processing. These have been delivered to R. Evans for implementation into the MODIS prototype processing algorithm for testing; (3) A method was developed for removal of the effects of the O2 'A' absorption band from SeaWiFS band 7 (745-785 nm). This is important in that SeaWiFS imagery will be used as a test data set for the MODIS atmospheric correction algorithm over the oceans; and (4) Construction of a radiometer, and associated deployment boom, for studying the spectral reflectance of oceanic whitecaps at sea was completed. The system was successfully tested on a cruise off Hawaii on which whitecaps were plentiful during October-November. This data set is now under analysis.

Description: There were two major developments in 1994 in our collaborative GPS experiment in the Tien Shan of the Former Soviet Union (FSU). Both were motivated by our expectation that we will ultimately obtain better science at lower cost if we involve our colleagues in the FSU more deeply in (1) the collection and (2) the analysis of data. As an experimental test of the concept of having our local collaborators carry out the field work semi-autonomously, we sent 6 MIT receivers to the Tien Shan for a period of 3 months. To enable our collaborators to have the capability for data analysis, we provided computers for two data analysis centers and organized a two-week training session. This report emphasizes the rationale for deeper involvement of FSU scientists, describes the training sessions, discusses the data collection, and presents the results. We also discuss future plans. More detailed discussion of background, general scientific objectives, discussions with collaborators, and results for the campaigns in 1992 and 1993 have been given in previous reports.

Description: Efforts continue under this contract to develop algorithms for the computation of sea surface temperature (SST) from MODIS infrared retrievals. This effort includes radiative transfer modeling, comparison of in situ and satellite observations, development and evaluation of processing and networking methodologies for algorithm computation and data accession, evaluation of surface validation approaches for IR radiances, and participation in MODIS (project) related activities. Efforts in this contract period have focused on radiative transfer modeling, evaluation of atmospheric correction methodologies, involvement in field studies, production and evaluation of new computer networking strategies, and objective analysis approaches.

Description: During this quadrennium, very long baseline interferometry (VLBI) data acquisition and system development has focused on improving the accuracy of the system sufficiently to allow the determination of reliable estimates of height variations. Associated with this aim has been improvements in the determinations of horizontal velocity fields, monitoring water vapor delay using interferometric methods, and improvements to Earth rotation measurements. The primary aims of the improvements to height measurement accuracy have been to directly measure the contemporary magnitudes of post glacial rebound, and to determine a height reference system for measuring global sea level rise. High frequency Earth rotation studies have been carried out to better define the transformation parameters from an inertial coordinate system to an Earth fixed one, and to better understand the coupling between the components of the atmosphere-ocean-solid Earth system. Two major VLBI campaigns were carried out in support of these studies: (1) Epoch-92 in July 1992 and (2) Cont-94 in January 1994. Each of these campaigns lasted approximately two weeks and involved multiple VLBI networks operating simultaneously in addition to other space geodetic systems operating during these periods. Two major compilations of the VLBI results (and results from other space geodetic systems) have been published during this quadrennium.

Description: We have classified four different images, under various levels of JPEG compression, using the following classification algorithms: minimum-distance, maximum-likelihood, and neural network. The training site accuracy and percent difference from the original classification were tabulated for each image compression level, with maximum-likelihood showing the poorest results. In general, as compression ratio increased, the classification retained its overall appearance, but much of the pixel-to-pixel detail was eliminated. We also examined the effect of compression on spatial pattern detection using a neural network.

Description: This document describes the progress of the task of the Millimeter-wave Imaging Radiometer (MIR) data processing and the development of water vapor retrieval algorithms, for the second six-month performing period. Aircraft MIR data from two 1995 field experiments were collected and processed with a revised data processing software. Two revised versions of water vapor retrieval algorithm were developed, one for the execution of retrieval on a supercomputer platform, and one for using pressure as the vertical coordinate. Two implementations of incorporating products from other sensors into the water vapor retrieval system, one from the Special Sensor Microwave Imager (SSM/I), the other from the High-resolution Interferometer Sounder (HIS). Water vapor retrievals were performed for both airborne MIR data and spaceborne SSM/T-2 data, during field experiments of TOGA/COARE, CAMEX-1, and CAMEX-2. The climatology of water vapor during TOGA/COARE was examined by SSM/T-2 soundings and conventional rawinsonde.

Description: It has become increasingly important to conduct interdisciplinary environmental research assessments, both nationally and internationally. For this reason, the Committee for a Pilot Study on Database Interfaces was charged to review and advise on data interfacing activities. The committee used six case studies (1) to identify and understand the most important problems associated with collecting, integrating, and analyzing environmental data from local to global spatial scales and over a very wide range of temporal scales; and (2) to elaborate the common barriers to interfacing data of disparate sources and types. Consistent with the committee's charge, the primary focus was the interfacing of geophysical and ecological data. The case studies used by the committee were: The Impact Assessment Project for Drought Early Warning in the Sahel, The National Acid Precipitation Assessment Program, The H.J. Andrew Experimental Forest Long-Term Ecological Research Site, The Carbon Dioxide Information Analysis Center, The First International Satellite Land Surface Climatology Project (ISLSCP) Field Experiment, and The California Cooperative Oceanic Fisheries Investigation. The committee derived a number of lessons from the case studies, and these lessons are summarized at the end of each case study and analyzed in the last chapter. Some are generic in nature; others are more specific to a discipline or project. The conclusions and recommendations are based on the committee's analysis of the case studies and additional research. They are organized according to four major areas of barriers or challenges to the effective interfacing of diverse environmental data. In the final section the committee offers a set of broadly applicable principles (Ten Keys to Success) that can be used by scientists and data managers in planning and conducting interfacing activities.

Description: Through the CaPE Hydrometeorology Project, we have developed an understanding of some of the unique data quality issues involved in assimilating data of disparate types for regional-scale hydrologic modeling within a GIS framework. Among others, the issues addressed here include the development of adequate validation of the surface water budget, implementation of the STATSGO soil data set, and implementation of a remote sensing-derived landcover data set to account for surface heterogeneity. A model of land surface processes has been developed and used in studies of the sensitivity of surface fluxes and runoff to soil and landcover characterization. Results of these experiments have raised many questions about how to treat the scale-dependence of land surface-atmosphere interactions on spatial and temporal variability. In light of these questions, additional modifications are being considered for the Marshall Land Surface Processes Model. It is anticipated that these techniques can be tested and applied in conjunction with GCIP activities over regional scales.

Description: A few of the major accomplishments during the second phase include: (1) all three Special Sensor Microwave Imagers (SSM/I's: F08, F10, and F11) have been cross-calibrated; (2) a very large, quality-controlled, collocated SSM/I and radiosonde data set has been produced; (3) the SSM/I-radiosonde and SSM/I-buoy data sets have been used to calibrate the SSM/I ocean retrieval algorithm; (4) ocean products have been produced for both F10 and F11 SSM/I for 1991-1993; (5) the SSM/I-buoy data set was used to better determine the variation of the ocean T(sub B) with wind direction; and (6) it was demonstrated that under high wind conditions, wind direction information can be obtained from individual SSM/I observations.

Description: A method has been presented for the representation of sub-grid scale variability of surface properties within a land surface processes model. The method uses remotely-sensed data to directly or indirectly estimate probability density functions (PDF's) or key surface variables. Application of this technique in a coupled land surface-atmosphere model requires only grid-scale values of the variables of interest, obtained from low-resolution satellite imagery or surface/remote sensing data assimilation. The PDF's of each controlling surface property are superimposed on the respective grid-scale values to simulate sub-grid scale heterogeneity. Sensitivity studies will be carried out to ascertain the relative importance of the heterogeneity of several variables, and the degree to which non-linear property-process interactions impact large-scale fluxes.

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.