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

Mantle Melting Beneath Mid-Ocean Ridges

none availableEarth and Planetary Science

Feedback Between Deglaciation, Volcanism, and Atmospheric CO2

An evaluation of the historical record of volcanic eruptions shows that subaerial volcanism increases globally by two to six times above background levels between 12 ka and 7 ka, during the last deglaciation. Increased volcanism occurs in deglaciating regions. Causal mechanisms could include an increase in magma production owing to the mantle decompression caused by ablation of glaciers and ice caps or a more general pacing of when eruptions occur by the glacial variability. A corollary is that ocean ridge volcanic production should decrease with the rising sea level during deglaciation, with the greatest effect at slow spreading ridges. CO2 output from the increased subaerial volcanism appears large enough to influence glacial/interglacial CO2 variations. We estimate subaerial emissions during deglaciation to be between 1000 and 5000 Gt of CO2 above the long term average background flux, assuming that emissions are proportional to the frequency of eruptions. After accounting for equilibration with the ocean, this additional CO2 flux is consistent in timing and magnitude with ice core observations of a 40 ppm increase in atmospheric CO2 concentration during the second half of the last deglaciation. Estimated decreases in CO2 output from ocean ridge volcanoes compensate for only 20% of the increased subaerial flux. If such a large volcanic output of CO2 occurs, then volcanism forges a positive feedback between glacial variability and atmospheric CO2 concentrations: deglaciation increases volcanic eruptions, raises atmospheric CO2, and causes more deglaciation. Such a positive feedback may contribute to the rapid passage from glacial to interglacial periods. Conversely, waning volcanic activity during an interglacial could lead to a reduction in CO2 and the onset of an ice age. Whereas glacial/interglacial variations in CO2 are generally attributed to oceanic mechanisms, it is suggested that the vast carbon reservoirs associated with the solid Earth may also play an important role.Earth and Planetary Science

Origins of chemical diversity of back-arc basin basalts: a segment-scale study of the Eastern Lau Spreading Center

We report major, trace, and volatile element data on basaltic glasses from the northernmost segment of the Eastern Lau Spreading Center (ELSC1) in the Lau back-arc basin to further test and constrain models of back-arc volcanism. The zero-age samples come from 47 precisely collected stations from an 85 km length spreading center. The chemical data covary similarly to other back-arc systems but with tighter correlations and well-developed spatial systematics. We confirm a correlation between volatile content and apparent extent of melting of the mantle source but also show that the data cannot be reproduced by the model of isobaric addition of water that has been broadly applied to back-arc basins. The new data also confirm that there is no relationship between mantle temperature and the wet melting productivity. Two distinct magmatic provinces can be identified along the ELSC1 axis, a southern province influenced by a “wet component” with strong affinities to arc volcanism and a northern province influenced by a “damp component” intermediate between enriched mid-ocean ridge basalts (E-MORB) and arc basalts. High–field strength elements and rare earth elements are all mobilized to some extent by the wet component, and the detailed composition of this component is determined. It differs in significant ways from the Mariana component reported by E. Stolper and S. Newman (1994), particularly by having lower abundances of most elements relative to H_(2)O. The differences can be explained if the slab temperature is higher for the Mariana and the source from which the fluid is derived is more enriched. The ELSC1 damp component is best explained by mixing between the wet component and an E-MORB-like component. We propose that mixing between water-rich fluids and low-degree silicate melts occurs at depth in the subduction zone to generate the chemical diversity of the ELSC1 subduction components. These modified sources then rise independently to the surface and melt, and these melts mix with melts of the background mantle from the ridge melting regime to generate the linear data arrays characteristic of back-arc basalts. The major and trace element framework for ELSC1, combined with different slab temperatures and compositions for difference convergent margins, may be able to be applied to other back-arc basins around the globe

Glacial cycles drive variations in the production of oceanic crust

Glacial cycles redistribute water between oceans and continents causing pressure changes in the upper mantle, with consequences for melting of Earth's interior. Using Plio-Pleistocene sea-level variations as a forcing function, theoretical models of mid-ocean ridge dynamics that include melt transport predict temporal variations in crustal thickness of hundreds of meters. New bathymetry from the Australian-Antarctic ridge shows significant spectral energy near the Milankovitch periods of 23, 41, and 100 ky, consistent with model predictions. These results suggest that abyssal hills, one of the most common bathymetric features on Earth, record the magmatic response to changes in sea level. The models and data support a link between glacial cycles at the surface and mantle melting at depth, recorded in the bathymetric fabric of the sea floor.Comment: 30 pages, 6 figures (including supplementary information). Resubmitted to Science on 12 December 201

A New Parameterization of Hydrous Mantle Melting

Modeling of melt formation and transport in all tectonic settings requires the inclusion of water, since water has large effects on mantle solidi as well as physical properties of liquids. To facilitate the inclusion of water in melting models this paper presents a new parameterization for melt fraction as a function of pressure, temperature, water content and modal cpx, based on knowledge gained from recent advances in the fields of thermodynamic modeling as well as experimental investigations of peridotite melting and hydrous equilibria. The parameterization is computationally efficient and can be modified easily as better experimental data become available. We compare it to other published parameterizations and test it insimple calculations of adiabatic decompression melting (mid-ocean ridge) and hydrous melting (subduction zone).Earth and Planetary Science

Adakitic Dacites Formed by Intracrustal Crystal Fractionation of Water-rich Parent Magmas at Nevado de Longaví Volcano (36·2°S; Andean Southern Volcanic Zone, Central Chile)

The mid-Holocene eruptive products of Nevado de Longaví volcano (36·2°S, Chile) are the only reported occurrence of adakitic volcanic rocks in the Quaternary Andean Southern Volcanic Zone (33-46°S). Dacites of this volcano are chemically distinct from other evolved magmas of the region in that they have high La/Yb (15-20) and Sr/Y (60-90) ratios and systematically lower incompatible element contents. An origin by partial melting of high-pressure crustal sources seems unlikely from isotopic and trace element considerations. Mafic enclaves quenched into one of the dacites, on the other hand, constitute plausible parental magmas. Dacites and mafic enclaves share several characteristics such as mineral chemistry, whole-rock isotope and trace element ratios, highly oxidizing conditions (NNO + 1·5 to >NNO + 2, where NNO is the nickel-nickel oxide buffer), and elevated boron contents. A two-stage mass-balance crystal fractionation model that matches both major and trace elements is proposed to explain magmatic evolution from the least evolved mafic enclave to the dacites. Amphibole is the main ferromagnesian phase in both stages of this model, in agreement with the mineralogy of the magmas. We also describe cumulate-textured xenoliths that correspond very closely to the solid assemblages predicted by the model. We conclude that Nevado de Longaví adakitic dacites are the products of polybaric fractional crystallization from exceptionally water-rich parent magmas. These basaltic magmas are inferred to be related to an exceptionally high, but transient input of slab-derived fluids released from serpentinite bodies hosted in the oceanic Mocha Fracture Zone, which projects beneath Nevado de Longaví. Fractional crystallization that is modally dominated by amphibole, with very minor garnet extraction, is a mechanism for generating adakitic magmas in cold subduction zones where a high flux of slab-derived fluids is presen

Insights into the Mechanism of Intermediate-Depth Earthquakes from Source Properties as Imaged by Back Projection of Multiple Seismic Phases

This study investigates the spatial and temporal distribution of energy release of large, intermediate-depth earthquakes using a modified back projection technique first used to study the 2004 Sumatra-Andaman megathrust event. Multiple seismic phases are included in the back projection analysis, which provides the capability to determine the energy distribution with respect to depth and time. A total of 22 intermediate-depth earthquakes with moment magnitudes greater than or equal to 6.5 are investigated with hypocentral depths between 100 and 300 km. For most of these events, the vertical extent of energy release is either below the resolution of this study $(\leq5 km)$ or slightly above $(\leq15 km)$. This observation agrees with previous studies that find large, intermediate-depth earthquakes have subhorizontal rupture planes. The results also show a significant portion of the events have multiple rupture planes that are well separated in depth. The closeness in time of the ruptures on separate planes and the distance between the planes suggest dynamic triggering where the P waves from the first rupture initiate rupture on the second plane. We propose that a dehydration embrittlement mechanism combined with preferentially hydrated subhorizontal faults can explain the observations of dominant subhorizontal rupture planes and the frequent occurrence of rupture complexity involving multiple subevents.Earth and Planetary Science

Origins of Chemical Diversity of Back-Arc Basin Basalts: A Segment-Scale Study of the Eastern Lau Spreading Center

We report major, trace, and volatile element data on basaltic glasses from the northernmost segment of the Eastern Lau Spreading Center (ELSC1) in the Lau back-arc basin to further test and constrain models of back-arc volcanism. The zero-age samples come from 47 precisely collected stations from an 85 km length spreading center. The chemical data covary similarly to other back-arc systems but with tighter correlations and well-developed spatial systematics. We confirm a correlation between volatile content and apparent extent of melting of the mantle source but also show that the data cannot be reproduced by the model of isobaric addition of water that has been broadly applied to back-arc basins. The new data also confirm that there is no relationship between mantle temperature and the wet melting productivity. Two distinct magmatic provinces can be identified along the ELSC1 axis, a southern province influenced by a wet component with strong affinities to arc volcanism and a northern province influenced by a damp component intermediate between enriched mid-ocean ridge basalts (E-MORB) and arc basalts. High field strength elements and rare earth elements are all mobilized to some extent by the wet component, and the detailed composition of this component is determined. It differs in significant ways from the Mariana component reported by E. Stolper and S. Newman (1994), particularly by having lower abundances of most elements relative to $H_2O$. The differences can be explained if the slab temperature is higher for the Mariana and the source from which the fluid is derived is more enriched. The ELSC1 damp component is best explained by mixing between the wet component and an E-MORB-like component. We propose that mixing between water-rich fluids and low-degree silicate melts occurs at depth in the subduction zone to generate the chemical diversity of the ELSC1 subduction components. These modified sources then rise independently to the surface and melt, and these melts mix with melts of the background mantle from the ridge melting regime to generate the linear data arrays characteristic of back-arc basalts. The major and trace element framework for ELSC1, combined with different slab temperatures and compositions for difference convergent margins, may be able to be applied to other back-arc basins around the globe.Earth and Planetary Science