基本信息中文题目2000年以来道路扩张导致西伯利亚原生森林景观破碎化加剧四倍英文题目Fourfold increase in fragmentation from road expansion in Siberian intact forests since 2000期刊名称Nature Ecology Evolution发表时间2026年7月27日接收2026年9月1日在线发表作者信息张紫荆Zi Jing Zhang共同一作、李仲玢Zhongbin B. Li共同一作兼通讯、王梦秋Mengqiu Wang共同通讯、Alexey Egorov、David P. Edwards、陈吉泉Jiquan Chen、伍棋坪Qiping Wu、冯耀伟Yaowei Feng、罗莉Li Luo、杨澜Lan Yang、樊彬Bin Fan、马宇乐Yule Ma、康葳Wei Kang、张汉奎Hankui K. Zhang、许振赐Zhenci Xu作者单位武汉大学遥感信息工程学院、香港大学地球与行星科学系/地理系、俄罗斯生物多样性保护中心、美国爱达荷大学、英国剑桥大学植物科学系、美国密歇根州立大学、加利福尼亚大学河滨分校、南达科他州立大学等DOIhttps://doi.org/10.1038/s41559-026-03162-7核心观点 道路网络的快速扩张是全球森林退化、生境破碎化和碳损失的关键诱因但其对高纬度寒带原生森林景观Intact Forest Landscapes, IFLs的长期生态扰动此前缺乏定量解析。本研究基于 2000–2023 年间 20,480 景 Landsat 遥感影像重建了覆盖俄罗斯中东部西伯利亚 581 万 km² 区域的长时序道路网络演变数据库。研究发现该区域道路面积与长度均增长近 1.6 倍道路扩张直接占用了约 1,046 km² 的 IFLs并导致原生林斑块数量激增 329.5%破碎化程度加剧 4.3 倍斑块平均面积从 4,460 km² 骤降至 1,038 km²30 m 道路边缘区扩大 1.2 倍。道路扩张主要受木材采伐与采矿活动驱动并在 2010–2023 年间直接引致 39.1 兆吨Mt的森林地上生物量AGB损失。研究表明即便在划定的原生林与自然保护区内道路侵蚀依然显著亟需出台具有刚性法律约束力的空间规划以平衡经济开发与寒带森林生态安全。️ 一、 引言部分简要总结 (Introduction Summary)现实背景与科学背景 原生森林景观作为地球上保存最完整、受人类工业活动干扰极小的连片森林生态系统在维持全球生物多样性、调控水文过程以及提供巨量碳汇方面发挥着不可替代的作用。西伯利亚森林覆盖着约 171 万 km² 的原生林是地球三大温寒带原生林景观之一。然而伴随全球及区域资源消费需求的快速攀升西伯利亚地区的大规模工业伐木、矿产开采以及化石能源开发急剧升温导致俄罗斯亚洲部分森林碳汇在 1990–2020 年间衰减了 42%尤其在 2010 年代末期呈现出加速萎缩态势。前人研究局限与核心科学瓶颈研究重心过度偏向热带寒带森林长期处于认知盲区过去数十年关于道路驱动毁林、生境破碎化及边缘效应的研究绝大多数集中于亚马孙盆地、刚果盆地及东南亚等热带雨林针对北方针叶林Taiga野生动植物栖息地完整性受道路网络侵蚀的定量证据严重匮乏。现有全球道路数据集的严重漏检与数据缺陷应用广泛的全球道路数据库如 OpenStreetMap、GRIP在高纬度寒冷地区存在更新滞后、漏检严重以及缺乏多年代际连续观测等问题此外现有数据集多为无宽度的单线矢量格式无法精确衡量道路施工走廊Road corridors及其路堤、路肩对地表的真实破坏面积。驱动机制与碳亏缺的综合归因缺失以往针对西伯利亚森林受损的研究多聚焦于野火和森林虫害在多年代际尺度上道路扩张如何直接穿透原生林、如何诱发斑块几何构型退化以及道路建设引致的地上生物量损失规模缺乏基于严格像元级遥感反演与地统计学的系统性量化。本文试图解决的核心科学问题 2000–2023 年间西伯利亚中东部道路网络的空间扩展速率与地理分异格局如何新建道路对原生森林景观的直接侵占、空间切割破碎化以及边缘效应产生了何种程度的影响采矿与伐木活动在何种程度上塑造了道路网络的空间格局道路扩张对寒带森林地上生物量造成了多少直接碳亏损 二、 研究方法详细介绍与分析 (Methodology Framework)技术路线与时空范围界定 研究区横跨俄罗斯中东部西伯利亚东经 75°02′–147°21′北纬 42°17′–64°17′涵盖萨哈共和国RS、克拉斯诺亚尔斯克边疆区KK、哈巴罗夫斯克边疆区KhK、伊尔库茨克州IO等行政单元总面积达 581 万 km²。研究选择 2000 年、2010 年和 2023 年三个基准时间节点开展全覆盖对比分析。多源遥感数据源与无云合成Landsat 影像预处理通过 Google Earth Engine 平台调用 20,480 景经几何与辐射校正的 Landsat 5/7/8/9 表面反射率产品Collection 2 Tier 1其中 2000 年 3,914 景、2010 年 4,699 景、2023 年 11,867 景。利用 QA_PIXEL 波段严格掩膜云与云阴影全部选取无积雪覆盖的植被生长季数据采用中位数合成算法Median compositing生成 30 m 分辨率年度无云真彩色合成影像。环境与地理辅助变量引入 30 m 分辨率全球数字表面模型 AW3D30 提取高程与坡度获取 GLC_FCS30 土地覆盖产品并重分类为森林、草甸灌丛、农田、湿地和其他 5 类整合世界原生森林景观数据库IFL基准年 2000 年下载欧空局气候变化倡议ESA Biomass_cci v5.01100 m 空间分辨率森林地上生物量AGB数据集2010–2021 年收集世界保护区数据库WDPA与全球陆地生态区TEOW矢量图层。道路网络人工解译与设计基精度检验道路走廊Road corridors数字化7 名遥感专业科研人员耗时约 5,040 小时在 30 m 合成影像上独立开展全境人工目视解译与矢量化包括路面及推土扰动的路堤、路缘区。在 ArcGIS Pro 中采用测地线算法Geodesic measurement精确量算道路面积与长度并将道路走廊平均宽度量化为面积与长度之比。设计基精度评估Design-based Approach依据 Olofsson et al. 准则构建分层随机抽样框架定义道路权重 W_road 0.003与非道路W_non-road 0.997两个图层抽取 402 个独立验证样点结合 Google Earth 历史超高分辨率影像~0.5 m进行盲审评定。误差矩阵计算表明道路图层总体精度达到 0.9991 ± 0.0003制图精度为 0.974–1.000漏分误差极低用户精度为 0.71 ± 0.09受 30 m 混合像元略微高估路宽影响。工业扰动提取与数学建模工作流采矿与伐木区提取采矿多边形基于 Maus et al. 全球矿业数据集并结合 2023 年 Landsat 影像人工增补伐木斑块基于 PyTorch 搭建深度学习语义分割模型 SegFormer利用 16,422 个 256 × 256 像元切片样本实施数据增强与训练最终模型平均交并比mIoU达 0.91F1 分数达 0.95。景观破碎化与边缘效应量化以 2000 年 IFL 数据单块面积 ≥ 500 km² 且距道路 1 km 以上为基线将各期道路叠加于原生林斑块提取斑块数量、平均面积、总周长并划分为 1、1–10、10–100、100–200 及 200 km² 五个规模等级。向斑块内部建立 30 m、60 m 和 100 m 缓冲区计算边缘受损面积。卡方分布与蒙特卡洛置换检验构建 Pearson 卡方统计量公式 4结合 5,000 次蒙特卡洛随机置换检验道路扩张在土地覆盖类型、高程带100 m, 100–600 m, 600 m及坡度带0–5°, 5–15°, 15–30°, 30°上的选择性偏好。地理加权回归GWR采用自适应双平方核函数Adaptive bisquare kernel与修正赤池信息准则AICc寻优带宽解析道路扩张与采矿、伐木区增量之间的空间非平稳关联。残差自协变量回归RAC与生物量误差传播以 10 km × 10 km 网格n 18,423统计归一化道路扩张与归一化 AGB 损失量构建 k 近邻空间权重矩阵k 50通过引入残差空间自协变量RAC剔除空间自相关偏误。生物量不确定性核算采用地统计学球状变异函数模型在公式6中完全计入像元误差的空间自相关传播效应。 三、 主要研究结果与充实数据分析 (Main Results Quantified Data)1. 2000–2023 年道路网络的大规模扩张特征在 2000 至 2023 年的 23 年间西伯利亚中东部的道路网络规模扩张了近 1.6 倍。详见 Figure 1a, b。规模激增道路占地总面积从 2000 年的 7,110 km² 增加到 2023 年的 11,056 km²净增长 55.5%道路总长度由 172,521 km 延伸至 275,288 km净增长 59.6%。在此期间由 Landsat 观测到的道路走廊平均宽度稳定保持在约 48 m。区域分异格局道路绝对扩张集中在四大重点行政区——萨哈共和国RS面积增加 663 km²增长 122.8%长度增加 11,558 km增长 86.5%、克拉斯诺亚尔斯克边疆区KK面积增加 639 km²增长 62.6%长度增加 20,104 km增长 80.9%、哈巴罗夫斯克边疆区KhK面积增加 574 km²增长 109.6%长度增加 15,800 km增长 126.2%以及伊尔库茨克州IO面积增加 493 km²增长 32.0%长度增加 22,303 km增长 56.0%。详见 Figure 1c, d。下垫面与地形选择性侵占82.7% 的新建道路直接开凿于原始森林内部侵占灌丛和草甸的占 10.3%占用耕地的占 4.3%。卡方检验表明道路空间配置具有极显著的地形选择性偏好chi^2 266.9, P 0.001。多达 70.35% 的道路扩张密集分布在海拔 100–600 m 的低地平原Figure 1e且绝大部分坐落于坡度小于 5° 的平缓地带。空间向北与向东推移在 IO 和 KK道路网由传统的北纬 58°N 显著向 63°N 的高纬度深处推进在滨海边疆区PK道路则向东经 140°E 沿海地带延伸直接服务于油气勘探开采Figure 1f与皆伐作业Figure 1g。图片名Figure 1图例说明Figure 1 | 2000–2023 年西伯利亚中东部卫星遥感观测的道路网络扩张特征。(a) 2000 年上层灰色与 2023 年底层蓝色道路网络空间分布图(b) 2000、2010 和 2023 年道路面积与长度的阶段变化统计(c, d) 各行政单元 2000–2023 年净增道路面积 (c) 与净增道路长度 (d)(e) 道路扩张面积在 1 至 2,500 m 高程梯度上的分布比例标明 100–600 m 主体高程带占 70.35%(f) 萨哈共和国境内油气开采道路在 2000 年勘探期与 2023 年开采期的 Landsat 遥感影像对比(g) 伊尔库茨克州皆伐林区道路网络扩张的遥感影像对比。2. 采矿与伐木工业活动对道路扩张的空间驱动卫星遥感图层空间叠加与 GWR 模型证实道路网络的延伸与高强度的自然资源掠夺存在直接的空间共轭关系详见 Figure 2a, b。采矿活动的驱动全区采矿用地的扩张与道路扩张呈现强烈的正向空间耦合GWR 局部回归平均系数为 0.2平均判定系数 R^2 为 0.3在赤塔州ZK及哈巴罗夫斯克南部尤为显著。见 Figure 2a。木材采伐的驱动工业采伐对道路的牵引效应更为突出GWR 平均回归系数达 0.3平均 R^2 为 0.4伊尔库茨克州与克拉斯诺亚尔斯克南部呈现密集的采伐运输道路网。见 Figure 2b。图片名Figure 2图例说明Figure 2 | 西伯利亚中东部卫星观测道路与采矿、伐木工业活动的空间叠加图。(a) 2000 年灰色与 2023 年蓝色道路网络叠加 2023 年卫星提取的采矿分布区红色(b) 道路网络叠加 2023 年通过 SegFormer 深度学习模型识别的木材采伐区绿色。3. 原生森林景观IFL的直接侵占与破碎化加剧四倍道路的纵横切削使西伯利亚原生森林景观的完整性遭受了结构性重创详见 Figure 3。直接与间接破坏面积2000–2023 年间新建道路直接占压清除了 1,046 km² 的 IFL 核心区域Figure 3a。在道路两侧建立空间缓冲区分析显示道路 1 km 间接影响带内损失的 IFL 面积高达 30,427 km²占全区该时期 IFL 损失总面积31.4 万 km²的 9.7%。其中萨哈共和国承受了最大的直接侵占294 km²而克拉斯诺亚尔斯克边疆区在 1 km 缓冲带内的间接破坏面积居首达 10,025 km²。如 Figure 3b 所示。斑块裂变与四倍破碎化原生森林斑块数量从 2000 年的 383 个攀升至 2010 年的 628 个增长 64.0%到 2023 年更是爆发式增长至 1,645 个较 2000 年基准增加了 329.5%意味着生境破碎化程度加剧了 4.3 倍。见 Figure 3c。斑块几何形态萎缩伴随斑块数量的倍增IFL 斑块的平均面积从 2000 年的 4,460 km² 急剧塌陷至 2010 年的 2,719 km²至 2023 年进一步萎缩至仅 1,038 km²平均斑块规模缩减了 76.7%。见 Figure 3d。微型孤立斑块激增斑块规模分级统计显示小于 1 km² 的极微型破碎斑块在 2000 年基本不存在而 2023 年激增至数百个昔日广袤的连续荒野被分割为高度孤立的生境孤岛。如 Figure 3e 所示。边缘效应急剧扩张IFL 斑块总周长从 2000 年的 257,872 km 持续拉长至 2023 年的 300,107 km。沿道路向林内的 30 m 边缘区面积从 2000 年的 7,686 km²占 IFL 总面积的 0.45%增至 2023 年的 8,934 km²占比升至 0.64%累计扩张 16.0%100 m 边缘区面积更是由 2000 年的约 25,600 km² 激增至 2023 年的 31,000 km² 以上占 IFL 的比例由 1.50% 攀升至 2.13%。见 Figure 3f。图片名Figure 3图例说明Figure 3 | 2000–2023 年道路扩张导致西伯利亚中东部原生森林景观IFL毁林、破碎化与边缘效应加剧特征。(a) 2023 年道路网络穿透切断 IFL 核心斑块的典型遥感局部案例绿色为 2020 年现存 IFL灰色为 2000–2020 年损失区蓝色为道路(b) 道路两侧不同距离30 m 至 1 km缓冲区内的 IFL 间接损失面积内插小图展示 30–90 m 近距细节(c) 2000、2010 和 2023 年 IFL 斑块总周长随斑块数量的演变(d) 斑块平均面积随斑块数量激增的骤降趋势(e) 不同面积等级1、1–10、10–100、100–200、200 km²森林斑块数量的时序变化(f) 道路 30 m、60 m 和 100 m 缓冲区内森林边缘受损面积左轴及其在 IFL 中的面积占比右轴。4. 道路开发导致的森林地上生物量亏缺利用 ESA Biomass_cci 遥感反演数据论文系统量化了 2010–2023 年间新建道路所直接引起的植被碳损失详见 Figure 4地上生物量AGB损失总量2010 至 2023 年间道路占地直接清除的活立木地上生物量累计达到 39.1 兆吨Mt以 2010 年生物量为基准附带空间自相关校正后的 3 倍标准差误差棒。行政区损失贡献克拉斯诺亚尔斯克边疆区KK的 AGB 损失最高达 8.5 Mt哈巴罗夫斯克边疆区KhK次之损失 6.1 Mt伊尔库茨克州IO损失 5.8 Mt萨哈共和国RS损失 5.6 Mt托木斯克州TO损失 4.6 Mt。上述五个行政区合计贡献了全区道路开发引致 AGB 损失总量的 78.3%。见 Figure 4b。空间自相关校正回归在 10 km × 10 km 网格尺度上经残差自协变量RAC模型过滤空间自相关干扰后归一化道路扩张量与归一化 AGB 损失量之间呈现高度显著的正向线性相关回归斜率 beta 0.46模型判定系数 R^2 0.56样本量 n 18,423P 0.001。如 Figure 4c 所示。边缘带生物量损失强度梯度在道路两侧伴随微气候退化与林木枯死单位面积 AGB 损失强度呈现出强烈的距离衰减规律紧邻道路的 30 m 边缘带内 AGB 损失强度接近 2,000 tons/km²至 90 m 处降至约 1,770 tons/km²而在 1 km 处下降至约 1,250 tons/km²。见 Figure 4d。图片名Figure 4图例说明Figure 4 | 2010–2023 年西伯利亚中东部道路网络直接侵占引致的森林地上生物量AGB损失估算。(a) 2010–2023 年新增道路叠加欧空局 2010 年森林 AGB 空间底图Mg ha⁻¹(b) 各行政区道路直接占用的 AGB 损失总量柱状图附带考虑像元空间自相关误差传播的 ±3 倍标准差误差棒(c) 基于残差自协变量模型RAC拟合的归一化道路扩张与归一化 AGB 损失的散点回归beta 0.46, R^2 0.56, n 18,423(d) 道路两侧不同缓冲距离内单位面积 AGB 损失强度的沿程衰减曲线tons/km²。 四、 讨论部分的不足与未来展望 (Discussion Limitations Perspectives)研究不足与不确定性来源30 m Landsat 空间分辨率的物理局限受限于 Landsat 像元尺寸细小林间便道容易发生混合像元效应可能导致高反射率简易路面宽度的系统性高估与此同时茂密针叶林树冠重叠遮蔽下的狭窄拖拉机道Canopy-obscured skid trails无法被中分辨率光学传感器捕捉可能导致局部林区实际道路长度存在漏估。光谱“绿化”与真实生态恢复的认知偏差时序追踪显示2000–2010 年建设的道路中有约 85% 在 2023 年影像中依然清晰可见道路扰动具有极强持久性即使少数废弃道路表面出现先锋草本入侵并在卫星图上呈现光谱闭合但寒带冻土区森林垂直结构、群落分层与碳储量的彻底恢复往往需要数十年乃至上百年地表植被复绿绝不代表生态系统完整性的即时恢复。生物量评估的时间窗口截断受欧空局高质量 Biomass_cci 数据起始年份限制本研究仅系统评估了 2010–2023 年间的 AGB 亏损2000–2010 年前期的碳损失未能纳入定量核算。道路与人为野火的级联复合效应道路修筑打破了高寒针叶林的天然地理阻隔大幅提升了人类可达性进而显著放大了因人为弃火、吸烟引发森林火灾的风险本研究侧重于道路本身的直接与物理间接破坏未将道路诱发火灾的衍生生态损害完全解耦分离。未来研究展望与保护治理策略融合高精度星载激光雷达与雷达遥感立体监测后续研究需引入米级高分辨率光学数据如 Sentinel-2、具备穿透树冠能力的合成孔径雷达Sentinel-1, ALOS/PALSAR并协同利用 GEDI 与 ICESat-2 空间激光测高激光雷达精准解构路网两侧微地形沉陷、冻土退化以及树冠三维高度与生物量垂直断层的空间微动态。消除保护区“纸面公园”现象设立刚性禁路红线空间叠加分析发现2000 年以来已有 260 km² 的新增道路公然入侵了法定自然保护区内部直接导致 40.9 km² 原生森林景观被蚕食。这表明仅有保护区或 IFL 名义 designation 无法阻挡资源开发必须由国家立法强制设立不可逾越的“无路区”Roadless areas红线严禁采矿与皆伐作业向核心区延伸。跨国基建与战略通道的前瞻性生态空间规划在俄罗斯“远东大开发”倡议、西伯利亚大铁路支线复线升级及东西伯利亚-太平洋ESPO油气输运管道布局中必须建立先导性环境影响评价EIA机制在交通走廊选线中推行“避让第一”原则优先整合已退化路网空间严厉限制新线路切割东西伯利亚针叶林核心生态区防止世界关键温寒带天然碳库退化为不可逆的碳源排放地。