Assessing Local Climate Change by Spatiotemporal Seasonal LST and Six Land Indices, and Their Interrelationships with SUHI and Hot–Spot Dynamics: A Case Study of Prayagraj City, India (1987–2018)

: LST has been ﬂuctuating more quickly, resulting in the degradation of the climate and human life on a local–global scale. The main aim of this study is to examine SUHI formation and hotspot identiﬁcation over Prayagraj city of India using seasonal Landsat imageries of 1987–2018. The interrelationship between six land indices (NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI) and LST (using a mono-window algorithm) was investigated by analyzing correlation coefﬁcients and directional proﬁling. NDVI dynamics showed that the forested area observed lower LST by 2.25–4.8 ◦ C than the rest of the city landscape. NDBI dynamics showed that the built-up area kept higher LST by 1.8–3.9 ◦ C than the rest of the city landscape (except sand/bare soils). SUHI was intensiﬁed in the city center to rural/suburban sites by 0.398–4.016 ◦ C in summer and 0.45–2.24 ◦ C in winter. Getis–Ord G i * statistics indicated a remarkable loss of areal coverage of very cold, cold, and cool classes in summer and winter. MODIS night-time LST data showed strong SUHI formation at night in summer and winter. This study is expected to assist in unfolding the composition of the landscape for mitigating thermal anomalies and restoring environmental viability.


Introduction
Globally, 55% of the total populace resided in urban areas in 2018, and prediction statistics show that if this trend continues, then the urban population will account for 68% of the total in 2050 [1].Recently, the IPCC revealed in its report that the global mean surface temperature (GMST) increased by 1.53 • C and the global mean air temperature (GMAT) [both land and ocean] increased by 0.87 • C during the preindustrial period (1850-1900) and recent postindustrial period (2006-2015) [2].The rise in land surface temperature (LST) has severe environmental consequences because unplanned urbanization will deteriorate climate equilibrium and hamper human life and health from microscale to macroscale [3,4].In its AR6 2021 report, IPCC revealed that we are on the way to reaching 1.5 • C more global warming in the next 20 years.The unprecedented changes in the recent past have been highly challenging and alarming, leading to uncertain precipitation, increased glacier melting, mean sea level rising, floods, droughts, damage to agricultural land, and food shortage affecting every region of the globe [5].
The presence of Earth surface objects, such as asphalt, stones, pebbles, and sand, over the city landscape has diverse electromagnetic behavior in terms of evaporation, absorption, and radiation.Longwave radiation, as well as prevailing winds, assimilates into massive heat discharge from the Earth's surface [6].These surface objects weaken evapotranspiration and accelerate sensitivity [7].Consequently, a difference in LST has been observed in city landscapes where core city space experiences higher LST than suburban/rural sites.This distinctive LST characteristic is defined by the surface urban heat island (SUHI).The SUHI has emerged because of the conversion of natural land into built-up space at the cost of water bodies, bare land, and forest [8,9].
LST intensification in urban setups is a perilous factor responsible for deteriorating urban climate and degrading human life and living [10,11].Much attention is now given to mitigating its severity and threats to varied aspects of the environment by policymakers, health authorities, urban planners, climatologists, and environmentalists [12][13][14].In the recent past, researchers have been engrossed in the interrelationship between LST and different land indices, such as the normalized difference built-up index (NDBI), enhanced built-up and bareness index (EBBI), normalized difference vegetation index (NDVI), and normalized difference moisture index (NDMI).Worldwide, scientists have been intensively focused on how, where, and what magnitude of land use/land cover (LULC) or land indices dynamics have been influencing the climatic conditions as a result of LST intensification and SUHI, which directly or indirectly make the environment uncomfortable and unhealthy for all animals and plants.Some studies rigorously found similar facts in Taipei city of Taiwan [15], Phoenix city of the United States of America (USA) [16], Singapore [10], Dhaka city of Bangladesh [17], Kathmandu valley of Nepal [18], Nanjing city of China [19], Beijing city of China [20], Tokyo city of Japan [21], Tehran city of Iran [13], 70 selected cities of Europe [22], Hong Kong [23], Baltimore-DC metropolitan area of the USA [24], and Cairo city of Egypt [25].At the same time, researchers have discussed how their changing aspects affect and transform the environment of the city landscape in various Indian cities, such as Kolkata [26], Delhi and Mumbai [27], Chandigarh [28], Hyderabad [29], Noida [30], Lucknow [31], and Raipur [32].Mal et al. (2020) conducted a study on the relationship of LST with LULC and elevation in the Ganga River basin, which includes major cities such as Kolkata, Patna, Allahabad (now Prayagraj), Varanasi, Lucknow, Kanpur, New Delhi, and Kathmandu during 2001-2019 using 1 km of MODIS Terra datasets [33].Other studies over the Ganga River basin were also carried out for different cities such as Delhi [34], Kolkata [35], Kanpur and Patna [36], and Lucknow [30] using the MODIS/Landsat database.However, this study lacked a city-level analysis of LST profiling and SUHI information, especially for Prayagraj city.Furthermore, this study lacked effective land indices such as the NDBI, EBBI, NDVI, NDMI, normalized difference water index (NDWI), and soil-adjusted vegetation index (SAVI) to unfold the land dynamics and their role in the LST increase.Furthermore, the analytical results did not show details for thermal state analysis, including land indices dynamics and SUHI information, on a long spatiotemporal scale for Prayagraj city of India.
UN-Habitat (2018) introduced 17 sustainable development goals (SDGs).SDG-11 has stresses the city's resilience and sustainability and focuses on the significance of greenery and open spaces in bringing environmental viability and prosperity by coping with adverse local climate change and intense landscape transformation [37].Various studies have shown that two main approaches help study the LST of urban climate, i.e., ground observations (GOBs) and satellite observations (SOBs).The GOB involves conventional data calculation of air temperature using urban and rural meteorological stations.The SOB has spatiotemporal resolutions with a mathematical background for estimating LST, and they are required to study the spatial variations of SUHI [3].Therefore, in this study, we plan to use spatiotemporal Landsat imageries  to derive land indices, LST, and SUHI information in Prayagraj city using summer season (May-June) and winter season (December-January) datasets.We selected these times for our study on the basis of the finest spectral signature availability and albedo because of haze-free and cloud-free skies [38].
The primary aim of this present work is to examine the interrelations with LST dynamics using directional profiling on the summer/winter seasons during 1987-2018 in Prayagraj city of India by investigating local level climate change through a long spatiotemporal analysis of six land indices (NDBI, EBBI, NDVI, NDMI, NDWI, and SAVI).In this connection, the tasks are (i) to assess the thermal state over the city landscape along with six land indices, (ii) to explore the interrelationship between six land indices and LST, (iii) to delineate the role of the six different land indices in LST increase using urban-rural directional profiling, (iv) to extract the SUHI state both at daytime and night-time, and (v) to extract the hotspots using Getis-Ord G i * statistics.The dynamics of six land indices, LST, and their correlation are investigated to achieve these tasks.Then, the scenario of SUHI formation is discussed using directional profiling of LST by assessing how these six land indices impact the dynamics of LST in eight directions, namely, west, east, north, south, southwest, northeast, northwest, and southeast (center of the city to periphery).Next, the interrelationships between land indices and LST dynamics are validated using Google Earth images.Lastly, we delineate the hotspots to find warming or cooling spaces scattered over the city landscape.

Study Area
We selected Prayagraj city as a study area because this city was selected as a smart city by the MoHUA, i.e., the Ministry of Housing and Urban Affairs, Government of India (GoI), in 2015 [39].This city is one of the biggest in terms of size, as well as historically enriched cities, in Uttar Pradesh state in India.Its location ranges from 25 • 23 7 N to 25 • 32 14 N latitude and 81 • 43 57 E to 81 • 53 59 E longitude, where the mean elevation is 93.77m (Figure 1).The study area covers 72.98 km 2 .This city is located over the holy place called Sangam (confluence of Ganga, Yamuna, and invisible Saraswati rivers) [40,41].The sides of these rivers (Ganga and Yamuna) have been enriched by eroded materials from the Vidhyan uplands and Himalayas mountains [42].The study area has a CWG-type climate, i.e., monsoon type with dry winters based on Koppen's climatic regions (KCR) scheme in India [43], and it has 744.1 mm of mean annual rainfall and 20-32.6 • C of mean annual temperature [44].This city is widely known for religious gatherings at 6 year intervals, Maha Kumbh Mela and Ardha Kumbh Mela, where >100 million pilgrims congregate to make it the largest congregation in Asia [40].

Data Used
This study uses Landsat 5 (TM) and Landsat 8 (OLI/TIRS) satellite imageries with a spatial resolution of 30 m.These are employed for four distinct time points for distinctive seasons in summer and winter.This study uses night-time LST derived from MODIS (Terra) satellite datasets to investigate summer/winter seasonal LST dynamics and SUHI state at night-times from 2007 to 2018 only as this satellite has been providing imagery since 2000.Before 2000, there were no available data on night-time LST.Table 1 shows the satellite datasets used in this study.The software, ERDAS IMAGINE 2014 was used for preprocessing these satellite images.The same dry summer and dry winter seasons were selected for obtaining the cloud-free data with the finest spectral information.

Data Used
This This study uses night-time LST derived from MODIS (Terra) satellite datasets to investigate summer/winter seasonal LST dynamics and SUHI state at night-times from 2007 to 2018 only as this satellite has been providing imagery since 2000.Before 2000, there were no available data on night-time LST.Table 1 shows the satellite datasets used in this study.The software, ERDAS IMAGINE 2014 was used for preprocessing these satellite images.The same dry summer and dry winter seasons were selected for obtaining the cloud-free data with the finest spectral information.

Methods
Figure 2 shows the overall methodological framework for the execution process.All the methods were discussed under subsequent heads.

Methods
Figure 2 shows the overall methodological framework for the execution process.All the methods were discussed under subsequent heads.

Land Indices NDBI
In the evaluation of urban climate, NDBI is an important indicator.It ranges from −1 to +1.A higher positive value specifies bare soils.Lower positive and negative values adjacent to 0 show the dominance of built-up space.Higher negative values specify vegetation and water bodies.NDBI is calculated using Equation (1).
where, in Landsat 5 TM, NIR Band is band 4, and MIR Band is band 5, whereas, in Landsat 8, NIR Band is band 5 and MIR Band is band 6.

EBBI
EBBI is widely used as a significant indicator in assessing urban climate.It ranges from 0 to +1.A higher positive value (perceived threshold > 0.1) specifies bare soils.A lower positive value (perceived threshold between 0.06 and 0.1) specifies built-up space.m positive value (perceived threshold < 0.06) adjacent to 0 specifies water bodies and vegetation.EBBI is calculated using Equation (2) [17].
where, in Landsat 5 TM, NIR Band is band 4, MIR Band is band 5, and TIR Band is band 6, whereas, in Landsat 8, NIR Band is band 5, MIR Band is band 6, and TIR Band is band 10.

NDMI
Another critical indicator in urban climate assessment is NDMI.It ranges from −1 to +1.A positive value specifies vegetation and water bodies.A negative value indicates bare soils and built-up areas.It provides information on the moisture present in the landscape.NDMI is calculated using Equation (3) [47].
where, in Landsat 5 TM, NIR Band is band 4, and MIR Band is band 5, whereas, in Landsat 8, NIR Band is band 5, and MIR Band is band 6.

NDVI
NDVI ranges from −1 to +1.A higher positive value specifies vegetation.A lower positive value indicates bare soils or built-up areas.The negative values specify water bodies.NDVI is calculated using Equation (4) [12].
where Landsat 5 TM defines NIR Band as band 4 and Red Band as band 3, whereas Landsat 8 OLI/TIRS defines NIR Band as band 5 and Red Band as band 4.

NDWI
NDWI ranges from −1 to +1.A positive value specifies water bodies.A lower positive value adjacent to 0 specifies vegetation space.Negative values specify bare soils and built-up areas.NDWI is calculated using Equation ( 5) [48].

NDWI =
Green Band − NIR Band Green Band + NIR Band (5) where, in Landsat 5 TM, Green Band is band 2, and NIR Band is band 4, whereas, in Landsat 8, Green Band is band 3, and NIR Band is band 5.

SAVI
SAVI ranges from −1 to +1.A positive value specifies vegetation.A lower positive value adjacent to 0 specifies water bodies.Negative values specify bare soils and built-up space.It is calculated using Equation (6) [49].
where Landsat 5 TM defines NIR Band as band 4 and Red Band as band 3, whereas Landsat 8 OLI/TIRS defines NIR Band as band 5 and Red Band as band 4. L is the soil brightness correction factor, and it is a constant 0.5.

LST Retrieval Landsat-Based LST Calculation
Several algorithms are available for retrieving LST for distinct satellite sensors.The most popular algorithms are the mono-window algorithm (MWA) [50], radiative transfer equation (RTE) algorithm [51], split-window algorithm (SWA) [52][53][54], and single-channel algorithm (SCA) [55,56] for retrieving LST using thermal bands of Landsat.RTE cannot be considered for use if the information of the atmospheric profile is not available on in situ parameters at the satellite pass [57].The SWA provides accurate results but was not selected for use in the study area as it is specific to band 10 of Landsat 8 (OLI/TIRS) data only for LST computation because of its better calibration.The MWA and SCA also give good results [33,58].In this study, we selected MWA for LST computation from multitemporal Landsat images because MWA shows significant accuracy for computing LST, with three essential parameters being indispensable to LST retrieval, i.e., the ground emissivity, the atmospheric transmittance (AT), and the effective mean temperature of the atmosphere [59], calculated using Equations ( 7)- (9), respectively.
where T s defines the LST (K), T a defines the mean atmospheric temperature (K), T b defines the at-sensor pixel brightness temperature (K), C and D algorithm parameters are estimated through land surface emissivity (LSE) and AT, ε defines LSE, τ defines AT, and a and b are constants of the algorithm (−67.355351 and 0.458606, respectively).The following steps are needed to calculate Equations ( 7)-( 9).
where L λ represents spectral radiance at the top of atmosphere (TOA) (W/(m 2 •sr•µm)), Λ represents the multiplicative rescaling factor for each specific band in the metadata, Γ represents the additive rescaling factor for each specific band in the metadata, and QCAL represents the quantized and calibrated digital number (DN) values of standard product.QCAL max and QCAL min represent the maximum and minimum DN values of the images, respectively.L max and L min represent the TIR band's spectral radiance at QCAL max and QCAL min , respectively.These values of the rescaling factor are available in the metadata of respective Landsat images.
Then, the at-sensor brightness temperature (BT) is computed using Equation (12).
where T b represents the at-sensor BT (K), and K 1 and K 2 (Wm −2 ) are thermal conversion constants (prelaunch calibration) mentioned in the metadata (Table 1) of the respective sensors of Landsat datasets.
Step 2: To compute LST, LSE is one of the indispensable parameters [60].The NDVI threshold (NDVI THR ) method was selected to calculate the LSE because of its significance in segregating pixels of vegetation, water, and soil [45].LSE can be computed using Equations ( 13)- (15).
where ε represents LSE, ε v represents vegetation emissivity, ε s represents soil emissivity, P V represents proportionate of vegetation, C represents constant of surface characteristics, NDV I s represents NDVI of pure soil, NDV I v represents NDVI of pure vegetation, and F represents a geometric factor commonly considered as 0.55 [45,61].
The constant values of LSE were calculated using Equation ( 16) for Landsat 5 (TM) and Equation (17) for Landsat 8 (OLI/TIRS) [61], where ρ Red represents the reflectance value of respective red bands of the imageries.
Step 3: The AT is another indispensable parameter to calculate LST.Before calculating AT, water vapor content should be calculated on the basis of the atmospheric profile (Table 2) using Equation (18) [59,61]  Step 4: The effective mean atmospheric temperature is another indispensable parameter for computing LST (Table 3).It can be computed using Equation (20).

MODIS-Based Night-Time LST Calculation
Night-time LST was retrieved using MODIS night-time datasets.First of all, the MODIS night-time dataset reprojection was changed.Then, cloud-affected areas were eliminated using preprocessed quality control.Each pixel DN value was then converted into LST ( • C) using Equation ( 22) [46].
where x i defines the values of land indices (NDBI, EBBI, NDVI, NDMI, NDWI, and SAVI), and y i defines the LST values.

Intensity of SUHI Calculation
The SUHI is defined by the observed difference of LST between the urban space and the suburban/rural space over a city landscape.It was previously well defined by Oke and East [64] and Oke [65] using Equation (24), which is very popular in the literature [9,13].
where T U−R is the intensity of the SUHI, T U represents the LST of urban space, and T R represents the LST of suburban/rural space.

Hotspot Analysis (Getis-Ord G i *)
The spatial LST distribution over the city landscape was examined using hotspot analysis (Getis-Ord G i *) to characterize both hot and cold spots over the city using each feature (LST value) on the basis of its neighboring features.Hotspots are the clustered areas of high values of the feature, whereas cold spots are the clustered areas of low values of the feature.The Getis-Ord G i * statistic is derived from Equations ( 25)-( 27) [66].
where x j is defined by the feature attribute value of j, w i,j is defined by the spatial weight between features i and j, and n is defined by the total number of features.

Seasonal Spatiotemporal LST Dynamics
The summer spatiotemporal distribution of LST dynamics for Prayagraj city is shown in Figure 3a for summer time points S1, S2, S3, and S4 with boxplots in Figure 3b and their statistics in Table 4.The mean LST witnessed was 38.20 • C in S1, which increased to 40.44 • C in S2, but declined to 37.49 • C in S3, before again inclining to 38.09 • C in S4.In S1, the foremost area of warm temperature was northeast, and that of cool temperature was southeast.In S2, the foremost area of warm temperature was west (except barren land in the northwest at 8-8.5 km), and that of cool temperature was southeast (except for the area of Ganga River flow in the east at 8-9 km).In S3, the foremost area of warm temperature was southwest, and that of cool temperature was northwest (except for the area of the Ganga River flow in the northeast at 8-9 km).However, in S4, the foremost area of warm temperature was the northeast.The foremost area of cool temperature was the northwest (except for vegetation coverage area in the northeast at 1.5-2.5 km).We recommend further research on the surface types or land use/land cover (LULC) classes of the city landscape in Sarif and Gupta (2022) to additionally determine their distribution over the city landscape [67].The winter spatiotemporal LST is mapped in Figure 4a showing the distribution of LST dynamics.Their statistics are presented in Table 4 and Figure 4b.The mean LST witnessed was 19.72 • C in W1, which declined to 19.41 • C in W2.It further declined to 18.06 • C in W3.However, again, it inclined to 19.84 • C in W4.In W1, the foremost area of warm temperature was the northwest, and that of cool temperature was North (except for the area of the Ganga River flow in the northeast at 8-9 km).In W2, the foremost area of warm temperature again was the northwest, and that of cool temperature was the east (except for the area of the Ganga River flow in the northeast at 8-9 km).In W3, the foremost area of warm temperature was the northwest, and that of cool temperature was the northeast.In W4, the foremost area of warm temperature was the southeast, and that of cool temperature was the northwest (except for the area of vegetation coverage in the northeast at 1.5-2.5 km).
temperature was the northeast.The foremost area of cool temperature was the north (except for vegetation coverage area in the northeast at 1.5-2.5 km).We recommend ther research on the surface types or land use/land cover (LULC) classes of the city scape in Sarif and Gupta (2022) to additionally determine their distribution over th landscape [67].

Seasonal Magnitude of LST Based on Multiple Ring Profiling
The magnitude of the summer mean LST difference between the periods based on multiple ring profiling was extracted at 0.5 km intervals from the city's center to the city's periphery (Table 5 and Figure 5).It was detected that each zone at each period witnessed

Seasonal Magnitude of LST Based on Multiple Ring Profiling
The magnitude of the summer mean LST difference between the periods based on multiple ring profiling was extracted at 0.5 km intervals from the city's center to the city's periphery (Table 5 and Figure 5).It was detected that each zone at each period witnessed a higher temperature in comparison to its preceding time points (except for the period between S2 and S3) by a significant amount.Substantial temperature intensification was experienced in each distinct zone from the city center to the periphery.However, a declining trend was observed in the period of S1-S2.The magnitude of the winter mean LST difference distributional pattern between distinctive periods on the basis of multiple ring profiling was also extracted at 0.5 km intervals from the center of the city to the periphery of the city (Table 5 and Figure 6).It was detected that each zone at each period witnessed lower temperatures in comparison to its preceding time points by a significant amount except for the period between W3 and W4.The magnitude of the winter mean LST difference distributional pattern between distinctive periods on the basis of multiple ring profiling was also extracted at 0.5 km intervals from the center of the city to the periphery of the city (Table 5 and Figure 6).It was detected that each zone at each period witnessed lower temperatures in comparison to its preceding time points by a significant amount except for the period between W3 and W4.This means that a substantial temperature reduction was experienced in each zone from the city center to the periphery.However, an amplifying trend was observed during the period of W3-W4.6.Furthermore, Figure A1 shows whisker boxplots of all the land indices for the summer/winter seasons.The analysis found that summer mean NDBI witnessed a decrease of −0.023 in S1 which amplified to −0.015, −0.02, and 0.02 in S2, S3 and S4, respectively.The The spatial NDBI distributional dynamics maps of Prayagraj are shown in Figure 7a,b for the summer season for S1, S2, S3, and S4 time points and the winter season for W1, W2, W3, and W4, respectively.The statistics of all the six land indices and their relationship with LST for all the summer and winter time points are presented in Table 6.Furthermore, Figure A1 shows whisker boxplots of all the land indices for the summer/winter seasons.The analysis found that summer mean NDBI witnessed a decrease of −0.023 in S1 which amplified to −0.015, −0.02, and 0.02 in S2, S3 and S4, respectively.The winter mean NDBI witnessed −0.036 in W1 which amplified to −0.012 in W2 but declined to −0.018 in W3 and further to −0.039 in W4. to −0.018 in W3 and further to −0.039 in W4.Most of the highly dense impervious/built-up area was intensely concentrated in the city center to the south in S1 and W1, which further expanded in the south in S2, W2, S3, and W3.In S4 and W4, the highly dense built-up area spread in the southwest and northeast in S4 and W4, respectively.Therefore, it is apparent that the built-up area is mostly concentrated in the southwest and northeast up to 8 km.Most of the highly dense impervious/built-up area was intensely concentrated in the city center to the south in S1 and W1, which further expanded in the south in S2, W2, S3, and W3.In S4 and W4, the highly dense built-up area spread in the southwest and northeast in S4 and W4, respectively.Therefore, it is apparent that the built-up area is mostly concentrated in the southwest and northeast up to 8 km.
The correlations between LST and all six land indices for the summer and winter seasons are presented in Figures 13 and 14, respectively.In the summer/winter seasons, the correlations between LST vs. NDBI were found to be positive at all distinctive time points.On the basis of the values of correlation coefficients, it can be found that the NDBI's role in LST intensification in the summer was higher than in the winter.Further, the built-up intensity effect was amplified over the city landscape, resulting in increased temperature growth because of the conversion of forests and water bodies into built-up land.

EBBI Dynamics and Its Connection with LST
The spatiotemporal maps of EBBI dynamics are shown in Figure 8a,b for the summer/winter seasons, respectively, for all the distinctive summer/winter time points (however, summer/winter seasonal statistics are shown in Figure A1).The summer mean EBBI witnessed an amplifying trend in the summer/winter seasons at all time points.The highly dense impervious/built-up area was mostly concentrated in the city center to 4 km of its periphery, and the bare land was mostly concentrated to 6-8 km and 9.5-11 km in both S1 and W1.Then, the highly dense impervious/ built-up area expanded in the south and north up to 6 km in S2 and W2, and bare land mostly remained concentrated in the northeast (6-8 km) and southwest (6-8 km).In S3 and W3, the highly-dense impervious/built-up area expanded in the south and north up to 7 km, and bare land mostly remained concentrated in northeast (6-8 km) and southwest (7-9 km).However, the highly dense impervious/built-up area expanded in the southwest, northeast, and northwest up to 8 km in S4 and W4, and bare land mostly remained concentrated on northeast (6-8 km) and southwest (8-9 km).Therefore, it is apparent that the impervious/built-up area mostly existed in the southwest and northeast up to 7 km.
north up to 6 km in S2 and W2, and bare land mostly remained concentrated in the northeast (6-8 km) and southwest (6-8 km).In S3 and W3, the highly-dense impervious/builtup area expanded in the south and north up to 7 km, and bare land mostly remained concentrated in northeast (6-8 km) and southwest (7-9 km).However, the highly dense impervious/built-up area expanded in the southwest, northeast, and northwest up to 8 km in S4 and W4, and bare land mostly remained concentrated on northeast (6-8 km) and southwest (8-9 km).Therefore, it is apparent that the impervious/built-up area mostly existed in the southwest and northeast up to 7 km.A positive correlation was detected between LST vs. EBBI at all distinctive summer/winter time points (Figures 13 and 14).The role of EBBI in LST intensification remained high in both the summer and the winter (Table 6).This indicates that the impervious/built-up and bare land intensity effect was amplified over the city landscape, resulting in increased temperature growth due to the conversion of forests, water bodies, and agricultural land into built-up land and bare land settings.

NDMI Dynamics and Its Relationship with LST
Seasonal NDMI dynamics are shown in Figure 9a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown in Figure A1).The summer mean NDMI witnessed a nonuniform pattern where it first amplified in S1 and S2 but declined in S3, before again amplifying in S4.In contrast, the winter mean NDMI shows an increasing pattern in W1 and W2 but a decreasing one in W3 and W4.
This may be attributed to the fact that the high moisture content area was concentrated in the northwest and northeast in S1 and W1.A decrease was observed in the moisture content area in the northwest in S2 and W2 which further diminished in S3 and W3.However, the high moisture content in the northwest and southeast was again amplified A positive correlation was detected between LST vs. EBBI at all distinctive summer/winter time points (Figures 13 and 14).The role of EBBI in LST intensification remained high in both the summer and the winter (Table 6).This indicates that the impervious/built-up and bare land intensity effect was amplified over the city landscape, resulting in increased temperature growth due to the conversion of forests, water bodies, and agricultural land into built-up land and bare land settings.

NDMI Dynamics and Its Relationship with LST
Seasonal NDMI dynamics are shown in Figure 9a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown in Figure A1).The summer mean NDMI witnessed a nonuniform pattern where it first amplified in S1 and S2 but declined in S3, before again amplifying in S4.In contrast, the winter mean NDMI shows an increasing pattern in W1 and W2 but a decreasing one in W3 and W4.
This may be attributed to the fact that the high moisture content area was concentrated in the northwest and northeast in S1 and W1.A decrease was observed in the moisture content area in the northwest in S2 and W2 which further diminished in S3 and W3.However, the high moisture content in the northwest and southeast was again amplified in S4 and W4.In fact, the forest, grassland, and Ganga River with high moisture content were primarily present in the city landscape.
In both seasons, the correlation between LST vs. NDMI was found to be negative at all time points (Figures 13 and 14), while the correlation coefficient values (Table 6) indicate that the role of NDMI in decreasing LST was higher in the summer season than in the winter season.
in S4 and W4.In fact, the forest, grassland, and Ganga River with high moisture content were primarily present in the city landscape.In both seasons, the correlation between LST vs. NDMI was found to be negative at all time points (Figures 13 and 14), while the correlation coefficient values (Table 6) indicate that the role of NDMI in decreasing LST was higher in the summer season than in the winter season.

NDVI Dynamics and Its Relationship with LST
The NDVI maps of Prayagraj are shown in Figure 10a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown in Figure A1).An increase was observed in the summer mean NDVI at S2, followed by a decrease in S3 but again an increase in S4.However, the winter mean NDVI witnessed a decline in W2 and W3 but an increase in W4.
The analysis shows that high-density vegetation was concentrated in the northwest in S1 and W1, showing a diminishing trend in S2, W2, S3, and W3 but an amplification in S4 and W4.It was found that the forest and grassland were mainly dominant in the northwest direction from the city center.

NDVI Dynamics and Its Relationship with LST
The NDVI maps of Prayagraj are shown in Figure 10a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown in Figure A1).An increase was observed in the summer mean NDVI at S2, followed by a decrease in S3 but again an increase in S4.However, the winter mean NDVI witnessed a decline in W2 and W3 but an increase in W4.In the summer season, the correlations between LST vs. NDVI were found to be negative at all summer time points (Figure 13), while, in the winter season, the correlation between LST vs. NDVI was also positive at all winter time points except for W2 (Figure 14).This reflects that NDVI played a significant role in decreasing LST in the summer season, but its role was very weak in the winter.The analysis shows that high-density vegetation was concentrated in the northwest in S1 and W1, showing a diminishing trend in S2, W2, S3, and W3 but an amplification in S4 and W4.It was found that the forest and grassland were mainly dominant in the northwest direction from the city center.
In the summer season, the correlations between LST vs. NDVI were found to be negative at all summer time points (Figure 13), while, in the winter season, the correlation between LST vs. NDVI was also positive at all winter time points except for W2 (Figure 14).This reflects that NDVI played a significant role in decreasing LST in the summer season, but its role was very weak in the winter.

NDWI Dynamics and Its Relationship with LST
The NDWI maps of Prayagraj are shown in Figure 11a,b for the summer and winter seasons for all time points, respectively (however, summer/winter seasonal statistics are shown in Figure A1).The summer and winter mean NDWI followed a similar pattern, first declining in S2 and W2, then amplifying in S3 and W3 before again declining in S4 and W4.In both summer and winter seasons, the correlations between LST vs. NDWI was positive at all times except in W2 but with relatively low correlation coefficient values (Figures 13 and 14).This means that NDWI played an insignificant role in decreasing LST in both summer and winter because the whole city exhibited a lack of water bodies except for 8-9.5 km in the northeast direction of the Ganga River flow.High-density water bodies were concentrated in the northeast (8-9.5 km) due to the Ganga River flow in S1 and W1.The high-density vegetation was concentrated in the northwest and southeast, which decreased in S2, W2, S3, and W3.Then, high-density vegetation increased in the northwest and southeast in S4 and W4, respectively.Therefore, it is apparent that the forest and grassland were mainly dominant in the northwest direction from the city center.
In both summer and winter seasons, the correlations between LST vs. NDWI was positive at all times except in W2 but with relatively low correlation coefficient values (Figures 13 and 14).This means that NDWI played an insignificant role in decreasing LST in both summer and winter because the whole city exhibited a lack of water bodies except for 8-9.5 km in the northeast direction of the Ganga River flow.

SAVI Dynamics and Its Relationship with LST
The distributional dynamics maps of SAVI are shown in Figure 12a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown in Figure A1).The summer mean SAVI amplified in S2 but declined in S3 before again amplifying in S4.The winter mean SAVI declined in W2 and W3 but amplified in W4.In the summer season, the correlation between LST vs. SAVI was found to be negative at each distinctive summer time point (Figure 13), while, in the winter season, the correlation between LST vs. SAVI was found to be positive at all winter time points except in W2 (Figure 14).This means that SAVI played a noteworthy role in decreasing LST in the summer season, but its role was very weak in the winter.It can be observed that high-density forest and grassland were concentrated in the northwest and southeast in S1 and W1, which was reduced in S2 and W2.The high-density forest and grassland concentrated in the northwest diminished in S3 and W3.However, the high-density forest and grassland concentrated within the northwest and southeast increased in S4 and W4.Therefore, it was observed that the forest and grassland mostly dominated in the northwest direction from the city center.
In the summer season, the correlation between LST vs. SAVI was found to be negative at each distinctive summer time point (Figure 13), while, in the winter season, the correlation between LST vs. SAVI was found to be positive at all winter time points except in W2 (Figure 14).This means that SAVI played a noteworthy role in decreasing LST in the summer season, but its role was very weak in the winter.

Effects of Land Indices on LST Distribution
The effect of spatiotemporal seasonal dynamics of all six land indices (i.e., NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI) on LST profiling for summer/winter seasons from the city center to the periphery in eight different directions, i.e., north to south, northeast to southwest, northwest to southeast, and west to east was extracted on the city landscape.The effect of these six land indices on the LST in the summer and winter seasons is presented in Figures 15 and 16

Effects of Land Indices on LST Distribution
The effect of spatiotemporal seasonal dynamics of all six land indices (i.e., NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI) on LST profiling for summer/winter seasons from the city center to the periphery in eight different directions, i.e., north to south, northeast to southwest, northwest to southeast, and west to east was extracted on the city landscape.The effect of these six land indices on the LST in the summer and winter seasons is pre-

North to South
In summer, from the city center in both the north and the south directions, the mean LST followed a declining trend (up to 11 km) of 0.2 • C-2.75 • C except for 6.5-11 km in the north direction and 9-11 km in the south direction (Figure 15).The higher mean LST at the city center was because of the impervious/built-up land (high mean NDBI by −0.04 to 0.04).The decline in mean LST witnessed at 1-2 km in the north and 4.5-6.5 km in the south was due to the existence of forest land/grassland (higher mean NDVI by 0.02-0.10),as well as the availability of higher moisture content (higher mean NDMI and mean SAVI).However, the increase in mean LST at 6.5-7.5 km and 10-11 km in the north and 6.5-11 km in the south was because of the existence of barren land/bare soils (higher mean EBBI by 0.01-0.08).The mean NDWI was found to be negative (<−0.1) at all consecutive summer time points, which means that its role was insignificant due to the unavailability of water bodies of large size (except rivers) to impact the reduction in mean LST.
In winter, from the city center in both the north and the south directions, the mean LST witnessed declining inclination up to 11 km by 0.2 • C-3.6 • C except for 6.5-9.5 km in the north, and 5-9 km and 10-11 km in the south (Figure 16).The higher mean LST in the city center was due to the built-up land presence (high mean NDBI by −0.04 to 0.04).In contrast, the decline in mean LST witnessed at 1-2 km and 9.5-11 km in the north and 9-10 km in the south was because of the presence existence of grassland/forest land (higher mean NDVI by 0.02-0.25)coupled with high moisture content availability, as reflected by high mean NDMI and mean SAVI values.The existence of barren land/bare soils (higher mean EBBI by 0.01-0.08)at 6.5-9.5 km in the north and at 5-9 km and 10-11 km in the south resulted in a spike in mean LST.However, as in the summer season, mean NDWI had negative values in the winter season also [<0 in the north and south (except at city center and 1.5-2 km in the south)] at all consecutive winter time points, signifying its insignificant role in the reduction in mean LST.

Northeast (NE) to Southwest (SW)
In summer, in both the NE and the SW directions from the city center, mean LST detected a declining pattern up to 11 km by 0.2 • C-5.8 • C except for 6-7.5 km in the NE and 6.5-11 km in the SW, where a reverse trend was witnessed.The decline in mean LST witnessed at 1-3 km in the NE and 6.5-7.5 km and 9.5-10.5 km in the SW was because of the availability of grassland/forest land (higher mean NDVI by 0.02-0.13)and high moisture, as the mean NDMI and mean SAVI also witnessed a similar pattern to mean NDVI.The availability of barren land/bare soils (higher mean EBBI by 0.01-0.1)resulted in the increase in mean LST at 6.5-7.5 km in the NE and 7.5-9.5 km in the SW.The mean NDWI showed negative values in the summer, indicating that its role was insignificant due to the unavailability of water bodies.However, 7.5-8 km in the NE showed the reverse trend due to the Ganga River flow.
In winter, in both the NE and the SW directions from the city center, the mean LST witnessed a declining pattern up to 11 km by 0.2 • C-4.5 • C except for 6-7.5 km in the NE direction and 10-11 km in the SW direction where the mean LST witnessed an increasing trend (Figure 16).The city center had high-density built-up land (higher mean NDBI by −0.04 to 0.04), leading to higher mean LST which witnessed a decline at 1-3 km in the NE and 4.5-6 km and 10-11 km in the SW because of the existence of grassland/forest land (higher mean NDVI by 0.02-0.2).The mean NDMI and mean SAVI also witnessed a similar pattern to that of mean NDVI, indicating the availability of higher moisture due to the presence of forest/grassland, which helped in the decline in mean LST.However, the mean LST high spikes at 6-7.5 km in the NE direction and at 10-11 km in the SW direction were due to barren land/bare soils (higher mean EBBI by 0.01-0.08)at these locations.However, in summer, the role of NDWI was insignificant except for 7.5-8.5 km in the NE due to Ganga River.

Northwest (NW) to Southeast (SE)
In summer, in both the NW and the SE directions from the city center, it was detected that the mean LST witnessed a declining pattern up to 9 km in the NW and 8 km in the SE by 0.1 • C-2.6 • C because of the existence of grassland/forest land (higher mean NDVI by 0.01-0.15)coupled with high moisture availability due to high mean NDMI and mean SAVI (Figure 15).However, an increasing pattern was witnessed in the zones of 3.5-4.5 km and 6-9 km in the NW and 6.5-8 km in the SE because of the presence of barren land/bare soils (higher mean EBBI by 0.01-0.2).The role of the mean NDWI was insignificant due to negative values (<0.1) at all consecutive summer time points.
In winter, in both the NW and the SE directions from the city center, the mean LST witnessed a declining pattern up to 11 km by 0.1 • C-3.5 • C except for 3-5 km and 6.5-9 in the NW direction and 5-8 km in the SE direction, where the mean LST witnessed an increasing pattern.The decline in mean LST at 1-3 km in the NW and 4.5-6.5 km and 3-5 km in the SE was because of the presence of grassland/forest land (higher mean NDVI by 0.02-0.23)and the high peaks of mean NDMI and mean SAVI (Figure 16).The barren land/bare soils resulted in an increase in mean LST inclination at 3-5 km and 6.5-9 km in the NE and at 5-8 km in the SE (higher mean EBBI by 0.01-0.06).Again, the role of NDWI was insignificant due to negative values of the mean NDWI [except the city center to 1.5 km in the SE at W1, where the mean LST was >0] at all consecutive time points.

West to East
In summer, the mean LST witnessed a higher peak at the city center in comparison to nearby surroundings because of the presence incidence of high-density impervious/builtup land (higher mean NDBI by −0.04 to 0.04).Then, the mean LST followed a declining pattern up to 11 km by 0.1 • C-1.75 • C in both the west (in particular, at 1-3 km and 5.5-6.5 km) and the east (in particular, 3-5 km and 8-9.5 km) due to forest land/grassland (higher mean NDVI by 0.02-0.10)coupled with high moisture availability (Figure 15).The mean LST witnessed an increasing pattern in 4-5.5 km and 6.5-11 km in the west and 6-7.5 and 9.5-11 km in the east because of barren land/bare soils (higher mean EBBI by 0.01-0.2).However, mean NDWI was found to have negative values (<-0.1) at all consecutive summer time points, which means that its role was insignificant due to the unavailability of many water bodies.
In winter, the city center again witnessed higher LST compared to its surroundings.However, the mean LST witnessed declining inclination up to 11 km by 0.1 • C-1.6 • C [except 4-5.5 km and 6.5-11 km in the west and 5.5-7.5 km and 9.5-11 km in the east] in both the west and the east directions.The decline in mean LST witnessed at 1.5-3 km and 5.5-6.5 km in the west and 4-6 km and 8.5-9.5 km in the east was because of the existence of grassland/forest land (higher mean NDVI by 0.02-0.22).The mean NDMI and mean SAVI also witnessed a similar pattern to the mean NDVI, which means that higher moisture was available where forest/grassland was present, consequently helping the decline in the mean LST (Figure 16).However, the mean LST inclination at 4.5-5 km and 6.5-11 km in the west direction and at 5.5-7.5 km and 9.5-11 km in the east direction was because of barren land/bare soils presence (higher mean EBBI by 0.01-0.05).However, the mean NDWI was found to have negative values [<0 in the north and south (except at city center to 0.5 km in the west)] at all consecutive winter time points, which means that its role was insignificant due to the unavailability of water bodies large in size, which may have impacted the reduction in the mean LST.The distribution of SUHI dynamics was evaluated in the summer season at all four time points (S1, S2, S3, and S4), as well as in the winter season at all four time points (W1, W2, W3, and W4) over Prayagraj city.For this purpose, five rural/suburban point locations, well distributed over the study area (shown in Figures 3 and 4) for sum-mer/winter seasons, were selected to study the actual variation in LST between urban and rural/suburban areas.The seasonal SUHI dynamics for all time points of the summer and winter seasons at five rural/suburban locations in Prayagraj city are shown in Figure 17.Furthermore, the statistics of point location-based SUHI are presented in Table 7.It can be observed that these five case points showed strong SUHI evidence for the summer and the winter.
3, 15, 179 30 of 42 locations, well distributed over the study area (shown in Figures 3 and 4) for summer/winter seasons, were selected to study the actual variation in LST between urban and rural/suburban areas.The seasonal SUHI dynamics for all time points of the summer and winter seasons at five rural/suburban locations in Prayagraj city are shown in Figure 17.Furthermore, the statistics of point location-based SUHI are presented in Table 7.It can be observed that these five case points showed strong SUHI evidence for the summer and the winter.In the present work, directional ring profiling was used to delineate the actual difference of LST from the city center to its periphery by taking eight-directional ring profiling, i.e., north to south, NE to SE, NW to SE, and west to east.The distribution of LST dynamics was evaluated to investigate the SUHI state in the summer and winter at all four time points over Prayagraj city.Figure A2 shows the seasonal LST profiling for SUHI formation in Prayagraj city  for both the summer and the winter seasons.

North to South
In summer, at S1, the mean LST was higher in the city center than in both the north and the south directions by ~0.75-1.5 °C [except for 5-11 km in the north and 6.5-11 km in the south because these areas were constituted mostly of sands and bare soil].At S2, the mean LST was higher in the city center than in both the north and the south directions by ~0.2-2.75 °C [except for 9.5-10.5 km in the north and 7-11 km in the south because  In the present work, directional ring profiling was used to delineate the actual difference of LST from the city center to its periphery by taking eight-directional ring profiling, i.e., north to south, NE to SE, NW to SE, and west to east.The distribution of LST dynamics was evaluated to investigate the SUHI state in the summer and winter at all four time points over Prayagraj city.Figure A2 shows the seasonal LST profiling for SUHI formation in Prayagraj city  for both the summer and the winter seasons.

North to South
In summer, at S1, the mean LST was higher in the city center than in both the north and the south directions by ~0.75-1.5 • C [except for 5-11 km in the north and 6.5-11 km in the south because these areas were constituted mostly of sands and bare soil].At S2, the mean LST was higher in the city center than in both the north and the south directions by ~0.2-2.75 • C [except for 9.5-10.5 km in the north and 7-11 km in the south because these areas were mostly sands and bare soils].At S3, the mean LST was higher at the city center than in both the north and the south directions by ~0.2-2 • C [except for 6.5-7 km, 8-8.5 km, and 9.5-10.5 km in the north, and 10-11 km in the south because these areas were mostly sands and bare soils].At S4, the mean LST was higher at the city center than in both the north and the south directions by ~0.2-2.2 • C except for 6.5-7.5 km and 9.5-10.5 km in the north, and 6.5 km and 10-11 km in the south because these areas were mostly sands and bare soils.
In winter, at W1, the mean LST was higher in the city center than in both the north and the south directions by ~0.2-1.7 • C [except for 7.5-11 km in the north, and 6.5-11 km in the south because these areas were mostly sands and bare soil].At W2, the mean LST was higher at the city center than in both the north and the south directions by ~0.2-1.5 • C [except at 6 km in the north and 10-11 km in the south because these areas were mostly sands and bare soils].At W3, the mean LST was higher in the city center than in both the north and the south directions by ~0.2-3.6 • C [except for 8-8.5 km and 9.5-10.5 km in the north, and 10-11 km in the south because these areas were covered mostly with sands and bare soils].At W4, the mean LST was higher in the city center than in both the north and the south directions by ~0.3-2.1 • C except for 10.5-11 km in the south because of the existence of barren land/bare soils.

NE to SE
In summer, the mean LST showed a higher peak at the city center than in both the NE and the SW directions at all four time points, i.e., S1, S2, S3, and S4, by ~0.1-1.2 • C, ~0.1-4.1 • C, ~0.12-5.8• C, and ~0.2-1.8,respectively.However, because of the existence of sand, barren land, and bare soil, a reverse pattern was observed in zones 5-7.5 km and 9-11 km in the NE, and 5-11 km in the SW at the S1 timepoint, in zones 6-6.5 km and 10-11 km in the NE, and 5.5 km and 6.5-11 km in the SW at the S2 timepoint, in zones 7 km and 9.5-10.5 km in the NE and 10-11 km in the SW at the S3 timepoint, and in zones 6.5-8 km in the NE and 6.5 km and 10-11 km in the SW at the S4 timepoint.
In winter, the mean LST in the city center also followed the same pattern as that of the summer, with higher temperatures observed as compared to suburban areas in both the NE and the SW directions.The mean LST in the city center was higher by ~0.1-3.25 • C, ~0.1-2.4 • C, ~0.2-4.5 • C, and ~0.3-2.1 • C at the W1, W2, W3, and W4 time points.However, in the NE direction, the zones 7-7.5 km at W1, 5.0 and 7.5 km at W2, 6.5-7 km, 8.5 km, and 10 km at W3, and 6-7 km at W4 exhibited a higher LST peak as compared to the city center due to the presence of sand and bare soil.Along the same line, in the SW direction, zones 4.5-11 at W1, 7-8.5 km and 9.5-11 km at W2, 4.5-5.5 km, 8 km, and 9.5-11 km at W3, and 10-11 km at W4 exhibited the same trend.

NW to SE
In summer, at S1, the mean LST was higher at the city center than in both the NW and the SE directions by ~0.15-1.8• C except for zones having predominantly bare soil at 3.5-5 km and 6.5-9 km in the NW and 4.5-8 in the SE.The same pattern was also observed at the S2, S3, and S4 time points.The mean LST at the city center was higher by ~0.1-2.6 • C, ~0.2-1.3 • C, and ~0.3-2 • C in both directions at S2, S3, and S4, respectively.The zones with a reverse trend, due to the presence of predominantly bare soil, were 3.5-4.5 km and 7.5-9 km in the NW and 6-8 km in the SE at S2, 3.5-4.5 km and 7.5-9 km in the NW and 5 km and 6-8 km in the SE at S3, and 3.5-4.5 km and 7.5-9 km in the NW and 5.5-8 km in the SE at S4.
In winter, the city center observed a higher LST peak than the suburban areas in both the NW and the SE directions, as observed in the summer season.The mean LST in the city center was higher by ~0.1-1 • C, ~0.1-1.9 • C, ~0.1-3.5 • C, and ~0.3-1.7 • C at the W1, W2, W3, and W4 time points, respectively.However, exceptions to this trend were also observed in certain zones in both the NW and the SE directions at each distinct timepoint, as can be clearly witnessed in Figure A2, because these areas were predominantly exposed to bare soil.

West to East
In summer, the city center observed higher LST than suburban areas by ~0.1-1.2 • C at S1, ~0.1-1.75 • C at S2, ~0.1-1.7 • C at S3 and ~0.1-1.5 • C at S4 in both the west and the east directions.However, because of the existence of mostly sand, barren land, and bare soil, a few zones observed a reverse trend, namely, 3.5-11 km in the west and 5-11 km in the east at S1, 4.5-5 km and 7-11 km in the west, and 7.5-8 km and 9.5-11 km in the east at S2, 4.5-5.5 km and 7-11 km in the west, and 6.5-7 km and 9.5-11 km in the east at S3, and 10-11 km in the west, and 7-8 km and 10.5-11 km in the east at S4.
In winter, the mean LST in the city center also followed the same pattern as that of the summer, with higher temperatures observed as compared to suburban areas in both the west and the east directions.The mean LST in the city center was higher by ~0.1-0.85 • C, ~0.1-1.6 • C, ~0.1-1.5 • C, and ~0.1-1.3 • C at the W1, W2, W3, and W4 time points.However, exceptions to this trend were also observed in certain zones in both the west and the east directions at each distinctive timepoint, as can be clearly seen in Figure A2, as these areas mostly had sand and bare soil.

Hotspot Identification
The hotspot analysis was performed using the Getis-Ord G i * approach to analyze the spatial distribution of LST over Prayagraj city.This approach uses LST values of neighboring features and delineates both hot and cold spots over the city landscape.Hotspots are the clusters of features of high values of LST, while cold spots aggregate the features of low LST values.On the basis of this analysis, the city landscape was categorized into seven classes: very cold, cold, cool, not significant, warm, hot, and very hot.
The summer spatiotemporal hotspot maps of Prayagraj city are shown in Figure 18a for summer time points S1, S2, S3, and S4 to present the clustering distribution of hotspot dynamics, and their statistics are compiled in Table 8.The very cold spot class experienced a high loss of 1.29 km 2 of areal coverage during S1-S4.The cold spot class experienced a loss of 0.07 km 2 of areal coverage during S1-S4.The cool spot class also experienced a loss of 0.44 km 2 of areal coverage during S1-S4.However, the not significant class experienced an enormous gain of 7.56 km 2 of areal coverage during S1-S4.The warm spot class experienced a loss of 1.81 km 2 of areal coverage during S1-S4.The hot spot class experienced a loss of 2.06 km 2 of areal coverage during S1-S4.The very hot spot class experienced a loss of 1.88 km 2 of areal coverage during S1-S4.This summer hotspot pattern indicates that the comfort level of living space intensively decreased in the city landscape as the areas of very cold, cold, and cool spots severely declined.
The winter spatiotemporal hotspot maps of Prayagraj city are shown in Figure 18b for the W1, W2, W3, and W4 winter time points to present the clustering distribution of hotspot dynamics.Their statistics are shown in Table 8.The very cold spot class experienced a severe loss of 6.49 km 2 of areal coverage during S1-S4.The cold spot class experienced a loss of 0.80 km 2 of areal coverage during S1-S4.The cool spot class also experienced a loss of 2.08 km 2 of areal coverage during S1-S4.However, the not significant class has experienced a considerable gain of 15 km 2 of areal coverage during S1-S4.The warm spot class experienced a loss of 1.46 km 2 of areal coverage during S1-S4.The hot spot class experienced a loss of 1.30 km 2 of areal coverage during S1-S4.The very hot spot class experienced a loss of 2.87 km 2 of areal coverage during S1-S4.This winter hotspot pattern also indicates that the comfort level of living space was on the decline in the city landscape as the areas of very cold, cold, and cool spots severely declined.This work assessed the seasonal (summer and winter) thermal state over the city landscape of Prayagraj city (India).The effects of land indices, namely, NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI, on the thermal state were extensively examined to investigate how water bodies, forest land, wetland, and barren soils control intensification and/or cooling of LST over the landscape of the study area.The study area was delineated using eight-directional ring profiling of land indices, including LST to explore how, where, and what magnitude the LST changed either in increasing or in decreasing patterns due to different land indices.This can help policymakers and planners to conduct sustainable planning and enrich the carrying capacities of the landscape.As per the IPCC AR6 report of 2021, at the local to the global level, an extreme transformation of landscape has occurred, especially in the postindustrial era, which needs to be monitored and mitigated to control the rise in global mean temperature, leading to a decrease in the adverse consequences of climate change [5].In connection to this issue, multiple cities (such as Taipei city of Taiwan [15], Phoenix city of the United States of America (USA) [16], Singapore [10], Dhaka city of Bangladesh [17], Kathmandu valley of Nepal [18], Nanjing city of China [19], Beijing city of China [20], Tokyo city of Japan [21], Tehran city of Iran [13], 70 selected cities of Europe [22], Hong Kong [23], Baltimore-DC metropolitan area of the USA [24], and Cairo city of Egypt [25]) have witnessed a similar amplifying pattern, which is very concerning and critical for our blue planet.UN-Habitat (2018), through SDG-11, categorically acknowledged the significance of spatial identification of land coverage and their possible effects on the safety, resilience, and sustainability of the city landscape, especially using greenery and open spaces [37].
Worldwide, to minimize the adverse effects of SUHI, including the rise in energy consumption, water scarcity, air pollution, and health problems, e.g., sunstroke, cardiac, and respiratory issues [68,69], populaces are converging on different mitigation strategies, such as the use of light paints and materials, planting of trees on streets, and cool and green roof creations [70,71].Other strategies, such as designing the size, orientation, and shape of buildings, could improve the wind flow in the city landscape [72].Moreover, distinctive mitigation approaches, such as preserving wetlands, small to large water bodies, and greenery plantations on barren spaces, can significantly improve the local climate and ecosystem of the city landscape [38,73,74].For this strategy and subsequent planning and creation, remote sensing and GIS-based information using a multitemporal spatial database can be applied to recognize the ground reality.Along with this, public awareness regarding the adoption and implementation of the above strategies is a prerequisite for their success in sustainable development and ecosystem restoration of the Prayagraj city landscape.

Conclusions
This study explored the dynamics of seasonal (summer and winter) land indice [namely, NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI] and LST dynamics using Landsa 5 (TM) and Landsat 8 (OLI/TIRS) imagery for Prayagraj city of India.The multitempora The spatiotemporal winter night-time LST maps are shown in Figure 19b.Winter time points, such as January 2007 and December 2018, were incorporated to present the night-time LST dynamics and their whisker boxplot statistics.The mean night-time LST was also intensified by 1.69 • C during the winter periods from 2007 to 2018.At night-time, the central urban core area intensively experienced higher LST than its periphery.It was further detected that the SUHI phenomenon at night-time was severely intensified by 0-2.98 • C and 0-4.56 • C in summer and winter, respectively, according to multiple ring profiling.Other major cities over the Indo-Gangetic plain, such as Delhi [34], Lucknow [30], Patna [36], and Kolkata [35] have also been experiencing a strong SUHI phenomenon, whereby the city center has exhibited a higher LST than the periphery in the last few decades; a similar pattern was observed in our study area.
The above findings of night-time LST dynamics indicate a severe SUHI occurrence in both the summer and the winter seasons at night-time.These findings further strengthen the results of the daytime LST of Landsat imagery from 2008 to 2018.In the daytime, impervious surfaces absorb albedo extensively, and a long time is taken to radiate it back to the atmosphere due to its physical properties.As a result of increasing urbanization, this phenomenon has been further intensified.Both daytime and night-time LST experience a strong SUHI phenomenon; however, at night-time, SUHI becomes more vulnerable.As per this study, it is suggested to take immediate attention to reduce the SUHI severity using an effective mitigation strategy after urban planners and policymakers consider the aforementioned spatial thermal anomalies in the city.

Conclusions
This study explored the dynamics of seasonal (summer and winter) land indices [namely, NDBI, EBBI, NDMI, NDVI, NDWI, and SAVI] and LST dynamics using Landsat 5 (TM) and Landsat 8 (OLI/TIRS) imagery for Prayagraj city of India.The multitemporal spatial pattern of land indices and LST and their correlation dynamics with directional ring profiling over the city landscape were investigated, including the formation of SUHI and its dynamics over the city landscape.It was found that summer periodical LST magnitudes were highly intensified during S1-S4 by 0.32-2.45• C (except for zone 11 km) at 0.5 km intervals from the city center to the periphery.Winter periodical LST magnitudes were also intensified by 0.02-1.06• C (except for zones 9, 11, and 15-19 km).It was witnessed that the northeast and southwest directions had a high growth of LST distribution, while the northwest direction had a low growth of LST distribution in both seasons.
The results based on directional ring profiling of the effect of land indices on LST found that most of the vegetation/forest land available, at 1-3 km in the northwest and 5-6 km in the southeast direction, were depleted during the selected period.The impervious/built-up land expanded from the city center to 8 km in all directions during the study time, whereas bare soils and sand were primarily present in the northeast and the northwest (6-11 km).The presence of different land covers significantly controlled the LST distribution as forested area decreased the LST distribution whereas built-up area, bare soils, and sands increased the LST distribution.
Forest cover played a crucial role in declining the LST by 2.25-4.8• C (except for water bodies), whereas bare soils and sand played a critical role in amplifying the LST by 1.9-5.6 • C. At the same time, built-up land amplified the LST by 1.8-3.9• C (except for sand and bare soils).These results were further strengthened by the findings that LST vs. NDVI, LST vs. SAVI, LST vs. NDMI, and LST vs. NDWI had a positive correlation in the summer season.However, the LST vs. NDWI relationship had a very weak positive correlation due to the unavailability of water bodies of a significant size, which may have reduced the LST.In contrast, LST vs. NDBI and LST vs. EBBI had a strong positive correlation.

Figure 3 .
Figure 3. Spatiotemporal dynamics of LST in the summer season over Prayagraj city: (a) LST during 1988-2018 (five rural/suburban areas shown here used for computing SUHI) and (b plots of LST dynamics during 1988-2018.

Figure 3 .
Figure 3. Spatiotemporal dynamics of LST in the summer season over Prayagraj city: (a) LST maps during 1988-2018 (five rural/suburban areas shown here used for computing SUHI) and (b) boxplots of LST dynamics during 1988-2018.

Figure 4 .
Figure 4. Spatiotemporal dynamics of LST (°C) in winter season over Prayagraj city: (a) LST maps during 1987-2018 (five rural/suburban areas shown here used for computing SUHI) and (b) boxplots of LST dynamics during 1987-2018.

Figure 4 .
Figure 4. Spatiotemporal dynamics of LST ( • C) in winter season over Prayagraj city: (a) LST maps during 1987-2018 (five rural/suburban areas shown here used for computing SUHI) and (b) boxplots of LST dynamics during 1987-2018.

Figure 5 .
Figure 5. Magnitude of the summer mean LST in different periods (1988-2018).Between S1 and S2, the mean LST observed amplifying pattern ranged from 0.94 • C to 3.62 • C, wherein the lowest intensified zone was 12 (5.5-6km) by 0.94 • C because there was existence of grass and forest land.The highest intensified zone was 0 (city center) by 3.62 • C because of the presence of highly dense, impervious land.Between S2 and S3, the difference in mean LST observed a declining pattern which ranged from −1.80 • C to −3.88 • C, wherein the lowest decline zone was 12 (5.5-6km) by −1.80 • C because of the dominance of grass and forest land, but the highest declined zone was 15 (7-7.5 km) by −3.88 • C because of the dominance of sand and barren land.Between S3 and S4, the difference in mean LST was found in an amplifying pattern that ranged from 0.33 • C to 2.65 • C, wherein the lowest intensified zone was 22 (10.5-11km) by 0.33 • C because of the presence of grassland.The highest intensified zone was again 0 (city center) by 2.65 • C because of the presence of highly dense, impervious land.Between S1 and S4, the difference in mean LST showed an increasing pattern which ranged from 0.12 • C to 2.45 • C except for zone 11 (5-5.5 km), in which the magnitude of mean LST declined by −0.10 • C because of forest and grassland, wherein the lowest intensified zone was 17 (8-8.5km) by 0.12 • C because of the presence of the Ganga River flow.The highest intensified zone again was 0 (city center) by 2.45 • C because of the presence of highly dense, impervious land.The magnitude of the winter mean LST difference distributional pattern between distinctive periods on the basis of multiple ring profiling was also extracted at 0.5 km intervals from the center of the city to the periphery of the city (Table5and Figure6).It was detected that each zone at each period witnessed lower temperatures in comparison to its preceding time points by a significant amount except for the period between W3 and W4.This means that a substantial temperature reduction was experienced in each zone from the city center to the periphery.However, an amplifying trend was observed during the period of W3-W4.
km) by 2.48 °C because of the dominant barren land.Between W1 and W4, the difference in mean LST exhibited a very interesting amplifying pattern [except for zone 9 (4-4.5 km), zone 11 (5-5.5 km), and zone 15 to zone 19 (7-9.5 km), which ranged from −0.89 °C to −0.21 °C as these zones had forest and grassland, as well as Ganga River flow] which ranged from 0.02 °C to 1.06 °C, wherein the lowest intensified zone was 8 (3.5-4 km) by 0.02 °C because there was the existence of impervious land.The highest intensified zone again was 22 (10.5-11km) by 1.06 °C due to the dominant barren land.

3. 3 .
Spatiotemporal Dynamics of Land Indices and LST, and Their Relationships 3.3.1.NDBI Dynamics and Its Connection with LST The spatial NDBI distributional dynamics maps of Prayagraj are shown in Figure 7a,b for the summer season for S1, S2, S3, and S4 time points and the winter season for W1, W2, W3, and W4, respectively.The statistics of all the six land indices and their relationship with LST for all the summer and winter time points are presented in Table

Figure 6 .
Figure 6.Magnitude of the winter mean LST in different periods (1987-2018).Between W1 and W2, the difference in mean LST showed a declining pattern (except zone 1 by 0.45 • C and zone 2 by 0.22 • C) which ranged from −0.76 • C to −0.03 • C, wherein the lowest decline zone was 2 (0.5-1 km) by −0.03 • C because there was the existence of high-density impervious land.The highest decline zone was 11 (5-5.5 km) by −0.76 • C because of the presence of forest and grassland.Between W2 and W3, the difference in mean LST was found to decline, which ranged from −2.24 • C to −0.48 • C, wherein the lowest decline zone was 20 (city center) by −0.48 • C because of the presence of bare soil with scattered grassland.The highest decline zone was 0 (city center) by −2.24 • C because of the high moisture content over scattered vegetation in the high-density impervious land.Between W3 and W4, the difference in mean LST was found amplified in the range 0.64 • C to 2.48 • C, wherein the lowest intensified zone was 17 (8-8.5km) by 0.64 • C because of the presence of Ganga River flow.The highest intensified zone was 22 (10.5-11km) by 2.48 • C because of the dominant barren land.Between W1 and W4, the difference in mean LST exhibited a very interesting amplifying pattern [except for zone 9 (4-4.5 km), zone 11 (5-5.5 km), and zone 15 to zone 19 (7-9.5 km), which ranged from −0.89 • C to −0.21 • C as these zones had forest and grassland, as well as Ganga River flow] which ranged from 0.02 • C to 1.06 • C, wherein the lowest intensified zone was 8 (3.5-4 km) by 0.02 • C because there was the existence of impervious land.The highest intensified zone again was 22 (10.5-11km) by 1.06 • C due to the dominant barren land.3.3.Spatiotemporal Dynamics of Land Indices and LST, and Their Relationships 3.3.1.NDBI Dynamics and Its Connection with LST

3. 3 . 6 .
SAVI Dynamics and Its Relationship with LST The distributional dynamics maps of SAVI are shown in Figure 12a,b for the summer and winter seasons, respectively (however, summer/winter seasonal statistics are shown

42 Figure 13 .
Figure 13.Correlation between the LST and the six land indices of Prayagraj city in the summer season (1988-2018).

Figure 13 . 42 Figure 14 .
Figure 13.Correlation between the LST and the six land indices of Prayagraj city in the summer season (1988-2018).

Figure 14 .
Figure 14.Correlation between the LST and the six land indices of Prayagraj city in the winter season (1987-2018).

Figure 15 .
Figure 15.Effects of the land indices on the distribution of LST dynamics in the summer season in Prayagraj city (1988-2018).

Figure 15 .
Figure 15.Effects of the land indices on the distribution of LST dynamics in the summer season in Prayagraj city (1988-2018).

Figure 16 .
Figure 16.Effects of the land indices on the distribution of LST dynamics in the winter season in Prayagraj city (1987-2018).

Figure 16 .
Figure 16.Effects of the land indices on the distribution of LST dynamics in the winter season in Prayagraj city (1987-2018).

Figure 17 .
Figure 17.Seasonal SUHI dynamics at five case point locations in Prayagraj city (1987-2018) for all time points: (a) summer and (b) winter.

Figure 17 .
Figure 17.Seasonal SUHI dynamics at five case point locations in Prayagraj city (1987-2018) for all time points: (a) summer and (b) winter.

Figure 19 .
Figure 19.Maps of night-time LST dynamics with whisker boxplots for exploring SUHI at night-time: (a) summer and (b) winter.

Table 1 .
Summary of the data used.

(Sensor)/Ancillary Data Path/Row Resolution/Scale Season Acquisition Date Time (GMT) Constants of Thermal Conversion Source K 1 K 2 Landsat-5 (TM)
. Then, AT is computed using Equation (19).

Table 2 .
Estimation equations of atmospheric transmittance.

Table 3 .
The estimation equation for effective mean atmospheric temperature in four standard atmospheres.

Table 4 .
Summer/winter seasonal LST dynamics of Prayagraj city.

Table 5 .
Periodical LST magnitudes based on multiple ring buffers from the city center to the periphery at 0.5 km of intervals in Prayagraj city.

Table 6 .
Statistics of the six land indices and their relationship with LST for all the summer and winter time points.

Table 6 .
Statistics of the six land indices and their relationship with LST for all the summer and winter time points.

Table 7 .
Seasonal SUHI magnitude based on rural/suburban point locations in Prayagraj city.

Table 7 .
Seasonal SUHI magnitude based on rural/suburban point locations in Prayagraj city.

Table 8 .
Seasonal hotspot classes based on Getis-Ord G i * statistics. Hot-