Recent production of land surface anisotropy, diffuse bihemispherical (white-sky) albedo, and direct-beam directional hemispherical (black-sky) albedo from observations acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) instruments aboard the National Aeronautics and Space Administration's Terra and Aqua satellite platforms have provided researchers with unprecedented spatial, spectral, and temporal information on the land surface's radiative characteristics. Cloud cover, which curtails retrievals, and the presence of ephemeral and seasonal snow limit the snow-free data to approximately half the global land surfaces on an annual equal-angle basis. This precludes the MOD43B3 albedo products from being used in some remote sensing and ground-based applications, climate models, and global change research projects. An ecosystem-dependent temporal interpolation technique is described that has been developed to fill missing or seasonally snow-covered data in the official MOD43B3 albedo product. The method imposes pixel-level and local regional ecosystem-dependent phenological behavior onto retrieved pixel temporal data in such a way as to maintain pixel-level spatial and spectral detail and integrity. The phenological curves are derived from statistics based on the MODIS MOD12Q1 IGBP land cover classification product geolocated with the MOD43B3 data. The resulting snow-free value-added products provide the scientific community with spatially and temporally complete global white- and black-sky surface albedo maps and statistics. These products are stored on 1-min and coarser resolution equal-angle grids and are computed for the first seven MODIS wavelengths, ranging from 0.47-2.1 μm and for three broadband wavelengths 0.3-0.7, 0.3-5.0, and 0.7-5.0 μm.