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

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

On the formation/dissolution of equilibrium droplets

We consider liquid-vapor systems in finite volume $V\subset\R^d$ at parameter values corresponding to phase coexistence and study droplet formation due to a fixed excess $\delta N$ of particles above the ambient gas density. We identify a dimensionless parameter $\Delta\sim(\delta N)^{(d+1)/d}/V$ and a \textrm{universal} value \Deltac=\Deltac(d), and show that a droplet of the dense phase occurs whenever \Delta>\Deltac, while, for \Delta<\Deltac, the excess is entirely absorbed into the gaseous background. When the droplet first forms, it comprises a non-trivial, \textrm{universal} fraction of excess particles. Similar reasoning applies to generic two-phase systems at phase coexistence including solid/gas--where the ``droplet'' is crystalline--and polymorphic systems. A sketch of a rigorous proof for the 2D Ising lattice gas is presented; generalizations are discussed heuristically.Comment: An announcement of a forthcoming rigorous work on the 2D Ising model; to appear in Europhys. Let

Timing of volcanism along the northern East Pacific Rise based on paleointensity experiments on basaltic glasses

Samples from two adjacent and contrasting ridge segments along the East Pacific Rise were measured for their magnetic paleointensity in order to further explore the possibilities of dating very young volcanic samples using secular variations in the Earth's magnetic field. The ridge segment north of the Orozco transform fault (15°22′-16°20′N) is the shallowest and broadest along more than 5000 km of the East Pacific Rise, whereas the adjacent segment to the north (16°16′-18°N) has a "typical" morphology for its intermediate spreading rate. Both ridge segments were densely sampled during the PANR01MV cruise and 36 samples of axial lava flows, consisting mainly of glasses from the rims of the flows and some fragments of lobate basalts, were selected from this collection for paleointensity experiments. The Coe version of the Thellier double-heating procedure (in air) was used. Twenty-seven units provide internally consistent paleointensity estimates leading to precise estimates of the paleofield, which range between 8 μT and 57 μT. Comparisons with reference paleointensity curves compiled from subaerial flows, archeomagnetic data and sedimentary records projected to the sampling site coordinates show that the measured values can be used to constrain the volcanic history of the ridge segments over the past few thousand years. A good agreement was found between apparent "freshness" of the glasses, the geochemistry of the lavas, and their magnetic paleointensity values. The inflated southern segment seems characterized by recent activities as indicated by numerous flows with paleointensities clustering around today's value (39 μT) or around the high values typical of 2000-3000 years ago (~55 μT). We interpret this distribution to indicate the flooding by effusive lava flows of the entire axial plateau some 2000-3000 years ago, followed by a volcanic phase producing smaller volume lava flows confined to the innermost 200 m of the ridge axis. The northern ridge segment is characterized by dispersed paleointensity values consistent with a series of small eruptions of diverse ages. Samples collected at the tips of both ridge segments across the 16°20′N axial discontinuity have the lowest paleointensities and are thus thought to be significantly older, consistent with models advocating reduced magmatism near ridge axis discontinuities. This study demonstrates the strong potential of paleointensity measurements as a tool to help constrain volcanic history at ridge axes

Influence of Hydrodynamic Interactions on the Kinetics of Colloidal Particle's Adsorption

The kinetics of irreversible adsorption of spherical particles onto a flat surface is theoretically studied. Previous models, in which hydrodynamic interactions were disregarded, predicted a power-law behavior $t^{-2/3}$ for the time dependence of the coverage of the surface near saturation. Experiments, however, are in agreement with a power-law behavior of the form $t^{-1/2}$. We outline that, when hydrodynamic interactions are considered, the assymptotic behavior is found to be compatible with the experimental results in a wide region near saturation.Comment: 4 pages, 1 figures, Phys. Rev. Lett. (in press

Geodynamics and Rate of Volcanism on Massive Earth-like Planets

We provide estimates of volcanism versus time for planets with Earth-like composition and masses from 0.25 to 25 times Earth, as a step toward predicting atmospheric mass on extrasolar rocky planets. Volcanism requires melting of the silicate mantle. We use a thermal evolution model, calibrated against Earth, in combination with standard melting models, to explore the dependence of convection-driven decompression mantle melting on planet mass. Here we show that (1) volcanism is likely to proceed on massive planets with plate tectonics over the main-sequence lifetime of the parent star; (2) crustal thickness (and melting rate normalized to planet mass) is weakly dependent on planet mass; (3) stagnant lid planets live fast (they have higher rates of melting than their plate tectonic counterparts early in their thermal evolution) but die young (melting shuts down after a few Gyr); (4) plate tectonics may not operate on high mass planets because of the production of buoyant crust which is difficult to subduct; and (5) melting is necessary but insufficient for efficient volcanic degassing - volatiles partition into the earliest, deepest melts, which may be denser than the residue and sink to the base of the mantle on young, massive planets. Magma must also crystallize at or near the surface, and the pressure of overlying volatiles must be fairly low, if volatiles are to reach the surface. If volcanism is detected in the Tau Ceti system, and tidal forcing can be shown to be weak, this would be evidence for plate tectonics.Comment: Revised version, accepted by Astrophysical Journa

Development of a thermal ionizer as ion catcher

An effective ion catcher is an important part of a radioactive beam facility that is based on in-flight production. The catcher stops fast radioactive products and emits them as singly charged slow ions. Current ion catchers are based on stopping in He and H$_2$ gas. However, with increasing intensity of the secondary beam the amount of ion-electron pairs created eventually prevents the electromagnetic extraction of the radioactive ions from the gas cell. In contrast, such limitations are not present in thermal ionizers used with the ISOL production technique. Therefore, at least for alkaline and alkaline earth elements, a thermal ionizer should then be preferred. An important use of the TRI$\mu$P facility will be for precision measurements using atom traps. Atom trapping is particularly possible for alkaline and alkaline earth isotopes. The facility can produce up to 10$^9$ s$^{-1}$ of various Na isotopes with the in-flight method. Therefore, we have built and tested a thermal ionizer. An overview of the operation, design, construction, and commissioning of the thermal ionizer for TRI$\mu$P will be presented along with first results for $^{20}$Na and $^{21}$Na.Comment: 10 pages, 4 figures, XVth International Conference on Electromagnetic Isotope Separators and Techniques Related to their Applications (EMIS 2007

Properties of Physical Systems: Transient Singularities on Borders and Surface Transitive Zones

Certain alternative properties of physical systems are describable by supports of arguments of response functions (e.g. light cone, borders of media) and expressed by projectors; corresponding equations of restraints lead to dispersion relations, theorems of counting, etc. As supports are measurable, their absolutely strict borders contradict the spirit of quantum theory and their quantum evolution leading to appearance of subtractions or certain needed flattening would be considered. Flattening of projectors introduce transitive zones that can be examined as a specification of adiabatic hypothesis or the Bogoliubov regulatory function in QED. For demonstration of their possibilities the phenomena of refraction and reflection of electromagnetic wave are considered; they show, in particular, the inevitable appearing of double electromagnetic layers on all surfaces that formerly were repeatedly postulated, etc. Quantum dynamics of projectors proves the neediness of subtractions that usually are artificially adding and express transient singularities and zones in squeezed forms.Comment: 12 p

Ewald methods for polarizable surfaces with application to hydroxylation and hydrogen bonding on the (012) and (001) surfaces of alpha-Fe2O3

We present a clear and rigorous derivation of the Ewald-like method for calculation of the electrostatic energy of the systems infinitely periodic in two-dimensions and of finite size in the third dimension (slabs) which is significantly faster than existing methods. Molecular dynamics simulations using the transferable/polarizable model by Rustad et al. were applied to study the surface relaxation of the nonhydroxylated, hydroxylated, and solvated surfaces of alpha-Fe2O3 (hematite). We find that our nonhydroxylated structures and energies are in good agreement with previous LDA calculations on alpha-alumina by Manassidis et al. [Surf. Sci. Lett. 285, L517, 1993]. Using the results of molecular dynamics simulations of solvated interfaces, we define end-member hydroxylated-hydrated states for the surfaces which are used in energy minimization calculations. We find that hydration has a small effect on the surface structure, but that hydroxylation has a significant effect. Our calculations, both for gas-phase and solution-phase adsorption, predict a greater amount of hydroxylation for the (012) surface than for the (001) surface. Our simulations also indicate the presence of four-fold coordinated iron ions on the (001) surface.Comment: 23 pages, REVTeX (LaTeX), 8 figures not included, e-mail to [email protected], paper accepted in Surface Scienc

Mechanical tuning of the evaporation rate of liquid on crossed fibers

We investigate experimentally the drying of a small volume of perfectly wetting liquid on two crossed fibers. We characterize the drying dynamics for the three liquid morphologies that are encountered in this geometry: drop, column and a mixed morphology, in which a drop and a column coexist. For each morphology, we rationalize our findings with theoretical models that capture the drying kinetics. We find that the evaporation rate depends significantly on the liquid morphology and that the drying of liquid column is faster than the evaporation of the drop and the mixed morphology for a given liquid volume. Finally, we illustrate that shearing a network of fibers reduces the angle between them, changes the morphology towards the column state, and so enhances the drying rate of a volatile liquid deposited on it