Melt mixing causes negative correlation of trace element enrichment and CO<inf>2</inf> content prior to an Icelandic eruption

Major elements, trace elements and volatiles were measured in 110 olivine-hosted melt inclusions from the subglacial Skuggafjöll eruption in the Eastern Volcanic Zone of Iceland. Variations in melt inclusion trace element concentrations can be accounted for by incomplete mixing of diverse mantle parental melts accompanied by variable extents of fractional crystallisation. Binary mixing between an incompatible trace element-enriched and depleted melts provides a good fit to observed variations in trace element ratios such as Ce/Y. Surprisingly, the CO2 contents of melt inclusions correlate negatively with their degree of trace element enrichment. Depleted, low-Ce/Y inclusions with ∼1200 ppm CO2 have high CO2/Nb contents (∼400), suggesting that melts experienced little or no CO2 exsolution before inclusion entrapment. Enriched, high-Ce/Y inclusions contain ∼300 ppm CO2, have low CO2/Nb (contents 50–100) and melts are likely to have exsolved much of their original CO2 contents prior to inclusion entrapment. The negative correlation between CO2 content and trace element enrichment may arise either from the more efficient exsolution of CO2 from enriched melts, or from the intrusion of CO2-supersaturated depleted melts into enriched melts that had already exsolved much of their original CO2 contents. Some inclusions have lower CO2 contents than predicted from binary mixing models, which suggests that at least some CO2 exsolution occurred concurrently with mixing. Enriched inclusions record entrapment pressures of ∼0.5 kbar. These pressures probably correspond to the depth of mixing. Higher pressures recorded in depleted inclusions may have resulted from the development of CO2 supersaturation during ascent from storage at ≥1.5 kbar. The presence of CO2 supersaturation in melt inclusions has the potential to constrain timescales of melt inclusion entrapment.This work was supported by a Natural Environment Research Council studentship to D.A.N. (NE/1528277/1) and a Natural Environment Research Council Ion Microprobe Facility award (IMF461/0512).This is the final published version of the article, which can also be found online here: http://www.sciencedirect.com/science/article/pii/S0012821X14003604

Magma mixing and high fountaining during the 1959 K īlauea Iki eruption, Hawai'i

The 1959 Kīlauea Iki eruption provides a unique opportunity to investigate the process of shallow magma mixing, its impact on the magmatic volatile budget and its role in triggering and driving episodes of Hawaiian fountaining. Melt inclusions hosted by olivine record a continuous decrease in H2O concentration through the 17 episodes of the eruption, while CO2 concentrations correlate with the degree of post-entrapment crystallization of olivine on the inclusion walls. Geochemical data, when combined with the magma budget and with contemporaneous eruption observations, show complex mixing between episodes involving hot, geochemically heterogeneous melts from depth, likely carrying exsolved vapour, and melts which had erupted at the surface, degassed and drained-back into the vent. The drained-back melts acted as a coolant, inducing rapid cooling of the more primitive melts and their olivines at shallow depths and inducing crystallization and vesiculation and triggering renewed fountaining. A consequence of the mixing is that the melts became vapor-undersaturated, so equilibration pressures cannot be inferred from them using saturation models. After the melt inclusions were trapped, continued growth of vapor bubbles, caused by enhanced post-entrapment crystallization, sequestered a large fraction of CO2 from the melt within the inclusions. This study, while cautioning against accepting melt inclusion CO2 concentrations “as measured” in mixed magmas, also illustrates that careful analysis and interpretation of post-entrapment modifications can turn this apparent challenge into a way to yield novel useful insights into the geochemical controls on eruption intensity.IS was supported by a NERC-funded studentship and a USGS Jack Kleinman Grant for Volcano Research. ME acknowledges NERC ion probe grant IMF376/0509. BH’s participation was funded by NSF EAR-1145159. We acknowledge the NERC Edinburgh Ion Microprobe facility, where we undertook the SIMS analyses.This version is the author accepted manuscript and will be under embargo until the 6th of June 2015. The final version has been published by Elsevier in Earth and Planetary Science Letters here: http://www.sciencedirect.com/science/article/pii/S0012821X14003264

Smoothed Analysis of the Minimum-Mean Cycle Canceling Algorithm and the Network Simplex Algorithm

The minimum-cost flow (MCF) problem is a fundamental optimization problem with many applications and seems to be well understood. Over the last half century many algorithms have been developed to solve the MCF problem and these algorithms have varying worst-case bounds on their running time. However, these worst-case bounds are not always a good indication of the algorithms' performance in practice. The Network Simplex (NS) algorithm needs an exponential number of iterations for some instances, but it is considered the best algorithm in practice and performs best in experimental studies. On the other hand, the Minimum-Mean Cycle Canceling (MMCC) algorithm is strongly polynomial, but performs badly in experimental studies. To explain these differences in performance in practice we apply the framework of smoothed analysis. We show an upper bound of $O(mn^2\log(n)\log(\phi))$ for the number of iterations of the MMCC algorithm. Here $n$ is the number of nodes, $m$ is the number of edges, and $\phi$ is a parameter limiting the degree to which the edge costs are perturbed. We also show a lower bound of $\Omega(m\log(\phi))$ for the number of iterations of the MMCC algorithm, which can be strengthened to $\Omega(mn)$ when $\phi=\Theta(n^2)$. For the number of iterations of the NS algorithm we show a smoothed lower bound of $\Omega(m \cdot \min \{ n, \phi \} \cdot \phi)$.Comment: Extended abstract to appear in the proceedings of COCOON 201

Stronger or longer: Discriminating between Hawaiian and Strombolian eruption styles

The weakest explosive volcanic eruptions globally, Strombolian explosions and Hawaiian fountaining, are also the most common. Yet, despite over a hundred years of observations, no classifications have offered a convincing, quantitative way of demarcating these two styles. New observations show that the two styles are distinct in their eruptive time scale, with the duration of Hawaiian fountaining exceeding Strombolian explosions by ∼300–10,000 s. This reflects the underlying process of whether shallow-exsolved gas remains trapped in the erupting magma or is decoupled from it. We propose here a classification scheme based on the duration of events (brief explosions versus prolonged fountains) with a cutoff at 300 s that separates transient Strombolian explosions from sustained Hawaiian fountains.The authors wish to acknowledge grants from NSF (EAR-0409303, 0810332, 1145159, 1427357) and ARRA (113153 via the Hawaiian Volcano Observatory), which funded this research.This is the accepted manuscript. The final version is available at http://dx.doi.org/10.1130/G37423.

Analytic descriptions for transitional nuclei near the critical point

Exact solutions of the Bohr Hamiltonian with a five-dimensional square well potential, in isolation or coupled to a fermion by the five-dimensional spin-orbit interaction, are considered as examples of a new class of dynamical symmetry or Bose-Fermi dynamical symmetry. The solutions provide baselines for experimental studies of even-even [E(5)] and odd-mass [E(5|4)] nuclei near the critical point of the spherical to deformed gamma-unstable phase transition.Comment: LaTeX (elsart), 53 pages; typographical correction to (3.15

Indirect exchange in GaMnAs bilayers via spin-polarized inhomogeneous hole gas: Monte Carlo simulation

The magnetic order resulting from an indirect exchange between magnetic moments provided by spin-polarized hole gas in the metallic phase of a GaMnAs double layer structure is studied via Monte Carlo simulation. The coupling mechanism involves a perturbative calculation in second order of the interaction between the magnetic moments and carriers (holes). We take into account a possible polarization of the hole gas due to the existence of an average magnetization in the magnetic layers, establishing, in this way, a self-consistency between the magnetic order and the electronic structure. That interaction leads to an internal ferromagnetic order inside each layer, and a parallel arrangement between their magnetizations, even in the case of thin layers. This fact is analyzed in terms of the inter- and intra-layer interactions.Comment: 17 pages and 14 figure

Slepian functions and their use in signal estimation and spectral analysis

It is a well-known fact that mathematical functions that are timelimited (or spacelimited) cannot be simultaneously bandlimited (in frequency). Yet the finite precision of measurement and computation unavoidably bandlimits our observation and modeling scientific data, and we often only have access to, or are only interested in, a study area that is temporally or spatially bounded. In the geosciences we may be interested in spectrally modeling a time series defined only on a certain interval, or we may want to characterize a specific geographical area observed using an effectively bandlimited measurement device. It is clear that analyzing and representing scientific data of this kind will be facilitated if a basis of functions can be found that are "spatiospectrally" concentrated, i.e. "localized" in both domains at the same time. Here, we give a theoretical overview of one particular approach to this "concentration" problem, as originally proposed for time series by Slepian and coworkers, in the 1960s. We show how this framework leads to practical algorithms and statistically performant methods for the analysis of signals and their power spectra in one and two dimensions, and on the surface of a sphere.Comment: Submitted to the Handbook of Geomathematics, edited by Willi Freeden, Zuhair M. Nashed and Thomas Sonar, and to be published by Springer Verla

Search for direct production of charginos and neutralinos in events with three leptons and missing transverse momentum in √s = 7 TeV pp collisions with the ATLAS detector

A search for the direct production of charginos and neutralinos in final states with three electrons or muons and missing transverse momentum is presented. The analysis is based on 4.7 fb−1 of proton–proton collision data delivered by the Large Hadron Collider and recorded with the ATLAS detector. Observations are consistent with Standard Model expectations in three signal regions that are either depleted or enriched in Z-boson decays. Upper limits at 95% confidence level are set in R-parity conserving phenomenological minimal supersymmetric models and in simplified models, significantly extending previous results

Measurement of the inclusive and dijet cross-sections of b-jets in pp collisions at sqrt(s) = 7 TeV with the ATLAS detector

The inclusive and dijet production cross-sections have been measured for jets containing b-hadrons (b-jets) in proton-proton collisions at a centre-of-mass energy of sqrt(s) = 7 TeV, using the ATLAS detector at the LHC. The measurements use data corresponding to an integrated luminosity of 34 pb^-1. The b-jets are identified using either a lifetime-based method, where secondary decay vertices of b-hadrons in jets are reconstructed using information from the tracking detectors, or a muon-based method where the presence of a muon is used to identify semileptonic decays of b-hadrons inside jets. The inclusive b-jet cross-section is measured as a function of transverse momentum in the range 20 < pT < 400 GeV and rapidity in the range |y| < 2.1. The bbbar-dijet cross-section is measured as a function of the dijet invariant mass in the range 110 < m_jj < 760 GeV, the azimuthal angle difference between the two jets and the angular variable chi in two dijet mass regions. The results are compared with next-to-leading-order QCD predictions. Good agreement is observed between the measured cross-sections and the predictions obtained using POWHEG + Pythia. MC@NLO + Herwig shows good agreement with the measured bbbar-dijet cross-section. However, it does not reproduce the measured inclusive cross-section well, particularly for central b-jets with large transverse momenta.Comment: 10 pages plus author list (21 pages total), 8 figures, 1 table, final version published in European Physical Journal

Inclusive search for same-sign dilepton signatures in pp collisions at root s=7 TeV with the ATLAS detector

An inclusive search is presented for new physics in events with two isolated leptons (e or mu) having the same electric charge. The data are selected from events collected from p p collisions at root s = 7 TeV by the ATLAS detector and correspond to an integrated luminosity of 34 pb(-1). The spectra in dilepton invariant mass, missing transverse momentum and jet multiplicity are presented and compared to Standard Model predictions. In this event sample, no evidence is found for contributions beyond those of the Standard Model. Limits are set on the cross-section in a fiducial region for new sources of same-sign high-mass dilepton events in the ee, e mu and mu mu channels. Four models predicting same-sign dilepton signals are constrained: two descriptions of Majorana neutrinos, a cascade topology similar to supersymmetry or universal extra dimensions, and fourth generation d-type quarks. Assuming a new physics scale of 1 TeV, Majorana neutrinos produced by an effective operator V with masses below 460 GeV are excluded at 95% confidence level. A lower limit of 290 GeV is set at 95% confidence level on the mass of fourth generation d-type quarks