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.     

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.     

The impact of summer temperature changes on vegetation development in Ireland during the Weichselian Lateglacial Interstadial

We analysed pollen from a sediment core from Fiddaun, a small Lateglacial lake basin in western Ireland. Results reflect the general Lateglacial vegetation development in Ireland, as reconstructed from other pollen records. The Fiddaun diagram shows a number of short‐lived regressive vegetation phases during the Interstadial. The close similarity between two pollen records from the same region (Fiddaun and Lurga) indicates that these fluctuations probably reflect regional rather than local changes. Comparison with a previously published climate reconstruction, based on a chironomid‐inferred mean July air temperature reconstruction, lithology, and oxygen and carbon isotopes of lake marl from the Fiddaun record, allowed us to establish the relationship between summer temperature and vegetation changes. Results reveal that two temporary regressive shifts in the pollen record correspond to cold oscillations, which have been correlated to Greenland Interstadial 1b and 1d. It seems that the first cold oscillation (GI‐1d) had the most distinct effect on vegetation in Ireland. In contrast, it appears that the transition from Juniperus shrubland and Empetrum heath to grassland, which is estimated at ～13.7 ka BP, was not caused by decreasing summer temperatures, as no substantial change is observed in the climate proxies. Copyright © 2012 John Wiley & Sons, Ltd.     

The Influence of Path Corrections and a Three-dimensional Global <Emphasis Type="Italic">P</Emphasis>-wave Velocity Model on Seismic Event Location in Kazakhstan

We relocate 81 large nuclear explosions that were detonated at the Balapan and Degelen Mountain subregions of the Semipalatinsk test site in Kazakhstan during the years 1978 to 1989. The absolute locations of these explosions are available, as well as very accurate estimates of their origin times. This ground truth information allows us to perform a detailed analysis of location capability. We use a sparse network of stations with highly accurate first arrival picks measured using a waveform cross-correlation method. These high quality data facilitate very accurate location estimates with only a few phases per event. We contrast two different approaches: 1) a calibration-based approach, where we achieve improved locations by using path corrections, and 2) a model-based approach, where we achieve improved locations by relocating in a recently published global 3-D P-velocity model. Both methods result in large improvements in accuracy of the obtained absolute locations, compared to locations obtained in a 1-D reference earth model (ak135). The calibration-based approach gives superior results for this test site, in particular when arrival times from regional stations are included. Estimated locations remain well within a 1000 km² region surrounding the ground truth locations when the path corrections for the Balapan and Degelen Mountain subregions are interchanged, but even for the short separation between these two regions, we find variations in the path corrections that cause systematic mislocations. The model-based approach also results in substantially reduced mislocation distances and has the distinct advantage that it is, in principle, transportable to other source regions around the world.     

Effects of an oil spill in outdoor model tidal flat ecosystems

For several reasons the fate and effects of oil pollution are not easily studied in natural ecosystems. Furthermore, the results of laboratory tests cannot easily be extrapolated to natural systems. A logical objective of ecotoxicological research is, therefore, to simulate natural systems in model ecosystem experiments. A one-year feasibility study was carried out with MOdel TIdal Flat ecosystems (MOTIFs), designed to represent tidal flat systems typical of the Wadden Sea and other temperate estuarine areas. During the first 6 months the 4 MOTIFs studied developed very similarly. After this period an oil spill was simulated by exposing two of the four MOTIFs to a floating oil mousse for one week. During this week part of the oil was buried in the sediment by bioturbation. This sediment-bound fraction caused a prolonged exposure of the MOTIFs to oil compounds. Short- and long-term effects resulted in large differences between oil-treated and control MOTIFs. These differences persisted to the end of the experiment, 7 months after the removal of most of the oil.     

A map-view restoration of the Alpine-Carpathian-Dinaridic system for the Early Miocene

A map-view palinspastic restoration of tectonic units in the Alps, Carpathians and Dinarides reveals the plate tectonic configuration before the onset of Miocene to recent deformations. Estimates of shortening and extension from the entire orogenic system allow for a semi-quantitative restoration of translations and rotations of tectonic units during the last 20 Ma. Our restoration yielded the following results: (1) The Balaton Fault and its eastern extension along the northern margin of the Mid-Hungarian Fault Zone align with the Periadriatic Fault, a geometry that allows for the eastward lateral extrusion of the Alpine-Carpathian-Pannonian (ALCAPA) Mega-Unit. The Mid-Hungarian Fault Zone accommodated simultaneous strike-perpendicular shortening and strike-slip movements, concomitant with strike-parallel extension. (2) The Mid-Hungarian Fault Zone is also the locus of a former plate boundary transforming opposed subduction polarities between Alps (including Western Carpathians) and Dinarides. (3) The ALCAPA Mega-Unit was affected by 290 km extension and fits into an area W of present-day Budapest in its restored position, while the Tisza-Dacia Mega-Unit was affected by up to 180 km extension during its emplacement into the Carpathian embayment. (4) The external Dinarides experienced Neogene shortening of over 200 km in the south, contemporaneous with dextral wrench movements in the internal Dinarides and the easterly adjacent Carpatho-Balkan orogen. (5) N–S convergence between the European and Adriatic plates amounts to some 200 km at a longitude of 14° E, in line with post-20 Ma subduction of Adriatic lithosphere underneath the Eastern Alps, corroborating the discussion of results based on high-resolution teleseismic tomography.The displacement of the Adriatic Plate indenter led to a change in subduction polarity along a transect through the easternmost Alps and to substantial Neogene shortening in the eastern Southern Alps and external Dinarides. While we confirm that slab-pull and rollback of oceanic lithosphere subducted beneath the Carpathians triggered back-arc extension in the Pannonian Basin and much of the concomitant folding and thrusting in the Carpathians, we propose that the rotational displacement of this indenter provided a second important driving force for the severe Neogene modifications of the Alpine-Carpathian-Dinaridic orogenic system.     

The CM carbonaceous chondrite regolith Diepenveen

A carbonaceous chondrite was recovered immediately after the fall near the village of Diepenveen in the Netherlands on October 27, 1873, but came to light only in 2012. Analysis of sodium and poly‐aromatic hydrocarbon content suggests little contamination from handling. Diepenveen is a regolith breccia with an overall petrology consistent with a CM classification. Unlike most other CM chondrites, the bulk oxygen isotopes are extremely ¹⁶O rich, apparently dominated by the signature of anhydrous minerals, distributed on a steep slope pointing to the domain of intrinsic CM water. A small subset plots closer to the normal CM regime, on a parallel line 2 ‰ lower in δ¹⁷O. Different lithologies in Diepenveen experienced varying levels of aqueous alteration processing, being less aqueously altered at places rather than more heated. The presence of an agglutinate grain and the properties of methanol‐soluble organic compounds point to active impact processing of some of the clasts. Diepenveen belongs to a CM clan with ~5 Ma CRE age, longer than most other CM chondrites, and has a relatively young K‐Ar resetting age of ~1.5 Ga. As a CM chondrite, Diepenveen may be representative of samples soon to be returned from the surface of asteroid (162173) Ryugu by the Hayabusa2 spacecraft.     

The impact of surface‐wave attenuation on 3‐D seismic survey design

Three‐dimensional seismic survey design should provide an acquisition geometry that enables imaging and amplitude‐versus‐offset applications of target reflectors with sufficient data quality under given economical and operational constraints. However, in land or shallow‐water environments, surface waves are often dominant in the seismic data. The effectiveness of surface‐wave separation or attenuation significantly affects the quality of the final result. Therefore, the need for surface‐wave attenuation imposes additional constraints on the acquisition geometry. Recently, we have proposed a method for surface‐wave attenuation that can better deal with aliased seismic data than classic methods such as slowness/velocity‐based filtering. Here, we investigate how surface‐wave attenuation affects the selection of survey parameters and the resulting data quality. To quantify the latter, we introduce a measure that represents the estimated signal‐to‐noise ratio between the desired subsurface signal and the surface waves that are deemed to be noise. In a case study, we applied surface‐wave attenuation and signal‐to‐noise ratio estimation to several data sets with different survey parameters. The spatial sampling intervals of the basic subset are the survey parameters that affect the performance of surface‐wave attenuation methods the most. Finer spatial sampling will reduce aliasing and make surface‐wave attenuation easier, resulting in better data quality until no further improvement is obtained. We observed this behaviour as a main trend that levels off at increasingly denser sampling. With our method, this trend curve lies at a considerably higher signal‐to‐noise ratio than with a classic filtering method. This means that we can obtain a much better data quality for given survey effort or the same data quality as with a conventional method at a lower cost.     

Power law or power flaw?

Since its introduction by Svensson in 1959, the power law curve y = ax^b (where x and y are horizontal and vertical direction, respectively) has been widely used in morphological analysis of glacial trough cross-profiles. The numerical constants a and b are obtained by a linear regression analysis of the logarithmic form of the power law curve (ln y = ln a + b ln x). The value b then gives a measure for the form of the cross-section. However, over the years this form of the power law has endured a lot of criticism. This criticism is well founded, since this particular form of the power law is not suitable for curve fitting in morphological analyses. In this paper a general power law is proposed, of the form y − y₀ = a|x − x₀|^b (where x₀, y₀ are the coordinates of the origin of the cross-profile). A unique and unbiased solution for this equation is obtained with a general least-squares method, thereby minimizing the error between the cross-profile data and the curve, and not between the logarithmic transform of the data and its regression line. This provides a robust way to characterize trough cross-section forms. © 1998 John Wiley & Sons, Ltd.