Evaluation of Sustainable Regional Development Combining Remote Sensing Data and Ecological Constraints: A Case Study of Chaohu Basin, China

The maintenance of sustainable development in light of climate change is presenting society with a serious challenge. This study analyzes the spatial environment and its coupled relationship with socio-economic and food security in the Chaohu Lake Basin in central Anhui Province, China. Previous studies have shown that this study area has good environmental conditions. However, the ecological security, economic development, and food security in the Chaohu Lake Basin was generally poor and poorly coordinated in 2018, especially in the center of local cities and surrounding areas. Nevertheless, many towns are in a well-organized state, meaning sustainable and coordinated regional development can be achieved through resource optimization. Generally, the economic conditions in the Chaohu Lake Basin should be improved in a way that ensures that the environment is not significantly degraded. The main functional zoning requirements of the counties’ land management planning designate most of the counties as parts of ecological conservation areas. The present study provides monitoring of important basin ecosystems and serves as a guide to future sustainable regional development by providing directions and paths to follow; it can also provide a reference for eco-functional zoning that is applicable to various natural basins worldwide.


Introduction
Global environmental change and sustainable regional development are topics of increasing relevance as human modification of the global landscape continues unabated. As stated in the 2019 United Nations Global Sustainable Development Report, the terrestrial ecosystem of China still faces significant challenges. Implementing sustainable land-use practices requires integrating agricultural and climatic interventions. Current land use patterns are not meeting people's needs while agriculture destroys forests and biodiversity, squanders water, and releases one-quarter of global greenhouse-gas emissions. Since the late 1980s, the economy of China has been developing rapidly [1], while at the same time economic development has also caused damage to the environment. The environment of China is fragile in many regions, and the trend of environmental deterioration has not been effectively addressed in recent years. In 1999, pilot projects for the "Grain for Green" program were implemented in some regions of China, and the program was introduced nationwide [2] starting in 20 January 2003. So far, more than 33.80 million hm 2 of farmland has been returned to forests or grasslands. Because the area of cultivated land has gradually decreased in recent years, food security is also an important topic. Economic development always conflicts with the development of sustainable agricultural production and an ecologically sound civilization to a certain extent; knowing how to coordinate the relationships in an "ecological-economic society" can provide a powerful way to promote sustainable regional development.
After China's reform and opening up processes, with the rapid increase in population and land [3][4][5] urbanization of China, a larger proportion of the rural populations migrated to urban areas for various reasons [6][7][8]. The effects of this change in population on ecosystem function have been considered in the analysis of the effects of sustainable development of human-natural systems. Attention needs to be paid to minimizing negative ecological impacts and protecting the visual quality of the environment [9,10]. The rational allocation of resources including water, land, and human resources integrates the social and ecological impact of these changes and is involved in the initiative to modernize the economy [11,12]. Rapid economic development is a double-edged sword for environmental security [13,14] because improving the environment requires investment [15] while economic development negatively affects ecological resources [16]. Increasing agricultural production involves the almost irreconcilable problem of supporting economic growth while conserving the environment and working toward the goal of sustainable development [17,18]. Changes in polices and administrative jurisdictions [19] can promote local economic development indirectly affect the environment. While developing its economy in recent years, China has also attached great importance to protecting the environment as one of the countries worldwide that has initiated policies that support ecological compensation. The principle of ecological compensation is based upon the concept of "those who develop are those who protect; those who destroy are those who restore; those who benefit are those who compensate; those who pollute are those who pay for pollution." Promoting local economic development will indirectly affect the environment. While developing its economy in recent years, China has also attached great importance to the protection of the environment, making it of the main countries worldwide that initiated ecological compensation. The principle of ecological compensation is based upon the idea that "who develops, who protects; who destroys, who restores; who benefits, who compensates; who pollutes, who pays" [20]. At the current stage of development, China has begun to implement a plan for the rehabilitation of cultivated land, grasslands, rivers, and lakes and has conducted pilot projects in some regions and watersheds. However, the system designed to support and protect cultivated land has lagged behind, which seriously affects the advancement of cultivated land protection. It is important to coordinate the relationship between the socio-economic (SE) development of the basin and the development and use of resources, and to manage the relationships involved in resource use between regions and departments.
Under the Sustainable Development Goals (SDGs) of the United Nations, conducting a multi-dimensional sustainability assessment is necessary during regional development. Many scholars use the coupling function to study the relationship between the SDGs, while the direction of sustainable development involves more than two subsystems [21][22][23]. Complexities in sustainability evaluations of subsystems will reflect multiple roles and these roles need to be integrated [24][25][26]. It is interconnected and holistic, and previous studies lack systematic and comprehensive holistic consideration. Scholars have explored ways of linking the SDGs, the degree of interaction [27], and the network of relationships [28]. System dynamics and evaluation methods have been used in these integration frameworks. Fu [29] connects geography and SDGs, presented a classification-coordination-collaboration approach for advancing the SDGs, and found that the quantitative output results of different research scales vary greatly. Also, remote sensing technology can be used during regional sustainability assessment and monitoring [30]. In this entire process, however, it is difficult to feedback the spatiotemporal changes of sustainable development simultaneously with these methods. A remote sensing ecological index (RSEI) can be applied to environmental monitoring in urban (impervious surface) [31,32], basin [33], and agricultural engineering areas [34]. It is necessary to explore the connections between subsystems, as well as to show their spatial differences.
Hefei, the central city of the Chaohu Lake Basin, has grown into a new first-tier city; most counties in the basin have been included in the Yangtze River Delta urban agglomeration and have experienced rapid economic development in recent years. Chaohu Lake is the fifth largest freshwater lake in China. In 2001, the administrative division of Chaohu City was revoked; the administrative divisions had been adjusted from county to city status in 1999. The entire water area was placed under the unified jurisdiction of Hefei City after 2001, and the area around the Chaohu Lake Basin experienced rapid urbanization [35].
The Chaohu Lake Basin had been a typical important grain producing area for the last several centuries. It is an important freshwater lake basin and is part of the New Yangtze River Delta urban agglomeration [36]. Studying the relationship between changes in the natural environment of the basin and the SE situation can help in monitoring important basin ecosystems and serve to provide direction and a path for future sustainable regional development. Managers of the Chaohu Lake Basin are implementing spatial ecological planning in the entire basin. Therefore, this area is a typical research object and the present study can also be used as a reference for eco-functional zoning that is applicable to various natural basins worldwide. This study constructs an eco-socio-economic sustainable development system coordination model, make an empirical study with the statistical data of Chaohu Basin, find the development status of regional subsystems, and study the relationships between subsystems. The aims of this research are based on developing an understanding of the status quo and relationships, to make the subsystems be optimized, and even promote regional sustainable development.

Study Area
The Chaohu Lake Basin is located in Jianghuai Hills of central Anhui Province, between the Huaihe and Yangtze rivers and is a typical traditional agricultural planting area. It ranges between 30 • 58'00"-32 • 58'00" N and 116 • 24'30"-118 • 30'00" E and covers an area of 14,203 km 2 (including the 770 km 2 of Chaohu Lake) [34], adjacent to the Yangtze River in the southeast, the Dabie Mountains in the west, the Jianghuai watershed in the north, and the Chuhe River valley in the northeast. Considering the availability of the data, the counties involved in the Chaohu Lake Basin were included in the study area (Yuexi County only contains one town, Yaohe town), covering five cities (Hefei, Lu'an, Anqing, Wuhu, and Ma'anshan), 19 counties and districts, and 161 towns and streets, with 33 rivers in the entire basin ( Figure 1). Chaohu Lake Basin was the first national tourism and leisure zone approved by the National Tourism Administration, and also one of the first pilot demonstration areas for the development of an ecological civilization in China in which construction is prohibited within 1 km of a red line around the Lake; this Basin supports a population of about 10.75 million people.

Data Sources
The present study employed a long time series of high resolution (30 m) Landsat images. As inputs into the spatiotemporal fusion models, Landsat 8 Operational Land Imager images based on the least amount of cloud cover (less than 8.92%) were used. Nine Landsat images (Path 120-121/Row 38) acquired in 2018 were downloaded from the US Geological Survey website (https://earthexplorer. usgs.gov); information related to the specific images is shown in Table 1. Because the downloaded Landsat data were level 1T, the data have been geometrically corrected. Using ENVI 5.3 software, first, radiant calibration of each phase of each image was conducted, so the DN value of each image was converted into the reflectivity of the sensor; next, the Fast Line-of-sight Atmospheric Analysis of Spectral Hypercubes (FLAASH) atmospheric correction tool was used to perform atmospheric correction on the visible and infrared bands of each phase to eliminate the effects of the atmosphere and the influence of factors such as reflections from some objects; then, the two images were merged by using the seamless mosaic and cubic convolution resampling method, while the annual mean values for each township were calculated using principal component analysis. Finally, ArcGIS 10.2 software was used to crop all of the data based on the boundary of the study area. The SE data was acquired from statistical yearbooks (http://tjj.ah.gov.cn/) and government websites (http://www.ah.gov.cn/).

The Eco-Economic-Social Sustainable Development System
The essence of human society is a complex ecosystem formed by human economic production or the activities of daily life, social development, and changes, resource development and use, and changes in the natural environment. The composition of the composite ecosystem has three main aspects: social, economic, and natural. The system can comprehensively reflect the status of regional development, population status, environmental status and changes, and resource status; it can reflect comprehensive development in a certain area and the impact of the status quo on the composite regional environment. The status of regional economic development, social development, resource status and use, and environmental conditions directly affect the regional economic, social, resource, and environmental systems. By analyzing the complex regional environment and integrated ecosystems, the concept of an Integrated Ecosystem can be employed for interpretation. The three-dimensional development of the environment, SE, and agricultural production need to be integrated and coordinated to achieve sustainable regional goals and requirements. The ecological subsystem helps to maintain regional ecological safety, the SE system mainly creates economic benefits to provide funding for regional development, and the agricultural production system helps to guarantee human food security (FS). The relationship between them is shown in Figure 2. The analysis framework shows that the interconnected and holistic subsystems have connections and interactions. The purple dashed lines subcategorize the structure of the eco-socio-economic sustainable development system; the green dashed lines indicate the dominant functions and sustainable development targets of various types of sub-systems in the regional development; the black solid lines indicate the interactions among the sub-systems; the red solid lines represent the interaction mechanism within the eco-socio-economic sustainable development systems.

Remote Sensing Ecological Index
A comprehensive index is more representative than a single index. The response of environment factors in the Chaohu Lake Basin can be specifically reflected by indicators such as humidity (x 1 ), greenness (x 2 ), heat (x 3 ), and dryness (x 4 ). Among them, the distribution of soil moisture and water resources causes changes in humidity. Vegetation types, soil fertility, and water quality, which are closely related to plant growth, cause changes in greenness. Soil temperature, surface cover, and other causes changes in heat; changes in dryness are caused by soil texture, soil temperature, land degradation, or off-farm construction activities. In the technical implementation of a comprehensive index, the enhancement of thematic information with technology can be used to extract the representation information of relevant indicators from remote sensing images, such as the humidity component of heat transformation, a vegetation index, and land surface temperature (LST) representing humidity, greenness, and heat, respectively. Because buildings are an important part of the human ecosystem, the impermeable surfaces of buildings and infrastructure replace the original natural ecosystem, leading to a "drying" of the earth's surface. Therefore, the normalized difference built-up and soil index (NDBSI) represents the "dryness" of the landscape. In this study, a remote sensing ecological index (RSEI) evaluation model coupled with the above indices was used to monitor and evaluate the environmental quality (EQ) of a land reclamation project area. Humidity, greenness, heat, and dryness were represented by a humidity component, a normalized vegetation index, surface temperature, and dryness index, respectively; these indicators reflect changes in EQ. An RSEI is not suitable for classifying large areas of water. If a large area of water exists in a research area, it is better to mask the area before applying an RSEI; a water area is the region with the best environment by default. Since the construction of an RSEI requires a thermal infrared image, RSEI is mainly suitable for mesoscale mapping. The specific calculation method is shown in Table 2. Table 2. Details of the evaluation indicator system of sustainable development.

Index Quantitative Method Units Indicator Definition References
Environmental Reflecting the humidity of surface water, soil, and vegetation [37,38] Greenness (x 2 ) NDVI = (ρ NIR -ρ Red )/(ρ NIR + ρ Red ) % Reflecting the parameters of plant biomass, leaf area index, and vegetation coverage [39,40] Heat (x 3 ) Reflecting growth and distribution of vegetation, crop yield, and evaporation of surface water resources [33,41]  Urban-rural community (x 13 ) Fiscal expenditures on urban and rural community development yuan Support of urban-rural community development [52] Population density (x 14 ) Household registered population/area of town Capita/km 2 Reflecting the level of regional people agglomeration [53] Notes: ρ Blue , ρ Green , ρ Red , ρ NIR , ρ SWIR1 , and ρ SWIR2 respectively represent the reflectance of the blue, green, red, near-infrared (NIR), short wave infrared 1 (SWIR1), and short wave infrared 2 bands corresponding to the thematic mapper (TM) and the operational land unit (OLI) images; gain and skewness are the gain and offset values of TM6 band, or 0.056 and 1.238, respectively; K 1 and K 2 are scaling coefficients, obtained by image metadata, K 1 = 606.09 W/(m 2 × sr × µm), K 2 = 1282.71 K; center wavelength λ is 11.48 µm, α is 1.438 × 10 -2 mK; ε 6 is the surface specific emissivity based on TM6; while using the landsat OLI data, L 10 is the radiance value at the sensor, obtained by radiometric calibration; τ 10 is the transmittance of the atmosphere in the thermal infrared band; ε 10 is the surface specific emissivity; T s is the surface temperature; B 10 (T s ) is the same black body as T s radiation brightness; I 10 ↑ and I 10 ↓ are the radiance of the atmosphere up and down, respectively; the surface specific emissivity of the water bodies is 0.995; the surface specific emissivity ε 6 and ε 10 were obtained by Sobrino's NDVI threshold method; τ 10 , I 10 ↑, and I 10 ↓ parameters, with reference to the mid-latitude summer standard atmospheric profile, were obtained by interpolating the atmospheric profile according to the central latitude and longitude of the image and imaging time set; 1 RMB (yuan) ≈ 0.150 $.

Food Security and Socio-Economic Subsystem Evaluation Index
The indicators of regional FS and SE subsystems reflect different aspects; however, the selected indicators are representative, accessible, and comparable. Also, the indicators in line with characteristics of sustainable regional development refer to existing literature; the details are shown in Table 3. Table 3. Dynamic deviation state and characteristics of "eco-socio-economic" coordination.

Quadrant
Angle Range

Vector Dynamic Deviation State
Ecological Subsystem

Socio-Economic Subsystem
Food Security Subsystem Description vector that is higher (or lower) than the set target; "--" indicates it could not be determined if it meets the set target.
Multifunctional land use and economic development directly affect regional sustainability; the cropping system of the study region largely involves the intensive use of cultivated land. The grain output of cultivated land used for two to three crops annual will definitely be better than single cropped land or even land abandonment, so a multiple cropping index (x 5 ) was used reflecting the intensive level of land use for cultivated land resources in the study area. Although a strong correlation exists between grain production and the normalized difference vegetation index (NDVI), regional grain output can still be calculated more accurately using statistics provided by townships; we used the grain output per unit of cultivated land (x 6 ) to characterize the capacity of the study area for grain production. To ensure the sustainable development of agricultural production and the economy, China has passed legislation designed to implement special protection measures for prime farmland; the protection of prime farmland area (x 7 ) is designed to guarantee regional food security. Finally, the guarantee of food security is inseparable from government input. The financial expenditures in regional agricultural production (x 8 ), including agricultural, forestry, and water conservation expenditures, will comprehensively promote food output.
The social and economic development of townships and towns involve multiple sectors. During urbanization, the most important sectors include industrial development (x 9 ), per capita income (x 10 ), social security and employment (x 11 ), health (x 12 ), community development (x 13 ), and others. At the same time, population density (x 14 ) reflects the level of regional population agglomeration and is prominent when urbanization is ongoing. The level of industrialization in the town varies greatly due to natural endowments, policy support, and development opportunities; as a result, the growth rate of large-scale industries can better reflect the level of local economic growth and restructuring. The per capita income of residents can better reflect the comprehensive social and economic level of the region. These indicators are used to measure the achievements in social development and the protection of people's lives and work.

Principal Component Analysis and Standardization of Indicators
Principal component analysis (PCA) was used to integrate the indicators; the variance of the principal components was used as a weight to evaluate the various evaluation indicators. The specific steps started with Equation (1): where x i is the i th indicator, and I also needs to be standardized. The value range normalized by the Equation (1) is [0, 1]; the I closer is to 1, the better the sustainable development of the subsystem. Because a large value in the first principal component represents good environmental conditions, the initial ecological index I 0 was obtained using Equation (2): where PC 1 is the first principal component. Because the indicators that were obtained from remote sensing surveys, statistical yearbooks, and government websites were not uniform, the indicators were standardized based on the standard deviation, and the ranges were standardized to [0, 1], as shown in Equation (3).
where NI is the normalized index value, I is the numerical value of the index, and I max and I min are the maximum and minimum values of the indicator in the image, respectively.

Multi-Dimensional Evaluation Model
The evaluation of the coordination of regional environment changes and SE development, that is, the equilibrium state between various elements in a certain area, allowed the development of the following model: Assuming other conditions are homogeneous, the relationship between "eco-social-economic" is abstracted into a flute, the forces in the three directions in a Cartesian coordinate system. If all three forces reach the expected goal, the combined force is 0, that is, F is at the equilibrium point O, and the regional elements achieve sustainable development; otherwise, if one or more of them do not reach or perhaps exceed the established goal, the force makes F lean toward the equilibrium point. At this time, the regional development is in an unbalanced state. The concept is shown in Figure 3. In the present study, OA, OB, and OC represent the positive development relationships of the ecological, SE, and FS subsystem, respectively. When the combined force is 0, the angle between the corresponding angles of OA, OB, and OC is 2π/3; when F deviates from the equilibrium point O, F, the magnitude and orientation classes quantitatively characterize the corresponding deviation relationships. In order to facilitate the calculations, polar coordinates (F, θ) are used to characterize the sustainable regional development ordination. The polar diameter F represents the distance of any deviation. The larger the value of F, the poorer the coordination of the research unit; the polar angle θ represents the direction of the deviation. The direction angles of OA, OB, and OC were defined as π/2, 7π/6, and 11π/6, and the specific calculations for each variable are shown in Equations (4)- (6): where OA is the ecological degree of sustainability, X 1i is the ecological efficiency of the research unit, X 1S is the standard value of ecological efficiency, OB is the SE degree of sustainability, X 2i is the SE efficiency of the research unit, X 2S is the standard value of SE efficiency, OC is the FS degree of sustainability, X 3i is the FS efficiency of the research unit, and X 3S is the standard value of FS efficiency. All the standard values are the averages of the indicators. Based on the rules of vector operation, the polar coordinates (F) were obtained according to Equations (7)-(9): OC and → F are the vector expressions of ecological, economic, and social sustainability and the combined force of the three; OA, OB, and OC are the magnitudes of the forces where a positive or negative value indicates that its action direction is the same or opposite as the predetermined direction, respectively; x A , x B , x C and y A , y B , y C are the x-axis and y-axis coordinates. Based on the direction of the force determined in Figure 1, the resultant coordinate (x', y') was accumulated, and then the polar coordinate (F', θ) could be obtained. The extension lines OA', OB', and OC' were made in the opposite direction of OA, OB and OC. As the direction of the ecological, economic, and social reverse development vectors, the results were divided into six quadrants, and the characteristics of the dynamic state of each quadrant vector were formed according to the positive and negative relationship of the three elements of the "eco-socio-economic" vector. Each quadrant had positive and negative vectors used to define the spatial position relationship of each research unit. Through the quadrant of the result and its relatives to the vector direction, the degree of coordination and matching between the elements of each research unit could be determined. For more details, see Table 3.

Spatial Differential Results for RSEI
The RSEI was calculated using PCA in ArcGIS software before evaluating coordination; the indices of humidity, greenness, dryness, and heat were synthesized by band in each period after the standardization of the indices. Based on the spatial difference of the four indicators, RSEI fluctuated in different periods and different regions. In this study, the four indicators in each image were calculated separately, and then the average values were calculated. Finally, the index images in the two images were mosaicked together, and finally PCA was used. According to the calculation results of PCA, the contribution rate of the first principal component (FPC) of RSEI in 2018 was 84.29%, indicating that the FPC contained most of the variation in the characteristics of the four indicators. In the FPC, WET (abbreviation for humidity), and NDVI had positive values, indicating that they had a promoting effect on the quality of the environment, while LST and NDBSI had negative values, indicating that the impact of these two indicators on the environment was negative.
In addition, the coefficient affected by LST was the largest one. During the process of the data, to calculate the result corresponding to the quality of the environment, the first PC grid value of the region was used in Equation (1) and then standardized to obtain the RSEI index. The standardized RSEI and the four indicator values are shown in Figure 4a,b, respectively. .95% of all towns, respectively. The overall EQ of the Chaohu Lake Basin in 2018 was slightly lower. Based on a digital elevation model of the Chao Lake Basin ( Figure 5) and land use/cover types, the basin was divided into five regions (excluding the lake) to analyze the regional spatial changes in EQ. These regions included the Dabie Mountains in the southwest, the northwestern Hefei urban area, the northwestern hills, the farmland in the northeast, and the alluvial plain in the south. The findings suggest that the regional EQ has been greatly affected differences in precipitation and temperature. The highly erodible soil in the Dabie Mountains resulted in an ecological security pattern is characterized by being low in the north and high in the south; the southern area is one of the areas with a poor EQ in Chaohu Lake Basin.

Coordination of Sustainable Development
When we considered the coordination of the eco-socio-economic system of the Chaohu Lake Basin, the county scale was too large and the region was too heterogeneous for this type of analysis. Therefore, we chose township scale and urban built-up areas such as downtown Hefei as research units and used PCA to evaluate the sustainable development of FS and SE subsystems. The three dimensions of each system were calculated using Equations (4)-(6), and the results are shown in Figure 6a-c.

Analysis of the State of Sustainable Development in Independent Systems
(1) Degree of ecological sustainability In 2018, the areas with an average degree of ecological degree of sustainability where OA > 0 were mainly in the east of the Chaohu Lake Basin, which specifically includes towns with large areas of water or paddy fields. A total of 80 towns accounted for 57.36% of the entire basin (Figure 6a), of which the most prominent areas lie in the east along the Yangtze River; meanwhile, the EQ of 23 towns around the Dabie Mountains and Hefei City was lower than the average value. Judging only from the characteristics of the image, the environment in the northwest is negatively affected by the radiation of urban expansion.
(2) Degree of SE sustainability During the study period, the areas with a high OB (OB > 0) with the potential for sustainable socio-economic development were concentrated in downtown Hefei and in the counties of Chaohu, Hanshan, and Wuwei (including 61 towns; see Figure 6b). The sustainability of SE development in these well-developed regions was significantly higher than in other regions and was mainly dependent on the rapid development of their industries and large-scale service industries; this high OB is also inseparable from the government's support for these regions. The overall rate of sustainable socio-economic development of Jin'an District and Shucheng County in Lu'an City and the southern part of Lujiang County was low (OB < −0.5), in which 54 towns are involved. The SE development in these regions was relatively slow and therefore was less sustainable.

(3) Degree of FS sustainability
During the study period, the area of sustainable FS development (OC > 0) was mainly distributed in the northeast and southwest away from the city center where the terrain was relatively flat (Figure 6c). This area includes a total of 85 towns and accounts for 71.95% of the entire basin area. The Chaohu Lake Basin, a traditional agricultural production area, can still play a major role in the country's food security. However, 22 townships that are distributed in the downtown areas of Hefei and Chaohu as well as in Wuwei County still have the lowest sustainability of food security (OC < −0.5). Overall, FS sustainability is approximately sufficient for regional development.

Multidimensional System Coordination Analysis
(1) Coordination characteristics According to the geometric interval classification method, the value of F' can be divided into four categories, which are in a coordinated state (0 ≤ F' < 0.5), semi-coordinated state (0.5 ≤ F' < 1.0), antagonistic state (1.0 ≤ F' < 1.5), and disordered state (F' > 1.5). According to the calculations, the average F' value was 0.887, with no difference in most of each subsystem. It can be seen that the degree of coordination in support of sustainable development of the ecological, social, and economic systems of the Chaohu Lake Basin in 2018 was in a semi-coordinated state, and the environment, SE, and FS subsystems were still in the process of rapid development, meaning resource allocation needs to be further optimized.
(2) Spatial differences in coordination In Figure 6d, in terms of the spatial pattern of coordination, 42 towns were in a coordinated state, accounting for 23.95% of the total basin area; these towns were decentralized and were mainly located in Lujiang, Hanshan, and Feidong counties, and the coordination degree of sustainable development in these areas were at the highest level; ecological protection, SE development, and guaranteed balanced FS development and mutual promotion, and the comprehensive matching effect of regional resources was the best here.
A total of 69 towns were in the semi-coordinated state, accounting for 42.43% of the basin area. They are mostly distributed around the coordination region and are widely distributed and without obvious clusters; this type of area also needs to be integrated in the use of resources. Although it is inferior to the coordination type, this type of an area can achieve positive resource investment effects through the integration of resource elements.
Thirty towns were in the antagonistic state, accounting for 19.65% of the entire basin. These towns are evenly distributed in the borders of the counties with the exception of Jin'an District and Changfeng County, which are also close to the imbalanced area. Compared with the semi-coordinated type, the degree of coordination of factor allocation in this region is skewed, and ecological, SE, and FS have been over-or under-developed in one aspect or another, which leads to the mutual resistance of the factors and affects the sustainable development of the area. If the effective deployment of resources is not carried out appropriately, this will gradually transition to an imbalanced state.
Twenty towns were in a state of imbalance, accounting for 13.97% of the total area, which was the smallest for the research units of the present study; these areas were concentrated in downtowns, county centers, and in some towns along the river or reservoirs. In terms of sustainable development coordination, the resources have not been effectively matched, and most of this imbalanced region has experienced rapid social and economic development, while the sustainability of the EQ and FS are relatively low. Under the comprehensive effects, this is a current problem area for eco-social-economic comprehensive planning and sustainable development of Chaohu Lake Basin and should be the object of an urgently needed adjustment of resources.

(3) Skewness characteristics of sustainable development
Identifying the coordinated skewness characteristics of each subsystem is helpful to finding problem areas and providing an effectively remedy. Coordination skewness θ was distributed in quadrants I-VI, with the area of each quadrant accounting for 30.58%, 19.94%, 8.65%, 15.73%, 12.20%, and 12.91% of the basin, respectively. In general, there are the most dysfunctional quadrants (I, VI), followed by the transitional quadrants (II, V), and the least coordinated quadrants (III, IV). To some extent, this indicates that the development of the study area had an imbalance phenomenon in 2018, and resource allocation needs to be optimized.
Under the spatial pattern of sustainable regional development (Figure 6e), the resource allocation skewness of the dysfunctional quadrants is relatively serious, which reflects that only food safety is guaranteed, but environment security and SE development are lower than the currently established goals, which restricts the comprehensive and sustainable development of the region. Quadrant I mostly includes the traditional agricultural production areas of the Chaohu Lake Basin, such as Shucheng, Jin'an, and Changfeng that are mainly in the northernmost and southwest of the Chaohu Lake Basin. The coordination skewness of Quadrant I is expressed by having a better than average food supply and poor socioeconomic conditions; meanwhile, FS was high according to the set goals in Quadrant VI, although the environment of such areas was relatively poor. Such areas are mostly located in ecologically vulnerable areas located in hilly areas. These areas should be appropriately addressed to develop short-board subsystems in the future under the premise of protecting the environment.
The skewness of the degree of sustainable development coordination of the transitional quadrant stems from the inconsistency between the regional environment and SE development; that is, compared with the sustainable development goals in Quadrant II, the environment subsystem is more efficient and the matching SE development is extremely low; meanwhile, the coordination characteristics of the subsystems in Quadrant V are just the opposite. Quadrant II is widely distributed in the southern edge of the Chaohu Lake Basin and some towns in the northeast; Quadrant V is mainly located in Wuwei County or scattered around the cities. Ecological protection in these areas does not fully match the pace of rapid economic development.
The characteristics of the resource allocation and coordination skewness of the coordination quadrant are manifested because FS cannot be guaranteed, environmental protection or SE development has reached the currently established goals, and deviations in resource allocation are small. This quadrant is located in the center of some cities and counties, as well as in some strong industrial counties. It is represented by the urban areas of Hefei and Chaohu. In the future, the entire Chaohu Lake Basin should be developed radiantly with this region as the center.

Identification of Regional Sustainable Development Issues
Identifying problem areas is a prerequisite for formulating and improving related regional related. Based on comprehensive, coordinated measures and their spatial characteristics, identifying problem areas that need to be optimized in support of the sustainable development of the Chaohu Lake Basin will help to formulate site-specific regional control measures for planned development. This starts with the idea that the requirements of sustainable development include maintaining, rationally using, and improving the natural resource base, which supports ecological resistance to stress as well as SE growth. As far as regional coordination is concerned, through an analysis of the spatial pattern of the research results, six types of problem areas were identified in the research unit of the study area: (1) Imbalanced between high agriculture yield and inefficient SE (Quadrant I) The area of Quadrant I is mainly distributed in 45 towns in the main agricultural production areas, and this area has the largest number of towns of all types. The average value of F' in this area is 0.67, which is generally characterized by a high FS factor and slow SE development. This pattern is closely related to regional characteristics. As a major food producing area, the area has high land use rate and limited space for land use structure optimization. To coordinate sustainable regional development, it will be necessary to increase government investment in fixed assets, focus on SE development, and promote agriculture in the premise of ensuring future ecological and FS based on the development of agricultural auxiliary industries designed to boost the regional SE development.
(2) Imbalanced between inefficient SE and high ecological security (Quadrant II) The region of Quadrant II is relatively scattered, involving 36 towns, with an average F' of 0.69. In the Chaohu Lake Basin, except for the administrative area of Hefei City, although these SE areas are generally underdeveloped, they have a valuable and good EQ, so this pattern of inconsistency occurs. From the perspective of the overall sustainable development of the region, the region can still maintain a high-quality environment in the future. In view of the inefficient SE development, the income of residents in the region can be increased through government transfer payments and ecological compensation, and thereby better maintain a desired EQ.
(3) Imbalance between high ecological safety and low agricultural production (Quadrant III) Regardless of SE conditions of the region, the region of Quadrant III has a better environment and a lower FS includes a total of 17 towns with an average F' of 0.77. Except for Jin'an, Changfeng, and Feixi, a few towns in each county fall into this category. Because of a scarcity of cultivated land or the neglect of FS, this area lacks agriculture production. It is necessary to pay attention to FS in this region and to improve the land use rate and the efficiency of cultivated land use to maximize the food guarantee rate in the region in the next period. This will help the region to adapt to natural disasters and promote sustainable regional development. However, land use efficiency needs to be improved based on the premise of not harming the environment.
(4) Imbalance between low agriculture production and high SE efficiency (Quadrant IV) The low level of agricultural production, rapid SE development, and the EQ is not too bad in the area of Quadrant IV, which mainly includes the city or county centers for a total of 28 towns with an average F' of 1.87. The current level of urban and related industrial development in such areas means that almost no more cultivated land is available for reclamation, so it is not practical to increase agricultural output in these areas. A food reserve system needs to be established and food needs to be transferred from other regions to adjust the balance of grain supply and demand across the county and to guarantee FS. If natural disasters occur in the county, the ability to transfer food from the outside will also be required. Meanwhile, part of the financial support for the region should be used to finance ecological compensation funds to compensate local residents for ecological conservation in other regions.
(5) Imbalance between high SE efficiency and low ecological security (Quadrant V) The area of Quadrant V is characterized by economic development that has occurred at the expense of the environment; most of these areas are well-known industrial towns. The 17 towns have an average F' of 0.62, which shows that the degree of coordination of sustainable development in such regions is mostly in a semi-coordinated state. Then future sustainable regional development can be coordinated by optimizing resource allocation. Sustainable development will require strengthening the environment, which means putting ecological protection first, and focusing on developing green and pollution-free enterprises.
(6) Imbalance between low ecological security and high agriculture production (Quadrant VI) The area in Quadrant VI includes 18 towns with an average F' of 0.76. This area has good regional coordination and is also the main agricultural production area in the Chaohu Lake Basin. However, the environment is poor because of geomorphological undulations and a high potential for soil erosion. Management should focus on controlling soil erosion and identifying fragile areas to prevent ecological and anthropogenic damage in the future to maintain efficient crop production. This area is suitable for inclusion in a restoration plan to protect the cultivated land; after the fertility of cultivated land is restored, food security would be guaranteed.

Multifunctional Zoning in Support of Regional Sustainable Development
The county area is the most important and basic space unit in the administrative management system used for land in China. The county unit reflects the system's own integrity, stability, and regional multi-functionality. Current China's spatial planning of land use at the provincial-city-county-level had been implemented with the goal of forming "a county, a plan, and a blueprint." In this paper, zones divided based on land space planning functions divide the study area into priority development, agglomeration development, agricultural development, ecological conservation, and prohibited development zones (Figure 7). Based on regional industrial strengths, the development of SE systems is prioritized in the priority development and agglomeration development zones. The main function of the comprehensive agricultural development zone is to realize the sustainable development of regional agriculture and ensure FS; development is restricted in this zone. The ecological conservation zone refers to areas where ecological functions are important or ecosystems are fragile, and basic conditions for ecological and economic development are good. Prohibited development zones include the key ecological function areas and resource areas requiring special protection.
Because regional development is multi-functional, regional linkages and compound functions should be considered during sustainability analysis. A comprehensive evaluation system used to analyze the multifunctional suitability of land use and its functional zoning should meet the requirements of a land spatial planning system. Based on the diagnosis of regional problems, social, ecological, and economic development should be integrated to support the three subsystems of SE, FS, and ES, to allow the integrated and sustainable development of a system run in an orderly way.

Conclusions
This paper analyzes the core indicators of the three subsystems including ecological security, SE development, and food security. First, the regional EQ was calculated so that an eco-socio-economic model could be constructed based on a mechanical equilibrium model. The study used a multidimensional evaluation model and takes the township as the research unit in 2018. Then, the coordination degree of ecological security, SE development, and food security in the Chaohu Lake Basin was comprehensively evaluated. The main conclusions are as follows: (1) The environment of the Chaohu Lake Basin is mainly affected by factors such as climate, urbanization, and topography. The EQ of the entire study area is slightly inferior to that of other areas, and a slight trend toward deterioration was identified. The spatial differences across the landscape are obvious. The urban centers radiate into the surrounding areas so that the environment in those area tends to be degraded, while the environment in other areas is expected to improve over time in the future. Ecological conservation should continue to be emphasized in the future in support of sustainable regional development. (2) The coordination of ecological security, SE development, and food security in the Chaohu Lake Basin is generally in the process of development and run-in, especially in towns and city surrounding areas. However, many towns are still in a coordinated state, meaning sustainable and coordinated development can still be achieved on a regional scale through resource optimization. The area was divided into six patterns according to the degree of coordination deviation. Most of the areas are mainly rich in grain production but poor in economic development conditions. It can be seen that the overall SE conditions in the Chaohu Lake Basin need to be improved while ensuring that the environment is not degraded. (3) Based on the requirements of the main functional zoning of a county's territorial and spatial planning, the 161 towns in the entire Chaohu Lake Basin were divided into five functional zones, of which most are considered to be comprehensive agricultural zones; after all, the basin is a traditional agricultural area. The second most common zones are ecological conservation areas in which development is restricted; these areas have good environmental conditions. These regions continue to be maintained, although the protection of ecologically fragile areas needs to be strengthened. In the future, ecological compensation and food security system development should be carried out in urban priority and agglomeration development areas to pave the way for sustainable regional development.
This study provides monitoring of important nature ecosystems and serves as a guide to future sustainable regional development by providing directions and paths to follow, it can be used as a basis for spatial planning and also as a reference for eco-functional zoning that is applicable to various natural regions worldwide. In the future, this research will continue to expand to explore 30-years of change using historical data so as to better study sustainable development.