Research on Land Use Versatility Evaluation, Spatiotemporal Coupling, and Influencing Factors Based on Multimethod Integration

: Land use multifunctionality research is important for the efficient use of land resources and the resolution of land use conflicts. With the use of methods such as the technique for order of preference by similarity to ideal solution (TOPSIS) model, coupling coordination model, and geographical detector, the land use multifunctionality level, spatiotemporal coupling, and influencing factors in the Sichuan Province of China from 2000 to 2020 were systematically analyzed in this paper. It was revealed that, from 2000 to 2020, the comprehensive land use functionality in Sichuan Province was continuously improved with increasing economic, social, and ecological functionality levels. The comprehensive land use functionality in each city (prefecture) exhibited a positive development trend. The coupling coordination degree of the land use multifunctionality in Sichuan Province has been continuously improved, undergoing an evolutionary process from the brink of disarray to barely coordinated, then to primary coordination, and finally to medium coordination. The spatial differentiation of land use multifunctionality coupling coordination among cities (pre-fectures) was notable, showing center–periphery spatial distribution characteristics. The average slope and employed population density exhibited the highest explanatory power for the spatial differences in land use multifunctionality coupling coordination. The interaction between any two factors exerted a greater impact than any single factor on the spatial differentiation of land use multi-functionality coupling coordination. Based on the regional development characteristics, region-specific strategies should be adopted to enhance the land use multifunctionality level in Sichuan Province.


Introduction
Land serves as a vital foundation for human production and living, offering various products and services during its usage, thereby demonstrating the multifunctionality of land use.The study of land use multifunctionality is beneficial for measuring and evaluating the impact on human well-being, as well as aiding in the harmonious development of the national economy, society, and the ecological environment.As societal and economic progress continues, with the acceleration of industrialization, urbanization, and the modernization of agriculture, along with the increasing intensity of human activities [1,2], the scarcity of land resources in underdeveloped areas has become increasingly promi-nent.This scarcity leads to competition, where social development encroaches upon economic development spaces, economic growth divides ecological environmental spaces, and social development areas lack ecological spaces [3].Consequently, investigating land use multifunctionality is important for enhancing the comprehensive benefits of land resources, optimizing the national territorial spatial pattern, facilitating ecological civilization construction, and contributing to the building of a beautiful China [2,4,5].
The study of land use multifunctionality traces its origins to the concept of multifunctionality in agriculture [6][7][8].This term refers to the state and manifestation of the land use functions of a given region, including environmental, economic, and social aspects.This represents the diverse products and services provided by different land use practices and is a crucial concept and methodological framework for assessing the impact of land use changes on sustainability [2,6,9].In recent years, land use multifunctionality has become a focal point of interest among scholars both domestically and internationally [2][3][4][10][11][12].International researchers, relying on the Sustainability Impact Assessment-Tools for Environmental, Social and Economic Effects of Multifunctional Land Use in European Regions (SENSOR) project under the Sixth Framework Programme of the European Union [13,14], have empirically studied land use multifunctionality [14].They have noted that changes in land use functionality can impact regional sustainable development [15].Qualitative and quantitative studies have been carried out in typical regions, including the Mediterranean area [16,17], central Europe [18], northern Spain [19], North America [20], the northeastern part of South Asia in the Ganges River Basin [21], and the Arctic region [22].Domestically, scholars have conducted research at various scales, ranging from national and provincial to municipal levels, as well as finer scales, such as land parcels and grid cells.Previous studies are generally based on typical geographic units such as ecologically fragile zones [2,3], the Xiangjiang River Basin [6], and the middle reaches of the Heihe River [7].Land use multifunctionality can be categorized into economic-social-environmental-cultural [23], economic-social-production-ecological [24], and production-living-ecological [25] dimensions.Various mathematical models have been employed in these studies, including an exploratory spatial data analysis [3], geographically weighted regression [26], the gray relational projection method [23], and a comprehensive diagram method of full permutation polygons [27].Overall, the existing research has achieved certain progress in the quantitative analysis of land use multifunctionality, thereby preliminarily establishing a research framework system.However, there are deficiencies in the scope of study regions and in the construction of theoretical frameworks.
Sichuan, characterized by its typical basin topography, exhibits significant regional differences in land use.Within this context, in this study, Sichuan Province was selected as the research subject, and an evaluative index system was constructed from three functional perspectives: economic, social, and ecological.Via the adoption of methodologies such as the entropy-weighted technique for order of preference by similarity to ideal solution (TOPSIS) model, coupling coordination degree, and a geographical detector and the utilization of analytical software such as ArcGIS (10.3) and SPSS (19.0), this research focuses on examining the spatiotemporal differentiation characteristics of land use multifunctionality and its coupling coordination in the region from 2000 to 2020.Furthermore, the aim was to identify the influencing factors of the spatial differentiation of land use multifunctionality and its coupling coordination.The objective was to provide insights for underdeveloped areas to enhance the efficacy of land resource utilization and optimize the national territorial spatial pattern.

Literature Review
The land use system is a composite ecosystem with specific attributes, and the functionality of land use intuitively reflects this system [15,16,28].The multifunctionality of land use refers to the state and manifestations of the functions of land use in a specific region, including its economic, social, and ecological aspects, that is, the products and services provided to human production and life under different land use methods [15,[29][30][31].Generally speaking, the economic function is the foundation for the healthy development of land use, the social function is the goal, i.e., meeting the needs of human survival and development, while the ecological function serves as the safeguard.All three impact the functionality of land use.Based on different objectives, the subsystems of land use multifunctionality interact and establish a coupled and coordinated relationship in their respective developmental processes.This is characterized by economic functions supporting the enhancement of social functions and the maintenance of ecological functions: social functions, relying on economic development and ecological improvement, promote both, serving as a directive for internal functional optimization and coordination; the ecological function is a key component of sustainable development.The quality of ecological functions directly influences the development of economic and social functions, and the development of socioeconomics also poses severe challenges to environmental protection.In light of this, this study, drawing upon systems theory and the concept of sustainable development, constructs an evaluation index system for economic, social, and ecological functions with three dimensions: economy, society, and ecology (Figure 1).This system is used to comprehensively assess the level of land use multifunctionality.Additionally, by employing a coupling coordination model and integrating a geographic detector, this study calculates the internal coordination state of land use functions and comprehensively identifies influencing factors.This approach aims to further promote regional sustainable development, enhancing land use efficiency, and resolving conflicts in land utilization.

Overview of the Research Area
Sichuan Province, situated in the inland southwest of China, exemplifies a typical underdeveloped region in the western part of the country, characterized by notable poverty in the Liangshan area to the west and the Qinba mountainous belt to the east.Located in the upper reaches of the Yangtze River, the province is characterized by significant topographical diversity and a complex geomorphology, featuring a west-to-east descending terrain.The western region comprises plateaus and mountains, the central area consists of plains, and the eastern region is characterized by basins and hills.The province is rich in river systems, predominantly the Yangtze River, and is largely covered with purple soils.The forest coverage rate in Sichuan reaches 38.83%, with an average annual temperature of approximately 17 °C, an annual precipitation of 1193 mm, and an average annual sunlight period of 1443.89h.Spanning an area of 48.61 km², the province encompasses 18 prefecture-level cities and 3 autonomous prefectures.It exhibits an average population density of 188 people/km², which is 1.27 times higher than the national average value of 148 people/km², and an average economic density of 999.87 million renminbi/km², equivalent to 89.56% of the national average value of 1058.32 million renminbi/km².Sichuan faces acute and complex challenges in terms of human-land relations.Therefore, selecting it as the subject for a study of land use multifunctionality in underdeveloped southwestern regions is both typical and representative (Figure 2).

Data Sources
The primary sources of raw data include the China Statistical Yearbooks (https://www.stats.gov.cn,accessed on 15 November 2022) and Statistical Yearbooks of Sichuan Province and its various cities (prefectures) (http://tjj.sc.gov.cn,accessed on 15 November 2022), the National Economic and Social Development Bulletins of Sichuan Province and its cities (prefectures) (https://www.sc.gov.cn,accessed on 15 November 2022), and the Water Resources Bulletins of Sichuan Province and its cities (prefectures) (http://slt.sc.gov.cn,accessed on 10 December 2022), all for the 2000-2020 period.Additional data were sourced from official websites and other relevant databases.The data needed to calculate the influencing factors such as the average elevation and slope were obtained from the Geospatial Data Cloud of the Chinese Academy of Sciences (https://www.gscloud.cn,accessed on 8 February 2023).Slope data were generated using ArcGIS (10.2) and the digital elevation model (DEM).

Construction of the Index System
Land use multifunctionality evaluation is a complex systemic endeavor, but a unified standard for assessing land use multifunctionality has not yet been established in the academic community.Drawing on existing research [2,5,9,32,33] and considering the unique characteristics of Sichuan Province, in this paper, an evaluation index system for land use multifunctionality was constructed.This system was developed in accordance with the principles of completeness, scientific rigor, systemic coherence, and data availability.It comprises three primary indicators, namely, economic, social, and ecological functions, thereby providing a comprehensive framework for assessing land use multifunctionality (Table 1).The economic function refers to the capacity to provide material resources and infrastructure for human production and living.It primarily includes agricultural production and economic development functions.Per capita grain possession, fruit yield per unit area, and the proportion of the total output value of agriculture, forestry, animal husbandry, and fishery in GDP were selected to represent the agricultural production function.The proportion of secondary and tertiary industries and per capita GDP were used to characterize the economic development function.The social function addresses the human pursuit of a better life and the needs of human production and living, encompassing spatial carrying capacity, social security, employment support, and transportation guarantee functions.Population density was chosen to reflect the spatial carrying capacity; the number of medical beds per 10,000 people, and the average wages of employees, and the urban-rural income balance index represents the social security function.Road network density illustrates the transportation guarantee function, while employment density demonstrates the employment support function.The ecological function refers to the ability to maintain the natural conditions and utilities essential for human survival, mainly including environmental purification and resource supply functions.Sewage treatment rate and fertilizer application per unit area were selected to represent the environmental purification function; the land reclamation rate and average meat production of pigs, cattle, and sheep per unit of land illustrate the resource supply function.

Resource Provision Functionality
Land reclamation rate (%) Total cultivated land area per regional total area Per capita pork, beef, and mutton production (ten thousand tons/km²) Pork, beef, and mutton production per regional total area

Determination of Indicator Weights
The range method was used to standardize the 18 original data indicators [34], and the entropy method [35] was employed to calculate the weights of the indicators.
(1) Data standardization processing and coordinate translation, where B represents the translation distance, which was 0.01 in this study.
Positive indicator: Negative indicator: (2) Calculate the proportion of the value of indicator j in city i, as follows: (3) Calculation of indicator information entropy and its diversity coefficient: In the formula, D = 1/LN n , where n represents the number of samples.
(4) Calculation of indicator weights: In the formula, Xij represents the value of the jth evaluation indicator for city (or state) i; min{Xj} and max{Xj} denote the minimum and maximum values of the jth indicator across all cities (or states).D = 1/LN n , where m is the number of cities (or states) being evaluated, and n is the number of indicators.

Measurement of Land Use Multifunctionality Level
This study employs the MATLAB R2016 software and utilizes the TOPSIS model to evaluate the merits and demerits of land use multifunctionality.The TOPSIS model is a comprehensive evaluation method for multiobjective decision-making [36].Its core concept involves establishing an ideal model and determining the proximity of the evaluation targets to the optimal and worst solutions.This approach can be used to assess the relative merits of the evaluation targets, allowing for both a horizontal evaluation and longitudinal comparison.This method is also suitable for evaluating small samples and can be applied in a comprehensive evaluation of large systems with multiple objectives [37].The specific implementation steps are as follows: (1) Constructing the decision matrix: (2) Determining the positive and negative ideal solutions: where  + and  − denote the positive and negative ideal solutions, respectively, and rij is the standardized value of the R matrix.
(3) Calculating the Euclidean distance: this refers to the distance of the ith (i = 1, 2, 3, …, n) evaluation object from the best and worst solutions: where   + and   − denote the distances of the evaluation object from the best and worst solutions, respectively, and   is the entropy weight.
(4) Calculating the relative closeness of each evaluation object to the optimal and suboptimal solutions: where   occurs within 0≤   ≤1.A   value closer to 1 indicates superior land use multifunctionality, while a value closer to 0 indicates inferior multifunctionality.

Analysis of Land Use Multifunctionality Coupling Coordination
The synergy among the various functions of land use can be measured using a coupling coordination model [5,23], which can be expressed as follows: where C is the coupling degree, and (), (), and ℎ() denote the indices of the economic, social, and ecological functions, respectively.The symbols α, β, and λ denote the respective weights.To provide a qualitative classification standard for the coupling coordination development coefficient, in this study, in alignment with existing research [23,38], the coupling coordination coefficient was categorized into three major classes and 10 subcategories (Table 2).

Geographical Detector
The geographical detector, proposed by Wang Jinfeng and Xu Chengdong [39], is a statistical method for detecting and utilizing spatial heterogeneity and revealing potential driving factors.This method is widely used to analyze the spatial pattern evolution and spatial differentiation of geographical elements [38,39].In light of the above, factor detection and interaction detection in the geographical detector toolkit were primarily utilized in this study.Factor detection was mainly used to analyze the explanatory degree of the various influencing factors of land use multifunctionality.
where  is the degree of explanation for the dependent variable, with a range of [0, 1], denotes the stratification of either the independent or dependent variable, and  ℎ and  denote the number of units in stratum h and the total number of units in the entire region, respectively.Furthermore,  ℎ 2 and  2 denote the variances in the Y values in stratum h and the entire region, respectively, and  and  denote the sums of squares of the within-stratum variance and the total variance in the entire region, respectively.
Interaction detection aims to identify the interaction effects of the different influencing factors, which involves detecting whether the interaction between X1 and X2 can enhance or diminish the explanatory power of the dependent variable.There are a total of five types of factor interactions: single-factor nonlinear attenuation, nonlinear attenuation, dual-factor enhancement, mutual independence, and nonlinear enhancement.

Construction of Impact Factor Indices
To investigate the factors influencing the spatial differentiation of land use multifunctionality coupling coordination in Sichuan Province and to further enhance the land use multifunctionality level, in this study, a geographical detector was utilized.An index system of the influencing factors for the spatial differentiation of land use multifunctionality coupling coordination was constructed, encompassing natural, social, and economic dimensions (Table 3).Gross domestic product/total area

Evaluation of the Land Use Multifunctionality in Sichuan Province
From 2000 to 2020, the economic functionality of the land use in the entire province continuously increased.The index increased from 0.3792 in 2000 to 0.6865 in 2020, with an average annual increase of 3.01% (Figure 3).The change process can be divided into two periods: an initial decline followed by a consistent increase.Specifically, the 2000-2005 period marked the initial decline stage, where the index decreased from 0.3792 in 2000 to 0.3138 in 2005, with an average annual decrease of 3.71%.The 2005-2020 period was characterized by a continuous increase, with the index increasing from 0.3138 in 2005 to 0.6865 in 2020, at an average annual growth rate of 5.36%.During the 10th Five-Year Plan of China, the province underwent adjustments in the agricultural industrial structure, with a continuous increase in the urbanization rate and a relative decrease in the share of the total grain output and the output value of agriculture, forestry, animal husbandry, and fishery in the GDP.After 2005, the province continued to promote the transformation of its economic development mode, adjusted its economic structure, and enhanced the economic efficiency, leading to rapid increases in the added value of the secondary and tertiary industries, per capita GDP, and economic agglomeration, thereby continuously enhancing the economic functionality.
The social functionality of the land use in Sichuan Province continued to improve, with the index increasing from 0.3027 in 2000 to 0.6833 in 2020, with an average annual increase of 3.95% (Figure 3).The process of change can be divided into two stages: slow growth and rapid growth.The 2000-2010 period represented the slow growth stage, with the index increasing from 0.3027 to 0.3888.The 2010-2020 period constituted the rapid growth stage, with the index increasing from 0.3888 to 0.6833.Since the implementation of the national 10th Five-Year Plan and the formal initiation of the Western Development Strategy, along with the implementation of the poverty alleviation campaign and the subsequent introduction of the Rural Revitalization Strategy, the provincial infrastructure continued to improve, the social security system was increasingly refined, medical and health services steadily advanced, and wages consistently increased, driving a continuous increase in the social functionality.
The ecological functionality of the land use in Sichuan Province exhibited a cyclically fluctuating upward trend, with the index increasing from 0.2656 in 2000 to 0.5378 in 2020, an average annual increase of 3.42% (Figure 3).The change process can be divided into three stages: increase, decline, and increase again.The 2000-2012 period was characterized by a rapid increase, with the index rapidly increasing from 0.2656 to 0.4534.A slight decline occurred between 2012 and 2013, during which the index decreased from 0.4534 to 0.4390.The 2013-2020 period exhibited slow growth, with the index increasing from 0.4390 to 0.5378.The ecological environment is a crucial marker of social and civilizational progress, and the advancement of green and ecological development is a key aspiration when aiming to improve the quality of life for the population.Following the proposal in 2010 to establish a long-term green development mechanism, the 12th and 13th Five-Year Plans introduced corresponding measures for ecological and environmental protection, thereby facilitating positive development in the ecological functionality.Overall, from 2000 to 2020, the comprehensive functionality of the land use in Sichuan Province demonstrated an upward trend, with the index increasing from 0.3245 in 2000 to 0.6324 in 2020, marking an average annual increase of 3.39% (Figure 3).This growth was primarily attributable to the rapid development of the economic functionality of land use, the continuous refinement of the social functionality, and the sustained improvement in the ecological functionality in the province.From a current perspective, the ranking of the various functionalities is economic > social > ecological, indicating clear differences.In particular, the ecological functionality of land use still exhibits considerable room for enhancement, which also suggests that the coordinated development among the various functionalities requires further strengthening.

Evaluation of the Land Use Multifunctionality in the Various Cities (Prefectures)
During the 2000-2020 period, the land use functionality level varied across different cities (prefectures).In 2020, high-value economic functionality areas were found in Chengdu, Ziyang, and Meishan, with indices of 0.5529, 0.5428, and 0.5316, respectively.Chengdu, as the core city of the Chengdu-Chongqing economic circle and an important central city in the western region, has experienced rapid economic development.Additionally, the region's flat terrain and fertile soil, coupled with the polarization-trickledown effect observed in the capital city, are particularly important.Ziyang and Meishan, both parts of the Chengdu Plain economic zone, are significantly influenced by the capital city's trickle-down effect, thereby positioning the economic functionality of land use at the forefront in the province.Low-value areas were identified in the Ganzi Tibetan Autonomous Prefecture and Aba Tibetan and Qiang Autonomous Prefecture, with indices of 0.3093 and 0.3299, respectively.Despite the abundant tourism resources, these regions exhibit a low gross regional product, remote geographical locations, a harsh natural environments, and an underdeveloped transportation infrastructure, all of which hinder economic development and result in the lagging economic functionality of land use (Figure 4).
In terms of the social functionality of land use, Chengdu city is far ahead of other cities (prefectures), with an index reaching 0.6719.As the provincial capital, Chengdu is the political, economic, and cultural center of Sichuan Province.Its high spatial carrying capacity, favorable employment and transportation conditions, and continuously improving social security system have significantly driven the rapid development of its land use social functionality.The social functionality of land use in Ganzi and Liangshan ranks at the bottom across the province.This is mainly because both areas are highly impoverished, characterized by a prevalence of mountains and scarcity of plains, which enhances the conflict between human and land resources.Additionally, the low economic level and difficulties implementing social security and other related infrastructure have resulted in a low social functionality of the land use in these regions (Figure 4).
High-value ecological functionality areas are located in Bazhong and Ya'an, with indices of 0.5630 and 0.5538, respectively.Bazhong serves as a crucial ecological barrier in the upper reaches of the Yangtze River and the Qinling Mountains, leading in the province with high forest coverage and wastewater treatment rates, among other ecological maintenance and environmental purification indicators.Ya'an, referred to as the Lung of Tianfu, benefits from its unique natural resources.Coupled with the government's continuous efforts in strengthening ecological and environmental management, these factors have positioned Bazhong and Ya'an as the top two provinces in terms of the ecological functionality of land use.A low-value area was identified in Liangshan Prefecture, with an index of only 0.4866.This is mainly attributed to the rapid population growth and severe soil erosion, exacerbating ecological degradation and resulting in the relatively low ecological functionality of land use in this region (Figure 4).
Overall, the economic, social, and ecological functionalities are continuously improving across the various cities (prefectures).The absolute differences in the comprehensive land use functionality among these cities (prefectures) are declining, with the range decreasing from 0.2707 in 2000 to 0.1887 in 2020.However, the comprehensive functionality indices of the land use in these cities (prefectures) remain relatively low, indicating that there is substantial room for further improvement (Figure 4).

Temporal Evolution Analysis of the Coupling Coordination Degree of the Land Use Multifunctionality
From 2000 to 2020, the coupling coordination degree of the land use multifunctionality in Sichuan Province was continuously improved, with the coordination coefficient increasing from 0.4415 in 2000 to 0.7416 in 2020, marking an average annual increase of 2.63%.The state of coupling coordination evolved from the brink of disarray to a medium level of coordination, transitioning from a nonideal type to an ideal type of coupling coordination state.The process encompassed three stages, namely, nonrational, transitional, and ideal stages, experiencing an evolutionary process from the brink of disarray to barely coordinated, then to primary coordination, and finally to medium coordination.The year 2010 represented a pivotal point, with the land use multifunctionality at a barely coordinated state.Prior to 2010, only the index of the economic functionality of land use was above 0.4, while the indices of the social and ecological functionalities were below 0.4, indicating a lower level of development and hindering the enhancement in the overall effectiveness of land use.After 2010, the comprehensive evaluation value was improved, propelling the land use multifunctionality to a state of primary coordinated development.The coupling of the economic, social, and ecological functionalities of land use was favorable, with the overall functionality trending toward coordinated development (Figure 5).

Spatial Evolution Analysis of the Land Use Multifunctionality Coupling Coordination Degree
Examining the changes in the coupling coordination status of the land use multifunctionality in the various cities (prefectures), it could be observed that the coupling coordination status was synchronously improved across these regions, with all indices increasing, although the extent of improvement varied, ranging from 0.05 to 0.33.Cities with an increase in the coupling coordination index between 0.05 and 0.10 include Luzhou, Deyang, Suining, Leshan, Meishan, Guang'an, Dazhou, Ya'an, Bazhong, and Ziyang, accounting for 47.62% of the province.This mainly occurs because these cities initially exhibited a relatively favorable coordinated development among the economic, social, and ecological functions, hence the smaller magnitude of the increase in their coupling coordination coefficients.Cities (prefectures) with an index increase between 0.10 and 0.20 include Chengdu, Zigong, Mianyang, Guangyuan, Neijiang, Nanchong, Yibin, and Liangshan, accounting for 38.10% of the province.In these areas, the economic functionality, represented by the economic density and per capita GDP, showed significant growth.In terms of the social security functionality, the average wage of employees increased from 7600 to 70,471 renminbi.Environmental purification indicators, such as the wastewater treatment rate and fertilizer usage, also showed positive trends, effectively enhancing the degree of coupling coordination of land use multifunctionality.Cities (prefectures) with an index increase between 0.20 and 0.33 include Panzhihua, Aba, and Ganzi, accounting for 14.29% of the province.In 2000, these areas were categorized as exhibiting a disordered state of coupling coordination.However, the implementation of poverty alleviation campaigns and rural revitalization strategies, continuous improvements in infrastructure, economic development, social security system improvements, and ecological environment protection in western Sichuan significantly enhanced the coupling coordination of the comprehensive land use functionality.The range of the coupling coordination coefficients of land use multifunctionality among the various cities (prefectures) continuously declined, decreasing from 0.4140 in 2000 to 0.1672 in 2020, with an average annual decline of 4.43%.This is primarily due to the synergistic advancement and steady improvements in the economic, social, and ecological functionalities in each city (prefecture), facilitating the evolution of land use multifunctionality coupling coordination across the province (Table 4).During the 2000-2020 period, the spatial differentiation of land use multifunctionality coupling coordination was notable across the various cities (prefectures).In 2000, the spatial pattern exhibited a structure comprising one core, two zones, and multiple points.The core was the Chengdu Plain economic zone, where Chengdu, Deyang, Meishan, Suining, and Ziyang all showed a state of primary coordination.The two zones are the Southern Sichuan economic zone and the northeastern Sichuan economic zone, with Zigong, Luzhou, Dazhou, Guang'an, and Bazhong in these zones also showing a state of primary coordination.The multiple points were distributed across 11 cities (prefectures) throughout the province, accounting for 52.38% of the total province, where the coupling coordination coefficients of the land use multifunctionality were relatively low.By 2020, all cities (prefectures) had reached coordination, forming a center-periphery spatial distribution structure that was primarily centered on the Chengdu Plain economic zone and decreased toward the periphery.In 2020, there were eight cities in a state of medium coordination, accounting for 38.10% of the province, exhibiting a 521 distribution pattern, with five cities in the Chengdu Plain economic zone, two cities in the Southern Sichuan economic zone, and one city in the northeastern Sichuan economic zone.There were 10 cities in a state of primary coordination, accounting for 47.62% of the province, mainly distributed in the Panxi economic zone and the northeastern Sichuan economic zone.Three cities (prefectures) showed a barely coordinated state, accounting for 14.29% of the province, namely, Aba, Ganzi, and Ya'an (Figure 6).

The interpretation of the Findings
In recent years, significant progress has been made in the field of land use multifunctionality research [40][41][42].Notably, there has been a significant increase in studies focusing on the driving factors of internal functional coordination within land use multifunctionality [43][44][45].This trend correlates with the intensity of China's policies over the past decade, aimed at poverty alleviation, rural revitalization, and sustainable rural development.With the rapid socioeconomic development, the internal differences in land use functions are becoming increasingly evident, with the hierarchy being economic functions > social functions > ecological functions [46].This is primarily due to the continuous implementation of economic development policies across the province, leading to a consistent elevation in economic levels, with the average economic density reaching 999.87 million CNY/km 2 , progressively narrowing the gap with the national average of 1058.32 million CNY/km 2 .The unique regional location and typical basin topography of Sichuan Province have resulted in a diverse range of geomorphological types.Generally, the conditions in the central basin differ from those in the surrounding mountainous areas, leading to a marked spatial differentiation in the coupling and coordination of land use multifunctionality, characterized by higher levels in the central region compared to the periphery [47,48].Owing to its location in Western China and significant geomorphological features, the primary factor influencing the spatial differentiation in the coupling and coordination of land use multifunctionality is the average slope [49].Moreover, apart from the distinct geomorphological characteristics, the western region of Sichuan, being an ethnic area with relatively recent economic development, contributes to the varied influencing factors across different municipalities and prefectures.This also indicates that the interaction between any two influencing factors has a greater impact on the coupling and coordination of land use multifunctionality than any individual factors [50,51].

Policy Recommendations
To comprehensively promote the coordinated development of land use multifunctionality and effectively enhance the efficiency of land resource utilization while reconciling human-land conflicts, Sichuan Province should rationally optimize its territorial spatial development pattern.In line with the high-quality development requirements oriented towards ecological priority and green development, the province should expedite the preparation of provincial territorial spatial planning based on the results of the "dual evaluation."It should adhere to a bottom-line thinking approach, and, considering regional developmental characteristics, promote economic development in a manner suited to local conditions.Underdeveloped areas should particularly integrate and coordinate the relationships among economic, social, and ecological functions.They should balance and coordinate the input of various elements within the economic-social-ecological system, thereby effectively enhancing the level of land use multifunctionality.This approach will continuously improve the benefits of land resource utilization.Moreover, underdeveloped areas should also comprehensively plan and manage these aspects, establishing land use planning and management mechanisms that incorporate factors like economic development, ecological protection, and social security, to ensure the coordinated development of various functions.Concurrently, strengthening land use monitoring and assessment is crucial to adjust planning and management strategies in a timely manner.

Future Directions for Research
This paper focuses on exploring the multifunctionality of land use at the provincial and municipal (prefecture) scales in Sichuan Province, although it acknowledges certain limitations.First, future research will aim to further perfect the evaluation indicator system.Based on regional spatial characteristics, the system will be refined to emphasize research on the multifunctionality of land use at the micro-scales of counties, townships, and villages throughout the province.This research will uncover the spatiotemporal characteristics and influencing factors of land use multifunctionality at these scales, contributing to the sustainable development of land resources in the province.Secondly, there will be an innovation in the methods used to evaluate multifunctionality.The research will fully utilize and apply novel evaluation methods, such as remote sensing technology, Geographic Information Systems (GIS), and model simulations, to enhance the accuracy and efficiency of evaluating land use multifunctionality.Thirdly, the study will focus on effectively linking the evaluation of multifunctionality with sustainable development goals.It will explore how to achieve sustainable development in economic, social, and ecological terms through the optimization of land use, effectively integrating land use multifunctionality evaluations with sustainable development objectives.

Conclusions
In this study, based on an exploration of an evaluation index system for the land use multifunctionality in Sichuan Province, the entropy-weighted TOPSIS model and coupling coordination degree were used to investigate the development level and the spatiotemporal differentiation characteristics of the coupling coordination degree of the land use multifunctionality in Sichuan Province and its 21 cities (prefectures).With the use of a geographical detector, the factors influencing the spatial differentiation in the coordinated development of land use multifunctionality were also detected.The conclusions are as follows: (1) The comprehensive functionality of the land use in Sichuan Province continues to improve, with significant progress in the economic, social, and ecological functionalities.Overall, the functionalities can be ranked as follows: economic > social > ecological.Among the various cities (prefectures), the level of each functionality varies, with clear differences in the extent of improvement.The absolute differences in the comprehensive land use functionality across cities (prefectures) are decreasing, with the range declining from 0.2707 in 2000 to 0.1887 in 2020.(2) The coupling coordination degree of the land use multifunctionality in Sichuan Province is continuously improving, with the coordination coefficient increasing from 0.4415 in 2000 to 0.7416 in 2020.There is clear spatial differentiation in the coupling coordination degree of the land use multifunctionality among cities (prefectures).In 2000, a spatial pattern comprising one core, two zones, and multiple points was evident, while in 2020, all cities (prefectures) entered a state of coordination, forming a center-periphery distribution pattern.(3) The main single factors influencing the existing spatial differentiation in the coupling coordination of land use multifunctionality are the average slope and the employed population density, with the explanatory power of both factors exceeding 0.9.The interaction between any two influencing factors exerts a greater impact on land use multifunctionality coupling coordination than a single factor, indicating that the spatial differentiation in the coupling coordination development of the land use multifunctionality in Sichuan Province is influenced by multiple factors.
Author Contributions: N.Z.designed the experiments and drafted the manuscript.Data processing and visualization were performed by X.X., G.J. and Y.G. revised the manuscript.All authors have read and agreed to the published version of the manuscript.

Figure 3 .
Figure 3. Variation trend of the land use versatility index in Sichuan Province.

Figure 4 .
Figure 4. Spatial and temporal evolution of the land use multifunctionality evaluation results for the cities (prefectures) in Sichuan Province.

Figure 5 .
Figure 5. Change trend of the land use multifunctional coupling coordination degree in Sichuan.

Figure 6 .
Figure 6.Evolution of the spatial pattern of land use multifunctional coupling coordination in Sichuan.

Table 1 .
Land use function categories and indicators.
Yield of fruit per unit area (tons/km 2 ) Fruit production per regional total area Proportion of the total output value of agriculture, Total output value of agriculture, forestry, animal husbandry, and fishery/GDP forestry, animal husbandry, and fishery in the gross domestic product (GDP) (%)

Table 2 .
Coupling coordination degree division criteria.

Table 3 .
Index system of the influencing factors.

Table 4 .
Coupling coordination degree of the land use multifunctionality in cities (prefectures) in Sichuan Province from 2000 to 2020.