一种高分辨率方法用于绘制和预测盐沼动态：整合LiDAR数据与时空变化潮汐基准面

《Ecological Informatics》：A high-resolution methodology for mapping and predicting salt marsh dynamics: integrating LiDAR data and spatiotemporal varying tidal datums

　　

准确评估盐沼动态对于扭转其广泛衰退至关重要。以往的建模研究考察了盐沼对平均海平面上升的恢复力，但忽略了在许多河口观测到的潮汐动态。基于盐沼出现由土壤高程相对于潮汐基准面（tidal datums）决定这一原理，本研究旨在通过将高分辨率LiDAR数据与空间变化的潮汐基准面相结合，绘制当前盐沼范围并预测其长期演变。以正在经历潮汐放大的潟湖Ria de Aveiro作为案例研究。研究人员开发了一个简化的、基于规则的长期盐沼演变模型，其中垂直加积（vertical accretion）作为潮汐基准面的函数内生推导得出。该模型通过重建1987年盐沼范围进行后报模式验证，并应用于预测至2100年的盐沼动态。蒙特卡洛分析（Monte Carlo analysis）量化了预测不确定性，并区分了加积与基准面变异性各自的影响。土壤高程-基准面关系被证明对制图有效，总体精度达到90.3%。该模型重现了1987年范围，并捕捉到1987年至2024年间盐沼面积26%的减少（从34.8 km²降至25.7 km²）。预计这一减少将持续至2100年，中位面积达到17.5 km²，部分原因在于阻碍向陆迁移的人为屏障。预测面积受加积控制，而基准面变异性在内陆区域变得决定性。未来研究应将本模型应用于其他河口系统，特别是那些盐沼受到潮汐放大威胁的系统。

**论文解读：基于LiDAR与时空变化潮汐基准面的高分辨率盐沼动态制图与预测方法**

**研究背景、问题与意义**

盐沼作为陆海过渡带的动态生态系统，提供栖息地、海岸防护、碳汇等关键服务，但全球范围内其面积持续衰退。以往研究多关注盐沼对平均海平面上升（mean sea level rise, MSLR）的响应，忽略了潮汐动态（如潮差变化、潮汐放大）对盐沼分布的影响。潮汐基准面（tidal datums，如平均高潮位MHW、平均大潮高潮位MHWS等）直接决定盐沼的淹没频率与空间分布，而许多河口因人工疏浚、航道加深导致潮汐传播改变，潮汐基准面呈现时空变异。现有模型如SLAMM（Sea Level Affecting Marshes Model）和SMRM（Simplified Marsh Response Model）或依赖固定加积率、或假设均匀潮汐基准面，难以适用于潮汐动态演变的系统。因此，亟需开发一种计算高效、能反映潮汐基准面时空变化的高分辨率盐沼动态预测模型。本研究以葡萄牙Ria de Aveiro潟湖为案例，该潟湖已观测到潮汐放大导致的盐沼退化，但缺乏高精度现状分布与未来演变评估。研究旨在：1）利用LiDAR（Light Detection and Ranging, 激光雷达）数据与空间变化潮汐基准面绘制当前盐沼范围（2 m分辨率）；2）开发内生加积模型预测至2100年的盐沼动态，并定量不确定性来源。论文发表在《Ecological Informatics》。

**关键技术方法**

研究人员采用以下主要技术方法（忽略具体培养、质粒构建等步骤）：1）高分辨率LiDAR数字地形模型（Digital Terrain Model, DTM，2 m分辨率，来自葡萄牙国土总局DGT）提供土壤高程；2）基于Delft3D流体动力学模型（Lopes等，2022校准验证）模拟2024年及1987年潟湖内各水体的潮汐基准面（MHWN、MHW、MHWS、HAT）；3）发展简化规则模型：垂直加积分为有机加积（基于生物量-高程抛物线关系）和无机加积（基于高程-加积率S型函数），加积率由潮汐基准面内生驱动；4）后报验证：基于1987年空间可变潮汐基准面与参数敏感性分析（576种组合），重建1987年盐沼范围并与航片对比；5）蒙特卡洛集合（Monte Carlo ensemble, 1500次模拟）预测2024-2100年演变，考虑潮汐基准面空间变异与加积参数不确定性。

**研究结果**

**4.1 盐沼制图**
**4.1.1 潮汐与土壤高程的局部特征**：潮汐基准面空间变异小（标准差≤0.10 m），但WB5（上游水体）的MHWN比WB1（入口）高0.07 m。LiDAR显示60%潮滩高程在0.6-1.8 m间，WB2潮滩面积最大（12.4 km2），WB1最小（1.1 km2）。各水体高差曲线（hypsometric curves）表明WB2和WB5在MHWS-HAT间坡度陡峭，限制高沼泽空间。

**4.1.2 盐沼范围评估与验证**：通过分层随机抽样（1000个验证点）与高分辨率正射影像（OrtoSat2023）验证，总体精度达90.3%，生产者精度（盐沼类）88.7%，用户精度81.2%。当前盐沼面积25.7±0.5 km2，中沼泽（mid marsh）占82.5%（21.2 km2），低沼泽（low marsh）7%，高沼泽（high marsh）10.5%。WB2面积最大（12.1 km2），WB1最小（1.4 km2）。

**4.2 盐沼对潮汐与海平面变化的恢复力**
**4.2.1 潮汐基准面与平均海平面的年际变化**：1976-2023年验潮站数据显示各潮汐基准面线性上升，HAT速率7.8 mm/yr，MHWS 6.7 mm/yr，MHW 5.9 mm/yr，MHWN 5.0 mm/yr，平均海平面（MSL）上升速率3.4 mm/yr。潮汐基准面上升速率高于MSL，表明潮汐放大。

**4.2.2 后报盐沼动态（1987-2024）**：模型对最大无机加积率（Δzi,max）敏感，最优参数组合（上四分位数，149组）总体精度≥74.7%。重建1987年总面积34.8 km2 [32.9-35.2]，中沼泽损失最大（-4.8 km2），低沼泽与高沼泽相对减少54%和40%。1987-2024年面积下降26%，WB5损失最大（50%），WB2损失18.2%，WB1、WB3、WB4变化较小。

**4.2.3 预测盐沼动态至2100年**：中位面积17.5 km2，较当前减少约32%。低沼泽和高沼泽几乎完全消失，中沼泽在WB1略有增加（0.6→0.8 km2），WB2和WB5下降。预测面积受垂直加积率主控，所有水体中加积率越高面积越大，但即使加积率超过MSLR，面积仍低于当前。潮汐基准面空间变异在WB5起决定性作用（MHWN变异0.4 m时面积接近0），其余水体影响较小。海岸挤压（coastal squeeze）——人为屏障阻碍向陆迁移——是限制盐沼存续的关键因素。

**讨论与结论**

**讨论总结**：制图方法精度90.3%，与同类研究一致，误分类主要发生在过渡带与盐田区。潮汐基准面空间变异揭示WB5潮汐放大最显著（1987-2024年MHWN上升35 cm，MHWS上升53 cm），与航道加深导致的潮波传播改变一致。后报模型虽低估某些区域面积，但成功再现了盐沼大面积丧失的空间模式（中沼泽损失为主）。预测显示即使加积率较高，盐沼面积仍下降，主要因高差曲线陡峭导致高海拔空间不足，且人为屏障限制横向迁移。模型不确定性来源包括：未来海平面与潮汐轨迹、加积率参数（缺乏本地实测）、模型未包含横向侵蚀（边缘侵蚀与潮沟扩张）以及植被动态（如物种竞争与盐度影响）。模型计算效率极高（75年模拟约15秒/19,575,996个网格），使得蒙特卡洛不确定性分析可行，并易于推广至其他系统。与SLAMM和SMRM相比，本模型能解析潮汐基准面的年际与亚区域变异，且仅需三种简单沼泽类型，更具可转移性。

**结论部分翻译**：本研究提供了Ria de Aveiro盐沼当前分布的高空间分辨率评估，并评估了其对海平面上升的恢复力。当前面积为25.7±0.5 km2，以中沼泽植被为主（82.5%），低沼泽和高沼泽分别占7%和10.5%。中央潟湖（WB2）承载45%的盐沼面积，而Mira水道（WB1）仅占5%。至2100年的预测表明，总盐沼面积将下降，中位损失约36%。这一下降伴随着沼泽类型的重新分布：低沼泽和高沼泽几乎完全消失，而中沼泽仍为主导类型。垂直加积成为盐沼存续的一阶控制因素，预测面积随加积率增加而增加。潮汐基准面的空间变异施加了次要的、空间异质的影响，在内陆水体（WB5和WB2）中更为重要。这些盐沼有限的适应能力因海岸挤压而加剧：低海拔沼泽丧失后，向陆迁移被人为屏障所阻碍。预测模型通过重建1987年盐沼分布进行了后报验证，垂直加积由高程与演变基准面之间的相互作用内生推导得出。其计算效率使得蒙特卡洛集合能够显式传播参数与基准面不确定性。这些特性使得该框架可复制并转移至其他河口系统，特别是那些正在经历潮汐放大的系统。此处提供的详细表征为设计旨在保护盐沼及其所支持生态系统服务的监测与管理计划提供了基线数据。鉴于预测损失集中在内陆水体，监测工作应优先部署在那里。在国家尺度上，该方法可应用于所有葡萄牙河口系统，利用DGT提供的LiDAR衍生DTM以及现有流体动力学模型，根据水框架指令评估过渡水体。