Multi-Temporal Analysis of the Impact of Summer Forest Dynamics on Urban Heat Island Effect in Yan’an City

: In this study, MODIS land products and China land cover datasets were used to extract normalized difference vegetation index, land surface temperature, and vegetation cover type in Yan’an City during the summers of 2017–2022. On this basis, analysis of spatial change and correlation were carried out as a way to study the mitigation effect on urban heat islands in Yan’an City with forest. The study showed that: (1) The coverage of normalized difference vegetation index over 0.4 in summer in Yan’an City increased from 59.38% to 69.12%, and the vegetation showed good growth conditions. It has a spatial distribution pattern of more in the south and less in the north. (2) The proportion of the urban heat island in Yan’an City increased from 15.51% to 16.86%. Urban heat island intensity fluctuated year by year, with the maximum urban heat island intensity of 6.26 ◦ C appearing in 2019. It has a spatial distribution pattern of less in the south and less in the north. The transition rate of temperature field grade from low to high is 73.32%, and the transition rate to low is only 0.31%. (3) There is a negative correlation between land surface temperature and normalized difference vegetation index in Yan’an City. Vegetation has a mitigating effect on the UHI and the best cooling effect among the vegetation is shown by forest. The cooling effect of forest in Yan’an City is attenuated by an increase in distance, and the effective range is greater than 1000 m. In this study, the regulation effect of forest on the urban heat island was obtained by digging deeper into the intrinsic connection between spatial change in vegetation cover and land surface temperature change in Yan’an City. It provides an important reference for the formulation of meteorological protection policy as well as the promotion of sustainable development of the urban ecological environment and is of guiding significance for future urban planning and ecological construction.


Introduction
The urban heat island (UHI) effect refers to the phenomenon where urban temperatures are notably higher compared to surrounding rural areas [1].Historical UHI data indicate London has experienced rising temperatures since the late 19th century, with hotter summers and a widening temperature gap with its surrounding areas.[2].After the Industrial Revolution, rapid urban population growth, extensive construction and road paving, industrial activities, and anthropogenic heat discharge have escalated.The accelerating urbanization and rapid industrialization have made UHIs a global environmental concern.In the context of global warming today [3], the impact of UHIs on the thermal environment in residential areas can sometimes be more direct than the effects of climate change [4], and there is a mutual influence between the two [5].Excessive temperatures in cities can imbalance the body's regulatory functions (especially in the elderly) and increase the incidence of diseases such as heart attacks and pyrexia.For every 1 K rise in temperature, the number of heat-related deaths in people over 65 years of age is expected to double [6].At the same time, the heat island effect reduces human living comfort, and residents will consequently consume 2-4% more energy for cooling for every 1 • C rise, creating a vicious cycle [7,8].
Since the 20th century, the link strategy of the heat island effect has received more and more attention.The mitigation methods of the heat island effect can reduce the temperature up to 1.5-2.0K on average, which is important for urban policy making and sustainable development [9,10].The denser the building complex in a city, the greater the warming effect, and the more open the urban design, the less likely it is to gather heat [11].In response to the challenge of the UHI effect, greening mitigation initiatives are also gaining attention.In recent years, London has introduced the All London Green Grid, which continuously integrates green space into the urban fabric to reduce urban temperature by linking green spaces [12].In Singapore, green spaces are integrated into buildings and cityscapes to create a "garden city" and a cooling urban microclimate [13].Recent studies have shown that the cooling effect of forests mainly works in two ways.The first is through utilizing the characteristics of the trees themselves to achieve cooling.For example, by using the shading effect of leaves to block shortwave radiation from the sun, reducing ground surface temperatures; by absorbing carbon dioxide through photosynthesis, reducing greenhouse gas emissions, and improving the microclimate in cities; and the transpiration of vegetation can also achieve cooling by increasing air humidity, reducing heat flow in the air, and reducing longwave radiation, among other mechanisms [14][15][16].Second is by slowing down the heat island effect through the reasonable allocation of green space across the landscape spatial layout.Existing studies have confirmed that the proportion and distribution characteristics of forests and shrubs in urban areas can significantly reduce the intensity of urban heat islands.The green space aggregation index and the largest patch index are strongly negatively correlated with ground surface temperature, and reasonable forest arrangements can have a certain cooling effect on urban heat islands.In contrast, the reduction of forest areas leads to a harder urban surface, which makes the urban ground more prone to absorb and store heat, resulting in an exacerbation of the urban heat island effect [17][18][19][20].
Around the world, the majority of research aimed at mitigating the UHI effect is concentrated in densely populated and economically advanced regions [21][22][23][24].Recent studies have infrequently addressed how aspects such as urban layout and green space planning impact UHIs, in which environmental conditions and socioeconomic development levels differ significantly.This paper focuses on Yan'an, a tier-3 city in China known for its significant ecological improvements over the past two decades.This city represents the kind of place that has not seen much research on how to relieve the heat island effect.There is a lack of research concerning the mitigation of UHI effects in similar cities.This research is significant for the development of small and mid-sized cities as it fills the gap in sustainable studies related to strategic land use planning and urban climate management for these regions.Through GIS technology, the study deeply analyzes the state of urban vegetation coverage and the spatial distribution changes of the UHI.It reveals how changes in vegetation affect the temporal and spatial characteristics of the UHI.Qualitative analysis of how vegetation eases the UHI effect in Yan'an fills the gap in research on forests' ability to counter UHIs in less developed areas.Additionally, it compares different types of vegetation to identify the most effective ones for cooling.This research highlights the potential value of various plantings in alleviating the UHI effect, providing theoretical support for greening strategies-like the pattern, area, and distribution of plants in underdeveloped regions.The forests in Yan'an are primarily artificial, dense, and lie at some distance from the city center, raising questions about their temperature regulatory capabilities.This study uses the NDVI and LST data for analysis and modeling.It shows how the forest's cooling effect relates to its distance from the city and identifies an effective range of this cooling.The paper not only explores the value and effectiveness of artificial forests in mitigating the UHI effect but also informs city planning and ecological strategies for less developed regions.

Study Area
The study area is Yan'an City (Figure 1), selected for its distinct climatic conditions and ecological significance.Yan'an City, a renowned historical city with a resident population of about 2,282,581 and an urbanization rate of 61.97%, represents an ideal urban case for studying heat island effects due to its typical temperate monsoon climate and rapid urbanization spells, positioned as a third-tier city in China.The geographical coordinates of Yan'an City range between 109 • 13 ′ ~110 • 49 ′ E and 36 • 25 ′ ~37 • 35 ′ N, with a total area of 37,031.3square kilometers.Because of its diverse vegetation types and distribution, it is known as the "one of the lungs of Shaanxi".According to the survey data released by the Forestry Bureau, the total area covered by forests in Yan'an City is about 19,652.5 square kilometers as of 2022, with a coverage rate of 53.07%, and the total volume of standing trees has reached 86,966,000 cubic meters, most of which are the forests on the mountain ranges, and there are also small urban green islands such as Ganquan Laoshan Forest Park, etc.During the period of 2017-2022, a total of 38,555.10square kilometers of afforestation area was completed in Yan'an City.Yan'an City is within the scope of the Three-North Project, which is known as "the world's largest ecological project".Since the implementation of various ecological restoration projects, the city has continuously realized the "double increase" in forest accumulation and resource area [25,26] and has been one of the cities with the fastest growth of vegetation cover.

Methodological Framework
This study aims to investigate the mitigating effects of forest on the UHI effect in Yan'an City, delineating the methodological approach into three distinct stages.First, LST and NDVI data from MODIS products were downloaded, as well as county and municipal administrative boundaries and annual land cover datasets with a spatial resolution of 30 m.The analysis phase consisted of examining vegetation changes, changes in heat island areas and their intensity, and correlations between NDVI and LST.This phase also included the creation of buffer zones around the affected areas to study the effect of the distance of forests in mitigating UHIs.Ultimately, the results of the study demonstrated the spatial and temporal effects of vegetation on UHIs, identified the most effective types of vegetation mitigation, and determined the effective range of forests in mitigating UHIs.The results of the study provide a basis for decision making on forest layout and ecological conservation strategies in Yan'an City. Figure 2 illustrates the flowchat of this paper.Yan'an City, located in the northern part of Shaanxi Province and in the middle of the Loess Plateau, represents the typical transitional zone between arid and semi-arid regions, making it particularly sensitive to UHIs.It is connected to Yulin in the north, Xianyang, Tongchuan, and Weinan cities in the south, and bordered by the Yellow River and Shanxi in the east and Gansu in the west.Yan'an has a temperate monsoon climate and belongs to the arid and semi-arid transition zone of China, which makes the city conducive to studying the interaction between vegetation dynamics and UHI mitigation.The main climatic characteristics include abundant light and heat resources, dryness, little rain, strong evaporation, and a large temperature difference between day and night, with average annual temperature of 9.9 • C and total annual precipitation of 507.7 mm.The study area is characterized by a wide range of round hills and valleys, a high terrain in the north and a low terrain in the south, and major mountain ranges such as the Baiyu Mountains, the Ziwuling Mountains, and the Huanglong Mountains, with a fragmented surface, large undulations and steep slopes, and significant spatial variations in temperature and precipitation, which have an important influence on the distribution of vegetation cover.

Methodological Framework
This study aims to investigate the mitigating effects of forest on the UHI effect in Yan'an City, delineating the methodological approach into three distinct stages.First, LST and NDVI data from MODIS products were downloaded, as well as county and municipal administrative boundaries and annual land cover datasets with a spatial resolution of 30 m.The analysis phase consisted of examining vegetation changes, changes in heat island areas and their intensity, and correlations between NDVI and LST.This phase also included the creation of buffer zones around the affected areas to study the effect of the distance of forests in mitigating UHIs.Ultimately, the results of the study demonstrated the spatial and temporal effects of vegetation on UHIs, identified the most effective types of vegetation mitigation, and determined the effective range of forests in mitigating UHIs.The results of the study provide a basis for decision making on forest layout and ecological conservation strategies in Yan'an City. Figure 2 illustrates the flowchat of this paper.

MODIS Land Products
MODIS is suitable for mesoscale urban heat island studies.A large volume of UHIrelated research uses it because of its free access, easy acquisition, and wide coverage [27][28][29][30][31]. MODIS is equipped with two types of sensors, Terra and Aqua.The Aqua sensor has a temporal resolution of 8 days and a spatial resolution of 1000 m, with overpass times at approximately 1:30 p.m. and 1:30 a.m.Compared to the Terra sensor, the daytime temperatures captured by Aqua are very close to the highest temperatures of the day, making it ideal for extracting heat island regions [31,32].This study uses the L3 MODIS LST product (MYD11A1) and the L3 MODIS NDVI product (MYD13Q1) from 2017 to 2022.Both are located within the same coordinate system (WGS 84) but they have different spatial resolutions.It is necessary to resample the LST product through spatial interpolation methods to match the spatial resolution of the NDVI product.Table 1 gives basic information about the MODIS data used in this research.

Data Sources 2.3.1. MODIS Land Products
MODIS is suitable for mesoscale urban heat island studies.A large volume of UHI-related research uses it because of its free access, easy acquisition, and wide coverage [27][28][29][30][31]. MODIS is equipped with two types of sensors, Terra and Aqua.The Aqua sensor has a temporal resolution of 8 days and a spatial resolution of 1000 m, with overpass times at approximately 1:30 p.m. and 1:30 a.m.Compared to the Terra sensor, the daytime temperatures captured by Aqua are very close to the highest temperatures of the day, making it ideal for extracting heat island regions [31,32].This study uses the L3 MODIS LST product (MYD11A1) and the L3 MODIS NDVI product (MYD13Q1) from 2017 to 2022.Both are located within the same coordinate system (WGS 84) but they have different spatial resolutions.It is necessary to resample the LST product through spatial interpolation methods to match the spatial resolution of the NDVI product.Table 1 gives basic information about the MODIS data used in this research.This article utilized the GEE platform to retrieve the MYD11A1 and MYD13Q1 datasets, selecting the LST and NDVI products from June to August (summer in Yan'an) for annual compounding, mosaicking, clipping, and temperature conversion before downloading.The specific steps are as follows: (1) Load the corresponding date's LST/NDVI dataset in the GEE platform.(2) Import the research area vector file and clip the imagery.(3) For the LST dataset, it is also necessary to convert the temperature values to degrees Celsius.The mean composite method ensures the spatial integrity of the data.The urban heat island effect is strongest and the vegetation growth condition is optimal in summer.This meets the requirements of this article to quantify the interannual variation of the urban heat island and to understand its association with vegetation.Figure 3 shows that June-August is the period of the year with the highest temperatures in the city of Yan'an, and it is summer in Yan'an.

China Land Cover Dataset (CLCD)
The CLCD used in this study [33] is based on a comprehensive assembly of 335,709 Landsat images on the GEE.This has led to the creation of the first Landsat-derived annual CLCD.The training samples were collected by combining stable samples extracted from the CLCD with visual interpretation samples from satellite time-series data, Google Earth, and Google Maps imagery.Several temporal indicators were constructed using all available Landsat satellite data, and these were fed into a random forest classifier to obtain classification results.To enhance the temporal and spatial consistency of the CLCD, a post-processing method incorporating spatiotemporal filtering and logical reasoning was proposed.This dataset was employed to delineate the land use types within Yan'an City by downloading six phases of land use type coverage data for the years 2017 to 2022.The land use types of forests, croplands, and grasslands in Yan'an City were extracted using the ArcGIS 10.8 software for mosaicking and clipping and then resampled to a 250 m resolution to ensure consistency with the spatial resolutions of NDVI and LST data.
Celsius.The mean composite method ensures the spatial integrity of the data.The urban heat island effect is strongest and the vegetation growth condition is optimal in summer.This meets the requirements of this article to quantify the interannual variation of the urban heat island and to understand its association with vegetation.Figure 3 shows that June-August is the period of the year with the highest temperatures in the city of Yanʹan, and it is summer in Yanʹan.

China Land Cover Dataset (CLCD)
The CLCD used in this study [33] is based on a comprehensive assembly of 335,709 Landsat images on the GEE.This has led to the creation of the first Landsat-derived annual CLCD.The training samples were collected by combining stable samples extracted from the CLCD with visual interpretation samples from satellite time-series data, Google Earth, and Google Maps imagery.Several temporal indicators were constructed using all available Landsat satellite data, and these were fed into a random forest classifier to obtain classification results.To enhance the temporal and spatial consistency of the CLCD, a postprocessing method incorporating spatiotemporal filtering and logical reasoning was proposed.This dataset was employed to delineate the land use types within Yan'an City by According to existing research, the peak values of NDVI occur during the summer each year, representing the maximum level of vegetation cover throughout the year [34].This study extracts NDVI data from the peak growing season for trend analysis, rather than aggregate NDVI over an entire year.A period from June to August is selected annually from 2017 to 2022, and an average value synthesis is performed for the NDVI of each pixel to determine the maximum NDVI value of Yan'an City within a year.This timeframe minimizes the impact of seasonal changes, as well as climatic and weather disturbances.During the summer, vegetation usually experiences its most vigorous growth; extended daylight hours and suitable temperature conditions create the ideal natural environment for plant growth, leading to the highest biomass and chlorophyll activity and, consequently, the highest annual NDVI values.The larger weather variations in spring and autumn may affect the precision of NDVI, while the relatively stable climate of the summer favors the acquisition of accurate data, thus more accurately reflecting the health status and coverage of vegetation.

Extraction of UHI Areas and UHI Intensity Calculation (1) Calculation of LST
The MODIS MYD11A1 dataset provides estimates of LST, which is particularly applicable to the study of the UHI effect.In this case, the ground temperature data in the MODIS product are mathematically linearly transformed before storage and need to be calculated and processed by Equation ( 2) to obtain the actual value ( • C): where LST is the land surface temperature, DN is the image pixel value, a is the scaling factor, and b is the offset.In this study, the scaling factor is 0.02 and the offset is −273.15.
The data thus processed may have negative values, and all negative values need to be processed as Nodata in ArcGIS 10.8 to avoid subsequent result errors.
(2) Extraction of UHI Region In this study, the mean-standard deviation method was selected to extract the UHI region [35].The study area of this paper spans five years, and the surface temperatures between different years are not comparable.It is necessary to consider the mean and standard deviation of different years comprehensively and utilize the combination of surface temperature mean and different standard deviation multiples to delineate the surface heat field.The mean value and standard deviation can reflect the variation of temperature in different features, with a more reasonable hierarchical structure and better ability to express details.This can avoid the error caused by the temporal phase difference to a certain extent, and the temperature field in the study area is divided into 4 levels using Table 2 as the standard.The high-temperature area is extracted as the UHI region in Yan'an City.surface temperature mean and different standard deviation multiples to delineate the surface heat field.The mean value and standard deviation can reflect the variation of temperature in different features, with a more reasonable hierarchical structure and better ability to express details.This can avoid the error caused by the temporal phase difference to a certain extent, and the temperature field in the study area is divided into 4 levels using Table 2 as the standard.The high-temperature area is extracted as the UHI region in Yan'an City.

Temperature Field Level Division Criteria Colors
High (3) Calculation of surface UHI intensity (SUHII) The high-temperature area is designated as the UHI region.The UHI intensity for Yan'an City from 2017 to 2022 is calculated as Equation ( 3): where SUHII represents the UHI intensity, Tu is the annual average LST of the UHI area, and Tr is the annual average temperature of other areas.surface temperature mean and different standard deviation multiples to delineate the surface heat field.The mean value and standard deviation can reflect the variation of temperature in different features, with a more reasonable hierarchical structure and better ability to express details.This can avoid the error caused by the temporal phase difference to a certain extent, and the temperature field in the study area is divided into 4 levels using Table 2 as the standard.The high-temperature area is extracted as the UHI region in Yan'an City.

Temperature Field Level Division Criteria Colors
High (3) Calculation of surface UHI intensity (SUHII) The high-temperature area is designated as the UHI region.The UHI intensity for Yan'an City from 2017 to 2022 is calculated as Equation (3):

SUHII Tu Tr
(3) where SUHII represents the UHI intensity, Tu is the annual average LST of the UHI area, and Tr is the annual average temperature of other areas.

Statistical Analysis
(1) Correlation Analysis The main purpose of correlation analysis is to assess whether there is a linear relationship between two variables and the strength of the relationship.The correlation analysis in this study was conducted in three main steps: surface temperature mean and different standard deviation multiples to delineate the surface heat field.The mean value and standard deviation can reflect the variation of temperature in different features, with a more reasonable hierarchical structure and better ability to express details.This can avoid the error caused by the temporal phase difference to a certain extent, and the temperature field in the study area is divided into 4 levels using Table 2 as the standard.The high-temperature area is extracted as the UHI region in Yan'an City.

Temperature Field Level Division Criteria Colors
High (3) Calculation of surface UHI intensity (SUHII) The high-temperature area is designated as the UHI region.The UHI intensity for Yan'an City from 2017 to 2022 is calculated as Equation (3):

SUHII Tu Tr
(3) where SUHII represents the UHI intensity, Tu is the annual average LST of the UHI area, and Tr is the annual average temperature of other areas.

Statistical Analysis (1) Correlation Analysis
The main purpose of correlation analysis is to assess whether there is a linear relationship between two variables and the strength of the relationship.The correlation analysis in this study was conducted in three main steps: First, "Create Random Points" was used in ArcGIS 10.8.One thousand different ran- surface temperature mean and different standard deviation multiples to delineate the surface heat field.The mean value and standard deviation can reflect the variation of temperature in different features, with a more reasonable hierarchical structure and better ability to express details.This can avoid the error caused by the temporal phase difference to a certain extent, and the temperature field in the study area is divided into 4 levels using Table 2 as the standard.The high-temperature area is extracted as the UHI region in Yan'an City.

Temperature Field Level Division Criteria Colors
High (3) Calculation of surface UHI intensity (SUHII) The high-temperature area is designated as the UHI region.The UHI intensity for Yan'an City from 2017 to 2022 is calculated as Equation (3):

SUHII Tu Tr
(3) where SUHII represents the UHI intensity, Tu is the annual average LST of the UHI area, and Tr is the annual average temperature of other areas.

Statistical Analysis (1) Correlation Analysis
The main purpose of correlation analysis is to assess whether there is a linear relationship between two variables and the strength of the relationship.The correlation analysis in this study was conducted in three main steps: First, "Create Random Points" was used in ArcGIS 10.8.One thousand different random points (six thousand points in total) were extracted from the NDVI and LST images for each year.Second, the ArcGIS 10.8 tool "Extract Multiple Values to Points" was used.
where Ts is the surface temperature; µ is the mean; and std is the standard deviation.
(3) Calculation of surface UHI intensity (SUHII) The high-temperature area is designated as the UHI region.The UHI intensity for Yan'an City from 2017 to 2022 is calculated as Equation (3): where SUHII represents the UHI intensity, Tu is the annual average LST of the UHI area, and Tr is the annual average temperature of other areas.

Statistical Analysis (1) Correlation Analysis
The main purpose of correlation analysis is to assess whether there is a linear relationship between two variables and the strength of the relationship.The correlation analysis in this study was conducted in three main steps: First, "Create Random Points" was used in ArcGIS 10.8.One thousand different random points (six thousand points in total) were extracted from the NDVI and LST images for each year.Second, the ArcGIS 10.8 tool "Extract Multiple Values to Points" was used.Based on the year as the grouping basis, we divided them into six groups and assigned the values calculated by NDVI and LST to the points and exported the results.Thirdly, the data with anomalies or null values were excluded from the results.Based on previous related studies [34], it can be assumed that there is a linear relationship between NDVI and LST.The Pearson correlation coefficient was used to measure the linear relationship between the variables.The value of this coefficient ranges from −1 to 1, where a positive value indicates positive correlation, a negative value indicates negative correlation, while 0 indicates no linear relationship, and a larger value indicates a more significant correlation, which is calculated as in Equation ( 4): where, r is the Pearson coefficient, n is the number of samples, i is the year ordering, X i is the standardized score of NDVI samples and is the sample mean, Y i is the standardized score of LST samples and is the sample mean.In this study, n = 6, 2017 data i = 1, plus 1 for each year, and 2022 data i = 6.
(2) Overlay Analysis Overlay analysis is often used to study the relationship between different data layers in geographical space.By combining multiple geographic data layers, the spatial association and influence relationship between various geographical elements can be obtained.In this study, vegetation type and temperature field information were overlaid to analyze changes in the temperature fields of various types of vegetation.

Interannual Variation of NDVI and Vegetation Coverage in Yan'an City
The NDVI is very sensitive to the biophysical characteristics of vegetation.It objectively discloses critical information like vegetation cover and growth conditions on a rather extensive spatiotemporal scale.In this study, we have extracted the NDVI value for Yan'an City from MODIS terrestrial products.NDVI values range from −1 to 1, with water bodies exhibiting a negative value.Bare land/wasteland and developed areas have lesser values and are collectively referred to as non-vegetation areas (NDVI ≤ 0.2) [36,37].Areas with an NDVI value greater than 0.2 are defined as vegetation areas.Substantial positive values usually symbolize healthy green vegetation.In this paper, we have integrated data of land use and cover type, defining areas with NDVI ≥ 0.4 as zones with healthy vegetation growth.Based on the above rules, reclassification was conducted and the area and ratio of each category were statistically analyzed.The result indicates that the overall vegetation coverage in Yan'an City is exceptionally high, surpassing 95%.Except for 2017, the proportion of regions with good vegetation growth has always exceeded 60%.Table 3 presents the areas of each value domain.
Spatial distribution maps of different numerical levels in Yan'an City from 2017 to 2022 were created by ArcGIS 10.8.As can be seen from Figure 4, the spatial distribution of vegetation in Yan'an City varies prominently between the six years, overall displaying a "higher in the south and lower in the north" pattern.Well-growing areas (represented in olive green and dark green) have been gradually spreading from the south to the north of Yan'an City, especially with improved vegetation growth around Baota District.Non-vegetation parts of Yan'an City are concentrated in Baota District and the south and east edges of the city.The distribution of built-up areas follows a pattern of "more in the north and less in the south", which is opposite to the distribution of vegetation (see Figures 1b and 4).This is consistent with the spatial and temporal variation rule of the Loess Plateau, which presents a pattern of being higher in the south than the north and gradually decreasing from the south to the north.It coincides with the shift of the ecological farming retreat's focus from the northeast to the southwest in Yan'an City [38,39].The northern part of Yan'an City is characterized by loess ridges and gullies, with a dry climate and scarce rainfall, which is not conducive to vegetation growth.In contrast, the southern part consists of low mountains, hills, and plains with a moist climate that is suitable for vegetation growth [40].As a result, the vegetation cover and quality in Huangling County and Huanglong County in the southern area have always been at a higher level.Meanwhile, due to the impact of ecological restoration projects, the vegetation cover levels and quality within the northern part of Yan'an City, including Zhidan County, Zichang County, and Yanchang County, have been steadily increasing in recent years [41].
In this study, the NDVI values ranged from −0.13 to 0.78 in 2017, −0.11 to 0.82 in 2018, −0.15 to 0.82 in 2019, −0.15 to 0.81 in 2020, −0.12 to 0.83 in 2021, and −0.06 to 0.82 in 2022.Peak values exceeded 0.8 except in 2017, with the average NDVI values increasing from 0.50 to 0.67 and the minimum values increasing from −0.13 to −0.06, as shown in Figure 5.This indicates overall good vegetation growth during the summer in Yan'an City, with improvements in vegetation coverage, growth conditions, and greenness.After analyzing remote sensing data, it was discovered that peak NDVI values are entirely located within the extensive southern forests along the city's borders as depicted in Figure 4.This demonstrates that vegetation growth conditions are superior in the southern part of Yan'an City The northern part of Yan'an City is characterized by loess ridges and gullies, with a dry climate and scarce rainfall, which is not conducive to vegetation growth.In contrast, the southern part consists of low mountains, hills, and plains with a moist climate that is suitable for vegetation growth [40].As a result, the vegetation cover and quality in Huangling County and Huanglong County in the southern area have always been at a higher level.Meanwhile, due to the impact of ecological restoration projects, the vegetation cover levels and quality within the northern part of Yan'an City, including Zhidan County, Zichang County, and Yanchang County, have been steadily increasing in recent years [41].
In this study, the NDVI values ranged from −0.13 to 0.78 in 2017, −0.11 to 0.82 in 2018, −0.15 to 0.82 in 2019, −0.15 to 0.81 in 2020, −0.12 to 0.83 in 2021, and −0.06 to 0.82 in 2022.Peak values exceeded 0.8 except in 2017, with the average NDVI values increasing from 0.50 to 0.67 and the minimum values increasing from −0.13 to −0.06, as shown in Figure 5.This indicates overall good vegetation growth during the summer in Yan'an City, with improvements in vegetation coverage, growth conditions, and greenness.After analyzing remote sensing data, it was discovered that peak NDVI values are entirely located within the extensive southern forests along the city's borders as depicted in Figure 4.This demonstrates that vegetation growth conditions are superior in the southern part of Yan'an City compared to the north, aligning with the general distribution pattern of more vegetation cover in the south than the north.Many forests are concentrated within the southern administrative regions of Fuxian County, Huangling County, and Huanglong County, far away from the dense urban construction areas like Baota District and Ansai County.Forest coverage is less within the central urban area, with nearby forest zones being fragmented and distributed among forest parks.According to official statistics, from 2017 to 2022, Yan'an City completed afforestation over an area of 3855.10 km 2 .The afforestation capacity has been consistently high (Table 4) with forest coverage area experiencing a dramatic increase.Combining CLCD and NDVI data, forest, grassland, and cultivated land areas in Yan'an City from 2017 to 2022 were extracted (CLCD represents forest, grassland, and cultivated land, and NDVI ≥ 0.2) [42][43][44].The proportion of each type of vegetation is shown in Figure 6.Forests are the land type with the most coverage area in Yan'an City, with a coverage rate that shows an increasing trend year by year, growing from 44.77% to 51.69%, an increase rate of 6.92% over six years.Cultivated land exhibits a slight upward trend, increasing from 8.30% to 8.49%.The proportion of grassland cover increased from 38.14% to 39.94%.Many forests are concentrated within the southern administrative regions of Fuxian County, Huangling County, and Huanglong County, far away from the dense urban construction areas like Baota District and Ansai County.Forest coverage is less within the central urban area, with nearby forest zones being fragmented and distributed among forest parks.According to official statistics, from 2017 to 2022, Yan'an City completed afforestation over an area of 3855.10 km 2 .The afforestation capacity has been consistently high (Table 4) with forest coverage area experiencing a dramatic increase.Combining CLCD and NDVI data, forest, grassland, and cultivated land areas in Yan'an City from 2017 to 2022 were extracted (CLCD represents forest, grassland, and cultivated land, and NDVI ≥ 0.2) [42][43][44].The proportion of each type of vegetation is shown in Figure 6.Forests are the land type with the most coverage area in Yan'an City, with a coverage rate that shows an increasing trend year by year, growing from 44.77% to 51.69%, an increase rate of 6.92% over six years.Cultivated land exhibits a slight upward trend, increasing from 8.30% to 8.49%.The proportion of grassland cover increased from 38.14% to 39.94%.The vegetation coverage ratio increased from 84.45% to 93.54%, with a growth rate of 19.08%/6a in Yan'an City.This indicates that Yan'an City had a relatively high level of vegetation coverage originally, and ecological restoration projects were one of the driving factors for the growth of various types of vegetation [45,46].

Interannual Variation of Temperature Field and UHI in Yan'an City
The temperature field in Yan'an City was divided into four temperature field classes of high-temperature zone, sub-high-temperature zone, medium-temperature zone, and low-temperature zone based on the mean-standard deviation method.Table 5 demonstrates the 2017-2022 classification criteria for the summer temperature field in Yan'an City.Between 2017 and 2021, the average surface temperature during the summer in Yan'an City fluctuated and decreased from 32.82 °C to 31.27 °C.The highest average temperature in this six-year span was recorded in 2021 at 32.96 °C, with an average annual change of −1.1%.From the perspective of standard deviation, during the summers from 2017 to 2021, the standard deviation of Yan'an City's surface temperature fluctuated within the range of 2.46-3.45,indicating that the gap between high and low temperatures was relatively stable.As indicated in Table 6, the average surface temperature and standard deviation declined in 2018 and 2022.According to the statistical data, Yan'an City had a greater afforestation effort in 2017 and 2021 (Table 4), suggesting that such changes may be due to anthropogenic factors.The vegetation coverage ratio increased from 84.45% to 93.54%, with a growth rate of 19.08%/6a in Yan'an City.This indicates that Yan'an City had a relatively high level of vegetation coverage originally, and ecological restoration projects were one of the driving factors for the growth of various types of vegetation [45,46].

Interannual Variation of Temperature Field and UHI in Yan'an City
The temperature field in Yan'an City was divided into four temperature field classes of high-temperature zone, sub-high-temperature zone, medium-temperature zone, and lowtemperature zone based on the mean-standard deviation method.Table 5 demonstrates the 2017-2022 classification criteria for the summer temperature field in Yan'an City.Between 2017 and 2021, the average surface temperature during the summer in Yan'an City fluctuated and decreased from 32.82 • C to 31.27 • C. The highest average temperature in this six-year span was recorded in 2021 at 32.96 • C, with an average annual change of −1.1%.From the perspective of standard deviation, during the summers from 2017 to 2021, the standard deviation of Yan'an City's surface temperature fluctuated within the range of 2.46-3.45,indicating that the gap between high and low temperatures was relatively stable.As indicated in Table 6, the average surface temperature and standard deviation declined in 2018 and 2022.According to the statistical data, Yan'an City had a greater afforestation effort in 2017 and 2021 (Table 4), suggesting that such changes may be due to anthropogenic factors.Table 7 shows the area and proportion of various thermal field grades in Yan'an City from 2017 to 2022.It can be seen that high-temperature zone areas and mediumtemperature zone areas show a fluctuating upward trend.The high-temperature zone areas increased by 1.35% over six years, and the medium-temperature zone areas by 6.59% over six years.At the same time, the proportion of sub-high-temperature areas zone decreased, and the proportion of low-temperature zone areas slightly decreased by 0.94% over six years, with the largest decline of 2.68% occurring in 2018.According to the mean-standard deviation method, high-temperature zone areas were extracted as UHI regions.The UHI area saw a significant increase over six years, peaking at 6603.64 square kilometers in 2018.The average temperature of the UHI overall decreased by 2.04 • C, a larger drop than in other areas (1.48 • C).The UHII showed its highest levels from 2019-2021, exceeding 6 • C, but it significantly decreased by 2022.Using ArcGIS 10.8 software, the LST images over the six years were reclassified for mapping (see Figure 7).From the maps, the spatial changes in the UHI regions of Yan'an City can be observed.Between 2017 and 2022, the UHI regions in Yan'an City were concentrated in the northern and eastern border areas, as well as in a very small part of Baota District and the south.According to the reclassification results, there was a certain degree of reduction in the UHI areas of Yan'an City.The medium-temperature and low-temperature zone areas are mainly distributed in the southwest and cover a large area, which might be related to the aforementioned conclusion (the vegetation-covered areas in Yan'an City continuously spread northward, and there are large areas of dense forests in the southwestern region).might be related to the aforementioned conclusion (the vegetation-covered areas in Yan'an City continuously spread northward, and there are large areas of dense forests in the southwestern region).Further spatial overlay analysis of the LST in Yan'an City from 2017 to 2022 in the summer was conducted to obtain the characteristics and intensity of changes in LST area during the period.As shown in Figure 8, the areas in Yan'an City where the temperature field grade increased cover a vast area.For example, 73.32% of the city shows a trend of transitioning from a lower temperature zone to a higher one, including the forest-dense areas in the southwest, where temperatures are gradually shifting from the low-temperature zone to the medium-temperature and sub-high-temperature zones.Among the six types of transitions from lower to higher temperature zones in the city, the shift from medium-temperature zone to sub-high-temperature zone has the highest proportion, reaching 25.08%, followed by the transition from sub-high-temperature zone to high-temperature zone, at 20.69%.The areas transitioning from higher-temperature zones to lower ones account for only 0.31%, noticeably less than the areas transitioning from lower to higher zones.Areas with no change make up 26.37%, mainly consisting of low-temperature and medium-temperature zones, accounting for 6.87% and 6.29%, respectively.The statistical results indicate that in Yan'an City, the high-temperature zone has the smallest area transitioning to other temperature zones, the medium-temperature zone has the largest area transitioning to other zones, with the smallest area transitioning to lower-temperature zones and the largest area changing to the sub-high-temperature zone.Further spatial overlay analysis of the LST in Yan'an City from 2017 to 2022 in the summer was conducted to obtain the characteristics and intensity of changes in LST area during the period.As shown in Figure 8, the areas in Yan'an City where the temperature field grade increased cover a vast area.For example, 73.32% of the city shows a trend of transitioning from a lower temperature zone to a higher one, including the forestdense areas in the southwest, where temperatures are gradually shifting from the lowtemperature zone to the medium-temperature and sub-high-temperature zones.Among the six types of transitions from lower to higher temperature zones in the city, the shift from medium-temperature zone to sub-high-temperature zone has the highest proportion, reaching 25.08%, followed by the transition from sub-high-temperature zone to hightemperature zone, at 20.69%.The areas transitioning from higher-temperature zones to lower ones account for only 0.31%, noticeably less than the areas transitioning from lower to higher zones.Areas with no change make up 26.37%, mainly consisting of lowtemperature and medium-temperature zones, accounting for 6.87% and 6.29%, respectively.The statistical results indicate that in Yan'an City, the high-temperature zone has the smallest area transitioning to other temperature zones, the medium-temperature zone has the largest area transitioning to other zones, with the smallest area transitioning to lower-temperature zones and the largest area changing to the sub-high-temperature zone.

Impact of Forests on the UHI in Yan'an City
Based on Figures 4 and 7, the spatial distribution patterns of densely vegetated zones and high-temperature zones exhibit a certain degree of correlation, while there is a certain overlap with the lower-temperature zones.Analyzing the spatial changes in the vegetation corresponding to the temperature field in Yan'an City from 2017 to 2022 (Figure 9), it is observed that the area maintaining a low-temperature forest zone remains unchanged at the highest proportion, accounting for 9.63% of all types of changes.Following this is the transition from a medium-temperature grassland zone to high-temperature grassland zone, accounting for 8.68%.The proportion of temperature fields transitioning from low to high is the lowest in forests, while for grasslands, the transition proportions from low to high and high to low are both high.In cultivated lands, transitions from low to high in temperature fields are the most frequent.These results indicate that grasslands and cultivated lands are more susceptible to changes in temperature field levels due to other factors, while forests exhibit a certain stability and are less susceptible to changes in temperature field levels due to other factors.

Impact of Forests on the UHI in Yan'an City
Based on Figures 4 and 7, the spatial distribution patterns of densely vegetated zones and high-temperature zones exhibit a certain degree of correlation, while there is a certain overlap with the lower-temperature zones.Analyzing the spatial changes in the vegetation corresponding to the temperature field in Yan'an City from 2017 to 2022 (Figure 9), it is observed that the area maintaining a low-temperature forest zone remains unchanged at the highest proportion, accounting for 9.63% of all types of changes.Following this is the transition from a medium-temperature grassland zone to high-temperature grassland zone, accounting for 8.68%.The proportion of temperature fields transitioning from low Forests are the most stable vegetation type in Yan'an City, and the cooling effect of forests is analyzed separately.Buffer zones of 500 m, 1000 m, and 2000 m are established around the forest.The correlation analysis of LST and NDVI is conducted within different ranges.Analysis results (Figure 10) show that the slopes of the regression equations within the three buffer zone ranges are all negative, consistent with the calculation results of the correlation coefficient (r), indicating a certain negative correlation between LST and NDVI, i.e., as LST increases, NDVI tends to decrease.For every 1-unit increase in LST, the corresponding NDVI decreases by an average of about 0.02 units.The data analysis within the 500 m buffer zone (Figure 10a) shows that in 2017, 2019, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.69, −0.67, and −0.71, respectively, and relatively weaker in 2018 (r = −0.43).The data analysis within the 1000 m buffer zone (Figure 10b) shows that in 2017, 2019, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.69, −0.6, and −0.65, respectively, and relatively weaker in 2018 (r = −0.39).The data analysis within the 2000 m buffer zone (Figure 10c) shows that in 2017, 2020, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.41, −0.35, and −0.35, respectively, and relatively weaker in 2018 (r = −0.17).Overall, the negative correlation between NDVI and LST in Yan'an City is significant, indicating that forests have a noticeable inhibitory effect on temperature rise and can mitigate the urban heat island effect to a certain extent.Forests are the most stable vegetation type in Yan'an City, and the cooling effect of forests is analyzed separately.Buffer zones of 500 m, 1000 m, and 2000 m are established around the forest.The correlation analysis of LST and NDVI is conducted within different ranges.Analysis results (Figure 10) show that the slopes of the regression equations within the three buffer zone ranges are all negative, consistent with the calculation results of the correlation coefficient (r), indicating a certain negative correlation between LST and NDVI, i.e., as LST increases, NDVI tends to decrease.For every 1-unit increase in LST, the corresponding NDVI decreases by an average of about 0.02 units.The data analysis within the 500 m buffer zone (Figure 10a) shows that in 2017, 2019, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.69, −0.67, and −0.71, respectively, and relatively weaker in 2018 (r = −0.43).The data analysis within the 1000 m buffer zone (Figure 10b) shows that in 2017, 2019, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.69, −0.6, and −0.65, respectively, and relatively weaker in 2018 (r = −0.39).The data analysis within the 2000 m buffer zone (Figure 10c) shows that in 2017, 2020, and 2021, the negative correlation between LST and NDVI is more significant, with correlation coefficients of −0.41, −0.35, and −0.35, respectively, and relatively weaker in 2018 (r = −0.17).Overall, the negative correlation between NDVI and LST in Yan'an City is significant, indicating that forests have a noticeable inhibitory effect on temperature rise and can mitigate the urban heat island effect to a certain extent.

Discussion
Since the ecological restoration project was carried out in Yan'an City in 1999, the ecological environment has changed dramatically.The studies that selected this area as the study area mostly focused on exploring its land use type changes and soil erosion improvement [47][48][49], and a few studied the impact of its new forest on UHIs, which is one of the innovations of this paper.As the development of megacities is gradually becoming saturated, more attention should be turned to cities with more development potential, and more reasonable urban planning should be carried out for underdeveloped cities with reference to past experiences.According to the conclusion from the results of this paper, the maximum UHI intensity of Yan'an City reaches 6.03 °C from 2017-2022, and it is necessary to prevent further intensification of the UHI.The study of vegetation dynamics and the spatial distribution of heat islands in Yan'an has revealed a high rate of vegetation cover within the city and intensification of the UHI over the past six years.Correlation analysis between NDVI and LST demonstrates that increasing greenery contributes positively to the mitigation of the UHI, and forests are particularly effective.Correlation results within various buffer zones indicate that afforestation in the urban periphery is effective in reducing temperature fluctuations, and the mitigation of vegetation is more pronounced when it is closer to green areas.Scientific greening strategies and methods are used to enhance the ecological quality and living environment of cities.In the field of urban planning, this study suggests that planners should plant more trees in areas seriously affected by the heat island effect and consider the design of green spaces based on the relationship between vegetation and the cooling effect distance, employing scientific greening strategies and methods to enhance the ecological quality and living environment of cities.Moreover, from a global perspective, the case of Yan'an contributes to a broader understanding of the role of green spaces in ecological restoration and their ability to improve environmental quality.This provides new perspectives and solutions for the heat island phenomenon and environmental challenges faced by cities at different development levels around the world.This paper focuses on the mitigation of the UHI effect in Yan'an City by the studied forest, which has been applied in Shanghai and Zhengzhou,

Discussion
Since the ecological restoration project was carried out in Yan'an City in 1999, the ecological environment has changed dramatically.The studies that selected this area as the study area mostly focused on exploring its land use type changes and soil erosion improvement [47][48][49], and a few studied the impact of its new forest on UHIs, which is one of the innovations of this paper.As the development of megacities is gradually becoming saturated, more attention should be turned to cities with more development potential, and more reasonable urban planning should be carried out for underdeveloped cities with reference to past experiences.According to the conclusion from the results of this paper, the maximum UHI intensity of Yan'an City reaches 6.03 • C from 2017-2022, and it is necessary to prevent further intensification of the UHI.The study of vegetation dynamics and the spatial distribution of heat islands in Yan'an has revealed a high rate of vegetation cover within the city and intensification of the UHI over the past six years.Correlation analysis between NDVI and LST demonstrates that increasing greenery contributes positively to the mitigation of the UHI, and forests are particularly effective.Correlation results within various buffer zones indicate that afforestation in the urban periphery is effective in reducing temperature fluctuations, and the mitigation of vegetation is more pronounced when it is closer to green areas.Scientific greening strategies and methods are used to enhance the ecological quality and living environment of cities.In the field of urban planning, this study suggests that planners should plant more trees in areas seriously affected by the heat island effect and consider the design of green spaces based on the relationship between vegetation and the cooling effect distance, employing scientific greening strategies and methods to enhance the ecological quality and living environment of cities.Moreover, from a global perspective, the case of Yan'an contributes to a broader understanding of the role of green spaces in ecological restoration and their ability to improve environmental quality.This provides new perspectives and solutions for the heat island phenomenon and environmental challenges faced by cities at different development levels around the world.This paper focuses on the mitigation of the UHI effect in Yan'an City by the studied forest, which has been applied in Shanghai and Zhengzhou, China and Medan, Indonesia [50][51][52].
There are some other measures to mitigate the UHI effect, which are summarized in the following three categories: 1.
These generally feature sun-tolerant plants [53].Covering some or all of the top floors of buildings with green plants reduces the warming caused by the reflection of sunlight on the buildings.Plant transpiration and soil evaporation cool down the building surface and improve the thermal and urban microclimate of the building, thus mitigating the UHI effect [54].

2.
Altering the surface materials of the building The principle is to reflect or absorb a portion of the thermal radiation.Cooling materials reduce thermal energy and regulate heat gain during convection, conduction, and radiation heat transfer [55].The cooling power is determined by the radiation emitted and absorbed by the material [56].Common materials include nanomaterials, photonic materials, plasma, phase change materials, etc.

3.
Rationalize the urban pattern This is performed in order to better evacuate the heat gathered in the city [57][58][59].
In this paper, we focus on the mitigating effect of vegetation on UHIs, mainly because of the vast area of forest in Yan'an City.Factors such as the construction location, area, and density of forest affect the mitigation effect on the UHI [60,61].Some studies have shown that increasing the green area near residential areas can increase living comfort [62], and lower temperatures can be capitalized on in terms of house prices [24].Based on the findings of this paper, it is innovative to consider the cooling distance of forest and the spatial heterogeneity of UHIs as some of the considerations for urban planning.The limited nature of land resources determines that green space cannot be increased indefinitely, while strengthening the rationalization of the spatial layout of forest is necessary.This helps to reduce the differences in the thermal environment of the city and can make the layout of urban green space more balanced, scientific, and sustainable.
This study has conducted a profound analysis of the mitigating effects of forests, yet it still has some limitations.Firstly, although this paper discusses the cooling potential of artificial forests on UHIs, it does not include other variables that might impact UHIs, such as environmental humidity, altitude, and economic factors.Future research should carry out a more diversified analysis to consider these variables.Secondly, due to constraints in data acquisition, the MODIS temperature product data have a large spatial resolution, which could deviate from the actual surface temperatures.Future studies may utilize more precise charged sources of temperature data to enhance the accuracy of temperature assessments.Moreover, the evaluation of the limiting distance for forest cooling effects is not exhaustive, as the analysis only considered three distance points.Future studies could refine the step size to yield more precise results.Lastly, further analysis could be conducted on canopy cover and vegetation health in future research.

Conclusions
In exploring the spatiotemporal variability of green spaces in Yan'an City, this study used NDVI calculations to effectively reveal changes in vegetation cover and growth conditions.Vegetation zones (NDVI ≥ 0.2) and well-growing zones (NDVI ≥ 0.4) identified in this study exhibit differentiated spatial distribution patterns and reflect the regional characteristic of "more in the south and less in the north", which aligns with the geographical pattern of Yan'an where the south is suitable for vegetation growth and the north is relatively arid.Overall, the vegetation coverage in Yan'an City is very high, exceeding 95%, and the proportion of well-growing vegetation areas is over 60% in all years except 2017.There is a trend of expansion towards the north in both the vegetation-covered areas and well-growing zones annually, with especially noticeable improvements in the vegetation growth in Baota District.Moreover, the coverage of all types of vegetation in Yan'an City has increased, with the most significant growth observed in forested areas.The overall vegetation coverage ratio has risen to 19.08% per year, and ecological restoration projects are one of the main drivers of this growth across various vegetation types.
Temperature field analysis indicates that the high-temperature zones of Yan'an City are mainly concentrated on the northern and eastern edges of the city, including Baota District, with some appearances in the southern areas.The area of UHI zones has expanded over six years, peaking at 6603.64 square kilometers in 2018, with an average growth rate of 1.35% per year.The UHII was highest from 2019 to 2021, exceeding 6 • C with a peak at 6.26 • C, followed by a decrease in 2022.The analysis of temperature field transitions shows that the majority of areas (73.32%) transitioned from cold to hot zones, with very few areas (0.31%) transitioning to even hotter zones, indicating a general trend of rising temperatures during this period.Additionally, the coverage of different temperature zones has also changed, with an increase in medium-temperature zones and a slight decrease in sub-high-temperature and low-temperature zones.The largest transition area from medium-temperature zones to other fields suggests that the medium-temperature zone is the primary source of change to other zones, representing the most unstable area that is easily affected by other factors and should therefore be given extra attention in urban planning.
By analyzing and comparing the spatial distribution of vegetation and UHIs in Yan'an City, a certain correlation in their spatial patterns is evident.Among different types of landscapes, such as forests, cultivated lands, and grasslands, forests exhibit the most stability against temperature field type changes and effectively mitigate the UHI effect.Further correlative analysis reveals a significant negative correlation between NDVI and LST, signifying that vegetation considerably suppresses a rise in temperature.The relationship between LST and NDVI also varies with spatial scale, where the negative correlation is more significant within 500 and 1000 m of forested regions and weaker within a 2000 m range, demonstrating that the cooling effect of forests is more pronounced when closer to them, with an effective range extending beyond 1000 m.This finding has implications for future urban planning-it is recommended to adjust the urban green space planning to enhance forest, which can help mitigate the urban heat island phenomenon.

Figure 1 .
Figure 1.Study Area.(a) Location of the study area in China; (b) Land Use and Land Cover of Yan'an City in 2017.

Figure 1 .
Figure 1.Study Area.(a) Location of the study area in China; (b) Land Use and Land Cover of Yan'an City in 2017.

Figure 2 .
Figure 2. Flowchart of the data used in the article.This figure shows how they were utilized in a step-by-step manner to ultimately arrive at the conclusions.

Figure 2 .
Figure 2. Flowchart of the data used in the article.This figure shows how they were utilized in a step-by-step manner to ultimately arrive at the conclusions.

Figure 3 .
Figure 3. Monthly average maximum and minimum temperatures in Yan'an City.

Figure 3 .
Figure 3. Monthly average maximum and minimum temperatures in Yan'an City.

2. 3
.3.Statistical Data This study uses new afforestation area data from 2017 to 2022 from the Shaanxi Provincial Statistical Yearbook issued by the Shaanxi Provincial Bureau of Statistics, obtained through the official website of the organization (http://tjj.shaanxi.gov.cn/tjsj/ndsj/tjnj/,accessed on 17 March 2024), which contains in its dataset the specific camping area of Yan'an City for the period of 2017-2022 and the camping area under different afforestation methods, which are credible and authoritative.2.4.Processing Methods 2.4.1.Calculate NDVI NDVI is widely used to assess the quantitative value of vegetation cover and is one of the most suitable indices for monitoring crop growth dynamics.The NDVI was calculated by the following Equation (1): NDVI -Temperature Zone Ts > µ + std Sub-High-Temperature Zone µ < Ts ≤ µ + std Medium-Temperature Zone µ − std ≤ Ts B ≤ µ Low-Temperature Zone Ts < µ − stdwhere Ts is the surface temperature; µ is the mean; and std is the standard deviation.
-Temperature Zone Ts > µ + std Sub-High-Temperature Zone µ < Ts ≤ µ + std Medium-Temperature Zone µ − std ≤ Ts B ≤ µ Low-Temperature Zone Ts < µ − stdwhere Ts is the surface temperature; µ is the mean; and std is the standard deviation.
-Temperature Zone Ts > µ + std Sub-High-Temperature Zone µ < Ts ≤ µ + std Medium-Temperature Zone µ − std ≤ Ts B ≤ µ Low-Temperature Zone Ts < µ − stdwhere Ts is the surface temperature; µ is the mean; and std is the standard deviation.

Figure 4 .
Figure 4. Spatial distribution maps of different numerical levels in Yan'an City from 2017 to 2022.

Figure 4 .
Figure 4. Spatial distribution maps of different numerical levels in Yan'an City from 2017 to 2022.

Figure 8 .
Figure 8. Yan'an City Land Surface Temperature Field Transition Types Map.

Figure 8 . 23 Figure 9 .
Figure 8. Yan'an City Land Surface Temperature Field Transition Types Map.

Figure 9 .
Figure 9.The 2017-2022 Spatial Changes in Vegetation Corresponding to Temperature Field Grades in Yan'an City.

Table 1 .
MODIS Land Products and Relevant Properties.

Table 1 .
MODIS Land Products and Relevant Properties.

Table 2 .
Method for Dividing Land Surface Thermal Field Levels Using Mean-Standard Deviation.

Table 2 .
Method for Dividing Land Surface Thermal Field Levels Using Mean-Standard Deviation.

Table 2 .
Method for Dividing Land Surface Thermal Field Levels Using Mean-Standard Deviation.

Table 2 .
Method for Dividing Land Surface Thermal Field Levels Using Mean-Standard Deviation.

Table 2 .
Method for Dividing Land Surface Thermal Field Levels Using Mean-Standard Deviation.

Table 3 .
Proportions and Area of different NDVI value domains in Yan'an City (2017-2022).

Table 5 .
Proportions of Vegetated Area and Non-vegetated Area in Yan'an City (2017-2022).

Table 5 .
Proportions of Vegetated Area and Non-vegetated Area in Yan'an City (2017-2022).

Table 6 .
Criteria for Division of Temperature Zones in Yan'an City (2017-2022).

Table 7 .
Information on temperature fields and UHI in Yan'an in different years (2017-2022).