Urban Heat Island and Its Regional Impacts Using Remotely Sensed Thermal Data—A Review of Recent Developments and Methodology

: 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 ﬁlling, 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 ﬁlling 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.


Introduction
Urbanization is known to have substantial impacts on landscapes and ecosystems [1][2][3][4], and urban inhabitants are expected to reach 70% of the world population by 2050 [5]. Moreover, the nature of urban development has been changing from a single city model to a group of cities (urban agglomeration) worldwide. Urban heat island (UHI), urbanization, and climate change are increasingly interconnected, resulting in several environmental consequences (such as heat stress, biodiversity loss, fire risk, warming water due to run off, and diminished air quality) at both local and regional levels [2,[6][7][8][9]. Such UHI related impacts are also called UHI regional impacts (UHIRIP). Generally, UHI research includes data from two major sources: air temperature data that are observed by weather or climate stations and remotely sensed data to observe UHI through land surface temperature. Before the availability of remotely sensed data, UHI was widely observed in the field, with the first scientific observation of UHI in 1833 [10]. Field observations of UHI continue to be a critical source of training and validation data [11,12]. These observations, along with modeling studies, continue to help unravel the factors that are responsible for UHI development, and are providing a basis for the development and application of sustainable adaptation strategies. Communicating scientific knowledge quickly and effectively of UHI and UHIRIP to architects, engineers, scientists, and planners could help inform urban Table 1. Example of main reviews, bibliographies, and summaries on UHI and UHIRIP using remotely sensed data.

Reference Topics Sensors Measurements
Hall et al. [11] Satellite remote sensing of surface energy balance success, failures, and unresolved issues in field experiment (FIFE)

Landsat, SPOT Thermal
Gallo et al. [13] Assessment of urban heat islands: A satellite perspective AVHRR, Landsat MSS Thermal Voogt and Oke [6] Thermal remote sensing of urban climates Multiple, review Thermal Weng and Larson [49] Satellite remote sensing of urban heat islands: current practice and prospects Multiple, review Thermal Jiang et al. [50] Land surface emissivity retrieval from combined mid-infrared and thermal infrared data of MSG-SEVIRI

Meteosat Second Generation (MSG)
Spinning Enhanced Visible and Infrared Imager (SEVIRI) Kalma et al. [51] Estimating land surface evaporation: A review of methods using remotely sensed surface temperature data Multiple, review Thermal Racoviteanu et al. [52] Optical remote sensing of glacier characteristics: a review focusing on the Himalaya

ASTER Indices
Rizwan et al. [53] A review on the generation, determination, and mitigation of urban heat island Review Determination of UHI Weng [7] Thermal infrared remote sensing for urban climate and environmental studies: Methods, applications, and trends Multiple, review Thermal Bowler et al. [31] Urban greening to cool towns and cities: A systematic review of the empirical evidence Review Synthesis analysis Sailor [54] A review of methods for estimating anthropogenic heat and moisture emissions in the urban environment Review Bibliometric profile Li et al. [55] Satellite-derived land surface temperature: current status and perspectives Multiple, review Thermal

Reference Topics Sensors Measurements
Ngie et al. [56] Assessment of urban heat island using satellite remotely sensed imagery: A review Multiple, review Thermal Rasul et al. [16] A review of remote sensing of urban heat and cool islands Multiple, review Thermal Huang and Lu [57] Urban heat island research from 1991 to 2015: A bibliometric analysis Review Bibliometric profile Zhang et al. [58] A bibliometric profile of the remote sensing open access journal published by MDPI between 2009 and 2018 Multiple, review Bibliometric profile Deilami et al. [59] Urban heat island effect: A systematic review of spatio-temporal factors, data, methods, and mitigation measures Multiple, review Thermal Zhou et al. [60] Satellite remote sensing of surface urban heat islands: Progress, challenges, and perspectives Multiple, review Thermal Becker and Zhao-Liang [61] Surface temperature and emissivity at various scales: definition, measurement, and related problems Multiple, review Thermal (surface temperature and emissivity) Dash et al. [62] Land surface temperature and emissivity estimation from passive sensor data: theory and practice current trends Multiple, review Thermal (surface temperature and emissivity) Table 2. Examples of research publications investigating UHI and UHIRIP using remotely sensed data.

Urban Heat Island
UHI studies have been conducted for over 200 years, since the first conceptualization by Luke Howard in 1818 [98]. Generally, an urban heat island (UHI) is an urban area or metropolitan area that is significantly warmer than its surrounding rural areas because of human activities. The temperature difference is usually greater at night than during the day and is most apparent when winds are weak. Some research [99,100] shows that the annual mean air temperature of a city with 1 million people or more can be 1-3 • C warmer than its surroundings. In the evening, the difference can be as high as 12 • C. Heat islands can affect Land 2021, 10, 867 5 of 30 communities by increasing the summertime peak energy demand (such as air conditioning costs), air pollution and greenhouse gas emissions, heat-related illness, and mortality, and decreasing water quality and ecosystem security. Higher temperature "domes" are created over an urban or industrial areas by hot air layers forming at building-top or chimney-top level. This dome is about 5 • C to 7 • C warmer than the air above it and the ground level temperature, and can trap all polluting emissions within its confines (see also temperature inversion [53,57]).
The large amount of heat generated from urban structures and pavements, as they absorb and re-radiate solar radiation, as well as the heat from other anthropogenic sources, are the main causes of UHI. These heat sources increase the temperatures of an urban area compared with its surroundings, which is known as UHI intensity (UHII). Traditionally, regardless of the methodology employed, whether it refers to (1) differences between two fixed observatories, one urban and another peripheral or non-urban; (2) mobile urban transects; or (3) remote sensing analysis, UHII provides a value of thermal differences between contrasted points, sectors, or areas, one urban and another that could be termed non-urban. Thus, the intensity of the UHI is seen in the temperature difference expressed at a given time between the hottest sector (areas) of the city and the surrounding non-urban space. The intensity of the heat island is the simplest and most quantitative indicator of the thermal modification imposed by the city upon the territory in which it is situated and of its relative warming in relation to the surrounding rural environment. The intensity could be defined for various time scales and geographical locations [101,102].

The Study of the Spatial Structure of Urban Thermal Patterns, Change Dynamics, and Their Relation to Urban Surface Characteristics
Based on the fractional theory of ecology [103,104], the spatial structure of urban thermal patterns and temporal change dynamics can be studied in two and three dimensions. Figure 1 shows an example of the UHI and UHIRIP profile in Sioux Falls, South Dakota, USA, and the surrounding area, derived from Landsat ARD LST over different land cover classes [97]. The study of temporal change in UHI can include multiple scales of change, including daily, day and night, monthly, seasonal, yearly, and long-term time series. The physical mechanisms driving UHI are well documented [28]. UHIRIP may be described in multiple ways with various methodological approaches to investigate each type; specifically, it can impact the ground, the surface, and various heights in the air [105,106] at a regional scale. Different pictures arise for each type of UHI when measured by different methods. Tam et al. [107] suggested that the magnitude of total change in day to day temperature variability can be used to decide a suitable urban/rural pair for any urbanization impact study. Generally, the UHI at a regional scale is best measured using remotely sensed data with one or multiple thermal bands. When explaining the character of remotely sensed UHI, Roth et al. [108] assert, "satellite-derived surface heat islands are in a separate class and it is not clear that they will match others measured by more conventional means in the urban canopy layer or the urban boundary layer". Their precautionary statement relates in part to the surface "seen" by remote sensing platforms that depend on altitude and the camera or sensor angle. Imagery collected at nadir and/or high altitude primarily consists of rooftops, streets, crop fields, and vegetation canopies. Observations from lower heights at oblique angles consist of items seen from a bird's-eye perspective plus varying degrees of vertical features in the landscape, such as the walls of buildings. As a result, angle can have a large influence on the urban surface temperatures recorded by airborne and spaceborne thermal infrared sensors [109]. Another concern regards mixed pixels (i.e., individual pixels containing surfaces having different physical properties, depending on the spatial resolution of the data), which can complicate image analysis. This is especially true for thermal sensors aboard satellites, because most have a spatial resolution that is coarser than the other spectral bands on the satellite. The typical variation of urban surface properties also complicates thermal sensors. A final consideration when using remotely sensed imagery involves correcting for atmospheric attenuation. For many applications, these issues are far outweighed by remote sensing's benefits. With high spatial resolution thermal data, these issues can typically be resolved. Additionally, from a macro research perspective, remotely sensed thermal data have the major advantage of investigating UHI and UHIRIP at a broad scale, permitting focus on environmental issues in urban agglomerations and surrounding areas, and at urban and non-urban interfaces.
Land 2021, 10, x FOR PEER REVIEW 6 of 31 when using remotely sensed imagery involves correcting for atmospheric attenuation. For many applications, these issues are far outweighed by remote sensing's benefits. With high spatial resolution thermal data, these issues can typically be resolved. Additionally, from a macro research perspective, remotely sensed thermal data have the major advantage of investigating UHI and UHIRIP at a broad scale, permitting focus on environmental issues in urban agglomerations and surrounding areas, and at urban and non-urban interfaces.

Simulation and Projection of UHI and UHIRIP
Applying theories of landscape ecology [104], UHI studies focus on moving from static spatial structures of urban thermal patterns to the change dynamic of spatial patterns and processes of urban thermal characteristics. The spatial structure of UHI patterns determines the processes of UHI impacts. Li et al. [110] simulated the urban climate of various generated cities under the same weather conditions. By studying various city shapes, they generalized and proposed a reduced form to estimate UHI intensities based only on the structure of urban sites, as well as their relative distances. They concluded that in addition to the size, the UHI intensity of a city is directly related to the density and the amplifying effect that urban sites have on each other. Their approach can serve as a UHI rule of thumb for the comparison of urban development scenarios. Ramírez-Aguilar and Lucas Souza [111] present a study based on the relationship between UHI and population size (p) by considering the population density (PD) and the urban form parameters of different neighborhoods in the city of Bogotá, Colombia. They concluded that urban form, expressed by land cover and urban morphology changes caused by population density, has a great effect on temperature differences within a city. Advances in computing technology have fostered the development of new and powerful deep learning techniques that have demonstrated promising results in a wide range of applications. In particular, deep learning methods have been successfully used to classify remotely sensed data collected by Earth observation instruments [112]. Deep learning algorithms, which learn the representative and discriminative features in a hierarchical manner from the data, have recently become a hotspot in the machine-learning area, and have been introduced into

Simulation and Projection of UHI and UHIRIP
Applying theories of landscape ecology [104], UHI studies focus on moving from static spatial structures of urban thermal patterns to the change dynamic of spatial patterns and processes of urban thermal characteristics. The spatial structure of UHI patterns determines the processes of UHI impacts. Li et al. [110] simulated the urban climate of various generated cities under the same weather conditions. By studying various city shapes, they generalized and proposed a reduced form to estimate UHI intensities based only on the structure of urban sites, as well as their relative distances. They concluded that in addition to the size, the UHI intensity of a city is directly related to the density and the amplifying effect that urban sites have on each other. Their approach can serve as a UHI rule of thumb for the comparison of urban development scenarios. Ramírez-Aguilar and Lucas Souza [111] present a study based on the relationship between UHI and population size (p) by considering the population density (PD) and the urban form parameters of different neighborhoods in the city of Bogotá, Colombia. They concluded that urban form, expressed by land cover and urban morphology changes caused by population density, has a great effect on temperature differences within a city. Advances in computing technology have fostered the development of new and powerful deep learning techniques that have demonstrated promising results in a wide range of applications. In particular, deep learning methods have been successfully used to classify remotely sensed data collected by Earth observation instruments [112]. Deep learning algorithms, which learn the representative and discriminative features in a hierarchical manner from the data, have recently become a hotspot in the machine-learning area, and have been introduced into the geoscience and remote sensing community for remotely sensed big data analysis [113]. With climate change, the simulation and projection of UHI and its regional impact by using computer technology (deep learning) and remotely sensed data are becoming more important for urban planning and policy makers.

Challenges for Land Surface Temperature and Emissivity Retrieval (Separation)
Land surface temperature and emissivity are two important surface parameters that can be derived from remotely sensed data after atmospheric correction [114][115][116]. Besides radiometric calibration and cloud detection, two main problems need to be resolved in order to obtain land surface temperature and emissivity values from various satellite sensors. These problems are often referred to as land surface temperature and emissivity separation (TES) from radiance at ground level, and as atmospheric corrections in the literature [117,118]. Reliable retrieval of urban and intra-urban thermal characteristics using satellite thermal data depends on accurate removal of the effects of atmospheric attenuations, as well as angular and land surface emissivity. In the thermal infrared of remotely sensed data, the emission of the targets is dominant when compared with the reflection, and this radiation is a function of two unknowns-the emissivity and the temperature of the target [119]. The temperature and emissivity separation is complex because of the existing non-linear relationship between temperature and radiance. The complex dynamics of these relationships within the target (atomic level) propagates in a cascade effect, reflecting variations in determining emissivity. Mohamed et al. [117] reviewed details of LST and land surface emissivity (LSE) retrieval methods and their potential for adoption in medium spatial resolution, including ASTER and Landsat. The review further comments on spatial and temporal prospects of effective intra-urban surface thermal mapping. They also suggested future development of land surface temperature and emissivity estimation for UHI assessment. Li et al. [120] described the theoretical basis of LSE measurements and reviewed the published methods. They also categorized these methods into (1) (semi-) empirical or theoretical methods, (2) multi-channel temperature emissivity separation (TES) methods, and (3) physically based methods (PBMs). Then, they discussed the validation methods that are important for verifying the uncertainty and accuracy of retrieved emissivity. Finally, the prospects for further developments are given. These studies provided a forum for assessing what had been achieved by the UHI community over four decades, and what needs to be done in the near future. It is clear that the observation, experiments, and algorithm development efforts, although completely worthwhile for scientists, need to deliver various datasets, especially from remotely sensed data to modelers working in the areas of UHI and UHIRIP at local, regional, and global levels. A lot of basic theoretical research and scientific verification work has been done on scale issues, as well as scaling issues including emissivity and temperature measurements related to remote sensing standards [121]. All of the methods described in Rolim et al. [119] represent the largest and main part of the existing methods of temperature and emissivity separation developed in the last four decades, but further research is necessary for more precise methods that are less susceptible to errors during the separation of these variables.

The Relationship between Atmospheric Heat Islands and Surface UHI through Combining Coincident Remote Sensing and Ground-Based Observations
Generally, UHI data are obtained from one of two sources-weather stations and remote sensing. Remotely sensed data have been used to observe how UHI impacts climate change in urban and non-urban areas for decades because of the multiple temporal and spatial resolutions of remotely sensed datasets. Hundreds of published studies explore UHI and its impacts by using these two data sources, but the relationship between air temperature obtained from field stations and surface temperature derived from remote sensing remains unclear. Wang et al. [122] investigated the relationship of canopy UHI (CUHI) and surface UHI (SUHI) using four observations per day, without temporal averaging, in four different cities in two different global regions, with 201 of 2232 CUHI-SUHI pairs exhibiting significant UHI differences in their spatial distributions and intensities. The results indicate that 81.09% of the UHI differences occurred during the daytime and were caused by local air advection related to wind speed ≥2 m/s and land surface conditions in the study areas. They concluded that a joint analysis of CUHIs and SUHIs should be conducted to characterize urban thermal environments, and that current urban planning procedures should integrate these UHI differences to develop effective mitigation strategies and adaptation measures. The combination of both types of UHI sub-components provides added value for quantifying urban thermal environments, which can assist in developing effective mitigation strategies and adaptation measures. A growing trend is to combine the two methods, both with their own advantages [59].

Develop Controlling Approaches for UHI and UHIRIP
UHIs occur when cities replace other land covers with dense concentrations of pavement, buildings, and other surfaces that absorb and retain heat. This effect increases energy costs (e.g., for air conditioning), air pollution levels, and heat-related illness and mortality. UHI results from increases in built-up surfaces in urban areas, whereas increasing vegetation cover and water surfaces within cities or urban agglomerations could improve the urban ecological function and thereby improve urban environments for humans [123]. The importance of optimizing urban LULC planning and the development of UHI mitigation methods is increasing. Progress has been made to this end [67,124,125], with the development of UHI mitigating technologies [126]. Ulpiani et al. [127] reviewed an infrared emissivity dynamic switch against overcooling, which is aimed at collecting state-of-theart technologies and techniques to dynamically control the heat transfer to and from the radiative emitter and to ultimately modulate its cooling capacity using spacecraft thermal control, thermal camouflage, and electronics. This work discussed prominent pathways toward technically and economically effective integration in the built environment for UHI and UHIRIP.

Impacts on Human Health
Climate change, increasing urbanization, and an aging population in much of the world are likely to increase the risks to health from UHI, particularly from heat exposure. Additionally, increased urbanization has resulted in a more extensive UHI effects, causing more frequent and intense heat waves in urban regions compared with rural locales [67,128,129]. In urban and surrounding areas, heat waves will be exacerbated by the UHI effect and will have the potential to negatively influence the health and welfare of residents. Heaviside et al. [130] suggest that UHI contributed around 50% of the total heat-related mortality during the 2003 heat wave in the West Midlands of the UK. Moon [131] concluded that the mortality and morbidity risks of diabetic patients under the heat wave were mildly increased by about 18% for mortality and 10% for overall morbidity. Li et al. [132] found that high temperature significantly increases the risk of mortality in the population of Jinan, China. Most research in this topic uses both weather station and in situ measurements in order to investigate the health effects of UHI [129]. Some results [133] show that different sites (city center or surroundings) have experienced different degrees of warming as a result of increasing urbanization [131]. Johnson et al. [134] suggest that thermal remote sensing data can be utilized to improve the understanding of intra-urban variations of risk from extreme heat. The refinement of the current risk assessment systems could increase the likelihood of survival during extreme heat events and assist emergency personnel in the delivery of vital resources during such disasters. The conclusion is that UHI is directly linked to adverse health effects from exposure to extreme thermal conditions.

UHI and UHIRIP on LULC Differences and Change
UHI is a result of continued urbanization, urban agglomeration, and associated increases in paved areas and buildings. Mitigation strategies have been developed to increase vegetation and water surface areas within urban areas to reduce the magnitude of the temperature. One measure of UHI's ecological footprint is estimated by calculating the increase of the cooling demand caused by the heat island over the urban area, and then translating the increased energy use to environmental cost [123,125,135]. Some research shows that the UHI effect has become more prominent in areas of rapid urbanization and in urban agglomerations [136,137]. The spatial distribution of UHI has changed from a mixed pattern, where bare land, semi-bare land, and land under development were warmer than other LULC types, to extensive UHI, as contiguous urbanized blocks grew larger [38,138]. Some analyses showed that the higher temperature in the UHI had a scattered pattern and was related to certain LULC types [97]. In order to analyze the relationship between UHI and LULC changes, some studies attempted to employ a quantitative approach for exploring the relationship between surface temperature and several indices, including the Normalized Difference Vegetation Index (NDVI), Normalized Difference Water Index (NDWI), Normalized Difference Bareness Index (NDBaI), and Normalized Difference Build-up Index (NDBI). It was found that correlations between NDVI, NDWI, NDBI, and temperature are negative when NDVI and NDWI are limited in range, but there is a positive correlation between NDBI and temperature [139][140][141][142].

Impacts on Regional Economics
Because UHI is related to a significant increase in surface temperature and changes in precipitation patterns, it can potentially affect local economies and the social systems of cities [143]. Some studies [144,145] show that the critical sectors of services, agriculture, and tourism may be strongly affected by future UHIs. To counterbalance the consequences of the increased urban surface temperatures, important research has been carried out resulting in the development of efficient mitigation technologies. In particular, some studies [102,146] have documented the development of highly reflective materials, cool and green roofs, cool pavements, urban greens, water surface, and other mitigation technologies. UHIRIP includes economic impacts, such as increases of energy consumption for cooling purposes, as well as an increase in the peak electricity load, which is a factor for planning maximum power source capacities [147]. Scientists from Australia reported that the total economic cost to the community due to hot weather is estimated to be approximately $1.8 billion in present value terms. Approximately one-third of these impacts are due to heat waves. Of the total heat impact, the UHI effect contributes approximately $300 million (AUD) in present value terms for the city of Melbourne, Australia [9]. Estrada, Botzen and Tol [144] provided a quantitative assessment of the economic costs of the joint impacts of local and global climate change for all main cities around the world. They estimated the UHI effect for the 1692 largest cities in the world for the period 1950-2015, and predicted that the percentage of city gross domestic product (GDP) that would be lost for the median city in 2050 due to global climate change alone would be relatively small: 0.9% and 0.7% for the RCP8.5 and RCP4.5 emission scenarios, respectively [144]. At the end of the century, these impacts will increase to 3.9% and 1.2%, respectively. Cost-benefit analyses are presented of urban heat island mitigation options, including green and cool roofs and cool pavements. It has been shown that local actions can be climate risk-reduction instruments. Furthermore, limiting the urban heat island through city adaptation plans can substantially amplify the benefits of international mitigation efforts.

Impact on Biodiversity
Besides UHI, urban development causes wildlife habitat loss and fragmentation, threatens wildlife populations, increases fire risk, and reduces biodiversity [2,148]. These problems are of particular concern in the wildland urban interface (WUI), where homes and associated structures are built among forests, shrublands, or grasslands [1,148,149]. The WUI has received considerable attention because of recent increases in both the number of structures destroyed and the area burned annually by wildland fire.Čeplová et al. [150] studied three habitats with different disturbance regimes in 45 central European settlements of three different sizes. Their results highlight the importance of urban size as a factor shaping the biodiversity of native and alien plant communities in individual urban habitats, Land 2021, 10, 867 10 of 30 and the important role of habitat mosaic for maintaining high species richness in city floras. The study of Coluzzi et al. [151] represented a first step to improve the description of relevant processes to protect natural habitats and quality agriculture, therefore combating land degradation and detrimental climate change effects. Kaiser et al. [152] monitored temperature and relative air humidity in wooded sites characterized by different levels of urbanization in the surroundings, and investigated the effect of urbanization at the local and landscape scale on two key traits of biological fitness in two closely related butterfly species that differ in thermal sensitivity.

Remotely Sensed Thermal Datasets
Remote sensing derived LST is effective for UHI and UHIRIP studies. Satellites can quickly obtain continuous information over a large geographic area that can be maintained in long-term archives. LST for large geographic areas can be derived from surface radiation of heat measured by satellite sensors. This is particularly attractive when investigating the surface UHI in multiple cities or urban agglomerations at various spatial extents. Along with the extensive spatial coverage, many satellites record multiple wavelengths of electromagnetic energy that can be used to decipher a wealth of information, in addition to thermal information ( Figure 2). Consequently, multispectral imagery allows for a comparative analysis between LST and other variables, such as land cover and vegetation indexes [50,153], specifically the interaction between UHI and LULCC [154]. Remote sensing can also be used to track the patterns of change in UHI over time through various time periods from a day, to years, and even a time series of decades [38,[155][156][157][158]. Because information is desired at a high spatial resolution and dense temporal frequency, data from multiple sensors can make more accurate and reliable quantitative assessments of UHIRIP studies [60]. Table 3 includes a list of the main remotely sensed datasets that have been recently used to derive LST and analyze UHI and UHIRIP. Voogt and Oke [6], and others [156,178] pointed out that improved spatial and spectral resolution of sensors and advances in digital image processing techniques increase the usefulness of remote sensing for UHI and UHIRIP studies. Forster [179] also stated that satellite, radar, and airborne sensors can provide spatially continuous information  The rapid development of remote sensing technology offers more potential for accurate and reliable quantitative assessments of UHI (Table 3 and Figure 3). Many researchers ( Table 2) have used remotely sensed LST to assess UHI over various geographic areas. However, for all of these studies, the 1 × 1 km spatial resolution of coarse datasets was found to be suitable only for broad-scale urban temperature mapping ( Table 4). The higher resolution of Landsat time series is suitable for UHIRIP at various scales (Table 4).
Land 2021, 10, x FOR PEER REVIEW 12 of 31 ECOSTRESS can accurately monitor UHI trends over time in vulnerable areas such as the urban and non-urban interface. With more available remotely sensed data (Figure 3), innovative studies like these hint at the potential for remote sensing to play an even more prominent role in research of urban climate, urban environment, urban ecological service, and urban planning in the future.

Algorithms for UHI and UHIRIP in Urban and Non-Urban Interface Studies Based on Remotely Sensed Data
Generally, the methods for evaluating UHI and UHIRIP can be summarized into four basic types: (1) historical weather station data, (2) field observation, (3) computer simulation, and (4) remote sensing technology. The limitations of the first three methods have been well documented [53,57,105,180]. In this paper, we only focus on the methods that use remote sensing technology. A number of algorithms (or methods) have been developed to estimate UHI and UHIRIP from remotely sensed data (Table 5), including simple empirical approaches to complex methods based on remotely sensed data assimilation using various models. The structure of the UHIRIP pattern centroid in three dimensions indicates the overall variation of the intensity and distribution of the UHI in space and time. The simplified relationship of thermal data and UHI has been applied from a local spatial scale using airborne very high-resolution images to a broad scale with AVHRR, MODIS, ASTER, and Landsat data at regional and continental levels. Assimilation procedures of UHI often require remotely sensed data over different spectral domains to retrieve input parameters that characterize surface properties such as thermal properties, albedo, NDVI, and other indices. A brief review of these approaches is presented in Table 5, with a discussion about the main physical bases and assumptions of various models.
Detailed knowledge of UHI and UHIRIP, especially latent and sensible heat flux components, is important for monitoring the climate change of the land surface. The main methods classically used to measure UHI are appropriate to field observations [24][25][26]189], but do not allow for an estimation of UHI at large spatial scales. For operational  Voogt and Oke [6], and others [156,178] pointed out that improved spatial and spectral resolution of sensors and advances in digital image processing techniques increase the usefulness of remote sensing for UHI and UHIRIP studies. Forster [179] also stated that satellite, radar, and airborne sensors can provide spatially continuous information pertaining to numerous variables in urban environments that complement field observations. An increasing number of studies directly relate remotely sensed data to in situ field data [180,181], and applications of remote sensing technology will expand UHI studies to various geographic extents. An exciting recent trend in UHI and UHIRIP research involves coupling remotely sensed data with ancillary and social economic datasets from multiple sources ( Table 2). The typical examples include (1) fractional vegetation cover derived from satellite data to improve model simulations of UHI [182], (2) incorporation of remotely sensed data into a model that partitioned various fluxes in the surface energy balance [183], (3) integrating high-resolution multispectral data with property tax records to investigate the contribution of residential land use to UHI formation [184], (4) studying the potential application of change in urban green space as an indicator of urban environmental quality change [185], (5) using parameters from thermal satellite data and three-dimensional virtual reality models to better understand the factors controlling urban environmental quality (UEQ) [186], (6) further advancing the use of remotely sensed imagery to evaluate UEQ by estimating ground-level particulate matter (PM) concentrations using satellite-based data [187], and (7) estimating the value of U.S. urban tree cover for reducing heat-related health impacts and electricity consumption [188]. In addition, NASA's Ecosystem Spaceborne Thermal Radiometer Experiment on the International Space Station (ECOSTRESS) was launched in June 2018, and is able to image fine-scale temperatures in cities at a 70 × 70 m resolution throughout different times of the day, every 3-5 days on average, over most of the globe [146]. With new algorithm development, ECOSTRESS can accurately monitor UHI trends over time in vulnerable areas such as the urban and non-urban interface. With more available remotely sensed data (Figure 3), innovative studies like these hint at the potential for remote sensing to play an even more prominent role in research of urban climate, urban environment, urban ecological service, and urban planning in the future.

Algorithms for UHI and UHIRIP in Urban and Non-Urban Interface Studies Based on Remotely Sensed Data
Generally, the methods for evaluating UHI and UHIRIP can be summarized into four basic types: (1) historical weather station data, (2) field observation, (3) computer simulation, and (4) remote sensing technology. The limitations of the first three methods have been well documented [53,57,105,180]. In this paper, we only focus on the methods that use remote sensing technology. A number of algorithms (or methods) have been developed to estimate UHI and UHIRIP from remotely sensed data (Table 5), including simple empirical approaches to complex methods based on remotely sensed data assimilation using various models. The structure of the UHIRIP pattern centroid in three dimensions indicates the overall variation of the intensity and distribution of the UHI in space and time. The simplified relationship of thermal data and UHI has been applied from a local spatial scale using airborne very high-resolution images to a broad scale with AVHRR, MODIS, ASTER, and Landsat data at regional and continental levels. Assimilation procedures of UHI often require remotely sensed data over different spectral domains to retrieve input parameters that characterize surface properties such as thermal properties, albedo, NDVI, and other indices. A brief review of these approaches is presented in Table 5, with a discussion about the main physical bases and assumptions of various models.
Detailed knowledge of UHI and UHIRIP, especially latent and sensible heat flux components, is important for monitoring the climate change of the land surface. The main methods classically used to measure UHI are appropriate to field observations [24][25][26]189], but do not allow for an estimation of UHI at large spatial scales. For operational applications to ecological conservation and city planning, managers and engineers need accurate estimates of land surface temperature and UHI at broad spatial scales. New algorithms based on remotely sensed data have been developed to use the imagery of various spatial resolutions and temporal frequency to evaluate UHI [190][191][192]. It is often difficult to classify these methods because their complexity depends on the balance between the empiricaland physical-based modules used. Nevertheless, we summarize some algorithm (model) categories in the following subsection. Table 5. Methods used to measure UHI and UHIRIP using remotely sensed data.

Method Sensor Period Example
Calculate LST All thermal bands 1970s-current Avdan and Jovanovska [193], and Peng et al. [194]  Compare multi-temporal LST images The normalization of the temperature based on the mean and standard deviation in high and low temperature areas.
Streutker [39] Common normalization of temperture based on min and max LST of the same image in the same way as for NDVI. A normalized ratio scale technique.
Chen et al. [38] Statistical analyses of UHI The relationship between LST, NDVI, ground vegetation (GV), and impervious surface area (ISA). Multiple linear regression. Geographically weighted regression.

LST and UHI Intensity Calculation
LST calculation, including empirical direct methods where remotely sensed data are introduced directly in semi-empirical models to estimate LST, is the simplified relationship between thermal infrared remotely sensed and meteorological data [14]. This method allows for the characterization of UHI intensity both at the local scale, using ground measurements, and over large areas, using satellite data, by calculating a cumulative temperature difference [55,92]. Most current operational models [60] use remote sensing directly to estimate the input parameters and LST.
Seasonal information captures the annual profile of LST and its trend over long time periods, and is essential to the study of UHI [207]. Therefore, remote sensing has been used to accurately monitor and compare the LST difference in the same season in different years and trends over long time periods. In the last 10-15 years, thermal sensor technology has been rapidly developing (Figure 3). Three types of methods have been developed to estimate LST with remotely sensed data: the single infrared channel method; the split window method; and a new day-night MODIS LST method, which is designed to take advantage of the unique capability of the MODIS instrument [55]. Recently, Peng et al. [194] proposed a wavelet coherence approach to exploring spatial heterogeneity and the scale-dependence of the relationship between LST and multiple influencing factors. The advantages, disadvantages, and applicability of these three types of algorithms are summarized in Table 6.

Comparing the Difference between Core Urban and Non-Urban Area
Many studies have documented the use of LST data to observe meso-scale temperature differences between urban and rural areas in cities worldwide [156,208,209]. The land surface temperature (LST) of core urban areas is generally higher than the surrounding rural areas, and has a strong correlation with land cover [153]. UHII analysis is the most common method to compute the magnitude and extent of UHI by evaluating the LST difference between urban and surrounding non-urban areas [162]. These analyses are often supported with auxiliary land surface information, such as land cover and impervious surface area (ISA). Deterministic models generally are generally based on more complex models that compute the intensity of UHIRIP in space and time. Remotely sensed data are used at different modeling levels, either as the input parameters to characterize the different surface covers, or in assimilation procedures, which aim to retrieve adequate parameters for the LST computation. Some examples of these studies are shown in Table 5. UHI intensity was typically quantified in two steps in these studies [60]. First, urban and non-urban areas were defined and delineated from land cover or ISA maps. Urban areas are usually defined as land with a relatively higher proportion of ISA [38,95], whereas nonurban areas have various definitions in different studies, but generally include non-urban land cover classes. Different sized rural and suburban zones have been used as reference areas. Other land covers, such as water bodies, cropland, forest, and low-intensity ISA, have also been used as references in the studies [101]. Second, the area-weighted mean urban-reference LST differences were calculated to reflect the UHI intensity [69,210] or magnitude. Some studies identified "hotspots" based on positive UHI intensity in certain time periods [211,212]. A positive value of UHI intensity indicated an urban heating effect, while a negative value represented a cooling effect. A few studies also quantified the UHI intensity using small numbers of representative pixels in urban and reference areas instead of the area-weighted mean value for the purpose of surface-air UHI comparison [122][123][124] or UHI attribution analysis [125,126]. The urban-reference difference method facilitates a comparative analysis of UHIs among cities and urban agglomerations, regions, and across the globe, but the validity of such comparisons can be limited by the large uncertainties associated with urban and reference definitions [68]. Recent research [97] performed a comprehensive and consistent analysis of surface UHI and UHIRIP using Landsat LST ARD time series and dynamic land cover datasets in the Sioux Falls, SD, area. It shows that the use of time series of LST and land change dynamic data provided a consistent and quantitative analysis for the distribution and change of UHI intensity and UHIRIP (Figure 4). We further discuss limitations in Section 5. magnitude. Some studies identified "hotspots" based on positive UHI intensity in certain time periods [211,212]. A positive value of UHI intensity indicated an urban heating effect, while a negative value represented a cooling effect. A few studies also quantified the UHI intensity using small numbers of representative pixels in urban and reference areas instead of the area-weighted mean value for the purpose of surface-air UHI comparison [122][123][124] or UHI attribution analysis [125,126]. The urban-reference difference method facilitates a comparative analysis of UHIs among cities and urban agglomerations, regions, and across the globe, but the validity of such comparisons can be limited by the large uncertainties associated with urban and reference definitions [68]. Recent research [97] performed a comprehensive and consistent analysis of surface UHI and UHIRIP using Landsat LST ARD time series and dynamic land cover datasets in the Sioux Falls, SD, area. It shows that the use of time series of LST and land change dynamic data provided a consistent and quantitative analysis for the distribution and change of UHI intensity and UHIRIP (Figure 4). We further discuss limitations in Section 5.

UHI and UHIRIP Analysis by Using Urban Ecological Indices
Many studies have compared UHIRIP to ecological indices [103,149,[213][214][215], vegetation fraction, and percent ISA, finding strong correlations with mean LST. Landscape metrics indicate that urban landscape configuration also influences the surface UHIRIP [216]. The latest vegetation index methods and inference methods use remote sensing to compute a reduction factor (such as Kc or Priestley Taylor-alpha parameters) for the estimation of the actual UHI [203]. Different papers deal with these approaches in the various journals, and these approaches use land cover [217,218], LST pattern [219,220], and a combination of land cover and LST pattern [221][222][223] as monitoring indicators of UHI.
Urban ecological status is closely related to the quality of human life and the development of urban economies. A timely and objective understanding of urban ecological status, particularly in urban and non-urban interface areas [224], has become an increasing important. Scientists have been developing a remote sensing-based ecological index for the measure of urban ecological status under UHI [213,225]. This urban ecological status index (UESI) aims to integrate four important ecological indicators that are frequently used in evaluating urban ecology. The four indicators include greenness, wetness, dryness, and heat, and can be represented by four remote sensing indices or components: NDVI, normalized difference built-up and soil index (NDBSI), wetness component of the tasseled cap transformation (Wet), and LST, respectively. Instead of a simple or weighted addition of the four indicators, a principal component analysis (PCA) can be utilized to compress the four indicators into one index in order to assess the overall urban ecological status under UHI. The calculation of the UESI can be fully automated, avoiding the need to assign threshold values or weights during the computing procedure. Therefore, the UESI can be used to easily and objectively assess urban ecological status. Combined with change detection, UESI can also be used to monitor the change of the ecological status of the core urban and surrounding non-urban areas between different years. In practice, the index was successfully applied in a multitemporal ecological status assessment [34]. Pan [221] used the G index spatial aggregation analysis to calculate the urban heat island ratio index, and the landscape metrics to quantify the changes of the spatial pattern of the UHI from the aspects of quantity, shape, and structure. Pan found that the heat island strength had a negative linear correlation with urban vegetation coverage, and a positive logarithmic correlation with urban impervious surface coverage. Bala et al. [226] developed the Urban Heat Intensity Ratio Index (UHIRI) to quantify urban heat intensity. This work analyzed the variation in LST with land cover changes in Varanasi, India, from 1989 to 2018, using Landsat images, and concluded that the replacement of vegetation with urban land cover has a severe impact on increasing UHI intensity.

Various Statistical Models
Statistical models and machine learning have also been proposed to measure UHI [227]. Among these studies, a Gaussian surface model has been utilized the most because it can provide not only the intensity, but also the spatial extent and the central location of the UHI. The kernel convolution method has also been proposed to study UHI effects because of its high efficiency in characterizing the temperature values over space in a continuous surface [227]. Chun and Guldmann [77] explored the urban determinants of UHI using two-dimensional (2D) and three-dimensional (3D) urban information as the input for spatial statistical models. The results show that solar radiations, open spaces, vegetation, building roof-top areas, and water strongly impact surface temperatures, and that spatial regressions are necessary in order to capture the neighboring effects. Recently, Li et al. [81] estimated UHI intensity by linear regression functions between LST and regionalized ISA. These statistical models could avoid the bias caused by the definitions of urban−rural areas or the choice of the representative pixels, and thus facilitate the comparison of UHI among cities. Szymanowski and Kryza [228] addressed the issue of the potential usefulness of remotely sensed data and their derivatives for UHI modeling. Statistically significant models explained 71% to 85% of the air temperature variance. It has been stressed that remotely sensed data are important sources to model urban air temperature heat islands. However, in all of these studies, such models worked less effectively in cities frequently covered by clouds, in arid landscapes, and in urban agglomerations, so they have only been applied in a few UHI studies to date. Recently, Lai et al. [79] published the statistical estimation of next-day nighttime surface urban heat islands of selected cities. Most previous studies modelled the SUHI variations for the past period, but rarely investigated the estimation for future UHIs, especially at the daily (i.e., day-to-day) scale. To address this issue, this study incorporated both meteorological and surface controls to estimate next-day nighttime UHIs using a support vector machine regression (SVR) model. Some uncertainties exist in terms of the Gaussian modelling and UHI estimators, which may limit estimation accuracy. Nevertheless, by providing a feasible yet simple approach for estimating next-day nighttime UHIs, this study fills a knowledge gap in the UHI estimation and is helpful for supporting adaptation to and mitigation of UHI and UHIRIP.

Spatial−Temporal Time-Series Algorithm
Advances in computing technology have fostered the development of new and powerful data fusion, gap filling, machine learning, and deep learning techniques that have demonstrated promising results in a wide range of applications [229]. Models to fuse data from multiple sensors and fill gaps can improve UHI monitoring using an ensemble of dense time series of thermal data with a high spatial resolution [199,230]. Zhou et al. [200] presented a new algorithm that focuses on data gap filling using clear observations from orbit overlap regions to obtain Landsat LST data. Multiple linear regression models were established for each pixel time series to estimate the stable predictions and uncertainties. Liu et al. [83] comprehensively quantified the spatial-temporal patterns of surface urban heat island by investigating the relationship between LST and the land cover types, and the associated landscape components. Such approaches have been used to generate temporally dense and high-resolution LST over long time periods by integrating the observations of Landsat, MODIS, AVHRR, VIIRS, and ECOSTRESS [191,231]. These datasets facilitate subtle analyses of monthly, seasonal, and yearly trends in UHI intensity at regional levels. Machine learning (ML) has become popular in UHI and UHIRIP, but its use has remained restricted to predicting, rather than understanding, the natural world. ML techniques may not be the solution to all the problems remotely sensed data might have. However, these techniques provide a powerful set of tools that deserve serious attention to deal with some relevant UHI and UHIRIP remotely sensed data problems [232]. Lucas [233] points out that ML differs from the broader field of statistics in two respects: (1) the estimation of parameters that relate to the real world is less emphasized, and (2) the driver of the predictions is expected to be the data rather than expert opinion and careful selection of plausible mechanistic models. The Google Earth engine (GEE) is a cloud-based platform for planetary-scale geospatial analysis that brings Google's massive computational capabilities to bear on a variety of high-impact societal and environmental issues [234]. GEE has many functions that could be used to analyze UHI and UHIRIP at local, regional, and global levels ( Figure 5). Some research has generated consistent large-scale UHI and UHIRIP analysis based on optimal data and ML algorithm selection using GEE [235,236]. The advanced GEE cloud-based platform and the large number of geosciences and remote sensing datasets archived in GEE were used to analyze land the cover dynamics (236), and the results showed the advantages of using GEE to analyze the spatiotemporal dynamics of the LULCC, vegetation cover, LST, and climate for a long time series, and highlighted the importance of environmental protection. The power here lies in the way a scientist defines their questions and uses machine learning alongside other methods. Techniques for data analysis and interpretation that fully incorporate the temporal dimension remain an area of intense research and represent an important challenge for operational UHI and UHIRIP monitoring.

Summary of UHI and UHIRIP Based on Remotely Sensed Data
This review provides an overview of research on UHI and UHIRIP based on remote sensing techniques, sensors, and algorithms, as listed in Tables 3-6, respectively. Much work has been completed on UHIRIP in the last four decades, and we have endeavored to keep updated with new methods and results. A significant research limitation still exists: the quantification of UHI and its regional impacts using high-resolution time-series

Summary of UHI and UHIRIP Based on Remotely Sensed Data
This review provides an overview of research on UHI and UHIRIP based on remote sensing techniques, sensors, and algorithms, as listed in Tables 3-6, respectively. Much work has been completed on UHIRIP in the last four decades, and we have endeavored to keep updated with new methods and results. A significant research limitation still exists: the quantification of UHI and its regional impacts using high-resolution timeseries remotely sensed thermal data in the urban and non-urban interface. Some of the algorithms listed in Table 5 may be the most practical approaches to assess UHI in core urban areas of cities and surrounding areas, but characterization of UHI across broad areas is necessary in order to inform monitoring, reporting, science, and policy. Being able to relate LST to UHI is especially important when such datasets are being used to inform policy decisions or to communicate outside of the scientific community. The increasing availability of remotely sensed data across a range of spatial resolutions and temporal frequencies, and technological improvements in image processing capacity and storage, have led to advances in the methods used to monitor UHI more frequently and accurately.
Assessing the uncertainty and accuracy of UHI data is important. A sensitivity analysis not only provides a framework for assessing the potential for bias and the extent of uncertainty in UHI estimates, but also reveals significant factors that determine the extent of UHIRIP in the urban and non-urban interface. Oleson et al. [237] developed an approach to evaluate the robustness of models used to simulate urban heat islands in different environments. The findings indicated that heat storage and sensible heat flux are most sensitive to uncertainties in the input parameters within the atmospheric and surface conditions considered. Sensitivity studies indicate that it is important to not only to accurately characterizing the structure of the urban area, but also to ensuring that the input data reflect the thermal admittance properties of each of the city surfaces.
Currently, a wide variety of methods are employed to characterize UHI for major cities worldwide (Table 5), although most of the applications cited were limited to small areas because of data availability and constraints of storage and computing resources. With the development of gap filling and data fusion models [238], advances in high-performance computing (HPC), and cheaper storage, applications based on high-resolution time series at larger or even regional scales will become the mainstream in the near future [199,231]. While much of the methodological variation described here will persist, future methods will evolve and adapt to greater data volumes and processing capabilities [239]. Legacy change mapping methods that rely on analyst interactions with individual scenes should decline over time given the improved ability to process and characterize time series of rich high-resolution thermal data. However, such spatial−temporal methods that are based on gap filling and data fusion should match the institutional requirements for accuracy. Near-term research objectives will require robust validation datasets in establishing which data-intensive methods are the most appropriate for quantifying UHI over large areas. Techniques for LST data analysis and interpretation that fully incorporate the temporal dimension still require intense research and represent an important challenge for operational UHI research in order to meet management needs.
Technological advances that include machining leaning and artificial intelligence in UHI and UHIRIP using remotely sensed data have led to an explosion of UHI and UHIRIP profiling data from large numbers of multiple data sources [233,240]. This rapid increase in the remotely sensed data dimension and acquisition rate is challenging conventional analysis strategies. Modern machine learning methods, such as deep learning, promise to leverage very large datasets for finding a hidden structure within them, and for making accurate predictions [232]. Deep learning methods are a powerful complement to classical machine learning tools and other analysis strategies, and have been used in a number of applications in UHII and remotely sensed image analyses [241]. The explainable artificial intelligence in UHII and UHIRIP modeling has become more and more important [242]. Interpretable machine learning methods either target a direct understanding of the model architecture (i.e., model-based interpretability) or interpret the model by analyzing the model behavior (post hoc interpretability) [242].
Currently, most of the time-series algorithms used to map UHIRIP include data from the temporal domain of AVHRR and MODIS, and the spatial domain of the data is almost entirely neglected. Although these datasets with a lower spatial resolution and higher temporal frequency can detect a change of UHI in real time, they often lack pertinent spatial detail. Even though many UHI analysis algorithms have been developed [60], most of the UHI monitoring data derived from the Landsat archive are provided in a time frame that is not near enough to real time to be relevant for specific management needs. With the advances in HPC and cheaper storage, applications based on Landsat time series at continental or even global scales will be the mainstream in the next few years.
To date, information from Landsat time-series thermal data has taken the form of statistical metrics, change metrics, pattern distribution, or trend components used in UHI impact applications [243]. Improvement of existing approaches, as well as the inclusion of novel techniques, often imported and adapted from other disciplines, are important to fully capitalize on the thermal data in order to produce monthly, seasonal, and annual LST results that meet a wide range of UHI and UHIRIP research needs. Landsat-9, which will be launched in September 2021, will continue collecting images of the Earth's surface in visible, near-infrared, and shortwave-infrared bands, as well as the thermal infrared radiation, or heat, of the Earth's surface from two thermal bands. The future European Space Agency's LSTM (Land Surface Temperature Monitoring) or Sentinel 8 mission will carry a high spatial−temporal resolution thermal infrared sensor to provide records of land-surface temperature. Land-surface temperature measurements are key variables to understand and respond to climate variability and natural hazards, such as urban heat island issues. The main objective of LSTM is to deliver global high spatial−temporal dayand night-time land surface temperature measurements. LSTM will operate from a low-Earth, polar orbit, to map both land-surface temperature and rates of evapotranspiration. It will be able to identify the temperatures of individual fields and image the Earth every three days at a 50 m resolution. Another future thermal sensor is Thermal infraRed Imaging Satellite for High-resolution Natural resource Assessment (TRISHNA), which is a future high-resolution space-time mission in the thermal infrared (TIR) led jointly by the French (CNES) and Indian (ISRO) space agencies. One of scientific objectives guiding the definition of the mission is the monitoring of the urban environment. TRISHNA will be positioned on a polar orbit and provide a revisit of three passages over 8 days with global coverage. The time of passage around 13:00 p.m. LST allows thermal data to be collected in the middle of the day, but also in the middle of the night. The instrument will offer four thermal channels (8.6 µm, 9.1 µm, 10.4 µm, and 11.6 µm) and six optical channels (485 nm, 555 nm, 650 nm, 860 nm, 1380 nm, and 1650 nm) with a spatial resolution between 50 m and 60 m for all channels. All of these observations acquired from thermal remote sensing will provide more valuable information for natural resource management, hazard monitoring, and scientific research and applications.

Future Research Directions
Remote sensing technology has been widely applied in the research of UHI and UHIRIP. The most important advantage of using remote sensing thermal data is the wallto-wall coverage of UHI patterns that can meet the needs of spatial and temporal analyses. Remotely sensed data can be used to investigate the surface temperatures of cities and urban agglomerations for various ecosystems with different climate conditions, for example tropical and sub-tropical, temperate and cold temperate, coastal and inland, and arid and semi-arid land at regional scales. These studies are needed to describe surface temperature characteristics in these specific environments and how climate change may be modulating UHI patterns. UHIRIP produces an aggregate impact on weather conditions, land use, human health, biodiversity, ecosystem security, economics, and urban planning [16,244].
Land surface temperature and emissivity retrieval (separation) has always been challenging. Generally, the LSE values needed to apply the method have been estimated from a procedure that uses the visible and near-infrared bands. The algorithm was created using the brightness temperature of the thermal and emissivity of different land cover types, derived from visible and near-infrared bands of various sensors. Compared with field-based observation, remote sensing offers the advantages of a harmonized, long-term, and spatially extensive record to observe LST change. The retrieved LSTs are verified using the near surface temperature of weather station datasets, which will help to improve the accuracy of LST derived from thermal bands. The difference between retrieved LST and Automatic Weather Station (AWS) data indicates that the technique works by giving an error of ±3 • C [245]. These differences can be because of the difference between the resolutions of thermal and visible bands, and a comparison was made between the point measurement (AWS data) 2 m above the surface and surface temperature (retrieved LST). Communicating the results of time-series LST studies that are based on both field weather station observations and remote-sensing time-series data to urban planners, policymakers, and the general public could help inform urban design and decision making.
Using temporally dense time series of remotely sensed data at a high spatial resolution is a growing trend in UHI and UHIRIP research, facilitated by increasing computer capabilities to handle big datasets, machine leaning, deep learning, and Google Earth Engine applications. Landsat ARD, in particular, has great potential to derive LST. Models used to fuse data from across multiple sensors will be developed to increase data temporal density and spatial resolution. Moreover, future sensor improvement on Landsat and aircraft thermal data are possible options. On the other hand, in order to determine the temporal variation of LST using satellite data with restricted overpass times, it appears necessary to use long-time weather station observations to investigate diurnal UHI in various ecosystems, although some new sensors (e.g., ECOSTRESS) can provide this information. Future research is anticipated to improve on methods to simultaneously derive LST and land surface emission (LSE) from hyperspectral TIR, multi spectral-temporal, and TIR-microwave data; additionally, future methods will consider aerosol and cirrus effects [18]. Another viable angle of potential future studies is urban development strategies for mitigating UHI, such as increasing vegetation and water surfaces in urban development.
Climate models are the only tools that account for the complex set of processes that will determine future climate change at both a global and regional level, and assessing regional impacts of climate change begins with the development of climate projections at relevant temporal and spatial scales [246]. The most current existing climate change modeling covers large geographic areas at regional and global levels with relatively low spatial resolutions (>10 km). In the future, LST that is derived from remotely sensed data will support climate change modeling (regional climate models and statistical downscaling models) in UHI and UHIRIP analyses in urban and surrounding areas.
Our analysis indicated that determination is still a central topic of UHI research. Modeling will continue to provide vital and useful results on the spatiotemporal assessment of UHI, especially when models more effectively combine thermal data from multiple sensors. ML (DL) and AI are continuing to grow in popularity in UHI and UHIRIP research. For time series analyses with remote sensing data, a cloud computing platform such as GEE could bring about a substantial change in UHI and UHIRIP analyses, as they have the capability to process big remote sensing datasets and assess the spatiotemporal dynamics of the area quickly. A better integration of remote sensing and station measurements into models is expected. This study also suggests that direct and indirect UHIRIP, especially human health issues, heat wave impacts, air pollution, and ecological security, will receive increasing scientific attention in the future. Research on controlling and adapting to UHI impacts may warrant special attention. The interaction of UHI and UHIRIP, and their changes to LULC based on urban planning, are actively being studied.

Acknowledgments:
We are grateful to the USGS National Land Imaging Program for supporting this research. We would like to thank Carol Deering for assistance in collecting and preparing the literary material. This manuscript was improved thanks to comments and suggestions from Norman Bliss, Matthew Rigge, and Thomas Adamson. Hua Shi and Qiang Zhou's work was performed under USGS contracts G13PC00028 and 140G0119C0001. For Kevin Gallo, the scientific results and conclusions, as well as any views or opinions expressed herein, are those of the author and do not necessarily reflect those of NOAA or the Department of Commerce. Any use of trade, firm, or product names is for descriptive purposes only and does not imply endorsement by the U.S. Government.