矿物岩石地球化学通报  2017, Vol. 36 Issue (1): 14-25   PDF    
硅酸盐岩浆液态不混溶作用的理论基础概述
侯通1,2     
1. 中国地质大学(北京), 地质过程与矿产资源国家重点实验室, 北京 100083;
2. Institute of Mineralogy, Leibniz University of Hannover 30167, Germany
摘要: 液态不混溶作用是岩浆演化过程中的重要作用之一。本文针对具有普遍意义的硅酸盐岩浆液态不混溶作用的研究历史和进展,对其理论基础,包括熔体结构、相平衡、动力学和元素分配以及同位素分馏等问题进行了阐述。同时还就岩浆混溶作用中涉及到的稳定态和亚稳态液态不混溶不同的相平衡关系,以及从热力学角度液态不混溶作用发生的驱动力,不混溶相从成核到长大再到最终相分离的动力学过程进行了详细介绍。此外,针对岩浆系统,还总结了发生不混溶的条件和不混溶作用过程中元素的地球化学行为。最后以攀枝花层状岩体中部岩相带的形成为例,说明粒间熔体液态不混溶作用在韵律层形成过程中发挥的重要作用。
关键词: 硅酸盐      熔体结构      液态不混溶      岩浆     
A Brief Introduction on the Principle of Silicate Liquid Immiscibility
HOU Tong1,2     
1. State Key Laboratory of Geological Process and Mineral Resources, China University of Geosciences, Beijing 100083, China;
2. Institute of Mineralogy, Leibniz University of Hannover 30167, Germany
Abstract: Silicate liquid immiscibility is one of the most important magma processes. This paper reviews theoretical principles for the separation of conjugate silicate melts, i.e. silica-rich vs. iron-rich, including the role of melt structure, kinetics, and the partitioning of elements and fractionation of stable isotopes between the two immiscible conjugates. The Coulombic expulsion between ions is the root reason that facilitates the liquid immiscibility at the atomic level. On basis of the recent progress of experiments, this paper presents phase equilibria relevant to stable and metastable immiscibility, and the driving force for immiscibility both thermodynamics and kinetics studies on the nucleation barrier in producing unmixing, coarsening and potential separation of equilibrated melts. Moreover, this paper reviews the geochemical behavior of elements and isotopes during the liquid immiscibility. At the end of the paper, taking the Middle zone of the Panzhihua layered intrusion as an example, this paper proposed a model advocating the interstitial liquid immiscibility probably played an important role in the formation of rhythmic layers.
Key words: silicate     melt structure     liquid immiscibility     magma    

硅酸盐熔体液态不混溶作用(silicate liquid immiscibility)是岩浆演化的重要方式之一。但是, 现代岩浆岩岩石学建立至今, 学术界对于这一概念从认识、理解到接受可谓一波三折——自20世纪初的新兴(Daly, 1914), 到20世纪20年代遭到了以Bowen为代表的科学家们的质疑(Greig, 1927Bowen, 1928), 直至上世纪后半叶起又逐步“复兴”(De, 1974McBirney and Nakamura, 1974Dixon and Rutherford, 1979Philpotts, 1979, 1982Roedder, 1979)。这不仅仅是因为早期测试手段落后且缺少必要的实验验证, 更多情况下是由于岩石学家对于自然样品所传递信息不同程度的忽视所导致的。例如, 月球玄武岩中的不混溶现象就是Roedder(1978)在橄榄石和斜长石中寻找流体包裹体的过程中偶然发现的。笔者认为, 从某种意义上说, 如果把岩浆液态不混溶作用这一理论比作岩浆岩岩石学中的“灰姑娘”, 它长期得不到国内外同行认可的原因就是因为一直没能“穿上”隐藏在自然样品形成过程中的不混溶作用这双“水晶鞋”。近年来, 围绕着岩浆液态不混溶作用的岩相学、地球化学和实验岩石学研究又掀起了新一波的热潮(Thompson et al., 2007Veksler et al., 2007Namur et al., 2012Charlier et al., 2013Kamenetsky et al., 2013Veksler and Charlier, 2015Hou and Veksler, 2015)。有理由相信, 在这一次研究热潮之后, 对于液态不混溶作用在岩浆演化过程中所起到的重要作用将会有更加清晰的认识, 再也不会被忽视。

一般来说, 岩浆液态不混溶作用可以导致: ① 硅酸盐-碳酸盐熔体的分离, 从而形成碳酸岩(Freestone and Hamilton, 1980Dasgupta et al., 2006);② Fe-S-C-Si系统的相分离, 尤其是富铁的金属相与硅酸盐熔体的分离, 对地球以及其他类地行星的形成有重要的启示意义(Corgne et al., 2008);③ 硫化物熔体和硅酸盐熔体的相分离, 可以形成铜镍(钴)硫化物矿床(Naldrett, 1999)等。值得指出的是, 不论是在实验中还是在自然界, 除上述3种液态不混溶作用以外, 更为普遍的是发生在2种硅酸盐熔体(例如富铁玄武质岩浆-富硅流纹质岩浆)之间的液态不混溶作用(Roedder and Weiblen, 1970Philpotts, 1979Veksler et al., 2007Jakobsen et al., 2005, 2011)。本文将主要讨论具有普遍意义的硅酸盐岩浆液态不混溶作用的研究历史和进展, 回顾其理论基础, 包括熔体结构、相平衡、动力学(kinetics)和元素分配以及同位素分馏等, 并在文末以攀枝花层状岩体韵律层最为发育的中部岩相带为例, 揭示液态不混溶作用在韵律层形成过程中的所扮演的重要角色。

1 硅酸盐液态不混溶作用 1.1 原子尺度上的液态不混溶

一般来说, 如果液相之间化学键(离子键、共价键和金属键以及分子间作用力)不同, 则彼此无法混溶。对于硅酸盐熔体, 它基本上是不同组分以共价键和离子键相互作用的产物。因此, 当硅酸盐熔体中出现了以金属键为主(如金属和金属硫化物熔体等)、解离的离子键(如熔融的盐类)或者分子态(如超临界的C-O-H流体)的液相, 则会通过不混溶作用与之分离。不过, 具体到本文讨论的硅酸盐熔体, 两相由同样的化学键组成, 不混溶作用成因则更加复杂, 下面将从熔体结构的角度来探讨不混溶作用发生的本质。

根据对熔体聚合程度(degree of polymerization)的影响, 阳离子可以分为2大类: 网络形成子(network former)和网络修饰子(network modifier)。网络形成子一般离子半径小而价态较高, 因此在熔体结构中, 配位数一般都不大于4。例如四次配位的Si4+是岩浆中最常见的网络形成子, 其他的还有P5+, B3+和某些微量元素(Ge4+、Sb3+和As3+)。网络修饰子一般离子半径要大于87.2 pm(1 pm=0.001 nm), 配位数一般不小于5(Hudon and Baker, 2002), 包括碱金属、碱土金属(但不包括Li+, Be2+, Mg2+)、稀土元素、Th4+和U4+等。离子半径大于Si4+(26 pm)但小于87.2 pm的离子, 属于两性行为离子, 其配位方式既可以是4次配位, 也可以是5次或6次配位, 在岩浆中主要的是Al3+, 其他的还包括Li+、Be2+、Mg2+、Mn2+、Zn2+、Sc3+、Ga3+、Ti4+、Zr4+、Hf4+、Nb5+、Ta5+和某些高场强元素。

SiO2或者富SiO2的熔体黏度很高, 如果向其中掺杂一定量的金属氧化物(如网络修饰子)基本都会出现不混溶现象。金属氧化物-SiO2二元系中, 不混溶域的大小和金属阳离子的电价和离子半径的比值(Z/r, 亦称离子势)密切相关(Hess, 1995)。如果金属元素为碱金属或者碱土金属, 这种效应会更加明显(图 1)。因此, 研究者普遍认为在硅酸盐熔体中发生液态不混溶作用的根本原因是网络修饰子之间的库仑排斥力(Coulombic repulsion;Hess, 1995, 1996Thompson et al., 2007)。不过也有例外, 当网络修饰子在硅酸盐熔体结构中同时被桥氧(bridging-oxygen;Si-O-Si)和非桥氧(non-bridging oxygen-NBO;Si-O-M)包围时, 这2类氧会“遮蔽”网络修饰子所携带的正电荷。但是, 当桥氧围绕网络修饰子时, 后者会与Si4+形成更强的共价键, 网络修饰子的正电荷则无法被有效“遮蔽”, 网络修饰子之间的库仑排斥力进而变得更强, 并最终导致不混溶作用的发生。实验表明离子势越高, 这种相分离的发生越容易。对于双性行为离子, 当其为4次配位时, 库仑排斥力则表现不明显, 不会发生液态不混溶作用。

二元系相图显示出不混溶域(据Levin et al., 1964) 图 1 碱金属、碱土金属氧化物与二氧化硅的 binaries with two liquid fields(after Levin et al., 1964) Figure 1 Silica-rich portions of alkali and alkaline-earth oxides

具体到硅酸盐岩浆液态不混溶作用, 多数都发生在岩浆富铁的情况下, 比如富铁玄武岩和流纹岩。因此, 这种液态不混溶作用很可能与铁离子在熔体结构中的特殊性有关。元素周期表中过渡元素列表中第一行除了铁以外的所有元素(Co2+、Ni2+、Cu2+、V3+、Cr3+), 它们的核外电子排布都比较类似, 在d亚层有5个轨道且均未充满电子, 所以对原子核的遮蔽效果都比较差。因此, 这些过渡元素在金属氧化物-SiO2二元系中的不混溶域要比根据离子势推测的区域要大, 如果是三价过渡元素(Cr3+、V3+和Fe3+)则不混溶域扩大的效应更加明显(Hudon and Baker, 2002)。

1.2 部分混溶双液系统

在岩浆演化过程中, 部分混溶双液系统的特征是随着温度降低, 单一液相将自发地变成双液相。主要表现为2种液相的溶解度随温度的升高而增加, 以致达到某一温度(最高会溶温度)时, 2种液相组成相同, 形成了单一的液相。继续升温时, 无论成分如何变化, 仅有单相区的存在。从热力学的角度来说, 如果不混溶作用是系统自发过程, 则在这一过程中必然伴随系统吉布斯自由能的降低(James, 1975Doremus, 1994Shelby, 2005)。能量-组成变化图(图 2a)显示, 系统组成和能量构成一光滑曲线, 和完全互溶的系统相比, 不混溶系统在曲线上出现一个局部的能量最大值和2个能量低值点, 这2个点就是不混溶两相a和b, 组成分别为L1和L2。以这两点为切点的直线与纵轴的交点就是不混溶两相各自的化学势(μaμb)。因此, 在a和b点之间的曲线上任意一点的组成在热力学上都是亚稳态的, 它们会倾向于发生液态不混溶, 从而使体系能量降低。此外, 由于L1和L2的能量和温度密切相关, 可将能量-组成曲线转换为温度-组成曲线(图 2b), 图中a和b两点之间的曲线被称为双节线(binodal)。一般来说, 这个曲线都是上凸(开口向下)的, 这样的体系又称为最高会溶温度(consolute temperature)部分互溶系统。一个处于a和b之间的单一组成随着温度的降低, 在逐渐接近最高会溶温度的过程中会发生临界现象, 又被称为临界乳光(critical opalescence;图 2b)。物理学的研究表明, 在临界点附近有许多有趣的现象, 样品的许多物理特性, 例如: 比热、压缩率等都有一些奇特的发散现象(Simmons et al., 1970Simmons and Macedo, 1971Mazurin and Porai-Koshits, 1984)。此外, 图 2a的曲线ab之间有2个拐点s1和s2, 其二阶导数(δ2GX2)T, P为0, 在温度-组成图(图 2b)中, 对应2点随着温度的变化曲线即为旋节线(spinodal)。图 2a中s1和s2之间的曲线上每一点的二阶导数(δ2GX2)T, P均大于0, 因此有局部最大值, 说明系统不稳定, 即任意一点的组成会自发发生液态不混溶来降低系统总能量, 这个能量差值就是发生不混溶作用的驱动力。对于曲线a-s1和s2-b来说, 曲线上每个点的二阶导数(δ2GX2)T, P均小于0, 因此有局部最小值, 系统相对稳定, 对于这2段曲线上的点, 如果要发生不混溶作用到达稳定的a和b, 则需要能量的涨落或者扰动(fluctuation)来克服势垒(详见下文)。

(a)吉布斯自由能-组成图, T1高于最高会溶温度, Tc代表最高会熔点温度, T2低于最高会熔点。切线与曲线相切的2个点a和b的组成分别为相互平衡的不混溶两相L1和L2μaμb分别为L1和L2的化学势。2个拐点s1和s2之间的区域为旋节线。(b)双节线和旋节线在组成-温度坐标系中的示意图 图 2 双节线和旋节线示意图(修改自James, 1975) Figure 2 Binodal and spinodal in a schematic binary A-B (modified after James, 1975)
1.3 稳定态和亚稳定态不混溶

对均一液相来说, 降温除了可能会导致不混溶以外, 更为常见的是导致出现结晶。以玄武质岩浆为例, 系统开始结晶矿物的温度称为液相线温度(liquidus temperature), 在岩浆液相线温度以上出现的不混溶都是稳定的, 而且在双节线和液相线相交时可以出现不混溶两相与矿物晶体的共存(图 3a)。此外, 在液相线以下岩浆处于过冷状态, 也会出现不混溶, 而且随着过冷度的增加, 不混溶的效应愈加明显(图 3b)。但是这种不混溶作用是亚稳态的, 会使系统中不混溶域扩大, 这种现象在多组分的玄武质岩浆系统(图 4)中也比较常见。

(a)液相线以上的稳定不混溶域和液相线下的亚稳态不混溶域(虚线区域);(b)液相线的亚稳态不混溶域(修改自Thompson et al., 2007) 图 3 温度和组成坐标系 Figure 3 Schematic binary temperature vs. composition phase diagrams(modified after Thompson et al., 2007)
1.4 形核能量势垒(nucleation barrier)

现代动力学研究表明, 液态不混溶作用发生相分离的初期还是以形核和生长为主要的作用机制(James, 1975Mazurin and Porai-Koshits, 1984Shelby, 2005及其引文), 因此, 要产生不混溶的新相则需要从旧的相中成核。那么新相如何成核以及成核速率的快慢则必然要同时受到系统热力学和动力学条件的制约。目前, 有关不混溶作用发生时的相成核过程还很难定量理解, 只能依据现有的理论进行定性描述。如图 5(图 2a已做介绍)所示, 如果有一个组成为x的熔体, 成分介于a(双节线)和s1(旋节线)之间, 系统会倾向于通过不混溶作用来形成目标产物——a和b相来达到稳定态。在曲线上, x的切线与b点的自由能差值ΔG就是此时发生相分离的驱动力, 会使其中一个不混溶相的组成向b点移动。不过, 从切点向b点移动的过程中, 是有能量的涨落的, 因此在成核过程中会形成一些成分上类似c的新相。如果要达到b点的组成, 则要克服这个能量障(能量势垒)Δg(图 5), 至少到达d点的组成, 系统才会继续自发生成不混溶的b相。因此, 在b相生成之前出现的不混溶现象在成分上都是不稳定的, 如果有能量的扰动成功克服形核势垒, f相就会在组成上变为b相。此外, 如果初始成分处在s1和s2之间(二阶导数大于0;旋节线以内), 在发生不混溶的过程中则不需要克服形核能量势垒。不过即使如此, 在相分离过程中, 动力学的制约作用也是普遍存在的, 因为除了能量势垒, 要想完成相分离还会受到不同元素扩散率和熔体黏度的影响。

Qz-石英;Fa-铁橄榄石;Or-正长石;Hd-钙铁辉石;An-钙长石 图 4 (a) Qz-Fa-Or伪三元系;(b)Or-Qz-Fa-Hd相图;(c)Or-Qz-Fa-An相图(显示稳定不混溶域和亚稳态不混溶域) (修改自Roedder, 1951, 1978Visser and van Groos, 1979Hoover and Irvine, 1977Bogaerts and Schmidt, 2006Irvine, 1976) Figure 4 Liquid immiscibility in the(a)fayalite Fa-orthoclase Or-silica Qz pseudo-ternary join;(b)join Fa-Hd-Or-Qz);(c)join Fa-An-Or-Qz(Solid grey curve shows the extent of sub-liquidus immiscibility)(modified after Roedder, 1951, 1978Visser and van Groos, 1979Hoover and Irvine, 1977Bogaerts and Schmidt, 2006Irvine, 1976)

图 5 低于最高会溶温度自由能-组成变化图解显示导致不混溶两相分离的热力学驱动力(据James, 1975) Figure 5 Schematic free energy-composition diagram for TTc showing graphical method of determining the thermodynamic driving force for liquid-phase separation(after James, 1975)
2 岩浆液态不混溶作用的实验研究进展 2.1 岩浆演化路径和不混溶作用发生

除了由陨石撞击产生的高温岩浆外, 常见的由地幔(包括地幔柱)或者地壳部分熔融产生的岩浆在离开源区后很快就会开始结晶矿物。因此, 控制残余岩浆成分变化, 即决定岩浆演化路径(liquid line of descent, 下降液体线)的关键因素是发生分离结晶的矿物的种类和含量。一般来说, 对于名义上无水的玄武岩浆体系, 发生液态不混溶作用的先决条件是残余岩浆富铁(FeO>18%~19%, wt.)且没有发生铁氧化物的分离结晶。不过Charlier和Grove(2012)的实验也表明即使有磁铁矿的结晶分离, 即使岩浆铁含量已经开始有所下降, 只要主要岩浆演化路径还在不混溶域内, 不混溶作用就仍然可以发生。因此, 在岩浆演化过程中, 能否发生液态不混溶的关键是岩浆演化路径能否穿过不混溶域(图 6)。例如, 以拉斑玄武岩为初始成分的2个经典实验(Juster et al., 1989Toplis and Carroll, 1995), 其演化路径都没有穿过这一区域, 因此没有在实验中观察到不混溶现象。此外, 如果玄武质岩浆由于富水(水活度aH2O=0.1~1;Botcharnikov et al., 2008)而导致氧逸度升高, 在演化过程中, 磁铁矿会成为早期液相线矿物, 这样岩浆演化路径较快地偏向富硅一方也无法穿过不混溶域。

(a)显示实验(Charlier and Grove, 2012)确定的不混溶域。玄武岩体系实验获得的下降液体线, [据Toplis and Carroll,1995(T & C; 黑色);Juster et al. 1989(J & G;灰色)];(b)显示液态不混溶作用的开始和结束都和液相的演化趋势直接相关。NBO-非桥氧数量, T=Si+Al+P+Ti(据Charlier and Grove, 2012) 图 6 温度-熔体聚合程度变化图解(a)和下降液体线穿过不混溶域示意图(b) Figure 6 Temperature as a function of the degree of polymerization(NBO/T)of experimental melts(a), and Schematic illustration of the evolution of the liquid line of descent through a two-liquid field(b)

在三元系的相图中, 不混溶域既可以在同结线(cotectic)上, 也可以在“热堤”(或热坝;thermal barrier)上(图 7)。这两种情况会导致液态不混溶作用发生后共轭两相(immiscible conjugate)演化路径出现明显的差别。换句话说, 当液态不混溶作用发生后直到结束, 岩浆演化路径会出现两种情况: 例如, 在Roedder系统(图 4a)中, 不混溶域覆盖部分同结线, 不混溶作用并不会改变岩浆的演化路径, 液相成分可以从Fa-Q同结线开始演化, 液态不混溶作用之后, 虽然两相各自演化, 但是都在生成钾长石的共结点(eutectic point)结束, 可谓“殊途同归”。不过在某些三元系(如CaO-B2O3-SiO2)中, 不混溶域则在“热堤”上。如果在岩浆演化过程中出现这样的情况, 就会导致发生不混溶作用之后的两相彼此之间由于“热堤”的阻拦无法互相演化, 从而各自沿着不同的路径演化。对于玄武岩体系, 尤其是富铁玄武岩-安山岩体系, 在低压环境下演化时, 斜长石一般都为早期的液相线矿物。因此, 斜长石在玄武岩体系中稳定存在时往往会形成“热堤”(例如低压Ol-Cpx-Pl热堤;Yoder and Tilley, 1962), 这样由不混溶作用形成的富铁和富硅会各自演化, 互不干扰, 在富铁和富硅各自的共结点结束演化。不难看出, 这2种演化路径分别与Fenner和Bowen演化趋势非常类似。

a-不混溶域在同结线上;b-不混溶域在“热堤”上 图 7 液相面和不混溶域之间的2种关系(据Veksler and Charlier, 2015) Figure 7 Two topological types of liquid immiscibility regions in relationship to liquidus surface(after Veksler and Charlier, 2015)

如果岩浆演化路径只要穿过不混溶域就会出现液态不混溶现象的话, 那么在自然界应该有很多的实例, 但通常在自然样品的研究过程中又很难从岩相学上识别出液态不混溶现象。主要原因是, 如果不混溶作用在演化过程中发生并在演化结束前消失的话, 那么它只是发挥“隐形”作用, 在演化的最终产物中, 是无法观察到任何不混溶作用的直接岩相学证据的。具体来说, 岩浆液相成分进入不混溶域后, 随着岩浆温度降低和两相中矿物结晶作用的进行, 残余熔体还会离开不混溶域(图 6b), 从而使不混溶两相消失。以伪三元系Qz-Fa-Or为例(图 8), 当液相成分在X点时(铁橄榄石首晶区), 此时铁橄榄石开始不断结晶使残余液相成分进入不混溶域, 然后形成不混溶两液相——富铁的L1和富硅的L2。不过富铁相还会继续和铁橄榄石反应, 使两相成分点沿着降温矢分别向着L1′和L2′点移动, 而使两相成分彼此接近。直到富硅相的成分点和铁橄榄石的连线穿过初始组分X点时, 所有的富铁相被消耗殆尽, 而此时系统只有单一的富硅相。因此, 当系统完全固结之后, 笔者未观察到任何不混溶作用发生的证据, 但这并不代表液态不混溶作用在岩石形成过程中没有发挥作用。

X点为起始液相成分, 随着温度的下降, 液态不混溶作用发生并随后消失。灰色区域为稳定不混溶域。Qz-石英, Or-正长石, Fa-铁橄榄石 图 8 Fa-Or-Qz伪三元系相图(据Roedder 1951, 1978) Figure 8 Illustration of a crystallization path for liquid composition X in the Fa-Or-Qz pseudo-ternary join (after Roedder, 1951, 1978)
2.2 不混溶两相之间的元素分配和同位素分馏

前文已介绍液态不混溶作用在原子尺度上和不同离子之间的相互作用力(库仑定律)有关, 而且不混溶两相必须在热力学上是彼此平衡的。因此, 在平衡两相之间的元素分配和不同元素的离子势有关, 即元素分配要遵循能斯特分配定律, 分配系数Di=CLFe/CLSi, 其中C为元素质量, 下标为元素种类, 上标表示不混溶的两相(LFe为富铁相, LSi为富硅相, 下同)。研究表明分配系数的大小和指定元素-SiO2二元系曲线开口的大小以及最高会溶温度之间存在密切的联系(Hudon and Baker, 2002Schmidt et al., 2006Veksler et al., 2006): 双节线向下开口越大(图 1), 最高会溶温度越高(图 9a), 这种元素在富铁相中的分配系数越大。在离子势和分配系数关系的图解上, 只有微量元素Cs、Rb和Sb, 还有主量元素K, Al和Si的分配系数小于1, 即倾向于进入到富硅相中, 其他元素都倾向于进入富铁相。其中Cs和Rb离子势最低, 属于网络修饰子, Sb属于网络形成子(图 9b)。分配系数最大的包括稀土元素和三价过渡元素Cr和V。虽然分配系数的绝对值会随着温度、压力、氧逸度和液相成分的变化而变化, 但不同元素分配系数之间的相对大小不变。

使用激光剥蚀电感耦合等离子体质谱仪测定实验样品(LFe=富铁液相;LSi=富硅液相)中元素含量 图 9 不同元素与SiO2二元系最高会溶温度和离子势的关系(a)(据Hudon and Baker 2002) 和不同元素在不混溶两相之间的能斯特分配系数和离子势关系(b)(Veksler私人通讯) Figure 9 Consolute(critical)temperatures of the miscibility fields in the metal oxide-silica binaries(a)(after Hudon and Baker 2002), and Liquid-liquid Nernst partition coefficients D plotted against ionic potential Z/r(b)

Kyser等(1998)在简单体系K2O-FeO-Fe2O3-Al2O3-SiO2的实验(1180℃, 105 Pa)中发现在不混溶的两相中, 富硅相的 δ18 O值相对富铁相要高, 出现了氧同位素的分馏。在类似的体系中, 加入了挥发分(水、磷和硫)之后(Lester et al., 2013b), 富硅相的氧同位素值依然高于富铁相, 总体上大约高0.4‰~1‰。铁和镁同位素在不混溶两相之间的分馏效应与氧同位素类似, 都在富硅相中富集重同位素。这说明稳定同位素在不混溶两相之间的分馏很可能和熔体的聚合程度有关, 聚合程度越高(富硅相), 越富集重同位素(Zhu et al., 2015)。

2.3 岩浆不混溶作用发生的动力学过程

实验表明, 在降温过程中, 当熔体成分演化到双节线(双节线以上)附近时淬火会出现成分不均一(微米级的富铁和富硅相共存)的玻璃(图 10a)。而这种现象在自然界也普遍存在, 只是很多时候被研究者直接忽视。例如, 快速冷凝的玄武岩中可见在斑晶斜长石边缘出现有2种成分构成不均一的玻璃, 相对于基质(富铁玄武岩), 一相更富铁, 一相更富硅(图 10b)。如何在动力学上合理的解释这种现象是区分稳定不混溶和亚稳态不混溶的关键。

Mt-磁铁矿; Fa-铁橄榄石; Pl-斜长石 图 10 实验样品与自然样品中成分不均一的玻璃 Figure 10 Heterogenized glass in both experimental and natural samples

在过去几十年, 围绕着岩浆不混溶作用发生的动力学过程这一主题, 科学家们完成了大量的实验(Visser and van Groos 1976, 1977Freestone and Hamilton 1977Roedder 1978Veksler et al., 2007, 2008a, 2008b), 但争议依旧。一派观点是认为这种微米级的不均一是在缓慢淬火过程中出现的亚稳态不混溶, 也就是说在淬火前熔体成分是均一的, 只是由于淬火缓慢, 导致过冷状态下出现相分离(图 3), 且淬火温度仍然高于玻璃转化温度, 所以这种不均一现象被保留了下来(Philpotts, 2008)。另一种观点则认为这些不均一的玻璃是稳定的不混溶两相刚刚开始相分离的产物, 由于形核缓慢, 所以形成了两相共存的胶体(emulsion;Veksler et al., 2007)。

实验表明, 只要硅酸盐熔体的组成进入不混溶域, 很快就会出现纳米级的不混溶两相液滴(droplet), 但是随后2种液滴的生长变大则属于Ostwald熟化过程, 即相对小的液滴消失进而使相对大液滴增大(Mazurin and Porai-Koshits, 1984), 这一过程一般会伴随界面自由能γ的降低。在熟化过程中, 理论模拟显示液滴半径r的增大速度和熟化时间的平方根或立方根线性相关(James, 1975)。虽然液滴增大的过程也受熔体过饱和度和扩散率等因素影响, 但是在理论上有一点可以明确, 那就是界面自由能γ越大, 液滴增大的速度就越快。

Veksler等(2010)在简单体系K2O-FeO-Fe2O3-Al2O3-SiO2中测定的不混溶两相之间界面自由能是截至目前针对岩浆体系开展的唯一一组实验, 虽然实验是在空气中进行的, 环境过于氧化, 但依然有重要的参考价值。实验表明, 不混溶两相之间的界面自由能值很低, 只有8~16 mN/m, 甚至低于室温下水和家用食用油之间的界面能(30~50 mN/m)。这说明, 即使不混溶作用产生的两相能够很快成核形成小液滴, 但是它们通过Ostwald熟化生长增大却注定会很慢, 因此不混溶两相共存的胶体可以稳定存在很长时间。例如, 在玄武岩体系的实验中, 由于存在转熔反应(peritectic reaction): 橄榄石+熔体=易变辉石+斜长石+普通辉石。与这4种矿物平衡的液相的成分应该比较均一, 但实际情况远非如此, 其变化范围较大(Veksler, 2009Namur et al., 2012)。这很可能也是由于不混溶两相分离较慢所导致的。此外, 如果把这个过程放到自然界岩浆房尺度来考虑的话, 那么岩浆不混溶作用在大多数情况下都会产生富铁和富硅相共存的胶体, 但却很难在重力作用下在岩浆房里形成下部的富铁岩浆和上部的富硅岩浆发生(岩浆房尺度的液态不混溶)。这也是全晶质的自然样品(尤其是深成侵入岩)中很难观察到清晰的不混溶作用的岩相学证据的主要原因之一。

2.4 挥发分在岩浆液态不混溶作用中的影响

Lester等(2013a)的实验表明, 在简单体系K2O-FeO-Fe2O3-Al2O3-SiO2±An50中, 在较为氧化的环境下(MH和NNO缓冲剂), 水的加入会导致系统整体液相线温度的降低, 扩大磁铁矿的稳定区域并抑制硅酸盐矿物的结晶, 因此会导致不混溶域在相同氧逸度条件下的扩大, 甚至出现富铁相-富硅相-铁的硫化物的三相不混溶;磷和硫同样可以帮助扩大不混溶域。在磁铁矿-铁橄榄石-流纹岩的多组分系统中, 在100 MPa和1010℃, 较为还原的条件下(FMQ-1.5缓冲剂), 不论水活度怎样变化(aH2O=0.05~1), 氟饱和的岩浆系统(萤石成为液相线矿物)都会出现明显的不混溶(图 11)。

水活度为aH2O=0.05~1, 100 MPa、1 010℃中, 当F过饱和之后, 磁铁矿-铁橄榄石-流纹岩系统会发生液态不混溶作用, 其中磁铁矿和萤石会被富铁相熔体包裹(未发表实验);Mt-磁铁矿, Flu-萤石, LFe-富铁不混溶熔体, LSi-富硅不混溶熔体 图 11 不同水活度条件下磁铁矿-铁橄榄石-流纹岩多组分系统电子背散射图 Figure 11 BSE image showing liquid immiscibility occurred after the saturation of fluorite in the join magnetite-fayalite-rhyolite system
3 粒间熔体液态不混溶与层状岩体韵律层成因

基于上述背景, 下面以攀枝花岩体为例, 简要分析一下液态不混溶作用在岩体形成过程中可能发挥的重要作用。王坤等(2013)在攀枝花岩体中部带浅色辉长岩磷灰石中发现深色和浅色2类熔融包裹体, 指示了在磷灰石结晶时周围熔体发生了液态不混溶作用, 并通过熔融包裹体均一化和质量平衡计算分别得到富硅和富铁熔融包裹体的平均组成。根据这样的组成, Hou和Veksler(2015)使用常压气氛电炉, 实验验证结果表明富硅和富铁不混溶两相可以稳定存在(图 12a)。在岩浆不混溶作用发生的动力学过程表明, 岩浆房尺度的液态不混溶一般很难出现, 而岩相学观察这些磷灰石结晶比单斜辉石和斜长石都要晚, 因此磷灰石中熔融包裹体记录的是晶粥层中粒间熔体尺度的液态不混溶(图 12b)。液态不混溶出现在这个阶段也和粒间熔体相对更演化有关(图 6)。

Ab-钠长石;An-钙长石;Di-透辉石 (a)实验证实攀枝花岩体中部岩相带磷灰石熔融包裹体的成分是可以稳定存在的液态不混溶; (b)攀枝花岩体中部岩相带粒间熔体液态不混溶的模式图;(c)Ab-An-Di相图(引自Veksler et al., 2015) 图 12 粒间溶体液态不混溶作用与韵律成因模式 Figure 12 An integrated model for the formation of rhythmic layers and interstitial liquid immiscibility

在粒间熔体的尺度上, 液态不混溶发生时一开始以出现富硅和富铁小液滴为特征, 但是富硅小液滴由于密度小于粒间熔体, 从而会在浮力作用下向上迁移, 进入到上部以单斜辉石(cpx)-斜长石(plag)-粒间熔体(L)三者为主要组分的体系中(在缓慢降温的深成条件下, 此体系可以看作是平衡的)。富硅的不混溶相, 除了富硅以外, 还富集铝和碱, 因此成分上类似钠长石Ab。这与开放体系的Di-Ab-An三元系非常类似, 在图 12c中, 假设富硅液滴加入之前, 上部晶粥层的总组成相当于m点, 单斜辉石的组成为C1, 斜长石的组成为P1, 粒间熔体的组成为L1。如果富硅液滴进入体系之后, 总组成会从m点向q点移动, 为了“缓冲”系统中加入的硅铝质成分, 富镁铁质的单斜辉石会溶解。当更多Ab成分加入使总组成移动到q点时, 系统中单斜辉石消失殆尽, 只剩下L1和P1两相, 如果再继续有富硅熔体进入使总组成移动到r点, 会有L2和P2重新平衡。不论总组成在q还是r, 都会导致系统中相对缺失暗色堆晶矿物, 从而导致浅色层的形成。而对于富铁液滴, 不论留在原地还是在重力作用下向下迁移, 伴随着粒间熔体液态不混溶作用的持续进行, 将会导致浅色矿物的消失和暗色层的形成。二者原理类似, 只是m点会在单斜辉石首晶区内(图 12c)。

4 存在问题和今后研究方向

在岩浆岩岩石学研究领域, 虽然越来越多的证据表明液态不混溶作用的发生以及其对岩浆演化路径(liquid line of descent)的重要性, 但是岩浆液态不混溶作用发生的具体机制还存在争议。两种不同密度和黏度的硅酸盐熔体能否在岩浆房中发生大规模的分离还没有得到实验证实, 不过正如本文所述, 小规模的不混溶作用(如粒间熔体尺度)还是比较有说服力的。

岩浆液态不混溶作用的规模不论大小, 由于其可以相对快速地富集某些元素, 所以对于矿床的形成, 尤其是铁、钛、磷等矿床的形成可能有着重要的意义。但值得注意的是在目前的相关研究中, 普遍存在将堆晶(cumulate)和熔体(或液相, liquid或melt)彼此混淆的现象。实际上, 2种成分截然的彼此相邻的堆晶岩石并不代表2种不混溶的熔体。

今后的研究方向可以针对岩石和矿床的形成机制, 集中在不混溶作用发生的诱因, 规模以及相分离等动力学过程。

致谢: 感谢中国地质大学(北京)张招崇教授, 德国汉诺威大学Francois Holtz教授, 德国GFZ波茨坦地球科学研究中心Ilya Veksler博士和比利时列日大学Bernard Charlier研究员长期的支持和指导, 以及匿名审稿人和编辑。

参考文献
[] Bogaerts M, Schmidt M W. 2006. Experiments on silicate melt immiscibility in the system Fe2SiO4-KAlSi3O8-SiO2-CaO-MgO-TiO2-P2O5 and implications for natural magmas. Contributions to Mineralogy and Petrology , 152 (3) : 257–274. DOI:10.1007/s00410-006-0111-6
[] Botcharnikov R E, Almeev R R, Koepke J, Holtz F. 2008. Phase relations and liquid lines of descent in hydrous ferrobasalt-implications for the Skaergaard intrusion and Columbia River flood basalts. Journal of Petrology , 49 (9) : 1687–1727. DOI:10.1093/petrology/egn043
[] Bowen N L. 1928. The evolution of the igneous rocks. Princeton: Princeton University Press: 332 .
[] Charlier B, Grove T L. 2012. Experiments on liquid immiscibility along tholeiitic liquid lines of descent. Contributions to Mineralogy and Petrology , 164 (1) : 27–44. DOI:10.1007/s00410-012-0723-y
[] Charlier B, Namur O, Grove T L. 2013. Compositional and kinetic controls on liquid immiscibility in ferrobasalt-rhyolite volcanic and plutonic series. Geochimica et Cosmochimica Acta , 113 : 79–93. DOI:10.1016/j.gca.2013.03.017
[] Corgne A, Keshav S, Wood B J, McDonough W F, Fei Y W. 2008. Metal-silicate partitioning and constraints on core composition and oxygen fugacity during Earth accretion. Geochimica et Cosmochimica Acta , 72 (2) : 574–589. DOI:10.1016/j.gca.2007.10.006
[] Daly R A. 1914. Igneous rocks and their origin. NewYork:McGraw-Hill Book Company, Inc., 563
[] Dasgupta R, Hirschmann M M, Stalker K. 2006. Immiscible transition from carbonate-rich to silicate-rich melts in the 3 GPa melting interval of eclogite+CO2 and genesis of silica-undersaturated ocean island lavas. Journal of Petrology , 47 (4) : 647–671.
[] De A. 1974. Silicate liquid immiscibility in the Deccan traps and its petrogenetic significance. Geological Society of America Bulletin , 85 (3) : 471–474. DOI:10.1130/0016-7606(1974)85<471:SLIITD>2.0.CO;2
[] Dixon S, Rutherford M J. 1979. Plagiogranites as late-stage immiscible liquids in ophiolite and mid-ocean ridge suites:An experimental study. Earth and Planetary Science Letters , 45 (1) : 45–60. DOI:10.1016/0012-821X(79)90106-7
[] Doremus R H. 1994. Glass science.2nd ed. New York: John Wiley & Sons Inc .
[] Freestone I C, Hamilton D L. 1980. The role of liquid immiscibility in the genesis of carbonatites-An experimental study. Contributions to Mineralogy and Petrology , 73 (2) : 105–117. DOI:10.1007/BF00371385
[] Greig J W. 1927. Immiscibility in silicate melts. American Journal of Science , 13 (73) : 133–154.
[] Hess P C. 1995. Thermodynamic mixing properties and the structure of silicate melts. In:Stebbins J F, McMillan P F, Dingwell D B, eds. Structure, Dynamics & Properties of Silicate Melts. Washington, D C:Mineralogical Society of America , 32 : 145–190.
[] Hess P C. 1996. Permian and lower critical points:Thermodynamic constraints on the solution properties of silicate melts. Geochimica et Cosmochimica Acta , 60 (13) : 2365–2377. DOI:10.1016/0016-7037(96)00095-6
[] Hoover J D, Irvine T N. 1977. Liquidus relations and Mg-Fe partitioning on part of the system Mg2SiO4-Fe2SiO4-CaMgSi2O6-CaFeSi2O6-KAlSi3O8-SiO2. Carnegie Institution of Washington Yearbook , 77 : 774–784.
[] Hou T, Veksler I V. 2015. Experimental confirmation of high-temperature silicate liquid immiscibility in multicomponent ferrobasaltic systems. American Mineralogist , 100 (5-6) : 1304–1307. DOI:10.2138/am-2015-5285
[] Hudon P, Baker D R. 2002. The nature of phase separation in binary oxide melts and glasses. I. Silicate systems. Journal of Non-Crystalline Solids , 303 (3) : 299–345. DOI:10.1016/S0022-3093(02)01043-8
[] Irvine T N. 1976. Metastable liquid immiscibility and MgO-FeO-SiO2 fractionation patterns in the system Mg2SiO4-Fe2SiO4-CaAl2Si2O8-KAlSi3O8-SiO2. Carnegie Institution of Washington Yearbook , 75 : 597–611.
[] Jakobsen J K, Veksler I V, Tegner C, Brooks C K. 2005. Immiscible iron-and silica-rich melts in basalt petrogenesis documented in the Skaergaard intrusion. Geology , 33 (1) : 885–888.
[] Jakobsen J K, Veksler I V, Tegner C, Brooks C K. 2011. Crystallization of the Skaergaard intrusion from an emulsion of immiscible iron-and silica-rich liquids:Evidence from melt inclusions in plagioclase. Journal of Petrology , 52 (2) : 345–373. DOI:10.1093/petrology/egq083
[] James P F. 1975. Liquid-phase separation in glass-forming systems. Journal of Materials Science , 10 (10) : 1802–1825. DOI:10.1007/BF00554944
[] Juster T C, Grove T L, Perfit M R. 1989. Experimental constraints on the generation of FeTi basalts, andesites, and rhyodacites at the Galapagos spreading center, 85°W and 95°W. Journal of Geophysical Research:Solid Earth(1978-2012) , 94 (B7) : 9251–9274. DOI:10.1029/JB094iB07p09251
[] Kamenetsky V S, Charlier B, Zhitova L, Sharygin V, Davidson P, Feig S. 2013. Magma chamber-scale liquid immiscibility in the Siberian Traps represented by melt pools in native iron. Geology , 41 (10) : 1091–1094. DOI:10.1130/G34638.1
[] Kyser T K, Lesher C E, Walker D. 1998. The effects of liquid immiscibility and thermal diffusion on oxygen isotopes in silicate liquids. Contributions to Mineralogy and Petrology , 133 (4) : 373–381. DOI:10.1007/s004100050459
[] Lester G W, Clark A H, Kyser T K, Naslund H R. 2013a. Experiments on liquid immiscibility in silicate melts with H2O, P, S, F and Cl:Implications for natural magmas. Contributions to Mineralogy and Petrology , 166 (1) : 329–349. DOI:10.1007/s00410-013-0878-1
[] Lester G W, Kyser T K, Clark A H. 2013b. Oxygen isotope partitioning between immiscible silicate melts with H2O, P and S. Geochimica et Cosmochimica Acta , 109 : 306–311. DOI:10.1016/j.gca.2013.01.037
[] Levin E M, Robbins C R, McMurdie H F. 1964. Phase equilibria diagrams, vol 1. Columbus:American Ceramic Society
[] Mazurin O V, Porai-Koshits E A. 1984. Phase separation in glass. Amsterdam, North-Holland, Oxford, New York, Tokyo 382
[] McBirney A R, Nakamura Y. 1974. Immiscibility in late stage magmas of the Skaergaard intrusion. Carnegie Institution of Washington Yearbook , 73 : 348–354.
[] Naldrett A J. 1999. World-class Ni-Cu-PGE deposits:Key factors in their genesis. Mineralium Deposita , 34 (3) : 227–240. DOI:10.1007/s001260050200
[] Namur O, Charlier B, Holness M B. 2012. Dual origin of Fe-Ti-P gabbros by immiscibility and fractional crystallization of evolved tholeiitic basalts in the Sept Iles layered intrusion. Lithos , 154 : 100–114. DOI:10.1016/j.lithos.2012.06.034
[] Philpotts A R. 1979. Silicate liquid immiscibility in tholeiitic basalts. Journal of Petrology , 20 (1) : 99–118. DOI:10.1093/petrology/20.1.99
[] Philpotts A R. 1982. Compositions of immiscible liquids in volcanic rocks. Contributions to Mineralogy and Petrology , 80 (3) : 201–218. DOI:10.1007/BF00371350
[] Philpotts A R. 2008. Comments on:Liquid immiscibility and the evolution of basaltic magma. Journal of Petrology , 49 (12) : 2171–2175. DOI:10.1093/petrology/egn061
[] Roedder E, Weiblen P W. 1970. Silicate liquid immiscibility in lunar magmas, evidenced by melt inclusions in lunar rocks. Science , 167 (3918) : 641–644. DOI:10.1126/science.167.3918.641
[] Roedder E. 1951. Low-temperature liquid immiscibility in the system K2O-FeO-Al2O3-SiO2. American Mineralogist , 36 : 282–286.
[] Roedder E. 1978. Silicate liquid immiscibility in magmas and in the system K2O-FeO-Al2O3-SiO2:An example of serendipity. Geochimica et Cosmochimica Acta , 42 : 1597–1617. DOI:10.1016/0016-7037(78)90250-8
[] Roedder E. 1979. Silicate liquid immiscibility in magmas. In:Yoder H S Jr, ed. The Evolution of Igneous Rocks. Fiftieth Anniversary Perspectives. Princeton:Princeton University Press, 15-58
[] Schmidt M W, Connolly J A D, Günther D, Bogaerts M. 2006. Element partitioning:The role of melt structure and composition. Science , 312 (5780) : 1646–1650. DOI:10.1126/science.1126690
[] Shelby J E. 2005. Introduction to glass science and technology.2nd ed. Cambridge: The Royal Society of Chemistry .
[] Simmons J H, Napolitano A, Macedo P B. 1970. Supercritical viscosity anomaly in oxide mixtures. Journal of Chemical Physics , 53 (3) : 1165–1170. DOI:10.1063/1.1674114
[] Simmons J H, Macedo P B. 1971. Analysis of viscous relaxation in critical oxide mixtures. Journal of Chemical Physics , 54 (3) : 1325–1331. DOI:10.1063/1.1674972
[] Thompson A B, Aerts M, Hack A C. 2007. Liquid immiscibility in silicate melts and related systems. Reviews in Mineralogy and Geochemistry , 65 (1) : 99–127. DOI:10.2138/rmg.2007.65.4
[] Toplis M J, Carroll M R. 1995. An experimental study of the influence of oxygen fugacity on Fe-Ti oxide stability, phase relations, and mineral-melt equilibria in ferro-basaltic systems. Journal of Petrology , 36 (5) : 1137–1170. DOI:10.1093/petrology/36.5.1137
[] Veksler I V, Dorfman A M, Rhede D, Wirth R, Borisov A A, Dingwell D B. 2008a. Liquid unmixing kinetics and the extent of immiscibility in the system K2O-CaO-FeO-Al2O3-SiO2. Chemical Geology , 256 (3-4) : 119–130. DOI:10.1016/j.chemgeo.2008.06.033
[] Veksler I V, Dorfman A M, Borisov A A, Wirth R, Dingwell D B. 2007. Liquid immiscibility and the evolution of basaltic magma. Journal of Petrology , 48 (11) : 2187–2210. DOI:10.1093/petrology/egm056
[] Veksler I V, Dorfman A M, Rhede D, Wirth R, Borisov A A, Dingwell D B. 2008a. Liquid unmixing kinetics and the extent of immiscibility in the system K2O-CaO-FeO-Al2O3-SiO2. Chemical Geology , 256 (3-4) : 119–130. DOI:10.1016/j.chemgeo.2008.06.033
[] Veksler I V, Dorfman A M, Borisov A A, Wirth R, Dingwell D B. 2008b. Liquid immiscibility and evolution of basaltic magma:Reply to S. A. Morse, A. R. McBirney and A. R. Philpotts. Journal of Petrology , 49 (12) : 2177–2186.
[] Veksler I V. 2009. Extreme iron enrichment and liquid immiscibility in mafic intrusions:Experimental evidence revisited. Lithos , 111 (1-2) : 72–82. DOI:10.1016/j.lithos.2008.10.003
[] Veksler I V, Kähn J, Franz G, Dingwell D B. 2010. Interfacial tension between immiscible liquids in the system K2O-FeO-Fe2O3-Al2O3-SiO2 and implications for the kinetics of silicate melt unmixing. American Mineralogist , 95 (11-12) : 1679–1685. DOI:10.2138/am.2010.3456
[] Veksler I V, Reid D L, Dulski P, Keiding J K, Schannor M, Hecht L, Trumbull R B. 2015. Electrochemical Processes in a Crystal Mush:Cyclic Units in the Permian Critical Zone of the Bushveld Complex, South Africa. Journal of Petrology , 56 (6) : 1229–1250. DOI:10.1093/petrology/egv036
[] Veksler I V, Charlier B. 2015. Silicate liquid immiscibility in layered intrusions. In:Charlier B, Namur O, Latypov R, Tegner C(eds). Layered Intrusions. Netherlands:Springer, 229-258
[] Veksler I V, Dorfman A M, Danyushevsky L V, Jakobsen J K, Dingwell D B. 2006. Immiscible silicate liquid partition coefficients:Implications for crystal-melt element partitioning and basalt petrogenesis. Contributions to Mineralogy and Petrology , 152 (6) : 685–702. DOI:10.1007/s00410-006-0127-y
[] Visser W, van Groos A F K. 1976. Liquid immiscibility in K2O-FeO-Al2O3-SiO2. Nature , 264 (5585) : 426–427. DOI:10.1038/264426a0
[] Visser W, van Groos A F K. 1977. Liquid immiscibility in K2O-FeO-Al2O3-SiO2(reply). Nature , 267 (5611) : 560. DOI:10.1038/267560a0
[] Visser W, van Groos A F K. 1979. Phase relations in the system K2O-FeO-Al2O3-SiO2 at 1 atmosphere with special emphasis on low temperature liquid immiscibility. American Journal of Science , 27 (1) : 70–91.
[] Yoder H S Jr, Tilley C E. 1962. Origin of basalt magmas:An experimental study of natural and synthetic rock systems. Journal of Petrology , 3 (3) : 342–532. DOI:10.1093/petrology/3.3.342
[] Zhu D, Bao H M, Liu Y. 2015. Non-traditional stable isotope behaviors in immiscible silica-melts in a mafic magma chamber. Scientific Reports , 5 : 17561. DOI:10.1038/srep17561
[] 王坤, 邢长明, 任钟元, 王焰. 2013. 攀枝花镁铁质层状岩体磷灰石中的熔融包裹体:岩浆不混熔的证据. 岩石学报 , 29 (10) : 3503–3518.