The Influence of Land Use Evolution on the Visitor Economy in Wuhan from the Perspective of Ecological Service Value

: This research used transfer matrix, dynamic attitude, and a linear regression model to investigate the characteristics of land-use change and evolution of ecological service values and their impacts on Wuhan’s visitor economy. The results showed that: (1) the land -use scale in the Wuhan metropolitan area changed significantly from 1990 to 2018. The area of arable land, forest land, and grassland decreased at a faster rate, whereas that of water and construction land continued to increase; (2) there were differences in the dynamic attitudes of land-use at different stages. The dynamic attitude of construction land-use changed the most with cultivated land, water area, forest land, unused land, and grassland. From 1990 to 2005, land-use change exhibited a relatively gentle trend, whereas from 2005 to 2020, it accelerated; (3) although land-use regulation service, support service, and cultural service values positively responded to tourism economic growth, their influences were dissimilar. This study clarifies the effects of urban land-use on tourism economic development and provides a reference for its effective control.


Introduction
Land is an important natural resource and a spatial carrier of human economic and social activities, and ecology. Land-use change reflects the impact of human activities on the natural environment, causing changes in surface structure (i.e., water bodies, climate, and ecology) and affect the ecosystem service value (ESV). With the rapid development of tourism in China, the urban ESV caused by changes in land-use forms has a crucial impact on the development of the visitor economy. With rapid urbanization in China, natural landscapes, such as forests, grasslands, and water bodies, have been gradually replaced by construction [1]. The contradictions between urban construction and tourism resource protection have become increasingly noticeable. Tourism is an important component of urban economies. The relationship between land use and the visitor economy has received significant academic research attention. McMurry (1930) raised the issue of the relationship between tourism and land-use in his book "The Relationship between Recreational Activities and Land-use" [2]. Since then, the land use issues associated with tourism development have gradually become a popular topic for tourism researchers.
International scholars have analyzed the impacts of different types of tourism activities, such as leisure experiences, resorts, and sightseeing tours on land use patterns [3,4]. With the rapid development of urban tourism, Weaver (1993) developed a spatial model of land use in the urban tourism belt after studying its spatial structure. He took the urban core area as the tourism center and formed the urban, suburban, rural, and remote zones in its periphery, thus, forming a spatial land layout with complementary functions [5]. Due to population growth, increased disposable time and income, im-proved transport accessibility, and being driven by new policies and governance systems, the land used for tourism and recreation is changing. As a result, land use demand has gradually shifted to urban, rural, and coastal areas [6]. Tyrväinen et al. (2014) found that the land use development mode of tourist destinations has a significant impact on the environmental quality and tourist satisfaction with destinations [7]. A good ecological environment provides many opportunities for tourism development. Saha and Paul (2020) constructed a land use change matrix to identify land use/land cover change and found that cultivated land, vegetation cover, aquaculture, water bodies, and dunes around coastal tourist attractions are rapidly decreasing [8]. Agricultural lands, fallow areas, dunes, and vegetation have been cleared to make way for the establishment of hotels, resorts, and other tourism-related enterprises, thus, gradually transforming these areas into tourist zones. Hjalager (2020) investigated 213 land use conflicts in coastal tourism in Denmark and found that the issues causing tensions included construction and reconstruction, infrastructure transformation, and landscape transformation. These are the main contradictions between land use and tourism development [9].
In contrast, domestic scholars in China have had more of a focus on the study of urban land use patterns, with tourism and recreation functions at the core. Tourism is typically regarded as an effective way to promote economic growth in the coastal areas of developing countries. Wang (2013) took Sanya, a coastal area in the south of Hainan Island, to analyze the changes in land use patterns caused by tourism development and the huge pressure brought to the local ecological environment [10]. Du (2020) studied the relationship between land use change and red tourism in Jinggangshan, China's red tourism base, and found that the demand for tourism infrastructure and service facilities increased, which promoted the conversion of forest and agricultural land to tourism. There are four tourism land use patterns, including integrated development, reach of scenic areas, traffic reach, and traditional town reach [11]. Wu (2001) provided a quantitative description of the spatial scope of urban tourism land use and explained the formation of land-use employing the extreme rent theory. He considered that the reason behind the formation of the ring-city recreation belt is that investors and tourists reach a compromise under the influence of land rent and travel costs [12]. Subsequently, more studies on the relationship between urban land use and tourism development have been conducted. All indicate that the relationship between urban land use and urban tourism is inseparable [11]. The development process of tourism activities is essentially a process of land use, which is also one of the major reasons for changes in land use patterns [13]. Tourism development has a direct or indirect impact on land use [14]. Especially, in cities where all economic activity revolves around tourism and tourism is the main industry, the demand for tourism land can become the driving force of urban land use change [15]. Chai (2021) studied the tourism driving force of land use change and found that tourism development was connected with regional land use change [16]. Campbell et al. (2005) studied the tourist land of different tourist destinations in Anhui Province and concluded that tourism development has led to evident changes in land use patterns in tourist areas and has become the main driving force of land use change in some [17]. Su (2011) found that tourism product development strategies led by eco-tourism, landscape leisure tourism, and agricultural tourism are conducive to local ecological environment protection and promote the continuous improvement of land ESVs [18]. Cao et al. (2021) determined that tourism economic development significantly impacts land use change [19]. The development of tourism and leisure has led to the rapid expansion of tourism facilities and construction land, reduction in cultivated land areas, and concentration of urban and tourism facilities land in areas rich in landscape resources [20].
Overall, it is indisputable that tourism development is one of the driving factors for urban land use change [21]. However, is the relationship between tourism development and land use only one-way? In particular, with the rapid advancement in urbanization and industrialization, will the changes in the spatial structure of urban land use also affect the development of the visitor economy? Will changes in urban land use directly determine or influence the levels of urban tourism economic development? Taking the Wuhan metropolitan area (the largest urban grouping in central China) as the example, this study analyzed scale change, dynamic attitude, and ecological value of land use using remote sensing images from 1990 to 2020 and explored the impact of land use changes on tourism economic development to provide a theoretical basis and practical guidance for urban tourism economic development and effective land use control.
The research contributions are as follows: (1) Most of the existing literature focuses on conflicts and contradictions between tourism and land use. There is a lack of research on the impact of land use on tourism economic development. Based on the perspective of ESV, this research examines the impacts of urban land use change on tourism development, representing an important supplement to the existing research literature; (2) a majority of the current research on land use change focuses on single cities, concentrating on the characteristics of change and influential factors, and the results are relatively rich. However, there are fewer studies on urban areas with river basins, and quantitative analyses of relationships between land use change and tourism development are minimal. Therefore, this research employed an econometric model to empirically test impacts, and the findings provides useful reference for the literature and future research.

Study Area
The Wuhan metropolitan area is located at 112°30′-116°07′ E, 29°05′-31°51′ N in the middle reaches of the Yangtze River in the central and eastern parts of Hubei Province. It consists of Wuhan and its surrounding cities of Huangshi, Ezhou, Xiaogan, Huanggang, Xianning, Xiantao, Tianmen, Qianjiang ( Figure 1). The total land area of the region is 5.79 × 104 km 2 . The Wuhan metropolitan area is a national pilot zone for comprehensive support reform of resource saving and environmentally friendliness. It is a significant economic engine driving the rise of central China and is one of the key hub points for the coordinated development of eastern, central, and western China with the Yangtze River Basin Economic Belt as the axis. The Wuhan metropolitan area is rich in water resources and has prominent geographical advantages. It is an important agricultural production area and is one of the most developed regions in central China. By 2020, the permanent resident population of the Wuhan metropolitan area was 98 million, and the regional GDP was 963.576 billion Chinese Yuan, accounting for 52.83% and 60.35% of Hubei Province, respectively. The Wuhan metropolitan area only accounts for one-third of Hubei Province. However, it has more than half of the population and more than 60% of the total GDP of the province. It is the core area for economic development in Hubei Province.

Data Sources and Processing
Land Remote-Sensing Satellite Thematic Mapper image data from 1990,1995,2000,2005,2010,2015, and 2020 in the Wuhan metropolitan area were selected. The spatial resolution of the data from data elevation model was 30 m. Earth Resources Data Analysis System 9.2 was used to carry out geometric accuracy correction, data fusion, and image registration for remote sensing image data. T, 2, 3, and 4, bands, which are more sensitive to green vegetation, were selected for RGB false color synthesis and equalization processing. Remote sensing images were clipped according to the administrative edge-vector graph. Supervised classification and human-computer interactive interpretation were used to comprehensively interpret remote sensing images. Then, artificial visual interpretation was used to correct the classification results with topographic maps, Google Earth images, land-use status, and other auxiliary data to obtain the land-use vector maps for corresponding years. It was found that the Kappa coefficients of each class were all over 0.82 after the classification accuracy test, indicating that the interpretation results met the research needs. According to the land use and cover classification system of the Chinese Academy of Sciences, the land use types were divided into six categories: cultivated land, forest land, grassland, construction land, water area, and unused land. In addition, the socioeconomic data for this study were from the Statistical Yearbook of Hubei Province and the Statistical Yearbooks of various cities from 1990 to 2020.

Dynamic Attitude of Single Land Use
The dynamic attitude of single land-use refers to the change in a certain land use type in a certain study period [22]. It is used to express the rate of change and range of different land use types over a certain period [23]. Its calculation formula is as follows: where K is the dynamic attitude of a certain land use type in the research period, is the year, and and represent the amount of land use at the beginning and end of the study period, respectively.

Dynamic Attitude of Comprehensive Land Use
The comprehensive land-use dynamic attitude refers to the rate of change of the entire land use type within a certain time range in the research area [24]. It is an indicator of regional differences in the rate of change of a land use type. It reflects the comprehensive impact of socioeconomic activities in the region on land use change [25]. Its formula is as follows: where is the total area of category land use type at the beginning of monitoring, ∆ − is the total area of category land use type converted to other types of land use type in the monitoring period, and is the monitoring period. The land use dynamic attitude reflects the rate of change for a land use type in the study sample area corresponding to period .

Land Use Transfer Matrix
The land use transfer matrix describes the conversion between different land use types. It not only reflects the land use type structure at the beginning and end of the period but also shows the transfer and change of various land use types during the research period. It is used to describe the direction of regional land use change and the source and composition of various land use types at the end of the study period [26]. The transition matrix generates the transition probability matrix of land use change to predict the changing trend of regional land use in a specific scenario [27]. Its mathematical form is as follows: where is the land area, is the number of land use types, and and are land use types at the beginning and end of the study period, respectively.

Measure of Ecosystem Service Value
Referring to the calculation methods of Costanza et al. (1998) [28] and Xie et al. (2008) [29], a standard ecosystem value is defined as 1/7 of the annual natural grain economic value of farmland with 1 hm 2 national average yield [30][31][32]. The economic value of one ecosystem service equivalent factor in China is 449.1 RMB/hm 2 [29]. The formula is as follows: where is the economic value of grain provided per unit area of farmland system, is the type of food crop, is the total sown area of the grain crop, is the average price of the th grain crop, is the average yield per unit area of the grain crop, and is the total sown area of all food crops.
The food production function of cultivated land was taken as the standard (value is 1) to assign the equivalent land use value of its type. The average ESV per unit area of the ecosystem in the Wuhan metropolitan area was obtained based on the survey results of different experts (Table 1). Some scholars set the ESV of construction land as 0 [33]. However, it is considered that construction land has a certain value for cultural and entertainment functions, i.e., construction land is dependent on human land use/cover, forming a unique artificial landscape pattern type [34] that functions as a component of the ecosystem and affects human well-being [35]. Therefore, this research retained the equivalent value of cultural services in construction land. The calculation formula of ESV was as follows: where ESV is the ecosystem service value, is the area (hm 2 ) of land-use type, is the ESV coefficient of the th land-use type, ESV represents the ESV in item , represents the ESV coefficient of item of the land-use type, is the value equivalent of item of ecosystem service of a certain land-use type, and is a standard equivalent ESV [36].

Linear Regression Model
Referring to the research of Macdonald (1980) [37] and Ploberger (1989) [38], this study constructed a linear regression model of the impact of land use change on the visitor economy. To reduce the heteroscedasticity of the original data and ensure that its changing trend remained unchanged, the logarithmic method was adopted for the panel data of nine cities in the Wuhan metropolitan area from 1990 to 2020. The model was as follows: where is coefficients; is the total tourism revenue of each city in the Wuhan metropolitan area; is the total ESV of each city; represents other control variables, including population (PO ), an education level (EL ), and the average wage of each city (AW ); and are the fixed effects of individual and year, respectively; and is the standard error. Table 2 shows that the land use scale of the Wuhan metropolitan area changed significantly from 1990 to 2020. First, from the perspective of various land use changes, in the past 30 years, the area of cultivated land, forest land, and grassland in the Wuhan metropolitan area decreased; however, the areas of water and construction land continued to increase. The cultivated land, forest land, and grassland areas decreased by 2668, 159, and 44 km 2 , respectively. The decrease in farmland area was 8.81%, followed by grassland at 3.09%, and woodland at 0.91%. Construction land and water area increased by 2075 (73.21%) and 770.84 (13.89%), respectively. The unutilized land only increased by 1.15 km 2 (0.67%). The economies of various cities in the Wuhan metropolitan area developed rapidly with the acceleration of China's urbanization. Urban land expansion is obvious, with large-scale occupation of cultivated land for built land, a large amount of cultivated land has been converted into lakes, wetlands, and land use has changed significantly.

Change of Land Use Scale
Second, land use in the Wuhan metropolitan area is mainly cultivated land, forest land, and water area. The area proportions of the three types of land were 52.34%, 30.38%, and 9.59%, respectively, in 1990 and 50.84%, 30.25%, and 10.35%, respectively, up to 2005. It can be seen that the cultivated land area decreased rapidly, the woodland area changed steadily, and the water area increased to a certain extent. By 2020, the land use in the Wuhan metropolitan area is still dominated by woodland, arable land, and water bodies, accounting for 47.75%, 30.12%, and 10.93%, respectively. However, it is worth noting that the expansion of construction land is more obvious, with the proportion of construction land increasing from 4.9% in 1990 to 8.49% in 2020. Thus, the demand for land use in Wuhan metropolitan area has entered a more active and prosperous period.
Finally, from the perspective of the development stage of land use change, after 2005, the expansion rate of construction land in the Wuhan metropolitan area accelerated. The rapid increase in the urbanization rate had a severe impact on land use types. Specifically, during 1990-2005 and 2005-2020, construction land and water area increased by 460 and 439, respectively, and then increased by 1614 and 331 km 2 , respectively; however, cultivated land, forest land, and grassland decreased by 882, 84, and 22 km 2 , respectively, and then decreased by 1785, 75, and 21 km 2 , respectively. It can thus be concluded that the degree of land use change from 2005 to 2020 was greater than that from 1990 to 2005. The main reason is that since the beginning of the 21st century, China's urbanization strategy has been adjusted to "actively yet prudently promote urbanization" and "coordinated development of large, medium and small cities, and small towns." Rapid urbanization accelerated the expansion of urban construction land, sharply decreased the arable land area, and slowed down the expansion rate of water bodies.
Overall, with the advancement of urbanization and the continuous improvement of tourism development, the urban scale of the central city continued to expand (Figure 2). The scale of construction land and water bodies in the Wuhan metropolitan area increased. In contrast, the share of forest land, grassland, and cultivated land decreased, and the structure of urban area land constantly adjusted. However, the latter is still the most important form of land use in the Wuhan metropolitan area. This situation will not change greatly in a short time.

Rate of Change of Land Use Type
The dynamic attitude of land-use reflects the quantitative changes of land-use types in a certain time range and is divided into single and comprehensive dynamic attitudes [22]. The dynamic attitude of single and comprehensive land use was used to measure the rate of change of land use in the Wuhan metropolitan area. The results were as follows: First, the dynamic attitude analysis of single land use, according to formula (1), in different periods, the rates of change of land-use types in the Wuhan metropolitan area showed significant differences. Specifically, the area of forest land and grassland gradually decreased; from 1990 to 2005 and from 2005 to 2020, the average annual rates of change of forest land were −0.03% and −0.1%, respectively. The average annual rate of change of the cultivated land area was −0.19% from 1990 to 2005 and −0.4% from 2005 to 2020. However, water and construction land areas increased yearly. From 1990 to 2005, the annual growth rates of water and construction land were 1.08% and 0.53%, respectively. The annual growth rates from 2005 to 2020 were 3.27% and 0.37%, respectively. In conclusion, among the six types of land, the dynamic attitude of construction land use changed the most, while the degree of dynamic change of grassland was the least. In addition, the unused land area initially showed an increasing trend and then a decreasing trend. The unused land area increased annually from 1990 to 2005 with a dynamic variation of 2.52%, while it decreased annually from 2005 to 2020 with a dynamic variation of 1.79% (Figure 3). Second, the dynamic attitude analysis of comprehensive land use reflects the overall change in all types of land in the research area during the research period [39]. According to Formula (2), the rate of land use change in the Wuhan metropolitan area increased, indicating that socioeconomic activities had a more severe impact on land patterns. Specifically, from 1990 to 2005, the comprehensive dynamic attitude of land use in the Wuhan metropolitan area was 2.50%, and the rate of land use change was relatively gentle. From 2005 to 2020, the comprehensive dynamic attitude of land use in the Wuhan metropolitan area was 9.78%, and compared with the previous period, the rate of change of land use type showed an accelerating trend, and the impact of human activities increased. In addition, from the perspective of regional differences, the dynamic attitudes of integrated land use varied greatly among different cities. From 1990 to 2005, the dynamic attitudes of soil use of Wuhan, Ezhou, and Huanggang were larger, while smaller for Huangshi, Xianning, and Xiaogan. From 2005 to 2020, the dynamic attitude of soil use of Tianmen, Huanggang, and Wuhan were larger, while smaller for Xianning, Xiantao, and Xiaogan (Figure 4). In recent years, the land use in Wuhan and Huanggang have fluctuated greatly.

Changes in Land Use Types
ArcGIS was used to classify land use in the Wuhan metropolitan area from 1990 to 2005 and from 2006 to 2020, as well as to carry out spatial superposition analysis to obtain the land use transfer matrices during these years (Table 3). From 1990 to 2020, the land use types in the Wuhan metropolitan area underwent an evident transformation: (1) From 1990 to 2005, based on maintaining an original water area of 5549 km 2 , an area of 636 km 2 underwent transformation from cultivated land to water bodies, making it the largest input; however, only 4.86 km 2 underwent transformation from construction land to water bodies, making it the smallest input. Based on maintaining an original construction land area of 2834.82 km 2 , an area of 410, 41 and 47 km 2 underwent transformation from cultivated land, water bodies, and forest land, respectively, to construction land. The three types of land mentioned above were the main sources of construction land, and cultivated land was its largest source. A total of 1186 km 2 of cultivated land was converted to another land use (such as water body, built land, grassland) of which 636 2 was converted to water bodies, accounting for 53.85%; 410 km 2 was converted to construction land, accounting for 34.6%, of which 122 km 2 was converted to forest land, accounting for 10.35%. The total area of forest land decreased annually. A total of 143 km 2 of forest land was reclaimed as cultivated land, 47 km 2 was converted into construction land, and a small part was converted into water bodies and grassland. Cultivated land was the largest input source of forest land, covering 122 km 2 , and the scale of conversion of unused land and construction land to forest land was the smallest, 2 and 3.6 km 2 , respectively. A small portion of grassland was converted to woodland, 25 km 2 ; arable land, 14 km 2 ; water area, 11 km 2 ; and construction land, 5.7 km 2 . The conversion of land in this period mainly occurred in the following areas: conversion of forest land and cultivated land to construction land, conversion of cultivated land to the water bodies, and conversion of forest grassland and cultivated land. This shows that woodland and arable land were the main sources of construction land.
(2) From 2006 to 2020, arable land was still the major contributor to the water bodies, with an input of 1038 km 2 . The input sources of construction land were scattered. The contributions of cultivated land, forest land, water bodies, and grassland to construction land were 2005, 281, 183, and 28 km 2 , respectively, indicating that construction land was an important land use type increase during this period. The increase in cultivated land mainly originated from forest land, water bodies, and construction land, with conversions of 1258, 678, and 746 km 2 , respectively. In addition, the improvement of grassland and unused land contributed 52 and 27 km 2 , respectively, to the expansion of agricultural land. The increase in forest and grassland area was mainly from cultivated land, with 1404 km 2 of cultivated land being converted to forest land, which was equal to the amount of forest and grassland that converted to agricultural land (1258 km 2 ). The changes in land use types in the two periods of 1990-2005 and 2006-2020 were mainly characterized by the conversion of forest land, grassland, and cultivated land to construction land and water bodies. Moreover, the scale of transformation of various land use types in the later stage was larger than that in the previous stage, indicating that the demand for construction land in the later stage of urban development was larger. The degree of transformation among various types of land was more intense.

Analysis of Changes in Ecological Service Value
Formulas (5)-(7) were used to calculate ESV. The analysis of changes in the total value of ecosystem services in the Wuhan metropolitan area increased annually. The ESV increased from 525.320 billion RMB in 1990 to 530.529 billion RMB in 2020, with an increase rate of 0.93% (Table 4). Specifically, the ESV of forest land and grassland decreased by 1.21 billion RMB and 226 million RMB, respectively, accounting for 8.81% and 0.91%, respectively. The rate of increase of the ESV of water bodies was 13.89%. The ESV of unused land first increased and then decreased. From 1990 to 2005, the ESV of unused land increased by 37.75% and decreased by 26.92% from 2006 to 2020. The growth rate of the ESV of construction land ranked first among the six types, accounting for 73.21%. In particular, during 2005-2020, the growth rate reached 49%. Second, the change in the ESV structure was analyzed. From the perspective of the evolution of ecological service function value, from 1990 to 2020, the value of hydrological regulation, waste treatment, aesthetic landscape provision, and cultural and entertainment ecosystem services showed an increasing trend yearly, with the rate of change of ESV being 5.48%, 4.51%, 5.2%, and 2.11%, respectively (Table 5). Other individual ecosystem services showed a declining trend. The value of food production and soil formation protection declined rapidly, with a changing rate of −6.24% and −3.73%. From 1990 to 2020, the ESV of raw material production, gas exchange, gas regulation, biodiversity conservation, natural heritage and diversity showed little change. The rates of change were −1.93%, −2.27%, −1.52%, −0.67%, and −1.96%, respectively. Regarding ecosystem services of food production, the ESV decreased by 1.69% from 1990 to 2005 and by 4.26% from 2006 to 2020. The main reason for the decrease in the ESV of food production was that a large number of rural laborers migrated to cities with the acceleration of urbanization. A large amount of cultivated land was abandoned. Meanwhile, the cultivated land in the suburbs was gradually transformed into construction land, which decreased the food production value year on year.
From the perspective of ESV composition, the main functions in the Wuhan metropolitan area are water conservation, waste treatment, soil formation and protection, and biodiversity protection. The sum of these four functions account for more than 50% of the ecosystem services. Wuhan is a typical water town, south of the Yangtze River, with a dense water network, reservoirs, lakes, and abundant water. Therefore, the hydrological regulation is the highest ecological function, accounting for nearly 20% of the total ecological value. Finally, from the perspective of the value scale of ecological services, the values of various ecological services during 1990-2020 are regulation, cultural, expenditure, and supply services. The value of a single ecological service function is hydrological regulation, waste treatment, biodiversity protection, soil formation and protection, climate regulation, gas exchange, landscape aesthetics, raw material production, food production, natural heritage and diversity, and cultural entertainment. It is worth noting that there are many rivers and lakes in the Wuhan metropolitan area, such as the Yangtze River, Han River, and East Lake [40], which significantly impact microclimate improvement and ecosystem regulation services and also make hydrological regulation services occupy the first place in various ecological service functions.

Effect of Land Use Change on the Visitor Economy
The development of the visitor economy in the Wuhan City circle is closely related to land-use change [41]. The number of visitors and total tourism income was significantly positively correlated with construction land [42]. Different types of land have dissimilar ESVs. For example, forest land provides tourists with beautiful natural environments, clean air, diverse flora and fauna, and plays a significant role in gas exchange, climate regulation, and biodiversity; grasslands are of great significance for animal and plant protection and biodiversity in scenic areas. In particular, the development of agricultural tourism converts cultivated land into a tourism function, and this plays an important role in the visitor economy. Although the contribution of building land in gas exchange, climate regulation, and food production is small, it has the service functions of landscape aesthetics, culture and entertainment, and natural heritage [43]. Land use mainly promotes tourism economic growth through the ESVs in three aspects: regulation, support, and cultural services ( Figure 5).

Impact of Land Use Change on the Visitor Economy
Different types of land are dissimilar in the value of supply, regulation, support, and cultural services. In addition, among the above four types of service values, supply services provide food production for residents, which is less related to the development of the visitor economy. In contrast, regulation, support, and cultural services are closely related to the development of the visitor economy. Therefore, this research used regulation service, supporting service, and cultural service values as independent variables, and the tourism income of each city as the dependent variables and constructed a linear regression model to explore the impact of land use change on the visitor economy in the Wuhan metropolitan area.
According to the model setting form of Formula (8), the parameters of each variable were estimated. As shown in Table 6: (1) Overall, the value of various ecological services had a significant promoting effect on visitor economy growth. The influence coefficients of overall value, regulation, support, and cultural service values were positive and pass the significance test i.e., these variables have a significant promoting effect on tourism revenues; (2) the value of various ecological services had significant heterogeneity in visitor economy growth. The regression results showed that the value of cultural services had the largest influence coefficient on tourism revenue (60.97), followed by the value of support services (21.72), the value of regulation services (11.94), and the total value (6.84). The land supply service is mainly reflected in the waste disposal, soil formation and conservation, biodiversity conservation, and landscape aesthetic provision functions. Land support services are the core function of land, which have a crucial impact on people's everyday lives and economic development. Cultural service is reflected in the cultural and entertainment functions of land. For example, heritage architecture and sites can be regarded as cultural tourism resources with cultural dissemination and educational value. Thus, different land types have dissimilar impacts on the visitor economy. Note: * represents the significance level of 10%, ** represents the significance level of 5% and *** represents the significance level of 1%.
Thus, different types of land have different ecological service values and have significant differences in tourism development; and (3) other control variables, such as population, education levels, and average salary of employees, all passed the significance test and were consistent with the actual situation.

Conclusions
As the core area of economic development in Hubei Province, the Wuhan metropolitan area is also an important strategic fulcrum for the rise of central China. Its land use transformation is of great significance for social and economic development. Based on the spatial evolution of land use, the dynamic attitude of land use, transfer matrix, and a linear regression model were used to study the evolution characteristics of land use and the impact on tourism development in the Wuhan metropolitan area. The main conclusions are as follows: First, the land use scale in the Wuhan metropolitan area changed significantly from 1990 to 2018. From the perspective of various land use ratios, the land use in the study area is mainly cultivated land, forest land, and water bodies. From the perspective of various land use changes, the area of cultivated land, forest land, and grassland in the study area has been decreasing, while the area of water bodies and construction land has been increasing in the past 30 years. From the perspective of the development stage of land use change, after 2005, the expansion rate of construction land significantly accelerated, and the urbanization rate increased rapidly, which had a severe impact on land use types. With the acceleration of urbanization in the world, the characteristics of land use evolution in the Wuhan metropolitan area provide a useful reference for urban land planning in other countries.
Second, there are great differences in the dynamic attitudes of land use at various stages. On the one hand, in terms of the dynamic attitude of single land use, the area of forest land, grassland, and cultivated land is decreasing annually, and the rate of decrease is increasing. However, the water bodies and construction land area increased annually, the rate of increase was smaller. The unused land area showed an increasing trend initially and then a decreasing trend. Among the six types of land, the dynamic change in construction land use was the largest, while the degree of dynamic change of grassland was the smallest. On the other hand, in terms of the dynamic attitude of integrated land use, from 1990 to 2005, the comprehensive dynamic attitude of land use in the Wuhan metropolitan area was 2.50%, and the rate of land use change was relatively gentle. From 2005 to 2020, the comprehensive dynamic attitude of land-use in the Wuhan metropolitan area was 9.78%, and when compared with the previous period, the rate of change of land use type showed an accelerating trend, and the impact of human activities increased.
Third, based on the perspective of ESV, this study focused on the impact of land use change on the development of the visitor economy by constructing a regression equation model. The results showed that although regulating, supporting, and cultural service values had significant promoting effects on visitor economy growth, their effects varied; cultural service value, support service value, and regulation service value had different influence coefficients of 60.97, 21.72, and 11.94, respectively. In addition, this study introduces the three dimensions of regulatory, support, and cultural services that can promote land ecological service value enhancement and tourism economic growth. This may be new ideas for other countries or regions in formulating tourism economic development policies.
Based on the above findings, the following policy implications were derived. (1) All planning of land use should be carried out to achieve sustainable development of land use. Land use change is the most direct reason for the change of the total value of regional ecosystem services, and whether the change direction is reasonable or not directly determines the total value of regional ecosystem services. Land use planning must give full consideration to the value of land supply services, adjustment services, support services and cultural services [44]. The government should arrange the quantity and spatial layout of all types of land according to local conditions, so as to make the land use in more scientifically rational, and achieve the optimal allocation of land use resources and sustainable use. (2) The comprehensive utilization efficiency of land should be used to solve the contradictions between land supply and demand. With continuous popluation growth and the acceleration of urbanization, the amount of cultivated land decreases rapidly, which leads to increasingly sharp contradictions between man and land [45]. Therefore, the government and land resource management departments should take measures to solve the contradictions between land supply and demand to ensure economic development and food security. On the one hand, efforts in land development and consolidation should be intensified, the construction of basic farmland protection areas stengthened, and the economical and intensive use of land resources promoted. On the other hand, wasteland should be developed and transformed such as sandy land, tidal flats and fields, as well as medium and low producing areas, and the intensity of land pollution control must be strengthened. (3) The value of land ecological services should be improved to lay a foundation for the development of the visitor economy. The value of ecosystem services in the Wuhan metropolitan area shows an increasing trend year on year, which indicates that the ecosystem structure is constantly optimizing and has achieved good ecological effects, which lays a foundation for the development of the visitor economy [46]. However, it is worth noting that at present, the increase of ecosystem service value in the Wuhan metropolitan area mainly comes from the reclamation and utilization of unused land, which is converted into the increment of the value of forest land and water land. In fact, the land use mode and technology are still at a low level. Therefore, the Wuhan metropolitan area should attach importance to the value of land ecological services, improve the ecological environmental benefits per unit of area, and realize the sustainable development of ecosystems and visitor economy.
This study had some limitations, (1) Lack of discussion of the spatial dimension. The impact of land use change on tourism economic development involves time and space dimensions. However, this study only discussed the impact of land use change on tourism from a time dimension, and lacked research on spatial autocorrelation and adjacency. The essence of the evolution of land use patterns is a comprehensive process of its spatial form, spatial functional structure, and spatial element distribution. Therefore, the influences on the visitor economy can be further analyzed in the future from these three aspects. (2) The calculation of land ecological service value needs further study. At present, there is no unified understanding of ecological service value coefficients of various types of land. This paper refers to the methods provided by IPCC and relevant scholars. However, the land ecological service value coefficient, which is more accurate and consistent with the provincial scale in China, still needs further study.  Data Availability Statement: All data required to support this research is already presented in this manuscript.