综述：在茎秆尺度上表示柔性盐沼植被：生物力学、建模方法及研究空白

《Coastal Engineering Journal》：Representing flexible salt marsh vegetation at the stem scale: biomechanics, modeling approaches, and research gaps

　　

在茎秆尺度上准确表示柔性盐沼植被对于提高研究人员对流体-植被相互作用（fluid–vegetation interaction）及其在海岸防护和其他生态系统服务中作用的理解至关重要。本文批判性地回顾了在物理实验和数值建模方法中用于复制柔性盐沼植被的现有方法。具体而言，该综述（i）评估了哪些生物力学茎秆属性（biomechanical stem properties）对于在水动力强迫（hydrodynamic forcing）下再现真实的植物行为最相关，（ii）考察了以往研究如何实现柔性茎秆的代用和数值表示，以及（iii）综合了当前与准确复制植被力学相关的研究空白（research gaps）。特别关注了诸如茎秆几何形状（stem geometry）、抗弯刚度（flexural stiffness）、材料密度（material density）及其在实验和计算框架中的缩放（scaling）等参数。通过识别不一致性、方法上的权衡和未解决的挑战，本文为工程师和建模者选择合适的表示策略提供了指导。同时，本文告知生态学家哪些现场测量的植物属性对于支持可靠的植被-流动建模（vegetation–flow modeling）工作最为关键，从而促进生态实地调查与水动力建模研究之间更强的整合。

**1. 引言**
本文指出，全球沿海人口持续增长，增加了对风暴潮和海平面上升等海岸灾害的暴露程度。传统灰色基础设施虽能提供短期防护，但常干扰自然海岸动态。盐沼作为基于自然的解决方案（NbS）的典型代表，全球面积约5.5 Mha，通过植被诱导的流动减缓促进沉积物积累。然而，盐沼植被的物种组成、形态和生物力学特性存在区域和季节差异，导致基于单一地点测量值的模型难以推广。本文旨在通过综述现有方法，识别茎秆尺度下柔性盐沼植被表示的关键生物力学属性，评估建模方法的可重复性，并指出研究空白。

**2. 盐沼在海岸防护中的作用**
盐沼植被通过茎秆和叶片增加水力粗糙度，降低流速并促进沉积，从而减弱波浪能量并增强土壤稳定性。波浪衰减能力取决于淹没比、迎风面积（frontal area A）、植被密度和入射波条件，例如在40 m长的实验段中波浪高度可减少高达60%。侵蚀减少方面，地上和地下生物量联合作用可降低侵蚀率高达80%，根系网络增强土壤抗剪强度。垂直沉积方面，单个茎秆改变局部流场，促进沉积，植被区域沉积速率比裸泥滩高42%。

**2.1 波浪衰减**
波浪衰减源于植被诱导的流动阻碍，耗散波浪能量并减小冠层内的轨道速度。衰减幅度受淹没比（d/LS）、迎风面积A、植被密度和入射波条件影响。即使在风暴潮条件下，水槽实验也显示波浪高度可减少高达60%，但早期生长季节的衰减效率降低，表明需要量化季节变化。

**2.2 侵蚀减少与垂直沉积**
地上和地下生物量联合作用可减少侵蚀率高达80%。根系网络增强土壤抗剪强度，提供全年抗侵蚀能力。植被茎秆通过改变局部流动模式促进沉积，例如两周内植被区域沉积量比裸泥滩高42%。

**3. 物理过程与建模考虑**
**3.1 茎秆灵活性**
盐生植物的茎秆灵活性允许其随波运动，降低流体动力学载荷。柔度影响波浪衰减效率，研究报告降低幅度在30%至70%之间。柔性植被表现为回避行为，随流弯曲，减少局部湍流；而刚性茎秆能承受更高流速，可能增强沉积但也会增加湍流和再悬浮。

**3.2 研究尺度**
盐沼植被动力学可在草甸、斑块或茎秆尺度上研究。波浪衰减通常在草甸或斑块尺度量化，而单个茎秆动力学影响局部沉积物移动。茎秆灵活性在草甸尺度模型中的整合仍具挑战性，例如有效长度（le）概念用于考虑植物重构。简化方法如有效植物高度和物种无关的阻力系数（CD）与雷诺数（Re）关系可用于预测波浪衰减，但无法解析茎秆摆动等过程。

**3.3 实验与数值模型中的植被表示**
现场研究因水柱内能见度有限而稀少，受控水槽实验需真实植被表示。活体植物采样存在生态和后勤挑战，人工代用或数值表示提供灵活替代，但许多简化方法依赖刚性元素，无法再现植物重构或生物力学属性的季节变化。

**4. 植被属性的复制**
**4.1 几何表示**
茎秆应复制自然植被的茎长（LS）和直径（dS）或宽度（WS）及厚度（TS）。茎秆横截面形状影响惯性矩（I），例如圆形截面I=π/64 dS⁴，矩形截面I=WS TS³/12。对于不规则形状，需使用校正系数（如Elymus spp.的cElymus=0.79）计算修正惯性矩Ic。迎风面积和抗弯刚度（EI）需与活体植被匹配。

**4.2 相似性参数**
当无法精确匹配所有属性时，需确保关键无量纲数的相似性，包括柯西数（Ca=ρ u² bS LS³/EI）、浮力参数（B=(ρ-ρS) g AS LS³/EI）、雷诺数（Re=ρ u dS/μ）和凯莱根-卡彭特数（KC=uw T/WS）。茎秆运动主要由Ca控制，浮力对挺水植被重要，Re确保流动条件可比。

**4.3 缩放考虑**
物理实验中常用弗劳德缩放，但会引入Re和Ca的尺度效应。弗劳德-柯西缩放通过调整杨氏模量E来保持Ca相似性，但Re可能不匹配。替代缩放法则（如Attili等提出的）可匹配结构宽度缩放。对于草甸，可通过迎风面积或固体体积分数缩放茎秆以克服Re差异。

**5. 物理植被建模中的代用方法**
**5.1 抗弯刚度的相似性**
代用材料的选择需复制活体植被的EI。例如，Markov等使用乳胶橡胶管匹配S. alterniflora的茎秆运动，但未达到精确一致。其他研究通过验证茎秆运动或比较EI来评估代用材料的适用性，但常存在偏差。

**5.2 复制挑战**
van Veelen等的研究中，硅密封胶代用茎秆的直径和EI分别匹配不同物种，但未完全复制所有特征。Keimer等使用3D树脂打印制造可变柔度代用茎秆，但杨氏模量低于现场测量值，且树脂性能受环境影响。整体而言，复制所有相关参数仍具挑战性，需探索3D打印等新技术。

**6. 植被茎秆建模中的数值方法**
**6.1 结构求解器**
常用梁或杆理论，包括欧拉-伯努利梁（仅弯曲）、铁木辛柯梁（考虑剪切）、基尔霍夫杆（弯曲和扭转）和科瑟拉杆（弯曲、剪切、扭转和轴向拉伸）。科瑟拉杆理论已成功应用于柔性植被建模，但计算成本高。铁木辛柯梁和基尔霍夫杆在植被建模中应用较少，后者被认为有潜力。

**6.2 受力方法**
半经验方法（如莫里森方程）基于参数化力表达式，计算效率高但依赖经验系数（如CD、CM），且未考虑复杂尾流相互作用。直接力方法从解析的流场计算力，准确性更高但计算成本高，且依赖网格分辨率和湍流模型。CD值受Re、KC和Ca影响，但缺乏针对物种和季节变化的综合数据。

**6.3 耦合考虑**
单向耦合仅考虑流体对结构的影响，双向耦合考虑相互反馈。空间耦合方法包括浸入边界法（IBM）、自适应网格法（如任意拉格朗日-欧拉方法ALE）和光滑粒子流体动力学（SPH）。IBM广泛应用于植被研究，但需注意插值误差；自适应网格法对网格变形敏感；SPH避免网格问题但计算成本高。现有FSI求解器验证范围有限，大多针对单向流或振荡流，且很少直接整合现场测量数据。

**7. 结论与研究需求**
本文强调茎秆尺度灵活性在捕获流体-植被相互作用中的关键作用。物理和数值模型需考虑茎秆几何、抗弯刚度、材料密度和季节变化。研究空白包括：标准化现场测量、改进FSI求解器效率、开发复杂几何代用材料、整合沉积物输运、以及从茎秆到草甸的升级方法。未来工作应聚焦于这些方面，以增强对盐沼生态系统服务的量化能力。