大型火山碎屑岩（ignimbrite）覆盖下埋藏地貌揭示喷发前岩石隆升速率

《SCIENCE ADVANCES》：Landscapes buried beneath large-volume ignimbrites reveal preeruptive uplift rates

　　

来自大规模爆炸式喷发的火山碎屑密度流（pyroclastic density current，PDC）可形成具有平面状、低角度表面的广阔火山碎屑岩（ignimbrite）沉积。研究人员证明了利用这些表面，可通过被埋藏古地貌的地形起伏来约束喷发前岩石隆升速率。研究人员使用由河流侵蚀驱动的景观演化模型，在给定表面坡度的火山碎屑岩之下探索了允许的景观几何形态范围。结果表明，河流河道陡度（channel steepness，ksn）与地形坡度作为岩石隆升速率的函数之间存在明显相关性。研究人员显示，要生成被火山碎屑岩沉积埋藏的平缓地貌，需要较低的岩石隆升速率；因此火山碎屑岩埋藏更倾向于记录造山带在构造-岩浆状态下的地貌，从而可能将长期地貌演化与大型爆炸式火山系统发育条件联系起来。将该模型应用于中央安第斯西翼早中新世Cardones火山碎屑岩，得到喷发前岩石隆升速率小于0.26千米每百万年（km Myr-1），与独立的热年代学（thermochronometric）和地温压力计（geothermobarometric）约束相当。

论文解读：《SCIENCE ADVANCES》刊发关于大型火山碎屑岩埋藏地貌与喷发前岩石隆升速率研究

研究背景方面，俯冲带汇聚板块边界分布着地球大型山脉与最具爆炸性火山。理解山脉演化需掌握岩石隆升历史，但长尺度岩石隆升速率难以获取。现有长期约束多依赖热年代学冷却史、地质压力计等间接代理，需数值模型转为地表参考速度，时间分辨率受数据封闭温度与空间分布限制；且缺乏独立检验，难在景观调整所需时标上评估隆升史。因此，研究人员提出一种地貌学方法：利用大型体积、空间广布火山碎屑岩之下埋藏地貌几何形态，约束喷发前岩石隆升速率，积分于景观响应时标。

主要关键技术方法方面，研究人员采用Landlab开源Python地表过程工具包构建二维景观演化模型，以河流动力侵蚀模型（stream-power incision model）描述河道下切；用最大高程剖面回归提取山脉尺度坡度（range-scale slope，S_R）；以χ图（χ-plot）与河道陡度（k_sn）量化河流地貌；并通过χ图随机化生成瞬态复杂河道剖面。研究区以中央安第斯西翼早中新世Cardones火山碎屑岩（>1260 km³，源自Lauca Caldera）为案例，其初始表面坡度约1.5°±0.3°。

研究结果部分，引言（INTRODUCTION）中，研究人员界定低长宽比火山碎屑岩（low aspect ratio ignimbrite，LARI）具0.5°至3°缓倾表面，并由源区向下倾方向减薄；提出完全埋藏条件为被埋景观山脉尺度坡度不大于火山碎屑岩初始表面坡度。研究人员以Cardones火山碎屑岩为安第斯案例，建立河流 Relief、河道陡度与山脉尺度坡度理论框架，公式给出地表抬升率 dz/dt = U(x,t) - K(x,t) A(x,t)^m S(x,t)ⁿ，稳态下 k_sn = (U/K)^1/n。

相对影响部分（Relative influence of erodibility and rock uplift rate on range-scale slope），560组模拟显示 S_R = k_sn/108.3；Cardones要求 S_R ≤ 1.5°，将U-K解限制于1.5°等值线之上。

现有可蚀性约束部分（Constraints from existing estimates of erodibility），归一化可蚀性 K 约1×10^-9 至1×10^-8 m^-1 yr^-1；结合中新世阿里卡（arid）条件，区域 K～1×10^-9 m^-1 yr^-1，得出 S_R≤1.5°时岩石隆升速率 <0.26 km Myr^-1。

景观响应时标部分（Constraining timescales of landscape response），大尺度100×50 km模拟表明低K、低U、高S_R景观达新稳态最慢；以95%稳态计，域长与初烈10 m时约14 Myr可达1.5°。

一维复杂剖面部分（Assessing complex topographies in 1D），χ图随机化显示复杂河道若位于稳态剖面之下，其平均 k_sn 必小于稳态；稳态代表最大陡度保守端元。

讨论部分（DISCUSSION）中，研究人员指出河流动力模型可刻画山脉尺度几何，m/n≈0.5时结论不依赖m、n绝对值；非稳态/非均匀地貌仅改变凹度，不推翻LARI坡度封顶逻辑。埋藏地貌低起伏类似现代挤压造山带低海拔段；Ignimbrite埋藏具地貌选择性，偏好保存缓慢或空间均匀隆升期形成的景观，将长期地貌演化与大型爆炸式火山系统发育的构造-岩浆状态联系。热年代学50 Ma至今 <0.2 km Myr^-1 与本文短中期 <0.26 km Myr^-1 一致，支持Incaic后构造较稳、Quenchua造山主要影响秘鲁而非北智利。

研究结论部分译为：利用中央安第斯西翼中新世Cardones低长宽比火山碎屑岩初始表面几何与河流动力模型，研究人员将渐新世—中新世平均岩石隆升速率上限约束为0.26 km Myr^-1；该结果提供与现有热年代学和地温压力计一致、稳健独立的中央安第斯发育约束，并与较慢俯冲带汇聚速率期相符。LARI作为长波长地貌标志，采样造山演化的特定阶段；区域广布火山碎屑岩埋藏优先记录缓慢或空间均匀隆升塑造的地貌，说明火山碎屑岩表面不是造山隆升的通用标志，而是筛选隆升史的地貌滤波器，并暗示长期地貌演化与大型爆炸式火山系统发育条件间存在联系。