We developed a global, 30-m resolution dataset of percent tree cover by rescaling the 250-m MOderate-resolution Imaging Spectroradiometer (MODIS) Vegetation Continuous Fields (VCF) Tree Cover layer using circa- 2000 and 2005 Landsat images, incorporating the MODIS Cropland Layer to improve accuracy in agricultural areas. Resulting Landsat-based estimates maintained consistency with the MODIS VCF in both epochs (RMSE =8.6% in 2000 and 11.9% in 2005), but showed improved accuracy in agricultural areas and increased discrimination of small forest patches. Against lidar measurements, the Landsat-based estimates exhibited accuracy slightly less than that of the MODIS VCF (RMSE=16.8% for MODIS-based vs. 17.4% for Landsat-based estimates), but RMSE of Landsat estimates was 3.3 percentage points lower than that of the MODIS data in an agricultural region. The Landsat data retained the saturation artifact of the MODIS VCF at greater than or equal to 80% tree cover but showed greater potential for removal of errors through calibration to lidar, with post-calibration RMSE of 9.4% compared to 13.5% in MODIS estimates. Provided for free download at the Global Land Cover Facility (GLCF) website (www.landcover.org), the 30-m resolution GLCF tree cover dataset is the highest-resolution multi-temporal depiction of Earth's tree cover available to the Earth science community.

We introduce and test a new method to detect annual forest cover loss from time series estimates of percent tree cover. Our approach is founded on two realistic assumptions: (1) land cover disturbances are rare events over large geographic areas that occur within a short time frame; and (2) spatially discrete land cover disturbances are continuous processes over time. Applying statistically rigorous algorithms, we first detect disturbance pixels as outliers of an underlying chi-square distribution. Then, we fit nonlinear, logistic curves for each identified change pixel to simultaneously characterize the magnitude and timing of the disturbance. Our method is applied using the yearly Vegetation Continuous Fields (VCF) tree cover product from Moderate Resolution Imaging Spectroradiometer (MODIS), and the resulting disturbance-year estimates are evaluated using a large sample of Landsat-based forest disturbance data. Temporal accuracy is ~65% at 250-m, annual resolution and increases to >85% when temporal resolution is relaxed to ±1 yr. The r² of MODIS VCF-based disturbance rates against Landsat ranges from 0.7 to 0.9 at 5-km spatial resolution. The general approach developed in this study can be potentially applied at a global scale and to other land cover types characterized as continuous variables from satellite data.

Monitoring current forest characteristics in the taiga−tundra ecotone (TTE) at multiple scales is critical for understanding its vulnerability to structural changes. A 30 m spatial resolution Landsat-based tree canopy cover map has been calibrated and validated in the TTE with reference tree cover data from airborne LiDAR and high resolution spaceborne images across the full range of boreal forest tree cover. This domain-specific calibration model used estimates of forest height to determine reference forest cover that best matched Landsat estimates. The model removed the systematic under-estimation of tree canopy cover >80% and indicated that Landsat estimates of tree canopy cover more closely matched canopies at least 2 m in height rather than 5 m. The validation improved estimates of uncertainty in tree canopy cover in discontinuous TTE forests for three temporal epochs (2000, 2005, and 2010) by reducing systematic errors, leading to increases in tree canopy cover uncertainty. Average pixel-level uncertainties in tree canopy cover were 29.0%, 27.1% and 31.1% for the 2000, 2005 and 2010 epochs, respectively. Maps from these calibrated data improve the uncertainty associated with Landsat tree canopy cover estimates in the discontinuous forests of the circumpolar TTE.

Measurements of solar radiation reflected from Earth’s surface are the basis for calculating albedo, vegetation indices, and other terrestrial attributes. However, the “bi-directional” geometry of illumination and viewing (i.e., the Bi-directional Reflectance Distribution Function (BRDF)) impacts reflectance and all variables derived or estimated based on these data. The recently released MODIS BRDF/Albedo Model Parameters (MCD43A1 Collection 6) dataset enables retrieval of directional reflectance at arbitrary solar and viewing angles, potentially increasing precision and comparability of data collected under different illumination and observation geometries. We quantified the ability of MCD43A1 Collection 6 for retrieving directional reflectance and compared the daily Collection 6 retrievals to those of MCD43A1 Collection 5, which are retrieved on an eight-day basis. Correcting MODIS-based estimates of surface reflectance from the illumination and viewing geometry of the Terra satellite (MOD09GA) to that of the MODIS Aqua (MYD09GA) overpass, as well as MCD43A4 Collection 6 and Landsat-5 TM images show that the BRDF correction of MCD43A1 Collection 6 results in greater consistency among datasets, with higher R² (0.63–0.955), regression slopes closer to unity (0.718–0.955), lower root mean squared difference (RMSD) (0.422–3.142), and lower mean absolute error (MAE) (0.282–1.735) compared to the Collection 5 data. Smaller levels of noise (observed as high-frequency variability within the time series) in MCD43A1 Collection 6 in comparison to Collection 5 corroborates the improvement of BRDF parameters time series. These results corroborates that the daily MCD43A1 Collection 6 product represents the anisotropy of surface features and results in more precise directional reflectance derivation at any solar and viewing geometry than did the previous Collection 5.