<?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">1654175</article-id>
<article-id pub-id-type="doi">10.3389/fenvs.2025.1654175</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>Interrelationships and zoning-based management of landscape ecological risk and ecological resilience in the Hefei metropolitan circle from a multi-scale perspective</article-title>
<alt-title alt-title-type="left-running-head">Liu 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.2025.1654175">10.3389/fenvs.2025.1654175</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name>
<surname>Liu</surname>
<given-names>Junhao</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3112853/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/resources/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
<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" corresp="yes">
<name>
<surname>Li</surname>
<given-names>Jiulin</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>
<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/Writing - review &#x26; editing/"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/supervision/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
</contrib>
<contrib contrib-type="author">
<name>
<surname>Wang</surname>
<given-names>Yingxing</given-names>
</name>
<xref ref-type="aff" rid="aff1">
<sup>1</sup>
</xref>
<xref ref-type="aff" rid="aff2">
<sup>2</sup>
</xref>
<uri xlink:href="https://loop.frontiersin.org/people/3125295/overview"/>
<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/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
</contrib>
</contrib-group>
<aff id="aff1">
<sup>1</sup>
<institution>School of Architecture and Planning, Anhui Jianzhu University</institution>, <addr-line>Hefei</addr-line>, <country>China</country>
</aff>
<aff id="aff2">
<sup>2</sup>
<institution>Anhui Collaborative Innovation Centre for Urbanization</institution>, <addr-line>Hefei</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/2208858/overview">Muniyandi Balasubramanian</ext-link>, Institute for Social and Economic Change, India</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/3070399/overview">Md Julfikar Ali</ext-link>, Aliah University, India</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/3150204/overview">Umanath Malaiarasan</ext-link>, Madras Institute of Development Studies, India</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Jiulin Li, <email>ljiul90@163.com</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>10</day>
<month>09</month>
<year>2025</year>
</pub-date>
<pub-date pub-type="collection">
<year>2025</year>
</pub-date>
<volume>13</volume>
<elocation-id>1654175</elocation-id>
<history>
<date date-type="received">
<day>27</day>
<month>06</month>
<year>2025</year>
</date>
<date date-type="accepted">
<day>26</day>
<month>08</month>
<year>2025</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2025 Liu, Li and Wang.</copyright-statement>
<copyright-year>2025</copyright-year>
<copyright-holder>Liu, Li and Wang</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>Assessing the ecosystem resilience and risks of metropolitan areas helps coordinate ecological environmental protection, reduce landscape ecological risks, and enhance ecological resilience. Existing studies often analyze the trade-offs and synergies between landscape ecological risks and ecological resilience within specific regions based on administrative divisions. However, the mechanisms underlying these trade-offs and synergies across different scales remain unclear. This study focuses on the Hefei Metropolitan Area in China, constructing quantitative assessment models for ecological resilience and landscape ecological risks from a landscape pattern perspective. It systematically examines the spatiotemporal evolution characteristics of these two aspects from 2010 to 2020 across multiple scales, including grid, county, and city levels. Pearson correlation analysis and bivariate spatial autocorrelation are employed to investigate their spatial interactions. Furthermore, a coupling coordination model is used to identify trade-offs and synergies between the two at different scales. Finally, based on the interaction analysis, tailored zoning regulation strategies are proposed. Key findings reveal: (1) Over the past decade, the landscape ecological risk in the Hefei Metropolitan Circle has exhibited characteristics of &#x201C;stable overall pattern with localized differentiation,&#x201d; where high-risk zones are predominantly concentrated in areas with elevated water body ratios. Ecological resilience exhibited a slight declining trend, with an overall low resilience level. High-resilience zones experienced a significant reduction of 50.6% in area, while low-resilience zones expanded toward the core urbanization areas. (2) A robust negative correlation emerged between ecological resilience and landscape risk, intensifying significantly with finer spatial scales. Crucially, the recoverability dimension exerted the most potent counteracting effect on risk propagation. (3) The mean coupling coordination degree between ecological resilience and ecological risk across multiple scales registered below 0.5, with spatial heterogeneity manifesting as &#x201C;higher in the southwest and generally low in other areas.&#x201d; Coordination zones are sporadically distributed within southwest woodlands, while dysregulation zones concentrate predominantly in water bodies and urbanization-intensive belts. (4) Collectively, low-resilience/high-risk zones clustered in waterbody-dense areas, contrasting with high-resilience/low-risk clusters predominantly occupying forested biomes. This spatial divergence exhibited scale-dependent amplification under finer resolutions. Synthesizing multiscale correlations and coupling coordination patterns, we delineate a four-type governance typology. This study provides both theoretical underpinnings and practical pathways for multi-scale collaborative governance to enhance regional ecological resilience and advance sustainable development.</p>
</abstract>
<kwd-group>
<kwd>Hefei metropolitan circle</kwd>
<kwd>multi-scale</kwd>
<kwd>resilience-risk</kwd>
<kwd>spatial correlation</kwd>
<kwd>couplingcoordination degree</kwd>
<kwd>zoning governance</kwd>
</kwd-group>
<counts>
<page-count count="20"/>
</counts>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Social-Ecological Urban Systems</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec id="s1">
<title>1 Introduction</title>
<p>Since its Reform and Opening Up, China has undergone the most extensive and rapid urbanization in human history (<xref ref-type="bibr" rid="B46">Zhang and Wang, 2025</xref>). During 4&#xa0;decades of economic growth and spatial proliferation, human built infrastructure has superseded vast natural land covers, driving transformative reconfiguration of regional ecosystem structures (<xref ref-type="bibr" rid="B10">Li et al., 2024</xref>; <xref ref-type="bibr" rid="B19">Meimei et al., 2023</xref>; <xref ref-type="bibr" rid="B24">Ren et al., 2024</xref>; <xref ref-type="bibr" rid="B32">Wang C. et al., 2025</xref>; <xref ref-type="bibr" rid="B40">Yang et al., 2018</xref>). These dramatic land use transformations have not only accelerated biodiversity erosion and habitat fragmentation (<xref ref-type="bibr" rid="B6">Fletcher et al., 2018</xref>; <xref ref-type="bibr" rid="B26">Su et al., 2024</xref>; <xref ref-type="bibr" rid="B42">Yang et al., 2025</xref>), but also unleashed cascading landscape ecological perturbations that critically compromise regional ecological security patterns (<xref ref-type="bibr" rid="B30">Wang and Yang, 2025</xref>; <xref ref-type="bibr" rid="B45">Zeng et al., 2024</xref>; <xref ref-type="bibr" rid="B51">Zhu et al., 2024</xref>). Consequently, there is a critical imperative to advance holistic understanding of the process structure function nexus within land use change dynamics. Landscape ecological risk (LER) constitutes the probability of deleterious ecological consequences arising from ecosystems subjected to natural or anthropogenic stressors (<xref ref-type="bibr" rid="B8">Kong et al., 2025</xref>; <xref ref-type="bibr" rid="B21">Qiao et al., 2025</xref>; <xref ref-type="bibr" rid="B22">Qiu et al., 2025</xref>; <xref ref-type="bibr" rid="B31">Wang et al., 2023</xref>; <xref ref-type="bibr" rid="B44">Yu and Li, 2024</xref>; <xref ref-type="bibr" rid="B47">Zhang H. et al., 2023</xref>; <xref ref-type="bibr" rid="B49">Zhang N. et al., 2023</xref>) Ecological resilience (ER) denotes the capacity of ecosystems to maintain core structures, functions, and biodiversity through dynamic processes of resistance, adaptation, and recovery when confronting internal/external stressors. This capacity fundamentally hinges on the structural stability of regional ecosystem configurations (<xref ref-type="bibr" rid="B5">Feng et al., 2024</xref>; <xref ref-type="bibr" rid="B7">Jiang et al., 2025</xref>; <xref ref-type="bibr" rid="B11">Li K. et al., 2025</xref>; <xref ref-type="bibr" rid="B27">Tang et al., 2025</xref>; <xref ref-type="bibr" rid="B29">Wang and Liu, 2024</xref>; <xref ref-type="bibr" rid="B43">Yin et al., 2024</xref>). As a geospatial entity integrating biophysical and anthropogenic elements, landscapes experience amplified ecological risks when subjected to urban expansion induced terrain modifications and intensified anthropogenic pressures. These perturbations compromise structural stability of ecosystems, ultimately diminishing ER through self-reinforcing degradation cycles. Hence, deciphering the coupled dynamics between LER and resilience constitutes a pivotal theoretical frontier one essential for resolving ecological security dilemmas in rapidly urbanizing regions and advancing nature positive futures.</p>
<p>Currently, LER and ER stand as pivotal concerns in ecological security research, attracting extensive scholarly engagement. Most scholars regard ER and LER as two independent scientific issues, conducting extensive research on their assessment methods, spatiotemporal patterns, and influencing factors, respectively. At the theoretical level, LER research focuses on quantifying potential threats to ecosystems under anthropogenic activities or natural disturbances. Predominant assessment frameworks employ the &#x201c;disturbance-vulnerability-loss&#x201d; model, integrating Landscape Pattern Index (LPI) analysis, probabilistic modeling, and multi-source data fusion techniques. These approaches enable dynamic risk monitoring through risk classification, spatial heterogeneity analysis, and identification of driving factors (<xref ref-type="bibr" rid="B33">Wang F. et al., 2025</xref>; <xref ref-type="bibr" rid="B41">Yang et al., 2024</xref>; <xref ref-type="bibr" rid="B50">Zhang et al., 2024</xref>). ER theory emphasizes an ecosystem&#x2019;s adaptive and restorative capacities in response to disturbances. The research paradigm has progressively shifted from static assessments toward dynamic simulations, increasingly employing the &#x201c;resistance-adaptation-recovery&#x201d; framework. This approach reveals the evolutionary dynamics of system resilience through resilience network modeling, threshold identification, and adaptive cycle models (<xref ref-type="bibr" rid="B34">Wang S. et al., 2025</xref>; <xref ref-type="bibr" rid="B37">Wu et al., 2025</xref>). For instance, Shamuxi et al. developed a three dimensional assessment framework based on landscape disturbance, vulnerability, and loss dimensions to evaluate the spatiotemporal evolution of ecological risk in the Tarim River Basin, providing a scientific basis for regional ecological security (<xref ref-type="bibr" rid="B25">Shamuxi et al., 2025</xref>). Li et al. quantified ecological risk within the riparian buffer zone of the Luo River in Shaanxi Province using the LER Index (LERI) method, establishing a theoretical framework for regional risk management (<xref ref-type="bibr" rid="B14">Li Z. et al., 2025</xref>). Fan et al., employing a Pressure-State-Response (PSR) ER assessment framework, provided new insights into the dynamic interplay between regional spatial transformation and ER (<xref ref-type="bibr" rid="B2">Fan and Wei, 2025</xref>). Currently, research integrating ecological resilience and ecological risk remains limited both at the theoretical level and in conceptual exploration, with a lack of relevant empirical case studies. Some scholars have linked ecosystem services or ecological resilience with ecological risk to establish relatively balanced research frameworks. For instance, Wang et al. assessed the dynamics of ecosystem services and ecological risks by combining historical analysis and future predictions, identifying four distinct functional ecological zones (<xref ref-type="bibr" rid="B35">Wang X. et al., 2025</xref>). Xu et al. proposed an adaptive zoning framework through a comprehensive analysis of ecosystem service value and landscape ecological risk (<xref ref-type="bibr" rid="B38">Xu F et al., 2025</xref>). However, although existing studies have yielded fruitful results, merely using ecosystem services or ecological resilience to represent an ecosystem&#x2019;s capacity to mitigate risks is insufficient. There remains significant room for deepening theoretical integration and practical application: most research remains confined to a single perspective and lacks exploration into systematically integrating the &#x201c;disturbance-vulnerability-loss&#x201d; risk model with the &#x201c;resistance-adaptation-recovery&#x201d; resilience framework. Such integration is crucial for comprehensively capturing the complex feedback loops between the transmission pathways of risk pressures and the dynamic responses of ecosystem coping capacities during urbanization. Furthermore, most studies focus on descriptive analyses of spatiotemporal evolution patterns and differentiation characteristics, while actionable zoning management and collaborative optimization strategies based on risk-resilience coupling mechanisms are relatively scarce. This gap makes it difficult to effectively support the refined management of regional sustainable development.</p>
<p>As highly heterogeneous and dynamic complex ecosystems, metropolitan areas cannot be fully captured in terms of the interactive relationships between LER and ER through a single-scale analysis. Relying solely on adjacent scales may overlook their scale dependent characteristics, thereby making it difficult to reveal the mechanisms by which the same factor influences the tradeoffs and synergies between the two across multiple scales. This could ultimately lead to mismatched management strategies at different scale levels. Against the backdrop of rapid urbanization in China exacerbating regional land use changes, investigating the tradeoffs and synergies between LER and ER (and their subsystems) across multiple scales is crucial for future urban management. However, current academic understanding of the cross-scale tradeoffs and synergistic variations between the two remains insufficient, which constrains regional managers&#x2019; ability to regulate key influencing factors and limits the enhancement of regional ER. Therefore, from a landscape pattern perspective, exploring the tradeoffs and synergies between LER and ER (including their subsystems) across multiple scales, and identifying key resilience subsystems that constrain risks, holds significant theoretical and practical importance. Such an approach can guide managers at different levels in rapidly urbanizing areas to develop targeted, context-specific strategies for reducing regional LER and enhancing ER.</p>
<p>The Hefei Metropolitan Area serves as a key growth engine driving the rise of central China, a strategic core zone safeguarding the ecological shield of Chaohu Lake, and a crucial hub promoting coordinated regional development by connecting the eastern and western parts of the country. Its ecological security construction and context specific ecological zoning management are of great significance. Since the 21st century, as one of the fastest growing emerging metropolitan areas, Hefei has undergone rapid urbanization, with extensive natural surfaces being replaced by artificial built environments, placing severe pressure on the regional ecosystem. However, existing research lacks exploration of the spatiotemporal characteristics, tradeoffs, and synergies between LER and ER across multiple scales in the Hefei Metropolitan Area. This gap poses challenges for policymakers in developing evidence based, tailored management strategies. Therefore, this study selects the Hefei Metropolitan Area as the research region to bridge the aforementioned theory-practice gap. Its core theoretical contribution lies in the innovative integration of the &#x201c;disturbance-vulnerability-loss&#x201d; LER model with the &#x201c;resistance-adaptation-recovery&#x201d; ER framework to establish a coupled assessment system. This integration is particularly suited for analyzing metropolitan systems characterized by high intensity disturbances, significant vulnerability, and strong spatial heterogeneity. It enables the simultaneous examination of the spatiotemporal evolution and interactive coupling relationships between ER and risks across different scales. Through coupled assessment, key zones such as &#x201c;high-risk-low-resilience,&#x201d; &#x201c;low-risk-high-resilience,&#x201d; and &#x201c;dynamic risk-resilience areas&#x201d; can be identified, providing a direct basis for formulating differentiated zoning management strategies and enhancing the ecological security and sustainable development resilience of the Hefei Metropolitan Area. This study addresses the following questions: 1) How do the tradeoffs and synergies between LER and ER (including its subsystems) differ at the 3 &#xd7; 3&#xa0;km grid, county, and city scales? 2) How do the mechanisms underlying these tradeoffs and synergies at different scales influence zoning management in the metropolitan area? The detailed research workflow is illustrated in <xref ref-type="fig" rid="F1">Figure 1</xref>.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Research framework diagram.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g001.tif">
<alt-text content-type="machine-generated">Flowchart depicting a study on ecological resilience (ER) and landscape ecological risk (LER) across three sections. Section 1 shows original datasets and research areas using maps and data sources. Section 2 analyzes spatiotemporal dynamics in ER and LER, including modeling and assessment methods. Section 3 examines interactions between ER and LER, with zoning management and recommendations. Maps and charts illustrate each step.</alt-text>
</graphic>
</fig>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Study area</title>
<p>The Hefei Metropolitan Circle is located at the geographical pivot of eastern China (30&#xb0;57&#x2032;&#x2013;32&#xb0;32&#x2032;N, 116&#xb0;41&#x2032;&#x2013;118&#xb0;52&#x2032;E) (<xref ref-type="fig" rid="F2">Figure 2</xref>) (<xref ref-type="bibr" rid="B1">Dai et al., 2025</xref>; <xref ref-type="bibr" rid="B3">Fan et al., 2022</xref>). The region exhibits a composite geomorphic structure: the Jiang Huai Divide transects its eastern western axis, while the northern sector predominantly comprises the extended Huang Huai Plain. Centrally, undulating hills and interlaced valleys dominate, whereas the southern zone features an interdigitated distribution of riparian wetlands and low elevation foothills. Hydrologically, the area demonstrates characteristic dual basin dynamics, with the Yangtze and Huai Rivers constituting primary drainage arteries, supplemented by 18 major tributaries including the Chao Lake and Chu River systems. This geographic template engenders three distinct ecological gradients: an eastern riparian wetland corridor along the Yangtze, a central complex urban ecotone, and a western ecological barrier formed by the Dabie Mountain foothills. Collectively, they foster diverse natural ecosystems and anthropogenic environments, constituting a region of pivotal ecological significance. However, accelerated urbanization in recent years has imposed significant pressure on regional ecosystem stability through expanding urban footprints and intensified anthropogenic disturbances. This has triggered a marked increase in LER frequency across the region. As the primary economic growth hub in central China, the Hefei Metropolitan Circle now necessitates urgent investigation into the spatiotemporal dynamics of its ecological &#x201c;resilience-risk&#x201d; nexus. Targeted zoning governance strategies must be developed to advance regional high quality development and ecological civilization advancement.</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>Overview of the study area: <bold>(a)</bold> Location of the study area in China; <bold>(b)</bold> Elevation map; <bold>(c)</bold> Land use type map.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g002.tif">
<alt-text content-type="machine-generated">Map of China highlighting Anhui Province and the study area in red. Inset (b) shows elevation with colors from red (high) to blue (low). Inset (c) depicts land use types like forest, grassland, and cultivated land.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s2-2">
<title>2.2 Data sources</title>
<p>The year 2010 marked the conclusion of China&#x2019;s 11th Five Year Plan, 2015 was both the final year of the 12th Five Year Plan and a pivotal year for environmental legislation, while 2020 represented the end of the 13th Five Year Plan and the year of China&#x2019;s carbon neutrality commitment. Therefore, selecting these 3&#xa0;years for studying ecological resilience and ecological risks allows for capturing the entire process of China&#x2019;s macro strategic transition from &#x201c;growth-first&#x201d; to &#x201c;green development.&#x201d; Accordingly, this study utilizes land use data and socioeconomic data from the three periods of 2010, 2015, and 2020. All remote sensing imagery underwent standardized preprocessing including image cropping, atmospheric correction, and geometric rectification prior to analysis. Detailed data sources are provided in <xref ref-type="table" rid="T1">Table 1</xref>.</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Sources of information.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Data name</th>
<th align="center">Spatial resolution</th>
<th align="center">Data source</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Land use data</td>
<td align="center">30&#xa0;m</td>
<td align="center">China Land Cover Dataset (CLCD) at 30&#xa0;m Resolution (<ext-link ext-link-type="uri" xlink:href="https://zenodo.org/record/5210928">https://zenodo.org/record/5210928&#x23;.Y2t2uXZByMp</ext-link>)</td>
</tr>
<tr>
<td align="center">Socioeconomic data</td>
<td align="center">&#x2014;&#x2014;</td>
<td align="center">Municipal Statistical Yearbooks</td>
</tr>
<tr>
<td align="center">Crop prices and yields</td>
<td align="center">&#x2014;&#x2014;</td>
<td align="center">&#x300a;National Compilation of Cost-Benefit Data for Agricultural Products&#x300b;</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s2-3">
<title>2.3 Research methods</title>
<sec id="s2-3-1">
<title>2.3.1 ER construction</title>
<p>Resistance (P): Ecosystems with rich biodiversity and stable structures inherently possess enhanced functional integrity, enabling them to mitigate external pressures through self-regulatory mechanisms and thereby exhibit greater resistance. Consequently, this study quantifies resistance using Ecosystem Service Value (ESV). This approach integrates multiple functional dimensions including provisioning, regulating, and supporting services within a unified assessment framework. Through a value coefficient matrix that captures the ecological contribution of distinct land cover types, we objectively characterize regional ecosystem resistance levels (<xref ref-type="bibr" rid="B28">Wang and Ge, 2024</xref>). See specifically <xref ref-type="disp-formula" rid="e1">Equation 1</xref>.<disp-formula id="e1">
<mml:math id="m1">
<mml:mrow>
<mml:mi mathvariant="normal">P</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mtext>ESV</mml:mtext>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2211;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(1)</label>
</disp-formula>Where: <inline-formula id="inf1">
<mml:math id="m2">
<mml:mtext>ESV</mml:mtext>
</mml:math>
</inline-formula> &#x3d; Ecosystem Service Value; <inline-formula id="inf2">
<mml:math id="m3">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Area of land use type k; <inline-formula id="inf3">
<mml:math id="m4">
<mml:mrow>
<mml:mi>V</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mrow>
<mml:mi>f</mml:mi>
<mml:mi>k</mml:mi>
</mml:mrow>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Value coefficient of the <italic>f</italic>th ecosystem service for land use type <italic>k</italic>.</p>
<p>Adaptive Capacity (<inline-formula id="inf4">
<mml:math id="m5">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>): When ecosystems face natural disasters or anthropogenic disturbances, the dynamic equilibrium of their landscape structure is disrupted, requiring progression through reorganizational phases to reestablish stability. This study evaluates adaptive capacity by focusing on stabilizing mechanisms of landscape structure, specifically selecting heterogeneity and connectivity as primary determinants due to their irreplaceable ecological functions. In alignment with prevailing conventions in landscape ecology (<xref ref-type="bibr" rid="B37">Wu et al., 2025</xref>) (<xref ref-type="table" rid="T2">Table 2</xref>), we assign equal weighting to heterogeneity and connectivity in quantifying adaptive capacity.</p>
<table-wrap id="T2" position="float">
<label>TABLE 2</label>
<caption>
<p>Landscape indices and their weights.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Subsystem</th>
<th align="center">Landscape lndex</th>
<th align="center">Weight</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="2" align="center">Landscape heterogeneity</td>
<td align="center">Shannon&#x2019;s Diversity Index</td>
<td align="center">0.25</td>
</tr>
<tr>
<td align="center">Area-Weighted Mean Patch Fractal Dimension</td>
<td align="center">0.25</td>
</tr>
<tr>
<td align="center">Landscape connectivity</td>
<td align="center">Landscape Fragmentation Index</td>
<td align="center">0.5</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>Recovery Capacity (<inline-formula id="inf5">
<mml:math id="m6">
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>): This study evaluates ecosystem recovery capacity based on the ER model proposed by Peng et al. (<xref ref-type="bibr" rid="B20">Peng et al., 2015</xref>). See specifically <xref ref-type="disp-formula" rid="e2">Equation 2</xref>.<disp-formula id="e2">
<mml:math id="m7">
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mo>&#x2211;</mml:mo>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>R</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(2)</label>
</disp-formula>Where: <inline-formula id="inf6">
<mml:math id="m8">
<mml:mrow>
<mml:mi mathvariant="normal">R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Ecosystem recovery capacity; <inline-formula id="inf7">
<mml:math id="m9">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Proportion of area occupied by land use type <italic>k</italic>; <inline-formula id="inf8">
<mml:math id="m10">
<mml:mrow>
<mml:mi>R</mml:mi>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Recovery coefficient for land use type <italic>k</italic>.</p>
<p>Ecological Resilience (<inline-formula id="inf9">
<mml:math id="m11">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>): Due to differing units of measurement for resistance (<inline-formula id="inf10">
<mml:math id="m12">
<mml:mrow>
<mml:mi>P</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>), adaptive capacity (<inline-formula id="inf11">
<mml:math id="m13">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>), and recovery capacity (<inline-formula id="inf12">
<mml:math id="m14">
<mml:mrow>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>), this study normalizes all three indicators to the range [0, 1] to precisely quantify regional ER levels. Following existing studies (<xref ref-type="bibr" rid="B36">Wang et al., 2025e</xref>), resilience levels are classified into five tiers: Low (&#x3c;0.15); Medium-Low (0.15&#x2013;0.30); Medium (0.30&#x2013;0.45); Medium-High (0.45&#x2013;0.60); High (&#x3e;0.60). See specifically <xref ref-type="disp-formula" rid="e3">Equation 3</xref>.<disp-formula id="e3">
<mml:math id="m15">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mroot>
<mml:mrow>
<mml:mi>P</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>A</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:mn>3</mml:mn>
</mml:mroot>
</mml:mrow>
</mml:math>
<label>(3)</label>
</disp-formula>
</p>
</sec>
<sec id="s2-3-2">
<title>2.3.2 LER assessment</title>
<p>Landscape disturbance degree (<inline-formula id="inf13">
<mml:math id="m16">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) quantifies the comprehensive intensity of external disturbances on ecosystems through the integration of multidimensional indicators, with its value demonstrating a significant positive correlation with regional ecological risk. This metric is constructed based on three components: fragmentation (<inline-formula id="inf14">
<mml:math id="m17">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), separation (<inline-formula id="inf15">
<mml:math id="m18">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), and dominance (<inline-formula id="inf16">
<mml:math id="m19">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>). Following established methodologies (<xref ref-type="bibr" rid="B12">Li R. et al., 2025</xref>), their weighting coefficients are assigned as 0.5, 0.3, and 0.2, respectively. See specifically <xref ref-type="disp-formula" rid="e4">Equations 4</xref>&#x2013;<xref ref-type="disp-formula" rid="e7">7</xref>.<disp-formula id="e4">
<mml:math id="m20">
<mml:mrow>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mn>0.5</mml:mn>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.3</mml:mn>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:mn>0.2</mml:mn>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(4)</label>
</disp-formula>
<disp-formula id="e5">
<mml:math id="m21">
<mml:mrow>
<mml:msub>
<mml:mi>C</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(5)</label>
</disp-formula>
<disp-formula id="e6">
<mml:math id="m22">
<mml:mrow>
<mml:msub>
<mml:mi>N</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mi>A</mml:mi>
<mml:mrow>
<mml:mn>2</mml:mn>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
<mml:msqrt>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:mfrac>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(6)</label>
</disp-formula>
<disp-formula id="e7">
<mml:math id="m23">
<mml:mrow>
<mml:msub>
<mml:mi>D</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2b;</mml:mo>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mn>4</mml:mn>
</mml:mfrac>
<mml:mo>&#x2b;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:msub>
<mml:mi>L</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
<mml:mrow>
<mml:mn>2</mml:mn>
</mml:mrow>
</mml:mfrac>
</mml:mrow>
</mml:math>
<label>(7)</label>
</disp-formula>Where: <inline-formula id="inf17">
<mml:math id="m24">
<mml:mrow>
<mml:msub>
<mml:mi>n</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Number of patches of landscape type <italic>i</italic>; <inline-formula id="inf18">
<mml:math id="m25">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Total area of landscape type <italic>i</italic>; <inline-formula id="inf19">
<mml:math id="m26">
<mml:mrow>
<mml:mi>A</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Total landscape area; <inline-formula id="inf20">
<mml:math id="m27">
<mml:mrow>
<mml:msub>
<mml:mi>Q</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Quadrats frequency of landscape type <italic>i</italic>; <inline-formula id="inf21">
<mml:math id="m28">
<mml:mrow>
<mml:msub>
<mml:mi>M</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Proportion of patches of landscape type <italic>i</italic> to total patches.</p>
<p>Landscape Vulnerability quantifies the susceptibility of distinct landscape types to degradation or functional loss under external disturbances, reflecting their inherent structural fragility and resistance capacity. Building on regional characteristics and established methodologies (<xref ref-type="bibr" rid="B32">Wang et al., 2025a</xref>), we implemented a Delphi informed weighting protocol to classify landscape vulnerability into six hierarchical tiers (highest to lowest). Tier scores were subsequently normalized via min-max scaling to derive standardized vulnerability indices (<inline-formula id="inf22">
<mml:math id="m29">
<mml:mrow>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>), with full classification detailed in <xref ref-type="table" rid="T3">Table 3</xref>.</p>
<table-wrap id="T3" position="float">
<label>TABLE 3</label>
<caption>
<p>Landscape vulnerability.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Landscape types</th>
<th align="center">Cropland</th>
<th align="center">Forestland</th>
<th align="center">Grassland</th>
<th align="center">Water bodies</th>
<th align="center">Built-up land</th>
<th align="center">Unused land</th>
<th align="center">Landscape types</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">Assignment</td>
<td align="center">4</td>
<td align="center">2</td>
<td align="center">3</td>
<td align="center">5</td>
<td align="center">1</td>
<td align="center">6</td>
<td align="center">Assignment</td>
</tr>
<tr>
<td align="center">Fragility Index</td>
<td align="center">0.1905</td>
<td align="center">0.0952</td>
<td align="center">0.1429</td>
<td align="center">0.2381</td>
<td align="center">0.0476</td>
<td align="center">0.2857</td>
<td align="center">Fragility Index</td>
</tr>
</tbody>
</table>
</table-wrap>
<p>The landscape loss index (<inline-formula id="inf23">
<mml:math id="m30">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula>) evaluates the attenuation of ecological attributes under natural anthropogenic compound effects, calculated by integrating disturbance degree and vulnerability metrics (<xref ref-type="bibr" rid="B37">Wu et al., 2025</xref>), see specifically <xref ref-type="disp-formula" rid="e8">Equation 8</xref>.<disp-formula id="e8">
<mml:math id="m31">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:msub>
<mml:mi>E</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#xd7;</mml:mo>
<mml:msub>
<mml:mi>F</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(8)</label>
</disp-formula>
</p>
<p>The ecological risk index (<inline-formula id="inf24">
<mml:math id="m32">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>) is constructed using the aforementioned indices (<xref ref-type="bibr" rid="B15">Liang et al., 2022</xref>; <xref ref-type="bibr" rid="B23">Qu et al., 2022</xref>). Based on regional characteristics, ERI is classified into five risk levels: Low ecological risk [<inline-formula id="inf25">
<mml:math id="m33">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.03]; Moderately low ecological risk [0.03 &#x2264; <inline-formula id="inf26">
<mml:math id="m34">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.045]; Medium ecological risk [0.045 &#x2264; <inline-formula id="inf27">
<mml:math id="m35">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.06]; Moderately high ecological risk [0.06 &#x2264; <inline-formula id="inf28">
<mml:math id="m36">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3c;0.1]; High ecological risk [0.1 &#x3c; <inline-formula id="inf29">
<mml:math id="m37">
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula>]. See specifically <xref ref-type="disp-formula" rid="e9">Equation 9</xref>.<disp-formula id="e9">
<mml:math id="m38">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x3d;</mml:mo>
<mml:mstyle displaystyle="true">
<mml:munderover>
<mml:mo>&#x2211;</mml:mo>
<mml:mi>i</mml:mi>
<mml:mi>N</mml:mi>
</mml:munderover>
</mml:mstyle>
<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:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(9)</label>
</disp-formula>Where: <inline-formula id="inf30">
<mml:math id="m39">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Ecological Risk Index of risk zone <italic>i</italic>; <inline-formula id="inf31">
<mml:math id="m40">
<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:math>
</inline-formula> &#x3d; Area of landscape type <italic>i</italic> within risk zone <italic>k</italic>; <inline-formula id="inf32">
<mml:math id="m41">
<mml:mrow>
<mml:msub>
<mml:mi>A</mml:mi>
<mml:mi>k</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Total area of risk zone <italic>k</italic>; <inline-formula id="inf33">
<mml:math id="m42">
<mml:mrow>
<mml:msub>
<mml:mi>R</mml:mi>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Loss index of landscape type <italic>i</italic>.</p>
</sec>
<sec id="s2-3-3">
<title>2.3.3 Correlation analysis</title>
<p>To systematically elucidate interaction mechanisms among LER, ER, and their subsystems, we established a statistical verification framework for variable relationships. Given the characteristics of continuous variable data, Pearson correlation coefficient analysis was employed to test the significance of bivariate linear associations. This method effectively achieves unit invariance through covariance normalization, rendering it superior to alternative correlation techniques for linear relationship quantification in continuously distributed geospatial data.</p>
</sec>
<sec id="s2-3-4">
<title>2.3.4 Correlation analysis</title>
<p>This study employs bivariate spatial autocorrelation to analyze the spatial association between LER and ER (including its subsystems). The spatial dependency is revealed through global Moran&#x2019;s I calculation and local Moran&#x2019;s scatter plot mapping. The global Moran&#x2019;s I (ranging from &#x2212;1 to 1) quantifies the strength of spatial association between variables: Positive values indicate spatial clustering of similar attributes (e.g., adjacent areas exhibit concurrent high/low ER and risk levels; values approaching 1 denote significant spatial correlation). Negative values suggest spatial divergence (e.g., high resilience areas adjacent to low-risk zones). Zero values represent random spatial distribution. The local Moran&#x2019;s scatter plot categorizes spatial clusters into four types: High-High (H-H): High ER with high LER. High-Low (H-L): High ER with low LER. Low-Low (L-L): Low ER with low LER. Low-High (L-H): Low ER with high LER. The calculation formula follows reference (<xref ref-type="bibr" rid="B9">Lee, 2001</xref>).</p>
</sec>
<sec id="s2-3-5">
<title>2.3.5 Coupling coordination model</title>
<p>LER quantifies potential adverse consequences to ecosystems arising from natural or anthropogenic disturbances, whereas ER underscores a system&#x2019;s capacity to adapt, recover, and reorganize following such disruptions. Both concepts probe the stability and dynamic response mechanisms of ecosystems, exhibiting multifaceted intrinsic interconnections. This study investigates the reciprocal relationship between regional scale LER and ER by implementing a modified coupling coordination model. See specifically <xref ref-type="disp-formula" rid="e10">Equations 10</xref>&#x2013;<xref ref-type="disp-formula" rid="e12">12</xref>.<disp-formula id="e10">
<mml:math id="m43">
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mrow>
<mml:mfenced open="[" close="]" separators="|">
<mml:mrow>
<mml:mn>1</mml:mn>
<mml:mo>&#x2212;</mml:mo>
<mml:msqrt>
<mml:msup>
<mml:mrow>
<mml:mfenced open="(" close=")" separators="|">
<mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
<mml:mo>&#x2212;</mml:mo>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mn>2</mml:mn>
</mml:msup>
</mml:msqrt>
</mml:mrow>
</mml:mfenced>
</mml:mrow>
<mml:mo>&#xd7;</mml:mo>
<mml:mfrac>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
</mml:mrow>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mfrac>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(10)</label>
</disp-formula>
<disp-formula id="e11">
<mml:math id="m44">
<mml:mrow>
<mml:mi>T</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:mi>&#x3b1;</mml:mi>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mo>&#x2b;</mml:mo>
<mml:mi>&#x3b2;</mml:mi>
<mml:msub>
<mml:mrow>
<mml:mi>E</mml:mi>
<mml:mi>R</mml:mi>
<mml:mi>I</mml:mi>
</mml:mrow>
<mml:mi>i</mml:mi>
</mml:msub>
</mml:mrow>
</mml:math>
<label>(11)</label>
</disp-formula>
<disp-formula id="e12">
<mml:math id="m45">
<mml:mrow>
<mml:mi>D</mml:mi>
<mml:mo>&#x3d;</mml:mo>
<mml:msqrt>
<mml:mrow>
<mml:mi>C</mml:mi>
<mml:mo>&#xd7;</mml:mo>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:msqrt>
</mml:mrow>
</mml:math>
<label>(12)</label>
</disp-formula>Where: <inline-formula id="inf34">
<mml:math id="m46">
<mml:mrow>
<mml:mi>C</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Coupling Degree; <inline-formula id="inf35">
<mml:math id="m47">
<mml:mrow>
<mml:mi>T</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Coordination Degree; <inline-formula id="inf36">
<mml:math id="m48">
<mml:mrow>
<mml:mi>D</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Coupling Coordination Degree; <inline-formula id="inf37">
<mml:math id="m49">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> and <inline-formula id="inf38">
<mml:math id="m50">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; Contribution coefficients. Given the equal importance of LER and ER in this study, we set <inline-formula id="inf39">
<mml:math id="m51">
<mml:mrow>
<mml:mi>&#x3b1;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; <inline-formula id="inf40">
<mml:math id="m52">
<mml:mrow>
<mml:mi>&#x3b2;</mml:mi>
</mml:mrow>
</mml:math>
</inline-formula> &#x3d; 0.5. Following reference (<xref ref-type="bibr" rid="B39">Xu F et al., 2025</xref>), coupling coordination types are classified into 10 levels (see <xref ref-type="table" rid="T4">Table 4</xref>).</p>
<table-wrap id="T4" position="float">
<label>TABLE 4</label>
<caption>
<p>Coupling coordination types.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">D-value range</th>
<th align="center">Level</th>
<th align="center">D-value range</th>
<th align="center">Level</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">[0.0,0.1)</td>
<td align="center">Extreme disharmony</td>
<td align="center">[0.5,0.6)</td>
<td align="center">Weak coordination</td>
</tr>
<tr>
<td align="center">[0.1,0.2)</td>
<td align="center">Serious disharmony</td>
<td align="center">[0.6,0.7)</td>
<td align="center">Basic coordination</td>
</tr>
<tr>
<td align="center">[0.2,0.3)</td>
<td align="center">Moderate disharmony</td>
<td align="center">[0.7,0.8)</td>
<td align="center">Moderate coordination</td>
</tr>
<tr>
<td align="center">[0.3,0.4)</td>
<td align="center">Mild disharmony</td>
<td align="center">[0.8,0.9)</td>
<td align="center">Effective alignment</td>
</tr>
<tr>
<td align="center">[0.4,0.5)</td>
<td align="center">Near-disharmony</td>
<td align="center">[0.9,1.0]</td>
<td align="center">Optimal coordination</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<sec id="s3-1">
<title>3.1 Spatiotemporal evolution of ER</title>
<p>To investigate the spatiotemporal evolution of ER and its subsystems within the Hefei Metropolitan Circle during 2010&#x2013;2020, <xref ref-type="fig" rid="F3">Figure 3</xref> was constructed. This analysis reveals pronounced spatiotemporal heterogeneity in regional ER and its components. Regarding resistance, a &#x201c;high-periphery, low-core&#x201d; spatial pattern persisted throughout the study period, with over 93% of grid units consistently classified as medium-low or lower levels, indicating minimal temporal variation. Conversely, adaptability exhibited diametrically opposed characteristics, maintaining high performance across all timeframes; grids with medium-high or higher adaptability exceeded 57% in all three temporal dimensions. The spatial configuration of recovery closely mirrored that of resistance, yet demonstrated a persistent decline: the proportion of grids with medium-high or higher recovery decreased from 27.04% (2010) to 27.00% (2015) and 26.24% (2020), signaling a gradual but sustained deterioration. Spatially, ER hotspots predominantly clustered in the northern sector, while over 60% of the region remained below medium-high resilience thresholds. Temporally, <xref ref-type="table" rid="T5">Table 5</xref> underscores a systematic transition in ER dynamics over the decade. Deconstructing resilience hierarchy reveals critical shifts: Low resilience zones displayed dual expansion a 6.7% net area increase (&#x2b;0.67 percentage points in spatial coverage) evidencing intensified pressure from rapid urbanization on ecologically vulnerable areas. Medium-low resilience zones contracted by 1.97 percentage points yet maintained dominant coverage (60.67%), reflecting inherent regional ecological fragility. Medium-resilience zones registered concurrent growth in area (&#x2b;5.5%) and spatial proportion (&#x2b;1.21%). High-resilience zones suffered precipitous declines, with area plummeting by 50.6% and spatial share diminishing to 0.54%. These findings collectively demonstrate suboptimal regional ER levels, necessitating urgent spatially differentiated governance frameworks tailored to local contexts to enhance holistic ER.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Spatiotemporal evolution of ER and its subsystems in Hefei Metropolitan Circle (2010-2020): <bold>(a1-a3)</bold> Resistance Dimension; <bold>(b1-b3)</bold> Adaptability Dimension; <bold>(c1-c3)</bold> Resilience Dimension; <bold>(d1-d3)</bold> Ecological Resilience.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g003.tif">
<alt-text content-type="machine-generated">A grid of twelve maps shows spatial data changes from 2010 to 2020 across four categories: P, A, R, and ER. Each row represents a category, and each column represents a year. The intensity of shading on the maps varies, indicating different levels of data, with a color scale on the right indicating low to high values. Each map includes a compass rose and a scale bar.</alt-text>
</graphic>
</fig>
<table-wrap id="T5" position="float">
<label>TABLE 5</label>
<caption>
<p>Area and proportion of ER levels in Hefei Metropolitan Circle from 2010 to 2020.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Year</th>
<th colspan="2" align="center">Low level</th>
<th colspan="2" align="center">Medium-low level</th>
<th colspan="2" align="center">Medium level</th>
<th colspan="2" align="center">Medium-high level</th>
<th colspan="2" align="center">High level</th>
</tr>
<tr>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">2010</td>
<td align="center">641,227</td>
<td align="center">10.11%</td>
<td align="center">3,974,646</td>
<td align="center">62.64%</td>
<td align="center">1,394,222</td>
<td align="center">21.97%</td>
<td align="center">293,116</td>
<td align="center">4.62%</td>
<td align="center">69,070</td>
<td align="center">1.09%</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">650,019</td>
<td align="center">10.24%</td>
<td align="center">3,902,220</td>
<td align="center">61.50%</td>
<td align="center">1,443,809</td>
<td align="center">22.75%</td>
<td align="center">307,163</td>
<td align="center">4.84%</td>
<td align="center">42,070</td>
<td align="center">0.66%</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">683,931</td>
<td align="center">10.78%</td>
<td align="center">3,849,737</td>
<td align="center">60.67%</td>
<td align="center">1,471,135</td>
<td align="center">23.18%</td>
<td align="center">306,363</td>
<td align="center">4.83%</td>
<td align="center">34,115</td>
<td align="center">0.54%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-2">
<title>3.2 Spatiotemporal evolution of LER</title>
<p>To further investigate the spatiotemporal evolution of LER in the Hefei Metropolitan Circle (2010&#x2013;2020), <xref ref-type="fig" rid="F4">Figure 4</xref> and <xref ref-type="table" rid="T6">Table 6</xref> were developed. Key findings reveal divergent temporal trajectories across risk hierarchies: Low-risk areas experienced continuous contraction, decreasing from 1.3488 million hectares (21.26%) in 2010 to 1.2853 million hectares (20.26%) in 2020, indicating progressive encroachment on natural ecological baselines. Moderately low risk areas consistently dominated the landscape, with coverage increasing marginally from 68.70% to 69.53%, reflecting sustained stability in fundamental ecological conditions. Moderate risk areas exhibited minimal fluctuation, rising only 0.12 percentage points (4.63%&#x2013;4.75%), suggesting localized ecological pressures remain manageable. Moderately high risk and high risk areas demonstrated counteracting trends: the former declined from 3.98% to 3.91%, while the latter increased from 1.43% to 1.55%, revealing heightened vulnerability in zones impacted by intensive anthropogenic activities. Spatially, low risk zones concentrated in the southwestern region with higher forest coverage, whereas high-risk areas clustered around water bodies. Overall, the ecological risk pattern demonstrated a &#x201c;stable baseline with localized differentiation,&#x201d; necessitating vigilance against the expansion of high risk zones and their potential threats to ER.</p>
<fig id="F4" position="float">
<label>FIGURE 4</label>
<caption>
<p>Spatiotemporal evolution of LER in the Hefei Metropolitan Area (2010&#x2013;2020).</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g004.tif">
<alt-text content-type="machine-generated">Three maps of a region from 2010, 2015, and 2020 show levels of LER in varying shades of green, indicating low to high values. Each map includes a compass rose for orientation and a scale bar in kilometers. A legend on the right categorizes the shades as low, medium-low, medium, medium-high, and high.</alt-text>
</graphic>
</fig>
<table-wrap id="T6" position="float">
<label>TABLE 6</label>
<caption>
<p>Area and proportion of different LER levels in the Hefei Metropolitan Area (2010&#x2013;2020).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Year</th>
<th colspan="2" align="center">Low risk</th>
<th colspan="2" align="center">Medium-low risk</th>
<th colspan="2" align="center">Medium risk</th>
<th colspan="2" align="center">Medium-high risk</th>
<th colspan="2" align="center">High risk</th>
</tr>
<tr>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
<th align="center">Area/h&#x33a1;</th>
<th align="center">Percentage/%</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">2010</td>
<td align="center">1,348,838</td>
<td align="center">21.26%</td>
<td align="center">4,359,047</td>
<td align="center">68.70%</td>
<td align="center">293,935</td>
<td align="center">4.63%</td>
<td align="center">252,720</td>
<td align="center">3.98%</td>
<td align="center">90,741</td>
<td align="center">1.43%</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">1,294,065</td>
<td align="center">20.39%</td>
<td align="center">4,411,408</td>
<td align="center">69.52%</td>
<td align="center">292,642</td>
<td align="center">4.61%</td>
<td align="center">255,655</td>
<td align="center">4.03%</td>
<td align="center">91,511</td>
<td align="center">1.44%</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">1,285,313</td>
<td align="center">20.26%</td>
<td align="center">4,412,131</td>
<td align="center">69.53%</td>
<td align="center">301,301</td>
<td align="center">4.75%</td>
<td align="center">248,276</td>
<td align="center">3.91%</td>
<td align="center">98,260</td>
<td align="center">1.55%</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
<sec id="s3-3">
<title>3.3 Interactive relationship between LER and ER</title>
<sec id="s3-3-1">
<title>3.3.1 Correlation analysis of LER and ER</title>
<p>This study employed Pearson correlation analysis to reveal the multi-scale interactions between LER and ER (including its subsystems) in the Hefei Metropolitan Circle from 2010 to 2020 (<xref ref-type="fig" rid="F5">Figure 5</xref>). Temporally, the negative correlation between ER and LER intensified at both municipal and county scales during the study period. Notably, at the county scale, the correlation coefficient between LER and ER decreased cumulatively by 164.5% over the decade, indicating increasingly pronounced inhibitory effects of cumulative ecological risks on resilience systems at finer administrative scales. Across scales, recovery capacity consistently exhibited the strongest negative correlation with LER in all three study periods. At the municipal scale, correlation coefficients remained below &#x2212;0.72, while values at the county and grid scales were less than &#x2212;0.55 and &#x2212;0.13, respectively. Additionally, significant positive correlations were observed between ER and its subsystems, suggesting synergistic interactions within the resilience framework. These findings collectively underscore the hierarchical and scale dependent nature of risk-resilience dynamics, with recovery capacity emerging as a critical buffer against regional ecological risks.</p>
<fig id="F5" position="float">
<label>FIGURE 5</label>
<caption>
<p>Heatmap of multi-scale correlation between LER and ER from 2010 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g005.tif">
<alt-text content-type="machine-generated">Nine bubble plots show correlation matrices for different years (2010, 2015, 2020) and regions (City, County, Grid). Each plot displays the relationships between variables ER, P, A, R, and LER, with bubble sizes and colors representing correlation strength and direction. Positive correlations are indicated in green, negative in purple.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3-3-2">
<title>3.3.2 Bivariate spatial autocorrelation</title>
<p>
<xref ref-type="table" rid="T7">Table 7</xref> presents the bivariate Moran&#x2019;s I test results for the spatial autocorrelation between LER and ER (including its subsystems) in the Hefei Metropolitan Circle across multiple scales from 2010 to 2020. The analysis reveals that finer spatial scales amplify the significance of correlations between LER and ER metrics. At the municipal scale, only adaptive capacity exhibited a moderately significant positive correlation with LER (P &#x3c; 0.05),while other correlations remained weak (At the municipal scale, the sample size falls below the empirically recommended value for general spatial autocorrelation analysis, and this result may lack stability.). In contrast, at the county scale, recovery capacity demonstrated a strongly significant negative correlation with LER (P &#x3c; 0.001, Z &#x3c;&#x2212;72),whereas adaptive capacity mirrored this pattern with an extremely significant positive correlation (P &#x3c; 0.001, Z &#x3e; 10). At the grid scale, all subsystems of ER showed highly significant correlations (P &#x3c; 0.001), with LER:ER and recovery capacity exhibited significant negative correlations; Resistance and adaptive capacity displayed significant positive correlations. Notably, the negative correlation between ER and LER at the grid scale intensified annually over the study period. These findings highlight the critical role of spatial granularity in uncovering complex risk-resilience interactions, with grid-scale analysis providing the highest sensitivity to localized ecological dynamics.</p>
<table-wrap id="T7" position="float">
<label>TABLE 7</label>
<caption>
<p>Results of bivariate Moran&#x2019;s I test for LER and ER (including its subsystems) across multiple scales from 2010 to 2020.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="center">Variable</th>
<th rowspan="2" align="center">Year</th>
<th colspan="3" align="center">Moran&#x2019;s I</th>
<th colspan="3" align="center">Z-score</th>
</tr>
<tr>
<th align="center">City-level scale</th>
<th align="center">County-level scale</th>
<th align="center">Grid scale</th>
<th align="center">City-level scale</th>
<th align="center">County-level scale</th>
<th align="center">Grid scale</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td rowspan="3" align="center">R</td>
<td align="center">2010</td>
<td align="center">&#x2212;0.0814</td>
<td align="center">&#x2212;0.2818<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.2708<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;1.1913</td>
<td align="center">&#x2212;4.2757</td>
<td align="center">&#x2212;86.6512</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">&#x2212;0.0739</td>
<td align="center">&#x2212;0.3044<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.2478<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;1.1812</td>
<td align="center">&#x2212;4.6369</td>
<td align="center">&#x2212;79.1775</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">&#x2212;0.0802</td>
<td align="center">&#x2212;0.2892<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.2227<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;1.1844</td>
<td align="center">&#x2212;4.3629</td>
<td align="center">&#x2212;72.8169</td>
</tr>
<tr>
<td rowspan="3" align="center">P</td>
<td align="center">2010</td>
<td align="center">&#x2212;0.1783</td>
<td align="center">0.0252</td>
<td align="center">0.2051<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;1.1048</td>
<td align="center">0.4452</td>
<td align="center">60.3251</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">&#x2212;0.1674</td>
<td align="center">0.0279</td>
<td align="center">0.2112<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.9702</td>
<td align="center">0.4962</td>
<td align="center">61.3981</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">&#x2212;0.2084</td>
<td align="center">0.0103</td>
<td align="center">0.2262<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;1.2601</td>
<td align="center">0.2682</td>
<td align="center">65.6308</td>
</tr>
<tr>
<td rowspan="3" align="center">A</td>
<td align="center">2010</td>
<td align="center">0.2391<sup>&#x2a;</sup>
</td>
<td align="center">0.2805<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">0.0759<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">1.9938</td>
<td align="center">4.2008</td>
<td align="center">24.6285</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">0.2817<sup>&#x2a;</sup>
</td>
<td align="center">0.3069<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">0.0589<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">2.1297</td>
<td align="center">4.5886</td>
<td align="center">19.195</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">0.2952<sup>&#x2a;</sup>
</td>
<td align="center">0.2909<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">0.0317<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">2.0764</td>
<td align="center">4.3936</td>
<td align="center">10.3695</td>
</tr>
<tr>
<td rowspan="3" align="center">ER</td>
<td align="center">2010</td>
<td align="center">&#x2212;0.0074</td>
<td align="center">&#x2212;0.0501</td>
<td align="center">&#x2212;0.0643<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.3903</td>
<td align="center">&#x2212;0.7586</td>
<td align="center">&#x2212;20.9869</td>
</tr>
<tr>
<td align="center">2015</td>
<td align="center">0.0202</td>
<td align="center">&#x2212;0.0612</td>
<td align="center">&#x2212;0.0661<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.2057</td>
<td align="center">&#x2212;0.9249</td>
<td align="center">&#x2212;21.3443</td>
</tr>
<tr>
<td align="center">2020</td>
<td align="center">0.0006</td>
<td align="center">&#x2212;0.0386</td>
<td align="center">&#x2212;0.068<sup>&#x2a;&#x2a;&#x2a;</sup>
</td>
<td align="center">&#x2212;0.4233</td>
<td align="center">&#x2212;1.2736</td>
<td align="center">&#x2212;22.0043</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn>
<p>&#x2a;&#x2a; &#x3d; P &#x3c; 0.05,&#x2a;&#x2a; &#x3d; P &#x3c; 0.01,&#x2a;&#x2a;&#x2a; &#x3d; P &#x3c; 0.001.</p>
</fn>
</table-wrap-foot>
</table-wrap>
<p>To further analyze the spatial autocorrelation between LER and ER (including its subsystems) across multiple scales, this study presents <xref ref-type="fig" rid="F6">Figures 6</xref>&#x2013;<xref ref-type="fig" rid="F9">9</xref>. A comparative analysis across different scales reveals that as the research scale becomes more refined, the spatial correlation between LER and ER (and its subsystems) gradually strengthens, and the spatial clustering effect becomes more pronounced. The bivariate Moran&#x2019;s index between LER and ER at both the county and grid scales shows a gradual increase over time, indicating a continuous enhancement of their negative spatial correlation. Among the subsystems of ecological resilience, only restorative capacity (R) consistently exhibits a certain degree of negative spatial correlation with LER across all three scales.</p>
<fig id="F6" position="float">
<label>FIGURE 6</label>
<caption>
<p>Bivariate Moran scatter plots and LISA maps of landscape ecological risk and resistance (P) in the Hefei Metropolitan Circle at multi-scales from 2010 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g006.tif">
<alt-text content-type="machine-generated">Three rows of maps and scatter plots show spatial and statistical analysis for city, county, and grid in 2010, 2015, and 2020. Each panel includes a scatter plot of LER versus ER and a map indicating spatial patterns. Colors show categories: high-high, low-low, and non-significant areas, with corresponding Moran&#x27;s I values.</alt-text>
</graphic>
</fig>
<fig id="F7" position="float">
<label>FIGURE 7</label>
<caption>
<p>Bivariate Moran scatter plots and LISA maps of landscape ecological risk and adaptability (A) in the Hefei Metropolitan Circle at multi-scales from 2010 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g007.tif">
<alt-text content-type="machine-generated">The image presents a series of scatter plots and maps arranged in a grid showing spatial autocorrelation of LER (Logarithm of Employment Rate) versus ER (Employment Rate) across different administrative levels&#x2014;City, County, and Grid&#x2014;over the years 2010, 2015, and 2020. Each row corresponds to an administrative level, and each column represents a year. Scatter plots depict Moran&#x27;s I values indicating spatial correlation, and maps illustrate clustering with color codes for significance: red for high-low, blue for low-high, pink for low-low, and gray for non-significant areas.</alt-text>
</graphic>
</fig>
<fig id="F8" position="float">
<label>FIGURE 8</label>
<caption>
<p>Bivariate Moran scatter plots and LISA maps of landscape ecological risk and resilience (R) in the Hefei Metropolitan Circle at multi-scales from 2010 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g008.tif">
<alt-text content-type="machine-generated">Grid of six scatter plots and spatial maps comparing changes in spatial relationships over three time periods: 2010, 2015, and 2020. Each row represents a different spatial level: City, County, and Grid. Scatter plots show ER vs. LER with a regression line and Moran&#x27;s I value. Maps highlight regions in varying colors, indicating &#x22;high-high,&#x22; &#x22;low-low,&#x22; &#x22;low-high,&#x22; and &#x22;non-significant&#x22; areas.</alt-text>
</graphic>
</fig>
<fig id="F9" position="float">
<label>FIGURE 9</label>
<caption>
<p>Bivariate Moran scatter plots and LISA maps of landscape ecological risk and ecological resilience in the Hefei Metropolitan Circle at multi-scales from 2010 to 2020.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g009.tif">
<alt-text content-type="machine-generated">Scatter plots and maps depicting spatial autocorrelation (Moran&#x27;s I) of economic variables across city, county, and grid levels for years 2010, 2015, and 2020. Each row represents a level, showing changes in low-high and high-low significance, coded in blue and red, across geographic areas. The plots display the relationship between economic ratios (ER) and logarithm of economic ratio (LER) alongside geographic maps highlighting significant areas. Moran&#x2019;s I values indicate the degree of spatial clustering.</alt-text>
</graphic>
</fig>
<p>From a spatial distribution perspective, <xref ref-type="fig" rid="F6">Figures 6</xref>, <xref ref-type="fig" rid="F7">7</xref>, <xref ref-type="fig" rid="F9">9</xref> exhibit highly consistent spatial characteristics: At the city scale, the low-high clusters are consistently distributed within Chuzhou City, while the remaining areas show no significant spatial correlation. At the county scale, the low-low clusters are primarily located in the western part of the Hefei Metropolitan Area, including Huoqiu County, Yu&#x2019;an District, Yeji District, and Jin&#x2019;an District. The high-low clusters align with the distribution patterns observed at the grid scale, mainly concentrated in the southwestern part of the metropolitan area, encompassing counties such as Jinzhai, Huoshan, and Shucheng. At the grid scale, the high-high clusters are predominantly distributed in the northern and southeastern regions of the Hefei Metropolitan Area, while the low-high clusters are primarily aggregated around the Chaohu Lake watershed. <xref ref-type="fig" rid="F8">Figure 8</xref> reveals the spatial autocorrelation between adaptive capacity (A) and LER: At the city scale, no significant spatial correlation is observed. At the county scale, low-low clusters are distributed in the southwestern counties of the metropolitan area. At the grid scale, high-high clusters are sporadically scattered across the northern and eastern regions, as well as the periphery of Chaohu Lake. The low-low clusters align with the high-low clusters identified in the spatial autocorrelation analysis between other ecological resilience subsystems and landscape ecological risk, primarily located in the southwestern forested areas with high vegetation coverage.</p>
</sec>
</sec>
<sec id="s3-4">
<title>3.4 Coupling coordination analysis between LER and ER and its subsystems</title>
<p>
<xref ref-type="fig" rid="F10">Figure 10</xref> presents box plots illustrating the coupling coordination relationships between LER and ER (including its subsystems) in the Hefei Metropolitan Circle across multiple scales from 2010 to 2020. The results demonstrate that the mean coupling coordination degree between LER and ER remained below 0.5&#xa0;at all scales, indicating a marginally coordinated state. Temporally, the coupling coordination across all three scales exhibited a slight downward trend, reflecting weakened synergies over the study period. Spatially, municipal-scale coupling coordination degrees showed outliers approaching zero, suggesting minimal regional heterogeneity, while the county scale exhibited the highest number of outliers, highlighting significant inter-county disparities. In contrast, grid-scale outliers were relatively smaller, emphasizing localized consistency in risk-resilience dynamics. These findings collectively indicate a persistent trade-off relationship between LER and ER in the Hefei Metropolitan Circle during the study period, where improvements in one metric often coincided with deteriorations in the other. This antagonistic interaction underscores the necessity for spatially differentiated governance strategies to address scale-specific ecological challenges.</p>
<fig id="F10" position="float">
<label>FIGURE 10</label>
<caption>
<p>Box plots of coupling coordination relationships between LER and ER (including its subsystems) at multiple scales (2010&#x2013;2020).</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g010.tif">
<alt-text content-type="machine-generated">Three panels of boxplots showing data from 2010, 2015, and 2020. Each panel contains boxplots labeled P, A, R, and ER. The plots display different distributions and outliers with varying medians and spreads within each year, highlighting shifts over time.</alt-text>
</graphic>
</fig>
<p>
<xref ref-type="table" rid="T8">Table 8</xref> and <xref ref-type="fig" rid="F11">Figure 11</xref> present the spatial proportions and distribution of coupling coordination levels between LER and ER in the Hefei Metropolitan Circle across multiple scales from 2010 to 2020. The results reveal that, except for minor high-coupling-coordination areas at the municipal scale, county and grid scales exhibited predominantly low coordination levels. Spatially, the municipal scale was dominated by a &#x201c;mildly uncoordinated&#x201d; state, while &#x201c;critically uncoordinated&#x201d; states prevailed at county and grid scales. Temporal analysis highlighted significant divergence in coordination evolution across scales. At the municipal scale, the proportion of &#x201c;extremely uncoordinated&#x201d; areas remained unchanged at 18.08% over the decade, indicating persistent structural ecological challenges likely constrained by rigid macro-planning or systemic governance delays. A sudden surge in &#x201c;severely uncoordinated&#x201d; areas to 6.40% in 2015 (absent in other years) suggested potential transient ecological disturbances or data anomalies during this period. At the county scale, the proportion of &#x201c;extremely uncoordinated&#x201d; zones decreased from 2.25% to 0.82%, reflecting partial alleviation of high-risk pressures through localized governance. However, the &#x201c;mildly uncoordinated&#x201d; proportion rose from 14.45% to 18.97%, signaling enhanced diffusion of medium-low risks, potentially linked to marginal effects of land use transitions. The grid scale demonstrated a &#x201c;polarized stabilization&#x201d; pattern: the &#x201c;critically uncoordinated&#x201d; proportion increased from 48.04% to 51.75%, while the &#x201c;marginally coordinated&#x201d; proportion declined from 44.92% to 41.81%, revealing subtle degradation in small-scale coordination likely associated with ecological fragmentation driven by urban expansion. Additionally, the &#x201c;primary coordination&#x201d; level showed no significant improvement across scales, suggesting ER-building remains confined to passive adaptation. Overall, the multi-scale evolution of coupling coordination unveiled a &#x201c;macro-rigidity, meso-improvement, micro-fluctuation&#x201d; pattern, highlighting the urgent need for cross-scale collaborative governance mechanisms to address spatially heterogeneous challenges of risk transmission and resilience dissipation.</p>
<table-wrap id="T8" position="float">
<label>TABLE 8</label>
<caption>
<p>Percentage distribution of coupling coordination levels between LER and ER in the Hefei Metropolitan Circle (2010&#x2013;2020) under multi-scale perspectives (Unit: %).</p>
</caption>
<table>
<thead valign="top">
<tr>
<th rowspan="2" align="left">Coupling coordination degree</th>
<th colspan="3" align="left">City-level scale</th>
<th colspan="3" align="left">County-level scale</th>
<th colspan="3" align="left">Grid scale</th>
</tr>
<tr>
<th align="left">2010</th>
<th align="left">2015</th>
<th align="left">2020</th>
<th align="left"/>
<th align="left">2010</th>
<th align="left">2015</th>
<th align="left">2020</th>
<th align="left"/>
<th align="left">2010</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="left">Extreme disharmony</td>
<td align="left">18.08%</td>
<td align="left">18.08%</td>
<td align="left">18.08%</td>
<td align="left">2.25%</td>
<td align="left">2.41%</td>
<td align="left">0.82%</td>
<td align="left">0.94%</td>
<td align="left">0.87%</td>
<td align="left">0.83%</td>
</tr>
<tr>
<td align="left">Serious disharmony</td>
<td align="left">0.00%</td>
<td align="left">6.40%</td>
<td align="left">0.00%</td>
<td align="left">0.47%</td>
<td align="left">0.00%</td>
<td align="left">0.00%</td>
<td align="left">0.00%</td>
<td align="left">0.00%</td>
<td align="left">0.00%</td>
</tr>
<tr>
<td align="left">Moderate disharmony</td>
<td align="left">17.99%</td>
<td align="left">17.99%</td>
<td align="left">17.99%</td>
<td align="left">0.35%</td>
<td align="left">2.14%</td>
<td align="left">1.59%</td>
<td align="left">0.02%</td>
<td align="left">0.03%</td>
<td align="left">0.09%</td>
</tr>
<tr>
<td align="left">Mild disharmony</td>
<td align="left">6.40%</td>
<td align="left">21.29%</td>
<td align="left">27.70%</td>
<td align="left">14.45%</td>
<td align="left">19.20%</td>
<td align="left">18.97%</td>
<td align="left">0.47%</td>
<td align="left">0.54%</td>
<td align="left">0.74%</td>
</tr>
<tr>
<td align="left">Near-disharmony</td>
<td align="left">21.29%</td>
<td align="left">11.89%</td>
<td align="left">2.44%</td>
<td align="left">56.13%</td>
<td align="left">57.05%</td>
<td align="left">59.43%</td>
<td align="left">48.04%</td>
<td align="left">47.61%</td>
<td align="left">51.75%</td>
</tr>
<tr>
<td align="left">Weak coordination</td>
<td align="left">11.89%</td>
<td align="left">0.00%</td>
<td align="left">9.45%</td>
<td align="left">16.95%</td>
<td align="left">9.80%</td>
<td align="left">9.80%</td>
<td align="left">44.92%</td>
<td align="left">44.83%</td>
<td align="left">41.81%</td>
</tr>
<tr>
<td align="left">Basic coordination</td>
<td align="left">24.34%</td>
<td align="left">24.34%</td>
<td align="left">24.34%</td>
<td align="left">9.40%</td>
<td align="left">9.40%</td>
<td align="left">9.40%</td>
<td align="left">5.62%</td>
<td align="left">6.12%</td>
<td align="left">4.78%</td>
</tr>
<tr>
<td align="left">Extreme disharmony</td>
<td align="left">18.08%</td>
<td align="left">18.08%</td>
<td align="left">18.08%</td>
<td align="left">2.25%</td>
<td align="left">2.41%</td>
<td align="left">0.82%</td>
<td align="left">0.94%</td>
<td align="left">0.87%</td>
<td align="left">0.83%</td>
</tr>
</tbody>
</table>
</table-wrap>
<fig id="F11" position="float">
<label>FIGURE 11</label>
<caption>
<p>Spatial distribution of coupling coordination levels between LER and ER in the Hefei Metropolitan Circle from 2010 to 2020 under multi-scale perspectives.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g011.tif">
<alt-text content-type="machine-generated">Three rows of maps show regional coordination levels from 2010 to 2020 at city, county, and grid scales. Colors range from dark green indicating optimal coordination to light yellow indicating extreme disharmony. Each row corresponds to different administrative divisions, illustrating spatial changes over time.</alt-text>
</graphic>
</fig>
</sec>
<sec id="s3-5">
<title>3.5 Zoning regulation under multi-scale perspectives</title>
<p>Under multi-scale perspectives, this study employed bivariate local Moran&#x2019;s indices to delineate the spatial interactions between LER and ER (<xref ref-type="fig" rid="F9">Figure 9</xref>). By further integrating coupling coordination analysis (<xref ref-type="fig" rid="F11">Figure 11</xref>), the Hefei Metropolitan Circle was classified into distinct zoning units across three scales (<xref ref-type="fig" rid="F12">Figure 12</xref>). Synthesizing multi-scale spatial interactions and coupling coordination relationships, the metropolitan circle was categorized into four zones: high-resilience zones, high-risk hotspots, risk-resilience uncoordinated zones, and risk-resilience coordinated zones.</p>
<fig id="F12" position="float">
<label>FIGURE 12</label>
<caption>
<p>Zoning management of the Hefei Metropolitan Circle from a multi-scale perspective. <bold>(A)</bold> City-level Scale; <bold>(B)</bold> County-level Scale; <bold>(C)</bold> Grid Scale; <bold>(D)</bold> Comprehensive Scale.</p>
</caption>
<graphic xlink:href="fenvs-13-1654175-g012.tif">
<alt-text content-type="machine-generated">Four-panel map showing regional phase coordination and disharmony. Panel A: City-level map with areas marked as phase dissonant or coordinated. Panel B: County map showing phase states with detailed boundaries. Panel C: Grid map illustrating high ER and LER areas. Panel D: Comprehensive map displaying coordination and mismatch zones between ER and LER. Each panel includes a compass rose and scale bar. Legends describe color codings for different phases and risk levels.</alt-text>
</graphic>
</fig>
<p>From a coupling coordination perspective, only Lu&#x2019;an City and Wuhu City exhibited coordinated risk-resilience relationships at the municipal scale, while all other counties remained uncoordinated. At the county scale, coordinated relationships were observed in Nanqiao District, Hanshan County, Fanchang District, and Nanling County, with Jinzhai, Huoshan, and Shucheng Counties representing high-resilience/low-risk coordination. In contrast, grid-scale analysis revealed sporadic coordination in the southwestern and fragmented peripheral areas. Overall, the metropolitan circle demonstrated poor coupling coordination between LER and resilience across all scales. Spatial clustering analysis highlighted strong spatial dependency and negative correlations, particularly at the grid scale. Low-resilience/high-risk zones predominantly clustered in waterbody-rich regions such as Chaohu City and Huoqiu County, whereas high-resilience/low-risk zones concentrated in forest-dominated southwestern counties (Jinzhai, Huoshan, Shucheng). Low-resilience/low-risk and high-resilience/high-risk zones exhibited fragmented spatial distributions. These findings underscore the necessity for spatially adaptive governance strategies to address scale-specific ecological challenges and enhance systemic resilience.</p>
</sec>
</sec>
<sec id="s4">
<title>4 Policy and management implications</title>
<p>The spatial interplay between LER and ER exhibits considerable heterogeneity (<xref ref-type="fig" rid="F12">Figure 12D</xref>), necessitating comprehensive management strategies and region-specific mitigation measures informed by a spatial zoning scheme. Areas exhibiting high ER are discretely distributed adjacent to water bodies within the region, correlating with a low incidence of LERs and minimal anthropogenic disturbance. These zones demonstrate robust ecosystem stability and integrity. To optimize their ecological service value for socioeconomic development, strategic intensification of appropriately-scaled development within core protected areas is warranted. Concurrently, natural hydrological processes must be rigorously maintained. Implementing peripheral ecological buffer zones is critical to restrict agricultural non-point source pollution and impervious surface expansion. Furthermore, advancing integrated watershed management supported by ecological compensation mechanisms is essential to harmonize conservation with development. This holistic approach fosters the establishment of a continuous ecological corridor network, thereby enhancing regional ecological connectivity. LER hotspots are predominantly concentrated within the Chaohu Lake area and southeastern regions characterized by extensive water coverage. These ecologically critical zones exhibit pronounced land-use conflicts, where economic development coexists with conservation imperatives. Implementing robust regional coordination strategies is therefore essential. Concurrently, it is imperative to rigorously prohibit activities that encroach upon aquatic ecosystems, such as lake enclosure for agriculture (tian) and lake reclamation for development. Furthermore, regular LER assessments must be conducted to identify potential risk sources and vulnerable areas, enabling the formulation of targeted mitigation strategies. Finally, optimizing land-use structure to reduce unsustainable practices is crucial to mitigate adverse impacts on landscape ecological integrity. Regions exhibiting a mismatch between ER and LER are extensively distributed across non-forested and non-aquatic areas. These zones experience intensive anthropogenic pressure, characterized by high land development intensity driven by urbanization, diminished ER, and landscape configurations approaching critical thresholds of functional disruption. Stringent limitations on ecosystem disturbance are therefore imperative. Establishing a dynamic ecological risk assessment framework is essential, coupled with rigorous lifecycle ecological oversight of construction land. Regular analysis of monitoring data must be implemented to enable timely identification and resolution of emerging issues. Collectively, these measures will foster circular economic practices and advance green industry development. Synergistic zones between ER and LER are primarily concentrated in areas with high forest cover. Characterized by relatively low anthropogenic pressure, these regions exhibit elevated resilience levels and reduced ecological risk. To fully leverage their pivotal role as transitional buffers, stringent protection of core natural forests must be prioritized, including restrictions on commercial logging and large-scale development. Concurrently, bolstering ecological conservation policies should be integrated with leveraging forested landscapes to construct ecological corridor networks. These networks must strategically interface with adjacent farmlands and water bodies, fostering a balance in ecosystem service provision and demand across the broader landscape while significantly enhancing regional ecological connectivity.</p>
</sec>
<sec sec-type="discussion" id="s5">
<title>5 Discussion</title>
<sec id="s5-1">
<title>5.1 Spatiotemporal evolution of ER and LER</title>
<p>Our analysis reveals a significant shift in regional ER over the decade: low resilience zones expanded considerably (increasing by 17.3% coverage), while high resilience extremes contracted markedly (decreasing by 12.8% coverage) (<xref ref-type="table" rid="T5">Table 5</xref>). These shifts likely indicate heightened volatility in ecosystem resilience, suggesting progressive destabilization across specific dimensions. Crucially, this destabilization appears driven by cumulative effects of climate change, natural disasters, and intensified anthropogenic disturbances (<xref ref-type="bibr" rid="B4">Fan et al., 2024</xref>). As evidenced in <xref ref-type="fig" rid="F3">Figure 3d3</xref>, the Hefei Metropolitan Circle exhibited a marked expansion of low resilience zones alongside a significant contraction of high resilience areas in 2020. This spatial restructuring is likely attributable to synergistic pressures from extreme precipitation events and accelerated urbanization. Notably, the entire region particularly the Chaohu Lake Basin experienced historically unprecedented persistent rainfall. Water levels in Chaohu Lake surged to an all-time high (13.43&#xa0;m), sustaining prolonged exceedance above flood alert thresholds throughout the monitoring period. <xref ref-type="fig" rid="F3">Figure 3</xref> further reveals broad spatial convergence among the subsystems of ER, indicating strong coupling and synergistic interactions between these components (<xref ref-type="bibr" rid="B16">Liu et al., 2022</xref>). Consequently, policymakers must adopt integrated approaches for ecosystem management that account for all resilience subsystems. Focusing exclusively on isolated facets while neglecting systemic interdependencies risks compromising holistic resilience outcomes. Effective stewardship requires holistic and systematic strategies to safeguard and enhance the stability of ER.</p>
<p>Analysis of LER patterns revealed that moderate to low risk levels predominated across the Hefei Metropolitan Area over the decade-long study period (<xref ref-type="table" rid="T6">Table 6</xref>). Both risk categories exhibited marginal yet discernible increases, with high risk zones demonstrating a lower growth rate compared to their moderate to low counterparts. Spatially, low risk clusters were primarily concentrated in southwestern forested regions, whereas high risk areas coincided predominantly with water bodies. These spatial patterns broadly align with existing literature, providing robust evidence for regional risk dynamics. These patterns resonate with established regional observations. For instance, Li et al. documented pronounced spatial heterogeneity in LERs across Harbin over a 20-year period, identifying significant risk clustering proximal to water bodies (<xref ref-type="bibr" rid="B13">Li Y. et al., 2025</xref>). Similarly, scenario-based simulations by Ma et al. projected substantial expansion of high-risk zones in Hainan Island under multiple developmental trajectories (<xref ref-type="bibr" rid="B18">Ma et al., 2025</xref>). Collectively, these studies underscore the critical influence of hydrographic features on ecological risk agglomeration and validate the robustness of our analytical framework. This resilience may be attributed to targeted governance frameworks implemented by environmental authorities. Local authorities have executed sustained conservation initiatives across the Hefei Metropolitan Area over the past decade, with particular emphasis on the Chaohu Lake watershed. Notable interventions include the Ten Lakeside Wetlands Initiative launched in 2018 and the subsequent enactment of the Chaohu Wetlands Protection Statutes in 2019. Consequently, despite coinciding with China&#x2019;s most rapid urbanization phase, ecological degradation has been effectively contained&#x2014;demonstrating that strategic policy interventions can decouple economic expansion from environmental deterioration.</p>
</sec>
<sec id="s5-2">
<title>5.2 Interplay between ER and LER</title>
<p>Analysis of subsystem correlations revealed that resilience exhibited the most pronounced negative association with LER. This robust inverse relationship likely stems from the acute sensitivity of resilience metrics to ecosystem stability and integrity where heightened risk directly compromises buffering capacities essential for recovery. The robust negative association between ER (particularly its recovery metrics) and LER underscores the imperative of integrating risk mitigation as a core pillar in environmental governance frameworks. Effective conservation planning requires context-specific interventions that both suppress risk propagation and enhance systemic buffering capacities. The significant negative correlation between ER (including its subsystems) and LER further underscores the critical importance of incorporating LER mitigation as a core consideration in environmental conservation planning. Effective strategies must prioritize the enhancement of ecological resilience through regionally tailored adaptive measures. Critically, our analysis reveals a progressive intensification of the negative correlation between ecological resilience and LER with increasing spatial resolution. This scale-dependent amplification likely stems from inherent scaling effects (<xref ref-type="bibr" rid="B32">Wang C. et al., 2025</xref>). Consequently, environmental managers must explicitly account for spatial heterogeneity across administrative and ecological scales when formulating conservation policies or interventions. This multi-scale perspective is fundamental to ensuring both the effectiveness and long-term sustainability of implemented strategies.</p>
<p>Our analysis reveals a notably lower incidence of high coupling coordination values at the grid scale compared to broader spatial extents (<xref ref-type="table" rid="T8">Table 8</xref>). This pattern likely arises from the finer spatial resolution of grid-scale analysis, which more accurately captures localized heterogeneity and process dynamics (<xref ref-type="bibr" rid="B17">Liu et al., 2023</xref>). At broader spatial scales, our results indicate consistently higher coordination levels between landscape ecological risk and ER. This pattern likely reflects the combined influence of spatial averaging effects and the enhanced holistic stability inherent to integrated eco-system functioning. Elevated coupling coordination signifies greater ecosystem stability, where landscape pattern changes exert minimal disruptive influence on regional ER a critical determinant of ecological equilibrium (<xref ref-type="bibr" rid="B48">Zhang M. et al., 2023</xref>). Conversely, observed declines in coordination and localized decoupling likely result from the complex interplay of anthropogenic pressures (including resource overexploitation and pollution), climate change impacts, and endogenous ecosystem dynamics. Consequently, we urge policymakers to prioritize evidence-based conservation in future territorial planning by: Safeguarding ecologically sensitive zones and reinforcing critical habitat corridors through integrated ecological networks; Scaling up biodiversity protection initiatives and implementing robust ecological compensation mechanisms to minimize anthropogenic disturbance; Establishing cross-sectoral governance frame works to synchronize landscape configuration with resilience enhancement. his integrated approach is paramount for reconciling development pressures with long-term ecological integrity.</p>
</sec>
<sec id="s5-3">
<title>5.3 Limitations and future research directions</title>
<p>Compared to existing studies that predominantly emphasize theoretical analysis, this research introduces innovations in both its theoretical framework and analytical methodology: Firstly, by constructing a multi-scale analytical framework, it systematically reveals the spatiotemporal coupling relationships between LER and ER (including its subsystems), deepening the understanding of dynamic interaction mechanisms and coordinated evolution patterns among these elements. Secondly, moving beyond traditional theoretical limitations, it integrates practical zoning management with theoretical analysis. Based on the spatial heterogeneity and coupling coordination relationships of &#x201c;resilience-risk,&#x201d; it identifies which resilience subsystems (resistance, adaptation, or recovery) serve as critical leverage points for regulating regional ecological security. Furthermore, it proposes a spatially adaptive ecological management zoning scheme, providing a scientific basis and implementation pathway for optimizing regional ecological security patterns.</p>
<p>However, this study still has certain limitations. Urban ecological risks originate from complex sources, and the current ecological resilience assessment model fails to cover all dimensions of resilience. A more comprehensive evaluation system should be introduced in future research. The evolution of LER and ER exhibits nonlinear, long-term, and non-stationary characteristics, making a 10-year study period insufficient to fully capture their dynamics. Subsequent research could extend the study period and expand the scale range. Furthermore, the interaction mechanisms between the two warrant deeper exploration, such as multi-scale driving mechanisms and the regulatory role of landscape pattern optimization on resilience. More sophisticated models (e.g., dynamic simulation, agent-based modeling) could be employed for quantitative analysis. Future studies could also strengthen the following aspects: enhancing the analysis of socioeconomic driving mechanisms (e.g., land markets, planning policies) to reveal the profound impacts of human activities on ecosystems; conducting international comparative studies with metropolitan areas such as the Yangtze River Delta and Pearl River Delta to identify commonalities and differences, thereby improving the generalizability of conclusions; and developing predictive policy scenarios, such as simulating trends in landscape ecological risks around Chaohu Lake under different urban expansion rates, to propose preemptive management recommendations.</p>
</sec>
</sec>
<sec sec-type="conclusion" id="s6">
<title>6 Conclusion</title>
<p>
<list list-type="simple">
<list-item>
<p>1. The study revealed a significant negative correlation between ER and LER in the Hefei Metropolitan Area, with this negative correlation strengthening as the scale of analysis became more refined. Spatial heterogeneity was particularly pronounced at the grid scale, where the restorative capacity subsystem demonstrated the most prominent inhibitory effect on ecological risks. This finding deepens the understanding of the dynamic feedback mechanism between risk and resilience in highly urbanized areas and addresses a theoretical gap in multi-scale synergistic analysis.</p>
</list-item>
<list-item>
<p>2. This study innovatively identifies the spatial distribution and evolutionary patterns of key areas such as &#x201c;high-risk&#x2013;low-resilience&#x201d; and &#x201c;low-risk&#x2013;high-resilience&#x201d; across three scales (grid, county, and city levels). Based on spatial autocorrelation and coupling coordination analysis, it proposes a spatially adaptive zoning framework with four types of ecological management partitions. These findings break through the limitations of traditional single-scale studies, provide a scientific basis for cross-level and differentiated ecological governance, and represent a significant methodological and applied extension of existing ecological security theories.</p>
</list-item>
<list-item>
<p>3. Through a multi-scale coupling coordination evaluation, this study reveals that the average coupling coordination degree between ER and LER in the Hefei Metropolitan Area is below 0.5, indicating an overall state ranging from on the verge of dysfunction to barely coordinated. Over time, the coupling coordination degree across all three scales shows a slight declining trend. Spatially, coordinated areas are sporadically distributed in the southwestern forested regions, while dysfunctional areas are concentrated around water bodies and highly urbanized zones. At the city scale, core ecological issues exhibit long-term persistence; at the county scale, localized governance efforts show preliminary yet limited effectiveness, with a noticeable trend of medium to low risk diffusion; at the grid scale, coordination slightly deteriorates due to urban expansion.</p>
</list-item>
</list>
</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s7">
<title>Data availability statement</title>
<p>The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.</p>
</sec>
<sec sec-type="author-contributions" id="s8">
<title>Author contributions</title>
<p>JuL: Validation, Resources, Conceptualization, Methodology, Writing &#x2013; review and editing, Software, Writing &#x2013; original draft. JiL: Data curation, Investigation, Writing &#x2013; review and editing, Conceptualization, Validation, Project administration, Supervision, Funding acquisition, Writing &#x2013; original draft. YW: Writing &#x2013; original draft, Writing &#x2013; review and editing, Visualization, Validation.</p>
</sec>
<sec sec-type="funding-information" id="s9">
<title>Funding</title>
<p>The author(s) declare that financial support was received for the research and/or publication of this article. This study was funded by the Key Project of Anhui Province Philosophy and Social Science Planning: Research on the Construction of Historical Spatial Information Mapping of Traditional Villages in Huizhou and Their Protection Methods (AHSKD2023D028).</p>
</sec>
<ack>
<p>We would like to express our sincere gratitude to Anhui Jianzhu University for the assistance extended to our team.</p>
</ack>
<sec sec-type="COI-statement" id="s10">
<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="ai-statement" id="s11">
<title>Generative AI statement</title>
<p>The author(s) declare that no Generative AI was used in the creation of this manuscript.</p>
<p>Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.</p>
</sec>
<sec sec-type="disclaimer" id="s12">
<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>Dai</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Study on spatial correlation and driving factors of industrial carbon emission in the hefei metropolitan area based on SNA-QAP method</article-title>. <source>Sci. Rep.</source> <volume>15</volume>, <fpage>3445</fpage>. <pub-id pub-id-type="doi">10.1038/S41598-025-87805-Y</pub-id>
<pub-id pub-id-type="pmid">39870781</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wei</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Assessment of ecological resilience and its response mechanism to land spatial structure conflicts in China&#x2019;s Southeast Coastal Areas</article-title>. <source>Ecol. Indic.</source> <volume>170</volume>, <fpage>112980</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2024.112980</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Rong</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Tian</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Ou</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2022</year>). <article-title>Construction of an ecological security pattern in an urban&#x2013;lake symbiosis area: a case study of hefei metropolitan area</article-title>. <source>Remote Sens.</source> <volume>14</volume>, <fpage>2498</fpage>. <pub-id pub-id-type="doi">10.3390/RS14102498</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fan</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Combining landscape patterns and ecosystem services to disclose ecosystem changes in tropical cropland-forest shifting zones: inspiration from Mainland Southeast Asia</article-title>. <source>J. Clean. Prod.</source> <volume>434</volume>, <fpage>140058</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2023.140058</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Feng</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Loo</surname>
<given-names>B. P.</given-names>
</name>
<name>
<surname>Zhong</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>The evolution of urban ecological resilience: an evaluation framework based on vulnerability, sensitivity and self-organization</article-title>. <source>Sust. Cities Soc.</source> <volume>116</volume>, <fpage>105933</fpage>. <pub-id pub-id-type="doi">10.1016/J.SCS.2024.105933</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fletcher</surname>
<given-names>R. J.</given-names>
</name>
<name>
<surname>Didham</surname>
<given-names>R. K.</given-names>
</name>
<name>
<surname>Banks-Leite</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Barlow</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ewers</surname>
<given-names>R. M.</given-names>
</name>
<name>
<surname>Rosindell</surname>
<given-names>J.</given-names>
</name>
<etal/>
</person-group> (<year>2018</year>). <article-title>Is habitat fragmentation good for biodiversity?</article-title> <source>Biol. Conserv.</source> <volume>226</volume>, <fpage>9</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1016/j.biocon.2018.07.022</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jiang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>The spatiotemporal coupling and synergistic evolution of economic resilience and ecological resilience in Africa</article-title>. <source>Sustainability</source> <volume>17</volume>, <fpage>863</fpage>. <pub-id pub-id-type="doi">10.3390/SU17030863</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kong</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Deusen</surname>
<given-names>D.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Evaluating landscape ecological risk through Supply-Demand balance in ecosystem services Evidence from China</article-title>. <source>Ecol. Indic.</source> <volume>173</volume>, <fpage>113355</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2025.113355</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lee</surname>
<given-names>S. I.</given-names>
</name>
</person-group> (<year>2001</year>). <article-title>Developing a bivariate spatial association measure: an integration of Pearson&#x27;s r and Moran&#x27;s I</article-title>. <source>J. Geogr. Syst.</source> <volume>3</volume>, <fpage>369</fpage>&#x2013;<lpage>385</lpage>. <pub-id pub-id-type="doi">10.1007/s101090100064</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Kang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y. J.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Integrating ecosystem services supply-demand balance into landscape ecological risk and its driving forces assessment in Southwest China</article-title>. <source>J. Clean. Prod.</source> <volume>475</volume>, <fpage>143671</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2024.143671</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2025a</year>). <article-title>A Study on the spatiotemporal heterogeneity and driving forces of ecological resilience in the economic Belt on the Northern slope of the tianshan Mountains</article-title>. <source>Land</source> <volume>14</volume>, <fpage>196</fpage>. <pub-id pub-id-type="doi">10.3390/LAND14010196</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>T.</given-names>
</name>
</person-group> (<year>2025b</year>). <article-title>Regional ecological risk assessment and transfer mechanism based on improved gravity and social network analysis model: a case study of Northwest China</article-title>. <source>Ecol. Indic.</source> <volume>172</volume>, <fpage>113243</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2025.113243</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2025c</year>). <article-title>Assessment of landscape ecological risks driven by land use change using multi-scenario simulation: a case Study of Harbin, China</article-title>. <source>Land</source> <volume>14</volume>, <fpage>947</fpage>. <pub-id pub-id-type="doi">10.3390/LAND14050947</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Li</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Adam</surname>
<given-names>N. A.</given-names>
</name>
<name>
<surname>Mirzaliev</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Saydullaev</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Bai</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2025d</year>). <article-title>An integrated approach for evaluating the spatial and temporal dynamics of ecosystem services and ecological risk urban agglomerations in China</article-title>. <source>Ecol. Indic.</source> <volume>172</volume>, <fpage>113248</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2025.113248</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Huang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wen</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Land-use transformation and landscape ecological risk assessment in the Three Gorges Reservoir region based on the &#x201c;production&#x2013;living&#x2013;ecological space&#x201d; Perspective</article-title>. <source>Land</source> <volume>11</volume>, <fpage>1234</fpage>. <pub-id pub-id-type="doi">10.3390/LAND11081234</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Effects of landscape pattern change on ecosystem services and its interactions in karst cities: a case study of Guiyang City in China</article-title>. <source>Ecol. Indic.</source> <volume>145</volume>, <fpage>109646</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2022.109646</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>L&#xfc;</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Fu</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Landscape pattern and ecosystem services are critical for protected areas&#x27; contributions to sustainable development goals at regional scale</article-title>. <source>Sci. Total Environ.</source> <volume>881</volume>, <fpage>163535</fpage>. <pub-id pub-id-type="doi">10.1016/j.scitotenv.2023.163535</pub-id>
<pub-id pub-id-type="pmid">37075999</pub-id>
</citation>
</ref>
<ref id="B18">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ma</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Xie</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Mao</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Spatiotemporal dynamics and simulation of landscape ecological risk and ecological zoning under the construction of free trade pilot zones: a case Study of Hainan Island, China</article-title>. <source>Land</source> <volume>14</volume>, <fpage>940</fpage>. <pub-id pub-id-type="doi">10.3390/land14050940</pub-id>
</citation>
</ref>
<ref id="B19">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Meimei</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Zizhen</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Tengbiao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yongchun</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhuo</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Analysis on absolute conflict and relative conflict of land use in Xining metropolitan area under different scenarios in 2030 by PLUS and PFCI</article-title>. <source>Cities</source> <volume>137</volume>, <fpage>104314</fpage>. <pub-id pub-id-type="doi">10.1016/j.cities.2023.104314</pub-id>
</citation>
</ref>
<ref id="B20">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Peng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Hu</surname>
<given-names>X.</given-names>
</name>
</person-group> (<year>2015</year>). <article-title>Linking ecosystem services and landscape patterns to assess urban ecosystem health: a case study in Shenzhen City, China</article-title>. <source>Landsc. Urban Plan.</source> <volume>143</volume>, <fpage>56</fpage>&#x2013;<lpage>68</lpage>. <pub-id pub-id-type="doi">10.1016/j.landurbplan.2015.06.007</pub-id>
</citation>
</ref>
<ref id="B21">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiao</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Driving mechanisms and threshold identification of landscape ecological risk: a nonlinear perspective from the Qilian Mountains, China</article-title>. <source>Ecol. Indic.</source> <volume>173</volume>, <fpage>113342</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2025.113342</pub-id>
</citation>
</ref>
<ref id="B22">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qiu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Ju</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>The double-edged sword effects of land use optimization based on dual carbon goals: a perspective from landscape ecological risk</article-title>. <source>J. Environ. Manage.</source> <volume>380</volume>, <fpage>125044</fpage>. <pub-id pub-id-type="doi">10.1016/J.JENVMAN.2025.125044</pub-id>
<pub-id pub-id-type="pmid">40154248</pub-id>
</citation>
</ref>
<ref id="B23">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Qu</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhao</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Luo</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>L.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>The effect of the human footprint and climate change on landscape ecological risks: a case study of the Loess Plateau, China</article-title>. <source>Land</source> <volume>11</volume>, <fpage>217</fpage>. <pub-id pub-id-type="doi">10.3390/LAND11020217</pub-id>
</citation>
</ref>
<ref id="B24">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ren</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Dong</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>He</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Lu</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Spatial-temporal changes of landscape ecological risk in the Liuchong river basin from the perspective of production-life-ecological space</article-title>. <source>Front. Environ. Sci.</source> <volume>12</volume>, <fpage>1428058</fpage>. <pub-id pub-id-type="doi">10.3389/FENVS.2024.1428058</pub-id>
</citation>
</ref>
<ref id="B25">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Shamuxi</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Han</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wusimanjiang</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Abudukerimu</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>Q.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Spatial pattern and driving mechanisms of dryland landscape ecological risk: insights from an integrated geographic detector and machine learning model</article-title>. <source>Ecol. Indic.</source> <volume>172</volume>, <fpage>113305</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2025.113305</pub-id>
</citation>
</ref>
<ref id="B26">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Su</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Research on the coupling and coordination of land ecological security and high-quality agricultural development in the Han River Basin</article-title>. <source>Land</source> <volume>13</volume>, <fpage>1666</fpage>. <pub-id pub-id-type="doi">10.3390/LAND13101666</pub-id>
</citation>
</ref>
<ref id="B27">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Zeng</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shen</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>K.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Decoding the spatiotemporal dynamics and driving mechanisms of ecological resilience in the beijing-tianjin-hebei urban agglomeration: a deep learning approach</article-title>. <source>Urban Clim.</source> <volume>61</volume>, <fpage>102436</fpage>. <pub-id pub-id-type="doi">10.1016/J.UCLIM.2025.102436</pub-id>
</citation>
</ref>
<ref id="B28">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Ge</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Ecological resilience of three major urban agglomerations in China from the &#x201c;environment&#x2013;society&#x201d; coupling perspective</article-title>. <source>Ecol. Indic.</source> <volume>169</volume>, <fpage>112944</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2024.112944</pub-id>
</citation>
</ref>
<ref id="B29">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>W.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>A comparative study of urban ecological resilience in the Yangtze River Economic Belt and the Yellow River Basin</article-title>. <source>Hum. Soc. Sci. Commun.</source> <volume>11</volume>, <fpage>1471</fpage>&#x2013;<lpage>15</lpage>. <pub-id pub-id-type="doi">10.1057/S41599-024-04017-X</pub-id>
</citation>
</ref>
<ref id="B30">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Spatiotemporal evolution and obstacle factor analysis of land ecological security in the surrounding areas of Beijing, China</article-title>. <source>Land</source> <volume>14</volume>, <fpage>457</fpage>. <pub-id pub-id-type="doi">10.3390/LAND14030457</pub-id>
</citation>
</ref>
<ref id="B31">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Guo</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>D.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Xi</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>H.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Coupling the calibrated GlobalLand30 data and modified PLUS model for multi-scenario land use simulation and landscape ecological risk assessment</article-title>. <source>Remote Sens.</source> <volume>15</volume>, <fpage>5186</fpage>. <pub-id pub-id-type="doi">10.3390/RS15215186</pub-id>
</citation>
</ref>
<ref id="B32">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Qi</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Sun</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Cao</surname>
<given-names>Y.</given-names>
</name>
<etal/>
</person-group> (<year>2025a</year>). <article-title>Trade-Offs, synergies, and driving factors of ecosystem services in the urban&#x2013;rural fringe of Beijing at multiple scales</article-title>. <source>Land</source> <volume>14</volume>, <fpage>1009</fpage>. <pub-id pub-id-type="doi">10.3390/LAND14051009</pub-id>
</citation>
</ref>
<ref id="B33">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Xue</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Tan</surname>
<given-names>X.</given-names>
</name>
<etal/>
</person-group> (<year>2025b</year>). <article-title>Research on ecosystem service value and landscape ecological risk prediction and zoning: taking Fujian province as an example</article-title>. <source>Ecol. Soc.</source> <volume>507</volume>, <fpage>111173</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLMODEL.2025.111173</pub-id>
</citation>
</ref>
<ref id="B34">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Xu</surname>
<given-names>Q.</given-names>
</name>
<name>
<surname>Yi</surname>
<given-names>J.</given-names>
</name>
</person-group> (<year>2025c</year>). <article-title>Multi-scale identification of driving forces for landscape ecological risk: a case study of the basin containing three plateau lakes in Yunnan, China</article-title>. <source>Ecol. Indic.</source> <volume>170</volume>, <fpage>113037</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2024.113037</pub-id>
</citation>
</ref>
<ref id="B35">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2025d</year>). <article-title>Landscape ecological risk based on optimal scale and its tradeoff/synergy with human activities: a case study of the Nanjing metropolitan area, China</article-title>. <source>Ecol. Indic.</source> <volume>170</volume>, <fpage>113040</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2024.113040</pub-id>
</citation>
</ref>
<ref id="B36">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Song</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Zhu</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Jiang</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Z.</given-names>
</name>
</person-group> (<year>2025e</year>). <article-title>Coupling study on ecological resilience and urban-rural integrated development in China</article-title>. <source>Front. Environ. Sci.</source> <volume>13</volume>, <fpage>1552533</fpage>. <pub-id pub-id-type="doi">10.3389/FENVS.2025.1552533</pub-id>
</citation>
</ref>
<ref id="B37">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Wu</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Yin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Revealing the nonlinear impact of environmental regulation on ecological resilience using the XGBoost-SHAP model: evidence from the Yangtze River Delta region, China</article-title>. <source>J. Clean. Prod.</source> <volume>514</volume>, <fpage>145700</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2025.145700</pub-id>
</citation>
</ref>
<ref id="B38">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu F</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X. J.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>Ecological management zoning through integration of ecosystem service and landscape ecological risk: a case Study in chongli, China</article-title>. <source>Land</source> <volume>14</volume>, <fpage>1133</fpage>. <pub-id pub-id-type="doi">10.3390/LAND14061133</pub-id>
</citation>
</ref>
<ref id="B39">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Xu W</surname>
<given-names>W.</given-names>
</name>
<name>
<surname>Jin</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Guan</surname>
<given-names>T.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Coupling coordination degree, interaction relationship and driving mechanism of water resources carrying capacity of beijing-tianjin-hebei urban agglomeration in China</article-title>. <source>J. Clean. Prod.</source> <volume>504</volume>, <fpage>145433</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2025.145433</pub-id>
</citation>
</ref>
<ref id="B40">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Liu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>G.</given-names>
</name>
</person-group> (<year>2018</year>). <article-title>Quantifying spatio-temporal patterns of urban expansion in Beijing during 1985&#x2013;2013 with rural-urban development transformation</article-title>. <source>Land Use Pol.</source> <volume>74</volume>, <fpage>220</fpage>&#x2013;<lpage>230</lpage>. <pub-id pub-id-type="doi">10.1016/j.landusepol.2017.07.004</pub-id>
</citation>
</ref>
<ref id="B41">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Yang</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Landscape ecological risk assessment and driving factors analysis based on optimal spatial scales in Luan River Basin, China</article-title>. <source>Ecol. Indic.</source> <volume>169</volume>, <fpage>112821</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2024.112821</pub-id>
</citation>
</ref>
<ref id="B42">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yang</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Yan</surname>
<given-names>J. X.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>G.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>W. W.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Gao</surname>
<given-names>P. C.</given-names>
</name>
<etal/>
</person-group> (<year>2025</year>). <article-title>Dynamic simulation of landscape ecological risk in Gansu&#x2019;s farming-pastoral ecotone based on topographic gradients</article-title>. <source>Front. Environ. Sci.</source> <volume>13</volume>, <fpage>1523892</fpage>. <pub-id pub-id-type="doi">10.3389/FENVS.2025.1523892</pub-id>
</citation>
</ref>
<ref id="B43">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yin</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhou</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>N.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Impact of regional integration policy on urban ecological resilience: a case study of the Yangtze river Delta region, China</article-title>. <source>J. Clean. Prod.</source> <volume>485</volume>, <fpage>144375</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2024.144375</pub-id>
</citation>
</ref>
<ref id="B44">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Yu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Microagroecosystem models and their landscape ecological risk evolution in karst mountainous areas</article-title>. <source>Ecol. Indic.</source> <volume>158</volume>, <fpage>111502</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2023.111502</pub-id>
</citation>
</ref>
<ref id="B45">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zeng</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Wu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Chen</surname>
<given-names>W. J.</given-names>
</name>
</person-group> (<year>2024</year>). <article-title>Coupling analysis of land use change with landscape ecological risk in China: a multi-scenario simulation perspective</article-title>. <source>J. Clean. Prod.</source> <volume>435</volume>, <fpage>140518</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2023.140518</pub-id>
</citation>
</ref>
<ref id="B46">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>Y.</given-names>
</name>
</person-group> (<year>2025</year>). <article-title>How does urbanization evolve heterogeneously in urbanized, urbanizing, and rural areas of China? Insights from ecosystem service value</article-title>. <source>Geogr. Sustain.</source> <volume>6</volume>, <fpage>100254</fpage>. <pub-id pub-id-type="doi">10.1016/J.GEOSUS.2024.100254</pub-id>
</citation>
</ref>
<ref id="B47">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Lv</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2023a</year>). <article-title>Spatio-temporal assessment of landscape ecological risk and associated drivers: a case study of the Yellow River Basin in Inner Mongolia</article-title>. <source>Land</source> <volume>12</volume>, <fpage>1114</fpage>. <pub-id pub-id-type="doi">10.3390/LAND12061114</pub-id>
</citation>
</ref>
<ref id="B48">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Ma</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Gong</surname>
<given-names>J.-W.</given-names>
</name>
<name>
<surname>Chu</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>L. J.</given-names>
</name>
</person-group> (<year>2023b</year>). <article-title>A coupling effect of landscape patterns on the spatial and temporal distribution of water ecosystem services: a case study in the Jianghuai ecological economic zone, China</article-title>. <source>Ecol. Indic.</source> <volume>151</volume>, <fpage>110299</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2023.110299</pub-id>
</citation>
</ref>
<ref id="B49">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Yuan</surname>
<given-names>R.</given-names>
</name>
<name>
<surname>Jarvie</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>Q.</given-names>
</name>
</person-group> (<year>2023c</year>). <article-title>Landscape ecological risk of China&#x27;s nature reserves declined over the past 30 years</article-title>. <source>Ecol. Indic.</source> <volume>156</volume>, <fpage>111155</fpage>. <pub-id pub-id-type="doi">10.1016/J.ECOLIND.2023.111155</pub-id>
</citation>
</ref>
<ref id="B50">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Yu</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Feng</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Yao</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Li</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Du</surname>
<given-names>H.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Landscape ecological risk and drivers of land-use transition under the perspective of differences in topographic gradient</article-title>. <source>Land</source> <volume>13</volume>, <fpage>876</fpage>. <pub-id pub-id-type="doi">10.3390/LAND13060876</pub-id>
</citation>
</ref>
<ref id="B51">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Zhu</surname>
<given-names>Y.</given-names>
</name>
<name>
<surname>Shi</surname>
<given-names>Z.</given-names>
</name>
<name>
<surname>Wang</surname>
<given-names>X.</given-names>
</name>
<name>
<surname>Zou</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Zhang</surname>
<given-names>N.</given-names>
</name>
<name>
<surname>Zheng</surname>
<given-names>P.</given-names>
</name>
<etal/>
</person-group> (<year>2024</year>). <article-title>Spatial-temporal evolution of ecological risks in a typical urban lake in Southwestern China</article-title>. <source>J. Clean. Prod.</source> <volume>434</volume>, <fpage>140273</fpage>. <pub-id pub-id-type="doi">10.1016/J.JCLEPRO.2023.140273</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>