Pressure, temperature and time constraints on tectonic models for southwestern Sweden

In this work, a number of key localities have been investigated in detail in order to provide precise constraints on models for the tectonic evolution of southern Sweden. The new data presented in this thesis show that there are large differences in terms of pre-Sveconorwegian tectonic evolution between the Eastern and Median-Western Segments situated on either side of the Mylonite Zone, a major shear zone as well as a structural and lithological boundary. This has a direct influence on the possible tectonic scenarios that can be suggested when reconstructing the formation of the south-western part of the Fennoscandian Shield. At Viared in the central Eastern Segment, Sveconorwegian eclogite facies metamorphism is dated at 0.97 Ga using mainly U-Pb on zircon. This is similar to other localities showing high-pressure granulite or eclogite facies metamorphism in the Eastern Segment and suggests that this high grade event was a regional feature east of the Mylonite Zone. On the well exposed Nordön Island in the Western Segment, both pre-Sveconorwegian and Sveconorwegian metamorphism and deformation was dated using several isotope methods, including U-Pb SIMS zircon and Sm-Nd garnet dating. In the Median Segment, veining was dated using U-Pb SIMS zircon methods. In addition to age determinations, thermobarometry was done on several samples and the results compared with published data. The results show that Sveconorwegian peak metamorphism reached amphibolite to granulite facies west of the Mylonite Zone at 1.02-1.04 Ga. This is 50-70 Ma before the orogenic activity started in the Eastern Segment. In the central Eastern Segment, Pre-Sveconorwegian veining is dated at 1.42 Ga, thus belonging to the 1.42-1.46 Ga Hallandian veining found elsewhere in the Eastern Segment. Pre-Sveconorwegian veining and isoclinal folding in the Western Segment is dated at 1.55 Ga. Neither the 1.55 Ga nor the 1.02-1.04 Ga events have been found in the Eastern Segment. Further, the 1.42-1.46 Ga veining event documented in the Eastern Segment has not been found in the Western or Median Segments. Only a few 0.97 Ga zircon U-Pb ages have been found in the Western and Median Segments west of the Mylonite Zone. In the Western Segment, Ar-Ar dating of muscovite related to late Sveconorwegian uplift gives 981±4 Ma, interpreted to reflect the crystallization of muscovite below 400 °C. This indicates that the Western Segment already was exhumed at the time when the Eastern Segment experienced its Sveconorwegian high-pressure peak metamorphism. Additionally this work shows that there is an age difference between the 1.34 Ga Askim, the 1.30 Ga Göta and 1.31 Ga Kärra granites west of the Mylonite Zone and the 1.38-1.40 Ga granites and monzonites east of the Mylonite Zone, thus disproving the previous concept of “stitching granites” that was used as an argument for a pre-Sveconorwegian correlation between the Eastern and Western Segments. The differences in pre-Sveconorwegian history between the segments east and west of the Mylonite zone suggest that the zone originally formed as a Sveconorwegian first order crustal suture. However, no ophiolites have been found along the zone and no calk-alkaline magmatism related to subduction of oceanic crust slightly before 0.97 Ma has yet been documented. If such features were found it would strengthen the idea that the Mylonite Zone is a crustal suture, however their lack does not preclude it. Geothermobarometry on retrograde eclogite facies rocks exposed at Viared indicates that those rocks experienced metamorphic conditions of 15.0–16.7 kbar at temperatures of 719 to 811°C. The equivalent burial depth of ~50 km is compatible with metamorphism in a subduction setting. The partial preservation of a high pressure paragenesis indicates rapid exhumation after burial. A two-dimensional model involving continental subduction of the Eastern Segment is proposed to explain the cycle of eclogite and high-pressure granulite facies that took place during the ~0.97 Ga Falkenberg phase of the Sveconorwegian Orogeny

Oxygen-isotope and trace element constraints on the origins of silica-rich melts in the subarc mantle

Peridotitic xenoliths in basaltic andesites from Batan island in the Luzon arc contain silica-rich (broadly dacitic) hydrous melt inclusions that were likely trapped when these rocks were within the upper mantle wedge underlying the arc. These melt inclusions have been previously interpreted to be slab-derived melts. We tested this hypothesis by analyzing the oxygen isotope compositions of these inclusions with an ion microprobe. The melt inclusions from Batan xenoliths have δ 18OVSMOW values of 6.45 ± 0.51‰. These values are consistent with the melts having been in oxygen isotope exchange equilibrium with average mantle peridotite at temperatures of ≥875°C. We suggest the δ 18O values of Batan inclusions, as well as their major and trace element compositions, can be explained if they are low-degree melts (or differentiation products of such melts) of peridotites in the mantle wedge that had previously undergone extensive melt extraction followed by metasomatism by small amounts (several percent or less) of slab-derived components. A model based on the trace element contents of Batan inclusions suggests that this metasomatic agent was an aqueous fluid extracted from subducted basalts and had many characteristics similar to slab-derived components of the sources of arc-related basalts at Batan and elsewhere. Batan inclusions bear similarities to “adakites,” a class of arc-related lava widely considered to be slab-derived melts. Our results suggest the alternative interpretation that at least some adakite-like liquids might be generated from low-degree melting of metasomatized peridotites

Sulfur loss from subducted altered oceanic crust and implications for mantle oxidation

© The Author(s), [year]. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Walters, J. B., Cruz-Uribe, A. M., & Marschall, H. R. Sulfur loss from subducted altered oceanic crust and implications for mantle oxidation. Geochemical Perspectives Letters, 13, (2020): 36-41, doi:10.7185/geochemlet.2011.Oxygen fugacity (fO2) is a controlling factor of the physics of Earth’s mantle; however, the mechanisms driving spatial and secular changes in fO2 associated with convergent margins are highly debated. We present new thermodynamic models and petrographic observations to predict that oxidised sulfur species are produced during the subduction of altered oceanic crust. Sulfur loss from the subducting slab is a function of the protolith Fe3+/ΣFe ratio and subduction zone thermal structure, with elevated sulfur fluxes predicted for oxidised slabs in cold subduction zones. We also predict bi-modal release of sulfur-bearing fluids, with a low volume shallow flux of reduced sulfur followed by an enhanced deep flux of sulfate and sulfite species, consistent with oxidised arc magmas and associated copper porphyry deposits. The variable SOx release predicted by our models both across and among active margins may introduce fO2 heterogeneity to the upper mantle.We thank James Connolly for modelling support and Peter van Keken for providing updated P–T paths for the Syracuse et al. (2010) models. The manuscript benefited from the editorial handling by Helen Williams and from constructive reviews of Maryjo Brounce, Katy Evans, and an anonymous reviewer. JBW acknowledges Fulbright and Chase Distinguished Research fellowships. This work was supported by NSF grant EAR1725301 awarded to AMC

Recycling Argon through Metamorphic Reactions: the Record in Symplectites

The 40Ar/39Ar ages of metamorphic micas that crystallized at high temperatures are commonly interpreted as cooling ages, with grains considered to have lost 40Ar via thermally-driven diffusion into the grain boundary network. Recently reported laser-ablation data suggest that the spatial distribution of Ar in metamorphic micas does not always conform to the patterns predicted by diffusion theory and that despite high metamorphic temperatures, argon was not removed efficiently from the local system during metamorphic evolution. In the Western Gneiss Region (WGR), Norway, felsic gneisses preserve microtextural evidence for the breakdown of phengite to biotite and plagioclase symplectites during near isothermal decompression from c. 20–25 to c. 8–12 kbar at ~700°C. These samples provide an ideal natural laboratory to assess whether the complete replacement of one K-bearing mineral by another at high temperatures completely ‘resets’ the Ar clock, or whether there is some inheritance of 40Ar in the neo-crystallized phase. The timing of the high-temperature portion of the WGR metamorphic cycle has been well constrained in previous studies. However, the timing of cooling following the overprint is still much debated. In-situ laser ablation spot dating in phengite, biotite-plagioclase symplectites and coarser, texturally later biotite yielded 40Ar/39Ar ages that span much of the metamorphic cycle. Together these data show that despite residence at temperatures of ~700°C, Ar is not completely removed by diffusive loss or during metamorphic recrystallization. Instead, Ar released during phengite breakdown appears to be partially reincorporated into the newly crystallizing biotite and plagioclase (or is trapped in fluid inclusions in those phases) within a close system. Our data show that the microtextural and petrographic evolution of the sample being dated provides a critical framework in which local 40Ar recycling can be tracked, thus potentially allowing 40Ar/39Ar dates to be linked more accurately to metamorphic history

Deeply subducted continental fragments - Part 2: Insight from petrochronology in the central Sesia Zone (western Italian Alps)

Subducted continental terranes commonly comprise an assembly of subunits that reflect the different tectono-metamorphic histories they experienced in the subduction zone. Our challenge is to unravel how, when, and in which part of the subduction zone these subunits were juxtaposed. Petrochronology offers powerful tools to decipher pressure–temperature–time (P–T–t) histories of metamorphic rocks that preserve a record of several stages of transformation. A major issue is that the driving forces for re-equilibration at high pressure are not well understood. For example, continental granulite terrains subducted to mantle depths frequently show only partial and localized eclogitization. The Sesia Zone (NW Italy) is exceptional because it comprises several continental subunits in which eclogitic rocks predominate and high-pressure (HP) assemblages almost completely replaced the Permian granulite protoliths. This field-based study comprises both main complexes of the Sesia terrane, covering some of the recently recognized tectonic subunits involved in its assembly; hence our data constrain the HP tectonics that formed the Sesia Zone. We used a petrochronological approach consisting of petrographic and microstructural analysis linked with thermodynamic modelling and U–Th–Pb age dating to reconstruct the P–T–t trajectories of these tectonic subunits. Our study documents when and under what conditions re-equilibration took place. Results constrain the main stages of mineral growth and deformation, associated with fluid influx that occurred in the subduction channel. In the Internal Complex (IC), pulses of fluid percolated at eclogite facies conditions between 77 and 55 Ma with the HP conditions reaching  ∼  2 GPa and 600–670 °C. By contrast, the External Complex (EC) records a lower pressure peak of  ∼  0.8 GPa for 500 °C at  ∼  63 Ma. The juxtaposition of the two complexes occurred during exhumation, probably at  ∼  0.8 GPa and 350 °C; the timing is constrained between 46 and 38 Ma. Mean vertical exhumation velocities are constrained between 0.9 and 5.1 mm year−1 for the IC, up to its juxtaposition with the EC. Exhumation to the surface occurred before 32 Ma, as constrained by the overlying Biella Volcanic Suite, at a mean vertical velocity between 1.6 and 4 mm year−1. These findings constrain the processes responsible for the assembly and exhumation of HP continental subunits, thus adding to our understanding of how continental terranes behave during subduction

Evidence for a reducing Archean ambient mantle and its effects on the carbon cycle

Chemical reduction-oxidation mechanisms within mantle rocks link to the terrestrial carbon cycle by influencing the depth at which magmas can form, their composition, and ultimately the chemistry of gases released into the atmosphere. The oxidation state of the uppermost mantle has been widely accepted to be unchanged over the past 3800 m.y., based on the abundance of redox-sensitive elements in greenstone belt–associated samples of different ages. However, the redox signal in those rocks may have been obscured by their complex origins and emplacement on continental margins. In contrast, the source and processes occurring during decompression melting at spreading ridges are relatively well constrained. We retrieve primary redox conditions from metamorphosed mid-oceanic ridge basalts (MORBs) and picrites of various ages (ca. 3000–550 Ma), using V/Sc as a broad redox proxy. Average V/Sc values for Proterozoic suites (7.0 ± 1.4, 2s, n = 6) are similar to those of modern MORB (6.8 ± 1.6), whereas Archean suites have lower V/Sc (5.2 ± 0.4, n = 5). The lower Archean V/Sc is interpreted to reflect both deeper melt extraction from the uppermost mantle, which becomes more reduced with depth, and an intrinsically lower redox state. The pressure-corrected oxygen fugacity (expressed relative to the fayalite-magnetite-quartz buffer, DFMQ, at 1 GPa) of Archean sample suites (DFMQ –1.19 ± 0.33, 2s) is significantly lower than that of post-Archean sample suites, including MORB (DFMQ –0.26 ± 0.44). Our results imply that the reducing Archean atmosphere was in equilibrium with Earth’s mantle, and further suggest that magmatic gases crossed the threshold that allowed a build-up in atmospheric O2 levels ca. 3000 Ma, accompanied by the first “whiffs” of oxygen in sediments of that age

The gabbro-eclogite phase transition and the elevation of mountain belts on Venus

Among the four mountain belts surrounding Lakshmi Planum, Maxwell Montes is the highest and stands up to 11 km above the mean planetary radius and 7 km above Lakshmi Planum. The bulk composition and radioactive heat production of the crust on Venus, where measured, are similar to those of terrestrial tholeiitic basalt. Because the thickness of the low-density crust may be limited by the gabbro-garnet granulite-eclogite phase transitions, the 7-11 km maximum elevation of Maxwell Montes is difficult to understand except in the unlikely situation that the crust contains a large volume of magma. A possible explanation is that the base of the crust is not in phase equilibrium. It has been suggested that under completely dry conditions, the gabbro-eclogite phase transition takes place by solid-state diffusion and may require a geologically significant time to run to completion. Solid-state diffusion is a strongly temperature-dependent process. In this paper we solve the thermal evolution of the mountain belt to attempt to constrain the depth of the gabbro-eclogite transition and thus to assess this hypothesis quantitatively. The one-dimensional heat equation is solved numerically by a finite difference approximation. The deformation of the horizontally shortening crustal and mantle portions of the thermal boundary layer is assumed to occur by pure shear, and therefore the vertical velocity is given by the product of the horizontal strain rate and depth

The effect of eclogitization of crustal rocks on the seismic properties on variable scales: Implications for geophysical imaging of eclogitization at depth

Plate tectonics shapes the face of the earth and subduction and collision zones are among the most important features on Earth. Here, crustal material is recycled into the mantle or integrated into growing orogens. However, the processes active at depth cannot be studied directly and we thus rely on geophysical imaging methods to visualize the geometries that result from the ongoing processes. Additionally, these processes can be studied in fossil subduction and collision zones. However, the scales at which observations from geophysical imaging are made are orders of magnitude larger than those made in field-based studies of fossil subduction and collision zones. This thesis provides insight into how eclogitization modifies the physical properties of deeply buried rocks and what influence the resulting lithologies and their geometrical configuration have on geophysical imaging. In an interdisciplinary approach, I show how structures that are likely representative for those present at depth in subduction and collision zones develop and what their geometries at depth will be. I then derive their petrophysical properties and show how these are modified on various scales, and how this influences the detectability of such associations using geophysical imaging techniques. To do so, the island of Holsnøy in western Norway serves as a natural laboratory that is ideal to study eclogitization of crustal material. Geological mapping on Holsnøy constrains the geometric framework of the constituting lithologies and the scales at which such structures could be expected to establish. Previously, several authors have shown that many of the eclogite occurrences on Holsnøy are produced contemporaneously with ductile deformation forming shear zones at various scales. Our geological mapping aided by photogrammetry using drone images reveals that large parts of this exposed continental sliver were eclogitized statically without associated ductile deformation. This shows that even in domains with ongoing regional deformation, low-strain domains develop within the descending crustal material. Nevertheless, even the major shear zones that are exposed are only a few hundred meters thick, and thus far below the scale that is detectable by geophysical imaging techniques. However, geological mapping of the area suggests that the exposed structures are, at least in a qualitative sense, scale independent, suggesting that the same structural framework could be present at a larger scale in active subduction and collision zones. Measurements of P and S wave velocities of the exposed granulitic protolith and eclogites suggest that eclogitization of the lower crust causes three major changes of the petrophysical properties: (1) increased P and S wave velocities, (2) an increase of the seismic anisotropy, and (3) a decrease of the VP/VS ratio, suggesting distinct variations in the geophysical signal when the descending material is partially eclogitized. Additionally, testing the signal that the exposed shear zones would give in reflection seismic and receiver function studies reveals that the variations in shear zone structure indeed produces variations in the retrieved waveforms. Nevertheless, as the exposed structures are too small for geophysical imaging, the finite element method is used to calculate the effective properties of representative structures acting as an effective medium. The results show that the geometrical configuration of the constituting lithologies only has a minor impact on the P wave velocities and anisotropies of the resulting effective medium. Furthermore, our effective medium calculations on the kilometer scale show that eclogitization of crustal material can indeed produce significant seismic anisotropy. In this case, the calculated anisotropy reaches ~5%, which would produce a dependence of the retrieved signal in, for example, receiver function studies on the backazimuth of the sampled rays. Such backazimuthal dependence is indeed observed in active collision zones such as the Himalaya-Tibet collision system and the results presented here can thus be used to constrain the lithologies at depth, suggesting that the lower crust of India below the Himalaya is partially eclogitized along shear zones similar to those exposed on Holsnøy

High-MgO lavas associated to CFB as indicators of plume-related thermochemical effects: the case of ultra-titaniferous picrite-basalt from the Northern Ethiopian-Yemeni plateau

A comprehensive petrological and geochemical dataset is reported in order to define the thermo-compositional characteristics of Ti (Fe)-enriched picrite-basalt lavas (HT2, TiO2 3-7 wt%), erupted close to the axial zone of the inferred Afar mantle plume, at the centre of the originally continuous Ethiopian-Yemeni CFB plateau (ca. 30Ma) which is zonally arranged with progressively lower Ti basalts (HT1, TiO2 2-4 wt%; LT, TiO2 1-3 wt%) toward the periphery. Integrated petrogenetic modelling based on major and trace element analyses of bulk rocks, minerals and melt inclusions in olivines, as well as Sr-Nd-Pb-He-O isotope compositional variations enables us to make several conclusions. 1) The phase equilibria constraints indicate that HT2 primary picrites were generated at ca. 1570°C mantle potential temperatures (Tp) in the pressure range 4-5 GPa whereas the HT1 and LT primary melts formed at shallower level (< 2 to 3 GPa, Tp 1530 °C for HT1 and 1430°C for LT). Thus the Afar plume head was a thermally and compositionally zoned melting region with maximum excess temperatures of 300-350°C with respect to the ambient mantle. 2) The HT2 primary melts upwelled nearly adiabatically to the base of the continental crust (ca. 1 GPa) where fractionation of olivine, followed by clinopyroxene, led to variably differentiated picritic and basaltic magmas. 3) Trace element modelling requires that the primary HT2 melts were generated - either by fractional or batch melting (F 9-10%) - from a mixed garnet peridotite source (85%) with 15% eclogite (derived from transitional MORB protoliths included in Panafrican terranes) that has to be considered a specific Ti-Fe and incompatible element enriched component entrained by the Afar plume. 4) The LT, HT1 and HT2 lavas have 143Nd/144Nd = 0.5131-0.5128, whereas Sr-Pb isotopes are positively correlated with TiO2, varying from 87 Sr/86Sr 0.7032 and 206Pb/204Pb 18.2 in LT basalts to 87Sr/86Sr 0.7044 and 206Pb/204Pb 19.4 in HT2 picrite-basalts. High 3He/4He (15-20 RA) ratios are exclusively observed in HT2 lavas, confirming earlier evidence that these magmas require a component of deep mantle in addition to eclogite, while the LT basalts may more effectively reflect the signature of the pre-existing mantle domains. The comparison between high-MgO (13-22%) lavas from several Phanerozoic CFB provinces (Karoo, Paranà-Etendeka, Emeishan, Siberia, Deccan, North Atlantic Province) shows that they share extremely high mantle potential temperatures (Tp 1550-1700°C) supporting the view that hot mantle plumes are favoured candidates for triggering many LIPs. However, the high incompatible element and isotopic variability of these high-MgO lavas (and associated CFB) suggest that plume thermal anomalies are not necessarily accompanied by significant and specific chemical effects, which depend on the nature of mantle materials recycled during the plume rise, as well as by the extent of related mantle enrichments (if any) on the pre-existing lithospheric section

Fundamental studies in geodynamics

Research in fundamental studies in geodynamics continued in a number of fields including seismic observations and analysis, synthesis of geochemical data, theoretical investigation of geoid anomalies, extensive numerical experiments in a number of geodynamical contexts, and a new field seismic volcanology. Summaries of work in progress or completed during this report period are given. Abstracts of publications submitted from work in progress during this report period are attached as an appendix