含水安山质熔体黏度模型：来自樱岛安山岩的约束与全球喷发前黏度重估

《Chemical Geology》：A viscosity model for hydrous andesitic melts: Insights from Sakurajima andesite

　　

岩浆黏度（viscosity）的研究对于理解火山过程至关重要，因为它控制着岩浆运移（magma transport）、变形响应（deformation response），并最终决定喷发方式（eruptive style）。这对安山质岩浆（andesitic magmas）尤为相关，其熔体黏度相对基性岩浆（mafic magmas）较高，且通常具有较高的喷发前含水量（pre-eruptive water contents），因而倾向于爆炸性喷发行为。熔体脆性（melt fragility, m），定义为Angell图中玻璃化转变温度（glass transition temperature, Tg）处黏度-温度关系的陡峭程度，是结合Tg在宽温度范围内描述黏度的关键参数。先前研究提出，m可通过纵波和剪切波声速经由体积-剪切模量比（bulk-to-shear modulus ratio, K/G）估算，或与每个四面体配位阳离子的非桥氧数（NBO/T）相关联。然而，近期研究表明，基于这些方法以及广泛使用的黏度模型得出的脆性估算，可能导致熔体黏度预测不准确，在无水和含水条件下对安山质、玄武质、碱玄质响岩（tephri-phonolitic）和响岩质碱玄岩（phono-tephritic）成分的计算值与实验测量值存在系统性偏差。这种预测失败源于常影响黏度测量的过程，包括熔体的非晶纳米结构重组、Fe-Ti氧化物纳米晶（nanolites）的结晶以及实验过程中的水损失。关键的是，纳米晶不仅限于实验室，其存在已在多种成分的天然火山产物中得到记录，包括2011年新燃岳（Shinmoedake）喷发的安山质产物。这些观察表明，纳米晶可在岩浆上升和喷发过程中形成，产生耦合的流变学和化学变化。然而，含水岩浆中这些流变学变化的幅度仍缺乏约束，代表了喷发建模和预测中的关键空白。研究人员通过开发一个含水安山质熔体黏度模型来解决这一空白，该模型考虑了温度、水含量和过渡金属浓度（如Fe、Ti和Mn）的影响。研究人员使用了代表樱岛（Sakurajima）火山的含过渡金属和不含过渡金属的模拟成分，樱岛火山是日本最活跃和最具危险性的火山系统之一，其喷发产物以富含纳米晶而闻名。研究人员结合了差示扫描量热法（differential scanning calorimetry）、黏度测定法（viscometry）、拉曼光谱（Raman spectroscopy）和布里渊光谱（Brillouin spectroscopy）来确定玻璃化转变和熔体脆性，而不会引发损害现有数据集的实验假象。基于已有方法，研究人员推导出将熔体脆性与布里渊光谱衍生的声速相关联的简化经验方程，分别针对含过渡金属和不含过渡金属成分进行校准，从而无需耗时的高温测量即可快速确定熔体脆性。由此产生的黏度模型为樱岛火山及一系列安山质成分提供了修正的喷发前黏度估算，对安山质火山岩浆上升和破碎模型具有潜在意义。

论文解读：含水安山质熔体黏度模型——来自樱岛安山岩的约束与全球喷发前黏度重估

一、研究背景与问题提出

岩浆黏度是控制火山过程中岩浆运移、变形响应乃至最终喷发方式的核心物理参数。对于安山质岩浆而言，其熔体黏度显著高于基性岩浆，加之喷发前通常富含水，使其倾向于爆炸性喷发，对环太平洋火山带等人口密集区域构成严重火山灾害威胁。因此，精确的黏度模型对于火山灾害评估和岩浆动力学数值模拟至关重要。

熔体脆性（melt fragility, m）作为描述黏度-温度关系陡峭程度的关键参数，结合玻璃化转变温度（glass transition temperature, T_g）可在宽温度范围内预测黏度。传统上，m的估算依赖于体积-剪切模量比（K/G）或非桥氧数（NBO/T）等指标。然而，近期研究揭示这些方法以及广泛使用的黏度模型在安山质、玄武质等成分上存在系统性偏差，计算值与实验测量值不符。

这种预测失败的根本原因在于黏度测量过程中常见的一系列复杂过程：熔体的非晶纳米结构重组、Fe-Ti氧化物纳米晶（nanolites）的结晶析出以及实验过程中的水损失。尤为关键的是，纳米晶不仅存在于实验室，在2011年新燃岳（Shinmoedake）喷发等天然火山产物中也有广泛记录，表明其可在岩浆上升和喷发过程中形成，从而产生耦合的流变学和化学变化。然而，含水岩浆中这些流变学变化的幅度长期缺乏定量约束，成为喷发建模和预测中的关键空白。

二、研究内容与意义

针对上述空白，研究人员以日本最活跃和最具危险性的火山系统之一——樱岛（Sakurajima）火山为对象，开发了一个含水安山质熔体黏度模型。该模型系统考虑了温度、水含量和过渡金属浓度（如Fe、Ti和Mn）的综合影响。研究采用含过渡金属（AND100）和不含过渡金属（AND0）的模拟成分，以及过渡金属含量为65%的中间成分（AND65），分别代表樱岛安山岩的天然化学组成及其简化类似物。

通过将黏度测定与光谱学方法相结合，研究人员成功规避了传统实验方法中易引入的假象，建立了成分依赖的K/G-脆性关系，并据此提出了修正的喷发前黏度估算。该研究的意义在于：首次系统量化了纳米晶化过程中Fe-Ti提取对含水岩浆流变学的影响，为安山质火山岩浆上升和破碎模型提供了关键流变学约束，并揭示了全球范围内安山质喷发前熔体黏度可能被系统性高估的现象。

三、关键技术方法

研究人员综合运用了多种实验与表征技术。样品合成方面，AND100含水玻璃采用内加热压力容器（IHPV）装置合成，AND65和AND0含水玻璃采用活塞-圆筒装置合成，均在德国拜罗伊特大学巴伐利亚地球科学研究所（BGI）完成。成分分析采用电子探针（EPMA）和扫描电子显微镜（SEM）；铁氧化态通过穆斯堡尔光谱（M?ssbauer spectroscopy）表征；水含量通过卡尔费休滴定（KFT）和傅里叶变换红外光谱（FTIR）双重测定。流变学测量方面，低黏度区采用微压入（MP）黏度计和同心圆柱（CC）黏度计，玻璃化转变区采用常规差示扫描量热法（C-DSC）和闪速差示扫描量热法（F-DSC）。此外，布里渊光谱（Brillouin spectroscopy）用于获取K/G比值，拉曼光谱用于监测实验前后样品纳米-微观尺度的结构变化。

四、研究结果

4.1 成分与结构表征

EPMA分析确认了各合成玻璃的化学组成与设计目标一致，AND100系列含有完整的过渡金属（FeO_tot约6.77 wt%，TiO₂约0.79 wt%），AND0系列几乎不含过渡金属（FeO_tot约0.04 wt%，TiO₂约0.04 wt%）。穆斯堡尔光谱显示，水饱和条件下合成的AND100样品具有较高的Fe³⁺/Fe_T比值（0.64），而含水样品的比值较低（0.22-0.42），反映了不同合成条件下的氧逸度差异。SEM背散射成像确认了起始玻璃中不存在微晶和气泡。

4.2 玻璃化转变温度与脆性测定

C-DSC和F-DSC测量揭示了水含量对T_g的显著影响：AND100系列从无水状态的654.0°C降至含水5.15 wt%时的380.3°C；AND0系列从737.0°C降至含水5.56 wt%时的383.4°C。基于K/G比值与脆性的经验关系，计算得到的脆性值m在AND100系列中为30.5-34.7，在AND0系列中为31.8-35.1，均随水含量增加而增大，表明水使熔体变得更“脆弱”（fragile）。

4.3 黏度模型建立

研究人员建立了成分依赖的K/G-脆性关系，分别针对含过渡金属和不含过渡金属成分进行校准。该关系使得仅通过布里渊光谱测量声速即可快速确定熔体脆性，无需耗时的高温黏度测量。新模型综合考虑了温度、水含量和过渡金属浓度的协同效应，能够更准确地预测含水安山质熔体在宽温度范围内的黏度行为。

4.4 喷发前黏度重估

应用新模型对樱岛火山喷发前条件进行估算，结果显示其喷发前熔体黏度比先前估计值低约9倍。更广泛地，研究人员发现全球安山质喷发前熔体黏度可能被系统性高估，最大偏差可达27倍。这一发现对岩浆上升速率、脱气行为以及喷发风格的数值模拟具有重要影响。

五、讨论与结论

讨论部分指出，传统黏度模型对安山质成分预测失效的根源在于未能充分考虑纳米晶化过程中Fe-Ti氧化物的提取效应以及实验过程中的水损失。纳米晶的形成不仅改变了残余熔体的化学组成，还通过耦合的流变学和化学变化影响岩浆的上升和喷发动力学。研究人员强调，新的K/G-脆性关系为快速评估熔体脆性提供了便捷途径，而新黏度模型则为安山质火山喷发前黏度的精确估算提供了可靠工具。

研究结论表明，含水安山质熔体的黏度受温度、水含量和过渡金属浓度的共同控制，纳米晶化过程中的Fe-Ti提取效应不可忽视。樱岛火山喷发前熔体黏度约为先前估计的九分之一，全球安山质喷发前黏度可能被高估多达27倍。这些修正的黏度估算为安山质火山岩浆上升和破碎模型提供了新的流变学约束，有助于更准确地评估火山喷发风险和灾害潜力。