Search NASA⌕ Search

SEARCH · Search NASA

Results for “terrain classification”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 163 records · Page 9

Real-Time Application of Multi-Satellite Precipitation Analysis for Floods and Landslides

Satellite data acquired and processed in real time now have the potential to provide the spacetime information on rainfall needed to monitor flood and landslide events around the world. This can be achieved by integrating the satellite-derived forcing data with hydrological models and landslide algorithms. Progress in using the TRMM Multi-satellite Precipitation Analysis (TMPA) as input to flood and landslide forecasts is outlined, with a focus on understanding limitations of the rainfall data and impacts of those limitations on flood/landslide analyses. Case studies of both successes and failures will be shown, as well as comparison with ground comparison data sets-- both in terms of rainfall and in terms of flood/landslide events. In addition to potential uses in real-time, the nearly ten years of TMPA data allow retrospective running of the models to examine variations in extreme events. The flood determination algorithm consists of four major components: 1) multi-satellite precipitation estimation; 2) characterization of land surface including digital elevation from NASA SRTM (Shuttle Radar Terrain Mission), topography-derived hydrologic parameters such as flow direction, flow accumulation, basin, and river network etc.; 3) a hydrological model to infiltrate rainfall and route overland runoff; and 4) an implementation interface to relay the input data to the models and display the flood inundation results to potential users and decision-makers, In terms of landslides, the satellite rainfall information is combined with a global landslide susceptibility map, derived from a combination of global surface characteristics (digital elevation topography, slope, soil types, soil texture, and land cover classification etc.) using a weighted linear combination approach. In those areas identified as "susceptible" (based on the surface characteristics), landslides are forecast where and when a rainfall intensity/duration threshold is exceeded. Results are described indicating general agreement with landslide occurrences.

Adler, Robert↗

Application of Multi-Satellite Precipitation Analysis to Floods and Landslides

Satellite data acquired and processed in real time now have the potential to provide the spacetime information on rainfall needed to monitor flood and landslide events around the world. This can be achieved by integrating the satellite-derived forcing data with hydrological models and landslide algorithms. Progress in using the TRMM Multi-satellite Precipitation Analysis (TMPA) as input to flood and landslide forecasts is outlined, with a focus on understanding limitations of the rainfall data and impacts of those limitations on flood/landslide analyses. Case studies of both successes and failures will be shown, as well as comparison with ground comparison data sets both in terms of rainfall and in terms of flood/landslide events. In addition to potential uses in real-time, the nearly ten years of TMPA data allow retrospective running of the models to examine variations in extreme events. The flood determination algorithm consists of four major components: 1) multi-satellite precipitation estimation; 2) characterization of land surface including digital elevation from NASA SRTM (Shuttle Radar Terrain Mission), topography-derived hydrologic parameters such as flow direction, flow accumulation, basin, and river network etc.; 3) a hydrological model to infiltrate rainfall and route overland runoff; and 4) an implementation interface to relay the input data to the models and display the flood inundation results to potential users and decision-makers. In terms of landslides, the satellite rainfall information is combined with a global landslide susceptibility map, derived from a combination of global surface characteristics (digital elevation topography, slope, soil types, soil texture, and land cover classification etc.) using a weighted linear combination approach. In those areas identified as "susceptible" (based on the surface characteristics), landslides are forecast where and when a rainfall intensity/duration threshold is exceeded. Results are described indicating general agreement with landslide occurrences. However, difficulties in comparing landslide event information (mostly from news reports) with the satellite-based forecasts are analyzed.

Adler, Robert↗

Spectral signature selection for mapping unvegetated soils

Airborne multispectral scanner data covering the wavelength interval from 0.40-2.60 microns were collected at an altitude of 1000 m above the terrain in southeastern Pennsylvania. Uniform training areas were selected within three sites from this flightline. Soil samples were collected from each site and a procedure developed to allow assignment of scan line and element number from the multispectral scanner data to each sampling location. These soil samples were analyzed on a spectrophotometer and laboratory spectral signatures were derived. After correcting for solar radiation and atmospheric attenuation, the laboratory signatures were compared to the spectral signatures derived from these same soils using multispectral scanner data. Both signatures were used in supervised and unsupervised classification routines. Computer-generated maps using the laboratory and multispectral scanner derived signatures resulted in maps that were similar to maps resulting from field surveys. Approximately 90% agreement was obtained between classification maps produced using multispectral scanner derived signatures and laboratory derived signatures.

May, G. A.↗

Characterizing Detailed Grain Shape and Size Distribution Properties of Lunar Regolith

Introduction: As the nation prepares to return to the Moon, there is an increasing need for testing tools, instruments, and equipment in simulated environments on Earth to ensure successful operations during lunar missions. Regolith will affect all aspects of future lunar missions, from plume interactions during landing to space suit and tool design [1]. Because of this, it is important to understand the grain shape and size properties of lunar regolith and how those influence regolith behavior in order to prepare for these missions. This knowledge is also vital to create more accurate lunar regolith simulants for testing equipment in a lunar environment. While particle size analyses have been performed on most Apollo soils using simple sieving, shape has only been crudely addressed [2]. This work analyzes 4 lunar regolith samples to provide a better understanding of these size and shape parameters and will provide a more accurate baseline of data to create high-fidelity lunar regolith simulants. Methods: New technologies exist today that are capable of measuring size and shape simultaneously for hundreds of thousands of particles in a single measurement. We conducted a rigorous analysis of the particle size distribution (PSD) as well as the size-dependent 2D and 3D shape parameters of lunar regolith samples of different compositions and maturity levels. This analysis was done using a Microtrac SYNC which provides a unique combination of tri-laser diffraction and Dynamic Image Analysis (DIA). Sample Selection. Four regolith samples were selected for analysis based on maturity level and lunar terrain type: • 10084 – Mature high-Ti mare regolith • 15601 – Immature low-Ti mare regolith • 64501 – Mature highland regolith • 67461 – Immature highland regolith Sample Analysis. After receiving the samples, each sample was imaged with an optical microscope (Figure 1). To obtain 2D and 3D particle size and shape, 0.1 g of each sample was then added to the SYNC for analysis, which outputs more than 30 size and shape parameters for each individual grain as well as the complete size distribution from 0.01-2000 μm by blending laser diffraction and DIA together. For the 0.1 g sample masses requested, ~100,000 grains per sample were captured by DIA. Results: From the PSD analysis, it can be seen that the average particle size of 10084 is ~24.5 µm which is much smaller than the average particle size of the Apollo sample collection (~72 µm) [3], while the average particle size of samples 15601, 64501, and 67461 are larger than the Apollo sample average at ~106 µm, ~103 µm, and ~118.5 µm respectively. The size and shape measurements for the samples output ~30 parameters, four of which were focused on for this study: sphericity, aspect ratio, roundness, and concavity (Figure 3). However, after analyzing the plots it was found that only sphericity and aspect ratio showed differences between the samples. These two parameters are measured on a scale from 0 to 1, with 1 being a perfect sphere with equal dimensions. The results show that the sphericity values are slightly high-er for the mature samples 10084 and 64501 (~0.96) than they are for samples 15601 and 67461 (~0.95) (Figure 4). The aspect ratio values are slightly lower for samples 64501 (~0.7) and 67461 (~0.74) than they are for samples 10084 (~0.76) and 15601 (~0.8) (Figure 5). Discussion: Mature regoliths are those that have been exposed to micrometeorite and solar wind bombardment for long periods of time, breaking up particles and causing them to become more rounded [4]. Therefore, the smaller particle sizes and higher sphericity values for samples 10084 and 64501 are expected due to their higher maturity level compared to 15601 and 67461. However, the aspect ratio values are not dependent on maturity and are instead dependent on terrain type. The lower aspect ratio values for the high-land samples is potentially due to the higher plagioclase content, which occurs in elongated particles and does not break down as easily as the pyroxenes or olivines that are present in the mare. Conclusion and Future Work: The results of this work provide a baseline of high-quality data that will contribute to the creation of high-fidelity lunar simulants and will greatly benefit NASA’s efforts of establishing a human presence on the Moon. Future work includes performing unsupervised image classification on the ~105 particle images per sample in order to identify different classes of grains. These grain classes can then be linked to detailed shape properties, and the relative abundance of each class in the samples can be compared. Acknowledgments: We would like to thank the Extraterrestrial Materials Analysis Group (ExMAG) and the Astromaterials Allocation Review Board (AARB) for allocating lunar regolith samples 10084,9010, 15601,365, 64501,249, and 67461,171 for this work. References: [1] Taylor, L. A. et al. (2005) AIAA #2510. [2] Katagiri et al. (2015) ASCE. [3] Carrier III, W. D. (2005) Lunar Geotech Institute Tech Report. [4] McKay, D. S. et al. (1991) The Lunar Sourcebook, Chapter 7.

S R Deitrick↗

Measurements from an Aerial Vehicle: A New Tool for Planetary Exploration

Aerial vehicles fill a unique planetary science measurement gap, that of regional-scale, near-surface observation, while providing a fresh perspective for potential discovery. Aerial vehicles used in planetary exploration bridge the scale and resolution measurement gaps between orbiters (global perspective with limited spatial resolution) and landers (local perspective with high spatial resolution) thus complementing and extending orbital and landed measurements. Planetary aerial vehicles can also survey scientifically interesting terrain that is inaccessible or hazardous to landed missions. The use of aerial assets for performing observations on Mars, Titan, or Venus will enable direct measurements and direct follow-ons to recent discoveries. Aerial vehicles can be used for remote sensing of the interior, surface and atmosphere of Mars, Venus and Titan. Types of aerial vehicles considered are airplane "heavier than air" and airships and balloons "lighter than air". Interdependencies between the science measurements, science goals and objectives, and platform implementation illustrate how the proper balance of science, engineering, and cost, can be achieved to allow for a successful mission. Classification of measurement types along with how those measurements resolve science questions and how these instruments are accommodated within the mission context are discussed.

Wright, Henry S.↗

An interdisciplinary analysis of Colorado Rocky Mountain environments using ADP techniques

The author has identified the following significant results. Good ecological, classification accuracy (90-95%) can be achieved in areas of rugged relief on a regional basis for Level 1 cover types (coniferous forest, deciduous forest, grassland, cropland, bare rock and soil, and water) using computer-aided analysis techniques on ERTS/MSS data. Cost comparisons showed that a Level 1 cover type map and a table of areal estimates could be obtained for the 443,000 hectare San Juan Mt. test site for less than 0.1 cent per acre, whereas photointerpretation techniques would cost more than 0.4 cent per acre. Results of snow cover mapping have conclusively proven that the areal extent of snow in mountainous terrain can be rapidly and economically mapped by using ERTS/MSS data and computer-aided analysis techniques. A distinct relationship between elevation and time of freeze or thaw was observed, during mountain lake mapping. Basic lithologic units such as igneous, sedimentary, and unconsolidated rock materials were successfully identified. Geomorphic form, which is exhibited through spatial and textual data, can only be inferred from ERTS data. Data collection platform systems can be utilized to produce satisfactory data from extremely inaccessible locations that encounter very adverse weather conditions, as indicated by results obtained from a DCP located at 3,536 meters elevation that encountered minimum temperatures of -25.5 C and wind speeds of up to 40.9m/sec (91 mph), but which still performed very reliably.

Hoffer, R. M.↗

ICESat-2 Photon Classification: Finding Signal Photons in the ATL03 Geolocated Photon Data Product

ICESat-2 carries NASA's next-generation laser altimeter, ATLAS, (Advanced Topographic Laser Altimeter System), designed to measure changes in ice sheet height, sea ice freeboard, and vegetation canopy height. ATLAS contains a photon-counting lidar which transmits green (532-nm) pulses at 10kHz. Each pulse is split into 3 pairs of beams (one strong and one weak). Approximately 1014 photons per pulse travel from ATLAS through the atmosphere to reflect off the Earth's surface. Some return back into the ATLAS telescope where they are recorded. Photons from sunlight and instrument noise at the same wavelength are also recorded. The flight software time tags all photons within a 500m to 6 km range window and generates histograms. Using the histograms, it selects a telemetry window which varies from 20m over flat surfaces to hundreds of meters over rougher terrain. ATL03 contains the time, height (relative to the WGS-84 ellipsoid), latitude and longitude of every photon within the telemetry window. The basic challenge is to determine which of these photons were reflected off the surface. We have developed an algorithm that identifies these signal photons and assigns a confidence level (low, medium, or high) to each signal photon based on the signal to noise ratio. We present an overview of the signal identification algorithm and show the results on actual ICESat-2 data over ice sheet, sea ice, vegetated, and water surfaces. Higher level ATLAS products work with aggregations of the photons in order to determine the ellipsoidal height of the Earth, canopy height and structure, and other quantities of geophysical interest.

Remote Sensing↗

Evaluation of the radar response to land surfaces and volumes: Examination of theoretical models, target statistics, and applications

Four areas of L-band radar remote sensing of terrain were examined: (1) the behavior of the radar backscatter coefficient of distributed surface and volumes as a function of the targets' dielectric and geometric parameters and as a fnction of their physical parameters; (2) the correspondence of the angular behavior of the relative backscatter coefficient as extracted from SIR-B digital imagery and truck mounted L-band scatterometer measurements for about 100 fields; (3) the statistical behavior of SIR-B image density for targets that appear homogeneous on Thematic Mapper (TM) optical imagery and/or color IR photography; and (4) the applicability of SIR-B imagery both alone and in conjunction with TM imagery for the classification and monitoring of land cover and renewable resources.

Ulaby, F. T.↗

Computer-aided analysis of Skylab multispectral scanner data in mountainous terrain for land use, forestry, water resource, and geologic applications

The author has identified the following significant results. One of the most significant results of this Skylab research involved the geometric correction and overlay of the Skylab multispectral scanner data with the LANDSAT multispectral scanner data, and also with a set of topographic data, including elevation, slope, and aspect. The Skylab S192 multispectral scanner data had distinct differences in noise level of the data in the various wavelength bands. Results of the temporal evaluation of the SL-2 and SL-3 photography were found to be particularly important for proper interpretation of the computer-aided analysis of the SL-2 and SL-3 multispectral scanner data. There was a quality problem involving the ringing effect introduced by digital filtering. The modified clustering technique was found valuable when working with multispectral scanner data involving many wavelength bands and covering large geographic areas. Analysis of the SL-2 scanner data involved classification of major cover types and also forest cover types. Comparison of the results obtained wth Skylab MSS data and LANDSAT MSS data indicated that the improved spectral resolution of the Skylab scanner system enabled a higher classification accuracy to be obtained for forest cover types, although the classification performance for major cover types was not significantly different.

Hoffer, R. M.↗

Lineament and polygon patterns on Europa

A classification scheme is presented for the lineaments and associated polygonal patterns observed on the surface of Europa, and the frequency distribution of the polygons is discussed in terms of the stress-relief fracturing of the surface. The lineaments are divided on the basis of albedo, morphology, orientation and characteristic geometry into eight groups based on Voyager 2 images taken at a best resolution of 4 km. The lineaments in turn define a system of polygons varying in size from small reticulate patterns the limit of resolution to 1,000,000 sq km individuals. Preliminary analysis of polygon side frequency distributions reveals a class of polygons with statistics similar to those found in complex terrestrial terrains, particularly in areas of well-oriented stresses, a class with similar statistics around the antijovian point, and a class with a distribution similar to those seen in terrestrial tensional fracture patterns. Speculations concerning the processes giving rise to the lineament patterns are presented.

Pieri, D. C.↗

Building a Bilingual Google Earth Engine Dashboard to Increase Accessibility to Long-term Time Series Remote Sensing Data for Monitoring Saline System Changes in Chile’s Atacama Desert

Saline systems, consisting of salt flats, ponds, and marshes, provide vital water resources to wildlife and communities in northern Chile’s Atacama Desert, one of the driest regions in the world. Mining is extensive in the Atacama, which contains 30% of the world’s lithium reserves and is abundant in potassium and boron. The groundwater that feeds into salt marshes and ponds is extracted in large volumes for mining operations, limiting the availability of water for ecosystems. However, identifying long-term and large-scale environmental impacts from local lithium mining on the saline systems is limited by region inaccessibility and terrain variability. Open access satellite imagery and cloud computing technology has made studying Atacama saline systems feasible and allowed for collaboration across different agencies and countries. The NASA DEVELOP Program partnered with Chile’s la Universidad de La Serena and Servicio Nacional de Geología y Minería (SERNAGEOMIN) to create the Saline Analysis Tool (SalT) in Google Earth Engine (GEE). SalT is used to analyze the extent and distribution of remote saline systems in the Atacama from 1986 to the present day. The tool filters Landsat 5 Thematic Mapper (TM) and Landsat 8 Operational Land Imager (OLI) data from GEE’s data catalog and creates a single composite image per year for analysis. Additional output analyses include land cover classification, Normalized Difference Vegetation Index (NDVI) and Normalized Difference Water Index (NDWI) raster images that can be displayed on the map interface or exported. The tool can also generate time-lapse videos and charts displaying NDVI, NDWI, and land cover over time. A key feature of the tool is the use of a bilingual graphical user interface to make analysis accessible and customizable to different users’ needs—SalT provides options to select an analysis area, analysis time period, and outputs to display or export. The tool also incorporates new Earth observations as they are added to GEE’s catalog. The ability to easily visualize and analyze long-term remote sensing imagery will enable SERNAGEOMIN and la Universidad de la Serena to continually monitor changes in these saline systems and inform future land management policy.

NASA DEVELOP↗

Interdisciplinary applications and interpretations of EREP data within the Susquehanna River Basin

The author has identified the following significant results. Photography from the S190A and S190B sensors was compared, for terrain analysis. The S190B photographs were used for terrain mapping of three test areas selected as being representative of major physiographic regions in Pennsylvania. Skylab photography was superior to both LANDSAT imagery and high altitude aircraft photography for purposes of accurate location of lineaments. Analysis of Skylab imagery has shown that long lineaments originally plotted on LANDSAT images are actually made up of shorter segments. Correlation of lineaments with ore deposits was determined following the preparation of a Pennsylvania mineral deposit map. Digital wave number analysis (spatial filtering) was attempted to determine if it can be used to enhance certain subtle features, and in particular, to locate and verify lineaments. Various spectral bands and channels of the MSS digital data were evaluated for their value in the classification and thematic mapping.

Mcmurtry, G. J.↗

Age discrimination among eruptives of Menengai Caldera, Kenya, using vegetation parameters from satellite imagery

Results are presented of an investigation to determine the degree to which digitally processed Landsat TM imagery can be used to discriminate among vegetated lava flows of different ages in the Menengai Caldera, Kenya. A selective series of five images, consisting of a color-coded Landsat 5 classification and four color composites, are compared with geologic maps. The most recent of more than 70 postcaldera flows within the caldera are trachytes, which are variably covered by shrubs and subsidiary grasses. Soil development evolves as a function of time, and as such supports a changing plant community. Progressively older flows exhibit the increasing dominance of grasses over bushes. The Landsat images correlated well with geologic maps, but the two mapped age classes could be further subdivided on the basis of different vegetation communities. It is concluded that field maps can be modified, and in some cases corrected by use of such imagery, and that digitally enhanced Landsat imagery can be a useful aid to field mapping in similar terrains.

Blodget, Herbert W.↗

Processing LiDAR Data to Predict Natural Hazards

ELF-Base and ELF-Hazards (wherein 'ELF' signifies 'Extract LiDAR Features' and 'LiDAR' signifies 'light detection and ranging') are developmental software modules for processing remote-sensing LiDAR data to identify past natural hazards (principally, landslides) and predict future ones. ELF-Base processes raw LiDAR data, including LiDAR intensity data that are often ignored in other software, to create digital terrain models (DTMs) and digital feature models (DFMs) with sub-meter accuracy. ELF-Hazards fuses raw LiDAR data, data from multispectral and hyperspectral optical images, and DTMs and DFMs generated by ELF-Base to generate hazard risk maps. Advanced algorithms in these software modules include line-enhancement and edge-detection algorithms, surface-characterization algorithms, and algorithms that implement innovative data-fusion techniques. The line-extraction and edge-detection algorithms enable users to locate such features as faults and landslide headwall scarps. Also implemented in this software are improved methodologies for identification and mapping of past landslide events by use of (1) accurate, ELF-derived surface characterizations and (2) three LiDAR/optical-data-fusion techniques: post-classification data fusion, maximum-likelihood estimation modeling, and hierarchical within-class discrimination. This software is expected to enable faster, more accurate forecasting of natural hazards than has previously been possible.

Fairweather, Ian↗

The Value of Long-Term (40 years) Airborne Gamma Radiation SWE Record for Evaluating Three Observation-Based Gridded SWE Data Sets by Seasonal Snow and Land Cover Classifications

Observation‐based long‐term gridded snow water equivalent (SWE) products are important assets for hydrological and climate research. However, an evaluation of the currently available SWE products has been limited due to the lack of independent SWE data that extend over a large range of environmental conditions. In this study, three daily long‐term SWE products (Special Sensor Microwave Imager and Sounder [SSMI/S] SWE, GlobSnow‐2 SWE, and University of Arizona [UA] SWE) we reevaluated by seasonal snow cover and land cover classifications over the conterminous United States from1982 to 2017, using the historical airborne gamma radiation SWE observations (20,738 measurements).We found that there are similar patterns in SSMI/S and GlobSnow‐2 SWE when compared against the gamma SWE. However, GlobSnow‐2 SWE had better agreement with gamma SWE than SSMI/S SWE in some forested‐type classes and maritime and prairie snow classes. As compared to SSMI/S and GlobSnow‐2SWE, UA SWE has much better agreement with gamma SWE in all land cover types and snow classes. Tree cover and topographic heterogeneity affect the agreement between the gamma and gridded SWE and accuracy of gamma SWE itself with the largest differences typically occurring when the percent tree cover was 80% or higher, the terrain slope was steeper than 2.5°, and the elevation range exceeded 100 m. The results demonstrate the reliability of the UA SWE products and the benefits of the gamma radiation approach to measure SWE, especially in forested regions.

Eunsang Cho↗

Field Tests With Trident Drill in Bishop Tuff Help Prepare for Future Missions to Moon and Mars

We performed drilling in volcanic deposits near Bishop California using an engineering model of the Honeybee Robotics TRIDENT (The Regolith and Ice Drill for Exploration of New Terrains) drill [1] a rotary percussive 1-meter class drill that is carried on the PRIME1and VIPER (Volatiles Investigating Polar Exploration Rover)[2] missions that launch in 2024. A similar drilling system was planned for the proposed Icebreaker Discovery class mission to Mars [3] and the Mars Life Explorer mission recommended by the 2020 Decadal Survey of planetary science [4]. The objectives of the project were (1) to use data collected by the drill for operational purposes as a probe of subsurface material properties in formations that are analogous to those that may be encountered on planetary surfaces; (2) correlate subsurface structures with those deduced from Ground Penetrating Radar (GPR); and (3) inspect the boreholes after they were drilled to test PERISCOPE (Probe for Exploring Regolith and Ice by Subsurface Classification of Organics, polycyclic aromatic hydrocarbons (PAHs), and Elements), a newly developed downhole UV fluorescence spectrometer [5].

Carol R. Stoker↗

Development of an information data base for watershed monitoring

Landsat multispectral scanner data, Defense Mapping Agency digital terrain data, conventional maps, and ground data were integrated to create a comprehensive information data base (the Image Based Information System), to monitor the water quality of the Lake Tahoe Basin. Landsat imagery was used as the planimetric base to which all other data were registered. A georeference image plane, which provided an interface between all data planes for the Lake Tahoe Basin data base, was created from the drainage basin map. The data base was used to extract each drainage basin for separate display. The Defense Mapping Agency-created elevation image was processed with VICAR software to produce a component representing slope magnitude, which was cross-tabulated with the drainage basin georeference table. Future applications of the data base include the development of precipitation modeling, surface runoff models, and classification of drainage basin cover types.

Smith, A. Y.↗

Multitemporal Snow Cover Mapping in Mountainous Terrain for Landsat Climate Data Record Development

A multitemporal method to map snow cover in mountainous terrain is proposed to guide Landsat climate data record (CDR) development. The Landsat image archive including MSS, TM, and ETM+ imagery was used to construct a prototype Landsat snow cover CDR for the interior northwestern United States. Landsat snow cover CDRs are designed to capture snow-covered area (SCA) variability at discrete bi-monthly intervals that correspond to ground-based snow telemetry (SNOTEL) snow-water-equivalent (SWE) measurements. The June 1 bi-monthly interval was selected for initial CDR development, and was based on peak snowmelt timing for this mountainous region. Fifty-four Landsat images from 1975 to 2011 were preprocessed that included image registration, top-of-the-atmosphere (TOA) reflectance conversion, cloud and shadow masking, and topographic normalization. Snow covered pixels were retrieved using the normalized difference snow index (NDSI) and unsupervised classification, and pixels having greater (less) than 50% snow cover were classified presence (absence). A normalized SCA equation was derived to independently estimate SCA given missing image coverage and cloud-shadow contamination. Relative frequency maps of missing pixels were assembled to assess whether systematic biases were embedded within this Landsat CDR. Our results suggest that it is possible to confidently estimate historical bi-monthly SCA from partially cloudy Landsat images. This multitemporal method is intended to guide Landsat CDR development for freshwaterscarce regions of the western US to monitor climate-driven changes in mountain snowpack extent.

Landsat↗