Mixing of magmas from enriched and depleted mantle sources in the northeast Pacific: West Valley segment,Juan de Fuca Ridge

The 50 km-long West Valley segment of the northern Juan de Fuca Ridge is a young, extension-dominated spreading centre, with volcanic activity concentrated in its southern half. A suite of basalts dredged from the West Valley floor, the adjacent Heck Seamount chain, and a small near-axis cone here named Southwest Seamount, includes a spectrum of geochemical compositions ranging from highly depleted normal (N-) MORB to enriched (E-) MORB. Heck Seamount lavas have chondrite-normalized La/Sm_cn0.3, ⁸⁷Sr/⁸⁶Sr=0.70235–0.70242, and ²⁰⁶Pb/²⁰⁴Pb=18.22–18.44, requiring a source which is highly depleted in trace elements both at the time of melt generation and over geologic time. The E-MORB from Southwest Seamount have La/Sm_cn1.8, ⁸⁷Sr/⁸⁶Sr=0.70245–0.70260, and ²⁰⁶Pb/²⁰⁴Pb=18.73–19.15, indicating a more enriched source. Basalts from the West Valley floor have chemical compositions intermediate between these two end-members. As a group, West Valley basalts from a two-component mixing array in element-element and element-isotope plots which is best explained by magma mixing. Evidence for crustal-level magma mixing in some basalts includes mineral-melt chemical and isotopic disequilibrium, but mixing of melts at depth (within the mantle) may also occur. The mantle beneath the northern Juan de Fuca Ridge is modelled as a plum-pudding, with plums of enriched, amphibole-bearing peridotite floating in a depleted matrix (DM). Low degrees of melting preferentially melt the plums, initially removing only the amphibole component and producing alkaline to transitional E-MORB. Higher degrees of melting tap both the plums and the depleted matrix to yield N-MORB. The subtly different isotopic compositions of the E-MORBs compared to the N-MORBs require that any enriched component in the upper mantle was derived from a depleted source. If the enriched component crystallized from fluids with a DM source, the plums could evolve to their more evolved isotopic composition after a period of 1.5–2.0 Ga. Alternatively, the enriched component could have formed recently from fluids with a less-depleted source than DM, such as subducted oceanic crust. A third possibility is that enriched material might be dispersed as plums throughout the upper mantle, transported from depth by mantle plumes.     

Production of lunar fragmental material by meteoroid impact

The rate of production of new fragmental lunar surface material is derived theoretically on the hypothesis that such material is excavated from a bedrock layer by meteoroid impacts. An overlaying regolith effectively shields the bedrock layer from small impacts, reducing the production rate of centimeter-sized and smaller blocks by a large factor. Logarithmic production rate curves for centimeter to meter-sized blocks are nonlinear for any regolith from centimeters to tens of meters in thickness, with small blocks relatively much less frequent for thicker (older) regoliths, suggesting the possibility of a statistical reverse bedding. Modest variations in the exponents of scaling laws for crater depth-diameter ratio and maximum block-diameter to crater diameter ratio are shown to have significant effects on the production rates. The production rate increases slowly with increasing size of the largest crater affecting the region.     

Multidimensional GPR array processing using Kirchhoff migration

We compare the ability of several practical ground-penetrating radar (GPR) array processing methods to improve signal-to-noise ratio (SNR), increase depth of signal penetration, and suppress out-of-plane arrivals for data with SNR of roughly 1. The methods include two-dimensional (2-D) monostatic, three-dimensional (3-D) monostatic, and 3-D bistatic Kirchhoff migration. The migration algorithm is modified to include the radiation pattern for interfacial dipoles. Results are discussed for synthetic and field data. The synthetic data model includes spatially coherent noise sources that yield nonstationary signal statistics like those observed in high noise GPR settings. Array results from the model data clearly indicate that resolution and noise suppression performance increases as array dimensionality increases. Using 50-MHz array data collected on a temperate glacier (Gulkana Glacier, AK), we compare 2-D and 3-D monostatic migration results. The data have low SNR and contain reflections from a complex, steeply dipping bed. We demonstrate that the glacier bed can only be accurately localized with the 3-D array. In addition, we show that the 3-D array increases SNR (relative to a 2-D array) by a factor of three.     

Using isotopes to understand the evolution of water ages in disturbed mixed land-use catchments

Tracer studies have been key to unravelling catchment hydrological processes, yet most insights have been gained in environments with relatively low human impact. We investigated the spatial variability of stream isotopes and water ages to infer dominant flow paths in a ~10-km² nested catchment in a disturbed, predominantly agricultural environment in Scotland. We collected long-term (>5 years) stable isotope data of precipitation, artificial drainage, and four streams with varying soil and land use types in their catchment areas. Using a gamma model, Mean Transit Times (MTTs) were then estimated in order to understand the spatial variability of controls on water ages. Despite contrasting catchment characteristics, we found that MTTs in the streams were generally very similar and short (<1 year). MTTs of water in artificial drains were even shorter, ranging between 1 to 10 months for a typical field drain and <0.5 to 1 month for a country road drain. At the catchment scale, lack of heterogeneity in the response could be explained by the extensive artificial surface and subsurface drainage, “short-circuiting” younger water to the streams during storms. Under such conditions, additional intense disturbance associated with highway construction during the study period had no major effect on the stream isotope dynamics. Supplementary short-term (~14 months) sampling of mobile soil water in dominant soil-land use units also revealed that agricultural practices (ploughing of poorly draining soils and soil compaction due to grazing on freely draining soils) resulted in subtle MTT variations in soil water in the upper profile. Overall, the isotope dynamics and inferred MTTs suggest that the evolution of stream water ages in such a complex human-influenced environment are largely related to near-surface soil processes and the dominant soil management practices. This has direct implications for understanding and managing flood risk and contaminant transport in such environments.     

Street Versus Rooftop Level Concentrations of Fine Particles in a Cambridge Street Canyon

Dispersion of particles, as evidenced by changes in their number distributions (PNDs) and concentrations (PNCs), in urban street canyons, is still not well understood. This study compares measurements by a fast-response particle spectrometer (DMS500) of the PNDs and the PNCs (5–1000 nm, sampled at 1 Hz) at street and rooftop levels in a Cambridge UK street canyon, and examines mixing, physical and chemical conversion processes, and the competing influences of traffic volume and rooftop wind speed on the PNDs and the PNCs in various size ranges. PNCs at street level were ≈6.5 times higher than at rooftop. Street-level PNCs followed the traffic volume and decreased with increasing wind speed, showing a larger influence of wind speed on 30–300 nm particles than on 5–30 nm particles. Conversely, rooftop PNCs in the 5–30 nm size range increased with wind speed, whereas those for particles between 30 and 300 nm did not vary with wind speed.     

Titanium in subduction zone fluids: First insights from abinitio molecular metadynamics simulations

The preferential incorporation of High-Field-Strength Elements (HFSE) in rutile (TiO₂), combined with its supposed stability in subduction zone settings, make it an ideal candidate to explain the low HFSE concentrations in subduction-derived magmas. The solubility behaviour of rutile is key to these arguments, but at present experimental and field-based evidence are contradictory.We have used abinitio molecular (meta)dynamics to investigate the coordination environment of Ti(IV) in pure water at 300 and 1000 K and densities ranging from 900-1260 kg m⁻³ (approximate pressures 0.9-3.6 GPa). In all high temperature simulations, the long-range structure of the solvent indicates a breakdown of the hydrogen bonding network as expected for supercritical water. The five-fold coordination of titanium to water is energetically most favourable in aqueous fluids at room temperature and pressure, separated from four and six-fold configurations by ∼175 and ∼200 kJ mol⁻¹, respectively. The average first shell Ti-O distance is 2.00 Å, in excellent agreement with bond lengths obtained from experiments. At similar densities and 1000 K, titanium is on average six-fold coordinate with water, and shows some degree of water dissociation in the first hydration shell. This coordination environment is remarkably persistent with increasing density from 1021 to 1260 kg m⁻³ at constant temperature (1000 K). At lower densities, however, (900 kg m⁻³ at 1000 K), the coordination with first shell water molecules is less than five. The observed coordination changes could promote association of titanium with peralkaline or peraluminous domains in the aqueous fluid, and thereby explain field-and laboratory based evidence of enhanced HFSE concentrations.This study demonstrates that abinitio molecular dynamics has considerable potential to access details of element behaviour in aqueous fluids at geologically relevant conditions that are impossible to examine otherwise. Changes in the solvent structure due to density variations lead to differences in solvent behaviour allowing access to new domains for fluid-solid interaction. Moreover, changes in the solvent structure are strongly linked to the effectiveness of element solvation.     

Architecture of the Oman–UAE ophiolite: evidence for a multi-phase magmatic history

The Oman–United Arab Emirates ophiolite is the world’s largest ophiolite. It is divided into 12 separate fault-bounded blocks, of which the northern three lie wholly or partly in the United Arab Emirates. Extensive mapping has shown that the United Arab Emirates blocks contain mantle and crustal sections which correspond to the classic ‘Penrose conference’ ophiolite definition but which are cut by a voluminous later magmatic sequence including ultramafic, mafic and felsic components. Samples from the later magmatic sequence are dated at 96.4?±?0.3, 95.74?±?0.3 and 95.2?±?0.3 Ma; the early crustal section, which has not been dated directly, is thus constrained to be older than c. 96.4 Ma. Petrological evidence shows that the early crustal section formed at a spreading ridge, but the later magmatic sequence was formed from hydrous magmas that produced different mineral crystallisation sequences to normal mid-ocean ridge basalt (MORB). Mineral and whole-rock geochemical analyses show that the early crustal rocks are chemically similar to MORB, but the later magmatic sequence has chemical features typically found in supra-subduction zone (SSZ) settings. The ophiolite in the United Arab Emirates thus preserves clear evidence for two stages of magmatism, an early episode formed at a spreading centre and a later episode associated with the onset of subduction. Similar two-stage magmatism has been recognised in the Oman sector, but the United Arab Emirates contains the most voluminous SSZ magmatism yet described from this ophiolite.