《Environmental and Sustainability Indicators》：Spatial insights on edge-driven degradation in protected tropical forests of Bangladesh
编辑推荐：
摘要：保护区内（PAs）的热带森林退化对生物多样性、碳储存和生态系统服务构成了一种隐蔽但持续的威胁。然而，退化在精细尺度上的格局，特别是沿保护区边界附近的退化动态，目前仍知之甚少。本研究利用空间显式方法评估了2001–2023年间孟加拉国三个保护区——巴瓦尔国
摘要：保护区内（PAs）的热带森林退化对生物多样性、碳储存和生态系统服务构成了一种隐蔽但持续的威胁。然而，退化在精细尺度上的格局，特别是沿保护区边界附近的退化动态，目前仍知之甚少。本研究利用空间显式方法评估了2001–2023年间孟加拉国三个保护区——巴瓦尔国家公园（BNP）、特克纳夫野生动物保护区（TWS）和雷马-卡伦加野生动物保护区（RKWS）——边缘驱动的森林损失。研究人员使用Hansen全球森林变化数据集（30 m分辨率），量化了同心缓冲带（0–500 m、500–1000 m、1000–1500 m）内的树冠覆盖损失，并按冠层密度进行分层，以识别结构性脆弱的森林斑块。研究发现，退化集中在边缘地带：RKWS和TWS在0–500 m区域内的损失最高（分别为12.4%和11.7%）。在所有研究地点，靠近道路加剧了干扰强度，但空间梯度因地而异。BNP在50 m缓冲区内损失最大，而RKWS和TWS在500–1000 m较宽缓冲区内出现主要损失峰值。依赖冠层的损失格局也因地点而异：在RKWS中，高冠层森林（覆盖率>75%）受影响最为严重，而BNP的损失集中在65-75%冠层等级，TWS则集中在>80%冠层等级。BNP表现为扩散型退化，而TWS和RKWS则呈现加剧的边缘破碎化。研究结果强调，仅靠法律保护不足以维持该国的森林完整性。将高分辨率遥感与空间分析相结合，可提供一种可扩展的方法来检测退化热点并为恢复优先事项提供指导。这种方法可为可持续森林管理提供信息，并支持实现全球恢复承诺的进展，包括波恩挑战和联合国生态系统恢复十年。
**论文解读：孟加拉国受保护热带森林边缘驱动退化的空间格局与启示**
**一、研究背景与问题**
热带森林是全球生物多样性最丰富、生态功能最重要的陆地生态系统之一，对气候调节、碳固存、水源涵养及当地民生具有不可替代的作用。然而，即使在法律划定的保护区内，热带森林仍面临选择性采伐、农业扩张、道路建设、非法采伐和火灾等多重人为干扰。森林退化往往以隐蔽形式发生，表现为冠层变薄、生物量选择性移除、破碎化和结构简化，虽未造成完全毁林，却严重损害生态系统功能。传统监测手段如土地利用/覆盖分类图和NDVI等植被指数，通常只能捕捉粗尺度上的毁林事件，难以检测冠层以下的渐进性退化过程。更关键的是，退化在空间上并非均匀分布，而是集中于人类可达性高、治理薄弱的森林边缘和道路廊道沿线。这种边缘驱动的退化格局在国家级或生物群系尺度的评估中往往被掩盖。
孟加拉国森林覆盖率仅为国土面积的17.08%，其中约10.31%为法律保护地。位于人口密集或快速城市化地区的巴瓦尔国家公园（简称BNP）、特克纳夫野生动物保护区（简称TWS）和雷马-卡伦加野生动物保护区（简称RKWS）面临非法土地利用、基础设施扩张和资源偷猎等压力。TWS还因罗兴亚难民危机而承受额外的生态压力。尽管已有全国尺度的毁林研究，但保护区内退化的空间分布——特别是与道路、边界和森林结构的关系——仍缺乏精细尺度的认识。本研究正是基于这一空白，利用高分辨率遥感数据与空间分析方法，揭示保护区内部边缘驱动的退化机制。
**二、研究设计与技术方法**
本研究采用Hansen全球森林变化数据集（GFC，30 m分辨率，版本1.11，2023年），基于Google Earth Engine提取2001–2023年森林损失信息。研究人员以2000年树冠覆盖为基线，沿保护区边界向内建立0-500 m、500-1000 m和1000-1500 m三个同心缓冲带，量化归一化树冠损失率。通过OpenStreetMap提取道路网络，生成50 m、200 m、500 m和1000 m道路缓冲区，分析道路邻近性对退化强度的影响。按1%间隔划分冠层覆盖等级，计算每个等级内的比例损失，识别易受干扰的结构类型。基于Landsat影像、使用随机森林分类器生成2000年、2010年和2023年森林/非森林二值图，采用连通分量标记算法计算斑块数量，指示破碎化程度。
**三、主要研究结果**
*3.1 边缘效应与道路邻近性*
所有三个保护区的森林损失均呈现强烈的年际波动，且干扰事件不成比例地集中在道路附近。BNP于2020年出现显著峰值，50 m缓冲区损失达约22-24%，随距离递减（500 m处约20%，1000 m处约18-19%），表现为典型的单调道路边缘梯度。RKWS于2022年出现最大峰值，但梯度呈倒置特征：1000 m缓冲区损失最高（约28-31%），50 m缓冲区最低（约13-16%），表明干扰从道路边缘向内部扩散并形成景观尺度的脉冲事件。TWS于2017年出现主要峰值，1000 m缓冲区损失约33-35%，同样高于近路区域。结果表明道路邻近性是驱动退化的持续因素，但不同地点的空间梯度模式存在明显差异。
*3.2 边缘效应与保护区边界距离*
三个保护区在0-500 m边界缓冲区内均呈现最高归一化树冠损失率：BNP为11.497%，TWS为11.713%，RKWS为12.4%。随着距离增加，损失率逐步下降：在500-1000 m区域，BNP降至8.695%，TWS降至9.822%，RKWS降至8.884%；在1000-1500 m区域，进一步降至6.370%、8.536%和7.149%。热力图分析显示边界附近存在局部损失热点，尤其在BNP和RKWS中。
*3.3 冠层覆盖依赖性损失*
BNP中低冠层（0-20%）像素占比高但相对损失低，中期冠层（10-40%）出现中等程度损失峰值，高冠层（>75%）损失率虽高（如66%等级达28.4%）但像素数量极少，代表性有限。TWS中损失集中在高冠层区域：85%等级损失57.8%（3451像素），97%等级损失57.4%（141像素），表明致密森林遭遇广泛退化。RKWS表现出最显著的模式：85-87%冠层等级的损失极高（分别为83.7%、86.5%和71.5%），呈现初始冠层密度越高、损失越大的清晰梯度。
*3.4 斑块数量与森林损失的时间动态*
2000年至2023年间，斑块数量在全部地点均大幅增加。TWS增幅最剧烈，从191个增至1022个（+435.1%）；RKWS从93个增至314个（+237.6%）；BNP从537个增至861个（+60.9%）但2010年后趋于稳定。三个地点合计从821个增至2197个（+167.6%）。森林损失面积方面：RKWS累计损失最高，于2004年（103.95 ha）、2007年（80.82 ha）和2022年（83.79 ha）出现明显峰值；TWS自2001年（85.5 ha）整体下降但于2017年反弹（47.17 ha）；BNP长期保持较低水平，仅在2010年（15.57 ha）和2020年（23.04 ha）有轻微增加。
**四、讨论与结论总结**
讨论部分指出，道路邻近性与保护区边界效应共同塑造了空间异质的退化格局。BNP的经典边缘梯度与其城郊区位、土地侵占和工业扩张相关；TWS的2017年损失峰值与罗兴亚难民涌入造成的燃料木材需求激增吻合；RKWS的高冠层损失集中于高价值林区，与土地权属争议和茶园扩张密切相关。研究人员强调，法律保护不足以确保森林完整性，高冠层森林的结构完好性并不等于生态韧性。恢复策略必须因地点而异：RKWS应优先保护残余高冠层森林，TWS应重点关注连接度恢复，BNP应维持现有森林完整性并防止新的破碎化。全球恢复目标如波恩挑战和联合国生态系统恢复十年，需要超越单纯面积指标，采纳空间诊断性与治理敏感的评估框架。
研究结论部分的核心表述为：本研究对孟加拉国三个保护区提供了空间精细尺度的森林退化评估，证明法律指定本身不足以防止森林损失。尽管拥有保护地位，TWS、RKWS和BNP在2001-2023年间仍遭受了由人口迁移、小农侵占和工业扩张等不同地点特定压力驱动的实质性森林损失。结构完整性不是韧性的可靠指标，因为高冠层区域遭受了最高程度的退化，凸显了森林面积本身不足以作为保护成功的度量标准。退化格局在空间和时间上呈现异质性，表明恢复和管理干预必须针对每个地点的特定生态和社会政治背景进行定制，而不是依赖标准化方法。持久的恢复成果需要支持执法、促进社区有意义参与并确保长期持续投入的治理体系。孟加拉国的案例表明，有效的森林恢复超越了植树指标或正式保护；它需要整合细尺度空间诊断、适应性管理和参与式治理，以解决森林退化的根本驱动因素。只有协调生态、社会和制度维度，恢复工作才能产生有意义且持久的影响，将退化景观转变为具有韧性的功能性森林。然而，该研究未纳入一手社会经济数据，未来研究可结合遥感数据与家庭调查以更严格地评估差异森林损失的因果归因。
