Quartz-in-garnet barometry constraints on formation pressures of eclogites from the Franciscan Complex, California

Determining pressure and temperature variations between high-pressure/low-temperature (HP–LT) eclogite blocks is crucial for constraining end-member exhumation models; however, it has historically been challenging to constrain eclogite pressures due to the high variance associated with this bulk-rock composition. In this work, we utilize quartz-in-garnet elastic barometry to constrain formation pressures of eclogites from the northern (Junction School, Ring Mountain, Jenner Beach) and southern Franciscan Complex (Santa Catalina Island). Multiple eclogite blocks from Jenner Beach are analyzed, and a single eclogite from the other localities. By comparing garnet growth conditions from within a single outcrop and between distinct outcrops, we evaluate the local and regional spatial distribution of P conditions recorded by eclogites. We compare the mean, median, and max pressures between different garnet zones and eclogites. Pressures sometimes exhibit systematic changes across garnet zones; however, some eclogites exhibit no systematic pressure variations across garnet zones. Pressures from northern Franciscan eclogites range from $$\sim $$ ∼ 1.4–1.8 GPa, at an estimated temperature of 500 $$^{\circ }$$ ∘ C; pressures from the Catalina eclogite range from $$\sim $$ ∼ 1.2–1.5 GPa, at an estimated temperature of 650 $$^{\circ }$$ ∘ C. Mean and maximum pressures of different eclogites from the northern Franciscan exhibit negligible differences (< 0.1 GPa). The results are inconsistent with models that propose exhumation of metamorphic blocks from different structural levels, and suggest that now exposed HP–LT eclogites from the northern Franciscan Complex may represent rocks that were coherently underplated, and exhumed from similar structural levels.

Ore Mineralogy, Geochemistry, and Dating (Re-Os, Ar-Ar) of the El Volcán Porphyry/Epithermal System, Maricunga Gold Belt, Northern Chile

The El Volcán gold project (8.9 Moz Au @ 0.71 g/t Au) is located in the Maricunga gold belt in northern Chile, on the flank of the large Cenozoic Copiapó Volcanic Complex. Precious metal mineralization is hosted in two zones (Dorado and Ojo de Agua) of (pervasively) altered Miocene porphyry intrusions and lava flows of andesitic to rhyolitic composition, and in breccias. The ore zones reflect an evolving magmatic-hydrothermal system with mineral assemblages of magnetite-ilmenite-pyrite-molybdenite (early), bornite-chalcopyrite-pyrite-rutile (stage I), chalcocite-chalcopyrite-enargite-fahlore-pyrite (stage II), and chalcopyrite-covellite-pyrite (stage III). Alteration is dominantly of Maricunga-style (illite-smectite-chlorite ± kaolinite), partly obscured by quartz-kaolinite-alunite ± illite ± smectite alteration. Powdery quartz-alunite-kaolinite alteration with native sulfur and cinnabar forms shallow steam-heated zones. Early K-feldspar ± biotite alteration is preserved only in small porphyry cores and in deep drill holes. Most gold is submicrometer size and is in banded quartz veinlets, which are characteristic of the Maricunga gold belt. However, some gold is disseminated in zones of pervasive quartz-kaolinite-alunite alteration, with and without banded quartz veinlets. Minor visible gold is related to disseminated chalcocite-chalcopyrite-enargite-fahlore-pyrite. The lithogeochemical database identifies a pronounced Au-Te-Re signature (&gt;100× bulk crust) of the hydrothermal system. Molybdenum-rich bulk rock (100–400 ppm Mo) has an Re-Os age of 10.94 ± 0.17 Ma (2σ). 40Ar-39Ar ages on deep K-feldspar alteration and on alunite altered rock have the same age within error and yield a combined age of 11.20 ± 0.25 Ma (2σ). The formation of the El Volcán gold deposit took place during the establishment of the Chilean flat-slab setting in a time of increasing crustal thickness when hydrous magmas were formed in a mature arc setting. The vigorous nature of the hydrothermal system is expressed by abundant one-phase vapor fluid inclusions recording magmatic vapor streaming through a large rock column with a vertical extent of ≥1,500 m.

Multi-stage serpentinization of ultramafic rocks in the Manlay Ophiolite, southern Mongolia

Serpentinization of ultramafic rocks in ophiolites is key to understanding the global cycle of elements and changes in the physical properties of lithospheric mantle. Mongolia, a central part of the Central Asian Orogenic Belt (CAOB), contains numerous ophiolite complexes, but the metamorphism of ultramafic rocks in these ophiolites has been little studied. Here we present the results of our study of the serpentinization of an ultramafic body in the Manlay Ophiolite, southern Mongolia. The ultramafic rocks were completely serpentinized, and no relics of olivine or orthopyroxene were found. The composition of Cr-spinels [Mg# = Mg/(Mg + Fe2+) = 0.54 and Cr# = Cr/(Cr + Al) = 0.56] and the bulk rock chemistry (Mg/Si = 1.21–1.24 and Al/Si < 0.018) of the serpentinites indicate their origin from a fore-arc setting. Lizardite occurs in the cores and rims of mesh texture (Mg# = 0.97) and chrysotile is found in various occurrences, including in bastite (Mg# = 0.95), mesh cores (Mg# = 0.92), mesh rims (Mg# = 0.96), and later-stage large veins (Mg# = 0.94). The presence of lizardite and chrysotile and the absence of antigorite suggests low-temperature serpentinization (<300 °C). The lack of brucite in the serpentinites implies infiltration of the ultramafic rocks of the Manlay Ophiolite by Si-rich fluids. Based on microtextures and mineral chemistry, the serpentinization of the ultramafic rocks in the Manlay Ophiolite took place in three stages: (1) replacement of olivine by lizardite, (2) chrysotile formation (bastite) after orthopyroxene and as a replacement of relics of olivine, and (3) the development of veins of chrysotile that cut across all previous textures. The complex texture of the serpentinites in the Manlay Ophiolite indicates multiple stages of fluid infiltration into the ultramafic parts of these ophiolites in southern Mongolia and the CAOB.

Garnet Chemical Zoning Based Thermobarometry: Method Evaluation and Applications in the Menderes Massif, Western Turkey

The garnet chemical zoning method (GZM) is a reliable thermodynamic approach for forward modeling pressure-temperature (P-T) paths using observed garnet and bulk rock compositions. However, intracrystalline diffusion is known to compromise the integrity of GZM modeled garnet-growth P-T paths. For this reason, extracting reliable metamorphic estimates from garnet-bearing schists in the Central Menderes Massif (CMM), western Turkey, has been difficult. To evaluate the impact of diffusion on GZM, we simulate garnet growth and diffusion for an average metapelite using the program Theria_G. Modeled garnet compositions from four simulations are used to estimate P-T conditions and paths by GZM, which are compared against Theria_G specified P-T-t trajectories. Factors influencing results are heating/cooling rate, grain size, and peak T. At a maximum T of 610 °C, both undiffused and diffused garnet compositions returned estimates comparable to prescribed conditions regardless of heating/cooling rate. Diffused profiles from simulations reaching a maximum T of 670 °C also reproduced prescribed P-T paths if tectonism occurred at high heating/cooling rates (50 °C/my). From these insights and additional Theria_G simulation-derived observations for CMM garnets, we deduce that metamorphism in the region exceeded 650 °C and achieved a maximum burial P between 8–10 kbar prior to Cenozoic exhumation.

A simple and robust method for calculating temperatures of granitoid magmas

Calculating the temperatures of magmas from which granitoid rocks solidify is a key task of studying their petrogenesis, but few geothermometers are satisfactory. Zircon saturation thermometry has been the most widely used because it is conceptually simple and practically convenient, and because it is based on experimental calibrations with significant correlation of the calculated zircon saturation temperature (TZr) with zirconium (Zr) content in the granitic melt (i.e., TZr ∝ ZrMELT). However, application of this thermometry to natural rocks can be misleading, resulting in the calculated TZr having no geological significance. This thermometry requires Zr content and a compound bulk compositional parameter M of the melt as input variables. As the Zr and M information of the melt is not available, petrologists simply use bulk-rock Zr content (ZrBULK-ROCK) and M to calculate TZr. In the experimental calibration, TZr shows no correlation with M, thus the calculated TZr is only a function of ZrMELT. Because granitoid rocks represent cumulates or mixtures of melt with crystals before magma solidification and because significant amount Zr in the bulk-rock sample reside in zircon crystals of varying origin (liquidus, captured or inherited crystals) with unknown modal abundance, ZrBULK-ROCK cannot be equated with ZrMELT that is unknown. Hence, the calculated magma temperatures TZr using ZrBULK-ROCK have no significance in both theory and practice. As an alternative, we propose to use the empirical equation $$T_{SiO_{2}}$$ T S i O 2 (°C) = -14.16 × SiO2 + 1723 for granitoid studies, not to rely on exact values for individual samples but focus on the similarities and differences between samples and sample suites for comparison. This simple and robust thermometry is based on experimentally determined phase equilibria with T ∝ 1/SiO2.

Geochemical and sedimentary constraints on the formation of the Venta Micena early Pleistocene site (Guadix-Baza Basin, Spain)

Despite the paleontological relevance of the terrestrial Early Pleistocene Venta Micena bonebed (Baza Basin, Spain), it lacks a comprehensive geochemical/sedimentological study. Here, we demonstrate that the 1.5-m-thick Venta Micena limestone formed in a relatively small freshwater wetland/pond located at the periphery of the large saline Baza paleolake. Two microfacies are observed, with high and low contents of invertebrate fossils, and which originated in the centre and margin of the wetland, respectively. X-ray diffraction (XRD) mineralogy and paleohydrological characterization based on ostracod and bulk-rock geochemistry (δ13C and δ18O) indicate that the limestone reflects a general lowstand of the Baza lake, permitting the differentiation of freshwater wetlands that were fed by adjacent sources. Conversely, during highstands, the Baza lake flooded the Venta Micena area and the freshwater fauna was replaced by a saline one. Bulk-rock isotopic data indicate that the lower interval C1 of the limestone (bone-rich in marginal settings) displays general negative values, while the upper interval C2 (bone free) displays less negative values. The bones of predated mammals accumulated in the marginal areas, which were flooded and buried by recurring water-table fluctuations. Lake dynamics played a critical role in bone accumulation, which was previously considered as representing a hyena den.

Metamorphism of metachert from the Southern Alps, New Zealand. A thermodynamic forward modelling study

<p>Scattered, scarce occurrences of garnet- and quartz-rich metamorphic rock, probably derived from Mn- and Fe-rich chert, occur within metamorphosed greywacke sequences worldwide. The metamorphism of such garnetiferous metacherts has not previously been investigated using modern thermodynamic forward modelling techniques due to the lack of appropriate, internally-consistent activity-composition (a–x) models for Mn-bearing minerals. The present study applies thermodynamic forward modelling using the recently-proposed a–x models of White et al. (2014) to investigate the metamorphism of garnetiferous metachert samples from the Southern Alps, New Zealand. Pressure-temperature (P–T) pseudosections are used in combination with results from petrography, element composition mapping using micro X-ray fluorescence (µXRF) and scanning electron microscope (SEM) methods, and garnet composition data from analytical transects by electron probe microanalysis (EPMA), to study metachert metamorphism. All the samples are compositionally layered, so the possibility exists that an input bulk rock composition might not match the effective bulk composition at the site of garnet growth. If a mineral assemblage stability field in a calculated P–T pseudosection matched the mineral assemblage in the rock, this was taken as an initial indication of a permissible input bulk rock composition. In that case, refined constraints on the P–T conditions were sought by comparing calculated and measured garnet compositions. The studied rocks include samples that are carbonate-bearing, which require consideration of the effects of fluid composition in mixed H₂O–CO₂ fluids, as well as a sample in which the garnet is strongly zoned, texturally-complex, and inferred to be of polymetamorphic origin. The effects of element fractionation by that garnet were investigated by recalculating the P–T pseudosection using a new bulk rock composition with the garnet core content removed. In none of the samples did the calculated and observed composition isopleths for the garnet cores match, suggesting that initial garnet nucleation in these Mn-rich rocks was locally controlled. For most samples in which the calculated and observed mineral assemblages matched, successful estimates of the peak metamorphic conditions were obtained. A garnet chert (A12E) from the mylonite zone of the Alpine Fault at Vine Creek, near Hokitika, gave a tight intersection of composition isopleths, indicating peak metamorphic conditions of 510 °C/5.5 kbar, after recalculation to correct for element fractionation by the strongly-zoned garnet. This tight, modern constraint is within error of previously-reported results from traditional geothermobarometry (420–600 °C/5.9–13 kbar) and Raman spectroscopy of carbonaceous material (RSCM T = 556 °C) from nearby sites. A peak metamorphic estimate of 520–550 °C/7–10 kbar was obtained from a dolomite-bearing sample from the garnet zone near Fox Glacier (J34), in good comparison with published temperatures from Raman spectroscopy of carbonaceous material in nearby metagreywacke samples (526–546 °C). The prograde metamorphic P–T path was probably steep, based on growth of the garnet core at ~475535 °C/5–9 kbar. The successful results for these garnet chert samples show that the new a-x models for Mn-bearing minerals extend the range of rock types that are amenable to pseudosection modelling. Results obtained in this study also serve to highlight several possible concerns: a) garnet nucleation and initial growth in very Mn-rich rocks may be subject to local compositional or kinetic controls; b) bulk rock compositions may not always mimic the effective bulk composition; c) the existing a–x models for Mn-bearing minerals and white micas may need refining; and d) some rocks may simply be ill-suited to thermodynamic forward modelling. Items a) and b) may be indicated by the common observation of a mismatch between predicted and measured garnet composition isopleths for garnet cores, and by a mismatch between garnet composition isopleths and the appropriate mineral assemblage field for sample AMS01, from the mylonite zone, Hari Hari, Southern Alps. For item c) every P–T pseudosection calculated using the new a–x models for Mn-bearing minerals showed garnet stable to very low temperatures below 300 °C. In addition, the P–T pseudosection for an oligoclase-zone metachporphyroblasts of Fe-Ti oxides (magnetitert (Sample J36) from Hari Mare stream, Franz Josef - Fox Glacier, indicated that the white mica margarite should be present instead of plagioclase (oligoclase), for a rock in which oligoclase is present and margarite is absent, a problem previously noted elsewhere. Item d) is exemplified by a very garnet-rich ferruginous metachert sample (J35, garnet zone, headwater region, Moeraki River, South Westland) which proved impossible to model successfully due to its complex mineral growth and deformation history. This sample contained multiple generations of carbonate with differing compositions, amphibole (not incorporated for modelling with the new a–x models for Mn-bearing minerals), large e associated with smaller, possibly later-formed ilmenite), and the garnet bands were offset by late deformation. The garnetiferous metachert samples studied here preserve in their textures and compositions clues to their growth mechanism and metamorphic history. The textures in at least two of the samples are consistent with the diffusion controlled nucleation and growth model for garnet. This research has successfully used state of the art thermodynamic modelling techniques in combination with the latest internally consistent a-x models on Mn-rich metachert, for the first time, extracting P–T conditions of the metamorphism of garnetiferous metachert from the Southern Alps.</p>

Metamorphism of metachert from the Southern Alps, New Zealand. A thermodynamic forward modelling study

<p>Scattered, scarce occurrences of garnet- and quartz-rich metamorphic rock, probably derived from Mn- and Fe-rich chert, occur within metamorphosed greywacke sequences worldwide. The metamorphism of such garnetiferous metacherts has not previously been investigated using modern thermodynamic forward modelling techniques due to the lack of appropriate, internally-consistent activity-composition (a–x) models for Mn-bearing minerals. The present study applies thermodynamic forward modelling using the recently-proposed a–x models of White et al. (2014) to investigate the metamorphism of garnetiferous metachert samples from the Southern Alps, New Zealand. Pressure-temperature (P–T) pseudosections are used in combination with results from petrography, element composition mapping using micro X-ray fluorescence (µXRF) and scanning electron microscope (SEM) methods, and garnet composition data from analytical transects by electron probe microanalysis (EPMA), to study metachert metamorphism. All the samples are compositionally layered, so the possibility exists that an input bulk rock composition might not match the effective bulk composition at the site of garnet growth. If a mineral assemblage stability field in a calculated P–T pseudosection matched the mineral assemblage in the rock, this was taken as an initial indication of a permissible input bulk rock composition. In that case, refined constraints on the P–T conditions were sought by comparing calculated and measured garnet compositions. The studied rocks include samples that are carbonate-bearing, which require consideration of the effects of fluid composition in mixed H₂O–CO₂ fluids, as well as a sample in which the garnet is strongly zoned, texturally-complex, and inferred to be of polymetamorphic origin. The effects of element fractionation by that garnet were investigated by recalculating the P–T pseudosection using a new bulk rock composition with the garnet core content removed. In none of the samples did the calculated and observed composition isopleths for the garnet cores match, suggesting that initial garnet nucleation in these Mn-rich rocks was locally controlled. For most samples in which the calculated and observed mineral assemblages matched, successful estimates of the peak metamorphic conditions were obtained. A garnet chert (A12E) from the mylonite zone of the Alpine Fault at Vine Creek, near Hokitika, gave a tight intersection of composition isopleths, indicating peak metamorphic conditions of 510 °C/5.5 kbar, after recalculation to correct for element fractionation by the strongly-zoned garnet. This tight, modern constraint is within error of previously-reported results from traditional geothermobarometry (420–600 °C/5.9–13 kbar) and Raman spectroscopy of carbonaceous material (RSCM T = 556 °C) from nearby sites. A peak metamorphic estimate of 520–550 °C/7–10 kbar was obtained from a dolomite-bearing sample from the garnet zone near Fox Glacier (J34), in good comparison with published temperatures from Raman spectroscopy of carbonaceous material in nearby metagreywacke samples (526–546 °C). The prograde metamorphic P–T path was probably steep, based on growth of the garnet core at ~475535 °C/5–9 kbar. The successful results for these garnet chert samples show that the new a-x models for Mn-bearing minerals extend the range of rock types that are amenable to pseudosection modelling. Results obtained in this study also serve to highlight several possible concerns: a) garnet nucleation and initial growth in very Mn-rich rocks may be subject to local compositional or kinetic controls; b) bulk rock compositions may not always mimic the effective bulk composition; c) the existing a–x models for Mn-bearing minerals and white micas may need refining; and d) some rocks may simply be ill-suited to thermodynamic forward modelling. Items a) and b) may be indicated by the common observation of a mismatch between predicted and measured garnet composition isopleths for garnet cores, and by a mismatch between garnet composition isopleths and the appropriate mineral assemblage field for sample AMS01, from the mylonite zone, Hari Hari, Southern Alps. For item c) every P–T pseudosection calculated using the new a–x models for Mn-bearing minerals showed garnet stable to very low temperatures below 300 °C. In addition, the P–T pseudosection for an oligoclase-zone metachporphyroblasts of Fe-Ti oxides (magnetitert (Sample J36) from Hari Mare stream, Franz Josef - Fox Glacier, indicated that the white mica margarite should be present instead of plagioclase (oligoclase), for a rock in which oligoclase is present and margarite is absent, a problem previously noted elsewhere. Item d) is exemplified by a very garnet-rich ferruginous metachert sample (J35, garnet zone, headwater region, Moeraki River, South Westland) which proved impossible to model successfully due to its complex mineral growth and deformation history. This sample contained multiple generations of carbonate with differing compositions, amphibole (not incorporated for modelling with the new a–x models for Mn-bearing minerals), large e associated with smaller, possibly later-formed ilmenite), and the garnet bands were offset by late deformation. The garnetiferous metachert samples studied here preserve in their textures and compositions clues to their growth mechanism and metamorphic history. The textures in at least two of the samples are consistent with the diffusion controlled nucleation and growth model for garnet. This research has successfully used state of the art thermodynamic modelling techniques in combination with the latest internally consistent a-x models on Mn-rich metachert, for the first time, extracting P–T conditions of the metamorphism of garnetiferous metachert from the Southern Alps.</p>

Petrogenesis of Neoproterozoic Ultramafic Rocks, Wadi Ibib–Wadi Shani, South Eastern Desert, Egypt: Constraints from Whole Rock and Mineral Chemistry

This contribution deals with new geology, petrography, and bulk-rock/mineral chemistry of the poorly studied ultramafics of Wadi Ibib–Wadi Shani (WI–WS) district, South Eastern Desert, Egypt. These ultramafics are dismembered ophiolitic rocks that can be subdivided into serpentinites and serpentinized peridotites. Primary minerals such as olivine and pyroxene are absent in serpentinites, but relics of them occur in serpentinized peridotites. Pseudomorph after olivine is indicated by common hourglass textures with less mesh, whilst schistose bastites reflect a pyroxene pseudomorph. Chromite can be subdivided into Cr-spinel and Al-spinel. Cr-spinel ranges from chromite to magnesochromite in composition, whereas Al-spinel belongs to the spinel field. Cr-spinel includes YCr (Cr/(Cr+Al+Fe+3), YAl (Al/(Al+Cr+Fe+3), and YFe+3 (Fe+3/(Fe+3+Al+Cr), similar to forearc peridotite, whilst Al-spinel is more similar to abyssal peridotite, and may be formed during inanition of subduction processes in proto forearc environments. The main secondary minerals are tremolite, talc, and chlorite—which is subdivided into pycnochlorite and diabantite—and their temperature ranges from 174 to 224 °C. The examined rocks had undergone high partial melting degrees (>25%), as indicated by the Cr# of their unaltered cores (Cr-spinel, >0.6), whole rocks (Al2O3, SiO2, CaO, and MgO), trace and REEs, depleted Na2O, Al2O3, and Cr2O3 of clinopyroxene, and high forsterite content ((Fo = 100 Mg/Mg + Fe), av. 95.23 mol%), consistent with forearc settings.

Geology and Petrology of Ruapehu Volcano and Related Vents

<p>Ruapehu Volcano is an active, multiple-vent, andesite composite volcano at the southern terminus of the Taupo Volcanic Zone, central North Island, New Zealand. The present-day volume of Ruapehu is estimated at 110 km3, and construction of the massif probably occurred during the past 0.5 m.y. Geologic mapping and stratigraphic studies have led to the recognition of four periods of cone construction, each occurring over 104-105 year time intervals. On the basis of lithologic/petrographic differences, and conspicuous unconformities which separate the deposits of each cone-building period, four new formations are defined, comprises the Ruapehu Group. Te Herenga formation (new formation name) comprises the oldest deposits of Ruapehu (upper lavas ca. 0.23 Ma) and is exposed as planeze surfaces and aretes on N and NW Ruapehu. The formation includes lava flows, tuff breccias, and small intrusive bodies surrounded by zones of hydrothermal alteration. There is little petrographic and compositional diversity; most lavas are porphyritic titanomagnetite- augite- hypersthene- plagioclase basic andesites. Wahiance Formation (new formation name) is younger than Te Herenga Fm,. but of unknown age. It is well exposed on SE Ruapehu, and comprises mostly lava flows and tuff breccias. The lavas comprise acid and basic andesites. Mangewhero Formation (new formation name) is well exposed everywhere except SE Ruapehu, and the upper lavas and pyroclastics (ca. 0.02 Ma) form the present high peeks and main cone of Ruapehu. The lavas are petrographically and geochemically diverse, ranging from basalt to decite in bulk composition. Some of the lower lavas are olivine-beering andesites of hybrid orgin. Whakapapa Formation (new formation name; ca 15,000 years to present) comprises conspicuously young lava flows, tuff breccias, airfall pyroclastics and minor pyroclastic flows of acid- and basic andesite. The deposits of these post-glacial summit and flank eruptions are subdivided into the lwikau, Rangataua, Tama and Crater Lake Members. 'Related vents' produced Heuhungatahi Andesite Fm. (> 0.5 Ma?), and Holocene deposits of basalt and basic andesite at isolated, monogenetic centres comprising Ohakune Andesite Fm., Pukeonake Andesite Fm., and Waimarino Basalt Fm. (new formation name). Most Ruapehu lavas are medium-K acid and basic andesites (mean of 144 bulk rock analyses is 57.8 wt % SiO2), but rare basalt and minor decite are present. Nearly all lavas are porphyritic in plagioclase, augite and hypersthene [plus or minus] olivine, with titanomagnetite micro- phenocrysts, and contain abundant metamorphic and igneous rock inclusions. Petrography, mineral chemistry and bulk rock chemistry indicate fractional crystallization series from parental basalts (52-53 % SiO2, Q-normative, low-alumina) to medium-K basic- and acid andesites (58-59 % SiO2). Early fractionating minerals are olivine and clinopyroxene with minor chrome spinel and plagioclase, followed by plagioclase, orthopyroxene, clinopyroxene and minor titanomagnetite in later stages of differentiation. Thus, basalt differentiation to produce andesites involves 'POAM-type' (Gill, 1981) fractional crystallization. Three second-order differentiation processes operate concurrently with frational crystallization: (1) Crystal accumulation involves addition of co-genetic plutonic rock fragments and crystals derived from them. These inclusions are common and few rocks represent liquid compositions. (2) Magma mixing involves mingling of magmas in repeatedly-occupied conduits. End members are as diverse as basalt and decite, yielding petrogaphically and chemically distinctive high-Mg andesites of the upper cone complex and parasitic centres. (3) Selective crustal assimilation is suggested by partially fused metamorphic inclusions, positive correlation of 87Sr/86Sr with SiO2, and failure of simple 'POAM' fractionation to explain decites (63-65 % SiO2). Petrogenesis of Ruapehu andesites takes place under open-system condition, involving production of parental Q-normative basalts in the mantle wedge, concurrent fractional crystallization and crustal contamination, entrainment of co-genetic plutonic rocks, and mixing of magmas in common conduits.</p>