The Ecosystem Services Value Change and Its Driving Forces Responding to Spatio-Temporal Process of Landscape Pattern in the Co-Urbanized Area

The development of urbanization is still expanding on the earth, and the rapid expansion of cities has changed the regional landscape pattern and significantly affected the value of regional ecosystem services in developing countries such as China. Zhongmu County, as the core area of Zhengzhou-Kaifeng integration, studying the temporal and spatial transformations of its landscape pattern and ecosystem service value (ESV) is of great significance to the region’s sustainable development. Based on remote sensing images and socioeconomic data, this study aims to explore the landscape pattern of Zhongmu County from 2005 to 2018 and its impact on ESV. Research methods include an ESV equivalent factor method, landscape pattern index, spatial autocorrelation, and other methods. The results showed that: (1) During the study period, the patch density and shape complexity, landscape diversity, and fragmentation of the overall landscape in the study area continued to increase while landscape connectivity decreased. (2) The total amount of ESV increased by 10.05 million USD; ESV had certain differences in spatial distribution: high-value clusters were mainly located at the boundary of the Yellow River in the north, while low-value clusters had a significant eastward expansion trend. (3) ESV increased significantly in areas where cultivated land was transferred to waters and forests, and ESV in areas where waters transferred to construction land and cultivated land was significantly reduced. (4) ESV had a significant positive spatial correlation with patch density, edge density, mean patch fractal dimension, mean patch size, and the landscape shape index, and a spatially negative correlation with the contagion index and cohesion index. The spatial and temporal changes in landscape patterns and ESV were all mainly driven by the regional development “Zhengzhou-Kaifeng integration” policy. Therefore, the ESV can be improved, and the ecological security of the urban integration area can be guaranteed through policy measures such as optimizing the layout of construction land and adjusting the uniform distribution of green areas by the land-use policy.


Introduction
Ecosystem services (ESs) refer to life support products and services that humans obtain directly or indirectly through the structure, process, and function of ecosystems [1][2][3]. The United Nations Millennium Ecosystem Assessment (MA) divides ES into four categories: support, regulation, supply, and cultural functions [4]. Ouyang et al. [5] found that ES forms the life support system on which human beings depend, such as the generation and maintenance of biodiversity, climate regulation, nutrient storage and cycling, soil renewal, Land 2021, 10, 1043 3 of 21 increased its urbanization development speed. From 2005 to 2018, Zhengkai Avenue and Zhengkai Intercity Railway in Zhongmu County were completed and opened to traffic. A comprehensive transportation system coordinated with the Zhengzhou-Kaifeng integration strategy was initially established. At the same time, the construction land for industrial parks and residential communities is increasing, and the urban development pattern of "three axes, three districts, and multiple clusters" was initially formed. It can be seen that the landscape pattern of Zhongmu County has undergone rapid and drastic changes [26]. At the same time, ES, which is closely related to it, will inevitably be affected. In order to prevent the ES of Zhongmu County from being negatively affected by urbanization, the stable supply and gradual increase in ESV must be ensured, and at the same time the balanced development of ecological environmental protection and the construction of Zhengzhou-Kaifeng integration must be realized.
This study was based on the remote sensing images and socioeconomic data of Zhongmu County from 2005 to 2018, and used multiple methods such as the ecosystem service value equivalent factor method, landscape pattern index, land-use transfer matrix, cold and hot spots analysis and spatial autocorrelation to analyze the spatio-temporal changes in Zhongmu County's landscape pattern and the ESV in the Zhengzhou-Kaifeng integration. The equivalent factor method was used to evaluate the value of the ecosystem in Zhongmu County. The landscape pattern index was used to express the dynamic changes in different landscape pattern characteristics in Zhongmu County, such as landscape composition, fragmentation, and connectivity. The land-use transfer matrix was used to reflect the temporal and spatial changes in the landscape pattern in Zhongmu County from 2005 to 2018. The cold and hot spots analysis was used to explore whether there are statistically significant hot and cold spot clusters in the ESV variation in Zhongmu County, and determine the location of the hot spot clusters and cold spot clusters in Zhongmu County. The spatial autocorrelation analysis method was used to reflect the spatial distribution of ESV itself and the characteristics of landscape patterns in Zhongmu County. The objectives of this study were: (1) in-depth exploration of the dynamic response between landscape pattern and ESV; (2) to explore how to optimize the landscape pattern of urban integration areas; (3) to discuss the dynamics of driving forces and their impact on the development and planning of urban integration areas.

Study Area
Zhongmu County is a county under the jurisdiction of Zhengzhou City, the capital of Henan Province (34 • 26 -34 • 56 N, 113 • 46 -114 • 12 E) [27]. Located in the hinterland of the Central Plains, on the shore of the Yellow River. It is a warm temperate continental monsoon climate with abundant water resources and more than 40 large and small rivers. Zhongmu County is adjacent to the central urban area of Zhengzhou City in the west and the Western New District of Kaifeng City in the east. It is the core area and strategic fulcrum of the Zhengzhou-Kaifeng Integration. At the same time, as the main battlefield of Zhengzhou City's eastward expansion, it is located in the "three districts" superimposed development area of the Central Plains urban agglomeration, Zhengzhou metropolitan area, and the Zhengzhou Airport Economic Comprehensive Experimental Zone. The Zhongmu County Urban-Rural Master Plan (2015)(2016)(2017)(2018)(2019)(2020)(2021)(2022)(2023)(2024)(2025)(2026)(2027)(2028)(2029)(2030) indicates that the administrative area of Zhongmu County will deduct the Baisha group that is classified into the Zhengdong New District, the Jiulong group that is classified into the Economic and Technological Development Zone, and the area outside the Airport Economic Comprehensive Experimental Zone, with a planned area of about 934 km 2 [26]. Since the above areas in Zhongmu County were not classified into other regions in 2005 and 2010, their development had an important relationship with and impacts on the landscape pattern and ecological environment of Zhongmu County. As well as to facilitate the comparison of ESV spatio-temporal characteristics between different years, the research chose to conduct ment had an important relationship with and impacts on the landscape pattern and ecological environment of Zhongmu County. As well as to facilitate the comparison of ESV spatio-temporal characteristics between different years, the research chose to conduct the analysis based on the extent of Zhongmu County in 2005, whose location and satellite images are shown in Figure 1.

Geospatial Data
The landscape classification data and administrative division data of the three phases of Zhongmu County in 2005, 2010 and 2018 are all sourced from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resde.cn, accessed on 24 December 2020). The landscape classification data is a raster map (resolution 30 m) generated by manual visual interpretation of remote sensing images, and the interpretation accuracy is more than 90%. The administrative division data of Zhongmu County was combined and the clipping tool of the ArcGIS 10.2 software was used to clip the 30 m raster map of China's land-use remote sensing monitoring, from which the landscape classification map of Zhongmu County can be obtained. The classification is based on the "China Land Use/Land Cover Remote Sensing Monitoring data classification system", and the specific classification description is shown in Table A1. Using ArcGIS 10.2, the study area was divided into 1 km×1 km fishing nets, and a total of 1511 grids were formed. Calculating the ESV and landscape pattern index in the grid separately can further realize the in-depth exploration of the spatial dynamics of ESV and landscape pattern characteristics and the spatial coupling relationship between the two at a smaller unit scale.

Geospatial Data
The landscape classification data and administrative division data of the three phases of Zhongmu County in 2005, 2010 and 2018 are all sourced from the Resource and Environmental Science Data Center of the Chinese Academy of Sciences (http://www.resde.cn, accessed on 24 December 2020). The landscape classification data is a raster map (resolution 30 m) generated by manual visual interpretation of remote sensing images, and the interpretation accuracy is more than 90%. The administrative division data of Zhongmu County was combined and the clipping tool of the ArcGIS 10.2 software was used to clip the 30 m raster map of China's land-use remote sensing monitoring, from which the landscape classification map of Zhongmu County can be obtained. The classification is based on the "China Land Use/Land Cover Remote Sensing Monitoring data classification system", and the specific classification description is shown in Table A1. Using ArcGIS 10.2, the study area was divided into 1 km×1 km fishing nets, and a total of 1511 grids were formed. Calculating the ESV and landscape pattern index in the grid separately can further realize the in-depth exploration of the spatial dynamics of ESV and landscape pattern characteristics and the spatial coupling relationship between the two at a smaller unit scale. fluctuations of different types of crops on the total value, the study used the main food types in Zhongmu County in 2018: wheat, corn, and rice as basic data. For the time being, the impact of climate on food prices was considered.

Landscape Pattern Index
In order to describe the dynamic changes in different landscape pattern features such as landscape composition, fragmentation, and connectivity in Zhongmu County from 2005 to 2018, and reduce the one-sidedness and limitations of a single index, this study selected ten representative landscape pattern indexes. The descriptions of the indexes are shown in Table A2. The comprehensive analysis of the ten indexes can more comprehensively reflect the characteristics of the landscape pattern of Zhongmu County from 2005 to 2018, and conduct landscape dynamic research. The index calculation was performed in the Fragstats 4.2.1 software.

Land-Use Transfer Matrix
The land-use transfer matrix can quantitatively describe the mutual transfer of landscape types in the initial and final stages of Zhongmu County, and reflect the transfer direction and nature of various landscapes, and then reveal the structural characteristics of multiple landscapes before and after the transfer. It can also effectively reflect the spatio-temporal changes in the landscape pattern in Zhongmu County during a certain period [28,29]. The formula is: where S is the area of various landscapes, n is the number of landscape types, and i and j are the landscape types at the beginning and end of the study, respectively.

Ecosystem Service Value Assessment
The equivalence factor method is the main method to evaluate the value of the ecosystem. Scholars have widely adopted it to quantify ESV because it is easy to use, requires less data, and has highly comparable results [30,31]. The method can evaluate the value of the ecosystem in Zhongmu County. The method was first proposed by Costanza and modified by Xie Gaodi and other scholars to derive the equivalent factor of the terrestrial ecosystem service value in China (Table A3), which is defined as the economic value of one identical factor equal to 1/7 of the market value of the average grain yield in the study area in the current year [32]. The formula is as follows: where Q 1 is the average grain output, F 1 is the average grain price. E x is the ESV per unit area of landscape type x, a 1 is the value equivalent of dry land, and a x is the value equivalent of other landscape types. Since this method does not involve construction land, this study does not consider the ESV of construction land [15]. In addition, according to Formulas (4) and (5), the ESV of different landscape types and different ES functions can be obtained, and the formulas are as follows: Land 2021, 10, 1043 6 of 21 where A k and VC k represent the area and ESV coefficient of landscape type k, respectively, ESV f represents the value of the fth ecosystem service function, and VC f k represents the ESV coefficient of the fth service function of landscape type k [33].

Analysis of Cold and Hot Spots of ESV
The cold and hot spot analysis was used to explore whether there are statistically significant hot and cold spot clusters in the ESV variation in Zhongmu County, and determine the hot spot clusters' location [34]. The hot spot indicates that the ESV of Zhongmu County increased significantly, and the increase is much higher than that of the surrounding area; the cold spot indicates that the ESV of Zhongmu County dropped significantly, and the decrease is far more than that of the surrounding area. The cold and hot spot analysis is based on the Getis-Ord Gi* index, the formula is: where X j is the attribute value of spatial unit j; W ij represents the spatial weight between spatial units i and j (adjacent is 1, non-adjacent is 0); n is the number of spatial units; X is the mean value.

Spatial Autocorrelation of ESV
The spatial autocorrelation analysis method can reflect the spatial distribution of ESV itself and the characteristics of landscape patterns in Zhongmu County. Spatial autocorrelation analysis measures whether the distribution of elements is spatially related to neighboring elements [35], including global spatial autocorrelation and local spatial autocorrelation. Global autocorrelation is used to describe the degree and significance of the spatial correlation between variables in a region. The local spatial autocorrelation is usually analyzed by Moran's I scatter plot and local indicators of spatial association (LISA). Moran's I scatter diagram can reflect the correlation between spatial unit elements and neighboring elements. At the same time, LISA decomposes global Moran's I into each spatial unit and forms a LISA clustering map, which can show the specific location of the spatial aggregation or differentiation of ESV and landscape pattern [36,37], the correlation mode of Moran's I scatter plot and LISA cluster plot is shown in Table 1. The following are the calculation formulas for the global Moran's I index and the local Moran's I index: where I is global Moran's I; n is the total number of spatial units; x is the element value; W ij is the spatial weight, and S 2 is the variance. This study builds W ij based on the queen adjacency criterion commonly used by the GeoDa software. When the spatial units are adjacent, the value of W ij is 1, otherwise, it is 0. The value of I is (−1, 1). The closer I is to 1, the stronger the positive spatial correlation of the elements. The more concentrated the distribution, the closer to −1, the stronger the negative spatial correlation and the greater the distribution difference; equal to 0 it shows no spatial autocorrelation, and the distribution is random [37,38].

Scatter Image Limit LISA Gathering Spatial Correlation Description
First quadrant High-High Cluster Positive spatial correlation The element values of the spatial unit and its forest are high, and the spatial difference is small.

Second quadrant Low-High Outlier Spatial negative correlation
The element value of the spatial unit is low, the element value of its neighborhood is high, and the spatial difference is large.
Third quadrant Low-Low Cluster Positive spatial correlation The feature value of the spatial unit and its neighborhood is low, and the spatial difference is small.

Fourth quadrant
High-Low Outlier Spatial negative correlation The higher the element value of a spatial unit, the lower the element value of its neighbors, and the greater the spatial difference.

Dynamic Analysis of the Overall Landscape
The landscape composition of Zhongmu County was dominated by dry land. Large paddy fields were distributed from east to west in the north, urban land was mainly distributed in the central and western regions, and rural residential areas were widely scattered in the central and southern regions (  (Table 2), and the values of CONTAG and CO-HESION showed a downward trend. It showed that the density of landscape edges in Zhongmu County increased during the study period. With the increase in patches per unit area, patch shapes became more and more complex, and landscape diversity increased. However, the degree of fragmentation of the landscape was also growing year by year, and the overall landscape connectivity was decreasing and the degree of patch dispersion was increasing year by year.  Table 3 showed that dry land was the main source of increase in other landscape types. The reduction in dry land from 2005 to 2018 was 206.80 km 2 , which was the largest reduction among all landscape types, and the landscape type with the largest transfer area of dry land was rural residential land and the deduction was 73.80 km 2 . In addition, the size of urban land, rural residential land and other construction land transferred to woodland, gardens, high-coverage grassland, canals and reservoir pond did not exceed 1.15 km 2 , indicated that less construction land was transformed into ecological landscapes. The input source of the ecological landscape was relatively single, and the input was much lower than the output.  (Table 2), and the values of CONTAG and COHE-SION showed a downward trend. It showed that the density of landscape edges in Zhongmu County increased during the study period. With the increase in patches per unit area, patch shapes became more and more complex, and landscape diversity increased. However, the degree of fragmentation of the landscape was also growing year by year, and the overall landscape connectivity was decreasing and the degree of patch dispersion was increasing year by year.  Table 3 showed that dry land was the main source of increase in other landscape types. The reduction in dry land from 2005 to 2018 was 206.80 km 2 , which was the largest reduction among all landscape types, and the landscape type with the largest transfer area of dry land was rural residential land and the deduction was 73.80 km 2 . In addition, the size of urban land, rural residential land and other construction land transferred to woodland, gardens, high-coverage grassland, canals and reservoir pond did not exceed 1.15 km 2 , indicated that less construction land was transformed into ecological landscapes. The input source of the ecological landscape was relatively single, and the input was much lower than the output.  Table 4). The increase rate was 3.08%. Among all types of landscape ESV, the ESV of dry land always remained the highest from 2005 to 2018, and the average contribution rate of dry land was as high as 37.46%; while the average contribution rate of high-coverage grassland ESV was the smallest, 0.76%; the ESV of dry land and paddy fields, and their contribution rate showed a declining trend; while the ESV contribution rate of reservoir ponds and canals was constantly increasing.
Among the individual service values in Zhongmu County, as shown in Table 5, the average contribution rate of adjustment services was the highest, accounting for 58.61% of the overall ESV. Among the 11 secondary ecosystem service functions, the ESV of hydrological regulation has always remained the largest, with an average contribution rate of 33.64% during the study period. The ESV of water resources supply always remained the smallest during the study period, with an average contribution rate of only 0.58%. From the perspective of the rate of change, only the rate of change in hydrological regulation ESV from 2005 to 2018 was positive and showed a trend of first decreasing and then increasing; while the rate of change in the other four types of service functions were all negative, including soil conservation; the rate of reduction in gas regulation and food production was increasing.   Figure 3 and Table 6 showed that the spatial distribution of total ESV in Zhongmu County from 2005 to 2018 was dominated by low-value areas. The average ESV of 80.01% of the spatial units was less than USD 0.24 million. ESV high-value areas were distributed in the Yellow River in the north, and the average proportion of spatial units was only 0.53%; ESV higher-value areas were mainly distributed in the Yellow River, Yanming Lake, northwest paddy fields and water areas in the north, the average ratio was only 2.36%; there were mainly ESV low-value areas in the middle and southern part of Zhongmu County.     From the perspective of the spatial distribution of ESV changes from 2005 to 2018, the ESV changes were generally between USD−0.22 and 0.21 million, and the areas with an increase of USD 1.00 to 2.91 million only existed in local areas of the Yellow River. In addition, the ESV variation in the Jialu River Basin in the central part was significantly higher than that of the surrounding areas, forming a northwest-southeast belt-shaped area in the entire Zhongmu County. At the same time, the ESV variation in the belt-shaped area had obvious spatial differentiation characteristics: the ESV in the upper reaches of the northwestern is generally in a state of reduction, and the reductions are all more than USD 0.22 million; the opposite is true in the middle and lower reaches of the central and southeastern regions, the growth of ESV both exceeded USD 0.21 million, mainly due to the influence of the "Zhengzhou-Kaifeng integration" policy.

Landscape Transfer of ESV Changes in Cold/Hot Spots
From 2005 to 2018, the distribution of ESV hot spots and their landscape shifts were as follows: (1) The northern Yanming Lake and Yellow River areas mainly shifted from cultivated land to waters. (2) There were a small number of scattered hot spots in the central and southern regions, mainly shifting cultivated land to waters and forests, respectively. The distribution of ESV cold spots and their landscape shifts were as follows: (1) Northeast of Zhongmu County, mainly due to the conversion of waters to cultivated land. (2) The Baisha cluster area in the west was mainly converted from waters to cultivated land and cultivated land to construction land ( Figure 4).
Synthesizing the landscape shifts in the cold and hot spots in the three periods, it was found that the change in cultivated land to waters and forests was the main area in the hot spots of ESV change. Waters and forests contributed to hydrological regulation, water supply, climate regulation, and other service functions. Therefore, ESV in areas where cultivated land is transferred to waters and forests increased significantly. In the cold spot area where ESV changes, waters were mainly assigned to construction land and cultivated land. This was primarily due to the weak ESV of construction land and food production and consumption, hydrological regulation, soil conservation, biodiversity, and other ecological service functions. Furthermore, the ecological contribution per unit area of cultivated land was much smaller than that of water, so ESV decreased significantly in areas where water was transferred to construction land and cultivated land. Table 7 showed that the global Moran's I values of PD, ED, LSI, SHAPE_MN, FRAC_MN, and ESV were all greater than 0, showing a positive spatial correlation, indicating that the higher the degree of fragmentation of the landscape, the higher the richness of the patch, and the more complex the shape of the patch, the higher the ESV of its neighborhood. The connectivity index, CONTAG, COHESION, and ESV global Moran's I values were all negative and relatively close to 0, indicating a weak spatial negative correlation with ESV. When the connectivity of the spatial unit is high, the ESV of its neighborhood is low. Land 2021, 10, x FOR PEER REVIEW 12 of 22  Table 7 showed that the global Moran's I values of PD, ED, LSI, SHAPE_MN, FRAC_MN, and ESV were all greater than 0, showing a positive spatial correlation, indicating that the higher the degree of fragmentation of the landscape, the higher the richness of the patch, and the more complex the shape of the patch, the higher the ESV of its neighborhood. The connectivity index, CONTAG, COHESION, and ESV global Moran's I values were all negative and relatively close to 0, indicating a weak spatial negative correlation with ESV. When the connectivity of the spatial unit is high, the ESV of its neighborhood is low.   The global Moran's I of ESV and landscape level pattern index in 2005, 2010, and 2018 were obtained through the GeoDa software. After 999 random permutations, all p-values were 0.001, which shows that the spatial autocorrelation between ESV and landscape index is significant under 99.9% confidence.

Spatial Correlation between ESV and Landscape-Level Index
The LISA clustering map of PD, ED, LSI, SHAPE_MN, FRAC_MN, and ESV was mainly high-high clusters and low-low clusters. Figure 5 showed that the high-high clusters of the above index and ESV were primarily located in the Yanming Lake, the Yellow River, and the paddy fields and reservoir pond areas in the northwest corner of Zhongmu County. It is shown that the higher the abundance of these three landscape patches, the higher the ESV of the neighborhood. Therefore, the number and shape complexity of landscape patches in Yanming Lake, the Yellow River, and the northwest corner of Zhongmu County can be increased to improve the ESV of Zhongmu County effectively.

Landscape Elements and ESV
From 2005 to 2018, as the areas of canal, reservoir pond, and woodland increased by 5.36 km 2 , 16.40 km 2 , and 8.03 km 2 , respectively, and the areas of dry land, paddy fields and gardens decreased by 84.24 km 2 , 43.98 km 2 , and 20.37 km 2 , respectively. The ESV of The LISA clustering map of CONTAG, COHESION, and ESV mainly focused on the spatial aggregation modes of low-high outliers and high-low outliers. The low-high outliers were mainly distributed in the north, and the high-low outliers were mainly distributed in the south. The distributions of the high-high clusters and low-low clusters of PD, ED, LSI, SHAPE_MN, FRAC_MN, and ESV were consistent, which means that the higher the connectivity of landscape patches in the Yellow River, Yanming Lake and the northwest region, the lower the overall ESV, which further proved that attention should be focused on increasing the richness and complexity of the three regional landscape patches to achieve an efficient increase in ESV.

Landscape Elements and ESV
From 2005 to 2018, as the areas of canal, reservoir pond, and woodland increased by 5.36 km 2 , 16.40 km 2 , and 8.03 km 2 , respectively, and the areas of dry land, paddy fields and gardens decreased by 84.24 km 2 , 43.98 km 2 , and 20.37 km 2 , respectively. The ESV of hydrological regulation was compared with soil conservation and gas regulation functions, the only performance was positive growth. At the same time, through the superposition analysis of the ESV variation cold, hot spots, and landscape transfer matrix maps, it is found that the ESV in the area where the cultivated land was transferred to the water area and the forest was significantly increased. The ESV, where the water was transferred to the construction land and the cultivated land, was reduced substantially. The above conclusions fully reflected that the changes in elements in the landscape pattern of Zhongmu County had a significant impact on ESV. This was more consistent with previous research conclusions: For example, Hou et al. [39] found that the conversion of forest land to construction land in Xi'an had the greatest negative impact on ESV, followed by the conversion of grassland to construction land; Assefa et al. [40] found that the area of Bahir Dar City and its surrounding wetlands and water bodies in Ethiopia decreased by 1618 hectares, resulting in a reduction in the total ESV of USD 8.9 million, and the reduction in ESV had a negative impact on the local climate, waste management and the livelihoods of poor communities.

Landscape Index and ESV
In the Yanming Lake, Yellow River in the north, the paddy fields, and reservoir pond areas in the northwestern corner in Zhongmu County, ESV was positively correlated with PD, ED, SHAPE_MN, FRAC_MN, and LSI, and negatively correlated with CONTAG and COHESION, indicating that the more landscape patches in the Yanming Lake, Yellow River and northwest corner, the higher the complexity and the higher the ESV. Lin et al. [41] discussed in the study of Guangdong-Hongkong-Macao Greater Bay Area that abundant patch types in the landscape were conducive to improving soil conservation ability, so ESV also increased correspondingly, while enhanced connectivity was not conducive to maintaining soil conservation ability, so ESV was reduced. In addition, other scholars also had the opposite conclusion. For example, through correlation analysis, Gu et al. [42] found that ESV in Lancang County was negatively correlated with PD and LSI, but positively correlated with CONTAG; Qiu et al. [43] found in Fuyang City that the total ESV negatively correlated with ED, PD, LSI, and SHDI, suggesting that landscape fragmentation had an adverse impact on the overall ESV in Fuyang City. Based on the above research conclusions, it can be concluded that the impact of landscape pattern index on ESV is complicated. Cen et al. [19] also proved a non-linear and complex dynamic relationship between ESV and landscape pattern index in the study of the south bank area of Hangzhou Bay. Therefore, the correlation and interaction degree of ESV and the landscape pattern index are different in different research areas or different spaces in the same research area. It is necessary to analyze other regions to propose targeted landscape pattern control strategies, which effectively improve the level of regional ES.

Driving Force Analysis
Most studies have shown that policy factors are the main driving forces for landscape patterns and ESV changes: Raviv et al. [44] found that in the Serekamer Mountain Biosphere Reserve, although wildfires and land cover changes affect the ecosystem, good land restoration programs increase the value of ecosystem services by 30-50%; Mendoza-González et al. [45] proved that the development of tourism policies in the cen-tral Gulf of Mexico led to the loss of the value of ecological services such as coastal protection. In this study, the implementation of urban planning policies and ecological policies at different levels during 2005-2018 was also the main driving force for the changes in landscape patterns and ES in Zhongmu County.
In terms of urban planning policies, the Zhengzhou-Kaifeng integration played an important role in promoting the development of Zhongmu County. In 2005, when the concept of Zhengzhou-Kaifeng integration was formed, Zhongmu County's arable land, mainly dry land and paddy fields accounted for 36.70%, and construction land accounted for only 4.54%. From 2005 to 2018, the Zhengzhou-Kaifeng integration was promoted by the state [25]. In 2016 the National Development and Reform Commission officially approved Zhengzhou as a national central city, promoting Zhengzhou to lead the construction of a metropolitan area, which further promoted the urbanization development process of Zhongmu County in the Zhengzhou-Kaifeng integration [46]. During this period, Zhongmu County initially established a comprehensive transportation system that coordinated with the Zhengzhou-Kaifeng Integration strategy. At the same time, construction land was continuously increasing, such as industrial parks and residential communities. The urban development pattern of "three axes, three districts, and multiple clusters" was initially formed [26]. This further promoted the changes in the landscape pattern in Zhongmu County in the Zhengzhou-Kaifeng integration. In 2018, the proportion of cultivated land decreased to 32.47%. The ratio of construction land rapidly increased to 8.64%, which led to a negative growth in soil conservation and gas regulation functions that have contributed much to the county's ecological value.
In terms of environmental policies, the concepts emphasized by the central government, such as "ecological civilization", "green development", "ecological city", and "ecological red line", provide national-level policy guarantees for the sustainable development of the ecological environment in Zhongmu County. From 2005 to 2018, the Outline of the Central Plains Urban Agglomeration Overall Development Plan proposed that the Zhengzhou-Kaifeng integration area build green corridors on both sides of the main arterial roads, construct forests and greenbelt ecological adjustment zones in feasible areas, and develop urban agriculture and sightseeing agriculture to realize the ecological docking between Zhengzhou and Kaifeng [47,48]. The Zhengzhou Municipal Government successively implemented the Zhengzhou City Master Plan (2010-2020), Zhengzhou City Green Space System Planning (2013-2030), Zhengzhou National Central City Ecological Construction Plan (2016-2025), Zhengzhou Forestry Ecological Construction Implementation Plan in 2017, and the Zhengzhou City Ecological Protection and Construction Planning (2017-2035) which put forward construction along the Yellow River ecological cultural tourism industry area, the Yellow River ecological zone in the north, the Yellow River ecological wetland landscape zone, the ecological barrier along the Yellow River zone in the north of Zhengzhou, and the Yellow River ecological conservation and ecological culture development belt, respectively, by returning farmland to forest, widening the drainage of the levee protection, ecological hydrating and other measures. This greatly improved the ESV in the Yellow River Basin and Yanming Lake region in the north, which is also the reason why the area with the increase in ecological value as high as USD 1.00 million to USD 2.91 million only existed in the Yellow River region, and the high-value accumulation area of the total ESV was mainly distributed in the north and the range was increasing.

Optimization Measures and Outlook
This study analyzed the spatial and temporal changes in landscape patterns and ecosystem service values in the study area in the context of integrated urban development. Based on the results of the study, the following optimization measures are proposed for the future, to provide recommendations for the future by means of planning policies and strategies.
(1) Enhance the development the quality of green infrastructure. Improve the development level of green infrastructures such as rivers, wetlands, and woodlands to enhance the function of ecosystem services such as hydrology, climate, and soil in the region and effectively increase the total amount of ESV; (2) Control the unlimited expansion of construction land and strengthen the core ecological infrastructure. We should strictly control the unlimited expansion of construction land, optimize the layout of construction land, strengthen the intensive use, and avoid the massive erosion of arable land, water and other high-quality landscapes. At the same time, provide full play to the core advantages of the central city and strengthen the construction of "ecological green core: Zhengzhou-Kaifeng Central Park" and "ecological green core: comprehensive parks, special parks, community parks, pocket parks, etc." Form a comprehensive park system to ensure the high quality of life for the high-density population in the central region and the stable growth of ESV; (3) Promote high-quality development of watershed ecological reserves. The Yellow River basin is the most important ecological source [49]. Further, promote the highquality development of the ecological reserve along the Yellow River. Gradually improve the landscape pattern of the Yellow River basin with "wetland-protection forest-agricultural land", and restore the value of the Yellow River ecosystem services.

Limitation
Several limitations in this paper should be acknowledged. Firstly, the resolution of the acquired remote sensing interpretation data were limited, and unnecessary errors may occur due to the vacancy of the data itself. Secondly, the consideration of social and human factors was relatively inadequate. In the future, the accuracy of remote sensing interpretation data and social and human factors should be considered as comprehensively as possible.

Conclusions
Based on the remote sensing interpretation data and socioeconomic data of Zhongmu County in 2005, 2010 and 2018, this study analyzed the spatial pattern dynamics of the overall and various types of landscapes in combination with the landscape index, and analyzed the dynamic process of ecosystem service values in Zhongmu County based on the equivalent factor method. Finally, the spatial relationship between landscape patterns and ecosystem services was analyzed comprehensively, and the influence of policy drivers was systematically sorted. Ecological problems and optimization measures in Zhongmu County were proposed in response to the research results. The conclusions of this paper are as follows: (1) Changes in the total value of ecosystem services. The total value of ecosystem services in Zhongmu County in 2005, 2010 and 2018 showed a fluctuating trend of first increasing and then decreasing, with the total value increasing by USD 10.05 million, and the ecological environment construction achieved certain results. Spatially, the low-value area was the main area (80.01%), the high-value area was only in the local area of the Yellow River (0.53%), and the higher-value area was mainly in the Yellow River, Yanming Lake, and the paddy fields and reservoir ponds (2.36%) in the northwest corner. The spatial aggregation pattern of ecological values has a border effect: the high-high clusters area was mainly at the northern Yellow River border, while the low-low clusters area was at the east-west border, and there is an obvious trend of expansion to the east; (2) Changes in the ecosystem service values of various landscapes. The ecosystem service value of dry lands always remains the highest, with an average contribution of 37.46% to the overall value, and was the most dominant landscape in Zhongmu County in maintaining ecosystem service value, followed by reservoir ponds (14.65%) and canals (10.51%), with changes mainly related to the increase or decrease in the corresponding landscape area and value coefficients; County were hydrological regulation (33.64%), soil conservation (12.71%) and climate regulation (10.16%), among which only the ecological value of hydrological regulation had positive growth; (4) Through the overlay analysis of the cold, hot spot and landscape transfer matrix mapping of the amount of ecosystem service change from 2005 to 2010, 2010 to 2018 and 2005 to 2018, it was found that the value of ecosystem services in the areas where cropland was transferred to waters and forests increased significantly, and the value of ecosystem services in the areas where waters were transferred to construction land and cropland decreased significantly. Therefore, we can focus on increasing watershed and forest in the landscape regulation to improve the ecosystem service value; (5) Through bivariate spatial autocorrelation analysis, it is found that the ecosystem service value can be effectively increased by enriching the landscape types and increasing the number and complexity of landscape patches in the region.
The study had important implications for the planning and management policies for Zhongmu County and other urban developing regions: (1) urban integration areas should strictly limit the expansion of construction land, optimize the layout of construction land to enhance its intensive use, and increase the plot ratio and use efficiency on the limited construction land; (2) improve the development level of green infrastructures such as rivers, wetlands and woodlands, so as to enhance the function of ecosystem services such as hydrology, climate and soil in the region and effectively increase the total amount of ESV.

Data Availability Statement:
The data presented in this study are available on request from the first author.

Acknowledgments:
The authors would like to thank the support of the International Joint Laboratory of Landscape Architecture, Henan Agricultural University, for their infinite help.

Conflicts of Interest:
The authors declare no conflict of interest. Table A1. Description of landscape classification in Zhongmu County.

Cultivated land
Paddy field Irrigated arable land for growing rice, lotus root, and other aquatic crops, including arable land where rice and dry land crops are rotated.
Dry land Cultivated land that relies on natural precipitation to grow crops and cultivated land that mainly grows vegetables.

First Classification Secondary Classification Description
Forest Woodland Natural forests and artificial forests (including timber forests, economic forests, shelterbelts, etc.) with a canopy closure greater than 30%, low forest land and shrubland with a canopy closure greater than 40%, and a height of fewer than 2 m.

Garden plot
Unformed forest afforestation land, site, nursery and orchard, mulberry garden, tea garden, hot cropping forest garden, etc.

Grassland
High-coverage grassland Natural grassland, improved grassland and mowing grass with a coverage >50%. Generally, the water condition is good and the growth is dense.

Water area
Canal Naturally formed or artificially excavated rivers and artificial canals.

Reservoir pond
The land below the perennial water level in the artificially constructed water storage area.

Bottomland
The land between the water level of rivers and lakes during the normal period and the water level of the flood period.

Construction land
Urban land Large, medium and small cities and built-up areas above counties and towns. Rural residential land Independent of rural settlements outside of cities and towns.

Other construction lands
Lands for factories and mines, large industrial areas, oil fields, saltworks, quarries, etc., as well as traffic roads, airports, and special land. Table A2. Landscape pattern index and its ecological significance.

Landscape Pattern Index Ecological Significance
Percentage of Landscape (PLAND) Reflects the relative proportion of a certain landscape type to the entire landscape area [50].
Largest Patch Index (LPI) It reflects the proportion of the largest patch in the landscape area and determines the dominant type in the landscape [51].
Edge Density (ED) Reflect the degree of differentiation or fragmentation of the overall landscape patch [52].

Mean Patch Fractal Dimension (FRAC_MN)
Reflect the complexity of the patch shape of the landscape type. The larger the value, the greater the space occupied by the landscape type, and the more complex its shape [53].

Mean Patch Size (SHAPE_MN)
Reflects the degree of deviation of the plaque shape from the standard shape (round or square). When all patches in the landscape are square, the value is 1; when the shape of the patches deviates from the square, the value increases [54].

Patch Density (PD)
Reflecting the number of patches per unit area and characterizing the complexity of landscape spatial structure, is an important indicator for evaluating landscape fragmentation [55].

Patch Cohesion Index (COHESION)
Reflecting the physical connectivity of patches, the smaller the value, the more scattered the landscape patches [56,57].
Landscape Shape Index (LSI) Reflecting the overall geometric complexity of the landscape, the larger the value, the longer and irregular the boundary of the patch, the more discrete the patch, that is, the higher the complexity or fragmentation of the landscape [58].

Contagion Index (CONTAG)
Reflects the degree of agglomeration or spreading trend of different plaque types.
To measure to what extent landscapes are aggregated or clumped as a percentage of the maximum possible [59].

Shannon's Diversity (SHDI)
A measure of patch diversity in a landscape, which is determined by both the number of different patch types and the proportional distribution of area among patch types [60].