Journal of Petrology Advance Access published online on July 24, 2008
Journal of Petrology, doi:10.1093/petrology/egn036
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
The Roles of Fractional Crystallization, Magma Mixing, Crystal Mush Remobilization and Volatile–Melt Interactions in the Genesis of a Young Basalt–Peralkaline Rhyolite Suite, the Greater Olkaria Volcanic Complex, Kenya Rift Valley
1Environment Centre, Lancaster University, Lancaster LA1 4YQ, UK
2Institute of Geochemistry, Mineralogy and Petrology, University of Warsaw, AL.
wirki i Wigury 93, 02-089 Warsaw, Poland
3US Geological Survey, 956 National Center, Reston, VA 20192, USA
4School of Geosciences, University of Edinburgh, West Mains Road, Edinburgh EH9 3JW, UK
5Department of Earth Sciences, Cespar, the Open University, Walton Hall, Milton Keynes MK7 6AA, UK
Received March 13, 2008; Revised typescript accepted July 2, 2008
| ABSTRACT |
|---|
The Greater Olkaria Volcanic Complex is a young (
20 ka) multi-centred lava and dome field dominated by the eruption of peralkaline rhyolites. Basaltic and trachytic magmas have been erupted peripherally to the complex and also form, with mugearites and benmoreites, an extensive suite of magmatic inclusions in the rhyolites. The eruptive rocks commonly represent mixed magmas and the magmatic inclusions are themselves two-, three- or four-component mixes. All rock types may carry xenocrysts of alkali feldspar, and less commonly plagioclase, derived from magma mixing and by remobilization of crystal mushes and/or plutonic rocks. Xenoliths in the range gabbro–syenite are common in the lavas and magmatic inclusions, the more salic varieties sometimes containing silicic glass representing partial melts and ranging in composition from anorthite ± corundum- to acmite-normative. The peralkaline varieties are broadly similar, in major element terms, to the eruptive peralkaline rhyolites. The basalt–trachyte suite formed by a combination of fractional crystallization, magma mixing and resorption of earlier-formed crystals. Matrix glass in metaluminous trachytes has a peralkaline rhyolitic composition, indicating that the eruptive rhyolites may have formed by fractional crystallization of trachyte. Anomalous trace element enrichments (e.g.
2000 ppm Y in a benmoreite) and negative Ce anomalies may have resulted from various Na- and K-enriched fluids evolving from melts of intermediate composition and either being lost from the system or enriched in other parts of the reservoirs. A small group of nepheline-normative, usually peralkaline, magmatic inclusions was formed by fluid transfer between peralkaline rhyolitic and benmoreitic magmas. The plumbing system of the complex consists of several independent reservoirs and conduits, repeatedly recharged by batches of mafic magma, with ubiquitous magma mixing. KEY WORDS: alkali basalt; geochemistry; magma mixing; trace element
| INTRODUCTION |
|---|
Many recent studies have demonstrated the complexity of silicic magma chambers, a complexity related, inter alia, to magma composition, the pressure–temperature conditions of formation and differentiation, the local tectonics and variable residence times within the upper crust. In addition to crystal fractionation, combined assimilation–fractional crystallization and magma mixing, many systems have evolved in part by the remobilization of crystal mushes and/or partly to completely solidified plutonic bodies that may or may not be cognate with their eruptive hosts (Nakada et al., 1994
The majority of recent studies have focused on metaluminous systems, in both continental and subduction-related settings. None has attempted to place the evolution of a young peralkaline rhyolite system into a time–composition framework, although the differences in bulk composition, with consequent effects on phase equilibria, density, viscosity and volatile behaviour, might be expected to generate evolutionary histories different from the metaluminous cases. Furthermore, most studies of silicic systems have focused on relatively large, commonly caldera-forming, centres. In this paper, we report on the Greater Olkaria Volcanic Complex (GOVC; Fig. 1), a young peralkaline rhyolite dome and lava field in the south–central Kenya rift valley that has been marked over the past
20 kyr by the development of at least 80 small centres from which
15 km3 of magma have been erupted in hundreds of eruptions. The eruptive rate can be contrasted with, for example, the highly productive Okataina centre in New Zealand, where 80 km3 of rhyolite magma were erupted in nine episodes during the last 26 kyr (Nairn et al., 2004
).
|
An additional issue concerns the origin of peralkaline rhyolite magmas and the long-standing debate as to whether silicic magmas form by extended fractional crystallization of mafic parents or by crustal anatexis. Previous studies of the GOVC used trace element and isotopic evidence to suggest that the rhyolites formed by volatile-induced melting of heterogeneous crustal rocks (Davies & Macdonald, 1987
In this paper, we address the complex temporal, structural and magmatic evolution of the GOVC. We show that trachyte and rhyolite magmas formed by a combination of fractional crystallization of basalts, ubiquitous magma mixing, assimilation of earlier crystal mushes and/or plutonic material, and interactions between melts and volatile phases. Marshall et al. (in prep.) have addressed the further histories of the rhyolite magmas and their significance for the development of the higher structural levels of the GOVC.
The specific aims of this paper are:
- to provide mineral chemical and geochemical data on the basalt and trachyte eruptive rocks and on the magmatic inclusions in the peralkaline rhyolites;
- to determine the roles of fractional crystallization, magma mixing and remobilization of pre-existing solid material in the evolution of the basalt–trachyte suite;
- to re-examine the petrogenesis of the peralkaline rhyolitic magmas;
- to place constraints on models of the GOVC plumbing system.
| GEOLOGICAL EVOLUTION OF THE GOVC |
|---|
The GOVC is part of the central Kenya peralkaline province (CKPP; Macdonald & Scaillet, 2006
|
The earliest phase (stage 3) of post-caldera magmatism, the Lower Comendite Member of the Olkaria Comendite Formation, is represented by rhyolite lavas (O2) and pyroclastics (Op2), dated by 14C at >9150 ± 110 BP. Stage 4, the Middle Comendite Member, was a period of ring dome formation (O3) and eruption of thick surge deposits (Op3). 14C dates bracket its age as between 9150 ± 110 and 3280 ± 120 BP. This was followed (stage 5) by a resurgence of the caldera floor and the formation of short, thick flows, the Upper Comendite Member (O4). The sixth stage saw the eruption of very thick flows of comendite from a north–south fissure system (O5), the youngest flow having yielded a 14C date of 180 ± 50 BP (Clarke et al., 1990
Mafic eruptive rocks occur both north and south of the dome field (Fig. 2). To the south, the Lolonito Basalt Formation (Ba1) forms a string of small outcrops and comprises flows and pyroclastic cones with eruption ages of <0·45 Ma. The Formation may represent an early phase of the generally younger Akira Basalt Formation (Ba2), which consists of lavas and scoria and spatter cones (Clarke et al., 1990
). The youngest cone is unvegetated and ungullied and may be only a few hundred years old. The formation is geographically associated with, and of the same general age as, the Tandamara Trachyte Formation (Tt), composed of welded scoriaceous lapilli and blocks and some lavas (or rheomorphic ignimbrites).
To the north of the dome field, mafic volcanism accompanied rhyolite emplacement on the Ndabibi plains, an 11 km wide, low-lying area between the GOVC and the Eburru Volcanic Complex (Fig. 2). Faulted lavas and pyroclastic cones of pre-caldera age comprise the Ndabibi Basalt Formation (Bn; Clarke et al., 1990
). Younger units include surtseyan tuff cones (kBt) with laterally equivalent fall tuffs (kBtm) and transitional surtseyan–strombolian ash cones (Bt). Judging from their degree of dissection, the latter units are broadly coeval with the post-caldera rhyolites (i.e. with ages of <20 ka BP).
Magmatic inclusions, ranging in size mainly from a few millimetres to 10 cm, but rarely up to
40 cm, occur in rhyolites from all stages of post-caldera activity. They are most easily studied in faulted and/or eroded outcrops, particularly in the Ol Njorowa erosion channel. Where present, they usually make up <1% of the outcrop but rarely can be 5–10%. Magmatic inclusions also occur in composite intrusions. For example, the 50 m high Central Tower plug of O4 age (Fig. 2) has an
20 cm thick skin of trachyte, which was still semi-molten when intruded by a core of comendite containing magmatic inclusions ranging in bulk composition from mugearite to benmoreite.
Two kinds of inclusion are broadly distinguishable in the field: a dark, sparsely porphyritic, more mafic variety, and purple–brown, more felsic types. The dark inclusions tend to be small (from <1 mm to
10 cm), vary in form from sub-rounded to lensoid, have irregular, crenulated margins and low porosities (
3%). Rare, highly vesiculated (50–60%) inclusions vary in shape but invariably have crenulated or streaky margins. Some inclusions have trails of mafic material incorporated in the comendite host, and patches and re-entrants of vesicular rhyolite occur within the inclusions. Some have pillow-like forms, sometimes with crenulate and/or chilled margins, and with inter-pillow spaces filled by rhyolite. Chilled margins are not always present, perhaps because many inclusions were small enough to chill en masse. Other samples show mingled zones and/or zones of vesiculation at the contacts between inclusion and host. The lighter-coloured inclusions are generally larger, up to
40 cm, and sub-rounded to slightly elongate. The margins are crenulate but less irregularly shaped than in the dark inclusions.
The textural evidence indicates that the inclusions were incorporated into the host rhyolites in a dominantly liquid state (Bacon & Metz, 1984
). Swirling of comendites and inclusions suggest mixing in eruptive conduits. More detailed textural evidence (below) shows that the inclusions themselves are very commonly mixed magma rocks.
The distribution of surface outcrops of the basalts and rhyolites indicates that the GOVC dome field displays a shadow zone, where mafic magma is prevented by its density from penetrating the overlying silicic magma, except as components of mixed rocks. The question arises as to what constitutes the Olkaria magmatic system, particularly as the Eburru Volcanic Complex and the Longonot caldera volcano lie immediately to the north and east, respectively, of the GOVC (Fig. 1). We include in the system the Ndabibi basalts, because basalts and comendites have co-erupted and mixed there, and the Akira basalts and Tandamara trachytes, because geochemically they are closely similar to the Ndabibi basalts and Olkaria trachytes.
| SAMPLING STRATEGY AND ANALYTICAL METHODS |
|---|
|
|
|---|
Our sampling covered, with the exception of the Maiella Pumice, the major units described above; namely, the Olkaria and Tandamara Trachyte Formations, the Lolonito, Akira and Ndabibi Basalt Formations, selected samples from the Olkaria Comendite Formation, and magmatic inclusions in the post-caldera comendites. The four available analyses of Maiella Pumice (Clarke et al., 1990
|
The first record of magmatic inclusions at Olkaria was by Bone (1987
Electron microprobe analysis (EMPA) of phenocryst and xenocryst phases in 18 samples, selected to cover the range of petrographic features in the suite, was carried out at the US Geological Survey, Reston, VA, and EMPA of a further set of phenocrysts, xenocrysts and glasses was carried out in the Department of Earth Sciences, The Open University. Details of the techniques employed in both laboratories are given in Supplementary Appendix 1 (available for downloading from the Journal of Petrology web site, at http://petrology.oxfordjournals.org). Representative analyses of phenocryst phases are given in Tables 2–6![]()
![]()
![]()
. The full dataset, including structural formulae and details of the analytical laboratory, is available as supplementary tables on the Journal of Petrology web site; these tables are referred to subsequently as ST1–8. Eighteen analyses of trace element abundances (Table 3; ST2) were made using laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) at the Open University, using techniques outlined in Supplementary Appendix 1.
|
|
|
|
|
Major element compositions of 59 samples were determined by X-ray fluorescence (XRF) on fused discs at Lancaster University. Results for standard AGV1 are given in ST7. Representative analyses are given in Table 7 and the full dataset in ST7. Wet-chemical analyses of a further three samples were taken from Bliss (1979
|
|
| PETROGRAPHY |
|---|
Phenocryst assemblages are given in the Appendix and modal analyses of 21 representative rocks in Table 1. We divide the rocks into three groups. Lacking petrographic information for the BGS samples, we have assigned them to groups on the basis of their chemical composition (see below). A main suite (Group 1) ranges in bulk composition from basalt through hawaiite, mugearite, benmoreite and trachyte to rhyolite. The group is represented by extrusive rocks and magmatic inclusions and includes various kinds of mixed magma rocks. Group 2 comprises rocks with bulk compositions in the hawaiite–benmoreite range, but compositionally different from Group 1, formed by the mixing of basaltic and trachytic or rhyolitic magmas. This group is dominated by extrusive rocks but includes a few magmatic inclusions. Five silica-undersaturated magmatic inclusions constitute Group 3. All three groups contain an extensive suite of feldspar xenocrysts and gabbroic to syenitic xenoliths in various states of disaggregation.
Basalts and hawaiites
The mafic rocks (Lolonito and Ndabibi Basalt Formations, mafic components of mixed rocks) are plagioclase-phyric, with subordinate olivine and clinopyroxene phenocrysts (Table 1; Appendix). There are two petrographic types, finely (more sparsely) and coarsely porphyritic, which do not differ significantly in chemical composition. The two types occasionally occur in the same flow. The less-phyric rocks (e.g. BL64) contain <5% plagioclase phenocrysts
2 mm long; these tend to be euhedral, although some internal resorption is common. Plagioclase crystals in more coarsely porphyritic varieties (e.g. BL25, BL158a) typically form 10–20% modally, and are elongate laths or blocky crystals up to 5 mm long. They are weakly zoned and some show poorly developed sieved interiors with partially resorbed rims. Olivine and clinopyroxene phenocrysts in both varieties tend to be small (<1 mm) and subhedral to euhedral. Glomeroporphyritic clusters of all three phenocryst phases are very common.
Mugearites and benmoreites
Phenocrysts (plagioclase + clinopyroxene + FeTi-oxides ±olivine) in the mugearites and benmoreites tend to be smaller and less abundant than in the basaltic rocks; some rocks are aphyric (e.g. BB84.24, BB84.39; Appendix). The plagioclase most commonly forms laths up to 0·5 mm long, some showing swallowtail terminations and, occasionally, hollow cores; these textures indicate rapid growth from the magma (Lofgren, 1980
). Olivines are scarce, colourless, small (<0·15 mm) and usually associated with clinopyroxene. The clinopyroxene phenocrysts are typically small (
0·2 mm), colourless and slightly rounded; they tend to occur in clusters, sometimes with plagioclase. The FeTi-oxides are equant, small (
0·2 mm) and also tend to cluster together.
Trachytes
Where fresh, the trachytes are massive and grey but some are microvesicular and some have a brownish colour as a result of oxidation. They contain up to 32% phenocrysts (Table 1), the assemblage being alkali feldspar + olivine + clinopyroxene + FeTi-oxide + apatite ±zircon ± chevkinite ± sulphide. The dominant phase is blocky alkali feldspar, up to 7 mm across, usually retaining euhedral or subhedral form but almost invariably strongly resorbed internally. Many crystals are simply zoned. Some rocks (e.g. SMN180) contain a generation of almost totally resorbed feldspars. Olivine phenocrysts are most commonly present as colourless prisms, up to 1·5 mm long, showing varying degrees of resorption. The clinopyroxene crystals are grey–green, small (< 1·5 mm) and fairly euhedral, although normally showing incipient resorption. FeTi-oxides are usually equant in form and < 1 mm, rarely up to 2 mm, in size; they can be discrete or included in clinopyroxene or alkali feldspar. Apatite microphenocrysts sometimes occur as discrete crystals but more commonly are associated with, or included in, clinopyroxene, FeTi-oxide, olivine or alkali feldspar. Crystals are prismatic, up to 300 µm long, frequently show simple or oscillatory zonation (with zones
10 µm across), and may have a central hollow. Zircon and chevkinite microphenocrysts are relatively scarce, typically only a few grains per thin section, and occur either as discrete crystals or in association with alkali feldspar, clinopyroxene and FeTi-oxide. Crystals are elongated, up to 350 µm long, and either euhedral or with a slightly rounded form.
Mixed magma rocks
We distinguish two types of mixed magma rocks; mingled, where two or more melt components have not completely blended, and hybridized, where the mixing process has produced a well-blended magma, the mixing stage being recorded mainly in disequilibrium phenocryst assemblages. Mixed magma rocks are found in all three groups as defined above. Mixing takes two forms in Group 1 : mixing of basaltic magmas and mixing of mugearites, benmoreites, trachytes and rhyolites in various proportions. The dominant expression of basalt–basalt mixing is a streaky admixture of melts of slightly different composition and/or texture. For example, the matrix in BB85.250 has two intimately intermixed components: one is slightly coarser and carries plagioclase microphenocrysts; the other is finer-grained, lacks plagioclase microphenocrysts, is rather more clinopyroxene-rich and tends to form spheres in the coarser component. A less common form of basalt–basalt mixing is exemplified in BL146a (Fig. 3a), where rounded inclusions of a titanaugite-bearing, oxide-rich basalt, up to 8 mm across and forming 12% modally, are embedded in a more typical olivine basalt. A high Ti content in pyroxene typically indicates high rates of crystallization (Feeley & Dungan, 1996
).
|
The second type of mixing in Group 1, found also in Group 3, occurs in the magmatic inclusions. They contain one or more components in the mugearite–rhyolite compositional range in various degrees of mingling and hybridization. Some representative examples are shown in Fig. 3. The mafic material has commonly chilled to pillows, 0·5–1 mm across; in some cases, the comendite host has filled the interstices between pillows and is often fluorite- and riebeckite-rich (Fig. 3b). Vesiculation of host and/or inclusion magmas at their contact is common (Fig. 3b). In Fig. 3c, the mafic magma has dispersed into rounded to crenulated blebs within the comendite host (light brown). This sample (SMN105) contains 22% of mafic inclusions but most samples have significantly smaller amounts (e.g. 3·4% in SMN168 and 1·3% in BB85.26). In some cases, as shown in Fig. 3d, the mafic magma has fingered into the salic host; the slightly greyer material on the right may be a third component or a mix of the other two. There is sometimes a zone of intermediate composition between the rhyolite host and the inclusion, formed either by physical commingling or diffusive mixing of the melts (Fig. 3g and h).
The petrographic and textural relationships in Group 1 clearly indicate that the inclusions range from single-component (unmixed) melts to, more commonly, mixtures of two, and more rarely three or four, magmatic components. In BB85.269N, for example, the components range from mugearitic to trachytic. The inclusions were then intruded into, and locally interacted with, the comendite melts. The extrusive trachytes show little evidence of mingling but it is clear from phenocryst compositions (below) that at least three flows (SMN179, SMN180, BB85.18) are hybrids.
Group 2 comprises rocks from the Ndabibi, Lolonito and Akira Basalt Formations and two magmatic inclusions (Appendix). They are mingled (BB85.38, SMN103, BB85.94, BL229, BL469, SMN158, BL136) or hybridized (BL120c, BB85.42, 860130, BB85.93, BL74, BL563a, BL574) mixtures of mafic and trachytic or rhyolitic (comendite and pantellerite) magmas. The two components normally form streaky or patchy mixtures at a scale <2 mm.
Perhaps the most notable feature of the Group 3 rocks, all magmatic inclusions, is the intimate character of the interaction between the inclusions and their rhyolitic hosts. In BB84.39 and BB85.185, the mafic component quenched into micro-spherules < 1 mm in diameter, enclosed by, and sometimes veined by, salic material, often rich in fluorite and sodic amphibole. In SMN153, the mafic inclusion was breaking up in, and was veined by, the host. The two components in SMN121 are streaked out and rolled together. Mixing at the junction between inclusion and host was common.
In summary, mafic, intermediate and salic magmas have been present contemporaneously in the Olkaria system for at least the past 20 kyr and have very frequently mixed at various stages during ascent and eruption.
Xenocrysts and xenoliths
A significant proportion (59%) of the mafic–intermediate extrusive rocks and magmatic inclusions carry large (up to 6 mm) grains of alkali feldspar, and less frequently plagioclase, forming up to 9% modally (Table 1; Appendix). Such feldpars invariably have a jacket of more mafic material against the rhyolitic host (Fig. 3g and h), a feature reported for large xenocrysts of alkali feldspar in intermediate enclaves at the Gedemsa volcano in Ethiopia by Peccerillo et al. (2003
). The feldspars are generally intensely embayed and resorbed, to the extent that they can be difficult at a cursory glance to distinguish from the host matrix (Fig. 3e). The patchy replacement of plagioclase by sanidine and the growth of calcic dendritic rims on sanidine are almost certainly a response to magma mixing (Hibbard, 1981
). The growth of dendritic plagioclase has been reproduced experimentally and is highly dependent on the degree of undercooling of the mafic melt (Lofgren, 1974
). Textural and compositional relationships are complex, depending, as we shall see below, on the number and proportions of mixing components.
Coarse-grained xenoliths ranging in composition from gabbro (SMN105; Fig. 3c) through syenodiorite–syenite (SMN166, BB85.40, SMN153, BB85.278) to syenite (SMN105, Fig. 3g; BB85.268) commonly form strings up to 7 mm long and have almost invariably been partially disaggregated by resorption. In some cases, the alkali feldspars show dendritic rims (Fig. 3g and h).
Glass-bearing feldspars and feldspar aggregates occur in samples SMN103, BL225a and SMN168I. In SMN103, the mugearitic host magma has infiltrated into a plagioclase crystal (An32·0Ab61·7Or6·4) (Table 2) along cleavage planes and has largely replaced the inner parts of the crystal. The residual feldspar is penetrated by glass (Fig. 4). A 5 mm x 3 mm xenolith (no. 3) in BL225a comprises residual alkali feldspar crystals slightly reverse zoned from An0·1Ab59·2Or40·7 to An1·7Ab60·9Or37·4 (Table 2) in a matrix with two components (Fig. 5): a clear glass showing perlitic cracks, and sub-rounded areas up to 300 µm across containing vermicular quartz and platy feldspar in grey glass. There was some mingling of the two textural varieties of melt along their contacts.
|
|
Xenolith 1 in BL225a is 10 mm x 6 mm and comprises alkali feldspar (An3·3–0·2Ab60–64Or33–40) in glass that is dominantly pale brown, with much less abundant, arcuate patches of deep brown glass. Xenolith 2 (6 mm x 3 mm) is similar except that the deep brown glass usually contains irregular aggregates of FeTi-oxides (Fig. 6). SMN168I carries a 5 mm long cluster of alkali feldspars that have been partially melted along their rims. The residual feldspars have the average composition An0·2Ab56·7Or43·1 and there are two types of glass, one deep brown, the other pale brown to colourless and containing minute oxide filaments. Following Bachmann et al. (2002
|
| MINERAL CHEMISTRY |
|---|
Feldspars
Phenocrysts
Plagioclase phenocrysts in the mafic rocks, up to 17% modally (Table 1), form subhedral to euhedral plates up to 4 mm x 3 mm, although incipient resorption features are common. Core compositions are in the range An89–38 (Fig. 7; Table 2; ST1a) and the extent and type of zoning (up to 8% An content) vary between rocks.
|
In the mugearites and benmoreites, core compositions are in the range An57–30 and zoning is modest (An
8). Alkali feldspar phenocrysts in the trachytes are mainly sanidine, but many crystals have anorthoclase cores (Table 2; ST1b). In SMN180, for example, the range of core compositions is An19–2. Both normal and reverse zoning occur, usually modest (An <4). There are several trends within the sanidine crystals; for example, the displacement of the trends for SMN180, SMN199 and SMN168 in Fig. 7a (inset) should be noted. The shift is probably related to compositional differences in the parental magmas.
Xenocrysts
Compositions are presented in ST1c. We exemplify the complexity using one rock, although we stress that all the mafic and intermediate rocks are to some degree unique in terms of their xenocryst content and composition. We have attempted to distinguish hybrid phenocryst assemblages from higher-temperature residuals by the presence of reverse zoning in the former but we note that such zoning can also be formed by reheating without mixing.
In the mixed magma rock BB85.40 (Fig. 7b), plagioclase phenocrysts (An52·2, An46·6), comprising euhedral laths some 80 µm x 200 µm, show modest reverse zoning (An <5). Phenocrysts of composition An40–30, forming partly resorbed, platy crystals
300 µm x 200 µm in size, show modest normal zoning (An <7). Xenocryst 1 is reversely zoned from An8·2 to An21·5. Xenocryst 2 (not shown) has darker and lighter core zones [on back-scattered electron (BSE) images] of An38·7 and An8·7, respectively, and a dendritic rim of An25·7. Xenocrysts 4 and 6 are patchily zoned between An27 and An21 and xenocryst 5 is normally zoned from An33·8 to An26. Overall, it seems that sodic plagioclase and sanidine phenocrysts from mugearitic and trachytic magmas, respectively, were attempting to re-equilibrate to
An25 after mixing, the composition in equilibrium with the new melt composition. The reverse zoning from andesine to labradorite requires a different explanation. It may have been caused by the host basalt mixing with, or being heated by, a more magnesian magma, or by loss of a volatile phase before mixing with the trachyte.
We have also made several electron microprobe analyses of a large (5 mm) crystal in BB85.40 (Z in Fig. 7b; Fig. 8), which shows significant resorption by the host magma. Although the textural relationships are complex, we distinguish an outer calcic zone (>An20), an interior zone patchily replaced by An8, and bright zones on the BSE image which tend to have the highest Ba contents in the crystal (BaO
5·61 wt%; Cn 11%). The substitution scheme in this crystal, and in all the xenocrysts, is Ba + Al
K + Si.
|
The compositional range in this crystal covers much of the range in the phenocrysts (Fig. 7a). Our interpretation is that an alkali feldspar xenocryst was incorporated into trachyte magma and was being intensely, patchily replaced before its host was admixed with a more mafic magma (benmoreite?) and the dendritic rim was formed.
A comparable situation is found in a 4 mm crystal in SMN111 (Fig. 9): a partially resorbed, patchily replaced feldspar (An21–16) has a narrow dendritic rim (An26). Our interpretation is that the crystal was incorporated into a melt (the lighter zone in the image) more mafic than that from which it crystallized and a more calcic, dendritic rim was formed (Hibbard, 1981
), followed by intense resorption. Before resorption was completed, the crystal plus surrounding melt were mixed into the host rhyolite (darker zone).
|
Anhedral, partially resorbed plagioclase xenocrysts, patchily zoned and up to 2 mm across, have been found in eight rocks. They occur singly or in glomerophyric clusters, often involving alkali feldspar. We have analyses from two rocks, both indicating reverse zoning, in BB85.268 from An42·7 to An46·1 and in BB85.278 from An35·3 to An43·3.
The question remains as to the source of the large xenocrysts. In many cases, the alkali feldspars were derived from a trachytic magma; their size and composition match the phenocrysts in the trachytes. However, the plagioclase and ternary feldspars, too calcic to have crystallized from trachyte, may have been part of coarse-grained, solidified rocks, such as the syenodiorites and syenites found in BB85.278, or from crystal mushes. Such recycled crystals have been termed antecrysts by W. Hildreth (see Bacon & Lowenstern, 2005
).
Trace elements in feldspars
Barium and Sr contents of feldspar phenocrysts and xenocrysts have been determined by EMPA (Table 2; ST1). In addition, Ba, Cs, Ga, Li, Pb, Rb, Sr and the rare earth elements (REE) have been determined by LA-ICP-MS in feldspar phenocrysts from the mafic portion of a mixed rock, a trachyte, a mixed dyke rock and a magmatic inclusion (Table 3; ST2).
Barium contents in the feldspar phenocrysts and xenocrysts show a peaked distribution against An content. In the more calcic (An >30) and least calcic (An <5) feldpars, Ba is generally low (<0.l0 a.p.f.u.). Between An30 and An5, and peaking at An10, many are notably Ba-rich, the maxima being 5·61 and 8·78 wt% BaO determined by EMPA and LA-ICP-MS, respectively, in BB85.40 (Tables 2 and 3). Barium contents thus peaked in trachytic melts somewhat less evolved than the eruptive trachytes. In contrast, all the feldspars are Sr-poor (SrO
0·16 wt%). Sr contents generally decrease, and Rb contents increase, with decreasing An content.
Chondrite-normalized REE patterns for feldspars are shown in Fig. 10. The plagioclase phenocryst (An87) in the basalt is light REE (LREE)-enriched relative to the heavy REE (HREE), which are all below detection, and shows a positive Eu anomaly (Eu/Eu* = 12). Feldspars in the mixed rocks have compositions in the range An36–20. The REE patterns are broadly intermediate between those for the basalt and trachyte but with very large, positive Eu anomalies (Eu/Eu* = 183–1137).
|
Alkali feldspar in the trachyte shows a steeper decrease towards the middle REE (MREE) ([La/Sm]CN = 72, as opposed to 5 in the basalt) but a fairly flat pattern between the MREE and HREE. The Eu anomaly is smaller than in the ternary feldspars (Eu/Eu* = 20). The maximum Eu anomalies are, therefore, generally found in the more Ba-rich feldspars. REE abundances and chondrite-normalized patterns vary both within crystals (e.g. Ce from 4 to 8 ppm in crystal 1, SMN199) and between crystals in the same rock (e.g. 4–11 ppm in crystal 1 and 11–16 ppm in crystal 6, SMN199), which may have been related to element diffusion during reheating and mixing events.
Olivine
Olivine phenocryst cores in the mafic rocks are in the range Fo88–65 (Table 4; ST3). Zoning, where present, is normal (Fo <6). Core compositions in the mugearites are more Fe-rich (Fo65–60). There are only three data points between Fo60 and Fo20 and none between Fo34 and Fo21, consistent with the disappearance of olivine from the phenocryst assemblage in the benmoreites. Core compositions in the trachytes range from Fo21·0–17·9 in SMN179 to Fo<1 in SMN168. MnO and CaO concentrations reach 2·71 and 0·57 wt%, respectively, and generally show negative correlations with mg-number. Concentrations of NiO are very low (
0·11 wt%), even in the most magnesian rocks.
In Fig. 11, we assess whether the phenocrysts were in equilibrium with their respective whole-rock compositions, using the technique of Reubi & Nicholls (2005
).
|
The distribution coefficients (Fe/MgKDmin/liq) used were: basalt 0·25–34; mugearite–benmoreite 0·29–0·40; and trachyte/rhyolite 0·40–0·69; and it was assumed that Fe2+ = 0·9 Fetotal in the melt. Olivines in the basalts and trachytes have compositions close to equilibrium compositions; cores of
Fo81 in SMN103 and BB85.42 (zoned to a rim of Fo57) probably represent relicts from a higher-temperature stage of evolution. In the mixed rocks, the phenocrysts have clearly been derived from one or more of the mix components.
Clinopyroxene
The compositional range of phenocryst cores is fairly continuous from augite through ferroaugite to ferrohedenbergite (En50–0), although there are relatively few analyses around En30 (Fig. 12; Table 5; ST4). Concentrations of Ti and Al are moderate and the two elements show a strong positive correlation. Maximum abundance levels (Ti
0·1 and Al
0·33 a.p.f.u., respectively) occur in pyroxenes with mg-number
55. Sodium abundances are low, even in the most Fe-rich pyroxenes (<0·15 a.p.f.u. =
5% acmite molecule).
|
In Fig. 13, we assess whether the phenocrysts were in equilibrium with their respective whole-rock compositions, using the technique of Reubi & Nicholls (2005
|
As for the olivines, phenocrysts in the mafic rocks, with core compositions in the range En50–40, approach equilibrium compositions. The mixed rocks and trachytes usually contain two or more phenocryst populations, the variable sense of zoning strongly suggesting that the clinopyroxenes from two mixing magmas were attempting to re-equilibrate with the new host composition (Fig. 12, inset).
FeTi-oxides
FeTi-oxides join the phenocryst assemblage in the benmoreites and occur as discrete crystals up to 1 mm long and as inclusions in fayalite, clinopyroxene and alkali feldspar. Both rhombohedral and spinel phases occur mainly as subhedral to heavily resorbed tabular grains, although the rhombohedral phase also forms rods and prisms. Ilmenite occasionally forms thin plates within spinels. Spinels sometimes show exsolution features and both phases can contain sulphide inclusions.
The rhombohedral phase is ilmenite (Xilm 0·97–0·91; Table 6; ST5). Magnesium contents can be moderately high (
0·09 a.p.f.u.) but Mn contents are invariably low (
0·03 a.p.f.u.). Niobium abundances range up to 0·02 a.p.f.u. and are variable in crystals in the same rock (e.g. 0·003–0·009 a.p.f.u. in SMN179). No compositional parameter varies systematically with whole-rock composition.
The spinels are compositionally more variable than the ilmenites (Xusp 0·55–0·74; Table 6; ST5). Niobium concentrations are low (Nb <0·005 p.f.u.) but Al abundances reach 0·12 a.p.f.u. and Mg contents 0·1 a.p.f.u.
Apatite and chevkinite
We will report on the GOVC apatites elsewhere. Here we note that the major coupled substitution is REE3+ + Si4+ = Ca2+ + P5+ and that there is significant solid substitution towards britholite-(Ce). The maximum britholite component, in apatites from SMN168, is 20%. The most REE-rich phases occur in the most evolved rocks (e.g. SMN168 and SMN199).
Chevkinite-(Ce) has been analysed in two trachytes and the salic component of the mixed magma rocks SMN168 and SMN105 (ST6). They conform closely to the general formula for the group; namely, A4BC2D2Si4O22, where A = (La, Ce, Ca, Sr, Th), B = Fe2+, C = (Fe2+, Fe3+, Ti, Al, Nb, Zr) and D = Ti. There is some variation in the proportions of Ce, La and Ca in the A site, although all are chevkinite-(Ce), and of Fe, Ti and Nb in the C site. Thus, Ce/La varies from 1·52 to 1·86 and is positively correlated with Fe/Ti, which ranges from 0·58 to 1·74. There is variable LREE enrichment relative to the MREE and HREE; La/Sm and La/Y range from 19 to 56 and La/Y from 10 to 63; both ratios tend to be higher in the chevkinites with lower Ce/La and Fe/Ti ratios.
| GEOTHERMOBAROMETRY |
|---|
Coexisting FeTi-oxides are common in the trachyte lavas (e.g. SMN179, SMN199, BB85.18) and in the salic components of mixed rocks (e.g. BB85.268, SMN168, BB85.40, SMN105). We have analytical data (Table 6; ST5) for three oxide pairs from trachytes; crystal sizes range from 100 to 700 µm and crystals have sharp contacts against glass. None of the pairs is in mutual contact but those in BB85-18 are attached to the same olivine crystal. All three pairs pass the Bacon & Hirschmann (1988
FMQ + 0·5 to –1·6. The higher temperatures for the Olkaria rocks may reflect their less peralkaline nature, lower volatile contents or reheating of the magmas by more mafic material. The latter process would be consistent with the strongly resorbed nature of the alkali feldspar phenocrysts in the trachytes. There would probably be no time constraints on the oxides recording the temperature increase; Pichavant et al. (2007| GEOCHEMISTRY |
|---|
Group 1
With the proviso that in many cases the terms refer to the bulk composition of mixed magma rocks, the main suite rocks (Group 1) vary from basalts through hawaiites, mugearites, benmoreites and trachytes to rhyolites in the TAS classification scheme (Fig. 14; Le Bas et al., 1986
2; hy
11); all the more evolved rocks are hy ± q-normative. The trachytes are mainly metaluminous, although five samples are mildly peralkaline. We have been unable to discern any systematic differences in composition between the mafic rocks of the Lolonito, Ndabibi and Akira Formations; accordingly, they are grouped together in subsequent discussion.
|
On a MgO–FeO* plot, Group 1 rocks show a peaked trend, with maximum FeO* (total Fe as FeO) abundances of
13 wt % at 3–4 wt % MgO (Fig. 15). In an attempt to identify the likely liquid line(s) of descent, the inset shows rocks for which there is no petrographic evidence of magma mixing (Appendix) or, with the exception of Ba, significant contamination by feldspar xenocrysts. Group 1 as a whole mirrors the trend shown by these rocks and we infer that the compositional ranges at given MgO values reflect either the presence of several liquid lines of descent and/or additional processes causing scatter.
|
MgO plots for major elements and selected trace elements are shown in Fig. 16a and b. In each case, the liquid line of descent, as defined above, is also shown. Overall, there are, with decreasing MgO abundance, increases in SiO2, Na2O, K2O, Ga, Rb, Nb, Ta, LREE, Mo, Sn, Zr, Hf and Th, and decreases in CaO, Ni, Cr, Co, V and Sc. Abundances of FeO*, TiO2, P2O5, Li and Sb peak at
3–4 wt % MgO, and of Ba, Cs and Y at 1–2 wt % MgO. For many elements, there is significant scatter between 3 and 1 wt % MgO. Incompatible trace element (ITE) ratios (e.g. Zr/Nb, Zr/Y, Ce/Y, Ce/Zr, Rb/Zr and Th/U) tend to be constant in the range basalt–benmoreite and then to show much greater ranges in the trachytes.
|
Part of the scatter at low MgO values is due to the fact at these low abundances, MgO is not an efficient index of differentiation. Using, as an alternative, the Thornton–Tuttle Differentiation Index (DI), there are, with increasing DI, overall increases in SiO2, Al2O3, K2O, Zr, Nb and the LREE, and decreases in TiO2, FeO*, MgO, CaO, P2O5, Ba, Sr, Y and Na2O/K2O (Fig. 17). Zr/Nb, Zr/Y and Ce/Y ratios initially rise gently but then increase more sharply at DI
85. There are no systematic compositional differences between the Olkaria trachytes, Tandamara trachytes and the trachytic melt inclusions, although the melt inclusions tend to be slightly less evolved (average DI = 79) than the Tandamara trachytes (86) and the Olkaria trachytes (87). The trend includes metaluminous and peralkaline types and we can infer that the trachytic magmas evolved along several lineages close to alumina-undersaturation. Sample SMN118 is anomalous among the trachytes in being ol-normative, having low MgO and TiO2 contents and in being aphyric. It appears to represent a separate lineage.
|
An important difference in REE distribution between the extrusive rocks and the magmatic inclusions of Group 1 is shown on a Ce–Y plot (Fig. 18). There are two trends in rocks with Ce >100 ppm: with increasing Ce, the majority of extrusive varieties show increasing Ce/Y, whereas the magmatic inclusions show a scattered trend of increasing Y and decreasing Ce/Y. The extrusive trend is essentially that defined by the inferred liquid line of descent and we note below that the Ce enrichment factor along the inferred liquid line of descent (
10) is comparable with that of Zr (
12). Thus rocks with low Ce/Y ratios are Y-enriched, not Ce-depleted. We now explore this further using chondrite-normalized REE patterns.
|
Davies & Macdonald (1987
|
The patterns in the magmatic inclusions are more diverse (Fig. 19). The rather flat patterns ([La/Yb]CN = 2·8–8·4) are a result of strong HREE enrichment (see Fig. 18) and they also show strong negative Eu anomalies (0·55–0·12), the patterns resembling those in the more evolved GOVC comendites (Macdonald et al., 1987
As well as the high abundances in the HREE noted above, some Group 1 rocks show anomalous enrichments in other elements. For example, BB85.01 and BB85.55 have high Zn, Cu and U contents (Zn 3553 and 3842 ppm; Cu 225 and 257 ppm; U 60·0 and 62·7 ppm; Table 8), BB85.01 and SMN121 have high Pb contents (233 and 243 ppm, respectively). The high-Y rocks (SMN171, BB85.01 and SMN118) are different from each other; for example, SMN171 has high Fe, Ti, Ba and Sr, whereas BB85.01 has normal abundances of those elements but the high Zn, Cu, U and Pb abundances noted above.
Groups 2 and 3
The location of the Group 2 rocks in the TAS (Fig. 14) and MgO–FeO* plots (Fig. 15), and all other compositional features, are consistent with their being mixtures of mafic and trachytic/rhyolitic magmas, although in many cases the presence of feldspar xenocrysts makes it a three-component calculation. In the straightforward case of the mingled rock BL229, which has neither unusually high Ba contents nor petrographic evidence of feldspar xenocrysts, the end-member proportions are 72% basalt (with 6·5 wt % MgO) and 28% pantellerite (
0·0 wt % MgO). The total range of mafic:salic components is from 85 to 33. Sample BB85.93 has a positive Eu anomaly (Eu/Eu* = 1·19).
Relative to Group 1 rocks of similar MgO content, the ne-normative rocks of Group 3 generally have higher abundances of Na2O and, with the exception of SMN121, high Rb, Th and the REE. Three of the five rocks are peralkaline, with Peralkalinity Indices [PI; mol. (Na2O + K2O)/Al2O3] up to 1·39; the other two (SMN153 and BB85.185) have PI = 0·94 and 0·95 (i.e. are close to peralkaline). Magmatic inclusions SMN153 and SMN169ii show modest LREE enrichment (LaCN 624, 717; [La/Yb]CN 9·3, 8·2) and negative Eu anomalies (Eu/Eu* 0·54, 0·39). SMN121, also a magmatic inclusion, is fairly similar (LaCN 418; [La/Yb]CN 9·4) but uniquely in our dataset shows a positive Ce anomaly (Ce/Ce* = 1·4).
Glass in partially melted xenoliths and feldspars
We have analyses of the glass in partially melted xenoliths and feldspars in SMN103, BL225a and SMN168I (Table 9; ST7). The glasses cover the range trachyte to rhyolite (SiO2 62·4–74·6 wt %) and vary from an ± c- to ac-normative. In SMN168I, the deep brown glass is higher in SiO2, FeO* and TiO2 and lower in Al2O3, CaO, Na2O, BaO, an and ab than the pale brown to colourless glass, whereas K2O is about constant. All are an ± c-normative but the PI increases with increasing SiO2. These trends are consistent with a partial melting sequence in which the earliest melts were enriched in quartz and femic components and later melts increasingly rich in feldspar. Xenolith 3 in BL225a contains both ac- and an-normative glass, the former types being richer in SiO2, TiO2, FeO*, MnO and CaO, and poorer in Al2O3, Na2O and K2O, than the latter. This is consistent with an increasing ratio of alkali feldspar to mafic components (including sodic pyroxene and/or amphibole) in the residual solids as melting proceeded.
|
Five glass analyses in xenolith 2 in BL225a are an + c-normative and one is ac-normative. The only significant compositional difference is in the Na2O contents (
2·62 as opposed to 4·94 wt %). Sodium may have been lost from originally peralkaline melts as a volatile (fluoride?) phase. This is consistent with glass data from xenolith 1; the compositional range is very similar to that in xenolith 2 except that Na2O contents are higher (3·64–4·50 wt %) and the glasses peralkaline. The glasses from SMN103, although trachytic, are less evolved than in the other xenoliths (e.g. SiO2 62·44–64·13 wt %; an
5·7). The residual feldspar is also more calcic (An32) and it is likely that this xenolith shows partial melting of a rock more mafic than syenite (syenodiorite?). Given that peralkaline glasses have been formed by partial melting of the xenoliths, it is instructive to compare their compositions with those of the eruptive comendites of the GOVC. The glasses plot, with the exception of two of pantelleritic composition, in the comendite field on an FeO*–Al2O3 plot (Fig. 20), forming an array parallel to, but displaced above, the eruptive comendites and pantellerites. The glass data from the xenoliths show, therefore, that smaller degrees of melting of syenites inferred to contain sodic pyriboles generated peralkaline rhyolitic melts similar, but not identical, to the GOVC comendites, whereas more advanced degrees of melting, and melting of more mafic xenoliths, generated an ± c-normative trachytes and rhyolites that had no extrusive equivalents.
|
| PETROGENESIS |
|---|
Causes of the scatter in Group 1
At a given MgO value, the Group 1 rocks show significant ranges of many elements (Fig. 16a and b) and normative characteristics. There were undoubtedly several liquid lines of descent during magma evolution. The variability was extended by basalt–basalt (or hawaiite) mixing in the mafic rocks and mixing between various combinations of mugearite, benmoreite, trachyte and rhyolite in the more evolved rocks. The amount of mixing is difficult to quantify because it caused variations similar to those from crystal–melt processes, although it is usually possible to identify the end-members via the phenocryst assemblages.
Post-eruptive modifications
Peralkaline rocks are prone to compositional modifications on crystallization and/or secondary hydration, particularly loss of Na, oxidation of Fe, and mobility of Rb and Sr. Some Olkaria trachytes are peralkaline and, without analyses of glassy analogues, we cannot preclude the possibility that those rocks lost some Na on post-eruptive crystallization, a possibility also suggested by the fact that some of the an-normative trachytes are also mildly c-normative (c.
2·5%). Plots of potentially mobile trace elements (e.g. Rb) against immobile elements (e.g. Nb, Zr) for the trachytes are linear, strongly suggesting that the trace elements have not been significantly mobilized during crystallization.
Feldspar accumulation
The trachytes contain up to 31% alkali feldspar phenocrysts (Table 1) and some rocks may be feldspar-acumulitic. However, there is no correlation in the trachytes between the modal abundance of phenocrysts and the whole-rock Al2O3 and FeO* contents, which would be expected to increase and decrease, respectively, as a result of alkali feldspar accumulation. We infer that feldspar accumulation has not been a significant factor in the trachytes and that the high modal proportions reflect large degrees of crystallization in relatively static magmas.
Resorption of feldspar xenocrysts
A large proportion of the Olkaria rocks contain up to
10% modally of variably resorbed alkali feldspar xenocrysts. Here we assess the compositional effects of the resorption.
Judging from rocks that do not contain feldspar xenocrysts and are not mixed, Ba abundances increase steadily from
300 ppm in the basalts to a maximum of
1000 ppm at
2·5 wt % MgO and then drop sharply to <500 ppm in the trachytes (Fig. 16b). However, 13 samples contain in excess of 2000 ppm Ba. Simple mass-balance calculations show that addition of 3·7 wt % of alkali feldspar with 3 wt % Ba (within the range of values of the Olkaria feldspars; Tables 2 and 3) would increase the Ba content of BB85.231 from 2630 ppm to 4000 ppm, close to the maximum in the suite (Fig. 16b). It is very likely that a significant proportion of the Ba budget of rocks with >1000 ppm has been derived from feldspar resorption, the actual abundances being determined by the amount and composition of the feldspar being resorbed. Anomalous Ba enrichment in alkaline rock suites has previously been ascribed to alkali feldspar resorption by Edgar et al. (2002
) and Sumner & Wolff (2003
). High Ba values can also be used as a guide to the likely former presence of feldspar xenocrysts, now completely resorbed (e.g. in samples SMN171 and BB85.231). Alkali feldspar resorption also increases the Al2O3, Na2O, K2O, and usually SiO2, whole-rock abundances, whereas the content of most other elements is diluted.
Davies & Macdonald (1987
) noted that many Ndabibi basalts show positive Eu and Sr anomalies when plotted on primitive mantle-normalized diagrams, which they ascribed to accumulation and/or resorption of plagioclase. Our new data allow us to add some detail. As the An contents of the feldspar phenocrysts and xenocrysts decrease, Eu and Sr abundances initially increase, reaching maxima of
28 ppm at An20 and 1877 ppm at An30, respectively. Abundances of both elements then sharply decrease. As the highest Eu and Sr contents in an Olkaria whole-rock are 6·7 ppm and 550 ppm, respectively (Tables 7 and 8), feldspar resorption has the potential to increase those values. For example, addition of 5% (wt) of feldspar An26·6 (analysis 4HC in BB85.40; ST2) to basalt SMN195 (Table 8) increases the Eu abundance from 1·53 to 2·51 ppm. Importantly, the Eu anomaly (Eu/Eu*) increases from 1·02 to 1·73. It is likely, therefore, that positive Eu (and Sr) anomalies in the Olkaria rocks (e.g. Eu/Eu* = 2·04 in SMN165b; Table 8) are a result of the resorption of ternary feldspars in the compositional range An30–20. Furthermore, feldspar resorption may have counteracted the development of negative Eu anomalies during fractional crystallization.
Significance of new analyses for the Daly Gap
The eruptive rocks of Group 1, taken to include lavas, high-level dykes, and blocks and bombs in pyroclastic deposits, show an absence of compositions in the range 49·4–62·9 SiO2 wt % (MgO 4·75–0·52 wt %). This gap is filled by the magmatic inclusions (50·8–64·9 SiO2 wt %; 4·39–0·03 MgO wt %). The majority of the inclusions are mixed; nevertheless, the compositional range shown by the non-mixed varieties suggests that there was a continuum of melt compositions from basalt to trachyte (Fig. 15, inset). This is consistent with the essentially continuous ranges shown by the cores of feldspar (Fig. 7) and clinopyroxene (Fig. 12) phenocrysts. In addition, the presence of trachytes with q > 10% and rhyolites in the Olkaria and Tandamara Trachyte Formations provides a transition to the comendites (Fig. 14). Thus the silica-gap between mafic and silicic rocks identified in the GOVC by Davies & Macdonald (1987
) and Macdonald et al. (1987
) does not exist among the magmatic products. A complete spectrum of compositions has existed in the wider plumbing system for the past 20 kyr, although basalts have been erupted largely outside the caldera and mugearites and benmoreites tended to reach the surface only as magmatic inclusions.
A similar situation, where the Daly gap has been filled by magmatic inclusions (enclaves), has been described from the Gedemsa volcano, Ethiopia (Peccerillo et al., 2003
), the Alid volcanic centre, Eritrea (Lowenstern et al., 2006
) and Pantelleria, Italy (Ferla & Meli, 2006
). Bacon & Metz (1984
) recorded that intermediate (andesitic) compositions occur as inclusions in high-silica rhyolites of the Coso volcanic field in California, and Bachmann et al. (2002
) recorded hybrid andesitic inclusions in the upper intracaldera Fish Canyon Tuff, Colorado. In both the latter cases, the andesites were formed by mixing of basaltic and silicic magmas, as in the GOVC where the majority of inclusions are mixes of two or more components.
Role of fractional crystallization
Major element variations in the Group 1 rocks (Fig. 16a) are qualitatively consistent with fractionation of olivine–clinopyroxene–plagioclase assemblages in the range 8–3·5 wt % MgO, clinopyroxene–plagioclase–FeTi-oxides–apatite from 3·5 to 1·5 wt % MgO, and sanidine–clinopyroxene–FeTi-oxides–apatite ± olivine in the trachytes. We do not present quantitative models of fractional crystallization, because the multi-lineage nature and the complexity of phenocryst or xenocryst relationships (e.g. Figs 7b, 11 and 13) make it virtually impossible to identify likely parent–daughter pairs. Instead, we estimate degrees of crystallization using ITE enrichment factors. The maximum enrichment factor along the inferred liquid line of descent are
12 for Zr, Rb and Th. Assuming that these elements have been fairly incompatible, the enrichment factors indicate that
8% of melt (of quartz trachytic to rhyolitic composition) remained after fractionation of parental basalts with 8 wt % MgO, as in the most magnesian members of the Ndabibi Basalt Formation. This value is comparable with those estimated for similar magmatic suites, such as the Gedemsa, Ethiopia, basalt–pantellerite sequence (
10%; Peccerillo et al., 2003
) and the basalt–rhyolite suite of the Alid centre, Eritrea (
10–20%; Lowenstern et al., 2006
).
Enrichment factors for Ce and Nb are slightly lower at
10. The fractionation of these elements relative to Zr and Rb occurs in rocks with MgO <3·5 wt % and can probably be ascribed to Ce entry into apatite and chevkinite (ST6) and Nb entry into ilmenite (Table 6). The enrichment factor for Y is
4 and it has clearly been fractionated relative to the other ITE; we discuss this below.
Relationship of comendites to trachytes
Compositional relationships between the trachytes and the least evolved post-caldera comendites are explored on DI-element/ratio (Fig. 17) and FeO*–Al2O3 (Fig. 20) plots. On the FeO*–Al2O3 diagram, the trachytes plot close to an extension of data for the trachytes from Longonot and Eburru. However, the tie-line from whole-rock to interstitial glass (residual melt) in SMN199 trends towards the comendite field and it is significant that PI increases from 0·94 to 1·01. Thus 26·5% (by volume) crystallization of an alkali feldspar–olivine–clinopyroxene–FeTi-oxide assemblage (SMN199; Table 1) generated a peralkaline residuum. On the DI plots (Fig. 17), the comendites generally lie on extensions of the trachyte field, consistent with a crystal fractionation relationship. The sharp drop in Al2O3 and Zr would require fractionation of alkali feldspar and zircon. It is noticeable that the SMN199 tie-line trends towards the comendites for Al2O3. Because Zr enters zircon preferentially over Hf, zircon fractionation results in lower Zr/Hf in residual melts (Lowery Claiborne et al., 2006
). Trachytes from the Olkaria and Tandamara Trachyte Formations have Zr/Hf in the range 46–50 (Table 8), whereas the least evolved comendites range from 24 to 33 (Macdonald et al., 1987
).
The fractional crystallization model for the origin of the comendites resolves a problem raised by the partial melting model (Davies & Macdonald, 1987
; Macdonald et al., 1987
; Black et al., 1997
; Heumann & Davies, 2002
). Low Ba–Sr–Eu rhyolites cannot be formed directly by partial melting. Extensive feldspar fractionation from intermediate magma is required, with the formation of cumulates enriched in Ba, Sr and Eu, as we have documented at Olkaria.
The role of accessory mineral fractionation
We noted above that ITE ratios tend to remain about constant in the basalt–benmoreite range but to show greater ranges in the trachytes. Some process, or processes, fractionated the ratios in the salic magmas. Judging from variations in ITE ratios against DI, the relative degrees of incompatibility were Zr > Ce + Nb > Rb + Th > Y. Rubidium was partitioned into alkali feldspar (Table 3), although the KD was <1. Y may have been fractionated by apatite; Y2O3 levels reach 1·24 wt % in that phase (unpublished data) and there is a positive correlation in the whole-rocks between P2O5 and Y. However, apatite fractionation should result in decreased Ce/Y ratios, whereas they increase with increasing DI and there is a negative correlation between P2O5 and Ce. Fractionation of LREE-enriched chevkinite would have similar results and chevkinite also concentrates Th (ST6). It seems more likely that Ce (and the other LREE) were partitioned into clinopyroxene, where KD values of up to three have been recorded from salic rocks (Nash & Crecraft, 1985
; Mahood & Stimac, 1990
; Sisson, 1991
). Either we have not recognized the Y-fractionating phase in the trachytes, or its distribution was not controlled simply by crystal–melt processes, as we suggest below.
Role for a volatile phase?
Certain aspects of GOVC geochemistry may require a role for a volatile phase or phases in magma evolution. We appreciate that the following discussion is speculative but, taken together, the geochemical data point to processes other than simple melt–solid interactions operating in the Olkaria magmas.
(1) Four magmatic inclusions and a trachyte lava have negative Ce anomalies, with Ce/Ce* in the range 0·91–0·33. Development of negative Ce anomalies is commonly ascribed to oxidation of Ce3+ to Ce4+ and its being taken up by a halogen-rich oxidizing vapour and lost from the magma. Such a situation is quite probable in the GOVC, where halogen-rich peralkaline magmas were commonly mixing with and/or being heated by more mafic melts, potentially promoting volatile loss. Direct evidence of volatile loss is seen in the vesiculation common at the borders between magmatic inclusions and comenditic hosts. Uniquely in our dataset, magmatic inclusion SMN121 shows a positive Ce anomaly (Ce/Ce* = 1·4). This may be a case of Ce addition through interaction with a volatile phase, a feature consistent with the rock's high Pb content (243 ppm).
(2) Many rocks in the Olkaria suite show anomalous element enrichments; for example, in Y, HREE, U, Cu, Pb and Zn. Yttrium has low solubility in CO2-rich fluids but, with U, may have partitioned into carbonate melts separating from the silicate magmas. We noted above that the enrichment factor for Y (
4) was much smaller than for the other ITE (
10). A possible explanation is that Y has been lost from the major parts of the system but concentrated in others. We note that Macdonald et al. (1993
) documented separation of a carbonate phase from trachytes of the neighbouring Suswa volcano, although no data for Y distribution between the two melts were presented. Enrichment of Pb and Zn may have been via a sulphur-rich fluid phase; sulphide blebs are common inclusions in the FeTi-oxide phenocrysts, indicating that the melts were S-bearing. Lowenstern et al. (1991
) found that Cu is several hundred times more concentrated in magmatic vapour than in the coexisting melt represented by a pantellerite from Pantelleria and suggested that Cu resides in a vapour phase in some shallow magma chambers, perhaps as CuCl.
(3) Relative to rocks of similar MgO content in Group 1, the ne-normative rocks of Group 3 generally have high abundances of Na2O, Rb, Th and the REE, especially perhaps the HREE (magmatic inclusions SMN153, SMN169ii and SMN121 show only modest LREE enrichment ([La/Yb]CN
9)). The rocks provide evidence of extensive interaction with the host comendites, including the common growth of fluorite. We suggest that there was selective diffusion into the inclusions of a Na–F-rich fluid, which imparted, inter alia, a ne-normative, and in some cases peralkaline, character to them.
| DISCUSSION |
|---|
Nature of the Olkaria plumbing system
We have suggested that the following mechanisms operated during the evolution of the Olkaria rocks: (1) crystal fractionation of basaltic magmas generated mugearitic magmas; (2) mixing between basalt magmas was common; (3) further crystal fractionation along several liquid lines of descent produced benmoreitic and trachytic magmas of Groups 1 and 3 and the post-caldera comendites; (4) two- and three-component mixing between mugearites, benmoreites and, less commonly, trachytes produced a range of mingled and hybrid magmas; (5) mixing between basalts and trachytes/rhyolites generated the rocks of Group 2; (6) magma types less evolved than trachytes commonly contain significant amounts of feldspar xenocrysts, presumably derived from the margins of the conduits and reservoirs through which the ascending magmas passed; (7) fluid phases of various composition escaped from, and interacted with, the magmas, possibly a result of heating during mixing; (8) eruption of the comendites was sometimes accompanied by mixing with small amounts of magmas, themselves mixed during an earlier phase.
These features have, to varying degrees, distinguished the GOVC plumbing system since eruption of the oldest Ndabibi basalts (0·5 Ma?) and have affected magmas erupted both inside and peripheral to the ring structure. The system is clearly very open. It is also a very dynamic environment, where crystals with a variety of growth and resorption histories have been juxtaposed shortly before eruption (see Bachmann et al., 2002
). The repeated occurrence of different magma types during very short periods (e.g. the <3 kyr of O4 activity) strongly suggests that the plumbing system comprises several independent reservoirs and conduits (compare Reubi & Nicholls, 2004
, 2005
). The reservoirs must be small enough to have been repeatedly pervaded by recharge batches of mafic magma (compare Mammoth Mountain, California; Hildreth, 2004
). The concept of the lower parts of the plumbing system being a focus for basaltic dykes is consistent with geophysical evidence. A series of short-wavelength gravity highs (–190 mGal) superimposed on the broad negative anomaly in the axial region of the rift and coincident with the central complexes of the Kenya Dome has been modelled by Swain (1992
) as resulting from dyke injection in the central 40 km of the rift. Perhaps 22–26% of the present crust, down to a depth of 22 km, consists of this intruded material.
Magma mixing has been ubiquitous in the development of the GOVC. Indeed, it is uncommon to find rocks in the mugearite–benmoreite range that are not mixed. The swirling textures common between comendites and magmatic inclusions and the intimate admixture of two or three components in some inclusions suggest that significant mixing took place in the conduits. However, the longer times required for hybridization may point to some mixing in small reservoirs. In the basalt–trachyte/rhyolite mixes of Group 2, mingling was more common than hybridization, perhaps pointing either to insufficient time for more complete mixing and/or a viscosity or density barrier preventing mixing. In Group 1 rocks, mugearites and benmoreites have been able to mix with trachytes only where the trachytes have been phenocryst-poor; the crystal-rich eruptive trachytes show relatively little evidence of mingling, although some rocks are hybrids (Appendix). Magmatic inclusions in the comendites show only local interactions with their hosts; more thorough mixing has presumably been prevented by differences in viscosity.
Rejuvenation of earlier mushes and plutonic bodies
Many GOVC rocks in the range basalt–benmoreite, including some erupted outside the ring structure, carry significant amounts of alkali feldspar ± plagioclase xenocrysts. Some also carry gabbro–syenite xenoliths, which have sometimes been partially melted to form silicic melts. The compositions of the xenocrysts are very similar to those of the phenocrysts in the mugearite–trachyte range (Fig. 7a and b) and it seems likely that the bulk of the xenocrysts represent recycled crystals from intrusive rocks or crystal mushes (i.e. are antecrysts). The presence of holocrystalline comagmatic xenoliths indicates that some parts of the plumbing system had crystallized completely. The xenocrysts occur both within and outside the ring structure, possibly indicating that the plutons and mushes were derived from a caldera-scale mass, possibly the residue of the system that generated the syn-caldera rocks. No post-caldera pantellerites have been erupted within the ring fracture and it seems possible that the pre-caldera magma reservoir was totally evacuated and that caldera formation resulted in collapse of the roof onto the floor, as envisaged by Macdonald et al. (1995
) for the basalt–trachyte caldera volcano Silali, in the north–central rift. Later dykes were then able to cut through the former chamber, incorporating xenocrystic and xenolithic material. Material erupted peripherally to the ring structure may have been carried there laterally; there is some evidence in the young volcanoes of the central rift valley of magmas moving between complexes along dykes (Macdonald et al., 1995
; Scott & Skilling, 1999
; Rogers et al., 2004
).
Thermal rejuvenation of near-solidus high-level intrusions shortly before eruption is common in systems of very different volumes (Bachmann et al., 2002
) and is probably induced by recurrent pulses of basaltic intrusion and crystallization. Mafic to intermediate magmas may also evolve H2O–F–S–Cl-rich fluids that flux melting in the crystal mush (Bachmann et al., 2002
; the gas sparging effect of Bachmann & Bergantz, 2006
) and there is some evidence in the GOVC for the release of such fluids.
Figure 21 is a schematic illustration of the development of the GOVC, based on the various constraints described above. Given the complexity of magma relationships, the scheme is necessarily much simplified and is certainly not definitive. However, the GOVC seems to have similarities to the plumbing system of Batur volcano, Bali, prior to the second caldera-forming eruption there at
20 ka, as drawn by Reubi & Nicholls (2005
).
|
The system is driven by basaltic magmas intruded as dykes at mid- to upper-crustal levels. Mingling between basalt magmas sometimes occurs. Basalt ponds at various levels, to form the crystal mushes and/or gabbros recorded as mafic xenoliths in more evolved mixed rocks. Crystal-poor melts of mugearitic–benmoreitic composition are expelled from the reservoirs and rise through dykes, mixing frequently. Incorporation of partially to completely solidified plutonic material is almost ubiquitous during magma ascent. The intermediate magmas pond and fractionate through to trachytes, with crystal residua of syenodioritic composition.
The expelled trachytes pond at
5–6 km depth, which is comparable with the depth estimated from seismic data to be the interface between the Pan-African basement and the Miocene–Holocene volcaniclastic infill beneath the GOVC (Mooney & Christensen, 1994
; Mechie et al., 1997
). This may also have been the level at which the pantelleritic magmas erupted during caldera collapse were stored. There is extensive crystallization within the trachytic reservoirs, possibly promoted by gas exsolution, and magma mixing is relatively restricted. The compositional range in the trachytes points to the need for several of these mush-rich reservoirs, which equate, on a smaller scale, with the viscosity barrier of Smith (1979
) and mush zone envisaged by Bachmann & Bergantz (2004
), Hildreth (2004
) and Hildreth & Wilson (2007
) for large silicic systems. At a critical value of crystallization, interstitial rhyolitic melts escape upwards as dykes and pods.
Eruption of the rhyolites is either promoted by, or entrains, magmas from deeper levels, ranging from mugearite to trachyte and normally forming two- to four-component mixes. Mafic magmas extruded outside the ring fracture have also incorporated plutonic material and have mixed with rhyolitic magmas.
Origin of the peralkaline rhyolites
Comenditic melts, now represented by glass, have formed by two mechanisms in the GOVC: fractional crystallization of trachyte and partial melting of syenites. The trachytes themselves formed by fractionation of basalts and the inferred liquid line of descent to the least evolved comendites is long, requiring >90% of crystallization. Given the generally phenocryst-poor nature of the comendites, extraction from the parental trachytic mushes must have been very efficient. Extraction would have been aided by the low viscosity of the volatile-rich, peralkaline melts and by gas-driven filter pressing related to volatile saturation and exsolution as crystallization proceded (Sisson & Bacon, 1999
; Bachmann & Bergantz, 2004
; Simon & Reid, 2005
).
Generation of the comendites by partial melting of syenites cannot be excluded on geochemical grounds and the process itself may have been promoted by fluxing of volatiles released from underlying magmas. However, the evidence from the rocks is that partial melting acted on a very small (millimetre) scale and produced a range of rhyolitic compositions, including an- and c-normative types that have not been erupted as discrete magma bodies. It may well be, however, that at least some of the comendites are mixtures of melts formed by fractional crystallization and partial melting.
| CONCLUSIONS |
|---|
(1) There is a continuum, in the GOVC, of compositions from basalt to comendite, the SiO2 range 50–63 wt % being represented only by magmatic inclusions.
(2) The trachytes were formed by fractional crystallization of basaltic magmas via mugearites and benmoreites.
(3) Magma mixing has been ubiquitous in all rock types, although hybridization has been less common in the trachytes and rhyolites.
(4) Xenocrystic alkali feldspar occurs in 51% of specimens, up to 9% modally; plagioclase xenocrysts are less common (21%) and less abundant (usually in trace amounts). Some of the xenocrysts have probably been derived through remobilization of partially or completely solidified plutonic material. Xenocryst resorption has had profound effects on certain geochemical features, such as Ba abundances and Eu/Eu*.
(5) The occurrence of peralkaline rhyolitic glasses indicates that the eruptive comendites could have formed by two mechanisms: (a) partial melting of syenites; (b) fractional crystallization of metaluminous trachyte.
(6) Separation of Na- or K-enriched fluids may have generated, inter alia, anomalous trace element enrichments in some rocks and a group of silica-undersaturated peralkaline rocks.
(7) The Olkaria system comprises a plexus of dykes, pods and partially to completely solidified intrusions, where ductile deformation promotes extraction, aggregation and mixing of a wide range of melts.
| SUPPLEMENTARY DATA |
|---|
Supplementary data for this paper are available at Journal of Petrology online.
| APPENDIX: SAMPLE DETAILS |
|---|
Appendix:
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Coordinates for locations refer to UTM grid numbers on East Africa 1:50 000 (Y503) topographic maps. Units refer to the stratigraphic subdivisions of Clarke et al. (1990
). Gp:TAS, Groups 1–3 as defined petrographically and geochemically in the text and classification according to the TAS system (T, trachyte; R, rhyolite; B, basalt; M, mugearite; BA, basaltic andesite; Bn, benmoreite; A, andesite; Pt, phonotephrite; Tp, tephriphonolite; P, phonolite). Mix, the number of components in mingled rocks; hyb, rock believed to be hybridized. —, not present; =, no information; ol, olivine; cpx, clinopyroxene; pl, plagioclase; ox, FeTi-oxides; af, alkali feldspar; ap, apatite; z, zircon; ch, chevkinite; fl, fluorite.
*In SMN105, more than one mafic component.
**BL469 may contain more than 2 mix components.
| ACKNOWLEDGEMENTS |
|---|
The Natural Environment Research Council (UK) is thanked for provision of a research studentship to A.S.M. and various research grants to R.M. R.M. gratefully acknowledges tenure of a Visiting Research Fellowship at the Open University during preparation of this paper. Charlie Bacon and John Wolff provided very detailed and helpful journal reviews and Alastair Lumsden for his careful and patient editorial work. We thank Nic Odling (Edinburgh) for help with XRF analyses, and Michelle Higgins and Kay Green (Open University) for the preparation of thin sections. Dr Luc Lepage kindly supplied a copy of the worksheets for his ILMAT program. Field work in Kenya was greatly assisted by Geoffrey Muchemi and Johnson Mungania at the Olkaria Geothermal Project. The Warden of Hell's Gate National Park and John Mackay and Tim Trent of Oserian Development Co. Ltd made access to parts of the field area possible.
*Corresponding author. E-mail: r.macdonald{at}lancaster.ac.uk
| REFERENCES |
|---|
Andersen DJ, Lindsley DH. New (and final!) models for the Ti-magnetite/ilmenite geothermometer and oxygen barometer. EOS Transactions, American Geophysical Union (1985) 66:416.
Bachmann O, Bergantz GW. Rejuvenation of the Fish canyon magma body: a window into the evolution of large-volume silicic magma systems. Geology (2003) 31:789–792.
Bachmann O, Bergantz GW. On the origin of crystal-poor rhyolites: extracted from batholitic crystal mushes. Journal of Petrology (2004) 45:1565–1582.
Bachmann O, Bergantz GW. Gas percolation in upper-crustal silicic crystal mushes as a mechanism for upward heat advection and rejuvenation of near-solidus magma bodies. Journal of Volcanology and Geothermal Research (2006) 149:85–102.[CrossRef][Web of Science]
Bachmann O, Dungan MA, Lipman PW. The Fish Canyon magma body, San Juan volcanic field, Colorado: rejuvenation and eruption of an upper-crustal batholith. Journal of Petrology (2002) 43:1469–1503.
Bacon CR, Hirschmann MC. Mg/Mn partitioning as a test for equilibrium between coexisting Fe–Ti oxides. American Mineralogist (1988) 73:57–61.[Abstract]
Bacon CR, Lowenstern JB. Late Pleistocene granodiorite source for recycled zircon and phenocrysts in rhyodacite lava at Crater Lake, Oregon. Earth and Planetary Science Letters (2005) 233:277–293.[CrossRef][Web of Science]
Bacon CR, Metz J. Magmatic inclusions in rhyolites, contaminated basalts, and compositional zonation beneath the Coso volcanic field, California. Contributions to Mineralogy and Petrology (1984) 85:346–365.[CrossRef][Web of Science]
Black S, Macdonald R, Kelly MR. Crustal origin for peralkaline rhyolites from Kenya: evidence from U-series disequilibria and Th-isotopes. Journal of Petrology (1997) 39:995–1008.[CrossRef][Web of Science]
Bliss CM. Geology and petrochemistry of the Naivasha area, Kenya Rift Valley. In: PhD thesis (1979) Lancaster University.
Bone BD. The geological evolution of the S. W. Naivasha volcanic complex, Kenya. In: PhD thesis (1987) Lancaster University.
Burt RM, Brown SJA, Cole JW, Shelley D, Waight TE. Glass-bearing plutonic fragments from ignimbrites of the Okataina caldera complex, Taupo Volcanic Zone, New Zealand: remnants of a partially molten intrusion associated with preceding eruptions. Journal of Volcanology and Geothermal Research (1998) 84:209–237.[CrossRef][Web of Science]
Charlier BLA, Wilson CJN, Lowenstern JB, Blake S, van Calsteren PW, Davidson JP. Magma generation at a large, hyperactive silicic volcano (Taupo, New Zealand) revealed by U–Th and U–Pb systematics in zircons. Journal of Petrology (2005) 46:3–32.
Clarke MC, Woodhall DG, Allen D, Darling G. Geological, volcanological and hydrogeological controls on the occurrence of geothermal activity in the area surrounding Lake Naivasha, Kenya (1990) Nairobi: Ministry of Energy. 138.
Davies GR, Macdonald R. Crustal influences in the petrogenesis of the Naivasha basalt–comendite complex: combined trace element and Sr–Nd–Pb isotope constraints. Journal of Petrology (1987) 28:1009–1031.
Edgar CJ, Wolff JA, Nichols HJ, Cas RAF, Marti J. A complex Quaternary ignimbrite-forming phonolitic eruption: the Poris Member of the Diego Hernández Formation (Tenerife, Canary Islands). Journal of Volcanology and Geothermal Research (2002) 118:99–130.[CrossRef][Web of Science]
Feeley TC, Dungan MA. Compositional and dynamic controls on mafic–silicic interactions at continental arc volcanoes: evidence from Cordon El Guadal, Tatara–San Pablo Complex, Chile. Journal of Petrology (1996) 37:1547–1577.
Ferla P, Meli C. Evidence of magma mixing in the Daly Gap of alkaline suites: a case study from the enclaves of Pantelleria (Italy). Journal of Petrology (2006) 47:1467–1507.
Fitton JG, Saunders AD, Larsen LM, Hardarson BS, Norry MJ. Volcanic rocks from the southeast Greenland margin at 63°N: composition, petrogenesis and mantle sources. In: Proceedings of the Ocean Drilling Program Scientific Results, 152—Saunders AD, Larsen HC, Wise SH, eds. (1998) College Station, TX: Ocean Drilling Program. 331–350.
Heumann A, Davies GR. U–Th disequilibrium and Rb–Sr age constraints on the magmatic evolution of peralkaline rhyolites from Kenya. Journal of Petrology (2002) 43:557–577.
Hibbard MJ. The magma mixing origin of mantled feldspars. Contributions to Mineralogy and Petrology (1981) 76:158–170.[CrossRef][Web of Science]
Hildreth W. Volcanological perspectives on Long Valley, Mammoth Mountain, and Mono Craters: several contiguous but discrete systems. Journal of Volcanology and Geothermal Research (2004) 136:169–198.[CrossRef][Web of Science]
Hildreth W, Wilson CJN. Compositional zoning of the Bishop Tuff. Journal of Petrology (2007) 48:951–999.
Le Bas MJ, Le Maitre RW, Streckeisen A, Zanettin B. A chemical classification of volcanic rocks based on the total alkali–silica diagram. Journal of Petrology (1986) 27:745–750.
Lepage LD. ILMAT: an Excel program for ilmenite–magnetite geothermometry and geobarometry. Computers and Geosciences (2003) 29:673–678.[CrossRef]
Lofgren GE. An experimental study of plagioclase crystal morphology: isothermal crystallization. American Journal of Science (1974) 274:243–273.[Abstract]
Lofgren GE. Experimental studies on the dynamic crystallization of silicate melts. In: Physics of Magmatic Processes—Hargraves RB, ed. (1980) Princeton, NJ: Princeton University Press. 487–511.
Lowenstern JB, Mahood GA, Rivers ML, Sutton SR. Evidence for extreme partitioning of copper into a magmatic vapour phase. Science (1991) 252:1405–1409.
Lowenstern JB, Charlier BLA, Clynne MA, Wooden JL. Extreme U–Th disequilibrium in rift-related basalts, rhyolites and granophyric granite and the timescale of rhyolite generation, intrusion and crystallization at Alid volcanic center, Eritrea. Journal of Petrology (2006) 47:2105–2122.
Lowery Claiborne L, Miller CF, Walker BA, Wooden JL, Mazdab FK, Bea F. Tracking magmatic processes through Zr/Hf ratios in rocks and Hf and Ti zoning in zircons: An example from the Spirit Mountain batholith, Nevada. Mineralogical Magazine (2006) 70:517–543.
Macdonald R. Nomenclature and petrochemistry of the peralkaline oversaturated extrusive rocks. Bulletin Volcanologique (1974) 38:398–516.
Macdonald R, Scaillet B. The central Kenya peralkaline province: Insights into the evolution of peralkaline salic magmas. Lithos (2006) 91:59–73.[CrossRef][Web of Science]
Macdonald R, Davies GR, Bliss CM, Leat PT, Bailey DK, Smith RL. Geochemistry of high-silica rhyolites, Naivasha, Kenya rift valley. Journal of Petrology (1987) 28:979–1008.
Macdonald R, Kjarsgaard BA, Skilling IP, Davies GR, Hamilton DL, Black S. Liquid immiscibility between trachyte and carbonate in ash flow tuffs from Kenya. Contributions to Mineralogy and Petrology (1993) 114:276–287.[CrossRef][Web of Science]
Macdonald R, Davies GR, Upton BGJ, Dunkley PN, Smith M, Leat PT. Petrogenesis of Silali volcano, Gregory Rift, Kenya. Journal of the Geological Society, London (1995) 152:703–720.
Macdonald R, Rogers NW, Fitton JG, Black S, Smith M. Plume–lithosphere interactions in the generation of basalts of the Kenya Rift, East Africa. Journal of Petrology (2001) 42:877–900.
Mahood GA, Stimac JA. Trace-element partitioning in pantellerites and trachytes. Geochimica et Cosmochimica Acta (1990) 54:2257–2276.[CrossRef][Web of Science]
Mechie J, Keller GR, Prodehl C, Khan MA, Gaciri SJ. A model for the structure, composition and evolution of the Kenya rift. Tectonophysics (1997) 278:95–118.[CrossRef][Web of Science]
Middlemost EAK. Iron oxidation ratios, norms and the classification of volcanic rocks. Chemical Geology (1989) 77:19–26.[CrossRef][Web of Science]
Mooney WD, Christensen NI. Composition of the crust beneath the Kenya rift. Tectonophysics (1994) 236:391–408.[CrossRef][Web of Science]
Nairn IA, Shane PR, Cole JW, Leonard GJ, Self S, Pearson N. Rhyolite magma processes of the
AD 1315 Kaharoa eruption episode, Tarawera volcano, New Zealand. Journal of Volcanology and Geothermal Research (2004) 131:265–294.[CrossRef][Web of Science]
Nakada S, Bacon CR, Gartner AE. Origin of phenocrysts and compositional diversity in pre-Mazama rhyodacite lavas, Crater Lake, Oregon. Journal of Petrology (1994) 35:127–162.
Nash WP, Crecraft HR. Partition coefficients for trace elements in silicic magmas. Geochimica et Cosmochimica Acta (1985) 49:2309–2322.[CrossRef][Web of Science]
Peccerillo A, Barberio MR, Yirgu G, Ayalew D, Barbieri M, Wu TW. Relationships between mafic and peralkaline silicic magmatism in continental rift settings: a petrological, geochemical and isotopic study of the Gedemsa volcano, Central Ethiopian Rift. Journal of Petrology (2003) 44:2003–2032.
Pichavant M, Costa F, Burgisser A, Scaillet B, Martel C, Poussineau S. Equilibration scales in silicic to intermediate magmas—implications for experimental studies. Journal of Petrology (2007) 48:1955–1972.
Ren M, Omenda PA, Anthony EY, White JC, Macdonald R, Bailey DK. Application of the QUILF thermobarometer to the peralkaline trachytes and pantellerites of the Eburru volcanic complex, East African Rift, Kenya. Lithos (2006) 91:109–124.[CrossRef][Web of Science]
Reubi O, Nicholls IA. Magmatic evolution of Batur volcanic field, Bali, Indonesia: petrological evidence for polybaric fractional crystallization and implications for caldera-forming eruptions. Journal of Volcanology and Geothermal Research (2004) 138:345–369.[CrossRef][Web of Science]
Reubi O, Nicholls IA. Structure and dynamics of a silicic magmatic system associated with caldera-forming eruptions at Batur volcanic field, Bali, Indonesia. Journal of Petrology (2005) 46:1367–1391.
Rogers NW, Evans PJ, Blake S, Scott SC, Hawkesworth CJ. Rates and timescales of fractional crystallization from 238U–230Th–226Ra disequilibria in trachyte lavas from Longonot volcano, Kenya. Journal of Petrology (2004) 45:1747–1776.
Rogers NW, Thomas LE, Macdonald R, Hawkesworth CJ, Mokadem F. 238U–230Th disequilibrium in recent basalts and dynamic melting beneath the Kenya rift. Chemical Geology (2006) 234:148–168.[CrossRef][Web of Science]
Scott SC, Skilling IP. The role of tephrochronology in recognising synchronous caldera-forming events at the Quaternary volcanoes Longonot and Suswa, south Kenya Rift. In: Volcanoes in the Quaternary. Geological Society, London, Special Publications—Firth CR, McGuire WJ, eds. (1999) 161:47–67.[CrossRef]
Simon JI, Reid MR. The pace of rhyolite differentiation and storage in an archetypical silicic magma system, Long Valley, California. Earth and Planetary Science Letters (2005) 235:123–140.[CrossRef][Web of Science]
Sisson TW. Pyroxene–high silica rhyolite trace element coefficients measured by ion microprobe. Geochimica et Cosmochimica Acta (1991) 55:1575–1585.[CrossRef][Web of Science]
Sisson TW, Bacon CR. Gas-driven filter pressing in magmas. Geology (1999) 27:613–616.
Smith RL. Ash-flow magmatism. In: Ash-flow tuffes. Geological Society of America, Special Papers—Chapin CE, Elston WE, eds. (1979) 180:5–27.
Smith VC, Shane P, Nairn IA. Reactivation of a rhyolitic magma body by new rhyolitic intrusion before the 15·8 ka Rotorua eruptive episode: implications for magma storage in the Okataina Volcanic Centre, New Zealand. Journal of the Geological Society, London (2004) 161:757–772.
Smith VC, Shane P, Nairn IA. Trends in rhyolite geochemistry, mineralogy, and magma storage during the last 50 kyr at Okataina and Taupo volcanic centres, Taupo Volcanic Zone, New Zealand. Journal of Volcanology and Geothermal Research (2005) 148:372–406.[CrossRef][Web of Science]
Sumner JM, Wolff J. Petrogenesis of mixed-magma, high-grade, peralkaline ignimbrite TL (Gran Canaria): diverse styles of mixing in a replenished, zoned magma chamber. Journal of Volcanology and Geothermal Research (2003) 126:109–126.[CrossRef][Web of Science]
Sun S-S, McDonough WF. Chemical and isotopic systematics of oceanic basalts: implications for mantle composition and processes. In: Magmatism in the Ocean Basins. Geological Society, London, Special Publications—Saunders AD, Norry MJ, eds. (1989) 42:313–345.
Swain CJ. The Kenya rift axial gravity high: a reinterpretation. Tectonophysics (1992) 204:59–70.[CrossRef][Web of Science]
Tamura Y, Yuhara M, Ishii T, Irino N, Shukuno H. Andesites and dacites from Daisen volcano, Japan; partial-to-total remelting of an andesite magma body. Journal of Petrology (2003) 44:2243–2260.
Wolff JA, Gardner JN. Is the Valles Caldera entering a new cycle of activity? Geology (1995) 23:411–414.
![]()
CiteULike
Connotea
Del.icio.us What's this?
This article has been cited by other articles:
![]() |
J. Biggs, E.Y. Anthony, and C.J. Ebinger Multiple inflation and deflation events at Kenyan volcanoes, East African Rift Geology, November 1, 2009; 37(11): 979 - 982. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Macdonald and B. Baginski The central Kenya peralkaline province: a unique assemblage of magmatic systems Mineralogical Magazine, June 15, 2009; 73(1): 1 - 16. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. Macdonald, B. Baginski, H.E. Belkin, P. Dzierzanowski, and L. Jezak REE partitioning between apatite and melt in a peralkaline volcanic suite, Kenya Rift Valley Mineralogical Magazine, April 17, 2009; 72(6): 1147 - 1161. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. S. Marshall, R. Macdonald, N. W. Rogers, J. G. Fitton, A. G. Tindle, K. Nejbert, and R. W. Hinton Fractionation of Peralkaline Silicic Magmas: the Greater Olkaria Volcanic Complex, Kenya Rift Valley J. Petrology, February 1, 2009; 50(2): 323 - 359. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
























