We describe and validate a model that retrieves fractional snow-covered area and the grain size and albedo of that snow from surface reflectance data (product MOD09GA) acquired by NASA's Moderate Resolution Imaging Spectroradiometer (MODIS). The model analyzes the MODIS visible, near infrared, and shortwave infrared bands with multiple endmember spectral mixtures from a library of snow, vegetation, rock, and soil. We derive snow spectral endmembers of varying grain size from a radiative transfer model specific to a scene's illumination geometry; spectra for vegetation, rock, and soil were collected in the field and laboratory. We validate the model with fractional snow cover estimates from Landsat Thematic Mapper data, at 30 m resolution, for the Sierra Nevada, Rocky Mountains, high plains of Colorado, and Himalaya. Grain size measurements are validated with field measurements during the Cold Land Processes Experiment, and albedo retrievals are validated with in situ measurements in the San Juan Mountains of Colorado. The pixelweighted average RMS error for snow-covered area across 31 scenes is 5%, ranging from 1% to 13%. The mean absolute error for grain size was 51 µm and the mean absolute error for albedo was 4.2%. Fractional snow cover errors are relatively insensitive to solar zenith angle. Because MODSCAG is a physically based algorithm that accounts for the spatial and temporal variation in surface reflectances of snow and other surfaces, it is capable of global snow cover mapping in its more computationally efficient, operational mode.

Forest canopies influence the proportion of the land surface that is visible from above, or the viewable gap fraction (VGF). The VGF limits the amount of information available in satellite data about the land surface, such as snow cover in forests. Efforts to recover fractional snow cover from the spectral mixture analysis model Moderate Resolution Imaging Spectroradiometer (MODIS) snow-covered area and grain size (MODSCAG) indicate the importance of view angle effects in forested landscapes. The VGF can be estimated using both hemispherical photos and forest canopy models. For a set of stands in the Cold Land Field Processes Experiment (CLPX) sites in the Fraser Experimental Forest in Colorado, the convergence of both measurements and models of the VGF as a function of view angle supports the idea that VGF can be characterized as a function of forest properties. A simple geometric optical (GO) model that includes only between-crown gaps can capture the basic shape of the VGF as a function of view zenith angle. However, the GO model tends to underestimate the VGF compared with estimates derived from hemispherical photos, particularly at high view angles. The use of a more complicated geometric opticalradiative transfer (GORT) model generally improves estimates of the VGF by taking into account withincrown gaps. Small footprint airborne lidar data are useful for mapping forest cover and height, which makes the parameterization of the GORT model possible over a landscape. Better knowledge of the angular distribution of gaps in forest canopies holds promise for improving remote sensing of snow cover fraction.