Ivrea-Verbano Zone Research Articles

Heavy halogen compositions of peridotite massifs in the Ivrea-Verbano Zone and implications for strong modification of mantle rocks

The lithospheric mantle affected by metasomatism is an important reservoir for halogens and often preserves significant heterogeneity in halogen abundances. However, the halogen compositions of the non-metasomatized mantle rocks with depleted geochemical features have not been well studied, impeding a comprehensive understanding of the halogen distribution in the mantle. The Balmuccia (BM) and Baldissero (BD) peridotite massifs in the Ivrea-Verbano Zone, Italian Alps, represent fragments of the sub-continental lithospheric mantle (SCLM) that are fresh and devoid of subduction-related mantle metasomatism. Here we examine the heavy halogen contents of well-characterized spinel-facies mantle peridotites (n = 20) and pyroxenites (n = 16) from the BM and BD massifs. The results indicate significant variation in whole rock halogen contents of peridotites and pyroxenites, with chlorine (Cl) ranging from <30 to 244 μg/g, bromine (Br) from <0.10 to 0.94 μg/g, and iodine (I) from <0.015 to 0.12 μg/g. Some peridotites and pyroxenites show rather low halogen contents (below the detection limit), consistent with their depleted features. However, a large number of samples display halogen contents far higher than the depleted MORB mantle (i.e., Cl: 3–10 μg/g) and peridotite xenoliths from the strongly metasomatized mantle, although ∼30 % of the halogens were leached out by ultrapure water or dilute nitric acid. The data suggest variable halogen enrichment in many mantle rocks. The halogen contents of different types of mantle rocks do not vary with major elements and incompatible trace elements (e.g., Nb), indicating that halogen compositions are not controlled by magmatic processes such as melting and melt-peridotite reaction. Instead, secondary fluid inclusions are abundant and mainly distributed as trails within olivine and pyroxene in both peridotites and pyroxenites, which were likely trapped during the exhumation of these mantle rocks. Their elevated halogen/Nb ratios compared with MORB mantle and strong correlations with Ba/Nb support the effect of fluids. Micro-Raman analyses identify hydrous serpentines within fluid inclusions. Moreover, mass balance calculation indicates that halogen enrichment is mainly caused by secondary fluids (> 90 %). Together, our findings highlight that the fluid percolation in the mantle rocks during exhumation can intensely modify the halogens, to a similar degree as mantle metasomatism, suggesting the sensitivity of halogens to secondary modification.

Active seismic surveys for drilling target characterization in Ossola Valley: International Continental Scientific Drilling Program (ICDP) project Drilling the Ivrea–Verbano zonE (DIVE) phase I

Abstract. Drilling target locations of the International Continental Scientific Drilling Program (ICDP) project Drilling the Ivrea–Verbano zonE (DIVE) have been initially proposed based on geological knowledge of surface outcrops and the structural context of the Ivrea–Verbano zone (IVZ) and of the Insubric Line. For the determination of the exact locations of drilling sites as well as for drilling geometry planning, we have carried out a series of active seismic experiments to image the subsurface at high resolution. The two drilling sites of project DIVE in Ossola Valley, one near Ornavasso and the other in Megolo di Mezzo, in the central part of the Ivrea–Verbano zone have been surveyed with site-specific velocity models and a seismic data processing chain. The findings have been interpreted in relation with the outcropping structures. These suggest a reasonable continuity from the surface. They also guide the planned borehole orientations: near-vertical at DT-1B (Ornavasso) into the tightly folded Massone Antiform and at 15–20° from the vertical in Megolo across a flank of the broad Proman Anticline. The seismic surveys indicate that the sedimentary overburden is up to 50 m deep at the specific drill sites and can be minimized by relocating the proposed locations. The seismic surveys also indicate that the center of the Ossola Valley contains about 550 m of sedimentary infill, defining the interface of bedrock and Quaternary glacial sediments at about 300 m below sea level.

Composite Melt–Rock Interactions in the Lowermost Continental Crust: Insights from a Dunite–Pyroxenite–Gabbronorite Association of the Mafic Complex from the Ivrea–Verbano Zone (Italian Alps)

Abstract The processes leading to the building of the continental crust through magmatic underplating are fundamentally unknown, mainly because of the rare accessibility to deep level sections of the continental crust. The Italian Alps expose the Permian Mafic Complex, an 8-km-thick gabbronorite–diorite batholith that intruded the lower continental crust during the post-Variscan transtensional tectonics. We present here a petrological and geochemical study of a concentric dunite–pyroxenite–gabbronorite association, called Monte Mazzucco sequence, enclosed at deep levels of the Mafic Complex, thereby allowing us to provide new insights into the magmatic processes driven by emplacement of mantle melts in deep crustal continental areas. The studied sequence includes a ~ 60-m-thick dunite lens, in which olivine (82 mol % forsterite) is associated with accessory Cr-spinel including blebs and lamellae made up of magnetite. The dunite lens is permeated by mm- to cm-scale thick magmatic veins, which range in composition from hornblende lherzolite to olivine hornblendite and hornblende websterite. The lens is mantled by a m-scale ring consisting of amphibole-bearing (≤1 vol %) websterite, and the websterite ring is in turn enclosed by amphibole-free gabbronorites. Both magmatic veins within the dunites and mantling websterites typically include an oxide association of Al-spinel and magnetite. Remarkably, the hornblende websterite veins and the mantling websterites are typically plagioclase-free and include clinopyroxene and amphibole with chondrite-normalized rare earth element patterns characterized by negative Eu anomaly. The mantling websterites display a subtle, gradual outward decrease of Mg# for orthopyroxene, clinopyroxene and accessory olivine, coupled with an increase of the negative Eu anomaly in clinopyroxene and amphibole. The enclosing gabbronorites are amphibole-free and have a chemically evolved signature. We propose a petrogenetic scenario including two major events of melt–dunite interaction. The first resulted from focused reactive melt infiltration and formed the magmatic veins within dunites. The hornblende lherzolite and the olivine hornblendite veins were produced by focused reactive melt migration through dunite grain boundaries, involving dissolution of olivine and recrystallization of Cr-spinel into Al-spinel and magnetite, whereas the hornblende websterite veins crystallized from melts penetrating through narrow fractures and recording earlier plagioclase fractionation. Most likely, the infiltrating melts were overall derived from an evolving H2O-rich magma emplaced below the dunite body. The second event of melt–dunite reactive interaction developed the websterite ring around dunites. We envision that the outermost domain of the dunite body was replaced by websterites in response to reaction with an invading H2O-poor melt, which had previously undergone plagioclase fractionation. The dunite replacement occurred under dynamic conditions, which promoted the reaction progress, thereby leading to total or almost total dissolution of precursor olivine, and started to form the lens shape of the Monte Mazzucco ultramafic association. The gabbronorites closely adjacent to the websterite ring represent the crystallization products of the invading melt.

The Age of Granulite-Facies Metamorphism in the Ivrea–Verbano Zone (NW Italy) Determined Through In Situ U–Pb Dating of Garnet

Abstract Rates of melt production, extraction and crystallization, as well as scales of melt transfer and interaction with their residuum change continuously in migmatite and granulite, affecting the behavior of monazite and zircon as time capsules. Therefore, accessory mineral chronometers may be ambiguous and incomplete in providing an overview of the temperature–time evolution of high-grade metamorphic rocks. In this study, we applied the novel technique of in situ U–Pb dating of garnet to the archetypal lower continental crust of the Ivrea–Verbano Zone (IVZ), NW Italy. In the IVZ, the temporal relationship between granulite-facies metamorphism and mafic underplating has long been debated, because of the interplay between tectonic, magmatic, metamorphic and metasomatic processes over a period of more than a hundred million years. Garnet from mafic and pelitic granulites yielded U–Pb ages between 287.4 ± 4.9 Ma and 280.1 ± 12.4 Ma, overlapping within uncertainty the time proposed for the emplacement of the Mafic Complex (286–282 Ma). These results indicate that the thermal climax in granulitic rocks was caused by mafic underplating and concomitant asthenospheric upwelling, rather than being inherited from the post-Variscan Carboniferous evolution. Providing robust dating of garnet with as low as 4 ng/g U, this study demonstrates the strength of garnet petrochronology in resolving complex tectono-metamorphic histories of high-grade terranes. It also represents a further step forward towards establishing garnet as part of the in situ U–Pb geochronology repertoire.

Seismic structure of the Balmuccia Peridotite from a high-resolution refraction and reflection survey

SUMMARY In anticipation of a forthcoming scientific deep drilling initiative within the Western Alps near Balmuccia, Italy, a high-resolution seismic survey is performed at the proposed drill site. This site is situated within the Ivrea Verbano Zone (IVZ), characterized by lower crustal materials and fragments of upper-mantle rocks exposed adjacent to the Insubric Line. The 2-km-long seismic survey crosses an isolated km-scale outcrop of peridotite near the town of Balmuccia. Applying P-wave traveltime tomography, a substantial contrast in seismic velocities is identified, with velocities in the range of 1–8 km s−1. The peridotite displays velocities ranging from 6 to 8 km s−1. The higher velocities near 8 km s−1 are consistent with laboratory measurements on small-scale samples, while the low-velocity areas within the peridotite body reflect the influence of fractures and faults. The mean velocity derived for the peridotite body is ca. 7 km s−1. The reflection seismic analysis reveals subvertical reflectors positioned at the peridotite boundaries mapped at the surface, converging at a depth of ca. 0.175 km b.s.l. which images a lens-like structure for the peridotite body. However, the area beneath the imaged lens and the deeper Ivrea Geophysical Body (IGB) suggested by earlier studies is not well imaged, which leaves room for other interpretations regarding the relationship of these two bodies. Prior geophysical investigations provide only approximate depth estimates for the top of the IGB, spanning between 1–3 km depth b.s.l. Although the reflection data do not exhibit a series of continuous reflectors beneath the peridotite, a prominent reflection at ca. 1.3 km depth may indicate the top of the IGB.

Effectiveness of Ti-in-amphibole thermometry and performance of different thermometers across lower continental crust up to UHT metamorphism

Metabasites are important constituents of deep crustal sections and are the favored rock type for studying lower crustal amphibolite to granulite transitions. However, metapelites may develop a larger number of temperature-sensitive mineral assemblages and are particular useful when extreme ultrahigh temperature (UHT) conditions are envisaged. A recent calibration of the Ti-in-amphibole thermometer by Liao et al. (2021) was supposed to make thermometry on metabasites quick and easy to apply. However, their calibration is based on experiments which were not originally designed to investigate in detail the temperature dependence of Ti in amphibole. In addition, a possible effect of aTiO2 and/or pressure on the Ti content of amphibole was not fully taken into account. This resulted in a calibration uncertainty of ± 70 °C (2σ), much higher than that of other single-mineral thermometers. In this study we firstly test the newly calibrated Ti-in-amphibole thermometer across the mid to lower crustal section of the Ivrea–Verbano Zone (IVZ; NW Italy) and compare the performance of different thermometric techniques across the sequence. Ti-in-amphibole thermometry records increasing peak temperatures from amphibolite (600–700 °C), transition (750–800 °C) and granulite (850–950 °C) zones. Titanium content of amphibole may be modified by retrograde fluid influx returning temperatures c. 200–300 °C lower than in non-altered domains. The comparison reveals that Zr-in-rutile thermometer in pelitic granulites seems to be more prone to post-peak resetting than the Ti-in-amphibole thermometry in nearby mafic rocks. This behavior is also confirmed by amphibole analyses from other UHT localities, where the performance of Ti-in-amphibole thermometry is comparable with that of Al-in-orthopyroxene in pelitic granulites. However, Ti-in-amphibole temperatures are underestimated in rutile-bearing samples and this limitation is not solely restricted to rocks containing high H2O contents as previously thought. Derived constraints on the diffusivity of Ti through amphibole demonstrate the robustness of the Ti-in-amphibole thermometer to later thermal disturbances. However, ad-hoc experiments are still necessary to improve the accuracy and precision of calibration and to extend its applicability. This advance will make mafic granulites routine targets for studies devoted to understanding the regional extent of UHT metamorphism.

Accretion of “young” Phanerozoic subcontinental lithospheric mantle triggered by back-arc extension—the case of the Ivrea-Verbano Zone

The subcontinental lithospheric mantle (SCLM) beneath Phanerozoic regions is mostly constituted by fertile lherzolites, which sharply contrast with cratonic mantle made of highly-depleted peridotites. The question of whether this chemical difference results from lower degrees of melting associated with the formation of Phanerozoic SCLM or from the refertilization of ancient depleted SCLM remains a subject of debate. Additionally, the timing and geodynamic environment of accretion of the fertile SCLM in many Phanerozoic regions are poorly constrained. We here document new geochemical and Nd-Hf isotopic data for orogenic lherzolite massifs from the Ivrea-Verbano Zone (IVZ), Southern Alps. Even though a few Proterozoic Re depletion ages are locally preserved in these mantle bodies, our data reveal that the IVZ lherzolitic massifs were “recently” accreted to the SCLM in the Upper Devonian (ca. 370 Ma) during Pangea amalgamation, with a petrochemical evolution characterized by low-degree (~ 5–12%) depletion and nearly contemporaneous pervasive to focused melt migration. The lithospheric accretion putatively took place through asthenospheric upwelling triggered by Variscan intra-continental extension in a back-arc setting related to the subduction of the Rheic Ocean. We thus conclude that the fertile sections of Phanerozoic SCLM can be accreted during “recent” events of back-arc continental extension, even where Os isotopes preserve memories of melting events in much older times.

Fracture network characterisation of the Balmuccia peridotite using drone-based photogrammetry, implications for active-seismic site survey for scientific drilling

The presence of discontinuities (e.g. faults, fractures, veins, layering) in crystalline rocks can be challenging for seismic interpretations because the wide range of their size, orientation, and intensity, which controls the mechanical properties of the rock and elastic wave propagation, resulting in equally varying seismic responses at different scales. The geometrical characterisation of adjacent outcrop discontinuity networks allows a better understanding of the nature of the subsurface rocks and aids seismic interpretation. In this study, we characterise the discontinuity network of the Balmuccia peridotite (BP) in the Ivrea–Verbano Zone (IVZ), northwestern Italy. This geological body is the focus of the Drilling the Ivrea–Verbano zonE (DIVE), an international continental scientific drilling project, and two active seismic surveys, SEismic imaging of the Ivrea ZonE (SEIZE) and high-resolution SEIZE (Hi-SEIZE), which aim to resolve the subsurface structure of the DIVE drilling target through high-resolution seismic imaging. For fracture characterisation, we developed two drone-based digital outcrop models (DOMs) at two different resolutions (10−3–10 m and 10−1–103 m), which allowed us to quantitatively characterise the orientation, size, and intensity of the main rock discontinuities. These properties affect the seismic velocity and consequently the interpretation of the seismic data. We found that (i) the outcropping BP discontinuity network is represented by three more sets of fractures with respect to those reported in the literature; (ii) the discontinuity sizes follow a power-law distribution, indicating similarity across scales, and (iii) discontinuity intensity is not uniformly distributed along the outcrop. Our results help to explain the seismic behaviour of the BP detected by the SEIZE survey, suggesting that the low P-wave velocities observed can be related to the discontinuity network, and provide the basic topological parameters (orientation, density, distribution, and aperture) of the fracture network unique to the BP. These, in turn, can be used for interpretation of the Hi-SEIZE seismic survey and forward modelling of the seismic response.

Underplated melts control sulfide segregation at the continental crust-mantle transition

Exposures of the Earth’s crust-mantle transition are scarce, thus, limiting our knowledge about the formation of subcontinental underplate cumulates, and their significance for metal storage and migration. Here, we investigated chalcophile metals to track sulfide crystallization within the Contact Series, an <150-m-thick pyroxenite-gabbronorite sequence, formed by mantle-derived melts, highlighting the boundary between the Balmuccia mantle peridotite and gabbronoritic Mafic Complex of the Ivrea-Verbano Zone. Within the Contact Series, numerous sulfides crystallized in response to the differentiation of mantle-derived underplated melts. Such sulfide-controlled metal differentiation resulted in anomalous Cu contents (up to ~380 ppm), compared to reference mantle (~19 ppm) and crustal samples (~1 ppm). We propose that the assimilation of continental crust material is a critical mechanism driving sulfide segregation and sulfide-controlled metal storage. Our results evidence that sulfides are trapped in the underplated mafic-ultramafic cumulates and that their enrichment in Cu may provide essential implications for crustal metallogeny.

Transition from orogenic-like to anorogenic magmatism in the Southern Alps during the Early Mesozoic: Evidence from elemental and Nd-Sr-Hf-Pb isotope geochemistry of alkali-rich dykes from the Finero Phlogopite Peridotite, Ivrea–Verbano Zone

The Ivrea-Verbano Zone (IVZ) in the westernmost sector of the Southern Alps is an iconic upper mantle to lower continental crust sequence of the Adriatic Plate and provides a geological window into the tectono-magmatic events that occurred at the Gondwana–Laurussia boundary from Late Paleozoic to Early Mesozoic. In this work, we document new geochemical and Nd-Sr-Hf-Pb isotopic data for Early Mesozoic alkali-rich dyke swarms which intruded the Finero Phlogopite Peridotite (northern IVZ) to provide geological constraints on the nature, origin and evolution of Early Mesozoic magmatism in the Southern Alps. The studied dykes are amphibole-phlogopite-bearing and show geochemical features varying between two end-member groups. A dyke group is characterized by HFSE-poor, Al-rich amphibole (Al2O3 up to 16 wt%) with high LILE and LREE contents, high radiogenic 87Sr/86Sr(i) (0.704732 to 0.704934) and low radiogenic Nd isotopes (εNd(i) from –0.1 to –0.7), which support the occurrence of significant amounts of recycled continental crust components in the parental mantle melts and impart an overall “orogenic-like” affinity. This dyke group was largely derived from metasomatized lithospheric mantle sources. The second group is HFSE-rich with Al-poorer amphibole enriched in LILE and LREE, low radiogenic 87Sr/86Sr(i) (0.703761–0.704103) and higher radiogenic Nd isotopes (εNd(i) from +3.4 to +5.4) pointing to an “anorogenic” alkaline affinity and asthenospheric to deep lithospheric mantle sources. Some dykes show both orogenic and anorogenic affinities, providing evidence that the orogenic-like magmatism in the IVZ predates the alkaline anorogenic magmatism. The Finero dyke swarms therefore record a geochemical change of the Early Mesozoic magmatism of the Southern Alps from orogenic-like magmatism, typical of post-collisional settings, to anorogenic alkaline magmatism, common in intraplate to extensional settings, and places a temporal correlation of Early Mesozoic magmatism in the IVZ to those in the eastern and central sectors of the Southern Alps.

New constraints on the shear wave velocity structure of the Ivrea geophysical body from seismic ambient noise tomography (Ivrea-Verbano Zone, Alps)

SUMMARY We performed seismic ambient noise tomography to investigate the shallow crustal structure around the Ivrea geophysical body (IGB) in the Ivrea-Verbano Zone (IVZ). We achieved higher resolution with respect to previous tomographic works covering the Western Alps, by processing seismic data collected by both permanent and temporary seismic networks (61 broad-band seismic stations in total). This included IvreaArray, a temporary, passive seismic experiment designed to investigate the IVZ crustal structure. Starting from continuous seismic ambient noise recordings, we measured and inverted the dispersion of the group velocity of surface Rayleigh waves (fundamental mode) in the period range 4–25 s. We obtained a new, 3-D vS model of the IVZ crust via the stochastic neighbourhood algorithm (NA), with the highest resolution between 3 to 40 km depth. The fast and shallow shear wave velocity anomaly associated with the IGB presents velocities of 3.6 km s−1 directly at the surface, in remarkable agreement with the location of the exposed lower-to-middle crustal and mantle outcrops. This suggests a continuity between the surface geological observations and the subsurface geophysical anomalies. The fast IGB structure reaches vS of 4 km s−1 at 20–25 km depth, at the boundary between the European and Adriatic tectonic plates, and in correspondence with the earlier identified Moho jump in the same area. The interpretation of a very shallow reaching IGB is further supported by the comparison of our new results with recent geophysical investigations, based on receiver functions and gravity anomaly data. By combining the new geophysical constraints and the geological observations at the surface, we provide a new structural interpretation of the IGB, which features lower crustal and mantle rocks at upper crustal depths. The comparison of the obtained vS values with the physical properties from laboratory analysis of local rock samples suggests that the bulk of the IGB consists of a combination of mantle peridotite, ultramafic and lower crustal rocks, bound in a heterogeneous structure. These new findings, based on vS tomography, corroborate the recent interpretation for which the Balmuccia peridotite outcrops are continuously linked to the IGB structure beneath. The new outcomes contribute to a multidisciplinary framework for the interpretation of the forthcoming results of the scientific drilling project DIVE. DIVE aims at probing the lower continental crust and its transition to the mantle, with two ongoing and one future boreholes (down to 4 km depth) in the IVZ area, providing new, complementary information on rock structure and composition across scales. In this framework, we constrain the upper crustal IGB geometries and lithology based on new evidence for vS, connecting prior crustal knowledge to recent active seismic investigations.

Structural, lithostratigraphic and thermal features of a Permian lower crust from the Western Italian Alps (Valpelline Series, Valle d’Aosta)

AbstractThe Valpelline Series (Dent-Blanche Tectonic System, Western Italian Alps) is a sector of lower continental crust, which consists of Permian migmatitic metapelite with different mineral assemblages (i.e., garnet-, cordierite- and orthopyroxene-bearing), minor amphibolite and marble, intruded by aplite and pegmatite. Widespread melt production in metapelite and locally in amphibolite occurred during the development of the regional foliation. The P–T conditions during migmatisation, estimated using conventional geothermobarometers, range between 800–900 °C and 0.5–0.8 GPa, with a difference of up to ∼50 °C between cordierite- and orthopyroxene-bearing migmatites, the latter reaching higher temperatures. The Valpelline Series shows rock types, metamorphic assemblages, P–T conditions and published ages of high-temperature regional metamorphism like the archetypal lower crust section of the Ivrea-Verbano Zone in the Southern Western Alps. The latter likely represents an external portion of the same extending lower crust, at the onset of the Tethyan rifting due to lithospheric extension and asthenospheric rising.

First Ce-Nd isotope measurements of middle and lower continental crust samples support massive lower crust recycling over Earth's history

Long-lived radioactive isotopes provide valuable information on the evolution of Earth's geological reservoirs. Coupled measurements of the 138La-138Ce and 147Sm-143Nd systems have attracted much interest, but some critical reservoirs, such as the deep continental crust, have yet to be investigated. To address this gap, we report Ce-Nd isotope measurements of 35 crustal samples from Canada, western Europe, and Siberia, with Archean to Phanerozoic ages. Most samples from western Europe and Siberia are interpreted to represent the deepest continental crust because of low SiO2 content (< 55 wt%) and generally mafic compositions. In contrast, Precambrian composites from the Canadian shield have highly felsic compositions (61–71 wt% SiO2) and represent the old upper continental crust. Xenoliths from the Massif Central (France) and samples from uplifted massifs in the Ivrea-Verbano zone (Italy) plot on the Ce-Nd mantle array. Precambrian Canadian composites and xenoliths from Udachnaya (Siberian Craton) plot to the left of the mantle array, with the Siberian samples farthest from the array. The La-Ce system yields a ca. 1.8 Ga age for the Siberian xenoliths, distinct from the ca. 2.7 Ga Sm-Nd errorchron age. The La-Ce system was most likely reset during a large-scale episode of delamination and rejuvenation of the Archean lower lithosphere established in Siberia. A change in the La/Ce ratio during this event explains an increasing deviation with time of the Ce-Nd isotopic composition of the Siberian deep crust from the mantle array. We introduce a new parameter (θ) to quantify the likelihood of any rock evolving away from the mantle array. It represents the angle between the evolution vector of a sample and the mantle array, and it is calculated from measured La/Ce and Sm/Nd ratios. Samples that plot to the right of the mantle array have θ > 0, whereas samples that plot to the left of the mantle array have θ < 0. The angle θ was used to compare our samples to a worldwide granulite compilation (1581 samples) because it is independent of the age and the initial isotopic composition of the samples. The majority of lower crust rocks have negative θ, suggesting they are prone to evolve to the left of the mantle array, but to a lesser extent than the Siberian samples that are characterized by the most negative θ values. Thus, the Siberian samples with the most unusual Ce-Nd isotopic compositions are the only group to plot close to the deep crust end-member calculated from classical mass-balance budget of the bulk silicate Earth. These results suggest that the lower crust sampled by Siberian xenoliths is only a minor component in the terrestrial deep crust. We use this database to estimate new parent/daughter ratios for the lower continental crust: 138La/142Ce = 0.00380 ±0.00007 and 147Sm/144Nd = 0.128 ±0.007. Measurements and models are finally reconciled when considering massive recycling of lower continental crust through Earth's history, totaling 3 to 4 present-day, continental crust masses.

Crystal plasticity and fluid availability govern the ability of titanite to record the age of deformation

Here, we study the relationships of titanite-hosting microdomains, intragrain chemical variations, microstructure and fluids with the aim of deciphering the reliability of titanite U–Pb dating to constrain the age of deformation in mylonitic rocks. We investigate these relationships in a post-Variscan amphibolite-facies shear zone developed in the mid-low continental crust (Ivrea-Verbano Zone, Southern Alps, Italy). Quantitative orientation analyses along with textural imaging of titanite are combined with trace-element analyses and U–Pb age dating. Titanite is studied in mm- to cm-scale layered rocks showing compositional variation consisting of alternating ‘amphibole-rich’ (i.e., amphibolites) and ‘clinopyroxene/plagioclase-rich’ domains (i.e., calc-silicates). Titanite from amphibole-rich domains shows predominance of crystal–plastic deformation features, as abrupt or progressive core-to-rim structures characterized by increasing lattice distortion and local dislocation density, associated with the development of abundant subgrains and rare newly nucleated grains. We suggest that these microstructures form while interacting with small amounts of fluids circulating along the grain boundaries. Consequently, locally the chemistry of titanite is changing. In the clinopyroxene/plagioclase-rich domains, titanite is mostly undeformed and rarely shows bending localized in discontinuous narrow rims/tips. In these domains, fluid-mediated replacement reactions are either rare or absent, as also indicated by weak chemical variations across and among grains. These observations suggest different reactivities with respect to the same P-T-fluid conditions of the two compositional domains, which coexist within the same sample at the thin section scale. U–Pb data show correlations with chemical and microstructural domains that differ as function of the composition of the microdomain. This correlation is more apparent within amphibole-rich domains where microstructures characterized by high lattice distortion/dislocations and/or subgrains show significant variations of REE, Zr, Y, Nb, U with respect to the low deformed portion of grains. These titanite domains define an isotopic population providing the youngest (Jurassic) lower intercept age. A less clear correlation between titanite chemistry and microstructures is observed in clinopyroxene/plagioclase-rich domains. Here, the rare titanites showing lattice distortion and minor Sr depletion define a population providing a similar Jurassic lower intercept age. Therefore, our results demonstrate that microstructurally and chemical calibrated U–Pb dating of titanite provides realistic ages of shear zone activity, only in case of predominance of crystal-plastic processes and of local interaction of titanite with small amounts of fluids focused along grain boundaries. Finally, the different footprints recorded by titanite grains strongly depend on the composition of cm- or mm-scale interlayered domains in which titanite occurs.

Tectono-metamorphic evolution of a post-variscan mid-crustal shear zone in relation to the Tethyan rifting (Ivrea-Verbano Zone, Southern Alps)

The Ivrea-Verbano Zone, in the Italian Southern Alps, samples a complete section of middle to lower continental crust, which records the transition between the Variscan and the Alpine Tethys rift-related tectonics. We present a structural-metamorphic characterization of one of the most important structures, the Forno-Rosarolo shear zone which drove the post-Variscan tectonic evolution of the Ivrea-Verbano Zone. The Forno-Rosarolo shear zone is a NE-SW-oriented, subvertical shear zone, mainly made of metapelites, amphibolites, calc-silicates and granulites involved in anastomosed proto- to ultra-mylonite layers enveloping weakly deformed lenses. Mylonitic deformation developed on an already folded and sheared metamorphic (Variscan) complex but predates brittle fracturing and faulting, locally associated with pseudotachylites.In present day orientation, the kinematic indicators point to a sinistral sense of shear. Removing the Alpine ∼90° tilt of the IVZ crustal section, this kinematic results in a former extensional shear zone. Investigations on the kinematic of the flow in the mylonites revealed a non-coaxial deformation with a major component of pure shear acting together with simple shear. Metamorphic conditions of the wall rocks vary from upper amphibolite (SE, footwall) to granulite facies (NW, hanging wall). Estimations from the mineral assemblages within the mylonites point to amphibolite facies conditions during deformation (∼650 °C and ∼5.5 kbar).In absence of precise geochronological constraints, we discuss the role of the Forno-Rosarolo shear zone within the Tethyan rifting system in comparison with recent reconstructions. Kinematic data and metamorphic conditions allow to constrain the development of the Forno-Rosarolo shear zone mylonitic deformation during the intermediate phase of the Tethyan rift known as “thinning mode”. This stage was characterized by general shear conditions (pure shear between 70% and 50%) suggesting a phase of transition from a symmetric to an asymmetric configuration of rift.

3-D imaging of the Balmuccia peridotite body (Ivrea–Verbano zone, NW-Italy) using controlled source seismic data

SUMMARY We provide new results from a controlled-source seismic experiment on the deepest part of the Val Sesia crust–mantle section of the Ivrea–Verbano zone (IVZ) in the Italian Alps. The IVZ is a tilted, almost complete section through the continental crust and exposes gabbros and peridotites in the structurally deepest level, coinciding with high-resolution gravity anomalies imaging the Ivrea geophysical body. The seismic experiment SEIZE (SEismic imaging of the Ivrea ZonE) was conducted along two crossing profiles: an NNE-SSW profile of ∼11 km length and an E-W profile of ∼16 km length. 432 vibration points were recorded with 110 receivers resulting in 24 392 traveltime picks. Inversion methods using Markov chain Monte Carlo techniques have been used to derive an isotropic 3-D P-wave velocity model based on first break traveltimes (refracted phases) from controlled source seismic data. Resulting seismic P-wave velocities (Vp ) range from 4.5 to 7.5 km s−1, with an expected general trend of increasing velocities with depth. A sharp velocity change from low Vp in the West to high Vp in the East marks the Insubric Zone (ISZ), the Europe–Adria plate boundary. The most prominent feature of the 3-D tomography model is a high-velocity body (Vp increases from 6 to 7.5 km s−1) that broadens downwards. Its pointy shape peaks the surface East of Balmuccia at a location coincident with the exposed Balmuccia peridotite. Considering rock physics, high-resolution gravity and other geophysical data, we interpret this high-velocity body as dominantly composed of peridotite. The dimension of this seismically imaged peridotite material is far bigger than interpreted from geological cross-sections and requires a revision of previous models. The interpretation of ultramafic bodies in the IVZ as fragments of mantle peridotites interfingered in the crust during pre-Permian accretion is not supported by the new data. Instead, we revive a model that the contact between the Balmuccia peridotite and the Permian mafic magmas might represent a fossil continental crust–mantle transition zone.

C–O–H fluid-melt-rock interaction in graphitic granulites and problems of quantifying carbon budget in the lower continental crust

Estimates on the geological carbon cycle are subject to large uncertainties that can be reduced by thorough observation of rocks. In this contribution, we focus specifically on C–O–H fluid-melt-rock interactions in graphitic metapelitic granulites and on their bearing to the carbon budget of granulitic roots of continents. We provide robust microstructural and thermometric constraints on the coexistence of anatectic silicate melts and C–O–H fluids up to near ultrahigh temperature conditions in the archetypal crustal section of Ivrea-Verbano Zone (IVZ, Italian Alps). Fluid inclusions in garnet are investigated before and after high-temperature experiments, and contain considerable proportions of CO2, CH4, N2, but lower H2O than predicted for graphitic systems at granulite facies. When comparing and contrasting the melt compositions obtained by Perple_X and rhyolite-MELTS with natural melts from IVZ, a much better match is obtained by the former, questioning the choice of rhyolite-MELTS for modelling melting equilibria of metasedimentary rocks and for quantifying carbon budget of the lower crust. Overall, data show that assuming only a limited extent of fluid-melt immiscibility in the deep crust contradicts the evidence from natural rocks and prompts to an incomplete view of actual carbon behavior and carbonic fluids. The available experimental dataset on CO2 solubility in felsic melts cannot be used to interpret the volatile budget of melt inclusions in graphitic migmatites and granulites, as most solubility experiments were conducted under carbonate-saturated (i.e. highly oxidizing) conditions which maximize CO2 content of melt, compared to graphitic (i.e. more reducing) protoliths. As a consequence, thermodynamic models still cannot account for all the complexities related with interactions among H2O–CO2–CH4 ternary fluids, H2O- and CO2-bearing anatectic melts and graphite-bearing residues in graphitic metapelites. Targeted experimental studies are therefore crucial to boost substantial computational efforts, before any precise estimates on carbon budget and fluxes in the lower anatectic crust can be made.

The peridotite-pyroxenite sequence of Rocca d'Argimonia (Ivrea-Verbano Zone, Italy): Evidence for reactive melt flow and slow cooling in the lowermost continental crust

The present study investigates the origin of a ∼ 400 m thick peridotite-pyroxenite sequence, locally known as Rocca d'Argimonia, which is encased within the lowest levels of the ∼8 km thick Mafic Complex, a gabbro-dioritic batholith from the lower continental crust section of the Ivrea-Verbano Zone (Southern Alps). The broad purpose of this work is to elucidate the evolution of mantle-derived magmas emplaced in the lowermost continental crust. We also wish to provide new geochronological constraints on the intrusion and the cooling evolution of the Lower Mafic Complex, which are still mostly unknown.Zircon grains separated from two gabbronorites enclosing the peridotite-pyroxenite sequence document a whole reset of the original UPb isotopic system at 286 ± 2 Ma. A Rocca d'Argimonia peridotite provided a consistent age of 296 ± 19 Ma, based on a LuHf clinopyroxene-amphibole alignment (initial εHf = +7.9). These results are interpreted to date the early cooling evolution of the Lower Mafic Complex. The clinopyroxene-amphibole peridotite pair also gave a RbSr alignment corresponding to an age of 250 ± 16 Ma (initial 87Sr/86Sr = 0.7045), which is reconciled with the subsequent, slow cooling evolution of the Lower Mafic Complex.The Rocca d'Argimonia peridotites (dunites to harzburgites and lherzolites) have whole-rock εNd(286 Ma) and εHf(286 Ma) values ranging from +0.4 to −1.1 and from +8.8 to +3.4, respectively, and 87Sr/86Sr(286 Ma) ratios varying from 0.7048 to 0.7060. The peridotites are crosscut by gabbronorite dykes and associated with olivine-free orthopyroxene-dominated pyroxenites. The gabbronorite dykes and the pyroxenites share the Nd-Hf-Sr signature of the peridotites. The gabbronorites enclosing the peridotite-pyroxenite sequence differ in the relatively enriched Nd-Hf-Sr isotopic fingerprint (εNd(286 Ma) = −4.5 to −5.6, εHf(286 Ma) = +0.2 to −2.5, 87Sr/86Sr(286 Ma) = 0.7073 to 0.7086). Remarkably, pyroxenes from the pyroxenites typically contain higher Cr2O3 than the peridotite counterparts (e.g., 0.23–0.37 wt% Cr2O3 in peridotite orthopyroxenes and 0.43–0.50 wt% Cr2O3 in pyroxenite orthopyroxenes).We propose that the peridotite-pyroxenite sequence formed by reactive melt percolation through an olivine-rich spinel-bearing matrix. In particular, the relatively high Cr2O3 content of the pyroxenitic pyroxenes is attributed to redistribution of Cr released by the breakdown of accessory spinel that was originally present in the olivine-dominated matrix, in response to a high melt/solid ratio. The melt residual after the crystallization of the pyroxenites continued to percolate the olivine-dominated matrix, thereby forming spinel-bearing lherzolites, harzburgites and dunites with the progression of the melt migration process. The chemically most primitive dunite records Nd-Hf-Sr isotopic disequilibrium conditions, thereby providing further information about the process of reactive melt flow. It is inferred that the original olivine-dominated matrix crystallized from a mantle magma having a slightly enriched Nd-Hf-Sr isotopic signature, whereas the percolating melts overall had a relatively depleted Nd-Hf-Sr isotopic fingerprint, despite being variably contaminated by the continental crust.

Metasomatism in the Finero Phlogopite Peridotite: New insights from C and N concentrations and δ13C - δ11B signatures

The interconnection between element recycling and mantle metasomatism may provide powerful insights to unravel the mechanisms governing the transfer, partitioning and residence times of volatiles among the different Earth's reservoirs. In this study, we provide new constraints on the incorporation of C and N in the rock-forming minerals of the Finero Phlogopite Peridotite (FPP) from the Ivrea-Verbano Zone (Southern Alps, Italy), representing a natural laboratory to shed lights on the amounts of volatiles that may be fixed in a deeply metasomatized sub-continental lithospheric mantle (SCLM). Constraints on the relationships between metasomatic events and volatile element enrichment were given by in-situ LA-ICP-MS and micro-Raman investigations together with C (δ13C) and B (δ11B) isotopes.All rock-forming minerals (olivine, clinopyroxene, amphibole and phlogopite) show enrichments in C (270–690 ppm) and N (7–20 ppm) compared to depleted mantle values. The high volatile contents in clinopyroxene and in amphibole are associated with enrichment in incompatible elements (e.g., Li, B, Rb, Sr, Pb, Th, U) and high LaN/LuN ratios (up to 42.0) thus confirming that the metasomatic agent affecting the FPP derive form a hydrous melt with crustal affinity enriched in C and N. Micro-Raman investigations document the occurrence of graphite, carbonate and N2 inclusions, specifically in olivine and clinopyroxene, that are likely responsible for the high CN concentrations reported here. The δ13C signatures measured in mineral separates range from −12.6 to −22.7 ‰, suggesting an overprinting of the C mantle signature by a 12C-rich crustal metasomatic agent. High temperature decarbonation processes, simulated by heating the samples up to 900 °C, indicate that up to10 ppm of C can be residually stored in the SCLM, whereas N is completely released even at lower temperature. The LA-MC-ICP-MS in-situ δ11B signature of amphibole is −5.4 ± 4.1 ‰ (2σ, n = 18), whereas bulk δ11B of olivine and clinopyroxene is −5.51 ± 0.01 ‰ (2σ) and − 1.82 ± 0.06 ‰ (2σ), respectively. Modelling shows that residual slab melts plus variable amounts of 11B-rich serpentinite-derived fluids may account for these B isotope imprints. A subduction-related metasomatizing event during the FPP's evolution is thus proposed to explain the geochemical features documented in our study, in agreement with previous works. At global scale, our results suggest that the metasomatized SCLM might significantly contribute to the so-called hidden C reservoir invoked to match the global C mass balance calculations.

Multistage tectono-stratigraphic evolution of the Canavese Intracontinental Suture Zone: New constraints on the tectonics of the Inner Western Alps

The Canavese Intracontinental Suture Zone (CISZ) within the Inner Western Alps represents the remnant of a long-lived minor subduction zone involving a narrow, thinned continental crust/oceanic lithosphere seaway between two continental domains of the Adria microplate (i.e., the Sesia Zone and the Ivrea-Verbano Zone). As opposed to many suture zones, the CISZ mostly escaped pervasive tectonic deformation and metamorphism, thus preserving the original stratigraphy and allowing the relationships between tectonics and sedimentation to be defined. Through detailed geological mapping (1:5000 scale), structural analysis, stratigraphic and petrographic observations, we document evidences for the late Paleozoic to late Cenozoic tectonic evolution of the CISZ, showing that it played a significant role in the context of the tectonic evolution of the Inner Western Alps region from the early to late Permian Pangea segmentation, to the Jurassic Tethyan rifting, and up to the subduction and collisional stages, forming the Western Alps. The site of localization/formation of the CISZ was not accidental but associated with the re-use of structures inherited from regional-scale wrench tectonics related to the segmentation of Pangea, and from the subsequent extensional tectonics related to the Mesozoic rifting, as documented by crosscutting relationships between stratigraphic unconformities and tectonic features. Our findings document that evidences derived from stratigraphy, facies indicators, and relationships between tectonics and sedimentation in the shallow crustal portions of suture zones, such in the CISZ, are important to better constrain the tectonic history of those metamorphic orogenic belts around the world in which evolutionary details are commonly complicated by high-strain deformation and metamorphic transformations.