<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Environ. Sci.</journal-id>
<journal-title>Frontiers in Environmental Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Environ. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-665X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1076344</article-id>
<article-id pub-id-type="doi">10.3389/fenvs.2023.1076344</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Environmental Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Assessment of ecosystem health based on landscape pattern in ecologically fragile regions at different spatial scales: A case study of Dianchi Lake basin, China</article-title>
<alt-title alt-title-type="left-running-head">Lei et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenvs.2023.1076344">10.3389/fenvs.2023.1076344</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Lei</surname>
<given-names>Dongmei</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2061890/overview"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhang</surname>
<given-names>Yuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Ran</surname>
<given-names>Yuju</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Gao</surname>
<given-names>Liping</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Jie</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Li</surname>
<given-names>Zhiguo</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Mo</surname>
<given-names>Jinxiao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Xuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>College of Finance and Public Administration</institution>, <institution>Yunnan University of Finance and Economics</institution>, <addr-line>Kunming</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Geography and Tourism</institution>, <institution>Zhaotong University</institution>, <addr-line>Zhaotong</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/89366/overview">Shuisen Chen</ext-link>, Guangzhou Institute of Geography, China</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1993494/overview">Xiang Zhang</ext-link>, Wuhan University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1569182/overview">Chao Chen</ext-link>, Zhejiang Ocean University, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Dongmei Lei, <email>dmLei@ynufe.edu.cn</email>
</corresp>
<fn fn-type="other">
<p>This article was submitted to Environmental Informatics and Remote Sensing, a section of the journal Frontiers in Environmental Science</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>23</day>
<month>01</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1076344</elocation-id>
<history>
<date date-type="received">
<day>21</day>
<month>10</month>
<year>2022</year>
</date>
<date date-type="accepted">
<day>12</day>
<month>01</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Lei, Zhang, Ran, Gao, Li, Li, Mo and Liu.</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Lei, Zhang, Ran, Gao, Li, Li, Mo and Liu</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>In the process of urbanization, maintaining ecosystem health in ecologically fragile regions is an inevitable requirement for sustainable development in the future, but the ecological health status varies greatly at different spatial scales. Determining the appropriate spatial scale of ecological health assessment is the basis and premise of ecological management and environmental protection policy making. Taking Dianchi Lake basin as the research area, based on the vigor-organization-resilience (VOR) model and landscape pattern index, this study established an assessment system of ecological health to evaluate the ecosystem health status from quadrat, ecological and administrative scales, with the aim to determine a reasonable evaluation scale and strengthen regional sustainability. The results show that: 1) The suitable spatial grain of landscape pattern analysis in Dianchi Lake basin is 150&#xa0;m. 2) The average values of ecosystem health indices at quadrat scale, ecological scale and administrative scale are 0.5466, 0.4977 and 0.5373, respectively, demonstrating a sub-health state. 3) The most suitable scale for ecological health assessment of Dianchi Lake basin is ecological scale, because taking sub-watershed as ecological unit can ensure the integrity of the evaluated ecosystem to the greatest extent. Multiscale assessment is helpful for eco-environmental management departments to understand the ecosystem health status at different scales and provide a scientific basis for regional eco-environmental management decisions.</p>
</abstract>
<kwd-group>
<kwd>region ecosystem health</kwd>
<kwd>multiscale assessment</kwd>
<kwd>VOR model</kwd>
<kwd>landscape pattern</kwd>
<kwd>Dianchi Lake basin</kwd>
</kwd-group>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Ecosystem health refers to the ability of an ecosystem to keep its organizational structure, respond to changing conditions, and recover from stressors (<xref ref-type="bibr" rid="B29">Rapport, 1989</xref>; <xref ref-type="bibr" rid="B6">Das et al., 2020</xref>). Watershed, as a special regional ecosystem and a unique geographical unit, is crucial for managing runoff, conserving water, purifying water quality, and preserving regional biological balance. A healthy watershed is very important to natural ecosystems and social systems (<xref ref-type="bibr" rid="B20">Luo et al., 2003</xref>; <xref ref-type="bibr" rid="B39">Wu et al., 2015</xref>). It can not only provide different ecological products or services for human beings, but also has certain self-adjustment and recovery ability to natural and human disturbance (<xref ref-type="bibr" rid="B39">Wu et al., 2015</xref>). Unfortunately, human activities directly or indirectly affect all aspects of the watershed, resulting to the continuous deterioration of the watershed ecosystem (<xref ref-type="bibr" rid="B3">Chen et al., 2016</xref>). In particular, watersheds in ecologically fragile regions are facing mounting problems of ecological protection and economic development. The effects of urbanization process, particularly the wasteful development, over-exploitation and utilization of water and soil resources, will result in large-scale changes in the watershed environment, and the structure and function of the watershed ecosystem, which poses a serious threat to the ecosystem health (<xref ref-type="bibr" rid="B52">Zhao et al., 2012</xref>; <xref ref-type="bibr" rid="B47">Zhang H. et al., 2013a</xref>). Therefore, accurately evaluating ecological health of watershed is critical for promoting regional sustainable development in ecologically fragile region, both theoretically and practically.</p>
<p>As a new goal of environmental management, ecosystem health has become a hot spot and a new direction of environmental research (<xref ref-type="bibr" rid="B10">Guo et al., 2002</xref>; <xref ref-type="bibr" rid="B41">Xiao and Ou, 2002</xref>; <xref ref-type="bibr" rid="B7">Fan et al., 2006</xref>; <xref ref-type="bibr" rid="B27">Peng et al., 2007</xref>; <xref ref-type="bibr" rid="B50">Zhang W. et al., 2020b</xref>; <xref ref-type="bibr" rid="B4">Cheng et al., 2020</xref>; <xref ref-type="bibr" rid="B6">Das et al., 2020</xref>; <xref ref-type="bibr" rid="B26">Pan et al., 2020</xref>; <xref ref-type="bibr" rid="B33">Shen et al., 2021</xref>). <xref ref-type="bibr" rid="B32">Schaeffer et al. (1988)</xref> first discussed the measurement of ecosystem health in 1988, and proposed principles and methods to evaluate ecosystem health. Since the 1980s, a great deal of research has been done on watershed ecosystem health, involving the method of biological integrity index evaluation (IBI) (<xref ref-type="bibr" rid="B1">Baek et al., 2014</xref>; <xref ref-type="bibr" rid="B18">Lin et al., 2021</xref>), planktonic index of biotic integrity (PIBI) (<xref ref-type="bibr" rid="B51">Zhang Y. et al., 2020c</xref>), landscape pattern index (<xref ref-type="bibr" rid="B8">Fang et al., 2012</xref>; <xref ref-type="bibr" rid="B42">Xie et al., 2015</xref>), press-state-response (PSR) model (<xref ref-type="bibr" rid="B11">He et al., 2015</xref>), fuzzy mathematical model (<xref ref-type="bibr" rid="B42">Xie et al., 2015</xref>), vitality-organization-resilience (VOR) model and so on (<xref ref-type="bibr" rid="B21">Ma et al., 2001</xref>; <xref ref-type="bibr" rid="B20">Luo et al., 2003</xref>; <xref ref-type="bibr" rid="B16">Lei et al., 2020</xref>). In particular, VOR model is a commonly used comprehensive evaluation method for ecosystem health, which contains three factors: ecosystem vitality, ecosystem organization and ecosystem resilience (<xref ref-type="bibr" rid="B5">Costanza et al., 1997</xref>; <xref ref-type="bibr" rid="B28">Rapport et al., 1998</xref>; <xref ref-type="bibr" rid="B46">Yuan et al., 2019</xref>; <xref ref-type="bibr" rid="B44">Yan et al., 2020</xref>). On the other hand, the process of watershed urbanization has significantly changed land use and landscape pattern (<xref ref-type="bibr" rid="B31">Robert et al., 1989</xref>). Accordingly, it is a trend that current scholars use landscape patterns to evaluate ecosystem health at the regional scale. Landscape pattern-based ecological health evaluation is generally carried out by constructing a comprehensive ecological health evaluation index system and using an analytic hierarchy process (<xref ref-type="bibr" rid="B25">Ou et al., 2018</xref>). As an example, <xref ref-type="bibr" rid="B8">Fang et al. (2012)</xref> calculated the ecological health index of the Luanhe River basin by using the landscape pattern index, and analyzed the ecosystem health status of various cities. <xref ref-type="bibr" rid="B42">Xie et al. (2015)</xref> clarified the spatial distribution and limiting factors of watershed ecosystem health by using the landscape pattern index and other evaluation indicators in the Tiaoxi watershed. Existing research provides an established research method and system for landscape pattern-based ecological health assessment. However, there are still some obvious shortcomings in the existing research. Firstly, the diversity of different ecosystems makes it difficult to establish a unified indicator system to evaluate all ecosystems (<xref ref-type="bibr" rid="B44">Yan et al., 2020</xref>). Secondly, ecological health assessment based on landscape pattern index is scale-dependent, but most studies fail to clarify the appropriate spatial granularity of landscape index calculation (<xref ref-type="bibr" rid="B35">Su et al., 2020</xref>). Additionally, in the regional ecosystem health assessment, a single-scale assessment of city, province (city, county) region, watershed, and other ecosystems is mainly used, while a multiscale assessment is rare. In terms of the study area, the assessment of watershed ecosystem health in the ecologically fragile regions is still insufficient.</p>
<p>In this study, taking the typical ecologically fragile region&#x2014;Dianchi Lake basin as the study area, we used the VOR model to determine the optimal size of a landscape pattern index to evaluate the distribution of ecosystem health index variables at three spatial levels: quadrat, ecological, and administrative scale. The aim was to deepen the understanding of ecosystem health assessment by examining the following related issues: 1) How does the landscape pattern index alter when grain size changes? 2) Are there any differences in ecosystem health status at different scales in the study area using the VOR model based on landscape pattern index? 3) Which scale is the most suitable scale for watershed ecosystem health assessment? Thus, it provides a scientific basis for making regional resource and environment management decisions and sustainable development strategies in the plateau lake basin.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Material and methods</title>
<sec id="s2-1">
<title>2.1 Research case</title>
<p>Dianchi Lake basin, located in Yunnan-Guizhou Plateau in western China, is a typical plateau basin with remarkable ecological sensitivity and vulnerability (<xref ref-type="bibr" rid="B47">Zhang H. et al., 2013a</xref>). Specifically, it has a subtropical monsoon humid climate with an annual temperature of 15&#xb0;C and an average annual rainfall of 935&#xa0;mm. The whole Dianchi Lake basin is within the administrative area of Kunming city, covering an area of about 2909.41&#xa0;km<sup>2</sup>. It is made up of nine county-level administrative districts, including the main metropolitan region of Kunming: Wuhua, Guandu, Panlong, Xishan, and Chenggong (<xref ref-type="bibr" rid="B52">Zhao et al., 2012</xref>; <xref ref-type="bibr" rid="B47">Zhang H. et al., 2013a</xref>) (<xref ref-type="fig" rid="F1">Figure 1</xref>). As the economic, political and cultural core of Yunnan Province, the Dianchi Lake basin only accounts for 14% of Kunming&#x2019;s land area, but it has 52% of this city&#x2019;s population and 70% of its economic aggregate (<xref ref-type="bibr" rid="B17">Li et al., 2012</xref>; <xref ref-type="bibr" rid="B36">Wang et al., 2020</xref>). The economic development of the whole basin is closely related to the economic development of Kunming city. The land use types of Dianchi Lake basin are complex and diverse, including farmland, forest land, grassland, water and built-up land, of which 27% are farmland, 26% are built-up land, and 11% are water and wetland (<xref ref-type="table" rid="T1">Table 1</xref>).</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Map showing the geographical location and land use of Dianchi Lake basin in 2018.</p>
</caption>
<graphic xlink:href="fenvs-11-1076344-g001.tif"/>
</fig>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Area of different categories of land use in Dianchi Lake basin in 2018.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Land use types</th>
<th align="center">Farmland</th>
<th align="center">Forest land</th>
<th align="center">Grassland</th>
<th align="center">Built-up land</th>
<th align="center">Water</th>
<th align="center">Unused land</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Land-use area (km<sup>2</sup>)</td>
<td align="center">785.26</td>
<td align="center">922.77</td>
<td align="center">59.00</td>
<td align="center">768.17</td>
<td align="center">306.61</td>
<td align="center">67.60</td>
</tr>
<tr>
<td align="center">Area proportion (%)</td>
<td align="center">27</td>
<td align="center">32</td>
<td align="center">2</td>
<td align="center">26</td>
<td align="center">11</td>
<td align="center">2</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>However, the development intensity of Dianchi Lake basin is more than 15%, and available land and water resources are limited and scarce, with low environmental capacity (<xref ref-type="bibr" rid="B36">Wang et al., 2020</xref>). As the region with the highest degree of urbanization in Yunnan Province, it cannot ignore the realistic damage and potential threat caused by urbanization to watershed ecosystem health (<xref ref-type="bibr" rid="B43">Xu, 2011</xref>). Since the 1980s, rapid urbanization of Kunming city has caused serious disturbance and damage to the natural hydrological cycle of Dianchi Lake basin, mainly in the following aspects: the degradation of river ecosystem function, water pollution, especially the increasingly serious eutrophication of Dianchi Lake water, the reduction of ecological land area, serious soil and water loss, and the sharp decrease of the density and curvature of major rivers (<xref ref-type="bibr" rid="B15">Lei and Xu., 2018</xref>). At present, the urbanization process has a very prominent impact on the Dianchi Lake basin, especially the water environment deterioration and eutrophication in Dianchi Lake that has threatened the ecosystem health status of the whole basin and restricted the sustainable development of the lakeside social economy (<xref ref-type="bibr" rid="B43">Xu, 2011</xref>; <xref ref-type="bibr" rid="B42">Xie et al., 2015</xref>). Therefore, Dianchi Lake basin has highly urgent and representative in regional development, outstanding ecological vulnerability, rapid urbanization, and dramatic land-use/land-cover change, which can offer practical guiding significance for multiscale ecosystem health assessment (<xref ref-type="bibr" rid="B15">Lei and Xu., 2018</xref>).</p>
</sec>
<sec id="s2-2">
<title>2.2 Data source</title>
<p>The basic data of this study include the following. 1) Land use data. We used Landsat TM/OLI remote sensing images (2018) with a spatial resolution of 30&#xa0;m. The images were taken on cloudless days from February to March. On this basis, according to the national standard Classification of Land Use Status, the landscape types of the study area were divided into six categories: forest land, farmland, water, grassland, built-up land and unused land. These classification results were further corrected by manual visual interpretation using background data from field investigations. Finally, it ensures that the Kappa coefficient of classification results is greater than 85%, which meets the application accuracy requirements. 2) DEM data with resolution of 30&#xa0;m. All these data were downloaded from the geospatial data cloud (<ext-link ext-link-type="uri" xlink:href="http://www.gscloud.cn/">www.gscloud.cn</ext-link>). All spatial data were processed and mapped using ENVI (version 5.1) and ArcGIS (version 10.2) software.</p>
</sec>
<sec id="s2-3">
<title>2.3 Methods</title>
<p>The purpose of this study was to determine the appropriate scale for ecosystem health assessment in Dianchi Lake basin based on landscape pattern index and VOR model. Three steps are included in this study. The first was to ensure the appropriate grain size for the landscape pattern index. The second was to evaluate ecosystem health at different scales based on VOR model. The third was to compare and identify the most appropriate evaluation scale based on the results of multiscale assessment. The specific technology roadmap was shown in <xref ref-type="fig" rid="F2">Figure 2</xref>.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Technology roadmap of this study.</p>
</caption>
<graphic xlink:href="fenvs-11-1076344-g002.tif"/>
</fig>
<sec id="s2-3-1">
<title>2.3.1 Selection of landscape metrics and suitable space granularity</title>
<p>The landscape pattern is the spatial arrangement and structural composition of landscape elements in a region (<xref ref-type="bibr" rid="B40">Xiao et al., 2003</xref>). Landscape index can effectively and comprehensively reflect landscape pattern, such as landscape structure, landscape diversity, spatial heterogeneity and fragmentation (<xref ref-type="bibr" rid="B37">Wu, 2004</xref>; <xref ref-type="bibr" rid="B24">Mairota et al., 2013</xref>; <xref ref-type="bibr" rid="B2">Chefaoui, 2014</xref>). Therefore, in this study, based on previous studies (<xref ref-type="bibr" rid="B48">Zhang L. L. et al., 2013b</xref>; <xref ref-type="bibr" rid="B49">Zhang Q. et al., 2020</xref>; <xref ref-type="bibr" rid="B12">Hou et al., 2020</xref>), nine landscape indices (LPIs) were selected to analyze the landscape pattern and the grain change of LPIs in Dianchi Lake basin (<xref ref-type="table" rid="T2">Table 2</xref>). Based on the land use/cover data, FRAGSTATS 4.2 software was used to calculate selected landscape indices.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Landscape metrics used in this study.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Landscape metrics</th>
<th align="left">Abbreviation</th>
<th align="left">Explanation</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Mean Patch Shape Index</td>
<td align="left">SHAPE-MN</td>
<td align="left">Mean patch shape index can directly represent the complexity of shape</td>
</tr>
<tr>
<td align="left">Mean Patch Fractal Dimension</td>
<td align="left">FRAC-MN</td>
<td align="left">This indicator is based on the fractal dimension of each patch, which is related to the circumference and area of the patch. The higher the FRAC-MN value, the more complex the patch shape</td>
</tr>
<tr>
<td align="left">Shape Coefficient of Variation</td>
<td align="left">SHAPE-CV</td>
<td align="left">The index of shape regularity of each patch in landscape</td>
</tr>
<tr>
<td align="left">Patch Richness</td>
<td align="left">PR</td>
<td align="left">Number of landscape types</td>
</tr>
<tr>
<td align="left">Shannon&#x2019;s Evenness Index</td>
<td align="left">SHEI</td>
<td align="left">SHEI values range from 0 to 1. When SHEI value is small, dominance degree is generally high, which can reflect that landscape is dominated by one or a few dominant patch types. When SHEI approaches 1, the dominance degree is low, indicating that there are no obvious dominant patches in the landscape, and each patch type is evenly distributed in the landscape</td>
</tr>
<tr>
<td align="left">Landscape Shape Index</td>
<td align="left">LSI</td>
<td align="left">The LSI value is greater than or equal to 1 and has an upper limit. LSI reflects the shape changes of the landscape. A larger value indicates that the shape is more complex and irregular</td>
</tr>
<tr>
<td align="left">Contagion Index</td>
<td align="left">CONTAG</td>
<td align="left">The CONTAG value ranges from 0 to 100. This indicator describes the degree of agglomeration or extension of different patch types in a landscape. Theoretically, a small CONTAG value indicates that there are many small patches in the landscape. When it approaches 100, it indicates that there are dominant patch types with high connectivity in the landscape</td>
</tr>
<tr>
<td align="left">Patch Cohesion Index</td>
<td align="left">COHESION</td>
<td align="left">The physical connectedness of the focal patch type. COHESION approaches 0 as the proportion of the landscape comprised of the focal class decreases and becomes increasingly subdivided and less physically connected</td>
</tr>
<tr>
<td align="left">Aggregation Index</td>
<td align="left">AI</td>
<td align="left">The degree of aggregation of cells of the focal class</td>
</tr>
<tr>
<td align="left">Shannon&#x2019;s Diversity Index</td>
<td align="left">SHDI</td>
<td align="left">It is an index of complexity and variability of various patch types in landscape</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Spatial granularity is the area, volume and length of features represented by the smallest identifiable units in the landscape (<xref ref-type="bibr" rid="B35">Su et al., 2020</xref>). In landscape ecology, the study of ecological health based on landscape pattern is scale dependent (<xref ref-type="bibr" rid="B38">Wu et al., 2002</xref>). When the spatial granularity of landscape analysis changes, the characteristics of ecosystem will change accordingly. Therefore, the selection of appropriate spatial granularity is the premise of regional ecosystem health assessment and can ensure the accuracy of the results (<xref ref-type="bibr" rid="B22">Ma et al., 2018</xref>; <xref ref-type="bibr" rid="B35">Su et al., 2020</xref>). This study adopted the traditional spatial granularity selection method, that was, the changing trend of multiple landscape pattern indices under different spatial granularity is analyzed to select the appropriate spatial granularity (<xref ref-type="bibr" rid="B48">Zhang L. L. et al., 2013b</xref>; <xref ref-type="bibr" rid="B13">Ji et al., 2013</xref>; <xref ref-type="bibr" rid="B22">Ma et al., 2018</xref>). First, the landscape types of Dianchi Lake basin were divided into grid images by albers coordinate projection method. Second, we used the nearest method in ArcGIS 10.2 to re-sample at 90, 150, 210, 270, 330, and 390&#xa0;m for assessing granularity.</p>
</sec>
<sec id="s2-3-2">
<title>2.3.2 The division of multiscale space</title>
<p>Referring to the existing relevant studies (<xref ref-type="bibr" rid="B22">Ma et al., 2018</xref>), this study divided the regional ecosystem into three different spatial scales, which were quadrat scale, ecological scale, and administrative scale. Each spatial scale is composed of different numbers of basic units. The units of these three scales are quadrat unit, ecological unit, and administrative unit, respectively.<list list-type="simple">
<list-item>
<p>(1) Quadrat scale. The study area was divided into different plots by setting different size quadrat and the ecosystem health index of different plots was calculated. In this scale, ecosystem health calculation is based on a single quadrat, so the size of quadrat must affect the spatial distribution of ecosystem health. In order to analyze and verify this problem, this study divided different sizes of quadrat units and analyzed the trend of ecosystem health changes at different quadrat sizes. The quadrat sizes were divided into 5&#xa0;km &#xd7; 5&#xa0;km, 10&#xa0;km &#xd7; 10&#xa0;km, 15&#xa0;km &#xd7; 15&#xa0;km, and 20&#xa0;km &#xd7; 20&#xa0;km respectively. And the most suitable quadrat size for ecosystem health assessment of Dianchi Lake basin was selected by comparison and evaluation.</p>
</list-item>
<list-item>
<p>(2) Ecological scale. The integrity of the structure and process of natural elements in the evaluation unit can be guaranteed by the boundary of natural surfaces such as watersheds (<xref ref-type="bibr" rid="B27">Peng et al., 2007</xref>). This paper used the hydrological analysis function of ArcGIS10.2 software to extract the sub-watersheds with relatively complete structure and ecological elements through the catchment area of the basin, which are regarded as the evaluation ecological units to assess and characterize the ecosystem health. The designed flow rate of this study was &#x3e;25000. Finally, 51 ecological units, namely sub-watersheds, were obtained for ecosystem health evaluation of Dianchi Lake basin on the ecological scale.</p>
</list-item>
<list-item>
<p>(3) Administrative scale. According to the administrative boundary of the region, the study area was divided into administrative units, and the ecosystem health values of each administrative unit were calculated and spatially characterized as a whole. The evaluation on the administrative scale is helpful to understand the situation of each administrative region and is conducive to decision-makers to formulate a variety of policies to deal with ecological risks (<xref ref-type="bibr" rid="B27">Peng et al., 2007</xref>). The study region was divided into eight administrative units based on the municipal districts and counties under the administration of Kunming city within the scope of the study.</p>
</list-item>
</list>
</p>
</sec>
<sec id="s2-3-3">
<title>2.3.3 Multiscale assessment of ecosystem health</title>
<sec id="s2-3-3-1">
<title>2.3.3.1 VOR model</title>
<p>In this study, the VOR model was used to evaluate ecosystem health at quadrat, administrative and ecological scales. The VOR model includes three factors: ecosystem vitality, ecosystem organization and ecosystem resilience (<xref ref-type="bibr" rid="B5">Costanza et al., 1997</xref>; <xref ref-type="bibr" rid="B46">Yuan et al., 2019</xref>). In detail, vitality can be defined as the power and nutrient cycling capacity of an ecosystem. Organization refers to the diversity of ecosystem composition and approaches, which is the foundation of vitality, resilience and service function. Resilience represents the ability of an ecosystem to keep the stability of its structure and pattern or restore its original health status under external disturbance (<xref ref-type="bibr" rid="B8">Fang et al., 2012</xref>). Based on the VOR model, the calculation <xref ref-type="disp-formula" rid="e1">formula (1)</xref> of ecosystem health index is as follows:<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>H</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf1">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>H</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem health index of the evaluation unit <italic>i</italic>; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>V</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem vitality index of the evaluation unit <italic>i,</italic> equal to the value of NDVI of the evaluation unit <italic>i</italic>; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem organization index of the evaluation unit <italic>i;</italic> <inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem resilience index of the evaluation unit <italic>i.</italic>
</p>
<p>
<inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is calculated by the following <xref ref-type="disp-formula" rid="e2">formula (2)</xref>:<disp-formula id="e2">
<mml:math id="m7">
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>Where <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem organization index of the evaluation unit <italic>i;</italic> n represents the number of indicators representing <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>O</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>; <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:msubsup>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the standardized value of <inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>;</mml:mo>
</mml:mrow>
</mml:math>
</inline-formula> <inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the weight of the indicator <italic>i</italic>.</p>
<p>
<inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is calculated by the following <xref ref-type="disp-formula" rid="e3">formula (3)</xref>:<disp-formula id="e3">
<mml:math id="m14">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2211;</mml:mo>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>Where <inline-formula id="inf12">
<mml:math id="m15">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the ecosystem resilience index of the evaluation unit <italic>i,</italic> namely comprehensive resilience index; <inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the area of land use type <italic>i</italic>; <inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:msubsup>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the resilience value of land use type <italic>i</italic>; <inline-formula id="inf15">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is Shannon&#x2019;s Diversity index of the evaluation unit <italic>i.</italic> <xref ref-type="table" rid="T3">Table 3</xref> showed the resilience value of six land use types in the Dianchi Lake basin.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Resilience value of land use types in Dianchi Lake basin.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Land use types</th>
<th align="center">Farmland</th>
<th align="center">Forest land</th>
<th align="center">Grassland</th>
<th align="center">Built-up land</th>
<th align="center">Water</th>
<th align="center">Unused land</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Resilience value</td>
<td align="center">0.5</td>
<td align="center">0.9</td>
<td align="center">0.7</td>
<td align="center">0.4</td>
<td align="center">0.9</td>
<td align="center">0.0</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-3-3-2">
<title>2.3.3.2 Determination of assessment indices and index weight</title>
<p>To make the selected index system reflect the health status of watershed ecosystem more completely and accurately, there are three principles for selecting indicators from the landscape pattern index. 1) Rationality. The selected indexes have integrity and spatial scale suitability. That is, informative and comprehensive indicators, which can explain the interrelation and influence within the basin ecosystem; The index should also consider the spatial size of the evaluation object, which is suitable for the spatial scale of the watershed evaluation. 2) Accessibility. The selected indicators are simple and operable, and have strong comparability. The parameters are easy to obtain, easy to calculate and analyze. 3) Representativeness. It can be identified and reflect the essential characteristics of Dianchi Lake basin.</p>
<p>Following the above principles, and combining with the existed literature (<xref ref-type="bibr" rid="B8">Fang et al., 2012</xref>; <xref ref-type="bibr" rid="B44">Yan et al., 2020</xref>), we set and screened the indicators of ecosystem health evaluation based on the VOR model and the landscape pattern index in Dianchi Lake basin, and the results were shown in <xref ref-type="table" rid="T4">Table 4</xref>. Specifically, the vitality index can be represented by the first potential productivity of organisms. The organization index can select Shannon&#x2019;s diversity index, Shannon&#x2019;s evenness index, Contagion index, and Interspersion and Juxtaposition index from two aspects of landscape level and patch type level. The resilience index can be represented by the comprehensive resilience index (<xref ref-type="bibr" rid="B5">Costanza et al., 1997</xref>; <xref ref-type="bibr" rid="B46">Yuan et al., 2019</xref>) (<xref ref-type="table" rid="T4">Table 4</xref>).</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Ecosystem health evaluation index and its weights in Dianchi Lake basin.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Target layer</th>
<th align="center">Element layer</th>
<th align="center">Weight</th>
<th align="center">Indicator layer</th>
<th align="center">Weight</th>
<th align="center">Property of indicator</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="7" align="center">Basin ecosystem health index</td>
<td align="center">Vitality</td>
<td align="center">0.2014</td>
<td align="center">Normalized difference vegetation index (NDVI)</td>
<td align="center">0.2014</td>
<td align="center">Positive</td>
</tr>
<tr>
<td rowspan="5" align="center">Organization</td>
<td rowspan="5" align="center">0.6796</td>
<td align="center">Forest Coverage rate</td>
<td align="center">0.0840</td>
<td align="center">Positive</td>
</tr>
<tr>
<td align="center">SHDI</td>
<td align="center">0.2208</td>
<td align="center">Positive</td>
</tr>
<tr>
<td align="center">COHESION</td>
<td align="center">0.2244</td>
<td align="center">Positive</td>
</tr>
<tr>
<td align="center">CONTAG</td>
<td align="center">0.0694</td>
<td align="center">Positive</td>
</tr>
<tr>
<td align="center">Interspersion and Juxtaposition index (IJI)</td>
<td align="center">0.0810</td>
<td align="center">Positive</td>
</tr>
<tr>
<td align="center">Resilience</td>
<td align="center">0.1190</td>
<td align="center">Comprehensive resilience index</td>
<td align="center">0.1190</td>
<td align="center">Positive</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The weight of each index was calculated using the entropy approach. <xref ref-type="table" rid="T4">Table 4</xref> showed the weighted results of the indices.</p>
</sec>
<sec id="s2-3-3-3">
<title>2.3.3.3 Standardization of assessment indices</title>
<p>Because the units of each evaluation index are different, it is necessary to standardize the evaluation index. We use <xref ref-type="disp-formula" rid="e4">formula (4)</xref> to standardize the data and eliminate the influence of dimension, amplitude and positive directions. In this way, all the index values will be ranged in [0,1]. The <xref ref-type="disp-formula" rid="e4">formula (4)</xref> is as follows:<disp-formula id="e4">
<mml:math id="m19">
<mml:mrow>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">s</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">v</mml:mi>
<mml:mi mathvariant="normal">e</mml:mi>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">n</mml:mi>
<mml:mi mathvariant="normal">d</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:mi mathvariant="normal">t</mml:mi>
<mml:mi mathvariant="normal">o</mml:mi>
<mml:mi mathvariant="normal">r</mml:mi>
<mml:mo>:</mml:mo>
<mml:mtext>&#x2009;</mml:mtext>
<mml:msubsup>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">X</mml:mi>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
</p>
<p>Where <inline-formula id="inf16">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> denotes the value of index <italic>i</italic>; <inline-formula id="inf17">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="italic">max</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf18">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi mathvariant="italic">min</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is the maximum and minimum value of index <italic>i</italic>, respectively; <inline-formula id="inf19">
<mml:math id="m23">
<mml:mrow>
<mml:msubsup>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> denotes the standardized value of indicator <italic>i</italic>.</p>
</sec>
<sec id="s2-3-3-4">
<title>2.3.3.4 Comprehensive evaluation of ecosystem health</title>
<p>Ecosystem health assessment is an indicator measuring the degree of ecological health, and its value is between 0 and 1. The standardized index value and index weight are brought into the comprehensive evaluation model (<xref ref-type="disp-formula" rid="e5">formula (5)</xref>).<disp-formula id="e5">
<mml:math id="m24">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:munderover>
<mml:mstyle displaystyle="true">
<mml:mo>&#x2211;</mml:mo>
</mml:mstyle>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msup>
<mml:msub>
<mml:mi>X</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2032;</mml:mo>
</mml:msup>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
</p>
<p>In the formula, <inline-formula id="inf20">
<mml:math id="m25">
<mml:mrow>
<mml:mi>E</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> is comprehensive level of ecosystem health; <inline-formula id="inf21">
<mml:math id="m26">
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> is weight of the indicator <italic>i</italic>; <inline-formula id="inf22">
<mml:math id="m27">
<mml:mrow>
<mml:msubsup>
<mml:mrow>
<mml:mtext>&#x2009;</mml:mtext>
<mml:mi>X</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2032;</mml:mo>
</mml:msubsup>
</mml:mrow>
</mml:math>
</inline-formula> is the standardized value of indicator <italic>i</italic>; n represents the number of indicators.</p>
</sec>
<sec id="s2-3-3-5">
<title>2.3.3.5 Classification standard of ecosystem health</title>
<p>At present, there is no unified classification standard of ecosystem health. Based on the results of ecosystem health research (<xref ref-type="bibr" rid="B34">Sheng et al., 2011</xref>) and the current situation of the Dianchi Lake basin, the comprehensive ecosystem health index was divided into 5 classes from high to low by using the equal-discontinuous method. Each level was assigned with different ecosystem health index value range and ecosystem characteristics. Specifically, the value ranges are level &#x2160; superiorly healthy (&#x3e;0.8), level &#x2161; healthy (0.6&#x2013;0.8), level &#x2162; sub-healthy (0.4&#x2013;0.6), level &#x2163; unhealthy (0.2&#x2013;0.4), and level &#x2164; disease (0.2&#x2013;0). The classification standard of watershed ecosystem health status is shown in <xref ref-type="table" rid="T5">Table 5</xref>.</p>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Classification of comprehensive ecosystem health.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Classification standard</th>
<th align="left">Health level</th>
<th align="left">Comprehensive index value</th>
<th align="left">Characteristics of watershed ecosystem</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Superiorly healthy</td>
<td align="left">&#x2160;</td>
<td align="left">1.0&#x2013;0.8</td>
<td align="left">The landscape of the watershed remains in good condition, and the ecosystem functions are very perfect with strong systemic vitality. The ecological structure is very reasonable and stable, with strong natural resilience and little influence from human activities</td>
</tr>
<tr>
<td align="left">Healthy</td>
<td align="left">&#x2161;</td>
<td align="left">0.8&#x2013;0.6</td>
<td align="left">The landscape of the watershed remains in the natural state, and the ecosystem functions are perfect and has strong system vitality. The ecological structure is reasonable and the natural resilience is strong, which is less affected by human activities</td>
</tr>
<tr>
<td align="left">Sub-healthy</td>
<td align="left">&#x2162;</td>
<td align="left">0.6&#x2013;0.4</td>
<td align="left">The landscape of the watershed is changed to some extent, and the ecosystem plays its basic ecological function and has certain system vitality. The ecological structure is complete, the natural resilience is normal, and the ecosystem is greatly affected by human activities</td>
</tr>
<tr>
<td align="left">Unhealthy</td>
<td align="left">&#x2163;</td>
<td align="left">0.4&#x2013;0.2</td>
<td align="left">The landscape of the watershed is destroyed; the ecosystem structure is missing and incomplete; the vitality is low; the ecosystem is unstable</td>
</tr>
<tr>
<td align="left">Disease</td>
<td align="left">&#x2164;</td>
<td align="left">0.2&#x2013;0</td>
<td align="left">The landscape of the watershed is destroyed completely; the ecosystem structure is seriously missing and incomplete; the vitality is very low; the ecosystem is seriously deteriorated</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Selection of suitable space granularity</title>
<p>As the grain size increased, the values of landscape pattern indices changed and showed different trends, which indicated that landscape pattern indices were affected by the change of spatial grain size to a certain extent (<xref ref-type="table" rid="T6">Table 6</xref>). Compared with SHAPE-CV and LSI, other landscape pattern indices varied more slightly. SHAPE-CV and LSI changed obviously and fluctuated greatly, indicating the existence of particle size effect. The reason was that with the increase of spatial grain size, smaller plaques merged into peripheral plaques, which led to the rapid reduction of small plaques, and the reduction of small plaques led to the regularity of plaque shape, which was manifested as the decrease of LSI and the fluctuation of SHAPE-CV. When the spatial grain size was 150&#xa0;m, the two landscape pattern indices, SHAPE-CV and LSI, fluctuated obviously. Therefore, 150&#xa0;m was chosen as the appropriate spatial grain size for this study.</p>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Change of landscape pattern indices at different spatial grains in Dianchi Lake basin.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Space granularity(m)</th>
<th align="center">SHAPE-MN</th>
<th align="center">FRAC-MN</th>
<th align="center">SHAPE-CV</th>
<th align="center">PR</th>
<th align="center">SHEI</th>
<th align="center">LSI</th>
<th align="center">Contag</th>
<th align="center">Cohesion</th>
<th align="center">AI</th>
<th align="center">SHDI</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">90</td>
<td align="center">1.23</td>
<td align="center">1.03</td>
<td align="center">64.72</td>
<td align="center">6.00</td>
<td align="center">0.82</td>
<td align="center">77.00</td>
<td align="center">37.33</td>
<td align="center">98.67</td>
<td align="center">74.98</td>
<td align="center">1.47</td>
</tr>
<tr>
<td align="center">150</td>
<td align="center">1.25</td>
<td align="center">1.03</td>
<td align="center">67.90</td>
<td align="center">6.00</td>
<td align="center">0.80</td>
<td align="center">47.82</td>
<td align="center">38.38</td>
<td align="center">98.14</td>
<td align="center">74.51</td>
<td align="center">1.44</td>
</tr>
<tr>
<td align="center">210</td>
<td align="center">1.26</td>
<td align="center">1.03</td>
<td align="center">67.93</td>
<td align="center">6.00</td>
<td align="center">0.79</td>
<td align="center">35.03</td>
<td align="center">39.00</td>
<td align="center">97.78</td>
<td align="center">74.20</td>
<td align="center">1.41</td>
</tr>
<tr>
<td align="center">270</td>
<td align="center">1.28</td>
<td align="center">1.03</td>
<td align="center">68.73</td>
<td align="center">6.00</td>
<td align="center">0.78</td>
<td align="center">28.27</td>
<td align="center">39.32</td>
<td align="center">97.27</td>
<td align="center">73.62</td>
<td align="center">1.40</td>
</tr>
<tr>
<td align="center">330</td>
<td align="center">1.29</td>
<td align="center">1.03</td>
<td align="center">65.26</td>
<td align="center">6.00</td>
<td align="center">0.77</td>
<td align="center">23.57</td>
<td align="center">39.53</td>
<td align="center">96.88</td>
<td align="center">73.42</td>
<td align="center">1.38</td>
</tr>
<tr>
<td align="center">390</td>
<td align="center">1.31</td>
<td align="center">1.03</td>
<td align="center">65.41</td>
<td align="center">6.00</td>
<td align="center">0.77</td>
<td align="center">20.39</td>
<td align="center">39.89</td>
<td align="center">96.12</td>
<td align="center">73.22</td>
<td align="center">1.37</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2">
<title>3.2 Characteristics of ecosystem health at different scales</title>
<sec id="s3-2-1">
<title>3.2.1 Quadrat scale</title>
<p>At four quadrat scales (5&#xa0;km &#xd7; 5&#xa0;km, 10&#xa0;km &#xd7; 10&#xa0;km, 15&#xa0;km &#xd7; 15&#xa0;km, and 20&#xa0;km &#xd7; 20&#xa0;km), the ecosystem health values of the study region were spatially represented (<xref ref-type="fig" rid="F3">Figure 3</xref>). The results of ecosystem health distribution at four quadrat scales showed that with the increase of quadrat units, the overall patch shape tended to be simple, the loss of ecosystem health information was serious, and the precision of ecosystem health representation decreased. Because the characteristics of ecosystem health at the 10&#xa0;km &#xd7; 10&#xa0;km quadrat scale were more detailed, and the integrity of ecosystem was higher than that at the 5&#xa0;km &#xd7; 5&#xa0;km quadrat scale and other quadrat scales, the spatial differentiation analysis was carried out at this scale (<xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Ecosystem health in Dianchi Lake basin at four quadrat scales.</p>
</caption>
<graphic xlink:href="fenvs-11-1076344-g003.tif"/>
</fig>
<p>At the quadrat scale (10&#xa0;km &#xd7; 10&#xa0;km), the areas of healthy, sub-healthy and unhealthy were 1308.50&#xa0;km<sup>2</sup>, 1499.21&#xa0;km<sup>2</sup> and 99.38&#xa0;km<sup>2</sup>, accounting for 45.01%, 51.57% and 3.42% of the total area of Dianchi Lake basin, respectively (<xref ref-type="fig" rid="F3">Figure 3</xref>). The average value of ecosystem health was 0.5466 on this scale, which indicated that it was in the sub-health state. In the study area, the unhealthy regions were mainly distributed in the central part of the Dianchi Lake basin, that was, Dianchi Lake and lakeside areas with high urbanization level. While the sub-healthy and healthy regions were mainly located in the mountainous regions with high forest coverage in the north and south. From the perspective of the whole basin, the spatial distribution of health index showed that the value of health index gradually increased from the Dianchi Lake as the center to the surrounding lakeside and mountain areas.</p>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Ecological scale</title>
<p>At the ecological scale (51 sub-watersheds), the areas of healthy, sub-healthy and unhealthy were 873.90&#xa0;km<sup>2</sup>, 1 843.35&#xa0;km<sup>2</sup> and 192.38&#xa0;km<sup>2</sup>, accounting for 30.03%, 63.35% and 6.62% of the total area of Dianchi Lake basin, respectively (<xref ref-type="fig" rid="F4">Figure 4A</xref>). The average value of ecosystem health was 0.4977 on this scale, indicating that it was in the sub-health state. As shown in <xref ref-type="fig" rid="F4">Figure 4A</xref>, the healthy areas in the study area were mainly distributed in the northern and southern edge areas of the Dianchi Lake basin, which were also areas with low urbanization level and low human disturbance degree. The sub-healthy areas were concentrated in the middle of the basin. The unhealthy area formed four landscape patches, two of which were distributed in the north of Dianchi Lake basin that was the main urban area of Kunming city.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Ecosystem health in Dianchi Lake basin at ecological scale and administrative scale.</p>
</caption>
<graphic xlink:href="fenvs-11-1076344-g004.tif"/>
</fig>
</sec>
<sec id="s3-2-3">
<title>3.2.3 Administrative scale</title>
<p>At the administrative scale (8 units of administrative district), the areas of healthy, sub-healthy and unhealthy were 1102.49&#xa0;km<sup>2</sup>, 1 772.62&#xa0;km<sup>2</sup>, and 8.81&#xa0;km<sup>2</sup>, accounting for 38.23%, 61.47%, and 0.30% of the total area of Dianchi Lake basin, respectively (<xref ref-type="fig" rid="F4">Figure 4B</xref>). The average value of ecosystem health was 0.5373 on this scale, indicating that it was in the sub-health state. As shown in <xref ref-type="fig" rid="F4">Figure 4B</xref>, the healthy areas in the study area were mainly distributed in two administrative regions in the south of Dianchi Lake region&#x2014;Jinning District with low urbanization level and Xishan District with high forest coverage rate. The sub-healthy areas were mainly concentrated in Wuhua District, Chenggong District, Guandu District, Panlong District and Songming District with high urbanization level in the middle and north of Dianchi Lake basin. The unhealthy area formed a landscape patch, which was located in Fuming District on the northern margin of Dianchi Lake basin. Overall, the distribution characteristics of health index showed that the ecosystem health status of the southern region of the Dianchi Lake basin was better than that of the northern region.</p>
</sec>
</sec>
<sec id="s3-3">
<title>3.3 Comparative analysis of multiscale ecosystem health assessment</title>
<p>The results of multiscale ecosystem health assessment are shown in <xref ref-type="table" rid="T7">Table 7</xref>. There were differences in ecosystem health index values, health level and area at three different scales in Dianchi Lake basin, and the order of ecosystem health index was as follows: quadrat scale &#x3e; administrative scale &#x3e; ecological scale (<xref ref-type="table" rid="T7">Table 7</xref>).</p>
<table-wrap id="T7" position="float">
<label>TABLE 7</label>
<caption>
<p>Ecosystem health area and its proportion at different scales.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="left">Scale</th>
<th align="left">Health level</th>
<th align="left">Area (km<sup>2</sup>)</th>
<th align="left">Area ratio (%)</th>
<th align="left">Mean of overall health status</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="5" align="left">Quadrat scale</td>
<td align="left">Superiorly healthy</td>
<td align="left">0</td>
<td align="left">0</td>
<td rowspan="5" align="center">0.5466</td>
</tr>
<tr>
<td align="left">Healthy</td>
<td align="left">1308.50</td>
<td align="left">40.01</td>
</tr>
<tr>
<td align="left">Sub-healthy</td>
<td align="left">1499.21</td>
<td align="left">51.57</td>
</tr>
<tr>
<td align="left">Unhealthy</td>
<td align="left">99.38</td>
<td align="left">3.42</td>
</tr>
<tr>
<td align="left">Disease</td>
<td align="left">0</td>
<td align="left">0</td>
</tr>
<tr>
<td rowspan="5" align="left">Ecological scale</td>
<td align="left">Superiorly healthy</td>
<td align="left">0</td>
<td align="left">0</td>
<td rowspan="5" align="center">0.4977</td>
</tr>
<tr>
<td align="left">Healthy</td>
<td align="left">873.90</td>
<td align="left">30.03</td>
</tr>
<tr>
<td align="left">Sub-healthy</td>
<td align="left">1843.35</td>
<td align="left">63.35</td>
</tr>
<tr>
<td align="left">Unhealthy</td>
<td align="left">192.38</td>
<td align="left">6.62</td>
</tr>
<tr>
<td align="left">Disease</td>
<td align="left">0</td>
<td align="left">0</td>
</tr>
<tr>
<td rowspan="5" align="left">Administrative scale</td>
<td align="left">Superiorly healthy</td>
<td align="left">0</td>
<td align="left">0</td>
<td rowspan="5" align="center">0.5373</td>
</tr>
<tr>
<td align="left">Healthy</td>
<td align="left">1102.69</td>
<td align="left">38.23</td>
</tr>
<tr>
<td align="left">Sub-healthy</td>
<td align="left">1772.62</td>
<td align="left">61.47</td>
</tr>
<tr>
<td align="left">Unhealthy</td>
<td align="left">8.81</td>
<td align="left">0.30</td>
</tr>
<tr>
<td align="left">Disease</td>
<td align="left">0</td>
<td align="left">0</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The overlay analysis of ecosystem health results at three scales (<xref ref-type="table" rid="T8">Table 8</xref>) showed that: 1) For the quadrat scale and administrative scale, there were two overlapping regions of healthy and sub-healthy grades, and the overlap area of sub-healthy accounted for the largest proportion (29.61%). 2) For the quadrat scale and ecological scale, there were three overlapping regions of healthy, sub-healthy and unhealthy grades, and the overlap area of sub-healthy accounted for the largest proportion (38.26%). 3) For the administrative scale and ecological scale, there were two overlapping regions of healthy and sub-healthy grades, and the overlap area of sub-healthy accounted for the largest proportion (43.03%). 4) For the three scales, there were two overlapping regions of healthy and sub-healthy grades, with the total overlapped area accounting for 34.59%.</p>
<table-wrap id="T8" position="float">
<label>TABLE 8</label>
<caption>
<p>Proportion of overlay area of ecosystem health at different scales (%).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Overlay of different scales</th>
<th colspan="5" align="center">Health level</th>
<th rowspan="2" align="left">Total</th>
</tr>
<tr>
<th align="left">Superiorly healthy</th>
<th align="left">Healthy</th>
<th align="left">Sub-healthy</th>
<th align="left">Unhealthy</th>
<th align="left">Disease</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Quadrat scale &#xd7; Ecological scale</td>
<td align="left">0</td>
<td align="left">20.39</td>
<td align="left">38.26</td>
<td align="left">1.51</td>
<td align="left">0</td>
<td align="left">60.16</td>
</tr>
<tr>
<td align="left">Administrative scale &#xd7; Quadrat scale</td>
<td align="left">0</td>
<td align="left">16.42</td>
<td align="left">29.61</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">46.03</td>
</tr>
<tr>
<td align="left">Administrative scale &#xd7; Ecological scale</td>
<td align="left">0</td>
<td align="left">15.13</td>
<td align="left">43.03</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">58.16</td>
</tr>
<tr>
<td align="left">Administrative scale &#xd7; Ecological scale &#xd7; Quadrat scale</td>
<td align="left">0</td>
<td align="left">10.67</td>
<td align="left">23.92</td>
<td align="left">0</td>
<td align="left">0</td>
<td align="left">34.59</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>
<xref ref-type="table" rid="T8">Table 8</xref> further showed that, the overlap area of ecological health regions at the ecological scale and quadrat scale was the largest based on the overlapping results of ecological health regions at all levels, which was 60.16%. The mean values of the overall ecosystem health of these two scales were in the range of 0.4&#x2013;0.6, with the most similar health evaluation results. Although it was easy to control the accuracy of ecological health analysis by dividing ecological health regions with quadrat scale, its disadvantage was that it could not ensure the ecological integrity of the study region, especially for the plateau watershed with a high degree of terrain fragmentation. On the contrary, taking the sub-watershed as the evaluation unit could ensure the ecosystem integrity of the study area to the greatest extent. Therefore, it was more reasonable to use the ecological scale than the quadrat scale to analyze the ecosystem health in this study area.</p>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<sec id="s4-1">
<title>4.1 Evaluation index system based on landscape pattern index</title>
<p>Establishing an adequate index system is the basis and key to evaluating ecosystem health objectively and accurately. However, there is no universally adopted index system in ecosystem health assessment (<xref ref-type="bibr" rid="B45">Yang et al., 2017</xref>). The main reasons are as follows. Firstly, the establishment of the index system involves many subjects, and it is necessary to have enough understanding of the evaluation system. Secondly, different indicators should be selected to evaluate ecosystem health at different scales for different types of ecosystems. Therefore, it is impossible to generalize the index system for evaluating ecosystem health (<xref ref-type="bibr" rid="B9">Gao et al., 2010</xref>). Comprehensively considering the actual situation of land use in the Dianchi Lake basin, this study used the landscape pattern theory of landscape ecology to understand the ecosystem health of the basin and constructed the ecosystem health evaluation index system consisting of ecosystem vitality, organization and resilience. The evaluation results of the three scales coincided well, which verified the scientific and operational nature of this evaluation system. The evaluation index system constructed in this study enriches the existing technical method system of ecosystem health evaluation.</p>
<p>Some studies also believe that the landscape pattern index quantitatively describes the quantity and pattern of land cover and the morphological structure of surface vegetation, which can reflect the impact of human activities on watershed ecosystem health and is an appropriate indicator for watershed ecosystem health evaluation (<xref ref-type="bibr" rid="B40">Xiao et al., 2003</xref>). Because the regional ecosystem is composed of multiple landscapes, the ecosystem service function of the region is maintained dynamically by landscape structure and function, and the different adjacent relationships will inevitably affect the service function. Accordingly, landscape pattern plays an important role in ecosystem health, and landscape pattern index can be used as a suitable indicator for regional ecosystem health assessment (<xref ref-type="bibr" rid="B30">Ren et al., 2000</xref>; <xref ref-type="bibr" rid="B8">Fang et al., 2012</xref>). In addition, some studies also consider other indicators for ecosystem health assessment (<xref ref-type="bibr" rid="B9">Gao et al., 2010</xref>). For example, <xref ref-type="bibr" rid="B9">Gao et al. (2010)</xref> chose the landscape pattern index but also selected social and economic statistics data and ecological and environmental factors data from surface monitoring to construct an indicator system of ecosystem health evaluation.</p>
</sec>
<sec id="s4-2">
<title>4.2 The suitable scale of ecosystem health assessment</title>
<p>There are three kinds of evaluation units for ecosystem health assessment. Firstly, the spatial grid is taken as the evaluation unit (<xref ref-type="bibr" rid="B14">Jiang et al., 2021</xref>). Its advantage is that regional spatial differences can be expressed in more detail, while its disadvantage is that the spatial resolutions of multi-source data are mostly different, and data resampling operation is required to obtain the same spatial resolution (<xref ref-type="bibr" rid="B14">Jiang et al., 2021</xref>). Secondly, as an evaluation unit, the administrative unit is more common and widely-used (<xref ref-type="bibr" rid="B19">Lu et al., 2010</xref>). Since the monitoring projects and data statistics of China&#x2019;s social economy, environment, forestry, agriculture, water conservancy and other departments are conducted based on administrative units, the selection of administrative areas for evaluation units will facilitate data acquisition, ensure the integrity of the evaluated ecosystem, and promote the application of evaluation results (<xref ref-type="bibr" rid="B9">Gao et al., 2010</xref>). Lastly, the sub-watershed is taken as the evaluation unit (<xref ref-type="bibr" rid="B39">Wu et al., 2015</xref>; <xref ref-type="bibr" rid="B23">Ma et al., 2019</xref>). The sub-watershed is a complete geomorphic unit with ecological integrity and can reduce the interference of some overly detailed information (<xref ref-type="bibr" rid="B42">Xie et al., 2015</xref>).</p>
<p>The results of this study indicated that the most suitable scale for evaluating ecosystem health in the plateau lake basin was the ecological scale, that is, the sub-watershed was the evaluation unit. This is consistent with the research of some scholars (<xref ref-type="bibr" rid="B20">Luo et al., 2003</xref>). For example, <xref ref-type="bibr" rid="B20">Luo et al. (2003)</xref> also believed that the study of watershed ecosystem health should take into special consideration the drainage system of the watershed. Due to the continuous characteristics of rivers, the health damage of tributaries often affects the health of the mainstream and even the whole watershed. Therefore, the continuity and integrity of tributaries should be fully considered when selecting the scale of watershed ecosystem health research (<xref ref-type="bibr" rid="B20">Luo et al., 2003</xref>).</p>
</sec>
<sec id="s4-3">
<title>4.3 Applicability of VOR model</title>
<p>The results of this study showed that using the VOR model is a good way to evaluating the ecosystem health status of the Dianchi Lake basin: at the quadrat, administrative and ecological scales, the spatial distribution of the comprehensive evaluation results of ecosystem health are relatively consistent, and the overall evaluation results are relatively similar. However, in this study, the ecosystem health evaluation index system based on the VOR model highlighted the characteristics of the natural ecosystem of the watershed and ignored the urbanization pressure and social and economic development pressure of the watershed as a typical socio-economic-natural complex ecosystem.</p>
<p>In the future, to further study the ecological health of watershed, the main exploration contents are: 1) To build a new interdisciplinary and multiscale health research framework, collect data by &#x201c;3S&#x201d; technology, and update the health indicators timely; 2) The evaluation units are classified according to the main driving forces affecting the ecological process of sub-watershed units, and then the evaluation index system is set up by zoning classification. Additionally, the connotation, evaluation process and application of watershed ecosystem health have not been unified, and need to be further analyzed and discussed.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>Based on the landscape pattern index, this study constructed a VOR model to evaluate the ecosystem health of the Dianchi Lake basin from three spatial scales of quadrat, ecological and administrative, and determined the appropriate spatial scale for ecological health evaluation. The conclusions were as follows.<list list-type="simple">
<list-item>
<p>(1) According to the change characteristics of landscape pattern indices with grain size change, the appropriate spatial grain size to calculate landscape pattern indices in Dianchi Lake basin is 150&#xa0;m.</p>
</list-item>
<list-item>
<p>(2) The ecosystem health status at quadrat scale, ecological scale and administrative scale are all sub-health, but there are differences in ecosystem health index and healthy area proportion. There are overlapping areas among these three scales, but the overlap area between the quadrat scale and the ecological scale is the largest, accounting for 60.16%.</p>
</list-item>
<list-item>
<p>(3) The health distribution results at the quadrat scale and the ecological scale had the highest similarity, and the sub-watershed as the evaluation unit could guarantee the integrity of the ecosystem to the greatest extent. Therefore, the ecological scale is the most suitable scale to evaluate the ecosystem health in the Dianchi Lake basin.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>DL: Conceptualization, Methodology, Formal analysis, Writing, Supervision, Project administration, Funding acquisition. YZ: Investigation, Original draft, Data curation, Visualization. YR: Data curation. LG: Data curation. JL: Review and editing. ZL: Review. JM: Resources. XL: Data curation.</p>
</sec>
<sec id="s8">
<title>Funding</title>
<p>This research was supported by the Humanities and Social Sciences Foundation of the Ministry of Education (Grant No. 20XJAZH005), and the Basic Research Project of Science and Technology Department of Yunnan Province (Grant No. 202001BB050073).</p>
</sec>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Baek</surname>
<given-names>S. H.</given-names>
</name>
<name>
<surname>Son</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Choi</surname>
<given-names>H. W.</given-names>
</name>
<name>
<surname>Kim</surname>
<given-names>Y. O.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Assessing the ecosystem health status of Korea Gwangyang and Jinhae bays based on a planktonic index of biotic integrity (P-IBI)</article-title>. <source>Ocean Sci. J.</source> <volume>49</volume> (<issue>3</issue>), <fpage>291</fpage>&#x2013;<lpage>311</lpage>. <pub-id pub-id-type="doi">10.1007/s12601-014-0029-2</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chefaoui</surname>
<given-names>R. M.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Landscape metrics as indicators of coastal morphology: A multi-scale approach</article-title>. <source>Ecol. Indic.</source> <volume>45</volume> (<issue>5</issue>), <fpage>139</fpage>&#x2013;<lpage>147</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2014.04.004</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>S. Z.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Z. F.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>L. H.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Influence of human intervention on ecosystem health in the Karst rock desertification controlling process</article-title>. <source>Res. Soil Water Conservation</source> <volume>23</volume> (<issue>2</issue>), <fpage>213</fpage>&#x2013;<lpage>219</lpage>.</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheng</surname>
<given-names>W. J.</given-names>
</name>
<name>
<surname>Xi</surname>
<given-names>H. Y.</given-names>
</name>
<name>
<surname>Sindikubwabo</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Si</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>C. G.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>T. F.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Ecosystem health assessment of desert nature reserve with entropy weight and fuzzy mathematics methods: A case study of badain jaran desert</article-title>. <source>Ecol. Indic.</source> <volume>119</volume>, <fpage>106843</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.106843</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Costanza</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>D&#x2019;Arge</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Groot</surname>
<given-names>R. D.</given-names>
</name>
<name>
<surname>Farber</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Grasso</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hannon</surname>
<given-names>B.</given-names>
</name>
<etal/>
</person-group> (<year>1997</year>). <article-title>The value of the world&#x2019;s ecosystem services and natural capital</article-title>. <source>Nature</source> <volume>25</volume> (<issue>1</issue>), <fpage>253</fpage>&#x2013;<lpage>260</lpage>. <pub-id pub-id-type="doi">10.1038/387253a0</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Das</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Das</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Mandal</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Research note: Ecosystem Health (EH) assessment of a rapidly urbanizing metropolitan city region of eastern India &#x2013; a study on Kolkata Metropolitan Area</article-title>. <source>Landsc. Urban Plan.</source> <volume>204</volume>, <fpage>103938</fpage>. <pub-id pub-id-type="doi">10.1016/j.landurbplan.2020.103938</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>R. L.</given-names>
</name>
<name>
<surname>Su</surname>
<given-names>W. C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z. J.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>A preliminary discussion on the influence factor and the assessment method of the ecosystem health</article-title>. <source>Res. Soil Water Conservation</source> <volume>13</volume> (<issue>6</issue>), <fpage>82</fpage>&#x2013;<lpage>86</lpage>.</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Z. C.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Shan</surname>
<given-names>C. J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Regional ecosystem health assessment based on landscape pattern&#x2014;a case study for the administrative region in the Luanhe River basin</article-title>. <source>South-to-North Water Divers. Water Sci. Technol.</source> <volume>10</volume> (<issue>6</issue>), <fpage>37</fpage>&#x2013;<lpage>41</lpage>.</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>Z. G.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Weng</surname>
<given-names>Y. B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Y. Q.</given-names>
</name>
<name>
<surname>D</surname>
<given-names>Y. F.</given-names>
</name>
<etal/>
</person-group> (<year>2010</year>). <article-title>Comprehensive assessment of ecosystem health on multiple scales: A case study for metropolitan area of ningbo</article-title>. <source>Acta Ecol. Sin.</source> <volume>30</volume> (<issue>7</issue>), <fpage>1706</fpage>&#x2013;<lpage>1717</lpage>.</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Guo</surname>
<given-names>X. R.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>X. Q.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Primary studies on urban ecosystem health assessment</article-title>. <source>China Environ. Sci.</source> <volume>22</volume> (<issue>6</issue>), <fpage>46</fpage>&#x2013;<lpage>50</lpage>.</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>He</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>G. H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>R. J.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>W. Q.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Temporal and spatial variation of land ecosystem health based on the pressure-state-response model: A case study of pinggu district, beijing</article-title>. <source>J. Nat. Resour.</source> <volume>30</volume> (<issue>12</issue>), <fpage>2057</fpage>&#x2013;<lpage>2068</lpage>.</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Hou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>F. Q.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>X. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The spatial characteristics and relationships between landscape pattern and ecosystem service value along an urban-rural gradient in Xi&#x27;an city, China</article-title>. <source>Ecol. Indic.</source> <volume>108</volume>, <fpage>105720</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2019.105720</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ji</surname>
<given-names>Y. Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X. L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H. B.</given-names>
</name>
</person-group> (<year>2013</year>). <article-title>Research on the differences of landscape patterns under several spatial granularities in Zhenjiang city</article-title>. <source>China Land Sci.</source> <volume>27</volume> (<issue>5</issue>), <fpage>54</fpage>&#x2013;<lpage>62</lpage>.</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>D. F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Ecological space zoning based on lake ecosystem health maintenance: Taking the hongfeng Lake basin as example</article-title>. <source>Acta Ecol. Sin.</source> <volume>41</volume> (<issue>7</issue>), <fpage>2571</fpage>&#x2013;<lpage>2581</lpage>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lei</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>X. Y.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Study on indicator system of ecosystem health assessment in Dianchi Lake basin on background of urbanization</article-title>. <source>Resour. Dev. Mark.</source> <volume>34</volume> (<issue>7</issue>), <fpage>902</fpage>&#x2013;<lpage>906</lpage>.</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lei</surname>
<given-names>J. R.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z. Z.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y. Q.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>X. H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y. L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>T. T.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Development and application of wetland ecosystem health evaluation system in Hainan province</article-title>. <source>Wetl. Sci.</source> <volume>18</volume> (<issue>5</issue>), <fpage>555</fpage>&#x2013;<lpage>563</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z. J.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D. W.</given-names>
</name>
<name>
<surname>Ni</surname>
<given-names>J. B.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Impacts of 20-year socio-economic development on aquatic environment of Lake Dianchi Basin</article-title>. <source>J. Lake Sci.</source> <volume>24</volume> (<issue>6</issue>), <fpage>875</fpage>&#x2013;<lpage>882</lpage>. <pub-id pub-id-type="doi">10.18307/2012.0610</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>F. F.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Fang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>T. Y.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>W. Z.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Ecological health assessments of rivers with multiple dams based on the biological integrity of phytoplankton: A case study of north creek of jiulong river</article-title>. <source>Ecol. Indic.</source> <volume>121</volume>, <fpage>106998</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.106998</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>L. Z.</given-names>
</name>
<name>
<surname>Zhan</surname>
<given-names>Y. Z.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>Y. M.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>J. Y.</given-names>
</name>
<name>
<surname>Mou</surname>
<given-names>Y. M.</given-names>
</name>
</person-group> (<year>2010</year>). <article-title>Regional ecosystem health assessment based on land use pattern: A case of study of zhoushan island</article-title>. <source>Acta Ecol. Sin.</source> <volume>30</volume> (<issue>1</issue>), <fpage>245</fpage>&#x2013;<lpage>252</lpage>.</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>Y. C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Z. X.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Deng</surname>
<given-names>H. B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Assessment methods of watershed ecosystem health</article-title>. <source>Acta Ecol. Sin.</source> <volume>23</volume> (<issue>8</issue>), <fpage>1606</fpage>&#x2013;<lpage>1614</lpage>.</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>K. M.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>H. M.</given-names>
</name>
<name>
<surname>Guan</surname>
<given-names>Q. B.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>B. J.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Ecosystem health assessment: Methods and directions</article-title>. <source>Acta Ecol. Sin.</source> <volume>21</volume> (<issue>12</issue>), <fpage>2106</fpage>&#x2013;<lpage>2116</lpage>.</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>X. Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Duan</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Multi-scale landscape ecological risk assessment in ecologically fragile regions: A case study in gaoqu town in mizhi county, shaanxi province</article-title>. <source>Chin. J. Ecol.</source> <volume>37</volume> (<issue>10</issue>), <fpage>3171</fpage>&#x2013;<lpage>3178</lpage>.</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>Z. G.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L. L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>R. T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>B. J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Assessment and countermeasures of ecosystem health in the Beijing-Zhangjiakou area</article-title>. <source>Arid Zone Res.</source> <volume>36</volume> (<issue>1</issue>), <fpage>228</fpage>&#x2013;<lpage>236</lpage>.</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mairota</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Cafarelli</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Boccaccio</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Leronni</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Labadessa</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Kosmidou</surname>
<given-names>V.</given-names>
</name>
<etal/>
</person-group> (<year>2013</year>). <article-title>Using landscape structure to develop quantitative baselines for protected area monitoring</article-title>. <source>Ecol. Indic.</source> <volume>33</volume> (<issue>10</issue>), <fpage>82</fpage>&#x2013;<lpage>95</lpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2012.08.017</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ou</surname>
<given-names>W. X.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L. J.</given-names>
</name>
<name>
<surname>Tao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>A land-cover-based approach to assessing the spatio-temporal dynamics of ecosystem health in the Yangtze River Delta region</article-title>. <source>China Popul. Resour. Environ.</source> <volume>28</volume> (<issue>5</issue>), <fpage>84</fpage>&#x2013;<lpage>92</lpage>.</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pan</surname>
<given-names>Z. Z.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>J. H.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>D. F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J. W.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Predicting the joint effects of future climate and land use change on ecosystem health in the middle reaches of the Yangtze River Economic Belt, China</article-title>. <source>Appl. Geogr.</source> <volume>124</volume>, <fpage>102293</fpage>. <pub-id pub-id-type="doi">10.1016/j.apgeog.2020.102293</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. L.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J. S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y. Q.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Evaluation for regional ecosystem health: Methodology and research progress</article-title>. <source>Acta Ecol. Sin.</source> <volume>27</volume> (<issue>11</issue>), <fpage>4877</fpage>&#x2013;<lpage>4885</lpage>. <pub-id pub-id-type="doi">10.1016/s1872-2032(08)60009-8</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rapport</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Constanza</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Mcmichaeld</surname>
<given-names>A. J.</given-names>
</name>
</person-group> (<year>1998</year>). <article-title>Assessing ecosystem health</article-title>. <source>Trends Ecol. Evol.</source> <volume>13</volume> (<issue>10</issue>), <fpage>397</fpage>&#x2013;<lpage>402</lpage>. <pub-id pub-id-type="doi">10.1016/s0169-5347(98)01449-9</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rapport</surname>
<given-names>D. J.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>What constitutes ecosystem health?</article-title> <source>Perspect. Biol. Med.</source> <volume>33</volume> (<issue>1</issue>), <fpage>120</fpage>&#x2013;<lpage>132</lpage>. <pub-id pub-id-type="doi">10.1353/pbm.1990.0004</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ren</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J. G.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>S. L.</given-names>
</name>
</person-group> (<year>2000</year>). <article-title>Evaluation and monitoring of ecosystem health</article-title>. <source>Trop. Geogr.</source> <volume>20</volume> (<issue>4</issue>), <fpage>310</fpage>&#x2013;<lpage>316</lpage>.</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Robert</surname>
<given-names>M. G.</given-names>
</name>
<name>
<surname>Robert</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Gardner</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Bruce</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Effects of changing spatial scale on the analysis of landscape pattern</article-title>. <source>Landsc. Ecol.</source> <volume>3</volume> (<issue>3-4</issue>), <fpage>153</fpage>&#x2013;<lpage>162</lpage>. <pub-id pub-id-type="doi">10.1007/bf00131534</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Schaeffer</surname>
<given-names>D. J.</given-names>
</name>
<name>
<surname>Henricks</surname>
<given-names>E. E.</given-names>
</name>
<name>
<surname>Kerster</surname>
<given-names>H. W.</given-names>
</name>
</person-group> (<year>1988</year>). <article-title>Ecosystem health: I. Measuring ecosystem health</article-title>. <source>Environ. Manag.</source> <volume>12</volume> (<issue>3</issue>), <fpage>445</fpage>&#x2013;<lpage>455</lpage>. <pub-id pub-id-type="doi">10.1007/bf01873258</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>Z. H.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>Y. C.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>A integrated method for assessing the urban ecosystem health of rapid urbanized area in China based on SFPHD framework</article-title>. <source>Ecol. Indic.</source> <volume>121</volume>, <fpage>107071</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.107071</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sheng</surname>
<given-names>Z. L.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>X. Q.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>P. Z.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Progress on watershed ecosystem health assessment research in China</article-title>. <source>Yunnan Geogr. Environ. Res.</source> <volume>23</volume> (<issue>2</issue>), <fpage>52</fpage>&#x2013;<lpage>58</lpage>.</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>Y. H.</given-names>
</name>
<name>
<surname>Di</surname>
<given-names>X. H.</given-names>
</name>
<name>
<surname>Ming</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>B. T.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Ecological risk assessment in Yongchuan District of Chongqing city based on landscape pattern</article-title>. <source>Bull. Soil Water Conservation</source> <volume>40</volume> (<issue>3</issue>), <fpage>195</fpage>&#x2013;<lpage>201&#x2b;215</lpage>.</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>X. C.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Ecological space control of Dianchi Lake basin based on ecological protection red line</article-title>. <source>Econ. Geogr.</source> <volume>40</volume> (<issue>5</issue>), <fpage>191</fpage>&#x2013;<lpage>197</lpage>.</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Effects of changing scale on landscape pattern analysis: Scaling relations</article-title>. <source>Landsc. Ecol.</source> <volume>19</volume> (<issue>2</issue>), <fpage>125</fpage>&#x2013;<lpage>138</lpage>. <pub-id pub-id-type="doi">10.1023/b:land.0000021711.40074.ae</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Tueller</surname>
<given-names>P. T.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Empirical patterns of the effects of changing scale on landscape metrics</article-title>. <source>Landsc. Ecol.</source> <volume>17</volume> (<issue>8</issue>), <fpage>761</fpage>&#x2013;<lpage>782</lpage>. <pub-id pub-id-type="doi">10.1023/a:1022995922992</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>X. F.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>G. J.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H. Y.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Ecosystem health assessment of dongyang River basin</article-title>. <source>Bull. Soil Water conservation</source> <volume>35</volume> (<issue>2</issue>), <fpage>326</fpage>&#x2013;<lpage>331</lpage>.</citation>
</ref>
<ref id="B40">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname>
<given-names>D. N.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X. Z.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>T. S.</given-names>
</name>
</person-group> (<year>2003</year>). <source>Landscape ecology</source>. <publisher-loc>Beijing</publisher-loc>: <publisher-name>Science Press</publisher-name>.</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiao</surname>
<given-names>F. J.</given-names>
</name>
<name>
<surname>Ou</surname>
<given-names>Y. H.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Ecosystem health and its evolution indicator and method</article-title>. <source>J. Nat. Resour.</source> <volume>17</volume> (<issue>2</issue>), <fpage>203</fpage>&#x2013;<lpage>209</lpage>.</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>X. F.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>G. J.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Bian</surname>
<given-names>H. J.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2015</year>). <article-title>Evaluation of ecosystem health in Western Tiaoxi River watershed based on matter element model</article-title>. <source>Acta Sci. Circumstantiae</source> <volume>35</volume> (<issue>4</issue>), <fpage>1250</fpage>&#x2013;<lpage>1258</lpage>.</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>S. L.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Re-thinking of reestablish eco-system health on watershed scale during the process of urbanization</article-title>. <source>Yunnan Geogr. Environ. Res.</source> <volume>23</volume> (<issue>6</issue>), <fpage>70</fpage>&#x2013;<lpage>73&#x2b;85</lpage>.</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yan</surname>
<given-names>S. G.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H. D.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Ecosystem health assessment method of eco-redline based on land use and landscape pattern in Nanjing</article-title>. <source>J. Nat. Resour.</source> <volume>35</volume> (<issue>5</issue>), <fpage>1109</fpage>&#x2013;<lpage>1118</lpage>. <pub-id pub-id-type="doi">10.31497/zrzyxb.20200508</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bi</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Jing</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhen</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Impact of mining on ecosystem health in loess plateau&#x2014;a case study on changhe basin in jincheng city</article-title>. <source>J. Shanxi Agric. Univ. Nat. Sci. Ed.</source> <volume>37</volume> (<issue>4</issue>), <fpage>268</fpage>&#x2013;<lpage>275</lpage>.</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yuan</surname>
<given-names>M. N.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y. X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Peng</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Ecosystem health assessment based on the framework of vigor, organization, resilience and contribution in Guangzhou city</article-title>. <source>Chin. J. Ecol.</source> <volume>38</volume> (<issue>4</issue>), <fpage>1249</fpage>&#x2013;<lpage>1257</lpage>.</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>D. M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>H. L.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2013a</year>). <article-title>Relationship between landscape pattern of construction land in Dianchi Lake basin and water quality in Dianchi Lake</article-title>. <source>Bull. Soil Water Conservation</source> <volume>33</volume> (<issue>4</issue>), <fpage>103</fpage>&#x2013;<lpage>107&#x2b;4</lpage>.</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>L. L.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>Y. F.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y. H.</given-names>
</name>
</person-group> (<year>2013b</year>). <article-title>Effects of spatial grain change on the landscape pattern indices in Yimeng Mountain area of Shandong Province, East China</article-title>. <source>Chin. J. Ecol.</source> <volume>32</volume> (<issue>2</issue>), <fpage>459</fpage>&#x2013;<lpage>464</lpage>.</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2020a</year>). <article-title>The spatial granularity effect, changing landscape patterns, and suitable landscape metrics in the three gorges reservoir area, 1995-2015</article-title>. <source>Ecol. Indic.</source> <volume>114</volume> (<issue>12</issue>), <fpage>106259</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.106259</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>M. S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>T. H.</given-names>
</name>
<name>
<surname>Cui</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>M. H.</given-names>
</name>
<etal/>
</person-group> (<year>2020b</year>). <article-title>Ecosystem health assessment and eco-restoration mode discussion for lakeside zone of Taihu lake</article-title>. <source>J. Hydroecology</source> <volume>41</volume> (<issue>4</issue>), <fpage>48</fpage>&#x2013;<lpage>54</lpage>.</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ban</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Xiao</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2020c</year>). <article-title>Evaluating ecological health in the middle-lower reaches of the Hanjiang River with cascade reservoirs using the Planktonic index of biotic integrity (P-IBI)</article-title>. <source>Ecol. Indic.</source> <volume>114</volume>, <fpage>106282</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2020.106282</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>Y. L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y. C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Examining land-use/land-cover change in the lake Dianchi watershed of the yunnan&#x2013;guizhou plateau of southwest China with remote sensing and GIS techniques: 1974-2008</article-title>. <source>Int. J. Environ. Res. Public Health</source> <volume>9</volume> (<issue>11</issue>), <fpage>3843</fpage>&#x2013;<lpage>3865</lpage>. <pub-id pub-id-type="doi">10.3390/ijerph9113843</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>