The Spatiotemporal Implications of Urbanization for Urban Heat Islands in Beijing: A Predictive Approach Based on CA–Markov Modeling (2004–2050)

The rapid increase in infrastructural development in populated areas has had numerous adverse impacts. The rise in land surface temperature (LST) and its associated damage to urban ecological systems result from urban development. Understanding the current and future LST phenomenon and its relationship to landscape composition and land use/cover (LUC) changes is critical to developing policies to mitigate the disastrous impacts of urban heat islands (UHIs) on urban ecosystems. Using remote sensing and GIS data, this study assessed the multi-scale relationship of LUCC and LST of the cosmopolitan exponentially growing area of Beijing, China. We investigated the impacts of LUC on LST in urban agglomeration for a time series (2004–2019) of Landsat data using Classification and Regression Trees (CART) and a single channel algorithm (SCA), respectively. We built a CA–Markov model to forecast future (2025 and 2050) LUCC and LST spatial patterns. Our results indicate that the cumulative changes in an urban area (UA) increased by about 908.15 km2 (5%), and 11% of vegetation area (VA) decreased from 2004 to 2019. The correlation coefficient of LUCC including vegetation, water bodies, and built-up areas with LST had values of r = −0.155 (p > 0.419), −0.809 (p = 0.000), and 0.526 (p = 0.003), respectively. The results surrounding future forecasts revealed an estimated 2309.55 km2 (14%) decrease in vegetation (urban and forest), while an expansion of 1194.78 km2 (8%) was predicted for a built-up area from 2019 to 2050. This decrease in vegetation cover and expansion of settlements would likely cause a rise of about ~5.74 ◦C to ~9.66 ◦C in temperature. These findings strongly support the hypothesis that LST is directly related to the vegetation index. In conclusion, the estimated overall increase of 7.5 ◦C in LST was predicted from 2019–2050, which is alarming for the urban community’s environmental health. The present results provide insight into sustainable environmental development through effective urban planning of Beijing and other urban hotspots.


Introduction
Land use/cover change (LUCC) most often surrounds natural vegetation alteration for logging, urbanization, and agriculture expansion [1,2]. LUCC generates many environmental problems at both local and global scales, including biodiversity loss due to greenhouse gas release [2][3][4], changes in land surface temperature (LST), and precipitation shifts. The negative environmental consequences of urbanization, which include population growth, large-scale industrial and infrastructure development, and rapidly changing landscapes, are a global concern [5,6]. LUCC associated with urbanization is one of the major causes of shifting LST. LUCC and LST hold great importance due to the consequential impacts of LST on the urban environment. China, the world's second-largest economy, is undergoing rapid urbanization with dramatic infrastructure and urban area growth [5,7]. The efficient development of urbanization policies aimed at the adjustment of land use structure and  (Table 1). Before LUCC classification, the Landsat images were pretreated to remove atmospheric effects [24,26,34,35]. For LUCC classification accuracy, a ground survey was conducted to collect 250 ground control points for each land cover class. The digital elevation model (DEM) was collected from SRTM by using GEE. Each spatial scene was enhanced using the histogram equalization approach to attain a higher image contrast [36,37]. A series of steps were completed in a step-wise workflow pattern ( Figure 2). First, information sets were processed in GEE to construct a False Color Composite (FCC). The study area was extracted from all spatial imaginaries via masking of Beijing's geo-referenced outline boundary map. The Supervised Classification method was used to improve classification results from Landsat images. Statistical inferences from the mean LST and the percentage proportion of different land cover types, vegetated and non-vegetated areas, were drawn through correlation analysis from 2004 to 2019. The CA-Markov model was applied to predict future (2025 and 2050) trends of LUCC and LST.  (Table 1). Before LUCC classification, the Landsat images were pretreated to remove atmospheric effects [24,26,34,35]. For LUCC classification accuracy, a ground survey was conducted to collect 250 ground control points for each land cover class. The digital elevation model (DEM) was collected from SRTM by using GEE. Each spatial scene was enhanced using the histogram equalization approach to attain a higher image contrast [36,37]. A series of steps were completed in a step-wise workflow pattern ( Figure 2). First, information sets were processed in GEE to construct a False Color Composite (FCC). The study area was extracted from all spatial imaginaries via masking of Beijing's georeferenced outline boundary map. The Supervised Classification method was used to improve classification results from Landsat images. Statistical inferences from the mean LST and the percentage proportion of different land cover types, vegetated and non-vegetated areas, were drawn through correlation analysis from 2004 to 2019. The CA-Markov model was applied to predict future (2025 and 2050) trends of LUCC and LST.

Land Use/Cover Change
The Supervised Classification method is an algorithmic probability program applied to land cover classification. Classification and Regression Trees (CART) [21,38,39] is a primer supervised classification algorithm used in remote sensing techniques of spectral images. Ground verification in uncertain areas was completed through Google Earth Pro and GEE, where misclassified areas were corrected by positioning and rearranging the GEE script samples. The ground truth point was used to estimate mapping accuracy. Finally, Landsat-5 (TM) and Landsat-8 OLI (ETM) were classified into six land cover types using the CART classification algorithm (Table 2). Producers and users assessed the remote sensing image; overall, a confusion matrix was calculated for the classification accuracy and kappa coefficient. Producer precision is the number of pixels that have been correctly identified as a percentage of the classifier's pixels as a training sample of each batch. User precision is the pixels that have been accurately categorized as percentages of the total number of pixels identified as the class. Overall classification accuracy is the ratio of the total number of pixels correctly assigned to the total number. Equation (1) defines the kappa coefficient. The number of random points used to determine accuracy for each image was 250.

Land Use/Cover Change
The Supervised Classification method is an algorithmic probability program applied to land cover classification. Classification and Regression Trees (CART) [21,38,39] is a primer supervised classification algorithm used in remote sensing techniques of spectral images. Ground verification in uncertain areas was completed through Google Earth Pro and GEE, where misclassified areas were corrected by positioning and rearranging the GEE script samples. The ground truth point was used to estimate mapping accuracy. Finally, Landsat-5 (TM) and Landsat-8 OLI (ETM) were classified into six land cover types using the CART classification algorithm (Table 2). Producers and users assessed the remote sensing image; overall, a confusion matrix was calculated for the classification accuracy and kappa coefficient. Producer precision is the number of pixels that have been correctly identified as a percentage of the classifier's pixels as a training sample of each batch. User precision is the pixels that have been accurately categorized as percentages of the total number of pixels identified as the class. Overall classification accuracy is the ratio of the total number of pixels correctly assigned to the total number. Equation (1) defines the kappa coefficient. The number of random points used to determine accuracy for each image was 250.

Calculation of Normalized Difference Vegetation Index (NDVI)
The Normalized Difference Vegetation Index (NDVI) is the most commonly used vegetation index to observe greenery globally. The Near-Infrared (NIR) and red band ratio are also used to map vegetation to test its condition using Equation (2) [40]. The magnitude of that index is between −1 and 1. For green plants, the standard range is between 0.2 and 0.8 [41].
where L λ is the spectral radiance in Wm −2 sr −1 mm −1 .
We converted the spectral radiance to at-satellite brightness temperature (TB) under the assumption of uniform emissivity. The conversion formula is given in Equation (4).
T B is the brightness temperature in Kelvin (K), L λ is the spectral radiance in Wm −2 sr −1 mm −1 , and K 2 and K 1 are calibration constants. For Landsat-8 OLI, K 1 is 774.89 and K 2 is 1321.08.
The fractional vegetation, F v , of each pixel was determined from the NDVI using Equation (5) [46,47].
NDVI min 's value (0.2) and pixels considered bare soil, NDVI max is the maximum NDVI value (0.5), and pixels are regarded as healthy vegetation.
dε is the effect of the geometrical distribution of natural surfaces and internal reflections calculated by Equation (6).
where ε v is vegetation emissivity, ε s is soil emissivity, F v is fractional vegetation, and F is a shape factor with a mean of 0.55 [16,38,46,48].

Relative LST Change Detection
For the years 2004 and 2019, the relative LST was calculated to compare the effects of LUCC on the thermal urban climate. RLST changes from LUCC (increase/decrease) are derived from the study region's mean LST and using Equation (10) via increasing pixel value [35,44,49].
RLST jk = LST jk − LST j mean (10) where RLST jk represents the relative temperature of pixel j of class k, LST jk is the temperature of cell j of class k, and LST j indicates the mean value of LST for urban landscape j. If RLST k j > 0, the pixel shows a positive contribution of LUCC conversion, and if RLST k j < 0, then it is a negative contribution to the thermal environment.

Cellular Automata-Markov Chain (CA-Markov) Model Analysis
This model is developed based on a Markov stochastic probability matrix for predicting the transition of one status to another [26,50,51]. Generally, the Markov chain model is used to simulate transitions, parameters, and trends. It generated probability transition matrices to predict and classify potential land use/cover change (LUCC) and urban development scenarios and examined land surface temperature (LST) simulation patterns [20,52,53]. Trends were estimated using Equations (11)-(13) based on the conditional probability formula. Moreover, where S(t) is the state of the system at time t, S(t + 1) is the state of the system at the time (t + 1), and P ij is the matrix of the transition probability in a state. The cellular automata (CA) and Markov chain model is used to calculate LUCC and LST's future scenario by projecting 2025 and 2050. Forecasting LUCC and LST for the projected period was made through the CA-Markov model using the land use change modeler (LCM) in Terrset (Clark Labs TerrSet 18.31). According to the statistics (Table 3), about 44% of the urban area remained unchanged from 2004 to 2019, while 23% loss and 33% gain were observed in LUCC types in the urban landscape. For vegetation overall, 755.10 km 2 (69%) of the area was lost, 288.93 km 2 (26%) was added, and 55.22 km 2 (5%) of sites remained unchanged, which clearly shows the massive changes in urban vegetation. Based on the change detection map ( Figure 4A), about 61.42% of the land use area remained unchanged from 2004 to 2019. Overall, 38.58% of changes were observed in LUCC types in the city landscape, in which 3.7% were observed in CL converted into an urban area, followed by BL (1.23%) into cropland, FA (7.09%) swapped to impervious surfaces, and CL into FA (6.24%) ( Figure 5). A good change was observed in the conversion of BL (1.22%) and WB (0.78%) into an urban area, cropland (0.92%) converted to urban vegetation, and 0.54% of agricultural land replaced by impervious land ( Figure 4A). Additionally, a slight change was found in the conversion of WB (0.4%) into the forest area, VA (0.39%) into impervious surfaces, and BL (0.15%) into FA. Rapid anthropogenic activities have significantly replaced the natural surface area with semi-natural or impervious surfaces and other land uses between 2004 and 2019. According to the statistics (Table 3), about 44% of the urban area remained unchanged from 2004 to 2019, while 23% loss and 33% gain were observed in LUCC types in the urban landscape. For vegetation overall, 755.10 km 2 (69%) of the area was lost, 288.93 km 2 (26%) was added, and 55.22 km 2 (5%) of sites remained unchanged, which clearly shows the massive changes in urban vegetation. Based on the change detection map ( Figure 4A), about 61.42% of the land use area remained unchanged from 2004 to 2019. Overall, 38.58% of changes were observed in LUCC types in the city landscape, in which 3.7% were observed in CL converted into an urban area, followed by BL (1.23%) into cropland, FA (7.09%) swapped to impervious surfaces, and CL into FA (6.24%) ( Figure 5). A good change was observed in the conversion of BL (1.22%) and WB (0.78%) into an urban area, cropland (0.92%) converted to urban vegetation, and 0.54% of agricultural land replaced by impervious land ( Figure 4A). Additionally, a slight change was found in the conversion of WB (0.4%) into the forest area, VA (0.39%) into impervious surfaces, and BL (0.15%) into FA. Rapid anthropogenic activities have significantly replaced the natural surface area with semi-natural or impervious surfaces and other land uses between 2004 and 2019.

Estimation of Land Surface Temperature (LST)
Land surface temperature (LST) for Beijing city (  (Table 4).      The results in Figures 5 and 6  Apart from global warming and climate change, the city's surface temperature was minimal in 2004 because of the proliferation of green space and a low percentage of impervious surfaces spread around the city center. In 2019, the minimum and maximum surface temperatures increased due to urbanization.

Relationship between LUCC and LST
A non-significant negative linear relationship between vegetation and land surface temperature (LST) was observed with an R-value of −0.155(p > 0.419) (Figure 7). Simultaneously, the water index showed a strong negative correlation with land surface temperature, resulting in an R-value of −0.809 (p = 0.000). A significant positive correlation between urban areas and LST was found with an R-value of 0.526 (p = 0.003). Apart from global warming and climate change, the city's surface temperature was minimal in 2004 because of the proliferation of green space and a low percentage of impervious surfaces spread around the city center. In 2019, the minimum and maximum surface temperatures increased due to urbanization.

Relationship between LUCC and LST
A non-significant negative linear relationship between vegetation and land surface temperature (LST) was observed with an R-value of −0.155(p > 0.419) (Figure 7). Simultaneously, the water index showed a strong negative correlation with land surface temperature, resulting in an R-value of −0.809 (p = 0.000). A significant positive correlation between urban areas and LST was found with an R-value of 0.526 (p = 0.003).

Warming and Cooling Impacts of LUCC from 2004 to 2019
LUCC had warming and cooling impacts on urban climate (Figures 8 and 9). The overall maximum warming effects of 5.10 • C, 4.92 • C, and 4.58 • C were observed in the urban area (UA) transformed from barren land (BL), the vegetation area (VA), and agriculture/cropland (CL). Continued urbanization replaced BL with UA and CL, respectively, which raised the RLST by 2.90 • C and 2.84 • C. In the forest area (FA) transition by vegetation, urban, and barren areas, the mild warming effects of 2.82 • C, 2.36 • C, and 1.86 • C were quantified. The minimum positive effect on LST was noticed in the conversion of UA and CL to water bodies.
On the other side, the optimum cooling impact of about −3.20 • C, −2.28 • C, and −1.05 • C was observed in converting the urban area, cropland, and vegetation area into water bodies. A restrained decline in LST of approximately −2.41 • C, −1.05 • C, and −0.66 • C was quantified in the forest area transformed from urban, cropland, and vegetation areas.
The significant positive influence on the RLST is correlated with the transformation of forest areas into impervious and urban regions, while the cooling impact is due to BL and UA being converted into FA. We also investigated moderate changes in surface temperature in the water bodies and grass/agriculture land converted into barren land and urban counties.
The minimum cooling impact of about −0.58 • C and −0.50 • C was observed in converting FA and VA to agriculture/cropland due to the minimum difference in deciduous plantation and urban vegetation with cultivated land.

Warming and Cooling Impacts of LUCC from 2004 to 2019
LUCC had warming and cooling impacts on urban climate (Figures 8 and 9). The overall maximum warming effects of 5.10 °C, 4.92 °C, and 4.58 °C were observed in the urban area (UA) transformed from barren land (BL), the vegetation area (VA), and agriculture/cropland (CL). Continued urbanization replaced BL with UA and CL, respectively, which raised the RLST by 2.90 °C and 2.84 °C. In the forest area (FA) transition by vegetation, urban, and barren areas, the mild warming effects of 2.82 °C, 2.36 °C, and 1.86 °C were quantified. The minimum positive effect on LST was noticed in the conversion of UA and CL to water bodies.

Cellular Automata-Markov Chain (CA-Markov) Model Analysis
The combination of cellular automata (CA) and the stochastic transition matrix of the Markov chain model resulted in LUCC and LST for the projected period of 2025 and 2050 ( Figure 10). Map accuracy for the projected land use/cover change for predictive years was classified by the sufficient kappa coefficient value of 0.97. A decrease of 6% and 11% of forest cover areas was estimated during 2019-2025 and 2019-2050, respectively ( Table 5). The urban area will expand by 5% (2019-2025) and 8% (2019-2050), adding a specific rise in relative temperature for 2025 and 2050. Urban vegetation will decrease to 482.27 km 2 (3%) and 436.57 km 2 (3%), and cropland will increase by 4% and 5%. Between 2019-2025 and 2019-2050, there would be little change in the increase in barren land and water bodies.
During 2019-2025, the area under the slightly warm (20-25 • C) and warm (25-30 • C) categories decreased by about 3.51% and 1.60% (Table 6), which shifted to a positive increase in the hot (30-35 • C) and very hot (>35 • C) categories up to 5.51% and 0.01%, respectively. However, this alarming shift continues under warm and slightly warm areas with a decrease of about 21.60% and 12.15%. Areas under the hot (30-35 • C) and very hot (>35 • C) categories increase by about 26.95% and 1.68%, respectively, during 2025-2050 (Table 7). The significant positive influence on the RLST is correlated with the transformation of forest areas into impervious and urban regions, while the cooling impact is due to BL and UA being converted into FA. We also investigated moderate changes in surface temperature in the water bodies and grass/agriculture land converted into barren land and urban counties.
The minimum cooling impact of about −0.58 °C and −0.50 °C was observed in converting FA and VA to agriculture/cropland due to the minimum difference in deciduous plantation and urban vegetation with cultivated land.

Cellular Automata-Markov Chain (CA-Markov) Model Analysis
The combination of cellular automata (CA) and the stochastic transition matrix of the Markov chain model resulted in LUCC and LST for the projected period of 2025 and 2050 ( Figure 10). Map accuracy for the projected land use/cover change for predictive years was classified by the sufficient kappa coefficient value of 0.97. A decrease of 6% and 11% of forest cover areas was estimated during 2019-2025 and 2019-2050, respectively ( Table  5). The urban area will expand by 5% (2019-2025) and 8% (2019-2050), adding a specific rise in relative temperature for 2025 and 2050. Urban vegetation will decrease to 482.27 km 2 (3%) and 436.57 km 2 (3%), and cropland will increase by 4% and 5%. Between 2019-2025 and 2019-2050, there would be little change in the increase in barren land and water bodies.     Figure 10. The predicted LULC and LST maps for the years 2025 and 2050. Figure 10. The predicted LULC and LST maps for the years 2025 and 2050.

Implication of Land Use/Land Cover Change for LST
Although relationships between land surface temperature (LST) and land use and cover change (LUCC) have been studied previously [26,28,54], they lack the application of remote sensing technology to infer patterns and relationship dynamics. Further, 2004 to 2019 provides a unique, nearly two-decade span of extreme growth and economic change in Beijing's large metropolitan area. Moreover, we also simulated these parameters for 2025 and 2050 to provide future urban development plans [5,26,55]. The results provide insight into testable hypotheses, the quantities of individual contributions of LUCC to LST in hotspot areas, and potential mitigation measures to combat consequential adverse effects. The findings have shown that urban sprawl is the primary driving force in land surface temperatures (LSTs). This radiation and heat exchange is essential in rising urban heat islands (UHIs) [5,56,57]. Land use and land cover change (LUCC) impact RLST, especially in urban areas [58]. LUCC is essential to active management for awareness of anthropogenic climate change and global warming [24,26,59].
Currently, the available literature has shown that land cover classification could estimate the existing relationship between LST and LUCC [18,22]. UHIs may have been the product of some essential factors such as macro/meso-climates, urban morphology, population growth, geographical and biophysical shifts in surface area, anthropogenic changes, wind corridors, population challenges, and human lifestyles. We noticed that there was a significant trend in LUCC between 2004 and 2050.
The massive changes seen in LUCC could be attributed to Beijing's rapid development in the 1990s and 2000s. This period was characterized by the intense deforestation and demotion of cropland for various developmental projects. Agricultural land shifted to use as impervious surfaces and built-up areas for housing and industry, characterized by a simultaneous decline in total vegetation cover (VC) [60,61]. This rapid depletion of VC has a wide range of impacts on natural cooling due to the shading and evapotranspiration provided by plants and shrubs [62,63]. To bolster this [24,64], the negative linear relationship between the Normalized Difference Vegetation Index (NDVI) and LST demonstrated that VC acts as a sink within a UHI because of its cooling effects. This change could eventually obliterate the processes of surface evaporation and transpiration in plants [65]. Previous studies have seconded this phenomenon, while its impact on LST reduction in grasslands and ornamental plants is less than in vegetation covered by forests and urban treebanks and gardens [36,66,67].
Our results elucidate that land cover has a dominant impact on LST in urban environments. The impact values of LUCC such as Urban/Built-up areas, urban vegetation, water bodies, and forests on their LST vary according to their proportional area [19,68]. Vegetation is shown to play an essential role in mitigating or controlling temperatures in urban areas [64]. Moisture is added to the surrounding air by evaporation from the top of water bodies. Previous studies have shown that water bodies play a significant role in regulating LST in residential areas [69,70]. Urban areas play an essential role in the development of urban intricate heat flows.
Our study revealed that the vegetation area was not significant since LST has been significantly influenced by growth in developed regions. This may be attributed to the upward extension of external expansion, thereby mitigating vegetation's impact on LST. Temperature values measured in densely vegetated areas were low. In contrast, the highest LST values were observed in barren/impervious land compared to other land cover areas in various urban districts of Beijing [71]. According to previous findings, our results are primarily that Beijing's built-up areas have a strong positive linear relationship with LST [19,61,72]. A higher temperature rise results from the increase in impervious rough, dark surfaces, including dirt, metal, and asphalt. Due to low reflection and intense solar radiation uptake, the building material increased LST, which often emits heat day and night [9,62]. This analysis shows that UHIs have a negative NDVI trajectory [73]. In addition to urban expansion, the downtown area of Beijing has undergone significant changes.
Government and private parties have revitalized a substantial amount of land for new residential, commercial, and industrial plants. This has destroyed traditional wooden houses and the architecture of rice straw roofs and tile roofs and replaced them with skyscrapers and tall buildings made with impermeable, anti-transpiring non-evaporated materials such as concrete, glass, and solid aluminum. These materials directly affect heat fluxes in urban air spaces. Previous studies have shown that LST in China is severely affected by urban rather than rural systems [73][74][75][76], similar to other countries [77,78]. The conversion of forest and agricultural land into urban territories has also contributed to the LST increase. The government has moved several factories and companies to the outskirts of cities to make them more successful in environmental health. The new factories and subsidiary infrastructures are often placed in well-conditioned agricultural or forest land, augmenting the LST of that particular zone. Historically, vegetation or forested areas were considered a fringe line between urban and rural areas to absorb excess heat generated by automobiles and factories harbored within city limits [6,78].
These results work in concert to establish urban sprawl as the primary factor producing an abnormal heat-flux and impacting LST. This radiation and heat exchange is considered a significant effect in UHIs, which causes a substantial contribution to climate change in the city canyon. While LUCC broadly has an intense impact on the relative LST, anthropogenically induced LUCC increases RLST within the urban micro-atmosphere [38].

Land Use Conversion and Its Contribution to UHIs
The land use/land cover changes and associated RLST demonstrate the urban climate's positive and negative contributions at the standardized scale (10%). For example, this study reveals that barren land's transition into metropolitan areas contributes 1.01% to UHIs. In comparison, the warming effect is 1.23% in converting the 0.59% area, which is higher than the warming effect. The reasons for such surprising findings may be the locations of the LUCC close to water bodies and forest land. On the other hand, the conversion of the same area of 0.15% barren land into forest areas contributes a warming effect of 0.41% to UHIs. Similarly, the transition of barren land into vegetation areas pays about 0.61% to UHIs, while the negative input is −0.07 in reverse change (Table 8). The overall contribution is higher than the contribution for ventilation, owing to low to high surface energy. The results show that the warming contribution of the natural surface to the impervious surface is much more significant than the cooling contribution. Overall, this study showed that the maximum contribution to warming is caused by vegetation change, particularly urban forests, to high LST impervious areas. Simultaneously, the minimum result is converting from an impervious surface to grass/agricultural land with minimal surface temperature variation. This reveals that increasing green spaces in urban regions could benefit the urban climate. Simultaneously, the developed sites and impervious parts have low heat transmission capacity and trapped solar energy that can boast the UHI/thermal climate.
Recent scientific literature [47] concluded that the overall contribution of up to~65% warming impact associated with LUCC is well-matched with the expected warming effect due to physical interaction across the Indian region, which supports our findings.
Indeed, the maximum cooling contribution is 0.23% in converting forest areas into water bodies in the land use dynamic process. In comparison, the minimum cooling effect of 0.04% was observed from vegetation change to forests. The purpose behind the least contribution of vegetation area is the minimal difference in surface temperature and proximity to the impervious surface. Meanwhile, the highest cooling is investigated in the conversion of forest areas to water bodies. At the same time, the most negligible impact is obtained from the transformation of vegetation areas to forest land. Some land cover's cooling contribution is changed into an urban area, forest area, water bodies, vegetation, and crop/agriculture land.
Similarly, the inadequate literature on LUCC, associated relative land surface temperature, and their percent contribution makes it hard to make the correct analogy. Recently, a study also seconded our findings by concluding the impact of LUCC on surface temperature and the percent contribution in LST related to LUCC in Donating Lake Area, China [48].

Conclusions
The present study evaluated LUCC's influence on the LST in a large urban area of Beijing, observing the developmental resolution and different socio-economic parameters using RS data. This research established a link between LST and the urban environment. Using coefficient and projections analysis, landscape dynamics' contribution to temperature was evaluated using the Pearson correlation and CA-Markov model. It was discovered that increasing the built-up area by 5% increased the temperature by 1%. Increases in vegetation cover of 10% had a negative correlation as well. For the study period of 2019 to 2050, the rise in land surface temperature (LST) was 7.5 • C (10.35%). Overall, the findings reveal the average warming effect of 2.80 • C and the average cooling effect of −1.41 • C of LUCC in Beijing from 2004 to 2019. Due to forest land transformation to impervious surfaces, the average warming contribution of LUCC to UHI is 0.5%. In comparison, 0.11% cooling was found in reverse direction shifts. The positive contribution in UHIs due to the conversion of forest land into impervious surfaces is higher than the negative contribution.
More targeted green efforts should be placed on urban design, infrastructure planning, and development to negate temperature increases. Enhancing water bodies such as lakes, canals, waterfalls, and fountains and a significant increase in green spaces such as artificial parks, gardens, the linear plantation of woody plants, and the promotion of open-concept areas are all part of this plan. This study indicates that more extensive research is urgently required to assess changing land use/cover on local and regional climate. Many areas are rapidly changing due to the expanding effects of modern climate change and progressive activity. The geographic environment and known weather patterns of Beijing should be harnessed to promote natural cooling processes. In amalgamation, environmental education should be accessible to encourage ecological development concerning resource planning and management. Hence, there is a greater need for green policy intervention and effective urban planning to control the soaring thermal environment. A quantitative examination of these parameters should be included in a prospective study.
Additionally, while we discovered that urbanization has a direct effect on LST, the inverse relationship will reveal the extent to which urbanization affects LST. Future research should focus on quantifying the impacts of urbanization on LST. This study suggests intensive research in the future, specifically investigating the proportion of green space and impervious surfaces and its effects on RLST for sustainable future cities. The present study provides practical implications for urban landscape planning that requires the rational use of landscape connectivity between green and impervious surfaces and their impact on LST. Future urban research could focus on the issue of public health and infrastructure burden associated with rapid urbanization. Funding: This work is funded by the National Natural Science Foundation of China's (Grant Numbers 51978447) study on the response mechanism between the evolution of urbanization policy and the development of rural space networks in the Beijing-Tianjin-Hebei region.
Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: The authors will provide the data of this research upon request.