Melting of the subcontinental lithospheric mantle by the Emeishan mantle plume; evidence from the basal alkaline basalts in Dongchuan, Yunnan, Southwestern China

The Emeishan continental flood basalt (ECFB) sequence in Dongchuan, SW China comprises a basal tephrite unit overlain by an upper tholeiitic basalt unit. The upper basalts have high TiO₂ contents (3.2–5.2 wt.%), relatively high rare-earth element (REE) concentrations (40 to 60 ppm La, 12.5 to 16.5 ppm Sm, and 3 to 4 ppm Yb), moderate Zr/Nb and Nb/La ratios (9.3–10.2 and 0.6–0.9, respectively) and relatively high _{Nd (t)} values, ranging from − 0.94 to 2.3, and are comparable to the high-Ti ECFB elsewhere. The tephrites have relatively high P₂O₅ (1.3–2.0 wt.%), low REE concentrations (e.g., 17 to 23 ppm La, 4 to 5.3 ppm Sm, and 2 to 3 ppm Yb), high Nb/La (2.0–3.9) ratios, low Zr/Nb ratios (2.3–4.2), and extremely low _{Nd (t)} values (mostly ranging from − 10.6 to − 11.1). The distinct compositional differences between the tephrites and the overlying tholeiitic basalts cannot be explained by either fractional crystallization or crustal contamination of a common parental magma. The tholeiitic basalts formed by partial melting of the Emeishan plume head at a depth where garnet was stable, perhaps > 80 km. We propose that the tephrites were derived from magmas formed when the base of the previously metasomatized, volatile-mineral bearing subcontinental lithospheric mantle was heated by the upwelling mantle plume.     

A late Neoproterozoic paleomagnetic pole for the Congo craton: Tectonic setting, paleomagnetism and geochronology of the Nola dike swarm (Central African Republic)

New structural, geochronological and paleomagnetic data were obtained on dolerite dikes of the Nola region (Central African Republic) at the northern border of the Congo craton. In this region, metavolcanic successions were thrust southward onto the craton during the Panafrican orogenic events. Our structural data reveal at least two structural klippes south of the present-day limits of the Panafrican nappe suggesting that it has once covered the whole Nola region, promoting the pervasive hydrothermal greenschist metamorphism observed in the underlying cratonic basement and also in the intrusive dolerite dikes. Paleomagnetic measurements revealed a stable dual-polarity low-inclination magnetization component in nine dikes (47 samples), carried by pyrrhotite and magnetite. This component corresponds to a paleopole at 304.8°E and 61.8°S (dp = 5.4, dm = 10.7) graded at Q = 6. Both metamorphism and magnetic resetting were dated by the ⁴⁰Ar/³⁹Ar method on amphibole grains separated from the dikes at 571 ± 6 Ma. The Nola pole is the first well-dated paleomagnetic pole for the Congo craton between 580 and 550 Ma. It marks a sudden change in direction of the Congo craton apparent polar wander path at the waning stages of the Panafrican orogenic events.     

Sulfur isotope measurement of sulfate and sulfide by high-resolution MC-ICP-MS

We have developed a technique for the accurate and precise determination of ³⁴S/³²S isotope ratios (δ³⁴S) in sulfur-bearing minerals using solution and laser ablation multiple-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). We have examined and determined rigorous corrections for analytical difficulties such as instrumental mass bias, unresolved isobaric interferences, blanks, and laser ablation- and matrix-induced isotopic fractionation. Use of high resolution sector-field mass spectrometry removes major isobaric interferences from O₂⁺. Standard-sample bracketing is used to correct for the instrumental mass bias of unknown samples. Background on sulfur masses arising from memory effects and residual oxygen-tailing are typically minor (< 0.2‰, within analytical error), and are mathematically removed by on-peak zero subtraction and by bracketing of samples with standards determined at the same signal intensity (within 20%). Matrix effects are significant (up to 0.7‰) for matrix compositions relevant to many natural sulfur-bearing minerals. For solution analysis, sulfur isotope compositions are best determined using purified (matrix-clean) sulfur standards and sample solutions using the chemical purification protocol we present. For in situ analysis, where the complex matrix cannot be removed prior to analysis, appropriately matrix-matching standards and samples removes matrix artifacts and yields sulfur isotope ratios consistent with conventional techniques using matrix-clean analytes. Our method enables solid samples to be calibrated against aqueous standards; a consideration that is important when certified, isotopically-homogeneous and appropriately matrix-matched solid standards do not exist. Further, bulk and in situ analyses can be performed interchangeably in a single analytical session because the instrumental setup is identical for both. We validated the robustness of our analytical method through multiple isotope analyses of a range of reference materials and have compared these with isotope ratios determined using independent techniques. Long-term reproducibility of S isotope compositions is typically 0.20‰ and 0.45‰ (2σ) for solution and laser analysis, respectively. Our method affords the opportunity to make accurate and relatively precise S isotope measurement for a wide range of sulfur-bearing materials, and is particularly appropriate for geologic samples with complex matrix and for which high-resolution in situ analysis is critical.     

Numerical modeling assessment of three-gate structures for capturing contaminated groundwater

This modeling study evaluated the capability of alternative funnel-and-gate structures with three gates for capturing contaminated groundwater in a hypothetical unconfined aquifer. Simulated interceptor structures were linear and 45 m wide, consisting of three gates and two funnels (walls). One gate occupied the center and two gates occupied the ends of the interceptor structures. The structures, positioned perpendicular to regional groundwater flow, traversed the entire thickness of the aquifer. A total of four structures were evaluated (numbers designate widths of end, center, and end gates, respectively, in meters): 3-3-3, 2-5-2, 1-7-1, and 4-1-4. Particle tracking and zonal water budgets identified shapes of capture zones and discharge patterns for each interceptor structure. A mass transport model, accounting for advection and hydrodynamic dispersion, tested the capability of each structure for capturing a contaminant plume. Results suggest that: time-dependent capture zones underestimate the amount of time to capture a contaminant plume, wide center gates facilitate plume capture, and wide end gates facilitate lateral containment of contaminants. Of the structures simulated, the 2-5-2 configuration was relatively efficient at processing and containing the simulated contaminant plume.     

Late Palaeozoic hydrocarbon migration through the Clair field, West of Shetland, UK Atlantic margin

Geochemical analysis of bitumen- and hydrocarbon-bearing fluid inclusions from the Devonian-Carboniferous Clair field indicates that the reservoirs contain a mixture of oils from different marine and lacustrine sources. Reconstruction of the Clair field oil-charge history using fluid inclusion petrography show that oil-charging occurred at times of K-feldspar, quartz and calcite cementation. Temperature–composition–time data yielded from the integration of fluid inclusion microthermometry with high-resolution Ar–Ar dating, date hydrocarbon-bearing K-feldspar overgrowths at 247 ± 3.3 Ma. These data show that in order for oil to be trapped within primary fluid inclusions in K-feldspar overgrowths, hydrocarbon migration throughout the UK Atlantic margin must have been taking place during the Late Palaeozoic and as such, current industry oil-play models based solely on oil charging from Jurassic-Cretaceous marine sources are clearly incomplete and need revision. Apatite fission track analysis and vitrinite reflectance data were used to reconstruct thermal burial histories and assess potential oil generation from Middle Devonian lacustrine source rocks. Thermal history data from wells along The Rona Ridge adjacent to the Clair field show that the Palaeozoic section was heated to greater than 100 °C at some time between 270 and 230 Ma, confirming that Devonian source rocks were mature and expelling oil during the Late Palaeozoic at the time that authigenic K-feldspar overgrowths were growing in the Clair field.     

An experimental study of the grain-scale processes of peridotite melting: implications for major and trace element distribution during equilibrium and disequilibrium melting

The grain-scale processes of peridotite melting were examined at 1,340°C and 1.5 GPa using reaction couples formed by juxtaposing pre-synthesized clinopyroxenite against pre-synthesized orthopyroxenite or harzburgite in graphite and platinum-lined molybdenum capsules. Reaction between the clinopyroxene and orthopyroxene-rich aggregates produces a melt-enriched, orthopyroxene-free, olivine + clinopyroxene reactive boundary layer. Major and trace element abundance in clinopyroxene vary systematically across the reactive boundary layer with compositional trends similar to the published clinopyroxene core-to-rim compositional variations in the bulk lherzolite partial melting studies conducted at similar P–T conditions. The growth of the reactive boundary layer takes place at the expense of the orthopyroxenite or harzburgite and is consistent with grain-scale processes that involve dissolution, precipitation, reprecipitation, and diffusive exchange between the interstitial melt and surrounding crystals. An important consequence of dissolution–reprecipitation during crystal-melt interaction is the dramatic decrease in diffusive reequilibration time between coexisting minerals and melt. This effect is especially important for high charged, slow diffusing cations during peridotite melting and melt-rock reaction. Apparent clinopyroxene-melt partition coefficients for REE, Sr, Y, Ti, and Zr, measured from reprecipitated clinopyroxene and coexisting melt in the reactive boundary layer, approach their equilibrium values reported in the literature. Disequilibrium melting models based on volume diffusion in solid limited mechanism are likely to significantly underestimate the rates at which major and trace elements in residual minerals reequilibrate with their surrounding melt. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Fast kimberlite ascent rates estimated from hydrogen diffusion profiles in xenolithic mantle olivines from southern Africa

Fourier transform infrared spectrometry (FTIR) analyses of olivines from peridotite xenoliths found in southern African kimberlites indicate 0 to 80 ppm H₂O concentrations. OH absorbance profiles across olivine grains show homogeneous H contents from core to edge for most samples. In one sample the olivines are H-free, while another has olivines characterized by lower H contents at the grain edges compared to the cores, indicating H loss during transport of the xenolith to the surface. Flat or near-flat H profiles place severe constraints on the duration of H loss from olivine grains, with implications for kimberlite magma ascent rates. Diffusion equations were used to estimate times of H loss of about 4 h for the sample with heterogeneous olivine H contents. Resulting kimberlite ascent rates are calculated to be 5-37 m s⁻¹ minimum, although these estimates are highly dependent on volatile contents and degassing behavior of the host kimberlite magma. Xenolithic olivines from alkali basalts generally have lower H contents and more pronounced H diffusion profiles than do those from kimberlites. This difference is likely caused by higher magma temperatures and lower ascent rates of alkali basalts compared to kimberlites.     

A depleted, not ideally chondritic bulk Earth: The explosive-volcanic basalt loss hypothesis

It has long been customary to assume that in the bulk composition of the Earth, all refractory-lithophile elements (including major oxides Al₂O₃ and CaO, all of the REE, and the heat-producing elements Th and U) occur in chondritic, bulk solar system, proportion to one another. Recently, however, Nd-isotopic studies (most notably Boyet M. and Carlson R. W. (2006) A new geochemical model for the Earth’s mantle inferred from ¹⁴⁶Sm-¹⁴²Nd systematics. Earth Planet. Sci. Lett.250, 254-268) have suggested that at least the outer portion of the planet features a Nd/Sm ratio depleted to ∼0.93 times the chondritic ratio. The primary reaction to this type of evidence has been to invoke a “hidden” reservoir of enriched matter, sequestered into the deepest mantle as a consequence of primordial differentiation. I propose a hypothesis that potentially explains the evidence for Nd/Sm depletion in a very different way. Among the handful of major types of differentiated asteroidal meteorites, two (ureilites and aubrites) are ultramafic restites so consistently devoid of plagioclase that meteoriticists were once mystified as to how all the complementary plagioclase-rich matter (basalt) was lost. The explanation appears to be basalt loss by graphite-fueled explosive volcanism on roughly 100-km sized planetesimals; with the dispersiveness of the process dramatically enhanced, relative to terrestrial experience, because the pyroclastic gases expand into vacuous space (Wilson L. and Keil K. (1991) Consequences of explosive eruptions on small Solar System bodies: the case of the missing basalts on the aubrite parent body. Earth Planet. Sci. Lett.104, 505-512). By analogy with lunar pyroclastic products, the typical size of pyroclastic melt/glass droplets under these circumstances will be roughly 0.1 mm. Once separated from an asteroidal or planetesimal gravitational field, droplets of this size will generally spiral toward the Sun, rather than reaccrete, because drag forces such the Poynting-Robertson effect quickly modify their orbits (the semimajor axis, in a typical scenario, is reduced by several hundred km during the first trip around the Sun). Assuming a similar process occurred on many of the Earth’s precursor planetesimals while they were still roughly 100 km in diameter, the net effect would be a depleted composition for the final Earth. I have modeled the process of trace-element depletion in the planetesimal mantles, assuming the partial melting was nonmodal and either batch or dynamic in terms of the melt-removal style. Assuming the process is moderately efficient, typical final-Earth Nd/Sm ratios are 0.93-0.96 times chondritic. Depletion is enhanced by a relatively low assumed residual porosity in batch-melting scenarios, but dampened by a relatively high value for “continuous” residue porosity in dynamic melting scenarios. Pigeonite in the source matter has a dampening effect on depletion. There are important side effects to the Nd/Sm depletion. The heat-producing elements, Th, U and K, might be severely depleted. The Eu/Eu^∗ ratio of the planet is unlikely to remain precisely chondritic. One of the most inevitable side effects, depletion of the Al/Ca ratio, is consistent with an otherwise puzzling aspect of the composition of the upper mantle. A perfectly undepleted composition for the bulk Earth is dubious.     

Thermal evolution and physics of melt extraction on the ureilite parent body

We develop a physical model of the thermal history of the ureilite parent body (UPB) that numerically tracks the history of its heating, hydration, dehydration, partial melting and smelting as a function of its formation time and the initial values of its composition, formation temperature and water ice content. Petrologic and chemical data from the main group (non-polymict) ureilite meteorites, which sample the interior of the UPB between depths corresponding to pressures in the range 3-10 MPa, are used to constrain the model. We find that to achieve the ∼30% melting inferred for ureilites from all sampled depths, the UPB must have had a radius between ∼80 and ∼130 km and must have accreted about 0.55 Ma after CAI formation. Melting began in the body at ∼1 Ma after CAI, and the time at which 30% melting was reached varied with depth in the asteroid but was always between ∼4.5 and ∼5.8 Ma after CAI. The total rate at which melt was produced in the UPB varied from more than 100 m³ s⁻¹ in the very early stages of melting at ∼1 Ma after CAI to ∼5 m³ s⁻¹ between 2 and 3 Ma after CAI, decreasing to extremely small values as the end of melting was approached beyond ∼5 Ma. Although the initial period of high melt production occupied only a short time around 1 Ma after CAI, it corresponded to ∼half (16%) of total silicate melting, and all strictly basaltic (i.e. plagioclase-saturated) melts must have been produced during this period.A very efficient melt transport network, consisting of a hierarchy of veins and larger pathways (dikes), developed quickly at the start of melting, ensuring rapid (timescales of months) transport of any single parcel of melt to shallow levels, thus ensuring that chemical interaction between melts and the rocks through which they subsequently passed was negligible. Volatile (mainly carbon monoxide) production due to smelting began at the start of silicate melting in the shallowest parts of the UPB and at later times at greater depths. Except at the very start and very end of melting, the volatile content of the melts produced was always high - generally between 15 and 35 mass % - and most of the melt produced was erupted at the surface of the UPB with speeds well in excess of the escape velocity and was lost into space. However, we show that 30% melting at the 3 MPa pressure level was only possible if ∼15% of the total melt produced in the asteroid was retained as a small number (∼5) of very extensive, sill-like intrusions centered at a depth of ∼7 km below the surface, near the base of the ∼8 km thick outer crust of the asteroid that was maintained at temperatures below the basalt solidus by conductive heat loss to the surface. The horizontal extents of these sills occupied about 75% of the surface area of the UPB, and the sills acted as buffers between the steady supply of melt from depth and the intermittent explosive eruption of the melt into space. We infer that samples from these intrusions are preserved as the rare feldspathic (loosely basaltic) clasts in polymict ureilites, and show that the cooling histories of the sills are consistent with these clasts reaching isotopic closure at ∼5 Ma after CAI, as given by ²⁶Al-²⁶Mg, ⁵³Mn-⁵³Cr and Pb-Pb age dates.