Experimental and computational study of trace element distribution between orthopyroxene and anhydrous silicate melt: substitution mechanisms and the effect of iron

Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 Although orthopyroxene (Opx) is present during a wide range of magmatic differentiation processes in the terrestrial and lunar mantle, its effect on melt trace element contents is not well quantified. We present results of a combined experimental and computational study of trace element partitioning between Opx and anhydrous silicate melts. Experiments were performed in air at atmospheric pressure and temperatures ranging from 1,326 to 1,420°C in the system CaO–MgO–Al₂O₃–SiO₂ and subsystem CaO–MgO–SiO₂. We provide experimental partition coefficients for a wide range of trace elements (large ion lithophile: Li, Be, B, K, Rb, Sr, Cs, Ba, Th, U; rare earth elements, REE: La, Ce, Nd, Sm, Y, Yb, Lu; high field strength: Zr, Nb, Hf, Ta, Ti; transition metals: Sc, V, Cr, Co) for use in petrogenetic modelling. REE partition coefficients increase from $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 $ D_{\text{La}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.0005 to $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 $ D_{\text{Lu}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.109 , D values for highly charged elements vary from $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 $ D_{\text{Th}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0026 through $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 $ D_{\text{Nb}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0033 and $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 $ D_{\text{U}}^{{{\text{Opx}} {\hbox{-}} {\text{melt}}}} \sim 0.0066 to $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 $ D_{\text{Ti}}^{{\text{Opx}} {\hbox{-}} {\text{melt}}} \sim 0.058 , and are all virtually independent of temperature. Cr and Co are the only compatible trace elements at the studied conditions. To elucidate charge-balancing mechanisms for incorporation of REE into Opx and to assess the possible influence of Fe on Opx-melt partitioning, we compare our experimental results with computer simulations. In these simulations, we examine major and minor trace element incorporation into the end-members enstatite (Mg₂Si₂O₆) and ferrosilite (Fe₂Si₂O₆). Calculated solution energies show that R²⁺ cations are more soluble in Opx than R³⁺ cations of similar size, consistent with experimental partitioning data. In addition, simulations show charge balancing of R³⁺ cations by coupled substitution with Li⁺ on the M1 site that is energetically favoured over coupled substitution involving Al–Si exchange on the tetrahedrally coordinated site. We derived best-fit values for ideal ionic radii r ₀, maximum partition coefficients D ₀, and apparent Young’s moduli E for substitutions onto the Opx M1 and M2 sites. Experimental r ₀ values for R³⁺ substitutions are 0.66–0.67 ? for M1 and 0.82–0.87 ? for M2. Simulations for enstatite result in r ₀ = 0.71–0.73 ? for M1 and ~0.79–0.87 ? for M2. Ferrosilite r ₀ values are systematically larger by ~0.05 ? for both M1 and M2. The latter is opposite to experimental literature data, which appear to show a slight decrease in $ r_{0}^{{{\text{M}}2}} $ r_{0}^{{{\text{M}}2}} in the presence of Fe. Additional systematic studies in Fe-bearing systems are required to resolve this inconsistency and to develop predictive Opx-melt partitioning models for use in terrestrial and lunar magmatic differentiation models.     

Penetrometry of granular and moist planetary surface materials: Application to the Huygens landing site on Titan

The Huygens probe landed on the then unknown surface of Titan in January 2005. A small, protruding penetrometer, part of the Surface Science Package (SSP), was pushed into the surface material measuring the mechanical resistance of the ground as the probe impacted the landing site. We present laboratory penetrometry into room temperature surface analogue materials using a replica penetrometer to investigate further the nature of Titan’s surface and examine the sensor’s capabilities. The results are then compared to the flight instrument’s signature and suggest the Titan surface substrate material consists of sand-sized particles with a mean grain size ～2 mm. A possible thin 7 mm coating with mechanical properties similar to terrestrial snow may overlie this substrate, although due to the limited data we are unable to detect any further layering or grading within the near-surface material. The unusual weakening with depth of the signature returned from Titan has, to date, only been reproduced using a damp sand target that becomes progressively wetter with depth, and supports the suggestion that the surface may consist of a damp and cohesive material with interstitial liquid contained between its grains. Comparison with terrestrial analogues highlights the unusual nature of the landing site material.     

Investigation of runoff generation in a pristine,poorly gauged catchment in the Chilean Andes II: Qualitative and quantitative use of tracers at three spatial scales

Understanding runoff generation processes is important for flood prediction, water management, erosion control, water quality, contaminant transport and the evaluation of impacts of land use change. However, little process research has been carried out in southern Chile. In particular the young volcanic ash soils, which are typical for this area, are not well understood in their hydrologic behaviour. To establish a ‘reference study’ which can then be used for comparison with other (disturbed) sites, this study focuses on the investigation of runoff generation processes in an undisturbed, forested catchment in the Chilean Andes. The paper reports on an investigation of these processes with different tracer methods at different spatial scales. Hydrograph separation with environmental isotopes and geochemical constituents was used on the catchment scale. Thermal energy was used as a tracer to investigate groundwater–surface water interactions at the local stream reach scale and dye tracers were used to study infiltration and percolation characteristics at the plot scale. It was found that pre‐event water dominates the storm hydrograph. In the lower reaches, however, water usually exfiltrates from the stream into the adjacent aquifer. The dye tracer experiments showed that while preferential vertical flow dominates under forest, water infiltrates as a straight horizontal front in the bare volcanic ashes (no vegetation) on the catchment rim. Subsurface flow patterns in the forest differ significantly from summer to winter. All three approaches used in this study suggest an important shift in dominant processes from dry to wet season. Copyright © 2008 John Wiley & Sons, Ltd.     

Boron and Oxygen Isotope Composition of Certified Reference Materials NIST SRM 610/612 and Reference Materials JB-2 and JR-2

We present data on the concentration, the isotope composition and the homogeneity of boron in NIST silicate glass reference materials SRM 610 and SRM 612, and in powders and glasses of geological reference materials JB-2 (basalt) and JR-2 (rhyolite). Our data are intended to serve as references for both microanalytical and wet-chemical techniques. The δ¹¹ B compositions determined by N-TIMS and P-TIMS agree within 0.5% and compare with SIMS data within 2.5%. SIMS profiles demonstrate boron isotope homogeneity to better than δ¹¹ B = 2% for both NIST glasses, however a slight boron depletion was detected towards the outermost 200 μm of the rim of each sample wafer. The boron isotope compositions of SRM 610 and SRM 612 were indistinguishable. Glasses produced in this study by fusing JB-2 and JR-2 powder also showed good boron isotope homogeneity, both within and between different glass fragments. Their major element abundance as well as boron isotope compositions and concentrations were identical to those of the starting composition. Hence, reference materials (glasses) for the in situ measurement of boron isotopes can be produced from already well-studied volcanic samples without significant isotope fractionation. Oxygen isotope ratios, both within and between wafers, of NIST reference glasses SRM 610 and SRM 612 are uniform. In contrast to boron, significant differences in oxygen isotope compositions were found between the two glasses, which may be due to the different amounts of trace element oxides added at ten-fold different concentration levels to the silicate matrix.     

Reviews: Monitoring the Comprehensive Nuclear-Test-Ban Treaty: Hydroacoustics. Edited by Catherine deGroot-Hedlin and John Orcutt. Birkhäuser Verlag Basel-Boston-Berlin, 2001, Pure and Applied Geophysics (PAGEOPH), Volume 158 (2001), No. 3, 212 pages. Softcover. ISBN 3-7643-6538-2, sFr. 68,- / DM 90,- / öS 657,-.

Studia Geophysica et Geodaetica -     

Earth orientation determinations by short duration VLBI observations

Summary In May 1989 and April 1990 the radio telescopes of the Wettzell Geodetic Fundamental Station in Germany and of the Shanghai Observatory near Seshan in China observed two series of daily VLBI experiments of short duration for precise determination of UT1. In 1990 a few experiments were complemented by the Hartebeesthoek Radio Astronomy Observatory in South Africa and the Medicina telescope of the Bologna Istituto di Radioastronomia in Italy. Employing the South African station together with the east-west baseline formed by the observatories of Seshan and Medicina permitted simultaneous determinations of UT1 and polar motion. Here we report on the results of these observations. Comparing the UT1 results with those of the IRIS Intensive series gives a clear indication of the absolute accuracy of such short duration VLBI measurements which is estimated to be of the order of ±60µs.     

ДИНАМИЧЕСКИЕ ПАРАМЕТРЫ СИЛЬНЫШ ЗЕМЛЕТРЯСЕНИЙ ПО ДАННЫМ ШИРОКОПОЛОСНОЙ ДЛИННОПЕРИОДНОЙ АППАРАТУРЫ

Summary The dynamic focal parameters of a series of strong earthquakes, recorded by the FBV broadband velocigraph at KHC (Kaperské Hory, Czechoslovakia), are estimated. The influence of specific wave pattern features on the estimation accuracy is studied. Polarization analysis is used to demonstrate the complicated interference character even of record portions corresponding to a single teleseismic wave group. The values of the dynamic focal parameters do not depend significantly on the length of the analysed record portions and are sufficiently consistent with those reported by Riznichenko [2]. The possibility of estimating dynamic focal parameters from recording of standard classes B (intermediate band) and C (long-period) seismographs is pointed out.     

The FBV broadband seismograph station profile KHC-KSP

nua naam u nuu amamuu un uu mau an () u (), umua nu ¶rt; ¶rt;a¶rt; nu 7 u aua ¶rt;a¶rt; m¶rt;um aa naa u . a ¶rt;a nu amua ¶rt;m un ¶rt; u¶rt;au mmu am aua u a au uu uu mu.