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.

Surface water is of great importance to ecosystems and economies. Crucial to understanding hydrological variability and its relationships to human activities at large scales, open-access satellite datasets and big-data computational methods are now enabling the global mapping of the distribution and changes of inland water over time. A machine-learning algorithm, previously used only to map water at single points in time, was applied over 16 years of the USGS Landsat archive to detect and map surface water over central Asia from 2000 to 2015 at a 30-m, monthly resolution. The resulting dataset had an overall classification accuracy of 99.59% (±0.32% standard error), 98.24% (±1.02%) user’s accuracy, and 87.12% (±3.21%) producer’s accuracy for water class. This study describes the temporal extension of the algorithm and the application of the dataset to present patterns of regional surface water cover and change. The findings indicate that smaller water bodies are dramatically changing in two specific ecological zones: the Kazakh Steppe and the Tian Shan Montane Steppe and Meadows. Both the maximum and minimum extent of water bodies have decreased over the 16-year period, but the rate of decrease of the maxima was double that of the minima. Coverage decreased in each month from April to October, and a significant decrease in water area was found in April and May. These results indicate that the dataset can provide insights into the behavior of surface water across central Asia through time, and that the method can be further developed for regional and global applications.

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.

Stand age distribution is an important descriptor of boreal forest structure, which is directly linked to many ecosystem processes including the carbon cycle, the land–atmosphere interaction and ecosystem services, among others. Almost half of the global boreal biome is located in Russia. The vast extent, remote location, and limited accessibility of Russian boreal forests make remote sensing the only feasible approach to characterize these forests to their full extent. A wide variety of satellite observations are currently available to monitor forest change and infer its structure; however, the period of observations is mostly limited to the 2000s era. Reconstruction of wall-to-wall maps of stand age distribution requires merging longer-term site observations of forest cover change available at the Landsat scale at a subset of locations in Russia with the wall-to-wall coverage available from coarse resolution satellites since 2000. This paper presents a dataset consisting of a suite of multi-year forest disturbance samples and samples of undisturbed forests across Russia derived from Landsat Thematic Mapper and Enhanced Thematic Mapper Plus images from 1985 to 2000. These samples provide crucial information regarding disturbance history in selected regions across the Russian boreal forest and are designed to serve as a training and/or validation dataset for coarse resolution data products. The overall accuracy and Kappa coefficient for the entire sample collection was found to be 83.98% and 0.83%, respectively. It is hoped that the presented dataset will benefit subsequent studies on a variety of aspects of the Russian boreal forest, especially in relation to the carbon budget and climate.