Abstract

The future of ecology lies on the promotion of eco-economy, a resource saving and environment-friendly society, and the research of ecosystem services. This article applied weighted Voronoi diagram to simulate the ecosystem service capacity in Zhongyuan Urban Agglomeration, China, using remote sensing and geographic information systems. The results show that (1) compared with the traditional Voronoi diagram, the weighted Voronoi diagram expresses more accurate on the ecosystem service range, (2) the total value of the ecosystem service in the study area was 575.7 billion yuan in 2015.The maximum value was about 253.8 billion yuan greater than the minimum value. The unit forest ecosystem has been improved after the coefficient was modified, and the value of ecosystem services reached 32,638.87 yuan/hm2, (3) the service value of construction ecosystem has a great effect on the total value of ecosystem services; the overall pulldown effect of Kaifeng City is most considerable, and (4) the region with strong ecosystem service capacity has strong effect in annexing the nearby areas with weak ecosystem service capacity. And the distance of the annexing is inversely proportional. The closer the distance is, the stronger the annex will be.

1. Introduction

The survival of human societies is determined by the carrying capacity of their ecosystems [1, 2]. The idea that natural resources and environments are free causes ecological resource development to exceed ecosystems’ basic carrying capacities. This, in turn, leads to the reduction or complete loss of ecosystem service functions [3, 4]. Now, more than two-thirds of services provided by ecosystems worldwide have already begun to decrease [5]; this trend may not be effectively reversed in the next 50 years [3, 4, 6]. Resource, ecological, and environmental crises have already seriously threatened China’s development; for China to develop in a sustainable manner, scientific studies must establish macroscopic ecological security patterns and define key ecological spaces for prior protection. Therefore, exploring a scientific and effective method for delineating the range of ecosystem services is of positive significance in promoting the balanced development of ecosystem services and establishing and improving the ecosystem protection policy.

Ecosystem service functions are defined as “the natural benefits which human survival relies on and which are created and maintained by ecosystems and ecological processes” [7]. To describe ecosystem service functions more intuitively, Costanza et al. [8] conducted an assessment of ecosystem service values, which has become a key ecological research topic. Estoque and Murayama [9] found that the area of urban in the city of Baguio, the summer capital of the Philippines, expanded continuously during the study period, while the ecosystem service value dropped to 40% of that at the beginning of the study. Gómez-Baggethun and Barton [10] discovered that traditional monotonous urban planning is one-sided, and that a single socioeconomic indicator cannot assure the balanced ecological and economic development of a society. Moreover, introducing ecosystem service values has a positive impact on the development of human societies. Su et al. [11] primarily studied the Hangjiahu Plain agroecological zone, the Tianmu Mountain forest ecological zone, the Hangzhou urban ecological zone, and the Qiandao Lake Basin forest and wetland ecological zone and analyzed the effects of urbanization on regional ecosystems. Using research from around the world [1217], Costanza et al. [6] summarized ecosystem service value assessments and further improved their evaluation. However, key research on ecosystem service values remains: for example, how should coverage of ecosystem services be delimited? How should a method of ecosystem service coverage quantification be developed?

Voronoi diagrams are partitioning that have been shown to be highly similar to some natural structures; as such, they have been studied extensively in fields including meteorology, geology, surveying and mapping, archaeology, molecular chemistry, ecology, and computer science [18, 19]. Voronoi diagrams are mainly applied to site selection and the determination of the scope of influence of different objects. Qin and Shi [20] utilized Voronoi diagrams to examine the spatial distribution of large-scale residential areas, ultimately making suggestions for planning in the region. Fan et al. [21] integrated Voronoi diagrams with spatial interaction models to investigate the Horqin Left Middle Banner, in Inner Mongolia, and proposed a method for improving this model by substituting town centrality for population size and integrating the model with weighted Voronoi diagrams. Traditional Voronoi diagrams are often used to describe the evenness of ground objects and are based on the principle of even distribution; however, they can better solve practical problems after some reasonable improvement. For example, introducing weights can compensate for the defects of traditional uniform Voronoi diagrams [2225]. Li et al. [26] and Miao et al. [27] used Voronoi diagrams to examine parks in Shenyang and the wetland functional zone in Beijing, respectively, confirming the possibility of employing this method in ecosystem research. Nevertheless, there is still room for development in studying ecosystem services and enhancing its assessment accuracy.

To explore a suitable method of studying ecosystem service coverage areas, we introduce Voronoi diagrams and weighted Voronoi diagrams into ecosystem service value research and analyze ecosystem services coverage based on actual conditions. This approach helps quantify ecosystem service research not only at fixed points but also across regions. The Zhongyuan Urban Agglomeration is one of seven major urban agglomerations in China; it is the key development area of the Zhongyuan Economic Zone. Recent rapid social and economic development in this area has changed the types of ecosystems in the region, which in turn has led to changes in ecosystem service values in the coverage areas. In order to explore a reasonable method for estimating ecosystem service coverage areas, we will value the ecosystem services, that is, use economic methods to describe the strength of ecosystem services. Furtherly, we introduce the Voronoi diagram and the weighted Voronoi diagram to simulate and analyze the ecosystem service areas. The Voronoi diagram, the Euclidean distance map, and the rank assignment map are compared at various levels to explore the practicality of the weighted Voronoi diagram, thereby improving the spatially estimation practices of ecosystem service and providing the scope for future ecosystem service research.

2. Study Area and Data Resources

2.1. Study Area

The Zhongyuan Urban Agglomeration (shown in Figure 1) is located in central Henan Province, China; it has a total area of 58,700 square kilometers and includes nine cities: Zhengzhou, Luoyang, Kaifeng, Xinxiang, Jiaozuo, Xuchang, Luohe, Pingdingshan, and Jiyuan. This area is an important growth pole, not only of Henan Province but also of the Zhongyuan Economic Zone. In 2015, the gross national product of the agglomeration was 2.19 × 109 yuan and its population was 4.15 × 108. It has a nice natural environment with a wide distribution of varied land resources, a temperate monsoon climate, high topography in the west and north, and low topography in the east and south. The eastern part of the Zhongyuan Urban Agglomeration is located in the North China Plain, which is an essential grain production base and consists mainly of farmland. The Taihang Mountains in the west and north are important production sites of forest resources. Given the premise that ecological and economic development must be constantly promoted, it is necessary to protect regional ecologies to maintain balance in regional ecosystem services. Thus, we investigate the ecosystem service coverage of different ecosystems in the region and make suggestions for regional planning.

2.2. Data Sources

We aim to verify the applicability of weighted Voronoi diagrams to ecosystem service value research and place little emphasis on spatial and temporal variations in ecosystem service values. Weighted Voronoi diagrams of the ecosystem service value of the Zhongyuan Urban Agglomeration in 2015 are created to examine their meaning in ecosystem service research. Land use types significantly control the functions of regional ecosystem services and are commonly used as indictors of the values of regional ecosystem services [28]. We acquired land use map of the study area by interpreting Landsat OLI images in April 2015, which is in the growing seasons and the images with low cloud cover percentages were selected. We used a manual interpretation method to derive land use maps from the remote sensing images. The land use is classified into four main types: farmland, forest, water, and construction area. To evaluate the accuracy of the land use map, we selected 110 random points from the map. The Land Use Atlas in Henan Province (1 : 1 : 100000) and high-resolution satellite images from Google Earth were referred to validate the map. As a result, the overall accuracy of the land use map is 91.3%, and the Kappa coefficient is 0.88.

3. Research Methods

Based on Costanza’s and Daly’s [13] measurement and calculation of global ecosystem service values and the “Ecological Service Value Equivalence per Ecosystem Unit Area in China” table established by Xie at al. [29] based on actual conditions in China, this paper considers the actual condition of the Zhongyuan Urban Agglomeration to propose a novel research method, which is introduced into the correction coefficient of the study area’s ecosystem service value. ArcGIS 10.1 and MATLAB are utilized to describe the spatial distribution of ecosystem service values in the Zhongyuan Urban Agglomeration. The results are compared in color gradings, traditional Voronoi diagrams, and weighted Voronoi diagrams to analyze these methods’ advantages and drawbacks in ecosystem service value assessment (Figure 2) and provide a new perspective on ecosystem service value research.

3.1. The Correction Coefficient Model of Ecosystem Service Values

Based on previous research [13, 29, 30] and real-world regional conditions, the assessment coefficient of ecosystem service values (as shown in Table 1) is improved. We propose an improved model based on grain yields (1) and a correction coefficient model of ecosystem service values (2). where Ea is the value of food supply services per unit agroecosystem per unit ecosystem (yuan), is the crop area (hm2), is the crop yield per unit area (kg/hm2), is the notional average price of crops (yuan/kg), is the total sown area (hm2), denotes the equivalent value of ecosystem for function , and is the function of ecosystem service in the th ecosystem type relative to the food supply services provided by the agroecosystem. The of the study area was calculated to be 1493.77 yuan/hm2. The correction coefficient of ecosystem service value is listed in Table 1.

3.2. The Ecosystem Service Value (ESV) Model

The ecosystem service value (ESV) model is effective for estimating ecosystem service values based on different ecological types [30] and is expressed as illustrated in (3) and (4). where is the ecosystem service value (yuan), is the area of the ith ecosystem type, and denotes the equivalent coefficient of function for ecological land .

3.3. Voronoi Diagrams

Voronoi diagrams (abbreviated as V diagrams, as shown in Figure 3) were proposed by Ukrainian mathematician Voronoi in 1908 and were applied to high-dimensional spaces. By assuming that for a set of discrete points on Euclidean plane , where any two arbitrary points are noncoincident (i.e., ) and any four arbitrary points are nonconcyclic, is calculated as where is the Euclidean distance between points and . Thus, the distance from each point that belongs to is closer to than any other points of discrete point set .

3.4. Weighted Voronoi Diagrams

Although Voronoi diagrams can depict the basic characteristics of some natural environments, they are unable to present the condition of some relatively complex ground objects. Based on traditional Voronoi diagrams, weighted Voronoi diagrams (as shown in Figure 4) introduce weight indices and convert traditional Euclidean distances into weighted Euclidean distances [31]. Given that , nonnegative real number weights are introduced; hence, the weighted distance between points and is

Then, should be revised to

When , weighted Voronoi diagrams become traditional Voronoi diagrams.

4. Results and Analysis

4.1. Ecosystem Service Value Assessment

The Zhongyuan Economic Zone has lately been an important economic development region in China; more importantly, it is also the key region for stimulating the development of the central region of China and promoting balance in China’s economic development. Because the Zhongyuan Urban Agglomeration is the key region of the Zhongyuan Economic Zone, its overall economic development acts as a model for the development of other regions. Therefore, proper ecosystem planning in the agglomeration and balanced ecological and economic development of the region enhances its competitiveness and encourages the development of other regions.

Before ecosystem service value assessment was completed, Costanza and Daly [13] and Xie et al. [29] performed coefficient correction of ecosystem service values according to real-world conditions in the study area. The results are listed in Table 2.

The coefficients of ecosystem service values in Table 2 and the ecosystem of different land use areas in the Zhongyuan Urban Agglomeration after satellite image interpretation were integrated to assess the ecosystem service value of each land use of the region. The results are shown in Figure 5.

The total ecosystem service value of the Zhongyuan Urban Agglomeration in 2015 was 575.7 billion yuan. The ecosystem service value of Luoyang was the highest, accounting for 43.93% of the regional total; in contrast, that of Luohe was only 1.42% of the regional total. Furthermore, the forest ecosystems of Luoyang contributed the most to the total ecosystem service value of the city, accounting for 89.09% of the total. The total ecosystem service value of Kaifeng was 28.259 billion yuan, of which the value of farm ecosystem was 48.216 billion yuan, whereas the loss in value caused by construction ecosystem expansion was up to 28.037 billion yuan. The ecosystem types combined with the administrative area could lead to the differences of ecosystem service values. Specifically, higher proportion of forest and water ecosystems in a region would increase the region’s ecosystem service values and furtherly have positive influences upon the neighbor ecosystems.

4.2. Simulation of Ecosystem Service Coverage

Similar to social services, ecological services have coverage directly related to their intensity. Past ecosystem service studies could not determine the areas of influence of multiple service functions. For example, the ecosystem service function of climate regulation by forest ecosystems is affected by several factors and it is not possible to accurately assess specific service coverage. Hence, selecting an appropriate angle for measuring and calculating ecosystem service coverage can compensate for shortcomings in previous research.

Ecosystem service values quantify all kinds of ecosystem service functions and present them intuitively. The appearance of this variable allows ecosystem service coverage to be simulated and analyzed using an economic approach. The resulting Voronoi diagrams are shown in Figure 6 in a map showing a uniform distribution; however, there are considerable differences in ecosystem service values between cities. A uniform distribution cannot correctly depict the actual coverage of ecosystem services. Therefore, ecosystem service values of different regions are presented as weighted Voronoi diagrams, which partly compensate for defects in the expression of service coverage by traditional Voronoi diagrams.

The ecosystem service values of all cities in the Zhongyuan Urban Agglomeration were used as weight indices to improve the original Voronoi diagrams, as illustrated in Figure 7. These results, together with the data in Figure 5, demonstrate that the ecosystem service value of Luoyang accounts for a significant portion of the total ecosystem service value of the Zhongyuan Urban Agglomeration. Its influence is considerably greater than those of other cities. Integration of the different service boundaries of Figures 5 and 7 shows that, taking Luoyang as an example, the ecosystem service value coverage of the city expands significantly, which is highly related to the larger ecosystem service value of the city. For regions with lower service values, such as Jiyuan and Luohe, except for certain preserved coverage, areas are impacted by regions with greater ecosystem service values. That is, if the ecosystem service value is greater, the service coverage is more extensive. To an extent, these areas affect and compensate for areas with smaller ecosystem service values.

To highlight the interactions and relationships between ecosystem services of different regions, we classified ecosystem service coverage using the Jenks natural breaks classification method into five types: high-, medium to high-, medium-, medium to low-, and low-impact areas. The traditional Voronoi diagrams are compared with the weighted ones in Figure 8.

As shown in Figure 8, regions may have various effects on each other, although the center of each region is mainly influenced by the region itself. In addition, regions with greater ecosystem service capabilities compensate for regions with lower capabilities, which conceal the originally insignificant ecosystem service capabilities of regions with lower service capacities. The ecosystem service interactions engulfing different regions can be monitored and controlled artificially. Rational distribution of ecosystems over the study area will considerably enhance a region’s ecosystem service capability and thus its ecosystem service quality. If the stepwise progressive decreases of different kinds of ecosystem service capabilities can be successfully determined, it would be possible to more accurately evaluate the influences and coverage of ecosystem services.

From the total ecosystem service value in Figure 8, the type of ecosystem that dominates the total ecosystem services can be seen and various types of ecosystem service coverage can be studied. Thus, the uneven distribution of land use types between regions can be obtained. We illustrate the ecosystem service capabilities and coverage of farm ecosystem, forest ecosystems, water ecosystems, and construction ecosystems of regions in the Zhongyuan Urban Agglomeration in Figures 912.

Excluding Jiyuan, Luohe, and Jiaozuo, the farmland ecosystem service values of the study area were 40–60 billion yuan. According to Figures 912, it is obvious that regions with less farmland ecosystem service capability are more considerably affected by regions with greater capabilities. For regional division, compared to traditional Voronoi diagrams, the weighted Voronoi diagrams also consider the actual spatial distribution of farmland. Additionally, compared to the traditional Voronoi diagrams, which use color grading to represent ecosystem service values, the weighted diagrams reflect the progressive reduction over space and the nonuniformity of ecosystem services, more closely representing reality. Traditional ecosystem service value evaluation only considers the ecosystem service capabilities of cities, but not the direct interactions between cities, whereas weighted Voronoi diagrams can better explain these interactions. Construction ecosystems cause negative ecosystem service values and lower overall ecosystem service capabilities. Therefore, during calculation, all construction was assigned positive values. High-impact areas are those that are more greatly affected by construction ecosystems. The ecosystem service value losses of these areas are significantly higher than those of other areas. Areas with higher levels of urbanization and more extensive construction experience stronger reduction. In the study area, there are more forest ecosystems in the west and fewer in the east. Thus, areas with higher forest ecosystem service capabilities impose greater influences on areas with lower capabilities. Luoyang has vast forest ecosystems. Strong forest ecosystem service capability is the basic characteristic of the city’s ecosystem services; Luoyang’s ecosystem service value is 6.91, 7.32, 15.06, 119.42, 10.15, 27.3, 27.3, and 323.29 times of those of Xinxiang, Jiyuan, Jiaozuo, Kaifeng, Zhengzhou, Xuchang, Pingdingshan, and Luohe, respectively. Engulfing effects and compensation due to the forest ecosystem services of Luoyang gradually weaken with distance. However, the low-impact forest ecosystem service area of the city still has relatively significant effects on other areas. The weighted Voronoi diagrams, compared to the traditional Voronoi diagrams, have the advantage of intuitively presenting ecosystem service capabilities and their coverage over space. The coverage of water ecosystem services and the actual distribution of water ecosystems are relatively strongly correlated. The Yellow River runs through the central and the central-north regions of the Zhongyuan Urban Agglomeration. Rivers in Luoyang are widely distributed, which places Luohe in a low-impact water ecosystem service area. This distribution of high-impact areas in the northwest and low-impact areas in the northeast matches the satellite image interpretation results. Comparison between the satellite image interpretation (Figure 13) and ecosystem service coverage map reveals that high-, medium-, and low-impact ecosystem service areas are directly related to the distinct types of ecosystems in different cities, further confirming that weighted Voronoi diagrams are suitable for simulation of ecosystem service coverage and provide a new method and direction for these studies.

5. Discussion and Conclusion

5.1. Discussion

This study employs cities as central points for investigating ecosystem service coverage. This method considers an entire region and thereby presents the spatial distribution of ecosystem services, providing a new perspective for ecosystem service coverage research. Nevertheless, this method neglects nonuniform ecosystem distributions, ecosystem distribution within the region, and the interactions between ecosystems, which reduce the accuracy of ecosystem service coverage research. Thus, in future studies, better results can be obtained if fixed point studies can be conducted for different land use types.

Ecosystem service capabilities were classified using the Jenks natural breaks classification method. Although this method can partly reveal the basic features of research topics and data, it neglects the gradual reduction in service distance of service providers and the acceptance capabilities of service. Proper investigation of the gradual decrease of different ecosystem services over distance, and the acceptance rates of different service providers, allows the new method to be better applied to the actual production activities.

5.2. Conclusions

The total ecosystem service value of the Zhongyuan Urban Agglomeration in 2015 was 575.7 billion yuan. The total ecosystem service value of Luoyang was the highest, at 267.665 billion yuan; that of Luohe was the lowest, at only 13.881 billion yuan. The differences between regions are huge. Forest ecosystems have the greatest ecosystem service value per unit area—up to 32.639 thousand yuan per hm. Meanwhile, construction ecosystem service values significantly reduce the service values of different ecosystem types. Continuous expansion of construction ecosystems has decreased the ecosystem service value of Zhengzhou by 37.012 billion yuan, with the most remarkable decline seen in Kaifeng.

Compared to traditional Voronoi diagrams, weighted Voronoi diagrams are more suitable for studying ecosystem service coverage. Comparison between the two indicates that the spatial coverage of ecosystem services presented by the weighted Voronoi diagrams better fit actual conditions. Traditional Voronoi diagrams describe uniform ecosystem services and neglect the discrepancies in actual service coverage caused by regional differences. In weighted Voronoi diagrams, regions with higher ecosystem service capabilities partly engulf and compensate for regions with lower capabilities. These engulfing and compensating actions are inversely proportional to the distance between two regions. When two ecosystems are closer to each other, their interactions are stronger. The region with lower capability is more easily engulfed by the region with higher capability, and the ecosystem services of the former are enhanced. These actions partly contribute to the spatial and temporal distribution of ecosystems in the region and promote the union of different regions, thereby enhancing overall regional ecosystem service capability.

The ecosystem service coverage of farm ecosystem, forest ecosystems, water ecosystems, and construction ecosystems was studied separately, and the results compared to actual conditions. We found that the two agree relatively well, again proving that weighted Voronoi diagrams perform well in assessing ecosystem service coverage. The service coverage of different ecosystem types was also examined. We discovered that, for larger differences in ecosystem type and size between different regions, weighted Voronoi diagrams differ more significantly than their traditional counterparts. This result demonstrates that, for two regions whose ecosystem service differences are relatively large, traditional Voronoi diagrams do not match actual service coverage.

Data Availability

The readers can contact the authors to access the data in the manuscript.

Disclosure

The content does not represent the official position of the Chinese government and is entirely the responsibility of the authors.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This paper was supported by the National Natural Science Foundation of China (41601175 and 41771380). The Key Scientific Research Project of Henan Province (16A610001), supported by the Program for Innovative Research Team (in Science and Technology) in the University of Henan Province (16IRTSTHN012), the GDAS’ Project of Science and Technology Development (2017GDASCX-0101 and 2018GDASCX-0904), the Guangdong Innovative and Entrepreneurial Research Team Program (2016ZT06D336), the Science and Technology Program of Guangdong Province (2017B010117008), and the Science and Technology Program of Guangzhou (201604016047).