<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" article-type="research-article" dtd-version="2.3" xml:lang="EN">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Ecol. Evol.</journal-id>
<journal-title>Frontiers in Ecology and Evolution</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Ecol. Evol.</abbrev-journal-title>
<issn pub-type="epub">2296-701X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fevo.2023.1200218</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Ecology and Evolution</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Ecological space resilience assessment of Baiyangdian Basin from the perspective of evolutionary resilience</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liang</surname>
<given-names>Xinyu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2198224"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Gao</surname>
<given-names>Xing</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<xref ref-type="author-notes" rid="fn001">
<sup>*</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Zhao</surname>
<given-names>Meiran</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Guo</surname>
<given-names>Zhongyuan</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Shen</surname>
<given-names>Weining</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Public Administration, Hebei University of Economics and Business</institution>, <addr-line>Shijiazhuang</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Hebei Collaborative Innovation Center for Urban-Rural Integrated Development, Hebei University of Economics and Business</institution>, <addr-line>Shijiazhuang</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>Hebei Province Economics and Social Development Geographic Information Big Data Platform, Hebei University of Economics and Business</institution>, <addr-line>Shijiazhuang</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Center for Urban Sustainability and Innovation Development (CUSID), Hebei University of Economics and Business</institution>, <addr-line>Shijiazhuang</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>Edited by: Hualin Xie, Jiangxi University of Finance and Economics, China</p>
</fn>
<fn fn-type="edited-by">
<p>Reviewed by: Zhe Feng, China University of Geosciences, China; Pattiyage Ishan Ayantha Gomes, Sri Lanka Institute of Information Technology, Sri Lanka</p>
</fn>
<fn fn-type="corresp" id="fn001">
<p>*Correspondence: Xing Gao, <email xlink:href="mailto:gaoxing@hueb.edu.cn">gaoxing@hueb.edu.cn</email>
</p>
</fn>
</author-notes>
<pub-date pub-type="epub">
<day>27</day>
<month>10</month>
<year>2023</year>
</pub-date>
<pub-date pub-type="collection">
<year>2023</year>
</pub-date>
<volume>11</volume>
<elocation-id>1200218</elocation-id>
<history>
<date date-type="received">
<day>04</day>
<month>04</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>17</day>
<month>10</month>
<year>2023</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2023 Liang, Gao, Zhao, Guo and Shen</copyright-statement>
<copyright-year>2023</copyright-year>
<copyright-holder>Liang, Gao, Zhao, Guo and Shen</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>The structure and function of ecological spaces play a crucial role in resisting risk interference and providing ecological services. Conducting resilience assessments of ecological spaces is of great significance for maintaining ecological security. Taking the Baiyangdian Basin as the research object, based on the evaluation results of ecosystem service function importance and ecological environmental sensitivity, the ecological space importance level of the Baiyangdian Basin is divided. Based on the three basic characteristics of &#x201c;resistance, adaptability, recoverability&#x201d;, the resilience of ecological space model was constructed by selecting ecosystem service value, landscape structure stability, ecological resilience model and coefficient. The research results show that: (1) The overall importance index of ecological space in the Baiyangdian Basin is relatively low. The non-critical area accounted for 43.00% of the total area, the general area accounted for 30.78%, and the core area accounted for 26.22%, which mainly distributed in the northern part of the study area and at the confluence of rivers. (2) The proportion of areas with extremely important ecosystem service functions in Baiyangdian Basin is 31.73%, concentrated in the central and southwestern parts of the basin. The highly sensitive ecological environment areas accounted for 21.96%, mainly distributed in high-altitude areas in the northern part of the basin and river convergence areas. (3) On the whole, the overall resilience level of the study area is not high and there is a significant difference in resilience between the north and south, with 68.60% of the regions in the low level of resilience. The area with high resilience level accounted for only 1.51%, and the highest resilience level was found at river confluence. The research results provide a theoretical basis for maintaining and improving the ecological environment of the basin, and provide a basis for implementing precise policies to optimize the ecological space of the basin.</p>
</abstract>
<kwd-group>
<kwd>ecological space</kwd>
<kwd>resilience</kwd>
<kwd>ecosystem service</kwd>
<kwd>landscape pattern</kwd>
<kwd>Xiong&#x2019;an New Area</kwd>
</kwd-group>
<counts>
<fig-count count="7"/>
<table-count count="7"/>
<equation-count count="11"/>
<ref-count count="48"/>
<page-count count="14"/>
<word-count count="6835"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-in-acceptance</meta-name>
<meta-value>Urban Ecology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1" sec-type="intro">
<label>1</label>
<title>Introduction</title>
<p>Ecological space is the basic space to maintain ecosystem balance and ensure regional ecological security (<xref ref-type="bibr" rid="B45">Zhang and Yue, 2019</xref>). At present, climate change, socio-economic development, construction land expansion and other phenomena have brought about many problems such as water resource shortage, reduction of biological and landscape diversity, and continuous degradation of ecological functions (<xref ref-type="bibr" rid="B22">Liu et&#xa0;al., 2016</xref>; <xref ref-type="bibr" rid="B17">Jin et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B11">Gao et&#xa0;al., 2021</xref>; <xref ref-type="bibr" rid="B20">Liu et&#xa0;al., 2022</xref>), resulting in unreasonable distribution of ecological space, affecting the ecological environment and the well-being of residents (<xref ref-type="bibr" rid="B23">Liu et&#xa0;al., 2023</xref>). Under this background, it is urgent to change the thinking of ecological space management and control, from the perspective of restoration and management of ecological space after damage to prevention and control of ecological space before damage, and enhance the ability of ecological space to cope with the impact of uncertain risks. Resilience theory, as a dynamic approach to analyzing how systems respond to risks, provides a new perspective for ecological spaces to deal with uncertain risks (<xref ref-type="bibr" rid="B25">Luo et&#xa0;al., 2018</xref>). &#x201c;Ecological space resilience&#x201d; is a new attempt to combine the concept of resilience with the study of ecological space. The concept can be defined as the ability of ecological space to resist and adapt to the impact of external shocks and recover or transform from them by optimizing its structure and function. Therefore, the resilience assessment of ecological space can provide a theoretical basis for improving its resilience level and achieving sustainable development.</p>
<p>Ecological space originates from the term &#x201c;green space&#x201d; (<xref ref-type="bibr" rid="B26">Ngom et&#xa0;al., 2016</xref>). It is an important space component to maintain biodiversity, ensure ecosystem service supply and improve residents&#x2019; quality of life (<xref ref-type="bibr" rid="B39">Xie et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B8">Fabien et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B6">Deng et&#xa0;al., 2022</xref>). The follow-up research of ecological space should be based on the quantitative identification of ecological space. In recent years, the research of ecological space identification has made some progress. From the perspective of research objects, most of them focus on meso-scale and macro-scale such as urban agglomeration and municipal level, while few studies are conducted from the watershed scale (<xref ref-type="bibr" rid="B15">Huang et&#xa0;al., 2020</xref>; <xref ref-type="bibr" rid="B35">Wang et&#xa0;al., 2022</xref>). From the perspective of research, the structure and function of ecological space are mainly studied from the perspectives of ecological functions, ecological elements and ecological protection red lines (<xref ref-type="bibr" rid="B24">Long et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B42">Yu et&#xa0;al., 2015</xref>; <xref ref-type="bibr" rid="B19">Kong et&#xa0;al., 2019</xref>; <xref ref-type="bibr" rid="B28">Pan et&#xa0;al., 2020</xref>). From the perspective of ecological space delineation method, the main methods are land use type merging method and spatial attribute evaluation method. Among them, land use type merging method is to identify the structure of ecological space by using different land use types (<xref ref-type="bibr" rid="B16">Jiang and Liu, 2020</xref>; <xref ref-type="bibr" rid="B47">Zhao et&#xa0;al., 2022</xref>), but this method can draw the amount of ecological space, but it is difficult to draw the function of ecological space. Spatial attribute method is mainly used to quantitatively identify ecological space through ecosystem service function evaluation, ecological sensitivity evaluation and other methods (<xref ref-type="bibr" rid="B18">Jin et&#xa0;al., 2020</xref>). This method can more intuitively reflect the structure and function of ecological space, and is currently an important way to identify ecological space.</p>
<p>Originally meant to &#x201c;restore vitality&#x201d;, resilience was introduced into ecology by Canadian scholar Holling in 1973 to characterize the characteristics of stable state of ecosystem (<xref ref-type="bibr" rid="B13">Holling, 1973</xref>). Since the concept of resilience was put forward, it has undergone step modification and development from engineering resilience to ecological resilience and then to evolution resilience. Engineering resilience is characterized by the ability of the system to recover to a balanced or stable state after external impact (<xref ref-type="bibr" rid="B34">Walker and Salt, 2006</xref>), focusing on a single ultimate equilibrium state (<xref ref-type="bibr" rid="B41">Yang et&#xa0;al., 2021</xref>). With the deepening of the research, the academic community gradually realized that resilience has a certain ability to absorb disturbance before it suffers from impact and changes its own state (<xref ref-type="bibr" rid="B9">Folke, 2006</xref>), and ecological resilience to mitigate impact, resist risks, restore balance and adapt to new paths has become the mainstream view (<xref ref-type="bibr" rid="B14">Holling, 1996</xref>). With the development of The Times and the frequent occurrence of various crises, scholars have new thoughts on resilience and gradually realize that the path of system development is not single or balanced, but complex and unbalanced. Therefore, evolution resilience, as a resilience theory that emphasizes system interaction and dynamic feedback, immediately becomes a mainstream paradigm widely recognized by the current academic community (<xref ref-type="bibr" rid="B33">Walker et&#xa0;al., 2004</xref>). In the perspective of evolution resilience, resilience is not the recovery to normal state, but the ability of social complex ecosystem to adjust, adapt and transform in response to pressure (<xref ref-type="bibr" rid="B1">Carpenter et&#xa0;al., 2005</xref>; <xref ref-type="bibr" rid="B48">Zhou, 2015</xref>). Evolutionary resilience focuses on the maintenance of system structure and function, focusing on three aspects of social ecosystems: resistance, adaptability, and recoverability. Under the concept of evolutionary resilience, the concept of resilience has been applied in multiple fields. In recent years, the research objects on resilience involve urban resilience, economic resilience, ecological resilience, organizational resilience and other aspects. As modern society is facing threats from natural disasters, public security and other aspects, the current research on urban resilience accounts for the majority. The research methods can be divided into qualitative research and quantitative research. Among them, the qualitative research mainly focuses on the concept and analysis framework of resilience (<xref ref-type="bibr" rid="B37">Wei and Xiu, 2020</xref>) and the formation mechanism of resilience (<xref ref-type="bibr" rid="B7">Du et&#xa0;al., 2022</xref>). Quantitative studies mostly focus on building index systems to assess the level of resilience and analyze its influencing factors. Most scholars explain the spatial distribution of resilience level from the three core characteristics of resilience (namely resistance, adaptability and recoverability) (<xref ref-type="bibr" rid="B38">Xia et&#xa0;al., 2022</xref>), among which, Resistance is the ability of the system to withstand disturbances without significantly deviating from the normal development path (<xref ref-type="bibr" rid="B32">Vugrin et&#xa0;al., 2011</xref>), and adaptability is the ability of the system to adapt to shocks in the normal development process (<xref ref-type="bibr" rid="B2">Chen et&#xa0;al., 2022</xref>). Recoverability is the ability of a system to quickly recover from a disaster with minimal self-damage (<xref ref-type="bibr" rid="B12">Haimes Yacov, 2009</xref>). The combination of resilience and ecological space gives rise to the concept of ecological space resilience. Resilience can be seen as an attribute of ecological space, which can measure the ability of ecological space to face future disturbances, optimize its structure and function by taking preemptive actions to reduce the impact of disturbances, and quickly respond and restore to a new equilibrium state when impacts occur (<xref ref-type="fig" rid="f1">
<bold>Figure&#xa0;1</bold>
</xref>). At present, although there is a high degree of recognition of the concept of resilience in the existing literature, the grasp and emphasis of resilience vary. There is a lack of research on the measurement and comprehensive evaluation of ecological spatial resilience, and a complete theoretical framework has not yet been formed. Further exploration is needed to incorporate the theory of evolutionary resilience into the practice of ecological spatial optimization layout.</p>
<fig id="f1" position="float">
<label>Figure&#xa0;1</label>
<caption>
<p>Understanding about Ecological Space Resilience Concepts.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g001.tif"/>
</fig>
<p>As the ecological hinterland of Xiong&#x2019;an New Area, the structure and function of the ecological space in the Baiyangdian Basin are crucial for maintaining and improving the regional ecological environment. With the improvement of urbanization level, the ecological environment of Baiyangdian Basin is facing increasing pressure, and the function of ecological space cannot be effectively played. It is necessary to transform the concept of ecological problem restoration into the concept of ecological problem prevention, and integrate resilience thinking into the construction of ecological space to enhance the adaptability of the basin to respond to uncertain risks. In this paper, 35 counties and cities in five prefectural cities in Hebei Province, China, through which Baiyangdian Basin flows, are studied. Based on two evaluation methods, the importance of ecosystem services and the sensitivity of ecological environment, the ecological spatial importance index system is established to quantitatively identify the ecological space of Baiyangdian Basin. Based on the three basic characteristics of &#x201c;resistance, adaptability and recoverability&#x201d;, the resilience evaluation index was selected to analyze the spatial distribution of resilience level of ecological space structure and function, and divide resilience zones, in order to provide quantitative basis for improving the ecological space structure and ecological space function of Baiyangdian Basin and realize the sustainable development of the basin.</p>
</sec>
<sec id="s2" sec-type="materials|methods">
<label>2</label>
<title>Materials and methods</title>
<sec id="s2_1">
<label>2.1</label>
<title>Overview of the study area</title>
<p>The Baiyangdian Basin in Hebei Province, China, as defined in this study, is a part of the Daqing River Basin (<xref ref-type="fig" rid="f2">
<bold>Figure&#xa0;2</bold>
</xref>). The study area is located at 113&#xb0;46&#x2019; -116&#xb0;25&#x2019; east longitude and 37&#xb0;51&#x2019; - 40&#xb0;29&#x2019; north latitude. The terrain presents a stepped downward trend from northwest to southeast, showing a mountain-hill-plain with a relative height of 2815 meters. The average temperature is 12.1&#xb0;C, the average annual water storage is 1.32&#xd7;10<sup>9</sup>m<sup>3</sup>. By 2020, the total area of the study area is 4.43&#xd7;10<sup>6</sup>hm<sup>2</sup>, among which the cultivated land area is 2.51&#xd7;10<sup>6</sup>hm<sup>2</sup>, accounting for 56.58% of the total land use type. The area of forest land was 1.01&#xd7;10<sup>6</sup>hm<sup>2</sup>, accounting for 22.84%. The grassland area was 2.47&#xd7;10<sup>5</sup>hm<sup>2</sup>, accounting for 5.56%. The water area was 4.56&#xd7;10<sup>4</sup>hm<sup>2</sup>, accounting for 1.03%; The construction land area is 6.19&#xd7;10<sup>5</sup>hm<sup>2</sup>, accounting for 13.96%. The unused land area is 1.1&#xd7;10<sup>3</sup>hm<sup>2</sup>, accounting for 0.02%. In general, cultivated land is the main study area, followed by woodland and construction land. In recent years, the structure and function of the ecological space in the Baiyangdian Basin have been affected to a certain extent. Studying the resilience level is important to ensure the ecological security of the Baiyangdian Basin.</p>
<fig id="f2" position="float">
<label>Figure&#xa0;2</label>
<caption>
<p>Spatial location of the study area (<bold>(A)</bold> China scope; <bold>(B)</bold> Hebei Province; <bold>(C)</bold> Baiyangdian Basin).</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g002.tif"/>
</fig>
</sec>
<sec id="s2_2">
<label>2.2</label>
<title>Data source and preprocessing</title>
<p>The data used in this study include remote sensing image data, land use data, DEM elevation data, slope data, NDVI data, soil texture data, meteorological data, NPP data and socio-economic data. Among them:</p>
<list list-type="order">
<list-item>
<p>The remote sensing image data of the study area in 2020 came from the geospatial data cloud (<ext-link ext-link-type="uri" xlink:href="https://www.gscloud.cn/">https://www.gscloud.cn/</ext-link>).</p>
</list-item>
<list-item>
<p>Land use data was interpreted from remote sensing images. ENVI software was used to preprocess, supervise and classify remote sensing images, and land use types were divided into cultivated land, forest land, grassland, water area, construction land and unused land. After testing, the accuracy of the image is higher than 85%, which can be used in this study.</p>
</list-item>
<list-item>
<p>DEM elevation data comes from geospatial data cloud (<ext-link ext-link-type="uri" xlink:href="https://www.gscloud.cn/">https://www.gscloud.cn/</ext-link>) with a resolution of 30m and a time of 2020. Slope data are extracted from DEM elevation data.</p>
</list-item>
<list-item>
<p>All NDVI data and soil data were obtained from Data Center for Resources and Environmental Sciences, Chinese Academy of Sciences (<ext-link ext-link-type="uri" xlink:href="https://www.resdc.cn/">https://www.resdc.cn/</ext-link>). ArcGIS software was used to carry out reprojection and mask processing on Chinese soil data, extract soil data in the study area, connect the &#x201c;Value&#x201d; field in its attribute table with the &#x201c;MU_GLOBAL&#x201d; field in the soil attribute table HWSD.mdb, and obtain the raster data of soil texture in the study area.</p>
</list-item>
<list-item>
<p>Meteorological source: National Data Center for Meteorological Sciences (<ext-link ext-link-type="uri" xlink:href="https://data.cma.cn/">https://data.cma.cn/</ext-link>). In the grid calculator, the maximum and minimum method is used to normalize the annual average rainfall, and the distribution map of annual average precipitation in the study area is obtained. The annual average temperature data of several meteorological stations around the study area was obtained by spatial interpolation calculation.</p>
</list-item>
<list-item>
<p>NPP data from MODIS data set (<ext-link ext-link-type="uri" xlink:href="https://ladsweb.modaps.eosdis.nasa.gov/">https://ladsweb.modaps.eosdis.nasa.gov/</ext-link>). The MRT tool was used to extract the NPP data. Referring to the MODIS17A3 data document, the conversion factor 0.0001 was used in the raster calculator to process the obtained NPP data, and the mask extraction was carried out based on the vector boundary of the study area in ArcGIS software. Because of the holes in the data, raster point transformation was carried out, and Kriging method was used to carry out spatial interpolation on the point data, and finally resampling the data after interpolation was processed to obtain data with a spatial resolution of 30m&#xd7;30m.</p>
</list-item>
<list-item>
<p>Socio-economic data are derived from statistical yearbooks.</p>
</list-item>
</list>
</sec>
<sec id="s2_3">
<label>2.3</label>
<title>Research methods</title>
<sec id="s2_3_1">
<label>2.3.1</label>
<title>Identification of ecological space</title>
<p>Ecological space is a space that provides ecosystem services (<xref ref-type="bibr" rid="B36">Wang et&#xa0;al., 2017</xref>; <xref ref-type="bibr" rid="B5">Cui et&#xa0;al., 2020</xref>), whose quantitative scale and spatial pattern have a great impact on the ecological security of the watershed (<xref ref-type="bibr" rid="B3">Chen et&#xa0;al., 2018</xref>). In view of a series of ecological problems in the study area, the study of ecological space should not only pay attention to the &#x201c;structure&#x201d;, but also to the &#x201c;function&#x201d; of ecological space. Therefore, the factors that can reflect the ecological spatial function of Baiyangdian basin are selected. Based on the two evaluation methods of the importance of ecosystem service function and the sensitivity of ecological environment, the ecological space is identified comprehensively and quantitatively, and the importance degree is graded.</p>
<sec id="s2_3_1_1">
<label>2.3.1.1</label>
<title>Evaluation of the importance of ecosystem service functions</title>
<p>(1) Importance index of water conservation function</p>
<p>The water conservation function is an important function in ecosystem regulation services, and its change directly affects the hydrology and vegetation of the basin (<xref ref-type="bibr" rid="B44">Zhang et&#xa0;al., 2023</xref>). The calculation formula is:</p>
<disp-formula>
<label>(1)</label>
<mml:math display="block" id="M1">
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>N</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>c</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>s</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>WR</italic> is the ecosystem water conservation function importance index, <italic>NPP<sub>mean</sub>
</italic> is the average annual net primary productivity of vegetation, <italic>F<sub>sic</sub>
</italic>is the soil seepage factor, <italic>F<sub>pre</sub>
</italic> is the average annual precipitation factor, <italic>F<sub>slo</sub>
</italic> is the slope factor.</p>
<p>(2) Importance index of soil and water conservation function</p>
<p>The function of soil and water conservation is the regulation service of ecosystem (such as forest, grassland, etc.) to reduce soil erosion caused by water erosion through its structure (<xref ref-type="bibr" rid="B46">Zhang et&#xa0;al., 2015</xref>). The modified general soil and water loss equation is used to measure the soil and water conservation function. The calculation formula is:</p>
<disp-formula>
<label>(2)</label>
<mml:math display="block" id="M2">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>R</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>K</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>L</mml:mi>
<mml:mi>S</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>C</mml:mi>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>A</italic> is the importance index of soil and water conservation function, <italic>R</italic> is rainfall erosivity factor, <italic>K</italic> is soil erodibility factor, <italic>LS</italic> is the topographic relief factor, <italic>C</italic> is the vegetation cover factor.</p>
<p>(3) Importance index of biodiversity conservation function</p>
<p>Biodiversity conservation function refers to the role played by an ecosystem in maintaining the diversity of genes, species and ecosystems, and is one of the most important functions provided by an ecosystem (<xref ref-type="bibr" rid="B21">Liu and Liu, 2021</xref>). The calculation formula is:</p>
<disp-formula>
<label>(3)</label>
<mml:math display="block" id="M3">
<mml:mrow>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mi>N</mml:mi>
<mml:mi>P</mml:mi>
<mml:msub>
<mml:mi>P</mml:mi>
<mml:mrow>
<mml:mi>m</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>p</mml:mi>
<mml:mi>r</mml:mi>
<mml:mi>e</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>t</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>m</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mo stretchy="false">(</mml:mo>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mrow>
<mml:mi>a</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>t</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mo stretchy="false">)</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>S<sub>bio</sub>
</italic>is the importance index of biodiversity conservation function, <italic>NPP<sub>mean</sub>
</italic> is the average annual net primary productivity of vegetation, <italic>F<sub>pre</sub>
</italic> is the average annual precipitation, <italic>F<sub>alt</sub>
</italic> is the average annual temperature, and is the altitude factor.</p>
<p>(4) Evaluation of the importance of ecosystem service functions</p>
<p>The ecosystem service value after the normalization of each single function is divided into three levels: unimportant, generally important, and extremely important. Then the ecosystem service function importance index (<inline-formula>
<mml:math display="inline" id="im1">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) is obtained by using the disjunction algorithm. Using the natural breakpoint method, the results of the importance of ecosystem service functions are divided into three levels: unimportant, generally important, and extremely important. The calculation formula is:</p>
<disp-formula>
<label>(4)</label>
<mml:math display="block" id="M4">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mi>W</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mrow>
<mml:mi>b</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>o</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
<sec id="s2_3_1_2">
<label>2.3.1.2</label>
<title>Ecological environment sensitivity assessment</title>
<p>The sensitivity of ecological environment can reflect the sensitivity of ecosystem to natural environment changes and interference from human activities (<xref ref-type="bibr" rid="B27">Ouyang et&#xa0;al., 2000</xref>). The terrain of the Baiyangdian Basin is complex, and its ecological environment sensitivity is influenced by various ecological factors. Therefore, dominant and representative factors are selected to comprehensively evaluate the ecological environment sensitivity (<xref ref-type="bibr" rid="B43">Zhan and Zhu, 2019</xref>), and the sensitivity degree is divided.</p>
<list list-type="order">
<list-item>
<p>Build a hierarchical model. The model is divided into three levels: target layer, criterion layer and reference layer (<xref ref-type="table" rid="T1"><bold>Table 1</bold></xref>). The target layer is the ecological sensitivity assessment of Baiyangdian Basin, and the criterion layer includes two judgment criteria: terrain factor and land use type. The reference layer can be determined by the following indicators:</p>
</list-item>
<list-item>
<p>Sensitivity level definition. Determine the weights of each factor using the Analytic Hierarchy Process and define the sensitivity level of each ecological factor (<xref ref-type="table" rid="T2"><bold>Table 2</bold></xref>).</p>
</list-item>
<list-item>
<p>The ecological environmental sensitivity (<inline-formula>
<mml:math display="inline" id="im2">
<mml:mrow>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) of the study area is obtained by weighted superposition of the selected factors. The natural breakpoint method is used to reclassify the results, resulting in three levels of non-sensitivity, medium sensitivity, and high sensitivity.</p>
</list-item>
</list>
<table-wrap id="T1" position="float">
<label>Table&#xa0;1</label>
<caption>
<p>Sensitivity information of ecological Factors in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Target layer</th>
<th valign="middle" align="center">Criterion layer</th>
<th valign="middle" align="center">Reference layer</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="5" align="center">Ecological sensitivity assessment of Baiyangdian Basin</td>
<td valign="middle" rowspan="3" align="center">terrain factor</td>
<td valign="middle" align="center">slope/&#xb0;</td>
</tr>
<tr>
<td valign="middle" align="center">elevation/m</td>
</tr>
<tr>
<td valign="middle" align="center">Slope direction</td>
</tr>
<tr>
<td valign="middle" rowspan="2" align="center">land use type</td>
<td valign="middle" align="center">vegetation index</td>
</tr>
<tr>
<td valign="middle" align="center">land type</td>
</tr>
</tbody>
</table>
</table-wrap>
<table-wrap id="T2" position="float">
<label>Table&#xa0;2</label>
<caption>
<p>Sensitivity information of various ecological factors in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center">Ecological factors</th>
<th valign="middle" align="center">Secondary factor</th>
<th valign="middle" align="center">Classification</th>
<th valign="middle" align="center">Sensitivity level</th>
<th valign="middle" align="center">Weight</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" rowspan="9" align="center">terrain factor</td>
<td valign="middle" rowspan="3" align="center">slope/&#xb0;</td>
<td valign="middle" align="center">0&#x2014;7.21</td>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" rowspan="3" align="center">0.159</td>
</tr>
<tr>
<td valign="middle" align="center">7.21&#x2014;19.61</td>
<td valign="middle" align="center">medium sensitivity</td>
</tr>
<tr>
<td valign="middle" align="center">19.61&#x2014;57.49</td>
<td valign="middle" align="center">non-sensitive</td>
</tr>
<tr>
<td valign="middle" rowspan="3" align="center">elevation/m</td>
<td valign="middle" align="center">-104&#x2014;441</td>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" rowspan="3" align="center">0.195</td>
</tr>
<tr>
<td valign="middle" align="center">441&#x2014;1134</td>
<td valign="middle" align="center">medium sensitivity</td>
</tr>
<tr>
<td valign="middle" align="center">1134&#x2014;2832</td>
<td valign="middle" align="center">non-sensitive</td>
</tr>
<tr>
<td valign="middle" rowspan="3" align="center">slope direction</td>
<td valign="middle" align="center">flat ground, due south</td>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" rowspan="3" align="center">0.194</td>
</tr>
<tr>
<td valign="middle" align="center">due east, due west, southeast, southwest</td>
<td valign="middle" align="center">medium sensitivity</td>
</tr>
<tr>
<td valign="middle" align="center">due north, northeast, northwest</td>
<td valign="middle" align="center">non-sensitive</td>
</tr>
<tr>
<td valign="middle" rowspan="6" align="center">land use type</td>
<td valign="middle" rowspan="3" align="center">vegetation index</td>
<td valign="middle" align="center">0.00&#x2014;0.44</td>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" rowspan="3" align="center">0.240</td>
</tr>
<tr>
<td valign="middle" align="center">0.44&#x2014;0.64</td>
<td valign="middle" align="center">medium sensitivity</td>
</tr>
<tr>
<td valign="middle" align="center">0.64&#x2014;0.90</td>
<td valign="middle" align="center">non-sensitive</td>
</tr>
<tr>
<td valign="middle" rowspan="3" align="center">land type</td>
<td valign="middle" align="center">forest land, water area</td>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" rowspan="3" align="center">0.212</td>
</tr>
<tr>
<td valign="middle" align="center">cultivated land, grassland</td>
<td valign="middle" align="center">medium sensitivity</td>
</tr>
<tr>
<td valign="middle" align="center">construction land, unused land</td>
<td valign="middle" align="center">non-sensitive</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2_3_1_3">
<label>2.3.1.3</label>
<title>Comprehensive identification of ecological space</title>
<p>In the ArcGIS software, ecological space (<inline-formula>
<mml:math display="inline" id="im3">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>L</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) is obtained by disjunctive algorithm, and it is divided into three levels: non-critical type, general type, and core type. The calculation formula is:</p>
<disp-formula>
<label>(5)</label>
<mml:math display="block" id="M5">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>L</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>m</mml:mi>
<mml:mi>a</mml:mi>
<mml:mi>x</mml:mi>
<mml:mo>{</mml:mo>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>,</mml:mo>
<mml:mo>&#xa0;</mml:mo>
<mml:mi>S</mml:mi>
<mml:mi>S</mml:mi>
</mml:mrow>
<mml:mo>}</mml:mo>
</mml:mrow>
</mml:math>
</disp-formula>
</sec>
</sec>
<sec id="s2_3_2">
<label>2.3.2</label>
<title>Horizontal measure of resilience</title>
<p>In view of the negative impact brought by the disturbance, the concept of resilience is introduced to improve the ability of the study area to cope with various risks. From the perspective of evolutionary resilience, resilience refers to the persistence, adaptability, and transferability of complex social ecosystems, focusing on the ability of ecological space to respond, adapt, and shape changes. Its emphasis is not on the pursuit of equilibrium, but on the ability of the system to continuously improve its own structure to cope with risks, thus enabling sustainable development. Therefore, constructing resilience evaluation indicators by considering the three basic characteristics of &#x201c;resistance, adaptability and recoverability&#x201d; in evolutionary resilience, and the horizontal spatial distribution of resilience in ecological space was analyzed.</p>
<sec id="s2_3_2_1">
<label>2.3.2.1</label>
<title>Resistance index</title>
<p>Resistance (<italic>P</italic>) refers to the ability of ecological space to resist external disturbance. Many studies have pointed out that the resistance characteristic of ecological space is closely related to ecosystem service function. This study is based on the modified ecological service equivalence table per unit area of Chinese ecosystems (<xref ref-type="bibr" rid="B40">Xie et&#xa0;al., 2015</xref>), combined with the socio-economic development status, to modify the economic value created by the unit area grain yield in the study area, and obtain the unit area ecological service value of the Baiyangdian Basin, which represents the resistance of ecological space resilience. The calculation formula is:</p>
<disp-formula>
<label>(6)</label>
<mml:math display="block" id="M6">
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>=</mml:mo>
<mml:mi>E</mml:mi>
<mml:mi>S</mml:mi>
<mml:mi>V</mml:mi>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(7)</label>
<mml:math display="block" id="M7">
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>k</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<disp-formula>
<label>(8)</label>
<mml:math display="block" id="M8">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>a</mml:mi>
</mml:msub>
<mml:mo>=</mml:mo>
<mml:mfrac>
<mml:mn>1</mml:mn>
<mml:mn>7</mml:mn>
</mml:mfrac>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>m</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>p</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mi>M</mml:mi>
</mml:mfrac>
</mml:mrow>
</mml:mstyle>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>P</italic> is ecological space resistance, <italic>ESV</italic> is ecosystem service value, <italic>A<sub>i</sub>
</italic> is the area of land use type <italic>i</italic>, <italic>VC<sub>i</sub>
</italic> is the ecosystem service value per unit area of class <italic>i</italic> land use type, <italic>EC<sub>j</sub>
</italic> is the equivalent of the ecosystem service value of item <italic>j</italic> of certain land use type, and <italic>E<sub>a</sub>
</italic>is the economic value (yuan/hectare) of providing food production service function for farmland ecosystem per unit area. <italic>i</italic> is the crop type, <italic>p<sub>i</sub>
</italic> is the national average price of a crop in a given year (yuan/ton), <italic>q<sub>i</sub>
</italic> is the planted area of a crop (hectare), <italic>M</italic> is the planted area of all crops (hectare).</p>
</sec>
<sec id="s2_3_2_2">
<label>2.3.2.2</label>
<title>Adaptability</title>
<p>The more stable the ecological space is, the more adaptability (<italic>A</italic>) it is. In this study, the adaptability is expressed through the relevant indicators of landscape structure stability (<xref ref-type="bibr" rid="B30">Peterson, 2002</xref>). Among them, the stability of landscape structure is measured by the relevant landscape index (<xref ref-type="bibr" rid="B31">Turner, 1989</xref>). The calculation formula is:</p>
<disp-formula>
<label>(9)</label>
<mml:math display="block" id="M9">
<mml:mrow>
<mml:mi>A</mml:mi>
<mml:mo>=</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:mi>C</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>0.25</mml:mn>
<mml:mi>S</mml:mi>
<mml:mi>H</mml:mi>
<mml:mi>D</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>+</mml:mo>
<mml:mn>0.25</mml:mn>
<mml:mi>A</mml:mi>
<mml:mi>W</mml:mi>
<mml:mi>M</mml:mi>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>A</italic> is ecological space adaptability index; <italic>C</italic> is landscape fragmentation degree; <italic>SHDI</italic> is Shannon&#x2019;s diversity index (SHDI); <italic>AWM</italic> is the the area-weighted mean patch fractal dimension (AWMPFD).</p>
</sec>
<sec id="s2_3_2_3">
<label>2.3.2.3</label>
<title>Recoverability</title>
<p>The resilience of ecological space (<italic>R</italic>) reflects the ability and potential of ecological space to recover from damage. Land use plays an important role in the concept of ecosystem recoverability (<xref ref-type="bibr" rid="B10">Foster et&#xa0;al., 2003</xref>; <xref ref-type="bibr" rid="B4">Colding, 2007</xref>), so the recoverability is quantified by summing the area-weighted ecosystem resilience coefficients of all land use types. The calculation formula is based on the ecological resilience model and coefficient proposed by <xref ref-type="bibr" rid="B29">Peng et&#xa0;al. (2015)</xref>. The calculation formula is:</p>
<disp-formula>
<label>(10)</label>
<mml:math display="block" id="M10">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:mo>=</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mtext>&#xa0;</mml:mtext>
<mml:mo>&#x2211;</mml:mo>
<mml:mtext>&#xa0;</mml:mtext>
</mml:munderover>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mstyle>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>R</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>R</italic> is ecological space recoverability; <italic>A<sub>i</sub>
</italic> is the area ratio of land use type; <italic>RC<sub>i</sub>
</italic> is the ecosystem elasticity coefficient of land use <inline-formula>
<mml:math display="inline" id="im4">
<mml:mi>i</mml:mi>
</mml:math>
</inline-formula> type.</p>
<p>The calculations of resistance, adaptability and recoverability have different units, so when calculating resilience levels, it is necessary to standardize various indicators to the [0,1] range. In addition, to avoid amplification calculations when multiplying indicators, a root sign is needed to neutralize the order of magnitude. According to the natural breakpoint method, the resilience level is divided into high level, medium level and low level.</p>
<p>The specific calculation formula is as follows:</p>
<disp-formula>
<label>(11)</label>
<mml:math display="block" id="M11">
<mml:mrow>
<mml:mi>Re</mml:mi>
<mml:mi>s</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>l</mml:mi>
<mml:mi>i</mml:mi>
<mml:mi>e</mml:mi>
<mml:mi>n</mml:mi>
<mml:mi>c</mml:mi>
<mml:mi>e</mml:mi>
<mml:mo>=</mml:mo>
<mml:mroot>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mroot>
</mml:mrow>
</mml:math>
</disp-formula>
<p>Where, <italic>P</italic> is resistance, <italic>A</italic> is adaptability, <italic>R</italic> is recoverability.</p>
</sec>
</sec>
</sec>
</sec>
<sec id="s3" sec-type="results">
<label>3</label>
<title>Results</title>
<sec id="s3_1">
<label>3.1</label>
<title>Evaluation results of ecosystem service function importance</title>
<sec id="s3_1_1">
<label>3.1.1</label>
<title>The importance of water conservation function</title>
<p>According to the extracted average annual net primary productivity of vegetation, soil seepage factor, annual average precipitation factor and slope factor, formula (1) is used to calculate and obtain the evaluation result of the importance of water conservation function in Baiyangdian Basin, which is normalized in ArcGIS software. According to the natural breakpoint method, the results are divided into the areas of water conservation function is unimportant, generally important and extremely important (<xref ref-type="table" rid="T3"><bold>Table 3</bold></xref>) (<xref ref-type="fig" rid="f3"><bold>Figure 3</bold></xref>). The results show that the water conservation function in Baiyangdian Basin is affected by landform, temperature and precipitation, land use conditions, vegetation and other factors. Among them, 23.31% of the water conservation function in Baiyangdian Basin is unimportant, distributed in Zhuolu County and Yu County in the northern part of the study area. The terrain in this area is relatively high with a large slope, and an annual precipitation of less than 400 millimeters. Therefore, the importance of water conservation is relatively low. Generally important areas accounted for 53.57%, mainly distributed in flat terrain areas such as urban construction and development areas. There is a large amount of arable land with high density, and rainfall and vegetation coverage are at a moderate level in the entire study area. Therefore, the water conservation capacity is average. The extremely important area accounts for 23.12%, mainly distributed in the southern part of Laishui County, the southeastern part of Yi County, Dingxing County and other areas. The terrain in this area is relatively flat and open, so the water conservation capacity is better.</p>
<fig id="f3" position="float">
<label>Figure&#xa0;3</label>
<caption>
<p>Each Single Factor Importance Evaluation Results [<bold>(A)</bold> Water conservation; <bold>(B)</bold> Soil and water conservation; <bold>(C)</bold> Biodiversity conservation].</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g003.tif"/>
</fig>
<table-wrap id="T3" position="float">
<label>Table&#xa0;3</label>
<caption>
<p>Statistics of importance levels of each function.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" rowspan="2" align="center">degree of importance</th>
<th valign="middle" colspan="2" align="center">water conservation</th>
<th valign="middle" colspan="2" align="center">soil and water conservation</th>
<th valign="middle" colspan="2" align="center">biodiversity conservation</th>
</tr>
<tr>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">unimportant</td>
<td valign="middle" align="center">0.00-0.08</td>
<td valign="middle" align="center">23.31%</td>
<td valign="middle" align="center">0.00-0.04</td>
<td valign="middle" align="center">83.90%</td>
<td valign="middle" align="center">0.00-0.08</td>
<td valign="middle" align="center">19.17%</td>
</tr>
<tr>
<td valign="middle" align="center">generally important</td>
<td valign="middle" align="center">0.08-0.16</td>
<td valign="middle" align="center">53.57%</td>
<td valign="middle" align="center">0.04-0.13</td>
<td valign="middle" align="center">12.95%</td>
<td valign="middle" align="center">0.08-0.18</td>
<td valign="middle" align="center">49.60%</td>
</tr>
<tr>
<td valign="middle" align="center">extremely important</td>
<td valign="middle" align="center">0.16-0.36</td>
<td valign="middle" align="center">23.12%</td>
<td valign="middle" align="center">0.13-0.52</td>
<td valign="middle" align="center">3.14%</td>
<td valign="middle" align="center">0.18-0.42</td>
<td valign="middle" align="center">31.23%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3_1_2">
<label>3.1.2</label>
<title>The importance of soil and water conservation</title>
<p>According to the extracted rainfall erosivity factor, soil erodibility factor, topographic relief factor and vegetation cover factor, formula (2) is used for calculation, and the evaluation results of soil and water conservation function importance in Baiyangdian Basin are obtained, which are normalized in ArcGIS software. According to the method of natural segment point, the results are divided into the areas of unimportant, generally important and extremely important soil and water conservation function (<xref ref-type="table" rid="T3"><bold>Table 3</bold></xref>) (<xref ref-type="fig" rid="f3"><bold>Figure 3</bold></xref>). The results show that the spatial distribution of water and soil conservation importance index in Baiyangdian Basin show a trend of higher in northwest and lower in southeast. Among them, 83.90% of the areas in Baiyangdian Basin are unimportant in soil and water conservation function, mainly distributed in the plain area in the southeast of the study area. The terrain in this area is relatively flat, and most of them are cultivated land and construction land. The economy in this area is relatively developed, and there are relatively many artificial soil and water conservation measures. Therefore, there is less soil erosion and the importance of soil and water conservation functions is relatively low. The general important area accounts for 12.95%, while the extremely important area accounts for 3.14%. It is mainly distributed in the northwest of the study area. This area is mostly low mountains and hills, mainly composed of forest land and grassland, with high altitude, significant terrain fluctuations, and more precipitation in the area. Therefore, soil and water loss is serious.</p>
</sec>
<sec id="s3_1_3">
<label>3.1.3</label>
<title>The importance of biodiversity conservation functions</title>
<p>According to the extracted average annual net primary productivity of vegetation, average annual precipitation, average annual temperature and altitude factors, formula (3) is used to calculate and obtain the assessment result of biodiversity conservation function importance in Baiyangdian Basin, which is normalized in ArcGIS software. According to the method of natural segment point, the results are divided into regions with unimportant, generally important and extremely important biodiversity conservation functions (<xref ref-type="table" rid="T3"><bold>Table 3</bold></xref>) (<xref ref-type="fig" rid="f3"><bold>Figure 3</bold></xref>). The results show that the biodiversity conservation function in Baiyangdian Basin is mainly of general importance. Among them, 19.17% of the areas in the Baiyangdian Basin have unimportant biodiversity conservation functions, mainly distributed in the northern part of the study area. Due to the high terrain and significant changes in temperature during the day and night, the importance of biodiversity conservation functions is relatively low in this area. The generally important areas account for 49.60% of the total, mainly distributed in the southern part of the study area where there is more cultivated land and construction land. This area has a large population, urban and agricultural spaces continue to expand, and the self repair ability of the ecosystem is weak, so the biodiversity conservation function is poor. The extremely important area accounts for 31.23%, mainly distributed in the south of Laishui County, the southeast of Yi County and Dingxing County. The terrain in this area is mostly low mountains, hills and plains. Compared with the northern part of the study area, the terrain is lower, with abundant rainfall, abundant light and heat, and rich vegetation, so the importance of biodiversity conservation function is high.</p>
</sec>
<sec id="s3_1_4">
<label>3.1.4</label>
<title>Evaluation of the importance of ecosystem service functions</title>
<p>The evaluation results of a single function often can only reflect the impact results of a single type of factor, rather than being holistic. Therefore, ArcGIS is used to overlay each importance index to calculate the importance of ecosystem service functions. The results are normalized in the ArcGIS software, and the results are divided into regions with unimportant ecosystem service functions, generally important ecosystem service functions and extremely important ecosystem service functions according to the natural breakpoint method (<xref ref-type="table" rid="T4"><bold>Table 4</bold></xref>) (<xref ref-type="fig" rid="f4"><bold>Figure 4</bold></xref>).</p>
<table-wrap id="T4" position="float">
<label>Table&#xa0;4</label>
<caption>
<p>Statistics of importance levels of ecosystem service functions in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center"/>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">unimportant</td>
<td valign="middle" align="center">0.00-0.08</td>
<td valign="middle" align="center">19.14%</td>
</tr>
<tr>
<td valign="middle" align="center">generally important</td>
<td valign="middle" align="center">0.08-0.18</td>
<td valign="middle" align="center">49.13%</td>
</tr>
<tr>
<td valign="middle" align="center">extremely important</td>
<td valign="middle" align="center">0.18-0.42</td>
<td valign="middle" align="center">31.73%</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f4" position="float">
<label>Figure&#xa0;4</label>
<caption>
<p>Evaluation on the Importance of Ecosystem Service Function in Baiyangdian Basin.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g004.tif"/>
</fig>
<p>The results indicate that 19.14% of the areas with unimportant ecosystem service functions in Baiyangdian Basin are mainly distributed in the northern part of the study area, where the net primary productivity of vegetation is low and the altitude is high. The general important areas account for 49.13%, mainly distributed in transition areas from high to low terrain, as well as flat terrain areas such as urban development areas. The extremely important area accounts for 31.73%, mainly distributed in the central and southwestern parts of the basin. This region is low in elevation and flat in terrain, with more rivers flowing through it and better water storage capacity and water regulation capacity. And the land use type of this area is mostly cultivated land, with good cultivation conditions. In addition, the importance of water conservation function is higher in the eastern part of the study area, the importance of water and soil conservation is lower, and the importance of biodiversity conservation function is higher, so the ecosystem in this region is better.</p>
</sec>
</sec>
<sec id="s3_2">
<label>3.2</label>
<title>Assessment result of eco-environmental sensitivity</title>
<p>The ecologically sensitive areas of Baiyangdian basin are mainly distributed in the north of the basin with higher elevations and the river pooling area, and the overall pattern is non-sensitive &gt; medium sensitive &gt; highly sensitive (<xref ref-type="table" rid="T5"><bold>Table 5</bold></xref>) (<xref ref-type="fig" rid="f5"><bold>Figure 5</bold></xref>). Among them, the non-sensitive areas accounted for 54.63%, mainly concentrated in the southern plain, where the terrain was flat and open, with distinct four seasons and obvious continental climate characteristics. However, because human production and living activities would bring certain pressure to the ecological environment, there were also some medium-sensitive areas in this region. The medium sensitive area accounted for 23.41%, and the highly sensitive area accounts for 21.96%. They are mainly distributed in areas such as Zhuolu County, Laishui County, Laiyuan County, and Fuping County with higher elevations in the northwest. The average temperature in this area is low, the temperature difference between day and night is large, there are many mountains, basins and hills in some areas, complicated terrain, poor living environment and frequent geological disasters, so the ecological sensitivity is relatively high. In addition, Baiyangdian District has a high ecological sensitivity due to the vast water area where rivers gather, rich aquatic animal and plant resources, and high ecosystem diversity.</p>
<table-wrap id="T5" position="float">
<label>Table&#xa0;5</label>
<caption>
<p>Statistics of ecological and environmental sensitivity levels in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center"/>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">non-sensitive</td>
<td valign="middle" align="center">0.00-0.36</td>
<td valign="middle" align="center">54.63%</td>
</tr>
<tr>
<td valign="middle" align="center">medium sensitive</td>
<td valign="middle" align="center">0.36-0.60</td>
<td valign="middle" align="center">23.41%</td>
</tr>
<tr>
<td valign="middle" align="center">highly sensitive</td>
<td valign="middle" align="center">0.60-1.00</td>
<td valign="middle" align="center">21.96%</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f5" position="float">
<label>Figure&#xa0;5</label>
<caption>
<p>Ecological Sensitivity Assessment of Baiyangdian Basin.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g005.tif"/>
</fig>
</sec>
<sec id="s3_3">
<label>3.3</label>
<title>Identification of ecological space</title>
<p>By using ArcGIS to overlay the importance layer of ecosystem service functions and the sensitivity layer of ecological environment, the ecological space of the study area is comprehensively identified. According to the natural breakpoint method, the ecological space is classified into importance levels (<xref ref-type="table" rid="T6"><bold>Table 6</bold></xref>) (<xref ref-type="fig" rid="f6"><bold>Figure 6</bold></xref>).</p>
<table-wrap id="T6" position="float">
<label>Table&#xa0;6</label>
<caption>
<p>Ecological spatial identification results in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center"/>
<th valign="middle" align="center">threshold</th>
<th valign="middle" align="center">proportion</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">non-critical area</td>
<td valign="middle" align="center">0.00-0.35</td>
<td valign="middle" align="center">43.00%</td>
</tr>
<tr>
<td valign="middle" align="center">general area</td>
<td valign="middle" align="center">0.35-0.59</td>
<td valign="middle" align="center">30.78%</td>
</tr>
<tr>
<td valign="middle" align="center">core area</td>
<td valign="middle" align="center">0.59-1.00</td>
<td valign="middle" align="center">26.22%</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f6" position="float">
<label>Figure&#xa0;6</label>
<caption>
<p>Distribution of Ecological Space Importance in Baiyangdian Basin.</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g006.tif"/>
</fig>
<p>The results show that 43.00% of the ecological space in the Baiyangdian Basin is non-critical area, mainly distributed in the southern part of the Baiyangdian Basin and the northern part of Yu County. Among them, Yu County has a relatively high altitude and typical loess landform development, while the southern region of the study area is mainly composed of construction land and cultivated land, which is a gathering area for human activities with less vegetation coverage. Therefore, the importance of ecological space is relatively low. The general type area accounts for 30.78%, and its distribution is relatively uniform but fragmented. The core type accounts for 26.22%, mainly in Laishui County, Layuan County, Fuping County, Tang County, Yi County and other central mountainous areas as well as Baiyangdian District where rivers come together. Among them, the land use types in the central part of the basin are mainly woodland and grassland, which is the main distribution area of forest ecosystem in the basin. There are a number of forest parks and natural ecological scenic spots, the terrain includes mountains, hills, plains, river gullies, is the main ecological space protection zone; Baiyangdian District is the gathering place of several rivers, which is the largest lake in Hebei Province and an important ecological water body, so it has become the core area of ecological space in the basin.</p>
</sec>
<sec id="s3_4">
<label>3.4</label>
<title>Horizontal spatial distribution of resilience</title>
<p>According to the calculation results of formula (11), the resilience level is divided into high level, medium level and low level according to the natural breakpoint method in ArcGIS  (<xref ref-type="table" rid="T7"><bold>Table 7</bold></xref>) (<xref ref-type="fig" rid="f7"><bold>Figure 7</bold></xref>).</p>
<table-wrap id="T7" position="float">
<label>Table&#xa0;7</label>
<caption>
<p>Statistics of resilience level in Baiyangdian Basin.</p>
</caption>
<table frame="hsides">
<thead>
<tr>
<th valign="middle" align="center"/>
<th valign="middle" align="center">resilience level value</th>
<th valign="middle" align="center">proportion</th>
</tr>
</thead>
<tbody>
<tr>
<td valign="middle" align="center">low level</td>
<td valign="middle" align="center">0.00-0.06</td>
<td valign="middle" align="center">68.60%</td>
</tr>
<tr>
<td valign="middle" align="center">medium level</td>
<td valign="middle" align="center">0.06-0.19</td>
<td valign="middle" align="center">29.89%</td>
</tr>
<tr>
<td valign="middle" align="center">high level</td>
<td valign="middle" align="center">0.19-0.57</td>
<td valign="middle" align="center">1.51%</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="f7" position="float">
<label>Figure&#xa0;7</label>
<caption>
<p>Distribution of Ecological Space Resilience Level in Baiyangdian Basin [<bold>(A)</bold> the result of resistance; <bold>(B)</bold> the result of adaptability; <bold>(C)</bold> the result of recoverability; <bold>(D)</bold> the result of resilience].</p>
</caption>
<graphic mimetype="image" mime-subtype="tiff" xlink:href="fevo-11-1200218-g007.tif"/>
</fig>
<p>The results show that the resilience level of Baiyangdian Basin is not high, and the area with low resilience level accounts for 68.60%, mainly in the southern part of the study area and the northern part of Yu County. Among them, there are a lot of cultivated land and construction land in the southern part of the study area, which occupy a large amount of natural resources. In addition, the flat terrain and convenient transportation have led to the continuous improvement of the urbanization level in these areas, which has led to the increase of landscape fragmentation and the increase of human interference on the landscape, thus reducing the level of ecological resilience in these areas; The northern part of Yu County is high in altitude, and the temperature difference between day and night is large, so that its ability to resist external interference is reduced. The medium level areas account for 29.89%, mainly distributed in Laishui - Laiyuan - Fuping area. These areas have variable topography, including mountains, hills, plains, rivers and ravines, etc. The land use types are mainly woodland and grassland, and there are a large number of forest parks and ecological scenic spots. However, due to the high altitude of this area, insufficient rainfall and great disparity in geomorphology, the resilience level is not high and is in the medium grade. The proportion of areas with high resilience levels is only 1.51%, mainly in the Baiyangdian District. This area has a high diversity index of aroma, high landscape heterogeneity, and abundant species. Therefore, its ecosystem has a good self-regulation ability and a high resilience level.</p>
<p>In contrast, the resilience level of core and important ecological space regions is higher, while that of non-critical ecological space regions is lower. This shows that the more complex the structure of ecological space is, the more complete the function is, the better the ecological environment is, the higher its own stability is, the less susceptible it is to external interference, and the higher the resilience level is. On the other hand, when the ecological space is not critical, the production efficiency of the ecosystem is low, the ability to resist natural disasters and external interference is also low, and the resilience level of the ecosystem is low due to the simple species composition of the ecosystem.</p>
</sec>
</sec>
<sec id="s4">
<label>4</label>
<title>Conclusion and discussion</title>
<sec id="s4_1" sec-type="conclusions">
<label>4.1</label>
<title>Conclusion</title>
<p>In this study, the ecological space of Baiyangdian Basin was identified by combining the evaluation of the importance of ecosystem service function and the evaluation of eco-environmental sensitivity. Secondly, based on the principles of integrity, orientation and operability, resilience evaluation indexes were selected from three aspects: resistance, adaptability and recoverability to cope with crises. Finally, the ecological resilience level of the study area was analyzed. It lays a theoretical foundation for the optimization of ecological space pattern in Baiyangdian Basin. Through the above research, this paper draws the following conclusions:</p>
<list list-type="order">
<list-item>
<p>There is a significant difference in the importance of ecological space in the study area. In general, 26.22% of the study area belongs to the core type of ecological space, 30.78% to the general type of ecological space, and 43.00% to the non-critical type of ecological space. In addition, it can be clearly found through the research that the core areas of ecological space are concentrated in the areas with good vegetation, superior water conditions and high biodiversity. The non-critical areas are mostly plain areas, where the terrain is flat, mainly cultivated land and intensive construction land, or the non-critical areas are the northern Yu County with higher elevation, poor ecological environment and relatively less vegetation.</p>
</list-item>
<list-item>
<p>The resilience level in the study area is not high, the threshold of low resilience level is 0.06, the threshold of medium resilience level is 0.19, and the threshold of high resilience level is 0.57. From the perspective of spatial distribution, regions with low resilience level accounted for 68.60%, regions with medium resilience level accounted for 29.98%, and regions with high resilience level accounted for 1.51%. The low resilience level area is concentrated in the southern part of the basin where the city is rapidly expanding, the medium level area is concentrated in the central and northern part of the basin with more forest vegetation and relatively uneven terrain, and the high resilience level area is less concentrated in the river convergence area.</p>
</list-item>
</list>
</sec>
<sec id="s4_2" sec-type="discussion">
<label>4.2</label>
<title>Discussion</title>
<p>At present, the study of ecological spatial resilience is still in the exploratory stage, and most of the existing studies are on the theoretical level of ecological resilience. This study considered the structure and function of ecological space, delineated the importance level of ecological space in Baiyangdian Basin, combined with the ecological environment characteristics of the study area, quantitatively assessed the resilience level of ecological space from the perspective of evolutionary resilience, and established an evaluation index system of ecological space resilience.</p>
<p>In terms of the method of identifying ecological space, from the perspective of ecological needs, this study selects the factors that can directly reflect the function of ecological space. Combined with the importance analysis of ecosystem service function and the sensitivity analysis of ecological environment, the ecological space of Baiyangdian Basin was more comprehensively identified and classified into important types of ecological space. However, due to the different sources of basic data, there may be some errors in data fusion and calculation. In the future, it is necessary to further improve the evaluation accuracy when identifying ecological space.</p>
<p>In terms of the method of evaluating resilience, a more comprehensive evolutionary resilience concept was chosen from a research perspective. One is to view ecological space as a complex dynamic system with structure and function, with a focus on considering its ability to cope with uncertain risks. Secondly, from the perspective of evolutionary resilience, resilience is seen as the ability of ecological spaces to absorb external interference, reduce their own risks, and accelerate self recovery, in order to enhance the resistance, adaptability and recoverability of ecological spaces. In order to ensure the ecological security of Baiyangdian Basin, a different evaluation system was established according to the current situation of land use and local characteristics, and a resilience measurement model was established to measure the resilience level on a spatial scale. However, the results obtained by this research method are greatly influenced by the process of land use change, and there is inevitable uncertainty in the interpretation of remote sensing images. In addition, this study only considers ecological factors and neglects the impact of economy and society on resilience levels. In the future, comprehensive consideration should be given from multiple dimensions to improve the ecological space resilience indicator system.</p>
</sec>
</sec>
<sec id="s5" sec-type="data-availability">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec id="s6" sec-type="author-contributions">
<title>Author contributions</title>
<p>XG, XL, WS contributed to the conception of the study. XG, XL performed the experiment. XG, XL, ZG, MZ, WS contributed significantly to analysis and manuscript preparation. XL performed the data analyses and wrote the manuscript. All authors contributed to the article and approved the submitted version.</p>
</sec>
</body>
<back>
<sec id="s7" sec-type="funding-information">
<title>Funding</title>
<p>This research is funded by the Social Science Fund of Hebei Province (HB22GL030).</p>
</sec>
<sec id="s8" sec-type="COI-statement">
<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 id="s9" sec-type="disclaimer">
<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>Carpenter</surname> <given-names>S. R.</given-names>
</name>
<name>
<surname>Westley</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Turner</surname> <given-names>M. G.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Surrogates for resilience of social-ecological systems</article-title>. <source>Ecosystems</source> <volume>8</volume>, <fpage>941</fpage>&#x2013;<lpage>944</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10021-005-0170-y</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>M.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>E. D.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Measuring urban infrastructure resilience via pressure-state-response framework in four chinese municipalities</article-title>. <source>Appl. Sci.</source> <volume>12</volume>, <fpage>2819</fpage>&#x2013;<lpage>2819</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/APP12062819</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname> <given-names>Y. L.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>B. G.</given-names>
</name>
<name>
<surname>Zhong</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>L. Q.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>A. M.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Predictive simulation of ecological space based on a particle swarm optimization-Markov composite model: a case study for Chang Zhu-Tan urban agglomcrations</article-title>. <source>Acta Ecol. Sinica</source> <volume>38</volume>, <fpage>55</fpage>&#x2013;<lpage>64</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb201612082530</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Colding</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>&#x2018;Ecological land-use complementation&#x2019; for building resilience in urban ecosystems</article-title>. <source>Landscape Urban Planning</source> <volume>81</volume>, <fpage>46</fpage>&#x2013;<lpage>55</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.landurbplan.2006.10.016</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cui</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Sun</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Lv</surname> <given-names>C. D.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Construction and optimization of green space ecological networks in urban fringe areas: a case study with the urban fringe area of Tongzhou district in Beijing</article-title>. <source>J. Cleaner Prod.</source> <volume>276</volume>, <elocation-id>124266</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.jclepro.2020.124266</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Deng</surname> <given-names>H. D.</given-names>
</name>
<name>
<surname>Ye</surname> <given-names>L. F.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>L. X.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Han</surname> <given-names>Y. W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Review of urban ecological security research</article-title>. <source>J. Environ. Eng. Technol.</source> <volume>12</volume>, <fpage>248</fpage>&#x2013;<lpage>259</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.12153/j.issn.1674-991X.20210179</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Du</surname> <given-names>W. X.</given-names>
</name>
<name>
<surname>Shi</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>L. H.</given-names>
</name>
<name>
<surname>Zhai</surname> <given-names>G. F.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>W.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>Z. W.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Multidimensional measurements and analysis of urban economic resilience under the disturbance of economic risks: A case study of the Yangtze River Delta, China</article-title>. <source>Prog. Geography</source> <volume>41</volume>, <fpage>956</fpage>&#x2013;<lpage>971</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18306/d1kxiz.2022.06.002</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fabien</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Schulp</surname> <given-names>C. J. E.</given-names>
</name>
<name>
<surname>Verburg</surname> <given-names>P. H.</given-names>
</name>
<name>
<surname>Teeffelen</surname> <given-names>A. J. A. v.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Testing the applicability of ecosystem services mapping methods for peri-urban contexts: A case study for Paris</article-title>. <source>Ecol. Indicators</source> <volume>83</volume>, <fpage>504</fpage>&#x2013;<lpage>514</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.ecolind.2017.07.046</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Folke</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Resilience: The emergence of a perspective for social- ecological systems analyses</article-title>. <source>Global Environ. Change</source> <volume>16</volume>, <fpage>253</fpage>&#x2013;<lpage>267</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.gloenvcha.2006.04.002</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Foster</surname> <given-names>D.</given-names>
</name>
<name>
<surname>Swanson</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Aber</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Burke</surname> <given-names>I.</given-names>
</name>
<name>
<surname>Brokaw</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Tilman</surname> <given-names>D.</given-names>
</name>
<etal/>
</person-group>. (<year>2003</year>). <article-title>The importance of land-use legacies to ecology and conservation</article-title>. <source>BioScience</source> <volume>53</volume>, <fpage>77</fpage>&#x2013;<lpage>88</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1641/0006-3568(2003)053[0077:TIOLUL]2.0.CO;2</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname> <given-names>J. X.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>X. M.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>Z.</given-names>
</name>
<name>
<surname>Hou</surname> <given-names>P.</given-names>
</name>
<etal/>
</person-group>. (<year>2021</year>). <article-title>Evaluating changes in ecological land and effect of protecting important ecological spaces in China</article-title>. <source>J. Geographical Sci.</source> <volume>31</volume>, <fpage>1245</fpage>&#x2013;<lpage>1260</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.11821/dlxb202107010</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Haimes Yacov</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>On the definition of resilience in systems</article-title>. <source>Risk Analysis</source> <volume>29</volume>, <fpage>498</fpage>&#x2013;<lpage>501</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1111/j.1539-6924.2009.01216.x</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Holling</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>1973</year>). <article-title>Resilience and stability of ecological systems</article-title>. <source>Annu. Rev. Ecol. Systematics</source> <volume>4</volume>, <fpage>1</fpage>&#x2013;<lpage>23</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/ANNUREV.ES.04.110173.000245</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Holling</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>1996</year>). &#x201c;<article-title>Engineering resilience versus ecological resilience</article-title>,&#x201d; in <source>Engineering within ecological constraints</source> (<publisher-loc>Washington</publisher-loc>: <publisher-name>National academies press</publisher-name>), <fpage>31</fpage>&#x2013;<lpage>34</lpage>.</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Huang</surname> <given-names>X. Y.</given-names>
</name>
<name>
<surname>Zhao</surname> <given-names>X. M.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Jiang</surname> <given-names>Y. F.</given-names>
</name>
<name>
<surname>Lai</surname> <given-names>X. H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>The natural ecological spatial management zoning based on ecosystem service function and ecological sensitivity</article-title>. <source>Acta Ecol. Sinica</source> <volume>40</volume>, <fpage>1065</fpage>&#x2013;<lpage>1076</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb201811152474</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname> <given-names>M. Q.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Discussion on the concept definition and spatial boundary classification of &#x201c;Production-Living-Ecological&#x201d; space</article-title>. <source>Urban Dev. Stud.</source> <volume>27</volume>, <fpage>43</fpage>&#x2013;<lpage>48</lpage>. doi: CNKI:SUN:CSFY.0.2020-04-007
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname> <given-names>G.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>K.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Guo</surname> <given-names>B. S.</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Trade-offs in land-use competition and sustainable land development in the north China plain</article-title>. <source>Technol. Forecasting Soc. Change</source> <volume>141</volume>, <fpage>36</fpage>&#x2013;<lpage>46</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.techfore.2019.01.004</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Wei</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Lu</surname> <given-names>Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Construction of ecological security pattern based on the importance of ecosystem service functions and ecological sensitivity assessment: a case study in Fengxian County of Jiangsu Province, China</article-title>. <source>Environment Dev. Sustainability</source> <volume>23</volume>, <fpage>1</fpage>&#x2013;<lpage>28</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10668-020-00596-2</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname> <given-names>L. Q.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>Y. Q.</given-names>
</name>
<name>
<surname>Zheng</surname> <given-names>H.</given-names>
</name>
<name>
<surname>Xiao</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L.</given-names>
</name>
<etal/>
</person-group>. (<year>2019</year>). <article-title>A method for evaluating ecological space and ecological conservation redlines in river basins: A case of the Yangtze River Basin</article-title>. <source>Acta Ecol. Sinica</source> <volume>39</volume>, <fpage>835</fpage>&#x2013;<lpage>843</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb201802060314</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>H. J.</given-names>
</name>
<name>
<surname>Niu</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Yu</surname> <given-names>Q.</given-names>
</name>
<name>
<surname>Yang</surname> <given-names>L. Z.</given-names>
</name>
<name>
<surname>Ma</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Qiu</surname> <given-names>S.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Evaluation of the spatiotemporal evolution of China&#x2019;s ecological spatial network function-structure and its pattern optimization</article-title>. <source>Remote Sensing</source> <volume>14</volume>, <fpage>4593</fpage>&#x2013;<lpage>4593</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3390/RS14184593</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>H. F.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Value evaluation of biodiversity conservation in Dongting lake wetland</article-title>. <source>J. Cent. South Univ. Forestry Technol.</source> <volume>41</volume>, <fpage>140</fpage>&#x2013;<lpage>147</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.14067/j.cnki.1673-923x.2021.10.016</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>J. Y.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>W. C.</given-names>
</name>
<name>
<surname>Kuang</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Ning</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Remote sensing-based analysis of the spatiotemporal characteristics of built-up area across China based on the Plan for Major Function-oriented Zones</article-title>. <source>Acta Geographica Sinica</source> <volume>71</volume>, <fpage>355</fpage>&#x2013;<lpage>369</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.11821/dlxb201603001</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname> <given-names>X. Q.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>K. L.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Simulation and prediction of multi-scenario evolution of ecological space based on FLUS model: A case study of the Yangtze River Economic Belt, China</article-title>. <source>J. Geographical Sci.</source> <volume>33</volume>, <fpage>373</fpage>&#x2013;<lpage>391</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/S11442-023-2087-9</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Long</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y. Q.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>T. T.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>A. X.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>A primary study on ecological land use classification</article-title>. <source>Ecol. Environ. Sci.</source> <volume>24</volume>, <fpage>1</fpage>&#x2013;<lpage>7</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.16258/j.cnki.1674-5906.2015.01.001</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname> <given-names>F. H.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y. X.</given-names>
</name>
<name>
<surname>Peng</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>J. S.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Assessing urban landscape ecological risk through an adaptive cycle framework</article-title>. <source>Landscape Urban Planning</source> <volume>180</volume>, <fpage>125</fpage>&#x2013;<lpage>134</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.landurbplan.2018.08.014</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ngom</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Gosselin</surname> <given-names>P.</given-names>
</name>
<name>
<surname>Blais</surname> <given-names>C.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Reduction of disparities in access to green spaces: Their geographic insertion and recreational functions matter</article-title>. <source>Appl. Geography</source> <volume>66</volume>, <fpage>35</fpage>&#x2013;<lpage>51</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.apgeog.2015.11.008</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ouyang</surname> <given-names>Z. Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>X. K.</given-names>
</name>
<name>
<surname>Miao</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>China&#x2019;s eco-environmental sensitivity and its spatial heterogeneity</article-title>. <source>Acta Ecol. Sinica</source> <volume>20</volume>, <fpage>9</fpage>&#x2013;<lpage>12</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3321/j.issn:1000-0933.2000.01.002</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Pan</surname> <given-names>F. J.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H. Z.</given-names>
</name>
<name>
<surname>Song</surname> <given-names>M. J.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>A.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L. Y.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Research on spatial identification of ecological space in national key ecological function regions based on GIS: a case study of Changyang County in Hubei Province, China</article-title>. <source>J. Cent. China Normal Univ. (Nat. Sci.)</source> <volume>54</volume>, <fpage>658</fpage>&#x2013;<lpage>669</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.19603/j.cnki.1000-1190.2020.04.017</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname> <given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>Y. X.</given-names>
</name>
<name>
<surname>Wu</surname> <given-names>J. S.</given-names>
</name>
<name>
<surname>Lv</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Hu</surname> <given-names>X. X.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Linking ecosystem services and landscape patterns to assess urban ecosystem health: A case study in Shenzhen City, China</article-title>. <source>Landscape Urban Planning</source> <volume>143</volume>, <fpage>56</fpage>&#x2013;<lpage>68</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1016/j.landurbplan.2015.06.007</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peterson</surname> <given-names>G. D.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Contagious disturbance, ecological memory, and the emergence of landscape pattern</article-title>. <source>Ecosystems</source> <volume>5</volume>, <fpage>329</fpage>&#x2013;<lpage>338</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s10021-001-0077-1</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Turner</surname> <given-names>M. G.</given-names>
</name>
</person-group> (<year>1989</year>). <article-title>Landscape ecology: The effect of pattern on process</article-title>. <source>Annu. Rev. Ecol. Systematics</source> <volume>20</volume>, <fpage>171</fpage>&#x2013;<lpage>197</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1146/annurev.es.20.110189.001131</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Vugrin</surname> <given-names>E. D.</given-names>
</name>
<name>
<surname>Warren</surname> <given-names>D. E.</given-names>
</name>
<name>
<surname>Ehlen</surname> <given-names>M. A.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>A resilience assessment framework for infrastructure and economic systems: Quantitative and qualitative resilience analysis of petrochemical supply chains to a hurricane</article-title>. <source>Process Saf. Progress</source> <volume>30</volume>, <fpage>280</fpage>&#x2013;<lpage>290</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1002/prs.10437</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Walker</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Holling</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Carpenter</surname> <given-names>S. R.</given-names>
</name>
<name>
<surname>Kinzig</surname> <given-names>A.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Resilience, adaptability and transformability in social-ecological systems</article-title>. <source>Ecol. Society</source> <volume>9</volume>, <elocation-id>5</elocation-id>. doi:&#xa0;<pub-id pub-id-type="doi">10.5751/ES-00650-090205</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="book">
<person-group person-group-type="author">
<name>
<surname>Walker</surname> <given-names>B.</given-names>
</name>
<name>
<surname>Salt</surname> <given-names>D.</given-names>
</name>
</person-group> (<year>2006</year>). <source>Resilience thinking: sustaining ecosystem and people in a changing world</source> (<publisher-loc>Washington, DC</publisher-loc>: <publisher-name>Island Press</publisher-name>).</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>S. H.</given-names>
</name>
<name>
<surname>Huang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>X. L.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>J. H.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Spatio-temporal variations in ecological spaces and their ecological carrying status in China&#x2019;s mega-urban agglomerations</article-title>. <source>J. Geographical Sci.</source> <volume>32</volume>, <fpage>1683</fpage>&#x2013;<lpage>1704</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.1007/s11442-022-2018-1</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname> <given-names>F. Y.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>K. Y.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>T.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>P.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Progress and prospect of research on urban ecological space</article-title>. <source>Prog. Geography</source> <volume>36</volume>, <fpage>207</fpage>&#x2013;<lpage>218</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18306/dlkxjz.2017.02.007</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Xiu</surname> <given-names>C. L.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Study on the concept and analytical framework of city network resilience</article-title>. <source>Prog. Geography</source> <volume>39</volume>, <fpage>488</fpage>&#x2013;<lpage>502</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.18306/dlkxiz.2020.03.013</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xia</surname> <given-names>C. Y.</given-names>
</name>
<name>
<surname>Dong</surname> <given-names>Z. Y. Z.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>B.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Spatio-temporal analysis and simulation of urban ecological resilience: A case study of Hangzhou</article-title>. <source>Acta Ecol. Sinica</source> <volume>42</volume>, <fpage>116</fpage>&#x2013;<lpage>126</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb202101030007</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>G. D.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Biodiversity and ecosystem services in city</article-title>. <source>Environ. Protection</source> <volume>43</volume>, <fpage>25</fpage>&#x2013;<lpage>28</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.14026/j.cnki.0253-9705.2015.05.004</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname> <given-names>G. D.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>C. X.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S. M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Improvement of the evaluation method for ecosystem service value based on per unit area</article-title>. <source>J. Natural Resources</source> <volume>30</volume>, <fpage>1243</fpage>&#x2013;<lpage>1254</lpage>. doi: <pub-id pub-id-type="doi">10.11849/zrzyxb.2015.08.001</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname> <given-names>X. P.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>L. K.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>Y. B.</given-names>
</name>
<name>
<surname>Hou</surname> <given-names>Y. J.</given-names>
</name>
<name>
<surname>Niu</surname> <given-names>J.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Review and prospects of resilient city theory</article-title>. <source>Geogr. Geo-Information Sci.</source> <volume>37</volume>, <fpage>78</fpage>&#x2013;<lpage>84</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.3969/j.issn.1672-0504.2021.06.012</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname> <given-names>F.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>X. B.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L. J.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>W. H.</given-names>
</name>
<name>
<surname>Fu</surname> <given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname> <given-names>H.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Study of ecological land in China: Conception, classification, and spatial-temporal pattern</article-title>. <source>Acta Ecol. Sinica</source> <volume>35</volume>, <fpage>4931</fpage>&#x2013;<lpage>4943</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb201311052672</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhan</surname> <given-names>M. S.</given-names>
</name>
<name>
<surname>Zhu</surname> <given-names>J. H.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Spatial planning of Qingyun Mountain scenic spot in Benxi, Northeast China based on ecological sensitivity assessment</article-title>. <source>Chin. J. Appl. Ecol.</source> <volume>30</volume>, <fpage>2352</fpage>&#x2013;<lpage>2360</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.13287/j.1001-9332.201907.037</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>C. S.</given-names>
</name>
<name>
<surname>Fan</surname> <given-names>N.</given-names>
</name>
<name>
<surname>Liu</surname> <given-names>C. L.</given-names>
</name>
<name>
<surname>Xie</surname> <given-names>G. D.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Spatio-temporal pattern and evolution of ecosystem water conservation in China from 1990 to 2018</article-title>. <source>Acta Ecol. Sinica</source> <volume>43</volume>, <fpage>1</fpage>&#x2013;<lpage>10</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb202107191937</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>L.</given-names>
</name>
<name>
<surname>Yue</surname> <given-names>W. Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Integrated recognition of urban ecological space: A case study on Hangzhou</article-title>. <source>Acta Ecol. Sinica</source> <volume>39</volume>, <fpage>6460</fpage>&#x2013;<lpage>6468</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5846/stxb201804100817</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname> <given-names>C. X.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>L. M.</given-names>
</name>
<name>
<surname>Li</surname> <given-names>S. M.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>C. S.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Soil conservation of national key ecological function areas</article-title>. <source>J. Resour. Ecol.</source> <volume>6</volume>, <fpage>397</fpage>&#x2013;<lpage>404</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.5814/j.issn.1674-764x.2015.06.007</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname> <given-names>H. L.</given-names>
</name>
<name>
<surname>Xu</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname> <given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname> <given-names>G. L.</given-names>
</name>
<name>
<surname>Chen</surname> <given-names>X.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Quantitative identification and analysis of spatiotemporal evolution characteristics of ecological space: A case study of Hulun Buir city</article-title>. <source>J. Natural Resources</source> <volume>37</volume>, <fpage>2153</fpage>&#x2013;<lpage>2168</lpage>. doi:&#xa0;<pub-id pub-id-type="doi">10.31497/zrzyxb.20220815</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname> <given-names>L. M.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>From social vulnerability to social ecological resilience: paradigm transformation of disaster social science research</article-title>. <source>Thinking</source> <volume>41</volume>, <fpage>50</fpage>&#x2013;<lpage>57</lpage>.</citation>
</ref>
</ref-list>
</back>
</article>