公众号名称：地刊速览
作者名称：地刊速览
发布时间：2026-07-07 12:00
这篇论文用碎屑锆石定年和拉曼测温示踪新西兰南阿尔卑斯山五个流域的沉积物来源，发现泥沙并非主要来自冰川，而更符合整个流域均匀侵蚀的稳定状态。
冰川侵蚀被认为是山地最重要的侵蚀过程之一，但冰川到底贡献了多少沉积物，在学术上一直有争议。要回答这个问题，关键是要搞清楚河流里的泥沙到底是从哪儿来的——是从冰川覆盖的高海拔源头区，还是从山坡崩塌和河流侵蚀等其他地方。
新西兰的南阿尔卑斯山（Southern Alps）是研究这个问题的理想地点。这座山脉沿着阿尔卑斯断层快速隆升，从断层往东，岩石的冷却年龄和变质温度呈现出非常规律的梯度变化。研究人员利用这个特点，把锆石裂变径迹热年代学和碳质物拉曼光谱这两种方法结合起来，像给沉积物贴“示踪剂”一样，来追溯泥沙的来源。
他们在西侧选取了五个部分被冰川覆盖的流域，冰川面积占流域的12%到50%不等。在河床沙坝上采集了现代河流沉积物，然后分离出细砂粒级的锆石和含碳质物的岩屑颗粒进行分析。同时，他们设定了七种不同的物源模型——包括均匀侵蚀、坡度控制、滑坡物质、冰川侵蚀、冻融过程、水流路径和同震抬升等——再把这些模型预测的年龄和温度分布与实测数据进行统计比较。
结果出乎意料：最能解释实测数据的模型是“均匀采样”和“坡度依赖”模型，而不是“冰川侵蚀”模型。也就是说，当前这些河流输出的沉积物，并不是主要来自冰川覆盖的高海拔源头区，而是来自整个流域较均匀的贡献，且稍微偏向陡峭地形。冰川侵蚀模型在所有七个模型中拟合效果最差。不过，研究也发现了一个微弱的趋势：冰川覆盖面积越大的流域，冰川侵蚀的信号会稍微强一点，说明在冰期冰川扩张的时候，高海拔来源的物质占比可能会更高。
研究人员还发现，沉积物的主要来源区域大约在距离断层10到15公里的中游地带，这里恰好是冰川刚退去不久、山坡被削得最陡的地方。也就是说，冰川可能并不直接产生大量泥沙，而是通过削陡山坡，间接促进了后续的崩塌和侵蚀过程。更值得注意的是，“均匀采样”模型的良好表现提示，这些流域可能已经接近或达到了侵蚀平衡状态——即整个流域的侵蚀速率大致均匀。考虑到末次冰盛期时这里几乎完全被冰川覆盖，而冰退至今不过一万多年，这种快速的转型本身就说明冰川侵蚀和河流侵蚀之间的主导权切换，可以发生得比预想中快得多。
图1：研究区位置和采样流域
这张图展示了新西兰南岛（South Island）西侧五个采样流域的位置（星号标记），流域范围用色块标出。图中的黑圈和白圈分别代表已发表的锆石裂变径迹（ZFT）和碳质物拉曼光谱（RSCM）基岩数据点，用于拟合年龄和温度与断层距离的关系。灰色区域是末次冰盛期（MIS 2）的最大冰范围，蓝色是现代的冰川覆盖范围。
图2：河流纵剖面和基岩数据的拟合关系
图a显示了各流域河流的纵剖面，箭头标出了现代冰川的末端位置。图b是锆石裂变径迹冷却年龄随距离阿尔卑斯断层（Alpine Fault）的变化及拟合曲线，内插图放大了断层15公里以内的部分。图c是碳质物拉曼光谱得到的最高变质温度（Tmax）随断层距离的变化，两个流域组的拟合曲线分开显示。
图3：七种物源模型的采样概率分布
这张图以库克（Cook）流域为例，展示了七种不同物源模型中每个格网被采样的概率，颜色越亮表示概率越高。模型分别对应均匀侵蚀、坡度控制、滑坡物质、水流路径、冰川侵蚀、冻融过程和同震抬升等不同的地表过程。
图4：锆石裂变径迹的年龄分布
这张图展示了五个流域碎屑锆石裂变径迹冷却年龄的概率密度分布和核密度估计，每种颜色代表一个流域，并标出了通过混合模型自动识别的年龄组分。
图5：碳质物拉曼光谱的温度分布
这张图展示了五个流域碎屑碳质物的最高变质温度（Tmax）分布，所有样品都包含至少一个400–500°C的主组分。
图6：不同物源模型的对比结果
这张图把所有物源模型的预测分布（虚线）和实际观测分布（实线）放在一起对比，便于直观看出哪个模型跟实测数据最接近。
图7：KS检验统计量与冰川覆盖面积的关系
这张图展示了冰川侵蚀模型（模型5）的KS检验统计量随流域内冰川覆盖面积比例的变化，显示出两者之间微弱的负相关关系，说明冰川覆盖越大的流域，冰川侵蚀的信号可能稍微强一点。此外，山坡坡角的中位数也与冰川面积呈正相关。
