Provenance discrimination of Lower Cretaceous synrift sandstones (eastern Iberian Chain Spain) Constraints from detrital modes heavy minerals and geochemistry

Petrography, geochemical whole-rock composition, and chemical analyses of tourmaline were performed in order to determine the source areas of Lower Cretaceous Mora, El Castellar, and uppermost Camarillas Formation sandstones from the Iberian Chain, Spain. Sandstones were deposited in intraplate subbasins, which are bound by plutonic and volcanic rocks of Permian, Triassic, and Jurassic age, Paleozoic metamorphic rocks, and Triassic sedimentary rocks. Modal analyses together with petrographic and cathodoluminescence observations allowed us to define three quartz-feldspathic petrofacies and recognize diagenetic processes that modified the original framework composition. Results from average restored petrofacies are: Mora petrofacies = P/F >1 and Q(r)70 F(r)22 R(r)9; El Castellar petrofacies = P/F >1 and Q(r)57 F(r)25 R(r)18; and Camarillas petrofacies = P/F ∼ zero and Q(r)64 F(r)28 R(r)7 (P—plagioclase; F—feldspar; Q—quartz; R—rock fragments; r—restored composition). Trace-element and rare earth element abundances of whole-rock analyses discriminate well between the three petrofacies based on: (1) the Rb concentration, which is indicative of the K content and reflects the amount of K-feldspar modal abundance, and (2) the relative modal abundance of heavy minerals (tourmaline, zircon, titanite, and apatite), which is reproduced by the elements hosted in the observed heavy mineral assemblage (i.e., B and Li for tourmaline; Zr, Hf, and Ta for zircon; Ti, Ta, Nb, and their rare earth elements for titanite; and P, Y, and their rare earth elements for apatite). Tourmaline chemical composition for the three petrofacies ranges from Fe-tourmaline of granitic to Mg-tourmaline of metamorphic origin. The three defined petrofacies suggest a mixed provenance from plutonic and metamorphic source rocks. However, a progressively major influence of granitic source rocks was detected from the lowermost Mora petrofacies toward the uppermost Camarillas petrofacies. This provenance trend is consistent with the uplift and erosion of the Iberian Massif, which coincided with the development of the latest Berriasian synrift regional unconformity and affected all of the Iberian intraplate basins. The uplifting stage of Iberian Massif pluton caused a significant dilution of Paleozoic metamorphic source areas, which were dominant during the sedimentation of the lowermost Mora and El Castellar petrofacies. The association of petrographic data with whole-rock geochemical compositions and tourmaline chemical analysis has proved to be useful for determining source area characteristics, their predominance, and the evolution of source rock types during the deposition of quartz-feldspathic sandstones in intraplate basins. This approach ensures that provenance interpretation is consistent with the geological context.

Kyzylkumite, Ti2V3+O5(OH): new structure type, modularity and revised formula

The crystal structure of kyzylkumite, ideally Ti2V3+O5(OH), from the Sludyanka complex in South Baikal, Russia was solved and refined (including the hydrogen atom position) to an agreement index, R1, of 2.34 using X-ray diffraction data collected on a twinned crystal. Kyzylkumite crystallizes in space group P21/c, with a = 8.4787(1), b = 4.5624(1), c = 10.0330(1) Å, β = 93.174(1)°, V = 387.51(1) Å3 and Z = 4. Tivanite, TiV3+O3OH, and kyzylkumite have modular structures based on hexagonal close packing of oxygen, which are made up of rutile TiO2 and montroseite V3+O(OH) slices. In tivanite the rutile:montroseite ratio is 1:1, in kyzylkumite the ratio is 2:1. The montroseite module may be replaced by the isotypic paramontroseite V4+O2 module, which produces a phase with the formula Ti2V4+O6. In the metamorphic rocks of the Sludyanka complex, vanadium can be present as V4+ and V3+ within the same mineral (e.g. in batisivite, schreyerite and berdesinskiite). Kyzylkumite has a flexible composition with respect to the M4+/M3+ ratio. The relationship between kyzylkumite and a closely related Be-bearing kyzylkumite-like mineral with an orthorhombic norbergite-type structure from Byrud mine, Norway is discussed. Both minerals have similar X-ray powder diffraction patterns.

Seismic anisotropy in the Morcles nappe shear zone: Implications for seismic imaging of crustal scale shear zones

Microstructures and textures of calcite mylonites from the Morcles nappe large-scale shear zone in southwestern Switzerland develop principally as a function of 1) extrinsic physical parameters including temperature, stress, strain, strain rate and 2) intrinsic parameters, such as mineral composition. We collected rock samples at a single location from this shear zone, on which laboratory ultrasonic velocities, texture and microstructures were investigated and quantified. The samples had different concentration of secondary mineral phases (< 5 up to 40 vol.%). Measured seismic P wave anisotropy ranges from 6.5% for polyphase mylonites (~ 40 vol.%) to 18.4% in mylonites with < 5 vol.% secondary phases. Texture strength of calcite is the main factor governing the seismic P wave anisotropy. Measured S wave splitting is generally highest in the foliation plane, but its origin is more difficult to explain solely by calcite texture. Additional texture measurements were made on calcite mylonites with low concentration of secondary phases (≤ 10 vol.%) along the metamorphic gradient of the shear zone (15 km distance). A systematic increase in texture strength is observed moving from the frontal part of the shear zone (anchimetamorphism; 280 °C) to the higher temperature, basal part (greenschist facies; 350–400 °C). Calculated P wave velocities become increasingly anisotropic towards the high-strain part of the nappe, from an average of 5.8% in the frontal part to 13.2% in the root of the basal part. Secondary phases raise an additional complexity, and may act either to increase or decrease seismic anisotropy of shear zone mylonites. In light of our findings we reinterpret the origin of some seismically reflective layers in the Grône–Zweisimmen line in southwestern Switzerland (PNR20 Swiss National Research Program). We hypothesize that reflections originate in part from the lateral variation in textural and microstructural arrangement of calcite mylonites in shear zones.

Pressure-temperature estimates of the lizardite/antigorite transition in high pressure serpentinites

Serpentine minerals in natural samples are dominated by lizardite and antigorite. In spite of numerous laboratory experiments, the stability fields of these species remain poorly constrained. This paper presents petrological observations and the Raman spectroscopy and XRD analyses of natural serpentinites from the Alpine paleo-accretionary wedge. Serpentine varieties were identified from a range of metamorphic pressure and temperature conditions from sub-greenschist (P < 4 kbar, T ~ 200–300 °C) to eclogite facies conditions (P > 20 kbar, T > 460 °C) along a subduction geothermal gradient. We use the observed mineral assemblage in natural serpentinite along with the Tmax estimated by Raman spectroscopy of the carbonaceous matter in associated metasediments to constrain the temperature of the lizardite to antigorite transition at high pressures. We show that below 300 °C, lizardite and locally chrysotile are the dominant species in the mesh texture. Between 320 and 390 °C, lizardite is progressively replaced by antigorite at the grain boundaries through dissolution–precipitation processes in the presence of SiO2 enriched fluids and in the cores of the lizardite mesh. Above 390 °C, under high-grade blueschist to eclogite facies conditions, antigorite is the sole stable serpentine mineral until the onset of secondary olivine crystallization at 460 °C.

Deciphering High-Pressure metamorphism in collisional context using microprobe-mapping methods : application to the Stak eclogitic massif (NW Himalaya)

The Stak massif, northern Pakistan, is a newly recognized occurrence of eclogite formed by the subduction of the northern margin of the Indian continent in the northwest Himalaya. Although this unit was extensively retrogressed during the Himalayan collision, records of the high-pressure (HP) event as well as a continuous pressure-temperature (P-T) path were assessed from a single thin section using a new multiequilibrium method. This method uses microprobe X-ray compositional maps of garnet and omphacitic pyroxene followed by calculations of ∼200,000 P-T estimates using appropriate thermobarometers. The Stak eclogite underwent prograde metamorphism, increasing from 650 °C and 2.4 GPa to the peak conditions of 750 °C and 2.5 GPa, then retrogressed to 700–650 °C and 1.6–0.9 GPa under amphibolite-facies conditions. The estimated peak metamorphic conditions and P-T path are similar to those of the Kaghan and Tso Morari high- to ultrahigh-pressure (HP-UHP) massifs. We propose that these three massifs define a large HP to UHP province in the northwest Himalaya, comparable to the Dabie-Sulu province in China and the Western Gneiss Region in Norway.

Neotethys closure history of Anatolia: insights from 40Ar–39Ar geochronology and P–T estimation in high pressure metasedimentary rocks

The multiple high-pressure (HP), low-temperature (LT) metamorphic units of Western and Central Anatolia offer a great opportunity to investigate the subduction- and continental accretion-related evolution of the eastern limb of the long-lived Aegean subduction system. Recent reports of the HP–LT index mineral Fe-Mg-carpholite in three metasedimentary units of the Gondwana-derived Anatolide–Tauride continental block (namely the Afyon Zone, the Ören Unit and the southern Menderes Massif) suggest a more complicated scenario than the single-continental accretion model generally put forward in previous studies. This study presents the first isotopic dates (white mica 40Ar–39Ar geochronology), and where possible are combined with P–T estimates (chlorite thermometry, phengite barometry, multi-equilibrium thermobarometry), on carpholite-bearing rocks from these three HP–LT metasedimentary units. It is shown that, in the Afyon Zone, carpholite-bearing assemblages were retrogressed through greenschist-facies conditions at c. 67–62 Ma. Early retrograde stages in the Ören Unit are dated to 63–59 Ma. In the Kurudere–Nebiler Unit (HP Mesozoic cover of the southern Menderes Massif), HP retrograde stages are dated to c. 45 Ma, and post-collisional cooling to c. 26 Ma. These new results support that the Ören Unit represents the westernmost continuation of the Afyon Zone, whereas the Kurudere–Nebiler Unit correlates with the Cycladic Blueschist Unit of the Aegean Domain. In Western Anatolia, three successive HP–LT metamorphic belts thus formed: the northernmost Tavşanlı Zone (c. 88–82 Ma), the Ören–Afyon Zone (between 70 and 65 Ma), and the Kurudere–Nebiler Unit (c. 52–45 Ma). The southward younging trend of the HP–LT metamorphism from the upper and internal to the deeper and more external structural units, as in the Aegean Domain, points to the persistence of subduction in Western Anatolia between 93–90 and c. 35 Ma. After the accretion of the Menderes–Tauride terrane, in Eocene times, subduction stopped, leading to continental collision and associated Barrovian-type metamorphism. Because, by contrast, the Aegean subduction did remain active due to slab roll-back and trench migration, the eastern limb (below Southwestern Anatolia) of the Hellenic slab was dramatically curved and consequently teared. It therefore is suggested that the possibility for subduction to continue after the accretion of buoyant (e.g. continental) terranes probably depends much on palaeogeography.

REE and Hf distribution among mineral phases in the CV-CK clan : a way to explain present-day Hf isotopic variations in chondrites

Chondrites are among the most primitive objects in the Solar System and constitute the main building blocks of telluric planets. Among the radiochronometers currently used for dating geological events, Sm–Nd and Lu–Hf are both composed of refractory, lithophile element. They are thought to behave similarly as the parent elements (Sm and Lu) are generally less incompatible than the daughter elements (Nd and Hf) during geological processes. As such, their respective average isotopic compositions for the solar system should be well defined by the average of chondrites, called Chondritic Uniform Reservoir (CHUR). However, while the Sm–Nd isotopic system shows an actual spread of less than 4% in the average chondritic record, the Lu–Hf system shows a larger variation range of 28% [Bouvier A., Vervoort J. D. and Patchett P. J. (2008) The Lu–Hf and Sm–Nd isotopic composition of CHUR: Constraints from unequilibrated chondrites and implications for the bulk composition of terrestrial planets. Earth Planet. Sci. Lett.273, 48–57]. To better understand the contrast between Sm–Nd and Lu–Hf systems, the REE and Hf distribution among mineral phases during metamorphism of Karoonda (CK) and Vigarano-type (CV) carbonaceous chondrites has been examined. Mineral modes were determined from elemental mapping on a set of five CK chondrites (from types 3–6) and one CV3 chondrite. Trace-element patterns are obtained for the first time in all the chondrite-forming minerals of a given class (CK chondrites) as well as one CV3 sample. This study reveals that REE are distributed among both phosphates and silicates. Only 30–50% of Sm and Nd are stored in phosphates (at least in chondrites types 3–5); as such, they are not mobilized during early stages of metamorphism. The remaining fraction of Sm and Nd is distributed among the same mineral phases; these elements are therefore not decoupled during metamorphism. Of the whole-rock total of Lu, the fraction held in phosphate decreases significantly as the degree of metamorphism increases (30% for types 3 and 4, less than 5% in type 6). In contrast to Lu, Hf is mainly hosted by silicates with little contribution from phosphates throughout the CK metamorphic sequence. A significant part of Sm and Nd are stored in phosphates in types 3–5, and these elements behave similarly during CK chondrite metamorphism. That explains the robustness of the Sm/Nd ratios in chondrites through metamorphism, and the slight discrepancies observed in the present-day isotopic Nd values in chondrites. On the contrary, Lu and Hf are borne by several different minerals and consequently they are redistributed during metamorphism–induced recrystallization. The Lu/Hf ratios are therefore significantly disturbed during chondrites metamorphism, leading to the high discrepancies observed in present-day Hf isotopic values in chondrites.

Calibration of cosmogenic noble gas production in ordinary chondrites based on 36 Cl- 36 Ar ages. Part 1: Refined produced rates for cosmogenic 21 Ne and 38 Ar

We measured the concentrations and isotopic compositions of He, Ne, and Ar in bulk samples and metal separates of 14 ordinary chondrite falls with long exposure ages and high metamorphic grades. In addition, we measured concentrations of the cosmogenic radionuclides 10Be, 26Al, and 36Cl in metal separates and in the nonmagnetic fractions of the selected meteorites. Using cosmogenic 36Cl and 36Ar measured in the metal separates, we determined 36Cl-36Ar cosmic-ray exposure (CRE) ages, which are shielding-independent and therefore particularly reliable. Using the cosmogenic noble gases and radionuclides, we are able to decipher the CRE history for the studied objects. Based on the correlation 3He/21Ne versus 22Ne/21Ne, we demonstrate that, among the meteorites studied, only one suffered significant diffusive losses (about 35%). The data confirm that the linear correlation 3He/21Ne versus 22Ne/21Ne breaks down at high shielding. Using 36Cl-36Ar exposure ages and measured noble gas concentrations, we determine 21Ne and 38Ar production rates as a function of 22Ne/21Ne. The new data agree with recent model calculations for the relationship between 21Ne and 38Ar production rates and the 22Ne/21Ne ratio, which does not always provide unique shielding information. Based on the model calculations, we determine a new correlation line for 21Ne and 38Ar production rates as a function of the shielding indicator 22Ne/21Ne for H, L, and LL chondrites with preatmospheric radii less than about 65 cm. We also calculated the 10Be/21Ne and 26Al/21Ne production rate ratios for the investigated samples, which show good agreement with recent model calculations.

Assessing Initial Conditions for Chloride Transport Across Low- permeability Argillaceous Rocks, Wellenberg, Switzerland

Information about fluid evolution and solute transport in a low-permeability metamorphic rock sequence has been obtained by comparing chloride concentrations and chlorine isotope ratios of pore water, groundwater, and fluid inclusions. The similarity of d37Cl values in fluid inclusions and groundwater suggests a closed-system evolution during the metamorphic overprint, and signatures established at this time appear to form the initial conditions for chloride transport after exhumation of the rock sequence.

Small quantity but large effect — How minor phases control strain localization in upper mantle shear zones

Low viscosity domains such as localized shear zones exert an important control on the geodynamics of the uppermost mantle. Grain size reduction and subsequent strain localization related to a switch from dislocation to diffusion creep is one mechanism to form low viscosity domains. To sustain strain localization, the grain size of mantle minerals needs to be kept small over geological timescales. One way to keep olivine grain sizes small is by pinning of mobile grain boundaries during grain growth by other minerals (second phases). Detailed microstructural studies based on natural samples from three shear zones formed at different geodynamic settings, allowed the derivation of the olivine grain-size dependence on the second-phase content. The polymineralic olivine grain-size evolution with increasing strain is similar in the three shear zones. If the second phases are to pin the mobile olivine grain boundary the phases need to be well mixed before grain growth. We suggest that melt-rock and metamorphic reactions are crucial for the initial phase mixing in mantle rocks. With ongoing deformation and increasing strain, grain boundary sliding combined with mass transfer processes and nucleation of grains promotes phase mixing resulting in fine-grained polymineralic mixtures that deform by diffusion creep. Strain localization due to the presence of volumetrically minor minerals in polymineralic mantle rocks is only important at high strain deformation (ultramylonites) at low temperatures (<~800°C). At smaller strain and stress conditions and/or higher temperatures other parameters like overall energy available to deform a given rock volume, the inheritance of mechanical anisotropies or the presence of water or melts needs to be considered to explain strain localization in the upper mantle.

Magnetic anisotropy in clinopyroxene and orthopyroxene single crystals

Pyroxenes constitute an important component in mafic igneous and metamorphic rocks. They often possess a prismatic habit, and their long axis, the crystallographic c axis, helps define a lineation in a textured rock. Anisotropy of magnetic susceptibility (AMS) serves as a fabric indicator in igneous and metamorphic rocks. If a rock’s AMS is carried by pyroxenes, it can be related to their crystallographic preferred orientation and degree of alignment. This requires knowing the intrinsic AMS of pyroxene single crystals. This study provides a comprehensive low-field and high-field AMS investigation of chemically diverse orthopyroxene and clinopyroxene crystals in relation to crystal structure, chemical composition, oxidation state of Fe, and the possible presence of ferromagnetic inclusions. The paramagnetic anisotropy, extracted from high-field data, shows clear relationships to crystallographic directions and Fe concentration both in clinopyroxene and orthopyroxene. In the diopside-augite series, the intermediate susceptibility is parallel to b, and the maximum is at 45° to the c axis. In aegirine, the intermediate axis remains parallel to b, while the maximum susceptibility is parallel to c. The AMS of spodumene depends on Fe concentration. In enstatite, the maximum susceptibility aligns with c and the minimum with b, and in the case of hypersthene, the maximum susceptibility is normal to the exsolution lamellae. Magnetite inclusions within augite possess a ferromagnetic anisotropy with consistent orientation of the principal susceptibilities, which dominates the low-field anisotropy. These results provide better understanding of magnetic anisotropy in pyroxenes and form a solid basis for interpretation of magnetic fabrics in pyroxene-bearing rocks.

Age of cleft monazites in the eastern Tauern Window: constraints on crystallization conditions of hydrothermal monazite

Monazite-bearing Alpine clefts located in the Sonnblick region of the eastern Tauern Window, Austria, are oriented perpendicular to the foliation and lineation. Ion probe (SIMS) Th–Pb and U–Pb dating of four cleft monazites yields crystallization ages of different growth domains and aggregate regions ranging from 18.99 ± 0.51 to 15.00 ± 0.51 Ma. The crystallization ages obtained are overlapping or slightly younger than zircon fission track ages but older than zircon (U–Th)/He cooling ages from the same area. This constrains cleft monazite crystallization in this area to *300–200 �C. LA-ICP-MS data of dated hydrothermal monazites indicate that in graphite-bearing, reduced host lithologies, cleft monazite is poor in As and has higher La/Yb values and U concentrations, whereas in oxidised host rocks opposite trends are observed. Monazites show negative Eu anomalies and variable La/Yb values ranging from 520 to 6050. The positive correlation between Ca and Sr concentration indicates dissolution of plagioclase or carbonates as the source of these elements. The data show that early exhumation and cleft formation in the Tauern is related to metamorphic dome formation caused by the collision of the Adriatic with the European plate and that monazite crystallization in the clefts occurred later. Our data also demonstrate that hydrothermal monazite ages offer great potential in helping to constrain the chronology of exhumation in collisional orogens.

XMapTools: a MATLAB©-based program for electron microprobe X-ray image processing and geothermobarometry

XMapTools is a MATLAB©-based graphical user interface program for electron microprobe X-ray image processing, which can be used to estimate the pressure–temperature conditions of crystallization of minerals in metamorphic rocks. This program (available online at http://www.xmaptools.com) provides a method to standardize raw electron microprobe data and includes functions to calculate the oxide weight percent compositions for various minerals. A set of external functions is provided to calculate structural formulae from the standardized analyses as well as to estimate pressure–temperature conditions of crystallization, using empirical and semi-empirical thermobarometers from the literature. Two graphical user interface modules, Chem2D and Triplot3D, are used to plot mineral compositions into binary and ternary diagrams. As an example, the software is used to study a high-pressure Himalayan eclogite sample from the Stak massif in Pakistan. The high-pressure paragenesis consisting of omphacite and garnet has been retrogressed to a symplectitic assemblage of amphibole, plagioclase and clinopyroxene. Mineral compositions corresponding to ~165,000 analyses yield estimates for the eclogitic pressure–temperature retrograde path from 25 kbar to 9 kbar. Corresponding pressure–temperature maps were plotted and used to interpret the link between the equilibrium conditions of crystallization and the symplectitic microstructures. This example illustrates the usefulness of XMapTools for studying variations of the chemical composition of minerals and for retrieving information on metamorphic conditions on a microscale, towards computation of continuous pressure–temperature-and relative time path in zoned metamorphic minerals not affected by post-crystallization diffusion.