<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="research-article" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Environ. Sci.</journal-id>
<journal-title>Frontiers in Environmental Science</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Environ. Sci.</abbrev-journal-title>
<issn pub-type="epub">2296-665X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1488439</article-id>
<article-id pub-id-type="doi">10.3389/fenvs.2024.1488439</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Environmental Science</subject>
<subj-group>
<subject>Original Research</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Landscape pattern evolution and ecological security assessment based on land use changes in mining cities: a case study of Heihe city</article-title>
<alt-title alt-title-type="left-running-head">Shi et al.</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fenvs.2024.1488439">10.3389/fenvs.2024.1488439</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Shi</surname>
<given-names>Yu</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2774012/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Fan</surname>
<given-names>Xiaoxiao</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="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2787598/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Chen</surname>
<given-names>Haonan</given-names>
</name>
<xref ref-type="aff" rid="aff3">
<sup>3</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2742572/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/software/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Kao</surname>
<given-names>Yuhan</given-names>
</name>
<xref ref-type="aff" rid="aff4">
<sup>4</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2891255/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Sun</surname>
<given-names>Meiqi</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/2891262/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/data-curation/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Architecture and Art Design</institution>, <institution>Hebei University of Technology</institution>, <addr-line>Tianjin</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>School of Architecture and Art Design</institution>, <institution>Hebei University of Technology</institution>, <institution>Key Laboratory of Healthy Human Settlements in Hebei Province</institution>, <addr-line>Tianjin</addr-line>, <country>China</country>
</aff>
<aff id="aff3">
<sup>3</sup>
<institution>School of Economics</institution>, <institution>Capital University of Economics and Business</institution>, <addr-line>Beijing</addr-line>, <country>China</country>
</aff>
<aff id="aff4">
<sup>4</sup>
<institution>Jiangxi University of Finance and Economics</institution>, <institution>School of Design and Art</institution>, <addr-line>Nanchang</addr-line>, <country>China</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2336688/overview">S&#xe9;rgio Ant&#xf3;nio Neves Lousada</ext-link>, University of Madeira, Portugal</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2709510/overview">Duanqiang Zhai</ext-link>, Tongji University, China</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2855828/overview">Biao Liu</ext-link>, China University of Geosciences, China</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Xiaoxiao Fan, <email>2013116@hebut.edu.cn</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>18</day>
<month>11</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1488439</elocation-id>
<history>
<date date-type="received">
<day>30</day>
<month>08</month>
<year>2024</year>
</date>
<date date-type="accepted">
<day>31</day>
<month>10</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Shi, Fan, Chen, Kao and Sun.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Shi, Fan, Chen, Kao and Sun</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 development of mining cities has led to changes in land use and the evolution of landscape patterns. Constructing an ecological security evaluation system can reflect the ecological security status of mining city areas and provide planning references for these cities. This study, based on Heihe City&#x2019;s land use data from 1980 to 2020, systematically analyzes the spatial and temporal evolution characteristics of land use, landscape patterns, and landscape ecological security levels by constructing a land use transfer matrix and calculating landscape pattern indices and landscape ecological security indices. The results show that: 1) Forest land is the main type of land use in Heihe City, accounting for over 50% of the total area. Land use changes primarily occurred between 2000&#x2013;2010 and 2010&#x2013;2020, with the spatial pattern characterized by overall stability and localized dramatic shifts, mainly involving the conversion of forest land to farmland and unused land. 2) From 1980 to 2020, the landscape ecological security pattern in Heihe City improved. Landscape diversity and landscape contagion increased, while landscape fragmentation, the largest patch area, and the average patch area decreased. Land use was optimized overall, but the trend of forest fragmentation became noticeable. 3) Between 1980 and 2020, the landscape ecological security level in Heihe City improved significantly. Driven by ecological restoration policies, the ecological security level in the southern region greatly improved, while the spatial pattern continued to show a trend of lower security in the south and higher security in the north. Specifically, the area proportions of low-security and relatively low-security areas increased by 6.23% and 9.55%, respectively. The spatial clustering of landscape ecological security levels is evident, with high-high value clustering mainly in the north and low-low value clustering mainly in the south. It is necessary to further promote ecological protection in the north to ensure the ecological barrier function, while strengthening ecological restoration in the south to improve ecological security levels. Additionally, continuous macro-policy regulation is needed to maintain long-term ecological security in Heihe City. The ecological security level of mining city landscapes is influenced by economic dependence, policy constraints, and environmental issues. This study can provide guidance for planning improvements in the study area and other mining cities.</p>
</abstract>
<kwd-group>
<kwd>ecological security assessment</kwd>
<kwd>Heihe city</kwd>
<kwd>land use</kwd>
<kwd>landscape pattern</kwd>
<kwd>mining city</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Land Use Dynamics</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<sec id="s1-1">
<title>1.1 Background</title>
<p>With the rise of mining cities, China introduced relevant policies to align with global standards. Since 1978, the global coal industry has gone through four stages: the first stage, from 1978 to 1992, was a period of deepening reform in the coal industry, during which developed countries began comprehensive adjustments and restructuring, promoting the construction of township coal mines and achieving rapid increases in coal production. During this time, a large number of mining cities emerged, while developing countries were still in the preparatory phase of resource development (<xref ref-type="bibr" rid="B18">Li et al., 2002</xref>). The second stage, from 1993 to 2001, was marked by coal industry reform aimed at poverty alleviation. In developing countries like China, the number of small coal mines decreased from over 80,000 in 1998 to around 22,000&#xa0;b&#xa0;y 2001. During this period, industrial and mining land use in China significantly decreased, leading to notable changes in land use patterns, while developed countries were experiencing the growth phase of mining city development (<xref ref-type="bibr" rid="B6">Fan et al., 2005</xref>). The third stage, from 2002 to 2011, focused on the sustainable and healthy development of the coal industry, with efforts to build a new coal industry system and transform the economic development model (<xref ref-type="bibr" rid="B56">Zhou and Tang, 2004</xref>). During this time, awareness grew about the limitations resource depletion posed to urban development, leading to explorations of sustainable development strategies for mining cities. As global coal production reached its designed capacity, the comprehensive regional development of mining cities improved, and their scale expanded. The fourth stage, from 2012 to the present, has been characterized by a focus on high-quality development in the coal industry, with policymakers worldwide emphasizing the ecological security of mining cities. Over these 4&#xa0;decades, the planning and development of mining cities have become a key focus of research in global land-use and spatial planning (<xref ref-type="bibr" rid="B22">Liu Y. et al., 2023</xref>).</p>
</sec>
<sec id="s1-2">
<title>1.2 Related work</title>
<p>Mining cities refer to cities or regions where the mining industry is the dominant economic driver (<xref ref-type="bibr" rid="B21">Liu X. et al., 2023</xref>). Mining activities have significant impacts on the environment, leading to ecological issues such as large-scale land subsidence and a decline in biodiversity due to long-term exploitation (<xref ref-type="bibr" rid="B31">Shi and Xiang, 2023</xref>). These environmental problems result in land-use changes that pose challenges to the sustainable development of these cities (<xref ref-type="bibr" rid="B27">Mardonova and Han, 2023</xref>). However, there are several gaps in the current research on mining cities: First, in terms of research content, scholars primarily focus on ecological restoration strategies in mining cities, particularly on water ecosystem management and vegetation restoration techniques, but there is a lack of systematic analysis of the evolution of landscape ecological security patterns in the context of land-use and spatial planning (<xref ref-type="bibr" rid="B8">Feng et al., 2023</xref>; <xref ref-type="bibr" rid="B17">Li Q. et al., 2023</xref>; <xref ref-type="bibr" rid="B41">Xu et al., 2023</xref>; <xref ref-type="bibr" rid="B54">Zhao et al., 2023</xref>). Second, from the research perspective, insufficient attention has been given to mining cities as resource-depleted urban areas, especially those with rich mining resources and significant ecological importance, such as Heihe City (<xref ref-type="bibr" rid="B36">Wei et al., 2023</xref>)., which remains under-researched. Exploring the long-term evolution of landscape ecological security patterns in mining cities can deepen our understanding of the ecological security dynamics of these cities and advance knowledge of the comprehensive management of ecological degradation. It also provides a scientific basis for assessing the effectiveness of ecological restoration efforts. Therefore, conducting long-term research on the landscape ecological security patterns of mining cities is essential for guiding ecological protection and restoration efforts.</p>
<p>Land-use change refers to alterations in the types, functions, or characteristics of land use within a specific region or area. Such changes directly influence the structure and patterns of landscapes (<xref ref-type="bibr" rid="B48">Zhang J. et al., 2023</xref>). The landscape pattern index method is commonly used to assess and describe the characteristics of landscape patterns in a given area or region. Previous research has often employed this method for quantitative studies on the landscape pattern changes of urban land-use types (<xref ref-type="bibr" rid="B46">Ye et al., 2022</xref>). However, in studies focusing on mining cities, which are resource-depleted urban areas, most research leans towards qualitative planning, and there is relatively little quantitative research using landscape indices (<xref ref-type="bibr" rid="B5">Dong et al., 2023</xref>; <xref ref-type="bibr" rid="B10">Gao et al., 2024</xref>). Li Baojie and others conducted a landscape pattern analysis of mining areas using the optimal analysis grain size, revealing that changes in regional land-use landscape patterns are the result of the combined effects of natural and human factors. Their analysis of the spatiotemporal evolution of landscape patterns in the Jiawang mining area based on the optimal grain size provides a new method for the quantitative assessment of land-use rationality and scientific validity (<xref ref-type="bibr" rid="B16">Li et al., 2016</xref>) Although some scholars have studied landscape patterns in mining cities, most of these studies focus on short-term dynamic landscape patterns over a few specific years, without addressing the dynamic changes over a period of nearly 40&#xa0;years (<xref ref-type="bibr" rid="B5">Dong et al., 2023</xref>; <xref ref-type="bibr" rid="B25">Luo et al., 2023</xref>; <xref ref-type="bibr" rid="B37">Wu et al., 2014</xref>; <xref ref-type="bibr" rid="B38">Wu et al., 2021</xref>). This paper selects a time frame that covers key periods of land policy reform and mining culture development in China, making the study of landscape pattern evolution in mining cities within this context more representative.</p>
<p>Ecological security assessment is the foundation and premise of ecological security research. It is the process of evaluating the health status of a specific region or ecosystem and its potential to withstand external impacts (<xref ref-type="bibr" rid="B26">Lv et al., 2024</xref>). This process reflects the degree to which human beings are protected from ecological damage and environmental pollution, playing a crucial role in maintaining the ecological balance of the Earth and supporting the development of mining cities. The ecological security of mining cities refers to the ability of their natural, artificial, and complex ecosystems to maintain their normal structure and function while also recovering from external pressures and stresses. As cities driven by mining activities, the ecosystems of mining cities are more fragile and complex than those of traditional cities. Currently, research on the ecological security of mining cities focuses more on establishing evaluation mechanism models and selecting evaluation indicators, with little attention paid to the relationship between ecological security and land-use change. This gap makes it difficult to guide the management and regulation of ecological security effectively (<xref ref-type="bibr" rid="B20">Li Z. et al., 2023</xref>; <xref ref-type="bibr" rid="B44">Yang and Ye, 2023</xref>; <xref ref-type="bibr" rid="B45">Yao et al., 2023</xref>; <xref ref-type="bibr" rid="B55">Zhao et al., 2024</xref>). By conducting ecological security assessments of mining cities, the impact of different land-use practices on the ecosystem can be evaluated, thus informing land-use planning and management to achieve sustainable utilization of land resources and protection of the ecological environment in mining cities (<xref ref-type="bibr" rid="B1">Chen and Ma, 2023</xref>; <xref ref-type="bibr" rid="B19">Li X. et al., 2023</xref>). This study combines the analysis of landscape pattern evolution with ecological security assessments in mining cities based on land-use change. The methodology primarily uses quantitative analysis, supplemented by qualitative analysis, to make the results of the ecological security assessment more convincing.</p>
</sec>
<sec id="s1-3">
<title>1.3 Objectives and contributions</title>
<p>Based on the above, this study selected Heihe City in Heilongjiang Province, China, as the research subject through data pre-selection, as it provided the most representative data results (<xref ref-type="bibr" rid="B43">Yang et al., 2022</xref>). Heihe City has a relatively low level of mining development and is still in its early stages, but it is rich in various mineral resources. As a typical mining city in northern China and a forestry city, Heihe boasts abundant forest resources and serves as a natural ecological barrier for the &#x201c;granary of Northeast China,&#x201d; giving it significant ecological importance. During its development, Heihe faces common challenges of resource-dependent cities, such as the over-exploitation of natural resources and ecological degradation. This study used land-use data from five periods&#x2014;1980, 1990, 2000, 2010, and 2020&#x2014;to conduct a detailed assessment of the landscape pattern evolution and ecological security of Heihe from 1980 to 2020. Methods employed include land-use transition matrices, landscape pattern indices, landscape ecological security assessments, and spatial autocorrelation analysis. The research aims to answer two specific questions.<list list-type="simple">
<list-item>
<p>(1) How did the landscape patterns of Heihe change across different spatial and temporal dimensions during the 4&#xa0;decades of policy turbulence from 1980 to 2020?</p>
</list-item>
<list-item>
<p>(2) How did these changes in landscape patterns affect the ecological security of the mining city across different spatial and temporal dimensions from 1980 to 2020?</p>
</list-item>
</list>
</p>
<p>Addressing these questions will provide policy insights for improving the ecological security of similar mining cities in China and offer a reference for studying the ecological security of mining cities in other countries and regions.</p>
</sec>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Study area overview</title>
<p>Heihe City is located in the northern part of Heilongjiang Province (47&#xb0;42&#x2032;&#x2013;51&#xb0;03&#x2032;N, 124&#xb0;45&#x2032;&#x2013;129&#xb0;18&#x2032;E), at the eastern end of the Greater Khingan Range and the northern part of the Lesser Khingan Range, with a total area of 6,872,600&#xa0;hm<sup>2</sup> (<xref ref-type="fig" rid="F1">Figure 1</xref>). The terrain within the city features rolling mountains and valleys, with higher elevations in the northwest gradually sloping down towards the southeast, forming a landscape that runs from the northwest to the southeast. Heihe has a cold temperate continental monsoon climate, characterized by high temperatures and strong winds in spring, simultaneous heat and rain in summer, sharp temperature drops in autumn, and cold, dry winters. The seasons are distinct, with long winters and short summers. According to the results of Heihe City&#x2019;s third national land survey, the current land-use types in Heihe include 2.1879 million hm<sup>2</sup> of arable land, 10,300&#xa0;hm<sup>2</sup> of orchards, 3.2125 million hm<sup>2</sup> of forest land, 160,800&#xa0;hm<sup>2</sup> of grassland, 858,400&#xa0;hm<sup>2</sup> of wetlands, 75,300&#xa0;hm<sup>2</sup> of urban, village, and industrial land, 51,900&#xa0;hm<sup>2</sup> of transportation land, and 116,300&#xa0;hm<sup>2</sup> of water areas and water conservancy facilities.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Location map of the study in Heihe City, Heilongjiang Province.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g001.tif"/>
</fig>
<p>Heihe City is located in the key metallogenic belt of the Greater and Lesser Khingan Ranges, with extremely rich mineral resources and vast development potential. A total of 95 types of minerals have been discovered within the city&#x2019;s jurisdiction, accounting for 68.35% of the 139 types found in the entire province. There are verified reserves of 38 types, which make up 40.86% of the 93 types identified in the province. These include one type of energy mineral, 20 types of metal minerals, and 17 types of non-metal minerals. Heihe holds the province&#x2019;s largest reserves for 26 types of minerals, including copper, gold, lead, zinc, tungsten, perlite, pumice, and shale used for cement production, while coal ranks sixth. Additionally, the city has over 140 million tons of marble reserves. The primary minerals being developed and utilized are coal, iron, copper, lead, zinc, molybdenum, gold, silver, and mineral water. However, Heihe City is currently facing several challenges. Some mines are approaching closure due to resource depletion, large-scale mine development is insufficient, and the costs of land and forest requisition are high, all of which have slowed the pace of mining development. As Heihe&#x2019;s mineral resources gradually become exhausted, problems such as environmental degradation and difficulties in sustaining economic development have become increasingly prominent (<xref ref-type="bibr" rid="B36">Wei et al., 2023</xref>).</p>
</sec>
<sec id="s2-2">
<title>2.2 Research methodology</title>
<sec id="s2-2-1">
<title>2.2.1 Data source and preprocessing</title>
<p>The land use data in this article is derived from the national land use type remote sensing monitoring spatial distribution data provided by the Resource and Environment Science and Data Center of the Chinese Academy of Sciences (<ext-link ext-link-type="uri" xlink:href="http://www.resdc.cn/">http://www.resdc.cn</ext-link>) (accessed on 3 February 2024), with a spatial resolution of 30 m &#xd7; 30&#xa0;m. The unified spatial reference coordinate system is the Krasovsky_1940_Albers coordinate system, including land use data for five periods: 1980, 1990, 2000, 2010, and 2020. The comprehensive accuracy of land use type evaluation in this dataset exceeds 90%.</p>
<p>Based on the &#x201c;Current land use classification&#x201d; (GB/T 21,010&#x2013;2017) and considering the landscape characteristics of Heihe City, the landscape types in Heihe City are classified into seven categories: arable land, forest land, grassland, water bodies, urban and rural residential land, industrial and mining land, and unused land. Urban and rural residential land includes commercial and service land, industrial and mining storage land, residential land, public management and public service land, special land, and transportation and transportation land. Considering the spatial resolution of the image data and the precision of information extraction, other construction land in urban and rural areas, industrial and mining land, and residential land, including factory areas, large industrial areas, quarries, etc., are uniformly classified as industrial and mining land.</p>
</sec>
<sec id="s2-2-2">
<title>2.2.2 Landscape pattern indices</title>
<p>Landscape pattern indices can highly condense landscape spatial pattern information, reflecting its structural composition and spatial configuration characteristics (<xref ref-type="bibr" rid="B11">Han et al., 2018</xref>; <xref ref-type="bibr" rid="B15">Jingwei and Haize, 2006</xref>; <xref ref-type="bibr" rid="B24">Lu et al., 2012</xref>). Based on principles of landscape ecology, nine landscape indices were selected for calculation: at the patch scale level, the indices include LPI, AREA_MN, FRAC_AM, and COHESION; at the landscape scale, the indices include N.P., CONTAG, LPI, FRAC_AM, and SHDI (<xref ref-type="table" rid="T1">Table 1</xref>).</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Landscape pattern index description.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<italic>Scale</italic>
</th>
<th align="center">
<italic>Landscape Pattern Index</italic>
</th>
<th align="center">
<italic>Index Significance</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="4" align="center">Patch Scale</td>
<td align="center">LPI</td>
<td align="left">- The proportion of the largest single patch area in the landscape is typically expressed as a percentage. The higher the value of LPI, the larger the patches of that type in the entire landscape, indicating a more significant impact on the landscape (<xref ref-type="bibr" rid="B3">Corry, 2005</xref>)</td>
</tr>
<tr>
<td align="center">AREA_MN</td>
<td align="left">- AREA_MN describes the average area of all landscape patches within a specific area. A higher AREA_MN index indicates that the average size of landscape patches is larger, resulting in lower landscape fragmentation (<xref ref-type="bibr" rid="B11">Han et al., 2018</xref>)</td>
</tr>
<tr>
<td align="center">FRAC_ AM</td>
<td align="left">- FRAC_AM is used to measure the complexity of landscape patch edges. When FRAC_AM approaches 1, it indicates that the shape of the patch is relatively simple, with smooth edges, such as circular or square patches. When FRAC_AM approaches 2, it indicates that the shape of the patch is very complex and convoluted, with long and irregular edges (<xref ref-type="bibr" rid="B15">Jingwei and Haize, 2006</xref>)</td>
</tr>
<tr>
<td align="center">COHESION</td>
<td align="left">- The COHESION index reflects the degree of spatial connectivity between individual patches in the landscape. A higher COHESION index indicates that the patches in the landscape are more likely to form a more continuous and cohesive whole (<xref ref-type="bibr" rid="B35">Wei et al., 2009</xref>)</td>
</tr>
<tr>
<td rowspan="5" align="center">Landscape Scale</td>
<td align="center">N.P.</td>
<td align="left">- The N.P. index measures the number of landscape patches (i.e., discrete areas of different landscape types) within a specific area. A higher N.P. index indicates higher landscape fragmentation (<xref ref-type="bibr" rid="B39">Xiaoping et al., 2009</xref>)</td>
</tr>
<tr>
<td align="center">LPI</td>
<td align="left">- Same as above</td>
</tr>
<tr>
<td align="center">CONTAG</td>
<td align="left">- The CONTAG index measures the degree of adjacency or connectivity between patches of different types within the landscape. A higher CONTAG index indicates better connectivity in the landscape (<xref ref-type="bibr" rid="B11">Han et al., 2018</xref>)</td>
</tr>
<tr>
<td align="center">FRAC_AM</td>
<td align="left">- Same as above</td>
</tr>
<tr>
<td align="center">SHDI</td>
<td align="left">- The SHDI index measures the biodiversity or landscape diversity within an area. A higher SHDI index indicates higher landscape diversity in that area (<xref ref-type="bibr" rid="B39">Xiaoping et al., 2009</xref>)</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2-3">
<title>2.2.3 Landscape ecological security index</title>
<p>
<list list-type="simple">
<list-item>
<p>(1) Selection of Landscape Ecological Security Index</p>
</list-item>
</list>
</p>
<p>The Landscape Ecological Security Index construction method based on landscape disturbance and vulnerability has been widely applied in many studies (<xref ref-type="bibr" rid="B7">Feng et al., 2016</xref>; <xref ref-type="bibr" rid="B28">Nitschke and Innes, 2008</xref>; <xref ref-type="bibr" rid="B30">Sheng et al., 2018</xref>; <xref ref-type="bibr" rid="B33">Sun and Chen, 2017</xref>; <xref ref-type="bibr" rid="B53">Zhang et al., 2019</xref>). Landscape disturbance represents the degree of disturbance to different landscape ecosystems (<xref ref-type="bibr" rid="B53">Zhang et al., 2019</xref>), calculated by weighting landscape fragmentation, isolation, and dominance (<xref ref-type="bibr" rid="B53">Zhang et al., 2019</xref>). Landscape fragmentation describes the continuity and integrity of natural or anthropogenic landscapes, typically referring to the degree of fragmentation or isolation between different landscape elements (such as forests, grasslands, water bodies, etc.) within an area (<xref ref-type="bibr" rid="B29">Qing et al., 2007</xref>). Landscape isolation refers to the spatial distance and degree of isolation between different landscape types within a region (<xref ref-type="bibr" rid="B14">Jin and Hu, 2003</xref>). Landscape dominance is usually determined by comparing the coverage area of different landscape types (<xref ref-type="bibr" rid="B2">Cheng et al., 2005</xref>; <xref ref-type="bibr" rid="B4">Daqing et al., 2005</xref>; <xref ref-type="bibr" rid="B13">Jiang et al., 2002</xref>; <xref ref-type="bibr" rid="B23">Lohbeck et al., 2016</xref>). Existing research considers landscape fragmentation as the most important factor in landscape disturbance, followed by isolation and dominance, thus assigning values of 0.5, 0.3, and 0.2 (<xref ref-type="bibr" rid="B9">Gao et al., 2021</xref>; <xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>), rrespectively. Landscape vulnerability represents the degree of change in landscape use types after disturbance (<xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>). Based on existing research (<xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>) and the actual situation of the study area, values were assigned as follows: urban and rural residential land is the most stable, assigned a value of 1; forests and grasslands are mostly located far from residential areas and relatively stable, assigned values of two and 3, respectively; farmland is susceptible to human disturbance, assigned a value of 4; water bodies, located close to residential and farmland areas, are greatly disturbed by human activities, assigned a value of 5; unused land and industrial and mining areas represent the most fragile ecosystems, assigned a value of 6. Using ArcGIS 10.3, the study area was divided into 1 km &#xd7; 1&#xa0;km grids, and each grid&#x2019;s Landscape Ecological Security Index was calculated. Finally, the Landscape Ecological Security Index for the study area was obtained through Kriging interpolation (<xref ref-type="table" rid="T2">Table 2</xref>).</p>
<p>
<list list-type="simple">
<list-item>
<p>(2) Determination of Landscape Ecological Security Standards.</p>
</list-item>
</list>
</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Calculation method of landscape ecological security index.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<italic>Name</italic>
</th>
<th align="center">
<italic>Calculation Method</italic>
</th>
<th align="center">
<italic>Explanation</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Landscape Fragmentation Degree</td>
<td align="center">
<inline-formula id="inf1">
<mml:math id="m1">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">N</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">
<italic>N</italic>
<sub>
<italic>i</italic>
</sub> and <italic>A</italic>
<sub>
<italic>i</italic>
</sub> represent the number and area of patches for landscape type <italic>i</italic>, respectively</td>
</tr>
<tr>
<td align="center">Landscape Separation Degree</td>
<td align="center">
<inline-formula id="inf2">
<mml:math id="m2">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">N</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mi mathvariant="normal">A</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">
<italic>A</italic> denotes the total landscape area (in hm<sup>2</sup>)</td>
</tr>
<tr>
<td align="center">Landscape Dominance Degree</td>
<td align="center">
<inline-formula id="inf3">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">D</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:mo>&#x2061;</mml:mo>
<mml:mi>ln</mml:mi>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>/</mml:mo>
<mml:mn>4</mml:mn>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mi>ln</mml:mi>
<mml:mo>&#x2061;</mml:mo>
<mml:msub>
<mml:mi mathvariant="normal">A</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">
<italic>P</italic>
<sub>
<italic>i</italic>
</sub> stands for the perimeter of landscape type I (in hm<sup>2</sup>)</td>
</tr>
<tr>
<td align="center">Landscape Disturbance Degree</td>
<td align="center">
<inline-formula id="inf4">
<mml:math id="m4">
<mml:mrow>
<mml:msub>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mi mathvariant="normal">a</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">C</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="normal">b</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mi mathvariant="normal">c</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">D</mml:mi>
<mml:mi mathvariant="normal">i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">
<italic>a</italic>, <italic>b</italic>, <italic>c</italic> are weights, assigned values of 0.5, 0.3, and 0.2, respectively (<xref ref-type="bibr" rid="B9">Gao et al., 2021</xref>; <xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>)</td>
</tr>
<tr>
<td align="center">Landscape Vulnerability</td>
<td align="center">Fi: The referencing method in literature</td>
<td align="left">Values are assigned based on the sensitivity of landscapes to disturbance, specifically: Unused land, mining land: 6; Water bodies: 5; Farmland: 4; Grassland: 3; Forest land: 2; Construction land: 1 (<xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>; <xref ref-type="bibr" rid="B40">Xie, 2011</xref>)</td>
</tr>
<tr>
<td align="center">Landscape Ecological Security Index</td>
<td align="center">
<inline-formula id="inf5">
<mml:math id="m5">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:msub>
<mml:mi>S</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>m</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mi>k</mml:mi>
<mml:mi>i</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:mn>10</mml:mn>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:math>
</inline-formula>
</td>
<td align="left">
<italic>m</italic> is the number of landscape types, and <italic>A</italic>
<sub>
<italic>ki</italic>
</sub> represents the area of landscape type <italic>i</italic> in the <italic>k</italic>th sampling unit</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>To reflect the spatiotemporal differences in the landscape ecological security index of the study area, reference to existing research results and combined with the actual distribution of the landscape ecological security index in the study area (<xref ref-type="bibr" rid="B9">Gao et al., 2021</xref>; <xref ref-type="bibr" rid="B34">Wang et al., 2021</xref>), the natural breakpoint method is adopted to divide the landscape ecological security level of the study area into five grades. The classification intervals for other years are standardized based on the classification intervals for the reference year (<xref ref-type="table" rid="T3">Table 3</xref>). To make the landscape ecological security pattern map more intuitive, each level of landscape ecological security is further divided into three equal parts while maintaining the same level intervals.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Landscape ecological security evaluation standards.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">
<italic>Landscape Ecological Security Index</italic>
</th>
<th align="center">
<italic>Landscape Ecological Security Level</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">0.151&#x2264; <inline-formula id="inf6">
<mml:math id="m6">
<mml:mrow>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.194</td>
<td align="center">Low Security</td>
</tr>
<tr>
<td align="center">0.194&#x2264; <inline-formula id="inf7">
<mml:math id="m7">
<mml:mrow>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.233</td>
<td align="center">Lower Security</td>
</tr>
<tr>
<td align="center">0.233&#x2264; <inline-formula id="inf8">
<mml:math id="m8">
<mml:mrow>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.280</td>
<td align="center">General Security</td>
</tr>
<tr>
<td align="center">0.280.&#x2264; <inline-formula id="inf9">
<mml:math id="m9">
<mml:mrow>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.341</td>
<td align="center">Safer</td>
</tr>
<tr>
<td align="center">0.341.&#x2264; <inline-formula id="inf10">
<mml:math id="m10">
<mml:mrow>
<mml:mi mathvariant="normal">E</mml:mi>
<mml:msub>
<mml:mi mathvariant="normal">S</mml:mi>
<mml:mi mathvariant="normal">k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.450</td>
<td align="center">Safe</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-2-4">
<title>2.2.4 Spatial autocorrelation analysis</title>
<p>Global spatial autocorrelation can reveal whether geographic phenomena exhibit spatial clustering and is often characterized using Moran&#x2019;s I (<xref ref-type="bibr" rid="B12">HuaLin, 2008</xref>). Local spatial autocorrelation reflects the degree of correlation between a geographic phenomenon or attribute in a local area unit within a larger region and the same phenomenon or attribute in adjacent local area units, and is often represented using LISA<sub>i</sub> indices (<xref ref-type="bibr" rid="B12">HuaLin, 2008</xref>). The calculation formulas are as follows (<xref ref-type="disp-formula" rid="e1">Equations 1</xref>, <xref ref-type="disp-formula" rid="e2">2</xref>).<list list-type="simple">
<list-item>
<p>(1) Global Spatial Autocorrelation:</p>
</list-item>
</list>
<disp-formula id="e1">
<mml:math id="m11">
<mml:mrow>
<mml:mtext>Moran</mml:mtext>
<mml:mo mathvariant="sans-serif">&#x2019;</mml:mo>
<mml:mi mathvariant="sans-serif">s</mml:mi>
<mml:mi>I</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>n</mml:mi>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>j</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:msub>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>1</mml:mn>
</mml:mrow>
<mml:mi>n</mml:mi>
</mml:munderover>
</mml:mstyle>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>
</p>
<p>In the equation, <italic>n</italic> represents the total number of sample points for variable <italic>x</italic>; <italic>x</italic>
<sub>
<italic>i</italic>
</sub> and <italic>x</italic>
<sub>
<italic>j</italic>
</sub> denote the observed values of a certain attribute in areas i and j, respectively; <inline-formula id="inf11">
<mml:math id="m12">
<mml:mrow>
<mml:mover accent="true">
<mml:mi mathvariant="normal">x</mml:mi>
<mml:mo>&#x20d0;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:math>
</inline-formula> denotes the average value of the attribute across all <italic>n</italic> areas; <italic>W</italic>
<sub>
<italic>ij</italic>
</sub> represents the spatial weight matrix. When the Moran&#x2019;s I index is positive, it indicates spatial clustering, meaning similar values are adjacent to each other; when the Moran&#x2019;s I index is negative, it indicates spatial dispersion, meaning similar values are far apart from each other; when the Moran&#x2019;s I index is close to 0, it indicates a random spatial distribution of patterns.<list list-type="simple">
<list-item>
<p>(2) Local Spatial Autocorrelation:</p>
</list-item>
</list>
<disp-formula id="e2">
<mml:math id="m13">
<mml:mrow>
<mml:msub>
<mml:mtext>LISA</mml:mtext>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mi>i</mml:mi>
</mml:munder>
</mml:mstyle>
<mml:mfrac>
<mml:mrow>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:mrow>
<mml:mrow>
<mml:mi>n</mml:mi>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:mfrac>
<mml:mrow>
<mml:mstyle displaystyle="true">
<mml:munder>
<mml:mo>&#x2211;</mml:mo>
<mml:mi>j</mml:mi>
</mml:munder>
</mml:mstyle>
<mml:mrow>
<mml:msub>
<mml:mi>W</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:msub>
<mml:mi>x</mml:mi>
<mml:mi>j</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mover accent="true">
<mml:mi>x</mml:mi>
<mml:mo>&#xaf;</mml:mo>
</mml:mover>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#x2001;</mml:mo>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="&#x7c;">
<mml:mrow>
<mml:mi>i</mml:mi>
<mml:mo>&#x2260;</mml:mo>
<mml:mi>j</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>
</p>
<p>The LISA<sub>i</sub> does not have a fixed range of values. A LISA<sub>i</sub> value greater than 0 indicates spatial clustering of similar values around the area unit, manifested as high-high or low-low clustering. A LISA<sub>i</sub> value less than 0 indicates spatial clustering of dissimilar values, manifested as high-low or low-high clustering.</p>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Land use change analysis</title>
<sec id="s3-1-1">
<title>3.1.1 Analysis of land use structural changes</title>
<p>Between 1980 and 2020, forest land, cultivated land, grassland, and unused land were the primary land-use types in Heihe City. Among them, forest land accounted for the largest area, though it steadily decreased over time, while the area of cultivated land and industrial mining land continuously expanded. The total area of industrial/mining land increased 32-fold, from 180.27 hm<sup>2</sup> in 1980 to 5,775.10 hm<sup>2</sup> in 2020. Spatially, forest land was widely distributed across various parts of Heihe City, with southern forest areas gradually being converted into farmland. Industrial and mining land was mainly concentrated in the southern and southwestern regions. Water bodies were primarily located in the southern part of the study area, while cultivated land was mainly distributed in the southwest and northeast. Unused land was more scattered, primarily dispersed in the southwestern part of the study area (<xref ref-type="table" rid="T4">Table 4</xref>; <xref ref-type="fig" rid="F2">Figure 2</xref>).</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Land use change in the study area from 1980 to 2020.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">
<break/>
<italic>Land Use Type</italic>
</th>
<th colspan="2" align="center">
<italic>1980</italic>
</th>
<th colspan="2" align="center">
<italic>1990</italic>
</th>
<th colspan="2" align="center">
<italic>2000</italic>
</th>
<th colspan="2" align="center">
<italic>2010</italic>
</th>
<th colspan="2" align="center">
<italic>2020</italic>
</th>
</tr>
<tr>
<th align="center">
<italic>Area/hm</italic>
<sup>
<italic>2</italic>
</sup>
</th>
<th align="center">
<italic>Percentage/%</italic>
</th>
<th align="center">
<italic>Area/hm</italic>
<sup>
<italic>2</italic>
</sup>
</th>
<th align="center">
<italic>Percentage/%</italic>
</th>
<th align="center">
<italic>Area/hm</italic>
<sup>
<italic>2</italic>
</sup>
</th>
<th align="center">
<italic>Percentage/%</italic>
</th>
<th align="center">
<italic>Area/hm</italic>
<sup>
<italic>2</italic>
</sup>
</th>
<th align="center">
<italic>Percentage/%</italic>
</th>
<th align="center">
<italic>Area/hm</italic>
<sup>
<italic>2</italic>
</sup>
</th>
<th align="center">
<italic>Percentage/%</italic>
</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Farmland</td>
<td align="center">1,042,380.35</td>
<td align="center">15.60</td>
<td align="center">1,191,385.92</td>
<td align="center">17.83</td>
<td align="center">1,611,517.76</td>
<td align="center">24.12</td>
<td align="center">1,766,664.60</td>
<td align="center">26.44</td>
<td align="center">1,626,910.82</td>
<td align="center">24.35</td>
</tr>
<tr>
<td align="center">Woodland</td>
<td align="center">4,050,119.12</td>
<td align="center">60.61</td>
<td align="center">4,014,045.89</td>
<td align="center">60.08</td>
<td align="center">3,592,790.19</td>
<td align="center">53.77</td>
<td align="center">3,458,221.23</td>
<td align="center">51.76</td>
<td align="center">3,569,805.22</td>
<td align="center">54.43</td>
</tr>
<tr>
<td align="center">Grassland</td>
<td align="center">872,989.17</td>
<td align="center">13.06</td>
<td align="center">664,838.34</td>
<td align="center">9.95</td>
<td align="center">673,624.75</td>
<td align="center">10.08</td>
<td align="center">719,621.61</td>
<td align="center">10.77</td>
<td align="center">657,851.47</td>
<td align="center">9.85</td>
</tr>
<tr>
<td align="center">Waters</td>
<td align="center">51,463.73</td>
<td align="center">0.70</td>
<td align="center">49,700.14</td>
<td align="center">0.74</td>
<td align="center">60,558.34</td>
<td align="center">0.91</td>
<td align="center">113,113.30</td>
<td align="center">1.70</td>
<td align="center">58,861.87</td>
<td align="center">0.88</td>
</tr>
<tr>
<td align="center">Urban and Rural Residential Land</td>
<td align="center">36,957.93</td>
<td align="center">0.55</td>
<td align="center">52,375.02</td>
<td align="center">0.78</td>
<td align="center">55,119.49</td>
<td align="center">0.82</td>
<td align="center">58,459.35</td>
<td align="center">0.87</td>
<td align="center">57,606.71</td>
<td align="center">0.86</td>
</tr>
<tr>
<td align="center">Industrial and Mining Land</td>
<td align="center">180.27</td>
<td align="center">0.09</td>
<td align="center">431.95</td>
<td align="center">0.01</td>
<td align="center">832.10</td>
<td align="center">0.01</td>
<td align="center">994.19</td>
<td align="center">0.01</td>
<td align="center">5,775.10</td>
<td align="center">0.09</td>
</tr>
<tr>
<td align="center">Unused land</td>
<td align="center">627,631.63</td>
<td align="center">9.39</td>
<td align="center">708,926.33</td>
<td align="center">10.61</td>
<td align="center">687,265.34</td>
<td align="center">10.29</td>
<td align="center">564,668.40</td>
<td align="center">8.45</td>
<td align="center">704,556.50</td>
<td align="center">10.54</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Land use changes in heihe city from 1980 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g002.tif"/>
</fig>
</sec>
<sec id="s3-1-2">
<title>3.1.2 Analysis of land use transfer changes</title>
<p>Between 1980 and 2020, the mutual conversion among cultivated land, forest land, grassland, and unused land was the main direction of land use type transitions in the study area. During these 40 years, 3,924.13 hm<sup>2</sup> of forest land were converted to cultivated land, primarily in the southern and northeastern parts of the study area. Additionally, 2,064.5 hm<sup>2</sup> of grassland were converted to cultivated land, distributed throughout various parts of Heihe City. This indicates that, alongside mining development, Heihe City has increasingly emphasized the sustainability of agriculture and the ecological environment. Faced with the driving force of urban construction and development, 154.09 hm<sup>2</sup> of cultivated land were converted to urban and rural residential construction land, along with 84.45 hm<sup>2</sup> of forest land and 55.04 hm<sup>2</sup> of grassland, with the main conversion areas located in the southwestern part of Heihe City. It was observed that the degree of land use transfer for industrial and mining land was relatively low, with only a small portion converted to cultivated land. The main reason for this is Heihe City&#x2019;s strong reliance on mineral resources, and the long-term environmental impacts caused by mining activities, which make it difficult for industrial and mining land to be repurposed for other uses in the short term after mining operations have ceased (Supplementary Table S1; <xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>The land use transfer patterns vary across different land categories in the study area over different years.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g003.tif"/>
</fig>
<p>From the perspective of different land use periods, the main land use transitions occurred between 2000&#x2013;2010 and 2010&#x2013;2020. During 2000&#x2013;2010, driven by the national policies of returning farmland to forest and grassland, along with the development of Heihe City&#x2019;s forestry industry, the primary land transitions involved forest land. A total of 2,560.98 hm<sup>2</sup> of forest land was converted to cultivated land, 2,534.68 hm<sup>2</sup> to grassland, and 1,192.94 hm<sup>2</sup> to unused land. During the same period, 1,156.12 hm<sup>2</sup> of cultivated land were converted to forest land, while grassland and unused land contributed 2,444.57 hm<sup>2</sup> and 1,465.97 hm<sup>2</sup>, respectively, to forest land (Supplementary Table S4; <xref ref-type="fig" rid="F3">Figure 3</xref>). In the period between 2010 and 2020, the most significant land transitions involved water bodies, with 1,124.71 hm<sup>2</sup> of water area being converted to other land types, while only 573.73 hm<sup>2</sup> of other land were converted to water bodies. During this period, the water area was reduced by half, primarily transitioning into cultivated land and unused land. This reduction suggests potential over-extraction and processing of water resources, leading to their depletion and conversion into cultivated land. Additionally, the failure of ecological restoration attempts in areas affected by mining development may have resulted in the transformation of water bodies into unused land (Supplementary Table S5; <xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
<p>Overall, from 1980 to 2020, the changes in cultivated and forest land were the most significant, with a scattered spatial distribution. The general trend shows forest land undergoing the most notable transformations. The land use pattern exhibited significant changes in the northeastern and southwestern regions, while the northwestern and southeastern areas remained relatively stable (Supplementary Table S1; <xref ref-type="fig" rid="F3">Figure 3</xref>).</p>
</sec>
</sec>
<sec id="s3-2">
<title>3.2 Landscape pattern change analysis</title>
<sec id="s3-2-1">
<title>3.2.1 Patch scale change analysis</title>
<p>Regarding patch contagion (COHESION), the patch indices for various land-use types were all close to 100, indicating a high degree of connectivity among the landscape patches, forming a relatively continuous and cohesive whole. However, the patch cohesion for urban and rural construction land and industrial and mining land was relatively low, though it gradually increased from 1980 to 2020. The primary reason for this increase is that Heihe City has vigorously developed its mining industry from 1980 to 2020, attracting a large number of mining-related industries and services, which led to the concentration of economic activities and, consequently, the focused distribution of industrial and mining land and urban-rural construction land. Additionally, the Heihe City government&#x2019;s &#x201c;General Land Use Plan for Heihe City (2006&#x2013;2020)&#x201d;(<ext-link ext-link-type="uri" xlink:href="https://www.heihe.gov.cn/hhs/c102620/202012/c11_225272.shtml">https://www.heihe.gov.cn/hhs/c102620/202012/c11_225272.shtml</ext-link>) (accessed on 8 October 2024) mandated adjustments in the scale of construction land, further enhancing the cohesion of urban-rural construction land and industrial and mining land (<xref ref-type="fig" rid="F4">Figure 4</xref>).</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Changes in patch scale landscape pattern indices from 1980 to 2020. <bold>(A)</bold> Changes in cohesion landscape pattern index from 1980 to 2020. <bold>(B)</bold> Changes in FRAC_AM landscape pattern index from 1980 to 2020. <bold>(C)</bold> Changes in the LPI landscape pattern index from 1980 to 2020 <bold>(D)</bold> changes in AREA_MN landscape pattern index from 1980 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g004.tif"/>
</fig>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Changes in Landscape Pattern Index at the landscape Scale from 1980 to 2020. <bold>(A)</bold> 1980&#x2013;2020 FRAC_AM Landscape Pattern Index Changes. <bold>(B)</bold> 1980&#x2013;2020 LPI Landscape Pattern Index Changes. <bold>(C)</bold> 1980&#x2013;2020 N.P. Landscape Pattern Index Changes. <bold>(D)</bold> 1980&#x2013;2020 SHDI Landscape Pattern Index Changes. <bold>(E)</bold> 1980&#x2013;2020 CONTAG Landscape Pattern Index Changes.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g005.tif"/>
</fig>
<p>In terms of area-weighted average patch fractal dimension (FRAC_AM), forest land and unused land exhibited relatively high values, close to 1.3. In contrast, the area-weighted average patch fractal dimensions for other land-use types were more balanced, hovering around 1.2 and 1.1. This indicates that the shapes of patches for forest land and unused land are relatively complex with jagged edges, while the patches for other land types have simpler shapes with smoother edges. By 2020, the area-weighted average patch fractal dimension for industrial and mining land had increased, mainly influenced by Heihe City&#x2019;s ecological restoration policies. Human interference has caused changes in the shapes of these industrial and mining areas. The &#x201c;General Planning for Mineral Resources in Heilongjiang Province (2016&#x2013;2020)&#x201d; (<ext-link ext-link-type="uri" xlink:href="https://www.hlj.gov.cn/hlj/c108411/202110/c00_31181663.shtml">https://www.hlj.gov.cn/hlj/c108411/202110/c00_31181663.shtml</ext-link>) (accessed on 8 October 2024) proposed land reclamation in mining areas, which effectively addresses ecological issues in these regions, resulting in more complex shapes for industrial/mining land. Furthermore, the rugged terrain of Heihe City, characterized by rolling mountains and intersecting valleys, has contributed to the increasing complexity of the boundaries of industrial and mining land as mineral extraction deepens, which is another reason for the rise in the fractal dimension of these patches (<xref ref-type="fig" rid="F4">Figure 4</xref>).</p>
<p>In terms of the largest patch index (LPI), forest land has consistently been the dominant landscape type in Heihe City. However, from 1990 to 2020, the LPI area of forest land gradually decreased, while the LPI area of cultivated land gradually increased. This shift can be attributed to Heihe City being primarily a forestry-based economy during this period, where economic development was closely tied to forest resources. As urban development progressed, activities such as deforestation and mining caused the originally continuous forest land to be fragmented into smaller patches of residential land, industrial and mining land, and unused land, resulting in spatial fragmentation (<xref ref-type="fig" rid="F4">Figure 4</xref>).</p>
<p>Regarding average patch area (AREA_MN), forest land had the largest AREA_MN, indicating the lowest degree of landscape fragmentation. This is because forest land accounted for 60% of the original land-use types in Heihe City, and resource extraction in mining areas was concentrated in regions planned by the Heihe City government. As a result, forest areas that were either undeveloped or minimally impacted by development were often retained on the periphery of mining zones or in difficult-to-develop areas, leading to relatively concentrated patches of forest land. The AREA_MN of forest land peaked in 1990 and subsequently began to decline, mirroring the trend of increasing fragmentation in unused land and cultivated land during the same period. This indicates that after 1990, the process of urban integration contributed to the increased fragmentation of forest land, cultivated land, and unused land.</p>
</sec>
<sec id="s3-2-2">
<title>3.2.2 Analysis of landscape-level changes</title>
<p>In terms of the area-weighted average patch fractal dimension (FRAC_AM), the study area showed a gradual decline from 1980 to 2020, decreasing by 0.02 over this period, indicating that the landscape in the study area has become simpler and the edges have become smoother. Regarding the largest patch index (LPI), the area of dominant landscapes in the study area gradually decreased from 1990 to 2010, but began to increase again from 2010 to 2020. This suggests a decline in the status of forest land as a dominant landscape, although it has somewhat recovered due to support from the return-to-forest policy. Concerning the number of patches (NP), there was a decline in the number of patches from 1980 to 1990, reaching a minimum value, indicating an increase in landscape connectivity in Heihe City in 1990, which enhanced the overall resilience of the landscape to disturbances. After 1990, the number of patches began to rise gradually, suggesting an increasing degree of landscape fragmentation in Heihe City. In terms of landscape diversity (SHDI), there has been an upward trend in landscape diversity in the study area since 1990. This increase is primarily associated with the rise in grassland area during the ecological restoration process, as well as the expansion of urbanization leading to an increase in the area of built-up and unused land. Regarding landscape contagion (CONTAG), the contagion index in the study area began to decline from 1990, dropping by 4% by 2010. This indicates that during this period, the degree of adjacency or connectivity between different patch types in the landscape of Heihe City weakened, while the extent of human activities affecting the natural landscape increased, coinciding with a gradual rise in urbanization levels.</p>
</sec>
</sec>
<sec id="s3-3">
<title>3.3 Landscape ecological security pattern evolution</title>
<sec id="s3-3-1">
<title>3.3.1 Overall analysis of landscape ecological security pattern</title>
<p>From a general perspective, between 1980 and 2020, the spatial distribution of landscape ecological security levels in the study area showed a trend of being lower in the south and higher in the north. Over time, the overall landscape ecological security pattern has shown improvement. Specifically, in 1980, the proportions of low-security and relatively low-security areas in the study area were 28.37% and 36.03% respectively, mainly concentrated in the southern part with a fragmented landscape of agricultural land, forests, urban and rural construction land, and unused land. Areas classified as relatively safe and safe were primarily located in the north and east, covering 12.08% of the area, predominantly composed of forests with more intact land use. Overall, the ecological security level of Heihe City was relatively low in 1980. By 2020, the proportions of low-security and relatively low-security areas had decreased to 22.14% and 26.48% respectively, showing increases of 6.23% and 9.55%. This improvement was particularly notable in areas such as Aihui District, Wudalianchi City, and Bei&#x2019;an City, likely due to integrated planning of agricultural land and controlled zoning for mineral development (<xref ref-type="fig" rid="F6">Figure 6</xref>).</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Evolution of ecological security landscape pattern in the study area from 1980 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g006.tif"/>
</fig>
<p>From 1980 to 1990, the ecological security levels in the study area continued to exhibit a pattern of lower levels in the south and higher levels in the north, but the ecological security in the low-security and relatively low-security areas in the southern and central regions improved. The proportion of low-security and relatively low-security areas decreased from 64.40% in 1980 to 41.11% in 1990, indicating an improvement in the ecological environment of Heihe City during this period, with an overall increase in landscape ecological security. From 1990 to 2000, the ecological security levels in the study area gradually shifted from the pattern of lower in the south and higher in the north to one of lower in the central region and higher in the surrounding areas. The environmental quality of Sunwu County improved, while that of Wudalianchi City and Nenjiang County deteriorated. From 2000 to 2010, the ecological security levels in the study area saw an overall improvement, forming a landscape ecological pattern of lower in the central region and higher in the surrounding areas, with particularly improved environmental quality in Wudalianchi City, where most areas were classified as relatively safe or safe. However, from 2010 to 2020, with the accelerated urbanization process, the ecological quality of Heihe City declined, and the overall ecological security level decreased. This decline was especially prominent in the central and southern regions, where areas classified as safe transitioned to relatively low-security areas, primarily due to an increase in unused land and a reduction in landscape restoration capacity (<xref ref-type="fig" rid="F6">Figure 6</xref>).</p>
</sec>
<sec id="s3-3-2">
<title>3.3.2 Spatial autocorrelation analysis of landscape ecological security index</title>
<p>From 1980 to 2020, the global spatial autocorrelation indices for the five periods were 0.892, 0.992, 1.029, 1.038, and 1.103, all of which passed the significance test at <italic>p</italic> &#x3d; 0.01. This indicates that the ecological security index exhibits spatial clustering. Spatially, two extremes in Heihe City&#x2019;s ecological security index were observed, with relatively few areas showing high-low or low-high clustering. Most areas exhibited no significant correlation, with local spatial autocorrelation primarily characterized by high-high value clustering and low-low value clustering.</p>
<p>In 1980, the clusters of ecological security indices were relatively scattered, with low-low value clusters appearing around the periphery. By 1990 and 2000, the spatial distribution of ecological security in Heihe City became more stable. Locally, high-high value clusters were mainly concentrated in the northern and eastern regions, indicating a relatively good level of ecological security in these areas. Meanwhile, the low-low value clusters around the edges gradually expanded southward, reflecting a decline in ecological security in the southern region, with a growing risk of ecological degradation. The spatial autocorrelation results for the ecological security index suggest that the northern and eastern regions, where forests are distributed, acted as ecological barriers, effectively maintaining the environment. In contrast, the central and southern regions, due to the development of industrial, mining, and urban construction land, experienced more frequent human disturbances, making the ecosystems more fragile. During the 2010 and 2020 periods, the area of high-high value clusters in the northern region became more stable, while the low-low value clusters in the central region shifted entirely to the south. This was due to the central region&#x2019;s unused land patches gradually forming larger, more regular patches under proper urban planning, reducing ecological damage in the surrounding areas. However, urban and rural construction in the southern region began to impact the environment, highlighting the need for future attention to ecosystem restoration and protection in this area (<xref ref-type="fig" rid="F7">Figure 7</xref>).</p>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Local spatial autocorrelation map of ecological security index in the study area from 1980 to 2020 (<italic>p</italic> &#x3d; 0.01)</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g007.tif"/>
</fig>
</sec>
</sec>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<sec id="s4-1">
<title>4.1 Nsights on landscape pattern index</title>
<sec id="s4-1-1">
<title>4.1.1 Screening and analysis of landscape pattern index</title>
<p>Evaluating the spatial pattern changes of mining cities using landscape pattern indices is meaningful and necessary because it can provide a reference for policymakers in planning mining cities (<xref ref-type="bibr" rid="B35">Wei et al., 2009</xref>). In selecting landscape pattern indices, covering every aspect in practical applications is optional. Landscape pattern indices serve repetitive purposes, and only a few indices are needed to fully meet the requirements of pattern analysis, reflecting the characteristics of landscape patterns in the study area. Therefore, this study screened landscape factors, and nine indices representing landscape fragmentation, aggregation, and diversity were selected for calculation.</p>
<p>Compared to the study by Ya et al., which only analyzed the spatial changes of the landscape from a temporal perspective (<xref ref-type="bibr" rid="B42">Ya-Hui et al., 2021</xref>), this study further divided landscape indices into patch-level and landscape-level analyses. It described the changes in landscape patterns of different land types in the study area over different time dimensions and the overall landscape pattern changes in the study area over different time dimensions, providing a comprehensive analysis from macro and micro perspectives, making the results more scientifically sound. From the results of landscape pattern analysis, it can be observed that human disturbances significantly impact the spatial pattern of a mining city.</p>
</sec>
<sec id="s4-1-2">
<title>4.1.2 Analysis of landscape pattern index results</title>
<p>Based on the analysis of landscape pattern indices, the landscape pattern of Heihe City has been primarily influenced by economic dependence, policy restrictions, and environmental issues.</p>
<p>From the patch scale perspective, significant changes were observed in the average patch area (AREA_MN) and the largest patch index (LPI) for forest land between 1980 and 2020. However, other landscape patterns showed little variation. This is mainly because the forestry industry holds an important position in the economic structure of Heihe City. The forestry industry is one of Heihe&#x2019;s traditional sectors, with abundant forest resources, particularly high-quality timber such as red pine and white birch. The city has actively developed industries related to timber processing, non-timber forest product development, and forest tourism. With the growing awareness of sustainable development and ecological protection, Heihe City has gradually promoted the sustainable use of forest resources and the conservation of its forest resources in recent years. Over the past 40 years, mining activities and land use policies in Heihe City have been subject to strict government regulation and oversight. The local government has controlled the area of forest land, effectively increasing its extent by implementing afforestation projects on reclaimed farmland and unused land in various forest farms. Consequently, the changes in the landscape pattern of forest patches have been the most significant.</p>
<p>From the landscape scale perspective, the landscape diversity index of Heihe City showed an upward trend from 1980 to 2020, while other indices displayed a downward trend. This pattern is closely related to the mining policies, ecological restoration policies, and urbanization processes implemented in Heihe City. The city is rich in mineral resources, with a diverse range of types distributed across various patches, which facilitates economic diversification. Under this model of mining development, Heihe City has experienced significant changes in its economic structure, which has also influenced the updating and alteration of the landscape pattern. Additionally, due to the transformation of the industrial structure, the local government has adopted policies for the redevelopment and reuse of abandoned mining areas or extraction sites. This involves converting mined-out land into green spaces, parks, wetlands, or other forms of urban land. For instance, the Xigangzi coal mining area in Heihe City has combined mine environmental restoration with farmland occupation and compensation projects to transform industrial wasteland into fertile land. This has enriched the types of landscape patches within the city, thereby increasing landscape diversity. At the same time, Heihe City has initiated ecological restoration projects aimed at restoring ecosystems in mining areas through vegetation recovery, soil remediation, and water body management. These efforts have diversified the landscape structure of damaged areas and increased the complexity of the ecosystem. Furthermore, with the acceleration of the urbanization process, the government has expanded land types for various uses, such as commercial zones, residential areas, industrial districts, and transportation facilities. This urbanization process has clarified the boundaries between different land uses, increasing the number and variety of landscape patches and enhancing overall landscape diversity.</p>
</sec>
</sec>
<sec id="s4-2">
<title>4.2 Ecological security evaluation insights</title>
<sec id="s4-2-1">
<title>4.2.1 Optimization insights for ecological security evaluation</title>
<p>Due to the adoption of different standards, the results of landscape ecological security assessment of mining cities obtained from different research objects can only present conclusions partially consistent with those of previous studies (<xref ref-type="bibr" rid="B1">Chen and Ma, 2023</xref>). Scholars typically combined ecological security assessment with evaluation model studies in previous research. For instance, Yin et al. combined landscape ecological security assessment with the PSR model (<xref ref-type="bibr" rid="B47">Yin et al., 2012</xref>); While Zhang et al. integrated ecological security assessment with the Pressure-State-Response framework (<xref ref-type="bibr" rid="B52">Zhang R. et al., 2023</xref>). However, due to differences in evaluation objects and evaluation models, the scientific validity of the assessment results cannot be guaranteed. Additionally, evaluation criteria have become one of the controversial issues in landscape ecological security assessment results. For example, Wang et al. constructed an ecological security evaluation system through the evaluation of ecosystem service importance and ecological sensitivity, mainly using the Analytic Hierarchy Process (AHP) combined with expert questionnaires to obtain evaluation criteria, leading to highly subjective conclusions (<xref ref-type="bibr" rid="B57">Zhou et al., 2022</xref>). In this study, however, theoretical models were not relied upon, and all research was based on the land use changes in Heihe City from 1980 to 2020. In constructing ecological evaluation criteria, quantitative research methods were adopted, and the landscape ecological security assessment results obtained by calculating landscape pattern indices reduced the subjectivity and errors introduced by expert scoring methods.</p>
<p>Moreover, this study did not use a single ecological security evaluation system. Still, it supplemented the landscape ecological security assessment results through spatial autocorrelation analysis of landscape ecological security spaces, helping to deepen the understanding of the spatial characteristics and dynamic changes of ecosystems, thereby better assessing the health status of ecosystems and guiding ecological protection and management efforts. Spatial autocorrelation analysis of landscape ecological security spaces was used to examine the overall analysis of landscape ecological security indices, validating the feasibility of the ecological security assessment results.</p>
</sec>
<sec id="s4-2-2">
<title>4.2.2 Ecological regulation measures in ecological security areas</title>
<p>Considering the economic development characteristics and natural environment conditions of Heihe City, and based on the 2020 ecological security evaluation results, corresponding management measures are proposed according to the different ecological security levels and characteristics, in order to reduce ecological security risks to a controllable range. The division of ecological functional areas from the perspective of sustainable development in mining cities has not yet been proposed for Heihe City (<xref ref-type="bibr" rid="B31">Shi and Xiang, 2023</xref>). This study, based on Heihe City&#x2019;s ecological security index and higher-level planning, uses methods such as overlaying current land use data, DEM data, and administrative boundaries, along with spatial interpolation analysis, to divide Heihe City into four functional areas: the ecological security control area for mining cities, ecological protection area, ecological optimization area, and ecological comprehensive development area (<xref ref-type="fig" rid="F8">Figure 8</xref>) (<xref ref-type="bibr" rid="B32">Song et al., 2023</xref>).</p>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Ecological security zoning of heihe city.</p>
</caption>
<graphic xlink:href="fenvs-12-1488439-g008.tif"/>
</fig>
<p>In the low and relatively low ecological security areas of Heihe City, there is a phenomenon of low-low value clustering, which is highly sensitive to human interference (<xref ref-type="bibr" rid="B32">Song et al., 2023</xref>). These clusters appear in agricultural and forested lands, where the ecological security status is poor, and serious ecological degradation and related issues are present. This region has been designated as the ecological security control area for mining cities. The total area of this control area is 4,308,956&#xa0;hm<sup>2</sup>, accounting for 64.52% of the total area. Key areas include the central and southern parts of Nenjiang County, Wudalianchi City, and Bei&#x2019;an City, as well as the northern part of Nenjiang County and the central region of Aihui District. In this area, ecological regulation efforts should focus on reversing the trend of ecological retrogression, with particular attention to post-development land quality protection and soil restoration. Conservation and rehabilitation measures should be implemented to protect agricultural and forested lands. By constructing ecological corridors and restoring ecological nodes and barrier points, all development activities that could cause pollution should be prohibited, minimizing the impact of human activities on the existing ecological environment (<xref ref-type="bibr" rid="B50">Zhang et al., 2017</xref>; <xref ref-type="bibr" rid="B49">Zhang et al., 2018</xref>; <xref ref-type="bibr" rid="B51">Zhang et al., 2024</xref>; <xref ref-type="bibr" rid="B17">Li Q. et al., 2023</xref>).</p>
<p>The general ecological security area of Heihe City is in a state of ecological alert, where both low-high and high-low value clustering phenomena occur, indicating that the area is sensitive to human interference. Although the ecological environment has been partially damaged, it still retains some capacity for self-recovery (<xref ref-type="bibr" rid="B8">Feng et al., 2023</xref>). This region has been classified as the ecological security protection area for mining cities, mainly distributed in the northern parts of Xunke County, Sunwu County, and the northern part of Aihui District, covering an area of 15,948&#xa0;hm<sup>2</sup>. The dominant landscape types in this area include degraded unused land (bare rock), large expanses of high-vegetation-cover paddy fields, and human-disturbed forested areas (shrubland) and farmland. Therefore, priority should be given to protecting the degraded unused land and forested areas in this area, preventing the expansion of land use driven by the pursuit of economic gains, and reducing unreasonable human development activities. For the degraded unused land in this area, measures such as land leveling, increasing green cover, and prohibiting further development should be taken to maintain ecological security levels (<xref ref-type="bibr" rid="B36">Wei et al., 2023</xref>). As for the human-disturbed forested areas, strict implementation of the &#x201c;Grain for Green&#x201d; project should be enforced to prevent the over-expansion of farmland for economic development purposes.</p>
<p>The relatively safe areas of Heihe City are regions with strong resilience to disturbances, exhibiting high-high or high-low value clustering. These areas have an intact ecological structure and a high level of ecological security. This region has been designated as the ecological optimization area for mining cities, primarily located in Aihui District, covering a total area of 6,419&#xa0;hm<sup>2</sup>. The dominant landscape type in this area is undisturbed forestland, which can serve as an ecological barrier and support sustainable regional developments (<xref ref-type="bibr" rid="B36">Wei et al., 2023</xref>). Therefore, it is essential to continue promoting active forestland protection policies, implementing the &#x201c;Grain for Green&#x201d; and &#x201c;Wetland Restoration&#x201d; policies, and intensifying efforts to protect the forested landscape in this area (<xref ref-type="bibr" rid="B8">Feng et al., 2023</xref>). Human activity should be minimized to reduce negative impacts, and unreasonable agricultural practices such as deforestation for farming must be strictly limited. These measures will help maintain the natural ecological succession of Heihe City while balancing social, ecological, and economic benefits.</p>
<p>The safe areas in Heihe City are areas with strong self-recovery capabilities, where high-high value clustering is evident. These areas have the highest level of ecological security and can accommodate moderate development. This region has been designated as the ecological comprehensive development area, located in the northern part of Aihui District, covering a total area of 1,326&#xa0;hm<sup>2</sup>. Although this area enjoys a high level of ecological security, it also faces significant environmental pressures. In the course of development, efforts should leverage the concentration of forestland to enhance infrastructure development, promote the clustering of the forestry industry, and prioritize the advancement of modern forestry. Additionally, it is important to vigorously develop the service sector, accommodate industrial transfers from central cities, and manage population transfers from restricted and prohibited development areas. This approach aims to prevent ecological degradation caused by blind development.</p>
</sec>
</sec>
<sec id="s4-3">
<title>4.3 Limitations and future directions of the study</title>
<p>This study provides important insights into Chinese mining cities&#x2019; ecological protection and planning development. The assessment results are a scientific basis for government agencies to formulate ecological protection policies and plans. For mining cities at different stages of development, this study can help determine development strategies and propose appropriate ecological protection measures to ensure the sustainable development of landscape ecological security in mining cities. Despite rigorous data processing and indicator selection, the choice of landscape pattern index, the accuracy of remote sensing data, and the determination of ecological evaluation criteria may all affect the assessment results. To further improve the accuracy of the assessment, future research could consider increasing the analysis of ecological risk assessment to validate the accuracy of the results. While this study analyzed the spatiotemporal dynamics of ecological security levels in the study area and explored their spatial correlations, it did not investigate the driving factors affecting ecological security levels. In the future, it is necessary to conduct a quantitative analysis of the driving factors behind ecological security changes, investigate land use conflicts, and perform multi-scenario simulations and predictions. This will help develop systematic strategies that focus on improving the ecological environment baseline, formulating reasonable spatial layouts, and planning for tourism development to enhance the ecological security status of the research area and promote the sustainable and coordinated development of geological parks and regional ecosystem safety. Furthermore, the ecological regulation measures and planning development strategies proposed in this paper are framed from a policy planning perspective; specific measures need to be implemented through the establishment of an ecological security framework.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s5">
<title>5 Conclusion</title>
<p>In this study, for the first time in its kind, the ecological security level of mining cities is integrated with the development strategies of mining cities.Based on the land use data from Heihe, Heilongjiang, from 1980 to 2020, this study analyzed the spatial pattern changes in Heihe using landscape pattern indices. By constructing a landscape ecological assessment system and conducting spatial autocorrelation analysis, the ecological security level of Heihe from 1980 to 2020 was determined. The results of this study are also applicable to mining cities worldwide that are in similar stages. Compared to other studies, this research adopts a &#x201c;dual evaluation index,&#x201d; making the constructed ecological security assessment system more comprehensive and the assessment results more objective. The main research conclusions are as follows.<list list-type="simple">
<list-item>
<p>(1)&#x2009;&#x2009;This article reveals the spatiotemporal evolution of land use in Heihe City from 1980 to 2020. Through land use data and land transfer matrices, the study found that over the 40&#xa0;years, the spatial pattern of land use exhibited overall stability with significant local changes. By 2020, the current land use in Heihe City primarily consists of arable land, forest land, and unused land. Arable land and forest land have gradually shifted towards unused land. This trend indicates that during the urbanization process, mineral exploitation can lead to land degradation.</p>
</list-item>
<list-item>
<p>(2)&#x2009;&#x2009;Using Fragstats 4.2 software, four patch-level indices (AREA_MN, COHESION, FRAC_AM, LPI) and five landscape-level indices (FRAC_AM, LPI, NP, SHDI, CONTAG) were calculated to analyze the spatiotemporal dynamics of the landscape pattern in Heihe City from 1980 to 2020. The study found that over these 40 years, at the patch scale, except for forest land, other patches experienced severe fragmentation, had simpler shapes, and smoother edges. The patches of various land use types showed a high degree of cohesion, with forest land and arable land being the dominant landscapes in the study area. At the landscape scale, the overall landscape diversity in Heihe City increased, while landscape fragmentation and contagion decreased, indicating an optimization of the landscape pattern. Therefore, in future planning, attention should continue to be paid to the connectivity of Heihe City&#x2019;s landscape spatial pattern.</p>
</list-item>
<list-item>
<p>(3)&#x2009;&#x2009;Spatial autocorrelation analysis was conducted to assess and analyze the ecological security level of Heihe City and its spatiotemporal evolution characteristics from 1980 to 2020. During this period, the overall ecological security index of Heihe City showed an upward trend, but the ecological security index in the southern region remained low. Additionally, the landscape ecological security index displayed spatial clustering, mainly characterized by high-high value clusters and low-low value clusters. High-high value clusters were primarily located in the northern forested areas of the study region, while low-low value clusters were mainly found in the southern mixed land use areas. Therefore, the government of Heihe City should focus on the ecological security of the southern region and take proactive measures to improve it.</p>
</list-item>
<list-item>
<p>(4)&#x2009;&#x2009;Based on the ecological security index, the urban ecological security space of Heihe City was regulated by dividing the city into mining urban ecological security control areas, ecological protection areas, ecological optimization areas, and comprehensive ecological development areas, corresponding to low and lower security areas, general ecological security areas, relatively secure areas, and secure areas, respectively. The ecological security control area focuses on reversing the trend of ecological retrogression. The ecological protection area emphasizes protecting the degraded unused land and forest land within the area. The ecological optimization area promotes active forest protection and ecological policies, while the comprehensive ecological development area aims to improve infrastructure construction and promote industrial diversification. The research results contribute to maximizing ecological and economic benefits in mining cities, promoting the sustainable development of such cities, enriching the theoretical framework, and guiding similar mining cities.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s6">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec sec-type="author-contributions" id="s7">
<title>Author contributions</title>
<p>YS: Conceptualization, Data curation, Investigation, Methodology, Software, Validation, Writing&#x2013;original draft, Writing&#x2013;review and editing. XF: Data curation, Methodology, Supervision, Writing&#x2013;original draft, Writing&#x2013;review and editing. HC: Conceptualization, Data curation, Investigation, Methodology, Software, Writing&#x2013;original draft. YK: Conceptualization, Formal Analysis, Methodology, Writing&#x2013;original draft. MS: Conceptualization, Data curation, Investigation, Methodology, Writing&#x2013;original draft.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This research was funded by &#x201c;The Art and Science Planning Foundation of Tianjin, grant number A16038&#x201d;, &#x201c;The Social Science Foundation of Hebei, grant number HB17YS015&#x201d;, &#x201c;The outcomes of the first-class curriculum&#x2019;s Innovative Design A&#x2019;s in Tianjin Municipality.&#x201d; and &#x201c;The Industry -Education Cooperation and Education project of the Ministry of Education grant number 221003711084739&#x201d;.</p>
</sec>
<ack>
<p>We sincerely thank the two reviewers for their careful review of this article, which has greatly improved its rigor. We are also very grateful to Mr. Zhang Qingtai from Guangxi Arts University for providing technical guidance for this paper. Additionally, on a personal note, I would like to deeply thank Mr. Xia Xusheng for his companionship and support, and for the encouragement he has given me throughout my three years of graduate study. I sincerely hope that he can continue to accompany me throughout my academic career.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Chen</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Ecological risk identification and ecological security pattern construction of productive wetland landscape</article-title>. <source>Water Resour. Manag.</source> <volume>37</volume>, <fpage>4709</fpage>&#x2013;<lpage>4731</lpage>. <pub-id pub-id-type="doi">10.1007/s11269-023-03574-1</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Cheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Study on the change of agriculture landscape pattern based on GIS--A case study from Maping Town</article-title>. <source>Chin. J. eco-agriculture</source> <volume>13</volume>, <fpage>184</fpage>&#x2013;<lpage>186</lpage>.</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Corry</surname>
<given-names>R. C.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Characterizing fine-scale patterns of alternative agricultural landscapes with landscape pattern indices</article-title>. <source>Landsc. Ecol.</source> <volume>20</volume>, <fpage>591</fpage>&#x2013;<lpage>608</lpage>. <pub-id pub-id-type="doi">10.1007/s10980-004-5036-8</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Daqing</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Junbin</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yushan</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Study on heterogeneity and diversity of landscape ecosystem for slopeland vegetation</article-title>. <source>Res. Soil Water Conservation</source> <volume>12</volume>, <fpage>163</fpage>&#x2013;<lpage>165</lpage>.</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Dong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Coupling coordination analysis of industrial mining land, landscape pattern and carbon storage in a mining city: a case study of Ordos, China</article-title>. <source>Geomatics Nat. Hazards and Risk</source> <volume>14</volume>. <pub-id pub-id-type="doi">10.1080/19475705.2023.2275539</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2005</year>). <article-title>Problems, reasons and strategies for sustainable development of mining cities in China</article-title>. <source>J. Nat. Resour.</source> <volume>20</volume>, <fpage>68</fpage>&#x2013;<lpage>77</lpage>.</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Analysis of ecological vulnerability in yuyang district based on landscape pattern</article-title>. <source>Res. Soil Water Conservation</source> <volume>23</volume>, <fpage>179</fpage>&#x2013;<lpage>184</lpage>. <pub-id pub-id-type="doi">10.13869/j.cnki.rswc.2016.06.023</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feng</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Strategies and practical paths for ecological restoration and comprehensive management in yulin coal mining area in the new era</article-title>. <source>Northwest. Geol.</source> <volume>56</volume>, <fpage>19</fpage>&#x2013;<lpage>29</lpage>. <pub-id pub-id-type="doi">10.12401/j.nwg.2023087</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>B.-P.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.-M.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>K.-J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Landscape ecological risk assessment and influencing factors in ecological conservation area in Sichuan-Yunnan provinces, China</article-title>. <source>Ying yong sheng tai xue bao &#x3d; J. Appl. Ecol.</source> <volume>32</volume>, <fpage>1603</fpage>&#x2013;<lpage>1613</lpage>. <pub-id pub-id-type="doi">10.13287/j.1001-9332.202105.018</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Gao</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lai</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Sensing the multi-scale landscape functions heterogeneity by big geodata from parcel to urban agglomerations-a case of the Greater Bay Area, China</article-title>. <source>Int. J. Digital Earth</source> <volume>17</volume>. <pub-id pub-id-type="doi">10.1080/17538947.2024.2304073</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Han</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Rao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Correlation analysis between cultivated land quality and landscape pattern index in south hilly areas</article-title>. <source>Jiangsu J. Agric. Sci.</source> <volume>34</volume>, <fpage>1057</fpage>&#x2013;<lpage>1065</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1000-4440.2018.05.013</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>HuaLin</surname>
<given-names>X. I. E.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Regional eco-risk analysis of based on landscape structure and spatial statistics</article-title>. <source>Acta Ecol. Sin.</source> <volume>28</volume>, <fpage>5020</fpage>&#x2013;<lpage>5026</lpage>.</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Analysis on the spatial pattern dynamic of forest landscape in Anji Mountainous District, Zhejiang province</article-title>. <source>J. Fujian Coll. For.</source> <volume>22</volume>, <fpage>150</fpage>&#x2013;<lpage>153</lpage>. <pub-id pub-id-type="doi">10.1324/j.cnki.jfcf.2002.02.015</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jin</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>B.</given-names>
</name>
</person-group> (<year>2003</year>). <article-title>Discussion on several indices assessing landscape dispersion</article-title>. <source>Ying yong sheng tai xue bao &#x3d; J. Appl. Ecol.</source> <volume>14</volume>, <fpage>314</fpage>&#x2013;<lpage>316</lpage>. <pub-id pub-id-type="doi">10.13287/j.1001-9332.2003.0071</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jingwei</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Haize</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2006</year>). <article-title>Application of landscape index in landscape pattern based on damaike wetland nature reserve in anshan</article-title>. <source>Res. Soil Water Conservation</source> <volume>13</volume>, <fpage>230</fpage>&#x2013;<lpage>233</lpage>. <pub-id pub-id-type="doi">10.13203/j.whugis20140269</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Qu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>Landscape pattern analysis for mining area based on optimal grain size</article-title>. <source>Geomatics Inf. Sci. Wuhan Univ.</source> <volume>41</volume>, <fpage>939</fpage>&#x2013;<lpage>945</lpage>.</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Ai</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Z.</given-names>
</name>
<etal/>
</person-group> (<year>2023a</year>). <article-title>Theoretical analysis and practical exploration on ecological restoration of mines with multi-source solid wastes: example from yulin city, shaanxi province</article-title>. <source>Northwest. Geol.</source> <volume>56</volume>, <fpage>70</fpage>&#x2013;<lpage>77</lpage>. <pub-id pub-id-type="doi">10.12401/j.nwg.2023093</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Pang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2002</year>). <article-title>Study on countermeasures to realize sustainable development for mining aress</article-title>. <source>China Soft Sci.</source> <volume>0</volume>, <fpage>16</fpage>&#x2013;<lpage>19</lpage>.</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Ai</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Ren</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>M.</given-names>
</name>
<etal/>
</person-group> (<year>2023b</year>). <article-title>Ecological risk assessment and restoration area identification of Pengyang County on the basis of the landscape pattern and function</article-title>. <source>Environ. Monit. Assess.</source> <volume>195</volume>, <fpage>998</fpage>. <pub-id pub-id-type="doi">10.1007/s10661-023-11596-w</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2023c</year>). <article-title>Study on the construction of ecological security pattern in kunming based on MSPA and MCR model</article-title>. <source>J. Ecol. Rural Environ.</source> <volume>39</volume>, <fpage>69</fpage>&#x2013;<lpage>79</lpage>. <pub-id pub-id-type="doi">10.19741/j.issn.1673-4831.2022.0211</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023a</year>). <article-title>Development status of coal mining in China</article-title>. <source>J. South. Afr. Inst. Min. Metall.</source>, <volume>123</volume>, <fpage>19</fpage>&#x2013;<lpage>27</lpage>. <pub-id pub-id-type="doi">10.17159/2411-9717/1506/2023</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2023b</year>). <article-title>The spatio-temporal patterns and driving forces of land use in the context of urbanization in China: evidence from nanchang city</article-title>. <source>Int. J. Environ. Res. public health</source> <volume>20</volume>, <fpage>2330</fpage>. <pub-id pub-id-type="doi">10.3390/ijerph20032330</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lohbeck</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Bongers</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Martinez-Ramos</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Poorter</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2016</year>). <article-title>The importance of biodiversity and dominance for multiple ecosystem functions in a human-modified tropical landscape</article-title>. <source>Ecology</source> <volume>97</volume>, <fpage>2772</fpage>&#x2013;<lpage>2779</lpage>. <pub-id pub-id-type="doi">10.1002/ecy.1499</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>T.-T.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>N.-N.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Applications of 2D and 3D landscape pattern indices in landscape pattern analysis of mountainous area at county level</article-title>. <source>Ying yong sheng tai xue bao &#x3d; J. Appl. Ecol.</source> <volume>23</volume>, <fpage>1351</fpage>&#x2013;<lpage>1358</lpage>. <pub-id pub-id-type="doi">10.13287/j.1001-9332.2012.0182</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Luo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Guan</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>P.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Evolution of green space pattern and its ecosystem services change in mining cities: a case study of xuzhou city</article-title>. <source>Resour. Environ. Yangtze Basin</source> <volume>32</volume>, <fpage>1686</fpage>&#x2013;<lpage>1697</lpage>. <pub-id pub-id-type="doi">10.11870/cjlyzyyhj202308012</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lv</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ji</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Chao</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Integrated assessment and prediction of ecological security in typical ecologically fragile areas</article-title>. <source>Environ. Monit. Assess.</source> <volume>196</volume>, <fpage>286</fpage>. <pub-id pub-id-type="doi">10.1007/s10661-024-12453-0</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Mardonova</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>Y. S.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Environmental, hydrological, and social impacts of coal and nonmetal minerals mining operations</article-title>. <source>J. Environ. Manag.</source> <volume>332</volume>, <fpage>117387</fpage>. <pub-id pub-id-type="doi">10.1016/j.jenvman.2023.117387</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Nitschke</surname>
<given-names>C. R.</given-names>
</name>
<name>
<surname>Innes</surname>
<given-names>J. L.</given-names>
</name>
</person-group> (<year>2008</year>). <article-title>Integrating climate change into forest management in South-Central British Columbia: an assessment of landscape vulnerability and development of a climate-smart framework</article-title>. <source>For. Ecol. Manag.</source> <volume>256</volume>, <fpage>313</fpage>&#x2013;<lpage>327</lpage>. <pub-id pub-id-type="doi">10.1016/j.foreco.2008.04.026</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qing</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ranghui</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Shixin</surname>
<given-names>W. U.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Spatial- temporal dynamic changes of landscape fragmentation in the upper reaches of tarim river</article-title>. <source>J. Arid Land Resour. Environ.</source> <volume>21</volume>, <fpage>73</fpage>&#x2013;<lpage>77</lpage>. <pub-id pub-id-type="doi">10.1348/j.cnki.jalre.2007.09.001</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sheng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Analysis of landscape vulnerability in youyu county based on landscape pattern and ecological sensibility</article-title>. <source>J. Northwest For. Univ.</source> <volume>33</volume>, <fpage>167</fpage>&#x2013;<lpage>174</lpage>. <pub-id pub-id-type="doi">10.3969/j.issn.1001-7461.2018.04.28</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shi</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Xiang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Suitability evaluation of underground space development in coal mining cities in Henan province</article-title>. <source>J. Min. Sci. Technol.</source> <volume>8</volume>, <fpage>714</fpage>&#x2013;<lpage>724</lpage>. <pub-id pub-id-type="doi">10.19606/j.cnki.jmst.2023.05.013</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Song</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Gu</surname>
<given-names>H.-H.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>F.-P.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>S.-P.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.-X.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Environmental assessments in dense mining areas using remote sensing information over Qian?an and Qianxi regions China</article-title>. <source>Ecol. Indic.</source> <volume>146</volume>, <fpage>109814</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2022.109814</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Sun</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>F.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Assessment on ecological health of landscape of the CoastalWetlands in the mouth of the yalu river</article-title>. <source>Wetl. Sci.</source> <volume>15</volume>, <fpage>40</fpage>&#x2013;<lpage>46</lpage>. <pub-id pub-id-type="doi">10.13248/j.cnki.wetlandsci.2017.01.006</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Spatio-temporal differentiation of landscape pattern evolution and its ecological security in karst areas: a case study of trough valley in northeast Guizhou Province</article-title>. <source>Acta Ecol. Sin.</source> <volume>41</volume>, <fpage>7273</fpage>&#x2013;<lpage>7291</lpage>. <pub-id pub-id-type="doi">10.5846/stxb202004160913</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>Q. I.</given-names>
</name>
<name>
<surname>Yan-bo</surname>
<given-names>Q. U.</given-names>
</name>
<name>
<surname>Hong-yi</surname>
<given-names>L. I. U.</given-names>
</name>
<name>
<surname>Le</surname>
<given-names>L. I.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Selection of representative landscape pattern indexes and land use division</article-title>. <source>China Land Sci.</source> <volume>23</volume>, <fpage>33</fpage>&#x2013;<lpage>37</lpage>. <pub-id pub-id-type="doi">10.13708/j.cnki.cn11-2640.2009.01.006</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wei</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Kong</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Construction of ecological security pattern and quantitative diagnosis of ecological problems in Shenmu City</article-title>. <source>Acta Ecol. Sin.</source> <volume>43</volume>, <fpage>82</fpage>&#x2013;<lpage>94</lpage>. <pub-id pub-id-type="doi">10.5846/stxb202204120969</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>C.-X.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>S.-F.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Q.-H.</given-names>
</name>
</person-group> (<year>2014</year>). <article-title>Ecological risk assessment of land use based on exploratory spatial data analysis (ESDA): a case study of Haitan Island, Fujian Province</article-title>. <source>Ying yong sheng tai xue bao &#x3d; J. Appl. Ecol.</source> <volume>25</volume>, <fpage>2056</fpage>&#x2013;<lpage>2062</lpage>. <pub-id pub-id-type="doi">10.13287/j.1001-9332.20140425.004</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lei</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Bian</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Landscape ecological network construction controlling surface coal mining effect on landscape ecology: a case study of a mining city in semi-arid steppe</article-title>. <source>Ecol. Indic.</source> <volume>133</volume>, <fpage>108403</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2021.108403</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xiaoping</surname>
<given-names>L. I. U.</given-names>
</name>
<name>
<surname>Xia</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Yimin</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Q. I. N.</given-names>
</name>
<name>
<surname>Shaoying</surname>
<given-names>L. I.</given-names>
</name>
<name>
<surname>Minghui</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2009</year>). <article-title>Landscape expansion index and its applications to quantitative analysis of urban expansion</article-title>. <source>Acta Geogr. Sin.</source> <volume>64</volume>, <fpage>1430</fpage>&#x2013;<lpage>1438</lpage>.</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xie</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2011</year>). <article-title>Spatial characteristic analysis of land use eco-risk based on landscape structure:a case study in the Xingguo County, Jiangxi Province</article-title>. <source>China Environ. Sci.</source> <volume>31</volume>, <fpage>688</fpage>&#x2013;<lpage>695</lpage>.</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Qin</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Observational study on the change of vegetation carbon sink in cities with depleted coal resources in transition:Xuzhou as an example</article-title>. <source>J. China Coal Soc.</source> <volume>48</volume>, <fpage>2916</fpage>&#x2013;<lpage>2924</lpage>. <pub-id pub-id-type="doi">10.13225/j.cnki.jccs.cn22.1674</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ya-Hui</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>W.-B.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Landscape pattern of built-up land in a typical mining city, Datong, Shanxi, China from 1986 to 2018</article-title>. <source>Ying yong sheng tai xue bao &#x3d; J. Appl. Ecol.</source> <volume>32</volume>, <fpage>1614</fpage>&#x2013;<lpage>1622</lpage>. <pub-id pub-id-type="doi">10.13287/j.1001-9332.202105.011</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Future urban waterlogging simulation based on LULC forecast model: a case study in Haining City, China</article-title>. <source>Sustain. Cities Soc.</source> <volume>87</volume>, <fpage>104167</fpage>. <pub-id pub-id-type="doi">10.1016/j.scs.2022.104167</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Ye</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Identification of ecological networks in the Guangdong-Hong Kong-Macao Greater Bay Area based on habitat quality assessment</article-title>. <source>Acta Ecol. Sin.</source> <volume>43</volume>, <fpage>10430</fpage>&#x2013;<lpage>10442</lpage>. <pub-id pub-id-type="doi">10.20103/j.stxb.202211173329</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yao</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Cheng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Evaluation of ecological security pattern and optimization suggestions in Minjiang River Basin based on MCR model and gravity model</article-title>. <source>Acta Ecol. Sin.</source> <volume>43</volume>, <fpage>7083</fpage>&#x2013;<lpage>7096</lpage>. <pub-id pub-id-type="doi">10.20103/j.stxb.202209062546</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ye</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Shang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Spatiotemporal evolution law and driving force of mining city patterns</article-title>. <source>Environ. Sci. Pollut. Res.</source> <volume>29</volume>, <fpage>10291</fpage>&#x2013;<lpage>10307</lpage>. <pub-id pub-id-type="doi">10.1007/s11356-021-16488-5</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Qiu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Pan</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2012</year>). <article-title>Ecological security assessment of land use in small towns based on PSR model A case study of 22 towns of tongnan</article-title>. <source>J. Southwest China Normal Univ. Nat. Sci. Ed.</source> <volume>37</volume>, <fpage>126</fpage>&#x2013;<lpage>130</lpage>. <pub-id pub-id-type="doi">10.13718/j.cnki.xsxb.2012.02.015</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liesch</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Goldscheider</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2023a</year>). <article-title>Global analysis of land-use changes in karst areas and the implications for water resources</article-title>. <source>Hydrogeology J.</source> <volume>31</volume>, <fpage>1197</fpage>&#x2013;<lpage>1208</lpage>. <pub-id pub-id-type="doi">10.1007/s10040-023-02650-5</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Lin</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Research advance on ecological security in China: assessment models and management measures</article-title>. <source>J. Ecol. Rural Environ.</source> <volume>34</volume>, <fpage>1057</fpage>&#x2013;<lpage>1063</lpage>. <pub-id pub-id-type="doi">10.11934/j.issn.1673-4831.2018.12.001</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Jing</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2017</year>). <article-title>Assessment of the ecological security of water environment in Henan based on PSR Model</article-title>. <source>Ecol. Sci.</source> <volume>36</volume>, <fpage>49</fpage>&#x2013;<lpage>54</lpage>. <pub-id pub-id-type="doi">10.14108/j.cnki.1008-8873.2017.05.007</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Jia</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Research on the sustainability of greening process in the Mu Us Sandy Land based on the spatiotemporal stability of ecological land</article-title>. <source>Plos One 19</source> <volume>19</volume>, <fpage>e0292469</fpage>. <pub-id pub-id-type="doi">10.1371/journal.pone.0292469</pub-id>
</citation>
</ref>
<ref id="B52">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Xiong</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2023b</year>). <article-title>Ecological security assessment of China based on the Pressure-State-Response framework</article-title>. <source>Ecol. Indic.</source> <volume>154</volume>, <fpage>110647</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2023.110647</pub-id>
</citation>
</ref>
<ref id="B53">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Qu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>New method of constructing landscape ecological vulnerability index and rationality verification</article-title>. <source>Sci. Surv. Mapp.</source> <volume>44</volume>, <fpage>129</fpage>&#x2013;<lpage>136</lpage>. <pub-id pub-id-type="doi">10.16251/j.cnki.1009-2307.2019.11.019</pub-id>
</citation>
</ref>
<ref id="B54">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>W.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Integrated approach for ecological restoration and ecological spatial network optimization with multiple ecosystem functions in mining areas</article-title>. <source>Ecol. Indic.</source> <volume>156</volume>, <fpage>111141</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2023.111141</pub-id>
</citation>
</ref>
<ref id="B55">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Prediction of ecological security patterns based on urban expansion: a case study of Chengdu</article-title>. <source>Ecol. Indic.</source> <volume>158</volume>, <fpage>111467</fpage>. <pub-id pub-id-type="doi">10.1016/j.ecolind.2023.111467</pub-id>
</citation>
</ref>
<ref id="B56">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Tang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Economy development of mining cities of China:data and analysis</article-title>. <source>China Ind. Econ.</source>, <fpage>43</fpage>&#x2013;<lpage>51</lpage>. <pub-id pub-id-type="doi">10.19581/j.cnki.ciejournal.2004.03.006</pub-id>
</citation>
</ref>
<ref id="B57">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Z.-C.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>S.-Y.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>G.-D.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Construction of ecological security pattern based on ecological sensitivity assessment in jining city, China</article-title>. <source>Pol. J. Environ. Stud.</source> <volume>31</volume>, <fpage>5383</fpage>&#x2013;<lpage>5404</lpage>. <pub-id pub-id-type="doi">10.15244/pjoes/151866</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>