Evidence for melt leakage from the Hawaiian plume above the mantle transition zone

Dehydration reactions at the top of the mantle transition zone (MTZ) can stabilize partial melt in a seismic low-velocity layer (LVL), but the seismic effects of temperature, melt and volatile content are difficult to distinguish. We invert P-to-S receiver function phases converted at the top and bottom of a LVL above the MTZ beneath Hawaii. To separate the thermal and melting related seismic anomalies, we carry out over 10 million rock physics inversions. These inversions account for variations arising from the Clapeyron slope of phase transition, bulk solid composition, dihedral angle, and mantle potential temperature. We use two independent seismic constraints to evaluate the temperature and shear wave speed within the LVL. The thermal anomalies reveal the presence of a hot and seismically slow plume stem surrounded by a “halo” of cold and fast mantle material. In contrast to this temperature distribution, the plume stem contains less than 0.5 vol% melt, while the surrounding LVL—within the coverage area—contains up to 1.7 vol% melt, indicating possible lateral transport of the melt. When compared to the melting temperatures of mantle rocks, the temperature within the LVL, calculated from seismic observations of MTZ thickness, suggests that the observed small degrees of melting are sustained by the presence of volatiles such as CO2 and H2O. We estimate the Hawaiian plume loses up to 1.9 Mt/yr H2O and 10.7 Mt/yr CO2 to the LVL, providing a crucial missing flux for global volatile cycles

Evidence of Volatile-Induced Melting in the Northeast Asian Upper Mantle

A seismic low velocity layer (LVL) above the mantle transition zone (MTZ), often thought to be caused by volatile-induced melting, can significantly modulate planetary volatile cycles. In this work, we show that an LVL observed beneath northeast Asia is characterized by small, 0.8 (Formula presented.) 0.5 vol%, average degrees of partial melting. Seismically derived P-T conditions of the LVL indicate that slab-derived (Formula presented.), possibly combined with small amounts of (Formula presented.) O, is necessary to induce melting. Modeling the reactive infiltration instability of the melt in a stationary mantle above a stalled slab, we demonstrate that the volatile-rich melt slowly rises above the stalled slab in the MTZ, with percolation velocities of 200–500 (Formula presented.) m/yr. The melt remains stable within the LVL for this geologically significant period of time, potentially transferring up to 52 Mt/yr of (Formula presented.) from the subducting slab to the mantle for an LVL similar in areal extent ((Formula presented.)) to the northeast Asian LVL. Reaction between the melt channels and the LVL mantle precipitates up to 200 ppmw solid C in localized zones. Using the inferred small melt volume fraction to model trace element abundances and isotopic signatures, we show that interaction between this melt and the surrounding mantle can over the long-term produce rocks bearing a HIMU like geochemical signature

Pore network analysis of Brae Formation sandstone, North Sea

In this work, we apply digital rock physics (DRP) to characterize the pore networks of the Brae Formation sandstones from two different wells in the Miller field area (North Sea, UK). Using X-ray micro-CT scans, we calculate the porosity and permeability and generate pore network models to assess pore shape characteristics. The porous samples are marked by macroporosities ranging from 4.9% to 15.2% with the effective porosities varying from 0 to 14.8%. The samples also contained some microporosity hosted in secondary and accessory mineral phases, varying between 2.6% and 10.7%. Pore network model results for total porosity indicate that the samples have median pore and throat radii ranging from 5.5 μm to 16.8 μm and 6.4 μm–12.9 μm, respectively. The throat length of all samples has a median value ranging between 36.3 μm and 82.4 μm. The ratio between effective porosity and total porosity (φ∗) varies with total porosity (φ) following the exponential relation φ∗ = 0.98 − e− (φ− 0.032)/0.028. Pore network connectivity is established at a porosity of 3% and full communication is achieved at porosities exceeding 10%. Permeability was found to vary with total porosity with an exponent of 3.67. Based on these observations and the results from our models, the connectivity of the pore network has important implications for predicting reservoir performance during large scale subsurface projects such as hydrocarbon production and CO2 storage

Data assimilation for plume models

International audienceWe use a four-dimensional variational data assimilation (4D-VAR) algorithm to observe the growth of 2-D plumes from a point heat source. In order to test the predictability of the 4D-VAR technique for 2-D plumes, we perturb the initial conditions and compare the resulting predictions to the predictions given by a direct numerical simulation (DNS) without any 4D-VAR correction. We have studied plumes in fluids with Rayleigh numbers between 106 and 107 and Prandtl numbers between 0.7 and 70, and we find the quality of the prediction to have a definite dependence on both the Rayleigh and Prandtl numbers. As the Rayleigh number is increased, so is the quality of the prediction, due to an increase of the inertial effects in the adjoint equations for momentum and energy. The horizon predictability time, or how far into the future the 4D-VAR method can predict, decreases as Rayleigh number increases. The quality of the prediction is decreased as Prandtl number increases, however. Quality also decreases with increased prediction time

Textures in experimentally deformed olivine aggregates: the effects of added water and melt

Abstract. The texture development in experimentally sheared aggregates of olivine was monitored as a function of increased water content and added melt. In dry samples, an alignment of {010} with the shear plane and <100> and <001> with the shear direction, respectively, was observed, consistent with intracrystalline glide on the (010) [100] an

Mass-Radius Relationships for Solid Exoplanets

We use new interior models of cold planets to investigate the mass-radius relationships of solid exoplanets, considering planets made primarily of iron, silicates, water, and carbon compounds. We find that the mass-radius relationships for cold terrestrial-mass planets of all compositions we considered follow a generic functional form that is not a simple power law: $\log_{10} R_s = k_1 + 1/3 \log_{10}(M_s) - k_2 M_s^{k_3}$ for up to $M_p \approx 20 M_{\oplus}$, where $M_s$ and $R_s$ are scaled mass and radius values. This functional form arises because the common building blocks of solid planets all have equations of state that are well approximated by a modified polytrope of the form $\rho = \rho_0 + c P^n$. We find that highly detailed planet interior models, including temperature structure and phase changes, are not necessary to derive solid exoplanet bulk composition from mass and radius measurements. For solid exoplanets with no substantial atmosphere we have also found that: with 5% fractional uncertainty in planet mass and radius it is possible to distinguish among planets composed predominantly of iron or silicates or water ice but not more detailed compositions; with $\sim$~5% uncertainty water ice planets with $\gtrsim 25$ water by mass may be identified; the minimum plausible planet size for a given mass is that of a pure iron planet; and carbon planet mass-radius relationships overlap with those of silicate and water planets due to similar zero-pressure densities and equations of state. We propose a definition of "super Earths'' based on the clear distinction in radii between planets with significant gas envelopes and those without.Comment: ApJ, in press, 33 pages including 16 figure

Processes controlling lithium isotopic distribution in contact aureoles: A case study of the Florence County pegmatites, Wisconsin

Li isotopes may be useful tracers of fluid flow in a number of geological environments and case studies of contact aureoles have highlighted the very large Li isotopic fractionation that can be generated in these settings. However, the amount of isotopic fractionation and the distance that Li travels into the country rocks vary greatly from place to place. Seeking to identify the parameters that govern Li distribution in contact aureoles, we apply a combination of Li isotope analyses, 1-D diffusion and 2-D advection-diffusion modeling to two country rock profiles adjacent to Li-rich pegmatite dikes from the Florence County pegmatite field, Wisconsin. Although less than ∼3 m thick, the pegmatite sheets have a large impact on the Li budget of the country rocks (amphibolites and schists); Li is enriched in adjacent country rocks by up to a factor of 20 over more distant amphibolites and schists. Li from the pegmatite has traveled more than 50 m into the country rocks, and Li isotopes are systematically fractionated with distance from the contacts (with δ7Li varying from +6 at the contact to-7 at 30 m from the contact in one case). These observations are consistent with diffusive fractionation of Li through an advecting grain-boundary fluid. Both one-dimensional diffusion and two-dimensional advection-diffusion models fail to reproduce the exact Li distribution in the profiles, suggesting that fluid advection, coupled with heterogeneous permeability, plays an important role in determining the final Li distribution within the contact aureoles

Digital rock physics analysis of core integrity using deep neural networks and computer vision: MR23C-0124

Core imaging and classification is an important step for generating a digital database for subsurface geology. The British Geological Survey collection contains cores from over 15,000 onshore and 8,000 offshore boreholes. Many cores are photographed at high-resolution creating an archive of over 100,000 core tray images containing between 1 m and 3 m of core per image. A crucial challenge in storage and classification of the digital data associated with these cores is the degree of damage or degradation which determines their suitability for further analysis by the scientific community. The large volume of core image data precludes manual phase segmentation and core quality determination from individuals. In this work, we present a new method using both deep learning and computer vision to automate the process. To test the feasibility of our technique, we use a small subset of 62 core tray images, captured with 3 light spectra (Red, Green, Blue). We use pre-trained neural networks to segment the image, which is followed by traditional computer vision techniques for edge detection. We also automate the process of calculating the number of fragments and area of each fragment present in each individual core image. Finally, we present an index for core integrity based on the output of these measurements. The work-flow demonstrates that deep neural networks and computer vision can be leveraged to quantify and non-intrusively assess geophysical properties at a large scale, using only a subset of the data, with open-source packages. This core quality index will allow users to quickly and consistently assess core condition, and in particularly degradation due to: transport, storage and previous sampling. By automating this process it is possible to quickly assess tens to hundreds of metres of core to identify areas suitable for sampling. It also provides semi quantitative information on how representative individual core samples are of bulk rock properties. This will improve integration between core analysis and other datasets, for example wireline logs

Electrical conductivity during incipient melting in the oceanic low-velocity zone

International audienceThe low-viscosity layer in the upper mantle, the asthenosphere, is a requirement for plate tectonics1. The seismic low velocities and the high electrical conductivities of the asthenosphere are attributed either to subsolidus, water-related defects in olivine minerals2, 3, 4 or to a few volume per cent of partial melt5, 6, 7, 8, but these two interpretations have two shortcomings. First, the amount of water stored in olivine is not expected to be higher than 50 parts per million owing to partitioning with other mantle phases9 (including pargasite amphibole at moderate temperatures10) and partial melting at high temperatures9. Second, elevated melt volume fractions are impeded by the temperatures prevailing in the asthenosphere, which are too low, and by the melt mobility, which is high and can lead to gravitational segregation11, 12. Here we determine the electrical conductivity of carbon-dioxide-rich and water-rich melts, typically produced at the onset of mantle melting. Electrical conductivity increases modestly with moderate amounts of water and carbon dioxide, but it increases drastically once the carbon dioxide content exceeds six weight per cent in the melt. Incipient melts, long-expected to prevail in the asthenosphere10, 13, 14, 15, can therefore produce high electrical conductivities there. Taking into account variable degrees of depletion of the mantle in water and carbon dioxide, and their effect on the petrology of incipient melting, we calculated conductivity profiles across the asthenosphere for various tectonic plate ages. Several electrical discontinuities are predicted and match geophysical observations in a consistent petrological and geochemical framework. In moderately aged plates (more than five million years old), incipient melts probably trigger both the seismic low velocities and the high electrical conductivities in the upper part of the asthenosphere, whereas in young plates4, where seamount volcanism occurs6, a higher degree of melting is expected