Feldspathic clast populations in polymict ureilites: Stalking the missing basalts from the ureilite parent body

Polymict ureilites DaG 164/165, DaG 319, DaG 665, and EET 83309 are regolith breccias composed mainly of monomict ureilite-like material, but containing ∼2 vol% of feldspathic components. We characterized 171 feldspathic clasts in these meteorites in terms of texture, mineralogy, and mineral compositions. Based on this characterization we identified three populations of clasts, each of which appears to represent a common igneous (generally basaltic) lithology and whose mafic minerals show a normal igneous fractionation trend of near-constant Fe/Mn ratio over a range of Fe/Mg ratios that extend to much higher values than those in monomict ureilites. The melts represented by these populations are unlikely to be impact melts, because the ubiquitous presence of carbon in polymict ureilites (the regolith of the ureilite parent body) implies that impact melts would have crystallized under conditions of carbon redox control and therefore have highly magnesian mafic mineral compositions with constant Mn/Mg ratio. Therefore, these melts appear to be indigenous products of igneous differentiation on the ureilite parent body (UPB), complementary to the olivine-pigeonite residues represented by the majority of monomict ureilites.The most abundant population is characterized by albitic plagioclase in association with pyroxenes, phosphates, ilmenite, silica, and incompatible-element enriched glass. Model calculations suggest that it formed by extensive fractional crystallization of the earliest melt(s) of precursor materials from which the most magnesian (shallowest) olivine-pigeonite ureilites formed. A less abundant population, characterized by labradoritic plagioclase, may have formed from melts complementary to more ferroan olivine-pigeonite ureilites, and derived from deeper in the UPB. The third population, characterized by the presence of olivine and augite, could only have formed from melts produced at greater depths in the UPB than the olivine-pigeonite ureilites. Many other feldspathic clasts cannot be positively associated with any of these three populations, because their mafic mineral compositions exhibit carbon redox control. However, they may be products of early crystallization of basaltic melts produced on the UPB, before carbon was exhausted by reduction.Partial melting on the ureilite parent body was a fractional (or incremental) process. Melts were produced early in UPB history, and most likely extracted rapidly, thus preserving primitive chemical and oxygen isotopic signatures in the residues.     

Prospects for the William Herschel Telescope

René Rutten, Director of the Isaac Newton Group of Telescopes, sets the scene for future developments in international astronomy on La Palma.     

Chemical and carbon isotopic evolution of hydrocarbons during prograde metamorphism from 100°C to 550°C: Case study in the Liassic black shale formation of Central Swiss Alps

Hydrocarbon distributions and stable isotope ratios of carbonates (δ¹³C_car, δ¹⁸O_car), kerogen (δ¹³C_ker), extractable organic matter (δ¹³C_EOM) and individual hydrocarbons of Liassic black shale samples from a prograde metamorphic sequence in the Swiss Alps were used to identify the major organic reactions with increasing metamorphic grade. The studied samples range from the diagenetic zone (<100°C) to amphibolite facies (∼550°C). The samples within the diagenetic zones (<100 and 150°C) are characterized by the dominance of C_<20n-alkanes, suggesting an origin related with marine and/or bacterial inputs. The metamorphic samples (200 to 550°C) have distributions significantly dominated by C₁₂ and C₁₃n-alkanes, C₁₄, C₁₆ and C₁₈n-alkylcyclopentanes and to a lesser extend C₁₅, C₁₇ and C₂₁n-alkylcyclohexanes. The progressive ¹³C-enrichment (up to 3.9‰) with metamorphism of the C_>17n-alkanes suggests the occurrence of cracking reactions of high molecular weight compounds. The isotopically heavier (up to 5.6‰) C_<17n-alkanes in metamorphic samples are likely originated by thermal degradation of long-chain homologous with preferential release of isotopically light C₁ and C₂ radicals. The dominance of specific even C-number n-alkylcyclopentanes suggests an origin related to direct cyclization mechanism (without decarboxylation step) of algal or bacterial fatty acids occurring in reducing aqueous metamorphic fluid conditions. The regular increase of the concentrations of n-alkylcycloalkanes vs. C_>13n-alkanes with metamorphism suggests progressive thermal release of kerogen-linked fatty acid precursors and degradation of n-alkanes. Changes of the steroid and terpenoid distributions are clearly related to increasing metamorphic temperatures. The absence of 18α(H)-22,29,30-trisnorneohopane (Ts), the occurrence of 17β(H)-trisnorhopane, 17β(H), 21α(H)-hopanes in the C₂₉ to C₃₁ range and 5α(H),14α(H),17α(H)-20R C₂₇, C₂₉ steranes in the low diagenetic samples (<100°C) are characteristic of immature bitumens. The higher thermal stress within the upper diagenetic zone (150°C) is marked by the presence of Ts, the disappearance of 17β(H)-trisnorhopane and thermodynamic equilibrium of the 22S/(22S + 22R) homohopane ratios. The increase of the ααα-sterane 20S/(20S + 20R) and 20R ββ/(ββ + αα) ratios (from 0.0 to 0.55 and from 0.0 to 0.40, respectively) in the upper diagenetic zone indicates the occurrence of isomerization reactions already at <150°C. However, the isomerization at C-20 (R → S) reaches thermodynamic equilibrium values already at the upper diagenesis (∼150°C) whereas the epimerisation at C-14 and C-17 (αα → ββ) arrives to constant values in the lower anchizone (∼200°C). The ratios Ts vs. 17α(H)-22,29,30-trisnorneohopane [(Ts/(Ts + Tm)] and 18α(H)-30-norneohopane (C₂₉Ts) vs. 17α(H),21β(H)-30-norhopane [C₂₉Ts/(C₂₉Ts + C₂₉)] increase until the medium anchizone (200 to 250°C) from 0.0 to 0.96 and from 0.0 to 0.44, respectively. An opposite trend towards lower values is observed in the higher metamorphic samples.The occurrence of specific hydrocarbons (e.g., n-alkylcyclopentanes, cadalene, hydrogenated aromatic compounds) in metamorphic samples points to kerogen degradation reactions most probably occurring in the presence of water and under reducing conditions. The changes of hydrocarbon distributions and carbon isotopic compositions of n-alkanes related to metamorphism suggest that the organic geochemistry may help to evaluate the lowest grades of prograde metamorphism.     

Numerical modeling of toxic nonaqueous phase liquid removal from contaminated groundwater systems: mesh effect and discretization error estimation

Numerical modeling has now become an indispensable tool for investigating the fundamental mechanisms of toxic nonaqueous phase liquid (NAPL) removal from contaminated groundwater systems. Because the domain of a contaminated groundwater system may involve irregular shapes in geometry, it is necessary to use general quadrilateral elements, in which two neighbor sides are no longer perpendicular to each other. This can cause numerical errors on the computational simulation results due to mesh discretization effect. After the dimensionless governing equations of NAPL dissolution problems are briefly described, the propagation theory of the mesh discretization error associated with a NAPL dissolution system is first presented for a rectangular domain and then extended to a trapezoidal domain. This leads to the establishment of the finger‐amplitude growing theory that is associated with both the corner effect that takes place just at the entrance of the flow in a trapezoidal domain and the mesh discretization effect that occurs in the whole NAPL dissolution system of the trapezoidal domain. This theory can be used to make the approximate error estimation of the corresponding computational simulation results. The related theoretical analysis and numerical results have demonstrated the following: (1) both the corner effect and the mesh discretization effect can be quantitatively viewed as a kind of small perturbation, which can grow in unstable NAPL dissolution systems, so that they can have some considerable effects on the computational results of such systems; (2) the proposed finger‐amplitude growing theory associated with the corner effect at the entrance of a trapezoidal domain is useful for correctly explaining why the finger at either the top or bottom boundary grows much faster than that within the interior of the trapezoidal domain; (3) the proposed finger‐amplitude growing theory associated with the mesh discretization error in the NAPL dissolution system of a trapezoidal domain can be used for quantitatively assessing the correctness of computational simulations of NAPL dissolution front instability problems in trapezoidal domains, so that we can ensure that the computational simulation results are controlled by the physics of the NAPL dissolution system, rather than by the numerical artifacts. Copyright © 2014 John Wiley & Sons, Ltd.