The Construction of Ecological Security Pattern under Rapid Urbanization in the Loess Plateau: A Case Study of Taiyuan City

: Taiyuan City in the eastern Loess Plateau has experienced severe ecological problems caused by urban expansion. For cities undergoing rapid urbanization, building an ecological security pattern (ESP) is an effective means to improve urban resilience. Here, geographic information systems (GIS) were used to analyze, manipulate, and visualize urban ecological multi-source information and remote sensing (RS) for the history of land use/land-cover (LULC) changes and the structure of the urban ecological system. Four important ecosystem service functions were estimated: soil conservation, habitat quality, water yield, and carbon storage. The minimum cumulative resistance (MCR) model was combined with the circuit theory method to determine the ecological corridors, pinch points, and barrier points. Our results showed that: (1) from 1980 to 2020, Taiyuan’s built-up area showed increased construction land and enhanced landscape fragmentation. The decline in cultivated land was mainly attributed to construction land. During the period from 2000 to 2010, a greater amount of land was changed in Taiyuan than in other periods; (2) The ecosystem service evaluation based on the LULC in 2020 revealed that the central urban area was lower than the other areas; (3) 38 ecological source sites were identiﬁed, accounting for 16% of the total study area. An area of 106 km 2 was allocated to construct 79 ecological corridors. We identiﬁed 31 ecological pinch points and 6 ecological barrier points; (4) an ESP optimization governance model of “two rings, four zones, and nine corridors” was proposed. Our study provides theoretical guidance for sustainable development and ecological design in Taiyuan City and other regions.


Introduction
Human activity directly affects the natural ecosystem of the Earth's land surface through land use/land-cover (LULC) change [1,2].Globally, LULC changes are as old as human activity [3].Recently, the rate of changing differs from that of earlier changes [4].Rapid urbanization has promoted the proliferation of urban space; most agricultural and ecological land areas have already been converted to construction land, and unnatural landscape elements have increased dramatically [5][6][7].As a result, environmental problems, including ecological patch fragmentation, air pollution, biodiversity reduction, water shortages [8,9], and heat islands [10] have affected the landscape pattern and urban ecological security [11][12][13].There is a complex interactive coupling relationship between urbanization and the ecological environment.How to regulate the relationship between urbanization and the ecological environment has become a worldwide scientific problem and a global 2 of 23 strategic problem.In the critical transition period of China's urbanization development from high speed to high quality, it is an urgent need to study the urbanization process and ecological environment response pattern to provide a scientific basis for regional high-quality development and ecological environment protection in the face of national strategic needs.
It is important to build ecological security patterns (ESPs) to ensure ecological security, maintain ecosystem function, and balance the relationship between economic development and ecological, environmental protection [14][15][16].ESPs originate from the coupling method of spatial patterns and ecological processes in the research of landscape, which focuses on ecological patterns of potential landscape [17].Some nodes, patches, corridors, and entire habitats that contribute to the maintenance and control of certain ecological processes are considered key ecological elements.ESPs construction can perform spatial identification, restoration, and reconstruction of these key ecological elements to improve landscape patterns, achieve effective regulation of specific ecological processes, and significantly enhance ecological benefits [18].For example, Chai et al. [19] explored the relationship between ecological security and urbanization in the Wuhan City Circle during 2005-2018 using Landsat-TM/ETM and Landsat-8 images; based on a time series of multi-source RS data, Wang et al. [20] used pressure-state-response (P-S-R) model to construct Beijing ecological security assessment during 1995-2015; using multi-source RS data, Xu et al. [21] analyzed the interactive coupling mechanism between eco-environmental quality and urban development in China's eastern regions.Through the rational allocation of regional natural resources, achieving point, line, and surface multi-category, multi-level spatial element layout optimization effectively regulates and controls the ecological process so as to maintain a sustainable ecological security state.
Currently, the construction of ESPs mainly includes ecological source identification, resistance surface construction, and ecological corridor identification [22].The basis of building an ESPs is to identify ecological sources and extract ecological corridors.Identifying the ecological sources is mainly based on the quantification of ecological importance, and its quantitative indicators are primarily focused on a certain ecological field; thus, the indicators must be improved.The minimum cumulative resistance (MCR) model is one of the technical methods determining the ecological corridor because it can systematically consider the relationship across internal land cover units [23].However, the MCR model is challenging to identify the key nodes in the corridor.With the support of circuit theory, the ecological flow has been compared to the current in physics, trying to simulate the movement of species in complex ecological networks.It has made certain progress in ESP research [22,24,25].Coupling with MCR Model, the circuit theory is used to identify areas that need to be protected and restored in the ecological corridor.It is a new way for the future construction of regional ESPs.
The Loess Plateau is an important part of the ESP of China's "two screens and three belts", and it is also the most prominent ecology of the contradiction between humans and land.This region has been affected by human activities and climate change, with serious ecological damage, and the coordinated development of urbanization and the ecological environment has been widely concerned.In addition, the Loess Plateau is not only the key area and demonstration area of the project of converting farmland to forests but also an important part of the high-quality development of the Yellow River basin [26].More research has been conducted on the impact of vegetation cover changes on ecosystem services [27][28][29], while less research has been conducted on the relationship between urban expansion and ESPs.Taiyuan City is located in the eastern part of the Loess Plateau.As the capital of Shanxi Province, the urban expansion of Taiyuan has severely impacted ESPs.The harsh ecological background and distinctive mountainous topography have influenced urban development trends, leading to the rapid expansion of the central area, disruption of ecological links between the east and west ends, and serious challenges for urban biodiversity conservation.Our study aimed to: (1) analyze the history of LULC changes in Taiyuan City based on RS data from 1980 to 2020; (2) identify ecological source areas based on ecosystem service evaluation results (soil conservation, habitat quality, water yield, and carbon storage); (3) perform a resistance surface based on slope, elevation, LULC, normalized difference vegetation index (NDVI), and habitat quality; and (4) build the ecological corridors using the MCR Model and circuit theory and determining the key areas of ecological security to protect.Our study provides theoretical guidance for sustainable development and ecological design in Taiyuan City and provides a theoretical reference for the study of urban landscape evolution and ecological security pattern optimization.

Study Area
Taiyuan City (111 1) and the eastern part of the Loess Plateau.Its east, west, and north are surrounded by mountains.Taiyuan City is the most important energy and heavy industry base in China, making substantial contributions to the modernization of construction in China with an important strategic position.The Fen River, as the second major tributary of the Yellow River, flows through Taiyuan City [30].There are more tributaries, albeit with a small flow, that is evenly distributed in each administrative district.The climate in Taiyuan City is affected by monsoons.The mean annual precipitation is 450 mm.There are six districts (Yingze, Xinghualing, Jiancaoping, Wanbolin, Jinyuan, and Xiaodian), three counties in the territory (Yangqu, Loufan, and Qingxu), and a county-level city (Gujiao).
disruption of ecological links between the east and west ends, and serious challenges for urban biodiversity conservation.Our study aimed to: (1) analyze the history of LULC changes in Taiyuan City based on RS data from 1980 to 2020; (2) identify ecological source areas based on ecosystem service evaluation results (soil conservation, habitat quality, water yield, and carbon storage); (3) perform a resistance surface based on slope, elevation, LULC, normalized difference vegetation index (NDVI), and habitat quality; and (4) build the ecological corridors using the MCR Model and circuit theory and determining the key areas of ecological security to protect.Our study provides theoretical guidance for sustainable development and ecological design in Taiyuan City and provides a theoretical reference for the study of urban landscape evolution and ecological security pattern optimization.

Study Area
Taiyuan City (111°30′-113°09′E, 37°27′-38°25′N) is distributed in the central part of Shanxi Province (Figure 1) and the eastern part of the Loess Plateau.Its east, west, and north are surrounded by mountains.Taiyuan City is the most important energy and heavy industry base in China, making substantial contributions to the modernization of construction in China with an important strategic position.The Fen River, as the second major tributary of the Yellow River, flows through Taiyuan City [30].There are more tributaries, albeit with a small flow, that is evenly distributed in each administrative district.The climate in Taiyuan City is affected by monsoons.The mean annual precipitation is 450 mm.There are six districts (Yingze, Xinghualing, Jiancaoping, Wanbolin, Jinyuan, and Xiaodian), three counties in the territory (Yangqu, Loufan, and Qingxu), and a countylevel city (Gujiao).

Data Description
The detailed data sources of this study are shown in Table 1.After processing by ArcGIS 10.8, the slope, elevation, and watershed data were obtained.The LULC data (Landsat-5 TM and Landsat-8 OLI) for 1980, 1990, 2000, 2010, and 2020 were classified into five categories (cultivated, woodland, grassland, water, and unused land).NDVI was calculated using Landsat bands in ArcGIS 10.8.

Analytical Steps
Figure 2 illustrates the technical flow of this study.We adopted an integrated workflow incorporating RS data to explore long-term LULC changes in Taiyuan City during 1980-2020, summarized the spatial evolution of LULC, and accurately identified the city's ecological and environmental problems.LULC and ecosystem services were evaluated to identify and determine ecological sources.A comprehensive resistance surface was constructed using GIS.Furthermore, ecological corridors were built using the MCR model and circuit theory to determine the pinch and obstacle areas.Finally, we constructed an ecological security pattern for Taiyuan.

Landscape Transition Matrix
It is composed of rows and columns representing land types for different periods, of which the behavior is Phase I and the column is Phase II in the land use transfer matrix.
where S stands for land area; n represents the number of LULC; i and j represent the LULC before and after the transfer, respectively; and S ij is the percentage of area converted from i to j in the total land area from Phase I to Phase II.

LULC Dynamic Degree
Within the study period, the single LULC dynamic degree (K) showed the amount of change in a particular type of LULC and reflected the speed of change.The higher the value is, the faster the type changes.The K is given by [31,32]: where K represents the dynamic degree of a certain LULC; U b and U a represent the quantity of a certain LULC at the end and beginning from 1980 to 2020, respectively; and T represents the time interval.

Landscape Transition Matrix
It is composed of rows and columns representing land types for different periods, of which the behavior is Phase I and the column is Phase II in the land use transfer matrix.
where S stands for land area; n represents the number of LULC; i and j represent the LULC before and after the transfer, respectively; and Sij is the percentage of area converted from i to j in the total land area from Phase I to Phase II.

LULC Dynamic Degree
Within the study period, the single LULC dynamic degree (K) showed the amount of change in a particular type of LULC and reflected the speed of change.The higher the value is, the faster the type changes.The K is given by [31,32]: LULC dynamics are characterized by a change in the overall LULC over a given period of time, reflecting the changes in regional LULC [33].The formula is: where LC represents the comprehensive degree of LULC dynamics; LU i represents the area of class i land during the initial year; ∆LU i−j represents the absolute value of the area converted from class i to j within the study period; and T represents the time interval.

Soil Conservation Assessment
A modified universal soil loss equation was applied to estimate the soil conservation service function and estimate potential and actual soil erosion.The formula is: where A c stands for soil loss per year; R is the effect of precipitation on soil erosion; K is the soil erodibility factor; LS is the slope length (L) and steepness (S); C is the effect of cropping and management; and P is the support practice factor.We obtained C and P values from [34].
The precipitation erosion factor (R) was calculated by Wischmeier's model.The formula used is as follows: R = 12 ∑ i=1 1.735 × 10 (1.5lg Pa 2 P −0.08188) (5) where P is the mean annual precipitation (mm) and P a is the mean monthly precipitation (mm).
In this study, calculations were performed according to the characteristics of soil particle size composition and organic carbon (SOC) content.The K is given by Equation ( 6), as follows: where SN = 1 − SAN/100; SAN, SIL, CLA, and C represent the composition percentages of sand, silt, clay, and organic carbon in the soil, respectively.
Soil erosion is directly influenced by topography (i.e., LS).The LS factor calculations are given in Equation ( 7) as follows: where ϕ represents the non-cumulative slope length (m); m represents the slope-length exponent; and θ represents the slope.

Habitat Quality Assessment
The habitat quality module considers the distance weights of threat factors (r) when connecting LULC with threats [35,36].According to the following equation, the stresses level (i) is calculated: where x and y are both grids; d is the linear distance; d rmax is the maximum impact distance on the habitat.Habitat degradation degree (D) is as follows: where R represents the number of r; Y r represents the total number of r grid cells; W r is the weight of the r; r y is the number of stress factors; β is the legal accessibility, and S j represents the sensitivity of land use type j, ranging from 0 to 1.The larger value represents the more sensitive value.The habitat quality index (Q xj ) is as follows: where x and j are the grid and LULC, respectively; H is the habitat suitability of LULC; D is the habitat degradation degree; K is the semi-saturation constant; and z is a normalization constant.Q xj range is 0-1, where a higher value indicates a better habitat quality.

Water Yield Assessment
Water yield is not only an important type of ecosystem service but also the basis of many ecological processes and ecological service functions.It affects ecological functions, maintains ecosystem stability, and improves human well-being [37].According to the climate, topography, LULC, and soil type, the water yield module in the InVEST model calculates the amount of water that will be produced in each grid cell within the watershed [38].In each grid cell, precipitation (P) is subtracted from the actual evapotranspiration (AET) to determine the amount of water produced.According to the following equation, the annual water yield Y(x) was calculated: The AET(x)/P(x) value was calculated as follows [39]: ) 16) where PET is the potential evapotranspiration (mm); x represents the plant evapotranspiration coefficient in pixels; ET 0 , which is calculated by Equation ( 18) [40], represents the reference evapotranspiration of pixels; ω(x) is the coefficient of water availability for plants in pixels; AWC(x) is the volume of water content available to plants (mm); P(x) is the annual precipitation in pixels (mm); and Z is the Zhang coefficient, which ranges from 1 to 10 and is used to reflect the seasonal distribution and hydrological characteristics of rainfall [41].When the value was 2.5, the experimental results were more consistent with the surface runoff in the Water Resources Bulletin of Taiyuan City; therefore, the Zhang coefficient of Taiyuan City was determined to be 2.5.PAWC is the available water capacity of the plants obtained from the difference between the wilting point and field capacity [39].
where Ra is the extra-terrestrial solar radiation (MJ•m −2 •day −1 ); T max and T min are the average maximum and minimum daily air temperature ( • C), respectively.

Carbon Storage Assessment
Based on LULC and soil carbon pools, dead organic matter, aboveground and underground biomass, and harvested wood products, carbon storage was simulated under the current scenario in the InVEST model [42].We included four carbon pools owing to the limited number of HWP measurements.As a result, carbon storage across a landscape can be estimated by aggregating the carbon in each pool across different LULC.The total carbon storage (C) is calculated below: where C soil , C dead , C above , and C below are the carbon storage of soil, soil dead, plant aboveground, and belowground, respectively.Considering that the proportion of the dead organic matter carbon pool in the entire evaluation was small, the change in its value did not affect the evaluation result; therefore, it was excluded from this study.

Extraction of Ecological Corridors
In the present study, the ecological corridors were extracted by the MCR model and circuit theory.The Linkage Mapper can identify ecological corridors and key nodes in heterogeneous landscapes based on circuit theory.Linkage Pathways directly create the environmental information necessary for tools such as Euclidean distance to run in ArcGIS software.

Information Extraction from RS Data
LULC involves complex processes and mechanisms driven by natural and human factors [43].This strongly influences the composition and structure of the regional landscape.It is important to understand that these changes directly affect the spatial distribution of ecosystem services and indirectly affect various ecological processes in the landscape.Land use changes can be effectively reflected in landscape patterns by analyzing their changes [44].
From 1980 to 2020, the landscape structure of Taiyuan City experienced the significant spatial change (Figure 3).There was a decrease in grassland, cultivated land, and unused land in Taiyuan but an increase in urban development, water, and woodland.Construction land has increased by 432.05 km 2 ; whereas a total of 344.82 km 2 of cultivated land and 146.1 km 2 of grassland decreased.The proportion of woodland basically remained unchanged.Although the area of unused land in Taiyuan is small, it generally increased and then decreased, with a total reduction of approximately 1.03 km 2 .There were notable numerical changes in construction land, water area, and unused land from 1980 to 1990 (Table 2): (1) the area of construction land increased by 29.73 km 2 , which was mainly transferred from cultivated land and grassland; (2) the total amount of cultivated land decreased by 28.32 km 2 after subtracting the transferred area in 1990; it was mainly transferred to construction land; (3) water area, woodland, and unused land exhibited an increasing trend, with the area increasing by 3.1 km 2 , 0.39 km 2 , and 0.21 km 2 .During this period, Taiyuan City was in a slow development stage and urban development was relatively balanced.From 1990 to 2000, the areas of cultivated land, woodland, and unused land showed decreasing trends, whereas the other land areas showed increasing trends.The specific performance is as follows: (1) the construction land area has increased by 66.0 km 2 ; (2) the cultivated land area continued to decrease.While the woodland area began to decline, with a cumulative decrease of 36.05 km 2 ; and (3) grassland and water area increased slightly, mainly from woodland and cultivated land.At this stage, the urbanization process began to accelerate, the area of cultivated land and ecological land began to decrease significantly, and the contradiction between urban development and ecological balance began.Wood and construction land expanded from 2000 to 2010, while other types of land decreased.The specific performance is as follows: (1) the construction land increased by 268.6 km 2 , with a significant growth rate; (2) 228.65 km 2 of cultivated land decreased, with most of the land being converted to woodland, grassland, and construction land; and 137.92 km 2 of grassland decreased, with most of the land being There were notable numerical changes in construction land, water area, and unused land from 1980 to 1990 (Table 2): (1) the area of construction land increased by 29.73 km 2 , which was mainly transferred from cultivated land and grassland; (2) the total amount of cultivated land decreased by 28.32 km 2 after subtracting the transferred area in 1990; it was mainly transferred to construction land; (3) water area, woodland, and unused land exhibited an increasing trend, with the area increasing by 3.1 km 2 , 0.39 km 2 , and 0.21 km 2 .During this period, Taiyuan City was in a slow development stage and urban development was relatively balanced.From 1990 to 2000, the areas of cultivated land, woodland, and unused land showed decreasing trends, whereas the other land areas showed increasing trends.The specific performance is as follows: (1) the construction land area has increased by 66.0 km 2 ; (2) the cultivated land area continued to decrease.While the woodland area began to decline, with a cumulative decrease of 36.05 km 2 ; and (3) grassland and water area increased slightly, mainly from woodland and cultivated land.At this stage, the urbanization process began to accelerate, the area of cultivated land and ecological land began to decrease significantly, and the contradiction between urban development and ecological balance began.Wood and construction land expanded from 2000 to 2010, while other types of land decreased.The specific performance is as follows: (1) the construction land increased by 268.6 km 2 , with a significant growth rate; (2) 228.65 km 2 of cultivated land decreased, with most of the land being converted to woodland, grassland, and construction land; and 137.92 km 2 of grassland decreased, with most of the land being converted to cultivated land, woodland, and construction sites; (3) woodland area increased significantly by 114.24 km 2 .However, since 2002, in response to the policy of "returning cultivated land to forest", Taiyuan has actively carried out the work of returning cultivated land to forest, and ecological land has been protected to a certain extent.From 2010 to 2020, the areas of grassland and cultivated land decreased, whereas those of woodland, construction land, and water increased.The specific performance was as follows: (1) construction land area continued to increase, and the growth rate was slower than that of the previous period; (2) although cultivated land and grassland were still declining, the trend was slowing down; and (3) woodland and water areas increased.At this stage, Taiyuan entered a stable development state, the urban expansion rate slowed, and the ecological land use increased.

Characteristics of Land Use Dynamic Change
In addition to revealing the change rate of all land types, the single land use dynamic degree also reflects the intensity of regional land use change [45].According to the single land use dynamic degree change (Table 3), the area of construction land changed the most from 1980 to 2020, followed by unused land, cultivated land, and grassland with dynamic degrees of 3.65%, −1.03%, −0.40%, and −0.2%, respectively.From 2000 to 2010, the single land use dynamic degree of all land types changed significantly, and the change rates of construction and woodland were 4.01% and 0.47%, respectively, which were much higher than those of the other periods.Unused land, water area, cultivated land, and grassland showed decreasing trends, with change rates of −6.86%, −2.08%, −1.19%, and −0.77%, respectively.This phenomenon indicates that the period was characterized by rapid economic development and urbanization in Taiyuan.

Identification of Ecological Source
Taiyuan City is surrounded by mountains on three sides, with a high coverage of natural forests and artificial forests.Moreover, it plays a key role in protecting the soil and water loss in the Yellow River basin.Therefore, in this study, four typical ecosystem service functions (soil conservation, water yield, carbon storage, and habitat quality) were selected for analysis.Then, the analytic hierarchy process (AHP) combined with expert scoring was used to assign corresponding weights to the evaluation results.The ecosystem service functions were evaluated based on their weighted importance.Finally, the evaluation results were graded, and high-value areas were selected as the basis for establishing the ecological source.

Soil Conservation, Water Yield, and Carbon Storage Assessment
In 2020, the total amount of soil conservation in Taiyuan was 2.18 × 10 8 t (Figure 4a), the overall soil conservation capacity was poor, and the soil conservation quantity displayed clear spatial differences.The six main city districts of Taiyuan and Qingxu County have relatively flat terrain, high urbanization levels, continuous areas, and high proportions of construction and cultivated land.Therefore, these areas were included in the overall zoning statistics in this study.Gujiao City, Loufan, and Yangqu counties had higher vegetation coverage and larger topographic fluctuations, and zoning statistics should be conducted.Soil conservation in the six main City Districts of Taiyuan and Qingxu County was 0.41 × 10 8 t, accounting for only 19% of the soil conservation in the region (Figure 4d).The total soil conservation in Yangqu and Loufan Counties and Gujiao City was 0.61 × 10 8 t, 0.53 × 10 8 t, and 0.63 × 10 8 t, respectively, accounting for 81% of the total region area.From the perspective of regional average soil conservation, the values in the central urban area of Taiyuan were significantly lower than those in Yangqu County, Gujiao City, and Loufan County, which had better vegetation coverage.Considering that most areas in the Yangqu and Loufan Counties and Gujiao City have high forest and grassland coverage and are less affected by human activities, the soil conservation capacity is correspondingly high.However, the soil conservation ability of built-up areas, with a high impact on human activities and low vegetation cover, is weak.Therefore, to enhance Taiyuan's development planning, it is necessary to further protect the areas to retain high forest and grass coverage, restore and construct green spaces in built-up areas, and achieve regional ecological security by improving regional soil conservation capacity.
In 2020, the total water yield of Taiyuan was 6.47 × 10 8 m 3 , which is close to the annual runoff of 7.15 × 10 8 m 3 in the Water Resources Bulletin of Taiyuan.Therefore, the InVEST model could better assess water yield in Taiyuan (Figure 4b).The total water production of Yangqu, Loufan Counties, and Gujiao City is relatively high, accounting for 66.8% of the total water production area of the city, with water production of 1.5 × 10 8 m 3 , 1.06 × 10 8 m 3 , and 1.76 × 10 8 m 3 , respectively, followed by Qingxu County, which accounts for 10.23% of the city's total water production, and Taiyuan's "six main urban districts", which only account for 23.01%(Figure 4d).Water yield was inversely proportional to vegetation evapotranspiration.Generally, construction land has the least vegetation and the lowest evapotranspiration; therefore, it has a higher water yield.However, most urban runoff flows through sewage pipes [46,47].Although the water yield was high, its water retention ability was weak.In 2020, the total water yield of Taiyuan was 6.47 × 10 8 m 3 , which is close to the annual runoff of 7.15 × 10 8 m 3 in the Water Resources Bulletin of Taiyuan.Therefore, the InVEST model could better assess water yield in Taiyuan (Figure 4b).The total water production of Yangqu, Loufan Counties, and Gujiao City is relatively high, accounting for 66.8% of the total water production area of the city, with water production of 1.5 × 10 8 m 3 , 1.06 × 10 8 m 3 , and 1.76 × 10 8 m 3 , respectively, followed by Qingxu County, which accounts for 10.23% of the city's total water production, and Taiyuan's "six main urban districts," which only account for 23.01%(Figure 4d).Water yield was inversely proportional to vegetation evapotranspiration.Generally, construction land has the least vegetation and the lowest evapotranspiration; therefore, it has a higher water yield.However, most urban runoff flows through sewage pipes [46,47].Although the water yield was high, its water retention ability was weak.
The carbon storage in Taiyuan ranges from 0.14 to 14.37 Mg/ha, with an average of 12.83 Mg/ha (Figure 4c).There is large spatial heterogeneity in carbon storage in Taiyuan.Carbon storage in urban and water areas was generally low, whereas that in cultivated land, woodland, and grassland was high.The total carbon storage of Taiyuan City was 9.81 × 10 7 t and that of Yangqu County was 3.13 × 10 7 t, which is 31.9% of the total carbon storage (Figure 4d).Loufan and Gujiao Counties have 2.60 × 10 7 t and 1.90 × 10 7 t of carbon storage, respectively, accounting for 23.10% and 19.40%, respectively.These three regions are composed of high vegetation coverage, scattered urban space, and a part of cultivated land; therefore, carbon storage is relatively high.Large areas of cultivated land and rural settlements were distributed in Qingxu County.Carbon storage reached 0.8 × 10 7 t, which is 8.16% of the total carbon storage.The order of carbon storage was as follows: Wanbailin District > Jinyuan District > Xiaodian District > Xinghualing District > Yingze District.The carbon storage in Taiyuan ranges from 0.14 to 14.37 Mg/ha, with an average of 12.83 Mg/ha (Figure 4c).There is large spatial heterogeneity in carbon storage in Taiyuan.Carbon storage in urban and water areas was generally low, whereas that in cultivated land, woodland, and grassland was high.The total carbon storage of Taiyuan City was 9.81 × 10 7 t and that of Yangqu County was 3.13 × 10 7 t, which is 31.9% of the total carbon storage (Figure 4d).Loufan and Gujiao Counties have 2.60 × 10 7 t and 1.90 × 10 7 t of carbon storage, respectively, accounting for 23.10% and 19.40%, respectively.These three regions are composed of high vegetation coverage, scattered urban space, and a part of cultivated land; therefore, carbon storage is relatively high.Large areas of cultivated land and rural settlements were distributed in Qingxu County.Carbon storage reached 0.8 × 10 7 t, which is 8.16% of the total carbon storage.The order of carbon storage was as follows: Wanbailin District > Jinyuan District > Xiaodian District > Xinghualing District > Yingze District.

Habitat Quality Assessment
In this study, the biodiversity module of the InVEST model was used to generate habitat quality using landscape pattern distribution and habitat threat factors.Considering the multiple objective conditions of the study area and related literature [48], we selected cultivated land, roads, rural and urban land use, and industrial and mining land as threat factors.Consulting experts investigated the threat source in Taiyuan City using the maximum distance, weight, and attenuation methods, as shown in Table 4. Based on the InVEST model manual and related literature [48], we defined the effects of different threat factors on land covers (Figure 5).Habitat suitability has a significant impact on land cover.Cultivated land and rural land had a lower impact on land cover, and their impact on river channels was the lowest.All threat factors had a significant effect on woodland, other woodlands, and high-cover grassland, but had a lower effect on open woodland, medium-cover grassland, low-cover grassland, other grasslands, and rivers.
In this study, the biodiversity module of the InVEST model was used to generate habitat quality using landscape pattern distribution and habitat threat factors.Considering the multiple objective conditions of the study area and related literature [48], we selected cultivated land, roads, rural and urban land use, and industrial and mining land as threat factors.Consulting experts investigated the threat source in Taiyuan City using the maximum distance, weight, and attenuation methods, as shown in Table 4. Based on the InVEST model manual and related literature [48], we defined the effects of different threat factors on land covers (Figure 5).Habitat suitability has a significant impact on land cover.Cultivated land and rural land had a lower impact on land cover, and their impact on river channels was the lowest.All threat factors had a significant effect on woodland, other woodlands, and high-cover grassland, but had a lower effect on open woodland, medium-cover grassland, low-cover grassland, other grasslands, and rivers.Ecological lands, such as grasslands, woodlands, and water, are generally subjected to human intervention to a lesser degree and can provide long-term stable habitats and home ranges for organisms.Therefore, the habitat suitability was relatively high.The natural breakpoint method in the ArcGIS is the most appropriate grouping of similar values and maximizes the differences between classes [49].In this study, the results of habitat quality assessments in Taiyuan were divided into five categories using the natural breakpoint method: 0.0-0.4,0.4-0.65,0.65-0.75,0.75-0.85,and 0.85-1.0(Figure 6a).The average habitat quality index for Taiyuan in 2020 is 0.56.The overall habitat quality was low, with the proportions of the five habitat quality types being 38.1, 2.5, 20.8, 13.6, and 24.8%, respectively.Areas with low habitat quality accounted for up to 38.1%, mainly concentrated Ecological lands, such as grasslands, woodlands, and water, are generally subjected to human intervention to a lesser degree and can provide long-term stable habitats and home ranges for organisms.Therefore, the habitat suitability was relatively high.The natural breakpoint method in the ArcGIS is the most appropriate grouping of similar values and maximizes the differences between classes [49].In this study, the results of habitat quality assessments in Taiyuan were divided into five categories using the natural breakpoint method: 0.0-0.4,0.4-0.65,0.65-0.75,0.75-0.85,and 0.85-1.0(Figure 6a).The average habitat quality index for Taiyuan in 2020 is 0.56.The overall habitat quality was low, with the proportions of the five habitat quality types being 38.1, 2.5, 20.8, 13.6, and 24.8%, respectively.Areas with low habitat quality accounted for up to 38.1%, mainly concentrated in the plains of the main urban areas of Taiyuan City (Xiaodian, mean: 0.16; Yingze, mean: 0.43; Xinghualing, mean: 0.41; Jiancaoping, mean: 0.33; Wanbolin, mean: 0.46; Jinyuan, mean: 0.41) and Qingxu County (mean: 0.34) (Figure 6b).The land use types in these regions are mainly urban construction, cultivated, industrial and mining land, and rural residential areas, which are strongly disturbed by regional economic development and frequent human social activity, resulting in poor habitat quality in these regions.Areas with high habitat quality, such as Yangqu (mean value = 0.68), Loufan (mean value = 0.64), and Gujiao (mean value = 0.64), were mainly distributed in mountainous areas owing to less human disturbance and higher vegetation cover.mean: 0.41) and Qingxu County (mean: 0.34) (Figure 6b).The land use types in these regions are mainly urban construction, cultivated, industrial and mining land, and rural residential areas, which are strongly disturbed by regional economic development and frequent human social activity, resulting in poor habitat quality in these regions.Areas with high habitat quality, such as Yangqu (mean value = 0.68), Loufan (mean value = 0.64), and Gujiao (mean value = 0.64), were mainly distributed in mountainous areas owing to less human disturbance and higher vegetation cover.

Ecological Source in Taiyuan City
Based on the above results, habitat quality, soil conservation, water yield, and carbon sequestration storage were weighted using ArcGIS 10.8 with weights of 0.3, 0.3, 0.15, and 0.15, respectively.The natural breakpoint method was used to divide the evaluation results of ecosystem services in Taiyuan into five levels (Figure 7a).Regions with the top 25% of the extracted ecosystem services were selected as the evaluation results [50].The areas with high ecosystem service value in Taiyuan were concentrated in the west and east of Yangqu County, north and south of Gujiao City, the southern part of Loufan County, Jinyuan District, and the southern part of Wanbolin District.These areas are generally covered by forests with large terrain fluctuations and low human intervention, thereby providing a good ecosystem service value for the region.In contrast, in the "six districts" of Taiyuan City, most areas of Qingxu County, the middle part of Yangqu County, the middle part of Gujiao City, and the middle part of Loufan County generally had low ecosystem service values.The ecological sources of Taiyuan were identified by removing patches with an area of less than 4 km 2 .There are 38 ecological source areas, with a total area of 1110.2 km 2 , accounting for 16.13% of the area of Taiyuan City (Figure 7b).Woodland and grassland were the main components of the ecological source areas, with a total area of 95%.The results showed that Yangqu County, Gujiao City, and Loufan County contributed 442.71 km 2 , 364.58 km 2 , and 153.54 km 2 , respectively, of ecological sources, accounting for 86.5% of the total ecological sources.

Ecological Source in Taiyuan City
Based on the above results, habitat quality, soil conservation, water yield, and carbon sequestration storage were weighted using ArcGIS 10.8 with weights of 0.3, 0.3, 0.15, and 0.15, respectively.The natural breakpoint method was used to divide the evaluation results of ecosystem services in Taiyuan into five levels (Figure 7a).Regions with the top 25% of the extracted ecosystem services were selected as the evaluation results [50].The areas with high ecosystem service value in Taiyuan were concentrated in the west and east of Yangqu County, north and south of Gujiao City, the southern part of Loufan County, Jinyuan District, and the southern part of Wanbolin District.These areas are generally covered by forests with large terrain fluctuations and low human intervention, thereby providing a good ecosystem service value for the region.In contrast, in the "six districts" of Taiyuan City, most areas of Qingxu County, the middle part of Yangqu County, the middle part of Gujiao City, and the middle part of Loufan County generally had low ecosystem service values.The ecological sources of Taiyuan were identified by removing patches with an area of less than 4 km 2 .There are 38 ecological source areas, with a total area of 1110.2 km 2 , accounting for 16.13% of the area of Taiyuan City (Figure 7b).Woodland and grassland were the main components of the ecological source areas, with a total area of 95%.The results showed that Yangqu County, Gujiao City, and Loufan County contributed 442.71 km 2 , 364.58 km 2 , and 153.54 km 2 , respectively, of ecological sources, accounting for 86.5% of the total ecological sources.

Construction of Resistance Surface
The spatial distribution of regional ecological resistance is defined as the resistance surface [51].Regions with perfect ecological functions and high ecosystem service values

Construction of Resistance Surface
The spatial distribution of regional ecological resistance is defined as the resistance surface [51].Regions with perfect ecological functions and high ecosystem service values have lower resistance to species migration.A reciprocal relationship between habitat quality and ecological resistance was also considered.In addition, land use type, elevation, slope, and NDVI were selected as resistance factors.Scoring standards and specific resistance values were formulated according to the research of professionals and predecessors.Weight was determined using the AHP (Table 5).A comprehensive resistance surface for Taiyuan City was constructed based on the basis of the single-resistance surface evaluation results.As shown in Figure 8, the resistance surface of each urban and rural construction area in Taiyuan City was large.By contrast, the resistance surface was smaller in areas with high relief and low human activity.Compared with mountainous areas, flat terrain is more conducive to animal migration; however, with the continuous increase in impervious pavement areas in the city, it is difficult to form large green patches in urban built-up areas.Thus, the resistance coefficient of the built-up area in Taiyuan remained very high, which increased the difficulty of constructing and restoring the ecological source area in Taiyuan.Compared to built-up areas, the mountainous area of Taiyuan City has become an important area for ecological protection and biodiversity maintenance.Lower ecological resistance also creates more suitable conditions for corridor construction and reduces the economic input and difficulty of corridor construction.

Identification of Ecological Corridors, Nodes, and Obstacle Points
We identified the length of the corridor and the construction areas according to the constructed source and resistance surface.In the width design, considering the local economic development level and the difficulty of corridor construction, the potential corridor was simulated using a threshold of 5000.This study identified 79 corridors with a total length of 626.7 km and an average length of 8.14 km (Figure 9a).Among the corridors, the shortest was 0.51 km, and the longest was 32.83 km.The corridor space was mainly dominated by woodlands and grasslands, and these areas also had low resistance values.The corridor construction area was 109.82 km 2 , of which cultivated land accounted for 8.5 km 2 , forest land accounted for 70.9 km 2 , grassland accounted for 27.4 km 2 , water area accounted for 1.1 km 2 , and construction land accounted for 1.92 km 2 .The spatial layouts of the corridors in Taiyuan show clear differences.Owing to the higher drag coefficient, the corridors in the six districts and Qingxu County are less distributed.There were four long corridors, two horizontal and two verticals, which were traversed by the low-resistance spaces in the city and effectively created a connection between the ecological sources around the "six urban areas".A number of ecological sources of varying sizes constituted the northern part of Yangqu County and the southern part of Gujiao City; therefore, the corridors were generally short.The ecological sources in the southern parts of Loufan and Yangqu Counties were far apart; the area was composed of many ecological patches with smaller areas, and the distribution of each corridor was "cobwebbed".A comprehensive resistance surface for Taiyuan City was constructed based on the basis of the single-resistance surface evaluation results.As shown in Figure 8, the resistance surface of each urban and rural construction area in Taiyuan City was large.By contrast, the resistance surface was smaller in areas with high relief and low human activity.Compared with mountainous areas, flat terrain is more conducive to animal migration; however, with the continuous increase in impervious pavement areas in the city, it is difficult to form large green patches in urban built-up areas.Thus, the resistance coefficient of the built-up area in Taiyuan remained very high, which increased the difficulty of constructing and restoring the ecological source area in Taiyuan.Compared to built-up areas, the mountainous area of Taiyuan City has become an important area for ecological protection and biodiversity maintenance.Lower ecological resistance also creates more suitable conditions for corridor construction and reduces the economic input and difficulty of corridor construction.

Identification of Ecological Corridors, Nodes, and Obstacle Points
We identified the length of the corridor and the construction areas according to the constructed source and resistance surface.In the width design, considering the local economic development level and the difficulty of corridor construction, the potential corridor was simulated using a threshold of 5000.This study identified 79 corridors with a total length of 626.7 km and an average length of 8.14 km (Figure 9a).Among the corridors, the shortest was 0.51 km, and the longest was 32.83 km.The corridor space was mainly dominated by woodlands and grasslands, and these areas also had low resistance values.The corridor construction area was 109.82 km 2 , of which cultivated land accounted for 8.5 km 2 , forest land accounted for 70.9 km 2 , grassland accounted for 27.4 km 2 , water area accounted for 1.1 km 2 , and construction land accounted for 1.92 km 2 .The spatial layouts of the corridors in Taiyuan show clear differences.Owing to the higher drag coefficient, the corridors in the six districts and Qingxu County are less distributed.There were four long corridors, two horizontal and two verticals, which were traversed by the low-resistance spaces in the city and effectively created a connection between the ecological sources around the "six urban areas".A number of ecological sources of varying sizes constituted Ecological lands and corridors were overlaid, and the nodes of a corridor segment that passed through the ecological land were considered ecological nodes.An ecological breakpoint is highly resistant to normal ecological flow.There are 31 pinch-point locations in the city of Taiyuan.Ecological pinch points in the study area were mostly located near corridors, near the junction of ecological sources, and in the middle of the long corridors (Figure 9b).In addition, corridors in the eastern and western regions of Taiyuan exhibited Ecological lands and corridors were overlaid, and the nodes of a corridor segment that passed through the ecological land were considered ecological nodes.An ecological breakpoint is highly resistant to normal ecological flow.There are 31 pinch-point locations in the city of Taiyuan.Ecological pinch points in the study area were mostly located near corridors, near the junction of ecological sources, and in the middle of the long corridors (Figure 9b).In addition, corridors in the eastern and western regions of Taiyuan exhibited a network distribution, and the number of pinch points increased substantially.According to the identification results of obstacle points, the largest area of obstacle points in Taiyuan was distributed in the central city of Taiyuan and north of TISCO Industrial Park.Land use in these areas is complex, mainly comprising industrial, cultivated, and urban construction land.The main urban road spans Taiyuan from east to west, notably affecting biological passage.The second is Yangqu County, located east of Taiyuan City.This area is at the eastern border of Taiyuan, where a large area of cultivated land and rural settlements is located, and several expressway networks interweave, resulting in a large area of obstacles under the influence of human activity.The obstacle points west of Taiyuan are mainly distributed near the second reservoir of the Fen River.The fragile water ecological environment and longitudinal ecological corridor network led to obstacles in this region.The obstacles and pinch points directly affect the connectivity of regional corridors.Therefore, effectively protecting the ecological environment at these points is key to promoting the quality of regional ecological governance and restoration, and it is also an important factor influencing the smooth connectivity, construction, and maintenance of corridors.

Ecological Source Assessment
The scientific and accurate identification of ecological source areas is an important goal for the construction of ecological security patterns and territorial ecological restoration, and the construction of ecological source area identification systems based on local conditions is the core technology.Owing to the different needs of ecological security, the identification system of ecological source areas varies significantly [22].In the present study, we selected soil conservation, habitat quality, water yield, and carbon storage to identify ecological sources.
The Loess Plateau is one of the most serious soil erosions and the most fragile ecological environment in our country.There are significant differences in the functions of soil and water conservation services between different ecosystems in different regions.Moreover, forest and grass vegetation measures are influenced by the warming and drying climate and human activities [52].The primary function of soil conservation is to alleviate the regional soil erosion caused by water erosion.A high level of soil conservation capacity can effectively improve the water conservation capacity in the region, play a vital role in reducing water and soil loss, and enhance wind and sand fixation [53,54].The results show that the soil conservation of different ecosystems in Taiyuan City in 2020 is in the order of cultivated land < grassland < woodland.The woodland ecosystem has a high function of soil and water conservation, which can be realized mainly through the interception and absorption of precipitation by vegetation canopy, ground cover, and root soil layer to reduce its splash and scour on the soil.The grassland ecosystem lacks the interception effect of the forest canopy on rainfall, and land cover is relatively low.Compared with forest and grass ecosystems, the vegetation coverage of the farmland ecosystem is relatively lower; therefore, soil conservation is the lowest.
Urbanization not only has a direct impact on habitat, but also has an indirect impact on land use through the transmission of ecological protection [55].Taiyuan City is located in the eastern part of the Loess Plateau and has a unique geographical location and diverse landform types.The degradation of habitat quality in typical landscape areas was significant.There are many natural landscapes, such as loess beams, tablelands, mounds, hilly landforms, erosional residual tablelands, and serious problems such as soil erosion and desertification.In addition, urban construction and energy and mineral exploitation have greatly changed the spatial pattern composition and structural characteristics of Taiyuan City, resulting in the reduction in regional greenspace area, the increase in landscape heterogeneity, and the degradation of natural ecology.The average habitat quality index for Taiyuan in 2020 was 0.56, which was higher than that of Tianjin City in 2020 (0.2311) [56].The better the habitat quality, the more the ecological flow can be promoted, corresponding to low resistance [57].Based on habitat suitability and threat factor, the habitat quality results in this study were representative to some extent; woodland and high-coverage grassland had the highest habitat quality, whereas bare land and urban land had lower habitat quality.However, habitat quality in Taiyuan has been declining over the past 40 years owing to land use changes.
Water yield can be evaluated quantitatively and visualized to understand trends in the supply of water in ecosystems and determine the impact of human activities on water resources [58][59][60].From the perspective of water yield service, compared with other land use types, the soil infiltration capacity of woodlands is strong, and the water yield capacity is low owing to the effect of surface runoff interception.Construction land has a high-water yield because of its impervious surface and low precipitation infiltration.The results for Shanghai and Beijing-Tianjin-Hebei in China clearly show that construction land is often covered with asphalt, cement, and concrete, creating an impervious layer that reduces the number of penetrations and concentrations, which leads to an increase in water yield [61].
Regions with high carbon storage capacity can play a positive role in promoting ecosystem production and climate regulation as key indicators of ecosystem services; therefore, they are increasingly used in the evaluation of regional ecosystem services [62][63][64].Research has suggested that LULC has a significant impact on the potential of carbon storage and that natural vegetation is positively correlated with carbon storage [65].According to the results of this study, cultivated land, woodland, and grassland are the most important carbon storage.
This study found that ecological lands, such as mountains, grasslands, and water surfaces, have high ecological suitability and are the main areas for ecological corridor construction.Areas with more serious human interference, such as agricultural land, industrial land, and urban construction land, have higher ecological resistance, which affects the connectivity of ecological corridors.Therefore, urban planning and construction should be prioritized.In general, plains are more suitable for biological migration.However, the rapid development of cities on plains poses great resistance to animal migration.To realize the harmonious development of humans and nature and build ecologically friendly and habitable cities, more ecological patches, such as wetland parks and forest botanical gardens, need to be built and preserved in the process of existing urban planning and renovation instead of moving mountains and filling lakes [66].Simultaneously, the restoration of ecological corridors must be accompanied by the avoidance of urban pollutants and urban flood risks [67,68].The safeguarding of biological migration is one of the main functions of ecological corridors.If organisms encounter a highly polluted environment during their migration to the corridor, the pollutants will have a more negative ecological impact; therefore, the discharge of pollutants near ecological corridors must be strictly controlled.In terms of flood prevention, the urban green corridors identified in this study, such as the Fen River corridor, can be coupled with ecological sponge measures for simultaneous construction [69,70].

Construction of Ecological Security Pattern in Taiyuan City
Human activities have a significant impact on ESP, which can easily damage the local ecological landscape.This results in an increase in landscape complexity, the landscape patch shape becomes more complex, and the connectivity is reduced [71].Similar results in arid areas of China can be found in Ganzhou District [72] and the Middle Reaches of the Heihe River [73].Therefore, improving the connectivity and stability of the ecological landscape has a significant effect on the response to disturbances in human activities.In the present study, the regional topographical features and migration habits of biological species were fully considered, and resistance surfaces and potential ecological corridors were constructed to protect biodiversity and promote river wetland environment construction to identify ESP.Taiyuan City's ecological security pattern was built by connecting the ecological corridors and ecological sources.In this study, owing to the large area of impervious pavement in built-up urban and rural areas, urban ecology was hindered as a result of a lack of ecological sources, corridors, and human settlement environments.In recent years, with the relevant departments of Taiyuan City continuing to strengthen the Fen River landscape transformation and restoration, the Fen River has effectively played the role of an urban green corridor, becoming an important place for wildlife migration and the recreation of citizens.Therefore, this study used the Fen River as an independent corridor in the construction of an urban corridor and constructed an ecological security pattern optimization governance model of "two ring roads, four zones, and nine corridors" (Figure 10).The two rings refer to the land corridor around the city of Taiyuan and the waterfront corridor along the bank of the Fen River.Through the construction of these two corridors, an ecological corridor network space connecting the north, south, east, and west was formed, effectively connecting the ecological sources of Taiyuan.The central part of Taiyuan is the primary construction and development area.Urban construction land has a large spatial area and lacks large green areas.In the long run, it will continue to divide the connection between Taiyuan City and ecological land at both ends and affect the migration of animals.It is urgent to build an ecological corridor in the east-west direction to repair the ecological connection between the eastern and western mountains.
"Four zones" refers to the ecological control zone of high coverage vegetation in Yangqu county and Gujiao City in the north of Taiyuan City; Tianlong, Meng Mountain, and Tai Mountain (scenic spots with restricted development zones) in Jinyuan District; the ecological barrier zone in the southwest of Gujiao City and Loufan County; and an ecological conservation zone in the northeast of Yangqu County.A common feature of these areas is that they have large ecological patches and are the main ecological source distribution areas.
The term "nine corridors" refers to the "two vertical and one horizontal" corridors which connect the ecological source areas in the main urban area of Taiyuan City; "one horizontal and three verticals" corridors connecting Yangqu County, Gujiao City, and Loufan County of Taiyuan City; and "one vertical and one horizontal" corridor in the northeast of Yangqu County.The "two vertical corridors and one horizontal" are located The two rings refer to the land corridor around the city of Taiyuan and the waterfront corridor along the bank of the Fen River.Through the construction of these two corridors, an ecological corridor network space connecting the north, south, east, and west was formed, effectively connecting the ecological sources of Taiyuan.The central part of Taiyuan is the primary construction and development area.Urban construction land has a large spatial area and lacks large green areas.In the long run, it will continue to divide the connection between Taiyuan City and ecological land at both ends and affect the migration of animals.It is urgent to build an ecological corridor in the east-west direction to repair the ecological connection between the eastern and western mountains.
"Four zones" refers to the ecological control zone of high coverage vegetation in Yangqu county and Gujiao City in the north of Taiyuan City; Tianlong, Meng Mountain, and Tai Mountain (scenic spots with restricted development zones) in Jinyuan District; the ecological barrier zone in the southwest of Gujiao City and Loufan County; and an ecological conservation zone in the northeast of Yangqu County.A common feature of these areas is that they have large ecological patches and are the main ecological source distribution areas.
The term "nine corridors" refers to the "two vertical and one horizontal" corridors which connect the ecological source areas in the main urban area of Taiyuan City; "one horizontal and three verticals" corridors connecting Yangqu County, Gujiao City, and Loufan County of Taiyuan City; and "one vertical and one horizontal" corridor in the northeast of Yangqu County.The "two vertical corridors and one horizontal" are located in the main urban area, with frequent human activity, resulting in high resistance to corridor construction.Located in the west of Taiyuan City, the "one horizontal and three verticals" corridors connecting Yangqu County, Gujiao City, and Loufan County are composed of small ecological patches.These areas encompass the water source protection area of Taiyuan City, the second reservoir of the Fen River.Therefore, they are more suitable for the living and breeding space of wild animals and plants.Located in northeast Yangqu County, Taiyuan City, the "one vertical and one horizontal" corridor passes through the important "ecological conservation area".Although there was high vegetation coverage in the study area, the ecosystem services value was low.Through the construction of the corridor, the ecological source space was effectively combined, and strict source protection and active corridor construction were implemented to promote regional ecological development and improve ecosystem service value.Collectively, these nine corridors constitute the construction of an ecological network in Taiyuan, which can effectively increase the links between ecological sources and become an important way to achieve sustainable development in the region.

Limitations and Future Research Directions
Firstly, limited by the length of the paper, this study only constructed ESP in 2020.In future studies, changes in ESPs during different periods should be considered.Secondly, the settings of some parameters in the InVEST model are subjective, and more objective methods will be integrated to set the model parameters in the future.Thirdly, in terms of selecting ecological nodes, more practical nodes can be proposed from the perspective of biodiversity protection.Only land use type and terrain were considered in the construction of the resistance surface.In the future, social resistance factors such as population and economy can be included to build a more comprehensive resistance surface.Finally, future studies can be strengthened from the aspects of the diachronic, cross-scale, and driving mechanism of ESPs, and the influencing factors, dynamic processes, and driving mechanisms of the evolution of ESPs.

Conclusions
In this study, we systematically analyzed the history of LULC changes in Taiyuan City from 1980 to 2020.We then assessed four ecosystem services in 2020, including soil conservation, habitat quality, water yield, and carbon storage.Furthermore, we identified ecological corridors, pinch points, and barrier points by combining the MCR model and circuit theory method.Finally, we constructed an optimization model of the ecological security pattern.The main conclusions are as follows: (1) During the period from 1980 to 2020, urban construction land has continuously increased, with a cumulative increase of 432.05 km 2 ; cultivated land area continued to decrease, with a cumulative decrease of 344.82 km 2 ; grassland and water area have also experienced a decrease of 146.1 km 2 and 1.98 km 2 , respectively.Significant LULC changes occurred from 2000 to 2010, which was also a decade of rapid development.
(2) The total amount of soil conservation in Taiyuan was 2.18 × 10 8 t.The total water yield was 6.47 × 10 8 m 3 , which is close to the annual runoff of 7.15 × 10 8 m 3 in the Water Resources Bulletin of Taiyuan.The carbon storage ranges from 0.14 to 14.37 Mg/ha, with an average of 12.83 Mg/ha.The average habitat quality index was 0.56.(3) A total of 38 ecological sources were identified, covering an area of 1124.16 km 2 and accounting for 16% of the total area of Taiyuan City.The ecological sources were mainly woodlands and grasslands.Comprehensive resistance surface analysis shows that areas with high resistance values were mainly concentrated in the central urban area of Taiyuan City and the built-up area of local counties.A total of 79 corridors and 31 ecological "pinch points" were identified in Taiyuan City, which were mainly distributed in the corridors near each ecological source.Six ecological barrier points were also identified, the largest of which was located near the northern corridor of Taiyuan City.There are significant spatial differences in ecological sources in Taiyuan, which are shown as "six districts" and Qingxu county are lower than Yangqu County, Gujiao City, and Loufan County.(4) Integrate the ecological source, corridor, pinch, and obstacle points, and bring the Fenhe River bank zone into the construction of regional ecological security patterns as an important river ecological corridor.The optimal management model of the ecological security pattern of "two rings, four regions, and nine corridors" was constructed.

Figure 1 .
Figure 1.Map of Taiyuan city location and elevation.

Figure 1 .
Figure 1.Map of Taiyuan city location and elevation.

Figure 2 .
Figure 2. Map of the ESPs construction.

Figure 2 .
Figure 2. Map of the ESPs construction.

25 Figure 3 .
Figure 3. Land use change from 1980 to 2020 in Taiyuan City.

Figure 3 .
Figure 3. Land use change from 1980 to 2020 in Taiyuan City.

Figure 4 .
Figure 4. Distribution of (a) soil conservation, (b) water yield, and (c) carbon storage amount in Taiyuan City; (d) amount of soil conservation, water yield, and carbon storage in each administrative region.

Figure 4 .
Figure 4. Distribution of (a) soil conservation, (b) water yield, and (c) carbon storage amount in Taiyuan City; (d) amount of soil conservation, water yield, and carbon storage in each administrative region.

Figure 5 .
Figure 5. Sensitivity of land covers to habitat threat factor.

Figure 5 .
Figure 5. Sensitivity of land covers to habitat threat factor.

Figure 6 .
Figure 6.(a) Habitat quality assessments.(b) Mean habitat quality value in each administrative region.

Figure 6 .
Figure 6.(a) Habitat quality assessments.(b) Mean habitat quality value in each administrative region.

25 Figure 7 .
Figure 7. (a) Ecological source in Taiyuan City.(b) Ecological source areas in each administrative region.

Figure 7 .
Figure 7. (a) Ecological source in Taiyuan City.(b) Ecological source areas in each administrative region.

Figure 9 .
Figure 9. (a) Ecological corridors in Taiyuan City and (b) ecological pinch and barrier point construction in Taiyuan City.

25 Figure 10 .
Figure 10.Optimization model of ecological security pattern in Taiyuan City.

Figure 10 .
Figure 10.Optimization model of ecological security pattern in Taiyuan City.

Table 1 .
The detailed data descriptions.

Table 2 .
Land use transfer matrix of Taiyuan City from 1980 to 2020.

Table 3 .
Single landscape dynamics in Taiyuan City from 1980 to 2020.

Table 5 .
Classification and weight of each resistance factor.