Chalcophile element geochemistry of the Boggy Plain Zoned Pluton, southeastern Australia : a S-saturated barren compositionally diverse magma

The behavior of the platinum group elements (PGE) and Re in felsic magmas is poorly understood due to scarcity of data. We report the concentrations of Ni, Cu, Re, and PGE in the compositionally diverse Boggy Plain zoned pluton (BPZP), which shows a variation of rock type from gabbro through granodiorite and granite to aplite with a SiO2 range from 52 to 74 wt %. In addition, major silicate and oxide minerals were analyzed for Ni, Cu, and Re, and a systematic sulfide study was carried out to investigate the role of silicate, oxide, and sulfide minerals on chalcophile element geochemistry of the BPZP. Mass balance calculation shows that the whole rock Cu budget hosted by silicate and oxide minerals is <13 wt % and that Cu is dominantly located in sulfide phases, whereas most of the whole rock Ni budget (>70 wt %) is held in major silicate and oxide minerals. Rhenium is dominantly hosted by magnetite and ilmenite. Ovoid-shaped sulfide blebs occur at the boundary between pyroxene phenocrysts and neighboring interstitial phases or within interstitial minerals in the gabbro and the granodiorite. The blebs are composed of pyrrhotite, pyrite, chalcopyrite, and S-bearing Fe-oxide, which contain total trace metals (Co, Ni, Cu, Ag, Pb) up to ~16 wt %. The mineral assemblage, occurrence, shape, and composition of the sulfide blebs are a typical of magmatic sulfide. PGE concentrations in the BPZP vary by more than two orders of magnitude from gabbro (2.7–7.8 ppb Pd, 0.025–0.116 ppb Ir) to aplite (0.05 ppb Pd, 0.001 ppb Ir). Nickel, Cu, Re, and PGE concentrations are positively correlated with MgO in all the rock types although there is a clear discontinuity between the granodiorite and the granite in the trends for Ni, Rh, and Ir when plotted against MgO. Cu/Pd values gradually increase from 6,100 to 52,600 as the MgO content decreases. The sulfide petrology and chalcophile element geochemistry of the BPZP show that sulfide saturation occurred in the late gabbroic stage of magma differentiation. Segregation and distribution of these sulfide blebs controlled Cu and PGE variations within the BPZP rocks although the magma of each rock type may have experienced a different magma evolution history in terms of crustal assimilation and crystal fractionation. The sulfide melt locked in the cumulate rocks must have sequestered a significant portion of the chalcophile elements, which restricted the availability of these metals to magmatic-hydrothermal ore fluids. Therefore, we suggest that the roof rocks that overlay the BPZP were not prospective for magmatic-hydrothermal Cu, Au, or Cu–Au deposits.

Generation and evolution of magma beneath the East Pacific Rise: Constraints from U-series disequilibrium and plagioclase-hosted melt inclusions

Mid-ocean ridge basalt (MORB) samples from the East Pacific Rise (EPR 12 degrees 50'N) were analyzed for U-series isotopes and compositions of plagioclase-hosted melt inclusions. The Ra-226 and Th-230 excesses are negatively correlated; the Ra-226 excess is positively correlated with Mg# and Sm/Nd, and is negatively correlated with La/Sm and Fe-8; the Th-230 excess is positively correlated with Fe-8 and La/Sm and is negatively correlated with Mg# and Sm/Nd. Interpretation of these correlations is critical for understanding the magmatic process. There are two models (the dynamic model and the "two-porosity" model) for interpreting these correlations, however, some crucial parameters used in these models are not ascertained. We propose instead a model to explain the U-series isotopic compositions based on the control of melt density variation. For melting either peridotite or the "marble-cake" mantle, the FeOt content, Th-230 excess and La/Sm ratio increases and Sm/Nd decreases with increasing pressure. A deep melt will evolve to a higher density and lower Mg# than a shallow melt, the former corresponds to a long residence time, which lowers the Ra-226 excess significantly. This model is supported by the existence of low Ra-226 excesses and high Th-230 excesses in MORBs having a high Fe-8 content and high density. The positive correlation of Ra-226 excess and magma liquidus temperature implies that the shallow melt is cooled less than the deep melt due to its low density and short residence time. The correlations among Fe-8, Ti-8 and Ca-8/Al-8 in plagioclase-hosted melt inclusions further prove that MORBs are formed from melts having a negative correlation in melting depths and degrees. The negative correlation of Ra-226 excess vs. chemical diversity index (standard deviation of Fe-8, Ti-8 and Ca-8/Al-8) of the melt inclusions is in accordance with the influence of a density-controlled magma residence time. We conclude that the magma density variation exerts significant control on residence time and U-series isotopic compositions. (c) 2010 Elsevier B.V. All rights reserved.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Petrogenesis of Basaltic Andesites at Tofua Volcano

Tofua volcano is situated midway along the Tonga oceanic arc and has undergone two phases of ignimbrite-forming activity. The eruptive products are almost entirely basaltic andesites (52 center dot 5-57 wt % SiO2) with the exception of a volumetrically minor pre-caldera dacite. The suite displays a strong tholeiitic trend with K2O <1 wt %. Phenocryst assemblages typically comprise plagioclase + clinopyroxene +/- orthopyroxene with microlites of Ti-magnetite. Olivine (Fo(83-88)) is rare and believed to be dominantly antecrystic. An increase in the extent and frequency of reverse zoning in phenocrysts, sieve-textured plagioclase and the occurrence of antecrystic phases in post-caldera lavas record a shift to dynamic conditions, allowing the interaction of magma batches that were previously distinct. Pyroxene thermobarometry suggests crystallization at 950-1200 degrees C and 0 center dot 8-1 center dot 8 kbar. Volatile measurements of glassy melt inclusions indicate a maximum H2O content of 4 center dot 16 wt % H2O, and CO2-H2O saturation curves indicate that crystallization occurred at two levels, at depths of 4-5 center dot 5 km and 1 center dot 5-2 center dot 5 km. Major and trace element models suggest that the compositions of the majority of the samples represent a differentiation trend whereby the dacite was produced by 65% fractional crystallization of the most primitive basaltic andesite. Trace element models suggest that the sub-arc mantle source is the residuum of depleted Indian mid-ocean ridge basalt mantle (IDMM-1% melt), whereas radiogenic isotope data imply addition of 0 center dot 2% average Tongan sediment melt and a fluid component derived from the subducted altered Pacific oceanic crust. A horizontal array on the U-Th equiline diagram and Ra excesses of up to 500% suggest fluid addition to the mantle wedge within the last few thousand years. Time-integrated (Ra-226/Th-230) vs Sr/Th and Ba/Th fractionation models imply differentiation timescales of up to 4500 years for the dacitic magma compositions at Tofua.

Magma Evolution in the Primitive, Intra-oceanic Tonga Arc: Rapid Petrogenesis of Dacites at Fonualei Volcano

Fonualei is unusual amongst subaerial volcanoes in the Tonga arc because it has erupted dacitic vesicular lavas, tuffs and phreomagmatic deposits for the last 165 years. The total volume of dacite may approach 5 km(3) and overlies basal basaltic andesite and andesite lavas that are constrained to be less than a few millennia in age. All of the products are crystal-poor and formed from relatively low-viscosity magmas inferred to have had temperatures of 1100-1000 degrees C, 2-4 wt % H2O and oxygen fugacities 1-2 log units above the quartz-fayalite-magnetite buffer. Major and trace element data, along with Sr-Nd-Pb and U-Th-Ra isotope data, are used to assess competing models for the origin of the dacites. Positive correlations between Sc and Zr and Sr rule out evolution of the within-dacite compositional array by closed-system crystal fractionation of a single magma batch. An origin by partial melting of lower crustal amphibolites cannot reproduce these data trends or, arguably, any of the dacites either. Instead, we develop a model in which the dacites reflect mixing between two dacitic magmas, each the product of fractional crystallization of basaltic andesite magmas formed by different degrees of partial melting. Mixing was efficient because the two magmas had similar temperatures and viscosities. This is inferred to have occurred at shallow (2-6 km) depths beneath the volcano. U-Th-Ra disequilibria in the basaltic andesite and andesite indicate that the parental magmas had fluids added to their mantle source regions less than 8 kyr ago and that fractionation to the dacitic compositions took less than a few millennia. The 165 year eruption period for the dacites implies that mixing occurred on a similar timescale, possibly during ascent in conduits. The composition of the dacites renders them unsuitable candidates as contributors to average continental crust.

Plume-lithosphere interaction at Santiago Island (Cape Verde) : implications for ocean island magma genesis models

Tese de doutoramento, Geologia (Geoquímica), Universidade de Lisboa, Faculdade de Ciências, 2014

The influence of viscous and latent heating on crystal-rich magma flow in a conduit

The flow dynamics of crystal-rich high-viscosity magma is likely to be strongly influenced by viscous and latent heat release. Viscous heating is observed to play an important role in the dynamics of fluids with temperature-dependent viscosities. The growth of microlite crystals and the accompanying release of latent heat should play a similar role in raising fluid temperatures. Earlier models of viscous heating in magmas have shown the potential for unstable (thermal runaway) flow as described by a Gruntfest number, using an Arrhenius temperature dependence for the viscosity, but have not considered crystal growth or latent heating. We present a theoretical model for magma flow in an axisymmetric conduit and consider both heating effects using Finite Element Method techniques. We consider a constant mass flux in a 1-D infinitesimal conduit segment with isothermal and adiabatic boundary conditions and Newtonian and non-Newtonian magma flow properties. We find that the growth of crystals acts to stabilize the flow field and make the magma less likely to experience a thermal runaway. The additional heating influences crystal growth and can counteract supercooling from degassing-induced crystallization and drive the residual melt composition back towards the liquidus temperature. We illustrate the models with results generated using parameters appropriate for the andesite lava dome-forming eruption at Soufriere Hills Volcano, Montserrat. These results emphasize the radial variability of the magma. Both viscous and latent heating effects are shown to be capable of playing a significant role in the eruption dynamics of Soufriere Hills Volcano. Latent heating is a factor in the top two kilometres of the conduit and may be responsible for relatively short-term (days) transients. Viscous heating is less restricted spatially, but because thermal runaway requires periods of hundreds of days to be achieved, the process is likely to be interrupted. Our models show that thermal evolution of the conduit walls could lead to an increase in the effective diameter of flow and an increase in flux at constant magma pressure.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

Tracing pre-eruptive magma degassing using (Pb-210/Ra-226) disequilibria in the volcanic deposits of the 1980-1986 eruption of Mount St. Helens

Disequilibria between Pb-210 and Ra-226 can be used to trace magma degassing, because the intermediate nuclides, particularly Rn-222, are volatile. Products of the 1980-1986 eruptions of Mount St. Helens have been analysed for (Pb-210/Ra-226). Both excesses and deficits of Pb-210 are encountered suggesting rapid gas transfer. The time scale of diffuse, non-eruptive gas escape prior to 1980 as documented by Pb-210 deficits is on the order of a decade using the model developed by Gauthier and Condomines (Earth Planet. Sci. Lett. 172 (1999) 111-126) for a non-renewed magma chamber and efficient Rn removal. The time required to build-up Pb-210 excess is much shorter (months) as can be observed from steady increases of (Pb-210/Ra-226) with time during 1980-1982. The formation of Pb-210 excess requires both rapid gas transport through the magma and periodic blocking of gas escape routes. Superposed on this time trend is the natural variability of (Pb-210/Ra-226) in a single eruption caused by tapping magma from various depths. The two time scales of gas transport, to create both Pb-210 deficits and Pb-210 excesses, cannot be reconciled in a single event. Rather Pb-210 deficits are associated with pre-eruptive diffuse degassing, while Pb-210 excesses document the more vigorous degassing associated with eruption and recharge of the system. (c) 2006 Elsevier B.V. All rights reserved.

Pb-210-Ra-226 and Ra-228-Th-232 systematics in young arc lavas: implications for magma degassing and ascent rates

New data show that island arc rocks have (Pb-210/Ra-226)(o) ratios which range from as low as 0.24 up to 2.88. In contrast, (Ra-22S/Th-232) appears always within error of I suggesting that the large Ra-226-excesses observed in arc rocks were generated more than 30 years ago. This places a maximum estimate on melt ascent velocities of around 4000 m/year and provides further confidence that the Ra-226 excesses reflect deep (source) processes rather than shallow level alteration or seawater contamination. Conversely, partial melting must have occurred more than 30 years prior to eruption. The Pb-210 deficits are most readily explained by protracted magma degassing. Using published numerical models, the data suggest that degassing occurred continuously for periods up to several decades just prior to eruption but no link with eruption periodicity was found. Longer periods are required if degassing is discontinuous, less than 100% efficient or if magma is recharged or stored after degassing. The long durations suggest much of this degassing occurs at depth with implications for the formation of hydrothermal and copper-porphyry systems. A suite of lavas erupted in 1985-1986 from Sangeang Api volcano in the Sunda arc are characterised by deficits of Pb-210 relative to Ra-226 from which 6-8 years of continuous Rn-222 degassing would be inferred from recent numerical models. These data also form a linear (Pb-210)/Pb-(Ra-226)/Pb array which might be interpreted as a 71-year isochron. However, the array passes through the origin suggesting displacement downwards from the equiline in response to degassing and so the slope of the array is inferred not to have any age significance. Simple modelling shows that the range of (Ra-226)/Pb ratios requires thousands of years to develop consistent with differentiation occurring in response to cooling at the base of the crust. Thus, degassing post-dated, and was not responsible for magma differentiation. The formation, migration and extraction of gas bubbles must be extremely efficient in mafic magma whereas the higher viscosity of more siliceous magmas retards the process and can lead to Pb-210 excesses. A possible negative correlation between (Pb-210/Ra-226)(o) and SO2 emission rate requires further testing but may have implications for future eruptions. (C) 2004 Elsevier B.V. All rights reserved.

The magma budget of Volcan Arenal, Costa Rica from 1980 to 2004

Using topographic data collected by radar interferometry, stereo-photogrammetry, and field survey we have measured the changing surface of Volcan Arenal in Costa Rica over the period from 1980 to 2004. During this time this young volcano has mainly effused basaltic andesite lava, continuing the activity that began in 1968. Explosive products form only a few percent of the volumetric output. We have calculated digital elevation models for the years 1961, 1988 and 1997 and modified existing models for 2000 and 2004. From these we have estimated the volume of lava effused and coupled this with the data presented by an earlier study for 1968-1980. We find that a dense rock equivalent volume of 551 M m(3) was effused from 1968 to 2004. The dense rock equivalent effusion rate fell from about 2 m(3) s(-1) to about 0.1-0.2 m(3) s(-1) over the same period, with an average rate of about 0.5 m(3) s(-1). Between 1980 and 2004, the average effusion rate was 0.36 m(3) s(-1), a similar rate to that measured between 1974 and 1980. There have been two significant deviations from this long-term rate. The effusion rate increased from 1984 to 1991, at the same time as explosivity increased. After a period of moderate effusion rates in the 1990s, the rate fell to lower levels around 1999. (c) 2006 Elsevier B.V. All rights reserved.

The Sphene-Centered Ocellar Texture: An Effect of Grain-Supported Flow and Melt Migration in a Hyperdense Magma Mush

The sphene-centered ocellar texture consists of leucocratic ocelli with sphene (titanite) crystals at the center, enclosed in a biotite-rich matrix. This texture has been recognized worldwide in hybrid intermediate rocks. On the basis of structural, petrological, and geochronological data from selected outcrops of the Variscan Ribadelago pluton (NW Iberian Massif), we propose that the ocelli were formed by migration and accumulation of a residual melt through a plagioclase- and biotite-dominated crystalline framework. At the late stage of crystallization, the magma acted as a hyperdense suspension and reacted to the pressure gradient caused by the regional stress field, entering the domain of grain-supported flow. Microstructures reveal that aligned crystal domains arose in the crystal framework from the shearing and compaction of the crystal mush and behaved as magmatic microshears. Relative displacement of adjacent crystal clusters along these microshears corresponded to the onset of Reynolds dilatancy that generated an expansion of the crystal mush, involving melt migration and pore aperture. The mineralogy of the ocelli, dominated by andesine and sphene, represents the composition of the migrating melt. The chemistry of this late, Ti-rich melt stems from the incongruent melting of biotite. Magmatic sphene from the ocelli yields a U-Pb age of 317 +/- 1 Ma, which represents the final crystallization of the hybridized magmatic system. Moreover, this texture offers an opportunity to better understand the rheological behavior of highly crystallized magmas.