Ivrea-Verbano Zone Research Articles

Role of inherited compositional and structural heterogeneity in shear zone development at mid-low levels of the continental crust (the Anzola shear zone; Ivrea-Verbano Zone, Southern Alps)

The formation of shear zones is crucial to understand the deformation of the crust and the evolution of rifted margins. However, despite their intrinsic importance, a detailed characterization of the compositional and structural patterns of shear zones is often lacking, resulting in poorly constrained models of shear initiation. In this contribution, we reconstruct the pre-shearing lithological, structural and metamorphic proprieties of rocks forming a major, rift-related shear zone with the aim to assess the conditions promoting the strain localization. We focus on the Anzola shear zone, a major extensional structure from one of the best-preserved cross-sections through the middle to lower continental crust of a fossil passive margin, the Ivrea-Verbano Zone (Southern Alps, Italy). Until now, the Anzola shear zone is believed to have developed within a rheologically hard and isotropic gabbro rather than in the surrounding weaker and anisotropic volcano-sedimentary metamorphic sequence. New geological mapping shows that significant pre-existing heterogeneities related to the composition and deformation structures, characterize the Anzola shear zone. Field relationships and geochemistry reveal that the (ultra-)mylonitic rocks overprinted a multi-lithological sequence that have already experienced Variscan deformation and late Variscan High-Temperature metamorphism, at the boundary between amphibolite and granulite facies. Our in-depth trace elements study is shown to be a powerful tool in reconstructing the pre-shearing relationships between wall rocks and mylonites and determining the protoliths of tectonites. Estimated P-T conditions indicate that mylonitic deformation started at high temperature (~820 °C) with presence of melt and continued as solid-state deformation down to amphibolite facies (~650 °C), following a retrograde path. We argue that strain localization was promoted by the combination of rheological boundaries derived from pre-existing conditions, including: i) compositional and structural anisotropies of the volcano-sedimentary metamorphic sequence contrasted by ii) the close intrusion of a nearly isotropic gabbro and iii) the presence of melt within the metamorphic boundary depicted by the transition between granulites (dominated by anhydrous minerals) and amphibolite facies (dominated by hydrous minerals). Our findings finally suggest that pre-existing significant heterogeneities relate to rock composition, deformation and metamorphism represent the preferential loci for strain localization controlling the initiation and development of rift-related structures in the mid to lower crust of passive margins. • Field mapping and chemical analysis of an extensional shear zone from middle/lower crust. • Trace element analysis reveals composite protoliths of the mylonites. • Anisotropic and/or hydrated rocks are folded with isotropic and/or anhydrous layers. • Mylonites postdates folding and upper-amphibolite-granulite facies metamorphism. • Pre-existing contrasted rheological boundaries promoted strain localization.

Thorium zoning in monazite: A case study from the Ivrea–Verbano zone, NW Italy

AbstractMetamorphism and partial melting of the lower crust is commonly assumed to cause depletion in heat producing elements (HPEs; K, U, Th). In the deep crust, volumetrically subordinate metasedimentary layers, which are source to crustal granites (sensu lato), host the majority of Th ± U, primarily within the REE + Th + U + Y phosphate mineral monazite. We examine the spatial and temporal distribution of Th within monazite grains in metasedimentary rocks from the lower crustal section of the Ivrea–Verbano Zone (Italy), using textural, compositional and geochronological data. We link this to outcrop and regional scale trends described by in‐field gamma‐ray spectrometry data (in‐field GRS) for the purpose of understanding how Th distribution is controlled by progressive metamorphism and partial melting. In‐field GRS data shows that the whole rock budget of Th does not change between granulite facies rocks and their unmelted equivalents but is significantly lower in rocks that have undergone more significant melt loss at ultra‐high temperature (UHT) conditions. Concurrently, the bulk Th budget of monazite increases with metamorphic grade to granulite facies conditions and is greatly reduced in UHT samples. Monazite geochronology returns dates mostly in the range 240–320 Ma with two main peaks at circa 290 and 270 Ma. Textural and chemical constraints indicate that these dates record the timing of pre‐peak to peak metamorphic conditions. Amphibolite facies monazite compositional zones are absent from granulite facies monazite, in contrast to examples from lower‐pressure terranes. This is consistent with the expanded stability of allanite relative to monazite with increasing pressure having an important role in determining the internal structure, composition and extent of inheritance of monazite in going from amphibolite facies to granulite facies rocks. We propose high‐pressure granulites should preserve less monazite inherited from amphibolite facies conditions than low‐pressure granulites. Monazite is preserved at all metamorphic grades and presents a mineralogical mechanism for retaining Th in residual deep crust during partial melting and after melt loss.

Assembly and Tectonic Evolution of Continental Lower Crust: Monazite Petrochronology of the Ivrea‐Verbano Zone (Val Strona di Omegna)

AbstractMetasediments are common constituents of exhumed lower‐to‐mid‐crustal granulite terranes; understanding their emplacement is significant for the assembly and tectonic evolution of deep continental crust. Here, we report a monazite U‐Th‐Pb petrochronological investigation of the Variscan Ivrea‐Verbano Zone (IVZ) (Val Strona di Omegna section)—an archetypal section of lower crust. Monazite Th‐Pb dates from 11 metapelitic samples decrease with structural depth from 310 to 285 Ma for amphibolite‐facies samples to <290 Ma for granulite‐facies samples. These dates exhibit a time‐resolved variation in monazite trace‐element composition, dominated by the effects of plagioclase and garnet partitioning. Monazite growth under prograde to peak metamorphic conditions began as early as 316 ± 2 Ma. Amphibolite‐facies monazite defines a trend consistent with progressively decreasing garnet modal abundances during decompression and cooling starting at ∼310 Ma; the timing of the onset of exhumation decreases to ∼290 Ma at the base of the amphibolite‐facies portion of the section. Structurally lower, granulite‐facies monazite equilibrated under garnet‐present pressure‐temperature conditions at <290 Ma, with monazite (re)crystallization persisting until at least ∼260 Ma. Combined with existing detrital zircon U‐Pb dates, the monazite data define a <30 Myr duration between deposition of clastic sediments and their burial and heating, potentially to peak amphibolite‐to‐granulite‐facies conditions. Similarly brief timescales for deposition, burial and prograde metamorphism of lower crustal sediments have been reported from continental magmatic arc terranes—supporting the interpretation that the IVZ represents sediments accreted to the base of a Variscan arc magmatic system >5 Myr prior to the onset of regional extension and mafic magmatism.

Role of Inherited Compositional and Structural Heterogeneity in Shear Zone Development at Mid-Low Levels of the Continental Crust (the Anzola Shear Zone; Ivrea-Verbano Zone, Southern Alps)

Role of Inherited Compositional and Structural Heterogeneity in Shear Zone Development at Mid-Low Levels of the Continental Crust (the Anzola Shear Zone; Ivrea-Verbano Zone, Southern Alps)

High-resolution seismic reflection survey crossing the Insubric Line into the Ivrea-Verbano Zone: Novel approaches for interpreting the seismic response of steeply dipping structures

A high-resolution seismic reflection survey has been conducted across the Insubric Line from the Sesia Zone into the Ivrea-Verbano Zone (IVZ), where a remarkably complete cross-section of lower continental crust is exposed. The survey was carried out in preparation for the DIVE (Drilling the Ivrea-Verbano zonE) project, which was recently approved by the International Continental Scientific Drilling Program (ICDP). DIVE aims to gain new insights into the characteristics of the lower continental crust through targeted drilling, sampling, and borehole logging. A key borehole is planned near the Insubric Line at Balmuccia, where the deepest parts of the lower continental crust are exposed. As such, the primary objective of this seismic survey was to explore whether the sub-vertical structures prevailing at the surface can be expected to continue at depth or whether there are any indications for major flattening or fault-related offsets. Correspondingly, the acquisition and processing of the seismic reflection data were geared towards revealing weak backscattered events from local heterogeneities associated with the prevailing sub-vertical structural grain. The migrated sections, contain coherent backscattered events to a depth of ~1 km, which form numerous short lineaments that seem to align sub-vertically. To substantiate this observation, we have generated synthetic seismic reflection surveys for canonical models of sub-vertical structures associated with Gaussian- and binary-distributed heterogeneities. Both the observed and synthetic seismic data were then subjected to energy-based attribute analysis as well as geostatistical estimations of the structural aspect ratios and the associated dips. The results of these quantitative interpretation approaches are indicative of the overall consistency between the synthetic and the observed seismic data and, hence, support the original qualitative interpretation of the latter in that the sub-vertical structural grain evident at the surface seems to prevail throughout the imaged part of the upper crust.

TRIASSIC-JURASSIC RIFT-RELATED DEFORMATION AND TEMPERATURE-TIME EVOLUTION OF THE FOSSIL ADRIATIC MARGIN: A REVIEW FROM OSSOLA AND STRONA DI OMEGNA VALLEYS (IVREA-VERBANO ZONE)

The direct observation and investigation of rift-related structures at the mesoscale is uncommon. Hence, detailed constraints on the evolution of the main faults and shear zones developed during crustal extension are not always available.The Ivrea-Verbano Zone, in the Italian Southern Alps, samples remnants of the former lower crust of the rifted margins surrounding the Alpine Tethys and therefore provides the opportunity to directly investigate rift-related tectonics. Here, several shear zones have been recognized and interpreted as related to Mesozoic rifting. However, even if there is a general agreement with this interpretation, the precise age of activity of many of those shear zones is not well constrained.In this paper we present a review of the Triassic-Jurassic geochronological and thermochronological data available for two sections of the Ivrea-Verbano Zone, the Strona di Omegna and Ossola valleys, where at least two extensional shear zones are exposed. Ductile deformation occurred under amphibolite-facies conditions and it was alternatively attributed to late Variscan deformation or to Triassic-Jurassic rifting-related tectonics.We discuss the available chronological data and the different interpretations provided for the shear zones considering also new geochronological studies on other lower crustal shear zones exposed in other sectors of the Ivrea-Verbano Zone. This review allows us to strengthen the more recent interpretations indicating that these shear zones are important tectonic structures related to Late Triassic-Jurassic deformation in the lower crust of the Adriatic margin.

Joint Seismic and Gravity Data Inversion to Image Intra-Crustal Structures: The Ivrea Geophysical Body Along the Val Sesia Profile (Piedmont, Italy)

We present results from a joint inversion of new seismic and recently compiled gravity data to constrain the structure of a prominent geophysical anomaly in the European Alps: the Ivrea Geophysical Body (IGB). We investigate the IGB structure along the West-East oriented Val Sesia profile at higher resolution than previous studies. We deployed 10 broadband seismic stations at 5 km spacing for 27 months, producing a new database of ∼1000 high-quality seismic receiver functions (RFs). The compiled gravity data yields 1 gravity point every 1–2 km along the profile. We set up an inversion scheme, in which RFs and gravity anomalies jointly constrain the shape and the physical properties of the IGB. We model the IGB’s top surface as a single density and shear-wave velocity discontinuity, whose geometry is defined by four, spatially variable nodes between far-field constraints. An iterative algorithm was implemented to efficiently explore the model space, directing the search toward better fitting areas. For each new candidate model, we use the velocity-model structures for both ray-tracing and observed-RFs migration, and for computation and migration of synthetic RFs: the two migrated images are then compared via cross-correlation. Similarly, forward gravity modeling for a 2D density distribution is implemented. The joint inversion performance is the product of the seismic and gravity misfits. The inversion results show the IGB protruding at shallow depths with a horizontal width of ∼30 km in the western part of the profile. Its shallowest segment reaches either 3–7 or 1–3 km depth below sea-level. The latter location fits better the outcropping lower crustal rocks at the western edge of the Ivrea-Verbano Zone. A prominent, steep eastward-deepening feature near the middle of the profile, coincident with the Pogallo Fault Zone, is interpreted as inherited crustal thickness variation. The found density and velocity contrasts of the IGB agree with physical properties of the main rock units observed in the field. Finally, by frequency-dependent analysis of RFs, we constrain the sharpness of the shallowest portion of the IGB velocity discontinuity as a vertical gradient of thickness between 0.8 km and 0.4 km.

Anatectic granites in their source region: A comparison between experiments, thermodynamic modelling and nanogranitoids

Anatectic granites represent the products of crustal maturation by partial melting, therefore, studies on their origin are crucial to understanding how continental crust forms and evolves. This contribution compares composition of melts obtained from three tools allowing the investigation of anatectic granites in their source region: experiments, phase equilibria modelling and nanogranitoids. Benchmark experimental runs are selected along with nanogranitoids from the Ivrea–Verbano Zone (IVZ; NW Italy) and compositions of melt are then calculated using the same bulk rock compositions considered in experimental and natural case studies. The possible effects of melt loss, diverse P–T paths and variable equilibration volumes on melt composition are modeled for IVZ rocks. Results are then compared with those of previous comparative studies to verify the presence of systematic compositional divergences and similarities. Collectively, all data reflect pressure and temperature conditions (≈700–950 °C and 0.3–1.5 GPa) typical of the suprasolidus orogenic crust, with source regions ranging from metapelite, metagreywacke to metagranite.Contents of FeO, MgO, K2O and, to a lesser extent, CaO in the melt show the largest divergences. FeO, MgO and CaO tend to be enriched in experimental melts and nanogranitoids, whereas K2O is mostly overestimated in calculated melts.After a consideration of the common pitfalls for each approach we discuss the potential causes for the main geochemical discrepancies. The most extreme deviations of FeO and MgO, often observed in low-T (≤ 750 °C) melts, seem to be related to the selected melt models, which have been mostly calibrated on the basis of high-T (≥800 °C) experiments. Some experimental melts tend to be strongly depleted in K2O due to the metastable persistence of micas at much higher temperatures in laboratory runs, whereas the compositions of nanogranitoids in some investigated rocks may be partially controlled by diffusion in the melt. Disequilibrium melting of plagioclase may explain the CaO enrichment observed in some experimental and natural melts. When local equilibration volumes are considered in the calculations, a better match is observed for some oxides (e.g. FeO and MgO) between nanogranitoids and their calculated counterparts, even though other significant divergences remain unaffected. The combination of diverse petrologic tools, with proper recognition of their shortcomings, represent the best approach to shed light into the complexities of melting processes in nature which control the composition of anatectic granites at the source.

Contaminating melt flow in magmatic peridotites from the lower continental crust (Rocca d'Argimonia sequence, Ivrea–Verbano Zone)

Abstract. The lower continental crust section of the Ivrea–Verbano Zone (Italian Alps) was intruded by a ∼ 8 km thick gabbroic–dioritic body (Ivrea Mafic Complex) in the Upper Carboniferous–Lower Permian, in conjunction with the post-collisional transtensional regime related to the Variscan orogeny. In the deepest levels of the Ivrea Mafic Complex, several peridotite–pyroxenite sequences considered of magmatic origin are exposed. We present here a petrological–geochemical investigation of the peridotites from the largest magmatic ultramafic sequence of the Ivrea Mafic Complex, locally called Rocca d'Argimonia. In spite of the widespread subsolidus re-equilibration under granulite facies conditions, most likely reflecting a slow cooling evolution in the lower continental crust, the Rocca d'Argimonia peridotites (dunites to harzburgites and minor clinopyroxene-poor lherzolites) typically retain structures and microstructures of magmatic origin. In particular, the harzburgites and the lherzolites typically show poikilitic orthopyroxenes enclosing partially dissolved olivine and minor spinel. Olivine has forsterite proportion diminishing from the dunites to the harzburgites and the lherzolites (90 mol % to 85 mol %) and negatively correlating with δ18O (+5.8 ‰ to +6.6 ‰). Gabbronorite dykes locally crosscut the peridotites and show millimetre-scale thick, orthopyroxenite to websterite reaction zones along the contact with host rocks. We propose that the Rocca d'Argimonia peridotites record a process of reactive melt flow through a melt-poor olivine-rich crystal mush or a pre-existing dunite. This process was most likely responsible for the olivine dissolution shown by the poikilitic orthopyroxenes in the harzburgites–lherzolites. We infer that the reactively migrating melts possessed a substantial crustal component and operated at least at the scale of ∼ 100 m.

On the formation of magmatic sulphide systems in the lower crust by long‐lived mass transfer through the lithosphere: Insights from the Valmaggia pipe, Ivrea Verbano Zone, Italy

AbstractThe lower crustal domain of the Ivrea‐Verbano Zone (NW Italy) hosts five ~300‐m‐wide pipe‐like ultramafic intrusions that are metasomatized and exhibit Ni‐Cu‐PGE sulphide mineralization. To better constrain the role of metasomatism in the ore genesis, we studied the best‐preserved pipe at Valmaggia which was emplaced 249 Myrs ago. Phlogopite 40Ar/39Ar analyses show that the pipe was infiltrated by metasomatic fluids derived from the subcontinental lithospheric mantle (SCLM) in two pulses at ~208 Ma and ~189 Ma which introduced sulphides into the pipe. Consequently, the pipe repeatedly acted as a preferred path for mass transfer from the SCLM into the lower crust over >60 Myrs (i.e., emplacement to second metasomatic pulse). Uplifted block margins, such as the Ivrea‐Verbano Zone, are potentially important exploration targets for magmatic sulphides. We argue that exploration strategies should focus on structures such as pipes that can focus metasomatic agents during ascent through the lithosphere.

Mantle-Derived Corundum-Bearing Felsic Dykes May Survive Only within the Lower (Refractory/Inert) Crust: Evidence from Zircon Geochemistry and Geochronology (Ivrea–Verbano Zone, Southern Alps, Italy)

Corundum-rich (up to 55 vol.%) felsic dykes formed with albite, ±K-feldspar, ±hercynite and ±biotite-siderophyllite cut the lower crustal rocks exposed in the Ivrea–Verbano Zone (NW Italy). Zircon is an abundant accessory mineral and its investigation through laser ablation-inductively coupled plasma (multi-collector)-mass spectrometer (LA-ICP-(MC)MS) has allowed results to directly constrain the timing of emplacement, as well as petrology and geochemistry of parental melts. Zircons are characterized by very large concentration in rare earth elements (REE), Th, U, Nb and Ta, and negative Eu anomaly. U–Pb analysis points to Norian emplacement ages (223 ± 7 Ma and 224 ± 6 Ma), whereas large positive εHf(t) values (+13 on average) indicate a derivation from depleted to mildly enriched mantle source. The mantle signature and the corundum oversaturation were preserved thanks to limited crustal contamination of the host, high-temperature refractory granulites and mafic intrusives. According to the geochemical data and to the evidence of the development of violent explosions into the conduits, it is proposed that dykes segregated from peraluminous melts produced by exsolution processes affecting volatile-rich differentiates during alkaline magmatism. This work provides robust evidence about the transition of the geochemical affinity of Southern Alps magmatism from orogenic-like to anorogenic during Norian time, linked to a regional uprising of the asthenosphere and change of tectonic regime.

High Temperature Reduced Granulite-Facies Nature of Garnetites in the Khabarny Mafic–Ultramafic Massif, Southern Urals: Evidence from Fluid and Mineral Analyses

Abstract High-grade metamorphic rocks underlying the intrusive layered dunite–pyroxenite–gabbronorite East-Khabarny Complex (EKC) are integrated in the complex Khabarny mafic–ultramafic Massif in the Sakmara Allochthon zone in the Southern Urals. These rocks are associated with high-temperature shear zones. Garnetites from the upper part of the metamorphic unit close to the contact with EKC gabbronorite are chemically and texturally analysed to estimate their formation conditions and fluid regime. Fluids provide crucial information of formation conditions and evolution of these garnetites during high-grade metamorphism, and are preserved in channel positions within Si6O1812- rings of cordierite, and in fluid inclusions in quartz and garnet. Minerals and fluid inclusions of the garnetites are studied with X-ray fluorescence spectrometry, electron microprobe analyses, Raman spectroscopy, and microthermometry. The garnetites mainly consist of garnet (up to 80 vol. %), cordierite and quartz. Accessory minerals are rutile, ilmenite, graphite, magnetite and cristobalite. Granulite-facies metamorphic conditions of the garnetites are estimated with the garnet–cordierite–sillimanite–quartz geothermobarometer: temperatures of 740 to 830 ˚C and pressures of 770–845 MPa. The average garnet composition in end-member concentrations is 48·5 mole % almandine (±3·9), 34·7 mole % pyrope (±3·3), 10·3 mole % spessartine (±1·1), 1·8 mole % grossular (±1·5), and 1·5 mole % andradite (±1·5). The cordierite electron microprobe analyses reveal an average Mg2+ fraction of 0·79 ± 0·01 in the octahedral site. Relicts of a strong positive temperature anomaly (up to 1000 ˚C) are evidenced by the preservation of cristobalite crystals in garnet and the high titanium content of quartz (0·031 ± 0·008 mass % TiO2) and garnet (0·31 ± 0·16 mole % end-member Schorlomite-Al). The fluid components H2O, CO2, N2 and H2S are detected in cordierite, which correspond to a relatively oxidized fluid environment that is common in granulites. In contrast, a highly reduced fluid environment is preserved in fluid inclusions in quartz nodules, which are mono-fluid phase at room temperature and composed of CH4 (>96 mole %) with locally minor amounts of C2H6, N2, H2S and graphite. The fluid inclusions occur in homogeneous assemblages with a density of 0·349 to 0·367 g·cm-3. The CH4-rich fluid may represent peak-temperature metamorphic conditions, and is consistent with temperature estimation (∼1000 ˚C) from Ti-in-garnet and Ti-in-quartz geothermometry. Tiny CH4-rich fluid inclusions (diameter 0·5 to 2 µm) are also detected by careful optical analyses in garnet and at the surface of quartz crystals that are included in garnet grains. Graphite in fluid inclusions precipitated at retrograde metamorphic conditions around 300–310 ± 27 ˚C. Aragonite was trapped simultaneously with CH4-rich fluids and is assumed to have crystallized at metastable conditions. The initial granulite facies conditions that led to the formation of a cordierite and garnet mineral assemblage must have occurred in a relative oxidized environment (QFM-buffered) with H2O–CO2-rich fluids. Abundant intrusions or tectonic emplacement of mafic to ultramafic melts from the upper mantle that were internally buffered at a WI-buffered (wüstite–iron) level must have released abundant hot CH4-rich fluids that flooded and subsequently dominated the system. The origin of the granulite-facies conditions is similar to peak-metamorphic conditions in the Salda complex (Central Urals) and the Ivrea–Verbano zone (Italian Alps) as a result of magmatic underplating that provided an appearance of a positive thermal anomaly, and further joint emplacement (magmatic and metamorphic rocks together) into upper crustal level as a high temperature plastic body (diapir).

New gravity data and 3-D density model constraints on the Ivrea Geophysical Body (Western Alps)

SUMMARY We provide a high-resolution image of the Ivrea Geophysical Body (IGB) in the Western Alps with new gravity data and 3-D density modelling, integrated with surface geological observations and laboratory analyses of rock properties. The IGB is a sliver of Adriatic lower lithosphere that is located at shallow depths along the inner arc of the Western Alps, and associated with dense rocks that are exposed in the Ivrea-Verbano Zone (IVZ). The IGB is known for its high seismic velocity anomaly at shallow crustal depths and a pronounced positive gravity anomaly. Here, we investigate the IGB at a finer spatial scale, merging geophysical and geological observations. We compile existing gravity data and we add 207 new relative gravity measurements, approaching an optimal spatial coverage of 1 data point per 4–9 km2 across the IVZ. A compilation of tectonic maps and rock laboratory analyses together with a mineral properties database is used to produce a novel surface rock-density map of the IVZ. The density map is incorporated into the gravity anomaly computation routine, from which we defined the Niggli gravity anomaly. This accounts for Bouguer Plate and terrain correction, both considering the in situ surface rock densities, deviating from the 2670 kg m–3 value commonly used in such computations. We then develop a 3-D single-interface crustal density model, which represents the density distribution of the IGB, including the above Niggli-correction. We retrieve an optimal fit to the observations by using a 400 kg m–3 density contrast across the model interface, which reaches as shallow as 1 km depth below sea level. The model sensitivity tests suggest that the ∼300–500 kg m–3 density contrast range is still plausible, and consequently locates the shallowest parts of the interface at 0 km and at 2 km depth below sea level, for the lowest and the highest density contrast, respectively. The former model requires a sharp density discontinuity, the latter may feature a vertical transition of densities on the order of few kilometres. Compared with previous studies, the model geometry reaches shallower depths and suggests that the width of the anomaly is larger, ∼20 km in west–east direction and steeply E–SE dipping. Regarding the possible rock types composing the IGB, both regional geology and standard background crustal structure considerations are taken into account. These exclude both felsic rocks and high-pressure metamorphic rocks as suitable candidates, and point towards ultramafic or mantle peridotite type rocks composing the bulk of the IGB.

Porphyroclasts: Source and Sink of Major and Trace Elements During Deformation-Induced Metasomatism (Finero, Ivrea-Verbano Zone, Italy)

Petrographic and geochemical data for mylonites from a metric-scale shear zone in mantle peridotites from the Finero massif (Southern Alps) record large mineralogical and geochemical modifications compared to surrounding coarse-grained ultramafic rocks, which were pervasively deformed in presence of hydrous melts. The mylonites are composed by olivine and orthopyroxene and, less frequently, clinopyroxene, phlogopite, and pargasite porphyroclasts enclosed in a fine-grained matrix of phlogopite and olivine, with subordinate amounts of orthopyroxene, clinopyroxene, pargasite, and chromite. P-T estimates indicate that deformation occurred under granulite- to upper-amphibolite facies conditions. Field relationships and U-Pb dating indicate that the shear zone was active during Lower Jurassic and/or later, in an extensional setting at the western margin of the Adria plate, which led to the opening of the Alpine Tethys. The major and trace element composition of the porphyroclasts in the mylonites significantly differ from those in the hosting coarse-grained ultramafics. Porphyroclasts were chemically active during deformation acting as source (diffusion-out) or sink (diffusion-in) for some trace elements. The chemical modifications were promoted by the interaction with aqueous fluids and the composition varied from mantle- (enriched in Ni, Co, Li, Na, REE, Y, and Sr) to crustal-derived (enriched in Zn, K, Al, Ti, and Fe).

Timescales and Mechanisms of Crystal-mush Rejuvenation and Melt Extraction Recorded in Permian Plutonic and Volcanic Rocks of the Sesia Magmatic System (Southern Alps, Italy)

AbstractSilicic calderas can evacuate 100 to >1000 km3 of rhyolitic products in a matter of days to months, leading to questions on pre-eruptive melt generation and accumulation. Whereas silicic plutonic units may provide information on the igneous evolution of crystal-mush bodies, their connection with volcanic units remains enigmatic. In the Ivrea–Verbano Zone of the southern Alps, the plumbing system of a Permian rhyolitic caldera is exposed to a depth of about 25 km in tilted crustal blocks. The upper-crustal segment of this magmatic system (also known as the Sesia Magmatic System) is represented by the Valle Mosso pluton (VMP). The VMP is an ∼260 km3 composite silicic intrusion ranging from quartz-monzonite to high-silica leucogranite (∼67–77 wt% SiO2), which intrudes into roughly coeval rhyolitic products of the >15 km diameter Sesia Caldera. In the caldera field, the emplacement of a large, crystal-rich rhyolite ignimbrite(s) (>400 km3) is followed by eruption of minor volumes (1–10 km3) of crystal-poor rhyolite. Here, we compare silicic plutonic and volcanic units of the Sesia Magmatic System through a combination of geochemical (X-ray fluorescence, inductively coupled plasma mass spectrometry and electron microprobe analyses) and petrological (rhyolite-MELTS, trace element and diffusion modeling) tools to explore their connection. Textural and compositional features shared by both VMP and crystal-rich ignimbrites imply thermal rejuvenation of crystal-mush as the mechanism to create large volumes of eruptible rhyolitic magma. Bulk-rock composition of crystal-rich rhyolite erupted during the caldera collapse overlaps that of the bulk VMP. Quartz and plagioclase from these two units show resorbed cores and inverse zoning, with Ti- and anorthite-rich rims, respectively. This indicates crystallization temperatures in rims >60 °C higher than in cores (780–820 versus ∼720 °C), if temperature is the sole parameter responsible for zonation, suggesting heating and partial dissolution of the crystal-framework. Decrease in crystallinity associated with thermal energy input was calculated through rhyolite-MELTS and indicates lowering of the mush crystal fraction below the rheological lock-up threshold, which probably promoted eruptive activity. Also, after the climatic eruption, Si-rich melts in the Sesia Magmatic System were produced by extraction of interstitial melt from un-erupted, largely crystalline mush. Regarding both textures and chemical variations, we interpret the deep quartz-monzonite unit of the VMP as a compacted silicic cumulate. Fractionated melts extracted from this unit were emplaced as a leucogranite cupola atop the VMP, generating the final internal architecture of the silicic intrusion, or alternatively erupted as minor post-caldera, crystal-poor rhyolite. Ti-in-quartz diffusion profiles in thermally rejuvenated units of the Sesia Magmatic System demonstrate that the process of reheating, mobilization and eruption of crystal-mush took place rapidly (c. 101–102 years). A protracted cooling history is instead recorded in the diffusion timescales of quartz from the silicic cumulate units (c. 104–106 years). These longer timescales encompass the duration of evolved melt extraction from the cumulate residue. We argue that the VMP preserves a complex record of pre-eruptive processes, which span mechanisms and timescales universally identified in volcanic systems and are consistent with recently proposed numerical models.

Crystallographic preferred orientation, magnetic and seismic anisotropy in rocks from the Finero peridotite, Ivrea-Verbano Zone, Northern Italy – Interplay of anisotropy contributions from different minerals

Mineral alignment can provide valuable information on a rock's geological history. Many structural, tectonic, or geodynamic studies rely on determining the preferred orientation of rock-forming minerals, either by direct texture determination or using anisotropy measurements as proxies. Robust interpretations of anisotropy require detailed understanding of its relationship with texture. This study investigates the texture, magnetic and seismic anisotropies of ultramafic rocks from the Ivrea-Verbano Zone (IVZ), Northern Italy, with a special focus on the interplay between anisotropy contributions carried by different minerals. Texture was obtained from electron backscatter diffraction, magnetic anisotropy from low- and high-field anisotropy of susceptibility measurements, and seismic anisotropy from P- and S-wave velocities along the three principal directions of the macroscopic fabric, using 22 mm diameter cores. Texture-based models of magnetic anisotropy agree well with corresponding measurements. Larger variability is observed for seismic anisotropy. Possible explanations are (1) velocities were measured in three directions, insufficient to determine the full elasticity tensor, and (2) seismic anisotropy depends on grain boundaries, pores, and cracks in addition to crystallographic preferred orientation. The contributions to magnetic and seismic anisotropy of each constituent mineral can interfere positively, leading to a larger overall anisotropy, or negatively, resulting in a weaker anisotropy. Seismic anisotropy is mostly controlled by olivine, whereas magnetic anisotropy can be dominated by olivine, pyroxene, or hornblende. This study illustrates that, although the relationship between texture and anisotropy can be complex, it is possible to quantitatively predict magnetic, and to a lesser degree seismic, anisotropy from texture data. Even though both magnetic and seismic anisotropies are common proxies for texture, they are not coaxial, and no simple relationship between the degree of magnetic and seismic anisotropy was found. These results underline the importance of understanding different minerals' contributions to anisotropy in upper mantle rocks.

3D rock fabric analysis using micro-tomography: An introduction to the open-source TomoFab MATLAB code

The study of rock fabric properties (orientation, planar, linear, anisotropy) is key to unravelling the geological processes that generated them. With advancements in data acquisition and treatment, X-ray micro-computed tomography (μXCT) represents a powerful method to analyse the shape preferred orientation (SPO) of rock-forming elements, including minerals, aggregates, and pores, in the three-dimensional space. After reconstruction and segmentation of μXCT images, we developed a novel protocol to construct and analyse the fabric tensor, a second-rank symmetric tensor constructed using the orientation and the length of the three characteristic axes of each grain (simplified to a best fit ellipsoid). The analysis of the fabric tensor permits calculation of mean principal directions and associated confidence ellipses, and quantifies the degree of anisotropy (P′) and the shape (T) of the fabric ellipsoid by eigenvalue and eigenvector analysis.We implement this method in the TomoFab open-source MATLAB package. The code integrates a graphical user interface (GUI) that allows the visualisation of the full set of ellipsoid orientation, shape, and size. Density plots and contouring can be utilised to identify fabrics graphically, and a full set of fabric parameters can be calculated based on the analysis of the fabric tensor and/or the analysis of each principal direction orientation tensor.We demonstrate the versatility of TomoFab with synthetic datasets and a field- and laboratory-based investigation of a sample presenting a magmatic foliation and lineation, collected in the Mafic Complex within the lower crustal section of the Ivrea-Verbano Zone (North Italy). In the light of these developments, we stress that μXCT represents a pertinent tool for rock fabric analysis to characterise the SPO of rock components. This approach can be performed parallel or complementary to other rock fabric quantification methods (e.g., AMS, EBSD) and applied to various rock types. TomoFab is freely available for download at https://github.com/benpetri/tomofab.

Evidence of subduction-related components in sapphirine-bearing gabbroic dykes (Finero phlogopite–peridotite): Insights into the source of the Triassic–Jurassic magmatism at the Europe–Africa boundary

A gabbroic dyke swarm containing magmatic sapphirine occurs in the Finero phlogopite–peridotite (FPP), one of the major mantle massifs in the Ivrea–Verbano Zone (IVZ; western Southern Alps). Sapphirine is part of a particular mineral assemblage, including plagioclase, titanian pargasite, titanian phlogopite, and Cl-rich apatite; the latter mineral hosts calcite inclusions. The dykes cut the mantle foliation at a high angle, are bounded by orthopyroxenite layers, and show symmetric internal banding, represented by two outer hornblendite selvages and an inner leucogabbro band. The sapphirine occurs in up to 3 cm-thick irregular patches in both hornblendite salvages, along with Al-rich amphibole and green spinel.We present major and trace elements of minerals and bulk rock, as well as mineral O, Sr, and Nd isotopic compositions of dykes and the host peridotite from two different outcrops in the FPP area.Our data show that early melt migration developed through porous flow within cm-thick channels and was characterised by orthopyroxene dissolution. Following progressive percolation and reaction, the melt became silica saturated with segregation of orthopyroxenite in the centres of the channels. The banded internal structure of the dykes was caused by three different evolutionary stages, involving opening and enlargement of the conduits. The sapphirine and green spinel segregation took place at T > 1000 °C, in the presence of melt with transient composition, which interstitially migrated and reacted with the cumulus minerals to form the hornblendite layers. The mineral chemistry of the newly-formed amphiboles indicates that the sapphirine parental melt was Al-rich, depleted to strongly depleted in Hf, Zr, Nb, Ta, Ti, Sc, V, and middle and heavy rare earth elements, and characterised by a positive Eu anomaly and (Zr/Hf)N < 1. These data suggest a parental melt with a significant amount of normative plagioclase. However, the studied veins do not show evidence of plagioclase assimilation, and we argue that this process could have occurred in magmatic bodies that are not outcropping today to the surface or in the melt source.The δ18O values of vein amphiboles and plagioclases vary from 6.9 to 8.6‰ SMOW, which is well above the mantle range, even when considering fractionation upon cooling. Given that orthopyroxene from the wall has “normal” mantle δ18O values (5.8‰), reaction with the host metasomatised peridotite cannot be responsible for the heavy δ18O signature, and the latter must have been imparted by crustal components deeper in the mantle.Our petrographic and geochemical evidence demonstrates that the northern IVZ records an extremely prolonged release, from the Variscan orogenic cycle to the Mesozoic exhumation, of K-H2O-rich mantle-derived melts, mixed with subduction-related components. This finding provides valuable insights into the Triassic–Jurassic magmatism and the geodynamic environment at the Europe–Africa boundary.

Evolution of mantle melts intruding the lowermost continental crust: constraints from the Monte Capio–Alpe Cevia mafic–ultramafic sequences (Ivrea–Verbano Zone, northern Italy)

This study presents a new petrological–geochemical data set for the Monte Capio and Alpe Cevia mafic–ultramafic sequences, which are exposed in the deepest levels of the Ivrea–Verbano Zone. These sequences are composed of a peridotite core, with dunite in the center, mantled by minor orthopyroxene-dominated pyroxenites and subordinate hornblende gabbronorites. Amphibole is ubiquitous in the peridotites and the pyroxenites (≤ 15 vol % and 10–40 vol %, respectively), and the peridotite–pyroxenite associations are frequently crosscut by amphibole-rich (45–90 vol %) veins/dykes showing sinuous-to-sharp planar boundaries towards host rocks. The whole-rock Mg# [100 × Mg/(Mg + Fetot2+)] decreases from the peridotites to the pyroxenites and the crosscutting amphibole-rich dykes (84–81, 80–77, and 73–66, respectively), consistently with the Mg# variations shown by included orthopyroxene, clinopyroxene, and amphibole. Olivine has relatively low forsterite and NiO amounts (84–78 mol % and ≤ 0.14 wt%), and spinel is characterized by low Cr# [100 × Cr/(Cr + Al)] of 7–24. The anorthite content of plagioclase varies from 91 to 88 mol% in plagioclase-bearing pyroxenites to 91–75 mol% in amphibole-rich dykes. The chondrite-normalized REE patterns of amphibole from peridotites and pyroxenites show nearly flat MREE–HREE, no evident Eu anomaly, and LREE that are slightly depleted to slightly enriched with respect to MREE. Amphibole from the amphibole-rich veins/dykes exhibits slight LREE depletion. Whole-rock and amphibole separates show substantial variations in initial Nd–Sr isotopic compositions (e.g., whole-rock eNd calculated at 290 Ma ranges from − 0.3 to − 4.7), irrespective of the rock-type and of incompatible element amphibole compositions. We propose that the Monte Capio–Alpe Cevia dunites formed by cooling of magma lenses that intruded the lowermost continental crust of the Ivrea–Verbano Zone. The chemically evolved signature of the dunites documents earlier crystallization of chemically primitive dunites at lower levels, or olivine fractionation within the dunites during melt ascent. Associated pyroxene-bearing peridotites show a magmatic evolution ruled by reaction of a melt-poor crystal mush with migrating melts relatively rich in SiO2 and H2O, which developed orthopyroxene and amphibole at the expenses of olivine ± clinopyroxene. These migrating melts may be reconciled with those feeding the crosscutting amphibole-rich veins/dykes, whose compositions suggest formation by chemically evolved H2O-rich basalts with an arc-type incompatible trace-element fingerprint. Unraveling the origin of the Monte Capio–Alpe Cevia pyroxenites is hampered by the complex open-system magmatic evolution, which also included assimilation of material released by basement metasediments and/or involvement of primary melt batches with different compositions.

Lead oxide nanospheres in seismically deformed zircon grains

The presence of Pb oxide nanospheres in zircon from mylonites of the Ivrea-Verbano Zone (IVZ), Northern Italy is described for the first time. Isotopic dating of zircon reveals detrital cores with ages scattered from 960 to 320 Ma, and metamorphic rims with a mean age of 280 ± 4 Ma. Zircon crystals, derived from samples that contain pseudotachylite formed during seismic events, have planar fractures (PFs) and planar deformation bands (PDBs). The PDBs are associated with straight dislocation arrays in glide configuration. Detrital zircon cores have mottled diffraction contrast in TEM bright field images, indicative of irradiation damage and/or annealing in radiation damaged zircon. Lead oxide nanospheres up to 9 nm in diameter, as recognized in TEM, occur in detrital cores but not in metamorphic rims of the zircon. No relationship was observed between the nanospheres and PDBs or PFs. This is the first report of lead oxide nanospheres in zircon from a Phanerozoic metamorphic rock, as well as being the first observed case of such nanoinclusions being composed of Pb oxide, rather than native Pb. The formation of the nanospheres is attributed to high-temperature metamorphism, with subsequent seismic activity having little effect on distribution or preservation.