Search NASA⌕ Search

SEARCH · Search NASA

Results for “CLOUD MASK”

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 181 records · Page 10

Operational Cloud Detection in GEOS Imagery

The bispectral spatial coherence approach for the detection and masking of clouds is revisited with modifications for operational applications to Geostationary Operational Environmental Satellite (GOES) Imager and Sounder data. The approach applies an edge detection algorithm to a "difference" image to identify the edges or contours of the clouds. The difference image is generated as the difference between the 3.9 and 11 micrometer channels of the GOES sensors. This difference image better highlights cloud edges than either image alone and draws upon the emissivity differences in the various channels. The cloud contours are "filled-in" by referencing the thermal gradients in the 11-micrometer channel or difference image associated with the cloud edge. This paper and associate poster will present a preliminary validation of the bispectral spatial coherence approach and will compare the cloud masks generated with this approach to the National Oceanic and Atmospheric Administration/National Environmental Satellite, Data, and Information Service (NOAA/NESDIS) operational cloud products.

Jedlovec, Gary J.↗

The Harmonized Landsat and Sentinel-2 Surface Reflectance Data Set

The Harmonized Landsat and Sentinel-2 (HLS) project is a NASA initiative aiming to produce a VirtualConstellation (VC) of surface reflectance (SR) data acquired by the Operational Land Imager (OLI) and MultiSpectral Instrument (MSI) aboard Landsat 8 and Sentinel-2 remote sensing satellites, respectively. The HLS products are based on a set of algorithms to obtain seamless products from both sensors (OLI and MSI): atmospheric correction, cloud and cloud-shadow masking, spatial co-registration and common gridding, bidirectional reflectance distribution function normalization and spectral bandpass adjustment. Three products are derivedfrom the HLS processing chain: (i) S10: full resolution MSI SR at 10 m, 20 m and 60 m spatial resolutions; (ii)S30: a 30 m MSI Nadir BRDF (Bidirectional Reflectance Distribution Function)-Adjusted Reflectance (NBAR);(iii) L30: a 30 m OLI NBAR. All three products are processed for every Level-1 input products from Landsat 8/OLI (L1T) and Sentinel-2/MSI (L1C). As of version 1.3, the HLS data set covers 10.35 million km2 and spans from first Landsat 8 data (2013); Sentinel-2 data spans from October 2015. The L30 and S30 show a good consistency with coarse spatial resolution products, in particular MODIS Collection 6 MCD09CMG products (overall deviations do not exceed 11%) that are used as a reference for quality assurance. The spatial co-registration of the HLS is improved compared to original Landsat 8 L1T and Sentinel 2A L1C products, for which misregistration issues between multi-temporal data are known. In particular, the resulting computed circular errors at 90% for the HLS product are 6.2 m and 18.8 m, for S10 and L30 products, respectively. The main known issue of the current data set remains the Sentinel-2 cloud mask with many cloud detection omissions. The cross-comparison with MODIS was used to flag products with most evident non-detected clouds. A time series outlier filtering approach is suggested to detect remaining clouds. Finally, several time series are presented to highlight the high potential of the HLS data set for crop monitoring.

Landsat Sentinel-2↗

Spatially Varying Spectrally Thresholds for MODIS Cloud Detection

The EOS science team has developed an elaborate global MODIS cloud detection procedure, and the resulting MODIS product (MOD35) is used in the retrieval process of several geophysical parameters to mask out clouds. While the global application of the cloud detection approach appears quite robust, the product has some shortcomings on the regional scale, often over determining clouds in a variety of settings, particularly at night. This over-determination of clouds can cause a reduction in the spatial coverage of MODIS derived clear-sky products. To minimize this problem, a new regional cloud detection method for use with MODIS data has been developed at NASA's Global Hydrology and Climate Center (GHCC). The approach is similar to that used by the GHCC for GOES data over the continental United States. Several spatially varying thresholds are applied to MODIS spectral data to produce a set of tests for detecting clouds. The thresholds are valid for each MODIS orbital pass, and are derived from 20-day composites of GOES channels with similar wavelengths to MODIS. This paper and accompanying poster will introduce the GHCC MODIS cloud mask, provide some examples, and present some preliminary validation.

Haines, S. L.↗

Overview and Statistical Analysis of Boundary Layer Clouds and Precipitation Over the Western North Atlantic Ocean

Due to their fast evolution and large natural variability in macro- and microphysical properties, the accurate representation of boundary layer clouds in current climate models remains a challenge. One of the regions with large intermodel spread in the Coupled Model Intercomparison Project Phase 6 ensemble is the western North Atlantic Ocean. Here, statistically representative in situ measurements can help to develop and constrain the parameterization of clouds in global models. To this end, we performed comprehensive measurements of boundary layer clouds, aerosol, trace gases, and radiation in the western North Atlantic Ocean during the NASA Aerosol Cloud meTeorology Interactions oVer the western ATlantic Experiment (ACTIVATE) mission. In total, 174 research flights with 574 flight hours for cloud and precipitation measurements were performed with the HU-25 Falcon during three winter (February–March 2020, January–April 2021, and November 2021–March 2022) and three summer seasons (August–September 2020, May–June 2021, and May–June 2022). Here we present a statistical evaluation of 16 140 individual cloud events probed by the fast cloud droplet probe and the two-dimensional stereo cloud probe during 155 research flights in a representative and repetitive flight strategy allowing for robust statistical data analyses. We show that the vertical profiles of distributions of the liquid water content and the cloud droplet effective diameter (ED) increase with altitude in the marine boundary layer. Due to higher updraft speeds, higher cloud droplet number concentrations (N liquid ) were measured in winter compared to summer despite lower cloud condensation nucleus abundance. Flight cloud cover derived from statistical analysis of in situ data is reduced in summer and shows large variability. This seasonal contrast in cloud coverage is consistent with a dominance of a synoptic pattern in winter that favors conditions for the formation of stratiform clouds at the western edge of cyclones (post-cyclonic). In contrast, a dominant summer anticyclone is concomitant with the occurrence of shallow cumulus clouds and lower cloud coverage. The evaluation of boundary layer clouds and precipitation in the N liquid ED phase space sheds light on liquid, mixed-phase, and ice cloud properties and helps to categorize the cloud data. Ice and liquid precipitation, often masked in cloud statistics by a high abundance of liquid clouds, is often observed throughout the cloud. The ACTIVATE in situ cloud measurements provide a wealth of cloud information useful for assessing airborne and satellite remote-sensing products, for global climate and weather model evaluations, and for dedicated process studies that address precipitation and aerosol–cloud interactions.

Aerosol↗

Changes in Tropical Clouds and Atmospheric Circulation Associated with Rapid Adjustment Induced by Increased Atmospheric CO2 – A Multiscale Modeling Framework Study

The radiative heating increase due to increased CO2 concentration is the primary source for the rapid adjustment of atmospheric circulation and clouds. In this study, we investigate the rapid adjustment resulting from doubling of CO2 and its physical mechanism using a multiscale modeling framework (MMF). The MMF includes an advanced higher-order turbulence closure in its cloud-resolving model component and simulates realistic shallow and deep cloud climatology and boundary layer turbulence. The rapid adjustment over the tropics is characterized by 1) reduced ascent and descent strengths over the ocean, 2) increased lower tropospheric stability (LTS) over the subsidence region, 3) shoaling of planetary boundary layers over the ocean, 4) increased deep convection over lands and shift of cloud coverage from the ocean to lands, and 5) reduced sensible (SH) and latent heat (LH) fluxes over the oceanic deep convective regions. Unlike conventional general circulation models and another MMF, a reduction in the global-mean shortwave cloud radiative cooling is not simulated, due to the increase in low clouds at lower altitudes over the ocean, resulting from reduced cloud-top entrainment due to strengthened inversion. Changes in regional circulation play a key role in cloud changes and shift of cloud coverage to lands. Weaker energy transport resulting from water vapor and cloud CO2 masking effects reduces the upward motion and convective clouds in the oceanic regions. The ocean-land transports are linked to the partitioning of surface SH and LH fluxes that increases humidity over lands and enhances deep convection over the tropical lands.

Xu, Kuan-Man↗

Cloud Motion in the GOCI COMS Ocean Colour Data

The Geostationary Ocean Colour Imager (GOCI) instrument, on Koreas Communications, Oceans, and Meteorological Satellite (COMS), can produce a spectral artefact arising from the motion of clouds the cloud is spatially shifted and the amount of shift varies by spectral band. The length of time it takes to acquire all eight GOCI bands for a given slot (portion of a scene) is sucient to require that cloud motion be taken into account to fully mask or correct the eects of clouds in all bands. Inter-band correlations can be used to measure the amount of cloud shift, which can then be used to adjust the cloud mask so that the union of all shifted masks can act as a mask for all bands. This approach reduces the amount of masking required versus a simple expansion of the mask in all directions away from clouds. Cloud motion can also aect regions with unidentied clouds thin or fractional clouds that evade the cloud identication process yielding degraded quality in retrieved ocean colour parameters. Areas with moving and unidentied clouds require more elaborate masking algo-rithms to remove these degraded retrievals. Correction for the eects of moving fractional clouds may also be possible. The cloud shift information can be used to determine cloud motion and thus wind at the cloud levels on sub-minute timescales. The benecial and negative eects of moving clouds should be con-sidered for any ocean colour instrument design and associated data processing plans.

Robinson, Wayne D.↗

Automated Detection of Clouds in Satellite Imagery

Many different approaches have been used to automatically detect clouds in satellite imagery. Most approaches are deterministic and provide a binary cloud - no cloud product used in a variety of applications. Some of these applications require the identification of cloudy pixels for cloud parameter retrieval, while others require only an ability to mask out clouds for the retrieval of surface or atmospheric parameters in the absence of clouds. A few approaches estimate a probability of the presence of a cloud at each point in an image. These probabilities allow a user to select cloud information based on the tolerance of the application to uncertainty in the estimate. Many automated cloud detection techniques develop sophisticated tests using a combination of visible and infrared channels to determine the presence of clouds in both day and night imagery. Visible channels are quite effective in detecting clouds during the day, as long as test thresholds properly account for variations in surface features and atmospheric scattering. Cloud detection at night is more challenging, since only courser resolution infrared measurements are available. A few schemes use just two infrared channels for day and night cloud detection. The most influential factor in the success of a particular technique is the determination of the thresholds for each cloud test. The techniques which perform the best usually have thresholds that are varied based on the geographic region, time of year, time of day and solar angle.

Jedlovec, Gary↗

How well is Martian ozone inferred with reflectance spectroscopy?

The Mariner 9 UV spectrometer scanned from 2100 to 3500 Angstroms in one of its two spectral channels every 3 seconds with a spectral resolution of 15 Angstroms and an effective field-of-view of approximately 300 km(exp 2). The only gaseous absorption in the 2000 to 3000 Angstrom region was assumed to come from the Hartley band system of ozone, and therefore the amount of ozone was inferred by fitting this absorption feature with laboratory data of ozone absorption. Mars O3 as inferred from these spectra shows strong seasonal and latitudinal variation, with column abundances ranging from 0.2 microns at equatorial latitudes to 60 microns over the northern winter polar latitudes. The detectability limit of the spectrometer was approximately 3 microns. I use a radiative transfer model based on the discrete ordinate method to calculate synthetic radiance spectra. When typical amounts of dust and cloud are present, significant underestimation of O3 occurs. A factor of 3 times as much O3 is needed to generate the same spectrum for cloudy, dusty atmospheres as for a clear atmosphere. If the scattering properties of Martian clouds and dust were well known, then their appearance would not be a problem, as a model would be capable of retrieving the O3 abundance. However, these properties are not well known, which raises doubts about the effectiveness of the current UV spectroscopy technique used to measure O3. Spatial and temporal variability in temperature and water vapor have been claimed to account for the scatter of the data points. However, water vapor is a small source of odd hydrogen in the winter polar atmosphere, and may not account for most of the variability. Masking by clouds and dust may also account for some of the observed O3 variability, because the nature and opacity of the clouds and dust in the polar hood change dramatically in latitude and even on a day-to-day basis. As the maximum O3 abundance resides near the surface, spacecraft must be able to observe through the entire cloud and dust abundance in order to actually see the total O3 column abundance. If reflectance spectroscopy is used, as on Mariner 9, then the cloud and the airborne dust must be traversed twice; first by the incoming solar flux down to the surface, and then once again upon reflection from the surface out to the spacecraft.

Lindner, Bernhard Lee↗

H31G-1596: DeepSAT's CloudCNN: A Deep Neural Network for Rapid Cloud Detection from Geostationary Satellites

Cloud and cloud shadow detection has important applications in weather and climate studies. It is even more crucial when we introduce geostationary satellites into the field of terrestrial remote sensing. With the challenges associated with data acquired in very high frequency (10-15 mins per scan), the ability to derive an accurate cloud shadow mask from geostationary satellite data is critical. The key to the success for most of the existing algorithms depends on spatially and temporally varying thresholds,which better capture local atmospheric and surface effects.However, the selection of proper threshold is difficult and may lead to erroneous results. In this work, we propose a deep neural network based approach called CloudCNN to classify cloudshadow from Himawari-8 AHI and GOES-16 ABI multispectral data. DeepSAT's CloudCNN consists of an encoderdecoder based architecture for binary-class pixel wise segmentation. We train CloudCNN on multi-GPU Nvidia Devbox cluster, and deploy the prediction pipeline on NASA Earth Exchange (NEX) Pleiades supercomputer. We achieved an overall accuracy of 93.29% on test samples. Since, the predictions take only a few seconds to segment a full multispectral GOES-16 or Himawari-8 Full Disk image, the developed framework can be used for real-time cloud detection, cyclone detection, or extreme weather event predictions.

GOES-1↗

In situ Mars ozone detector

We propose sending a balloon-borne UV photometer sensor package to measure atmospheric ozone on Mars, and this package could be a Discovery Program sensor candidate. Past measurements of ozone on Mars are highly uncertain, perhaps a factor of 3 or so uncertain, due primarily to interference and masking by cloud and dust. In-situ balloon measurements would avoid these problems, and would provide 'ground truth' which remote sensing techniques cannot. We have explored this approach to measure ozone abundance in the terrestrial stratosphere with a balloon-borne UV absorption photometer. Atmospheric pressures and temperatures and ozone concentrations near the surface of Mars are similar to those in the terrestrial stratosphere.

Lindner, Bernhard Lee↗

Arctic sea ice leads from advanced very high resolution radiometer images

A large number of advanced very high resolution radiometer (AVHRR) images from throughout 1989 are analyzed to determine lead characteristics. The units of analysis are square 200-km cells, and there are 270 such cells in the data set. Clouds are masked manually. Leads determine from images of the potential open water delta, a scaled version of the surface temperature or albedo that weights thin ice by its thermal or brightness impact. The lead fraction is determined as the mean delta, the monthly mean lead fraction ranges from 0.02 in winter to 0.06 in summer in the central Arctic and is near 0.08 in the winter in the peripheral seas. A method of accounting for lead width sampling errors due to the finite sample areas is introduced. In the central Arctic the observed mean lead width for a threshold of delta = 0.1 ranges from 2 or 3 km (near the resolution of the instrument) in the winter to 6 km in the summer. In the peripheral seas it is about 5 km in the winter. Width distributions are often more heavily weighted in the tail than exponential distributions and are well approximated by a power law. The along-track, number density power law N = aw(exp -6) has a mean exponent of b = 1.60 (standard deviation 0.18) and shows some seasonal variability. Mean floe widths in the central Arctic are 40 to 50 km in the winter, dropping to about 10 km in the summer. For floes the power law has a mean exponent of 0.93 and exhibits a clearer annual cycle. Lead orientation is determined with a method based on the direction of maximum extent.

Lindsay, R. W.↗

Harmonized Landsat/Sentinel-2 Products for Land Monitoring

The Harmonized Landsat-8 and Sentinel-2 (HLS) project is a NASA initiative aiming to produce a seamless, harmonized surface reflectance record from the Operational Land Imager (OLI) and Multi-Spectral Instrument (MSI) aboard Landsat-8 and Sentinel-2 remote sensing satellites, respectively. The HLS products are based on a set of algorithms to obtain seamless products from both sensors (OLI and MSI): atmospheric correction, cloud and cloud-shadow masking, geographic co-registration and common gridding, bidirectional reflectance distribution function normalization and bandpass adjustment. As of version 1.3, the HLS v1.3 data set covers 9.12 million km2 and spans from first Landsat-8 data (2013) to present. HLS products provide near-daily surface reflectance information with a common geometric framework, and are suitable for a variety of agricultural and vegetation monitoring tasks, including analysis of crop type, condition, and phenology.

Masek, Jeffrey↗

Overview of SatCORPS Satellite-Derived Cloud and Radiation Retrieval Coverage of ARM Domains

The ARM Climate Research Facility program maintains a number of climatically representative sites, which provide long-term cloud- and climate-monitoring records via ground-based instrumentation. These measurements provide a valuable data record over a localized region, but can be greatly enhanced by use of satellite monitoring. Satellite analyses over larger domains can provide parameters helpful for monitoring climate and evaluating models. The NASA/Langley Cloud group routinely derives such cloud and radiative parameters, from various geostationary and polar-orbiting satellite coverage over ARM sites; the group provides near-realtime analyses covering the 3 ARM fixed sites, as well as the GO-Amazon AMF site. This is accomplished by employing a suite of algorithms including VISST (Visible Infrared Solar Split-Window Technique), SIST (Solar Infrared Split-Window Technique), and SINT (Solar-infrared Infrared Near-Infrared Technique), now collectively called SatCORPS (Satellite Cloud Observations and Radiative Property retrieval System). An overview and catalog of SatCORPS-derived datasets processed for ASR, and available from both the ARM archive and the NASA/Langley Cloud group website, is provided. Specific improvements included in recently added datasets such as GO-Amazon and Azores are highlighted, including an improved cloud-detection mask, as well as improvements in derived Top-of-Atmosphere (TOA) SW albedoes and LW fluxes. New narrowband-to-broadband (NB-BB) fits and corrections for improved TOA fluxes are illustrated, including MTSAT-1/CERES Aqua NB-BB fits for the TWPICE field campaign, as well as new fits covering the Azores region which incorporate GERB TOA fluxes (Geostationery Earth Radiation Budget). Finally, validation of the reprocessed SatCORPS datasets is shown.

M M Khader↗

Adaptive Data Screening for Multi-Angle Polarimetric Aerosol and Ocean Color Remote Sensing Accelerated by Automatic Differentiation

Remote sensing measurements from multi-angle polarimeters (MAPs) contain rich aerosol microphysical property information, and these sensors have been used to perform retrievals in optically complex atmosphere and ocean systems. Previous studies have concluded that, generally, five moderately separated viewing angles in each spectral band provide sufficient accuracy for aerosol property retrievals, with performance gradually saturating as angles are added above that threshold. The Hyper-Angular Rainbow Polarimeter (HARP) instruments provide high angular sampling with a total of 90-120 unique angles across four bands, a capability developed mainly for liquid cloud retrievals. In practice, not all view angles are optimal for aerosol retrievals due to impacts of clouds, sun glint, and other impediments. The many viewing angles of HARP can provide resilience to these effects, if the impacted views are screened from the dataset, as the remaining views may be sufficient for successful analysis. In this study, we discuss how the number of available viewing angles impacts aerosol and ocean color retrieval uncertainties, as applied to two versions of the HARP instrument. AirHARP is an airborne prototype that was deployed in the ACEPOL field campaign, while HARP2 is an instrument in development for the upcoming NASA Plankton, Aerosol, Cloud, ocean Ecosystem (PACE) mission. Based on synthetic data, we find that a total of 20-30 angles across all bands (i.e. five to eight viewing angles per band) are sufficient to achieve good retrieval performance. Following from this result, we develop an adaptive multi-angle polarimetric data screening (MAPDS) approach to evaluate data quality by comparing measurements with their best-fitted forward model. The FastMAPOL retrieval algorithm is used to retrieve scene geophysical values, by matching an efficient, deep learning-based, radiative transfer emulator to observations. The data screening method effectively identifies and removes viewing angles affected by thin cirrus clouds and other anomalies, improving retrieval performance. This was tested with AirHARP data, and we found agreement with the High Spectral Resolution Lidar-2 (HSRL-2) aerosol data. The data screening approach can be applied to modern satellite remote sensing missions, such as PACE, where a large amount of multi-angle, hyperspectral, polarimetric measurements will be collected.

multi-angle polarimeter↗

First retrieval of AOD at fine-resolution over shallow and turbid coastal waters from MODIS

The widely used Moderate Resolution Imaging Spectroradiometer (MODIS) Dark-Target (DT) aerosol product fails to accurately retrieve Aerosol Optical Depth (AOD) over shallow and turbid Coastal Waters (CWs). To fill in gaps, and to improve land to ocean AOD continuity, we developed a coastal water retrieval algorithm at a spatial resolution of 1 km (CW-1km). CW-1km relies on observed top-of-atmosphere reflectance at 2.1 μm (ρ2.1), both to derive AOD and to perform a spatial variation test that enhances the existing DT masks for clouds and land. We show that the CW-1km improves spatial continuity of AOD between land, coast, and open ocean, while also increasing AOD product availability by 47.0%. Comparing with 15 years of marine aerosol network measurements, CW-1km AODs are validated to have a normalized mean bias of 1.0%, which is much smaller than 17.6% for the original DT product.

MODIS AOD↗

Cloud Tolerance of Remote-Sensing Technologies to Measure Land Surface Temperature

Conventional methods to estimate land surface temperature (LST) from space rely on the thermal infrared(TIR) spectral window and is limited to cloud-free scenes. To also provide LST estimates during periods with clouds, a new method was developed to estimate LST based on passive microwave(MW) observations. The MW-LST product is informed by six polar-orbiting satellites to create a global record with up to eight observations per day for each 0.25resolution grid box. For days with sufficient observations, a continuous diurnal temperature cycle (DTC) was fitted. The main characteristics of the DTC were scaled to match those of a geostationary TIR-LST product. This paper tests the cloud tolerance of the MW-LST product. In particular, we demonstrate its stable performance with respect to flux tower observation sites (four in Europe and nine in the United States), over a range of cloudiness conditions up to heavily overcast skies. The results show that TIR based LST has slightly better performance than MW-LST for clear-sky observations but suffers an increasing negative bias as cloud cover increases. This negative bias is caused by incomplete masking of cloud-covered areas within the TIR scene that affects many applications of TIR-LST. In contrast, for MW-LST we find no direct impact of clouds on its accuracy and bias. MW-LST can therefore be used to improve TIR cloud screening. Moreover, the ability to provide LST estimates for cloud-covered surfaces can help expand current clear-sky-only satellite retrieval products to all-weather applications.

polar-orbiting↗

High Ice Water Concentrations in the 19 August 2015 Coastal Mesoconvective System

During August 2015, NASA's DC-8 research aircraft was flown into High Ice Water Content (HIWC) events as part of a three-week campaign to collect airborne radar data and to obtain measurements from microphysical probes. Goals for this flight campaign included improved characterization of HIWC events, especially from an airborne radar perspective. This paper focuses on one of the flight days, in which a coastal mesoscale convective system (MCS) was investigated for HIWC conditions. The system appears to have been maintained by bands of convection flowing in from the Gulf of Mexico. These convective bands were capped by a large cloud canopy, which masks the underlying structure if viewed from an infrared sensing satellite. The DC-8 was equipped with an IsoKinetic Probe that measured ice concentrations of up to 2.3 g m(exp -3) within the cloud canopy of this system. Sustained measurements of ice crystals with concentrations exceeding 1 g m(exp -3) were encountered for up to ten minutes of flight time. Airborne Radar reflectivity factors were found to be weak within these regions of high ice water concentrations, suggesting that Radar detection of HIWC would be a challenging endeavor. This case is then investigated using a three-dimensional numerical cloud model. Profiles of ice water concentrations and radar reflectivity factor demonstrate similar magnitudes and scales between the flight measurements and model simulation. Also discussed are recent modifications to the numerical model's ice-microphysics that are based on measurements during the flight campaign. The numerical model and its updated ice-microphysics are further validated with a simulation of a well-known case of a supercell hailstorm measured during the Cooperative Convective Precipitation Experiment. Differences in HIWC between the continental supercell and the coastal MCS are discussed.

Proctor, Fred H.↗

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↗