Next Article in Journal
Embedding Sustainable Mathematics Higher Education in the Fourth Industrial Revolution Era Post-COVID-19: Exploring Technology-Based Teaching Methods
Next Article in Special Issue
The Influence of Urbanization and Environmental Factors on the Financial Performance of Retail and Automotive Industries in China
Previous Article in Journal
The Role of the Forest Recreation Industry in China’s National Economy: An Input–Output Analysis
Previous Article in Special Issue
Co-Application of Sewage Sludge, Chinese Medicinal Herbal Residue and Biochar Attenuated Accumulation and Translocation of Antibiotics in Soils and Crops
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial and Temporal Evolution of Ecosystem Service Values and Topography-Driven Effects Based on Land Use Change: A Case Study of the Guangdong–Hong Kong–Macao Greater Bay Area

1
College of Forestry and Landscape Architecture, South China Agricultural University, Guangzhou 510642, China
2
Guangdong Rural Construction Research Institute, Guangzhou 510642, China
3
College of Architecture and Planning, Yunnan University, Kunming 650300, China
*
Authors to whom correspondence should be addressed.
Sustainability 2023, 15(12), 9691; https://doi.org/10.3390/su15129691
Submission received: 24 April 2023 / Revised: 12 June 2023 / Accepted: 14 June 2023 / Published: 16 June 2023
(This article belongs to the Special Issue Advances in Ecosystem Services and Urban Sustainability)

Abstract

:
The Guangdong–Hong Kong–Macao Greater Bay Area (GBA) is rich in natural and marine resources, and it is scientifically valuable to study the evolution patterns and driving mechanisms of the ecosystem service values (ESVs) of the GBA for the governance and conservation of its ecosystems. Based on the land use changes in the GBA from 2000 to 2020, the ESVs in the GBA were measured at the grid scale, and the Markov model was used to predict the ESVs in 2030; the calculated results were used to analyze the spatial and temporal variation characteristics of the ESVs during the 30-year period, while the driving role of the topographic factors on the ESVs is revealed through the construction of the geographically weighted regression model (GWR). The results show the following: (1) During the 20-year period, the area of arable land and water in the GBA fluctuated greatly, with the area decreasing year by year and shifting mainly into construction land; in terms of shifting the center of gravity of the land, and the center of gravity of the grassland and unused land shifted the greatest distance due to the expansion of construction land, with the center of gravity shifting westward as a whole. (2) The ecosystem services (ESs) in the GBA show obvious aggregation in the spatial distribution, with the total ESVs decreasing year by year. Among them, the areas with an increasing total value are mainly located in the cities of Zhaoqing and Huizhou in the GBA, accounting for 27%, and the areas with a decreasing total value year by year are mainly located in the dense urban areas in the central part of the GBA, accounting for 35%, and the area is increasing, indicating that the habitat quality is deteriorating, and the model prediction shows that the value of ecosystem services in 2030 have a decreasing trend under the development of the natural state. (3) Topographic factors have a significant influence on the ESVs, and in terms of spatial distribution, the areas with the strongest effect are distributed in the northwestern and northeastern parts of the GBA, and the main uses for the land are wood land, arable land, water and the area of the water–land intersection near the sea.

1. Introduction

As a strategic national development region, the Guangdong–Hong Kong–Macao Greater Bay Area (GBA) has a key geographical location and a wealth of natural resources. The sustainable development of the GBA plays a crucial role in the enhancement of regional ESVs as well as the construction and optimization of ecological patterns [1]. For the coastal areas with complex and fragile environments, impaired habitat quality will lead to a decline in the ESVs and negatively affect environmental security and socioeconomic development [2]. Therefore, strengthening the accounting of the ESVs and analyzing the influencing factors that affect the evolution of the ESVs will help with the governance and protection of ecosystems in the GBA [3].
The concept of ecosystems was first proposed by Tasley in the 1930s, marking the gradual refinement of the study of ecosystems; on this basis, SCEP first proposed the concept of ESs in the 1970s [4], and listed a number of ESs provided by the natural environment. Subsequently, Costanza deepened his research on ESVs and classified them into 17 categories, including gas regulation, climate regulation, water regulation, soil formation, nutrient cycling and recreation, based on four dimensions including production, basic functions, environmental benefits and recreation, which became the basis for the classification of subsequent scholars and greatly enriched the research results in this field. At the beginning of the 21st century, the UN Millennium Ecosystem Assessment redefined the concept of ESs and classified them as provisioning, regulating, cultural and supporting services in relation to human benefits, leading to a boom in research on the value of ecosystem services. In China, the study of ESs began in the 1990s when Liu Xiaodi introduced the concept of ESs to the country based on the research results of Daily [5]; subsequently, Ouyang Zhiyun conducted in-depth research on this concept, continuously enriching the theoretical findings on ecosystem services [6]. In the same era, Xie Gaoqi defined ES functions as products and services obtained directly or indirectly through ecosystem functions and classified them into three major categories, namely, production functions, ecological functions and recreational and leisure functions [7]; since then, domestic research on ES functions has been gradually enriched and diversified [8,9].
Early research on ESs at home and abroad focused on the valuation of ecosystem services. There are two main approaches to accounting; the first is the physical quality assessment method, which, in turn, includes the functional volume assessment method and the energy value assessment method. The material quality refers to the value of the final products or services obtained by humans directly or indirectly from the ecosystem [10]. Its principle is to obtain the total value of an ecosystem service from its functional volume and the unit price of that volume [11,12]; the functional volume of the relevant indicator is usually calculated using models such as InVEST and ARIES, which are commonly used. The advantages of the physical quality assessment method are that the results can objectively reflect the structural functions and ecological processes of the ecosystem, but the limitations are that there are uncertainties in the data acquisition and the calculation process is cumbersome, and there are differences in the units of measurement for the individual functional indicators, which makes it difficult to measure the value of ecosystem services with multiple functions. The second approach is the value–volume approach, which includes the monetary and value-equivalent approaches. The principle of this approach is to quantify ecosystem services using economic algorithms, including the direct market approach, the substitution market approach and the simulated market approach [13]. The analysis of the different ecosystem service valuation methods shows that each method has its own advantages and limitations, so the appropriate method can be weighed according to the purpose and focus of the study [1]. Based on changes in the social environment and the different methods of accounting for value, scholars measured the ESVs at different environmental scales and gradually built up a systematic research method for quantifying ESs [14,15,16,17], of which a highly representative research result was the principle and method of valuing ESs proposed by Costanza et al. in 1997 [18]; Costanza first accounted the ESVs at the global scale using this method. Subsequently, Chinese scholars such as Xie Gaodi combined the actual situation of ESs in China and, after a series of revisions, developed a table of value equivalent factors applicable to terrestrial ecosystems in China [19]. Later research directions integrate spatial techniques in value assessment to explore the coupling between the land use types and the ESVs. The number and rate of land type changes can directly reflect the intensity of the land type conversion [20,21,22], thus exploring the temporal and spatial evolution and variation of ESVs in the region from multiple dimensions.
At present, research on the factors influencing the evolution of ESVs has made great progress, mainly focusing on the relevance of physical geographic factors to ESVs [23,24,25,26], while regional topography, as one of the determinants of landscape patterns, has an important influence on the evolution of ESVs [27,28]. There are two main methods for coupling ESVs with topographic factors. The first is the direct coupling method, which is mainly calculated through the topographic position index; this index is a composite analysis of the elevation and slope, and is often used to quantify the spatial effects of land use on topographic gradients. The second is the model coupling method, which is mainly carried out through a spatial autocorrelation analysis, such as the Moran index and the geographically weighted models, among which the GWR is more advantageous in the analysis of coupling relationships, because it can express the data characteristics of spatial data at different locations and can reflect the weight size of local areas in the data [29,30]. Therefore, it is important to use spatial autocorrelation analysis to study the driving effect of ESVs in order to optimize the regional ESVs and build an ecological security pattern [31,32].
From the current research results, it can be seen that studies on ESVs in the GBA mostly focus on the relationship between urbanization and ESVs, but the depth of the studies are still insufficient [33,34,35], mainly in two aspects. Firstly, current studies on the spatial and temporal variation of regional ESVs seldom involve future-year scenarios, and at the same time, when analyzing the spatial and temporal variation of ESVs in regions with rapid urbanization, the coupling study of ESVs and land use changes is seldom involved; secondly, studies of the driving effects of ESVs within the GBA are less likely to involve topographic factors. Therefore, based on the remote sensing image interpretation maps of the GBA from 2000 to 2020, this study measures and predicts ESVs in the GBA in 2030 at the grid scale, analyzes its spatial and temporal variation characteristics, and introduces a model called GWR to reveal the driving effect of topographic factors on ESVs. The aim is to explore the following issues: (1) the spatial and temporal evolutionary characteristics of ESVs in the GBA; (2) the impact of land use on changes in the ESVs; and (3) the relationship between topographic factors and ESVs. These issues are explored to deepen the content of the study and to provide theoretical guidance for the GBA ecosystem to achieve regional ecological governance and integrated urban development.

2. Materials and Methods

2.1. Study Area Overview

The GBA belongs to one of the four major bay areas in the world, with a total surface area of about 56,000 square kilometers. In terms of administrative area, it mainly consists of the nine cities of Guangzhou, Shenzhen, Zhuhai, Foshan, Huizhou, Dongguan, Zhongshan, Jiangmen and Zhaoqing in the Guangdong Province and the special administrative regions of Hong Kong and Macau (Figure 1). The study area has a mild and humid climate, with a dense network of rivers and water, and is distributed in the mountainous forests of Mafeng Mountain–Baiyun Mountain, Gudou Mountain–Wuguishan Mountain–Phoenix Mountain, Daling Mountain–Yangtai Mountain–Tanglang Mountain, etc., which form an inter-urban ecological transition zone. The study area also has a wide distribution of forests and agricultural land, with arable land scattered throughout the bay area; wetland resources are abundant, with water entering the rivers located in the four cities of Foshan, Zhuhai, Shenzhen and Zhongshan [36,37].

2.2. Data Sources and Initial Data Processing

The remote sensing image source data selected for this study were obtained from the Resource and Environment Science and Data Centre of the Chinese Academy of Sciences (http://www.resdc.cn/, accessed on 1 July 2021), with three periods including 2000, 2010 and 2020, and the spatial resolution of the data is 30 m. The topographic data were obtained from the geospatial data cloud database (http://www.gscloud.cn/, accessed on 1 July 2021), with the type ASTER GDEM. The relevant statistical data were mainly obtained from the Guangdong Statistical Yearbook, the Hong Kong Statistical Yearbook and the Macao Statistical Yearbook during the study period. The collected remote sensing data were corrected and classified using ArcGIS software, and according to the previous research results, the land types were classified into six categories, such as construction land, arable land, grassland, woodland, water and unused land [38]; the relevant yearbook data were collated, summarized and analyzed, and processed using SPSS software.

2.3. Research Methods

The research idea and process are shown in Figure 2.

2.3.1. Analysis of Land Use Change

The land use transfer matrix can concretely reflect the characteristics of the land use structure. It reflects the transition relationship between a certain land type within a certain time interval, and calculates the transition process of land state from moment t to moment t + 1 using Equation (1) [39,40]; the direction and distance of the land center of gravity shift can reflect the natural economic conditions as well as the quality of land in the region [41,42], and the coordinate change of the land use type center of gravity can be calculated using Equation (2) as follows:
S ab = S 11 S 12 S 1 n S 21 S 22 S 2 n S n 1 S n 2 S nn
In the equation, Sab represents the land use status at the beginning and end of the study, and n represents the number of land use types.
X a = i = 1 n ( C a i × X i ) i = 1 n C a i Y a = i = 1 n ( C a i × Y i ) i = 1 n C a i
In the equation, Xa and Ya represent the latitude and longitude coordinates of the center of gravity of a land type in year a, Cai represents the area of the i patch in a land type in year a, and Xi and Yi represent the latitude and longitude coordinates of the center of gravity of the i patch in a land type.

2.3.2. Ecosystem Service Value at the Grid Scale

The grid scale can quantify the spatial heterogeneity of the land use types in the study area and express their information adequately, which is an effective evaluation unit for quantifying the spatial and temporal evolution of the land use [43]. In this study, 1 km × 1 km, 3 km × 3 km, 5 km × 5 km, 10 km × 10 km, 15 km × 15 km and 20 km × 20 km grids were constructed as pre-selected evaluation units based on previous research results. Considering the scale and shape characteristics of the study area, in order to make the grid cells not only have certain regional characteristics, but also contain sufficient research information, and at the same time, be able to reflect the spatial variation of ESVs in the study area more freshly, we referred to a large amount of the literature relevant to our study and compared the spatial variation effect of different grid calculation results, and finally chose 3 km × 3 km as the base measurement cell (Figure 3). The calculation of ESVs referred to the equivalence factor table modified by Xie Gaodi according to the actual situation in China; in addition, considering the differences in grain prices in different regions, the average grain price of farmland in the Guangdong Province was chosen as the average grain price in this study, while the ESV coefficients for construction land referred to the coefficients proposed in the table for the value of each ecosystem service per unit area of the different types of terrestrial ecosystems in China to make corrections [44,45]. In terms of the value per unit area, ESVs at the grid scale were corrected using the coefficient proposed by Xie Gaodi (1.40 for Guangdong Province) and using the vegetation cover [46], resulting in ESVs per unit area in the GBA (Table 1), with the following revised equation:
F V = N D V I N D V I min N D V I max N D V I min F a = f ai / f E = E a × F a
In the equation, FV is the vegetation index, NDVI is the standard value of vegetation cover, fa denotes the sum of the FV values of the ath grid cell, f is the mean value of FV, Ea and E refer to ESVs before and after the revision and Fa denotes the vegetation cover revision factor of the ath grid cell.

2.3.3. Markov Model

Markov chain represents the sequential state of the development process of things, which is discrete in time and state, and the development process of things is only linked to adjacent things. Markov model can be used on this basis to predict the future state of things in the same time interval according to the development process and state of things [47,48,49]. In this study, the Markov model is used to predict and analyze the change in land use types in the GBA in 2030 and estimate the ESVs in the study area based on the prediction results. The relevant calculation equation is as follows:
S ( a + 1 ) = P ij × S ( a )
In the equation, S(a), S(a + 1) represent the states of the ground class at moments a, a + 1; Pij represents the state transfer matrix.

2.3.4. Geographically Weighted Regression Model

Geographically weighted regression model (GWR) is a local regression model that assists the original spatial data in determining the coordinate location parameters [50], which is based on the principle of independent linear regression calculations at all point locations, thereby allowing for the expression of data characteristics of spatial data at different locations, perfectly compensating for the shortcomings of the global regression model, while reflecting the magnitude of the weights of local areas in the data [51,52]. In this study, the relationship between topography and slope data and ESVs is analyzed comprehensively by constructing a GWR model, in which topography and slope data are used as independent variables, and the ESVs are used as dependent variables, from which the corresponding explanatory variable coefficients (predicted) and local R2 data can be calculated. The Y value represents the predicted estimate of the dependent variable, and the closer the Y value is to the ESVs, the better the fit between the two; the local R2 means the fit of the local regression model to the Y value, with larger values indicating a higher model accuracy [53]. The relevant calculation equation is as follows:
Y i = β 0 ( μ a , ν a ) + a = 1 k β k ( μ a , ν a ) X ak + ε a
In the equation, (ua, va) represents the spatial coordinate point a; Ya and Kak are the dependent variable Y and the set of sub variables Xk of the measured values of the spatial position; k is the number of independent variables; β0 (ua, va) is the spatial position of the constant term (ua, va); βk (ua, va) is the value of the continuous function βk (ua, va) at point a and εa represents the random error term.

3. Results

3.1. Analysis of Spatial and Temporal Land Use Change Characteristics

3.1.1. Land Use Transfer Matrix

From the land use transfer matrix of the study area (Table 2), the urban boundaries within the GBA expanded continuously from 2000 to 2010, encroaching on the surrounding arable land, woodland and water. The areas of arable land, water and woodland transferred out of the GBA during the 10-year period were 2636 km2, 1098 km2 and 992 km2, respectively, while the largest amount of construction land was transferred in, accounting for 43% of the original area, which is mainly due to the fact that the urban master plan expanded the use of reserve land resources such as mudflats and coastal zones, while the construction of infrastructure also accelerated the change in the construction land. From 2010 to 2020, arable land, waters and woodland remained the largest land types transferred out, with areas of 3459 km2, 2141 km2 and 1937 km2, respectively, while construction land was transferred in at 4205 km2. In general, the area of arable land in the GBA fluctuated due to the national policy of returning farmland to forests and the continuous reduction in the primary industry over the past 20 years; the area of woodland changed less, while the area of land transferred for construction reached 51% of the original area over the past 20 years, which shows that the rapid development of the GBA has put the ecological land under pressure (Figure 4).

3.1.2. Shift in Land Use Focus

From 2000 to 2020, the center of gravity of the grassland, unused land and water areas in the GBA shifted significantly. The main reason for this is the expansion of the urban land area and the excessive artificial reclamation, which led to the occupation of ecological land such as ponds, dike-ponds and lakes; the center of gravity of the grassland moved 22,897 m to the northwest, and the center of gravity of the urban construction land moved 3602 m to the north and then 936 m to the southeast. The results of the analysis show the impact of the reorientation of urban development in different periods on various types of land, and that the shift in the center of gravity of the land use indicates a change in the structure of the land use in the region, as the environmental carrying capacity and self-regulating ability of the GBA are affected, resulting in a change in the ESVs (Table 3).

3.2. Analysis of the Spatial and Temporal Evolution of Ecosystem Service Values

3.2.1. Time Series Changes in Ecosystem Service Values

In terms of the changes in the total ESVs, there is an overall decreasing trend from 2000 to 2020; in terms of the magnitude of change, there is a greater change in the waters, woodlands and croplands, thus indicating a continuous decline in the environmental quality of the region (Table 4).
In terms of the change in the ESVs of the different land types, hydrological regulation, biodiversity and soil conservation have the highest contribution ratings, with all three contributing more than 44.8%, while the least contribution is made by recreation and food production, which account for less than 10%. Specifically, soil conservation, food production and biodiversity maintenance are the main service functions of arable land; the main service function of woodland is biodiversity; the main service function of grassland is soil conservation, but the overall contribution of grassland is low due to its small footprint and the main service function of the watershed area is hydrological regulation and waste disposal. It is worth mentioning that the ESVs of the watershed area declined significantly by 139.37 × 108 CNY from 2010 to 2020, which is mainly due to the destruction of watershed habitats as a result of the dike-pond reclamation project and urban expansion, with the watershed area decreasing year by year. Overall, the ESVs in the study area show a decreasing trend due to the expansion of urban construction sites and the increasing negative impacts of human disturbances such as industrial waste emissions, domestic waste and vehicle exhaust (Figure 5).

3.2.2. Spatial Changes in Ecosystem Service Values

The results of ESVs in the GBA are classified into the following five standard classes: very low(ESV < 15 × 106 CNY/hm2), low (ESV of 15 × 106–30 × 106 CNY/hm2), medium (ESV of 30 × 106–45 × 106 CNY/hm2), high (ESV of 40 × 106–60 × 106 CNY/hm2) and very high (ESV > 60 × 106 CNY/hm2). In terms of the spatial distribution of the regional ESV class (Figure 6), the areas with a very high ESV class in the GBA are mainly located in the northeastern part of Zhaoqing, the northeastern part of Huizhou and the western and southern part of Jiangmen, which is mainly due to the rich woodland resources and good ecological substrate in the area. The very low ESV classes are found in the central part of the region, which is mainly due to the concentration of urban agglomerations in the region and the high level of disturbance from human activities, which accelerated changes in vegetation, climate and other related factors and led to a continuous decrease in the quality of water areas and water resources, with a large number of ponds and other water areas being encroached upon, as well as the negative impact of urban sewage on rivers and marine ecosystems, ultimately leading to a continuous decrease in the quality of habitats in the central urban agglomerations. Specifically, from 2000 to 2010, the area of low ESVs in the GBA expanded significantly. By 2020, low ESV areas became the dominant class type, accounting for 27.24% of the area, and still show an increasing trend, indicating that the ESVs in the GBA have been decreasing over the past 20 years. From the grid cells, the ESV cells in the GBA are mainly medium and low grade, and the high ESV area shows a small increase between 2000 and 2020, from 1534 cells to 1740 cells, which is mainly due to the related cities in the northwest of the GBA to increase the protection and restoration of forest land, slowing down the trend of reducing the forest land area. Among them, the proportion of high-grade cells in the Hong Kong and Macao special administrative regions is 3%, and it shows a small increase over 20 years, which is mainly due to the fact that the forest area in the Hong Kong and Macao special administrative region has increased by about 85 km2 in 20 years, which contributed to the rebound of ESVs to a certain extent (Table 5).
In terms of the spatial changes in the ESV class, the overall ESV rank in the GBA remained stable from 2000 to 2010, with 34% of the regional ranks increasing, and 29% of them decreasing, mainly because the economic growth in the GBA was still slow at the beginning of the 21st century, with less disturbance to the ecological environment and land. Among them, the nine cities in the Pearl River Delta have a high proportion of areas with no significant changes and small increases in the ES levels, while the Hong Kong and Macao special administrative regions have an increasing trend in the ES function levels, mainly because the ecological restoration of Mai Po Nature Reserve in Hong Kong was effective, and the restored mangrove and mudflat ecosystems can provide a better ecological benefit. During 2010–2020, the overall ES level of the GBA declined significantly, among which the ES level of nine cities in the Pearl River Delta declined significantly, which is mainly due to the continuous adjustment in urban construction and industrial structure, and the accelerated growth rate of economic development in Zhaoqing, Huizhou, Zhongshan and other cities. A total of 33% of the regional ESV levels in the Hong Kong and Macau special administrative regions showed a declining trend, among which the level in the coastal areas declined significantly, which is mainly due to the prosperity of land reclamation projects in coastal areas, the decline of the retention rate of the natural shoreline, and the large changes in the regional ecological environment quality and ecological pattern (Table 6). Overall, 38% of the regional ESV class in the GBA did not change significantly over the 20-year period, the proportion of decreasing regional cells was 35%, the proportion of increasing regional cells was 27% and the number of decreasing regional cells continuously increased, indicating that the regional habitat quality is deteriorating (Figure 7).
In terms of the regional spatial correlation, the ESVs in the GBA show significant spatial aggregation, with “high-high” and “low-low” indicating the proximity of high-value and low-value areas of ESVs. The global Moran index of the study area shows an increasing trend from 2000 to 2020, indicating the increasing aggregation of ESVs in the GBA. Specifically, from 2000 to 2010, the “high-high” area continued to expand southward, which is mainly due to the obvious transfer of woodland to the south of the GBA, the high vegetation cover of wetlands and river corridors, the better maintained landscape pattern and the more stable ecological land function structure. The “low-low” region is shifting from scattered distribution to patchy development, such as in Foshan, Zhongshan, Jiangmen and other cities, mainly because the concentrated development of construction land increased the intensity of land use and increased pressure on the ecological environment within the city. The Hong Kong and Macao special administrative regions also belong to the “low-low” region, mainly because of the limited space for urban development in the region, so a large number of land reclamation projects were carried out, encroaching on natural ecological land such as coastal shorelines and water areas, and the increase in the proportion of human-made land surface reduced the ESVs. From 2010 to 2020, the “high-high” zone continued to develop to the south, and the “low-low” zone is dominated by areas of concentrated urban agglomerations, which generate low ESVs due to a low intra-urban landscape diversity (Figure 8).

3.2.3. Forecasting the Ecosystem Service Values

By running the Markov model under the natural development scenario, the predicted ESV in 2030 reaches 1853.921 × 108 CNY, a decrease of 64.839 × 108 CNY compared to 2020, which is mainly due to the decreasing areas of watershed, woodland and grassland from 2000 to 2020, while the area of arable land shows stability in the context of policy protection. According to this pattern of development, the areas of watershed, woodland and grassland will keep decreasing in the future, resulting in a decreasing ESV in the forecast year 2030. At the same time, it can be observed that the future decline of water resources will remain a serious problem, and therefore, scientific and rational planning are needed to reduce damage to water resources and the near-coastal environment (Table 7).

3.3. Analysis of the Topographic Factor-Driven Effects of Ecosystem Service Value

The elevation and slope were selected as the topographic characterization factors of the study area with reference to the relevant literature. The results of the geographically weighted regression model (GWR) show that R2 > 0.8 and local R2 > 0.8, indicating that the model operation results meet the accuracy requirements. The results show that the local R2 of the data from 2000 to 2020 changed significantly, showing a decreasing trend and then an increasing trend, indicating that the influence of topographic factors on ESVs generally show an increasing trend. The regression parameters of the topographic factors (Table 8) show that the local R2 values are above 0.85 in all years, indicating a significant correlation between the topographic factors and the ESVs.
In terms of spatial distribution, the strongest influence of topography on ESV is found in the northwestern and northeastern parts of the GBA, where the correlation with the land use data shows that woodland, arable land and water are the main land types, while the influence of the land–water interface in the south is also evident. Elevation and slope have a direct influence on plant growth, which determines the vegetation cover of the area and, thus, has a further influence on the ESVs provided (Figure 9).
From the perspective of construction land, apart from parks and green belts within the city, most areas have poor green infrastructure, and grey infrastructure, human housing and other construction land directly change the ecological landscape pattern of the land, which shows that human factors have a very strong influence on the ESVs. From 2000 to 2020, the trend of the overall pattern of the GBA clusters into patches, which shows that the influence of topographic factors on ESVs is gradually increasing.

4. Discussion and Conclusions

4.1. Discussion

4.1.1. Land Use Change Significantly Affects Ecosystem Service Values

Land is fundamental; the state and structure of the ecosystem are dependent on the land use approach, and the excessive repurposing of ecological land areas will cause the center of gravity of the land to continuously migrate, changing the landscape pattern of the region to a large extent, while affecting the supply of ESs. To achieve the sustainable development of the regional ecosystem, we must first pay attention to the adjustment of the land use structure; the regulation capacity and carrying capacity of the natural environment need reasonable planning of land use [54]. For example, in the areas around construction sites, the construction of nature reserves, the ecological restoration of rivers and lakes and the management of land reclamation projects can be strengthened to safeguard biodiversity, soil conservation and other ecological functions so as to safeguard the carrying capacity and self-regulating ability of the natural environment as a whole and enhance ESVs [55,56].

4.1.2. Spatial and Temporal Variation in Ecosystem Service Values and Strategies for Optimization

The spatial and temporal variations in ESVs in the GBA were significant from 2000 to 2020. Overall, in terms of the spatial distribution, the high-value areas are concentrated in the northwest, northeast and southwest regions where ecological land is widely distributed, while the low-value areas are mainly concentrated in the urban contiguous development areas and some coastal areas in the central part of the GBA, and the area of the low-value areas is still increasing. Specifically, in order to slow down the rate of diminishing ESVs, we should take the initiative to find entry points and explore optimization strategies in terms of human social life [57,58,59]. In terms of water resources, while the deep treatment and re-circulation of domestic fecal water and industrial wastewater should be carried out to achieve less or no discharge and improve the utilization rate of water resources, the rational use of water resources should be strengthened, water supply and drainage systems should be improved, land use planning should be carried out scientifically and reasonably and damage to water resources and the near-shore environment should be reduced [60,61]. In terms of energy use, the characteristics of each city in the GBA should be addressed, the upgrading of equipment for renewable and clean energy and the development of new energy sources should be enhanced to provide a reasonable and effective supplement [62]. In terms of solid waste treatment, while improving the harmless treatment of waste, the whole process of management and control should be implemented from all aspects of waste generation, collection, transportation, reuse and treatment to avoid causing adverse effects on the ecological environment and achieve control of waste generation.

4.1.3. Topographic Factors and Ecosystem Service Values

Topography has a direct impact on the growth state and diversity of plants, while the type, growth and distribution of plants have a critical impact on key regional ESs, such as soil and water conservation and water harvesting functions. The above analysis shows that there is a significant correlation between topographic factors and ESVs in the GBA, with an increase from 2000 to 2020. In the northwestern and northeastern parts of the region, the influence of topography on ES is particularly significant in areas such as woodland, arable land and the land–water interface, so that the protection and optimization of these areas is conducive to the sustainable development of the region, to the optimization of ESVs and to the construction of ecological security patterns across the region [63].

4.1.4. Sustainability in the Guangdong–Hong Kong–Macao Greater Bay Area

The results of the valuation of ESs in the GBA show that the region is currently facing significant ecological degradation and intensified resource and environmental constraints. At the same time, due to the GBA being located in the subtropical monsoon climate zone, its special geographical location and climatic conditions expose it to ecological risks such as flooding, storm surges and sea level rise, which pose challenges to regional ecological restoration and the construction of ecological security patterns. In the face of the continuous decline of ecological and environmental quality in the GBA, it is necessary to gradually promote the formation of green and low-carbon production and lifestyle and urban construction and operation mode to realize the sustainable development of the GBA. Therefore, the future development planning of the GBA should pay attention to the delineation of ecological protection zones and the control of urban expansion boundaries, and resist ecological risks such as storm surges and sea level rise by building a resilient multi-objective, cross-habitat and land–sea integrated regional ecological security pattern, while promoting the ecological restoration of coastal wetlands and coastal zone maintenance to construct a coastal security barrier in the GBA to achieve sustainable development [41].

4.2. Conclusions

Based on the land use changes in the GBA from 2000 to 2020, this study measures the ESVs in the GBA at the grid scale, predicts the ESVs in 2030 and analyzes its spatial and temporal evolution characteristics and the influence of topographic factors in ESVs. This will provide scientific guidance for optimizing ESVs.
(1)
Over the past 20 years, in terms of the shift in the land area, the most obvious fluctuations are in the areas of water, arable land and construction land; the northwestern and northeastern areas with high vegetation cover gradually decreased, while the impervious area of the urban agglomerations in the central area is increasing, and the ecological land around the urban agglomerations, such as arable land and wetlands, is decreasing. In terms of the shift in the center of gravity of the land, the center of gravity of all types of land has shifted to different degrees, with the center of gravity of construction land shifting northward, while the center of gravity of grassland and unused land shifted westward, and the center of gravity of forest land has shifted northwestward, with significant contradictions between the survival of ecological land and the expansion of urban construction land.
(2)
From an overall perspective, the ESVs in the GBA from highest to lowest are the following: regulating services > supporting services > provisioning services > cultural services. Hydrological regulation, biodiversity and soil conservation are the three items that contribute the most to the values of individual ESs. In terms of temporal changes, the ESs in the GBA showed a decreasing trend in the overall and individual values over the 20-year period. In terms of the spatial changes in rank, the area with the lowest ESV rank is located in the dense urban area in the central part of the GBA, accounting for 35% of the total area and increasing in size, indicating the deteriorating quality of the habitat; the area with a relatively high ESV rank is located in the city of Zhaoqing, northeastern Huizhou and Jiangmen in the GBA, accounting for 27% of the total area. The Markov model predicts that the ecosystem service value in 2030 shows a decreasing trend under the development of a natural state.
(3)
The spatial distribution of ESVs in the GBA is aggregated, and there is a regional adjacency between the high-value and low-value areas of ESVs. From 2010 to 2020, the “high” area continues to develop to the south, and the scope of the low-value aggregation area is also expanding. The construction land in the low-value area is increasing, and the population and traffic pressure are becoming bigger and bigger, which has a great impact on the urban green space system and reduces the anti-disturbance ability of the inner-city green space ecosystem.
(4)
There is a significant correlation between topographic factors and ESV, where topographic factors have a strong influence on ESVs mainly in the northwestern and northeastern parts of the GBA, and the coupling analysis with the land use data show that the main land types in this area are forest land, arable land and water. At the same time, the influence of the water–land intersection zone in the south is also obvious; the weaker influence is mainly located in the central part of the study area and the border area, and the main land types are construction land, dike-ponds and forest land.

Author Contributions

Conceptualization, H.L., Y.H. and Y.Z.; methodology, H.L., Y.H. and Y.Z.; software, Y.H. and Y.Z.; validation, Y.H. and S.W.; formal analysis, H.L., Y.H. and S.W.; investigation, W.G. and Y.L.; resources, J.W.; data curation, Y.H., Y.Z., J.W. and Y.L.; writing—original draft preparation, H.L. and Y.Z.; writing—review and editing, H.L., Y.H., S.W., Y.Z. and Q.X.; visualization, Y.H., Q.X., X.Z., K.Y., Q.H. and L.L.; supervision, H.L., Y.H. and W.L.; project administration, H.L.; funding acquisition, H.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a National Natural Science Foundation Project, grant number 52078222, and the Guangdong Provincial Education Department’s Key Scientific Research Projects of the year 2020, grant number 2020ZDZX1033.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Liu, Y.; Yan, Y.; Li, X. An Empirical Analysis of an Integrated Accounting Method to Assess the Non-Monetary and Monetary Value of Ecosystem Services. Sustainability 2020, 12, 8296. [Google Scholar] [CrossRef]
  2. Li, X.; Zhou, J.; Huang, Y.; Huang, M. Understanding the Guangdong-Hong Kong-Macao Greater Bay Area from the perspective of megacity region. Prog. Geogr. 2018, 37, 1609–1622. [Google Scholar]
  3. Feng, S.; Fan, F. Spatiotemporal changes of landscape pattern using impervious surface in Guangdong-Hong Kong-Macao Greater Bay Area, China Hong Kong-Macao Greater Bay Area, China. J. Appl. Ecol. 2018, 29, 2907–2914. [Google Scholar]
  4. SCEP. Man’s Impact on the Global Environment: Study of Critical Environmental Problems; MIT Press: Cambridge, MA, USA, 1970. [Google Scholar]
  5. Liu, X. Ecosystem Services. Environ. Her. 1998, 1, 44. [Google Scholar]
  6. Ouyang, Z.; Zheng, H. Ecological mechanisms of ecosystem services. Acta Ecol. Sin. 2009, 29, 6183–6188. [Google Scholar]
  7. Xie, G.; Xiao, Y.; Zhen, L.; Lu, C. Study on ecosystem services value of food production in China. Chin. J. Eco-Agric. 2005, 13, 10–13. [Google Scholar]
  8. Yan, S.; Wang, X. Advances of Research on Ecosystem Service. Chin. J. Ecol. 2002, 21, 61–68. [Google Scholar]
  9. Aodeng, G.; Liang, Y.; Han, G. Overview of grassland ecosystem service functions and their eco-economic value. Inn. Mong. Prataculture 2004, 16, 46–47. [Google Scholar]
  10. Li, L.; Wang, X.; Luo, L.; Yi, X.; Zhao, Y.; Zhao, Y.; Bachagha, N. A systematic review on the methods of ecosystem services value assessment. Chin. J. Ecol. 2018, 37, 1233–1245. [Google Scholar] [CrossRef]
  11. Kareiva, P.; Marvier, M. Conserving biodiversity coldspots-Recent calls to direct conservation funding to the world’s biodiversity hotspots may be bad investment advice. Am. Sci. 2003, 4, 220–224. [Google Scholar]
  12. Robertson, G.P.; Swinton, S.M. Reconciling agricultural productivity and environmental integrity: A grand challenge for agriculture. Front. Ecol. Environ. 2005, 1, 38–46. [Google Scholar] [CrossRef]
  13. Li, L.; Lin, W.; Ren, C.; Xv, D. Comparative Study on Two Assessment Methods of Ecosystem Service Value—A Case Study of Hangzhou Bay Area. Res. Soil Water Conserv. 2022, 29, 228–234. [Google Scholar]
  14. Zhou, Y.; Yi, L.; Kuang, H.; He, Y.; Yang, D.; Zhou, L. Study on the Changes and Ecological System Service Values in Guangdong-Hong Kong-Macao Greater Bay Area. Environ. Prot. 2019, 47, 56–60. [Google Scholar]
  15. Wang, S.; Huang, L.; Xv, X.; Xv, S. Spatial and temporal evolution of ecosystem services and its trade-offs and synergies in Guangdong-Hong KongMacao Greater Bay Area. Acta Ecol. Sin. 2020, 40, 8403–8416. [Google Scholar]
  16. Lin, M.; Zhou, R.; Zhong, L. Research on the changes of ecosystem services in Guangdong-Hong Kong-Macao Greater Bay Area based on the change of landscape pattern. J. Guangzhou Univ. (Nat. Sci. Ed.) 2019, 18, 87–95. [Google Scholar]
  17. Li, J.; Wang, J. Identification, classification, and mapping of coastal ecosystem services of the Guangdong, Hong Kong, and Macao Great Bay Area. Acta Ecol. Sin. 2019, 39, 6393–6403. [Google Scholar]
  18. Costanza, R.; d’Arge, R.; De Groot, R.; Farber, S.; Grasso, M.; Hannon, B.; Limburg, K.; Naeem, S.; O’neill, R.V.; Paruelo, J.; et al. The value of the world’s ecosystem services and natural capital. Nature 1997, 387, 253–260. [Google Scholar] [CrossRef]
  19. Xie, G.; Zhang, C.; Zhang, L.; Chen, W.; Li, S. Improvement of the Evaluation Method for Ecosystem Service Value Based on Per Unit Area. J. Nat. Resour. 2015, 30, 1243–1254. [Google Scholar]
  20. Wang, W.; Guo, H.; Chuai, X.; Dai, C.; Lai, L.; Zhang, M. The impact of land use change on the temporospatial variations of ecosystems servicesvalue in China and an optimized land use solution. Environ. Sci. Policy 2014, 44, 62–72. [Google Scholar] [CrossRef]
  21. Xie, W.; Fu, W.; Yang, D.; Liu, J. Spatio-temporal Evolution and Simulation Forecast of Ecosystem Service Value in Jiangsu Province Based on Land Use Changes. Areal Res. Dev. 2022, 41, 126–158. [Google Scholar]
  22. Ma, S.; Shi, M.; Yang, G.; Xv, X.; Yin, J. Analysis on spatiotemporal Change of Land Use Based on GIS Technology—Taking Xingjiang Tarim Basin as an Example. Res. Soil Water Conserv. 2013, 20, 177–181. [Google Scholar]
  23. Hu, X.; Hong, W.; Wu, C.; Qi, X. Spatiotemporal differentiation ofland eco-system service in Fuzhou City. China Environ. Sci. 2013, 33, 881–888. [Google Scholar]
  24. Wang, X.; Yan, F.; Su, F. Impacts of Urbanization on the Ecosystem Services in the Guangdong-Hong Kong-Macao Greater Bay Area, China. Remote Sens. 2020, 12, 3269. [Google Scholar] [CrossRef]
  25. Wu, L.; Sun, C.; Fan, F. Estimating the Characteristic Spatiotemporal Variation in Habitat Quality Using the InVEST Model—A Case Study from Guangdong-Hong Kong-Macao Greater Bay Area. Remote Sens. 2021, 13, 1008. [Google Scholar] [CrossRef]
  26. Liu, Y.; Hou, X.; Li, X.; Song, B.; Wang, C. Assessing and predicting changes in ecosystem service values based on land use/cover change in the Bohai Rim coastal zone. Ecol. Indic. 2020, 111, 106004. [Google Scholar] [CrossRef]
  27. Jiang, T. Thesis for Master Degree of Science Change of Farmland Ecosystem Services Values and Its Driving Forces in Northwestern Arid Area—A Case Study on Minqin County Tongchang JIANG. Master’s Thesis, Lanzhou University, Lanzhou, China, 2010. [Google Scholar]
  28. Yao, C.; Zhu, H.; Lv, X.; Liu, Y. Study on theImpact of Socio-economic Driving Factors of Land Use Change on the Ecosystem Service Values in Fujian Province. J. Nat. Resour. 2009, 24, 225–233. [Google Scholar]
  29. Chen, Y.; Li, J.; Xv, J. The influence of socio-economic factors on the value of ecological services in Hubei Province based on GWR. China Land Sci. 2015, 29, 89–96. [Google Scholar]
  30. Li, Q.; Li, T.; Yang, M.; Ying, W. Spatiotemporal variation of ecosystem services value based on gradient analysis in Wuhan. Acta Ecol. Sin. 2017, 37, 2118–2125. [Google Scholar]
  31. Yao, X.; Zhou, D.; Huang, W.; Liu, X. Spatial and temporal effects of ecosystem service values on topographic gradients in northwest Qianxi. Agric. Technol. 2022, 42, 93–100. [Google Scholar]
  32. Zhang, T.; Qu, Y.; Liu, Y.; Yan, G.; Foliente, G. Spatiotemporal Response of Ecosystem Service Values to Land Use Change in Xiamen, China. Sustainability 2022, 14, 12532. [Google Scholar] [CrossRef]
  33. Liu, Z.; Wang, S.; Fang, C. Spatiotemporal evolution and influencing mechanism of ecosystem service value in the Guangdong-Hong Kong-Macao Greater Bay Area. Acta Geogr. Sin. 2021, 76, 2797–2813. [Google Scholar]
  34. Wu, J.; Jiang, W.; Huang, Y. Spatial-Temporal Changes of Ecosystem Service Value in the Guangdong-Hong Kong-Macao Greater Bay Area. South Archit. 2022, 6, 52–61. [Google Scholar]
  35. Hou, J.; Qin, T.; Liu, S.; Wang, J.; Dong, B.; Yan, S.; Nie, H. Analysis and Prediction of Ecosystem Service Values Based on Land Use/Cover Change in the Yiluo River Basin. Sustainability 2021, 13, 6432. [Google Scholar] [CrossRef]
  36. Chunshan, Z.; Luo, L.; Shi, C.; Wang, J. Spatio-temporal Evolutionary Characteristics of the Economic Development in the Guangdong-Hong Kong-Macao Greater Bay Area and Its Influencing Factors. Trop. Geogr. 2017, 37, 802–813. [Google Scholar]
  37. Gu, S.; Cao, D.; Yang, M. The Strategic Thinking of Building Innovation Ecosystem in Guangdong-Hong Kong-Macao Greater Bay Area. Chin. Soft Sci. 2018, 4, 1–9. [Google Scholar]
  38. Zhu, M.; Fu, X. Spatial-temporal Evolution of Urban Land Use Efficiency in the Guangdong-Hong Kong-Macao Greater Bay Area. Trop. Geogr. 2017, 37, 814–823. [Google Scholar]
  39. Wang, F.; Gao, J.; Shao, H.; Zhang, T.; Zhang, Y.; Xv, X.; Zhao, C.; Wang, H. Response of ecosystem service values to land use change based on GIS and ecological compensation in Loess Plateau. Sci. Soil Water Conserv. 2013, 11, 25–31. [Google Scholar]
  40. Niu, X. Research on Distributed Assessment and Quantitative Compensation of Forest Ecological Benefit—Case Studies with Guangdong and Liaoning Provinces. Master’s Thesis, Beijing Forestry University, Beijing, China, 2012. [Google Scholar]
  41. Zhang, J.; Li, Q.; Wu, X.; Zhang, C.; Wang, J.; Wu, G. Evolution of the ecosystem services value and carrying capacity in the Guangdong-Hong KongMacao Greater Bay Area based on land use changes. Acta Ecol. Sin. 2021, 41, 8375–8386. [Google Scholar]
  42. Chang, X.; Li, X.; Diao, H. Analysis of spatial and temporal changes in land use and shift of focus in mining areas. J. Inn. Mong. Agric. Univ. (Nat. Sci. Ed.) 2021, 42, 41–48. [Google Scholar]
  43. Yang, G.; Lv, K.; Li, F. Spatial and Temporal Correlation Analysis of Land Use Change and Ecosystem Service Value in Nanchang City Based on Grid Scale. China Land Sci. 2022, 36, 121–130. [Google Scholar]
  44. Shi, Y.; Wang, R.; Huang, J.; Yang, W. Analysis of spatial and temporal variation in terrestrial ecosystem service functions in China. Sci. Bull. 2012, 57, 720–731. [Google Scholar] [CrossRef] [Green Version]
  45. Zhang, L.; Yu, X.; Jiang, M.; Xue, Z.; Lu, X.; Zou, Y. A consistent ecosystem services valuation method based on Total Economic Value and Equivalent Value Factors: A case study in the Sanjiang Plain, Northeast China. Ecol. Complex. 2017, 29, 40–48. [Google Scholar] [CrossRef]
  46. Xie, G.; Lu, C.; Cheng, S. Advances in global ecosystem service valuation research. Resour. Sci. 2001, 6, 5–9. [Google Scholar]
  47. Xv, Y.; Zhou, Y.; He, Y.; Kuang, H.; Yang, D. Study on the prediction of ecological system service values based on the changes of land use in Guangdong-Hongkong-Macau Greater Bay Area. Ecol. Sci. 2020, 39, 207–212. [Google Scholar]
  48. Luo, D.; Zhang, W. A comparison of Markov model-based methods for predicting the ecosystem service value of land use in Wuhan, central China. Ecosyst. Serv. 2014, 7, 57–65. [Google Scholar] [CrossRef]
  49. Gao, X.; Wang, J.; Li, C.; Shen, W.; Song, Z.; Nie, C.; Zhang, X. Land use change simulation and spatial analysis of ecosystem service value in Shijiazhuang under multi-scenarios. Environ. Sci. Pollut. Res. 2021, 28, 31043–31058. [Google Scholar] [CrossRef] [PubMed]
  50. Brunsdon, C.; Fotheringham, A.S.; Charlton, M.E. Geographically weighted regression: A method for exploring spatial nonstationarity. Geogr. Anal. 1996, 28, 281–298. [Google Scholar] [CrossRef]
  51. Noresah, M.S.; Ruslan, R. Modelling urban spatil structure using Geographically weighed Ression. In Proceedings of the 18th World IMACS Congress and MODMS09 International Congress on Modelling and Simulation, Cairns, Australia, 13–17 July 2009; Australian National University: Canberra, Australia, 2009; Volume 9, pp. 13–17. [Google Scholar]
  52. Platt, R.V. Global and local analysis of fragmentation in a mountain region of Colorado. Agric. Ecosyt. Environ. 2004, 2, 207–218. [Google Scholar] [CrossRef]
  53. Chen, X. Comparison and analysis of spatial prediction and variability of soil pH in Anhui Province based on three kinds of geographically weighted regression. Sci. Geogr. Sin. 2023, 5, 5–6. [Google Scholar]
  54. Hasan, S.S.; Zhen, L.; Miah, M.G.; Ahamed, T.; Samie, A. Impact of land use change on ecosystem services: A review. Environ. Dev. 2020, 34, 100527. [Google Scholar] [CrossRef]
  55. Fu, B.; Zhang, L. Land-use change and ecosystem services: Concepts, methods and progress. Prog. Geogr. 2014, 33, 441–446. [Google Scholar]
  56. Talukdar, S.; Singha, P.; Shahfahad; Mahato, S.; Praveen, B.; Rahman, A. Dynamics of ecosystem services (ESs) in response to land use land cover (LU/LC) changes in the lower Gangetic plain of India. Ecol. Indic. 2020, 112, 106121. [Google Scholar] [CrossRef]
  57. Wang, Y.; Shen, J.; Xiang, W. Ecosystem service of green infrastructure for adaptation to urban growth: Function and configuration. Ecosyst. Health Sustain. 2018, 4, 132–143. [Google Scholar] [CrossRef] [Green Version]
  58. Zheng, W.; Ke, X.; Xiao, B.; Zhou, T. Optimising land use allocation to balance ecosystem services and economic benefits-A case study in Wuhan, China. J. Environ. Manag. 2019, 248, 109306. [Google Scholar] [CrossRef]
  59. Jin, G.; Deng, X.; Chu, X.; Li, Z.; Wang, Y. Optimization of land-use management for ecosystem service improvement: A review. Phys. Chem. Earth Parts A/B/C 2017, 101, 70–77. [Google Scholar] [CrossRef]
  60. Jin, G.; Wang, P.; Zhao, T.; Bai, Y.; Zhao, C.; Chen, D. Reviews on land use change induced effects on regional hydrological ecosystem services for integrated water resources management. Phys. Chem. Earth Parts A/B/C 2015, 89–90, 33–39. [Google Scholar] [CrossRef]
  61. Liu, J.; Li, J.; Gao, Z.; Yang, M.; Qin, K.; Yang, X. Ecosystem Services Insights into Water Resources Management in China: A Case of Xi’an City. Int. J. Environ. Res. Public Health 2016, 13, 1169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Picchi, P.; van Lierop, M.; Geneletti, D.; Stremke, S. Advancing the relationship between renewable energy and ecosystem services for landscape planning and design: A literature review. Ecosyst. Serv. 2019, 35, 241–259. [Google Scholar] [CrossRef]
  63. Xie, B.; Meng, S.; Zhang, M. Evolution of Habitat Quality and Its Response to Topographic Gradient Effect in a Karst Plateau: A Case Study of the Key Biodiversity Conservation Project Area of Wuling Mountains. Int. J. Environ. Res. Public Health 2023, 20, 331. [Google Scholar] [CrossRef]
Figure 1. Map of the Guangdong–Hong Kong–Macao Greater Bay Area, China.
Figure 1. Map of the Guangdong–Hong Kong–Macao Greater Bay Area, China.
Sustainability 15 09691 g001
Figure 2. Article flow chart.
Figure 2. Article flow chart.
Sustainability 15 09691 g002
Figure 3. Schematic diagram of grid scale in the research area.
Figure 3. Schematic diagram of grid scale in the research area.
Sustainability 15 09691 g003
Figure 4. Land use interpretation map of Guangdong, Hong Kong and Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Figure 4. Land use interpretation map of Guangdong, Hong Kong and Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Sustainability 15 09691 g004
Figure 5. Individual ecosystem service values in the Guangdong–Hong Kong–Macao Greater Bay Area.
Figure 5. Individual ecosystem service values in the Guangdong–Hong Kong–Macao Greater Bay Area.
Sustainability 15 09691 g005
Figure 6. Ecosystem service values class for the Guangdong–Hong Kong–Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Figure 6. Ecosystem service values class for the Guangdong–Hong Kong–Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Sustainability 15 09691 g006
Figure 7. Areas of change in ecosystem service values class in 2000 (a), 2010 (b) and 2020 (c).
Figure 7. Areas of change in ecosystem service values class in 2000 (a), 2010 (b) and 2020 (c).
Sustainability 15 09691 g007
Figure 8. Local spatial autocorrelation of ecosystem service values in the Guangdong–Hong Kong–Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Figure 8. Local spatial autocorrelation of ecosystem service values in the Guangdong–Hong Kong–Macao Greater Bay Area in 2000 (a), 2010 (b) and 2020 (c).
Sustainability 15 09691 g008
Figure 9. Elevation local R2 and predicted map from 2000 to 2020 (a) and slope local R2 and predicted map from 2000 to 2020 (b).
Figure 9. Elevation local R2 and predicted map from 2000 to 2020 (a) and slope local R2 and predicted map from 2000 to 2020 (b).
Sustainability 15 09691 g009
Table 1. Ecosystem service values (1 × 108 CNY) per unit area in the Guangdong–Hong Kong–Macao Greater Bay Area.
Table 1. Ecosystem service values (1 × 108 CNY) per unit area in the Guangdong–Hong Kong–Macao Greater Bay Area.
Ecosystem ServicesArable LandWoodlandGrasslandWaterUnused LandConstruction Land
Food production1833.63605.10788.46971.8236.67-
Raw material production715.125464.22660.11641.7773.35-
Atmospheric regulation1320.217921.282750.44935.15110.02-
Climate regulation1778.627462.872860.463777.28238.37-
Hydrological regulation1411.907499.552805.4534,417.24128.35−6678.00
Waste disposal2548.753153.842420.3927,229.41476.74−2174.10
Soil maintenance2695.447371.194107.33751.79311.723480.00
Biodiversity maintenance1870.308269.673428.896289.35733.45-
Recreation311.723813.951595.268141.32440.07-
Total14,485.6851,561.6821,416.8083,155.122548.75−5372.10
Table 2. Land Use Transfer Matrix for the Greater Bay Area of Guangdong, Hong Kong and Macau (km2).
Table 2. Land Use Transfer Matrix for the Greater Bay Area of Guangdong, Hong Kong and Macau (km2).
20002010
Arable LandWoodlandGrasslandWaterConstruction LandUnused Land
Arable land11,794.79200.476.57640.161788.170.14
Woodland146.5429,617.0523.7465.25755.930.39
Grassland10.1476.231063.2311.5264.400.03
Water542.5639.284.263152.38511.800.03
Construction land130.0392.713.7971.544135.000.05
Unused land2.100.910.023.236.6110.63
20102020
Arable LandWoodlandGrasslandWaterConstruction LandUnused Land
Arable land10,971.15760.7336.61244.382407.709.73
Woodland617.7728,671.65221.40115.43973.958.69
Grassland42.13270.49820.7311.5979.910.71
Water1169.43144.2511.692108.89736.7179.34
Construction land232.96122.175.6334.434026.3011.63
Unused land3.441.850.140.186.7411.17
Table 3. Land use focus shift in the Guangdong, Hong Kong and Macao Greater Bay Area (m).
Table 3. Land use focus shift in the Guangdong, Hong Kong and Macao Greater Bay Area (m).
Type of Land Use2000–20102010–2020
Transfer DistanceTransfer DirectionTransfer DistanceTransfer Direction
Arable land3733Northwest3625West
Woodland1171Northwest591Southeast
Grassland2935Northwest19,962Northwest
Water3448South4276Southeast
Construction land3602North936Southeast
Unused land11,086East24,404Southwest
Table 4. Changes of ecosystem service values (1 × 108 CNY) in the Guangdong–Hong Kong–Macao Greater Bay Area.
Table 4. Changes of ecosystem service values (1 × 108 CNY) in the Guangdong–Hong Kong–Macao Greater Bay Area.
Ecosystem Service Functions200020102020
ESVContribution Rate (%)ClassESVContribution Rate (%)ClassESVContribution Rate (%)Class
Supply ServicesFood production50.552.33946.692.25945.292.369
Raw material production182.238.417177.958.567175.069.127
Reconciliation ServicesAtmospheric regulation270.5012.485263.5912.685259.4413.524
Climate regulation275.7512.724267.6012.874258.9713.505
Hydrological regulation377.0517.401344.9016.591277.1614.443
Support ServicesWaste disposal247.6811.436229.2811.036180.659.416
Soil maintenance289.9213.383289.8913.943290.5815.142
Biodiversity maintenance313.7114.482304.3414.642291.2815.181
Cultural ServicesRecreation159.837.378155.277.478140.347.318
Total2167.22100-2079.50100-1918.76100-
Table 5. Area statistics for different classes of ecosystem service values (number of cells).
Table 5. Area statistics for different classes of ecosystem service values (number of cells).
YearVery Low Low MediumHighVery High
200010921577109111571534
20101543115993310551761
2020175799291610461740
Table 6. Area statistics on changes in ecosystem service value class (number of cells).
Table 6. Area statistics on changes in ecosystem service value class (number of cells).
YearESV Reduction Zone (<−2.4 × 106 CNY/hm2)ESV No Significant Change Zone (−2.4–2.4 × 106 CNY/hm2)ESV Increase Zone (>2.4 × 106 CNY/hm2)
2000–2010187523872189
Proportion293734
2010–202020983373980
Proportion335215
2000–2020228624361729
Proportion353827
Table 7. Projections of ecosystem service values (1 × 108 CNY) in the Guangdong–Hong Kong–Macao Greater Bay Area in 2030.
Table 7. Projections of ecosystem service values (1 × 108 CNY) in the Guangdong–Hong Kong–Macao Greater Bay Area in 2030.
Ecosystem Service TypesEcosystem ServicesArable LandWoodlandGrasslandWaterUnused LandConstruction LandTotal
Supply ServicesFood production23.93018.1210.8611.6970.0060.00044.616
Raw material production9.333163.6410.7211.1210.0120.000174.828
Reconciliation ServicesAtmospheric regulation17.230237.2243.0051.6330.0180.000259.110
Climate regulation23.212223.4963.1256.5980.0390.000256.470
Hydrological regulation18.426224.5943.06560.1150.021−60.090246.131
Support ServicesWaste disposal33.26394.4502.64447.5600.078−19.563158.433
Soil maintenance35.177220.7504.4871.3130.05131.314293.093
Biodiversity maintenance24.409247.6583.74610.9850.1210.000286.918
Cultural ServicesRecreation4.068114.2191.74314.2200.0720.000134.322
Total189.0481544.15423.395145.2440.419−48.3391853.921
Table 8. Topographic Factor Regression Parameter Table for Guangdong, Hong Kong and Macau Greater Bay Area.
Table 8. Topographic Factor Regression Parameter Table for Guangdong, Hong Kong and Macau Greater Bay Area.
Year200020102020
Local R20–0.85410–0.85320–0.8584
R20.80220.82670.8317
Adjusted R20.78830.81460.8200
AICC256,353.8256,973.8256,797.1
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Li, H.; Huang, Y.; Zhou, Y.; Wang, S.; Guo, W.; Liu, Y.; Wang, J.; Xu, Q.; Zhou, X.; Yi, K.; et al. Spatial and Temporal Evolution of Ecosystem Service Values and Topography-Driven Effects Based on Land Use Change: A Case Study of the Guangdong–Hong Kong–Macao Greater Bay Area. Sustainability 2023, 15, 9691. https://doi.org/10.3390/su15129691

AMA Style

Li H, Huang Y, Zhou Y, Wang S, Guo W, Liu Y, Wang J, Xu Q, Zhou X, Yi K, et al. Spatial and Temporal Evolution of Ecosystem Service Values and Topography-Driven Effects Based on Land Use Change: A Case Study of the Guangdong–Hong Kong–Macao Greater Bay Area. Sustainability. 2023; 15(12):9691. https://doi.org/10.3390/su15129691

Chicago/Turabian Style

Li, Hui, Yilin Huang, Yilu Zhou, Shuntao Wang, Wanqi Guo, Yan Liu, Junzhi Wang, Qing Xu, Xiaokang Zhou, Kexin Yi, and et al. 2023. "Spatial and Temporal Evolution of Ecosystem Service Values and Topography-Driven Effects Based on Land Use Change: A Case Study of the Guangdong–Hong Kong–Macao Greater Bay Area" Sustainability 15, no. 12: 9691. https://doi.org/10.3390/su15129691

APA Style

Li, H., Huang, Y., Zhou, Y., Wang, S., Guo, W., Liu, Y., Wang, J., Xu, Q., Zhou, X., Yi, K., Hou, Q., Liao, L., & Lin, W. (2023). Spatial and Temporal Evolution of Ecosystem Service Values and Topography-Driven Effects Based on Land Use Change: A Case Study of the Guangdong–Hong Kong–Macao Greater Bay Area. Sustainability, 15(12), 9691. https://doi.org/10.3390/su15129691

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop