Petrology of some oceanic island basalts: PRIMELT2.XLS software for primary magma calculation

PRIMELT2.XLS software is introduced for calculating primary magma composition and mantle potential temperature (TP) from an observed lava composition. It is an upgrade over a previous version in that it includes garnet peridotite melting and it detects complexities that can lead to overestimates in TP by >100°C. These are variations in source lithology, source volatile content, source oxidation state, and clinopyroxene fractionation. Nevertheless, application of PRIMELT2.XLS to lavas from a wide range of oceanic islands reveals no evidence that volatile-enrichment and source fertility are sufficient to produce them. All are associated with thermal anomalies, and this appears to be a prerequisite for their formation. For the ocean islands considered in this work, TP maxima are typically ~1450–1500°C in the Atlantic and 1500–1600°C in the Pacific, substantially greater than ~1350°C for ambient mantle. Lavas from the Galápagos Islands and Hawaii record in their geochemistry high TP maxima and large ranges in both TP and melt fraction over short horizontal distances, a result that is predicted by the mantle plume model

The Significance of Multiple Saturation Points in the Context of Polybaric Near-fractional Melting

Experimental petrologists have successfully located basaltic liquid compositions parental to mid-ocean ridge basalt that are, within experimental resolution, multiply saturated with three-phase harzburgite or four-phase lherzolite assemblages on their liquidus at some elevated pressure. Such an experimental result is a necessary consequence of any paradigm in which erupted basalts derive from single-batch primary liquids that equilibrate with a mantle residue and undergo no subsequent magma mixing before differentiation and eruption. Here we investigate whether, conversely, such evidence of multiple saturation is sufficient to exclude dynamic melting models wherein increments of melt are mixed after segregation from residues, during melt transport or in magma chambers. Using two independent models of crystal–liquid equilibria to simulate polybaric near-fractional peridotite melting, we find that aggregate liquids from such melting processes can display near-intersections of liquidus surfaces too close to distinguish experimentally from exact multiple saturation points. Given uncertainties in glass compositions, fractionation corrections, experimental temperature and pressure conditions, and achievement of equilibrium, these results suggest that polybaric mixtures can in fact masquerade as mantle-equilibrated single-batch primary liquids. Multiple saturation points on the liquidus surfaces of primitive basalts do, however, preserve information about the average pressure of extraction of their constituent increments of liquid

Fluid Outflows From Venus Impact Craters: Analysis From Magellan Data

Many impact craters on Venus have unusual outflow features originating in or under the continuous ejecta blankets and continuing downhill into the surrounding terrain. These features clearly resulted from flow of low-viscosity fluids, but the identity of those fluids is not clear. In particular, it should not be assumed a priori that the fluid is an impact melt. A number of candidate processes by which impact events might generate the observed features are considered, and predictions are made concerning the rheological character of flows produced by each mechanism. A sample of outflows was analyzed using Magellan images and a model of unconstrained Bingham plastic flow on inclined planes, leading to estimates of viscosity and yield strength for the flow materials. It is argued that at least two different mechanisms have produced outflows on Venus: an erosive, channel-forming process and a depositional process. The erosive fluid is probably an impact melt, but the depositional fluid may consist of fluidized solid debris, vaporized material, and/or melt

Adiabat_1ph: A new public front-end to the MELTS, pMELTS, and pHMELTS models

The program adiabat_1ph is a simple text-menu driver for subroutine versions of the algorithms MELTS, pMELTS, and pHMELTS [Asimow et al., 2004; Ghiorso et al., 2002; Ghiorso and Sack, 1995]. It may be used to calculate equilibrium assemblages along a thermodynamic path set by the user and can simultaneously calculate trace element distributions. The MELTS family of algorithms is suitable for multicomponent systems, which may be anhydrous, water-undersaturated, or water-saturated, with the options of buffering oxygen fugacity and/or water activity. A wide variety of calculations can be performed either subsolidus or with liquid(s) present; melting and crystallization may be batch, fractional, or continuous. The software is suitable for Linux, MacOS X, and Windows, and many aspects of program execution are controlled by environment variables. Perl scripts are also provided that may be used to invoke adiabat_1ph with some command line options and to produce output that may be easily imported into spreadsheet programs, such as Microsoft Excel. Benefits include a batch mode, which allows almost complete automation of the calculation process when suitable input files are written. This technical brief describes version 1.04, which is provided as ancillary material. Binaries, scripts, documentation, and example files for this and future releases may be downloaded at http://www.gps.caltech.edu/~asimow/adiabat. On a networked computer, adiabat_1ph automatically checks whether a newer version is available

Tests of random density models of terrestrial planets

Random density models are analyzed to determine the low degree harmonics of the gravity field of a planet, and therefrom two properties: an axiality P_l , the percent of the degree variance in the zonal term referred to an axis through the maximum for degree l; and an angularity E_(ln) , the angle between the maxima for two degrees l, n. The random density distributions give solutions reasonably consistent with the axialities and angularities for the low degrees, l < 5, of Earth, Venus, and Moon, but not for Mars, which has improbably large axialities and small angularities. Hence the random density model is an unreliable predictor for the non‐hydrostatic second‐degree gravity of Mars, and thus for the moment‐of‐inertia, which is more plausibly close to 0.365MR^2

A lesson in defining "extinct"

Toni Feder’s Issues and Events news item about the High-Altitude Water Cherenkov Gamma-Ray Observatory under construction in Mexico ( Physics Today, October 2013, page 22) begins by describing the site as “nestled at 4100 m on the slopes of Sierra Negra, an extinct volcano.” However, as the picture accompanying her piece or a cursory examination in Google Earth makes clear, the site is not on the slopes of Sierra Negra. It sits in the saddle between Sierra Negra and its much larger companion Pico de Orizaba, also known as Citlaltépetl. Sierra Negra, a minor flank cone of that larger volcanic system, may now be extinct, but Pico de Orizaba (http://www.volcano.si.edu/volcano.cfm?vnum=1401-10=) absolutely is not. It is dormant at best, having erupted as recently as 1846. To compound the issue, the photograph shows that the observatory site is built squarely in front of an obviously young lava flow. Now, I’m sure that site selection was done with due attention to natural hazards, but as long as astronomers insist on building expensive observatories on top of volcanoes, there needs to be clear understanding and common vocabulary between astronomers and geologists. “Extinct” is a troublesome word when applied to a volcano

PetroPlot: A plotting and data management tool set for Microsoft Excel

PetroPlot is a 4000-line software code written in Visual Basic for the spreadsheet program Excel that automates plotting and data management tasks for large amount of data. The major plotting functions include: automation of large numbers of multiseries XY plots; normalized diagrams (e.g., spider diagrams); replotting of any complex formatted diagram with multiple series for any other axis parameters; addition of customized labels for individual data points; and labeling flexible log scale axes. Other functions include: assignment of groups for samples based on multiple customized criteria; removal of nonnumeric values; calculation of averages/standard deviations; calculation of correlation matrices; deletion of nonconsecutive rows; and compilation of multiple rows of data for a single sample to single rows appropriate for plotting. A cubic spline function permits curve fitting to complex time series, and comparison of data to the fits. For users of Excel, PetroPlot increases efficiency of data manipulation and visualization by orders of magnitude and allows exploration of large data sets that would not be possible making plots individually. The source codes are open to all users

Grain boundary partitioning of Ar and He

An experimental procedure has been developed that permits measurement of the partitioning of Ar and He between crystal interiors and the intergranular medium (ITM) that surrounds them in synthetic melt-free polycrystalline diopside aggregates. ^(37)Ar and ^(4)He are introduced into the samples via neutron irradiation. As samples are crystallized under sub-solidus conditions from a pure diopside glass in a piston cylinder apparatus, noble gases diffusively equilibrate between the evolving crystal and intergranular reservoirs. After equilibration, ITM Ar and He is distinguished from that incorporated within the crystals by means of step heating analysis. An apparent equilibrium state (i.e., constant partitioning) is reached after about 20 h in the 1450 °C experiments. Data for longer durations show a systematic trend of decreasing ITM Ar (and He) with decreasing grain boundary (GB) interfacial area as would be predicted for partitioning controlled by the network of planar grain boundaries (as opposed to ITM gases distributed in discrete micro-bubbles or melt). These data yield values of GB-area-normalized partitioning, K¯^(Ar)_(ITM), with units of (Ar/m^3 of solid)/(Ar/m^2 of GB) of 6.8 x 10^3 – 2.4 x 104 m^(-1). Combined petrographic microscope, SEM, and limited TEM observation showed no evidence that a residual glass phase or grain boundary micro-bubbles dominated the ITM, though they may represent minor components. If a nominal GB thickness (δ) is assumed, and if the density of crystals and the grain boundaries are assumed equal, then a true grain boundary partition coefficient (K^(Ar)_(GB) = X^(Ar)_(crystals)/X^(Ar)_(GB) may be determined. For reasonable values of δ, K^(Ar)_(GB) is at least an order of magnitude lower than the Ar partition coefficient between diopside and melt. Helium partitioning data provide a less robust constraint with K¯^(He)_(ITM) between 4 x 10^3 and 4 x 10^4 cm^(-1), similar to the Ar partitioning data. These data suggest that an ITM consisting of nominally melt free, bubble free, tight grain boundaries can constitute a significant but not infinite reservoir, and therefore bulk transport pathway, for noble gases in fine grained portions of the crust and mantle where aqueous or melt fluids are non-wetting and of very low abundance (i.e., <0.1% fluid). Heterogeneities in grain size within dry equilibrated systems will correspond to significant differences in bulk rock noble gas content

A hydrous melting and fractionation model for mid-ocean ridge basalts: Application to the Mid-Atlantic Ridge near the Azores

The major element, trace element, and isotopic composition of mid-ocean ridge basalt glasses affected by the Azores hotspot are strongly correlated with H2O content of the glass. Distinguishing the relative importance of source chemistry and potential temperature in ridge-hotspot interaction therefore requires a comprehensive model that accounts for the effect of H2O in the source on melting behavior and for the effect of H2O in primitive liquids on the fractionation path. We develop such a model by coupling the latest version of the MELTS algorithm to a model for partitioning of water among silicate melts and nominally anhydrous minerals. We find that much of the variation in all major oxides except TiO2 and a significant fraction of the crustal thickness anomaly at the Azores platform are explained by the combined effects on melting and fractionation of up to ~700 ppm H2O in the source with only a small thermal anomaly, particularly if there is a small component of buoyantly driven active flow associated with the more H2O-rich melting regimes. An on-axis thermal anomaly of ~35°C in potential temperature explains the full crustal thickness increase of ~4 km approaching the Azores platform, whereas a ≥75°C thermal anomaly would be required in the absence of water or active flow. The polybaric hydrous melting and fractionation model allows us to solve for the TiO2, trace element and isotopic composition of the H2O-rich component in a way that self-consistently accounts for the changes in the melting and fractionation regimes resulting from enrichment, although the presence and concentration in the enriched component of elements more compatible than Dy cannot be resolved