ALBERT

Description: This report describes the transitional activities of the JPL Analysis Center.

Description: In August and September of 1995 the Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) was deployed to Brazil as part of the NASA Smoke Cloud Aerosol and Radiation experiment in Brazil (SCAR-B). AVIRIS measures spectra from 400 to 2500 nm at 10-nm intervals. These spectra are acquired as images with dimensions of 11 by up to 800 km with 20-m spatial resolution. Spectral images measured by AVIRIS are spectrally, radiometrically, and spatially calibrated. During the SCAR-B deployment, AVIRIS measured more than 300 million spectra of regions of Brazil. A portion of these spectra were acquired over areas of actively burning fires. Actively burning fires emit radiance in the AVIRIS spectral range as a function of temperature. This emitted radiance is expressed from the 2500-nm end of the AVIRIS spectrum to shorter wavelengths as a function of intensity and modeled by the Planck function.. The objective of this research and analysis was to use spectroscopic methods to determine the minimum high temperature of the most intense fires measured in the SCAR-B AVIRIS data set. Spectra measured by AVIRIS with hot sources have been previously examined for volcanic lava and fires in Brazil.

Description: Various locations in the southwestern U.S. are used to calibrate remote sensing instruments. This study shows how some of these targets compare in terms of albedo and homogeneity, and records the variation of these factors for a single location (Ivanpah Playa) over a period of one year. Results indicate that there is a great deal of variation among these targets in albedo, spectral flatness, and surface uniformity, and that these factors can change throughout the year.

Description: The purpose of the chemistry component of the model comparison is to assess to what extent differences in the formulation of chemical processes explain the variance between model results. Observed concentrations of chemical compounds are used to estimate to what degree the various models represent realistic situations. For readability, the materials for the chemistry experiment are reported in three separate sections. This section discussed the data used to evaluate the models in their simulation of the source gases and the Nitrogen compounds (NO(y)) and Chlorine compounds (Cl(y)) species.

Description: Space-based and airborne coherent Doppler lidars designed for measuring global tropospheric wind profiles in cloud-free air rely on backscatter, beta from aerosols acting as passive wind tracers. Aerosol beta distribution in the vertical can vary over as much as 5-6 orders of magnitude. Thus, the design of a wave length-specific, space-borne or airborne lidar must account for the magnitude of 8 in the region or features of interest. The SPAce Readiness Coherent Lidar Experiment under development by the National Aeronautics and Space Administration (NASA) and scheduled for launch on the Space Shuttle in 2001, will demonstrate wind measurements from space using a solid-state 2 micrometer coherent Doppler lidar. Consequently, there is a critical need to understand variability of aerosol beta at 2.1 micrometers, to evaluate signal detection under varying aerosol loading conditions. Although few direct measurements of beta at 2.1 micrometers exist, extensive datasets, including climatologies in widely-separated locations, do exist for other wavelengths based on CO2 and Nd:YAG lidars. Datasets also exist for the associated microphysical and chemical properties. An example of a multi-parametric dataset is that of the NASA GLObal Backscatter Experiment (GLOBE) in 1990 in which aerosol chemistry and size distributions were measured concurrently with multi-wavelength lidar backscatter observations. More recently, continuous-wave (CW) lidar backscatter measurements at mid-infrared wavelengths have been made during the Multicenter Airborne Coherent Atmospheric Wind Sensor (MACAWS) experiment in 1995. Using Lorenz-Mie theory, these datasets have been used to develop a method to convert lidar backscatter to the 2.1 micrometer wavelength. This paper presents comparison of modeled backscatter at wavelengths for which backscatter measurements exist including converted beta (sub 2.1).

Description: The typical fair weather electric field at the ground is between -100 and -300 V/m. At the NASA Kennedy Space Center and US Air Force Cape Canaveral Air Station (KSC) the electric field at the ground sometimes reaches -400 to -1200 V/m within an hour or two after sunrise on days that otherwise seem to be fair weather. We refer to the enhanced negative electric fields as the "sunrise enhancement." To investigate the sunrise enhancement at KSC we measured the electric field (E) in the first few hundred meters above the ground before and during several sunrise enhancements. From these E soundings we can infer the presence of charge layers and determine their thickness and charge density.

Description: Early observations with ERS-1 SAR image data revealed a large ice stream in North East Greenland (Fahnestock 1993). The ice stream has a number of the characteristics of the more closely studied ice streams in Antarctica, including its large size and gross geometry. The onset of rapid flow close to the ice divide and the evolution of its flow pattern, however, make this ice stream unique. These features can be seen in the balance velocities for the ice stream (Joughin 1997) and its outlets. The ice stream is identifiable for more than 700 km, making it much longer than any other flow feature in Greenland. Our research goals are to gain a greater understanding of the ice flow in the northeast Greenland ice stream and its outlet glaciers in order to assess their impact on the past, present, and future mass balance of the ice sheet. We will accomplish these goals using a combination of remotely sensed data and ice sheet models. We are using satellite radar interferometry data to produce a complete maps of velocity and topography over the entire ice stream. We are in the process of developing methods to use these data in conjunction with existing ice sheet models similar to those that have been used to improve understanding of the mechanics of flow in Antarctic ice streams.

Description: This research is focusing on two related areas that are fundamental to the NASA PARCA (Program for Arctic Regional Climate Assessment) program. The primary research area is the determination of the amount, rate, and timing of accumulation at distributed sites in the dry snow zone of Greenland and evaluation of these results in light of accumulation modeling results. The secondary research area is the calibration of the isotope "thermometer" at these ice sheet sites as well as the determination of long-term temperature trends in Greenland.

Description: The TARFOX (Tropospheric Aerosol Radiative Forcing Observational Experiment) intensive field campaign was designed to reduce uncertainties in estimates of the effects of anthropogenic aerosols on climate by measuring direct radiative effects and the optical, physical, and chemical properties of aerosols [1]. TARFOX was conducted off the East Coast of the United States between July 10-31, 1996. Ground, aircraft, and satellite-based sensors measured the sensitivity of radiative fields at various atmospheric levels to aerosol optical properties (i.e., optical thickness, phase function, single-scattering albedo) and to the vertical profile of aerosols. The LASE (Lidar Atmospheric Sensing Experiment) instrument, which was flown on the NASA ER-2 aircraft, measured vertical profiles of total scattering ratio and water vapor during a series of 9 flights. These profiles were used in real-time to help direct the other aircraft to the appropriate altitudes for intensive sampling of aerosol layers. We have subsequently used the LASE aerosol data to derive aerosol backscattering and extinction profiles. Using these aerosol extinction profiles, we derived estimates of aerosol optical thickness (AOT) and compared these with measurements of AOT from both ground and airborne sun photometers and derived from the ATSR-2 (Along Track and Scanning Radiometer 2) sensor on ERS-2 (European Remote Sensing Satellite-2). We also used the water vapor mixing ratio profiles measured simultaneously by LASE to derive precipitable water vapor and compare these to ground based measurements.

Description: The Jovian magnetosphere with its strong magnetic field and the rapid rotation of the planet present new opportunities and challenges for the use of electrodynamic tethers. An overview of the basic plasma physics properties of an electrodynamic tether moving through the Jovian magnetosphere is examined. Tether use for both propulsion and power generation are considered. Close to the planet, tether propulsive forces are found to be as high as 50 Newtons and power levels as high as 1 million Watts.

Description: The Automated Fluid Interface System (AFIS) is an advanced development prototype satellite servicer. The device was designed to transfer consumables from one spacecraft to another. An engineering model was built and underwent development testing at Marshall Space Flight Center. While the current AFIS is not suitable for spaceflight, testing and evaluation of the AFIS provided significant experience which would be beneficial in building a flight unit.

Description: Among the most important short-term dynamic biological processes are diurnal changes in canopy water relations. Plant regulation of water transport through stomatal openings affects other gaseous transport processes, often dramatically decreasing photosynthetic fixation of carbon dioxide during periods of water stress. Water stress reduces stomatal conductance of water vapor through the leaf surface and alters the diurnal timing of stomatal opening. Under non-water stressed conditions, stomates typically open soon after dawn and transpire water vapor throughout the daylight period. During stress periods, stomates may close for part of the day, generally near mid-day. Under prolonged stress conditions, stomatal closure shifts to earlier times during the day; stomates may close by mid-morning and remain closed until the following morning - or remain closed entirely. Under these conditions the relationship between canopy greenness (e.g., measured with a vegetation index or by spectral mixture analysis) and photosynthetic fixation of carbon is lost and the remotely sensed vegetation metric is a poor predictor of gas exchange. Prediction of stomatal regulation and exchange of water and trace gases is critical for ecosystem and climate models to correctly estimate budgets of these gases and understand or predict other processes like gross and net ecosystem primary production. Plant gas exchange has been extensively studied by physiologists at the leaf and whole plant level and by biometeorologists at somewhat larger scales. While these energy driven processes follow a predictable if somewhat asymmetric diurnal cycle dependent on soil water availability and the constraints imposed by the solar energy budget, they are nonetheless difficult to measure at the tree and stand levels using conventional methods. Ecologists have long been interested in the potential of remote sensing for monitoring physiological changes using multi-temporal images. Much of this research has focused on day-to-day changes in water use, especially for agricultural applications. Ustin et al. showed seasonal changes in canopy water content in chaparral shrub could be estimated using optical methods. Vanderbilt et al. followed asymmetric diurnal changes in the reflectance of a walnut orchard, but could not attribute specific reflectance changes to specific changes in canopy architecture or physiology. Forests and shrub lands in California experience prolonged periods of drought, sometimes extending six months without precipitation. The conifer and evergreen chaparral communities common to the foothill region around the central valley of California retain their foliage throughout the summer and have low transpiration rates despite high net radiation and temperature conditions. In contrast, grasslands and drought resistant deciduous species in the same habitat are seasonally dormant in summer. Because of differences in the mechanisms of drought tolerance, rooting depth and physiology between different plant communities in the region, it is likely that they display differences in diurnal water relations. The presence of diverse plant communities provides an opportunity to investigate possible diurnal landscape patterns in water relations that could be observed by an airborne hyperspectral scanner. This investigation of AVIRIS data collected over forest and shrub land represents the continuation of a prior investigation involving spectral mixture analysis of diurnal effects in the same AVIRIS data set.

Description: Spectroscopy is used in the laboratory to measure the molecular components and concentrations of plant constituents to answer questions about the plant type, status, and health. Imaging spectrometers measure the upwelling spectral radiance above the Earth's surface as images. Ideally, imaging spectrometer data sets should be used to understand plant type, plant status, and health of plants in an agricultural setting. An Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data set was acquired over agricultural fields near Wallula, Washington on July 23rd, 1997. AVIRIS measures upwelling radiance spectra through 224 spectral channels with contiguous 10-nm sampling from 400 to 2500 nm in the solar-reflected spectrum. The spectra are measured as images of 11 by up to 800 km with 20-m spatial resolution. The spectral images measured by AVIRIS represent the integrated signal resulting from: the solar irradiance; two way transmittance and scattering of the atmosphere; the absorptions and scattering of surface materials; as well as the spectral, radiometric and spatial response functions of AVIRIS. This paper presents initial research to derive properties of the agricultural fields near Wallula from the calibrated spectral images measured by AVIRIS near the top of the atmosphere.

Description: This report summarizes the current and future plans of the Crustal Dynamics Data Information System (CDDIS) with respect to the International VLBI Service for Geodesy and Astrometry (IVS). Included are background information about the CDDIS, the computer architecture, staffing supporting the system, archive contents, and future plans for the CDDIS within the IVS.

Description: This document presents a system controlling the motion of a spherical air bearing used in the modeling of spacecraft dynamics and controls in a laboratory environment. The system is part of the Spinning Rocket Simulator (SRS), used to simulate the coning of spacecraft during a thrusting stage. The reaction force at the spherical air bearing supporting the spacecraft model must coincide with the thrust axis of the model for proper simulation. Therefore, the bearing is translated in a circular path to introduce a centrifugal force. This horizontal force along with the gravitational reaction force at the bearing combines to simulate the direction of the spacecraft's thrust force. The control system receives attitude information from the spacecraft model via a laser beam embedded in the model that impinges on a photosensitive array. The non-linear system is controlled using high-speed lookup tables and digital techniques. A vector-controlled motor and a stepper motor are given the necessary signals to accurately control the turntable and platform supporting the air bearing. Preliminary performance data is presented. Mechanical elements of the table and platform are described in detail. A wireless (RF) data path for all devices on the spacecraft model to an off-table command computer is also described.

Description: NASA's Cross-Cutting Technology Development Program identified formation flying as a key enabler for the next generation Earth and Sciences campaign. It is hoped that this technology will allow a distributed network of autonomous satellites to act collaboratively as a single collective unit paving the way for extensive co-observing campaigns, coordinated multi-point observing programs, improved space-based interferometry, and entirely new approaches to conducting science. APL as a team member with GSFC, funded by the Earth Sciences and Technology Organization (ESTO), investigated formation deployment and initialization concepts which is central to the formation flying concept. This paper presents the analytical approach and preliminary results of the study. The study investigated a simple mission involving the deployment of six micro-satellites, one at a time, from a bus. At the initialization state, the satellites fly in an along-track trajectory separated by nominal spacing. The study entailed the development of a two-body (bus and satellite) relative motion propagator based on Clohessy-Wiltshire (C-W) equations with drag from which the relative motion of the micro-satellites is deduced. This code was used to investigate cluster development characteristics subject to "tip-off' (ejection) conditions. Results indicate that cluster development is very sensitive to the ballistic coefficients of the bus and satellites, and to relative ejection velocity. This information can be used to identify optimum deployment parameters, along with accuracy bounds for a particular mission, and to develop a cluster control strategy minimizing global fuel and cost. A suitable control strategy concept has been identified, however, it needs to be developed further.

Description: Breakthrough technology development is critical to securing the future of our space industry. The National Aeronautics and Space Administration (NASA) Cross-Enterprise Technology Development Program (CETDP) is developing critical space technologies that enable innovative and less costly missions, and spawn new mission opportunities through revolutionary, long-term, high-risk, high-payoff technology advances. The CETDP is a NASA-wide activity managed by the Advanced Technology and Mission Studies Division (AT&MS) at Headquarters Office of Space Science. Program management for CETDP is distributed across the multiple NASA Centers and draws on expertise throughout the Agency. The technology research activities are organized along Project-level divisions called thrust areas that are directly linked to the Agency's goals and objectives of the Enterprises: Earth Science, Space Science, Human Exploration and Development of Space; and the Office of the Chief Technologist's (OCT) strategic technology areas. Cross-Enterprise technology is defined as long-range strategic technologies that have broad potential to span the needs of more than one Enterprise. Technology needs are identified and prioritized by each of the primary customers. The thrust area manager (TAM) for each division is responsible for the ultimate success of technologies within their area, and can draw from industry, academia, other government agencies, other CETDP thrust areas, and other NASA Centers to accomplish the goals of the thrust area. An overview of the CETDP and description of the future directions of the thrust area called Distributed Spacecraft are presented in this paper. Revolutionary technologies developed within this thrust area will enable the implementation of a spatially distributed network of individual vehicles, or assets, collaborating as a single collective unit, and exhibiting a common system-wide capability to accomplish a shared objective. With such a capability, new Earth and space science measurement concepts become a reality.

Description: The overview of the International Space Station (ISS) is comprised of the program vision and mission; Space Station uses; definition of program phases; as well as descriptions and status of several scheduled International Space Station Overview assembly flights.

Description: In recent years, the role of PSCs in the ozone depletion process has become better understood. PSCs provide the surfaces upon which heterogeneous reactions take place that affect the gas phase partitioning between active and reservoir chlorine and nitrogen species. Present methods of PSC detection include in situ measurements by lidar and various satellite-borne instruments such as the Stratospheric Aerosol Measurement II (SAM II) on the Nimbus 7 spacecraft, which produced PSC measurements from 1978 to 1994 and several instruments onboard the Upper Atmospheric Research Satellite (UARS) such as the Cryogenic Limb Array Etalon Spectrometer (CLAES) which provided measurements for 1991-1993. All of the PSC-detection methods devised so far have been hampered by incomplete sampling of the places and times in which PSCs are likely to form. There is a need to understand the climatology of PSCS, in particular the timing of their onset and duration, their vertical distribution, geographic extent, annual variability and responses to volcanic aerosol forcing. Poole and Pitts [1994] assembled a PSC climatology based on SAM II data, but this climatology is incomplete, as it is limited to the edge of the polar night due to the limitations of the solar occultation scan geometry. The Advanced Very High Resolution Radiometer (AVHRR) five- channel sensors onboard the NOAA polar-orbiting satellites have been collecting data over the polar regions continuously since 1979. These operational satellites provide unmatched coverage in space and time of both polar regions, but were not designed for the detection of optically-thin PSCS. However, the AVHRR data archive would be an invaluable source for the construction of a long-term climatology of PSCs if techniques can be developed and tested to detect PSCs in AVHRR data. In the last few years, the members of our group at San Francisco State University and NASA Ames Research Center have been engaged in the development of various PSC detection methods using AVHRR data. There is strong evidence that a subset of PSCS, those that are optically thick, can be readily identified in the AVHRR data set. Our group has also made significant progress in the identification of optically thinner PSCs using a variety of techniques.

Description: The Southern Great Plains 1997 (SGP97) field experiment was conducted in Oklahoma during June-July 1997 to validate the models used for computing remote soil moisture using measurements by microwave radiometers. One of the objectives of SGP97 was to examine the effect of soil moisture on the evolution of the Atmospheric Boundary Layer (ABL) and clouds over the Southern Great Plains (SGP) during the warm season. The LASE (Lidar Atmospheric Sensing Experiment) airborne DIAL (Differential Absorption Lidar) system, which was flown autonomously on the NASA ER-2 aircraft during previous missions, was reconfigured to fly on the NASA P3 research aircraft. During SGP97 LASE was used to study the morning evolution of the ABL, particularly as manifested in the development of the convective boundary layer, and to study the influence of soil moisture variations on the development of ABL. The ABL development is strongly influenced by the surface energy budget, which is in turn influenced by soil moisture, mesoscale meteorology, clouds, and solar insolation. LASE data acquired during this mission are being used to study the ABL water vapor budget, the development of the ABL, spatial and temporal variabilities in the ABL, and the meteorological factors that influence the ABL development. This field experiment also permitted comparisons of LASE water vapor measurements with water vapor profiles acquired by radiosondes launched at the DOE (Department of Energy) Atmospheric Radiation Measurement (ARM) Southern Great Plain (SGP) site and at NASA/Wallops Flight Facility, as well as with measurements from other SGP97 aircraft.

Description: Lidar backscattering profiles available from the LITE data set have been used to estimate the optical depths of the Saharan dust layer over West Africa and E. Atlantic regions, in the context of validating the 3-D conceptual model of the Saharan dust plume proposed by Karyampudi and Carlson. The aerosol extinction profiles and optical depths were retrieved from LITE using the Fernald et al. (1972) method. An extinction-to-backscattering ratio, S(sub a), of 25 was selected for optical depth calculations. The spatial analysis of total column and Saharan dust layer optical depths show higher optical depths over W. Africa that decrease westward over E. Atlantic. The higher optical depths over W. Africa, in general, are associated with heavy dust being raised from the surface in dust source regions. Rapid depletion of these heavy dust particles, perhaps due to sedimentation, appear to decrease the dust loading within the dust layer as the plume leaves the west African continent. Higher optical depths are generally confined to the southern edge of the dust layer, where the middle level jet appears to transport the heavy dust concentrations that tend to mix downward from vertical mixing associated with the strong vertical shears underneath the middle jet. Thus, LITE measurements although, in general, validate the Saharan dust plume conceptual model, show maximum values of optical depths near the southern edge of the dust plume over the E. Atlantic region instead of near the center of the dust plume as described in the conceptual model.

Description: Water vapor plays an important role in the energetics of the boundary layer processes which in turn play a key role in regulating regional and global climate. It plays a primary role in Earth's hydrological cycle, in radiation balance as a direct absorber of infrared radiation, and in atmospheric circulation as a latent heat energy source, as well as in determining cloud development and atmospheric stability. Water vapor concentration, expressed as a mass mixing ratio (g kg(exp -l)), is conserved in all meteorological processes except condensation and evaporation. This property makes it an ideal choice for studying many of the atmosphere's dynamic features. Raman scattering measurements from lidar also allow retrieval of water vapor mixing ratio profiles at high temporal and vertical resolution. Raman lidars sense water vapor to altitudes not achievable with towers and surface systems, sample the atmosphere at much higher temporal resolution than radiosondes or satellites, and do not require strong vertical gradients or turbulent fluctuations in temperature that is required by acoustic sounders and radars. Analysis of highly-resolved water vapor profiles are used here to characterize two important mesoscale flows: thunderstorm outflows and a cold front passage. The data were obtained at the Atmospheric Radiation Measurement Site (CART) by the groundbased Department of Energy/Sandia National Laboratories lidar (CART Raman lidar or CARL) and Goddard Space Flight Center Scanning Raman Lidar (SRL). A detailed discussion of the SRL and CARL performance during the IOPs is given by others in this meeting.

Description: The NASA Langley Research Center's airborne UV Differential Absorption Lidar (DIAL) system participated in the Subsonic Assessment, Ozone and Nitrogen Oxide Experiment (SONEX) mission from October 13 to November 12, 1997. The purpose of the mission was to study the upper troposphere/lower stratosphere in and near the North Atlantic flight corridor to better understand this region of the atmosphere and how civilian air travel in the corridor might be affecting the atmospheric chemistry. Bases of operations included NASA Ames, California (37.4 deg N, 122.1 deg W); Bangor, Maine (44.8 deg N, 68.8 deg W); Shannon, Ireland (52.7 deg N, 8.9 deg W); and Lajes, Terceira Island, Azores (38.8 deg N, 27.1 deg W). Since the UV DIAL system observes in the nadir as well as the zenith, aerosol and ozone data were obtained from near the Earth's surface to the lower stratosphere. A number of interesting features were noted relating to both chemistry and dynamics of the troposphere, which are reported here.

Description: Worldwide, about ten Differential Absorption Lidars are used for long-term monitoring of stratospheric ozone. These systems are an important component of the Network for the Detection of Stratospheric Change. Although DIALs are self-calibrating in principle, regular intercomparisons with other ozone-lidars, microwave radiometers or ozone-sondes are highly desirable to ensure high data quality at a well known level. The Network for the Detection of Stratospheric Change (NDSC) validation policy suggests that such intercomparisons be "blind", meaning all participants submit their data to an impartial referee, without seeing results from the other participants. Here we report on the "blind" intercomparison taking place from January 20th to February 10th 1998 at Ny-Alesund, Spitsbergen (78.92 deg N, 11.95 deg E). Participating groups were from the Alfred Wegener Institute, Potsdam, operating the NDSC DIAL system at Ny-Alesund, from the University of Bremen operating the NDSC microwave radiometer for ozone profiling at Ny-Alesund, and the NASA Goddard Space Flight Center group with the "NDSC travelling standard" STROZ-LITE. The first author acted as the impartial referee. Also used for the intercomparison were data from ECC-6A/Vaisala RS80 ozone sondes routinely launched at Ny-Alesund by the AWI group. A 1% KI solution (3 ml) and the 1986 ECC pump correction (1.092 at 5 hPa) are used. The ECC-data were available to all participants during the campaign and thus were not "blind". Table 1 summarizes the expected performance of the instruments participating in the ozone intercomparison reported in this paper.

Description: Orbiter towing provides a backup reboost capability for the International Space Station (ISS). Results from recent studies are presented, showing performance, system configuration, mission operations, and programmatics. A proposed flight demonstration to mitigate risks is also discussed.

Description: Imaging Spectroscopy enables the identification and mapping of surface mineralogy over large areas. This study focused on assessing the utility of Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) data for environmental impact analysis over the Environmental Protection Agency's (EPA) high priority Superfund site Ray Mine, AZ. Using the Spectral Angle Mapper (SAM) algorithm to analyze AVIRIS data makes it possible to map surface materials that are indicative of acid generating minerals. The improved performance of the AVIRIS sensor since 1996 provides data with sufficient signal to noise ratio to characterize up to 8 image endmembers. Specifically we employed SAM to map minerals associated with mine generated acid waste, namely jarositc, goethite, and hematite, in the presence of a complex mineralogical background.

Description: The primary objective of the Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) Program at the Southern Great Plains (SGP) Cloud and Radiation Testbed (CART) is to acquire in situ and remote sensing data to improve cloud and atmospheric radiative models and parameterizations. As a consequence of this program, a large number of atmosphere and surface measurements are being acquired at the ARM SGP CART central site. NASA's Airborne Visible/Infrared Imaging Spectrometer (AVIRIS) overflew this site on August 1, 1997. AVIRIS measures the upwelling spectral solar radiance from 400 to 2500 nm at 10-nm intervals. From 20 km altitude, these calibrated spectra are acquired as images of 11 by up to 800 km with 20-by-20 m spatial resolution. These data were acquired at the ARM SGP CART Central Site to first investigate derivation of atmospheric parameters from the measured spectra, second study the variation of these parameters, and third demonstrate the inversion of the calibrated radiance spectra to apparent surface reflectance. These objectives have been pursued with AVIRIS data at other sites for atmospheric water vapor and derivation of apparent surface reflectance.

Description: Mineral maps generated for the Ray Mine, Arizona were analyzed to determine if imaging spectroscopy can provide accurate information for environmental management of active and abandoned mine regions. The Ray Mine, owned by the ASARCO Corporation, covers an area of 5700 acres and is situated in Pinal County, Arizona about 70 miles north of Tucson near Hayden, Arizona. This open-pit mine has been a major source of copper since 1911, producing an estimated 4.5 million tons of copper since its inception. Until 1955 mining was accomplished by underground block caving and shrinkage stope methods. (excavation by working in stepped series usually employed in a vertical or steeply inclined orebody) In 1955, the mine was completely converted to open pit method mining with the bulk of the production from sulfide ore using recovery by concentrating and smelting. Beginning in 1969 a significant production contribution has been from the leaching and solvent extraction-electrowinnowing method of silicate and oxide ores. Published reserves in the deposit as of 1992 are 1.1 billion tons at 0.6 percent copper. The Environmental Protection Agency, in conjunction with ASARCO, and NASA/JPL obtained AVIRIS data over the mine in 1997 as part of the EPA Advanced Measurement Initiative (AMI) (Tom Mace, Principal Investigator). This AVIRIS data set is being used to compare and contrast the accuracy and environmental monitoring capabilities of remote sensing technologies: visible-near-IR imaging spectroscopy, multispectral visible and, near-IR sensors, thermal instruments, and radar platforms. The goal of this effort is to determine if these various technologies provide useful information for envirorunental management of active and abandoned mine sites in the arid western United States. This paper focuses on the analysis of AVIRIS data for assessing the impact of the Ray Mine on Mineral Creek. Mineral Creek flows to the Gila River. This paper discusses our preliminary AVIRIS mineral mapping and environmental findings.

Description: This section contains a number of special diagnostics that are designed to examine certain mechanisms. Section 1 reports on the method used to test the photochemical partitioning in the models. Sections 2 and 3 represent efforts to examine the model calculated production and removal rates for ozone and how the values are combined with transport rates in the models to produce the simulated ozone distributions. Sections 4 and 5 concentrate on polar processes including the dynamics aspect of vortex confinement and the chemical aspects of chlorine activation.

Description: MM II defined a series of experiments to better understand and characterize model transport and to assess the realism of this transport by comparison to observations. Measurements from aircraft, balloon, and satellite, not yet available at the time of MM I [Prather and Remsberg, 1993], provide new and stringent constraints on model transport, and address the limits of our transport modeling abilities. Simulations of the idealized tracers the age spectrum, and propagating boundary conditions, and conserved HSCT-like emissions probe the relative roles of different model transport mechanisms, while simulations of SF6 and C02 make the connection to observations. Some of the tracers are related, and transport diagnostics such as the mean age can be derived from more than one of the experiments for comparison to observations. The goals of the transport experiments are: (1) To isolate the effects of transport in models from other processes; (2) To assess model transport for realistic tracers (such as SF6 and C02) for comparison to observations; (3) To use certain idealized tracers to isolate model mechanisms and relationships to atmospheric chemical perturbations; (4) To identify strengths and weaknesses of the treatment of transport processes in the models; (5) To relate evaluated shortcomings to aspects of model formulation. The following section are included:Executive Summary, Introduction, Age Spectrum, Observation, Tropical Transport in Models, Global Mean Age in Models, Source-Transport Covariance, HSCT "ANOY" Tracer Distributions, and Summary and Conclusions.

Description: A global lightning model that includes diurnal and annual lightning variation, and total flash density versus latitude for each major land and ocean, has been used as the basis for simulating the global electric circuit charging rate. A particular objective has been to reconcile the difference in amplitude ratios [AR=(max-min)/mean] between global lightning diurnal variation (AR approx. = 0.8) and the diurnal variation of typical atmospheric potential gradient curves (AR approx. = 0.35). A constraint on the simulation is that the annual mean charging current should be about 1000 A. The global lightning model shows that negative ground flashes can contribute, at most, about 10-15% of the required current. For the purpose of the charging rate simulation, it was assumed that each ground flash contributes 5 C to the charging process. It was necessary to assume that all electrified clouds contribute to charging by means other than lightning, that the total flash rate can serve as an indirect indicator of the rate of charge transfer, and that oceanic electrified clouds contribute to charging even though they are relatively inefficient in producing lightning. It was also found necessary to add a diurnally invariant charging current component. By trial and error it was found that charging rate diurnal variation curves in Universal time (UT) could be produced with amplitude ratios and general shapes similar to those of the potential gradient diurnal variation curves measured over ocean and arctic regions during voyages of the Carnegie Institute research vessels.

Description: This year repeat elevation surveys in the southern half of Greenland were made using the Airborne Topographic Mapper (ATM). The intent of these surveys is to compare present elevations to those measured in 1993 and determine the magnitude and spatial distributions of thickening and thinning rates. In order to effectively interpret any observed changes, it is important to understand the processes that affect these changes. Moreover, because the surveys are made over a brief period (2-4 weeks) during the spring or summer, it is also important to understand the effects of seasonal and interannual elevation variability, in relation to the timing of these surveys. Toward that end we are examining data from weather stations along the coast of Greenland along with data from GC-Net automatic weather stations (AWS's) on the ice sheet. The objectives are to assess: a) the importance of the timing of the flights in relation to natural processes that affect surface heights, namely accumulation and melt, and b) the temperature characteristics of the region in the five years that separated the two sets of surveys (1993-1998), in relation to the past 19 years.

Description: A distributed satellite formation, modeled as an arbitrary number of fully connected nodes in a network, could be controlled using a decentralized controller framework that distributes operations in parallel over the network. For such problems, a solution that minimizes data transmission requirements, in the context of linear-quadratic-Gaussian (LQG) control theory, was given by Speyer. This approach is advantageous because it is non-hierarchical, detected failures gracefully degrade system performance, fewer local computations are required than for a centralized controller, and it is optimal with respect to the standard LQG cost function. Disadvantages of the approach are the need for a fully connected communications network, the total operations performed over all the nodes are greater than for a centralized controller, and the approach is formulated for linear time-invariant systems. To investigate the feasibility of the decentralized approach to satellite formation flying, a simple centralized LQG design for a spacecraft orbit control problem is adapted to the decentralized framework. The simple design uses a fixed reference trajectory (an equatorial, Keplerian, circular orbit), and by appropriate choice of coordinates and measurements is formulated as a linear time-invariant system.

Description: The Microwave Anisotropy Probe (MAP) is a follow-on to the Differential Microwave Radiometer (DMR) instrument on the Cosmic Background Explorer (COBE) spacecraft. The MAP spacecraft will perform its mission, studying the early origins of the universe, in a Lissajous orbit around the Earth-Sun L(sub 2) Lagrange point. Due to limited mass, power, and financial resources, a traditional reliability concept involving fully redundant components was not feasible. This paper will discuss the redundancy philosophy used on MAP, describe the hardware redundancy selected (and why), and present backup modes and algorithms that were designed in lieu of additional attitude control hardware redundancy to improve the odds of mission success. Three of these modes have been implemented in the spacecraft flight software. The first onboard mode allows the MAP Kalman filter to be used with digital sun sensor (DSS) derived rates, in case of the failure of one of MAP's two two-axis inertial reference units. Similarly, the second onboard mode allows a star tracker only mode, using attitude and derived rate from one or both of MAP's star trackers for onboard attitude determination and control. The last backup mode onboard allows a sun-line angle offset to be commanded that will allow solar radiation pressure to be used for momentum management and orbit stationkeeping. In addition to the backup modes implemented on the spacecraft, two backup algorithms have been developed in the event of less likely contingencies. One of these is an algorithm for implementing an alternative scan pattern to MAP's nominal dual-spin science mode using only one or two reaction wheels and thrusters. Finally, an algorithm has been developed that uses thruster one shots while in science mode for momentum management. This algorithm has been developed in case system momentum builds up faster than anticipated, to allow adequate momentum management while minimizing interruptions to science. In this paper, each mode and algorithm will be discussed, and simulation results presented.

Description: The National Aeronautics and Space Administration (NASA) Goddard Space Flight Center (GSFC) has proposed a set of spacecraft flying in close formation around the Earth in order to measure the behavior of the auroras. The mission, named Auroral Lites, consists of four spacecraft configured to start at the vertices of a tetrahedron, flying over three mission phases. During the first phase, the distance between any two spacecraft in the formation is targeted at 10 kilometers (km). The second mission phase is much tighter, requiring satellite interrange spacing targeted at 500 meters. During the final phase of the mission, the formation opens to a nominal 100-km interrange spacing. In this paper, we present the strategy employed to initialize and model such a close formation during each of these phases. The analysis performed to date provides the design and characteristics of the reference orbit, the evolution of the formation during Phases I and II, and an estimate of the total mission delta-V budget. AI Solutions' mission design tool, FreeFlyer(R), was used to generate each of these analysis elements. The tool contains full force models, including both impulsive and finite duration maneuvers. Orbital maintenance can be fully modeled in the system using a flexible, natural scripting language built into the system. In addition, AI Solutions is in the process of adding formation extensions to the system facilitating mission analysis for formations like Auroral Lites. We will discuss how FreeFlyer(R) is used for these analyses.

Description: The small expendable deployable system and tether satellite system programs did not have a uniform written criteria for tethers. The JSC safety panel asked what criteria was used to design the tethers. Since none existed, a criteria was written based on past experience for future tether programs.

Description: The Propulsive Small Expendable Deployer System (ProSEDS) space experiment will demonstrate the use of an electrodynamic tether propulsion system. The flight experiment is a precursor to the more ambitious electrodynamic tether upper stage demonstration mission which will be capable of orbit raising, lowering and inclination changes-all using electrodynamic thrust. ProSEDS, which is planned to fly in 2000, will use the flight proven Small Expendable Deployer System (SEDS) to deploy a tether (5 km bare wire plus 15 km spectra) from a Delta II upper stage to achieve approx. 0.4N drag thrust, thus deorbiting the stage. The experiment will use a predominantly "bare" tether for current collection in lieu of the endmass collector and insulated tether approach used on previous missions. ProSEDS will utilize tether-generated current to provide limited spacecraft power. In addition to the use of this technology for orbit transfer and upper stages, it may also be an attractive option for future missions to Jupiter and any other planetary body with a magnetosphere.

Description: This report gives an overall view of the CORE program at Goddard Space Flight Center (GSFC). It summarizes the different CORE sessions and gives information about the technical staff. The outlook summarizes the evolution of the different CORE programs.

Description: Airborne laser-altimeter flight lines from 1993 over southern parts of the ice sheet were resurveyed with almost complete repeat coverage. In 1993 and 1994, NASA surveyed the entire Greenland ice sheet by airborne laser altimeter, obtaining surface-elevation profiles with root mean square (rms) accuracies of 10 cm or better (Krabill 1995) along flight lines that crossed all the major catchment basins. In 1998, the ten flight lines flown in 1993 in the south of Greenland were resurveyed with about 99% repeat coverage; flight lines in the north will be resurveyed in 1999. Additional flights in 1998 were over glaciers, identified by E. Rignot, where existing SAR data give information on ice motion.

Description: Many modem and ancient carbonate deposits around the world have been recognized as microbial buildups or microbialites. Ancient microbialite structures have been divided into two basic categories based on their internal fabric or texture. They include stromatolites which have a predominantly laminated internal fabric and thrombolites which have an open-porous clotted fabric, that lacks laminae. The origin of these two basic microbial fabrics is still being debated in the literature. Understanding the origin and the various microorganisms involved in forming these modem fabrics is the key to the interpretation of similar fabrics in ancient and possibly Martian rocks. Therefore, detailed studies are needed on the microbiological makeup and origin of the fabrics in modem microbialites. Such studies may serve as analogs for ancient and Martian microbialites in the future. The purpose of this study is to examine the textures and carbon isotopic signatures of the following modem microbialites from the Bahamas: 1) a modem subtidal microbialite from Iguana Cay, Bahamas and 2) a modem microbial mat (stromatolite) from a hypersaline pond on Lee Stocking Island, Bahamas.

Description: Life at the surface of the Earth, over the last 400 m.y., evolved under conditions of decreased short-wave radiation (i.e., ultraviolet) relative to solar output due to absorption and scattering by constituents (e.g., ozone, water vapor, aerosols) in the upper atmosphere. However, a significant amount of ultraviolet radiation in the range from 280-320 nm, known as ultraviolet-B radiation, reaches the Earth's surface and has sufficient energy to be damaging to biologic tissue. Natural fluctuations in atmospheric constituents (seasonal variation, volcanic eruptions, etc.), changes in the orbital attitude of the Earth (precession, axial tilt, orbital eccentricity), and long-term solar variability contribute to changes in the total amount of ultraviolet radiation reaching the surface of the Earth, and thus, the biosphere. More recently, the atmospheric release of commercial propellants and refrigerants, known as chlorofluorocarbons (CFCs), has contributed to a significant depletion in naturally occurring ozone in the stratosphere. Thus, decreased stratospheric ozone has resulted in an increased UV-B flux at the Earth's surface which may have profound effects on terrestrial and marine organisms. In this study, we are investigating the effects of differing solar UV-B fluxes on alfalfa (Medicago sativa L.), an important agricultural crop. A long-term goal of this research is to develop spectral signatures to detect plant response to increased UV-B radiation from remote sensor platforms.

Description: The Environmental Research Aircraft and Sensor Technology (ERAST) program is a NASA initiative that seeks to demonstrate the application of cost-effective aircraft and sensor technology to private commercial ventures. In 1997-98, a series of flight-demonstrations and image acquisition efforts were conducted over the Hawaiian Islands using a remotely-piloted solar- powered platform (Pathfinder) and a fixed-wing piloted aircraft (Navajo) equipped with a Kodak DCS450 CIR (color infrared) digital camera. As an ERAST Science Team Member, I defined a set of flight lines over the largest coffee plantation in Hawaii: the Kauai Coffee Company's 4,000 acre Koloa Estate. Past studies have demonstrated the applications of airborne digital imaging to agricultural management. Few studies have examined the usefulness of high resolution airborne multispectral imagery with 10 cm pixel sizes. The Kodak digital camera integrated with ERAST's Airborne Real Time Imaging System (ARTIS) which generated multiband CCD images consisting of 6 x 106 pixel elements. At the designated flight altitude of 1,000 feet over the coffee plantation, pixel size was 10 cm. The study involved the analysis of imagery acquired on 5 March 1998 for the detection of anomalous reflectance values and for the definition of spectral signatures as indicators of tree vigor and treatment effectiveness (e.g., drip irrigation; fertilizer application).

Description: Research has established the importance of global tropospheric wind measurements for large scale improvements in numerical weather prediction. In addition, global wind measurements provide data that are fundamental to the understanding and prediction of global climate change. These tasks are closely linked with the goals of the NASA Earth Science Enterprise and Global Climate Change programs. NASA Goddard has been actively involved in the development of direct detection Doppler lidar methods and technologies to meet the wind observing needs of the atmospheric science community. A variety of direct detection Doppler wind lidar measurements have recently been reported indicating the growing interest in this area. Our program at Goddard has concentrated on the development of the edge technique for lidar wind measurements. Implementations of the edge technique using either the aerosol or molecular backscatter for the Doppler wind measurement have been described. The basic principles have been verified in lab and atmospheric lidar wind experiments. The lidar measurements were obtained with an aerosol edge technique lidar operating at 1064 nm. These measurements demonstrated high spatial resolution (22 m) and high velocity sensitivity (rms variances of 0.1 m/s) in the planetary boundary layer (PBL). The aerosol backscatter is typically high in the PBL and the effects of the molecular backscatter can often be neglected. However, as was discussed in the original edge technique paper, the molecular contribution to the signal is significant above the boundary layer and a correction for the effects of molecular backscatter is required to make wind measurements. In addition, the molecular signal is a dominant source of noise in regions where the molecular to aerosol ratio is large since the energy monitor channel used in the single edge technique measures the sum of the aerosol and molecular signals. To extend the operation of the edge technique into the free troposphere we have developed a variation of the edge technique called the double edge technique. In this paper a ground based aerosol double edge lidar is described and the first measurements of wind profiles in the free troposphere obtained with this lidar will be presented.

Description: Earth's thin layer of soil is a fragile resource, made up of minerals, organic materials, air, water, and billions of living organisms. Soils plays a variety of critical roles that sustain life on Earth. If we think about soil, we tend to see it first as the source of most of the food we eat and the fibers we use, such as wood and cotton. Few students realize that soils also provide the key ingredients to many of the medicines (including antibiotics), cosmetics, and dyes that we use. Fewer still understand the importance of soils in integrating, controlling, and regulating the movement of air, water, materials, and energy between the hydrosphere, lithosphere, atmosphere, and biosphere. Because soil sustains life, it offers both a context and a natural laboratory for investigating these interactions. The enclosed poster, which integrates soil profiles with typical landscapes in which soils form, can also help students explore the interrelationships of Earth systems and gain an understanding of our soil resources. The poster, produced jointly by the American Geological Institute and the Soil Science Society of America, aims to increase awareness of the importance of soil, as does the GLOBE (Global Learning and Observations To Benefit the Environment) Program. Vice President Al Gore instituted the GLOBE Program on Earth Day of 1993 to increase environmental awareness of individuals throughout the world, contribute to a better scientific understanding of the Earth, and help all students reach higher levels of achievement in science and mathematics. GLOBE functions as a partnership between scientists, students, and teachers in which scientists design protocols for specific measurements they need for their research that can be performed by K-12 students. Teachers are trained in the GLOBE protocols and teach them to their students. Students make the measurements, enter data via the Internet to a central data archive, and the data becomes available to scientists and the general community. Students benefit by having a "hands-on"experience in science, math, and technology, using their local environment as a learning laboratory, as well as contact with scientists and other students around the world. Soil investigations have become an essential component of GLOBE. The protocols that have been developed so far within the GLOBE program include GPS Location, Atmosphere/Climate, Soil Characterization, Soil Moisture and Temperature, Land Cover/Biometry, Hydrology, and Satellite Image Classification. For the GLOBE Soil Characterization Protocol, students explore the physical. chemical, and morphological properties of the soil at their study site. They are asked to dig a pit or use an auger to about 1 meter at at least 2 sites.

Description: Geosynchronous Synthetic Aperture Radar (GeoSAR) is a consortium project consisting of The Jet Propulsion Laboratory (JPL), Calgis (a small GIS company based in Fresno, CA) and the California Department of Conservation with funding provided by Defense Research Projects Agency (DARPA) started in November 1996. The two main objectives of the GeoSAR Program are: 1) To develop a state of the art dual frequency interferometric radar mapping instrument capable of mapping the true ground surface height beneath the vegetation canopy; and 2) To transition this mapping technology to a commercial company, Calgis. JPL, the technical lead, has the following program deliverables at program completion in November 1999 include radar design and radar hardware for X-band (3 cm) and P-band (83 cm) radars, processor software, hardware and documentation, and calibrated X-band radar.

Description: Data in the wavelength range 0.545 - 1.652 microns from the Moderate Resolution Imaging Spectroradiometer (MODIS), to be launched aboard the Earth Observing System (EOS) Terra in the fall of 1999, will be used to map daily global snow cover at 500m resolution. However, during darkness, or when the satellite's view of the surface is obscured by cloud, snow cover cannot be mapped using MODIS data. We show that during these conditions, it is possible to supplement the MODIS product by mapping the snow cover using passive microwave data from the Special Sensor Microwave Imager (SSM/I), albeit with much poorer resolution. For a 7-day time period in March 1999, a prototype MODIS snow-cover product was compared with a prototype MODIS-SSM/I product for the same area in the mid-western United States. The combined MODIS-SSM/I product mapped 9% more snow cover than the MODIS-only product.

Description: The Laser Vegetation Imaging Sensor (LVIS) is an airborne, scanning laser altimeter designed and developed at NASA's Goddard Space Flight Center. LVIS operates at altitudes up to 10 km above ground, and is capable of producing a data swath up to 1000 m wide nominally with 25 m wide footprints. The entire time history of the outgoing and return pulses is digitized, allowing unambiguous determination of range and return pulse structure. Combined with aircraft position and attitude knowledge, this instrument produces topographic maps with decimeter accuracy and vertical height and structure measurements of vegetation. The laser transmitter is a diode-pumped Nd:YAG oscillator producing 1064 nm, 10 nsec, 5 mJ pulses at repetition rates up to 500 Hz. LVIS has recently demonstrated its ability to determine topography (including sub-canopy) and vegetation height and structure on flight missions to various forested regions in the U.S. and Central America. The LVIS system is the airborne simulator for the Vegetation Canopy Lidar (VCL) mission (a NASA Earth remote sensing satellite due for launch in 2000), providing simulated data sets and a platform for instrument proof-of-concept studies. The topography maps and return waveforms produced by LVIS provide Earth scientists with a unique data set allowing studies of topography, hydrology, and vegetation with unmatched accuracy and coverage.

Description: The authors propose a grouped threshold method for scene identification in Advanced Very High Resolution Radiometer imagery that may contain clouds, fire, smoke, or snow. The philosophy of the approach is to build modules that contain groups of spectral threshold tests that are applied concurrently, not sequentially, to each pixel in an image. The purpose of each group of tests is to identify uniquely a specific class in the image, such as smoke. A strength of this approach is that insight into the limits used in the threshold tests may be gained through the use of radiative transfer theory. Methodology and examples are provided for two different scenes, one containing clouds, forest fires, and smoke; and the other containing clouds over snow in the central United States. For both scenes, a limited amount of supporting information is provided by surface observers.

Description: Estimates of the effect of pulse stretching on satellite laser altimetry in particular the Geoscience Laser Altimeter System (GLAS), by cloud multiple scattering were made from an analytical method and from Monte Carlo calculations. The path delay of the return pulse was found to be largest for low-level clouds with particle radii (3-20 microns). The magnitude of the path delay was affected by several factors including cloud height, cloud optical depth, cloud particle size, particle shape, and receiver field of view. Polar aerosol and Rayleigh scattering usually added less than 1 cm to the overall path delay. Path delay estimates for all cloud conditions would be less if a simple Gaussian fit of the return pulse peak were used to measure the pulse's centroid. However, care must be taken in determining the centroid as factors such as pulse width, surface slope and the fitting method used will affect the estimate. A planned application for laser altimetry is high precision monitoring of the height change of polar ice sheets. In the absence of a correction algorithm, the required GLAS altimetry accuracies will not be achieved unless atmospheric channel information is used to remove profiles that are likely to be heavily contaminated by multiple scattering.

Description: Following the launch of the Earth Observing System first morning (EOS-AM1) satellite, daily, global snow-cover mapping will be performed automatically at a spatial resolution of 500 m, cloud-cover permitting, using Moderate Resolution Imaging Spectroradiometer (MODIS) data. A technique to calculate theoretical accuracy of the MODIS-derived snow maps is presented. Field studies demonstrate that under cloud-free conditions when snow cover is complete, snow-mapping errors are small (less than 1%) in all land covers studied except forests where errors are greater and more variable. The theoretical accuracy of MODIS snow-cover maps is largely determined by percent forest cover north of the snowline. Using the 17-class International Geosphere-Biosphere Program (IGBP) land-cover maps of North America and Eurasia, the Northern Hemisphere is classified into seven land-cover classes and water. Snow-mapping errors estimated for each of the seven land-cover classes are extrapolated to the entire Northern Hemisphere for areas north of the average continental snowline for each month. Average monthly errors for the Northern Hemisphere are expected to range from 5 - 10%, and the theoretical accuracy of the future global snow-cover maps is 92% or higher. Error estimates will be refined after the first full year that MODIS data are available.

Description: As increasingly complex scientific and environmental observation spacecraft are deployed, the burden on the downlink assets, and ground-based systems complexity and cost is becoming a major problem. Already, the limitations of communications bandwidth and processing throughput limit the science data gathering, both in volume and in rate. This poses a dilemma to the scientist experimenter forcing choices between data collection and bandwidth/processing/archiving. Advances in ground based processing and space-to-Earth links have fallen behind the requirements for observation data, at increasing rates, over the last few decades. As NASA achieves its 40th anniversary, the ability to observe and capture phenomena of theoretical and practical interest to life on Earth far outstrips the ability to transfer, process, or store these data. NASA recognizes the need to invest on technological advancements that will enable both the space and ground systems to address the limitations. Spacecraft onboard computing power is a clear one. The capability of creating data products onboard the spacecraft adds a new level of flexibility to address the more demanding observation needs. Current spacecraft computing power is limited and incapable of addressing the needs of the new generation of observation satellites because extensive onboard data processing is required. Traditional spacecraft architectures only collect, package, and transmit to Earth the data acquired by multiple instruments. Conversely, the experience on developing ground data systems shows the need for high performance computing systems to process and create information from the instrumentation data. The expectation is that supercomputing technology is required to enable spacecraft to create information onboard. Moving supercomputing capability onboard spacecraft requires an approach that considers an integrated data architecture. Otherwise, it may simply convert a compute-bound problem into a communications bound problem, as has been shown numerous times in the context of massively parallel architectures. What is left to determine are the technologies that will enable spacecraft high performance computing.

Description: The upcoming generation of laser altimeters record the interaction of emitted laser radiation with terrestrial surfaces in the form of a digitized waveform. We model these laser altimeter return waveforms as the sum of the reflections from individual surfaces within laser footprints, accounting for instrument-specific properties. We compare over 1000 modeled and recorded waveform pairs using the Pearson correlation. We show that we reliably synthesize the vertical structure information for vegetation canopies contained in a medium-large diameter laser footprint from a high-resolution elevation data set.

Description: This study investigates the use of H2, p-synthesis, and mixed H2/mu methods to construct full-order controllers and optimized controllers of fixed dimensions. The benchmark problem definition is first extended to include uncertainty within the controller bandwidth in the form of parametric uncertainty representative of uncertainty in the natural frequencies of the design model. The sensitivity of H2 design to unmodelled dynamics and parametric uncertainty is evaluated for a range of controller levels of authority. Next, mu-synthesis methods are applied to design full-order compensators that are robust to both unmodelled dynamics and to parametric uncertainty. Finally, a set of mixed H2/mu compensators are designed which are optimized for a fixed compensator dimension. These mixed norm designs recover the H, design performance levels while providing the same levels of robust stability as the u designs. It is shown that designing with the mixed norm approach permits higher levels of controller authority for which the H, designs are destabilizing. The benchmark problem is that of an active tendon system. The controller designs are all based on the use of acceleration feedback.

Description: The zonal mean eddy heat flux is directly proportional to the wave activity that propagates from the troposphere into the stratosphere. This quantity is a simple eddy diagnostic which is easily calculated from conventional meteorological analyses. Because this "wave driving" of the stratosphere has a strong impact on the stratospheric temperature, it is necessary to compare the impact of the flux with respect to stratospheric radiative changes caused by greenhouse gas changes. Hence, we must understand the precision and accuracy of the heat flux derived from our global meteorological analyses. Herein, we quantify the stratospheric heat flux using five different meteorological analyses, and show that there are 30% differences between these analyses during the disturbed conditions of the northern hemisphere winter. Such large differences result from the planetary differences in the stationary temperature and meridional wind fields. In contrast, planetary transient waves show excellent agreement amongst these five analyses, and this transient heat flux appears to have a long term downward trend.

Description: We examine concurrent upper tropospheric measurements of CN (diameter greater than 4 nm). NO, and NO(Y) during the SONEX Experiment over the North Atlantic (Oct.-Nov., 1997). Elevated CN and NO(Y) concentrations observed in the upper troposphere are attributed largely to enhancements in convective outflows. We estimate that less than 7% of observed high-CN plumes (greater than 10000 /cc) may be attributed to aircraft emissions. Dilution of high-CN convective and aircraft plumes appears to be much more rapid than losses of NO(X) and CN by oxidation and coagulation, respectively, and accounts for much of observed CN concentrations. When taking into account of different time scales against dilution for observable aircraft and convective high-CN plumes (estimated to be 1:4), the contribution by aircraft emissions to CN concentrations is significant, about 20% of the convective source. We find no evidence that particle formation in convective plumes is limited by OH oxidation of SO2.

Description: This study examines a unique data set returned by IMP8 and Geotail on January 29, 1995 during a substorm which resulted in the ejection of a plasmoid. The two spacecraft (s/c) were situated in the north lobe of the tail and both observed a traveling compression region (TCR). From single s/c observations only the length of the plasmoid in X and an estimate of its height in Z can be determined. However, we show that dual s/c measurements of TCRs can be used to model all three dimensions of the underlying plasmoid and to estimate of its rate of expansion or contraction. For this event plasmoid dimensions of Delta(X) approximates 18, Delta(Y) approximates 30, and Delta(Z) approximates 10 R(sub e) are inferred from the IMP8 and Geotail lobe magnetic field measurements. The earthward end of the plasmoid was inferred to be near the mean location of the near-earth neutral line, X approximates -26 R(sub e). Its center was underneath IMP 8 at X approximates -34 R(sub e) and its tailward end appeared to be near X approximates -44 R(sub e). Furthermore, a factor of approximately 2 increase in the amplitude of the TCR occurred in the 1.5 min it took to move from IMP 8 to Geotail. Modeled using conservation of the magnetic flux, this increase in lobe compression implies that the underlying plasmoid was expanding at a rate of approximately 140 km/s. Such an expansion is comparable to recently reported V(sub y) speeds in "young" plasmoids in this region of the tail. Finally, the Geotail measurements indicate that a reconfiguration of the lobe magnetic field closely followed the ejection of the plasmoid which moved magnetic flux tubes into the wake behind the plasmoid where they would convect into the near-earth neutral line and reconnect.

Description: Vertical mixing ratio profiles of four relatively long-lives gases, HCN, C2H2, CO, and C2H6, have been retrieved from 0.01/cm resolution infrared solar occultation spectra recorded between latitudes of 5.3degN and 31.4degN. The observations were obtained by the Atmospheric Trace Molecule Spectroscopy (ATMOS) Fourier transform spectrometer during the Atmospheric Laboratory for Applications and Science (ATLAS) 3 shuttle flight, 3-12 November 1994. Elevated mixing ratios below the tropopause were measured for these gases during several of the occultations. The positive correlations obtained between the simultaneously measured mixing ratios suggest that the enhancements are likely the result of surface emissions, most likely biomass burning and/or urban industrial activities, followed by common injection via deep convective transport of the gases to the upper troposphere. The elevated levels of HCN may account for at least part of the "missing NO," in the upper troposphere. Comparisons of the observations with values measured during a recent aircraft campaign are presented.

Description: Prior to 1991, major warmings (defined by increasing zonal mean temperatures and zonal mean easterly winds from 60degN to the pole at 10 hPa) typically occurred approximately once every two Arctic winters; a major warming in mid-Dec. 1998 was the first since Feb. 1991. The Dec. 1998 warming was also the second earliest on record. The earliest, and the only other major warming on record before the end of Dec. was in early Dec 1987; prior to that, the earliest was in late Dec./early Jan. 1984-85. The 1984-85 and 1987 warmings resulted in the warmest and weakest lower stratospheric polar vortices in the 20 years before 1998-99. Fig. 1 compares temperatures and vortex strength in 1998-99 with those in the previous 20 years, using the US National Center for Environmental Prediction (NCEP) record; 1987-88 and 1984-85 are also highlighted. The Dec. 1998 warming had a more pronounced effect on mid-stratospheric temperatures than the Dec. 1987 warming (Fig. 1a), although smaller than that of warmings later in winter (e.g., 1984-85). 10-hPa temperatures fell well below average again in late Jan. 1999 and remained unusually low until an early final warming began in late Feb. 840 K PV gradients (Fig. 1c) set a record minimum in Jan. 1999, but were near average in Feb before the final warming. The effect of the Dec. 1998 warming on lower stratospheric temperatures was comparable to that of other major warmings; there was a brief period of record-high minimum 46-hPa temperatures in early Jan 1999 (Fig. 1b), and temperatures then fell to near average for a short period in mid-Feb. Lower stratospheric PV gradients were the weakest on record during the 1998-99 winter (Fig. 1d). The evolution of the vortex and minimum temperatures during 1998-99 was remarkably similar to that during 1987-88, the only previous year when a major warming was observed before the end of Dec.

Description: Analysis of a time series of European Remote Sensing Satellite (ERS)-1 and -2, RADARSAT ScanSAR synthetic aperture radar (SAR) and Landsat images from 1973 to 1998, shows daily to interannual changes in Hofsjokull, a 923 sq km ice cap in central Iceland. A digital elevation model of Hofsjokull was constructed using interferometry, and then SAR backscatter coefficient (d) was plotted with elevation, and air temperature along a transect across the ice cap. Most of the a' changes measured along the transect are caused by a change in the state (frozen or thawed) of the surficial snow or ice when air temperature rises above or below about -5 to O C. Seasonal (sigma)deg patterns are identified in a 4-year time series of 57 ERS-1 and -2 images. In addition, June 1997 ScanSAR images display rapid changes in brightness that are tied closely to daily meteorological events. SAR and Landsat data were also used to measure changes in the areal extent of Hofsjokull, from 1973 to 1997, and to locate (sigma)deg and reflectance boundaries that relate to the glacier facies. Late-summer 1997 (sigma)deg and reflectance boundaries agree and are coincident with the approximate location of the fim line, and the January 1998 position of the equilibrium line as determined from ERS-2 data.

Description: Aerosol particles originate from man-made sources such as urban/industrial activities,rurning associated with land use processes, wind-blown dust, and natural sources. Their interaction with sunlight and their effect on cloud microphysics forms a major uncertainty in predicting climate change. Furthermore, the lifetime of only a few days causes high spatial variability in aerosol optical and radiative properties that requires global observations from space. Remote sensing of aerosol properties from space is reviewed both for present and planned national and international satellite sensors. Techniques that are being used to enhance our ability to characterize the global distribution of aerosol properties include well-calibrated multispectral radiometers, multispectral polarimeters, and multi-angle spectroradiometers. Though most of these sensor systems rely primarily on visible to mid-infrared spectral channels, the availability of of thermal channels to aid in cloud screening is an important additional piece of information that is not always incorporated into the sensor design. In this paper we describe the various satellite sensor systems being developed by Europe, Japan, and the U.S., and highlight the advantages and disadvantages of each of these systems for aerosol applications. An important underlying theme is that the remote sensing of aerosol properties, especially aerosol size distribution and single scattering albedo, is exceedingly difficult. As a consequence, no one sensor system is capable of providing totally unambiguous information, and hence a careful intercomparison of derived products from different sensors, together with a comprehensive network of ground-based sun-photometer and sky radiometer systems, are required to advance our quantitative understanding of global aerosol characteristics.

Description: Various configurations of a scanning satellite instrument are simulated by sampling realistic fields of nitrous oxide. Synoptic grids are computed from the resulting simulated orbital data and compared to the original sampled data fields. Results are compared with those obtained by flying a simulated satellite over low-resolution fields and fields that are static in time. Although increasing the number of instrument scan positions does provide more information along an orbital swath, using more than three to five scan positions does not significantly increase the accuracy of global synoptic grids using the gridding techniques described here.

Description: Radar data from vegetated land surfaces depend on many structural and compositional parameters describing the terrain. Because early, noninterferometric radar systems usually constituted an insufficient observation set from which to estimate parameters of the terrain, statistical regression techniques were used which incorporated some level of apriori knowledge or field measurements. With the advent of radar interferometry and polarimetric interferometry, potentially at multiple baselines, the observation set is now approaching that required to quantitatively estimate the parameters describing a vegetated land surface. Quantitative estimation entails formulating a physical scattering model relating the radar observations to the vegetation and surface parameters on which they depend. This paper describes the physics of candidate scattering models, and shows how the models determine the estimable parameter set. It also indicates the measurement accuracy of parameters such as vegetation height, height-to-base-of-live-crown, and surface topography with multibaseline polarimetric interferometry.

Description: This viewgraph presentation focuses on the past, present, and future space parts environment. The past environment was characterized by long lead time flagship missions having substantial support from NASA and DOD. The future environment is characterized by many BFC missions, short development cycles, smaller projects and shorter part delivery schedules. These ideas are elaborated upon.

Description: In this paper a fuzzy classification procedure is applied to polarimetric radar measurements, and street pixels are detected. These data are successively grouped into consistent roads by means of a dynamic programming approach based on the fuzzy membership function values. Further fusion of the 2D road network extracted and 3D TOPSAR measurements provides a powerful way to analyze urban infrastructures.

Description: An autonomous GPS array is being implemented in the Transantarctic Mountains, sponsored by NSF and NASA, for the purpose of measuring uplift resulting from post-glacial rebound (PGR). The rebound of the solid earth due to unloading of ice since the Last Glacial Maximum is expected to dominate the measured uplift for most of West Antarctica, dwarfing the signals due to present-day ice sheet mass balance changes and tectonic motion, as long as mantle viscosity is greater than about 10(exp 20) Pa-s. Predicted uplift patterns have been calculated for a range of model scenarios, which illustrate how the uplift pattern might distinguish between different-sized ice sheets and deglaciation histories as represented by the competing models. The scenarios considered by James and Ivins (1998) include ICE-3G, CLIMAP and a variation of the CLIMAP model by Denton et al. For these models, peak uplift rates occur in the Transantarctic Mountains, and differences between models is often large there. Thus, the Transantarctic Mountains are an ideal place to obtain uplift measurements to constrain deglaciation models.

Description: An autonomous spacecraft must balance long-term and short-term considerations. It must perform purposeful activities that ensure long-term science and engineering goals are achieved and ensure that it maintains positive resource margins. This requires planning in advance to avoid a series of shortsighted decisions that can lead to failure, However, it must also respond in a timely fashion to a somewhat dynamic and unpredictable environment. Thus, spacecraft plans must often be modified due to fortuitous events such as early completion of observations and setbacks such as failure to acquire a guidestar for a science observation. This paper describes the use of iterative repair to support continuous modification and updating of a current working plan in light of changing operating context.

Description: In this paper, a multiscale approach is introduced to classify the Japanese Research Satellite-1 (JERS-1) mosaic image over the Central African rainforest. A series of texture maps are generated from the 100 m mosaic image at various scales. Using a quadtree model and relating classes at each scale by a Markovian relationship, the multiscale images are classified from course to finer scale. The results are verified at various scales and the evolution of classification is monitored by calculating the error at each stage.

Description: In the execution of this proposal, we will first examine current and developing spacecraft materials and evaluate their ability to attenuate adverse biological mutational events in mammalian cell systems and reduce the rate of cancer induction in mice harderian glands as a measure of their protective qualities. The HZETRN code system will be used to generate a database on GCR attenuation in each material. If a third year of funding is granted, the most promising and mission-specific materials will be used to study the impact on mission cost for a typical Mars mission scenario as was planned in our original two year proposal at the original funding level. The most promising candidate materials will be further tested as to their transmission characteristics in Fe and Si ion beams to evaluate the accuracy of the HZETRN transmission factors. Materials deemed critical to mission success may also require testing as well as materials developed by industry for their radiation protective qualities (e.g., Physical Sciences Inc.) A study will be made of designing polymeric materials and composite materials with improved radiation shielding properties as well as the possible improvement of mission-specific materials.

Description: In this paper, we report the results of our recent research on the application of a multiprocessor Cray T916 supercomputer in modeling super-thermal electron transport in the earth's magnetic field. In general, this mathematical model requires numerical solution of a system of partial differential equations. The code we use for this model is moderately vectorized. By using Amdahl's Law for vector processors, it can be verified that the code is about 60% vectorized on a Cray computer. Speedup factors on the order of 2.5 were obtained compared to the unvectorized code. In the following sections, we discuss the methodology of improving the code. In addition to our goal of optimizing the code for solution on the Cray computer, we had the goal of scalability in mind. Scalability combines the concepts of portabilty with near-linear speedup. Specifically, a scalable program is one whose performance is portable across many different architectures with differing numbers of processors for many different problem sizes. Though we have access to a Cray at this time, the goal was to also have code which would run well on a variety of architectures.

Description: This project focuses on the adaptation of human populations to their environments from prehistoric times to the present. It emphasizes interdisciplinary research to develop ecological baselines through the use of remotely sensed imagery, in situ field work, and the modeling of human population dynamics. It utilizes cultural and biological data from dated archaeological sites to assess the subsistence and settlement patterns of human societies in response to changing climatic and environmental conditions. The utilization of remote sensing techniques in archaeology is relatively new, exciting, and opens many doors.

Description: The International Space Station has been in development since 1984, and has recently begun on orbit assembly. Most of the hardware for the Space Station has been manufactured and the rest is well along in design. The major sets of hardware that are still to be developed for Space Station are the pallets and interfacing hardware for resupply of unpressurized spares and scientific payloads. Over the last ten years, there have been numerous starts, stops, difficulties and challenges encountered in this effort. The Space Station program is now entering the beginning of orbital operations. The Program is only now addressing plans to design and build the carriers that will be needed to carry the unpressurized cargo for the Space Station lifetime. Unpressurized carrier development has been stalled due to a broad range of problems that occurred over the years. These problems were not in any single area, but encompassed budgetary, programmatic, and technical difficulties. Some lessons of hindsight can be applied to developing carriers for the Space Station. Space Station teams are now attempting to incorporate the knowledge gained into the current development efforts for external carriers. In some cases, the impacts of these lessons are unrecoverable for Space Station, but can and should be applied to future programs. This paper examines the progress and problems to date with unpressurized carrier development identifies the lessons to be learned, and charts the course for finally accomplishing the delivery of these critical hardware sets.

Description: Two different methods for retrieving Upper Tropospheric Humidities (UTH) from the TOVS (TIROS Operational Vertical Sounder) instruments aboard NOAA polar orbiting satellites are presented and compared. The first one, from the Environmental Technology Laboratory, computed by J. Bates and D. Jackson (hereafter BJ method), estimates UTH from a simplified radiative transfer analysis of the upper tropospheric infrared water vapor channel at wavelength measured by HIRS (6.3 micrometer). The second one results from a neural network analysis of the TOVS (HIRS and MSU) data developed at, the Laboratoire de Meteorologie Dynamique (hereafter the 3I (Improved Initialization Inversion) method). Although the two methods give very similar retrievals in temperate regions (30-60 N and S), an absolute bias up to 16% appears in the convective zone of the tropics. The two datasets have also been compared with UTH retrievals from infrared radiance measurements in the 6.3 micrometer channel from the geostationary satellite METEOSAT (hereafter MET method). The METEOSAT retrievals are systematically drier than the TOVS-based results by an absolute bias between 5 and 25%. Despite the biases, the spatial and temporal correlations are very good. The purpose of this study is to explain the deviations observed between the three datasets. The sensitivity of UTH to air temperature and humidity profiles is analysed as are the clouds effects. Overall, the comparison of the three retrievals gives an assessment of the current uncertainties in water vapor amounts in the upper troposphere as determined from NOAA and METEOSAT satellites.

Description: Carbonados are porous polycrystalline (with crystal sizes up to 100 micrometer) diamonds. Carbonado is found only in alluvial deposits in Bahia, Brazil and in the Central African Republic (CAR). Alluvial deposit host is 1.5 Ga while the carbonados are between 2.6 - 3.8 Ga. The process of fusing the carbonado microcrystals together is not fully understood, partly due to fact that the origin of these carbonado, is not known. Several modes of origins are proposed for carbonado. First, a crustal origin. Carbonados have a light carbon and helium isotopic signature. They contain an enrichment of the rare-earth elements (REE). Carbonados have tightly trapped atmospheric noble gases and contain an evidence of high He content despite the carbonado expected depletion of He at mantle temperatures. Carbonados have high porosity incompatible with high pressure mantle conditions. Second, a mantle origin is proposed based on similar REE pattern to kimberlites. The presence of nitrogen platelet (by IR spectra) indicates high temperature origin and syngenetic inclusions of rutile, ilmenite, and magnetite indicates high pressure and high temperature conditions consistent with mantle origin as well. Third, it is proposed that carbonado diamonds are a result of early impacts into crustal rocks. This is indicated by the rare and controversial occurrence of high pressure diamond polymorph, londsdaleite, frequently found in diamonds from meteorite impact sites, and by observation of planar deformation features, possibly indicating shock events. Finally, it is suggested that carbonados have an extraterrestrial origin, as indicated by a long term annealing based on observation of a zero-phonon line, attributed to paired nitrogen atoms in association with a vacancy.

Description: This study extends a previous investigation on estimating surface soil moisture using the Special Sensor Microwave/Imager (SSM/I) over a grassland region. Although SSM/I is not optimal for soil moisture retrieval, it can under some conditions provide information. Rigorous analyses over land have been difficult due to the lack of good validation data sets. A scientific objective of the Southern Great Plains 1997 (SGP97) Hydrology Experiment was to investigate whether the retrieval algorithms for surface soil moisture developed at higher spatial resolution using truck-and aircraft-based passive microwave sensors can be extended to the coarser resolutions expected from satellite platform. With the data collected for the SGP97, the objective of this study is to compare the surface soil moisture estimated from the SSM/I data with those retrieved from the L-band Electronically Scanned Thinned Array Radiometer (ESTAR) data, the core sensor for the experiment, using the same retrieval algorithm. The results indicated that an error of estimate of 7.81% could be achieved with SSM/I data as contrasted to 2.82% with ESTAR data over three intensive sampling areas of different vegetation regimes. It confirms the results of previous study that SSM/I data can be used to retrieve surface soil moisture information at a regional scale under certain conditions.

Description: Two ascending European Space Agency (ESA) Earth Resources Satellites (ERS)-1/-2 tandem-mode, synthetic aperture radar (SAR) pairs are used to calculate the surface elevation of Hofsjokull, an ice cap in central Iceland. The motion component of the interferometric phase is calculated using the 30 arc-second resolution USGS GTOPO30 global digital elevation product and one of the ERS tandem pairs. The topography is then derived by subtracting the motion component from the other tandem pair. In order to assess the accuracy of the resultant digital elevation model (DEM), a geodetic airborne laser-altimetry swath is compared with the elevations derived from the interferometry. The DEM is also compared with elevations derived from a digitized topographic map of the ice cap from the University of Iceland Science Institute. Results show that low temporal correlation is a significant problem for the application of interferometry to small, low-elevation ice caps, even over a one-day repeat interval, and especially at the higher elevations. Results also show that an uncompensated error in the phase, ramping from northwest to southeast, present after tying the DEM to ground-control points, has resulted in a systematic error across the DEM.

Description: Vertical electron-density (N (sub e)) profiles, deduced from newly-available ISIS-II digital ionospheric topside-sounder data, are used to investigate the "polar-hole" region within the winter, nighttime polar cap ionosphere during solar minimum. The hole region is located around 0200 MLT near the poleward side of the auroral oval. Earlier investigations had revealed very low N (sub e) values in this region (down to 200/cu cm near 300 km). In the present study, such low N, values (approx. 100/cu cm) were only found near the ISIS (International Satellite for Ionospheric Study)-II altitude of 1400 km. The peak ionospheric concentration below the spacecraft remained fairly constant (approx. 10 (exp 5)/cu cm across the hole region but the altitude of the peak dropped dramatically. This peak dropped, surprisingly, to the vicinity of 100 km. These observations suggest that the earlier satellite in situ measurements, interpreted as deep holes in the ionospheric F-region concentration, could have been made during conditions of an extreme decrease in the altitude of the ionospheric N (sub e) peak. The observations, in combination with other data, indicate that the absence of an F-layer peak may be a frequent occurrence at high latitudes.

Description: Lidar altimeter observations of vegetated landscapes provide a time-resolved measure of laser pulse backscatter energy from canopy surfaces and the underlying ground. Airborne lidar altimeter data was acquired using the Scanning Lidar Imager of Canopies by Echo Recovery (SLICER) for a successional sequence of four, closed-canopy, deciduous forest stands in eastern Maryland. The four stands were selected so as to include a range of canopy structures of importance to forest ecosystem function, including variation in the height and roughness of the outer-most canopy surface and the vertical organization of canopy stories and gaps. The character of the SLICER backscatter signal is described and a method is developed that accounts for occlusion of the laser energy by canopy surfaces, transforming the backscatter signal to a canopy height profile (CHP) that quantitatively represents the relative vertical distribution of canopy surface area. The transformation applies an increased weighting to the backscatter amplitude as a function of closure through the canopy and assumes a horizontally random distribution of the canopy components. SLICER CHPs, averaged over areas of overlap where lidar ground tracks intersect, are shown to be highly reproducible. CHP transects across the four stands reveal spatial variations in vegetation, at the scale of the individual 10 m diameter laser footprints, within and between stands. Averaged SLICER CHPs are compared to analogous height profile results derived from ground-based sightings to plant intercepts measured on plots within the four stands. Tbe plots were located on the segments of the lidar ground tracks from which averaged SLICER CHPs were derived, and the ground observations were acquired within two weeks of the SLICER data acquisition to minimize temporal change. The differences in canopy structure between the four stands is similarly described by the SLICER and ground-based CHP results, however a Chi-square test of similarity documents differences that are statistically significant. The differences are discussed in terms of measurement properties that define the smoothness of the resulting CHPs and Lidar Altimeter Measurements of Canopy Structure - Harding et al. canopy properties that may vertically bias the CHP representations of canopy structure. The statistical differences are most likely due to the more noisy character of the ground-based CHPs, especially high in the canopy where ground-based sightings are rare resulting in an underestimate of canopy surface area and height, and to departures from the assumption of horizontal randomness which bias the CHPs toward the observer (upward for SLICER and downward for ground-based CHPs). The results demonstrate that the SLICER observations reliably provide a measure of canopy structure that reveals ecologically interesting structural variations such as those characterizing a successional sequence of closed-canopy, broadleaf forest stands.

Description: We have examined the sea level height tide records at seven tide gauge sites in the region of southcentral Alaska that were affected by the 1964 Prince William Sound earthquake to determine the history of crustal uplift subsequent to the earthquake. There is considerable variation in the behavior depending on the location of the site relative to the 1964 rupture. At Seward, on the eastern side of the Kenai Peninsula we find a slow uplift that is consistent with elastic strain accumulation while at Seldovia and Nikiski on the western side of the Kenai we find a persistent rapid uplift of about 1 cm/yr that most likely represents a long term transient response to the earthquake, but which cannot be sustained over the expected recurrence interval for a great earthquake of several hundred years. Further to the southwest, at Kodiak, we find evidence that the rate of uplift, which is still several mm/yr, has slowed significantly over the past three and a half decades. To the east of the Kenai Peninsula we find subsidence at Cordova and an uncertain behavior at Valdez. At both of these sites there is a mathematically significant time-dependence to the uplift behavior, but the data confirming this time dependence are not as convincing as at Kodiak. At Anchorage, to the north there is little evidence of vertical motion since the earthquake. We compare these long term tide gauge records to recent GPS observations. In general there is reasonable consistency except at Anchorage and Cordova where the GPS measurement indicate somewhat more rapid uplift and subsidence, respectively.

Description: The geographic limits of the South Atlantic Anomaly, as defined by radiation damage, are compared to contours of geomagnetic total field intensity, as defined by the 1995 IGRF, for the present and recent past. The most likely secular variation of the geomagnetic field is estimated and used to extrapolate the geomagnetic field to the year 2100. This indicates that radiation damage to spacecraft and humans in space will likely increase and to cover a much larger geographic area.

Description: We have adopted the transport scenarios used in Part 1 to examine the sensitivity of stratospheric aircraft perturbations to transport changes in our 2-D model. Changes to the strength of the residual circulation in the upper troposphere and stratosphere and changes to the lower stratospheric K(sub zz) had similar effects in that increasing the transport rates decreased the overall stratospheric residence time and reduced the magnitude of the negative perturbation response in total ozone. Increasing the stratospheric K(sub yy) increased the residence time and enhanced the global scale negative total ozone response. However, increasing K(sub yy) along with self-consistent increases in the corresponding planetary wave drive, which leads to a stronger residual circulation, more than compensates for the K(sub yy)-effect, and results in a significantly weaker perturbation response, relative to the base case, throughout the stratosphere. We found a relatively minor model perturbation response sensitivity to the magnitude of K(sub yy) in the tropical stratosphere, and only a very small sensitivity to the magnitude of the horizontal mixing across the tropopause and to the strength of the mesospheric gravity wave drag and diffusion. These transport simulations also revealed a generally strong correlation between passive NO(sub y) accumulation and age of air throughout the stratosphere, such that faster transport rates resulted in a younger mean age and a smaller NO(y) mass accumulation. However, specific variations in K(sub yy) and mesospheric gravity wave strength exhibited very little NO(sub y)-age correlation in the lower stratosphere, similar to 3-D model simulations performed in the recent NASA "Models and Measurements" II analysis. The base model transport, which gives the most favorable overall comparison with inert tracer observations, simulated a global/annual mean total ozone response of -0.59%, with only a slightly larger response in the northern compared to the southern hemisphere. For transport scenarios which gave tracer simulations within some agreement with measurements, the annual/globally averaged total ozone response ranged from -0.45% to -0.70%. Our previous 1995 model exhibited overly fast transport rates, resulting in a global/annually averaged perturbation total ozone response of -0.25%, which is significantly weaker compared to the 1999 model. This illustrates how transport deficiencies can bias model simulations of stratospheric aircraft.

Description: According to the slip partitioning concept, the trench parallel component of relative plate motion in regions of oblique convergence is accommodated by strike-slip faulting in the overriding continental lithosphere. The pattern of postseismic surface deformation due to viscoelastic flow in the lower crust and asthenosphere following a major earthquake on such a fault is modified from that predicted from the conventual elastic layer over viscoelastic halfspace model by the presence of the subducting slab. The predicted effects, such as a partial suppression of the postseismic velocities by 1 cm/yr or more immediately following a moderate to great earthquake, are potentially detectable using contemporary geodetic techniques.

Description: Many spacecraft attitude determination methods use exactly two vector measurements. The two vectors are typically the unit vector to the Sun and the Earth's magnetic field vector for coarse "sun-mag" attitude determination or unit vectors to two stars tracked by two star trackers for fine attitude determination. TRIAD, the earliest published algorithm for determining spacecraft attitude from two vector measurements, has been widely used in both ground-based and onboard attitude determination. Later attitude determination methods have been based on Wahba's optimality criterion for n arbitrarily weighted observations. The solution of Wahba's problem is somewhat difficult in the general case, but there is a simple closed-form solution in the two-observation case. This solution reduces to the TRIAD solution for certain choices of measurement weights. This paper presents and compares these algorithms as well as sub-optimal algorithms proposed by Bar-Itzhack, Harman, and Reynolds. Some new results will be presented, but the paper is primarily a review and tutorial.

Description: The subject of neotectonics, covering the structures and structural activity of the last 5 million years (i.e., post-Miocene) is a well-recognized field, including "active tectonics," focussed on the last 500,000 years in a 1986 National Research Council report of that title. However, there is a cartographic gap between tectonic maps, generally showing all features regardless of age, and maps of current seismic or volcanic activity. We have compiled a map intended to bridge this gap, using modern data bases and computer-aided cartographic techniques. The maps presented here are conceptually descended from an earlier map showing tectonic and volcanic activity of the last one million years. Drawn by hand with the National Geographic Society's 1975 "The Physical World" map as a base, the 1981 map in various revisions has been widely reproduced in textbooks and various technical publications. However, two decades of progress call for a completely new map that can take advantage of new knowledge and cartographic techniques. The digital tectonic activity map (DTM), presented in shaded relief (Fig. 1) and schematic (Fig. 2) versions, is the result. The DTM is intended to show tectonism and volcanism of the last one million years, a period long enough to be representative of global activity, but short enough that features such as fault scarps and volcanos are still geomorphically recognizable. Data Sources and Cartographic Methods The DTM is based on a wide range of sources, summarized in Table 1. The most important is the digital elevation model, used to construct a shaded relief map. The bathymetry is largely from satellite altimetry, specifically the marine gravity compilations by Smith and Sandwell (1996). The shaded relief map was designed to match the new National Geographic Society world physical map (1992), although drawn independently, from the digital elevation model. The Robinson Projection is used instead of the earlier Van der Grinten one. Although neither conformal nor equal-area, the Robinson Projection provides a reasonable compromise and retains useful detail at high latitudes.

Description: Vegetational changes are primary indicators of the present and future ecological status of the globe. These are changes which not only impact upon the primary productivity, but the total of the biogeochemical processes occurring on the planet. The impacts of global climatic and other environmental changes on vegetation must be monitored by some means in order to develop models which will allow us to predict long term effects. Large scale monitoring is now possible only with remote sensing systems, primarily passive reflectance, obtained by the use of satellite and aircraft platforms. However, passive reflectance techniques at this time are limited in their ability to detect subtle changes in the concentration and oxidation states of the many compounds involved in the light reactions of photosynthesis. Knowledge of these changes we consider to be fundamental in the remote assessment of both the rate and efficiency of photosynthesis and also the early detection of stress damage. The above factors pointed to the desirability of a sensing technique with the sensitivity and specificity necessary for detecting and quantifying those biological entities involved in photosynthesis. Another optical technique for vegetation monitoring is fluorescence. Previously, the lack of adequate excitation light sources and detector technologies have limited the use of fluorescence on intact plant leaves in the field. It is only recently with the advent of lasers with short pulse duration and advanced detector technologies that fluorescence measurements in the remote mode have become possible in the presence of ambient light.

Description: The age of secondary forests in the Amazon will become more critical with respect to the estimation of biomass and carbon budgets as tropical forest conversion continues. Multitemporal Thematic Mapper data were used to develop land cover histories for a 33,000 Square kM area near Ariquemes, Rondonia over a 7 year period from 1989-1995. The age of the secondary forest, a surrogate for the amount of biomass (or carbon) stored above-ground, was found to be unimportant in terms of biomass budget error rates in a forested TM scene which had undergone a 20% conversion to nonforest/agricultural cover types. In such a situation, the 80% of the scene still covered by primary forest accounted for over 98% of the scene biomass. The difference between secondary forest biomass estimates developed with and without age information were inconsequential relative to the estimate of biomass for the entire scene. However, in futuristic scenarios where all of the primary forest has been converted to agriculture and secondary forest (55% and 42% respectively), the ability to age secondary forest becomes critical. Depending on biomass accumulation rate assumptions, scene biomass budget errors on the order of -10% to +30% are likely if the age of the secondary forests are not taken into account. Single-date TM imagery cannot be used to accurately age secondary forests into single-year classes. A neural network utilizing TM band 2 and three TM spectral-texture measures (bands 3 and 5) predicted secondary forest age over a range of 0-7 years with an RMSE of 1.59 years and an R(Squared) (sub actual vs predicted) = 0.37. A proposal is made, based on a literature review, to use satellite imagery to identify general secondary forest age groups which, within group, exhibit relatively constant biomass accumulation rates.

Description: The use of TOMS (Total Ozone Mapping Spectrometer) satellite data demonstrates the recently developed technique of using satellite UV radiance measurements to detect absorbing tropospheric aerosols is effective over snow/ice surfaces. Instead of the traditional single wavelength (visible or infrared) method of measuring tropospheric aerosols, this method takes advantage of the wavelength dependent reduction in the backscattered radiance due to the presence of absorbing aerosols over snow/ice surfaces. An example of the resulting aerosol distribution derived from TOMS data is shown for an August 1998 event in which smoke generated by Canadian forest fires drifts over and across Greenland. As the smoke plume moved over Greenland, the TOMS observed 380 nm reflectivity over the snow/ice surface dropped drastically from 90-100% down to 30-40%. To study the effects of this smoke plume in both the UV and visible regions of the spectrum, we compared a smoke-laden spectrum taken over Greenland by the high spectral resolution (300 to 800 nm) GOME instrument with one that is aerosol-free. We also discuss the results of modeling the darkening effects of various types of absorbing aerosols over snow/ice surfaces using a radiative transfer code. Finally, we investigated the history of such events by looking at the nearly twenty year record of TOMS aerosol index measurements and found that there is a large interannual variability in the amount of smoke aerosols observed over Greenland. This information will be available for studies of radiation and transport properties in the Arctic.

Description: As a result of photochemistry, some relationship between the stratospheric age-of-air and the amount of tracer contained within an air sample is expected. The existence of such a relationship allows inferences about transport history to be made from observations of chemical tracers. This paper lays down the conceptual foundations for the relationship between age and tracer amount, developed within a Lagrangian framework. In general, the photochemical loss depends not only on the age of the parcel but also on its path. We show that under the "average path approximation" that the path variations are less important than parcel age. The average path approximation then allows us to develop a formal relationship between the age spectrum and the tracer spectrum. Using the relation between the tracer and age spectra, tracer-tracer correlations can be interpreted as resulting from mixing which connects parts of the single path photochemistry curve, which is formed purely from the action of photochemistry on an irreducible parcel. This geometric interpretation of mixing gives rise to constraints on trace gas correlations, and explains why some observations are do not fall on rapid mixing curves. This effect is seen in the ATMOS observations.

Description: The Global Modeling Initiative (GMI) science team is developing a three dimensional chemistry and transport model (CTM) to be used in assessment of the atmospheric effects of aviation. Requirements are that this model be documented, be validated against observations, use a realistic atmospheric circulation, and contain numerical transport and photochemical modules representing atmospheric processes. The model must also retain computational efficiency to be tractable to use for multiple scenarios and sensitivity studies. To meet these requirements, a facility model concept was developed in which the different components of the CTM are evaluated separately. The first use of the GMI model will be to evaluate the impact of the exhaust of supersonic aircraft on the stratosphere. The assessment calculations will depend strongly on the wind and temperature fields used by the CTM. Three meteorological data sets for the stratosphere are available to GMI: the National Center for Atmospheric Research Community Climate Model (CCM2), the Goddard Earth Observing System Data Assimilation System (GEOS DAS), and the Goddard Institute for Space Studies general circulation model (GISS). Objective criteria were established by the GMI team to identify the data set which provides the best representation of the stratosphere. Simulations of gases with simple chemical control were chosen to test various aspects of model transport. The three meteorological data sets were evaluated and graded based on their ability to simulate these aspects of stratospheric measurements. This paper describes the criteria used in grading the meteorological fields. The meteorological data set which has the highest score and therefore was selected for GMI is CCM2. This type of objective model evaluation establishes a physical basis for interpretation of differences between models and observations. Further, the method provides a quantitative basis for defining model errors, for discriminating between different models, and for ready re-evaluation of improved models. These in turn will lead to a higher level of confidence in assessment calculations.

Description: Supported by numerical experiment results, the abrupt change of the location of the intertropical convergence zone (ITCZ), from the equatorial trough flow regime to the monsoon trough flow regime is interpreted as a subcritical instability. The existence of these multiple quasi-equilibria is due to the balance of two "forces" on the ITCZ: one toward the equator, due to the earth's rotation, has a nonlinear latitudinal dependence; and the other toward the latitude of the sea surface (or ground) temperature peak has a relatively linear latitudinal dependence. This work pivots on the finding that the ITCZ and Hadley circulation can still exist without the pole-to-equator gradient of radiative-convective equilibrium temperature.

Description: Solar radiation is the primary source of energy driving atmospheric and oceanic circulations. Concerned with the huge computing time required for computing radiative transfer in weather and climate models, solar heating in minor absorption bands has often been neglected. The individual contributions of these minor bands to the atmospheric heating is small, but collectively they are not negligible. The solar heating in minor bands includes the absorption due to water vapor in the photosynthetically active radiation (PAR) spectral region from 14284/cm to 25000/cm, the ozone absorption and Rayleigh scattering in the near infrared, as well as the O2 and CO2 absorption in a number of weak bands. Detailed high spectral- and angular-resolution calculations show that the total effect of these minor absorption is to enhance the atmospheric solar heating by approximately 10%. Depending upon the strength of the absorption and the overlapping among gaseous absorption, different approaches are applied to parameterize these minor absorption. The parameterizations are accurate and require little extra time for computing radiative fluxes. They have been efficiently implemented in the various atmospheric models at NASA/Goddard Space Flight Center, including cloud ensemble, mesoscale, and climate models.

Description: There are procedures and methods for verification of coding algebra and for validations of models and calculations that are in use in the aerospace computational fluid dynamics (CFD) community. These methods would be efficacious if used by the glacier dynamics modelling community. This paper is a presentation of some of those methods, and how they might be applied to uncertainty management supporting code verification and model validation for glacier dynamics. The similarities and differences between their use in CFD analysis and the proposed application of these methods to glacier modelling are discussed. After establishing sources of uncertainty and methods for code verification, the paper looks at a representative sampling of verification and validation efforts that are underway in the glacier modelling community, and establishes a context for these within overall solution quality assessment. Finally, an information architecture and interactive interface is introduced and advocated. This Integrated Cryospheric Exploration (ICE) Environment is proposed for exploring and managing sources of uncertainty in glacier modelling codes and methods, and for supporting scientific numerical exploration and verification. The details and functionality of this Environment are described based on modifications of a system already developed for CFD modelling and analysis.

Description: We obtain valuable information on the angular and seasonal variability of surface reflectance using a hand-held spectrometer from a light aircraft. The data is used to test a procedure that allows us to estimate visible surface reflectance from the longer wavelength 2.1 micrometer channel (mid-IR). Estimating or avoiding surface reflectance in the visible is a vital first step in most algorithms that retrieve aerosol optical thickness over land targets. The data indicate that specular reflection found when viewing targets from the forward direction can severely corrupt the relationships between the visible and 2.1 micrometer reflectance that were derived from nadir data. There is a month by month variation in the ratios between the visible and the mid-IR, weakly correlated to the Normalized Difference Vegetation Index (NDVI). If specular reflection is not avoided, the errors resulting from estimating surface reflectance from the mid-IR exceed the acceptable limit of DELTA-rho approximately 0.01 in roughly 40% of the cases, using the current algorithm. This is reduced to 25% of the cases if specular reflection is avoided. An alternative method that uses path radiance rather than explicitly estimating visible surface reflectance results in similar errors. The two methods have different strengths and weaknesses that require further study.

Description: Space Flight hardware and software designers are increasingly turning to Commercial-Off-the-Shelf (COTS) products in hopes of meeting the demands imposed on them by projects with short development cycle times. The Technology Validation Assurance (TVA) team at NASA GSFC has embarked on applying a method for inserting COTS hardware into the Spartan 251 spacecraft. This method includes Procurement, Characterization, Ruggedization/Remediation and Verification Testing process steps which are intended to increase the user's confidence in the hardware's ability to function in the intended application for the required duration. As this method is refined with use, it has the potential for becoming a benchmark for industry-wide use of COTS in high reliability systems.

Description: Volcanic eruptions loft gases and ash particles into the atmosphere and produce effects that are both short term (aircraft hazards, interference with satellite measurements) and long term (atmospheric chemistry, climate). Large (greater than 0.5mm) ash particles fall out in minutes [Rose et al, 1995], but fine ash particles can remain in the atmosphere for many days. This fine volcanic ash is a hazard to modem jet aircraft because the operating temperatures of jet engines are above the solidus temperature of volcanic ash, and because ash causes abrasion of windows and airframe, and disruption of avionics. At large distances(10(exp 2)-10(exp 4) km or more) from their source, drifting ash clouds are increasingly difficult to distinguish from meteorological clouds, both visually and on radar [Rose et al., 1995]. Satellites above the atmosphere are unique platforms for viewing volcanic clouds on a global basis and measuring their constituents and total mass. Until recently, only polar AVHRR and geostationary GOES instruments could be used to determine characteristics of drifting volcanic ash clouds using the 10-12 micron window [Prata 1989; Wen and Rose 1994; Rose and Schneider 1996]. The NASA Total Ozone Mapping Spectrometer (TOMS) instruments aboard the Nimbus-7, Meteor3, ADEOS, and Earth Probe satellites have produced a unique data set of global SO2 volcanic emissions since 1978 (Krueger et al., 1995). Besides SO2, a new technique has been developed which uses the measured spectral contrast of the backscattered radiances in the 330-380nm spectral region (where gaseous absorption is negligible) in conjunction with radiative transfer models to retrieve properties of volcanic ash (Krotkov et al., 1997) and other types of absorbing aerosols (Torres et al., 1998).

Description: The Support Vector Machine provides a new way to design classification algorithms which learn from examples (supervised learning) and generalize when applied to new data. We demonstrate its success on a difficult classification problem from hyperspectral remote sensing, where we obtain performances of 96%, and 87% correct for a 4 class problem, and a 16 class problem respectively. These results are somewhat better than other recent results on the same data. A key feature of this classifier is its ability to use high-dimensional data without the usual recourse to a feature selection step to reduce the dimensionality of the data. For this application, this is important, as hyperspectral data consists of several hundred contiguous spectral channels for each exemplar. We provide an introduction to this new approach, and demonstrate its application to classification of an agriculture scene.

Description: The behavior of mesospheric odd nitrogen species during and following relativistic and diffuse auroral precipitation events is simulated, Below 75 km nitric oxide is enhanced in proportion to the ion pair production function associated with the electron precipitation and the length of the event. Nitrogen dioxide and nitric acid are also enhanced. At 65 km the percentage of odd nitrogen for N is 0.1%, HNO3 is 1.6%, NO2 is 15%, and NO is 83.3%. Between 75 and 85 km NO is depleted during particle events due to the faster destruction of NO by N relative to the production of NO by N reacting with O2. Recovery of NO depends on transport from the lower thermosphere, where NO is produced in abundant amounts during particle events.

Description: Long-lived tropospheric source gases, such as nitrous oxide, enter the stratosphere through the tropical tropopause, are transported throughout the stratosphere by the Brewer-Dobson circulation, and are photochemically destroyed in the upper stratosphere. These chemical constituents, or "tracers" can be used to track mixing and transport by the stratospheric winds. Much of our understanding about the stratospheric circulation is based on large scale gradients and other spatial features in tracer fields constructed from satellite measurements. The point of view presented in this paper is different, but complementary, in that transport is described in terms of tracer probability distribution functions (PDFs). The PDF is computed from the measurements, and is proportional to the area occupied by tracer values in a given range. The flavor of this paper is tutorial, and the ideas are illustrated with several examples of transport-related phenomena, annotated with remarks that summarize the main point or suggest new directions. One example shows how the multimodal shape of the PDF gives information about the different branches of the circulation. Another example shows how the statistics of fluctuations from the most probable tracer value give insight into mixing between different regions of the atmosphere. Also included is an analysis of the time-dependence of the PDF during the onset and decline of the winter circulation, and a study of how "bursts" in the circulation are reflected in transient periods of rapid evolution of the PDF. The dependence of the statistics on location and time are also shown to be important for practical problems related to statistical robustness and satellite sampling. The examples illustrate how physically-based statistical analysis can shed some light on aspects of stratospheric transport that may not be obvious or quantifiable with other types of analyses. An important motivation for the work presented here is the need for synthesis of the large and growing database of observations of the atmosphere and the vast quantities of output generated by atmospheric models.

Description: Magnetic remanence of crustal rocks can reside in three common rock-forming magnetic minerals: magnetite, pyrrhotite, and hematite. Thermoremanent magnetization (TRM) of magnetite and pyrrhotite is carried mostly by single domain (SD) grains. The TRM of hematite grains, however, is carried mostly by multidomain (NM) grains. This characteristic is illustrated by TRM acquisition curves for hematite of variable grainsizes. The transition between truly NM behavior and tendency towards SD behavior his been established between hematite grainsizes of 0. 1 and 0.05 mm. Coarse grainsize of lower crustal rocks and the large sensitivity of MD hematite grains to acquire TRM indicates that hematite could be a significant contributor to long-wavelength magnetic anomalies.

Description: We present an extension for the 2D (zonal mean) version of our Numerical Spectral Mode (NSM) that incorporates Hines' Doppler spread parameterization (DSP) for small scale gravity waves (GW). This model is applied to describe the seasonal variations and the semi-annual and quasi-biennial oscillations (SAO and QBO). Our earlier model reproduced the salient features of the mean zonal circulation in the middle atmosphere, including the QBO extension into the upper mesosphere inferred from UARS measurements. In the present model we incorporate also tropospheric heating to reproduce the upwelling at equatorial latitudes associated with the Brewer-Dobson circulation that affects significantly the dynamics of the stratosphere as Dunkerton had pointed out. Upward vertical winds increase the period of the QBO observed from the ground. To compensate for that, one needs to increase the eddy diffusivity and the GW momentum flux, bringing the latter closer to values recommended in the DSP. The QBO period in the model is 30 months (mo), which is conducive to synchronize this oscillation with the seasonal cycle of solar forcing. Multi-year interannual oscillations are generated through wave filtering by the solar driven annual oscillation in the zonal circulation. Quadratic non-linearities generate interseasonal variations to produce a complicated pattern of variability associated with the QBO. The computed temperature amplitudes for the SAO and QBO are in substantial agreement with observations at equatorial and extratropical latitudes. At high latitudes, however, the observed QBO amplitudes are significantly larger, which may be a signature of propagating planetary waves not included in the present model. The assumption of hydrostatic equilibrium not being imposed, we find that the effects from the vertical Coriolis force associated with the equatorial oscillations are large for the vertical winds and significant for the temperature variations even outside the tropics but are relatively small for the zonal winds.

Description: Snow-cover maps generated from surface data are based on direct measurements, however they are prone to interpolation errors where climate stations are sparsely distributed. Snow cover is clearly discernable using satellite-attained optical data because of the high albedo of snow, yet the surface is often obscured by cloud cover. Passive microwave (PM) data is unaffected by clouds, however, the snow-cover signature is significantly affected by melting snow and the microwaves may be transparent to thin snow (less than 3cm). Both optical and microwave sensors have problems discerning snow beneath forest canopies. This paper describes a method that combines ground and satellite data to produce a Multiple-Dataset Snow-Cover Product (MDSCP). Comparisons with current snow-cover products show that the MDSCP draws together the advantages of each of its component products while minimizing their potential errors. Improved estimates of the snow-covered area are derived through the addition of two snow-cover classes ("thin or patchy" and "high elevation" snow cover) and from the analysis of the climate station data within each class. The compatibility of this method for use with Moderate Resolution Imaging Spectroradiometer (MODIS) data, which will be available in 2000, is also discussed. With the assimilation of these data, the resolution of the MDSCP would be improved both spatially and temporally and the analysis would become completely automated.