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

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

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

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

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

Petrological systematics of mid-ocean ridge basalts: Constraints on melt generation beneath ocean ridges

Mid-ocean ridge basalts (MORB) are a consequence of pressure-release melting beneath ocean ridges, and contain much information concerning melt formation, melt migration and heterogeneity within the upper mantle. MORB major element chemical systematics can be divided into global and local aspects, once they have been corrected for low pressure fractionation and interlaboratory biases. Regional average compositions for ridges unaffected by hot spots (“normal” ridges) can be used to define the global correlations among normalized Na2O, FeO, TiO2 and SiO2 contents, CaO/Al2O3 ratios, axial depth and crustal thickness. Back-arc basins show similar correlations, but are offset to lower FeO and TiO2 contents. Some hot spots, such as the Azores and Galapagos, disrupt the systematics of nearby ridges and have the opposite relationships between FeO, Na2O and depth over distances of 1000 km.Local variations in basalt chemistry from slow- and fast-spreading ridges are distinct from one another. On slow-spreading ridges, correlations among the elements cross the global vector of variability at a high angle. On the fast-spreading East Pacific Rise (EPR), correlations among the elements are distinct from both global and slow-spreading compositional vectors, and involve two components of variation. Spreading rate does not control the global correlations, but influences the standard deviations of axial depth, crustal thickness, and MgO contents of basalts.Global correlations are not found in very incompatible trace elements, even for samples far from hot spots. Moderately compatible trace elements for normal ridges, however, correlate with the major elements. Trace element systematics are significantly different for the EPR and the mid-Atlantic Ridge (MAR). Normal portions of the MAR are very depleted in REE, with little variability; hot spots cause large long wavelength variations in REE abundances. Normal EPR basalts are significantly more enriched than MAR basalts from normal ridges, and still more enriched basalts can erupt sporadically along the entire length of the EPR. This leads to very different histograms of distribution for the data sets as a whole, and a very different distribution of chemistry along strike for the two ridges. Despite these differences, the mean Ce/Sm ratios from the two ridges are identical.Existing methods for calculating the major element compositions of mantle melts [Klein and Langmuir, 1987; McKenzie and Bickle, 1988; Niu and Batiza, 1991] are critically examined. New quantitative methods for mantle melting and high pressure fractionation are developed to evaluate the chemical consequences of melting and fractionation processes and mantle heterogeneity. The new methods rely on new equations for partition coefficients for the major elements between mantle minerals and melts. The melting calculations can be used to investigate the chemical compositions produced by small extents of melting or high pressures of melting that cannot yet be determined experimentally. Application of the new models to the observations described above leads to two major conclusions: (1) The global correlations for normal ridges are caused by variations in mantle temperature, as suggested by Klein and Langmuir [1987] and not by mantle heterogeneity. (2) Local variations are caused by melting processes, but are not yet quantitatively accounted for. On slower spreading ridges, local variations are controlled by the melting regime in the mantle. On the EPR, local variations are predominantly controlled by ubiquitous, small scale heterogeneites. Volatile content may be an important and as yet undetermined factor in affecting the observed variations in major elements.We propose a hypothesis, similar to one proposed by Allegre et al [1984] for isotopic data, to explain the differences between the Atlantic and Pacific local trends, and the trace element systematics of the two ocean basins, as consequences of spreading rate and a different distribution of enriched components from hot spots in the two ocean basins. In the Atlantic, the hot spot influence is in discrete areas, and produces clear depth and chemical anomalies. Ridge segments far from hot spots do not contain enriched basalts. Melting processes associated with slow-spreading ridges vary substantially over short distances along strike and lead to the local trends discussed above, irrespective of hot spot influence. In the Pacific, enriched components appear to have been more thoroughly mixed into the mantle, leading to ubiquitous small scale heterogeneities. Melting processes do not vary appreciably along strike, so local chemical variations are dominated by the relative contribution of enriched component on short time and length scales. Thus the extent of mixing and distribution of enriched components influences strongly the contrasting local major element trends. Despite the difference in the distribution of enriched components, the mean compositions of each data set are equivalent. This suggests that the hot spot influence is similar in the two ocean basins, but its distribution in the upper mantle is different. These contrasting relationships between hot spots and ridges may result from differences in both spreading rate and tectonic history. Unrecognized hot spots may play an important role in diverse aspects of EPR volcanism, and in the chemical systematics of the erupted basalts.The observations and successful models have consequences for melt formation and segregation. (1) The melting process must be closer to fractional melting than equilibrium melting. This result is in accord with inferences from abyssal peridotites [Johnson et al., 1990]. (2) Small melt fractions generated over a range of pressures must be extracted rapidly and efficiently from high pressures within the mantle without experiencing low pressure equilibration during ascent. This requires movement in large channels, and possibly more efficient extraction mechanisms than nonnally envisaged in porous flow models with small residual porosity. (3) Diverse melts from the melting regime produce variations in basalts that are observable at the surface. (4) Basalt data can be used to constrain the melting process (e.g. active vs. passive upwelling) and its relationship to segmentation. The data cannot be used to constrain the shape of the melting regime, however, for many shapes lead to similar chemical results. (5) Highly incompatible elements and U-series disequilibria results appear not yet to be explained by melting models, and may require additional processes not yet clearly envisaged

Space Charge Limited 2-d Electron Flow between Two Flat Electrodes in a Strong Magnetic Field

An approximate analytic solution is constructed for the 2-d space charge limited emission by a cathode surrounded by non emitting conducting ledges of width Lambda. An essentially exact solution (via conformal mapping) of the electrostatic problem in vacuum is matched to the solution of a linearized problem in the space charge region whose boundaries are sharp due to the presence of a strong magnetic field. The current density growth in a narrow interval near the edges of the cathode depends strongly on Lambda. We obtain an empirical formula for the total current as a function of Lambda which extends to more general cathode geometries.Comment: 4 pages, LaTex, e-mail addresses: [email protected], [email protected]

Overcoming the Child-Langmuir law via the magnetic mirror effect

The maximum current in a vacuum tube prescribed by the classical Child-Langmuir law can be overcome, when the space-charge effect of the induced potential is mitigated by the mirror effect in a spatially varying magnetic field. The current could exceed the Child-Langmuir value by as much as a few factors. The regime of practical interest is examined