Geothermal palaeogradients and metamorphic zonation from the conodont colour alteration index (CAI)

The conodont colour alteration index (CAI) is potentially valuable in thermal history investigations, but there are problems in its use. An approach is proposed herein to overcome the problem of establishing reliable geothermal palaeogradients in areas with burial metamorphism, based on knowledge of the stratigraphic succession thickness, the boundaries of chronostratigraphic units, and the age and duration of deformation episodes. Although also problematical that a correlation between CAI and Kübler index (KI) values does not exist, there are areas where CAI and KI values are notably consistent for the anchizone boundaries; however, anomalous correlations between CAI and KI values have also been found in some areas. A CAI anchizone (or ancaizone) is defined in order to enable CAI values to be used independently to establish a metamorphic zonation.     

Lead isotope systematics of Late Cretaceous – Tertiary Andean arc magmas and associated ores between 8°N and 40°S: evidence for latitudinal mantle heterogeneity beneath the Andes

ABSTRACT Lead isotope variability of magmatic arc rocks and associated mineralization of the Central Andes is usually considered to be the result of mixing between a homogeneous mantle and heterogeneous continental crust. About 230 new lead isotope data on the Northern and Central Andes allow us to compare for the first time lead isotope systematics of the Late Cretaceous – Tertiary arc magmatism and associated mineralization along the Andean chain between 8°N and 40°S. Lead isotope compositions indicate mixing between mantle and upper crustal rocks along the whole Andean chain. Additionally, we have found that mantle end-members of the Late Cretaceous – Tertiary magmatism are heterogeneous and systematically shifted towards less radiogenic ²⁰⁶Pb/²⁰⁴Pb compositions from north to south along the Andes. This heterogeneity most likely results from mixing between a low radiogenic mantle, possibly carrying a DMM or EM I component, and a more radiogenic mantle, possibly carrying an HIMU component. Thus, our results imply that lead isotope variability of Andean magmas at the continental scale is caused not only by crustal but also by mantle heterogeneity.     

Thermodynamic controls on element partitioning between titanomagnetite and andesitic–dacitic silicate melts

Titanomagnetite–melt partitioning of Mg, Mn, Al, Ti, Sc, V, Co, Ni, Cu, Zn, Ga, Zr, Nb, Mo, Hf and Ta was investigated experimentally as a function of oxygen fugacity (fO₂) and temperature (T) in an andesitic–dacitic bulk-chemical compositional range. In these bulk systems, at constant T, there are strong increases in the titanomagnetite–melt partitioning of the divalent cations (Mg²⁺, Mn²⁺, Co²⁺, Ni²⁺, Zn²⁺) and Cu²⁺/Cu⁺ with increasing fO₂ between 0.2 and 3.7 log units above the fayalite–magnetite–quartz buffer. This is attributed to a coupling between magnetite crystallisation and melt composition. Although melt structure has been invoked to explain the patterns of mineral–melt partitioning of divalent cations, a more rigorous justification of magnetite–melt partitioning can be derived from thermodynamic principles, which accounts for much of the supposed influence ascribed to melt structure. The presence of magnetite-rich spinel in equilibrium with melt over a range of fO₂ implies a reciprocal relationship between a(Fe²⁺O) and a(Fe³⁺O_1.5) in the melt. We show that this relationship accounts for the observed dependence of titanomagnetite–melt partitioning of divalent cations with fO₂ in magnetite-rich spinel. As a result of this, titanomagnetite–melt partitioning of divalent cations is indirectly sensitive to changes in fO₂ in silicic, but less so in mafic bulk systems.     

A new scheme for numerical modelling of flow and transport processes in 3D fractured porous media

The objective of this work is to develop a new numerical approach for the three-dimensional modelling of flow and transient solute transport in fractured porous media which would provide an accurate and efficient treatment of 3D complex geometries and inhomogeneities. For this reason, and in order to eliminate as much as possible the number of degrees of freedom, the fracture network, fractures and their intersections, are solved with a coupled 2D–1D model while the porous matrix is solved independently with a 3D model. The interaction between both models is accounted for by a coupling iterative technique. In this way it is possible to improve efficiency and reduce CPU usage by avoiding 3D mesh refinements of the fractures. The approach is based on the discrete-fracture model in which the exact geometry and location of each fracture in the network must be provided as an input. The formulation is based on a multidimensional coupling of the boundary element method-multidomain (BEM-MD) scheme for the flow and boundary element dual reciprocity method-multidomain (BE-DRM-MD) scheme for the transport. Accurate results and high efficiency have been obtained and are reported in this paper.     

Study of a rockfall in a limy conglomerate canyon (Covarrubias, Burgos, N. Spain)

In this paper, a rockfall of some tenths of cubic metres of rock, occurred in a conglomerate canyon in the northern Castilian plateau (N-Central Spain), is described in detail. This includes the basic data (geology, mapping, meteorological data and rock involved), the mechanical mechanism of the instability, the propagation and reach of the fallen rock boulders and the final location of the pieces. The study is based on detailed geological, aerial photography and geotechnical studies and advance topographical measurements. Geological engineering techniques are used to analyze the instability initiation and rockfall propagation. This study contributes to a better understanding of the geological processes leading to the occurrence of calcareous canyons and cliffs in the drainage network of the upper Castilian plateau. It is also helpful to improve our knowledge on natural rockfall phenomena for risk engineering purposes aiming to quantify, control and mitigate associated hazards.     

The Alleret maar (Massif Central, France): A new lacustrine sequence of the early Middle Pleistocene in western Europe

The Alleret maar (Massif Central, France) provides a long lacustrine sequence (40.6 m) attributed to the early Middle Pleistocene. Sediment, pollen and diatoms analysis of its upper part (AL2 core, 14.6 m) indicates two temperate phases marked by high lake levels, forest development and vegetation expansion. They are separated by a cold period during which lake level drops, coarse sediment input increases and steppic and xerophilous plants develop. Pollen data suggests that this sequence belongs to the upper part of the Cromerian complex. These results are in agreement with the 557 ± 3 ka (±12 ka, including all errors) ⁴⁰Ar/³⁹Ar age obtained from an interbedded tephra layer emitted by the Mont-Dore/Sancy strato-volcano and establish that this sequence probably covers the MIS 15 substages.     

Olivine in ultramafic lamprophyres: chemistry,crystallisation, and melt sources of Siberian Pre- and post-trap aillikites

We studied olivines from the Devonian pre-trap (the Ilbokich occurrence) and the Triassic post-trap (the Chadobets occurrence) carbonate-rich ultramafic lamprophyres (UMLs) in the southwestern portion of the Siberian craton. On the basis of detailed investigations of major, minor, and trace-element distributions, we have reconstructed the main processes that control the origins of these olivines. These include fractional crystallisation from melt, assimilation, and fractional crystallisation processes with orthopyroxene assimilation, melt-reaction diffusive re-equilibration, alkali enrichment, and CO₂ degassing of the melt. Furthermore, we inferred the composition of the sources of the primary UML melt and their possible correlations with proto-kimberlitic melts, as well as the influence of the Triassic Siberian plume on the composition of the lithospheric mantle. The main differences between olivines from the Ilbokich and the Chadobets aillikites were that the olivines from the former had more magnesium-rich cores (Mg# = 89.2?±?0.2), had Mg- and Cr-rich transition zones (Mg# = 89.7?±?0.2 and 300–500 ppm Cr), had lower Ni (up to 3100 ppm) and Li (1.4–1.5 ppm), and had higher B (0.8–2.6 ppm) contents, all at higher Fo values (90–86), relative to the olivines from the latter (Mg# = 88–75; 200–300 ppm Cr; up to 3400 ppm Ni; 1.4–2.4 ppm Li; 0.4–2.2 ppm B). The Siberian aillikite sources contained a significant amount of metasomatic material. Phlogopite-rich MARID-type veins provided the likely metasomatic component in the pre-trap Devonian Ilbokich aillikite source, whereas the Triassic Chadobets aillikitic post-trap melts were derived from a source with a significant carbonate component. A comparison of UML olivines with olivines from the pre-trap and post-trap Siberian kimberlites shows a striking similarity. This suggests that the carbonate component in the aillikitic source could have been produced by evolved kimberlite melts. The differences in the lithospheric metasomatic component that contributed to pre-trap and post-trap aillikitic melts can be interpreted as reflections of the thermal impact of the Siberian Traps, which reduced phlogopite-bearing metasomes within the southwestern Siberian sub-continental lithospheric mantle.