Response of Ecosystem Service Value to Spatio-Temporal Pattern Evolution of Land Use in Typical Heavy Industry Cities: A Case Study of Taiyuan City, China

: Ecosystem services value (ESV) has been one index of quantitative evaluation for the ecological livability of heavy industry cities in the new era, which is intimately relevant to patterns of spatio-temporal changes in land use. This study aims to reveal the response of ecosystem service value in heavy industrial cities to the spatial-temporal evolution structure of land use and to analyze the cold and hot spots and sensitivity. In this study, Taiyuan was taken as an example


Introduction
Ecosystem service (ES) means the environmental conditions and utilities formed and maintained by an ecosystem that support the survival and development of humans [1]. It could be divided into the supply of resources, environmental regulation, cultural entertainment, production support, and other functions [2,3]. Ecosystem service value (ESV) is a valid way to quantitatively assess the magnitude of ES. ESV is directly driven by land use changes caused by human activities [4]. The area of each type of land use, land uses, and patterns of spatial variation could affect ecological processes such as material cycle, energy conversion, biological production, and biodiversity, thereby affecting ESV changes [5][6][7][8].
Resource-based heavy industrial cities are areas where economic and environmental conflicts gather, and how to take into account the protection of the ecological environment while the urban economy is sustainable is one of the important issues to be considered Therefore, it is urgent to evaluate heavy industry cities' ESV to determine the impact of human positive and negative ecological behaviors on the regional ecological environment.
Taiyuan is located in the center of the coal mines in Shanxi Province, an industrial base with concentrated energy and heavy chemical industry. Intense human activities (mining, afforestation, ecological restoration, etc.) have changed the quantitative structure and spatial pattern of land use in Taiyuan, which have impacted the natural ecological environment. In this study, we relied on the Landsat remote sensing image data of 2003,2008,2013, and 2018 to make land use type maps and applied the intensity analysis method to explain the spatial and temporal change in land use and land type conversion in Taiyuan in the past 16 years. Meanwhile, we revised the "Ecosystem Service Value Per Unit Area Scale" based on climate productivity and characterized the spatial and temporal evolution of ESVs by constructing fishing nets. Moreover, hotspot analysis was combined with a spatial overlay method to explore the response mechanism between LUCC and ESV. This study can provide the scientific basis for land use management, spatial pattern optimization, and ecological governance and protection in heavy industry cities in the new era.

Study Area
Taiyuan (111 • 30 -113 • 09 E, 37 • 27 -38 • 25 N) is located in the central part of Shanxi Province, covering an area of 6988 km 2 ( Figure 1). It governs six municipal districts, three counties, and one county-level city, of which six municipal districts are the main urban areas of Taiyuan City. Fenhe River valley plain is located in the middle and south, surrounded by mountains in the west, north, and east. The overall terrain is high in the north and low in the south, in the shape of a dustpan, with an average elevation of about 800 m. Taiyuan has a continental monsoon climate in the north temperate zone. The climate is dry, with low and concentrated precipitation. The Fenhe River, the mainstream of the Yellow River, and the Hutuo River, the Haihe River system, are the two major water systems in Taiyuan. The Fenhe River runs through the whole city from north to south. Taiyuan is rich in coal and iron, among which coal mines are mainly distributed in Dongshan, Xishan, and the northern part of Loufan. With abundant mineral resources, Taiyuan's heavy industry has developed rapidly, which has promoted the urbanization process of Taiyuan but caused many problems, such as soil erosion, forest vegetation degradation, environmental pollution, etc., at the same time. In recent years, Taiyuan has carried out a series of afforestation, mine ecological restoration, urban greening, and other measures to improve the environment. In 2018, the afforestation area was 17.28 km 2 , and the green coverage rate in the built-up area reached 42.78% [40]. These data show that the environmental situation has gradually improved.

Data Sources and Processing
The remote sensing data were collected from the geospatial data cloud platform (http: //www.gscloud.cn/, accessed on 21 August 2020). In 2003 and 2008, landsat5 TM digital products were selected for remote sensing images, and in 2013 and 2018, Landsat8 OLI digital products were selected for remote sensing images, which had been systematically corrected for radiation, geometry, and elevation. The fourth phase of image data was processed by GIS software, such as band synthesis, mask extraction, and clipping. The combination of supervised classification and unsupervised classification was applied to classify land use in Taiyuan, and the accuracy of classification results was evaluated by the AOI District method. The overall accuracy of phase 4 image interpretation is 0.95, 0.95, 0.94, and 0.90.
The temperature and precipitation data of Taiyuan are from the daily data set of surface climate data in China provided (V3.0) by the National Meteorological Information Center of the China Meteorological Administration (https://data.cma.cn/, accessed on 21 October 2021). For the date and month of partial data vacancy, the data of the same site in the same period in the previous and next two years of the year were applied, and the average value was taken for interpolation.

Data Sources and Processing
The remote sensing data were collected from the geospatial data cloud platform (http://www.gscloud.cn/, accessed on 21 August 2020). In 2003 and 2008, landsat5 TM digital products were selected for remote sensing images, and in 2013 and 2018, Landsat8 OLI digital products were selected for remote sensing images, which had been systematically corrected for radiation, geometry, and elevation. The fourth phase of image data was processed by GIS software, such as band synthesis, mask extraction, and clipping. The combination of supervised classification and unsupervised classification was applied to classify land use in Taiyuan, and the accuracy of classification results was evaluated by the AOI District method. The overall accuracy of phase 4 image interpretation is 0.95, 0.95, 0.94, and 0.90.
The temperature and precipitation data of Taiyuan are from the daily data set of surface climate data in China provided (V3.0) by the National Meteorological Information Center of the China Meteorological Administration (https://data.cma.cn/, accessed on 21 October 2021). For the date and month of partial data vacancy, the data of the same site in the same period in the previous and next two years of the year were applied, and the average value was taken for interpolation.

Intensity Analysis
The method of intensity analysis was first introduced and applied by Aldwaik et al. [24]. Based on the cross matrix, it explains changes in the intensity and area of land use at time points at a certain location from the levels of time, category, and transition. The interval level explains the rate of land use change in different time periods, the category level elucidates the changes in loss and gain of different land use types, and the transition level indicates the inter-conversion between land use categories. The program is

Intensity Analysis
The method of intensity analysis was first introduced and applied by Aldwaik et al. [24]. Based on the cross matrix, it explains changes in the intensity and area of land use at time points at a certain location from the levels of time, category, and transition. The interval level explains the rate of land use change in different time periods, the category level elucidates the changes in loss and gain of different land use types, and the transition level indicates the inter-conversion between land use categories. The program is encapsulated in Excel 2007 software (https://sites.google.com/site/intensityanalysis/, accessed on 15 October 2021).
The C tij variable is used in all three levels of strength analysis and C tij is a number of pixels that transition from category i at a time Y t to category j at a time Y t+1 .
At the interval level, S t (the annual change percentage during each time interval) and U (the uniform annual change percentage) were compared to reflect which time interval the land use change is rapid (S t > U) or slow (S t < U) at all time intervals.
At the category level, Equation (1) is defined as the uniform intensity at each time interval. G tj (the intensity of the annual gain during a time interval [Y t , Y t+1 ]) and L ti The analysis on the transition level reflects whether the transition between different land use types is strong in a certain time interval. Equations (5) and (6) correspond to the transition-in pattern for the category n, and Equations (7) and (8) correspond to the transition-out pattern for the category m. When the annual intensity of transition R tin exceeds the uniform intensity W tn , it means that the increase of category n is due to the occupation of the category i. On the contrary, it means that the increase of category n avoids the occupation of the category i. On the transition-out pattern of the category mQ tmj m, Q tmj and V tm represents the annual intensity of transition and uniform intensity during the time interval [Y t , Y t+1 ], respectively; moreover, the analysis of the transition-out pattern is the same as that on the transition-in pattern.

Ecosystem Services Value Calculation
The size of biomass reflects the differences in ecosystem service functions to some extent. Climate productivity reflects the plant growth rate and the degree of biomass accommodation by the environment under certain climatic and soil conditions and is an important factor affecting the size of biomass, and the differences in regional environmental conditions of the ecosystem itself also directly affect the type and size of the service function [41,42]. Therefore, based on the research of Xie et al. [32], we adopted two factors of temperature and precipitation in calculating climate productivity to conduct the regional revision for the ecological service value equivalent [34,35]. Moreover, the six types of land use in this paper were sorted as each ecosystem type, among which farmland corresponds to farmland, forestland corresponds to forest, grassland corresponds to grassland, and water body represents the average value of wetlands, rivers, and lakes, unused land corresponds to the desert, and construction land was not taken into account [2]. The calculation process is as follows: In this formula, t is the annual average temperature, • C; r is annual average precipitation, mm; L is the annual average evaporation, mm; V is actual annual evaporation, mm; NPP is the climate productivity; S k is the coefficient of climate difference; NPP a is the climate productivity in the study area; NPP q is the national climate productivity; E a is the economic value of unit equivalent factor in the study area, RMB·hm −2 ; E a is the economic value of the national unit equivalent factor, RMB·hm −2 ; ESV is the value of ecosystem services in the study area, RMB; i is the type of land use; n is the number of types; A i is the area of the land type i.
After calculation, the climate difference coefficients S k of Taiyuan (Table 1).

Sensitivity Analysis
In economics, elasticity is a measure of how buyers and sellers react to changes in market conditions [43]. The elasticity in economics was wielded to quantify the sensitivity of regional ESV to each 1% change in the coefficient of service value of different ecosystem types, which is defined as the sensitivity index (CS). The calculation formula is as follows: Among them, VC ik VC ij , VC ik and VC ij are the ecosystem service value coefficient before and after the adjustment of land types k; ESV i and ESV j are before and after the adjustment of the coefficient of value. When CS > 1, it means that ESV is elastic to CS; if CS is less than 1, it expresses that ESV lacks elasticity to CS, which means that the value coefficient changes in the same direction as the ecosystem service value. Meanwhile, if the value coefficient changes by 1%, the change of ecosystem service value will be less than 1%, then the result is credible.

Analysis of Spatio-Temporal Characteristics of Land Use
From 2003 to 2018, the farmland, grassland, water body, and unused land both decreased in Taiyuan, while the area of forestland and construction land increased. Among these lands, the reduction rate of farmland area was the largest, 0.12%, and the increase rate of construction land area was the most prominent, 0.71%.
The construction land of Taiyuan is mainly concentrated in the east. From 2003 to 2018, it continued to expand to the periphery. However, due to the limitation of terrain, its expansion is mainly along the river valley to the south to occupy the area of farmland and water body ( Figure 2). The northern and western regions are mountains and hilly, mainly grassland and woodland. In addition, the forestland area of Xinghualing District had a significant increase from 2003 to 2018, primarily as a result of the greening and forestry construction area of the East and West Mountains, and it belongs to the security control area of ecological environment in the overall land use planning of Taiyuan (2006-2020). The spatial pattern of land use in Taiyuan is consistent with its natural characteristics, such as topography, topography, and hydrology [44].
Land 2022, 11, x FOR PEER REVIEW and water body ( Figure 2). The northern and western regions are mountains and mainly grassland and woodland. In addition, the forestland area of Xinghualing D had a significant increase from 2003 to 2018, primarily as a result of the greening an estry construction area of the East and West Mountains, and it belongs to the se control area of ecological environment in the overall land use planning of Taiyuan ( 2020). The spatial pattern of land use in Taiyuan is consistent with its natural chara tics, such as topography, topography, and hydrology [44].

Analysis of Spatio-Temporal Characteristics of Land Use
The rapid changes in land use were mainly focused on the two periods of 2008 and 2013-2018, with the uniform annual changing areas accounting for 1.95% and of the study area, respectively, which are greater than the average intensity of 1.28

Interval Level
The rapid changes in land use were mainly focused on the two periods of 2008-2013 and 2013-2018, with the uniform annual changing areas accounting for 1.95% and 1.72% of the study area, respectively, which are greater than the average intensity of 1.28%, illustrating that Taiyuan is in a period of rapid development in two intervals ( Figure 3). However, the average annual changing areas from 2003 to 2018 accounted for only 0.18% of the study area, which was a slowly increasing trend. The economic development and urban expansion of Taiyuan were dilatory and in the early stage of development.

Category Level
The category hierarchy further explains the interval level ( Figure 4). From 2003 to 2008, except for unused land (small base of unused land area, large fluctuation in intensity), the loss and gain intensity of each type were slight (Figure 4a). Among them, the loss intensity of unused land and farmland exceeds the average intensity, and farmland is the type of land with the largest loss (right side of Figure 4a), while the gain intensity of construction land and water bodies exceeds the average intensity, but the change is relatively small. Compared with the period from 2003 to 2008 (Figure 4a), the intensity and area of the transformation of each land type in Taiyuan changed dramatically from 2008 to 2018 (Figure 4b,c), and the urbanization process began to accelerate gradually. From 2008 to 2018,

Category Level
The category hierarchy further explains the interval level ( Figure 4). From 2003 to 2008, except for unused land (small base of unused land area, large fluctuation in intensity), the loss and gain intensity of each type were slight (Figure 4a). Among them, the loss intensity of unused land and farmland exceeds the average intensity, and farmland is the type of land with the largest loss (right side of Figure 4a), while the gain intensity of construction land and water bodies exceeds the average intensity, but the change is relatively small.

Category Level
The category hierarchy further explains the interval level ( Figure 4). From 2003 to 2008, except for unused land (small base of unused land area, large fluctuation in intensity), the loss and gain intensity of each type were slight (Figure 4a). Among them, the loss intensity of unused land and farmland exceeds the average intensity, and farmland is the type of land with the largest loss (right side of Figure 4a), while the gain intensity of construction land and water bodies exceeds the average intensity, but the change is relatively small.  (Figure 4a), the intensity and area of the transformation of each land type in Taiyuan changed dramatically from 2008 to 2018 (Figure 4b,c), and the urbanization process began to accelerate gradually. From 2008 to 2018, the gain intensity of construction land exceeded the average change intensity and was   (Figure 4b,c), and the urbanization process began to accelerate gradually. From 2008 to 2018, the gain intensity of construction land exceeded the average change intensity and was greater than the loss intensity, and the area continued to increase. The loss intensity of farmland and grassland was more significant than the gain intensity from 2008 to 2018, and the area continued to decrease. However, the imbalance difference in gain and loss of each type mainly occurred in the periods from 2008 to 2013. From 2013-2018, only two types of land were concentrated, construction land and farmland. The overall changes in the number of land types tended to balance, showing a spatial transfer. From 2003 to 2018, the loss intensity of farmland and unused land and the gain intensity of construction land and water body were in excess of uniform intensity, which fit with the definition of stability [22,45].

Transition Level
Due to the small proportion of unused land and water body, its changes have little impact on the overall pattern of land use, so only farmland, forestland, grassland, and construction land are analyzed at the transition level.
The intensity that farmland transformed into construction land in three periods was above the uniform intensity, and the stable transformation reached the largest, with a transition area of 28.  (Figure 5b), which resulted from the grain-for-green policy, the implementation of forestry ecological construction, and the further consolidation of future achievements. From 2003-2008, the increased area of grassland stemmed from land. In the same period, the intensity of grassland transfer to construction land overtook the average intensity. From 2013-2018, although there was an area conversion between land types, the area remained relatively stable (Figure 5c). In the whole period, construction land was mainly transferred from various land types, and farmland was the main source, but the average reduction intensity in each period gradually enhanced (from 0 to 0.13 to 0.16) (Figure 5d).
The transfer of various land types in Taiyuan mainly occurred from 2008 to 2013, while the conversion of land types from 2003 to 2008 behaved an evident one-way increase (construction land) or decrease (farmland) in a certain land type. Moreover, two-way transitions among forestland, grassland, and farmland were gradually flat from 2013 to 2018 because 2003-2008 is the initial stage of urban development, and urban development relies solely on the increase of construction land. From 2013 to 2018, the urban development of Taiyuan was gradually stable, the construction land area was charged, and the policy of Farmland Protection, ecological environment construction, and the adjustment of land layout optimized the urban spatial pattern step by step. In combination with the spatial comparability and visualization effect of ESV, researchers selected the grid of 1 × 1 km to figure the sum of ESV of each land type in the grid, which is the spatial distribution of ESV ( Figure 6). The low-value zones are mainly located in the Fenhe Valley, with flat terrain. It is the main urban area of Taiyuan, which mainly consists of construction land and farmland. The high-value zones are distributed in mountains and hilly regions in the west and north, and the land types are mainly forestland and grassland. Due to the difference in ESV of various land types, the spatial trend of ESV in Taiyuan is roughly consistent with the distribution of ESV in a single grid, demonstrating that the ESV of grids with mixed land types is transitory in space.  From 2003 to 2008, the minor change in ESV of Taiyuan showed no significant regional spatial difference, and the conversion intensity between land uses was also tiny during this period. In the same term, due to the increase of forestland area in Xinghualing District and Yangqu County, the ESV in Xinghualing District and the east of Yangqu County increased. The forestland area in the two counties increased by 1396.31 km 2 and 649.00 km 2 . Otherwise, in the east of Qingxu County and the center of Gujiao County, there is an obvious transition from the high-value zone to the low-value zone. From 2013 to 2018, the area of construction land in Taiyuan increased by 43.16 km 2 , which lowered its overall ESV. Meanwhile, the In combination with the spatial comparability and visualization effect of ESV, researchers selected the grid of 1 × 1 km to figure the sum of ESV of each land type in the grid, which is the spatial distribution of ESV ( Figure 6). The low-value zones are mainly located in the Fenhe Valley, with flat terrain. It is the main urban area of Taiyuan, which mainly consists of construction land and farmland. The high-value zones are distributed in mountains and hilly regions in the west and north, and the land types are mainly forestland and grassland. Due to the difference in ESV of various land types, the spatial trend of ESV in Taiyuan is roughly consistent with the distribution of ESV in a single grid, demonstrating that the ESV of grids with mixed land types is transitory in space. From 2003 to 2008, the minor change in ESV of Taiyuan showed no significant regional spatial difference, and the conversion intensity between land uses was also tiny during this period. In the same term, due to the increase of forestland area in Xinghualing District and Yangqu County, the ESV in Xinghualing District and the east of Yangqu County increased. The forestland area in the two counties increased by 1396.31 km 2 and 649.00 km 2 . Otherwise, in the east of Qingxu County and the center of Gujiao County, there is an obvious transition from the high-value zone to the low-value zone. From 2013 to 2018, the area of construction land in Taiyuan increased by 43.16 km 2 , which lowered its overall ESV. Meanwhile, the increase in construction land area (303.89 km 2 from 2003

Sensitivity Analysis
We adjusted the ESV coefficient of each land type in Taiyuan by 50% up and down. The sensitivity index of the ESV to the values coefficient of Taiyuan in four years is less than ( Table 2), indicating that the ESV of Taiyuan is inelastic for value change. Furthermore, the research results are reliable. In addition, the CS of farmland, grassland, and water bodies showed a decreasing trend as a whole. Meanwhile, the CS value of forestland gradually increased and was the highest in all years, which was closely associated with the large proportion of forestland and annual growth in the area in Taiyuan. In 2018, the CS value of the forestland in Taiyuan was 0.631, stating that the ESV in the study area increased by 0.631 percentage points with the 1% change in the ESC coefficient of the forestland.

Analysis of ESV Cold and Hot Spots Based on Land Use Transfer
To discuss the response relationship between ESV and land use change in Taiyuan, we overlaid the results of the hotspot analysis of ESVs with layers of land use transfer from 2003 to 2018. The hot and cold spots of ESV in Taiyuan were mostly concentrated in the eastern main urban area and its surrounding areas and expanded to the south and north (Figure 7).  From 2003 to 2008, the hot spots of ESV mainly gathered in the central part of Yangqu County, which was mainly aroused by the transfer of farmland to forestland and grassland. Cold spots were distributed in the main urban area in the east of Taiyuan City, and the amount of farmland transferred to construction land is up to 22.35 km 2 , which led to a significant reduction of the ESV in this area. From 2008 to 2013, the number of cold and hot spots was more, the distribution scope was wider, and the spatial distribution was relatively diffuse. The hot spots were distributed in Xinghualing District and Yangqu County in the east, and the spatial distribution of ESV cold spots was more dispersed. Large areas of cold spots appeared on both sides of the Fenhe River in Gujiao County and in the west of Qingxu County due to the increase of construction land. From 2013 to 2018, the cold and hot spots of ESV were small and scattered, and the patches of the change in land use were relatively broken, mainly due to the transfer of farmland, forestland, and grassland around the construction land. The urban development tended to be gentle and steadily expanded to the periphery (Figure 7).

Discussion
The change in land use reflects the relationship between people and land [46]. It could announce not only spatial allocation and development of land resources but also the From 2003 to 2008, the hot spots of ESV mainly gathered in the central part of Yangqu County, which was mainly aroused by the transfer of farmland to forestland and grassland. Cold spots were distributed in the main urban area in the east of Taiyuan City, and the amount of farmland transferred to construction land is up to 22.35 km 2 , which led to a significant reduction of the ESV in this area. From 2008 to 2013, the number of cold and hot spots was more, the distribution scope was wider, and the spatial distribution was relatively diffuse. The hot spots were distributed in Xinghualing District and Yangqu County in the east, and the spatial distribution of ESV cold spots was more dispersed. Large areas of cold spots appeared on both sides of the Fenhe River in Gujiao County and in the west of Qingxu County due to the increase of construction land. From 2013 to 2018, the cold and hot spots of ESV were small and scattered, and the patches of the change in land use were relatively broken, mainly due to the transfer of farmland, forestland, and grassland around the construction land. The urban development tended to be gentle and steadily expanded to the periphery (Figure 7).

Discussion
The change in land use reflects the relationship between people and land [46]. It could announce not only spatial allocation and development of land resources but also the concentrated expression of the changes in spatial layout and structure during urban development [47]. Economic development, the grain-for-green policy, ecological governance, and resto-ration policies make the dynamic changes in land use in Taiyuan from 2003 to 2018 mainly occur among farmland, forestland, and construction land. For instance, the southern expansion of the farmland reflects the urban spatial planning of transferring from the north to the south in Taiyuan. According to the analysis of land use intensity at the three levels of time, category, and transition, the farmland has been stably transferred to construction land within three periods, indicating that the progress of urbanization is based on the reduction of farmland. In addition, from 2003 to 2013, the main source of forestland was farmland, which further proved the actual effect of the national policies, such as the grain-for-green policy, forestry ecological construction, and the ecological restoration of the mining area. However, the fact that the transition between forestland and farmland was gradually balanced from 2013 to 2018 indicates the equilibrium between economic development and the ecological environment. The focus of ecological environment protection transformed from the increase in the quantity of forestland to the structural balance between the economy and the environment.
Through the comparison between the transfer matrix, dynamic change of land use, and comprehensive index of land degree, the application of intensity analysis in this paper could quantify the impact of human activities and the internal structural change in regional urban land use, which has been affirmed in the test of Sun et al. [48] and Wang et al. [49]. However, the intensity analysis failed to express the intensity at all levels in the spatial dimension, and it is difficult to prove the impact of spatial adjacency on the change in land use. Hence to this, how to solve this problem is the emphasis of future research [50]. Moreover, the intensity analysis also provides the interpretation for the error in the classification of remote sensing images [51,52], which could better precisely its changing process. This interpretation will be the next step of this paper.
According to the report of the Millennium Ecosystem Assessment (MA), the change in land cover is one of the two most crucial direct driving forces leading to changes in terrestrial ecosystem services in the past 50 years [3]. Changes in the type of land use, the pattern of land use, and the intensity of land use will lead to linkage changes in ES [7]. Based on the revision of climate productivity, researchers obtained the ecosystem service of Taiyuan. From 2003-2018, the ESV of Taiyuan showed a declining trend year by year. During this period, the ESV in 2003 and 2008 had little difference, which was consistent with the changing trend of the ESV in the same period calculated by Liu et al. [53] and Song et al. [37], and the sensitivity analysis verified that the change depends on the forest ecosystem. However, it is difficult to conduct a comprehensive assessment of the regional ecosystem service due to the lack of consideration towards the positive and negative effects of construction land and the single climatic factor in the revised indicator. Otherwise, the heterogeneity of the natural environment, complex multi-scale effects of ecosystem services [7], market failure, and price vacancy [54] all make it more difficult to calculate the ESV. Despite that, it is still of great practical significance to explore the characteristics of changes in ecosystem service value over time [53] and the mechanism process [7]. Meanwhile, according to the sensitivity and hotspots analysis, the dynamic change in forestland, farmland and construction land is an important driving force for the ESV change in Taiyuan. The change in land use is a specific manifestation of urban expansion, so clarifying the impact of urbanization (spatial expansion, expansion forms, etc.) on ESV [18,55] and the response relationship between the two [56] is also one of the key directions in the future study.

Conclusions
Heavy industry cities need to optimize the allocation of land use from the response of land use structure to the spatial and temporal evolution of ESV to improve the ecological environment and the quality of livable ecology. Optimizing the spatial pattern of land use and strengthening the construction of the ecological environment has been urgent issues in heavy industry cities. Meanwhile, the yearly decrease in the area of farmland, grassland, and water body should also be valued. In the later stage of the development of the heavy industry cities, the industrial layout is gradually concentrated, and the environmental problems left are prominent after the relocation of industrial and mining land. Through land reclamation, ecological restoration and other measures to optimize the land use structure, so as to achieve the purpose of increasing the value of ecosystem services. For example, the ecological restoration of the East Mountain in Taiyuan City resulted in the transfer of grassland to woodland and thus increased ESV. According to the research, the intensity of land use in Tai-yuan increases first and then decreases with time. Particularly, farmland has become the main source of the stable increase in construction land. Therefore, it is necessary to coordinate the relationship between urban development and food security. Moreover, the construction of regional ecological civilization is also significant. The above changes in the intensity of land and structure directly affect the temporal and spatial evolution of its ESV. The most obvious is that the ESV in Taiyuan is decreasing year by year, with the greatest change from 2008 to 2013. The cold and hot spots are mainly concentrated in the eastern region, expanding to the south and the north. The research results can provide some references for the early warning of farmland and water area control, the optimization and adjustment of land use structure, and the spatial layout of ecological environment restoration measures in heavy industrial cities.

Data Availability Statement:
The data used in this study are from the geospatial data cloud platform (http://www.gscloud.cn/, accessed on 21 August 2020), the National Meteorological Information Center of the China Meteorological Administration (https://data.cma.cn/, accessed on 21 October 2021).

Conflicts of Interest:
The authors declare no conflict of interest.