解读《美国科学院院刊》2025 年研究：六十年的一道石笋形态方程，终于解出来了
进洞看石笋，第一件让人纳闷的事不是它们有多老，是它们长得太不一样了。
有的细得像一根蜡烛，顶上还是尖的。有的粗成一截柱子，顶圆滚滚的。还有的顶上平平的，像有人拿刀切过。同一座洞里，隔着三五米，就能同时站着这三种。好看的那些还被人起了绰号，烛台、女巫手指、婚礼蛋糕，还有罗密欧与朱丽叶。
图1 十二根真实石笋。A 烛台（捷克 Sloupsko-šošůvské 洞）、D 尖塔（澳大利亚 Jenolan 洞）、G、K 婚礼蛋糕（美国卢雷洞）是锥形；B、C 女巫手指（美国卡尔斯巴德洞窟）、E、H、L 是柱状；F、J 是平顶（斯洛文尼亚 Postojna 洞）；I 是捷克 Punkva 洞里上下相对的「罗密欧与朱丽叶」。引自原文 Figure 1。
这些形状不是随便长的。它们背后有一套相当确定的物理：水从洞顶滴下来，落在石笋顶上，然后沿着表面往下流，一路流一路把碳酸钙沉在石头上。既然每一步都能写方程，那形状就应该能算出来。
一、六十年前，有人猜了一个答案
1960 年代，一位叫 Franke 的研究者给石笋生长写了一个模型，简单到只有两句话。第一句，石头长得有多快，正比于水里的钙有多「过饱和」。第二句，水从顶上流下来，一路沉淀，越往下钙越少。写完他顺手猜了个结果：条件不变的话，石笋迟早会长成一个固定的形状，之后只管往上平移，轮廓再也不变。
这个猜测后来被计算机模拟一遍遍证实了。模拟出来的确会收敛到一个稳定的柱状。但那个形状的数学表达式，六十年里没人写出来过。
为什么难？因为这是一个自己咬自己尾巴的问题。表面是什么形状，决定了水怎么流；水怎么流，决定了钙在哪里被消耗；钙在哪里被消耗，又决定了下一层石头长在哪里。三个环节扣成一个环，方程是非线性的，硬解解不动。
二、把方程解开，只剩一个数
这篇论文的破题办法不是硬算，是换一个问题问。
既然要找的是「形状不变、只往上平移」的那种解，那就先把这个条件当成前提用起来。这一步看着像偷懒，其实是这类问题唯一走得通的路。黏性指进里的 Saffman-Taylor 指形、凝固里的 Ivantsov 抛物面、冰柱、石灰华锥，都是这么解出来的。
前提一代进去，方程立刻松了。钙浓度沿半径是一条抛物线，从顶上最高，往外一路降，降到某个半径 R 的时候正好等于饱和浓度。再往外，水里的钙已经沉不出东西了，石头也就不再往外长。所以石笋有一个确定的半径。
把这个半径、这个浓度、这个生长速度凑到一起，所有参数会自己塌成一个无量纲数：
Da = kA / Q。
k 是方解石沉淀的速率常数，A 是石笋的底面积，Q 是滴下来的水的体积流量。这个东西叫 Damköhler （达姆科勒）数，化工里常用，意思是「反应有多快」比上「送料有多快」。
换成大白话：水给得少，一滴水落在顶上，钙还没流到边上就沉完了，顶就摊平；水给得多，钙来不及沉，被水一路带到边上才沉，形状就被冲成了尖的。中间那个不多不少的点，是柱状。
图2 三种理想形状。左，平顶（Da 大于 1）；中，柱状（Da 等于 1）；右，锥状（Da 小于 1）。表面的网格是等间距的母线和纬圈，能看出顶部曲率的差别。底下的 U 是往上长的速度，U0 是顶面水平时的生长速率。引自原文 Figure 3。
三、三种形状，洞里都找得到
把 Da 代回方程，得到的不是一条曲线，是一族曲线。按 Da 的大小分成三段：
Da 大于 1，平顶。方程在中心那一圈是没有解的，因为水里的钙一落地就沉光了，中心反而长不出斜面，只能是一块平盖。
Da 等于 1，柱状。顶是半圆的，身子是直的。这是唯一一个前人算出来过的形状，而且只是数值解。
Da 小于 1，锥状。顶是尖的。
三种形状往上长的速度也不一样：平顶慢于 U0，柱状正好等于 U0，锥状快于 U0。也就是说，尖的石笋不只是长得尖，它长得也快。
图1 里那十二根，A、D、G、K 是锥，B、C、E、H、L 是柱，F、J 是平顶。三类都在真实洞穴里，不用去凑。
四、可方程没告诉你，该落在哪一个
到这里有个麻烦。方程给出一整族解，从很尖到很平，全都能满足方程。但真实的一根石笋只会是其中一种。那么，是什么在挑？
这个问题在物理学里有个名字，叫形状选择。凝固、黏性指进里都有，而且都知道：光靠方程本身选不出来，必须再往模型里塞一点新的物理。
平顶的选择规则，作者找到了，而且很具体：水滴是从多高掉下来的。
一滴水从十米高的洞顶落下来，落到石笋顶上的时候早就不是一滴水了。它会溅开，铺成一片薄薄的水膜，覆盖的面积比水滴本身大得多。落得越高，溅得越开。所以真实的情况不是「一个点源」，是「一小片均匀洒水的圆盘」。
把点源换成圆盘，重新解一遍，答案变了：Da 不再等于 1，而是 Da = 1 + k/P，P 是单位面积的洒水速率。平顶的半径、平顶石笋往上长的速度，全都由这一个量定死。
这解释了一件事。斯洛文尼亚的 Postojna 洞里，矮的石笋大多是平顶，高的石笋顶多是圆的。看着奇怪，其实是落距：矮的那些，顶离洞顶远，水滴落得久、溅得开，摊成了平顶；高的那些，顶快顶到洞顶了，水滴几乎是直接落上去的，溅不开，就成了柱状。
柱状石笋是这套解释的极限情形：湿区小到可以忽略，等于回到点源。这也顺带说明了柱状为什么最常见。
锥形就难了。湿区大小解释不了尖顶。作者换了个方向，去看洞顶那一位。
五、锥形石笋的答案，藏在钟乳石身上
钟乳石和石笋常常是一对，同一滴水，上面沉一点，下面沉一点。可它们长得完全不像。
钟乳石细长，像冰柱。石笋粗壮，像个墩子。如果是同一套物理的上下镜像，那形状应该对称才对。不对称，说明两边沉淀的机制根本不是一回事。
作者的说法是：差别在水里的二氧化碳。
水在土里的时候溶了很多生物呼吸产生的二氧化碳，偏酸，能溶下很多钙。进了洞，洞里空气的二氧化碳分压低，水里的二氧化碳开始往外跑，钙就沉出来。
钟乳石那边，限制速度的是「碳酸变成二氧化碳再跑掉」这一步。膜薄的地方，总水量少，表面张力又牢牢拽着水分子，二氧化碳的逃逸效率低，钟乳石长得就慢。冰柱那套理论就是这么算的。
石笋那边不一样。水已经挂在洞口滴了半天，落下来又溅开，二氧化碳早就跑得差不多了。这时候限制速度的是石头表面的反应本身，正比于过饱和度，跟水膜的厚薄基本无关。水一溅开，钙就全方位往外扩张。
切换点在哪里？就在水滴悬挂、坠落、溅开那一瞬间。
这个说法有一个可以当场验的例子：捷克 Punkva 洞里那对「罗密欧与朱丽叶」（图1 的 I）。上面是钟乳石，下面是石笋。石笋的上半段，轮廓确实像是把钟乳石倒过来照了一遍；再往下，才慢慢变粗，变成典型的锥形。中间那个转折，就是水从「还没平衡」变成「已经平衡」的位置。
还有一个间接的旁证。按这个说法，Da 越小，锥越细长，上半段「像倒过来的钟乳石」的部分应该越长。论文拿两根锥形石笋去拟合，Da 等于 0.23 的那根烛台，像钟乳石的那一段明显比 Da 等于 0.59 的那根长。
图3 两根锥形石笋的拟合。红线是本文的解析解，白虚线是把钟乳石的理论倒过来用的结果。尖端一段两者贴得很紧，越往下越分家：钟乳石那条算不出石笋底部那种突然变陡的走势。左图那根烛台 Da 等于 0.23，右图那根 Da 等于 0.59。引自原文 Figure 5。
六、拿三根真石头验一验
解析解再漂亮，也得跟真石头对一遍。
作者从 Postojna 洞拿了三根小石笋。这三根原本是别人采去做古气候研究、已经离开洞穴的标本。给它们做 X 射线断层扫描，把里面的三维形状重建出来，再把解析解画出的曲面直接叠上去。
A 和 B 是柱状，只有一个可调参数，就是半径。C 是平顶，两个：底半径和平顶半径。除此之外没有任何自由度。
吻合得有点过分。要知道解析解假设的是滴水速率恒定、水流是层流、形状完全轴对称，这几条在真洞穴里没有一条严格成立。
这套东西还能倒过来用。既然形状由 Da 和半径决定，那反过来，量出形状，就能反推滴水量。三根石笋算下来，流量分别是每秒 0.002、0.0045、0.001 立方厘米。按一滴水 0.1 毫升算，就是 50 秒一滴、22 秒一滴、100 秒一滴。
这个数以前只能靠人在洞里掐表数，现在可以从石头本身读出来。
七、最要紧的一条：形状会污染同位素信号
石笋之所以金贵，是因为它能读古气候。读的办法是测氧和碳的同位素，也就是 δ18O 和 δ13C。
但这两个数不是纯净的气候信号。水滴在石笋表面流动、二氧化碳往外跑的过程中，轻重同位素跑的速度不一样，会额外产生一个偏移，叫动力分馏。要还原气候，就得先把这个偏移扣掉。
以前的数值模型算过这个偏移，但只算柱状石笋。默认的结论是：偏移沿半径是一条抛物线，中心最小，边上最大。
这篇说，不对。形状不一样，偏移的分布也不一样，而且差在最关键的地方，中心线。
平顶石笋的中心那一圈，水的浓度是恒定的，所以 δ13C 在那里也是平的，一条水平线，一直平到湿区的边缘，出了那个圈才开始往上抬。柱状和锥状才是抛物线。
图4 同位素。横轴是离石笋中心的距离，纵轴是 δ13C。黑点是塞尔维亚 Ceremosnia 洞一根平顶石笋的实测值；红线是平顶的预测，中心一段是平的，到 9 厘米左右才抬起来；蓝线是柱状和锥状的预测，是一条抛物线。引自原文 Figure 6。
图4 里那根石笋的实测点，几乎全贴在红线上。中心那 9 厘米，确实是一条水平线。如果按老办法当柱状处理，用蓝线去读，中心会差出 1.4‰。
1.4‰ 在古气候里是什么概念？够写一篇文章了。
所以以后的规矩得改一条：读同位素之前，先看一眼这根石笋是什么形状。如果它是平顶，中心那一段不能按抛物线反演。反过来，在不同高度切几个径向剖面，还能看出形状什么时候变过。
八、这套理论管不到的地方
作者自己列了几条，都挺实在。
第一，滴得特别慢的时候，水不成连续的水膜。它是一滴一滴打在顶上，散成几道细流往下走，走几厘米就饱和了，不再沉淀。这意味着石笋的直径有一个下限，不可能无限细。这个「最小石笋」的问题几十年前就有人提过，连续膜模型天生答不了。
第二，真实洞穴的条件一直在漂。季节在变，滴水量在变，Da 就在变。所以一根石笋的轮廓记的不是当下的 Da，是它这一辈子 Da 走过的整条路。
第三，轴对称是理想化。洞里有风，水滴也不总落在正中间。
我自己的疑问是另一个：这套解假定石笋一直往上长、轮廓不变。可真实石笋会分叉，会长歪，会被水滴砸出一个坑。这些「不理想」的部分，恰恰是最好看的那些。
看完这篇，我回头又翻了一遍图1。
以前看这些照片，只觉得好看。现在知道，那些轮廓不是长出来的，是算出来的。一个数，管着它顶上那块是平的、圆的，还是尖的。
下次再进洞，大概会先抬头看看洞顶有多高，再低头看看石笋顶是什么样。
原文信息
Szymczak P., Ladd A.J.C., Lipar M., Pekarovič D. Shapes of ideal stalagmites. PNAS 2025, 122(42): e2513263122. DOI: 10.1073/pnas.2513263122.
数据说明
一，Damköhler 数定义为 Da = kA/Q，其中 k = 1.3×10⁻⁵ 厘米每秒是方解石沉淀速率常数，A 是石笋底面积，Q 是滴水的体积流量。柱状石笋 Da 等于 1，半径 R 等于根号下 Q 除以 πk。
二，平顶石笋的 Da 等于 1 加 k/P，P 是单位面积的滴水速率；平顶半径 Rc 等于 R 乘根号下 1 减 1/Da。湿区越小，Rc 越小，形状越接近柱状。
三，锥状石笋在 Da 很小时可近似为 z/R = C 减 Da⁻¹ 乘反双曲正切（r/R）。论文用完整解拟合两根锥形石笋，得到 Da 分别为 0.23 和 0.59。
四，三根断层扫描石笋反推的滴水流量：A 为 0.002 立方厘米每秒，B 为 0.0045，C 为 0.001；按一滴 0.1 立方厘米估算，滴间隔分别为 50 秒、22 秒、100 秒。
五，同位素预测采用论文的典型洞穴条件，平顶与柱状在中心处差约 1.4‰；实测样品来自塞尔维亚 Ceremosnia 洞，恒定区半径约 9 厘米。六，图片：图1 至图4 分别引自原文 Figure 1、3、5、6。图内英文由本文译为中文，数据与图形未作改动；面板字母、数学符号、单位及照片版权标记保留原样。