Many novel research algorithms have been developed to analyze urban heat island (UHI) and UHI regional impacts (UHIRIP) with remotely sensed thermal data tables. We present a comprehensive review of some important aspects of UHI and UHIRIP studies that use remotely sensed thermal data, including concepts, datasets, methodologies, and applications. We focus on reviewing progress on multi-sensor image selection, preprocessing, computing, gap filling, image fusion, deep learning, and developing new metrics. This literature review shows that new satellite sensors and valuable methods have been developed for calculating land surface temperature (LST) and UHI intensity, and for assessing UHIRIP. Additionally, some of the limitations of using remotely sensed data to analyze the LST, UHI, and UHI intensity are discussed. Finally, we review a variety of applications in UHI and UHIRIP analyses. The assimilation of time-series remotely sensed data with the application of data fusion, gap filling models, and deep learning using the Google Cloud platform and Google Earth Engine platform also has the potential to improve the estimation accuracy of change patterns of UHI and UHIRIP over long time periods.

Remotely monitoring changes in central U.S. grasslands is challenging because these landscapes tend to respond quickly to disturbances and changes in weather. Such dynamic responses influence nutrient cycling, greenhouse gas contributions, habitat availability for wildlife, and other ecosystem processes and services. Traditionally, coarse-resolution satellite data acquired at daily intervals have been used for monitoring. Recently, the harmonized Landsat-8 and Sentinel-2 (HLS) data increased the temporal frequency of the data. Here we investigated if the increased data frequency provided adequate observations to characterize highly dynamic grassland processes. We evaluated HLS data available for 2016 to (1) determine if data from Sentinel-2 contributed to an improvement in characterizing landscape processes over Landsat-8 data alone, and (2) quantify how observation frequency impacted results. Specifically, we investigated into estimating annual vegetation phenology, detecting burn scars from fire, and modeling within-season wetland hydroperiod and growth of aquatic vegetation. We observed increased sensitivity to the start of the growing season (SOST) with the HLS data. Our estimates of the grassland SOST compared well with ground estimates collected at a phenological camera site. We used the Continuous Change Detection and Classification (CCDC) algorithm to assess if the HLS data improved our detection of burn scars following grassland fires and found that detection was considerably influenced by the seasonal timing of the fires. The grassland burned in early spring recovered too quickly to be detected as change events by CCDC; instead, the spectral characteristics following these fires were incorporated as part of the ongoing time-series models. In contrast, the spectral effects from late-season fires were detected both by Landsat-8 data and HLS data. For wetland-rich areas, we used a modified version of the CCDC algorithm to track within-season dynamics of water and aquatic vegetation. The addition of Sentinel-2 data provided the potential to build full time series models to better distinguish different wetland types, suggesting that the temporal density of data was sufficient for within-season characterization of wetland dynamics. Although the different data frequency, in both the spatial and temporal dimensions, could cause inconsistent model estimation or sensitivity sometimes; overall, the temporal frequency of the HLS data improved our ability to track within-season grassland dynamics and improved results for areas prone to cloud contamination. The results suggest a greater frequency of observations, such as from harmonizing data across all comparable Landsat and Sentinel sensors, is still needed. For our study areas, at least a 3-day revisit interval during the early growing season (weeks 14–17) is required to provide a >50% probability of obtaining weekly clear observations.

The U.S. Geological Survey’s Land Change Monitoring, Assessment, and Projection (LCMAP) initiative involves detecting changes in land cover, use, and condition with the goal of producing land change information to improve the understanding of the Earth system and provide insights on the impacts of land surface change on society. The change detection method ingests all available high-quality data from the Landsat archive in a time series approach to identify the timing and location of land surface change. Annual thematic land cover maps are then produced by classifying time series models. In this paper, we describe the optimization of the classification method used to derive the thematic land cover product. We investigated the influences of auxiliary data, sample size, and training from different sources such as the U.S. Geological Survey’s Land Cover Trends project and National Land Cover Database (NLCD 2001 and NLCD 2011). The results were evaluated and validated based on independent data from the training dataset. We found that refining the auxiliary data effectively reduced artifacts in the thematic land cover map that are related to data availability. We improved the classification accuracy and stability considerably by using a total of 20 million training pixels with a minimum of 600,000 and a maximum of 8 million training pixels per class within geographic windows consisting of nine Analysis Ready Data tiles (450 by 450 km²). Comparisons revealed that the NLCD 2001 training data delivered the best classification accuracy. Compared to the original LCMAP classification strategy used for early evaluation (e.g., Trends training data, 20,000 samples), the optimized classification strategy improved the annual land cover map accuracy by an average of 10%.

The recently released Landsat analysis ready data (ARD) over the United States provides the opportunity to investigate landscape dynamics using dense time series observations at 30-m resolution. However, the dataset often contains data gaps (or missing data) because of cloud contamination or data acquisition strategy, which result in different capabilities for seasonality modeling. We present a new algorithm that focuses on data gap filling using clear observations from orbit overlap regions. Multiple linear regression models were established for each pixel time series to estimate stable predictions and uncertainties. The model’s training data came from stratified random samples based on the time series similarity between the pixel and data from the overlap regions. The algorithm was first evaluated using four tiles (5000 × 5000 30-m pixels for each tile) from 2018 land surface temperature data (LST) in Atlanta, Georgia. The accuracy was assessed using randomly masked clear observations with an average Root Mean Square Error (RMSE) of 3.88 and an average bias of −0.37, which were comparable to the product accuracy. We also applied the method on ARD surface reflectance bands at Fairbanks, Alaska. The accuracy assessment suggested a majority RMSE of less than 0.04 and a bias of less than 0.0023. The gap-filled time series can be of help for reliable seasonal modeling and reducing artifacts related to data availability. This approach can also be applied to other datasets, vegetation indexes, or spectral reflectance bands of other sensors.

The land carbon sink of the conterminous United States was evaluated using a bottom-up modeling framework and 30-meter land change data from 1985 to 2020. This cross-scale, cross-landscape, and cross-system approach tracked fractional land cover changes and applied regional model calibration. Results show average terrestrial and aquatic carbon sinks of +110 ± 37 and +19 ± 0.5 teragrams of carbon per year, respectively. The terrestrial carbon sink, showing no clear trend, peaked in the 1990s, with more years as a carbon source since 2000, contradicting recent national and global studies. Land change had the largest impact (−70 ± 5.5 teragrams of carbon per year), exceeding impacts of climate (−33 ± 48 teragrams of carbon per year), wildfire (−7.7 ± 2.4 teragrams of carbon per year), and erosion transport (−1.9 ± 0.13 teragrams of carbon per year). The positive CO ₂ fertilization effect (+69 ± 12 teragrams of carbon per year) was insufficient to maintain the carbon sink strength. Our framework reveals key paths of carbon loss, with implications for carbon budget and energy policies in the United States and beyond.

Global observations at 30-m ground sampling distance (GSD) are now possible at a cadence of 1-3 days by combining Landsat 8 and 9 with Sentinel-2A and -2B satellites. Previous studies characterizing pixel-level Landsat-class measurement frequency used data from different sources but offered little information on observation availability after rigorous quality screening. This study examined the coverage frequency of HLS V2.0 data for 2022, the first year all four satellites data were available. These data have had quality control filtering and harmonization, and therefore reflect the spatial-temporal distribution of usable observations. On average, HLS data provide observations every 1.6 days at the global scale, and 2.2 days in the data-scarce tropical regions, regardless of cloud cover. The global mean and median cloud-free observations were 69 and 64, respectively. The frequency of good-quality observations varies geographically and seasonally due to changes in satellite swath overlap, cloud frequency, and solar illumination. High latitudes (>~75°N) exhibit the highest number of cloud-free observations between March and September. However, data are unavailable during winter months due to low solar elevation angles and boreal regions have a lower number of clear observations in the summer months. The tropical regions have the lowest number of clear observations. More frequent HLS observations could improve terrestrial monitoring. We mapped the monthly and weekly number of clear observations globally to show where HLS data could support monthly or sub-weekly time series applications.