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.

Lake extent is an indicator of water capacity as well as the aquatic ecological and environmental conditions. Due to the small sizes and rapid water dynamics, monitoring the extent of small lakes fluctuating between 2.5 and 30 km² require observations with both high spatial and temporal resolutions. The paper applied an improved surface reflectance (SR) downscaling method (i.e., IMAR (Improved Modified Adaptive Regression model)) to downscale the daily SR acquired by the Moderate Resolution Imaging Spectroradiometer (MODIS) Terra platform to a consistent 250-m resolution, and derived monthly water extent of four small lakes in the Tibetan Plateau (Longre Co, Ayonggongma Co, Ayonggama Co, and Ayongwama Co)) from 2000 to 2014. Using Landsat ETM+ acquired on the same date, the downscaled MODIS SR and identified water extent were compared to the original MODIS, observations downscaled using an early SR downscaling method (MAR (Modified Adaptive Regression model)) and Wavelet fusion. The results showed IMAR achieved the highest correlation coefficients (R²) (0.89–0.957 for SR and 0.79–0.933 for water extent). The errors in the derived water extents were significantly decreased comparing to the results of MAR and Wavelet fusion, and lakes morphometry of IMAR is more comparable to Landsat results. The detected lake extents dynamic between 2000 and 2014 were analyzed using the trend and season decomposition model (BFAST), indicating an increasing trend after 2005, and it likely had higher correlations with temperature and precipitation variation in the Tibetan region (R²: 0.598–0.728 and 0.61–0.735, respectively).

Inland surface water is essential to terrestrial ecosystems and human civilization. Accurate mapping of surface water dynamic is vital for both scientific research and policy-driven applications. MODIS provides twice observation per day, making it perfect for monitoring temporal water dynamic. Although MODIS provides two bands at 250 m resolution, accurately deriving water area always depends on observations from the spectral bands with 500 m resolution, which limits its discrimination ability over small lakes and rivers. The paper presents an automated method for downscaling the 500 m MODIS surface reflectance (SR) to 250 m to improve the spatial discrimination of water body extraction. The method has been tested at Co Ngoin and Co Bangkog in Qinghai-Tibet plateau. The downscaled SR and the derived water bodies were compared to SR and water body mapped from Landsat-7 ETM+ images were acquired on the same date. Consistency metrics were calculated to measure their agreement and disagreement. The comparisons indicated that the downscaled MODIS SR showed significant improvement over the original 500 m observations when compared with Landsat-7 ETM+ SR, and both commission and omission errors were reduced in the derived 250 m water bodies.