华东传统矿业城市生态网络重构研究

传统矿业城市生境碎片化严重,重构生态网络对提升区域生态系统服务和维护生物多样性及景观多样性具有重要意义。以徐州市中心城区为对象,采用形态学空间格局分析(MSPA)和景观指数方法,提取对维持城市生态系统健康具有重要价值的生境作为源地,利用最小累积阻力模型(MCR)构建阻力面并提取潜在生态廊道,并结合电路理论、图论和网络结构分析等方法定量评价和优化生态网络,重构徐州市中心城区生态安全网络。结果表明,徐州市中心城区生态源地共237个,面积为6 421.49 hm~2,主要分布于西南部;尚需添加70块小生境和192条新的生态廊道,才能弥补现有生态网络的不足;重构后的生态网络使徐州市中心城区景观破碎化程度和分离度大幅度降低,斑块连通性和斑块结合度提高,区域景观格局和生态安全显著改善。该研究可为传统矿业城市转型、低效用地再开发以及国土空间规划提供科学依据。     

城乡建设用地景观格局多尺度差异

以郑州市域内1989、2000和2009年的典型城市与农村样地（25 km×25 km）为例,应用空间自相关、半方差、以及九项轨迹方差分析方法,对城市与农村建设用地景观多尺度上的空间格局特征与动态差异进行研究,结果表明,（1）城市建设用地景观聚集强度较高,斑块规模较大,而农村建设用地景观在大部分尺度上呈随机或均匀分布,斑块规模较小。（2）随着尺度的增加,城市及农村建设用地景观分布状态都由聚集向均匀过渡,但前者过渡较为平缓,而后者在小尺度范围（0~2km）变化速度较快,在较大尺度（3~16 km）上其分布状态的波动较小。（3）1989-2009年,城市与农村建设用地景观动态特征有较大差异,并表现出尺度依赖性。城市较小尺度上（0~4 km）建设用地景观的分布状态表现出聚集-分散的趋势,在较大尺度上（8~14 km）表现出分散-聚集的变化趋势。农村小尺度上（0~4 km）建设用地景观的分布状态也表现出聚集-分散的趋势,但大尺度（4~16 km）上状态变化很小且没有规律。     

快速城市化进程中城市扩张对景观格局分异特征的影响

在中国快速城市化背景下,大部分研究仅是单一的利用景观指数分析城市景观格局,往往只侧重于对景观斑块或格局的几何特征进行简单分析和描述,忽略了城市扩张与城乡景观格局动态变化过程之间的内部联系,研究城市扩张与城乡景观格局动态变化过程之间的内部联系对揭示城市化规律有着重要意义。文章以江苏省昆山市（县）范围为研究对象,基于全市1985—2008年间多期TM影像图,在采用遥感与地理信息系统相结合的技术方法获取昆山市景观格局动态变化以及相关景观生态学指数分析的基础上,探讨城市扩张对城乡景观格局的影响。结果表明：23年间自然景观类型的面积基本保持不变;人工景观面积显著增涨了4倍。斑块密度、景观斑块数量破碎度指数受人工景观类型影响基本呈上升趋势;斑块密度、景观破碎度、最大斑块、边缘密度等指数与建设用地增长之间存在明显相关性,表明城市扩张是城乡景观格局特征变化的重要驱动力。     

哈尔滨市绿地景观格局与过程的连通性和完整性

城市绿地景观的空间结构和生态过程的研究是城市景观生态系统研究的重要内容和基本特色之一．应用景观生态学原理，在遥感和地理信息系统(GIS)技术支持下对哈尔滨城市现有景观生态格局进行了研究和分析，认为主要存在以下景观生态问题：景观生态连通性低；景观结构单一；绿地面积少，且绿地空间分布不均衡．在现有景观格局基础上，哈尔滨市可望通过以下几个方面改善景观生态过程和格局的连续性：在景观生态战略点开辟新的绿色斑块；建立和完善生态化的水系廊道网络；建立城效绿化网络，连接绿地系统的点、线、面；增加城市绿地系统中的景观异质性和多样性；在城市扩展中维护景观生态过程与格局的连续性．图4表4参19     

基于Landsat时间序列数据的重庆市热力景观格局演变分析

城市热岛效应是伴随城市化而产生的生态环境问题,研究城市热力景观格局的演变有助于掌握城市热岛效应变化的机制与规律,为城市产业合理布局、城市生态环境改善提供科学的决策支持。以重庆市主城九区为研究对象,基于2001年Landsat7 ETM+、2007年Landsat5 TM以及2014年Landsat8 OLI_TIRS三期遥感影像数据,利用辐射传输方程法和Jiménez-Mu?oz et al.(2014)的分裂窗算法反演地表温度,并在此基础上,计算热力景观格局指数,研究重庆市热力景观格局的演变过程。结果表明,(1)2001—2014年,重庆市主城九区热岛和强热岛景观类型范围不断扩大,热岛效应明显增强,渝中区热岛和强热岛所占百分比最高,江北区、南岸区、九龙坡区和大渡口区的增长速度较快。(2)热力景观类型中热岛面积的增加主要由正常区斑块转化而来,而强热岛面积的增加主要由正常区和热岛斑块转化而来,跨越正常区的斑块转化较难。(3)在斑块类型水平上,热岛与强热岛斑块优势度增大。2001—2014年热岛斑块密度减小3.33,平均斑块面积增大6倍;强热岛斑块密度减少0.65,平均斑块面积增大5倍。热岛和强热岛斑块变得大而集中,破碎度减小。热岛间的连通性与强热岛斑块间的连通性越来越高,连通性指数分别增大了3.51和8.41,强热岛斑块形状的复杂程度逐期变大。在景观水平上,2001—2014年重庆市的斑块数量和斑块密度减小,平均斑块面积增大,热力景观破碎化程度逐渐降低,斑块连通性指数高,均大于99.5。聚合度和均匀度指数分别增大14.85和0.09,像素间聚合成斑块的程度变大且斑块类型面积越来越均匀。由此可见,随着城市的发展,重庆市的热环境问题越来越严重,利用热力景观格局指数分析城市热环境,可了解城市热力景观格局的演变趋势,为热岛的缓解提供理论依据。     

天津滨海新区湿地环境演化与景观格局动态

从历史演化和景观格局动态变化两个方面对天津滨海新区湿地生态系统的演变进行了研究。历史上天津地区为海洋,自滨海新区成陆以来,经历了多次的海浸和洪水淹没,海浸期间发生陆地湿地化,海退和洪水退却时发生水体湿地化。本文总结了自地质时期以来的历史演化模型。同时在GIS和RS技术支持下,运用FRAG-STATS软件,研究了1979-2008年滨海新区湿地景观格局动态变化特征。结果表明,随着滨海新区的发展和土地利用强度的增大,研究区景观斑块数量增加、平均斑块面积减少的趋势比较明显,景观呈现破碎化的趋势;总体上,近30年中,滨海新区湿地总面积变化不大,并没有大面积持续减少。景观破碎程度从1979年开始持续上升,中间略有波动,但总体上升较快。     

珠江三角洲地区高速公路沿线景观格局变化研究

根据1995年和2005年土地利用数据,借助地理信息系统软件和景观生态学方法研究珠江三角洲区域发展轴线——广深高速公路和广珠高速公路运营前后沿线区域的土地利用和景观格局变化,并进一步对比分析相同功能等级道路的生态效应差异。结果表明：（1）高速公路运营后,沿线区域建设用地增加迅速,耕地大量流失且迅速向建设用地转化;（2）公路的建设运营使得沿线区域人工干扰强度加大,景观斑块密度增加、面积减少、斑块形状趋向复杂,景观破碎化加剧;（3）道路途径区域的社会经济条件以及原有的景观格局特征是导致道路生态破碎效应差异的重要原因;（4）公路对自然景观的影响效应远大于其对人工景观的影响效应。     

西南天然林保护工程区生态成效评估

随着大尺度生态保护和恢复工程的开展和不断深化,工程有效性评估逐渐成为备受关注的热点问题。开展大尺度生态保护和恢复的成效评估对于相应的宏观决策与管理具有重要科学意义。基于年均植被覆盖度、景观格局指数和人类威胁因子评估指标体系,对西南4省市天然林保护工程区的生态保护成效进行定量评估。结果表明:研究区生态保护成效总体上明显改善,地类景观斑块分布更加合理规整,连通性提高,人类干扰强度逐步减小,尤其是4省市交界林区以及重庆与湖北、四川与陕西、四川与三江源地区的交界区域生态成效显著提高;而汶川周边林区、横断山脉部分林区以及云南文山林区生态保护成效较弱。因此,在天然林保护工程的进一步实施过程中有必要在分析保护成效空间分异特征的基础上采取保护和恢复并重的策略,将植被、景观格局和生境质量与人类活动的干扰和压力综合起来,以促进大尺度生态保护和恢复的可持续有效管理。     

山西中部城市群景观格局演变对其热环境的影响研究

城市群化使得城市热环境问题变得更为复杂。以2010-2020年4期Landsat影像数据为基础,运用空间分析、景观格局指数与数理统计等方法,分析山西中部城市群景观格局动态演变和热环境分布及变化特征,并定量探究城市群景观格局演变对其热环境的影响作用。研究表明,（1）近10年来山西省中部城市群的建设用地逐年增加,在空间分布上均有向太原市偏移的趋势。（2）研究区地表温度总体呈上升趋势,城市群热岛效应逐渐显现,2013-2016年,高温区重心向东北方移动30.03 km,转移至太原市中轴线偏北方位置,各城市间热环境的连通性加强。（3）城市群的景观格局演变对其热环境有着明显的影响：(1)各景观类型中建设用地、耕地及草地热环境贡献指数较高,林地与水体贡献指数为负,其中林地贡献值为-0.975,贡献指数绝对值最大;(2)从景观格局指数来看,在类型水平上,林地的聚集度指数、平均分布斑块面积、最大斑块所占景观面积比例、景观类型比例与地表温度（LST）呈负相关,斑块密度、边缘密度与LST呈正相关,而在景观水平上,景观形状指数、斑块密度、香浓均匀度指数等与LST呈正相关,2013年各相关系数最大,不同景观类...     

中型景观尺度下斑块类型和地理特征影响杨树人工林干部病害的发生——以河南省清丰县为例

景观病理学在开展病原扩散、病害发生及其严重程度方面以其全新的视角,为森林病害的区域控制提供了新的研究技术及理论支持。首次利用景观病理学原理和方法对河南省清丰县一个中型景观下杨树人工林干部病害发生特征开展了研究,目的是解析在大尺度下斑块类型和地理特征对杨树人工林病害发生的影响。在100 km2的调查区域,以种植方式和林分类型划分斑块类型,分析显示发病株率在不同斑块间差异显著：农田间作斑块的林木发病株率显著低于孤立斑块、纯林斑块、混交林等斑块的发病株率;但发病株率在孤立斑块、纯林斑块及混交林等斑块间无显著差异。抚育管理措施对预防和减轻杨树人工林干部病害的发生起到关键左右：精细管理林分（有修枝、施肥和锄草）的林木发病株率（p=0.001）和发病指数（p〈0.001）均显著低于粗放管理林分（无修枝、无施肥和锄草等）。人类活动,如无序修剪和放牧很可能是造成村落附近林分发病率显著高于其他地点林分的主要原因。采用logistic回归,以品种编号、树龄、树高、林分密度、林分郁闭度、林分类型、斑块类型、地理特征,等为自变量建立病害发生预测模型。方程拟合达到极显著水平（Wald=71.248,p〈0.001）。方程总的预测正确率为68.2%,发病的预测正确率为79.8%。     

Assessing effects of land use on landscape connectivity: loss and fragmentation of western U.S. forests

Effects of land-use change on the conservation of biodiversity have become a concern to conservation scientists and land managers, who have identified loss and fragmentation of natural areas as a high-priority issue. Despite urgent calls to inform national, regional, and state planning efforts, there remains a critical need to develop practical approaches to identify where important lands are for landscape connectivity (i.e., linkages), where land use constrains connectivity, and which linkages are most important to maintain network-wide connectivity extents. Our overall goal in this paper was to develop an approach that provides comprehensive, quantitative estimates of the effects of land-use change on landscape connectivity and illustrate its use on a broad, regional expanse of the western United States. We quantified loss of habitat and landscape connectivity for western forested systems due to land uses associated with residential development, roads, and highway traffic. We examined how these land-use changes likely increase the resistance to movement of forest species in non-forested land cover types and, therefore, reduce the connectivity among forested habitat patches. To do so, we applied a graph-theoretic approach that incorporates ecological aspects within a geographic representation of a network. We found that roughly one-quarter of the forested lands in the western United States were integral to a network of forested patches, though the lands outside of patches remain critical for habitat and overall connectivity. Using remotely sensed land cover data (ca. 2000), we found 1.7 million km2 of forested lands. We estimate that land uses associated with residential development, roads, and highway traffic have caused roughly a 4.5% loss in area (20 000 km2) of these forested patches, and continued expansion of residential land will likely reduce forested patches by another 1.2% by 2030. We also identify linkages among forest patches that are critical for landscape connectivity. Our approach can be readily modified to examine connectivity for other habitats/ecological systems and for other geographic areas, as well as to address more specific requirements for particular conservation planning applications.     

Coupling urban 3-D information and circuit theory to advance the development of urban ecological networks

Ongoing, rapid urban growth accompanied by habitat fragmentation and loss challenges biodiversity conservation and leads to decreases in ecosystem services. Application of the concept of ecological networks in the preservation and restoration of connections among isolated patches of natural areas is a powerful conservation strategy. However, previous approaches often failed to objectively consider the impacts of complex 3-D city environments on ecological niches. We used airborne lidar-derived information on the 3-D structure of the built environment and vegetation and detailed land use and cover data to characterize habitat quality, niche diversity, and human disturbance and to predict habitat connectivity among 38 identified habitat core areas (HCAs) in Nanjing, China. We used circuit theory and Linkage Mapper to create a landscape resistance layer, simulate habitat connectivity, and identify and prioritize important corridors. We mapped 64 links by using current flow centrality to evaluate each HCA's contribution and the links that facilitate intact connectivity. Values were highest for HCA links located in the west, south, and northeast of the study area, where natural forests with complex 3-D structures predominate. Two smaller HCA areas had high centrality scores relative to their extents, which means they could act as important stepping stones in connectivity planning. The mapped pinch-point regions had narrow and fragile links among the HCAs, suggesting they require special protection. The barriers with the highest impact scores were mainly located at the HCA connections to Purple Mountain and, based on these high scores, are more likely to indicate important locations that can be restored to improve potential connections. Our novel framework allowed us to sufficiently convey spatially explicit information to identify targets for habitat restoration and potential pathways for species movement and dispersal. Such information is critical for assessing existing or potential habitats and corridors and developing strategic plans to balance habitat conservation and other land uses based on scientifically informed connectivity planning and implementation.     

Connecting science,policy, and implementation for landscape‐scale habitat connectivity

We examined the links between the science and policy of habitat corridors to better understand how corridors can be implemented effectively. As a case study, we focused on a suite of landscape‐scale connectivity plans in tropical and subtropical Asia (Malaysia, Singapore, and Bhutan). The process of corridor designation may be more efficient if the scientific determination of optimal corridor locations and arrangement is synchronized in time with political buy‐in and establishment of policies to create corridors. Land tenure and the intactness of existing habitat in the region are also important to consider because optimal connectivity strategies may be very different if there are few, versus many, political jurisdictions (including commercial and traditional land tenures) and intact versus degraded habitat between patches. Novel financing mechanisms for corridors include bed taxes, payments for ecosystem services, and strategic forest certifications. Gaps in knowledge of effective corridor design include an understanding of how corridors, particularly those managed by local communities, can be protected from degradation and unsustainable hunting. There is a critical need for quantitative, data‐driven models that can be used to prioritize potential corridors or multicorridor networks based on their relative contributions to long‐term metacommunity persistence.     

Invasive species impacts on ecosystem structure and function: A comparison of the Bay of Quinte, Canada, and Oneida Lake, USA, before and after zebra mussel invasion

As invasion rates of exotic species increase, an ecosystem level understanding of their impacts is imperative for predicting future spread and consequences. We have previously shown that network analyses are powerful tools for understanding the effects of exotic species perturbation on ecosystems. We now use the network analysis approach to compare how the same perturbation affects another ecosystem of similar trophic status. We compared food web characteristics of the Bay of Quinte, Lake Ontario (Canada), to previous research on Oneida Lake, New York (USA) before and after zebra mussel (Dreissena polymorpha) invasion. We used ecological network analysis (ENA) to rigorously quantify ecosystem function through an analysis of direct and indirect food web transfers. We used a social network analysis method, cohesion analysis (CA), to assess ecosystem structure by organizing food web members into subgroups of strongly interacting predators and prey. Together, ENA and CA allowed us to understand how food web structure and function respond simultaneously to perturbation. In general, zebra mussel effects on the Bay of Quinte, when compared to Oneida Lake, were similar in direction, but greater in magnitude. Both systems underwent functional changes involving focused flow through a small number of taxa and increased use of benthic sources of production; additionally, both systems structurally changed with subgroup membership changing considerably (33% in Oneida Lake) or being disrupted entirely (in the Bay of Quinte). However, the response of total ecosystem activity (as measured by carbon flow) differed between both systems, with increasing activity in the Bay of Quinte, and decreasing activity in Oneida Lake. Thus, these analyses revealed parallel effects of zebra mussel invasion in ecosystems of similar trophic status, yet they also suggested that important differences may exist. As exotic species continue to disrupt the structure and function of our native ecosystems, food web network analyses will be useful for understanding their far-reaching effects.     

Ranking individual habitat patches as connectivity providers: Integrating network analysis and patch removal experiments

Here we propose an integrated framework for modeling connectivity that can help ecologists, conservation planners and managers to identify patches that, more than others, contribute to uphold species dispersal and other ecological flows in a landscape context. We elaborate, extend and partly integrate recent network-based approaches for modeling and supporting the management of fragmented landscapes. In doing so, experimental patch removal techniques and network analytical approaches are merged into one integrated modeling framework for assessing the role of individual patches as connectivity providers. In particular, we focus the analyses on the habitat availability metrics PC and IIC and on the network metric Betweenness Centrality. The combination and extension of these metrics jointly assess both the immediate connectivity impacts of the loss of a particular patch and the resulting increased vulnerability of the network to subsequent disruptions. In using the framework to analyze the connectivity of two real landscapes in Madagascar and Catalonia (NE Spain), we suggest a procedure that can be used to rank individual habitat patches and show that the combined metrics reveal relevant and non-redundant information valuable to assert and quantify distinctive connectivity aspects of any given patch in the landscape. Hence, we argue that the proposed framework could facilitate more ecologically informed decision-making in managing fragmented landscapes. Finally, we discuss and highlight some of the advantages, limitations and key differences between the considered metrics.     

Using network centrality measures to manage landscape connectivity

We use a graph-theoretical landscape modeling approach to investigate how to identify central patches in the landscape as well as how these central patches influence (1) organism movement within the local neighborhood and (2) the dispersal of organisms beyond the local neighborhood. Organism movements were theoretically estimated based on the spatial configuration of the habitat patches in the studied landscape. We find that centrality depends on the way the graph-theoretical model of habitat patches is constructed, although even the simplest network representation, not taking strength and directionality of potential organisms flows into account, still provides a coarse-grained assessment of the most important patches according to their contribution to landscape connectivity. Moreover, we identify (at least) two general classes of centrality. One accounts for the local flow of organisms in the neighborhood of a patch, and the other accounts for the ability to maintain connectivity beyond the scale of the local neighborhood. Finally, we study how habitat patches with high scores on different network centrality measures are distributed in a fragmented agricultural landscape in Madagascar. Results show that patches with high degree and betweenness centrality are widely spread, while patches with high subgraph and closeness centrality are clumped together in dense clusters. This finding may enable multispecies analyses of single-species network models.     

Multiple structural modifications to dendritic ecological networks produce simple responses

Structural modifications to landscapes affect the ability of organisms to access different habitat patches. There exist, however, very few general methods by which to relate modifications to expectations of effects, and even fewer that enable understanding of how multiple modifications may interact. In the absence of any guiding principles, ecologists have assumed that interactions will result in complex landscape-scale effects. One way of understanding such effects is through rendering a landscape as a graph or network, among the simplest of which are dendritic networks typified by stream systems. Yet even for stream networks, there are no known general principles concerning the nature of interactions between multiple modifications. We developed a model to describe the ability of fish to access and use different habitat patches within dendritic networks. We used mathematical and numerical analyses of the model to investigate how the habitat value of a network is affected by changes in connectivity and habitat quality, and then to examine interactions between multiple modifications. Rather than showing complex interactions, our analytic and simulation-based results show that the combined effect of multiple modifications approximately equals the sum of individually predicted effects. Dendritic networks thus appear to respond far more simply to multiple modifications than has previously been assumed. These results have implications for stream management planning, and offer a firm foundation from which to better understand population processes within dendritic networks.     

Multi-temporal land use and cover changing analysis: the environmental impact in Mediterranean area

ABSTRACT

Human activity shapes the levels of anthropogenic pressure that depend on the land management method adopted. This has a fundamental role in the transformation of traditional landscapes. This study focuses on a representative region of the Mediterranean area with the objective to analyse the landscape’s dynamics, to detect the spatial arrangement of class patches, to identify the main agroecosystem characters and to provide a framework to assess ecosystems services. In order to assess land use/land cover changes and landscape persistence, the period between 1960 and 2012 was analysed, taking into consideration the years 1960, 2000 and 2012 using comparable land use maps. Land use and land cover analysis show an urban area growth of 24% during 2000–2012 and of 523% over between 1960 and 2012. The very high levels of land abandonment up to the year 2000 (+7216%) have reversed their trend between 2000 and 2012 (?95%). The orchards showed a relevant increase, particularly after 2000, while the vineyards were linked to the highest value of surface erosion (?74%). The outcomes showed that urban settlements can damage the ecological network with negative effects on the landscape’s environmental sustainability in proximity of significant urban centres. Instead, the ecological network is well preserved and highly associated to the agricultural areas when there is the persistence of many land uses and low urban density, despite the presence of dynamic changes.     

Validating graph-based connectivity models with independent presence–absence and genetic data sets

Habitat connectivity is a key objective of current conservation policies and is commonly modeled by landscape graphs (i.e., sets of habitat patches [nodes] connected by potential dispersal paths [links]). These graphs are often built based on expert opinion or species distribution models (SDMs) and therefore lack empirical validation from data more closely reflecting functional connectivity. Accordingly, we tested whether landscape graphs reflect how habitat connectivity influences gene flow, which is one of the main ecoevolutionary processes. To that purpose, we modeled the habitat network of a forest bird (plumbeous warbler [Setophaga plumbea]) on Guadeloupe with graphs based on expert opinion, Jacobs’ specialization indices, and an SDM. We used genetic data (712 birds from 27 populations) to compute local genetic indices and pairwise genetic distances. Finally, we assessed the relationships between genetic distances or indices and cost distances or connectivity metrics with maximum-likelihood population-effects distance models and Spearman correlations between metrics. Overall, the landscape graphs reliably reflected the influence of connectivity on population genetic structure; validation R² was up to 0.30 and correlation coefficients were up to 0.71. Yet, the relationship among graph ecological relevance, data requirements, and construction and analysis methods was not straightforward because the graph based on the most complex construction method (species distribution modeling) sometimes had less ecological relevance than the others. Cross-validation methods and sensitivity analyzes allowed us to make the advantages and limitations of each construction method spatially explicit. We confirmed the relevance of landscape graphs for conservation modeling but recommend a case-specific consideration of the cost-effectiveness of their construction methods. We hope the replication of independent validation approaches across species and landscapes will strengthen the ecological relevance of connectivity models.     

Ecological network and emergy analysis of urban metabolic systems: Model development,and a case study of four Chinese cities

Analysis of the structure and function of urban metabolic systems is an important goal of urban research. We used network pathways and network utility analysis to analyze the basic network structure of the urban metabolic system and the complex ecological relationships within the system, providing a new way to perform such research. Using four Chinese cities as examples, we developed an ecological network model of the urban metabolic system. By using network pathway analysis, we studied the changing relationships between metabolic length and the number of metabolic pathways, and between metabolic length and reachability. Based on the distribution of the number of metabolic pathways, we describe the basic structure and intercompartment relationships of the system. By using the sign distribution in the network utility matrix, we determined the ecological relationships and degree of mutualism between the compartments of the system. The basic components of the system consisted of the internal environment, the external environment, and the agricultural, industrial, and domestic sectors. With increasing metabolic length, the ecological relationships among the components of the system became more diverse, and the numbers of metabolic paths and their reachability improved. Although the basic network structure of the four cities was identical, the mutualism index differed. Beijing's mutualism index was superior to that of Shanghai, and much higher than those of Tianjin and Chongqing. By analyzing the structure and function of the urban metabolic system, we provide suggestions for optimizing the structure and adjusting the relationships, and propose methods for the application of ecological network analysis in future urban system research.