《Geographies》：Sediment-Driven Expansion of Tropical Mangroves in the Bengawan Solo Delta Revealed by Multi-Decadal Google Earth Engine Analysis
编辑推荐：
## 摘要
红树林砍伐是一个全球性问题，但以沉积物为主的河口可以通过其自身的形态动力学过程抵御这一问题。本研究重建了Ujung Pangkah河口（1995-2025年）的时空演变轨迹，以权衡自然向海推进与人为压力之间的作用。研究人员在Google Eart
## 摘要
红树林砍伐是一个全球性问题，但以沉积物为主的河口可以通过其自身的形态动力学过程抵御这一问题。本研究重建了Ujung Pangkah河口（1995-2025年）的时空演变轨迹，以权衡自然向海推进与人为压力之间的作用。研究人员在Google Earth Engine（GEE）平台中利用Landsat影像档案和随机森林（RF）分类器，结合四种光谱指数（归一化差异植被指数（NDVI）、修正归一化差异水体指数（mNDWI）、增强型植被指数（EVI）和红树林植被指数（MVI））对河口进行制图，并对全部四个时期进行了独立验证（总体精度92.25-95.00%；Kappa系数0.845-0.900）。误差校正的变化检测分析显示了一条非单调的演变轨迹：红树林面积从1995年的898.88公顷增长至2015年的2410.44公顷峰值，随后在2025年缩减至1816.87公顷（误差校正后：从1261.11公顷增至2626.03公顷，再缩减至2213.71公顷）。在整个30年记录中，这一变化具有统计学意义的净扩张（z = 3.60，p < 0.001）。空间归因分析表明，2015-2025年的缩减主要来自陆向人为转化（78.1%），而非海向侵蚀（21.9%）。悬浮沉积物代理指标与十年尺度净变化之间的滞后相关性（r = 0.92）为持续的沉积驱动型海岸向海推进提供了定量支持。因此，Ujung Pangkah的韧性与其真实的、可定位的人为压力并存。保护这一蓝碳生态系统需要将政策针对已在进行中的内部转化区域，同时持续保护海岸前沿。
# 沉积驱动的热带红树林扩张与人为收缩：Ujung Pangkah河口三十年时空动态研究解读
## 一、研究背景与科学问题
红树林是全球生物地球化学活性最高、生态功能最重要的生态系统之一。红树林位于陆海动态交错带，其复杂根系结构能够衰减水动力能量、缓解海岸侵蚀，并为海岸带提供抵御风暴潮和海平面上升的天然屏障。同时，三角洲和河口红树林被公认为重要的蓝碳汇——其在淹水河口沉积物中的长期碳封存能力甚至超过热带陆地森林，在全球气候变化减缓议程中具有重要地位。
然而，过去三十年间，以海岸城市化和潮间带大规模改造为养殖塘为主的人为压力，导致全球红树林面积大幅减少。尽管印度尼西亚拥有全球最大面积的红树林，但同时其红树林丧失速率也居高不下，主要驱动因素为虾塘养殖扩张。值得注意的是，某些以沉积物为主的河口水系展现出独特生态轨迹。Ujung Pangkah河口位于印度尼西亚东爪哇Bengawan Solo河入海口，该区域的红树林前沿非但在人为压力下退缩，反而在持续河流沉积物通量的促进下经历向海推进。Bengawan Solo三角洲的形态演变历来受大量沉积物排放塑造。
历史研究中，理解复杂多年代形态动力学变化一直受方法论局限制约。以往研究多依赖时间受限的实地调查或局部卫星影像分析，难以捕捉长期地貌转变。此外，全球遥感文献主要聚焦于红树林砍伐和丧失的量化，对热带三角洲系统中自然沉积驱动的多年代向海推进机制研究仍存在关键空白。基于云端的对地观测计算平台——Google Earth Engine（GEE）结合随机森林（RF）等机器学习算法，为突破上述局限提供了技术路径。
本研究的目标是：(1) 利用多时相Landsat影像和GEE平台上的随机森林分类，重建Ujung Pangkah红树林生态系统1995-2025年的三十年时空演变轨迹；(2) 基于设计推断的精度和不确定性估计，量化长期变化动态（面积增加、减少和稳定）；(3) 检验河流沉积物通量是否与观测到的河口动态存在定量关联。
## 二、主要技术方法
本研究基于GEE云平台，使用Landsat 5专题制图仪（TM，1995年）、Landsat 7增强型专题制图仪加（ETM+，2005年）和Landsat 8陆地成像仪（OLI，2015和2025年）的Level-2地表反射率（SR）产品。采用年度中值合成技术克服热带云覆盖并标准化潮汐波动，并用QA_PIXEL位掩码去除云和云阴影。特征工程方面，除原生多光谱波段外，引入归一化差异植被指数（NDVI）、修正归一化差异水体指数（mNDWI）、增强型植被指数（EVI）和红树林植被指数（MVI）四个光谱指数作为随机森林分类器的预测变量。使用50棵决策树的RF分类器进行“红树林/非红树林”二元分类。精度评估采用分层随机抽样，每期独立生成400个验证点（每类200个），利用混淆矩阵报告总体精度（OA）、生产者精度（PA）、用户精度（UA）和Kappa系数，并采用Olofsson等提出的误差校正面积估计法计算面积及其95%置信区间。变化检测限于两期均有效分类的像素交集，以避免数据空缺被误判为土地覆盖变化。2015-2025年损失像素的空间归因采用JRC全球地表水数据集中≥90公顷连通水域掩膜，以200米欧氏距离区分海向侵蚀和陆向人为转化。沉积物驱动分析以红波段地表反射率作为总悬浮固体（TSS）浓度代理指标。
## 三、主要研究结果
### 1. 分类性能与精度评估
基于2025年模型的随机森林特征重要性分析表明，NDVI（贡献值8.94）和EVI（7.17）是主要判别参数，其次为MVI（6.21）、近红外波段（5.64）、mNDWI（4.42）和短波红外波段（4.22），可见光波段贡献最低。消融实验显示，剔除MVI后分类性能低于仅用原始波段的基线模型（OA = 87.65%，κ = 0.751 vs. 88.00%，κ = 0.758），表明MVI捕获了植被指数无法提供的判别信息。四个时期独立验证的总体精度均超过92%，Kappa均高于0.84。各时期红树林类别的生产者精度显著低于用户精度（2005年最为显著：51.96% vs. 92.46%），表明系统性遗漏误差使误差校正面积在每个时期均高于原始制图面积9-78%。
### 2. 红树林面积时空动态
红树林面积呈现非单调演变轨迹：1995年基线面积为898.88公顷（误差校正1261.11 ± 374.24公顷）；2005年增至998.97公顷（校正1777.52 ± 518.42公顷）；2005-2015年间经历最大规模扩张，2015年达到峰值2410.44公顷（校正2626.03 ± 352.03公顷）；2015-2025年缩减至1816.87公顷（校正2213.71 ± 358.46公顷）。整个30年记录中，生态系统仍实现从898.88公顷到1816.87公顷（制图面积）的显著净扩张。
### 3. 分类后变化轨迹
2005-2015年是记录中净增加最大的时期，而2015-2025年是唯一出现净损失的时期。空间归因分析显示，在可归因的损失面积（782.84公顷，占总损失的57.6%）中，仅21.9%（171.35公顷）发生在距大型水体200米范围内，78.1%（611.49公顷）发生在更内陆的区域。这表明2015-2025年的缩减主要归因于内部/陆向转化而非海向侵蚀。内陆损失区域呈现出规则的矩形几何网格图案，与人工养殖塘特征一致。
### 4. 面积变化的幅度与统计显著性
Mann-Kendall检验得到Kendall's τ = 0.667（S = 4），双侧精确p值为0.333，未达到常规显著性水平。Sen斜率估计为+35.75公顷/年。基于误差校正面积估计的独立z检验显示：2005-2015年扩张（z = 2.47，p = 0.014）和整个30年净增加（z = 3.60，p < 0.001）具有统计稳健性；1995-2005年增加和2015-2025年缩减虽方向一致，但未达到统计显著性，应谨慎解读。
### 5. 水动力驱动趋势
总悬浮固体（TSS）代理指标（距河口3公里内开阔水域红波段反射率）的均值分别为0.0895（1995）、0.0975（2005）、0.0723（2015）和0.0601（2025）。同期相关分析显示TSS代理与制图红树林面积之间的相关性较弱且不显著（r = -0.77，p = 0.23）。然而，十年初始TSS代理与该十年内发生的净面积变化之间存在强正相关（r = 0.92，n = 3）：2005-2015年初期升高的沉积物代理预示了该时期的大规模净增加（+1395.16公顷），而2015-2025年初期的相对降低的沉积物代理则预示了记录中唯一的净损失（-565.49公顷）。受限于仅三个十年转换，该相关性无法作为统计确认证据（p = 0.25），但为河流沉积物沉积与红树林定殖之间的多年滞后效应提供了定量支持。
## 四、讨论
### 长期向海推进与新出现的人为收缩并存
Ujung Pangkah生态系统的30年净扩张相对于全球热带红树林衰退的主流模式而言是一个异常案例，支持局地地貌过程可以维持红树林韧性这一更广泛的论点。然而，净扩张掩盖了时间复杂性的存在。两个独立检验支持将2015-2025年的缩减视为真实生态信号而非数据伪影：其一，将变化检测限制在两年均有效分类的像素中，仅使估计损失变化了2.2%；其二，潮汐阶段代理指标显示2025年的平均潮位是四个时期中最低的，这反而可能高估而非低估2025年面积。
### 系统内的双重动态模式
Ujung Pangkah的海向边缘与亚马逊河口和孙德尔本斯部分地区等以沉积物为主的系统存在功能相似性，河流沉积物供给持续超过海岸边缘的侵蚀和沉降压力。然而系统内部日益类似于人为驱动的三角洲。独立研究证实，2006年至2016年间仅Ujung Pangkah分区就有101.70公顷红树林被转为鱼塘。在国家尺度上，水产养殖转化是印度尼西亚红树林丧失的最大单一驱动因素，仅2009年至2019年间就造成超过80,000公顷的森林损失。
### 河流沉积物通量与向海推进
Bengawan Solo河是悬浮颗粒物的持续来源，河流到达Madura海峡后流速减缓与潮汐混合引发絮凝和沉积。沉积形成的潮间滩涂在达到临界高程后成为先锋物种（如Avicennia属和Sonneratia属）的定殖基质。滞后相关分析（r = 0.92）为这一机制提供了定量支持。研究局限在于红波段地表反射率是悬浮沉积物浓度的间接光学代理，可能受有色溶解有机物、叶绿素a和太阳耀斑等其他水体成分影响。
### 理论贡献与蓝碳意义
维持自然水文连通性和河流沉积物输送是红树林净存续的必要条件，但显然并非充分条件。Ujung Pangkah的政策窗口应从预期性（防范假设的未来扰动）转向纠正性（应对正在进行的、空间可识别的损失模式）。地方和国家空间规划应优先保障积极向海推进的海岸前沿的划定和保护，同时将执法注意力集中于当前养殖塘和聚落转化集中的内部区域。
## 五、研究结论
本研究利用Google Earth Engine云计算框架成功重建了Ujung Pangkah红树林生态系统1995-2025年的三十年时空演变轨迹。通过多指数特征工程方法与Olofsson误差校正面积估计的整合，随机森林分类在全部四个独立验证时期实现了92.25-95.00%的总体精度（Kappa 0.845-0.900）。与全球河口退化趋势相反，Ujung Pangkah生态系统在完整30年记录中呈现统计学显著的净扩张，从1995年的898.88公顷增至2025年的1816.87公顷（误差校正后从1261.11增至2213.71公顷；z = 3.60，p < 0.001）。然而该轨迹并非单调：面积至2005年温和增长，至2015年大幅扩张至2410.44公顷峰值，随后至2025年净缩减565.49公顷——这是序列中唯一的十年净损失。2015-2025年缩减的空间归因表明其以陆向成因为主：78.1%的可归因损失发生在与养殖塘和聚落转化一致的内部区域，而21.9%发生在开阔海域或主要河道附近。同时，沉积物（TSS）代理与十年净变化之间的滞后相关（r = 0.92，n = 3）为Bengawan Solo河沉积物输送持续维持向海推进的假说提供了定量支持。沉积驱动的海岸前沿 accretion 与内部新兴的人为 attrition 两个过程并行运作，二者的净平衡将决定生态系统未来轨迹及其对区域蓝碳封存的贡献。维持生态韧性需要海岸管理政策既保护新生的向海推进前沿，又直接应对河口内部已在进行的人为转化——这是一个纠正性而非纯预期性的政策优先事项。
