The Coupling and Coordinated Development from Urban Land Using Benefits and Urbanization Level: Case Study from Fujian Province (China)

In recent years, urbanization has been developing rapidly. However, it is also accompanied by land management problems, such as low land use efficiency. In this research, we manage to explore the temporal and spatial evolution laws as well as characteristics of the coupling and coordinated development between urbanization and land use benefits. Through this, it is possible for us to provide policy recommendations for the sustainable development of the urbanization in Fujian Province. In this study, we take prefecture-level municipal districts and county-level cities in Fujian as the research subject. We construct an index system, based on data in 2002, 2005, 2010, 2015, and 2017, to evaluate the urban land use benefits and urbanization. Besides, we leverage the Gini coefficient weighting method to give weight to each index and calculate the value of its benefits. Moreover, it is the relative development degree and the coupling coordination degree model that we comprehensively leverage to study the spatiotemporal evolution law of the coupling coordination degree (CCD). The results show that: (1) Urban land use benefits and urbanization level are positively correlated with the regional administrative level and economic development status; (2) The CCD of urban land use benefits and urbanization level in various regions of Fujian is still low. However, the overall development direction is good; (3) From the perspective of spatial distribution, the CCD owns a “center-periphery” pattern that is based on the law of diminishing CCD power from three central cities of Fuzhou, Xiamen, and Sanming. Consequently, it requires governments to take action. Firstly, they should promote the intensive land use in the urbanization process. Meanwhile, they should also pay attention to ecological environment protection. Besides, it is recommendable to give full play to the radiating and leading effect of central cities on surrounding ones. Finally, they are required to provide appropriate policies and resource support to peripheral cities.


Introduction
As a basic carrier of human activities and urban operations, land resources are a basic guarantee of urban socioeconomic stability. Their effects on urbanization have been gradually strengthening. The land resources are not only closely related to the development of urbanization itself, but also to the food security and social stability [1][2][3]. In recent years, the rapid progress of China's urbanization has brought huge economic benefits. Nevertheless, it has developed at the expense of the ecological 3 of 18 regions along the southeast coast. Fujian is the core area of the 21st Century Maritime Silk Road supported by the Party Central Committee [50]. Excellent geographical location and policy support in Fujian promote its development greatly. Nevertheless, its unique geographical feature of "eight mountains, one river, and one field" demonstrates the lack of available land resources. This defect severely limits the socioeconomic development and urbanization process of Fujian. Consequently, the intensive utilization of land, which contributes to the effective coordinated development of land use benefits and urbanization, is particularly significant for promoting the socioeconomic advancement of Fujian.
In view of the above theoretical and practical background, this study, taking prefecture-level municipal districts and county-level cities in Fujian Province as the research object, builds an index system. This system is based on statistical data in 2002, 2005, 2010, 2015, and 2017, to evaluate the urban land use benefits and urbanization level. We also comprehensively leverage the Gini coefficient weighting method and the CCD model to study the spatial and temporal evolution of the coupling coordination degree (CCD). The contributions of this paper are as follows. Firstly, the utilization of the new weight determination method, the Gini coefficient weighting method (GCWM), can reflect the information more fully. Secondly, this paper, taking county-level cities and prefecture-level municipal districts as the research object, reflects the local situation more accurately.

Study Area
Fujian Province is situated at 23 • 33 N-28 • 20 N and 115 • 50 E-120 • 40 E southeastern China. It is on the coast of the East China Sea, across the Taiwan Strait, and opposite to Taiwan. It is an important estuary in mainland China (see Figure 1), adjacent to Zhejiang Province in the northeast, crossing the Wuyi Mountains and bordering with Jiangxi Province in the northwest, connecting with Guangdong Province in the southwest, as well as connecting the Yangtze River Delta and the Pearl River Delta. As of the end of 2017, Fujian Province has jurisdiction over 9 prefecture-level cities of Fuzhou, Xiamen, Quanzhou, Zhangzhou, Putian, Longyan, Sanming, Nanping, and Ningde (for convenience, we all call the districts of prefecture-level cities below according to the names of prefecture-level cities. For example, we call Fuzhou city as Fuzhou). There are a total of 29 districts. Fujian spans about 124,000 square kilometers of land, with 14 county-level cities and 44 counties. It has a permanent population of 39.11 million [51]. The territory of the regional city mainly consists of mountains and hills, accounting for more than 80% of the total area. River valleys and basins are covered with dense forests [52]. The sea area of Fujian, reaching 136,000 square kilometers, is slightly larger than the land. The coastline is long and winding, with unique harbor resources. The weather in Fujian is mostly the subtropical marine monsoon climate, which is warm and humid. Fujian's forest overage rate has ranked first in the country for 40 years consecutively. It reached as high as 65.95% in 2017, thanks to its superior geographic location and natural conditions [51].
In recent years, the economy in Fujian has developed rapidly. As shown in Table 1, Gross Domestic Product (GDP) increased from 0.45 trillion yuan in 2002 to 3.22 trillion yuan in 2017, with an annual growth rate of 14% [51,53]. Economic development has promoted the improvement of people's living standards. From 2002 to 2017, the per capita disposable income of urban residents increased from 9189 yuan to 39,001 yuan, multiplied by nearly three times [51,53]. The per capita disposable income of rural residents increased from 3539 yuan to 16,335 yuan, multiplied by 3.6 times [51,53]. With economic development, the urbanization process is also steadily advancing. In 2002, the urban developed area was about 502.58 square kilometers, while increasing to 1516.88 square kilometers by 2017 [51,53]. The urban area has been swallowing huge amounts of land for agricultural and environmental uses. The scarcity of available land resources has gradually limited the development of urbanization. With the implementation of the "One Belt One Road" initiative in China and the deepening of economic globalization, Fujian has gained an increasingly prominent status. Consequently, it is of great significance to strengthen the research on the relationship between the urbanization in Fujian and the land use benefits. It is also necessary to explore the temporal and spatial evolution of its coordinated development proposing feasible plans. In this way, it is possible for us to promote the coordinated development and comprehensive management of cities in Fujian. Besides, we are also likely to avoid human-land conflicts and problems brought by urbanization. Finally, we can create a city cluster that is reasonably and scientifically organized. the urbanization in Fujian and the land use benefits. It is also necessary to explore the temporal and spatial evolution of its coordinated development proposing feasible plans. In this way, it is possible for us to promote the coordinated development and comprehensive management of cities in Fujian. Besides, we are also likely to avoid human-land conflicts and problems brought by urbanization. Finally, we can create a city cluster that is reasonably and scientifically organized.

Data Source and Processing
In view of data use, previous studies often used prefecture-level municipal districts data. However, it was difficult to accurately reflect the relationship between local land use benefits and the urbanization level. Moreover, compared with data of county-level cities and prefecture-level municipal districts, those of counties are incomplete. Consequently, we finally leverage the data of county-level cities and prefecture-level municipal districts in Fujian for this research. In terms of time, due to the missing data in some years, we

Data Source and Processing
In view of data use, previous studies often used prefecture-level municipal districts data. However, it was difficult to accurately reflect the relationship between local land use benefits and the urbanization level. Moreover, compared with data of county-level cities and prefecture-level municipal districts, those of counties are incomplete. Consequently, we finally leverage the data of county-level cities and prefecture-level municipal districts in Fujian for this research. In terms of time, due to the missing data in some years To mitigate the impact of different dimensions or magnitude orders of indices, we first perform dimensionless processing to the data. We divide the indices into two categories, positive index and negative index according to their effect on the system [54]. Then, we standardize the data by Equation (1).
The X stk is the value and x stk is the standardized value. s is the number of regions, ranging from 1 to 23 (

Building the Index System of Urban Land Use Benefits and Urbanization Level
Constructing a scientific index evaluation system is a prerequisite for evaluating the coupling degree of urban land use benefits and urbanization [55]. To ensure scientificity, integrity, hierarchy, and operability, this research, based on the idea of synergy in physics and previous research, constructs a comprehensive evaluation index system of urban land use benefits and urbanization in Fujian. Urban land benefits are the four in one of the economic, social, ecological, and environmental benefits [14]. To comprehensively reflect this connotation, instead of only focusing on economic benefits, we evaluate urban land use benefits from the four dimensions of economic, social, ecological, and environmental benefits. Urbanization has a multidimensional meaning. It mainly includes the four interacting aspects of population migration, economic development, spatial expansion, and improvement of living standards [56]. Therefore, we perform a comprehensive evaluation of the four dimensions of population, economic, social, and spatial urbanization. In this way, we can avoid the one-dimensional view of the spatial transfer of the rural population that the government often thinks.
Specifically, to the index layer, the economic benefits of urban land use refer to the value of products and services produced per unit area of land. They can be measured by GDP, industrial output value, and fixed asset input [40]. Social benefits are mainly measured by population density, urban road area per capita, and developed area per capita [3]. For ecological benefits, the per capita park area, as well as the green coverage and area rate of the developed area, can fully reflect the urban ecological status. For environmental benefits, we usually use the pollutant compliance and removal rate for reflection. Due to the defects of statistical work in the statistical department, a large number of indices, such as industrial waste-water compliance rate and industrial solid comprehensive utilization rate, are missing. Consequently, we can only leverage the sewage and harmless treatment rate of domestic garbage to measure environmental benefits [57]. In the urbanization system, population urbanization refers to the process of transforming the agricultural population into the urban ones. Therefore, we leverage the rate of population urbanization and the number of nonagricultural populations to reflect the process of population urbanization [40]. In the process of urbanization, the economy is accumulating in the secondary industry and the tertiary one. Besides, the proportion of the tertiary industry keeps expanding. Therefore, we select per capita GDP, per capita industrial production value, and the proportion of tertiary industry to reflect economic urbanization [58]. In addition, social urbanization is a process that public service facilities continue to improve. Consequently, we choose indices, such as the number of hospital beds per 10,000 people, buses per 10,000 people, ordinary teachers per 10,000 people, and the total wages of urban employees, to reflect the process of social urbanization. We also reflect the process from multiple angles such as medical treatment, transportation, education, and wage level. The spatial urbanization is represented by the expansion of urban construction land, which is measured by urban construction land area and proportion of construction land [59].
Finally, we determined 8 primary indices and 22 secondary indices, as shown in Table 2.

Index Weighting
The weighting of indices is the core of multiattribute decision-making. The weighting method mainly includes subjective and objective types. The subjective weighting method mainly weights by expert experience. However, this method is likely to be interfered with by subjective factors. Besides, it also does not fully utilize the information in the data. Comparatively, the core idea of the objective one is to give indices weight by comparing the content or differentiation of the data in indices. The more obvious the change in a certain index is, the richer the information the index contains, and the heavier the weight is [60]. Overall, the objective weighting method is more appropriate. Nevertheless, it is also a problem to choose the most suitable one from various objective weighting methods.
The Gini coefficient is an important analysis index in economics. It is internationally leveraged to comprehensively measure the difference in the distribution of income among residents. Its method of calculation coincides with the core idea of the objective weighting method. For both, the greater the degree of data differentiation is, the larger the value is. The Gini coefficient weighting method draws on this method, weighting each index by deriving the Gini coefficient value. It comprehensively reflects the difference between any two data in the same index. It fully utilizes the value of information. It scientifically and objectively reflects the differentiation (distinction) of the data in a certain index. In addition, the definition of the Gini coefficient itself eliminates the effect of dimensions. Consequently, there is no need for us to standardize data in advance. As a result, it avoids the loss of information in data processing. It simplifies the calculation. Besides, its applicability is stronger and order preservation is better, compared to other objective weighting methods [61]. Therefore, we leverage the Gini coefficient weighting method in our index system. The steps are as follows [61]: Calculate the Gini coefficient value of the evaluation indices, as shown in Equations (2) and (3): where G k is the Gini coefficient value of the k-th index. m is the number of evaluation indices. n is the number of samples of the indices. X ki refers to the i-th sample of the k-th index. µ k refers to the sample of the k-th index. In particular, when the average value of the index data a not 0, the Gini coefficient value is calculated by Equation (2). When it is 0, the value is calculated by Equation (3).

2.
Calculate the Gini coefficient weight of the evaluation index: After calculating the Gini coefficient values G k of m indices in the evaluation system by Equation (2) or (3), and then normalizing them, the Gini coefficient weight g k of the k-th index can be derived, as shown in Equation (4): where g k is the weight of the Gini coefficient of the k-th index. G k is the value of the Gini coefficient of the k-th index. m is the number of evaluation indices. The sum of all weights is 1. The specific weight values are shown in Table 2.

Relative Development Model
The relative development degree is a specific index that reflects the urban land use benefits as well as the relative development degree and status of urbanization in a certain period. The calculation steps are as follows [58].

1.
Calculate the overall benefits of urban land use benefits and urbanization system, as shown in Equations (5) and (6): a k x stk (s = 1, 2, . . . , n 1 and t = 1, 2, . . . , n 2 ), b k y stk (s = 1, 2, . . . , n 1 and t = 1, 2, . . . , n 2 ), Among them, f st (x) and g st (y) are the overall benefits of the urban land use benefits and urbanization system in the t-th year in the s area. n 1 and n 2 refer to the sample sizes in areas and years, respectively. x stk and y stk refer to the standardized values of the k-th index in the t-th year in the s-th area, respectively. a k and b k refer to the weights of the corresponding indices, respectively. m 1 and m 2 refer to the numbers of indices in the urban land use benefits and urbanization system, respectively. The two are added as the indices in the evaluation system, that is, m 1 + m 2 = m.

2.
Calculate the relative development degree, as shown in Equation (7): . . , n 1 and t = 1, 2, . . . , n 2 ), where E st is the relative development degree of urban land use benefits and urbanization in the s-th area in the t-th year. When E st > 1, we call it an advanced city; when E st = 1, we call it a synchronous city; when E st < 1, we call it a lagging city.

Coupling Coordination Degree (CCD) Model
Refer to Tang's studies, we build the CCD model by the following steps [62], as shown in Equations (8)- (10): . . , n 1 and t = 1, 2, . . . , n 2 ), . . , n 1 and t = 1, 2, . . . , n 2 ), D st (x, y) = C st × T st (s = 1, 2, . . . , n 1 and t = 1, 2, . . . , n 2 ), Among them, C st is the coupling degree of urban land use benefits and urbanization of the s-th area in the t-th year, and C st ∈ [0, 1]. T st is the coordination degree of urban land use benefits and urbanization of the s-th area in the t-th year. D st is the coupling coordination degree of urban land use benefits and urbanization of the s-th area in the t-th year, and D st ∈ [0, 1]. α and β are the contributions of the urban land use benefits system and the urbanization system, respectively. According to the existing literature, we cannot conclude whether the urban land use benefits and the urbanization are more important, so we take 0.5 for both α and β. According to the existing research, we divide the CCD into three stages and 10 categories [63], as shown in Table 3.

Weight Analysis
It can be seen from Table 2 that the difference between the weights of the secondary indices is huge. The largest is 0.1084, while the smallest is 0.0062, 17.48 times as large as the latter. Among the 22 secondary indices, there are three indices, concentrated on economic benefits and spatial urbanization, which are above 0.1. There are five indices between 0.05 and 0.099. Fourteen are below 0.05. As shown in Table 4, among the primary indices, the weight of the economic benefit index is the largest, close to 1/3 of the total. Besides, the difference between the weight of the ecological benefits index and the environmental benefits index is not large, far lower than the weights of other primary indices.
To conclude, it shows that the gap between ecological and environmental benefits among regions is narrow. However, the gap between economic benefits and spatial urbanization is huge, comparatively.

Relative Development Analysis
From Figure   The relative development degree is derived by the ratio of the urban land use benefits and the urbanization level. Considering this, we cannot help asking whether urban land use benefits are lower in areas with a low relative development degree. Taking 2017 as an example, we make a comparison chart of urban land use benefits and urbanization level in various regions, as shown in Figure 3. The relative development degree is derived by the ratio of the urban land use benefits and the urbanization level. Considering this, we cannot help asking whether urban land use benefits are lower in areas with a low relative development degree. Taking 2017 as an example, we make a comparison chart of urban land use benefits and urbanization level in various regions, as shown in Figure 3. The relative development degree is derived by the ratio of the urban land use benefits and the urbanization level. Considering this, we cannot help asking whether urban land use benefits are lower in areas with a low relative development degree. Taking 2017 as an example, we make a comparison chart of urban land use benefits and urbanization level in various regions, as shown in Figure 3. From Figure 3, the urban land use benefits in areas with relatively low development are higher. Correspondingly, the urbanization level is also higher. From the perspective of the administrative level, the urban land use benefits and urbanization level of prefecture-level municipal districts, such as Fuzhou, Xiamen, Zhangzhou, Quanzhou, Sanming, and Longyan, are generally higher than their county-level cities. From the perspective of the regional economic development level, Fuzhou, Xiamen, Quanzhou, Zhangzhou, and other well-developed regions have higher levels of From Figure 3, the urban land use benefits in areas with relatively low development are higher. Correspondingly, the urbanization level is also higher. From the perspective of the administrative level, the urban land use benefits and urbanization level of prefecture-level municipal districts, such as Fuzhou, Xiamen, Zhangzhou, Quanzhou, Sanming, and Longyan, are generally higher than their county-level cities. From the perspective of the regional economic development level, Fuzhou, Xiamen, Quanzhou, Zhangzhou, and other well-developed regions have higher levels of urbanization. However, the discrepancy of land use benefits between well-developed regions and others is much narrower than that of urbanization levels.

Analysis of the Coupling Coordination Degree (CCD) of Urban Land Use Benefits and Urbanization Level
From Table 5, we can see that in the five years of 2002, 2005, 2010, 2015, and 2017, regions whose coupling degree C of urban land use benefits and urbanization level were above 0.9 accounted for 100%, 100%, 95.6%, 95%, and 100%, respectively. Only ratios of Xiamen in 2010 and Shaowu in 2015 were below 0.9, while still between 0.8 and 0.9. This variable is relatively stable. The coupling coordination degree D was between 0.2 and 0.6. The change is relatively obvious. From Figure 4, it can be seen that the coupling coordination degree is mainly in four stages. They are moderately uncoordinated, slightly uncoordinated, at the edge of being uncoordinated, and barely coordinated. It is moderately uncoordinated and slightly uncoordinated that dominate in them. These two stages accounted for 74%, 74%, 82.6%, 68.2%, and 66.7% of the previous years, respectively. No region has reached the stage of coordinated development. However, most have moved from moderately uncoordinated to slightly uncoordinated, with an overall trend to develop well. Considering the stage of each region, there are obvious trends and laws of distribution.
In view of the time series, from 2010 to 2015, there were eight cities that achieved the improvement of the coupling coordination phase. The number far exceeded the total of other adjacent years. In addition, only the coupling coordination phase in one city had a reduction. This number of reductions was also far lower than that in other adjacent years. From the perspective of spatial distribution, regions with relatively high coupling coordination, such as Xiamen, Fuzhou, Quanzhou, and Shishi, are all coastal regions. Except for Sanming, it is an inland city. Moreover, Xiamen and Fuzhou become the coastal central cities while Sanming is the inland central city. Consequently, there grows a sort of "central-peripheral" development pattern. That is, for the surrounding cities of these three central cities, the closer to the central city they are, the higher their coordination degree is, and vice versa.   In view of the time series, from 2010 to 2015, there were eight cities that achieved the improvement of the coupling coordination phase. The number far exceeded the total of other adjacent years. In addition, only the coupling coordination phase in one city had a reduction. This number of reductions was also far lower than that in other adjacent years. From the perspective of spatial distribution, regions with relatively high coupling coordination, such as Xiamen, Fuzhou, Quanzhou, and Shishi, are all coastal regions. Except for Sanming, it is an inland city. Moreover, Xiamen and Fuzhou become the coastal central cities while Sanming is the inland central city. Consequently, there grows a sort of "central-peripheral" development pattern. That is, for the surrounding cities of these three central cities, the closer to the central city they are, the higher their coordination degree is, and vice versa.

Discussion
The type of mainstream cities in Fujian Province has changed from leading cities to lagging ones. It reflects that in the urbanization process, the development of urban land use benefits lags behind urbanization. This also confirms that the rapid development of urbanization has brought many defects, such as the spreading development of urbanization and inefficient land use [64,65]. The main reason is that, in the process of urbanization, local governments blindly expand urban areas. They overpursue urbanization while neglecting the efficient utilization of land. The "BIG PIE" policy has resulted in low urban land use benefits, far behind the development of urbanization. Moreover, the deeper reason is related to the assessment mechanism of government officials. The assessment mainly considers economic indices. Since urban land expansion has a motivating effect on economic development [66,67], government officials are encouraged to expand urban areas. It causes that officials only focus on socioeconomic benefits in land planning, while ignoring the ecological and

Discussion
The type of mainstream cities in Fujian Province has changed from leading cities to lagging ones. It reflects that in the urbanization process, the development of urban land use benefits lags behind urbanization. This also confirms that the rapid development of urbanization has brought many defects, such as the spreading development of urbanization and inefficient land use [64,65]. The main reason is that, in the process of urbanization, local governments blindly expand urban areas. They overpursue urbanization while neglecting the efficient utilization of land. The "BIG PIE" policy has resulted in low urban land use benefits, far behind the development of urbanization. Moreover, the deeper reason is related to the assessment mechanism of government officials. The assessment mainly considers economic indices. Since urban land expansion has a motivating effect on economic development [66,67], government officials are encouraged to expand urban areas. It causes that officials only focus on socioeconomic benefits in land planning, while ignoring the ecological and environmental benefits. As a result, it finally leads to the fact that urban land use benefits are lower than the urbanization level.
Besides, we find that the urban land use benefits and urbanization level are related to the administrative level and economic development of the city. It is manifested that areas with higher administrative and economic development levels tend to own higher urban land use efficiency and urbanization level. However, the impact of economic conditions on urbanization is greater than that of the urban land use benefits. The reason is not complex to explain. One the one hand, areas with high administrative level will receive more policy support and resource tilt [68]. As a result, it is more beneficial to urban construction and economic development. Comparatively, economically well-developed areas tend to pay more attention to economic and social benefits. Nevertheless, economic development often comes at the price of the environment. This partially hinders the improvement of overall benefits. Consequently, the difference in urban land use benefits among cities is insignificant.
Moreover, it is also found that the CCD is still at a relatively low level in various regions of Fujian. However, it develops in a good direction. There are situations of high coupling and low coupling coordination. This is mainly because the coupling degree only indicates the strength of the effect between urban land use benefits and the urbanization level. However, we do not know whether they promote or inhibit each other. Comparatively, the coupling coordination degree fully considers the coordination degree between them [69]. Early urban land use benefits limited to urbanization. However, the two gradually promote each other as time goes by. This situation does not only occur in Fujian exclusively. Jia et al. have studied the three major urban agglomerations, including the Yangtze River Delta, the Pearl River Delta, and the Beijing-Tianjin-Hebei area, which are relatively better-developed [70]. They obtained the same result. Similarly, Zuo et al. also got the same result when studying the Shaanxi-Gansu-Ningxia region, which is of lower development level [40]. It requires further study of whether this situation is suitable for the whole country. However, the increase in the CCD took place in large numbers from 2010 to 2015. It far exceeded the sum of other years. This may partly benefit from the establishment of the Western Taiwan Straits Economic Zone in 2009 and the peaceful cross-strait relations.
Finally, we find that each city generally exhibits a development pattern according to its distance from the three central cities of Fuzhou, Xiamen, and Sanming. That is, the closer to the central city it is, the higher the coordination degree it has; vice versa. What is more, cities with developed coastal transportation have a higher coupling and coordination degree than inland cities. Zhang et al. [3] and Wang et al. [58] studied the three provinces in Northeast China and the Bohai rim area, respectively, to obtain this "center-periphery" pattern. Zhang et al. found that cities with relatively high coupling and coordination degrees are distributed along the Harbin-Dalian Transportation Economic Belt in strips [3]. This finding is so interesting. It shows that the development of regions not only depends on their own influencing factors, but also on the radiation and leading role of the central cities. The farther away the distance is, the weaker this radiation and driving effect is. Meanwhile, developed traffic can also strengthen this driving effect. This provides a way for the government to give full play to the leading role of central cities. However, there are still several questions that require us to perform further research and practice. For example, we should consider how to leverage the central city's radiation and driving role more efficiently. This is significant to promote the coordinated development of land use benefits and urbanization in other cities.
When it comes to research methods, the research methods of this study are applicable to various areas with similar issues, such as urban agglomerations, provinces, and cities. However, the completeness of the data in each region, the difference of regional conditions, and the different perspectives of the target connotation will all affect the construction of the index system. The universal system still needs to be tested in practice. Besides, the selection of the index weighting method will also have an impact on the results. Weights, computed by the objective weighting method, are based on the information of the data itself [61]. Due to the discrepancy in the data of various regions, weights will change accordingly. Consequently, weights obtained in this study cannot be directly applied to other regions.
In addition, this research also has several drawbacks, which are likely to be solved in future research. (1) This study only selects 2002, 2005, 2010, 2015, and 2017 as time nodes, causing the time series data to be insufficient. Consequently, we can only tentatively study the spatial and temporal evolution of urban land use benefits and urbanization in various regions. Therefore, future research can leverage continuous-time data to analyze its evolution in more detail. (2) Considering a large number of county statistical data is missing, we only take prefecture-level municipal districts and county-level cities as research samples. However, the county is also an important part of China's urbanization process. As a result, there are limitations in the spatial analysis of urban land use benefits and urbanization in regions of Fujian; (3) In the construction of the index system, some data are missing. For example, compared with other studies, some important indices, such as the proportion of the tertiary industry's employed population [59], urban residents per capita use area [57], and the three-waste treatment rate [54], are not included. Consequently, it requires us to add more comprehensive data in future relevant studies to improve the index system, reflecting the spatiotemporal evolution of the CCD of each city more accurately.

Conclusions
It is of practical significance to study the law of spatiotemporal evolution of the coupling and coordination relation between urban land use benefits and urbanization level in Fujian. It is likely to promote the healthy development of the urbanization process. Therefore, this study takes county-level cities and prefecture-level municipal districts in Fujian Province as the research object, leveraging the data of 2002, 2005, 2010, 2015, and 2017 to construct the evaluation index system of urban land use benefits and urbanization system. We derive the CCD of the index system and analyze the temporal and spatial evolution law. The conclusions of this study are as follows: (1) Urban land use benefits and urbanization levels are positively correlated with the regional administrative level and economic development status. (2) The CCD of urban land use benefits and urbanization levels in various regions of Fujian is still low. However, the overall development direction is good. (3) In terms of spatial distribution, the CCD has a "center-periphery" pattern. That is, the closer to the three central cities of Fuzhou, Xiamen, and Sanming it is, the higher the coordination degree it has; vice versa.
The above findings can provide certain implications for policy formulation. First of all, in the process of urbanization, Fujian Province cannot blindly pursue the speed. It is necessary to plan urban expansion rationally and strengthen the intensive use of land, avoiding the development model of "BIG PIE". It is also significant to add some indices on the ecological environment to the assessment of government officials. It can avoid only concentrating on the economy while neglecting the ecological environment. It can finally promote the overall improvement of urban land use benefits. Secondly, it is necessary to give targeted priority to the improvement of the urban land use benefits or urbanization level, based on the reality of different regions. Finally, it is recommendable to strengthen the radiation and the leading role of central cities to surrounding ones. Appropriate policy and resource support for peripheral cities should be provided, such as financial subsidies for regional construction.