草没有连成一片，影响却可以连起来
在一片刚被草占据的海滩上，地面可能仍大部分裸露。每一丛草只占很小的位置，却不只影响自己脚下：风绕过茎叶，颗粒撞上植物、落在背风侧，附近裸地的输沙也随之改变。相邻草丛即使没有接触，也可能通过同一股风、同一批流动的沙相互影响。
经过许多次输沙事件，原来分散的小沙堆可能逐渐并成较大的沙体。地面抬高又改变风、海水接触和埋沙条件，植物随之调整生长。这里的“共同建造”指这些物理与生长过程相互作用，不需要草有共同目标。
2026 年的一项研究追踪了荷兰 Texel 岛 Hors 海岸的一块约 12 公顷地段。当地有宽阔的前方沙滩，研究区受到的直接海浪侵蚀较少。研究者识别出 2013 年的 4,349 个草丛，利用后续八个年份的高分辨率地形资料，比较早期草丛配置与 2014—2023 年的地形变化。它为这个问题提供了很好的长期窗口，但并不代表所有海岸都有相同条件。1
先把问题拆成三步：草怎样让沙留下，空间排列怎样改变这个作用，短时沉积又怎样累积为多年地形。
沙丘长高，要看沙在哪里少流走了
风沙输运中的沙粒可以跳跃、滚动或被悬起。植物一方面直接拦截颗粒，另一方面承受部分空气阻力，并改变近地风和湍动。地面原本承受的风力、颗粒再次起动的机会，以及下游能够继续携带的沙量，都可能变化。4
但“风变慢”还不是完整的沉积解释。要看一小块地面有没有积沙，必须比较流进和流出的沙量。在忽略压实变化、只讨论水平方向输沙的简化条件下：
这里 $z$ 是地面高程，$\rho_b$ 是包含孔隙影响的沙层体积密度，$\mathbf q$ 是每单位宽度的质量输沙通量。大白话是：流进来的比流出去的多，差额才能留下来抬高地面。 有草、有弱风，却没有沙源，不会凭空长出一座沙丘；一个位置截留得更多，也可能使更下游的草少分到沙。
还有两个容易混在一起的调整过程。风进入草地后，要经过一定距离才充分受到植物阻力影响；沙粒群体的输运量也不会在每一个位置立即等于当地风对应的平衡值。
2023 年的一项研究结合美国俄勒冈海岸的现场测量与已有风洞资料，专门分析了这种空间滞后。它指出，稀疏植被中的沉积区不一定从第一排草立即开始；模型如果默认植被一出现就实现全部遮蔽效果，可能把沉积位置放错。该研究的三条野外样带、短时观测与模型标定，支持认真处理这种调整长度，但不能直接提供任意草种、任意地形的统一长度。4
本库用一个守恒的小例子说明这件事。先指定风在草地中逐渐减弱，再让输沙通量 $q$ 以长度 $\ell_s$ 向当地的平衡通量 $q_{\rm eq}$ 调整：
示意中把近地摩擦速度写成 $u_*(x)=0.4[0.55+0.45e^{-x/\ell_c}]$ 米/秒，并采用 $q_{\rm eq}\propto\max(u_*-0.25,0)^3$ 的阈值关系。摩擦速度是从地面剪切应力换算的尺度，不是风速计直接读到的风速。$\ell_c$ 控制预设的风调整长度，$\ell_s=0.8$ 米控制输沙调整。这些是用于解释的参数，没有拟合 Hors 的现场风场。10
风调整长度由 0.5 米改成 2 米、8 米，计算中的沉积峰位置由约 0.29 米移到 0.67 米、1.29 米。它们并不等于各自的风调整长度，因为沉积同时受风、输运阈值和沙粒调整控制。网格加密后，归一化通量的最大差小于 $5.1\times10^{-5}$，并且逐格保持收支一致。
这个例子只解释沉积位置怎样改变。它没有证明真实草地会截留全部来沙，也没有用一种风速公式替代完整的湍流和颗粒碰撞过程。
同样多的草，为什么排列会改变结果
一大团草可能形成很强的局部遮蔽，但后排草面对的是已经被前排改变过的风沙；分散的草可以接触更多不同的输沙路径，却可能过于稀疏，使很多沙从间隙继续通过。多个小簇则可能同时兼有局部遮蔽和较宽的影响范围。哪一种更合适，取决于风向、草形、供沙和评价指标。
有一个不需要复杂流体计算的收支例子。设两条并行输沙路径各流来 1 单位沙，有两个相同的截沙单元，每个能留下经过它的沙的一半。如果两个单元串在同一路径上，留下的是 $0.5+0.25=0.75$ 单位；若分别放在两条路径上，共留下 1 单位。截沙总比例由 37.5% 变成 50%，没有增加单元数量，也没有改变每个单元的截留规则。10
这只是“共同使用同一股来沙”的示意，不是实际草丛固定有 50% 的截沙率。它说明：不能把每丛草单独测得的作用直接相加，也不能仅凭草丛数量预测总效果。
2019 年的一项野外实验提供了更直接的证据。研究者使用相同数量的人工模拟草，在 4 平方米样地中安排为分散、多个小簇和集中一簇三种构型，以随机区块和无模拟草的对照处理当地的地面变化，并在三个月份重复测量。它改变了空间排列，因此比仅比较自然界里“草多的地方和草少的地方”更接近检验排列的因果作用。5
公开的图表数据给出：以分散构型的截沙量为 1，多个小簇约为 0.672，集中一簇约为 0.307。截沙总量最高的是分散构型。 但作者还定义了一个效率指标：把截沙量除以平均相邻连接距离，作为扩展投入的几何代理。在这个指标下，多个小簇相对分散约为 2.34 倍，相对集中一簇约为 1.91 倍。6
这里的效率不是直接测得的植物能耗、存活率或繁殖成功率。人工模拟草也不会长根和改变生长方式。实验较直接支持的是“排列能够改变截沙，改变分母又会改变优劣排序”，不能由此宣布某个空间图样在一切环境中都最优。
多年的地形里，早期配置留下了什么
回到 Hors 的长期资料，研究者用草丛中心周围 7 米范围内的数量描述邻域，再除以圆形面积，得到局部密度。7 米是在多个候选邻域尺度中，根据与地形的关联及其增益变化选出的分析尺度；它不是草丛作用范围的直接测量。1
本库读取论文随附的逐丛资料，抽出 2013 年、7 米邻域的 4,349 条记录，重算得到以下关系。这里列的是 Spearman 秩相关，表示排序上的关联强弱，不是因果效应大小。3
| 2013 年的预测量 | 与 2014 年高程的相关 | 与 2023 年高程的相关 |
|---|---|---|
| 单丛面积 | 0.179 | 0.102 |
| 邻域草丛密度 | 0.573 | 0.529 |
这些结果支持：在这块研究地段，只知道一丛草有多大，远不如知道它周围的草如何分布有用。 它不等于单丛面积完全没有作用。研究者把初始地形、距海岸距离和单丛面积一起加入分析后，邻域密度仍提供更多解释；在近似相同的邻域总植被面积内，分成更多草丛也与更高的后期地形相关。这些是对竞争解释有价值的检查。2
但研究区不是随机分配的沙滩。微地形、实际来沙方向、局部水分与建立历史不可能被几个变量全部代表；相邻记录还会共享邻域和地形，因此 4,349 条记录不能当成 4,349 个相互独立的实验。这里的“预测”指早期资料与后来数据的关联，尚不等于在另一个海岸验证过的预测能力。本文重算相关和曲线，没有重新完成空间区块统计或航测处理。
还有一个很有意思的历史信号：用 2013、2017、2021 年各自的草丛密度预测 2023 年高程，本库复算的相关分别约为 0.529、0.276、−0.007。三批草丛数量分别为 4,349、2,908、3,188，后来的草丛边界和数量已发生变化。3
这支持早期结构留下了持续痕迹。但“后来数草丛不再好用”不等于后来植物生长不重要。草丛扩展、合并后，原来许多小丛可能被识别成少数大片；同一个“丛数”指标在不同阶段代表的结构已经不同，地形也保留了过去多年累积的沙。
一条 S 形曲线里，哪些是数据，哪些是解释
这项长期研究发现，低密度区的平均高程增加有限，中间一段密度对应较快的高度变化，再往上逐渐趋平。作者用 Hill 函数描述它：
其中 $x$ 是局部草丛密度，$y$ 是高程，两个平台值和形状参数都来自拟合。它是一条描述关系的曲线，方程本身没有风、沙或植物生长。
本库按每个密度箱等权重拟合八个观测年份，重建了公开的曲线；在原图数据位置，预测高程与发布文件的最大差约为 0.000013 米。这个很小的差只说明曲线计算对上了，不是航测具有这种测量精度。若改成每条草丛记录等权，2023 年的 $k$ 由约 0.0541 变为 0.0571，$h$ 由约 4.32 变为 3.72。稀少密度箱获得多少权重，会影响阈值样貌。310
参数的数学含义也需要准确。由方程直接可得，$x=k$ 时，高程处于上下平台的正中间；在线性密度横轴上，真正最陡的拐点是：
2023 年的公开参数给出 $k\approx0.05406$，但 $x_{\rm inf}\approx0.04847$ 丛/平方米。正文把 $k$ 称作拐点并不严格；在这里将它称为“半高密度”更准确。这不消除 S 形关系，却提醒我们：从一个拟合参数到一个真实的物理门槛，还需要额外论证。
本库还没复现补充表里的 AIC 原值。以公开的 31 个密度箱均值、四参数 Hill 曲线及常规高斯残差似然计算，2014 年和 2023 年的 AIC 约为 −24.98 和 34.65，补充表为 −48.52 和 4.66。公开曲线的残差也给出同样的差异；原 R 脚本未能取得，原因保持未决。因此这里只把 S 形当作已能复算的描述，不宣称独立验证了整张模型优劣表。210
约 4.5 米，是怎样得到的
作者提出：每个草丛的影响可近似画成一个圆，两个圆重叠就可以相互作用；当许多圆连接成跨越区域的群组时，整体建造能力增强。这个思路借用了连续渗流理论。1
论文先把约 0.05 丛/平方米的半高密度解释成连通转变附近，再结合圆盘模型和草丛聚集程度，推得影响半径约 2.27 米；两个半径相加，得到约 4.54 米的中心间作用距离。之后，它检查了最近邻距离与高程、按该距离构图的跨越组比例，以及分析窗口改变时拟合密度的变化。12
这些检查提供了支持，却不是三次独立测得了同一个风场长度：它们仍来自同一片草丛与地形资料，其中构图还使用了已经推定的距离。将这种距离当成现场流场仪器量出的结果，证据就被放大了。
圆盘理论中还有两个需要分清的量。对于均匀泊松分布、允许重叠且半径相同的圆盘，填充参数为 $\eta=\rho\pi r^2$。文献给出的临界值约为 1.128；但一个圆与另一个圆重叠，要求两者中心相距不超过 $2r$，所以平均相连邻居数为 $4\eta$，在这一模型的临界点约为 4.51，不是 1.128 个邻居。8
2026 年论文在这一处把填充参数写成平均相连邻居数；其聚集修正也需要谨慎。如果 $G(r)$ 指通常的局部对关联函数，半径 $R$ 内的平均邻居数涉及 $2\pi\rho\int_0^R sG(s)\,ds$，不能一般地用边界处的 $G(R)$ 乘圆面积替代。相关完整实现没有取得，这里指出的是书面定义与推断仍待澄清，不能据此认定作者程序一定采用了哪种错误实现。
本库另做了几何例子：40 米见方的区域放入同样 100 个点、使用同样的 2.27 米影响半径。规则排列能全部连接；分成彼此较远的四团时，每团内部连接很好，最大群组却只有全部点的四分之一。密度完全相同，整体连接可以非常不同。10
真实风场还具有方向、强弱、植被孔隙和高度差，圆形、等半径、二值连接只能是一个近似。沙粒也不需要沿一条“草丛接草丛”的网络才能前进；它能飞过裸地。几何上连起来，必须再接上改变风沙通量的动力学，才算解释了建造过程。
对窗口尺度的检查也应如此理解。按论文选择的 $k(L)=k_c+aL^{-0.75}$ 重算，$R^2\approx0.944$；换一个示例指数 1.5，同样两参数拟合也可达到约 0.973。这个比较不证明后者是正确的物理规律，只说明“看起来很直、拟合很好”不足以单独识别二维渗流的普适指数。10
因此，较稳妥的判断是：连通式的相互作用是有物理基础、也有数据支持的候选解释；约 4.5 米的精确含义和渗流是否为主导原因，尚未被这些资料唯一确定。 同样不能反过来说，所有邻丛效应都只是巧合。
地形会改变草，最初的排列不会永远固定
如果把每丛草当成不会变化的挡板，就只能解释第一段过程。活的草会生长、扩展和衰退，而埋沙本身就是改变这些行为的因素。
2021 年的一项研究比较了不同供沙地点，并进一步在同源植物的室外实验中改变沉积条件。它发现，缺沙条件下的扩展更集中，有沙沉积时形成的空间布局更分散为多个簇。较高的埋沙水平却没有简单带来更多枝条：强埋沙处理的整体枝条数量明显减少，部分个体生长很差，而只分析枝条足够多、能够拟合空间分布的个体会漏掉这部分受损者。7
这为反馈增加了两层内容。第一，植物配置可以是环境作用的结果，因此不能把每一种配置都当成先于环境、固定不变的原因。第二，埋沙可能促进向上或侧向扩展，但超过耐受能力也会压制生长；正反馈有条件，不能写成“沙越多—草越好—沙丘永远越高”。
地形继续抬高后，迎风坡和背风区会变化，草丛可能合并，沙源也可能被前方的新沙体截走。海浪侵蚀、风暴和植被受损又能重新打开输沙通道。最终形态因此来自一段历史，而不是把某一年草的密度代入曲线就永远确定。
怎样进一步判断一片真实沙丘由什么控制
最有区分力的证据应把空间结构与实际通量接起来：在总植被量、初始地形和供沙近似相同的地段改变排列，同时跟踪近地风、进入和离开的沙量，以及草丛之间的地形变化。若草丛间的影响确实重叠，能否在风和通量中先看到这种重叠，再看到共同沙体形成，会比仅从最终高度反推更有说服力。
其次要改变方向和尺度。旋转同一构型、改变主要来沙方向，或者使用孔隙与高度不同的草丛，若推定的相互作用距离随之变化，就能判断“一个固定圆半径”在哪些条件下适用。比较相同总面积的单大丛、多小丛与不同间距，还能分开总遮蔽、邻丛相互作用和上游抢先截沙。
长期观测则需要持续记录植物的扩展与合并，不能只追踪“丛数”；还应区分高程、相对初始地面的增高和总沙体积。本篇数据里的变化列与直接用两个年度高程做差，在部分记录上并不完全相同，本库没有取得足以解释差异的处理链。两种算法与早期密度的正关联接近，但这不允许把它们当作逐条相同的测量。310
现有证据已经足以说明：草丛能通过近地流动和输沙影响比自身更大的区域；空间排列会改变截沙效果；地形和埋沙又会改变植物布局。把这三步接起来，可以理解分散的草怎样参与建造一片沙丘。尚未收紧到同一答案的是：在某个具体海岸，供沙、单丛形态、邻丛相互作用和长期生长各占多大分量，以及哪一种空间模型能可靠预测它们随条件变化的结果。