Skip Navigation



Journal of Petrology Advance Access published online on October 6, 2007

Journal of Petrology, doi:10.1093/petrology/egm056
This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
48/11/2187    most recent
egm056v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Veksler, I. V.
Right arrow Articles by Dingwell, D. B.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?

© The Author 2007. Published by Oxford University Press. All rights reserved. For Permissions, please e-mail: journals.permissions@oxfordjournals.org

Liquid Immiscibility and the Evolution of Basaltic Magma

Ilya V. Veksler1,*, Alexander M. Dorfman2, Alexander A. Borisov3,4, Richard Wirth1 and Donald B. Dingwell2,5

1Geoforschungszentrum Potsdam, Department 4.1, Telegrafenberg, 14473 Potsdam, Germany
2Earth and Environmental Sciences, University of Munich, Theresienstrasse 41, 80333 Munich, Germany
3Institute of Geology and Mineralogy, University of Cologne, Zülpicherstrasse, 49 B, 50674 Cologne, Germany
4Institute of Geology of Ore Deposits, Petrography, Mineralogy and Geochemistry, Russian Academy of Sciences, Staromonetny 35, 109017 Moscow, Russia
5Geological and Environmental Sciences, Stanford University, Stanford, CA, 94305-2115 USA

Received April 2, 2007; Revised typescript accepted August 17, 2007


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
This experimental study examines relationships between alternative evolution paths of basaltic liquids (the so-called Bowen and Fenner trends), and silicate liquid immiscibility. Synthetic analogues of natural immiscible systems exhibited in volcanic glasses and melt inclusions were used as starting mixtures. Conventional quench experiments in 1 atm gas mixing furnaces proved unable to reproduce unmixing of ferrobasaltic melts, yielding instead either turbid, opalescent glasses, or crystallization of tridymite and pyroxenes. In contrast, experiments involving in situ high-temperature centrifugation at 1000g (g = 9·8 m/s2) did yield macroscopic unmixing and phase separation. Centrifugation for 3–4 h was insufficient to complete phase segregation, and resulted in sub-micron immiscible emulsions in quenched glasses. For a model liquid composition of the Middle Zone of the Skaergaard intrusion at super-liquidus temperatures of 1110–1120°C, centrifugation produced a thin, silicic layer (64·5 wt% SiO2 and 7·4 wt% FeO) at the top of the main Fe-rich glass (46 wt% SiO2 and 21 wt% FeO). The divergent compositions at the top and bottom were shown in a series of static runs to crystallize very similar crystal assemblages of plagioclase, pyroxene, olivine, and Fe–Ti oxides. We infer from these results that unmixing of complex aluminosilicate liquids may be seriously kinetically hampered (presumably by a nucleation barrier), and thus conventional static experiments may not correctly reproduce it. In the light of our centrifuge experiments, immiscibility in the Skaergaard intrusion could have started already at the transition from the Lower to the Middle Zone. Thus, magma unmixing might be an important factor in the development of the Fe-enrichment trend documented in the cumulates of the Skaergaard Layered Series.

KEY WORDS: liquid immiscibility; Skaergaard; layered intrusions; experimental petrology


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Very few ideas in the history of igneous petrology have experienced such dramatic turns in acceptability as the concept of igneous petrogenesis by liquid immiscibility. A popular idea at the dawn of modern petrology (e.g. Daly, 1914Go), it fell into disrepute in late 1920s (Greig, 1927Go; Bowen, 1928Go), was revived by new evidence in the 1970s and 1990s (Roedder & Weiblen, 1970Go, 1971Go; De, 1974Go; McBirney & Nakamura, 1974Go; Philpotts, 1976Go, 1982Go; Dixon & Rutherford, 1979Go; Roedder, 1979Go), and has remained largely ignored during the last two to three decades. Nevertheless, a modern re-evaluation of petrogenetic and geochemical aspects of silicate immiscibility has been attempted recently (e.g. Jakobsen et al., 2005Go; Bogaerts & Schmidt, 2006Go; Schmidt et al. 2006Go; Veksler et al., 2006Go). Observations reported in some of the recent papers imply that the phenomenon of melt unmixing may hold a key to a number of major, long-standing problems of igneous petrology.

One early recognized and hotly debated problem of igneous petrology has been the Bowen–Fenner controversy. The controversy dates back to the classical works by Bowen (1928Go) and Fenner (1929Go), and it deals with the general chemical trend of evolution of sub-alkaline basaltic magma. Simply put, Fenner and his followers since that time have advocated a liquid differentiation trend towards strong enrichment in Fe-oxides at constant or slightly decreasing silica content. Bowen led proponents of the alternative claim that the content of Fe-oxides in basaltic liquids never achieves the concentrations inferred from the Fenner trend, and the liquid evolution is characterized rather by silica enrichment. The trends, which are also called tholeiitic and calc-alkaline, both exist in nature and result in the formation of distinct volcanic rock associations and plutonic complexes (e.g. Wager & Deer, 1939Go; Yoder & Tilley, 1962Go; Carmichael, 1964Go). The processes behind these opposing trends have been traditionally explained in terms of fractional crystallization, with particular emphasis on (1) the onset of magnetite crystallization, which depends on the oxidation state of magma (e.g. Osborn, 1979Go), and (2) cotectic proportions of olivine, plagioclase and clinopyroxene (Grove & Baker, 1984Go). Others have stressed the effects of melt chemistry; for example, the role of alkalis (Irvine, 1976Go; Shi, 1993Go).

Silica and Fe-oxides, the definitive components of the Bowen and Fenner trends, are also key players in the occurrence of silicate liquid immiscibility (Roedder, 1951Go). During the last surge of interest in silicate melt unmixing 20–30 years ago, it was established that certain ferrobasaltic liquids were unstable, unmixing to form two contrasting liquid compositions, one with a SiO2 content of less than 45 wt% and very high total iron [FeO(t) = FeO + Fe2O3] at 25–30 wt% or more, and the other with SiO2 contents at or exceeding 55–60 wt%, and FeO(t) at about 5–10 wt% (e.g. Philpotts, 1976Go, 1979Go, 1982Go). Observations both in natural rocks and in laboratory experiments have confirmed the stable coexistence of certain silicic and Fe-rich liquids with gabbroic mineral assemblages (Roedder & Weiblen, 1970Go; De, 1974Go; McBirney & Nakamura, 1974Go; Dixon & Rutherford, 1979Go; Philpotts, 1979Go; Roedder, 1979Go; Longhi, 1990Go). Silicate liquid immiscibility in ferrobasaltic melts results in the formation of (1) extremely Fe-rich, and (2) silicic conjugate liquids. Thus, it has the potential to place the Bowen vs Fenner debate on a new basis, possibly linking the divergent liquid evolution trends. In the presence of immiscibility, the trends of FeO(t) and silica enrichment do not appear as mutually exclusive alternatives, but rather as complementary evolution paths of immiscible conjugate liquids in a heterogeneous magma. In this debate experimental evidence is crucial. It should be pointed out that reliable experimental demonstrations of immiscibility in basaltic compositions have been restricted so far to temperatures below 1050°C, and very advanced stages of crystallization. As a consequence, immiscibility has been relegated to a minor role in igneous petrogenesis. Unless it is proven that immiscibility in natural ferrobasaltic melts can start earlier (i.e. at about 50–60% crystallization of a common tholeiitic liquid, and temperatures around 1070–1100°C), unmixing will remain insignificant in the Bowen vs Fenner debate.

In this study, we examine the experimental evidence for and against earlier unmixing of ferrobasaltic liquids. We revisit the thermodynamic and phase equilibria constraints on the Bowen and Fenner types of liquid evolution. We present new results for static and dynamic (centrifuge) experiments, which demonstrate that unmixing of multicomponent basaltic liquids appears to be seriously experimentally hampered by a nucleation barrier, formation of colloidal emulsions, and metastable crystallization. Our results show that conventional static experiments on small melt droplets or rock chips in metal wire loops may be unable to correctly reproduce the onset of liquid immiscibility in evolved, Fe-rich basaltic or dacitic liquids. The onset of macroscopic immiscibility may be pushed to lower temperatures by metastable crystallization, especially in charges with high initial volume proportion of crystals (e.g. in experiments on rock chips). In some cases, a second immiscible liquid phase may not nucleate at all, or not grow beyond the sub-micron, colloidal stage, and after hours and days it may remain unrecognizable in final glasses from fine exsolutions formed during quenching. High-temperature in situ centrifugation is, however, capable of separating immiscible emulsions, and thus revealing stable, super-liquidus immiscibility even when the result is a very fine emulsion. In light of our centrifuge experiments, we propose that unmixing of natural ferrobasaltic liquids may in fact start at a much earlier stage of crystallization than was previously thought. Our starting mixtures were based on compositions of natural immiscible glasses, and so we begin with a brief overview of those natural melt inclusions that appear to record the early stages of unmixing. In conclusion, we examine the petrogenetic implications for tholeiitic magmas, mostly using the classical example of the Skaergaard intrusion, East Greenland, and propose that the existence of the divergent trends primarily stems from the strong non-ideality of ferrobasaltic and ferrodacitic liquids.


    SILICATE LIQUID IMMISCIBILITY IN MAGMAS AND SELECTION OF EXPERIMENTAL COMPOSITIONS
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Volcanic glasses and melt inclusions
Droplet exsolution textures in natural volcanic glasses have been traditionally regarded as the strongest evidence supporting silicate liquid immiscibility. Classical examples were presented by Philpotts (1982Go), who carried out a detailed study of the mesostasis of basalts and andesites from numerous localities worldwide, and described glassy, or partly crystallized spherical droplets dispersed in a second glass of a different composition. The compositions of the droplets and matrix glasses were shown to compare well with experimental Fe-rich and silica-rich immiscible conjugate melts, consistent with an origin by liquid immiscibility. The average compositions of the coexisting phases in the mesostasis of tholeiitic lavas according to Philpotts (1982Go) are listed in the first two columns of Table 1. Both phases are notably depleted in MgO, and imply equilibration at low temperature. Experiments on similar basaltic compositions (Roedder & Weiblen, 1970Go; McBirney & Nakamura, 1974Go; Dixon & Rutherford, 1979Go; Philpotts, 1979Go; Philpotts & Doyle, 1983Go; Longhi, 1990Go), confirmed the late-stage, low-temperature origin of the groundmass phase assemblages forming the immiscible globules. According to those experiments, immiscibility in terrestrial basaltic lavas probably occurs below 1020°C.


View this table:
[in this window]
[in a new window]

 
Table 1: Examples of natural immiscible liquids (quenched glasses) in mesostasis of basaltic lavas, and in melt inclusions (wt% oxides)

 
Significantly less evolved, more refractory compositions have been reported, however, in glassy melt inclusions hosted by olivine and plagioclase phenocrysts (Roedder & Weiblen, 1970Go; Krasov & Clocchiatti, 1979Go; Roedder, 1979Go; Fujii et al., 1980Go; Philpotts, 1981Go, 1982Go; Shearer et al., 2001Go). Philpotts (1981Go, 1982Go) argued that the more primitive glass compositions in plagioclase-hosted inclusions may be due to metastable unmixing. He proposed that early exsolution of pyroxenitic, Fe-rich liquid in plagioclase-hosted melt inclusions could be a consequence of arrested crystallization of clinopyroxene in small inclusion cavities because of the absence of suitable sites for nucleation. Such an explanation is, however, difficult to apply to the case studied by Krasov & Clocchiatti (1979Go). Those workers attempted to homogenize inclusions in An70 plagioclase phenocrysts from subduction-related andesite–dacitic lava of the Karymsky volcano, Kamchatka, Russia, and observed two conjugate silicate liquids at temperatures of up to 1250–1280°C. Importantly, some of the immiscible inclusions had been originally crystallized to fine-grained daughter crystal assemblages, and immiscibility was shown to be reversible in heating–cooling experiments around 1280°C. The compositions of the quenched immiscible glass phases reported by Krasov & Clocchiatti (1979Go) are also presented in Table 1. The glasses are characterized by elevated MgO contents, and a lower contrast in SiO2 and Fe-oxide contents. These features are consistent with the equilibration of immiscible liquids at higher temperatures.

Perhaps the most convincing example of high-temperature silicate liquid immiscibility is provided by melt inclusions in native iron, which, although exotic in terrestrial environments, do occur in a few basaltic sills and lava flows. Grains, nuggets, and large blocks of native iron and minor cohenite (Fe3C) have been reported in basaltic lavas on Disco island, West Greenland (Bird et al., 1981Go), and in the roof endocontact zone of the Khutungunsky sill of the Siberian trap province (Ryabov, 1988Go, 1989Go). The occurrences of native iron are believed to result from reduction of basaltic magma by carbon at the contact with organic-bearing sediments. Some iron nuggets contain numerous inclusions of rock fragments, single silicate crystals, and spherical glass droplets showing clear, classical immiscibility textures, as can be seen in Fig. 1. The photograph shows that the Fe-rich and silica-rich glasses are very distinct, and, in those inclusions where low-temperature devitrification is minor, the interfaces between the glasses are razor-sharp. Representative electron microprobe data on the glasses in the Siberian samples from Ryabov (1989Go) are also presented in Table 1. Very similar compositions have been reported from Disco island (Bird et al., 1981Go). The elevated MgO contents of the glasses are similar to those of the plagioclase-hosted melt inclusions studied by Krasov & Clocchiatti (1979Go; see also Table 1), and imply a high-temperature origin.


Figure 1
View larger version (66K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 1. Natural inclusions of immiscible Fe-rich (Lfe) and silica-rich (Lsi) glasses in native iron, Siberian Traps; reflected light. The Fe-rich glass forms a rim around the silica-rich glass, and there are also a few smaller spherical droplets of Lfe inside the Lsi. Glass compositions reported by Ryabov (1989Go) are given in Table 1.

 
Two of our starting mixtures were based on the average compositions of the immiscible glasses of the iron-hosted melt inclusions. The compositions RY40 and RY20 (Table 2) were calculated as mixtures between the average Fe-rich (Lfe) and silica-rich (Lsi) iron-hosted glasses published by Ryabov (1989Go) (Table 1) in mass proportions of 3:2 and 4:1. In other words, the composition RY40 upon complete melting should theoretically yield 40 wt% of the conjugate liquid Lsi, and 60 wt% Lfe, whereas the composition RY20 should give 20 and 80 wt% respectively of the same liquids.


View this table:
[in this window]
[in a new window]

 
Table 2: Compositions of synthetic starting mixtures normalized to 100 wt%

 
Three other starting mixtures were designed in a similar way from compositions of other immiscible glasses listed in Table 1. The composition SF-1 is based on the analyses of low-temperature immiscible groundmass globules and glasses published by Philpotts (1982Go), and compiled in the two leftmost columns of Table 1. Accordingly, the composition has the lowest MgO content among the mixtures used in this study. The mixture SF-1 has been used in previous centrifuge experiments (Veksler et al., 2006Go), where it produced roughly equal proportions of Lfe and Lsi immiscible phases. The composition KC is based on electron microprobe analyses of the Fe-rich and silica-rich glasses of the plagioclase-hosted melt inclusions published by Krasov & Clocchiatti (1979Go). The mass proportion of the Lfe and Lsi (see the analyses in Table 1) is 7:3. The mixture FJ is based on electron microprobe analyses of immiscible glassy inclusions in plagioclase phenocrysts from an olivine-bearing andesite of the Akita-Komagatake volcano, Japan, which were published by Fujii et al. (1980Go). The analyses were carried out using defocused beam, and the composition FJ represents an average mixture of the Fe-rich and silica-rich immiscible liquids in unidentified proportions. Notably, the MgO content in the bulk of the inclusions is the highest ever reported among immiscible liquids.

By their very nature, plagioclase-hosted melt inclusions must be saturated in their host plagioclase. The compositions KC and FJ are, however, surprisingly low in alumina, and, as we learned from our experiments on these mixtures (see below), do not produce plagioclase on the liquidus. To compensate for possible post-entrapment modification as a result of precipitation of plagioclase on inclusion walls [see Veksler (2006Go) for discussion], we prepared another mixture FJ20, which has been derived from the composition FJ by the addition of 20 wt% of pure anorthite.

Liquids of the Skaergaard intrusion
Since the study by Wager & Deer (1939Go), the Skaergaard intrusion of layered gabbros has been the prime example of the Fenner, tholeiitic, or Fe-enrichment trend. Thus, Skaergaard has remained for decades at the centre of the Bowen–Fenner controversy. The intrusion is located in East Greenland, and after many years of comprehensive research was once termed ‘the most intensely studied body of rock on Earth’ (McBirney & Naslund, 1990Go). Numerous attempts have been made to reconstruct liquid lines of descent in the intrusion by: (1) mass-balance calculations (Wager & Brown, 1968Go; Hunter & Sparks, 1987Go; Nielsen, 2004Go); (2) experimental simulation of crystallization or partial melting of Skaergaard rocks (McBirney & Naslund, 1990Go; Toplis & Carroll, 1995Go; Thy et al., 2006Go); (3) tracing specific elements using mineral–melt partition coefficients (Tegner, 1997Go; Jang et al., 2001Go); (4) geochemical thermometry of cumulus mineral assemblages (Ariskin, 1999Go, 2002Go); (5) studies of crystallized melt inclusions in cumulus minerals (Hanghøj et al., 1995Go; Jakobsen et al., 2005Go, 2007Go); (6) comparisons with compositions of coeval and spatially associated dykes (Brooks & Nielsen, 1978Go, 1990Go). The disparate results of these reconstructions are illustrated in Fig. 2. In general, one may conclude that the disagreement among researchers is the strongest when it comes to the late stages of magma evolution, especially after the onset of ilmenite and magnetite crystallization. The disagreement between the models is, however, not so dramatic regarding the starting composition of the Skaergaard parental magma (see also Nielsen, 2004Go), and the models also agree reasonably well at the early stages of magma evolution, up to about 50–60% crystallization. This stage corresponds roughly to the transition from the Lower to the Middle Zone of the Layered Series marked by the disappearance of olivine from the liquidus mineral association and the appearance of liquidus pigeonite [see Wager & Brown (1968Go) and McBirney (1989Go) for a detailed description of the Skaergaard stratigraphy]. The liquid composition at this stage has been well constrained by detailed experimental studies (McBirney & Naslund, 1990Go; Toplis & Carroll, 1995Go; Thy et al., 2006Go), thermodynamic modelling (Ariskin, 2002Go), and mass-balance calculations (Nielsen, 2004Go). Our starting mixtures MZ-1 and MZ-2 (Table 2) are based on estimates of the liquid composition at the top of the Lower Zone. The composition MZ-1 is very similar to the glasses produced in experiments by Toplis & Carroll (1995Go) at temperatures close to 1100°C. The composition MZ-2 differs from MZ-1 only by its somewhat higher FeO content, which brings it closer to the mass-balance estimation of the Middle Zone liquid by Nielsen (2004Go). The remaining two starting mixtures in Table 2, the Fe-rich composition MZF and silica-rich composition MZS, are explained below, where we describe the results of our centrifuge experiments.


Figure 2
View larger version (24K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 2. Liquid lines of descent in the Skaergaard intrusion in terms of total iron (FeO) vs silica (SiO2) variations. Average compositions of immiscible Fe-rich and Si-rich melt inclusions in apatite reported by Jakobsen et al. (2005Go) are shown by filled and open stars.

 

    EXPERIMENTAL AND ANALYTICAL METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Preparation of starting mixtures
All the starting mixtures listed in Table 2 were synthesized by fusion in Pt crucibles of carefully weighed and mixed reagent-grade chemicals [SiO2, Al2O3, MgO, TiO2, CaCO3, Na2CO3, K2CO3 and Ca3(PO4)2]. To avoid strong oxidation and Fe losses to the Pt crucibles, Fe-free glasses were first prepared by repeated fusions in an electric furnace once at 900°C, and then twice at 1350°C. The Fe-free glasses were then crushed to a fine-grained powder (grain size less than 3 µm) and mixed and ground with reagent-grade FeO in an agate mortar under acetone. Fine powders of the glass–FeO mixtures were used for the experimental charges.

Static experiments
Static experiments were conducted in a 1 atm vertical tube furnace in which fO2 was controlled by CO–CO2 gas mixtures. A PtRh6–PtRh30 (type B) thermocouple was employed for temperature measurements. The uncertainties of cited log fO2 values and temperature do not exceed ±0·2 and ±2°C, respectively. In most cases, silicate melts were suspended from Re loops, prepared from commercially available Re ribbon (Rhenium Alloys, Inc.). Rhenium has been shown to be an excellent loop material, preventing iron loss from a charge even at very reducing conditions (Borisov & Jones, 1999Go). The samples were either (1) directly placed in a pre-heated furnace at the required temperature or (2) first over-heated by 50–150°C and then cooled to run temperature at a rate of 10–50°C per hour. After sufficient dwell exposure at constant temperature, the samples were quenched in air within a few seconds.

Centrifuge experiments
The experimental equipment and Fe–Pt double containers for the centrifuge experiments have been described in detail previously (e.g. Dorfman et al., 1996Go; Veksler et al., 2006Go). Briefly, the equipment represents a small cylindrical wire-wound electric furnace mounted on a centrifuge, which provides acceleration of up to 1000g (g = 9·8 m/s2) over maximal run durations of 5 h. Temperature is controlled by three Pt–PtRh10 (type S) thermocouples at the top, bottom and the middle of the container, respectively. Comparisons between the thermocouple readings imply that ‘top to bottom’ temperature gradients in the samples did not exceed ±2°C. Inner containers made of ‘iron’ (more precisely, low-carbon steel with C concentration of less than 0·6 wt%), wrapped and sealed inside platinum foil (Veksler et al., 2006Go) were used in most of the runs. However, a few runs were carried out in nickel or graphite containers. At the start of centrifuge runs, starting mixtures were first completely melted and homogenized at high, super-liquidus temperatures and the slowest rotation speed (50g) for at least 30 min, but usually for 1 h or more. The purpose of the rotation at this stage was to prevent melts from creeping up the container walls by capillary forces. Complete melting and homogenization was checked by a few test runs in which containers were opened and the glasses were examined after the heating at slow rotation. The main centrifugation stage at 1000g and experimental temperature followed after the homogenization and lasted for a few hours. Quenching was performed in air by turning off the heating, while the centrifuge rotation continued at a constant speed. During quenching, temperature inside the sample dropped from 1200 to 1100°C in 10 s, from 1100 to 1000°C in 15 s, and further cooling of the sample to 800°C took about 1·5–2 min.

Electron microprobe analyses
Run products were mounted in epoxy, ground and polished and studied by electron microprobe at the GeoForschungsZentrum Potsdam using a Cameca SX50 instrument in wavelength-dispersive spectrometry (WDS) mode at 15 kV accelerating voltage and 15 nA beam current with spot sizes ranging from 1 to 3 µm depending on the properties and compositions of the analyzed materials. Counting time for all the elements was set to 20 s on peak and 10 s on background. The following synthetic and natural standards were used for the calibration: orthoclase (Al and K), rutile (Ti), wollastonite (Si and Ca), albite and jadeite (Na), apatite (P), hematite (Fe), diopside and periclase (Mg).

Transmission electron microscopy (TEM)
TEM was carried out using a FEI TecnaiTMG2 F20 X-Twin electron microscope operated at 200 kV. Focused ion beam milling (FIB) was applied to cut site-specific TEM-foils from the area of interest. FIB preparation was performed under ultrahigh vacuum conditions in an oil-free system FEI FIB200TEM. TEM-ready foils with dimensions 20 µm x 10 µm x 0·1 µm were cut directly from microprobe epoxy mounts using a gallium-ion beam. The TEM foils were placed on a perforated carbon film on a copper grid. Carbon coating to prevent charging of the TEM samples was not applied.


    RESULTS
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Static experiments on synthetic analogues of natural immiscible glasses
The results of the static experiments are summarized in Tables 3 and 4. The most important result in the context of this study is that none of the products of the conventional static runs showed unambiguous, macroscopic signs of liquid immiscibility. The composition SF-1 that had repeatedly unmixed at almost identical T and fO2 conditions in our previous centrifuge runs (Veksler et al., 2006Go) quenched in the loops to beads of turbid, milky glass, occasionally with trace amounts of small (3–10 µm), dendritic tridymite crystals. Over-heating to 1400°C (run L-34) did not result in any visible changes in the optical appearance of the glass. TEM images of foils cut from one of the turbid glasses (Fig. 3b) revealed sub-micron heterogeneity comprising silica- and Fe-rich amorphous phases. Compositional TEM scans and profiles showed sharp co-variations in Fe, Ca and K at the sub-micron scale with areas enriched in Fe and Ca being depleted in K, but a more uniform distribution of Al and Si. The Fe-rich and Si-rich phases form 3D sponge-like interconnected textures that are usually interpreted as products of spinodal phase separation (Vogel, 1985Go; Shelby, 2005Go). Thus, it appears that although the composition SF-1 did not produce macroscopic immiscible liquid (glass) segregations in static runs, it falls within the spinodal region of the liquid miscibility gap, and consequently cannot be quenched to a homogeneous glass.


Figure 3
View larger version (91K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 3. Products of a static experiment at atmospheric pressure, 1123°C and log fO2 = –10·6 on starting composition SF-1 (sample L-6). The charge quenched to a bead of turbid glass. (a) General view of the polished epoxy mount in reflected light. The location of the fragment studied by TEM is indicated by an arrow. (b) TEM image showing sub-micron heterogeneity in the glass fragment. The high-absorption, Fe-rich amorphous phase is darker. The lightest areas are those where the low-absorption silica-rich phase goes through the whole thickness of the foil.

 

View this table:
[in this window]
[in a new window]

 
Table 3: Conditions of static runs and phase composition of the run products

 

View this table:
[in this window]
[in a new window]

 
Table 4: Electron microprobe analyses of glasses and crystal phases in products of static loop experiments

 
The composition FJ quenched at 1249°C to clear, dark brown, optically homogeneous glass. At 1203°C, it produced a similar glass with a few elongated, euhedral crystals of low-Ca pyroxene corresponding to pigeonite (sample L-9; see analyses in Table 4). As the temperature was brought further down to 1156 and 1123°C, the size and the bulk volume proportion of the pyroxene crystals increased. A strong compositional zoning developed in the crystals, with low-Ca, pigeonitic compositions in the cores, and progressively Ca-rich, augitic compositions at the rims (samples L-11 and L-5). Glass compositions retained very high FeO(t) contents of about 24–25 wt%, but no signs of liquid unmixing have been observed in the run products. The liquid clearly remained under-saturated in plagioclase down to at least 1123°C. The addition of 20 wt% of plagioclase components (mixture FJ20) also failed to stabilize plagioclase. Pigeonite remained the only liquidus phase in the samples L-14 and L-25 produced from the mixture FJ20 at 1200°C, and the mixture completely melted at and quenched to clear, optically homogeneous brownish glass from 1249 and 1280°C (samples L-17 and L-23).

Static experiments on starting mixture KC produced homogeneous, very dark brown glass at 1200 and 1249°C, and liquidus assemblages comprising low-Ca pyroxene at 1169 and 1098°C (samples L-29 and L-33). The pyroxene crystals are similar to those produced in the experiments on the mixtures FJ and FJ20, and they also showed a strong compositional zoning from pigeonitic cores to sub-calcic and high-Ca rims (Table 4). A few smaller tridymite crystals were found in sample L-33 quenched at 1098°C. Plagioclase did not crystallize, and no signs of liquid unmixing were observed in the samples. The absence of plagioclase from the products of mixtures FJ, FJ20 and KC is not consistent with the fact that the original natural glasses occur as inclusions in plagioclase, but it is not entirely unexpected in view of the low Al2O3 contents in the bulk compositions as well as quench partial liquids. As noted above, the original natural plagioclase-hosted melt inclusions may have been modified by post-entrapment diffusion and crystallization of host plagioclase on inclusion walls (see also Veksler, 2006Go).

Eight static experiments on the mixture RY40 at temperatures 1098–1319°C and log fO2 ranging from –7·8 to –12·1 all produced tridymite crystals suspended in clear, dark brown glass, and no signs of liquid immiscibility. Unlike the sample L-33 (see above), where tridymite formed tiny, micron-sized needles, and small (2–10 µm) elongated euhedral crystals, in the products of mixture RY40 it mostly formed strange-looking oval blebs or crystal aggregates measuring up to 0·4 mm in diameter. The bulk mass fraction of tridymite in the samples varied from 2% at 1319°C (sample L-30) to 13% at 1098°C (sample L-31), and tridymite showed no visible resorption at run durations of up to 35 h. Crystallization of tridymite as a single liquidus phase over such a broad temperature interval from a melt with SiO2 content of only 53·3 wt% is remarkable, and is discussed in more detail below. The absence of liquid immiscibility is no less remarkable, as the failure of melt droplets to unmix starkly contrasts with spectacular immiscibility textures of the natural iron-hosted melt inclusions, on which the mixture RY40 was based (Fig. 1).

Centrifuge experiments
Optical and electron microscope images of centrifuge run products are presented in Figs 4–8,GoGoGoGo and experimental conditions and results are summarized in Tables 5 and 6. As mentioned above, the mixture SF-1 has previously been studied and described in detail by Veksler et al. (2006Go). In this study, we present the results of centrifuge experiments on new starting compositions KC, FJ, FJ20, RY20, RY40, MZ, and MZ-2.


Figure 4
View larger version (111K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 4. Examples of centrifuged run products. (a) Sample C-117, mixture KC at 1150°C. Pyroxene crystals (cpx) in a turbid glass. (b) Sample C-106, mixture FJ and 1200°C. Large clinopyroxene crystals in a dark brown glass. (c) Sample C-114, mixture RY20 in iron container at 1150°C. Emulsion of two liquids Lfe and Lsi and suspended tridymite crystals (tr) separated at the top of the clear Fe-rich glass Lfe. (d) Sample C-116, mixture RY20 in nickel container at 1170°C. Poorly separated emulsion of two liquids and suspended tridymite crystals.

 

Figure 5
View larger version (54K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 5. Phase separation and vertical compositional gradients in products of centrifuge experiments on mixtures RY40 and RY20 analogous to the bulk compositions of the iron-hosted melt inclusions (Fig. 1). (a)–(c) represent compositional plots and a photograph of sample C-104; (d)–(f) represent analogous plots and a photograph of sample C-112. Run conditions and electron microprobe analyses of the phases are given in Tables 5 and 6. Abbreviations for phases are as in Fig. 4.

 

Figure 6
View larger version (110K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 6. Back-scattered electron image of a small area at the upper right corner of sample C-104 (see Fig. 5) showing tridymite crystals (tr) and the diffuse transition between the Fe-rich (Lfe) and silicic (Lsi) glasses.

 

Figure 7
View larger version (57K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 7. Centrifuged products of the starting mixture RY20 at 1200°C (sample C-110). (a) General view of the polished epoxy mount in reflected light. The location of the fragment studied by TEM is indicated by an arrow. (b) TEM image of the glass fragment showing sub-micron globular heterogeneity. High-absorption, Fe-rich amorphous phase is darker.

 

Figure 8
View larger version (51K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 8. Vertical compositional gradients in centrifuged products of mixtures MZ-1 and MZ-2. (a)–(c) represent compositional plots and a photograph of sample C-115; (d)–(f) represent analogous plots and a photograph of sample C-120.

 

View this table:
[in this window]
[in a new window]

 
Table 5: Run conditions and phase assemblages in products of centrifuge experiments

 

View this table:
[in this window]
[in a new window]

 
Table 6: Electron microprobe analyses of the products of centrifuge experiments

 
Our first pilot centrifuge run on starting mixture KC lasted for a total duration of 90 min (sample C-101 in Table 5), and was probably not long enough to equilibrate the charge. Olivine, which was never observed in static runs, formed at the bottom of the container, whereas tridymite crystallized at the top in a turbid, brownish glass. Glass at the bottom, just above the olivine crystals, was darker and more transparent than at the top. Vertical electron microprobe (EMP) profiles of the glass over the height of the sample revealed strong, semi-linear compositional gradients with top-to-bottom variations in SiO2 content from 72·6 to 45·5 wt%. However, when the total duration of the centrifuge runs was increased to 180 min, and initial heating to 1180°C was performed at a slower rotation speed, the mixture KC quenched to a uniformly turbid, non-transparent brownish glass, which showed no detectable compositional gradients in vertical electron microprobe profiles (see the analyses at the top and bottom of sample C-105 in Table 6). At 1150°C (sample C-117), the mixture produced abundant crystals of low-Ca pyroxene, and tridymite grains fractionated to the top of the sample in a brownish–greenish, turbid glass (Fig. 4c). Vertical electron microprobe profiles revealed significant variations in glass composition, with higher FeO contents at the bottom and more silicic compositions at the top (Table 6). Thus, we conclude that the composition KC showed some signs of unmixing (probably metastable) in the centrifuge runs, but the results were obscured by crystallization and were not reproducible at run durations shorter than 2 h.

Centrifuge experiments on mixture FJ at 1200°C produced clear, brown glasses with large, elongated crystals of low-Ca pyroxene (samples C-106 and 108). Mixture FJ20 quenched at 1200°C to turbid, greenish glass with no detectable compositional gradients in the vertical cross-section (see analyses in Table 6). In a static run (sample L-25), the mixture quenched at the same temperature to clear glass, and small amounts of low-Ca pyroxene. Thus, for mixtures FJ and FJ20 the results of centrifuged runs are generally in good agreement with static experiments at similar conditions (Tables 3 and 4). No signs of liquid immiscibility have been detected by either method.

At different run durations and 1200°C, the centrifuged products of the mixture RY40 retained, in agreement with isothermal static runs, variable amounts of tridymite, but showed also unequivocal signs of liquid immiscibility and phase separation. Longer run durations resulted in a better phase separation. However, immiscibility led to formation of LfeLsi emulsions, and at the exposures of 2 and 3 h liquid segregation was still incomplete. Vertical cross-sections of the quenched samples featured three distinct zones (Figs 5–7GoGo): (1) a bottom zone of dark brown Fe-rich glass; (2) a middle creaming zone of tridymite crystals suspended in semi-transparent, milky glass with plumes and clouds representing quenched emulsion of silica-rich and Fe-rich liquids (Fig. 7); (3) a thin (2–3 mm) layer of yellowish, honey-coloured, transparent silica-rich glass at the very top of the sample (Figs 5Go and 6). According to the TEM image of a glass foil cut from the creaming zone of the sample C-110 (Fig. 7b), the LsiLfe emulsions are characterized by drop-in-matrix morphology, which implies immiscibility in the nucleation regime (Shelby, 2005Go), and a mean droplet size of about 100–150 nm. Some droplets appear to merge by flocculation and coalescence. In addition to the silicate unmixing, a small amount of Fe–P liquid metal alloy (FeP) accumulated at the bottom of the container in most of the runs. The formation of the alloy liquid probably resulted from P2O5 reduction by C residing in the inner low-carbon steel container. The formation of the Fe–P alloy had a negligible effect on the FeO content of the silicate liquids, but resulted in a significant drop of P2O5 concentrations from 4·7 wt% in the starting charge to an average of less than 4 wt% in some of the centrifuged glasses. The increase of run duration, and the decrease in the proportion of silica-rich components (mixture RY20) visibly improved phase separation. Tridymite crystals are virtually absent from some of the products of the RY20 mixture (e.g. sample C-112 in Fig. 5), and a narrow but clear and distinct layer of silicic glass formed at the top. The compositions of the Fe-rich and silica-rich experimental glasses at 1200°C are, however, not so contrasting as those in the natural iron-hosted melt inclusions (Table 1), upon which the starting mixtures RY40 and RY20 were originally based. To increase fO2 to the level of the Ni–NiO buffer, and completely block the formation of the Fe–P alloy, we carried our a few runs in nickel containers. Two experiments on mixture RY20 (samples C-116 and C-118) produced poorly separated LfeLsi emulsions, which quenched to turbid, non-transparent glasses, and small amounts of tridymite crystals that floated to the top.

The model Skaergaard mixture MZ-1 at 1110–1120 °C was, as expected, slightly super-liquidus, and quenched in iron containers to very dark brown glass, which looked optically homogeneous. Electron microprobe profiles revealed, however, steep vertical compositional gradients for all the major oxides in a narrow layer at the very top of the samples (Fig. 8) The profiles of the components are consistent with element partitioning between Fe-rich and silica-rich immiscible liquids (e.g. Schmidt et al., 2006Go; Veksler et al., 2006Go, and references therein), and imply stratification by liquid immiscibility. As in the other ferrobasaltic compositions, liquid unmixing apparently started at a sub-micron, colloidal scale, and emulsion coalescence and creaming after 2 or 3 h of centrifugation remained incomplete. TEM line scans of foils cut at the top of the sample C-111 immediately below the thin silica-rich layer revealed fine-grained chemical heterogeneity of the glass at a sub-micron scale. Peaks in the Fe scan, with a characteristic width of 20–30 nm and a magnitude at least two times greater than background noise, corresponded to depressions in the Si scan, implying the nucleation of Fe-rich and silica-rich amorphous phases. The total proportion of silica-rich liquid Lsi in centrifuged products of the mixture MZ-1 was small, and the consolidated Lsi layer, which formed on the top, is a mere 0·2 mm thick. Longer durations apparently enhance phase separation. Experiments were also performed in nickel and graphite containers (samples C-113 and C-122). Unfortunately, both containers reacted with the melt, and changed the liquid composition. Glass in the nickel container was contaminated at 1120°C by 0·8–1·2 wt% NiO, and the charge precipitated small amounts of Ni-rich olivine and magnetite (see Tables 5 and 6). Vertical compositional gradients in the quenched glass were not so pronounced as in the iron containers, but still significant (Table 6). In the graphite container, FeO reduction resulted in crystallization of metallic iron on the container walls, and an overall drop in the FeO concentration in the melt to below 15 wt%. A slight increase in the silica content in the glass at the very top (Table 6) may indicate incipient liquid phase separation, but the gradient is even less pronounced than in the Ni container, and probably insignificant.

Centrifuge experiments on the second, more Fe-rich Skaergaard model mixture MZ-2 did not show significant compositional gradients or any other signs of liquid unmixing. In an iron container at 1110°C, melt quenched to dark brown homogeneous glass (sample C-120). The homogeneity of the glass at a sub-micron scale was confirmed by TEM study, which did not reveal any compositional heterogeneity above the background noise. In the nickel container (sample C-121), the liquid dissolved 0·8–1 wt% NiO and precipitated Ni-rich magnetite at the bottom and on the container walls. The glass showed no systematic compositional gradients, and appeared to be homogeneous to sub-micron scale.

In Fig. 9 the compositions of the glasses from the centrifuge experiments (Table 6) are plotted on the olivine–feldspar–silica projection, and compared with the position of the liquid miscibility gap in the Fe2SiO4 – KAlSi3O8 – SiO2 system (Roedder, 1951Go, 1979Go). For the projection, molar fractions of oxides (XMO) calculated from the analyses in Table 6 were combined into normative components: Fo = XMgO/2; Fa = (XFeOXTiO2)/2; Ab = 2XNa2O; Or = 2XK2O; An = XAl2O3 – 2(XNa2O + XK2O); Qz = (XSiO2 + XTiO2)/2 + 3XP2O5 – (XMgO + XFeO)/2 – 5XNa2O 5XK2OXAl2O3XCaO. Although the projection excludes important normative components such as wollastonite, apatite and ilmenite, the compositions of the immiscible Lsi and Lfe glasses from the centrifuge experiments plot well inside the immiscibility gap of the simplified synthetic system. The compositions of homogeneous glasses tend to lie to the left of the two-liquid field.


Figure 9
View larger version (16K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 9. Glass compositions from centrifuge experiments plotted on the (Fo + Fa)–(An + Ab + Or)–Qz projection. Homogeneous glasses from samples without phase separation are plotted in black circles. Immiscible glasses are shown by open circles (mixtures RY40 and RY20), and black triangles (mixture MZ-1). Pairs of the conjugate Lsi and Lfe compositions are connected by tie-lines. The gray area and dashed curves represent the miscibility gap and crystallization fields in the system Fa–Or–Qz (after Roedder, 1951Go). (See text for explanation of the projection method.)

 
Static mixing experiments
The conflicting results of static and centrifuge experiments on mixtures RY40, RY20 and SF-1, as well as the unexpected unmixing of the Skaergaard model composition MZ-1, call for a critical examination of the experimental techniques, and a careful check of equilibrium. Equilibrium is normally checked by reversal experiments, and to perform those, we prepared mixtures and synthesized glasses of the compositions MZF and MZS (Table 2), which correspond to the end-member glass compositions produced at the bottom and at the top in centrifuge runs on mixture MZ-1 and believed to represent conjugate immiscible liquids Lfe and Lsi. Because the duration of our centrifuge experiments was limited to a few hours, we carried out a series of control reversal experiments in a static setup. For that purpose, pairs of MZF and MZS melt droplets were loaded in Re loops, and Fe–Pt double containers identical to those that were used in centrifuge experiments. The idea was to see whether the liquids, which had been produced in the centrifuge, would mix back when kept in contact in a loop for a sufficient time.

In the loop set-up, approximately half of the loop was filled first with the silica-rich mixture MZS. The powdered mixture was melted for a few minutes at 1350°C and reducing conditions in a 1 atm gas mixing furnace, and quenched to glass. Then the rest of the loop was filled with the Fe-rich mixture MZF, and the charge was once again rapidly heated to 1350°C, held at that temperature for a few minutes, and quenched. Finally, the loop with two glasses was placed into gas-mixing furnace pre-stabilized at 1090 or 1114°C, and fO2 conditions indicated in Table 3, annealed for 1 or 3 days, and quenched. Run durations longer than 3 days were deemed not practical, because the losses of more volatile components, such as alkalis or P2O5, were expected to become significant, and affect bulk chemistry of the two-liquid charges.

In containers (run S-141 in Tables 3 and 4), the mixture MZF was first fused at the bottom of the inner iron container to a layer of bubble-free, homogeneous Fe-rich glass in the centrifuge furnace at 1110°C. The upper part of the container was then filled with the silica-rich mixture MZS; the outer Pt container was welded shut; the container assembly was placed in a vertical position into a static furnace, heated for 5 h at 1270°C to fuse the MZS mixture, and then kept for 72 h at 1110°C.

Despite the simple experimental design, the results appeared not as simple and straightforward as anticipated. In four loop runs L-51, 52, 56 and 58 the Lfe and Lsi liquid droplets did not mix back (Fig. 10); however, they did mix completely in the run L-57 (Table 3). Furthermore, electron microprobe analyses of quenched glasses (Table 4) revealed a significant partial mixing in all the loop runs. Notably, it mostly affected the Fe-rich liquid. Regardless of the run durations, and initial mass proportions of the two starting droplets, the Fe-rich melts dissolved roughly equal amounts of Lsi, whereas the composition of the silicic liquid (provided it did not dissolve completely as in the run L-57) changed only slightly. In all the experiments, including the run L-57, mixing stopped on the Fe-rich side at the same final melt composition characterized by a SiO2 content of 54 ± 0·5 wt%. Electron microprobe profiles across the samples demonstrated that the final Lfe glasses were chemically homogeneous, probably because of effective mixing by convection. The Lsi glasses were also homogeneous, and the transition zone between the glasses in most cases was narrow and sharp (Fig. 10b). In a single run at a lower temperature of 1090°C, groups of thin, needle-like, Ca-rich pyroxene crystals nucleated in the silica-rich glass, and the portion of glass near the crystals was further depleted in CaO, MgO, and FeO, and enriched in SiO2 by about 2 wt%. We interpret the crystals as the first liquidus phase forming in the LfeLsi mixtures.


Figure 10
View larger version (29K):
[in this window]
[in a new window]
[Download PowerPoint slide]
 
Fig. 10. Products of a static mixing run. (a) General view of sample L-56 in the Re ribbon loop showing incomplete mixing between Fe-rich glass (Lfe) and silicic glass (Lsi). (b) Electron microprobe profiles across the sample. The interface between the Fe-rich and silicic glasses is manifested by a step in the elemental profiles. (See text for further discussion.)

 
It appears that the outcome of the mixing experiments in loops depended on the mass proportion of the Fe-rich and silica-rich (MZF and MZS) liquid droplets. The masses of the MZF and MZS mixtures were carefully measured in three runs, L-56, L-57 and L-58, and the MZS/MZF mass ratios in those runs were 1·35, 0·77 and 1·2, respectively. Thus, it is not surprising that the Lsi glass completely dissolved in the run L-57 in which the initial mass proportion of the silica-rich MZS mixture was the lowest. The presence of a massive iron container also seems to play a role in the mixing process, as demonstrated by the static experiment S-141. In that experiment, the MZF and MZS liquids layers did not mix back after 3 days at 1110°C. However, microprobe profiles show that the transition between optically distinct Fe-rich and Si-rich glasses is not razor-sharp. The width of the diffusion zone between the glasses is about 0·2 mm. The diffusion zone may have formed during the brief heating to 1270°C, and it seems that no further mixing took place during the long exposure at 1110°C. The composition of the bottom Fe-rich layer remained virtually intact (Table 4). In contrast to the centrifuge experiments, minor crystallization took place in the charge. The bottom Fe-rich liquid precipitated small amounts of olivine crystals (Fo50), and a few pyroxenes crystallized in the silicic glass at the container walls, and at the MZF–MZS contact. As discussed in the next section, the iron container may function as a buffer of FeO activity and thus affect the equilibrium between two immiscible silicate liquids.

Static mixing experiment S-141 is analogous in design to some experiments performed by Longhi (1990Go) on Fe-rich and silica-rich immiscible lunar compositions. Those experiments were also performed in iron containers, and different loading methods were used. In some runs, powders of Fe-rich and silica-rich conjugate compositions were loaded in two separate layers, and this is very similar to what we did in the run S-141. Longhi (1990Go) may well have encountered the same metastability problems as we did in our static runs. He reported extensive crystallization of silicates at the interfaces between glasses of contrasting composition, and noted the absence of a clear sharp meniscus in some of his experimental products. Crystallization of pyroxene in the mixing zone between the MZF and MZS glasses in our run S-141 may be of the same nature.

The reasons why mixing repeatedly stopped (or dramatically slowed down) in loops when SiO2 concentration in Lfe reached 54 wt% are not clear. First, the composition may represent a sort of kinetic barrier at which the viscosity of Lfe becomes too high for convective mixing, forcing further mixing to proceed at a much slower rate by diffusion only. Alternatively, the composition may represent a true conjugate liquid. If so, the difference from glass compositions formed in the centrifuge may be attributed to incomplete equilibration in shorter centrifuge runs, especially for slowly diffusing network-forming components (see also Veksler et al., 2006Go). Finally, the composition may correspond to the borders of the spinodal region of the miscibility gap, where further mixing becomes energetically unfavourable. In our view, the first, kinetic explanation is less likely because the run durations were apparently long enough for effective diffusional transport in even more viscous Lsi, and no compositional gradients, which would arise from incomplete diffusion, were observed in the quenched Lfe glasses. With some reservations, we conclude that immiscibility in the MZ-1 liquid is likely to be stable, although chemical equilibration in the centrifuge was probably incomplete, and true conjugate liquids may somewhat differ in composition from the glasses that were analysed at the top and the bottom of centrifuged samples.

Crystallization experiments on the Fe-rich and silica-rich conjugate liquids
Thermodynamics requires that equilibrium immiscible liquids should crystallize identical mineral assemblages. To test this requirement of our putative immiscible liquids, we performed a series of static crystallization experiments on pure compositions MZF and MZS. Run conditions and the results are summarized in Tables 3 and 4. Plagioclase and augite are the main phases crystallizing in both mixtures at 1081 and 1062°C, and compositions of the solid solutions in the products of the MZF and MZS runs are similar (Table 4). The MZF mixture crystallized at 1081°C, plagioclase An56 and pyroxene Wo42En34Fs24 (sample L-62), whereas the mixture MZS at the same temperature crystallized zoned plagioclase An75–65 (average composition An72) and augite Wo32En42Fs26 (sample L-63). The crystalline products also include minor amounts of Fe–Ti oxides, and olivine, which appeared in products of both mixtures between 1062 and 1043°C. All the products of the MZS mixture also contain variable amounts of a pure silica phase, probably tridymite. However, this tridymite, as discussed below, may be partly or fully metastable. With falling temperature, plagioclase generally becomes, as expected, more albitic. However, in the low-temperature MZF products at 1022°C (sample L-64) plagioclase crystals showed broad compositional variations clustering around two modes (plag-1 and plag-2 in Table 4). It should be noted that equilibrium is hardly achievable at temperatures below 1050°C, especially in the silicic, viscous liquids of the mixture MZS. Despite these complications at lower temperatures, the main conclusion from the crystallization experiments is that the crystal assemblages formed from the MZF and MZS mixtures above 1050°C are similar, and so the liquids may indeed represent immiscible conjugate phases.


    DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Kinetic and thermodynamic limitations
Our results imply that reaching thermodynamic equilibrium in experiments on multicomponent ferrobasaltic compositions at temperatures below 1100°C is a very challenging task. The liquids are strongly non-ideal, and show a clear tendency towards liquid immiscibility. We believe that our results pose serious doubts regarding the reliability of the conventional static loop technique when it is applied to dry ferrobasaltic, andesitic and dacitic compositions at temperatures below 1100°C. Here, none of the unambiguous cases of immiscibility revealed in natural glasses was reproducible in static experiments.

Centrifuge experiments differed from the static runs not only in their higher acceleration, but also in a significantly greater sample mass, and the presence of massive containers. We used inner containers made of several materials. None of them proved to be perfectly inert. However, iron appears to be the best choice. The main disadvantage of iron is that it sets fO2 values somewhat below the iron–wüstite buffer [see Veksler et al. (2006Go) for further discussion], and such conditions are too reducing for the vast majority of terrestrial magmatic environments (with the rare exception of the magmas described above that precipitate native iron). Equilibration with an iron container may slightly change the FeO(t) of the silicate charge. However, on the positive side, the dissolution reaction between silicate melt and metal iron regulates not only fO2, but also the activity of FeO in the melt. Thus, the FeO component, which is crucial for liquid immiscibility, becomes effectively constrained.

Graphite containers at low pressures are even more reducing than iron and proved to be unsuitable. On the other hand, our experiments in nickel containers and in the presence of the more oxidizing Ni–NiO buffer were hampered, as described above, by Ni contamination and crystallization of Ni-bearing oxides and silicates. It is not clear, however, why liquid immiscibility and phase separation in mixtures RY20 and MZ-1 appear to be slower in nickel than in iron, as evidenced by less pronounced vertical compositional gradients in the centrifuged products. Additions of less than 1 wt% of NiO to the liquid phase(s) are not likely to cause large effects on unmixing, and a greater proportion of the Fe2O3 component in more oxidized melts is expected, if anything, to broaden the miscibility gap and enhance the unmixing (Naslund, 1983Go). However, there may be other, less obvious and indirect effects, of a kinetic nature. If liquid phase separation develops via diffusion-controlled nucleation and growth (and TEM images of quenched glasses, such as that in Fig. 7b, support this view), the process of nucleation and droplet growth may, for example, significantly slow down in liquids with a higher ferric–ferrous ratio, as amphoteric, network-forming Fe3+ is likely to diffuse much more slowly than the network-modifying Fe2+.

Metastability and liquid immiscibility
The classical, best-studied example of low-temperature silicate immiscibility in the SiO2–Al2O3–FeO–K2O system is characterized by a flat solvus emerging by maximally 170°C above the liquidus surfaces of tridymite and fayalite (Roedder, 1951Go; 1979Go; Schmidt et al., 2006Go). The stable miscibility gap along the fayalite–ferrobustamite cotectic in the Fe2SiO4 – CaFeSi2O6 – KAlSi3O8 – SiO2 system was reported to be no more than 10°C above the liquidus (Hoover & Irvine, 1978Go). Numerous other aluminosilicate systems present examples of metastable immiscibility characterized by S-shaped liquidi and liquid–liquid immiscibility domes submerged just below the liquidus surfaces (e.g. Irvine, 1976Go). Such topological relationships imply that the Gibbs free energy of two immiscible liquids and alternative phase assemblages of a single homogeneous liquid and crystal phase(s) are very close for a broad range of petrogenetically relevant basaltic and dacitic compositions. This situation should be favourable for metastability.

Tridymite as a metastable phase
Unusual crystallization of tridymite in mixtures SF-1, RY20 and RY40 has precedents in other previously published experimental studies. For example, Irvine (1976Go) reported an unexpected appearance of primary liquidus tridymite inside the plagioclase crystallization field on two pseudo-ternary joins of the Mg2SiO4 – Fe2SiO4 – CaAl2Si2O8 – KAlSi3O8 – SiO2 system. Tridymite crystallization was also reported in the mesostasis of natural tholeiitic basalt at the brink of stable liquid immiscibility (Philpotts & Doyle, 1983Go). Notably, in the latter case tridymite crystallized at temperatures below 1030–1060°C from a homogeneous liquid with only 57–59 wt% SiO2. We propose that precipitation of tridymite may in some cases represent a metastable response of silicate liquids to strong interactions between FeO and SiO2 components at ionic level (Hudon & Baker, 2002Go), which in a normal, equilibrium state would have led to liquid immiscibility. One could speculate that for complex, multicomponent liquids it may be sometimes easier, from the kinetic point of view, to decrease the instability of melt structure by precipitation of pure crystalline silica than by nucleation of a complex immiscible liquid. Our centrifuge experiments on mixtures RY20 and RY40 demonstrate that better phase separation and a greater volume of the continuous Lsi layer at longer run durations also correlate with decreasing proportion of tridymite crystals. In our view, occurrences of sporadic tridymite and quartz precipitation from melts that are high in Fe-oxides and low in SiO2 deserve a thorough re-examination from the point of view of potential liquid immiscibility.

Immiscible emulsions
Formation of colloidal emulsions is another form of metastability, which was often observed in the products of our centrifuge and static runs. Sub-micron heterogeneities are widespread in silicate glasses, and they have been usually thought to form during quench or to result from metastable liquid immiscibility (Vogel, 1985Go; Shelby, 2005Go). Plumes and clouds of colloidal droplets, as well as blurred and diffuse interfaces between immiscible liquid layers, are very common in our centrifuged run products. Although some turbid glasses in our static and centrifuge experiments may have formed by metastable unmixing during quench (e.g. samples C-107 and C-108), those centrifuged glasses that developed pronounced vertical compositional gradients probably formed from stable immiscible emulsions. The textures of quenched samples, which can be viewed as snapshots of droplet distribution after different centrifugation times, demonstrate that sub-micron droplets were slowly moving during centrifugation over macroscopic distances, and in some samples (e.g. products of the mixtures RY20 and RY40) formed dense clouds, and eventually merged into compact layers of homogeneous liquids (Figs 4c, 5 and 6). Such features cannot form in one or two minutes of quench time, and imply that those fine emulsions of silicate liquids formed by super-liquidus liquid immiscibility, and remained stable for hours. TEM images of quenched glasses allow us to see single droplets at high magnifications, measure their size, and reveal the morphology of emulsions. According to the TEM images, some turbid glasses produced in static runs may result from spinodal decomposition (Fig. 3), but droplets-in-matrix dispersed morphologies, which are typical for centrifuged glasses (Fig. 7b), imply formation and growth in a nucleation regime (Shelby, 2005Go).

The Bowen and Fenner trends: physicochemical background
Oxygen fugacity, melt non-ideality and magnetite stability
Textbook explanations of the Bowen and Fenner trends (e.g. Morse, 1994Go) are usually built upon the equation


Formula 1

(1)
when the reaction takes place between pure crystalline phases, it defines the fayalite–magnetite–quartz (FMQ) oxygen buffer, and establishes a functional dependence between T and fO2. Although natural basaltic systems do not crystallize assemblages of pure fayalite, magnetite and silica minerals but rather solid solutions of olivine and Fe–Ti oxides, with silica remaining in a residual melt, the equation is still useful as a simple, shorthand formulation of two key principles: (1) magnetite stability in magmatic systems is a function of fO2; (2) crystallization of magnetite should lead to a strong increase in the SiO2 content of the liquid.

The link between magnetite crystallization and the Bowen (silica enrichment) trend as a result of fundamental mass-balance principles is undeniable. However, two things should be kept in mind when equation (1) is applied to complex natural basaltic systems. First, fO2 is surely not the only factor controlling the magnetite stability in multicomponent magmas. In absence of liquidus tridymite or quartz, silica activity (aSiO2) is not fixed at unity, and becomes a variable, which must be raised to the third power and incorporated in the equilibrium constant for reaction (2). A detailed analysis of the relationship between silica activity and oxygen fugacity in layered gabbroic intrusions has been carried out by Morse (1980Go, 1990Go). Because of the strong non-ideality of silicate liquids, aSiO2 is a complex function of melt composition. The ferrous–ferric ratio in basaltic liquids certainly depends on aSiO2 and the activities of other components, especially alkalis (e.g. Kilinc et al., 1983Go). In view of the important role of melt chemistry, it is hardly surprising that alkalis have been shown to have a strong effect on magnetite stability, and the general trend of melt evolution in petrologically relevant aluminosilicate systems (Irvine, 1976Go; Shi, 1993Go).

The second consideration is also related to the strong non-ideality of ferrobasaltic liquids. Some experimental models show that magnetite crystallization does not necessarily prevent the formation of Fenner-type (Fe-rich) liquids. Furthermore, in those cases where non-ideality takes the extreme form of stable immiscibility, Fe-rich liquids coexist with silicic, Bowen-type liquids, and various assemblages of Fe–Mg silicates, Fe–Ti oxides, and silica minerals. Such relationships have been most vividly demonstrated in simplified synthetic systems where liquidus phases comprise pure fayalite and magnetite, and aSiO2 is fixed at unity by liquidus tridymite (e.g. Naslund, 1983Go). Naslund (1983Go) showed that fayalite–tridymite and magnetite–tridymite cotectics in the system Na2O(K2O)–FeO–Fe2O3–Al2O3–SiO2 intersected a region of stable liquid immiscibility over a broad range of fO2 from 10–12 to 0·2. As a result, a series of Bowen- and Fenner-type aluminosilicate liquids were shown to coexist regardless of the direction of the magnetite–fayalite reaction (2). This work also showed that the miscibility gap greatly expanded, and the compositions of coexisting liquids diverged at oxidizing conditions and in equilibrium with magnetite. Apparently, this means that the Fe2O3 component causes greater non-ideality in aluminosilicate melts than FeO. The main conclusion, however, is that the effect of magnetite crystallization upon the direction of melt evolution may be secondary and subordinate to strong chemical interactions between components and species in the liquid itself. Strong intrinsic non-ideality of ferrobasaltic liquids may be a more important factor in magma evolution than the composition of liquidus mineral assemblages or fO2 conditions imposed on the system.

Concentration–activity relationships in complex, multicomponent liquids
Unlike the immiscible synthetic systems referred to above, experiments on natural basaltic compositions have seldom produced stable liquid unmixing. Above, we have already expressed our doubts regarding the reliability of conventional static techniques in the vicinity of or inside stable miscibility gaps. Here we would like to point out that non-ideality arising from anomalously strong interactions between FeO and SiO2 components in silicate melts [see Hudon & Baker (2002Go) for a theoretical discussion] should be universal and reveal itself also in liquid compositions lying outside miscibility gaps. In the case that non-ideality does not result in stable unmixing or metastable phase separation, it reveals itself in broad compositional variations of activity coefficients. Activities in the liquid phase can be assessed from crystal–liquid equilibria. A few examples illustrating this approach are presented in Table 7. Here we compare our samples L-69 and L-70, the products of the MZF and MZS starting mixtures, with two samples published earlier by McBirney & Naslund (1990Go) and Toplis & Carroll (1995Go). The data show that very different liquid compositions with SiO2 concentrations from 46 to 72 wt% can be equilibrated at similar T and fO2 with almost identical crystal phase assemblages comprising olivine (Fo27–35), plagioclase (An43–48), clinopyroxene, and Fe–Ti oxides. Notably, the liquids appear to show a strong negative correlation between the SiO2 and FeO components, and positive correlations between FeO and CaO, SiO2 and Al2O3, and SiO2 and alkalis. Also notable is that the spectrum of liquids from silicic to extremely Fe-rich forms in a narrow range of redox conditions, and the presence or absence of magnetite does not seem to be a crucial factor. Importantly, Fenner-type liquids at the Fe-rich side of the spectrum are characterized not only by low silica content, but also by low alkalis and high Ca/Al values. These chemical characteristics appear to arise from activity–composition relationships in complex aluminosilicate liquids where FeO, on the one hand, increases the activity coefficient of SiO2, whereas alkalis decrease it (e.g. Shi, 1993Go). A very similar pattern of the activity–composition relationships is revealed by element partitioning between the Lsi and Lfe immiscible liquids [e.g. Tables 1and 6; see also Veksler et al. (2006Go) and references therein]. In the process of the formation of immiscible liquids, alkalis partition to the Bowen-type silicic liquids, whereas FeO and alkaline earths partition to the Fenner-type liquids, which are also often characterized by very low alumina contents.


View this table:
[in this window]
[in a new window]

 
Table 7: Liquid compositions and crystal assemblages produced by melting of natural Skaergaard gabbros (sample UB-4145) and crystallization of synthetic glasses

 
It was also demonstrated that variations of silica content in melts along the pseudo-binary join AnDi–SiO2 (where AnDi is the eutectic composition in the anorthite–diopside system) have complex, non-linear effects on Fe solubility (Borisov, 2007Go). The maximum of Fe solubility (that is, a minimum activity coefficient of FeO) was observed in melts with about 55–57 wt% SiO2, regardless of temperature and oxygen fugacity. Melts on both sides of the maximum show significantly lower Fe solubility and, correspondingly, higher FeO activity coefficients. It is clear that Lsi melts should lie on the silica-rich side of the maximum, where increasing silica contents result in decreasing FeO concentrations.

Implications for the Skaergaard intrusion and tholeiitic magmas in general
The most recent major dispute about the Fenner and Bowen trends in the Skaergaard intrusion took place in the late 1980s. It was started by Hunter & Sparks (1987Go), who argued against the strong Fe-oxide enrichment in the Skaergaard liquid on the basis of mass-balance considerations. A series of papers followed (Brooks & Nielsen, 1990Go; Hunter & Sparks, 1990Go; McBirney & Naslund, 1990Go; Morse, 1990Go), in which the controversy was examined from different angles, and on the basis of different types of evidence. In the context of the present discussion, we would like to emphasize the thermodynamic and experimental aspects of the debate, and point out that liquid evolution paths at the high- and low-SiO2 ends of Fig. 2 are directly based on experimental studies. Ferrobasaltic glasses with FeO(t) well above 26 wt% were produced in a series of partial melting experiments on natural Skaergaard gabbros (McBirney & Nakamura, 1974Go; McBirney & Naslund, 1990Go). On the other hand, numerous crystallization experiments on synthetic glasses and rock powders, which modeled the starting composition of the Skaergaard magma, produced liquids with FeO(t) concentrations of no more than 18–22 wt% (e.g. Toplis & Carroll, 1995Go; Thy et al., 2006Go), near the compositional path inflections where Fe–Ti oxides start to crystallize, after which experimental liquids evolve towards silica enrichment. This apparent discrepancy demands an explanation.

Our static experiments on compositions MZF and MZS (Tables 3, 4 and 7) confirm that very different liquid compositions from ferrobasaltic to rhyolitic can be equilibrated under identical conditions with crystal assemblages typical for cumulates of the Middle Zone (MZ) and Upper Zone (UZ) of Skaergaard. It is also true that very similar mineral assemblages can be found in felsic volcanic rocks (Hunter & Sparks, 1987Go, 1990Go; Sparks, 1988Go) characterized by rhyolitic groundmass, bulk-rock SiO2 concentrations well above 60 wt% and low FeO(t). Because the compositions MZF and MZS were produced by immiscibility in centrifuge experiments we propose that the diverging liquid evolution paths in Fig. 2 skirt a region of stable liquid immiscibility. Alternatively, immiscibility may be metastable, and the paths may fan out across a flat and broad liquidus plateau. Flat liquidus surfaces are typical for regions above metastable immiscibility domes (e.g. Irvine, 1976Go). In any case, strong non-ideality rooted in the anomalous incompatibility between Fe-oxide and SiO2 components in silicate melt structures appears to be the principal reason for the divergent liquid evolution paths.

Silicate liquid immiscibility has been mentioned but previously never as a central issue in debates about the Bowen and Fenner trends (e.g. McBirney, 1975Go). Immiscibility certainly operated in the Upper Zone of the Skaergaard Layered Series, close to the end of magma crystallization. It has been reproduced in melting experiments on the UZ cumulates (McBirney & Nakamura, 1974Go; McBirney & Naslund, 1990Go), confirmed by the discovery of silicic and Fe-rich melt inclusions in cumulus apatite from UZ layered gabbros (Jakobsen et al., 2005Go; see also Fig. 2) and indicated also by the common occurrence of large-scale, outcrop-sized segregations of melanogranophyre. Our results imply that the process may have started earlier, and demand that the role of silicate liquid immiscibility should be reconsidered.

We propose that the Fe-enrichment trend in Skaergaard was driven at first by fractional crystallization of troctolitic and gabbroic cumulates that formed the lower parts of the Layered Series (the LZa and LZb of the conventional Skaergaerd stratigraphy, e.g. McBirney, 1989Go). This is a classical tholeiitic trend, primarily caused by a high cotectic proportion of plagioclase in the cumulus assemblages, along the lines explained, for example, by Grove & Baker (1984Go). Such Fe-enrichment continued until about 60 wt% of crystallization (Nielsen, 2004Go), which corresponds to crystallization of the upper part of the Lower Zone (LZc), when the magma became saturated in Fe–Ti oxides. Any further significant Fe-enrichment by fractional crystallization would have been impossible, and the liquid would have started to evolve towards silica enrichment. However, in view of our centrifuge experiments on the composition MZ-1, liquid immiscibility may have started in the Skaergaard magma at about the same time. The initial proportion of the immiscible silica-rich liquid Lsi in the Fe-rich Skaergaard magma was probably small, at around 15 wt%, but the unmixing and gravitational separation of the conjugate liquids provided a mechanism for continuous Fe enrichment in the denser liquid at the bottom of the magma chamber despite the crystallization of Fe–Ti oxides. In our currently preferred scenario, fractional crystallization coupled with liquid immiscibility continued in Skaergraard until the end of crystallization in the Sandwich Horizon. We believe that the combined effects of crystallization and unmixing may explain many enigmatic features of the Skaergaard rocks and other layered gabbroic intrusions. Studies of plagioclase-hosted melt inclusions (Jakobsen et al., 2007Go) and pyroxene–plagioclase symplectites (Stripp et al., 2007Go) have already produced circumstantial evidence in favour of early liquid immiscibility in the Skaergaard intrusion. However, a detailed discussion of the numerous geological implications of this hypothesis is beyond the scope of this paper and will be presented elsewhere.

Immiscible emulsions of silicate melts, which were routinely formed in our experiments, are of special interest for future research. Kinetic and thermodynamic factors regulating the formation and stability of such emulsions in complex aluminosilicate melts are at the moment poorly understood. It is not clear whether such emulsions can be sufficiently stable at the timescale of natural magmatic processes to seriously affect magma dynamics. As discussed above, our results hint at significant effects of melt composition upon the kinetics of unmixing (see also Veksler et al., 2007Go). Viscosity, interfacial energy and diffusion are likely to be the key factors, and diffusion-driven coarsening of silicate emulsions is expected to be slow.

In general, it is envisaged that the proposed immiscibility between Bowen- and Fenner-type liquids is most likely, and should be most effective in large, slowly cooling basaltic magma chambers, where the proportion of the Fe-rich, low-viscosity liquid is large; where there is enough time for emulsion coarsening as a result of coalescence and/or Ostwald ripening; and where other conditions are met for large-scale gravitational phase separation. The Skaergaard intrusion may represent an example of such a magma chamber. Perhaps the long-standing Bowen–Fenner controversy will be finally resolved at the location where it has been debated most, the Skaergaard intrusion in East Greenland.


    ACKNOWLEDGEMENTS
 
This study was inspired by fruitful collaboration and discussions of Skaergaard with Alexey Ariskin, Charles Kent Brooks, Marian Holness, Jakob Jakobsen, Christian Tegner, Troels Nielsen, and Gemma Stripp. We are grateful to Dieter Rhede and Oona Appelt for the help with microprobe analyses. Reviews by Peter Thy, Tony Morse and two anonymous reviewers were thoughtful, friendly, and most helpful. I.V.V. has been supported by grants from the German Science Foundation (DFG).


*Corresponding author. E-mail: veksler{at}gfz-potsdam.de


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 SILICATE LIQUID IMMISCIBILITY IN...
 EXPERIMENTAL AND ANALYTICAL...
 RESULTS
 DISCUSSION
 REFERENCES
 
Ariskin AA. Phase equilibria modelling in igneous petrology: use of COMAGMAT model for simulating fractionation of ferrobasaltic magmas and the genesis of high-alumina basalts. Journal of Volcanology and Geothermal Research (1999) 90:115–162.[CrossRef][Web of Science]

Ariskin AA. Geochemical thermometry of layered series rocks of the Skaergaard intrusion. Petrology (2002) 10:495–518.[Web of Science]

Bird JM, Goodrich CA, Weathers MS. Petrogenesis of Uifaq iron, Disko island, Greenland. Journal of Geophysical Research (1981) 86:11787–11805.[CrossRef]

Bogaerts M, Schmidt MW. 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 (2006) 152:257–274.[CrossRef][Web of Science]

Borisov AA. Experimental study of the effect of SiO2 on cobalt and iron solubilities in silicate melts. Petrology (2007) 15. (in press).

Borisov A, Jones JH. An evaluation of Re, as an alternative to Pt, for 1-bar loop technique: an experimental study at 1400°C. American Mineralogist (1999) 84:1528–1534.[Abstract]

Bowen NL. The Evolution of the Igneous Rocks (1928) Princeton, NJ: Princeton University Press.

Brooks CK, Nielsen TFD. Early stages in the differentiation of the Skaergaard magma as revealed by a closely related suite of dike rocks. Lithos (1978) 11:1–14.[CrossRef][Web of Science]

Brooks CK, Nielsen TFD. The differentiation of the Skaergaard intrusion. A discussion of Hunter and Sparks (Contrib Mineral Petrol 95: 451–461). Contributions to Mineralogy and Petrology (1990) 104:244–247.[CrossRef][Web of Science]

Carmichael ISE. The petrology of Thingmuli, a Tertiary volcano in eastern Iceland. Journal of Petrology (1964) 5:435–460.[Abstract/Free Full Text]

Daly RA. Igneous Rocks and their Origin (1914) New York: McGraw–Hill.

De A. Silicate liquid immiscibility in the Deccan Traps and its petrogenetic significance. Geological Society of America Bulletin (1974) 85:471–474.[Abstract/Free Full Text]

Dixon S, Rutherford MJ. Plagiogranites as late-stage immiscible liquids in ophiolite and mid-ocean ridge suites: an experimental study. Earth and Planetary Science Letters (1979) 45:45–60.[CrossRef][Web of Science]

Dorfman AM, Hess K.-U, Dingwell DB. Centrifuge-assisted falling-sphere viscosimetry. European Journal of Mineralogy (1996) 8:507–514.[Abstract/Free Full Text]

Fenner CN. The crystallization of basalts. American Journal of Science (1929) 18:225–253.[Web of Science]

Fujii T, Kushiro I, Nakamura Y, Koyaguchi T. A note on silicate liquid immiscibility in Japanese volcanic rocks. Journal of the Geological Society of Japan (1980) 86:409–412.

Greig JW. Immiscibility in silicate melts. American Journal of Science (1927) 13:133–154.[Web of Science]

Grove TL, Baker MB. Phase equilibrium controls on the tholeiitic versus calc-alkaline differentiation trends. Journal of Geophysical Research (1984) 89(B5):3253–3274.

Hanghøj K, Rosing M, Brooks CK. Evolution of the Skærgaard magma: evidence from crystallized melt inclusions. Contributions to Mineralogy and Petrology (1995) 120:265–269.[Web of Science]

Hoover JD, Irvine TN. Liquidus relations and Mg–Fe partitioning on part of the system Mg2SiO4–Fe2SiO4–CaMgSi2O6–CaFeSi2O6 –KAlSi3O8–SiO2. Carnegie Institution of Washington Yearbook (1978) 77:774–784.

Hudon P, Baker DR. The nature of phase separation in binary oxide melts and glasses. I. Silicate systems. Journal of Non-Crystalline Solids (2002) 303:299–345.[CrossRef]

Hunter RH, Sparks RSJ. The differentiation of the Skaergaard Intrusion. Contributions to Mineralogy and Petrology (1987) 95:451–461.[CrossRef][Web of Science]

Hunter RH, Sparks RSJ. The differentiation of the Skaergaard Intrusion. A reply. Contributions to Mineralogy and Petrology (1990) 104:248–254.[CrossRef][Web of Science]

Irvine TN. Metastable liquid immiscibility and MgO–FeO–SiO2 fractionation patterns in the system Mg2SiO4–Fe2SiO4–CaAl2Si2O8–KAlSi3O8–SiO2. Carnegie Institution of Washington Yearbook (1976) 75:597–611.

Jakobsen JK, Veksler IV, Tegner C, Brooks CK. Immiscible iron- and silica-rich melts in basalt petrogenesis documented in the Skaergaard intrusion. Geology (2005) 33:885–888.[Abstract/Free Full Text]

Jakobsen JK, Veksler IV, Tegner C, Brooks CK. Crystallization of the Skaergaard Intrusion from an emulsion of immiscible iron- and silica-rich liquids: Constraints from the melt inclusions in plagioclase. In: Frontiers in Mineral Sciences 2007, Programme and Abstracts—Harte B, Carpenter M, eds. (2007) Cambridge: The Mineralogical Society of Great Britain and Ireland. 230–231.

Jang YD, Naslund HR, McBirney AR. The differentiation trend of the Skaergaard intrusion and the timing of magnetite crystallization: Iron enrichment revisited. Earth and Planetary Science Letters (2001) 189:189–196.[CrossRef][Web of Science]

Kilinc A, Carmichael ISE, Rivers ML, Sack RO. The ferrous–ferric ratio in natural silicate liquids equilibrated in air. Contributions to Mineralogy and Petrology (1983) 83:136–140.[CrossRef][Web of Science]

Krasov NF, Clocchiatti R. Immiscibility in silicate melts and its possible petrogenetic importance, as shown by study of melt inclusions. Transactions (Doklady) of the USSR Academy of Sciences (1979) 248:92–95.

Longhi J. Silicate liquid immiscibility in isothermal crystallization experiments. Proceedings of the 20th Lunar and Planetary Science Conference—Sharpton VL, Ryder G, eds. (1990) New York: Pergamon. 13–24.

McBirney AR. Differentiation of the Skaergaard intrusion. Nature (1975) 253:691–694.[CrossRef]

McBirney AR. The Skaergaard layered series: I. Structure and average compositions. Journal of Petrology (1989) 30:363–397.[Abstract/Free Full Text]

McBirney AR, Nakamura Y. Immiscibility in late-stage magmas of the Skaergaard intrusion. Carnegie Institution of Washington Yearbook (1974) 73:348–352.

McBirney AR, Naslund HR. The differentiation of the Skaergaard intrusion. A discussion of Hunter and Sparks (Contrib Mineral Petrol 95: 451–461). Contributions to Mineralogy and Petrology (1990) 104:235–240.[CrossRef][Web of Science]

Morse SA. Kiglapait mineralogy II: Fe–Ti oxide minerals and the activities of oxygen and silica. Journal of Petrology (1980) 21:685–719.[Abstract/Free Full Text]

Morse SA. The differentiation of the Skaergaard intrusion. A discussion of Hunter and Sparks (Contrib Mineral Petrol 95: 451–461). Contributions to Mineralogy and Petrology (1990) 104:240–244.[CrossRef][Web of Science]

Morse SA. Basalts and Phase Diagrams (1994) Melbourne, FL: Krieger.

Naslund HR. The effect of oxygen fugacity on liquid immiscibility in iron-bearing silicate melts. American Journal of Science (1983) 283:1034–1059.[Abstract/Free Full Text]

Nielsen TFD. The shape and volume of the Skaergaard intrusion, Greenland: implications for mass balance and bulk composition. Journal of Petrology (2004) 45:507–530.[Abstract/Free Full Text]

Osborn EF. The reaction principle. In: The Evolution of Igneous Rocks. Fiftieth Anniversary Perspectives—Yoder H. S. Jr, ed. (1979) Princeton, NJ: Princeton University Press. 133–169.

Philpotts AR. Silicate liquid immiscibility: its probable extent and petrogenetic significance. American Journal of Science (1976) 276:1147–1177.[Abstract/Free Full Text]

Philpotts AR. Silicate liquid immiscibility in tholeiitic basalts. Journal of Petrology (1979) 20:99–118.[Abstract/Free Full Text]

Philpotts AR. Liquid immiscibility in silicate melt inclusions in plagioclase phenocrysts. Bulletin de Minéralogie (1981) 104:317–324.

Philpotts AR. Compositions of immiscible liquids in volcanic rocks. Contributions to Mineralogy and Petrology (1982) 80:201–218.[CrossRef][Web of Science]

Philpotts AR, Doyle CD. Effects of magma oxidation state on the extent of silicate liquid immiscibility in a tholeiitic basalt. American Journal of Science (1983) 283:967–985.[Abstract/Free Full Text]

Roedder E. Low-temperature liquid immiscibility in the system K2O–FeO–Al2O3–SiO2. American Mineralogist (1951) 36:282–286.[Web of Science]

Roedder E. Silicate liquid immiscibility in magmas. In: The Evolution of Igneous Rocks. Fiftieth Anniversary Perspectives—Yoder H. S. Jr, ed. (1979) Princeton, NJ: Princeton University Press. 15–58.

Roedder E, Weiblen PW. Lunar petrology of silicate melt inclusions, Apollo 11 rocks. In:. Proceedings of the Apollo 11 Lunar Science Conference. Geochimica et Cosmochimica Acta, Supplement 1 (1970) 1:507–528.

Roedder E, Weiblen PW. Petrology of silicate melt inclusions, Apollo 11 and Apollo 12 and terrestrial equivalents. In:. Proceedings of the 2nd Lunar Science Conference. Geochimica et Cosmochimica Acta, Supplement 2 (1971) 1:507–528.

Ryabov VV. Chemical composition of immiscible liquids in natural glasses from traps. Soviet Geology and Geophysics (1988) 29(11):62–70.

Ryabov VV. Liquation in Natural Glasses: the Example of Traps. (1989) Novosibirsk: Nauka (in Russian).

Schmidt MW, Connolly JAD, Günther D, Bogaerts M. Element partitioning—the role of melt structure and composition. Science (2006) 312:1646–1650.[Abstract/Free Full Text]

Shearer CK, Papike JJ, Spilde MN. Trace element partitioning between immiscible lunar melts: An example from naturally occurring lunar melt inclusions. American Mineralogist (2001) 86:238–241.[Abstract/Free Full Text]

Shelby JE. edn. In: Introduction to Glass Science and Technology (2005) 2nd. Cambridge: The Royal Society of Chemistry.

Shi P. Low-pressure phase relationships in the system Na2O–CaO–FeO–MgO–Al2O3–SiO2 at 1100°C, with implications for the differentiation of basaltic magmas. Journal of Petrology (1993) 34:743–762.[Abstract/Free Full Text]

Sparks RSJ. Petrology and geochemistry of the Loch Ba ring-dyke, Mull (N. W. Scotland): an example of extreme differentiation of tholeiitic magma. Contributions to Mineralogy and Petrology (1988) 100:446–461.[CrossRef][Web of Science]

Stripp G, Holness M, Nielsen T, Veksler I, Tegner C. Reactive symplectite formation as tracers of late-stage liquids: the Skaergaard Intrusion. In: Frontiers in Mineral Sciences 2007, Programme and Abstracts—Harte B, Carpenter M, eds. (2007) Cambridge: The Mineralogical Society of Great Britain and Ireland. 227.

Tegner C. Iron in plagioclase as a monitor of the differentiation of the Skaergaard intrusion. Contributions to Mineralogy and Petrology (1997) 128:45–51.[CrossRef][Web of Science]

Toplis MJ, Carroll MR. 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 (1995) 36:1137–1171.[Abstract/Free Full Text]

Thy P, Lesher CE, Nielsen TFD, Brooks CK. Experimental constraints on the Skaergaard liquid line of descent. Lithos (2006) 92:154–180.[CrossRef][Web of Science]

Veksler IV. Crystallized melt inclusions in gabbroic rocks. Melt Inclusions in Plutonic Rocks. Mineralogical Association of Canada Short Courses—Webster JD, ed. (2006) 36:100–122.

Veksler IV, Dorfman AM, Danyushevsky LM, Jakobsen JK, Dingwell DB. Immiscible silicate liquid partition coefficients: implications for crystal–melt element partitioning and basalt petrogenesis. Contributions to Mineralogy and Petrology (2006) 152:685–702.[CrossRef][Web of Science]

Veksler IV, Dorfman AM, Wirth R, Dingwell DB. Unmixing of ferrobasaltic liquids: a study of nucleation, settling and coalescence of emulsion droplets in centrifuged glasses. In: Frontiers in Mineral Sciences 2007, Programme and Abstracts—Harte B, Carpenter M, eds. (2007) Cambridge: The Mineralogical Society of Great Britain and Ireland. 213.

Vogel W. Chemistry of Glass (1985) Columbus, OH: American Ceramic Society.

Wager LR, Brown GM. Layered Igneous Rocks (1968) Edinburgh: Oliver & Boyd.

Wager LR, Deer WA. Geological investigation in East Greenland, Part III, The petrology of the Skaergaard intrusion, Kangerdluqssuaq, East Greenland. Meddelser om Grønland (1939) 105((4)).

Yoder HS, Jr, Tilley CE. Origin of basaltic magmas: an experimental study of natural and synthetic rock systems. Journal of Petrology (1962) 3:342–532.[Abstract/Free Full Text]


Add to CiteULike CiteULike   Add to Connotea Connotea   Add to Del.icio.us Del.icio.us    What's this?


This article has been cited by other articles:


Home page
J PetrologyHome page
R. G. Cawthorn and L. D. Ashwal
Origin of Anorthosite and Magnetitite Layers in the Bushveld Complex, Constrained by Major Element Compositions of Plagioclase
J. Petrology, September 1, 2009; 50(9): 1607 - 1637.
[Abstract] [Full Text] [PDF]


Home page
Eur J MineralHome page
T. Martin Crespo, R. P. Lozano Fernandez, and R. Gonzalez Laguna
The fulgurite of Torre de Moncorvo (Portugal): description and analysis of the glass
European Journal of Mineralogy, August 1, 2009; 21(4): 783 - 794.
[Abstract] [Full Text] [PDF]


Home page
J PetrologyHome page
A. R. Philpotts
Comments on: Liquid Immiscibility and the Evolution of Basaltic Magma
J. Petrology, December 9, 2008; (2008) egn061v1.
[Full Text] [PDF]


Home page
J PetrologyHome page
A. R. McBirney
Comments on: 'Liquid Immiscibility and the Evolution of Basaltic Magma' Journal of Petrology 48, 2187-2210
J. Petrology, December 9, 2008; (2008) egn062v1.
[Full Text] [PDF]


Home page
J PetrologyHome page
S. A. Morse
Compositional Convection Trumps Silicate Liquid Immisciblity in Layered Intrusions: a Discussion of 'Liquid immiscibility and the evolution of basaltic magma' by Veksler et al., Journal of Petrology 48, 2187-2210
J. Petrology, December 9, 2008; (2008) egn063v1.
[Full Text] [PDF]


Home page
J PetrologyHome page
I. V. Veksler, A. M. Dorfman, A. A. Borisov, R. Wirth, and D. B. Dingwell
Liquid Immiscibility and Evolution of Basaltic Magma: Reply to S. A. Morse, A. R. McBirney and A. R. Philpotts
J. Petrology, December 9, 2008; (2008) egn064v1.
[Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow FREE Full Text (PDF) Freely available
Right arrow All Versions of this Article:
48/11/2187    most recent
egm056v1
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Email this article to a friend
Right arrow Similar articles in this journal
Right arrow Alert me to new issues of the journal
Right arrow Add to My Personal Archive
Right arrow Download to citation manager
Right arrowRequest Permissions
Google Scholar
Right arrow Articles by Veksler, I. V.
Right arrow Articles by Dingwell, D. B.
Right arrow Search for Related Content
GeoRef
Right arrow GeoRef Citation
Social Bookmarking
 Add to CiteULike   Add to Connotea   Add to Del.icio.us  
What's this?