The multi-decadal Landsat data record is a unique tool for global land cover and land use change analysis. However, the large volume of the Landsat image archive and inconsistent coverage of clear-sky observations hamper land cover monitoring at large geographic extent. Here, we present a consistently processed and temporally aggregated Landsat Analysis Ready Data produced by the Global Land Analysis and Discovery team at the University of Maryland (GLAD ARD) suitable for national to global empirical land cover mapping and change detection. The GLAD ARD represent a 16-day time-series of tiled Landsat normalized surface reflectance from 1997 to present, updated annually, and designed for land cover monitoring at global to local scales. A set of tools for multi-temporal data processing and characterization using machine learning provided with GLAD ARD serves as an end-to-end solution for Landsat-based natural resource assessment and monitoring. The GLAD ARD data and tools have been implemented at the national, regional, and global extent for water, forest, and crop mapping. The GLAD ARD data and tools are available at the GLAD website for free access.

A comprehensive record of the land use history of South America reveals that human impact on land increased by 60% since 1985.

Lack of national data on water-related ecosystems is a major challenge to achieving the Sustainable Development Goal (SDG) 6 targets by 2030. Monitoring surface water extent, wetlands, and water quality from space can be an important asset for many countries in support of SDG 6 reporting. We demonstrate the potential for Earth observation (EO) data to support country reporting for SDG Indicator 6.6.1, ‘Change in the extent of water-related ecosystems over time’ and identify important considerations for countries using these data for SDG reporting. The spatial extent of water-related ecosystems, and the partial quality of water within these ecosystems is investigated for seven countries. Data from the Moderate Resolution Imaging Spectroradiometer (MODIS) and Landsat 5, 7, and 8 with Shuttle Radar Topography Mission (SRTM) are used to measure surface water extent at 250 m and 30 m spatial resolution, respectively, in Cambodia, Jamaica, Peru, the Philippines, Senegal, Uganda, and Zambia. The extent of mangroves is mapped at 30 m spatial resolution using Landsat 8 Operational Land Imager (OLI), Sentinel-1, and SRTM data for Jamaica, Peru, and Senegal. Using Landsat 8 and Sentinel 2A imagery, total suspended solids and chlorophyll-a are mapped over time for a select number of large surface water bodies in Peru, Senegal, and Zambia. All of the EO datasets used are of global coverage and publicly available at no cost. The temporal consistency and long time-series of many of the datasets enable replicability over time, making reporting of change from baseline values consistent and systematic. We find that statistical comparisons between different surface water data products can help provide some degree of confidence for countries during their validation process and highlight the need for accuracy assessments when using EO-based land change data for SDG reporting. We also raise concern that EO data in the context of SDG Indicator 6.6.1 reporting may be more challenging for some countries, such as small island nations, than others to use in assessing the extent of water-related ecosystems due to scale limitations and climate variability. Country-driven validation of the EO data products remains a priority to ensure successful data integration in support of SDG Indicator 6.6.1 reporting. Multi-country studies such as this one can be valuable tools for helping to guide the evolution of SDG monitoring methodologies and provide a useful resource for countries reporting on water-related ecosystems. The EO data analyses and statistical methods used in this study can be easily replicated for country-driven validation of EO data products in the future.

While publicly available, cost-free coarse and medium spatial resolution satellite data such as MODIS and Landsat perform well in characterizing industrial cropping systems, commercial high spatial resolution satellite data are often preferred alternative for fine scale land tenure agricultural systems such as found in Pakistan. In this article, we integrated commercial 5 m spatial resolution RapidEye and free 30 m Landsat imagery in characterizing winter wheat in Punjab province, Pakistan. Specifically, we used 5 m spatial resolution RapidEye imagery from peak of the winter wheat growing season to derive training data for the characterization of time-series Landsat data. After co-registration, each RapidEye image was classified into wheat/no wheat labels at the 5 m resolution and then aggregated as percent cover to 30 m Landsat grid cells. We produced four maps, two using RapidEye derived continuous training data (of percent wheat cover) as input to a regression tree model, and two using RapidEye derived categorical training data as input to a classification tree model. From the RapidEye-derived 30 m continuous training data, we derived Map 1 as percent wheat per pixel, and Map 2 as binary wheat/no wheat classification derived using a 50% threshold applied to Map 1. To create the categorical wheat/no wheat training data, we first converted the continuous training data to a wheat/no wheat classification, and then used these categorical RapidEye training data to produce a categorical wheat map from the Landsat data. Two methods for categorizing the training data were used. The first method used a 50% wheat/no wheat threshold to produce Map 3, and the second method used only pure wheat (≥75% cover) and no wheat (≤25% cover) training pixels to produce Map 4. The approach of Map 4 is analogous to a standard method in which whole, pure, high-confidence training pixels are delineated. We validated the wheat maps with field data collected using a stratified, two-stage cluster design. Accuracy of the maps produced from the percent cover training data (Map 1 and Map 2) was not substantially better than the accuracy of the maps produced from the categorical training data as all methods yielded similar overall accuracies (±standard error): 88% (±4%) for Map 1, 90% (±4%) for Map 2, 90% (±4%) for Map 3, and 87% (±4%) for Map 4. Because the percent cover training data did not produce significantly higher accuracies, sub-pixel training data are not required for winter wheat mapping in Punjab. Given sufficient expertise in supervised classification model calibration, freely available Landsat data are sufficient for crop mapping in the fine-scale land tenure system of Punjab. For winter wheat mapping in Punjab and other like landscapes, training data for supervised classification may be collected directly from Landsat images without the need for high resolution reference imagery.

Direct human action, principally land use expansion, and natural dynamics, such as fire and drought, drive global land change. Here we present a global land change monitoring system, DIST-ALERT, that rapidly tracks vegetation loss anomalies with 30 m resolution using imagery from Landsat 8/9 and Sentinel-2A/B/C satellites. The alerts capture agricultural expansion, urbanization, logging, mining, fire, drought, landslides, and other dynamics, but without attribution. Identified through a probability sample, 2023 anthropogenic land use conversions totaled 28.6 ± 7.6 Mha (±standard error), half of which replaced long-lived or secondary natural vegetation. Fires resulting in land cover conversion totaled 14.9 ± 4.3 Mha (±standard error). Combined, these dynamics equal 0.3% of the global land surface, equivalent to the area of the state of California. Annual DIST-ALERT summaries of land use expansion and climate-driven land change can serve as a future long-term global environmental data record.