《Sustainability》：Spatiotemporal Dynamics of Landscape Ecological Risk Under Vegetation Loss and Urban Expansion in Dhaka
编辑推荐：
摘要：景观生态风险(Landscape Ecological Risk, LER)反映景观变化对生态系统结构、功能及稳定性的潜在不利影响。在达卡等快速城市化特大城市中，植被丧失与建设用地扩张近几十年来加剧了环境压力。本研究考察了2004—2024年达卡在植被变
摘要：景观生态风险(Landscape Ecological Risk, LER)反映景观变化对生态系统结构、功能及稳定性的潜在不利影响。在达卡等快速城市化特大城市中，植被丧失与建设用地扩张近几十年来加剧了环境压力。本研究考察了2004—2024年达卡在植被变化与城市扩张共同影响下LER的时空动态。研究人员利用多时相遥感数据制作土地覆被图，基于归一化差值植被指数(NDVI)估算 Fractional Vegetation Cover（FVC，植被覆盖度），并使用夜间灯光(Nighttime Light, NTL)数据量化城市化强度。景观生态风险指数(Landscape Ecological Risk Index, LERI)由景观格局指标计算获得，采用双变量空间自相关和地理加权回归(Geographically Weighted Regression, GWR)检验空间关联性与局部空间异质性。结果表明，2004—2024年间研究区34.39%的区域发生植被退化，高风险区占比从2004年的24.36%上升至2024年的42.95%。土地覆被分析进一步显示建设用地大幅扩展，植被、农用地及低地类呈萎缩与破碎化趋势。空间分析揭示植被覆盖、城市化强度与生态风险的关系在城市内部存在差异且随时间推移空间分异性增强。研究发现植被丧失与城市扩张在空间上与达卡生态风险升高相关联。受遥感数据不确定性、无监督分类、重采样及地面验证有限等因素影响，结论应谨慎解读。尽管存在局限，本研究提供了理解生态风险动态的空间显式框架，并为快速城市化地区的绿地保护、生态修复与可持续城市规划提供了依据。
论文解读：《Sustainability》刊载——达卡地区植被丧失与城市扩张背景下景观生态风险(LER)的时空动态研究
一、研究背景与立项依据
快速城市化与无计划土地利用转换加剧了全球特大城市的环境挑战，南亚孟加拉国首都达卡(Dhaka)人口近二十年迅猛增长，建成区持续向外侵蚀自然与农用地，导致植被丧失、景观破碎化及生态脆弱性上升。景观生态风险(Landscape Ecological Risk, LER)是评估景观变化对生态系统结构、功能及稳定性潜在不利影响的有效框架，Fractional Vegetation Cover（FVC，植被覆盖度）可表征植被状况，Nighttime Light（NTL，夜间灯光）数据可作为城市化强度的代理指标，Landscape Ecological Risk Index（LERI，景观生态风险指数）则基于景观格局指标量化区域生态风险。已有达卡相关研究多单独关注土地利用变化或生态风险，常采用传统统计方法而未能捕捉空间异质性与动态交互作用；植被恢复(Vegetation Restoration, VR)与城市扩张对LER时空演化的联合效应尚缺乏深入探讨，也少有将生态风险评估转化为空间显式规划建议的研究。为此，研究人员开展本次综合研究，旨在厘清2004—2024年达卡FVC、NTL、土地覆被与LERI的时空变化，揭示植被条件、城市化强度与LER的空间关系及空间异质性影响，并为生态保护与城市增长管理提供分区规划建议。该论文发表于国际期刊《Sustainability》。
二、主要关键技术方法概述
研究人员选取达卡市建成区约306.38 km2为研究区，获取2004年Landsat 5 TM、2014年Landsat 8 OLI及2024年Landsat 9 OLI-2 Collection 2 Level-2地表反射率影像，在Google Earth Engine(GEE)中对每年3—7月影像进行大气校正并生成中值合成影像，采用Weka K-means无监督聚类生成五类土地覆被（建设用地、农用地、低地、植被、水体）；基于NDVI线性像元分解模型估算30 m分辨率FVC；选用VIIRS月度NTL数据经双线性插值重采样至30 m网格并作归一化处理。以30×30像元渔网网格为评价单元，计算每个单元内各地类面积占比，结合景观破碎化指数(Fragmentation Index, Fi)、分离度指数(Si)及优势度指数(Di)（权重分别取0.5、0.3、0.2）构建景观干扰指数(Landscape Disturbance Index, LDi)，依据专家打分法赋各土地覆被类型景观脆弱指数(Landscape Vulnerability Index, LFi)并归一化，进而计算LERI并由Kriging插值表达空间分布，按Natural Breaks将LERI划为五级生态风险等级。采用FRAGSTATS 4.2计算各土地覆被类型的Patch Density（PD，斑块密度）、Largest Patch Index（LPI，最大斑块指数）与Splitting Index（SPLIT，分裂指数）以量化景观结构变化。对FVC与LERI做双变量全局及局部空间自相关分析(Local Indicators of Spatial Association, LISA)，并通过Geographically Weighted Regression（GWR，地理加权回归）以行政区街道面汇总栅格变量，自适应带宽、高斯核函数探索FVC、NTL与LERI间局部空间变异性关系。
三、研究结果
3.1 植被恢复(Vegetation Restoration, VR)与退化动态
FVC分级统计显示2004—2024年极低与低FVC合计占比由46.83%升至54.61%，高FVC占比下降7.83个百分点；植被变化转移分析表明34.39%区域发生植被退化，优于29.23%区域有植被恢复(VR)，净退化面积15.69 km2。研究人员由此得出结论：达卡整体植被条件趋于下降，植被退化范围超过恢复范围。
3.2 景观生态风险(LER)等级分布
2004年高风险与极高风险区合计占44.59%，2014年升至66.38%，2024年为59.70%；其中高风险类由2004年24.36%增至2024年42.95%，极低与低风险合计由31.66%降至15.98%。研究人员指出高风险区显著扩张、低风险区收缩，说明生态风险总体升高并在2004—2014年变化最为剧烈。
3.3 LERI时序变化(LERI Change over Time)
2004—2014年44%区域LERI升高，32%降低；2014—2024年63.48%区域LERI降低，30.99%升高；全时段(2004—2024)升、降风险区分别占44.45%与35.42%。研究人员认为早期生态风险普遍加剧，近期部分区域风险有所回落，反映城市成熟或景观发育过程中的差异响应。
3.4 夜间灯光指数(Nighttime Light Index, NTL)
NTL空间格局由2004年连片高强度照明转为2014、2024年以极低和低光强为主、高强度呈局部簇团状分散分布，均值在后期有所回升。研究人员解读此变化体现达卡近年城市增长以既有城区填充式开发与集约化为主，而非均质蔓延。
3.5 土地覆被(Land Cover, LC)变化
建设用地由中心城区及交通廊道向周边大幅扩展，农用地、植被及低地面积递减且破碎化明显，水体范围缩小且斑块更零散。研究人员总结土地覆被变化特征为建设用地持续侵占生态用地并引致自然景观组分萎缩与分割。
3.6 LERI与FVC的双变量局部空间自相关分析(Local Autocorrelation Analysis Between LERI and FVC)
High–Low（高生态风险–低植被覆盖）聚类占比由2004年33%增至2024年42%，High–High（高生态风险–高植被覆盖）由9%升至13%，Low–High（低生态风险–高植被覆盖）由28%降至18%。研究人员认为高生态风险区空间集聚性增强，FVC与LER间的简单负向空间关联随城市化推进趋于减弱和非均匀化。
3.7 地理加权回归(GWR)模型结果
FVC–LERI的GWR显著负相关区（Vegetation Negative Zones, VNZ）由2004年44.37%降至2024年9.49%，显著正相关区（Vegetation Positive Zones, VPZ）由2.88%升至33.98%；NTL–LERI显著相关区由2004年50.69%升至2024年66.53%，Urban Positive Zones（UPZ，NTL与LERI正相关区）由6.70%升至16.61%。研究人员提出FVC与LERI的反向局部关系随景观破碎化被削弱，NTL与LERI的正向局部关联范围扩大，印证城市化强度与生态风险的空间非平稳关系。
四、讨论与结论总结（翻译结论部分核心内容）
本研究以2004—2024年多时相遥感数据、景观格局指标及空间统计分析集成框架，考察了达卡在植被丧失与快速城市扩张共同作用下景观生态风险(LER)的时空动态。研究发现：研究期内达卡发生显著植被退化及持续性建设用地扩张，虽局部存在植被恢复(VR)但退化范围超过恢复，整体植被状况下滑；同时生态高风险区明显扩增、低风险区缩减，反映城市系统生态压力增大及景观破碎化加剧；空间分析表明植被覆盖、城市化强度(NTL)与生态风险的关系于城市内部并不均一——高强度城市扩张与高NTL区通常伴生更高生态风险，而植被覆盖与LER的空间关系随时间推移分化愈加明显，说明达卡生态风险受植被变化、土地覆被转型及空间不均衡城市增长的共同塑造。规划层面提示需保护并连通城市绿地、限制无计划建成区蔓延，并将生态风险考量纳入城市发展政策，空间显式分析有助于识别局部风险格局并支撑在地化规划干预。受土地覆被分类未经混淆矩阵定量验证及GWR结果属探索性空间关联而非因果推断等限制，结论宜在方法学范畴内审慎解读。未来研究结合野外生态观测、更高空间分辨率数据集及因果建模方法，可深化对城市生态动态及城市化–植被变化–生态风险三者关系的认识。
---研究分区规划建议简述（原文Section 5浓缩）：依LERI、FVC、NTL及LC空间分异，研究人员将达卡划分为四类空间导控区——西部/中部/南部划为城市增长管理区(Urban Growth Management Zones)，严控现存植被湿地转不透水面，提倡紧凑竖向开发并配建绿色基础设施；西、西南及临河带划为生态修复优先区(Ecological Restoration Priority Zones)，实施城市造林、湿地低地修复及河岸缓冲带建设；东部及东北保留较低生态风险及较稳定植被农地区划为生态保育与连通区(Ecological Conservation and Connectivity Zones)，防止无序开发以维持生境连续性与生态系统服务；北部城乡过渡带（Uttarkhan、Dakshinkhan等）为管理过渡区(Managed Transition Zones)，需在进一步退化前以生态敏感分区管控并预留生态廊道，最终强调通过连接东部生态源地与残存湿地、河道及破碎绿地来强化城市生态网络。
